File Exchange

image thumbnail

LatLon distance

version 1.0 (1.64 KB) by

Find distance between two points based on latlon coordinates



View License

This functions finds distance (in km) between two points on Earth using latitude-longitude coordinates of the two points. Distance is calculated using two distance funstions: Haversine and Pythagoran. The first output is based on Haversine function, which is more accurate especially for longer distances.

latlon1: latlon of origin point [lat lon]
latlon2: latlon of destination point [lat lon]

d1km: distance calculated by Haversine formula
d2km: distance calculated based on Pythagoran theorem

--Example 1, short distance:
latlon1=[-43 172];
latlon2=[-44 171];
[d1km d2km]=distance(latlon1,latlon2)
d1km =
137.365669065197 (km)
d2km =
137.368179013869 (km)
%d1km approximately equal to d2km

--Example 2, longer distance:
latlon1=[-43 172];
latlon2=[20 -108];
[d1km d2km]=distance(latlon1,latlon2)
d1km =
10734.8931427602 (km)
d2km =
31303.4535270825 (km)
d1km is significantly different from d2km (d2km is not able to work for longer distances).

Comments and Ratings (15)


Works perfectly.

Manoj Pai

How to apply inputs

Maninder Ade

Aysu Özel

M Sohrabinia

M Sohrabinia (view profile)

@Seb: decimal degrees as shown in the examples. Thanks for rating the code.

M Sohrabinia

M Sohrabinia (view profile)

@Justin: thanks for the suggestions Justin, I will try to get back to this, have been away from Matlab for the last couple of years. I acknowledge I have not spent a lot work on handling the exceptions, the aim was to calculate the correct distance which I believe it does.

Rafael D

Seb Biass


Just a question, are inputs in deg min sec or decimal degrees?



Straight forward. Handy function for GPS users

This code does what it says it should but there are two 2-minute changes to make it much better:

1) make all operations vectorized so you can put a list of lat,lon coordinates. Basically, add a "." before everything. Super easy and much better

2) Add a check for NaN's. If two points are too close, you get NaN when it should be zero. Very easy to do:

>> d1km(isnan(d1km)) = 0;

I hope this helps

Tianyang Li

Ian Shapiro

Simple but effective. Does just what it needs to do.

MATLAB Release
MATLAB 7.7 (R2008b)

Inspired: pathdist

Download apps, toolboxes, and other File Exchange content using Add-On Explorer in MATLAB.

» Watch video

Win prizes and improve your MATLAB skills

Play today