Python根据shp裁剪nc

您可以使用以下代码利用Python中的xarray和geopandas库来裁剪netCDF文件：

import xarray as xr
import geopandas as gpd

# Load the shapefile as a GeoDataFrame
gdf = gpd.read_file('path/to/shapefile.shp')

# Load the netCDF file as an xarray dataset
ds = xr.open_dataset('path/to/netcdffile.nc')

# Crop the netCDF file to the shapefile extent
ds_cropped = ds.sel(lat=slice(gdf.bounds.miny.max(), gdf.bounds.maxy.min()),
                    lon=slice(gdf.bounds.minx.min(), gdf.bounds.maxx.max()))

# Save the cropped netCDF file to a new file
ds_cropped.to_netcdf('path/to/croppedfile.nc')

在这里，我们首先将形状文件加载为GeoDataFrame对象。然后，我们将netCDF文件加载为xarray数据集对象。接下来，我们将netCDF数据集对象裁剪为形状文件的范围，并将其保存到新文件中。

请注意，此代码假定您的netCDF文件是具有经度和纬度坐标的二维数据数组。如果您正在处理三维数据（例如，带有时间维度的数据），则可能需要进一步裁剪。