На главную... Продукты | Технологии | Классификаторы | Проекты | Скачать | Цены| Форум | Статьи | Обучение | Контакты

Функция mapPlaneToGeo() для проекции Меркатора

Поиск  Пользователи  Правила  Войти
Форум » Настольные приложения » GIS ToolKit, GIS ToolKit Active, ГИС Конструктор для Windows
Страницы: 1
RSS
Функция mapPlaneToGeo() для проекции Меркатора
 
Здравствуйте!

Функция mapPlaneToGeo() возвращает неправильные геодезические координаты.
Карта векторная map:
-  обобщенный тип карты   7  (моpскaя нaвигaционнaя)
-  вид рамки  3 (прямоугольная)
-  система координат  1 (системa кооpдинaт 42 годa, прямоугольная)
-  пpоекция исходного мaтеpиaлa 8 (цилиндpическaя пpямaя paвноугольнaя (Меpкaтоpa))
-  системa высот   1 (бaлтийскaя системa высот, ноль Кpонштaдского футштокa)
-  вид эллипсоидa   1  (Кpaсовского 1942 г.)
-  флаг соответствия проекции 1
-  метрика объектов пpедстaвленa в pеaльных кооpдинaтaх нa местности (метры), высоты задаются в метрах.

Углы (из паспорта, если перевести в градусы):
 ЮЗ широта= 47, долгота= 84
 СЗ широта= 53, долгота= 84
 СВ широта= 53, долгота= 90
 ЮВ широта= 47, долгота= 90
Углы в проекции (из паспорта в метрах) получены из mapCheckAndUpdate():
 ЮЗ x=3719858.8  y=-210170.4
 СЗ x=4373129.8  y=-210170.4
 СВ x=4373129.8  y=210170.4
 ЮВ x=3719858.8  y=210170.4

Параметры проекции:
 p1 := 51*pi/180;  // Первая главная параллель, радианы
 p2 := -p1;        // Вторая главная параллель, радианы
 f0 := 87*pi/180;  // Осевой меридиан, радианы
 p0 := 0;  // Параллель главной точки, радианы

Прямоугольные координаты определяются правильно (стандартными функциями и методами),
а при переводе в геодезические - правильно только справа от осевого меридиана,
слева получается отрицательная долгота, которую можно преобразовать в правильную по формуле:
L_прав := L_расчетн + f0 + L_юз (=L_расчет + 90 + 84), где L_расчет получена из mapPlaneToGeo().

Посоветуйте, как исправить ситуацию:
- что можно изменить в паспорте;
- или может имеется другая функция перевода для проекции Меркатора на эллипсоиде Красовского.

С Уважением, Bolotoff.
 
Доброго времени суток господин Bolotoff !
Судя по всему наши модераторы еще не прочитали Ваш вопрос, и я приняла инициативу на себя.
Отсылаю Ван на ящик карту с Вашими параметрами и код (на Паскале) с функцией которая Вас беспокоит.
Я не знаю что у Вас с библиотеками GisTool, но мой Вам совет - обновите.
В центре карты с параметрами которые Вы указали выше
для точки  X=(4373129.8+3719858.8)/2       Y=0  
получается B=50.094318, L=87.00000 см.код.
Не тот глуп кто не знает, а тот, кто не знает где искать.
 
Здравствуйте!
Поздравляю всех с праздником Рождества Христова!

Библиотеки обновил до версии 9.15 от 24.12.08.
Ситуация не изменилась (отрицательная долгота слева от осевого меридиана).

С Уважением, Bolotoff.
 
Полностью с Вами согласна!
Проблема присутсвует для западной части карты от осевого мередиана !
Не тот глуп кто не знает, а тот, кто не знает где искать.
Страницы: 1
Читают тему (гостей: 1)



© КБ Панорама, 1991-2024

Регистрируясь или авторизуясь на форуме, Вы соглашаетесь с Политикой конфиденциальности