Engineering:Project:Bayesian Segmentation

• [#BAYESIAN_CLASSIFIER_IMAGE_FILTER 1 BAYESIAN CLASSIFIER IMAGE FILTER]
  • [#Introduction 1.1 Introduction]
  • [#Use_Case 1.2 Use Case]
  • [#Data 1.3 Data]
  • [#Algorithm 1.4 Algorithm]
  • [#The_ITK_filter_design 1.5 The ITK filter design]
  • [#Project_Status 1.6 Project Status]
  • [#Contacts 1.7 Contacts]

BAYESIAN CLASSIFIER IMAGE FILTER

File:BayesianProgWeekProject.ppt

Introduction

This ITK filter is a segmentation algorithm which utilizes Bayes's Rule along with an affine-invarient anisotropic smoothing filter.

Use Case

I'd like to segment a volume or sub-volume into 'N' classes in a very general manner. I will provide the data and the number of classes that I expect and the algorithm will output a labelmap with 'N' classes.

Data

We are using the Harvard structural datasets.

Algorithm

1. The user sets the number of distinct classes for segmentation: 'N'
2. Generate 'N' prior images (default, 'N' uniform prior images)
   
3. Generate 'N' statistical distributions (default, 'N' normal distributions)
   
4. Generate 'N' membership images by applying the statistical distributions to the raw data
   
5. Generate 'N' posterior images by applying Bayes' rule to the prior and membership images
   
6. Smooth the posterior images for 'm' iterations using an affine-invarient anisotropic smoothing filter and renormalize after each iteration (default, m = 5)
   
7. Apply maximum a posteriori rule to apply labeling and finalize segmentation

The ITK filter design

Flowchart

Project Status

• Fully incorporated into itkBayesianClassificationImageFilter and itkBayesianClassificationInitializationImageFilter in ITK CVS.
  
• Fully wrapped in VTK for use in Slicer.
  
• The working ITK code has been committed to the SandBox.

Contacts

• John Melonakos @ Georgia Tech
• Luis Ibanez @ Kitware
• Karthik Krishnan @ Kitware