Loading the geosphere package you can use a number of different functions library(geosphere) distm(c(lon1, lat1), c(lon2, lat2), fun = distHaversine) Also: distHaversine() distMeeus() distRhumb() distVincentyEllipsoid() distVincentySphere() …
Thanks! Here is a C# version of OP’s solutions using degrees. It utilises the System.Device.Location.GeoCoordinate class public static GeoCoordinate GetCentralGeoCoordinate( IList<GeoCoordinate> geoCoordinates) { if (geoCoordinates.Count == 1) { return geoCoordinates.Single(); } double x = 0; double y = 0; double z = 0; foreach (var geoCoordinate in geoCoordinates) { var latitude = geoCoordinate.Latitude * Math.PI … Read more
Here are two options, one using sf and one using sp package functions. sf is the more modern (and, here in 2020, recommended) package for analyzing spatial data, but in case it’s still useful, I am leaving my original 2012 answer showing how to do this with sp-related functions. Method 1 (using sf): library(sf) library(spData) … Read more
Needed to convert answers from radians back to degrees. Working code below: import math R = 6378.1 #Radius of the Earth brng = 1.57 #Bearing is 90 degrees converted to radians. d = 15 #Distance in km #lat2 52.20444 – the lat result I’m hoping for #lon2 0.36056 – the long result I’m hoping for. … Read more
Based on another question on stackoverflow, I got this code.. This calculates the result in meters, not in miles 🙂 public static float distFrom(float lat1, float lng1, float lat2, float lng2) { double earthRadius = 6371000; //meters double dLat = Math.toRadians(lat2-lat1); double dLng = Math.toRadians(lng2-lng1); double a = Math.sin(dLat/2) * Math.sin(dLat/2) + Math.cos(Math.toRadians(lat1)) * Math.cos(Math.toRadians(lat2)) … Read more
Best option is to use libraries sp and rgeos, which enable you to construct spatial classes and perform geoprocessing. library(sp) library(rgeos) Read the data and transform them to spatial objects: mydata <- read.delim(‘d:/temp/testfile.txt’, header=T) sp.mydata <- mydata coordinates(sp.mydata) <- ~long+lat class(sp.mydata) [1] “SpatialPointsDataFrame” attr(,”package”) [1] “sp” Now calculate pairwise distances between points d <- gDistance(sp.mydata, … Read more
See this example, initializes the map to “San Diego, CA”. Uses the Google Maps Javascript API v3 Geocoder to translate the address into coordinates that can be displayed on the map. <html> <head> <meta name=”viewport” content=”initial-scale=1.0, user-scalable=no”/> <meta http-equiv=”content-type” content=”text/html; charset=UTF-8″/> <title>Google Maps JavaScript API v3 Example: Geocoding Simple</title> <script type=”text/javascript” src=”http://maps.google.com/maps/api/js?sensor=false”></script> <script type=”text/javascript”> var … Read more
The GeoCoordinate class (.NET Framework 4 and higher) already has GetDistanceTo method. var sCoord = new GeoCoordinate(sLatitude, sLongitude); var eCoord = new GeoCoordinate(eLatitude, eLongitude); return sCoord.GetDistanceTo(eCoord); The distance is in meters. You need to reference System.Device.
I think your question says you have the city values for the two cities between which you wish to compute the distance. This query will do the job for you, yielding the distance in km. It uses the spherical cosine law formula. Notice that you join the table to itself so you can retrieve two … Read more
The Java code given by Dommer above gives slightly incorrect results but the small errors add up if you are processing say a GPS track. Here is an implementation of the Haversine method in Java which also takes into account height differences between two points. /** * Calculate distance between two points in latitude and … Read more