[Insight-developers] Making cuberille polygonization work with QuadEdgeMesh

[Alexandre GOUAILLARD]

hi dan,

Excellent paper. and thanks for the well documented way of reproducing the code.

I will take a look this week end.

Internal point means that the point has already a full one-ring (which
means, it is already locally a 2 manifold and you cannot had anymore
face / edge without breaking the 2-manifoldness). it is possible that
you either produce non 2 manifolds or during the process at one stage
are not 2 manifold (add a new face before deleting a previous one).
If I understand this should not be as you are tessellating an implicit
surface, which gives me high hopes to be able to find a way to do it
in a QE (always 2 manifold at any point, any time) way.

alex.

On Fri, Jul 23, 2010 at 2:02 AM, Dan Mueller <dan.muel at gmail.com> wrote:
> 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
> _______________________________________________
> Powered by www.kitware.com
>
> Visit other Kitware open-source projects at
> http://www.kitware.com/opensource/opensource.html
>
> Kitware offers ITK Training Courses, for more information visit:
> http://kitware.com/products/protraining.html
>
> Please keep messages on-topic and check the ITK FAQ at:
> http://www.itk.org/Wiki/ITK_FAQ
>
> Follow this link to subscribe/unsubscribe:
> http://www.itk.org/mailman/listinfo/insight-developers
>