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.

As my previous questions have indicated, I am working on an air quality problem where I have a dataset of industrial emissions that consists of a location (long, lat), and an annual volume of emissions in pounds or tonnes.

The end result that I am looking for is a 'dispersion raster' that will provide an aggregate value of emission dispersion (see below) based on overlapping 11KM buffers, where the dispersion value of each buffer is aggregated wherever the buffers (2, up to 57!) overlap.

Details

I am trying to understand how I can model the dispersion of emissions, from the emissions source, through a simple linear decay model where the dispersion at a point (away from the emissions source) is ((1-d)/(11000m)) (from Hystad 2011 - "Environmental Health Perspectives") where d = the distance to the emissions source in meters.

I envision creating a buffer around each emissions point (and have done so, see below) but do not know how to apply the above formula to the buffer. Once the buffer has that formula, thus providing a range of values that degrade to zero as you move from the emissions source, I would like to aggregate all the values (through the vertical or Z) of overlapping buffers.

For instance, the attached image provides a view of the emissions sources (red dots) and the buffers (the grey buffers are 'low values' and the red 'high'). How would I create buffers that decay from the emissions point by ((1-d)/(11000m))? Once these are created, how could I calculate the aggregate (i.e. the exposure from all the buffers that overlap a point) exposure value at a point, as represented by the blue dot.

This is hard to explain, so please ask questions! Thanks!

Red= Emissions Sources. Buffers = Range of dispersion. Blue = How do I calculate an aggregate dispersion volume for here?

share|improve this question
2  
If you have access to any raster-based GIS or image program, your solution is easily found by convolving the point data with the dispersion kernel. For an account of what this is and how it works, please see the three-part article starting at directionsmag.com/features/convolution/129753 and continued at directionsmag.com/features/convolution-part-2-of-3/129751 and directionsmag.com/features/convolution-part-3-of-3/129750 (the links on the site are broken, unfortunately). –  whuber Nov 6 '12 at 22:44
    
Note, that calculating pollutant dispersion by a linear model with buffers is an extremely simplified method. If you want to produce real dispersion maps you can investigate various atmospheric dispersion modeling methods and specific softwares. –  nadya Nov 7 '12 at 2:47
    
Nadya, I appreciate this comment. At this point, a simplified model is a simpler model. Just curious, and I hope you don't mind my asking, where are you studying geomatics? –  Matt Dance Nov 8 '12 at 4:41
    
Geomatics in Portugal, Ecology was in Ukraine. Good luck in this modeling! –  nadya Nov 12 '12 at 1:06

1 Answer 1

So, you want to create raster with values equal to (1-d)/11000 where d is the distance to the nearest point. I think for this you need to pass the following steps.

Firstly you need to create a distance raster where each cell will have just an euclidean distance to the nearest point. It was discussed in this question. Then you can use raster calculator to compute (1-d)/11000. After, you can merge this raster with that one where you aggregated your 11km buffers with the emission values, by sum or multiple, what you prefer. (You don't have the emission component in your formula (1-d)/11000 and I can't suggest something more precise).

share|improve this answer
    
And thanks for this. I will have some time to work thought this tomorrow, and will report back here. –  Matt Dance Nov 8 '12 at 4:42

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.