监督分类和面向对象分类流程_第1页
监督分类和面向对象分类流程_第2页
监督分类和面向对象分类流程_第3页
监督分类和面向对象分类流程_第4页
监督分类和面向对象分类流程_第5页
已阅读5页,还剩20页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

高分一号城市绿地现状调查与分析实现教程本文将介绍基于高分一号影像数据的城市绿地信息提取的实现步骤,下图是主要的操作流程(图一)外调版证切割及资和蟆(年鉴■匡廉标绿地现状分析与评价矢量化处理影像预处理绿地信息提取总结面向对象分类本文将介绍基于高分一号影像数据的城市绿地信息提取的实现步骤,下图是主要的操作流程(图一)外调版证切割及资和蟆(年鉴■匡廉标绿地现状分析与评价矢量化处理影像预处理绿地信息提取总结面向对象分类图一首先对高分影像进行预处理,其次使用监督分类法和面向对象分类法对城市绿地进行分类,然后对分类出来的影像进行矢量化处理,最后另其在arcGIS中进行统计分析,得出武汉市城市绿地的现状,下面是具体步骤。第一章数据预处理因为处理数据是高分一号影像,本文处理软件为ENVI5.1,因为ENVI5.2以下版本不能对高分一号直接进行处理,所以需要安装r6补丁,将下面两个文件直接粘贴到软件所在位置(图二),然后就可以打开高分影像了(图三)ENV151IDL83CARTOSAT-1CRESDADigitsElevationDMCDMSP(NOAA)EditDisplayViewsHelpOpen...Ctrl+OOpenAsENV151IDL83CARTOSAT-1CRESDADigitsElevationDMCDMSP(NOAA)EditDisplayViewsHelpOpen...Ctrl+OOpenAs►直|L»!OpenRecentOpenWorldDataOpenRemoteDateset...RemoteConnectionManagerCtrl-FWNew►Views&Layers►ADS40ALOSATSR►软件(D:JtenvikenM5,l去装所需►envi51idIB3r&?&H7T►力工具CD群助出)刻录新建文角夹名称,修改日期类型K小2017/5/4星期四…文胖夫2017/5/4星期四…文件夫图三为了加快数据处理的速度,本文是选择先进行辐射定标然后将图像裁剪在进行后续的操作,预处理流程如下图(图四):影像裁剪

~^T~

辐射校正影像裁剪

~^T~

辐射校正正射校正图像融合影像镶嵌影像增强图四1.1辐射校正分为辐射定标和大气校正

3侬m3侬GF1P1MS.1Ei耳口NJ0.3XH独工场LLAmJULJb&lPiiMMlaGFl_™Sl^ElliC!_IMaO5_3aLM214_LLAflCTldL3i5WMSLaGFl_PM51_Elli5_IMaDJ_ZaLWZ2G_LLAfl(ffll42JmF^l0GFl_™&i_Eil4e_N3Ll5JKiL6D22t&_LLAOMi4JJDE41.ML£L"GF1PMS1E1U6N如万30160Z2SLigOOlK加沁3侬m3侬GF1P1MS.1Ei耳口NJ0.3XH独工场LLAmJULJb&lPiiMMlaGFl_™Sl^ElliC!_IMaO5_3aLM214_LLAflCTldL3i5WMSLaGFl_PM51_Elli5_IMaDJ_ZaLWZ2G_LLAfl(ffll42JmF^l0GFl_™&i_Eil4e_N3Ll5JKiL6D22t&_LLAOMi4JJDE41.ML£L"GF1PMS1E1U6N如万30160Z2SLigOOlK加沁二GFi_™&i_Eil<[l_NaDLl5j!CiL6D2i5_LLAflMi4L36fiOPANi2LGFl_OM£lJll^l_«a04_2O16021i_Ll^m)Ul-36S-MSl*GFl_™&l_Ell^l_IM3DLa_2aL6n216_LLA™]UL36&2.FArJl兰GF1_AMS1_£11^2_IM3LL_2OL6O^2a_LLAfH»1^2M2MWL&GFl.m&l.ElKZ-HSLLjnLM3J3_LLAQHll-^2[KSFAHl上GF1FM5LE1地工时3L4_2OL轮迎£LL切0耻后9199&"M珞L三GFl.™Sl_El]A3_H3L4._ZOLM3Ja_LLAfiC(]ML9SSFANl三日FLPM&LEli甩5/^flUL2OLW22O-LLWmi超205SMSSL工GF1加5LE11jraCLC-ZOL如ZZQJJJM1CT13口砂FAN1■IQF1加&1E114.5NJ0132OL和功LUWtii42J0«QM^L1GFl_™51_E115>5_M3a3_Z01Wnfl_LLW:^lil3CnMMlXWLOT却LfllQld至Q1工雕3峪4事五^停沿制jK*加―如阳GafefFI^CDSx»Q闾鼻JbTH.@Op*ni廿区机■KIE0II聿脂4生E吊型Sliktt+HE3AVM卜F”函―呼1JD|融立壮左岩,13④aHtt--刈姻2日嘴丸小标坦主"关■期三1^23侬9K6M通3K63侬3尊3K&3侬mK&3侬9QML朝目LE・期二1至布MLfl/Z/LflIM二13?«2DLWLE!1翼二13:4E加1的。S二】和H2DL加LE百二13c4T州代闻魂耽一M23013/3/23ft*B—]©*&L阴冢期耽一]^103侬9K6M通3K63侬3尊3K&3侬mK&3侬9Q(2)辐射定标:选择Toolbox->RadiometricCorrection->RadiometricCalibration,选择待处理的高分数据ToolboxSearchthe-toQlbvK/Radi5netrieCiTredlc-ri/Radi5netriiCalibratiin回回曜;j阅±口InageSharpening百口LIDARE-OlHc>said£inB土口RadarB-tSRai±.onetrieC^rreciion:--4SApplyCalnandOffset;|+t!_|AlnosphericCoirectionHodule!CalibrateAVBfLR:-4«CalibrateTIJ13:--4mCross-TiackIlliminationCOTrectlon4s-DarkSubtracticn:T®EFFORTPiLiahlng;--4SEnissivity氏LphdEesiduilE4fi-EmissivityNornallsatlonEnissivltyReferenceChanneL42EnpiricalLineCalihrateEsis-ting:--4£'EnpiriealLineComputeFactors:-TgmtFieldCalibration踹I£RReflectanceCalibration:--4mL^BResiduiilsCalibralionRadLicietrlcCaLibxa-tiQn4fi-ThermalAticsphericCorrect1onEl-I^ZIRasterIsnagdb.ent±L3KegLi>neofInterestT口SPEAREQ)Spectralt口StatisticsEl-lZ3TargetDetectionJ-C3THORH--L3Terrain£-Q|TxansfomT口Vector:--SlEiiensionE弹出RadiometricCalibration对话框,进行如图设置。对于多光谱影像,点击ApplyFLAASHSetting设置成默认值;如果是对全色影像进行辐射定标,那么Calibration则是Reflectance,OutPutType为UInt,ScaleFactor为1000,如下图:

RadiometricCaibiratioriCalibrationTypeRadianceOutputImerleaveElOutputDataTypeFloatScaleFactor①IDApplyFLAASHSettings-lOutputFilename:-lF:\JELMa51cAGFT顺处理文档\WFULdatV]DisplayresultOKCancelOKCancel高分一号多光谱影像参数设置RadiomedricCslibrartionCalibratiTypeReflectanceOutputIn+erleaveOutputDataTjcp高分一号多光谱影像参数设置RadiomedricCslibrartionCalibratiTypeReflectanceOutputIn+erleaveOutputDataTjcp虺UIntSealeFactcr1UUU.UUApplyFLAASHSettingEOutputFilerL:=iiriH:HlsplayresuitOKCMLCtslOKCMLCtsl高分一号全色影像参数设置(3)大气校正:选择Toolbox->RadiometricCorrection->AtmosphericCorrectionModule->FLAASHAtmosphericCorrection,弹出FLAASHAtmosphericCorrectionModelInputParameters对话框。要注意,全色影像不做大气校正,多光谱影像则需要做大气校正处理。

国Toolboxq-LtJ回Searchthet>>uLbc-i/E-adiorietricCc-rrection/AtnosphcricCorrectionNock,田曰E|।LDFavxsritesi:_lAnona1yDetection守口BandRatioE-jQChfltigtDetestimj_lLJClassificatiQn!_|FeslureExtraction再向FilterE+口Geonetil<Correclionj_lLJImaseSharrrninsSLlLiDAE?i!_lflossickLng由口Rachrj'_3GadiqictriCc?Tjec:tiOT।MApplyCainandOffset•FjAtnospherlcCorreetl&TiNodu.1e--•&QUlckktnosphtricCorrection(QTIAC):1-4JICalibrateAVKER•;-4HCalibrateTItfS■■Cr*st-TraekIlluiLLtjatiottCarrestIon•卜41DarkSubtiacti^n::■<«EFFORTPolishinE|:-EnisEivityAlphaFiesidiu1s4};EnlssivltyMornaliiatitnEnissivityReferen<eChannel451EnpiTicalLineCalibrateE^istinE

EmpiricalLineCimputeFacl^rs;;-4}iFlatFieldCalIbratlon\「事W工占RRtileetinceCaiibxatloti■■/LosResidualsCilitirati<>nRadicrjietricCalibration4iiThernalATnc-sphErlcCcirreellon由口RasterManaEeientddORegionsInterest币USFEAEiI_jSpectral;l_|Statistics-i-ir~iTaxsctDeteci1on国FLAA&HAtmoipJ-'i&rc国FLAA&HAtmoipJ-'i&rcCcrieiztianModel[rtputParairiMeT(4)点击InputRadianceImage,选择前面处理好的数据,在RadianceScaleFactors面板中选择Usesinglescalefactorforallbands,由于定标的辐射量数据与FLAASH的辐射亮度的单位相差10倍,所以在此Singlescalefactor选择:1,单击OK;■RadianeeScaleFactorsCReadarrayofscalefactors(1perband)fromASCIIfileQUsesinglescalefactorforallbandsSinglescalefactor1-OOODOdOKCancel(5)设置文件输出路径。(6)传感器基本信息设置:StJ-ne%ntenL.atE汁LKOttcfl[fl土)DJIS|阿口北fypeGF-1FllshrtBatetpr-14-州4SPlish-rT〔-eCHT]0.30Scns-srAltitudeUn)5C6.Lin31,泥GrcnmiE1edtlxgOn。0,016gd-11#5641F1IB151se(il)SceneCenterLocation从影像中自动获取;SensorType为GF-1;GroundElevation通过统计DEM数据获得;PixelSize根据相机选择,PMS相机全色2m,多光谱8m,WFV相机16m;FlightDate从影像xml头文件中读取,减去8换算成GMT时间;anii^-^uC.i_t■*ur■<二mrterTiiine>2014-04-1411:11:56^/CenterTime>(7)大气模型和气溶胶模型:AtmosphericModel根据经纬度和影像区域选择:

latrtLde(°N)Jan.h.brchh.^yJdi/SeptNdv.eoS/WS/WSAA/MLWNLWS/W70S/Ws/wMLWMLWMLWS顺tiLlWLWMLWMLWSASSASMLW50MLWMLWSASSASSAS5A540SASSASSASMLSMLSSAS30MLSMLSMLSITMLSZDTTTTTT1DITTTTTDTTTTTT-1DTTTITT-2DTTThilLSMLST-3DMLSMLSMLS帆L5MLSMLSSASS用SASSASSASSAS-5DSAS5ASSASMLWMLW5A5表2数据经纬度与获取时间决定选用的大气模型AerosolModel根据实际情况选择;AerosolRetrieval选择None;其他默认;(8)单击MultispectralSetting按钮,在FilterFunctionFile导入光谱响应曲线"gf-lpms.sli”,单击OK。地MultispectralSettingsIndextofirstband1SelectChannelDefinitionsbyFile©jGUICancelY/aterRetrievalKaufJian-TanreAerosolRetrievalWaterRefErencEChannelUndrfinedWaterAbsorptionChannel地MultispectralSettingsIndextofirstband1SelectChannelDefinitionsbyFile©jGUICancelY/aterRetrievalKaufJian-TanreAerosolRetrievalWaterRefErencEChannelUndrfinedWaterAbsorptionChannelF'iL1erFurLCtlonFile上:\」;:;]_1|11,在51:5、翻据臼化(勤:1居1遮普库-»-1\11111,1.511一]Undefined(9)单击AdvancedSettings,在高级设置中ModtranResolution选择5cm-1。■FLAASHAcv.ancedSettinge■FLAASHAcv.ancedSettingeSpectr&ErapTiDefinitic'nFileForDfcrrnadirLoginsInstrumentsDD<->I'¥S।1Zen3.1hAnsle]aQ口口-0口ierosolS^aleHeLsh.1(kn)1-50A2irmihAngle口。必00C02ffixinsEallo(ppm)390.00IJseScjji-arc3111FunEtL。(!口。|Use占比巧cun有yC<irrcctLoiifss[lt|ReusenoDTRAHCaleuLattansRq|,1TiledProcessing^es-tt|TileSize(Kb)100Radi或k。InageSpatialSubsetFullSceneRe-deflre舐式日Fa。totsF-jrRadianceIxageCh>>e-sbCutpjtReflectanceScaleFactor10000AuTenaticallySaveTenplateFile?es|ltOutput时羽nqwtiGFiles口中1W|Flcdtr-aaEezaluiian5ci-l-FlcdtraaerHoi=LScaltc口5berciDISORTStreans[.(10)设置好后,在大气校正模块面板中,单击Apply。1.2影像裁剪(1)File->Open,这里使用的是不规则裁剪,具体的理论知识可以在网上搜索,首先打开你矢量的"矢量.shp"数据,如图;(2)在Toolbox中,打开RegionsofInterest/SubsetDatafromROIs。SelectInputFile选择Beijing_TM.dat,点击OK,打开SubsetDatafromROIsParameters面板;(3)在SubsetDatafromROIsParameters面板中,设置以下参数:SelectInputROIs:选择EVF:矢量.shpMaskpixelsoutputofROI?:YesMaskBackgroundValue背景值:0

(4)选择输出路径和文件名,单击OK执行图像裁剪。1.3正射校正使用RPC文件进行正射校正,点开GeometricCorrection—Orthorectification—RPCOrthorectificationWorkflow,如下图所示:自•百ornetricCorrectirmfASCIIC»:iordinateC口工旧已工fMillGLTBuildG电口mornetricCorrectirmfASCIIC»:iordinateC口工旧已工fMillGLTBuildG电口m巳tryFiL^by;BuildSup"GLT,jeoreferencebySensorCjeLireferencefromGLTCjeLireferencefromKjfilfilapCoordirLSttfConverte:Orthorectification\--SfEuildRFCs-…-EFCUfthcixcctiEicwtiCEiguTiiiUEUrthorertit(±|"QjRegistrationRejifcijectGLTwithBowl3Super(jLTGeoreterenceSupei-IcjfilG吃口refere工四曰在该工作页面中将要处理的影像输入,然后点击next在弹出的页面的Advanced选项中选择输出像素大小,多光谱选择8m,全色选择2m,然后在该工作页面中将要处理的影像输入,然后点击next在弹出的页面的Advanced选项中选择输出像素大小,多光谱选择8m,全色选择2m,然后再export页面选择输出的位置,点击finish;RPCOrthorectificaticnRPCRefinementImportGZ?s团G再export页面选择输出的位置,点击finish;RPCOrthorectificaticnRPCRefinementImportGZ?s团G•莉dCirriageRezairiplirigBilinearGridSpacing10Freviev1.4图像融合m:ageSharpemTLg二JC1ISpectr:=ilSh:u:_jieri.iTLg二4C口1ororm:=J.izedl.BrovaHSVSharpEningFCSpTcMHLSharpEni工喧LiIlWLG&orr出+ricCorrsctiunFilterFestnj-h!Ektractionm:ageSharpemTLg二JC1ISpectr:=ilSh:u:_jieri.iTLg二4C口1ororm:=J.izedl.BrovaHSVSharpEningFCSpTcMHLSharpEni工喧LiIlWLG&orr出+ricCorrsctiunFilterFestnj-h!EktractionClaEsificationChangeDetectionB:±ndRatio十rlrin-rGr:=jti-Schm1ItP:itlSh:d2-peni2-e:=Lfchthe10I:1±uxw-CJo-cl-aaaa-Qh"oc■■■■■・・----…:0...s...s.:s.!s.!s:.s:.bfflfflAilqti:=lLvDetectiQn1.5图像镶嵌图像镶嵌也是所谓的羽化操作,一般来说一景高分影像能覆盖的区域是有限的,所以需要进行镶嵌操作,对其进行拼接到所需要的研究范围,点击工具箱中的Mosaicking-seamlessMosaicIPeaVji-eExtr3cti。工lS'C]FiL+.3S"C-ILimetricCoi_i_ectiorL由“匚二I工面演二亡Sharpening:由£jLiDAUS"ti5lHo互aicking丁・/仔FiselEasedflosaiekingSe:urileEEMosaic由.匚|Ej&itSe:urileEEMosaic£]■■QjJiadioirie+ricCcrre-2tioilJiasterM:driagenierut由“|Eegionso£InterestEFE虹:s-rispectr:il镶嵌过程最好两张两张镶嵌,这样电脑可以承受,先将要镶嵌的两张图通过绿色的加号加载进来,然后使用右上角的seamlines中的自动接边线进行自动拼接,再观察自动拼接的效果,如果有不好的地方再对其进行修改。接边线的修改主要就是画三角形,将不想要的部分去除,这个很简单,在网上也可以搜到。如何觉得对接边线进行修改比较麻烦,可以加我qq545698595

接边效果:1.6影像增强在遥感集市网站上有植被增强的插件,如果要做绿地信息提取的话可以去上面下载安装,然后对影像进行处理,但是很多情况下不需要这一步,对影像做一个拉伸处理便可。第二章绿地信息提取监督分类这一步实习课上都有,不需要演示面向对象绿地信息提取面向对象信息提取是目前最流行的信息提取方法,现在可以实现的软件有econition8.7ERDAS还有envi的EX模块,总体来说易康软件的面向对象提取方法比较完善,分割方法也比较多,但是因为软件需要收费,所以本文使用envi来进行处理。

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

图切割效果预览基于规则特征提取本文需要提取绿地,分类主要依靠的是绿地的植被覆盖指数。正常情况下,当NDVI>0便可以认为有植被存在,本文发现当令NDVI在0-1之间一些建筑物的房顶也被提了出来,所以将最低值提高到了0.1,虽然发现还是会有一些建筑物被提取,但是情况已经好了很多,其次因为建筑物大多是矩形的,所以可以通过设置矩形度对房屋进行剔除,因此对绿地的描述设置为:NDVI在0.1-1.0之间,rectangle在0.0182-0.8之间,两个属性的分类结果如下图(同一区域),光谱值和纹理提取出来的效果良好,将结果进行输出:图两种属性提取信息预览(左基于光谱,右基于几何)试验区绿地提取从分类结果可以看出分类的效果良好,将肉眼不能分辨的绿地提取了出来,适用于高分辨遥感影像的绿地信息提取,接着将其应用于全图:

图分类结果图精度评价精度评价有两种方法,在classic中的classification-postclassicification-confusionMatrix中的两种,一种是从高分辨率提取的图像验证低分辨率的图像,另一种是勾选正确的样本来验证用面向对象方法提取的影像信息,我使用了8m分辨率的多光谱影像提取结果和融合后影像的提取结果进行对比,来得出此次面向对象分类的准确度。..F1-291G10Uri9jperaisedGene-rafte-Random^mpleDe-d-skriTreeEnc^nernberCcJlertonCne-ra0上雪fromHC4ePast口号"市匚加0aA^?k|nClanColorsRd#ClacGlITArChangeDwectiDn0ta出th岫djcir匈..F1-291G10Uri9jperaisedGene-rafte-Random^mpleDe-d-skriTreeEnc^nernberCcJlertonCne-ra0上雪fromHC4ePast口号"市匚加0aA^?k|nClanColorsRd#ClacGlITArChangeDwectiDn0ta出th岫djcir匈/MinorhyMd洋is口UEp匚Ig£E«ECambreCla-siseEBijftrrZeroCorvfucioriHOCCunras.SrgFntfflEn]™gp□"fHkai口nId^Ktaf涂如出■>rvFileIBbmcTtoIe「击一布皿后口TErrafonFiFfawr即EraIMapVEkxTcjpogiraphkRadarWind™Hdp-Ll,ingGrwnd『力由limiago>UnjngGroundTruthROhLHM第三章矢量化处理3.1矢量切割在arcgis中打开.dat文件,然后把提取出来的绿地信息进行分级,我的论文中将绿地分成了:公园绿地、生产绿地、防护绿地、附属绿地和其他绿地5种绿地类型,首先将在ArcGIS中新建一个和.dat文件一样轮廓的面文件命名为五类,如图:tLJK■a5*CLUSSHIM:Ury蛆■EJir:D.XF2DD然后打开编辑器选择切割按钮将影像切割成几个大的部分,再对每个部分进行细部分割,通过与原影像叠加的目视解译和实地外调确定绿地类型。以上为五类绿地分割的情况,新建CLASS字段来对每一块区域进行赋值。3.2矢量处理将处理的.dat影像转换成shp文件:点开小红箱子-T中的数据管理工具datamanagementtools-features-featuretopolygon:-T中的数据管理工具data臼骸DataManagementTools田田田田田日Attachments田田田田田日DataComparisonDistributedGe-odatabaseDomainsFeatureClassFeatures区禽AddX¥Coordinates飞Adjust3DZ声飞BearingDistanceToLine气.CheckGeometry典%CopyFeatures典飞DeleteFeatures「Dice州飞FeatureEnvelopeToPolygonK'FeatureToLine“飞FeatureToPoint*飞、FeatureToPolygonFeatureVerticesToPoints声飞、Mini

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论