U …7hã.ã@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 polygoncCs~t ||¡\}}t |¡}t |¡}|dkrV|dkr@|ksRn|dkrV|dkrVdS|dkrz|dkrr|krznndSdS)z/Count crossings of prime meridian for AddPoint.rééÿÿÿÿ)rZAngDiffZ AngNormalize)Úlon1Úlon2Zlon12Ú_©r úS/var/www/formularioweb/env/lib/python3.8/site-packages/geographiclib/polygonarea.pyÚ_transit7s  ÿÿ&þzPolygonArea._transitcCsXt |d¡}t |d¡}d|kr,dkr4nndndd|krJdkrRnndndS)z.Count crossings of prime meridian for AddEdge.g€†@rihr©rÚ remainder)rr r r r Ú_transitdirectDs  ÿzPolygonArea._transitdirectcCs¼| |¡|d@r4| | ¡dkr&dnd|d¡|s@| ¡|r€| ¡|dkrb| | ¡q°| ¡| dkr°| |¡n0| ¡|krš| | ¡n| ¡dkr°| |¡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@r0||dkr"dnd|d7}|s<|d9}|rn||dkrV||8}q|| dkr||7}n"||kr€||8}n|dkr||7}d|S)z$Reduce double area to allowed range.rrrrrrrr r r Ú _areareduceBds     zPolygonArea._areareduceBFcCsŽ||_||_dtj|j|_tjtjBtj B|jr:tj n tj tj BB|_ |jsXtƒ|_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$||_|_||_|_np|j |j|j|||j¡\ }}}}}}}}}}|j  |¡|j sˆ|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Œ|jdkrˆ|j |j|j|d||j¡\ }}}}}}}}}|j |¡|jsn|j  |¡|j t   |j|¡7_ ||_||_|jd7_dS)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Ô|jr tj}|jdkr0d}|js$d}|j||fS|jrL|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Ïs4    ÿ    ÿzPolygonArea.Computec CsN|jr tj}|jdkr.d}|js$d}d||fS|j ¡}|jrBdn|j ¡}|j}|jd} |jrhdgnddgD]¬} |j  | dkrˆ|j n|| dkr˜|j n|| dkr¨|j n|| dkr¸|j n||j¡\ } } } } } } } } } } || 7}|jsp|| 7}|t | dkr|j n|| dkr|j n|¡7}qp|jr0| ||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ús@     ýÿ  ÿzPolygonArea.TestPointc Cs|jdkrdtjtjfS|jd}|j ¡|}|jrB||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#sB    ÿÿ ÿzPolygonArea.TestEdgeN)F)FT)FT)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 .