买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

申请/专利权人：无锡学院

摘要：本发明公开了一种基于非对称多层分解的ISAR与VIS图像融合方法，载入空间分辨率相同的逆合成孔径雷达图像与可见光图像，比较逆合成孔径雷达图像与可见光图像的加权空间频率方差的大小，将两幅图像分为细节图像Ia和粗糙图像Ib；使用多层高斯边窗滤波器分解框架对Ia和Ib分别进行分解，获得Ia的细节保留层Sda、边缘保留层Sea、基本能量层Sga、Ib的细节保留层Sdb、边缘保留层Seb和基本能量层Sgb；通过获得的Sda对Sdb进行引导融合策略获得Ib最终的非对称细节保留融合层Sfb；使用局部方差与空间频率构造判别标准对Sda和Sfb进行融合，获得最终的细节保留融合层Sfd；利用ω对Sea和Seb进行融合，获得最终的边缘保留融合层Sfe；将Sga与Sgb融合，获得最终的基本能量层Sfg；将Sfd，Sfe和Sfg相加获得最终的融合图像If。

主权项：1.一种基于非对称多层分解的ISAR与VIS图像融合方法，其特征在于，该方法为：步骤一、载入空间分辨率相同的逆合成孔径雷达图像ISAR：I1和可见光图像VIS：I2；步骤二、确定I1的空间频率F1，并且通过F1确定I1的加权空间频率方差σ1；确定I2的空间频率F2，并且通过F2确定I2的加权空间频率方差σ2；具体通过以下步骤实现：201按照下式确定I1的行频率R1和I1的列频率C1为 其中，M和N分别表示I1的长度和宽度，I1x,y表示I1在第x行第y列像素点的灰度值，x∈{1,...,M},y∈{1,...,N}，∑·表示求和操作；202按照下式确定I2的行频率R2和I2的列频率C2为 其中，M和N分别表示I2的长度和宽度，I2x,y表示I2在第x行第y列像素点的灰度值，x∈{1,...,M},y∈{1,...,N}；203按照下式确定I1的空间频率F1和I2的空间频率F2为 利用步骤201～203，即可完成一幅图像的求空间频率操作；204按照下式确定I1的灰度平均值A1和I2的灰度平均值A2为 205按照下式确定I1的灰度值方差V1和I2的灰度值方差V2为 206按照下式确定I1的加权空间频率方差σ1和I2的加权空间频率方差σ2为 步骤三、比较σ1与σ2的大小，如果σ1大于等于σ2，则将I1记为细节图像Ia，I2记为粗糙图像Ib；如果σ1小于σ2，则将I2记为细节图像Ia，I1记为粗糙图像Ib；步骤四、通过多层高斯边窗滤波器分解框架对Ia和Ib分别进行多层分解，获得Ia的细节保留层Sda、Ia的边缘保留层Sea、Ia的基本能量层Sga、Ib的细节保留层Sdb、Ib的边缘保留层Seb和Ib的基本能量层Sgb；具体通过以下步骤实现：301使用多层高斯边窗滤波器分解框架中的高斯高频提取方法对Ia和Ib分别进行分解，按照下式确定Ia通过高斯滤波器的第i个低频信息Igi，Ib通过高斯滤波器的第m个低频信息Ifm为 其中，i∈{1,2,3}，m∈{1,2,3}，Ig0＝Ia，If0＝Ib，GF·表示进行高斯滤波操作；302使用多层高斯边窗滤波器分解框架中的边窗高频提取方法对Ia和Ib分别进行分解，按照下式确定Ia通过边窗滤波器的第j个低频信息Idj，和Ib通过边窗滤波器的第n个低频信息Ihn为 其中，j∈{1,2,3}，n∈{1,2,3}，Id0＝Ia，Ih0＝Ib，SWF·表示进行边窗滤波操作，r表示滤波窗口的半径，e表示滤波器迭代的次数；303按照下式确定Ia的细节保留层Sda和Ib的细节保留层Sdb为 304按照下式确定Ia的边缘保留层Sea和Ib的边缘保留层Seb为 305按照下式确定Ia的基本能量层Sga和Ib的基本能量层Sgb为 步骤五、通过Sda对Sdb进行引导融合，获得Ib的细节保留融合层Sfb；步骤六、确定Sda的局部方差Ldak,o和空间频率Fda，确定Sfb的局部方差Lfbk,o和空间频率Ffb，通过Ldak,o、Lfbk,o、Fda和Ffb构造判别标准对Sda和Sfb进行融合，获得最终的细节保留融合层Sfd；具体通过以下步骤实现：501按照下式确定Sda的局部方差Ldak,o和Sfb的局部方差Lfbk,o为 其中，在Sda中随机生成一个长度为P，宽度为Q的大小的区域，在Sfb的相同位置处也生成一个长度为P，宽度为Q的大小的区域，k表示区域中心位置的横坐标，o表示区域中心位置的纵坐标，w表示区域中的横坐标像素索引，z表示区域中的纵坐标像素索引，μda表示Sda在此区域中像素灰度值的平均值，μfb表示Sfb在此区域中像素灰度值的平均值；502按照下式确定Sda的空间频率Fda和Sfb的空间频率Ffb为 其中，Cda表示Sda的列频率，Rda表示Sda的行频率，Cfb表示Sfb的列频率，Rfb表示Sfb的行频率；503按照下式确定最终的细节保留融合层Sfd为 步骤七、确定Sea的加权空间频率方差σea和Seb的加权空间频率方差σeb，通过σea和σeb确定权重系数ω，通过ω对Sea和Seb进行加权融合，获得最终的边缘保留融合层Sfe；步骤八、将Sga与Sgb融合，获得最终的基本能量融合层Sfg；步骤九、通过Sfd、Sfe和Sfg相加的和确定最终的融合图像If。

全文数据：

权利要求：

百度查询： 无锡学院 基于非对称多层分解的ISAR与VIS图像融合方法