From 8f171396cb1cd6d958f539495836a43543ab1c42 Mon Sep 17 00:00:00 2001 From: Jonas Hoersch Date: Wed, 19 Dec 2018 15:06:59 +0100 Subject: [PATCH] build_natura_raster: Fix snakemake specifics --- scripts/build_natura_raster.py | 7 ++++--- 1 file changed, 4 insertions(+), 3 deletions(-) diff --git a/scripts/build_natura_raster.py b/scripts/build_natura_raster.py index b3a01f98..d0ad7cdf 100644 --- a/scripts/build_natura_raster.py +++ b/scripts/build_natura_raster.py @@ -1,16 +1,17 @@ +import numpy as np import atlite from osgeo import gdal import geokit as gk def determine_cutout_xXyY(cutout_name): - cutout = atlite.Cutout(cutout_name, cutout_dir="../cutouts") + cutout = atlite.Cutout(cutout_name, cutout_dir="cutouts") x, X, y, Y = cutout.extent dx = (X - x) / (cutout.shape[1] - 1) dy = (Y - y) / (cutout.shape[0] - 1) return [x - dx/2., X + dx/2., y - dy/2., Y + dy/2.] -cutout_names = np.unique([res['cutout'] for res in config['renewable'].values()]) -xs, Xs, ys, Ys = zip(*(determine_cutout_xyXY(cutout) for cutout in cutout_names)) +cutout_names = np.unique([res['cutout'] for res in snakemake.config['renewable'].values()]) +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])