Процедура реализации полиномиального преобразования в GDAL

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

Вот бы посмотреть на исходники, чтобы сравнить с теми, что использует QGIS. Кто-то мог бы вытащить ее из исходников?

[Lynx]

Вот что я нашел в последней версии gdal (скачать можно здесь):

Модуль достаточно хорошо описан в комментариях в нем.

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

спасибо, это не совсем то, самое интересное в этом модуле вот эта часть:

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

/* -------------------------------------------------------------------- */
/*      Compute the forward and reverse polynomials.                    */
/* -------------------------------------------------------------------- */
nCRSresult = CRS_compute_georef_equations( &sPoints,
                                      psInfo->adfToGeoX, psInfo->adfToGeoY,
                                      psInfo->adfFromGeoX, psInfo->adfFromGeoY,
                                      nReqOrder );

То есть как я понимаю, используется процедура CRS_compute_georef_equations, у которой на входе точки и порядок полинома и которая выдает прямые и обратные коэффициенты. Можно ли привести еще ее?

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

пардон, невнимательно посмотрел, искомая процедура начинается вот отсюда:

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

/***************************************************************************/
/*
    COMPUTE THE GEOREFFERENCING COEFFICIENTS BASED ON A SET OF CONTROL POINTS
*/
/***************************************************************************/

static int 
CRS_compute_georef_equations (struct Control_Points *cp, 
                                      double E12[], double N12[], 
                                      double E21[], double N21[], 
                                      int order)

Интересно, что она заимствована полностью из GRASS