Skip to content

Commit 353caca

Browse files
authored
Update caso_estudio.md
Typos y alguna corrección a capitalización de títulos
1 parent 7d0bd7d commit 353caca

File tree

1 file changed

+30
-29
lines changed

1 file changed

+30
-29
lines changed

episodes/caso_estudio.md

Lines changed: 30 additions & 29 deletions
Original file line numberDiff line numberDiff line change
@@ -25,7 +25,7 @@ exercises: 0
2525

2626
# Deforestación en Maranhão
2727

28-
La Amazonía es una de las regiones más biodiversas del planeta y un componente clave del sistema climático global que además sostiene múltiples comunidades indígenas. En particular el estado de Maranhão, en Brasil, es uno de los focos más críticos de deforestación en el país. Se estima que el 76 % de la cobertura original de bosque amazónico en este estado ha sido destruida. Según Global Forest Watch, Maranhão ha registrado una de las tasas más altas de pérdida de cobertura boscosa en Brasil en los últimos años, impulsada por incendios, expansión agropecuaria y tala ilegal. Estos procesos están estrechamente ligados a la fragmentación ecológica, la pérdida de biodiversidad y la violencia hacia comunidades indígenas. Frente a este escenario, el monitoreo sistemático de los cambios en la cobertura vegetal es fundamental. Los productos OPERA DIST-HLS, derivados principalmente de Landsat (NASA/USGS) y Sentinel-2 (ESA), ofrecen una herramienta poderosa para detectar disturbios recientes y aportar evidencia clave para la conservación y la formulación de políticas públicas basadas en datos.
28+
La Amazonía es una de las regiones más biodiversas del planeta y un componente clave del sistema climático global, que además sostiene múltiples comunidades indígenas. En particular, el estado de Maranhão, en Brasil, es uno de los focos más críticos de deforestación en el país. Se estima que el 76 % de la cobertura original de bosque amazónico en este estado ha sido destruida. Según Global Forest Watch, Maranhão ha registrado una de las tasas más altas de pérdida de cobertura boscosa en Brasil en los últimos años, impulsada por incendios, expansión agropecuaria y tala ilegal. Estos procesos están estrechamente ligados a la fragmentación ecológica, la pérdida de biodiversidad y la violencia hacia comunidades indígenas. Frente a este escenario, el monitoreo sistemático de los cambios en la cobertura vegetal es fundamental. Los productos OPERA DIST-HLS, derivados principalmente de Landsat (NASA/USGS) y Sentinel-2 (ESA), ofrecen una herramienta poderosa para detectar disturbios recientes y aportar evidencia clave para la conservación y la formulación de políticas públicas basadas en datos.
2929

3030
<figure style="text-align: center;">
3131
<img src="https://tse3.mm.bing.net/th/id/OIP.sSfdF5nBFUWbh3UImdbyVgHaE7?pid=Api" alt="Buriticupu - erosión" style="max-width: 100%; height: auto;">
@@ -36,21 +36,21 @@ La Amazonía es una de las regiones más biodiversas del planeta y un componente
3636

3737
## Ruta de trabajo
3838

39-
Nuestro objetivo es evaluar la deforestación en un area cercana a la ciudad de Buriticupu en el estado Maranhao.
40-
Para eso en esta notebook vamos a :
39+
Nuestro objetivo es evaluar la deforestación en un área cercana a la ciudad de Buriticupu en el estado de Maranhão.
40+
Para eso en esta notebook vamos a:
4141

42-
1. Filtrar y seleccionar los productos OPERA DIST-ALERT desde la nube
43-
2. Visualizar y explorar los subproductos VEG_DIST_STATUS
42+
1. Filtrar y seleccionar los productos OPERA DIST-ALERT desde la nube.
43+
2. Visualizar y explorar los subproductos VEG_DIST_STATUS.
4444
3. Gráficar la evolución del disturbio a lo largo del tiempo.
45-
4. Generar un mapa de disturbios
46-
5. Explorar subproducto VEG_DIST_DATE
45+
4. Generar un mapa de disturbios.
46+
5. Explorar el subproducto VEG_DIST_DATE.
4747

4848
### Antes de empezar - Importar librerías que vamos a utilizar
4949

5050
````python
5151

5252
#librerias para manipulación de datos
53-
from warnings import filterwarnings #suprime los warning
53+
from warnings import filterwarnings #suprime los warnings
5454
filterwarnings('ignore')
5555
import numpy as np, pandas as pd, xarray as xr
5656
import rioxarray as rio
@@ -72,32 +72,32 @@ gdal.SetConfigOption('CPL_VSIL_CURL_ALLOWED_EXTENSIONS','TIF, TIFF')
7272

7373
````
7474

75-
### FILTRAR Y SELECCIONAR LOS PRODUCTOS OPERA DESDE LA NUBE
75+
### Filtrar y seleccionar los productos OPERA desde la nube
7676

7777

78-
#### **1.a. Seleccionar el area de estudio**
78+
#### **1.a. Seleccionar el área de estudio**
7979

8080
Usa la herramienta [Bounding Box](https://boundingbox.klokantech.com/) para obtener las coordenadas geográficas (latitud y longitud) del área seleccionada.
8181

8282
Bounding Box es un selector visual de cajas geográficas (bounding boxes) que permite:
8383

84-
- Dibujar un rectángulo o polígono directamente sobre un mapa interactivo para delimitar un área específica-
84+
- Dibujar un rectángulo o polígono directamente sobre un mapa interactivo para delimitar un área específica.
8585
- Obtener las coordenadas geográficas (latitud y longitud) del área seleccionada y copiarlas fácilmente en diversos formatos útiles como MARC, DublinCore, KML, GeoJSON, OGC WKT, CSV, FGDC, entre otros.
8686

87-
Sigue los siguientes pasos:
87+
Sigue estos pasos:
8888

8989
1. Navega hasta la página https://boundingbox.klokantech.com/
9090
2. Busca la zona de interés y dibuja un rectángulo sobre el mapa.
9191
3. En la sección "Copy & Paste", selecciona el formato "CSV".
92-
4. Copia las coordenadas
92+
4. Copia las coordenadas.
9393

9494
Estas coordenadas están en el orden correcto requerido por STAC:
9595

9696
`bbox = [xmin, ymin, xmax, ymax] = [long_oeste, lat_sur, long_este, lat_norte]`
9797

9898
El siguiente ejemplo supone que las coordenadas copiadas de boundingbox son **-46.52993,-4.383815,-43.363075,-4.243793**
9999

100-
5. Define el area utilizando las coordenadas copiadas y define el rango de fechas de interés.
100+
5. Define el área de interés (AOI, por Area of Interest en inglés) utilizando las coordenadas copiadas y define el rango de fechas de interés.
101101

102102
````python
103103

@@ -126,7 +126,7 @@ search_params = {
126126
client = Client.open("https://cmr.earthdata.nasa.gov/stac/LPCLOUD/")
127127
items = list(client.search(**search_params).get_items())
128128

129-
# Extraemos fechas disponibles
129+
# Extraemos las fechas disponibles
130130
fechas = sorted({item.datetime.date() for item in items})
131131
print(f"Fechas disponibles ({len(fechas)}):")
132132
print(fechas)
@@ -188,8 +188,9 @@ Output is truncated. View as a scrollable element or open in a text editor. Adju
188188

189189
````python
190190

191-
# Recorremos todos los productos encontrados y seleccionamos solo los archivos .tif
192-
# correspondientes al asset 'VEG-DIST-STATUS'
191+
# Recorremos todos los productos encontrados
192+
# y seleccionamos solo los archivos .tif
193+
# correspondientes al asset 'VEG-DIST-STATUS'.
193194
# Guardamos la fecha del producto y el link al archivo
194195

195196
veg_status_assets = []
@@ -299,14 +300,14 @@ Se encontraron 12 archivos con menos de 40% de nubes.
299300
...
300301

301302

302-
### VISUALIZAR Y EXPLORAR Productos VEG_DIST_STATUS
303+
### Visualizar y explorar productos VEG_DIST_STATUS
303304

304305

305306
````python
306307

307-
#Visualizar productos VEG_DIST_STATUS en el area de interes.
308+
#Visualizar productos VEG_DIST_STATUS en el área de interés.
308309

309-
#convertir a shp las coordenadas del area de interes (AOI)
310+
#convertir a shp las coordenadas del AOI
310311
from shapely.geometry import box
311312
import geopandas as gpd
312313

@@ -315,7 +316,7 @@ aoi_coords = [-46.78, -4.61, -46.58, -4.41] # xmin, ymin, xmax, ymax
315316
aoi_geom = box(*aoi_coords)
316317
AOI = gpd.GeoDataFrame(geometry=[aoi_geom], crs="EPSG:4326")
317318

318-
#Visualizar la primer y ultima fecha del los productos filtrados
319+
#Visualizar primera y ultima fechas de los productos filtrados
319320
import matplotlib.pyplot as plt
320321
import matplotlib.colors as mcolors
321322
import geopandas as gpd
@@ -364,7 +365,7 @@ Valores del producto `VEG-DIST-STATUS`:
364365
- **0:** Sin alteración
365366
- **1:** Primera detección de alteraciones con cambios en la cobertura vegetal <50%
366367
- **2:** Detección provisional de alteraciones con cambios en la cobertura vegetal <50%
367-
- **3:** Detección confirmada de alteraciones con cambios en la cobertura vegetal < 50%
368+
- **3:** Detección confirmada de alteraciones con cambios en la cobertura vegetal <50%
368369
- **4:** Primera detección de alteraciones con cambios en la cobertura vegetal ≥50%
369370
- **5:** Detección provisional de alteraciones con cambios en la cobertura vegetal ≥50%
370371
- **6:** Detección confirmada de alteraciones con cambios en la cobertura vegetal ≥50%
@@ -418,7 +419,7 @@ plt.show()
418419
![](fig/output2.png)
419420

420421

421-
### EVOLUCIÓN DEL DISTURBIO A LO LARGO DEL TIEMPO
422+
### Evolución del disturbio a lo largo del tiempo
422423

423424

424425
````python
@@ -510,8 +511,8 @@ plt.show()
510511

511512
**¿Hay algo raro en los gráficos?**
512513

513-
Tal vez notes que una fecha tiene dos puntos. Eso es porque el producto OPERA usa imágenes SENTINEL y LANDSAT, y puede ocurrir que en alguna fecha haya dos productos. En ese caso, es necesario seleccionar alguno de ellos en base a algún criterio.
514-
Debajo, seleccionamos la imagen con menos nubes y, de esa forma, cuando volvemos a ejecutar el código que genera el gráfico encontramos solo un producto por fecha.
514+
Tal vez notes que una fecha tiene dos puntos. Eso es porque el producto OPERA usa imágenes SENTINEL y LANDSAT, y puede ocurrir que en alguna fecha haya dos productos. En ese caso, es necesario seleccionar alguno de ellos sobre la base de algún criterio.
515+
Debajo, seleccionamos la imagen con menos nubes y, de esa forma, cuando volvemos a ejecutar el código que genera el gráfico, encontramos solo un producto por fecha.
515516

516517

517518
:::
@@ -532,7 +533,7 @@ df_filtrados = df_filtrados.drop_duplicates(subset="fecha", keep="first")
532533
# Reconstruir la lista filtrada
533534
filtrados_unicos = df_filtrados.to_dict(orient="records")
534535

535-
#VOLVEMOS A CORRER EL STACK USANDO filtrados_unicos
536+
# VOLVEMOS A CORRER EL STACK USANDO filtrados_unicos
536537

537538
# Stack de los subproductos VEG-DIST-STATUS
538539

@@ -559,7 +560,7 @@ stack["time"] = fechas
559560

560561
````
561562

562-
### GENERAR UN MAPA DE DISTUBIOS
563+
### Generar un mapa de disturbios
563564

564565

565566

@@ -603,7 +604,7 @@ hvplot_map
603604
![](fig/output4.png)
604605

605606

606-
### EXPLORAR SUBPRODUCTO VEG_DIST-DATE
607+
### Explorar subproducto VEG_DIST-DATE
607608

608609
````python
609610

@@ -620,7 +621,7 @@ search_params = {
620621
"collections": ["OPERA_L3_DIST-ALERT-HLS_V1_1"]
621622
}
622623

623-
# Buscar items en el catálogo
624+
# Buscar ítems en el catálogo
624625
items = list(catalog.search(**search_params).get_items())
625626

626627
# Filtrar solo los assets VEG-DIST-DATE accesibles por HTTPS

0 commit comments

Comments
 (0)