第五章、vtk图像处理

vtk图像处理

5.1vtk图像创建

5.1.1vtk图像数据结构

数字图像文件内容有两个部分组成：图像头信息和数据。

图像头信息包括起点位置（Origin）、像素间隔（Space）和维数（Dimension）。这三个参数即可确定图像空间位置和大小。图像可看作空间中的一个规则的网格，网格中的每个最小单元称为像素（二维）或体素（三维），网格在每个方向上的像素或者体素个数即为图像在该方向上的维数。

下图表示一个4 $ \times $ 2 $ \times $ 3的体素的图像，即图像的位数，每个小球表示一个像素。图像原点为（5.1，10.0，6.5），两个像素之间的间隔表示像素间隔，每个方向上的像素间隔分别为1.5、1.5和1.8。

图像数据既为图像像素的像素值，一般采用一维数组来表示和存储。已知像素索引和图像维数的情况下，即可计算每个像素对应的像素值。医学图像中经常使用灰度图像。灰度图像的像素或体素中不需要考虑，其默认范围是0~255，可以采用一个unsigned char类型来表示。医学图像的会不范围要大于255，一般的范围为0~65536

5.1.2vtk图像创建

1.图像源Source

vtk中内置了多个创建图像的Source类，利用这些Source类可以快速创建图像。例如，用vtkImageCanvasSource2D为代表。该Source类的功能是创建画布，并提供多种几何图形的绘制填充功能。

vtkSmartPointer<vtkImageCanvasSource2D> canvas =
		vtkSmartPointer<vtkImageCanvasSource2D>::New();//定义一个vtkImageCanvasSource2D对象
	canvas->SetScalarTypeToUnsignedChar();
	canvas->SetNumberOfScalarComponents(1);
	canvas->SetExtent(0, 100, 0, 100, 0, 0);
	canvas->SetDrawColor(0, 0, 0, 0);
	canvas->FillBox(0,100,0,100);//绘制矩形
	canvas->SetDrawColor(255, 0, 0, 0);//设置矩形颜色
	canvas->FillBox(20,40,20,40);
	canvas->Update();

// Create actors
	vtkSmartPointer<vtkImageActor> redActor =
		vtkSmartPointer<vtkImageActor>::New();
	redActor->SetInputData(canvas->GetOutput());

	// Define viewport ranges
	// (xmin, ymin, xmax, ymax)
	double redViewport[4] = { 0.0, 0.0, 1.0, 1.0 };
	// Setup renderers
	vtkSmartPointer<vtkRenderer> redRenderer =
		vtkSmartPointer<vtkRenderer>::New();
	redRenderer->SetViewport(redViewport);
	redRenderer->AddActor(redActor);
	redRenderer->ResetCamera();
	redRenderer->SetBackground(1.0, 1.0, 1.0);

	// Setup render window
	vtkSmartPointer<vtkRenderWindow> renderWindow =
		vtkSmartPointer<vtkRenderWindow>::New();
	renderWindow->AddRenderer(redRenderer);
	renderWindow->SetSize(640, 480);
	renderWindow->Render();
	renderWindow->SetWindowName("ImageCanvasSource2D");

	// Setup render window interactor
	vtkSmartPointer<vtkRenderWindowInteractor> renderWindowInteractor =
		vtkSmartPointer<vtkRenderWindowInteractor>::New();
	vtkSmartPointer<vtkInteractorStyleImage> style =
		vtkSmartPointer<vtkInteractorStyleImage>::New();

	renderWindowInteractor->SetInteractorStyle(style);
	// Render and start interaction
	renderWindowInteractor->SetRenderWindow(renderWindow);
	renderWindowInteractor->Initialize();
	renderWindowInteractor->Start();

	return EXIT_SUCCESS;

效果：

其他vtk快速生成特定图像：

vtk类	功能
vtkImageEllipsoidSource	根据指定中心以及各个轴的半径来生成一个前景为椭圆的二值图像
vtkImageGuassianSource	生成一幅像素值服从高斯分布的图像
vtkImageGridSource	生成网格线图像
vtkImageNoiseSource	生成一个像素值为随机数的噪声图像
vtkImageSinusoidSource	生成图像像素值由正弦函数决定

5.2vtk图像显示

5.2.1vtkImageViewer2

vtkImageViewer2中封装了vtk图像显示的可视化渲染引擎。该类提供的主要交互操作有图像放缩、窗宽窗位调节，并提供切片选择及方向设置接口，尤其适合三维图像的切片显示。

vtkSmartPointer<vtkMetaImageReader> reader =
		vtkSmartPointer<vtkMetaImageReader>::New();
	reader->SetFileName(argv[1]);
	reader->Update();

	vtkSmartPointer<vtkImageViewer2> imageViewer =
		vtkSmartPointer<vtkImageViewer2>::New();
	imageViewer->SetInputConnection(reader->GetOutputPort());

	vtkSmartPointer<vtkRenderWindowInteractor> renderWindowInteractor =
		vtkSmartPointer<vtkRenderWindowInteractor>::New();
	imageViewer->SetupInteractor(renderWindowInteractor);

	imageViewer->SetColorLevel(500);
	imageViewer->SetColorWindow(2000);
	imageViewer->SetSlice(40);
	imageViewer->SetSliceOrientationToXY();
	imageViewer->Render();

	imageViewer->GetRenderer()->SetBackground(1.0, 1.0, 1.0);
	imageViewer->SetSize(640, 480);
	imageViewer->GetRenderWindow()->SetWindowName("DisplayImageExample");

	renderWindowInteractor->Start();

效果：

1.窗宽、窗位的概念

窗宽是图像显示的灰度范围。一般显示器的灰度范围为256级，而医学图像的灰度范围则远大于该范围，因此显示器不能显示所有灰度级，需要使用窗宽来定义欲现实的灰度范围。当灰度值高于该范围的最大值，均以白影显示。当低于该范围时，均以黑色显示。

2.医学图像二维视图

切片或切面时三维图像比较常见的概念，不同方向的切面都有特定的名字，分别是：矢状面，沿着身体前后径所做的与地面垂直的切面；冠状面，沿着身体左右径做的与地面垂直的切面；横断面，是指横断身体与地面平行的切面。

5.2.2 vtkImageActor

vtkImageActor是一个三维图像渲染Actor，通过纹理映射将图像映射到多边形上进行显示。

5.2.3 图像融合

使用vtkImageBlend实现两个图像融合：

	vtkSmartPointer<vtkJPEGReader> reader = vtkSmartPointer<vtkJPEGReader>::New();
	reader->SetFileName("E:\\PeopleImage\\wallhaven-newke8.jpg");
	reader->Update();

	vtkSmartPointer<vtkImageCanvasSource2D> imageSource = vtkSmartPointer<vtkImageCanvasSource2D>::New();
	imageSource->SetNumberOfScalarComponents(1);
	imageSource->SetScalarTypeToUnsignedChar();
	imageSource->SetExtent(0,512,0,512,0,0);
	imageSource->SetDrawColor(0.0);
	imageSource->FillBox(0, 512, 0, 512);
	imageSource->SetDrawColor(255.0);
	imageSource->FillBox(100, 400, 100, 400);
	imageSource->Update();

	vtkSmartPointer<vtkImageBlend> imageBlend = vtkSmartPointer<vtkImageBlend>::New();
	imageBlend->SetInputData(0, reader->GetOutput());
	imageBlend->SetInputData(1, imageSource->GetOutput());
	imageBlend->SetOpacity(0, 0.4);
	imageBlend->SetOpacity(1, 0.6);
	imageBlend->Update();

效果：

5.3 vtk图像基本操作

5.3.1图像信息的访问与修改

1.利用vtkImageData的方法

	vtkSmartPointer<vtkJPEGReader> reader = vtkSmartPointer<vtkJPEGReader>::New();
	reader->SetFileName("E:\\PeopleImage\\wallhaven-newke8.jpg");
	reader->Update();

	int dims[3];
	reader->GetOutput()->GetDimensions(dims);
	std::cout << "图像维数" << dims[0] << " " << dims[1] << " " << dims[2] << std::endl;

	double origin[3];
	reader->GetOutput()->GetOrigin(origin);
	std::cout << "图像原点" << origin[0] << " " << origin[1] << " " << origin[2] << std::endl;

	double spaceing[3];
	reader->GetOutput()->GetSpacing(spaceing);
	std::cout << "像素间隔：" << spaceing[0] << " " << spaceing[1] << " " << spaceing[2] << std::endl;