Blind Deconvolution for Ultrasound Sequences Using a Noninverse Greedy Algorithm

Algorithm 1

Noninverse greedy deconvolution.

Input: signal envelope 𝑦(𝑛), PSF ℎ(𝑛), threshold 𝑘.

Output: reflectivity function 𝑥(𝑛).

Initialisation:

 𝑥(𝑛)←0; 𝑅(𝑛)←𝑦(𝑛); 𝑖←0;

Repeat

 (i) 𝑛_𝑖=argmax_𝑛(𝑅_𝑖(𝑛))

 (ii) 𝑥(𝑛_𝑖)←𝑅_𝑖(𝑛_𝑖);

 (iii) 𝑅_𝑖+1←𝑅_𝑖−ℎ⊗𝑅_𝑖(𝑛_𝑖);

 (iv) 𝑖←𝑖+1;

Until stop criterion (max(𝑅_𝑖(𝑛))<𝑘)