注:此文知识学习笔记,仅记录完整程序和实现结果,具体原理参见:
https://blog.csdn.net/www_doling_net/article/details/8541534
1.vtkImageData实现图像信息的访问与修改
vtkImageData提供了多个函数用于访问或者获取图像的基本信息,这些函数通常使用Set或者Get加上相应 的信息名的形式,列如获取图像维数的方法定义为GetDimensions().
#include <vtkSmartPointer.h>
#include <vtkBMPReader.h>
#include <vtkImageData.h>
#include <vtkImageViewer2.h>
#include <vtkRenderer.h>
#include <vtkRenderWindow.h>
#include <vtkRenderWindowInteractor.h>
int main(){
//读取数据
vtkSmartPointer<vtkBMPReader> reader = vtkSmartPointer<vtkBMPReader>::New();
reader->SetFileName("test.bmp");
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 spacing[3];
reader->GetOutput()->GetSpacing(spacing);
std::cout << "图像间距:" << spacing[0] << " " << spacing[1] << " " << spacing[2] << std::endl;
//显示图像
vtkSmartPointer<vtkImageViewer2> imgViewer = vtkSmartPointer<vtkImageViewer2>::New();
imgViewer->SetInputConnection(reader->GetOutputPort());
vtkSmartPointer<vtkRenderWindowInteractor> rwi = vtkSmartPointer<vtkRenderWindowInteractor>::New();
imgViewer->SetupInteractor(rwi);
imgViewer->Render();
imgViewer->GetRenderer()->ResetCamera();
imgViewer->Render();
imgViewer->GetRenderer()->SetBackground(1.0, 1.0, 1.0);
imgViewer->SetSize(640, 480);
imgViewer->GetRenderWindow()->SetWindowName("GetImageInfo");
rwi->Start();
return 0;
}
运行结果:
2. 类vtkImageChangeInformation用于图像修改
vtkImageData中提供了多个Set函数用于设置图像的基本信息。当对一个管线的输出修改图像信息后,如果管线重新Update,那么这些修改都会恢复回原来的值。而vtkChangeImageInformation可以作为管线中的一个filter来修改图像信息。利用这个filter可以修改图像的原点,像素间隔以及范围起点(extent),另外还可以对图像平移缩放等操作。
#include"vtkSmartPointer.h"
#include"vtkBMPReader.h"
#include"vtkImageData.h"
#include"vtkImageChangeInformation.h"
#include"vtkImageViewer2.h"
#include"vtkRenderer.h"
#include"vtkRenderWindow.h"
#include"vtkRenderWindowInteractor.h"
using namespace std;
int main(){
vtkSmartPointer<vtkBMPReader> reader = vtkSmartPointer<vtkBMPReader>::New();
reader->SetFileName("test.bmp");
reader->Update();
int dims[3];
double origin[3];
double spacing[3];
reader->GetOutput()->GetDimensions(dims);
cout << "图像维度" << dims[0] << " " << dims[1] << " " << dims[2] << endl;
reader->GetOutput()->GetOrigin(origin);
cout << "图像中心" <<origin[0] << " " << origin[1] << " " << origin[2] << endl;
reader->GetOutput()->GetSpacing(spacing);
cout << "图像间隔" << spacing[0] << " " << spacing[1] << " " << spacing[2] << endl;
vtkSmartPointer<vtkImageChangeInformation> changer = vtkSmartPointer<vtkImageChangeInformation>::New();
changer->SetInputData(reader->GetOutput());
changer->SetOutputOrigin(100, 100, 0);
changer->SetOutputSpacing(10,10,1);
changer->SetCenterImage(1);
changer->Update();
changer->GetOutput()->GetDimensions(dims);
cout << "图像维度" << dims[0] << " " << dims[1] << " " << dims[2] << endl;
changer->GetOutput()->GetOrigin(origin);
cout << "图像中心" << origin[0] << " " << origin[1] << " " << origin[2] << endl;
changer->GetOutput()->GetSpacing(spacing);
cout << "图像间隔" << spacing[0] << " " << spacing[1] << " " << spacing[2] << endl;
vtkSmartPointer<vtkImageViewer2> imgViewer = vtkSmartPointer<vtkImageViewer2>::New();
imgViewer->SetInputConnection(changer->GetOutputPort());
vtkSmartPointer<vtkRenderWindowInteractor> rwi = vtkSmartPointer<vtkRenderWindowInteractor>::New();
imgViewer->SetupInteractor(rwi);
imgViewer->Render();
imgViewer->GetRenderer()->ResetCamera();
imgViewer->GetRenderer()->SetBackground(1.0,1.0,1.0);
imgViewer->Render();
imgViewer->GetRenderWindow()->SetWindowName("ChangeImageInfo");
imgViewer->SetSize(640,480);
rwi->Start();
return EXIT_SUCCESS;
}
运行结果:
这里的解释是参考《VTK图形图像开发进阶》,用lena图像跑出来的图像中心是(-2555,-2555, 0)。
首先读入图像,由vtkImageData提供函数接口获取图像的维数,原点和像素间隔。然后定义vtkImageChangeInformation指针,并设置输出图像原点为(100, 100, 0),输出图像像素间隔为(10, 10, 1),然后调用CenterImage()函数将图像的原点置于图像的中心。
操作后的结果使得图像的原点位于(-2555, -2555, 0),SetOutputOrigin(100, 100, 0)并没有起作用。
原因在哪里呢?
如果看下CenterImage()函数的注释,可以发现该函数的作用是将(0, 0, 0)点置于图像的中心。当CenterImage该函数执行时会重写SetOutputOrigin(),所以SetOutputOrigin函数不会产生任何作用。那(-2555, -2555, 0)又是如何计算出来的呢?
这个原理如下图所示:
根据图像的维数和像素间隔计算得到新的图像的宽度和高度为(512-1)*10,初始图像的原点位于(0, 0,0),现在将图像的中心平移至原点,平移量为(-(512-1)*10/2,(512-1)*5/2, 0) = (-2555, -2555,0)。
参考资料:
1.《The Visualization Toolkit – AnObject-Oriented Approach To 3D Graphics (4th Edition)》
2. 张晓东, 罗火灵. VTK图形图像开发进阶[M]. 机械工业出版社, 2015.