def _rays(obj: Rays) -> dict:
if not isinstance(obj, gpd.GeoSeries):
return {'Not a GeoSeries': True}
else:
return {
'Missing geometries': obj.is_empty.any(),
'Expecting Linestrings': not obj.geom_type.eq("LineString").all(),
'Missing z coordinates': not obj.has_z.all(),
'More than 2 points in Linestring': np.not_equal(
pygeos.count_coordinates(obj.array.data),
2 * len(obj)).any(),
}
def rays(rays: gpd.GeoSeries) -> None:
# errors
tests = {
'Missing geometries':
rays.is_empty.any(),
'Expecting Linestrings':
not rays.geom_type.eq("LineString").all(),
'Missing z coordinates':
not pygeos.has_z(pygeos.io.from_shapely(rays)).all(),
'More than 2 points in Linestring':
np.not_equal(pygeos.count_coordinates(pygeos.io.from_shapely(rays)),
2 * len(rays)).any(),
}
_raise(tests)
return None
def test_count_coords(geoms, count):
actual = count_coordinates(np.array(geoms, np.object_))
assert actual == count