В геодезии для моделирования гравитационного поля Земли (ГПЗ) традиционно используются сферические функции. Но ГПЗ пространственно нестационарно и может довольно сильно меняться в различных направлениях. Сферические функции по своей природе не приспособлены для отображения локальных особенностей поля. Поэтому для этих задач целесообразно использовать пространственно и частотно локализованные базисные функции. При этом удобно предварительно моделируемый регион разбить на сегменты с почти стационарным полем. Сложность поля в сегменте может быть охарактеризована с помощью матрицы анизотропности, получаемой в результате ковариационного анализа поля. При таком подходе к моделированию может возникнуть проблема плохой согласованности локальных моделей на границах сегментов. Для преодоления этой проблемы в статье предлагается использовать новые базисные функции с метрикой Махаланобиса вместо обычного евклидова расстояния. Использование метрики Махаланобиса, а также квадратичной формы, обобщающей эту метрику, позволяет учитывать структуру поля при определении отстояния точек друг от друга и сделать процесс моделирования непрерывным.