annageo Можно, ну нужно иметь минимум 3-4 точки в обоих системах координат. Если этого нет, то нужно будет определять в МСК-50 от ближайших пунктов ГГС или ОМС. От точности определения координат в двух системах будет зависеть и точность пересчета.
Я бы перевел из МГГТ в 63г. (т.к. у меня есть программка перевода) пару точек, потом в МСК-50. Забил бы эти точки в Акад и посадил бы на них весь массив.
есть такие точки, но как перевести все сразу? пробовала в 63 потом в мск-50... переводимый участок ложится примерно (естественно с расхождением) а остальной массив очень сильно улетает
На самом деле пересчет из МГГТ в МСК50 (СК42,СК63) и обратно не составляет труда. Если задача еще актуальна могу рассказать как оно делается. Или просто пересчитать чего надо.
Для начала нужно понять с чем имеем дело. А дело мы имеем, в случае с МГГТ, с проекцией Гаусса-Крюгера на эллипсоид Бесселя 1841, и случае с МСК50, СК42, СК63 с той же проекцией на эллипсоид Красовского. Начало координат МГГТ широта 55º40', долгота 37º30'. Значит при вычислении плоских координат осевой меридиан уже известен 37º30' ну и дельта по Х легко вычисляется по формулам Бугаевкого, это -6170810.1136 м. С минусом потому, что эллипсоид Бесселя -африканский, т.е. плюсы у него в южном полушарии. С СК42 все еще проще осевой меридиан на Москву 39º дельта по Y стандартная -7500000.0000. МСК50 и СК63 рассматривать пока не станем, поскольку их параметры проекций законспирированы, но если Вы работаете с этими СК, они должны быть Вам известны. Теперь при желании можно получить прямоугольные координаты и для МГГТ и для СК42, и имея не менее 3-х общих точек в той и иной системе вычислить 7 элементов перехода от одной системы к другой. Можно по формулам аффинного преобразования, можно по формулам (23) из ГОСТ Р 51794-2001. Однако, в жизни все еще проще. Народ создал Leica geomatic office, Trimble geomatic office, и много другого софта для обработки GPS. Расскажу на примере Leica, просто я сам именно с ее софтом все это и проделываю, но аналогичные возможности существуют и у других производителей GPS. Например, нужно пересчитать координаты из МГГТ в СК42. Закачиваем координаты МГГТ в проект, создаем проекцию с ранее названными параметрами, а именно False Northing -6170810.1136, Central Meridian 37º30' и Zone Width (ширина зоны)3º, на элипсоиде Бесселя 1841. Присваиваем проекцию проекту, софт сразу выдает нам координаты в WGS84. Экспортируем полученные координаты WGS84 в отдельный файл. Далее, в другой проект импортируем координаты в СК42 и соответственно создаем проекцию False Easting -7500000.0000,Central Meridian 39º и Zone Width 6º, на эллипсоиде Красовского. Соответственно назначаем проекцию проекту. Далее создаем еще один проект, в который импортируем координаты WGS84 из ранее созданного отдельного файла. А потом совсем все просто, в модуле Datum Map (Leica geomatic office) выбираем проект с WGS84 как трансформируемый, а проект с СК42 как опорный, по общим точкам, в автомате, вычисляются параметры перехода между эллипсоидами и соответственно перевычисляются координаты определяемых точек. И все. Точность перевычислений в пределах 5 мм. Я использую около 40 точек со всей Москвы, координаты которых известны в обеих системах. С МСК50 и СК63 все аналогично, только вместо параметров СК42 используются параметры этих проекций. (Добавление) Да кстати, пересчитать координаты из СК63 в МСК50 или СК42 или в любом другом направлении, можно используя всего лишь EXCEL, просто опять требуется знать параметры проекций исходной и конечной систем координат.
Вижу, что некоторые выражают сдержанное недоверие. И тем не менее, так оно все и есть. Это все кажется сложным до тех пор, пока сам не увидишь и не попробуешь. А затем становится непонятным, как из-за такой ерунды вообще могли возникать какие-то проблемы.
::gramercy:: ::shock1.gif:: :o Аффтар, йад не пей, в бабруйск ни еть. лучше пешы исчо.© Надеюсь, после моего эмоционального всплеска вы больше не станете писать о сдержанном недоверии. З.Ы. Что есть формулы Бугаевcкого в вашем понимании?
Л.М. Бугаевский (профессор) "Математическая картография" Москва "Златоуст" 1998 г. стр. 246. Учите мат. часть, и желательно без эмоций.
2 Yuriy Честно говоря, я надеялся что вы всё-таки откажитесь от более чем странного объяснения с привлечением "африканского" эллипсоида Бесселя. Потому что всё гораздо проще. Смещение - длина дуги меридиана от экватора до параллели 55°40'. Я, также как и вы, вычислил эти смещение. Но тремя разными способами. 1. Координатный калькулятор Pinnacle. Получилось: -6170810.069609 2. По формулам Крюгера, которые вы по недоразумению назвали формулами Бугаевского. ( L. Krueger,Konforme Abbildung des Erdellipsoids in der Ebene, Royal Prussian Geodetic Institute, New Series 52,172 pp. 1912 см. рисунок внизу). Точность этих формул чуть выше, чем тех, которые приведены в книге Морозова "Курс сферической геодезии" и, соответственно, "Математическая картография". Если воспользоваться какой-нибудь системой символьных вычислений (мне давно нравиться MuPAD, сейчас это составная часть MATLAB), то можно уточнить выражения: R = a/(n + 1) * (1 + 1/4*n^2 + 1/64*n^4 + 1/256*n^6 + 25/16384*n^8+…) a2 = 1/2*n - 2/3*n^2 + 5/16*n^3 + 41/180*n^4 - 127/288*n^5 + 7891/37800*n^6 + 72161/387072*n^7 - 18975107/50803200*n^8 … a4 = 13/48*n^2 - 3/5*n^3 + 557/1440*n^4 + 281/630*n^5 - 1983433/1935360*n^6 + 13769/28800*n^7 + 148003883/174182400*n^8 a6 = 61/240*n^3 - 103/140*n^4 + 15061/26880*n^5 + 167603/181440*n^6 - 67102379/29030400*n^7 + 79682431/79833600*n^8 a8 = 49561/161280*n^4 - 179/168*n^5 + 6601661/7257600*n^6 + 97445/49896*n^7 - 40176129013/7664025600*n^8 a10 = 34729/80640*n^5 - 3418889/1995840*n^6 + 14644087/9123840*n^7 + 2605413599/622702080*n^8 a12 = 212378941/319334400*n^6 - 30705481/10378368*n^7 + 175214326799/58118860800*n^8 a14 = 1522256789/1383782400*n^7 - 16759934899/3113510400*n^8 a16 = 1424729850961/743921418240*n^8 Получилось: -6170810.06964166 3. Вычисление длины дуги, как эллиптического интеграла 2 рода (формулы и эллиптические интегралы Карлсона) Получилось: -6170810.06964165 Для этих вычислений я использовал следующие параметры эллипсоида Бесселя 1841г a = 6377397.155 1/f = 299.1528128 Это не единственные "официальные" значения для эллипсоида Бесселя. Дело в том, что Бессель определял не сжатие, а малую полуось. Но и с ними результат не совпадает с вашим.
Я и не утверждал, что самый умный. Насколько было понятно из предыдущих сообщений, народ просто был не в теме. И к пересчетам координат относился как к священной корове. А это совсем не так. Если речь лично обо мне, то мне нравится Delpfi, но справляюсь и Visual Basic(ом). Насчет авторства формул и точности вычислений, то попались они мне у Бугаевского, я и обозвал их Бугаевскими, хотя они наверняка встречаются и у Красовского и иных авторов, абсолютная точность при определении смещения начала координат нужна если пересчитывать координаты непосредственно, без опорных точек, а в предложенном мной варианте свободный WGS84 блок, трансформируется по опорным точкам, и в этом случае неприемлемая ошибка за положение блока на эллипсоиде возникает если смещение превышает 20 м, и в эти 20 м, как Вы очевидно заметили, я попал. Кстати, итоговые невязки в действительности не превышают 2 мм, насчет 5 мм я просто подстраховался, случаи разные бывают. Насчет африканскости эллипсоида, акцентировать на этом внимание пришлось из за особенностей Лейковского софта, например Ски Про просто выдает чушь если забыть про африку, а Лейка Геоматик Оффисе уже нормально отрабатывает зеркалку. Да, опять к формулам, я заметил, что в тот полиномчик не помешал бы еще один членчик, так что не исключено, что Ваш результат точнее моего. Но руководствуясь вышеприведенными соображениями про 20 м уточнять результат не стал.
2. По формулам Крюгера, которые вы по недоразумению назвали формулами Бугаевского. ( L. Krueger,Konforme Abbildung des Erdellipsoids in der Ebene, Royal Prussian Geodetic Institute, New Series 52,172 pp. 1912 см. рисунок внизу). Точность этих формул чуть выше, чем тех, которые приведены в книге Морозова "Курс сферической геодезии" и, соответственно, "Математическая картография". Пришлось делать програмульку по непосредственному пересчету МСК50-МГГТ-МСК50, т.к. тяжко народу со SKI морочиться, и попутно выяснилось, что у Бугаевского все же иные формулы, нежли у Крюгера с Морозовым. И поскольку ссылок на иное авторство не обнаружено, то видать они все же Бугаевские. Однако грубоваты. Но так и назначение у них иное.
Yuriy, вы хороший человек, но истина мне дороже. Поэтому ещё раз, в разрядку, ф о р м у л ы____о д и н а к о в ы е. Если считаете иначе, то покажите, в чем видите разницу. Фраза имеет право на существование, только в том случае, если подтверждается цифрой. Иначе это пустое сотрясание воздуха. По моим прикидкам на пальцах, ошибка формул, приведённых в книге Льва Моисеевича не превышает сотых долей миллиметра для стандартной зоны в 6°. Так как похоже, что тема вам интересна, то посмотрите ссылки: A HIGHLY ACCURATE WORLD WIDE ALGORITHM FOR THE TRANSVERSE MERCATOR MAPPING - обратите внимание на то, что для суммирования рекомендутся использовать алгоритм Кленшоу. И это в высшей степени правильно! Великолепная библиотечка, я уже давал на неё ссылку - GeographicLib is a small set of C++ classes for performing conversions between geographic, UTM, UPS, MGRS, geocentric, and local cartesian coordinates Успехов в написании программулек! Пусть их будет много и разных.
На самом деле, все видится немного не так. Считать можно как угодно, но результаты будут скорее всего неправильными. Есть официальные пораметры связи МСК с СК-42 и СК-42 с СК-63. Этими параметрами пользуются ТИГГН, АГП и другие. Мы их не знаем, нам их не открывают. А вот параметров связи МСК с СК-63 нет. Их, конечно, можно вычислить, но на каждом участке они получаются разными и координаты намного отличаются от официальных МГГТ. В этом-то и проблема. Ну, а связь СК-63, как региональной в МО до 2009 г. с МСК-50 совершенно официально известна и на этом сайте неоднократно публиковалась. Т.е. обязательно надо ччитать через 42-й год. Есть совершенно "правильные" DOS программы, но они не позволяют пакетно пересчитывать координаты. А те, что позволяют, ужасно критичны к кодировке и мороки с ними много. А официальной АРМ ТИГГН Вам скорее всего не найти.
Уважаемый, вы бы такое жирное IMHO (In My Humble Opinion - по моему скромному мнению) к своему посту добавили, дабы не вводить людей в заблуждение. З.Ы. Если параметров мы не знаем, то хотя бы перечислите, что является или может являться параметрами связи. Зело интересно взглянуть, хоть одним глазком, именно на параметры связи, а не на параметры самих систем координат.
