Cross-Scale Pansharpening via ScaleFormer and the PanScale Benchmark¶
会议: CVPR 2026
arXiv: 2603.00543
代码: GitHub
领域: 遥感
关键词: 遥感图像融合, 跨尺度泛化, Transformer, 旋转位置编码, Pansharpening
一句话总结¶
提出首个跨尺度全色锐化数据集PanScale和评测基准PanScale-Bench,以及ScaleFormer框架——将分辨率变化重新解释为序列长度变化,通过Scale-Aware Patchify分桶采样+解耦空间-序列建模+RoPE实现跨尺度泛化。
研究背景与动机¶
领域现状:全色锐化(Pansharpening)利用高分辨率全色(PAN)图像和低分辨率多光谱(LRMS)图像融合生成高分辨率多光谱图像(HRMS),是遥感图像处理的核心任务。CNN/Transformer方法(MSDCNN、HFIN、ARConv等)已取得长足进步。
现有痛点:(i) 计算与内存瓶颈——从训练裁剪尺寸(200-256px)推理到800/1600/2000px时,Transformer显存剧增,常规GPU常在800px就OOM;(ii) 分块推理伪影——被迫分块推理引入边界不连续和明显块状伪影;(iii) 跨尺度泛化弱——单一低分辨率训练导致尺度诱导的分布偏移,亮度分布随分辨率增大显著偏移。
核心矛盾:现有数据集(PanCollection、NBU、PAirMax)仅提供有限尺度多样性和分辨率,缺乏标准化多尺度+高分辨率评估协议。
本文目标 在数据、算法、计算三个维度系统解决跨尺度全色锐化的挑战。
切入角度:将分辨率变化重新表述为序列长度变化——固定空间大小的patch作为token,仅序列长度随图像尺度线性增长。
核心 idea:用Scale-Aware Patchify引入序列轴,将空间建模与尺度建模解耦,配合RoPE实现对未见尺度的外推泛化。
方法详解¶
整体框架¶
ScaleFormer 要解决的是「训练在 200–256px 小图、推理却要在 800–2000px 大图上跑」的跨尺度泛化问题。它的核心思路是把分辨率变化重新理解成序列长度变化:固定每个 patch 的空间大小,图像变大时只是 token 序列变长。输入一张 PAN 图 \(\mathbf{P} \in \mathbb{R}^{H \times W \times 1}\) 和上采样后的 MS 图 \(\mathbf{L} \in \mathbb{R}^{H \times W \times C}\),先由 Scale-Aware Patchify 切成 5D 张量 \(\mathbf{P}_{5d} \in \mathbb{R}^{B \times T \times C \times h \times w}\)(\(T\) 是序列长度),再依次过 Single Transformer(空间域与序列域分别建模)和 Cross Transformer(PAN-MS 跨模态融合),最后回归出高分辨率多光谱图。
%%{init: {'flowchart': {'rankSpacing': 24, 'nodeSpacing': 28, 'padding': 6, 'wrappingWidth': 400, 'subGraphTitleMargin': {'top': 8, 'bottom': 16}}}}%%
flowchart TD
P["PAN 图 P (H×W×1)"] --> SAP
L["上采样 MS 图 L (H×W×C)"] --> SAP
SAP["Scale-Aware Patchify<br/>分桶采样切成 5D token 序列 (B×T×C×h×w)"] --> ST
subgraph ST["解耦空间-序列建模(Single Transformer)"]
direction TB
SPA["Spatial 自注意力<br/>patch 内建模空间关系"] --> SEQ["Sequence 自注意力 + RoPE<br/>跨 patch 建模、向未见序列长度外推"]
end
ST --> CT["Cross Transformer<br/>MS 交叉查询 PAN,注入高频空间细节"]
CT --> HEAD["回归头"]
HEAD --> OUT["HRMS 高分辨率多光谱图"]
关键设计¶
1. Scale-Aware Patchify:把分辨率泛化变成序列长度泛化
直接拿小图训练、大图推理会撞上尺度诱导的分布偏移——亮度统计量随分辨率漂移,模型没见过的尺度就崩。SAP 的做法是训练时随机采样分桶索引 \(t\) 来决定窗口大小 \(w(t)\),用 Patch-to-Sequence Tokenizer 把输入切成不同长度的 token 序列,让模型在训练阶段就见过多种有效序列长度;推理时固定窗口大小,高分辨率仅靠延长序列来承接。因为每个 token 的空间尺寸始终不变,它的均值和方差就稳定下来,不会随整图变大而漂移,这正是跨尺度外推能站住的前提。
2. 解耦空间-序列建模:让尺度建模独立于空间建模
如果空间关系和尺度变化耦在同一套注意力里,序列一长就难以泛化。这里把两者拆开:Spatial Transformer 只在每个 patch 内部建模空间关系, $\(\mathbf{f}_{i,1} = \mathbf{f}_i + SA_{spa}(LN(\mathbf{f}_i))\)$ Sequence Transformer 则在序列维度建模跨 patch 的相关性, $\(\mathbf{f}_{i+1,1} = \mathbf{f}_{i+1} + SA_{seq}(LN(\mathbf{f}_{i+1}))\)$ \(SA_{seq}\) 计算时把 batch 和空间维度合并,并注入 RoPE 编码连续的相对位置。RoPE 的好处是相对位置能平滑外推到训练时没出现过的序列长度,于是模型在 1600/2000px 这种长序列上仍能保持位置感知。
3. Cross Transformer:用交叉注意力做 PAN-MS 融合
PAN 提供高频空间细节、MS 提供光谱信息,二者要融合而不是简单相加。Cross Transformer 沿用前面解耦的结构,但把自注意力换成交叉注意力,让 MS 特征去查询 PAN 特征: $\(\mathbf{f}_{i,1}^{ms} = \mathbf{f}_i^{ms} + CA_{spa}(LN(\mathbf{f}_i^{ms}), LN(\mathbf{f}^{pan}))\)$ 这样空间细节从 PAN 注入 MS,同时整套融合仍保持逐 patch、可变序列长度的处理方式,于是也具备跨尺度一致性。
损失函数 / 训练策略¶
使用L1损失 \(\mathbf{L} = \|\mathbf{H}_{out} - \mathbf{G}\|_1\)。Adam优化器,初始学习率 \(5 \times 10^{-4}\),余弦退火衰减到 \(5 \times 10^{-8}\),500 epochs,NVIDIA 3090,32通道。
实验关键数据¶
主实验:PanScale数据集跨三个子集的平均结果¶
| 方法 | Jilin PSNR/SSIM | Landsat PSNR/SSIM | Skysat PSNR/SSIM |
|---|---|---|---|
| HFIN | 38.00/0.9698 | 40.21/0.9666 | 43.96/0.9658 |
| ARConv | 38.23/0.9697 | 39.66/0.9638 | 43.40/0.9797 |
| Pan-mamba | 35.55/0.9480 | 36.73/0.9206 | 41.39/0.9493 |
| ScaleFormer | 39.29/0.9761 | 41.04/0.9711 | 44.65/0.9827 |
ScaleFormer在所有数据集上全面领先SOTA,且在分辨率增大时性能保持稳定。
消融实验:Landsat数据集¶
| 消融配置 | 200px PSNR | 400px PSNR | 800px PSNR | 1600px PSNR |
|---|---|---|---|---|
| w/o RoPE | 40.46 | 40.95 | 40.76 | 40.69 |
| SeqT→SpaT | 40.91 | 41.30 | 40.72 | 40.51 |
| w/o SAP | 40.53 | 40.93 | 40.62 | 40.39 |
| Full Model | 40.61 | 41.37 | 41.13 | 41.03 |
各消融变体均在大分辨率上出现明显性能下降,证实每个组件对跨尺度泛化不可或缺。
关键发现¶
- 模型参数量仅0.52M(HFIN的1/4,ARConv的1/9),计算效率显著优势
- 随分辨率增大,ScaleFormer的GFLOPs和显存增长远慢于HFIN/ARConv
- ARConv在分块推理时出现严重块伪影(DDC-IoU显著下降)
- 全分辨实世界场景评估(无GT)中ScaleFormer同样保持竞争力
亮点与洞察¶
- 问题重构巧妙:将分辨率泛化转化为序列长度泛化,借用NLP/视频模型中序列建模的思想
- 计算效率突出:在参数量和GFLOPs上大幅领先SOTA,且优势随分辨率增大而扩大
- 数据集贡献:PanScale是首个覆盖3种卫星平台(0.5~15m分辨率)的跨尺度全色锐化数据集
- RoPE应用创新:将文本/视频领域的RoPE引入遥感融合任务实现尺度外推
局限与展望¶
- 仅关注全色锐化任务,对其他遥感融合任务(超光谱融合、SAR-光学融合)的泛化未验证
- SAP的分桶策略是预定义的固定窗口大小集合,自适应策略可能更优
- 仅使用L1损失,感知损失或GAN损失可能进一步提升视觉质量
- Sequence Transformer的自注意力仍为 \(O(T^2)\),超大规模输入时仍有瓶颈
相关工作与启发¶
- 传统方法(GS、IHS、GFPCA)在跨尺度场景表现差(PSNR低10+dB)
- CNN方法(MSDCNN、SFINet、MSDDN)跨尺度泛化有限
- HFIN/ARConv是当前SOTA但显存和计算瓶颈严重
- Pan-mamba用Mamba架构但性能不及Transformer方案
- FlexViT的多分辨率训练和视频生成的分桶训练策略是SAP的灵感来源
PanScale数据集详情¶
- 三个子数据集:Jilin(吉林一号,0.5~1m分辨率)、Landsat(Landsat-8,15m分辨率)、Skysat(Planet SkySat,~1m分辨率)
- 测试集设计:每个子数据集包含reduced-resolution(200×200到2000×2000)和full-resolution多尺度测试集
- 数据来源:通过Google Earth Engine (GEE)系统获取和预处理
- 评估指标:PanScale-Bench整合参考指标(PSNR/SSIM/ERGAS/Q)和无参考指标(\(D_\lambda\)/\(D_S\)/QNR)
效率优势¶
| 方法 | 参数量(M) | GFLOPs(G) |
|---|---|---|
| ARConv | 4.4147 | 38.32 |
| HFIN | 1.9836 | 46.21 |
| ScaleFormer | 0.5151 | 20.57 |
评分 ⭐¶
- 新颖性: ⭐⭐⭐⭐ — 分辨率→序列长度的重构视角新颖,SAP+RoPE组合有效
- 实验充分度: ⭐⭐⭐⭐⭐ — 三数据集+多尺度+全分辨率+消融+效率分析+可视化全覆盖
- 写作质量: ⭐⭐⭐⭐ — 图表设计优秀,Fig 1/2清晰展示问题和方案对比
- 价值: ⭐⭐⭐⭐⭐ — 数据集+基准+方法三位一体贡献,推动遥感融合领域发展