版权声明:本文为博主原创文章,转载请注明出处,谢谢 https://blog.csdn.net/rabbitbride/article/details/82422973
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;
}