Computing Vincenty Sphere Distance Matrix in R (Example Code)
This post illustrates how to compute the geospatial Vincenty sphere distance in R.
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 |
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 Sphere Great Circle Distance Matrix of a Set of Geographical Points
install.packages("geosphere") # Install geosphere package library("geosphere") # Load geosphere |
install.packages("geosphere") # Install geosphere package library("geosphere") # Load geosphere
distm(mat_ll, fun = distVincentySphere) # Return geospatial distance # [,1] [,2] # [1,] 0.0 437554.5 # [2,] 437554.5 0.0 |
distm(mat_ll, fun = distVincentySphere) # Return geospatial distance # [,1] [,2] # [1,] 0.0 437554.5 # [2,] 437554.5 0.0