shp文件裁剪netcdf python

code•2022-4-26•java•阅读54

shp文件裁剪netcdf

import xarray
import rioxarray
import geopandas
from shapely.geometry import mapping
#先导入必要的库

xds = xarray.open_dataset(".nc")
#使用xarray打开你的nc文件
print(xds)
#查看变量字段
xds = xds[['tp']].transpose('time', 'latitude', 'longitude')
#根据你的字段调整这里的内容
xds.rio.set_spatial_dims(x_dim="longitude", y_dim="latitude", inplace=True)
xds.rio.write_crs("EPSG:4326", inplace=True)


#加载shp文件
geodf = geopandas.read_file(".shp")
#利用shp裁剪
clipped = xds.rio.clip(geodf.geometry.apply(mapping), geodf.crs)
#将文件保存为nc格式
clipped.to_netcdf('.nc', mode='w', format="NETCDF4")