Snap to grid with sf in R

Question

Is there a snap to grid function available in sf for R similar to ST_SnapToGrid(geometry geomA, float size)in PostGIS? It is mentioned as a solution to non-noded intersection problem for PostGIS and am looking for a solution in R. Using other snap options with just the polygon lead to geometry errors (even with small tolerance).

Perhaps all coordinates can be rounded to certain precision, for the same effect as ST_SnapToGrid — Michael Dorman, Feb 03 '19 at 13:17
I'm not sure I understand your suggestion. I ran into errors when I tried to simplify my features as bits of some polygons were missing bits of other polygons or st_snap. — user3386170, Feb 03 '19 at 19:07

score 2 · Answer 1 · answered Feb 04 '19 at 06:19

2

Rounding coordinates may be equivalent to snapping all of the shapes coordinates to a regular grid.

Before -

> pnt = st_point(c(0,0))
> pol = st_buffer(pnt, 1)
> plot(pol)

After -

> pol[[1]] = round(pol[[1]], 1)
> plot(pol)

answered Feb 04 '19 at 06:19

Michael Dorman

783
5
13

I'll try this later and mark it as accepted once I've checked that it works, but it looks promising! – user3386170 Feb 04 '19 at 14:31

score 2 · Answer 2 · answered Mar 29 '19 at 22:35

2

The R equivalent of the ST_SnapToGrid in PostGIS is in the lwgeom package:

# Snap to grid of 5000 m
lwgeom::st_snap_to_grid(x, 5000)

Works well to solve the non-noded intersection problem, and is quicker than applying a buffer of the same tolerance.

answered Mar 29 '19 at 22:35

avianattackarmada

123
4

Snap to grid with sf in R

2 Answers2