8

I am working with Latitude / Longitude coordinates in a google map.

I have two lines :

  • Line A : 48.31508162629726, -2.591741396838972 to 48.40216156645915, -2.2218462112093404
  • Line B : 48.383816077371215, -2.274292940053768 to 48.66103546935337, -1.7066197241571377

I then use the following formula to find the point where they cross.

var XAsum = A.LngStart - A.LngEnd;
var XBsum = B.LngStart - B.LngEnd;
var YAsum = A.LatStart - A.LatEnd;
var YBsum = B.LatStart - B.LatEnd;

var LineDenominator = XAsum * YBsum - YAsum * XBsum;
if(LineDenominator == 0.0)
    return false;

var a = A.LngStart * A.LatEnd - A.LatStart * A.LngEnd;
var b = B.LngStart * B.LatEnd - B.LatStart * B.LngEnd;

var x = (a * XBsum - b * XAsum) / LineDenominator;
var y = (a * YBsum - b * YAsum) / LineDenominator;

This tells me that the lines do indeed cross and returns the x and y values.

However, when I plot the returned point, it is offset (not much) from the real intersection.

Is there a better and just as fast algorithm I could use which will return me the correct intersection point ?

It needs to be fast as I am iterating over a large number of lines (~1000).

EDIT : Note this is giving me an error offset of around 7.5 meters

Simon
  • 2,056
  • 4
  • 31
  • 41
  • What language are you programming in? Are there any libraries you can use (or are using) to help you solve GIS/spatial-analysis problems? – Jeff B Jul 25 '12 at 16:39
  • Did you solve your problem? I have the same challenge to find the intersection of two lines (arcs) on google maps. The Cartesian intersection produces errors on the map. For each arc I have a pair of latlons and I need to find their intersection on google maps. – Ash Apr 20 '15 at 01:34
  • No, I had to live with the error. Sorry – Simon Apr 20 '15 at 05:55

1 Answers1

4

I'm assuming the algorithm you're using is the one for finding line intersections on a Cartesian coordinate system (i.e. a plane). Unfortunately, the earth is not a plane (or even a sphere) so using that algorithm will introduce error. Google Maps uses an ellipsoid (specifically WGS84) to approximate the surface of the earth, so you'll need an algorithm for finding the intersections of arcs on an ellipsoid.

This page might contains some helpful information: http://mathhelpforum.com/calculus/90196-point-intersection-two-lines.html

Jeff B
  • 8,123
  • 17
  • 63
  • 134
  • thank you for taking the trouble to answer. You are right in your assumptions. I looked at the page you mention, but I am not sufficiently advanced (!) in maths to be able to translate what is indicated to my use. Would you have any further suggestions ? – Simon Jul 25 '12 at 18:28
  • Any particular language you care to use? – Jeff B Jul 26 '12 at 00:57
  • Didn't have time to do much research, but this might help? http://geotools.org/ Or maybe this? http://geospatialmethods.org/spheres/ – Jeff B Jul 26 '12 at 17:40