geoWeight {remote}

Geographic weighting

Description

The function performs geographic weighting of non-projected long/lat data. By default it uses the cosine of latitude to compensate for the area distortion, though the user can supply other functions via f.

Usage

geoWeight(x, f = function(x) cos(x), ...)

Arguments

x

a Raster* object

f

a function to be used to the weighting. Defaults to cos(x)

...

additional arguments to be passed to f

Value

a weighted Raster* object

Examples

data(vdendool)

wgtd <- geoWeight(vdendool)

opar <- par(mfrow = c(1,2))
plot(vdendool[[1]], main = "original")
plot(wgtd[[1]], main = "weighted")
par(opar)