I grabbed a database of the zip codes and their langitudes/latitudes, etc from this This page. It has got the following fields:
ZIP, LATITUDE, LONGIT
It can be done with just maths...
function calc_distance($point1, $point2)
{
$distance = (3958 * 3.1415926 * sqrt(
($point1['lat'] - $point2['lat'])
* ($point1['lat'] - $point2['lat'])
+ cos($point1['lat'] / 57.29578)
* cos($point2['lat'] / 57.29578)
* ($point1['long'] - $point2['long'])
* ($point1['long'] - $point2['long'])
) / 180);
return $distance;
}
This is mike's answer with some annotations for the magic numbers. It seemed to work fine for me for some test data:
function calc_distance($point1, $point2)
{
$radius = 3958; // Earth's radius (miles)
$deg_per_rad = 57.29578; // Number of degrees/radian (for conversion)
$distance = ($radius * pi() * sqrt(
($point1['lat'] - $point2['lat'])
* ($point1['lat'] - $point2['lat'])
+ cos($point1['lat'] / $deg_per_rad) // Convert these to
* cos($point2['lat'] / $deg_per_rad) // radians for cos()
* ($point1['long'] - $point2['long'])
* ($point1['long'] - $point2['long'])
) / 180);
return $distance; // Returned using the units used for $radius.
}
You can also try hitting a web service to calc the distance. Let someone else do the heavy lifting.
https://www.zipcodeapi.com/API#distance
In your zip table you need to grab the coordinates (lat,long) for the two points for which you want to get the distance.
You can then either calculate the distance right in the SQL or with PHP. Both methods are outlined in this post:
http://dev.strategystar.net/2011/10/mysql-php-get-distance-between-two-coordinates-in-miles-kilometers/
The calculation in the example is based on the formula already discussed in this thread. It provides the radius for the earth in both miles and kilometers so it will allow you to get the distance between two points in both units.
The link above is great because the method for calculation does not include any magic numbers, just the radius of the earth!
Check out the Haversine formula for calculating great circle distances between two points. Some more samples can be found here
Haversine formula:
(Note that angles need to be in radians to pass to trig functions).