def test_getStat(self):
reference = Raster('../../examples/data.tif')
simulated = Raster('../../examples/data1.tif')
reference.resetMask([2])
simulated.resetMask([2])
eb = EBudget(reference, simulated)
stat = eb.getStat(nIter=3)
ans0 = {
'NoNo': 0.5,
'NoMed': (5.0 * 5 / 8 + 3.0 * 3 / 8) / 8,
'MedMed': 4.0 / 8,
'MedPer': 1.0,
'PerPer': 1.0
}
for k in stat[0].keys():
np.testing.assert_almost_equal(stat[0][k], ans0[k])
ans1 = {
'NoNo': 0.5,
'NoMed': (5.0 / 8 + 5.0 / 32 + 3.0 / 16 + 3.0 / 32) / 2,
'MedMed': 4.0 / 8,
'MedPer': 1.0,
'PerPer': 1.0
}
for k in stat[1].keys():
np.testing.assert_almost_equal(stat[0][k], ans1[k])
def test_getStat(self):
reference = Raster('../../examples/data.tif')
simulated = Raster('../../examples/data1.tif')
reference.resetMask([2])
simulated.resetMask([2])
eb = EBudget(reference, simulated)
stat = eb.getStat(nIter=3)
ans0 = {'NoNo': 0.5, 'NoMed': (5.0*5/8 + 3.0*3/8)/8, 'MedMed': 4.0/8, 'MedPer': 1.0, 'PerPer': 1.0}
for k in stat[0].keys():
np.testing.assert_almost_equal(stat[0][k],ans0[k])
ans1 = {'NoNo': 0.5, 'NoMed': (5.0/8+5.0/32 + 3.0/16 + 3.0/32)/2, 'MedMed': 4.0/8, 'MedPer': 1.0, 'PerPer': 1.0}
for k in stat[1].keys():
np.testing.assert_almost_equal(stat[0][k],ans1[k])
