Wei_1144 wrote:
分頁二 以 差異距離(相同名稱基點之間的最大距離) 與 基點名稱 將資料作群組排列...(恕刪)
後面的差異距離怎麼會變成 398,812.34 m ? 真的差了這麼大嗎? 是比對三筆資料中的最大與最小的數據出來的差距嗎?
個人積分：339分
文章編號：5350561
個人積分：5分
文章編號：5360658
Wei_1144 wrote:
分頁二 以 差異距離(相同名稱基點之間的最大距離) 與 基點名稱 將資料作群組排列
我下載了您的excel檔看了一下，有一點問題想請教，就是這裡面的"距離差異" 是從別的地方計算後再貼進來的嗎?
我舉個例：
第386列：(距離差異 點名 內X 內Y 7777橫 7777縱)
39.07 卑南社 262598 2519398 261770 2519606
我用成大的轉換公式把7777的TWD67換成TWD97 TM2：
結果為：(262598 2519398)
取四捨五入到TM2的米，和內政部的X,Y並沒有差異。
不知是否有何處我誤解了?
Wei_1144 wrote:
小花的資料，我取得時已經是WGS84DMS，應該已經做過轉換，轉換工具與準確度不明
小花是用Garmin的Hu-Tzu-Shan/UTM 的3+2參數轉換的
個人積分：540分
文章編號：5364702
差異距離是把 7777的 (TWD67TM2_X, TWD67TM2_Y)、內政部的 (TWD97TM2_X, TWD97TM2_Y) 都轉成 TWD97DMS，然後再用DMS的測距公式去算距離，Miller問為什麼有的基點差異距離會是398公里，是因為我是計算同名基點之間的最大距離，差了398公里是因為7777資料中同名基點實在很多
關於 卑南社，我用Minstrelsy的程式把7777的TWD67TM2轉成TWD97TM2，結果是:
262600.45, 2519400.23 ，
差異距離約是3公尺左右，不是如我Excel中的29公尺，可能是我用的DMS測距公式有問題
但此Excel主要是為了把不同來源的資料做比較，並以差異大小與名稱作群組排列，期望能幫Minstrelsy省下一些整理原始資料的工夫
DMS測距公式的問題，有時間再debug
個人積分：540分
文章編號：5364888
再出一份只含7777與內政部的資料: 比較資料.zip
卑南社的差異距離是 2.97公尺
個人積分：339分
文章編號：5365395
個人積分：5分
文章編號：5369229
個人積分：5分
文章編號：5369390
ptX1: 第一點的經度
ptY1: 第一點的緯度
ptX2: 第二點的經度
ptY2: 第二點的緯度
double getSurfaceDistance(double ptX1, double ptY1, double ptX2, double ptY2)
{
const double dEquatorialRadius = 6378137.0; // TWD97, Semi-major axis of ellipsoid in meters
const double dFlattening = 298.257222101; // TWD97, Inverse flattening, Equator/(Equator-Polar)
double dX;
double dY1, dY2;
double dVal, dVal1, dVal2;
double dValue, dDist;
double dSinX, dCosX, dSinY1, dSinY2, dCosY1, dCosY2, dSinValue, dCosValue;
double dDegtoRad = 0.01745329251994329572; // pi/180
dX = (ptX2 - ptX1)*dDegtoRad;
dY1 = ptY1*dDegtoRad;
dY2 = ptY2*dDegtoRad;
dSinX = sin(dX);
dCosX = cos(dX);
dSinY1 = sin(dY1);
dSinY2 = sin(dY2);
dCosY1 = cos(dY1);
dCosY2 = cos(dY2);
dVal = dSinY1*dSinY2+dCosY1*dCosY2*dCosX;
if(abs(dVal)>1.0)
{
return 0.0f;
}
dValue = acos(dVal);
dSinValue = sin(dValue);
dCosValue = cos(dValue);
dVal1 = dSinY1+dSinY2;
dVal2 = dSinY1-dSinY2;
dVal = (((3*dSinValue-dValue)*dVal1*dVal1/(1+dCosValue)) - ((3*dSinValue+dValue)*dVal2*dVal2)/(1-dCosValue));
dDist = dValue*dEquatorialRadius + (dEquatorialRadius/(4*dFlattening))*dVal;
return dDist;
}
個人積分：540分
文章編號：5374527
minstrelsy wrote:
附上我自己用的測距程...(恕刪)
我是用這個: http://www.spatialtech.com.tw/pic/DGPS/DIST-C.JPG
看它的結構，若兩點距離近 應該誤差不大，若兩點距離遠 可能就會有較大的誤差 (我沒有仔細研究，純猜測)
曾經把出遊的航跡檔用這個公式計算航跡的長度，和Garmin算出來的長度約有 20公尺/69公里(差異/總長) 的差異
個人積分：5分
文章編號：5374802
Wei_1144 wrote:
剛剛看了一下，DMS...(恕刪)
了解了，謝謝您的說明和整理，辛苦了
昨天我看到了積丹尼先生也有三角點的整理，稍微看了一下，TWD67和7777有小數點三位的座標是一樣的，如此來看，此資料應該就是官方測量的資料，如7777說明所述。
millerliu wrote:
真是讓人興奮的消...(恕刪)
這86年內政部七參數，我試驗了一下，在台北的四十張，7777TWD67和內政部資料，差別在1米內，但台南的大坪山和牛息山，差別在4~5米左右。
以我的理解，七參數轉換，在有網形變形的座標系統中，不是每個點都可以達到小誤差的轉換，要用更多參數的多項式多次方回歸公式，才能作到滿足最多點的微小誤差。
個人積分：5分
文章編號：5375679
millerliu wrote:
4. 虎子山datum的XYZ參數是多少呢? 目前我蒐集到的是 -637, -547, -203 和 -634, -549, -201兩個版本, 因為沒有早期座標資料驗證, 無法知道哪一組是誤差小的?
參考本串第88樓
在登山補給站7777的留言板，小花和他的討論中有Hu-tzu-shan UTM的座標，對比同一點的7777的TWD67座標，我把它們都轉成TWD97 DMS：
7777 TWD67：(309681.949, 2667403.999) ==> TWD97DMS (121.595255, 24.108781)
小花 UTM：(356397.12, 2667230.37)
所用da,df皆不變 (International 1924橢球)
小花的UTM轉換後和7777相比：
參數一：(-637,-547,-203) ==> (121.595229, 24.108799) 差別 (0.000026, 0.000018)
參數二：(-634,-549,-201) ==> (121.595215, 24.108828) 差別 (0.000040, 0.000047)
參數三：(-637,-549,-203) ==> (121.595240, 24.108806) 差別 (0.000015, 0.000025)
網上這份資料
寫的是和參數三相同
我的筆記上寫Garmin的UTM用的是參數二，但我忘了記下出處
要想拿到更多UTM的座標來驗証，可能得和小花要
提供參考
