|
|
| (One intermediate revision by one other user not shown) |
| Line 1: |
Line 1: |
| <div class="floatcenter">[[File:ITK_Examples_Baseline_ImageProcessing_TestRelabelComponentImageFilter.png]]</div>
| | {{warning|1=The media wiki content on this page is no longer maintained. The examples presented on the https://itk.org/Wiki/* pages likely require ITK version 4.13 or earlier releases. In many cases, the examples on this page no longer conform to the best practices for modern ITK versions.}} |
| ==Description==
| |
| Relabel the components (blobs) in an image such that consecutive labels are used. This example creates a black image with one rectangle consisting of pixels with value 100 and a second rectangle consisting of pixels with value 200. There are three resulting components (the background and the two rectangles). The relabel filter labels these three components 0, 1, and 2.
| |
| | |
| ==Displayed Output==
| |
| The original image and an intensity scaled version of the relabeled image. | |
| | |
| ==RelabelComponentImageFilter.cxx==
| |
| <source lang="cpp">
| |
| #include "itkImage.h"
| |
| #include "itkImageFileReader.h"
| |
| #include "itkCustomColormapFunction.h"
| |
| #include "itkScalarToRGBColormapImageFilter.h"
| |
| #include "itkRGBPixel.h"
| |
| #include "itkMersenneTwisterRandomVariateGenerator.h"
| |
| | |
| #include "itkRelabelComponentImageFilter.h"
| |
| | |
| #include "QuickView.h"
| |
| | |
| typedef itk::RGBPixel<unsigned char> RGBPixelType;
| |
| typedef itk::Image<RGBPixelType, 2> RGBImageType;
| |
| typedef itk::Image<unsigned char, 2> ImageType;
| |
| typedef itk::Function::CustomColormapFunction<
| |
| ImageType::PixelType, RGBImageType::PixelType> ColormapType;
| |
| | |
| static void CreateImage(
| |
| ImageType::Pointer image);
| |
| static void CreateRandomColormap(\
| |
| unsigned int size, ColormapType::Pointer colormap);
| |
| | |
| int main(int, char *[])
| |
| {
| |
| ImageType::Pointer image = ImageType::New();
| |
| | |
| CreateImage(image);
| |
|
| |
| typedef itk::RelabelComponentImageFilter<ImageType, ImageType> FilterType;
| |
| FilterType::Pointer relabelFilter =
| |
| FilterType::New();
| |
| relabelFilter->SetInput(image);
| |
| | |
| typedef itk::ScalarToRGBColormapImageFilter<ImageType, RGBImageType> ColormapFilterType;
| |
| ColormapFilterType::Pointer colormapFilter1 =
| |
| ColormapFilterType::New();
| |
| | |
| ColormapType::Pointer largeColormap
| |
| = ColormapType::New();
| |
| CreateRandomColormap(255, largeColormap);
| |
| | |
| colormapFilter1->SetInput (image);
| |
| colormapFilter1->SetColormap(largeColormap);
| |
| | |
| ColormapFilterType::Pointer colormapFilter2 =
| |
| ColormapFilterType::New();
| |
| colormapFilter2->SetInput (relabelFilter->GetOutput());
| |
| colormapFilter2->SetColormap(largeColormap);
| |
| | |
| QuickView viewer;
| |
| viewer.AddRGBImage(
| |
| colormapFilter1->GetOutput(),
| |
| true,
| |
| "Original");
| |
| | |
| viewer.AddRGBImage(
| |
| colormapFilter2->GetOutput(),
| |
| true,
| |
| "Relabeled");
| |
| | |
| viewer.Visualize();
| |
| | |
| return EXIT_SUCCESS;
| |
| }
| |
| | |
| void CreateImage(ImageType::Pointer image)
| |
| {
| |
| // Create an image with 2 connected components
| |
| ImageType::RegionType region;
| |
| ImageType::IndexType start;
| |
| start[0] = 0;
| |
| start[1] = 0;
| |
| | |
| ImageType::SizeType size;
| |
| size[0] = 200;
| |
| size[1] = 300;
| |
| | |
| region.SetSize(size);
| |
| region.SetIndex(start);
| |
| | |
| image->SetRegions(region);
| |
| image->Allocate();
| |
| | |
| itk::ImageRegionIterator<ImageType> imageIterator(image,region); | |
| | |
| while(!imageIterator.IsAtEnd())
| |
| {
| |
| if(imageIterator.GetIndex()[0] > 100 &&
| |
| imageIterator.GetIndex()[0] < 150 &&
| |
| imageIterator.GetIndex()[1] > 100 &&
| |
| imageIterator.GetIndex()[1] < 150)
| |
| {
| |
| imageIterator.Set(200);
| |
| }
| |
| else if(imageIterator.GetIndex()[0] > 50 &&
| |
| imageIterator.GetIndex()[0] < 70 &&
| |
| imageIterator.GetIndex()[1] > 50 &&
| |
| imageIterator.GetIndex()[1] < 70)
| |
| {
| |
| imageIterator.Set(100);
| |
| }
| |
| else
| |
| {
| |
| imageIterator.Set(0);
| |
| }
| |
| | |
| ++imageIterator;
| |
| }
| |
| }
| |
| | |
| void CreateRandomColormap(unsigned int size, ColormapType::Pointer colormap)
| |
| {
| |
| ColormapType::ChannelType redChannel;
| |
| ColormapType::ChannelType greenChannel;
| |
| ColormapType::ChannelType blueChannel;
| |
| itk::Statistics::MersenneTwisterRandomVariateGenerator::Pointer random =
| |
| itk::Statistics::MersenneTwisterRandomVariateGenerator::New();
| |
| | |
| random->SetSeed ( 8775070 );
| |
| for (unsigned int i = 0; i < size; ++i)
| |
| {
| |
| redChannel.push_back(static_cast<ColormapType::RealType>
| |
| (random->GetUniformVariate(.3, 1.0)));
| |
| greenChannel.push_back(static_cast<ColormapType::RealType>
| |
| (random->GetUniformVariate(.3, 1.0)));
| |
| blueChannel.push_back(static_cast<ColormapType::RealType>
| |
| (random->GetUniformVariate(.3, 1.0)));
| |
| }
| |
| colormap->SetRedChannel(redChannel);
| |
| colormap->SetGreenChannel(greenChannel);
| |
| colormap->SetBlueChannel(blueChannel);
| |
| }
| |
| </source>
| |
| | |
| {{ITKVTKCMakeLists|RelabelComponentImageFilter|}}
| |
