Simple estimator of the irregular open surface area in 3-D
Project description
PyOpenSurfaceArea
Estimates the area of an open surface represented as a point cloud in 3-D space.
Installation
Install it simply by
pip install pyosa
Usage
A simple example of estimating the surface area composed of all points within some radius around a targeted point on the surface of the Stanford bunny is shown below. Run example.py
to reproduce this example.
from matplotlib.patches import Circle
import numpy as np
from scipy import spatial
import pyosa
# laod the data
xyz = np.genfromtxt('bunny100k.xyz', delimiter=' ')
# downsampling the point cloud just a bit for visualization purposes
N = xyz.shape[0]
num = 0.05 * N
mask = np.arange(0, N, int(N/num))
xyz_ds = xyz[mask]
# extract an arbitrary, open surface for which we want to extract the area
point = xyz_ds[np.argmax(xyz_ds[:, 2]), :]
tree = spatial.KDTree(xyz_ds)
r = 0.5
ind = tree.query_ball_point(point, r)
# extraction of the surface area is as simple as
surf = xyz_ds[ind]
area = pyosa.estimate(surf)
# show the resulting figure
circle = patches.Circle(point, r, fc='none', ec='k')
fig, ax = show_pcd(xyz_ds, point, circle) # impementation in example.py
ax.set_title(f'surface area: {area:.2f}')
Shortcomings
The current implementation is limited to open surfaces in 3-D space. The surface is assumed to have the least amount of variance in its "height" direction, that is, the surface is observed from the point of view on the normal at the center of the surface.
License
Project details
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
pyosa-0.0.1.tar.gz
(1.8 MB
view hashes)
Built Distribution
pyosa-0.0.1-py3-none-any.whl
(5.8 kB
view hashes)