Как подключить ITK к VTK с C++?

Question

Я новичок в ITK, VTK и Qt. Я использую Visual Studio 9. Я пытаюсь прочитать серию DICOM с помощью ITK и отобразить с помощью VTK в QVTKWidget (Qt). Я основывался на этом коде http://www.itk.org/Wiki/VTK/Examples/Cxx/IO/ReadDICOMSeries и модифицировал в соответствии со своими потребностями.

когда я читаю серию DICOM с помощью VTK и отображаю ее в QVTKWidget, она работает, но когда я хочу прочитать эту серию с помощью ITK и отображать с помощью VTK в QVTKWidget, программа отображает первое изображение серии, а когда я перехожу к следующему изображению с колесо мыши, программа вылетает.

когда я отладил, я получил эту ошибку:

void VTKImageExportBase::UpdateInformationCallbackFunction(void *userData)
{
static_cast< VTKImageExportBase * >
( userData )->UpdateInformationCallback();//the error is here
}

Я пробовал ausssi ImageToVTKImageFilter класс, но с той же проблемой.

вот мой код:

void essaieAppQtVTK::drawDCMSeries(std::string folderDCM)
{
typedef unsigned short    PixelType;
const unsigned int      Dimension = 3;

typedef itk::Image< PixelType, Dimension >         ImageType;

typedef itk::VTKImageExport<ImageType> ImageExportType; 
typedef itk::ImageSeriesReader< ImageType >        ReaderType;
ReaderType::Pointer reader = ReaderType::New();
typedef itk::GDCMImageIO       ImageIOType;
ImageIOType::Pointer dicomIO = ImageIOType::New(); 

reader->SetImageIO( dicomIO );
typedef itk::GDCMSeriesFileNames NamesGeneratorType;
NamesGeneratorType::Pointer nameGenerator = NamesGeneratorType::New();

nameGenerator->SetUseSeriesDetails( true );
nameGenerator->AddSeriesRestriction("0008|0021" );  
nameGenerator->SetDirectory( folderDCM);

typedef std::vector< std::string >    SeriesIdContainer;    
const SeriesIdContainer & seriesUID = nameGenerator->GetSeriesUIDs();   
std::cout << seriesUID.size() << std::endl;
SeriesIdContainer::const_iterator seriesItr = seriesUID.begin();
SeriesIdContainer::const_iterator seriesEnd = seriesUID.end();
while( seriesItr != seriesEnd )
{
std::cout << seriesItr->c_str() << std::endl;
seriesItr++;
}
std::string seriesIdentifier;
seriesIdentifier = seriesUID.begin()->c_str();
std::cout << seriesIdentifier.c_str() << std::endl;

typedef std::vector< std::string >   FileNamesContainer;
FileNamesContainer fileNames; 
fileNames = nameGenerator->GetFileNames( seriesIdentifier ); 

reader->SetFileNames( fileNames );
try
{
reader->Update();
}
catch (itk::ExceptionObject &ex)
{
std::cout << ex << std::endl;
}

//------------------------------------------------------------------------
// ITK to VTK pipeline connection.
//------------------------------------------------------------------------

// Create the itk::VTKImageExport instance and connect it to the
// itk::CurvatureFlowImageFilter.
ImageExportType::Pointer exporter = ImageExportType::New();
exporter->SetInput(reader->GetOutput());

// Create the vtkImageImport and connect it to the
// itk::VTKImageExport instance.
vtkImageImport* importer = vtkImageImport::New();  
ConnectPipelines(exporter, importer);

//------------------------------------------------------------------------
// VTK pipeline.
//------------------------------------------------------------------------

this->imageViewer= vtkImageViewer2::New();
imageViewer->SetInput(importer->GetOutput());


// slice status message
//******same code *****//

// usage hint message
//******same code *****//

// create an interactor with our own style (inherit from vtkInteractorStyleImage)
// in order to catch mousewheel and key events
//******same code *****//

// add slice status message and usage hint message to the renderer
//******same code *****//

//to display the result: 
ui.qvtkWidget->SetRenderWindow(imageViewer->GetRenderWindow()); 
ui.qvtkWidget->GetRenderWindow()->GetInteractor()-   >SetInteractorStyle(myInteractorStyle); 
imageViewer->Render(); 
ui.qvtkWidget->update();
}

Необязательно: мой экспортер и импортер указаны ниже:

template <typename ITK_Exporter, typename VTK_Importer>
void ConnectPipelines(ITK_Exporter exporter, VTK_Importer* importer)
{
importer->SetUpdateInformationCallback(exporter->GetUpdateInformationCallback());
importer->SetPipelineModifiedCallback(exporter->GetPipelineModifiedCallback());
importer->SetWholeExtentCallback(exporter->GetWholeExtentCallback());
importer->SetSpacingCallback(exporter->GetSpacingCallback());
importer->SetOriginCallback(exporter->GetOriginCallback());
importer->SetScalarTypeCallback(exporter->GetScalarTypeCallback());
importer->SetNumberOfComponentsCallback(exporter->GetNumberOfComponentsCallback());
importer->SetPropagateUpdateExtentCallback(exporter- >GetPropagateUpdateExtentCallback());
importer->SetUpdateDataCallback(exporter->GetUpdateDataCallback());
importer->SetDataExtentCallback(exporter->GetDataExtentCallback());
importer->SetBufferPointerCallback(exporter->GetBufferPointerCallback());
importer->SetCallbackUserData(exporter->GetCallbackUserData());
}
/**
* This function will connect the given vtkImageExport filter to
* the given itk::VTKImageImport filter.
*/
template <typename VTK_Exporter, typename ITK_Importer>
void ConnectPipelines(VTK_Exporter* exporter, ITK_Importer importer)
{
importer->SetUpdateInformationCallback(exporter->GetUpdateInformationCallback());
importer->SetPipelineModifiedCallback(exporter->GetPipelineModifiedCallback());
importer->SetWholeExtentCallback(exporter->GetWholeExtentCallback());
importer->SetSpacingCallback(exporter->GetSpacingCallback());
importer->SetOriginCallback(exporter->GetOriginCallback());
importer->SetScalarTypeCallback(exporter->GetScalarTypeCallback());
importer->SetNumberOfComponentsCallback(exporter->GetNumberOfComponentsCallback());
importer->SetPropagateUpdateExtentCallback(exporter->GetPropagateUpdateExtentCallback());
importer->SetUpdateDataCallback(exporter->GetUpdateDataCallback());
importer->SetDataExtentCallback(exporter->GetDataExtentCallback());
importer->SetBufferPointerCallback(exporter->GetBufferPointerCallback());
importer->SetCallbackUserData(exporter->GetCallbackUserData());
}

Возможно, в трубопроводе возникла ошибка между ITK и VTK. Пожалуйста, помогите мне найти решение этой проблемы. Я потратил две недели на поиск решения, но все методы не сработали, возможно, есть другой способ связать ITK и VTK. кроме класса ImageToVTKImageFilter. Я рассчитываю на вашу помощь. заранее спасибо.

Answer 1

Я столкнулся с той же ошибкой в ​​тот же момент, что вы показали в пустоте VTKImageExportBase::UpdateInformationCallbackFunction функция, которую я использовал ImageToVTKImageFilter для конвертации. Чего мне не хватало, так это обновления конвейера ITK перед его подключением к средству просмотра VTK. В ITK to VTK pipeline connection часть вашего кода также отсутствует то же самое. Так что используйте:

exporter->SetInput(reader->GetOutput());
exporter->Update();

тогда иди дальше.