一种基于稀疏表示分类的遥感图像变化检测方法
【专利摘要】本发明涉及一种基于稀疏表示分类的遥感图像变化检测方法，包括以下步骤：1)输入经过校正处理的变化前和变化后的遥感图像；2)获取变化前后图像的差值图像；3)在差值图像中对变化区域选取少量样本点；4)对每个像元计算利用样本点对其进行稀疏表示得到的估计值；5)计算各样本点实际值与估计值的差值；6)若差值小于给定阀值，则该样本点属于变化区，否则，属于无变化区。本发明的有益效果为：提出了一种判断遥感图像中有无发生变化的快速检测方法，仅需少量样本点，即可得到准确的变化检测区域；方法仅考虑变化区域特性，对复杂背景有强的抗干扰能力。
【专利说明】
-种基于稀疏表示分类的遥感图像变化检测方法
技术领域
[0001] 本发明设及遥感图像处理领域，尤其设及一种多时相遥感图像的变化检测方法。
【背景技术】
[0002] 利用遥感图像对地表的变化进行识别和范围检测，对于灾害范围统计、损失评估 等具有重要的实际意义。由于不同的地物在遥感图像中呈现出不同的光谱特征，可W通过 多时相遥感图像中的差异对发生变化的区域进行识别和检测。在进行变化检测时，大多数 情况下，如对火灾、洪灾、滑坡等进行检测时，最重要的是对变化区域进行快速的识别，对于 未发生变化的背景部分并不关注。因而，在此提出了一种针对变化区域进行快速高精度识 别的方法。
【发明内容】
[0003] 本发明的目的是提出了一种判断遥感图像中有无发生变化的快速检测方法。本方 法基于稀疏表示理论，对遥感图像中的变化区域进行半自动检测。在变化检测中，着眼于对 发生变化的区域进行分析和精细识别，仅需要对变化区域选取少量样本，而无需考虑发生 变化地物所处的复杂的、多样化的背景，具有强的抗干扰能力。
[0004] 为达到上述目的，本发明的技术方案提供一种基于稀疏表示分类的变化检测方 法，所述方法包括W下步骤：
[0005] 1)输入经过校正处理的变化前和变化后的遥感图像Xi和X2。前后两个时相的图像 必须经过几何校正和福射校正，保证前后时相图像的几何误差不超过一个像元，且地物福 射特性较为一致；
[0006] 2)获取变化前后图像的差值图像Xt, Xt = Xi-X2,并对差值做归一化处理到0-1值域 范围；
[0007] 3)在差值图像中，对变化区域选取n个样本点，获取其图像坐标值[i，j]，并得到差 值图像中该坐标的数值Xt[i，j]作为训练样本；
[000引4)对每个像元y，利用Xt对其进行稀疏线性表示得到估计值;;=JTiS:, W估计结果与
原值相似程度和系数的稀疏情况为约束，对W上稀疏表示的系数a进行求解，其目标方程 为：
[0009]
[0010] 其中，A是正则系数，用于平衡估计值与原始值的相似性和系数的稀疏程度。用 Matlab中的SPAMS工具箱对稀疏表示系数进行求解；
[0011] 5)计算各样本点实际值与估计值的残差值r(y):
[0012]
[0013] 6)对残差进行误差曲线分析确定阀值。对残差值
扣曽量，取 得10000个阀值t，对每个阀值计算检测率和误检率，并做出误差曲线。检测率(纵轴)为实际 检测出的变化占总变化的比率，误检率(横轴)为检测出的错误变化数占无变化区的比率。 误差曲线趋势类似于上升的对数函数曲线，误检率为加寸检测率为0,误检率为1时检测率趋 近于1。随误检率增加检测率先大幅增加，到一定程度时斜率降低。取误差曲线拐点处的检 测率和误检率所对应的为t值为阀值，即误差曲线由睹变平处的值，取该值作为阀值能在保 证较高检测率的情况下获得最低的误检率。
[0014] 7)若差值r(y)小于给定阀值，则该样本点属于变化区，否则，属于无变化区。
[0015] 本发明的有益效果为：
[0016] 1)本发明仅需对变化区域选取少量样本即可得到准确的变化检测区域，无需获取 变化区域之外的背景地物样本，对方法的训练较为简单。
[0017] 2)本发明仅对变化区域的光谱特性进行分析和描述，并寻找类别界限，无需考虑 背景区域的特性，对复杂背景有强的抗干扰能力。
[0018] 3)本发明采用稀疏表示分类技术对变化区域进行描述，发挥稀疏表示技术利用少 数样本点就可W进行完善信息描述的优势，将发生变化的区域作为一个类别进行精确识 别。
【附图说明】
[0019] 图1是本发明实施例所述的基于稀疏表示分类的遥感图像变化检测方法的流程示 意图。
[0020] 图2是误差曲线图。
【具体实施方式】
[0021 ]下面根据附图对本发明作进一步详细说明。
[0022] 如图1所示，本发明实施例所述的一种遥感图像变化检测方法，所述方法包括W下 步骤：
[0023] 1)输入经过校正处理的变化前和变化后的遥感图像Xi和X2。前后两个时相的图像 必须经过几何校正和福射校正，保证前后时相图像的几何误差不超过一个像元，且地物福 射特性较为一致；
[0024] 2)获取变化前后图像的差值图像Xi2 = Xi-X2,并对差值Xi2做归一化处理到0-1值域 范围；
[0025]
[0026] 3)在差值图像中，对变化区域选取不少于50个样本点，获取其图像坐标值[i，j]， 并得到差值图像中该坐标的数值Xt[i，j]作为训练样本；
[0027] 4)对每个像元y，利用稀疏表示的方式计算利用变化区域训练样本对其进行稀疏 表示得到的估计值y = Xta，求解的目标方程为：
[002引
[0029]其中，A是正则系数，用于平衡估计值与原始值的相似性，和系数的稀疏程度，此处 取入=I。M 112为取平方和，I I I I I为取绝对值之和。
[0030] 用Matlab中的SPAMS工具箱对稀疏表示系数进行求解。
[0031] 5)计算各样本点实际值与估计值的残差值：
[0032]
[0033] 6)对残差进行误差曲线分析确定阀值。对残差值W
J增量，取 得10000个阀值t，对每个阀值计算检测率和误检率，并做出误差曲线。检测率(纵轴)为实际 检测出的变化占总变化的比率，误检率(横轴)为检测出的错误变化数占无变化区的比率。 误差曲线趋势类似于上升的对数函数曲线，误检率为加寸检测率为0,误检率为1时检测率趋 近于1。随误检率增加检测率先大幅增加，到一定程度时斜率降低(实施本方法得到的误差 曲线案例如图2所示）。取误差曲线拐点处的检测率和误检率所对应的为t值为阀值，即误差 曲线由睹变平处的值，取该值作为阀值能在保证较高检测率的情况下获得最低的误检率。
[0034] 7)若差值r(y)小于给定阀值，则该样本点属于变化区，否则，属于无变化区。
[0035] 本发明所述的遥感图像变化检测方法具有W下特点：
[0036] 1、采用单类的稀疏表示分类技术，只需对变化区域取少量样本即可进行变化检 测，步骤3)中，仅需对变化区域进行样本点选取。
[0037] 2、基于稀疏表示技术得到各像元与变化区的相似性。在所述步骤4)中，对每个像 元得到其用变化区样本表示的估计值，该值可用于判定该像元是否属于变化区。若该像元 属于变化区，则该像元值能够通过变化区的样本进行高精度的稀疏表示;若不属于变化区， 则估计值与实际值会有较大的残差。
[0038] 3、通过误差分析的方式得到变化检测的阀值:在步骤6)中，对残差进行阀值探测， 对每个待选的阀值计算误差曲线，并通过误差曲线的趋势分析得到最优阀值，W实现高精 度的变化检测。
[0039] 由W上实施例可W看出，本发明实施例基于稀疏表示的技术特点，通过对单类别 的特性进行分析，仅需对变化区域选取少量样本即可进行变化区域的精确探测，需要输入 少，且能获得很好的检测效果。
[0040] 本发明通过上面的实施例进行举例说明，但是，本发明并不限于运里所描述的特 殊实例和实施方案。任何本领域中的技术人员很容易在不脱离分发明精神和范围的情况下 进行进一步的改进和完善，因此分发明只受到本发明权利要求的内容和范围的限制，其意 图涵盖所有包括在由附录权利要求所限定的本发明精神和范围内的备选方案和等同方案。
【主权项】
1. 一种基于稀疏表示分类的遥感图像变化检测方法，其特征在于：该方法包括以下步 骤： 1) 输入经过校正处理的变化前和变化后的遥感图像XdPX2;前后两个时相的图像必须 经过几何校正和辐射校正，保证前后时相图像的几何误差不超过一个像元，且地物辐射特 性较为一致； 2) 获取变化前后图像的差值图像XhXtiXi-Xs，并对差值做归一化处理到0-1值域范围； 3) 在差值图像中，对变化区域选取η个样本点，获取其图像坐标值[i，j]，并得到差值图 像中该坐标的数值Xt[i，j]作为训练样本； 4) 对每个像元y，利用Xt对其进行稀疏线性表示得到估计值;^ ，.以估计结果与原值 相似程度和系数的稀疏情况为约束，对以上稀疏表示的系数α进行求解，其目标方程为：其中，λ是正则系数，用于平衡估计值与原始值的相似性和系数的稀疏程度；用Matlab 中的SPAMS工具箱对稀疏表示系数进行求解； 5) 计算各样本点实际值与估计值的残差值r(y):6) 对残差进行误差曲线分析确定阀值；对残差值以为增量，取得 10000个阀值t，对每个阀值计算检测率和误检率，并做出误差曲线;检测率(纵轴)为实际检 测出的变化占总变化的比率，误检率(横轴)为检测出的错误变化数占无变化区的比率;误 差曲线趋势类似于上升的对数函数曲线，误检率为〇时检测率为〇,误检率为1时检测率趋近 于1;随误检率增加检测率先大幅增加，到一定程度时斜率降低;取误差曲线拐点处的检测 率和误检率所对应的为t值为阀值，即误差曲线由陡变平处的值，取该值作为阀值能在保证 较高检测率的情况下获得最低的误检率； 7) 若差值r(y)小于给定阀值，则该样本点属于变化区，否则，属于无变化区。
【文档编号】G06T7/00GK105957049SQ201610077917
【公开日】2016年9月21日
【申请日】2016年2月3日
【发明人】冉琼, 张蒙蒙, 李伟
【申请人】北京化工大学