PostGIS  2.3.8dev-r@@SVN_REVISION@@

◆ edge_distance_to_edge()

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 1234 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().

1235 {
1236  double d;
1237  GEOGRAPHIC_POINT gcp1s, gcp1e, gcp2s, gcp2e, c1, c2;
1238  double d1s = edge_distance_to_point(e1, &(e2->start), &gcp1s);
1239  double d1e = edge_distance_to_point(e1, &(e2->end), &gcp1e);
1240  double d2s = edge_distance_to_point(e2, &(e1->start), &gcp2s);
1241  double d2e = edge_distance_to_point(e2, &(e1->end), &gcp2e);
1242 
1243  d = d1s;
1244  c1 = gcp1s;
1245  c2 = e2->start;
1246 
1247  if ( d1e < d )
1248  {
1249  d = d1e;
1250  c1 = gcp1e;
1251  c2 = e2->end;
1252  }
1253 
1254  if ( d2s < d )
1255  {
1256  d = d2s;
1257  c1 = e1->start;
1258  c2 = gcp2s;
1259  }
1260 
1261  if ( d2e < d )
1262  {
1263  d = d2e;
1264  c1 = e1->end;
1265  c2 = gcp2e;
1266  }
1267 
1268  if ( closest1 ) *closest1 = c1;
1269  if ( closest2 ) *closest2 = c2;
1270 
1271  return d;
1272 }
Point in spherical coordinates on the world.
Definition: lwgeodetic.h:47
GEOGRAPHIC_POINT start
Definition: lwgeodetic.h:58
double edge_distance_to_point(const GEOGRAPHIC_EDGE *e, const GEOGRAPHIC_POINT *gp, GEOGRAPHIC_POINT *closest)
Definition: lwgeodetic.c:1183
GEOGRAPHIC_POINT end
Definition: lwgeodetic.h:59
Here is the call graph for this function:
Here is the caller graph for this function: