bjoern/pypsa-eur
1
0
Fork 0
Code Issues Pull Requests Actions 2 Packages Projects Releases Wiki Activity

build_natura: add cutouts as input (closes #95)

Browse Source
This commit is contained in:
Fabian Neumann 2019-11-11 16:29:27 +01:00
parent 7cd3b8a42b
commit a35dbc067c
2 changed files with 7 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
Snakefile
View File

@ -117,13 +117,17 @@ if config['enable']['build_cutout']:
script: "scripts/build_cutout.py" script: "scripts/build_cutout.py"
else: else:
rule retrieve_cutout: rule retrieve_cutout:
output: directory("cutouts/{cutout}") output:
directory("cutouts/europe-2013-era5"),
directory("cutouts/europe-2013-sarah")
script: 'scripts/retrieve_cutout.py' script: 'scripts/retrieve_cutout.py'
if config['enable']['build_natura_raster']: if config['enable']['build_natura_raster']:
rule build_natura_raster: rule build_natura_raster:
input: "data/bundle/natura/Natura2000_end2015.shp" input:
natura="data/bundle/natura/Natura2000_end2015.shp",
cutouts=expand("cutouts/{cutouts}", **config['atlite'])
output: "resources/natura.tiff" output: "resources/natura.tiff"
script: "scripts/build_natura_raster.py" script: "scripts/build_natura_raster.py"
else: else:

2
scripts/build_natura_raster.py
View File

@ -52,6 +52,6 @@ if __name__ == "__main__":
xs, Xs, ys, Ys = zip(*(determine_cutout_xXyY(cutout) for cutout in cutout_names)) xs, Xs, ys, Ys = zip(*(determine_cutout_xXyY(cutout) for cutout in cutout_names))
xXyY = min(xs), max(Xs), min(ys), max(Ys) xXyY = min(xs), max(Xs), min(ys), max(Ys)
natura = gk.vector.loadVector(snakemake.input[0]) natura = gk.vector.loadVector(snakemake.input.natura)
extent = gk.Extent.from_xXyY(xXyY).castTo(3035).fit(100) extent = gk.Extent.from_xXyY(xXyY).castTo(3035).fit(100)
extent.rasterize(natura, pixelWidth=100, pixelHeight=100, output=snakemake.output[0]) extent.rasterize(natura, pixelWidth=100, pixelHeight=100, output=snakemake.output[0])