ncl画图少了一半

刺娓娓

本人是新手，用ncl画图，出错了，就是读取txt格式1324行，3列（经度，纬度，亮温），得到了一半的图咋回事，大神帮忙啊
脚本如下：
;***********************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"


begin
ascii_filename = "/nuist/p/work/***/ncl/yunyu.txt"
argu  = asciiread(ascii_filename,(/1324,3/),"float")

obs = argu(:,2)
lat=argu(:,1)
lon=argu(:,0)
olon = fspan(120,146,256)
olat = fspan(10,30,153)

olon!0          = "lon"
olon@long_name  = "lon"
olon@units      = "degrees-east"
olon&lon        = olon
olat!0          = "lat"
olat@long_name  = "lat"
olat@units      = "degrees_north"
olat&lat        = olat
obs@_FillValue = 999999.000000
rscan = (/1/)   
dataobs = obj_anal_ic_Wrap(lon,lat,obs,olon,olat,rscan,False)
dataobs@units="K"


wks = gsn_open_wks ("png","obs")  
gsn_define_colormap( wks ,"Rainbow")         
res                         = True            
res@gsnAddCyclic            = False
res@tiMainString            ="obs"
res@gsnMaximize             = True
;res@gsnDraw                 = False
;res@gsnFrame                = False

;>--------------------------------------------<
;            set for the map
;>--------------------------------------------<
res@mpMinLatF               = 10.                        
res@mpMaxLatF               = 30.
res@mpMinLonF               = 120.
res@mpMaxLonF               = 146.

res@mpFillOn                = True
res@mpOutlineOn             = True  ; Use outlines from shapefile
res@mpDataBaseVersion       = "MediumRes"
res@mpDataSetName           = "Earth..4"
res@mpOutlineSpecifiers     = (/"China:states","Taiwan"/)
;>--------------------------------------------<
; set for the plot

res@cnFillOn                = True               
res@cnLinesOn               = False            
;res@cnLevelSpacingF         = 1.            
res@gsnSpreadColors         = True         
res@lbLabelAutoStride       = True
res@cnSmoothingOn           = True               
plot=gsn_csm_contour_map(wks,dataobs,res)

end

刺娓娓

画的图如下

刺娓娓

以下是贴的数据

andrewsoong

刺娓娓 发表于 2014-11-25 17:25
画的图如下

我看正常啊~~~极轨卫星数据？

刺娓娓

andrewsoong 发表于 2014-11-25 17:39
我看正常啊~~~极轨卫星数据？

amsu数据,不是的，右半边应该是有数据的，但是不知道为什么，画图只出来一半

andrewsoong

刺娓娓 发表于 2014-11-25 18:37
amsu数据,不是的，右半边应该是有数据的，但是不知道为什么，画图只出来一半

我看你这右边又锯齿状，所以猜测是卫星扫描的边缘。

刺娓娓

我不知道为什么，修改rscan = (/1/) 为rscan = (/10,5,3/)结果图是这样的。 我之前用其他工具画图，右边是有数据的，所以很奇怪不知道怎么回事啊

刺娓娓

andrewsoong 发表于 2014-11-25 19:08
我看你这右边又锯齿状，所以猜测是卫星扫描的边缘。

我不知道为什么，修改rscan = (/1/) 为rscan = (/10,5,3/)结果图是这样的。 我之前用其他工具画图，右边是有数据的，所以很奇怪不知道怎么回事啊

longlivehj

刺娓娓 发表于 2014-11-25 18:37
amsu数据,不是的，右半边应该是有数据的，但是不知道为什么，画图只出来一半

数据的确只有一半！
先不要插值，把obs直接画出来，一看就知道了。

刺娓娓

longlivehj 发表于 2014-11-26 09:50
数据的确只有一半！
先不要插值，把obs直接画出来，一看就知道了。

不好意思，我还是不太懂，要怎么改呢？