ó ø`9Sc@s,ddlZdgZddd„ƒYZdS(iÿÿÿÿNt GeodeticPointcBs\eZdZdZdZdd„Zd„Zd„Zd„Zd„Z d „Z d „Z RS( s@Geodetic coordinates of a point with respect to WGS84 ellipsoid.g@¦TXAg—Ä?XAgcCs||_||_||_dS(ssParameters: - 'latitide': geodetic latitude, in degrees (angle between the equatorial plane and a normal to the reference ellipsoid passing through the point, positive direction northward). - 'longitude': geodetic longitude, in degrees (angle between the meridional plane of the reference meridian and a normal to the reference ellipsoid passing through the point, positive direction eastward). - 'altitude': geodetic altitude, in meters (normal distance to the point from the reference ellipsoid, positive direction upward).N(tlatitudet longitudetaltitude(tselfRRR((scode\common\GeodeticPoint.pyt__init__ s  cCst|j|jƒS(sXReturn a copy of 'self' with the same latitude and longitude, but zero altitude.(tGeoPointRR(R((scode\common\GeodeticPoint.pytzeroAltitudePoint0scCs tj|jƒ}tj|jƒ}tj|jdtj|ƒd|jdtj|ƒdƒ}|jdtj|ƒ|}|jdtj|ƒ|}||j tj|ƒ7}||j tj|ƒ7}|tj|ƒ}|tj|ƒ}|||fS(s!Calculate geocentric cartesian coordinates of the point, in meters, and return them as a 3-tuple '(x, y, z)'. 'z' axis points to the North pole, 'x' axis is aligned with the reference meridian, whereas 'y' axis complements the system to a right-handed one.i( tmathtradiansRRtsqrttellipsoidMajorRadiustcostellipsoidMinorRadiustsinR(Rt latRadianst longRadianst denominatortxytztxty((scode\common\GeodeticPoint.pytcartesianCoords6s"csD|jƒ‰|jƒ‰tjt‡‡fd†tdƒDƒƒƒS(s¬Return the distance between 'self' and 'other' points, measured along a straight line, in meters. Warning: Altitudes of both points are taken into account.c3s%|]}ˆ|ˆ|dVqdS(iN((t.0ti(t selfCoordst otherCoords(scode\common\GeodeticPoint.pys `si(RRR tsumtrange(Rtother((RRscode\common\GeodeticPoint.pyt distToPointWs  cCs®|jƒ\}}}|jƒ\}}}||||} ||||} ||||} |jƒ\} } }| | | | || tj| d| d| dƒS(s×Return the distance between 'self' point and a plane passing through 'pt1', 'pt2', and Earth's center of mass, measured along the normal to the plane, in meters. The return value is positive if 'self' belongs to right-hand hemisphere with respect to the vector 'pt1'->'pt2' and negative otherwise. Exception is raised if points 'pt1' and 'pt2' coincide or are opposite to each other with respect to Earth's center of mass.i(RRR (Rtpt1tpt2tx1ty1tz1tx2ty2tz2txNormtyNormtzNormtx0ty0tz0((scode\common\GeodeticPoint.pyt distFromLinecs cCs4|j|jko3|j|jko3|j|jkS(N(tlatitideRRR(RR((scode\common\GeodeticPoint.pyt__eq__scCs't|jƒt|jƒAt|jƒAS(N(thashRRR(R((scode\common\GeodeticPoint.pyt__hash__„s( t__name__t __module__t__doc__R R RRRRR-R/R1(((scode\common\GeodeticPoint.pyRs   !  ((Rt__all__R(((scode\common\GeodeticPoint.pyts