ITK分割，VTK+VC怎么显示

2009-03-18 07:09:44 411 阅读 AI 写同款 · GEO 优化AI 选题数据监测一键分发 ›

31套VTK3D图像体绘制/VTK光线投射法/VTK三维重建程序源码 31套VTK3D图像体绘制/VTK光线投射法/VTK三维重建程序源码 1、基于VTK的3D图像体绘制（实现了以VTK为基础的3D图像重建） 2、基于VTK的人头骨3D图像（用表面重建的方法。vs2005，程序比较简单） 3、读取并可视化D1O0M图像（VIK &OpenCV）（将DIOXM格式的CT图用VTK读取，转给OpenCV处理并可视化） 4、VTK三维可视化读取RAM数据的c++源代码VS2010 5、VTK8.1VS2017编译好的VTK库c++（vtk8.1编译好的库R/D版本，使用vs 阅读详情

请问，有谁知道ITK分割后的图像通过VTK和VC显示出来？帮帮忙谢谢！我写了一个代码，编译能通过，但是老是在Viewer->Render()出错。
CVTK::CVTK()
{
lowerThreshold = 150.0;
upperThreshold = 180.0;
index[0] =60;
index[1] =116;
ReaderType1::Pointer reader = ReaderType1::New();
reader->SetFileName("G:\\ITKVTK\\test\\test\\myProject\\Debug\\BrainProtonDensitySlice.png");
WriterType1::Pointer writer = WriterType1::New();
writer->SetFileName("G:\\ITKVTK\\test\\test\\myProject\\Debug\\aa.png");
CurvatureFlowImageFilterType::Pointer smoothing = CurvatureFlowImageFilterType::New();
ConnectedFilterType::Pointer connectedThreshold = ConnectedFilterType::New();
CastingFilterType::Pointer caster = CastingFilterType::New();
smoothing->SetInput( reader->GetOutput() );
connectedThreshold->SetInput( smoothing->GetOutput() );
caster->SetInput( connectedThreshold->GetOutput() );
writer->SetInput( caster->GetOutput() );

smoothing->SetNumberOfIterations( 5 );
smoothing->SetTimeStep( 0.125 );
connectedThreshold->SetLower( lowerThreshold );
connectedThreshold->SetUpper( upperThreshold );
connectedThreshold->SetReplaceValue( 255 );
connectedThreshold->SetSeed( index );
ConnectorType::Pointer connector = ConnectorType::New();
connector->SetInput( caster->GetOutput());

// writer->Update();

viewer = vtkImageViewer::New();
viewer->SetInput(connector->GetOutput());
viewer->SetColorWindow(256);
viewer->SetColorLevel(128);

iren = vtkWin32RenderWindowInteractor::New();
viewer->SetupInteractor(iren);
// connector->Update();
picker=vtkCellPicker::New();
//picker->SetTolerance(0.00001);

vtkCallbackCommand * cmd = vtkCallbackCommand::New();
cmd->SetCallback(CVTK::MyExecute);
cmd->SetClientData( this );
picker->AddObserver(vtkCommand::EndPickEvent, cmd);
iren->SetPicker(picker);
cmd->Delete();

}

void CVTK::BeginRenderOn(CStatic * aStatic)
{
CRect rect;
aStatic->GetClientRect(&rect);
viewer->SetSize(rect.Width(),rect.Height());
viewer->SetParentId(aStatic->m_hWnd); //关键代码一：设置父窗口
viewer->Render(); // 关键代码二：开始绘制，启动交互器
// iren->Initialize();
}
void CVTK::MyExecute(vtkObject *caller, unsigned long eid, void *clientdata, void *calldata)
{
CVTK *m_viewer = reinterpret_cast<CVTK*>(clientdata );
if (m_viewer->picker->GetCellId()<0)
{
AfxMessageBox("CellID < 0");
}
else
{
AfxMessageBox("CellID > 0");
double ScreenCoordinate[3],globalCoordinate[3]; //拾取点的屏幕坐标和全局坐标
m_viewer->picker->GetSelectionPoint(ScreenCoordinate); //得到屏幕坐标
m_viewer->picker->GetPickPosition(globalCoordinate); //得到全局坐标
};
}

...全文