wav-dwt 0.2.1

An experimental implementation of the A Trous Wavelet Decomposition algorithms for audio instead of images
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

use crate::kernels::{B3SplineKernel, Kernel, LinearInterpolationKernel};
use convolve_nd::exports::ndarray::Array1;
use convolve_nd::Convolution;

use crate::layer::WaveletLayerBuffer;

pub trait WaveletDecompose {
    fn wavelet_decompose(&mut self, kernel: Kernel, scale: usize) -> WaveletLayerBuffer;
}

impl WaveletDecompose for Array1<f32> {
    fn wavelet_decompose(&mut self, kernel: Kernel, scale: usize) -> WaveletLayerBuffer {
        let stride = 2_usize.pow(scale as u32);
        let mut current_data: Array1<f32> = self.clone();

        match kernel {
            Kernel::LinearInterpolationKernel => {
                current_data.convolve(LinearInterpolationKernel::new().into(), stride);
            }
            Kernel::LowScaleKernel => {
                unimplemented!("Low scale is not a separable kernel");
            }
            Kernel::B3SplineKernel => current_data.convolve(B3SplineKernel::new().into(), stride),
        }

        let final_data = self.clone() - &current_data;
        *self = current_data;

        WaveletLayerBuffer { data: final_data }
    }
}