0

I have a folder containing 100 DEM rasters. I have a point file containing thousands of points. I need to collect raster values at all of these points.

Only some of the rasters intersect with the points, so I need to somehow pull out rasters that intersect with points and then extract values to the points that intersect with the raster.

Improve this question
3
  • 1
    Can you use Spatial Analyst tools?
    – Erica
    Mar 27, 2015 at 16:40
  • I need either an idea to use in model builder or python script. The tools would work great if there was one raster, but to go through 100s of them would take too long.
    – user3407191
    Mar 27, 2015 at 17:07
  • 1
    Have you got the tool working on one raster and need advice on writing it into a loop? Do you need the model/script from scratch? Are you hung up on just the "does the raster intersect" part, or where?
    – Erica
    Mar 27, 2015 at 17:11

3 Answers 3

Reset to default
2

Extract multiple values to points from Spatial Analyst can deal with multiple raster data as input. But I dont know if it is working with 100 Datasets.

Improve this answer
1
  • It would if you batch it or loop it in a script :)
    – Erica
    Mar 27, 2015 at 22:59
2

Your rasters, I guess from your description, are not overlapping. Therfore you should create a raster mosaic (using a gdalbuildvrt or the ArcGIS mosaicing tools) that would behave like a single image. Thenyou can use your tool (extract value to point or extract multivalue to point) only once.

Improve this answer
1
  • 1 up for better understanding of the question.
    – EikeMike
    Apr 1, 2015 at 4:38
0

EMVTP (http://desktop.arcgis.com/en/arcmap/10.3/tools/spatial-analyst-toolbox/extract-multi-values-to-points.htm ExtractMultiValuesToPoints) can extract several attributes, from several rasters at a time. This is how it works:

in_point_feature     = "ResearchPlots"
in_raster_and_fields = [["SoilData", "Texture"],      \
                        ["SoilData", "OrganicMatter"],  \
                        ["DEM", "Value"],               \
                        ["RainRaster", "PredictedRain"] \
                       ]

ExtractMultiValuesToPoints(in_point_feature, in_raster_and_fields, "NONE")

Improve this answer
2
  • Actually, this is part of my question related to EMVTP, posted last Friday.
    – Delonix R.
    Jun 20, 2016 at 12:35
  • The problem (or unexpected return) is that the tool doesn't give you the actual values, but rather the object ID associated to that cell.
    – Delonix R.
    Jun 20, 2016 at 12:37

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.