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 <[email protected]>
>Para: [email protected] 
>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
>[email protected]
>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.
>
>
>
_______________________________________________
R-br mailing list
[email protected]
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.
  • [R-br] maptools João Augusto P Batista
    • Re: [R-br] maptools Elias T. Krainski

Responder a