Получил и в Delphi значения совпадающие с cs2cs (используя функцию pj_transform):
> 5444142,83145087
> 4864854,17866784
Так что, уже легче.
Upd: Показания функции pj_transform изменяются на аналогичные pj_fwd, если в строке инициализации добавить параметр: +nadgrids=@null - в faq пишут, что этот параметр нужен, если у нас отличаются целевой и исходный эллипс.