CDO —maskregion 提取任意边界区域

下载了欧洲再分析资料ERA5后，想从中提取某个区域范围内的值。网上辗转寻找方法，发现资料很少，在此记录下来，希望可以帮到大家。
CDO——Climate Data Operator，气象中很好用的处理工具，是各种常用函数的集合，功能强大，想要实现的目的，基本几条语句就可以解决！

提取矩形边界区域

常用的cdo函数包括：
1.maskregion，myregion ifile ofile
2.sellonlatbox，lon1，lon2，lat1，lat2 ifile ofile
3.两者区别在于，maskregion函数，原文件范围不变，mask以内区域保留原值，mask之外设为缺失值；
sellonlatbox函数，原文件范围缩小为截取的矩形范围，给定矩形四角经纬度值后，只保留矩形内的值。
4.maskregion后跟随的myregion是一个ascii文件，如果剪裁经度120E到90W，纬度20N到20S，那么myregion可以是：
120 20
120 -20
270 -20
270 20

提取任意边界区域

如果所提取的边界范围不是矩形的，而是不规则区域，可以用maskregion 实现。
步骤如下：
1.利用不规则区域shp图，提取边界点；
可以参考