I need the global SOC stocks map at 0-30cm from the new SoilGrids 2.0: https://files.isric.org/soilgrids/latest/data/ocs/
I have no previous experience with VRT files. I have simply downloaded "ocs_0-30cm_mean.vrt" and "ocs_0-30cm_mean.vrt.ovr". Honestly, I don't know if I need one, the other, or both.
In R, I load it with the raster package.
soc <- raster("ocs_0-30cm_mean.vrt.ovr") # t/ha soc class : RasterLayer dimensions : 14509, 39812, 577632308 (nrow, ncol, ncell) resolution : 1, 1 (x, y) extent : 0, 39812, 0, 14509 (xmin, xmax, ymin, ymax) crs : NA source : ocs_0-30cm_mean.vrt.ovr names : ocs_0.30cm_mean.vrt values : -32768, 32767 (min, max)
However, the resolution appears to be 1 degree, when it should be 250m. The extent makes no sense (I guess it should be c(-180, 180, -90, 90), and there is no CRS.
I can plot (soc), though I don't know how to change the resolution to e.g. 0.25.
I have also tried the same with "ocs_0-30cm_mean.vrt" (instead of "ocs_0-30cm_mean.vrt.ovr"), but the file does not
plot(soc) Error in rgdal::getRasterData(con, offset = offs, region.dim = c(1, nc), : Failure during raster IO
How do I make this raster work, with 0.25 resolution and appropriate extent?