(完整版)遥感地质蚀变异常信息提取实验

1、基础地质学”创新性实验遥感地质蚀变异常信息提取实验基础地质学”国家级教学实验示范中心二O二年二月二十二日49目录1 实验目的32 实验内容32.1 熟悉遥感影像的辐射定标的方法与流程32.2 掌握遥感影像的波段合成、投影转换、影像裁剪的方法32.3 掌握ETM+遥感影像的Flaash大气校正、掩膜的应用方法32.4 掌握ETM+遥感影像羟基和铁染异常信息提取的方法与流程32.5 掌握ENVI与Surfer软件协同制图的方法。33 实验要求44 实验条件44.1 软件平台：ENVI4.6、Surfer944.2 遥感数据源：金川地区Landsat7ETM+遥感影像45 实验原理45.1 蚀变异常

2、提取的地质依据45.2 蚀变异常提取的物理依据45.2.1 利用主成分分析方法提取矿化蚀变信息55.2.2 铁染蚀变异常分析55.2.3含羟基类矿物和含CO32-矿物蚀变异常分析66 实验步骤76.1 金川地区ETM+遥感影像辐射定标76.2 金川地区ETM+遥感影像不同波段的合成96.3 定义金川地区ETM+遥感影像的地理坐标116.4 对金川地区ETM+遥感影像的地理坐标进行投影转换136.5 对金川地区ETM+遥感影像进行裁剪156.6 对金川地区ETM+遥感影像进行FLAASH大气校正176.6.1 数据转换176.6.2 编辑头文件信息186.6.3 进行FLAASH大气校正206.

3、7 简易去除ETM+遥感影像的干扰信息236.7.1 建立ROI（感兴趣区）236.7.2 建立掩膜256.7.3 应用掩膜266.7.4 掩膜的反选286.7.5 掩膜反选后的应用296.8 主成分分析306.9 提取蚀变异常信息336.10 用Surfer软件修饰铁染蚀变异常信息356.11 总结417 课后练习及作业421 实验目的通过学习基于地学软件的矿化蚀变异常信息提取的流程，了解利用主成分分析方法提取矿化蚀变异常信息提取的原理，掌握ENVI4.6以及Surfer9软件的基本使用方法。2 实验内容2.1 熟悉遥感影像的辐射定标的方法与流程通过学习遥感影像的辐射定标的方法与流程，掌握用

4、ENVI软件进行辐射定标的具体操作加深对遥感图像处理知识的理解。2.2 掌握遥感影像的波段合成、投影转换、影像裁剪的方法通过学习遥感影像的波段合成、投影转换、影像裁剪的方法与流程，掌握用ENVI软件进行这些操作的具体步骤，加深对遥感图像处理知识的理解。2.3 掌握ETM+遥感影像的Flaash大气校正、掩膜的应用方法通过学习ETM+遥感影像的Flaash大气校正，掩膜的建立与应用方法，熟练掌握ENVI软件进行这些操作的具体步骤，加深对大气校正知识的理解。2.4 掌握ETM+遥感影像羟基和铁染异常信息提取的方法与流程通过学习利用主成分分析方法对ETM+遥感影像进行羟基和铁染异常信息提取的方法与流

5、程，熟练掌握ENVI软件进行这些操作的具体步骤，加深对矿化蚀变信息提取知识的理解。2.5 掌握ENVI与Surfer软件协同制图的方法。通过学习使用ENVI软件提取ETM+遥感影像矿化蚀变信息与Surfer软件绘制矿化蚀变异常空间分布图，综合掌握ENVI与Surfer软件协同制图的方法，加深对所学知识的理解。3 实验要求上课时，需认真听取实验老师的细致讲解；练习时，按照实验指导书的步骤，自己亲自动手使用ENVI4.6软件提取矿化蚀变异常信息；以及在Surfer软件中制作矿化蚀变异常信息空间分布图。4 实验条件4.1 软件平台：ENVI4.6、Surfer9ENVI4.6软件为遥感图像处理软件，

6、用于提取矿化蚀变异常信息；Surfer9软件为一款以画三维图（等高线,imagemap,3dsurface）的软件，具有的强大插值功能和绘制图件能力，使它成为用来处理XYZ数据的首选软件，是地质工作者必备的专业成图软件。4.2 遥感数据源：金川地区Landsat7ETM+遥感影像Landsat7于1999年4月15日发射升空后，由于其优越的数据质量，以及与以前的Landsat系列卫星保持了在数据上的延续性，现在已成为中国遥感卫星地面站的主要产品之一，2003.5月出现故障退役。5 实验原理5.1 蚀变异常提取的地质依据目前认为中等强度以上的蚀变带对于TM（ETM+）蚀变信息提取是十分有利的。有

7、时尽管有蚀变岩存在，但不一定有矿，然而围岩的剧烈而较大范围的蚀变常常与大矿及富矿石的生成互为隶属，大型特大型内生热液矿床一般均有强烈且较大范围的围岩蚀变，并且具有分带现象（如斑岩铜矿）。这便是以找矿（首先是大矿、富矿）为最终目的蚀变遥感异常提取的地质依据。5.2 蚀变异常提取的物理依据从HuntG.RC（1978）和他领导的实验室的研究成果以及阎积惠等（1995）依据矿物反射波谱特征吸收谱带特点的定性分类研究中可以知道：主要造岩矿物在可见光近红外光谱（0.35-0.25um）并不产生具有鉴定意义的反射谱带，其光谱特征主要由岩石中为数不多的次要矿物决定：是含铁（Fe2+，Fe3+）基团产生，含铁

8、矿物主要有角闪石、赤铁矿、褐铁矿、针铁矿、磁铁矿、黄钾铁矶等，他们在TM1、TM4波段有强的吸收带，若岩石中含多量的Fe3+，而含Fe2+较少,这类岩石的主要吸收谱带位于TM4和TM1波段，反射波长相当于TM3波段的电磁波。若含大量的Fe2+、含Fe3+较少，则主要吸眩Q*in3：Z3-W5*T¥诫配f心=-_畫记收谱带位于1波段，对于波长相当于TM2波段的电磁波有某种程度的反射。二是含羟基（0H-）、水（H2O）或碳酸根（CO32-）基团产生，羟基的吸收谱带主要有二处：2.2um,2.3um，由于OH-在2.22.3um附近存在强吸收谷（称为羟基谱带），使得TM7产生低值，TM5产

9、生高值，含羟基矿物大多为次级蚀变矿物，如高岭土、叶腊石、云母类矿物、绿泥石、绿帘石等，水在1.4um和1.9um处有特征吸收带。含碳酸根矿物主要有五个特征吸收谱带（1.9-2.55um），较强的两个在2.35um和2.55um波长处（称为碳酸根谱带），相对较弱的在1.9um，2.0um，2.16um三处。5.2.1 利用主成分分析方法提取矿化蚀变信息遥感图像各波段间存有一定的相关性，为了减少相关性对分类的影响，常使用主分量分析法去相关。主分量分析基于变量之间的相互关系，在信息总量守恒的前提下，利用线性变换的方法来实现去相关性。由于所获各主分量之间不相关，故各主分量之间信息没有重复或冗余。故主分

10、量分析这一基本性质在蚀变异常信息提取中被充分利用并很有成效。5.2.2 铁染蚀变异常分析针对各个波段波谱所反映的特征，选取TM1、TM3、TM4、TM5这4个波段作为组合波段做主成分分析。这是因为铁氧化物的特征光谱信息集中在了TM14波段，在TM4和TM1波段有吸收峰，在TM3波段无特征吸收而呈高反射。同时为了避免含羟基和碳酸根矿物的干扰，在选取波段组合时舍弃了TM7波段。对TM1、TM3、TM4、TM5应用掩膜做主成分分析，统计分析如下：特征向量(Eigenvectors)主成分Band1Band3Band4Band5PC10.4153040.5278860.4966230.549750PC

11、2-0.204294-0.5514550.807515-0.045623PC3-0.636068-0.097958-0.1858670.742479PC4-0.6174170.6384660.258339-0.380024从上表可以看出，PC1主要反映了TM3和TM5波段的信息；PC2主要反映了TM4波段的信息；PC3反映TM5波段的加信息和TM1波段的减信息；PC4反映了TM1波段的减信息和TM3波段的加信息。根据铁染类蚀变矿物的波谱特征，包含这类蚀变信息的图像应该具有TM3与TM1或TM3与TM4具有相反的贡献值且绝对值较大，从而选择PC2分量，故PC4分量中暗色调部分表征了铁染信息。5.

12、2.3含羟基类矿物和含CO32-矿物蚀变异常分析选取TM1、TM4、TM5、TM7这4个波段作为组合波段做主成分分析【4】。这是因为粘土类矿物(含羟基矿物)和含CO32-矿物的特征光谱信息集中在TM5和TM7波段，在TM7波段为特征吸收带，在TM5相对高反射。主成分变换的波段组合的选择也是以此为依据的，由于可见光波段对铁氧化物敏感，为了避免铁氧化物信息的干扰，故只选择了一个可见光波段参与运算。对TM1、TM4、TM5、TM7应用掩膜做主成分分析，统计分析如下：特征向量(Eigenvectors)主成分Band1Band4Band5Band7PC10.4250780.5101330.565080

13、0.489650PC2-0.1433960.809782-0.146893-0.549649PC3-0.8701090.0487390.4575770.176519PC40.204102-0.2857310.670619-0.653429从上表可以看出，PC1主要反映的是TM4和TM5波段的信息；PC2主要反映的是TM4波段的信息；PC3反映的是TM1波段的信息；PC4反映的是TM5和TM7波段的信息，且符号相反。根据含羟基类矿物和含CO32-矿物的波谱特征，PC4图像中的高亮度信息就是表征含羟基类矿物和含CO32-矿物的蚀变异常信息。6 实验步骤6.1金川地区ETM+遥感影像辐射定标在ENV

14、I4.6软件中，单击菜单栏中的fileopenimagefile工具，打开金川地区Landsat7ETM+遥感影像中bandi影像文件。然后，单击菜单栏中的BasicToolspreprocessingCalibrationUtilitiesLandsatTM工具，为金川地区ETM+遥感影像中bandi影像文件进行辐射定标，如图6-1所示。Ifile喙EokC7AvailableBandiLktFileOpTionsSalctedEsrdBand：1ipL3Zt：i：i：33_7t2KtZtiTS4_xs|S_nj-iL0.tLISUTlstlcsChangeDetectionMe*5urer

15、nentT口IGene-rjlPurposeUtihtiie&DandlMaithSpectralMarthH-egiomOfMowkkingAVH1RRLardmtLand/tTMFleslzeDatdSpailal/Spctral)SubsetDataviaROIsRjotate/Flipl>staLayer£tadkingCcnvi-rtDaria(KQrBJL,PIP、ScrttchDat白Gray滋虹亡ChobDolor曰L3pl3ar033_Tt3002W2OiS_ixd0lllfl口1H囂t-ipId£o臼內址Tool盯ClB5tTfk-ai：i

16、c<!fransforTiFilterSpectralMapVectorTopographicRadarWindowHe*pr-=praceg春=一_Ir1.4FL._II,I图6-1打开辐射定标工具对话框在TMCalibrationInputFile弹窗下的SelectInputFile栏中选择Landsat7的bandi文件，然后单击OK，如图6-2所示。TMOlilxationInputFile2011bpatlalSubs亡tSenserTyre;UnkrrnByteCig;Elvst(Intel；)Ei-ojecrticn:0TN?ZFixel:28,5ffEtin:KS-&a

17、mp;L阳velength;ffcce图6-2TMCalibrationInputFile对话框接来下，在TMCalibrationParameters弹窗中，单击GetCalibrationParametersfromWeb选I1:i:I1'11SunE1d¥?ition(.des)"TMI匚alibrationParametersDataAcquisitionffonthDataAcquisitionDayDataAcquisitlonYearScale:Kin:0-COGetCalibrationParanetersfrn.VebQueuieCaneElDat

18、aAcquisitionBand:存EnleiOcilputFilenan亡ChooseCompressCallbrall<inTypeCRadiance(*ReflectanceOutputResulttofl-FileCNenoryLandsa!SatellUe厂4C5i-ETH+?项，如图6-3所示。图6-3TMCalibrationParameters对话框在TMWebCalibrationParameters弹窗中，选择Date/Path/Row/Band选项，根据金川地区Landsat7ETM+遥感影像的头文件信息，依次填写下面红色标注处的信息，如图6-4所示。注意：每一景影

19、像的信息都不一样，一定要根据头文件信息填写。图6-4TMWebCalibrationParameters对话框填写完成后，在TMCalibrationParameters中依次填写下图红色标注部分的信息，如图6-5所示。图6-5TMCalibrationParameters对话框单击0K,即完成了金川地区Landsat7ETM+遥感影像中band1影像文件的辐射定标操作;接下来，依次对band2、band3、band4、band5、band7波段进行辐射定标操作。6.2金川地区ETM+遥感影像不同波段的合成首先，在ENVI4.6软件中，单击菜单栏中的fileopenimagefile工具，打开

20、经过辐射定标的金川地区ETM+遥感影像的band1、band2、band3、band4、band5、band7数据文件，如图6-6所示。图6-6AvailableBandList对话框然后，单击菜单栏中的BasicToolsLayerStacking工具；如图6-7所示。|Bi鬲kTgkClassificationTransfcrmFResizeData(Spatial/Spectral)SubsetDataviaROIsRotate/FlipDataLayrStackingConvertData(B£Q,BIL,BIP)StretchDataStatistics.#SpatialS

21、tatisticsk匚hangeDetectionMeasurementToolBandMathSpectralMathSegmentationImageRegionOfInterest卜MosaickingMaskingPrepracessing卜图6-7打开波段合成工具在LayerStackingParameters弹窗中，单击ImportFile按钮，将经过辐射定标的金川地区ETM+遥感影像的bandl、band2、band3、band4、band5、band7数据文件全部添加进来。之后，单击ReorderFiles按钮，在ReorderFiles弹窗中将顺序从上到下调整为：band1、

22、band2、band3、band4、band5、band7。选择好保存的路径，单击OK即完成了金川地区ETM+遥感影像不同波段的波段合成操作。400KQKLayerStackingParametersSElRUt聞PilETMrLay已tStxKlng:Outputtfabandl-ra.dBandljband2-riidBand.1batjd3-radBand1band4-radBand1Iip<?rtFile,R4;OutputFileRang亡：bandS-radEnii1Units.Inclusive:IranE1?allthefiles广Exclusive:rHEetn.c&l

23、t;apa55esfileoverlapCarK-rlOutputResultW檻File厂ffemvryVPixelMeters|CsDaeraAdii1nistratorADealitdpJIIfttlZETlCancel11111111EnterOutputFilensueChoosebandl-radbarulZ-radbaiid3-xadband屯一TadbandS-iad甘iindT-TMdBand1Earid_BandBandBandfgiind£:band7-ra.dBand1Size2B.SOOOOOO?ReorderFilesArbitraryGeofiraphi

24、<StatePlatStatePlatArgenlinaArgcnlinaArseniina_Resimpling|Hrarcs：tRcighbpr二ReorderFilesDelete图6-8LayerStackingParameters对话框6.3定义金川地区ETM+遥感影像的地理坐标首先，在ENVI4.6软件中打开金川地区ETM+影像（波段合成后），并加载彩色合成影像;在主菜单中单击Map中选择CustomizeMapProjections选项。h吨JVectorTopographicRachrWindowRegistraticinRigercxKOrthorectilFicati

25、onOnhorectifncatloob-Momldklngi-GeorefererxeftcniinputGeornetryNeorefereinceSPOTh-Georeference5eaWiF5p-GeoreferebceA5TE.R-GcDrefercnceAVHRRp-GwirftrcircuENV1SATCieoreferenceMODISeoreferenceHADAR5ATBuildRPCiCiKtomizeMspProjecticwisConrtMap叶屋ktiooLayerStatklmgiMpCoorcfinatcrConverterASCHCoordmateConve

26、rsionMergeOld"map-projuixCFileGPS-Link图6-9打开定义地理坐标工具然后，在新弹窗的ProjectionName中输入投影的名字，如jinchuantouying；在Falseeasting中输入500000.00；Falsenorthing中输入0.00；在ProjectionType中选择TransverseMercator；在ProjectionDatum中选择xi'an-80；在Latitude中分别输入0，0,0.00；在Longitude中输入102，0，0.00。Scalefactor中输入1。订：WTIFL1Pt-DJHCt

27、ionDpt.iDE1.X5OTOOT.mt'rdJctionITarnaFi"c-j*ft.ian.Typaitud*Epr«cti«z-nv-ri.不in.I3I3I31-aztia£am*vid.3*ztLrdXL丑屯7汪0IaScalefactorU1distok60LtlIIEllip5oia|lAC-7G鼻母|D_jmudxi.Z'n-urnle£4tVWWWYZBKimF口即加图6-10定义地理坐标对话框接来下，输入完以上数据后，在菜单栏中点击Projection,选择AddNewProjection,载这个投影信

28、息。rrojectuiLHatonJll:illjOLtiLevn1916IYalxe-En.iyetok、GOWGS-72恼-54facsreZajidsrijEEi_jing54GrWagner3YWagnerVTIEllipsold175Froj皀匚tLQTlTyp电SpaceObliineHercatorASpaceObliqum1'ilercatorIEtor4jgt-aphic（sphere：）Stereographicnorthng.CODDOO图6-11定义地理坐标对话框最后，点击File中SaveProjections,保存这个投影信息。CiisToSieaTapPr

29、ojectinTiDefiEWToriyaltPrajfic七。!！Optnar>zFrojectionTsreSpaceObliQueNercatnrAtge-SpaceOliliqaeMercatorESt0&£1*phic1.zphc：r4)Stereo:-:沁hic(ellipse)kIranv«rseMercatorVsrid.srGriatenflagnerIVflaguerVI工<1"1l>tjprojeclicinDalwriitLevu1916Waic-Enifteioh'BOW&E-72WK-34Tac

30、ireZ-utdQrij血L_jing54nnElLiiEoidlyj-7'5图6-12定义地理坐标对话框6.4对金川地区ETM+遥感影像的地理坐标进行投影转换OtrlIlDsaickliLgSTOTsWiFSASTERZJfVISWIflflDISMDARSK7i,&&F5-LinkBtl辽dRTCsLayerStackingASCHCDardin.B.teConversidilCnstomizeMapTrojeclionsCooteConvtrterRagi匚七旷ntionEnl'i-mctTools.WiiLdcwI'iVap¥±

31、;ciarTopogrELphia%indo、Kigorou5Orthorectificati&ilDrthorscti£i&.1ion2.U8-2a3bS'1S1(KsBeEize(UbzidCinvertNepProjcation首先，在主菜单栏中选择ConvertMapProjection,进行投影转换。图6-13打开投影转换工具其次，在新弹窗的SelectInputFile里选择要进行转换的影像，如甘肃金川jubu.img；然后,单击OK。图6-14投影转换输出影像窗口接着，在新弹窗中Resampling里选择Bilinear方法；同时，选择好要保存的

32、路径。然后，单击ChangeProj。OutputProjection,and.E玄toutUpparLeftCorriarLo»rd.3n-at«Proj:jLnckarigtOTiyLngD县tun：II'*as.-EO1|Vidts:netersKiLSTfllSiiapT0403396.BOODCOuDnBt-cre启SaOCEdODOgiersOutputXSix=pixelsOatpntEe=ultio<*Flit矗niDtryOutputYSis*1393图6-15投影转换参数设置对话框再在新弹窗的SelectNewProjection中选择前面

33、保存好的投影jinchuantouying；然后，单击OK。New.Jtfbitr：±i_yGei：iaphicLat/LonUTMStateFl：me(W27)EtalePl=me(WE：3：liinchaiLEtoiLi123yang：匸IdengHatiaiTi.图6-16投影选择对话框下图为投影后的彩色影像。图6-17投影转换后的影像6.5对金川地区ETM+遥感影像进行裁剪首先，在ENVI4.6软件中打开波段合成后的金川地区ETM+遥感影像；在菜单栏中单击BasicToolsResizeData(Spatial/Spectral)工具。BasicToolsClassifica

34、tionTansfbrrrFResizeData(Spatisl/Spectrl)SubsetDataviaROIsRotate/HipDataLayerStackingConvertData(BSQ,BILBIP)StretchSLdtistic5SpatialStdtisticE*ChangeDetection*MeasurenientTodBandMathSpectralMathSegmentationImageRegionOfInterest卜Mcaickin*Masking»Preprocessing卜图6-18打开裁剪工具在ResizeDataInputFile弹窗中选择

35、要进行裁剪的影像，单击SpatialSubset按钮；ResizeDatwIi叩ut卜eI'li-.11:'iri.:i“ii：Fii-：Dins:0922x8045/6ESQSizt:Fl：'itins?:i：71.722.3:?.':vteFileType:EtlVIStandardSensc-rType:UnkncnvnByteOrder:HostQmel)ProJedion.:UTII,Zone4HNorthPlsei:28,!5HetergDatiia:VC-S-Q4VavElengtli:IToneUpperLefiComer:Lj1Descrlpl

36、lon:CreateLayerFileR亡suitTurJun20LG:S3:i52011图6-19输入裁剪数据对话框在SelectSpatialSubset弹窗中，可以通过Image、Map、File、ROI/EVF四种方法进行裁剪，根据需要裁剪出相应的研究区域。最后，单击OK,即完成裁剪操作。SelectSpatidfSub-setIBPile:金川地区辐射定标Dims:B922x0045(FloatingPoint)FullSize:287j10%960bytesSubsetSize:287r109j960bytesSubsetUsingI皿ag已|JU口|Fl巳|ROI/EVF丸&qu

37、ot;11|SubsetbyImageInputBand卩ResetPrevio-usOpenOKCancel图6-20选择空间裁剪区域对话框6.6对金川地区ETM+遥感影像进行FLAASH大气校正6.6.1数据转换首先，在ENVI软件中打开裁剪后的金川地区ETM+遥感影像；再单击BasicToolsConvertData(BSQ,BIL，BIP)工具，对需要进行FLAASH大气校正的数据进行数据转换操作;日融引fTook|Classiflcadon1ransformRResizeData(Spatial/SpectraDSubsetDataviaRO阪Rotate/FlipDataLayer

38、由匚kingConvertData怛辺ETL.ETP)StretchDataStatisticsSpatialStatisticsChangeDetectionMeasuremeniToolBandMathSpectralMathSegmentationIrngeRegionOfInterestMosaickingMasking图6-21打开数据转换工具在ConvertFileInputFile弹窗中，选择要进行数据转换的文件，再单击OK；ConcertFite1叩血FileFileInformation:faet*r!l.iXiXOO.1T;15;53ZCd.1邱处匸F£ll|&a

39、mp;/6EiJjdSe八州州丈件罐阳定硕(反射1!、金IIDins:1531X1911x:6LBSQSize:tFlatingP«imES,384tytca.FilsType;EtIVI飢址加PSensorTyre:llrkrCTnByteCiiJtr;Elvst(Inttl)Eroji-crtico:OT!I?Zco=4SNorthFiirel;38,5ffeleraCatua:VCE-Sl/avelength;0,485WZ,ZSpatialSubsetPiu1ScKnt图6-22输入数据转换文件窗口再在ConvertFileParameters弹窗中选择BIL格式，选择好存储

40、路径后，单击OK,即完成了数据格式的转换。ConvertFileParametersInputIrErlmmBSC|3_i7putInturlctvc:<BILEIP!iniFFrt.r.Pl'丿Eijli：rCuLullIFilci.a.iiT：Choose|C:UsersAdmInLs1ratorDesk1op金川地区ET0EQueueCancel图6-23转换参数设置窗口6.6.2 编辑头文件信息然后,在AvailableBandsList窗口里选择刚才转换格式的文件,右击该文件,点击EditHeader选项。AviableBan.FiJeOptionsGrayScaleC

41、RGBOpenFileinENVIZoomOpenFileinArcMapAddselectedHietoMemoryCloseSelectedFileResize(Layer(T1CalEdi?Hedder.,QuickStats.LoadTrueColorto-<neiv>LoadC1Rto<new>LoadTrue匚olorto<current>LoadORtocurrent*LoadBanddisplay"FoldAlfFilesUnfoldAllFiles图6-24编辑头文件信息在弹出窗口中，EditAttributes，在下拉菜单中选择

42、Wavelengths进行下步编辑。Header山归壬酬州交忡刑镐射圭顾反射）谧川地医子区<梏.，FileSize:SS?5fi3?3BibytesInputHeaderInfo-FromSajiplesI»LinesOffsetA=sstartBandNames.DefaultBandstoLoad.SpectralLibraryNames-FileType|EWVIStandardOKCancelWavelengths,.BadBandsList.DataType(FloatingFoiniFileResizeResultjkreresizefactor:_Q00000.L

43、TucJun2817:16:31201Maplnfo._.RPCorRSMProjectionEmulation.AssociateDEMFilerGeographicCorn已rs.“DlwpI7丁戸盂.FWHWGains.,图6-25编辑头文件中的中心波长信息在弹出的对话框中，依次将band1、band2、band3、band4、band5、band7波段的中心波长设置为:0.4825、0.565、0.66、0.825、1.65、2.22；之后，单击OK。EditWavelengthvaluesRe?etHm'rer.tWa-zelenjlh'Zal'.ies:区韬

44、至疋标:):band!-rad)I0:tiandZ-rad):金HD:tiandd-rad):金丿I'):bandi-rad):金JI:):band5-rad):金JI0.495000.565000.660000.82500i.esonoi7F(.lyei'C|M41：.Kr.r.1：7UK.：-.P.ImportASCII.ClEar图6-26中心波长编辑对话框6.6.3进行FLAASH大气校正最后，再单击菜单栏中BasicToolsPreprocessingCalibrationUtilitiesFLAASH工BasicTochI肉叮庵闭TriformFilterSpect

45、ralM即VectorlopogrspkRwarwncHeipnslzedCSiMtldl/SpecUal)Sub&e-tDaIjviaROEsRotdtf/Fkp禹酋LayrStackingCorwertDatil(BSQFBJL,BIPjStreicihDataSratislkSpatialtatistic-sGh*n学OetKtlcwMeasur&mefitToolBandMaifiSpeechMarhSegmental:kmknngjciRegtonOfinterestMoidkiflrigIvlaskiiig#LScroll:D.l.FLAASHCorrectionC

46、afirdtkjnUtStisOeoerdlPurposeUtlllwsCuta-pec帝匚UtifitiesLad能tMS5LndktTMQuklBirdRdidnce典oddViLHR.fldiwce*Gr卿scaleri>jEcoiwStleiedE-mdffesiifl(Laynr<nCaltBanrSaidov*illUCErEckiSElEe：(LayertTKCal口E«S.xe<La：ys-iI；T1：CaI=口Re3lie(layertilCal口K«lie!：1_吧也fl：nCal 口Ejh1xe>HqyvrCHCal* III&

47、quot;'Di«|TCLxL911(tlaatijiBfA；HHRR具，进行FLAASH大气校正；图6-27打开FLAASH大气校正工具在弹出的窗口中，单击InputRadianceImage按钮，选择前步得到的影像数据后，单击OK再在RadianceScaleFactors弹窗中，按下图参数进行设置，单击OK即可。回申tNVIij.S-LFlAASHAtEm少创itCorn子工AMly駅競歯er化Model|t血归iinxitBlit!*llelnrc-T>atlc-Ti：IjqiwtKaflladicrIaigeOutizytcfletuncicflic：Quiu

48、ux:frlrecKivrsfaBo&tn.weforFLAISHFLlci貼ddrweScdifcFatafsR«ajdaxTiys:fmmliefaci-f-z-s1perbmd!froi.hSCItfileTirMlngF£<«Lefaejserf*r屯11bandjFileBMicf>ofcClti：严“wFLAASHIhpulFileScentCentctLcn-stlon卫Flit：Zi1Ci*>；131kUliz&BIL53he:F«atln(Ftini:SH,bfJU*tortes.FileIyp&

49、;S3TIStaler3SeriFrTyj-e£CnlnwnEytsOrdsr:Clntel)Fr>J«<tlLfrn:TTljZ>ne43IcrthPixel:邮弋Dtua,:初£-£4Hyperspeetr-sLSettings.AdvancedS?r11ngs.-Sa.w.REfft6TC.SLnELes:defa«:tor|l'£时叩|Spat5nl£iih"tWatctRetrievilp>7tt|ffstejAUstbjtiFestweQKCarseOpeti-图6-28

50、FLAASH大气校正数据输入、影像选择及参数设置图之后，在SensorType中选择JNKNOWN-MSMultispectralLandsatTM7工具；FLAA5HAtmosphericCorrectionrdodellnputParametersAVHRR-7In（pu.tRadiaoceInaEc-E：出'I出'I克件ENTI辐克JjilS1（反财）'金川地IK子12格式蒔AVHRK-8Outnui:ReflectiiKDFileC:Fn>granFilesITTIDLT08llbbookAVHRIR-9AVHRIR-10OutputDirectoryfo

51、rFLAASBFiLesCAProgranFile3ITTIDL709libhciQkAVHRR-11SootnanLEfeirFLAASHFlLesAVHRR-12AVHRR-14Scene-CenterLcaticnDD<->DNSSensorTypeNKNOWN-M5FliAVHRR-15AVHRK-16Lat0-ooooooooSensorAlti+iLon|0-DOOOOOOOGT<>undElevaHyperspectralMukispsctralTzrrrr广*、p、zPiselSize(n)28.500AVHRft-171KONOSIRSLISSIirs

52、ussmAtMQsphericMdrl|TrQpicalAerQgglN>?delRuralLandsatTM4ffaterEetrievalSo史JffaterC&luinnul-tiplier|'l-|0吕AcrQsglRetrieval|2_Eajid(K_T)ILandlsatTM5LandsatTM7.-.ii：.TlslhllityCkn)|40.00LandsatMSS4LandsatMSS5ITultispectralSettings.-Ad¥:MODIS图6-29SensorType参数设置图接着,根据金川地区ETM+遥感影像文件的头文件信息依

53、次将下图红色标记部分的信息填完,并设置好保存的路径；然后，单击MultispectralSettings按钮，进行参数设置操作；之后，在单击AdvancedSettings按钮，进行参数设置操作。图6-30FLAASH大气校正工具对话框A.在单击MultispectralSettings按钮，进行参数设置操作时；在新的弹出窗口中，安全按照下面红色标注部分的参数进行设置；图6-31MultispectralSettings对话框按照红色标注信息继续设置其他的选项；之后，单击0K即可。图6-32MultispectralSettings对话框B.在单击AdvaneedSettings按钮，进行参数

54、设置操作时，完全按照下面红色标注的信息进行填写。填写完后，单击0K。在图？中单击应用，即完成了FLAASH大气校正操作。图6-33AdvancedSettings对话框6.7简易去除ETM+遥感影像的干扰信息6.7.1建立ROI（感兴趣区）首先，在IMAGE窗口中点击Overlay，选择RegionofInterest选项。图6-34打开ROI工具然后，在弹出的#1ROITool窗口中的Window里选择在Scroll中建立ROI。左图红色部分为建立的ROI。图6-35建立ROI区域图最后，将建立的ROI保存。Got0FileR0I_TypOptioreHelproffRestoreEOIs.

55、StitsGrowFiKelDele+eRlJl5ubee七DataviaFUJIeUutpu七RUlstoASCII.CidTLCylExportKOIztoEVFExpor七RUIstoEhapefile.ExportROIston-DVi5n：ilizerSelfictMlHiieROIsFolYHl='rLS5/1,530,i0图6-36保存ROI区域6.7.2 建立掩膜首先，在主菜单栏中单击BasicTools，然后单击Masking,选择BuildMask，建立掩膜。F:L*Eai-att/flLpBatIfleasureaientToolT)dtc>clion匚aavflrtBati(BSQ,BIL,ELF)Resile£l.laISpa.tial/Spsclral)t=ic1(?Le|n