Skip to main content
Log in

An algorithm for maximum inscribed circle based on Voronoi diagrams and geometrical properties

  • Original Article
  • Published:
Journal of Geographical Systems Aims and scope Submit manuscript

Abstract

The aim of this study is to formulate an algorithm for the calculation of maximum inscribed circle (MIC) that can be placed within a polygon and to implement it by using free and open source software (FOSS) for GIS. MIC is used in a wide range of fields, ranging from cartography, planning, agriculture, forestry and geology to medicine, biology, astronomy, security, and engineering applications. Due to the complexity of the problem, there is no single and simple algorithm for the computation of MIC for arbitrary polygons. The algorithm developed in this study (MICGIS) for the computation of MIC can be applied to both convex and concave polygons represented in vector data format. MICGIS makes use of the Voronoi diagrams and geometrical properties by benefiting from the solutions proposed for the special cases of Apollonius’ Problem. Thanks to the employment of Voronoi diagrams and FOSS for GIS, MICGIS also works successfully for polygons with holes. For the implementation of MICGIS, FOSS libraries written in Java are used. What is evident from the various runs of the script produced on the base of MICGIS for a set of arbitrary polygons is that it is both faster and more accurate in finding MIC compared with the alternative algorithms and software.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Fig. 1
Fig. 2
Fig. 3
Fig. 4
Fig. 5
Fig. 6

Similar content being viewed by others

References

  • Agafonkin V (2016) A new algorithm for finding a visual center of a polygon. Available: https://blog.mapbox.com/a-new-algorithm-for-finding-a-visual-center-of-a-polygon-7c77e6492fbc Accessed 15 Dec 2019

  • Alonso R, Ojeda J, Hitschfeld N, Hervías C, Campusano LE (2018) Delaunay based algorithm for finding polygonal voids in planar point sets. Astron Comput 22:48–62

    Google Scholar 

  • Amenta N, Bern M, Eppstein D (1998) The Crust and the beta-skeleton: combinatorial curve reconstruction. Gr Models Image Process 60:125–135

    Google Scholar 

  • Amenta N, Choi S, Kolluri RK (2001) The power crust, unions of balls, and the medial axis transform. Comput Geom 19(2–3):127–153

    Google Scholar 

  • Arcelli C, Frucci M (1992) Reversible skeletonization by (5, 7, 11)-erosion. In: Arcelli C, Cordella LP, di Sanniti Baja G (eds) Visual form, analysis and recognition. Springer, Boston, pp 21–28

    Google Scholar 

  • Aschoff T, Spiecker H (2004) Algorithms for the automatic detection of trees in laser scanner data. In: Thies M, Koch B, Spiecker H, Weinacker H (eds) WG VIII/2 laser–scanners for forest and landscape assessment, 3–6 October 2004 Freiburg, Germany. Albert Ludwigs University, Freiburg, pp 71–75

    Google Scholar 

  • Attali D, Montanvert A (1997) Computing and simplifying 2D and 3D continuous skeletons. Comput Vis Image Underst 67(3):261–273

    Google Scholar 

  • Birdal T (2011) Maximum inscribed circle using distance transform. Available: https://www.mathworks.com/matlabcentral/fileexchange/30805-maximum-inscribed-circle-using-distance-transform Accessed 29 Jan 2019

  • Birdal T (2013) Maximum inscribed circle using voronoi diagram. Available: https://www.mathworks.com/matlabcentral/fileexchange/32543-maximum-inscribed-circle-using-voronoi-diagram Accessed 15 Dec 2019

  • Blum H (1967) A transformation for extracting new descriptors of shape. In: Walthen-Dunn W (ed) Models for the perception of speech and visual form. MIT Press, Boston, pp 362–380

    Google Scholar 

  • Boissonnat JD, Wormser C, Yvinec M (2006) Curved voronoi diagrams. In: Boissonnat JD, Teillaud M (eds) Effective computational geometry for curves and surfaces. Springer, Berlin, pp 67–116

    Google Scholar 

  • Borgefors G (1993) Centres of maximal discs in the 5-7-11 distance transform. In: Proceedings of the Scandinavian conference on image analysis, vol 1, p 105

  • Brezina T, Graser A, Leth U (2017) Geometric methods for estimating representative sidewalk widths applied to Vienna’s streetscape surfaces database. J Geogr Syst 19(2):157–174

    Google Scholar 

  • Brune C, Vogt S, Peiper C, Brinker K, Trzewik J (2018) A new universal pore measurement and clustering approach for surgical meshes. bioRxiv: 446450. Available: http://biorxiv.org/content/early/2018/10/17/446450.abstract

  • Burri O, Guiet R (2016) Max inscribed circles. Available: https://imagej.net/Max_Inscribed_Circles Accessed 30 Jan 2019

  • Butsch A (2016) Districting problems—new geometrically motivated approaches. Karlsruhe: Doctoral dissertation, Karlsruhe Institut für Technologie

  • Cavic L, Šileryte R, Beirão JN (2017) 3D-informed convex spaces: the automated generation of convex representation for open public space analysis. In: Heitor et al. (eds) Proceedings of the 11th international space syntax symposium, Instituto Superior Técnico, Departamento de Engenharia Civil, Arquitetura e Georrecursos, Portugal, pp 2433–2447

  • Cheng T, Li P, Zhu S (2012) An algorithm for jammer localization in wireless sensor networks. In: Barolli L, Enokido T, Xhafa F, Takizawa M (eds) 26th IEEE international conference on advanced information networking and applications, 26–29 Mar 2012 Fukuoka-shi, Japan. IEEE Computer Society Conference Publishing Services (CPS), Tokyo, pp 724–731

  • Chin F, Snoeyink J, Wang CA (1999) Finding the medial axis of a simple polygon in linear time. Discre Comput Geom 21(3):405–420

    Google Scholar 

  • Court NA (1961) The problem of Apollonius. Math Teach 54(6):444–452

    Google Scholar 

  • Davis M (2012) JTS topology suite. Available: http://www.tsusiatsoftware.net/jts/main.html Accessed 30 Jan 2019

  • Davis LG, Bean DW, Nyers AJ, Brauner DR (2015) GLiMR: a GIS-based method for the geometric morphometric analysis of artifacts. Lithic Technol 40(3):199–217

    Google Scholar 

  • ESRI (Environmental Systems Research Institute) (2015) ArcMap Version 10.3.1. 380 New York Street, Redlands, 92373-8100, USA

  • Fanwu M, Chunguang X, Haiming L, Juan H, Dingguo X (2011) A quick algorithm of maximum inscribed circle method for roundness evaluation. In: 2011 international conference on system science, engineering design and manufacturing informatization, Guiyang, The Institute of Electrical and Electronics Engineers, 22–23 Oct 2011, vol 2, pp 348–351

  • Flowmap (2013) Educational Flowmap Release 7.4.2. Available: http://flowmap.geo.uu.nl/ Accessed 30 Jan 2019

  • Garcia-Castellanos D, Lombardo U (2007) Poles of inaccessibility: a calculation algorithm for the remotest places on Earth. Scott Geogr J 123(3):227–233

    Google Scholar 

  • Gass SI, Witzgall C, Harary HH (1998) Fitting circles and spheres to coordinate measuring machine data. Int J Flex Manuf Syst 10(1):5–25

    Google Scholar 

  • ET GeoWizards (2018) ET GeoWizards Version 12.1. Available: https://www.ian-ko.com/ETGeoWizards.html Accessed 10 Feb 2018

  • Hadaś E (2015) An overview of methods for tree geometric parameter estimation from ALS data in the context of their application for agricultural trees. Acta Sci Pol Geodesia et Descriptio Terr 14(3–4):5–28

    Google Scholar 

  • Jones DR (2014) A fully general, exact algorithm for nesting irregular shapes. J Global Optim 59(2):367–404

    Google Scholar 

  • Karimipour F, Ghandehari M (2013) Voronoi-based medial axis approximation from samples: issues and solutions. In: Gavrilova ML, Tan CJK, Kalantari B (eds) Transactions on computational science XX, vol 8110. Lecture notes in computer science. Springer, Berlin

    Google Scholar 

  • Kimmel R, Shaked D, Kiryati N, Bruckstein AM (1995) Skeletonization via distance maps and level sets. Comput Vis Image Underst 62(3):382–391

    Google Scholar 

  • Koch B, Heyder U, Weinacker H (2006) Detection of individual tree crowns in airborne Lidar Data. Photogramm Eng Remote Sens 72(4):357–363

    Google Scholar 

  • Lam L, Lee SW, Suen CY (1992) Thinning methodologies—a comprehensive survey. IEEE Trans Pattern Anal Mach Intell 14(9):869–885

    Google Scholar 

  • Lee DT (1982) Medial axis transformation of a planar shape. IEEE Trans Pattern Anal Mach Intell 4:363–369

    Google Scholar 

  • Li W, Goodchild MF, Church R (2013) An efficient measure of compactness for two-dimensional shapes and its application in regionalization problems. Int J Geogr Inf Sci 27(6):1227–1250

    Google Scholar 

  • Li YH, Feng QG, Hou ZS (2017) Classification of Diospyros Lotus seeds based on MLP neural network. In: 2nd international conference on software, multimedia and communication engineering (SMCE 2017), DEStech Publications, Pennsylvania, Shanghai, 23–24 Apr 2017, pp 498–503

  • Liu F, Xu G, Liang L, Zhang Q, Liu D (2016) Minimum circumscribed circle and maximum inscribed circle of roundness deviation evaluation with intersecting chord method. IEEE Trans Instrum Meas 65(12):2787–2796

    Google Scholar 

  • MacEachren AM (1985) Compactness of geographic shape: comparison and evaluation of measures. Geogr Ann Ser B Hum Geogr 67(1):53–67

    Google Scholar 

  • Masset L (2002) Cvoronoi version 1.0.00. Computes the best approaching circle of a set of 2D points. Available: https://uk.mathworks.com/matlabcentral/fileexchange/2794-cvoronoi Accessed 30 Jan 2019

  • OpenJUMP (2019) Version 1.14.1 release rev.6147. Available: http://www.openjump.org/ Accessed 2 Aug 2019

  • Petrík M, Kováč J, Kaťuch P, Bednarčíková L, Hudák R, Živčák J (2009) Roundness: determining the reference circle for MCCI and MICI system. In: Tyšler M, Maňka J, Witkovský V (eds) Proceedings of the 7th international conference measurement 2009, Smolenice, Slovakia. Komprint, Bratislava, 20–23 May 2009, pp 352–355

  • Powers MC (1953) A new roundness scale for sedimentary particles. J Sediment Res 23(2):117–119

    Google Scholar 

  • QGIS (2019) QGIS version 3.8.3-Zanzibar. Available: https://qgis.org/tr/site/ Accessed 15 Aug 2019

  • Ramamurthy R, Farouki RT (1999) Voronoi diagram and medial axis algorithm for planar domains with curved boundaries-I: theoretical foundations. J Comput Appl Math 102(1):119–141

    Google Scholar 

  • Reddy JM, Turkiyyah GM (1995) Computation of 3D skeletons using a generalized Delaunay triangulation technique. Comput Aided Des 27(9):677–694

    Google Scholar 

  • SAGA (2018) SAGA Version 6.3.0. Available: http://www.saga-gis.org/en/index.html Accessed 25 Feb 2018

  • Saxon J (2020) Reviving legislative avenues for gerrymandering reform with a flexible, automated tool. Political Analysis, pp 1–23

  • Shen Z, Yu X, Sheng Y, Li J, Luo J (2015) A fast algorithm to estimate the deepest points of lakes for regional lake registration. PLoS ONE 10(12):e0144700. https://doi.org/10.1371/journal.pone.0144700

    Google Scholar 

  • Siddiqi K, Bouix S, Tannenbaum A, Zucker SW (2002) Hamilton-jacobi skeletons. Int J Comput Vision 48(3):215–231

    Google Scholar 

  • Sun S (2016) Symbolize map distortion with inscribed circles in polygons. Int J Cartogr 2(2):166–185

    Google Scholar 

  • Sun Y, Che R (2003) Novel method for solving maximum inscribed circle. Opt Precis Eng 11(2):181–187

    Google Scholar 

  • Taghanaki SA, Liu Y, Miles B, Hamarneh G (2017) Geometry-based pectoral muscle segmentation from mlo mammogram views. IEEE Trans Biomed Eng 64(11):2662–2671

    Google Scholar 

  • Talbot H, Vincent LM (1992) Euclidean skeletons and conditional bisectors. In: Visual communications and image processing’ 92, vol 1818, pp 862–876

  • Tauscher S, Neumann K (2016) A displacement method for maps showing dense sets of points of interest. In: Gartner G, Jobst M, Huang H (eds) Progress in cartography. Lecture notes in geoinformation and cartography. Springer, Cham, pp 3–16. https://doi.org/10.1007/978-3-319-19602-2_1

    Google Scholar 

  • Thies J (2018) LayoutEditor Free Version. Available: http://www.layouteditor.net/ Accessed 10 Jan 2018

  • Tsygankov D, Bilancia CG, Vitriol EA, Hahn KM, Peifer M, Elston TC (2014) Cell Geo: a computational platform for the analysis of shape changes in cells with complex geometries. J Cell Biol 204(3):443–460. https://doi.org/10.1083/jcb.201306067

    Google Scholar 

  • Vogt S, Ritter C, Trzewik J, Brinker K (2017) Measuring, clustering and classifying pores of surgical meshes with an ImageJ plug-in. In: 2nd YRA MedTech symposium, Hochschule Ruhr-West, Germany. DuEPublico online, 8–9 June 2017, pp 5–8. http://dx.doi.org/10.17185/duepublico/43984

  • Wadell H (1932) Volume, shape and roundness of rock particles. J Geol 40(5):443–451

    Google Scholar 

  • Wadell H (1933) Sphericity and roundness of rock particles. J Geol 41(3):310–331

    Google Scholar 

  • Wadell H (1935) Volume, shape, and roundness of quartz particles. J Geol 43(3):250–280

    Google Scholar 

  • Walz U, Schumacher U (2005) Landscape fragmentation in the Free State of Saxony and the surrounding border areas. In: Hřebíček J, Ráček J (eds) Proceedings of the 19th international conference on informatics for environmental protection networking environmental information, Brno, Czech Republic. Masaryk University, Brno, 7–9 September 2005, pp 754–758

  • Weinacker H, Koch B, Heyder U, Weinacker R (2004) Development of filtering, segmentation and modelling modules for LIDAR and multispectral data as a fundament of an automatic forest inventory system. In: Thies M, Koch B, Spiecker H, Weinacker H (eds) WG VIII/2 laser-scanners for forest and landscape assessment. Albert Ludwigs University, Freiburg, pp 50–55

    Google Scholar 

  • Zheng J, Hryciw RD (2015) Traditional soil particle sphericity, roundness and surface roughness by computational geometry. Géotechnique 65(6):494–506

    Google Scholar 

Download references

Acknowledgements

Birsen F. Biçgel is thanked for her assistance in preparation of technical drawings of the first benchmark dataset (named Regular_Geopolys.shp) composed of 28 regular geometric shapes. Whereas, the second benchmark dataset (named Complex_Geopolys.shp) is composed of 491 complex polygon features representing a real-world dataset adapted from Güler et al. (2012) [https://doi.org/10.1016/j.jhydrol.2011.11.021]. The authors would like to thank them for sharing this dataset with us. The authors are also grateful to anonymous referees for their useful comments, suggestions and constructive critiques. This manuscript in draft was presented at the 7th International Conference on Cartography & GIS (7 ICC&GIS 2018) held by Bulgarian Cartographic Association in Sozopol, Bulgaria (June 18–23, 2018). The authors would like to acknowledge the participants who made insightful comments and suggestions at the presentation. This work received no specific grant from any funding agency in the public, commercial, or not-for-profit sectors.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Burak Beyhan.

Ethics declarations

Conflict of interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary material 1 (RAR 201 kb)

Appendix

Appendix

See Fig. 7.

Fig. 7
figure 7

MICGIS algorithm developed in the study

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Beyhan, B., Güler, C. & Tağa, H. An algorithm for maximum inscribed circle based on Voronoi diagrams and geometrical properties. J Geogr Syst 22, 391–418 (2020). https://doi.org/10.1007/s10109-020-00325-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s10109-020-00325-3

Keywords

JEL Classification

Navigation