1.直接访问图像像素（索引法）
1 #include <vtkAutoInit.h>
2 VTK_MODULE_INIT(vtkRenderingOpenGL);
3
4 #include <vtkSmartPointer.h>
5 #include <vtkImageData.h>
6 #include <vtkBMPReader.h>
7 #include <vtkImageViewer2.h>
8 #include <vtkRenderer.h>
9 #include <vtkRenderWindow.h>
10 #include <vtkRenderWindowInteractor.h>
11
12 int main()
13 {
14 vtkSmartPointer<vtkBMPReader> reader =
15 vtkSmartPointer<vtkBMPReader>::New();
16 reader->SetFileName("lena.bmp");
17 reader->Update();
18
19 int dims[3];
20 reader->GetOutput()->GetDimensions(dims);
21
22 int nbofComp;
23 nbofComp = reader->GetOutput()->GetNumberOfScalarComponents();
24
25 for (int k = 0; k < dims[2]; k++)
26 {
27 for (int j = 0; j < dims[1]; j++)
28 {
29 for (int i = 0; i < dims[0]; i++)
30 {
31 if (i < 384 && i > 128 && j > 128 && j < 384)
32 {
33 unsigned char *pixel = (unsigned char *)(reader->GetOutput()->GetScalarPointer(i, j, k));
34 *pixel = 255 - *pixel;
35 *(pixel + 1) = 255 - *(pixel + 1);
36 *(pixel + 2) = 255 - *(pixel + 2);
37 }
38 }
39 }
40 }
41
42 vtkSmartPointer<vtkImageViewer2> imgViewer =
43 vtkSmartPointer<vtkImageViewer2>::New();
44 imgViewer->SetInputData(reader->GetOutput());
45
46 vtkSmartPointer<vtkRenderWindowInteractor> rwi =
47 vtkSmartPointer<vtkRenderWindowInteractor>::New();
48 imgViewer->SetupInteractor(rwi);
49 imgViewer->Render();
50 imgViewer->GetRenderer()->ResetCamera();
51 imgViewer->Render();
52 imgViewer->GetRenderer()->SetBackground(1.0, 1.0, 1.0);
53 imgViewer->SetSize(640, 480);
54 imgViewer->GetRenderWindow()->SetWindowName("VisitImagePixelDirectly");
55
56 rwi->Start();
57
58 return 0;
59 }
输出结果：
上述案例实现了将图像的100*100大小的区域设置为反色。首先定义一个reader读取一副bmp图像，通过vtkImageData函数GetDimensions()获取图像的大小。建立三次循环，通过GetScalarPointer(i, j,k)函数获取访问图像像素值。需要注意的是，GetScalarPointer()函数返回的是void*类型，因此需要根据图像的实际类型进行强制转换。如上面代码中将像素值数组的头指针类型转换为unsigned char *。如果对于数据类型不确定的话，还可以先通过vtkImageCast将图像数据类型强制转换为特定的数据类型，再进行遍历。
在这里，我们需要注意的一点是VTK中彩色图以及矢量图的存储方式，具体如下：
因此在修改RGB图像以及向量图像像素时，需要根据像素的元组的组分数目来访问。上例中，需要修改每个像素的RGB值时，首先获得第(i, j, k)个像素的地址也就是R值的地址，然后将地址加1来访问后续G值以及B值。如果对于像素的元组组分不确定时，可以通过函数GetNumberOfScalarComponents()来获取。如下所示：
int nbOfComp = reader->GetOutput()->GetNumberOfScalarComponents();
2.迭代器方法访问图像像素
另外VTK中提供了vtkImageIterator类来利用迭代器方法访问图像像素。该类是一个模板类，使用时，需要提供迭代的图像像素类型以及迭代的区域大小。
1 #include <vtkAutoInit.h>
2 VTK_MODULE_INIT(vtkRenderingOpenGL);
3
4 #include <vtkSmartPointer.h>
5 #include <vtkBMPReader.h>
6 #include <vtkImageData.h>
7 #include <vtkImageIterator.h>
8 #include <vtkImageViewer2.h>
9 #include <vtkRenderer.h>
10 #include <vtkRenderWindow.h>
11 #include <vtkRenderWindowInteractor.h>
12
13 int main()
14 {
15 vtkSmartPointer<vtkBMPReader> reader =
16 vtkSmartPointer<vtkBMPReader>::New();
17 reader->SetFileName("lena.bmp");
18 reader->Update();
19
20 int subRegion[6] = { 64, 448, 64, 448, 0, 0 };
21 vtkImageIterator<unsigned char> iter(reader->GetOutput(),subRegion);
22
23 while (!iter.IsAtEnd())
24 {
25 unsigned char *inSI = iter.BeginSpan();
26 unsigned char *inSIEnd = iter.EndSpan();
27
28 while ( inSI != inSIEnd )
29 {
30 *inSI = 255 - *inSI;
31 ++inSI;
32 }
33 iter.NextSpan();
34 }
35
36 vtkSmartPointer<vtkImageViewer2> imgViewer =
37 vtkSmartPointer<vtkImageViewer2>::New();
38 imgViewer->SetInputConnection(reader->GetOutputPort());
39
40 vtkSmartPointer<vtkRenderWindowInteractor> rwi =
41 vtkSmartPointer<vtkRenderWindowInteractor>::New();
42 imgViewer->SetupInteractor(rwi);
43 imgViewer->Render();
44 imgViewer->GetRenderer()->ResetCamera();
45 imgViewer->Render();
46 imgViewer->GetRenderer()->SetBackground(1.0, 1.0, 1.0);
47 imgViewer->SetSize(640, 480);
48 imgViewer->GetRenderWindow()->SetWindowName("VisitImagePixelIteratively");
49
50 rwi->Start();
51
52 return 0;
53 }
输出结果：
如果对于ITK图像区域迭代器熟悉的话，可能会对上面代码存在疑问。上面代码中首先读取了一副bmp图像，然后定义了一个子区域。注意在定义子区域的时候，不要超过图像的大小范围。subRegion的六个值分别表示区域中x的最小最大值，y的最小最大值，z的最小最大值。由于处理的图像为二维图像，因此z的取值范围为[0,0]。然后根据图像类型unsigned char定义实例化一个图像迭代器it，定义it时有两个参数：一个是要访问的图像，另外一个是访问的图像区域。设置完毕后，迭代器开始工作。注意，上面代码中有两个while循环。
首先看第一个while循环，这里判断迭代器是否结束。进入循环后，对于每个迭代器it，又存在第二个循环。这个循环判断的是当前像素的组分是否迭代完毕。由于vtk中所有类型的图像格式都是vtkImageData，因此每个像素可能是标量，也可能是向量。因此，每当访问到一个像素时，需要迭代当前像素的组分。组分迭代时，inSI = it.BeginSpan()获取第一个组分，inSIEnd = it.EndSpan()表示组分迭代完毕，通过++inSI不断迭代组分，并对像素的组分值进行处理，当inSI与inSIEnd相等时组分迭代完毕。然后继续迭代像素it，直至迭代完毕所有像素。