VTK作为医学显示库,得到较多使用。作为较为上层的设计,对OpenGL进行了封装,并且有Windows、Linux、安卓等开发版本,可移植性较强。

          不过VES暂时没有编译成功。

          以下是嵌入MFC-IDC代码,有个缺陷,暂时没能修改。

代码:

showPointsInVtk(short* &depthImage, int HeightDepth, int WidthDepth,int IDCX)
{
	
	//double arr[3];
	vtkPoints * points = vtkPoints::New();
	int n = 0; int idx = 0; int Length = 3 * HeightDepth*WidthDepth;
	//首先读取点云数据到点表points同时指定点对应的id:
	for (int i = 0; i < HeightDepth; ++i){
		int LastPixel = i* WidthDepth;
		for (int j = 0; j < WidthDepth; ++j){
			int y = depthImage[3 * (j + LastPixel)];
			int x = depthImage[3 * (j + LastPixel) +1];
			int z = depthImage[3 * (j + LastPixel) +2];
			points->InsertPoint(n, x, y, z);
			n++;
		}
	}
	printf("%d\n", n);

	vtkPolyVertex * polyvertex = vtkPolyVertex::New();
	polyvertex->GetPointIds()->SetNumberOfIds(n);

        int i = 0; //建立拓扑关系
        for (i = 0; i < n; i++){
		polyvertex->GetPointIds()->SetId(i, i);
	}

	vtkUnstructuredGrid * grid = vtkUnstructuredGrid::New();
	grid->SetPoints(points);
	grid->InsertNextCell(polyvertex->GetCellType(),
		polyvertex->GetPointIds());

	vtkDataSetMapper *map1 = vtkDataSetMapper::New();
	map1->SetInputData(grid);

	vtkActor *actor1 = vtkActor::New();
	actor1->SetMapper(map1);
	actor1->GetProperty()->SetColor(0.194, 0.562, 0.75);

	vtkRenderer *ren = vtkRenderer::New();
	ren->AddActor(actor1);
	ren->SetBackground(0, 0, 0);

	//renWin和iren分别是vtkWin32OpenGLRenderWindow, vtkWin32RenderWindowInteractor
	vtkWin32RenderWindowInteractor        *iren = vtkWin32RenderWindowInteractor::New();
	vtkWin32OpenGLRenderWindow             *win = vtkWin32OpenGLRenderWindow::New();
	vtkInteractorStyleTrackballCamera    *style = vtkInteractorStyleTrackballCamera::New();
	win->Clean();
	win->SetSize(480, 360);
	win->BordersOn();
	win->AddRenderer(ren);
	iren->SetRenderWindow(win);
	iren->SetInteractorStyle(style);

	//搭上一个显示用的MFC窗口框架。
         CWnd *viewer_Wnd;
        viewer_Wnd = this->GetDlgItem(IDCX);   
         HWND hWnd = (HWND)viewer_Wnd->m_hWnd;
         HWND hParent = ::GetParent(hWnd);
        //::ShowWindow(hParent, SW_HIDE);
        ::ShowWindow(hParent, 1);
        win->SetParentId(this->m_hWnd);

        iren->Start();

        map1->Delete();
        grid->Delete();
       actor1->Delete();
        polyvertex->Delete();

        ren->Delete();
        win->Delete();
        iren->Delete();

        return 1;
}
 
 

显示效果:

 


 注意事项:

使用

iren->Start();
导致renter不能更行,修改为以下代码可以完成动态更新:

把 render和win放在初始对话框初始化,其他代码修改为:

this->win->Render();

//在 对话框初始化时初始化Render
inline int Ciisu2vrMfcsDlg::showPointsInVtkG(short* &depthImage, int HeightDepth, int WidthDepth, int IDCX)
{

	//double arr[3];
	vtkPoints * points = vtkPoints::New();
	int n = 0; int idx = 0; int Length = 3 * HeightDepth*WidthDepth;
	//首先读取点云数据到点表points同时指定点对应的id:
	for (int i = 0; i < HeightDepth; ++i){
		int LastPixel = i* WidthDepth;
		for (int j = 0; j < WidthDepth; ++j){
			int y = depthImage[3 * (j + LastPixel)];
			int x = depthImage[3 * (j + LastPixel) + 1];
			int z = depthImage[3 * (j + LastPixel) + 2];//这里压缩四倍是比较好的方法//int z = (depthImage[(j + LastPixel)]) ;
			points->InsertPoint(n, x, y, z);
			n++;
		}
	}
	vtkPolyVertex * polyvertex = vtkPolyVertex::New();
	polyvertex->GetPointIds()->SetNumberOfIds(n);
	int i = 0;
	//建立拓扑关系
	for (i = 0; i < n; i++){
		polyvertex->GetPointIds()->SetId(i, i);
	}

	vtkUnstructuredGrid * grid = vtkUnstructuredGrid::New();
	grid->SetPoints(points);
	grid->InsertNextCell(polyvertex->GetCellType(),
		polyvertex->GetPointIds() );

	vtkDataSetMapper *map = vtkDataSetMapper::New();
	map->SetInputData(grid);

	vtkActor *actorG = vtkActor::New();
	actorG->SetMapper(map);
	actorG->GetProperty()->SetColor(0.194, 0.562, 0.75);

	vtkRenderer *ren = vtkRenderer::New();
	ren->AddActor(actorG);
	ren->SetBackground(0.13, 0.40, 0.15);

	this->win->AddRenderer(ren);
	this->win->Render();//自己写一个重载-更新函数!

	map->Delete();
	grid->Delete();
	actorG->Delete();
	polyvertex->Delete();

	return 1;
}

后记:

        效果还不错,就是和OpenCV结合时有点麻烦,命名空间冲突。

Logo

更多推荐