Geo-Information Science · Lunar Microwave Remote Sensing

基于STARFM改进的月球微波亮温数据时空融合研究

Spatio-Temporal Fusion of Lunar Microwave Brightness Temperature Data Based on An Improved STARFM Model
朱红春 黄术成 张勇 朱立德 地球信息科学学报 2026, 28(9): 2884–2897 DOI 10.12082/dqxxkx.2026.260247
Abstract · 摘要 【目的】月球凭借其稳定的表面环境、特殊的空间位置及丰富的能源潜力,已成为深空探测与空间资源开发的重要目标。月球表面处于准真空环境,太阳辐射是其主要热源,而亮温(Brightness Temperature,BT)作为关键热物理参数,能够反映月表物质的热学与介电特性。针对当前月球微波与红外亮温数据存在高空间分辨率与高时间分辨率难以兼顾、观测盲区与轨道间隙明显等问题,本研究结合月球无大气、强辐射、高真空等环境特征,对传统时空融合模型进行改进,构建适用于月球球面几何与热物理特征的时空自适应融合框架,实现多源、多分辨率、多时相亮温数据的有效融合,生成兼具高空间分辨率与高时间一致性的全月面无缝月球微波亮温数据集,为月球热环境监测、月壤介电特性反演及水冰探测等研究提供可靠的数据支撑。【方法】本文以 CE-2 MRM 月球微波亮温数据与 LRO Diviner 红外亮温数据为基础,先结合时角、太阳几何、地形等时空特征构建数理统计模型,生成高时空一致性的模拟月球微波亮温数据集;再针对月球无大气、球面几何的环境特征,对 STARFM 模型进行改进,剔除光谱差异权重、保留时空权重,并引入月球球面距离公式替代平面距离,最终实现模拟月球微波亮温与 CE-2 实测月球微波亮温的时空融合。【结果】以改进 STARFM 模型融合得到的月球微波亮温数据为估计值、以 CE-2 MRM 实测月球微波亮温数据为参考值进行定量验证,结果显示融合数据决定系数 R²≈0.90,与实测值变化趋势高度吻合;均方根误差 RMSE≈5 K,数值偏差小、反演精度高;结构相似性指数 SSIM > 0.95,空间纹理与热分布结构保持良好。在阿波罗、月球号及嫦娥系列探测器着陆区开展独立验证,融合亮温时序变化与实测亮温日变化规律之间的时序一致性误差小于 5 K,时间序列连续稳定,能够真实反映月表亮温的昼夜动态演变特征。【结论】该模型有效解决了月球微波亮温数据时空分辨率不足、覆盖不连续的问题,成功生成全月面无缝覆盖、高时空分辨率的月球微波亮温数据集。该数据集可靠性高、精度良好,可为月球热环境动态分析、月壤参数反演、水冰探测等研究提供高精度数据支撑,对月球地质勘探、着陆点优选及月球基地建设具有重要参考价值。
月球 月球微波亮温 月球红外亮温 时空融合 改进STARFM模型
Scroll · 向下滚动
Section 01 · Background & Challenges

月球亮温观测的「高空间」与「高时间」两难

月球亮温蕴含月表及浅层月壤的热学与介电信息,是月球热环境研究的关键参数。然而受卫星成像系统固有权衡的制约,任何单一传感器都难以同时实现高时间分辨率与高空间分辨率:CE-2 MRM 具备全月面微波穿透探测能力,却空间分辨率粗糙、混合像元显著、时空分布不连续;LRO Diviner 红外亮温空间分辨率高、覆盖完整,但仅响应月表极浅层辐射且存在设备周期性故障与高纬可靠性问题。与此同时,成熟的时空融合方法在地球陆地、大气场景已形成完整技术体系,在月球遥感领域的落地却十分有限。

CHALLENGE 01

微波亮温时空分辨率与覆盖不足

MRM 受窄带工作频带、有限天线孔径与仪器设计限制等固有约束,获取的数据空间分辨率粗糙、混合像元效应明显,空间分布与时间采集间隔均不连续。CE-2 单时相微波亮温数据平均覆盖度仅约 30%,观测时段分布不均,无法实现全天时均匀覆盖,观测盲区与轨道间隙明显,难以支撑连续、高精度的月表热物理行为监测。

CHALLENGE 02

既有建模策略缺乏充分验证

过往研究主要通过插值技术或经验校准方法开展全球月球亮温制图,强调辐射校正与基本作图流程,往往忽略太阳入射角、当地时间与季节效应对亮温变率的影响。自时间-角度变换概念引入以来,多数建模策略假设给定纬度上亮温变化模式在时空上一致,全球亮温模拟常依赖理想化的日亮温周期模型,这类假设缺乏地面真实值或高分辨率参考数据集的充分验证,其泛化性与准确性受到质疑。

CHALLENGE 03

地球时空融合模型难以适配月球

亮温数据通常为单通道数据,不存在高光谱意义上的光谱变化,沿用地球场景的光谱差异权重并不适用;月球无大气、强辐射、高真空,且无源微波传感器空间分辨率常在千米尺度,必须顾及月球表面曲率,用平面距离计算像元空间权重会带来偏差。具备多源协同与高精度重构能力的先进融合方法在月球遥感中的成熟方案相对匮乏。

15 kmMRM 37.0 GHz 空间分辨率
3.0 GHz 频段为 25 km
约 30%CE-2 MRM 单时相平均覆盖度
观测时段分布不均
4 频段MRM 工作频率
3.0 / 7.8 / 19.35 / 37.0 GHz
251 dCE-2 在轨飞行时长
累计采集 2 401 道亮温数据

Section 02 · Framework & Method

改进 STARFM:三步递进式时空自适应融合框架

本研究以构建兼顾月球球面几何特性与表面热物理时空异质性的多源亮温时空自适应融合框架为目标,遵循「多源数据统一处理与特征分析—时空关系建模与微波亮温模拟—改进 STARFM 时空融合与精度验证」的三步递进式技术路线。核心改进在于:针对亮温单通道特性剔除光谱差异权重、保留时空权重,并引入月球球面距离公式(月球半径 r = 1 737.1 km)替代平面距离,将像元间平面距离转换为大圆距离,从而对滑动窗口内相邻像元实现更精确的空间加权。

01

多源亮温数据统一处理与基础分析

针对 CE-2 MRM 微波亮温与 LRO Diviner 红外亮温原始数据,开展时角转换、辐射归一化处理及空间插值等预处理,统一时间基准与空间尺度;计算多源多波段月球亮温分布结果,并统计分析微波与红外亮温在时空维度上的相关关系,为后续关系建模提供基础。

时角转换辐射归一化空间插值时空基准统一
02

顾及时空特征的亮温关系建模与模拟

以 LRO Diviner Channel 7 红外亮温与 CE-2 37 GHz 微波亮温为输入,通过时角转换得到太阳方位角与入射角并与红外亮温统计匹配;采用高分辨率 LOLA DEM 提取高程、坡度、坡向等空间变量,结合月球当地时间、太阳几何等时间特征以及地表成分,建立顾及纬度、当地时间、季节效应的数理统计模型,反演生成高时空一致性的模拟微波亮温数据集,并通过定性与定量分析评估模拟精度。

LOLA DEM太阳入射角当地时间季节效应数理统计模型
03

STARFM 模型的月球化改进

在 STARFM 框架内引入多项优化:利用已知的模拟 BT(x, y, t_k) 与 t_k 时刻实测 BT,对目标时相 t_0 的预测模拟亮温进行校正,得到与 t_0 观测值一致的高分辨率亮温分布;由于亮温数据为单通道、不存在光谱变化,仅计算像元的时空权重,合并权重贡献 C = T × d;空间距离权重 d 采用月球球面距离公式计算,顾及月球表面曲率,将平面距离转换为大圆距离。

剔除光谱权重时空权重球面大圆距离r = 1 737.1 km滑动窗口
04

双输入时空融合与精度验证

将模拟微波亮温与 CE-2 实测微波亮温作为双输入,生成高时空分辨率、全月无缝覆盖的月球微波亮温融合数据集;采用交叉验证、误差分析等方法,以融合亮温为估计值、CE-2 MRM 实测亮温为参考值,通过决定系数、均方根误差与结构相似性指数进行精度验证与质量评价,并在多个历史着陆区开展时序一致性独立检验。

双输入融合交叉验证误差分析R²/RMSE/SSIM着陆区检验
SUPPLEMENT A · 数据集

CE-2 MRM 与 LRO Diviner 多源数据

CE-2 于 2010 年 10 月 1 日发射,沿距月面约 100 km 的极轨椭圆轨道飞行,轨道周期约 118 分钟,相邻轨道赤道间距约 32.7 km,回归周期接近一个月;MRM 工作于 3.0、7.8、19.35、37.0 GHz,各波段灵敏度 0.5 K,本研究采用穿透深度约 30 cm 的 37.0 GHz 通道 2c 级亮温数据(已完成几何校正与系统定标,2010 年 10 月—2011 年 5 月共 2 401 道)。Diviner 为 LRO 九通道推扫辐射计,覆盖约 0.3~400 μm,测温范围 40~400 K,沿轨空间分辨率约 320 m、垂轨约 160 m;本研究选用光谱范围 25~41 μm、信噪比高的通道 7,数据取自 2009—2016 年全球累积产品 GCP。另使用 LOLA DEM 提取地形变量。

SUPPLEMENT B · 实验设置

虹湾局部验证与全月重建配置

选取虹湾地区作为局部时空融合结果的重点展示与验证区域:该区位于月球雨海西北部,地势平缓开阔、月壤覆盖均匀,玄武岩地形可有效降低地形辐射干扰,且已有嫦娥三号着陆器实地探测数据与长时序 CE-2 微波、LRO Diviner 红外多源观测资料。模拟数据由空间分辨率 0.5° 的月球亮温模型生成,为每个 0.5°×0.5° 网格单元提供近乎完整的时序亮温信息;融合时每隔 2 h 选取一个时间点作为目标时相,虹湾区域与全月面均按该配置重建,着陆区时序验证的融合亮温每 0.5 h 输出一个数据点。


Section 03 · Results & Validation

实验结果:R²≈0.90、RMSE≈5 K 的全月无缝重建

以改进 STARFM 模型融合得到的月球微波亮温为估计值、CE-2 MRM 实测微波亮温为参考值开展定量验证,并在虹湾局部区域、全月面尺度与历史着陆区三个层次上检验模型的精度、空间可迁移性与时序一致性。

R² ≈ 0.90决定系数
多时相散点紧密聚集于 1:1 线
RMSE ≈ 5 K均方根误差
数值偏差小、反演精度高
SSIM > 0.95结构相似性指数
空间纹理与热分布结构保持良好
100%全月面覆盖率
较原始条带数据提升约 70%
数据源 空间分辨率 时间与覆盖特性 探测深度 精度表现与局限
CE-2 MRM
37.0 GHz 实测微波亮温
15 km(3.0 GHz 为 25 km) 2010 年 10 月—2011 年 5 月共 2 401 道;单时相平均覆盖度约 30%,观测时段分布不均 约 30 cm 浅层月壤 灵敏度 0.5 K;混合像元效应显著,空间分布与时间采集间隔不连续
LRO Diviner
Channel 7 红外亮温
沿轨约 320 m,垂轨约 160 m 2009—2016 年全球累积产品 GCP;全球覆盖更完整、条带范围更广 月表极浅层(毫米级)热辐射 测温范围 40~400 K;存在设备周期性故障、高纬数据可靠性差、噪声异常与时空空缺
改进 STARFM
融合微波亮温数据集
0.5°(每 0.5°×0.5° 网格单元) 时间分辨率 0.25 h,每 2 h 重建一个目标时相;全月面无缝覆盖率 100% 继承 37.0 GHz 约 30 cm 浅层探测能力 R²≈0.90、RMSE≈5 K、SSIM≥0.95;着陆区时序一致性误差小于 5 K
FINDING 01 · 虹湾局部

时空不连续性显著降低

以虹湾地区为例,每隔 2 h 选取一个目标时相重建 CE-2 实测微波亮温:融合结果在所有时间阶段均实现高分辨率的全空间覆盖,显著降低了原始 CE-2 观测中存在的时空不连续性。融合亮温的空间分布与月球地形吻合良好,在融合影像中可清晰识别虹湾地区的多个大型撞击坑,表明模型对月表热非均质细节具备良好还原能力。

FINDING 02 · 全月尺度

覆盖率提升约 70%,日循环特征清晰

将模型扩展至全球并每 2 h 重建一次,融合数据集以 0.5° 空间分辨率实现 100% 全球覆盖,较原始 CE-2 MRM 条带数据(覆盖率不足 30%)提升约 70%。空间上呈赤道高、两极低的纬度梯度,正午赤道附近亮温达 270 K,夜间下降约 50 K;时间维度上有效捕获日循环特征,正午与夜间亮温差异达 220 K,与 Diviner 红外亮温及 CE-2 实测结果一致。

FINDING 03 · 着陆区验证

绝对误差限制在 ±5 K 以内

选取 Apollo 系列(美国)、Luna 系列(前苏联)与嫦娥系列(中国)等地理与地质条件多样的历史着陆点开展独立评估,融合亮温每 0.5 h 输出一个数据点,可详细描述整个日热循环。时间剖面与经验测量的相位对准与幅度一致,即使在存在微小差异的时段,绝对误差通常也被限制在 ±5 K 以内,完全在月球微波热遥感数据的可接受范围内。

改进的 STARFM 模型通过融合多源数据,将月球微波亮温的空间覆盖率从不足 30% 提升至 100%,同时保持 R²≈0.90、RMSE≈5 K 的高精度,为月球热环境分析、月壤参数反演及未来任务规划提供了高时空分辨率的数据支撑。

— 核心发现 · 基于STARFM改进的月球微波亮温数据时空融合研究

Section 04 · Conclusion & Prospect

三大贡献、应用价值与未来议程

本研究以月球微波亮温数据的时空融合为目标,基于 LRO Diviner 模拟得到的月球微波亮温数据并结合 CE-2 MRM 实测数据,在 STARFM 方法改进的基础上融合生成了高时空分辨率的全月月球微波亮温数据集,突破了实测微波亮温在时空分辨率与覆盖范围上的限制。

CONTRIBUTION 01

构建月球化的时空融合框架

针对现有遥感融合模型无法适配月球无大气、球面几何的特殊环境问题,剔除光谱差异权重、保留时空权重,并以月球球面距离公式替代平面距离,构建面向月球亮温重建的时空自适应融合框架,充分挖掘微波与红外观测数据的互补优势。

CONTRIBUTION 02

生成全月无缝亮温数据集

实现月球微波亮温数据在 0.5° 空间分辨率与 0.25 h 时间分辨率下的全月无缝覆盖,融合数据集决定系数 R²≈0.90、均方根误差 RMSE≈5 K、结构相似性指数 SSIM≥0.95,均表现优异,时序连贯、空间无缝。

CONTRIBUTION 03

完成多区域独立精度验证

依托多区域实测与模拟数据完成精度校验:在 Apollo、Luna 及 CE 系列探测器着陆区的验证中,时序一致性验证误差小于 5 K,融合值与实测值的相关性与结构一致性均较强,凸显模型在长期月球热环境监测中的实际适用性与科学可靠性。

APPLICATION VALUE

应用价值

该数据集可靠性高、精度良好,可为月球热环境动态分析、月壤参数反演、水冰探测等研究提供高精度数据支撑,并支持月表热建模、地下结构分析等下游应用;对月球地质勘探、着陆点优选及月球基地建设具有重要参考价值,也可服务于未来着陆器或漫游车的任务规划。研究团队来自山东科技大学测绘与空间信息学院、国家卫星气象中心、许健民气象卫星创新中心等单位。

LIMITATIONS & FUTURE

局限与未来

当前工作主要基于 CE-2 MRM 37.0 GHz 单一通道与 Diviner 通道 7 数据,模型普适性仍有提升空间。未来将把融合模型进一步扩展至多源异构数据(如红外-微波联合反演),并推广至 CE-2 MRM 的其他波段,推动月球科学研究的精细化与多源化发展,为反演月壤厚度、成分等地表地质特征提供更丰富的数据支持,对加深对地球和其他天体微波遥感特征的认识具有重要意义。