Skip to Content.
Sympa Menu

cgal-discuss - Re: [cgal-discuss] Smooth surfaces at image boundary

Subject: CGAL users discussion list

List archive

Re: [cgal-discuss] Smooth surfaces at image boundary


Chronological Thread 
  • From: Sebastien Loriot <>
  • To:
  • Subject: Re: [cgal-discuss] Smooth surfaces at image boundary
  • Date: Mon, 9 May 2022 09:23:27 +0200
  • Authentication-results: mail3-smtp-sop.national.inria.fr; spf=None ; spf=Pass ; spf=None
  • Ironport-data: A9a23:2IR2maPSztX26H7vrR2WkcFynXyQoLVcMsEvi/4bfWQNrUoq3z0Pz TNMXWmAPvnca2HxctB/b9zg8RgBvsTTyddiGXM5pCpnJ55ogZqcVI7Bdi8cHAvLc5adFBo/h yk6QoOdRCzhZiaE/n9BCpC48T8kk/vgqoPUUIYoAAgoLeNfYHpn2EsLd9IR2NYy24DkW1PV4 7senuWGULOb824sWo4rw/nbwP9flKyaVOQw4zTSzdgS1LPvvyF94KA3fcldHFOkKmVgJdNWc s6YpF2PEsw1yD92Yj+tuu6TnkTn2dc+NyDW4pZdc/DKbhSvOkXe345jXMfwZ3u7hB24sOFr+ vF/8qa5WCA7YbHOxNkyD0NxRnQW0a1uoNcrIFC6uM2XilLDKj7in6koA0YxMokVvO1wBAmi9 9RCcGFLPk3F3rjmhu7nIgVvrpxLwM3DJ4cYoHBn0XfcCd4pRJnCR+PB4towMDIY35sTR6mCP 5ZxhTxHYh/yICAIal0tUI9vt9+J21PDQQJ6pwfAzUYwyzGLkFYZPKLWGNHac9jPScROlVuDv Urd7mHhC1cbMsaewHyL6BqRavTnmCr6XMcNEeT9+KI2xlKUwWMXBVsdUl7TTeSFZlCWav1uM X4wyG0UpLkV5hO3TIXAcEHiiSvR1vIDYOZ4H+o/4QCL76Pb5QeFG2QJJgKtjvR25KfaohR6h je0c8PV6S9H6+LKFCrMnluAhXbjZnhPdD5qiTosFFNdu7HeTJcPYgUjp+uP/Yawh9zxXC71m nWE8HR4iLIUgsoGka68+DgrYg5ARLCZH2bZBS2NBgpJCz+Vgqb7PuREDnCFt559wH6xFAXpg ZT9s5H2ABoyJZ+MjjeRZ+4GAauk4f2IWBWF3wM1R8R4p278qif8FWy13N2YDBc5WirjUW+5C HI/RSsMjHOuFCD3MvcoO9zZ5zoClPK+To6Nug/ogipmO8AtLmdrDQlhYkmf222FraTfufBXB HtvSu71VSxyIf0/klKeHr5BuZd2mH1W7T6MHfjTkkX/uZLDNSX9YepUYDOmMLplhIva+1692 4gEZ6O3J+B3CrKWjt//qt5Nczjn7BETWfjLliCgXrPdclc+RDF9Wpc8A9oJIuRYokicrc+Ql lnVZ6OS4AOXaaTvJVrYZ3Z9RqnoWJoj/3s3MTZ9b1mt0nknJ42o6f5HJZcweLAm8s1lzOJ1H 6FVIZXeXqwXR2SV4SkZYLn8sJdmK0amiAeICCyvP2oycptmcArW94K2ZQDo7iQPUnG6uJJm8 b2t3w/WW7QZQAFmAJqEYf6j1QLjsn0UmeY0VEzNe4EBdELp+YlsCirwkv5ne5FWeUufnmOXj l/EDw0ZqO/Bp54O3OPI3a3U/Z20F+ZeH1ZBGzaJ4LuzMx7c9DXxzIJFVtGOYmmBBm75/aOVZ dJVwev5B/sJkQsYqIF7Cbtqkfsz6taz9b9XygNoQCfCY1ixUOgyJ3CH2YxQs/QIyOMF/wSxX U2L959RPrDQYJHpF1sYJQwEaOWf1KFLxmODs6xteEiqtjVq+LenUFlJO0XegiJqKrYoYpgux v0suZJL5gHj2AAmNM2K0nJd+2iWdCdSVqwmstQFBdauhFZ3jF5FZpPYB2n955TWM4dANUwjI zm1gqvehuQDmhCTLSJrTXWdj/BAgZkuuQxRyANQLVq+nNeY1OQ82wdc8GhqQwlYpvmdPzmf5 oS224xJyaSyE/NAgcFCWyWzHlgECkHFvEP2zFQNmSvSSEzAuqkh6oEiEb7lwazb2zs0kvtnE HWwx2PsUDKsd8b0tsf3cVAwsOTtFLSd6SWb8P1K3K25810Sbj/sg6vobm0Nw/ci7QXdm2Wfz dRXEC1MhWEX+MLeT2DXy2VX6FjIdC25GQ==
  • Ironport-hdrordr: A9a23:3Tqi+6NIxpIHK8BcTuKjsMiBIKoaSvp037BZ7TEWdfUzSL3lqy nOpoVn6faQslwssR4b6LK90cW7MAjhHNtOkPMs1NSZLWzbUQmTXeJfBOLZqlWKexEWtNQtrZ uIG5IeNDSaNykcsS+V2njALz9t+qjizElqv4vjJrVWID2Cp5sO0+6xMGimLnE=
  • Ironport-phdr: A9a23:d683IR194WXdaFmvsmDOZw4yDhhOgF0UFjAc5pdvsb9SaKPrp82kY BaEo649xwaQA82bs6sC17CM9fi4GCQp2tWoiDg6aptCVhsI2409vjcLJ4q7M3D9N+PgdCcgH c5PBxdP9nC/NlVJSo6lPwWB6nK94iQPFRrhKAF7Ovr6GpLIj8Swyuu+54Dfbx9HiTajb75+N hu7oAbeusULnIdvK7s6xh/VrndVYehbyn1kKUyJkxrg+su8+YNo/jhNtf4m68NOS7jxcb4iT bxfAjQmMmQ169PuuBLeUwaB5WYSX3sPnBZQDAfL8B/1XpHqsivnreV9wzWVPdf3Tb8vRzuv6 bpgRQLyhycGMz4593zXitB1galGrh+tuwBzzojJa4yTKfFwfL7SfckCSGRCUMheWCJODJ6iY YsBD+QOP/pXoYbmqlsSrxazHxWgCP/hxzJKgHL9wK000/4mEQHDxAEuAswBv2/Oo9XuKKcSU P66zK/VxjvFaPNWwzb96I/Uchs8pvyMRrZwcc3LyUkuEAPFkk+cqYnhPz6O2eQNtnKU7+tkV e61l2EnrARxryGpy8wxhYbHmpgbxUrY9SVl3ok1P9u4RVZ4bNCqDpddtC6UOpVoTs4jXmxkp ik3xqEbtZOnfCYHyZoqygDRZvGGcIWG7Q7vWeeRLDtmmn9rdr2yigq2/ES+zOD3S8e60FFPr iVfk9nMsGgA2ADc6sedTvtx5ECh2SyA1w3V9+pKIlg0mLLFJ5I9xrM8jJkevETZEiPrhkn7j 7Wae0or9+Sw6OnqY6/qq52dOoNulA3zPLoiltaiDek4PQUDWXWQ9/6m273550L5Ra1Hjv0on andt5DXPcEbqbS4Aw9Ry4oi5RO/Ayq/3NQWknQKI0hJeB2Aj4juNFHOJO73Ae2jjFSrlTdn3 /HGPrv/DZXRNnXPjqvtcLJn50NfyAc/181T64xXB70bL//+VVf9tNnCAR84Nwy0zfznCNJ41 o4GQmKPBqiZMLnTsV+P+uIjOeyMZIsOtTbyLvgq/f/ujXsjlVABeqmp2IMbaHajEft+OUWWf WLsgssdEWcNpgcyUfTmiEeeXj5Le3ayQ6U86yknB4KpF4jDQpmhj6GA3CegApJWe3tGCkuXH Hfzd4SEXu8MZziILs9glDwET7mhRJU72RGgrg+pg4Zge+Hb8ylduZP438Vu/MXSkwsz/Hp6F ZezyWaIGlp5lGoTWz472uhbplZ8zUvLhbN8hOZZEsAV4vdhXQIzNJqaxOt/XYOhEjndd8uEH Q71Cu6tBis8G4pZK74mZk98H47nlRXfx2+xBKdTkbWXBZsy+6aa3n7rJs87xWyVnLI5gQwAR c1CfXajmrY57xLaUpXNlF+YkLrscKA00yvE9WPFxm2L7wlDSAAlaazeRjgEY1fO69Hw50fMV birXK8jNRFAztLEL69iZdjgjFEAT/Dma5zFe2zkvWC2CF6Tw6+UKornf2JIxCLGFE0NiBwe5 16DPAk6Qzik+ifQUGAoGlXobEfht+J5rRtXV2cSyAeHJw1k3ruxoVsOgOCEDugUxvQCsTsgr DN9GBC82cjXApyOvVgpeqIUet477FpdsACR/wVgIpytKbxjjV8CYkx2uU3pzRB+FoRHl4Ajs noryAN4La/Q3klGcnuU2pX5O7ufLWeXnljncKrRwFDZzJCT/o8A7f05rxPouwToXksu/nN70 sVEhmOG78aCBw4TXJTtF0cvokIi9veKP29kvdOSiSE/VMv8+iXP0N8oGuY/nxOpftMEdbiBC Be3CMoCQc6nNO0tnVGtKBMCJuFbsqAubKbEP7OL3rCmOOF4kXeol2NCtcpm1kWW9i1gDOvM9 5kAyvCcmACAUn2v6TXp+tCygo1CaTwISyCkzS/+BYlNIKh2VYkOAGaqZcaww58t4vylE24d/ 1mlCVQc3cavchfHdF3x0zpb0kEPqGCmky+1p9BtuwkgtbHXnCnHwuC5MQECJnYOX256y1HlP YmzidkeGkmudQkg0hW/tw72wK1SpaI3KGe2Iw8AZCzxNWBlTu21spKNZsdO7NUjtiAfXOmnY F+cQ6LwuFNAi3KlTzYYnWprMW328pzi1wR3kmecMGp+oB+7MYlryBHT6cadDf9d0zwaRTVp3 DzeB1yyJd6srp2fk5bOtPz7Vnr0DMUCN3m2i9ra5G3nuT4PY1X3hf24l9z5HBJv1Cb604IvT iDUtFPnZZGt0a2mMOVhd00uBVnm6sM8FJst9+l4zJwWx3UegY2YuHQdlmKme85f3rj/a2ZLQ DojzNvc4Qyj00pmZCHspcqxRjCGz81tasPvKHgS3Tg86NwMDaO847lNnC8zqV2951G0A7A1j nIWzv0g72QfiucCtV82zymTNbsVGFFRIS3mkxnbp8D7tqhcY3yjNKShzEcr1870F6mM+0sPP RSxModnByJ76d9zdU7BwGGmoJ+xY8HeNJoSrkHGyEqG1rkNbshtybxSwnA7cWPl4S97l6hh1 kcohM/i+tDAcjQInurxAwYEZGOrIZpLoHe1y/4ZxJ7e3pjzTMs/XG9XDd24FbTwV2hK/fX/a 1TRSntl9jHCSOCZRUjGuCIE5zrOC8z5aC3RfSNEi40kHF7EegRemFxGBW1q2MdmSUb6gpSmK h4x5yhNtAen8V0Vm74ub1+nFT6BwWXgIjYsFMrFdEsQvlwEvh2Fd5TZt707HjkErMf48krQe irCNl4OVSZQCwSFHwyxZODwo4Oboq7DXKzmaKKfBNfG4fpXU/PCrX62+q1h+TvEdsCGP305S uY+xlIGRndhXcLQhzQITSUT0SPLdc+S4hmmqGVxqYik/fLnVRiKh8PHAqZOMdhp5xG9gLuSf ++WiiFjLD9E15QKjXbWwbkb1VQWhmlgbT6oWbgHsCfMSurXlMo1R1YDbDhvMcJT868m9gxEO Mqeldask7Ak0rg6DFBKUVGnkcasJIQLL2y7KFLbFROLObCBdlipi4n8ZaKxT6EVjf0B7UXh/ 2bGVRa5bnLfy2qMNVjnK+xHgSCFMQYLvYi8dk0oEm3/VJf9bQX9NtZrjDowyLlyh3XQNGdaP yIvFiEF5rCW8y5che1yXmJb6X8wZ/KAnDyY6PWeL5I+vv5iAyAynOVfqidfqfMd/GRfSfp5l TGH5MZpuE2jm/KTxyBPVRNPrnNUgdvOsxk+f6re8ZZEVDDP+xdHvgDyQ1wa4tBiDNPooaVZz NPCwbnyJDl1+NXR5cIABsLQJapv1VIuNBPoHHjfCw5XFVZD1Enajk1Z1e6Xrzib8sR8pZ/rl 54DDLRcUQ5tfhv1IktgFd0GZpxwW2F9+YM=
  • Ironport-sdr: Q3JVxgmgar7rSxzTKjD/Qz2sKxRO6COa8oRGt1shFBL6ND4z1z1AgygglJv3gy+aF5sCt05QLw TNDe/SBZR4cgmDdeHqRjKk8/Skc8ZE5/OG9P6DYdLPI7B/Yb822AtmAWuO4ozhYtcW76LAQO30 2ujVTic08vQYnIbQFLUXasfTwzy4FGgVpdY7RRNg1RJYMz2uAeOEhCfOilCnPwuz7wtT9b2I03 85Sv9rv9clnsw+mpysHAhHQg1zsVX52pgU4jj3g/n53sJ4kDdhKg5TuAVoI4tHEcixqXjDvSrH T2H5Mk9W9N4xCCcu4RSlAhPI

Hello,

could you share a minimal example showing the problem?

Thanks,

Sebastien.

On 5/7/22 13:37, Sebastian Bachmann ( via cgal-discuss Mailing List) wrote:
Hello,
No one has an idea? Is it not somehow possible to restrict nodes onto a surface? I.e. if I have the polyline of the surface's boundary, can I specify that all nodes must be on that surface?

For now I have written a function to go over all nodes and check if they are near such a boundary face and then "pull" the nodes onto that surface. However, there are certain edge cases were this then yields collapsed triangles etc. Thus, solving this during meshing would be a great enhancement.

Best,
Sebastian

On 25.04.2022 07:16, Sebastian Bachmann ( via cgal-discuss Mailing List) wrote:
Hello,

I created an image meshing tool, similar to the one given here in this example code: https://doc.cgal.org/latest/Mesh_3/Mesh_3_2mesh_3D_image_with_features_8cpp-example.html with the only difference, that I extract the surface elements and not write the full 3D mesh.
I noticed, that on the image boundary, the vertices would be not exactly on the face but roughly +- the relative error bound given in create_labeled_image_mesh_domain. For examples, nodes on the xy plane at z=0 are not on z=0 but for example -0.01 to 0.01. While the protected polylines are exact and have always z=0.
I created for example a mesh from the 420.inr example with the same meshing criteria as given in the example and can see that behaviour there.

While you can control the maximal deviation using the relative error bound, there will always be some deviation. And if I set the relative error bound too low, it seems the meshing never finishes.
I wonder now, if it is possible to constraint the mesh (similar to the 1d polylines), such that the vertices for facets on the image boundary, are always exactly on the image boundary?

Best,
Sebastian





Archive powered by MHonArc 2.6.19+.

Top of Page