Take the 2-minute tour ×
Geographic Information Systems Stack Exchange is a question and answer site for cartographers, geographers and GIS professionals. It's 100% free, no registration required.

I've a DTM and a shapefile of built-up areas. For a fire risk model, I've to find (within a defined buffer zone) two zones around built-up areas: one that is above/higher and lateral and the other below (downstream) the built-up area. Then I'll classify them according to model weights.

I'm using ArcGIS 10. I've made some attempts with spatial analyst but still haven't found a solution as it could be that in the same (dissolved) buffer zone an area is at the same time above a built-up area and below another built-up area I need to keep the "below" information and discard the "above" one.

share|improve this question

1 Answer 1

You were right about the spatial analyst extension. Depending on some assumptions (see below) with any luck there shouldn't be too much analysis.

  1. Take your shapefile of built areas and turn it into a grid - make sure that it's the same origin and cell size as the original DTM.
  2. Run a Combine between your built areas raster and your DTM - note that the built areas will have null regions so you'll only get data back in the regions where you have built areas.
  3. (WARNING: Big assumption!) Assuming you are only interested in classifying your raster based on height in comparison to the built areas, note the max height in the combine output (in the DTM column), and classify your original DTM based on this - any area of DTM >= max height is upstream or lateral, otherwise downstream.

This doesn't take into account that the region of the DTM might not be directly downstream of the built area, just that it's altitude is lower. If you are interested in directly downstream it becomes a bit more difficult (the following is just one way you could try for your result):

  1. Based on your grid of built areas (from step one above) run a Region Group to separate built areas into contiguous objects. Alternately we can just use the polygons.
  2. Make a new raster as the -ve DTM as we're looking at downstream.
  3. Run Flow direction on your -ve DTM - you may need to Fill or Sink the dataset to smooth it.
  4. Run a Watershed with the new flow direction and the grouped built areas as sources. Watershed selects everything upstream, but we've reversed this, so rather everything downstream is selected.
  5. (WARNING: Another assumption) We can assume that everywhere on the DTM that is not in a watershed is therefor Upstream of our built areas, and can classify as such.
  6. Any built area that edges onto the watershed of another built area is Downstream, and can be classified as such.
  7. Any remaining buildings can be classified as Upstream
  8. Any remaining areas of watershed can be classified as Downstream.

You might want to experiment with expanding to get lateral areas, and depending on the size of the area and the slop angle, firstly classifying your DTM based on a relative flatness. Hope this has given you something to try at least!

share|improve this answer

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

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