Linux and CDO处理

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:

利用gdal+cdo+shp将全球范围的nc文件裁剪为中国大陆区域 - 知乎 (zhihu.com)

你可能感兴趣的:(linux,运维,服务器)