Clase Electiva II

Codigo html:

#instalamos paquetes raster, sp y rgdal.
#install.packages("raster")
#install.packages("rgdal")

#los llamamos.
library("raster")

## Loading required package: sp

library("rgdal")

## Please note that rgdal will be retired by the end of 2023,
## plan transition to sf/stars/terra functions using GDAL and PROJ
## at your earliest convenience.
## 
## rgdal: version: 1.5-25, (SVN revision 1143)
## Geospatial Data Abstraction Library extensions to R successfully loaded
## Loaded GDAL runtime: GDAL 3.2.1, released 2020/12/29
## Path to GDAL shared files: C:/Users/felip/Documents/R/win-library/4.1/rgdal/gdal
## GDAL binary built with GEOS: TRUE 
## Loaded PROJ runtime: Rel. 7.2.1, January 1st, 2021, [PJ_VERSION: 721]
## Path to PROJ shared files: C:/Users/felip/Documents/R/win-library/4.1/rgdal/proj
## PROJ CDN enabled: FALSE
## Linking to sp version:1.4-5
## To mute warnings of possible GDAL/OSR exportToProj4() degradation,
## use options("rgdal_show_exportToProj4_warnings"="none") before loading sp or rgdal.
## Overwritten PROJ_LIB was C:/Users/felip/Documents/R/win-library/4.1/rgdal/proj

#asignamos las variables.
b1 <-raster("C:/Users/felip/Desktop/electiva2/practica2/LC08_L2SP_217064_20210723_20210729_02_T1_SR_B1.tif")
b2 <-raster("C:/Users/felip/Desktop/electiva2/practica2/LC08_L2SP_217064_20210723_20210729_02_T1_SR_B2.tif")
b3 <-raster("C:/Users/felip/Desktop/electiva2/practica2/LC08_L2SP_217064_20210723_20210729_02_T1_SR_B3.tif")
b4 <-raster("C:/Users/felip/Desktop/electiva2/practica2/LC08_L2SP_217064_20210723_20210729_02_T1_SR_B4.tif")
b5 <-raster("C:/Users/felip/Desktop/electiva2/practica2/LC08_L2SP_217064_20210723_20210729_02_T1_SR_B5.tif")
b6 <-raster("C:/Users/felip/Desktop/electiva2/practica2/LC08_L2SP_217064_20210723_20210729_02_T1_SR_B6.tif")
b7 <-raster("C:/Users/felip/Desktop/electiva2/practica2/LC08_L2SP_217064_20210723_20210729_02_T1_SR_B7.tif")

#probamos los datos de las bandas.
#banda b1 hasta la b7

#sistema de referencia
crs(b1)

## CRS arguments:
##  +proj=utm +zone=24 +datum=WGS84 +units=m +no_defs

crs(b2)

## CRS arguments:
##  +proj=utm +zone=24 +datum=WGS84 +units=m +no_defs

crs(b3)

## CRS arguments:
##  +proj=utm +zone=24 +datum=WGS84 +units=m +no_defs

crs(b4)

## CRS arguments:
##  +proj=utm +zone=24 +datum=WGS84 +units=m +no_defs

crs(b5)

## CRS arguments:
##  +proj=utm +zone=24 +datum=WGS84 +units=m +no_defs

crs(b6)

## CRS arguments:
##  +proj=utm +zone=24 +datum=WGS84 +units=m +no_defs

crs(b7)

## CRS arguments:
##  +proj=utm +zone=24 +datum=WGS84 +units=m +no_defs

#tambien podemos ver las filas por culumnas #dim(#banda)

#comparacion de bandas raster
compareRaster(b1,b7)

## [1] TRUE

#creamos raster (objeto con varias capas) (stack)
s <- stack(b4, b3, b2, b1)


#ploteo de bandas
#probamos para generar la banda
par(mfrow= c(3,3))
plot(b1, main="coastal",col =gray (0:100/100))
plot(b2, main="blue",col =gray (0:100/100))
plot(b3, main="green",col =gray (0:100/100))
plot(b4, main="red",col =gray (0:100/100))
plot(b5, main="Near infrared",col =gray (0:100/100))
plot(b6, main="Short infrared",col =gray (0:100/100))
plot(b7, main="Short infrared",col =gray (0:100/100))


#ploteamos a un color verdadero

landsatRGB <- stack (b7, b6, b5, b4 ,b3, b2, b1)

plotRGB(landsatRGB, axes = TRUE, stretch = "lin", main = "color verdadero")

Practica2

Andres felipe Rivas

14/9/2021

UNIVERSIDAD COOPEARTIVA DE COLOMBIA

Facultad de Ingenieria

practica2

Clase Electiva II