/FALCON

FALCON is a Python-based software application designed to facilitate PET motion correction, both for head and total-body scans. Our program is built around the fast 'greedy' registration toolkit, which serves as the registration engine. With FALCON, users can enjoy a streamlined experience for implementing motion correction.

Primary LanguagePythonGNU General Public License v3.0GPL-3.0

🦅 About FALCON 2.0

PyPI version Documentation Status Monthly Downloads Daily DownloadsDiscord YouTube

FALCON V2 (Fast Algorithms for motion correction) is an advanced, fully-automatic tool for motion correction in dynamic total-body or whole-body PET imaging. Designed with flexibility and reliability at its core, it's now even more versatile, capable of operating across various operating systems and architectures. 🚀

Shows an illustrated MOOSE story board adopted to different themes

🌟 Major Features

  • 🐍 Now Python-Powered for Effortless Integration: Say goodbye to cumbersome shell scripts! FALCON V2 is now a Python package, offering seamless compatibility with Python 3.9 and beyond. Integrate it effortlessly into your modern workflows and enjoy the unparalleled convenience.
  • 💻🖥️ Unveiling Cross-Platform Freedom: Break free from the Linux-only limitations! FALCON V2 now extends its robust support to Windows and Mac as well. A seamless, uniform experience across all platforms is no longer a wish—it's a reality!
  • 🛠️ Universal Architecture Compatibility with a Twist: Think you've seen FALCON V2 at its best? Wait until you experience it on ARM architecture! While designed to operate seamlessly on x86 and the latest M1 Silicon, it's on ARM where FALCON V2 truly soars. Prepare to be blown away by unparalleled performance!
  • 🚀 Say Goodbye to Memory Woes and Server Farms: Forget about expensive, beefy servers and memory constraints. With FALCON V2's optimized out-of-core computing, powered by Dask, we've redefined efficiency. Get ready to experience unparalleled speed and performance without breaking the bank or your system's RAM!
  • 🌍🔬 The Ultimate Flex: Versatility Reimagined: For the first time ever, a tool that knows no boundaries—meet FALCON V2. Whether you're dealing with any region, tracer, or modality, FALCON V2 is the only tool you'll ever need for your diagnostic and analytical adventures. Say goodbye to specialized solutions; this is the new standard - thanks to the amazing 'greedy' registration library!

Star History 🤩

Star History Chart

🚀 FALCON's motion correction in action

Shows an illustrated MOOSE story board adopted to different themes

In this analysis, we are examining the mean image of 20 dynamic frames of a 68Ga-PSMA study both before and after motion correction. By comparing the two mean images, we can clearly see the significant improvement that results from motion correction. The mean image after motion correction appears noticeably sharper and more defined than the one before correction.

🛠 Installation Guide

Virtual Environment Setup

Creating a virtual environment is highly recommended to avoid any potential conflicts with other Python packages.

  • Windows:

    python -m venv falconz_env
.\falconz_env\Scripts\activate

  • Linux/Mac:

    python3 -m venv falconz_env
source falconz_env/bin/activate

Installing FALCON V2

With your virtual environment activated, install FALCON V2 using pip:

pip install falconz==2.0.25 # stable recommended version

🚀 Usage

FALCON supports DICOM, Nifti, Analyze, and Metaimage file formats, whether it's a single 4D image or multiple 3D images. Simply specify the directory where the files are located and indicate the registration type. FALCON will take care of the rest.

To use FALCON, use the following syntax:

falconz -d path_to_4d_images -r <rigid | affine | deformable> -i <number_of_iterations_per_level> -sf <starting_frame_from_which_moco_should_be_performed> -rf <reference_frame>

Here's an example of using FALCON in Pro mode:

falconz -d /Documents/Sub001 -r deformable -i 100x50x25 -sf 0 -rf -1

In this example, FALCON is performing deformable registration with 100, 50, and 25 iterations at each level of the multi-scale registration. The registration starts from the 1st frame and uses the last frame as the reference.

Here's an example of using FALCON in lazy mode:

falconz -d /Documents/Sub001 -r deformable # for whole-body registration
falconz -d /Documents/Sub001 -r rigid # for brain only studies (much faster processing)

We also offer a specialized 🚀 Dash mode, engineered for rapid motion correction across total-body datasets. Execute complex whole-body registration tasks at unprecedented speeds ⚡ with a simple command! 👩‍💻👨‍💻

falconz -d /Documents/Sub001 -r deformable -m dash # for high-velocity whole-body registration

As shown above, you don't need to specify many additional parameters. The rest of the parameters are either inferred or set automatically based on common standards.

⚠️ Note: If you're not satisfied with the 'inferred' start frame, you can always set it on your own (the internal threshold is set to be quite safe). Refer to the manuscript for more information.

If you need help with FALCON or want to review the command-line options, use:

falconz --help

Please note that the number of iterations is specified as a string of values separated by 'x' in the -i option. For example, to perform 50 iterations at each level, you would use -i 50x50x50.

🗂 Required Folder Structure:

Required Folder Structure

FALCON accepts the following file formats for dynamic PET images: .dcm, .ima, .nii, .nii.gz, and .img/hdr.

Here are some examples to illustrate the accepted folder structures:

  1. For a bunch of DICOM (.dcm) or IMA (.ima) files:

    └── PET_WB_DYNAMIC_(QC)_0005  
    ├── XXX_1.dcm
    ├── XXX_2.dcm
    └── ...

  2. For a single 4D Nifti (.nii or .nii.gz) file:

    └── PET_WB_DYNAMIC_(QC)_0005  
    └── XXX.nii.gz

  3. For a bunch of 3D Nifti (.nii or .nii.gz) files:

    └── PET_WB_DYNAMIC_(QC)_0005  
    ├── XXX_1.nii.gz
    ├── XXX_2.nii.gz
    └── ...

  4. For Analyze (.img/hdr) files:

    └── PET_WB_DYNAMIC_(QC)_0005  
    ├── XXX.img
    ├── XXX.hdr
    └── ...

The main folder, PET_WB_DYNAMIC_(QC)_0005, should contain the dynamic PET images to be motion corrected.

🗂 Resultant Folder Structure

Upon successful execution, FALCON auto-generates an organized output directory, positioned at the same hierarchical level as your original dynamic PET image folder. This dedicated directory carries a unique naming schema that incorporates 'FALCONZ', the version number, and a timestamp for easy identification.

Here's a snapshot of the output folder structure:

FALCONZ-V02-2023-09-03-17-28-17/  # Automatically generated results folder
├── Motion-corrected-images       # Corrected dynamic PET images
├── ncc-images                    # Normalized Cross-Correlation images for start frame identification
├── Split-Nifti-files             # Individual 3D Nifti files
└── transforms                    # Transformation data for motion correction

Folder Components

🖼️ Motion-corrected-images

This is where you'll find the final dynamic PET images, now refined through motion correction procedures. 🌟

📊 ncc-images

A collection of Normalized Cross-Correlation images—these serve as essential tools 🛠️ for determining the most appropriate start frame for motion correction.

📁 Split-Nifti-files

This folder contains individual 3D Nifti files, which are crucial 🗝️ for conducting the motion correction operations.

🔄 transforms

This section archives the warp fields in cases of deformable registration and the transformation matrices 📐 for rigid or affine registrations, allowing for transparency and potential reusability of these parameters.

FALCON doesn't just deliver high-precision motion-corrected images; it also provides a comprehensive, organized output structure 🗂️ designed for immediate utility and future analysis. 🚀

📊 FalconZ Performance Metrics

The table below outlines FalconZ's performance metrics across various operating systems, architectures, dynamic PET datasets, and system hardware:

🖥️ OS / Arch Hardware Dataset Voxel Size (mm³) Matrix Size Time (Minutes)
🪟 Windows x86 ... 🧠 Brain ... ... ...
❤️ Cardiac ... ... ...
🕺 Total-Body 4x4x4 150x150x486 ...
🐧 Linux x86 Intel 20 cores, 60 GB RAM 🧠 Brain ... ... ...
❤️ Cardiac ... ... ...
🕺 Total-Body 4x4x4 150x150x486 12 min
🍏 macOS x86 Intel 4 cores, 32 GB RAM 🧠 Brain ... ... ...
❤️ Cardiac ... ... ...
🕺 Total-Body 4x4x4 150x150x486 8 min
🍏 macOS ARM64 M2 8 Cores, 8GB RAM 🧠 Brain ... ... ...
❤️ Cardiac ... ... ...
🕺 Total-Body 4x4x4 150x150x486 3 min

❤️ Citations

If you find FALCON useful, please consider citing the following publications:

1. Lalith Kumar Shiyam Sundar, et al. “Fully-automated, fast motion correction of dynamic whole and total-body PET/CT imaging studies.” JNM, 2023.
2. Gutschmayer S, et al. “A scale space theory based motion correction approach for dynamic PET brain imaging studies.” Frontiers in Physics, vol. 10, 2022.
3. Venet L, et al. “Accurate and Robust Alignment of Differently Stained Histologic Images Based on Greedy Diffeomorphic Registration.” Appl Sci, vol. 11, 2021.

Thank you for choosing FALCON V2. Let's soar to new heights together! 🦅