《itk实用demo》-PointsToImage和ImageToPoints

最新推荐文章于 2021-11-04 11:01:45 发布

原创最新推荐文章于 2021-11-04 11:01:45 发布 · 322 阅读

1 ·

CC 4.0 BY-SA版权

ITK医学图像处理专栏收录该内容

94 篇文章

订阅专栏

PointsToImage和ImageToPoints

ImageType::Pointer 
    PointsToImage(ImageType::Pointer image,PointSetType::Pointer pointset_temp)
{
    ImageType::RegionType inputRegion = image->GetLargestPossibleRegion();
    ImageType::SizeType size = inputRegion.GetSize();
    ImageType::IndexType start  = inputRegion.GetIndex();
    itk::ImageRegion<2> region(start, size);
    ImageType::Pointer outputImage = ImageType::New();
    outputImage->SetRegions(region);
    outputImage->Allocate();
    outputImage->FillBuffer(0);

    for (int i=0;i < pointset_temp->GetNumberOfPoints();i++)
    {
        PointType pos_last;
        pointset_temp->GetPoint( i, &pos_last );

        ImageType2D::IndexType index;
        index[0] = pos_last[0];
        index[1] = pos_last[1];

        outputImage->SetPixel(index, 1);
    }

    return outputImage;
}

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