<div>Hi Vero,<br></div><div><br></div><div>in alternative you can use the WCS. A simple tool like wget is enough to query the server:<br></div><div><br></div><div>$ wget "<a href="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" rel="noopener noreferrer" target="_blank">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</a>" -O Uruguai_nit_5-15.tiff<br></div><div><br></div><div>--2021-12-13 12:42:56-- <a href="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" rel="noopener noreferrer" target="_blank">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</a><br></div><div>Resolving maps.isric.org (maps.isric.org)... 137.224.206.3<br></div><div>Connecting to maps.isric.org (maps.isric.org)|137.224.206.3|:443... connected.<br></div><div>HTTP request sent, awaiting response... 200 OK<br></div><div>Length: unspecified [image/tiff]<br></div><div>Saving to: ‘Uruguai_nit_5-15.tiff’<br></div><div><br></div><div>Uruguai_nit_5-15.tiff [ <=> ] 3,16M 5,62MB/s in 0,6s<br></div><div><br></div><div>2021-12-13 12:43:03 (5,62 MB/s) - ‘Uruguai_nit_5-15.tiff’ saved [3314780]<br></div><div><br></div><div>Then apply r.in.gdal directly:<br></div><div><br></div><div>> r.in.gdal input=Uruguai_nit_5-15.tiff output=uruguai_nit_5_15 -o<br></div><div>Over-riding projection check<br></div><div>Importing raster map <uruguai_nit_5_15>...<br></div><div>100%<br></div><div><br></div><div>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. <br></div><div class="protonmail_signature_block"><div class="protonmail_signature_block-user"><div><br></div><div>Best.<br></div><div><br></div></div></div><div><br></div><div class="protonmail_signature_block"><div class="protonmail_signature_block-user"><div><br></div><div>-- <br></div><div>Luís<br></div><div><br></div><div>Sent with <a rel="noopener noreferrer" href="https://protonmail.com/" target="_blank">ProtonMail</a> Secure Email.<br></div><div><br></div></div></div><div><br></div><div class="protonmail_quote">
‐‐‐‐‐‐‐ Original Message ‐‐‐‐‐‐‐<br>
On Tuesday, December 7th, 2021 at 6:11 PM, Veronica Andreo <veroandreo@gmail.com> wrote:<br>
<blockquote class="protonmail_quote" type="cite">
<div dir="ltr"><div>Dear all, <br></div><div><br></div><div>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:<br></div><div><br></div><div>Fist attempt: <br></div><div><br></div><div>gdalinfo /vsicurl/<a href="https://files.isric.org/soilgrids/latest/data/ocd/ocd_0-5cm_mean.vrt" rel="noreferrer nofollow noopener" target="_blank">https://files.isric.org/soilgrids/latest/data/ocd/ocd_0-5cm_mean.vrt</a><br>Driver: VRT/Virtual Raster<br>Files: /vsicurl/<a href="https://files.isric.org/soilgrids/latest/data/ocd/ocd_0-5cm_mean.vrt" rel="noreferrer nofollow noopener" target="_blank">https://files.isric.org/soilgrids/latest/data/ocd/ocd_0-5cm_mean.vrt</a><br> /vsicurl/<a href="https://files.isric.org/soilgrids/latest/data/ocd/ocd_0-5cm_mean.vrt.ovr" rel="noreferrer nofollow noopener" target="_blank">https://files.isric.org/soilgrids/latest/data/ocd/ocd_0-5cm_mean.vrt.ovr</a><br> /vsicurl/<a href="https://files.isric.org/soilgrids/latest/data/ocd/./ocd_0-5cm_mean/tileSG-018-027/tileSG-018-027_2-4.tif" rel="noreferrer nofollow noopener" target="_blank">https://files.isric.org/soilgrids/latest/data/ocd/./ocd_0-5cm_mean/tileSG-018-027/tileSG-018-027_2-4.tif</a><br> /vsicurl/<a href="https://files.isric.org/soilgrids/latest/data/ocd/./ocd_0-5cm_mean/tileSG-018-027/tileSG-018-027_1-1.tif" rel="noreferrer nofollow noopener" target="_blank">https://files.isric.org/soilgrids/latest/data/ocd/./ocd_0-5cm_mean/tileSG-018-027/tileSG-018-027_1-1.tif</a><br>...<br>Size is 159246, 58034<br>Coordinate System is:<br>PROJCRS["Interrupted_Goode_Homolosine", <br> BASEGEOGCRS["GCS_WGS_1984 ellipse",<br> DATUM["World Geodetic System 1984",<br> ELLIPSOID["WGS 84",6378137,298.257223563,<br> LENGTHUNIT["metre",1]],<br> ID["EPSG",6326]],<br> PRIMEM["Greenwich",0,<br> ANGLEUNIT["Degree",0.0174532925199433]]],<br> CONVERSION["unnamed",<br> METHOD["Interrupted Goode Homolosine"]],<br> CS[Cartesian,2],<br> AXIS["(E)",east,<br> ORDER[1],<br> LENGTHUNIT["metre",1,<br> ID["EPSG",9001]]],<br> AXIS["(N)",north,<br> ORDER[2],<br> LENGTHUNIT["metre",1,<br> ID["EPSG",9001]]]]</div><div><br></div><div>r.import input="/vsicurl/<a href="https://files.isric.org/soilgrids/latest/data/ocd/ocd_0-5cm_mean.vrt" rel="noreferrer nofollow noopener" target="_blank">https://files.isric.org/soilgrids/latest/data/ocd/ocd_0-5cm_mean.vrt</a>" output=ocd_0_5cm_mean extent=region</div><div>...<br>Selected PROJ pipeline:<br>+proj=pipeline +step +proj=unitconvert +xy_in=deg +xy_out=rad +step<br>+proj=igh +lon_0=0 +x_0=0 +y_0=0 +ellps=WGS84<br>************************<br>Reproyectado primitivas...<br>Construyendo topología para el mapa vectorial<br><tmp_r_import_region_localhost_localdomain_918305@PERMANENT>...<br>Registrando primitivas...<br>Construyendo áreas.<br> 100%<br>Adjuntando islas...<br> 100%<br>Adjuntando centróides...<br> 100%<br>Importando mapa ráster <ocd_0_5cm_mean>...<br>ERROR 4: `/vsicurl/<a href="https://files.isric.org/soilgrids/latest/data/ocd/./ocd_0-5cm_mean/tileSG-022-030/tileSG-022-030_2-3.tif" rel="noreferrer nofollow noopener" target="_blank">https://files.isric.org/soilgrids/latest/data/ocd/./ocd_0-5cm_mean/tileSG-022-030/tileSG-022-030_2-3.tif</a>' not recognized as a supported file format.<br>ERROR: Unable to read row 11740<br>ERROR: Unable to import GDAL dataset<br> </vsicurl/<a href="https://files.isric.org/soilgrids/latest/data/ocd/ocd_0-5cm_mean.vrt" rel="noreferrer nofollow noopener" target="_blank">https://files.isric.org/soilgrids/latest/data/ocd/ocd_0-5cm_mean.vrt</a>></div><div><br></div><div>Second attempt, newer dataset soilgrids 2.0 pre-release<br></div><div><br></div><div>r.import input="/vsicurl/<a href="https://storage.googleapis.com/isric-share-soilgrids/pre-release/wv1500/wv1500_0-5cm_mean.tif" rel="noreferrer nofollow noopener" target="_blank">https://storage.googleapis.com/isric-share-soilgrids/pre-release/wv1500/wv1500_0-5cm_mean.tif</a>" output=wv1500_0_5cm_mean extent=region<br>Warning 1: /vsicurl/<a href="https://storage.googleapis.com/isric-share-soilgrids/pre-release/wv1500/wv1500_0-5cm_mean.tif" rel="noreferrer nofollow noopener" target="_blank">https://storage.googleapis.com/isric-share-soilgrids/pre-release/wv1500/wv1500_0-5cm_mean.tif</a>: This file used to have optimizations in its layout, but those have been, at least partly, invalidated by later changes<br>ADVERTENCIA: Actualizando el sistema de referencia espacial con la<br> definición proj4 incrustada.<br>...<br>Selected PROJ pipeline:<br>+proj=pipeline +step +proj=unitconvert +xy_in=deg +xy_out=rad +step<br>+proj=igh +lon_0=0 +x_0=0 +y_0=0 +ellps=WGS84<br>************************<br>Reproyectado primitivas...<br>Construyendo topología para el mapa vectorial<br><tmp_r_import_region_localhost_localdomain_918176@PERMANENT>...<br>Registrando primitivas...<br>Construyendo áreas.<br> 100%<br>Adjuntando islas...<br> 100%<br>Adjuntando centróides...<br> 100%<br>Warning 1: /vsicurl/<a href="https://storage.googleapis.com/isric-share-soilgrids/pre-release/wv1500/wv1500_0-5cm_mean.tif" rel="noreferrer nofollow noopener" target="_blank">https://storage.googleapis.com/isric-share-soilgrids/pre-release/wv1500/wv1500_0-5cm_mean.tif</a>: This file used to have optimizations in its layout, but those have been, at least partly, invalidated by later changes<br>ADVERTENCIA: Actualizando el sistema de referencia espacial con la<br> definición proj4 incrustada.<br>Importando mapa ráster <wv1500_0_5cm_mean>...<br> 100%<br>Removiendo vector <tmp_r_import_region_localhost_localdomain_918176><br>Removiendo vector <tmp_r_import_region_localhost_localdomain_918176><br>Estimated target resolution for input band <wv1500_0_5cm_mean>:<br>0.0022003864773300253<br>Reprojecting <wv1500_0_5cm_mean>...<br>ADVERTENCIA: proj_create() failed for '(null)'<br>ERROR: Unable to initialize coordinate transformation<br>ERROR: Unable to to reproject raster <wv1500_0_5cm_mean></div><div><br></div><div>It imports the raster but fails at re-projection... CRS is unknown in the tif</div><div><br></div><div>gdalinfo /vsicurl/<a href="https://storage.googleapis.com/isric-share-soilgrids/pre-release/wv1500/wv1500_0-5cm_mean.tif" rel="noreferrer nofollow noopener" target="_blank">https://storage.googleapis.com/isric-share-soilgrids/pre-release/wv1500/wv1500_0-5cm_mean.tif</a><br>Warning 1: /vsicurl/<a href="https://storage.googleapis.com/isric-share-soilgrids/pre-release/wv1500/wv1500_0-5cm_mean.tif" rel="noreferrer nofollow noopener" target="_blank">https://storage.googleapis.com/isric-share-soilgrids/pre-release/wv1500/wv1500_0-5cm_mean.tif</a>: This file used to have optimizations in its layout, but those have been, at least partly, invalidated by later changes<br>Driver: GTiff/GeoTIFF<br>Files: /vsicurl/<a href="https://storage.googleapis.com/isric-share-soilgrids/pre-release/wv1500/wv1500_0-5cm_mean.tif" rel="noreferrer nofollow noopener" target="_blank">https://storage.googleapis.com/isric-share-soilgrids/pre-release/wv1500/wv1500_0-5cm_mean.tif</a><br>Size is 159246, 58034<br>Coordinate System is:<br>PROJCRS["unknown",<br> BASEGEOGCRS["unknown",<br> DATUM["World Geodetic System 1984",<br> ELLIPSOID["WGS 84",6378137,298.257223563,<br> LENGTHUNIT["metre",1]],<br> ID["EPSG",6326]],<br> PRIMEM["Greenwich",0,<br> ANGLEUNIT["degree",0.0174532925199433],<br> ID["EPSG",8901]]],<br> CONVERSION["unknown",<br> METHOD["Interrupted Goode Homolosine"],<br> PARAMETER["Longitude of natural origin",0,<br> ANGLEUNIT["degree",0.0174532925199433],<br> ID["EPSG",8802]],</div><div>...<br></div><div><br></div><div>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?</div><div><br></div><div>Thanks much in advance!</div><div>Vero<br></div><div><br></div><div>My system info:</div><div>Versión de GRASS: <a href="http://8.0.dev" rel="noreferrer nofollow noopener" target="_blank">8.0.dev</a> <br>Fecha de compilación: 2021-10-13 <br>GDAL: 3.2.2 <br>PROJ: 7.2.1 <br>GEOS: 3.9.0 <br>SQLite: 3.34.1 <br>Python: 3.9.7 <br>wxPython: 4.0.7 <br>Plataforma: Linux-5.14.9-200.fc34.x86_64-x86_64-with-glibc2.33 </div></div>
</blockquote><br>
</div>