> On Thu, 22 Jun 2006, Andy Bunn wrote: > > > Howdy all: I have a simple question. Given two sp objects, a > > SpatialGridDataFrame tdmean and SpatialLinesDataFrame states I > want to plot > > tdmean as an image and draw the lines on top. Can this be done with sp? > > > > > > > class(tdmean) > > [1] "SpatialGridDataFrame" > > attr(,"package") > > [1] "sp" > > > is.projected(tdmean) > > [1] TRUE > > > class(states) > > [1] "SpatialLinesDataFrame" > > attr(,"package") > > [1] "sp" > > > is.projected(states) > > [1] TRUE > > > # plot states and put the image on top: works but is unsatisfying > > > plot(states,axes=T) > > > image(tdmean["z"],add=T) > > > > > > # plot image and put the state lines on top: doesn't work > > > image(tdmean["z"],axes=T) > > > lines(states) > > > > > Cany you please try: > > image(tdmean, "z", axes=T) > plot(states, add=TRUE) > > I don't think there is a lines() method for SpatialLines* objects.
This works. Thanks. Andy > > Roger > > > > > > _______________________________________________ > > R-sig-Geo mailing list > > [email protected] > > https://stat.ethz.ch/mailman/listinfo/r-sig-geo > > > > -- > Roger Bivand > Economic Geography Section, Department of Economics, Norwegian School of > Economics and Business Administration, Helleveien 30, N-5045 Bergen, > Norway. voice: +47 55 95 93 55; fax +47 55 95 95 43 > e-mail: [EMAIL PROTECTED] > > _______________________________________________ R-sig-Geo mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-sig-geo
