Поворот векторного объекта

uses Math;

var P1, P2 : TDoublePoint;
    angle  : double;
begin
  mapGetPlanePoint(Obj, @p1, 1);
  mapGetPlanePoint(Obj, @p2, 1);
  angle := arctan2(p2.y-p1.y, p2.x-p1.x);  
end;

 // Вычисление азимута участка объекта (стороны)
 // Возвращает величину угла в радианах
 // Если геодезия не поддерживается, то вычисляется дирекционный угол
 // Для последней точки вычисляет направление на первую точку
 // У замкнутых объектов первая и последняя точки совпадают
 // info  - идентификатор объекта карты в памяти
 // number - номер точки, начиная с 1
 // subject - номер подобъекта (если = 0, обрабатывается объект)
 // При ошибке возвращает 0 (при совпадении точек также)

_MAPIMP double _MAPAPI mapSideAzimuth(HOBJ info,long int number,      //16/01/01
long int subject);
_MAPIMP long int _MAPAPI mapSideAzimuthEx(HOBJ info,long int number,
long int subject, double * azimuth);

 // Вычисление дирекционного угла участка объекта (стороны)
 // Возвращает величину угла в радианах
 // Для последней точки вычислет направление на первую точку
 // У замкнутых объектов перва и последн точки совпадают
 // info  - идентификатор объекта карты в памяти
 // number - номер точки, начиная с 1
 // subject - номер подобъекта (если = 0, обрабатывается объект)
 // При ошибке возвращает 0 (при совпадении точек также)

_MAPIMP double _MAPAPI mapSideDirection(HOBJ info,long int number,
long int subject/* = 0*/);

 // Определение угла наклона биссектрисы угла, заданного
 // точками p1,p2,p3 с вершиной в точке p2
 // Определение угла наклона биссектрисы угла, заданного
 // точками p1,p2,p3 с вершиной в точке p2
 // Возвращаемый угол задан относительно оси X, его положительное
 // направление соответствует положительному направлению оси Y

_MAPIMP double _MAPAPI mapBisectorAngle(DOUBLEPOINT * p1, DOUBLEPOINT * p2,
DOUBLEPOINT * p3);

 // Прямая геодезическая задача на эллипсоиде                       // 03/04/12
 // Для расстояния не более 250 км координаты определяются с ошибкой до 0,0001",
 // а обратный азимут - до 0,001", что соответствует триангуляции 1 класса
 // Способ вспомогательной точки по методу Красовского
 // Метод предназначен для расстояний меньше радиуса Земли
 // Вычисления выполняются на текущем эллипсоиде, установленном
 // в документе - mapSetDocProjection
 // hmap     - идентификатор открытых данных
 // b1,l1    - геодезические координаты исходной точки
 // angle1   - азимут на вторую точку
 // distance - расстояние до второй точки
 // b2,l2    - рассчитанные координаты второй точки
 // angle2   - рассчитанный азимут со второй точки на первую
 //            (если angle2, то обратный азимут не вычисляется)
 // При ошибке в параметрах возвращает 0

_MAPIMP  long int _MAPAPI mapDirectPositionComputation(HMAP hmap,
double b1, double l1, double angle1,
double distance,
double * b2, double * l2, double * angle2);

 // Обратная геодезическая задача на эллипсоиде                     // 03/04/12
 // Для расстояния не более 180 градусов по широте
 // Выполняется построение ортодромии функцией mapOrthodromeObject
 // и запрос длины объекта и азимута первого отрезка
 // Точность порядка точности триангуляции 1 класса
 // Вычисления выполняются на текущем эллипсоиде, установленном
 // в документе - mapSetDocProjection
 // hmap     - идентификатор открытых данных
 // b1,l1    - геодезические координаты первой точки
 // b2,l2    - геодезические координаты второй точки
 // angle    - рассчитанный азимут с первой точки на вторую
 // Возвращает расстояние между заданными точками на текущем эллипсоиде
 // При ошибке в параметрах возвращает 0

_MAPIMP double _MAPAPI mapInversePositionComputation(HMAP hmap,
double b1, double l1,
double b2, double l2,
double * angle);

 // Повернуть объект вокруг заданной в прямоугольной системе точки   // 15/01/07
 // на заданный угол
 // info   - идентификатор объекта карты в памяти
 // center - координаты точки, вокруг которой поворачивается объект (метры)
 // angle  - угол поворота по часовой стрелке (радианы, от -PI до +PI)
 // При ошибке возвращает 0

_MAPIMP long int _MAPAPI mapRotateObject(HOBJ info, DOUBLEPOINT * center, double * angle);

  /*********************************************************
  **
  *     ФУНКЦИИ НАСТРОЙКИ ОТОБРАЖЕНИЯ ЭЛЕКТРОННОЙ КАРТЫ    *
  *                     С ПОВОРОТОМ                        *
  **
  *********************************************************/

 // Функция настройки отображения карты с поворотом
 //  hmap     - идентификатор открытых данных
 //  angle    - угол поворота карты в плане с вершиной в
 //             юго-западном углу карты (от -Pi до Pi)
 //  fixation - угол сектора фиксации поворота отображения карты
 //             относительно предыдущего положения (от 0 до Pi/6).
 //             По умолчанию = Pi/18 (10 градусов)
 // Угол fixation используется для минимизации дрожания изображения 
 // при движении по повернутой карте по прямой (или почти по прямой), 
 // когда при последовательном вызове функции подаются близкие 
 // значения угла поворота (angle). В случае, если разность между 
 // текущим углом поворота и требуемым будет меньше fixation, 
 // то новый угол поворота не устанавливается.
 // Возвращает значение установленного угла поворота.
 // При ошибке возвращает 0

_MAPIMP double _MAPAPI mapSetupTurn(HMAP hmap, double angle,
double fixation);

 // Активен ли поворот ?
 // hmap - идентификатор открытых данных
 // Возвращает (1 - активен, 0 - нет)

_MAPIMP long int _MAPAPI mapTurnIsActive(HMAP hmap);

 // Запросить угол поворота
 // hmap - идентификатор открытых данных
 // Возвращает значения от -Pi до Pi

_MAPIMP double _MAPAPI mapGetTurnAngle(HMAP hmap);

  // Найти точку на контурах объекта и подобъектов,
  // ближайшую к заданной
  // hMap    - идентификатор открытой карты,
  // info    - идентификатор объекта в памяти
  // pointin - координаты точки в прямоугольной
  //           системе координат , в метрах на местности
  // Возвращает номер точки метрики за которой расположена
  // точка на контуре или ноль при ошибке
  // Координаты точки (в метрах) помещаются по адресу pointout

_MAPIMP  long int _MAPAPI mapSeekNearVirtualPoint(HMAP hMap, HOBJ info,
DOUBLEPOINT * pointin, DOUBLEPOINT * pointout);
 

    const double x1 = 7400229.988;
    const double y1 = 4131918.366;
    const double x2 = 7415961.944;
    const double y2 = 4253383.378;
    HMAP hMap = mapOpenDataUn(WTEXT("/home/okbadmin/GIS/tmp/cmake/map/map.sit"), GENERIC_WRITE);
    cout << "Is map opened: " << hMap << endl;
    HOBJ object = mapCreateObject(hMap);
    auto findObject = mapSeekObject(hMap, object, mapListName(object), 4);
    cout << " Is found: " << findObject << " listname " << object << endl;
    double angle = 0;
    double Bx1 = x1;
    double Ly1 = y1;
    double Bx2 = x2;
    double Ly2 = y2;
    mapPlaneToGeo(hMap, &Bx1, &Ly1);
    mapPlaneToGeo(hMap, &Bx2, &Ly2);
    auto errorCode = mapInversePositionComputation(hMap, Bx1, Ly1, Bx2, Ly2, &angle);
    auto currentAzimuth = mapSideAzimuth(findObject, 1, 0);
    double centerX = Bx2;
    double centerY = Ly2;
    double directionX = 0;
    double directionY = 0;
    mapDirectPositionComputation(hMap, centerX, centerY, angle, 100, &directionX, &directionY, 0);
    mapGeoToPlane(hMap, &centerX, &centerY);
    mapGeoToPlane(hMap, &directionX, &directionY);
    cout.precision(std::numeric_limits<double>::digits10);
    cout << " plane " << centerX << " " << centerY << " " << directionX << " " << directionY << endl;
    cout << " angle: " << angle << " azumuth: " << currentAzimuth << std::endl;
    mapUpdatePointPlane(findObject, centerX, centerY, 1);
    mapUpdatePointPlane(findObject, directionX, directionY, 2);
    auto isCommited = mapCommitObject(findObject);
    cout << " isCommited " << isCommited << std::endl;