请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 4213|回复: 5

给txt文件加上经纬度信息

[复制链接]

新浪微博达人勋

发表于 2017-1-11 20:31:09 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x

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

脚本:
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”。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-1-16 10:51:41 | 显示全部楼层
不好意思,我也弄过
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-1-17 14:00:18 | 显示全部楼层
txt格式是不行,模仿NETCDF没得用,你非得使用txt,建议先读成一维数据,再转二维,再给坐标属性
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-1-18 20:29:36 | 显示全部楼层
权威 发表于 2017-1-17 14:00
txt格式是不行,模仿NETCDF没得用,你非得使用txt,建议先读成一维数据,再转二维,再给坐标属性

谢谢您的回复,我试试!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-1-19 16:54:47 | 显示全部楼层
txt,不是netcdf数据
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-1-19 17:05:09 | 显示全部楼层
f1=asciiread("E:/NCL/a/ri1number.txt",-1,"float")
dim = dimsizes(fi)

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

data = onedtond(f1,(/71,144/))
print(data);打印数据信息
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表