发明名称 基于深度稀疏自编码的多尺度几何遥感图像融合方法
摘要 本发明公开了一种基于深度稀疏自编码的多尺度几何遥感图像融合。其方案为:1)输入低分辨率多光谱图像M和高分辨率全色图像P,提取M第一主成分C1;2)获得C1的低通系数L<sub>M</sub>、带通系数H<sub>M</sub>和P的低通系数L<sub>P</sub>、带通系数H<sub>P</sub>;3)融合L<sub>M</sub>和L<sub>P</sub>得到低通系数L<sub>N</sub>;4)构造空间自相似字典D<sub>S</sub>,并在字典D<sub>S</sub>下融合H<sub>M</sub>和H<sub>P</sub>获得带通系数H<sub>N</sub>;5)更新H<sub>N</sub>得到带通系数<img file="DDA0001044169120000011.GIF" wi="98" he="62" />6)对L<sub>N</sub>和<img file="DDA0001044169120000012.GIF" wi="74" he="63" />逆变换得到融合第一主成分C2,并更新C2得到更新第一主成分C<sub>S</sub>;7)对C<sub>S</sub>逆变换得到高分辨多光谱图像。本发明减少了不匹配细节的注入,改善了融合后多光谱图像空间扭曲,可用于目标识别、地物分类和遥感监测。
申请公布号 CN106204450A 申请公布日期 2016.12.07
申请号 CN201610537055.5 申请日期 2016.07.08
申请人 咸阳师范学院 发明人 李红;苏晓萌;雷亮;宋笑雪;吴粉侠;刘小豫;段群;韩丽娜
分类号 G06T3/40(2006.01)I;G06T5/50(2006.01)I 主分类号 G06T3/40(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 基于深度稀疏自编码的多尺度几何遥感图像融合,包括:(1)分别输入低分辨率多光谱图像和高分辨率全色图像,并提取低分辨多光谱图像的第一主成分C1;(2)通过非采样轮廓波变换,分别获得低分辨多光谱图像的第一主成分低通系数L<sub>M</sub>与带通系数H<sub>M</sub>,及高分辨全色图像的低通系数L<sub>P</sub>与带通系数H<sub>P</sub>;(3)将低分辨多光谱图像的第一主成分低通系数L<sub>M</sub>与全色图像的低通系数L<sub>P</sub>融合,获得融合后的低通系数L<sub>N</sub>;(4)获得空间自相似字典D<sub>S</sub>:4a)将多光谱图像第一主成分的带通系数H<sub>M</sub>和全色图像的带通系数H<sub>P</sub>进行分块,得到多光谱图像第一主成分的带通系数小块M1和全色图像的带通系数小块P1;4a)分别对每一个光谱图像第一主成分的带通系数小块M1和全色图像的带通系数小块P1进行近邻扩充,得到扩充后的第一主成分带通系数I<sub>M</sub>和扩充后的全色图像的带通系数I<sub>P</sub>,该I<sub>M</sub>和I<sub>P</sub>的中心分别为H<sub>M</sub>和H<sub>P</sub>;4b)分别对扩充后的多光谱图像的第一主成分带通系数I<sub>M</sub>和扩充后的全色图像的带通系数I<sub>P</sub>以1为步长划窗取块,再将块列向量化作为字典原子,将扩充后的第一主成分带通系数I<sub>M</sub>的字典原子进行组合得到多光谱图像第一主成分字典D<sub>M</sub>,将扩充后的全色图像的带通系数I<sub>P</sub>的字典原子进行组合得到全色字典D<sub>P</sub>;将多光谱图像第一主成分字典D<sub>M</sub>与全色字典D<sub>P</sub>进行组合,获得空间自相似字典D<sub>S</sub>;(5)将低分辨多光谱图像的第一主成分带通系数H<sub>M</sub>与全色图像的带通系数H<sub>P</sub>进行融合,获得融合后的带通系数H<sub>N</sub>;(6)根据融合后的带通系数H<sub>N</sub>获得多级残差稀疏自编码的带通系数<img file="FDA0001044169090000011.GIF" wi="99" he="64" />6a)计算融合后的带通系数H<sub>N</sub>分别与低分辨多光谱图像的第一主成分带通系数H<sub>M</sub>的残差R1和与全色图像的带通系数H<sub>P</sub>的残差R2;6b)将步骤6a)获得的两个残差R1和R2按照区域能量最大的规则进行融合,获得带通残差系数<img file="FDA0001044169090000021.GIF" wi="87" he="64" />6c)将带通残差系数<img file="FDA0001044169090000022.GIF" wi="64" he="64" />与融合后的带通系数H<sub>N</sub>相加,获得残差稀疏自编码带通系数<img file="FDA0001044169090000023.GIF" wi="98" he="64" />6d)重复6a)‑6c)共S次,获得多级残差稀疏自编码带通系数<img file="FDA0001044169090000024.GIF" wi="99" he="64" />1≤S≤4;(7)对融合后的低通系数L<sub>N</sub>和多级残差稀疏自编码带通系数<img file="FDA0001044169090000025.GIF" wi="75" he="64" />同时进行非采样轮廓波变换的逆变换,获得融合后的第一主成分C2;(8)获得空域多级深度稀疏自编码主成分C<sub>S</sub>:8a)计算融合后的第一主成分C2与低分辨多光谱图像第一主成分C1的残差R3,计算融合后的第一主成分C2与全色图像的残差R4;8b)将8a)中获得的两个残差R3和R4按照深度稀疏自编码的融合方法得到融合后的空域残差<img file="FDA0001044169090000026.GIF" wi="87" he="64" />8c)将融合后的空域残差<img file="FDA0001044169090000027.GIF" wi="64" he="64" />与融合后的第一主成分C2相加获得空域深度稀疏自编码主成分C<sub>1</sub>;8d)重复8a)‑8c)共S次,获得空域多级深度稀疏自编码主成分C<sub>S</sub>,1<S≤4;(9)将空域多级深度稀疏自编码主成分C<sub>S</sub>进行逆主成分分析变换,获得高分辨多光谱图像。
地址 712000 陕西省咸阳市渭城区文林路东段