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.

  • 1
    Can you use Spatial Analyst tools? – Erica Mar 27 '15 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 '15 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 '15 at 17:11
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.

  • It would if you batch it or loop it in a script :) – Erica Mar 27 '15 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.

  • 1 up for better understanding of the question. – EikeMike Apr 1 '15 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")

  • Actually, this is part of my question related to EMVTP, posted last Friday. – Delonix R. Jun 20 '16 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 '16 at 12:37

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

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