[R-br] Conversão de matriz em raster georeferenciado

[Éder Comunello]

Alexandre, bom dia!

Com seu código eu obtenho:
> res(r)
[1] 1.978261 2.823529

Seria esperado 2 e 3. O problema parece estar na forma de associar seu grid
com os limites do raster.

Sugiro algo assim:

###
## Dados simulados
dCol <- 3; nCol <- 10
dLin <- 2; nLin <- 10
injuria <- matrix(rnorm(nCol*nLin, mean=12), nLin, nCol)
dim(injuria)

## Coordenadas do grid espaçadas 3 x 2
p0 <- c(680200,7803300)
grid0 <- expand.grid(x=seq(p0[1], by=dCol, len=nCol), y=seq(p0[2], by=dLin,
len=nLin))
grid1 <- expand.grid(x=seq(p0[1], by=dCol, len=nCol+1), y=seq(p0[2],
by=dLin, len=nLin+1))
plot(grid0, asp=T, pch=21)

require(raster)
r <-raster(injuria, xmn=p0[1], xmx=p0[1]+(nCol*dCol),
                    ymn=p0[2], ymx=p0[2]+(nLin*dLin),
                    crs=CRS("+proj=utm +zone=23 +south +datum=WGS84
+units=m +no_defs"))

plot(r); points(grid0, pch=21)
plot(r); points(grid1, pch=21)
r
###





​
================================================
Éder Comunello
PhD Student in Agricultural Systems Engineering (USP/Esalq)
Brazilian Agricultural Research Corporation (Embrapa)
Dourados, MS, Brazil [22 16.5'S, 54 49.0'W]




Em 18 de janeiro de 2016 15:44, ASANTOS <alexandresantosbr em yahoo.com.br>
escreveu:

> Prezados Listeiros,
>
>
>              Gostaria de realizar a conversão de uma matriz em raster
> georeferenciado com duas informações, a representação de plantas em uma
> matriz (B) e juntar as informações geográficas conhecendo-se o espaçamento
> entre plantas (d1), no entanto, gostaria de saber se conforme o meu CRM
> abaixo essa conversão esta correta, pois não estou bem certo de que
> realização de um cbind() entre os objetos B e d1 correspondam a real
> posição desses elementos. Alguém poderia me ajudar?
>
>
> ##Dados artificiais
> injuria<-rnorm(1550, mean=12)
>
> ##Crio a matriz com os dados artificiais
> B = matrix(
>     data=injuria,
>     nrow=25,
>     ncol=62)
>
> # Inserindo coordenadas espaciais
> p0<-c(680234.482625,7803361.55128)# Ponto inicial
> c0 <- seq(p0[1],p0[1]+182,by=2)
> d0 <- seq(p0[2],p0[2]+50,by=3)
> d1 <- expand.grid(x=c0, y=d0) ## Coordenadas do grid espaçadas 3 x 2
> d1<-d1[-(1:14),] ## Removo algumas corrdenadas para ficar com o mesmo
> número de dados
> plot(d1[,1],d1[,2])
> length(d1[,1])
>
> require(raster)
> FORM<-cbind(d1,B)
> colnames(FORM)<-c('x','y','area')
> r <- rasterFromXYZ(as.data.frame(FORM)[, c("x", "y", "area")])
> extent(r) <- c(680234.5,680416.5,7803362,7803410)
> projection(r) <- CRS("+proj=utm +zone=23 +south +datum=WGS84 +units=m
> +no_defs")
> plot(r)
> #
>
>
>
> Obrigado,
>
> --
> ======================================================================
> Alexandre dos Santos
> Proteção Florestal
> IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso
> Campus Cáceres
> Caixa Postal 244
> Avenida dos Ramires, s/n
> Bairro: Distrito Industrial
> Cáceres - MT                      CEP: 78.200-000
> Fone: (+55) 65 8132-8112 (TIM)   (+55) 65 9686-6970 (VIVO)
> e-mails:alexandresantosbr em yahoo.com.br
>         alexandre.santos em cas.ifmt.edu.br
> Lattes: http://lattes.cnpq.br/1360403201088680
> OrcID: orcid.org/0000-0001-8232-6722
> Researchgate: https://www.researchgate.net/profile/Alexandre_Santos10
> LinkedIn: https://br.linkedin.com/in/alexandre-dos-santos-87961635
> ======================================================================
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e fornea cdigo
> mnimo reproduzvel.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160119/8b940a68/attachment.html>