Alison,
there's two ways for this:
1. provide krige() with regression coordinates, by
ObsSat.rk <- krige(Obs_rain~sat_rain, ObsSat_UTM, SatDat_UTM,
ObsSat.rvgm, beta = c(a,b))
with a and b your "manual" intercept and slope, or
2. by kriging manually computed residuals and adding them to the
Hello list,
I am trying to use universal kriging (gstat, R) to create a gridded
daily rainfall surface. I have rain gauge data and satellite grid data
and have applied the following script successfully. I would now like to
manually define the regression model that defines the relationship
between
Roger Bivand wrote:
> On Thu, 13 Nov 2008, Patrick Giraudoux wrote:
>
>> Roger Bivand a écrit :
>>> On Thu, 13 Nov 2008, Patrick Giraudoux wrote:
>>>
Dear listers,
I am using GRASS6 embedded within QGIS under Windows for teaching
purpose. This prevents reading GRASS objects fro
On Thu, 13 Nov 2008, Patrick Giraudoux wrote:
Roger Bivand a écrit :
On Thu, 13 Nov 2008, Patrick Giraudoux wrote:
Dear listers,
I am using GRASS6 embedded within QGIS under Windows for teaching purpose.
This prevents reading GRASS objects from R within GRASS (R cannot be
launched within G
Roger Bivand a écrit :
On Thu, 13 Nov 2008, Patrick Giraudoux wrote:
Dear listers,
I am using GRASS6 embedded within QGIS under Windows for teaching
purpose. This prevents reading GRASS objects from R within GRASS (R
cannot be launched within GRASS in this configuration).
Patrick,
Could y
Tomislav Hengl wrote:
Alessandro,
You need to separate the put the list of input grids as a character vector e.g.:
GRIDS=c("DSM1.sgrd,DSM2.sgrd,DSM3.sgrd")
...not when calling the rsaga.geoprocessor, if I recall correctly. The
high-level wrapper functions would accept a vector, but the
rsa
Alessandro,
You need to separate the put the list of input grids as a character vector e.g.:
GRIDS=c("DSM1.sgrd,DSM2.sgrd,DSM3.sgrd")
See also:
http://spatial-analyst.net/wiki/index.php?title=Software#SAGA_GIS
Maybe you should also study this posting guide:
http://spatial-analyst.net/wiki/ind
On Thu, 13 Nov 2008, Patrick Giraudoux wrote:
Dear listers,
I am using GRASS6 embedded within QGIS under Windows for teaching purpose.
This prevents reading GRASS objects from R within GRASS (R cannot be launched
within GRASS in this configuration).
Patrick,
Could you explain why not? I ca
Hi All
when I use this function with R and SAGA:
rsaga.get.usage("grid_tools", 3)
rsaga.geoprocessor(lib="grid_tools", module=3,
param=list(GRIDS="DSM1.sgrd,DSM2.sgrd,DSM3.sgrd",GRID_TARGET="DSM.sgrd",TYPE=7, INTERPOL=2,
OVERLAP=0))
I have this problem
SAGA CMD 2.0.3
library path: C:\Wi
