-1

尝试向具有浮点像素类型的 ikt::Image显示或应用任何过滤器时出现分段错误。我可以使用相同的代码（交换像素类型）打开一个具有unsigned char 像素类型的图像，它工作得很好。我提到的两个图像（浮点像素图像和无符号字符图像）都是 Nifti 类型。

我还使用 ImageIOBase 打开了图像以检查像素类型及其确实是浮动的。

#include "itkImage.h"
#include "itkImageFileReader.h"
#include "itkExtractImageFilter.h"
#include "itkImageToVTKImageFilter.h"
#include "vtkImageViewer.h"
#include "vtkImageMapper3D.h"
#include "vtkRenderWindowInteractor.h"
#include "vtkSmartPointer.h"
#include "vtkImageActor.h"
#include "vtkInteractorStyleImage.h"
#include "vtkRenderer.h"
#include "itkRGBPixel.h"
#include <unistd.h>
#include "QuickView.h"
#include <armadillo>
#include "itkBinaryThresholdImageFilter.h"
#include "itkMedianImageFilter.h"



typedef float FloatPixelType;
typedef itk::Image<FloatPixelType,  3> Input3DFloatImageType;
typedef itk::Image<FloatPixelType, 2> FloatImageType;
typedef itk::ImageFileReader< Input3DFloatImageType> Float3DReaderType;
typedef itk::ImageToVTKImageFilter<FloatImageType> FloatConnectorType;
typedef itk::ExtractImageFilter< Input3DFloatImageType, FloatImageType > FloatFilterType;

int main(int argc, char *argv[]){

const char * grayMatterFile = "../../I3T/I3TGM.img";
Float3DReaderType::Pointer floatReader = Float3DReaderType::New();
floatReader->SetFileName( grayMatterFile );

FloatFilterType::Pointer floatFilter = FloatFilterType::New();
floatFilter->SetDirectionCollapseToSubmatrix();
floatFilter->InPlaceOn();
floatReader->UpdateOutputInformation();
Input3DFloatImageType::RegionType floatRegion = floatReader->GetOutput()->GetLargestPossibleRegion();

Input3DFloatImageType::SizeType floatSize = floatRegion.GetSize();
floatSize[1] = 0;
Input3DFloatImageType::IndexType floatIndexStart = floatRegion.GetIndex();

for(int i = 0; floatRegion.GetSize()[1]; i++){
        
        floatIndexStart[1] = i;
        Input3DFloatImageType::RegionType floatDesiredRegion;
        
        floatDesiredRegion.SetSize( floatSize);
        floatDesiredRegion.SetIndex( floatIndexStart);
    
        // Float image
        floatFilter->SetInput( floatReader->GetOutput() );
        floatFilter->SetExtractionRegion( floatDesiredRegion ); 
        floatFilter->Update();
        ShowImage<FloatImageType>(floatFilter->GetOutput());
    
    
    }
}

template<typename TImageType>
void ShowImage(typename TImageType::Pointer image){

    typedef itk::ImageToVTKImageFilter<TImageType> ThisConnectorType;
    typename ThisConnectorType::Pointer connector = ThisConnectorType::New();
    connector->SetInput(image);

    try{
        connector->Update();
     } catch( itk::ExceptionObject & error ) {
        std::cerr << "Error: " << error << std::endl;
        return;
     }


    vtkSmartPointer<vtkImageActor> actor =  vtkSmartPointer<vtkImageActor>::New();

    actor->GetMapper()->SetInputData(connector->GetOutput());

    vtkSmartPointer<vtkRenderer> renderer = vtkSmartPointer<vtkRenderer>::New();
    renderer->AddActor(actor);
    renderer->ResetCamera();

    vtkSmartPointer<vtkRenderWindow> renderWindow = vtkSmartPointer<vtkRenderWindow>::New();
    renderWindow->AddRenderer(renderer);

    vtkSmartPointer<vtkRenderWindowInteractor> renderWindowInteractor = vtkSmartPointer<vtkRenderWindowInteractor>::New();
    vtkSmartPointer<vtkInteractorStyleImage> style = vtkSmartPointer<vtkInteractorStyleImage>::New();

    renderWindowInteractor->SetInteractorStyle(style);

    renderWindowInteractor->SetRenderWindow(renderWindow);
    renderWindowInteractor->Initialize();

    renderWindowInteractor->Start();
}

行中出现分段错误：floatFilter->Update();

4

1 回答 1

1

你缺乏floatFilter->Update();对前ShowImage。

编辑：

您更新的代码（稍作修改）在我的计算机上运行。您安装的 ITK 或 VTK 可能已损坏。尝试使用 IDE 在调试模式下运行它，这将显示故障发生在哪里。

于 2017-11-30T14:20:00.000 回答