天天看點

《itk實用demo》-itkMesh

itk中的網格Mesh
#include "itkMesh.h"

MeshType
    ::Pointer PointsToMesh(ImageType::Pointer sliceimage)
{
    MeshType::Pointer  mesh = MeshType::New();

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

    int temp_index = ;
    for(int i=; i<size[]; i++)
    {
        for(int j=; j<size[]; j++)
        {
            ImageType::IndexType tmp_point;
            tmp_point[] = i;
            tmp_point[] = j;
            float gray_value = sliceimage->GetPixel(tmp_point);
            if (gray_value!=)
            {
                MeshType::PointType p0;
                p0[] = tmp_point[]; 
                p0[] = tmp_point[]; 
                p0[] = sliceIndex;
                mesh->SetPoint(temp_index, p0 );
                temp_index++;
            }
        }
    }
    return mesh;
}