經建版地形圖數值資料檔匯入工作

最後編輯：2017-10-01 建立：2016-09-02 歷史紀錄

登山

REX T內政部部國土測繪中心的「經建版地形圖數值資料檔(比例尺為二萬五千分之一、五萬分之一及十萬分之一)」前經「105年行政院資料開放諮詢小組第2次會議」列為甲類資料，並經內政部105年7月26日台內地字第1051306149號令修正發布「國土測繪成果資料收費標準」第2條附表附件2，開放資料供免費下載使用[1]，授權條款採用行政院「政府資料開放授權條款－第1版」[2]

這批 2016/07/28 釋出的圖檔包括二萬五千分之一經建版地形圖計262幅、五萬分之一經建版地形圖計80幅及十萬分之一經建版地形圖計7幅，共計349幅。

雖然不含等高線圖層，但是包含水系、道路、行政界線、鐵道、高壓線、建築區等圖層，及圖例、中文註記等向量圖層。以 QCad/LibreCAD 等開源軟體確認，

其格式為

AC1027 - AutoCAD 2013/2014
TWD67 二度分帶坐標

目前需要將其轉為 DXF，整理出可以匯入的資料後，再轉成 SHP/OSM 格式上傳。如能匯入資料，會大幅提供 OSM 在水系資訊、偏鄉道路的資訊量。由於資料龐大，需要志工來轉檔、整理這批資料。歡迎加入 Trello 工作事項[3]討論。

[1] http://www.moi.gov.tw/chi/chi_news/news_detail.aspx?sn=10850&type_code=01

[2] 政府資料開放授權條款－第1版 | 政府資料開放平臺 - http://data.gov.tw/node/18250

[3] OpenStreetMap Backlog 挖洞填坑的匯入國土測繪中心新增政府開放資料「經建版地形圖數值資料檔」 | Trello - https://trello.com/c/SEPcEDyM/196--

REX T工作流程

REX T1. AC1027 轉檔成開放格式如 SHP
2. 對比現有圖資。

3. 將可用資料轉出成 osm 格式。
4. 人工以 josm 整合 osm 資料。

REX T技術進展

REX T圖檔格式轉換

QCad/LibreCad 可以轉換 AC1027 成 DXF，並從 QGIS 進行後續處理。

HAPPYMAN C轉換好的 DXF 在此: https://goo.gl/eGBv7N

REX T已知問題

QCad 點位資料的編碼錯誤，TEXT 變成亂碼。

座標系統轉換

請參考 psilotum (mutolisp) 的轉換參數

台灣大地座標系統的轉換 « mutolisp 技術筆記 - http://mutolisp.logdown.com/posts/207563-taiwan-geodetic-coordinate-system-conversion

基本地形圖資料庫圖層與定義

圖層編碼與圖示可以參考李萬凱先生的教材 - AutoCAD基礎_10_圖層管理(經建版地形圖) - http://blog.xuite.net/lwkntu/blog/443706365-%2812%29AutoCAD%E5%9F%BA%E7%A4%8E_10_%E5%9C%96%E5%B1%A4%E7%AE%A1%E7%90%86%28%E7%B6%93%E5%BB%BA%E7%89%88%E5%9C%B0%E5%BD%A2%E5%9C%96%29

基本地形圖資料庫專有名詞定義 http://ngis2.moi.gov.tw/Storage/MOI_NGIS/Files/20100825104857644.pdf
基本地形圖資料庫地形資料分類編碼表 http://ngis2.moi.gov.tw/Storage/MOI_NGIS/Files/20100825104618317.pdf

REX T資料爬梳

REX T山頂高度

happyman: 本想從裡面擷取山頭標高的點位資料，但事實上這個數值資料已經不是經建三版的，而且資料雜亂無章，我放棄~~~ 如果看紀錄會需要某某峰(沒山名,只有標高)的資料，還是自己去對照經建版地圖吧。 - https://www.facebook.com/photo.php?fbid=10154799301732269&set=p.10154799301732269&type=3&theater

查「基本地形圖資料庫地形資料分類編碼表」，

獨立標高點 - 98120
標高點註記 - 90360

工作流程

1. 以 QGIS 取出以下兩個圖層

a. "Layer" = '98120' 標高打叉。Layer = '98120' 可以抓出獨立標高點的具體位置。但是他是以兩條 LineString 打叉的中心來標注標高點。

b. "Layer" = '98120_TXT' 標高文字。Layer = '98120_TXT' 可以抓出獨立標高點的文字。但是他的位置是經建版的文字位置，不是山頭的位置。

c. 用以下 CSR 可以正確標出位置

+proj=tmerc +lat_0=0 +lon_0=121 +k=0.9999 +x_0=250000 +y_0=0 +ellps=aust_SA +towgs84=-750.739,-359.515,-180.510,0.00003863,0.00001721,0.00000197,0.99998180 +units=m +no_defs

2. 以 QGIS Line Intersection 演算法取出標高打叉點之交叉點座標

3. 將上述 1,2 之結果存入 SpatiaLite，以交叉點座標查周圍 0.005 的標高文字。將結果輸出為新的 SpatiaLite Table.

4. 透過 org2osm 將資料轉為 OSM 資料格式，可匯入 josm 中檢查。

5. 利用 JOSM Script 檢查 OSM 資料庫中否有重複資料。