Re: [R-sig-Geo] accumulate gridcells to a specific value

2008-07-31 Thread Alexander.Herr
Hi Barry, > Can you explain a bit more about your problem? > * What are the inputs? Just a grid of values or are there some starting points? Inputs is a floating grid. Starting point is 1) the highest values in the grid. Once first accumulation achieved, exclude this area and 2) find highest valu

Re: [R-sig-Geo] ncdf file with non-equally spaced grid

2008-07-31 Thread Matt Oliver
Sorry, dlat is what changes an example file can be found at http://www.ocean.udel.edu/CMS/moliver/20071003.276.0237.n17.nc use require(ncdf) f <- open.ncdf("20071003.276.0237.n17.nc" ) lon <- get.var.ncdf(f, "lon") lat <- get.v

Re: [R-sig-Geo] ncdf file with non-equally spaced grid

2008-07-31 Thread Michael Sumner
You can use gdalwarp in FWTools to re-grid the data with control points (perhaps with with some effort), but as you say that will interpolate to a regular grid. Is it definitely dlon that increases? Can you post a link to an example file? Some netCDF files I've seen are actually using Mercat

Re: [R-sig-Geo] specifying different colors for vector map as background

2008-07-31 Thread Roger Bivand
On Thu, 31 Jul 2008, Jaime Carrera wrote: Hi, I'm trying to create a spplot with a vector map on the background. It works fine if I don't want the polygon map to be filled, but what I do want to is to have the polygons filled with different colors. The following command plots the polygon map

[R-sig-Geo] specifying different colors for vector map as background

2008-07-31 Thread Jaime Carrera
Hi, I'm trying to create a spplot with a vector map on the background. It works fine if I don't want the polygon map to be filled, but what I do want to is to have the polygons filled with different colors. The following command plots the polygon map fine, filling each polygon with its respect

[R-sig-Geo] ncdf file with non-equally spaced grid

2008-07-31 Thread Matt Oliver
Hi All, I have CF compliant ncdf files of sea surface temperature. The lon, lat grid is known but not equally spaced (ie dlon increases). Is there a clean way to get this matrix ported over to a geo-tiff or some other file that Arc will read? I have talked with some arc users that tell me that the

Re: [R-sig-Geo] accumulate gridcells to a specific value

2008-07-31 Thread Barry Rowlingson
2008/7/31 <[EMAIL PROTECTED]>: > I am trying to accumulate gridcells to a specific value. This means I > need to extent the nearest neighbours (or search radius) to a size so > that the sum of gridcell values reaches the desired value (or stops if > sum value of gridcells/radius goes beyond a spe

Re: [R-sig-Geo] reading WKB geometry

2008-07-31 Thread Edzer Pebesma
Jaime, I've been successful in reading points from a PostGIS data base through OGR; see http://wiki.intamap.org/index.php/PostGIS -- Edzer Jaime Carrera wrote: Hi list, I've been looking around trying to solve this problem; unfortunately I haven't found a way to solve it. I have a database