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

[Lom Navanyo]

I had the errors with rtree using R 3.6.3. I have since changed to R 4.0.0
but I got the same error.

And  yes, for Roger's example, I have the objects wd1, ... wd4, all with
length 101. I think my difficulty is my inability to output the list
detailing the point IDs t50_fid.

---------
Lom

On Tue, Jun 2, 2020 at 8:02 PM Kent Johnson <kent3737 using gmail.com> wrote:

> 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]]