VeinsOfTheEarth / rabpro

Delineating watershed basins and computing attribute statistics using Google Earth Engine
https://VeinsOfTheEarth.github.io/rabpro
BSD 3-Clause "New" or "Revised" License
35 stars 8 forks source link

Soil thickness layer(s) #105

Closed jonschwenk closed 2 years ago

jonschwenk commented 2 years ago

Title

Pelletier 2016 Soil Thickness

Description

This is a geotiff from a collection of 6. The full dataset description can be found here . This ingestion is only for the 6th layer called average_soil_and_sedimentary-deposit_thickness.tif. This layer contains "a grid that averages soil and sedimentary deposit thicknesses in meters for users who want a single thickness value that averages across upland hillslopes and valley bottoms." It was made by combining estimates of upland hillslope regolith and soil thicknesses with valley bottom/lowland deposit thicknesses.

Spatial Data Properties (all files) 

Spatial Representation Type: Raster
Number of Bands: 1

Raster Format: GeoTIFF

Number Columns: 43,200
Number Rows: 18,000

Properties specific to each file:

File name | Variable(units) | Data type | NoData value | Minimum valid value | Maximum valid value -- | -- | -- | -- | -- | -- 1) upland_hill-slope_regolith_thickness.tif | regolith  thickness(meters) | Byte | -1 | 0 | 50 2) upland_hill-slope_soil_thickness.tif | soil thickness(meters) | Float | -1 | 0 | 4.2 3) hill-slope_valley-bottom.tif | area within each pixel composed of hillslopes(fraction) | Float | none | 0 | 1 4) upland_valley-bottom_and_lowland_sedimentary_deposit_thickness.tif | sedimentary deposit thickness(meters) | Byte | -1 | 0 | 50 5) land_cover_mask.tif | land cover class(0 - 4) | Byte | none | 0 | 4 6) average_soil_and_sedimentary-deposit_thickness.tif | soil and sedimentary deposit thicknesses(meters) | Byte | -1 | 0 | 50

Citation

Pelletier, J. D., P. D. Broxton, P. Hazenberg, X. Zeng, P. A. Troch,G.-Y. Niu, Z. Williams, M. A. Brunke, and D. Gochis (2016), A gridded global data set of soil, immobile regolith, and sedimentary deposit thicknesses for regional and global land surface modeling, J. Adv. Model. Earth Syst.,8, doi:10.1002/2015MS000526

Projection

4326

Time start

No response

Time end

No response

URL

https://daac.ornl.gov/SOILS/guides/Global_Soil_Regolith_Sediment.html

github-actions[bot] commented 2 years ago

Thank you for your request! Maintainers can take the following steps to accept it:

jonschwenk commented 2 years ago

@jsta While I presently only need the 6th layer specified above, this would make a good test case for ingestion of multiple images into the same imageCollection.

jonschwenk commented 2 years ago

Data had to be converted to int8 to accommodate -1 as a nodata type (I have no idea how -1 was nodata in their original Byte stoarge scheme). Have not tested as 5/6 of the datasets are currently being ingested by GEE.

jonschwenk commented 2 years ago

Test successful! There's some weirdness in the exported id names (many 0s appended):

subbasin_stats.main(datasets, gee_feature_path='users/jonschwenk/galena_box',test=True)
Out[33]: 
([{'type': 'FeatureCollection',
   'columns': {},
   'features': [{'type': 'Feature',
     'geometry': None,
     'id': 'average_soil_and_sedimentary-deposit_thickness_00000000000000000000',
     'properties': {'count': 20286, 'mean': 27.495010165260872}},
    {'type': 'Feature',
     'geometry': None,
     'id': 'hill-slope_valley-bottom_00000000000000000000',
     'properties': {'count': 20286, 'mean': 0.9862760621817886}},
    {'type': 'Feature',
     'geometry': None,
     'id': 'land_cover_mask_00000000000000000000',
     'properties': {'count': 20286, 'mean': 1.722342989921635}},
    {'type': 'Feature',
     'geometry': None,
     'id': 'upland_hill-slope_regolith_thickness_00000000000000000000',
     'properties': {'count': 20286, 'mean': 183.23965342368015}},
    {'type': 'Feature',
     'geometry': None,
     'id': 'upland_hill-slope_soil_thickness_00000000000000000000',
     'properties': {'count': 5985, 'mean': 0.6192558818250185}},
    {'type': 'Feature',
     'geometry': None,
     'id': 'upland_valley-bottom_and_lowland_sedimentary_deposit_thickness_00000000000000000000',
     'properties': {'count': 20286, 'mean': 36.46651460225319}}]}],
 [<Task NSKNH2DJZITSDKZV2E4KD2X4 EXPORT_FEATURES: users-jonschwenk-pelletier2016__b1 (UNSUBMITTED)>])