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.