From a35dbc067ce20ecab47389de50aab753625879fb Mon Sep 17 00:00:00 2001 From: Fabian Neumann Date: Mon, 11 Nov 2019 16:29:27 +0100 Subject: [PATCH] build_natura: add cutouts as input (closes #95) --- Snakefile | 8 ++++++-- scripts/build_natura_raster.py | 2 +- 2 files changed, 7 insertions(+), 3 deletions(-) diff --git a/Snakefile b/Snakefile index 2474d516..ed5c86cb 100644 --- a/Snakefile +++ b/Snakefile @@ -117,13 +117,17 @@ if config['enable']['build_cutout']: script: "scripts/build_cutout.py" else: rule retrieve_cutout: - output: directory("cutouts/{cutout}") + output: + directory("cutouts/europe-2013-era5"), + directory("cutouts/europe-2013-sarah") script: 'scripts/retrieve_cutout.py' if config['enable']['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" script: "scripts/build_natura_raster.py" else: diff --git a/scripts/build_natura_raster.py b/scripts/build_natura_raster.py index 78c7e804..e85c561e 100644 --- a/scripts/build_natura_raster.py +++ b/scripts/build_natura_raster.py @@ -52,6 +52,6 @@ if __name__ == "__main__": xs, Xs, ys, Ys = zip(*(determine_cutout_xXyY(cutout) for cutout in cutout_names)) 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.rasterize(natura, pixelWidth=100, pixelHeight=100, output=snakemake.output[0])