Здравствуйте, nrwl, Вы писали:
N>Насколько я знаю, GDAL с XYZ не работает для растра. Привязка — двумерная, там и параметров-то нет для третьей координаты

Такой код нашел, можно через GCP высоту задать. Но если задать к примеру 1млн. GCP не грохнется ли алгоритм.

    GDAL_GCP *pasGCPs;
    ...
//1st point
    //Unique identifier, often numeric
    pasGCPs->pszId = "0";
    //Informational message or ""
    pasGCPs->pszInfo = NULL;
    //Pixel (x) location of GCP on raster
    pasGCPs->dfGCPPixel = 0;
    //Line (y) location of GCP on raster
    pasGCPs->dfGCPLine = 0;
    //X position of GCP in georeferenced space
    pasGCPs->dfGCPX = 590000.000;
    //Y position of GCP in georeferenced space
    pasGCPs->dfGCPY = 4928000.000;
    //Elevation of GCP, or zero if not known
    pasGCPs->dfGCPZ = 0;
    ...
    GDALSetGCPs(dst, GCP_count,pasGCPs,"WGS84");