Struct Ist 

pub struct Ist<'s> { /* private fields */ }

Iterative Soft Thresholding; an algorithm for reconstructing NUS data.

The Reconstructor::reconstruct implementation returns the fourier transform after interpolating unsampled data.

Interpolating is done by setting the unsampled positions in the data to zero, then iteratively thresholding the spectrum and setting the unsampled points to the IFT of what is above the threshold. On each iteration, the threshold is lowered until it reaches zero.

WARNING: IST is difficult to get right and it is unlikely that this does. These results should be trusted less than those of other, better tested, implementations. Feel free to contribute if you know any secret sauce.

A. Shchukina, P. Kasprzak, R. Dass et al. J Biomol NMR 68, 79–98 (2017). https://doi.org/10.1007/s10858-016-0068-3 (This implements IST-S with Virtual Echo)

Implementations

impl<'s> Ist<'s>

pub fn new( schedule: &'s Schedule<Ix1>, threshold_function: impl Fn(f64) -> f64 + 's, iterations: usize, mode: DisplayMode, ) -> Ist<'s>

Create a new Ist instance

schedule is the sampling schedule used. threshold_function is a function that returns the relative threshold for a given iteration. The function is called with values between zero and one, where each iteration i ∈ [0, iterations) is mapped to the range [0, 1]. The value returned must be between zero and one. iterations is the number of iterations to perform. mode is whether the data will be displayed in absolute value or real mode.

Trait Implementations

impl Debug for Ist<'_>

fn fmt(&self, f: &mut Formatter<'_>) -> Result

Formats the value using the given formatter.

impl Reconstructor for Ist<'_>

fn reconstruct(&mut self, data: &mut [Complex<f64>])

Reconstruct a complex signal in place.

fn reconstruct_from_reals(&mut self, data: &[f64], out: &mut [Complex<f64>])

Reconstruct data represented as a list of floating point values. Converts to complex numbers by setting the imaginary component to zero. Your data will have mirrored quadrature peaks.

fn quadrature_reconstruct( &mut self, cos: &[f64], sin: &[f64], out: &mut [Complex<f64>], )

Reconstruct a spectrum from cos and sin signals using quadrature.