小白求助~~~

gloria394083

刚开始学ncl，现在是相出一张等值线图，各种错误，我把一部分代码写上来，请大家帮帮忙！谢谢了！txt文件里是每个格点值，没有经纬度


begin
lon=ispan(100,145,1)
lon@units="degree_east"
lat=ispan(20,45,1)
lat@units="degree_north"

data=asciiread("./1b.txt",(/26,26/),"float")

wks=gsn_open_wks("png","test")

gsn_define_colormap(wks,"rainbow")

res=True
res@gsnAddCyclic=False
res@mpDataSetName="Earth..4"
res@mpDataBaseVersion="MediumRes"
res@mpOutlineOn=True

res@mpMinLatF=20
res@mpMaxLatF=45
res@mpMinLon=100
res@mpMaxLon=125

res@mpProjection ="LambertConformal"  
res@mpLambertMeridianF =110.0
res@mpLimitMode = "LatLon"
res@mpLambertParallel1F =.001  
res@mpLambertParallel2F =89.999

res@mpAreaMaskingOn =True

res@cnFillOn     = True
res@cnLinesOn    =False      
res@cnLineLabelsOn = False
res@cnFillDrawOrder ="PreDraw"  

res@cnSmoothingDistanceF=True

res@tiMainString="test"

plot=gsn_csm_contour_map(wks,data,res)

end

1649518749

站点资料要插值

gloria394083

1649518749 发表于 2017-1-3 21:05
站点资料要插值

我改成这样：

lon=ispan(100,145,1)
lat=ispan(20,45,1)
data=asciiread("./aaa.txt",(/26,26/),"float")

olon=new(26,"float")
olat=new(26,"float")
data1=new((/26,26/),"float")
do i=0,25
   olon=100+i
   end do
do j=0,25
   olat=20+j
   end do
   
olon!0="lon"
olon@long_name="lon"
olon@units="degree_east"
olon&lon=olon

olat!0="lat"
olat@long_name="lat"
olat@units="degree_north"
olat&lat=olat

data@_FillValun=-9999
rscan =(/10,5,3/)  
data1 =obj_anal_ic_deprecated(lon,lat,data,olon,olat,rscan,False)

它的报错是:Number of dimensions in parameter (2) of (obj_anal_ic_deprecated) is (2), (1) dimensions were expected
fatal:["Execute.c":8578]:Execute: Error occurred at or near line 34 in file 1bb.ncl

gloria394083

请问还要怎么改呢

1649518749

gloria394083 发表于 2017-1-3 22:16
请问还要怎么改呢

olon=fspan(100,145,45)   
olat=fspan(20,45,25)
olon@units="degrees_east"
olat@units="degrees_north"
加上这个试试？

1649518749

正常情况不是用站点的经纬度去插值吗？

gloria394083

1649518749 发表于 2017-1-3 23:03
正常情况不是用站点的经纬度去插值吗？

Number of dimensions in parameter (2) of (obj_anal_ic_deprecated) is (2), (1) dimensions were expected这个报错是什么意思呢？

1649518749

gloria394083 发表于 2017-1-3 23:17
Number of dimensions in parameter (2) of (obj_anal_ic_deprecated) is (2), (1) dimensions were expe ...

这本来是二维的，但现在只是一维吧，大概这意思，插值没成功

beacher

我也是这个问题，请问你怎么解决的呀