UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS

(Universidad del Perú, Decana de América)

FACULTAD DE CIENCIAS SOCIALES
ESCUELA PROFESIONAL DE GEOGRAFÍA

TALLER DE GEOMATICA
Base 19
TRABAJO: “Aplicación de herramientas de la geomática para el análisis
cualitativo y cuantitativo de los Cambios de la Cobertura y Uso de la Tierra
en tres periodos de análisis”.
DOCENTE:
Varillas Arquiñigo María Luisa
INTEGRANTES:
Flores Castillo Ayrton (19150328)
Inga Basurto Eduardo Alonso (19150325)
CORREO:
[email protected]
[email protected]

Lima – Perú
AGOSTO 2022

1
ÍNDICE
1. INTRODUCCIÓN..............................................................................................................3
2. OBJETIVOS...........................................................................................................................3
2.1 Objetivo General................................................................................................................3
2.2 Objetivos específicos..........................................................................................................3
3. UBICACIÓN Y LOCALIZACIÓN DEL ÁMBITO DE AREA DE ESTUDIO................4
4. MATERIALES Y MÉTODOS..............................................................................................6
5. PROCESO METODOLÓGICO PAUTADO SOBRE EN CAMPO Y GABINETE
SOBRE EL USO DE LAS HERRAMIENTAS DE LA GEOMÁTICA.................................6
Proceso Ráster..........................................................................................................................7
Selección de imágenes de satélite. Años y meses.....................................................................7
Proceso de corrección de imágenes de satélite.........................................................................7
Combinación de bandas. Selección de combinación................................................................8
Uso de filtros y mejoras de las imágenes..................................................................................9
Herramienta Pan-Sharpening...............................................................................................9
Paso de Ráster a Vectorial......................................................................................................10
Proceso Vectorial...................................................................................................................11
Elaboración de unidades de la leyenda en base a metodología Corine Land Cover hasta el
nivel de detalle 2. Utilizar la misma leyenda para los tres mapas...........................................12
Elaboración de la base de datos de los CCUT........................................................................12
Elaboración de los mapas.......................................................................................................12
6. Resultados y Análisis............................................................................................................15
Cuantificación en hectáreas los cambios de la cobertura y uso de la tierra de un periodo a otro
por unidades de la leyenda Corine Land Cover......................................................................15
Presentación de gráficos estadísticos comparativos de los resultados....................................16
7. CONCLUSIONES DE LO TRABAJADO EN CAMPO Y GABINETE..........................17
8. RECOMENDACIONES......................................................................................................17
9. BIBLIOGRAFÍA..................................................................................................................18
Bibliografía.................................................................................................................................18
ANEXOS...................................................................................................................................18

2
 1. INTRODUCCIÓN
Los mapas han servido desde épocas inmemorables para ubicarnos espacialmente,
tienen un papel importante en la historia, pues sin ellos muchos territorios nunca
hubieran sido descubiertos, actualmente con el avance de la tecnología, el mundo se ha
digitalizado y los mapas se encuentran en todos lado (GPS, Google Maps, entre otros),
además contamos con muchos softwares centrados en los sistemas de información
geográfica, gracias a los avances de la ciencia desde el siglo XX, tenemos satélites que
nos brindan imágenes de nuestro planeta de una gran calidad, es así que contamos con
poderosas herramientas que hace algunas décadas se veían muy lejanas. En el presente
curso de Taller de Geomática, se nos ha encargado usar tales herramientas para hallar
los cambios de uso y cobertura de la tierra en tres periodos de tiempo, es así que se
demostrar el poder y la confiabilidad de estas herramientas, además de que se podrá
conocer el impacto que ha tenido el hombre en el aumento o disminución de cada
cobertura y uso de la tierra por año. A continuación, se pasará a desarrollar los pasos
para obtener tanto los mapas como los resultados de la variación de la cobertura y uso
de la tierra por años.

2. OBJETIVOS
2.1 Objetivo General.

Emplear las herramientas de geomática para detectar los cambios en la cobertura y uso
de la tierra en 3 periodos de tiempo.

2.2 Objetivos específicos.

● Conocer las características principales de la Subcuenca del Río Pariac

● Aprender a utilizar las herramientas del Software ArcGis
● Elaborar mapas de cobertura y uso de la tierra año por año
● Conocer la variación de la cobertura vegetal de un año a otro
● Entender la utilidad de las herramientas de los sistemas de información
geográfica

3
3. UBICACIÓN Y LOCALIZACIÓN DEL ÁMBITO DE AREA DE
ESTUDIO
Siendo la ubicación una de las características más importantes para la realización del
presente trabajo, pues conocer en qué parte de nuestro territorio se encuentra esta laguna
nos ayuda a situarnos espacialmente.

SUBCUENCA DEL RÍO PARIAC

Geográficamente se encuentra entre las coordenadas UTM (WGS84 Zona L-18 Sur):
222951.5595 E – 8939594.865 N (Desembocadura del río Pariac_P1) y 246344.358 E –
8947502.7151 N (Punto extremo al Este de la subcuenca - Nevado Huantsán_P2).

La quebrada Pariac – laguna Rajucolta, se ubica en, la provincia de Huaraz,

departamento de Ancash – Perú.

Se puede decir que la altitud de la subcuenca del río Pariac varía entre 6331 m.s.n.m. y
3125 m.s.n.m. y presenta un área total de 107,35661 km2. Por otra parte, tenemos al
clima, el cual es variable resaltando la precipitación entre 700 mm en el extremo
occidental, hasta 1000 mm hacia la parte oriental sobre los 5000 m.s.n.m., las
temperaturas son variables a lo largo de la subcuenca, alcanzando hasta 20°C de
diferencia entre el día y la noche.

4
 FIG.01 UBICACIÓN POLÍTICA DE LA SUBCUENCA DEL RÍO PARIAC

Ubicación Hidrográfica

Con respecto a este tema el INAIGEM nos menciona que:

Hidrográficamente, la subcuenca se localiza en la Cuenca del río Santa,

perteneciente a la Vertiente hidrográfica del Pacífico, que drena sus aguas
por la margen derecha del río Santa. Se encuentra ubicada al Sur de la
ciudad de Huaraz. Tiene un área de 107,35661 km2 y 63.720546 km de
perímetro. (INAIGEM, 2018, pág.11)

LAGUNA RAJUCOLTA

Con respecto a la Ubicación y la Localización de la laguna, el INAIGEM nos menciona

en otro informe que: “La laguna Rajucolta se ubica al pie del Nevado Huantsán (6366
m.s.n.m.). Hidrológicamente se encuentra en la vertiente del Pacifico, cuenca del río
Santa y subcuenca del río Pariac. Políticamente pertenece al distrito y provincia de
Huaraz, departamento de Ancash” (INAIGEM, 2016, pág.8).

Las coordenadas UTM WGS 84 son:

Norte 8946459

Este 242820 Altitud

4280 m.s.n.m.

5
 FIG.02 UBICACIÓN GEOGRÁFICA DE LA SUBCUENCA DEL RÍO PARIAC Y DE LA LAGUNA RAJUCOLTA.

4. MATERIALES Y MÉTODOS
Con respecto a los materiales utilizados en el presente trabajo, estos fueron mayormente
digitales, entre ellos tenemos a la página web EarthExplorer que nos proporcionó las
imágenes Satelitales, además se revisó distintos shapes de cobertura y uso de la tierra
para identificar estas unidades en nuestro mapa, sin embargo, optamos por realizar la
clasificación no supervisada como ayuda fundamental para determinar las áreas de
cambio.

 Software ArcGIS Pro.

 Página Web USGS herramienta EarthExplorer. (Imágenes Satelitales 2013,
2017 y 2021 junio)
 Ayuda de Shapes de Cobertura y uso de la Tierra (Cobertura Vegetal, Uso de
la Tierra de Ancash, Shape Geomorfológico de Ancash).
 Ayuda de Mapas de Clasificación Supervisada
 Utilización de la Clasificación No Supervisada
 Metodología Corine Land Cover

Con respecto a los métodos, estos se describirán más adelante, sin embargo, se puede
decir que con el objetivo de identificar zonas de variación o cambio en las áreas de
cobertura y uso de la tierra. se ha utilizado mayormente el Software ArcGIS Pro,
además procesos con el pan-sharpening que permitan aumentar la resolución de las
imágenes, luego se prosiguió con la determinación de las unidades de cobertura y uso de
la tierra con la ayuda del procedimiento de Clasificación No Supervisada ,
seguidamente se dio la reclasificación de los polígonos en 7 categorías poniéndoles
nombres y colores en base a la metodología Corine Land Cover, finalmente se calculó el
área de cada clase en cada periodo de tiempo (Junio de 2013, 2017 y 2021), en donde se
pudo evidenciar una variación en algunos casos. esto se explicará en los resultados

5. PROCESO METODOLÓGICO PAUTADO SOBRE EN CAMPO Y

GABINETE SOBRE EL USO DE LAS HERRAMIENTAS DE LA
GEOMÁTICA.

6
Proceso Ráster

Con respecto a los ráster, estos constan de una matriz de celdas (o pixeles) organizadas
en filas y columnas (o una cuadrícula) en la que cada celda contiene un valor que
representa información, como la temperatura. Los ráster pueden ser fotografías aéreas
digitales, imágenes satelitales, imágenes digitales o incluso mapas escaneados

Selección de imágenes de satélite. Años y meses.

El proceso de selección de las imágenes se realizó en la página Web USGS con la

herramienta EarthExplorer, este sitio web nos brinda de manera gratuita gran cantidad
de imágenes satelitales, de todas las opciones que había para escoger, se eligió las
imágenes del Satelite Landsat 8, LandSat Collection 1 – Level 1, luego Landsat 8
OLI/TIRS C1 Level – 1, luego se selección la ubicación de la laguna Rajucolta.

FIG.03
PROCESO DE BÚSQUEDA DE LAS IMÁGENES SATELITALES.

Con respecto a la fecha de las imágenes, estas fueron tomadas del mes de junio para los
años 2021, 2017 y 2013.

7
Proceso de corrección de imágenes de satélite

Con respecto a este punto, consideramos que no era necesario, debido a que no había
mayor inconveniente en la visualización de las imágenes, aunque sí sabemos de la
importancia de la corrección atmosférica, radiométrica, entre otras.

Combinación de bandas. Selección de combinación.

Esta temática es muy interesante con respecto a este el sitio Web MappingGis nos dice
que:

Una de las potencialidades de la Teledetección aplicadas al estudio del

medio que nos rodea, es la capacidad de discriminar diferentes cubiertas
vegetales, usos de suelo, masas de agua, o la detección de fenómenos
naturales o provocados por la actividad humana. Esto puede analizarse
gracias a la existencia de las diferentes bandas multiespectrales con las que
cuentan los satélites de hoy día. (MappingGis, 2019, párr.1).

Además, los datos captados por los satélites que se encargan de la Teledetección se
llegan a registrar en diferentes bandas del espectro electromagnético. Lo que a su vez
genera una imagen monocroma que podemos visualizar en escala de grises con una
paleta de 256 tonos. Por lo que cada pixel de la imagen puede contener un valor que
oscila entre el negro (valor 0) y el blanco (valor 256).

Podremos distribuir las bandas de las imágenes de satélite a través de tres canales: rojo,
verde y azul. El paso de cada banda por estos canales dotará de diferentes tonalidades
los elementos dependiendo de la respuesta de los mismos frente a las longitudes de
onda. (MappingGis, 2019, párr.3).

Para el caso de nuestras imágenes Landsat, se eligió la combinación que suele llamarse
“color natural”, esta involucra a las tres bandas visibles y se le asigna a cada una de
ellas su verdadero color, resultando una combinación que se aproxima a los colores
naturales de la escena. Podemos generarla utilizando la combinación de bandas 4,3,2
en Landsat 8 y 4,3,2 en Sentinel 2.

8
 FIG.04 IMAGEN DE LA ZONA ELEGIDA (RAJUCOLTA) EN COLOR NATURAL

Uso de filtros y mejoras de las imágenes.

El problema con las imágenes Landsat 8, es que estas tienen solo una resolución o
tamaño de celdas de 30, 30, por lo cual es muy difícil diferenciar los cambios de la
cobertura y uso de la tierra.

Herramienta Pan-Sharpening

Para solucionar esto, se utilizó la herramienta Create Pan – Sharpened Ráster Dataset en
el ArcToolBox, Geoprocesing, este método consiste en combinar por asi decirlo la
imagen multiespectral de baja resolución (30,30) con la banda 8 pancromática de una
mayor resolución, (15,15), es de esta manera que se realizó el proceso y se obtuvo una
imagen con un leve aumento de resolución (15,15).

9
FIG.05 DETALLES DE LA IMAGEN LANDSAT CON LA RESOLUCIÓN MEJORADA (15,15) GRACIAS AL
PAN – SHARPENING.

Paso de Ráster a Vectorial.

En este punto se decidió realizar un mapa de clasificación no supervisada de las

imágenes satelitales que obtuvimos, con la herramienta, ArctoolBox, luego Spatial
Analysis Tool, seleccionando la Opción Iso Cluster Unsupervised Classification, se
eleigió 12 clasificaciones debido a que no siempre los resultados serán los deseados
debido a ser una clasificación no supervisada.

FIG.06 12 CLASES AUTOMÁTICAS CON LA HERRAMIENTA ISO CLUSTER UNSUPERVISED

CLASSIFICATION

Una vez obtenidos los resultados se obtuvieron 12 clases, es así que se usó la
herramienta Reclassify para unir clases que presentan similitudes, en donde se clasificó
obteniéndose 7 clases.

Se transformó la clasificación a vector con la herramienta Conversion Tools, Raster to

Polygon con la opción simplificar polígonos activada.

10
 FIG.07 CONVERSIÓN DE RÁSTER A POLÍGONO

Proceso Vectorial.

Después de crear este polígono o vector se decidió reducir el número de polígonos

creados en el campo grid code, se creó un nuevo campo llamado cobertura y uso corine
land cover.

Para reducir el número de polígonos fue útil la herramienta Dissolve hallada en

Geoprocessing, seleccionando el campo Grid Code y el campo creado cobertura y uso
corine land cover, finalmente se insertó la ruta de salida.

Resultando en un nuevo vector con menos polígonos o 7 clases

FIG.08 REDUCCIÓN DE NÚMERO DE POLÍGONOS A SOLO 7

11
Elaboración de unidades de la leyenda en base a metodología Corine Land Cover hasta
el nivel de detalle 2. Utilizar la misma leyenda para los tres mapas.

Primero se nombró a nuestras 7 clases obtenidas de la siguiente manera: Clase 1:

Bosque Denso, Clase 2: Áreas Agrícolas, Clase 3: Zonas Húmedas, Clase 4: Zonas
Abiertas sin o con poca Vegetación (Glaciar), Clase 5: Zonas Abiertas sin o con poca
vegetación y Clase 6: Áreas Urbanas y 7 Aguas Continentales. El siguiente paso fue
elegir el color de cada uno de estos polígonos, para lo cual se descargó un paquete de
colores adaptado a la metodología corine land cover, de tal manera que se obtuvo el
siguiente resultado como leyenda.

FIG.09 LEYENDA CORINE LAND COVER ELEGIDA

Elaboración de la base de datos de los CCUT

La base de datos de cambios de cobertura y uso de tierra se realizará en un Excel donde

se anotarán los principales cambios de acuerdo a la clasificación realizada para cada
periodo de tiempo (2013, 2017, 2021).

Elaboración de los mapas

Antes de la realización de los mapas finales, se decidió realizar mapas intermedios de

áreas cambio y no cambio de los periodos (2013 - 2017) y (2017 - 2021).

Esto se realizó con la herramienta geoprocessing, intersect, insertando la capa de

cobertura 2013 y la de 2017, dándonos como resultado un nuevo mapa vectorial de
intersección, en la tabla de atributos de esta se calculó el área de los polígonos de
intersección en hectáreas, luego o bien con un código o manualmente se calculó cuáles
elementos de cobertura y uso de la tierra cambiaron de un periodo a otro, entre la

12
columna de 2017 y la de 2018. En una nueva Columna llamada cambio se evidenció los
cambios entre un periodo a otro asignándoles o bien un SI o un NO, y de esta columna
de SI o NO se creó el mapa de cambios para el periodo 2013 - 2017 y 2017 y 2021. Las
areas con cambio estan en color Rojo y las de no cambio de color Amarillo.

Periodo 2013 a 2017 presenta un área de NO CAMBIO de 4018.342 Hectáreas y de SI

CAMBIO de 1371.315 Hectáreas

13
Periodo 2017 a 2021 presenta un área de NO CAMBIO de 3993.331 Hectáreas y de SI
CAMBIO de 1396.818 Hectáreas

FIG.10 MAPAS DE CAMBIO Y NO CAMBIO (2013 - 2017) Y (2017 - 2021)

La Elaboración de los mapas se realizó con los datos de las áreas o polígonos de la
clasificación no supervisada hecha anteriormente de tal manera que se insertó este mapa
vectorial con los nombres y colores de acuerdo a la metodología corine land cover.

Luego se realizó la generación de los títulos, la generación del membrete, la escala, el

norte geográfico, la leyenda y un mapa de ubicación, Quedando como resultado final el
siguiente mapa que representa la cobertura y uso de la Tierra para la Subcuenca del Río
Pariac o la Laguna Rajucolta.

FIG.11 MAPA FINAL DE COBERTURA Y USO DE LA TIERRA (2017).

6.

Resultados y Análisis.
Con respecto a los resultados, de nuestros 7 clases: Clase 1: Bosque Denso, Clase 2:
Áreas Agrícolas, Clase 3: Zonas Húmedas, Clase 4: Zonas Abiertas sin o con poca
Vegetación (Glaciar), Clase 5: Zonas Abiertas sin o con poca vegetación, Clase 6: Áreas
Urbanas y 7 Aguas Continentales, la que más varió en decremento fue la categoría 4
Zonas Abiertas sin o con poca Vegetación (Glaciar), siendo que en cada año disminuye

14
su área en hectáreas, por otro lado también tenemos a la categoría 7 o aguas
continentales que también fue en disminución con el paso de los años. con respecto al
Bosque Denso, Áreas Urbanas, Zonas Abiertas sin o con poca Vegetación, estas
categorías variaron muy poco o casi nada, con respecto a las categorías de zonas
húmedas y áreas agrícolas hubo también muy poca variación, pero no se evidenció un
patrón de disminución año por año

Cuantificación en hectáreas los cambios de la cobertura y uso de la tierra de un periodo

a otro por unidades de la leyenda Corine Land Cover.

FIG.12 CUADRO DE UNIDADES O 7 CLASES EN HECTÁREAS Y CUADRO DE VARIACIÓN PORCENTUAL.

Presentación de gráficos estadísticos comparativos de los resultados.

15
 FIG.13 GRÁFICO COMPARATIVO DE LAS HECTÁREAS DE CADA CLASE ENTRE LOS 3 AÑOS (2013,
2017 Y 2021) Y GRÁFICO DE DISPERSIÓN.

7. CONCLUSIONES DE LO TRABAJADO EN CAMPO Y

GABINETE
Luego de concluir nuestros resultados a través de la metodología respectiva podemos
mencionar las siguientes conclusiones con respecto al trabajo de campo y  en gabinete:

- Con respecto a lo visualizado en el campo es notable a simple vista la disminución de 

las aguas continentales (laguna del Rajucolta),

- Los resultados de los mapas temáticos de cobertura y uso actual de suelos con la
metodología CLC sirven para la generación de información que apoyará a la toma de
decisiones orientadas a la planificación y ordenamiento del territorio.

- Ya en gabinete y con la ayuda del software GIS como sus herramientas y de las
imágenes satelitales podemos concluir que se denota la disminución de las Áreas
abiertas sin o con poca vegetación (Glaciares) y de las aguas continentales.  Con
respecto a lo último mencionado podemos agregar que el aumento de la temperatura
terrestre (calentamiento global) es el principal causante de la disminución de los
glaciares.

-Se observa un pequeño aumento en las Áreas agrícolas de nuestra tabla.

- La baja pendiente que existe en el Nevado del Huantsan ayuda de cierta manera a la
desaceleración del derretimiento del nevado en su zona de acumulación. Todo lo
contrario pasa en los flancos izquierdo y derecho de nuestra imagen Landsat en donde
también se observa masa glaciar, pero en mucho menor proporción y con una pendiente
más pronunciada.

16
8. RECOMENDACIONES

 Elegir una imagen satelital de acorde a un mes específico para todas las
imágenes (ejemplo: junio).
 Utilizar adecuadamente los nombres y los colores de la metodología Corine
Land Cover
 Calcular el área de cada cobertura de manera adecuada
 Evitar realizar actividades potencialmente dañinas para el medio ambiente

9. BIBLIOGRAFÍA

Bibliografía
INAIGEM. (Abril de 2016). RECONOCIMIENTO DEL NIVEL DE PELIGRO DE LA LAGUNA
RAJUCOLTA CON FINES DE SEGURIDAD Y APROVECHAMIENTO. Obtenido de
repositorio.inaigem.gob.pe:
https://repositorio.inaigem.gob.pe/server/api/core/bitstreams/cc700e7b-deb2-480d-
a136-4a3c2922cbba/content

INAIGEM. (Abril de 2018). INVESTIGACIÓN PARA LA ELABORACIÓN DEL MAPA DE PELIGROS, SU

REDUCCIÓN DE RIESGOS Y LA ADAPTACIÓN AL CAMBIO CLIMÁTICO EN AMBIENTES DE
SUBCUENCAS GLACIARES EN EL PERÚ EN LA SUBCUENCA DEL RÍO PARIAC, HUARAZ-
ANCASH. Obtenido de CDN:
https://cdn.www.gob.pe/uploads/document/file/914465/Informe_T
%C3%A9cnico_N__06_Subcuenca_Pariac.pdf

MappingGis. (Mayo de 2019). Combinación de bandas en imágenes de satélite Landsat y

Sentinel. Obtenido de mappinggis.com:
https://mappinggis.com/2019/05/combinaciones-de-bandas-en-imagenes-de-satelite-
landsat-y-sentinel/

17
ANEXOS.

18
Fig.14 Mapas intermedios ráster y vectorial: Mapas áreas de cambio y no cambio,
Mapas finales vectoriales: Mapas de Cobertura y Uso de la Tierra.

19
20
Fig.15 (MAPAS (GEOMOROFOLGICO, GEOLOGICO, DE PENDIENTES ,
CLASIFICACIÓN SUPERVISADA).

21