I have a netcdf 700MB, dp(time=169, lat=2250, lon=1670) and I want to clip
it with a shapefile.
data = xarray.open_dataset('mynetcdf.nc')
data.rio.set_spatial_dims(x_dim="lon", y_dim="lat", inplace=True)
data.rio.write_crs("epsg:4326", inplace=True)
Shapefile = geopandas.read_file('mask1.shp', crs="epsg:4326")
data_clipped = data.rio.clip(Shapefile.geometry.apply(mapping),
Shapefile.crs, drop=False)
data_clipped.to_netcdf('new.nc')
and the clipped one becomes 2gb.
_______________________________________________
Geoserver-users mailing list
Please make sure you read the following two resources before posting to this
list:
- Earning your support instead of buying it, but Ian Turton:
http://www.ianturton.com/talks/foss4g.html#/
- The GeoServer user list posting guidelines:
http://geoserver.org/comm/userlist-guidelines.html
If you want to request a feature or an improvement, also see this:
https://github.com/geoserver/geoserver/wiki/Successfully-requesting-and-integrating-new-features-and-improvements-in-GeoServer
Geoserver-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/geoserver-users
