Настройте обратное предсказание расстояния в растровых классах в R

Я хочу изменить следующий код, чтобы для всех значений за пределами определенного максимального расстояния я мог использовать средние значения точек, которые попадают в одну и ту же категорию порядка почвы. Есть рекомендации?

# load packages and data
library(gstat)
library(raster)
data(meuse)
data(meuse.grid)
################## 
meuse <- meuse[110:155,]
meuse <- na.omit(meuse)
coordinates(meuse) <- ~x+y
coordinates(meuse.grid) <- ~x+y
gridded(meuse.grid) <- TRUE
spplot(meuse.grid[4], main="soil class")
# inverse distance prediction for maximum distance of 2km.
ca.idw = idw(cadmium ~ 1, meuse, meuse.grid, omax =5, maxdist = 2000)