K清风影像信息提取之--面向对象特征提取流程

1、影像信息提取之-面向对象特征提取流程“同物异谱，同谱异物会对影像分类产生的影响，加上高分辨率影像的光谱信息不是很丰富， 还有经常伴有光谱相互影响的现象， 这对基于像素的分类方法提出了一种挑战， 面向对象的影像分类技术可以一定程度减少上述影响。本专题以ENVI中的面向对象的特征提取FX工具为例，对这种技术和处理流程做一个简单的介绍。本专题包括以下内容：面向对象分类技术概述ENVIFX简介ENVIFX操作说明1、面向对象分类技术概述面向对象分类技术集合临近像元为对象用来识别感兴趣的光谱要素， 充分利用高分辨率的全色和多光谱数据的空间，纹理，和光谱信息来分割和分类的特点，以高精度的分类结果或者矢量输

2、出。它主要分成两局部过程：对象构建和对象的分类。影像对象构建主要用了影像分割技术，常用分割方法包括基于多尺度的、基于灰度的、纹理的、基于知识的及基于分水岭的等分割算法。比拟常用的就是多尺度分割算法， 这种方法综合遥感图像的光谱特征和形状特征， 计算图像中每个波段的光谱异质性与形状异质性的综合特征值，然后根据各个波段所占的权重，计算图像所有波段的加权值， 当分割出对象或基元的光谱和形状综合加权值小于某个指定的阈值时，进行重复迭代运算， 直到所有分割对象的综合加权值大于指定阈值即完成图像的多尺度分割操作。影像对象的分类，目前常用的方法是“监督分类和“基于规那么知识分类。这里的监督分类和我们常说的监

3、督分类是有区别的， 它分类时和样本的比照参数更多，不仅仅是光谱信息，还包括空间、纹理等对象属性信息。基于规那么知识分类也是根据影像对象的属性和阈值来设定规那么进行分类。表1为三大类分类方法的一个大概的比照。传统地物的单个的中低分辨丰富的空类型根本原理影像的最小单元适用数据源缺陷表1传统基于光谱、基于专家知识决策树与基于面向对象的影像分类比照表基 于光 谱的 分类 方法基 于专 家知 识决 策树面 向对 象的 分类 方法光谱信息特征根据光谱特征、空间关系和其他上下文关系归类像元几何信息、结构信息以及光谱信息影像像元单 个 的影 像 像元一 个 个影 像 对象率 多 光 谱和 高 光 谱影像多源数

4、据中 高 分 辨率 多 光 谱和 全 色 影像间 信 息 利用 率 几 乎为零知识获取比拟复杂速度比拟慢2、ENVIFX简介全名叫“面向对象空间特征提取模块一FeatureExtraction，基于影像空间以及影像光谱特征，即面向对象，从高分辨率全色或者多光谱数据中提取信息，该模块可以提取各种特征地物如车辆、建筑、道路、桥、河流、湖泊以及田地等。该模块可以在操作过程中随时预览影像分割效果。 该项技术对于高光谱数据有很好的处理效果， 对全色数据一样适用。 对于高分辨率全色数据，这种基于目标的提取方法能更好的提取各种具有特征类型的地物。一个目标物体是一个关于大小、光谱以及纹理亮度、颜色等的感兴趣区

5、域。可应用于：从影像中尤其是大幅影像中查找和提取特征。添加新的矢量层到地理数据库输出用于分析的分类影像替代手工数字化过程一具有易于操作向导操作流程，随时预览效果和修改参数，保存参数易于下次使用和与同事共享，可以将不同数据源参加ENVIFX中DEMs、LiDARdatasets、shapefiles、地面实测数据以提高精度、交互式计算和评估输出的特征要素、提供注记工具可以标识结果中感兴趣的特征要素和对象等特点。3、ENVIFX操作说明ENVIFX的操作可分为两个局部：发现对象FindObject和特征提取Extractfeatures,如图1所示。图1FX操作流程示意图是，特征提取定义要素，规则

6、分类第监督分类“完成小rj|fjiidj是,要出导查看报告和统计下面在下的FX工具，以米的QB图像为例，介绍向对象信息提取的流程。下面我们。基于规那么的面向对象信息提取该方法的工具为Toolbox/FeatureExtraction/RuleBasedFeatureExtractionWorkflow第一步：准备工作根据数据源和特征提取类型等情况，可以有选择地对数据做一些预处理工作。空间分辨率的调整如果您的数据空间分辨率非常高，覆盖范围非常大，而提取的特征地物面积较大如云、大片林地等。可以降低分 辨 率 ， 提 供 精 度 和 运 算 速 度 。 可 禾 （J用Toolbox/RasterMa

7、nagement/ResizeData工具实现。光谱分辨率的调整如果您处理的是高光谱数据，可以将不用的波段除去。可利用Toolbox/RasterManagement/LayerStacking工具实现。多源数据组合当您有其他辅助数据时候，可以将这些数据和待处理数据组合成新的多波段数据文件，这些辅助数据可以是DEM,lidar影像，和SAR影像。当计算对象属性时候，会生成这些辅助数据的属性信息，可以提高信息提取精度。可利用Toolbox/RasterManagement/LayerStacking工具实现。空间滤波如果您的数据包含一些噪声，可以选择ENVI的滤波功能做一些预处理。这里直接在EN

8、VI中翻开图像文件。第二步：发现对象1启动RuleBasedFX工具在Toolbox中 ， 找 至 刊FeatureExtraction, 选 择/FeatureExtraction/RuleBasedFeatureExtractionWorkflow，翻开工作流的面板，选择待分类的影像，此外还有三个面板可切换：在InputMask面板可输入掩膜文件，在AncillaryData面板可输入其他多源数据文件，切换到CustomBands面板，有两个自定义波段，包括归一化植被指数或者波段比值、HSI颜色空间，这些辅助波段可以提高图像分割的精度，如植被信息的提取等自定义的属性，在Normalized

9、Difference和ColorSpace属性上打钩，如下列图所示，点击Next;图2输入数据和属性参数选择2影像分割、合并FX根据临近像素亮度、纹理、颜色等对影像进行分割，它使用了一种基于边缘的分割算法，这种算法计算很快，并且只需一个输入参数，就能产生多尺度分割结果。通过不同尺度上边界的差异控制，从而产生从细到粗的多尺度分割。选择高尺度影像分割将会分出很少的图斑，选择一个低尺度影像分割将会分割出更多的图斑，分割效果的好坏一定程度决定了分类效果的也triti+仆4(DataSdecUonDataSdecUonSdocinputHSdocinputH*s*s*HvnwladMtwnce*Hvnw

10、ladMtwnceLDIDF期ICB精确度，我们可以通过preview预览分割效果，选择一个理想的分割阀值，尽可能好地分割出边缘特征。有两个图像分割算法供选择：Edge,基于边缘检测，需要结合合并算法可以到达最正确效果；Intensity:基于亮度，这种算法非常适合于微小梯度变化如DEM、电磁场图像等，不需要合并算法即可到达较好的效果。调整滑块阀值对影像进行分割，这里设定阈值为40。注：按钮词是用来选择分割波段的，默认为BaseImage所有波段。影像分割时，由于阈值过低，一些特征会被错分，一个特征也有可能被分成很多局部。我们可以通过合并来解决这些问题。合并算法也有两个供选择：FullLamb

11、daSchedule,合并存在于大块、 纹理性较强的区域，如树林、云等，该方法在结合光谱和空间信息的根底上迭代合并邻近的小斑块；FastLambda:合并具有类似的颜色和边界大小相邻节段。设定一定阈值，预览效果,这里我们设置的阈值为90,点Next进入下一步。TextureKernalSize:纹理内核的大小）如果数据区域较大而纹理差异较小， 可以把这个参数设置大一点。 默认是3,最大是19。注：这一步是可选项，如果不需要可以按照默认的0直接跳过。图3图像分割、合并这时候FX生成一个RegionMeans影像自动加载图层列表中，并在窗口中显示，它是分割后的结果，每一块被填充上该块影像的平均光谱

12、值。接着进行下一步操作。目前，已经完成了发现对象的操作过程，接下来是特征的提取。第三步：根据规那么进行特征提取在规那么分类界面。每一个分类有假设干个规那么Rule组成，每一个规那么有假设干个属性表达式来描述。规那么与规那么之间是与的关系，属性表达式之间是并的关系。同一类地物可以由不同规那么来描述，比方水体，水体可以是人工池塘、湖泊、河流，也可以是自然湖泊、河流等，描述规那么就不一样，需要多条规那么来描述。每条规那么又有假设干个属性来描述，如下是对水的一个描述：面积大于500像素延长线小于NDVI小于对道路的描述：延长线大于紧密度小于标准差小于20这里以提取居住房屋为例来说明规那么分类的操作过程

13、。首先分析影像中容易跟居住房屋错分的地物有：道路、森林、草地以及房屋旁边的水泥地。点击遁按钮，新建一个类别，在右侧Classproperties下修改好类别的相应属性。图4规那么分类面板1）第一条属性描述，划分植被覆盖和非覆盖区。在默认白属性SpectralMean上单击，激活属性，右边出现属性选择面板， 如下图。 选择Spectral,Band下面选择NormalizedDifferenceo在第一步自定义波段中选择的波段是红色和近红外波段，所以在此计算的是NDVIc把ShowAttributeImage勾上，可以看到计算的属性图像。通过拖动滑条或者手动输入确定阈值。 在阈值范围内的在预览窗

14、口里显示为红色，在Advanced面板，有三个类别归属的算法：算法有二进制、线性和二次多项式。选择二进制方法时，权重为0或者1,即完全不匹配和完全匹配两个选项；中选择线性和二次多项式时，可通过Tolerance设置匹配程度，值越大，其他分割块归属这一类的可能性就越大。这里选择类别归属算法为Liner,分类阈值Tolerance为默认的5,如下列图图5对象属性面板图6归属类别算法和阈值设置2）第二条属性描述，剔除道路干扰居住房屋和道路的最大区别是房屋是近似矩形）我们可以设置Rectangularfit属性。在Rule上右键选择AddAttibute按钮，新建一个规那么）在右侧Type中选择Spa

15、tial）在Name中选择Rectangularfit。设置值的范围是0.51）其他参数为默认值。注：预览窗口默认是该属性的结果，点击AllClasses,可预览几个属性共同作用的结果。同样的方法设置Type:Spatial;Name:AreaArea45Type:Spatial;Name:ElongationElongation33）第三条属性描述，剔除水泥地干扰水泥地反射率比拟高，居住房屋反射率较低，所以我们可以设置波段的象元值。Type:spectral;Name:SpectralMean）Band:GREENSpectralMean（GREEN）RadialBasisSigmoid,需

16、要设置GammainKernelFunction参数。这个值是一个大于零的浮点型数据。 默认值是输入图像波段数的倒数。b为SVM规那么指定thePenalty参数，这个值是一个大于零的浮点型数据。这个参数控制了样本错误与分类刚性延伸之间的平衡，默认值是100。AllowUnclassified是允许有未分类这一个类别，将不满足条件的斑块分到该类，默认是允许有未分类的类别。Threshold为分类设置概率域值，如果一个像素计算得到所有的规那么概率小于该值，该像素将不被分类，范围是0100,默认是5。主成分分析是比拟在主成分空间的每个分割对象和样本，将得分最高的归为这一类。这里我们选择K邻近法，K

17、参数设置为5,点击Next,输出结果。最终结果的输出方法和基于规那么的一样。直接输出矢量该方法的工具为Toolbox/FeatureExtraction/SegmentOnlyFeatureExtractionWorkflow。操作方法参考前面的第一和第二步骤， 第三步直接选择路径输出分割栅格结果和矢量结果。从以上的实际操作可以看到，ENVIFX扩展模块操作具有易于操作向导操作流程，随时预览效果和修改参数。4、小结基于像元的分类方法，依据主要是利用像元的光谱特征，大多应用在中低分辨率遥感图像。而高分辨率遥感图像的细节信息丰富，图像的局部异质性大，传统的基于像元的分类方法易受高分辨率影像局部异质

18、性大的影响和干扰。而面向对象分类方法可以高分辨率图像丰富的光谱、形状、结构、纹理、相关布局以及图像中地物之间的上下文信息，可以结合专家知识进行分类，可以显著提高分类精度，而且使分类后的图像含有丰富的语义信息，便于解译和理解。对高分辨率影像来说，还是一种非常有效的信息提取方法，具有很好的应用前景。5、ENVIFX操作说明补充ENVIFX的操作可分为两个局部：发现对象FindObject和特征提取Extractfeatures。3.1准备工作根据数据源和特征提取类型等情况，可以有选择的对数据做一些预处理工作。l空间分辨率的调整如果您的数据空间分辨率非常高， 覆盖范围非常大，而提取的特征地物面积较大

19、如云、大片林地等。可以降低分辨率，提供精度和运算速度。可利用ENVI主界面-BasicTool-ResizeData工具实现。l光谱分辨率的调整如果您处理的是高光谱数据， 可以将不用的波段除去。可利用ENVI主界面-BasicTool-layerstacking工具实现。l多源数据组合当您有其他辅助数据时候， 可以将这些数据和待处理数据组合成新的多波段数据文件，这些辅助数据可以是DEM,lidar影像,和SAR影像。 当计算对象属性时候，会生成这些辅助数据的属性信息，可以提高信息提取精度。可利用ENVI主界面-BasicTool-layerstacking工具实现。l空间滤波如果您的数据包含一

20、些噪声，可以选择ENVI的滤波功能做一些预处理。3.2发现对象（一）翻开数据在ENVIZoom中翻开ProcessingFeatureExtraction。如图2所示，BaseImage必须要选择，辅助数据AncillaryData和掩膜文件MaskFile是可选。这里选择ENVI自带数据envidatafeature_extractionqb_colorado,它是米的快鸟数据。（二）影像分割FX根据临近像素亮度、纹理、颜色等对影像进行分割，它使用了一种基于边缘的分割算法，这种算法计算很快，并且只需一个输入参数，就能产生多尺度分割结果。通过不同尺度上边界的差异控制，从而产生从细到粗的多尺度分

21、割。选择高尺度影像分割将会分出很少的图斑，选择一个低尺度影像分割将会分割出更多的图斑，分割效果的好坏一定程度决定了分类效果的精确度，我们可以通过预览分割效果，选择一个理想的分割阀值，尽可能好地分割出边缘特征。调整滑块阀值对影像进行分割，这里设定阈值为30,点击Next按钮，这时候FX生成一个RegionMeans影像自动加载图层列表中，并在窗口中显示，它是分割后的结果，每一块被填充上该块影像的平均光谱值。接着进行下一步操作。注：按钮是用来选择分割波段的，默认为BaseImage所有波段。（三）合并分块影像分割时，由于阈值过低，一些特征会被错分，一个特征也有可能被分成很多局部。我们可以通过合并来

22、解决这些问题。FX利用了FullLambda-Schedule算法。这一步是可选项，如果不需要可以直接跳过。设定一定阈值，预览效果。这里我们设置的阈值为95,点Next进入下一步。（四）分块精炼FX提供了一种阈值法Thresholding进一步精炼分块的方法。对于具有高比照度背景的特征非常有效例如，明亮的飞机对黑暗的停机坪。可以将精炼结果生成掩膜图层Mask，按钮可以修改基于哪个波段。这里我们就直接选择NoThresholdingdefault,点击Next进入下一步操作。（五）计算对象属性计算4个类别的属性：光谱、空间、纹理、自定义颜色空间和波段比。其中颜色空间选择三个RGB波段转换为HSI

23、颜色空间，波段比选择两个波段用于计算波段比常用红色和近红外波段。各个属性的详细描述参考ENVI/IDL提供的文档。这里我们按照默认全选择，ColorSpace选择RGB,BandRatio选择红色和近红外波段，点击Next按钮进行下一步操作。目前，已经完成了发现对象的操作过程，接下来是特征的提取。特征提取如图7所示，有三种特征提取方法供选择，分别是监督分类、规那么分类和直接矢量输出。（一）输出矢量选择ExportVectors,进入图8界面，选择保存路径，属性信息也可选择输出输出完成会出来一个报表。不关闭FX浮动面板，在ENVIZoom中将得到的矢量特征加载显示。点击Previous按钮，回到

24、图7界面。（二）监督分类在图7界面中选择Classifybyselectionexamples,下步至！如图8所示界面。1）选择样本在ENVIZoom中，切换到Select方式，双击Feature，翻开一个类别的属性，如图10所示，修改显示颜色、名称等信息。在分割图上选择一些样本，为了方便样本的选择，可以在ENVIZoom的图层管理中将原图移到最上层，选择一定数量的样本，如果错选样本，可以在这个样本上点击左键删除。一个类别的样本选择完成之后，新增类别，用同样的方法修改类别属性和选择样本。在选择样本的过程中，可以随时预览结果。可以把样本保存为xml文件以备下次使用。2）设置样本属性在图9中，切换

25、到Attributes选项。默认是所有的属性都被选择，可以根据提取的实际地物特性选择一定的属性。这里我们按照默认全部选择。3）选择分类方法在图9中，切换到Algorithm选项。FX提供了两种分类方法：K邻近法KNearestNeighbor和支持向量机SupportVectorMachine,SVM，如图12所示。这里我们选择K邻近法，K参数设置为5,点击下一步，输出结果。4）输出结果特征提取结果可以以两种格式输出，矢量和图像，如图12所示。矢量可以是所有分类以单个文件输出或者每一个类别分别输出；图像可以把分类结果和规那么结果分布输出。这里我们选择单个文件以及属性数据一块输出，分类图像和规那

26、么图像一块输出。点击Next按钮完成输出， 同时可以看到整个操作的参数和结果统计报表。（三）规那么分类在图7界面中选择Classifybycreatingrules,点击Next,到图15规那么分类界面。每一个分类有假设干个规那么Rule组成，每一个规那么有假设干个属性表达式来描述。规那么与规那么直接是与的关系，属性表达式之间是并的关系。同一类地物可以由不同规那么来描述，比方水体，水体可以是人工池塘、湖泊、河流，也可以是自然湖泊、河流等，描述规那么就不一样，需要多条规那么来描述。每条规那么又有假设干个属性来描述，如下是对水的一个描述：l面积大于500像素l延长线小于lNDVI小于对道路的描述：

27、l延长线大于l紧密度小于l标准差小于20这里以提取居住房屋为例来说明规那么分类的操作过程。首先分析影像中容易跟居住房屋错分的地物有：道路、森林、草地以及房屋旁边的水泥地。双击Feature_1图标，修改好类别的相应属性。1）第一条属性描述，划分植被覆盖和非覆盖区双击rule,翻开对象属性选择面板，如图16所示。选择Customized-bandratio。FX会根据选择的波段情况技术波段比值，比方这里在属性计算步骤中选择的RatioBand是红色和近红外波段，所以此时计算的是NDVI。把ShowAttributeImage勾上， 可以看到计算的结果， 通过ENVIZoom工具查看各个分割块对应

28、的值。点击Next按钮，或者双击bandratio,进入bandratio属性设置对话框.通过拖动滑条或者手动输入确定阈值。FuzzyTolerance是设置模糊分类阈值，值越大，其他分割块归属这一类的可能性就越大。归类函数有线性和S-type两种。 这里设置模糊分类阈值为默认的5,归属类别为S-type,值的范围为， 勾选ShowRuleConfidenceImage可以预览规那么图像。点击Ok完成此条属性描述。2）第二条属性描述，去除道路影响居住房屋和道路的最大区别是房屋是近似矩形， 我们可以设置Rect_fit属性。 点击按钮或者双击rule,选择Spatial-rect_fit。设置值

29、的范围是0.51,其他参数为默认值。同样的方法设置Spatial-Area:FuzzyTolerance=0,90Areaelongation（延长）:elongationavgband_2:avgband_250。类似的思路可以提取道路、林地、草地等分类，这里就不一一例举。最终结果的输出方式和监督分类一样。批处理操作ENVI提供了ENVI_FX函数，具体语法如下：ENVI_DOIT,ENVI_FX_DOIT,A_FID=array,A_POS=array,BR_BANDS=array,CENTERLINE_OPTIONS=array,CONF_THRESHOLD=floatingpoint,

30、CS_BANDS=array,DIMS=array,FID=fileID,/EXPORT_ATTRIBUTES,/EXPORT_RASTER,/INVERSE_MASK,KERNEL_SIZE=longinteger,M_FID=fileID,MERGE_LEVEL=floatingpoint,POS=array,R_FID=variable,RASTER_FILENAME=stringorstringarray,/RAW_ATTRIBUTES儿RAW_FILENAME=string,REFINE_BAND=integer,/REFINE_INVERSE,/REFINE_MASK,RULESE

31、T_FILENAME=string,SCALE_LEVEL=floatingpoint,SEGMENT_BANDS=array,SMOOTHING_THRESHOLD=floatingpoint,TD_FILENAME=string,THRESHOLD_LOWER=floatingpoint,THRESHOLD_UPPER=floatingpoint,VECTOR_FILENAME=string,VECTOR_OPTIONS=stringarray详细的变量说明请参考ENVI/IDL提供的。利用此函数可以对根据分割参数和分类参数编写批处理，或者将此功能应用到别的系统上。从以上的实际操作可以看到，ENVIFX扩展模块操作具有易于操作向导操作流程，随时预览效果和修改参数；功能比拟强大，具有多种特征提取方法，包括监督分类和知识规那么分类，等等。附录对象属性说明：1Spatial属性属性描述AREA