自适应光学在双光子贝塞尔光片荧光显微镜中的应用【附代码】

✨ 长期致力于双光子贝塞尔光片荧光显微镜、自适应光学、直接波前探测与校正、哈特曼探测器、液晶空间光调制器研究工作，擅长数据搜集与处理、建模仿真、程序编写、仿真设计。
✅ 专业定制毕设、代码
✅如需沟通交流，点击《获取方式》

（1）基于菲涅耳衍射的双光子贝塞尔光束像差模拟与校正算法：

为了消除生物组织不均匀折射率引起的像差对贝塞尔光片质量的影响，建立了一种基于菲涅耳衍射原理的像差校正算法。将贝塞尔光束的产生与波前像差校正整合到一个液晶空间光调制器上，通过加载复合相位图同时实现轴棱锥相位和像差校正相位。利用哈特曼探测器测量样品中等晕区内的双光子激发荧光焦点（通过空间光调制器在样品中聚焦产生），获得前向像差。像差用前50项Zernike多项式拟合，校正相位为拟合系数的负值。模拟计算表明，当像差大小超过0.3λ时，必须校正至少50项才能将斯特列尔比恢复到0.8以上。在Rhodamine 6G溶液中，校正后双光子贝塞尔光束的强度恢复到了无像差时的89%，厚度从2.1μm压缩回1.0μm。

（2）基于荧光导星和哈特曼探测的直接波前探测方法：

针对贝塞尔光片线状照明导致传统哈特曼无法有效探测的问题，提出了一种荧光导星产生方法。空间光调制器在轴棱锥相位的基础上叠加一个聚焦透镜相位，使光束在样品中汇聚成一个高斯聚焦点，激发产生双光子荧光作为点状导星。该荧光被同一物镜收集后进入哈特曼探测器。为了匹配等晕区大小（约100μm），贝塞尔光束长度通过组合轴棱锥相位和离焦相位动态调节，长度可在100μm至300μm之间调节，同时光束厚度变化控制在40%以内。实验测得，经此法探测并校正后，成像深度从150μm提升到400μm，分辨率在200μm深度处仍保持0.6μm。

（3）逐等晕区扫描校正的大视场成像策略：

针对生物样品像差随空间变化的问题，设计了一套先探测后校正的成像流程。将成像视野划分为多个等晕区（每个区大小100μm×100μm），每个区内分别探测像差并保存校正相位。成像时，扫描振镜驱动贝塞尔光片快速扫描，同时每切换一个等晕区，空间光调制器加载对应的校正相位。整个流程中，每个等晕区的像差探测耗时约0.5秒，校正相位加载约10毫秒。使用高浓度琼脂固定荧光微珠样品验证，校正后整个300μm×300μm视野内的横向分辨率一致性从之前的差异40%改善到差异小于12%，信噪比平均提高6dB。对小鼠肾脏组织成像显示，自适应光学校正后可以清晰分辨肾小球足细胞突起结构，而未校正时只能看到模糊轮廓。

import numpy as np from scipy.special import jv from scipy.fft import fft2, ifft2, fftshift class FresnelPropagator: def __init__(self, wavelength=0.8e-6, pixel_size=0.5e-6, N=512): self.lam = wavelength self.dx = pixel_size self.N = N self.k = 2*np.pi / self.lam def angular_spectrum(self, U0, z): fx = np.fft.fftfreq(self.N, self.dx) Fx, Fy = np.meshgrid(fx, fx) H = np.exp(1j * self.k * z * np.sqrt(1 - (self.lam*Fx)**2 - (self.lam*Fy)**2)) U_fft = fft2(U0) Uz = ifft2(U_fft * H) return Uz class ZernikeAberration: def __init__(self, max_order=50): self.coeffs = np.zeros(max_order) def phase_screen(self, rho, theta): # 简化的Zernike多项式求和 (仅示意) phase = np.zeros_like(rho) for i in range(len(self.coeffs)): n = int(np.ceil((-3 + np.sqrt(9+8*i))/2)) m = 2*i - n*(n+2) # 近似 if m == 0: Z = np.sqrt(n+1) * self.radial_poly(rho, n, 0) else: Z = np.sqrt(2*(n+1)) * self.radial_poly(rho, n, m) * (np.cos(m*theta) if m>0 else np.sin(-m*theta)) phase += self.coeffs[i] * Z return phase @staticmethod def radial_poly(rho, n, m): from sympy import binomial R = 0 for k in range((n-abs(m))//2 + 1): R += (-1)**k * binomial(n-k, k) * binomial(n-2*k, (n-abs(m))//2 - k) * rho**(n-2*k) return R class SLMCompositePhase: def __init__(self, axicon_angle=0.01, focal_length=0.5): self.alpha = axicon_angle self.f = focal_length def axicon_phase(self, x, y): r = np.sqrt(x**2 + y**2) return -self.k * self.alpha * r def lens_phase(self, x, y): return -self.k * (x**2 + y**2) / (2 * self.f) def correction_phase(self, zernike_phase): return -zernike_phase def composite(self, x, y, zernike_phase, mode='bessel'): base = self.axicon_phase(x, y) if mode=='bessel' else self.lens_phase(x, y) return base + self.correction_phase(zernike_phase) # 荧光导星聚焦模拟 def create_guide_star(slm_phase, propagator, z_focus): U0 = np.exp(1j * slm_phase) # 入射场 U_focus = propagator.angular_spectrum(U0, z_focus) intensity = np.abs(U_focus)**2 return intensity