Какой азимут был исходным не так принципиально, как тот факт, что геодезический азимут Саблино-Бугры в СК-42 равен 317°02'54.29", а в СК-32 был равен 317°02'50.62" (это все открытые "публикабельные" данные). В статье "Саблинская базисная сеть" был вычислен геодезический азимут Бугры-Саблино в WGS84 равный 136°54'48,0", его надо только развернуть. --- Сообщения объединены, 14 янв 2019, Оригинальное время сообщения: 14 янв 2019 ---
Короче говоря, надо пробовать пересчитывать высоты по тем параметрам, что получились и, смотреть разницу. Благо данных для контроля теперь стало больше. Может какую-то зависимость можно будет увидеть. Потом попробовать с поэкспериментировать с Сигналом А. А, как Вы думаете, сильно ли разворот в 4" может повлиять на результат? Просто интересно мнение, как специалиста.
Историческое значение высоты на довоенной топокарте равно 74.9 метра и это дожна была быть высота "старого сигнал А" (новый восстановлен в 1981 году). Я сейчас не помню его высоту в саженях 19 века :) Астростолб это новодел 1973 года. --- Сообщения объединены, 14 янв 2019, Оригинальное время сообщения: 14 янв 2019 --- У меня есть такое мнение что СК-32 должна прекрасно соответствовать всем современным требованиям для кадастровых работ в Ленинградской, Новгородской и Псковской областях. Так как граница с Норвегией (и северный кусок с Финляндией) была (до самого последнего времени) демаркирована в СК-32, то можно оценить ее точность в Карелии и Мурманской области. Для Чехии надо считать.
Уважаемые cr2, andrew_klikunov ! Коль скоро здесь обсуждается некий "инструмент" для отождествления геодезической высоты и нормальной, неужели нельзя обойтись без азимутов исходных дат. А учебники говорят вот о чем: при установлении ИГД СК-42 эллипсоид был выбран из нескольких условий, одно из которых равенство нормального и возмущающего потенциалов в точке счета высот. Разве имеет богатый смысл исходная отметка? По-моему смысл находится в точной передаче этого значения дальше. пф.. и здесь кувырк, значения ошибок растут, не только на восток (но с прогрессией), сколько на северо-восток.
Уважаемые stout и Yuri V., не могли бы Вы разъяснить, для чего в WGS 84 используются каждая из двух разных GM? Для GRS-80 всего-то одна... Просто есть любопытные на мой взгляд результаты по расчету. Я их обязательно приложу, но в начале хотелось бы узнать по GM.
Объяснение содержится в основополагающем документе: Department of Defense World Geodetic System 1984 Its Definition and Relationships with Local Geodetic Systems "3.2.3.2 Special Considerations for GPS Based on a recommendation in a DMA letter to the Air Force [21], the refined WGS 84 GM value (3986004.418 x 10 8 m 3 /s 2 ) was implemented in the GPS Operational Control Segment (OCS) during the fall of 1994. This improvement removed a 1.3 meter radial bias from the OCS orbit estimates. The process that generates the predicted broadcast navigation messages in the OCS also uses a GM value to create the quasi-Keplerian elements from the predicted Cartesian state vectors. The broadcast elements are then interpolated by a GPS receiver to obtain the satellite position at a given epoch. To avoid any loss of accuracy, the GPS receiver’s interpolation process must use the same GM value that was used to generate the fitted parameters of the broadcast message. Note that this fitting process is somewhat arbitrary but must be commensurate with the algorithm in the receiver. Because there are many thousands of GPS receivers in use around the world and because proposed, coordinated software modifications to these receivers would be a costly, unmanageable endeavor, Aerospace Corporation [25] suggested that the original WGS 84 GM value (3986005.0 x 10 8 m 3 /s 2 ) be retained in GPS receivers and in the OCS process which fits a set of broadcast elements to the Cartesian vectors. This approach takes advantage of the improved orbit accuracy for both the estimated and predicted states facilitated by the refined GM value and avoids the expense of software modifications to all GPS receivers. For the above reasons, the GPS interface control document (ICD-GPS-200) which defines the space segment to user segment interface should retain the original WGS 84 GM value. The refined WGS 84 GM value should continue to be used in the OCS orbit estimation process. Most importantly, this approach avoids the introduction of any error to a GPS user. "
Как я понял, значение 3986004.418 x 10 8 m 3 /s 2 используют для высокоточных определений орбит. Это правильно? Но, и в программе h_synth EGM2008 оно тоже используется. А вот зачем нужно значение 3986005.0 x 10 8 m 3 /s 2 мне так и непонятно. Оффтоп (Move your mouse to the spoiler area to reveal the content) Пытался перевести через google, но там как-то вообще невнятно получилось
Грубо говоря, уже были тысячи приёмников, которые использовали для расчёта координат верха старое значение геогравитационной постоянной. Введение нового значения превратило бы их в кирпич. Странно, у меня google вполне понятно перевёл (если не считать сокращений) "3.2.3.2 Особые соображения для GPS, основанные на рекомендации в письме dma военно-воздушным силам [21], уточненное значение 84 мг / кг (3986004,418 x 10 8 м 3 / с 2) было реализовано в сегменте оперативного управления GPS (ocs). ) осенью 1994 г. это усовершенствование устранило радиальное смещение в 1,3 метра из оценок орбиты окс. процесс, который генерирует предсказанные широковещательные навигационные сообщения в ocs, также использует значение gm для создания квазипеплеровских элементов из предсказанных векторов декартовых состояний. элементы широковещания затем интерполируются GPS-приемником для получения положения спутника в заданную эпоху. во избежание потери точности процесс интерполяции GPS-приемника должен использовать то же значение gm, которое использовалось для генерации подогнанных параметров широковещательного сообщения. обратите внимание, что этот процесс подбора несколько произвольный, но должен быть соизмерим с алгоритмом в приемнике. Поскольку во всем мире используется много тысяч GPS-приемников и поскольку предложенные согласованные модификации программного обеспечения этих приемников были бы дорогостоящим, неуправляемым предприятием, аэрокосмическая корпорация [25] предположила, что исходное значение wgs 84 г (3986005,0 x 10 8 м). 3 / s 2) сохраняться в приемниках gps и в процессе ocs, который подгоняет набор широковещательных элементов к декартовым векторам. этот подход использует преимущества улучшенной точности орбиты как для расчетного, так и для прогнозируемого состояний, чему способствует уточненное значение gm, и позволяет избежать затрат на модификацию программного обеспечения для всех GPS-приемников. по вышеуказанным причинам в документе управления интерфейсом gps (icd-gps-200), в котором определяется интерфейс космического сегмента к пользовательскому сегменту, должно сохраняться исходное значение wgs 84 г. уточненное значение 84 г / кг в кг следует продолжать использовать в процессе оценки орбиты окс. самое главное, этот подход позволяет избежать появления любой ошибки пользователю GPS. " Не совсем. Должно использоваться то значение, которое используется при описании гравитационного поля, будь то разложение по сферическим функциям или модели точечных масс. Вспомнился один забавный эпизод при выводе ПЗ-90. На вход программы интегрирования уравнений движения ИСЗ величина GM должна была передаваться в размерности км3/с2. Один товарищ, из вполне здравых рассуждений, что работаем в системе СИ передал значение в м3/с2. Как следствие, спутник летал под поверхностью Земли на глубине в 1000 км. Невязки O-C (observantum minus calculantum) были "прекрасны"
Оффтоп (Move your mouse to the spoiler area to reveal the content) почему не понятно, гугл нормально перевёл: «3.2.3.2 Особые соображения для GPS На основе рекомендации в письме DMA ВВС [21] уточненное значение GM WGS 84 (3986004,418 x 10 8 м 3 / с 2) было реализовано в сегменте оперативного управления GPS ( OCS) осенью 1994 г. Это усовершенствование устранило радиальное смещение 1,3 м из оценок орбиты OCS. Процесс, который генерирует прогнозируемые широковещательные навигационные сообщения в OCS, также использует значение GM для создания квазикеплеровых элементов из прогнозируемого декартова векторы состояния. Затем элементы широковещания интерполируются приемником GPS для получения положения спутника в данной эпохе. Чтобы избежать потери точности, процесс интерполяции приемника GPS должен использовать то же значение GM, которое использовалось для генерации подходящих параметров широковещательное сообщение. Обратите внимание, что этот процесс подбора несколько произвольный, но должен быть соизмерим с алгоритм в приемнике. Поскольку в мире используется много тысяч GPS-приемников и поскольку предложенные согласованные модификации программного обеспечения этих приемников будут дорогостоящим, неуправляемым мероприятием, Аэрокосмическая корпорация [25] предположила, что первоначальный WGS 84 значение GM (3986005,0 x 10 8 м 3 / с 2) должно сохраняться в приемниках GPS и в процессе OCS, который подгоняет набор широковещательных элементов к декартовым векторам. Этот подход использует преимущества улучшенной точности орбиты как для расчетного, так и для прогнозируемого состояний, чему способствует уточненное значение GM, и позволяет избежать затрат на модификацию программного обеспечения для всех приемников GPS. По вышеуказанным причинам управление интерфейсом GPS документ (ICD-GPS-200), в котором определяется интерфейс космического сегмента к пользовательскому сегменту, должен сохранять первоначальное значение WGS 84 GM. Уточненное значение GM в WGS 84 следует продолжать использовать в процессе оценки орбиты OCS. Что наиболее важно, этот подход позволяет избежать появления ошибок для пользователя GPS." --- Сообщения объединены, 23 янв 2019, Оригинальное время сообщения: 23 янв 2019 --- Оффтоп (Move your mouse to the spoiler area to reveal the content) Уважаемый andrew_klikunov, вы слишком преувеличиваете мой авторитет в данных вопросах, а, спрашивая меня наравне со stout так обоготворяете. Я даже не знал о том, что показал в цитате stout - оказывается начальную эпоху траектории ИСЗ можно задать не только временем, но и гравит. постоянной. Но это цветочки. Ягодки в самом смысле произведения GM, размерность км3/с2 - ускорение объёма без массы . Этим ученые ушли от проблемы точного определения массы Земли, которая необходима для составления уравнений движения ИСЗ, по закону Всемирного тяготения. Короче (если дилетанту позволительно что-то укорачивать), это такой коэффициент/константа, регулирующая метрический (в декартовых видимо) масштаб описания ГПЗ. Итеративно, путем экспериментов, подбирается. Теперь рассказывайте вы.
Точно в яблочко! Один из возможных физических способов задания линейного масштаба. Другой – скорость света.
Выкладываю результаты. В таблицах приведены пересчет с GRS-80 и WGS 84 (маркировал цветами). points - сводная табличка, вторая табличка - с полным расчетом (точки приведены в том же порядке, что и в сводной).
Вот две интересных ссылочки https://www.svf.stuba.sk/en/departm...ence-and-research/downloads.html?page_id=4996 http://imi.cas.sc.edu/gravimetric-quantities/ Статьи с описанием Soft'a есть на Генезисе. Программки по второй ссылке зело охочие до памяти да и подготовительный этап неспешный. Но считает именно аномалию высоты, а не высоту геоида. И быстро считает ============================================================================ Подготовительный этап на моём компе HTML: Microsoft Windows [Version 10.0.17134.523] (c) Корпорация Майкрософт (Microsoft Corporation), 2018. Все права защищены. [font=courier new]E:\WORK\hsynth_source>hsynth_init.exe Gravimetric quantity? (T, zeta, Dg, dg, xi, eta, Trr): zeta number of rows; number of columns; number of ellipsoids 4040 8053 38 ellipsoid number; height above base line; norm on ellipsoid 1 0.0000000000000000 106.50031083776784 2 0.33550536649285584 106.50027846944580 3 5.3680901151875524 106.49979296257357 4 27.176049440019522 106.49768948772272 5 85.890529840350638 106.49202941857754 6 209.69776048845213 106.48010968857686 7 434.84467696367381 106.45848613311860 8 805.65043924020995 106.42302027652319 9 1374.5270658704494 106.36895977425463 10 2202.0131338176916 106.29215881966186 11 3356.8252343634376 106.18673315142873 12 4915.9326381545452 106.04666152874260 13 6964.6614182648018 105.86623013783019 14 9596.8351239046951 105.63990092143938 15 12914.960008583430 105.36433239835331 16 17030.463819035613 105.03474708958333 17 22063.998275486876 104.64631090239089 18 28145.816656867741 104.19580416401541 19 35416.239392223004 103.67829678148465 20 44026.222308038967 103.09256851906240 21 54138.044259390241 102.43417169984326 22 65926.133365949383 101.70229346717974 23 79578.054086928110 100.89530191597001 24 95295.681033063869 100.01374407470124 25 113296.58989219430 99.057174782105363 26 133815.70134400437 98.026694274434348 27 157107.22062007335 96.920553633108852 28 183446.92376069381 95.733389954994863 29 213134.85205786114 94.456340733640843 30 246498.48920829478 93.079710376308739 31 283896.51205292338 91.597110389964087 32 325723.22639739211 90.005862265621573 33 372413.82554462133 88.305308166881829 34 424450.64249549306 86.496034868945642 35 482370.60953074053 84.579236527274489 36 546774.19411615888 82.556494342869129 37 618336.15192932333 80.429992800451785 38 697818.53202118853 78.202865475063788 --------------------------------------------------------------- program hsynth_init --------------------------------------------------------------- Gravimetric quantity: height anomaly in meters Total time = 2338.56615 CPU seconds Coefficients load time = 9.28370 CPU seconds Coefficients formation time = 1.44343 CPU seconds Synthesis time = 2317.08255 CPU seconds Values write time = 9.68915 CPU seconds --------------------------------------------------------------- E:\WORK\hsynth_source>hsynth_init.exe Gravimetric quantity? (T, zeta, Dg, dg, xi, eta, Trr): T number of rows; number of columns; number of ellipsoids 4040 8053 38 ellipsoid number; height above base line; norm on ellipsoid 1 0.0000000000000000 1041.6852537336774 2 0.33550536649285584 1041.6848271681149 3 5.3680901151875524 1041.6784288686742 4 27.176049440019522 1041.6507068904837 5 85.890529840350638 1041.5761029546120 6 209.69776048845213 1041.4189485209313 7 434.84467696367381 1041.1337045231355 8 805.65043924020995 1040.6654421629637 9 1374.5270658704494 1039.9506704995863 10 2202.0131338176916 1038.9293306036463 11 3356.8252343634376 1037.5219602150917 12 4915.9326381545452 1035.6454865345693 13 6964.6614182648018 1033.2177379650727 14 9596.8351239046951 1030.1563675447517 15 12914.960008583430 1026.3981067201403 16 17030.463819035613 1021.8663715929290 17 22063.998275486876 1016.4809451461633 18 28145.816656867741 1010.1773702723312 19 35416.239392223004 1002.8737343167808 20 44026.222308038967 994.52636339918161 21 54138.044259390241 985.05848024122110 22 65926.133365949383 974.43173018325774 23 79578.054086928110 962.60102210282639 24 95295.681033063869 949.54446490777150 25 113296.58989219430 935.23295947149973 26 133815.70134400437 919.65673772325397 27 157107.22062007335 902.78275164348634 28 183446.92376069381 884.54991817135317 29 213134.85205786114 864.87253866652691 30 246498.48920829478 843.66731239377077 31 283896.51205292338 820.89345991569871 32 325723.22639739211 796.55330283574892 33 372413.82554462133 770.68028623997293 34 424450.64249549306 743.32971011652728 35 482370.60953074053 714.57530087128350 36 546774.19411615888 684.50579036218619 37 618336.15192932333 653.22882371548565 38 697818.53202118853 620.87130401768889 --------------------------------------------------------------- program hsynth_init --------------------------------------------------------------- Gravimetric quantity: disturbing potential T in m^2/s^2 Total time = 2361.72216 CPU seconds Coefficients load time = 9.28566 CPU seconds Coefficients formation time = 1.42114 CPU seconds Synthesis time = 2341.07810 CPU seconds Values write time = 8.78912 CPU seconds --------------------------------------------------------------- E:\WORK\hsynth_source>[/font]
Не совсем понятен Ваш замысел. Вы хотите, чтобы я посчитал аномалию высоты с помощью этого ПО и потом сделал пересчет на элл. Красовского? P.S.: забыл написать, что в тех таблицах использованы значения высоты геоида из источников, приведенных выше (в статье из журнала), т.е. я не использовал EGM2008.
Началось всё с того, что нашлись контрольные данные по высотам квазигеоида над эллипсоидом Красовского. Они приведены тут. Это территория бывшей Чехословакии. Там, так же приведены и источники. Формулы приведены тут. Результаты ( в т.ч. и дельты по осям, которые использовались в формуле расчета разницы) приведены здесь. На эти данные я и опирался. Вроде ничего не упустил.
Вы имеете ввиду пересчет по ГОСТу 2008? Я тут исходил вот из чего: для расчета силы тяжести над эллипсоидом +- 50м по широте - не критично. Кстати ГОСТ выдал +- 2м по высоте.
stout, Вы не могли бы поделиться exe-шниками? А то у меня не получается скомпилировать фортрановские файлы.
У меня hsynth правильно работает только в версии на sparc/solaris (полное совпадение тестов с NGA). Сейчас пытаюсь разобраться с порядком следования коэффициентов в бинарном файле GAO98.jff для pinnacle. Нет ли у кого завалявшегося файла EGM96.jff для сравнения ?
