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 the gauge and satellite at the collocated points so I can apply
the same model to a series of days. Is it possible to manually define
the regression equation and have it applied within the krige function? 

 

Datasets

ObsSat_UTM = collocated gauge and satellite data.

SatDat_UTM = Satellite grid. 

 

# Define the regression model to define relationships between dependent
(gauge) and independent (sat) variables at collocated points. 

ObsSat.lm <-lm(formula = Obs_rain~sat_rain, ObsSat_UTM)

 

#Variogram of residuals of regression between obs and sat data a gauge
sites

ObsSat.rev <- variogram(Obs_rain~sat_rain, ObsSat_UTM)

ObsSat.rvgm <- fit.variogram(ObsSat.rev, vgm(20,"Exp", 50,0))

 

#predict onto the Sat grid

ObsSat.rk <- krige(Obs_rain~sat_rain, ObsSat_UTM, SatDat_UTM,
ObsSat.rvgm)

 

If this question has already been answered I apologise. I have had
trouble finding any examples of how to include this into my script.

Thanks very much 

Alison Oke

 

Hydrologist

Water Division

Bureau of Meteorology

GPO Box 1289

Melbourne, VIC 3001

 

Phone:  +61 (0) 3 8638 8212

Fax:      +61 (0) 3 8638 8200

E-mail:   a.oke@bom.gov.au

 


	[[alternative HTML version deleted]]


