adaptivegrain/src/mask.rs

use super::PLUGIN_NAME;
use failure::Error;
use num::{NumCast, ToPrimitive};
use std::fmt::Debug;
use std::ops::Shl;
use std::ptr;
use vapoursynth::core::CoreRef;
use vapoursynth::format::ColorFamily;
use vapoursynth::plugins::{Filter, FrameContext};
use vapoursynth::prelude::*;
use vapoursynth::video_info::{Property, VideoInfo};

pub struct Mask<'core> {
    pub source: Node<'core>,
    pub luma_scaling: f32,
}

lazy_static! {
    static ref FLOAT_RANGE: Vec<f32> = {
        [0f32; 256]
            .iter()
            .enumerate()
            .map(|(i, _f)| (i as f32) * (1.0 / 256.0))
            .collect::<Vec<_>>()
    };
}

#[inline]
fn get_mask_value(x: f32, y: f32, luma_scaling: f32) -> f32 {
    f32::powf(
        1.0 - (x * (1.124 + x * (-9.466 + x * (36.624 + x * (-45.47 + x * 18.188))))),
        (y * y) * luma_scaling,
    )
}

#[inline]
fn from_property<T: Debug + Clone + Copy + Eq + PartialEq>(prop: Property<T>) -> T {
    match prop {
        Property::Constant(p) => p,
        Property::Variable => unreachable!(),
    }
}

fn filter_for_int<T>(
    frame: &mut FrameRefMut,
    src_frame: FrameRef,
    depth: u8,
    average: f32,
    luma_scaling: f32,
) where
    T: Component + NumCast + Clone + Shl<u8> + Debug,
    <T as Shl<u8>>::Output: ToPrimitive,
{
    let max = ((1 << depth) - 1) as f32;
    let lut: Vec<T> = FLOAT_RANGE
        .iter()
        .map(|x| T::from(get_mask_value(*x, average, luma_scaling) * max).unwrap())
        .collect();
    for row in 0..frame.height(0) {
        for (pixel, src_pixel) in frame
            .plane_row_mut::<T>(0, row)
            .iter_mut()
            .zip(src_frame.plane_row::<T>(0, row).iter())
        {
            let i = <usize as NumCast>::from(src_pixel.clone()).unwrap() >> (depth - 8);
            unsafe {
                ptr::write(pixel, lut[i].clone());
            }
        }
    }
}

fn filter_for_float(frame: &mut FrameRefMut, src_frame: FrameRef, average: f32, luma_scaling: f32) {
    let lut: Vec<f32> = FLOAT_RANGE
        .iter()
        .map(|x| get_mask_value(*x, average, luma_scaling))
        .collect();
    for row in 0..frame.height(0) {
        for (pixel, src_pixel) in frame
            .plane_row_mut::<f32>(0, row)
            .iter_mut()
            .zip(src_frame.plane_row::<f32>(0, row).iter())
        {
            unsafe {
                ptr::write(pixel, lut[(src_pixel * 255f32) as usize]);
            }
        }
    }
}

impl<'core> Filter<'core> for Mask<'core> {
    fn video_info(&self, _api: API, _core: CoreRef<'core>) -> Vec<VideoInfo<'core>> {
        let info = self.source.info();
        let format = match info.format {
            Property::Variable => unreachable!(),
            Property::Constant(format) => format,
        };
        vec![VideoInfo {
            format: Property::Constant(
                _core
                    .register_format(
                        ColorFamily::Gray,
                        format.sample_type(),
                        format.bits_per_sample(),
                        0,
                        0,
                    )
                    .unwrap(),
            ),
            flags: info.flags,
            framerate: info.framerate,
            num_frames: info.num_frames,
            resolution: info.resolution,
        }]
    }

    fn get_frame_initial(
        &self,
        _api: API,
        _core: CoreRef<'core>,
        context: FrameContext,
        n: usize,
    ) -> Result<Option<FrameRef<'core>>, Error> {
        self.source.request_frame_filter(context, n);
        Ok(None)
    }

    fn get_frame(
        &self,
        _api: API,
        core: CoreRef<'core>,
        context: FrameContext,
        n: usize,
    ) -> Result<FrameRef<'core>, Error> {
        let new_format = from_property(self.video_info(_api, core)[0].format);
        let mut frame = unsafe {
            FrameRefMut::new_uninitialized(
                core,
                None,
                new_format,
                from_property(self.source.info().resolution),
            )
        };
        let src_frame = self.source.get_frame_filter(context, n).ok_or_else(|| {
            format_err!("Could not retrieve source frame. This shouldn’t happen.")
        })?;
        let average = match src_frame.props().get::<f64>("PlaneStatsAverage") {
            Ok(average) => average as f32,
            Err(_) => bail!(format!(
                "{}: you need to run std.PlaneStats on the clip before calling this function.",
                PLUGIN_NAME
            )),
        };

        match from_property(self.source.info().format).sample_type() {
            SampleType::Integer => {
                let depth = from_property(self.source.info().format).bits_per_sample();
                match depth {
                    0..=8 => filter_for_int::<u8>(
                        &mut frame,
                        src_frame,
                        depth,
                        average,
                        self.luma_scaling,
                    ),
                    9..=16 => filter_for_int::<u16>(
                        &mut frame,
                        src_frame,
                        depth,
                        average,
                        self.luma_scaling,
                    ),
                    17..=32 => filter_for_int::<u32>(
                        &mut frame,
                        src_frame,
                        depth,
                        average,
                        self.luma_scaling,
                    ),
                    _ => bail!(format!(
                        "{}: input depth {} not supported",
                        PLUGIN_NAME, depth
                    )),
                }
            }
            SampleType::Float => {
                filter_for_float(&mut frame, src_frame, average, self.luma_scaling);
            }
        }
        Ok(frame.into())
    }
}