請問有什麼程式可以批次轉換虎子山TM2至WGS84經緯度呢?
若沒有的話,我想自己試著寫寫看這樣的程式,不曉得有沒有
類似的程式碼可用,或者那裡可以找到相關的技術參考資訊?
參考:
http://jidanni.org/geo/taiwan_datums/index.html
(
嗯! 虎子山跟TWD67是不一樣的.
http://www.keepon.com.tw/htmldata/2003/07/Message3266.html
發表:cygnet
Hu-tzu-shan 使用 International 1924參考橢球體
69年那一版是 TWD67 使用 GRS67 參考橢球體
)
TWD97 的 參考球體, 跟 WGS84 差不多,
(http://jidanni.org/geo/taiwan_datums/pshung.html
差幾公分到數十公分,
因此, 未來使用 TWD97 的地圖與 GPS 定位的坐標, 在導航及登山的需求上應可視為一致.
)
而 GPS 誤差是數公尺, 可以把 TWD97 當 WGS84.
轉換方法,
先把 TWD67 轉成 TWD97,
再把 TWD97 TM2 轉成參考球體的經緯度, 就當成是 WGS84 的經緯度.
TWD67 轉 TWD97 的方程式(誤差應在2M內):
"
另 Also john at thl.ncku.edu.tw Aug 2003 捐的 contributed:
A = 0.00001549
B = 0.000006521
X67 = X97 - 807.8 - A * X97 - B * Y97
Y67 = Y97 + 248.6 - A * Y97 - B * X97
X97 = X67 + 807.8 + A * X67 + B * Y67
Y97 = Y67 - 248.6 + A * Y67 + B * X67
「以四參數法(平移加旋轉因子)來統計, 此算法比七參數法來得簡單! 結果誤差應在2公尺內!」
"
轉換成球體經緯度座標,
利用
./proj -I +proj=tmerc +ellps=GRS80 +lon_0=121 +x_0=250000 +k=0.9999
就可以轉成 WGS84 經緯度.
proj 程式有 source code, 可以在
http://www.remotesensing.org/proj/
下載, 有 unix makefile, win32 makefile 等,
我是在 freebsd 測的...
我把 TWD67 轉 TWD97 的座標, 寫個小程式 t67_to_t97,
(306000,2771000), (311000,2767000) 是兩個 twd67 座標,
而轉出來的
(121d33'47.555"E 25d2'39.277"N), (121d36'45.29"E 25d0'28.569"N) 是 WGS84 座標.
./t67_to_t97 306000 2771000 > t97
cat t97 | ./proj -I +proj=tmerc +ellps=GRS80 +lon_0=121 +x_0=250000 +k=0.9999
121d33'47.555"E 25d2'39.277"N
./t67_to_t97 311000 2767000 > t97
cat t97 | ./proj -I +proj=tmerc +ellps=GRS80 +lon_0=121 +x_0=250000 +k=0.9999
121d36'45.29"E 25d0'28.569"N
<!--引言開始 (您引用了全部文章,請適度引言)-->
sm91kimo wrote:<!--引言結束 -->
請問有什麼程式可以批次轉換虎子山TM2至WGS84經緯度呢?
若沒有的話,我想自己試著寫寫看這樣的程式,不曉得有沒有
類似的程式碼可用,或者那裡可以找到相關的技術參考資訊?
ec 於 2003-09-01 10:49 修改過
原來 proj 這麼複雜,在freebsd make後立即打消port至PPC的
念頭 :p (若硬要port的話可能用perlce的celib將整個UNIX-like
發展環境移過來)
cygent前輩曾提到虎子山和TWD67有近30公尺的差距,待會我再用
MMP/VistaT/PPG/MEM/GoGoMap/FormosaGuide等軟體查查看。
請問ec兄的t67_to_t97就是用四參數法嗎?
另偷懶的請問一下 proj 可否直接輸出120.312546E,23.306712N這樣的格式?
sm91kimo 於 2003-09-01 11:02 修改過
> 念頭 :p (若硬要port的話可能用perlce的celib將整個UNIX-like
> 發展環境移過來)
對啊! 看起來是很麻煩...
我現在是在想,
能不能用 vc++.net 的 managed c++, build 出來,
知道可不可以直接拿到 ppc 上 run :)
.................
看來是不行....
> cygent前輩曾提到虎子山和TWD67有近30公尺的差距,待會我再用
> MMP/VistaT/PPG/MEM/GGM等軟體查查看。
啊! 我弄錯了...
虎子山跟TWD67 是不一樣的...
> 請問ec兄的t67_to_t97就是用四參數法嗎?
是的.
> 另偷懶的請問一下 proj 可否直接輸出120.312546E,23.306712N這樣的格式?
不知道,
可能要看一下 man page 吧!
ec 於 2003-09-01 14:28 修改過
TRG-pro 兄
您也來了!歡迎,歡迎!
請參考一下這些…
WGS84坐標與TM2度分帶坐標轉換公式
http://www.colorado.edu/geography/gcraft/notes/datum/edlist.html
http://www.colorado.edu/geography/gcraft/notes/datum/datum.html
http://www.ccit.edu.tw/~ccchang/datm_PE.pdf
http://www.environment.sa.gov.au/mapland/sicom/sicom/tp_scs.html#tag6
可惜『汽車導航 DIY』一文已失連,是值得參考的好文章!
另外GPSfun 有些轉換程式連結,可以參考一下!
一般而言,如果你有GPS手持機或 MapSource 之類的程式,
它的坐標轉換是使五參數,也就是說除了其上提供的各種 DATUM 外,
可由使用者自定 user datum 、user grid ,
你高興的話,甚至可自定一個以你家為基準點的坐標顯示格式,
以後所有資料全是以你家為方位角及距離…
In the GPS World, everyday is fun.
> 虎子山跟TWD67 是不一樣的...
........
http://www.ccit.edu.tw/~ccchang/datm_PE.pdf
http://www.google.com.tw/search?q=cache:MuUOfUNojNcJ:www.ccit.edu.tw/~ccchang/datm_PE.pdf+%EF%BC%A8u-tzu-shan+twd67&hl=zh-TW&ie=UTF-8
This set of coordinates used in Taiwan was called the Hu-Tzu-Shan coordinate system,
also named the TWD67 (Taiwan Geodetic Datum based on the GRS67).
http://jidanni.org/geo/taiwan_datums/
(note: "Hu-tzu-shan datum" is Garmin's name for TWD67)
http://www.ncc.up.pt/gpsman/gpsmanhtml/manual/html/GPSMandoc_23.html
This grid is usually employed with either the "Hu-Tzu-Shan" datum (also known as "TWD67")
這些似乎說 TWD67 就是 虎子山(Hu-Tzu-Shan).
http://etrex.webz.cz/datumindex.html
Hu-Tzu-Shan International 1924 -637 -549 -203 Taiwan
應該就是指 Hu-Tzu-Shan 是用 International 1924參考橢球體
http://ngis.moi.gov.tw/ngis-net/document2/6/news/26/car.html
汽車導航系統DIY
交通部運輸研究所 史習平
這篇?
<!--引言開始 (您引用了全部文章,請適度引言)-->
cygnet wrote:<!--引言結束 -->
可惜『汽車導航 DIY』一文已失連,是值得參考的好文章!
複製連結




























































































