fix raster format error due to import
This commit is contained in:
parent
6d070a19e0
commit
8af09d9bec
@ -54,10 +54,10 @@ if __name__ == "__main__":
|
|||||||
|
|
||||||
with zipfile.ZipFile(snakemake.input.ship_density) as zip_f:
|
with zipfile.ZipFile(snakemake.input.ship_density) as zip_f:
|
||||||
zip_f.extract("shipdensity_global.tif")
|
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")
|
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])
|
ship_density.to_netcdf(snakemake.output[0])
|
||||||
|
|
||||||
|
Loading…
Reference in New Issue
Block a user