From 6d7eae8ad28169f73f399197bc46eaad85e15c62 Mon Sep 17 00:00:00 2001 From: Fabian Date: Thu, 20 Apr 2023 19:56:23 +0200 Subject: [PATCH] build_ship_raster: save reduced ship raster as tif file --- rules/build_electricity.smk | 4 ++-- scripts/build_ship_raster.py | 6 +++--- 2 files changed, 5 insertions(+), 5 deletions(-) diff --git a/rules/build_electricity.smk b/rules/build_electricity.smk index 0074a8a9..133ceb6e 100644 --- a/rules/build_electricity.smk +++ b/rules/build_electricity.smk @@ -172,7 +172,7 @@ rule build_ship_raster: ], ), output: - RESOURCES + "shipdensity_raster.nc", + RESOURCES + "shipdensity_raster.tif", log: LOGS + "build_ship_raster.log", resources: @@ -202,7 +202,7 @@ rule build_renewable_profiles: ) ), ship_density=lambda w: ( - RESOURCES + "shipdensity_raster.nc" + RESOURCES + "shipdensity_raster.tif" if "ship_threshold" in config["renewable"][w.technology].keys() else [] ), diff --git a/scripts/build_ship_raster.py b/scripts/build_ship_raster.py index 59610ade..9d9dda0b 100644 --- a/scripts/build_ship_raster.py +++ b/scripts/build_ship_raster.py @@ -46,7 +46,7 @@ import logging import os import zipfile -import rioxarray as rio +import rioxarray from _helpers import configure_logging from build_natura_raster import determine_cutout_xXyY @@ -64,10 +64,10 @@ if __name__ == "__main__": with zipfile.ZipFile(snakemake.input.ship_density) as zip_f: zip_f.extract("shipdensity_global.tif") - with rio.open_rasterio("shipdensity_global.tif") as ship_density: + with rioxarray.open_rasterio("shipdensity_global.tif") as ship_density: ship_density = ship_density.drop(["band"]).sel( x=slice(min(xs), max(Xs)), y=slice(max(Ys), min(ys)) ) - ship_density.to_netcdf(snakemake.output[0]) + ship_density.rio.to_raster(snakemake.output[0]) os.remove("shipdensity_global.tif")