eb08d1ea38ef34b757b7a5b0736ce92e2a4b0fb0,pykrige/core.py,,great_circle_distance,#Any#Any#Any#Any#,58

Before Change

    lat1 = np.array(lat1)*np.pi/180.0
    lat2 = np.array(lat2)*np.pi/180.0
    d = np.cos((lon1-lon2)*np.pi/180.0)*np.cos(lat1)*np.cos(lat2) \
        + np.sin(lat1)*np.sin(lat2)
    // Angle is arccos of euclidean dot product. Avoid errors caused
    // by numerics (possibly d>1.0 or d<-1.0, resulting in NAN. This
    // can, however, only be caused by numerical errors for real
    // lat/lon):
    d[d>1.0] = 1.0
    d[d<-1.0]=-1.0
    return 180.0/np.pi*np.arccos(d)

def euclid3_to_great_circle(euclid3_distance):
    
    Convert euclidean distance between points on a unit sphere to

After Change

    // Formula can be obtained from [2] combining eqns. (14)-(16)
    // for spherical geometry (f=0).
    
    return 180.0/np.pi*np.arctan2(
               np.sqrt((c2*np.sin(dlon))**2 + 
                       (c1*s2-s1*c2*cd)**2),
               s1*s2+c1*c2*cd)

def euclid3_to_great_circle(euclid3_distance):
    
    Convert euclidean distance between points on a unit sphere to

Instances