[Insight-developers] Making cuberille polygonization work with QuadEdgeMesh

[Dan Mueller]

Hi Insight Developers,

I recently submitted an Insight Journal article describing the
"cuberille" image to mesh filter:
    http://www.insight-journal.org/browse/publication/740

This filter seems to be working fine with itk::Mesh, but does not
always work with itk::QuadEdgeMesh (for some datasets it works fine;
others cause warnings/errors).

To reproduce the issue, follow these steps:
  1. Download the article source code:
      http://www.insight-journal.org/download/sourcecode/740/3
  2. Open Testing/CuberilleTest01.cxx
      Change the line "#define USE_QUAD_EDGE_MESH 0"
      to "#define USE_QUAD_EDGE_MESH 1"
  3. Configure, generate and build the "Cuberille" project
  4. Run the tests.

With quad edge mesh enabled I get the following test results:
===============
1>------ Build started: Project: RUN_TESTS, Configuration: Release Win32 ------
1>
1>Performing Post-Build Event...
1>Start processing tests
1>Test project C:/Work/Research/Polygonization/Cuberille/Build
1>  1/ 17 Testing Cuberille_Blob0_00               Passed
1>  2/ 17 Testing Cuberille_Blob1_01               Passed
1>  3/ 17 Testing Cuberille_Blob2_01               Passed
1>  4/ 17 Testing Cuberille_MarschnerLobb_01    ***Failed
1>  5/ 17 Testing Cuberille_Fuel_01                Passed
1>  6/ 17 Testing Cuberille_Fuel_02                Passed
1>  7/ 17 Testing Cuberille_Fuel_03                Passed
1>  8/ 17 Testing Cuberille_HydrogenAtom_01        Passed
1>  9/ 17 Testing Cuberille_Neghip_01           ***Failed
1> 10/ 17 Testing Cuberille_Neghip_02           ***Failed
1> 11/ 17 Testing Cuberille_Neghip_03           ***Failed
1> 12/ 17 Testing Cuberille_Nucleon_01             Passed
1> 13/ 17 Testing Cuberille_Nucleon_02             Passed
1> 14/ 17 Testing Cuberille_Nucleon_03             Passed
1> 15/ 17 Testing Cuberille_Silicium_01            Passed
1> 16/ 17 Testing Cuberille_Silicium_02            Passed
1> 17/ 17 Testing Cuberille_Silicium_03            Passed
1>76% tests passed, 4 tests failed out of 17
1>The following tests FAILED:
1>   4 - Cuberille_MarschnerLobb_01 (Failed)
1>   9 - Cuberille_Neghip_01 (Failed)
1>  10 - Cuberille_Neghip_02 (Failed)
1>  11 - Cuberille_Neghip_03 (Failed)
===============

I get the following warnings/errors:
===============
Debug: In c:\insighttoolkit-3.18.0\code\review\itkGeometricalQuadEdge.txx,
line 312
 (016192E0): Internal point.

Debug: In c:\insighttoolkit-3.18.0\code\review\itkGeometricalQuadEdge.txx,
line 399
 (016192E0): This point is yet surrounded by faces.

Debug: In c:\insighttoolkit-3.18.0\code\review\itkGeometricalQuadEdge.txx,
line 312
 (016192E0): Internal point.

Debug: In c:\insighttoolkit-3.18.0\code\review\itkGeometricalQuadEdge.txx,
line 399
 (016192E0): This point is yet surrounded by faces.

Debug: In c:\insighttoolkit-3.18.0\code\review\itkGeometricalQuadEdge.txx,
line 312
 (017470A0): Internal point.
...
===============

Obviously I am violating one or more QuadEdgeMesh requirements. Can
anyone help me to determine how to solve these errors? Specifically:
what is an internal point? Why is the point "yet surrounded by
faces?".

Any help would be appreciated.

Cheers, Dan
dan dot muel at gmail dot com