ITK数据与C++数据转换

 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;
}