Este material está disponível em: http://rpubs.com/leonardoreffatti.
library(leaflet)
library(leaflet.extras)
setwd("C:/R/heatmap")
data <- read.table("teste3.txt", h=T)
attach(data)
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)