网上找了一些关于基于VTK的任意平面切割,但是代码都是那一个一样的基于OpenCV
的代码,而假如我不想使用OpenCV,(而且我还找不到头文件"VTKReBuild.h")那么这个程序改怎么写呢?而且我想切割出来的是一个有一定厚度的立体图,相当于一个圆环类似(以人体模型为例,若选择从腹部切割开来的话,保留的便是以切割平面为基准大概1cm上下的切割图)恳求大佬指点!
https://blog.csdn.net/u014003644/article/details/76559800
十分感谢你的回答,但是我想请问下头文件vtkTestUtilities.h和vtkRegressionTestImage.h是哪来的呢?请问是在安装VTK的时候我漏装了环境而使我没有这两个头文件吗?
#include "vtkAutoInit.h"
VTK_MODULE_INIT(vtkRenderingOpenGL2);
VTK_MODULE_INIT(vtkInteractionStyle);
#include "vtkSmartPointer.h"
#include "vtkRegressionTestImage.h"
#include "vtkActor.h"
#include "vtkCamera.h"
#include "vtkCellPicker.h"
#include "vtkCommand.h"
#include "vtkImageActor.h"
#include "vtkImageReslice.h"
#include "vtkInteractorStyleImage.h"
#include "vtkImageMapToColors.h"
#include "vtkImagePlaneWidget.h"
#include "vtkImageReader.h"
#include "vtkInteractorEventRecorder.h"
#include "vtkLookupTable.h"
#include "vtkOutlineFilter.h"
#include "vtkDICOMImageReader.h"
#include "vtkPolyDataMapper.h"
#include "vtkProperty.h"
#include "vtkRenderWindow.h"
#include "vtkRenderWindowInteractor.h"
#include "vtkRenderer.h"
#include "vtkVolume16Reader.h"
#include "vtkImageData.h"
#include "vtkPointData.h"
#include "vtkPlaneSource.h"
#include "vtkPlane.h"
#include "vtkResliceCursorActor.h"
#include "vtkResliceCursorPolyDataAlgorithm.h"
#include "vtkResliceCursor.h"
#include "vtkResliceCursorWidget.h"
#include "vtkResliceCursorLineRepresentation.h"
#include "vtkBiDimensionalWidget.h"
#include "vtkTestUtilities.h"
//----------------------------------------------------------------------------
class vtkResliceCursorCallback : public vtkCommand
{
public:
static vtkResliceCursorCallback *New()
{
return new vtkResliceCursorCallback;
}
void Execute(vtkObject *caller, unsigned long /*ev*/,
void *callData) override
{
vtkImagePlaneWidget* ipw =
dynamic_cast< vtkImagePlaneWidget* >(caller);
if (ipw)
{
double* wl = static_cast<double*>(callData);
if (ipw == this->IPW[0])
{
this->IPW[1]->SetWindowLevel(wl[0], wl[1], 1);
this->IPW[2]->SetWindowLevel(wl[0], wl[1], 1);
}
else if (ipw == this->IPW[1])
{
this->IPW[0]->SetWindowLevel(wl[0], wl[1], 1);
this->IPW[2]->SetWindowLevel(wl[0], wl[1], 1);
}
else if (ipw == this->IPW[2])
{
this->IPW[0]->SetWindowLevel(wl[0], wl[1], 1);
this->IPW[1]->SetWindowLevel(wl[0], wl[1], 1);
}
}
vtkResliceCursorWidget *rcw = dynamic_cast<
vtkResliceCursorWidget * >(caller);
if (rcw)
{
vtkResliceCursorLineRepresentation *rep = dynamic_cast<
vtkResliceCursorLineRepresentation * >(rcw->GetRepresentation());
vtkResliceCursor *rc = rep->GetResliceCursorActor()->
GetCursorAlgorithm()->GetResliceCursor();
for (int i = 0; i < 3; i++)
{
vtkPlaneSource *ps = static_cast< vtkPlaneSource * >(
this->IPW[i]->GetPolyDataAlgorithm());
ps->SetNormal(rc->GetPlane(i)->GetNormal());
ps->SetCenter(rc->GetPlane(i)->GetOrigin());
// If the reslice plane has modified, update it on the 3D widget
this->IPW[i]->UpdatePlacement();
//std::cout << "Updating placement of plane: " << i << " " <<
// rc->GetPlane(i)->GetNormal()[0] << " " <<
// rc->GetPlane(i)->GetNormal()[1] << " " <<
// rc->GetPlane(i)->GetNormal()[2] << std::endl;
//this->IPW[i]->GetReslice()->Print(cout);
//rep->GetReslice()->Print(cout);
//std::cout << "---------------------" << std::endl;
}
}
// Render everything
this->RCW[0]->Render();
}
vtkResliceCursorCallback() = default;
vtkImagePlaneWidget* IPW[3];
vtkResliceCursorWidget *RCW[3];
};
//----------------------------------------------------------------------------
int main(int argc, char *argv[])
{
vtkSmartPointer reader =
vtkSmartPointer::New();
reader->SetDataDimensions(64, 64);
reader->SetDataByteOrderToLittleEndian();
reader->SetImageRange(1, 93);
reader->SetDataSpacing(3.2, 3.2, 1.5);
reader->SetFilePrefix("D:/programe/VTK/mhd/肝脏/1/062d1");
reader->SetFilePattern("%s%30a.dcm");
reader->ReleaseDataFlagOn();
reader->SetDataMask(0x7fff);
reader->Update();
vtkSmartPointer<vtkOutlineFilter> outline =
vtkSmartPointer<vtkOutlineFilter>::New();
outline->SetInputConnection(reader->GetOutputPort());
vtkSmartPointer<vtkPolyDataMapper> outlineMapper =
vtkSmartPointer<vtkPolyDataMapper>::New();
outlineMapper->SetInputConnection(outline->GetOutputPort());
vtkSmartPointer<vtkActor> outlineActor =
vtkSmartPointer<vtkActor>::New();
outlineActor->SetMapper(outlineMapper);
vtkSmartPointer<vtkRenderer> ren[4];
vtkSmartPointer<vtkRenderWindow> renWin =
vtkSmartPointer<vtkRenderWindow>::New();
renWin->SetMultiSamples(0);
for (int i = 0; i < 4; i++)
{
ren[i] = vtkSmartPointer<vtkRenderer>::New();
renWin->AddRenderer(ren[i]);
}
vtkSmartPointer<vtkRenderWindowInteractor> iren =
vtkSmartPointer<vtkRenderWindowInteractor>::New();
iren->SetRenderWindow(renWin);
vtkSmartPointer<vtkCellPicker> picker =
vtkSmartPointer<vtkCellPicker>::New();
picker->SetTolerance(0.005);
vtkSmartPointer<vtkProperty> ipwProp =
vtkSmartPointer<vtkProperty>::New();
//assign default props to the ipw's texture plane actor
vtkSmartPointer<vtkImagePlaneWidget> planeWidget[3];
int imageDims[3];
reader->GetOutput()->GetDimensions(imageDims);
for (int i = 0; i < 3; i++)
{
planeWidget[i] = vtkSmartPointer<vtkImagePlaneWidget>::New();
planeWidget[i]->SetInteractor(iren);
planeWidget[i]->SetPicker(picker);
planeWidget[i]->RestrictPlaneToVolumeOn();
double color[3] = { 0, 0, 0 };
color[i] = 1;
planeWidget[i]->GetPlaneProperty()->SetColor(color);
planeWidget[i]->SetTexturePlaneProperty(ipwProp);
planeWidget[i]->TextureInterpolateOff();
planeWidget[i]->SetResliceInterpolateToLinear();
planeWidget[i]->SetInputConnection(reader->GetOutputPort());
planeWidget[i]->SetPlaneOrientation(i);
planeWidget[i]->SetSliceIndex(imageDims[i] / 2);
planeWidget[i]->DisplayTextOn();
planeWidget[i]->SetDefaultRenderer(ren[3]);
planeWidget[i]->SetWindowLevel(1358, -27);
planeWidget[i]->On();
planeWidget[i]->InteractionOn();
}
planeWidget[1]->SetLookupTable(planeWidget[0]->GetLookupTable());
planeWidget[2]->SetLookupTable(planeWidget[0]->GetLookupTable());
vtkSmartPointer<vtkResliceCursorCallback> cbk =
vtkSmartPointer<vtkResliceCursorCallback>::New();
// Create the reslice cursor, widget and rep
vtkSmartPointer< vtkResliceCursor > resliceCursor =
vtkSmartPointer< vtkResliceCursor >::New();
resliceCursor->SetCenter(reader->GetOutput()->GetCenter());
resliceCursor->SetThickMode(0);
resliceCursor->SetThickness(10, 10, 10);
resliceCursor->SetImage(reader->GetOutput());
vtkSmartPointer< vtkResliceCursorWidget > resliceCursorWidget[3];
vtkSmartPointer< vtkResliceCursorLineRepresentation > resliceCursorRep[3];
double viewUp[3][3] = { { 0, 0, -1 },{ 0, 0, 1 },{ 0, 1, 0 } };
for (int i = 0; i < 3; i++)
{
resliceCursorWidget[i] = vtkSmartPointer< vtkResliceCursorWidget >::New();
resliceCursorWidget[i]->SetInteractor(iren);
resliceCursorRep[i] =
vtkSmartPointer< vtkResliceCursorLineRepresentation >::New();
resliceCursorWidget[i]->SetRepresentation(resliceCursorRep[i]);
resliceCursorRep[i]->GetResliceCursorActor()->
GetCursorAlgorithm()->SetResliceCursor(resliceCursor);
resliceCursorRep[i]->GetResliceCursorActor()->
GetCursorAlgorithm()->SetReslicePlaneNormal(i);
const double minVal = reader->GetOutput()->GetScalarRange()[0];
if (vtkImageReslice *reslice =
vtkImageReslice::SafeDownCast(resliceCursorRep[i]->GetReslice()))
{
reslice->SetBackgroundColor(minVal, minVal, minVal, minVal);
}
resliceCursorWidget[i]->SetDefaultRenderer(ren[i]);
resliceCursorWidget[i]->SetEnabled(1);
ren[i]->GetActiveCamera()->SetFocalPoint(0, 0, 0);
double camPos[3] = { 0,0,0 };
camPos[i] = 1;
ren[i]->GetActiveCamera()->SetPosition(camPos);
ren[i]->GetActiveCamera()->ParallelProjectionOn();
ren[i]->GetActiveCamera()->SetViewUp(
viewUp[i][0], viewUp[i][1], viewUp[i][2]);
ren[i]->ResetCamera();
//ren[i]->ResetCameraClippingRange();
// Tie the Image plane widget and the reslice cursor widget together
cbk->IPW[i] = planeWidget[i];
cbk->RCW[i] = resliceCursorWidget[i];
resliceCursorWidget[i]->AddObserver(
vtkResliceCursorWidget::ResliceAxesChangedEvent, cbk);
// Initialize the window level to a sensible value
double range[2];
reader->GetOutput()->GetScalarRange(range);
resliceCursorRep[i]->SetWindowLevel(range[1] - range[0], (range[0] + range[1]) / 2.0);
planeWidget[i]->SetWindowLevel(range[1] - range[0], (range[0] + range[1]) / 2.0);
// Make them all share the same color map.
resliceCursorRep[i]->SetLookupTable(resliceCursorRep[0]->GetLookupTable());
planeWidget[i]->GetColorMap()->SetLookupTable(resliceCursorRep[0]->GetLookupTable());
}
// Add the actors
//
ren[0]->SetBackground(0.3, 0.1, 0.1);
ren[1]->SetBackground(0.1, 0.3, 0.1);
ren[2]->SetBackground(0.1, 0.1, 0.3);
ren[3]->AddActor(outlineActor);
ren[3]->SetBackground(0.1, 0.1, 0.1);
renWin->SetSize(600, 600);
//renWin->SetFullScreen(1);
ren[0]->SetViewport(0, 0, 0.5, 0.5);
ren[1]->SetViewport(0.5, 0, 1, 0.5);
ren[2]->SetViewport(0, 0.5, 0.5, 1);
ren[3]->SetViewport(0.5, 0.5, 1, 1);
// Set the actors' positions
//
renWin->Render();
ren[3]->GetActiveCamera()->Elevation(110);
ren[3]->GetActiveCamera()->SetViewUp(0, 0, -1);
ren[3]->GetActiveCamera()->Azimuth(45);
ren[3]->GetActiveCamera()->Dolly(1.15);
ren[3]->ResetCameraClippingRange();
vtkSmartPointer< vtkInteractorStyleImage > style =
vtkSmartPointer< vtkInteractorStyleImage >::New();
iren->SetInteractorStyle(style);
iren->Initialize();
iren->Start();
return 0;
}
按照上面的指示,我只是轻微的修改了细微的几个地方而已,
http://www.vtk.org/gitweb?p=VTK.git;a=blob_plain;f=Interaction/Widgets/Testing/Cxx/TestResliceCursorWidget2.cxx
这是官网的代码,可是我的代码却运行不了呢?麻烦大家解答。