ncl这么写了，怎么还会出现地理坐标不对的？求指点

海豚∩_∩浩

f="00.grd"
g=fbindirread(f,-1,(/73,50/),"float")
wks  = gsn_open_wks("pdf","wind3")
  


lon           = fspan (73.,145.,73)
lon@units     ="degrees_east"


lat           = fspan (4.,53.,50)
lat@units     ="degrees_north "



g!0           ="longitude"
g!1           ="latitude"
g&longitude   = lon
g&latitude    = lat




check_for_y_lat_coord: Warning: Data either does not contain a valid latitude coordinate array or doesn't contain one at all.
(0)     A valid latitude coordinate array should have a 'units' attribute equal to one of the following values:
(0)         'degrees_north' 'degrees-north' 'degree_north' 'degrees north' 'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN' 'deg north'
(0)     check_for_lon_coord: Warning: Data either does not contain a valid longitude coordinate array or doesn't contain one at all.
(0)     A valid longitude coordinate array should have a 'units' attribute equal to one of the following values:
(0)         'degrees_east' 'degrees-east' 'degree_east' 'degrees east' 'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE' 'deg east'

freekiller

你绘图的变量没有地理坐标属性

海豚∩_∩浩

freekiller 发表于 2012-11-18 18:28

                               
登录/注册后可看大图

你绘图的变量没有地理坐标属性

没有地理坐标属性，不影响在等值线图中插入地图吧？

angel1006

感觉很高深的问题

海豚∩_∩浩

海豚∩_∩浩 发表于 2012-11-18 19:05

                               
登录/注册后可看大图

没有地理坐标属性，不影响在等值线图中插入地图吧？

f="00.grd"
g=fbindirread(f,-1,(/73,50/),"float")
wks  = gsn_open_wks("pdf","wind3")
  

lon           = fspan (73.,145.,73)
lon@units     ="degrees_east"

lat           = fspan (4.,53.,50)
lat@units     ="degrees_north "


g!0           ="longitude"
g!1           ="latitude"
g&longitude   = lon
g&latitude    = lat


为什么这么写了，还是会出现这个错误啊？求指点

hzwjy

海豚∩_∩浩 发表于 2012-11-19 21:10

                               
登录/注册后可看大图

f="00.grd"
g=fbindirread(f,-1,(/73,50/),"float")
wks  = gsn_open_wks("pdf","wind3")

变量g的维数写反了。g(lat, lon) 才对。从你代码看，你第一步读入数据的时候也是写反了。
NCL中多维数组的用法和C类似。

下面附加经纬度信息应该是：
g!0 = "lat"
g!1 = "lon"

好久不见!

您好  请问您如何解决的