关键词:
电磁逆散射问题
不适定性
非线性
混合正则化
收缩积分方程
无相位电磁逆散射问题
摘要:
电磁逆散射问题的目标是利用测量得到的散射场数据或无相位总场数据重建未知散射体的物理属性。该问题广泛应用于生物医学成像、无损检测、地下勘探、微波遥感等领域。然而,逆散射问题存在病态性和非线性两大挑战,现有的一些成像方法面临着如成像精度不高以及无法重构电大尺寸或者高对比度的散射体等问题。本文研究如何有效地缓解逆散射问题中的不适定性、降低散射模型的非线性,并且探讨如何利用无相位总场数据建立优化算法。通过解决这些关键问题,以提高逆散射问题的成像精度,并能更好地重构电大尺寸或高对比度的散射体。
首先,对于二维逆散射问题,提出了一个新的混合正则化方法(非凸正则化和修正的傅里叶基扩展正则化)用于缓解不适定性。受非凸正则化(nonconvex regularization,NR)在许多领域成功应用的启发,尝试将其应用于逆散射问题的建模中。具体而言,在重建任务中引入描述散射体稀疏性的先验信息,旨在提高重建结果的准确性和质量。进一步通过分析矩阵算子奇异向量的离散傅里叶变换系数,提出了一个修正的傅里叶基扩展(modified Fourier bases-expansion,MFBE)正则化,其引入了新形式的傅里叶系数张量。然后,将MFBE正则化和NR结合,构建了一种新的混合正则化方法。基于有效降低逆散射问题非线性的收缩积分方程(contraction integral equation,CIE)模型,我们展示了如何在建模中利用混合正则化方法以获得稀疏解,以及如何求解新的损失函数。对合成和实测数据进行测试,实验结果表明该算法能有效地反演电大尺寸或高对比度散射体,并显著提高了成像精度。特别是在处理无损散射体时,在未使用无损先验信息的情况下,反演结果的虚部全为零。
其次,对于三维电磁逆散射问题,提出了不同于对比源反演型算法的邻近交替方向乘子法:将非凸L1/2正则化和傅里叶基扩展正则化相结合以求解三维逆散射问题。傅里叶基扩展(Fourier bases-expansion,FBE)正则化直接应用于建模,L1/2正则化应用于未知量,以获得更稀疏的解和更好的反演效果。此外,在反演建模过程中使用了权重调整嵌套方案(weight adjustment nested scheme,WA-NS)。基于CIE模型、WA-NS和上述混合正则化技术,构造了新的重建模型。由于该函数是非凸、非光滑、非Lipschitz的,很难直接求解。因此,我们设计了一个有效的邻近交替方向乘子法(proximal alternating direction method of multipliers,PADMM)以求解相应优化问题,所有子问题都很容易求解或者有闭式解。合成数据和Fresnel实测数据的反演实验表明所提出的算法可以有效去除伪影,恢复散射体的轮廓。
最后,探讨了利用全变分(total variation,TV)和非凸L1/2正则化求解无相位电磁逆散射问题。使用有效降低逆散射问题非线性的CIE模型,建立了一个新的非凸、非光滑、非Lipschitz的损失函数。由于模型包含非凸L1/2正则化和TV正则化,能够实现对散射体的稀疏表示和平滑化,从而提高问题的可行性和计算效率。所建立的模型难以直接求解,因此通过添加等式约束分离目标函数中的变量,然后设计交替方向乘子法(alternating direction method of multipliers,ADMM)以最小化其等价的增广拉格朗日问题。该研究分别对二维和三维情况下无相位逆散射问题进行了数值实验,结果显示新算法的反演结果具有更好的边界保留效果和成像质量。