一种基于光强传输方程的散射成像相位恢复仿真研究

一、前言

当光线经过云雾、烟尘等散射介质时,会引起光子的随机传输问题,导致光波无法聚焦成像,只能得到由散射光产生的散斑图像,散射成像技术就是从杂乱无章的散斑图样中恢复出目标信息的成像技术。目前散射成像技术主要有:利用光学相位共轭思想的波前调制技术,通过反馈调制算法控制空间光调制器等设备对入射光波前进行相位补偿,达到透过散射介质重新聚焦目的;利用测量散射介质的传输矩阵,对散斑场进行解码,实现透过散射介质成像;利用光学记忆效应(Optical Memory Effect,OME),根据单、多帧散斑场图像的相关特性,结合相位恢复算法实现成像[1]。

二、基于光强传输方程的散斑相位恢复仿真实例

2.1 基于双平面的光强轴向微分估计

光强传输方程法需要采集不同离焦面上的光强信息,可通过移动物平面或者移动图像平面(相机)获得物体离焦面上的光强信息。光路结构如图所示,由于加入了散射介质(如毛玻璃,散射片等),因而相机接收到的是散斑图像。

首先,仿真得到焦面处一幅散斑图,并通过角谱衍射得到离焦距离为1微米的离焦散斑图。通过光强传输方程,即可通过两幅散斑图像恢复出被测物体的真实相位,结果如图所示。

在这里插入图片描述
求解光强传输方程需要预先获得光强的轴向微分,该量不可直接测量,需要通过数值有限差分得到。Teague于1983年提出通过采集两幅轻微离焦图像,使两幅图像相对于中心聚焦图像的离焦距离相等且方向相反,利用中心有限差分法估计获得光强的轴向微分,即[2],
在这里插入图片描述如下图所示,分别获取相对于中心聚焦图像的离焦距离相等且方向相反散斑图,通过光强传输方程,恢复的相位如图所示。
在这里插入图片描述

2.2 基于多平面的光强轴向微分估计

考虑到双平面的轴向微分估计方法中仅采用两个平面的光强信息,其中唯一可控的参数为离焦距离。为了解决上述问题,许多研究人员提出采用多面(多于两个平面)强度测量估计轴向微分,更灵活地校正非线性误差或者降低噪声的影响。给定2n+1个测量平面的光强,基于多平面的光强轴向微分的差分公式均可表示为[2]

在这里插入图片描述
如下图所示,在仿真中,分别获取多幅离焦散斑图,通过光强传输方程,恢复的相位如图所示。

在这里插入图片描述

参考文献

[1] 邵晓鹏, 苏云, 刘金鹏, et al. 计算成像内涵与体系(特邀) [J]. 光子学报, 2021, 50(05): 9-31.
[2] 左超, 陈钱, 孙佳嵩, et al. 基于光强传输方程的非干涉相位恢复与定量相位显微成像:文献综述与最新进展 [J]. 中国激光, 2016, v.43;No.474(06): 227-57.


⭐️⭐️⭐️⭐️⭐️ ···· **博主简介** ··· ⭐️⭐️⭐️⭐️⭐️

博士研究生生 ,研究方向主要涉及定量相位成像领域,具体包括干涉相位成像技术(如**全息干涉**、散斑干涉等)、非干涉法相位成像技术(如波前传感技术相位恢复技术)、此外,还对各种相位解包裹算法相干噪声去除算法等开展过深入的研究。

  • 4
    点赞
  • 24
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 1
    评论
### 回答1: 好的,以下是一个简单的 MATLAB 程序,实现了利用光强传输方程恢复相位的功能: ``` % 定义光强传输方程(理想情况下) function I = intensity_transfer_equation(phase) I = abs(exp(1i * phase)); end % 定义一个模拟的相位场 true_phase = linspace(0, 2 * pi, 100); true_intensity = intensity_transfer_equation(true_phase); % 加入一些噪声 observed_intensity = true_intensity + 0.1 * randn(size(true_intensity)); % 利用光强传输方程恢复相位 recovered_phase = angle(observed_intensity); % 画出恢复相位与真实的相位的对比图 plot(true_phase, true_intensity, 'b', ... true_phase, observed_intensity, 'ro', ... true_phase, intensity_transfer_equation(recovered_phase), 'g--'); legend('真实光强', '观察到的光强', '恢复光强'); xlabel('相位'); ylabel('光强'); ``` 这段程序模拟了一个理想的光强传输方程,并且生成了一个模拟的相位场。接着,它加入了一些随机噪声,模拟了相位场的观察。最后,它利用光强传输方程恢复相位场,并且画出了真实的相位场、观察到的光强以及恢复光强的对比图。 ### 回答2: 利用光强传输方程恢复相位的MATLAB程序如下: ```matlab % 假设原始图像为I,相位信息为phi,光传播距离为d I = imread('原始图像.jpg'); phi = 2 * pi * rand(size(I)); % 随机生成相位信息 d = 10; % 光传播距离,单位为像素 % 进行光传输 J = I .* exp(1i * phi * d); % 恢复相位 phi_recover = angle(J) / d; % 显示结果 subplot(1, 2, 1); imshow(I); title('原始图像'); subplot(1, 2, 2); imshow(phi_recover, []); title('恢复相位'); % 保存结果 imwrite(phi_recover, '恢复后的相位.jpg'); ``` 在这个程序中,我们假设原始图像为`I`,使用`rand`函数生成一个与原始图像相同大小的随机相位信息`phi`。然后,通过将原始图像乘以一个复数指数以及相位信息,实现光传输。光传输的距离由变量`d`表示。接下来,我们恢复相位信息,通过求解复数图像`J`的角度,并将结果除以光传播距离。最后,我们将原始图像恢复相位信息进行显示,并将恢复相位信息保存为图片。 请注意,在实际应用中,要使用光强传输方程进行相位恢复,可能需要更复杂的处理步骤和采样策略。此处提供的代码仅是一个简单的示例,需要根据具体应用进行适当的修改和调整。 ### 回答3: 可以通过以下步骤编写一个利用光强传输方程恢复相位的MATLAB程序: 1. 首先,导入图像。使用MATLAB的imread函数读取一个图像,将其转换为灰度图像。 2. 然后,为图像添加一个已知的相位畸变。可以通过创建一个二维网格矩阵,根据所需的相位畸变模式在每个像素位置上添加一个相位值来实现。 3. 接下来,计算输入图像的傅里叶变换。使用MATLAB的fft2函数对灰度图像进行二维快速傅里叶变换,得到输入图像的频谱。 4. 然后,通过光强传输方程计算频谱的相位。利用光强传输方程,将频谱的光强和已知的相位畸变相乘,得到修正后的频谱。 5. 接下来,计算修正后频谱的逆傅里叶变换。使用MATLAB的ifft2函数对修正后的频谱进行二维逆傅里叶变换,得到修正后的图像。 6. 最后,显示原始图像和修正后的图像。使用MATLAB的imshow函数分别显示原始图像和修正后的图像,以进行比较。 整个MATLAB程序的代码如下: ```MATLAB % 步骤1:导入图像 img = imread('image.jpg'); gray_img = rgb2gray(img); % 步骤2:添加相位畸变 [rows, cols] = size(gray_img); [X, Y] = meshgrid(1:cols, 1:rows); phase_distortion = sin(X/10) + sin(Y/10); % 根据需要调整相位畸变模式 % 步骤3:计算输入图像的傅里叶变换 spectrum = fft2(double(gray_img)); % 步骤4:计算频谱的相位 corrected_spectrum = abs(spectrum) .* exp(1i * phase_distortion); % 步骤5:计算修正后频谱的逆傅里叶变换 corrected_img = abs(ifft2(corrected_spectrum)); % 步骤6:显示原始图像和修正后图像 subplot(1,2,1); imshow(gray_img); title('原始图像'); subplot(1,2,2); imshow(corrected_img, []); title('修正后的图像'); ``` 请注意,上述程序只是一个简单的示例,具体的相位恢复方法可能因应用而异。您可以根据自己的需求以及图像的特点进行相应的改进。
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

简单光学

您的鼓励是我创作的最大动力!

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值