Working, but slow implementation of dhb boundaries

Note layer issues!

GeoProcessing.R

@@ -0,0 +1,9 @@
+td <- tempdir()
+
+unzip("statsnzdistrict-health-board-2015-SHP.zip", exdir = td)
+dhb2015 <- readOGR(td, layer="district-health-board-2015")
+dhbsub <- subset(dhb2015, dhb2015$DHB2015_Co != "99")
+#dhbsub <- gSimplify(dhbsub, tol = 0.01, topologyPreserve = TRUE)
+#
+#saveRDS(dhbsub, "MedicalCentres/dhb2015.rds")
+saveRDS(dhb2015, "MedicalCentres/dhb2015.rds")

MedicalCentres/app.R

@@ -4,6 +4,9 @@ library(leaflet)
 library(dplyr)
 library(shinyWidgets)
 library(htmltools)
+library(rgdal)
+library(rgeos)
+library(ggfortify)
 
 leafcols <- c("blue", "red", "darkred", "lightred",
           "orange", "beige", "green", "darkgreen",
@@ -33,6 +36,7 @@ snargcol <- function(snarg) {
     sdf
 }
 
+dhb2015 <- readRDS("dhb2015.rds")
 
 popups <- function(data) {
 sprintf("<strong>%s</strong><br/>%s<br/>DHB: %s<br/>%s",
@@ -90,12 +94,16 @@ ui <- fillPage(title = "New Zealand GP Practice Locations 2017",
 
 server <- function(input, output, session) {
     output$medmap <- renderLeaflet({
-        gl <- leaflet(Geo2017)
+        gl <- leaflet(Geo2017) %>%
+          addMapPane("polygons", zIndex = 420) %>%
+          addMapPane("markers", zIndex = 430)
         gl <- addTiles(gl)
         gl <- addCircleMarkers(gl, lng = Geo2017$lon,
                    lat = Geo2017$lat,
                    popup = popups(Geo2017),
-                   col= "blue")
+                   col= "blue",
+                   #pathOptions(pane = "markers")
+                   )
         gl
     })
 
@@ -111,20 +119,24 @@ server <- function(input, output, session) {
                      lat = Geo2017$lat,
                      popup = popups(Geo2017),
                      # color = snargcols[as.character(Geo2017$dhb_name)]
-                     col = snargcol(Geo2017$Snargle)
+                     col = snargcol(Geo2017$Snargle),
+                     #pathOptions(pane = "markers")
                      )
       } else if (input$selcol == "DHB") {
           gl <- addCircleMarkers(gl, lng = Geo2017$lon,
                      lat = Geo2017$lat,
                      popup = popups(Geo2017),
                      # color = dhbcols[Geo2017$dhb_name]
-                     col = dhbcol(Geo2017$dhb_name)
+                     col = dhbcol(Geo2017$dhb_name),
+                     #pathOptions(pane = "markers")
                      )
       } else {
           gl <- addCircleMarkers(gl, lng = Geo2017$lon,
                      lat = Geo2017$lat,
                      popup = popups(Geo2017),
-                     col= "blue")
+                     col= "blue",
+                     #pathOptions(pane = "markers")
+                     )
       }
     })
 
@@ -143,6 +155,17 @@ server <- function(input, output, session) {
                                          zoom = 13)
       }
     })
+    observeEvent(input$dhboverlay, {
+      if (input$dhboverlay) {
+        leafletProxy("medmap") %>% addPolygons(data = dhbsub, fill = FALSE,
+                                               color = "#000000",
+                                               weight = 2,
+                                               #pathOptions(pane = "polygons")
+                                               )
+      } else {
+        leafletProxy("medmap") %>% clearShapes()
+      }
+    })
 
 }