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

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 260|回复: 3

[作图] ncl为什么画图出来像扫描一样?

[复制链接]

新浪微博达人勋

发表于 2020-3-19 13:35:37 | 显示全部楼层 |阅读模式

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

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

x
begin
data = fbindirread("/home/ramsey/Documents/r/rr11.grd", 0, (/61,71/), "float")

;------设定经度坐标---------
lon = ispan(70, 140, 1)
lon!0 = "lon"
lon@long_name = "longtitude"
lon@units = "degrees_east"
;-------设定纬度坐标--------
lat = ispan(0, 60, 1)
lat!0 = "lat"
lat@long_name = "latitude"
lat@units = "degrees_north"
;-----设置变量元数据-------
data!0 = "lat"
data!1 = "lon"
data&lat = lat
data&lon = lon
data@long_name = "Correlation Coefficient"

;------plot------
wks = gsn_open_wks("png", "rr1")

res = True
res@gsnAddCyclic = False
res@mpMinLatF = 0
res@mpMaxLatF = 60
res@mpMinLonF = 70
res@mpMaxLonF = 140

;res@cnLevelSelectionMode = "ManualLevels"
;res@cnMinLevelValF = -1
;res@cnMaxLevelValF = 1
;res@cnLevelSpacingF = 0.1

res@cnLinesOn = False
res@cnFillOn = True
res@cnFillMode = "RasterFill"
res@cnRasterSmoothingOn = True

res@tiMainString = "~F9~ Correlation Coefficient to Nanjing"
res@tiMainOffsetYF = -0.035
res@gsnLeftString = ""

res@pmTickMarkDisplayMode = "Always"
res@lbBoxLinesOn = False
plot = gsn_csm_contour_map(wks, data, res)
end
经常出现这样的问题,数据应该是没问题的,我感觉是fbindirread这个函数的读取问题?数据就是71*61的二维数组

rr1.png
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2020-3-19 13:45:31 | 显示全部楼层
有可能是数组大小反了
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2020-3-19 15:29:36 | 显示全部楼层
轩-瑾 发表于 2020-3-19 13:45
有可能是数组大小反了

谢谢,确实是这个问题我把这个和fbinrecread函数搞混了,以为fbindirread也是要按右边维变化
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2020-3-19 15:32:07 | 显示全部楼层
问题解决,但是不知道为什么按(71*61)读入,NCL还是报错,可能涉及到NCL右边维变化快的原理?所以用了transpose转置后变为(61*71)成功画图。
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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