[R-sig-Geo] gDistance on a LARGE Spatial Object?

Peter Versteegen pitiur at gmail.com
Sun Oct 28 13:22:14 CET 2012


Wouter Steenbeek <WSteenbeek <at> nscr.nl> writes:
 
Your problem suggests that the spatial extent is not very large.
I have a similar problem with over 600,000 data points.
They describe bathymetric measurements of the lake.
I'm working on making contours..

Project all of your coordinates on a plane.
Means that you can work with meters instead of lat/lon.
Euclidian distances are easily calculated, d = sqrt((x-x1)^2 + (y-y1)^2).
You have to write a few lines of basic code in R.

Setup the problem with a database or write/read a textfile.

Routine computations like this 600,000 times
looking at each of the 600,000 points,
takes under a minute on my aging Mac.


I hope this helps.



More information about the R-sig-Geo mailing list