Wrap a denoiser with a Rician VST (variance-stabilizing transform)

Converts magnitude data to approximately homoscedastic Gaussian via VST, applies `denoise_fun`, then inverts the VST.

Usage

vst_wrap(x, sigma = NULL, denoise_fun, ...)

Arguments

x

3D/4D magnitude array or neuroim2 volume/vec.

sigma

Optional noise sigma. If `NULL` and a 4D input is provided, sigma is estimated with `estimate_sigma_rician`.

denoise_fun

Function that takes an array and returns an array of the same dimensions.

...

Additional arguments forwarded to `denoise_fun`.

Value

The denoised data, same shape/type as `x` (wrapped via `nv_wrap_like`).