ITK image获取像素值的多种方法

对于图像处理中，对像素进行遍历的处理是基本操作，所以在使用第三方库的时候，如ITK，需要知道如何获取像素的值，才能更好的进行扩展，如添加自定义的滤波等；另一方面，方便将ITK的数据类转换为自己软件定义的数据封装类或者转换为其他库的数据类。

 

根据使用过程中，我总结出4种方法如下：

• 调用GetPixel 函数(低效率)
• 采用迭代器(高效率)
• 采用重载()操作符(低效率)
• 直接提取内存buffer(高效率)

 

调用GetPixel 函数(低效率)

    using ImageType = itk::Image<int, 3>;
    using ReaderType = itk::ImageFileReader<ImageType>;
    const char* src_path = "H:\\tmp\\a.mhd";
    auto reader = ReaderType::New();
    reader->SetFileName(src_path);
    reader->Update();

    auto image = reader->GetOutput();

    ImageType::IndexType idx{92,207,10};
    std::cout<<image->GetPixel(idx)<<std::endl;

采用迭代器(高效率)

    using ImageType = itk::Image<int, 3>;
    using ReaderType = itk::ImageFileReader<ImageType>;
    const char* src_path = "H:\\tmp\\a.mhd";
    auto reader = ReaderType::New();
    reader->SetFileName(src_path);
    reader->Update();

    auto image = reader->GetOutput();
    
    using IteratorType = itk::ImageRegionIterator<ImageType>;
    IteratorType it(image, image->GetLargestPossibleRegion());
    it.GoToBegin();
    while(!it.IsAtEnd()){
        std::cout<<it.Get()<<std::endl;
        ++it;
    }

采用重载()操作符(低效率)

    using ImageType = itk::Image<int, 3>;
    using ReaderType = itk::ImageFileReader<ImageType>;
    const char* src_path = "H:\\tmp\\a.mhd";
    auto reader = ReaderType::New();
    reader->SetFileName(src_path);
    reader->Update();

    auto image = reader->GetOutput();
    
    ImageType::IndexType idx{92,207,10};
    std::cout<<(*image)[idx]<<std::endl;

直接提取内存buffer(高效率)

	using ImageType = itk::Image<int, 3>;
    using ReaderType = itk::ImageFileReader<ImageType>;
    const char* src_path = "H:\\tmp\\a.mhd";
    auto reader = ReaderType::New();
    reader->SetFileName(src_path);
    reader->Update();

    auto image = reader->GetOutput();

    int * buffer = image->GetBufferPointer(); // 这是一块连续的内存
    ImageType::SizeType size = image->GetLargestPossibleRegion().GetSize();
    std::cout<<buffer[10 * size[0] * size[1] + 207 * size[0] + 92]<<std::endl;