ihist 0.1.1

Fast integer image histograms

ihist

Fast histogram computation for image data.

Currently experimental and API may still change.

Python API

Quick Start

import numpy as np
import ihist

# Grayscale image
image = np.random.randint(0, 256, (100, 100), dtype=np.uint8)
hist = ihist.histogram(image)  # Shape: (256,)

# RGB image
rgb = np.random.randint(0, 256, (100, 100, 3), dtype=np.uint8)
hist = ihist.histogram(rgb)  # Shape: (3, 256)

# With masking
mask = np.ones((100, 100), dtype=np.uint8)
hist = ihist.histogram(image, mask=mask)

Function

The Python API provides a single function for computing histograms:

import ihist

histogram = ihist.histogram(image, bits=None, mask=None,
                            components=None, out=None,
                            accumulate=False, parallel=True)

Parameters

image : array_like Input image data. Must be uint8 or uint16, and 1D, 2D, or 3D.

• 1D arrays (W,) are interpreted as (1, W, 1)
• 2D arrays (H, W) are interpreted as (H, W, 1)
• 3D arrays (H, W, C) use C as number of components

Total pixel count must not exceed 2^32-1.

bits : int, optional Number of significant bits per sample. If not specified, defaults to full depth (8 for uint8, 16 for uint16). Valid range: [0, 8] for uint8, [0, 16] for uint16.

mask : array_like, optional Per-pixel mask. Must be uint8. Shape must match image dimensions:

• For 1D images (W,): mask must be 1D, shape (W,)
• For 2D/3D images: mask must be 2D, shape (H, W)

Only pixels with non-zero mask values are included. If not specified, all pixels are included.

components : sequence of int, optional Indices of components to histogram. If not specified, all components are histogrammed. Each index must be in range [0, n_components).

Examples:

• [0] - histogram only the first component (e.g., red in RGB)
• [0, 1, 2] - histogram first three components (e.g., RGB in RGBA, skipping alpha)

out : array_like, optional Pre-allocated output buffer. Must be uint32, and either:

• 1D with shape (2^bits,) for single-component histogram, or
• 2D with shape (n_hist_components, 2^bits)

If not specified, a new array is allocated and returned.

accumulate : bool, optional If False (default), the output buffer is zeroed before computing the histogram. If True, histogram values are accumulated into the existing buffer values. No effect if out is not given.

parallel : bool, optional If True (default), allows automatic multi-threaded execution for large images. If False, guarantees single-threaded execution.

Returns

histogram : ndarray Histogram(s) as uint32 array.

• If the image is 1D or 2D and components is not specified, returns 1D array of shape (2^bits,)
• If the image is 3D or components is explicitly specified, returns 2D array of shape (n_hist_components, 2^bits)
• If out was provided, returns the same array after filling

Java API

Quick Start

import io.github.marktsuchida.ihist.HistogramRequest;
import java.nio.IntBuffer;

// Grayscale image
byte[] image = new byte[100 * 100];
IntBuffer hist = HistogramRequest.forImage(image, 100, 100).compute(); // 256 bins
// hist.remaining() == 256, access with hist.get(i)

// RGB image
byte[] rgb = new byte[100 * 100 * 3];
IntBuffer hist = HistogramRequest.forImage(rgb, 100, 100, 3).compute(); // 3 * 256 bins

// With advanced options
IntBuffer hist = HistogramRequest.forImage(image, 100, 100)
    .roi(10, 10, 80, 80)       // Region of interest
    .mask(maskData, 80, 80)    // Per-pixel mask
    .bits(8)                   // Significant bits
    .parallel(true)            // Allow multi-threading
    .compute();

Classes

HistogramRequest - Builder-style interface:

// Multi-component image (e.g., RGB)
IntBuffer hist = HistogramRequest.forImage(image, width, height, 3)
    .selectComponents(0, 1, 2)      // Which components to histogram
    .roi(x, y, roiWidth, roiHeight) // Region of interest
    .mask(maskData, maskWidth, maskHeight)  // Per-pixel mask with dimensions
    .maskOffset(offsetX, offsetY)   // Mask offset for ROI alignment
    .bits(sampleBits)               // Significant bits per sample
    .output(preallocatedBuffer)     // Pre-allocated output (size must be exact)
    .accumulate(true)               // Add to existing values
    .parallel(true)                 // Allow multi-threading
    .compute();

// All methods except for forImage() and compute() are optional.

// The returned IntBuffer has position/limit set to cover the histogram data.
// If output(IntBuffer) was used, the returned buffer is a duplicate that shares
// storage with the original; the original's position/limit are not modified.
// For heap buffers (default or when output(int[]) was used), get the array:
int[] array = hist.array();

IHistNative - Low-level JNI wrapper (advanced):

// Buffer-based (direct or array-backed buffers)
IHistNative.histogram8(sampleBits, imageBuffer, maskBuffer,
    height, width, imageStride, maskStride, nComponents, componentIndices,
    histogramBuffer, parallel);

(And, similarly, histogram16() for 16-bit images.)

Input Types

The Java API supports both arrays and NIO buffers:

• 8-bit images: byte[] or ByteBuffer
• 16-bit images: short[] or ShortBuffer
• Mask: byte[] or ByteBuffer
• Histogram output: int[] or IntBuffer

Performance: Zero-Copy vs Copy

Input Type	IHistNative	HistogramRequest	Notes
Array (byte[], short[], int[])	N/A (wrap first)	Zero-copy	Wrapped as heap buffer
Direct buffer	Zero-copy	Zero-copy	GetDirectBufferAddress
Array-backed buffer (ByteBuffer.wrap())	Zero-copy	Zero-copy	GetPrimitiveArrayCritical
View buffer (e.g., asReadOnlyBuffer())	Rejected	Copy	Copied to temp direct buffer
Read-only heap buffer	Rejected	Copy (input only)	Rejected for histogram output

IHistNative requires buffers to be either direct or array-backed. View buffers and other buffer types are rejected with IllegalArgumentException.

HistogramRequest handles all buffer types automatically. Arrays are wrapped in heap buffers (zero-copy). Unsupported buffer types are silently copied to temporary direct buffers, which incurs a performance overhead but ensures all input types work.

Notes

• Signed types: Java byte is signed (-128 to 127), but pixel values are interpreted as unsigned (0 to 255). The native code correctly handles this. Similarly for short with 16-bit images.

• Accumulation: Like the C API, histogram values are accumulated. Use .accumulate(false) (default) to zero the output first.

• Exact buffer sizes: Buffer and array parameters must have exactly the required size (not just at least the required size). For the image buffer, this is width * height * nComponents; for the mask, width * height; for the histogram output, nHistComponents * 2^bits. For buffers, the size that must match is the value of remaining() (distance from position() to limit()), not capacity().

• Thread safety: The histogram functions are thread-safe for independent calls. Multiple threads can compute histograms simultaneously.

C API

The C API provides two functions for computing histograms of 2D image data:

• ihist_hist8_2d() - for 8-bit samples (uint8_t)
• ihist_hist16_2d() - for 16-bit samples (uint16_t)

Both functions have identical behavior except for the sample data type.

Function Signatures

#include <ihist/ihist.h>

void ihist_hist8_2d(
    size_t sample_bits,
    uint8_t const *restrict image,
    uint8_t const *restrict mask,
    size_t height,
    size_t width,
    size_t image_stride,
    size_t mask_stride,
    size_t n_components,
    size_t n_hist_components,
    size_t const *restrict component_indices,
    uint32_t *restrict histogram,
    bool maybe_parallel);

void ihist_hist16_2d(
    size_t sample_bits,
    uint16_t const *restrict image,
    uint8_t const *restrict mask,
    size_t height,
    size_t width,
    size_t image_stride,
    size_t mask_stride,
    size_t n_components,
    size_t n_hist_components,
    size_t const *restrict component_indices,
    uint32_t *restrict histogram,
    bool maybe_parallel);

Overview

These functions compute histograms for one or more components (stored as interleaved multi-sample pixels) from image data. They support:

• Multi-component images (e.g., grayscale, RGB, RGBA)
• Selective histogramming of specific components
• Optional per-pixel masking
• Region of interest (ROI) processing via stride and pointer offset
• Automatic parallelization for large images
• Arbitrary bit depths (not just full 8 or 16 bits)

Parameters

sample_bits Number of significant bits per sample. Valid range: 1-8 for ihist_hist8_2d(), 1-16 for ihist_hist16_2d(). The histogram will contain 2^sample_bits bins per sample.

Values with bits set beyond sample_bits are discarded and not counted in any bin.

image Pointer to image data. Samples are interleaved in row-major order:

• Row 0, pixel 0: all samples
• Row 0, pixel 1: all samples
• ...
• Row 1, pixel 0: all samples
• ...

May be NULL if height or width is 0.

mask (optional) Per-pixel mask for selective histogramming. If non-NULL, must point to height * mask_stride uint8_t values. Only pixels where the corresponding mask value is non-zero are included in the histogram.

Pass NULL to histogram all pixels.

height Image height in pixels. May be 0 for empty input.

width Image width in pixels. May be 0 for empty input.

image_stride Row stride for the image in pixels (not bytes). Must be ≥ width.

When image_stride equals width, the image is treated as contiguous. Use image_stride > width together with an offset image pointer to process a rectangular region of interest (ROI) within a larger image.

mask_stride Row stride for the mask in pixels (not bytes). Must be ≥ width.

When mask_stride equals width, the mask is treated as contiguous. This parameter allows the mask to have a different stride from the image, which is useful, for example, if you have a mask that covers only a rectanglular ROI of the image.

n_components Number of interleaved per pixel. Examples:

• 1 for grayscale
• 3 for RGB
• 4 for RGBA

Must be > 0.

n_hist_components Number of components to histogram. Must be > 0.

This allows histogramming a subset of components, such as skipping the alpha component in RGBA images.

component_indices Array of n_hist_components indices specifying which components to histogram. Each index must be in the range [0, n_components).

Examples:

• {0} - histogram only the first component (e.g., red in RGB)
• {0, 1, 2} - histogram first three components (e.g., RGB in RGBA, skipping alpha)
• {1, 2, 3} - histogram last three components (e.g., skip first component in ARGB)

Must not be NULL.

histogram (output, accumulated) Output buffer for histogram data. Must point to n_hist_components * 2^sample_bits uint32_t values.

Histograms for each component are stored consecutively:

• Bins for component_indices[0]: histogram[0] to histogram[2^sample_bits - 1]
• Bins for component_indices[1]: histogram[2^sample_bits] to histogram[2 * 2^sample_bits - 1]
• ...

Important: The histogram is accumulated into this buffer. Existing values are added to, not replaced. To obtain a fresh histogram, zero-initialize the buffer before calling the function.

maybe_parallel Controls parallelization.

• true - Allows automatic multi-threaded execution for large images, if ihist was built with parallelization support (TBB).
• false - Guarantees single-threaded execution.