taigw / GeodisTK

geodesic distance transform of 2d/3d images
MIT License
132 stars 35 forks source link

GeodisTK: Geodesic Distance Transform Toolkit for 2D and 3D Images

This repository provides source codes for geodesic distance transforms used in the following papers. If you use our code, please cite them.

@article{Wang2018Deepigeos,
    author={Wang, Guotai and Zuluaga, Maria A. and Li, Wenqi and Pratt, Rosalind and Patel, Premal A. and Aertsen, Michael and Doel, Tom and David, Anna L. and Deprest, Jan and Ourselin, Sébastien and Vercauteren, Tom},
    journal={IEEE Transactions on Pattern Analysis and Machine Intelligence}, 
    title={DeepIGeoS: A Deep Interactive Geodesic Framework for Medical Image Segmentation}, 
    year={2019},
    volume={41},
    number={7},
    pages={1559-1572}}

@article{LUO2021102102,
    author = {Luo, Xiangde and Wang, Guotai and Song, Tao and Zhang, Jingyang and Aertsen,Michael  and  Deprest, Jan and  Ourselin, Sebastien and Vercauteren, Tom and Zhang, Shaoting},
    title = {MIDeepSeg: Minimally interactive segmentation of unseen objects from medical images using deep learning},
    journal = {Medical Image Analysis},
    volume = {72},
    pages = {102102},
    year = {2021}}

Introduction

Geodesic transformation of images can be implementated with two approaches: fast marching and raster scan. Fast marching is based on the iterative propagation of a pixel front with velocity F [1]. Raster scan is based on kernel operations that are sequentially applied over the image in multiple passes [2][3]. In GeoS [4], the authors proposed to use a 3x3 kernel for forward and backward passes for efficient geodesic distance transform, which was used for image segmentation.

ranster scan Raster scan for geodesic distance transform. Image from [4]

DeepIGeoS [5] proposed to combine geodesic distance transforms with convolutional neural networks for efficient interactive segmentation of 2D and 3D images.

2D example A comparison of fast marching and ranster scan for 2D geodesic distance transform. (d) shows the Euclidean distance and (e) is a mixture of Geodesic and Euclidean distance.

This repository provides a cpp implementation of fast marching and raster scan for 2D/3D geodesic and Euclidean distance transforms and a mixture of them, and proivdes a python interface to use it.

How to install

  1. Install this toolkit easily by typing pip install GeodisTK

  2. Alternatively, if you want to build from source files, download this package and run the following:

    python setup.py build
    python setup.py install

How to use

  1. See a 2D example, run python demo2d.py

  2. See a 3D example, run python demo3d.py