VTK嵌入MFC同步显示

         使用VTK嵌入MFC，实现四视图更新，机制和细节参考原文。

         原文链接：http://blog.csdn.net/www_doling_net/article/details/8939115

 

原文代码：

#pragma once

#include <vtkResliceCursor.h>
#include <vtkResliceCursorWidget.h>

#include <vtkPlane.h>
#include <vtkPlaneSource.h>
#include <vtkPlaneWidget.h>

#include <vtkImagePlaneWidget.h>
#include <vtkResliceCursorThickLineRepresentation.h>
#include <vtkResliceCursor.h>

#include <vtkCommand.h>
#include <vtkViewport.h>
#include <vtkViewDependentErrorMetric.h>

#include <vtkSmartPointer.h>

#include <vtkRenderer.h>
#include <vtkRendererSource.h>
#include <vtkRenderingOpenGLModule.h>
#include <vtkRenderWindow.h>
#include <vtkWin32OpenGLRenderWindow.h>
#include <vtkWin32RenderWindowInteractor.h>
// CVtkView

class CVtkView : public CStatic
{
DECLARE_DYNAMIC(CVtkView)

public:
CVtkView();
virtual ~CVtkView();

public:
//3.2 重载CvtkView类PreSubclassWindow（）函数和OnPaint()函数
//PreSubclassWindow函数负责创建VTK可视化管线，OnPaint()函数负责客户区内场景渲染。
//vtkAcor,vtkRenderer,vtkRenderWindow,vtkRenderWindowInteractor四个部分。当然根据需要还可以设置vtkRenderWindowInteractorStyle,以及光照，材质，颜色等。
//在CvtkView类头文件中定义相关对象，并在PreSubclassWindow函数中实例化和构建可视化管线
void PreSubclassWindow();
void SetImageData(vtkSmartPointer<vtkImageData> ImageData);
void SetupReslice();

private:
vtkSmartPointer< vtkImagePlaneWidget >   m_ImagePlaneWidget;
vtkSmartPointer< vtkResliceCursorWidget> m_ResliceCursorWidget;
vtkSmartPointer< vtkResliceCursor >      m_ResliceCursor;
vtkSmartPointer< vtkResliceCursorThickLineRepresentation > m_ResliceCursorRep;

vtkSmartPointer<vtkRenderer>         m_Renderer;
vtkSmartPointer<vtkRenderWindow> m_RenderWindow;
vtkSmartPointer<vtkImageData>       m_ImageData;

//m_Direction为方向标志，取值分别为0,1和2，分别代表X轴，Y轴和Z轴方向，
int m_Direction;

protected:
DECLARE_MESSAGE_MAP()
};

 

原代码：

// VtkView.cpp : implementation file
//

#include "stdafx.h"
#include "iisu2vrMfcs.h"
#include "VtkView.h"


// CVtkView

IMPLEMENT_DYNAMIC(CVtkView, CStatic)

CVtkView::CVtkView()
{
//在实例化时需要注意，该视图类在默认情况下渲染的是vtkResliceCursorWidget对象的输出，
//因此需要为vtkResliceCursorWidget对象指定相应的vtkRenderer对象，
//m_ResliceCursorWidget->SetInteractor(m_RenderWindow->GetInteractor());
//m_ResliceCursorWidget->SetDefaultRenderer(m_Renderer);
}

CVtkView::~CVtkView()
{
}


void CVtkView::PreSubclassWindow()
{
// TODO: Add your specialized code here and/or call the base class  
CRect rect;
GetClientRect(rect);

m_Renderer = vtkSmartPointer<vtkRenderer>::New();
m_RenderWindow = vtkSmartPointer<vtkRenderWindow>::New();
m_RenderWindow->SetParentId(this->m_hWnd);
m_RenderWindow->SetSize(rect.Width(), rect.Height());
m_RenderWindow->AddRenderer(m_Renderer);

if (m_RenderWindow->GetInteractor() == NULL)
{
vtkSmartPointer<vtkRenderWindowInteractor> RenderWindowInteractor =
vtkSmartPointer<vtkRenderWindowInteractor>::New();
RenderWindowInteractor->SetRenderWindow(m_RenderWindow);
RenderWindowInteractor->Initialize();
}

m_RenderWindow->Start();
CStatic::PreSubclassWindow();
}

void CVtkView::SetImageData(vtkSmartPointer<vtkImageData> ImageData)
{
if (ImageData == NULL) return;

m_ImageData = ImageData;
SetupReslice();
}
void CVtkView::SetupReslice()
{
if (m_ImageData == NULL) return;
int dims[3];
m_ImageData->GetDimensions(dims);

//////////////////////////////////////////////////////////////////////////  
m_ImagePlaneWidget->SetInputData(m_ImageData);
m_ImagePlaneWidget->SetPlaneOrientation(m_Direction);
m_ImagePlaneWidget->SetSliceIndex(dims[m_Direction] / 2);
m_ImagePlaneWidget->On();
m_ImagePlaneWidget->InteractionOn();

//////////////////////////////////////////////////////////////////////////  
m_ResliceCursor->SetCenter(m_ImageData->GetCenter());
m_ResliceCursor->SetImage(m_ImageData);
m_ResliceCursor->SetThickMode(0);

//m_ResliceCursorRep->GetResliceCursorActor()->GetCursorAlgorithm()->SetResliceCursor(m_ResliceCursor);
//m_ResliceCursorRep->GetResliceCursorActor()->GetCursorAlgorithm()->SetReslicePlaneNormal(m_Direction);

m_ResliceCursorWidget->SetEnabled(1);
m_Renderer->ResetCamera();

//////////////////////////////////////////////////////////////////////////  
double range[2];
m_ImageData->GetScalarRange(range);
m_ResliceCursorWidget->GetResliceCursorRepresentation()->
SetWindowLevel(range[1] - range[0], (range[0] + range[1]) / 2.0);
m_ImagePlaneWidget->SetWindowLevel(range[1] - range[0], (range[0] + range[1]) / 2.0);
}

BEGIN_MESSAGE_MAP(CVtkView, CStatic)
END_MESSAGE_MAP()

添加类成员：

public:
afx_msg void OnBnClickedButtonViewshow();
//CStatic m_ShowLU;
//CStatic m_ShowRU;
//CStatic m_ShowLD;
//CStatic m_ShowRD;
CVtkView m_ShowLU;
CVtkView m_ShowRU;
CVtkView m_ShowLD;
CVtkView m_ShowRD;

CVtkView m_AxialView;
CVtkView m_SagittalView;
CVtkView m_CoronalView;
CVtkView m_3DView;

四视图更新代码：

//当用户改变图像切分的坐标轴时（旋转坐标轴或者平移坐标系），图像切分平面会产生相应的改变，
//如果将新的切分平面更新到二维视图的vtkImagePlaneWidget对象中，即可实现三维视图的同步更新操作。
///基于以上设计，实现一个vtkCommand子类，来监听vtkResliceCursorWidget::ResliceAxesChangedEvent消息，并实现相应的更新操作。
class vtkResliceCursorCallback : public vtkCommand
{
public:
static vtkResliceCursorCallback *New()
{
return new vtkResliceCursorCallback;
}

CVtkView* view[4];

public:
void Execute(vtkObject *caller, unsigned long /*ev*/,
void *callData)
{
vtkResliceCursorWidget *rcw = dynamic_cast<vtkResliceCursorWidget *>(caller);
if (rcw)
{
for (int i = 0; i < 3; i++)
{
vtkPlaneSource *ps =
static_cast<vtkPlaneSource *>(view[i]->GetImagePlaneWidget()->GetPolyDataAlgorithm());

ps->SetOrigin(view[i]->GetResliceCursorWidget()->
GetResliceCursorRepresentation()->GetPlaneSource()->GetOrigin());
ps->SetPoint1(view[i]->GetResliceCursorWidget()->
GetResliceCursorRepresentation()->GetPlaneSource()->GetPoint1());
ps->SetPoint2(view[i]->GetResliceCursorWidget()->
GetResliceCursorRepresentation()->GetPlaneSource()->GetPoint2());

view[i]->GetImagePlaneWidget()->UpdatePlacement();
view[i]->Render();
}
view[3]->Render();
}

}

vtkResliceCursorCallback() {}


};

最终效果图：