精读笔记
Problem Setting
论文《On the stability of scale-space metrics》（arXiv preprint / 2026）关注的不是提出一个全新图像距离，而是证明一类老的 Gaussian scale-space 距离在图像比较中到底有多稳定、可算、可用于旋转不确定投影。真正问题是：在形变、投影角变化、加性噪声和有限采样同时存在时，能否有一个距离既不像 L^2 那样对局部错位过敏，又不像过强低通那样丢失判别信息。以前方法的卡点在于这些性质通常分开处理：Wasserstein 有几何稳定性但计算与噪声处理复杂；sliced 版本可算但受投影选择和噪声影响；经典 scale-space/Besov norm 有数学基础但缺少针对旋转投影比较的稳定性和算法闭环。关键矛盾是稳定性来自低频化，而判别性来自高频细节；这篇论文把这个矛盾显式交给 α、p 和尺度权重来调节。
Motivation
作者的核心观察是：Gaussian scale-space 距离其实已经位于 Besov dual norm、Wasserstein-like distance 和图像金字塔之间，只是过去没有被系统地用于“可旋转、可离散、抗噪”的图像度量。已有路线缺的是一套统一解释：为什么这个多尺度距离对小形变稳定，为什么在 r=2 时能快速做旋转对齐，为什么离散采样和噪声不会把距离破坏掉。尤其在 cryo-EM 投影比较里，图像天然有 in-plane rotation，投影方向有小角度扰动，观测又有噪声；单纯用 Euclidean 或直接 OT 都不是干净解法。论文的动机不是替代 OT 的理论地位，而是给出一个更工程可用、理论可控的 Wasserstein-adjacent metric。
Core Idea
核心思想是把函数差异不直接放在原始像素空间比较，而是放在一族 Gaussian-smoothed scale-space 表示中比较：每个尺度看一次差异，再用 τ_k^{α/2} 加权并用 ℓ^p 聚合。这样做等价于给频率轴施加一种幂律多尺度权重：低频项天然稳定，高频项只在权重允许时参与。α 越大，度量越偏低频、越稳但越可能欠分辨；α 越小，越保留细节但更像对高频噪声开放。
和 prior 的本质区别不是“用了 Gaussian smoothing”本身，而是把 Gaussian scale-space 当作一个可证明稳定的 metric family，并进一步把旋转不变计算写成频域角向相关问题。它引入的 inductive bias 很明确：图像差异应按尺度解释，小形变首先影响细尺度而不应在全局距离中被无限放大。这个 bias 比 raw L^2 更适合几何扰动，比 full Wasserstein 更容易 FFT/NUFFT 化，也比单一低通距离保留更多跨尺度信息。
Method
第一，定义 D_{α,p,r}：对 f-g 在不同 Gaussian scale τ_k 下的 L^r norm 加权求和。它解决原始距离对小位移不连续或过敏的问题；核心变化是把“点对点差异”改成“尺度相关差异”。
第二，证明与 difference-of-Gaussians 形式等价。这个部分不是主要创新，但重要，因为它说明该距离不是重复低通滤波的任意堆叠，而是与 band-pass 多尺度分解/Besov dual norm 同构；这解释了为什么 α 控制高频贡献。
第三，形变稳定性用 maximum displacement ε 控制。证明机制很直接：Gaussian kernel 在尺度 τ 上对位移的响应约为 ε/√τ；粗尺度稳定，细尺度用权重衰减截断。最终得到 α<1、α=1、α>1 三种速率，反映尺度求和的临界行为。
第四，针对 tomographic projection，引入 mixed norm 和 projection 版本的稳定性。这里真正作用是把 3D 体的局部形变/旋转投影到 2D 图像后仍可控，而不是只证明平面图像平移。
第五，旋转不变计算限定在 r=2。因为 Parseval 和 Fourier inner product 可用，旋转后的距离可以通过极坐标频域中的角向 Fourier 系数相关来一次性计算多个角度；这是算法可行性的关键，不是任意 r 都有的性质。
第六，噪声鲁棒性来自 Gaussian smoothing 后频域能量的集中界。它说明噪声项在 scale-space norm 下随 n 衰减，但衰减率由 α 决定：大 α 更强低通，噪声更快被压下。
Key Insight / Why It Works
最核心的 insight 是：几何稳定性和噪声稳定性在这个框架里来自同一个机制——高频被尺度权重软抑制。小形变主要制造高频差异；白噪声也主要在高频上均匀铺开。Gaussian scale-space 的加权求和等价于给频率差异一个可调的 low-frequency bias，因此同一个 α 同时控制 deformation sensitivity 和 noise sensitivity。
真正有效的部分很可能是频率加权形状，而不是 Gaussian 本身。Gaussian 的优势是半群性质、Fourier 形式简单、旋转下行为干净、证明容易；若换成 wavelet/Laplacian pyramid，类似机制仍可能成立。换言之，核心贡献是把“多尺度频率衰减 = 稳定 metric”这件事在图像投影和旋转对齐场景中闭合，而不是发现了新的距离现象。
α 的三段稳定率是这篇论文最值得记住的理论信号：α<1 时形变误差按 ε^α，α=1 有 log，α>1 接近线性 ε。这不是技术细节，而是说明该 metric family 的上限：想要更强稳定性就必须接受更强低频偏置。论文没有也不可能免费绕过 stability-discrimination trade-off。
旋转不变算法的增益主要是 representation alignment + FFT/NUFFT engineering，而不是新的数学稳定性。r=2 下把旋转搜索转成频域角向相关，这是经典 harmonic-analysis 思路的有效重组；实际优势很可能来自把所有候选角度的内积共享计算，而不是 metric 本身神奇。
实验中 D_2 表现更好，技术归因应谨慎：它可能主要因为更强低通抑制噪声，未必说明它在真实结构判别上更优。若任务需要细粒度高频差异，较大 α 的优势可能反转。文中未充分说明参数选择原则，因此实验增益存在 tuning/inductive-bias matching 的可能。
Relation To Prior Work
这篇论文最接近三条线：Besov/Triebel 型函数空间范数，scale-space / DoG / Laplacian pyramid 图像处理，以及 Wasserstein / sliced Wasserstein / Cramér 类几何距离。它不是从零发明 metric，而是把这些谱系连接起来，并补上稳定性、旋转不变计算和噪声离散化分析。
与 Wasserstein 的关系是：p=r=1、0<α<1 的特殊情形与某类 snowflake ground metric 的 Wasserstein 距离 metrically equivalent。但本文没有求解 OT，也不依赖 transport plan；它用 Gaussian semigroup 诱导的 multiscale norm 获得类似几何稳定性。因此它更像 Wasserstein-adjacent surrogate，而不是 OT 方法本身。
与 wavelet multiscale EMD 近似相比，本质差别是 Gaussian scale-space 更自然适配连续图像、旋转和 Fourier 计算；wavelet 更像空间-频率局部化，Gaussian 更像 isotropic smoothing semigroup。新意不在“多尺度”，而在选择了一个对旋转和投影特别友好的多尺度基底。
与 sliced Wasserstein/Cramér 相比，它不是通过随机/多方向一维投影降低维度，而是在 2D Fourier domain 中保留完整径向-角向结构，然后只对 in-plane rotation 做群不变化。这个差异很实质：sliced 方法的 inductive bias 是投影平均，本文的 inductive bias 是尺度平滑。
与经典 image pyramid/SIFT scale-space 相比，本文不是提取特征点或局部描述子，而是把整个 scale-space 作为 metric embedding。看似熟悉的 Gaussian stack 在这里被重新解释为带稳定性证明的函数距离。
Dataset / Evaluation
实验覆盖了三个合成场景：旋转/平移敏感性、加性噪声下的旋转比较、3D Gaussian mixture 的投影分类。它们确实对应论文的理论 claim：形变稳定、噪声稳定、投影角扰动下旋转不变距离可用。实验没有大段真实数据，也没有真实 cryo-EM pipeline，因此更像 controlled illustration 而不是 deployment-level validation。
benchmark 支持“该距离在合成几何扰动和噪声下行为合理”，但还不能支持“适合真实 cryo-EM class averaging”这个更强判断。真实 cryo-EM 中 CTF、非均匀噪声、mask、heterogeneity、unknown centering、particle contamination 都可能改变频率权重的最优性。文中实验中 D_2 优势明显，但它可能只是因为合成目标本身较平滑且噪声主要可由低通抑制。
没有看到系统 ablation 来分离 α、p、τ0、grid resolution、alignment resolution、NUFFT interpolation error 的贡献。增益来源不清：可能是 metric family 的结构优势，也可能是更细的角度搜索、更强低通或任务分布匹配。评价整体是合理的 proof-of-concept，但不是决定性实证证据。
Limitation
第一，方法成立依赖形变用 maximum displacement 控制，且 deformation 是质量保持型 f(φ(x))|Jφ|。这很适合 probability-density/transport 视角，但不覆盖所有图像变形模型；例如强 intensity deformation、occlusion、non-mass-preserving contrast change 不在核心理论内。
第二，稳定性上限由 α 决定，无法免费获得。α 大时噪声和形变稳定，但高频结构被弱化；α 小时更敏感但更接近 raw high-frequency comparison。论文承认参数选择是开放问题，这实际上是方法落地的核心问题，不是小缺陷。
第三，快速旋转不变算法基本锁定 r=2。r=1 与 Wasserstein 关系更直接，但旋转计算不再同样简单；因此理论上最 OT-like 的情形和算法上最方便的情形并不完全重合。
第四，离散化理论假设函数光滑、紧支撑、采样规则，且用零填充/截断控制误差。真实图像的边界、mask、非均匀采样和插值误差会让这些界变成指导性而非直接保证。
第五，实验泛化不足。当前结果不能排除“核心能力主要来自更强低通”和“合成数据频谱与 D_2 偏好匹配”。如果真实任务的判别信息在高频，D_2 的优势可能消失。
第六，旋转不变距离只处理 in-plane rotation，没有同时处理 translation 和完整 rigid motion。论文提到未来做刚体运动最小化，但这会显著增加计算和优化复杂度，也可能破坏当前漂亮的频域分解。
Takeaway
- 第一，这篇论文真正推动的是把 Gaussian scale-space metric 从“已有但分散的数学对象”变成一个可分析、可离散、可旋转对齐的图像比较工具。
- 第二，最可迁移的 insight 是：如果任务同时要求几何稳定和噪声稳定，可以直接设计频率/尺度加权的 metric，而不必总是诉诸完整 OT 或学习表示。
- 第三，α 是核心旋钮，不是普通超参数；它定义了稳定性-判别性的 Pareto 点。
- 未来真正重要的是数据自适应或任务自适应地选择尺度权重。
一句话总结
这篇论文把 Gaussian scale-space / Besov-dual 型多尺度距离系统化为一种可证明稳定、可 FFT/NUFFT 加速并适配旋转投影比较的 Wasserstein-adjacent 图像度量，真正贡献是明确了尺度权重如何控制几何稳定、噪声鲁棒与判别能力之间的 trade-off。