传统与深度学习遥感变化监测遥感技术路线与方法.docx

《传统与深度学习遥感变化监测遥感技术路线与方法.docx》由会员分享，可在线阅读，更多相关《传统与深度学习遥感变化监测遥感技术路线与方法.docx（35页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、传统与深度学习遥感变化监测遥感技术路线与方法基于多时相影像进行地表覆盖变化监测是遥感的一个重要应用，在城市违法 建设监测、水政违法建设监测、森林砍伐监测等领域得到广泛应用。由于遥感算 法、数据源质量等原因，遥感变化信息提取一般采用目视解译方式进行，但是目 视解译方式费时费力，大区域工作效率很低。而深度学习可以在很短时间内按照 模型训练要求，快速全面的进行解译，可以大大提高遥感解译的效效率,因此一 般采用传统解译方法与深度学习相结合的方法进行。一、传统多时相影像动态检测技术遥感影像具有周期性特点，利用这一特性就可以用于监测一个地区的地表变 化情况动态监测概述 动态检测的关键技术基于ENVI的动态

2、检测1.1动态监测概述很多人对动态监测和动态检测两个名词有疑惑。动态监测是一个广义的名词, 泛指数据预处理、变化信息发现与提取、变化信息挖掘与应用等，以对整个流程 的表达。动态检测是一个狭义的名词，主要指局部数据预处理、变化信息发现与 提取。在一定的意义和场合上讲，两个词的意思是一样的。当然很多地方把这两 个词说成变化监测和变化检测。遥感动态检测就是从不同时期的遥感数据中，定量地分析和确定地表变化 的特征与过程。它涉及到变化的类型、分布状况与变化量，即需要确定变化前后 的地面类型、界线及变化趋势，能提供地物的空间分布及其变化的定性和定量信 息。目前，遥感变化检测技术大多是针对两个时相的遥感影像

3、进行操作。根据处 理过程来分，遥感变化检测方法可分为三类：(1)图像直接比拟法构筑物构筑物（新修道路）填河动土动土图：变化目标特征2.2深度学习变化监测数据源数据源使用高分辨率卫星数据，现在以高分一号PMS传感器融合校正后的 正射影像DOM产品为例进行说明。研究区如下列图所示：表：数据源分辨卫星数据级别波段 率 时间（米）高分一正射镶嵌影R、G、B f 桐-52 两时期号PMS 像DOM 3波段图：研究区技术流程数据源n预处理深度学习变化地物提提取结果图：技术流程预处理预处理主要包括两时期影像的波段合成。在ENVI工具箱中选择，Raster Management/Build Layer Sta

4、ck,选择两时相3波段影像数据，其他参数按照默认， 设置文件输出路径，点击0K进行波段合成，得到两时期6波段合成结果。Build Layer Stack XInput Rasters：fir regional. dat.regi on_T2. dat-2 rasters.Drag and drop rasters to reorder.Drag and drop rasters to reorder.Grid Definition (opti onel)S S B口 CCoord, system + extents + pixel sizevCoordinate System ExtentsN

5、orthWestEastSouthPixel SizeXYResampling MethodNearest Neighbor Output Raster File J.) Virtual RasterD:Tl_T2_LayerStack. date0 Di splay resultOK | P5anceM图：两时期波段合成深度学习深度学习操作流程主要包括三个局部：首先，创立标签图像。标签栅格的生 成一般先选择典型子区绘制样本后生成标签栅格；第二，训练深度学习模型。有 了标签栅格就可以初始化深度学习模型并进行模型训练，推荐使用随机化参数组 的方法进行模型训练；第三，使用训练好的模型进行目标提取。

6、训练好模型之后 就可以使用训练好的模型进行目标提取。这三步操作是深度学习的主要内容，之 后可进一步将分类结果转换为矢量ROI再次进行编辑生成新的标签栅格重复训 练模型，直到得到较好的结果。创立标签图像影像收据源影像收据源选取影像子区域I fri图：深度学习操作流程创立标签图：1）样本获取创立标签图像之前首先需要获取样本。通过人工目视解译的方法绘制研究区 范围内ROI样本，在ENVI上方工具栏点击 按钮，创立变化图斑ROL沿着变 化图斑轮廓绘制变化图斑样本。如下列图所示：22T31d Vern 3g* Re92c o4 rfteresi flp图：绘制ROI样本如果已有Shapefile矢量样本

7、，也可在ROI工具面板选择FileImportImport Vector导入已有的Shapefile矢量作为ROI样本文件。图：研究区矢量样本关于样本的获取，通常有如下三个原那么：样本区域范围内，全面覆盖多种地物类型。绘制样本尽可能的多。样本轮廓尽可能精确，不要多余地物。2）构建样本库构建样本库可以分为如下四步：样本矢量生成缓冲区将缓冲区合并通过合并后的缓冲区矢量生成最小外接矩形使用最小外接矩形裁剪栅格沿矢量周围扩展1000*1000像素生成缓冲区并将缓冲区合并，之后生成最小 外接矩形（最小外接矩形生成工具可参考244节3）内容），使用外接矩形将影 像切片创立模型训练样本库。矢量裁剪在Tool

8、box中选择Extensions/Frame Subset via Shapefile矢量分幅裁剪工具对影像进行分幅裁剪。图：矩形裁剪创立样本库3）生成标签栅格库标签栅格用于深度学习模型训练，需要通过样本ROI和栅格影像构建标签 栅格。对于大批量标签栅格的生成，可使用ENVI Modeler建模工具进行批量标 签栅格的生成。在ENVI工具箱选择Task Processing/ENVI Modeler翻开建模工 具，构建如下批处理模型，点击Run按钮运行模型，生成批量标签栅格。注：Generate Output Filename 节点可在 App Store 中下载。Build L&bel Ra

9、ster rcm ROIMi ：=Build L&bel Raster rcm ROIMi ：=图：ENVI Modeler批处理生成标签栅格图：标签栅格库模型训I练对于深度学习模型的训练，推荐使用随机化参数组的方法训练模型。在进行 随机化参数训练之前，需要先初始化一个深度学习模型。在Toolbox工具箱选择 Deep Learning/Train TensorFlow Mask Model,点击 Input Model 下方的 New Model 弹出初始化模型参数设置对话框，Number of Bands为波段合成后影像的波段数 量，这里两时期波段合成后为6波段影像。Number of C

10、lasses为变化类别，这里 变化类别设置为1类即为变化地物类，Output Model选择初始化模型输出路径点 击0K,生成初始化模型。Initialize ENVINet5 Multiclass ModelOK Cancel图：初始化模型接下来使用随机化参数训练工具进行深度学习模型训练。在Toolbox中，选 择Deep Learning/Deep Learning Guide Map翻开深度学习向导工具，在工具对话 框选择 Tools Randomize Training Parameters Using the ENVI Modeler 翻开深度 学习模型随机化参数训练工具。直接点击上

11、方工具栏的Run按钮，弹出随机化 参数训练对话框。图像直接比拟法是最为常见的方法，它是对经过配准的两个时相遥感影像 中像元值直接进行运算和变换处理，找出变化的区域。目前常用的光谱数据直 接比拟法包括图像差值法、图像比值法、植被指数比拟法、主成分分析法、光 谱特征变异法、假彩色合成法、波段替换法、变化矢量分析法、波段交叉相关 分析以及混合检测法等。图像差值法图像差值法就是将两个时相的遥感图像相减。其原理是：图像中未发生变 化的地类在两个时相的遥感图像上一般具有相等或相近的灰度值，而当地类发 生变化时，对应位置的灰度值将有较大差异。因此在差值图像上发生地类变化 区域的灰度值会与背景值有较大差异，从

12、而使变化信息从背景影像中显现出来。 光谱特征变异法同一地物反映在一时相影像上的信息与其反映在另外时相影像上的光谱 信息是一一对应的。当将不同时相的影像进行融合时，如同一地物在两者上的 信息表现不一致时，那么融合后的影像中此地物的光谱就表现得与正常地物的 光谱有所差异，此时称地物发生了光谱特征变异，我们就可以根据发生变异的 光谱特征确定变化信息。 假彩色合成法由于地表的变化，相同传感器对同一地点所获取的不同时相的影像在灰度 上有较大的区别。在进行变化信息的发现时，将前、后两时相的数据精确配准， 再利用假彩色合成的方法，将后一时相的一个波段数据赋予红色通道，前一时 相的同一波段赋予蓝色和绿色通道。

13、利用三原色原理，形成假彩色影像。其中， 地表未发生变化的区域，合成后影像灰度值接近，而土地利用发生变化的区域 那么呈现出红色，即判定为变化区域。 波段替换法在RGB假彩色合成中，G和B分量用前时相的两个波段，用后一时相的 一个波段影像组成R分量，在合成的RGB假彩色图像上能够很容易地发现红 色区域即为变化区域。(2)分类后比拟法分类后结果比拟法是将经过配准的两个时相遥感影像分别进行分类，然后比拟晶 Deep_Learning_Randomize_TrainingIterationsIterations16Input Model,ini t_change_detecti on_model. h5

14、Name： ENVI Deep LearningBands： 6 Classes： 1 Full info.Training Rasters?k_dat_Area274_1atel_raster, dat 4 ?k_dat_Area275_1atel_raster, dat 改 dat Area276 label raster, dat da+ At a 77 1 ah al yay + ay da +278 files.Validation Rasters?k_dat_Ar e a2 74_1ab el_r as t er. dat 人 ;k_dat_Area275_1atel_raster

15、, dat :k dat Area276 label raster, dat *da+ Ay a 7 7 1 ah al yas + ay da +278 files.Augment ScaleOYes NoOVes N。OVes N。Augment Rotati onOKCancel图：随机化参数模型训练Iterations：为迭代训练的次数，这里默认16次，表示一共训练得到16组 模型。Input Model：选择初始化模型或者选择需要再次训练的模型。Training Rasters：选择上一步生成的标签栅格。Validation Rasters：同样选择上一步生成的标签栅格，工具会自动选

16、择80% 数据用于训练，20%数据作为验证。也可以将标签栅格库分为两局部，一局部 用于训练一局部用于验证。Augment Scale：是否使用缩放增强。增强是深度学习常用的一种技术， ENVI深度学习使用缩放和旋转来补充原始训练数据。它可以减少你所需的标 记数量，特别是捕捉各种方向不同和大小不一的同种要素，增强还可以提高模 型的能力，将它所学到的知识推广到新的图像中。但在样本充足的前提下使用 增强会使模型的泛化能力过高，影响提取精度。这里选择No,不进行缩放增 强。Augment Rotation：选No,不进行旋转增强。Raster to Classify：选择要进行变化地物提取的栅格影像，

17、这里选择波段 合成后的研究区影像。Output Directory：选择提取结果存放路径。2.243模型分类通过查看随机化参数训练工具的提取效果，选择泛化能力最强的模型作为深 度学习变化监测提取模型。如下列图显示了不同泛化能力模型的提取效果，图中泛 化能力强的模型不仅没有“漏提”，还发现了人工没有发现的变化图斑，但泛化能 力强的模型也有“过提”的缺点，但相比“漏提”，我们更能接受“过提”。图：泛化能力弱图：泛化能力强选择好模型之后就可以使用已经训练好的模型进行变化地物提取。在 Toolbox 中,选择 Deep Learning/TensorFlow Mask Classification 翻

18、开深度学习模 型分类工具，Input Raster选择研究区数据，Input Trained Model选择训练好的最 优模型，Output Classification Raster选择提取结果输出路径和文件名，Output Class Activation Raster选择类激活栅格输出路径，点击OK。魅 TensorFlow Mask ClassificationInput RasterInput Trained ModelOutput Classification Raster (op ti on b1)Output Class Activation Raster (optional)

19、曲研究区.datiined_change_detecti on_model. h5Name： ENVI Deep LearningBands： 6 Classes： 1 Full info.0 Di splay resultD：CAM. dat2 Di splay resultOK Cancel图：使用已训练好的模型进行分类得到深度学习提取的变化图斑结果，如下列图所示图中红色图斑表示发生变化 的区域，黑色为背景区域。图：深度学习变化监测结果生成矢量结果1)分类后处理分类后处理主要包括小图斑处理(Classification Aggregation)、聚类处理(Clump)、 主次要分析(Maj

20、ority/Minority Analysis)。通过分类后处理可对一些细小图斑进 行过滤。小图斑处理在 Toolbox 中，选择 Classification/Post Classification/Classification Aggregation, 在弹出的参数对话框中Input Raster选择分类后结果，Minimum Size和Aggregate Unclassified Pixels按照默认设置，设置输出路径点击0K,得到小图斑处理结果。聚类处理在 Toolbox 中，选择 Classification/Post Classification/Clump Classes,选择上

21、一步 的输出结果，在弹出的对话框中，Class Order删除Unclassified这一类，其他参 数按照默认，设置输出路径点击OK,得到聚类处理结果。主次要分析在 Toolbox 中，选择 Classification/Post Classification/Majority/Minority Analysis, 选择上一步的输出文件，在弹出的参数设置对话框中选择变化类别，Analysis Method默认Majority进行主要分析，Kernel Size设置为27x27,其他参数按照 默认，设置结果输出路径点击0K,得到主要分析处理结果。2）分类栅格转矢量在 Toolbox 中, 选择

22、 Classification/Post Classification/Classification to Vector 将分类栅格转为矢量。选择上一步的输出结果，在弹出的对话框中，Export Classes 选择变化类别，其他参数按照默认，设置输出路径点击0K,得到矢量结果。Convert Classification to Vecto.Input RasterOes NoApply Symbologyo 2)best model class, datOutput VectorOKOKCancel图：分类栅格转矢量3）生成最小外接矩形通过输出的矢量结果生成最小外接矩形，在Toolbox中，

23、选择Vector/Vector to Bounding Box,在弹出的对话框中选择上一步输出的矢量，Oriented Bounding Box选择Yes,根据图斑边界生成最小外接矩形，设置文件输出路径点击0K, 得到最终结果。画 Vector Records to Bounding BoxInput Vector圾 best_model_class. shpUnit： MetersOriented Bounding BoxMinimum Value (optional)Maximum Value (optional/Yes I NoD:best_model_class_box. shpMi

24、rDFrevi ew1 Di splay resultOKOutput Vector图：生成最小外接矩形2.3变化监测结果结果统计深度学习变化监测在中完成，计算机环 境：Dell Precision T3640 图形工作站（CPU： Intel i9-9900K, GPU： NVIDIA RTX2080ti显存11GB,内存：64GB,硬盘：SSD固态），模型训练时间：35 分钟，目标提取时间：5小时30分钟，22135平方公里范围内共发现变化图斑 42255 个。图：变化监测提取结果使用训练样本对深度学习训练模型进行精度验证，结果如下：Loss （损失）：0.111Accuracy （正确率

25、）：0.946Precision （用户精度）：0.783Recall （生产者精度）：0.953Fl （F1 分数）：0.8532.3.1 监测结果区域内人工目视绘制变化图斑共计653个，深度学习模型全部识别发现，无 “漏提”现象。同时深度学习还发现了大量人工没有发下的遗漏变化图斑。监测结 果如下列图所示，图中红色矢量为人工目视解译绘制的变化图斑，黄色矩形框为深 度学习发现的变化区域。I 二分类结果得到变化检测信息。虽然该方法的精度依赖于分别分类时的精度和分类 标准的一致性，但在实际应用中仍然非常有效，该方法的核心是基于分类基础上 发现变化信息。（3）直接分类法结合了图像直接比拟法和分类后结

26、果比拟法的思想，常见的方法有：多时相主 成分分析后分类法、多时相组合后分类法等。多时相主成分分析后分类法当地物属性发生变化时，必将导致其在影像某几个波段上的值发生变化，所 以只要找出两时相影像中对应波段值的差异并确定这些差异的范围，便可发现变 化信息。在具体工作中将两时相的影像各波段组合成一个两倍于原影像波段数的 新影像，并对该影像作PC变换。由于变换结果中前几个分量上集中了两个影像 的主要信息，而后几个分量那么反映出了两影像的差异信息，因此可以试着抽取后 几个分量进行波段组合来发现变化信息。以上变换检测方法大多只是变化信息的发现，还有一个过程是变化信息的提 取，由以下方法供选择：手工数字化法

27、屏幕数字化区域生长法图像自动分类 监督分类 非监督分类 面向对象的特征提取法图像分割 手工阈值分割自动阈值分割图：深度学习变化监测结果组合法上述检测方法和信息提取方法不能说哪个绝对的好与坏，只能是根据不同的数 据源和不同的应用需求选用适合的方法。我们针对工程的实际情况会选择最合适 的方法进行检测。1.2 动态检测关键技术数据源的选择由于动态检测需要多时相数据，数据源选择时候，除了考虑检测范围、检测精 度外，重点还需要考虑物候和影像成像条件。物候和成像条件都会给变化信息带 来很大的噪声，当然了，经过一些图像处理过程可以一定程度上减少这方面的噪 声。1.2.1 数据预处理动态检测过程中，非常重要的

28、两个预处理过程是影像配准和大气校正影像配准影像的配准可以有单个文件的精确几何校正来保证，也可以有一个文件作为基 准来配置另一个时相的文件几何精校正是在影像和地形图上找相同的点，找出来的点是具有相同地理坐 标的地面控制点.影像有畸变，地形图（认为）它是准确的，因此控制点多一些，能更 好地建立畸变模型,从而对整幅影像进行校正。相对大气校正在多时相遥感图像中，除了地物的变化会引起图像中辐射值的变化外，不变的 地物在不同时相图像中的辐射值也会有差异O辐射校正是消除非地物变化所造成 的图像辐射值改变的有效方法，按照校正后的结果可以分为2种，绝对辐射校正 方法和相对辐射校正方法。绝对辐射校正方法是将遥感图

29、像的DN（DigitalNumber）值转换为真实地表反射率的方法，它需要获取影像过境时的地表测量数据，并考虑地形起伏等因素来校正大气和传感器的影响，因此这类方法一般都很复杂。相对辐射校正是将一图像作为参考（或基准）图像，调整另一图像的DN使得两时相影像上同名的地物具有相同的DN值,这个过程也叫多时相遥感图像的光谱归一化。在动态检测中，相对大气校正用的还是比拟多。1.2.2 变化检测方法选择地表变化信息可分为两种，一是转化(Conversion),另一是改变(Modification), 前者是土地从一种土地覆盖类型向另一种类型的转化，如草地转变为农田、森林 转变为牧场，后者是一种土地覆盖类型

30、内部条件(结构和功能)的变化，如森林由 密变疏或由一种树种组成变成另外一种组成的改变、植物群落生物量、生产力、 物候现象变化。所以我们在选用变换检测方法时候应考虑这点，比方做植被生物 量的监测，就属于改变，如果选用分类后处理的方法就不管用了。反过来基于图像直接比拟法，不能够直接确定变化区域是属于何种变化类型, 需要实地调查与图像比照等一系列步骤才能确定变化类型，这对于土地利用变化 监测时候，那么需要慎重考虑用此类方法。1.2.3 变化信息提取方法的选择有些变化检测方法得到的变化信息是单波段，如图像直接比拟法，监督分类法作用就不是很大，可以考虑用影像分割和面向对象特征提取方法。而 且图像直接比拟

31、法在确定是否变化时，往往需要确定一定的阈值，而阈值确实定 那么会直接影响变化信息提取的准确性，阈值确实定常见的如直方图法，样本挖掘 法(C4.5算法)等。1.3 基于ENVI的动态检测ENVI集成了局部动态检测方法，包括图像直接比拟法、分类后比拟法、PCA 变换、Two-Color Multiview、MNF变换法、ICA变换法、波谱角检测方法。当 然还可以自定义方法。1.3.1 图像直接比拟法ENVI中的图像直接比拟法就是对两时相影像做差值或者比值运算，整合了 一些预处理功能，如数据值归一化和单位的统一。下面以经过配准的两时相影像 为例介绍这个功能的使用。(1)单击主菜单-File-Open

32、 Image File,将两时相影像同时翻开，并用 Display Link功能杳看变化区域。(2)单击主菜单-Basic Tools-Change Detection-Compute Difference Map, 分别选择前一时相影像一个波段和后一时相影像的一个波段。(3)在 Compute Difference Map Input Parameters 面板中，可以选择计算方 法(差值或者比值)、归一化(0-1)和单位统一，设置变化等级以及设置变化 等级划分阈值。单击Define Class Thresholds按钮，可以对每一个变化范围进行 划分，如图1所示。选择一个路径输出。(4)结

33、果查看和统计：在Display中将结果显示，用Display Link链接功能 将前后时相和结果链接查看。单击Image-Tools-Color M叩ping-Class Color Mapping,查看各个变化等级的颜色和代码。单击主菜单-Classification-Post Classification-Class Statistics,统计各个变化。1.3.2 分类后比拟法ENVI中的分类后比拟法是通过比拟两时相影像分类结果，获得变化类型、面 积、百分比等。下面介绍这个工具的使用。(1)单击 主菜单-File-Open Image File,将两时相的分类图翻开。(2)单击 主菜单-B

34、asic Tools-Change Detection-Change Detection Statistics,选择前后时相的分类图。(3)在Define Equivalent Class面板中，如果两个时相的分类图命名规那么一 直，那么会自动将两时相上的类别关联；否那么需要在Initial State Class和Final State Class列表中手动选择相对应的类别，如图3所示，点击0k按 钮。(4)在结果输出面板中，选择统计类型：像素(Pixels) 百分比(Percent) 和面积(Area),选择路径输出结果。(5)结果以二维表格和图像现实展现。1.3.3 流程化图像处理工具动

35、态监测ENVI的流程图像处理工具中集成了动态监测功能，它集成了影像配准和变化信 息发现过程等，下面介绍这个工具的使用。(1)单击主菜单-Spectral-SPEAR Tools-Change Detection,翻开流程化 图像处理工具中的动态监测功能，分别将两时相影像文件输入，单击Next按钮。(2)如果两个时相影像未精确配准，这一步可以通过手动或者自动选择同 名点的方式对影像进行配准，这里选择第三项，单击Next按钮。(3)这一步是目视查看两时相影像的变化情况，点击Next按钮到下一步。(4)这一步是选择变化检测方法，提供了四大类方法，包括Two-ColorMultiview图像变换（PC

36、A、MNF和ICA变换）、图像直接比拟法和波谱角 检测方法。每一种方法都有对应的高级设置，比方选择图像直接比拟法可以选择 基于黑暗象元的相对大气校正和光谱归一化。这几类方法可以多项选择，这里我们全 部选择，单击Next按钮。（5）这一步是查看变化信息监测结果，选择使用不同检测方法的结果查看, 点击Finish按钮完成。（6）对检测到的变化信息，可以通过图像分割、图像分类、决策树分类或 者面向对象特征提取等方法提取变化信息。二、利用ENVI深度学习进行遥感变化监测技术流程工程影像在 ENVI5.6+Deep Learning 1.1.3 中完成。ENVI Deep Learning 1.1.3

37、适配ENVI 5.6,要求计算机必须具备NVIDIA显卡，对显卡环境要求如下：显卡驱动版本要求36.06或更高版本。要求显卡CUDA计算能力（Compute Capability）在58.6之间。推荐8G以上显存的显卡，如性价比拟高的NVIDIA GeoForce系列显卡 GeoForce RTX 3080、3090 等。2.1变化地物类型在使用深度学习方法提取变化目标之前，首先要明确变化地物类型，变化地 物类别划分为如下几类：建筑物（主要为变化厂房），构筑物（新修道路等）， 填河，动土，堆土等。其影像特征如下列图所示：建筑物构筑物动土填河图：变化目标如下列图所示，左侧为第一时相影像右侧为第二时相影像，通过两时期比照, 展示更多变化目标特征。建筑物建筑物（厂房）