objscale 0.1.6

Object-based analysis functions for fractal dimensions and size distributions

objscale

Object-based analysis functions for fractal dimensions and size distributions in atmospheric sciences and beyond. Optimized for large datasets.

Description

objscale provides computational tools for analyzing the scaling properties of objects in 2D binary arrays. The package consolidates methods for calculating size distributions and fractal dimensions that account for finite domain effects and complex boundary conditions. Originally developed for atmospheric science applications, these methods apply broadly to any field where object scaling properties matter.

The package implements methods from two main papers:

• DeWitt & Garrett (2024) - finite domain effects in size distributions
• DeWitt et al. (in prep) - fractal dimensions for cloud field characterization

Key Functions

finite_array_powerlaw_exponent

Calculate power-law exponents for size distributions while accounting for finite domain truncation effects. Essential for accurate scaling analysis in bounded domains.

individual_fractal_dimension

Fractal dimension of individual objects using the perimeter-area relationship, with proper handling of interior holes and resolution effects.

ensemble_correlation_dimension

Correlation dimension for characterizing the collective scaling properties of object ensembles. Uses point-pair correlation analysis across multiple length scales.

ensemble_box_dimension

Box-counting dimension for object ensembles. New analyses should prefer ensemble_correlation_dimension. Counts boxes containing object boundaries at varying spatial scales.

ensemble_coarsening_dimension

Novel fractal dimension based on how total object perimeter changes under spatial coarsening operations. New analyses should prefer ensemble_correlation_dimension. Used in Rees et. al., 2024

Installation

pip install objscale

Documentation

📖 Full Documentation

Complete API reference, detailed examples, and usage guides are available at objscale.readthedocs.io.

Quick Example

import objscale
import numpy as np

# Create binary array (e.g., cloud mask, percolation lattice)
arrays = [(np.random.random((1000, 1000)) < 0.3).astype(int) for _ in range(4)]

# Size distribution with finite domain corrections
(exponent, error), (log10_sizes, log10_counts) = objscale.finite_array_powerlaw_exponent(
    arrays, 'area', return_counts=True
)

# Ensemble fractal dimensions
corr_dim, corr_error = objscale.ensemble_correlation_dimension(arrays)
box_dim, box_error = objscale.ensemble_box_dimension(arrays)

# Individual object analysis  
ind_dim, ind_error = objscale.individual_fractal_dimension(arrays)

Features

• Finite domain corrections: Proper handling of truncation effects at domain boundaries as recommended by [DeWitt & Garrett (2024)](ACP - Finite domains cause bias in measured and modeled distributions of cloud sizes)
• Multiple size metrics: Area, perimeter, width, height, nested perimeter
• Arbitrary boundaries: Support for NaN-demarcated non-rectangular domains
• Individual and Ensemble methods: Characterize both individual and collective properties of object fields
• Performance optimized: Numba acceleration for computational efficiency. Can handle billions of individual objects on a mid-range laptop.

Requirements

• Python ≥ 3.8
• NumPy ≥ 1.20.0
• SciPy ≥ 1.7.0
• scikit-image ≥ 0.18.0
• Numba ≥ 0.56.0

Available Functions

Fractal Dimensions

• individual_fractal_dimension - Fractal dimension of individual objects
• ensemble_correlation_dimension - Correlation dimension for object ensembles
• ensemble_box_dimension - Box-counting dimension for object ensembles
• ensemble_coarsening_dimension - Coarsening-based fractal dimension

Size Distributions

• finite_array_powerlaw_exponent - Power-law exponents with finite domain corrections
• finite_array_size_distribution - Size distributions with truncation analysis
• array_size_distribution - Basic size distribution for single arrays

Object Analysis

• get_structure_props - Calculate perimeter, area, width, height of structures
• total_perimeter - Total perimeter of all objects
• total_number - Count number of structures
• isolate_largest_structure - Extract the largest connected structure
• remove_structures_touching_border_nan - Remove border-touching structures
• remove_structure_holes - Fill holes in structures
• clear_border_adjacent - Clear structures touching array edges

Utilities

• coarsen_array - Coarsen array resolution by averaging
• linear_regression - Linear regression with error estimates
• encase_in_value - Add border of specified value around array

For detailed parameter descriptions and usage examples, see the full documentation or use help(objscale.function_name) or objscale.function_name? in IPython/Jupyter.

Support Statement

This package consolidates research code developed over several years. While functional and tested, it should be considered research software with limited ongoing support. Users are encouraged to understand the underlying methods through the referenced papers before applying to their data.

References

If you use this package, please cite:

DeWitt, T. D. and Garrett, T. J.: Finite domains cause bias in measured and modeled distributions of cloud sizes, Atmos. Chem. Phys., 24, 8457–8472, https://doi.org/10.5194/acp-24-8457-2024, 2024.

Author

Thomas D. DeWitt University of Utah Department of Atmospheric Sciences

Sonnet 4 with Claude Code Anthropic

License

MIT License