2019-04-11 01:30:40 8 Comments

The first ever black hole was "pictured" recently, per an announcement made on 10th April, 2019. See for example: https://www.bbc.com/news/science-environment-47873592 .

It has been claimed that state-of-the-art imaging algorithms were an enabler for this historic success. Does anybody care to describe the difficulties, and (quite certainly non-trivial) mathematics that went into this effort ?

### Related Questions

#### Sponsored Content

#### 2 Answered Questions

### [SOLVED] When is a mapping the proximity operator of some convex function?

**2016-05-17 21:21:14****dohmatob****595**View**9**Score**2**Answer- Tags: reference-request fa.functional-analysis convex-analysis signal-analysis image-processing

#### 0 Answered Questions

### Strategies to make the sum operation bijective in a contribution problem

**2018-10-21 12:03:04****Wilem2****23**View**0**Score**0**Answer- Tags: inverse-problems

#### 1 Answered Questions

### [SOLVED] Is there any paper which summarizes the mathematical foundation of deep learning?

**2018-10-11 08:28:05****almnagako****1339**View**15**Score**1**Answer- Tags: applied-mathematics machine-learning

#### 0 Answered Questions

### Can the wavelet bispectrum be normalised so that its integral "gives the right answer"?

**2018-02-23 03:08:57****Julian Newman****197**View**4**Score**0**Answer- Tags: signal-analysis wavelets

#### 2 Answered Questions

### [SOLVED] Choosing the order of Tikhonov regularization of an inverse problem

**2013-10-27 23:14:49****shva****2420**View**2**Score**2**Answer- Tags: regularization inverse-problems

#### 4 Answered Questions

### [SOLVED] Mathematics TV clips

**2011-02-23 19:54:36****Daniel Moskovich****1380**View**23**Score**4**Answer- Tags: soft-question popularization

#### 0 Answered Questions

### The jump set of $SBV$ function for different value of parameter in image denoising problem

**2016-06-27 16:47:38****JumpJump****55**View**1**Score**0**Answer- Tags: ap.analysis-of-pdes real-analysis na.numerical-analysis inverse-problems image-processing

#### 0 Answered Questions

### Sample based inversion of the Radon transform

**2015-02-17 09:38:24****user45183****73**View**4**Score**0**Answer- Tags: pr.probability inverse-problems

#### 3 Answered Questions

### [SOLVED] Estimating the derivative of a noisy, non-uniformly sampled function

**2011-02-13 17:10:33****Grönwall****1841**View**1**Score**3**Answer- Tags: signal-analysis

#### 2 Answered Questions

### [SOLVED] Switching from pure mathematics (e.g. geometry) to more applied areas (e.g imaging) after Ph.D., as postdoc and chance of getting such a postdoc?

**2013-06-28 16:35:06****Let's talk math****3021**View**14**Score**2**Answer- Tags: dg.differential-geometry computer-science teichmuller-theory mathjobs image-processing

## 1 comments

## @Carlo Beenakker 2019-04-11 07:24:58

Essential elements$^*$ of the reconstruction algorithm were developed at MIT under the name CHIRP =

Continuous High-resolution Image Reconstruction using Patch priors, as described in Computational Imaging for VLBI Image Reconstruction (2015).The difficulty of VLBI (

Very Long Baseline Interferometry Image) reconstruction is that the inversion problem is highly-ill posed, there are many images that explain the data. The challenge is to find an explanation that respects our prior assumptions about the “visual” universe while still satisfying the observed data. Bayesian approaches are generally employed for that purpose, in CHIRP machine learning is used to automatically identify visual patterns --- obviating the need for hand training of the algorithm.A key technical innovation is a way to correct for the delays in the signal received from the various telescopes. The delays are difficult to predict, since they depend the local variations in the speed of the radio waves through the noisy atmosphere. CHIRP adopts an algebraic solution known as phase closure to this problem: If the measurements from three telescopes are multiplied, the extra delays caused by atmospheric noise cancel each other out.

One test case that shows the resolving power of CHIRP, compared to a competing algorithm (BU) is shown below (taken from the MIT paper). Notice how CHIRP is able to resolve 2 separate, previously unresolved, bright emissions in the blazar OJ287.

$^*$ UPDATE:This statement must be qualified, as stated by the lead author, Katie Bouman:“No one person or algorithm made the image. It was actually made by combining the images produced by 3 separate imaging pipelines, leading to a more powerful result than we could ever achieve with any single method.”