ncl 使用shape文件图

zhangqh

我在用shp 文件的时候，shapefile_mask_data知道这个函数变量必须是二维，lat 和lon
然后下面是脚本。用了一下循环
lt   = new((/10,120/),"float")
do i=0,119
    do j=0,9
f     = addfile("/TPlake.clm2.h0.2003-2012.nc","r")
tlak =f->TLAKE
tlake=tlak-273.15
copy_VarCoords(tlake,tlak)
shp_f            = "lake.shp"
mask_clm     := shapefile_mask_data(tlake(i,j,:,:),shp_f,True)
printVarSummary(mask_clm)
tl                  = wgt_areaave_Wrap(mask_clm,1.0,1.0,0)
lt(j,i)             = tl
delete(tl)
delete(mask_clm)
delete(tlake)
delete(tlak)
end do
end do
报错提示： (0) shapefile_mask_data: Error: not a valid rectilinear, curvilinear, or unstructured grid

画了好几次，也改过shp文件，换过读入的数据。都是一样的错误。不知道为什么 。
各位帮忙看看这个错误吧！！

xhw0734

我在用shapefile_mask_data的时候也报了这个错误，一样的报错，现在找到原因了么？是啥原因？

xhw0734

xhw0734 发表于 2016-9-9 16:12
我在用shapefile_mask_data的时候也报了这个错误，一样的报错，现在找到原因了么？是啥原因？

原因已找到，有兴趣可以看一下shapefile_utils.ncl的注解

zhangqh

  谢谢

zongqing123

学习了，谢谢

subtropical

xhw0734 发表于 2016-9-9 18:09
原因已找到，有兴趣可以看一下shapefile_utils.ncl的注解

您好 这个shapefile的ncl脚本可以分享一下吗谢谢