def test_geodataframe(self):
"""Test that constructing a valid Centroids instance from gdf works."""
gdf = gpd.GeoDataFrame(pd.read_excel(HAZ_TEMPLATE_XLS))
gdf.geometry = gpd.points_from_xy(gdf['longitude'], gdf['latitude'])
gdf['elevation'] = np.random.rand(gdf.geometry.size)
gdf['region_id'] = np.zeros(gdf.geometry.size)
gdf['region_id'][0] = np.NaN
gdf['geom'] = gdf.geometry # this should have no effect on centroids
centroids = Centroids.from_geodataframe(gdf)
centroids.check()
self.assertEqual(centroids.geometry.size, 45)
self.assertEqual(centroids.lon[0], 32.57)
self.assertEqual(centroids.lat[0], -25.95)
self.assertEqual(centroids.elevation.size, 45)
self.assertEqual(centroids.on_land.sum(), 44)
self.assertIsInstance(centroids.geometry, gpd.GeoSeries)
self.assertIsInstance(centroids.geometry.total_bounds, np.ndarray)
