Skip to main content

主辅影像干涉

影像干涉公式推导

复数影像的两种表现形式

复数影像的数字存储形式

a+bia+bi

复数影像的强度相位表现形式

u=uejϕu=\lvert u\rvert e^{j\phi}

其中, u\lvert u\rvert为强度, ϕ\phi为相位。

两种表达方式的相互关系为:

u=u.real2+u.imag2=a2+b2\lvert u\rvert={u.real}^2+{u.imag}^2=a^2+b^2 ϕ=arctanu.imagu.real=arctanba\phi=\arctan\frac{u.imag}{u.real}=\arctan\frac{b}{a}

两种表现形式的干涉计算

u1=u1ejϕ1=a+biu_1=\lvert u_1\rvert e^{j\phi_1}=a+bi u2=u2ejϕ2=c+diu_2=\lvert u_2\rvert e^{j\phi_2}=c+di

强度相位表现形式

uint=u1u2=u1ejϕ1u2ejϕ2=u1u2ej(ϕ1ϕ2)=u1u2ej(ϕ2ϕ1)\begin{aligned} u_{int} &=u_1\cdot u_2^* \\ &=\lvert u_1\rvert e^{j\phi_1}\cdot \lvert u_2\rvert e^{-j\phi_2} \\ &= \lvert u_1\rvert \lvert u_2\rvert e^{j(\phi_1-\phi_2)} \\ &= \lvert u_1\rvert \lvert u_2\rvert e^{-j(\phi_2-\phi_1)} \end{aligned}

注,*表示共轭

进一步计算,可得:

强度

u1u2=(a2+b2)(c2+d2)=a2c2+b2c2+a2d2+b2d2\begin{aligned} \lvert u_1\rvert \lvert u_2\rvert&=(a^2+b^2)(c^2+d^2) \\ &=a^2c^2+b^2c^2+a^2d^2+b^2d^2 \end{aligned}

相位

ϕ1ϕ2=arctanbaarctandc+cπ=arctanbadc1+bdac+cπ=arctanbcadac+bd+cπ\begin{aligned} \phi_1-\phi_2&=\arctan\frac{b}{a}-\arctan\frac{d}{c}+c\pi \\ &=\arctan\frac{\frac{b}{a}-\frac{d}{c}}{1+\frac{bd}{ac}}+c\pi \\ &=\arctan\frac{bc-ad}{ac+bd}+c\pi \end{aligned}

其中,使用了反正切函数(arctan)相加定律

arctanα±arctanβ=arctanα±β1αβ+cπ\arctan\alpha\pm\arctan\beta=\arctan\frac{\alpha\pm\beta}{1\mp\alpha\beta}+c\pi

数字存储形式

uint=(a+bi)(c+di)=ac+bd+(bcad)i\begin{aligned} u_{int} &=(a+bi){(c+di)}^* \\ &=ac+bd+(bc-ad)i \end{aligned}

强度

uint=(ac+bd)2(bcad)2=a2c2+b2d2+2abcd+b2c2+a2d22abcd=a2c2+b2d2+b2c2+a2d2=u1u2\begin{aligned} \lvert u_{int}\rvert&={(ac+bd)}^2{(bc-ad)}^2 \\ &=a^2c^2+b^2d^2+2abcd+b^2c^2+a^2d^2-2abcd \\ &=a^2c^2+b^2d^2+b^2c^2+a^2d^2 \\ &=\lvert u_1\rvert \lvert u_2\rvert \end{aligned}

相位

ϕint=arctanbcadac+bd=ϕ1ϕ2\begin{aligned} \phi_{int}&=\arctan\frac{bc-ad}{ac+bd} \\ &=\phi_1-\phi_2 \end{aligned}

方位向滤波

距离向滤波