bounding-shape 0.1.0

Creates shapely shapes around given points.

Bounding Shapes

This library creates Shapely shapes designed to fit closely around a set of given points, and provides functionality to create a grid of points inside those shapes. This simplifies the problem of trying to find a shape that fits around all points, whilst allowing for gaps/multiple shapes to accomodate for situations where points may be clustered around a few areas, or handling points far outside a cluster.

import geopandas as gpd
from shapely import MultiPoint

# Example dataset of rainfall stations covering England
fn = "AllRainfallStations"
imported_gdf = gpd.read_file(f"zip://{fn}.zip!{fn}/{fn}.shp")
imported_gdf.head()

<style scoped> .dataframe tbody tr th:only-of-type { vertical-align: middle; }

.dataframe tbody tr th {
    vertical-align: top;
}

.dataframe thead th {
    text-align: right;
}

</style>

	sample_N	geometry
0	1	POINT (-1.17133 52.18628)
1	1	POINT (-1.1637 52.98624)
2	1	POINT (-1.67298 51.69289)
3	1	POINT (-0.64929 52.49389)
4	1	POINT (-0.99517 52.73152)

# Convert improted data into Multipoint object
reprojected_gdf = imported_gdf.to_crs(crs="EPSG:32631")
points = MultiPoint(list(zip(reprojected_gdf.geometry.x, reprojected_gdf.geometry.y)))
points

from bounding_shape.bounding_shape import create_bounding_area

bounding_area = create_bounding_area(points)
bounding_area

Using the libary to create a set of points inside an area:

from bounding_shape.bounding_shape import create_spaced_points
spaced_points = create_spaced_points(points=points) 
spaced_points

Targeting a number of grid points

If creating a grid and there is a target number of points for the shape to be split into, the get_grid_spacing_for_points method allows this to be specified with the magnitude parameter.

from bounding_shape.bounding_shape import create_spaced_points, get_grid_spacing_for_points
from shapely import MultiPoint

points = MultiPoint([[0.0, 0.0], [0.0, 100.0], [100.0, 100.0], [100.0, 0.0]])

distance = get_grid_spacing_for_points(points=points, magnitude=500)
spaced_points = create_spaced_points(points=points, grid_spacing=distance, buffer=0.0)
spaced_points