ncl站点资料做EOF

qxlong

我想请问，ncl使用站点资料来做EOF分析，思路步骤应该是这怎样？ncl小白有点懵，毫无头绪，请各位大佬指点，就是简单说一下最开始应该做什么，然后是要做什么之类的，这问题困惑一段时间了，查了网站也不太明白，拜托各位了。

MgLily

qxlong 发表于 2019-10-24 11:22
得到evecv的结果之后，是可以直接画图，不需要插值了是吗？还有就是，后面我画图出现了这些问题，可以帮 ...

不可以的，得到的结果是站点数据，你要是想用gsn_csm_contour_map函数的话，需要的数组是格点经纬度的二维数组。

;---------------------使用natgrid进行插值---------------------------
latN = 55.                          ;中国区域（你要插值的范围）
latS = 15.
lonW = 70.
lonE = 140.

dlt = 0.5                           ;插值格点精度
latnum =floattoint(180/dlt)+1       ;纬度格点数
lonnum =floattoint(360/dlt)         ;经度格点数
latitude =latGlobeF(latnum,"lat","latitude","degrees_north")
longitude=lonGlobeF(lonnum,"lon","longitude","degrees_east")

lat=latitude({latS:latN})    ;截取插值范围格点
lon=longitude({lonW:lonE})

nlat=dimsizes(lat)
nlon=dimsizes(lon)

griddata=new((/nlat,nlon/),"float")   

griddata(:,:,:)=natgrid(lats,lons,evecv(:,:),lat,lon) ;站点数据插值,lats,lons是要插值的站点的经纬度信息，你得到的是griddata(模态,lat,lon)这样的数组
griddata!0="lat"
griddata&lat=lat
griddata!1="lon"
griddata&lon=lon

MgLily

思路的话，先分析，从EOF分析的这个函数eofunc_Wrap (x,neval,False)出发，将你想分析的数据变成这个函数需要的data形式，即时间维总是在最右边，data(sta | :,time | :)，设置好参数，用函数出结果evecv(evn|:,sta|:)，然后再考虑怎么表示这个结果（画图）。

MDAmike

一样的呀，都是用eofunc_Wrap函数

MgLily

这个是官网的站点EOF例子http://www.ncl.ucar.edu/Applications/Scripts/eof_0.ncl
EOF分析的函数一般用eofunc_Wrap，你可以去官网看看这个函数的解释
以及论坛上也有很多例子的，可以搜搜看

qxlong

MgLily 发表于 2019-10-22 12:47
思路的话，先分析，从EOF分析的这个函数eofunc_Wrap (x,neval,False)出发，将你想分析的数据变成这个函数需 ...

好的好的，谢谢你

qxlong

MDAmike 发表于 2019-10-21 20:58
一样的呀，都是用eofunc_Wrap函数

好的，我先试试

qxlong

MgLily 发表于 2019-10-22 12:47
思路的话，先分析，从EOF分析的这个函数eofunc_Wrap (x,neval,False)出发，将你想分析的数据变成这个函数需 ...

得到evecv的结果之后，是可以直接画图，不需要插值了是吗？还有就是，后面我画图出现了这些问题，可以帮我看看吗？
图一是报错，图二是printVarSummary（evecv）的结果，图三图四是部分程序img]file:///C:\Users\Administrator\AppData\Roaming\Tencent\Users\1095312682\QQ\WinTemp\RichOle\I9NO2@QITD@{D`}LDM[MPXV.png[/img][img]file:///C:\Users\Administrator\AppData\Roaming\Tencent\Users\1095312682\QQ\WinTemp\RichOle\B%2}N45ENL80F3)RH09Y5@N.png[/img]file:///C:\Users\Administrator\AppData\Roaming\Tencent\Users\1095312682\QQ\WinTemp\RichOle\82I7WQBU()XUKG%ARGXU4QS.png

ANNA0217

第一个要设置变量维度名称和属性

小羊趴在云朵上

感谢感谢，非常有用