D. C. Homan,J. S. Roth,A. B. Pushkarev
Abstract:We propose the use of entropy, measured from the spatial and flux distribution of pixels in the residual image, as a potential diagnostic and stopping metric for the CLEAN algorithm. Despite its broad success as the standard deconvolution approach in radio interferometry, finding the optimum stopping point for the iterative CLEAN algorithm is still a challenge. We show that the entropy of the residual image, measured during the final stages of CLEAN, can be computed without prior knowledge of the source structure or expected noise levels, and that finding the point of maximum entropy as a measure of randomness in the residual image serves as a robust stopping criterion. We also find that, when compared to the expected thermal noise in the image, the maximum entropy of the residuals is a useful diagnostic that can reveal the presence of data editing, calibration, or deconvolution issues that may limit the fidelity of the final CLEAN map.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: in radio interferometry, how to find an optimal stopping point for the CLEAN algorithm. Although the CLEAN algorithm has been widely successful as a standard deconvolution method in radio astronomy, determining the best stopping point in its iterative process remains a challenge.
Specifically, the paper proposes using the entropy in the residual image as a diagnostic and stopping indicator. The authors point out that calculating the entropy of the residual image in the final stage of the CLEAN algorithm does not require prior knowledge of the source structure or the expected noise level, and using the maximum - entropy point as a standard for measuring the randomness of the residual image is a robust stopping criterion. In addition, compared with the expected thermal noise, the maximum entropy of the residual can also reveal data editing, calibration, or deconvolution problems, which may affect the fidelity of the final CLEAN image.
### Main contributions of the paper
1. **Entropy as a stopping criterion**: The paper proposes a new method, that is, to monitor the end of CLEAN by measuring the overall entropy of the residual image during the CLEAN process. When the entropy reaches its maximum value, it indicates that all significant source structures have been removed, and the CLEAN algorithm can be stopped at this time.
2. **No dependence on prior knowledge**: This method does not need to know the thermal noise or other sources of random uncertainty in advance, making it more general and automated.
3. **Detecting potential problems**: By comparing the maximum entropy of the residual with the expected thermal noise, potential problems in the data processing process, such as data editing errors, calibration problems, or deconvolution errors, can be found.
### Calculation of entropy
The authors use Shannon's information - theoretic entropy formula to calculate entropy:
\[ S = -\sum_{j = 1}^{N_s}\sum_{k = 1}^{N_f}p_{j,k}\ln p_{j,k} \]
where:
- \( N_s \) is the number of spatial bins,
- \( N_f \) is the number of flux bins,
- \( p_{j,k} \) is the probability that a pixel falls into a specific spatial and flux bin.
### Experimental verification
The paper has carried out multiple tests by simulating VLBA data, and the results show that the maximum value of entropy can indeed be used as a reasonable stopping criterion, and good results can be achieved whether in a large - range cleaning window or in a tight cleaning window.
In conclusion, this paper provides a novel and effective method to optimize the stopping point of the CLEAN algorithm, thereby improving the quality and reliability of radio interferometry images.