Envirometrix / LandGISmaps

Processing of global environmental layers at various resolutions
93 stars 25 forks source link

How to read the data? Help needed #14

Closed Noonehere11 closed 3 years ago

Noonehere11 commented 5 years ago

https://landgis.opengeohub.org/#/?base=BingMaps%20(Aerial)&center=40.9268,36.0293&zoom=4&opacity=84&layer=sol_texture.class_usda.tt_m&depth=0 .

I am very interested in this map and want to have a look at the worldwide data for the map, but I failed to find it on the website.

Could you please, if it is possible, send me the data for the maps? I want the data for the theme ‘soil texture class’, ‘sand, clay content’ and preferably the ‘carbon content’.

I don't quite understand codes, is it required to use codes and R to read the data?? Thank you in advance

thengl commented 5 years ago

Soil texture class map you can download from: https://doi.org/10.5281/zenodo.1475451 The soil texture classes/codes are explained in: https://github.com/Envirometrix/LandGISmaps/blob/master/tables/sol_texture.class_usda.tt_m_250m_b_1950..2017_v0.1.tif.csv Soil carbon content map you can download from: https://doi.org/10.5281/zenodo.1475457 If you need to download data for smaller areas please use the WCS functionality (see #7 section on WCS). We are working on an R package to make the use of LandGIS grids easier. Let me know if I can still help.

alkalifly commented 5 years ago

Are there any plans to aggregate the data to coarser resolutions like SoilGrids?

thengl commented 5 years ago

If you carefully follow the WCS request builder you can produce on-the-fly any resolution you wish. How to do this? Open the LandGIS Geoserver (https://geoserver.opengeohub.org/landgisgeoserver/), then select "Demos" -> "WCS request builder" and follow these instructions.

Even faster way to produce aggregated GeoTIFFs is to request them via the URL command e.g. to resample clay content at 30 cm depth to 1 km resolution I would use (see instructions on how to use the WCS in URL):

https://geoserver.opengeohub.org/landgisgeoserver/ows?service=WCS&version=2.0.1
&request=GetCoverage
&coverageId=predicted250m:sol_clay.wfraction_usda.3a1a1a_m_250m_b30..30cm_1950..2017_v0.2
&subset=Lat(41,45)&subset=Long(32,35)
&scalefactor=0.25
thengl commented 5 years ago

To produce the whole world at e.g. 5 km you can use:

https://geoserver.opengeohub.org/landgisgeoserver/ows?service=WCS&version=2.0.1
&request=GetCoverage
&coverageId=predicted250m:sol_clay.wfraction_usda.3a1a1a_m_250m_b30..30cm_1950..2017_v0.2
&scalefactor=0.05

Note, you can also convert to other projection systems on-the-fly.

alkalifly commented 5 years ago

Brilliant; thank you for the helpful examples.