https://wiki.osgeo.org/w/index.php?action=history&feed=atom&title=Taiwan_datums%2Fcad2twd67
Taiwan datums/cad2twd67 - Revision history
2026-04-19T09:48:34Z
Revision history for this page on the wiki
MediaWiki 1.35.9
https://wiki.osgeo.org/w/index.php?title=Taiwan_datums/cad2twd67&diff=41913&oldid=prev
Jidanni: establish new file
2009-10-15T08:16:56Z
<p>establish new file</p>
<p><b>New page</b></p><div><pre><br />
#!/usr/bin/perl<br />
# cad2twd67 -- # 地籍座標轉換成TWD67<br />
# Taiwan "Cadastral" to TWD67 coordinate transformation<br />
# Copyright : http://www.fsf.org/copyleft/gpl.html<br />
# Created On : Fri May 13 23:28:16 2005<br />
<br />
# 地籍座標即日據時代沿用至今地籍圖之常見座標。<br />
# 單位為「間」，等於 20/11 米，即六台尺。<br />
# 原點為台中公園的 89 號三角點。<br />
<br />
# 只是抄李弘洲﹑朱杏修﹑曾德福 在<br />
# 1999年第十八屆測量學術及應用研討會論文集之公式。<br />
use strict;<br />
use warnings;<br />
<br />
# 「六十九年成果」之二度分帶坐標：<br />
my $XCtm69 = 227361.634 + 0.0;<br />
my $YCtm69 = 2632574.582 + 0.0;<br />
<br />
# 大約本程式是針對全省。<br />
# 故若只針對單一鄉鎮不遠處，<br />
# 亦已有一已知點，<br />
# 可以試細調以上的 0.0，<br />
# 否則誤差可達幾十公尺。<br />
<br />
# 轉換參數<br />
my $XCcad = 5750;<br />
my $YCcad = -21300;<br />
my $A = 1.8182516286522;<br />
my $B = -0.004167109289753;<br />
<br />
# 每行輸入之地籍坐標（單位為「間」）<br />
my $Xcad;<br />
my $Ycad;<br />
<br />
while (<>) {<br />
die "Enter Xcad, Ycad" if ( ( $Xcad, $Ycad, undef ) = split ) != 3;<br />
<br />
# $Xtmtrn, $Ytmtrn 是由地籍坐標($Xcad, $Ycad)經四參數轉換套合至<br />
#「六十九年成果」而成之近似二度分帶坐標。<br />
my $Xtmtrn = $A * ( $Xcad - $XCcad ) - $B * ( $Ycad - $YCcad ) + $XCtm69;<br />
my $Ytmtrn = $B * ( $Xcad - $XCcad ) + $A * ( $Ycad - $YCcad ) + $YCtm69;<br />
printf "%.3f\t %.3f\n", $Xtmtrn, $Ytmtrn;<br />
}<br />
<br />
# 以上 "69" 大概類似 TWD67 我猜。<br />
# 轉回去？ 要深讀他們的文章...<br />
<br />
# Usage example:<br />
# $ more cadfile<br />
# -5000 1250<br />
# -14100 -23500<br />
# $ perl cad2twd67 < cadfile<br />
# 207909.397 2673620.953<br />
# 191260.172 2628657.146<br />
</pre><br />
<br />
[[Category:Taiwan]]</div>
Jidanni