o 2c.@sBdZddlZddlmZddlmZddlmZGdddZdS)aDefine 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@seZdZdZeddZeddZeddZedd Zdd d Z d dZ ddZ ddZ dddZ dddZdddZdS) 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)rAngDiff AngNormalize)lon1lon2lon12_rO/var/www/html/gps/gps/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@rihrr remainder)r r rrr_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) RemainderAddSumNegateareaarea0 crossingsreversesignrrr _areareduceAMs "      zPolygonArea._areareduceAcCst||}|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.rrrrrrrrrr _areareduceBds"    zPolygonArea._areareduceBFcCs||_ ||_ dtj|j|_ tjtjBtj B|jr tj ntj tj BB|_ |js/t|_t|_d|_ tj|_ tj|_ d|_tj|_|_dS)zConstruct 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)earthpolylinemathpi_c2rrLATITUDE LONGITUDEDISTANCEEMPTYAREA LONG_UNROLL_maskr_areasum _perimetersumnumnanlat1r _crossings_lat0_lon0)selfr#r$rrr__init__{s,    zPolygonArea.__init__cCsHd|_d|_|js|jd|jdtj|_|_ |_ |_ dS)zReset to empty polygon.rN) r1r4r$r/Setr0r%r2r5r6r3r )r7rrrClears  zPolygonArea.Clearc Cs|jdkr||_|_||_|_n8|j|j|j|||j\ }}}}}}}}}}|j ||j sD|j ||j t |j|7_ ||_||_|jd7_dS)zAdd 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)r1r5r3r6r r# _GenInverser.r0rr$r/r4rr)r7latlonr s12S12rrrAddPoints   zPolygonArea.AddPointc Cs|jdkrF|j|j|j|d||j\ }}}}}}}}}|j||js7|j ||j t |j|7_ ||_||_|jd7_dSdS)zAdd 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) r1r# _GenDirectr3r r.r0rr$r/r4rr)r7azisr r<r=r?rrrAddEdges   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)aCompute 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)r$r%r2r1r0rr#r;r3r r5r6r.rr/rr4rrr r) r7rrr perimeterr r>r?tempsumrrrrComputes&           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)r$r%r2r1r0rr/r4r#r;r3r r5r6r.rrr!r)r7r<r=rrrrErFrr1ir r>r?rrr TestPoints6         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 | |j7}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)r1r%r2r0rr$r/r4r#rAr3r r.rrr;r5r6rr!r)r7rBrCrrr1rErFrr r<r=r?r>rrrrTestEdge#s.      zPolygonArea.TestEdgeN)F)FT)__name__ __module__ __qualname____doc__ staticmethodrrr r!r8r:r@rDrGrIrJrrrrr4s"       +)r) rNr%geographiclib.geomathrgeographiclib.accumulatorrgeographiclib.geodesicrrrrrrs .