[论文解读] Visibility Interpolation in Solar Hard X-Ray Imaging: Application to RHESSI and STIX

Emma Perracchione, Paolo Massa+2 more|arXiv (Cornell University)|Jan 1, 2021

Solar and Space Plasma Dynamics参考文献 58被引用 11

一句话总结

本文提出了一种基于可变尺度核（VSKs）的增强可见度插值方法，用于太阳硬X射线成像，以提升从稀疏(u,v)平面数据中重建图像的性能。通过利用先验源信息进行自适应核缩放，VSKs在STIX的稀疏采样和RHESSI的高分辨率数据中显著优于标准uv平滑方法，有效减少了伪影，并提高了对具有窄带或振荡性可见度模式的耀斑源的重建准确性。

ABSTRACT

Space telescopes for solar hard X-ray imaging provide observations made of sampled Fourier components of the incoming photon flux. The aim of this study is to design an image reconstruction method relying on enhanced visibility interpolation in the Fourier domain. The interpolation-based method is applied to synthetic visibilities generated by means of the simulation software implemented within the framework of the Spectrometer/Telescope for Imaging X-rays (STIX) mission on board Solar Orbiter. An application to experimental visibilities observed by the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) is also considered. In order to interpolate these visibility data, we have utilized an approach based on Variably Scaled Kernels (VSKs), which are able to realize feature augmentation by exploiting prior information on the flaring source and which are used here, for the first time, in the context of inverse problems. When compared to an interpolation-based reconstruction algorithm previously introduced for RHESSI, VSKs offer significantly better performance, particularly in the case of STIX imaging, which is characterized by a notably sparse sampling of the Fourier domain. In the case of RHESSI data, this novel approach is particularly reliable when the flaring sources are either characterized by narrow, ribbon-like shapes or high-resolution detectors are utilized for observations. The use of VSKs for interpolating hard X-ray visibilities allows remarkable image reconstruction accuracy when the information on the flaring source is encoded by a small set of scattered Fourier data and when the visibility surface is affected by significant oscillations in the frequency domain.

研究动机与目标

• 解决从(u,v)平面上稀疏采样的可见度数据中重建高保真太阳硬X射线图像的挑战。
• 克服现有插值方法（如uv平滑）的局限性，这些方法在可见度数据因窄带或分离良好的耀斑源而出现显著振荡时失效。
• 开发一种新颖的插值框架，通过特征感知的核缩放整合关于源形态的先验知识。
• 在合成STIX数据和真实RHESSI观测数据上验证该方法，证明其在性能上优于标准技术。

提出的方法

• 该方法采用可变尺度核（VSKs）在(u,v)平面上执行插值，其中核缩放由粗略的初始图像重建结果引导。
• VSK方法使用Matérn C0核以确保数值稳定性和低正则性，从而更好地处理振荡性可见度数据。
• 插值被表述为基函数的加权和，其系数通过求解一个线性系统确定，以在采样点处强制满足数据保真度。
• 在图像域应用软阈值化步骤，以抑制反傅里叶变换引起的振铃伪影。
• 该方法以两种变体应用：使用反投影（uv平滑BP）和CLEAN分量图（uv平滑CC）作为初始源估计。
• 该方法在合成STIX可见度数据和真实RHESSI数据上进行了测试，并与uv平滑和CLEAN方法进行了性能比较。

实验结果

研究问题

• RQ1 在STIX中(u,v)平面采样稀疏的情况下，基于VSK的可见度插值是否能显著提升图像重建精度？
• RQ2 在重建具有窄带、带状结构的耀斑源时，VSK插值与uv平滑和CLEAN方法相比表现如何？
• RQ3 通过核缩放整合先验源信息，在高振荡可见度区域在多大程度上能减少成像伪影？
• RQ4 对于具有细网格探测器的真实RHESSI数据，VSK方法是否在提升重建保真度的同时保持计算可行性？

主要发现

• 在STIX模拟中，基于VSK的插值显著优于标准uv平滑方法，尤其在存在强可见度振荡的配置中，得益于对稀疏采样的更好处理。
• 在具有细网格探测器的RHESSI数据中，VSK方法减少了源的碎片化并提高了重建保真度，尤其在使用探测器2–9时表现更优。
• uv平滑CC变体（使用CLEAN分量作为初始估计）达到的χ²值与CLEAN相当或更低，同时消除了手动PSF卷积的需要。
• uv平滑BP变体在源定位和通量估计方面最为准确，对于包含两个足点的配置，总通量误差控制在模拟值的5%以内。
• 计算成本相比uv平滑增加了约4–7倍，但仍与其他硬X射线成像方法保持竞争力。
• VSK方法被证明是CLEAN的用户无关替代方案，将启发式的PSF卷积步骤替换为自动化的、特征增强的插值过程。