Abstract: Many natural objects exhibit radial or axial symmetry in a single plane. However, a universal tool for simulating and fitting the shapes of such objects is lacking. Herein, we present an R package called 'biogeom' that simulates and fits many shapes found in nature. The package incorporates novel universal parametric equations that generate the profiles of bird eggs, flowers, linear and lanceolate leaves, seeds, starfish, and tree-rings, and three growth-rate equations that generate the profiles of ovate leaves and the ontogenetic growth curves of animals and plants. 'biogeom' includes several empirical datasets comprising the boundary coordinates of bird eggs, fruits, lanceolate and ovate leaves, tree rings, seeds, and sea stars. The package can also be applied to other kinds of natural shapes similar to those in the datasets. In addition, the package includes sigmoid curves derived from the three growth-rate equations, which can be used to model animal and plant growth trajectories and predict the times associated with maximum growth rate. 'biogeom' can quantify the intra- or interspecific similarity of natural outlines, and it provides quantitative information of shape and ontogenetic modification of shape with important ecological and evolutionary implications for the growth and form of the living world.

Abstract: The shape of leaf laminae exhibits considerable diversity and complexity that reflects adaptations to environmental factors such as ambient light and precipitation as well as phyletic legacy. Many leaves appear to be elliptical which may represent a ‘default’ developmental condition. However, whether their geometry truly conforms to the ellipse equation (EE), i.e., (x/a)2 + (y/b)2 = 1, remains conjectural. One alternative is described by the superellipse equation (SE), a generalized version of EE, i.e., |x/a|n +|y/b|n = 1. To test the efficacy of EE versus SE to describe leaf geometry, the leaf shapes of two Michelia species (i.e., M. cavaleriei var. platypetala, and M. maudiae), were investigated using 60 leaves from each species. Analysis shows that the majority of leaves (118 out of 120) had adjusted root-mean-square errors of < 0.05 for the nonlinear fitting of SE to leaf geometry, i.e., the mean absolute deviation from the polar point to leaf marginal points was smaller than 5% of the radius of a hypothesized circle with its area equaling leaf area. The estimates of n for the two species were ˂ 2, indicating that all sampled leaves conformed to SE and not to EE. This study confirms the existence of SE in leaves, linking this to its potential functional advantages, particularly the possible influence of leaf shape on hydraulic conductance.

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: 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: Bamboo is an important component in subtropical and tropical forest communities. The plant has characteristic long lanceolate leaves with parallel venation. Prior studies have shown that the leaf shapes of this plant group can be well described by a simplified version (referred to as SGE-1) of the Gielis equation, a polar coordinate equation extended from the superellipse equation. SGE-1 with only two model parameters is less complex than the original Gielis equation with six parameters. Previous studies have seldom tested whether other simplified versions of the Gielis equation are superior to SGE-1 in fitting empirical leaf shape data. In the present study, we compared a three-parameter Gielis equation (referred to as SGE-2) with the two-parameter SGE-1 using the leaf boundary coordinate data of six bamboo species within the same genus that have representative long lanceolate leaves, with >300 leaves for each species. We sampled 2000 data points at approximately equidistant locations on the boundary of each leaf, and estimated the parameters for the two models. The root–mean–square error (RMSE) between the observed and predicted radii from the polar point to data points on the boundary of each leaf was used as a measure of the model goodness of fit, and the mean percent error between the RMSEs from fitting SGE-1 and SGE-2 was used to examine whether the introduction of an additional parameter in SGE-1 remarkably improves the model’s fitting. We found that the RMSE value of SGE-2 was always smaller than that of SGE-1. The mean percent errors among the two models ranged from 7.5% to 20% across the six species. These results indicate that SGE-2 is superior to SGE-1 and should be used in fitting leaf shapes. We argue that the results of the current study can be potentially extended to other lanceolate leaf shapes.

Abstract: Leaf shape is closely related to economics of leaf support and leaf functions, including light interception, water use, and CO2 uptake, so correct quantification of leaf shape is helpful for studies of leaf structure/function relationships. There are some extant indices for quantifying leaf shape, including the leaf width/length ratio (W/L), leaf shape fractal dimension (FD), leaf dissection index, leaf roundness index, standardized bilateral symmetrical index, etc. W/L ratio is the simplest to calculate, and recent studies have shown the importance of the W/L ratio in explaining the scaling exponent of leaf dry mass vs. leaf surface area and that of leaf surface area vs. leaf length. Nevertheless, whether the W/L ratio could reflect sufficient geometrical information of leaf shape has been not tested. The FD might be the most accurate measure for the complexity of leaf shape because it can characterize the extent of the self-similarity and other planar geometrical features of leaf shape. However, it is unknown how strongly different indices of leaf shape complexity correlate with each other, especially whether W/L ratio and FD are highly correlated. In this study, the leaves of nine Magnoliaceae species (>140 leaves for each species) were chosen for the study. We calculated the FD value for each leaf using the box-counting approach, and measured leaf fresh mass, surface area, perimeter, length, and width. We found that FD is significantly correlated to the W/L ratio and leaf length. However, the correlation between FD and the W/L ratio was far stronger than that between FD and leaf length for each of the nine species. There were no strong correlations between FD and other leaf characteristics, including leaf area, ratio of leaf perimeter to area, fresh mass, ratio of leaf fresh mass to area, and leaf roundness index. Given the strong correlation between FD and W/L, we suggest that the simpler index, W/L ratio, can provide sufficient information of leaf shape for similarly-shaped leaves. Future studies are needed to characterize the relationships among FD and W/L in leaves with strongly varying shape, e.g., in highly dissected leaves.

Abstract: Möbius bands have been studied extensively, mainly in topology. Generalized Möbius-Listing surfaces and bodies providing a full geometrical generalization, is a quite new field, motivated originally by solutions of boundary value problems. Analogous to cutting of the original Möbius band, for this class of surfaces and bodies, results have been obtained when cutting such bodies or surfaces. In general, cutting leads to interlinked and intertwined different surfaces or bodies, resulting in very complex systems. However, under certain conditions, the result of cutting can be a single surface or body, which reduces complexity considerably. Our research is motivated by this reduction of complexity. In the study of cutting Generalized Möbius-Listing bodies with polygons as cross section, the conditions under which a single body results, displaying the Möbius phenomenon of a one-sided body, have been determined for even and odd polygons. These conditions are based on congruence and rotational symmetry of the resulting cross sections after cutting, and on the knife cutting the origin. The Möbius phenomenon is important, since the process of cutting (or separation of zones in a GML body in general) then results in a single body, not in different, intertwined domains. In all previous works it was assumed that the cross section of the GML bodies is constant, but the main result of this paper is that it is sufficient that only one cross section on the whole GML structure meets the conditions for the Möbius phenomenon to occur. Several examples are given to illustrate this.

Abstract: Egg geometry can be described using Preston's equation, which has seldom been used to calculate egg volume (V) and surface area (S) to explore S versus V scaling relationships. Herein, we provide an explicit re-expression of Preston's equation (designated as EPE) to calculate V and S, assuming that an egg is a solid of revolution. The side (longitudinal) profiles of 2221 eggs of six avian species were digitized, and the EPE was used to describe each egg profile. The volumes of 486 eggs from two avian species predicted by the EPE were compared with those obtained using water displacement in graduated cylinders. There was no significant difference in V using the two methods, which verified the utility of the EPE and the hypothesis that eggs are solids of revolution. The data also indicated that V is proportional to the product of egg length (L) and maximum width (W) squared. A 2/3-power scaling relationship between S and V for each species was observed, that is, S is proportional to (LW2)(2/3). These results can be extended to describe the shapes of the eggs of other species to study the evolution of avian (and perhaps reptilian) eggs.

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: The goal of this study is to provide a rigorous tool to quantify the inequality of the leaf size distribution of an individual plant, thereby serving as a reference trait for quantifying plant adaptations to local environmental conditions. The tool to be presented and tested employs three components: (1) a performance equation (PE), which can produce flexible asymmetrical and symmetrical bell-shaped curves, (2) the Lorenz curve (i.e., the cumulative proportion of leaf size vs. the cumulative proportion of number of leaves), which is the basis for calculating, and (3) the Gini index, which measures the inequality of leaf size distribution. We sampled 12 individual plants of a dwarf bamboo and measured the area and dry mass of each leaf of each plant. We then developed a generalized performance equation (GPE) of which the PE is a special case and fitted the Lorenz curve to leaf size distribution using the GPE and PE. The GPE performed better than the PE in fitting the Lorenz curve. We compared the Gini index of leaf area distribution with that of leaf dry mass distribution and found that there was a significant difference between the two indices that might emerge from the scaling relationship between leaf dry mass and area. Nevertheless, there was a strong correlation between the two Gini indices (r2 = 0.9846). This study provides a promising tool based on the GPE for quantifying the inequality of leaf size distributions across individual plants and can be used to quantify plant adaptations to local environmental conditions.

Abstract: Stomata are essential for the exchange of water vapour and atmospheric gases between vascular plants and their external environments. The stomatal geometries of many plants appear to be elliptical. However, prior studies have not tested whether this is a mathematical reality, particularly since many natural shapes that appear to be ellipses are superellipses with greater or smaller edge curvature than predicted for an ellipse. Compared with the ellipse equation, the superellipse equation includes an additional parameter that allows generation of a larger range of shapes. We randomly selected 240 stomata from each of four Magnoliaceae species to test whether the stomatal geometries are superellipses or ellipses. The stomatal geometries for most stomata (943/960) were found to be described better using the superellipse equation. The traditional “elliptical stomata hypothesis” resulted in an underestimation of the area of stomata, whereas the superellipse equation accurately predicted stomatal area. This finding has important implications for the estimation of stomatal area in studies looking at stomatal shape, geometry, and function.

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 number and composition of species in a community can be quantified with alpha-diversity indices, including species richness (R), Simpson's index (D), and the Shannon-Wiener index (HGREEK TONOS). In forest communities, there are large variations in tree size among species and individuals of the same species, which result in differences in ecological processes and ecosystem functions. However, tree size inequality (TSI) has been largely neglected in studies using the available diversity indices. The TSI in the diameter at breast height (DBH) data for each of 999 20 m x 20 m forest census quadrats was quantified using the Gini index (GI), a measure of the inequality of size distribution. The generalized performance equation was used to describe the rotated and right-shifted Lorenz curve of the cumulative proportion of DBH and the cumulative proportion of number of trees per quadrat. We also examined the relationships of alpha-diversity indices with the GI using correlation tests. The generalized performance equation effectively described the rotated and right-shifted Lorenz curve of DBH distributions, with most root-mean-square errors (990 out of 999 quadrats) being < 0.0030. There were significant positive correlations between each of three alpha-diversity indices (i.e., R, D, and H') and the GI. Nevertheless, the total abundance of trees in each quadrat did not significantly influence the GI. This means that the TSI increased with increasing species diversity. Thus, two new indices are proposed that can balance alpha-diversity against the extent of TSI in the community: (1 – GI) x D, and (1 – GI) x H'. These new indices were significantly correlated with the original D and HGREEK TONOS, and did not increase the extent of variation within each group of indices. This study presents a useful tool for quantifying both species diversity and the variation in tree sizes in forest communities, especially in the face of cumulative species loss under global climate change.

Abstract: Many natural shapes exhibit surprising symmetry and can be described by the Gielis equation, which has several classical geometric equations (for example, the circle, ellipse and superellipse) as special cases. However, the original Gielis equation cannot reflect some diverse shapes due to limitations of its power-law hypothesis. In the present study, we propose a generalized version by introducing a link function. Thus, the original Gielis equation can be deemed to be a special case of the generalized Gielis equation (GGE) with a power-law link function. The link function can be based on the morphological features of different objects so that the GGE is more flexible in fitting the data of the shape than its original version. The GGE is shown to be valid in depicting the shapes of some starfish and plant leaves.