Lionel Hertzog shows off spatial capabilities in R:

All of these operations follow the same logic, st_

operation(A, B) checks for each combinations of the geometries in A and B whether AoperationB is true or false. For instance`st_within(A, B)`

checks whether the geometries in A arewithinB, this is similar to`st_contains(B, A)`

, the difference between the two being the shape of the returned object. If A has n geometries and B has m,`st_contains(B, A)`

returns a list of length m where each elements contains the row IDs (numbers between 1 and n) of the geometries in A satisfying the operation. By using`sparse=FALSE`

the functions returns matrices, like`st_within(A, B, sparse=FALSE)`

returns a n x m matrix,`st_within(B, A, sparse=FALSE)`

returns a m x n matrix. Note that running st_operation(A, A) checks the operation between all geometries of the object, so returning a n x n matrix.

Click through for part 1 of the series.