millerliu wrote:
金馬地區座標轉換好了, 請注意, 礙於軟體限制計算經緯度的關係, 馬祖地區座標的北緯一律改為25/30/00來計算TW97, TW67, 特別注意了!! 金門地區則是精確的沒有問題!
轉換座標後才發現原來北緯都是25.5度的座標, 居然連TW97/67 XY值的Y是不一樣的!
請問您是怎麼驗証林務局轉換出來的金門的TWD67是精確的? 有其他的TWD67座標資料來源嗎?
如我之前所問的,金門馬組有多少POI是用TWD67的呢?
因為現在是想取得POI的WGS84座標,要是手上有一群舊datum座標資料(不管是虎子山或TWD67),但並沒有WGS84座標,我們就可以藉由找到轉換方式來轉換這群資料。
但要是現在並沒有這一群資料,找轉換方式也就沒有用了。
轉換的計算X,Y值和舊的X,Y值都有相關,所以雖然您設的舊的Y都一樣,但舊的X不同,新的X,Y也都會不同
512筆真的太為難我了....我今天還要因為release BIOS工作....加班到十點以後呢....
我提供這三十二顆資料散落台灣各地,盡量是東南西北都有分佈....也盡量取測量等級高的
也為了怕都是軍測原點所以裡面有一部分是地籍原點....希望能有這少少三十二筆資料瞭解
黃皮書裡虎子山/經緯度跟虎子山/UTM的關聯
整本黃皮書還蠻厚的....整理你表列所有點我可能要六日假期都用上去才有辦法....但講真的
原點跟森林點部份如我發表GDB,已經精進南北座標數據到半米以內....對於現在手持式gramin
機子來說已經足夠了,如果你是要探討日據時代留下來測量資料,轉換的問題,這三十二筆資料應該
可以表達地籍原點跟軍測原點的一些轉換的問題
另外黃皮書裡除了原點跟森林點之外,還有我不熟悉的很多種點位....最起碼我昨晚還發現有陸補點位
基本上這些點我都不瞭解出處,因為日本時代根本沒這些點位....我不知道真實黃皮書記載這類型點位
準確度是如何....也沒花時間找這種點位.
如果是要研究日治時代圖根點,森林點的測量數據轉換問題...我覺得倒是個可以研究的目標,講到這裡....
也還是老調重彈.....玩基石的這群人已經算是爬山界少數人,這群人當中有極高比例只會關心日據時代
古蹟基石,譬如說原點,森林點,圖根點,水準點等.....對於國民政府立的基石....真的沒啥興趣,我相信這點
就聯涉獵基石甚廣的張大鵬大哥也是一樣的情況.用古老日劇座標資料找到沒人找到的基石....這才是我當
初整理資料最初的用意...哈哈....好像太自私了
PS.金門跟馬祖沒有日據時代基石.
powerspot wrote:
512筆真的太為難我了....我提供這三十二顆資料散落台灣各地,儘量是東南西北都有分佈....也儘量取測量等級高的也為了怕都是軍測原點所以裡面有一部分是地籍原點....希望能有這少少三十二筆資料瞭解黃皮書裡虎子山/經緯度跟虎子山/UTM的關聯
如果是要研究日治時代圖根點,森林點的測量數據轉換問題...我覺得倒是個可以研究的目標,...對於國民政府立的基石....真的沒啥興趣,我相信這點就聯涉獵基石甚廣的張大鵬大哥也是一樣的情況.用古老日劇座標資料找到沒人找到的基石....這才是我當初整理資料最初的用意...哈哈....(恕刪)
其實你填入我指定的32個三角點就可以啦, 我解釋為何是這32個指定點, 因為座標轉換最基本資料是共同點, 而且這共同點必須有新舊座標資料, 目前這些都是內政部衛星測量點(WGS84), TWD67座標, 同時還是一等三角點(精確性最高), 加上又都有虎子山系統與Bessel系統座標, 這樣就可以推算出GPS需要的轉換參數DX, DY, DZ. 你提供的若干點位甚至不在512個共同點內, 這樣資料是有應用的困難的!
基石迷無意找現在三角點我可以理解, 但是計算座標轉換參數時卻必須貫穿古今, 因此目前我想探究與改進的是:
1. Bessel 1841系統, 虎子山系統的最精細轉換參數為何?
2. 7777資料中的TWD67 TM2的計算來源有無轉換產生的誤差? (各類基石)
3. 基石資料中台灣澎湖地區的.000基石校正回應有的定點.
4. 澎湖地區所有基石座標如何精確定位?
hwsu315 wrote:
有,我有把s190n...(恕刪)
不知道您的程式運作正常了嗎? 你前次計算澎湖轉換參數是用36點(我不知道是哪些?), 可否有空時使用下面提供的共同點座標來計算澎湖地區的TWD67, Hu-Tzu-Shan系統轉換三參數? 驗證後發現該區域三角點幾乎沒有什麼三角點移動的問題, 不超過20cm. 今年又加入了三個共同點, 因此TWD67共同點有15個, Hu-Tzu-Shan則有7個, 可以直接轉化為WGS84.
小花兄, 這份資料我又修正了若干內容, 不妨下載保留起來!
不同基準共同點檔案: 200804/mobile01-0f5736c92bc85111a878a65b96686333.zip
或者教教我如何用Excel計算Multiple Regression, 因為我不太理解僅有TWD67-X. TWD67-Y, TWD97-X, TWD97-Y的二維座標卻能算出三維的變數DX, DY, DZ? 如果學通了也許就能獨自發掘各種基準的轉換參數.
millerliu wrote:
不知道您的程式運作正常了嗎? 你前次計算澎湖轉換參數是用36點(我不知道是哪些?)
還是一樣的結果,您林務局的程式,是在"大地基準轉換"的Tag那兒轉換澎湖的TWD97-TWD67嗎? 我再看的結果好像只有此處有轉換的功能,有其他人也和您拿了林務局轉換程式嗎?
我相信您所用的林務局程式在澎湖地區有使用上澎湖基準點來計算的參數,因為您列的這幾個轉換後的座標,和當時TWD67的座標相差很小,但奇怪的就是我的轉換結果和您的就是不同。
我前次計算澎湖的Abridged Moledensky轉換參數 是用內政部裡有的12個點,另24點是從7777實測的gdb檔中取得的(業餘手持GPSr誤差)
millerliu wrote:
或者教教我如何用Excel計算Multiple Regression, 因為我不太理解僅有TWD67-X. TWD67-Y, TWD97-X, TWD97-Y的二維座標卻能算出三維的變數DX, DY, DZ? 如果學通了也許就能獨自發掘各種基準的轉換參數.
Excel似乎沒有算Multiple Regression的功能 (或我沒找到),要在Excel中算,可能需要找相關的VB plug-in或自己寫計算。
我是用Origin軟體,我也試用過SPSS軟體,它也有此功能,我相信其他的統計軟體應該都有此功能,之前我也找到過線上的Multiple Regression計算,我未記下網址,您可以google一下。
要注意的是,在把數值代入Multiple Regression計算前,需要把經緯度代入Abridged Molodensky公式,取得多項式的形式,才可來計算。
而在參考點不多時(兩三百以下),套裝軟體的Multiple Regression會算不出正確的三參數,需要先預設某一個參數,再算出其他兩個;此後再微調預設的參數,求出此範圍最佳的三參數。
點多時,三個參數可直接由套裝軟體的Multiple Regression算出來,Garmin/Mapsource中所用的Datum參數我就是用這方法反推得到的。
以上是我計算的方式
ps.
澎湖的虎子山,印象中,我在算了澎湖的TWD67後有稍微算一下,平均誤差好像不太好,有七八公尺,不確定,我會再算一次確定一下。
hwsu315 wrote:
還是一樣的結果,您林務局的程式,是在"大地基準轉換"的Tag那兒轉換澎湖的TWD97-TWD67嗎? 我再看的結果好像只有此處有轉換的功能,有其他人也和您拿了林務局轉換程式嗎?...(恕刪)
Excel似乎沒有算Multiple Regression的功能 (或我沒找到),要在Excel中算,可能需要找相關的VB plug-in或自己寫計算。
我是用Origin軟體,我也試用過SPSS軟體,它也有此功能,我相信其他的統計軟體應該都有此功能,之前我也找到過線上的Multiple Regression計算,我未記下網址,您可以google一下。
要注意的是,在把數值代入Multiple Regression計算前,需要把經緯度代入Abridged Molodensky公式,取得多項式的形式,才可來計算。
而在參考點不多時(兩三百以下),套裝軟體的Multiple Regression會算不出正確的三參數,需要先預設某一個參數,再算出其他兩個;此後再微調預設的參數,求出此範圍最佳的三參數。
只有你用了該程式, 另外7777大大早有此軟體, 因此你是唯一的user. 的確就在大地基準轉換那個label項下, 真是頭痛, 有軟體卻不能用!
Excel算複回歸是在書局看到在統計上的應用時提到的, 還有實際例子, 只是你提出要轉為多項式形式計算, 我就不理解了! 是要計算出 http://webhelp.esri.com/arcgisdesktop/9.2/index.cfm?TopicName=Equation-based_methods
最下方公式中的哪些數字?
如果無法參透這些計算方式, 看來只有回頭研讀 http://mail.thhs.tpc.edu.tw/~jueping/personal/paper/4.pdf提到的嘗試錯誤法來找出三參數.
參考點居然要兩三百以下來計算, 果真是統計運算, 那對參考點已屬有限的澎湖地區來說真是天方夜譚了, 難怪你要另加7777資料. 上述論文還可以找到中正理工的教授寫七參數計算程式, 不知道可否
請託開教授幫忙寫三參數計算程式? 目前想不到該程式如何運算的, 難道只要輸入參考點的經緯度即可? 如果這樣就可以找出Bessel, Hu-Tzu-Shan的參數了.
關閉廣告




























































































