1. 在win系统下将shape文件转为nc文件,精度为0.1X0.1 (anaconda shell)
gdal_rasterize -of netCDF -burn 1 -tr 0.1 0.1 ALA.shp ALA_maskout.nc
在linux环境下处理(WSL)
2. 激活环境py39
source activate py39
3. 跳转访问win下的文件夹
cd /mnt
4. 利用cdo将ALA.nc文件线性插值到目标网格中,经纬度lon/lat为3600/1801,与coup文件保持一致。
cdo -f nc remapbil,r3600x1801 ALA_maskout.nc remap_ALA_maskout.nc
5. 将coup.nc也进行一下插值,为了保持纬度相同
cdo -f nc remapbil,r3600x1801 global/Coupling/coup_1980.nc remap_coup_1980.nc
6.将得到的remap_ALA_maskout.nc文件用于裁剪数据
cdo ifthen remap_ALA_maskout.nc remap_coup_1980.nc ALA_coup_1980.nc
完成!
reference: