超表面全息在仿真里和实验台上的差距估计只有被散斑噪声折磨过的人才会深有体会。仿真中重建图像锐利到能看清每一个线条实验台上同样的相位分布却总像蒙了一层细沙边缘还有明暗交替的鬼影。我最早把这些全记在加工误差头上直到一次偶然把单模激光换成经过毛玻璃扩散的照明图像里的颗粒噪声明显变淡才意识到问题有很大一部分来自激光空间相干性。超表面全息的技术路径其实不复杂用亚波长结构的相位调制代替传统空间光调制器让入射光在透过的瞬间被重排从而在观察面上叠加出目标图像。但激光作为照明光源时单模光束在空间上几乎是完全相干的这意味着超表面上哪怕只有几十纳米的尺寸波动都会在观察面上通过干涉被放大成强度起伏。这篇文章想复盘我从伪影成因分析到空间相干性调控落地的全过程聊聊怎样通过主动控制激光的空间相干性来抑制超表面全息重建中的伪影以及在真正动手时有哪些参数需要权衡。内容既包含定性和定量的原理分析也包含可以直接复现的建模流程和实验注意事项适合正在跟相干噪声较劲的研究生也适合想把超表面全息推向产品的光学工程师参考。1. 伪影问题不是超表面的错是激光太完美了超表面全息的基本结构是排布在界面上的亚波长散射单元每个单元可以理解成一个像素相位器。设计阶段我们通常用GS类算法计算目标图像对应的相位分布再用TiO₂纳米柱、非晶硅纳米柱等结构去近似。理论上只要每个单元的透射相位准确输出面上就会相干叠加出预期图像。但是实际样品总有加工误差、侧壁粗糙、几何形状偏差这些偏差相当于在理想波前上叠了一个随机相位和随机振幅。如果照明光是完全相干光理想信号波前和误差波前会以确定的相位关系干涉误差部分不会被平均掉而是会稳定叠加在目标图像上形成像被针刺过的噪声纹理。这就是散斑伪影的第一个来源。第二个来源来自超表面本身的周期性。不管目标相位如何纳米柱阵列在物理上就是一个周期为p的光栅结构。哪怕每个柱的相位设计得很理想这个周期也会在衍射面上产生周期性的衍射级次。当目标图像的频谱与高阶衍射级次发生交叠时图像边缘会出现复制像或者莫尔条纹。高相干光会让这些级次之间的干涉变得尖锐视觉上就像罩了一层网格状的纱。第三个来源是系统内部的寄生反射超表面基底背面、光学透镜表面、探测器窗口都可能反射部分光这些反射光和主光束来自同一个相干光源会在全息图上产生牛顿环状或条纹状的干涉图样。这类伪影尤其阴险因为它的周期和超表面本身无关纯属光路质量问题。上述三类伪影的共同点是它们都依赖稳定的相位关系。只要降低光源的空间相干性让不同位置的照明光不能保持固定的位相锁定这些干涉项就会被平均掉。但也别高兴太早——超表面全息的本质是干涉成像完全去掉空间相干性会让重建图像直接消失。所以真正的难点在于如何保留用来重建波前的那部分相干性同时抹掉放大误差的那部分相干性。1.1 相干照明下的散斑与周期性伪影散斑是相干光在随机散射体上发生多重散射后各散射分量在观察面上随机干涉的结果。超表面加工后必然存在纳米级的粗糙度这些粗糙度相当于无数个随机散射源。对于完全相干照明观察面上每一点的强度是所有散射源的振幅叠加相位差从0到2π随机分布于是强度呈指数型概率分布。常用的衡量指标是散斑对比度C σ_I / ⟨I⟩完全相干照明下C接近1也就是说噪声标准差和信号均值几乎一样大重建图像里散斑的能量是相当可观的。如果用LED或者宽带光源照明由于不同波长和不同横向位置的发光点不相干散斑图案会被加权平均对比度迅速降到0.1以下。超表面全息不能直接用LED因为需要近单色光来控制相位所以只能在保留部分相干的同时尽量增大独立散斑图样的数量。周期性伪影则不一样它由超表面固定周期导致在时间上并不随机。传统空间光调制器的像素周期通常在微米量级超表面的周期则缩短到了几百纳米这就把高阶衍射级次推到了更大的衍射角外理论上避免了大部分高阶混叠。但实验里目标图像频谱很宽总会有一部分能量落在高阶衍射级附近周期性鬼像就会冒出来。降低空间相干性并不会消除这些周期性级次但会让级次边缘的干涉条纹变柔和让视觉上的网格罩子不再扎眼。1.2 从相干长度到伪影的空间尺度映射空间相干性可以用横向相干长度r_c来描述。直观理解如果在波前上取两个点间距小于r_c它们的光场相位关系能保持稳定间距远大于r_c两点之间的相位就是随机的。对于完全相干光r_c→∞对于自然光r_c可能只有微米甚至更小。把空间相干性和散斑颗粒联系起来可以用一个粗估公式。观察面上的散斑颗粒横向尺寸大约为 s ≈ λL/D其中λ是波长L是超表面到探测器的距离D是照明光斑在超表面上的有效直径。比如λ633nm、L300mm、D10mm则s≈19μm这和常见CMOS相机的像元尺寸5~10μm很接近所以图像上的散斑是一粒一粒清晰可见的。当我们降低空间相干性时原本一个相干照明区域被分割成N个近似独立的相干子区域N ≈ (D/r_c)²。每个子区域自己产生一幅散斑图案由于子区域之间互不相干探测器上积累的是N幅随机图案的非相干叠加散斑对比度近似按1/√N下降。在同样的D和λ下如果r_c从10mm降到2mmN从1变到25C的理论值从1降到0.2如果r_c降到0.5mmN400C约0.05散斑基本平滑了。这个估计当然很粗糙但足以说明数量级之间的关系空间相干长度越小独立模式数越多伪影的强度就越低。状态空间相干长度 r_c独立模式数 N理论散斑对比度 C≈1/√N单模激光远大于D1≈1弱散射2mm250.20中等散射1mm1000.10强散射0.5mm4000.05不过这个去相干过程也会影响全息重建本身。超表面的相位编码依赖相邻单元的干涉如果r_c减小到与单元周期p相当那么相邻单元之间的相干性就不复存在高频信息会被彻底洗掉。所以合适的r_c必须同时满足两件事一是相对于目标光斑D要足够小才能提供足够的统计平均次数二是相对于超表面单元周期p要足够大才能保住波前重构所需的局部相干性。这两个约束之间就是我们可以做文章的区域。2. 空间相干性调控把理想光源拉回现实但要有度既然问题源于过于完美的空间相干性最直接的思路就是主动去打破它但绝不是把激光换成LED就能了事。超表面全息本质上是靠干涉编码和重建波前如果入射光完全没有空间相干性每个单元都只和自身耦合多单元干涉无法形成重建图像会彻底消失。因此调控的目标是找到一个中间地带既能保留宏观波前重构所需的相干耦合又能抹掉微观随机相位误差带来的干涉伪影。这个中间地带可以用r_c来标定。r_c太小等效于把超表面像素之间的对话切断r_c足够大等效于错误波前之间仍然产生稳定干涉伪影原形毕露。实际操作中r_c的选择范围通常在微米量级但具体数值高度依赖系统设计。2.1 部分相干光的常用工程实现在实验室里想把激光改造成可调的部分相干光有几个现成的方案各有各的脾气。旋转毛玻璃可能是最经典的办法。把激光先聚焦到旋转的毛玻璃上光透过随机粗糙表面后被随机相位扰动探测器在积分时间内看到的是很多幅不同散斑图案的平均等效于空间相干性下降。通过调整光斑在毛玻璃上的直径、毛玻璃的颗粒度、转速可以连续调节等效r_c。优点是结构简单成本低缺点也很明显如果需要极短的曝光时间旋转毛玻璃来不及做时间平均瞬时相干性仍然很高。空间光调制器SLM加载动态随机相位屏是更灵活的方案。SLM按帧切换随机相位分布配合相机采集可以在单帧内改变空间相干性。它的核心优点是可以用算法精细控制相位统计特性比如指定相干函数的高斯宽度或指数宽度。但市面上多数SLM的刷新率只有几十到几百赫兹如果曝光时间短于帧周期SLM本身的像素栅格反而会成为新的衍射源这一点特别容易踩坑。多模光纤本质上也是一个去相干器件。不同横模的光在光纤中走的光程不同在输出端叠加时横向相干度天然下降。把激光耦合进一段数米长的多模光纤输出端的空间相干长度可以小到微米量级。优点是全光纤化和紧凑化缺点是模式对光纤弯曲敏感输出强度分布不稳定。静态相位漫射体比如工程漫射体是另一个常见选项。它没有运动部件稳定性好而且可以通过设计随机相位谱来指定散射角分布。缺点是不可调更麻烦的是它会引入一个固定的散斑图样如果不继续做时间平均这个固定图样会和超表面本身的误差勾结起来形成新的伪影。所以静态漫射体更适合配合扫描振镜或移动平台使用让它相对光束运动起来。方案工作原理优点局限适合场景旋转毛玻璃随机相位随时间变化积分平均简单、低成本、连续可调短曝光失效有机械振动静态成像、长时间曝光SLM动态相位逐帧加载随机相位屏可编程、精确控制刷新率限制、效率低需要快速切换相干度多模光纤多横模非相干叠加紧凑、全光纤化模式不稳定集成化光源静态漫射体固定随机相位去相关无运动部件不可调、有固定散斑配合扫描或移动使用我们实验里最终选择的是旋转毛玻璃因为它能连续扫过r_c的整个取值范围方便先做完一整条伪影抑制—分辨率曲线。为了减少机械振动我们把毛玻璃安装在精密轴系上并且让聚焦光斑尽量小这样线速度高、振动影响相对小。2.2 调控尺度与图像分辨率的权衡关于r_c到底取多少我的习惯是先把目标图像想清楚。超表面的单元周期p决定了对波前采样的粒度。如果希望重建图像里有足够的细节那么至少要有几个连续的单元在相干体积内协同工作。经验上r_c取p的5~10倍是一个比较稳的起点。举个例子超表面p600nm那么r_c≈3~6μm。这些数值看起来很小但已经足够形成局部的波前重构。r_c和散射角大致满足 r_c ≈ λ/θ其中θ是照明光束经漫射体后的新增发散角。假设λ633nm想要r_c5μmθ大约是0.13 rad约7.5°。这个角度并不夸张但它会把激光光斑扩散到很大的面积。好在超表面本身尺寸有限我们可以在超表面平面前加一个中继透镜把扩散的光重新聚回有效口径同时基本保持r_c不变。这一步很关键否则功率密度下降得厉害。还有一点容易被轻视目标图像的频谱结构决定了r_c的敏感性。对于灰度均匀的大色块比如简单的字母F即使r_c只有1μm重建结果可能还能看对于高密度条纹或复杂图案r_c降到2μm就会明显发糊。所以不要指望一个固定r_c走天下最好在实验前用仿真先把目标图像对r_c的容忍度跑出来。后面第3节会展开具体建模流程。3. 完整设计链路从目标图到部分相干重建超表面设计和相干照明下的设计流程基本一致但部分相干照明会影响重建质量所以需要在每一个环节都留出相干度的考量。下面是我比较推荐的一条完整链路。3.1 超表面相位分布设计流程第一步是预处理目标图像。把待显示的灰度图转成振幅分布记为目标振幅A_t。如果后面打算用部分相干照明可以在这一步先做一个可选的频域预加重把目标图像的频谱乘以一个逆光学传递函数补偿后续相干性下降带来的高频衰减。但预加重幅度不能太大否则会放大噪声建议首次设计时不做等仿真和实验结果出来后再回来调整。第二步是用加权GS算法恢复相位。初始在观察面随机生成相位然后做傅里叶变换到超表面平面在超表面平面固定振幅为1或按实际孔径切趾保留相位再傅里叶逆变换回观察面把振幅替换成A_t保留相位如此循环。一般迭代200~500次就能收敛。加权GS可以改进光能均匀度比如在亮区缩小振幅误差权重在暗区扩大从而抑制零级亮斑。下面是一个极简的GS循环示意# target_amp: 目标振幅, N: 迭代次数 phase np.random.uniform(-np.pi, np.pi, (Nx, N_y)) for _ in range(N): slm_field np.exp(1j * phase) obs_field np.fft.fft2(slm_field) obs_field target_amp * np.exp(1j * np.angle(obs_field)) slm_field np.fft.ifft2(obs_field) amplitude np.ones((Nx,N_y)) # 均匀振幅约束 phase np.angle(amplitude * np.exp(1j * np.angle(slm_field)))第三步是把连续相位图按超表面单元周期p离散化。对于每个网格挑选一个纳米柱结构使其在波长λ下的透射相位尽可能接近目标相位同时透射率尽量高。这一步一般需要建立单元库用RCWA或FDTD扫描纳米柱的半径、高度、晶格周期得到相位-几何尺寸-透射率数据库。第四步是版图生成。把选定的纳米柱排布坐标输出成GDS文件交给工艺线。这里的一个细节是单元库扫描时要用与最终照明条件相近的入射角如果照明光经中继透镜后带有一定会聚角模拟中的边界条件要相应设置成平面波斜入射否则相位偏大实验会差很多。3.2 部分相干光源的数值建模方法仿真里不能简单地把目标图像做一次高斯模糊来模拟部分相干效应因为部分相干光的本质是不同横向模式之间的非相干叠加而不是单纯的频率滤波。一个实用的建模思路是随机相位屏系综方法生成M个独立的随机相位屏每一幅屏都作为一个人射场样本独立传播到观察面最后把M幅强度图像做平均。具体流程如下生成M个二维随机相位分布φ_m(x,y)统计特性与毛玻璃一致。如果模拟旋转毛玻璃φ_m一般取高斯随机场并带一个横向相关长度对应毛玻璃颗粒尺寸。每个样本的入射场 U_m A_in(x,y) * exp(iφ_m(x,y))其中A_in是入射高斯振幅或者平顶振幅。用角谱法或瑞利-索末菲衍射把U_m传播到超表面平面。乘以超表面的复透射系数T(x,y)由第3.1节的单元库给出包含相位和透射率。再传播到重建面得到强度I_m |ψ_m|²。最终重建强度 I (1/M)∑I_m。Python伪代码大致长这样M 200 I_total np.zeros((Nx,Ny)) for m in range(M): phi gaussian_random_field(Nx, Ny, corr_len) # 毛玻璃相位 U_in input_beam * np.exp(1j * phi) U_ms angular_spectrum_propagate(U_in, z1) # 到超表面 U_out U_ms * meta_transfer # 超表面复透过率 U_img angular_spectrum_propagate(U_out, z2) # 到重建面 I_total np.abs(U_img)**2 I_recon I_total / M有几个容易犯的错一是M取得太小比如只有10~20幅屏平均后的强度还会带有明显的样本噪声建议至少200幅。二是随机相位屏的相关长度要和实际毛玻璃对应相关长度应该接近光斑直径的1/10到1/100否则等效r_c不准。三是角谱法的采样窗口要足够大避免衍射场在边界折叠。我的经验是初始网格设为光斑尺寸的4倍以上边缘用零填充。3.3 重建算法与伪影量化指标仿真跑完以后需要量化伪影抑制效果。我一般同时报三组指标散斑对比度C、PSNR和SSIM。散斑对比度C的计算要特别注意区域选择。如果把整个重建图都纳入统计零级亮斑或者信号本身的高强度区域会拉高均值导致C被低估。我通常把重建区域分成信号区和背景区分别计算C。背景区的C更接近真实散斑水平信号区的C则能反映图像本身的信噪比。PSNR是对比重建强度和目标图之间的均方误差。使用前需要把重建图做灰度归一化并和目标图对齐。SSIM则能反映结构相似性对边缘模糊更敏感。两者的组合可以区分噪声被抹平和图像被抹糊两种完全不同的效果。还有一个实用的伪影抑制率定义ASR (C_相干 − C_部分) / C_相干 × 100%其中C_相干是用完全相干照明重建时的散斑对比度C_部分是用部分相干照明时的对比度。ASR能直观地告诉你去相干手段把噪声压下去了多少。需要注意ASR只反映噪声量不能反映细节损失所以一定要和PSNR、SSIM一起看。下面是一组典型仿真结果示例展示了r_c扫描时的性能变化条件相干长度 r_c (μm)散斑对比度 CPSNR (dB)SSIM主观效果完全相干∞0.8718.60.60散斑明显边缘鬼影重弱调控100.4023.40.73噪声减少细节仍可辨中等调控50.1827.10.82伪影抑制好图像清晰较强调控20.0824.80.66细节变虚边缘钝化过度调控0.50.0517.90.30图像基本失效从这组数据能明显看到甜点出现在r_c5μm附近散斑对比度降到了0.18PSNR和SSIM都是最高。继续把r_c降更小散斑确实更干净但全息重建本身的细节损失反而让PSNR跳水。所以所谓调控真的不是越弱越好而是要在重建能力和噪声抑制之间找平衡点。4. 实测中的关键参数与踩坑记录仿真跑得再漂亮上了实验台该踩的坑一个都躲不掉。下面这几个问题都是我实际摸过一遍之后总结出来的建议动手前先看。4.1 空间相干长度怎么测才算准先说测量方法。实验室里最经典的就是杨氏双缝在重建光路中放一簇不同缝间距的双缝板观察每个缝距下的条纹可见度V (I_max−I_min)/(I_maxI_min)。在Schell模型下V(d)大致按exp(−d²/2r_c²)衰减拟合出r_c。这个方法原理简单但有两个隐性坑一是必须保证光斑均匀照亮两条缝如果光束中心不对称两缝的光强不等可见度会被系统性拉低测出的r_c偏小二是缝间距范围要覆盖从几十微米到毫米的量级只测几个间距很容易拟合出一根歪曲线。如果不想花时间做双缝也可以用剪切干涉仪或商用波前探测器直接测量复相干度。我们实验室没有专用设备就用了一套横向剪切板给两个剪切像做条纹分析速度比双缝快精度也够用。重要的是测量位置尽量放在超表面实际放置的平面上因为光路中的透镜会改变r_c在光源端测到的值不能直接代表超表面入射端的值。还要额外注意时间尺度。旋转毛玻璃产生的扰动有时间特征尺度大致等于毛玻璃上光斑线速度除以散射颗粒尺度。如果相机的曝光时间远大于这个特征时间看到的是良好的时间平均如果曝光时间接近甚至小于特征时间看到的就是瞬时散斑伪影依旧。我有个简单的经验调整光斑直径和转速让扰动时间小于曝光时间的1/10。可以做一个快速验证把相机曝光时间减小一半如果散斑对比度明显变高说明时间平均不充分需要提转速或增大线速度。4.2 降低相干性后图像对比度下降怎么办这是最多人问的问题。把毛玻璃装进去以后重建图像确实干净了但整体亮度也肉眼可见地掉下来。原因很直接毛玻璃把能量从一个集中小斑扩散成了大角分布的散射光进入超表面口径的只有一小部分再加上随机相位在不同模式下分栏每个模式的有效能量都被稀释了。要应对可以从几个方向同时使劲。第一选散射角尽量小的毛玻璃并用透镜组把散射光重新耦合进超表面的口径。我们用的毛玻璃散射角在6°~10°之间经过中继透镜后功率损失可以控制在3dB以内。如果散射角太大功率损失会非常快而且r_c会偏小重建质量反而崩。第二适当提高激光输出功率但要盯着安全问题尤其是在会聚光路上。第三通过超表面设计补偿。在GS算法中加入部分相干照明下的光学传递函数补偿可以让中高频信息更强。实际操作时我先用第3.2节的仿真算出当前r_c下的系统OTF再把目标图像频谱除以这个OTF限制最大增益再跑GS。实测下来中等调控下PSNR能再提高2~3dB边缘也更锐利。但增益限制设得太高会引入振铃需要盯着看。第四不要把注意力全放在降噪上。有时牺牲一点散斑抑制换取亮度和对比度反而整体观感更好。我在实验里做了r_c5μm和r_c8μm两组后者的散斑对比度只从0.18涨到0.26但图像亮度和边缘锐利度明显更好最后选了8μm作为系统工作点。所以说仿真给的是参考曲线最终取舍还是得用眼睛做主观评价。4.3 一组对照实验数据我们用的样品是400×400像素的超表面周期p600nm工作波长633nm重建距离300mm照明光斑经透镜压缩到直径约8mm目标图像是一个带高密度线条纹理的字符图案。旋转毛玻璃放在光源之后通过调整光斑聚焦大小和转速来改变r_c。下面是一组实际采集数据的记录条件转速 (rpm)估算 r_c (μm)背景散斑对比度PSNR (dB)SSIM静态漫射体0—0.7219.80.58低速旋转50120.4622.60.70中速旋转30060.2126.90.81高速旋转100030.1225.20.71数据里有个细节很有意思静态漫射体的背景散斑对比度是0.72远低于完全相干下的理论值1但PSNR仍然不高。因为静态漫射体确实把空间相干性降下来了但它产生的散斑图案是固定的没有时间平均无法彻底抹平。只有让毛玻璃转起来让散斑图案不断变化才能在曝光时间内实现真正的非相干叠加。这就是为什么不能拿一张静态漫散射片去冒充部分相干调控的原因。另一个细节是转速从300rpm提到1000rpm后背景对比度从0.21降到0.12但PSNR从26.9掉到25.2。原因是高转速带来一定机械振动图像在探测器上有微米级的位移等效成额外的模糊。所以这台装置里300rpm反而是甜点而不是转得越快越好。5. 这项策略能用在哪些场景以及下一步思路激光空间相干性调控这条思路并不仅仅是为了解决我手头这一个样品的问题。超表面全息要走向实际应用散斑伪影几乎都是绕不开的坎而空间相干性调控恰好是普适性很强的手段。第一个典型场景是裸眼三维显示。超表面全息显示如果直接用激光照明观察者会看到明显的颗粒噪声动态画面的散斑还会随时间闪烁非常影响体验。在维持全息图像轮廓的同时降低空间相干性可以让显示画面接近自然光照下的观感。这里有个限制是刷新率旋转毛玻璃对高速显示不友好比较合适的做法是用SLM在每一帧全息图案之间加载不同的随机相位屏再配合人眼或者相机的时间积分等效降低相干度。第二个场景是激光加工中的多焦点整形。超表面可以把一束激光分成几十个焦点用于并行微加工。但相干照明下焦点之间会出现复杂的干涉条纹加工面能量分布不均匀。把照明光改成部分相干后每个焦点的旁瓣和干涉背景都能被压低能量更集中加工一致性更好。当前的部分相干光配合超表面做多焦点已有一些文献报道但多数是在空间光调制器系统上验证的换成超表面之后需要注意衍射效率和热效应。第三个场景是显微照明。超表面可以做平顶光束整形或者结构光照明但激光散斑一直是显微成像的痛点。空间相干性调控可以与现有的减相干散斑抑制方案结合起来使用例如在照明端加入动态漫射体在探测端配合快速采集获得更干净的照明场。再往下一步想我在琢磨能不能把去相干功能直接做进超表面本身。传统方案需要外置毛玻璃或SLM那就有体积、功耗、可靠性问题。超表面既然能任意调控相位理论上可以把目标相位和一个随机相位叠加在一起输出波前本身的空间相干性就被设计值降低省掉外部运动部件。但这会牺牲部分衍射效率而且随机相位图样一旦固定时间平均带来的平滑就没有了可能更适合那些对单帧散斑不敏感的应用。另一个思路是结合深度学习用神经网络从部分相干重建图像中恢复高频细节这样可以把r_c推得更低让散斑抑制到几乎为零再由算法把图像拉回来。我们初步试过一个小网络在较强去相干条件下能把SSIM从0.66拉回到0.8以上效果可观但需要针对不同目标图像重训工程上还要再打磨。最后分享一个实际过程中的小经验判断空间相干性调控效果不能只看全局PSNR。我第一次统计整个重建区域的PSNR总是不高因为背景噪声把均值拉得很低后来把背景和信号分开统计才看清中速条件下的甜点其实比仿真预测的要宽。所以实验上先把指标定义清楚比调参数本身更重要。搞超表面全息的人往往盯着纳米结构不放但光源这一端的变量同样值得花时间去调。如果这篇经验能帮你少走一点弯路那就挺好。