发明名称 基于构造边的高精度快速相位解缠方法
摘要 本发明公开了一种基于构造边的高精度快速相位解缠方法。包括以下几个步骤:利用三步相移法投射三幅光栅图像,用相机分别捕获,得到由包裹相位图;计算产生质量图;在质量图上,由相邻两像素质量值构造水平与竖直权值边;初始化各个像素点的状态;删除冗余的权值边;进行像素点的相位解缠。本发明通过优化权值边的数量,降低了数据处理的规模,节省了解缠时间。本发明的相位解缠方法具有精度高、效率高、抗干扰能力强的特点。
申请公布号 CN103968781B 申请公布日期 2016.09.14
申请号 CN201410216216.1 申请日期 2014.05.21
申请人 哈尔滨工程大学 发明人 陆军;李积江;王成成;苏航;蔡成涛;夏桂华;朱齐丹;韩吉瑞
分类号 G01B11/25(2006.01)I 主分类号 G01B11/25(2006.01)I
代理机构 代理人
主权项 基于构造边的高精度快速相位解缠方法,其特征在于:包括以下几个步骤,步骤一:利用三步相移法投射三幅光栅图像,用相机分别捕获,得到包裹相位图;步骤二:由包裹相位图计算产生质量图;步骤三:在质量图上,由相邻两像素点质量值构造水平与竖直权值边;步骤四:初始化各个像素点的状态;步骤五:删除冗余的权值边;步骤六:进行像素点的相位解缠;解缠过程为:第一步:选择具有最小权值的权值边关联的两像素点,比较两个像素点的质量值,质量值大的像素点相对质量值小的像素点解缠,并将质量值大的像素点设置为已解缠标志,将两点加入到一个集合中,记录该集合的秩;第二步:将剩余的权值边按其权值从小到大的排序作为路径向导进行解缠,根据权值边关联两像素点的状态,分三种情况进行解缠:a)有一像素点已解缠而另一像素点未解缠,解缠两像素点中未解缠的像素点,合并到已解缠像素点集合中;b)两像素点均未解缠,解缠这两像素点,并标记为一个新集合;c)两像素点均已解缠且分别属于不同集合,比较两个集合的秩的大小得到小秩集合和大秩集合,属于小秩集合中的像素点相对大秩集合的像素点解缠,合并到大秩集合中;解缠的方法为:<maths num="0001"><math><![CDATA[<mrow><mi>&phi;</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>m</mi></msub><mo>,</mo><msub><mi>y</mi><mi>n</mi></msub><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>&phi;</mi><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>y</mi><mi>n</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>&Delta;&phi;</mi><mi>r</mi></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>m</mi></msub><mo>,</mo><msub><mi>y</mi><mi>n</mi></msub><mo>)</mo></mrow><mo>+</mo><mn>2</mn><mi>&pi;</mi></mrow></mtd><mtd><mrow><mi>i</mi><mi>f</mi></mrow></mtd><mtd><mrow><msub><mi>&Delta;&phi;</mi><mi>r</mi></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>m</mi></msub><mo>,</mo><msub><mi>y</mi><mi>n</mi></msub><mo>)</mo></mrow><mo>&le;</mo><mo>-</mo><mi>&pi;</mi></mrow></mtd></mtr><mtr><mtd><mrow><mi>&phi;</mi><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>y</mi><mi>n</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>&Delta;&phi;</mi><mi>r</mi></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>m</mi></msub><mo>,</mo><msub><mi>y</mi><mi>n</mi></msub><mo>)</mo></mrow></mrow></mtd><mtd><mrow><mi>i</mi><mi>f</mi></mrow></mtd><mtd><mrow><mo>-</mo><mi>&pi;</mi><mo>&lt;</mo><msub><mi>&Delta;&phi;</mi><mi>r</mi></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>m</mi></msub><mo>,</mo><msub><mi>y</mi><mi>n</mi></msub><mo>)</mo></mrow><mo>&lt;</mo><mi>&pi;</mi></mrow></mtd></mtr><mtr><mtd><mrow><mi>&phi;</mi><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>y</mi><mi>n</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>&Delta;&phi;</mi><mi>r</mi></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>m</mi></msub><mo>,</mo><msub><mi>y</mi><mi>n</mi></msub><mo>)</mo></mrow><mo>-</mo><mn>2</mn><mi>&pi;</mi></mrow></mtd><mtd><mrow><mi>i</mi><mi>f</mi></mrow></mtd><mtd><mrow><msub><mi>&Delta;&phi;</mi><mi>r</mi></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>m</mi></msub><mo>,</mo><msub><mi>y</mi><mi>n</mi></msub><mo>)</mo></mrow><mo>&GreaterEqual;</mo><mi>&pi;</mi></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000977593780000011.GIF" wi="1414" he="223" /></maths>其中,φ(x<sub>m</sub>,y<sub>n</sub>)为待解缠像素点的绝对相位值,φ(x<sub>m‑1</sub>,y<sub>n</sub>)为相邻已解缠的像素点的绝对相位值,Δφ<sub>r</sub>(x<sub>m</sub>,y<sub>n</sub>)=φ<sub>r</sub>(x<sub>m</sub>,y<sub>n</sub>)‑φ<sub>r</sub>(x<sub>m‑1</sub>,y<sub>n</sub>)是相邻两像素点包裹相位值的差值,φ<sub>r</sub>(x<sub>m</sub>,y<sub>n</sub>)是像素点(x<sub>m</sub>,y<sub>n</sub>)的包裹相位值,φ<sub>r</sub>(x<sub>m‑1</sub>,y<sub>n</sub>)是像素点(x<sub>m‑1</sub>,y<sub>n</sub>)的包裹相位值。
地址 150001 黑龙江省哈尔滨市南岗区南通大街145号哈尔滨工程大学科技处知识产权办公室