-
Notifications
You must be signed in to change notification settings - Fork 1
/
index.Rmd
75 lines (48 loc) · 3.69 KB
/
index.Rmd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
---
title: "Clasificación supervisada de cultivos- Campaña 2019-2020"
author: "Nazarena Rojas y Pamela Pairo"
output: html_document
---
```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE, message = F, warning = F)
```
```{r, echo=FALSE}
library(sf)
library(tidyverse)
library(mapview)
library(leafsync)
```
## Metodología general para las clasificaciones supervisadas
Se utilizó información elaborada por el INTA de las campañas Invierno de 2019 y Verano 2019/2020. La información está disponible en formato tiff [aquí](http://geointa.inta.gob.ar/publico/MNC/).
El informe generado de dicha campaña 2019/2020 que describe la metodología y los resultados obtenidos se puede consultar [aquí](http://www.geointa.inta.gob.ar/2020/10/06/mapa-nacional-de-cultivos-campana-2019-2020/)
Para la generación de los productos se utilizó la colección de imágenes Sentinel-2 disponibles en la plataforma de Google Earth Engine para el Departamento de Pergamino de la Provincia de Buenos Aires. En la clasificación de las mismas se utilizó el algoritmo Random Forest, 500 puntos por clase^[Para cultivos de invierno: suelo, arbustos, trigo y agua. Para cultivos de verano: suelo, soja, trigo, agua y arbustos] y se generaron los vectores con sus correspondientes áreas de cultivo. De todos los datos, se usaron el 70% para entrenar el modelo y el 30% para la validación.
Se calcularon las áreas sembradas para los cultivos de invierno (trigo, cebada y centeno) y verano (maíz y soja).
Además se calcularon las áreas sembradas según una clasificación generada a partir de tomar 500 puntos por clase tomando como verdades de campo lo generado por INTA.
## Cultivo de invierno 2019
A continuación se muestran los cultivos de trigo para la el 29-30 de septiembre por parte del INTA (violeta) y la clasificación realizada con imágenes Sentinel 2 (rojo)
```{r, echo=FALSE}
#Trigo_INTA <- read_sf("Pergamino/maiz_2020.shp")
Trigo_sep <- read_sf("Pergamino/interseccion_29-30_sep_2019.shp")
Trigo_INTA <- read_sf("Pergamino/Pergamino_Trigo.shp")
mapview(Trigo_sep, col.regions="red", color="red", label=FALSE)+ mapview(Trigo_INTA, col.regions="purple",color="purple", label=FALSE)
```
## Cultivos de verano 2019-2020
El gráfico que aparece a la izquierda corresponde a las clasificaciones realizadas para el cultivo de soja y a la derecha las clasificaciones para el cultivo de maiz. La clasificación realizada por el INTA con imagenes LANDSAT para la campaña de verano 2019-2020 se representa en violeta. La clasificación supervisada realizada con imágenes Sentinel-2 del 23-28 de febrero del 2020 está representada en rojo.
```{r, echo=FALSE}
Maiz_INTA <- read_sf("Pergamino/maiz_2020.shp")
Soja_INTA <- read_sf("Pergamino/soja_2020_f.shp")
Soja_Feb <- read_sf("Pergamino/soja_punto_febrero.shp")
Maiz_Feb <- read_sf("Pergamino/maiz_punto_febrero.shp")
p1 <- mapview(Soja_Feb, col.regions="red", color="red", labael= FALSE)+
mapview(Soja_INTA, col.regions="purple", color="purple", label=FALSE)
p2<-mapview(Maiz_Feb, col.regions="red", color="red", label=FALSE)+
mapview(Maiz_INTA, col.regions="purple", color="purple", label=FALSE)
sync(p1, p2)
```
Con el fin de identificar los cultivos de soja de segunda, se superpusieron las capas de TRIGO y de SOJA. La intersección entre ambas capas corresponde a la SOJA DE SEGUNDA y el área no superpuesta pertenecería a la SOJA DE PRIMERA.
```{r, echo=FALSE}
#Trigo_INTA <- read_sf("Pergamino/maiz_2020.shp")
Trigo_sep <- read_sf("Pergamino/interseccion_29-30_sep_2019.shp")
Soja_Feb <- read_sf("Pergamino/soja_punto_febrero.shp")
mapview(Trigo_sep, col.regions="orange", color="orange", label=FALSE)+ mapview(Soja_Feb, col.regions="purple",color="purple", label=FALSE)
```