swh:1:snp:a4c99a50dc49f82b591f268001b320f8c3ca0041
Tip revision: ce60f670aac0a708d3a9af5f0cf46a752d46ba8a authored by John M Chambers on 28 October 2020, 07:59:48 UTC
version 1.0-6.1
version 1.0-6.1
Tip revision: ce60f67
geoXY.R
"geoXY" <-
function(latitude, longitude,
lat0 = min(latitude, na.rm=TRUE),
lon0 = min(longitude, na.rm=TRUE),
unit = 1.) {
lat0 <- rep(lat0, length(latitude))
lon0 <- rep(lon0, length(longitude))
yDist <- geoDist(lat0, lon0, latitude, lon0)
xDist <- geoDist(lat0, lon0, lat0, longitude)
cbind(X= xDist * sign(longitude - lon0)/unit,
Y = yDist * sign(latitude - lat0)/unit)
}