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

build_natura_raster: Fix snakemake specifics

Browse Source
This commit is contained in:
Jonas Hoersch 2018-12-19 15:06:59 +01:00
parent 70857f955c
commit 8f171396cb
1 changed files with 4 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
scripts/build_natura_raster.py
View File

@ -1,16 +1,17 @@
import numpy as np
import atlite import atlite
from osgeo import gdal from osgeo import gdal
import geokit as gk import geokit as gk
def determine_cutout_xXyY(cutout_name): 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 x, X, y, Y = cutout.extent
dx = (X - x) / (cutout.shape[1] - 1) dx = (X - x) / (cutout.shape[1] - 1)
dy = (Y - y) / (cutout.shape[0] - 1) dy = (Y - y) / (cutout.shape[0] - 1)
return [x - dx/2., X + dx/2., y - dy/2., Y + dy/2.] 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()]) cutout_names = np.unique([res['cutout'] for res in snakemake.config['renewable'].values()])
xs, Xs, ys, Ys = zip(*(determine_cutout_xyXY(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[0])