itk中的花式数据切割（二）

上篇切割的特征是越切越小，这次换个不会变小的方法，当然不同的方法有利有弊，中间取舍，看实际情况。

0.先构建一个与原始数据同等大小的图像，这个是前提条件

略。。。（不清楚的请转《itk中的基本图像操作》一文）

1.沿着Z轴按照预设范围切出方盒子（其实随便XYZ那一边都一样）

ImageType::SizeType size = input_data->GetLargestPossibleRegion().GetSize();int min_temp = point_z[2] - min;int max_temp = point_z[2] + max;for(int k= min_temp; k<max_temp; 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 = input_data->GetPixel(point_temp);output_data->SetPixel(point_temp,1);}}}

思路：构建同等大小的容器，然后将需要的像素点放进去。

优点：结果图像和原始图像一样大小，可以直接做数学运算

缺点：由于像素点相等，所以在遍历时候运算速度受到一定影响

2.切出序列图像的某一层

ImageType::RegionType inputRegion = input_data->GetLargestPossibleRegion();ImageType::SizeType size = inputRegion.GetSize();if (factor_mode == X){int x = num_temp;for (int y=0;y<size[1];y++){for (int z=0;z<size[2];z++){ImageType::IndexType  point_temp;point_temp[0] = x;point_temp[1] = y;point_temp[2] = z;ImageType::PixelType rad_value = input_data->GetPixel(point_temp);output_data->SetPixel(point_temp,rad_value);}}}else if (factor_mode == Y){int y = num_temp;for (int x=0;x<size[0];x++){for (int z=0;z<size[2];z++){ImageType::IndexType  point_temp;point_temp[0] = x;point_temp[1] = y;point_temp[2] = z;ImageType::PixelType rad_value = input_data->GetPixel(point_temp);output_data->SetPixel(point_temp,rad_value);}}}else if(factor_mode == Z){int z = num_temp;for (int x=0;x<size[0];x++){for (int y=0;y<size[1];y++){ImageType::IndexType  point_temp;point_temp[0] = x;point_temp[1] = y;point_temp[2] = z;ImageType::PixelType rad_value = input_data->GetPixel(point_temp);output_data->SetPixel(point_temp,rad_value);}}}

思路：和1中的方法如出一辙。只是提取了某一层。其实按照这方法，你可以把任何想提取的像素点放到对应的容器中。比如接下来的3。

3.切出一个球形（其实随便什么形状）

ImageType::RegionType inputRegion = input_data->GetLargestPossibleRegion();ImageType::SizeType size = inputRegion.GetSize();int radius2=radius*radius;for(int k=0;k<size[2];k++){for(int j=0;j<size[1];j++){for(int i=0;i<size[0];i++){if((i-x)*(i-x)+(j-y)*(j-y)+(k-z)*(k-z)<=radius2){ImageType::IndexType  point_temp2;point_temp2[0] = i;point_temp2[1] = j;point_temp2[2] = k;output_data->SetPixel(point_temp2,density);}}}}}

对比上篇《一》的方法，思路刚好是相反的，《一》中是在原始数据的基础上，把不需要的删除，而这篇是从原始数据中把需要的提取出来。还是那句话，方法没有优劣，选择要看实际情况。

道生一，一生二，二生三，三生万物。
人法地，地法天，天法道，道法自然。

代码之余，时刻要提醒自己跳出程序之外，不要太忙，不要把编码时间安排的太满，走一走，停一停，回头看一看，想一想。代码撸的再好，始终只是用来糊口的术，用到什么，就去学什么好了，随时拿起，随时放下。真正需要思考的，是怎么才能过好这一生，这才是道。