Mapped Variably Scaled Kernels: Applications to Solar Imaging

Francesco Marchetti,Emma Perracchione,Anna Volpara,Anna Maria Massone,Stefano De Marchi,Michele Piana
2023-04-03
Abstract:Variably scaled kernels and mapped bases constructed via the so-called fake nodes approach are two different strategies to provide adaptive bases for function interpolation. In this paper, we focus on kernel-based interpolation and we present what we call mapped variably scaled kernels, which take advantage of both strategies. We present some theoretical analysis and then we show their efficacy via numerical experiments. Moreover, we test such a new basis for image reconstruction tasks in the framework of hard X-ray astronomical imaging.
Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to improve the accuracy and stability of kernel - based interpolation methods when dealing with data with discontinuities and complex distributions, especially for the application of solar hard X - ray imaging. Specifically, the author introduced a new method - Mapped Variably Scaled Kernels (MVSKs), which combines the advantages of Variably Scaled Kernels (VSKs) and the False Node Approach (FNA) to improve the effects of function interpolation and image reconstruction. ### Problem Background 1. **Limitations of Traditional Kernel Interpolation**: - In many cases, it is difficult to construct a basis function that is suitable for the data distribution and the target function only by adjusting the shape parameter. - For functions with jump discontinuities, traditional interpolation methods are prone to Gibbs and Runge phenomena, resulting in inaccurate interpolation results. 2. **Deficiencies of Existing Methods**: - Although Variably Scaled Kernels (VSKs) can improve the flexibility of interpolation to a certain extent, there are still challenges when dealing with complex discontinuous functions. - The False Node Approach (FNA) can improve the condition number of interpolation, but it may not be efficient enough when dealing with large - scale data. ### Solutions Proposed in the Paper The author proposed a new method - Mapped Variably Scaled Kernels (MVSKs), aiming to combine the advantages of VSKs and FNA to improve the performance of interpolation and image reconstruction. Specific improvements include: 1. **Introducing the Mapping Function S**: Map the original data to a new space through the mapping function S, thereby improving the node distribution, reducing the fill distance and increasing the separation distance, and further improving the stability of interpolation. 2. **Defining the Shape Function ψ**: The shape function ψ is used to simulate the jump discontinuities of the target function and prevent the occurrence of Gibbs phenomenon, especially when dealing with discontinuous functions. 3. **Theoretical Analysis and Numerical Experiments**: The paper provides a theoretical analysis of MVSKs and verifies its effectiveness in dealing with solar hard X - ray imaging data through numerical experiments. ### Application Scenarios This method is especially suitable for solar hard X - ray imaging, especially the data of ESA's STIX telescope. By using MVSKs, the author shows that when dealing with these data, the image can be reconstructed more accurately, artifacts can be reduced, and the image quality can be improved. ### Summary The main contribution of the paper is to propose a new interpolation method - Mapped Variably Scaled Kernels (MVSKs), which combines the advantages of VSKs and FNA and effectively solves the limitations of traditional interpolation methods when dealing with complex data. Through theoretical analysis and practical application verification, the superior performance of MVSKs in solar hard X - ray imaging has been proven.