ENVI中自定義坐標(biāo)(以北京54和西安80為例

1、ENVI中自定義坐標(biāo)（以北京54和西安80為例）更新至ENVI4.7下本帖最后由 lianyi_dyq 于 2010-10-6 11:56 編輯 ENVI中自定義坐標(biāo)系說明以北京54和西安80為例1 地理投影的基本原理常用到的地圖坐標(biāo)系有2種，即地理坐標(biāo)系和投影坐標(biāo)系。地理坐標(biāo)系是以經(jīng)緯度為單位的地球坐標(biāo)系統(tǒng)，地理坐標(biāo)系中有2個重要部分，即地球橢球體（spheroid）和大地基準(zhǔn)面（datum）。由于地球表面的不規(guī)則性，它不能用數(shù)學(xué)公式來表達(dá)，也就無法實施運算，所以必須找一個形狀和大小都很接近地球的橢球體來代替地球，這個橢球體被稱為地球橢球體，我國常用的橢球體如下：橢球體名稱年代長半軸（米）短

2、半軸（米）扁率WGS8419846378137.06356752.31：298.257克拉索夫斯基（Krasovsky）19406378245.06356863.01：298.3IAG-7519756378140.06356755.31：298.257表1 我國常用橢球體 大地基準(zhǔn)面指目前參考橢球與WGS84參考橢球間的相對位置關(guān)系（3個平移，3個旋轉(zhuǎn)，1個縮放），可以用其中3個、4個或者7個參數(shù)來描述它們之間的關(guān)系，每個橢球體都對應(yīng)一個或多個大地基準(zhǔn)面。 投影坐標(biāo)系是利用一定的數(shù)學(xué)法則把地球表面上的經(jīng)緯線網(wǎng)表示到平面上，屬于平面坐標(biāo)系。數(shù)學(xué)法則指的是投影類型，目前我國普遍采用的是高斯克呂格投

3、影，在英美國家稱為橫軸墨卡托投影（Transverse Mercator）。高斯克呂格投影的中央經(jīng)線和赤道為互相垂直，分帶標(biāo)準(zhǔn)分為3度帶和6度帶。美國編制世界各地軍用地圖和地球資源衛(wèi)星像片所采用的全球橫軸墨卡托投影（UTM）是橫軸墨卡托投影的一種變型。高斯克呂格投影的中央經(jīng)線長度比等于1，UTM投影規(guī)定中央經(jīng)線長度比為0.9996。 我國規(guī)定1：1萬、1：2.5萬、1：5萬、1：10萬、1：25萬、1：50萬比例尺地形圖，均采用高斯克呂格投影。1：2.5至1：50萬比例尺地形圖采用經(jīng)差6度分帶，1：1和1：2.5萬比例尺地形圖采用經(jīng)差3度分帶。2 北京54與西安80坐標(biāo)系 先了解大地坐標(biāo)的概念

4、。 大地坐標(biāo)，在地面上建立一系列相連接的三角形，量取一段精確的距離作為起算邊，在這個邊的兩端點，采用天文觀測的方法確定其點位（經(jīng)度、緯度和方位角）， 用精密測角儀器測定各三角形的角值，根據(jù)起算邊的邊長和點位，就可以推算出其他各點的坐標(biāo)。這樣推算出的坐標(biāo)，稱為大地坐標(biāo)。 我國1954年在北京設(shè)立了大地坐標(biāo)原點，由此計算出來的各大地控制點的坐標(biāo)，稱為1954年北京坐標(biāo)系。為了適應(yīng)大地測量的發(fā)展，我國于1978年采用國際大地測量協(xié)會推薦的IAG-75地球橢球體建立了我國新的大地坐標(biāo)系，并在1986年宣布在陜西省涇陽縣設(shè)立了新的大地坐標(biāo)原點，由此計算出來的各大地控制點坐標(biāo)，稱為1980年大地坐標(biāo)系。

5、我們經(jīng)常給影像投影時用到的北京54或者西安80坐標(biāo)系是投影直角坐標(biāo)系，如下表所示為北京54和西安80坐標(biāo)系采用的主要參數(shù)，坐標(biāo)名稱投影類型橢球體基準(zhǔn)面北京54Gauss Kruger（Transverse Mercator）Krasovsky北京54西安80Gauss Kruger（Transverse Mercator）IAG75西安80表2 我國常用坐標(biāo)系參數(shù)列表 從中可以看到我們通常稱謂的北京54坐標(biāo)系、西安80坐標(biāo)系實際上指的是我國的兩個大地基準(zhǔn)面。3 參數(shù)的獲取 對于地理坐標(biāo)，只需要確定兩個參數(shù)，即橢球體和大地基準(zhǔn)面。對于投影坐標(biāo)，投影類型為Gauss Kruger（Transver

6、se Mercator），除了確定橢球體和大地基準(zhǔn)面外，還需要確定中央經(jīng)線。 大地基準(zhǔn)面的確定關(guān)鍵是確定7個參數(shù)（或者其中幾個參數(shù)），北京54基準(zhǔn)面可以用三個平移參數(shù)來確定，即“-12,-113,-41,0,0,0,0”，很多軟件近似為Krasovsky(0,0,0,0,0,0,0)基準(zhǔn)面；西安80的7參數(shù)比較特殊，各個區(qū)域不一樣。一般有兩個途徑：一是直接從測繪部門獲?。欢歉鶕?jù)三個以上具有西安80坐標(biāo)系與其他坐標(biāo)系的同名點坐標(biāo)值，利用軟件來推算，有一些綠色軟件具有這個功能，如Coord MG。中央經(jīng)線獲取可有以下兩種方法，第一種根據(jù)已知帶號計算，6度帶用6*N-3，3度帶用3*N；第二種方法

7、是根據(jù)經(jīng)度從圖1中查找。下載 (28.33 KB)2008-12-2 17:47圖1高斯克呂格投影的分帶4 ENVI中定義坐標(biāo) ENVI中的坐標(biāo)定義文件存放在HOMEITTIDL70productsenvi45map_proj 文件夾下，三個文件記錄了坐標(biāo)信息： ellipse.txt橢球體參數(shù)文件 datum.txt基準(zhǔn)面參數(shù)文件 map_proj.txt坐標(biāo)系參數(shù)文件 在ENVI中自定義坐標(biāo)系分三步：定義橢球體、基準(zhǔn)面和定義坐標(biāo)參數(shù)第一步、添加橢球體語法為 ,。 這里將“Krasovsky，6378245.0，6356863.0”和“IAG-75，6378140.0，6356755.3”加

8、入ellipse.txt末端。注：ellipse.txt文件中已經(jīng)有了克拉索夫斯基橢球，由于翻譯原因，這里的英文名稱是Krassovsky，為了讓其他軟件平臺識別，這里新建一個Krasovsky橢球體。第二步、添加基準(zhǔn)面 語法為,。 這里將“D_BEIJING_1954, Krasovsky, -12, -113, -41”和“D_Xian_1980,IAG-75,0,0,0”加入datum.txt 末端。注：為了更好的與ArcGIS系列產(chǎn)品兼容，從ENVI 4.7開始，所有產(chǎn)品包括 ENVI、ENVI+IDL、ENVI Zoom 和ENVI Ex，全部采用ArcGIS投影轉(zhuǎn)換引擎（ENVI4

9、.7之前的版本用的是GCTP常規(guī)制圖轉(zhuǎn)換包），對用戶來說，ENVI菜單中所有的投影操作不變，同時還直接支持ArcGIS中的投影類型。但是自定義北京54及西安80坐標(biāo)系有一些改變，即定義兩個坐標(biāo)系的基準(zhǔn)面（datum）時候使用統(tǒng)一的名稱：D_BEIJING_1954和D_Xian_1980。第三步、定義坐標(biāo) 在ENVI任何用到投影坐標(biāo)的功能模塊中都可以新建坐標(biāo)系(在任何地圖投影選擇對話框中，點擊“New”按鈕。)，這里我們選擇Map-Customize Map Projection，如圖1所示，添加一個20度帶（6度）的坐標(biāo)，添加的參數(shù)如圖2所示。 注：投影類型選擇Transverse Merc

10、ator，Scale factor填寫0.9996，與Gauss-Kruger等同。False easting中如果把帶號，即39500000，得到的坐標(biāo)就帶有帶號。 下載 (117.54 KB)2010-9-30 14:15圖2 自定義坐標(biāo)系選擇Projection-Add New Projection, 將投影添加到 ENVI 所用的投影列表中。選擇 File Save Projections，存儲新的或更改過的投影信息。這樣一個新的投影坐標(biāo)就新建完成。 相同方法，添加一個20度帶（6度）的西安80坐標(biāo)系（如圖3所示）。下載 (122.52 KB)2010-9-30 14:15圖3自定義西

11、安80坐標(biāo)系 打開map_proj.txt，可以看到新建的坐標(biāo)信息已經(jīng)自動加入。下載 (130.92 KB)2010-9-30 14:15圖4 自定義坐標(biāo)列表5 使用自定義坐標(biāo)系 下面將利用自定義坐標(biāo)系將一副北京54坐標(biāo)系轉(zhuǎn)化為西安80坐標(biāo)系。試驗的柵格數(shù)據(jù)情況為：一幅北京坐標(biāo)系的柵格數(shù)據(jù)，投影參數(shù)如下：投影類型：Transverse Mercator橢球：Krassovsky基準(zhǔn)面：Krassovsky（自定義）中央經(jīng)線：117東向偏移：500000m 由于數(shù)據(jù)的投影信息不是國際標(biāo)準(zhǔn)或者說其參數(shù)名稱不是標(biāo)準(zhǔn)的，所以在ENVI中有可能不能讀取數(shù)據(jù)的投影信息（如圖5），這個時候就需要重新設(shè)定投影

12、信息。下載 (70.47 KB)2010-9-30 14:15圖5查看坐標(biāo)信息 打開數(shù)據(jù)文件，在Available Bands List中選擇文件，點擊右鍵打開Header Info，在Edit Attributes中選擇Map Info。選擇Change Projection，將前面定義好的北京54坐標(biāo)選上。下載 (124.45 KB)2010-9-30 14:15圖6添加坐標(biāo)系 這樣在影像的map info里面就可以看到投影信息了。下載 (64.96 KB)2010-9-30 14:15圖7 更改坐標(biāo)信息后的效果 選擇Map-Convert Map Projecton,進(jìn)行西安80的坐標(biāo)轉(zhuǎn)換。下載 (141.6 KB)2010-9-30 14:15圖8 轉(zhuǎn)換參數(shù) 打開轉(zhuǎn)換后的結(jié)果，如圖8所示，可以看到Map Info中投影信息已經(jīng)更