Skip to content

Commit

Permalink
Adding the variable "classname" to deter function. (#245)
Browse files Browse the repository at this point in the history
  • Loading branch information
carolinamoura2000 authored Apr 17, 2024
1 parent 5aae3d2 commit f61461f
Showing 1 changed file with 30 additions and 29 deletions.
59 changes: 30 additions & 29 deletions R/deter.R
Original file line number Diff line number Diff line change
Expand Up @@ -23,43 +23,42 @@ load_deter <- function(dataset, raw_data = FALSE,
###########################
## Bind Global Variables ##
###########################

.data <- view_date <- name_muni <- code_muni <- sensor <- satellite <- abbrev_state <- NULL
uc <- classname <- path_row <- area <- quadrant <- geometry <- id_alerta <- NULL

#############################
## Define Basic Parameters ##
#############################

param <- list()
param$source <- "deter"
param$dataset <- dataset
param$language <- language
param$raw_data <- raw_data

# check if dataset is valid

check_params(param)



#################
## Downloading ##
#################

dat <- external_download(
dataset = param$dataset,
source = param$source
)

## Return Raw Data

if (param$raw_data) {
return(dat)
}

######################
## Data Engineering ##
######################

dat <- dat %>%
janitor::clean_names() %>%
dplyr::mutate(
Expand All @@ -68,70 +67,72 @@ load_deter <- function(dataset, raw_data = FALSE,
~ stringi::stri_trans_general(., id = "Latin-ASCII")
)
)

# Loading municipal map data
geo_br <- external_download(
dataset = "geo_municipalities",
source = "internal"
)

# Adding alert_id variable to preserve the information of which rows belong to the same alert

dat <- dat %>%
dplyr::mutate(id_alerta = dplyr::row_number())

###################
## Harmonize CRS ##
###################

# The crs that will be used to overlap maps below
operation_crs <- sf::st_crs("+proj=poly +lat_0=0 +lon_0=-54 +x_0=5000000 +y_0=10000000 +ellps=aust_SA +units=m +no_defs")

# Changing crs of both data to the common crs chosen above
dat$geometry <- sf::st_make_valid(sf::st_transform(dat$geometry, operation_crs))
geo_br$geom <- sf::st_transform(geo_br$geom, operation_crs)

# Overlaps shapefiles
sf::st_geometry(dat) <- dat$geometry
sf::st_geometry(geo_br) <- geo_br$geom

dat <- suppressWarnings(sf::st_intersection(dat, geo_br)) %>%
dplyr::mutate(area = sf::st_area(.data$geometry))

dat <- dat %>%
dplyr::select(
view_date, name_muni, code_muni, abbrev_state,
view_date, name_muni, code_muni, abbrev_state, classname,
area, geometry, id_alerta
)

###################
## Renaming Data ##
###################

if (param$language == "pt") {
dat_mod <- dat %>%
dplyr::rename(
"data" = view_date,
"municipio" = name_muni,
"cod_municipio" = code_muni,
"uf" = abbrev_state
"uf" = abbrev_state,
"tipo_de_alerta" = classname
)
}

if (param$language == "eng") {
dat_mod <- dat %>%
dplyr::rename(
"date" = view_date,
"municipality" = name_muni,
"municipality_code" = code_muni,
"state" = abbrev_state,
"alert_id" = id_alerta
"alert_id" = id_alerta,
"alert_type" = classname
)
}

#################
## Return Data ##
#################

return(dat_mod)
}

0 comments on commit f61461f

Please sign in to comment.