
20 Jun
2013
20 Jun
'13
17:21
Oi Alex, Veja este exemplo aplicado ao raster do exemplo da funcao citada. require(raster) example(rasterToPolygons) pol@data$i <- col(matrix(1L,36,18))[which(r@data@values>6)] pol@data$j <- row(matrix(1L,36,18))[which(r@data@values>6)] head(pol@data, 3) r[1, head(pol@data, 3)$j] tail(pol@data,3) r[18, tail(pol@data,3)$j]