Computing Vincenty Ellipsoid Distance in R (Example Code)

In this tutorial, I’ll show how to calculate a Vincenty ellipsoid distance in the R programming language.

Example Data

mat_ll <- matrix(c(75.13135, 79.06473,        # Define longitude & latitude points
                   15.28734, 14.29367),
                 ncol = 2)
colnames(mat_ll) <- c("long", "lati")
rownames(mat_ll) <- c("p1", "p2")
mat_ll                                        # Show longitude/latitude matrix in RStudio console
#        long     lati
# p1 75.13135 15.28734
# p2 79.06473 14.29367

Example: Create Vincenty Ellipsoid Great Circle Distance Matrix of a Set of Geographical Points

install.packages("geosphere")                 # Install & load geosphere
library("geosphere")
distm(mat_ll, fun = distVincentyEllipsoid)    # Return geospatial distance
#          [,1]     [,2]
# [1,]      0.0 437475.4
# [2,] 437475.4      0.0

Leave a Reply

Your email address will not be published.

Fill out this field
Fill out this field
Please enter a valid email address.
You need to agree with the terms to proceed

Menu
Top