Hi Veronica,
thanks for pointing this out. The projection should be defined in the VRTs
of Soilgrids 2.0 and normally it should be recognised in GRASS-GIS. I have
been using the same files you mentioned. We will look into the issue and
let you know.

Thanks
Laura

On Wed, 8 Dec 2021 at 01:31, Veronica Andreo <veroand...@gmail.com> wrote:

> Hi Markus,
>
> thanks for your answer :)
>
> El mar., 7 dic. 2021 14:43, Markus Neteler <nete...@osgeo.org> escribió:
>
>> 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.
>>
>
> Any idea why it's not recognized? I see it is in PROJ and as you show
> gdalwarp works. Is it something GRASS related?
>
> 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
>>
>
> Anyway replacing ocd/ocd_0-5cm_mean.vrt by /vsicurl/
> https://files.isric.org/soilgrids/latest/data/ocd/ocd_0-5cm_mean.vrt in
> your workaround does the trick, takes a while, but does it. Thanks a LOT
> for the hint!!
>
> @Luí­s Moreira de Sousa <luis.de.so...@protonmail.ch> would it be
> possible to provide the version 2.0 with the projection defined? Or
> additionally in the easy to use epsg 4326 as an alternative for lazy people
> like me? :)
>
> Vero
> _______________________________________________
> grass-user mailing list
> grass-user@lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/grass-user
>
_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Reply via email to