简介资源为一份关于基于光学相位共轭的数字化波前整形技术DOPC的完整技术文档主要面向生物光子学、光学成像和生物医学工程等领域的科研人员与高年级学生帮助理解如何利用波前整形克服生物组织散射、突破传统光学约1 mm的聚焦深度限制。文档先对比基于反馈、传输矩阵测量和光学相位共轭三条技术路线的优劣随后着重剖析DOPC系统中的相位共轭镜构成以及光场测量和光场调控两阶段工作机理并包含四步相移全息、相机与空间光调制器像素匹配、焦点PBR提升等实现细节可作为该方向的导览性综述或课题入门笔记。资源仅此1份docx文档压缩包大小约212KB便于快速下载阅读。目前已有123人学习下载。1. 光子无法穿透的 1 mm为什么需要数字化波前整形生物组织对光并不友好。折射率的不均匀分布会把一束准直光打散成散斑传统光学成像在组织里的工作深度被压在 1 mm 附近限制因素不是探测器灵敏度而是弹道光在出射前就已经消失殆尽。早年学界默认散射是随机破坏过程直到波前整形技术出现大家才意识到组织内的散射事件其实是确定性的光走了哪些路径、每条路径累积了多少相位延迟都可以被测量和反推。基于光学相位共轭的数字化波前整形DOPC是三条技术路线中效率最高的一种它用相位共轭镜对散射光场做时间反演一次测量就能拿到预补偿所需的最优相位图无需像反馈迭代法那样反复逼近因此在血流、心跳导致散射快速变化的活体场景中最有竞争力。本文按测量—调控—对准—引导星—速度这条链路拆开讲最后给出高速化路线的工程判断。2. 四步相移全息与相位共轭镜的构建2.1 两种相位共轭镜的取舍相位共轭镜负责把发散散射光的波前“翻过来”再送回散射介质。早期实现依赖非线性晶体的四波混频响应速度极快但产生的共轭波能量很低生物组织内聚焦需要足够光强来激发荧光或驱动光遗传晶体方案很难满足。DOPC 的替代思路是把相位共轭拆成两个数字化步骤相机记录散射光场空间光调制器SLM重新发一个共轭波前。SLM 上的每个像素可独立调制百万像素的模式数意味着对散射波前的刻画能做到非常细调制效率也远高于晶体方案。代价是系统复杂度上移——相机传感器面和 SLM 面必须严格光学共轭这一步做不好后续所有补偿都白费。2.2 光场测量四步相移恢复散射相位SLM 只能加载相位相机只能记录强度所以相位信息必须先编码到强度里才能被测量。做法是在散射光路里引入一束平面参考光与散射光干涉形成散斑然后对参考光依次附加 0、π/2、π、3π/2 四步相移记录四帧干涉图。第 k 帧强度分布为[ I_k A_R^2 A_S^2 2 A_R A_S \cos(\phi_S - \phi_R k\pi/2) \quad k0,1,2,3 ]式中 (A_R)、(\phi_R) 是参考光振幅和相位(A_S)、(\phi_S) 是散射光的振幅和相位。参考光相位在空间上近似均匀因此用两对差分帧就能把散射光相位从余弦项里解出来。import numpy as np def four_step_phase_shift(I0, I1, I2, I3): 四步相移法恢复散射光场相位分布。 I0..I3 : 参考光相位为 0, pi/2, pi, 3pi/2 时相机记录的干涉强度图 返回 : 相位图 phi, 弧度制, 范围 [-pi, pi] # 差分对消掉 A_R^2 A_S^2 直流背景 cos_term I0 - I2 # 2*A_R*A_S*cos(phi_S - phi_R) sin_term I3 - I1 # 2*A_R*A_S*sin(phi_S - phi_R) phi np.arctan2(sin_term, cos_term) # 取相位, 振幅项被自动约去 return phi这段代码的要点在arctan2它同时利用 sin 和 cos 两个正交分量相位输出范围是完整的 ([-\pi,\pi])不会出现单用 arccos 时正负象限混叠的问题。两帧差分把直流背景和 SLM 反射率不均匀性一并消掉这是四步相移比单帧离轴全息抗噪能力更强的原因。实际采集时四帧图必须在散斑相关时间内拍完否则相位信息已经被生物运动抹平这一点在高速系统里会成为主要瓶颈。得到散射相位分布 (\phi_S(r)) 后SLM 需要加载的是 (-\phi_S(r))也就是对每个像素做共轭。2.3 光场调控为什么共轭相位能让光“倒流”把 (-\phi_S(r)) 加载到 SLM 上让一束与参考光同源的平面波经过 SLM 反射等于给平面波逐点叠加一个与散射过程相反的相位。依据波动方程在时间反演下的对称性这个相位共轭波重新进入散射介质后会沿着散射路径原路返回在目标位置相干叠加成焦点。注意这里不是让光真的倒流而是用事先测得的相位延迟做预补偿使所有路径的光在目标位置重新同相。这个过程只需要一次测量加一次加载对比反馈迭代法动辄数百轮优化速度优势明显代价是测量阶段必须保证参考光的相位稳定性。2.4 三种波前整形路线的横向对比维度反馈迭代法传输矩阵法DOPC 相位共轭法最优相位图获取迭代优化需数百至数万次光强采样一次测量复传输矩阵矩阵规模随像素数平方增长单次测量散射光场一次共轭加载系统复杂度较低SLM 加反馈探测器即可较高参考光路需长期稳定中高相机与 SLM 需严格像素匹配对动态散射介质适应性差迭代期间散斑已变化差矩阵需频繁重测好单帧测量后可快速加载典型场景静态散射体、教学验证光纤模式控制、光通信生物活体聚焦、时变组织三条路线最终都在求解同一个问题预补偿相位图。区别只在于相位图从哪来、拿得多快。3. 像素级精度的系统校准从 PBR 61 到 52000 的对准工程3.1 PBR 的上限由 SLM 像素数决定焦点峰值与背景功率之比PBR是衡量聚焦质量的核心指标。若 SLM 有效像素数为 (N)理论上限是 (\pi N/4)一百万像素对应约 78 万 PBR。实际 DOPC 系统长期停留在数十到数百 PBR问题出在相机传感器面与 SLM 面之间存在 6 个自由度的对准误差三个平动 (\Delta x,\Delta y,\Delta z) 和三个转动 (\Delta\theta_x,\Delta\theta_y,\Delta\theta_z)。任意一个误差都会把单个 SLM 像素的共轭相位送到错误的位置散斑统计叠加后焦点被抹进背景里。3.2 物理掩模时代的粗对准2010 年第一套 DOPC 系统用物理掩模辅助像素匹配先用一相机标定掩模位置再对照 SLM 上加载的掩模图样逐像素找对应关系。这个方案能工作但光路复杂双相机切换引入额外像差对准精度完全依赖手动调节最终 PBR 约 600离理论值差三个数量级。这类方法的问题本质是“用实物标定去补偿光学误差”标定本身的误差无法自举消除。3.3 角谱传播补偿与 6 自由度搜索2014 年 Jang 的方案把问题从硬件搬到软件里先做一次粗略的机械对准再用角谱分析法在 SLM 上加载补偿相位把相机平面与 SLM 平面之间残余的传播误差数值化地补掉。import numpy as np from numpy.fft import fft2, ifft2, fftshift, ifftshift def angular_spectrum_propagate(field, dz, wavelength, pixel_size): 角谱法计算复振幅场自由传播 dz 距离后的结果。 field : 源平面复振幅 (对应 SLM 平面) dz : 传播距离, 正数表示沿光轴前进 wavelength : 工作波长 pixel_size : 像素物理尺寸, 与 wavelength 单位一致 ny, nx field.shape fx fftshift(np.fft.fftfreq(nx, dpixel_size)) fy fftshift(np.fft.fftfreq(ny, dpixel_size)) FX, FY np.meshgrid(fx, fy) # 衍射传递函数 exp(i * kz * dz), kz 由色散关系决定 kz_sq (1.0 / wavelength) ** 2 - FX**2 - FY**2 H np.exp(1j * 2 * np.pi * dz * np.sqrt(np.maximum(kz_sq, 0))) H np.where(kz_sq 0, H, 0) # 倏逝波分量直接置零 return ifft2(fft2(field) * H)这个函数解决的是 (\Delta z) 和部分离焦误差但 (\Delta x,\Delta y) 对应像素平移(\Delta\theta_x,\Delta\theta_y,\Delta\theta_z) 对应平面旋转需要在外层包装坐标变换。常见的做法是把 6 个参数依次扫描先用粗步长扫 (\Delta x,\Delta y) 找到焦点的大致位置然后扫 (\Delta z) 确定离焦量最后旋转三个角度做精细搜索每一步都以焦点强度作为反馈值保留最优参数。整个优化过程耗时约 10 分钟用贪婪策略从粗到细推进最终把 PBR 从粗对准的 61 提升到 52000。3.4 泽尼克多项式修正与 SLM 曲率补偿角谱方案假设系统像差是纯粹的刚体偏差实际 SLM 自身有表面曲率成像透镜也有高阶像差。2016 年 Azimipour 引入自适应光学里成熟的爬山算法在 SLM 上叠加矩形泽尼克多项式作为修正面形以焦点亮度为反馈逐阶扫描多项式系数修正转动误差和系统像差后 PBR 达到 121000。2018 年 Hemphill 更进一步用四步相移全息直接测量 SLM 表面曲率把曲率相位预先加载掉再用干涉法定标 SLM 上生成的十字符号完成空间映射在 5 mm 厚散射样品上拿到近一万 PBR达到理论极限的 14.6%。方案年份对准/补偿手段实现 PBRCui 物理掩模2010掩模 双相机标定~600Jang 角谱自动校准2014角谱传播补偿 6 参数贪婪搜索52000Azimipour 泽尼克修正2016爬山算法 泽尼克多项式系数扫描121000Hemphill 综合校正2018SLM 曲率预补偿 十字符干涉标定~10000提示角谱补偿代码里dz的单位必须和pixel_size一致否则相位因子差几个数量级。先验证空场传播一个已知距离的菲涅尔衍射结果再接入 DOPC 搜索流程。从 61 到 121000 的跨越说明 DOPC 系统到达理论极限之前校准精度是第一约束。即便到了 2018 年也还有 85% 的 PBR 空间没有榨干误差来自校准残余和 SLM 相位灰阶量化的固有损失。4. 超声引导星与 TRUE 聚焦外差干涉的时频分析4.1 组织内聚焦需要“坐标”前面讨论的 DOPC 解决的是穿透散射介质后的聚焦真正有生物医学价值的是在散射组织内部聚焦。问题变得棘手光在组织内部没有清晰的参考点必须引入一个“引导星”来标记目标位置。侵入式引导星包括荧光分子、二次谐波纳米粒子、磁性粒子和基因编码蛋白对比度高、尺寸小但都要向组织内注入异物安全性受限。内源性散射体如血流中的红细胞虽然无侵入但成片出现无法单独操控。理想的引导星要同时满足非侵入、尺寸可控、可自由扫描三个条件。4.2 AOM 频移与声光标记聚焦超声天然满足这三个条件。超声在散射组织内受散射影响远小于光可以通过换能器精确聚焦到目标位置。当散射光穿过超声焦点时声光效应会给一小部分光子叠加一个频移这部分光子被称为标记光其余散射光频率不变。关键在于分离标记光与背景散射光用声光调制器AOM工作频率 40~60 MHz把参考光频率搬到标记光附近让参考光与标记光外差干涉形成拍频而未被标记的散射光因为与参考光频率差太大拍频超出相机帧率表现为稳定背景。相机记录到的强度为[ I(t) A_R^2 A_T^2 A_S^2 2A_R A_T \cos(2\pi \Delta f t \phi_T - \phi_R) ]其中 (\Delta f) 是参考光与标记光的频差(\phi_T) 正是超声焦点处散射光携带的相位提取出它相位共轭就能把光送回声学焦点。这套方案被称为 TRUETime-reversed ultrasonically encoded optical focusing。4.3 外差拍频信号的锁相提取四步相移在低频静态场景里够用但拍频信号的相位随时间变化需要锁相检测来同时提取振幅和相位。锁相 CMOS 相机每个像素内置高速光电转换和移相电路并行输出信号的振幅和相位数值等价于下面这段正交解调import numpy as np def lockin_demodulate(signal, tvec, f_beat): 正交锁相解调: 从外差干涉时间序列中提取振幅与相位。 signal : 相机单像素记录的拍频强度序列 tvec : 对应的时间序列 (秒) f_beat : 参考光与标记光的频差 (Hz), 由 AOM 驱动频率决定 # 本地振荡器: 同频余弦和正弦, 系数 2 用来抵消解调增益 ref_cos 2.0 * np.cos(2 * np.pi * f_beat * tvec) ref_sin 2.0 * np.sin(2 * np.pi * f_beat * tvec) # 混频得到 I/Q 分量, 再用滑动平均做低通滤掉 2*f_beat 分量 win int(np.round(len(tvec) / (f_beat * (tvec[-1] - tvec[0])))) # 约一个拍频周期 kernel np.ones(win) / win I np.convolve(signal * ref_cos, kernel, modesame) Q np.convolve(signal * ref_sin, kernel, modesame) amp np.sqrt(I**2 Q**2) phase np.arctan2(Q, I) return amp, phase说明几点win取一个拍频周期做平均窗口低于这个窗口滤不干净高频项高于这个窗口会抹掉超声焦点处散射光的快速变化ref_cos和ref_sin必须严格正交频率偏差会导致相位随测量时间线性漂移实验上通常用与 AOM 同一时钟源驱动的信号发生器来保证相干性。锁相相机的 0.3 ms 测量速度就是这么来的——每个像素的模拟电路并行完成解调不需要把四帧强度图传回计算机再算。代价是像素规模受限早期锁相 CMOS 只有 300×300 像素。4.4 引导星性质决定焦点质量引导星的位置和尺寸直接决定焦点质量。超声的波长比光长几个数量级声焦点的体积远大于光学衍射极限因此 TRUE 单次聚焦的焦点偏大。迭代 TRUE 可以逐轮收缩焦点第一轮聚焦到声焦点第二轮以第一轮焦点为新的引导源继续优化实验上观察到 2~3 倍分辨率和约 20 倍 PBR 的提升理论上可以收敛到单个散斑大小。2013 年 Judkewitz 用散斑涨落编码做引导星通过数千次光场测量在组织内实现了约 5 μm 的单散斑聚焦代价是耗时数小时。2015 年 Ruan 用聚焦超声击破血管内微气泡来提高信号对比度焦点亮度和尺寸都有明显改善。引导星类型侵入性尺寸控制生物安全性可操控性荧光分子 / 纳米粒侵入亚微米需注射存在代谢负担只能随注入位置固定磁性粒子侵入微米级有安全争议可用磁场粗略引导内源散射体动态非侵入成片出现高几乎不可操控聚焦超声非侵入由声透镜焦深决定高可自由扫描5. 散斑相关时间约束下的高速 DOPC 路线5.1 动态散射介质的时间窗活体组织的散斑相关时间在数十微秒到数毫秒量级主要贡献来自流动红细胞。波前整形的完整流程必须在一个相关时间内完成测量、传输、计算、加载相位、聚焦。传统四步相移加 SLM 的流程需要几十毫秒散斑早就变了共轭相位打在已经失效的散射状态上焦点自然出不来。5.2 四条高速化路线对比方案测量方式调控器件系统响应主要代价DMD FPGAWang 2015单帧二值化相位DMD刷新 22 kHz5.3 ms调制效率降至 (1/2\pi)不能兼容超声引导星锁相相机Liu 2016并行锁相解调0.3 ms未完成调控测量级 0.3 ms像素仅 300×300传输速率低离轴全息Hemphill 2017单帧强度图提取光场普通 SLM约 9 ms牺牲像素分辨率计算繁琐铁电液晶 SLMLiu 2017两帧曝光二值化相位铁电液晶毫秒级响应约 6 ms完整 TRUE 聚焦相位调制能力有限器件成本高5.3 验证方法用系统响应时间做可行性判据判断一套 DOPC/TRUE 方案能否用于活体不需要先做动物实验直接比较两个时间系统从“开始测量”到“焦点能量到达目标位置”的总响应时间 (\tau_{sys})和该组织区域的散斑相关时间 (\tau_{corr})。(\tau_{sys} \tau_{corr}) 时焦点 PBR 会显著下跌数值上可以通过把相关时间设为变量拟合 PBR 衰减曲线来提取 (\tau_{sys})。读文献时注意区分“光场测量时间”和“系统响应时间”——锁相相机的 0.3 ms 只是测量端调控端没跟上整链路依然是断的。Zhang 等人在 2017 年用两帧曝光二值化结合铁电液晶 SLM 完成了第一套毫秒级完整 TRUE 系统测量和调控都压缩到了相关时间窗口内。工程上最快的提速路径是把 Lockin 相机的高并行测量与铁电液晶的高速刷新拼在一起同时把相位图从浮点灰度量化到二值化或者三值化用调制深度换时间。本文还有配套的精品资源点击获取 数字化波前整形DOPC：深层生物组织光学聚焦的关键技术
简介资源为一份关于基于光学相位共轭的数字化波前整形技术DOPC的完整技术文档主要面向生物光子学、光学成像和生物医学工程等领域的科研人员与高年级学生帮助理解如何利用波前整形克服生物组织散射、突破传统光学约1 mm的聚焦深度限制。文档先对比基于反馈、传输矩阵测量和光学相位共轭三条技术路线的优劣随后着重剖析DOPC系统中的相位共轭镜构成以及光场测量和光场调控两阶段工作机理并包含四步相移全息、相机与空间光调制器像素匹配、焦点PBR提升等实现细节可作为该方向的导览性综述或课题入门笔记。资源仅此1份docx文档压缩包大小约212KB便于快速下载阅读。目前已有123人学习下载。1. 光子无法穿透的 1 mm为什么需要数字化波前整形生物组织对光并不友好。折射率的不均匀分布会把一束准直光打散成散斑传统光学成像在组织里的工作深度被压在 1 mm 附近限制因素不是探测器灵敏度而是弹道光在出射前就已经消失殆尽。早年学界默认散射是随机破坏过程直到波前整形技术出现大家才意识到组织内的散射事件其实是确定性的光走了哪些路径、每条路径累积了多少相位延迟都可以被测量和反推。基于光学相位共轭的数字化波前整形DOPC是三条技术路线中效率最高的一种它用相位共轭镜对散射光场做时间反演一次测量就能拿到预补偿所需的最优相位图无需像反馈迭代法那样反复逼近因此在血流、心跳导致散射快速变化的活体场景中最有竞争力。本文按测量—调控—对准—引导星—速度这条链路拆开讲最后给出高速化路线的工程判断。2. 四步相移全息与相位共轭镜的构建2.1 两种相位共轭镜的取舍相位共轭镜负责把发散散射光的波前“翻过来”再送回散射介质。早期实现依赖非线性晶体的四波混频响应速度极快但产生的共轭波能量很低生物组织内聚焦需要足够光强来激发荧光或驱动光遗传晶体方案很难满足。DOPC 的替代思路是把相位共轭拆成两个数字化步骤相机记录散射光场空间光调制器SLM重新发一个共轭波前。SLM 上的每个像素可独立调制百万像素的模式数意味着对散射波前的刻画能做到非常细调制效率也远高于晶体方案。代价是系统复杂度上移——相机传感器面和 SLM 面必须严格光学共轭这一步做不好后续所有补偿都白费。2.2 光场测量四步相移恢复散射相位SLM 只能加载相位相机只能记录强度所以相位信息必须先编码到强度里才能被测量。做法是在散射光路里引入一束平面参考光与散射光干涉形成散斑然后对参考光依次附加 0、π/2、π、3π/2 四步相移记录四帧干涉图。第 k 帧强度分布为[ I_k A_R^2 A_S^2 2 A_R A_S \cos(\phi_S - \phi_R k\pi/2) \quad k0,1,2,3 ]式中 (A_R)、(\phi_R) 是参考光振幅和相位(A_S)、(\phi_S) 是散射光的振幅和相位。参考光相位在空间上近似均匀因此用两对差分帧就能把散射光相位从余弦项里解出来。import numpy as np def four_step_phase_shift(I0, I1, I2, I3): 四步相移法恢复散射光场相位分布。 I0..I3 : 参考光相位为 0, pi/2, pi, 3pi/2 时相机记录的干涉强度图 返回 : 相位图 phi, 弧度制, 范围 [-pi, pi] # 差分对消掉 A_R^2 A_S^2 直流背景 cos_term I0 - I2 # 2*A_R*A_S*cos(phi_S - phi_R) sin_term I3 - I1 # 2*A_R*A_S*sin(phi_S - phi_R) phi np.arctan2(sin_term, cos_term) # 取相位, 振幅项被自动约去 return phi这段代码的要点在arctan2它同时利用 sin 和 cos 两个正交分量相位输出范围是完整的 ([-\pi,\pi])不会出现单用 arccos 时正负象限混叠的问题。两帧差分把直流背景和 SLM 反射率不均匀性一并消掉这是四步相移比单帧离轴全息抗噪能力更强的原因。实际采集时四帧图必须在散斑相关时间内拍完否则相位信息已经被生物运动抹平这一点在高速系统里会成为主要瓶颈。得到散射相位分布 (\phi_S(r)) 后SLM 需要加载的是 (-\phi_S(r))也就是对每个像素做共轭。2.3 光场调控为什么共轭相位能让光“倒流”把 (-\phi_S(r)) 加载到 SLM 上让一束与参考光同源的平面波经过 SLM 反射等于给平面波逐点叠加一个与散射过程相反的相位。依据波动方程在时间反演下的对称性这个相位共轭波重新进入散射介质后会沿着散射路径原路返回在目标位置相干叠加成焦点。注意这里不是让光真的倒流而是用事先测得的相位延迟做预补偿使所有路径的光在目标位置重新同相。这个过程只需要一次测量加一次加载对比反馈迭代法动辄数百轮优化速度优势明显代价是测量阶段必须保证参考光的相位稳定性。2.4 三种波前整形路线的横向对比维度反馈迭代法传输矩阵法DOPC 相位共轭法最优相位图获取迭代优化需数百至数万次光强采样一次测量复传输矩阵矩阵规模随像素数平方增长单次测量散射光场一次共轭加载系统复杂度较低SLM 加反馈探测器即可较高参考光路需长期稳定中高相机与 SLM 需严格像素匹配对动态散射介质适应性差迭代期间散斑已变化差矩阵需频繁重测好单帧测量后可快速加载典型场景静态散射体、教学验证光纤模式控制、光通信生物活体聚焦、时变组织三条路线最终都在求解同一个问题预补偿相位图。区别只在于相位图从哪来、拿得多快。3. 像素级精度的系统校准从 PBR 61 到 52000 的对准工程3.1 PBR 的上限由 SLM 像素数决定焦点峰值与背景功率之比PBR是衡量聚焦质量的核心指标。若 SLM 有效像素数为 (N)理论上限是 (\pi N/4)一百万像素对应约 78 万 PBR。实际 DOPC 系统长期停留在数十到数百 PBR问题出在相机传感器面与 SLM 面之间存在 6 个自由度的对准误差三个平动 (\Delta x,\Delta y,\Delta z) 和三个转动 (\Delta\theta_x,\Delta\theta_y,\Delta\theta_z)。任意一个误差都会把单个 SLM 像素的共轭相位送到错误的位置散斑统计叠加后焦点被抹进背景里。3.2 物理掩模时代的粗对准2010 年第一套 DOPC 系统用物理掩模辅助像素匹配先用一相机标定掩模位置再对照 SLM 上加载的掩模图样逐像素找对应关系。这个方案能工作但光路复杂双相机切换引入额外像差对准精度完全依赖手动调节最终 PBR 约 600离理论值差三个数量级。这类方法的问题本质是“用实物标定去补偿光学误差”标定本身的误差无法自举消除。3.3 角谱传播补偿与 6 自由度搜索2014 年 Jang 的方案把问题从硬件搬到软件里先做一次粗略的机械对准再用角谱分析法在 SLM 上加载补偿相位把相机平面与 SLM 平面之间残余的传播误差数值化地补掉。import numpy as np from numpy.fft import fft2, ifft2, fftshift, ifftshift def angular_spectrum_propagate(field, dz, wavelength, pixel_size): 角谱法计算复振幅场自由传播 dz 距离后的结果。 field : 源平面复振幅 (对应 SLM 平面) dz : 传播距离, 正数表示沿光轴前进 wavelength : 工作波长 pixel_size : 像素物理尺寸, 与 wavelength 单位一致 ny, nx field.shape fx fftshift(np.fft.fftfreq(nx, dpixel_size)) fy fftshift(np.fft.fftfreq(ny, dpixel_size)) FX, FY np.meshgrid(fx, fy) # 衍射传递函数 exp(i * kz * dz), kz 由色散关系决定 kz_sq (1.0 / wavelength) ** 2 - FX**2 - FY**2 H np.exp(1j * 2 * np.pi * dz * np.sqrt(np.maximum(kz_sq, 0))) H np.where(kz_sq 0, H, 0) # 倏逝波分量直接置零 return ifft2(fft2(field) * H)这个函数解决的是 (\Delta z) 和部分离焦误差但 (\Delta x,\Delta y) 对应像素平移(\Delta\theta_x,\Delta\theta_y,\Delta\theta_z) 对应平面旋转需要在外层包装坐标变换。常见的做法是把 6 个参数依次扫描先用粗步长扫 (\Delta x,\Delta y) 找到焦点的大致位置然后扫 (\Delta z) 确定离焦量最后旋转三个角度做精细搜索每一步都以焦点强度作为反馈值保留最优参数。整个优化过程耗时约 10 分钟用贪婪策略从粗到细推进最终把 PBR 从粗对准的 61 提升到 52000。3.4 泽尼克多项式修正与 SLM 曲率补偿角谱方案假设系统像差是纯粹的刚体偏差实际 SLM 自身有表面曲率成像透镜也有高阶像差。2016 年 Azimipour 引入自适应光学里成熟的爬山算法在 SLM 上叠加矩形泽尼克多项式作为修正面形以焦点亮度为反馈逐阶扫描多项式系数修正转动误差和系统像差后 PBR 达到 121000。2018 年 Hemphill 更进一步用四步相移全息直接测量 SLM 表面曲率把曲率相位预先加载掉再用干涉法定标 SLM 上生成的十字符号完成空间映射在 5 mm 厚散射样品上拿到近一万 PBR达到理论极限的 14.6%。方案年份对准/补偿手段实现 PBRCui 物理掩模2010掩模 双相机标定~600Jang 角谱自动校准2014角谱传播补偿 6 参数贪婪搜索52000Azimipour 泽尼克修正2016爬山算法 泽尼克多项式系数扫描121000Hemphill 综合校正2018SLM 曲率预补偿 十字符干涉标定~10000提示角谱补偿代码里dz的单位必须和pixel_size一致否则相位因子差几个数量级。先验证空场传播一个已知距离的菲涅尔衍射结果再接入 DOPC 搜索流程。从 61 到 121000 的跨越说明 DOPC 系统到达理论极限之前校准精度是第一约束。即便到了 2018 年也还有 85% 的 PBR 空间没有榨干误差来自校准残余和 SLM 相位灰阶量化的固有损失。4. 超声引导星与 TRUE 聚焦外差干涉的时频分析4.1 组织内聚焦需要“坐标”前面讨论的 DOPC 解决的是穿透散射介质后的聚焦真正有生物医学价值的是在散射组织内部聚焦。问题变得棘手光在组织内部没有清晰的参考点必须引入一个“引导星”来标记目标位置。侵入式引导星包括荧光分子、二次谐波纳米粒子、磁性粒子和基因编码蛋白对比度高、尺寸小但都要向组织内注入异物安全性受限。内源性散射体如血流中的红细胞虽然无侵入但成片出现无法单独操控。理想的引导星要同时满足非侵入、尺寸可控、可自由扫描三个条件。4.2 AOM 频移与声光标记聚焦超声天然满足这三个条件。超声在散射组织内受散射影响远小于光可以通过换能器精确聚焦到目标位置。当散射光穿过超声焦点时声光效应会给一小部分光子叠加一个频移这部分光子被称为标记光其余散射光频率不变。关键在于分离标记光与背景散射光用声光调制器AOM工作频率 40~60 MHz把参考光频率搬到标记光附近让参考光与标记光外差干涉形成拍频而未被标记的散射光因为与参考光频率差太大拍频超出相机帧率表现为稳定背景。相机记录到的强度为[ I(t) A_R^2 A_T^2 A_S^2 2A_R A_T \cos(2\pi \Delta f t \phi_T - \phi_R) ]其中 (\Delta f) 是参考光与标记光的频差(\phi_T) 正是超声焦点处散射光携带的相位提取出它相位共轭就能把光送回声学焦点。这套方案被称为 TRUETime-reversed ultrasonically encoded optical focusing。4.3 外差拍频信号的锁相提取四步相移在低频静态场景里够用但拍频信号的相位随时间变化需要锁相检测来同时提取振幅和相位。锁相 CMOS 相机每个像素内置高速光电转换和移相电路并行输出信号的振幅和相位数值等价于下面这段正交解调import numpy as np def lockin_demodulate(signal, tvec, f_beat): 正交锁相解调: 从外差干涉时间序列中提取振幅与相位。 signal : 相机单像素记录的拍频强度序列 tvec : 对应的时间序列 (秒) f_beat : 参考光与标记光的频差 (Hz), 由 AOM 驱动频率决定 # 本地振荡器: 同频余弦和正弦, 系数 2 用来抵消解调增益 ref_cos 2.0 * np.cos(2 * np.pi * f_beat * tvec) ref_sin 2.0 * np.sin(2 * np.pi * f_beat * tvec) # 混频得到 I/Q 分量, 再用滑动平均做低通滤掉 2*f_beat 分量 win int(np.round(len(tvec) / (f_beat * (tvec[-1] - tvec[0])))) # 约一个拍频周期 kernel np.ones(win) / win I np.convolve(signal * ref_cos, kernel, modesame) Q np.convolve(signal * ref_sin, kernel, modesame) amp np.sqrt(I**2 Q**2) phase np.arctan2(Q, I) return amp, phase说明几点win取一个拍频周期做平均窗口低于这个窗口滤不干净高频项高于这个窗口会抹掉超声焦点处散射光的快速变化ref_cos和ref_sin必须严格正交频率偏差会导致相位随测量时间线性漂移实验上通常用与 AOM 同一时钟源驱动的信号发生器来保证相干性。锁相相机的 0.3 ms 测量速度就是这么来的——每个像素的模拟电路并行完成解调不需要把四帧强度图传回计算机再算。代价是像素规模受限早期锁相 CMOS 只有 300×300 像素。4.4 引导星性质决定焦点质量引导星的位置和尺寸直接决定焦点质量。超声的波长比光长几个数量级声焦点的体积远大于光学衍射极限因此 TRUE 单次聚焦的焦点偏大。迭代 TRUE 可以逐轮收缩焦点第一轮聚焦到声焦点第二轮以第一轮焦点为新的引导源继续优化实验上观察到 2~3 倍分辨率和约 20 倍 PBR 的提升理论上可以收敛到单个散斑大小。2013 年 Judkewitz 用散斑涨落编码做引导星通过数千次光场测量在组织内实现了约 5 μm 的单散斑聚焦代价是耗时数小时。2015 年 Ruan 用聚焦超声击破血管内微气泡来提高信号对比度焦点亮度和尺寸都有明显改善。引导星类型侵入性尺寸控制生物安全性可操控性荧光分子 / 纳米粒侵入亚微米需注射存在代谢负担只能随注入位置固定磁性粒子侵入微米级有安全争议可用磁场粗略引导内源散射体动态非侵入成片出现高几乎不可操控聚焦超声非侵入由声透镜焦深决定高可自由扫描5. 散斑相关时间约束下的高速 DOPC 路线5.1 动态散射介质的时间窗活体组织的散斑相关时间在数十微秒到数毫秒量级主要贡献来自流动红细胞。波前整形的完整流程必须在一个相关时间内完成测量、传输、计算、加载相位、聚焦。传统四步相移加 SLM 的流程需要几十毫秒散斑早就变了共轭相位打在已经失效的散射状态上焦点自然出不来。5.2 四条高速化路线对比方案测量方式调控器件系统响应主要代价DMD FPGAWang 2015单帧二值化相位DMD刷新 22 kHz5.3 ms调制效率降至 (1/2\pi)不能兼容超声引导星锁相相机Liu 2016并行锁相解调0.3 ms未完成调控测量级 0.3 ms像素仅 300×300传输速率低离轴全息Hemphill 2017单帧强度图提取光场普通 SLM约 9 ms牺牲像素分辨率计算繁琐铁电液晶 SLMLiu 2017两帧曝光二值化相位铁电液晶毫秒级响应约 6 ms完整 TRUE 聚焦相位调制能力有限器件成本高5.3 验证方法用系统响应时间做可行性判据判断一套 DOPC/TRUE 方案能否用于活体不需要先做动物实验直接比较两个时间系统从“开始测量”到“焦点能量到达目标位置”的总响应时间 (\tau_{sys})和该组织区域的散斑相关时间 (\tau_{corr})。(\tau_{sys} \tau_{corr}) 时焦点 PBR 会显著下跌数值上可以通过把相关时间设为变量拟合 PBR 衰减曲线来提取 (\tau_{sys})。读文献时注意区分“光场测量时间”和“系统响应时间”——锁相相机的 0.3 ms 只是测量端调控端没跟上整链路依然是断的。Zhang 等人在 2017 年用两帧曝光二值化结合铁电液晶 SLM 完成了第一套毫秒级完整 TRUE 系统测量和调控都压缩到了相关时间窗口内。工程上最快的提速路径是把 Lockin 相机的高并行测量与铁电液晶的高速刷新拼在一起同时把相位图从浮点灰度量化到二值化或者三值化用调制深度换时间。本文还有配套的精品资源点击获取