Система координат MTW-файла

//--------------------------------------------------------------------
// Инициализация
//--------------------------------------------------------------------
void _fastcall TTranslate::InitCylindricalEqualSpaced()
{
  // Вычисление постоянных по FirstMainParal
  va = ParallelCurvatureRadius(FirstMainParallel) * BigAxis;
  if (va == 0)
    va = 1;
}

//--------------------------------------------------------------------
// Перевод прямоугольных координат в геодезические
// Входные данные:
// x,y - координаты точки в метрах
// Выходные данные:
// b,l - геодезические координаты точки в радианах
//--------------------------------------------------------------------
void _fastcall TTranslate::XY2BL_CylindricalEqualSpaced(double x, double y, double& b, double& l)
{
  y -= FalseEasting;
  x -= FalseNorthing;

  l = y / va + AxisMeridian;  // долготу получили
  b = M_PI_4;               // начальное значение   45

  int sign = 0;
  if (x < 0)
  {
    x = -x;
    sign = 1;
  }

  double abc, abegin, aend;
  abegin = 0; aend = M_PI_2;  //начальный интервал 0--.--90
  for (int i = 2; i < 18; i++)
  {
    // метод половинного деления
    abc = x - ArcLength(b);
    if (abc < 0)
    {
      aend = b;
      b = (b + abegin) / 2.;
    }
    else
    {
      abegin = b;
      b = (b + aend) / 2.;
    }
    if (fabs(abegin - aend) < TRANSLATEPRECISION)
      break;
  }

  if (sign == 0)
    b = (abegin + aend) / 2.0;
  else
    b = -(abegin + aend) / 2.0;
}

//--------------------------------------------------------------------
// Перевод геодезических координат в прямоугольные
// Входные данные:
// b,l - геодезические координаты точки в радианах
// Выходные данные:
// double x,y - прямоугольные координаты точки в метрах
//--------------------------------------------------------------------
void _fastcall TTranslate::BL2XY_CylindricalEqualSpaced(double b, double l, double& x, double& y)
{
  if (b >= 0)
    x = ArcLength(b);                                   // прямоугольные
  else
    x = -ArcLength(-b);                                 // прямоугольные

  x += FalseNorthing;

  y = va * (l - AxisMeridian) + FalseEasting;;          // координаты в метрах
}

//---------------------------------------------------------------------
// Вычисление геодезических координат точки B,L по координатам
// точки X,Y в цилиндрической проекции "Широта/Долгота" на шаре
// Входные данные:
// x,y - координаты точки в метрах
// Выходные данные:
// b,l - широта и долгота точки в радианах
//---------------------------------------------------------------------
void _fastcall TTranslate::XY2BL_LATLON(double x, double y, double& b, double& l)
{
  b = (x - FalseNorthing) / BigAxis + MainPointParallel;
  l = (y - FalseEasting) / BigAxis + AxisMeridian;
}

//---------------------------------------------------------------------
// Вычисление координат X,Y точки по геодезическим координатам
// точки B,L в цилиндрической проекции "Широта/Долгота" на шаре
// Входные данные:
// b,l - геодезические координаты точки в радианах
// Выходные данные:
// x,y - координаты точки в метрах
//---------------------------------------------------------------------
void _fastcall TTranslate::BL2XY_LATLON(double b, double l, double& x, double& y)
{
  x = (b - MainPointParallel) * BigAxis + FalseNorthing;
  y = (l - AxisMeridian) * BigAxis + FalseEasting;
}

  // Запросить данные о проекции матрицы
  // hMap -  идентификатор открытых данных
  // number - номер матрицы в списке открытых матриц
  // mapregister, datumparam, ellipsoidparam - адреса структур, содержащих данные о проекции
  // Структуры MAPREGISTEREX, DATUMPARAM, ELLIPSOIDPARAM описаны в mapcreat.h
  // ttype  - тип локального преобразования координат (см. TRANSFORMTYPE в mapcreat.h) или 0
  // tparm - параметры локального преобразования координат (см. mapcreat.h)
  // При ошибке возвращает ноль

_MAPIMP long int _MAPAPI mapGetMtrProjectionDataPro(HMAP hMap, long int number, MAPREGISTEREX* mapregister,
                                                    DATUMPARAM  *datumparam, ELLIPSOIDPARAM *ellipsoidparam, long int * ttype, LOCALTRANSFORM * tparm);