Roger,
This is what the summary() command gives me before/after I use
fullgrid(dem) <- FALSE
dem<-readFLOAT6sp("ohrfc.dem")
> summary(dem)
Object of class SpatialGridDataFrame
Coordinates:
min max
coords.x1 -351000 675000
coords.x2 -486000 486000
Is projected: TRUE
proj4string : [+proj=lcc +la
On Fri, 28 Apr 2006, Thomas Adams wrote:
> Roger,
>
> Your suggestion:
>
> fullgrid(dem) <- FALSE
>
> did turn dem into class type SpatialGridDataFrame, but when I tried:
>
> z <- predict(UK_fit,newdata=dem)
>
> I got an error:
>
> Error in model.frame(... :
> invalid variable type.
>
> I t
Roger,
Your suggestion:
fullgrid(dem) <- FALSE
did turn dem into class type SpatialGridDataFrame, but when I tried:
z <- predict(UK_fit,newdata=dem)
I got an error:
Error in model.frame(... :
invalid variable type.
I think I should restate the problem:
I have a file 'temps' which has class
On Fri, 28 Apr 2006, Thomas Adams wrote:
> Roger,
>
> This got me further along, but I am encountering a problem with:
>
> z <- predict(UK_fit, newdata=BMcD_SPx)
>
> The gstat step works for me, where I have:
>
> UK_fit<-gstat(formula=temps$temp~dem,data=temps,model=efitted)
>
> temps has cla
Roger,
This got me further along, but I am encountering a problem with:
z <- predict(UK_fit, newdata=BMcD_SPx)
The gstat step works for me, where I have:
UK_fit<-gstat(formula=temps$temp~dem,data=temps,model=efitted)
temps has class SpatialPointsDataFrame:
coordinates cat x
age -
From: Roger Bivand <[EMAIL PROTECTED]>
Date: Thursday, April 27, 2006 4:21 pm
Subject: Re: [R-sig-Geo] Example of universal kriging with R/gstat in GRASS
needed
> On Thu, 27 Apr 2006, Thomas Adams wrote:
>
> > List:
> >
> > I can not seem to work out the
On Thu, 27 Apr 2006, Thomas Adams wrote:
> List:
>
> I can not seem to work out the syntax for using R/gstat within a GRASS
> 6.1 session to do universal kriging. I have a DEM (elevation data on a
> grid) and point data for temperature; theoretically, the temperatures
> should relate to elevat
List:
I can not seem to work out the syntax for using R/gstat within a GRASS
6.1 session to do universal kriging. I have a DEM (elevation data on a
grid) and point data for temperature; theoretically, the temperatures
should relate to elevation. So, I am trying to spatially interpolate the
tem