Shiro Ikeda,Takeshi Nakazato,Takashi Tsukagoshi,Tsutomu T. Takeuchi,Masayuki Yamaguchi
Abstract:We reformulate the gain correction problem of the radio interferometry as an optimization problem with regularization, which is solved efficiently with an iterative algorithm. Combining this new method with our previously proposed imaging method, PRIISM, the whole process of the self-calibration of radio interferometry is redefined as a single optimization problem with regularization. As a result, the gains are corrected, and an image is estimated. We tested the new approach with ALMA observation data and found it provides promising results.
What problem does this paper attempt to address?
### What problem does this paper attempt to solve?
This paper aims to solve the gain calibration problem in radio interferometry data, specifically for the data of the Atacama Large Millimeter/sub - millimeter Array (ALMA). The paper reformulates the gain calibration problem as an optimization problem with regularization and proposes an efficient iterative algorithm to solve this problem.
#### Main problem background
1. **Radio interferometry**:
- Radio interferometry is an important astronomical observation method used to obtain high - resolution images of celestial bodies. For example, the ALMA telescope has made important contributions in astronomy and astrophysics.
- In radio interferometry, the observation data from multiple stations are combined into complex values, called visibility, which correspond to the two - dimensional Fourier transform of the image.
2. **Importance of gain calibration**:
- The gain of each telescope changes over time, causing changes in the phase and amplitude of visibility. Calibration is the process of estimating and eliminating the influence of gain changes on the observation data, which is crucial for improving the quality of the reconstructed image.
- Self - calibration is a widely used method to reconstruct images through iterative image estimation and gain calibration. However, this method needs to be carefully handled because it is a non - convex optimization problem and is prone to getting trapped in local minima.
#### Solutions proposed in the paper
1. **Redefine gain calibration**:
- The authors reformulate the gain calibration problem as an optimization problem with regularization and propose an efficient iterative algorithm to solve it.
- By combining the previously proposed imaging method PRIISM, the entire self - calibration process is redefined as a single optimization problem with regularization.
2. **Specific form of the optimization problem**:
- Image estimation: Estimate the image \(\hat{x}\) by minimizing the negative log - likelihood function with regularization:
\[
\hat{x}=\arg\min_x\left[L_{\tilde{v}}(x,\hat{g}) + R_{\lambda_1,\lambda_2}(x)\right]\quad\text{subject to}\quad x\geq0
\]
where \(R_{\lambda_1,\lambda_2}(x)=\lambda_1\|x\|_1+\lambda_2\text{TSV}(x)\), TSV represents Total Squared Variation, which is used to promote the sparsity and smoothness of the image.
- Gain estimation: Estimate the gain \(\hat{g}\) by minimizing the loss function with regularization:
\[
\hat{g}=\arg\min_g\left[L_{\tilde{v}}(\hat{x},g)+S_{\mu_1,\mu_2}(g)\right]\quad\text{subject to}\quad\sum_l|g_{\alpha l}| = N_\alpha
\]
where \(S_{\mu_1,\mu_2}(g)=\mu_1\sum_{\alpha,l}w_{\alpha l}\|g_{\alpha l}-g_{\alpha l - 1}\|^2+\mu_2\sum_{\alpha,l}w_{\alpha l}(|g_{\alpha l}|-|g_{\alpha l - 1}|)^2\), which is used to control the time - varying of the gain.
3. **Experimental verification**:
- The authors tested the new method using ALMA observation data and found that it provided promising results, especially when dealing with data sets such as HL Tau, SDP.81 and HD 142527.
In summary, this paper, by introducing a new optimization framework, solves...