分辨率与傅里叶光学

Learning objectives / 学习目标

1. 在明确 circular aperture 与 Rayleigh convention 的前提下计算 telescope、camera 与 microscope 的 diffraction-limited resolution。
2. 把一维或二维图样分解为 spatial-frequency components，并由频率矢量判断周期与条纹方向。
3. 推导 paraxial lens 后焦面坐标与 object spatial frequency 的映射。
4. 分清 coherent field convolution 与 incoherent intensity convolution，正确使用 ATF、OTF 与 MTF。
5. 设计 4-f spatial filter，预测 low-pass、high-pass 与 notch filtering 的作用和代价。

Prerequisites / 先修知识

• Single-aperture diffraction：Airy pattern 与 Fraunhofer diffraction。
• Multiple slits and gratings：diffraction orders 与 grating equation。
• 复振幅 \(U=Ae^{i\phi}\)、强度 \(I\propto\lvert U\rvert^2\)。
• Fourier transform、inverse transform 与 convolution theorem。

Scope and conventions / 范围与约定

• 默认 scalar、monochromatic、paraxial optics；若无特别说明，取空气中 \(\lambda=\lambda_0\)。
• 孔径直径记为 \(D_a\)，焦距记为 \(f\)，f-number 为 \(N_f=f/D_a\)。
• Fourier transform 约定为

\[\mathcal F\{u\}(f_x,f_y) =\iint u(x,y)e^{-i2\pi(f_xx+f_yy)}\,dx\,dy.\]

• 透镜后焦面坐标 obey

\[f_x=\frac{X}{\lambda f}, \qquad f_y=\frac{Y}{\lambda f}.\]

• coherent imaging 线性作用于 complex field；incoherent imaging 线性作用于 intensity。二者的 PSF 与 transfer function 不可混用。
• “diffraction limit”必须连同 aperture shape、wavelength、criterion、coherence 与 medium 一起陈述；Rayleigh criterion 是实用约定，不是信息论的绝对界限。
• 所有器件图均为 schematic，not to scale。

Notation and modeling contract / 符号与建模契约

Knowledge 1 — Rayleigh criterion and diffraction limit / 瑞利判据与衍射极限

Knowledge ID: P05-U05-L03-K01

理想 circular unobstructed aperture 对 point source 产生 Airy pattern。第一暗环的小角度位置为

\[\theta_1\approx1.22\frac{\lambda}{D_a}.\]

Rayleigh criterion 规定：一个 central maximum 落在另一个 pattern 的第一 minimum 处时，两点“刚好可分辨”。因此

\[\boxed{\theta_R=1.22\frac{\lambda}{D_a}}.\]

焦平面中的第一零点半径为

\[r_1=f\theta_R =1.22\frac{\lambda f}{D_a} =1.22\lambda N_f.\]

图 1 仅展示“峰—第一零点”的判据关系；曲线不是精确 Airy profile。

图 2 给出极限检查：\(D_a\to\infty\) 时 \(\theta_R\to0\)；\(\lambda\) 增大时曲线整体上移。

显微镜常用 object-space numerical aperture

\[\mathrm{NA}=n\sin\alpha.\]

incoherent point imaging 的 lateral Rayleigh separation 常写成

\[\boxed{\delta_R\approx0.61\frac{\lambda_0}{\mathrm{NA}}}.\]

\(\lambda_0\) 是 vacuum wavelength，介质折射率已经包含在 NA 中。aberration、sampling、noise 与 specimen contrast 只会使实际任务更复杂。

Knowledge 2 — Spatial-frequency representation / 空间频率表示

Knowledge ID: P05-U05-L03-K02

一维 sinusoidal object 可写为

\[g(x)=g_0\left[1+m\cos(2\pi f_xx+\phi)\right],\]

其中 period \(\Lambda=1/\lvert f_x\rvert\)。二维分量写成

\[g(\boldsymbol r) =A\cos\left(2\pi\boldsymbol f\cdot\boldsymbol r+\phi\right), \qquad \boldsymbol f=(f_x,f_y).\]

\(\boldsymbol f\) 指向 phase 增长最快方向，因而垂直于等相位条纹；\(\lvert\boldsymbol f\rvert\) 决定条纹间距。

图 3 表明 high spatial frequency 表示更细密的变化，不是“更亮”或“更高的 optical temporal frequency”。

图 4 同时编码 orientation 与 period。\(\boldsymbol f=\boldsymbol0\) 是 DC，给出平均亮度；低频描述缓慢轮廓，高频描述细节与 edge。real-valued object 的 spectrum 具有 \(\tilde g(-\boldsymbol f)=\tilde g^*(\boldsymbol f)\)。

Knowledge 3 — Lens as a Fourier transformer / 透镜作为傅里叶变换器

Knowledge ID: P05-U05-L03-K03

把 field \(U_o(x,y)\) 放在 ideal thin lens 的 front focal plane；在 scalar paraxial Fresnel approximation 下，back focal plane field 为

\[U_f(X,Y) =C\,e^{i\Phi(X,Y)} \tilde U_o\left( \frac{X}{\lambda f}, \frac{Y}{\lambda f} \right),\]

其中 \(C\) 是整体因子，\(\Phi\) 是已知 phase factor。可测的核心映射为

\[\boxed{X=\lambda f f_x,\qquad Y=\lambda f f_y}.\]

图 5 中 front/back focal-plane placement 是关键；输入场整体映射到频谱，不是“一个物点对应一个 Fourier 点”的几何成像。任意物距下仍可能有 scaled transform，但伴随额外 propagation 和 phase。

若输入分量 \(U\propto e^{i2\pi f_xx}\)，它等价于小角度 plane wave：

\[f_x=\frac{\sin\theta_x}{\lambda} \approx\frac{\theta_x}{\lambda}, \qquad X=f\tan\theta_x\approx\lambda f f_x.\]

图 6 将 grating orders、angular spectrum 和 focal-plane spots 统一：每个方向分量在后焦面占据一个位置。

Knowledge 4 — Transfer functions, convolution, and image formation / 传递函数、卷积与成像

Knowledge ID: P05-U05-L03-K04

在 isoplanatic region 内，coherent system 对 complex field 线性：

\[U_i=U_o*h_c, \qquad \tilde U_i=H_c\,\tilde U_o.\]

\(h_c\) 是 amplitude impulse response，\(H_c=\mathcal F\{h_c\}\) 是 coherent ATF。理想 pupil 的 \(H_c\) 在尺度映射后与 pupil function 对应。

互不相干 object points 在 intensity 上相加：

\[I_i=I_o*h_I, \qquad \tilde I_i=\mathrm{OTF}\,\tilde I_o,\]

\[\mathrm{MTF}(\boldsymbol f) =\frac{\lvert\mathrm{OTF}(\boldsymbol f)\rvert} {\mathrm{OTF}(\boldsymbol0)}.\]

理想 incoherent OTF 是 pupil 的 normalized autocorrelation，而不是把 coherent pupil 直接改名。

图 7 是 real-space view；frequency domain 中同一过程成为逐频率乘法。

图 8 以同一 1D pupil 为基准：coherent amplitude support 截至 \(f_{c,\mathrm{coh}}=\mathrm{NA}/\lambda\)，其 autocorrelation 对应的 incoherent OTF support 延伸到 \(2f_{c,\mathrm{coh}}\)。circular 2D pupil 的 exact incoherent MTF 不是三角线。

Knowledge 5 — Spatial filtering and coherent imaging / 空间滤波与相干成像

Knowledge ID: P05-U05-L03-K05

两片 focal length 均为 \(f\) 的 lenses 可组成 4-f processor。第一片在 common focal plane 形成 input spectrum；mask \(H_F(f_x,f_y)\) 乘上该 spectrum；第二片 inverse transforms：

\[U_{\mathrm{out}}(x,y) \propto\mathcal F^{-1} \left[ H_F(f_x,f_y)\tilde U_{\mathrm{in}}(f_x,f_y) \right].\]

图 9 的 filter plane 是 physical Fourier plane；上方箭头表示 field transform 与 spectrum multiplication，不是会在 mask 处折转的几何光线。mask 坐标要由 \(X=\lambda f f_x\) 换算。

• circular low-pass aperture：保留 DC/低频，去除细节并拓宽 PSF；
• central stop/high-pass：抑制均匀背景、强调 edges，但会出现 ringing；
• directional notch：去除一对特定频率，可抑制 periodic artifact；
• phase mask：改变 phase，可把 phase variation 转成 measurable intensity variation。

图 10 强调 filter 是 spectrum multiplier。频谱删去后通常不可逆；hard mask edge 会在 real space 产生 ringing。coherent imaging 保留 phase，因此 dust、surface error 与 parasitic reflections 还会产生 interference fringes 和 speckle。

Worked Examples / 完整例题

Worked Example 1 — Telescope Rayleigh separation

\(D_a=0.150\ \mathrm{m}\) 的 circular telescope 观测 \(\lambda=550\ \mathrm{nm}\) 的两个远处点源。求 Rayleigh angle；若距离 \(L=1.00\ \mathrm{km}\)，求 transverse separation。

已知： \(D_a=0.150\ \mathrm{m}\)，\(\lambda=5.50\times10^{-7}\ \mathrm{m}\)。
研究对象与系统： entrance pupil 与两个 incoherent point sources；忽略 atmosphere、aberration 与 detector。
参考系与正向： optical axis 为 \(+z\)，transverse separation 取正 magnitude。
模型与理由： far field、small angle、ideal circular aperture，使用 Rayleigh criterion。

\[\theta_R=1.22\frac{5.50\times10^{-7}}{0.150} =4.47\times10^{-6}\ \mathrm{rad},\]

\[s_R\approx L\theta_R =(1.00\times10^3)(4.47\times10^{-6}) =4.47\times10^{-3}\ \mathrm{m}.\]

结果： \(\theta_R=4.47\ \mathrm{\mu rad}\)，\(s_R=4.47\ \mathrm{mm}\)（3 s.f.）。
检查与解释： \(\lambda/D_a\) 无量纲，\(L\theta\) 为 m；diameter 翻倍时两值减半。这只是 ideal diffraction baseline。

Worked Example 2 — Microscope Rayleigh limit

oil-immersion objective 的 \(\mathrm{NA}=1.30\)，\(\lambda_0=520\ \mathrm{nm}\)。估算 lateral Rayleigh separation。

已知： \(\lambda_0=5.20\times10^{-7}\ \mathrm{m}\)，\(\mathrm{NA}=1.30\)。
研究对象与系统： specimen 中两个 incoherent point emitters 与 corrected objective；NA 已包含 medium index。
参考系与正向： lateral distance 只取 magnitude。
模型与理由： scalar diffraction 与 Rayleigh convention。

\[\delta_R=0.61\frac{5.20\times10^{-7}}{1.30} =2.44\times10^{-7}\ \mathrm{m}.\]

结果： \(\delta_R=244\ \mathrm{nm}\)（3 s.f.）。
检查与解释： NA 无量纲；NA 增大则 distance 下降。不可把 refractive index 再除一次；noise、contrast、sampling 可能恶化实测值。

Worked Example 3 — Grating frequency to Fourier-plane position

\(\Lambda=20.0\ \mathrm{\mu m}\) grating 放在 \(f=0.200\ \mathrm{m}\) lens 前焦面，以 \(\lambda=633\ \mathrm{nm}\) coherent light 正入射。求 \(\pm1\) orders 的位置。

已知： \(f_x=\pm1/\Lambda\)，\(f=0.200\ \mathrm{m}\)，\(\lambda=6.33\times10^{-7}\ \mathrm{m}\)。
研究对象与系统： grating—lens—back focal plane。
参考系与正向： variation 沿 \(x\)，焦面右侧为 \(+X\)。
模型与理由： scalar paraxial Fourier transform；\(\lambda/\Lambda=0.0317\ll1\)。

\[f_x=\pm5.00\times10^4\ \mathrm{m^{-1}},\]

\[X=\lambda f f_x =\pm6.33\times10^{-3}\ \mathrm{m}.\]

结果： \(X=\pm6.33\ \mathrm{mm}\)（3 s.f.）。
检查与解释： \((\mathrm m)(\mathrm m)(\mathrm m^{-1})=\mathrm m\)；real grating 给出 symmetric conjugate orders。更高 spatial frequency 聚焦得更远。

Worked Example 4 — Fourier-plane low-pass cutoff

4-f coherent processor 采用 \(f=0.100\ \mathrm{m}\) lenses 和 \(\lambda=500\ \mathrm{nm}\) light。Fourier plane 的 centered circular aperture 半径为 \(a_F=0.500\ \mathrm{mm}\)。求 cutoff radial frequency 与 minimum transmitted period。

已知： \(a_F=5.00\times10^{-4}\ \mathrm{m}\)、\(\lambda=5.00\times10^{-7}\ \mathrm{m}\)、\(f=0.100\ \mathrm{m}\)。
研究对象与系统： coherent input field spectrum 与 hard low-pass mask。
参考系与正向： mask center 在 optical axis，radial frequency 非负。
模型与理由： ideal 4-f mapping \(R_F=\lambda f\rho_f\)。

\[\rho_{f,c} =\frac{a_F}{\lambda f} =1.00\times10^4\ \mathrm{m^{-1}},\]

\[\Lambda_{\min} =\frac{1}{\rho_{f,c}} =1.00\times10^{-4}\ \mathrm{m}.\]

结果： \(10.0\ \mathrm{cycles\,mm^{-1}}\)，\(\Lambda_{\min}=0.100\ \mathrm{mm}\)（3 s.f.）。
检查与解释： \(a_F/(\lambda f)\) 为 m\(^{-1}\)；\(a_F\to\infty\) 时 cutoff 增大。更细的周期被压制，hard edge 还会引入 ringing。

Lesson synthesis / 全课综合

aperture 决定 point-source PSF；PSF 重叠形成 resolution criterion；arbitrary object 又可分解为空间频率。lens 把 angle/spatial frequency 映射到 physical focal-plane position；pupil 与 aberration 对每个分量施加 transfer magnitude/phase；inverse transform 或 convolution 得到 image。Fourier-plane mask 可以有意修改 spectrum，但会同步改变 PSF 并丢失信息。

从系统设计角度看，“分辨率”不是某一个元件的孤立属性。入口孔径先限制可收集的角谱范围，像差再改变孔径上的相位，探测器则以有限像素对最终强度采样。即使理想光学系统仍能传递某一空间频率，若该频率接近探测器奈奎斯特频率，或者其调制度已经低于噪声造成的不确定度，实际图像仍可能无法可靠辨认。因此完整的分辨率声明应同时给出波长、孔径、相干性、视场位置、判据、采样和信噪比。

实空间和频率空间并不是两套互相竞争的描述。实空间中的一个窄点扩散函数，对应频率空间中较宽的传递带宽；频率空间的锐利截断，则对应实空间中具有长尾和旁瓣的脉冲响应。选择在哪个空间计算，应由问题结构决定：局部点源、边缘扩散和卷积更适合实空间；周期结构、孔径截止和滤波器设计更适合频率空间。最终结果必须能在两种描述之间互相核对。

还要区分“观察到频谱”与“重建出物体”。后焦面上某个亮点说明相应方向分量到达了该平面，但其亮度只给出强度信息；若物体的相位未知，仅测量焦面强度通常不足以唯一重建复振幅。相位恢复需要额外传播距离、参考波或先验约束。这个信息边界正是相干成像比几何光学多出的关键内容。

Formula-and-decision table / 公式与决策表

Model boundaries / 模型边界

1. high-NA/vector polarization coupling 需要 vector diffraction。
2. large angle 下不能使用 \(\sin\theta\approx\tan\theta\approx\theta\)。
3. broadband system 的 PSF、Fourier-plane scale 与 focus 都随 \(\lambda\) 改变。
4. wide field、vignetting 或 field-dependent aberration 下，PSF 依赖 object position，LSI 失效。
5. partially coherent illumination 不能强套 coherent/incoherent 任一极限。
6. obstruction、apodization 与 segmented pupil 会改变 Airy rings 与 transfer function。
7. optical cutoff 之外还要检查 detector Nyquist、shot noise 与 read noise。

Misconception clinic / 误区诊所

Practice and retrieval / 练习与检索

1. [K01 概念] circular aperture diameter 加倍，固定 wavelength 时 Rayleigh angle 如何变化？

2. [K01 定量] \(D_a=2.40\ \mathrm{m}\)、\(\lambda=600\ \mathrm{nm}\) 时求 \(\theta_R\)。

3. [K02 图像] 一组竖直条纹的 spatial-frequency vector 指向哪里？

4. [K02 定量] period 为 \(0.250\ \mathrm{mm}\) 的条纹频率是多少？

5. [K03 推导] 从 \(f_x=\sin\theta/\lambda\) 与 \(X=f\tan\theta\) 推出 paraxial mapping。

6. [K03 定量] \(\lambda=532\ \mathrm{nm}\)、\(f=150\ \mathrm{mm}\) 时，\(20.0\ \mathrm{cycles\,mm^{-1}}\) 分量在哪里？

7. [K04 判断] coherent illumination 下可否直接写 \(I_i=I_o*h_c\)？

8. [K04 综合] 两级 incoherent systems 在某频率的 MTF 为 0.80、0.75，detector MTF 为 0.90。忽略 phase reversal，total MTF 为何？

9. [K05 设计] coherent 4-f system 中，一对离散 spectrum spots 对应 periodic noise。应选何种 filter？

10. [K05 边界] hard low-pass stop 为什么在 edge 附近产生 ringing？

Exit ticket / 离堂检验

1. 设计实验：用已知-period grating 标定 lens Fourier-plane scale，写出测量量、方程与误差来源。
2. 对“细节消失”的 image，区分 diffraction cutoff、aberration MTF、detector sampling 与 noise，并说明诊断 measurement。

References and source notes / 参考资料与来源说明

所有链接于 2026-08-28 实际打开并核对；正文独立组织，图均依据方程原创绘制。

1. OpenStax University Physics Vol. 3, §4.5 Circular Apertures and Resolution。
2. MIT OCW 2.71, Supplement to Lecture 18。
3. MIT OCW 2.71, Lecture 18 transcript。
4. MIT OCW 2.71 (2014), Lecture 16 Wave Optics。
5. MIT OCW 2.71, Lecture 22。
6. MIT OCW 2.71, Lecture 23 transcript。
7. University of Iowa, Abbe Theory and Optical Fourier Transform。

Source-to-content map

复习与拓展

• 与 Imaging Quality and Aberrations 联结：pupil phase error 如何改变 PSF、OTF 与 Strehl ratio。
• 与 Polarization States 联结：high-NA/vector imaging 中 polarization 与 resolution 的耦合。
• 后续可扩展：Abbe coherent imaging、Hopkins partially coherent theory、phase contrast、synthetic aperture 与 regularized deconvolution。