本章探讨非线性光学现象在显微成像中的应用，以及如何将这些原理整合到计算机图形学的体积渲染框架中。我们将看到，显微镜成像的许多概念——从点扩散函数到光学切片——都可以用我们已经建立的体积渲染方程来描述，而非线性效应则提供了新的对比度机制和三维分辨能力。
完成本章后，您将能够：
共聚焦显微镜通过点扫描和针孔滤波实现光学切片。考虑照明点扩散函数 $h_{ill}(\mathbf{r})$ 和探测点扩散函数 $h_{det}(\mathbf{r})$，共聚焦响应为：
\[I_{conf}(\mathbf{r}_0) = \int_V \rho(\mathbf{r}) |h_{ill}(\mathbf{r} - \mathbf{r}_0)|^2 |h_{det}(\mathbf{r} - \mathbf{r}_0)|^2 d^3\mathbf{r}\]其中 $\rho(\mathbf{r})$ 是样品的反射率或荧光分布。这个方程表明共聚焦成像是照明和探测PSF的乘积，导致改善的分辨率。
这种改善可以从频域角度理解。设 $\tilde{h}{ill}(\mathbf{k})$ 和 $\tilde{h}{det}(\mathbf{k})$ 为PSF的傅里叶变换，共聚焦光学传递函数（OTF）为：
\[\text{OTF}_{conf}(\mathbf{k}) = \text{OTF}_{ill}(\mathbf{k}) \otimes \text{OTF}_{det}(\mathbf{k})\]这种卷积扩展了系统的频率响应，使得截止频率从 $2\text{NA}/\lambda$ 提高到接近 $4\text{NA}/\lambda$。
对于荧光共聚焦显微镜，需要考虑激发和发射波长的差异：
\[I_{conf}(\mathbf{r}_0) = \int_V c(\mathbf{r}) |h_{exc}(\mathbf{r} - \mathbf{r}_0; \lambda_{exc})|^2 |h_{em}(\mathbf{r} - \mathbf{r}_0; \lambda_{em})|^2 d^3\mathbf{r}\]其中 $c(\mathbf{r})$ 是荧光团浓度，$h_{exc}$ 和 $h_{em}$ 分别是激发和发射波长下的PSF。
由于Stokes位移，通常 $\lambda_{em} > \lambda_{exc}$，导致发射PSF略宽：
\[\text{FWHM}_{em} = \text{FWHM}_{exc} \cdot \frac{\lambda_{em}}{\lambda_{exc}}\]这种波长差异对最终分辨率的影响可以通过有效PSF来量化：
\[h_{eff}(\mathbf{r}) = |h_{exc}(\mathbf{r})|^2 |h_{em}(\mathbf{r})|^2\]探测针孔在共轭焦平面上的传输函数为：
\[P(\mathbf{r}_d) = \begin{cases} 1, & |\mathbf{r}_d| \leq r_{pinhole} \\ 0, & \text{否则} \end{cases}\]经过针孔后的探测PSF变为：
\[h_{det}^{(p)}(\mathbf{r}) = \mathcal{F}^{-1}\left\{ \mathcal{F}\{h_{det}(\mathbf{r})\} \cdot \mathcal{F}\{P(\mathbf{r}/M)\} \right\}\]其中 $M$ 是系统放大率。
针孔大小的优化涉及分辨率和信号强度的权衡。定义归一化针孔半径：
\[v_p = \frac{2\pi}{\lambda} \cdot \text{NA} \cdot r_{pinhole}/M\]针孔的作用可以通过其对不同轴向位置信号的响应来量化。对于位于 $z$ 处的点源，通过针孔的积分强度为：
\[T_{pinhole}(z) = \int_{|\mathbf{r}_d| \leq r_{pinhole}} |h_{det}(\mathbf{r}_d, z)|^2 d^2\mathbf{r}_d\]归一化传输函数：
\[T_{norm}(u) = \frac{T_{pinhole}(z)}{T_{pinhole}(0)}\]其中 $u = 8\pi z \sin^2(\alpha/2)/\lambda$ 是归一化轴向坐标，$\alpha = \arcsin(\text{NA}/n)$。
最优值通常为 $v_p \approx 3.8$（对应1个Airy单位），此时：
对于不同应用，针孔大小的选择策略：
共聚焦系统的轴向响应函数为：
\[h_{conf}(z) = |h_{ill}(z)|^2 |h_{det}(z)|^2\]对于高数值孔径物镜，轴向PSF近似为：
\[h(z) \approx \text{sinc}^2\left(\frac{k n z \text{NA}^2}{2}\right)\]其中 $k = 2\pi/\lambda$，$n$ 是折射率，NA 是数值孔径。
共聚焦检测将轴向响应提升为四次方，显著改善光学切片：
\[h_{conf}(z) = \text{sinc}^4\left(\frac{k n z \text{NA}^2}{2}\right)\]更精确的分析需要考虑针孔的影响。定义归一化轴向坐标：
\[u = \frac{8\pi z}{\lambda} \sin^2\left(\frac{\alpha}{2}\right) = \frac{8\pi z}{\lambda} \cdot \frac{\text{NA}^2}{2n(n + \sqrt{n^2 - \text{NA}^2})}\]包含针孔效应的轴向响应为：
\[I_{conf}(u, v_p) = \left| \int_0^1 J_0(v_p\rho) \exp\left(\frac{iu\rho^2}{2}\right) \rho d\rho \right|^4\]其中 $J_0$ 是零阶贝塞尔函数。
定义光学切片厚度为强度降至峰值一半的全宽（FWHM）：
\[\Delta z_{conf} = \frac{0.64\lambda}{n - \sqrt{n^2 - \text{NA}^2}} \approx \frac{0.64\lambda n}{\text{NA}^2} \quad (\text{NA} \ll n)\]光学切片能力还可以通过轴向积分强度响应（ISR）来表征：
\[\text{ISR}(z) = \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} I_{conf}(x, y, z) dx dy\]对于均匀荧光平面，共聚焦系统的背景抑制比为：
\[\text{BSR} = \frac{\text{ISR}_{conf}(0)}{\int_{-\infty}^{\infty} \text{ISR}_{conf}(z) dz} \approx 3\sqrt{\pi} \cdot \frac{\text{NA}^2}{\lambda n}\]这表明背景抑制能力与NA²成正比，解释了高NA物镜在共聚焦成像中的重要性。
激光扫描通过振镜系统实现，扫描位置 $\mathbf{r}_s(t)$ 的轨迹通常为：
\[\mathbf{r}_s(t) = (A_x \sin(2\pi f_x t), A_y \sin(2\pi f_y t + \phi))\]对于光栅扫描，更常见的是锯齿波模式：
\[x_s(t) = A_x \cdot \text{sawtooth}(2\pi f_x t), \quad y_s(t) = A_y \cdot \text{step}(t, \Delta y)\]完整的3D图像通过逐层扫描获得：
\[I_{3D}(x, y, z) = \int_0^T I_{conf}(\mathbf{r}_s(t), z) \delta(\mathbf{r} - \mathbf{r}_s(t)) dt\]实际系统中需要考虑：
扫描非线性：振镜的机械响应导致边缘失真 \(\mathbf{r}_{actual}(t) = \mathbf{r}_s(t) + \delta\mathbf{r}_{nonlinear}(t)\)
振镜的传递函数通常表现为二阶系统： \(H_{mirror}(s) = \frac{\omega_0^2}{s^2 + 2\zeta\omega_0 s + \omega_0^2}\)
其中 $\omega_0$ 是共振频率，$\zeta$ 是阻尼系数。
像素驻留时间：影响信噪比 \(\text{SNR} \propto \sqrt{\tau_{dwell} \cdot I_{signal}}\)
对于光子计数检测，信号遵循泊松统计： \(P(n) = \frac{(\lambda \tau_{dwell})^n e^{-\lambda \tau_{dwell}}}{n!}\)
其中 $\lambda$ 是光子到达率。
光漂白效应：累积光剂量 \(\frac{dc}{dt} = -k_{bleach} \cdot I_{exc} \cdot c\)
考虑扫描速度的影响，有效漂白率为： \(k_{eff} = k_{bleach} \cdot \frac{\tau_{dwell}}{\tau_{pixel}} \cdot \frac{I_{exc}}{I_{sat}}\)
其中 $I_{sat}$ 是饱和强度。
扫描模式优化：
扫描效率定义为： \(\eta_{scan} = \frac{\text{有效成像时间}}{\text{总扫描时间}} = \frac{N_x N_y \tau_{dwell}}{T_{frame}}\)
共聚焦检测通过互相关提取焦平面信息：
\[I_{conf} = \langle I_{ill}(\mathbf{r}) \cdot I_{det}(\mathbf{r}) \rangle\]这在频域表现为传递函数的扩展：
\[\text{OTF}_{conf}(\mathbf{k}) = \text{OTF}_{ill}(\mathbf{k}) \otimes \text{OTF}_{det}(\mathbf{k})\]导致有效截止频率提高：
\[k_{cutoff}^{conf} = 2 \cdot k_{cutoff}^{conv} = \frac{4\pi \cdot \text{NA}}{\lambda}\]共聚焦成像的信息理论分析揭示了其优势的本质。定义Fisher信息矩阵：
\[F_{ij} = \sum_\mathbf{r} \frac{1}{\sigma^2(\mathbf{r})} \frac{\partial \mu(\mathbf{r})}{\partial \theta_i} \frac{\partial \mu(\mathbf{r})}{\partial \theta_j}\]其中 $\mu(\mathbf{r})$ 是期望信号，$\theta_i$ 是待估参数（如位置、强度）。
对于定位精度，Cramér-Rao下界给出：
\[\sigma_{loc}^2 \geq \frac{1}{F_{pos}} = \frac{\sigma^2}{N} \cdot \frac{\int |h(\mathbf{r})|^2 d^3\mathbf{r}}{\int |\nabla h(\mathbf{r})|^2 d^3\mathbf{r}}\]| 共聚焦系统由于PSF更尖锐（$h_{conf} = | h | ^4$），梯度更陡： |
导致更好的定位精度。
从体积渲染角度，共聚焦成像实现了空间-角度耦合的解耦：
\[L_{conf}(\mathbf{x}, \boldsymbol{\omega}) \approx L(\mathbf{x}) \cdot \delta(\boldsymbol{\omega} - \boldsymbol{\omega}_0)\]这种方向选择性源于：
这解释了为什么共聚焦系统能有效抑制散射光——散射改变了光的传播方向，从而被针孔拒绝。
根据标量衍射理论，显微镜的3D PSF为：
\[h(\mathbf{r}) = \left| \int_{\Omega} A(\mathbf{k}_\perp) e^{i(\mathbf{k}_\perp \cdot \mathbf{r}_\perp + k_z z)} d^2\mathbf{k}_\perp \right|^2\]| 其中孔径函数 $A(\mathbf{k}_\perp)$ 由物镜的数值孔径决定，$k_z = \sqrt{k^2 - | \mathbf{k}_\perp | ^2}$ 保证波矢量满足 $ | \mathbf{k} | = k = 2\pi n/\lambda$。 |
对于圆形孔径：
\[A(\mathbf{k}_\perp) = \begin{cases} 1, & |\mathbf{k}_\perp| \leq k \cdot \text{NA} \\ 0, & \text{否则} \end{cases}\]这个孔径函数在实空间对应于：
\[a(\mathbf{r}_p) = \frac{2\pi k \text{NA}}{|\mathbf{r}_p|} J_1(k \text{NA} |\mathbf{r}_p|)\]其中 $\mathbf{r}_p$ 是瞳孔平面坐标，$J_1$ 是一阶贝塞尔函数。
使用Debye-Wolf积分，PSF的精确表达式为：
\[h(\mathbf{r}) = \left| \int_0^{\alpha} \int_0^{2\pi} \sqrt{\cos\theta} \, e^{ikr\sin\theta\cos(\phi-\phi_0)} e^{ikz\cos\theta} \sin\theta \, d\phi \, d\theta \right|^2\]其中 $\alpha = \arcsin(\text{NA}/n)$ 是最大收集角。
这个积分可以分解为三个分量（Richards-Wolf公式）：
\[\mathbf{E}(\mathbf{r}) = -\frac{ikf}{2\pi} e^{ikf} \int_0^{\alpha} \int_0^{2\pi} \mathbf{a}(\theta, \phi) \sqrt{\cos\theta} e^{i\Phi(\mathbf{r}, \theta, \phi)} \sin\theta d\phi d\theta\]其中相位项： \(\Phi(\mathbf{r}, \theta, \phi) = k[r\sin\theta\cos(\phi-\phi_0) + z\cos\theta]\)
对于高NA物镜（NA > 0.7），必须考虑矢量特性：
\[\mathbf{h}(\mathbf{r}) = \int_{\Omega} \mathbf{A}(\mathbf{k}) e^{i\mathbf{k} \cdot \mathbf{r}} d^2\mathbf{k}_\perp\]其中 $\mathbf{A}(\mathbf{k})$ 包含偏振信息。
矢量PSF的三个偏振分量为：
\[I_x = |E_x|^2, \quad I_y = |E_y|^2, \quad I_z = |E_z|^2\]总强度： \(I_{total} = I_x + I_y + I_z\)
对于线偏振入射光，轴向分量 $I_z$ 在高NA时变得显著，导致PSF的各向异性。
根据Rayleigh准则，分辨率极限为：
对于共聚焦系统，这些值改善为：
这些改善因子来自于PSF的平方关系。更精确的分析给出：
横向PSF形状（在焦平面）： \(h_\perp(r) = \left[ \frac{2J_1(v)}{v} \right]^2, \quad v = \frac{2\pi}{\lambda} \text{NA} \cdot r\)
其中 $J_1$ 是一阶贝塞尔函数。
Airy斑的第一个零点位于 $v = 3.83$，对应： \(r_{Airy} = \frac{1.22\lambda}{2\text{NA}}\)
能量分布：中心Airy斑包含总能量的84%。
轴向PSF形状（沿光轴）： \(h_\parallel(z) = \left[ \frac{\sin(u/4)}{u/4} \right]^2, \quad u = \frac{2\pi}{\lambda} n z \left(1 - \sqrt{1 - (\text{NA}/n)^2}\right)\)
对于小NA近似： \(u \approx \frac{\pi n z \text{NA}^2}{\lambda}\)
第一个零点位于 $u = 4\pi$，给出： \(z_{min} = \frac{4\lambda}{n\text{NA}^2}\)
完整3D PSF的近似（Born & Wolf）： \(h(r,z) \approx h_\perp(r) \cdot h_\parallel(z) \cdot \cos\left(\frac{vw}{4u}\right)\)
其中 $w = u \cdot (r/z)^2$，这个余弦项描述了横向和轴向的耦合。
分辨率的实用定义：
Sparrow准则（更宽松）：
FWHM准则（最常用）：
Houston准则（更严格）：
将PSF视为体积渲染中的重建核，渲染方程变为：
\[I(\mathbf{x}, \boldsymbol{\omega}) = \int_V \int_{4\pi} f(\mathbf{x}', \boldsymbol{\omega}', \boldsymbol{\omega}) h(\mathbf{x} - \mathbf{x}') L(\mathbf{x}', \boldsymbol{\omega}') d\boldsymbol{\omega}' d^3\mathbf{x}'\]这表明显微成像可以理解为带有空间变化核的体积渲染过程。
在频域中，这种关系更加清晰：
\[\tilde{I}(\mathbf{k}) = \tilde{h}(\mathbf{k}) \cdot \tilde{O}(\mathbf{k})\]其中 $\tilde{h}(\mathbf{k})$ 是光学传递函数（OTF），即PSF的傅里叶变换。
对于非相干成像，OTF是孔径函数的自相关：
\[\text{OTF}(\mathbf{k}) = \frac{\int A^*(\mathbf{u}) A(\mathbf{u} + \lambda f \mathbf{k}) d^2\mathbf{u}}{\int |A(\mathbf{u})|^2 d^2\mathbf{u}}\]其中 $f$ 是焦距。对于圆形孔径，归一化OTF为：
\[\text{OTF}(k_r) = \begin{cases} \frac{2}{\pi}\left[\cos^{-1}\left(\frac{k_r}{2k_0}\right) - \frac{k_r}{2k_0}\sqrt{1-\left(\frac{k_r}{2k_0}\right)^2}\right], & k_r \leq 2k_0 \\ 0, & k_r > 2k_0 \end{cases}\]其中 $k_0 = \text{NA}/\lambda$ 是截止频率。
体积渲染的等价性：
这种对应关系使得计算机图形学中的许多技术可以应用于显微成像：
光学切片能力可以通过积分强度响应（ISR）定量描述：
\[\text{ISR}(z) = \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} h(x, y, z) \, dx \, dy\]对于理想的光学切片，ISR应该是一个尖锐的峰。共聚焦系统的ISR为：
\[\text{ISR}_{conf}(z) = \left[ \frac{\sin(kz\text{NA}^2/4n)}{kz\text{NA}^2/4n} \right]^4\]定义切片选择性参数：
\[S = \frac{\text{ISR}(0)}{\int_{-\infty}^{\infty} \text{ISR}(z) \, dz}\]共聚焦系统的 $S_{conf} \approx 2S_{conv}$，表明其优越的背景抑制能力。
定量分析框架：
轴向响应宽度： 定义包含90%能量的轴向范围： \(\Delta z_{90\%} = 2z_0 \text{ where } \int_{-z_0}^{z_0} \text{ISR}(z) dz = 0.9 \int_{-\infty}^{\infty} \text{ISR}(z) dz\)
对于共聚焦：$\Delta z_{90\%}^{conf} \approx 2.3\lambda n/\text{NA}^2$
对比度传递函数（CTF）： 对于周期性轴向结构： \(\text{CTF}(k_z) = \frac{|\text{OTF}(0, 0, k_z)|}{\text{OTF}(0, 0, 0)}\)
共聚焦系统的轴向CTF延伸到更高频率。
切片强度（Sectioning Strength）： \(\text{SS} = \frac{d^2I/dz^2|_{z=0}}{I(0)}\)
量化焦点附近的强度变化率：
背景贡献函数： 来自深度 $z’$ 的均匀荧光层对焦点的贡献： \(B(z') = 2\pi \int_0^{\infty} h(r, z') r dr\)
总背景： \(B_{total} = \int_{-\infty}^{\infty} c(z') B(z') dz'\)
其中 $c(z’)$ 是荧光团浓度分布。
实际系统中的像差通过波前畸变影响PSF：
\[\Phi(\mathbf{k}_\perp) = \sum_{n,m} W_{nm} Z_{nm}(\rho, \theta)\]其中 $Z_{nm}$ 是Zernike多项式，$W_{nm}$ 是像差系数。常见像差包括：
球差（$W_{040}$）：深度成像时最重要 \(\Phi_{spher} = W_{040} \left( 6\rho^4 - 6\rho^2 + 1 \right)\)
折射率失配引起的球差： \(W_{040} = \frac{2\pi z}{\lambda} \left[ n_1 \sqrt{1 - (\text{NA}/n_1)^2} - n_2 \sqrt{1 - (\text{NA}/n_2)^2} \right]\)
其中 $n_1$ 是浸油折射率，$n_2$ 是样品折射率，$z$ 是成像深度。
彗差（$W_{131}$）：倾斜入射导致 \(\Phi_{coma} = W_{131} \left( 3\rho^3 - 2\rho \right) \cos\theta\)
倾斜角 $\alpha$ 引起的彗差系数： \(W_{131} = \frac{2\pi}{\lambda} \cdot \text{NA} \cdot t \cdot \sin\alpha\)
其中 $t$ 是盖玻片厚度。
像散（$W_{222}$）：非圆对称系统 \(\Phi_{astig} = W_{222} \rho^2 \cos(2\theta)\)
圆柱形表面引起的像散： \(W_{222} = \frac{\pi}{\lambda} \cdot \text{NA}^2 \cdot \Delta f\)
其中 $\Delta f$ 是两个主曲率焦距之差。
含像差的PSF为：
\[h_{aberr}(\mathbf{r}) = \left| \int_{\Omega} A(\mathbf{k}_\perp) e^{i\Phi(\mathbf{k}_\perp)} e^{i\mathbf{k} \cdot \mathbf{r}} d^2\mathbf{k}_\perp \right|^2\]像差的影响分析：
Strehl比（峰值强度比）： \(S = \frac{I_{aberr}(0,0,0)}{I_{ideal}(0,0,0)} \approx \exp(-\sigma_\Phi^2)\)
其中 $\sigma_\Phi^2$ 是波前RMS误差的方差。
Maréchal准则：$S > 0.8$ 时系统接近衍射极限。
PSF展宽： 球差导致的轴向展宽： \(\Delta z_{aberr} \approx \Delta z_{ideal} \cdot \left(1 + \frac{|W_{040}|}{\lambda}\right)\)
不对称性： 彗差引入的质心偏移： \(\Delta r_{centroid} \approx 0.28 \cdot W_{131} \cdot \frac{f}{\text{NA}}\)
自适应光学校正： 使用变形镜产生共轭相位： \(\Phi_{mirror}(\rho, \theta) = -\Phi_{aberr}(\rho, \theta)\)
校正后的残余误差： \(\sigma_{residual}^2 = \sigma_{measurement}^2 + \sigma_{fitting}^2 + \sigma_{temporal}^2\)
双光子吸收截面 $\delta_{2\gamma}$ 通过二阶微扰理论给出：
\[\delta_{2\gamma} = \frac{8\pi^2 e^4 \omega^2}{n^2 c^2 \hbar} \left| \sum_n \frac{\langle f | \mathbf{r} | n \rangle \langle n | \mathbf{r} | i \rangle}{E_n - E_i - \hbar\omega} \right|^2 g(\omega)\]其中 $g(\omega)$ 是光谱线型函数，描述两个光子的频率分布。
对于简并双光子吸收（两个相同频率的光子），跃迁速率为：
\[W_{2\gamma} = \frac{2\pi}{\hbar} \left| \sum_n \frac{\langle f | \hat{\boldsymbol{\mu}} \cdot \mathbf{E} | n \rangle \langle n | \hat{\boldsymbol{\mu}} \cdot \mathbf{E} | i \rangle}{E_n - E_i - \hbar\omega + i\Gamma_n/2} \right|^2 \rho(\omega)\]其中 $\hat{\boldsymbol{\mu}}$ 是偶极矩算符，$\Gamma_n$ 是中间态线宽。
双光子吸收截面的典型值：
双光子激发速率正比于光强的平方：
\[R_{2\gamma}(\mathbf{r}) = \frac{1}{2} \delta_{2\gamma} N(\mathbf{r}) \left(\frac{I(\mathbf{r})}{\hbar\omega}\right)^2\]对于脉冲激光，需要考虑时间分布：
\[R_{2\gamma}^{pulse} = \frac{\delta_{2\gamma} N}{2} \frac{\langle I^2 \rangle}{(\hbar\omega)^2} = \frac{\delta_{2\gamma} N}{2} \frac{g^{(2)} P_{avg}^2}{(\hbar\omega)^2 f_{rep}^2 \tau_p A^2}\]其中：
这导致固有的光学切片能力，因为激发主要发生在焦点附近。激发体积可定义为：
\[V_{2\gamma} = \frac{\left(\int I^2 dV\right)^2}{\int I^4 dV}\]对于高斯光束，$V_{2\gamma} \approx 0.3 \lambda^3/\text{NA}^3$。
双光子系统的有效PSF为：
\[h_{2\gamma}(\mathbf{r}) = |h_{exc}(\mathbf{r})|^4\]其中 $h_{exc}$ 是激发光的PSF。这导致更尖锐的轴向限制：
\[\text{FWHM}_{2\gamma}^{(z)} = \frac{\text{FWHM}_{1\gamma}^{(z)}}{\sqrt{2}}\]具体地，双光子PSF的解析形式：
横向分布： \(h_{2\gamma}^{\perp}(r) = \left[ \frac{2J_1(v)}{v} \right]^4, \quad v = \frac{2\pi}{\lambda_{exc}} \text{NA} \cdot r\)
轴向分布： \(h_{2\gamma}^{\parallel}(z) = \left[ \frac{\sin(u/4)}{u/4} \right]^4, \quad u = \frac{2\pi}{\lambda_{exc}} n z \text{NA}^2\)
分辨率改善：
双光子激发效率依赖于脉冲参数：
\[\eta_{2\gamma} \propto \frac{P_{peak}^2}{\tau_p f_{rep}} = \frac{P_{avg}^2}{\tau_p f_{rep}^2}\]优化策略：
功率控制： \(P_{sample} = P_{laser} \cdot T_{optics} \cdot \exp(-\tau_{tissue})\)
典型值：1-50 mW at sample
三光子和更高阶过程：
\[R_{n\gamma} \propto I^n\]优势：
挑战：
考虑焦深效应的体积渲染方程：
\[L(\mathbf{x}, \boldsymbol{\omega}, z_f) = \int_0^s \tau(0, t) \sigma(t) \int_V h_{DOF}(\mathbf{x}(t) - \mathbf{x}', z(t) - z_f) L_e(\mathbf{x}', \boldsymbol{\omega}) d^3\mathbf{x}' dt\]其中 $h_{DOF}$ 是依赖于离焦距离 $(z - z_f)$ 的模糊核。
对于显微镜物镜，PSF随深度变化：
\[h(\mathbf{r}, z) = \left| \int_{\Omega} A(\mathbf{k}_\perp) e^{i\Phi_{ab}(\mathbf{k}_\perp, z)} e^{i\mathbf{k} \cdot \mathbf{r}} d^2\mathbf{k}_\perp \right|^2\]其中 $\Phi_{ab}$ 包含球差等像差。
通过多个焦平面的图像融合实现扩展景深：
\[I_{EDF}(\mathbf{r}_\perp) = \sum_i w_i(\mathbf{r}_\perp, z_i) I(\mathbf{r}_\perp, z_i)\]权重函数 $w_i$ 基于局部对比度或频率内容。
景深 $\text{DOF}$ 与数值孔径的关系：
\[\text{DOF} = \frac{n\lambda}{\text{NA}^2} + \frac{n \cdot e}{M \cdot \text{NA}}\]其中 $e$ 是探测器像素尺寸，$M$ 是放大率。
薄透镜近似的弥散圆：
\[r_{CoC} = \frac{|z - z_f|}{z} \frac{f}{N(z_f - f)}\]而波动光学给出的艾里斑：
\[r_{Airy} = 1.22 \frac{\lambda f}{D} = 1.22 \lambda N\]两者在 $r_{CoC} \approx r_{Airy}$ 时过渡。
完整的显微镜成像模拟需要考虑：
相干传递函数： \(H(\mathbf{k}) = A^*(-\mathbf{k}) \otimes A(\mathbf{k})\)
部分相干照明： \(I(\mathbf{r}) = \int S(\mathbf{k}_s) |H_{coh}(\mathbf{k}, \mathbf{k}_s)|^2 |\tilde{O}(\mathbf{k})|^2 d^2\mathbf{k}_s\)
像差效应：通过Zernike多项式展开
预积分PSF查找表： \(\text{PSF}_{LUT}[z, r] = h(r, z)\)
可分离近似： \(h(\mathbf{r}, z) \approx h_\perp(r_\perp) h_\parallel(z)\)
高斯近似： \(h(\mathbf{r}, z) \approx \exp\left(-\frac{r_\perp^2}{2\sigma_\perp^2(z)} - \frac{z^2}{2\sigma_z^2}\right)\)
将显微镜像差应用于图形渲染：
球差散景： \(h_{spher}(r, \theta) = J_0(kr\sin\theta) \exp(iW_{040}k\sin^4\theta)\)
彗差效果： \(h_{coma}(r, \theta, \phi) = h_0(r, \theta) \exp(iW_{131}kr\sin\theta\cos\phi)\)
本章建立了非线性光学显微成像与计算机图形学体积渲染之间的数学联系。关键要点：
关键公式：
| 共聚焦响应：$I_{conf} = \int \rho | h_{ill} | ^2 | h_{det} | ^2 d^3\mathbf{r}$ |
| 双光子PSF：$h_{2\gamma} = | h_{exc} | ^4$ |
20.1 推导共聚焦显微镜的轴向分辨率改善因子。从单光子PSF开始： \(h(z) = \text{sinc}^2\left(\frac{k n z \text{NA}^2}{2}\right)\)
证明共聚焦PSF的FWHM约为单光子的$1/\sqrt{2}$倍。
| *提示：考虑$ | h(z) | ^4$的半高全宽。* |
20.2 计算双光子激发的局域化体积。假设高斯光束： \(I(r,z) = I_0 \frac{w_0^2}{w^2(z)} \exp\left(-\frac{2r^2}{w^2(z)}\right)\)
其中$w(z) = w_0\sqrt{1 + (z/z_R)^2}$，$z_R = \pi w_0^2/\lambda$。
提示：计算$I^2$的有效体积$V_{eff} = [\int I^2 dV]^2 / \int I^4 dV$。
20.3 证明在弱聚焦条件下（$\text{NA} < 0.5$），标量衍射理论给出的PSF可以分解为横向和轴向部分的乘积。
提示：使用Debye积分的小角度近似。
20.4 考虑具有球差的显微镜系统。球差相位为： \(\Phi_{spher}(\rho) = W_{040} \rho^4\)
| 其中$\rho = | \mathbf{k}_\perp | /(k\cdot\text{NA})$是归一化孔径坐标。推导球差对轴向PSF的影响，并分析如何通过补偿改善成像质量。 |
提示：计算含球差的PSF并展开到$W_{040}$的一阶。
20.5 推导多光子过程的一般选择定则。对于$n$光子吸收，证明跃迁必须满足： \(\Delta l = 0, \pm 2, \pm 4, ..., \pm n \text{ (偶宇称)}\) \(\Delta l = \pm 1, \pm 3, ..., \pm n \text{ (奇宇称)}\)
提示：使用角动量守恒和宇称选择定则。
20.6 设计一个算法，将显微镜的3D PSF测量数据转换为体积渲染中可用的空间变化模糊核。考虑：
提示：考虑主成分分析（PCA）或张量分解。
20.7 实现一个将普通体积渲染转换为共聚焦显微镜成像的算法框架。要求：
提示：利用FFT加速卷积运算。
20.8 开发一个模拟双光子激发过程的蒙特卡洛算法，包括：
提示：追踪光子对的时空重合。
| 正确：强度PSF = | 振幅PSF | ² |