library(leaflet.extras)
library(dplyr)
library(plyr)
## Warning: package 'plyr' was built under R version 3.6.3
library(data.table)
library(tidyverse)
## Warning: package 'tidyverse' was built under R version 3.6.3
## Warning: package 'forcats' was built under R version 3.6.3
#setwd("C:\\xampp\\htdocs\\fire\\2019")
file<-read.csv("fire_archive_M6_98797.csv",header=TRUE)
ct01<-filter(file,acq_date %like% '2019-01')
ct02<-filter(file,acq_date %like% '2019-02')
ct03<-filter(file,acq_date %like% '2019-03')
ct04<-filter(file,acq_date %like% '2019-04')
ct05<-filter(file,acq_date %like% '2019-05')
ct06<-filter(file,acq_date %like% '2019-06')
ct07<-filter(file,acq_date %like% '2019-07')
ct08<-filter(file,acq_date %like% '2019-08')
ct09<-filter(file,acq_date %like% '2019-09')
ct10<-filter(file,acq_date %like% '2019-10')
ct11<-filter(file,acq_date %like% '2019-11')
ct12<-filter(file,acq_date %like% '2019-12')
leaflet() %>%
addTiles(group = "OSM (default)") %>%
addProviderTiles(providers$Esri.WorldImagery, group = "Esri.WorldImagery") %>%
setView( -8.8018 , 37.9565, 6 ) %>%
addHeatmap(data= ct01,lng = ~longitude, lat = ~latitude, intensity = ~frp,
blur = 20, max = 0.05, radius = 15,group="2019-01")%>%
addHeatmap(data= ct02,lng = ~longitude, lat = ~latitude, intensity = ~frp,
blur = 20, max = 0.05, radius = 15,group="2019-02")%>%
addHeatmap(data= ct03,lng = ~longitude, lat = ~latitude, intensity = ~frp,
blur = 20, max = 0.05, radius = 15,group="2019-03")%>%
addHeatmap(data= ct04,lng = ~longitude, lat = ~latitude, intensity = ~frp,
blur = 20, max = 0.05, radius = 15,group="2019-04")%>%
addHeatmap(data= ct05,lng = ~longitude, lat = ~latitude, intensity = ~frp,
blur = 20, max = 0.05, radius = 15,group="2019-05")%>%
addHeatmap(data= ct06,lng = ~longitude, lat = ~latitude, intensity = ~frp,
blur = 20, max = 0.05, radius = 15,group="2019-06")%>%
addHeatmap(data= ct07,lng = ~longitude, lat = ~latitude, intensity = ~frp,
blur = 20, max = 0.05, radius = 15,group="2019-07")%>%
addHeatmap(data= ct08,lng = ~longitude, lat = ~latitude, intensity = ~frp,
blur = 20, max = 0.05, radius = 15,group="2019-08")%>%
addHeatmap(data= ct09,lng = ~longitude, lat = ~latitude, intensity = ~frp,
blur = 20, max = 0.05, radius = 15,group="2019-09")%>%
addHeatmap(data= ct10,lng = ~longitude, lat = ~latitude, intensity = ~frp,
blur = 20, max = 0.05, radius = 15,group="2019-10")%>%
addHeatmap(data= ct11,lng = ~longitude, lat = ~latitude, intensity = ~frp,
blur = 20, max = 0.05, radius = 15,group="2019-11")%>%
addHeatmap(data= ct12,lng = ~longitude, lat = ~latitude, intensity = ~frp,
blur = 20, max = 0.05, radius = 15,group="2019-12")%>%
addLayersControl(
baseGroups = c("OSM (default)", "Esri.WorldImagery"),
overlayGroups = c("2019-01", "2019-02" ,"2019-03","2019-04","2019-05","2019-06",
"2019-07", "2019-08" ,"2019-09","2019-10","2019-11","2019-12"),
options = layersControlOptions(collapsed = FALSE)
)
Font:https://earthdata.nasa.gov/earth-observation-data/near-real-time/firms