Abstract: According to Thompson’s principle of similarity, the area of an object should be proportional to its length squared. However, leaf area–length data of some plants have been demonstrated not to follow the principle of similarity. We explore the reasons why the leaf area–length allometry deviates from the principle of similarity and examine whether there is a general model describing the relationship among leaf area, width and length. We sampled more than 11,800 leaves from six classes of woody and herbaceous plants and tested the leaf area–length allometry. We compared six mathematical models based on root-mean-square error as the measure of goodness-of-fit. The best supported model described a proportional relationship between leaf area and the product of leaf width and length (i.e., the Montgomery model). We found that the extent to which the leaf area–length allometry deviates from the principle of similarity depends upon the extent of variation of the ratio of leaf width to length. Estimates of the parameter of the Montgomery model ranged between 1/2, which corresponds to a triangular leaf with leaf length as its height and leaf width as its base, and π/4, which corresponds to an elliptical leaf with leaf length as its major axis and leaf width as its minor axis, for the six classes of plants. The narrow range in practice of the Montgomery parameter implies an evolutionary stability for the leaf area of large-leaved plants despite the fact that leaf shapes of these plants are rather different.

Abstract: Many cross-sectional shapes of plants have been found to approximate a superellipse rather than an ellipse. Square bamboos, belonging to the genus Chimonobambusa (Poaceae), are a group of plants with round-edged square-like culm cross sections. The initial application of superellipses to model these culm cross sections has focused on Chimonobambusa quadrangularis (Franceschi) Makino. However, there is a need for large scale empirical data to confirm this hypothesis. In this study, approximately 750 cross sections from 30 culms of C. utilis were scanned to obtain cross-sectional boundary coordinates. A superellipse exhibits a centrosymmetry, but in nature the cross sections of culms usually deviate from a standard circle, ellipse, or superellipse because of the influences of the environment and terrain, resulting in different bending and torsion forces during growth. Thus, more natural cross-sectional shapes appear to have the form of a deformed superellipse. The superellipse equation with a deformation parameter (SEDP) was used to fit boundary data. We find that the cross-sectional shapes (including outer and inner rings) of C. utilis can be well described by SEDP. The adjusted root-mean-square error of SEDP is smaller than that of the superellipse equation without a deformation parameter. A major finding is that the cross-sectional shapes can be divided into two types of superellipse curves: hyperellipses and hypoellipses, even for cross sections from the same culm. There are two proportional relationships between ring area and the product of ring length and width for both the outer and inner rings. The proportionality coefficients are significantly different, as a consequence of the two different superellipse types (i.e., hyperellipses and hypoellipses). The difference in the proportionality coefficients between hyperellipses and hypoellipses for outer rings is greater than that for inner rings. This work informs our understanding and quantifying of the longitudinal deformation of plant stems for future studies to assess the influences of the environment on stem development. This work is also informative for understanding the deviation of natural shapes from a strict rotational symmetry.

Abstract: This publication introduces the audiovisual geometry database with Gielis transformations as initial records for a prototype of the database. A concise overview is given of the rationale behind the database and studying wave phenomena with Gielis transformations. First results on a form of timbral polyphony observed in Gielis curves and future work are briefly discussed.

Abstract: The richness, diversity, connection, depth and pleasure of studying symmetry continue to open doors. Here we report a connection between Coxeter's Apeirogon and the geometry associated with pictorial space, parabolic rotation and dual numbers.

Abstract: Many natural radial symmetrical shapes (e.g., sea stars) follow the Gielis equation (GE) or its twin equation (TGE). A supertriangle (three triangles arranged around a central polygon) represents such a shape, but no study has tested whether natural shapes can be represented as/are supertriangles or whether the GE or TGE can describe their shape. We collected 100 pieces of Koelreuteria paniculata fruit, which have a supertriangular shape, extracted the boundary coordinates for their vertical projections, and then fitted them with the GE and TGE. The adjusted root mean square errors (RMSEadj) of the two equations were always less than 0.08, and >70% were less than 0.05. For 57/100 fruit projections, the GE had a lower RMSEadj than the TGE, although overall differences in the goodness of fit were non-significant. However, the TGE produces more symmetrical shapes than the GE as the two parameters controlling the extent of symmetry in it are approximately equal. This work demonstrates that natural supertriangles exist, validates the use of the GE and TGE to model their shapes, and suggests that different complex radially symmetrical shapes can be generated by the same equation, implying that different types of biological symmetry may result from the same biophysical mechanisms.

Abstract: The Lamé curve is an extension of an ellipse, the latter being a special case. Dr. Johan Gielis further extended the Lamé curve in the polar coordinate system by introducing additional parameters (n1, n2, n3; m): rφ=1Acosm4φn2+1Bsinm4φn3−1/n1, which can be applied to model natural geometries. Here, r is the polar radius corresponding to the polar angle φ; A, B, n1, n2 and n3 are parameters to be estimated; m is the positive real number that determines the number of angles of the Gielis curve. Most prior studies on the Gielis equation focused mainly on its applications. However, the Gielis equation can also generate a large number of shapes that are rotationally symmetric and axisymmetric when A = B and n2 = n3, interrelated with the parameter m, with the parameters n1 and n2 determining the shapes of the curves. In this paper, we prove the relationship between m and the rotational symmetry and axial symmetry of the Gielis curve from a theoretical point of view with the condition A = B, n2 = n3. We also set n1 and n2 to take negative real numbers rather than only taking positive real numbers, then classify the curves based on extremal properties of r(φ) at φ = 0, π/m when n1 and n2 are in different intervals, and analyze how n1, n2 precisely affect the shapes of Gielis curves.

Abstract: Gielis transformations are a generalization of Lame curves. To combine domains, we can make use of the natural alliance between Lame's work and Rvachev's R-functions. A logical next step is the extension to n-valued logic dening dierent partitions.

Abstract: A new class of irregularly shaped dielectric lens antennas with a supershaped microstrip antenna feeder is presented and detailed in this work. The surface of the lens antenna and the feeder shape have been modelled by using the three and two-dimensional Gielis formula, respectively. The antenna design has been carried out by integrating an home-made software tool with the CST Microwave Studio®. The radiation properties of the whole antenna system have been evaluated using a dedicated high-frequency technique based on the tube tracing approximation. Moreover, the effects due to the multiple internal reflections have been properly modeled. The proposed model was applied to study unusual and complex lens antenna systems with the aim to design special radiation characteristics.

Abstract: In this work, a nonlinear dispersive multiphysic model based on Maxwell and asymptotic Smoluchowsky equations has been developed to analyze the electroporation phenomenon induced by pulsed electric field on biological cells. The irregular plasma membrane geometry has been modeled by incorporating in the numerical algorithm the Gielis superformula as well as the dielectric dispersion of the plasma membrane has been modeled using the multi-relaxation Debye-based relationship. The study has been carried out with the aim to compare our model implementing a thin plasma membrane with the simplified model in which the plasma membrane is modeled as a distributed impedance boundary condition. The numerical analysis has been performed exposing the cell to external electric pulses having rectangular shapes. By an inspection of the obtained results, significant differences can be highlighted between the two models confirming the need to incorporate the effective thin membrane into the numerical algorithm to well predict the cell response to the pulsed electric fields in terms of transmembrane voltages and pore densities, especially when the cell is exposed to external nanosecond pulses.

Abstract: Many geometries of plant organs can be described by the Gielis equation, a polar coordinate equation extended from the superellipse equation, . Here, r is the polar radius corresponding to the polar angle φ; m is a positive integer that determines the number of angles of the Gielis curve when φ ∈ [0 to 2π); and the rest of the symbols are parameters to be estimated. The pentagonal radial symmetry of calyxes and corolla tubes in top view is a common feature in the flowers of many eudicots. However, prior studies have not tested whether the Gielis equation can depict the shapes of corolla tubes. We sampled randomly 366 flowers of Vinca major L., among which 360 had five petals and pentagonal corolla tubes, and six had four petals and quadrangular corolla tubes. We extracted the planar coordinates of the outer rims of corolla tubes (in top view) (ORCTs), and then fitted the data with two simplified versions of the Gielis equation with k = 1 and m = 5: (Model 1), and (Model 2). The adjusted root mean square error (RMSEadj) was used to evaluate the goodness of fit of each model. In addition, to test whether ORCTs are radially symmetrical, we correlated the estimates of n2 and n3 in Model 1 on a log-log scale. The results validated the two simplified Gielis equations. The RMSEadj values for all corolla tubes were smaller than 0.05 for both models. The numerical values of n2 and n3 were demonstrated to be statistically equal based on the regression analysis, which suggested that the ORCTs of V. major are radially symmetrical. It suggests that Model 1 can be replaced by the simpler Model 2 for fitting the ORCT in this species. This work indicates that the pentagonal or quadrangular corolla tubes (in top view) can both be modeled by the Gielis equation and demonstrates that the pentagonal or quadrangular corolla tubes of plants tend to form radial symmetrical geometries during their development and growth.

Abstract: Measuring the inequality of leaf area distribution per plant (ILAD) can provide a useful tool for quantifying the influences of intra- and interspecific competition, foraging behavior of herbivores, and environmental stress on plants’ above-ground architectural structures and survival strategies. Despite its importance, there has been limited research on this issue. This paper aims to fill this gap by comparing four inequality indices to measure ILAD, using indices for quantifying household income that are commonly used in economics, including the Gini index (which is based on the Lorenz curve), the coefficient of variation, the Theil index, and the mean log deviation index. We measured the area of all leaves for 240 individual plants of the species Shibataea chinensis Nakai, a drought-tolerant landscape plant found in southern China. A three-parameter performance equation was fitted to observations of the cumulative proportion of leaf area vs. the cumulative proportion of leaves per plant to calculate the Gini index for each individual specimen of S. chinensis. The performance equation was demonstrated to be valid in describing the rotated and right shifted Lorenz curve, given that >96% of root-mean-square error values were smaller than 0.004 for 240 individual plants. By examining the correlation between any of the six possible pairs of indices among the Gini index, the coefficient of variation, the Theil index, and the mean log deviation index, the data show that these indices are closely related and can be used interchangeably to quantify ILAD.

Abstract: Gielis curves (GC) can represent a wide range of shapes and patterns ranging from star shapes to symmetric and asymmetric polygons, and even self intersecting curves. Such patterns appear in natural objects or phenomena, such as flowers, crystals, pollen structures, animals, or even wave propagation. Gielis curves and surfaces are an extension of Lamé curves and surfaces (superquadrics) which have benefited in the last two decades of extensive researches to retrieve their parameters from various data types, such as range images, 2D and 3D point clouds, etc. Unfortunately, the most efficient techniques for superquadrics recovery, based on deterministic methods, cannot directly be adapted to Gielis curves. Indeed, the different nature of their parameters forbids the use of a unified gradient descent approach, which requires initial pre-processings, such as the symmetry detection, and a reliable pose and scale estimation. Furthermore, even the most recent algorithms in the literature remain extremely sensitive to initialization and often fall into local minima in the presence of large missing data. We present a simple evolutionary algorithm which overcomes most of these issues and unifies all of the required operations into a single though efficient approach. The key ideas in this paper are the replacement of the potential fields used for the cost function (closed form) by the shortest Euclidean distance (SED, iterative approach), the construction of cost functions which minimize the shortest distance as well as the curve length using R-functions, and slight modifications of the evolutionary operators. We show that the proposed cost function based on SED and R-function offers the best compromise in terms of accuracy, robustness to noise, and missing data.