La fórmula a emplear para obtener el índice de vegetación, usando el comando r.mapcalc de GRASS,
sería r.mapcalc
'ndvi=100*(r4-r3)/(r4+r3)'
Posteriormente se deben reclasificar los valores menores de cero,
pues no tienen sentido. Esta operación puede realizarse con el comando r.mapcalc
'ndvi2=if(ndvi<0,0,ndvi)' con lo que obtendríamos la capa ndvi2
que nos muestra el índice de vegetación con valores iguales o
superiores a 0.
Índice de vegetación
|
Índice de vegetación modificado
|
2. Análisis de Componentes Principales. Haz el ACP de tu zona y evalúa los resultados.
acp.1 |
acp.2 |
acp.3 |
acp.4 |
acp.5 |
acp.6 que corresponde a la b7 |
Resultado de ejecutar i.pac input=b1,b2,b3,b4,b5,b7
output=acp , lo que genera estas seis capas con el prefijo
acp y como extensión el número del componente analizado. Se observa una mayor definición e información en la capa acp.2 lo que puede denotar un mayor diferencial en la varianza absorbida lo que destaca la importancia relativa de cada componente. |