ГАО2012

Какой азимут был исходным не так принципиально, как тот факт, что геодезический азимут Саблино-Бугры в СК-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 мне так и непонятно.

Оффтоп

Пытался перевести через 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 должна была передаваться в размерности км³/с². Один товарищ, из вполне здравых рассуждений, что работаем в системе СИ передал значение в м³/с². Как следствие, спутник летал под поверхностью Земли на глубине в 1000 км. Невязки O-C (observantum minus calculantum) были "прекрасны"

 

Оффтоп

почему не понятно, гугл нормально перевёл:
«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 ---

Оффтоп

Уважаемый andrew_klikunov, вы слишком преувеличиваете мой авторитет в данных вопросах, а, спрашивая меня наравне со stout так обоготворяете.
Я даже не знал о том, что показал в цитате stout - оказывается начальную эпоху траектории ИСЗ можно задать не только временем, но и гравит. постоянной. Но это цветочки. Ягодки в самом смысле произведения GM, размерность км³/с² - ускорение объёма без массы . Этим ученые ушли от проблемы точного определения массы Земли, которая необходима для составления уравнений движения ИСЗ, по закону Всемирного тяготения. Короче (если дилетанту позволительно что-то укорачивать), это такой коэффициент/константа, регулирующая метрический (в декартовых видимо) масштаб описания ГПЗ. Итеративно, путем экспериментов, подбирается.
Теперь рассказывайте вы.

 

Точно в яблочко!
Один из возможных физических способов задания линейного масштаба. Другой – скорость света.

 

Выкладываю результаты. В таблицах приведены пересчет с 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-шниками? А то у меня не получается скомпилировать фортрановские файлы.

 

GGMCalc - исправленная версия.

 

У меня hsynth правильно работает только в версии на sparc/solaris (полное совпадение тестов с NGA).
Сейчас пытаюсь разобраться с порядком следования коэффициентов в бинарном файле GAO98.jff для pinnacle.
Нет ли у кого завалявшегося файла EGM96.jff для сравнения ?