title: Projection issues point over multipolygon data sf objects
body: <p>I'm trying to plot point data (coördinates) on a multipolygoon map. But I can't get them on top of each other and in the right shape.</p>
  
  <p>Underlying projection data of the Netherlands <code>multipolygon</code> data:</p>
  
  <pre><code>Simple feature collection with 1 feature and 3 fields
  geometry type:  MULTIPOLYGON
  dimension:      XY
  bbox:           xmin: 13565.4 ymin: 306846.2 xmax: 278026.1 ymax: 619232.6
  epsg (SRID):    28992
  proj4string:    +proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel +towgs84=565.4171,50.3319,465.5524,-0.398957,0.343988,-1.87740,4.0725 +units=m +no_defs
  </code></pre>
  
  <p><a href="https://www.imergis.nl/gpkg/2019_voorlopige_bestuurlijke_grenzen_imergis_gpkg.zip" rel="nofollow noreferrer">Downloadlink to a <code>.zip</code> with 4 Geopackage</a> files from the Netherlands</p>
  
  <p>Underlying <code>point</code> projection data of coordinates:</p>
  
  <pre><code>Simple feature collection with 1466 features and 4 fields
      geometry type:  POINT
      dimension:      XY
      bbox:           xmin: 3.389737 ymin: 50.77106 xmax: 7.187327 ymax: 53.44402
      epsg (SRID):    28992
      proj4string:    +proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel +towgs84=565.4171,50.3319,465.5524,-0.398957,0.343988,-1.87740,4.0725 +units=m +no_defs
  </code></pre>
  
  <p>some rows of <code>df</code> dataframe:</p>
  
  <pre><code>df &lt;- structure(list(brouwerijnaam = c("Almelosche", "HertogJan", "HertogJan", 
  "Bavaria", "Bavaria", "Brand"), begin_datum_jaar = c(1987, 1981, 
  1981, 1719, 1719, 1340), eind_datum_jaar = c(1990, 0, 0, 0, 0, 
  0), plaats = c("Almelo", "Arcen", "Arcen", "Lieshout", "Lieshout", 
  "Wijlre"), lon = c(52.35192, 51.473427, 51.473427, 51.518607, 
  51.518607, 50.832695), lat = c(6.660562, 6.181219, 6.181219, 
  5.597416, 5.597416, 5.895593)), row.names = c(NA, -6L), class = c("tbl_df", 
  "tbl", "data.frame"))
  </code></pre>
  
  <p>some rows of <code>points</code> dataframe:</p>
  
  <pre><code>points &lt;- structure(list(brouwerijnaam = c("Almelosche", "HertogJan", "HertogJan", 
  "Bavaria", "Bavaria", "Brand"), begin_datum_jaar = c(1987, 1981, 
  1981, 1719, 1719, 1340), eind_datum_jaar = c(1990, 0, 0, 0, 0, 
  0), plaats = c("Almelo", "Arcen", "Arcen", "Lieshout", "Lieshout", 
  "Wijlre"), geometry = structure(list(structure(c(52.35192, 6.660562
  ), class = c("XY", "POINT", "sfg")), structure(c(51.473427, 6.181219
  ), class = c("XY", "POINT", "sfg")), structure(c(51.473427, 6.181219
  ), class = c("XY", "POINT", "sfg")), structure(c(51.518607, 5.597416
  ), class = c("XY", "POINT", "sfg")), structure(c(51.518607, 5.597416
  ), class = c("XY", "POINT", "sfg")), structure(c(50.832695, 5.895593
  ), class = c("XY", "POINT", "sfg"))), class = c("sfc_POINT", 
  "sfc"), precision = 0, bbox = structure(c(xmin = 50.832695, ymin = 5.597416, 
  xmax = 52.35192, ymax = 6.660562), class = "bbox"), crs = structure(list(
      epsg = NA_integer_, proj4string = NA_character_), class = "crs"), n_empty = 0L)), row.names = c(NA, 
  -6L), sf_column = "geometry", agr = structure(c(brouwerijnaam = NA_integer_, 
  begin_datum_jaar = NA_integer_, eind_datum_jaar = NA_integer_, 
  plaats = NA_integer_), .Label = c("constant", "aggregate", "identity"
  ), class = "factor"), class = c("sf", "tbl_df", "tbl", "data.frame"
  ))
  </code></pre>
  
  <p>Code I used to shape the x and y coördinates in a sf object with the same crs of multipolygon data.</p>
  
  <pre><code>library(sf)
  library(tidyverse)
  
  points &lt;- df %&gt;%
    filter(!is.na(lon)) %&gt;%
    st_as_sf(coords = c("lat", "lon")) 
  
  points &lt;- st_set_crs(brouwerijen_sf, 28992)
  
  ggplot() +
    geom_sf(data = multipolygon) +
    geom_sf(data = points) +
    coord_sf(datum = NA)
  </code></pre>
  
  <p><a href="https://i.stack.imgur.com/SrYrs.png" rel="nofollow noreferrer"><img src="https://i.stack.imgur.com/SrYrs.png" alt="multipolygon"></a>
  <a href="https://i.stack.imgur.com/OnbQl.png" rel="nofollow noreferrer"><img src="https://i.stack.imgur.com/OnbQl.png" alt="Point data"></a>
  <a href="https://i.stack.imgur.com/peQHk.png" rel="nofollow noreferrer"><img src="https://i.stack.imgur.com/peQHk.png" alt="enter image description here"></a></p>