//script for success rate calculation
// one parameter %1 = weight map resulting from the statistical analysis

del active%1.* -force

// Cross Final with Map: active
Active%1.tbt := TableCross(%1,active,IgnoreUndefs)

//In the cross table, calculate
tabcalc Active%1 npixact:=iff(active=1,npix,0)
tabcalc Active%1 Npcumactive = ColumnCumulative(npixact)
tabcalc Active%1 totalslide = ColumnAggregateSum(npixact,,1)
tabcalc Active%1 totalarea = ColumnAggregateSum(npix,,1)
tabcalc Active%1 percentage:=100*(Npcumactive / totalslide)
tabcalc Active%1 Percentlandslide:=100-percentage
tabcalc Active%1 Npixcumul:= cum(NPix)
tabcalc Active%1 reverse = totalarea -npixcumul
tabcalc Active%1 percentmap = 100*(reverse/totalarea)
// after this display a graph with Percentlandslide az y-axis and Percentmap as x-axis
 