create_chauffert_density

create_chauffert_density#

mrinufft.trajectories.sampling.create_chauffert_density(shape, wavelet_basis, nb_wavelet_scales, *, verbose=False)[source]#

Create a density based on Chauffert’s method.

This is a reproduction of the proposition from [CCW13]. A sampling density is derived from compressed sensing equations to maximize guarantees of exact image recovery for a specified sparse wavelet domain decomposition.

Parameters:

  • shape (tuple[int, ...]) – The shape of the density grid.

  • wavelet_basis (str, pywt.Wavelet) – The wavelet basis to use for wavelet decomposition, either as a built-in wavelet name from the PyWavelets package or as a custom pywt.Wavelet object.

  • nb_wavelet_scales (int) – The number of wavelet scales to use in decomposition.

  • verbose (bool, optional) – If True, displays a progress bar. Default to False.

Returns:

A density array created based on wavelet transform coefficients.

Return type:

NDArray

See also

pywt.wavelist

A list of wavelet decompositions available in the PyWavelets package used inside the function.

pywt.Wavelet

A wavelet object accepted to generate Chauffert densities.

References

[CCW13]

Chauffert, Nicolas, Philippe Ciuciu, and Pierre Weiss. “Variable density compressed sensing in MRI. Theoretical vs heuristic sampling strategies.” In 2013 IEEE 10th International Symposium on Biomedical Imaging, pp. 298-301. IEEE, 2013.