itkConnectedComponentImageFilter对输入图像大小的限制

来源：互联网发布：学霸君怎么样知乎编辑：程序博客网时间：2024/06/07 02:33

版本信息：Win10(64bit)，VS2013(64bit)，ITK4.8.2，VTK6.3.0，Qt5.4.2

1 问题描述

使用VTK实现对3D CT数据进行切片重新提取，然后使用ITK对切片进行二维分析，这是需要使用到itkConnectedComponentImageFilter将灰度图转换为Label Image。当切片的尺寸非常小时（如2x2），会出现错误：

 Debug Assertion Failed!  Program: C:\WINDOWS\SYSTEM32\MSVCP120D.dll File: C:\Program Files (x86)\Microsoft Visual Studio 12.0\VC\INCLUDE\vector Line: 1201  Expression: vector subscript out of range  For information on how your program can cause an assertion failure, see the Visual C++ documentation on asserts.  (Press Retry to debug the application)

上述提示vector下标越界，而比较令人疑惑的是，我并没有使用std::vector。继续使用VS2013进行调试发现错误原因为：

std::vector<std::vector<itk::ConnectedComponentImageFilter<itk::Image<unsigned char,2>,itk::Image<unsigned short,2>,itk::Image<unsigned char,2> >::runLength,std::allocator<itk::ConnectedComponentImageFilter<itk::Image<unsigned char,2>,itk::Image<unsigned short,2>,itk::Image<unsigned char,2> >::runLength> >,std::allocator<std::vector<itk::ConnectedComponentImageFilter<itk::Image<unsigned char,2>,itk::Image<unsigned short,2>,itk::Image<unsigned char,2> >::runLength,std::allocator<itk::ConnectedComponentImageFilter<itk::Image<unsigned char,2>,itk::Image<unsigned short,2>,itk::Image<unsigned char,2> >::runLength> > > >::operator[](unsigned __int64 _Pos) 行 1202C++

也即是说itkConnectedComponentImageFilter里面使用了std::vector，并且产生了下标越界。

2 解决方法

最好的方法就是不要使用itkConnectedComponentImageFilter，而改用其它方法生成Label Image：

//![0] -- UCharImageTypetypedef unsigned char UCharPixelType;typedef itk::Image<UCharPixelType, 2> ImageType;typedef ImageType itkImageData;typedef ImageType UCharImageType;// [0]//![1] -- UShortImageTypetypedef unsigned short UShortPixelType;typedef itk::Image<UShortPixelType, 2> UShortImageType;// [1]typedef itk::BinaryImageToLabelMapFilter<UCharImageType>    BinaryImageToLabelMapFilterType;typedef itk::LabelMapToLabelImageFilter        <BinaryImageToLabelMapFilterType::OutputImageType,        UShortImageType>    LabelMapToLabelImageFilterType;BinaryImageToLabelMapFilterType::Pointer binary2Label = BinaryImageToLabelMapFilterType::New();LabelMapToLabelImageFilterType::Pointer  labelMap2LabelImage = LabelMapToLabelImageFilterType::New();binary2Label->SetInput( image ); // image的类型是itkImageData，在上面有定义，此处略去imagebinary2Label->Update(); // 的产生方法labelMap2LabelImage->SetInput( binary2Label->GetOutput() );labelMap2LabelImage->UpdateLargestPossibleRegion();labelMap2LabelImage->Update();

0 0