Superresolution Interferometric Imaging with Sparse Modeling Using Total Squared Variation --- Application to Imaging the Black Hole Shadow

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-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*.

arXiv comment18 pages, 7 figures, ApJ in Press

References in corpus (55)

• Astropy: A Community Python Package for Astronomy
• Coevolution (Or Not) of Supermassive Black Holes and Host Galaxies
• Measuring Distance and Properties of the Milky Way's Central Supermassive Black Hole with Stellar Orbits
• Event-horizon-scale structure in the supermassive black hole candidate at the Galactic Centre
• The Black-Hole Mass in M87 from Gemini/NIFS Adaptive Optics Observations
• Jet Launching Structure Resolved Near the Supermassive Black Hole in M87
• The M87 Black Hole Mass from Gas-dynamical Models of Space Telescope Imaging Spectrograph Observations
• RadioAstron -- a Telescope with a Size of 300 000 km: Main Parameters and First Observational Results
• Compressed sensing imaging techniques for radio interferometry
• High Resolution Linear Polarimetric Imaging for the Event Horizon Telescope
• Imaging Optically-Thin Hot Spots Near the Black Hole Horizon of Sgr A* at Radio and Near-Infrared Wavelengths
• Resolved Magnetic-Field Structure and Variability Near the Event Horizon of Sagittarius A*
• 1.3 mm Wavelength VLBI of Sagittarius A*: Detection of Time-Variable Emission on Event Horizon Scales
• The old nuclear star cluster in the Milky Way: dynamics, mass, statistical parallax, and black hole mass
• Evidence for Low Black Hole Spin and Physically Motivated Accretion Models from Millimeter VLBI Observations of Sagittarius A*
• Imaging the Schwarzschild-radius-scale Structure of M87 with the Event Horizon Telescope using Sparse Modeling
• The Intrinsic Size of Sagittarius A* from 0.35 cm to 6 cm
• Methods for detecting flaring structures in Sagittarius A* with high frequency VLBI
• The size of the jet launching region in M87
• Sparsity Averaging Reweighted Analysis (SARA): a novel algorithm for radio-interferometric imaging
• Observational appearance of inefficient accretion flows and jets in 3D GRMHD simulations: Application to Sgr~A*
• The Power of Imaging: Constraining the Plasma Properties of GRMHD Simulations using EHT Observations of Sgr A*
• Plasma diagnostic potential of 2p4f in N${}^{+}$ -- accurate wavelengths and oscillator strengths
• A General Relativistic Null Hypothesis Test with Event Horizon Telescope Observations of the black-hole shadow in Sgr A*
• PURIFY: a new approach to radio-interferometric imaging
• 230 GHz VLBI observations of M87: event-horizon-scale structure at the enhanced very-high-energy $\gamma$-ray state in 2012
• The application of compressive sampling to radio astronomy I: Deconvolution
• Imaging the Supermassive Black Hole Shadow and Jet Base of M87 with the Event Horizon Telescope
• LOFAR Sparse Image Reconstruction
• Super-resolution Full Polarimetric Imaging for Radio Interferometry with Sparse Modeling
• Faraday rotation in GRMHD simulations of the jet launching zone of M87
• Modeling Seven Years of Event Horizon Telescope Observations with Radiatively Inefficient Accretion Flow Models
• Scalable splitting algorithms for big-data interferometric imaging in the SKA era
• Theory and Simulations of Refractive Substructure in Resolved Scatter-Broadened Images
• Persistent Asymmetric Structure of Sagittarius A* on Event Horizon Scales
• Imaging an Event Horizon: Mitigation of Scattering Toward Sagittarius A*
• Fast Variability and mm/IR flares in GRMHD Models of Sgr A* from Strong-Field Gravitational Lensing
• Dynamical Imaging with Interferometry
• Interstellar Scintillation and the Radio Counterpart of the Fast Radio Burst FRB150418
• MORESANE: MOdel REconstruction by Synthesis-ANalysis Estimators. A sparse deconvolution algorithm for radio interferometric imaging
• Compressed sensing for wide-field radio interferometric imaging
• Spread spectrum for imaging techniques in radio interferometry
• Constraining the Structure of Sagittarius A*'s Accretion Flow with Millimeter-VLBI Closure Phases
• Imaging an Event Horizon: Mitigation of Source Variability of Sagittarius A*
• The Effects of Accretion Flow Dynamics on the Black Hole Shadow of Sagittarius A${}^{\ast }$
• Stochastic Optics: A Scattering Mitigation Framework for Radio Interferometric Imaging
• Extreme Brightness Temperatures and Refractive Substructure in 3C273 with RadioAstron
• Fine-scale structure of the quasar 3C 279 Measured with 1.3 mm very long baseline interferometry
• The Optics of Refractive Substructure
• An accelerated splitting algorithm for radio-interferometric imaging: when natural and uniform weighting meet
• Compressed sensing reconstruction of a string signal from interferometric observations of the cosmic microwave background
• PRECL: A new method for interferometry imaging from closure phase
• Resolving the inner jet structure of 1924-292 with the EVENT HORIZON TELESCOPE
• Doppler Tomography by Total Variation Minimization
• Accelerating cross-validation with total variation and its application to super-resolution imaging

Cited by in corpus (34)

• First M87 Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole
• First M87 Event Horizon Telescope Results. IV. Imaging the Central Supermassive Black Hole
• First M87 Event Horizon Telescope Results. II. Array and Instrumentation
• Discriminating Accretion States via Rotational Symmetry in Simulated Polarimetric Images of M87
• Black Hole Shadow as a Test of General Relativity: Quadratic Gravity
• Forming intracluster gas in a galaxy protocluster at a redshift of 2.16
• Hybrid Very Long Baseline Interferometry Imaging and Modeling with Themis
• Imaging Black Holes and Jets with a VLBI Array Including Multiple Space-Based Telescopes
• Metrics and Motivations for Earth-Space VLBI: Time-Resolving Sgr A* with the Event Horizon Telescope
• Multi-frequency Black Hole Imaging for the Next-Generation Event Horizon Telescope
• Collimation of the relativistic jet in the quasar 3C 273
• Selective Dynamical Imaging of Interferometric Data
• Black Hole Spin Signature in the Black Hole Shadow of M87 in the Flaring State
• Unravelling the Innermost Jet Structure of OJ 287 with the First GMVA+ALMA Observations
• Testing General Relativity with the Black Hole Shadow Size and Asymmetry of Sagittarius A*: Limitations from Interstellar Scattering
• Variational Image Feature Extraction for the Event Horizon Telescope
• Global Mapping of an Exo-Earth using Sparse Modeling
• Polarization images of accretion flow around supermassive black holes: imprints of toroidal field structure
• ALMA Super-resolution Imaging of T Tau: r = 12 au Gap in the Compact Dust Disk around T Tau N
• The Evolving Radio Photospheres of Long-Period Variable Stars
• alpha-Deep Probabilistic Inference (alpha-DPI): efficient uncertainty quantification from exoplanet astrometry to black hole feature extraction
• Super-resolution Imaging of the Protoplanetary Disk HD 142527 Using Sparse Modeling
• VLBI imaging of black holes via second moment regularization
• Measuring the Ellipticity of M 87* Images
• Regularized Maximum Likelihood Image Synthesis and Validation for ALMA Continuum Observations of Protoplanetary Disks
• Search for Alignment of Disk Orientations in Nearby Star-Forming Regions: Lupus, Taurus, Upper Scorpius, $\rho$ Ophiuchi, and Orion
• CLEAN imaging systematics of M87 radio jet
• Event-horizon-scale Imaging of M87* under Different Assumptions via Deep Generative Image Priors
• 3D Detection and Characterisation of ALMA Sources through Deep Learning
• Imaging the event horizon of M87* from space on different timescales
• Global Mapping of Surface Composition on an Exo-Earth Using Sparse Modeling
• Resolving the inner jet of PKS 1749+096 with super-resolution VLBA images at 7 mm
• Deblurring galaxy images with Tikhonov regularization on magnitude domain
• Very-Long Baseline Interferometry Imaging with Closure Invariants using Conditional Image Diffusion