[R-br] maptools

[Elias T. Krainski]

Oi João Augusto, 

A função readShapePoly() possui uma limitação, segundo: 
   https://stat.ethz.ch/pipermail/r-sig-geo/2011-July/012358.html

Com readOGR() eu consegui. Se vc baixou os mapas de RJ? Se foi do diretorio 
     ftp://geoftp.ibge.gov.br/malhas_digitais/censo_2010/setores_censitarios/ 
então não há '33MEE250GC_SIR', há estes quatro conjuntos: 
  listf <- unique(substr(dir("mapas/rj"), 1, 14))
  listf
[1] "33DSE250GC_SIR" "33MUE250GC_SIR" "33SDE250GC_SIR"
[4] "33SEE250GC_SIR"
(talvez vc trocou S por M)

veja o resultado em cada um
> rj <- readOGR("mapas/rj", listf[1])
OGR data source with driver: ESRI Shapefile 
Source: "mapas/rj", layer: "33DSE250GC_SIR"
with 293 features and 3 fields
Feature type: wkbPolygon with 2 dimensions
> rj <- readOGR("mapas/rj", listf[2])
OGR data source with driver: ESRI Shapefile 
Source: "mapas/rj", layer: "33MUE250GC_SIR"
with 92 features and 3 fields
Feature type: wkbPolygon with 2 dimensions
> rj <- readOGR("mapas/rj", listf[3])
OGR data source with driver: ESRI Shapefile 
Source: "mapas/rj", layer: "33SDE250GC_SIR"
with 381 features and 3 fields
Feature type: wkbPolygon with 2 dimensions
> rj <- readOGR("mapas/rj", listf[4])
OGR data source with driver: ESRI Shapefile 
Source: "mapas/rj", layer: "33SEE250GC_SIR"
with 28318 features and 13 fields
Feature type: wkbPolygon with 2 dimensions

Att. 
Elias T. Krainski


>________________________________
> De: João Augusto P Batista <joao.augustopb em gmail.com>
>Para: r-br em listas.c3sl.ufpr.br 
>Enviadas: Sexta-feira, 12 de Outubro de 2012 21:06
>Assunto: [R-br] maptools
> 
>Olá!
>
>Esta é minha primeira mensagem, vou tentar fazê-la o mais inteligível possível.
>
>Estou tentando utilizar os mapas do IBGE, através da função readShapePoly(), mas sempre dá erro. Os arquivos na pasta dos mapas tem extensão .dbf, .shp, .shx, .prj.
>Segue a rotina e o  output:
>
>> library(descr)
>> library(Hmisc)
>> library(foreign)
>> library(maptools)
>> mapa <- readShapePoints("Mapas IBGE/Municípios/33MEE250GC_SIR")
>Shapefile type: PolygonM (25), # of Shapes: 6
>
>Error in read.shape(filen = fn, verbose = verbose, repair = repair) : 
>  Shapefile type not (yet) handled by this function
>
>Já procurei nas mensagens anteriores, mas nenhuma respondeu o meu problema, e já tentei readOGR() e readShapePoints(), sem sucesso.
>
>O que fazer?
>
>Muito obrigado, 
>
>João Augusto.
>_______________________________________________
>R-br mailing list
>R-br em listas.c3sl.ufpr.br
>https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
>
>
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121013/312ef792/attachment.html>