[R-br] Conversão de matriz em raster georeferenciado
Éder Comunello
comunello.eder em gmail.com
Terça Janeiro 19 11:41:32 BRST 2016
Alterando grid para representar centro dos pixels...
###
## 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) ### Lower Left
p1 <- p0 + c(dCol/2, dLin/2)
grid <- expand.grid(x=seq(p1[1], by=dCol, len=nCol),
y=seq(p1[2], by=dLin, len=nLin))
plot(grid, 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(grid, 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 19 de janeiro de 2016 09:43, ASANTOS <alexandresantosbr em yahoo.com.br>
escreveu:
> Muito obrigado Éder,
>
> Ficou muito bom mas era quase isso, pois a coordenada do grid1
> deveria marcar o centro do pixel e não os vértices
>
> --
> ======================================================================
> 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
> ======================================================================
>
> Em 19/01/2016 09:28, Éder Comunello escreveu:
>
> 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.
>
>
>
>
> _______________________________________________
> R-br mailing listR-br em listas.c3sl.ufpr.brhttps://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo m�nimo reproduz�vel.
>
>
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160119/4557f4d5/attachment.html>
Mais detalhes sobre a lista de discussão R-br