2022年2022年监督分类和面向对象分类流程 .pdf

《2022年2022年监督分类和面向对象分类流程 .pdf》由会员分享，可在线阅读，更多相关《2022年2022年监督分类和面向对象分类流程 .pdf（21页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

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

2、载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 21 页 - - - - - - - - - 图二图三为了加快数据处理的速度，本文是选择先进行辐射定标然后将图像裁剪在进行后续的操作，预处理流程如下图（图四）：图四1.1 辐射校正分为辐射定标和大气校正名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 21 页 - - - - - - - - - （1）打开数据： ENVI-Open As-C

3、RESDA-GF-1,选择处理的影像，打开XML 后缀文件；（2）辐射定标：选择Toolbox-Radiometric Correction- Radiometric Calibration ，选择待处理的高分数据弹出 Radiometric Calibration对话框， 进行如图设置。 对于多光谱影像，点击 Apply FLAASH Setting 设 置 成 默 认 值 ； 如 果 是 对 全 色 影 像 进 行 辐 射 定 标 ， 那 么Calibration则 是Reflectance,Out Put Type 为 UInt,Scale Factor 为 1000，如下图：名师资料总

4、结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 21 页 - - - - - - - - - 高分一号多光谱影像参数设置高分一号全色影像参数设置（ 3 ） 大 气 校 正 ： 选 择Toolbox-Radiometric Correction-Atmospheric Correction Module-FLAASH Atmospheric Correction， 弹出 FLAASH Atmospheric Correction Model Input Parameters对话框。

5、要注意，全色影像不做大气校正，多光谱影像则需要做大气校正处理。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 21 页 - - - - - - - - - （4）点击 Input Radiance Image，选择前面处理好的数据，在Radiance Scale Factors面板中选择 Use single scale factor for all bands ， 由于定标的辐射量数据与FLAASH 的辐射亮度的单位相差 10 倍，所以在此Single scale fa

6、ctor 选择： 1，单击 OK；名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 21 页 - - - - - - - - - （5）设置文件输出路径。（6）传感器基本信息设置：Scene Center Location 从影像中自动获取；Sensor Type为 GF-1；Ground Elevation 通过统计 DEM 数据获得；Pixel Size 根据相机选择，PMS 相机全色2m，多光谱8m，WFV 相机 16m；Flight Date 从影像 xml 头文件

7、中读取，减去8 换算成 GMT 时间；（7）大气模型和气溶胶模型：Atmospheric Model 根据经纬度和影像区域选择：名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页，共 21 页 - - - - - - - - - Aerosol Model 根据实际情况选择；Aerosol Retrieval 选择 None；其他默认；（8） 单击 Multispectral Setting 按钮， 在 Filter Function File 导入光谱响应曲线 “gf-1pms

8、.sli ” ，单击 OK 。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 7 页，共 21 页 - - - - - - - - - （9）单击 Advanced Settings，在高级设置中Modtran Resolution 选择 5 cm-1。（10）设置好后，在大气校正模块面板中，单击Apply 。1.2 影像裁剪（1）File - Open ，这里使用的是不规则裁剪，具体的理论知识可以在网上搜索，首先打开你矢量的 矢量 .shp 数据，如图；（2）在 Toolbox 中

9、，打开 Regions of Interest /Subset Data from ROIs。Select Input File选择 Beijing_TM.dat，点击 OK ，打开 Subset Data from ROIs Parameters面板；（3）在 Subset Data from ROIs Parameters面板中，设置以下参数：Select Input ROIs：选择 EVF:矢量 .shp Mask pixels output of ROI？： Yes Mask Background Value背景值： 0 名师资料总结 - - -精品资料欢迎下载 - - - - - -

10、 - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 8 页，共 21 页 - - - - - - - - - （4）选择输出路径和文件名，单击OK执行图像裁剪。1.3 正射校正使 用RPC 文 件 进 行 正 射 校 正 ， 点 开Geometric Correction Orthorectification RPC Orthorectification Workflow，如下图所示：在该工作页面中将要处理的影像输入，然后点击next 在弹出的页面的Advanced 选项中选择输出像素大小，多光谱选择8m，全色选择2m，然后名师资料总结 - - -精

11、品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 9 页，共 21 页 - - - - - - - - - 再 export 页面选择输出的位置，点击finish；1.4 图像融合对每一景的多光谱和全色影像进行正射校正之后就对其进行融合，一般来说高分一号影像使用 GS融合， 点击工具箱中的Image SharpeningGram-Schmidt Pan Sharpening 分别将多光谱和全色影像输入，然后设置输出位置，finish；名师资料总结 - - -精品资料欢迎下载 - - - - - - - - -

12、- - - - - - - - - 名师精心整理 - - - - - - - 第 10 页，共 21 页 - - - - - - - - - 1.5 图像镶嵌图像镶嵌也是所谓的羽化操作，一般来说一景高分影像能覆盖的区域是有限的，所以需要进行镶嵌操作， 对其进行拼接到所需要的研究范围，点击工具箱中的Mosaickingseamless Mosaic 镶嵌过程最好两张两张镶嵌，这样电脑可以承受，先将要镶嵌的两张图通过绿色的加号加载进来，然后使用右上角的seamlines 中的自动接边线进行自动拼接，再观察自动拼接的效果，如果有不好的地方再对其进行修改。接边线的修改主要就是画三角形，将不想要的部分去

13、除，这个很简单， 在网上也可以搜到。如何觉得对接边线进行修改比较麻烦，可以加我qq545698595 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 11 页，共 21 页 - - - - - - - - - 接边效果 : 1.6 影像增强在遥感集市网站上有植被增强的插件，如果要做绿地信息提取的话可以去上面下载安装，然后对影像进行处理，但是很多情况下不需要这一步，对影像做一个拉伸处理便可。第二章绿地信息提取2.1 监督分类这一步实习课上都有，不需要演示2.2 面向对象绿地信息提取面向

14、对象信息提取是目前最流行的信息提取方法，现在可以实现的软件有econition8.7 ERDAS还有 envi 的 EX模块， 总体来说易康软件的面向对象提取方法比较完善，分割方法也比较多，但是因为软件需要收费，所以本文使用envi 来进行处理。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 12 页，共 21 页 - - - - - - - - - 2.2.1 分类步骤2.2.2 实现对象的创建影像对象的创建主要进行了图像分割处理，本文采用多尺度分割算法对影像进行分割。首先要对数据

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

16、- - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 13 页，共 21 页 - - - - - - - - - 图 切割效果预览基于规则特征提取本文需要提取绿地，分类主要依靠的是绿地的植被覆盖指数。正常情况下，当NDVI0便可以认为有植被存在，本文发现当令NDVI 在 0-1 之间一些建筑物的房顶也被提了出来，所以将最低值提高到了0.1， 虽然发现还是会有一些建筑物被提取，但是情况已经好了很多，其次因为建筑物大多是矩形的，所以可以通过设置矩形度对房屋进行剔除，因此对绿地的描述设置为： NDVI 在 0.1-1.0 之间， recta

17、ngle 在 0.0182-0.8 之间，两个属性的分类结果如下图（同一区域） ，光谱值和纹理提取出来的效果良好，将结果进行输出：图两种属性提取信息预览（左基于光谱，右基于几何）试验区绿地提取从分类结果可以看出分类的效果良好，将肉眼不能分辨的绿地提取了出来，适用于高分辨遥感影像的绿地信息提取，接着将其应用于全图：名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 14 页，共 21 页 - - - - - - - - - 图 分类结果图2.2.3 精度评价精度评价有两种方法，在class

18、ic 中的 classification-post classicification-confusion Matrix中的两种， 一种是从高分辨率提取的图像验证低分辨率的图像，另一种是勾选正确的样本来验证用面向对象方法提取的影像信息，我使用了 8m 分辨率的多光谱影像提取结果和融合后影像的提取结果进行对比，来得出此次面向对象分类的准确度。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 15 页，共 21 页 - - - - - - - - - 第三章矢量化处理3.1 矢量切割在 ar

19、cgis 中打开 .dat 文件，然后把提取出来的绿地信息进行分级，我的论文中将绿地分成了：公园绿地、 生产绿地、 防护绿地、 附属绿地和其他绿地5 种绿地类型， 首先将在 ArcGIS中新建一个和 .dat 文件一样轮廓的面文件命名为五类，如图：然后打开编辑器选择切割按钮将影像切割成几个大的部分，再对每个部分进行细部分割，通过与原影像叠加的目视解译和实地外调确定绿地类型。以上为五类绿地分割的情况，新建CLASS 字段来对每一块区域进行赋值。3.2 矢量处理将处理的 .dat 影像转换成shp 文件：点开小红箱子中的数据管理工具data management tools-features-fe

20、ature to polygon:将面向对象提取的影像转成矢量文件名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 16 页，共 21 页 - - - - - - - - - 如下图，是转换出来的全部面文件，需要打开编辑器，将最外面的面文件删除，剩下才是我们真正需要的面，删除后的矢量文件在下图：名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 17 页，共 21 页 - - -

21、 - - - - - - 点击 Arcmap10.1 界面的第一行selection 的 select by attributes ， 打开页面之后选择CLASS字段，以生产绿地为例，其类别为1，以下高亮区域便是生产绿地的范围：接着点击 selection 中的 select by location ，目标图层选择自己用影像转出来的面文件，源图层选择你切割分类的图层，点击apply，分类结果导入了目标图层，以下高亮部分便是选择出来的生产绿地图斑，名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - -

22、 - 第 18 页，共 21 页 - - - - - - - - - 打开目标图层litongchangshi 的属性表， 选择右下方右边的图标，添加一个类别字段，把选择的图斑赋值为生产绿地，依次操作下去便可把所有图斑的绿地类型赋值完毕，如下图名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 19 页，共 21 页 - - - - - - - - - 最后再利用CLASS_ID 进行显示，最后可以分出五种绿地类型的区域分布名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 20 页，共 21 页 - - - - - - - - - 第四章绿地信息分析与评价使用arcgis 的属性选择器可以对绿地进行武汉市的绿地空间分布计算、绿量计、各城区绿地率计算等等，在此不讨论，如有其它阐述不足之处，请联系qq:545698595, 欢迎批评指正！名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 21 页，共 21 页 - - - - - - - - -