Если ваши точки находятся в пределах одного полушария, вы можете выполнить гномоническую проекцию от сферических до плоских координат, а затем триангулировать, поскольку большие круги становятся прямыми наименьшего расстояния.
простой (если существует) алгоритм, чтобы найти диаграмму Вороного для набора точек на поверхности сферы. Исходный код был бы великолепен. Я человек Delphi (да, я знаю ...), но я тоже ем C-код.