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

use crate::kernels::{B3SplineKernel, Kernel, LinearInterpolationKernel};
use convolve_image::Convolution;
use ndarray::{Array2, Array3};

use crate::layer::WaveletLayerBuffer;

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

impl WaveletDecompose for Array2<f32> {
    fn wavelet_decompose(&mut self, kernel: Kernel, pixel_scale: usize) -> WaveletLayerBuffer {
        let stride = 2_usize.pow(pixel_scale as u32);
        let mut current_data = 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::Grayscale { data: final_data }
    }
}

impl WaveletDecompose for Array3<f32> {
    fn wavelet_decompose(&mut self, kernel: Kernel, pixel_scale: usize) -> WaveletLayerBuffer {
        let stride = 2_usize.pow(pixel_scale as u32);
        let mut current_data = 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::Rgb { data: final_data }
    }
}