.f.HdZddlZddlmZddlmZddlmZGddZy)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 dZ d Z dd Z dd Zdd Zy ) PolygonAreazArea of a geodesic polygonctj||\}}tj|}tj|}|dkDr|dcxkr|ks n|dkDr|dk(rdS|dkr|dcxkr|krdSdS)z/Count crossings of prime meridian for AddPoint.r)rAngDiff AngNormalize)lon1lon2lon12_s XC:\Users\user\Documents\project_loop\venv\Lib\site-packages\geographiclib/polygonarea.py_transitzPolygonArea._transit7s ||D$'HE1   T "D   T "D!1T!1!%dai @ dQ&6$&6"@&6A@ctj|d}tj|d}d|cxkrdkrnndndd|cxkrdkrdz Sdz S)z.Count crossings of prime meridian for AddEdge.g@rihr r remainder)r rs r_transitdirectzPolygonArea._transitdirectDsX >>$ &Dt~~dE/J?s?a?s?a+-"+-rc |j||dzr,|j|jdkrdnd|zdz |s|j|rR|j|dz kDr|j| ns|j| dz kr\|j|nJ|j|k\r|j| n$|jdkr|j|d|jzS)z)Reduce accumulator area to allowed range.r rr ) RemainderAddSumNegateareaarea0 crossingsreversesigns r _areareduceAzPolygonArea._areareduceAMs NN51} hhdhhj1n"5a79 DKKM eAg  5& 88:% ! 5 u  5& 88:> 5  rctj||}|dzr||dkrdnd|zdz z }|s|dz}|r%||dz kDr ||z}d|zS|| dz kr||z }d|zS||k\r ||z}d|zS|dkr||z }d|zS)z$Reduce double area to allowed range.r rr rrrrs r _areareduceBzPolygonArea._areareduceBds >>$ &D1} D1Hq"-a//d DBJD a   : E6!8    :    : !8   :rc\||_ ||_ dtjz|jz|_ t jt jzt jz|jrt jn t jt jzz|_ |jst|_t|_d|_ tj$|_ tj$|_ d|_tj$x|_|_y)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(TWWuyy(DJ?##h&8&88##$%)]]8>>==8#7#779DJ ==+-$-$DDH@DI4DI5DO"hh&DJrcd|_d|_|js|jj d|j j dt jx|_x|_ x|_ |_ y)zReset to empty polygon.rN) r8r;r+r6Setr7r,r9r<r=r:r )r>s rClearzPolygonArea.Clears\DHDO ==$--++A.16:hh>DJ>>di$)rc |jdk(r|x|_|_|x|_|_n|j j |j|j|||j\ }}}}}}}}}}|jj||jsN|jj||xjtj|j|z c_ ||_||_|xjdz c_y)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 xx1}""dj49""dj49)-)?)? 499c3 *4&aaAq!Q3 S! ]] # ;// 3??didiHHMHrc |jdk7r|jj|j|j|d||j \ }}}}}}}}}|j j||jsN|jj||xjtj|j|z c_ ||_||_|xjdz c_yy)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 xx1}(, (=(= 499c5!TZZ)9%ac1aAq# Q ]] # ;55diiEEdidi hh!mhrc  |jrtj}|jdkrd}|jsd}|j|fS|jr)|jj }|j|fS|j j|j|j|j|j|j\ }}}}}}}}}}|jj |}t|j}|j||j t"j%|j|jz} t"j'||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,r9r8r7rr*rDr:r r<r=r5rr6rr;rrr%r!) r>r#r$r perimeterrrGrHtempsumr"s rComputezPolygonArea.Computes7, }}TXXd xx!|i ]]3D XXy$ && }}$$((*i XXy$ &&'+zz'='= iiDJJ DJJ(@$AsAq!Q1a""&&s+I$--(G KK+"6"6tyy$**"MMI  # #GTZZ$+T 3D 88Y $$rc r|jrtj}|jdk(rd}|jsd}d|fS|jj }|jrdn|j j }|j}|jdz} |jrdgnddgD]} |jj| dk(r |jn|| dk(r |jn|| dk7r |jn|| dk7r |jn||j\ } } } } } } } } } } || z }|jr|| z }|tj!| dk(r |jn|| dk7r |jn|z }|jr| |fStj#||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,r9r8r7rr6r;r*rDr:r r<r=r5rrr'r!)r>rErFr#r$r rPrQr"r8irrGrHs r TestPointzPolygonArea.TestPoints }}TXXd xx1}i ]]3D  4 ""&&(I]]c (9(9(;GItxx!|]]qcA H)-)?)?!V 16diis1f #Q!Vtzz *&aaAq!Q33i ]]3[))qAv$))378Av3H H  H }} )T !!  # #GTZZ$+T 3D  4 rc ^|jdk(r!dtjtjfS|jdz}|jj |z}|j r||tjfS|j j }|j}|jj|j|j|d||j\ } } } } } } } } } || z }|tj|j| z }|jj| | |j |j"|j\ } } } } } } } } } } || z }|| z }|tj%| |j"z }tj'||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#s{ xx1} $(( "" ((Q,C""&&(1,I }} )TXX %%}}  "GI&*jj&;&; iiC4::'7#AsCAq!Q sNG ++DIIs;;I'+zz'='= 3 DJJ (4$AsAq!Q1a I sNG %%c4::66I  # #GTZZ$+T 3D  4 rN)F)FT)__name__ __module__ __qualname____doc__ staticmethodrrr%r'r?rBrIrNrRrUrWrrrr4sz" @ @--,,'<?..(%V& R# rr) r[r,geographiclib.geomathrgeographiclib.accumulatorrgeographiclib.geodesicrrr]rrras%\ &1+R R r