Differences
This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision Last revision Both sides next revision | ||
tutorial:using_local_saturation [2013/08/14 19:55] admin |
tutorial:using_local_saturation [2013/08/14 19:56] 127.0.0.1 external edit |
||
---|---|---|---|
Line 4: | Line 4: | ||
Load the model ''simulate_deforestation_using_local_saturation.egoml'' from ''\Examples\advanced\local_saturation''. Local saturation prevents a change from occurring within a specific region, in which a class area is greater than an established threshold (Fig.1). This feature is useful to simulate diffusion process as well as to establish a minimum forest remaining area (as established by the Brazilian forest code for private properties). Local saturation can be implemented by reducing the probability of a transition using an asymptotic function as follows: | Load the model ''simulate_deforestation_using_local_saturation.egoml'' from ''\Examples\advanced\local_saturation''. Local saturation prevents a change from occurring within a specific region, in which a class area is greater than an established threshold (Fig.1). This feature is useful to simulate diffusion process as well as to establish a minimum forest remaining area (as established by the Brazilian forest code for private properties). Local saturation can be implemented by reducing the probability of a transition using an asymptotic function as follows: | ||
- | **P2 = P1 * (Li – Oc) / (Li + Oc), case Li >= Oc\\ | + | **P<sub>2</sub> = P<sub>1</sub> * (Li – Oc) / (Li + Oc), case Li >= Oc\\ |
- | P2 = 0, case Li < Oc | + | P<sub>2</sub> = 0, case Li < Oc |
** | ** | ||
Line 12: | Line 12: | ||
For each model step, the amount of deforested cells is calculated in every 3x3 window of the landscape map. In this case, the threshold consists of five cells as follows: | For each model step, the amount of deforested cells is calculated in every 3x3 window of the landscape map. In this case, the threshold consists of five cells as follows: | ||
- | **P2 = P1 * (5 - Oc) / (5 + Oc), case Li >= Oc\\ | + | **P<sub>2</sub> = P<sub>1</sub> * (5 - Oc) / (5 + Oc), case Li >= Oc\\ |
- | P2 = 0, case 5 < Oc** | + | P<sub>2</sub> = 0, case 5 < Oc** |
Line 22: | Line 22: | ||
**nbCount(i1, 3, 3)** | **nbCount(i1, 3, 3)** | ||
- | where nbCount is the neighborhood counting operator, i1 is map # 1, and 3,3 is the window size in cells. You could easily increase the neighborhood size, changing these values. Then, a third //[[:Calculate Map]]// is added to apply the local saturation rule on the probability map as follows: | + | where nbCount is the neighborhood counting operator, i1 is map #1, and 3,3 is the window size in cells. You could easily increase the neighborhood size, changing these values. Then, a third //[[:Calculate Map]]// is added to apply the local saturation rule on the probability map as follows: |
**if v1 - i1 >= 0 then i2 * (v1 - i1) / (v1 + i1) else 0** | **if v1 - i1 >= 0 then i2 * (v1 - i1) / (v1 + i1) else 0** |