ITK/Examples/Statistics/MembershipSample

From KitwarePublic
< ITK‎ | Examples
Revision as of 23:22, 2 March 2011 by Lorensen (talk | contribs)
Jump to navigationJump to search

A MembershipSample is a list ListSample with an additional attribute indicating which class the sample belongs to. This is the output type in many of the clustering and segmentation algorithms in the Statistics framework.

MembershipSample.cxx

<source lang="cpp">

  1. include "itkListSample.h"
  2. include "itkMembershipSample.h"
  3. include "itkVector.h"

int main(int, char *[]) {

 typedef itk::Vector< float, 3 > MeasurementVectorType;
 typedef itk::Statistics::ListSample< MeasurementVectorType > SampleType;
 SampleType::Pointer sample = SampleType::New();
 MeasurementVectorType mv;
 mv[0] = 1.0;
 mv[1] = 2.0;
 mv[2] = 4.0;
 sample->PushBack(mv);
 mv[0] = 2.0;
 mv[1] = 4.0;
 mv[2] = 5.0;
 sample->PushBack(mv);
 mv[0] = 3.0;
 mv[1] = 8.0;
 mv[2] = 6.0;
 sample->PushBack(mv);
 typedef itk::Statistics::MembershipSample< SampleType >
   MembershipSampleType;
 MembershipSampleType::Pointer membershipSample =
   MembershipSampleType::New();
 membershipSample->SetSample(sample);
 membershipSample->SetNumberOfClasses(2);
 membershipSample->AddInstance(0U, 0UL );
 membershipSample->AddInstance(0U, 1UL );
 membershipSample->AddInstance(1U, 2UL );
 MembershipSampleType::ConstIterator iter = membershipSample->Begin();
 while ( iter != membershipSample->End() )
   {
   std::cout << "instance identifier = " << iter.GetInstanceIdentifier()
             << "\t measurement vector = "
             << iter.GetMeasurementVector()
             << "\t frequency = "
             << iter.GetFrequency()
             << "\t class label = "
             << iter.GetClassLabel()
             << std::endl;
   ++iter;
   }


 MembershipSampleType::ClassSampleType::ConstPointer classSample =
   membershipSample->GetClassSample( 0 );
 MembershipSampleType::ClassSampleType::ConstIterator c_iter =
   classSample->Begin();
 while ( c_iter != classSample->End() )
   {
   std::cout << "instance identifier = " << c_iter.GetInstanceIdentifier()
             << "\t measurement vector = "
             << c_iter.GetMeasurementVector()
             << "\t frequency = "
             << c_iter.GetFrequency() << std::endl;
   ++c_iter;
   }
 return EXIT_SUCCESS;

} </source>


{{ITK CMakeLists|MembershipSample}}