[閒聊] 如何用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