r - 在 LatLong 中工作时，gridDistance 是否计算大地距离？（R中的光栅包）

Question

R（光栅包）中的 gridDistance 函数生成具有给定值的单元格的距离图。该函数以当前投影的单位给出距离，如果在 LatLong 中工作，则以米为单位。我想知道使用 LatLong 距离是否被“校正”（大地测量）。否则我可能需要重新投影到等距投影。谢谢

Answer 1

调查真正优秀的raster软件包的来源，看起来你的直觉是正确的......

# From gridDistance() when lonlat is true
if (lonlat) {
    distance <- pointDistance(

# From pointDistance()
if (! longlat ) {
    return( .planedist(p1[,1], p1[,2], p2[,1], p2[,2]) )
} else { 
    return( .haversine(p1[,1], p1[,2], p2[,1], p2[,2], r=6378137) )
}

# Finally from .haversine()
.haversine <- function(x1, y1, x2, y2, r=6378137) {
adj <- pi / 180
x1 <- x1 * adj
y1 <- y1 * adj
x2 <- x2 * adj
y2 <- y2 * adj
x <- sqrt((cos(y2) * sin(x1-x2))^2 + (cos(y1) * sin(y2) - sin(y1) * cos(y2) * cos(x1-x2))^2)
y <- sin(y1) * sin(y2) + cos(y1) * cos(y2) * cos(x1-x2)
return ( r * atan2(x, y) )
}

所以简而言之，是的，你得到了超精确的大地距离。