博文
地形数据网格化方法
||
最近做项目，需要把散乱的地形数据(经度，纬度，高度)插值网格化。
https://blog.sciencenet.cn/blog-873750-667762.html
上一篇：第一次被拒稿后的感想
下一篇：读《地震概论》笔记
不会使用GIS软件，只是去年学习过使用GMT软件画图，知道里面有些命令可以gridding data
翻出教程来看了看，实验了一下，确实可用，记在这里。
由于原始数据的分布比较均匀，大概计算了一下精确度，基本上每10minutes可以有一个地形控制点
所以第一步使用blockmedian命令如下：
blockmedian file.xyz -R119/126/31/41 -I10m -V>file10x10.xyz
把原文件中的数据的求中值，然后用surface命令如下：
surface file10x10.xyz -Gfile2x2.grd -I2m -R119/126/31/41 -T0.35
其中-T选项的参数0.35，是因为地形数据的曲率通常取0.35
reference：
Wessel P. and W. H. F. Smith. 2010. The Generic Mapping Tools Technical Reference and Cookbook(Version 4.5.3).
https://blog.sciencenet.cn/blog-873750-667762.html
上一篇：第一次被拒稿后的感想
下一篇：读《地震概论》笔记
