NCL如何提取特定不规则区域的数据？

Sunshining

我想问一下啊，我怎样从NCEP/NCAR下载的再分析资料（nc格式）中只提取出特定区域（比如大西洋）的数据来做其他的计算？
大洋是不规则的，用经纬网格来限制行不通啊。。
我从官网上找到一张类似的图，看了一下它的脚本，屏蔽其他区域一段的代码是这样写的：
; Mask out all regions but that sprcified by the user =>'region'
;=================================================;
   rmask     = f[iStrt]->REGION_MASK       ; read region info from 1st file
  ;printVarSummary(rmask)                  ; (nlat,mlon)   ... 2D
   rmask_3D  = conform(x, rmask, (/1,2/) ) ; broadcast the mask to 3D
   x         = mask(x, rmask_3D.ne.region, False)     ; data for region  
   delete ( rmask_3D )                     ; not needed
；================================================；
里面region=6。
我不太懂这里rmask具体是什么，我在我的代码里，就把变量的lat和lon（全球范围）读到了rmask里面，其余一样，但是不能画图，提示警告是区域中没有数据，我把False换成True也不行。
有没有大神可以指教一下，这里这个rmask变量里边是什么啊？是全球的经纬网格？还是什么意思？

Sunshining

warning:ContourPlotInitialize: no valid values in scalar field; ContourPlot not possible:[errno=1101]
这是弹出来的警告

香草拿铁

请问官网的链接可否贴一下？
另，请问问题解决了么

zhouxiaosky

多谢分享~~~~

Sunshining

香草拿铁 发表于 2017-3-31 22:33
请问官网的链接可否贴一下？
另，请问问题解决了么

http://www.ncl.ucar.edu/Applications/eof.shtml
这是官网的链接，里面的第三个例子。问题还没有解决啊啊啊啊
帮帮我

香草拿铁

Sunshining 发表于 2017-4-6 17:55
http://www.ncl.ucar.edu/Applications/eof.shtml
这是官网的链接，里面的第三个例子。问题还没有解决啊 ...

region=6
我看到这个了
初学NCL。说下自己的见解。我觉得这里的region编号是和另外一处的region得来的。出处这里：
https://www.ncl.ucar.edu/Applications/popmask.shtml
[size=13.3333px]popmask_2.ncl
这个6指的是大西洋。
ncl的部分代码为：
  printVarSummary(rmask)  printMinMax(rmask, True)       ; -14 ==> 11  do REGION=-14,11   
n = num(rmask.eq.REGION)   
  print("REGION="+REGION+"   n="+n)  
end do

我想我能理解楼主的意图，应该和我的也差不多。
我现在的想法是利用NCL读取shapefile的功能，进行区域的popmask来判定某一特定区域。对该区域内的数据进行计算处理
共勉

mon

香草拿铁 发表于 2017-4-6 22:00
region=6
我看到这个了
初学NCL。说下自己的见解。我觉得这里的region编号是和另外一处的region得来的 ...

对啊，我也是想用shp文件来选特定区域的数据，比如中国的一个省，但不知道具体怎么搞，请问大牛已经找到办法了吗

流浪的心

学习了，谢谢！

mwy021

“region=X”代表特定海域，它的NC数据比你多一个维度，这个维度就是各个海域

Sunshining

哇原来是这个样子啊，这种数据哪里可以下载呀