I _think_ raster may be solving a different problem. I am looking for a
voxel-type raster. Let me know if I am wrong.
plot3d(raster(kriged_sp))
Error in as.double(x) :
cannot coerce type 'S4' to vector of type 'double'
Let me describe the SpatialGridDataFrame a bit better. The coordinates
are position in volumetric space (x,y,depth_cm). The value at each
position is the predicted value from IDW (var1.pred). So the structure is
summary(kriged_sp)
Object of class SpatialGridDataFrame
Coordinates:
min max
x 740317.2 740747.2
y 181326.1 181646.1
depth_cm -90.0 0.0
Is projected: TRUE
proj4string :
[+init=epsg:32119 +proj=lcc +lat_1=36.16666666666666
+lat_2=34.33333333333334 +lat_0=33.75 +lon_0=-79 +x_0=609601.22 +y_0=0
+ellps=GRS80 +datum=NAD83 +units=m +no_defs +towgs84=0,0,0]
Number of points: 2
Grid attributes:
cellcentre.offset cellsize cells.dim
x 740322.2 10 43
y 181331.1 10 32
depth_cm -87.5 5 18
Data attributes:
var1.pred var1.var
Min. : 3.66 Min. : NA
1st Qu.:19.92 1st Qu.: NA
Median :22.97 Median : NA
Mean :22.08 Mean : NaN
3rd Qu.:24.62 3rd Qu.: NA
Max. :44.26 Max. : NA
NA's :24768
When I create a raster from this object, I get something that looks
two-dimensional:
class : RasterLayer
filename :
nrow : 32
ncol : 43
ncells : 1376
min value :
max value :
projection : +init=epsg:32119 +proj=lcc +lat_1=36.16666666666666 +lat_2=34.33333333333334 +lat_0=33.75 +lon_0=-79 +x_0=609601.22 +y_0=0 +ellps=GRS80 +datum=NAD83 +units=m +no_defs +towgs84=0,0,0
xmin : 740317.2
xmax : 740747.2
ymin : 181326.1
ymax : 181646.1
xres : 10
yres : 10
Robert J. Hijmans wrote:
Mark,
If spgdf is a SpatialGridDataFrame object
# install.packages("raster", repos="http://R-Forge.R-project.org";)
library(raster)
plot3D(raster(spgdf))
Or you can have a look at the code:
getMethod("plot3D", 'RasterLayer')
Robert
On Wed, Dec 9, 2009 at 5:10 AM, Mark Connolly <mark_conno...@acm.org> wrote:
I have a 3d SpatialGridDataFrame I created using krige (IDW). I am using
rgl facilities to create an interactive plot of the points, but I'd really
like to have something more volumetric than points. My copy of the truly
wonderful Applied Spatial Data Analysis with R (shameless pander) is pretty
mute on 3d. Any hints on where I might find something on transforming the
points into a mesh? Something I need to code for myself?
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
