Heatmap produção de Uva

Este material está disponível em: http://rpubs.com/leonardoreffatti.

Pacotes necessários e leitura dos dados:

library(leaflet)
library(leaflet.extras)
setwd("C:/R/heatmap")
data <- read.table("teste3.txt", h=T)
attach(data)

Heatmap:

data <- subset(data, !is.na(lat))
leaflet(data) %>%
  addProviderTiles(providers$Esri.WorldImagery) %>%
  addHeatmap(group="heat", lng=~lon, lat=~lat, blur=40, radius = 30, 
             intensity = ~producao_2/1000, max=10)