def test_append_diff_pass(self):
""" Append raster """
centr_ras = Centroids()
centr_ras.set_raster_file(HAZ_DEMO_FL, window= Window(0, 0, 50, 60))
centr_bis = Centroids()
centr_bis.set_raster_file(HAZ_DEMO_FL, window= Window(51, 61, 10, 10))
centr_bis.append(centr_ras)
self.assertAlmostEqual(centr_bis.meta['crs'], DEF_CRS)
self.assertAlmostEqual(centr_bis.meta['transform'].c, -69.33714959699981)
self.assertAlmostEqual(centr_bis.meta['transform'].a, 0.009000000000000341)
self.assertAlmostEqual(centr_bis.meta['transform'].b, 0.0)
self.assertAlmostEqual(centr_bis.meta['transform'].f, 10.42822096697894)
self.assertAlmostEqual(centr_bis.meta['transform'].d, 0.0)
self.assertAlmostEqual(centr_bis.meta['transform'].e, -0.009000000000000341)
self.assertEqual(centr_bis.meta['height'], 71)
self.assertEqual(centr_bis.meta['width'], 61)
def test_append_pass(self):
"""Append points"""
centr = Centroids()
centr.lat, centr.lon, centr.geometry = self.data_vector()
centr_bis = Centroids()
centr_bis.set_lat_lon(np.array([1, 2, 3]), np.array([4, 5, 6]))
with self.assertRaises(ValueError):
centr_bis.append(centr)
centr.geometry.crs = {'init': 'epsg:4326'}
centr_bis.append(centr)
self.assertAlmostEqual(centr_bis.lat[0], 1)
self.assertAlmostEqual(centr_bis.lat[1], 2)
self.assertAlmostEqual(centr_bis.lat[2], 3)
self.assertAlmostEqual(centr_bis.lon[0], 4)
self.assertAlmostEqual(centr_bis.lon[1], 5)
self.assertAlmostEqual(centr_bis.lon[2], 6)
self.assertTrue(np.array_equal(centr_bis.lat[3:], centr.lat))
self.assertTrue(np.array_equal(centr_bis.lon[3:], centr.lon))
