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() {}

};

最终效果图：