Nothing
# Author: Robert J. Hijmans
# Date : June 2008
# Version 0.8 (taken from Raster package)
# Licence GPL v3
# Vincenty formula for a sphere
# http://en.wikipedia.org/wiki/Great_circle_distance
distVincentySphere <- function(p1, p2, r=6378137) {
toRad <- pi / 180
p1 <- .pointsToMatrix(p1) * toRad
if (missing(p2)) {
p2 <- p1[-1, ,drop=FALSE]
p1 <- p1[-nrow(p1), ,drop=FALSE]
} else {
p2 <- .pointsToMatrix(p2) * toRad
}
p <- cbind(p1[,1], p1[,2], p2[,1], p2[,2], as.vector(r))
lon1 <- p[,1]
lat1 <- p[,2]
lon2 <- p[,3]
lat2 <- p[,4]
x <- sqrt((cos(lat2) * sin(lon1-lon2))^2 + (cos(lat1) * sin(lat2) - sin(lat1) * cos(lat2) * cos(lon1-lon2))^2)
y <- sin(lat1) * sin(lat2) + cos(lat1) * cos(lat2) * cos(lon1-lon2)
dist <- p[,5] * atan2(x, y)
return ( as.vector(dist ))
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.