R version 4.2.2 Patched (2022-11-10 r83330) -- "Innocent and Trusting"
Copyright (C) 2022 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> options("rgdal_show_exportToProj4_warnings"="none")
> library(sp)
> data(meuse)
> data(meuse.grid)
> coordinates(meuse) = ~x+y
> coordinates(meuse.grid) = ~x+y
> x = coordinates(meuse)
> y = coordinates(meuse.grid)
> out = spDists(meuse,meuse.grid)
> out2 = as.matrix(dist(rbind(coordinates(meuse),coordinates(meuse.grid))))
> out2 = out2[1:155,155+1:nrow(y)]
> # should be equal:
> sum(out2 - out)
[1] 3.666401e-12
> summary(as.vector(out2 - out))
Min. 1st Qu. Median Mean 3rd Qu. Max.
-9.095e-13 0.000e+00 0.000e+00 7.600e-18 0.000e+00 9.095e-13
>
> out = spDists(meuse.grid,meuse)
> sum(out2 - t(out))
[1] 3.666401e-12
> summary(as.vector(out2 - t(out)))
Min. 1st Qu. Median Mean 3rd Qu. Max.
-9.095e-13 0.000e+00 0.000e+00 7.600e-18 0.000e+00 9.095e-13
>
> proc.time()
user system elapsed
1.044 0.862 1.024
Generated by dwww version 1.15 on Sun Jun 30 10:38:48 CEST 2024.