相位展开

某些信号只能以模 2*pi 观察到，这同样也适用于二维和三维图像。在这些情况下，需要进行相位展开以恢复底层的未展开信号。在本例中，我们将演示 [1] 在 skimage 中实现的算法，用于解决此类问题。可以使用 skimage 对一维、二维和三维图像进行展开。在这里，我们将演示二维情况下的相位展开。

import numpy as np
from matplotlib import pyplot as plt
from skimage import data, img_as_float, color, exposure
from skimage.restoration import unwrap_phase


# Load an image as a floating-point grayscale
image = color.rgb2gray(img_as_float(data.chelsea()))
# Scale the image to [0, 4*pi]
image = exposure.rescale_intensity(image, out_range=(0, 4 * np.pi))
# Create a phase-wrapped image in the interval [-pi, pi)
image_wrapped = np.angle(np.exp(1j * image))
# Perform phase unwrapping
image_unwrapped = unwrap_phase(image_wrapped)

fig, ax = plt.subplots(2, 2, sharex=True, sharey=True)
ax1, ax2, ax3, ax4 = ax.ravel()

fig.colorbar(ax1.imshow(image, cmap='gray', vmin=0, vmax=4 * np.pi), ax=ax1)
ax1.set_title('Original')

fig.colorbar(ax2.imshow(image_wrapped, cmap='gray', vmin=-np.pi, vmax=np.pi), ax=ax2)
ax2.set_title('Wrapped phase')

fig.colorbar(ax3.imshow(image_unwrapped, cmap='gray'), ax=ax3)
ax3.set_title('After phase unwrapping')

fig.colorbar(ax4.imshow(image_unwrapped - image, cmap='gray'), ax=ax4)
ax4.set_title('Unwrapped minus original')

Text(0.5, 1.0, 'Unwrapped minus original')

展开过程接受掩蔽数组，并且还可以选择性地假设循环边界以连接图像的边缘。在下面的示例中，我们研究了一个简单的相位斜坡，它通过掩蔽图像的一行被分成两部分。

# Create a simple ramp
image = np.ones((100, 100)) * np.linspace(0, 8 * np.pi, 100).reshape((-1, 1))
# Mask the image to split it in two horizontally
mask = np.zeros_like(image, dtype=bool)
mask[image.shape[0] // 2, :] = True

image_wrapped = np.ma.array(np.angle(np.exp(1j * image)), mask=mask)
# Unwrap image without wrap around
image_unwrapped_no_wrap_around = unwrap_phase(image_wrapped, wrap_around=(False, False))
# Unwrap with wrap around enabled for the 0th dimension
image_unwrapped_wrap_around = unwrap_phase(image_wrapped, wrap_around=(True, False))

fig, ax = plt.subplots(2, 2)
ax1, ax2, ax3, ax4 = ax.ravel()

fig.colorbar(ax1.imshow(np.ma.array(image, mask=mask), cmap='rainbow'), ax=ax1)
ax1.set_title('Original')

fig.colorbar(ax2.imshow(image_wrapped, cmap='rainbow', vmin=-np.pi, vmax=np.pi), ax=ax2)
ax2.set_title('Wrapped phase')

fig.colorbar(ax3.imshow(image_unwrapped_no_wrap_around, cmap='rainbow'), ax=ax3)
ax3.set_title('Unwrapped without wrap_around')

fig.colorbar(ax4.imshow(image_unwrapped_wrap_around, cmap='rainbow'), ax=ax4)
ax4.set_title('Unwrapped with wrap_around')

plt.tight_layout()
plt.show()

在上图中，掩蔽行可以看作是穿过图像的白色线条。底行中两个展开图像之间的差异很明显：在不进行展开（左下）的情况下，掩蔽边界上下的区域根本不交互，导致两个区域之间存在任意整数倍于 2 pi 的偏移。我们也可能将这些区域作为两张单独的图像进行展开。如果为垂直方向启用环绕（右下），情况就会发生变化：展开路径现在允许从图像底部到顶部以及反过来传递，从而提供了一种确定两个区域之间偏移量的方法。

参考文献

[1]

Miguel Arevallilo Herraez、David R. Burton、Michael J. Lalor 和 Munther A. Gdeisat，“基于排序可靠性并在非连续路径后进行的快速二维相位展开算法”，应用光学杂志，第 41 卷，第 35 期，第 7437 页，2002 年