[GRASS-user] error when importing soil grids data into GRASS

[Markus Neteler]

Hi Vero,

On Tue, Dec 7, 2021 at 6:12 PM Veronica Andreo <veroandreo at gmail.com> wrote:
>
> Dear all,
>
> I'm trying to import soil grids data for south america into grass gis by means of r.import and the vsicurl driver and I'm facing different errors according to the soilgrid version I use. See below:
>
> Fist attempt:
>
> gdalinfo /vsicurl/https://files.isric.org/soilgrids/latest/data/ocd/ocd_0-5cm_mean.vrt
> Driver: VRT/Virtual Raster
> Files: /vsicurl/https://files.isric.org/soilgrids/latest/data/ocd/ocd_0-5cm_mean.vrt
>        /vsicurl/https://files.isric.org/soilgrids/latest/data/ocd/ocd_0-5cm_mean.vrt.ovr
>        /vsicurl/https://files.isric.org/soilgrids/latest/data/ocd/./ocd_0-5cm_mean/tileSG-018-027/tileSG-018-027_2-4.tif
>        /vsicurl/https://files.isric.org/soilgrids/latest/data/ocd/./ocd_0-5cm_mean/tileSG-018-027/tileSG-018-027_1-1.tif
> ...
> Size is 159246, 58034
> Coordinate System is:
> PROJCRS["Interrupted_Goode_Homolosine",


This is the problem.

Here my workaround:

### ocd
# PROJCRS["Interrupted_Goode_Homolosine" ... :-(
# requires a trick - r.import still fails
eval `g.region -g`
gdalwarp -t_srs epsg:4326 -te $w $s $e $n ocd/ocd_0-5cm_mean.vrt
ocd_0_5cm_mean_epsg4326.tif -co COMPRESS=LZW

r.import input=ocd_0_5cm_mean_epsg4326.tif output=ocd_0_5cm_mean extent=region
r.colors ocd_0_5cm_mean color=bgyr

HTH,
Markus