《itk实用demo》-ImageToPoints

版权声明：本文为博主原创文章，转载请注明出处，谢谢 https://blog.csdn.net/rabbitbride/article/details/82422766

ImageToPoints

PointSetType
    ::Pointer ImageToPoints(ImageType::Pointer sliceimage)
{    
    PointSetType::Pointer  PointSet = PointSetType::New();
    PointsContainerPointer  points = PointSet->GetPoints();

    int size[2];
    size[0]=sliceimage->GetLargestPossibleRegion().GetSize()[0];
    size[1]=sliceimage->GetLargestPossibleRegion().GetSize()[1];

    int temp_index = 0;
    for(int i=0; i<size[0]; i++)
    {
        for(int j=0; j<size[1]; j++)
        {
            ImageType::IndexType tmp_point;
            tmp_point[0] = i;
            tmp_point[1] = j;
            float gray_value = sliceimage->GetPixel(tmp_point);
            if (gray_value!=0)
            {
                PointType p0;
                p0[0] = tmp_point[0];
                p0[1] = tmp_point[1];
                p0[2] = sliceIndex;
                points->InsertElement(temp_index, p0);
                temp_index++;
            }
        }
    }
    return PointSet;
}