rem ILWIS Script //Created by T. Hengl (http://spatial-analyst.net) //Input: %1 - prediction map (z), %2 - z1 lower value, %3 - z2 upper value, %4 - prediction variance map (e), %5 - e1 lower value (default=0.4), %6 - e2 upper value (default=0.8), %7 - original variance, %8 - name of the output file //Purpose: Visualisation of uncertainty using whiteness. //Last update: 27. February 2006. //derive the normalized error err=sqrt(%4)/sqrt(%7) //stretch the values (z) to the inspection range: tmp_z=(%1-%2)/(%3-%2) tmp_zc.mpr{dom=Value;vr=0.0000:1.0000:0.0001}=iff(tmp_z<=0, 0, iff(tmp_z>1, 1, tmp_z)) //make the hue image tmp_f1=-90-tmp_zc*300 tmp_f2=iff(tmp_f1<=-360, tmp_f1+360, tmp_f1) tmp_H=iff(tmp_f2>=0, tmp_f2*240/360, (tmp_f2+360)*240/360) //stretch error values (e) to the inspection range tmp_e=(err-%5)/(%6-%5) tmp_ec.mpr{dom=Value;vr=0.0000:1.0000:0.0001}=iff(tmp_e<=0, 0, iff(tmp_e>1, 1, tmp_e)) //make the saturation and intensity images tmp_S=240*(1-tmp_ec) tmp_I=125*(1+tmp_ec) calc tmp_?.mpr //visualize the HSI images together %8=colorhsi(tmp_H,tmp_S,tmp_I) //*due to the systematic rounding of numbers in ILWIS, the intensity needs to be rescaled to 255 calc %8.mpr breakdep %8.mpr -force setdescr %8.mpr Predictions visualized with uncertainty z1=%2, z2=%3, e1=%5, e2=%6 del tmp_??.mpr -force del tmp_?.mpr -force del err.mpr -force open %8.mpr -noask