Misaligned sf object in ggplot

Question

I have two sf objects, which when plotting are slightly misaligned.

While initially the two objects have different CRS:

I proceed, as usual with:

sf2 <- st_transform(sf2, crs = st_crs(sf1))

to align the CRS. However still, when plotting this:

ggplot() 
  + geom_sf(sf1)
  + geom_sf(sf2, color = "green", fill = "transparent")

I get the attached output. Any ideas what I am doing wrong?

Datasources are:

Likely red herring, but set ` sf_use_s2(FALSE)` prior to st_transform. — Chris, Commented Jul 2 at 20:44
Where/how did you source your data? Some example data will be helpful here. Subset sf1 and sf2 to only the single most northwestern polygon, run dput() on both subsets, and add the results to a codeblock in your question. The results will look something like structure(list(...)). — L Tyrone, Commented Jul 2 at 23:21
Links to the source/s of your data will also be useful in order to rule out any processing errors you may have accidentally made e.g. the data may have already been corrupted. — L Tyrone, Commented Jul 2 at 23:30
many thanks for the comments: I found a "workaround" when exporting a subsample (see Answer). I however still do not understand the problem. I have updated the question by adding the data sources. I couldn't however attach the dput() output, as it was too long. — Tom, Commented Jul 3 at 8:34

Tom · Accepted Answer · 2024-07-03 08:33:12Z

0

The "solution" to the problem was, that:

all of a sudden the polygons align, even though running identical(st_crs(sf1), st_crs(sf2))gives me FALSE

I still do not understand the issue - but it worked.

answered Jul 3 at 8:33

Tom

235 bronze badges

Your "solution" still gives the wrong map, as in your question, for me.
– Edzer Pebesma
Commented Jul 16 at 14:50

Add a comment |

1 Answer 1