轉貼：清大彭明輝的部落格: 學術文獻回顧與分析的程序與技巧
清大彭明輝的部落格: 學術文獻回顧與分析的程序與技巧: "學術文獻回顧與分析的程序與技巧"
匯入TWD67圖層但卻不能套疊
在使用QGIS中，如果匯入TWD67圖層，也確定選擇了TWD67 / TM2 zone 121 的坐標系統(EPSG: 3828)，但是匯入的圖層卻不能跟TWD97或WGS84等圖層套疊成功，有平移的狀況。 原因可能是，原始的圖檔不小心在指定TWD67坐標系統時，採用了原點位於地心(geocentric)的狀況。而事實上，TWD67坐標系統當初在定義的時候卻不是原點位於地心的，當時採用的是最貼合台灣的locally best-fit橢球。所以修正的方法如下： 定義一個TWD67 Geocentric的坐標系統： 在QGIS中新增坐標系統 名稱：TWD67 Geocentric Format選擇Proj String 在Parameters輸入下列定義(連成一行不空格不換行)： +proj=tmerc +lat_0=0 +lon_0=121 +k=0.9999 +x_0=250000 +y_0=0 +ellps=GRS67 + towgs84=-752,-358,-179,-.0000011698,.0000018398,.0000009822,.00002329 +units=m +no_defs 注意 towgs84定義的7個參數，修正了原點的平儀、旋轉、尺度等，這組參數，反應了TWD67系統定義的時候採用GRS67橢球並非原點位於地心的事實。 最後將圖層坐標系統改成這個新的坐標系統，應該就可以正確套合了 將這個圖層另存到新圖層，儲存的時候選定為正確的TWD67 / TM2 zone 121 的坐標系統(EPSG: 3828)，TWD97 / TM2 zone 121 (EPSG: 3826)或WGS84 EPSG: 4326)，以後就可以正常使用。 2021補充： +proj=tmerc 為Transverse Mercator的代稱，代表用以公尺為單位的縱坐標、橫坐標 如果是經緯度，+proj=longlat，如下(連成一行不空格不換行)： +proj=longlat +ellps=GRS67 +no_defs +towgs84=-752,-358,-179,-.0000011698,.0000018398,.0000009822,.00002329 參考資料： GIS TECH
留言