0

我正在尝试使用 VTK 和 C++ 编写一个程序来读取 MRI 数据。但我无法获得主要的 MRI 原始数据的间距。“GetSpacing”仅适用于“ReadImageData”函数。我想我在 C++ 编程中犯了一些错误。但我不知道它在哪里。

vtkImageData* ReadImageData(string mri_imagedata_file)
{
    vtkSmartPointer<vtkMetaImageReader> reader =
        vtkSmartPointer<vtkMetaImageReader>::New();
    reader->SetFileName(mri_imagedata_file.c_str());
    reader->Update();
    vtkImageData* metaimage = reader->GetOutput();
    double sp[3];
    metaimage->GetSpacing(sp);
    cout << sp[0] << " " << sp[1] << " " << sp[2] <<endl; //<----------It works here.
    return metaimage;
}

int main (int argc, char *argv[])
{
    if(argc != 2)
    {
        cerr << "Usage: " << argv[0] << " MRI image data" <<endl;
        return EXIT_FAILURE;
    }
    string mri_imagedata_file = argv[1];// Input "prost00.mhd"

    vtkImageData* metaimage = ReadImageData(mri_imagedata_file);
    double sp2[3];
    metaimage->GetSpacing(sp2);
    cout << sp2[0] << " " << sp2[1] << " " << sp2[2] << endl; //<-----It doesn't work here

}

感谢您的关注。

4

3 回答 3

0

我的假设是vtkMetaImageReader::GetOutput()返回指向一些vtkMetaImageReader内部数据的指针，所以当你退出时你ReadImageData的reader被破坏并且你返回的指针变得无效。

于 2013-04-23T11:08:44.970 回答
0

假设vtkSmartPointer<vtkMetaImageReader>是一种智能指针， 函数中的指针在返回reader 时被破坏。ReadImageData这包括metaimage指向的内容。

    return metaimage;
    // reader->~();  // Destructed here. Including what is pointed to by metaimage. 
}

一个好的解决方案是返回智能指针而不是metaimage. 这样，返回reader时指向的内容不会被破坏，ReadImageData并且在 main 中可用。

于 2013-04-23T11:10:02.407 回答
0

看起来您忘记将数组作为参数传递给GetSpacingin main，因此它调用了返回 a 的重载double*，而您的数组保持不变。

您还打印了一个名为的数组，spacing尽管看起来您希望在sp2.

于 2013-04-23T11:10:25.463 回答