Hi Vero,
in alternative you can use the WCS. A simple tool like wget is enough to query
the server:
$ wget
"https://maps.isric.org/mapserv?map=/map/nitrogen.map&SERVICE=WCS&VERSION=2.0.1&REQUEST=GetCoverage&COVERAGEID=nitrogen_5-15cm_Q0.5&FORMAT=GEOTIFF_INT16&SUBSET=X(-6535883,-6012080)&SUBSET=Y(-3921059,-3304244)&SUBSETTINGCRS=http://www.opengis.net/def/crs/EPSG/0/152160&OUTPUTCRS=http://www.opengis.net/def/crs/EPSG/0/152160";
-O Uruguai_nit_5-15.tiff
--2021-12-13 12:42:56--
https://maps.isric.org/mapserv?map=/map/nitrogen.map&SERVICE=WCS&VERSION=2.0.1&REQUEST=GetCoverage&COVERAGEID=nitrogen_5-15cm_Q0.5&FORMAT=GEOTIFF_INT16&SUBSET=X(-6535883,-6012080)&SUBSET=Y(-3921059,-3304244)&SUBSETTINGCRS=http://www.opengis.net/def/crs/EPSG/0/152160&OUTPUTCRS=http://www.opengis.net/def/crs/EPSG/0/152160
Resolving maps.isric.org (maps.isric.org)... 137.224.206.3
Connecting to maps.isric.org (maps.isric.org)|137.224.206.3|:443... connected.
HTTP request sent, awaiting response... 200 OK
Length: unspecified [image/tiff]
Saving to: ‘Uruguai_nit_5-15.tiff’
Uruguai_nit_5-15.tiff [ <=> ] 3,16M 5,62MB/s in 0,6s
2021-12-13 12:43:03 (5,62 MB/s) - ‘Uruguai_nit_5-15.tiff’ saved [3314780]
Then apply r.in.gdal directly:
> r.in.gdal input=Uruguai_nit_5-15.tiff output=uruguai_nit_5_15 -o
Over-riding projection check
Importing raster map <uruguai_nit_5_15>...
100%
You can request the result in a different CRS directly from the WCS. The
GetCapabilities request informs you which are available. Although I would
advise to work either on the original CRS or with another composed by an
equal-area projection.
Best.
--
Luís
Sent with [ProtonMail](https://protonmail.com/) Secure Email.
‐‐‐‐‐‐‐ Original Message ‐‐‐‐‐‐‐
On Tuesday, December 7th, 2021 at 6:11 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",
> BASEGEOGCRS["GCS_WGS_1984 ellipse",
> DATUM["World Geodetic System 1984",
> ELLIPSOID["WGS 84",6378137,298.257223563,
> LENGTHUNIT["metre",1]],
> ID["EPSG",6326]],
> PRIMEM["Greenwich",0,
> ANGLEUNIT["Degree",0.0174532925199433]]],
> CONVERSION["unnamed",
> METHOD["Interrupted Goode Homolosine"]],
> CS[Cartesian,2],
> AXIS["(E)",east,
> ORDER[1],
> LENGTHUNIT["metre",1,
> ID["EPSG",9001]]],
> AXIS["(N)",north,
> ORDER[2],
> LENGTHUNIT["metre",1,
> ID["EPSG",9001]]]]
>
> r.import
> input="/vsicurl/https://files.isric.org/soilgrids/latest/data/ocd/ocd_0-5cm_mean.vrt";
> output=ocd_0_5cm_mean extent=region
> ...
> Selected PROJ pipeline:
> +proj=pipeline +step +proj=unitconvert +xy_in=deg +xy_out=rad +step
> +proj=igh +lon_0=0 +x_0=0 +y_0=0 +ellps=WGS84
> ************************
> Reproyectado primitivas...
> Construyendo topología para el mapa vectorial
> <tmp_r_import_region_localhost_localdomain_918305@PERMANENT>...
> Registrando primitivas...
> Construyendo áreas.
> 100%
> Adjuntando islas...
> 100%
> Adjuntando centróides...
> 100%
> Importando mapa ráster <ocd_0_5cm_mean>...
> ERROR 4:
> `/vsicurl/https://files.isric.org/soilgrids/latest/data/ocd/./ocd_0-5cm_mean/tileSG-022-030/tileSG-022-030_2-3.tif'
> not recognized as a supported file format.
> ERROR: Unable to read row 11740
> ERROR: Unable to import GDAL dataset
> </vsicurl/https://files.isric.org/soilgrids/latest/data/ocd/ocd_0-5cm_mean.vrt>
>
> Second attempt, newer dataset soilgrids 2.0 pre-release
>
> r.import
> input="/vsicurl/https://storage.googleapis.com/isric-share-soilgrids/pre-release/wv1500/wv1500_0-5cm_mean.tif";
> output=wv1500_0_5cm_mean extent=region
> Warning 1:
> /vsicurl/https://storage.googleapis.com/isric-share-soilgrids/pre-release/wv1500/wv1500_0-5cm_mean.tif:
> This file used to have optimizations in its layout, but those have been, at
> least partly, invalidated by later changes
> ADVERTENCIA: Actualizando el sistema de referencia espacial con la
> definición proj4 incrustada.
> ...
> Selected PROJ pipeline:
> +proj=pipeline +step +proj=unitconvert +xy_in=deg +xy_out=rad +step
> +proj=igh +lon_0=0 +x_0=0 +y_0=0 +ellps=WGS84
> ************************
> Reproyectado primitivas...
> Construyendo topología para el mapa vectorial
> <tmp_r_import_region_localhost_localdomain_918176@PERMANENT>...
> Registrando primitivas...
> Construyendo áreas.
> 100%
> Adjuntando islas...
> 100%
> Adjuntando centróides...
> 100%
> Warning 1:
> /vsicurl/https://storage.googleapis.com/isric-share-soilgrids/pre-release/wv1500/wv1500_0-5cm_mean.tif:
> This file used to have optimizations in its layout, but those have been, at
> least partly, invalidated by later changes
> ADVERTENCIA: Actualizando el sistema de referencia espacial con la
> definición proj4 incrustada.
> Importando mapa ráster <wv1500_0_5cm_mean>...
> 100%
> Removiendo vector <tmp_r_import_region_localhost_localdomain_918176>
> Removiendo vector <tmp_r_import_region_localhost_localdomain_918176>
> Estimated target resolution for input band <wv1500_0_5cm_mean>:
> 0.0022003864773300253
> Reprojecting <wv1500_0_5cm_mean>...
> ADVERTENCIA: proj_create() failed for '(null)'
> ERROR: Unable to initialize coordinate transformation
> ERROR: Unable to to reproject raster <wv1500_0_5cm_mean>
>
> It imports the raster but fails at re-projection... CRS is unknown in the tif
>
> gdalinfo
> /vsicurl/https://storage.googleapis.com/isric-share-soilgrids/pre-release/wv1500/wv1500_0-5cm_mean.tif
> Warning 1:
> /vsicurl/https://storage.googleapis.com/isric-share-soilgrids/pre-release/wv1500/wv1500_0-5cm_mean.tif:
> This file used to have optimizations in its layout, but those have been, at
> least partly, invalidated by later changes
> Driver: GTiff/GeoTIFF
> Files:
> /vsicurl/https://storage.googleapis.com/isric-share-soilgrids/pre-release/wv1500/wv1500_0-5cm_mean.tif
> Size is 159246, 58034
> Coordinate System is:
> PROJCRS["unknown",
> BASEGEOGCRS["unknown",
> DATUM["World Geodetic System 1984",
> ELLIPSOID["WGS 84",6378137,298.257223563,
> LENGTHUNIT["metre",1]],
> ID["EPSG",6326]],
> PRIMEM["Greenwich",0,
> ANGLEUNIT["degree",0.0174532925199433],
> ID["EPSG",8901]]],
> CONVERSION["unknown",
> METHOD["Interrupted Goode Homolosine"],
> PARAMETER["Longitude of natural origin",0,
> ANGLEUNIT["degree",0.0174532925199433],
> ID["EPSG",8802]],
> ...
>
> Any suggestion? The second attempt seems a problem of the data provider not
> defining the CRS properly. But what about the first case, where the CRS is
> defined, any ideas on why the projection is not recognized? Any experience
> importing these data?
>
> Thanks much in advance!
> Vero
>
> My system info:
> Versión de GRASS: 8.0.dev
> Fecha de compilación: 2021-10-13
> GDAL: 3.2.2
> PROJ: 7.2.1
> GEOS: 3.9.0
> SQLite: 3.34.1
> Python: 3.9.7
> wxPython: 4.0.7
> Plataforma: Linux-5.14.9-200.fc34.x86_64-x86_64-with-glibc2.33
_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user
