RikFerreira

library(tidyverse)
library(sf)
library(units)

ponto <- st_point(c(1, 2)) %>%
    st_sfc(crs = 31985) %>%
    st_sf()

offset <- 5 # 5 unidades de mapa (EPSG 31985: metros)

RikFerreira

RikFerreira

RikFerreira

RikFerreira

# Source: https://minato.sip21c.org/swtips/factor-in-R.pdf

bartlett.sphericity.test <- function(x) {
    method <- "Bartlett’s test of sphericity"
    
    data.name <- deparse(substitute(x))
    
    x <- subset(x, complete.cases(x)) # Omit missing values
    n <- nrow(x)
    p <- ncol(x)

RikFerreira

library(tidyverse)

dataset <- tibble(
  position = 0:9,
  x0 = 1:10,
  y0 = 1:10,
  x1 = 1:10,
  y1 = 1:10
)

RikFerreira

library(raster)
library(tidyverse)
library(sf)
library(caret)
library(RStoolbox)

source("scripts/functions.R")

# Generate cropped landsat ----------------------------------------------------

RikFerreira

library(tidyverse)
library(sidrar)
library(survey)
library(srvyr)
library(PNADcIBGE)
library(janitor)

info_sidra(1288)

pop_censos <- get_sidra(

RikFerreira

library(tidyverse)
library(randomForest)

train <- data.frame(
    x = 1:10,
    y = 1:10 * 100
)

lm <- lm(y ~ x, data = train)
rf <- randomForest(y ~ x, data = train)

RikFerreira

# diff LT05_L1TP_218072_20100716_20161014_01_T1_MTL.txt LT05_L1TP_218072_20100801_20161015_01_T1_MTL.txt

4,6c4,6
<     REQUEST_ID = "0501610125656_06672"
<     LANDSAT_SCENE_ID = "LT52180722010197COA00"
<     LANDSAT_PRODUCT_ID = "LT05_L1TP_218072_20100716_20161014_01_T1"
---
>     REQUEST_ID = "0501610125560_19672"
>     LANDSAT_SCENE_ID = "LT52180722010213CUB00"
>     LANDSAT_PRODUCT_ID = "LT05_L1TP_218072_20100801_20161015_01_T1"