首页 / 专利库 / 人工智能 / 深层模型 / 一种基于深层神经网络的遥感图像覆被变化检测方法

一种基于深层神经网络的遥感图像覆被变化检测方法

阅读:779发布:2020-05-12

专利汇可以提供一种基于深层神经网络的遥感图像覆被变化检测方法专利检索,专利查询,专利分析的服务。并且本 发明 公开了一种基于深层神经网络的遥感图像覆被变化检测方法,先建立深层神经网络模型的训练集,再将训练集中每个训练样本前M张图像作为输入,通过深层神经网络模型预测输出最后一张图像,并与 图像序列 中的第M+1张图像进行对比,获得的相似度差异添加演化场平滑惩罚项作为损失函数,使用 梯度下降法 优化损失函数,遍历所有训练样本多轮,直到获得满意的损失函数值,完成网络训练,获得训练好的模型,最后将待预测地区最新的时间间隔相近的M张遥感图像输入训练好的模型中,得到的输出图像即为预测的地表覆被图像。,下面是一种基于深层神经网络的遥感图像覆被变化检测方法专利的具体信息内容。

1.一种基于深层神经网络的遥感图像覆被变化检测方法,其特征在于,包括以下步骤:
(1)、建立深层网络模型的训练样本集
(1.1)、图像采集:获取待预测地区近期的时间间隔相等的N张同类型遥感图像,按时间顺序编号为第1~N张;
(1.2)、图像预处理:对每张遥感图像进行几何校正、缩放、剪切和增强处理;
(1.3)、构建训练样本集:将预处理后的N张遥感图像分成N-M个训练样本,N远大于M;
其中,第1个训练样本由第1~(M+1)张遥感图像构成,第2个训练样本由第2~(M+2)张遥感图像构成,然后以此类推,直到第N-M个训练样本,由第(N-M)~N张遥感图像构成;
(2)、构建深层网络模型
卷积神经网络CNN和空间变换网络STN级联组成深层网络模型,深层网络模型的输入为某一个训练样本中的前M张图像,输出为预测的遥感图像;
(3)、训练深层网络模型
(3.1)、深层网络模型初始化,设定阈值θ以及学习速率η;
(3.2)、随机选取一个训练样本,将训练样本中的前M图像输入至深层网络模型中;
在CNN中,先利用编码器的第第一层提取输入图像的特征,再从第二层开始依次对图像特征进行降维处理,得到降维后的特征图;然后利用解码器对特征图进行升维处理,最终输出覆被演化场Φ;
在STN中,基于演化场Φ,利用STN对输入的第M张图像Pm转化成预测输出图像Pn;
(3.3)、计算预测输出图像Pn与该样本中第M+1张图像Pm+1之间的像素误差平方和,记为Σ;计算演化场Φ的梯度值ΔΦ;再将Σ与ΔΦ的和作为该训练样本的损失函数值;
(3.4)、重复步骤(3.2)~(3.3),直至计算出所有训练样本的损失函数值,将所有训练样本损失函数值相加,作为深层网络模型的损失函数值,记为Ω;
(3.5)、比较Ω与θ的大小,如果Ω<θ,则终止训练,得到训练好的深层网络模型,否则,进入步骤(3.6);
(3.6)、利用批量梯度下降算法在损失函数负梯度方向,根据负梯度方向以学习速率η更新深层神经网络模型中的权重参数,然后返回步骤(3.2),再次利用训练样本集进行下一轮训练;
(4)、地表覆被预测
采集M张待预测地区的遥感图像并按时间顺序编号,然后将编号后的遥感图像序列输入训练好的深层网络模型,其输出图像为该预测地区未来的地表覆被预测图像。
2.根据权利要求1所述的一种基于深层神经网络的遥感图像覆被变化检测方法,其特征在于,所述深层网络模型的具体结构为:STN级联在CNN之后;其中,CNN采用UNet的结构,包括编码器和解码器,级数均为N;编码器第n层卷积层提取的图像特征Cn与解码器第N-n层反卷积层的输出Rn同尺寸,Cn与Rn进行通道拼接后,输入至第N-n+1层反卷积层,输出得到升维的特征图Rn-1,然后以此类推,直到第N层反卷积层,最后以C1与R1拼接后的多通道特征图为输入,经反卷积得到遥感图像的演化场Φ;
在STN中,基于演化场Φ,将输入的第M张图像Pm转化成预测输出图像Pn,具体转换的公式为:
Pn(p)=Pm·Φ(p)=∑q∈Z(p')Pm(q)Πd∈{x,y}(1-|p'd-qd|)
其中,·表示空间变换方法,基于演化场Φ通过线性插值将图像Pm变换成输出预测图像Pn;p为Pn中的任意一个像素点的像素坐标,p'=p+Φ(p)为根据演化场Φ计算出的p在Pm中所对应的亚像素点坐标;Z(p')为亚像素点p'对应的4个邻域像素点坐标的集合,d表示图像的像素点坐标维度,p'd、qd分别表示p'、q在维度d下的坐标值。

说明书全文

一种基于深层神经网络的遥感图像覆被变化检测方法

技术领域

[0001] 本发明属于遥感图像变化检测技术领域,更为具体地讲,涉及一种基于深层神经网络的遥感图像覆被变化检测方法。

背景技术

[0002] 随着遥感技术的发展和观测数据的积累,变化检测技术已成为了遥感图像处理与应用的一个研究热点。遥感图像变化检测技术通常通过分析比较同一区域、不同时期获取的两幅或多幅遥感图像,根据图像之间的差异来获取所感兴趣地物、场景或目标的变化情况,已在环境监测、军事侦察和灾情监测等方面得到了广泛的应用。
[0003] 现有的遥感变化检测方法通常分为两类:1、分类后比较法:首先对两幅原图像分别进行独立分类,再对两幅分类好的图像进行逐像素地比较,进而得到最终的变化检测结果;2、比较后分类法,首先通过差值法、比值法或者对数比值法等得到两幅图像的差异图,再对该差异图使用变换、概率分布等方法进行分析,得到最终的变化检测结果。这两种方法都只能从时序影像中获取一些统计特征,例如未来覆被的面积大小。而无法提供变化的位置信息,也无法将变化检测的结果可视化

发明内容

[0004] 本发明的目的在于克服现有技术的不足,提供一种基于深层神经网络的遥感图像覆被变化检测方法,通过预测遥感图像来直观的显示该地区未来全方位的地表覆被情况。
[0005] 为实现上述发明目的,本发明一种基于深层神经网络的遥感图像覆被变化检测方法,其特征在于,包括以下步骤:
[0006] (1)、建立深层网络模型的训练样本集
[0007] (1.1)、图像采集:获取待预测地区近期的时间间隔相等的N张同类型遥感图像,按时间顺序编号为第1~N张;
[0008] (1.2)、图像预处理:对每张遥感图像进行几何校正、缩放、剪切和增强处理;
[0009] (1.3)、构建训练样本集:将预处理后的N张遥感图像分成N-M个训练样本,N远大于M;
[0010] 其中,第1个训练样本由第1~(M+1)张遥感图像构成,第2个训练样本由第2~(M+2)张遥感图像构成,然后以此类推,直到第N-M个训练样本,由第(N-M)~N张遥感图像构成;
[0011] (2)、构建深层网络模型
[0012] 将卷积神经网络CNN和空间变换网络STN级联组成深层网络模型,深层网络模型的输入为某一个训练样本中的前M张图像,输出为预测的遥感图像;
[0013] (3)、训练深层网络模型
[0014] (3.1)、深层网络模型初始化,设定阈值θ以及学习速率η;
[0015] (3.2)、随机选取一个训练样本,将训练样本中的前M图像输入至深层网络模型中;
[0016] 在CNN中,先利用编码器的第一层提取输入图像的特征,再从第二层开始依次对图像特征进行降维处理,得到降维后的特征图;然后利用解码器对特征图进行升维处理,最终输出覆被演化场Φ;
[0017] 在STN中,基于演化场Φ,利用STN对输入的第M张图像Pm转化成预测输出图像Pn;
[0018] (3.3)、计算预测输出图像Pn与该样本中第M+1张图像Pm+1之间的像素误差平方和,记为Σ;计算演化场Φ的梯度值ΔΦ;再将Σ与ΔΦ的和作为该训练样本的损失函数值;
[0019] (3.4)、重复步骤(3.2)~(3.3),直至计算出所有训练样本的损失函数值,将所有训练样本损失函数值相加,作为深层网络模型的损失函数值,记为Ω;
[0020] (3.5)、比较Ω与θ的大小,如果Ω<θ,则终止训练,得到训练好的深层网络模型,否则,进入步骤(3.6);
[0021] (3.6)、利用批量梯度下降算法在损失函数负梯度方向,根据负梯度方向以学习速率η更新深层神经网络模型中的权重参数,然后返回步骤(3.2),再次利用训练样本集进行下一轮训练;
[0022] (4)、地表覆被预测
[0023] 采集M张待预测地区的遥感图像并按时间顺序编号,然后将编号后的遥感图像序列输入训练好的深层网络模型,其输出图像为该预测地区未来的地表覆被预测图像。
[0024] 本发明的发明目的是这样实现的:
[0025] 本发明一种基于深层神经网络的遥感图像覆被变化检测方法,先建立深层神经网络模型的训练集,再将训练集中每个训练样本前M张图像作为输入,通过深层神经网络模型预测输出最后一张图像,并与图像序列中的第M+1张图像进行对比,获得的相似度差异添加演化场平滑惩罚项作为损失函数,使用梯度下降法优化损失函数,遍历所有训练样本多轮,直到获得满意的损失函数值,完成网络训练,获得训练好的模型,最后将待预测地区最新的时间间隔相近的M张遥感图像输入训练好的模型中,得到的输出图像即为预测的地表覆被图像。
[0026] 本发明一种基于深层神经网络的遥感图像覆被变化检测方法还具有以下有益效果:
[0027] 1)、本发明可以直接快速获取预测的地表覆被图像,提供变化的位置信息以及统计特征,可用于科研或管理决策。
[0028] 2)、模型计算简单,无需过多的人工交互,实现了端到端的计算。
[0029] 3)、该模型无需人工标注数据,可以根据已有的无标注的历史遥感影像数据获取预测结果。附图说明
[0030] 图1是本发明一种基于深层神经网络的遥感图像覆被变化检测方法流程图
[0031] 图2是CNN模的结构图;
[0032] 图3是训练深层网络模型流程图。

具体实施方式

[0033] 下面结合附图对本发明的具体实施方式进行描述,以便本领域的技术人员更好地理解本发明。需要特别提醒注意的是,在以下的描述中,当已知功能和设计的详细描述也许会淡化本发明的主要内容时,这些描述在这里将被忽略。
[0034] 实施例
[0035] 图1是本发明一种基于深层神经网络的遥感图像覆被变化检测方法流程图。
[0036] 在本实施例中,如图1所示,本发明一种基于深层神经网络的遥感图像覆被变化检测方法,包括以下步骤:
[0037] S1、建立深层网络模型的训练样本集
[0038] S1.1、图像采集:获取待预测地区近期的时间间隔相等的N张同类型遥感图像,按时间顺序编号为第1~N张;
[0039] S1.2、图像预处理:对每张遥感图像进行几何校正、缩放、剪切和增强处理;本实例中,使用ENVI软件对所有遥感图像进行几何校正,消除传感器本身、大气、地形、地球自转引起的遥感图像畸变;使用灰度级校正增强;最后统一所有图像尺寸
[0040] S1.3、构建训练样本集:将预处理后的N张遥感图像分成N-M个训练样本,N远大于M;
[0041] 其中,第1个训练样本由第1~(M+1)张遥感图像构成,第2个训练样本由第2~(M+2)张遥感图像构成,然后以此类推,直到第N-M个训练样本,由第(N-M)~N张遥感图像构成;
[0042] 在本实施例中,收集待预测地区近20年的单色遥感图像20张,按时间顺序排列并编号为第1~20张,每张图像与前一张图像拍摄时间间隔都近似为1年。对获取的20张遥感图像进行预处理:通过几何校正、缩放、剪切、图像增强等处理完成图形序列的规范化,统一的地理坐标系分辨率和尺寸。将预处理后的图像序列拆分成4帧一组的训练样本,第1组为第1~4张,第2组为第2~5张,直至第17组为第17~20张。
[0043] S2、构建深层网络模型
[0044] 将卷积神经网络CNN和空间变换网络STN级联组成深层网络模型,深层网络模型的输入为某一个训练样本中的前M张图像,输出为预测的遥感图像;
[0045] 在本实施例中,深层网络模型的具体结构为:STN级联在CNN之后;其中,CNN采用UNet的结构,包括编码器和解码器,级数均为N;编码器第n层卷积层提取的图像特征Cn与解码器第N-n层反卷积层的输出Rn同尺寸,Cn与Rn进行通道拼接后,输入至第N-n+1层反卷积层,输出得到升维的特征图Rn-1,然后以此类推,直到第N层反卷积层,最后以C1与R1拼接后的多通道特征图为输入,经反卷积得到遥感图像的演化场Φ;
[0046] 在STN中,基于演化场Φ,将输入的第M张图像Pm转化成预测输出图像Pn,具体转换的公式为:
[0047] Pn(p)=Pm·Φ(p)=∑q∈Z(p')Pm(q)Πd∈{x,y}(1-|p'd-qd|)
[0048] 其中,·表示空间变换方法,基于演化场Φ通过线性插值将图像Pm变换成输出预测图像Pn;p为Pn中的任意一个像素点的像素坐标,p'=p+Φ(p)为根据演化场Φ计算出的p在Pm中所对应的亚像素点坐标;Z(p')为亚像素点p'对应的4个邻域像素点坐标的集合,q为其中任意一个像素点,d表示图像的像素点坐标维度,p'd、qd分别表示p'、q在维度d下的坐标值。
[0049] 进一步的,如图2所示,CNN模块的结构如下:
[0050] 1)第一层为输入层,输入层的神经元大小为串联的输入图像的大小。
[0051] 2)第二层为卷积层,卷积核的维度与输入图像的维度相同,输出特征图为原图像分辨率的1/2。
[0052] 3)第三层为卷积层,卷积核的个数是第二层卷积核个数的一倍,该层输出特征图为原图像分辨率的1/4。
[0053] 4)第四层为卷积层,卷积核的个数与第三层卷积核个数相同,该层输出特征图为原图像分辨率的1/8。
[0054] 5)第五层为卷积层,卷积核的个数与第三层卷积核个数相同,该层输出特征图为原图像分辨率的1/16。
[0055] 6)第六层为反卷积层,输入为第五层输出,输出特征图为原图像分辨率的1/8。
[0056] 7)第七层为反卷积层,第六层输出特征图与第四层输出特征图串联后输入到该层中,输出特征图为原图像分辨率的1/4。
[0057] 8)第八层为反卷积层,第七层输出特征图与第三层输出特征图串联后输入到该层中,输出特征图为原图像分辨率的1/2。
[0058] 9)第九层为反卷积层,第八输出特征图与第二层输出特征图串联后输入到该层中,输出为预测的覆被演化场。
[0059] S3、如图3所示,训练深层网络模型
[0060] S3.1、深层网络模型初始化,设定阈值θ以及学习速率η;
[0061] S3.2、随机选取一个训练样本,将训练样本中的前M帧图像输入至深层网络模型中;
[0062] 在CNN中,先利用编码器的第一层提取输入图像的特征,再从第二层开始依次对图像特征进行降维处理,得到降维后的特征图;然后利用解码器对特征图进行升维处理,最终输出覆被演化场Φ;
[0063] 在STN中,基于演化场Φ,利用STN对输入的第M张图像Pm转化成预测输出图像Pn;
[0064] S3.3、计算预测输出图像Pn与该样本中第M+1张图像Pm+1之间的像素误差平方和,记为Σ;计算演化场Φ的梯度值ΔΦ;再将Σ与ΔΦ的和作为该训练样本的损失函数值;
[0065] S3.4、重复步骤S3.2~S3.3,直至计算出所有训练样本的损失函数值,将所有训练样本损失函数值相加,作为深层网络模型的损失函数值,记为Ω;
[0066] S3.5、比较Ω与θ的大小,如果Ω<θ,则终止训练,得到训练好的深层网络模型,否则,进入步骤S3.6;
[0067] S3.6、利用批量梯度下降算法在损失函数负梯度方向,根据负梯度方向以学习速率η更新深层神经网络模型中的权重参数,然后返回步骤S3.2,再次利用训练样本集进行下一轮训练;
[0068] 在本实施例中,在训练集中抽取1组训练样本,4帧图像重新编号为P1-P4;将P1-P3以3通道图像的形式输入到CNN模块中,输出预测P4的演化场Φ;基于演化场Φ中的坐标变换量,空间变换模块从P3中通过复制或插值获取像素信息得到预测的第4张遥感图像 将预测的第4张遥感图像与真实的第4张图像进行比对,基于损失函数,应用随机梯度下降算法对模型参数进行更新,完成该样本的训练;最后重复以上步骤,遍历训练集中的所有样本完成一轮训练,重复多轮,直到获得满意的损失函数值,完成网络训练,获得训练好的模型。
[0069] S4、地表覆被预测
[0070] 采集M张待预测地区的遥感图像并按时间顺序编号,然后将编号后的遥感图像序列输入训练好的深层网络模型,其输出图像为该预测地区未来的地表覆被预测图像。
[0071] 尽管上面对本发明说明性的具体实施方式进行了描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈