D'ailleurs, l'écriture de ma formule peut être simplifiée en:
d((lat1,lon1), (lat2,lon2)) = min(π-|lat1+lat2|, max(|lat2-lat1|, min(|lon2-lon1|, |lon2-lon1+2π|, |lon2-lon1-2π|)))
(min(a-x,a+x)=a+min(-x,x)=a-|x| )