读取txt文件并且画图的问题

墨青衫

用的是风云3b的MWRI数据，在做完主成分分析之后输出到了txt文本里面。图画出来经纬度出现问题。求大神指点一下
下面是程序
txt文件是第一列是lat 第二列是lon 后面5列是数据 只画第三列的。
begin
fi = "d:/data/try.txt"
nrow = numAsciiRow(fi)
ncol = numAsciiCol(fi)
data = asciiread(fi, (/nrow, ncol/), "float")
;print(data)
lat=data(:,0)
lon=data(:,1)
RI=data(:,2)
;data@_FillValue = -999
wks = gsn_open_wks("ps","d:/RI")
res = True ;
res@sfXArray         = lon
res@sfYArray         = lat
;res@xyLineColor = "Blue" ; line color
res@tiMainString = "RI"
res@gsnMaximize = True
res@cnFillOn = True
res@cnLinesOn = False
res@cnLineLabelsOn = False
res@cnInfoLabelOn = False
res@mpMinLatF = 10
res@mpMaxLatF = 70
res@mpMinLonF = -90
res@mpMaxLonF = -30
res@pmTickMarkDisplayMode = "Always"
res@gsnSpreadColors = True
res@cnMonoLineColor = False
res@cnLevelSelectionMode = "ExplicitLevels"
res@cnLevels = (/-3,-2,-1,0,1,2/)
res@gsnMaximize = False
plot = gsn_csm_contour_map(wks,RI,res)
end
这个是画出来的图。。。

MeteoInfo

为什么你会觉得有问题呢？

墨青衫

图形是重复的啊，而且左边边缘很明显的一个不对的拉伸。

MeteoInfo

墨青衫 发表于 2015-11-6 16:56
图形是重复的啊，而且左边边缘很明显的一个不对的拉伸。

Swath的形状是很完整的，应该不是经纬度的问题，如果你觉得图形有问题应该是数据本身的问题。

南国的佳人

不过大神可不可以指点一下数据是在哪里下载到的呢？

墨青衫

MeteoInfo 发表于 2015-11-6 17:03
Swath的形状是很完整的，应该不是经纬度的问题，如果你觉得图形有问题应该是数据本身的问题。

国家卫星中心

墨青衫

墨青衫 发表于 2015-11-6 16:56
图形是重复的啊，而且左边边缘很明显的一个不对的拉伸。

已解决，由于是卫星数据，所以应该用2d的坐标
data@lat2d=lat
data@lon2d=lon
坐标是从原始下载数据里面读取的。

慧仔333

这个数据是格点的吧？请问ncl做站点等值线填充图有吗？