- 积分
- 590
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-9-13
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
假如你有一个二维水平分布的格点数据,每个值都对应地图上的某个经纬度坐标,当你想把这组数据画在地图背景上时,用ncl该怎么做呢?
首先需要做的是告诉ncl 你的数据的经纬度信息。 比如数据维度为 data(/M,N/),那么就要用attribute和units来指定它的经纬度信息:
data!0 = "lat" ; 指定第一个维度的名称。
data!1 = "lon" ; 指定第二个维度的名称。
data&lat = fspan(lat_start, lat_end, lat_num) ; 指定纵坐标分别对应的纬度,这里只是个例子,你可以指定你的纬度值。
data&lon = fspan(lon_start, lon_end, lon_num) ; 指定横坐标分别对应的经度。
data&lat@units = "degrees_north" ; 赋予lat方向信息,这个是必须有的。
data&lon@units = "degrees_east" ; 赋予lon方向信息,这个是必须有的。
最后需要注意的是,在画图时,应给resource设置如下属性:
res@gsnAddCyclic = False
这是告诉ncl 数据的东西两端不闭合,不用为形成经圈而弥补缺口了。默认的是True,会造成数据投射错误。ř
Searching...
|
评分
-
查看全部评分
|