etm是什麼數據

⑴ 空間數據

所有與地理位置有關的圖形數據都稱為空間數據。按照數據的種類，將塔里木河流域的空間數據分成遙感影像、空間基礎地理圖形、生態環境專題圖形三大類。

(一)遙感影像數據

塔里木河生態環境遙感影像數據是整個監測系統可持續運行的基礎，為形成多層次、運行化的動態監測體系，需採用多種衛星圖像數據，數據類型包括MODIS、LANDSATTM/ETM⁺、ASTER、CBERS、SPOT和QUICKBIRD等衛星圖像。

1.MODIS數據

MODIS衛星圖像數據，覆蓋范圍大，更新周期短，是宏觀動態監測的主要信息源。掃描寬度2330km，單景即可覆蓋塔里木河整個流域，時相為自2003年6月開始，數據更新周期為15d。MODIS數據共有36個波段，CH1、CH2空間解析度為250m，CH3至CH7為500m，CH8至CH36為1000m。

2.LANDSATTM/ETM⁺數據

TM圖像數據在可見光－熱紅外(0.45～12.5μm)共有7個波段，TM1～TM5、TM7地面解析度為30m，TM6為120m。ETM⁺1～ETM⁺5、ETM⁺7地面解析度為30m，ETM⁺6為60m，ETM⁺Pan(0.52～0.90μm)為15m。TM/ETM⁺數據單景幅寬185km×185km。覆蓋塔里木河全流域需要TM/ETM⁺50景，其中「四源一干」范圍29景。生態環境本底調查以2002年為准，動態監測每兩年一次。另外，為開展塔里木河幹流上游土壤鹽漬化監測，每監測周期採集2月下旬至3月初的TM/ETM⁺圖像5景(ASTER數據45景)。

3.ASTER數據

LANDSAT－7失效後選用ASTER圖像替代ETM⁺數據。其Band1～Band3相當於TM/ETM⁺2～TM/ETM⁺4，但地面解析度達15m;Band4～Band9對應TM/ETM⁺5、TM/ETM⁺7，地面解析度也是30m;Band10～Band14對應TM/ETM⁺6，但其地面解析度為90m。ASTER數據單景幅寬60km×60km。共收集2004年7～9月覆蓋塔里木河流域「四源一干」范圍ASTER數據120景。

4.CBERS數據

中巴地球資源衛星(CBERS)是我國和巴西共同發射的國土資源勘查衛星，其圖像數據可作為美國LANDSATTM/ETM⁺、ASTER圖像數據的補充。使用CBERSB1～B5數據，地面解析度19.5m，B1～B4對應於TM/ETM⁺1～TM/ETM⁺4，B5波段相當於ETM⁺Pan數據。單景幅寬113km×113km，共收集塔里木河全流域所需數據119景，時相為2002年7～9月。

5.SPOT－5數據

SPOT－5有5個較窄的、更適宜的光譜段。全色波段影像空間解析度達2.5m，多光譜波段為10m，數據單景幅寬60km×60km。覆蓋烏斯曼(中游)和喀爾達依(下游)兩個重點區PA+多光譜各一景，從2003年7～9月開始以1年為一個周期進行生態監測。

6.QUICKBIRD數據

QUICKBIRD衛星掃描寬度16.5km，星下點全色波段空間解析度0.61m，4個多光譜波段2.44m。系統在塔里木河下游末端台特瑪湖230km²范圍利用QUICKBIRD圖像監測調水方案實施後水面及植被恢復狀況，自2003年9月起，以兩年為一個周期進行監測。

(二)空間基礎地理圖形數據

本系統按照不同的管理范圍採集了不同尺度的基礎地理圖形數據，這些基礎圖形數據為分級管理提供了基礎數據資料，採集的數據包括1∶50萬、1∶10萬、1∶5萬、1∶1萬4種尺度。

1∶50萬地形圖主要用於信息系統整體圖面的顯示控制，是區內1990～1991年最新數據，覆蓋塔里木河全流域，面積約125×10⁴km²，數據成果為數字線劃圖(DLG)。1∶10萬地形圖主要用於遙感影像幾何校正和系統信息可視化的基礎地理數據，數據范圍為塔里木河流域「四源一干」區域，面積23.63×10⁴km²，最終數據成果包括數字柵格圖(DRG)、數字線劃圖(DLG)、數字高程模型(DEM)三種產品，為20世紀七八十年代測制。1∶10萬DLG數據主要作為系統查詢、分析的背景數據，為快捷方便地進行信息的可視化查詢管理提供底圖數據;DEM數據主要用於虛擬模擬和其他子系統的分析應用;DRG數據主要用於遙感影像幾何校正和快速配准。1∶5萬地形圖是對1∶10萬數據的補充，1995～1997年測制，范圍覆蓋塔河幹流區域，最終成果為DLG數據。塔里木河管理的重點是幹流河道兩岸，2002年沿幹流河道範圍利用航攝方法採集了1∶1萬基礎地形數據，數據成果包括數字正射影像圖(DOM)、數字高程模型(DEM)、數字線劃圖(DLG)3種產品。

(三)生態環境專題圖形數據

生態環境專題圖形數據包括塔里木河流域生態環境調查、監測所形成的圖件和由本系統遙感專題信息提取所產生的圖形資料。

塔里木河幹流區域是生態保護和監測的重點，對這一區域生態環境的變化分析不僅需要現實數據資料，也需要大量的歷史數據。本次系統收集了幹流地區1980年以來的八大類歷史專題圖形數據，包括草場、荒漠化、土地類型、土地資源評價、植被、土地利用、鹽漬化、地貌，數據比例尺為1∶10萬。另外，對在幹流河道1∶1萬地形測繪中所形成的1∶1萬和1∶5萬水系圖、植被分布圖和土地利用圖也已入庫。「四源一干」范圍主要收集了1990年和2000年1∶10萬土地利用專題圖，20世紀80年代和90年代1∶20萬水文地質圖。全流域范圍採集了1995年1∶50萬土壤專題圖、1995年1∶100萬草地資源專題圖、2005年塔里木河流域水文站網分布圖、重要水利工程分布圖等。

本系統運行後，還將不斷產生和生成新的專題圖數據。利用遙感技術對塔里木河流域進行植被、土地利用、沙質荒漠化、土壤鹽漬化的專題信息提取，提取的專題圖形數據是生態環境監測、評估的基礎。由遙感圖像提取的專題信息主要有:

塔里木河幹流中下游植被類型與蓋度圖(1∶10萬)，從2002年開始，每兩年一期;

塔里木河「四源一干」土地利用圖(1∶10萬)，從2002年開始，每兩年一期;

塔里木河幹流下游土地沙質荒漠化圖(1∶10萬)，從2002年開始，每兩年一期;

塔里木河幹流上游土壤鹽漬化圖(1∶10萬)，從2002年開始，每兩年一期;

塔里木河幹流烏斯滿地區植被類型與蓋度圖(1∶1萬)，從2003年開始，每年一期;

塔里木河幹流喀爾達依地區植被類型與蓋度圖(1∶1萬)，從2003年開始，每年一期;

塔里木河幹流台特瑪湖地區植被類型、蓋度與水體面積分布圖(1∶1萬)，從2003年開始，每年一期。

(四)GPS控制點數據

2004～2005年在四源流范圍共布設了404個GPS控制點，主要用於水利工程設計、施工，河道變形監測以及遙感影像圖校正。

⑵ 遙感數據類型及數據處理

開展遙感地質找礦工作的關鍵之一是精心選擇和獲取遙感數據。不同構造單元具有相異的地理地貌特徵、不同成礦條件和各自的成礦地質背景。所以選擇與研究區成礦特徵相適應的遙感數據，是遙感地質找礦取得良好效果的保證。

一、TM與ETM⁺數據特徵

陸地衛星TM和ETM+分別屬於美國陸地資源衛星的第二代和第三代感測器系統，具有一定的繼承性。兩者之間的多光譜波段數相同，都為7個波段，ETM+增加了一個地面解析度為15m的全色（PAN）波段，以及將TM6波段的空間解析度提高到60m。多光譜數據空間分辨滿足1∶20萬的制圖要求。對重點礦區（帶）進行遙感地質解譯，可以通過8波段與多光譜數據融合方法將影像空間解析度提高到15m，能夠達到1∶5萬地圖草測精度。因此該數據能夠滿足本次項目中對遙感地質信息提取的要求。

表5-1列出了TM、ETM+遙感數據的主要性能指標。可以看出，ETM（TM）1～4波段為可見光及近紅外光，含鐵礦物在此波段有明顯的光譜吸收特徵，對於探測乾旱、半乾旱環境中含褐鐵礦的岩石、土壤有較好效果；ETM7對於出露地表的粘土與碳酸鹽礦物敏感，將ETM7和ETM2～5組合分析，對於識別含鐵粘土礦物非常有效，這類礦物通常是熱液蝕變的標志；TM6可以探測地表的熱輻射強度，對於識別與金礦化密切相關的石英脈本身及含石英較多的岩石是非常有用的。總之，能從TM圖像中提取三種基本的光譜信息，即鐵染強度、岩石土壤中的羥基和碳酸鹽根礦物，是地質找礦的十分重要的資料。

表5-1 TM和ETM+的主要性能參數及光譜識別標志

二、工作使用的數據

研究區分布較廣，主要涉及研究區的圖像有3景，數據景號與採集日期見表5-2。

表5-2 研究區內遙感數據景號與採集日期

三、多波段遙感數據最佳波段組合優選

由於人們對彩色敏感程度比對全色大得多，遙感圖像應用研究中多採用RGB彩色合成圖像及彩色空間變換圖像來進行構造、地層岩性和岩體等解譯。然而ETM+圖像有7個波段，各波段物理性質不同，且各波段之間信息又有一定重疊和冗餘。一般圖像的方差越大，包含的信息越多；相關系數越小，波段間信息冗餘度越小。所以要求圖像方差要大而相關性要小這樣兩個條件。基於上述依據，目前最佳波段選擇的主要方法有：熵與聯合熵、OIF指數法、方差-協方差矩陣特徵值法等。其中OIF指數法，該方法是美國查維茨提出，理論依據是：圖像數據的標准差越大，所包含的信息量也越多，波段的相關系數越小，表明各波段的圖像數據獨立性也就越高，信息的冗餘度也就越小。由於計算方法簡單，易於操作，所以最常用。但這些波段選擇方法只考慮到兩個方面：信息量和相關性，而沒有考慮到研究對象的光譜特徵和衛星感測器的用途，即要使波段組合後研究地物的光譜特徵差異最大。綜合各方面因素，在本文中採用「最佳指數（OIF）+光譜特徵分析」法。

四、遙感圖像處理

遙感圖像處理過程中始終以區域控礦地質理論為基礎，結合工作區的自然地理地貌環境，在充分總結成礦規律的基礎上，從遙感圖像或數據中提取不同層次、不同內容的與成礦有關的控礦要素，圈定成礦遠景區。

數據處理過程包括：

1）數據預處理：包括圖像數據分析，校正，配准，子區裁剪等操作。

2）數據處理：包括圖像增強、信息提取等。主要有兩方面工作，即圖像分類、解譯和成礦信息提取。

3）生成專題圖層：研究區構造格架、影像構造單元劃分，蝕變遙感異常信息以及成礦位場等圖層，為多元信息統計分析提供數據源。

遙感圖像處理流程（圖5-1）。

五、遙感蝕變異常提取

（一）遙感技術應用於斑岩銅礦勘查的理論基礎

1.蝕變信息提取的物理學基礎

近三十年來，中外學者進行了岩石和礦物波譜特性的大量研究工作，這些研究涉及晶體場理論、礦物學、固體物理學、量子力學、遙感岩石學等眾多領域。在多年岩石和礦物波譜特徵研究的基礎上，特別是80年代發射的陸地衛星LandsatV增設了兩個短波紅外波段：TM5（1.55～1.75μm）和TM7（2.08～2.35μm），為找礦提供了可以提取具有找礦標志意義的熱液蝕變遙感信息。美國國家宇航局（NASA）1998年發射的主力遙感衛星LandsatⅦ號進一步將全色波段的空間解析度由LandsatV衛星的30m提高到15m。

圖5-1 遙感圖像處理流程

從HuntG.R（1978）和他領導的實驗室的研究成果以及閻積惠等（1995）依據礦物反射波譜特徵吸收譜帶特點的定性分類研究中可以知道：主要造岩礦物在可見光—近紅外光譜（0.35～2.5μm）並不產生具有鑒定意義的反射譜帶，其光譜特徵主要由岩石中為數不多的次要礦物決定：

一是含鐵（Fe²⁺、Fe³⁺）基團產生，含鐵礦物主要有角閃石、赤鐵礦、褐鐵礦、針鐵礦、磁鐵礦、黃鉀鐵礬等，他們在TM1—TM4波段有強的吸收帶，若岩石中含多量的Fe³⁺，而含Fe²⁺很少，這類岩石的主要吸收譜帶位於TM4和TM1波段，反射波長相當於TM3波段的電磁波。若含大量的Fe²⁺、含Fe³⁺很少，則主要吸收譜帶位於1波段，對於波長相當於TM2波段的電磁波有某種程度的反射。

二是含羥基（OH^-）、水（H₂0）或碳酸根基團產生，羥基的吸收譜帶主要有二處：2.2μm，2.3μm，由於OH^-在2.2～2.3μm附近存在強吸收谷（稱為羥基譜帶），使得TM7產生低值，TM5產生高值，含羥基礦物大多為次級蝕變礦物，如高嶺土、葉蠟石、雲母類礦物、綠泥石、綠簾石等，水在1.4μm和1.9μm處有特徵吸收帶。含碳酸根礦物主要有五個特徵吸收譜帶（1.9～2.55μm），較強的兩個在2.35μm和2.55μm波長處（稱為碳酸根譜帶），相對較弱的在1.9μm、2.0μm、2.16μm三處。常見礦物包括方解石、白雲石、石膏、菱鎂礦等。

2.遙感應用的地質學基礎

從控礦因素的分析得知，成礦規律的研究必須從分析控制和影響礦床形成的各種地質因素著手，通過對地層、岩體、斷裂、蝕變等地質因素分析，確定控礦地質因素，分析其對成礦有利的程度。

圍岩蝕變現象作為找礦標志已有數百年歷史，發現的大型金屬、非金屬礦床更是不勝枚舉：北美、俄羅斯的大部分斑岩銅礦、我國的銅官山銅礦、猶他州的大鋁礦、西澳大利亞的大型金礦、墨西哥的大鉑礦、美國許多白鎢礦、世界大多數錫礦、哈薩克的剛玉礦等。有用元素的逐步富集是形成礦床的根本，而這種成礦物質通常由成礦熱液進行遷移搬運和卸載沉澱。近礦圍岩蝕變主要是不同類型的熱液與周圍岩石相互作用的產物，是成礦物質逐步富集成礦過程中留下的印跡。熱液來源主要有：岩漿作用有關的熱液、變質作用有關的熱液以及地下水環流熱液等。熱液流體在運移過程中由於其溫度、物理化學條件以及圍岩組分差異，形成不同的蝕變礦物，按照圍岩組分可以分為三類：中酸性岩的蝕變（雲英岩化、絹雲母化、鈉長石化、鉀長石化）、基性、超基性岩的蝕變（蛇紋石化、綠泥石化、青磐岩化、碳酸鹽化）、石灰岩及其他碳酸鹽類的蝕變（矽卡岩化、硅化、重晶石化、白雲岩化）。並且不同類型的圍岩蝕變指示礦床存在的意義不同，這是由於各種蝕變和礦石同樣是成礦作用的產物，從而在時間空間上礦和有關的蝕變存在密切的聯系，並且不同類型的蝕變及其組合專屬一定的礦床（礦種或礦床類型），因而可以用蝕變來預測礦化礦產資源（Курекнн.1954）。近礦圍岩蝕變的研究和熱液蝕變岩石的發現，可以指示找礦的空間范圍，增加找到礦床的機會。而斑岩銅礦的形成過程中蝕變特徵尤其明顯。

（二）蝕變遙感異常信息提取方法

本次工作中主要採用的是主分量分析法、比值法增強蝕變遙感異常信息，通過閾值處理（決策樹技術）分級提取。工作過程為：預處理（去干擾）→信息提取→異常優化。

1.預處理

（1）高端切割去雲及鹽鹼地的影響

根據圖像采樣統計結果，雲的光譜特徵在TM（ETM）1、2、3的灰度值相對較高，TM（ETM）5、7的灰度值相對較低；鹽鹼地一般表現為在TM（ETM）3的灰度值相對較高，其次為TM（ETM）5、7、2、1。本次工作採用ETM1高端切割來去除雲的干擾，採用ETM3的高端切割來去除鹽鹼地的影響，效果較好。

（2）比值法消除植被影響

植被是蝕變異常提取過程中常見的干擾因素。工作區雖然植被覆蓋較少，但為了能盡可能完全地提取致礦異常，還是做了此項處理。根據植物的光譜曲線特徵綠色植被在TM（ETM）4（0.76～0.90μm）的反射率最高，可以認為，只有當有其他因素影響時，TM（ETM）5才可能大於TM（ETM）4，因此，可以選取ETM5/4≤1來消除植被干擾。

2.信息提取

一般常見的圖像處理方法有：主分量分析、光譜角法、比值法等。下面簡要介紹這幾種方法的基本原理：

（1）主分量分析（PCA：PrinCipal Component Analysis）又稱主成分分析，在計算機處理中稱K—L變換。TM數據在圖像處理系統中經K-L變換將TM圖像轉變為一組互不相關的表徵函數序列，目的在於壓縮TM的波譜維數、突出地物類別、提取與礦化有關的蝕變信息。K-L變換在數學含義上，它是一種基於圖像統計特徵的多維正交線性變換。經過這種變換後生成一組新的組分圖像（數目等於或小於原波段數），是輸入的若干圖像的線性組合，即：

新疆北部主要斑岩銅礦帶

其中，X是原多波段圖像的數據矩陣，矩陣元素為P個波段的像元值向量；Y是輸出的主組分矩陣，即q個組分的像元值向量，一般q≤p；T為變換核矩，通常為由變換波段之間的協方差矩陣所產生的特徵向量矩陣。

新疆北部主要斑岩銅礦帶

y₁=t₁₁x₁+t₁₂x₂+…+t_1px_p 第一組分

…

y_q=t_q1x₁+t_q2x₂+…+t_qpx_p 第q組分

y₁，y₂…y_q按協方差矩陣的特徵值大小依次排序。

K—L變換後，第一組分（y₁）取得最大信息量（可達90%左右），其餘依次減少，一、二、三組分基本是已集中了絕大部分信息，後面組分包含的信息量往往已非常小。

（2）光譜角法

把每一個多維空間點以其空間特徵向量來表徵，並以空間向量角的相似性作為判據。它是一種監督分類，要求每一類別有一個已知參考譜。此參考譜可以是地面實測入庫光譜，也可以是已知條件的圖面單元的統計入庫結果（又稱圖像采樣）。為了直觀，設三維空間點P在彩色坐標系中的特徵向量為OP，以此向量為軸作小角錐（圖5-2），凡位於此小角錐內的空間點都視為相似的。

圖5-2 光譜角法原理

根據線性代數理論，向量α、β間夾角θ為

新疆北部主要斑岩銅礦帶

式中，（d，β）為n維向量d，β的內積，|α|、|β|分別為向量d，β的長度，當存在已知礦點或礦床時，可以利用光譜角法圈定與其有相似譜特徵的成礦遠景區，以減少主分量分析所獲異常中的非礦異常；當存在兩種以上已知礦點時，可以用光譜角法對主分量分析異常進行類別區分。這兩點是光譜角法在異常信息提取中對主分量分析法可以起的輔助作用。

（3）比值法

大量的研究成果表明，蝕變礦物在不同TM波段之間存在光譜反差，其中含羥基的粘土礦物和碳酸鹽礦物，在TM7波段具有強吸收，在TM5波段為強反射，而褐鐵礦在TM3表現為高反射，在TM1、TM2和TM4則具不同程度的吸收特徵，故TM5/TM7、TM3/TM1、TM5/TM4和TM4/TM3通常可用於增強提取上述特定的熱液蝕變信息。

⑶ 遙感數據種類

在區域地質調查中，應用遙感技術經歷了從黑白航空相片目視解譯到廣泛應用多平台多感測器航空航天遙感各種信息的過程。目前隨著感測器系統的探測能力、質量、品種和解析度的大大提高，可供遙感地質應用的航空航天遙感數據越來越多，其應用領域也不斷擴大，遙感數據處理、解譯、成果和交流也正在逐漸向數字化和自動化方向發展。在中小比例尺的區域遙感地質調查中，已形成了以航天遙感數據為主，地面解析度高的航空遙感為重要補充的技術格局；同時，遙感技術與地質、物化探等地學學科，以及其他領域的某些新技術、新方法（如GIS、GPS等）緊密融合，使區域地質調查中的遙感應用成為多學科、多技術互相滲透的綜合應用體系。

目前國內使用的星載多光譜遙感系統主要為表1-1所列。

表1-1 目前常用的星載遙感技術系統