基于面向对象-随机森林分类法和中等分辨率遥感影像的沼泽湿地制图的方法技术

基于面向对象-随机森林分类法和中等分辨率遥感影像的沼泽湿地制图的方法，本发明专利技术涉及一种利用遥感影像进行沼泽湿地制图的方法。本发明专利技术的目的是要解决目前遥感分类方法，仅利用光谱信息，不能准确的提取草丛沼泽、灌丛沼泽和森林沼泽空间信息的技术问题。本发明专利技术利用对Landsat-8OLI数据进行辐射定标、大气校正、正射纠正、几何精纠正和投影转换；并利用eCognition 8.64软件对数据进行处理，将非沼泽湿地对象去除保留沼泽湿地对象，通过设定指数阈值区分出草丛沼泽、灌丛沼泽和森林沼泽并在ArcGIS 10.0软件中制作沼泽湿地专题地图。利用本发明专利技术的制图方法精度高。本发明专利技术用于沼泽湿地制图。

【技术实现步骤摘要】

本专利技术涉及一种利用遥感影像进行沼泽湿地制图的方法。
技术介绍
湿地是地球上单位面积生态服务价值最高、固碳能力最强、生物多样性保护意义 最大的生态系统，被誉为"地球之肾"。沼泽湿地（MarshWetland)是最主要的湿地类型， 包括草丛沼泽、灌丛沼泽和森林沼泽。它具有涵养水源、蓄洪防旱、降解污染、维系生物多样 性、调节河川径流和局地气候、补充地下水、净化环境和沉积营养物质、防止土壤侵蚀等多 种功能；沼泽湿地还可以为人类提供食品、水果、药材、能源和多种工业原料；沼泽湿地具 有的食物链及其所支撑的生物多样性，为许多野生动物提供独特的生存空间。因此保护沼 泽湿地的生态环境，对保护区域生态平衡具有重要的作用。 面向对象-随机森林的遥感图像解译方法是相对于传统遥感影像处理软件主要 针对单个像元的解译算法而言的。该方法在分类时不仅考虑地物的光谱特征，还主要利用 其几何特征和结构特征，图像中的最小单元不再是单个的像元，而是一个个对象。该方法突 破了传统遥感影像分类方法以像元为基本分类和处理单元的局限性，以含有更多语义信息 的多个相邻像元组成的对象为处理单元，可以实现较高层次的遥感图像分类和目标地物提 取。该方法是基于认知模型的遥感信息提取方法，更贴近人类的认知过程，已成为遥感信息 提取领域主要的研究方向之一。 2〇13年2月11日，美国航空航天局（NASA)成功发射Landsat-8卫星。Landsat-8 卫星上携带两个传感器，分别是0LI陆地成像仪（OperationalLandImager)和TIRS 热红外传感器（ThermalInfraredSensor)。0LI陆地成像仪有9个波段，成像宽幅为 185Xl85km。本专利技术选取Landsat-80LI陆地成像仪的前8个波段（0.43~0·45μπι， 0. 45 ~0. 51μm，0. 53 ~0. 59μm，0. 64 ~0. 67μm，0. 85 ~0. 88μm，1. 57 ~1. 65μm， 2. 11~2. 29μm，0. 50~0. 68μm)，其中波段1-7分辨率为30米，波段8为15米分辨率的 全色波段，卫星每16天可以实现一次全球覆盖。已被证实非常适合于沼泽湿地空间信息提 取研究。 遥感图像为沼泽湿地分布信息的获取提供了数据基础，为采取有效措施加强沼泽 湿地资源的科学保护和管理，合理开展沼泽湿地资源的开发和利用，以及有效维护和加强 沼泽湿地的生态环境及其功能，发挥其生态价值等方面提供理论和实践基础。然而准确地 利用遥感数据提取草丛沼泽、灌丛沼泽和森林沼泽信息仍是一个难以克服的问题。
技术实现思路
本专利技术的目的是要解决目前遥感分类方法，仅利用光谱信息，不能准确的提取草 丛沼泽、灌丛沼泽和森林沼泽空间信息的技术问题，提供了一种基于面向对象-随机森林 分类方法和中等分辨率遥感影像的草丛沼泽、灌丛沼泽和森林沼泽区分方法。 按照以下步骤进行： -、对Landsat-80LI数据进行辐射定标、大气校正、正射纠正、几何精纠正和投影 转换； 二、收集基于数字高程模型（DEM)获取的坡度、地形湿度指数和基于缨帽变换获 取的亮度、绿度和湿度辅助数据； 三、在eCognition8. 64软件支持下，对Landsat-80LI影像和辅助分类数据进行 多层多尺度分割，eCognition8. 64软件自带多尺度分割算法，可将遥感影像分割为一系列 均质单元，即对象，得到一系列分割单元，每个分割单元由空间上相邻、同质性达到80%以 上的像元组成，将每个分割单元作为一个对象； 多尺度分割设置参数 四、提取步骤三中所得到的每个对象的光谱特征、纹理特征、形状特征和结构特 征，并且利用归一化植被指数（NDVI)、归一化水体指数（NDWI)、比值植被指数（RVI)、差值 植被指数（DVI)、坡度（Slope)、地形湿度指数、HIS变换的色调（Hue)特征以及缨帽变换的 亮度（BI)、绿度（GVI)和湿度（WI)等计算对象的算数特征； 主要植被指数计算公式如（1)-(4): DVI = NIR-R (4) 式中，NIR为近红外波段，波长为0.85~0.88μm;R为红色波段，波长为0.64~ 0. 67μm;G为绿色波段，波长为0. 53~0. 59μm。 五、利用步骤四中得到的所有特征信息，确定区分沼泽湿地与非沼泽湿地的特征 值，从Landsat_80LI影像中去除非沼泽湿地对象； 提取过程和参数如下： 1、设定水体指数阈值，将水体剔除，参数设置为：BK=0.24； 2、设定林地指数阈值，将林地剔除，参数设置为：RVI> = 9. 5 或者HueR:G:B=SWIR1:NIR:R> =0·25; 3、设定建筑用地指数阈值，将建筑用地剔除，参数设置为： GVI<=0.3； 4、设定沼泽湿地指数阈值，将非沼泽湿地剔除，参数设置为： Slope< = 5 ； 六、在步骤五所提取的沼泽湿地对象基础上，利用步骤四所涉及的草丛沼泽、灌丛 沼泽和森林沼泽的特征逐步判定并提取草丛沼泽、灌丛沼泽和森林沼泽对象； 提取过程和参数如下： 1、设定森林沼泽指数阈值，提取森林沼泽，参数设置为： NDVI/R) = 0. 0021 ； 2、设定灌丛沼泽指数阈值，提取灌丛沼泽，参数设置为： StandarddeviationBK= 0. 0055 ； 3、设定草丛沼泽指数阈值，提取草丛沼泽，参数设置为： StandarddeviationBI> = 0. 0055 ； 七、导出草丛沼泽、灌丛沼泽和森林沼泽对象，获得沼泽湿地矢量数据； 八、在ArcGIS10. 0软件中制作沼泽湿地专题地图。 本专利技术的优点： 本专利技术将面向对象-随机森林的遥感图像分类方法应用于草丛沼泽、灌丛沼泽和 森林沼泽的提取中，将独立的像元合并成为同质的对象，提取草丛沼泽、灌丛沼泽和森林沼 泽信息的过程中不仅利用光谱特征，还利用了对于区分草丛沼泽、灌丛沼泽和森林沼泽十 分重要的纹理特征和拓扑特征，逐步得到草丛沼泽、灌丛沼泽和森林沼泽空间分布信息。所 得分类结果精度高，且具有明确的地理意义。本专利技术克服了传统分类技术提取草丛沼泽、灌 丛沼泽和森林沼泽信息时存在的"漏分"和"误分"现象严重的困难，同时也解决了分类得到 的草丛沼泽、灌丛沼泽和森林沼泽空间信息存在"椒盐现象"，不具有明确地理意义等问题。 本专利技术能够简捷、快速、准确地提取草丛沼泽、灌丛沼泽和森林沼泽信息以及制作沼泽湿地 专题地图。本专利技术对草丛沼泽、灌丛沼泽和森林沼泽空间信息提取结果的生产精度、用户精 度和Kappa系数分别为：90. 05 %，93. 30 %，0. 89,能够满足沼泽湿地制图需求。本专利技术克服 了传统沼泽湿地制图过程中存在的费时、费力等缺点，对中低分辨率遥感影像提取草丛沼 泽、灌丛沼泽和森林沼泽具有实践意义。【附图说明】 图1案例区域的面向对象-随机森林方法的遥感影像分类流程图 图2案例区域的标准假彩色合成图像，依据加色法彩色合成原理，选择 Landsat-80LI影像的NIR、R、G三个波段，分别赋予红、绿、蓝三种原色，合成假彩色图像。 图3案例区域的沼泽湿地专题图。其中1为草丛沼泽，2为灌丛沼泽，3为本文档来自技高网...

【技术保护点】
基于面向对象‑随机森林分类法和中等分辨率遥感影像的沼泽湿地制图的方法，其特征在于该制图方法按照以下步骤进行：一、对Landsat‑8 OLI数据进行辐射定标、大气校正、正射纠正、几何精纠正和投影转换；二、收集基于数字高程模型(DEM)获取的坡度、地形湿度指数和基于缨帽变换获取的亮度、绿度和湿度辅助数据；三、在eCognition 8.64软件支持下，对Landsat‑8 OLI影像和辅助分类数据进行多层多尺度分割，得到一系列分割单元，每个分割单元由空间上相邻和同质性达到80％以上的像元组成，将每个分割单元作为一个对象；多尺度分割设置参数分割尺度彩色因子形状因子光滑度紧致度500.80.20.60.4四、提取步骤三中所得到的每个对象的光谱特征、纹理特征、形状特征和结构特征，并且利用归一化植被指数(NDVI)、归一化水体指数(NDWI)、比值植被指数(RVI)、差值植被指数(DVI)、坡度(Slope)、地形湿度指数、HIS变换的色调(Hue)特征以及缨帽变换的亮度(BI)、绿度(GVI)和湿度(WI)计算对象的算数特征；主要植被指数计算公式如(1)‑(4)：NDVI=NIR-RNIR+R---(1)]]>NDWI=NIR-GNIR+G---(2)]]>RVI=NIRR---(3)]]>DVI＝NIR‑R                 (4)式中，NIR为近红外波段，波长为0.85～0.88μm；R为红色波段，波长为0.64～0.67μm；G为绿色波段，波长为0.53～0.59μm；五、利用步骤四中得到的所有特征信息，确定区分沼泽湿地与非沼泽湿地的特征值，从Landsat‑8 OLI影像中去除非沼泽湿地对象；提取过程和参数如下：1、设定水体指数阈值，将水体剔除，参数设置为：BI<＝0.24；2、设定林地指数阈值，将林地剔除，参数设置为：RVI>＝9.5或者Hue R:G:B＝SWIR1:NIR:R>＝0.25；3、设定建筑用地指数阈值，将建筑用地剔除，参数设置为：GVI<＝0.3；4、设定沼泽湿地指数阈值，将非沼泽湿地剔除，参数设置为：Slope<＝5；六、在步骤五所提取的沼泽湿地对象基础上，利用步骤四所涉及的草丛沼泽、灌丛沼泽和森林沼泽的特征逐步判定并提取草丛沼泽、灌丛沼泽和森林沼泽对象；提取过程和参数如下：1、设定森林沼泽指数阈值，提取森林沼泽，参数设置为：NDVI/R>＝0.0021；2、设定灌丛沼泽指数阈值，提取灌丛沼泽，参数设置为：Standard deviation BI<＝0.0055；3、设定草丛沼泽指数阈值，提取草丛沼泽，参数设置为：Standard deviation BI>＝0.0055；七、导出草丛沼泽、灌丛沼泽和森林沼泽对象，获得沼泽湿地矢量数据；八、在ArcGIS 10.0软件中制作沼泽湿地专题地图。...

【技术特征摘要】

【专利技术属性】
技术研发人员：满卫东，王宗明，刘明月，贾明明，毛德华，任春颖，路春燕，欧阳玲，于皓，
申请(专利权)人：中国科学院东北地理与农业生态研究所，
类型：发明
国别省市：黑龙江;23