From 8af09d9bec1e3ed6268047d01cd188192edade40 Mon Sep 17 00:00:00 2001 From: Philipp Glaum Date: Thu, 28 Jul 2022 15:21:21 +0200 Subject: [PATCH] fix raster format error due to import --- scripts/build_ship_raster.py | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/scripts/build_ship_raster.py b/scripts/build_ship_raster.py index 12ad64a6..d34a1c10 100644 --- a/scripts/build_ship_raster.py +++ b/scripts/build_ship_raster.py @@ -51,13 +51,13 @@ if __name__ == "__main__": cutouts = snakemake.input.cutouts xs, Xs, ys, Ys = zip(*(determine_cutout_xXyY(cutout) for cutout in cutouts)) - + with zipfile.ZipFile(snakemake.input.ship_density) as zip_f: zip_f.extract("shipdensity_global.tif") - ship_density=xarray.open_dataarray("shipdensity_global.tif", engine="rasterio") + ship_density=xarray.open_rasterio("shipdensity_global.tif") os.remove("shipdensity_global.tif") - ship_density=ship_density.drop("band").sel(x=slice(min(xs),max(Xs)), y=slice(max(Ys),min(ys))) - + 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])