cdmurph32/fast_image_resize

Rust library for fast image resizing with using of SIMD instructions.

fast_image_resize

Rust library for fast image resizing with using of SIMD instructions.

CHANGELOG

Supported pixel formats and available optimisations:

Format	Description	Native Rust	SSE4.1	AVX2	Neon
U8	One u8 component per pixel (e.g. L)	+	+	+	+
U8x2	Two u8 components per pixel (e.g. LA)	+	+	+	+
U8x3	Three u8 components per pixel (e.g. RGB)	+	+	+	+
U8x4	Four u8 components per pixel (e.g. RGBA, RGBx, CMYK)	+	+	+	+
U16	One u16 components per pixel (e.g. L16)	+	+	+	+
U16x2	Two u16 components per pixel (e.g. LA16)	+	+	+	+
U16x3	Three u16 components per pixel (e.g. RGB16)	+	+	+	+
U16x4	Four u16 components per pixel (e.g. RGBA16, RGBx16, CMYK16)	+	+	+	+
I32	One i32 component per pixel	+	-	-	-
F32	One f32 component per pixel	+	-	-	-

Colorspace

Resizer from this crate does not convert image into linear colorspace during resize process. If it is important for you to resize images with a non-linear color space (e.g. sRGB) correctly, then you need to convert it to a linear color space before resizing and convert back a color space of result image. Read more about resizing with respect to color space.

This crate provides the PixelComponentMapper structure that allows you to create colorspace converters for images whose pixels based on u8 and u16 components.

In addition, the crate contains functions create_gamma_22_mapper() and create_srgb_mapper() to create instance of PixelComponentMapper that converts images from sRGB or gamma 2.2 into linear colorspace and back.

Some benchmarks

• All x86-64 benchmarks.
• All arm64 benchmarks.

Rust libraries used to compare of resizing speed:

• image (https://crates.io/crates/image)
• resize (https://crates.io/crates/resize)

Resize RGB8 image (U8x3) 4928x3279 => 852x567

Pipeline:

src_image => resize => dst_image

• Source image nasa-4928x3279.png
• Numbers in table is mean duration of image resizing in milliseconds.

	Nearest	Bilinear	CatmullRom	Lanczos3
image	20.66	82.74	141.56	199.96
resize	-	48.75	97.42	145.70
fir rust	0.28	39.57	67.00	98.50
fir sse4.1	0.28	9.63	14.13	19.84
fir avx2	0.28	7.73	9.67	14.66

Resize RGBA8 image (U8x4) 4928x3279 => 852x567

Pipeline:

src_image => multiply by alpha => resize => divide by alpha => dst_image

• Source image nasa-4928x3279-rgba.png
• Numbers in table is mean duration of image resizing in milliseconds.
• The image crate does not support multiplying and dividing by alpha channel.

	Nearest	Bilinear	CatmullRom	Lanczos3
resize	-	65.74	130.18	194.13
fir rust	0.19	35.61	52.50	75.76
fir sse4.1	0.19	13.19	17.25	22.62
fir avx2	0.19	9.57	11.97	16.37

Resize L8 image (U8) 4928x3279 => 852x567

Pipeline:

src_image => resize => dst_image

• Source image nasa-4928x3279.png has converted into grayscale image with one byte per pixel.
• Numbers in table is mean duration of image resizing in milliseconds.

	Nearest	Bilinear	CatmullRom	Lanczos3
image	16.97	48.47	77.55	107.48
resize	-	17.10	35.46	60.94
fir rust	0.15	13.44	14.73	22.58
fir sse4.1	0.15	4.99	5.28	8.05
fir avx2	0.15	7.09	5.18	8.60

Examples

Resize RGBA8 image

use std::io::BufWriter;
use std::num::NonZeroU32;

use image::codecs::png::PngEncoder;
use image::io::Reader as ImageReader;
use image::{ColorType, ImageEncoder};

use fast_image_resize as fr;

fn main() {
    // Read source image from file
    let img = ImageReader::open("./data/nasa-4928x3279.png")
        .unwrap()
        .decode()
        .unwrap();
    let width = NonZeroU32::new(img.width()).unwrap();
    let height = NonZeroU32::new(img.height()).unwrap();
    let mut src_image = fr::Image::from_vec_u8(
        width,
        height,
        img.to_rgba8().into_raw(),
        fr::PixelType::U8x4,
    ).unwrap();

    // Multiple RGB channels of source image by alpha channel 
    // (not required for the Nearest algorithm)
    let alpha_mul_div = fr::MulDiv::default();
    alpha_mul_div
        .multiply_alpha_inplace(&mut src_image.view_mut())
        .unwrap();

    // Create container for data of destination image
    let dst_width = NonZeroU32::new(1024).unwrap();
    let dst_height = NonZeroU32::new(768).unwrap();
    let mut dst_image = fr::Image::new(
        dst_width,
        dst_height,
        src_image.pixel_type(),
    );

    // Get mutable view of destination image data
    let mut dst_view = dst_image.view_mut();

    // Create Resizer instance and resize source image
    // into buffer of destination image
    let mut resizer = fr::Resizer::new(
        fr::ResizeAlg::Convolution(fr::FilterType::Lanczos3),
    );
    resizer.resize(&src_image.view(), &mut dst_view).unwrap();

    // Divide RGB channels of destination image by alpha
    alpha_mul_div.divide_alpha_inplace(&mut dst_view).unwrap();

    // Write destination image as PNG-file
    let mut result_buf = BufWriter::new(Vec::new());
    PngEncoder::new(&mut result_buf)
        .write_image(
            dst_image.buffer(),
            dst_width.get(),
            dst_height.get(),
            ColorType::Rgba8,
        )
        .unwrap();
}

Resize with cropping

use std::num::NonZeroU32;

use image::codecs::png::PngEncoder;
use image::io::Reader as ImageReader;
use image::{ColorType, GenericImageView};

use fast_image_resize as fr;

fn resize_image_with_cropping(
    mut src_view: fr::DynamicImageView,
    dst_width: NonZeroU32,
    dst_height: NonZeroU32
) -> fr::Image {
    // Set cropping parameters
    src_view.set_crop_box_to_fit_dst_size(dst_width, dst_height, None);

    // Create container for data of destination image
    let mut dst_image = fr::Image::new(
        dst_width,
        dst_height,
        src_view.pixel_type(),
    );
    // Get mutable view of destination image data
    let mut dst_view = dst_image.view_mut();

    // Create Resizer instance and resize source image
    // into buffer of destination image
    let mut resizer = fr::Resizer::new(
        fr::ResizeAlg::Convolution(fr::FilterType::Lanczos3)
    );
    resizer.resize(&src_view, &mut dst_view).unwrap();

    dst_image
}

fn main() {
    let img = ImageReader::open("./data/nasa-4928x3279.png")
        .unwrap()
        .decode()
        .unwrap();
    let width = NonZeroU32::new(img.width()).unwrap();
    let height = NonZeroU32::new(img.height()).unwrap();
    let src_image = fr::Image::from_vec_u8(
        width,
        height,
        img.to_rgb8().into_raw(),
        fr::PixelType::U8x3,
    ).unwrap();
    resize_image_with_cropping(
        src_image.view(),
        NonZeroU32::new(1024).unwrap(),
        NonZeroU32::new(768).unwrap(),
    );
}

Change CPU extensions used by resizer

use fast_image_resize as fr;

fn main() {
    let mut resizer = fr::Resizer::new(
        fr::ResizeAlg::Convolution(fr::FilterType::Lanczos3),
    );
    unsafe {
        resizer.set_cpu_extensions(fr::CpuExtensions::Sse4_1);
    }
}