View Full Version : Finding the nearest road given a lat long.

Anonymous

09-11-2002, 10:37 AM

Hi,

I want to find the nearest road from a lat/long the reason I want this is because of the inherant inaccuracy in GPS's sometimes my car appears off road (only slightly) but when I try to calculate a route based on these waypoints it goes mad and totaly gets it wrong.

Thanks

John Meyer

09-11-2002, 08:16 PM

Check out this article: Reverse Geocoding, Pt. III by Walt Cygan

http://www.mp2kmag.com/articles.asp?ArticleID=50

dzhong

10-10-2002, 01:06 PM

Hi Scarr,

I encontered the same problem. Could you give me an update once you have a fix?

Thanks in advance!

larry_llama

12-22-2002, 10:09 PM

I have seen this issue pop up a lot, and I have an idea that I am going to try (though i am not quite there yet).. if you are tracking GPS then you know your compass heading. if you are moving along a road, and your gps is returning an invalid point, it must be either to the left or right of you. so if you take a line perpendicular to your compass heading, you can start looking along that line in both directions until you can get a valid location. if you check like 100 feet to either side of you, you should hit the road you want (and not the wrong one). just a thought, but it cuts down on the number of iterations required in the method where you spiral outward 360 degrees around your gps point.

hope this makes sense. post here if you get this working and i will do the same!

ps is there any way to get the car icon up on the map screen (like in streets and trips application during gps tracking)

cwest

04-04-2005, 03:33 PM

I’m new to MapPoint, and I want to drag a line from a point that shows distance and direction. I thought pointing a line in a give direction—say 286 degrees—would be an obvious feature similar to “Measure Distance.” Have I missed something or isn’t this doable?

Wilfried

04-05-2005, 04:22 AM

Hi,

I just encountered a similar problem (with other purpose, but my solution may be working for you). What I do is running in a circle around the point where the car 'pretends' to be. I start with a radius of 5 meter and check the circumfence of the circle every 5 meter. If not found I increment the radius by 5 meter till a street is found. In your case I only should do a small maximum of 10..15 meter because you can end up with a wrong street if car is on a non digitized street.

This is the code (in Delphi). DistanceRes is the resollution (I set it to 5 meter), and MaxDistance is the radius where the routine should give up. FindAddr is the routine that ask Mappoint for streetresults (not listed here).

function TMapPoint.CalcInCircle(CenterLat, CenterLon: double; Streets: TStrings; out Distance, Angle: integer): boolean;

var

Lat, Lon: double;

Radius: double; // radius of the circle in degree

Circum: integer; // circumfence of the circle in Distance

Count: integer; // the amount of points to calculate on the circumfence

AngleInc: integer; // angle increments in degree per calculation

begin

Distance := DistanceRes;

while Distance <= MaxDistance do begin

Radius := 1 / 60 / 1852 * Distance;

Circum := Trunc(Distance * 2 * PI);

Count := Circum div DistanceRes;

AngleInc := 360 div Count;

if AngleInc <= 0 then begin

Result := False;

Exit;

end;

Angle := 0;

while Angle < 360 do begin

CalcXY(CenterLat, CenterLon, Angle, Radius, Lat, Lon);

Result := FindAddr(Lat, Lon, Streets);

if Result then

Exit;

Inc(Angle, AngleInc);

end;

Inc(Distance, DistanceRes);

end;

Result := False;

end;

And this is the actulally working horse:

procedure TMapPoint.CalcXY(centerLat, centerLon, Angle, Radius: double; out Lat, lon: double);

var

a, b, c, r: double;

begin

{ B

/ |

/ |

/ |

/ |

c / | a

/ |

/ |

/ |

/_________________|

A b C

given A en c

a = c * sin(A)

b = c * cos(A)

A == 0 is direction East in this drawing

a = Latitude offset

b = Longitude offset

c = Radius of circle }

r := DegToRad(Angle);

c := Radius;

a := c * Sin(r);

b := c * Cos(r) / Cos(DegToRad(centerLat)); // correct longitude for circle

Lat := CenterLat + a;

Lon := CenterLon + b;

end;

micheln

12-04-2008, 05:54 PM

Hi Wilfried

I have a similar need. Actually, I need to show the nearest POI to the actual position of the vehicle. This info should be updated every 2 minutes or so. My POI table has some 1000s records.

What I'm doing so far is, in a separate thread, to load the POI data in an in-memory dataset (delphi kbmMemTable) and calculate the distance between every POI and the given position (for loop). I store the results (actual position - POI infos - distance) in an array, and finally select the closest POI (min distance).

I know what I'm doing is REALLY awful :oops2: but since I didn't have much time when I did it, I didn't really have the choice (optimizing and taking time to search for the better solution when my boss is crying about not meeting an impossible deadline isn't the easiest thing to do :p

Anyway, I would like to optimize this now, so Could you give me a hand in adapting your code to my need?

Best Regards,

Michel

Powered by vBulletin® Version 4.2.2 Copyright © 2017 vBulletin Solutions, Inc. All rights reserved.

Search Engine Friendly URLs by

vBSEO 3.6.1