def distance2To(self, other, radius=R_M, adjust=False, wrap=False): '''Compute the distance between this and an other geohash using the U{Equirectangular Approximation / Projection <https://www.Movable-Type.co.UK/scripts/latlong.html>}. @arg other: The other geohash (L{Geohash}). @kwarg radius: Mean earth radius (C{meter}) or C{None}. @kwarg adjust: Adjust the wrapped, unrolled longitudinal delta by the cosine of the mean latitude C{bool}). @kwarg wrap: Wrap and unroll longitudes (C{bool}). @return: Approximate distance (C{meter}, same units as I{radius}) or distance squared (C{degrees} squared) if I{radius} is C{None} or 0. @raise TypeError: The I{other} is not a L{Geohash}, C{LatLon} or C{str}. @see: U{Local, flat earth approximation <https://www.EdWilliams.org/avform.htm#flat>}, functions ''' other = _2Geohash(other) a1, b1 = self.latlon a2, b2 = other.latlon if radius: return equirectangular(a1, b1, a2, b2, radius=radius, adjust=adjust, limit=None, wrap=wrap) else: return equirectangular_(a1, b1, a2, b2, adjust=adjust, limit=None, wrap=wrap)[0]
def equirectangular3(lat1, lon1, lat2, lon2, **options): '''DEPRECATED, use function C{equirectangular_}. @return: 3-Tuple C{(distance2, delta_lat, delta_lon)}. ''' from pygeodesy.formy import equirectangular_ return tuple(equirectangular_(lat1, lon1, lat2, lon2, **options)[:3])
def d2yxu(self, i, j): '''Return distance squared, points[i] to -[j] deltas and longitudinal unrollment. ''' p1 = self.pts[i] p2 = self.pts[j] return equirectangular_(p1.lat, p1.lon, p2.lat, p2.lon, **self.options)