3
$\begingroup$

A common approach to image reconstruction is to solve the convex optimization problem \begin{equation} \text{minimize} \quad \frac12 \| Ax - b \|^2 + \gamma \| Dx \|_1 \end{equation} where $b$ is a vector of measurements (perhaps a blurry, noisy image), $A$ is a linear transformation that maps the reconstructed image to the vector of measurements (often $A$ is a convolution operator), and $D$ is a wavelet or tight frame analysis operator. The optimization variable is $x$, an image stored as a rectangular array of numbers.

What are the most popular choices for $D$ (and what are the most popular and easiest to use software implementations of $D$)? I've experimented with both curvelet and shearlet tight frames.

Is there a good choice of $D$ that is included in Matlab's wavelet toolbox?

$\endgroup$
2
  • $\begingroup$ What did you end up doing? $\endgroup$
    – Royi
    Commented Jul 14, 2023 at 8:45
  • $\begingroup$ @Royi Although I’ve experimented with curvelets and shearlets and common types of wavelets (like Daubechies, I believe), I still don’t have a clear answer in my mind as to which of these options is best or most popular, or if some other option that I’m not aware of has been found to be superior. I’m most interested in MR image reconstruction, if anyone else would like to chime in. $\endgroup$
    – littleO
    Commented Jul 14, 2023 at 9:08

1 Answer 1

2
$\begingroup$

EDIT: some code made available] From a personal experience, and in the work of colleagues, we observed that multi-band dual-tree wavelets, a type of low-redundancy oriented decompositions based on Hilbert pair of wavelets, were an interesting alternative to higher redundancy transformations (curvelets, shearlets). You can have a two-fold or four fold redundancy. Even if you can have fewer directions than with other transforms (a survey on 2D geometric wavelet-like representations is given in A panorama on multiscale geometric representations, 2011, L. Jacques et al.), the Hilbert properties yields useful statistical properties with Gaussian noise.

Dual tree of an 8-band dual-tree wavelet

Caveat: such design sometimes requires to relax the tight frame condition (depending on the optimization used), see Relaxing Tight Frame Condition in Parallel Proximal Methods for Signal Restoration, IEEE TSP, 2012. Another type of image restoration with the same wavelets and no issue with the tight frame condition is discussed in Sparse adaptive template matching and filtering for 2D seismic images with dual-tree wavelets and proximal methods, 2015, ICIP.

Those dual-tree wavelets originate from Image Analysis Using a Dual-Tree $M$-Band Wavelet Transform, 2006, and the improvement over curvelets was shown, for instance (for standard denoising), in A Nonlinear Stein Based Estimator for Multichannel Image Denoising, 2008.

[EDIT] The 2D code is now available, embedded in a denoising toolbox: $M$-band 2D dual-tree (Hilbert) wavelet multicomponent image denoising, so you might need a little work. I hope I can work on a cleaner implementation. Recently related discussions in How to implement a $j$-level $M$-band wavelet transform of an image? and Code for a wavelet based hilbert transform?

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.