<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
Prof. Paulo,<br>
<br>
Seguindo seus conselhos ate consigo criar o objeto
SpatialPixelsDataFrame, mas o data não aparece e se faço
img.kc<-raster(kc_output) e ploto, consigo visualizar o grid
krigado, sendo CRM:<br>
<br>
#<br>
require(geoR)<br>
require(sp)<br>
require(raster)<br>
plot(s100$coords, xlim=c(0,1.2), ylim=c(0,1.2), xlab="Coord X",
ylab="Coord Y")<br>
loci <- matrix(c(0.2, 0.6, 0.2, 1.1, 0.2, 0.3, 1.0, 1.1), ncol=2)<br>
text(loci, as.character(1:4), col="red")<br>
polygon(x=c(0,1,1,0), y=c(0,0,1,1), lty=2)<br>
<br>
# Fitting models with nugget fixed to zero<br>
ml <- likfit(s100, ini = c(1,0.5), fix.nugget = T)<br>
<br>
# defining the grid<br>
pred.grid <- expand.grid(seq(0,1, l=51), seq(0,1, l=51))<br>
# kriging calculations<br>
kc <- krige.conv(s100, loc = pred.grid, krige =
krige.control(obj.m = ml))<br>
# displaying predicted values<br>
image(kc, loc = pred.grid, col=gray(seq(1,0.1,l=30)), xlab="Coord
X", ylab="Coord Y")<br>
#<br>
<br>
#Centro das 25 parcela 0.2 x 0.2 <br>
p0<-c(0.1,0.1)<br>
#<br>
c0 <- seq(p0[1],p0[1]+0.9,by=0.2)<br>
d0 <- seq(p0[2],p0[2]+0.9,by=0.2)<br>
d1 <- expand.grid(x=c0, y=d0)<br>
points(d1[,1],d1[,2],col="red")<br>
ver<-c(0.0,0.2,0.4,0.6,0.8,1.0)<br>
hor<-c(0.0,0.2,0.4,0.6,0.8,1.0)<br>
abline(v=ver, col="red")<br>
abline(h=hor, col="red")<br>
#<br>
grid=GridTopology(c(0.1,0.1), c(0.2,0.2), c(5,5))<br>
grid=SpatialGrid(grid=grid)<br>
class(grid) ## SpatialPixels<br>
gridPol <- as(grid, "SpatialPolygons") <br>
options(digits=5)<br>
pixelsize <- ((1.0-0.0)/30)*((1.0-0.0)/30)#Tamanho do pixel é
dado por seq Xmaior - Xmenor/l * Ymaior-Ymenor/l<br>
#<br>
coordinates(pred.grid) <- ~ Var1+Var2## Transformando as
coordenadas do grid<br>
kc$predict<-as.data.frame(kc$predict)### Valores dos pixels
preditos<br>
kc_output<-SpatialPixelsDataFrame(points=pred.grid,
data=kc$predict)### Criando o objeto SpatialPixelsDataFrame<br>
summary(kc_output)## Não aparece o data mas:<br>
<br>
#<br>
<br>
<br>
Obrigado,<br>
<br>
<br>
Alexandre<br>
<div class="moz-cite-prefix">Em 08/04/2013 15:55, Paulo Justiniano
escreveu:<br>
</div>
<blockquote
cite="mid:alpine.DEB.2.00.1304081654110.13740@pataxo.est.ufpr.br"
type="cite">como nao?< é um veror numerico de valores!
<br>
<br>
a melhor solicao é fazer uma funcao que converta a saida de
kriging para o formato sp como SpatialPoitsDataFrame (ou Pixel ou
Grid - dependendo da malha de predicao)
<br>
<br>
<br>
vou ver se tenho aogo aqui, mas isto é simples de qq forma
<br>
<br>
<br>
On Mon, 8 Apr 2013, ASANTOS wrote:
<br>
<br>
<blockquote type="cite">Thiago,
<br>
<br>
Vou continuar tentando, mas extract() não funciona, pois o
objeto kc$predict pertence a classe numeric.
<br>
<br>
Abraço,
<br>
<br>
Alexandre
<br>
<br>
<br>
Em 07/04/2013 02:54, Thiago V. dos Santos escreveu:
<br>
Alexandre,
<br>
<br>
Você já tem o polígono gridPol de onde você quer criar a média
da krigagem. Não conheço bem as classes do pacote geoR, mas as
<br>
coordenadas do kc$predict parecem estar no objeto pred.grid. A
abordagem que eu usaria seria criar um objeto raster da krigagem
e em
<br>
seguida contar os valores usando a função extract (raster)
usando como argumentos o raster da krigagem e o gridPol.
<br>
<br>
Além da ajuda para a função extract no manual do raster, talvez
o último exemplo dessa página possa te
<br>
ajudar: <a class="moz-txt-link-freetext" href="http://stackoverflow.com/questions/14682606/extract-value-from-raster-stack-from-spatialpolygondataframe">http://stackoverflow.com/questions/14682606/extract-value-from-raster-stack-from-spatialpolygondataframe</a>
<br>
<br>
Mas não estou certo se era isso que você estava buscando. Acho
que o Paulo vai te ajudar melhor.
<br>
</blockquote>
<br>
<br>
<blockquote type="cite">
<br>
Saudações,
<br>
--
<br>
Thiago V. dos Santos
<br>
PhD student
<br>
Land and Atmospheric Science
<br>
University of Minnesota
<br>
<a class="moz-txt-link-freetext" href="http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/index.htm">http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/index.htm</a>
<br>
Phone: (612) 323 9898
<br>
<br>
____________________________________________________________________________________________________________________________________________
<br>
From: ASANTOS <a class="moz-txt-link-rfc2396E" href="mailto:alexandresantosbr@yahoo.com.br"><alexandresantosbr@yahoo.com.br></a>
<br>
To: <a class="moz-txt-link-abbreviated" href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>; Thiago V. dos Santos
<a class="moz-txt-link-rfc2396E" href="mailto:thi_veloso@yahoo.com.br"><thi_veloso@yahoo.com.br></a>
<br>
Sent: Saturday, April 6, 2013 10:26 PM
<br>
Subject: Re: [R-br] Obter valor médio dos pixels de uma imagem
gerada por krigagem
<br>
<br>
Thiago,
<br>
<br>
Já pensei isto, mas para SpatialPoints(vp), preciso fazer um
list() das coordenadas dos pixels (vp$x e vp$y) e isto não tenho
e
<br>
nem sei como recuperar o valor de um pixel à partir de uma
coordenada conhecida.
<br>
<br>
Abraço,
<br>
<br>
Alexandre
<br>
<br>
<br>
Em 06/04/2013 21:32, Thiago V. dos Santos escreveu:
<br>
Nao conheco bem todas as funcoes do sp, mas para pela
mensagem de erro eu transformaria a tabela vp em SpatialPoints,
que
<br>
eh o argumento esperado pela funcao over.
<br>
<br>
Saudações,
<br>
--
<br>
Thiago V. dos Santos
<br>
PhD student
<br>
Land and Atmospheric Science
<br>
University of Minnesota
<br>
<a class="moz-txt-link-freetext" href="http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/index.htm">http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/index.htm</a>
<br>
Phone: (612) 323 9898
<br>
<br>
____________________________________________________________________________________________________________________________________________
<br>
From: ASANTOS <a class="moz-txt-link-rfc2396E" href="mailto:alexandresantosbr@yahoo.com.br"><alexandresantosbr@yahoo.com.br></a>
<br>
To: <a class="moz-txt-link-abbreviated" href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>
<br>
Sent: Saturday, April 6, 2013 8:12 PM
<br>
Subject: Re: [R-br] Obter valor médio dos pixels de uma imagem
gerada por krigagem
<br>
<br>
Boa noite pessoal,
<br>
<br>
Acho evoluí um pouco na questão, mas agora a dificuldade
esta em contar os pixels contidos no interior das malhas, quando
<br>
tento um over(), ocorre o seguinte erro:
<br>
<br>
Erro em function (classes, fdef, mtable) :
<br>
unable to find an inherited method for function "over", for
signature "table", "SpatialPolygons"
<br>
<br>
<br>
segue CRM:
<br>
#
<br>
require(geoR)
<br>
plot(s100$coords, xlim=c(0,1.2), ylim=c(0,1.2), xlab="Coord X",
ylab="Coord Y")
<br>
loci <- matrix(c(0.2, 0.6, 0.2, 1.1, 0.2, 0.3, 1.0, 1.1),
ncol=2)
<br>
text(loci, as.character(1:4), col="red")
<br>
polygon(x=c(0,1,1,0), y=c(0,0,1,1), lty=2)
<br>
<br>
# Fitting models with nugget fixed to zero
<br>
ml <- likfit(s100, ini = c(1,0.5), fix.nugget = T)
<br>
<br>
# defining the grid
<br>
pred.grid <- expand.grid(seq(0,1, l=51), seq(0,1, l=51))
<br>
# kriging calculations
<br>
kc <- krige.conv(s100, loc = pred.grid, krige =
krige.control(obj.m = ml))
<br>
# displaying predicted values
<br>
image(kc, loc = pred.grid, col=gray(seq(1,0.1,l=30)),
xlab="Coord X", ylab="Coord Y")
<br>
#
<br>
<br>
#Centro das 25 parcela 0.2 x 0.2
<br>
p0<-c(0.1,0.1)
<br>
#
<br>
c0 <- seq(p0[1],p0[1]+0.9,by=0.2)
<br>
d0 <- seq(p0[2],p0[2]+0.9,by=0.2)
<br>
d1 <- expand.grid(x=c0, y=d0)
<br>
points(d1[,1],d1[,2],col="red")
<br>
ver<-c(0.0,0.2,0.4,0.6,0.8,1.0)
<br>
hor<-c(0.0,0.2,0.4,0.6,0.8,1.0)
<br>
abline(v=ver, col="red")
<br>
abline(h=hor, col="red")
<br>
#
<br>
require(sp)
<br>
grid=GridTopology(c(0.1,0.1), c(0.2,0.2), c(5,5))### Cria as
área
<br>
grid=SpatialGrid(grid=grid)
<br>
class(grid) ## SpatialPixels
<br>
gridPol <- as(grid, "SpatialPolygons") ## Transfoma as areas
em polígonos espaciais
<br>
options(digits=5)
<br>
pixelsize <- (1.0-0.0/30)*(1.0-0.0/30)#Tamanho do pixel é
dado por seq Xmaior - Xmenor/l * Ymaior-Ymenor/l
<br>
vp <- table(area <- kc$predict)## Valor dos pixels
<br>
count.pix<-table(over(vp, gridPol)) #contar quantos pontos
tem em cada area
<br>
<br>
<br>
Obrigado,
<br>
<br>
Alexandre
<br>
<br>
Em 06/04/2013 17:11, ASANTOS escreveu:
<br>
> Boa noite membros da lista,
<br>
>
<br>
> Gostaria de receber ideias para calcular o valor
médio dos pixel em uma imagem gerada por krigagem em sub-regiões
de
<br>
interesse na imagem, por exemplo usando um exemplo do Prof.
Paulo, onde ajusta-se um semivariograma e faz a krigagem
ordinária
<br>
dos dados:
<br>
>
<br>
> #
<br>
> require(geoR)
<br>
> plot(s100$coords, xlim=c(0,1.2), ylim=c(0,1.2), xlab="Coord
X", ylab="Coord Y")
<br>
> loci <- matrix(c(0.2, 0.6, 0.2, 1.1, 0.2, 0.3, 1.0,
1.1), ncol=2)
<br>
> text(loci, as.character(1:4), col="red")
<br>
> polygon(x=c(0,1,1,0), y=c(0,0,1,1), lty=2)
<br>
>
<br>
> # Fitting models with nugget fixed to zero
<br>
> ml <- likfit(s100, ini = c(1,0.5), fix.nugget = T)
<br>
>
<br>
> # defining the grid
<br>
> pred.grid <- expand.grid(seq(0,1, l=51), seq(0,1,
l=51))
<br>
> # kriging calculations
<br>
> kc <- krige.conv(s100, loc = pred.grid, krige =
krige.control(obj.m = ml))
<br>
> # displaying predicted values
<br>
> image(kc, loc = pred.grid, col=gray(seq(1,0.1,l=30)),
xlab="Coord X", ylab="Coord Y")
<br>
> #
<br>
>
<br>
> #Agora vou dividir esta área em 25 partes e lançar os
centroides destas subáreas:
<br>
>
<br>
> #Centro das 25 parcela 0.2 x 0.2
<br>
> p0<-c(0.1,0.1)
<br>
> #
<br>
> c0 <- seq(p0[1],p0[1]+0.9,by=0.2)
<br>
> d0 <- seq(p0[2],p0[2]+0.9,by=0.2)
<br>
> d1 <- expand.grid(x=c0, y=d0)
<br>
> points(d1[,1],d1[,2],col="red")
<br>
> ver<-c(0.0,0.2,0.4,0.6,0.8,1.0)
<br>
> hor<-c(0.0,0.2,0.4,0.6,0.8,1.0)
<br>
> abline(v=ver, col="red")
<br>
> abline(h=hor, col="red")
<br>
> #
<br>
>
<br>
> E aqui vem minha pergunta, como calcular o valor médio dos
pixels contidos em cada subárea de 0.2 x 0.2?
<br>
>
<br>
> Obrigado,
<br>
>
<br>
<br>
--
======================================================================
<br>
Alexandre dos Santos
<br>
Proteção Florestal
<br>
Coordenador do curso Técnico em Florestas
<br>
Vice Coordenador do curso de Engenharia Florestal
<br>
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de
Mato Grosso
<br>
Campus Cáceres
<br>
Caixa Postal 244
<br>
Avenida dos Ramires, s/n
<br>
Bairro: Distrito Industrial
<br>
Cáceres - MT CEP: 78.200-000
<br>
Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO)
<br>
<a class="moz-txt-link-abbreviated" href="mailto:e-mails:alexandresantosbr@yahoo.com.br">e-mails:alexandresantosbr@yahoo.com.br</a>
<br>
<a class="moz-txt-link-abbreviated" href="mailto:alexandre.santos@cas.ifmt.edu.br">alexandre.santos@cas.ifmt.edu.br</a>
<br>
======================================================================
<br>
<br>
_______________________________________________
<br>
R-br mailing list
<br>
<a class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<br>
<a class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
<br>
Leia o guia de postagem (<a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e
forneça código mínimo reproduzível.
<br>
<br>
<br>
<br>
<br>
_______________________________________________
<br>
R-br mailing list
<br>
<a class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<br>
<a class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
<br>
Leia o guia de postagem (<a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e
forneça código mínimo reproduzível.
<br>
<br>
<br>
-- <br>
======================================================================
<br>
Alexandre dos Santos
<br>
Proteção Florestal
<br>
Coordenador do curso Técnico em Florestas
<br>
Vice Coordenador do curso de Engenharia Florestal IFMT -
Instituto Federal de Educação, Ciência e Tecnologia de Mato
Grosso
<br>
Campus Cáceres
<br>
Caixa Postal 244
<br>
Avenida dos Ramires, s/n
<br>
Bairro: Distrito Industrial Cáceres - MT
CEP: 78.200-000
<br>
Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO)
<br>
<a class="moz-txt-link-abbreviated" href="mailto:e-mails:alexandresantosbr@yahoo.com.br">e-mails:alexandresantosbr@yahoo.com.br</a>
<br>
<a class="moz-txt-link-abbreviated" href="mailto:alexandre.santos@cas.ifmt.edu.br">alexandre.santos@cas.ifmt.edu.br</a>
======================================================================
<br>
<br>
<br>
<br>
-- <br>
======================================================================
<br>
Alexandre dos Santos
<br>
Proteção Florestal
<br>
Coordenador do curso Técnico em Florestas
<br>
Vice Coordenador do curso de Engenharia Florestal IFMT -
Instituto Federal de Educação, Ciência e Tecnologia de Mato
Grosso
<br>
Campus Cáceres
<br>
Caixa Postal 244
<br>
Avenida dos Ramires, s/n
<br>
Bairro: Distrito Industrial Cáceres - MT
CEP: 78.200-000
<br>
Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO)
<br>
<a class="moz-txt-link-abbreviated" href="mailto:e-mails:alexandresantosbr@yahoo.com.br">e-mails:alexandresantosbr@yahoo.com.br</a>
<br>
<a class="moz-txt-link-abbreviated" href="mailto:alexandre.santos@cas.ifmt.edu.br">alexandre.santos@cas.ifmt.edu.br</a>
======================================================================
<br>
<br>
<br>
</blockquote>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
R-br mailing list
<a class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<a class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</blockquote>
<br>
<pre class="moz-signature" cols="72">--
======================================================================
Alexandre dos Santos
Proteção Florestal
Coordenador do curso Técnico em Florestas
Vice Coordenador do curso de Engenharia 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)
<a class="moz-txt-link-abbreviated" href="mailto:e-mails:alexandresantosbr@yahoo.com.br">e-mails:alexandresantosbr@yahoo.com.br</a>
<a class="moz-txt-link-abbreviated" href="mailto:alexandre.santos@cas.ifmt.edu.br">alexandre.santos@cas.ifmt.edu.br</a>
======================================================================
</pre>
</body>
</html>