itk中的基本图像操作

自己会慢慢把itk的一些简单用法放在博客来，就从创建一个简单图像开始吧。

//1、创建图像
ImageType::Pointer mask_data= ImageType::New();
//2、创建一个与image_data等大的mask图像
ImageType::PointType origin_temp = image_data->GetOrigin();
ImageType::SpacingType spacing_temp = image_data->GetSpacing();
ImageType::RegionType inputRegion = image_data->GetLargestPossibleRegion();
ImageType::SizeType size = inputRegion.GetSize();
ImageType::IndexType start  = inputRegion.GetIndex();
itk::ImageRegion<3> region(start, size);
mask_data->SetOrigin(origin_temp);
mask_data->SetSpacing(spacing_temp);
mask_data->SetRegions(region);
mask_data->Allocate();
mask_data->FillBuffer(itk::NumericTraits<PixelType>::Zero);
//3、复制一个图像
mask_data->Graft(image_data);
//4、给图像赋值
mask_data = image_data；
//5、遍历访问图像中的像素点
ImageType::RegionType inputRegion = mask_data->GetLargestPossibleRegion();
ImageType::SizeType size = inputRegion.GetSize();
for(int k=0;k<size[2];k++)
	{
		for(int j=0;j<size[1];j++)
		{
			for(int i=0;i<size[0];i++)
			{
				ImageType::IndexType  point_temp;
				point_temp[0] = i;
				point_temp[1] = j;
				point_temp[2] = k;
				ImageType::PixelType value = mask_data->GetPixel(point_temp);					
			}
		}
	}
//6、给像素点赋值
int density；
ImageType::IndexType  point_temp;
point_temp[0] = i;
point_temp[1] = j;
point_temp[2] = k;
mask_data->SetPixel(point_temp,density);	

itk是一个开源的图像算法包，任何算法都是在定义了一个基本结构的基础之上的，了解清楚基本结构，要比熟悉接口和用法更重要！

「紫薇软剑」三十岁前所用，误伤义士不祥，悔恨无已，乃弃之深谷。 重剑无锋，大巧不工。四十岁前恃之横行天下 。 四十岁后，不滞于物，草木竹石均可为剑。自此精修，渐进于无剑胜有剑之境。

——————《神雕侠侣》