PostGIS  2.1.10dev-r@@SVN_REVISION@@
double edge_distance_to_edge ( const GEOGRAPHIC_EDGE e1,
const GEOGRAPHIC_EDGE e2,
GEOGRAPHIC_POINT closest1,
GEOGRAPHIC_POINT closest2 
)

Calculate the distance between two edges.

IMPORTANT: this test does not check for edge intersection!!! (distance == 0) You have to check for intersection before calling this function.

Definition at line 1216 of file lwgeodetic.c.

References edge_distance_to_point(), GEOGRAPHIC_EDGE::end, and GEOGRAPHIC_EDGE::start.

Referenced by circ_tree_distance_tree_internal(), ptarray_distance_spheroid(), and test_edge_distance_to_edge().

1217 {
1218  double d;
1219  GEOGRAPHIC_POINT gcp1s, gcp1e, gcp2s, gcp2e, c1, c2;
1220  double d1s = edge_distance_to_point(e1, &(e2->start), &gcp1s);
1221  double d1e = edge_distance_to_point(e1, &(e2->end), &gcp1e);
1222  double d2s = edge_distance_to_point(e2, &(e1->start), &gcp2s);
1223  double d2e = edge_distance_to_point(e2, &(e1->end), &gcp2e);
1224 
1225  d = d1s;
1226  c1 = gcp1s;
1227  c2 = e2->start;
1228 
1229  if ( d1e < d )
1230  {
1231  d = d1e;
1232  c1 = gcp1e;
1233  c2 = e2->end;
1234  }
1235 
1236  if ( d2s < d )
1237  {
1238  d = d2s;
1239  c1 = e1->start;
1240  c2 = gcp2s;
1241  }
1242 
1243  if ( d2e < d )
1244  {
1245  d = d2e;
1246  c1 = e1->end;
1247  c2 = gcp2e;
1248  }
1249 
1250  if ( closest1 ) *closest1 = c1;
1251  if ( closest2 ) *closest2 = c2;
1252 
1253  return d;
1254 }
Point in spherical coordinates on the world.
Definition: lwgeodetic.h:33
GEOGRAPHIC_POINT start
Definition: lwgeodetic.h:44
double edge_distance_to_point(const GEOGRAPHIC_EDGE *e, const GEOGRAPHIC_POINT *gp, GEOGRAPHIC_POINT *closest)
Definition: lwgeodetic.c:1165
GEOGRAPHIC_POINT end
Definition: lwgeodetic.h:45

Here is the call graph for this function:

Here is the caller graph for this function: