O i.PdZddlZddlmZddlmZddlmZGddZdS)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)GeodesicceZdZdZedZedZedZedZddZ dZ d Z d Z dd Z dd ZddZdS) PolygonAreazArea of a geodesic polygonctj||\}}tj|}tj|}|dkr|dcxkr|ksn|dkr|dkrdn|dkr|dcxkr|krnndndS)z/Count crossings of prime meridian for AddPoint.r)rAngDiff AngNormalize)lon1lon2lon12_s Q/srv/django_bis/venv311/lib/python3.11/site-packages/geographiclib/polygonarea.py_transitzPolygonArea._transit7s |D$''HE1  T " "D  T " "D!1!1!1!1T!1!1!1!1!%daii A dQ&6&6&6&6$&6&6&6&6&6""A@ctj|d}tj|d}d|cxkrdkrnndndd|cxkrdkrnndndz S)z.Count crossings of prime meridian for AddEdge.g@rihr r remainder)r rs r_transitdirectzPolygonArea._transitdirectDsm >$ & &Dt~dE/J/J????s?????aa????s?????aa+-rc~|||dzr5||dkrdnd|zdz |s||rd||dz kr|| n|| dz kr||n\||kr|| n-|dkr||d|zS)z)Reduce accumulator area to allowed range.r rr ) RemainderAddSumNegateareaarea0 crossingsreversesigns r _areareduceAzPolygonArea._areareduceAMs( NN51}: hhdhhjj1nn"5a7999 !DKKMMM  eAg   5& 88::% ! ! 5 u   5& 88::>> 5  rctj||}|dzr||dkrdnd|zdz z }|s|dz}|r||dz kr||z}n'|| dz kr||z }n||kr||z}n |dkr||z }d|zS)z$Reduce double area to allowed range.r rr rrrrs r _areareduceBzPolygonArea._areareduceBds >$ & &D1}0 D1HHqq"-a//d DBJD  a   E6!8        !88   :rFc||_ ||_ dtjz|jz|_ t jt jzt j z|jr t j nt j t j zz|_ |jst|_t|_d|_ tj|_ tj|_ d|_tjx|_|_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_c2r!rLATITUDE LONGITUDEDISTANCEEMPTYAREA LONG_UNROLL_maskr_areasum _perimetersumnumnanlat1r _crossings_lat0_lon0)selfr*r+s r__init__zPolygonArea.__init__{sDJ(DM(TWuy(DJ?#h&88#$%)]88>>=8#779DJ =7+--$-$DDH@DI4DI5DO"h&DJrcd|_d|_|js|jd|jdt jx|_x|_ x|_ |_ dS)zReset to empty polygon.rN) r8r;r+r6Setr7r,r9r<r=r:r )r>s rClearzPolygonArea.ClearsfDHDO =.$-++A...16:h>DJ>>di$)))rc |jdkr|x|_|_|x|_|_n|j|j|j|||j\ }}}}}}}}}}|j ||j sH|j ||xj t|j|z c_ ||_||_|xjdz c_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. rr N)r8r<r:r=r r* _GenInverser5r7rr+r6r;rr)r>latlonrs12S12s rAddPointzPolygonArea.AddPoints x1}}""dj49""dj499)-)?)? 49c3 *4*4&aaAq!Q3 S!!! ]@ # ;// 3???didiHHMHHHHrc |jdkr|j|j|j|d||j\ }}}}}}}}}|j||jsH|j ||xj t |j|z c_ ||_||_|xjdz c_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. rFr N) r8r* _GenDirectr:r r5r7rr+r6r;rr)r>azisrrErFrHs rAddEdgezPolygonArea.AddEdges x1}}(, (=(= 49c5!TZ)9)9%ac1aAq# Q ]F # ;55diEEEdidi hh!mhhhh}rTc ltj}|jdkrd}|jsd}|j||fS|jr#|j}|j||fS|j|j|j |j |j |j \ }}}}}}}}}}|j|}t|j}|||jt"|j |j z} 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,r9r8r+r7rr*rDr:r r<r=r5rr6rr;rrr%r!) r>r#r$r perimeterrrGrHtempsumr"s rComputezPolygonArea.Computes., 8D x!||i ]&3D Xy$ && }'$((**i Xy$ &&'+z'='= iDJ DJ(@(@$AsAq!Q1a"&&s++I$-((G KK+"6"6ty$*"M"MMI  # #GTZ$+T 3 3D 8Y $$rc tj}|jdkrd}|jsd}d||fS|j}|jrdn|j}|j}|jdz} |jrdgnddgD]} |j | dkr|j n|| dkr|j n|| dkr|j n|| dkr|j n||j\ } } } } } } } } } } || z }|js=|| z }|t| dkr|j n|| dkr|j n|z }|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,r9r8r+r7rr6r;r*rDr:r r<r=r5rrr'r!)r>rErFr#r$r rPrQr"r8irrGrHs r TestPointzPolygonArea.TestPoints 8D x1}}i ]&3D  4 "&&((I];cc (9(9(;(;GItx!|].qccA H H)-)?)?!VV 166diis1ff #Q!VVtzz **&aaAq!Q33i ]H3[))qAvv$))378Avv3HH H  }" )T !!  # #GTZ$+T 3 3D  4 rc |jdkrdtjtjfS|jdz}|j|z}|jr||tjfS|j}|j}|j |j |j |d||j \ } } } } } } } } } || z }|t|j | z }|j| | |j|j|j \ } } } } } } } } } } || z }|| z }|t| |jz }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) rr F)r8r,r9r7rr+r6r;r*rKr:r r5rrrDr<r=rr'r!)r>rLrMr#r$r8rPrQr"rrErFrHrGr s rTestEdgezPolygonArea.TestEdge#st x1}} $( "" (Q,C"&&((1,I }& )TX %%}  ""GI&*j&;&; iC4:'7'7#AsCAq!Q sNG ++DIs;;;I'+z'='= 3 DJ (4(4$AsAq!Q1a I sNG %%c4:666I  # #GTZ$+T 3 3D  4 rN)F)FT)__name__ __module__ __qualname____doc__ staticmethodrrr%r'r?rBrIrNrRrUrWrrrr4s"" @ @< @--<-<,<,''''<???..(%(%(%(%V& & & & R# # # # # # rr) r[r,geographiclib.geomathrgeographiclib.accumulatorrgeographiclib.geodesicrrr]rrras\ &&&&&&111111++++++R R R R R R R R R R r