openkspace-recon 0.2.0

K-space -> image reconstruction: IFFT, coil combination, corrections.
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59

//! OpenKSpace-Recon: K-space -> image reconstruction math.
//!
//! Core pipeline:
//!   1. (optional) noise pre-whitening      -- `prewhiten`
//!   2. (optional) navigator phase correction -- `phasecorr`
//!   3. fftshift -> IFFT -> fftshift           -- `fft` + `shift`
//!   4. RSS coil combine                     -- `coil`
//!   5. (optional) centre crop to recon matrix -- `crop`
//!
//! Reconstructions are composed behind the [`ReconStrategy`] trait so that
//! future parallel-imaging / compressed-sensing back-ends can slot in
//! without touching this crate's public shape.
//!
//! # Quick start
//!
//! ```rust,no_run
//! use openkspace_recon::{IfftRss, ReconStrategy};
//! use openkspace_io::ismrmrd::IsmrmrdFile;
//!
//! let file = IsmrmrdFile::open("scan.h5").unwrap();
//! let mut strategy = IfftRss::default();
//! strategy.prewhiten = true;
//! strategy.phase_correct = true;
//! let volume = strategy.reconstruct(&file).unwrap();
//! println!("Reconstructed {} slices", volume.data.shape()[0]);
//! ```

pub mod coil;
pub mod crop;
pub mod cs;
pub mod espirit;
pub mod fft;
pub mod grappa;
pub mod oversampling;
pub mod partial_fourier;
pub mod phasecorr;
pub mod prewhiten;
pub mod sense;
pub mod sensitivity;
pub mod shift;
pub mod strategy;
pub mod wavelet;

pub use coil::{rss_combine, rss_combine_4d};
pub use crop::center_crop_3d;
pub use cs::{fista_cs_single_coil, CsError};
pub use espirit::espirit_sensitivity_maps;
pub use fft::{ifft1_inplace, ifft2_inplace, ifft3_inplace};
pub use grappa::{GrappaError, GrappaKernel, SamplingPattern};
pub use oversampling::OversamplingRemover;
pub use partial_fourier::{homodyne_reconstruct, PartialFourierPlan};
pub use phasecorr::PhaseCorrector;
pub use prewhiten::NoisePrewhitener;
pub use sense::{sense_gfactor_1d, sense_unfold_1d, SenseError};
pub use sensitivity::walsh_sensitivity_maps;
pub use shift::{fftshift_axis, ifftshift_axis};
pub use strategy::{
    CsRss, FftMode, GrappaRss, IfftRss, ImageVolume, ReconStrategy, SenseMapSource, SenseRss,
};