Overview

The resampler library enables fast differentiable resizing and warping of arbitrary grids. It supports:

grids of any dimension (e.g., 1D, 2D images, 3D video, 4D batches of videos), containing
samples of any shape (e.g., scalars, colors, motion vectors, Jacobian matrices) and
any numeric type (e.g., uint8, float64, complex128)
within several array libraries (numpy, torch, and jax);
either 'dual' ("half-integer") or 'primal' grid-type for each dimension;
many boundary rules, specified per dimension, extensible via subclassing;
an extensible set of filter kernels, selectable per dimension;
optional gamma transfer functions for correct linear-space filtering;
prefiltering for accurate antialiasing when resize downsampling;
efficient backpropagation of gradients for torch and jax;
few dependencies (only numpy and scipy) and no C extension code, yet
faster resizing than C++ implementation in torch.nn.

A key strategy is to leverage existing sparse matrix representations and operations.

Example usage

!pip install -q mediapy numpy resampler

import mediapy as media
import numpy as np
import resampler

array = np.random.default_rng(1).random((4, 6, 3))  # 4x6 RGB image.
upsampled = resampler.resize(array, (128, 192))  # To 128x192 resolution.
media.show_images({'4x6': array, '128x192': upsampled}, height=128)

image = media.read_image('https://github.com/hhoppe/data/raw/main/image.png')
downsampled = resampler.resize(image, (32, 32))
media.show_images({'128x128': image, '32x32': downsampled}, height=128)

import matplotlib.pyplot as plt

array = np.array([3.0, 5.0, 8.0, 7.0])  # 4 source samples in 1D.
new_dual = resampler.resize(array, (32,))  # (default gridtype='dual') 8x resolution.
new_primal = resampler.resize(array, (25,), gridtype='primal')  # 8x resolution.

_, axs = plt.subplots(1, 2, figsize=(7, 1.5))
axs[0].set_title("gridtype='dual'")
axs[0].plot((np.arange(len(array)) + 0.5) / len(array), array, 'o')
axs[0].plot((np.arange(len(new_dual)) + 0.5) / len(new_dual), new_dual, '.')
axs[1].set_title("gridtype='primal'")
axs[1].plot(np.arange(len(array)) / (len(array) - 1), array, 'o')
axs[1].plot(np.arange(len(new_primal)) / (len(new_primal) - 1), new_primal, '.')
plt.show()

batch_size = 4
batch_of_images = media.moving_circle((16, 16), batch_size)
upsampled = resampler.resize(batch_of_images, (batch_size, 64, 64))
media.show_videos({'original': batch_of_images, 'upsampled': upsampled}, fps=1)

original

upsampled

Most examples above use the default resize() settings:

gridtype='dual' for both source and destination arrays,
boundary='auto' which uses 'reflect' for upsampling and 'clamp' for downsampling,
filter='lanczos3' (a Lanczos kernel with radius 3),
gamma=None which by default uses the 'power2' transfer function for the uint8 image in the second example,
scale=1.0, translate=0.0 (no domain transformation),
default precision and output dtype.

Advanced usage:

Map an image to a wider grid using custom scale and translate vectors, with horizontal 'reflect' and vertical 'natural' boundary rules, providing a constant value for the exterior, using different filters (Lanczos and O-MOMS) in the two dimensions, disabling gamma correction, performing computations in double-precision, and returning an output array in single-precision:

new = resampler.resize(
    image, (128, 512), boundary=('natural', 'reflect'), cval=(0.2, 0.7, 0.3),
    filter=('lanczos3', 'omoms5'), gamma='identity', scale=(0.8, 0.25),
    translate=(0.1, 0.35), precision='float64', dtype='float32')
media.show_images({'image': image, 'new': new})

Warp an image by transforming it using polar coordinates:

shape = image.shape[:2]
yx = ((np.indices(shape).T + 0.5) / shape - 0.5).T  # [-0.5, 0.5]^2
radius, angle = np.linalg.norm(yx, axis=0), np.arctan2(*yx)
angle += (0.8 - radius).clip(0, 1) * 2.0 - 0.6
coords = np.dstack((np.sin(angle) * radius, np.cos(angle) * radius)) + 0.5
resampled = resampler.resample(image, coords, boundary='constant')
media.show_images({'image': image, 'resampled': resampled})

Limitations:

Filters are assumed to be separable.
Although resize implements prefiltering, resample does not yet have it (and therefore may have aliased results if downsampling).
Differentiability is only with respect to the grid values, not wrt the resize shape, scale, translation, or the resampling coordinates.

Release files for resampler 1.1.0

For a detailed explanation of source distributions (sdists) and built distributions (wheels), please see the package formats documentation.

Source distribution (sdist)

Source distribution for resampler 1.1.0
File	Size	Uploaded
resampler-1.1.0.tar.gz	48.3 kB	Aug 6, 2026	Details

Built distribution (wheel)

Table of built distributions (wheels) for resampler 1.1.0
File	Interpreter	ABI	Platform	Reset
resampler-1.1.0-py3-none-any.whl 45.8 kB Aug 6, 2026	Python 3	none	any	Details

Total release size: 94.0 kB

Release files / resampler-1.1.0.tar.gz

Download URL	resampler-1.1.0.tar.gz
Size	48.3 kB
Tags	Source
SHA-256 checksum How to use checksums	`141d4cf152c20e364525e9b41d51842a8357354ddf214cdad1b784f1cf44dbfd`
BLAKE2b-256 checksum How to use checksums	`101b48be63ac1a431b4923d98f935412ea6ac844fb571fbb0d8e320dcb59226d`
Upload date	Aug 6, 2026
Uploaded using Trusted Publishing? What is trusted publishing?	No
Uploaded via	`twine/7.0.0 CPython/3.13.14`

Release files / resampler-1.1.0-py3-none-any.whl

Download URL	resampler-1.1.0-py3-none-any.whl
Size	45.8 kB
Tags	Python 3
SHA-256 checksum How to use checksums	`626775d90c2d57782c866dc17be67b0e612c4b8095ffaeeee98faf50372451f7`
BLAKE2b-256 checksum How to use checksums	`c73710331bfe73315c6c1e7616e3a773f9f5014e45ab00782954133e43ddc1c5`
Upload date	Aug 6, 2026
Uploaded using Trusted Publishing? What is trusted publishing?	No
Uploaded via	`twine/7.0.0 CPython/3.13.14`