给txt文件加上经纬度信息

小宝19950118

文本文件ri1number.txt是2.5*2.5的数据，维度为(/71,144/),用ncl画图时出现下面的警告，我想是因为我的数据没有经纬度的信息。想请教各位如何加上经纬度的信息？


脚本：
begin
f1=asciiread("E:/NCL/a/ri1number.txt",(/71,144/),"float")

wks=gsn_open_wks("pdf","aa")
res=True
res@cnFillOn=True
res@lbOrientation="Vertical"

plot=gsn_csm_contour_map(wks,f1,res)

end

我试着加上
f1!0 = "lat" ;第一维的名称
f1&lat =fspan(-87.5,87.5,71) ;第一维的值
;lat@units = "degrees_north"

f1!1 = "lon" ;第二维的名称
f1&lon =fspan(0,357.5,144)  ;第二维的值
;lon@units = "degrees_east"

但是还是会出现警告，我在网上找了一些解决的办法，但是他们都是从现成的nc文件中提取经纬度的信息，就是脚本中我用红色标注的部分写成“f1&lat =v&lat”。

tianqh

不好意思，我也弄过

权威

txt格式是不行，模仿NETCDF没得用，你非得使用txt,建议先读成一维数据，再转二维，再给坐标属性

小宝19950118

权威 发表于 2017-1-17 14:00
txt格式是不行，模仿NETCDF没得用，你非得使用txt,建议先读成一维数据，再转二维，再给坐标属性

谢谢您的回复，我试试！

权威

txt,不是netcdf数据

权威

f1=asciiread("E:/NCL/a/ri1number.txt",-1,"float")
dim = dimsizes(fi)

print(dim(0:30)) ;打印开始31字节

data = onedtond(f1,(/71,144/))
print(data);打印数据信息