def test_coarse(self):
reference = Raster('../../examples/data.tif')
simulated = Raster('../../examples/data1.tif')
reference.resetMask([2])
simulated.resetMask([2])
eb = EBudget(reference, simulated)
eb.coarse(2)
# W
answer = np.array([[1.0, 0.25], [0.5, 0.25]])
np.testing.assert_array_equal(eb.W, answer)
# Rj
answer1 = np.array([[1.0, 1.0], [0, 0]])
answer3 = np.array([[0, 0], [1.0, 1.0]])
ans = {1: answer1, 3: answer3}
np.testing.assert_equal(eb.Rj, ans)
# Sj
answer1 = np.array([[3.0 / 4, 0.0], [1.0, 0]])
answer3 = np.array([[1.0 / 4, 1.0], [0, 1.0]])
ans = {1: answer1, 3: answer3}
np.testing.assert_equal(eb.Sj, ans)
eb.coarse(2)
# W
answer = np.array([[0.5]])
np.testing.assert_array_equal(eb.W, answer)
# Rj
answer1 = np.array([[(1 + 1.0 / 4) / 2]])
answer3 = np.array([[(1.0 / 2 + 1.0 / 4) / 2]])
ans = {1: answer1, 3: answer3}
np.testing.assert_equal(eb.Rj, ans)
# Sj
answer1 = np.array([[(3.0 / 4 + 0.5) / 2]])
answer3 = np.array([[(1.0 / 4 + 1.0 / 4 + 1.0 / 4) / 2]])
ans = {1: answer1, 3: answer3}
np.testing.assert_equal(eb.Sj, ans)
def test_coarse(self):
reference = Raster('../../examples/data.tif')
simulated = Raster('../../examples/data1.tif')
reference.resetMask([2])
simulated.resetMask([2])
eb = EBudget(reference, simulated)
eb.coarse(2)
# W
answer = np.array([[1.0, 0.25], [0.5, 0.25]])
np.testing.assert_array_equal(eb.W, answer)
# Rj
answer1 = np.array([[1.0, 1.0], [0, 0]])
answer3 = np.array([[0, 0], [1.0, 1.0]])
ans = {1: answer1, 3: answer3}
np.testing.assert_equal(eb.Rj, ans)
# Sj
answer1 = np.array([[3.0/4, 0.0], [1.0, 0]])
answer3 = np.array([[1.0/4, 1.0], [0, 1.0]])
ans = {1: answer1, 3: answer3}
np.testing.assert_equal(eb.Sj, ans)
eb.coarse(2)
# W
answer = np.array([[0.5]])
np.testing.assert_array_equal(eb.W, answer)
# Rj
answer1 = np.array([[(1+1.0/4)/2]])
answer3 = np.array([[(1.0/2 + 1.0/4)/2]])
ans = {1: answer1, 3: answer3}
np.testing.assert_equal(eb.Rj, ans)
# Sj
answer1 = np.array([[(3.0/4 + 0.5)/2]])
answer3 = np.array([[(1.0/4 + 1.0/4 + 1.0/4 )/2]])
ans = {1: answer1, 3: answer3}
np.testing.assert_equal(eb.Sj, ans)
