首页 > 分享 > 相位解包裹（五）枝切法（Goldstein’s branch cut algorithm）

相位解包裹（五）枝切法（Goldstein’s branch cut algorithm）

花匠小妙招
2024-12-18 06:17

上一篇文章讲了相位解包裹算法的分类，那从这一篇文章开始，就可以开始介绍具体的相位解包裹算法了。

因为介绍过了残点的运算，所以这里首先介绍利用残点的特性去进行相位解包裹算法的一种经典算法，枝切法（Goldstein’s branch cut algorithm）。枝切法的解包裹是和路径有关的，所以其属于空间相位解包裹中的一种。

残点的极性（polarity/charge）

从残点的计算我们得知，在一个2*2的小区域里面求相位差的环路积分，结果可能为0，±1（所有结果都除2pi），+1成为正残点（positive），-1成为负残点（negative），这一属性成为残点的极性（polarity）或者charge（实在不知道怎么翻译）

可以得到，相位差的环路积分，实际就是经过多少个正负残点charge的和的2pi倍

oint{triangledown phi left( r right)dr=2pi times left( text{sum of enclosed phase residue charges} right)}

之前一直说，相位解包裹因为残点存在，是一个与路径有关的过程，那怎么可以让他变成一个与路径无关的过程呢？

从上面的式子可以看出来，要是相位差的环路积分等于0，就需要让charge为0，那怎么才能为0呢，就是即路径经过的区域里，残点的charge被平衡（balance）了，意思是区域中正负残点的数量相等。只要所有区域的正负残点都是平衡的，那相位解包裹的过程自然就和路径没有关系了。

空间相位解包裹的基础就是残点的理论，通过设置一些障碍，例如连接残点的branch，让解包裹的路径绕开这些障碍，解包裹就能顺利进行。

以上讨论都是正常情况，假如认为假如噪声，导致计算的残点不等于±1，这些特殊情况需要另外讨论。

枝切法

枝切法就是利用连接残点的branch，让解包裹路径躲开这些区域的相位解包裹的算法。其最重要的贡献还是如何连接残点。

算法的核心：通过枝节（branch）连接距离最近的正负残点，并使它们达到平衡（balance）

残点的连接

对于用branch连接残点是有讲究的，如下图（a）从A到B解包裹，有3条路径，明显中间的路径经过了残点的区域，自然得到的结果和另外两条路径会不一样。假如像图（b）一样把正负残点连接使它们平衡，解包裹的路径需要躲开branch，只能走另外两条路的其中一条，不管哪一条，它们的结果都一样。但如果像图（c）一样，通过吧残点和相位图边界连接起来以平衡它们，那么解包裹路径还是会选择中间错误的路。

所以如何通过branch连接正负残点，是需要一定的法则，如下图，连接残点的branch应该尽可能短。Goldstein等提出的方法，也只是能解决部分残点连接的问题，不是一个完美的方案，在他以后还有很多研究者提出各种各样的方法，要根据情况来决定算法的选取。这里也只讨论Goldstein枝切法一种。

Goldstein枝切法残点连接

1、首先当然是需要扫描检测出所有的残点。

2、将一个3*3的搜索窗放在残点的位置，寻找别的残点。如果找到了，将它们用branch连接起来。

（1）如果另一个残点的极性与当前残点相反，那么就认为这一个branch是已经平衡了，打上“balanced”标记。

（2）如果另一个残点的极性和当前残点相同，那搜索窗就放在新检测到的这个残点的位置，继续寻找残点，直到找到极性相反的残点并将它们用branch连接起来且他们charge的总和为0（也就是平衡了balanced）；或者在搜索窗中再无残点被检测到，如果是后面这种情况，那么搜索窗会增大2，也就是变成5*5的窗，继续从开始的那个残点开始检测别的残点。

肥肠重要，一定要注意：

1、当在搜索窗中检测到了残点，都要将它们用branch连接起来。如果检测到的这个残点之前还未被用branch连接过，那它的极性需要加到charge里面来，如果连接过了的就不需要加。（目的就是要用branch上残点的charge总和为0，使它平衡）

2、如果在搜索窗遇到了相位图的边界，那直接把branch和边界连接，同时认定这一branch的charge总和为0。

3、之前讲过，残点实际是一个2*2的小区域，我们标记这个小区域左上角的那个像素为残点，但在残点的连接中，不需要连接被标记的两个像素，而是选择两个小区域最短的路径将他们连接起来，就像下面的图。