Patch area and uniform sampling on the surface of any ellipsoid

Marples, Callum Robert; Williams, Philip Michael

Abstract

Algorithms for generating uniform random points on a triaxial ellipsoid are non-trivial to verify because of the non-analytical form of the surface area. In this paper, a formula for the surface area of an ellipsoidal patch is derived in the form of a one-dimensional numerical integration problem, where the integrand is expressed using elliptic integrals. In addition, analytical formulae were obtained for the special case of a spheroid. The triaxial ellipsoid formula was used to calculate patch areas to investigate a set of surface sampling algorithms. Particular attention was paid to the efficiency of these methods. The results of this investigation show that the most efficient algorithm depends on the required coordinate system. For Cartesian coordinates, the gradient rejection sampling algorithm of Chen and Glotzer is best suited to this task, when paired with Marsaglia’s method for generating points on a unit sphere. For outputs in polar coordinates, it was found that a surface area rejection sampler is preferable.

Citation

Marples, C. R., & Williams, P. M. (2024). Patch area and uniform sampling on the surface of any ellipsoid. Numerical Algorithms, 95(4), 1801-1827. https://doi.org/10.1007/s11075-023-01628-4

• Journal Article Type: Article
• Acceptance Date: Jul 17, 2023
• Online Publication Date: Aug 14, 2023
• Publication Date: 2024-04
• Deposit Date: Aug 15, 2023
• Publicly Available Date: Aug 16, 2023
• Journal: Numerical Algorithms
• Print ISSN: 1017-1398
• Electronic ISSN: 1572-9265
• Publisher: Springer Verlag
• Peer Reviewed: Peer Reviewed
• Volume: 95
• Issue: 4
• Pages: 1801-1827
• DOI: https://doi.org/10.1007/s11075-023-01628-4
• Keywords: Ellipsoid · Surface area · Random sampling
• Public URL: https://nottingham-repository.worktribe.com/output/24404676
• Publisher URL: https://link.springer.com/article/10.1007/s11075-023-01628-4