Skip to Content.
Sympa Menu

cgal-discuss - Re: [cgal-discuss] Create a coarse surface mesh produces areas of very small elements

Subject: CGAL users discussion list

List archive

Re: [cgal-discuss] Create a coarse surface mesh produces areas of very small elements


Chronological Thread 
  • From: Sebastien Loriot <>
  • To:
  • Subject: Re: [cgal-discuss] Create a coarse surface mesh produces areas of very small elements
  • Date: Wed, 24 May 2023 14:34:11 +0200
  • Authentication-results: mail3-smtp-sop.national.inria.fr; spf=None ; spf=Pass ; spf=None
  • Ironport-data: A9a23:YUMseKu/BmYp6X7QqLj3dXWy7OfnVJ5aMUV32f8akzHdYApBsoF/q tZmKWvSOq3bZTH3c9AjPoywoxkEuMDSmIAwT1BkrC8zHyMRgMeUXt7xwmXYb3rDdJWbJK5Ex 5xDMYeYdJhcolv0/ErF3m3J9CEkvU2wbuOgTrSCYEidfCc8IA85kxVvhuUltYBhhNm9Emult Mj7yyHlEAbNNwVcbCRMs8pvlDs15K6p4G5C4gRlDRx2lAa2e0c9XMp3yZ6ZdCOQrrl8RoaSW +vFxbelyWLVlz9F5gSNz94X2mVTKlLjFVDmZkh+A8BOsTAezsAG6ZvXAdJHAathZ5plqPgqo DlFncTYpQ7EpcQgksxFO/VTO3kW0aGrZNYriJVw2CCe5xSuTpfi/xlhJHAnI4hB6uImOGxLy b8mCxc8RTrS2f3jldpXSsE07igiBMziPYdapWs5iD+FUrApRpfMR6iM7thdtNsyrpoWTLCOO oxDM2MpMEuojx5nYj/7DLolmOCyh37jNThcgF2QrKszpWPUyWSd1ZCwboqKIofWHa25mG69q WDd2yPzCy0GJdCdzmWewi6crNXAyHaTtIU6TeXkrJaGmma7zWMaDFgaVECwvOKiolWvXspWb U0S4Csn66YonHFHVfH4Vhy85WGe51sSB4UWHOo95wWAjKHT5m51G1ToUBZwecIj7d8yFAUn0 2PKofDONBdF4eOKHCf1GqivkRu+Pi0cLGknbCACTBcY79SLnG3VpkKfJjqEOP7l5uAZCQ0c0 BjR83dj3+R7YdojkvTkrQqe0lpAs7CQFlZtjjg7SF5J+e+QWWJIT4mh6Fye8uwZaYjFER+Ou 38Ln8XY5+cLZX1sqMBvaLRUdF1Kz6zdWNE5vbKJN8d9n9hK0yD4Fb28GBkkeC9U3j8sIFcFm nP7twJL/4N0N3C3d6JxaI/ZI510nfSwS4q7Da6INIsmjn1NmOmvrHEGiam4jzCFraTQufxX1 WqzK5f3US5BVcyLMhLqHL1AuVPU+szO7TqLGcqTI+WP3r2Za3qYIYrpw3PfBt3VGJis+V2Pm /4Gb5Xi40wGDIXWPHOLmaZNdgtiBSZgVfjLRzl/LLHrzvxOQzF/VZc8ANoJJ+RYokiivr6Yr i3jARcEkQaXaL+uAVziV02PoYjHBf5XxU/X9wR1Vbpx8yl4OdSc/+0EeoEpfLIq0uVmwLQmB 7MGYsiMSLAHADjO5z1XP9G3oZ1AZSabo1uEHxOkRzwjIL9mZQjCoeH/ciXVqSIhMyuQtOkFm YOG6D/1e5Q4aj5ZPJ7kU873l1KVlloBqd12RHrNc4Vyel2z0Y1EKB7Rr/4QIuMKIyrt3jGxi gScW08ZgcLvoIYF1sbDqo7ZjoWuEsp4RlF7GUuC55mIFCDqxEiR6q4eb/SpJBf2S3HR1Ji5Q 9ls39XQEaEijUlblYhRCJNpxv8O3MTurLpk0QhUJnXHQFC1ALdGIHPd/81wmoBS571eqy2kc 1mu/4REBLC3J8/VKl4dCw47ZOCl1/tPuD3z7+wwEXrq9h1M47uLfkVDDSai0BUHAuNOD7ok5 uM9tOo9yQ+10EMqO+nbqBFkzT2HK3hYXpg3sp0fPpTQtTMq7VN8MLj8ES795a+dZ+pcancKJ iCmv4ucprB+6HebTV8NOyns59dNvbUPpxFA83EaLXurhNfup6E6zT9RwxsNXyVXyRRN7MxrM EMybEZ3CLm8/Q4zoM1PQWr2FxpzP0Cb8Bapyn8itm7Qf2+3XEPjcUw/Pues+hgC0mR+JzJ0w pCR+FzHYx3LIv7j+zQUWFF0jcDjQfha1BzwqOr+E+urR5AFMCfY2Imwbm82mj7bKMIWhnyfg 9J1/ex1OJbJBQRJr4IVU4ClhKktEjaaL2l/QNZkzqMDPUfYXBqQgTGuCUSASvlhFszw032TK pJRf5pUdhGEyiyxgCgRBvcML59KjfcZ3oc+VY2xF1EWkYm0j2RPgMrL+znclV0bZYxksfwAJ 7P7cxOAFW2thkVopVLdkfkcBE2GZYgrWQ6t+sG07+QDKLwbus5Ob0wZ8+W5rleVAiRd7jOWu wLxPfbWxtN91LU2zpfNE7pCNSqwO9jcROSFyyHtktVsPPfkE9bCiBMRkXbjZz9pBLo2X89ms 4iCq/vl9RrhkItudlvGir6tMrJs5/SieMZ2aeXJdGJ7mwmGU+/SuyozwXiyc8F1oYkM9/uZS BudQ+ruUMweRPN25mBfMgpaGDYjU5XHVL/q/36Bnq7dGyom8FL1Kf28/iXUdkBdTCgDPqP+B iLSu/qD4tN5rpxGNCQbBsNJUoNJH1v+ZZQIL9HBlyGULm2NsGOwvrHPkRkB6zaSBEfdQYy+q djASwPlfRu/hLDQwZsL+8ZutxkQFzBmjfN2YksZ/MVshiunCHIda94QKogCFooehxmaOEsUv 90RRDBK5eTBsTV4ndHU5d3iWkKAHLVLNIumYDMu+EyQZmG9A4bo7H6NMMt/yy8eR9ch5LjPx RIiFrnYMR24w5UvTuEWjhB+qfky3evUnxrk5mik+/EfwH8i7XEi23loHQ4LXivCey0IeIMnO kBtLV15rIqHpYId3Cqul7O52P3Ugd838wgVUA==
  • Ironport-hdrordr: A9a23:/yBdhalighUD9a5K3RCwq1xYl1LpDfIP3DAbv31ZSRFFG/Fw5P re+sjzsiWE9Ar5OUtQ5OxoXZPrfZqyz+8W3WB8B8bYYOCighrTEGgA1/qG/9SDIVycygc178 4JGMhD4Z/LfCBHZK7BgTVQeOxQp+VvnprY/Ns2g00dKT2CqJsN0+66MGum+4FNKzWuzKBWKK ah
  • Ironport-phdr: A9a23:6TYyKB8bRZsRTf9uWdS2ngc9DxPPW53KNwIYoqAql6hJOvz6uci4Y QqEuqgm0QeBdL6YwswHotKei7rnV20E7MTJm1E5W7sIaSU4j94LlRcrGs+PBB6zBvfraysnA JYKDwc9rDm0PkdPBcnxeUDZrGGs4j4OABX/Mhd+KvjoFoLIgMm7yeC/94fdbghKhDexbrN/I AmrpgjNq8cahpdvJak2xhbVuHVDZv5YxXlvJVKdnhb84tm/8Zt++ClOuPwv6tBNX7zic6s3U bJXAjImM3so5MLwrhnMURGP5noHXWoIlBdDHhXI4wv7Xpf1tSv6q/Z91SyHNsD4Ubw4RTKv5 LptRRT1iikIKiQ5/XnUhMJ+g61Urg+vqQJxw4DUYo6aKPVwc7jBfd4YX2dNQtpdWiJDD466c oABD/ABPeFdr4TlulUBswCxChO3BOPzyz9HmHn21rA83eQmCwHG3QsgH84Tu3nTt9X1MKASU fq0zKnP0TXPdfZW2TPm5YjHdxAuu/CMXbZqfcXNzkkvEhrIg1ONooPqIz2bzP4Cs3SH7+V+T +KvjXYqpg5wrzWuwsoiiofEiIIJxl3H9St3wIU4KcO3RUJnb9OpFJReuS+EOodoTc4uXnxlt ig0xLAFt5O2cjYHxZQoyRPZdveJfY+I4hf5W+aQJzd1nH1leLOjhxay7Eiv0ffwWdWz0FZPq CdOj9rCtmgV2hDN9sSKTuFx80Sh1DqVyQzf9PxILVo7mKfZM5Iswb89moARvEnCBCP6hkv7g LGMekk5+uWl5OLqaaj8qJCGLY97kAT+P7wumsOhBeQ4NRADX22B9uS90L3v5EP4QK9XgvEvn KnUv5LXKd4Upq6+BA9V3YIj5AilAzi619QYmGELLFNDeB2Zk4jkI0/CLOz8APulgFmhkC1ny +7YMrH8GJnAIXrOnK/kfbln6k5czAQzzcpY55JRErwBO/bzWkjrtNzZCR85Lxa4w+PjCNhm2 YMeXHiAAq6dMK/IrVCI4ecvL/GKZIAOoDn9MeQq5+byjX8lnl8QZbSl0YMNaH+kBvRmP1mZY X30j9gdHmcFpA4+QPX3h12DSj5ce2uyX7kn5jwgE4KnDYLDRpi3j7Cb3Se7GIdWZmFcBVyWH 3fobdbMZ/BZYy2bJoptkycPSKO6Y44nzxCn8gHgmJR9Ke+BwSASvIn/1dVzr8nUjxA16XQgF MCayWCKUydxmksHQjY32OZ0pkkrmQTL6rRxn/ENTY8b3PhOSApvbfY0rsR/AtH2AUfae8uRD U2hWpOgCC0wSdQ4x5kPZVx8EpOslEOLxDKkVpkSkbHDH5ko6uTExXGkPMJ60XfByO8ki3EpR 8JOMSutgassvxPLCdvxml6C372vabxa2SfM8GmZym/bp0VfSg93TePAWVgQY0LXqZLy4UaRB 6S2B+EBNQ1MgdWHNrMMatDtigBeQ+z/Pd3Ff2+rs2K5BBLN2anVKYS2JCMS2yLSDEVCmAcWl ZqfHS45ACrp42fXDTg1UEnqf1up6u5m7nWyUk4zyQiOKUxnzbu8vBAP17SaTLsI07QItT1Ey X08FUuh39/QF9uLphZwNKRab9Qn5V5b1GXf/wVjN52kJqpmixYQaQNy90/p0hx2DM1Hn61I5 DsxyA1oKKWEllZFXzyd1JH0fLbQLyi6/RyibbLXxkCLyMyfqe8E7PU1rUmmvRn8TBJztSU6l YMPij3AucavbkJaS5/6X0cp+gIvorjbZnJ4/IbIzTh3NrHytDbe2tUvDe9jyxC6ft4ZPrnXc W26W8AcGcWqL/Qn3lazaRdRdvtW87Q1ON/gcv+u16uiPeImlzWjxzcigsg1wgeX+ixwR/Sdl Y0BxOuZ2RfBUjPUg1Kos8SxkodBL2J3fCL32W3vA4hfYbd3dIAABDK1IsG58d55gobkR39S8 FPL60ou4MayYlLSalX82VYVzkELuTm8njP+yTVokjYvp67Z3SrUwu2kegBVcmJMQWBjix/rL +3Wx5gBWE+ybg810h6hzUn/zqlf4q94KiHfTFxJcC7/M2x5GvHo5/zSPogVsspu7XofWf/0e V2AT7/hvxYWtkGrV3BTwjw2bXDiu5n0mQB7lHPIKX9yqHTDfsQjjRzb5dHaWbtQxm9cHHg+2 WSRXALleYX2rrD239/Zv+uzVnysTMhWeCjvlsabsTejoHZtGVu5luyyndvuFU471zX63p9kT 3atzl60b4/12qC9Ke8icFNvAQq28M59AIB5js01gLkf3HEbgtOe+n9NwgKReZ1LnLnzanYAX 2tB2NrS+gnixAtmKlqGwov4UjOWxc4rNLzYKisGnyk66c5NEqKd6rdJyDB0rlSPpgXUefFhn z0ZxKhm+DsAjuoOogZo0jSFD+VYAxxDJSK13UftjZj2vOBNaW2oa7T1yEdugYXrEuSZug8FE HfhJsV5QGkpv505ag6TliW0sN2secGMP4xP8EfPyFGZ0bATcNVox59ozWJmIT6v4yNjkrZhy 0Qohdbg5MCGMzk/ovz/WEIJcG2tIZtUoGmljL4CzJnMmdnzWMwwQHNTG8K4KJDgWDMK6aa4a 0DXSmB68jHDXuOBVQ6HtBU/9yKJSs/0cSHRfD5DlJ1jXEXPfRMExllFAHNi2MZ+T17PpoSpc V8ltGpJtxio910VkLIub16mDS/evFv6MG5qDsXPakMHtEcaoB6Eec2GsrApRn8epMbw6lfXb DTcPlUtbylBTEWAAxqL0qCGw97G/qDYA+O/K6CLerCSsalFUO/Ow5uz04xg9jLKN8OVP3AkA edpkkxEFWt0HcjUgVBtA2Qeij7NYsiHpRy95jw/r8ax9+7uUR7u4o3HAqVbMNFm8RS7yamZM OvYiCF8IDdenpQCoB2AgKAYx0IXgjpyeiOFFL0BsWvSVvuVlPMISREcbCx3OY1D6Kd9lghBN MjHi8/kg751ivlmbjUNHVflm8yveYkLOzTnbAKBVBvNbu7fY2GUkKSVKeumRLZdjftZrUi1s DefSAr4Oyib0iLuT1apOP1NiyeSOFpfvpu8e1BjEzuGLpquZxulPdtwlTBzz6czgyaALmobK zl7bwVIqpWf6CpZhrN0HGkLvR8HZaGU3j2U6eXVMMNcqfxwHiF9jP5X+lw/wrpRqT5eHbl7x HSUodlprFWr1OKIz3A0NXgG4iYOj4WNs0J4PKzf/ZQVQnfI8iUG6mCIAggLrd9oYjUOk6VZy 9yKi7iqbTkerJTb+swTA8WSI8WCYiJJ2fXBFzvdDQ9DRjmuZzi3b6l1n/Sb93nTpZ8/+MCEp Q==
  • Ironport-sdr: 646e0445_zEnC0dsWFS+jKhSG1cD8mgyR14QEbVv9kM8QR535yAA7yvc u55dsJS22M+Y+dw6GLuxpACfL5dX9nl7d577iMg==

Do you observe this behavior on all images or only a specific one?
For example, when using 420.inr do you have such an issue?

What could case such an over refinement (provided you are not using master version of CGAL) is if the mesher is trying to respect some
topological criteria (a very tiny handle for example, but that would be surprising in an image).

Best,

Sebastien.

On 5/24/23 14:14, Sebastian Bachmann ( via cgal-discuss Mailing List) wrote:
Hi everyone,

I have an image from which I want to create a coarse surface mesh.
The code I use for that is based upon the example mesh_3D_image_with_features.cpp, however it is adopted to then extract the surface and write an OFF file: https://gist.github.com/reox/0df1f446f0ff9d46e0688b60e8e858b8

As you can see, the Mesh_criteria is set to a really coarse mesh. The input image (attached) I use has a voxel size of around 0.5mm and represents a proximal femur.
Now, for some reason, I get extremely fine triangle elements in certain places (see attached image).

I find that peculiar, as the facet_distance is already quite large, and I thought that this might be the main issue here. I remembered, that also the relative_error_bound in the create_labeled_image_mesh_domain can have an influence, but reducing this rather creates a "spiky" mesh than removing those small elements. (Eventually when it is set to 0.1, the small elements vanish, but then the vertices are all over the place and not on the isosurface anymore...)

Is there anything I missed here, maybe did I use wrong parameters? Or is there another way how I could produce such coarse meshes?

Best,
Sebastian




Archive powered by MHonArc 2.6.19+.

Top of Page