Skip to main content

Tools for topological processing of polygons vectorized from rasters

Project description

RastaChimp

RastaChimp is a utility library to perform topological processing of polygons vectorized from rasters, as returned by the rasterio.features.shapes function:

  • Simplification, using Douglas-Peucker: Shapely (and GEOS) can be used to simplify polygons but will create gaps and overlaps between polygons that share a boundary (even with preserve_topology=true)
  • Smoothing, using the Chaikin method
  • Densify

RastaChimp decomposes the polygons into edges (using Shapely), applies an operation to them then rebuilds transformed polygons. It works well as long as the output of rasterio.features.shapes is in pixel coordinates (default) or the transform to world coordinates does not have a rotation component, which is usually the case. Because of limited floating precision, it may not work well in the more general case, although it could if the polygons are noded (for example, by using the bpol tool of v.clean in GRASS).

Install

pip install rastachimp

Usage

The function to perform simplification is simplify_dp. As inputs, it takes the tuples (geometry, value) as output from rasterio.features.shapes and the tolerance (in the same unit as the input geometries). It returns an iterator of tuples (geometry, value), just like rasterio.features.shapes.

from rasterio import features
from rastachimp import as_shapely, simplify_dp

# 2D array with classification values
image = np.array(...., dtype=np.uint8)

# no transform => in pixel coordinates
shapes = features.shapes(image)

# convert to Shapely geometry (features.shapes returns a GeoJSON dict)
shapes = as_shapely(shapes)

tolerance = 2
# simpl is an interator of tuples (geometry, value)
simpl = simplify_dp(shapes, tolerance)

Issues

Report issues at https://github.com/gvellut/rastachimp/issues

TODO

  • Better detection and handling of bad inputs
  • Other generalization algorithms (Visvalingam)
  • Other smoothing algorithm (Bezier Curves)

Project details


Release history Release notifications

This version
History Node

0.3

Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Filename, size & hash SHA256 hash help File type Python version Upload date
rastachimp-0.3-py3-none-any.whl (8.2 kB) Copy SHA256 hash SHA256 Wheel py3
rastachimp-0.3.tar.gz (7.1 kB) Copy SHA256 hash SHA256 Source None

Supported by

Elastic Elastic Search Pingdom Pingdom Monitoring Google Google BigQuery Sentry Sentry Error logging AWS AWS Cloud computing DataDog DataDog Monitoring Fastly Fastly CDN SignalFx SignalFx Supporter DigiCert DigiCert EV certificate StatusPage StatusPage Status page