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