<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META content="text/html; charset=iso-8859-1" http-equiv=Content-Type>
<META name=GENERATOR content="MSHTML 8.00.6001.23536">
<STYLE></STYLE>
</HEAD>
<BODY bgColor=#ffffff>
<DIV><FONT size=2 face=Arial></FONT><FONT size=2 face=Arial></FONT><BR><FONT
color=#0000ff face=Arial>Eu não consegui baixar a bibliotéca sptsatat, diz que
não está disponível na versão R 3.0.1 do Windows.</FONT></DIV>
<BLOCKQUOTE
style="BORDER-LEFT: #000000 2px solid; PADDING-LEFT: 5px; PADDING-RIGHT: 0px; MARGIN-LEFT: 5px; MARGIN-RIGHT: 0px">
<DIV dir=ltr>Alexandre, bom dia!
<DIV><BR></DIV>
<DIV>Não tenho experiência com a {spatstat}, mas tentei trabalhar seu código e
acredito que possa lhe ser útil.<BR>
<DIV><BR></DIV>
<DIV>
<DIV><FONT face="courier new, monospace">### <BEGIN></FONT></DIV>
<DIV><FONT face="courier new, monospace">require(raster); require(spatstat);
require(sp) </FONT></DIV>
<DIV><FONT face="courier new, monospace"><BR></FONT></DIV>
<DIV><FONT face="courier new, monospace">## Raster </FONT></DIV>
<DIV><FONT
face="courier new, monospace">----------------------------------------------------------- </FONT></DIV>
<DIV><FONT face="courier new, monospace">set.seed(20130622)</FONT></DIV>
<DIV><FONT face="courier new, monospace">r <- raster(nc=100,
nr=100) </FONT></DIV>
<DIV><FONT face="courier new, monospace">r <- setValues(r,
round(runif(ncell(r))* 255))</FONT></DIV>
<DIV><FONT face="courier new, monospace">plot(r)</FONT></DIV>
<DIV><FONT face="courier new, monospace"><BR></FONT></DIV>
<DIV><FONT face="courier new, monospace">r.spdf <-
data.frame(coordinates(r), px=getValues(r)); names(r.spdf)</FONT></DIV>
<DIV><FONT face="courier new, monospace">gridded(r.spdf) <-
~x+y</FONT></DIV>
<DIV><FONT face="courier new, monospace">spplot(r.spdf)</FONT></DIV>
<DIV><FONT face="courier new, monospace"><BR></FONT></DIV>
<DIV><FONT face="courier new, monospace">## Amostras</FONT></DIV>
<DIV><FONT
face="courier new, monospace">------------------------------------------------- </FONT></DIV>
<DIV><FONT face="courier new, monospace">parc.p <- spsample(r.spdf, 15,
type="random")</FONT></DIV>
<DIV><FONT face="courier new, monospace">parc.c <-
vector("list",length(parc.p)) </FONT></DIV>
<DIV><FONT face="courier new, monospace">for(i in 1:length(parc.p))
{</FONT></DIV>
<DIV><FONT face="courier new, monospace">
parc.c[[i]] <- disc(radius=5, centre=c(parc.p$x[i],
parc.p$y[i]))}</FONT></DIV>
<DIV><FONT face="courier new, monospace"><BR></FONT></DIV>
<DIV><FONT face="courier new, monospace">W <-
do.call(union.owin,parc.c)</FONT></DIV>
<DIV><FONT face="courier new, monospace">W.df <- as.data.frame(<A
href="http://as.im">as.im</A>(W))</FONT></DIV>
<DIV><FONT face="courier new, monospace"><BR></FONT></DIV>
<DIV><FONT face="courier new, monospace">W.sp <- W.df; coordinates(W.sp)
<- ~x+y</FONT></DIV>
<DIV><FONT face="courier new, monospace">plot(W.sp)</FONT></DIV>
<DIV><FONT face="courier new, monospace"><BR></FONT></DIV>
<DIV><FONT face="courier new, monospace">W.df$px <- over(W.sp,
r.spdf)</FONT></DIV>
<DIV><FONT face="courier new, monospace">head(W.df)</FONT></DIV></DIV>
<DIV><FONT face="courier new, monospace">#</FONT></DIV>
<DIV><FONT face="courier new, monospace">
<DIV># x y value
px</DIV>
<DIV>#1 -179.9856 -39.83543 1 11</DIV>
<DIV>#2 -179.9856 -39.29907 1 76</DIV>
<DIV>#3 -179.9856 -38.76271 1 76</DIV>
<DIV>#4 -179.9856 -38.22635 1 76</DIV>
<DIV>#5 -179.9856 -37.68998 1 34</DIV>
<DIV>#6 -178.7263 -41.98087 1 157</DIV>
<DIV>#</DIV></FONT></DIV>
<DIV><FONT face="courier new, monospace">###
<END></FONT></DIV></DIV></DIV>
<DIV class=gmail_extra><BR clear=all>
<DIV>.</DIV></DIV></BLOCKQUOTE></BODY></HTML>