Fact-checked by Grok 2 weeks ago

Superellipsoid

A superellipsoid is a three-dimensional geometric solid that generalizes the ellipsoid, defined implicitly by the equation \left| \frac{x}{a} \right|^r + \left| \frac{y}{b} \right|^r + \left| \frac{z}{c} \right|^r = 1, where a, b, and c are the semi-axes lengths along the respective coordinate directions, and r > 0 is a parameter controlling the shape's "squareness."^[1] When r = 2, the superellipsoid reduces to a standard ellipsoid (including the sphere if a = b = c); for r > 2, it approaches a rectangular box with rounded edges, while for $0 < r < 2, it becomes more pinched and star-like at the poles.^[1] The surface can also be parameterized using spherical coordinates, facilitating its use in computational modeling:

\begin{align*} x &= a \cdot |\cos \eta|^{2/r} \cdot \operatorname{sgn}(\cos \eta) \cdot \cos \theta, \\ y &= b \cdot |\cos \eta|^{2/r} \cdot \operatorname{sgn}(\cos \eta) \cdot \sin \theta, \\ z &= c \cdot |\sin \eta|^{2/r} \cdot \operatorname{sgn}(\sin \eta), \end{align*}

where \eta \in [-\pi/2, \pi/2] and \theta \in [0, 2\pi).^[2] The concept of the superellipsoid originated in the mid-20th century as an extension of the superellipse, a two-dimensional curve independently described by French mathematician Gabriel Lamé in the 19th century^[1] but popularized in design contexts by Danish polymath Piet Hein in the 1950s and 1960s. Hein applied superellipses to practical designs, such as furniture and tabletops, to blend the smoothness of circles with the utility of rectangles, and extended the idea to three dimensions with his 1965 "superegg"—a balanced, egg-shaped object with r \approx 2.5 that neither rolls nor tips over easily.^[2] In mathematics, superellipsoids belong to the broader family of superquadrics, which encompass both surfaces and volumes with adjustable exponents for flexible shape generation.^[3] Superellipsoids gained prominence in computer graphics and vision through the work of Alan H. Barr, who in 1981 formalized them as superquadrics and introduced angle-preserving transformations to deform them while maintaining desirable properties like uniformity.^[3] These primitives are valued for their ability to model a wide range of natural and man-made objects—from rounded polyhedra to organic forms—with fewer parameters than polygonal meshes, enabling efficient rendering and animation.^[3] Applications extend to computer vision for object recognition and pose estimation via moment invariants, as well as in physics for approximating nanoparticle geometries and in engineering for optimizing designs like electromagnetic sources.^[4]^[5]^[6] Despite their simplicity, superellipsoids pose computational challenges in areas like ray tracing due to non-quadratic equations, though approximations and hierarchical bounding volumes mitigate this in modern graphics pipelines.^[7]

Definition and Formulation

Implicit Equation

The superellipsoid is defined implicitly as the level set consisting of all points (x, y, z) in \mathbb{R}^3 satisfying the equation

\left| \frac{x}{a} \right|^r + \left| \frac{y}{b} \right|^r + \left| \frac{z}{c} \right|^r = 1,

where a > 0, b > 0, and c > 0 denote the lengths of the semi-axes along the x-, y-, and z-directions, respectively, and r > 0 is the squareness parameter that governs the deviation from quadratic curvature.^[8]^[9] The parameter r determines the geometric character of the surface: when r = 2, the equation describes an ellipsoid; as r increases beyond 2, the surface approaches a rectangular cuboid with sharp edges; conversely, for $0 < r < 2, the surface becomes more pinched and diamond-like along the axes.^[8]^[2] In the normalized case where a = b = c = 1, the equation simplifies to the unit superellipsoid

\left| x \right|^r + \left| y \right|^r + \left| z \right|^r = 1,

which serves as a standard reference for symmetric shapes centered at the origin with unit extent along each axis.^[8]^[10] This implicit formulation arises as a direct extension of the two-dimensional superellipse, given by \left| \frac{x}{a} \right|^r + \left| \frac{y}{b} \right|^r = 1, into three dimensions by appending the corresponding z-term while retaining the same exponent r across all axes to preserve consistent squareness.^[8]^[9]

Parametric Representation

The parametric representation of a superellipsoid provides a method to generate points on its surface using two angular parameters, facilitating rendering, sampling, and geometric computations in fields such as computer graphics. For an axis-aligned superellipsoid centered at the origin with semi-axes lengths a, b, and c along the x, y, and z directions, respectively, and shape parameter r > 0, the surface points (x, y, z) are given by

\begin{align*} x &= a \cdot |\cos \eta|^{2/r} \cdot \operatorname{sgn}(\cos \eta) \cdot \cos \theta, \\ y &= b \cdot |\cos \eta|^{2/r} \cdot \operatorname{sgn}(\cos \eta) \cdot \sin \theta, \\ z &= c \cdot |\sin \eta|^{2/r} \cdot \operatorname{sgn}(\sin \eta), \end{align*}

where \eta \in [-\pi/2, \pi/2] and \theta \in [0, 2\pi).^[2] This formulation assumes a single squareness parameter r for all directions, corresponding to an isotropic superellipsoid. The signum function \operatorname{sgn}(u) = 1 if u > 0, -1 if u < 0, and $0 if u = 0 is incorporated to preserve the correct orientation and sign of the trigonometric terms, ensuring the parameterization covers all quadrants without discontinuities or reflections. Without it, raising negative values to non-integer powers (common for r \neq 2) could yield complex results or incorrect symmetries; the absolute value handles the magnitude, while \operatorname{sgn} restores the direction. Here, \eta serves as an analog to latitude, controlling the meridional (pole-to-pole) variation, while \theta acts like longitude, governing the azimuthal (equatorial) distribution around the surface. Unlike the standard sphere (r = 2), where the mapping is uniform and bijective, deviations in r introduce distortions: for small r < 2/3, the parameterization can exhibit singularities at the poles (\eta = \pm \pi/2), where the Jacobian determinant vanishes, leading to degenerate tangent planes and challenges in uniform sampling or normal computation.^[2] To position, orient, and scale the superellipsoid arbitrarily, the basic parametric point (x, y, z) can be transformed via an affine mapping: (x', y', z')^\top = R (x, y, z)^\top + t, where R is a $3 \times 3 rotation matrix and t is the translation vector. This generalization allows flexible posing while maintaining the underlying shape for modeling complex objects.

Special Cases

Sphere and Ellipsoid

A superellipsoid with exponent r = 2 and equal semi-axes a = b = c reduces to a sphere of radius a, as the implicit equation simplifies to \left( \frac{x}{a} \right)^2 + \left( \frac{y}{a} \right)^2 + \left( \frac{z}{a} \right)^2 = 1.^[1]^[11] For arbitrary positive semi-axes a, b, c > 0 and r = 2, the superellipsoid becomes a standard ellipsoid, governed by the quadratic equation \left( \frac{x}{a} \right)^2 + \left( \frac{y}{b} \right)^2 + \left( \frac{z}{c} \right)^2 = 1.^[1]^[11] This case connects superellipsoids directly to classical quadratic surfaces in geometry.^[1] Ellipsoids with r = 2 exhibit unique properties, including smoothness everywhere due to their polynomial nature, positive Gaussian curvature at all points, which classifies them as elliptic surfaces, and preservation under affine transformations, meaning the image of an ellipsoid under an affine map remains an ellipsoid.^[12]^[13]^[14] The case r = 2 corresponds exactly to the ellipsoid, where the exponent unifies generalizations of circles and ellipses in two dimensions to their three-dimensional counterparts, providing a foundational smooth, bounded quadratic form.^[1]^[11]

Cuboid and Other Limits

As the exponent r in the superellipsoid formulation approaches infinity, the shape converges to a cuboid, or rectangular box, aligned with the coordinate axes and having dimensions $2a \times 2b \times 2c, where a, b, and c are the semi-axes parameters.^[15]^[16] This limit reflects the increasing "squareness" of the bounding superellipses in each cross-section, with the surface area approaching that of the cuboid's six faces and the volume stabilizing at $8abc.^[15] Geometrically, the convergence occurs in the Hausdorff metric, meaning the superellipsoid Hausdorff distance to the cuboid diminishes to zero as r \to \infty, enabling precise approximations of polyhedral forms in modeling applications.^[15] At the opposite extreme, as r \to 0^+, the superellipsoid degenerates into a cross-shaped structure that collapses onto the coordinate axes, forming a spiky, star-like form while remaining star-convex with respect to the origin; for $0 < r < 1, the shapes are non-convex.^[15]^[17] The volume approaches zero for dimensions greater than one, and the shape emphasizes extensions along the axes, with surface features pinching inward dramatically.^[15] This limit highlights the superellipsoid's ability to model highly angular, non-convex forms without self-intersections, though the surface becomes increasingly singular near the axes.^[17] For intermediate values, when r = 1 and a = b = c, the superellipsoid reduces to a diamond-like regular octahedron, a Platonic solid with eight equilateral triangular faces.^[18] More generally, superellipsoids exhibit non-smoothness for r \neq 2, as the curvature vanishes or becomes undefined at the axial intersections due to singularities in the parametric representation.^[16] Specifically, for r < 2/3, sharp cusps form at the axes, where the Gaussian curvature becomes infinite, leading to pointed features that enhance the polyhedral character while preserving the implicit equation's continuity.^[15]^[16] These boundary behaviors underscore the superellipsoid's versatility in transitioning between smooth ellipsoids and faceted polyhedra, with visual implications for rendering angular objects in computational geometry.^[18]

Properties

Volume Calculation

The volume V enclosed by a superellipsoid defined by the implicit equation \left| \frac{x}{a} \right|^r + \left| \frac{y}{b} \right|^r + \left| \frac{z}{c} \right|^r \leq 1, where a, b, c > 0 are the semi-axes lengths and r > 0 is the squareness parameter, is given by

V = 8abc \frac{\left[ \Gamma\left(1 + \frac{1}{r}\right) \right]^3}{\Gamma\left(1 + \frac{3}{r}\right)},

with \Gamma denoting the gamma function.^[19] This formula arises from evaluating the triple integral over the region bounded by the implicit equation, leveraging the eightfold symmetry across the octants to restrict computation to the positive octant and multiply by 8. Through substitutions such as u = (x/a)^r, v = (y/b)^r, w = (z/c)^r with u + v + w \leq 1 and u, v, w \geq 0, the integral transforms into a form expressible via the multivariate beta function, which factors into products of univariate beta integrals B(m, n) = \int_0^1 t^{m-1} (1-t)^{n-1} \, dt = \frac{\Gamma(m) \Gamma(n)}{\Gamma(m+n)} along each dimension, ultimately yielding the gamma function expression above.^[19]^[2] For the special case r = 2, the superellipsoid reduces to a standard ellipsoid, and the formula simplifies to the well-known V = \frac{4}{3} \pi a b c, since \Gamma\left(\frac{3}{2}\right) = \frac{\sqrt{\pi}}{2} and \Gamma\left(\frac{5}{2}\right) = \frac{3 \sqrt{\pi}}{4}, confirming \frac{\left[ 2 \Gamma\left(\frac{3}{2}\right) \right]^3}{\Gamma\left(\frac{5}{2}\right)} = \frac{4 \pi}{3}.^[19] In the limit as r \to \infty, the superellipsoid approaches a cuboid (rectangular box) with side lengths $2a, $2b, $2c, and the volume tends to $8abc, as \Gamma\left(1 + \frac{1}{r}\right) \to 1 and \Gamma\left(1 + \frac{3}{r}\right) \to 1.^[19] Numerical evaluation of the volume for arbitrary r > 0 relies on computing the gamma function values, which are well-defined via analytic continuation for non-integer arguments but can pose challenges in precision for very small r (where arguments approach 1 from above) or large r (requiring asymptotic expansions). Standard libraries implement gamma via methods like the Lanczos approximation, ensuring accuracy to machine precision; for specific non-integer r, series expansions related to the hypergeometric function {}_1F_1 (the confluent hypergeometric function underlying certain gamma representations) provide efficient approximations when direct evaluation is unstable.^[19]

Surface Area and Curvature

The surface area S of a superellipsoid lacks a closed-form expression for the general squareness parameter r > 0. Computations typically rely on numerical quadrature over the parametric representation or by segmenting the surface into regions, such as the projections in the first octant (S_{xy}, S_{xz}, S_{yz}), with the total area obtained by multiplying the sum by 8. Approximations, such as those based on series expansions or empirical relations like S \approx 2 A_n(1,1) [a^{s(n)} b^{s(n)} c^{s(n)}]^{1/s(n)} where s(n) is a fitted exponent, achieve low errors (e.g., mean absolute percentage error of 0.177%). For the special case r=2, reducing to an ellipsoid, the surface area involves elliptic integrals of the first and second kinds; it simplifies to $4\pi a^2 for the sphere when a = b = c. The Gaussian curvature K and mean curvature H of the superellipsoid boundary, defined implicitly by F(x,y,z) = \left|\frac{x}{a}\right|^r + \left|\frac{y}{b}\right|^r + \left|\frac{z}{c}\right|^r - 1 = 0, are derived via formulas for implicit surfaces. Specifically, K = \frac{\nabla F \cdot \adj(H(F)) \cdot \nabla F^T}{|\nabla F|^4}, where \nabla F is the gradient of F and H(F) is its Hessian matrix, with \adj denoting the adjugate; the mean curvature follows as H = \frac{\nabla F \cdot H(F) \cdot \nabla F^T - |\nabla F|^2 \trace(H(F))}{2 |\nabla F|^3}. These expressions incorporate terms scaled by powers of r through the partial derivatives of F, such as \frac{\partial F}{\partial x} = r \left|\frac{x}{a}\right|^{r-1} \frac{\operatorname{sgn}(x)}{a} and second derivatives involving (r-1) factors. Curvature properties vary with r: for r=2 and a=b=c, the sphere exhibits constant positive Gaussian curvature K = 1/a^2. As r increases, curvature becomes more varied, approaching zero on the nearly flat faces of the limiting cuboid as r \to \infty. For r < 1, the parametric representation develops singularities at the coordinate axes, where curvature concentrates sharply due to the pinched geometry. The local differential geometry of the superellipsoid is captured by the first fundamental form from its parametric equations \mathbf{r}(\eta, \omega) = \left( a \left| \cos \eta \right|^{2/r} \operatorname{sgn}(\cos \eta) \left| \cos \omega \right|^{2/r} \operatorname{sgn}(\cos \omega), \, b \left| \cos \eta \right|^{2/r} \operatorname{sgn}(\cos \eta) \left| \sin \omega \right|^{2/r} \operatorname{sgn}(\sin \omega), \, c \left| \sin \eta \right|^{2/r} \operatorname{sgn}(\sin \eta) \right), with -\pi/2 \leq \eta \leq \pi/2 and -\pi \leq \omega \leq \pi. The coefficients are E = \mathbf{r}_\eta \cdot \mathbf{r}_\eta, F = \mathbf{r}_\eta \cdot \mathbf{r}_\omega, and G = \mathbf{r}_\omega \cdot \mathbf{r}_\omega, yielding the metric tensor for arc length and area elements \sqrt{EG - F^2} \, d\eta \, d\omega.

Historical Development

Origins in Lamé Curves

The mathematical foundations of the superellipsoid trace back to the superellipse, a curve introduced by French mathematician Gabriel Lamé in 1818 to analyze specific plane figures such as the astroid and hypocycloids. In his work Examen des différentes méthodes employées pour résoudre les problèmes de géométrie, Lamé generalized the equation of the ellipse by incorporating a variable exponent, allowing for shapes that interpolate between quadratic forms and more angular profiles. This form, known as the Lamé curve, is expressed as

\left| \frac{x}{a} \right|^n + \left| \frac{y}{b} \right|^n = 1,

where a and b are semi-axes lengths, and n > 0 is the exponent (for instance, n=2 recovers the standard ellipse, while n=2/3 yields the astroid as a hypocycloid). Lamé's innovation facilitated the study of these curves in theoretical geometry, particularly for understanding properties like tangency and rectification without relying on transcendental functions.^[20] The transition from two-dimensional superellipses to three-dimensional superellipsoids is a natural extension by adding a third term with the same exponent, generating a solid of consistent "squareness." Early treatments in analytical geometry recognized superellipsoids as versatile generalizations of ellipsoids and quadrics, capable of approximating a range of convex shapes through exponent variation, though without attributing a single inventor to the 3D form.^[2] In the 20th century, Danish designer Piet Hein popularized the concept through practical applications, such as his 1965 "superegg"—a balanced, egg-shaped object based on a superellipsoid with exponent around 2.5. Prior to widespread computation, superellipsoids remained primarily theoretical constructs in pure geometry, serving as tools for interpolating between familiar quadrics like spheres and cuboids.^[1]

Adoption in Computer Graphics

The adoption of superellipsoids in computer graphics began in 1981 with Alan H. Barr's introduction of superquadrics as a family of parametric shapes extending traditional quadric surfaces, enabling efficient modeling of organic forms through few control parameters such as squareness exponents. Published in IEEE Computer Graphics and Applications, this work emphasized their angle-preserving transformations, which preserved surface detail during deformations, making them suitable for parametric solid modeling in early graphics pipelines.^[3] By the 1990s, superellipsoids became a standard primitive in ray-tracing software, particularly POV-Ray, where they supported constructive solid geometry (CSG) operations for combining shapes into complex scenes with rounded features like pseudo-cylinders or squished spheres. This integration facilitated realistic rendering of smooth, non-spherical objects without excessive polygonal approximation, contributing to their popularity in hobbyist and professional visualization tools. Key developments in the 1980s and 1990s further solidified their role, including Franc Solina and Ruzena Bajcsy's 1990 method for recovering superquadric models with global deformations from range images, advancing 3D object reconstruction in vision-integrated graphics.^[21] Concurrently, superellipsoids were incorporated into CAD systems to approximate rounded polyhedra, allowing compact representation of manufacturable forms like blended corners on mechanical parts through combinations with convex hulls.^[22] Post-2000 extensions have leveraged superellipsoids in GPU-accelerated rendering for real-time applications and deformable models in virtual and augmented reality, where their parametric efficiency supports interactive shape manipulation. Ongoing research continues to explore integrations with machine learning for shape representation in computer vision and robotics.

Applications

Geometric Modeling

Superellipsoids serve as versatile primitives in computer-aided design (CAD) and computer animation for compactly representing rounded three-dimensional objects, such as furniture components and vehicle bodies, where the squareness parameter r adjusts the degree of roundness from spherical to more angular forms.^[3] This parametric control enables efficient shape synthesis in design workflows, allowing designers to generate smooth, convex solids with minimal adjustments to axes lengths and exponents.^[23] A key advantage of superellipsoids lies in their simplicity, requiring far fewer parameters than polygonal meshes or non-uniform rational B-splines (NURBS) to model convex bodies, which reduces computational overhead in rendering and manipulation while maintaining high-fidelity approximations of organic curves.^[24] Furthermore, they integrate seamlessly with constructive solid geometry (CSG) techniques, enabling blending and boolean operations to construct complex assemblies from basic primitives without excessive geometric complexity. In graphics applications, superellipsoids excel at modeling natural forms like fruits and eggs, particularly when r < 2, which yields pinched, asymmetric shapes that capture the tapered contours of such objects.^[3] They also find use in reverse engineering, where scanned point clouds of physical artifacts are approximated by fitted superellipsoids to recreate parametric CAD models for manufacturing or analysis.^[25] Despite these strengths, superellipsoids are inherently limited to convex geometries and perform poorly for concave shapes; this drawback was addressed in 1980s and 1990s research through hierarchical compositions and CSG-based assemblies of multiple superellipsoids to approximate non-convex forms.^[3]

Image Processing and Fitting

Superellipsoids serve as effective primitives in computer vision for fitting parametric models to range data or image silhouettes, enabling shape extraction from visual inputs. In human body pose estimation, for instance, superquadrics—encompassing superellipsoids—are fitted to segmented point clouds from multi-camera systems to model parts such as arms, forearms, and the torso, where the exponent parameter r (often denoted as \epsilon) adjusts the squareness to capture deviations from smooth ellipsoidal forms, achieving high inlier rates (e.g., over 65% of points within 2 cm) on real range data.^[26] In medical imaging, superellipsoids approximate organ geometries for segmentation and volume estimation, facilitating simulations and diagnostics. Studies from the late 1980s and 1990s onward have applied them to structures like the heart; for example, iterative superellipsoid fitting to CT slices segments the heart's 3D structure, incorporating partial data to enhance robustness and enable accurate volume computation even in occluded regions.^[27] Recent integrations with deep learning leverage superellipsoids as shape priors in neural networks for 3D reconstruction from 2D images, promoting compact representations that improve generalization on sparse datasets. Methods like learnable superquadrics decompose scenes into primitives via differentiable pipelines, enabling unsupervised parsing and reconstruction with fewer parameters than voxel-based alternatives, as demonstrated in 2024 frameworks that fit superquadrics end-to-end for object representation.^[28] Fitting superellipsoids to visual data faces challenges from occlusion and noise, which degrade parameter recovery in partial or corrupted inputs; robust approaches combine gradient-curvature constraints with algebraic distances to handle incomplete silhouettes, often extending to full superquadrics for explicit orientation modeling to mitigate alignment ambiguities.^[29] These techniques build on geometric modeling principles to adapt superellipsoids for data-driven analysis.

Parameter Recovery

Optimization Techniques

Optimization of superellipsoid parameters from observed data, such as range images or point clouds, typically involves nonlinear least-squares minimization to fit the implicit surface equation to the data points. The Levenberg-Marquardt algorithm is widely used for this purpose due to its robustness in handling the nonlinearities inherent in the superellipsoid model, blending gradient descent and Gauss-Newton methods to ensure convergence even from imperfect initial guesses.^[25]^[30] This approach minimizes the sum of squared residuals between the data points and the model surface, enabling recovery of shape parameters a, b, c, and r along with pose. The objective function can be formulated using either algebraic or geometric distance measures. For algebraic fitting, the residuals are based on the inside-outside function of the superellipsoid, defined implicitly as

f(\mathbf{x}) = \left| \frac{x'}{a} \right|^r + \left| \frac{y'}{b} \right|^r + \left| \frac{z'}{c} \right|^r - 1 = 0,

where \mathbf{x} = (x', y', z') are coordinates transformed by rotation and translation, and the cost is \sum_i f(\mathbf{x}_i)^2.^[31]^[32] This method is computationally efficient but sensitive to scale differences. Alternatively, geometric distance minimizes the Euclidean distance from each point to the nearest surface location, often approximated via iterative projection or radial measures like \sum_i \| \mathbf{x}_i - \mu_i \|^2, where \mu_i is the projected point on the surface; this provides a more accurate fit but increases computational cost.^[30]^[25] To account for object orientation and position, pose parameters—rotation matrix and translation vector—are incorporated into the cost function, transforming input points before evaluation. Rotation is commonly parameterized using unit quaternions to prevent gimbal lock and ensure smooth optimization over the SO(3) manifold, with the transformation applied as \mathbf{x}' = R(\mathbf{q}) (\mathbf{x} - \mathbf{t}), where \mathbf{q} is the quaternion and \mathbf{t} the translation.^[33] This joint optimization over shape and pose parameters totals up to 10 degrees of freedom for a full superellipsoid model. Effective initialization is crucial to avoid local minima in the non-convex optimization landscape. Common strategies include deriving initial scale parameters a, b, c from the bounding box of the data points, capturing the extents along principal axes, or using principal component analysis (PCA) on the point cloud to estimate orientation via eigenvectors of the covariance matrix and seed the size parameters from the spread along those directions. The exponent r is often initialized to 2 (ellipsoid) or based on aspect ratios. These seeds guide the Levenberg-Marquardt iterations toward a global fit.^[25]^[30]

Practical Algorithms

Gradient-based methods for superellipsoid parameter recovery typically employ nonlinear least squares optimization, such as the Levenberg-Marquardt algorithm, to minimize the error between observed data points and the implicit superellipsoid equation.^[31] This approach accelerates convergence in iterative solvers by computing the Jacobian matrix of the implicit form, which captures the partial derivatives of the superellipsoid function with respect to its parameters, enabling efficient updates via Gauss-Newton steps.^[34] The Jacobian is particularly useful for orthogonal distance fitting, where it helps locate points on the surface closest to the data while accounting for the nonlinearity of the exponents.^[34] Robust variants address noise and outliers common in real-world data, such as laser scans. The random sample consensus (RANSAC) algorithm, applied hierarchically, samples subsets of points to hypothesize superellipsoid parameters and refits using least squares on inliers, rejecting outliers to ensure reliable recovery from cluttered range images.^[35] For uncertainty quantification, probabilistic approaches like expectation-maximization with Student-t distributions model noise as heavy-tailed, providing posterior estimates of parameters and their variances, which outperform deterministic methods in corrupted point clouds. Software implementations facilitate practical recovery. In MATLAB, custom toolboxes and scripts, such as the EMS-superquadric_fitting library, support probabilistic fitting of superquadrics to point clouds using optimization routines from the Curve Fitting Toolbox.^[36] While OpenCV provides general tools for contour fitting and ellipse approximation, superellipsoid recovery often requires extensions via custom C++ modules integrating least-squares solvers.^[37] For validation, POV-Ray's built-in superellipsoid primitive allows rendering fitted models to visually assess parameter accuracy against input data.^[38] Evaluation relies on goodness-of-fit metrics like chi-squared statistics from least-squares residuals, measuring statistical consistency between fitted and observed points.^[31] Surface accuracy is quantified using Hausdorff distance, which captures maximum deviations between fitted superellipsoid boundaries and data contours.^[39] Benchmarks on synthetic data demonstrate convergence for squareness parameter r in [0.5, 10], with stable formulations achieving low errors (e.g., <1% deviation) across this range using numerically robust exponent evaluations.^[40]

References

[1]
Generalization of the super ellipsoid concept and its application in ...
In three-dimensional case, super ellipsoid (or super egg) is defined by the equation: (2) | x a | n + | y b | n + | z c | n = 1 representing a solid whose ...
[2]
[PDF] Chapter 2 SUPERQUADRICS AND THEIR GEOMETRIC ...
2.1). On the other hand, when m → 0 the curve takes up the shape of a cross. Superellipses are special cases of curves which are known in analytical geometry ...
[3]
Superquadrics and Angle-Preserving Transformations - IEEE Xplore
Abstract: A new and powerful family of parametric shapes extends the basic quadric surfaces and solids, yielding a variety of useful forms.Missing: paper | Show results with:paper
[4]
[PDF] Moments of Superellipsoids and their Application to Range Image ...
Abstract— Cartesian moments are frequently used global geo- metrical features in computer vision for object pose estimation and recognition.
[5]
A unified analytical description of the geometry of nanoscale second ...
Sep 26, 2006 · Superellipsoids: A unified analytical description of the geometry of nanoscale second-phase particles of any shape | Applied Physics Letters | ...
[6]
Electromagnetic source transformations using superellipse equations
May 11, 2009 · The superellipse equations can be used to describe different shapes with a single set of mathematical functions and are thus ideal to ...
[7]
[PDF] Faster Calculation of Superquadric Shapes - W. Randolph Franklin
Jul 3, 1981 · A. H. Barr, "Superquadrics and Angle-Preserving Trans- formations," IEEE Computer Graphics and Applications,. Vol. 1, No. 1, Jan. 1981, pp.Missing: paper | Show results with:paper
[8]
Super ellipsoids - SingSurf.org
Superellipsoids are three dimensional versions of superellipses which in turn are a cross between a square and a circle.
[9]
Generalization of the Super Ellipsoid Concept and Its Application in ...
Aug 10, 2025 · Super-ellipsoids, defined by replacing the exponent 2 in the tri-axial ellipsoid equation with any positive exponent n, offer a flexible and ...
[10]
Superellipsoid - Maple Help - Maplesoft
A superellipsoid is a 3-dimensional solid whose horizontal cross-sections are superellipses with the same exponent r, and whose vertical cross-sections through ...Missing: definition | Show results with:definition
[11]
Superellipse -- from Wolfram MathWorld
The generalization to a three-dimensional surface is known as a superellipsoid. Superellipses with a=b are also known as Lamé curves or Lamé ovals, and the case ...
[12]
[PDF] Superquadrics and Angle- Preserving Transformations
Superquadrics and Angle-. Preserving Transformations. Alan H. Barr. Rensselaer Polytechnic Institute. Over the past 20 years, a great deal of interest has de ...
[13]
Sampling Superquadric Point Clouds with Normals - arXiv
Feb 14, 2018 · This is the first complete approach for close-to-uniform sampling of superellipsoids and superparaboloids in one single framework.
[14]
[PDF] Generalization of the super ellipsoid concept and its application in ...
The super ellipse and super ellipsoids were apparently independently uncovered by Piet Hein (1905-1996), Danish mathematician, inventor, designer and poet by ...Missing: origin | Show results with:origin
[15]
Characterization of asteroid shapes and stability on their surface ...
Aug 24, 2025 · Super-ellipsoids, defined by replacing the exponent 2 in the tri-axial ellipsoid equation with any positive exponent n, offer a flexible and ...
[16]
[PDF] Gaussian and Mean Curvatures∗
Because the Gaussian curvature is everywhere positive on the ellipsoid, the area of the unit sphere, 4π, is the total Gaussian curvature of the ellipsoid.
[17]
[PDF] Geometry of Curves and Surfaces - CS@Purdue
▷ The Gaussian curvature is K = k1k2. ▷ It is the determinant of the shape ... ▷ Likewise for a smooth or piecewise smooth surface. ▷ χ = 2 for a ...
[18]
Ellipsoidal Calculus - GitHub Pages
Thus, ellipsoids are preserved under affine transformation. If the rows of A are linearly independent (which implies m\leqslant n ), and b=0 , the affine ...
[19]
Superellipsoids in Higher Dimensions
### Summary of Superellipsoid Limits and Properties
[20]
[PDF] Scattering from Superquadric Surfaces - DTIC
16. 3.3 Superellipse curvature cusp with for v = 2.0,2.1 ....... ... Note that in three dimensions the superellipsoid can have two different ... r(2/3). = ...
[21]
[PDF] On the Recovery of Superellipsoids - Columbia Academic Commons
When the 5 parameters are all unity, the superellipsoid is the unit sphere. Superellipsoids can also be defined implicitly as the the vol- ume where the SE ...
[22]
Characterization of asteroid shapes and stability on their surface ...
Sep 22, 2025 · In this paper, the shapes of asteroids and the stability on their surfaces are characterized using super-ellipsoids, i.e. geometric shapes ...Missing: cusps | Show results with:cusps
[23]
Volumes of Generalized Unit Balls - jstor
ball is a circular cylinder in R3. balls shown in FIGURE 1. Here we use Stirling's formula: 1-(x) ~ V2-rxx-1/2e-x, where - m of the quantity on the left and ...Missing: l_p | Show results with:l_p
[24]
(PDF) A generic geometric transformation that unifies a wide range ...
Aug 5, 2025 · Using deformed superellipse equations, 2D objects and their 3D orientations are described with a new parametric method and extended to ...
[25]
From Pythagoras to Fourier and From Geometry to Nature
... 1818 by Gabriel Lamé [66]:. (9.4) ... The superellipse solved the problem. It is neither round nor rectangular, but in between. Yet it is fixed, it is ...
[26]
Recovery of parametric models from range images - IEEE Xplore
Abstract: A method for recovery of compact volumetric models for shape representation of single-part objects in computer vision is introduced.
[27]
[PDF] CAD REPRESENTATIONS OF 3D SHAPES WITH ...
. A superellipsoid is defined as the solution of the general form of the implicit equation [11]:. 1. a z a y a x. 1. 1/2. 2. 2. 2/. 3. 2/. 2. 2/. 1. = ⎟. ⎟. ⎠.
[28]
SuperQuadrics - McGill University
[Barr81] Barr, A.H. (1981) Superquadrics and Angle-Preserving Transformations. IEEE Computer Graphics and Applications, 1, 11-22. [Ferr93] Ferrie, F.P. ...
[29]
[PDF] SHAPE RECOVERY USING EXTENDED SUPERQUADRICS
A superellipse is a closed curve defined by the following simple equation: ... +. = , the following implicit equation can be obtained. (3.12). This function ...
[30]
[PDF] 3D HUMAN BODY POSE ESTIMATION BY SUPERQUADRICS
Abstract: This paper presents a method for 3D Human Body pose estimation. 3D real data of the searched object is acquired by a multi-camera system and ...
[31]
3-D Heart Segmentation on CT Images Using Superellipsoid Model ...
The approach consists of fitting superellipse to each slice and making a 3-D superellipsoid model of the heart. We used an iterative method over a set of given ...
[32]
Superellipse fitting to partial data - ScienceDirect.com
This paper attempts to improve the process of fitting to partial data by combining gradient and curvature information with the standard algebraic distance.
[33]
[PDF] Robust and Accurate Superquadric Recovery: A Probabilistic ...
The al- gorithm first estimates the latent variables from a posterior perspective given the current estimation of the superquadric parameter (E-step). Next, the ...<|control11|><|separator|>
[34]
[PDF] Recovery of parametric models from range images - e-Prints.FRI
Bajcsy and F. Solina, "Three dimensional shape representation revisited," in Proc. First Int. Computer Vision Conf., London, En- gland ...
[35]
https://ieeexplore.ieee.org/document/4209233
[36]
Reconstructing Superquadrics from Intensity and Color Images - PMC
Jul 16, 2022 · Recent research has shown that deep learning methods can be used to accurately reconstruct random superquadrics from both 3D point cloud data and simple depth ...
[37]
(PDF) Orthogonal distance fitting of implicit curves and surfaces
Aug 7, 2025 · ... Jacobian matrix (29) and the linear system (28). can be completed and ... 9d, Table 6). 4.4 Superquadric Fitting. A superquadric (superellipsoid) ...
[38]
https://www.povray.org/documentation/view/3.60/285/
[39]
EMS: A probabilistic approach for superquadric recovery - GitHub
We propose an algorithm to recover a superquadric surface/primitive from a given point cloud, with good robustness, accuracy and efficiency.Missing: toolbox | Show results with:toolbox
[40]
Creating Bounding rotated boxes and ellipses for contours - OpenCV
In this tutorial you will learn how to: Use the OpenCV function cv::minAreaRect; Use the OpenCV function cv::fitEllipse. Theory. Code C++
[41]
2.4.1.11 Superquadric Ellipsoid - Documentation - POV-Ray
It is an extension of the quadric ellipsoid. It can be used to create boxes and cylinders with round edges and other interesting shapes.
[42]
(PDF) Parametric Shape Modeling Using Deformable Superellipses ...
Aug 5, 2025 · We found that the superellipse with simple parametric deformations can efficiently model the prostate shape with the Hausdorff distance error ( ...
[43]
Revisiting Superquadric Fitting: A Numerically Stable Formulation
Sep 9, 2019 · One of their core applications is the estimation of object shape characteristics from a set of discrete samples from the surface of an object.