Skip to Content.
Sympa Menu

cgal-discuss - Re: [cgal-discuss] Compute radius of inscribed circle

Subject: CGAL users discussion list

List archive

Re: [cgal-discuss] Compute radius of inscribed circle


Chronological Thread 
  • From: Sana Habeeb <>
  • To:
  • Subject: Re: [cgal-discuss] Compute radius of inscribed circle
  • Date: Tue, 12 Mar 2024 05:09:51 -0400
  • Authentication-results: mail3-smtp-sop.national.inria.fr; spf=None ; spf=Pass ; spf=None
  • Ironport-data: A9a23:6HHIuKsZuq5acUspv8S0W7GIe+fnVPZaMUV32f8akzHdYApBsoF/q tZmKTjTOv6OMGCjfdolPt/k8kIHuJeAnYMwQVFlq3pjFyoWgMeUXt7xwmXYb3rDdJWbJK5Ex 5xDMYeYdJhcolv0/ErF3m3J9CEkvU2wbuOgTrSCYEidfCc8IA85kxVvhuUltYBhhNm9Emult Mj7yyHlEAbNNwVcbCRMtMpvlDs15K6u4GhB5ARkDRx2lAa2e0c9XMp3yZ6ZdCOQrrl8RoaSW +vFxbelyWLVlz9F5gSNz94X2mVTKlLjFVDmZkh+A8BOsTAezsAG6ZvXAdJHAathZ5plqPgqo DlFncTYpQ7EpcQgksxFO/VTO3kW0aGrZNYriJVw2CCe5xSuTpfi/xlhJEY8Mq425uQqODlP+ vU0FwoncgKB2P3jldpXSsE07igiBMziPYdat305iD+EVLApRpfMR6iM7thdtNsyrpoWTLCOO oxDMWopMU2ojx5nYj/7DLolleGhwHD7cjQeqFufqKc9y2fWxQ11lrPqNbI5f/TaHJUExx7A/ Aoq+UzXEBYDD+e+wgGF/yiPqbPCoAH4e4UdQejQGvlC2wDKnjNCVnX6T2CTqvawjguyWslUN lcP0jE/qLA7sk2tVNj0GROiyENopTYZUttUVuA2sUSDl/OS7AGeCWwJCDVGbbTKqfPaWxQt3 FCqsc/rWAU3sbKcZU2n3KuGr2ytbH19wXA5WQcISg4M4t/GqY41jw7SQtsLLEJTpo2lcd0X6 2DaxBXSl4kuYdg3O7JXFG0rbhqpr5nNCwM5v0DZAzvj4QR+a4qoIYev7DA3DMqszq7IETFtX 1BdxKByCdzi67nTyURhp81TQNmUCw6tamG0vLKWN8BJG86R03CiZ5tMxzp1OV1kNM0JERewP xaL4FMItMEOZCD6BUOSX25XI5R6pUQHPYS1Ps04kvIXOfCdiSfeoH80Px/PgQgBbmB3yf1nZ v93jvpA/V5BVP0/k2voLwvs+bAswS86yCvSQ5u9pylLIpLPDEN5vYwtaQPUBshgtP3siFyMr 753aZHWoz0BC7aWSneMoeYuwaUidydT6Wbe8JwJKIZu42NORAkcNhMm6e94I9I4wfkEybagE 7PUchYw9WcTTEbvcW2iAk2PopuyNXqmhStjZ3J+DkXiwHU5f4ek4YEWcpZ9L/Ft9/VuwbQwB 7MJctmJSKYHADnW2SUvXb+kpqxbdTOvmV2vOQiha2MBZJJOfVHC1eLlWQrNzxMwKBSLm/Ewm YD96TODc6E/H1xjKO30dMOQy0iAuClBueBqAGrNDNphWGTt14lILSbOoOc9CJwOI0+bxx+x9 QWfMTEHr8bj/q4399jog/ifjoGLSuFRIGtTL1P5352XaxbI3zOE6pBSdcq1ZhbhbXPQ1IT+Q PRK3tf+HeYinl0Xg7FjEr1u870y1+Hvq5Be0A5gOnfBNHavNZ9NPViE2ttppIRW57oEpzazZ F2DyuNaNZqNJsnhNlwbfyghT+ab0MAriivg1us0LGr69R1I0uK+C2sKBCa1iQtZMLdRG6Emy 714uMcptiqOuiBzOdOC1i1p52CAK0IbaJoet7YYPZTKjzQ6wVQTcL3eDS7LuKu0UetuCXVzA DGoh/vlvY9+l27iaHs4EEbf0dVN3aouvA94935cBlCrtOecuNoJ8kxwyxoVQD5R7C176MNoG 20yN0RKNaSEpDhpo85YXlGTIQJKBTzH20n912o2kHb9Sm+2XFfsN0w4A/6Gp2oCwlJffx9a3 bCW83nkWjDUZ/PM3jM+dEpmiv77R/lzy1HmtOW4Of+aRr8WTCHAgKC8QUYp8T7cHtIXlknLg cJI7dRAQ/T3GgBIqpJqFrTA86obTS61AVBrQNZj2fgsNn7dcjTj4gq+ARm9Ve0VLsOb7HLiL dJlI/9OcBGM1CyuiDQ/LoxUKp9WmM8Z3vYzSonJF0Un7YTG9iFItajO/BfQnGUoGtVisfgsI 7PrKg6tLDaiumt2qUTs8u9/JWuKUfsVblbd3ce00tkzObAtjeVOSXw2g5yI5yi7EQ0/5B+Fn hLxV4mPxcxY9IlctY/NEKJCOgaKFe3ORNm4qACdj/kebPfkE9v/iAcOm1y2YyVUJeQwXvp0p 5Sssfn2/l/0g7Itd1/nmLy6Tqxv2eivbrAGLPCtPH1+mA2cUvTN+DoGwXizcrZSoeNe5+6mZ gq2U9SxftgrQOVgxGVZRixdMhQFAYHlR/3EiQLkiNrUETkb8wjMDO3/xE/TdWsBKxM5YczvO DH7q9OFx45+rr0VIDQmGvs/IZtzAGG7aJscb9er6AWpVDi5sGij5InnuwErswzQK3++F83/3 5LJaz7+eDm2u4DK1NtpiJNzjDJGEEdChfQMQWxF9+5UkzyaCEs0HdYZO7gCCbBWlXXW/7P8b zfvcmAjKHvcWRJpTBbC2+ngDzyvXrE2BtTEJzIXphLeL2/8AY6bG7Ju+xtx+3o8KHOp0OijL spY4XHqeAS4xpZyX+sI+/inmqFdy+jHwm4Ts1XI+yAo783y3Z1RvJCgIOZMacADO8TElUGOK GRsAG4ZGQe0Tkn+FcsmcHlQcP3cUPUD0B1wBRpjAv6G02lY8AGE4PL6MuD3lLYEaazm4ZYQE GjvSTLlD3++gxQuVGhAhz7tqaBxAPOPWMO9KccPgOHUc76YsgwaAi/JocbDoAzONuKS/5MxW wRAO0QDOXk=
  • Ironport-hdrordr: A9a23:HFIVRqtnA+5Mh3YvE2iALSY47skDT9V00zEX/kB9WHVpm62j5r mTdZEgvyMc5wxhPU3I9erwWpVoBEmslqKdgrNxAV7BZniDhILAFugLhrcKgQeBJ8SUzJ876U 4PSdkZNDQyNzRHZATBjTVQ3+xO/DBPys6Vuds=
  • Ironport-phdr: A9a23:SK+ItRzmu/zXa/DXCzKawlBlVkEcU1XcAAcZ59Idhq5Udez7ptK+Z h2Zuq8m1QWZFcWDsrQY0bqQ6/ihEUU7or+/81k6M6ZwHycfjssXmwFySOWkMmbcaMDQUiohA c5ZX0Vk9XzoeWJcGcL5ekGA6ibqtW1aFRrwLxd6KfroEYDOkcu3y/qy+5rOaAlUmTaxe7x/I AuyoAnLtMQan4RuJrsvxhbKv3BEf/hayX5yKV+cgRrx5d288IJ//yhVpvks69NOXaLmcqs3S rBVEC4oP2cp6sP1qBLNVxGP5nwSUmUXlhpHHQ3I5wzkU5nyryX3qPNz1DGVMsPqQ780Xy+i7 7pwRx/zlCgHLT85/3rJhcF2kalWvQiupx17w47TfYGVKP9zdb7TcN8GWWZMWNtaWipcCY2+c oQPFfIMMulWr4b/p1UAoxiwCxSyCuzz0TJHnGP60Lcg3ug9DQ3L3gotFM8OvnTOq9X1Mb8fX /2vw6nO1zrIcvFY1i3m54fWbB8uuuyHULVrccXLz0kgDRnJgU+UqYzkJTyZzOQNvHKA4OV+S +2gkHUqqw50oje1x8csjpPFiZ4SylDB7Ch0xps+KtKkRkBhe9GkDIdQuD+AN4twWs4sTH1lt Tg1xLACuZC2YDYGxZs6yhPca/GJfImF7x39WeuNJTp0mHFodryxihuv8kWt1/PwW8ep3FtJq idJj8fAu3YQ3BLd7ciHT+Fy/kan2TuX2ADT6/1EIUAplaXAMZIu3rkwlp8LvUTCGC/2gkT2g 7WMekUq4Oil7fnsbLb+ppKEKYN4lgXzPr4tl8G/G+g0LBQCUmuB9em8ybHu+1DyTatOjv0ri anZrI7VJd4GqK64Ag5azJ4u5wq5Ajy7ytoXh2MHI0hAeB+fj4jmJVXOIPfgAPe6mVuskTNry +nYMbL8H5nBN3bDnbf8cbpn5E5czw0zzd9b551KEL0OPPXzWkrpuNzZCB82LRC0zv75BNlh0 o4SQ2GCD6+DPK/MrVOF5fgjLumRaIMNvTbyMfkl5/rgjX8jnl8deLGk3YANZ3C+BPtqOV+WY X7tgtgbC2sKvhAxQ/f0h1KcSjFefHmyX6cm6TE6DIKqF5vMRoeogLCZ2ie0BYVZZnpaBVCUD Xfoa4KEVu8RZC6KOM9ujiQEVaS9S48mzRyhqAD6xKBjLubN5yIYtInj28Ro5+3IjhE/7id0D sSY02GVVW54hGIIRzks3KB+u0Nx0FmD0bIry8FeDsFZsvZAUwMmMszVyeBgBtnpHAzIdM2IY FmrRdCrHSsgQNs639gUcgB2HNDxsxbYwiD/A6MJj6fZQ9sv46fE1j7wIdx8wjDIzu46nlw+S 4xOM2OhwaVw/gyWC4/SmFiCjPWXc78B1g7R8WPWyGOPuF1fAgd2Xb/IVGtZaELct9XR6U7FS rKyE6UpOwBdztSTbKBNb47Ui05bTtf/PdCLY36tg3ziQlGT17aUZczrfX8c1WPTEg8fggUL9 DGHMwY5QSyuqmabAD11Hk/0eBDR9/Jjok+2XlNhzx2Wd1YzkP2u6xsNjLqdTekS1/QKomA6u jBsFRG82dzRTNGPrg4kcKRHas4m+wR70jfSuAV5e5ChNKt/nUU2cgJtvkqo2Q8kJJ9HlJ0Rp XUjhCp7I6+V2hsVbDqe1Nb5PrzebGn79h2iYYbZ31jf1JCd/aJZu6dwkEnqoAz8ThlqyH5gy dQAixN0h73PBQsWCtfqV1ovsgN9v/fcazU84IXd0TttN7O1u3nMwYFhH/MrnzCnedoXK6aYD EnqCcRPHMGoIaojn1mgKBINNu1f/oY7Osqnc72N36v4dP14kmeehH9cqJt4zlrK8iN9TuDS2 JNQ2Pie102OUzb1yleouM36k6hLYDgTGiy0zi23TJVJaPhUeoAGQXyrP9Xxxth6gMv1XGVE8 Ve4G14c8MqgeB7XYlCkmAMNjAIYpnuonSb+xDtx+90whoyY2iGGg+HrdR5cf3VOWHEnllD0Z 469k9EdWkGsKQkvjhqso0jgleBdo+xkIm/fTF0tHWC+Jnx+Uqa2qruJYtJeoJIuvyJNVe2gY FecArfjqhoe2innEiNQ3jc+Pz2tv5z4mVR9hgf/ZD5opXzQP8h0xB6Z5NXWSPNb9jUDTSh8z zLQAxn0PtWk+8mVi4aWqvq3BAfDHtVYdSjmy5/FtTPuvzU7R03i2arrypu6SFtfs2ezzdRhW CTWoQypZ4Dq0///Kud7ZgxzA0e67cNmG4Z4m492hZcK2HFciI/GmBhP2Wr1L9hf3rrzKXQXQ jteicXc5gGj0UliJzSJzoj0W3+1zc5oZt38aWQTkHFYjYgCGOKP4bpIkDEg6EGxqwuXaPx8l 3EdwPEj5XQyjOQAuQ5rxSKYSON3fwEQLWnnkBKG6Mq7paNcaTO0cLS+40F5mMioELCIpgw0t G/RQp45Bmcw68x+NAmJy3jv8sT/f8GWa9sPtxqSmhOGju5PKZt3mOBYzSZgPGv8uzUixYtZx VR11Jez+oGDKmwr/KuzDxNSHjLwbsIXvDrqiO5Sk92X0IamAph6Um9TDd24ELTySGJU6K2vP h3rcnV0snqBHLvDAQKToFxrqX7CCdHjNn2aImUY0cQ3QRCcIEJFhwVHFD4+n5M/CkWr3Jm7K BY/tm1XvAao7EYVmYcKf1HlX2zSpRmlcGIxQZmbd19N6x1aolzSKYqY5/5yGCdR+tugqhaMI yqVfVctbylBV0qaClTkJrTr68PH9r3SHeO4Jr3HZbSH7+xXXfqFz7qg14Jn+3CHMcDFbRwAR 7Urn1FOW3x0AZGTgDwJQmoekyPEKc6Srhqx8AV4q8m+9LLgXweltu7tQ/NCdN5o/R6xm6KKM eWd0T14JThv3ZQJ3XbUybIb0Q1am2R0ejKqC7hFqT/VQfeahPpMFxBCIXAWVoMA/+cm0wJKI 8Kelt7lyusykKsuE1kcHV35xpPyOIpTcjn7bg+YQh7MbujOJCWXkZ+rJ/nnEvsJ0r0S71rp6 FP5WwfiJmjRyWevDkj1d7kK1GbBZFRfoN3vLEgrUzSyCoK+LEX8aoc/jCVqk+JuwCqWczdNa 34kNBody9/YpSJA3qcgRyoYtCcjda/c3H/Hp+jAdsRP6astW3sr0bIcuDNgkvNU9H0WHaMu3 nKD8pg25Qnhy7fqqHIvUQIS+GwT1cTb4AM7Y/Wfrt4ZBj7F5E5ftzzOTUla4YI0UJu3/PkBg tnXyPCpcWkEqYmFu5BGQZCTcZPiUjJpJxPtHHS84BItazmtOCmfgkVclKrX7XiJtt0hrYCqn pMSS7hdXVhzF/UAC00jEsZQaJFwFigpl7KWlqtqrTK3sQXRScNGv5vGSuPaAPPhLyychKVFY B1AyK3xLIAaPIn2k0J4bVwykIPPEkvWFddDx08pJhczu1lI+WNiQ3cb3kvkbkag7CZWG6Pt2 BExjQR6bKIm8zKtq1Y7K1zWpTchxUk8ndK25FLZODX1Laq2QcRXE3+u7xl3Ys69GV4lK1Drz ikGfH/eSrlcjqVtbzVugQ7Y4t5UHOJECLZDeFkWzO2WYPMh1RJdrD+mzAlJ/7iga9MqmQ00f Jqrt38F1RhkaYt/PqXWLexAyllZwKKHtCKl3Mg+xQYfIwAG92bYK0tq8AQYc6IrISal5Lkm8 QuZhz5KY3QBTdIvq/Nus04xYqGOknum3LlEJUS8ceeYKunK3gqI3d7NSVQ22EQSkkBD9rUjy sYvfX2fUEU3xaeQHRAEXSInAQ5QZstWsnPUeHTW2Q0s6Zd8PoH4G+yxCOHT7eAbhUWrGAtvF IMJvJxp9naE30TRLMOhJ7kAm0xF2Q==
  • Ironport-sdr: 65f01bec_vqekN+iSA195vN5cZzWuyc+Fw9BUddazC54vCB0YNdfTg1k Lj6cPKZzu+g3fIz3dC1hmWkC2QYgyT5RjuaaR7w==

Hi all,

How can I unsubscribe from this list?

Thank you.

Kind regards,
Sana

On Tue, Mar 12, 2024 at 4:16 AM, Andreas Fabri <> wrote:

You are right.

I fixed the example:  https://github.com/afabri/cgal/blob/SVD-inscribed_circle-GF/Segment_Delaunay_graph_2/examples/Segment_Delaunay_graph_2/sdg-inscribed-circle.cpp

You can also see some input files in the directory data/

But again, there is room for improvement.

Best,

Andreas

On 3/12/2024 12:54 AM, Bob Bill ( via cgal-discuss Mailing List) wrote:
Hello,

I tried your snippet (thanks a lot!) on the following polygon (it's the square [-1,0]x[-1,0], I've attached a mwe):


boundary_polygon.push_back(Point_2(-1., -1.));
boundary_polygon.push_back(Point_2(-0.5, -1));
boundary_polygon.push_back(Point_2(0., -1.));
boundary_polygon.push_back(Point_2(0., -0.5));
boundary_polygon.push_back(Point_2(0., 0.));
boundary_polygon.push_back(Point_2(-0.5, 0.));
boundary_polygon.push_back(Point_2(-1., 0));
boundary_polygon.push_back(Point_2(-1., -0.5));

and it gives as result 0 instead of 0.5. After debugging, I realised the check:

if (fit->vertex(0)->site().is_segment() &&
fit->vertex(1)->site().is_segment() &&
fit->vertex(2)->site().is_segment()))

is always false for this example, but I can't see the reason why it fails. Is there an obvious fix for this?


Best,
Bob

Il lunedì 11 marzo 2024 alle ore 15:37:01 CET, Andreas Fabri ha scritto:


Hello,

The idea is here :  https://gist.github.com/afabri/da5f00ff59eaadad3f00a7873670b897

There is place for improvement as the test if the point is inside the polygon is each time O(n log n).

This should evolve in a function provided by the library, but maybe this is better than nothing for you.

Best,

Andreas


On 3/9/2024 12:58 PM, Bob Bill ( via cgal-discuss Mailing List) wrote:
Hi,

I've been trying for a while to compute **the radius** of the inscribed circle of a general polygonal shape in CGAL. As suggested in another post, I should compute the segment Voronoi diagram. 

Unfortunately, I am not able to understand how to compute the "segment voronoi diagram" in CGAL. Is  Segment_Delaunay_graph_2 the right class?
 
Moreover, once I have the Voronoi diagram, how can I compute such a radius? Is there any code snippet already available for such a task? 


Best regards,
Bob

--
You are currently subscribed to cgal-discuss.
To unsubscribe or access the archives, go to
https://sympa.inria.fr/sympa/info/cgal-discuss

-- 
Andreas Fabri, PhD
Chief Officer, GeometryFactory
Editor, The CGAL Project

--
You are currently subscribed to cgal-discuss.
To unsubscribe or access the archives, go to
https://sympa.inria.fr/sympa/info/cgal-discuss


--
You are currently subscribed to cgal-discuss.
To unsubscribe or access the archives, go to
https://sympa.inria.fr/sympa/info/cgal-discuss

-- 
Andreas Fabri, PhD
Chief Officer, GeometryFactory
Editor, The CGAL Project

--
You are currently subscribed to cgal-discuss.
To unsubscribe or access the archives, go to
https://sympa.inria.fr/sympa/info/cgal-discuss




Archive powered by MHonArc 2.6.19+.

Top of Page