More circular statistics, and use those for Distributions of Longitude (and Angle?) #14427

What is the problem this feature will solve?

In discussions on fixing Distribution so that LatitudeDistribution and LongitudeDistribution can be created, @eerovaher noted that the standard pdf_mean, pdf_std, etc., are not really useful for Longitude, given the wrapping. We probably should use those in stats.circstat instead (or at least give the option to), but those do not have median, mad, and various sortings (including histograms) are missing as well (and are obviously tricky!).

Describe the desired outcome

Ideally we implement a full complement of circular statistics. In #14421 (comment), @nstarman linked to an introduction of circular statistics that includes a median, https://insidebigdata.com/2021/02/12/circular-statistics-in-python-an-intuitive-intro/

Additional context

No response