vtk（5.0以上版本）5.3_ImageResliceExample2.cpp 最终成功运行版本

《VTK图形图像开发进阶》随书代码\Examples\Chap05  中的这个例子5.3_ImageResliceExample2.cpp有些错误，致使程序无法运行。

此篇博客经过多个网页查询，最终调试好的版本。填坑。

#include
VTK_MODULE_INIT(vtkRenderingOpenGL2);
VTK_MODULE_INIT(vtkInteractionStyle);
VTK_MODULE_INIT(vtkRenderingFreeType)

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include "vtkVolume16Reader.h"


class vtkImageInteractionCallback : public vtkCommand
{
public:
static vtkImageInteractionCallback *New()
{
return new vtkImageInteractionCallback;
}


vtkImageInteractionCallback()
{
this->Slicing = 0;
this->ImageReslice = 0;
this->Interactor = 0;
}


void SetImageReslice(vtkImageReslice *reslice)
{
this->ImageReslice = reslice;
}


void SetImageMapToColors(vtkImageMapToColors *mapToColors)
{
this->mapToColors = mapToColors;
}


vtkImageReslice *GetImageReslice()
{
return this->ImageReslice;
}


void SetInteractor(vtkRenderWindowInteractor *interactor)
{
this->Interactor = interactor;
}


vtkRenderWindowInteractor *GetInteractor()
{
return this->Interactor;
}


virtual void Execute(vtkObject *, unsigned long event, void *)
{
vtkRenderWindowInteractor *interactor = this->GetInteractor();


int lastPos[2];
interactor->GetLastEventPosition(lastPos);
int currPos[2];
interactor->GetEventPosition(currPos);


if (event == vtkCommand::LeftButtonPressEvent)
{
this->Slicing = 1;
}
else if (event == vtkCommand::LeftButtonReleaseEvent)
{
this->Slicing = 0;
}
else if (event == vtkCommand::MouseMoveEvent)
{
if (this->Slicing)
{
vtkImageReslice *reslice = this->ImageReslice;


// Increment slice position by deltaY of mouse
int deltaY = lastPos[1] - currPos[1];


reslice->Update();
double sliceSpacing = reslice->GetOutput()->GetSpacing()[2];
vtkMatrix4x4 *matrix = reslice->GetResliceAxes();
// move the center point that we are slicing through
double point[4];
double center[4];
point[0] = 0.0;
point[1] = 0.0;
point[2] = sliceSpacing * deltaY;
point[3] = 1.0;
matrix->MultiplyPoint(point, center);
matrix->SetElement(0, 3, center[0]);
matrix->SetElement(1, 3, center[1]);
matrix->SetElement(2, 3, center[2]);
mapToColors->Update();
interactor->Render();
}
else
{
vtkInteractorStyle *style = vtkInteractorStyle::SafeDownCast(
interactor->GetInteractorStyle());
if (style)
{
style->OnMouseMove();
}
}
}
}


private:
int Slicing;
vtkImageReslice *ImageReslice;
vtkRenderWindowInteractor *Interactor;
vtkImageMapToColors *mapToColors;
};


int main()
{
vtkSmartPointer reader =
vtkSmartPointer::New();
reader->SetFileName("~\\data\\brain.mhd");
reader->Update();

int extent[6];
double spacing[3];
double origin[3];

reader->GetOutput()->GetExtent(extent);
reader->GetOutput()->GetSpacing(spacing);
reader->GetOutput()->GetOrigin(origin);
double center[3];
center[0] = origin[0] + spacing[0] * 0.5 * (extent[0] + extent[1]);
center[1] = origin[1] + spacing[1] * 0.5 * (extent[2] + extent[3]);
center[2] = origin[2] + spacing[2] * 0.5 * (extent[4] + extent[5]);


static double axialElements[16] = {
1, 0, 0, 0,
0, 1, 0, 0,
0, 0, 1, 0,
0, 0, 0, 1
};


vtkSmartPointer resliceAxes =
vtkSmartPointer::New();
resliceAxes->DeepCopy(axialElements);


resliceAxes->SetElement(0, 3, center[0]);
resliceAxes->SetElement(1, 3, center[1]);
resliceAxes->SetElement(2, 3, center[2]);


vtkSmartPointer reslice =
vtkSmartPointer::New();
reslice->SetInputData(reader->GetOutput());//这里与vtk5以前的版本不同，也是出错的主要原因
reslice->SetOutputDimensionality(2);
reslice->SetResliceAxes(resliceAxes);
reslice->SetInterpolationModeToLinear();
reslice->Update();//这里不要漏了


vtkSmartPointer colorTable =
vtkSmartPointer::New();
colorTable->SetRange(0, 1000);
colorTable->SetValueRange(0.0, 1.0);
colorTable->SetSaturationRange(0.0, 0.0);
colorTable->SetRampToLinear();
colorTable->Build();


vtkSmartPointer colorMap =
vtkSmartPointer::New();
colorMap->SetLookupTable(colorTable);
colorMap->SetInputConnection(reslice->GetOutputPort());
colorMap->Update();//这里不要漏了


vtkSmartPointer imgActor =
vtkSmartPointer::New();
imgActor->SetInputData(colorMap->GetOutput());


vtkSmartPointer renderer =
vtkSmartPointer::New();
renderer->AddActor(imgActor);
renderer->SetBackground(1, 1, 1);


vtkSmartPointer renderWindow =
vtkSmartPointer::New();
renderWindow->SetSize(500, 500);
renderWindow->AddRenderer(renderer);


vtkSmartPointer renderWindowInteractor =
vtkSmartPointer::New();
vtkSmartPointer imagestyle =
vtkSmartPointer::New();


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


vtkSmartPointer callback =
vtkSmartPointer::New();
callback->SetImageReslice(reslice);
callback->SetInteractor(renderWindowInteractor);
callback->SetImageMapToColors(colorMap);//这里是新加的部分，参考其他的网页


imagestyle->AddObserver(vtkCommand::MouseMoveEvent, callback);
imagestyle->AddObserver(vtkCommand::LeftButtonPressEvent, callback);
imagestyle->AddObserver(vtkCommand::LeftButtonReleaseEvent, callback);


renderWindowInteractor->Start();
return EXIT_SUCCESS;
}