Probably only works mapping-wise for years with limited number and control cells.
Hint: you can create also centroids for the polygons which are point data and much more efficient to plot. This might also allow you to plot several thousand observations in one map.use the sf package and st_centroid function to convert polygons to points, if necessary.
Definition of Done:
complete data-base with treatment and control cells for all years is exported as geopackage (use write_sf) for later mapping in QGIS
example maps for selected years are created in the Rmd script that contains the matching descriptives with mapView. You can add several years to the map by concetenating several mapviews -> mapView(cuz) + mapview(abc). Coloring options are available with zcol option.
I have saved .gpkg files containing the matched cells (separate files for CEM and PSM) in the following directory: "../../datalake/mapme.protectedareas/output/polygon/matched_cells/"
I created a markdown file with maps for the years 2008 to 2011 (both CEM and PSM). As mentioned by Johannes, I have converted the polygons to points for the maps in the markdown files. Haven't pushed the file to Github yet but I'll post the link here as soon as the file is pushed.
Screenshot from QGIS with CEM (left) and PSM (right):
Probably only works mapping-wise for years with limited number and control cells.
Hint: you can create also centroids for the polygons which are point data and much more efficient to plot. This might also allow you to plot several thousand observations in one map.use the
sf
package andst_centroid
function to convert polygons to points, if necessary.Definition of Done:
write_sf
) for later mapping in QGISmapView
. You can add several years to the map by concetenating several mapviews -> mapView(cuz) + mapview(abc). Coloring options are available withzcol
option.