nifti 0.17.0

Rust implementation of the NIfTI file format
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
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121

#[cfg(feature = "nalgebra_affine")]
extern crate nalgebra;
#[cfg(feature = "nalgebra_affine")]
extern crate nifti;

#[cfg(feature = "nalgebra_affine")]
mod nalgebra_affine {
    use nalgebra::Vector4;
    use nifti::{affine::Affine4, NiftiHeader};

    #[test]
    fn affine() {
        let mut header = NiftiHeader::default();
        let affine = Affine4::from_diagonal(&Vector4::new(2.0, 2.0, 2.0, 1.0));
        header.set_affine(&affine);
        assert_eq!(affine, header.affine());
        assert_eq!(header.sform_code, 2);
        assert_eq!(header.qform_code, 0);
    }

    #[test]
    #[rustfmt::skip]
    fn sform() {
        let header = NiftiHeader {
            sform_code: 1,
            qform_code: 0,
            srow_x: [2.4, -0.0008, -0.0411765, -114.766396],
            srow_y: [0.1, 2.4995277, 0.0485984, -97.420204],
            srow_z: [0.4, -0.0485, 2.4991884, -89.12282],
            .. NiftiHeader::default()
        };
        let real_affine = Affine4::new(
            2.4, -0.0008,    -0.0411765, -114.766396,
            0.1,  2.4995277,  0.0485984,  -97.420204,
            0.4, -0.0485,     2.4991884,  -89.12282,
            0.0,  0.0,        0.0,          1.0
        );
        assert_eq!(header.affine(), real_affine);
    }

    #[test]
    #[rustfmt::skip]
    fn qform() {
        let header = NiftiHeader {
            sform_code: 0,
            qform_code: 1,
            pixdim: [-1.0, 0.9375, 0.9375, 3.0, 0.0, 0.0, 0.0, 0.0],
            quatern_b: 0.0,
            quatern_c: 1.0,
            quatern_d: 0.0,
            quatern_x: 59.557503,
            quatern_y: 73.172,
            quatern_z: 43.4291,
            .. NiftiHeader::default()
        };
        let real_affine = Affine4::new(
            -0.9375, 0.0,    0.0, 59.557503,
            0.0,     0.9375, 0.0, 73.172,
            0.0,     0.0,    3.0, 43.4291,
            0.0,     0.0,    0.0, 1.0
        );
        assert_eq!(header.affine(), real_affine);
    }

    #[test]
    #[rustfmt::skip]
    fn both_valid() {
        let header = NiftiHeader {
            sform_code: 1,
            srow_x: [2.4, 0.0, 0.0, -114.766396],
            srow_y: [0.1, 2.4, 0.0, -97.420204],
            srow_z: [0.4, 0.4, 2.4, -89.12282],
            // All this should be ignored
            qform_code: 1,
            pixdim: [-1.0, 0.9375, 0.9375, 3.0, 0.0, 0.0, 0.0, 0.0],
            quatern_b: 0.0,
            quatern_c: 1.0,
            quatern_d: 0.0,
            quatern_x: 59.0,
            quatern_y: 73.0,
            quatern_z: 43.0,
            .. NiftiHeader::default()
        };
        let real_affine = Affine4::new(
            2.4, 0.0, 0.0, -114.766396,
            0.1, 2.4, 0.0, -97.420204,
            0.4, 0.4, 2.4, -89.12282,
            0.0, 0.0, 0.0, 1.0
        );
        assert_eq!(header.affine(), real_affine);
    }

    #[test]
    #[rustfmt::skip]
    fn none_valid() {
        let header = NiftiHeader {
            dim: [3, 100, 100, 100, 0, 0, 0, 0],
            pixdim: [-1.0, 0.9, 0.9, 3.0, 0.0, 0.0, 0.0, 0.0],
            // All this should be ignored, because only `shape_zoom_affine` will be called.
            sform_code: 0,
            srow_x: [1.0, 0.0, 0.0, 1.0],
            srow_y: [0.0, 1.0, 0.0, 1.0],
            srow_z: [0.0, 0.0, 1.0, 1.0],
            qform_code: 0,
            quatern_b: 0.0,
            quatern_c: 1.0,
            quatern_d: 0.0,
            quatern_x: 59.0,
            quatern_y: 73.0,
            quatern_z: 43.0,
            .. NiftiHeader::default()
        };
        let real_affine = Affine4::new(
            -0.9, 0.0, 0.0,   44.55,
            0.0,  0.9, 0.0,  -44.55,
            0.0,  0.0, 3.0, -148.5,
            0.0, 0.0, 0.0, 1.0
        );
        assert_eq!(header.affine(), real_affine);
    }
}