0 votos
581 vistas
por Novato (106 puntos)   en ArcGIS Desktop

 

Buen día:

Tengo dos imágenes de cierta zona para dos periodos de tiempo distintos. Cada imagen tiene cierto número de coberturas, por ejemplo bosques, cultivos y zonas descubiertas. Lo que deseo es determinar el cambio de coberturas entre los dos periodos de tiempo, especificando que tipo de cobertura aumento, disminuyo, o se transformó  (por ejemplo el número de hectáreas de zona desnuda del periodo de tiempo 1 que paso a ser Bosques en el periodo de tiempo 2), para ello he realizado los siguientes procesos:

1. Corte las imágenes con un shapefile que contiene la zona de estudio, posteriormente realice una resta  (a través de la calculadora raster) entre las dos imágenes resultantes, con lo cual buscaba identificar el número de pixeles que cambiaron o se distribuyeron entre cada cobertura, lo que ocurre es que el resultado que me arroja el proceso me parece un poco incoherente, debido que  al sumar el número de pixeles de una cobertura en particular esperaría que me diera el número total de pixeles de esa categoría para cualquier periodo de tiempo analizado, lo cual no ocurre y no sé si es algún error en el proceso que realizo o si estoy interpretando mal los resultados.

2. Convertí las dos imágenes de las coberturas a shapefiles (CobA y CobB) y cree los siguientes campos para cada shapefile correspondiente: CoberturaPeriodoA, AreaPeriodoA, CoberturaPeriodoB y AreaPeriodoB.

3. Intersecte los shapefiles, con lo cual podía interpretar que si los campos CoberturaPeriodoA y CoberturaPeriodoB de la capa resultante tenían el mismo tipo, por ejemplo bosques, no existía cambio de cobertura  para los polígonos que presentaran  esa característica (pero si podía existir aumento o disminución del mismo tipo de cobertura, al comparar los campos de Área), si presentaban un tipo diferente (por ejemplo bosque y cultivos) había existido un cambio de cobertura y a través de los campos de área podía ver de cuanto había sido. Siguiendo este proceso podía cuantificar los cambios de cobertura entre los dos periodos de tiempo, pero al igual que con el resultado de la resta de las imágenes me ocurría que al seleccionar todos los registros de una cobertura en particular (Bosque del campo CoberturaPeriodoA) el área no coincidía con el área total presente en el shapefile original (CobA).

Algo que puede estar influyendo en esto puede ser que los shapefiles CobA y CobB no presentan la misma área, el proceso que se llevó a cabo para su obtención fue tomar las imágenes y cortarlas con respecto a una zona de interés, utilizando las herramientas extract by mask y Clip, con ambas herramientas se obtienen áreas totales diferentes, lo cual no debería ser así si se están cortando con el mismo shapefile.

Me gustaría saber si existe una mejor forma de determinar y cuantificar los cambios de cobertura con arcgis, o si los procedimientos que estoy realizando son los correctos para obtener lo que necesito y de ser así que errores estoy cometiendo.

Gracias.

1 Respuesta

+1 voto
por Administrador (7.4k puntos)  

Para realizar un análisis multitemporal, personalmente prefiero trabajar en imágenes ráster, si trabajas con vectores, puedes crear tres campos: Antes, Ahora y Cambio, en este último concatenas los dos primeros campos, categorizas y posteriormente realizas las operaciones que requieras, ejemplo:

Cambio de cobertura

  1. Cultivos
  2. Pastizales
  3. Bosque

Las categorías serían:

11, 22, 33. Sin cambios

21. Pasto a Cultivo

13. Cultivo a Bosque

32. Bosque a Pasto

12. Cultivo a Pasto

Si usas Idrisi Selva, tiene el módulo Land Change Modeler, es perfecto para detectar cambios en imágenes ráster (con vectores puedes analizar dos imágenes, en ráster múltiples), y realizar proyecciones, más que todo puedes agregar otras variables, también te dejo un enlace de un curso de uso de la tierra.

http://adaptacion.org/Educacion/otros-cursos/Curso-Intensivo-en-determinacion-de-Uso-de-la-Tierra-utilizando-imagenes-satelitales

por Novato (106 puntos)  
Muchas gracias, voy a probar las recomendaciones dadas y posteriormente comentare mis resultados
Ofertas de trabajo SIG

Un buen usuario

  1. Comparte lo que sabe
  2. Demuestra su amabilidad
  3. Hace preguntas claras
  4. Escribe correctamente
  5. No hace daño a nadie
Gana 500 puntos y obtén un Certificado SIG validado por ArcGeek.

 

...