I want to do kriging with external drift(universal kriging) in R. The covariate I use is in WGS84 refrence system. Because of nature of kriging which calculates distance in planar surface I have to do transformation from WGS84 to ETRS89. what I want out of kriging is a raster map. Because of the reference system transformation the prediction points are not regular grids and what I get as a result is a vector map.
I used 'gstat' package in and spplot()
function to create plots. But When grids are regular (wich I can not define because of nature of the covariate data) with equal distance I have a raster map.