Ahora, se hará lo mismo para conocer la producción de aguacate en Valle del Cauca:
azucar_valle <- eva_valle %>% filter(CULTIVO == "CANA AZUCARERA") %>% dplyr::select(MUNICIPIO, COD_MUN, YEAR, PERIODO, ton_Prod, RENDIM) azucar_valleunique(azucar_valle$YEAR) [1] 2007 2008 2009 2010 2011 2012 2013 2014 2015 2016 2017 2018azucar_valle %>% replace(is.na(.), 0) -> azucar_valle2azucar_valle %>% group_by(MUNICIPIO, COD_MUN, YEAR) %>%
summarize(ton_Prod=sum(ton_Prod)) -> azucar_valle2`summarise()` has grouped output by 'MUNICIPIO', 'COD_MUN'. You can override using the `.groups` argument.azucar_valle2 %>%
group_by(COD_MUN) %>%
gather("ton_Prod", key = variable, value = number) %>%
unite(combi, variable, YEAR) %>%
pivot_wider(names_from = combi, values_from = number, values_fill = 0) -> azucar_valle3mun_valle_azucar = left_join(mun_valle, azucar_valle3, by="COD_MUN")Se hace el mapa:
bins <- c(0, 5000, 20000, 50000, 100000, 250000, 500000, 1000000, 1500000)
pal <- colorBin("PRGn", domain = mun_valle_azucar$ton_Prod_2018, bins = bins)
mapa6 <- leaflet(data = mun_valle_azucar) %>%
addTiles() %>%
addPolygons(label = ~ton_Prod_2018,
popup = ~MPIO_CNMBR,
fillColor = ~pal(ton_Prod_2018),
color = "#444444",
weight = 1,
smoothFactor = 0.5,
opacity = 1.0,
fillOpacity = 0.5,
highlightOptions = highlightOptions(color = "white", weight = 2, bringToFront = TRUE)
) %>%
addProviderTiles(providers$OpenStreetMap) %>%
addLegend("bottomright", pal = pal, values = ~ton_Prod_2018,
title = "Producción de caña de azúcar en Valle del Cauca [Ton] (2018)",
opacity = 1
)Some values were outside the color scale and will be treated as NAmapa6Se hace un mapa para conocer el área cosechada de caña de Azúcar:
azucar_valle <- eva_valle %>% filter(CULTIVO == "CANA AZUCARERA") %>% dplyr::select(MUNICIPIO, COD_MUN, YEAR, PERIODO, Ha_cosecha) azucar_valle %>% group_by(MUNICIPIO, COD_MUN, YEAR) %>%
summarize(Ha_cosecha=sum(Ha_cosecha)) -> azucar_valle2`summarise()` has grouped output by 'MUNICIPIO', 'COD_MUN'. You can override using the `.groups` argument.azucar_valle2 %>%
group_by(COD_MUN) %>%
gather("Ha_cosecha", key = variable, value = number) %>%
unite(combi, variable, YEAR) %>%
pivot_wider(names_from = combi, values_from = number, values_fill = 0) -> azucar_valle3mun_azucar_valle = left_join(mun_valle, azucar_valle3, by="COD_MUN")mun_azucar_valleSimple feature collection with 42 features and 24 fields
geometry type: MULTIPOLYGON
dimension: XY
bbox: xmin: -77.54977 ymin: 3.091239 xmax: -75.70724 ymax: 5.047394
geographic CRS: WGS 84
First 10 features:
DPTO_CCDGO MPIO_CCDGO MPIO_CNMBR MPIO_CRSLC
1 76 76001 CALI 1536
2 76 76036 ANDALUCÍA Ordenanza 38 de Abril 25 de 1921
3 76 76041 ANSERMANUEVO Ordenanza 29 de 1925
4 76 76054 ARGELIA Ordenanza 15 de 1956
5 76 76111 BUGA 1555
6 76 76113 BUGALAGRANDE 1791
7 76 76122 CAICEDONIA Ordenanza 21 del 20 de Abril de 1923
8 76 76126 CALIMA Ordenanza 49 de Junio 23 de 1939
9 76 76130 CANDELARIA 1797
10 76 76147 CARTAGO 1863
MPIO_NAREA MPIO_NANO DPTO_CNMBR Shape_Leng Shape_Area KM2
1 563.04276 2017 VALLE DEL CAUCA 1.1636697 0.045733211 563.044
2 110.46042 2017 VALLE DEL CAUCA 0.6651004 0.008984851 110.460
3 305.45118 2017 VALLE DEL CAUCA 0.9460442 0.024871077 305.451
4 90.79604 2017 VALLE DEL CAUCA 0.4538261 0.007391017 90.796
5 825.86513 2017 VALLE DEL CAUCA 2.0063716 0.067157337 825.865
6 396.78132 2017 VALLE DEL CAUCA 1.0336063 0.032279294 396.781
7 166.98369 2017 VALLE DEL CAUCA 0.6465900 0.013590500 166.984
8 793.49323 2017 VALLE DEL CAUCA 1.5441977 0.064484981 793.497
9 296.46056 2017 VALLE DEL CAUCA 0.8709866 0.024086066 296.456
10 248.16005 2017 VALLE DEL CAUCA 0.8955580 0.020211481 248.224
COD_MUN MUNICIPIO Ha_cosecha_2007 Ha_cosecha_2008 Ha_cosecha_2009
1 76001 CALI 3438 2857 3437
2 76036 ANDALUCIA 163 1568 2223
3 76041 ANSERMANUEVO 1631 1286 1597
4 76054 <NA> NA NA NA
5 76111 BUGA 5352 5029 585
6 76113 BUGALAGRANDE 7259 6099 8069
7 76122 CAICEDONIA 0 0 106
8 76126 <NA> NA NA NA
9 76130 <NA> NA NA NA
10 76147 CARTAGO 2693 2233 3069
Ha_cosecha_2010 Ha_cosecha_2011 Ha_cosecha_2012 Ha_cosecha_2013
1 3051 3775 3886 3628
2 1805 1832 2151 2008
3 1391 1553 1823 1703
4 NA NA NA NA
5 5331 5351 6281 5865
6 7093 6117 718 6704
7 92 65 77 72
8 NA NA NA NA
9 NA NA NA NA
10 2825 3069 3603 3364
Ha_cosecha_2014 Ha_cosecha_2015 Ha_cosecha_2016 Ha_cosecha_2017
1 3932 4004 332 3492
2 2349 2219 1849 22
3 182 1409 148 1375
4 NA NA NA NA
5 5704 5698 5063 5046
6 7118 6604 5712 6405
7 82 94 55 4
8 NA NA NA NA
9 NA NA NA NA
10 3673 3703 3774 2946
Ha_cosecha_2018 geometry
1 3379 MULTIPOLYGON (((-76.59175 3...
2 1889 MULTIPOLYGON (((-76.22406 4...
3 1368 MULTIPOLYGON (((-76.01558 4...
4 NA MULTIPOLYGON (((-76.14316 4...
5 4703 MULTIPOLYGON (((-76.31608 3...
6 666 MULTIPOLYGON (((-76.15131 4...
7 88 MULTIPOLYGON (((-75.8539 4....
8 NA MULTIPOLYGON (((-76.51747 4...
9 NA MULTIPOLYGON (((-76.30455 3...
10 3287 MULTIPOLYGON (((-75.94518 4...bins <- c(0, 100, 500, 1000, 1500, 2000, 2500, 3000, 5000, 8000, 10000, 13000)
pal <- colorBin("Paired", domain = mun_azucar_valle$Ha_cosecha_2018, bins = bins)
mapa7 <- leaflet(data = mun_azucar_valle) %>%
addTiles() %>%
addPolygons(label = ~Ha_cosecha_2018,
popup = ~MPIO_CNMBR,
fillColor = ~pal(Ha_cosecha_2018),
color = "#444444",
weight = 1,
smoothFactor = 0.5,
opacity = 1.0,
fillOpacity = 0.5,
highlightOptions = highlightOptions(color = "white", weight = 2, bringToFront = TRUE)
) %>%
addProviderTiles(providers$OpenStreetMap) %>%
addLegend("bottomright", pal = pal, values = ~Ha_cosecha_2018,
title = "Área cosechada de caña de azúcar en Valle del Cauca [Ha] (2018)",
opacity = 1
)mapa7