itk中的网格Mesh
#include "itkMesh.h"
MeshType
::Pointer PointsToMesh(ImageType::Pointer sliceimage)
{
MeshType::Pointer mesh = MeshType::New();
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)
{
MeshType::PointType p0;
p0[0] = tmp_point[0];
p0[1] = tmp_point[1];
p0[2] = sliceIndex;
mesh->SetPoint(temp_index, p0 );
temp_index++;
}
}
}
return mesh;
}