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.