Спасибо. Всё получилось!
Теперь ГИС Конструктор может строить азимутальные проекции не хуже ГИС Панорама мини.
Остаётся решить следующий вопрос - устранения наложения объектов невидимого полушария на объекты видимого полушария (со стороны полюса проекции). Ранее для решения этой проблемы Вы предложили следующее:
"Чтобы ограничить область отображения Вы можете на исходной карте (Карта мира цилиндрической проекции) нанести окружность с центром будущей азимутальной проекции. Затем скопировать объекты с обрезанием по заданному объекту (окружности) на другую карту. На промежуточную цилиндрическую или азимутальную. В этом случае не будет заворачиваний объектов за горизонт. Выполнение этой процедуры на лету при отображении в заданной проекции требует существенных временных затрат."
Согласен, что временные затраты на формирование проекции при выполнении фильтрации объектов неизбежно возрастут. Но, ввиду невозможности обхода данной процедуры, с этим придется мириться. Сократить эти затраты можно лишь упрощением алгоритма фильтрации.
Т.е., на настоящий момент для формирования азимутальной проекции имеется следующий код:
<code>
/////////////////////////
hMap=mapOpenData("k0.map",0); // открываем исходную карту
MAPREGISTEREX mapreg;
LISTREGISTER listreg;
memset((void*)&mapreg, 0, sizeof(mapreg));
memset((void*)&listreg, 0, sizeof(listreg));
mapreg.Length = sizeof(MAPREGISTEREX);
listreg.Length = sizeof(LISTREGISTER);
mapRegisterFromMapType(GEOGRAPHIC, &mapreg);
mapreg.EllipsoideKind = SPHERE_WGS_84;
mapreg.MaterialProjection = AZIMUTALEQUIDISTANTGUAM; // выбираем новый тип проекции
mapreg.MainPointParallel = 89 * M_PI / 180.0; // широта полюса проекции
mapreg.AxisMeridian = 0 * M_PI / 180.0; // долгота полюса проекции
mapreg.Scale = 5000000;
mapreg.FlagRealPlace = 1;
mapreg.DeviceCapability = -1;
mapreg.DataProjection = 1
// создаём карту в новой проекции
HMAP hMapNew = mapCreateMapEx("MyMap.sit", "Karta_mira.rsc", &mapreg, &listreg);
HOBJ hObj = mapCreateSiteObject(hMap, hMap);
HOBJ hTemp = mapCreateSiteObject(hMap, hMap);
HSELECT hSelect = mapCreateSiteSelectContext(hMap, hMap);
//// здесь следует задать условия фильтрации объектов ////
int flag = WO_FIRST;
while (mapSeekSiteSelectObject(hMap, hMap, hObj, hSelect, flag) != 0)
{
flag = WO_NEXT;
// Сделать копию объекта, чтобы не нарушить
// последовательность поиска объектов в mapSeekSiteSelectObject
mapReadCopyObject(hTemp, hObj);
// Перенос объекта на другую карту
if (mapChangeObjectMap(hTemp, hMapNew, hMapNew) != 0)
{
// Запись на новую карту
mapCommitObject(hTemp);
}
}
mapFreeObject(hTemp);
mapFreeObject(hObj);
mapDeleteSelectContext(hSelect);
MapSortingWithEvent("MyMap.sit", 0, 0, 0, "./", 0);
/////////////////////////
</code>
Затруднение вызывает задание структуры hSelekt, описывающей область фильтрации объектов.
При этом указанная область может иметь форму окружности (как было предложено Вами) с центром в точке, совпадающей с полюсом проекции, или же можно реализовывать упрощённый алгоритм фильтрации по геодезическим координатам, например, такой:
1.
Зная геодезические координаты полюса проекции (Bo,Lo), вычисляем предельные значения области проецирования (в градусах):
предельная широта с севера Bn = Bo + 90,
предельная широта с юга Bs = Bo - 90,
предельная долгота слева Ll = Lo - 90,
предельная долгота справа Lr = Lo + 90.
Полученные величины приводим к допустимому диапазону их значений:
если Bn > 90, то Bn = 180 - Bn,
если Bn < -90, то Bn = -180 - Bn,
если Ll < -180, то Ll = Ll + 360,
если Lr > 180, то Lr = Lr - 360.
2.
Выполнение, собственно, фильтрации.
На новую проекцию переносятся только объекты, удовлетворяющие хотя бы одному из двух условий.
Условие А:
если Bo < 0, широта объекта должна быть меньше Bs, иначе (если Bo > 0) широта объекта должна быть больше Bn;
при Bo = 0 условие А не рассматривается и работает только условие Б.
Условие Б:
если Bo = 90 или Bo = -90, условие Б не учитывается и работает только условие А,
иначе (если Bo лежит в диапазоне ]-90,90[ ) геодезические координаты объекта должны попадать в диапазон [Bn,Bs] по широте, а по долготе быть
больше Ll и в то же время меньше Lr.
Иными словами, при нахождении полюса проекции на Экваторе, отсечение "невидимого" полушария происходит только по долготе за счет работы условия Б.
При совпадении полюса проекции с одним из географических полюсов отсечение "невидимого" полушария осуществляется только по широте за счет работы условия А.
В иных случаях на новую карту проецируются и зона А, и зона Б за счёт работы обоих условий:
<img src="http://i58.fastpic.ru/big/2015/0122/8a/b72ba0120060dd118d0536c78cd1188a.jpg">
Предполагаю, что описанный алгоритм будет работать несколько быстрее предложенного Вами за счёт того, что в нём нет необходимости вычислять расстояния от каждого из анализируемых объектов до полюса проекции, а достаточно лишь сравнить геодезические координаты объекта с их предельными значениями для выбранного положения полюса проекции. Это лишь незначительно увеличило бы временные затраты на формирование косой азимутальной проекции.
Ввиду всего вышеизложенного, прошу ответить на следующие вопросы:
1.
На Ваш взгляд, способна ли реализация описанного упрощённого алгоритма фильтрации существенно сократить временные затраты на формирование проекции по сравнению с предложенным Вами ранее вариантом построения окружности с центром будущей азимутальной проекции и последующей выборке объектов из заданной области?
Иными словами, какой из алгоритмов видится более предпочтительным с позиции сокращения временных затрат на вычисления?
2.
Каким образом описать в структуре hSelect границы области селекции объектов для предпочитаемого метода фильтрации?
В частности, если использовать для этого метод mapSelectSeekArea(...), то каким образом задавать границы используемого в нём объекта-области поиска (особенно если это окружность)? Обязательно ли при задании границ объекта указывать прямоугольные координаты или существуют методы создания объектов по геодезическим координатам (что в данном случае гораздо удобнее)?