typedef itk::LabelMap< ShapeLabelObjectType > LabelMapType;
typedef itk::ConnectedComponentImageFilter <MaskInputImageType, OutputImageType > ConnectedComponentImageFilterType;
typedef itk::LabelImageToShapeLabelMapFilter< OutputImageType, LabelMapType> I2LType;
typedef itk::LabelMap< ShapeLabelObjectType > LabelMapType;
typedef itk::ConnectedComponentImageFilter <MaskInputImageType, OutputImageType > ConnectedComponentImageFilterType;
typedef itk::LabelImageToShapeLabelMapFilter< OutputImageType, LabelMapType> I2LType;
typename ConnectedComponentImageFilterType::Pointer connected = ConnectedComponentImageFilterType::New ();
connected->SetInput(mask);
connected->FullyConnectedOn();
connected->Update();
typename ConnectedComponentImageFilterType::Pointer connected = ConnectedComponentImageFilterType::New ();
connected->SetInput(mask);
connected->FullyConnectedOn();
connected->Update();
//Create a map to contain all connectedComponent (even a little pixel)
typename I2LType::Pointer i2l = I2LType::New();
i2l->SetInput( connected->GetOutput() );
i2l->SetComputePerimeter(true);
i2l->Update();
//Create a map to contain all connectedComponent (even a little pixel)
typename I2LType::Pointer i2l = I2LType::New();
i2l->SetInput( connected->GetOutput() );
i2l->SetComputePerimeter(true);
i2l->Update();
//Write statitistics on the largest component
ShapeLabelObjectType *labelObject = labelMap->GetNthLabelObject(largestComponent);
std::cout << " Centroid: " << std::endl;
//Write statitistics on the largest component
ShapeLabelObjectType *labelObject = labelMap->GetNthLabelObject(largestComponent);
std::cout << " Centroid: " << std::endl;