Добро пожаловать!

Войдите или зарегистрируйтесь сейчас!

Войти

ГАО2012

Тема в разделе "Геодезия как наука", создана пользователем Родичкин, 18 июл 2017.

  1. cr2

    cr2
    Форумчанин

    Регистрация:
    23 сен 2014
    Сообщения:
    318
    Симпатии:
    180
    Какой азимут был исходным не так принципиально, как тот факт, что геодезический азимут Саблино-Бугры в СК-42 равен 317°02'54.29", а в СК-32 был равен 317°02'50.62" (это все открытые "публикабельные" данные).
    В статье "Саблинская базисная сеть" был вычислен геодезический азимут Бугры-Саблино в WGS84
    равный 136°54'48,0", его надо только развернуть.
    --- Сообщения объединены, 14 янв 2019, Оригинальное время сообщения: 14 янв 2019 ---
    pul-map-en.jpg
     
    #401
    Yuri V. и andrew_klikunov нравится это.
  2. andrew_klikunov

    Форумчанин

    Регистрация:
    30 сен 2015
    Сообщения:
    189
    Симпатии:
    41
    Короче говоря, надо пробовать пересчитывать высоты по тем параметрам, что получились и, смотреть разницу. Благо данных для контроля теперь стало больше. Может какую-то зависимость можно будет увидеть. Потом попробовать с поэкспериментировать с Сигналом А.

    А, как Вы думаете, сильно ли разворот в 4" может повлиять на результат? Просто интересно мнение, как специалиста.
     
    #402
  3. cr2

    cr2
    Форумчанин

    Регистрация:
    23 сен 2014
    Сообщения:
    318
    Симпатии:
    180
    Историческое значение высоты на довоенной топокарте равно 74.9 метра и это дожна была быть высота "старого сигнал А" (новый восстановлен в 1981 году). Я сейчас не помню его высоту в саженях 19 века :)
    Астростолб это новодел 1973 года.
    --- Сообщения объединены, 14 янв 2019, Оригинальное время сообщения: 14 янв 2019 ---
    У меня есть такое мнение что СК-32 должна прекрасно соответствовать всем современным требованиям
    для кадастровых работ в Ленинградской, Новгородской и Псковской областях.
    Так как граница с Норвегией (и северный кусок с Финляндией) была
    (до самого последнего времени) демаркирована в СК-32,
    то можно оценить ее точность в Карелии и Мурманской области.
    Для Чехии надо считать.
     
    #403
    andrew_klikunov нравится это.
  4. Yuri V.

    Форумчанин

    Регистрация:
    31 мар 2009
    Сообщения:
    1.467
    Симпатии:
    1.616
    Адрес:
    Ивантеевка, РФ
    Уважаемые cr2, andrew_klikunov !
    Коль скоро здесь обсуждается некий "инструмент" для отождествления геодезической высоты и нормальной, неужели нельзя обойтись без азимутов исходных дат. А учебники говорят вот о чем: при установлении ИГД СК-42 эллипсоид был выбран из нескольких условий, одно из которых равенство нормального и возмущающего потенциалов в точке счета высот.
    Разве имеет богатый смысл исходная отметка? По-моему смысл находится в точной передаче этого значения дальше. пф.. и здесь кувырк, значения ошибок растут, не только на восток (но с прогрессией), сколько на северо-восток.
     
    #404
  5. andrew_klikunov

    Форумчанин

    Регистрация:
    30 сен 2015
    Сообщения:
    189
    Симпатии:
    41
    Уважаемые stout и Yuri V., не могли бы Вы разъяснить, для чего в WGS 84 используются каждая из двух разных GM? Для GRS-80 всего-то одна... Просто есть любопытные на мой взгляд результаты по расчету. Я их обязательно приложу, но в начале хотелось бы узнать по GM.
     
    #405
  6. stout

    Форумчанин Форумчанин

    Регистрация:
    5 янв 2008
    Сообщения:
    3.563
    Симпатии:
    11.337
    Адрес:
    Златоглавая и Белокаменная
    Объяснение содержится в
    основополагающем документе:
    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.
    "
     
    #406
    Yuri V. и andrew_klikunov нравится это.
  7. andrew_klikunov

    Форумчанин

    Регистрация:
    30 сен 2015
    Сообщения:
    189
    Симпатии:
    41
    Как я понял, значение 3986004.418 x 10 8 m 3 /s 2 используют для высокоточных определений орбит. Это правильно? Но, и в программе h_synth EGM2008 оно тоже используется. А вот зачем нужно значение 3986005.0 x 10 8 m 3 /s 2 мне так и непонятно.

    Оффтоп
    Пытался перевести через google, но там как-то вообще невнятно получилось
     
    #407
  8. stout

    Форумчанин Форумчанин

    Регистрация:
    5 янв 2008
    Сообщения:
    3.563
    Симпатии:
    11.337
    Адрес:
    Златоглавая и Белокаменная

    Грубо говоря, уже были тысячи приёмников, которые использовали для расчёта координат верха старое значение геогравитационной постоянной. Введение нового значения превратило бы их в кирпич.
    Странно, у меня 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 должна была передаваться в размерности км32. Один товарищ, из вполне здравых рассуждений, что работаем в системе СИ передал значение в м32. Как следствие, спутник летал под поверхностью Земли на глубине в 1000 км. Невязки O-C (observantum minus calculantum) были "прекрасны"
     
    #408
  9. Yuri V.

    Форумчанин

    Регистрация:
    31 мар 2009
    Сообщения:
    1.467
    Симпатии:
    1.616
    Адрес:
    Ивантеевка, РФ
    Оффтоп
    почему не понятно, гугл нормально перевёл:
    «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, размерность км32 - ускорение объёма без массы ::crazy::. Этим ученые ушли от проблемы точного определения массы Земли, которая необходима для составления уравнений движения ИСЗ, по закону Всемирного тяготения. Короче (если дилетанту позволительно что-то укорачивать), это такой коэффициент/константа, регулирующая метрический (в декартовых видимо) масштаб описания ГПЗ. Итеративно, путем экспериментов, подбирается.
    Теперь рассказывайте вы.
     
    #409
    stout нравится это.
  10. stout

    Форумчанин Форумчанин

    Регистрация:
    5 янв 2008
    Сообщения:
    3.563
    Симпатии:
    11.337
    Адрес:
    Златоглавая и Белокаменная
    Точно в яблочко!
    Один из возможных физических способов задания линейного масштаба. Другой – скорость света.
     
    #410
  11. andrew_klikunov

    Форумчанин

    Регистрация:
    30 сен 2015
    Сообщения:
    189
    Симпатии:
    41
    Выкладываю результаты. В таблицах приведены пересчет с GRS-80 и WGS 84 (маркировал цветами). points - сводная табличка, вторая табличка - с полным расчетом (точки приведены в том же порядке, что и в сводной).
     

    Вложения:

    #411
  12. stout

    Форумчанин Форумчанин

    Регистрация:
    5 янв 2008
    Сообщения:
    3.563
    Симпатии:
    11.337
    Адрес:
    Златоглавая и Белокаменная
    Вот две интересных ссылочки
    https://www.svf.stuba.sk/en/departm...ence-and-research/downloads.html?page_id=4996
    http://imi.cas.sc.edu/gravimetric-quantities/
    Статьи с описанием Soft'a есть на Генезисе.
    Программки по второй ссылке зело охочие до памяти да и подготовительный этап неспешный.
    Но считает именно аномалию высоты, а не высоту геоида. И быстро считает
    Values per second.png
    ============================================================================
    Подготовительный этап на моём компе
    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]
     
    #412
  13. andrew_klikunov

    Форумчанин

    Регистрация:
    30 сен 2015
    Сообщения:
    189
    Симпатии:
    41
    Не совсем понятен Ваш замысел. Вы хотите, чтобы я посчитал аномалию высоты с помощью этого ПО и потом сделал пересчет на элл. Красовского? P.S.: забыл написать, что в тех таблицах использованы значения высоты геоида из источников, приведенных выше (в статье из журнала), т.е. я не использовал EGM2008.
     
    #413
  14. stout

    Форумчанин Форумчанин

    Регистрация:
    5 янв 2008
    Сообщения:
    3.563
    Симпатии:
    11.337
    Адрес:
    Златоглавая и Белокаменная
    Если коротко, то – да.
    Я потерял нить обсуждения, надо перечитать и восстановить.
     
    #414
  15. andrew_klikunov

    Форумчанин

    Регистрация:
    30 сен 2015
    Сообщения:
    189
    Симпатии:
    41
    Началось всё с того, что нашлись контрольные данные по высотам квазигеоида над эллипсоидом Красовского. Они приведены тут. Это территория бывшей Чехословакии. Там, так же приведены и источники. Формулы приведены тут. Результаты ( в т.ч. и дельты по осям, которые использовались в формуле расчета разницы) приведены здесь. На эти данные я и опирался. Вроде ничего не упустил.
     
    #415
  16. X-Y-H

    X-Y-H Администратор
    Команда форума Форумчанин Форумчанин

    Регистрация:
    18 май 2007
    Сообщения:
    17.530
    Симпатии:
    5.113
    Адрес:
    Россия
    #416
  17. andrew_klikunov

    Форумчанин

    Регистрация:
    30 сен 2015
    Сообщения:
    189
    Симпатии:
    41
    Вы имеете ввиду пересчет по ГОСТу 2008? Я тут исходил вот из чего: для расчета силы тяжести над эллипсоидом +- 50м по широте - не критично. Кстати ГОСТ выдал +- 2м по высоте.
     
    #417
  18. andrew_klikunov

    Форумчанин

    Регистрация:
    30 сен 2015
    Сообщения:
    189
    Симпатии:
    41
    stout, Вы не могли бы поделиться exe-шниками? А то у меня не получается скомпилировать фортрановские файлы.
     
    #418
    Последнее редактирование: 29 янв 2019
  19. stout

    Форумчанин Форумчанин

    Регистрация:
    5 янв 2008
    Сообщения:
    3.563
    Симпатии:
    11.337
    Адрес:
    Златоглавая и Белокаменная
    GGMCalc - исправленная версия.
     

    Вложения:

    #419
    andrew_klikunov нравится это.
  20. cr2

    cr2
    Форумчанин

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

Поделиться этой страницей

  1. Этот сайт использует файлы cookie. Продолжая пользоваться данным сайтом, Вы соглашаетесь на использование нами Ваших файлов cookie.
    Скрыть объявление