Determinación de la matriz de varianza-covarianza para bandas ráster: procedimiento resumido en lenguajeR

En un artículo precendente se señaló un procedimiento comprensivo para determinar la matriz de varianza-covarianza para bandas ráster. No obstante, en aras del tiempo, se va a resumir en el menor número de líneas posible y compararlo con el método directo.

El procedimiento comprensivo resumido es el siguiente:

>setwd('Rproject')
>library(rgdal)
>b3<-readGDAL('b3.ND.tif')
b3.ND.tif has GDAL driver GTiff 
and has 2493 rows and 4613 columns
> b4<-readGDAL('b4.ND.tif')
b4.ND.tif has GDAL driver GTiff 
and has 2493 rows and 4613 columns
>#se asignan sus values a dos vectores
>z3<-b3$band1
>z4<-b4$band1
>M<-cbind(z3,z4)  #matriz de vectores
>n <- nrow(M) #numero de filas de M
> M_mean <- matrix(data=1, nrow=n) %*% cbind(mean(z3),mean(z4))  #matriz de medias
>D<-M-M_mean  #matriz de valores menos matriz de medias
>E<- (t(D)%*%D)/n  #matriz de covarianza
> E
          z3        z4
z3 371.10435  29.85423
z4  29.85423 381.56172

El procedimiento directo (usa el método ‘cov’) es el siguiente:

> cov(M)
          z3        z4
z3 371.10438  29.85423
z4  29.85423 381.56175

Como se puede apreciar, los resultados son prácticamente idénticos.

Anuncios
Esta entrada fue publicada en Lenguaje R, SIG. Guarda el enlace permanente.

Responder

Por favor, inicia sesión con uno de estos métodos para publicar tu comentario:

Logo de WordPress.com

Estás comentando usando tu cuenta de WordPress.com. Cerrar sesión / Cambiar )

Imagen de Twitter

Estás comentando usando tu cuenta de Twitter. Cerrar sesión / Cambiar )

Foto de Facebook

Estás comentando usando tu cuenta de Facebook. Cerrar sesión / Cambiar )

Google+ photo

Estás comentando usando tu cuenta de Google+. Cerrar sesión / Cambiar )

Conectando a %s