[R-sig-Geo] How to efficiently generate data of neighboring points

Kent Johnson kent3737 @end|ng |rom gm@||@com
Wed Jun 3 03:02:00 CEST 2020


Roger's example works for me and gives a list of length 101. I did have
some issues that were resolved by updating packages. I'm using R 3.6.3 on
macOS 10.15.4. I also use rtree successfully on Windows 10 with R 3.6.3.

Kent

On Tue, Jun 2, 2020 at 12:29 PM Roger Bivand <Roger.Bivand using nhh.no> wrote:

> On Tue, 2 Jun 2020, Kent Johnson wrote:
>
> > rtree uses Euclidean distance so the points should be in a coordinate
> > system where this makes sense at least as a reasonable approximation.
>
> I tried the original example:
>
> remotes::install_github("hunzikp/rtree")
> library(spData)
> library(sf)
> projdata<-st_transform(nz_height, 32759)
> library(rtree)
> pts <- st_coordinates(projdata)
> rt <- RTree(st_coordinates(projdata))
> bufferR <- c(402.336, 1609.34, 3218.69, 4828.03, 6437.38)
> wd1 <- withinDistance(rt, pts, bufferR[1])
>
> but unfortunately failed (maybe newer Boost headers than yours?):
>
> Error in UseMethod("withinDistance", rTree) :
>    no applicable method for 'withinDistance' applied to an object of class
> "c('list', 'RTree')"
>
> >
> > Kent
> >
> > On Tue, Jun 2, 2020 at 9:59 AM Roger Bivand <Roger.Bivand using nhh.no> wrote:
> >
> >> On Tue, 2 Jun 2020, Kent Johnson wrote:
> >>
> >>>> Date: Tue, 2 Jun 2020 02:44:17 -0500
> >>>> From: Lom Navanyo <lomnavasia using gmail.com>
> >>>> To: r-sig-geo using r-project.org
> >>>> Subject: [R-sig-Geo] How to efficiently generate data of neighboring
> >>>>         points within specified radii (distances) for each point in a
> >> given
> >>>>         points data set.
> >>>>
> >>>
> >>>
> >>>> Hello,
> >>>> I have data set of about 3400 location points with which I am trying
> to
> >>>> generate data of each point and their neighbors within defined radii
> >> (eg,
> >>>> 0.25, 1, and 3 miles).
> >>>>
> >>>
> >>> The rtree package is very fast and memory-efficient for within-distance
> >>> calculations.
> >>> https://github.com/hunzikp/rtree
> >>
> >> Thanks! Does this also apply when the input points are in geographical
> >> coordinates?
> >>
> >> Roger
> >>
> >>>
> >>> Kent Johnson
> >>> Cambridge, MA
> >>>
> >>>       [[alternative HTML version deleted]]
> >>>
> >>> _______________________________________________
> >>> R-sig-Geo mailing list
> >>> R-sig-Geo using r-project.org
> >>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
> >>>
> >>
> >> --
> >> Roger Bivand
> >> Department of Economics, Norwegian School of Economics,
> >> Helleveien 30, N-5045 Bergen, Norway.
> >> voice: +47 55 95 93 55; e-mail: Roger.Bivand using nhh.no
> >> https://orcid.org/0000-0003-2392-6140
> >> https://scholar.google.no/citations?user=AWeghB0AAAAJ&hl=en
> >>
> >
>
> --
> Roger Bivand
> Department of Economics, Norwegian School of Economics,
> Helleveien 30, N-5045 Bergen, Norway.
> voice: +47 55 95 93 55; e-mail: Roger.Bivand using nhh.no
> https://orcid.org/0000-0003-2392-6140
> https://scholar.google.no/citations?user=AWeghB0AAAAJ&hl=en
>

	[[alternative HTML version deleted]]



More information about the R-sig-Geo mailing list