Есть карта изолиний с контрольными точками (например, карта высот). Необходимо проинтерполировать значения функции в узлах прямоугольной решётки (не обязательно параллельной осям координат). Есть ли готовые библиотеки или что можно почитать по этой теме?

Нашёл GSLib — но там нужно делать преобразование координат, потому что решётка во входном должна быть параллельна осям координат. Или изменять код, но не хочу это делать с фортрановским кодом.
А в HPGL — не разобрался с вводом данных. Может кто знает, как задать прямоугольную решётку, в узлах которой нужно считать значения, и набор точек с известными значениями ?