一、原理

如图1所示是灰度级为256的相位图,对于相位图中的每一个像素点,它的光强与其相应的相位是有联系的。结合式(1)可知,用于相位计算的反正切值函数的主值位于区间(−π,π](-\pi ,\pi ](−π,π]或[0,2π)[0,2\pi )[0,2π)之间。

图1 相位图[1]

将折叠在[−π,π][-\pi ,\pi ][−π,π]之间的相位值还原成其真实相位值的过程称为相位展开的过程,也称为相位解包裹的过程。由采样定理可知,一个周期内至少存在两个采样点,它们相邻两像素间的真实相位差满足:

ϕ(i,j)\phi \left( i,j \right)ϕ(i,j)为各点的相位值。相位展开的原理是:相邻两像素的相位主值进行比较,如果其相位差超出[−π,π][-\pi ,\pi ][−π,π]范围,则将后面像素对应的相位值加上或减去的整数倍,直到得到的相位是连续的。


其中,ϕ(x,y)\phi \left( x,y \right)ϕ(x,y)是解包后的相位,N是条纹数目且为实数,n是条纹级次,为N的整数部分。

相位图中的任意一点均可以被设置为0从而作为相位解包的参考点,每一点处2π的不确定性可以通过沿着从参考点到当前点的连续路线来计算2π跳跃的次数和公式(2)中的条纹级次n来解决。解包后的相位可以被认为是绝对相位。在图1中,如果选点A作为参考点,B点的绝对条纹级次是15,因为选择不同的参考点会导致其他点处的条纹级次不同,所以应该根据实际应用来选择参考点。

在理想的条件下,相位图中不存在相位不一致(如图2),因此,解包过程可以沿着相位图中的任何一条路线进行。如公式(2)所示,原始相位φ(x,y)\varphi \left( x,y \right)φ(x,y)加上2π的n倍等于解包后的相位ϕ(x,y)\phi \left( x,y \right)ϕ(x,y),其中n为条纹级次。相位解包的关键是确定条纹级次n。通常,条纹级次可以通过公式(3)得到,利用相位图中已解包的像素来计算其相邻像素的n值。

图2 一致相位的解包裹

二、实例分析

在matlab中模拟一个一维信号,如图3(a)所示。并进一步求得其包裹信号,如图3(b),随后,对该包裹相位进行解包裹(解缠)处理,恢复得到的相位,如图3©所示。图3(d)中,将原始信号,与解包裹得到的信号绘制在同一图中,两者完全重合。

图3 相位解包裹实例

上述演示实例matlab程序可私信咨询获取:

参考文献

  1. 马蛟, “基于数字光栅投影的形貌测量技术研究,” 硕士 (合肥工业大学, 2014).

如有问题,可私信交流讨论,博主可提供相关的技术支持。

【相位解缠】一维信号相位解包裹原理及其演示实例相关推荐

  1. 一维信号小波去噪原理及python实现示例

    信号去噪是经常用到的信号预处理过程,以达到在保留原有信号真是信息的基础上尽可能低降低或者消除噪声,获得更高质量的信号,从而为下一步的处理奠定基础. 去噪方法可分为时域方法与频域方法.时域方法是指直接在 ...

  2. matlab 解缠原理,相位解缠算法matlab

    3.2 结合预置滤波的二维UKF相位解缠算法(AUKFPU) 3.2.1 全方位局部估计与二维UKF相位解缠算法[17] 把不敏卡尔曼滤波与传统路径跟踪策略结合起来,利用相位质量...... 龙源期刊 ...

  3. ENVI教程:InSAR技术,相位解缠

    干涉相位只能以 2π为模,所以只要相位变化超过了 2π,就会重新开始和循环.相位解缠是对去平和滤波后的位相进行相位解缠,解决 2π 模糊的问题. (1) SARscape->Interferom ...

  4. 102-RTKLIB中的相位解缠

    这是一篇未解决问题的博文 rtklib手册中的相位解缠公式: 首先ErE_rEr​就很不明白,不知道为什么程序中要那样求解. 另外对于EsE^sEs,大小与前面卫星pco改正求得的各方向单位矢量大小相 ...

  5. 相位 unwrap 与 wrap 算法详解(附代码)

    相位 unwrap 与 wrap 算法详解(附代码) 最近接手了一个项目,光通信方面的,我负责编写初测结果的数据处理算法,其中有一个算法叫做 unwrap 与 wrap,之前没有听说过.通过询问同事与 ...

  6. 详解量子计算:相位反冲与相位反转

    前言 本文需要对量子计算有一定的了解.需要的请翻阅我的量子专栏,这里不再涉及基础知识的科普. 量子相位反冲是什么? 相位反转(phase kickback)是量子计算中的一种现象,通常在量子算法中使用 ...

  7. StaMPS中使用GACOS数据进行大气校正,解缠、计算形变速率

    StaMPS中使用GACOS数据进行大气校正 1 GACOS数据下载 2 查看GACOS用法 3 在StaMPS中使用GACOS进行大气校正 4 stamps6 相位解缠 stamps7时间形变速率 ...

  8. 相位变化对于真实信号的影响是什么?

    一个信号经过傅里叶分解后变成Phase+Magnitude两部分 形象的说, Phase代表信号的形状,Magnitude代表对应形状的大小比例 一维信号的确不太直观,我们来看二维的信号(图像)吧,上 ...

  9. 北斗信号服务器解算,GPS/北斗定位解算算法的研究

    摘要: 卫星导航是一种通过全球卫星导航系统(Global Navigation Satellite System,GNSS)精确的测定地球上任何一点的位置和时间的方法.目前,卫星导航接收机可提供个人定 ...

最新文章

  1. ios图像处理第2部分:核心图形,核心图像,GPUImage
  2. 2.12.ECMAScript--运算符
  3. 《Android插件化技术——原理篇》
  4. latex 波浪线_湖熟镇月牙刀波浪刀带哪家好厂家
  5. 了解css中伪元素 before和after的用法
  6. 火星上首座虚拟房屋以超过50万美元价格售出
  7. 基于JAVA+SpringBoot+Vue+Mybatis+MYSQL的排课管理系统
  8. springboot项目启动rabbitmq报错org.springframework.amqp.AmqpIOException: java.io.IOException
  9. 阻塞IO, 非阻塞IO, 同步IO,异步IO
  10. linux下载TCGA数据,好用的新版TCGA数据下载方法
  11. Java 程序连接 Informix 数据库方法实例介绍
  12. python动态调用自定义模块_python importlib动态导入模块 reload重载模块
  13. Tomcat源码下载地址
  14. python经典书籍推荐-7本经典的Python书籍,你都读过了么?
  15. android wear iphone7,多款 Android Wear 智能手表不兼容iPhone 7
  16. 计算机附件计算器的用法,计算机附件中的计算器使用方法.doc
  17. css html 实现刻度,用CSS划一把刻度尺
  18. LS-DYNA模拟炸药在无限空域中的爆炸
  19. 开源驱动 大黄蜂驱动_客户驱动的开源是软件的未来
  20. 我的京东阅读App书架-211本电子书

热门文章

  1. Java OFD文件转换,OFD转PDF、图片、SVG、HTML工具类(2.0.0版本之前)
  2. Cocos2d-x 3.10 使用BMFont 制作自定义位图字体.fnt文件
  3. linux下安装ffmpeg的教程
  4. 我的团长我的团第二十五集
  5. STM32 网页服务器 LWIP websever ajax实时刷新 (一)
  6. 人工智能中的运筹学与最优化就该这样学之学习路线
  7. 关于springboot持久层框架mybatis+jpa组合使用
  8. 调用百度“车牌识别”API接口(go语言实现)
  9. Matlab深度学习入门实例:基于AlexNet的红绿灯识别(附完整代码)
  10. 【内网安全】域信息收集应用网络凭据CS插件AdfindBloodHound