factpal <- colorFactor(topo.colors(15), stands_prj$Forest)
pal <- colorFactor(palette = "orange",
                    domain = plots_prj$BasalArea)

leaflet(plots_prj) %>% addTiles() %>%  
addPolygons(data=stands_prj,weight=0,fillOpacity = 0.5, smoothFactor = 0.5, color = ~factpal(stands_prj$Forest),
              label = stands_prj$Forest) %>%

addCircles(lng = ~Long, lat = ~Lat, weight = 3,
    radius = ~ BasalArea ,label= ~Plot, fillColor = ~pal(BasalArea),
                 stroke = TRUE, fillOpacity = 0.5)