一种基于可靠性排序的条纹方向角快速解包方法技术

本发明专利技术公开了一种基于可靠性排序的条纹方向角快速解包方法，该方法包括：首先利用干涉图计算梯度，再计算密度图，由密度图计算可靠性值，对可靠性值进行排序，然后按照可靠性值大小顺序依次对其相邻的两个像素点方向角相对解包，以此得到指向图。本发明专利技术采用的展开路径为非连续路径，解决了传统算法在处理大尺寸图像时受到连续路径限制处理速度慢的问题。寸图像时受到连续路径限制处理速度慢的问题。寸图像时受到连续路径限制处理速度慢的问题。

【技术实现步骤摘要】
一种基于可靠性排序的条纹方向角快速解包方法


[0001]本专利技术属于图像处理
，特别是一种基于可靠性排序的条纹方向角快速解包算法。

技术介绍

[0002]条纹方向包含着干涉条纹图的重要信息，是许多现有去噪算法和单帧解调技术的重要组成部分。条纹方向角一般通过条纹强度梯度获得，由于使用条纹强度梯度替代实际相位梯度求取方向角，导致方向角包裹在[
‑
π/2,π/2]之间，方向图解包就是把包裹在[
‑
π/2,π/2]之间的方向角还原成真实的指向图的过程。
[0003]为了很好的实现方向角解包运算，国内外学者已经提出了诸如正则化相位跟踪的解包方法，质量图引导的解包方法。上述两种算法都是基于先前解包的最可靠像素的解包算法，这一类算法归于连续路径算法。上述算法在处理大尺寸图像时受到连续路径的限制，解包时间长，无法满足单帧干涉图处理的实时性需求，因此，实有必要提供一种非连续路径的快速解包算法。

技术实现思路

[0004]本专利技术的目的在于针对上述现有技术存在的缺陷和不足，提供一种解包速度快，抗干扰能力强的方向图快速解包方法。
[0005]实现本专利技术目的的技术解决方案为：一方面，提供了一种基于可靠性排序的条纹方向角快速解包方法，包括以下步骤：
[0006]步骤1，计算干涉图条纹梯度值；
[0007]步骤2，根据条纹梯度值计算条纹密度图；
[0008]步骤3，利用条纹密度图计算垂直可靠性值和水平可靠性值；
[0009]步骤4，通过可靠性值位置引导其相邻两个像素点进行相对解包；
[0010]步骤5，重复执行步骤4直至解包结束,得到指向图。
[0011]进一步地，步骤1所述计算干涉图条纹梯度值，具体为：
[0012]计算干涉图像的水平梯度其中f(i,j)、f(i
‑
1,j)分别为干涉图像在(i,j)、(i
‑
1,j)位置处的像素值；
[0013]计算干涉图像的垂直梯度
[0014]进一步地，步骤2所述根据条纹梯度值计算条纹密度图，具体为：
[0015][0016]式中，D
q
(i,j)表示干涉图像在(i,j)位置处的条纹密度值，符号表示卷积运算符，K
ρ
表示核大小为ρ的高斯滤波器。
[0017]进一步地，步骤3所述根据条纹密度图计算水平可靠性值和垂直可靠性值，具体为：
[0018]计算水平可靠性值E
h
：E
h
(i,j)＝D
q
(i,j)+D
q
(i+1,j)；
[0019]计算垂直可靠性值E
v
：E
v
(i,j)＝D
q
(i,j)+D
q
(i,j+1)；
[0020]式中，E
h
(i,j)表示干涉图像在(i,j)位置处的水平可靠性值，E
v
(i,j)表示图像在(i,j)位置处的垂直可靠性值，D
q
(i+1,j)表示干涉图像在(i+1,j)位置处的条纹密度值。
[0021]进一步地，步骤4所述通过可靠性值位置引导其相邻两个像素点进行相对解包，具体包括：
[0022]解包前所有像素点单独分组：查找可靠性值相邻的两个像素点A、B，假设A属于G
A
组、B属于G
B
组，计算两点方向角差值Δ
β
，判断A、B分组是否相同，如果相同则跳过此可靠性值，否则判断cos(Δ
β
)是否大于零，若大于零则令α＝0，否则令α＝1；
[0023]判断G
A
组元素数量是否大于G
B
组元素数量，若是，则G
B
组内所有方向角加上α
×
π并加入至G
A
组，否则G
A
组内所有方向角加上α
×
π并加入至G
B
组。
[0024]进一步地，步骤5所述重复执行步骤4直至解包结束,得到指向图，具体为：依次循环所有边缘值对全部方向角解包，并将解包结果包裹在[
‑
π,π]之间，即为指向图。
[0025]进一步地，在步骤5执行之前，还包括执行：
[0026]对步骤3获得的可靠性值进行降序排列。
[0027]另一方面，提供了一种基于可靠性排序的条纹方向角快速解包系统，所述系统包括顺次执行的：
[0028]第一模块，用于计算干涉图条纹梯度值；
[0029]第二模块，用于根据条纹梯度值计算条纹密度图；
[0030]第三模块，用于利用条纹密度图计算垂直可靠性值和水平可靠性值；
[0031]第四模块，用于通过可靠性值位置引导其相邻两个像素点进行相对解包；
[0032]第五模块，用于重复运行第四模块直至解包结束,得到指向图。
[0033]进一步地，所述系统还包括在第四模块执行前的：
[0034]第六模块，用于对第三模块获得的可靠性值进行降序排列。
[0035]本专利技术与现有技术相比，其显著优点为：
[0036]1)先进行高质量的像素点相对展开，后进行低质量的像素点相对展开，避免了展开路径先经过低质量区域产生的传播误差。
[0037]2)该算法为非连续路径算法，减少了相对已展开像素点选取下一展开像素点的操作，降低了运算时间。
[0038]下面结合附图对本专利技术作进一步详细描述。
附图说明
[0039]图1为本专利技术基于可靠性排序的条纹方向角快速解包算法流程图。
[0040]图2为通过可靠性值位置引导解包的算法流程图。
[0041]图3为一个实施例中去除背景分量和噪声分量的干涉条纹图。
[0042]图4为一个实施例中采用梯度算法对图3计算的条纹方向图。
[0043]图5为一个实施例中对图3求得的条纹密度图。
[0044]图6为图5局部的可靠性值计算示例图。
[0045]图7为一个实施例中最终解包后的指向图。
具体实施方式
[0046]为了使本申请的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本申请进行进一步详细说明。应当理解，此处描述的具体实施例仅仅用以解释本申请，并不用于限定本申请。
[0047]需要说明的是，本专利技术的说明书和权利要求书及上述附图中的术语“包括”和“具体”以及他们的任何变形，意图在于覆盖不排他的包括，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
[0048]需要说明的是，若本专利技术实施例中有涉及方向性指示(诸如上、下、左、右、前、后
……
)，则该方向性指示仅用于解释本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于可靠性排序的条纹方向角快速解包方法，其特征在于，包括以下步骤：步骤1，计算干涉图条纹梯度值；步骤2，根据条纹梯度值计算条纹密度图；步骤3，利用条纹密度图计算垂直可靠性值和水平可靠性值；步骤4，通过可靠性值位置引导其相邻两个像素点进行相对解包；步骤5，重复执行步骤4直至解包结束,得到指向图。2.根据权利要求1所述的基于可靠性排序的条纹方向角快速解包方法，其特征在于，步骤1所述计算干涉图条纹梯度值，具体为：计算干涉图像的水平梯度其中f(i,j)、f(i
‑
1,j)分别为干涉图像在(i,j)、(i
‑
1,j)位置处的像素值；计算干涉图像的垂直梯度3.根据权利要求2所述的基于可靠性排序的条纹方向角快速解包方法，其特征在于，步骤2所述根据条纹梯度值计算条纹密度图，具体为：式中，D
q
(i,j)表示干涉图像在(i,j)位置处的条纹密度值，符号表示卷积运算符，K
ρ
表示核大小为ρ的高斯滤波器。4.根据权利要求3所述的基于可靠性排序的条纹方向角快速解包方法，其特征在于，步骤3所述根据条纹密度图计算水平可靠性值和垂直可靠性值，具体为：计算水平可靠性值E
h
：E
h
(i,j)＝D
q
(i,j)+D
q
(i+1,j)；计算垂直可靠性值E
v
：E
v
(i,j)＝D
q
(i,j)+D
q
(i,j+1)；式中，E
h
(i,j)表示干涉图像在(i,j)位置处的水平可靠性值，E
v
(i,j)表示图像在(i,j)位置处的垂直可靠性值，D
q
(i+1,j)表示干涉图像在(i+1,j)位置处的条纹密度值。5.根据权利要求4所述的基于可靠性排序的条纹方向角快速解包方法，其特征在于，步骤4所述通过可靠性值位置引导其相邻两个像素点进行相对解包，具体包括：...

【专利技术属性】
技术研发人员：李建欣，李志宇，
申请(专利权)人：南京理工大学，
类型：发明
国别省市：