I use ncdf V1.6 on Windows XP Sp2, I hava attached the sample file (tab
delim) which contains 2 columns for coordinates long/lat and 10 columns for
weekly averages of radon. Sp would be the best solution to produce map I
guess, but for time series analysis I think netcdf would perform better.
Tha
e.g.
# x is a gpc.poly object
as.Polygons.gpc.poly <- function(x, ID) {
thisPolys <- lapply(get.pts(x), function(p) {
Polygon(rbind(as.matrix(cbind(p$x,p$y)),c(p$x[1],p$y[1])),
hole=p$hole)
})
Polygons(thisPolys, ID)
}
This works exc
Hmm
I'm trying to work with gpclib, converting to and from the sp classes.
I don't know what is/was in the package spgpc, but it is no longer
available:
> install.packages("spgpc", repos=rSpatial)
Warning: unable to access index for repository
http://r-spatial.sourceforge.net/R/bin/windows/contrib
Hi, can you provide an example file?
I'd recommend avoiding NetCDF for various reasons, why don't you simply add
each value column from the separate files to a new column in a
SpatialGridDataFrame?
That would be the recommend sp way to proceed here. If you really need to
create a NetCDF file, ple
Hi!
I have a table with 3 columns, 2 for long/lat coordinates and 1 for values
(radon concentration). I have data for every week of the year
2006,distributed in 52 tables/files. I want to create a NetCDF file, which
is much easier to handle than extracting the values from 52 tables. I tried
the fo
On Wed, 8 Aug 2007, Felix Andrews wrote:
Thanks Michael
The extension is "flt" (float type). This resulted in a sensible image:
foo <-
readBin("G:/Projects/Tuross/surface_files/rainfall_surface/year2000/rainGrid_20001.flt",
"double", n=200*200, size=4, endian="little")
foo[foo==-] <- NA
r
Are there R routines that do not request that x,y,z data be increasing
before plotting? Trying to generate basic surfaces from
three-dimensional real data. even if I put an ascending sort on one
axis the data along projected remaining coordinates may not be
increasing, or monotonic. What routines
