[闲聊] 如何用QGIS拼最新的DEM

楼主: sharkbig (bigshark)   2022-09-19 21:12:21
以下文章只是一些技术分享,如果写得还不够亲民还请见谅。
首先你需要的是一台计算速度还不错快的电脑,然后并在这台电脑中安装QGIS或gdal程式
然后找到你的qgis安装位置,确认是否存在一个叫做bin的资料夹,并且复制那个资料夹
的路径(e.g. C:\Program Files\QGIS 3.16\bin)
step 2:
然后就可以去下载全台湾的20m DEM,并且把他们全部都解压缩,你就会得到一大堆的
grd档跟hdr档。因为台湾的模式是依照县市分类,但是这在应用上非常的麻烦,所以我会
把他们解压缩到同一个资料夹内。
step 3:
在"外一层"的资料夹中开启命令提示字符(通常在直接在路径列输入cmd就会
跳出一个黑色视窗)
输入:
set PATH=C:\Program Files\QGIS 3.16\bin;%PATH%
^^^^ 就是上面提到的安装路径
接下来就是拼图时间:
gdalwarp -t_srs EPSG:3826 -srcnodata -999 解压缩的资料夹/xxxx*.grd \
merge_xxxx.tif
^^^^ ^^^^
因为全台的DEM数目超过了gdalwarp
可以处理的数目,因此我这边采用
五万的编号方式来先做分组,请将xxxx依序替换成:
9623、9723、
9522、9622、9722
9421、9521、9621、9721
9420、9520、9620、9720
9419、9519、9619、
9418、9518、9618
9417、9517、9717
9516
(*上面的编号不含澎湖)
最后用相同的指令
gdalwarp -t_srs EPSG:3826 merge_*.tif TW20mDEM.tif
然后一个小提醒,因为图资释出的时候就是采用TWD97的座标系统,因此不能在过程中
对资料进行座标转换,要不然都会导致格点的偏移,使地形出现不自然的线条!
还有在qgis的接口中有一个网格的工具叫做merge,这个指令是没办法顺利拼图的!
原因在于资料的原始格式不是由下向上作排列
如果还有问题也欢迎提出,至于上面的我都看不懂,我还是很想要新版的DEM的人,就
自己连络我吧XD
作者: apieceofslow (缓慢)   2022-09-25 00:46:00
推推

Links booklink

Contact Us: admin [ a t ] ucptt.com