Subject: CGAL users discussion list
List archive
- From: Jehyun Cha <>
- To:
- Subject: [cgal-discuss] CGAL 3D mesh generation prerequisite condition
- Date: Tue, 15 Mar 2016 19:22:25 -0700 (PDT)
- Authentication-results: mail2-smtp-roc.national.inria.fr; spf=None ; spf=SoftFail ; spf=None
- Ironport-phdr: 9a23:mTFtJxwt46DTnlXXCy+O+j09IxM/srCxBDY+r6Qd0e8SIJqq85mqBkHD//Il1AaPBtWLraodwLaH+4nbGkU+or+5+EgYd5JNUxJXwe43pCcHRPC/NEvgMfTxZDY7FskRHHVs/nW8LFQHUJ2mPw6anHS+4HYoFwnlMkItf6KuSt6U1538i7H60qaQSjsLrQL1Wal1IhSyoFeZnegtqqwmFJwMzADUqGBDYeVcyDAgD1uSmxHh+pX4p8Y7oGxmkql/posZCe26JPxgDO8QMDNzOG884IjntALIUBCUznoaSGQf1BRSUCbf6xSvfJrxozb8/sF90WHOJ8ryTLM9BW/7sY9kTRbpjGEMMDtvozKfsdB5kK8O+UHpnBd42YOBO9nNOQ==
I'm trying to triangulate given coronary artery model(please refer image and
file).
At first, I've tried to triangulate them using 3D constrained Delaunay
triangulation in TetGen engine, but it appears that TetGen didn't generate
them in all time. I've tried about 40 models with closed boundary, but only
half of them was successful.
As an alternative, I found that CGAL 3D mesh generation will generate
similar mesh based on Delaunay triangulation(of course, it's different from
3D constrained Delaunay triangulation).
I also tested it for 40 models which is same dataset used in TetGen test,
but it appears that only 1/4 of them were successful. It is weird because
even less models were processed than in TetGen test.
Is there are any condition for CGAL mesh generation except closed manifold
condition(no boundary & manifold)? Here is the code I've used in my test
case. It is almost same to example code from CGAL website.
========================================================
// Create input polyhedron
Polyhedron polyhedron;
std::ifstream input(fileName.str());
input >> polyhedron;
// Create domain
Mesh_domain domain(polyhedron);
// Mesh criteria (no cell_size set)
Mesh_criteria criteria(facet_angle = 25, facet_size = 0.15, facet_distance =
0.008,
cell_radius_edge_ratio = 3);
// Mesh generation
C3t3 c3t3 = CGAL::make_mesh_3<C3t3>(domain, criteria, no_perturb(),
no_exude());
cout << "Polygon finish: " << c3t3.number_of_cells_in_complex() << endl;
=========================================================
Here is one of CA model
<https://drive.google.com/file/d/0B81xvhk118U8VHdJTnM5amVmUVk/view?usp=sharing>
which was used in test case.
Also, I want to preserve given model triangles in generated mesh like
constrained Delaunay triangulation. Is there are any way that generate mesh
without specific criteria?
Please let me know if you want to know more.
--
View this message in context:
http://cgal-discuss.949826.n4.nabble.com/CGAL-3D-mesh-generation-prerequisite-condition-tp4661693.html
Sent from the cgal-discuss mailing list archive at Nabble.com.
- [cgal-discuss] CGAL 3D mesh generation prerequisite condition, Jehyun Cha, 03/16/2016
- Re: [cgal-discuss] CGAL 3D mesh generation prerequisite condition, Andreas Fabri, 03/16/2016
- Re: [cgal-discuss] CGAL 3D mesh generation prerequisite condition, Jehyun Cha, 03/16/2016
- Re: [cgal-discuss] CGAL 3D mesh generation prerequisite condition, Andreas Fabri, 03/16/2016
- Re: [cgal-discuss] CGAL 3D mesh generation prerequisite condition, Jehyun Cha, 03/16/2016
- Re: [cgal-discuss] CGAL 3D mesh generation prerequisite condition, Andreas Fabri, 03/16/2016
Archive powered by MHonArc 2.6.18.