need (maybe) imports from sp so that raster() works with CRS, not sure
indexing methods for coercion to raster