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
use std::sync::Arc;
use crate::file_reader::AsyncFileReader;
use crate::model::{EDFHeader, EDF_HEADER_BYTE_SIZE};
use futures::future::{err};
use futures::Future;
use std::io::Error;
pub struct AsyncEDFReader<T: AsyncFileReader> {
pub edf_header: Arc<EDFHeader>,
file_reader: T,
}
impl<T: 'static + AsyncFileReader> AsyncEDFReader<T> {
pub fn init_with_file_reader(
file_reader: T,
) -> Box<Future<Item = AsyncEDFReader<T>, Error = std::io::Error>> {
Box::new(
file_reader
.read_async(0, 256)
.map(|general_header_raw: Vec<u8>| {
let mut edf_header = EDFHeader::build_general_header(general_header_raw);
file_reader
.read_async(
256,
edf_header.number_of_signals * EDF_HEADER_BYTE_SIZE as u64,
)
.map(|channel_headers_raw| {
edf_header.build_channel_headers(channel_headers_raw);
AsyncEDFReader {
edf_header : Arc::new(edf_header),
file_reader,
}
})
})
.flatten(),
)
}
pub fn read_data_window(
&self,
start_time_ms: u64,
duration_ms: u64,
) -> Box<Future<Item = Vec<Vec<f32>>, Error = std::io::Error>> {
if let Err(e) = super::check_bounds(start_time_ms, duration_ms, &self.edf_header) {
return Box::new(err::<Vec<Vec<f32>>, Error>(e));
}
let first_block_start_time = start_time_ms - start_time_ms % self.edf_header.block_duration;
let first_block_index = first_block_start_time / self.edf_header.block_duration;
let number_of_blocks_to_get =
(duration_ms as f64 / self.edf_header.block_duration as f64).ceil() as u64;
let offset = self.edf_header.byte_size_header
+ first_block_index * self.edf_header.get_size_of_data_block();
let header = self.edf_header.clone();
Box::new(
self.file_reader
.read_async(
offset,
number_of_blocks_to_get * self.edf_header.get_size_of_data_block(),
)
.map(move |data: Vec<u8>| {
let mut result: Vec<Vec<f32>> = Vec::new();
for _ in 0..header.number_of_signals {
result.push(Vec::new());
}
let mut index = 0;
for _ in 0..number_of_blocks_to_get {
for (j, channel) in header.channels.iter().enumerate() {
for _ in 0..channel.number_of_samples_in_data_record {
let sample = super::get_sample(&data, index) as f32;
result[j].push(
(sample - channel.digital_minimum as f32)
* channel.scale_factor
+ channel.physical_minimum,
);
index += 1;
}
}
}
result
})
)
}
}