Accurate assessment of canopy structure is crucial in studying plant-environment interactions. The advancement of functional-structural plant models (FSPM), which incorporate the 3D structure of individual plants, increases the need for a method for accurate mathematical descriptions of leaf shape. A model was developed as an improvement of an existing leaf shape algorithm to describe a large variety of leaf shapes. Modelling accuracy was evaluated using a spatial segmentation method and shape differences were assessed using principal component analysis (PCA) on the optimised parameters. Furthermore, a method is presented to calculate the mean shape of a dataset, intended for obtaining a representative shape for modelling purposes. The presented model is able to accurately capture a large range of single, entire leaf shapes. PCA illustrated the interpretability of the parameter values and allowed evaluation of shape differences. The model parameters allow straightforward digital reconstruction of leaf shapes for modelling purposes such as FSPMs.