Struct AxialGrid 

pub struct AxialGrid {
    pub dims: (usize, usize, usize),
    pub voxel_size: (f64, f64, f64),
    pub affine: [f64; 16],
    src_dims: (usize, usize, usize),
    world_min: [f64; 3],
    inv_r: [[f64; 3]; 3],
    t: [f64; 3],
}

The cardinal-aligned grid an oblique volume resamples onto, plus the mapping back to the original voxel space. Build once with axial_grid_for and reuse for every volume that shares the geometry (magnitude, phase, mask), so they land on identical grids.

Fields

dims: (usize, usize, usize)

voxel_size: (f64, f64, f64)

affine: [f64; 16]

Diagonal voxel→world affine of the new grid (row-major 4×4).

src_dims: (usize, usize, usize)

Source dimensions this grid was built for.

world_min: [f64; 3]

inv_r: [[f64; 3]; 3]

Inverse of the source 3×3, for world→source-voxel.

t: [f64; 3]

Implementations

impl AxialGrid

fn source_coord( &self, ni: usize, nj: usize, nk: usize, ) -> Option<(f64, f64, f64)>

Source-voxel coordinate of a target voxel, or None if it falls outside the source.

pub fn sample(&self, data: &[f64]) -> Vec<f64>

Trilinearly resample one continuous volume onto this grid. Voxels outside the source are 0.

pub fn sample_nearest(&self, data: &[u8]) -> Vec<u8>

Nearest-neighbour resample, for labels and binary masks.