#include "itkBayesianClassifierInitializationImageFilter.h"
#include "itkBayesianClassifierImageFilter.h"
//Bayesian Segmentation Approach
typedef itk::BayesianClassifierInitializationImageFilter< InputImageType > BayesianInitializerType;
typename BayesianInitializerType::Pointer bayesianInitializer = BayesianInitializerType::New();
bayesianInitializer->SetInput(reader->GetOutput());
bayesianInitializer->SetNumberOfClasses(numClass);// Background, WM, GM and CSF
bayesianInitializer->Update();
typedef float PriorType;
typedef float PosteriorType;
typedef itk::VectorImage< InputPixelType, Dimension > VectorInputImageType;
typedef itk::BayesianClassifierImageFilter< VectorInputImageType, OutputPixelType, PosteriorType, PriorType > ClassifierFilterType;
typename ClassifierFilterType::Pointer bayesClassifier = ClassifierFilterType::New();
bayesClassifier->SetInput(bayesianInitializer->GetOutput());
bayesClassifier->Update();
体会:相比kmean分类,比较慢