CN109472834B - 一种基于小波变换的卡尔曼滤波相位展开方法 - Google Patents
一种基于小波变换的卡尔曼滤波相位展开方法 Download PDFInfo
- Publication number
- CN109472834B CN109472834B CN201811234118.5A CN201811234118A CN109472834B CN 109472834 B CN109472834 B CN 109472834B CN 201811234118 A CN201811234118 A CN 201811234118A CN 109472834 B CN109472834 B CN 109472834B
- Authority
- CN
- China
- Prior art keywords
- phase
- pixel
- pixels
- order
- interference
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 68
- 238000001914 filtration Methods 0.000 title claims abstract description 36
- 230000009466 transformation Effects 0.000 title claims abstract description 9
- 239000011159 matrix material Substances 0.000 claims abstract description 24
- 238000010587 phase diagram Methods 0.000 claims abstract description 20
- 238000004804 winding Methods 0.000 claims description 55
- 238000000354 decomposition reaction Methods 0.000 claims description 21
- 238000010586 diagram Methods 0.000 claims description 18
- 230000008569 process Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 2
- 230000007480 spreading Effects 0.000 claims description 2
- 238000012545 processing Methods 0.000 description 5
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000009965 odorless effect Effects 0.000 description 2
- CCEKAJIANROZEO-UHFFFAOYSA-N sulfluramid Chemical group CCNS(=O)(=O)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)F CCEKAJIANROZEO-UHFFFAOYSA-N 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T9/00—Image coding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/10—Image enhancement or restoration using non-spatial domain filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20024—Filtering details
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20048—Transform domain processing
- G06T2207/20064—Wavelet transform [DWT]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Multimedia (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种基于小波变换的卡尔曼滤波相位展开方法,涉及相位展开技术领域,解决的技术问题是解决卡尔曼滤波相位展开方法中运算成本高以及对计算机内存要求较高问题,该方法包括以下步骤:(1)对干涉图进行N次分解;(2)根据高频相位梯度信息创建掩膜矩阵A;(3)展开最高阶低频干涉相位图sN;(4)判断相位展开是否完成,检查当前阶是否为最低阶s0,是,结束;否则,转入下一步,继续展开较低阶缠绕相位图;(5)展开较低阶缠绕相位图;(6)判断相位展开是否完成,检查当前阶是否为最低阶s0,是,结束;否则,N=N‑1,转第(5)步。本发明在相位展开中运算成本低,在维持高效率的情况下,依然保证了精度,稳健性较强。
Description
技术领域
本发明涉及相位展开技术领域,尤其涉及一种基于小波变换的卡尔曼滤波相位展开方法。
背景技术
相位展开是指把干涉纹图中的缠绕相位加上2π的整数倍还原出真实相位的过程。传统经典相位展开方法总体分为三类,第一类是基于路径跟踪策略,包含了Goldstein枝切法、质量图引导法和最小不连续相位展开方法;第二类是最小范数方法,代表是最小二乘相位展开法;第三类是基于网络流的各种方法。还有一类方法可以展开相位的同时进行滤波,如卡尔曼滤波类相位展开方法。非线性卡尔曼滤波能在状态估计的同时实现噪声滤波,具有较强噪声抑制能力,将传统质量引导策略与无味卡尔曼滤波方法(UKF)相结合,可以利用相位质量图来引导相位展开路径,使无味卡尔曼滤波方法沿高质量干涉图向低质量干涉图进行相位展开,既可以准确的恢复真实相位,又可以避免直接穿过噪声较大的像元区域导致精度下降。但是卡尔曼滤波方法比传统相位展开方法拥有较高精度的同时,也需要比其他大部分方法更多的运行时间,随着干涉图规模的增大,尤其是针对大规模低信噪比干涉相位图进行相位展开时,质量引导策略的无味卡尔曼滤波相位展开方法其计算成本将急剧增大,故难以直接应用到大规模干涉图相位展开中。
发明内容
针对现有技术的不足,本发明所解决的技术问题是解决卡尔曼滤波相位展开方法中运算成本高以及对计算机内存要求较高问题。
为解决上述技术问题,本发明采用的技术方案是一种基于小波变换的卡尔曼滤波相位展开方法,包括以下步骤:
(1)对干涉图进行N次分解,具体过程如下:
根据干涉图的规模与噪声选取合适的小波分解阶数,然后依次对干涉图的行方向和列方向进行小波分解,高频细节部分为干涉相位梯度信息,之后的分解只需对低频缠绕相位图处理,分解公式如下:
di+1(n)=W(si(2n+1)-si(2n))
其中,si+1(n)为第i+1层分解后低频信号,di+1(n)为第i+1层分解后高频细节,W(·)为缠绕算子,保证每一层分解结果在相位主值区间。
(2)根据高频相位梯度信息创建掩膜矩阵A,具体过程如下:
干涉相位图经过小波分解后得到每一阶的低频缠绕相位与高频梯度信息,低频缠绕相位保存为需要进行相位展开的待展开相位,而由于相位梯度需满足|Δi+1(n)|≤π的条件,所以对于噪声较大的缠绕相位需要在更低阶进行相位展开,计算公式如下:
di+1(n)=W(si(2n+1)-si(2n))
由公式可以看出:低阶缠绕相位梯度不大于π/2时,高一阶低频缠绕相位梯度可满足约束条件|Δi+1(n+1)|≤π,其对应缠绕相位可在高阶展开,否则需在较低阶进行相位展开,即掩膜矩阵A由高频部分梯度信息是否大于π/2来生成,缠绕相位展开情况与掩膜矩阵元素对应值如下所示:
0:当前相位已经在更高阶被展开;
1:当前相位应在当前阶被展开;
2:当前相位应在更低阶被展开;
低频缠绕相位是否展开取决于对应高频梯度信息,即取决于掩膜矩阵,掩膜矩阵的正确与否直接影响到最终相位展开的精确度与效率。
(3)展开最高阶低频干涉相位图sN,具体分步骤如下:
1)获取最高阶低频缠绕干涉相位图sN的相位质量图,指导相位展开方法展开路径;
2)根据所得干涉相位质量图,选取干涉图中相位质量最高的非边界像元作为起始像元,其展开相位为缠绕相位,并设定其估计误差方差,估计误差方差值的范围为(0,1),根据经验设定;
3)在干涉图中以起始像元为中心,标记与起始像元邻接的上、下、左和右4个缠绕像元为邻接像元,将所有在掩膜矩阵AN中对应位置值‘1’的邻接像元嵌入待展开像元组;
4)比较所得待展开像元组中所有像元对应的像元质量值,选取最佳待展开像元标记为像元x,然后利用卡尔曼滤波相位展开方法展开此缠绕像元,具体展开步骤如下;
设所述最佳待展开像元x为干涉缠绕相位图中的(m,n)像元,卡尔曼滤波相位展开方法利用干涉图中邻接待展开像元(m,n)的上、下、左、右、上左、下左、上右和下右8个像元中的已展开像元的信息递推估计待展开像元像元(m,n)的干涉相位;设Ψ为待展开像元(m,n)的8个邻接像元中已展开像元的集合,设像元(k,l)是干涉图中待展开像元像元(m,n)的8个邻接像元中的已展开像元,即(k,l)∈Ψ,像元(k,l)的状态估计及误差方差分别为和
首先,对干涉图中的像元(m,n)干涉相位值按如下卡尔曼滤波方法进行预测估计:
上式中,(m,n)像元为干涉图中的待展开像元,χj,(k,l)是(k,l)像元状态估计的Sigma point,为干涉图中(m,n)像元与(k,l)像元之间的相位梯度估计值,SNR(k,l)是干涉图中(k,l)像元信噪比,为干涉图中(m,n)像元Sigma point预测值,Q(m,n)|(k,l)为干涉图中的(m,n)像元与(k,l)像元之间的相位梯度估计误差方差,和为相应调节权值系数;
其次,对干涉图中的(m,n)像元状态进行估计:
其中,y(m,n)和v(m,n)分别为干涉图中(m,n)像元含噪声的观测矢量及其附加噪声矢量,u(m,n)和v(m,n)分别为附加在复干涉信号实部和虚部的噪声,h[x(m,n)]为干涉图中(m,n)像元真实观测矢量,表示干涉图中像元(m,n)观测矢量预测值,Π(m,n)表示干涉图中(m,n)像元增益矩阵,R(m,n)表示干涉图中的(m,n)像元观测误差方差,和分别为干涉图中的(m,n)像元状态估计及误差方差;
5)从分步骤3)步获得的待展开像元组中删除像元x,标记sN中相邻x的非边界缠绕像元为邻接像元,将所有在掩膜矩阵中对应位置值为‘1’的邻接像元嵌入待展开像元组;
6)检查待展开像元组是否为空,若没有待展开像元,干涉图sN相位展开完成;否则,转入分步骤4)步。
(4)判断相位展开是否完成
检查当前阶是否为最低阶s0,是,结束;否则.转入下一步,继续展开较低阶缠绕相位图;
(5)展开较低阶缠绕相位图,具体分步骤如下:
1)对sN中每一个展开像元和dN中对应的相位梯度进行逆小波变换,得到sN-1中两个对应像元的展开相位,实现对干涉图的sN-1更新;更新之后的干涉相位图sN-1含有已在较高阶展开的像元、应在当前阶展开的像元、应在较低阶展开的像元(最低阶s0只含有前两者);
2)标记中sN-1所有展开像元的非边界相邻缠绕像元为邻接像元组,将所有在AN-1中对应位置为‘1’的邻接像元嵌入待展开像元组;
3)对待展开像元组中质量最高的像元x利用卡尔曼滤波相位展开方法展开;
4)从第2)步获得的待展开像元组中删除像元x,标记sN-1中相邻x的非边界缠绕像元为邻接像元,将所有在掩膜矩阵AN-1中对应位置值为‘1’的邻接像元嵌入待展开像元组;
5)检查待展开像元组是否为空,若没有待展开像元,相位干涉图sN-1展开完成;否则,转入第3)步。
(6)判断相位展开是否完成,检查当前阶是否为最低阶s0,是,结束;否则,N=N-1,转第(5)步。
与现技术相比本发明有益效果:
(1)通过对较大的原始相位干涉图进行小波分解,把对大尺寸原始干涉图的展开任务主要转移到了对干涉图最高阶的低频缠绕相位图的展开中,利用质量图引导策略的卡尔曼滤波方法对最高阶缠绕相位展开,既大大提高了效率,也维持了较高精度;
(2)与基于质量引导策略的卡尔曼滤波相位展开方法比较,数据处理展开方法略低,运行时间远远小于基于质量引导策略的卡尔曼相位展开方法,表明本发明在维持高效率的情况下,依然保证了精度,稳健性较强。
附图说明
图1为本发明流程示意图;
图2为实施例1所针对的多山地形场景模拟相位图;
图3为对图2模拟相位图进行缠绕得到的缠绕相位图;
图4为实施例1对最高阶低频缠绕相位展开结果;
图5为实施例1相位展开结果;
图6为实施例1相位展开误差直方图;
图7为基于质量引导策略的卡尔曼滤波相位展开方法相位展开结果;
图8为基于质量引导策略的卡尔曼滤波相位展开方法相位展开误差直方图;
图9为实施例2针对的部分Etna火山实测干涉图;
图10为实施例2对部分Etna火山实测数据相位展开结果。
具体实施方式
下面结合附图和实施例对本发明的具体实施方式作进一步的说明,但不是对本发明的限定。
实施例1:
采用多山地形场景干涉图进行相位展开,图2为多山地形场景的模拟真实相位图,其横坐标为距离向像素,纵坐标为方位向像素,规模为256×256,灰度表示干涉相位值,单位为弧度;图3为多山地形场景的含噪声缠绕相位图,横、纵坐标与图2相同,其信噪比为6.1dB。
图1示出了一种基于小波变换的卡尔曼滤波相位展开方法,包括以下步骤:
(1)对干涉图进行分解,具体过程如下:
根据噪声与干涉图规模选取小波变换分解阶数为1阶。对图3二维干涉图进行模小波变换分解。先对干涉图行方向进行小波分解,高频信息为相位梯度,不需进行处理,然后对低频缠绕相位进行列方向小波分解,即完成对初始干涉图的1阶小波分解。
(2)根据高频相位梯度信息创建掩膜矩阵A:
经小波分解后的高频系数,也就是梯度,根据梯度是否大于π/2生成掩膜矩阵。对于地形较为陡峭或者噪声较大的区域,将被放在较低阶干涉相位图进行相位展开。
(3)展开最高阶低频干涉相位图s1,具体分步骤如下:
1)获取最高阶干涉相位质量图,反映二维最高阶低频干涉相位图s1像元质量,引导相位展开路径。本实施例采用的干涉相位质量图为相位微分偏差图;
2)根据所得干涉相位质量图,选取干涉图中相位质量最高的非边界像元作为起始像元,其展开相位为缠绕相位,并设定其估计误差方差,本例中设为0.6;
3)在干涉图s1中以起始像元为中心,标记与起始像元邻接的上、下、左和右4个缠绕像元为邻接像元,将所有在掩膜矩阵AN中对应位置值为‘1’的邻接像元嵌入待展开像元组;
4)比较所得待展开像元对应相位质量值,选取最佳待展开像元标记为像元x,然后利用卡尔曼滤波相位展开方法(UKF算法)展开此缠绕像元;
5)从分步骤3)步获得的待展开像元组中删除像元x,标记s1中相邻x的非边界缠绕像元为邻接像元,将所有在掩膜矩阵中对应位置值为‘1’的邻接像元嵌入待展开像元组;
6)检查待展开像元组是否为空,若没有待展开像元,干涉图s1相位展开完成;否则,转入分步骤4)步。干涉相位图s1展开图如图4所示。
(4)判断相位展开是否完成
当前展开的为低频干涉相位图s1,,转入下一步,继续展开较低阶缠绕相位图s0;
(5)展开较低阶缠绕相位图,具体分步骤如下:
1)对s1中每一个展开像元和d1中对应的相位梯度进行行方向逆小波变换,得的s0中两个对应像元的展开相位;
2)标记中s0所有展开像元的非边界相邻缠绕像元为邻接像元组,将所有在A0中对应位置为‘1’的邻接像元嵌入待展开像元组;
3)对待展开像元组中质量最高的像元x利用卡尔曼滤波相位展开方法展开;
4)从第2)步获得的待展开像元组中删除像元x,标记s0中相邻x的非边界缠绕像元为邻接像元,将所有在掩膜矩阵A0中对应位置值为‘1’的邻接像元嵌入待展开像元组;
5)检查待展开像元组是否为空,若没有待展开像元,低频缠绕干涉相位图展开完成;否则,转入第3)步。
(6)判断相位展开是否完成,检查当前阶是否为最低阶s0,是,结束;否则,N=N-1,转第(5)步。
图5为本实施例基于小波变换的卡尔曼滤波相位展开方法对图3展开的结果示意图,示意图横纵坐标与图2相同;图6为本实施例相位展开结果的误差直方图,其横坐标为相位展开误差,纵坐标为误差统计频数。
图7为质量引导策略的卡尔曼滤波相位展开方法对图3展开的结果,示意图横纵坐标与图5相同;图8为质量引导策略的卡尔曼滤波相位展开结果的误差直方图,示意图横纵坐标与图6相同。
由图5~图8可以看出本发明实施例方法对干涉相位展开效果良好,具有一定的抗噪抑制能力,且精度较高;表1列出了本发明与质量引导策略的无味卡尔曼滤波相位展开方法在同一MATALB计算环境下展开较大规模干涉相位图所消耗时间,可以看出本发明展开大规模干涉相位图的运行时间远远小于质量引导策略的无味卡尔曼滤波相位展开方法。
表1相位展开方法运行时间
实施例2:
针对Etna火山干涉图进行相位展开,图9为预滤波后的部分Etna火山干涉图。
采用与实施例1所述的步骤(1)至步骤(5)步展开该缠绕干涉图,本例相位展开结果见图10,可以看出本发明方法展开后相位连续且平滑,反映本发明相位展开方法的有效性。将图10重缠绕后相位图条纹与原始缠绕相位图条纹一致,说明本发明方法在实测数据处理中获得有效结果。
本发明是质量引导策略卡尔曼滤波方法与小波变换的相结合的结果,在实测与模拟数据处理实验中其有效性明显。
与现技术相比本发明有益效果:
(1)通过对较大的原始相位干涉图进行小波分解,把对大尺寸原始干涉图的展开任务主要转移到了对干涉图最高阶的低频缠绕相位图的展开中,利用质量图引导策略的卡尔曼滤波方法对最高阶缠绕相位展开,既大大提高了效率,也维持了较高精度;
(2)与基于质量引导策略的卡尔曼滤波相位展开方法比较,数据处理展开方法略低,运行时间远远小于基于质量引导策略的卡尔曼相位展开方法,表明本发明在维持高效率的情况下,依然保证了精度,稳健性较强。
以上结合附图和实施例对本发明的实施方式做出了详细说明,但本发明不局限于所描述的实施方式。对于本领域技术人员而言,在不脱离本发明的原理和精神的情况下,对这些实施方式进行各种变化、修改、替换和变型仍落入本发明的保护范围。
Claims (3)
1.一种基于小波变换的卡尔曼滤波相位展开方法,其特征在于,包括以下步骤:
(1)对干涉图进行N次分解:
根据干涉图的规模与噪声选取合适的小波分解阶数,然后依次对干涉图的行方向和列方向进行小波分解,高频细节部分为干涉相位梯度信息,之后的分解只需对低频缠绕相位图处理,分解公式如下:
di+1(n)=W(si(2n+1)-si(2n))
其中,si+1(n)为第i+1层分解后低频信号,di+1(n)为第i+1层分解后高频细节,为缠绕算子,保证每一层分解结果在相位主值区间。
(2)根据高频相位梯度信息创建掩膜矩阵A:
干涉相位图经过小波分解后得到每一阶的低频缠绕相位与高频梯度信息,低频缠绕相位保存为需要进行相位展开的待展开相位,而由于相位梯度需满足|Δi+1(n)|≤π的条件,所以对于噪声较大的缠绕相位需要在更低阶进行相位展开,计算公式如下:
di+1(n)=W(si(2n+1)-si(2n))
由公式可以看出:低阶缠绕相位梯度不大于π/2时,高一阶低频缠绕相位梯度可满足约束条件|Δi+1(n+1)|≤π,其对应缠绕相位可在高阶展开,否则需在较低阶进行相位展开,即掩膜矩阵A由高频部分梯度信息是否大于π/2来生成,缠绕相位展开情况与掩膜矩阵元素对应值如下所示:
0:当前相位已经在更高阶被展开;
1:当前相位应在当前阶被展开;
2:当前相位应在更低阶被展开;
低频缠绕相位是否展开取决于对应高频梯度信息,即取决于掩膜矩阵,掩膜矩阵的正确与否直接影响到最终相位展开的精确度与效率;
(3)展开最高阶低频干涉相位图sN;
(4)判断相位展开是否完成,检查当前阶是否为最低阶s0,是,结束;否则,转入下一步,继续展开较低阶缠绕相位图;
(5)展开较低阶缠绕相位图;
(6)判断相位展开是否完成,检查当前阶是否为最低阶s0,是,结束;否则,N=N-1,转第(5)步。
2.如权利要求1所述的基于小波变换的卡尔曼滤波相位展开方法,其特征在于,步骤(3)具体分步骤如下:
1)获取最高阶低频缠绕干涉相位图sN的相位质量图,指导相位展开方法展开路径;
2)根据所得干涉相位质量图,选取干涉图中相位质量最高的非边界像元作为起始像元,其展开相位为缠绕相位,并设定其估计误差方差,估计误差方差值的范围为(0,1),根据经验设定;
3)在干涉图中以起始像元为中心,标记与起始像元邻接的上、下、左和右4个缠绕像元为邻接像元,将所有在掩膜矩阵AN中对应位置值‘1’的邻接像元嵌入待展开像元组;
4)比较所得待展开像元组中所有像元对应的像元质量值,选取最佳待展开像元标记为像元x,然后利用卡尔曼滤波相位展开方法展开此缠绕像元,具体展开步骤如下;
设所述最佳待展开像元x为干涉缠绕相位图中的(m,n)像元,卡尔曼滤波相位展开方法利用干涉图中邻接待展开像元(m,n)的上、下、左、右、上左、下左、上右和下右8个像元中的已展开像元的信息递推估计待展开像元像元(m,n)的干涉相位;设Ψ为待展开像元(m,n)的8个邻接像元中已展开像元的集合,设像元(k,l)是干涉图中待展开像元(m,n)的8个邻接像元中的已展开像元,即(k,l)∈Ψ,像元(k,l)的状态估计及误差方差分别为和
首先,对干涉图中的像元(m,n)干涉相位值按如下卡尔曼滤波方法进行预测估计:
上式中,(m,n)像元为干涉图中的待展开像元,χj,(k,l)是(k,l)像元状态估计的Sigmapoint,为干涉图中(m,n)像元与(k,l)像元之间的相位梯度估计值,SNR(k,l)是干涉图中(k,l)像元信噪比,为干涉图中(m,n)像元Sigma point预测值,Q(m,n)|(k,l)为干涉图中的(m,n)像元与(k,l)像元之间的相位梯度估计误差方差,和为相应调节权值系数;
其次,对干涉图中的(m,n)像元状态进行估计:
其中,y(m,n)和ν(m,n)分别为干涉图中(m,n)像元含噪声的观测矢量及其附加噪声矢量,u(m,n)和v(m,n)分别为附加在复干涉信号实部和虚部的噪声,h[x(m,n)]为干涉图中(m,n)像元真实观测矢量,表示干涉图中像元(m,n)观测矢量预测值,Π(m,n)表示干涉图中(m,n)像元增益矩阵,R(m,n)表示干涉图中的(m,n)像元观测误差方差,和分别为干涉图中的(m,n)像元状态估计及误差方差;
5)从分步骤3)步获得的待展开像元组中删除像元x,标记sN中相邻x的非边界缠绕像元为邻接像元,将所有在掩膜矩阵中对应位置值为‘1’的邻接像元嵌入待展开像元组;
6)检查待展开像元组是否为空,若没有待展开像元,干涉图sN相位展开完成;否则,转入分步骤4)步。
3.如权利要求1所述的基于小波变换的卡尔曼滤波相位展开方法,其特征在于,步骤(5)具体分步骤如下:
1)对sN中每一个展开像元和dN中对应的相位梯度进行逆小波变换,得到sN-1中两个对应像元的展开相位,实现对干涉图的sN-1更新;更新之后的干涉相位图sN-1含有已在较高阶展开的像元、应在当前阶展开的像元、应在较低阶展开的像元,最低阶s0只含有前两者;
2)标记中sN-1所有展开像元的非边界相邻缠绕像元为邻接像元组,将所有在AN-1中对应位置为‘1’的邻接像元嵌入待展开像元组;
3)对待展开像元组中质量最高的像元x利用卡尔曼滤波相位展开方法展开;
4)从第2)步获得的待展开像元组中删除像元x,标记sN-1中相邻x的非边界缠绕像元为邻接像元,将所有在掩膜矩阵AN-1中对应位置值为‘1’的邻接像元嵌入待展开像元组;
5)检查待展开像元组是否为空,若没有待展开像元,相位干涉图sN-1展开完成;否则,转入第3)步。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811234118.5A CN109472834B (zh) | 2018-10-23 | 2018-10-23 | 一种基于小波变换的卡尔曼滤波相位展开方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811234118.5A CN109472834B (zh) | 2018-10-23 | 2018-10-23 | 一种基于小波变换的卡尔曼滤波相位展开方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109472834A CN109472834A (zh) | 2019-03-15 |
CN109472834B true CN109472834B (zh) | 2023-04-14 |
Family
ID=65664322
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811234118.5A Active CN109472834B (zh) | 2018-10-23 | 2018-10-23 | 一种基于小波变换的卡尔曼滤波相位展开方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109472834B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111122507B (zh) * | 2019-11-06 | 2022-03-15 | 中国工程物理研究院激光聚变研究中心 | 一种基于离散小波变换的光学平板波前的降噪方法和装置 |
WO2022103587A2 (en) * | 2020-11-09 | 2022-05-19 | Arizona Board Of Regents On Behalf Of The University Of Arizona | Determination of a true shape of an object based on transformation of its optical image |
Citations (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101339245A (zh) * | 2008-08-08 | 2009-01-07 | 西安电子科技大学 | 多基线干涉合成孔径雷达干涉相位展开方法 |
CN102032877A (zh) * | 2010-11-30 | 2011-04-27 | 东南大学 | 基于小波变换的三维测量方法 |
CN102621549A (zh) * | 2011-10-14 | 2012-08-01 | 中国人民解放军国防科学技术大学 | 多基线/多频段干涉相位解缠频域快速算法 |
CN102663695A (zh) * | 2012-03-31 | 2012-09-12 | 重庆大学 | 基于小波变换的dr图像去噪方法及系统 |
CN103530857A (zh) * | 2013-10-31 | 2014-01-22 | 清华大学深圳研究生院 | 基于多尺度的卡尔曼滤波图像去噪方法 |
WO2014118559A1 (en) * | 2013-02-01 | 2014-08-07 | Ucl Business Plc | Apparatus and method for correcting susceptibility artefacts in a magnetic resonance image |
CN104459633A (zh) * | 2014-12-01 | 2015-03-25 | 中国科学院电子学研究所 | 结合局部频率估计的小波域InSAR干涉相位滤波方法 |
CN104766280A (zh) * | 2015-03-26 | 2015-07-08 | 电子科技大学 | 一种基于堆排序的质量图相位解缠方法 |
CN104849714A (zh) * | 2015-05-29 | 2015-08-19 | 桂林电子科技大学 | 一种高效稳健的多基线ukf相位展开方法 |
CN105719253A (zh) * | 2016-01-20 | 2016-06-29 | 桂林电子科技大学 | 一种嵌入堆排序的卡尔曼滤波相位展开方法 |
CN106067163A (zh) * | 2016-05-24 | 2016-11-02 | 中国科学院深圳先进技术研究院 | 一种基于小波分析的图像去雨方法及系统 |
CN106097404A (zh) * | 2016-05-27 | 2016-11-09 | 山东科技大学 | 利用非线性矢量曲面构建InSAR相位图像模型的方法 |
CN106932773A (zh) * | 2017-01-12 | 2017-07-07 | 桂林电子科技大学 | 基于修正嵌入式容积卡尔曼滤波的相位展开算法 |
CN107167087A (zh) * | 2017-05-12 | 2017-09-15 | 天津大学 | 一种基于经验小波变换的物体变形测量方法 |
CN107193005A (zh) * | 2017-06-16 | 2017-09-22 | 桂林电子科技大学 | 一种无损卡尔曼滤波与粒子滤波相结合的相位展开算法 |
CN107798690A (zh) * | 2017-10-18 | 2018-03-13 | 中国科学院生物物理研究所 | 用于活细胞中囊泡运动追踪的方法 |
CN107942331A (zh) * | 2017-11-28 | 2018-04-20 | 西安电子科技大学 | 基于谱分析的多通道sar系统通道偏差估计方法 |
CN107977939A (zh) * | 2017-11-29 | 2018-05-01 | 浙江理工大学 | 一种基于可靠度的加权最小二乘相位展开计算方法 |
CN107977944A (zh) * | 2017-12-12 | 2018-05-01 | 华南农业大学 | 一种估测npp遥感图像数据生成方法 |
CN108052978A (zh) * | 2017-10-17 | 2018-05-18 | 桂林电子科技大学 | 一种基于支持向量机的ukf相位展开算法 |
CN108548502A (zh) * | 2018-03-28 | 2018-09-18 | 西安理工大学 | 一种动态物体三维测量方法 |
CN108596851A (zh) * | 2018-04-26 | 2018-09-28 | 安徽理工大学 | 一种基于小波变化的压缩感知图像处理算法 |
CN108680119A (zh) * | 2018-04-19 | 2018-10-19 | 东南大学 | 一种分区域的单幅快速相位展开方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9933519B2 (en) * | 2010-02-22 | 2018-04-03 | Elbit Systems Ltd. | Three dimensional radar system |
US9020293B2 (en) * | 2012-02-07 | 2015-04-28 | National Cheung Kung University | Integration of filters and phase unwrapping algorithms for removing noise in image reconstruction |
-
2018
- 2018-10-23 CN CN201811234118.5A patent/CN109472834B/zh active Active
Patent Citations (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101339245A (zh) * | 2008-08-08 | 2009-01-07 | 西安电子科技大学 | 多基线干涉合成孔径雷达干涉相位展开方法 |
CN102032877A (zh) * | 2010-11-30 | 2011-04-27 | 东南大学 | 基于小波变换的三维测量方法 |
CN102621549A (zh) * | 2011-10-14 | 2012-08-01 | 中国人民解放军国防科学技术大学 | 多基线/多频段干涉相位解缠频域快速算法 |
CN102663695A (zh) * | 2012-03-31 | 2012-09-12 | 重庆大学 | 基于小波变换的dr图像去噪方法及系统 |
WO2014118559A1 (en) * | 2013-02-01 | 2014-08-07 | Ucl Business Plc | Apparatus and method for correcting susceptibility artefacts in a magnetic resonance image |
CN103530857A (zh) * | 2013-10-31 | 2014-01-22 | 清华大学深圳研究生院 | 基于多尺度的卡尔曼滤波图像去噪方法 |
CN104459633A (zh) * | 2014-12-01 | 2015-03-25 | 中国科学院电子学研究所 | 结合局部频率估计的小波域InSAR干涉相位滤波方法 |
CN104766280A (zh) * | 2015-03-26 | 2015-07-08 | 电子科技大学 | 一种基于堆排序的质量图相位解缠方法 |
CN104849714A (zh) * | 2015-05-29 | 2015-08-19 | 桂林电子科技大学 | 一种高效稳健的多基线ukf相位展开方法 |
CN105719253A (zh) * | 2016-01-20 | 2016-06-29 | 桂林电子科技大学 | 一种嵌入堆排序的卡尔曼滤波相位展开方法 |
CN106067163A (zh) * | 2016-05-24 | 2016-11-02 | 中国科学院深圳先进技术研究院 | 一种基于小波分析的图像去雨方法及系统 |
CN106097404A (zh) * | 2016-05-27 | 2016-11-09 | 山东科技大学 | 利用非线性矢量曲面构建InSAR相位图像模型的方法 |
CN106932773A (zh) * | 2017-01-12 | 2017-07-07 | 桂林电子科技大学 | 基于修正嵌入式容积卡尔曼滤波的相位展开算法 |
CN107167087A (zh) * | 2017-05-12 | 2017-09-15 | 天津大学 | 一种基于经验小波变换的物体变形测量方法 |
CN107193005A (zh) * | 2017-06-16 | 2017-09-22 | 桂林电子科技大学 | 一种无损卡尔曼滤波与粒子滤波相结合的相位展开算法 |
CN108052978A (zh) * | 2017-10-17 | 2018-05-18 | 桂林电子科技大学 | 一种基于支持向量机的ukf相位展开算法 |
CN107798690A (zh) * | 2017-10-18 | 2018-03-13 | 中国科学院生物物理研究所 | 用于活细胞中囊泡运动追踪的方法 |
CN107942331A (zh) * | 2017-11-28 | 2018-04-20 | 西安电子科技大学 | 基于谱分析的多通道sar系统通道偏差估计方法 |
CN107977939A (zh) * | 2017-11-29 | 2018-05-01 | 浙江理工大学 | 一种基于可靠度的加权最小二乘相位展开计算方法 |
CN107977944A (zh) * | 2017-12-12 | 2018-05-01 | 华南农业大学 | 一种估测npp遥感图像数据生成方法 |
CN108548502A (zh) * | 2018-03-28 | 2018-09-18 | 西安理工大学 | 一种动态物体三维测量方法 |
CN108680119A (zh) * | 2018-04-19 | 2018-10-19 | 东南大学 | 一种分区域的单幅快速相位展开方法 |
CN108596851A (zh) * | 2018-04-26 | 2018-09-28 | 安徽理工大学 | 一种基于小波变化的压缩感知图像处理算法 |
Non-Patent Citations (3)
Title |
---|
Efficient multiscale phase unwrapping methodology with modulo wavelet transform;Blinder David;《Optics Express》;20160926;1-16 * |
改进的扩展卡尔曼滤波相位解缠算法;谢先明;《测绘科学》;20181020;16-21 * |
结合局部频率估计的小波域InSAR相位滤波新方法;李芳芳;《系统工程与电子技术》;20150707;2719-2724 * |
Also Published As
Publication number | Publication date |
---|---|
CN109472834A (zh) | 2019-03-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109633648B (zh) | 一种基于似然估计的多基线相位估计装置及方法 | |
CN107193005B (zh) | 一种无损卡尔曼滤波与粒子滤波相结合的相位展开算法 | |
CN109472834B (zh) | 一种基于小波变换的卡尔曼滤波相位展开方法 | |
CN109752715B (zh) | 一种sar数据全散射体探测方法及装置 | |
CN113589286B (zh) | 基于D-LinkNet的无迹卡尔曼滤波相位解缠方法 | |
CN113311433B (zh) | 一种质量图和最小费用流结合的InSAR干涉相位两步解缠方法 | |
CN108362200B (zh) | 一种快速更新InSAR变形序列结果的方法 | |
CN110187609A (zh) | 一种计算光刻的深度学习方法 | |
CN103823219B (zh) | 自适应迭代的非局部干涉合成孔径雷达干涉相位滤波方法 | |
CN108629746B (zh) | 一种基于相关损失卷积神经网络的雷达图像斑噪抑制方法 | |
CN115205360A (zh) | 复合条纹投影钢管的三维外轮廓在线测量与缺陷检测方法及应用 | |
Li | Using radial basis functions in airborne gravimetry for local geoid improvement | |
CN112859077B (zh) | 一种多级合成孔径雷达干涉相位解缠方法 | |
Gao et al. | Two-dimensional phase unwrapping method using a refined D-LinkNet-based unscented Kalman filter | |
CN106932773A (zh) | 基于修正嵌入式容积卡尔曼滤波的相位展开算法 | |
CN104714372B (zh) | 一种用于光刻机的光源掩模优化方法 | |
Kemao | On window size selection in the windowed Fourier ridges algorithm | |
CN112665529B (zh) | 基于条纹密度区域分割和校正的物体三维形貌测量方法 | |
CN110440935B (zh) | 一种基于扩展信息滤波的相位展开方法 | |
CN106707283B (zh) | 基于无味信息滤波的相位展开算法 | |
CN108052978A (zh) | 一种基于支持向量机的ukf相位展开算法 | |
CN116609779A (zh) | 两阶段InSAR多观测处理方法,系统及相关设备 | |
CN114280608A (zh) | 一种DInSAR高程相关大气效应去除方法及系统 | |
CN112835041A (zh) | 结合UKF和AMPM的多基线InSAR高程重建方法 | |
CN110988876B (zh) | 闭合鲁棒双基线InSAR相位解缠方法、系统及可读存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |