On Thursday 20 September 2007, Rick Reeves wrote: > Hello: > > I would like to superimpose vector polygons (state outlines) from a > Shape file on top of a satellite image, > imported into a SpatialGridDataFrame from GEOTIFF via gdal_translate and > readGDAL. > > When I plot polygon and point shape files in R, into > SpatialPointDataFrame and SpatialPolygonDataFrame, > the two feature sets line up geographically, so it seems logical that a > SpatialGridDataFrame should behave > in the same way. > > From my initial research, the spplot function is the correct function > for plotting grids/images with axes and > annotation. > > The big question is, how do I incorporate a Spatial(Points or > Polygons)DataFrame into the spplot display list? > > It seems as though many scientists would like to create such plots > without resorting to GRASS or > another GIS.
While GRASS and GMT are excellent tools, here are some ideas on how to make composite maps in R: http://casoilresource.lawr.ucdavis.edu/drupal/node/442 cheers, Dylan > Thanks for any advice, > Rick Reeves > > Rick Reeves > Scientific Programmer / Analyst > National Center for Ecological Analysis and Synthesis > UC Santa Barbara > [EMAIL PROTECTED] > www.nceas.ucsb.edu > 805 892 2533 -- Dylan Beaudette Soils and Biogeochemistry Graduate Group University of California at Davis 530.754.7341 ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.