A Energia Radiativa do Fogo (FRP) é uma técnica para quantificar a biomassa queimada usando dados de sensoriamento remoto.A FRP mede a energia radiante emitida por unidade de tempo pela vegetação queimada.
A estimativa do total de biomassa queimada é realizada a partir da Energia Radiativa do Fogo (FRE), que é definida como a energia emitida pelo fogo como radiação eletromagnética durante o ciclo de vida da queimada e pode ser obtida a partir da integração temporal da FRP (WOOSTER et al., 2003). interessante:https://revistapesquisa.fapesp.br/como-monitorar-o-fogo/
dados_queimadas_v1 <- read.csv("Focos_2021-06-26_2021-06-27.csv")
View(dados_queimadas_v1)
str(dados_queimadas_v1)
## 'data.frame': 627 obs. of 12 variables:
## $ datahora : chr "2021/06/26 17:00:00" "2021/06/26 17:00:00" "2021/06/26 17:00:00" "2021/06/26 17:00:00" ...
## $ satelite : chr "AQUA_M-T" "AQUA_M-T" "AQUA_M-T" "AQUA_M-T" ...
## $ pais : chr "Brasil" "Brasil" "Brasil" "Brasil" ...
## $ estado : chr "RIO DE JANEIRO" "SAO PAULO" "MATO GROSSO DO SUL" "SAO PAULO" ...
## $ municipio : chr "SANTO ANTONIO DE PADUA" "IBITINGA" "ANTONIO JOAO" "SANTA BARBARA D'OESTE" ...
## $ bioma : chr "Mata Atlantica" "Cerrado" "Cerrado" "Mata Atlantica" ...
## $ diasemchuva : int 17 15 9 4 -999 3 3 3 3 2 ...
## $ precipitacao: num 0 0 0 0.1 0 0.3 0.3 0.3 0.3 0 ...
## $ riscofogo : num -999 1 0.5 1 1 1 1 1 0.9 1 ...
## $ latitude : num -21.5 -21.9 -22.3 -22.8 -22.9 ...
## $ longitude : num -42.2 -48.8 -56 -47.5 -43.4 ...
## $ frp : num 11.4 4.7 28.9 15.3 9.3 ...
summary(dados_queimadas_v1)
## datahora satelite pais estado
## Length:627 Length:627 Length:627 Length:627
## Class :character Class :character Class :character Class :character
## Mode :character Mode :character Mode :character Mode :character
##
##
##
##
## municipio bioma diasemchuva precipitacao
## Length:627 Length:627 Min. :-999.0 Min. :0.00000
## Class :character Class :character 1st Qu.: 14.0 1st Qu.:0.00000
## Mode :character Mode :character Median : 24.0 Median :0.00000
## Mean : 21.9 Mean :0.06109
## 3rd Qu.: 37.0 3rd Qu.:0.00000
## Max. : 59.0 Max. :3.20000
## NA's :185 NA's :185
## riscofogo latitude longitude frp
## Min. :-999.00 Min. :-25.865 Min. :-72.32 Min. : 4.60
## 1st Qu.: 0.80 1st Qu.:-14.695 1st Qu.:-53.85 1st Qu.: 14.45
## Median : 1.00 Median :-11.747 Median :-49.72 Median : 24.20
## Mean : -10.43 Mean :-12.304 Mean :-50.26 Mean : 74.76
## 3rd Qu.: 1.00 3rd Qu.: -9.633 3rd Qu.:-46.94 3rd Qu.: 54.80
## Max. : 1.00 Max. : 3.875 Max. :-35.59 Max. :5013.60
## NA's :185
dados_queimadas_v1[1:5,]
## datahora satelite pais estado municipio
## 1 2021/06/26 17:00:00 AQUA_M-T Brasil RIO DE JANEIRO SANTO ANTONIO DE PADUA
## 2 2021/06/26 17:00:00 AQUA_M-T Brasil SAO PAULO IBITINGA
## 3 2021/06/26 17:00:00 AQUA_M-T Brasil MATO GROSSO DO SUL ANTONIO JOAO
## 4 2021/06/26 17:00:00 AQUA_M-T Brasil SAO PAULO SANTA BARBARA D'OESTE
## 5 2021/06/26 17:00:00 AQUA_M-T Brasil RIO DE JANEIRO RIO DE JANEIRO
## bioma diasemchuva precipitacao riscofogo latitude longitude frp
## 1 Mata Atlantica 17 0.0 -999.0 -21.522 -42.201 11.4
## 2 Cerrado 15 0.0 1.0 -21.870 -48.799 4.7
## 3 Cerrado 9 0.0 0.5 -22.263 -56.016 28.9
## 4 Mata Atlantica 4 0.1 1.0 -22.844 -47.454 15.3
## 5 Mata Atlantica -999 0.0 1.0 -22.851 -43.442 9.3
Verificamos que existem dados faltantes e outliers
Análise dias sem chuva
hist(dados_queimadas_v1$diasemchuva,
main = "Histograma dias sem chuva",
xlab = "Quantidade de Dias",
ylab = "Frequência",
col = "#6600cc",
labels = TRUE)
#Removendo outliers negativos
dados_queimadas_v2 =filter(dados_queimadas_v1, diasemchuva >=0 | is.na(diasemchuva),riscofogo >=0 | is.na(riscofogo))
summary(dados_queimadas_v2)
## datahora satelite pais estado
## Length:621 Length:621 Length:621 Length:621
## Class :character Class :character Class :character Class :character
## Mode :character Mode :character Mode :character Mode :character
##
##
##
##
## municipio bioma diasemchuva precipitacao
## Length:621 Length:621 Min. : 0.00 Min. :0.00000
## Class :character Class :character 1st Qu.:14.00 1st Qu.:0.00000
## Mode :character Mode :character Median :24.00 Median :0.00000
## Mean :24.27 Mean :0.06193
## 3rd Qu.:37.00 3rd Qu.:0.00000
## Max. :59.00 Max. :3.20000
## NA's :185 NA's :185
## riscofogo latitude longitude frp
## Min. :0.0000 Min. :-25.865 Min. :-72.32 Min. : 4.60
## 1st Qu.:0.8000 1st Qu.:-14.645 1st Qu.:-54.00 1st Qu.: 14.50
## Median :1.0000 Median :-11.721 Median :-49.77 Median : 24.30
## Mean :0.8782 Mean :-12.241 Mean :-50.30 Mean : 75.34
## 3rd Qu.:1.0000 3rd Qu.: -9.568 3rd Qu.:-46.94 3rd Qu.: 55.00
## Max. :1.0000 Max. : 3.875 Max. :-35.59 Max. :5013.60
## NA's :185
View(dados_queimadas_v1)
View(dados_queimadas_v2)
Note that the echo = FALSE parameter was added to the code chunk to prevent printing of the R code that generated the plot.