3D数据在用ITK读入后,需要在C++中做其他处理,这就需要数据转化的工作;
起先我先是试着转化到Mat数据,但是遇到很多问题,加上自己是初学者就没有什么人能问的,就自己摸索,这是我在另外的博客看到的用法,涉及数据转化、取值、赋值操作:不是完整的,用到哪些就取相应步骤好了
//读数据就不说了
seedImage = reader->GetOutput();
//ITK数据转换为C++数据
const unsigned short* seedData = seedImage->GetBufferPointer();//ITK数据转换为C++数据
//memset初始化结果图指针
//unsigned short* resultData = new unsigned //short[imgSize[0]*imgSize[1]*imgSize[2]];
//memset(resultData, 0, sizeof(unsigned short) * imgSize[0] * imgSize[1] * //imgSize[2]);
ImageType::IndexType point_tempx;
for(int k=0; k < seedImgSize[2]; k++)
for(int j = 0; j < seedImgSize[1]; j++)
for(int i = 0; i < seedImgSize[0]; i++)
{
if(seedData[k*seedImgSize[0]*seedImgSize[1]+j*seedImgSize[0]+i] > 阈值)
{ point.x = i;
point.y = j;
point.z = k;
ImageType::PixelType valuex = seedImage->GetPixel(point_tempx);
}
}
//C++数据转ITK数据
OutputImageType::IndexType newvoxelIndex;
for(int z = 0; z < imgSize[2]; ++z)
for(int y = 0; y < imgSize[1]; ++y)
for(int x = 0; x < imgSize[0]; ++x)
{
newvoxelIndex[0] = x;
newvoxelIndex[1] = y;
newvoxelIndex[2] = z;
if(y < sur1[z][x] || y > sur12[z][x])
dstImage->SetPixel(newvoxelIndex, 0);
else
dstImage->SetPixel(newvoxelIndex, resultData[z*size[0]*size[1]+y*size[0]+x]);
}
//输出结果图
typedef itk::ImageFileWriter<OutputImageType>WriterType;
WriterType::Pointer Writer = WriterType::New();
Writer->SetFileName(argv[3]);
Writer->SetInput(dstImage);
Writer->Update();
return 0;
}