王平平， 田淑芳

(中國(guó)地質(zhì)大學(xué)(北京)地球科學(xué)與資源學(xué)院,北京 100083)

?

基于WorldView2數(shù)據(jù)的巖性信息增強(qiáng)方法
——以?xún)?nèi)蒙古扎嘎烏蘇地區(qū)為例

王平平， 田淑芳

(中國(guó)地質(zhì)大學(xué)(北京)地球科學(xué)與資源學(xué)院,北京 100083)

為了進(jìn)一步探討巖性信息增強(qiáng)方法在巖性遙感解譯中的應(yīng)用效果，選擇內(nèi)蒙古索倫山扎嘎烏蘇地區(qū)進(jìn)行巖性遙感信息增強(qiáng)方法研究。針對(duì)因異物同譜、同物異譜及表層淺覆蓋等因素導(dǎo)致的巖石影像特征復(fù)雜、色彩反差較弱及細(xì)節(jié)信息不夠豐富等一系列問(wèn)題，對(duì)研究區(qū)巖石波譜曲線(xiàn)特征、WorldView2數(shù)據(jù)光譜特征及空間特征進(jìn)行了綜合分析。在對(duì)WorldView2數(shù)據(jù)進(jìn)行幾何糾正和圖像融合等預(yù)處理的基礎(chǔ)上，應(yīng)用基于光譜特征和空間特征的一系列有效的巖性信息增強(qiáng)方法，對(duì)研究區(qū)進(jìn)行巖性解譯。將解譯結(jié)果與研究區(qū)1： 50 000比例尺地質(zhì)圖對(duì)比，劃分出了更多的巖性種類(lèi)，提高了巖性解譯的精度，為更加精確的巖性識(shí)別提供了依據(jù)。

WorldView2數(shù)據(jù)； 巖石波譜特征； 遙感巖性信息增強(qiáng)； 光譜特征； 空間特征

0 引言

WorldView2衛(wèi)星提供的高分辨率遙感數(shù)據(jù)包括0.5 m空間分辨率的全色圖像和2 m空間分辨率的8波段多光譜圖像[1-2]。近年來(lái)，遙感地質(zhì)調(diào)查技術(shù)在礦產(chǎn)資源勘查及區(qū)域地質(zhì)調(diào)查中發(fā)揮著越來(lái)越重要的作用，遙感地質(zhì)解譯是其重要內(nèi)容。以往的遙感地質(zhì)解譯工作大多以中、低分辨率遙感數(shù)據(jù)(如TM，ETM及ASTER等)為主，這類(lèi)遙感數(shù)據(jù)因受光譜分辨率和空間分辨率的限制，在遙感巖性解譯中常常只能區(qū)分巖性大類(lèi)，不能對(duì)巖性進(jìn)行詳細(xì)劃分； 而WorldView2數(shù)據(jù)在設(shè)置有較多多光譜波段的同時(shí)，大幅度提高了空間分辨率，可以很好地彌補(bǔ)這一缺陷。遙感巖性信息增強(qiáng)對(duì)于遙感巖性解譯十分重要，它可以解決因異物同譜、同物異譜、表層淺覆蓋等因素引起的巖石影像特征復(fù)雜、色彩反差較弱及細(xì)節(jié)信息不夠豐富等一些列問(wèn)題，從而提高巖性解譯的精度與準(zhǔn)確度。

本文以“內(nèi)蒙古自治區(qū)烏拉特中旗索倫山地區(qū)1∶5萬(wàn)比例尺區(qū)域地質(zhì)礦產(chǎn)調(diào)查”子項(xiàng)目為依托，以扎嘎烏蘇地區(qū)為研究區(qū)，進(jìn)行了基于WorldView2數(shù)據(jù)的巖性信息增強(qiáng)方法研究，為開(kāi)展進(jìn)一步的巖性解譯工作奠定基礎(chǔ)。

1 研究區(qū)概況及數(shù)據(jù)處理

1.1 研究區(qū)概況

扎嘎烏蘇地區(qū)位于內(nèi)蒙古自治區(qū)烏拉特中旗索倫山地區(qū)西部，北部與蒙古國(guó)接壤，在E108°00′～108°15′，N42°20′～42°30′之間。區(qū)內(nèi)有EW向邊防公路橫貫全區(qū)，交通便利(圖1)。區(qū)內(nèi)海拔高度在1 100～1 400 m，屬低山丘陵地貌。水系不發(fā)育，無(wú)永久性和季節(jié)性河流，稀少的溝谷均為干溝，零星分布的少量季節(jié)性湖泊也因連年干旱少雨而干涸。夏季干燥炎熱，春季和冬季風(fēng)大寒冷，屬典型的大陸性氣候。植被以草本植物為主，零星分布有小叢灌木，植被覆蓋度較低。巖石風(fēng)化和破碎嚴(yán)重[4]。

圖1 索倫山扎嘎烏蘇地區(qū)交通位置圖[3]

研究區(qū)處于華北板塊北緣與西伯利亞板塊南緣的結(jié)合地帶，一級(jí)構(gòu)造單元為天山―興蒙造山帶，二級(jí)構(gòu)造單元為索倫山―西拉木倫結(jié)合帶和包爾汗圖―溫都爾廟弧盆系[5]。區(qū)內(nèi)出露地層主要包括古元古界寶音圖群(Pt1B)、中元古界桑達(dá)來(lái)呼都格組(Pt2s)、奧陶系烏賓敖包組(O1-2w)、二疊系包特格組(P2b)和白堊系二連組(K2e)。主要巖類(lèi)為沉積巖和變質(zhì)巖。其中沉積巖主要包括長(zhǎng)石石英砂巖、泥巖、硅質(zhì)巖、灰?guī)r等； 變質(zhì)巖主要包括板巖、紫紅色石英巖、大理巖和片巖。區(qū)內(nèi)出露巖漿巖較少，主要為二長(zhǎng)花崗巖和閃長(zhǎng)巖。

1.2 數(shù)據(jù)源及預(yù)處理

本文使用的WorldView2多光譜圖像數(shù)據(jù)，在4個(gè)傳統(tǒng)波段(藍(lán)、綠、紅和近紅外1波段)的基礎(chǔ)上增加了4個(gè)新的波段(海岸、黃、紅邊和近紅外2波段)(表1)[6-7]。

表1 WorldView2多光譜圖像波段設(shè)置

在進(jìn)行巖性信息增強(qiáng)處理前，對(duì)WorldView2數(shù)據(jù)進(jìn)行了必要的預(yù)處理。本文利用地面控制點(diǎn)和DEM對(duì)WorldView2數(shù)據(jù)進(jìn)行了幾何糾正[8]; 采用基于正交化算法的光譜銳化高保真影像融合方法(Gram-Schmidt,GS)對(duì)研究區(qū)WorldView2的多光譜數(shù)據(jù)(圖2(a))與其全色波段數(shù)據(jù)進(jìn)行融合處理[9]。融合后的多光譜圖像(圖2(b))，影像分辨率得到提升，可解譯能力有了提高，紅圈區(qū)域內(nèi)的細(xì)節(jié)信息得到明顯增強(qiáng)，可為后續(xù)的巖性信息增強(qiáng)及巖性解譯打下了基礎(chǔ)。

(a) WorldView2原始圖像 (b) GS融合效果

圖2 GS融合假彩色合成圖像

Fig.2 False color composition image with GS fusion

2 巖石波譜特征分析

地物波譜特征的差異是用遙感方法區(qū)分地物的主要依據(jù)[10]。由于巖石的成分和結(jié)構(gòu)較為復(fù)雜，光譜特征變化復(fù)雜多樣，因此巖石本身的光譜特征不像單一礦物那樣具有特殊的可鑒定的光譜特征[11]，無(wú)法根據(jù)其光譜特征準(zhǔn)確地確定具體巖性； 但巖石的光譜特征可以反映巖石基本組成的物質(zhì)成分和結(jié)構(gòu)特點(diǎn)，可用于區(qū)分巖性的大類(lèi)。為了更好地了解研究區(qū)內(nèi)的巖石波譜特征，選取研究區(qū)內(nèi)的巖石波譜曲線(xiàn)進(jìn)行對(duì)比分析，研究不同巖石的波譜特性及其規(guī)律性，為后續(xù)的巖性信息增強(qiáng)提供依據(jù)。觀察研究區(qū)內(nèi)的泥巖、長(zhǎng)石石英砂巖、板巖、紫紅色石英巖及白色大理巖的波譜曲線(xiàn)特征(圖3)。

分析圖3，可以得到如下規(guī)律：

1)對(duì)于研究區(qū)內(nèi)的變質(zhì)巖(即板巖、紫紅色石英巖及白色大理巖)，在0.4～1.04 μm波長(zhǎng)范圍內(nèi)，白色大理巖的反射率明顯高于板巖、紫紅色石英巖及沉積巖的反射率，這與大理巖的礦物成分有很大關(guān)系——大理巖主要由淺色礦物組成，淺色礦物的反射率較高，因此大理巖的反射率最大[12]，這一特性可以作為區(qū)分大理巖與其他巖石的重要標(biāo)志。在0.74 μm附近，紫紅色石英巖出現(xiàn)反射峰，在0.875 μm附近存在反射率低谷，之后隨波長(zhǎng)的增加反射率逐漸增大。在0.4～1.04 μm波長(zhǎng)范圍內(nèi)，板巖的反射率一直處于上升趨勢(shì)，并不斷趨近于長(zhǎng)石石英砂巖。

圖3 研究區(qū)巖石波譜曲線(xiàn)與WorldView2圖像波段設(shè)置關(guān)系

2)對(duì)于研究區(qū)內(nèi)的沉積巖(如泥巖和長(zhǎng)石石英砂巖)，反射率總體較低。在0.45 μm附近，泥巖和長(zhǎng)石石英砂巖出現(xiàn)交叉現(xiàn)象，交叉點(diǎn)前后反射率的大小出現(xiàn)明顯變化； 并且隨著波長(zhǎng)的增加，二者反射率差值不斷增大。

以上規(guī)律可為后續(xù)的巖性信息增強(qiáng)提供波譜特征方面的依據(jù)，特別是在B6和B7處存在的特征反射峰和吸收谷對(duì)巖性信息增強(qiáng)具有重要的指導(dǎo)意義。為了更好地解決因異物同譜、同物異譜、表層淺覆蓋等因素引起的巖石影像特征復(fù)雜、色彩反差較弱、細(xì)節(jié)信息不夠豐富等一些列問(wèn)題，本文通過(guò)巖性信息增強(qiáng)處理，充分利用遙感影像的光譜特征和空間特征，增強(qiáng)不同巖性之間的差異，提取巖性界線(xiàn)，并結(jié)合巖石光譜特征及地質(zhì)資料進(jìn)行巖性解譯。

3 巖性信息增強(qiáng)

3.1 最佳波段組合法

波段組合的優(yōu)化是一種增強(qiáng)不同地物類(lèi)別之間差異、增強(qiáng)目視解譯效果的方法。選取最佳的波段組合，可以更加充分地體現(xiàn)地物之間的差異，使不同地物之間的差異最大化，從而達(dá)到區(qū)分不同地物的目的[13]。選取最佳波段組合時(shí)，通常綜合考慮3個(gè)方面： ①波段組合內(nèi)的波段信息量較大； ②波段組合內(nèi)波段間的相關(guān)性較??； ③該波段組合的彩色圖像中，研究區(qū)內(nèi)所關(guān)注地物類(lèi)型之間的光譜差異要大、可分性要好[14]。

3.1.1 研究區(qū)WorldView2數(shù)據(jù)統(tǒng)計(jì)特征

數(shù)據(jù)統(tǒng)計(jì)特征可以反映數(shù)據(jù)的信息量、離散程度和相關(guān)性等。波段信息量的大小可通過(guò)標(biāo)準(zhǔn)差反映，波段的標(biāo)準(zhǔn)差越大，說(shuō)明該波段像元灰度值之間的離散程度越大，像元灰度值分布的動(dòng)態(tài)范圍越大，圖像反差越明顯，不同地物之間差異越大，所包含的信息量越大； 而波段的標(biāo)準(zhǔn)差越小，圖像反差就越低，不同地物之間的差異也隨之減弱，不利于不同地物類(lèi)別之間的區(qū)分，所包含的信息量也越少。波段間的相關(guān)系數(shù)反映了波段間的相關(guān)性，相關(guān)系數(shù)越大，相關(guān)性越大，波段間信息冗余度越大； 相關(guān)系數(shù)越小，相關(guān)性越低，波段間信息冗余度越小[15-16]。

本文對(duì)內(nèi)蒙古索倫山地區(qū)扎嘎烏蘇幅WorldView2數(shù)據(jù)進(jìn)行統(tǒng)計(jì)特征分析的結(jié)果表明： ①B7的標(biāo)準(zhǔn)差最大，所含信息量最大(表2)； ②8個(gè)波段間的相關(guān)性總體上都較大； ③相鄰波段間相關(guān)性較大，特別是B1,B2,B3之間的相關(guān)性最大，但B1與其他波段的相關(guān)性低于與B2,B3的相關(guān)性(表3)。因此，在進(jìn)行波段組合時(shí)，應(yīng)選取B1； 而B(niǎo)4,B5,B6之間相關(guān)性較大，可選取其中的1個(gè)波段； B7,B8之間相關(guān)性較大，選取其中的1個(gè)波段。

表2 研究區(qū)光譜信息統(tǒng)計(jì)

表3 相關(guān)系數(shù)統(tǒng)計(jì)

3.1.2 最佳波段因子法選取波段組合

目前，較為常用的選取最佳波段組合的方法為最佳波段因子法(optimum index facter,OIF),這是一種選擇最優(yōu)RGB假彩色合成方案的方法，它基于波段組合內(nèi)波段的標(biāo)準(zhǔn)差及波段間的相關(guān)系數(shù)選取最佳波段組合[10]，即

(1)

式中：Si為第i波段的標(biāo)準(zhǔn)差；Rij為3個(gè)波段中任意2個(gè)波段間的相關(guān)性[17]。OIF越大，波段組合內(nèi)波段的標(biāo)準(zhǔn)差越大，波段間的相關(guān)性越小，這樣的波段組合形成的圖像具有較大的信息量和較小的相關(guān)性，目視效果較好。

對(duì)于不同的地區(qū)，應(yīng)該從其實(shí)際情況出發(fā)，綜合考慮該區(qū)圖像各波段的標(biāo)準(zhǔn)差、波段之間的相關(guān)系數(shù)以及巖石波譜特征分析結(jié)果來(lái)確定最佳波段組合。根據(jù)本文研究區(qū)WorldView2數(shù)據(jù)統(tǒng)計(jì)特征結(jié)果，將8個(gè)波段分為3組： ①B1； ②B4,B5,B6； ③B7,B8。分別計(jì)算其OIF，并由大到小排列，結(jié)果如表4所示。

表4 波段組合及OIF

由研究區(qū)巖石波譜特征分析結(jié)果(表2)及WorldView2數(shù)據(jù)統(tǒng)計(jì)特征分析結(jié)果(表3)可知，B6和B7對(duì)區(qū)分巖性具有重要作用； 再結(jié)合各種波段組合圖像的目視效果對(duì)比分析，發(fā)現(xiàn)B6,B7,B1(圖4(a))和B4,B7,B1(圖4(b))組合圖像包含的信息量最大，巖性差異最明顯，可作為最佳波段組合圖像用于巖性解譯。但總體來(lái)看巖性間的差異還比較弱，且由于受到異物同譜、同物異譜、表層淺覆蓋等因素影響，需要進(jìn)一步加以區(qū)分。

(a) B6(R)B7(G)B1(B)組合 (b) B4(R)B7(G)B1(B)組合

圖4 最佳波段組合假彩色合成圖像

Fig.4 False color images composed with optimum band combinations

3.2 比值法

在巖性信息增強(qiáng)中，通過(guò)選擇2個(gè)亮度差異較大的波段進(jìn)行比值運(yùn)算，可以突出地物并增大地物之間的差異。由圖3可以看出，紫紅色石英巖在波長(zhǎng)0.74 μm附近出現(xiàn)反射峰，在波長(zhǎng)0.875 μm附近存在反射率低谷，即反射峰出現(xiàn)在B6，吸收谷出現(xiàn)在B7。因此本文利用這2個(gè)波段的反射率差異進(jìn)行比值運(yùn)算(B6/B7)，并將比值結(jié)果作為1個(gè)波段與其他波段進(jìn)行假彩色合成，取得了較好的巖性信息增強(qiáng)效果(圖5)。

(a) B6(R)B7(G)B1(B)組合 (b) B4(R)B7(G)B6/B7(B)組合

圖5 比值法增強(qiáng)假彩色合成圖像

Fig.5 False color composition images after ratio method enhancement

圖5(b)中，經(jīng)過(guò)比值增強(qiáng)后，紫紅石英巖(紅色調(diào)影像)較為突出，與其他巖性之間的差異更加明顯(黃藍(lán)色調(diào)的影像為灰白色大理巖)。

3.3 HSV變換

本文對(duì)B4(R)B7(G)B1(B)假彩色合成圖像(圖6(a))進(jìn)行色度-飽和度-亮度(hue saturation value,HSV)變換； 在HSV彩色空間中對(duì)S進(jìn)行反差擴(kuò)展； 然后再反變換到RGB空間，得到增強(qiáng)后的圖像[18-20](圖6(b))。

(a) B4(R)B7(G)B1(B)組合(b) HSV變換+S反差擴(kuò)展

圖6 HSV變換增強(qiáng)假彩色合成圖像

Fig.6 False color composition images after HSV transformation enhancement

從圖6可以看出，經(jīng)HSV變換處理后，不同巖性之間的色調(diào)差異得到顯著增強(qiáng)，巖性界線(xiàn)更加清晰明顯，可用于準(zhǔn)確地圈定巖性界線(xiàn)，提高了巖性的目視解譯效果。

3.4 主成分分析

根據(jù)表3中的統(tǒng)計(jì)結(jié)果可知，各波段間的相關(guān)性較大，冗余度較高，對(duì)巖性識(shí)別造成了一定的影響。因此，可通過(guò)主成分分析法[21-24]對(duì)WorldView2圖像的8個(gè)波段數(shù)據(jù)進(jìn)行處理，將信息壓縮到少數(shù)幾個(gè)不相關(guān)的主成分中，減弱波段間的相關(guān)性。對(duì)研究區(qū)的WorldView2數(shù)據(jù)進(jìn)行主成分分析的結(jié)果表明，信息主要集中在前4個(gè)主成分中，其他成分主要包含大量噪聲。PC1主要包含大量的地貌、構(gòu)造信息[12]，巖性信息主要集中在PC3，最終選擇了B6(R)PC3(G)B8(B)假彩色合成方案(圖7(b))。

(a) B6(R)B7(G)B1(B)組合 (b) B6(R)PC3(G)B8(B)組合

圖7 主成分分析增強(qiáng)假彩色合成圖像

Fig.7 False color composition images after principal components analysis enhancement

由圖7(b)可以看出，主成分分析處理后，巖性之間的色調(diào)差異得到了很好的增強(qiáng)，青灰色長(zhǎng)石石英砂巖與灰白色塊狀大理巖之間的界線(xiàn)更加清晰。結(jié)合地質(zhì)資料可知，粉紅色調(diào)影像為青灰色長(zhǎng)石石英砂巖，綠色調(diào)影像為灰白色塊狀大理巖。

4 基于空間特征的巖性信息增強(qiáng)

紋理信息對(duì)于巖性的識(shí)別至關(guān)重要，它可以宏觀地反映大面積出露的某一種地物，是解譯細(xì)小地物(特別是巖性、植被)的重要標(biāo)志[12]。巖石的紋理信息不受其表面色調(diào)的影響，它反映的是巖石表面結(jié)構(gòu)及其粗細(xì)程度，能夠反映出如地質(zhì)構(gòu)造、巖性界線(xiàn)等特征。遙感技術(shù)對(duì)于巖性的識(shí)別主要依據(jù)對(duì)影像光譜信息的分析，但如果僅用光譜特征進(jìn)行識(shí)別，因異物同譜、同物異譜及表層淺覆蓋等因素引起的巖石影像特征復(fù)雜、色彩反差較弱、細(xì)節(jié)信息不夠豐富等一系列問(wèn)題就難以解決，解譯結(jié)果的準(zhǔn)確性也會(huì)受到很大影響。為了解決這類(lèi)問(wèn)題，提高解譯結(jié)果的可信度，還需對(duì)影像紋理信息進(jìn)行分析[25]。

本文使用的WorldView2數(shù)據(jù)多光譜波段的空間分辨率為2 m，擁有豐富的細(xì)節(jié)信息，可充分利用這一特點(diǎn)來(lái)增強(qiáng)巖性的紋理及邊緣信息。本文通過(guò)對(duì)B4(R)B7(G)B1(B)假彩色合成圖像(圖8(a))進(jìn)行高通濾波處理[25-26]，在保持圖像高頻信息的同時(shí)，消除了圖像中的低頻成分，從而增強(qiáng)了紋理和邊緣等信息。通過(guò)紋理的一致性及規(guī)律性可以有效地排除由異物同譜、同物異譜、風(fēng)化因素及地表覆蓋物對(duì)影像色彩的干擾，將不同的巖性區(qū)域歸并，準(zhǔn)確圈定巖性界線(xiàn)(圖8(b))。

(a) Band4(R)7(G)1(B)組合 (b) 高通濾波紋理增強(qiáng)

圖8 高通濾波增強(qiáng)假彩色合成圖像

Fig.8 False color composition images after high-pass filtering enhancement

從圖8可以看出，經(jīng)高通濾波處理后，不同巖性之間的紋理及邊界信息差異更加明顯(其中灰白色鈣質(zhì)板巖夾青灰色泥巖、石英巖的紋理呈點(diǎn)狀分布，灰白色塊狀大理巖的紋理呈線(xiàn)狀分布)，對(duì)識(shí)別巖性起到重要作用。

5 巖性解譯結(jié)果

采用上述方法對(duì)研究區(qū)WorldView2圖像進(jìn)行1∶25 000比例尺的巖性解譯，結(jié)果如圖9(a)所示。

(a) 1∶25 000比例尺巖性解譯圖 (b) 1∶50 000比例尺地質(zhì)圖

圖9 研究區(qū)局部巖性解譯圖與地質(zhì)圖對(duì)比

與研究區(qū)相同區(qū)域的地質(zhì)圖(圖9(b))對(duì)比可以看出，通過(guò)巖性信息增強(qiáng)處理，巖性得到了細(xì)分，劃分出了更多的巖性類(lèi)別； 而以往的地質(zhì)填圖工作常常受到野外工作條件的限制，只能對(duì)巖性單元進(jìn)行粗略的劃分。從圖9(a)可以看出，通過(guò)巖性解譯劃分出了Pt1B(大理巖)、Pt1B(大理巖、構(gòu)造角礫巖)、O1-2w(青灰色板巖、粉砂巖夾泥巖)、K2e(磚紅、棕紅色泥巖、泥質(zhì)砂巖、泥質(zhì)砂礫巖)和Qhl(粉砂、泥、細(xì)砂)； 而在相同區(qū)域的1∶50 000比例尺的地質(zhì)圖(圖9(b))中，只劃分出了Pt1B(大理巖)、O1-2w(青灰色板巖、粉砂巖夾泥巖)和K2e(磚紅、棕紅色泥巖、泥質(zhì)砂巖、泥質(zhì)砂礫巖)。因此，巖性信息增強(qiáng)在巖性解譯中有著至關(guān)重要的作用。

6 結(jié)論

本文選擇內(nèi)蒙古烏索倫山扎嘎烏蘇地區(qū)進(jìn)行巖性遙感信息增強(qiáng)方法研究。在對(duì)研究區(qū)巖石波譜曲線(xiàn)特征、WorldView2數(shù)據(jù)光譜特征及空間特征進(jìn)行分析的基礎(chǔ)上，應(yīng)用了一系列的巖性信息增強(qiáng)方法，較好地解決了因異物同譜、同物異譜、表層淺覆蓋等導(dǎo)致的巖石影像特征復(fù)雜、色彩反差較弱、細(xì)節(jié)信息不夠豐富等一系列問(wèn)題，從而降低了巖性識(shí)別的難度，使巖性界線(xiàn)得到更準(zhǔn)確圈定，巖性種類(lèi)得到細(xì)分，巖性解譯的精度得到提高。得出以下結(jié)論：

1)巖性色調(diào)對(duì)于識(shí)別巖性具有重要作用。本研究基于WorldView2數(shù)據(jù)的光譜特征和研究區(qū)典型巖石波譜特征，分別采用波段組合法、比值法、HSV變換和主成分分析方法對(duì)研究區(qū)進(jìn)行了巖性信息增強(qiáng)。這些方法在不同程度上增強(qiáng)了不同巖性之間的色調(diào)差異，有利于準(zhǔn)確圈定巖性界線(xiàn)，識(shí)別出更多的巖性種類(lèi)，提高了巖性解譯的精度。

2)紋理信息對(duì)于巖性識(shí)別同樣有著不可忽視的作用，特別是對(duì)于WorldView2數(shù)據(jù)，其空間分辨率較高，可突出地物的細(xì)節(jié)信息。本研究采用高通濾波的方法對(duì)巖性的紋理信息進(jìn)行增強(qiáng)，突出了巖性的細(xì)節(jié)特征，有助于對(duì)巖性進(jìn)行更好的識(shí)別。

3)巖性信息增強(qiáng)在巖性解譯中起到了重要作用，但是，這些方法容易因研究區(qū)的改變而對(duì)增強(qiáng)變換后的結(jié)果產(chǎn)生影響，因此，在今后的工作中應(yīng)著重加強(qiáng)這方面的研究。

[1] 唐煥麗,劉凱,艾彬,等.WorldView-2遙感影像融合方法對(duì)比研究[J].北京測(cè)繪,2013(5):1-7. Tang H L,Liu K,Ai B,et al.Comparison analysis between different fusion methods for the case of WorldView-2 images[J].Beijing Surveying and Mapping,2013(5):1-7.

[2] 王曉鵬,楊志強(qiáng),康高峰,等.WorldView-2高分辨率衛(wèi)星數(shù)據(jù)在西昆侖塔什庫(kù)爾干地區(qū)遙感地質(zhì)調(diào)查中的應(yīng)用[J].地質(zhì)找礦論叢,2014,29(3):428-432. Wang X P,Yang Z Q,Kang G F,et al.Application of WorldView-2 data to remote sensing geological survey in Tashkurgan area at west Kunlun[J].Contributions to Geology and Mineral Resources Research,2014,29(3):428-432.

[3] 李鋼柱,羅文華,侯萬(wàn)榮,等.內(nèi)蒙古自治區(qū)烏拉特中旗索倫山地區(qū)扎嘎烏蘇、浩仁呼都格、哈日格那、索倫、沙布特、烏珠爾少布特(K49E010001、K49E010002、K49E010003、K49E010004、K49E010005、K49E010006)幅1∶5萬(wàn)區(qū)域地質(zhì)礦產(chǎn)調(diào)查項(xiàng)目總體設(shè)計(jì)書(shū)[R].北京:中國(guó)人民武裝警察部隊(duì)黃金指揮部,2013. Li G Z,Luo W H,Hou W R,et al.The Overall Design about the 1∶50 000 Regional Geology and Mineral Resources Survey Project in Zhagawusu,Haorenhuduge,Harigena,Suolun,Shabute and Wuzhuershabute of the Wulate Middle Banner Sauron Mountain Areas in Inner Mongolia Autonomous Region[R].Beijing:Chinese People’s Armed Police Force Gold Headquarters,2013.

[4] 李朋武,高銳,管燁,等.內(nèi)蒙古中部索倫-林西縫合帶封閉時(shí)代的古地磁分析[J].吉林大學(xué)學(xué)報(bào):地球科學(xué)版,2006,36(5):744-758. Li P W,Gao R,Guan Y,et al.Palaeomagnetic constraints on the final closure time of Solonker-Linxi suture[J].Journal of Jilin University:Earth Science Edition,2006,36(5):744-758.

[5] 陶繼雄,蘇茂榮,寶音烏力吉,等.內(nèi)蒙古達(dá)爾罕茂明安聯(lián)合旗滿(mǎn)都拉地區(qū)索倫山蛇綠混雜巖的特征及構(gòu)造意義[J].地質(zhì)通報(bào),2004,23(12):1238-1242. Tao J X,Su M R,Baoyin W L J,et al.Characteristics and tectonic significance of the Solon Mountain ophiolitic melange in the Mandula area,Darhan Muminggan,Inner Mongolia[J].Geological Bulletin of China,2004,23(12):1238-1242.

[6] 任夢(mèng)依,陳建平.ASTER與WorldView-2結(jié)合提取巖性信息流程——以西藏物瑪?shù)貐^(qū)為例[J].地質(zhì)學(xué)刊,2013,37(4):585-592. Ren M Y,Chen J P.On process of lithologic information extraction by ASTER and WorldView-2 data:A case study of Wuma area in Tibet[J].Journal of Geology,2013,37(4):585-592.

[7] 金謀順,王輝,張微,等.高分辨率遙感數(shù)據(jù)鐵染異常提取方法及其應(yīng)用[J].國(guó)土資源遙感,2015,27(3):122-127.doi:10.6046/gtzyyg.2015.03.20. Jin M S,Wang H,Zhang W,et al.Method for extraction of ferric contamination anomaly from high-resolution remote sensing data and its applications[J].Remote Sensing for Land and Resources,2015,27(3):122-127.doi:10.6046/gtzyyg.2015.03.20.

[8] 趙英時(shí).遙感應(yīng)用分析原理與方法[M].北京:科學(xué)出版社,2003. Zhao Y S.Analysis Principle and Method of Remote Sensing Applications[M].Beijing:Science Press,2003.

[9] 劉新星,陳建平,曾敏,等.基于多源遙感數(shù)據(jù)的西藏羌多地區(qū)地質(zhì)構(gòu)造解譯[J].國(guó)土資源遙感,2015,27(3):154-160.doi:10.6046/gtzyyg.2015.03.24. Liu X X,Chen J P,Zeng M,et al.Geological structural interpretation of Qiangduo area in Tibet based on multi-source remote sensing data[J].Remote Sensing for Land and Resources,2015,27(3):154-160.doi:10.6046/gtzyyg.2015.03.24.

[10]吳德文,朱谷昌,吳健生,等.青海芒崖地區(qū)巖石光譜特征分析及應(yīng)用[J].國(guó)土資源遙感,2001,13(4):28-34.doi:10.6046/gtzyyg.2001.04.05. Wu D W,Zhu G C,Wu J S,et al.The analysis and application of spectral characteristics of rock samples from Mangya area,Qinghai province[J].Remote Sensing for Land and Resources,2001,13(4):28-34.doi:10.6046/gtzyyg.2001.04.05.

[11]余健,張志,李閔佳,等.基于ASTER遙感影像的西昆侖巖性信息提取方法研究[J].國(guó)土資源遙感,2012,24(1):22-27.doi:10.6046/gtzyyg.2012.01.05. Yu J,Zhang Z,Li M J,et al.The Methodology of lithologic information extraction by using ASTER data in West Kunlun Mountains[J].Remote Sensing for Land and Resources,2012,24(1):22-27.doi:10.6046/gtzyyg.2012.01.05.

[12]田淑芳,詹騫.遙感地質(zhì)學(xué)[M].2版.北京:地質(zhì)出版社,2013. Tian S F,Zhan Q.Remote Sensing of Geology[M].2nd ed.Beijing:Geological Publishing House,2013.

[13]俞樂(lè).多源遙感信息快速處理與巖性信息自動(dòng)提取方法研究[D].杭州:浙江大學(xué),2010. Yu L.Towards Multi-Source Remote Sensing Information Fast Processing and Automatic Lithological Information Extraction[D].Hangzhou:Zhejiang University,2010.

[14]劉建平,趙英時(shí).高光譜遙感數(shù)據(jù)解譯的最佳波段選擇方法研究[J].中國(guó)科學(xué)院研究生院學(xué)報(bào),1999,16(2):153-161. Liu J P,Zhao Y S.Methods on optimal bands selection in hyperspectral remote sensing data interpretation[J].Journal of the Graduate School,Academia Sinica,1999,16(2):153-161.

[15]焦?jié)櫝?秦彥平,張淑云,等.WorldView-2數(shù)據(jù)在沉積巖地區(qū)的遙感巖性增強(qiáng)方法初探——以新疆喀什阿爾塔什地區(qū)為例[J].西北地質(zhì),2014,47(4):277-283. Jiao R C,Qin Y P,Zhang S Y,et al.A preliminary exploration of methods to enhance remote sensing lithology with WorldView-2 data in sedimentary rocks area:A case study in Kashi Altas district of Xinjiang[J].Northwestern Geology,2014,47(4):277-283.

[16]陳玲,梁樹(shù)能,周艷,等.國(guó)產(chǎn)高分衛(wèi)星數(shù)據(jù)在高海拔地區(qū)地質(zhì)調(diào)查中的應(yīng)用潛力分析[J].國(guó)土資源遙感,2015,27(1):140-145.doi:10.6046/gtzyyg.2015.01.22. Chen L,Liang S N,Zhou Y,et al.Potential of applying domestic high-resolution remote sensing data to geological survey in high altitudes[J].Remote Sensing for Land and Resources,2015,27(1):140-145.doi:10.6046/gtzyyg.2015.01.22.

[17]Pournamdari M,Hashim M,Pour A B.Spectral transformation of ASTER and Landsat TM bands for lithological mapping of Soghan ophiolite complex,south Iran[J].Advances in Space Research,2014,54(4):694-709.

[18]周云,符思濤.遙感圖像色彩增強(qiáng)處理方法探討[J].測(cè)繪與空間地理信息,2010,33(4):153-156. Zhou Y,Fu S T.Discussion on the processing method of color enhancement on remote sensing image[J].Geomatics & Spatial Information Technology,2010,33(4):153-156.

[19]曹建芳,陳俊杰,趙青杉.一種改進(jìn)的HSV顏色空間量化方法及其應(yīng)用[J].南京師范大學(xué)學(xué)報(bào):工程技術(shù)版,2014,14(2):68-73. Cao J F,Chen J J,Zhao Q S.An improved method on color space quantization and application[J].Journal of Nanjing Normal University:Engineering and Technology Edition,2014,14(2):68-73.[20]焦竹青,徐保國(guó).HSV變換和同態(tài)濾波的彩色圖像光照補(bǔ)償[J].計(jì)算機(jī)工程與應(yīng)用,2010,46(30):142-144. Jiao Z Q,Xu B G.Color image illumination compensation based on HSV transform and homomorphic filtering[J].Computer Engineering and Applications,2010,46(30):142-144.

[21]Yang J,Cheng Q M.A comparative study of independent component analysis with principal component analysis in geological objects identification,Part I:Simulations[J].Journal of Geochemical Exploration,2015,149:127-135.

[22]Yang J,Cheng Q M.A comparative study of independent component analysis with principal component analysis in geological objects identification.Part II:A case study of Pinghe District,Fujian,China[J].Journal of Geochemical Exploration,2015,149:136-146.

[23]衛(wèi)亞星,王莉雯.遙感圖像增強(qiáng)方法分析[J].測(cè)繪與空間地理信息,2006,29(2):4-7. Wei Y X,Wang L W.Analysis of enhancement methods about satellite images[J].Geomatics & Spatial Information Technology,2006,29(2):4-7.

[24]鄧書(shū)斌.ENVI遙感圖像處理方法[M].北京:科學(xué)出版社,2010. Deng S B.ENVI Remote Sensing Image Processing Methods[M].Beijing:Science Press,2010.

[25]金劍,田淑芳,焦?jié)櫝?等.基于地物光譜分析的WorldView-2數(shù)據(jù)巖性識(shí)別——以新疆烏魯克薩依地區(qū)為例[J].現(xiàn)代地質(zhì),2013,27(2):489-496. Jin J,Tian S F,Jiao R C,et al.Lithology identification with WorldView-2 data based on spectral analysis of surface features:A case study of Wulukesayi District in Xinjiang[J].Geoscience,2013,27(2):489-496.

[26]金劍.多元數(shù)據(jù)在策勒縣玉龍地區(qū)遙感成礦預(yù)測(cè)中的應(yīng)用[D].北京:中國(guó)地質(zhì)大學(xué)(北京),2013. Jin J.Application of Multivariate Data in Remote Sensing Metallogenic Prediction in the Yulong Region of Cele[D].Beijing:China University of Geosciences(Beijing),2013.

(責(zé)任編輯： 邢宇)

Research on lithological information enhancement method based on WorldView2 data： A case study of Zhagawusu district in Inner Mongolia

WANG Pingping， TIAN Shufang

(SchoolofEarthSciencesandResources，ChinaUniversityofGeosciences(Beijing)，Beijing100083，China)

In order to further explore and discuss the application of lithologic information enhancement method in the lithologic remote sensing interpretation, the authors selected the Zhagawusu district of Sauron Mountain in Inner Mongolia as the research area for studying lithologic remote sensing information enhancement method. To tackle the problems such as complex rock image characteristics, weak color contrast and a few details caused by the factor of the same spectra with different objects, the same object with different spectra, and the coverage of the surface of the rock, the authors analyzed the rock spectral curve features, the WorldView2 spectral features and spatial features in the study area. On the basis of geometric correction and image fusion of WorldView2 data preprocessing, the authors used a series of effective methods for the remote sensing lithologic information enhancement based on spectral characteristics and spatial characteristics. Based on these methods, the authors made the lithologic interpretation of the study area. Comparing the results of interpretation with 1∶50 000 geological map in the study area shows that the results of interpretation have more lithologic types, improve the accuracy of the lithologic interpretation and thus provide the basis for a more accurate remote sensing lithologic identification.

WorldView2 data； rock spectral characteristic； remote sensing lithological information enhancement； spectral characteristic； spatial characteristics

10.6046/gtzyyg.2016.04.27

王平平，田淑芳.基于WorldView2數(shù)據(jù)的巖性信息增強(qiáng)方法——以?xún)?nèi)蒙古扎嘎烏蘇地區(qū)為例[J].國(guó)土資源遙感,2016,28(4)：176-184.(Wang P P,Tian S F.Research on lithological information enhancement method based on WorldView2 data：A case study of Zhagawusu district in Inner Mongolia[J].Remote Sensing for Land and Resources,2016,28(4)：176-184.)

2015-05-15；

2015-06-21

中國(guó)地質(zhì)調(diào)查局地質(zhì)調(diào)查項(xiàng)目“西部重點(diǎn)地區(qū)地質(zhì)礦產(chǎn)調(diào)查評(píng)價(jià)”(編號(hào)： 12120113071800)資助。

TP 753; P 627

A

1001-070X(2016)04-0176-09

王平平(1990-),女,碩士研究生,主要研究方向?yàn)檫b感地質(zhì)應(yīng)用。Email： 770697872@qq.com。

田淑芳(1963-),女,副教授,主要從事遙感和GIS方面的教學(xué)和科研工作。Email： sftian@cugb.edu.cn。