监督分类和面向对象分类流程

1、高分一号城市绿地现状调查与分析实现教程本文将介绍基于高分一号影像数据的城市绿地信息提取的实现步骤，下图是主要的操作流程（图一）图一首先对高分影像进行预处理，其次使用监督分类法和面向对象分类法对城市绿地进行分类，然后对分类出来的影像进行矢量化处理，最后另其在arcGIS中进行统计分析，得出武汉市城市绿地的现状，下面是具体步骤。第一章数据预处理因为处理数据是高分一号影像，本文处理软件为ENVI5.1,因为ENVI5.2以下版本不能对高分一号直接进行处理，所以需要安装r6补丁，将下面两个文件直接粘贴到软件所在位置（图二），然后就可以打开高分影像了（图三）名称KhENV1512017/5/4四文件夹I

2、DL832017/5/4四文件夹悸改曰期ADS40ALCSATSRCRESDADigitalElevattiorDMCDMSP(NOAABinaryCARTOSAT-1EEEditDisplayViewsHelp*11Open.Ctrl+OOpenAskOpenRe-centOpenWorldData1lF_OpenRemotsDataset.RemoteConnectionManageitCtrl+WNewViewsfitLayers图三为了加快数据处理的速度，本文是选择先进行辐射定标然后将图像裁剪在进行后续的操作，预处理流程如下图（图四）：图四辐射校正分为辐射定标和大气校正OBn-Q皿a-

3、Tfth話&rsi?跆on一2问：Q霹F!|曲CD盘栢耳R3(1)打开数据：ENVI-OpenAs-CRESDA-GF-1选择处理的影像，打开XML后缀文件;Op*n*计=机卜耳|EI学曲h生p月匹litr1gdaaflSDSk-l*tII心EXsas-初立廿电兰GFl_FM51_E115_N3aalWnfl_IJJMMJlil3ff!-lMtSl苣GFi_klMll_Eii3ia_HaCl3_JOLia2i*_LLAeCfiULa6&i-Pi5MSIi二GFl_l51Eliaa_IM3DLW214J_LiMMliL35WMSLjlGFi_&i_Eiln_NaLl5JiOL6D21&_LLAO

4、COUL36fiOPArJi二地I圈犯S_2O16021J_L0422IW9FHjlGFi_&i_Eil4_NaL?JOL6D22G_LLAOM12JJDE41MtSL-GF1PMfilE1U6N30530160沏Li山网1也兀汕卩駅辐射定标：选择Toolbox-RadiometricCorrection-RadiometricCalibration选择待处理的高分数据ToolboxSearchthutoQlbvK/RadidiiEtricC*Trecllm/EtadicinetriciZalibration00j*lEliLIi-y-!i：i-ji-i-i*_:门圧iftcisaidsinB

5、TLJEsdarRsi.onetrieCorrection:-4SArijjlyCalnandOffsetjk!_|AtmosphericCorrectionModuleiCalibrateAVBFLR-4sCalibrateTIJ13Cross-TxackIllminationCTrealLon4DarkSubtract!qhEFFORTPolishing；EjiLEsivityALphEcsiduilsEmissivityHornalizatlonEnissivityReferenceChanneLEnpiricalLineCalibratsEsis-tins:-4EnpiTiealLin

6、eCmputmFactorsiFlatFieldCalibration4fi-IASReflectanceCalibrationRumlduaJs：匚3.11107屮1.血RHdLiuti7i.c：CaLibra-tiQnThermalAtiospherlcCorrect1unOlRasterManagBJiEtitQiKeglnsofInterestSFEARSpectralStatisticsTargetDetecticitiTHOR丁亡rrainTxansfornVectorEiimsioHESelectInputFile.卜FileInformationSpatialSabset.Sp

7、ecSubeet.FullExtent4of4BandsQK|Cancel弹出RadiometricCalibration对话框，进行如图设置。对于多光谱影像，点击ApplyFLAASHSetting设置成默认值;如果是对全色影像进行辐射定标，那么Calibration则是Reflectance,OutPutType为UInt,ScaleFactor为1000，如下图：高分一号全色影像参数设置（3）大气校正：选择Toolbox-RadiometricCorrection-AtmosphericCorrectionModule-FLAASHAtmosphericCorrection,弹出FLAA

8、SHAtmosphericCorrectionModelInputParameters对话框。要注意，全色影像不做大气校正，多光谱影像则需要做大气校正处理。Searchthe/Radioriet；ricCorrection/AtnosptiericCorrectionHodL,l-JFavnariteci!_|AnonalyDetectioni口EindRatio_J口ChfltiseDetetiiiJLiClassificatLQn!_|FeatureExtractioniClFilterHOCeoretiiCorreclion-I_JIiastSharpen!nsil_JLiDAE?:!_

9、|FlGSslckLng山口RadarjJRadiQietrlcCTTectLon；:-4HApplyCainandOffset:F-jAtdDEphericCorrectIon11e-djiile|EmSHAhgphEriv：hrrmicinl:OVickfttnc-sphrricCorrEctinn(QTIAC):：-l3iCalibrateAVrtRE?；./,：.ilil:il.-.IY-Crose-TraekIlliULLnatisnCarrestIonLarkSubtractiri:：-4HEFFORTPc-lishinEI:-#EmissivityAlphaRiesi1s翻Eni

10、ssivltyWornallzatlonEnissivityRtfermetChannelEripiTicalLineCalibrateEistlnEEnpiri:：alLineCnputcFaclvrEj-FlatFieldCalItirstl&n：-4ilIARRe：ilestajiceCalibratlutiLogR亡sidualsCalibrationRadioinetrLCCalibrationTherrialATnspherleCdrreeTl&n由RasterHansE&iLgrtAORegionsInteiest申CJSFEAE?1I_jSpectraLil_|Statisti

11、cskOTaisetDete匚tonRsatnaneforTLAftGEiFilesEyerieCanter|MMELat00Cl00Lon00Cl00SensorTypeJliKNOWW-HSGerixorftltitudeOsn)山WOGroLOjlfilevaTi-sn(kii0-凶0FiirlSizetiJ0-MOFllshtDate匚”二=FlightTiie37O:I:SS)o*；a_i_Jilxcsphrri匚Tode-1TropicalAefcqIHode-lturalEticctrilPlishins日111呛terRetrievalYcjKAefczo1Retrieval2

12、-Eand(XTFidth(ihErcfhaDife)9TaxerAbsirpTionFeaxure1133mInit5ilLsihLlity(kn)40.TavcleiifttfciR&callbratlon.防ITAd=e-dSettLng-.Save.Jestore.Hyptrsp&zLralSettings.点击InputRadianceImage,选择前面处理好的数据，在RadianceScaleFactors面板中选择Usesinglescalefactorforallbands,由于定标的辐射量数据与FLAASH的辐射亮度的单位相差10倍，所以在此Singlescalefacto

13、r选择：1,单击OK；5）设置文件输出路径。6）传感器基本信息设置：I-Bi-IT2014.I.：r.-.n-：.l,：l.:EmeuCenter|H)昭口血gf-111：T矽V卜-I：H-th-n-fi.OdO-：1：l932.-:沏.:-lb-ie-O)SceneCenterLocation从影像中自动获取；SensorType为GF-1；GroundElevation通过统计DEM数据获得；PixelSize根据相机选择，PMS相机全色2m,多光谱8m,WFV相机16m；FlightDate从影像xml头文件中读取，减去8换算成GMT时间；anii-uc.i_r*urai_il.bjla

14、uiimiZerterTimo2014-04-1411:11:567）大气模型和气溶胶模型：AtmosphericModel根据经纬度和影像区域选择：IjtrtLdeN)Jan.h.tarchh.(3yJlIi/SeptNdv.GDswG/WGAA/MLWMLW70SAWS/WMLWMLWMLWSAJlV60HLWMLWMLWSASSASMLW50ULWMLWSASSASSASSAS40SASSASSASMLSMLSSAS30MLSMLSMLSTTMLSZDTTTTTT10TTTTTTDrTTTTT-1DTTTTTT-2DTTTMLSMLST-3DMLSMLSMLSh-1L&MLS-4iiSA

15、S5ASSASSASSAS5A5-50SAS泗SAShdLWMLWSAS表2数据经纬度与获取时问决是选用的大气模型AerosolModel根据实际情况选择；AerosolRetrieval选择None；其他默认；(8)单击MultispectralSetting按钮，在FilterFunctionFile导入光谱响应曲线“gf-lpms.sli”单击OK。Indeslafirstband1越MultispectralSettingsSelectCh：srmelDefinitionsbyFile9GULWaterRetrievalKautJian-TaixreAerosolRetrievaLWa

16、terAbsorptionChannel|UndefinedWaterReferEn匚亡ChannelOKCancelHelpFiLterFunctionFiletJHL.astcrM冋化数推液普库-吕丄i-iUndrfined(9)单击AdvancedSettings,在高级设置中ModtranResolution选择5cm-1。(10)设置好后，在大气校正模块面板中，单击Apply。影像裁剪File-Open,这里使用的是不规则裁剪，具体的理论知识可以在网上搜索，首先打开你矢量的矢量.shp数据，如图；在Toolbox中，打开RegionsofInterest/SubsetDatafrom

17、ROIsSelectInputFile选择Beijing_TM.dat,点击OK,打开SubsetDatafromROIsParameters面板；在SubsetDatafromROIsParameters面板中，设置以下参数：SelectInputROIs:选择EVF:矢量.shpMaskpixelsoutputofROI?：YesMaskBackgroundValue背景值：0选择输出路径和文件名，单击0K执行图像裁剪。正射校正使用RPC文件进行正射校正，点开GeometricCorrectionOrthorectificationRPCOrthorectificationWorkflow

18、，如下图所示:ASCIIBuildDuildBuildFfCOrthorsctiFiciqtiuilDrthoi-ec11MapL-uijrdirL：ateLunverty：Orthorectification砖恥ildRFG在该工作页面中将要处理的影像输入，然后点击next在弹出的页面的Advanced选项中选择输出像素大小，多光谱选择8m,全色选择2m，然后再export页面选择输出的位置，点击finish；RPCOrthorectificaticnRPCRefinementImportGCPsAdvancedStszExpertVIC-fioid.CarractionOutpui.Pix

19、elSmtersImageResairiplirigBilinearGridSpacing10FrevievCBack|Finish|Cancel1.4图像融合对每一景的多光谱和全色影像进行正射校正之后就对其进行融合，一般来说高分一号影像使用GS融合，点击工具箱中的ImageSharpeningGram-SchmidtPanSharpening分别将多光谱和全色影像输入，然后设置输出位置，finish；._J:S-.4-oe:ELt-chtheto口丄ijKjImageSh：irpemrLgii1Gr：=jtiSclurilAtP：tl土.1.-_J.-_J-_J-_J-_J-_J.;-s-:

20、s-:s-:s-0-;.0-:s-0-FavQrilesEandRatioCh：rLeEetectio口CliSEificaticinFeatirft!EktractionFilterGecrriHtricCorFeetic-n：=LgeEh：=Lt_iiernngC1I5pactr：=ilSh：=Lrjiernn耳iColorfform：lize1l.)Brovey)Gr：=jtiSchniid_tP：itlSh：=LrpaniHSVSh：iXF已口Jtl巨FCSg已ctral5h：u-peniiLgL1IIWLrlritirFlieSeiectionGF2PMS1E92.&H28.2201

21、61113L1A0001959S01-HSS1ortLo.JatGF2_FMS1_E日彳2_庇也3_2O15L2O9_L1ADOIGF2PMS1E92.&20161113LIADO(1.5图像镶嵌图像镶嵌也是所谓的羽化操作，一般来说一景高分影像能覆盖的区域是有限的，所以需要进行镶嵌操作，对其进行拼接到所需要的研究范围，点击工具箱中的MosaickingseamlessMosaicLfSLlF旨:孟七低uExt匸wetiqilE肝/=%,B-|PiLterE-ClGecimetricCorrectionJ-(JIfTiag已ShiiLfpenirLt由國LiDAHS-tr3Mosaicking卜

22、毡FiselEasedfilosaiekingTiyitllbiE|W|E-f|E：刃:辽|JiadiCiffie+rieC心rreetioilAr|EasterM:dTLagemertt-T|-4*|JiegirmsoInterestE-;_JSPEAJ一4-1Sp已(：tral镶嵌过程最好两张两张镶嵌，这样电脑可以承受，先将要镶嵌的两张图通过绿色的加号加载进来，然后使用右上角的seamlines中的自动接边线进行自动拼接，再观察自动拼接的效果，如果有不好的地方再对其进行修改。接边线的修改主要就是画三角形，将不想要的部分去除，这个很简单，在网上也可以搜到。如何觉得对接边线进行修改比较麻烦，可

23、以加我qq545698595接边效果:1.6影像增强在遥感集市网站上有植被增强的插件，如果要做绿地信息提取的话可以去上面下载安装然后对影像进行处理，但是很多情况下不需要这一步，对影像做一个拉伸处理便可。第二章绿地信息提取2.1监分类这一步实习课上都有，不需要演示2.2面向对象绿地信息提取面向对象信息提取是目前最流行的信息提取方法，现在可以实现的软件有econition8.7ERDAS还有envi的EX模块，总体来说易康软件的面向对象提取方法比较完善，分割方法也比较多，但是因为软件需要收费，所以本文使用envi来进行处理。2.2.1分类步骤是1发黑愎分割合并分决特征提取I定矣要秦规则分抡管分类耳

24、出要盍查看报告和缄+2.2.2实现对象的创建影像对象的创建主要进行了图像分割处理，本文采用多尺度分割算法对影像进行分割。首先要对数据源进行一些预处理操作，本文将RGB调整成4:3:2显示模式，这样能够更好的判别出绿地和非绿地的区别，这个时候在图像上可以看到有些红色的建筑物顶和植被覆盖区域发生混淆，这个是之后需要解决的问题。然后打开基于规则面向对象分类模块，将影像输入，考虑到房屋和绿地之间的颜色混淆，在定制波段一栏选择归一化差值，将红波段和近红波段输入，可以在后期的规则特征提取中提供额外的属性数据。然后对图像进行影像分割和合并，经过反复试验，本文采用使用Edge和FullLambdaSchedu

25、e方法，阈值分别为39.9和89.7，纹理内核为3。分割效果如下：图切割效果预览基于规则特征提取本文需要提取绿地，分类主要依靠的是绿地的植被覆盖指数。正常情况下，当NDVI0便可以认为有植被存在，本文发现当令NDVI在0-1之间一些建筑物的房顶也被提了出来，所以将最低值提高到了0.1，虽然发现还是会有一些建筑物被提取，但是情况已经好了很多，其次因为建筑物大多是矩形的，所以可以通过设置矩形度对房屋进行剔除，因此对绿地的描述设置为：NDVI在0.1-1.0之间rectangle在0.0182-0.8之间，两个属性的分类结果如下图（同一区域），光谱值和纹理提取出来的效果良好，将结果进行输出：图两种属

26、性提取信息预览（左基于光谱，右基于几何）试验区绿地提取从分类结果可以看出分类的效果良好，将肉眼不能分辨的绿地提取了出来，适用于高分辨遥感影像的绿地信息提取，接着将其应用于全图：图分类结果图2.2.3精度评价精度评价有两种方法，在classic中的classification-postclassicification-confusionMatrix中的两种，一种是从高分辨率提取的图像验证低分辨率的图像，另一种是勾选正确的样本来验证用面向对象方法提取的影像信息，我使用了8m分辨率的多光谱影像提取结果和融合后影像的提取结果进行对比，来得出此次面向对象分类的准确度。,1望船TlIffi我Ex旳fCld

27、-IGChangeDectic-nttackksGene-rarte-Randamsampledllimp匚lfl-；E0PastCJais-srfi匚MionA?k|n匚ItsColorsRdClaccilTifiirOh5StabbcsCanfutiori冏之山H.0Curws-CambrieCh-sses巧daicacBjftcirZcEm-smrS?gn-i?nt3rlionInwgu百ifkelim2Vectoi1|Clar-sificathnITpnTcirmFit?rSpiralMapVertexTopographirRadarWindowH*口isicTookSupcrrsed

28、Lhjpc-nii-sedDed-sicnTreeEn-nernberCcJlHcbonCrt-araClasmSm*g-from!RO9si-匚兀-.-.也怒心-UcingGrcijndTuihliiaga*UnjngGrwndTnjthRDh第三章矢量化处理3.1矢量切割在arcgis中打开.dat文件，然后把提取出来的绿地信息进行分级，我的论文中将绿地分成了：公园绿地、生产绿地、防护绿地、附属绿地和其他绿地5种绿地类型，首先将在ArcGIS中新建一个和.dat文件一样轮廓的面文件命名为五类，如图:二耐會iiELJESXULCCLUEMUKKta?旳加曲丁出“击上丄Rd淖知严Ufctaii

29、nHi3f逊MUmnCdiTM-EJir:IXXKDD然后打开编辑器选择切割按钮将影像切割成几个大的部分，再对每个部分进行细部分割，通过与原影像叠加的目视解译和实地外调确定绿地类型。以上为五类绿地分割的情况，新建CLASS字段来对每一块区域进行赋值。3.2矢量处理将处理的.dat影像转换成shp文件：点开小红箱子中的数据管理工具datamanagementtools-features-featuretopolygon将面向对象提取的影像转成矢量文件日閤j|DataManagementToolsS屉AttachmentsE屉Data匚omparison裁DistributedGeodatabHe

30、1I1辱DomainsS屉FeatureClassB裁Features气AddXYCoordinates气Adjust3DZBearingDistanceToLine气CheckGeometry气CopyFeatures甄*DeleteFeatures气Dice込、*FeatureEnvelopeToPolygonFeatureTo-LineFtaiiureToPointFeatureToPolygon%、FeatureVerticesToPoints艺.：、MinimumBoundingGe-ometry气MultipartToSinglepart如下图，是转换出来的全部面文件，需要打开编辑器，将最外面的面文件删除，剩下才是我们真正需要的面，删除后的矢量文件在下图：点击ArcmaplO.1界面的第一行selection的selectbyattributes,打开页面之后选择CLASS字段，以生产绿地为例，其类别为1，以下高亮区域便是生产绿地的范围：接着点击selection中的selectbylocation,目标