bvreader 0.1.5

Reader for the BrainVision data 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


use crate::bv_reader::header::parser::parse_chan_info::ChannelInfo;

//use std::fmt::Error;

use crate::bv_reader::bv_error::Error;

pub fn scale_channels(data: &mut Vec<Vec<f32>>, info: &Vec<ChannelInfo>) -> Result<(), Error> {
    if data.len() != info.len() {return Err(Error::ChannelDataMismatch(data.len(), info.len()))}; // not the same amount of channels in data and info
    let chan_len = data[0].len();

    for i in 0..data.len() {
        if data[i].len() != chan_len {return Err(Error::ChannelDataMalformed(chan_len, data[i].len()))} // data malformed, unequal sample length per channel
        // iterate through channels
        let resolution = info[i].resolution;
        if resolution.is_none() {continue} // dont need to scale if there is no scale available, should probably error or log this?
        let resolution = f32::from(resolution.unwrap());
        for j in 0..chan_len {
            data[i][j] *= resolution; 
        }
    }
    Ok(())
}


#[cfg(test)]
mod tests {
    use super::*;

    #[test]
    fn test_scale_channels() {
        let mut in_data: Vec<Vec<f32>> = vec![vec![
            110., 20., 110., 20.,
        ],vec![
            80., 30., 80., 30.,
        ]];
        let mut in_info = ChannelInfo::default();
        in_info.resolution = Some(0.5);


        scale_channels(&mut in_data, &vec![in_info.clone(), in_info]).unwrap();

        let expected: Vec<Vec<f32>> = vec![vec![
            55., 10., 55., 10.,
        ],vec![
            40., 15., 40., 15.,
        ]];

        assert_eq!(in_data, expected)

    }
}