Dear all,

 

I'm trying to use maptools to produce some figures using the Brazilian
municipalities shape files (5,500 municipalities and ~9mb).

 

I've ran some examples suggested by Roger which worked fine with the
sids.shp demo file.

 

library(maptools)

 

xsp <- readShapePoly(system.file("shapes/sids.shp",package="maptools")[1],
proj4string=CRS("+proj=longlat"))

plot(xsp)

 

brks <- seq(min(xsp$BIR74), max(xsp$BIR74), length=7) 

brks

cols <- colorRampPalette(c("lightpink", "red"))(length(brks)-1) 

cols

 

plot(xsp, col=cols[findInterval(xsp$BIR74, brks, all.inside=TRUE)])

 

 

However, whenever I try to load the Brazilian data, R just stalls, and
nothing happens.

 

brasilsp <-
readShapePoly(system.file("shapes/Brasil_municipios_2001.shp",package="mapto
ols")[1], proj4string=CRS("+proj=longlat +ellps=clrk66"))

 

plot(brasilsp)

 

I've also tried to use the old  read.shape, which work fine the overall
plot.

 

x <- read.shape(system.file("shapes/Brasil_municipios_2001.shp",
package="maptools")[1])

 

plot(x, axes=FALSE, xlab = "", ylab = "")

 

plot(x, axes=FALSE, xlab = "", ylab = "", ol = "transparent")

 

However, if I try to use the auxvar option, I get an error message.

 

plot(x, auxvar=x$att.data$TOTAL, type="e",ol = "transparent", prbg=NULL)

 

Error in plot.Map(x, auxvar = x$att.data$TOTAL, type = "e", ol =
"transparent", :

 

 

Any help will be much appreciated.

 

Cheers,

 

JP

 


	[[alternative HTML version deleted]]


