Vamos a importar un archivo .csv de texto que contiene datos de temperatura promedio y indice NDVI (vegetacion) para las 720 colonias de hermosillo par el mes de Julio de 2015, datos obtenidos de Landsat 8 de los sensores TIRS (termico) y OLI (multiespectral)
setwd("~/Documents/ICIAM2021")
library(readr)
hermosillo <- read_csv("hermosillo.csv")
##
## ── Column specification ────────────────────────────────────────────────────────
## cols(
## NDVI = col_double(),
## TEMPERATURE = col_double()
## )
head(hermosillo)
## # A tibble: 6 x 2
## NDVI TEMPERATURE
## <dbl> <dbl>
## 1 0.236 28.1
## 2 0.196 28.3
## 3 0.241 26.6
## 4 0.460 26.6
## 5 0.194 27.3
## 6 0.151 28.0
cor(hermosillo)
## NDVI TEMPERATURE
## NDVI 1.0000000 -0.5731116
## TEMPERATURE -0.5731116 1.0000000
pairs(hermosillo)