Superquadrics

From HandWiki
Short description: Family of geometric shapes
Some superquadrics.

In mathematics, the superquadrics or super-quadrics (also superquadratics) are a family of geometric shapes defined by formulas that resemble those of ellipsoids and other quadrics, except that the squaring operations are replaced by arbitrary powers. They can be seen as the three-dimensional relatives of the superellipses. The term may refer to the solid object or to its surface, depending on the context. The equations below specify the surface; the solid is specified by replacing the equality signs by less-than-or-equal signs.

The superquadrics include many shapes that resemble cubes, octahedra, cylinders, lozenges and spindles, with rounded or sharp corners.[1] Because of their flexibility and relative simplicity, they are popular geometric modeling tools, especially in computer graphics. It becomes an important geometric primitive widely used in computer vision,[2][3] robotics,[4] and physical simulation.[5]

Some authors, such as Alan Barr, define "superquadrics" as including both the superellipsoids and the supertoroids.[1][6] In modern computer vision literatures, superquadrics and superellipsoids are used interchangeably, since superellipsoids are the most representative and widely utilized shape among all the superquadrics.[2][3] Comprehensive coverage of geometrical properties of superquadrics and methods of their recovery from range images and point clouds are covered in several computer vision literatures.[1][3][7][8]

Formulas

Implicit equation

The surface of the basic superquadric is given by

[math]\displaystyle{ \left|x\right|^r + \left|y\right|^s + \left|z\right|^t =1 }[/math]

where r, s, and t are positive real numbers that determine the main features of the superquadric. Namely:

  • less than 1: a pointy octahedron modified to have concave faces and sharp edges.
  • exactly 1: a regular octahedron.
  • between 1 and 2: an octahedron modified to have convex faces, blunt edges and blunt corners.
  • exactly 2: a sphere
  • greater than 2: a cube modified to have rounded edges and corners.
  • infinite (in the limit): a cube

Each exponent can be varied independently to obtain combined shapes. For example, if r=s=2, and t=4, one obtains a solid of revolution which resembles an ellipsoid with round cross-section but flattened ends. This formula is a special case of the superellipsoid's formula if (and only if) r = s.

If any exponent is allowed to be negative, the shape extends to infinity. Such shapes are sometimes called super-hyperboloids.

The basic shape above spans from -1 to +1 along each coordinate axis. The general superquadric is the result of scaling this basic shape by different amounts A, B, C along each axis. Its general equation is

[math]\displaystyle{ \left|\frac{x}{A}\right|^r + \left|\frac{y}{B}\right|^s + \left|\frac{z}{C}\right|^t = 1. }[/math]

Parametric description

Parametric equations in terms of surface parameters u and v (equivalent to longitude and latitude if m equals 2) are

[math]\displaystyle{ \begin{align} x(u,v) &{}= A g\left(v,\frac{2}{r}\right) g\left(u,\frac{2}{r}\right) \\ y(u,v) &{}= B g\left(v,\frac{2}{s}\right) f\left(u,\frac{2}{s}\right) \\ z(u,v) &{}= C f\left(v,\frac{2}{t}\right) \\ & -\frac{\pi}{2} \le v \le \frac{\pi}{2}, \quad -\pi \le u \lt \pi , \end{align} }[/math]

where the auxiliary functions are

[math]\displaystyle{ \begin{align} f(\omega,m) &{}= \sgn(\sin \omega) \left|\sin \omega \right|^m \\ g(\omega,m) &{}= \sgn(\cos \omega) \left|\cos \omega \right|^m \end{align} }[/math]

and the sign function sgn(x) is

[math]\displaystyle{ \sgn(x) = \begin{cases} -1, & x \lt 0 \\ 0, & x = 0 \\ +1, & x \gt 0 . \end{cases} }[/math]

Spherical product

Barr introduces the spherical product which given two plane curves produces a 3D surface. If [math]\displaystyle{ f(\mu)=\begin{pmatrix}f_1(\mu) \\ f_2(\mu)\end{pmatrix},\quad g(\nu)=\begin{pmatrix}g_1(\nu)\\g_2(\nu)\end{pmatrix} }[/math] are two plane curves then the spherical product is [math]\displaystyle{ h(\mu,\nu) = f(\mu)\otimes g(\nu) = \begin{pmatrix} g_1(\nu)\ f_1(\mu) \\ g_1(\nu)\ f_2(\mu) \\ g_2(\nu) \end{pmatrix} }[/math] This is similar to the typical parametric equation of a sphere: [math]\displaystyle{ \begin{align} x&=x_{0}+r\sin \theta \;\cos \varphi \\ y&=y_{0}+r\sin \theta \;\sin \varphi \qquad (0\leq \theta \leq \pi ,\;0\leq \varphi \lt 2\pi )\\ z&=z_{0}+r\cos \theta \end{align} }[/math] which give rise to the name spherical product.

Barr uses the spherical product to define quadric surfaces, like ellipsoids, and hyperboloids as well as the torus, superellipsoid, superquadric hyperboloids of one and two sheets, and supertoroids.[1]

Plotting code

The following GNU Octave code generates a mesh approximation of a superquadric:

function superquadric(epsilon,a)
  n = 50;
  etamax = pi/2;
  etamin = -pi/2;
  wmax = pi;
  wmin = -pi;
  deta = (etamax-etamin)/n;
  dw = (wmax-wmin)/n;
  [i,j] = meshgrid(1:n+1,1:n+1)
  eta = etamin + (i-1) * deta;
  w   = wmin + (j-1) * dw;
  x = a(1) .* sign(cos(eta)) .* abs(cos(eta)).^epsilon(1) .* sign(cos(w)) .* abs(cos(w)).^epsilon(1);
  y = a(2) .* sign(cos(eta)) .* abs(cos(eta)).^epsilon(2) .* sign(sin(w)) .* abs(sin(w)).^epsilon(2);
  z = a(3) .* sign(sin(eta)) .* abs(sin(eta)).^epsilon(3);

  mesh(x,y,z);
end

See also

References

  1. 1.0 1.1 1.2 1.3 Barr (1 January 1981). "Superquadrics and Angle-Preserving Transformations". IEEE Computer Graphics and Applications 1 (1): 11–23. doi:10.1109/MCG.1981.1673799. ISSN 0272-1716. https://ieeexplore.ieee.org/document/1673799. 
  2. 2.0 2.1 Paschalidou, Despoina; Ulusoy, Ali Osman; Geiger, Andreas (2019). "Superquadrics Revisited: Learning 3D Shape Parsing Beyond Cuboids". 2019 IEEE/CVF Conference on Computer Vision and Pattern Recognition (CVPR). pp. 10336–10345. doi:10.1109/CVPR.2019.01059. ISBN 978-1-7281-3293-8. https://ieeexplore.ieee.org/document/8953499. 
  3. 3.0 3.1 3.2 Liu, Weixiao; Wu, Yuwei; Ruan, Sipu; Chirikjian, Gregory S. (2022). "Robust and Accurate Superquadric Recovery: A Probabilistic Approach". 2022 IEEE/CVF Conference on Computer Vision and Pattern Recognition (CVPR). pp. 2666–2675. doi:10.1109/CVPR52688.2022.00270. ISBN 978-1-6654-6946-3. https://ieeexplore.ieee.org/document/9878948. 
  4. Ruan, Sipu; Wang, Xiaoli; Chirikjian, Gregory S. (2022). "Collision Detection for Unions of Convex Bodies With Smooth Boundaries Using Closed-Form Contact Space Parameterization". IEEE Robotics and Automation Letters 7 (4): 9485–9492. doi:10.1109/LRA.2022.3190629. ISSN 2377-3766. 
  5. Lu, G.; Third, J. R.; Müller, C. R. (2012-08-20). "Critical assessment of two approaches for evaluating contacts between super-quadric shaped particles in DEM simulations" (in en). Chemical Engineering Science 78: 226–235. doi:10.1016/j.ces.2012.05.041. ISSN 0009-2509. Bibcode2012ChEnS..78..226L. https://www.sciencedirect.com/science/article/pii/S0009250912003223. 
  6. Alan H. Barr (1992), Rigid Physically Based Superquadrics. Chapter III.8 of Graphics Gems III, edited by D. Kirk, pp. 137–159
  7. Aleš Jaklič, Aleš Leonardis, Franc Solina (2000) Segmentation and Recovery of Superquadrics. Kluwer Academic Publishers, Dordrecht
  8. Wu, Yuwei; Liu, Weixiao; Ruan, Sipu; Chirikjian, Gregory S. (2022). "Primitive-Based Shape Abstraction via Nonparametric Bayesian Inference". in Avidan, Shai; Brostow, Gabriel; Cissé, Moustapha et al. (in en). Computer Vision – ECCV 2022. Lecture Notes in Computer Science. 13687. Cham: Springer Nature Switzerland. pp. 479–495. doi:10.1007/978-3-031-19812-0_28. ISBN 978-3-031-19812-0. https://link.springer.com/chapter/10.1007/978-3-031-19812-0_28. 

External links