def test_create_unit_GGB_using_corners(self):
# create small GGB centred on (150.00025,-34.00025)
expectedShape = (1, 1)
scale = 0.00025
origin = (150.0, -34.0)
corner = (150.0 + scale, -34.0 - scale)
ggb = GriddedGeoBox.from_corners(origin, corner)
self.assertEqual(expectedShape, ggb.shape)
self.assertEqual(corner, ggb.corner)
def test_real_world(self):
# Flinders Islet, NSW
flindersOrigin = (150.927659, -34.453309)
flindersCorner = (150.931697, -34.457915)
originShouldBe = flindersOrigin
shapeShouldBe = (19, 17)
cornerShouldBe = (flindersOrigin[0] + shapeShouldBe[1] * 0.00025,
flindersOrigin[1] - shapeShouldBe[0] * 0.00025)
ggb = GriddedGeoBox.from_corners(flindersOrigin, flindersCorner)
assert ggb is not None
self.assertEqual(shapeShouldBe, ggb.shape)
self.assertAlmostEqual(originShouldBe[0], ggb.origin[0])
self.assertAlmostEqual(originShouldBe[1], ggb.origin[1])
self.assertAlmostEqual(cornerShouldBe[0], ggb.corner[0])
self.assertAlmostEqual(cornerShouldBe[1], ggb.corner[1])
def getFlindersIsletGGB():
flindersOrigin = (150.927659, -34.453309)
flindersCorner = (150.931697, -34.457915)
return GriddedGeoBox.from_corners(flindersOrigin, flindersCorner)
