【正文】 而定,DN格值不變，僅移動(dòng)坐標(biāo)零點(diǎn)。Ⅰ、7通道與1通道“AND”后放到8通道中，這樣就把原先的“去邊框”數(shù)據(jù)沖掉。Ⅱ、同理做6通道的掩模與“AND”處理，所得數(shù)據(jù)分別放在1113通道中。②準(zhǔn)歸一化校正從頭文件信息中讀取各波段屬低增益還是高增益狀態(tài)，查取列表得到對(duì)應(yīng)波段的GK 值，代入公式：K=GKd2/sinG，即可得到各波段的準(zhǔn)準(zhǔn)歸一化系數(shù)，如下：表2 各波段的準(zhǔn)歸一校正系數(shù)波段123457GK各拉伸后的波段對(duì)應(yīng)乘以準(zhǔn)準(zhǔn)歸一化系數(shù)。Ⅰ、如下圖彈出“ARI Image Arithmetic”對(duì)話框，“CNST”參數(shù)處輸入K值，“RVAL”的值寫“0，255，0，255”，然后點(diǎn)擊 “Status” “Run”。Ⅱ、同理做1113通道的準(zhǔn)準(zhǔn)歸一化校正，輸出的數(shù)據(jù)分別存在1113通道中。（注）當(dāng)因某種原因(頭文件缺失,或僅須做單景異常提取)不能或無(wú)須做準(zhǔn)準(zhǔn)歸一化處理時(shí)，本節(jié)的(準(zhǔn)準(zhǔn)歸一化校正)省略，而(各波段與掩模“AND”同時(shí)做徑輻射校正)改為各波段與掩?！癆ND”同時(shí)做徑輻射校正及無(wú)損失拉伸。同理做1113通道的徑輻射校正及無(wú)損失拉伸，輸出的數(shù)據(jù)分別存在1113通道中。①羥基異常（HA）利用BBB5及B7四個(gè)波段的組合進(jìn)行羥基異常的提取，輸入的應(yīng)是對(duì)應(yīng)波段的準(zhǔn)準(zhǔn)歸一化校正后的數(shù)據(jù)，該示例中是112及13通道組合。點(diǎn)擊“XPACE”按鈕，找到“PCA Principal Analysis”命令，點(diǎn)擊后彈出“PCA Principal Analysis”對(duì)話框，填好參數(shù)。其中：“DBIC”是參與主成分分析的通道（圖示1113）；“EIGN”是輸出哪些主分量，示例中輸出PCPCPCPC4共四個(gè)主分量；“DBOC”是輸出的四個(gè)主分量放于哪些通道（圖示為11114通道）；“DEVRANGE”是標(biāo)準(zhǔn)差的范圍（四倍）；“MASK”是掩模（輸入掩模在哪一“Segment”）。點(diǎn)擊 “Status” “Run”后，在DOS界面出現(xiàn)統(tǒng)計(jì)數(shù)據(jù)，從中可以得到各分量的特征值及分析羥基異常表現(xiàn)在哪個(gè)分量。選擇羥基異常分量的規(guī)定：在特征向量的相關(guān)矩陣中，B5與BB7的符號(hào)相反，而且BB7的值越大越好。分析特征向量的相關(guān)矩陣可以得出羥基異常集中在PC4中。②鐵染異常（FA）利用BBB4及B5四個(gè)波段的組合進(jìn)行鐵染異常的提取，輸入的應(yīng)是對(duì)應(yīng)波段的準(zhǔn)準(zhǔn)歸一化校正后的數(shù)據(jù)，該示例中是11及12通道組合。對(duì)112通道進(jìn)行PCA，得到的各分量分別放于11115通道中（1115已放羥基的PCPCPC3，因羥基異常存于PC4中，所以PCPCPC3的數(shù)據(jù)可以舍棄），其方法與羥基異常的提取方法一樣，統(tǒng)計(jì)的數(shù)據(jù)如下：選擇鐵染異常分量的規(guī)定：B1與B3的符號(hào)相反，由此可以分析出異常也是集中在PC4中。異常的分級(jí)利用標(biāo)準(zhǔn)離差的倍數(shù)Kσ（見附表3），采用不同的Kσ得到的異常范圍就不同。①羥基異常（HA）HA放在14通道中，打開“Edit”下拉菜單中的“Edit LUT”并選擇14通道的直方圖，在“Min X”“Max X”中輸入對(duì)應(yīng)不同的值，則該值就大約對(duì)應(yīng)某個(gè)Kσ，輸入“Min X”“Max X”后點(diǎn)擊“Apply”，顯示窗口即可顯示出異常的分布范圍。一般異常分為三級(jí)，Kσ值最大的為一級(jí)，次之為二級(jí)，最小的為三級(jí)。下面是HA的分級(jí)Kσ取值及異常分布情況： Ⅰ 一級(jí)異常。 Ⅱ 二級(jí)異常， Ⅲ 三級(jí)異常，以上幾幅圖是異常分布的顯示效果，下面的方法把分級(jí)的三層數(shù)據(jù)保存在三個(gè)通道中，其中：“RVAL”參數(shù)中的前兩個(gè)數(shù)即是不同Kσ值所對(duì)應(yīng)的上下限。經(jīng)過(guò)上述的分級(jí)處理，三級(jí)異常（HAHAHA3）分別放于1118通道中。②鐵染異常（FA）打開“View”下拉菜單中的“Histograms…”，顯示15通道的直方圖，發(fā)現(xiàn)直方圖的最右邊有一豎條，這可能與干擾有關(guān)，因此先將該層與“MASK”層（7通道）進(jìn)行“AND”處理。與HA分級(jí)的原理一樣，通過(guò)在“Edit LUT”顯示窗口的顯示效果，、得到的FA三級(jí)異常分別放于121通道中。①羥基異常（HA）點(diǎn)擊“FMO Mode Filter”，彈出“FMO Mode Filter”對(duì)話框，其中：“FLSZ”輸入濾波窗口的大?。▓D示為“33”）。分別對(duì)118通道進(jìn)行“33”、“55”濾波，濾波后的數(shù)據(jù)仍放于原通道。②鐵染異常（FA）用羥基異常濾波的方法對(duì)鐵染異常進(jìn)行濾波，19通道進(jìn)行“33”濾波，21通道進(jìn)行“55”濾波，濾波后的數(shù)據(jù)仍放于原通道。TM數(shù)據(jù)常以BBB7為背景圖，為了仍保持原有的地物地貌特征（原有的植被、鹽堿地等），因此需要未經(jīng)干擾預(yù)處理的數(shù)據(jù)，但需要進(jìn)行拉伸以及準(zhǔn)準(zhǔn)歸一化校正，以便增強(qiáng)地物信息，有助于目視解譯。以下三幅圖示是對(duì)原始數(shù)據(jù)進(jìn)行拉伸處理，處理后的數(shù)據(jù)放在2224通道中。 以上是拉伸處理，還需對(duì)拉伸后的數(shù)據(jù)進(jìn)行日地校正，其方法與前面的校正方法一致，在此不再贅述，經(jīng)校正后的數(shù)據(jù)對(duì)應(yīng)存在校正前的通道中，即2224通道中。:25萬(wàn)標(biāo)準(zhǔn)分幅異常圖制作在對(duì)浙江省12景ETM遙感數(shù)據(jù)分別進(jìn)行異常提取的基礎(chǔ)上，本項(xiàng)目還需要編制浙江省1：25萬(wàn)標(biāo)準(zhǔn)羥基、鐵染遙感異常圖15幅。具體工作過(guò)程如下。:25萬(wàn)標(biāo)準(zhǔn)分幅國(guó)家基本地形圖，按統(tǒng)一規(guī)定的經(jīng)差和緯差進(jìn)行分幅，每幅圖的內(nèi)圖廓皆由經(jīng)線和緯線構(gòu)成，并在國(guó)際百萬(wàn)分一地圖分幅編號(hào)的基礎(chǔ)上，建立了各種比例尺地形圖的圖幅編號(hào)系統(tǒng)。下表為浙江省1：25萬(wàn)分幅地圖的基本情況：表：1：25萬(wàn)浙江省分幅圖 (資料來(lái)源于國(guó)家基礎(chǔ)地理信息中心)分幅號(hào)名稱東邊經(jīng)度南邊緯度西邊經(jīng)度北邊緯度G51C001001瑞安縣G51C001004H50C001004常州市H50C002004宣城縣H50C003003景德鎮(zhèn)市H50C003004金華市H50C004003上饒市H50C004004衢州市H51C001001上海市H51C002001杭州市H51C002002定?？hH51C002003定海縣H51C003001嵊縣H51C003002寧波市H51C004001溫州市H51C004002一江山島在ENVI主菜單中，選擇Basic Tools——Mosaicking——Georeferenced，出現(xiàn)Map Based Mosaic窗口，點(diǎn)擊菜單Import——Import Files，出現(xiàn)Mosaic Input Files對(duì)話框，點(diǎn)擊Open，把前面提取好的浙江省ETM遙感異常數(shù)據(jù)中的相鄰2景先加載進(jìn)來(lái)，點(diǎn)擊OK，等大約3分鐘，在Map Based Mosaic窗口中就出現(xiàn)了2景數(shù)據(jù)鑲嵌一起的樣子，這時(shí)可以點(diǎn)擊下方選擇哪景數(shù)據(jù)放置在上面或下面，或者點(diǎn)擊Edit Entry，在出現(xiàn)的對(duì)話框中，可以進(jìn)行鑲嵌的色彩均衡化設(shè)置，和羽化距離的設(shè)置，這些設(shè)置都是影響鑲嵌的效果的。都設(shè)置好以后，點(diǎn)擊OK，然后點(diǎn)擊File——Apply，接下來(lái)就開始對(duì)2景數(shù)據(jù)進(jìn)行鑲嵌了，大概要30分鐘完成。按前面的方法，對(duì)12景數(shù)據(jù)依次進(jìn)行鑲嵌。最后，把浙江省12景ETM數(shù)據(jù)全部鑲嵌在一起。：25萬(wàn)標(biāo)準(zhǔn)分幅 在ENVI主菜單中，選擇Basic Tools——Resize Data——Open——加載12景鑲嵌數(shù)據(jù)——Spatial Subset，出現(xiàn)Select Spatial Subset對(duì)話框，在對(duì)話框中點(diǎn)選map按鈕，在出現(xiàn)的對(duì)話框中輸入1：25萬(wàn)標(biāo)準(zhǔn)分幅的左上角和右下角的經(jīng)緯度，點(diǎn)OK，然后重復(fù)前面的步驟，按標(biāo)準(zhǔn)分幅依次分出15幅異常圖來(lái)。①現(xiàn)在已經(jīng)有15幅異常的標(biāo)準(zhǔn)分幅圖了。接下來(lái)要做的是，首先在ENVI中，把15幅ENVI格式的異常數(shù)據(jù)轉(zhuǎn)換成PCI格式的。File——Save File As——PCI——Open，分別輸入各幅數(shù)據(jù)，點(diǎn)OK，然后確定路徑，再點(diǎn)OK，就生成了各自的PCI格式的文件了。②接著，在PCI中，對(duì)15幅數(shù)據(jù)分別用tif格式輸出以下這些圖，即背景圖、1級(jí)羥基異常、2級(jí)羥基異常、1級(jí)鐵染異常、2級(jí)鐵染異常。點(diǎn)擊ImageWorks模塊，點(diǎn)擊左上角的Use Image File，選擇要處理的分幅異常數(shù)據(jù)，點(diǎn)擊Accept，點(diǎn)擊File——Utility，在彈出的對(duì)話框中，點(diǎn)擊File——Export to…，出現(xiàn)Translate對(duì)話框，在Destination file輸入輸出文件的路徑和名字及格式，這里應(yīng)該按照輸出的5個(gè)層分別輸入，格式為tif；Source layers中點(diǎn)選相應(yīng)的通道，然后點(diǎn)Add加載到右邊，也就是說(shuō)當(dāng)要輸出背景圖的tif時(shí)候，按順序點(diǎn)選1通道到右邊，還有四個(gè)異常層的輸出也是同理。最后點(diǎn)一下左下角的Translate，就完成了相應(yīng)層的輸出工作。③第三步，就是在Photoshop中，把剛才輸出的5個(gè)層疊加起來(lái)。一共有15幅圖，按順序1幅1幅圖分別來(lái)做。首先，，由于此時(shí)的異常圖，背景是黑的，異常為白色亮塊，所以可以按下面方法做。點(diǎn)擊工具欄里的魔棒工具，這樣就選中了所有黑色的部分，然后按Delete去除黑色部分，接著點(diǎn)擊菜單的選擇——反選，這時(shí)候異常部分被選中，再點(diǎn)圖像——模式——RGB顏色，這樣做的目的是使得異常能被附上彩色；然后按照如下規(guī)則：一級(jí)羥基異常采用FF0000顏色值，二級(jí)羥基異常采用FFFF00顏色值，一級(jí)鐵染異常采用0000FF顏色值，二級(jí)鐵染異常采用00FFFF顏色值，用油漆桶工具給異常附上相應(yīng)的顏色。最后，這樣背景圖層在下面，異常圖層在上面。按照這個(gè)做法。至此，一幅標(biāo)準(zhǔn)分幅的遙感異常圖就完成了。然后把剩余的14幅也這樣完成。此項(xiàng)工作在Mapgis軟件中完成。需要完成如下工作：①配上與同比例尺地理底圖相應(yīng)的坐標(biāo)系、注記及圖廓整飾。②根據(jù)技術(shù)要求的規(guī)定進(jìn)行用色，內(nèi)容主要包括圖例、比例尺、圖名、責(zé)任表以及編圖說(shuō)明等。③圖廓外配置圖名、圖例、比例尺、技術(shù)說(shuō)明、責(zé)任表等，圖名放在圖廓正上方，比例尺放在圖廓正下方，文字說(shuō)明和技術(shù)說(shuō)明放在圖廓的左下方，編繪單位、負(fù)責(zé)人、時(shí)間放在圖廓的右下方。用專門的軟件進(jìn)行異常信息的入庫(kù)工作。35 /