И ЕЩЕ РАЗ О ПРОЕКЦИЯХ. Теоретические вопросы.

Системы координат, проекции, преобразования, привязка
Boris
Гуру
Сообщения: 4231
Зарегистрирован: 10 апр 2006, 22:34
Репутация: -344969098
Откуда: Париж

Re: И ЕЩЕ РАЗ О ПРОЕКЦИЯХ. Теоретические вопросы.

Сообщение Boris » 19 окт 2011, 01:35

Ну, что ж я готов представить алгоритм действий и результаты своих вычислений в MapInfo для СК-63 на Литву(?), приведенную выше, поскольку она прописана стандартных проекциях Мапинфо (со смещением в метрах). В архиве есть код на языке MapBasic, что бы желающие могли проверить и/или повторить результат в чем-то более точном:
Но в начале код на языке командной строки Window:

Код: Выделить всё

for /l %i in (755000,5000,1745000) do @for /l %j in (5505000,5000, 6495000) do @echo %i; %j >> i_j.txt 
цель - создать набор координат точек в пределах квадрата в 1000 км = +-500 км от цетрального меридиана СК и 6000 км от экватора. Причина именно таких данных - а) чистый волюнтаризм б) мапинфо оперирует границами карты - это повышает точность позиционирования
2.7z
код на мапбасике, файл-таблица с точками, результирующие таблицы прилагаются:
(794.1 КБ) 374 скачивания
Краткий поряд действий:
  • создаем таблицы, в которых смещение по Y (мапинфо) или X (геодезический гаусс-крюггер) задано либо метрами, либо угловыми величинами. Тут есть момент, котрый я не понял, смещение в градусах задано как южное, т.е. в нормальной записи должно быть отрицательным и ниже экватора, ложное смещение задано то же отрицательным, обычно это в мапинфо означает сдвиг на север. Честно скажу - долго не думал и задал смещение в градусах как положительное. Точки в результате совпали так, что бы не было сомнений в тождественности. Хотя может просто удачно совпало?
  • создали набор точек с приращением в 5 км и отстуспом от краев таблицы в те же 5 км. Что бы не было заморочек на границе таблицы.
  • Дальше скопировали координаты X,Y в каждую из таблиц. Колонки 4-5
  • выполнили кодирование координат в точки (points). для каждой таблицы расчет координаты -> точки проводился при установленной ее родной СК
  • установили для расчетов СК от одной из таблиц
  • вычислили координаты кодированных точек точек в обеих таблицах в одной и той же СК. При тождественности СК, то какая из них установлена - должно быть без разницы.
  • В конце файла mb.txt (кода на мапбасике) выписаны два запроса, которые считают разницу между парными точками. Сперва как сумму квадратов, затем как расстояние.
А вот статистика всего файла расстояний:
res.jpg
res.jpg (14.78 КБ) 7374 просмотра
Мой итог как никакого геодезиста:
  • средняя погрешность 3.5 см. Не слишком большая, что бы заметить в любой работе. Может возникать как набранный остаток округления при вычислениях переходов из одной СК в другую? Может, но для этого великовата.
  • Может 3.5 см быть ошибкой точности самой целочисленной мапинфо? Не может. Для квадрата в 1000 км, погрешность дожна быть от 0.5Е-06 до 0.5Е-04.
  • где-то есть погрешность погрешность в 50 см. скорее всего на границах +-500 км. Не смертельно, но много. С другой стороны на широте Москвы - это уже 17 градусов по меридиану.

stout
Участник
Сообщения: 91
Зарегистрирован: 13 апр 2012, 17:09
Репутация: 52

Re: И ЕЩЕ РАЗ О ПРОЕКЦИЯХ. Теоретические вопросы.

Сообщение stout » 14 апр 2012, 22:46

По просторам интернета гуляет анонимно рифмованая притча:
Шесть седовласых мудрецов сошлись из разных стран.
К несчастью, каждый был незряч, зато умом блистал.

Они исследовать слона явились в Индостан.
Один погладил бок слона... Довольный тем с полна
Сказал он: Истинна теперь как божий день видна.
Предмет, что мы зовем слоном - отвесная стена!

Другой, наткнувшийся рукой на бивня остриё,
Сказал: Коль интерестно вам, вот мнение моё:
Ни что иное этот слон, как острое копьё!

А третий хобот в руки взял и закричал: Друзья!
Гораздо проще наш вопрос, уверен в этом я!
Сей слон - живое существо, а именно — змея!

Мудрец четвертый обхватил одну из ног слона
И важно молвил: Это ствол. Картина мне ясна!
Слон - дерево, что зацветет, когда придет весна!

Однако пятый мудрый муж был, как гора, высок.
И потому достать рукой слоновье ухо смог.
Слон, без сомненья - веер, - он медленно изрек.

Тем временем шестой из них добрался до хвоста.
И рассмеялся от того, как истина проста.
Ваш слон - веревка. Если ж нет - зашейте мне уста!

А как известно,мудрецам присущ упрямый нрав.
Спор развязав, они дошли едва ль не до расправ.
Но правды ни один не знал, хотя был в чём-то прав.

Когда религиозный спор заходит где-нибудь,
Как верит каждый, что лишь он на правильный стал путь,
И слишком от гордыни слеп, чтобы увидеть суть.


Вполне вероятно, что ответы уже найдены, как-никак уже больше 3 лет прошло.
Но, всё же…
Boris писал(а):Имеем в качестве глобуса - эллипсоид Красовского.
Globus (лат.) - шар. ;)
Boris писал(а):Что задает центральный меридиан - понятно. Это место где цилиндр прижимают к эллипсоиду.

2. И связанный вопрос, является ли фигурой отображения именно цилиндр, или это более сложная фигура, повторяющая поверхность эллипсоида?
На этот вопрос ответ очень простой: А нет никакого цилиндра.
Один известный американский популяризатор науки как-то цинично заметил - основная задача популязатора при написании коммерчески успешной книжки - добится такого эффекта, когда читатель считал бы, что понял суть вопроса. Т.е. задача не в том, чтобы он понял, а в том, чтобы он только думал, что понял!
По моему, причиной ваших заблуждений является целый сонм неверных иллюстраций, коими пестрят не только популярные издания, но и некоторые монографии и более или менее солидные учебники. В интернете с этим просто беда. Постоянно натыкаешся на фразы типа "Cylindrical Projections are the result of projecting the globe onto a cylinder."
Повеселил рисуночек из "ArcGIS® 9 Understanding Map Projections" с лампочкой в центре Земли и обёртывающим цилиндром.
Очень редко можно найти правильные толкования, например такое, как это:
"The Mercator projection is often illustrated as a projection from the center of a globe onto a cylinder wrapped around the globe. However, that cylindrical projection is not the Mercator projection. If you see them side by side, you will see they dont look the same. The cylindrical projection does not have the property that Mercator was seeking; the rhumb lines are not straight lines. They do have the general form in common: meridians of longitude and parallels of latitude are straight, parallel lines. Also, since the meridians dont get closer together towards the poles (as they do on a globe), that means the scale must increase as the latitude increases." (http://babylon.acad.cai.cam.ac.uk/stude ... pirate.htm)

Вы имеете право спросить - а как же быть со книгой Снайдера (JOHN P. SNYDER) Map Projections — A Working Manual? Ведь там тоже нарисован и цилиндер и конус? Но посмотрите внимательно на подпись к этим ресункам:
Figure 1.—Projection of the Earth onto the three major surfaces. In a few cases, projection is geometric, but in most cases the projection is mathematical to achieve certain features.

Идея с цилиндром, на который происходит проектирование, это всего лишь иллюстрация, пытающаяся помочь понять на интуитивном уровне, что происходит на самом деле. Посмотрите, как люди пытаются представить себе, например, четырехмерный куб - https://www.google.com/search?q=%D1%87% ... 79&bih=665

Есть очень интересный (и незаслуженно забытый) метод Dr. David E. Wallis (существуют его отчеты по сходной тематике http://ntrs.nasa.gov/search.jsp?N=0&Ntk ... 20matchall). К сожалению электронных версий этих работ нет, да и сайта автора к огромному сожалению уже нет. Так вот, он утверждал, что получил проекцию ТМ путем проектирования эллипсоида на единичный диск(я надеюсь никто не усомнится в том, что этот диск покоится на трех китах) в комплЕксной плоскости, а затем развернул(спроектировал) этот диск на плоскость. По поводу последнего (по времени) метода Dr. Wallis'a была довольно жесткая пикировка в почтовой рассылке Proj.4 ещё в 2003 году.

Лучше всего по поводу цилиндра сказал Морозов в своей книге "Курс сфероидической геодезии":
В геодезической и картографической литературе проекцию Гаусса — Крюгера часто называют поперечно-цилиндрической равноугольной проекцией эллипсоида, однако без пояснения условности такого названия. Больше того, исходя из этого названия, иногда пытаются объяснить получение этой проекции посредством изображения эллипсоида на цилиндре, касающемся осевого меридиана. Конечно, таким способом получить проекцию Гаусса — Крюгера невозможно ни для эллипсоида, ни для шара. Дело в том, что в математической картографии-рассматривается целая группа поперечно-цилиндрических проекций, куда входит и проекция Гаусса — Крюгера. Главным признаком, объединяющим эти проекции, является постоянство масштаба на прямолинейном изображении одного из меридианов эллипсоида. Название этой группы проекций связано с одной-единственной проекцией, обладающей этим же признаком, а именно перспективной поперечно-цилиндрической проекцией. Эту проекцию действительно можно получить геометрическим построением: спроектировать поверхность эллипсоида на касательный по меридиану цилиндр прямолинейными лучами, проведенными из центра эллипсоида. Конечно, для получения любого перспективного изображения безусловно необходимо установить одну поверхность в определенное положение относительно другой и выбрать центр перспективы. Но во всех других видах изображений взаимное положение поверхностей не играет никакой роли и не может помочь найти уравнения изображения. Поэтому и для получения проекции Гаусса — Крюгера нет никакой необходимости «обертывать» эллипсоид цилиндром. Кстати, такие довольно распространенные в картографической литературе термины, как «касательный или секущий цилиндр», «касательный или секущий конус», «параллель касания», «меридиан касания» и т. п., надо понимать как чисто условные и не связывать с ними геометрического представления для объяснения или получения проекции.
geologic писал(а):Вообще-то Меркаторская проекция - чисто цилиндрическая, формулы-расчеты это подтверждают.
Не могли бы пояснить свою мысль? И какие формулы это подтверждают?
Boris писал(а):И еще пара вопросов:
3. Смещение начала координат, связанное со смещением по оси X (восток), как говорится, оно используется исключительно для удобства записи координат, что бы не было отрицательных значений. Имеет ли оно еще какой-нибудь смысл?
4. И собственно, основной вопрос, имеющий уже практический интерес. Смещение по оси Y (север) и смещение начальной параллели. Являются ли эти величины связанными? Можно ли одно выразить одно через другое в виде константы в пределах выбранной "дольки" проекции?
Ну, здесь ответ однозначный - да, являются. Чтобы там gamm не думал или говорил.
Точное аналитическое выражение(одно из возможных), связывающие геодезическую широту(B) и приращение долготы (удаление от осевого мередиана L) с прямоуголными координатами X и Y в проекции ТМ выглядит довольно просто:
X+iY = k*a(1-e^2)*П(e^2,e,invLam(Lam(B)+iL)) +(X_false + i(Y_false - k*M(B_origin)));
где k - масштабный коэффициент на осевом меридиане(=0.9996 для картографической системы координат UTM), а - большая полуось, е -эксцентриситет, В - геодезическая широта, L - удаление от осевого меридиана,i - мнимая единица; Lam - Ламбертиан (термин введен Lee, по аналогии с Гудерманианом):
Lam(e,B)=Arth(sin(B))-e*Arth(e*sin(B))),
invLam - функция обратная к Ламбертиану,
П - эллиптический интеграл 3-го рода в нормальной форме Лежандра (в данном частном случае его можно выразить через эллиптический интеграл второго рода плюс небольшой довесок из комбинации тригонометрических функций), B_origin - начальная широта, M(B_origin) - длина дуги меридиана от экватора до B_origin.

Если не верите мне, смотрите Map Projections — A Working Manual стр. 61 формула (8-10)
y = k0*{M-Mo+…}
… M is the true distance along the central meridian from the Equator to <B>.
… M0 = M calculated for <B_origin>, the latitude crossing the central meridian at the origin of the x, y coordinates.
Или раздел 5.9.14 Transverse Mercator CS specification; Table 5.20 — Transverse Mercator CS
в ISO/IEC 18026:2009(E) Information technology — Spatial Reference Model (SRM)
http://standards.sedris.org/18026_Ed2/t ... ACT_CS.PDF
там написано следующие (как и в предыдущем случае, я слегка изменил обозначения, не меняя сути)
N = N_false + k(a*V(B,L) - S(B_origin));
Вообще, весь этот стандарт заслуживает более пристального изучения. Очень полезный документ для тех, кто хочет разобраться в хитросплетениях систем координат.
Каминский Вадим писал(а):Поправка, в проекции ГК оси имеют обратное название чем принято в Декартовой СК, т.е. Х-север, Y-восток.
Сразу видно, что мы в разное время учились. Ещё в школе, сначала на уроках геометрии, а затем и физики нам рассказывали о правых и левых декартовых системах координат.
Boris писал(а):Кроме того, что я оси неправильно именовал. Это я признаю. Тем вы раскрыли мне давно мучившую меня тайну, почему у всех на север идет Y, а у геодезистов X. Теперь я понял - это проверка "свой-чужой".
Однако! Как споро вы Иоганна Генриха Луиса Крюгера в "свои"(или "чужие"?) записали. Потому, как в его фундаментальной работе "Konforme abbildung des erdellipsoids in der ebene" на странице под номером 3 в правом нижнем углу именно левая система координат изображена, и оси подписаны - Х - вверх, Y - вправо, на восток.
Boris писал(а):Дорогие коллеги,
ну ответьте все же мне, если "начальная параллель" НЕ входит в число этих самых параметров, то почему EPSG и ArcGIS, на областях (не субъектах РФ, а областях действия), приводят это самое смещение, а False Norting не приводят. Весь вопрос в том и соcтоит, сколько параметров надо учитывать в своей эмпирической работе.
Почему EPSG и ArcGIS используют B_origin - это вопрос лучше задавать им самим. Рискну выдвинуть рабочую гипотезу. Зная, что в наших "Х-файлах" фигурируют N(E)false, прекрасно понимая взаимно однозначное соответствие N_false и S(B_origin) и не желая подставлятся под наши "доблестные органы" используют те параметры, которые не прописаны в "Х-файлах".
Boris писал(а):И какова погрешность, если вместо False Norting берешь смещение параллели, или наоборот.
Теоретически, должна равняться нулю.

si jeunesse savait, si vieillesse pouvait.
best regards,
Alexander Starostin aka stout

Донецков
Гуру
Сообщения: 3058
Зарегистрирован: 19 май 2010, 19:44
Репутация: 189

Re: И ЕЩЕ РАЗ О ПРОЕКЦИЯХ. Теоретические вопросы.

Сообщение Донецков » 30 апр 2012, 17:58

Тогда вопрос для сведущих - кто видел "Х-файлы", в них N_false и S(B_origin) приведено?

Petruxin
Гуру
Сообщения: 1695
Зарегистрирован: 14 июн 2011, 16:47
Репутация: 133
Ваше звание: Завсегдатай
Откуда: Череповец

Re: И ЕЩЕ РАЗ О ПРОЕКЦИЯХ. Теоретические вопросы.

Сообщение Petruxin » 30 апр 2012, 20:25

Предлагаю это сообщение засчитать за 7-8 :shock:

Boris
Гуру
Сообщения: 4231
Зарегистрирован: 10 апр 2006, 22:34
Репутация: -344969098
Откуда: Париж

Re: И ЕЩЕ РАЗ О ПРОЕКЦИЯХ. Теоретические вопросы.

Сообщение Boris » 30 апр 2012, 20:39

"я не видел, но скажу!" :wink:
раз меня цитировали, то не вдаваясь в весь ответ, замечу, что EPSG и ESRI НЕ приводят параметров на территории, где они составляют какую-то там тайну. Они их приводят только на грузию, прибалтику, м.б. казахстан. К стати на прибалтику не только они одни. Сами прибалты на каком-то виденном мною эстонском сайте приводят смещение в градусах по широте. С чего бы это им опасаться каких-то мифических проблем на территории не входящие в СНГ? И очевидно, что эти параметры официально утекли из государственных органов свободных (теперь) стран.
Не очень понял про какие именно Х-файлы мы говорим. Про каталоги координат? Там вроде только цифры в столбик должны быть. Про руководство по введению или установлению СК-63? Официально заявлено, что их все про... (потеряли).
А в целом, в интернете я видел и так и так. Мой опыт говорит, что для мапинфо при разном задании погрешность - есть. Не существенная, но есть. Теперь бы вспомнить какая :) Помню, что очень несущественная. Если это действительно кого-то еще волнует - я создам MIF с шагом в 5 км на север в ширину соответствующей зоны. 3-х градусной или 6-ти это как скажите.

Донецков
Гуру
Сообщения: 3058
Зарегистрирован: 19 май 2010, 19:44
Репутация: 189

Re: И ЕЩЕ РАЗ О ПРОЕКЦИЯХ. Теоретические вопросы.

Сообщение Донецков » 01 май 2012, 10:10

В том то и вопрос, что у "соседей": Прибалтика (блок C), Закавказье (блок А), Казахстан (блок К) - приводятся градусы, все доморощенные российские и украинские программы (особенно интересны созданные в конце 80-х - начале 90-х, т.к. их явно создавали "знающие" люди без "иностранного" влияния, т.к. его еще просто не было) - используют смещение метры... На моем опыте расхождения в том и другом случае не превышает нескольких сантиметров. Как-то было дело, в АГП заказчики заказали пересчет имеющихся каталогов из СК-42 в СК-63, получили и стали работать и нам передали по нашей части объектов, а потом пыяснилось, что в АГП схалтурили, нет чтобы спуститься в хранилище и выписать каталожные данные, в чем-то сами пересчитали, расхождения составляли до нескольких метров...

Ответить

Вернуться в «Координаты и привязка»

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 2 гостя