Hi Vero,

On Tue, Dec 7, 2021 at 6:12 PM Veronica Andreo <veroand...@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
_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Reply via email to