Function intensity_statistics 

pub fn intensity_statistics(
    values: &[f64],
    saturation_value: Option<f64>,
) -> Result<IntensityStatistics>