画任意剖面问题

weiguangjian

手中有某个地区的DEM数据（即为海拔的栅格数据），已转成grads能画图的数据格式，并已做出平面图。但是这个数据只有x，y维是变化的，z维就是海拔高度值，t维固定。想要做地形的剖面是不可能实现的。现在就想要沿任意直线做一剖面，纵坐标表示高度值就可以，请问各路高人，该如何实现。...................

'set x 1'
'set y 1'
'set t 1'
lon1 = 120.03
lon2 = 120.13
lat1 = 35.91
lat2 = 35.99
lon = lon1
'collect 1 free'
while (lon <= lon2)
  lat = lat1 + (lat2-lat1)*(lon-lon1) / (lon2-lon1)
  'collect 1 gr2stn(h,'lon','lat')'
  lon = lon + 0.0083333
endwhile

.............
'd coll2gr(1,-u)'
................
按以上脚本执行提示如下
Collect Error: Z or T must be the only varying dimension
Error from COLL2GR:  Collection 1 empty
Operation Error:  Error from coll2gr function
  Error ocurred at column 1
DISPLAY error:  Invalid expression
  Expression = coll2gr(1,-u)
沿任意经线或纬线的这种剖面图倒是可以画出来

kongfeng0824

没用grads画过剖面图，知道surfer可以。但是批量画的画估计还是grads

善人/jw

http://bbs.06climate.com/forum.php?mod=viewthread&tid=10966LZ可以参考一下这个帖子

weiguangjian

我也没用过surfer，这个软件功能似乎更强大。

weiguangjian

kongfeng0824 发表于 2013-1-22 22:55

                               
登录/注册后可看大图

没用grads画过剖面图，知道surfer可以。但是批量画的画估计还是grads

我也没用过surfer，这个软件功能似乎更强大。

weiguangjian

善人/jw 发表于 2013-1-23 00:16

                               
登录/注册后可看大图

http://bbs.06climate.com/forum.php?mod=viewthread&tid=10966LZ可以参考一下这个帖子

我就是参看的这个帖子，由于我的数据的z维没有变化，只有高度值，通过上面的方法想把高度值插值到任意的直线上，但是就出现如上提示的错误信息。
沿经线或纬线的"剖面"能画出来，如下

kongfeng0824

weiguangjian 发表于 2013-1-23 15:18

                               
登录/注册后可看大图

我也没用过surfer，这个软件功能似乎更强大。

嗯是的。只要用好一款软件就可以的。

流浪的心

DEM数据

scwjx520

麻烦问下楼主，怎么得到一个省的分辨率为100米的DEM数据？希望能详细告诉我，谢谢！

weiguangjian

scwjx520 发表于 2014-2-7 20:53

                               
登录/注册后可看大图

麻烦问下楼主，怎么得到一个省的分辨率为100米的DEM数据？希望能详细告诉我，谢谢！

某地的高分辨率地形数据一般都是保密的，这只能向相关部门求解，我也没办法。