Hi,
I'm working on interpolating spatial point data to a grid and as I was
troubleshooting I found that the gridded locations I am predicting to are
not what I expected. In particularly, when I extracted the locations from
the raster grid to get the cell coordinates in which the locations lay, the
Dominik,
The last part of your example does not run. If I understand you well,
this is what you are after (but it all looks good to me):
library(raster)
# raster data with long and lat values
r - raster(xmn=-112.25, xmx=-104.125, ymn=33, ymx=43.75, res=1/240)
lat - init(r, v='y')
lon - init(r,