尝试

df <- read.table(sep=",", col.names=c("lat1", "lon1", "lat2", "lon2"), text="
52,4,52,13 
39,116,52,13")
library(geosphere)
distHaversine(df[, 2:1], df[, 4:3]) / 1000 # Haversine distance in km