《itk实用demo》-彩色连通域分块

最新推荐文章于 2021-12-28 16:16:50 发布

原创最新推荐文章于 2021-12-28 16:16:50 发布 · 599 阅读

2 ·

CC 4.0 BY-SA版权

ITK医学图像处理专栏收录该内容

94 篇文章

订阅专栏

本文介绍了一种使用ITK库标记图像中不同连通域的方法，并通过颜色区分这些区域。该程序读取灰度图像，利用连通组件过滤器识别独立对象，再将结果转换为彩色图像并保存。

将图像中不同的连通域用颜色标记出来。

#include "itkImage.h"
#include "itkImageFileReader.h"
#include "itkConnectedComponentImageFilter.h"
#include "itkLabelToRGBImageFilter.h"

#include "itksys/SystemTools.hxx"
#include <sstream>

#include "itkImageToVTKImageFilter.h"
#include "vtkImageViewer.h"
#include "vtkWin32RenderWindowInteractor.h"

#include "itkImageFileWriter.h"
#include <iostream>
#include <string>
using namespace std;
int main( int argc, char *argv[])
{
 const unsigned int Dimension = 2;
 typedef unsigned char                       PixelType;
 typedef itk::RGBPixel<unsigned char>         RGBPixelType;
 typedef itk::Image<PixelType, Dimension>     ImageType;
 typedef itk::Image<RGBPixelType, Dimension>  RGBImageType;
 std::string bmp = ".bmp";
 ImageType::Pointer image;
 typedef itk::ImageFileReader<ImageType> ReaderType;

 for (int num= 0;num<214;num++)
 {
  std::string  inputFilename = "D:\\se6\\bmp\\IM";
  std::string  outputFilename = "D:\\se6\\rgb\\IM";

  char t[256];
  string s;
  sprintf(t, "%d", num);
  s = t;
  inputFilename = inputFilename + s + bmp;
  outputFilename = outputFilename + s + bmp;

  ReaderType::Pointer reader = ReaderType::New();
  reader->SetFileName(inputFilename);
  reader->Update();

  image = reader->GetOutput();

  typedef itk::Image< unsigned short, Dimension > OutputImageType;

  typedef itk::ConnectedComponentImageFilter <ImageType, OutputImageType >
   ConnectedComponentImageFilterType;

  ConnectedComponentImageFilterType::Pointer connected =
   ConnectedComponentImageFilterType::New ();
  connected->SetInput(image);
  connected->Update();

  std::cout << "Number of objects: " << connected->GetObjectCount() << std::endl;

  typedef itk::LabelToRGBImageFilter<OutputImageType, RGBImageType> RGBFilterType;
  RGBFilterType::Pointer rgbFilter = RGBFilterType::New();
  rgbFilter->SetInput( connected->GetOutput() );

  typedef  itk::ImageFileWriter< RGBImageType  > WriterType;
  WriterType::Pointer writer = WriterType::New();
  writer->SetFileName(outputFilename);
  writer->SetInput(rgbFilter->GetOutput());
  writer->Update();
 }


 return EXIT_SUCCESS;
}