山友分享的全國基石及衛星控制點資料(98.11.18更新檔案)

Wei_1144 wrote:
我是用這個: http://www.spatialtech.com.tw/pic/DGPS/DIST-C.JPG...(恕刪)

hwsu315 wrote:
所用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用的是參數二，但我忘了記下出處

Wei_1144 wrote:
再出一份只含7777與內政部的資料: 比較資料.zip

Wei_1144 wrote:
表中沒有給定內政部資...(恕刪)

hwsu315 wrote:
另外，我加上了七參數的轉換比較，七參數是用86年內政部的平差工作報告的參數
我只取了成大轉換誤差小的五百多個點來計算，再加上平均的差別量
也許米兄可以算出平均差別更小的七參數

minstrelsy wrote:
我用內政部的七參數略...(恕刪)

lon = *x*DEG_RAD;
lat = *y*DEG_RAD;
height = *z*DEG_RAD;

if(lat<-HALF_PI&&lat>-1.001*HALF_PI)
lat = -HALF_PI;
else if(lat>HALF_PI&&lat<1.001*HALF_PI)
lat = HALF_PI;
else if((lat<-HALF_PI)||(lat>HALF_PI))
return;

if(lon>PI)
lon -= (2*PI);

sin_lat = sin(lat);
cos_lat = cos(lat);
r = TWD67_A/(sqrt(1.0-TWD67_ECC*sin_lat*sin_lat));
x1 = (r+height)*cos_lat*cos(lon);
y1 = (r+height)*cos_lat*sin(lon);
z1 = ((r*(1-TWD67_ECC))+height)*sin_lat;

x2 = x1+TWD67_DX+TWD67_S*(lon+TWD67_RZ*lat-TWD67_RY*height);
y2 = y1+TWD67_DY+TWD67_S*(-TWD67_RZ*lon+lat+TWD67_RX*height);
z2 = z1+TWD67_DZ+TWD67_S*(TWD67_RY*lon-TWD67_RX*lat+height);

hwsu315 wrote:
但第三行為何要乘上 DEG_RAD 呢?

hwsu315 wrote:
我所用的公式(Bursa-Wolf method/ Helmert Transformation)在此連結

hwsu315 wrote:
上述的差別，我把您的公式x2,y2,z2的式子改寫到我的Excel裡，分別用您blog的參數和86年內政部的參數計算，

hwsu315 wrote:
奇了, 這樣跟我算出來的數據還是不同.

以第一筆資料松嶺為例, 我用內政部參數轉出來的內政部TWD67_TM2分別是 (272068.232828, 2676476.019918),