o JÉfTã@sDddlmZmZddlmZddlmZGdd„dƒZd dd„Zd S) é)ÚheappopÚheappush)ÚTopologicalError)ÚPointc@sPeZdZdZdd„Zdd„Zdd„Zdd „Zd d „Zd d „Z dd„Z dd„Z dS)ÚCellaNA `Cell`'s centroid property is a potential solution to finding the pole of inaccessibility for a given polygon. Rich comparison operators are used for sorting `Cell` objects in a priority queue based on the potential maximum distance of any theoretical point within a cell to a given polygon's exterior boundary. cCs>||_||_||_t||ƒ|_| |¡|_|j|d|_dS)NgÍ;fž ö?)ÚxÚyÚhrÚcentroidÚ_distÚdistanceÚ max_distance)Úselfrrr Úpolygon©rúh/home/deployuser/azure_apps/autowriter/venv/lib/python3.10/site-packages/shapely/algorithms/polylabel.pyÚ__init__s   z Cell.__init__cCs |j|jkS©N©r ©rÚotherrrrÚ__lt__ó z Cell.__lt__cCs |j|jkSrrrrrrÚ__le__rz Cell.__le__cCs |j|jkSrrrrrrÚ__eq__"rz Cell.__eq__cCs |j|jkSrrrrrrÚ__ne__%rz Cell.__ne__cCs |j|jkSrrrrrrÚ__gt__(rz Cell.__gt__cCs |j|jkSrrrrrrÚ__ge__+rz Cell.__ge__cCsF| |j¡}|j |j¡}|jD] }t||j |¡ƒ}q|r |S| S)z¨Signed distance from Cell centroid to polygon outline. The returned value is negative if the point is outside of the polygon exterior boundary. )Úcontainsr r ÚexteriorÚ interiorsÚmin)rrÚinsider Úinteriorrrrr .s  z Cell._distN) Ú__name__Ú __module__Ú __qualname__Ú__doc__rrrrrrrr rrrrrs  rçð?cCs¤|jstdƒ‚|j\}}}}||}||}t||ƒ}|d} g} |jjd\} } t| | d|ƒ} t||d||dd|ƒ}|j| jkrG|} |} | |krq|} | |krit| t| | | | | |ƒƒ| |7} | |ksS| |7} | |ksM| rÏt | ƒ}|j| jkr|} |j | j|krˆqq|j d} t| t|j | |j | | |ƒƒt| t|j | |j | | |ƒƒt| t|j | |j | | |ƒƒt| t|j | |j | | |ƒƒ| ss| jS)aÁFinds pole of inaccessibility for a given polygon. Based on Vladimir Agafonkin's https://github.com/mapbox/polylabel Parameters ---------- polygon : shapely.geometry.Polygon tolerance : int or float, optional `tolerance` represents the highest resolution in units of the input geometry that will be considered for a solution. (default value is 1.0). Returns ------- shapely.geometry.Point A point representing the pole of inaccessibility for the given input polygon. Raises ------ shapely.errors.TopologicalError If the input polygon is not a valid geometry. Example ------- >>> from shapely import LineString >>> polygon = LineString([(0, 0), (50, 200), (100, 100), (20, 50), ... (-100, -20), (-150, -200)]).buffer(100) >>> polylabel(polygon, tolerance=10).wkt 'POINT (59.35615556364569 121.83919629746435)' zInvalid polygong@ré)Úis_validrÚboundsr!r Úcoordsrr rrr r rr)rÚ toleranceÚminxÚminyÚmaxxÚmaxyÚwidthÚheightÚ cell_sizer Ú cell_queuerrÚ best_cellÚ bbox_cellÚcellrrrÚ polylabel<sF  þû      ïr9N)r() ÚheapqrrÚshapely.errorsrÚshapely.geometryrrr9rrrrÚs   5