This is the same reason astronomers like HEALPix tiling for skymaps [0]. Equal areas-per-pixel (important for integrating over areas) with straightforward spherical harmonics calculations and hierarchical extensions that can store images at multiple resolutions by subdividing base cells (LIGO uses one such hierarchical strategy for our low-latency gravitational wave source direction estimates [1]).

[0] https://healpix.sourceforge.io

[1] https://arxiv.org/abs/1508.03634

Thanks for sharing. What do you use it for? Did you use S2 previously?