TWD67與TWD97


TWD67與TWD97恐怕是台灣使用與設計GPS軟硬體都會碰到的兩種座標系統. 依據內政部地政司衛星測量中心的說明, TWD67與TWD97的定義如下:


一、TWD67:





(一) 、參考橢球體採用1967年新國際地球原子如下:



1.長半徑: a =6378160公尺



2.短半徑: b =6356774.7192公尺



3.扁率: f =(a-b)/a=1/298.25





(二)、 大地基準點以南投埔里之虎子山起算:



經度λ=120° 58′ 25. 975〞



緯度φ=23° 58′ 32. 340〞



對頭拒山之方位角α=323° 57′ 23. 135〞





(三)、高程基準面:



台灣本島以基隆平均海水面起算。



澎湖以馬公平均海水面起算。





(四)、地圖投影:


     有關地籍測量及大比例尺測圖所應用之坐標系統,係採用橫麥卡托投影經差二度分帶,台灣本島之中央子午線為121。,坐標原點為中央子午線與與赤道交點,且橫坐標西移250,000公尺,中央子午線之尺度比率為0.9999。




 二、TWD97(TAIWAN DATUM 97,1997台灣大地基準)係民國八十六年內政部採用國際地球參考框架而定義更適合於全國大地坐標系統之新基準,其內容如下:






(一)、新 國家坐標系統之名稱命名為1997台灣大地基準(TWD97),其建構係採用國際地球參考框架(International Terrestrial Reference Frame簡稱為ITRF)。ITRF為利用全球測站網之觀測資料成果推算所得之地心坐標系統,其方位採國際時間局(Bureau International de l`Heure,簡稱為BIH)定義在1984.0時刻之方位。





(二)、新國家坐標系統之參考橢球體採用1980年國際大地測量學與地球物理學協會(International Union of Geodesy and geophysics 簡稱為IUGG)公布之參考橢球體 (GRS80),其橢球參數如下:
 長半徑 a=6378137公尺
 扁率 f=1/298.257222101





(三)、高程基準面:


  內政部已完成台灣一等水準網,計二、○六五個一等水準點測量工作,並於基隆設置水準原點及副點,高程系統以基隆港平均海水面為高程基準面,據此訂定二○○一年臺灣高程基準(簡稱TWVD2001),作為台灣高程測量控制系統之基準。






(四)、台 灣、琉球嶼、綠島、蘭嶼及龜山島等地區之投影方式採用橫麥卡托投影經差二度分帶,其中央子午線為東經121 度,投影原點向西平移250,000 公尺,中央子午線尺度比為0.9999;另澎湖、金門及馬祖等地區之投影方式,亦採用橫麥卡托投影經差二度分帶,其中央子午線定於東經119 度,投影原點向西平移250,000公尺,中央子午線尺度比為0.9999。

不 過TWD67與TWD97在座標應用上似乎還是要轉換成橫麥卡托投影的二度分帶. 官方對於其相關轉換的公式與參數亦無多加說明. TWD97是使用GRS80橢球, WGS84亦是使用GRS80橢球, 但據許多人的說法, 似乎TWD97與WGS84雖是極相近, 但座標亦非完全吻合.


意見

匿名 提到...

看過米兄在mobile01參與的討論串
今天下午上班不專心
抓了MTWDCON.exe搭配的
TWDCON.GRD(2003/4/28,1039kB)
用一些hex editor去看了一下
發現了一些可能有趣的資料
其中第一群數字在0x000038~0x0568EB=>88621=401*221=401*17*13個float
(FD 87 4E 44 as 32 bit IEEE single: 826.124816894531)
(5A 05 4F 44 as 32 bit IEEE single: 828.083618164063)
其中第二群數字在0x0568EC~0x0AD19F=>88621=401*221=401*17*13個float
(B5 14 54 C3 as 32 bit IEEE single: -212.08088684082)
(12 6E 4B C3 as 32 bit IEEE single: -203.42996215820)
其中第三群數字在0x0AD1A0~0x103A53=>88621=401*221=401*17*13個float
(26 B3 70 41 as 32 bit IEEE single: 15.043737411499)
(B8 F5 BE 41 as 32 bit IEEE single: 23.869979858398)
另 0x000004~0x000005=>DD 00 as integer: 221~~int nx;/* Number of columns */
另 0x000006~0x000007=>91 01 as integer: 401~~int ny;/* Number of rows */
另 0x000008~0x00000F=>00 00 00 00 00 17 01 41 as 64 bit IEEE double: 140000~~double x_min;/* Minimum x coordinate */
另 0x000010~0x000007=>00 00 00 00 00 F9 15 41 as 64 bit IEEE double: 360000~~double x_max;/* Maximum x coordinate */
另 0x000018~0x00001F=>00 00 00 00 08 63 42 41 as 64 bit IEEE double: 2410000~~double y_min;/* Minimum y coordinate */
另 0x000020~0x000027=>00 00 00 00 48 70 45 41 as 64 bit IEEE double: 2810000~~double y_max;/* Maximum y coordinate */
另 0x000028~0x00002F=>00 00 00 60 17 BE 89 40 as 64 bit IEEE double: 823.761413574219~~double x_inc;/* x increment */
另 0x000030~0x000037=>00 00 00 20 2B 15 8A 40 as 64 bit IEEE double: 834.646057128906~~double y_inc;/* y increment */
不過看起來不太對勁~

匿名 提到...

補充~
米兄在mobile01的討論串內似乎有提到MTWDCONV是用fortran寫的?
這是猜測的?或是米兄有看過原始檔?
如果有看過原始檔?可否惠示讓我們參研一下?
謝謝~

匿名 提到...

WGS84用的橢球模型應該不是GRS80吧,
1/f扁率不一樣