Oleg Belenkov (Все сообщения пользователя)

#include "mapapi.h"   

// --------------------------------------------------------------
//  Обязательная подпрограмма
// --------------------------------------------------------------
extern "C"
{
 __declspec(dllexport) long int WINAPI GetInterfaceVersion()
 {
    return 0;
 }
}

// --------------------------------------------------------------
//  Точка входа в задачу
// --------------------------------------------------------------
extern "C"
{
__declspec(dllexport) long int WINAPI CallMapFunctionEx(HMAP hmap, TASKPARMEX * parm)
{
  HOBJ info = mapCreateObject(hmap);

  mapDescribeObject(info, SQUARESERVINCODE);

  DOUBLEPOINT center = { 6197830, 7442044 };

  mapBuildEllipse(info, & center, 2000, 1000, 20.*M_PI/180., 32);

  mapCommitObject(info);

  mapFreeObject(info);

  return 0;
}
}

 #define PANASSERT(condition, testname) if (!(condition)) { mapErrorMessage(IDS_DOERROR, testname); return 0; }
 #define PANASSERTCOMPLETED(testname) return 1;
 #define TESTNAME "GeodeticEllipse"

// --------------------------------------------------------------
//  Check building ellipse
// --------------------------------------------------------------
int test_GeodeticEllipse()
{
  MAPREGISTEREX  mapreg;
  DATUMPARAM     datum;
  ELLIPSOIDPARAM ellipsoid;

  PANASSERT(mapGetParametersForEPSG(3395, &mapreg, &datum, &ellipsoid),TESTNAME);

  TEMPHMAP hmap = mapCreateTempSitePro(WS("service.rsc"), &mapreg, &datum, &ellipsoid, 1);
  PANASSERT(hmap, TESTNAME);

  TEMPHOBJ hobj = mapCreateSiteObject(hmap, hmap);
  PANASSERT(hobj, TESTNAME);

  PANASSERT(mapDescribeObject(hobj, SQUARESERVINCODE), TESTNAME);

  DOUBLEPOINT center(6197830, 7442044);

  PANASSERT(mapBuildEllipse(hobj, & center, 2000, 1000, 20.*M_PI/180., 32), TESTNAME);

  PANASSERT(mapCommitObject(hobj), TESTNAME);

  int count = mapPointCount(hobj, 0);
  PANASSERT(count == 33, TESTNAME);

  // Correct result: distance = 9 672.8923824823
  double distance = mapLengthInMap(hobj);
  PANASSERT(Compare(distance, 9672.8923824823) == 0, TESTNAME);

  // Correct result on ellipsoid: distance = 6 391.82228430545
  distance = mapLength(hobj);
  PANASSERT(Compare(distance, 6391.82228430545) == 0, TESTNAME);

  // Correct result: distance = 6 242 890.30451456
  double square = mapSquareInMap(hobj);
  PANASSERT(Compare(square, 6242890.30451456, 0.00001) == 0, TESTNAME);

  // Correct result: distance = 2 725 978.21426187
  square = mapSquare(hobj);
  PANASSERT(Compare(square, 2725978.21426187, 0.00001) == 0, TESTNAME);

  PANASSERTCOMPLETED(TESTNAME);
}

// -----------------------------------------------------------------
// Установить угол поворота карты
// -----------------------------------------------------------------
double TMapMetric::SetTurnAngle(double angle, double fixation)
{
  // Угол близок к нулю (установить 0)
  if (fabs(angle) < DELTANULL)
  {
    TurnActive = 0;
    Angle      = 0;
    return 0;
  }

  // Проверить угол на вхождение в диапазон
  if (angle < 0)
    while (angle < -M_PI)
      angle += M_PI * 2.;
  else
    while (angle > M_PI)
      angle -= M_PI * 2.;

  // Предыдущий угол не равен нулю
  if (fabs(Angle) > DELTANULL)
  {
    double delta;

    // Проверить угол сектора фиксации на вхождение в диапазон
    fixation = max(0.0, min(M_PI/6., fixation));
    delta    = fabs(Angle-angle);

    // Угол близок к предыдущему (оставить предыдущий)
    if (delta < fixation)
      return Angle;

    // Значения близких углов могут иметь разные знаки (около 180 градусов)
    if ((delta -= M_PI*2.) > 0)
    {
      // Угол близок к предыдущему (оставить предыдущий)
      if (delta < fixation)
        return Angle;
    }
  }

  // Установить новое значение угла
  SinTurn    = sin(Angle = angle);
  CosTurn    = cos(Angle);
  TurnActive = 1;

  return Angle;
}

// -----------------------------------------------------------------
// Повернуть точку
// -----------------------------------------------------------------
void _fastcall TMapMetric::Turn(DOUBLEPOINT& point, DOUBLEPOINT& turnpoint)
{
  turnpoint.X = point.X * CosTurn - point.Y * SinTurn;
  turnpoint.Y = point.X * SinTurn + point.Y * CosTurn;
}

// -----------------------------------------------------------------
// Повернуть точку обратно
// -----------------------------------------------------------------
void _fastcall TMapMetric::Inverse(DOUBLEPOINT& point, DOUBLEPOINT& turnpoint)
{
  turnpoint.X =  point.X * CosTurn + point.Y * SinTurn;
  turnpoint.Y = -point.X * SinTurn + point.Y * CosTurn;
}