快速地再網路上找尋一陣子,沒有看到很明確的轉換公式?不知是不是已有一些常用的 library ?暫時只找到中研院 GIS 小組,該網站上呈列出眾多開發的小工具程式 - GIS應用支援工具集,其中就看到了兩樣有關:

WGS84_TM2

Web版坐標轉換程式 - http://webgis.sinica.edu.tw/cgi-tran/webtrans.htm

但前者都是 Windows 程式,甚至連 Web 版呼叫的 CGI 都是 exe 檔案,一時之間,就先用 Web 版頂著用 XD 過幾天再來找找或寫一下 Python 版本好了

Python:

def WebQueryGPSToGIS(lat,lon):
        x = None 
        y = None 
        err = None 
        if lat <> None and lon <> None:
                try: 
                        raw_data = urllib.urlopen( 'http://webgis.sinica.edu.tw/cgi-tran/wgstrans.exe?module=5&Lat='+str(lat)+'&Lot='+str(lon) ).read()
                        m = re.search( r"([0-9\.]+),([0-9\.]+)" , raw_data )
                        x = m.group(1)
                        y = m.group(2)
                except Exception, e:
                        err = str(e)
        return { 'x': x, 'y': y, 'lat': lat, 'lon' : lon, }
        #return { 'x': x, 'y': y, 'lat': lat, 'lon' : lon, 'e': err }

後來有找到 Taiwan datums - OSGeo Wiki,上頭有提到 Perl - Geography-NationalGrid-TW-0.08,有空再來改寫成 Python 好了。


changyy 發表在 痞客邦 PIXNET 留言(0) 人氣()