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;
}全部评论 (0)
还没有任何评论哟~