LaminoCam

A C++ library for tomographic reconstruction of laminographic data, with missing cone correction and Model-Based Iterative Reconstruction (MBIR) capabilities. Developed at Lawrence Berkeley National Laboratory.

Overview

LaminoCam is a high-performance library developed at Lawrence Berkeley National Laboratory for advanced reconstruction of laminographic data. It provides forward and backward projection operators, iterative reconstruction algorithms, and Model-Based Iterative Reconstruction. The library is optimized for performance using:

• OpenMP parallelization
• Intel TBB (Threading Building Blocks)
• FFTW for fast Fourier transforms
• FINUFFT for non-uniform FFTs
• Optional GPU acceleration via CUDA

Features

• Forward/Backward Projection: Efficient projection operators on polar grids
• Iterative Reconstruction: Conjugate gradient and Nesterov accelerated gradient methods
• MBIR: Model-Based Iterative Reconstruction L1 regularization via Split Bregman method.
• TIFF I/O: Read and write reconstruction data

Requirements

Build Dependencies

• CMake ≥ 3.20
• C++ compiler with C++20 support (llvm recommended)
• OpenMP
• Intel TBB
• FFTW3 (single and double precision)
• libtiff
• FINUFFT
• CUDA for GPU acceleration

Supported Platforms

• Linux (tested on Arch Linux)
• macOS (15.5+)

Building

Configure and Build

# Linux
cmake --preset linux-release
cmake --build --preset release
ctest --test-dir release --output-on-failure

# macOS
cmake --preset macos
cmake --build --preset macos

Usage

Running Reconstruction

The reconstruction tool uses a TOML configuration file to specify input data and parameters:

./build/recon_scalar <config.toml>

TOML Configuration File

Create a TOML file with the following structure:

[input]
filename = "/path/to/projections.tif"
angles = "/path/to/angles.txt" 

[output]
filename = "path/to/output.tiff"

[recon_params]
max_iters = 100
tol = 1e-5
xtol = 1e-5
recon_dims = [21, 511, 511]

[recon_params.regularizer]
method = "split_bregman"

# Parameters for split_bregman method
[recon_params.regularizer.split_bregman]
inner_iters = 1
lambda = 0.1
mu = 10

Notes:

• Angles are expected be in degrees
• The angles file should be is ASCII format and contain one angle per line
• A template configuration file (config_template.toml) is generated automatically if no input is provided

Documentation

The latest version is available in readthedocs: https://laminocam.readthedocs.io/en/latest/

License

Copyright (c) 2018, The Regents of the University of California, through Lawrence Berkeley National Laboratory (subject to receipt of any required approvals from the U.S. Dept. of Energy). All rights reserved.

If you have questions about your rights to use or distribute this software, please contact Berkeley Lab's Innovation & Partnerships Office at IPO@lbl.gov.

This Software was developed under funding from the U.S. Department of Energy and the U.S. Government consequently retains certain rights. As such, the U.S. Government has been granted for itself and others acting on its behalf a paid-up, nonexclusive, irrevocable, worldwide license in the Software to reproduce, distribute copies to the public, prepare derivative works, and perform publicly and display publicly, and to permit other to do so.

Contributing

Contact

For questions or issues, please contact Berkeley Lab's Innovation & Partnerships Office at IPO@lbl.gov.