o ýòTiã.ã@sBdZddlZddlmZddlmZddlmZGdd„dƒZdS)a‡Define the :class:`~geographiclib.polygonarea.PolygonArea` class The constructor initializes a empty polygon. The available methods are * :meth:`~geographiclib.polygonarea.PolygonArea.Clear` reset the polygon * :meth:`~geographiclib.polygonarea.PolygonArea.AddPoint` add a vertex to the polygon * :meth:`~geographiclib.polygonarea.PolygonArea.AddEdge` add an edge to the polygon * :meth:`~geographiclib.polygonarea.PolygonArea.Compute` compute the properties of the polygon * :meth:`~geographiclib.polygonarea.PolygonArea.TestPoint` compute the properties of the polygon with a tentative additional vertex * :meth:`~geographiclib.polygonarea.PolygonArea.TestEdge` compute the properties of the polygon with a tentative additional edge The public attributes for this class are * :attr:`~geographiclib.polygonarea.PolygonArea.earth` :attr:`~geographiclib.polygonarea.PolygonArea.polyline` :attr:`~geographiclib.polygonarea.PolygonArea.area0` :attr:`~geographiclib.polygonarea.PolygonArea.num` :attr:`~geographiclib.polygonarea.PolygonArea.lat1` :attr:`~geographiclib.polygonarea.PolygonArea.lon1` éN)ÚMath)Ú Accumulator)ÚGeodesicc@s€eZdZdZedd„ƒZedd„ƒZedd„ƒZedd „ƒZdd d „Z d d„Z dd„Z dd„Z ddd„Z ddd„Zddd„ZdS)Ú PolygonAreazArea of a geodesic polygoncCszt ||¡\}}t |¡}t |¡}|dkr+|dkr |ks)n|dkr+|dkr+dS|dkr;|dkr:|kr;dSdS)z/Count crossings of prime meridian for AddPoint.rééÿÿÿÿ)rZAngDiffZ AngNormalize)Úlon1Úlon2Zlon12Ú_©r úY/var/www/html/evchargy.com/venv/lib/python3.10/site-packages/geographiclib/polygonarea.pyÚ_transit7s  ÿþþzPolygonArea._transitcCsVt |d¡}t |d¡}d|krdkrnndndd|kr'dkr(dSdS)z.Count crossings of prime meridian for AddEdge.g€†@rihr©rÚ remainder)rr r r r Ú_transitdirectDs ÿÿzPolygonArea._transitdirectcCs¼| |¡|d@r| | ¡dkrdnd|d¡|s | ¡|r@| ¡|dkr1| | ¡n'| ¡| dkr?| |¡n| ¡|krM| | ¡n | ¡dkrX| |¡d| ¡S)z)Reduce accumulator area to allowed range.rrréç)Z RemainderÚAddÚSumZNegate©ÚareaÚarea0Ú crossingsÚreverseÚsignr r r Ú _areareduceAMs "  €    zPolygonArea._areareduceAcCsªt ||¡}|d@r||dkrdnd|d7}|s|d9}|r=||dkr.||8}d|S|| dkr9||7}d|S||krI||8}d|S|dkrQ||7}d|S)z$Reduce double area to allowed range.rrrrrrrr r r Ú _areareduceBds"    øûýzPolygonArea._areareduceBFcCsš||_ ||_ dtj|j|_ tjtjBtj B|jr tj ntj tj BB|_ |js/tƒ|_tƒ|_d|_ tj|_ tj|_ d|_tj|_|_dS)zãConstruct a PolygonArea object :param earth: a :class:`~geographiclib.geodesic.Geodesic` object :param polyline: if true, treat object as a polyline instead of a polygon Initially the polygon has no vertices. érN)ÚearthÚpolylineÚmathÚpiZ_c2rrZLATITUDEZ LONGITUDEZDISTANCEZEMPTYZAREAZ LONG_UNROLLÚ_maskrÚ_areasumÚ _perimetersumÚnumÚnanÚlat1rÚ _crossingsÚ_lat0Ú_lon0)Úselfrrr r r Ú__init__{s,  ÿ  ýzPolygonArea.__init__cCsHd|_d|_|js|j d¡|j d¡tj|_|_ |_ |_ dS)zReset to empty polygon.rN) r%r(rr#ÚSetr$r r&r)r*r'r)r+r r r ÚClear™s  zPolygonArea.Clearc Cs¦|jdkr||_|_||_|_n8|j |j|j|||j¡\ }}}}}}}}}}|j  |¡|j sD|j   |¡|j t  |j|¡7_ ||_||_|jd7_dS)zÕAdd the next vertex to the polygon :param lat: the latitude of the point in degrees :param lon: the longitude of the point in degrees This adds an edge from the current vertex to the new vertex. rrN)r%r)r'r*rrÚ _GenInverser"r$rrr#r(rr )r+ÚlatÚlonr Ús12ÚS12r r r ÚAddPoint¡s ÿ  zPolygonArea.AddPointc Cs|jdkrF|j |j|j|d||j¡\ }}}}}}}}}|j |¡|js7|j  |¡|j t   |j|¡7_ ||_||_|jd7_dSdS)zëAdd the next edge to the polygon :param azi: the azimuth at the current the point in degrees :param s: the length of the edge in meters This specifies the new vertex in terms of the edge from the current vertex. rFrN) r%rÚ _GenDirectr'rr"r$rrr#r(rr)r+ÚaziÚsr r0r1r3r r r ÚAddEdge¸s ÿ  ÷zPolygonArea.AddEdgeTc CsÔ|jrtj}|jdkrd}|jsd}|j||fS|jr&|j ¡}|j||fS|j |j|j |j |j |j ¡\ }}}}}}}}}}|j |¡}t |jƒ}| |¡|jt |j |j ¡} t ||j| ||¡}|j||fS)a‹Compute the properties of the polygon :param reverse: if true then clockwise (instead of counter-clockwise) traversal counts as a positive area :param sign: if true then return a signed result for the area if the polygon is traversed in the "wrong" direction instead of returning the area for the rest of the earth :return: a tuple of number, perimeter (meters), area (meters^2) Arbitrarily complex polygons are allowed. In the case of self-intersecting polygons the area is accumulated "algebraically", e.g., the areas of the 2 loops in a figure-8 polygon will partially cancel. If the object is a polygon (and not a polyline), the perimeter includes the length of a final edge connecting the current point to the initial point. If the object is a polyline, then area is nan. More points can be added to the polygon after this call. rr)rr r&r%r$rrr/r'rr)r*r"rr#rr(rr rr) r+rrrÚ perimeterr r2r3Útempsumrr r r ÚComputeÏs&      ÿ    ÿ zPolygonArea.Computec CsH|jrtj}|jdkrd}|jsd}d||fS|j ¡}|jr!dn|j ¡}|j}|jd} |jr4dgnddgD]T} |j  | dkrD|j n|| dkrL|j n|| dkrT|j n|| dkr\|j n||j¡\ } } } } } } } } } } || 7}|jsŒ|| 7}|t | dkr€|j n|| dkrˆ|j n|¡7}q8|jr•| ||fSt ||j|||¡}| ||fS)aCompute the properties for a tentative additional vertex :param lat: the latitude of the point in degrees :param lon: the longitude of the point in degrees :param reverse: if true then clockwise (instead of counter-clockwise) traversal counts as a positive area :param sign: if true then return a signed result for the area if the polygon is traversed in the "wrong" direction instead of returning the area for the rest of the earth :return: a tuple of number, perimeter (meters), area (meters^2) rrr)rr r&r%r$rr#r(rr/r'rr)r*r"rr rr)r+r0r1rrrr9r:rr%Úir r2r3r r r Ú TestPointús6      ýÿ€  ÿ zPolygonArea.TestPointc Cs|jdkr dtjtjfS|jd}|j ¡|}|jr!||tjfS|j ¡}|j}|j  |j |j |d||j ¡\ } } } } } } } } } || 7}|t  |j | ¡7}|j | | |j|j|j ¡\ } } } } } } } } } } || 7}|| 7}|t  | |j¡7}t  ||j|||¡}|||fS)aCompute the properties for a tentative additional edge :param azi: the azimuth at the current the point in degrees :param s: the length of the edge in meters :param reverse: if true then clockwise (instead of counter-clockwise) traversal counts as a positive area :param sign: if true then return a signed result for the area if the polygon is traversed in the "wrong" direction instead of returning the area for the rest of the earth :return: a tuple of number, perimeter (meters), area (meters^2) rrF)r%r r&r$rrr#r(rr5r'rr"rrr/r)r*r rr)r+r6r7rrr%r9r:rr r0r1r3r2rr r r ÚTestEdge#s.    ÿÿ ÿ zPolygonArea.TestEdgeN)F)FT)Ú__name__Ú __module__Ú __qualname__Ú__doc__Ú staticmethodr rrrr,r.r4r8r;r=r>r r r r r4s"       +)r) rBr Zgeographiclib.geomathrZgeographiclib.accumulatorrZgeographiclib.geodesicrrr r r r Ús .