[R-sig-Geo] Dramatically slow plotting
MacQueen, Don
macqueen1 at llnl.gov
Thu Mar 20 16:44:34 CET 2014
I'm not sure this will help, but try
usePolypath=FALSE
in your call to spplot().
For further information, see
help.search('polypath')
I get:
Help files with alias or concept or title matching 'polypath' using fuzzy
matching:
graphics::polypath Path Drawing
sp::SpatialPolygons-class Class "SpatialPolygons"
-Don
--
Don MacQueen
Lawrence Livermore National Laboratory
7000 East Ave., L-627
Livermore, CA 94550
925-423-1062
On 3/17/14 3:23 PM, "Rolando Valdez" <rvaldezr at gmail.com> wrote:
>Hello,
>
>Recently, I acquired a MacBook Pro, Core i7, 8 GB ram. I Installed the
>newest R version, 3.0.3 from the web page. The problem is when I¹m
>plotting maps, because is going very, very slow, about 3 or 4 minutes
>just for a single map, while I¹ve done this in a few seconds in Windows
>with Core i5 and 4 GB ram.
>
>This is what I have:
>
>R version 3.0.3 (2014-03-06) -- "Warm Puppy"
>Copyright (C) 2014 The R Foundation for Statistical Computing
>Platform: x86_64-apple-darwin10.8.0 (64-bit)
>
>[R.app GUI 1.63 (6660) x86_64-apple-darwin10.8.0]
>
>I found a reproducible example in web and I took time with proc.time()
>
>ptm <- proc.time()
>library(sp)
>library(lattice) # required for trellis.par.set():
>trellis.par.set(sp.theme()) # sets color ramp to bpy.colors()
>
># prepare nc sids data set:
>library(maptools)
>nc <- readShapePoly(system.file("shapes/sids.shp",
>package="maptools")[1], proj4string=CRS("+proj=longlat +datum=NAD27"))
>arrow = list("SpatialPolygonsRescale", layout.north.arrow(),
> offset = c(-76,34), scale = 0.5, which = 2)
>#scale = list("SpatialPolygonsRescale", layout.scale.bar(),
># offset = c(-77.5,34), scale = 1, fill=c("transparent","black"),
>which = 2)
>#text1 = list("sp.text", c(-77.5,34.15), "0", which = 2)
>#text2 = list("sp.text", c(-76.5,34.15), "1 degree", which = 2)
>## multi-panel plot with filled polygons: North Carolina SIDS
>spplot(nc, c("SID74", "SID79"), names.attr = c("1974","1979"),
> colorkey=list(space="bottom"), scales = list(draw = TRUE),
> main = "SIDS (sudden infant death syndrome) in North Carolina",
> sp.layout = list(arrow), as.table = TRUE)
>
># sp.layout = list(arrow, scale, text1, text2), as.table = TRUE)
>proc.time() - ptm
>
>user system elapsed
> 2.408 0.064 2.616
>
>It was quick.
>
>Then I did a single plot with my shape:
>
>mapa <- readShapePoly(³Entidades_2013.shp²)
>ptm <- proc.time()
>spplot(mapa[1]); proc.time() - ptm
>
>user system elapsed
> 87.575 0.786 88.068
>
>Why it take a lot of time? I worked with same shapes in Windows and never
>took that time.
>
>Hope you can help me,
>
>Regards,
>
>Rolando Valdez
>
>_______________________________________________
>R-sig-Geo mailing list
>R-sig-Geo at r-project.org
>https://stat.ethz.ch/mailman/listinfo/r-sig-geo
More information about the R-sig-Geo
mailing list