很多人学傅里叶光学时最痛苦的点在于公式能背下来积分也能算但完全不知道它在描述什么物理图景。我当年也是这样——夫琅禾费衍射积分写了无数遍知道矩孔的衍射花样是sinc函数但这只是“会做题”。直到某天我把衍射积分和傅里叶变换在纸上对着看了一遍才发现光学里最值钱的东西根本不是那套新物理而是换了一种语言重新描述光场。这篇《傅里叶光学1》就是把第一课里最绕的几个概念用尽量朴素的方式拆开空间频率到底是什么夫琅禾费衍射为什么就是傅里叶变换透镜凭什么能做傅里叶变换以及4F系统和空间滤波这个经典应用怎么从频域去理解。适正在学光学的学生、刚入行的光学工程师以及做计算成像和机器视觉、想搞清楚镜头和成像本质的人。1. 光场为什么要“拆开”看平面波是光学的积木块1.1 平面波为什么是基元函数信号处理里有个很成熟的做法任何一段信号都可以看成不同频率正弦波的线性叠加。你分析一个音频不会去逐点看波形而是看它的频谱因为频谱把“这个信号里有什么成分”一次性说清了。光学里完全一样。任意一个复杂的单色光场都可以写成无穷多个平面波的叠加。平面波是最简单的光场振幅处处均匀等相位面是一个平面传播方向固定。它就像一个“积木块”所有复杂光场都能用这种积木拼出来。为什么必须是平面波因为波动方程是线性的叠加原理成立。一个平面波是波动方程最基本的解而线性系统的任何一个解都可以分解成基本解的叠加。这个逻辑和傅里叶级数展开一模一样只是把基函数从正弦波换成了平面波。你看一个复杂光场觉得它很乱说到底只是因为你在用一个固定的视角看它换到平面波视角一切就清楚多了。1.2 拆完之后的传播问题角谱法平面波当积木块的最大好处是传播过程极其简单。一个平面波沿某个方向传播从z0平面走到zd平面只需要乘一个相位因子exp(jk_z·d)振幅不变方向不变。这就是整个角谱理论的核心。所谓角谱就是“按传播方向展开的谱”。一个光场在z0处是什么样我先把它分解成一堆平面波每个平面波各自传一段距离分别乘上各自的传播相位最后再把它们叠加起来就得到任意距离处的光场。这个过程可以类比信号处理里的线性时不变系统平面波是特征函数自由空间传播是线性系统系统对每个特征函数的作用就是乘一个复增益。信号处理里你会把信号分解成不同频率的复指数看系统对每个频率的响应光学里你把光场分解成不同方向的平面波看每个方向在空间里怎么传。两者本质上是同一个套路只是从时间换成了空间。1.3 衍射的本质是频谱演化不是“边缘效应”教材里讲衍射经常画那种边缘绕过的示意图给人感觉衍射是“光碰到边之后拐弯了”。这个图像不能说错但很容易让人忽略更本质的东西。一个无限大的理想平面波在自由空间里传多久都不会衍射因为它的角谱只有一个方向分量没有任何展宽。但一旦遇到孔径光场被空间截断角谱就展宽了——出现了大量偏离原传播方向的高空间频率分量。这些分量传播速度在z方向上的分速度不一样互相干涉光场花样就会随距离变化。这就是衍射。所以衍射的实质是空间截断导致频谱展宽随后各频谱分量在传播中相位演化不一致叠加结果随距离改变。孔径越小频谱展宽越厉害衍射越明显。这正是傅里叶变换尺度定理的光学版本空域越窄频域越宽。想清楚这一点后面所有衍射现象都会变得非常顺理成章。2. 空间频率才是傅里叶光学的第一语言2.1 空间频率和时间频率是亲兄弟时间信号x(t)A·cos(2πf·t)f是每秒振动多少次单位Hz。空间信号里一束平面波斜着照到一个平面上它在x方向的投影是一个周期性分布E(x)A·cos(2πu·x)u就是每毫米有多少个明暗周期单位常用“线对每毫米”lp/mm。时间频率描述的是信号随时间变化的快慢空间频率描述的是光场在空间上变化的快慢。一个光场里有大片均匀亮区也有密集细纹理低空间频率对应那些变化缓慢的区域高空间频率对应细密的结构和边缘。这个类比特别有用做信号处理的工程师对“频率”有本能直觉但一到光学就开始晕其实就是没把这两个频率对应上。时间上频率越高信号变化越快空间上频率越高图案越密。脑子里有了这个对应关系傅里叶光学的门槛直接降一半。2.2 二维空间频率一个方向对应一个频率光场是二维分布所以空间频率也是二维的写成(u,v)。u对应x方向的周期变化快慢v对应y方向的。一束平面波沿某个方向传播在x方向和y方向上都有相位梯度所以它对应一个二维频率点(u,v)。这里最关键的一个关系是空间频率和传播方向是直接对应的。波长为λ的平面波如果传播方向与z轴夹角为θ那么它在垂直于z轴的平面上的空间频率就是usinθ/λ。方向越斜频率越高。换句话说频谱面上的一个点就代表了一束特定方向的平面波。这个视角会颠覆你理解“成像”的方式物体反射的光进入镜头后并不是简单地被“映射”到像面上而是先被分解成各个方向的平面波再由光学系统决定哪些方向能通过、哪些方向被拦掉。每一个方向的平面波就是一个空间频率分量。2.3 系统通频带为什么看不清就是丢高频所有光学系统的口径都是有限的所以能收集的平面波方向范围也有限。数值孔径NAn·sinθ定义了系统能收集的最大角度。对应的系统能通过的最大空间频率是[ v_c \frac{NA}{λ} ]这说明了什么说明分辨率极限的本质不是“像素不够”而是系统从物理上丢掉了高频空间频率。显微镜看不清楚细节不是因为探测器分辨率差而是物镜只收集了一个有限锥角内的光那些对应更细结构的更大角度衍射光根本没进到系统里。给做相机的读者一个对应F数≈F的镜头出射光锥角近似对应NA≈1/(2F)。550nm绿光下F/8的镜头截止频率大约227lp/mm。一个像元尺寸3.45μm的传感器奈奎斯特频率约145lp/mm。也就是说镜头本身的能力可能还没被传感器完全用满。这种匹配计算在日常光学系统设计里非常常用。3. 夫琅禾费衍射就是一次傅里叶变换3.1 从菲涅耳积分到傅里叶变换到底近似掉了什么标量衍射理论里菲涅耳衍射积分长这样[ U(x,y) \frac{e^{jkz}}{jλz} \iint U_0(ξ,η) \exp\left[j\frac{k}{2z}\left((x-ξ)^2(y-η)^2\right)\right] dξ dη ]这个公式猛一看很吓人但你只需要盯住平方项展开后的样子[ (x-ξ)^2 x^2 - 2xξ ξ^2 ]三项里x²在远场条件下近似对观察面只贡献一个整体相位因子不影响强度ξ²是物面自身的二次相位因子当距离足够远时可以忽略剩下的-2xξ刚好构成傅里叶变换的核exp(-j2πu·ξ)。夫琅禾费条件的判据是[ z \gg \frac{a^2}{λ} ]其中a是孔径的特征尺寸。满足这个条件物面上每个点的自相位近似一致衍射积分就退化成透射函数的一个傅里叶变换只不过自变量是空间频率坐标ux/(λz)。所以夫琅禾费衍射从来不是什么“额外近似”它就是傅里叶变换的物理实现方式之一。反过来想傅里叶变换是一个数学操作而光的远场传播天然就是做这个操作这就是傅里叶光学最底层的立足点。3.2 双缝、光栅、圆孔的频谱怎么一眼看穿一旦接受夫琅禾费衍射就是傅里叶变换很多衍射花样都可以直接“看”出来双缝的透射函数是两个δ函数的卷积乘以单缝宽度它的傅里叶变换是一个宽的sinc包络乘上细的余弦条纹。所以你在屏上看到的双缝花样其实是“单缝衍射调制的双缝干涉”——这句话在频域里根本不需要背卷积定理一说就通。光栅是梳状函数梳状函数的傅里叶变换还是梳状函数所以光栅衍射是离散的级次级次间距对应空间频率间隔1/d。光栅方程sinθnλ/d本质就是空间频率分量un/d在自由空间传播时对应的方向角。矩孔的衍射花样是二维sinc函数圆孔因为圆对称性傅里叶变换是贝塞尔函数第一暗环半径1.22λz/D就是艾里斑边界。这些经典结论全部可以统一到“衍射花样透射函数的傅里叶变换”这一句话里。以后看到任何衍射花样第一反应应该是这个是哪个函数的傅里叶变换而不是去翻积分公式。3.3 实验里怎么“看到”傅里叶变换用透镜拉近远场夫琅禾费条件z≫a²/λ在实验室里经常是不现实的。1mm的孔径532nm绿光下条件要求z远大于1.9米。在桌面上做实验这个距离很难摆得开。透镜解决了这个问题。一束平行光经过透镜在后焦面上会聚如果把孔径样品放在透镜前方后焦面上的复振幅分布正好是孔径透射函数的傅里叶变换。无穷远处的夫琅禾费衍射花样被透镜“拉”到了它的后焦面上。这就是为什么几乎所有光学信息处理实验都用透镜它把抽象的频谱面变成了物理上可触及的一个平面。实际操作时你可以直接用激光扩束后照一个光栅透镜后焦面放白纸或CMOS几秒钟就能看到清晰的衍射级次。那个面上的花样就是光栅透射函数的二维傅里叶变换的幅度分布。这比在远处架屏幕靠谱得多也清晰得多。4. 透镜凭什么是傅里叶变换器后焦面的秘密4.1 薄透镜的相位变换函数要理解透镜为什么能做傅里叶变换得先知道透镜对光场做了什么。理想薄透镜可以建模为一个相位变换[ t(x,y) \exp\left(-j\frac{k}{2f}(x^2y^2)\right) ]这个式子看着抽象但物理图像很简单凸透镜中间厚、边缘薄对离轴位置的光线相位延迟更小等效于给光场乘一个会聚的二次相位因子。平面波经过透镜后变成了向焦点会聚的球面波正是这个二次相位在起作用。平时几何光学说“透镜把平行光会聚到焦点”你已经听得太多了但在傅里叶光学里你需要记住的是它的另一面透镜对入射光场施加了一个二次相位变换。这个相位变换和自由空间传播中的二次相位因子刚好符号相反所以传播过程中积累的相位弯曲被透镜抵消最终在傅里叶面上得到一个干净的变换结果。4.2 前焦面放物后焦面就是严格频谱面物体放在透镜前焦面时物光先传播一个焦距距离到透镜被透镜乘上二次相位再传播一个焦距距离到后焦面。推导出来的结果是后焦面上的复振幅分布就是物透射函数的严格傅里叶变换连二次相位因子都没有残留是真正的“纯”频谱。物体不放前焦面时后焦面仍然近似是频谱面但频谱会多一个与物位置有关的二次相位因子。如果你只是探测光强这个相位不影响结果但如果你要做相干叠加、滤波后再成像这个相位就会改变后续的衍射行为。所以严格的光学信息处理实验都要求物放在前焦面。还有一个经常被忽略的细节频谱面的坐标换算关系是[ u \frac{x}{λf} ]波长不同同一个空间频率对应的位置就不同。宽光谱光源照明的时候频谱面会看到明显的彩色错位零频点甚至会呈现彩虹色。这也是很多频谱面实验要用单色光或激光的原因。4.3 4F系统光学信息处理的“标准间”4F系统是傅里叶光学里最经典的架构物面到第一个透镜距离为f第一个透镜后焦面就是频谱面频谱面到第二个透镜也是f第二个透镜的后焦面就是输出像面。整个过程“变换-滤波-逆变换”输入图像经过一次光学傅里叶变换在频谱面上被人为修改再做一次逆变换得到处理后的像。在这个架构下滤波变成了一个极其直观的操作在频谱面插一根针挡掉零频像面的背景就被去掉插一个针孔只让低频通过图像立刻变平滑挡一束高频方向的频谱对应的纹理方向就会消失。频谱面上的“点”就是“方向”也就是“频率”这就是4F系统最有价值的地方。这个架构不只是教学实验里玩玩。空间相干度测量、光学相关识别、图像边缘增强、暗场/相衬显微的硬件实现很多都能在4F系统的框架下理解。哪怕你后面做的是数字图像处理这个“频域滤波”的图像也值得刻在脑子里。5. 阿贝成像与空间滤波从频域重新理解成像5.1 像不是直接生成的是频谱“干涉”出来的阿贝1873年研究显微镜分辨率时提出了一个在今天看来非常“傅里叶光学”的观点成像不是光线简单地从物点映射到像点而是经历了两次衍射过程。第一次衍射光通过物结构后形成各级频谱第二次衍射这些频谱分量作为有效波源继续传播在像面干涉叠加重建出物体的像。换句话说成像的本质是频谱的合成。像的质量取决于系统到底让多少频谱分量参与了最终干涉。高频频谱被系统切得越多像的边缘就越模糊低频频谱缺失像的整体轮廓和背景亮度就会异常。这个观点对现代光学特别重要。现在做计算成像的人天天说的“编码孔径”“相位恢复”“叠层成像”本质上都是在跟频谱打交道要么想方设法恢复被系统丢掉的高频要么故意对频谱做调制来获取更多信息。5.2 挡频谱等于改图像阿贝-波特实验阿贝-波特实验是理解空间滤波的最佳方式实验布置就是一个标准的4F系统物面用光栅频谱面放各种滤波器件。做完这个实验几乎所有频域结论都能变成直观记忆只让零频通过像面变成均匀亮场。因为零频表示全场平均值它不能携带任何空间变化信息。挡掉零频像的背景变暗边缘变成亮线。因为高频对应边缘和细节滤掉零频相当于做了一次高通/微分。只让水平方向的频谱通过原本的二维图像输出只保留竖直方向的结构。频谱方向对应图像纹理方向这一点体现得淋漓尽致。这个实验为什么值得亲手做一次因为做完以后你看任何光学成像问题脑子里都会自动浮现频谱面的样子一个图像的内容落在低频中心边缘和纹理分布在高频纹理的方向对应频谱的方向。这个直觉比背一百条成像公式都有用。5.3 从频域看暗场、相衬和全息暗场显微的原理就是频谱面滤波普通显微镜下透明样品几乎看不见因为它不吸收光只改变相位暗场方法用挡光板挡住直接透射的零级光只让被样品散射的高频光进入成像系统结果背景变暗而样品细节变亮。本质上就是“挡零频”这个操作的显微镜版本。相衬显微镜更进一步。透明样品不改变振幅只改变相位而探测器只能测强度泽尼克的方法是在频谱面的零频位置放一个移相片让零级光相位移动约π/2。这样零级光和衍射光在像面干涉时相位差就转化成了强度差。这也是频谱面滤波直接改善成像能力的经典案例。全息术就更不用说了记录的时候用参考光和物光干涉本质上就是在频域或空域对复振幅做编码用强度记录下来重现时靠衍射把相位信息从编码中解出来。理解了“傅里叶光学频域视角”再去看全息、相干衍射成像甚至最新的傅里叶叠层显微FPM都不会觉得它们是独立的新技术只是同一套思想的延伸。6. 实验笔记找频谱面、估极限分辨、绕开常见翻车点6.1 桌面实验三步走想亲眼验证“频谱面就是傅里叶变换面”不需要复杂设备。一套最简配置一台激光笔或小型激光器一个扩束准直透镜组一片透射光栅或针孔样品一个焦距50~200mm的凸透镜一块白卡纸或任意CMOS相机。第一步把激光扩束并准直成直径1~2cm的平行光。准直质量非常重要入射光如果不是平面波波前带曲率就会给物光额外引入一个二次相位频谱面位置和花样都会受影响。第二步样品放在透镜前焦面附近。样品可以是一维光栅、单缝、双缝或者打印在透明片上的任意图案。光线经过样品后紧接着进入透镜。第三步在透镜后焦面放屏幕或相机。怎么确认那是后焦面很简单前后移动屏幕找到光斑最锐利、衍射级次最清晰的那个位置。这个位置就是频谱面。别嫌这点时间省不掉我见过太多人把屏幕放在“大约一个焦距”的位置就开拍结果花样边缘发虚还找不到原因其实只是焦面偏了两厘米。6.2 频率估算公式与一组常用换算傅里叶光学里最常用的三个估算公式值得贴在实验本上系统截止频率(v_c NA/λ)数值孔径越大、波长越短能通过的空间频率越高。传感器奈奎斯特频率像元尺寸p奈奎斯特频率 (f_N 1/(2p))。频谱面坐标换算(x u·λ·f)反过来 (u x/(λf))。举个例子532nm绿光100mm焦距透镜你在频谱面上把针孔移到离光轴5mm处对应的空间频率是多少算一下(u 0.005/(532×10^{-9}×0.1) ≈ 94lp/mm)。这就是你亲手滤掉的频率。如果你换200mm焦距透镜同样是5mm位移对应只有47lp/mm频谱被摊得更开更容易精确操作。这就是长焦透镜在4F系统里受青睐的原因。下面这组常见数值可以临时查系统/参数数值截止或奈奎斯特频率显微镜物镜 NA0.65绿光550nm截止频率约1182 lp/mmF/8镜头550nm截止频率约227 lp/mm像元3.45μm的传感器奈奎斯特频率约145 lp/mm像元1.4μm的传感器奈奎斯特频率约357 lp/mm把镜头截止频率和传感器奈奎斯特频率放在一张表里对比很多“为什么这个镜头配这个相机合适/不合适”的问题会立刻有答案。6.3 四个真实踩坑记录第一个坑照明光没准直。我早期做4F系统实验时图省事直接用点光源照明结果频谱面位置会随样品轴向移动而漂移换上准直光之后所有现象都变得干净利落。原因是波前曲率会给物光加一个二次相位等于额外引入了一个等效透镜。第二个坑透镜口径不够。透镜本身口径有限相当于一个低通滤波器高频分量还没到频谱面就被切掉了。这时候你在频谱面上看到的频谱和理论频谱不一样滤波效果也混杂了系统本身的响应。做实验最好选口径明显大于样品尺寸的透镜。第三个坑频谱面光强反差极大。零级亮斑常常刺眼而高频细节暗得看不见。这是真实频谱的样子不是设备问题。处理办法是用中性密度衰减片把零级压暗或者分多次曝光。强行一次拍全不是过曝就是欠曝。第四个坑把频谱面强度当成完整频谱。探测器只能测强度测不到相位。只凭频谱面的强度分布无法唯一反推出物光场这也是相位恢复领域存在的原因之一。做实验时千万别忽略相位信息很多现象解释不了就是因为只看了一半。6.4 先养出“方向频率”的直觉如果你只从这篇文章里带走一个东西我希望是这句空间频率等于传播方向频谱面里的一个点就是一束平面波。你什么时候看到衍射花样第一反应是“这是一次傅里叶变换”而不是去翻衍射积分公式这个学科就算正式入门了。这一篇先把平面波拆解、空间频率、夫琅禾费衍射和4F系统这条主线立住。下一篇会深入菲涅耳衍射不做远场近似时怎么处理、PSF与OTF的关系以及部分相干光场里的傅里叶光学。系列内容不追求堆公式重点是把物理图像讲顺。毕竟公式可以随时翻书图景一旦建立起来就再也丢不掉了。