透過ChatGPT協助修正QGIS中TWD67_TM2 CRS設定錯誤

QGIS軟體是利用Proj.4（地圖投影函式庫）來進行坐標系統轉換參考，在QGIS 1.9版本以上已經內建臺灣地區常用的EPSG代碼，然而目前Proj.4 492版中對於EPSG:3828（即TWD67 TM2 中央經線121度，適用於臺灣本島）以及EPSG:3827（即TWD67 TM2 中央經線119度，適用於澎湖群島）的定義不夠完整，因此會造成TWD67 TM2 <-> TWD97 TM2坐標轉換會產生問題。目前仍有部分政府機關（例如：地政事務所、林務局等單位）在政府開放資料平台上提供GIS資料時，仍會提供TWD67 TM2坐標系統的GIS圖層，因此這個問題仍然困擾很多人。

我們可以透過ChatGPT，直接提示” 請提供 橫麥卡托二度分帶坐標系統（TM2）以 TWD67 測地基準、中央經線為 121 度 的 Proj 4.0 格式定義 ” ，就會得到正確 Proj String格式定義（含完整七參數）：

+proj=tmerc +lat_0=0 +lon_0=121 +k=0.9999 +x_0=250000 +y_0=0 +ellps=GRS67 +towgs84=-752,-358,-179,-1.1698e-6,1.8398e-6,9.822e-7,2.329e-5 +units=m +no_defs

不過，目前最新版的QGIS軟體也不再支援 以PROJ 4.0 格式來自訂CRS，需改用WKT格式，在ChatGPT，進一步提示” 請輸出 WKT 格式 “，就會得到 WKT 格式 （符合ESRI 定義）：

PROJCS["TWD67_TM2_121",
    GEOGCS["GCS_TWD67",
        DATUM["D_Taiwan_Datum_1967",
            SPHEROID["GRS_1967",6378160,298.247167427,
                AUTHORITY["EPSG","7036"]],
            TOWGS84[-752,-358,-179,-0.0000011698,0.0000018398,0.0000009822,0.00002329]],
        PRIMEM["Greenwich",0],
        UNIT["Degree",0.0174532925199433]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["False_Easting",250000.0],
    PARAMETER["False_Northing",0.0],
    PARAMETER["Central_Meridian",121.0],
    PARAMETER["Scale_Factor",0.9999],
    PARAMETER["Latitude_Of_Origin",0.0],
    UNIT["Meter",1.0],
    AUTHORITY["EPSG","3828"]]

就可以在QGIS軟體中，[設定]->[自訂坐標參考系統]，新增一個 TWD67_TM2_121 坐標系統定義（CRS）。

亦可以直接將WKT內容儲存為副檔名為 .prj 的純文字檔，名稱要與.shp, .shx, .dbf相同，且放在相同的資料夾中，日後即可自動套用坐標系統，無須再重新指定。

TWD67 TM2，中央經線119度的WKT格式如下：

PROJCS["TWD67_TM2_119",
    GEOGCS["GCS_TWD67",
        DATUM["D_Taiwan_Datum_1967",
            SPHEROID["GRS_1967",6378160,298.247167427,
                AUTHORITY["EPSG","7036"]],
            TOWGS84[-752,-358,-179,-0.0000011698,0.0000018398,0.0000009822,0.00002329]],
        PRIMEM["Greenwich",0],
        UNIT["Degree",0.0174532925199433]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["False_Easting",250000.0],
    PARAMETER["False_Northing",0.0],
    PARAMETER["Central_Meridian",119.0],
    PARAMETER["Scale_Factor",0.9999],
    PARAMETER["Latitude_Of_Origin",0.0],
    UNIT["Meter",1.0],
    AUTHORITY["EPSG","3827"]]

參考資料：