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
use nalgebra::Point3;
use crate::ArraySequence;
pub type Point = Point3<f32>;
pub type Points = Vec<Point>;
pub type Streamlines = ArraySequence<Point>;
pub type TractogramItem = (Points, ArraySequence<f32>, Vec<f32>);
pub type RefTractogramItem<'data> = (&'data [Point], &'data [f32], &'data [f32]);
#[derive(Clone, PartialEq)]
pub struct Tractogram {
pub streamlines: Streamlines,
pub scalars: ArraySequence<f32>,
pub properties: ArraySequence<f32>,
}
impl Tractogram {
pub fn new(
streamlines: Streamlines,
scalars: ArraySequence<f32>,
properties: ArraySequence<f32>,
) -> Tractogram {
Tractogram { streamlines, scalars, properties }
}
pub fn item(&self, idx: usize) -> RefTractogramItem {
let scalars = if self.scalars.is_empty() { &[] } else { &self.scalars[idx] };
let properties = if self.properties.is_empty() { &[] } else { &self.properties[idx] };
(&self.streamlines[idx], scalars, properties)
}
}
impl<'data> IntoIterator for &'data Tractogram {
type Item = RefTractogramItem<'data>;
type IntoIter = TractogramIterator<'data>;
fn into_iter(self) -> Self::IntoIter {
TractogramIterator { tractogram: self, index: 0..self.streamlines.len() }
}
}
pub struct TractogramIterator<'data> {
tractogram: &'data Tractogram,
index: std::ops::Range<usize>,
}
impl<'data> Iterator for TractogramIterator<'data> {
type Item = RefTractogramItem<'data>;
fn next(&mut self) -> Option<Self::Item> {
let idx = self.index.next()?;
Some(self.tractogram.item(idx))
}
fn size_hint(&self) -> (usize, Option<usize>) {
(0, Some(self.tractogram.streamlines.len()))
}
}
impl<'data> ExactSizeIterator for TractogramIterator<'data> {}
impl<'data> DoubleEndedIterator for TractogramIterator<'data> {
fn next_back(&mut self) -> Option<Self::Item> {
let idx = self.index.next_back()?;
Some(self.tractogram.item(idx))
}
}