ClarkCGA / UDef-ARP

UDef-ARP was developed by Clark Labs, in collaboration with TerraCarbon, to facilitate implementation of the Verra tool, VT0007 Unplanned Deforestation Allocation (UDef-A).
GNU General Public License v3.0
25 stars 14 forks source link

fix(vulnerability_map.py): alternative geometric classification #37

Closed YaoTingYao closed 1 month ago

YaoTingYao commented 1 month ago

Refined alternative geometric classification algorithm to correct zone width assignments.

The original algorithm assigned narrow zones in areas of lowest vulnerability and wider zones in areas of highest vulnerability. This update reverses that logic, ensuring narrow zones are applied to areas of highest vulnerability, with progressively wider zones for lower-vulnerability areas.

According to UDef-A, in the alternative geometric classification, Common Ratio is as follow:

Common Ratio(r) = (LLmax/LLmin)^1/n_classes

Where:
r = Common ratio (unitless)
LLmax = Higher limit of the highest vulnerability class (m) = 2
LLmin = Lower limit of the lowest vulnerability class = 1
n_classes = Number of vulnerability classes = 30

r=(2/1)^1/30=~1.02337389199677

The alternative geometric classification differs from the original geometric classification, as the class ID increases, the bin width decreases (in contrast to the original geometric classification, where both the class ID and bin width increase). The formula for LLc in alternative geometric classification is as follow :

LLc = LLmin + (LLmax - LLmin) * r^c

Where:
LLc = Lower limit of the cth vulnerability class (m)
LLmax = Higher limit of the highest vulnerability class (m) = 2
LLmin = Lower limit of the lowest vulnerability class = 1
c = Class number (unitless)
r = Common ratio (unitless)

Class Limits of the Final Vulnerability Map

Class LLmin LLmax
1 1 1.0456800631315044
2 1.0456800631315044 1.0903167921791634
3 1.0903167921791634 1.1339340169263818
4 1.1339340169263818 1.1765550228835633
5 1.1765550228835633 1.2182025637193186
6 1.2182025637193186 1.258898873407749
7 1.258898873407749 1.2986656780982861
8 1.2986656780982861 1.337524207714422
9 1.337524207714422 1.3754952072875268
10 1.3754952072875268 1.4125989480317984
11 1.4125989480317984 1.4488552381662634
12 1.4488552381662634 1.4842834334896002
13 1.4842834334896002 1.518902447713434
14 1.518902447713434 1.5527307625596205
15 1.5527307625596205 1.5857864376269035
16 1.5857864376269035 1.6180871200322227
17 1.6180871200322227 1.6496500538318088
18 1.6496500538318088 1.6804920892271047
19 1.6804920892271047 1.71062969156042
20 1.71062969156042 1.740078950105126
21 1.740078950105126 1.768855586655083
22 1.768855586655083 1.7969749639178827
23 1.7969749639178827 1.8244520937163908
24 1.8244520937163908 1.8513016450029645
25 1.8513016450029645 1.8775379516906268
26 1.8775379516906268 1.9031750203053737
27 1.9031750203053737 1.9282265374637066
28 1.9282265374637066 1.952705877179373
29 1.952705877179373 1.976626108003225
30 1.976626108003225 2.0