Посчитать коэффициенты в GDAL 2-м полиномом

[Максим Дубинин]

Помогите подсчитать какие коэффициенты матрицы преобразования выдаст GDAL, если использовать 2й полином, используя процедуру изложенную вот здесь. Хочу сравнить со своими расчетами в R.

Данные:

Код: Выделить всё

x,y,lat,long
448,942,103,53.333333333
1761,938,103.5,53.3333333333
444,3871,103,52.66666667
3072,928,104,53.33333333
1774,3870,103.5,52.66666667
3104,3858,104 52.66666667

[Lynx]

на сколько мне удалось разобраться в магическом gdb , я получил следующие значения коэффициентов E (easting) и N (northing) во время отладки:

Код: Выделить всё

/* CALCULATE THE FORWARD TRANSFORMATION COEFFICIENTS */

Код: Выделить всё

(gdb) print *E@6
$46 = {102.88556183227459, 0.00038188556741073574, -7.4209626248944338e-05,
  1.291739396167488e-10, -1.6013767259965199e-09, 1.567426932570247e-08}
(gdb) print *N@6
$47 = {53.627320721400444, -1.2430362200880587e-08, -0.00033262232776502663,
  -5.0787758153134746e-10, 2.6246732378126636e-10, 2.1794416731445137e-08}

и:

Код: Выделить всё

/* CALCULATE THE BACKWARD TRANSFORMATION COEFFICIENTS */

Код: Выделить всё

(gdb) print *E@6
$53 = {2340449148192.2378, 5372.3697173522496, -88322344969.844818,
  -1.1563388041795769e-08, -51.500690814790815, 833229719.60298669}
(gdb) print *N@6
$54 = {-5851105321839.6221, 4949.0757496084861, 220805156122.53735,
  -21.999999971091633, -7.7482737706318456, -2083067544.4501576}

[Максим Дубинин]

вот что выдает для прямого преобразования R:

Код: Выделить всё

#an
1.028856e+02  3.818856e-04 -7.420963e-05  1.291739e-10 -1.601377e-09 1.567427e-08
#bn
5.362732e+01 -1.243106e-08 -3.326224e-04 -5.078775e-10  2.624674e-10 2.179443e-08

результат получается - один в один