itk数据表达之 创建图像

ITK 支持任何像素类型和空间维的图像。

创建图像：（源代码在 Examples/DataRepresentation/Image/Image.cxx ）

这个例子阐述了如何人为地创建一个 itk::Image 类,下面是对图像类进行实例化、声明和创建的最简单程序代码:

实际上,很少直接给图像分配内存和对图像进行初始化,图像通常都是从一个源文件直接读取的。

 在 ITK 中,图像以一个或多个区域组合的形式存在。一个区域是图像的一个子集,并有可能是系统中其他类所占有的图像的一部分。一个比较普遍的区域是 LargePossibleRegion它是一个以完整部分定义的图像。其他重要的区域还有 BufferedRegion 它是内存中图像的一部分,和 RequestedRegion 它是在对图像进行操作时被滤波器或其他类要求的一部分。
如上所述,在 ITK 中人为的创建图像需要对图像进行实例化,并将描述图像的区域与图像结合起来。一个区域是有两个类来定义的:itk::Index 和 itk::Size 类。与图像结合的图像中的原始区域是有 Index 来定义的。区域的延伸或大小是有 Size 来定义的。Index 是有一个 n 维数列来表示的,在拓扑图像结构中这些表示图像最初的像素的数列的元素都是整数。当人为创建图像时,用户就需要定义图像的大小和图像的起始位置。有了这两个参数,就可以选择处理的区域。图像的起始点是有一个 Index 类定义的,这个类中存放了一个 n 维数列,数列中的元素都是整数,表示图像中各维上最初的像素值。

#include “itkImage.h”
int main(int, char *argv[])
{
typedef itk::Image< unsigned short, 3 > ImageType;
ImageType::Pointer image = ImageType::New();//然后就可以调用 New()操作创建图像并将结果分配到 itk::SmartPointer.

ImageType::IndexType start;
start[0] = 0; // first index on X
start[1] = 0; // first index on Y
start[2] = 0; // first index on Z
//区域大小是有一个相同大小的图像数列来表示的(使用 Size 类),数列中的元素是无符号整数表示图像像素在各个方向上的延伸。
ImageType::SizeType size;
size[0] = 200; // size along X
size[1] = 200; // size along Y
size[2] = 200; // size along Z
//定义了起始地址和图像大小这两个参数就可以创建一个 ImageRegion 对象,这个区域是有图像的起始地址和大小来初始化的。
ImageType::RegionType region;
region.SetSize( size );
region.SetIndex( start );
//最后,这个区域传递给图像对象来定义他的延伸和初始地址。SetRegion 方法同时设定了 LargePossibleRegion、BufferedRegion 和 RequestedRegion。注意到并未执行
//任何操作来给图像像素数据分配内存,调用 Allocate()来进行分配内存。直到给区域分配了足够的内存来存放信息分配都不需要任何指令。
image->SetRegions( region );
image->Allocate();
 return   0;        
                                                             }

 

 

#include <iostream>
#include "itkImage.h"
#include "itkImageFileWriter.h"

int main()
{
	typedef itk::Image< unsigned char, 2 > ImageType;
	ImageType::Pointer image = ImageType::New();

	ImageType::IndexType start;
	start[0] =   0;  // first index on X
	start[1] =   0;  // first index on Y

	ImageType::SizeType  size;
	size[0]  = 100;  // size along X
	size[1]  = 100;  // size along Y

	ImageType::RegionType region;  
	region.SetSize( size ); 
	region.SetIndex( start );
	image->SetRegions( region );
	image->Allocate();

	ImageType::PixelType  initialValue = 0;
	image->FillBuffer( initialValue );

	for (int i = 0; i < size[0]; i++)
	{
		for (int j = 0; j< size[1]; j++)
		{
			ImageType::IndexType pixelIndex;
			pixelIndex[0] = i;
			pixelIndex[1] = j;

			ImageType::PixelType   pixelValue = 	image->GetPixel( pixelIndex );
			if (pixelIndex[0]*pixelIndex[0] + pixelIndex[1]*pixelIndex[1] < 90 * 90 )
			{
				image->SetPixel( pixelIndex,  100); 
			}  
		}
	}

	typedef itk::ImageFileWriter< ImageType > WriterType;
	WriterType::Pointer writer = WriterType::New();

	// Or other file format, such as mhd/mha, png, jpg etc.
	writer->SetFileName("output.jpg"); 
	writer->SetInput(image);
	writer->Update();
	return 0;
}