[R-br] Fwd: Mover coordenadas em shapefile

[Éder Comunello]

Alexandre, bom dia!

Você pode manipular o objeto 'sp', resultante da importação, com a função
elide() do pacote {maptools}. Se precisar de um novo shapefile, utilize a
função writeOGR() ou similar.

### <code r>
setwd("C:/LAB/RGIS/coords/elide")
lapply(c('rgdal', 'maptools', 'shapefiles'), require, character.only=TRUE)

#Download do controno da área em shapefile
links <- c(
"http://dl.dropbox.com/s/1oz5lhnk96jih3w/Tamandua288.dbf",
"http://dl.dropbox.com/s/6b3uygyrjhq1bby/Tamandua288.shp",
"http://dl.dropbox.com/s/wntluhn8bi7us2f/Tamandua288.shx")

if (sum(sapply(basename(links), file.exists))>3) {
          lapply(links, function(a) tryCatch(download.file(a,
dest=basename(a), mode='wb'),
                                              error=function(...)
message("ERRO!")))} else message('OK!')

## Abrir e visualizar contorno
CRS.new <- CRS("+proj=utm +zone=23 +south +datum=WGS84 +units=m +no_defs
+ellps=WGS84 +towgs84=0,0,0")
cont <- readOGR(getwd(), "Tamandua288")
cont em bbox
proj4string(cont) <- CRS.new
plot(cont)

cont2 <- elide(cont, shift=c(1032,280))
plot(cont2, col=3, add=T)

ogrDrivers()
writeOGR(cont2, getwd(), 'cont2', 'ESRI Shapefile')
dir(patt='.*shp')
### </code>

Éder Comunello <c <comunello.eder em gmail.com>omunello.eder em gmail.com>
Dourados, MS - [22 16.5'S, 54 49'W]
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140430/c1532bc7/attachment.html>