InteractiveMaps.RmdInteractive maps can be used to visually inspect the polygons. Although there is no way to directly edit the polygon geometry from them, it makes it easy to to quickly assess polygons without the need to load the data into GIS software such as ArcMap.
This tutorial explains how interactive maps can be generated within the polyCheck package.
The following analysis requires several packages to be loaded:
library(tidyverse)
library(raster)
library(polyCheck)
library(rgeos)As before, we must load the sample data:
# Load Data
dataPath <- system.file("extdata", package = "polyCheck")
polygons <- list_polygons(dataPath, "polygons")
survey_points <- shapefile(file.path(dataPath, "survey_points_sample.shp"))
settlement_points <- shapefile(file.path(dataPath, "settlement_points.shp"))The function display_polygons can be used to display the polygon. The map can be used to change the zoom and pan around the map to visually inspect the polygon boundaries. Note, that if the zoom level is too low, the basemap may not be displayed.
# Just displays the polygon
display_polygon(polygons[2])The maps can be made more informative by supplying the survey points dataset. The function will automatically find the point which matches the ID of the polygon, and plot both. This is demonstrated below, and helps in identifying whether the polygon is plotted in the right context.
# This is a slight workaround to make the code loop work
maps <- display_polygon(polygons[3], survey_points = survey_points)
htmltools::tagList(maps)# This is a slight workaround to make the code loop work
maps <- display_polygon(polygons[3], survey_points, settlement_points)
htmltools::tagList(maps)