Abstract: Preston's equation is a general model describing the egg shape of birds. The parameters of Preston's equation are usually estimated after re-expressing it as the Todd-Smart equation and scaling the egg's actual length to two. This method assumes that the straight line through the two points on an egg's profile separated by the maximum distance (i.e., the longest axis of an egg's profile) is the mid-line. It hypothesizes that the photographed egg's profile is perfectly bilaterally symmetrical, which seldom holds true because of photographic errors and placement errors. The existing parameter estimation method for Preston's equation considers an angle of deviation for the longest axis of an egg's profile from the mid-line, which decreases prediction errors to a certain degree. Nevertheless, this method cannot provide an accurate estimate of the coordinates of the egg's center, and it leads to sub-optimal parameter estimation. Thus, it is better to account for the possible asymmetry between the two sides of an egg's profile along its mid-line when fitting egg-shape data. In this paper, we propose a method based on the optimization algorithm (optimPE) to fit egg-shape data and better estimate the parameters of Preston's equation by automatically searching for the optimal mid-line of an egg's profile and testing its validity using profiles of 59 bird eggs spanning a wide range of existing egg shapes. We further compared this method with the existing one based on multiple linear regression (lmPE). This study demonstrated the ability of the optimPE method to estimate numerical values of the parameters of Preston's equation and provide the theoretical egg length (i.e., the distance between two ends of the mid-line of an egg's profile) and the egg's maximum breadth. This provides a valuable approach for comparing egg shapes among conspecifics or across different species, or even different classes (e.g., birds and reptiles), in future investigations.

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.