王寧，趙新生，陳川

(1.新疆大學(xué)地質(zhì)與礦業(yè)工程學(xué)院，烏魯木齊 830049；2.新疆大學(xué)新疆中亞造山帶大陸動(dòng)力學(xué)與成礦預(yù)測(cè)重點(diǎn)實(shí)驗(yàn)室，烏魯木齊 830049)

1 引 言

遙感技術(shù)因其具有區(qū)域性、經(jīng)濟(jì)性和快速性的特點(diǎn)，在地質(zhì)工作中起到越來(lái)越重要的作用，隨著地質(zhì)工作的深入發(fā)展和科學(xué)技術(shù)的不斷進(jìn)步，我國(guó)地質(zhì)填圖工作也在越來(lái)越普遍地使用遙感方法和技術(shù)。

遙感圖像具有視域廣闊、信息豐富、地貌形態(tài)逼真、地質(zhì)概括力強(qiáng)等特點(diǎn)，利用它進(jìn)行地質(zhì)解釋，能夠快速、準(zhǔn)確、全面的提供超前編圖信息。同時(shí)遙感影像的比例尺較大，因而可以反映出更豐富的地質(zhì)信息，能夠極大減少工作量，縮短工期，起到事半功倍的作用。同時(shí)，遙感技術(shù)的應(yīng)用受到客觀因素的限制，如遙感圖像受氣候及植被等自然條件影響，造成解譯準(zhǔn)確度降低[1]。

因此如何提高遙感數(shù)據(jù)的應(yīng)用效率成為關(guān)鍵問(wèn)題。由于ASTER數(shù)據(jù)的光譜和空間分辨率均具有一定優(yōu)勢(shì)，適合于中等比例尺的地質(zhì)填圖，因此本文以ASTER數(shù)據(jù)在新疆吉木乃縣托斯特鄉(xiāng)附近1∶5萬(wàn)區(qū)域地質(zhì)調(diào)查為例，綜合說(shuō)明遙感數(shù)據(jù)應(yīng)用效率在地質(zhì)填圖工作中的優(yōu)勢(shì)。

2 區(qū)域地質(zhì)背景及研究方法評(píng)述

2.1 區(qū)域地質(zhì)背景

研究區(qū)位于吉木乃縣托斯特盆地北部，烏倫古湖以西，福海盆地西部的山間凹陷，是在古生代基底上發(fā)育起來(lái)的由齋桑泊南-烏倫古湖斷裂和莫合爾臺(tái)南斷裂夾持的新生代斷陷盆地[2]。區(qū)域大地構(gòu)造位置為準(zhǔn)噶爾古板塊薩吾爾晚古生代島弧區(qū)中西段。區(qū)內(nèi)石炭系、二疊系、新近系及第四系地層均有發(fā)育，出露較好且分布廣泛；該地區(qū)巖漿活動(dòng)頻繁，侵入巖種類多，分布范圍廣，形成于后造山運(yùn)動(dòng)和活動(dòng)大陸邊緣?。粎^(qū)內(nèi)構(gòu)造較為復(fù)雜，褶皺、斷裂較發(fā)育，斷裂以北東向壓性、壓扭性為主，這些斷裂主要形成于華力西中晚期和燕山一喜馬拉雅期；巖石類型包括超基性巖、基性巖、中性巖、酸性巖和堿性巖等，如正長(zhǎng)花崗巖、花崗斑巖、閃長(zhǎng)巖、輝長(zhǎng)巖等。

2.2 ASTER數(shù)據(jù)來(lái)源

ASTER是搭載在美國(guó)NASA Terra衛(wèi)星上的一種高級(jí)光學(xué)傳感器，可以采集從可見光到熱紅外共3個(gè)通道的14個(gè)波段的遙感數(shù)據(jù)[3]，可見光-近紅外光譜域(0.52um～0.86um)內(nèi)設(shè)有3個(gè)波段(空間分辨率為為15m)，短波紅外譜域(1.60um～2.43um)內(nèi)有6個(gè)波段(空間分辨率為30m)，熱紅外譜域范圍(8.125um～11.65um)內(nèi)設(shè)有5個(gè)波段(空間分辨率為90m)，本研究中只利用了可見光-近紅外、短波紅外共9個(gè)波段(空間分辨率為15m)，沒有采用熱紅外波段[4]。

與ETM+數(shù)據(jù)相比，ASTER 數(shù)據(jù)的空間分辨率和光譜分辨率都有一定程度的提高，刈幅面積為60km×60km，適宜進(jìn)行制作1∶5萬(wàn)～1∶10萬(wàn)比例尺的遙感影像圖，這也為本文選擇ASTER數(shù)據(jù)進(jìn)行研究區(qū)的1∶5萬(wàn)遙感地質(zhì)填圖提供了一定依據(jù)。

2.3 遙感地質(zhì)填圖方法及應(yīng)用

本文所選擇的1∶5萬(wàn)遙感地質(zhì)填圖的方法包括影像單元法，影像巖石單元法和單元-剖面法等方法技術(shù)，這些方法技術(shù)貫穿于整個(gè)遙感填圖過(guò)程[5]。

(1)影像單元法

影像單元法是根據(jù)遙感圖像的信息進(jìn)行解譯，分類命名，最后形成遙感影像單元解譯圖的一種方法。影像單元可劃分為定性類影像單元，半定性類影像單元和非定性類影像單元，分類以影像單元建立劃分為基礎(chǔ)而實(shí)施，其將影像單元特征的表現(xiàn)形式和地質(zhì)屬性分析判斷的粗略作為依據(jù)。

定性類影像單元是指通過(guò)已建立的影像單元直接實(shí)現(xiàn)地質(zhì)體屬性的定性判斷。在本文中主要為斷裂構(gòu)造。

半定性類影像單元是指通過(guò)影像特征分析、對(duì)比，可基本判斷單元地質(zhì)屬性的一類影像單元?？蓜澐譃榄h(huán)形影像單元，層型影像單元和非層型影像單元。本研究為了增強(qiáng)解譯效果，對(duì)圖像進(jìn)行了增強(qiáng)和濾波處理，更清晰地表現(xiàn)影像特征，達(dá)到定性的目的。

非定性類影像單元是指根據(jù)影像單元特征難以實(shí)現(xiàn)地質(zhì)屬性歸類的影像單元。在本研究中，經(jīng)過(guò)影像處理之后，基本達(dá)到半定性目的，提高解譯精度。

(2)影像巖石單元法

影像巖石單元法是以影像巖石單元為基本填圖單位，進(jìn)行地質(zhì)圖初步解譯，遙感地質(zhì)填圖制圖的方法。為了更好的解譯巖性，本研究在一些典型地區(qū)進(jìn)行野外采樣，增強(qiáng)可靠度。

(3)單元-剖面法

單元-剖面法是以影像單元為調(diào)查單位，實(shí)施野外踏勘剖面選擇、測(cè)制和踏勘路線布置和調(diào)查的一種方法。在野外地質(zhì)調(diào)查過(guò)程中，實(shí)測(cè)主干剖面的布置以影像單元齊全，構(gòu)造簡(jiǎn)單，地層發(fā)育齊全地區(qū)為原則，通過(guò)全面分析研究區(qū)的地質(zhì)特征，在研究區(qū)布置了2條剖面(Ⅰ剖面和Ⅱ剖面)，基本上控制研究區(qū)內(nèi)所有的巖性類型。重點(diǎn)調(diào)查巖石類型、結(jié)構(gòu)構(gòu)造、巖石類型組合特征和影像單元界線性質(zhì)等特征，從而研究影像單元與火山巖，沉積巖體之間的對(duì)應(yīng)關(guān)系。

3 巖性(地層)遙感解譯

3.1 圖像預(yù)處理

首先對(duì)研究區(qū)遙感影像圖進(jìn)行大氣校正、幾何精校正、圖像裁剪和波段融合等圖像預(yù)處理，預(yù)處理之后則要進(jìn)行波段組合選擇，波段組合的選擇是多波段遙感影像合成中十分重要的過(guò)程。本文所選擇的ASTER遙感數(shù)據(jù)共14個(gè)波段，除了熱紅外波段不參與波段組合外，從其余9個(gè)波段中選取3個(gè)波段進(jìn)行組合，來(lái)更好的表達(dá)圖像信息，本文根據(jù)協(xié)方差矩陣，對(duì)每個(gè)波段進(jìn)行分析，選擇最佳波段組合。

因?yàn)閰f(xié)方差既反映了不同變量之間的相關(guān)程度，又反映了各個(gè)分量各自取值的離散程度，因此選擇波段的主要要求是該波段的協(xié)方差應(yīng)盡可能大。為了提高彩色合成圖像的視覺效果以及地物之間的可辨程度，現(xiàn)從研究區(qū)ASTER遙感數(shù)據(jù)9個(gè)波段之間的協(xié)方差(表 1)的8波段與VNIR(可見光-熱紅外)的4波段和1波段進(jìn)行假彩色合成(圖1)，相比其他波段合成，其圖像在空間分辨率及色彩豐富程度方面優(yōu)勢(shì)顯著。

表1 協(xié)方差矩陣

圖1 研究區(qū)遙感影像圖(ASTER數(shù)據(jù)8、4、1波段假彩色合成)

3.2 圖像增強(qiáng)處理和圖像分類

3.2.1 圖像增強(qiáng)處理

研究區(qū)遙感影像主要運(yùn)用主成分分析(PCA)去除波段之間多余信息，將多波段的圖像信息壓縮到比原波段更有效的少數(shù)幾個(gè)轉(zhuǎn)換波段，去相關(guān)拉伸消除ASTER數(shù)據(jù)中波段間的高度相關(guān)性，以期達(dá)到圖像增強(qiáng)的目的[6-7]。

3.2.2 圖像分類

圖像分類的目的是將圖像中每個(gè)像元根據(jù)其在不同波段的光譜亮度、空間結(jié)構(gòu)特征或者其他信息，按照某種規(guī)則或算法劃分為不同的類別[8]，是遙感圖像處理的重要部分。遙感圖像分類一般方法包括監(jiān)督分類和非監(jiān)督分類。監(jiān)督分類是通過(guò)選擇訓(xùn)練區(qū)樣本對(duì)整個(gè)圖像所有像元進(jìn)行分類的方法；非監(jiān)督分類是在沒有先驗(yàn)類別知識(shí)的情況下，根據(jù)圖像本身的統(tǒng)計(jì)特征及自然點(diǎn)群的分布情況來(lái)劃分地物類別的分類處理[9]。鑒于二者的區(qū)別，結(jié)合研究區(qū)的實(shí)際情況，本文主要運(yùn)用ENVI軟件中的監(jiān)督分類方法對(duì)圖像進(jìn)行處理。

(1)定義訓(xùn)練樣本

定義訓(xùn)練樣本一般用目視解譯，洪積層在遙感圖像RGB(8，4，1)上顯示淡紫色色調(diào)。沖積層分布于山麓帶，河流出口處，呈灰色色調(diào)，扇形形態(tài)，灰色均勻影紋圖案，部分有植被覆蓋，為綠色色調(diào)，且本研究區(qū)有呈辮狀的河流穿過(guò)沖積層；火山碎屑巖類則為略顯凹凸的紫灰色色調(diào)，其上有輕微植被覆蓋；凝灰礫巖與火山碎屑巖色調(diào)較像，只是顏色上略顯深色；石英砂巖在影像上顯示為亮白色色調(diào)，易于區(qū)分；花崗巖邊界明顯，在影像上部分顯團(tuán)塊狀，有少量植被覆蓋；粗砂巖在影像上凹凸明顯，深暗灰色；泥質(zhì)巖附近有斷裂發(fā)育，呈條帶狀分布；湖泊在影像上最好區(qū)分，且邊界明顯。

為了保證定義訓(xùn)練樣本的真實(shí)性，本文以遙感影像單元為調(diào)查單位，以研究現(xiàn)有資料為參考，實(shí)施野外踏勘布置剖面(圖2)，剖面Ⅰ位于研究區(qū)東部，剖面方向?yàn)榻媳毕?，全長(zhǎng)9.3km，剖面Ⅱ位于研究區(qū)西部，剖面方向?yàn)楸睎|南西方向，全長(zhǎng)5km；主要巖性單元特征總結(jié)見表2。

圖2 Ⅰ、Ⅱ?yàn)樵谘芯繀^(qū)的兩條實(shí)測(cè)剖面圖

表2 實(shí)測(cè)剖面所涉及的地層巖性單元及其特征

(2)執(zhí)行監(jiān)督分類

以校園建筑設(shè)計(jì)風(fēng)格為依據(jù)，校園文化景觀建設(shè)風(fēng)格定位為現(xiàn)代簡(jiǎn)約。以基地的地址特征、城市文脈、建筑肌理為基本依據(jù)，盡量采用原生材料及植物，硬化廣場(chǎng)景觀與自然生態(tài)景觀相結(jié)合，微地形成的多層次景觀與相應(yīng)軟景塑造相結(jié)合。通過(guò)現(xiàn)代簡(jiǎn)約景觀設(shè)計(jì)手法組織相應(yīng)的景觀設(shè)計(jì)元素，賦予整體景觀設(shè)計(jì)故事性和藝術(shù)性，運(yùn)用對(duì)構(gòu)、障景等手法，營(yíng)造景觀焦點(diǎn)，結(jié)合自然的草坡、溝渠小路，原生水體等豐富景觀層次，凸顯現(xiàn)代校園與自然生態(tài)完美結(jié)合的綠色生態(tài)校園景觀。

執(zhí)行監(jiān)督分類，需要進(jìn)行一定的計(jì)算。監(jiān)督分類的算法一般分為平行六面體分類法、馬氏距離分類法、最小距離分類法、最大似然法、神經(jīng)網(wǎng)絡(luò)法、Parallelpipe等，經(jīng)對(duì)效果圖的對(duì)比分析，本文采取最大似然法。

最大似然法分類也稱為貝葉斯(Bayes)分類，是基于圖像統(tǒng)計(jì)的監(jiān)督分類法，是典型的和應(yīng)用最廣的監(jiān)督分類方法。它建立在Bayes準(zhǔn)則的基礎(chǔ)上，偏重于集群分布的統(tǒng)計(jì)特性，并假定訓(xùn)練樣本數(shù)據(jù)在光譜空間服從高斯正態(tài)分布，利用遙感數(shù)據(jù)的統(tǒng)計(jì)特征，計(jì)算給定像元屬于某一訓(xùn)練樣本的似然度，像元最終被歸并到似然度最大的一類當(dāng)中，最終實(shí)現(xiàn)像元的分類，由于其分類精度較高，受到廣泛關(guān)注。

由于遙感圖像的分類工作是按照?qǐng)D像光譜特征，執(zhí)行監(jiān)督分類之后，根據(jù)目視解譯，去除干擾樣本，得到處理后的訓(xùn)練區(qū)(圖3)。

(3)檢驗(yàn)

當(dāng)選取訓(xùn)練樣本時(shí)，一般要求訓(xùn)練區(qū)域具有典型性和代表性，這就要對(duì)所選樣本進(jìn)行檢驗(yàn)。因此，使用ROI可分離性工具來(lái)計(jì)算任意所選訓(xùn)練樣本間的統(tǒng)計(jì)距離，對(duì)訓(xùn)練樣本進(jìn)行評(píng)價(jià)[10]。根據(jù)可分離性值的大小，從小到大列出感興趣區(qū)組合，如果兩個(gè)參數(shù)的可分離性小于1，考慮將兩類樣本合成一類樣本；小于1.8，需要重新選擇樣本；大于1.9說(shuō)明樣本之間的可分離性好，屬于合格樣本[7]。本文研究中結(jié)合1∶20萬(wàn)區(qū)域地質(zhì)礦產(chǎn)圖及其本區(qū)遙感影像圖進(jìn)行野外取樣，野外取得的樣本作為進(jìn)行監(jiān)督分類的訓(xùn)練樣本依據(jù)，經(jīng)分類處理，符合要求，樣本合格(表3)。

圖3 監(jiān)督分類(最大似然法)處理后的圖像

(4)分類后處理

監(jiān)督分類可以為遙感解譯提供一個(gè)可靠的依據(jù)，同時(shí)結(jié)合目視解譯，可更精確地解譯地層和巖性。

表3 本研究區(qū)訓(xùn)練樣本可分離性計(jì)算報(bào)表

續(xù)表1

4 構(gòu)造特征遙感解譯

區(qū)內(nèi)地質(zhì)構(gòu)造比較發(fā)育，在遙感影像上比較明顯，易于解譯。其中，褶皺呈淺灰色環(huán)狀色調(diào)，橢圓形影像形態(tài)，橢圓環(huán)帶狀影紋圖案；斷層構(gòu)造地貌呈線狀淺灰色色調(diào)，沿?cái)嗔褞Р糠职l(fā)育斷層三角面。在ASTER數(shù)據(jù)的RGB(8，4，1)假彩色合成影像上以線性，弧形，直線狀溝谷等標(biāo)志顯示，部分構(gòu)造兩側(cè)色調(diào)差異明顯。為了增強(qiáng)解譯效果，本文對(duì)圖像進(jìn)行濾波處理，增強(qiáng)圖像構(gòu)造信息表達(dá)，在遙感圖上更容易解譯構(gòu)造。

遙感影像經(jīng)過(guò)濾波處理后，使所獲得的線性體能最佳地反映構(gòu)造形跡，然后通過(guò)目視解譯進(jìn)行歸納整理，將遙感影像與1∶20萬(wàn)區(qū)域地質(zhì)礦產(chǎn)圖套合，結(jié)合野外調(diào)查資料，運(yùn)用ArcGIS軟件進(jìn)行解譯修正[12]，最終對(duì)主要線性體賦予地質(zhì)屬性，達(dá)到構(gòu)造解譯的目的(圖4)。

圖4 研究區(qū)構(gòu)造解譯圖

從研究區(qū)的解譯圖可以看出研究區(qū)的斷層以北東向、近東西向和北西向?yàn)橹鳎c已有的地質(zhì)資料記錄基本相符。結(jié)合研究區(qū)的巖性和構(gòu)造解譯結(jié)果，以ArcGIS為平臺(tái)，經(jīng)過(guò)數(shù)據(jù)圖層的疊加處理，得到本區(qū)的遙感解譯圖(圖5)。通過(guò)對(duì)具有代表性的區(qū)域進(jìn)行野外實(shí)地檢驗(yàn)，基本滿足了1∶5萬(wàn)地質(zhì)填圖的要求，確定了解譯精度的可靠性。

圖5 研究區(qū)遙感解譯地質(zhì)圖(1∶5萬(wàn))

5 結(jié)束語(yǔ)

以ASTER遙感圖像為數(shù)據(jù)源，ArcGIS和ENVI軟件為平臺(tái)，對(duì)吉木乃縣托斯特地區(qū)進(jìn)行了1∶5萬(wàn)的遙感區(qū)域地質(zhì)填圖，針對(duì)填圖結(jié)果，得出以下幾點(diǎn)結(jié)論：

①利用ASTER多光譜數(shù)據(jù)，進(jìn)行假彩色合成變換，可增強(qiáng)目視解譯效果，在1∶5萬(wàn)地質(zhì)填圖、室內(nèi)綜合研究和制圖工作中能收到事半功倍的效果。

②人工目視解譯和計(jì)算機(jī)監(jiān)督分類兩者結(jié)合，是一種提高分類精度的有效方法。充分利用現(xiàn)有遙感數(shù)據(jù)，在GIS技術(shù)支持下，提高遙感數(shù)據(jù)分類精度，保證遙感解譯圖的質(zhì)量。

③利用遙感地質(zhì)解譯圖，可以對(duì)地質(zhì)填圖路線之間巖石地層單元和地質(zhì)構(gòu)造界線進(jìn)行準(zhǔn)確連接，可以對(duì)地質(zhì)路線穿越到的和沒有穿越到的地質(zhì)體進(jìn)行圈定，避免工作中對(duì)地質(zhì)體的遺漏，特別是基巖出露地區(qū)以及第四紀(jì)沉積地層的邊界可以進(jìn)行準(zhǔn)確圈定，因此既節(jié)約填圖時(shí)間，又提高區(qū)域地質(zhì)填圖的精度和質(zhì)量，同時(shí)與常規(guī)填圖相比，利用遙感圖像連圖能減少人為因素帶來(lái)的誤差。

參考文獻(xiàn)：

[1] 楊云龍，趙培榮.遙感技術(shù)在地質(zhì)編圖中的應(yīng)用——以庫(kù)車盆地1∶10萬(wàn)遙感地質(zhì)填圖為例[J].云南地質(zhì)，1999(3)：283-290.

[2] 王廣瑞.中國(guó)新疆北部及鄰區(qū)構(gòu)造-建造圖及說(shuō)明書[M].北京：中國(guó)地質(zhì)大學(xué)出版社，1996.

[3] 李海濤，田慶久.ASTER數(shù)據(jù)產(chǎn)品的特性及其計(jì)劃介紹[J].遙感信息，2004(3)：53-55.

[4] FUJISADA H.Design and performance of ASTER instrument[C].Proc SPIE Int Soc Opt Eng，1995 (2583)：16-25.

[5] 方洪賓，趙福岳.1∶25萬(wàn)遙感地質(zhì)填圖方法和技術(shù)[M].北京：地質(zhì)出版社，2002：32-75.

[6] 李聰聰，孫順新，王新民，等.遙感技術(shù)在云南昭通地區(qū)煤炭資源調(diào)查評(píng)價(jià)中的應(yīng)用[J].中國(guó)煤炭地質(zhì)，2010，(10)：17-21.

[7] 呂錄仕，孫順新，馮富成，等.西部煤炭資源調(diào)查評(píng)價(jià)中遙感技術(shù)的應(yīng)用[J].中國(guó)煤炭地質(zhì)，2005，17(5)：32-33，43.

[8] 張永生，鞏丹超，等.高分辨率遙感衛(wèi)星應(yīng)用：成礦模型處理算法及應(yīng)用技術(shù)[M].北京：科學(xué)出版社.

[9] 鄧書斌.ENVI遙感圖像處理方法[M].北京：科學(xué)出版社，2010：128-129.

[10] 阿衣古力·喀德爾.ASTER數(shù)據(jù)地質(zhì)填圖關(guān)鍵技術(shù)研究及應(yīng)用—以彩溝地區(qū)為例[D].新疆大學(xué)，2012.

[11] 郭文娟，張佳華.利用ASTER遙感資料提取南京城郊土地利用信息的研究[J].農(nóng)業(yè)工程學(xué)報(bào)，2005，21(9)：62-66.

[12] 田勤虎.新疆哈拉喬拉地區(qū)1：5萬(wàn)遙感地質(zhì)填圖[D].長(zhǎng)安大學(xué)，2007.