momento 2
ANDRES FELIPE GUZMAN NARVAEZ
25/4/2021
library(raster)## Warning: package 'raster' was built under R version 4.0.5## 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/Usuario/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/Usuario/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/Usuario/Documents/R/win-library/4.0/rgdal/proj#concatenamos los nombres de las bandas
bandas<- paste0("C:/Users/Usuario/Desktop/electiva II/PROYECTO2/IMAGEN/LC08_L1TP_008058_20180317_20200901_02_T1_B",1:7,".tif")
bandas## [1] "C:/Users/Usuario/Desktop/electiva II/PROYECTO2/IMAGEN/LC08_L1TP_008058_20180317_20200901_02_T1_B1.tif"
## [2] "C:/Users/Usuario/Desktop/electiva II/PROYECTO2/IMAGEN/LC08_L1TP_008058_20180317_20200901_02_T1_B2.tif"
## [3] "C:/Users/Usuario/Desktop/electiva II/PROYECTO2/IMAGEN/LC08_L1TP_008058_20180317_20200901_02_T1_B3.tif"
## [4] "C:/Users/Usuario/Desktop/electiva II/PROYECTO2/IMAGEN/LC08_L1TP_008058_20180317_20200901_02_T1_B4.tif"
## [5] "C:/Users/Usuario/Desktop/electiva II/PROYECTO2/IMAGEN/LC08_L1TP_008058_20180317_20200901_02_T1_B5.tif"
## [6] "C:/Users/Usuario/Desktop/electiva II/PROYECTO2/IMAGEN/LC08_L1TP_008058_20180317_20200901_02_T1_B6.tif"
## [7] "C:/Users/Usuario/Desktop/electiva II/PROYECTO2/IMAGEN/LC08_L1TP_008058_20180317_20200901_02_T1_B7.tif"landsat <- stack(bandas)
landsat## class : RasterStack
## dimensions : 7731, 7581, 58608711, 7 (nrow, ncol, ncell, nlayers)
## resolution : 30, 30 (x, y)
## extent : 413385, 640815, 203385, 435315 (xmin, xmax, ymin, ymax)
## crs : +proj=utm +zone=18 +datum=WGS84 +units=m +no_defs
## names : LC08_L1TP//1_02_T1_B1, LC08_L1TP//1_02_T1_B2, LC08_L1TP//1_02_T1_B3, LC08_L1TP//1_02_T1_B4, LC08_L1TP//1_02_T1_B5, LC08_L1TP//1_02_T1_B6, LC08_L1TP//1_02_T1_B7
## min values : 0, 0, 0, 0, 0, 0, 0
## max values : 65535, 65535, 65535, 65535, 65535, 65535, 65535#install.packages("sf")
library(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.1library(readxl)## Warning: package 'readxl' was built under R version 4.0.5# LLAMAMOS EL ARCHIVO EXCEL CON PUNTOS
puntos <- read_xlsx("C:/Users/Usuario/Desktop/electiva II/PROYECTO2/IMAGEN/puntos.xlsx")
#puntos
#CONCVERTIR EXCEL EN BDE
puntos_geo <- puntos %>%
st_as_sf(coords = c("x","y"))
plot(puntos_geo)
df1 <- extract(landsat,puntos_geo)
df1## LC08_L1TP_008058_20180317_20200901_02_T1_B1
## [1,] 11282
## [2,] 11327
## [3,] 10908
## [4,] 12428
## [5,] 11086
## [6,] 11324
## [7,] 11206
## [8,] 11093
## [9,] 10887
## [10,] 10908
## [11,] 10808
## [12,] 10883
## [13,] 10954
## [14,] 11185
## [15,] 11494
## [16,] 10947
## [17,] 11387
## [18,] 11051
## [19,] 12980
## [20,] 11150
## [21,] 11043
## [22,] 10971
## [23,] 11048
## [24,] 11453
## [25,] 11111
## [26,] 11040
## [27,] 11858
## [28,] 10856
## [29,] 10757
## [30,] 10797
## [31,] 10838
## [32,] 10906
## [33,] 10967
## [34,] 10702
## [35,] 11046
## [36,] 12480
## [37,] 13271
## [38,] 12891
## [39,] 12044
## [40,] 13450
## [41,] 12274
## [42,] 12431
## [43,] 12596
## [44,] 12470
## [45,] 12721
## [46,] 12010
## [47,] 12880
## LC08_L1TP_008058_20180317_20200901_02_T1_B2
## [1,] 10470
## [2,] 10672
## [3,] 9909
## [4,] 11475
## [5,] 10242
## [6,] 10455
## [7,] 10316
## [8,] 10205
## [9,] 9923
## [10,] 9889
## [11,] 9798
## [12,] 9880
## [13,] 10022
## [14,] 10327
## [15,] 10795
## [16,] 10101
## [17,] 10727
## [18,] 10149
## [19,] 12839
## [20,] 10261
## [21,] 10093
## [22,] 9952
## [23,] 10128
## [24,] 10588
## [25,] 10197
## [26,] 10069
## [27,] 11158
## [28,] 9863
## [29,] 9765
## [30,] 9816
## [31,] 9867
## [32,] 9925
## [33,] 10016
## [34,] 9739
## [35,] 10135
## [36,] 11837
## [37,] 12610
## [38,] 12520
## [39,] 11163
## [40,] 12920
## [41,] 11538
## [42,] 11842
## [43,] 11957
## [44,] 11895
## [45,] 12080
## [46,] 11255
## [47,] 12146
## LC08_L1TP_008058_20180317_20200901_02_T1_B3
## [1,] 9889
## [2,] 10243
## [3,] 9372
## [4,] 11211
## [5,] 9549
## [6,] 9792
## [7,] 9627
## [8,] 9438
## [9,] 9167
## [10,] 8928
## [11,] 8744
## [12,] 8968
## [13,] 9518
## [14,] 9850
## [15,] 10276
## [16,] 9717
## [17,] 10453
## [18,] 9800
## [19,] 14143
## [20,] 9997
## [21,] 9684
## [22,] 9325
## [23,] 9767
## [24,] 9776
## [25,] 9489
## [26,] 9492
## [27,] 10692
## [28,] 9357
## [29,] 9095
## [30,] 9249
## [31,] 9346
## [32,] 9444
## [33,] 9436
## [34,] 9192
## [35,] 9800
## [36,] 11489
## [37,] 12086
## [38,] 11889
## [39,] 10674
## [40,] 12285
## [41,] 11181
## [42,] 11519
## [43,] 11618
## [44,] 11776
## [45,] 11450
## [46,] 10684
## [47,] 11693
## LC08_L1TP_008058_20180317_20200901_02_T1_B4
## [1,] 9153
## [2,] 10253
## [3,] 8148
## [4,] 10342
## [5,] 8992
## [6,] 9321
## [7,] 9026
## [8,] 8811
## [9,] 8142
## [10,] 7992
## [11,] 7780
## [12,] 7995
## [13,] 8563
## [14,] 9151
## [15,] 10469
## [16,] 8842
## [17,] 10623
## [18,] 8790
## [19,] 16160
## [20,] 8903
## [21,] 8381
## [22,] 7974
## [23,] 8630
## [24,] 9365
## [25,] 8860
## [26,] 8408
## [27,] 10910
## [28,] 8092
## [29,] 7908
## [30,] 7970
## [31,] 8164
## [32,] 8150
## [33,] 8297
## [34,] 8006
## [35,] 8702
## [36,] 11363
## [37,] 11627
## [38,] 11610
## [39,] 9936
## [40,] 11525
## [41,] 10981
## [42,] 10998
## [43,] 11333
## [44,] 11608
## [45,] 11009
## [46,] 10283
## [47,] 11313
## LC08_L1TP_008058_20180317_20200901_02_T1_B5
## [1,] 8474
## [2,] 17093
## [3,] 21435
## [4,] 18856
## [5,] 8987
## [6,] 9489
## [7,] 9281
## [8,] 8989
## [9,] 9754
## [10,] 11895
## [11,] 9616
## [12,] 14464
## [13,] 12398
## [14,] 15907
## [15,] 16478
## [16,] 20141
## [17,] 18387
## [18,] 19661
## [19,] 22161
## [20,] 19423
## [21,] 16681
## [22,] 20737
## [23,] 20584
## [24,] 13591
## [25,] 11705
## [26,] 21176
## [27,] 14783
## [28,] 20070
## [29,] 21016
## [30,] 20098
## [31,] 21264
## [32,] 20571
## [33,] 19764
## [34,] 21173
## [35,] 21743
## [36,] 17551
## [37,] 14824
## [38,] 16697
## [39,] 18263
## [40,] 16167
## [41,] 17035
## [42,] 18858
## [43,] 17871
## [44,] 18968
## [45,] 16528
## [46,] 16286
## [47,] 16347
## LC08_L1TP_008058_20180317_20200901_02_T1_B6
## [1,] 6981
## [2,] 17329
## [3,] 12932
## [4,] 14954
## [5,] 7127
## [6,] 7805
## [7,] 7503
## [8,] 7182
## [9,] 7655
## [10,] 8136
## [11,] 7435
## [12,] 9416
## [13,] 8994
## [14,] 13131
## [15,] 20033
## [16,] 15119
## [17,] 16353
## [18,] 15220
## [19,] 22362
## [20,] 15563
## [21,] 11500
## [22,] 12393
## [23,] 14741
## [24,] 13820
## [25,] 10085
## [26,] 13513
## [27,] 14996
## [28,] 12467
## [29,] 11689
## [30,] 11551
## [31,] 13408
## [32,] 12746
## [33,] 12827
## [34,] 13249
## [35,] 14945
## [36,] 16044
## [37,] 15164
## [38,] 15081
## [39,] 14263
## [40,] 15612
## [41,] 14980
## [42,] 15990
## [43,] 16423
## [44,] 15831
## [45,] 15790
## [46,] 15124
## [47,] 15365
## LC08_L1TP_008058_20180317_20200901_02_T1_B7
## [1,] 6250
## [2,] 12662
## [3,] 8667
## [4,] 11842
## [5,] 6322
## [6,] 6777
## [7,] 6555
## [8,] 6316
## [9,] 6530
## [10,] 6706
## [11,] 6340
## [12,] 7215
## [13,] 7029
## [14,] 9319
## [15,] 15837
## [16,] 10105
## [17,] 11515
## [18,] 10062
## [19,] 18823
## [20,] 10289
## [21,] 8281
## [22,] 8413
## [23,] 9648
## [24,] 11721
## [25,] 7944
## [26,] 8840
## [27,] 11729
## [28,] 8246
## [29,] 7843
## [30,] 7815
## [31,] 8713
## [32,] 8365
## [33,] 8582
## [34,] 8663
## [35,] 9793
## [36,] 12952
## [37,] 13570
## [38,] 12769
## [39,] 11214
## [40,] 13202
## [41,] 11914
## [42,] 12670
## [43,] 13274
## [44,] 11562
## [45,] 13338
## [46,] 11994
## [47,] 12836df2 <- aggregate(df1,list(puntos_geo$clase), mean)
df2## Group.1 LC08_L1TP_008058_20180317_20200901_02_T1_B1
## 1 agua 11055.17
## 2 constucciones 12611.23
## 3 suelos 11350.09
## 4 vegetacion 10976.36
## LC08_L1TP_008058_20180317_20200901_02_T1_B2
## 1 10138.50
## 2 11941.38
## 3 10573.18
## 4 10029.45
## LC08_L1TP_008058_20180317_20200901_02_T1_B3
## 1 9439.167
## 2 11504.231
## 3 10289.182
## 4 9469.455
## LC08_L1TP_008058_20180317_20200901_02_T1_B4
## 1 8635.667
## 2 11071.385
## 3 9853.636
## 4 8446.636
## LC08_L1TP_008058_20180317_20200901_02_T1_B5
## 1 11749.75
## 2 17250.08
## 3 18630.64
## 4 19368.64
## LC08_L1TP_008058_20180317_20200901_02_T1_B6
## 1 8859.167
## 2 15432.385
## 3 15857.545
## 4 12678.455
## LC08_L1TP_008058_20180317_20200901_02_T1_B7
## 1 7096.000
## 2 12549.000
## 3 11577.818
## 4 8673.364# 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/Usuario/Desktop/electiva II/PROYECTO2/IMAGEN/dataframe.xlsx")grafica <- raster("C:/Users/Usuario/Desktop/electiva II/PROYECTO2/curva espectral.png")
plot(grafica)
#install.packages("Rcpp")
library(Rcpp)
library(leaflet)## Warning: package 'leaflet' was built under R version 4.0.5longitude_astudy <- c(-75.350937, -75.231709, -75.236316, -75.347638)
latitud_astudy <- c( 2.835484, 2.837110, 2.767743, 2.765785)
plot(longitude_astudy, latitud_astudy)
mapa_uno <- leaflet() %>%
addTiles() %>%
addMarkers(lng = -75.257864, lat = 2.779560, popup="casa")
addPolygons(mapa_uno,longitude_astudy,latitud_astudy)