爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 5812|回复: 2

[作图] 格点数据画图求助

[复制链接]

新浪微博达人勋

发表于 2015-11-30 14:43:04 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 墨青衫 于 2015-11-30 14:46 编辑

QQ截图.jpg 有txt文本的格点,360*720*17的  0.5°*0.5°的数据
程序如下
begin
finame = "d:/222.txt"
nrow = numAsciiRow(finame)
data = asciiread(finame, (/nrow/), "integer")
data@_FillValue = -999
data2=onedtond(data,(/360,720,17/))
data2@_FillValue = -999
data2=data2/100
data2!0="lat"
data2!1="lon"
lat = fspan(-90,90,360)
lon = fspan(0,360,720)
lon@units = "degrees_east"
lat@units = "degrees_north"
data2&lat = lat
data2&lon = lon
wks = gsn_open_wks("ps", "d:/1")
gsn_define_colormap(wks,"WhBlGrYeRe") ; 设置色表
res = True
res@sfXArray=lon
res@sfYArray=lat
res@cnFillOn = True ; 填色
res@cnLinesOn = False ; 去掉等值线
res@cnLineLabelsOn = False ; 去掉等值线的数字
res@cnInfoLabelOn = False ; 去掉右下角说明
res@gsnSpreadColors = True ; 6.1.0 后默认即为 True
res@lbBoxLinesOn = False ; Label bar 上不标记黑线
res@gsnRightString = ""
res@gsnAddCyclic=False
plot = gsn_csm_contour_map(wks, data2(:,:,0), res)
end
求大神指导一下应该怎么做,是哪里出错了呢
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-12-3 08:14:50 | 显示全部楼层
好图  学习了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-12-4 23:18:57 | 显示全部楼层
很强大,可以的
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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