PRACTICA 05
JHON JAIRO CUENCA Y ANGIE DANNIELLE ZEDIEL
24/4/2021
#install.packages("leaflet")
#install.packages("Rcpp")
library(leaflet)## Warning: package 'leaflet' was built under R version 4.0.5longitude_astudy <- c(-75.194022, -75.145162, -75.098085, -75.099645, -75.182693)
latitud_astudy <- c(3.303722, 3.302725, 3.272956, 3.209103, 3.224697)
#plot(longitude_astudy, latitud_astudy)
map_one <- leaflet() %>%
addTiles() %>%
addMarkers(lng = -75.144727, lat = 3.255859)
addPolygons(map_one,longitude_astudy,latitud_astudy)library(raster)## Loading required package: splibrary(sp)
library(rgdal)## rgdal: version: 1.5-23, (SVN revision 1121)
## 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/Jairo Cuenca/Documents/R/win-library/4.0/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/Jairo Cuenca/Documents/R/win-library/4.0/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 rgdal.
## Overwritten PROJ_LIB was C:/Users/Jairo Cuenca/Documents/R/win-library/4.0/rgdal/projbandas <- paste0("C:/Users/Jairo Cuenca/Downloads/PRACTICA 05/LC08_L2SP_008058_20210410_20210416_02_T1_SR_B",1:7,".tif")
bandas## [1] "C:/Users/Jairo Cuenca/Downloads/PRACTICA 05/LC08_L2SP_008058_20210410_20210416_02_T1_SR_B1.tif"
## [2] "C:/Users/Jairo Cuenca/Downloads/PRACTICA 05/LC08_L2SP_008058_20210410_20210416_02_T1_SR_B2.tif"
## [3] "C:/Users/Jairo Cuenca/Downloads/PRACTICA 05/LC08_L2SP_008058_20210410_20210416_02_T1_SR_B3.tif"
## [4] "C:/Users/Jairo Cuenca/Downloads/PRACTICA 05/LC08_L2SP_008058_20210410_20210416_02_T1_SR_B4.tif"
## [5] "C:/Users/Jairo Cuenca/Downloads/PRACTICA 05/LC08_L2SP_008058_20210410_20210416_02_T1_SR_B5.tif"
## [6] "C:/Users/Jairo Cuenca/Downloads/PRACTICA 05/LC08_L2SP_008058_20210410_20210416_02_T1_SR_B6.tif"
## [7] "C:/Users/Jairo Cuenca/Downloads/PRACTICA 05/LC08_L2SP_008058_20210410_20210416_02_T1_SR_B7.tif"landsat <- stack(bandas)
landsat## class : RasterStack
## dimensions : 7731, 7591, 58686021, 7 (nrow, ncol, ncell, nlayers)
## resolution : 30, 30 (x, y)
## extent : 410385, 638115, 203385, 435315 (xmin, xmax, ymin, ymax)
## crs : +proj=utm +zone=18 +datum=WGS84 +units=m +no_defs
## names : LC08_L2SP//2_T1_SR_B1, LC08_L2SP//2_T1_SR_B2, LC08_L2SP//2_T1_SR_B3, LC08_L2SP//2_T1_SR_B4, LC08_L2SP//2_T1_SR_B5, LC08_L2SP//2_T1_SR_B6, LC08_L2SP//2_T1_SR_B7
## min values : 0, 0, 0, 0, 0, 0, 0
## max values : 65535, 65535, 65535, 65535, 65535, 65535, 65535#install.packages("readxl")
#install.packages("sf")
library(readxl)## Warning: package 'readxl' was built under R version 4.0.5library(sf)## Warning: package 'sf' was built under R version 4.0.5## Linking to GEOS 3.9.0, GDAL 3.2.1, PROJ 7.2.1# llamar al archivo de excel
PUNTOS <- read_xlsx("C:/Users/Jairo Cuenca/Downloads/PRACTICA 05/COORDENADAS.xlsx")
#PUNTOS
#pasar de excel a BDE
PUNTOS_geo <- PUNTOS %>%
st_as_sf(coords = c("X","Y"))
plot(PUNTOS_geo)
DF1 <- extract(landsat,PUNTOS_geo)
DF1## LC08_L2SP_008058_20210410_20210416_02_T1_SR_B1
## [1,] 9376
## [2,] 8717
## [3,] 8843
## [4,] 9208
## [5,] 8805
## [6,] 10494
## [7,] 13441
## [8,] 9682
## [9,] 15290
## [10,] 12098
## [11,] 8948
## [12,] 15182
## [13,] 20215
## [14,] 9331
## [15,] 11671
## [16,] 8851
## [17,] 12300
## [18,] 13568
## [19,] 9362
## [20,] 16356
## [21,] 9129
## [22,] 10465
## [23,] 9898
## [24,] 10939
## [25,] 8728
## [26,] 11418
## [27,] 8493
## [28,] 9351
## [29,] 10243
## [30,] 10377
## [31,] 10103
## [32,] 7857
## [33,] 9421
## [34,] 9572
## [35,] 8918
## [36,] 9145
## [37,] 9159
## [38,] 8104
## [39,] 14172
## [40,] 10678
## [41,] 7631
## [42,] 8860
## [43,] 8040
## [44,] 8956
## [45,] 8259
## [46,] 8487
## [47,] 9845
## [48,] 14658
## [49,] 14524
## [50,] 9228
## [51,] 15519
## [52,] 9134
## LC08_L2SP_008058_20210410_20210416_02_T1_SR_B2
## [1,] 9965
## [2,] 9187
## [3,] 9268
## [4,] 9704
## [5,] 8864
## [6,] 10870
## [7,] 13563
## [8,] 10004
## [9,] 15485
## [10,] 12844
## [11,] 9266
## [12,] 15349
## [13,] 20069
## [14,] 10347
## [15,] 11940
## [16,] 9193
## [17,] 12322
## [18,] 13130
## [19,] 9663
## [20,] 16160
## [21,] 9513
## [22,] 11236
## [23,] 10424
## [24,] 11474
## [25,] 9160
## [26,] 11487
## [27,] 8843
## [28,] 9501
## [29,] 11008
## [30,] 10897
## [31,] 10688
## [32,] 8363
## [33,] 9614
## [34,] 9726
## [35,] 9376
## [36,] 9571
## [37,] 9793
## [38,] 9456
## [39,] 14381
## [40,] 11767
## [41,] 7590
## [42,] 9343
## [43,] 8159
## [44,] 8890
## [45,] 8508
## [46,] 8734
## [47,] 10318
## [48,] 14564
## [49,] 14518
## [50,] 9606
## [51,] 15523
## [52,] 10054
## LC08_L2SP_008058_20210410_20210416_02_T1_SR_B3
## [1,] 11824
## [2,] 10320
## [3,] 10576
## [4,] 10965
## [5,] 10383
## [6,] 13245
## [7,] 14937
## [8,] 11755
## [9,] 16252
## [10,] 14313
## [11,] 11482
## [12,] 16740
## [13,] 21557
## [14,] 14219
## [15,] 13750
## [16,] 10689
## [17,] 15067
## [18,] 16430
## [19,] 11133
## [20,] 19445
## [21,] 11113
## [22,] 13188
## [23,] 12449
## [24,] 13621
## [25,] 10801
## [26,] 12378
## [27,] 10369
## [28,] 11756
## [29,] 13522
## [30,] 13044
## [31,] 12865
## [32,] 10478
## [33,] 11527
## [34,] 11473
## [35,] 11173
## [36,] 11032
## [37,] 13055
## [38,] 12358
## [39,] 15593
## [40,] 14105
## [41,] 8041
## [42,] 10488
## [43,] 9689
## [44,] 10116
## [45,] 9627
## [46,] 9643
## [47,] 11863
## [48,] 15982
## [49,] 16696
## [50,] 12294
## [51,] 17355
## [52,] 11614
## LC08_L2SP_008058_20210410_20210416_02_T1_SR_B4
## [1,] 11922
## [2,] 10092
## [3,] 10452
## [4,] 10685
## [5,] 9887
## [6,] 12800
## [7,] 15355
## [8,] 12396
## [9,] 16755
## [10,] 14801
## [11,] 12147
## [12,] 16565
## [13,] 21348
## [14,] 14210
## [15,] 14356
## [16,] 10454
## [17,] 15333
## [18,] 15960
## [19,] 10889
## [20,] 19165
## [21,] 10938
## [22,] 13454
## [23,] 13211
## [24,] 14729
## [25,] 10805
## [26,] 12594
## [27,] 10929
## [28,] 12361
## [29,] 13889
## [30,] 14417
## [31,] 13904
## [32,] 11774
## [33,] 11228
## [34,] 11223
## [35,] 10909
## [36,] 10923
## [37,] 13190
## [38,] 12366
## [39,] 16061
## [40,] 14606
## [41,] 7525
## [42,] 9866
## [43,] 8551
## [44,] 9048
## [45,] 9051
## [46,] 9266
## [47,] 11598
## [48,] 16194
## [49,] 16857
## [50,] 12802
## [51,] 17548
## [52,] 11448
## LC08_L2SP_008058_20210410_20210416_02_T1_SR_B5
## [1,] 16929
## [2,] 15171
## [3,] 17365
## [4,] 15301
## [5,] 12847
## [6,] 18426
## [7,] 21565
## [8,] 18874
## [9,] 22504
## [10,] 20876
## [11,] 19587
## [12,] 22026
## [13,] 26077
## [14,] 19744
## [15,] 20276
## [16,] 19466
## [17,] 21301
## [18,] 21487
## [19,] 17410
## [20,] 24291
## [21,] 18393
## [22,] 19879
## [23,] 19848
## [24,] 19479
## [25,] 16373
## [26,] 15546
## [27,] 15907
## [28,] 19581
## [29,] 19906
## [30,] 18943
## [31,] 18101
## [32,] 15091
## [33,] 15952
## [34,] 15336
## [35,] 14390
## [36,] 14025
## [37,] 18956
## [38,] 16893
## [39,] 21510
## [40,] 19470
## [41,] 12338
## [42,] 22114
## [43,] 24421
## [44,] 24124
## [45,] 18638
## [46,] 15364
## [47,] 19509
## [48,] 22711
## [49,] 22560
## [50,] 20588
## [51,] 22911
## [52,] 18111
## LC08_L2SP_008058_20210410_20210416_02_T1_SR_B6
## [1,] 16166
## [2,] 13853
## [3,] 14345
## [4,] 14574
## [5,] 12825
## [6,] 16816
## [7,] 19141
## [8,] 17289
## [9,] 19998
## [10,] 19355
## [11,] 17677
## [12,] 19689
## [13,] 23372
## [14,] 18302
## [15,] 18766
## [16,] 16177
## [17,] 19818
## [18,] 19842
## [19,] 15702
## [20,] 22807
## [21,] 16192
## [22,] 19020
## [23,] 17806
## [24,] 18475
## [25,] 14441
## [26,] 15183
## [27,] 15028
## [28,] 18044
## [29,] 18540
## [30,] 18194
## [31,] 17395
## [32,] 14662
## [33,] 16348
## [34,] 15269
## [35,] 14432
## [36,] 15345
## [37,] 18022
## [38,] 16761
## [39,] 19528
## [40,] 18625
## [41,] 9481
## [42,] 16104
## [43,] 14757
## [44,] 15731
## [45,] 13281
## [46,] 12572
## [47,] 16516
## [48,] 19860
## [49,] 20670
## [50,] 18516
## [51,] 20554
## [52,] 15781
## LC08_L2SP_008058_20210410_20210416_02_T1_SR_B7
## [1,] 13660
## [2,] 11645
## [3,] 12036
## [4,] 12478
## [5,] 11425
## [6,] 14912
## [7,] 16829
## [8,] 15185
## [9,] 17361
## [10,] 16867
## [11,] 15189
## [12,] 16677
## [13,] 20197
## [14,] 15521
## [15,] 16523
## [16,] 12351
## [17,] 16693
## [18,] 16706
## [19,] 12632
## [20,] 19769
## [21,] 12917
## [22,] 15225
## [23,] 14311
## [24,] 14966
## [25,] 11851
## [26,] 13068
## [27,] 12657
## [28,] 15653
## [29,] 14856
## [30,] 14727
## [31,] 13955
## [32,] 11772
## [33,] 13649
## [34,] 13454
## [35,] 12313
## [36,] 13108
## [37,] 15706
## [38,] 14404
## [39,] 16728
## [40,] 15240
## [41,] 8135
## [42,] 11873
## [43,] 10397
## [44,] 10996
## [45,] 10330
## [46,] 10348
## [47,] 13796
## [48,] 17502
## [49,] 18012
## [50,] 15662
## [51,] 18007
## [52,] 13135DF2 <- aggregate(DF1,list(PUNTOS_geo$CLASE), mean)
DF2## Group.1 LC08_L2SP_008058_20210410_20210416_02_T1_SR_B1
## 1 AGUA 11561.462
## 2 CASAS 10657.333
## 3 SUELOS 9812.417
## 4 VEGETACION 10261.750
## LC08_L2SP_008058_20210410_20210416_02_T1_SR_B2
## 1 11879.85
## 2 10959.53
## 3 10386.67
## 4 10483.92
## LC08_L2SP_008058_20210410_20210416_02_T1_SR_B3
## 1 13411.46
## 2 13093.87
## 3 12518.75
## 4 11950.67
## LC08_L2SP_008058_20210410_20210416_02_T1_SR_B4
## 1 13477.31
## 2 13292.53
## 3 12874.17
## 4 11646.17
## LC08_L2SP_008058_20210410_20210416_02_T1_SR_B5
## 1 19042.15
## 2 19265.40
## 3 17381.08
## 4 20282.42
## LC08_L2SP_008058_20210410_20210416_02_T1_SR_B6
## 1 17315.38
## 2 17706.87
## 3 16926.75
## 4 16151.92
## LC08_L2SP_008058_20210410_20210416_02_T1_SR_B7
## 1 14958.54
## 2 14722.87
## 3 14159.33
## 4 13182.75# EXPORTAR EL DF A UN EXCEL
#install.packages("writexl")
library(writexl)## Warning: package 'writexl' was built under R version 4.0.5write_xlsx(DF2,"C:/Users/Jairo Cuenca/Downloads/PRACTICA 05/INFORMACION P05.xlsx")GRAFICAS <- raster("C:/Users/Jairo Cuenca/Downloads/PRACTICA 05/CURVAS ESPECTRALES.png")
plot(GRAFICAS)