itkquadedgemeshdeleteedgetest.cxx

DTMK软件开发包,此为开源软件,是一款很好的医学图像开发资源.

/*=========================================================================

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkQuadEdgeMeshDeleteEdgeTest.cxx,v $
  Language:  C++
  Date:      $Date: 2007-12-25 12:57:00 $
  Version:   $Revision: 1.6 $

  Copyright (c) Insight Software Consortium. All rights reserved.
  See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.

     This software is distributed WITHOUT ANY WARRANTY; without even 
     the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR 
     PURPOSE.  See the above copyright notices for more information.

=========================================================================*/
#if defined(_MSC_VER)
#pragma warning ( disable : 4786 )
#endif

#include "itkQuadEdgeMesh.h"

int itkQuadEdgeMeshDeleteEdgeTest( int , char* [] )
{
  typedef double                            PixelType;
  typedef itk::QuadEdgeMesh< PixelType, 3 > MeshType;
  std::string indent = "    ";
    
  MeshType::Pointer mesh = MeshType::New( );

  // Points
  MeshType::PointType p0, p1, p2, p3, p4, p5;
  p0[ 0 ] =  0.00000000000000; p0[ 1 ] =  0.00000000000000; p0[ 2 ] = 5.0;
  p1[ 0 ] =  0.00000000000000; p1[ 1 ] = 10.00000000000000; p1[ 2 ] = 0.0;
  p2[ 0 ] = -9.51056516295153; p2[ 1 ] =  3.09016994374947; p2[ 2 ] = 0.0;
  p3[ 0 ] = -5.87785252292473; p3[ 1 ] = -8.09016994374947; p3[ 2 ] = 0.0;
  p4[ 0 ] =  5.87785252292473; p4[ 1 ] = -8.09016994374948; p4[ 2 ] = 0.0;
  p5[ 0 ] =  9.51056516295154; p5[ 1 ] =  3.09016994374947; p5[ 2 ] = 0.0;

  MeshType::PointIdentifier pid0 = mesh->AddPoint( p0 );
  MeshType::PointIdentifier pid1 = mesh->AddPoint( p1 );
  MeshType::PointIdentifier pid2 = mesh->AddPoint( p2 );
  MeshType::PointIdentifier pid3 = mesh->AddPoint( p3 );
  MeshType::PointIdentifier pid4 = mesh->AddPoint( p4 );
  MeshType::PointIdentifier pid5 = mesh->AddPoint( p5 );

  // Cells in a proper way
  mesh->AddEdge( pid3, pid4 );
  mesh->AddEdge( pid4, pid0 );
  mesh->AddEdge( pid0, pid3 );
  mesh->AddFaceTriangle( pid3, pid4, pid0 );

  mesh->AddEdge( pid4, pid5 );
  mesh->AddEdge( pid5, pid0 );
  mesh->AddFaceTriangle( pid4, pid5, pid0 );

  mesh->AddEdge( pid5, pid1 );
  mesh->AddEdge( pid1, pid0 );
  mesh->AddFaceTriangle( pid5, pid1, pid0 );

  mesh->AddEdge( pid1, pid2 );
  mesh->AddEdge( pid2, pid0 );
  mesh->AddEdge( pid2, pid3 );

  int EdgesBefore = mesh->ComputeNumberOfEdges();

  // Deleting two arbitrary edges: 
  mesh->DeleteEdge( pid3, pid4 );
  mesh->DeleteEdge( pid0, pid5 );

  std::cout << indent << "Trying to remove only two edges...";
  if ( EdgesBefore - mesh->ComputeNumberOfEdges() == 2 )
    {
    std::cout << "OK." << std::endl;
    return( EXIT_SUCCESS );
    }
  else
    {
    std::cout << "FAILED." << std::endl;
    return( EXIT_FAILURE );
    }

  return( EXIT_SUCCESS );
}