ncl用自己的格点数据画图，怎么赋值经纬度？

huangying

我有一个25*81的格点数据（从matlab处理来的csv格式），数据表示x方向81个格点（93.5E~105.5E，间隔0.125°）；y方向25个格点（36.5N~39.5N，间隔0.125°）。
用ncl画图，出来了图像，但是xy轴现在是显示了0~24和0~80，我想将xy轴改成我要的坐标轴36.5N~39.5N，93.5E~105.5E。不知道如何修改？
我想的解决方法有2个：
1.直接修改xy轴显示的内容
2.给我的变量赋值，再画图

由于我没有用过ncl ，所以以上两个解决方法都没有整出来。花了3天时间，希望能有人帮帮我。


以下是画图脚本：

begin
;=================================================;
  nrows = 25
  ncols = 81
  data  = asciiread("ydata.csv",(/nrows,ncols/),"float")
londata = asciiread("londata.csv",(/1,81/),"float")
latdata = asciiread("latdata.csv",(/1,25/),"float")
;=================================================;
  wks  = gsn_open_wks("x11","conwomap")   ; send graphics to PNG file
  res                     = True          ; plot mods desired
  res@cnFillOn            = True          ; turn on color
  res@cnFillPalette       = "gui_default" ; set color map
  res@vpWidthF            = 0.5           ; change aspect ratio of plot
  res@vpHeightF           = 0.3
  res@gsnMaximize         = True          ; maximize plot size
  plot = gsn_csm_contour(wks,data,res)       ; contour the variable
  end

huangying

已经解决了，给变量经纬度赋值：

begin
;=================================================;
; open file and read in data
;=================================================;
  nrows = 25
  ncols = 81
prec= asciiread("ydata.csv",(/nrows,ncols/),"float")

lon = fspan(93.5, 103.5, 81)
lat = fspan(36.5, 39.5, 25)
;print(lat)
;---Create lat and long coordinate variables
   lon!0          = "lon"
   lon@long_name  = "lon"
   lon@units      = "degrees-east"
   lon&lon        = lon
   lat!0          = "lat"
   lat@long_name  = "lat"
   lat@units      = "degrees_north"
   lat&lat        = lat

;---Name dimensions of t and assign coordinate variables
    prec!0    = "lat"
   prec!1    = "lon"
   prec&lat  = lat
   prec&lon  = lon
   prec@long_name = "precipitation"
   prec@units     = "mm"

;----------------------------------------------------------------------
  wks  = gsn_open_wks("x11","conwomap")   ; send graphics to PNG file
  res                     = True          ; plot mods desired
  res@cnFillOn            = True          ; turn on color
  res@cnFillPalette       = "gui_default" ; set color map
  res@vpWidthF            = 0.5           ; change aspect ratio of plot
  res@vpHeightF           = 0.3
  res@gsnMaximize         = True          ; maximize plot size
  plot = gsn_csm_contour(wks, prec,res)       ; contour the variable
  end

瑶池仙子121

大咖云集，菜鸟学习

杨豆

太棒啦 嘿嘿嘿菜鸟学到了

YQY优优

楼主太棒了

来自苹果APP客户端

ouyang

感谢楼主，学习了

Liujiaxin

这个经纬度太稀疏了，请问有什么解决办法？

小羊趴在云朵上

ncl绘图：感谢楼主