Re: [R-sig-Geo] how to plot different rows of a SpatialPolygonsDataFrame in trellis panels

2018-05-31 Thread Waichler, Scott R
, Washington, USA > -Original Message- > From: Kent Johnson [mailto:kent3...@gmail.com] > Sent: Thursday, May 24, 2018 6:05 AM > To: r-sig-geo@r-project.org; Waichler, Scott R > Subject: Re: how to plot different rows of a SpatialPolygonsDataFrame in > trellis panels &

[R-sig-Geo] how to plot different rows of a SpatialPolygonsDataFrame in trellis panels

2018-05-23 Thread Waichler, Scott R
Hello, I have a SpatialPolygonsDataFrame. I would like to do a trellis plot on one of the attributes, so that in the panel for a given attribute value, only those polygons with that value are plotted. So, each panel has different polygons plotted in it. I can't figure out how to do this.

[R-sig-Geo] stConstruct with multiple SpatialPolygonsDataFrame objects

2018-05-21 Thread Waichler, Scott R
Hi, Can anyone provide an example of using stConstruct() with multiple SpatialPolygonsDataFrame objects? I have created the latter by reading in multiple shapefiles, one shapefile per year. Each shapefile contains polygons at several levels of an attribute. I want to plot these polygons for

[R-sig-Geo] best way to find contiguous chunks in a sparse space-time matrix

2017-04-11 Thread Waichler, Scott R
Hi, I have a space-time sparse grid layout where columns are locations and rows are time. In general, a location is missing data for some of its times. I need to subset the layout such that I can represent contiguous blocks of time, wherein for each block of time (rows), all of the selected

[R-sig-Geo] how to get weights for inverse distance weighting

2017-04-07 Thread Waichler, Scott R
Hello, I am trying to get inverse distance weights to estimate values on a regular grid from a set of data points, over a sequence of times. The locations of the data points don't vary with time, but their values will with each instance, and in general some are NA. I want to determine the

Re: [R-sig-Geo] Problem installing rgeos

2015-08-06 Thread Waichler, Scott R
a full test of its state. Following R CMD INSTALL, this should then work. Roger HTH Chris On Mon, Aug 3, 2015 at 10:34 AM, Waichler, Scott R scott.waich...@pnnl.gov wrote: Thanks for the suggestions, but I think what I need is help tracking down what the exact error

Re: [R-sig-Geo] Problem installing rgeos

2015-08-06 Thread Waichler, Scott R
, and run R CMD check on it, you will get a full test of its state. Following R CMD INSTALL, this should then work. Roger HTH Chris On Mon, Aug 3, 2015 at 10:34 AM, Waichler, Scott R scott.waich...@pnnl.gov wrote: Thanks for the suggestions, but I think what I need is help tracking

[R-sig-Geo] Problem installing rgeos

2015-07-28 Thread Waichler, Scott R
Hi, I'm trying to install rgeos and I don't know how to respond to the error message: configure: error: cannot run C compiled programs install.packages(rgeos, dependencies = Depends, repos = http://cran.fhcrc.org/;) trying URL 'http://cran.fhcrc.org/src/contrib/rgeos_0.3-11.tar.gz' Content

Re: [R-sig-Geo] image background for spplot of gridded data

2013-11-19 Thread Waichler, Scott R
I was told, among other things, to look for previous messages by Oscar that would explain how to do this. Here is Oscar's solution that I found to work for the case where you want to plot two rasters that have different extents. In my case, gmap is a satellite image obtained with RgoogleMaps,

[R-sig-Geo] image background for spplot of gridded data

2013-11-15 Thread Waichler, Scott R
Hi, I have a map tile from RgoogleMaps that I want to use as a background image for an spplot of gridded data, where the data will be plotted with an alpha value 1. I am unable to plot them together. I was going to set up an example with meuse.grid, but I couldn't get CRS() to work with

Re: [R-sig-Geo] working with a rotated grid of SpatialPolygons

2013-11-14 Thread Waichler, Scott R
How about just using high-level maptools? Following on directly from your code, the function elide can rotate the grid polygons: library(maptools) plot(elide(grdp, rotate = -45), add = TRUE) Yes, elide should work for my needs. Thank you for pointing it out! Is your data in a real

[R-sig-Geo] working with a rotated grid of SpatialPolygons

2013-11-13 Thread Waichler, Scott R
I have an environmental process model whose domain is a regular grid that is not oriented in the cardinal directions. I want to show an outline of the active part of the domain, as well as other subgroupings of cells within the domain, on a map figure. These will be blocky polygons oriented

Re: [R-sig-Geo] Use a PNG image as background for lattice plotting with rasterVis

2013-10-21 Thread Waichler, Scott R
Waichler -Original Message- From: Oscar Perpiñan [mailto:oscar.perpi...@gmail.com] Sent: Saturday, October 19, 2013 1:06 AM To: Waichler, Scott R Cc: r-sig-geo@r-project.org Subject: Re: [R-sig-Geo] Use a PNG image as background for lattice plotting with rasterVis Hello, Let's

[R-sig-Geo] Use a PNG image as background for lattice plotting with rasterVis

2013-10-18 Thread Waichler, Scott R
Hi, I've learned the basics of using RgoogleMaps and rasterVis, and now I want to combine their capabilities to use levelplot() with a georeferenced PNG image as a background in the lattice panel. And, I want to plot to the normal R graphics devices, not the web as package plotGoogleMaps

Re: [R-sig-Geo] function to convert to/from U.S. State Plane Coordinates?

2013-09-17 Thread Waichler, Scott R
Erin, Thanks, now I see that I was trying to use spTransform() on a map object instead of a Spatial one. I followed your logic and came up with the following to plot the state of Washington using state plane coordinates instead of the original geographic values. But I get an extraneous

Re: [R-sig-Geo] function to convert to/from U.S. State Plane Coordinates?

2013-09-16 Thread Waichler, Scott R
? From: r-sig-geo-boun...@r-project.org [r-sig-geo-boun...@r-project.org] on behalf of Waichler, Scott R [scott.waich...@pnnl.gov] Sent: Monday, September 16, 2013 4:33 PM To: r-sig-geo@r-project.org Subject: [R-sig-Geo] function to convert to/from U.S. State Plane Coordinates

[R-sig-Geo] function to convert to/from U.S. State Plane Coordinates?

2013-09-16 Thread Waichler, Scott R
Hi, Does an R function exist to convert to and from the State Plane Coordinates system used in the U.S.? My search with Rseek was fruitless. Scott Waichler Hydrology Group, Energy Environment Directorate Pacific Northwest National Laboratory Richland, WA, USA scott.waich...@pnnl.gov

Re: [R-sig-Geo] function to convert to/from U.S. State Plane Coordinates?

2013-09-16 Thread Waichler, Scott R
, Waichler, Scott R scott.waich...@pnnl.gov wrote: Erin, thanks, I'm getting closer.  Rgdal's projInfo() lists UTM, which I believe is closely related to the State Plane Coordinate system, but I still don't see a direct transformation to SPC in R.  --Scott Waichler Could you mean spTransform from

Re: [R-sig-Geo] function to convert to/from U.S. State Plane Coordinates?

2013-09-16 Thread Waichler, Scott R
-geo-bounces@r- project.org] On Behalf Of Waichler, Scott R Sent: Monday, September 16, 2013 6:49 PM To: Rob Robinson Cc: r-sig-geo@r-project.org Subject: Re: [R-sig-Geo] function to convert to/from U.S. State Plane Coordinates? Ok, now I know about the existence of EPSG codes, but apparently

Re: [R-sig-Geo] making an inset in a map figure

2013-09-13 Thread Waichler, Scott R
I found some problems with this code, but no worries, I got subplot() working and finished my figure. The resolution and sharpness of the figure is definitely superior to a version I had been trying to make with ArcMap. 1.R-3.0.1 didn't find subplot() in TeachingDemos, but it did in

[R-sig-Geo] making an inset in a map figure

2013-09-10 Thread Waichler, Scott R
Hi, Does anyone have an example of making a map figure with an inset, where a portion of the map in the large figure is shown at a larger scale in the inset? Scott Waichler Hydrology Group Pacific Northwest National Laboratory Richland, WA, USA ___

Re: [R-sig-Geo] how to use levelplot() with a geographic projection

2013-09-06 Thread Waichler, Scott R
For the archive, here is a complete example of using rasterVis to do the following: -- use levelplot() to draw multiple panels with maps being geographically correct -- create raster objects and supply values for some of the cells, selected from a matrix of longitude, latitude coordinates --

Re: [R-sig-Geo] how to use levelplot() with a geographic projection

2013-09-04 Thread Waichler, Scott R
Hi again, In working on this some more I discovered that I need to transpose a data matrix, then flip the raster object created from the matrix to get the proper rendering in a plot. Can someone explain why that is? My example has to do with the state of California, with 2918 1/8-degree

[R-sig-Geo] how to use levelplot() with a geographic projection

2013-08-26 Thread Waichler, Scott R
Hi, I am making lattice plots with levelplot(), where my data are on a regular grid in longitude and latitude, and I want to plot state boundaries as an overlay for reference. How can I do this so that the plot looks more geographically true, with states having the appearance we are used to?