TransWikia.com

Downloading GeoTIFF files from Soilgrids 2.0 using Linux yields files without spatial information

Geographic Information Systems Asked by Emily H on February 9, 2021

I’ve been trying to download global GeoTIFF files from this Soilgrids database in Linux using this tutorial. The process seems to work, but when I try to use the files they have no associated spatial information and are essentially just pictures. I’m an undergrad who is new to coding and does not have access to ArcGIS Desktop.

Here’s the output when I run the gdalinfo function on one of the faulty GeoTIFF files:

Driver: GTiff/GeoTIFF

Files: ocs_0-5cm_mean_5000.tif

Size is 7962, 2902

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]]]]

Data axis to CRS axis mapping: 1,2

Origin = (-19949750.000000000000000,8361000.000000000000000)

Pixel Size = (5000.000000000000000,-5000.000000000000000)

Metadata:

  AREA_OR_POINT=Area

  Code_version=v2.0.0
              Covariates=clm_wcl_srdyrsum,clm_wcl_s04rad,clm_wcl_p10tot,veg_mod_eviyravg,clm_wcl_p09tot,veg_mod_nppy15,    mor_env_demm,clm_wcl_s08rad,clm_wcl_srdyrstd,clm_wcl_s07rad,clm_wcl_bio19,clm_wcl_p05tot,clm_wcl_bio17,mo    r_mrg_twi,clm_mod_lstd03std,clm_mod_lstd04std,clm_wcl_bio16,clm_wcl_bio08,clm_mod_lstd09std,clm_mod_lstd1    1std,clm_wcl_bio18,mor_mrg_vdp,clm_wcl_p04tot,clm_mod_lstd05std,clm_mod_lstd12std,clm_mod_lstd10std,clm_m    od_lstd02std,clm_wcl_p12tot,clm_mod_lstd01std,veg_mod_evient,veg_mod_evirng,veg_mod_evimax,veg_mod_evievn    ,luc_gfc_trely10,mor_mrg_vbf

Litter_layers=FALSE

Model=Quantile Regression Forests

Model_type=ranger

Mtry=15

Number_trees=200

Outputs_version=RUN03

WoSIS_version=Data stream 7

Image Structure Metadata:

COMPRESSION=DEFLATE

INTERLEAVE=BAND

Corner Coordinates:

Upper Left  (-19949750.000, 8361000.000)

Lower Left  (-19949750.000,-6149000.000)

Upper Right (19860250.000, 8361000.000)

Lower Right (19860250.000,-6149000.000)

Center      (  -44750.000, 1106000.000) (  0d51'53.55"W,  9d56' 7.32"N)

Band 1 Block=256x256 Type=Int16, ColorInterp=Gray

NoData Value=-32768

2 Answers

Using the url below.

https://maps.isric.org/mapserv?map=/map/ocs.map&SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&BBOX=-19948500.0,-6147449.0951699,19688885.757062882,8361000.0&CRS=EPSG:152160&WIDTH=1426&HEIGHT=895&LAYERS=ocs_0-30cm_mean&STYLES=&FORMAT=application/openlayers

This can be easily loaded into QGIS (Windows/Linux/Mac) enter image description here

Then export the area of interest you require as a GeoTiff

enter image description here

Answered by Mapperz on February 9, 2021

You have do have spatial information on your file, on your post you have the following information:


Coordinate System is:
PROJCRS["Interrupted_Goode_Homolosine",
BASEGEOGCRS["GCS_WGS_1984_ellipse",
:
:
Origin = (-19949750.000000000000000,8361000.000000000000000)
Pixel Size = (5000.000000000000000,-5000.000000000000000) 

This is the minimal spatial information information necessary for you to work with the file. The file is in Goode Homolosine projection, as indicated on the Soilgrids FAQ - How can I use the Homolosine projection.

You can then transform the geotiff from the Goode Homolosine projection into EPSG:4326 using gdalwrap or gdal_translate

gdalwarp -t_srs EPSG:4326 input.tif output.tif    

This information concern projection transformation is indicated on the soilgrids tutorial.

Hope this helps, the proposed solution of Mapperz is also a good way to get the maps, with the extra bonus of not requiring ArcGIS

Answered by Jorge Mendes on February 9, 2021

Add your own answers!

Ask a Question

Get help from others!

© 2024 TransWikia.com. All rights reserved. Sites we Love: PCI Database, UKBizDB, Menu Kuliner, Sharing RPP