From 350af4a273d8f6a813c70a1e949af20ed4701762 Mon Sep 17 00:00:00 2001 From: Jonas Hoersch Date: Wed, 19 Dec 2018 10:27:46 +0100 Subject: [PATCH] build_renewable_profiles: Compute potentials in EPSG:3035 to ensure a metric of meters --- scripts/build_renewable_profiles.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/scripts/build_renewable_profiles.py b/scripts/build_renewable_profiles.py index 6c906e13..53ca3475 100644 --- a/scripts/build_renewable_profiles.py +++ b/scripts/build_renewable_profiles.py @@ -50,7 +50,7 @@ def downsample_to_coarse_grid(bounds, dx, dy, mask, data): def calculate_potential(gid): feature = gk.vector.extractFeature(snakemake.input.regions, where=gid) - ec = gl.ExclusionCalculator(feature.geom, srs=4326, pixelRes=1e-3) # about 100m in EU, it also works in LAEA 3035 + ec = gl.ExclusionCalculator(feature.geom) corine = config.get("corine", {}) if isinstance(corine, list): @@ -99,7 +99,7 @@ if __name__ == '__main__': with Pool(initializer=init_globals, initargs=(bounds, dx, dy), maxtasksperchild=20, processes=snakemake.config['atlite'].get('nprocesses', 2)) as pool: - regions = gk.vector.extractFeatures(snakemake.input.regions, onlyAttr=True).iloc[:10] + regions = gk.vector.extractFeatures(snakemake.input.regions, onlyAttr=True) buses = pd.Index(regions['name'], name="bus") widgets = [ pgb.widgets.Percentage(),