ncl使用shapefile_utils.ncl出现gc_inout error 解决方法

Adiligirl

shapefile_utils.ncl 解决方法
http://www.ncl.ucar.edu/Applications/Scripts/shapefile_utils.ncl
cnmap来自
http://bbs.06climate.com/forum.p ... &fromuid=119514
1、问题
作图需要mask中国区域以外的数据，在使用
  cnmap_shp = "./NCL-Chinamap-master/cnmap/cnmap.shp"
  maskopt             = True
  maskopt@return_mask  = True
;注意使用2D数据（lat,lon)
  cnmapLand_mask = shapefile_mask_data(precip_ano_rate(0,:,:) ,cnmap_shp,maskopt)

会出现 fatal:gc_inout:the rightmost dimension of lat/lon must be at least three.
2、查错
设置  maskopt@debug       = True
地图数据无错
在arcmap中打开cnmap.shp, 察觉应该是地图九段线导致gc_inout函数运行出错
3、解决方法
copy cnmap.shp ,paste and rename it as cnmapLand.shp
删除九段线，重新运行，成功！
4、总结
当出现此类错误时，可检查是否是shp文件线段不闭合，即只有两个端点。

Adiligirl

还有一个bug，mask得到的区域只能有中国大陆的，海南岛和台湾都没加进来

angue666

你好，非常感谢您的分享，刚好遇到这个问题。请问您在解决方法中提到的“删除九段线”，具体是怎么做的呢？或者可以分享您删除后的shp文件吗？

Adiligirl

angue666 发表于 2020-7-16 18:12
你好，非常感谢您的分享，刚好遇到这个问题。请问您在解决方法中提到的“删除九段线”，具体是怎么做的呢？ ...

你好，我是在arcgis里面删除的。后面发现这样做只能提取大陆地区的数据，我就没有使用这个方法了。

angue666

哦，这样啊。非常感谢您的回复，那您还有其他方法来屏蔽中国以外的数据吗？

Adiligirl

angue666 发表于 2020-7-16 22:22
哦，这样啊。非常感谢您的回复，那您还有其他方法来屏蔽中国以外的数据吗？

我现在直接使用
res@mpAreaMaskingOn            = True
res@mpMaskAreaSpecifiers       = (/"China"/)
res@mpOutlineSpecifiers        = (/"China","China:Provinces"/)

angue666

Adiligirl 发表于 2020-7-17 16:42
我现在直接使用
res@mpAreaMaskingOn            = True
res@mpMaskAreaSpecifiers       = (/"China"/ ...

好的，非常感谢您！