[R-br] Selecionar coordenadas e valor de pixels contidos no interior de parcelas

[Mauro Sznelwar]

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.
  Alexandre, bom dia!


  Não tenho experiência com a {spatstat}, mas tentei trabalhar seu código e acredito que possa lhe ser útil.



  ### <BEGIN>
  require(raster); require(spatstat); require(sp) 


  ## Raster 
  ----------------------------------------------------------- 
  set.seed(20130622)
  r <- raster(nc=100, nr=100) 
  r <- setValues(r, round(runif(ncell(r))* 255))
  plot(r)


  r.spdf <- data.frame(coordinates(r), px=getValues(r)); names(r.spdf)
  gridded(r.spdf) <- ~x+y
  spplot(r.spdf)


  ## Amostras
  ------------------------------------------------- 
  parc.p <- spsample(r.spdf, 15, type="random")
  parc.c <- vector("list",length(parc.p)) 
  for(i in 1:length(parc.p)) {
            parc.c[[i]] <- disc(radius=5, centre=c(parc.p$x[i], parc.p$y[i]))}


  W    <- do.call(union.owin,parc.c)
  W.df <- as.data.frame(as.im(W))


  W.sp <- W.df; coordinates(W.sp) <- ~x+y
  plot(W.sp)


  W.df$px <- over(W.sp, r.spdf)
  head(W.df)
  #
  #          x         y value  px
  #1 -179.9856 -39.83543     1  11
  #2 -179.9856 -39.29907     1  76
  #3 -179.9856 -38.76271     1  76
  #4 -179.9856 -38.22635     1  76
  #5 -179.9856 -37.68998     1  34
  #6 -178.7263 -41.98087     1 157
  #
  ### <END>


  .
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20131226/b1c160d7/attachment.html>