任意分辨率三维地形mask制作

气候-谭志强

      最近打算用ERA Interim的0.125*0.125高分辨率资料绘制垂直螺旋度的剖面图，遇到了地形绘制问题，这2天在论坛上找了好久，发现没有满足这种高分辨率的三维mask文件，打算自己手动做一份地形mask文件。
      我是这么考虑的：
      ERA Interim有地面气压资料，我将逐格点地面气压和气压层（1000,925,850,700等）做比较，如果地面气压值小于某一气压层，说明存在地形，即该层在地下，将mask文件中该格点该层赋值为-1，无地形则赋值为1，如果地形在2层之间（如850和700），则该格点该层地形数据应为-1+(850-地面气压)/(850-700)；这样做后又遇到一个问题，即由于数据分辨率太高，导致地形轮廓锯齿，我进行了水平17*17格点平滑，最终出图如下：

沙颖凯

图中的地形是青藏高原东侧吗？考虑到楼主的地形在0.125*0.125的分辨率下也只是这种台阶形状，我感觉至少从结果上楼主的做法并不合理。

我记得论坛上这种东西挺多的，方法也不少，楼主可以多看看。

zihan6688

m文件和其他的是独立的吗？

气候-谭志强

沙颖凯 发表于 2014-11-26 12:04
图中的地形是青藏高原东侧吗？考虑到楼主的地形在0.125*0.125的分辨率下也只是这种台阶形状，我感觉至少从 ...

平滑过得结果

气候-谭志强

zihan6688 发表于 2014-11-26 13:25
m文件和其他的是独立的吗？

不是独立的

mofangbao

楼主是想说分辨率太低所以插值了吗，貌似这个地形分辨率是越高越好哈~另外多谢分享~

气候-谭志强

mofangbao 发表于 2014-11-26 14:08
楼主是想说分辨率太低所以插值了吗，貌似这个地形分辨率是越高越好哈~另外多谢分享~

哇，大神出现了！！不是分辨率太低插值，而是分辨率太高（0.125*0.125），画出的剖面图地形轮廓不平滑，所以滑动平滑处理了

mofangbao

气候-谭志强 发表于 2014-11-26 14:15
哇，大神出现了！！不是分辨率太低插值，而是分辨率太高（0.125*0.125），画出的剖面图地形轮廓不平滑， ...

楼主有请教过师兄师姐或者老师，问问这样做合理吗，美观固然重要，实事可能更重要吧。。。

气候-谭志强

mofangbao 发表于 2014-11-27 08:43
楼主有请教过师兄师姐或者老师，问问这样做合理吗，美观固然重要，实事可能更重要吧。。。

这个木有请教别人，不过觉得平滑到什么程度可以自己试，
D:\Z)URPJMQC%5BMC~6I42BG336N.jpg
9*9平滑
是不是感觉更接近实况

气候-谭志强

气候-谭志强 发表于 2014-11-27 17:07
这个木有请教别人，不过觉得平滑到什么程度可以自己试，

9*9平滑

补图。。。。。。。。。