高光譜遙感數(shù)據(jù)的處理及其在鈾資源勘查的應(yīng)用

—以廣西苗兒山地區(qū)為例

劉德長 謝紅接 趙英俊 張進(jìn)業(yè)

0引言

高光譜遙感數(shù)據(jù)是80年代以來發(fā)展起來的新一代(第三代)遙感數(shù)據(jù)。它使遙感技術(shù)步入可以同時(shí)獲取地球表面物質(zhì)成分信息和空間分布特征信息的新階段。本文以廣西苗兒山地區(qū)機(jī)載成像細(xì)分紅外高光譜數(shù)據(jù)為例，探討了高光譜遙感圖像的預(yù)處理和深化處理，野外實(shí)測光譜曲線庫建立與典型地物的波譜特征，以及高光譜數(shù)據(jù)在鈾資源勘查中的應(yīng)用等。

1 高光譜遙感數(shù)據(jù)的預(yù)處理和深化處理

1.1 苗兒山地區(qū)高光譜遙感數(shù)據(jù)特點(diǎn)

上海技術(shù)物理所航空紅外細(xì)分光譜掃描儀是以地質(zhì)遙感為主要應(yīng)用目的開發(fā)的，在2～2.5um光譜區(qū)間細(xì)分為12個(gè)工作波段。廣西苗兒山地區(qū)高光譜遙感數(shù)據(jù)，由中南地質(zhì)局230所和上海技術(shù)物理所于1988年11月2日、4日獲取。平臺(tái)為伊爾14型飛機(jī)，航速為250km/h，航高（絕對航高）達(dá)4500m,沿NE30。方向飛行。所選8個(gè)通道的中心波長分別為1.600μm、2.143μm、2.200μm、2.250μm、2.300μm、2.330μm、2.370μm和2.450μm。通過分析，苗兒山地區(qū)高光譜遙感數(shù)據(jù)有如下特點(diǎn)：

（1）由于儀器本身的原因，所有數(shù)據(jù)沿航向具周期性條帶噪聲；

（2）與TM和MSS等圖像不同，高光譜圖像在同一航帶的各波段之間，沿航向和掃描方向有幾個(gè)至十幾個(gè)象元之位移；

（3）由于航向分辯率（8.68m）和掃描方向分辯率（27m）的不一致，象元呈長方形，其長（航向）寬（掃描向）比為27:8.68=3.11:1;

（4）由于飛行方向、掃描角度、太陽方位角和太陽高度的幾何關(guān)系的變化和大氣的影響，掃描圖像存在邊緣輻射差異；

（5）由于飛行速度不等和掃描角度的變化引起圖像幾何差異；

（6）第7航帶數(shù)據(jù)由于記錄方式的錯(cuò)亂，引起圖像與正確方式所成的圖像呈鏡向?qū)ΨQ關(guān)系；

（7）高光譜圖像的一個(gè)顯著特點(diǎn)是能夠從圖像中提取任意象元的光譜率曲線，但由于成像光譜儀獲取的圖像象元亮度值并不能代表地物的反射率值，需要進(jìn)行反射率轉(zhuǎn)換。

1.2 高光譜數(shù)據(jù)的預(yù)處理

在研究上述所有特點(diǎn)的基礎(chǔ)上，開發(fā)了一系列數(shù)據(jù)預(yù)處理軟件，并進(jìn)行了條帶去除、波段間配準(zhǔn)、航向壓縮、鏡向變換、輻射校正、相對反射率轉(zhuǎn)換等預(yù)處理，最終得到清晰的圖像。

1.2.1去條帶方法研究

（1）原始數(shù)據(jù)分析 原始數(shù)據(jù)存在嚴(yán)重的條帶干擾，垂直于航線方向分布，且從第一通道到第八通道條帶越來越增強(qiáng)。經(jīng)分析，這種條帶是由于在成像過程中機(jī)械或光學(xué)部件不穩(wěn)定造成的，具有明暗交替的特征，并具有一定的周期性，以9個(gè)象素為周期，呈正弦波狀，部分以8個(gè)象素為周期。

（2）干擾條帶的消除方法 可以將干擾條帶作為加性噪聲處理以消除條帶，要求盡量減少信息損失，一般有兩種處理方法：

①平滑法，可以采用9點(diǎn)平滑或9×9點(diǎn)，但9×9窗口平滑效果不好，在航線方向上以9點(diǎn)平滑效果好，但也會(huì)損失有效信息。

②濾波法，從前面的分析可知干擾條帶的周期為9，通過不斷自動(dòng)調(diào)整波形振幅，可以有效地消除干擾，研究過程中使用了濾波法。

（3）效果分析圖 圖1左圖為原始圖像，右圖為去條帶后的圖像，對比左右兩側(cè)的圖像可以明顯地看出干擾條帶被有效地消除了，而原始數(shù)據(jù)信息的損失比較小。條帶消除后圖像的列數(shù)保持不變，但行數(shù)減少，如第5航帶行數(shù)由5712減為5632。

圖1去條帶影像效果對比圖

（左側(cè)為原圖像，右側(cè)為處理后圖像）

1.2.2 波段之間配準(zhǔn)

一般情況下，TM圖像的波段與波段之間都是完全重合的，不存在象元位移，而苗兒山高光譜遙感圖像存在著明顯的象元位移，一般為幾個(gè)→十幾個(gè)象元。從5、6、7、9條航帶來看，都是第2波段（2.143μm）和第3波段（2.2μm）之間象元位移量最小，如第5航帶的第2，3波段則完全重合，這可能說明該儀器的第2，3通道工作狀況最好。下面以第9航帶為例來說明波段之間的位移大小及配準(zhǔn)方法。以第2波段為基準(zhǔn)，其他波段則向左右或上下位移。

以第2波段為基準(zhǔn)

第1波段 右移1個(gè) 下移 7個(gè)

第3波段 左移2個(gè) 上移4個(gè)

第4波段 左移4個(gè) 下移14個(gè)

第5波段 左移4個(gè) 下移7個(gè)

第6波段 左移5個(gè) 下移8個(gè)

第7波段 右移2個(gè) 下移3個(gè)

根據(jù)以上位移量截取完全重合的圖像部分如表1，得到第9航帶最高配準(zhǔn)后圖像大小為列240，行5738。

表1 第9航帶各波段截取部分

	列數(shù)	行數(shù)
第1波段	9+2：256-6	8：6125-11
第2波段	9+3：256-5	15：6125-4
第3波段	9+5：256-3	19：6125-0
第4波段	9+4：256-4	1：6125-18
第5波段	9+7：256-1	8：6125-11
第6波段	9+8：256-17	7：6125-12
第7波段	9+1：256-17	13：6125-7

（說明列數(shù)的9是由于每幅圖像都有9象元的亮邊，所以都被裁掉）

1.2.3 航向壓縮和掃描方向拉伸處理

機(jī)載成像光譜儀沿飛行方向的分辯率由飛機(jī)的飛行速度決定，其速度為250公以/小時(shí)（即250000m/3600s），掃描率為每秒8線，2條線之間的距離為25000/3600/8=8.68，即航向分辨率為8.68m。掃描方向分辨率可由瞬時(shí)視場（6毫弧度）和飛行高度（4500米）給出，即（6/1000）×4500=27m。由此可見，圖像上每個(gè)象元不是呈正方形，而是呈長方形。長、寬分別為航向和掃描方向，長：寬=27：8.68≈3.11:1。從視覺效果來看，顯示出沿航向拉長或沿掃描方向壓縮（扁）的現(xiàn)象。通地沿航向的壓縮或沿掃描方向的拉伸處理（數(shù)據(jù)重采樣方法用三次樣條法），可得到27×27或8.68×8.68的正方形象元。 1.2.4 邊緣輻射校正處理

由于飛行方向、掃描角度、太陽方位角和太陽高度的幾何關(guān)系的變化及大氣的影像，掃描圖像存在邊緣輻射差異。這種差異可以通過幾何光學(xué)模型和大氣模型予以校正，但由于模型十分復(fù)雜，影響因素很多不易操作，因此，通常采用統(tǒng)計(jì)方法進(jìn)行圖像的邊緣輻射校正。設(shè)有圖像矩陣Anm

其中n為行數(shù)（沿飛行方向），m為列數(shù)（沿掃描方向），對圖像每列求均值得到向量：

反映了整個(gè)圖像邊緣輻射差異的平均趨勢，作者以圖像中機(jī)下點(diǎn)的輻射量為標(biāo)準(zhǔn)對其余各列像元進(jìn)行歸一化處理，得差值向量：

校正后得圖像為：

1.2.5 正切校正處理

由于掃描角度的變化，圖像中每個(gè)象元對應(yīng)的地面分辨率隨著變化，這種幾何差異的校正，通常使用正切校正方法。

h NP1 P1 NP2 P2 O 掃描角度變化如圖2：

圖2 掃描角度變化示意圖

其中0為傳感器，NP1，NP2為實(shí)際掃描面，P1、P2為地面。傳感器面向正下方等角速地收集地面數(shù)據(jù)， 為瞬時(shí)視場角，h為傳感器與地面之間的距離。

正切校正算法如下：

設(shè)圖像象元為等間隔網(wǎng)格上（即正方形象元）的點(diǎn)，則在掃描方向圖像網(wǎng)格，寬度h× 。從機(jī)下點(diǎn)開始計(jì)算，第i點(diǎn)網(wǎng)格與機(jī)下點(diǎn)的距離為i×h× 。實(shí)際上，地面網(wǎng)格點(diǎn)與機(jī)下點(diǎn)的距離為h×tan（i× ）。作者按照網(wǎng)格寬度h× 對圖像進(jìn)行重采樣，設(shè)重采樣后圖像網(wǎng)格點(diǎn)i對應(yīng)圖像網(wǎng)格點(diǎn)K，則有：

i·h·＝h·tan（k·）

求出K：

通常K不會(huì)落在網(wǎng)格點(diǎn)上，即K不是整數(shù)，可以得出i，滿足：

≤K≤ ＋1

采用距離倒數(shù)插值求出K點(diǎn)的象元值：

校正后的圖像象元值DN（i）=DN（K）

顯然通過上述方法校正后，圖像的列數(shù)將增加（圖3）。

圖3 正切校正前后對比影像圖

（左圖為原圖像；右圖為校正后圖像）

1.2.6 鏡像變換處理

從己經(jīng)處理過的4條航帶數(shù)據(jù)來看,只有第7航帶記錄錯(cuò)亂,顯示出來的圖像與實(shí)際圖像呈鏡向?qū)��?yīng)關(guān)系�？捎孟率奖磉_(dá)：

f（m，n）= f（M-m，n）

m=0,1,…,M-1； n=0,1,…,N-1

其中m為列數(shù)，n為行數(shù)，M、N為該圖像總的列數(shù)和總的代數(shù)，插值方法用的是三次樣條法（圖4）。

圖4 鏡像轉(zhuǎn)換前后對比影像圖

（左圖為原圖像：右圖為轉(zhuǎn)換后圖像）

1.2.7 相對反射率轉(zhuǎn)換處理

高光譜遙感圖像不同于多光譜遙感圖像的一個(gè)根本特點(diǎn)在于從高光譜圖像立方體里能夠提取出精確的地物光譜反射率曲線。據(jù)此，可以直接識(shí)別地物。因此，準(zhǔn)確地提取影像中的光譜信息是圖像分析的基礎(chǔ)，也是成像光譜技術(shù)研究的一項(xiàng)重要內(nèi)容。

成像光譜儀在0.4-2.5μm工作波長范圍內(nèi)主要光譜能量來自地球表面對太陽光的反射。由于太陽光受大氣吸收影響，太陽光的地面輻照度隨波長變化大。因此，成像光譜儀所獲得的圖像DN值不能代表地物的反射率值，需要進(jìn)行反射率轉(zhuǎn)換。

采用了2種方法進(jìn)行相對反射率轉(zhuǎn)換，下面介紹其原理：

（1）統(tǒng)計(jì)法 這是針對沒有進(jìn)行同步地面反射率測量采用的一種方便、快捷的方法，常用的有對數(shù)殘差轉(zhuǎn)換法和內(nèi)在平均相對反射率法（IARR）。研究過程編制了對數(shù)殘差轉(zhuǎn)換法軟件。

設(shè)有高光譜圖像DN（I,j,b）,i=1,2,…M為行數(shù)，j=1,2….N為列數(shù)，b=1,2,…k為波段數(shù)。

在i行j列有光譜量均值：

波段b的灰度均值：

現(xiàn)假設(shè)：

①整個(gè)圖像區(qū)域中各類成份足夠雜亂；

②對征一波段的整個(gè)區(qū)域取平均，所得的光譜與地面吸收特征無關(guān)，則可用對數(shù)殘差轉(zhuǎn)換，得相對反射率；

logRy＝logDN（i,j,b）一logSij一logDb

（2）地面定標(biāo)法 在飛行的同時(shí)，用野外光譜儀實(shí)地測量典型地物的光譜，與成像光譜圖像中對應(yīng)地點(diǎn)象元的DN值相對照，求出圖像DN值到地面反射率的轉(zhuǎn)換公式。理論與實(shí)踐都表明圖像DN值與地面反射率有以下近似的關(guān)系：

其中DN為波段b給定象元的數(shù)值，p（λ）為波段b的波長λ對應(yīng)的地面反射率，Ab為乘積項(xiàng)，反映了大氣傳輸及儀器設(shè)備的放大比；Bb為偏移項(xiàng)，反映了大氣輻射值及儀器的零點(diǎn)偏移。Ab、Bb可以通過地面定標(biāo)數(shù)據(jù)及相應(yīng)位置的圖像DN值來確定。

設(shè)有N類地面定標(biāo)點(diǎn)的反射率，則有一組方程：

 i=1，2,…，n

求此方程組的最小二乘方解得：

本次使用的數(shù)據(jù)是1988年11月初的飛行數(shù)據(jù)，但當(dāng)時(shí)沒有做同步地面測量。作者于1995年11月初進(jìn)行了補(bǔ)充測量，一方面在于了解工作區(qū)的典型地物的光譜特征，另一方面也想做一次似同步野外定標(biāo)實(shí)驗(yàn)。全部測量數(shù)據(jù)于1995.11.4-11日獲得，這期間天氣晴朗，采樣時(shí)間11-16點(diǎn)，大體與1988.11.2、4日的時(shí)相相當(dāng)。以第五航帶為例，當(dāng)時(shí)選取了資江水，張家礦區(qū)7號帶、煙竹河及公路交叉處的路邊白沙及附近正�；◢弾r作為定標(biāo)點(diǎn)。

表2 野外定標(biāo)點(diǎn)及其圖像坐標(biāo)表

表3 各點(diǎn)在對應(yīng)波段上的反射率值（%）和圖像DN值（0-255）表

這樣依此3點(diǎn)可以得到一條回歸曲線,可將圖像上任一點(diǎn)的DN值轉(zhuǎn)化為相對反射率值,并編制了相應(yīng)的軟件。上述3點(diǎn)轉(zhuǎn)換后的相對反射率值（表4）如下：

表4 轉(zhuǎn)換后的相對反射率值表

坐標(biāo)點(diǎn)	1600 nm	2143 nm	2200 nm	2250 nm	2300 nm	2330 nm	2450 nm
1061 1153	3.568	3.278	4.441	4.092	2.828	0.050	7.351
981 956	69	65	67	66	54	54	77
979 957	31	29	22	30	39	39	32

2. 波譜特征

野外地物波譜測量是遙感應(yīng)用的基礎(chǔ)，而高光譜分辯率數(shù)據(jù)具有圖譜合一的特征，因此波譜測量對航空高光譜數(shù)據(jù)的處理和解釋尤其重要。一般情況下，在飛行時(shí)應(yīng)同步進(jìn)行野外地物波譜測量，為野外定標(biāo)及數(shù)據(jù)處理和解釋提供依據(jù)。但由于條件限制，如前所述飛行時(shí)沒有同步進(jìn)行野外地物波譜測量。為了彌補(bǔ)這方面的不足，研究過程中擬同步性地補(bǔ)充進(jìn)行了野外地物波譜測量。在EVNI軟件上建立了野外實(shí)測光譜曲線庫，并研究了該區(qū)典型地物的波譜特征。

2.1 地物波譜測量及影響巖石反射光譜特征的主要因素

2.1.1 野外地物譜測量

利用美國地球物理環(huán)境研究公司（GER）研制的IRIS紅外智能波譜儀對苗兒山地區(qū)的典型地物進(jìn)行了波譜測試，共獲取了150多條波譜曲線，基本上反映了研究區(qū)各典型地物的波譜特征，并為野外定標(biāo)提供了基本數(shù)據(jù)。

野外光譜測量是在1995年11月4日至11月11日進(jìn)行的。如前所述，這個(gè)時(shí)間段與當(dāng)年獲取數(shù)據(jù)的時(shí)間相當(dāng)，便于地面所測光譜信息與航測的光譜信息做對應(yīng)分析。在測錄中用GPS在1：5的地形圖上做了較準(zhǔn)確的定位，并且對每一個(gè)測點(diǎn)的地物類型做了細(xì)致的描述，使以地面波譜測量為基礎(chǔ)的處理比較可靠。為了分析測得的光譜曲線，先分析一下影響巖石反射光譜特征的主要因素。

2.1.2 影響巖石反射光譜特征的主要因素

⑴巖石成分的影響 由于巖石的成因不同，巖石的成分也隨之發(fā)生變化，巖石呈現(xiàn)出不同顏色。例如，酸性中花崗巖以長英質(zhì)礦物為主，巖石呈淺色調(diào)；基性巖中的輝長巖以鐵鎂質(zhì)礦物為主，巖石呈暗色調(diào)。前者的光譜反射率高于后者。從基性巖到酸性巖，它們之間的反射率可在70%范圍內(nèi)變化，深色巖石的反射率要低于淺色巖�？梢妿r石光譜的反射率的高低與其顏色有關(guān)，更確切地說應(yīng)與巖石的成份有關(guān)。

⑵巖石結(jié)構(gòu)、構(gòu)造的影響 巖石在生成過程中由于物理化學(xué)條件不同、生成不同程度、不同結(jié)晶度的巖石。巖石的結(jié)構(gòu)、構(gòu)造反映巖石中不同礦物顆粒大小及它們之間的關(guān)系。通常，巖石的結(jié)晶度和粒度不同造成同類巖石光譜反射率的不均勻性，如粗�；◢弾r的光譜反射率變化要大于細(xì)�；◢弾r。

⑶蝕變和風(fēng)化作用的影響 在內(nèi)營力和外營力作用下，巖石發(fā)生蝕變或風(fēng)化，生成蝕變和風(fēng)化礦物，并引起巖石結(jié)構(gòu)、構(gòu)造和色調(diào)等的變化。

2.2 苗兒山地區(qū)典型地物波譜庫的建立及地物特征分析

2.2.1 典型地物波譜庫的建立

從野外獲取的150多條光譜曲線按類別可分成11個(gè)大類，并以此建成11個(gè)庫。它們分別是正常花崗巖、花崗巖風(fēng)化而成的沙土、花崗巖弱蝕變破碎物、花崗巖強(qiáng)水云母化、花崗巖與草混合物、礦石、硅化帶、稻田、水體、各類植被、混合植被等。

為了使實(shí)測曲線（0.4-2.5um）能與從高光譜圖像立方體中提取的曲線（1.6-2.45um）進(jìn)個(gè)對比分析和匹配研究,須將原曲線按高光譜圖像的七個(gè)波段及相應(yīng)的帶寬（100nm或50nm）進(jìn)行重采樣。 2.2.2 地物波譜特征分析

根據(jù)建成的庫，對其波譜特征進(jìn)行了分析。取各庫的平均光譜曲線，如S1、S2、S3、S4，分別代表正�；◢弾r、花崗巖風(fēng)化而成的砂土、花崗巖弱蝕變破碎物、花崗巖強(qiáng)水云母化等。

S1代表正�；◢弾r 總體反射率較低，在20%左右，曲線較平坦。

S2代表風(fēng)化后成的沙土 總體反射率在40-500%，在2200處有較明顯的吸收峰。

S3弱蝕變花崗巖 總體反射率比S2稍低，在2200nm處有較明顯的吸收峰。

S4強(qiáng)蝕變花崗巖 總體反射率明顯升高，特別是在低波段區(qū)間反射率升高更為明顯，在2200nm處有很明顯的吸收峰，880nm處有較明顯的吸收峰。

S5花崗巖型鈾礦石 總體反射率同正�；◢弾rS1，在20%左右，但不同的是在2200nm處有明顯吸收峰。

S6花崗巖與各種草混合后的曲線 在670-760nm之間有明顯的陡坡，反映植被的典型特征，即所謂的“葉綠素紅邊”。植被所占比例越多，這條斜邊越長，反之，越短。

S7硅化帶的特征曲線 總體反射率<40%，在670nm、2200nm有吸收峰（圖5-6）。

圖5 硅化帶的光譜特征曲線

圖6 重采樣后硅化帶的特征曲線

S8稻田反射率特征曲線 其中，S81與植被曲線（S10）有非常近似的特點(diǎn)，在670-760nm之間有明顯的陡坡，在1500nm處有明顯吸收峰，在1650nm和2250nm處有明顯反射峰，反映了田中稻草梗的反射率特點(diǎn)；S82反映的主要是田中泥土的反射率特征，由于它們源于花崗巖風(fēng)化破碎物，故近似于花崗巖特征曲線。

S9水體 資江水，梅溪水的特征曲線，反射率低。

S10各種植被及混合植被的特征 由于葉綠素在670nm紅波段具強(qiáng)吸收，在760nm近紅外波段具高反射的特點(diǎn)，便得各種植被在670-760nm之間具有明顯的陡坡，即所謂的“葉綠素紅邊”。它是健康植被的一個(gè)特征，如果植被或農(nóng)作物由于受到如重金屬、碳?xì)浠�合物、干旱、病蟲害、缺氮等因素引起脅迫（stress），則會(huì)使這些植被或農(nóng)作物的特征曲線產(chǎn)生微小的變化，即所謂的“紅邊藍(lán)移”現(xiàn)象。利用“紅邊藍(lán)移”現(xiàn)象可以預(yù)測病蟲害等的發(fā)生，并及時(shí)采取補(bǔ)救措施。本區(qū)的某些植物，如乍子樹、竹葉、松樹、雜草等都具有此種現(xiàn)象（圖7），說明上述植物受到脅迫。

圖 7 植物受到脅迫后的光譜特征曲線

（上圖為乍子樹、竹葉、松樹；下圖為雜草）

3 高光譜數(shù)據(jù)的地質(zhì)應(yīng)用及效果

3.1 高光譜分辨率圖像上硅化帶的識(shí)別及效果

高光譜分辨率數(shù)據(jù)的一個(gè)最重要用途是識(shí)別蝕變礦物。盡管由于原始數(shù)據(jù)和研究區(qū)植被覆蓋方面的原因，使應(yīng)用效果受到了一定的限制，但在硅化帶的識(shí)別方面仍取得了成功。

圖8為第5航帶第2波段（中心波長2.143μm）圖像（已經(jīng)各種處理和精校正），從該圖像上可以清楚地看到許多近東西向的暗色調(diào)線性體，這些線性體在對應(yīng)區(qū)域的TM圖像上沒有或僅略有顯示，在原有的大比例尺地質(zhì)圖上也只有少數(shù)幾條被填繪出來。通過對比圖像上直接提取的已知帶（硅化斷裂帶）和未知帶的光譜曲線，發(fā)現(xiàn)兩者具有一致的曲線特征。因此，推測這些未知的線性體可能是硅化斷裂帶。

作者對分析結(jié)果進(jìn)行了野外檢驗(yàn)，證實(shí)它們大多數(shù)確系硅化斷裂帶，其中馬家南邊的那條硅化帶規(guī)模很大，寬10～20m，主要為白色石英脈帶（鈾含量為26 x10－6～45義10—6），帶內(nèi)有殘留碎裂花崗巖（鈾含量40x10－6，具紅化、高嶺土化和弱硅化），圍巖為碎裂花崗巖（鈾含量為60 x10－6一70 x10－6）。硅化帶內(nèi)殘留碎裂花崗巖的鈾有丟失，而在圍巖中有富集現(xiàn)象。在豆詐山巖體及其周圍的高光譜圖中，解譯出的線性構(gòu)造包括了該區(qū)地質(zhì)圖上填出的主要硅化斷裂帶，其中有些硅化斷裂帶直接控制了鈾礦化。

圖8 第5航第2波段（2.143μm）圖像

3.2 高光譜分辨率圖像與其他多源信息的綜合分析及效果

將高光譜數(shù)據(jù)、TM數(shù)據(jù)與已有地質(zhì)圖等多源信息進(jìn)行融合，可形成一幅集多種信息于一體的新類型圖像（圖9）。

圖9 綜合解譯圖

在該圖像上，地質(zhì)圖上已有的主要斷裂構(gòu)造和硅化帶都已被解譯出來，同時(shí)，還解譯出很多線狀構(gòu)造和環(huán)狀構(gòu)造，而這些構(gòu)造在地質(zhì)圖上是沒有的。如fl—沙子江-茶坪斷裂構(gòu)造，f2一斷裂構(gòu)造，f3一斷裂構(gòu)造，f4一牛大江硅化帶等；環(huán)狀構(gòu)造有4個(gè)，分別以C1、C2、C3、C4表示。 圖9 綜合解譯圖 其中以Cl規(guī)模最大。本區(qū)是巖漿巖分布區(qū)，環(huán)狀構(gòu)造的形成一般與巖漿的侵人活動(dòng)有關(guān)。4個(gè)環(huán)狀構(gòu)造都發(fā)育在兩期巖體的交界部位，其中C2和C3是豆詐山巖體與香草坪巖體之接觸部位，C4是牛塘坳巖體（γ53－1）、香草坪巖體和苗兒山巖體之接觸部位，可能為更晚期的巖漿侵人活動(dòng)所致。C1不僅規(guī)模最大，且結(jié)構(gòu)復(fù)雜，該環(huán)形構(gòu)造西南部處于豆詐山巖體與香草坪巖體之接觸帶；東北部受香草坪巖體、苗兒山巖體和牛塘傲巖體三者的影響；東南部基本上是香草坪巖體和苗兒山巖體的分界線；西北部分還同時(shí)與小環(huán)C3相交。此環(huán)形構(gòu)造還是本區(qū)四條大硅化斷裂帶（F1, F2, F31, F32）的夾持部位。因此，該環(huán)形構(gòu)造區(qū)具有很好的構(gòu)造、鈾源（γ51、γ52－2）和成礦的先決條件，雙滑江礦床和白毛沖礦床就在此環(huán)形構(gòu)造范圍內(nèi)。前者產(chǎn)于該環(huán)形構(gòu)造的環(huán)邊上。本次研究預(yù)測的一處鈾礦化遠(yuǎn)景區(qū)（P1）就在該環(huán)的西北部分。由此可以認(rèn)為，該環(huán)形構(gòu)造的展布范圍應(yīng)該是很有成礦前景的地區(qū)。

經(jīng)過上述各種處理及綜合解譯，結(jié)合前人對本區(qū)成礦機(jī)理的認(rèn)識(shí)，本次研究預(yù)測了三個(gè)成礦有利區(qū)，通過驗(yàn)證，其中的兩個(gè)地區(qū)放射性測量值偏高，并有鈾礦化顯示。

參考文獻(xiàn)

1．鄭蘭芬 王晉年，1992，成像光譜遙感技術(shù)及成像光譜信息提取的分析與研究，《環(huán)境遙感》，Vol. 7，No. 10。

2．王晉年鄭蘭芬，1996，成像光譜圖像光譜吸收鑒別模型與礦物填圖研究，《環(huán)境遙感》》，Vol．11，No. 10。

3．王向軍，1995，成像光譜信息可視化表達(dá)系統(tǒng)，中國科學(xué)院遙感應(yīng)用研究所碩士學(xué)位論文。 4．王欽敏，1995，衛(wèi)星遙感技術(shù)發(fā)展和應(yīng)用動(dòng)態(tài)，《遙感信息》，Vol．40, No. 4。

5．R. J. P. Lyon, 1996，風(fēng)化及其它類荒漠漆表面層對高光譜分辨率反射率的影響（一），《環(huán)境遙感》，Vol. 11，No. 2。

6．Joseph C.H, Chein-I Chang, 1994, Hyperspectral Image，Classification and Dimenesionality Reduction：An Orthogonal Subspace Projection Approach,《IEEE Transaction on Geoscience and Remote Sensing》，Vo1.32, No.4（779-784）。

7．Geotz, A.，and Herring, M:，1989, The High Resolusion Imaging Spectrometer （HIRIS）for EOS，《IEEE Transaction on Geoscience and Remote Sensing》，Vol.27, No.2。

8．Mazer, A.S.,et al, 1988, Image Processing Software for Imaging Spectrometry Data Analysis.《Remote Sensing of Environment》,Vol.24，No. l（201-210）。

9．William H.F.，et al，1994, Retrieval of Apparent Surface Reflectance from AVRIS data：A Comparison of Empirical line, Radiative Transfer, and Spectral Mixture Methods,《Remote Sensing of Environment》》，Vo1.47（311-321）。

10．ENVI User's Guide, 1996,Research Systems, Inc.,USA

11．A Hyperspectral Remote Sensing Primer, 1996, NASA, USA。