Hi everyone,

I am just starting to use Spatial Data Analysis with R and I tumbled right away into a problem: I am not able to produce a plot of a SpatialPolygonsDataFrame on my MacOsX which is no problem on a Windows computer. It seems that the problem only arises with large .shp files. For example, I have obtained .shp files from Luxembourg and the USA from http://www.diva-gis.org. With the following code I am able to load .shp files into R, but I can only produce a plot of Luxembourg when I try to plot the US R crashes. I know that with map{maps} I can easily produce a map of the USA, however I need the USA and its administrative units as a polygon to further work with (extract values from raster object for example) and I would also like to plot my results....so

is there a way to solve the Quartz problem on my Mac?
or is there another way to obtain a polygon of the USA that I can work with?

library(maptools)
setwd("/")
lux<-readShapePoly("LUX_adm0.shp")
plot(lux)

setwd("/")
library(maptools)
us<-readShapePoly("USA_adm0.shp")
plot(us)

thank you so much in advance :))

Heike



	[[alternative HTML version deleted]]

