c++ - 使用 ITK 类的 C++ 上的奇怪段错误

Question

自 1 或 2 周前以来，我一直在尝试解决这个非常奇怪的问题。

我正在尝试使用 ITK 读取 DICOM 图像，然后进行一些处理。

问题是，我想独立于像素类型读取图像。为此，我必须找出像素类型并在运行时实例化 itkImage。

因此，当我读取图像文件时，我必须将它从 itkImage 转换为 itkImageBase(在 itkImage.h 中它被声明为 class itkImage : public itkImageBase)。

不幸的是，当我在 if block 内将它从一种类型转换为另一种类型后，我无法在该 if block 之外的任何地方访问存储在该变量上的值。这是代码:

*抱歉葡萄牙语的变量和方法名称。这里有一点帮助:

leitor 表示读者

tipoPixel 表示像素类型

abreImagem(templated method)表示openImage

iniciaAbrirImagem 表示 initiateOpenImage

//here is io.h
#ifndef IO_H
#define IO_H

#include "itkImageSeriesReader.h"
#include "itkImageFileWriter.h"
#include "itkGDCMSeriesFileNames.h"
#include "itkGDCMImageIO.h"
#include "itkMetaDataDictionary.h"
#include "itkImageBase.h"
#include "iobase.h"

typedef itk::GDCMImageIO GDCM;
typedef itk::Image<float, 3> TipoImagemIO;
typedef itk::MetaDataDictionary TipoDicionario;
typedef itk::MetaDataObject<std::string> MetaDataString;
typedef itk::ImageBase<3> ImageBase; //Typedef for the imagebase

class IO : public IOBase
{
private:

    ImageBase* imagemB; //This is the variable giving me some trouble
    ImageBase* imagemB2;

    GDCM::Pointer gdcm;

    std::string tipoPixel;
public:
    IO();
    template<typename TLeitor>void abreImagem(std::string s, TLeitor leitor);
    void abreDiretorio(std::string s);
    void iniciaGravarImagem(std::string s, ImageBase* imgNova);
    ImageBase* GetOutput();
    void extraiTags();
    void iniciaAbrirImagem(std::string s);
    template<typename TGravador, typename TImagem>void gravaImagem(TGravador gravador, TImagem img, std::string s);
};
#endif  


//here is io.cpp
#include "io.h"
#include "vtkKWImage.h"
#include "vtkKWImageIO.h"
#include "imagem.h"
#include <typeinfo>

IO::IO()
{
imagemB = 0;
imagemB2 = 0;
gdcm = GDCM::New();
}

template<typename TLeitor>
void IO::abreImagem(std::string s, TLeitor leitor)
{
leitor->SetImageIO(gdcm);
leitor->SetFileName(s);

try
{
    leitor->Update();
}
catch(itk::ExceptionObject &ex)
{
    std::cerr<<ex<<std::endl;
    exit(1);
}

try
{
    imagemB = static_cast<ImageBase*>(leitor->GetOutput());
}
catch(std::bad_cast &ex)
{
    std::cout<<"EXCECAO";
    exit(1);
}
std::cout<<imagemB->GetNameOfClass(); //This one is ok
}

void IO::iniciaAbrirImagem(std::string s)
{
extraiTipoDimensoes(s); 
tipoPixel = GetTipo();  

if(tipoPixel.compare("short") == 0)
{
    typedef itk::Image<short, 3> itkImg;
    typedef itk::ImageSeriesReader<itkImg> TipoLeitor;
    TipoLeitor::Pointer leitor = TipoLeitor::New();

    abreImagem<TipoLeitor::Pointer>(s, leitor);

std::cerr<<"Trying to call a method inside IF: "<<imagemB->GetNameOfClass()<<std::endl;  //This one is ok too
}
//...doing the same things for unsigned short, int, float, double, etc...
std::cerr<<"trying to call a method outside IF: ";  
std::cerr<<imagemB->GetNameOfClass()<<std::endl; //Here i get this weird segmentation fault
}

ImageBase* IO::GetOutput()
{
std::cerr<<"JUST ENTERED GETOUTPUT()";  
std::cout<<imagemB->GetNameOfClass()<<std::endl; //If i remove the last method call of this variable, i get a segfault here
std::cerr<<std::endl<<"HOORAY, ITS GETTING OUT OF GETOUPUT"<<std::endl;
return imagemB;
}

如您所见(或者至少我希望代码足够清晰以便您能看到)，我声明了一个指向 itk::ImageBase<3> 的指针并在类构造函数上用值 0 启动它(已经尝试过没有初始化，用 null 初始化，没有任何作用...)。

然后，我调用其他类的一些方法来找出输入图像的像素类型。然后，我进入并在 if block 中实例化一个 itk::ImageFileReader 和一个与输入图像具有相同像素类型的 itkImage。

然后，我将类型为 TipoLeitor(或 ReaderType，如您所愿)的变量发送到模板函数，我实际上在 try/catch block 中读取图像文件。然后，在下一次尝试中，我将读取器的输出(一个 itkImage)转换为 itkImageBase，并将其放入 imagemB 变量中。

如果我尝试调用该类的一些随机方法，例如模板化方法中的 GetNameOfClass，没关系，一切都按预期进行。如果我在调用模板化方法的 if block 内做同样的事情，那也没关系。但是，如果我在 if block 之后尝试在函数内部做同样的事情，它就不起作用了。如果我尝试在 GetOutput() 方法上执行此操作，也会发生同样的情况。

这就像变量 imagemB 在 if block 中初始化，当它退出 block 时它被设置为未初始化(这个词是否存在)。

没有详细的错误信息，只有“Segmentation fault”。

感谢您的帮助，对于糟糕的英语和冗长的文本感到抱歉。

Answer 1

当你说“独立于像素类型读取图像”时，如果你是在谈论磁盘上图像的像素类型，你不需要通过它: ImageSeriesReader 将在内部从磁盘上的像素类型转换为您作为 ImageSeriesReader 模板放置的图像的像素类型。 请注意，ImageSeriesReader 的模板是output 中的图像类型，而不是input。

如果你真的想在运行时控制输出类型，你可以使用模板函数:

template<typename OutputPixelType> void my_process() {
}

并将其用于:

switch (type)
  {
  case 1:
    my_process<unsigned char>(parseResult);
    break;
  case 2:
    my_process<short int>(parseResult);
    break;
  ...
  }

这里有完整的例子 http://hg.orfeo-toolbox.org/OTB-Applications/file/1c193d45e483/Utils/otbConvert.cxx