用CSV格式画站点图，怎么变换经纬度

指甲钳

用CSV格式的数据画图，第二列时经度，第三列是纬度，第四例时数据，我需要画站点的降雨量，之前是用下面的方式转换经纬度的，能够画出图，但是现在我需要设置颜色。用了cnlevels和CNFILLCOLORS后，就画不出图了。我想问一下时什么原因。一下时代码：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
wks = gsn_open_wks("x11","r20")     ; Open a workstation and.
gsn_define_colormap(wks,"precip_11lev") ; 设置色表


;--------------data-----------------------------
fpath = "e:\data\r20.csv"

lines = asciiread(fpath,-1,"string")

delim = ";"

  lon  = tofloat(str_get_field(lines,1,delim))
  lat  = tofloat(str_get_field(lines,2,delim))
  data = tofloat(str_get_field(lines,3,delim))
data@_FillValue = -9999

res = True

res@gsnMaximize   = True
res@gsnDraw       = False
res@gsnFrame      = False

res@sfXArray = lon
res@sfYArray = lat

res@trYMinF   = 22.75
res@trYMaxF   = 23.35
res@trXMinF   = 112.8
res@trXMaxF   = 113.3
;设置等值线参数

res@cnFillOn                = True
res@cnLinesOn               = True           ; turn off contour lines
res@gsnSpreadColors         = True           ; use full range of colormap
res@cnSmoothingOn           = True
res@cnSmoothingDistanceF    = 0.001
res@cnSmoothingTensionF     = -2.5
;res@cnLevels  = (/0.1,9.9,24.9,49.9,99.9,249.9,250.0/)
;res@cnFillColors= (/2,4,6,8,11,12,13/)
map = gsn_csm_contour(wks,data,res)

;-------shp
filename="gis/nhreverse/nhall.shp"

resp                  = True
resp@gsFillColor = (/"White"/)
polygon = gsn_add_shapefile_polygons(wks,map,filename,resp)


lnfilename="gis/nanhai/nanhai.shp"
lnres                  = True
lnres@gsEdgesOn = True
lnres@gsEdgeColor      = "white"
lnres@gsLineThicknessF = 1.0

polyline = gsn_add_shapefile_polylines(wks,map,lnfilename,lnres)


draw(map)
frame(wks)

end

这是之前的

这是之后的

指甲钳

没人知道咩？

指甲钳

我已经帮你改正了。

zsumore

粗心大意。。。。。。。。。。。。。。。

wsdeqx

顶顶顶顶顶

kongfeng0824

不错的探究     

kolrsong

try to set res@cnLevelSelectionMode="ExplicitLevels"

931404656

能不能插值呢？