0

I have several raster layers which represent particular geological formations.

Each layer is a raster dataset in a grey scale. Each cell has a pixel value which is an elevation of the bottom of a particular geological formation. Each cell also has a stretched value which is a thickness of a particular geological formation.

For each location, I would like to have the type of geological formation at each depths. Let's say, at 1 m depth step.

How can this be done? Can I somehow extract the information about layers into a csv file? Or can some manipulations be done in ArcMap?

Improve this question
2
  • Can be done. You also need dem though
    – FelixIP
    Commented Oct 29, 2015 at 18:34
  • @FelixIP Thanks! Very encouraging! Can you please let me know a bit more details?
    – Yule
    Commented Oct 29, 2015 at 21:19

1 Answer 1

Reset to default
1

I don’t have your grey layers. So I generated them: enter image description here

Bottom of Sandstone increases downwards from 0 to 342, similar with basalt (East-West)

For all 3 bottoms compute vert. distance to slice level (DEM-depth), set as flat surface=100 m in this example: enter image description here

Find minimum of distances

enter image description here

Calculate if bottom of Clay is closest to slice

enter image description here

Do the same for sandstone and basalt, using 2 and 3 accordingly

Mosaic results:

enter image description here

Result:

enter image description here

Improve this answer
7
  • thank you so much for the answer! I will try it and let you know if it worked for me.
    – Yule
    Commented Oct 30, 2015 at 10:04
  • Thank you very much again for this example! I am still working though it... Here is how my data looks like: imgur.com/qWIEYjK and imgur.com/eSlTzoB. I think that that the pixel value is an elevation about the sea level. Do you know what stretched value might represent? Or can it be just not useful information.
    – Yule
    Commented Nov 2, 2015 at 13:44
  • I have quite a few different materials, around 64. All materials have the elevation of the bottom value. Hence, I cannot figure out how to find the level of the surface at each point.
    – Yule
    Commented Nov 2, 2015 at 13:48
  • I did some reading and downloaded dem_europe layer from ERSI online. This suppose to be surface elevation layer
    – Yule
    Commented Nov 2, 2015 at 15:23
  • FelixIP, your answer taught me a lot. Thanks again! I am just wondering if you had 64 different materials and some of them are not present in a particular spot, would you still use ArcGIS to do the calculations? Or would you choose to extract data into e.g. Python, than do the required calculations?
    – Yule
    Commented Nov 2, 2015 at 15:52

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.