An efficient algorithm of recursive calculation of the probability of intersection of an unlimited number of vector polygons is created. The developed algorithm can be used for statistical analysis of any data presented as polygonal objects. Such objects are used in many Earth sciences, including geology, hydrometeorology, and glaciology, for example, to describe the edges of ice, the boundary of glaciers, permafrost, snow cover, and the position of certain isolines (isotherms, isobars, etc.). In addition to the algorithm for calculating the statistical characteristics of vector polygons, the article describes the adaptation of the method for calculating spatial density and interquantile zones for linear and point objects. As an example, the article describes the results of analyzing the position of landfast ice (polygonal objects), icebergs (point objects), and navigation routes of ships in the Arctic according to AIS data (linear objects). © The Author(s), under exclusive license to Springer Nature Switzerland AG 2025.