-
Superresolution Interferometric Imaging with Sparse Modeling Using Total Squared Variation --- Application to Imaging the Black Hole Shadow
Authors:
Kazuki Kuramochi,
Kazunori Akiyama,
Shiro Ikeda,
Fumie Tazaki,
Vincent L. Fish,
Hung-Yi Pu,
Keiichi Asada,
Mareki Honma
Abstract:
We propose a new superresolution imaging technique for interferometry using sparse modeling, utilizing two regularization terms: the $\ell_1$-norm and a new function named Total Squared Variation (TSV) of the brightness distribution. TSV is an edge-smoothing variant of Total Variation (TV), leading to reducing the sum of squared gradients. First, we demonstrate that our technique may achieve super…
▽ More
We propose a new superresolution imaging technique for interferometry using sparse modeling, utilizing two regularization terms: the $\ell_1$-norm and a new function named Total Squared Variation (TSV) of the brightness distribution. TSV is an edge-smoothing variant of Total Variation (TV), leading to reducing the sum of squared gradients. First, we demonstrate that our technique may achieve super-resolution of $\sim 30$% compared to the traditional CLEAN beam size using synthetic observations of two point sources. Second, we present simulated observations of three physically motivated static models of Sgr A* with the Event Horizon Telescope (EHT) to show the performance of proposed techniques in greater detail. We find that $\ell_1$+TSV regularization outperforms $\ell_1$+TV regularization with the popular isotropic TV term and the Cotton-Schwab CLEAN algorithm, demonstrating that TSV is well-matched to the expected physical properties of the astronomical images, which are often nebulous. Remarkably, in both the image and gradient domains, the optimal beam size minimizing root-mean-squared errors is $\lesssim 10$% of the traditional CLEAN beam size for $\ell_1$+TSV regularization, and non-convolved reconstructed images have smaller errors than beam-convolved reconstructed images. This indicates that the traditional post-processing technique of Gaussian convolution in interferometric imaging may not be required for the $\ell_1$+TSV regularization. We also propose a feature extraction method to detect circular features from the image of a black hole shadow with the circle Hough transform (CHT) and use it to evaluate the performance of the image reconstruction. With our imaging technique and the CHT, the EHT can constrain the radius of the black hole shadow with an accuracy of $\sim 10-20$% in present simulations for Sgr A*.
△ Less
Submitted 13 March, 2018; v1 submitted 15 February, 2018;
originally announced February 2018.
-
Imaging the Schwarzschild-radius-scale Structure of M87 with the Event Horizon Telescope using Sparse Modeling
Authors:
Kazunori Akiyama,
Kazuki Kuramochi,
Shiro Ikeda,
Vincent L. Fish,
Fumie Tazaki,
Mareki Honma,
Sheperd S. Doeleman,
Avery Broderick,
Jason Dexter,
Monika Mościbrodzka,
Katherine L. Bouman,
Andrew Chael,
Masamichi Zaizen
Abstract:
We propose a new imaging technique for radio and optical/infrared interferometry. The proposed technique reconstructs the image from the visibility amplitude and closure phase, which are standard data products of short-millimeter very long baseline interferometers such as the Event Horizon Telescope (EHT) and optical/infrared interferometers, by utilizing two regularization functions: the…
▽ More
We propose a new imaging technique for radio and optical/infrared interferometry. The proposed technique reconstructs the image from the visibility amplitude and closure phase, which are standard data products of short-millimeter very long baseline interferometers such as the Event Horizon Telescope (EHT) and optical/infrared interferometers, by utilizing two regularization functions: the $\ell_1$-norm and total variation (TV) of the brightness distribution. In the proposed method, optimal regularization parameters, which represent the sparseness and effective spatial resolution of the image, are derived from data themselves using cross validation (CV). As an application of this technique, we present simulated observations of M87 with the EHT based on four physically motivated models. We confirm that $\ell_1$+TV regularization can achieve an optimal resolution of $\sim 20-30$% of the diffraction limit $λ/D_{\rm max}$, which is the nominal spatial resolution of a radio interferometer. With the proposed technique, the EHT can robustly and reasonably achieve super-resolution sufficient to clearly resolve the black hole shadow. These results make it promising for the EHT to provide an unprecedented view of the event-horizon-scale structure in the vicinity of the super-massive black hole in M87 and also the Galactic center Sgr A*.
△ Less
Submitted 23 February, 2017;
originally announced February 2017.
-
Super-resolution Full Polarimetric Imaging for Radio Interferometry with Sparse Modeling
Authors:
Kazunori Akiyama,
Shiro Ikeda,
Mollie Pleau,
Vincent L. Fish,
Fumie Tazaki,
Kazuki Kuramochi,
Avery Broderick,
Jason Dexter,
Monika Mościbrodzka,
Michael Gowanlock,
Mareki Honma,
Sheperd S. Doeleman
Abstract:
We propose a new technique for radio interferometry to obtain super-resolution full polarization images in all four Stokes parameters using sparse modeling. The proposed technique reconstructs the image in each Stokes parameter from the corresponding full-complex Stokes visibilities by utilizing two regularization functions: the $\ell _1$-norm and total variation (TV) of the brightness distributio…
▽ More
We propose a new technique for radio interferometry to obtain super-resolution full polarization images in all four Stokes parameters using sparse modeling. The proposed technique reconstructs the image in each Stokes parameter from the corresponding full-complex Stokes visibilities by utilizing two regularization functions: the $\ell _1$-norm and total variation (TV) of the brightness distribution. As an application of this technique, we present simulated linear polarization observations of two physically motivated models of M87 with the Event Horizon Telescope (EHT). We confirm that $\ell _1$+TV regularization can achieve an optimal resolution of $\sim 25-30$\% of the diffraction limit $λ/D_{\rm max}$, which is the nominal spatial resolution of a radio interferometer for both the total intensity (i.e. Stokes $I$) and linear polarizations (i.e. Stokes $Q$ and $U$). This optimal resolution is better than that obtained from the widely used Cotton-Schwab CLEAN algorithm or from using $\ell _1$ or TV regularizations alone. Furthermore, we find that $\ell _1$+TV regularization can achieve much better image fidelity in linear polarization than other techniques over a wide range of spatial scales, not only in the super-resolution regime, but also on scales larger than the diffraction limit. Our results clearly demonstrate that sparse reconstruction is a useful choice for high-fidelity full-polarimetric interferometric imaging.
△ Less
Submitted 3 March, 2017; v1 submitted 1 February, 2017;
originally announced February 2017.
-
The Detection of TeV Gamma Rays from Crab using the Telescope Array Prototype
Authors:
S. Aiso,
M. Chikawa,
Y. Hayashi,
N. Hayashida,
K. Hibino,
H. Hirasawa,
K. Honda,
N. Hotta,
N. Inoue,
F. Ishikawa,
N. Ito,
S. Kabe,
F. Kajino,
T. Kashiwagi,
S. Kawakami,
Y. Kawasaki,
N. Kawasumi,
H. Kitamura,
K. Kuramochi,
K. Kurata,
E. Kusano,
H. Lafoux,
E. C. Loh,
Y. Matsubara,
T. Matsuyama
, et al. (29 additional authors not shown)
Abstract:
The Telescope Array prototype detectors were installed at Akeno Observatory and at the Utah Fly's Eye site. Using these detectors, we have observed the Crab Nebula and AGN's since the end of 1995. The successful detections of TeV gamma rays from Crab Nebula and Mkn501 are reported.
The Telescope Array prototype detectors were installed at Akeno Observatory and at the Utah Fly's Eye site. Using these detectors, we have observed the Crab Nebula and AGN's since the end of 1995. The successful detections of TeV gamma rays from Crab Nebula and Mkn501 are reported.
△ Less
Submitted 8 June, 1997;
originally announced June 1997.
