如何通过Shapefile与CDO从NetCDF数据提取巴基斯坦区域？

你之前的remap命令报错，核心原因是**remap无法直接将Shapefile作为参数**，它需要的是网格定义（如匹配的NetCDF/GRIB文件网格、手动指定的网格字符串）。要通过Shapefile提取目标区域，正确流程是先生成与输入数据匹配的掩码文件，再用掩码过滤数据。

正确操作步骤 #

1. 生成匹配输入数据的掩码文件 #

先将Shapefile转换为和你的temp.nc同网格的掩码文件，巴基斯坦区域标记为1，其余区域为0：

cdo remapcon,temp.nc -shape2grid,pakistan.shp pak_mask.nc

• shape2grid,pakistan.shp：将Shapefile转换为网格格式数据
• remapcon,temp.nc：把Shapefile的网格重映射到和temp.nc完全一致的网格
• 输出的pak_mask.nc就是和输入数据匹配的掩码文件

2. 用掩码过滤数据并替换缺失值 #

将原数据与掩码相乘保留目标区域值，再把区域外的缺失值替换为0（不需要可省略setmisstoc,0）：

cdo setmisstoc,0 -mul temp.nc pak_mask.nc masked_pak.nc

• mul temp.nc pak_mask.nc：原数据与掩码相乘，区域外值变为缺失值
• setmisstoc,0：将所有缺失值替换为0

一步式简化命令 #

不想单独生成掩码文件的话，可合并两步为一条命令：

cdo setmisstoc,0 -mul temp.nc -remapcon,temp.nc -shape2grid,pakistan.shp - masked_pak.nc

其中-表示将前一步输出作为后一步输入。

额外注意事项 #

• 确保pakistan.shp的配套文件（.shx、.dbf等）和主文件在同一目录
• 如果NetCDF数据与Shapefile坐标系不一致，先用GDAL转换Shapefile坐标系（以WGS84为例）：

ogr2ogr -t_srs EPSG:4326 pakistan_wgs84.shp pakistan.shp

之后用转换后的pakistan_wgs84.shp生成掩码即可。

内容的提问来源于stack exchange，提问作者Sabina Munir