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

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

Войти

  1. С 13.10.2020 снова работает авторизация/регистрация через социальные сети: VK (Вконтакте) и Facebook.
    Скрыть объявление

Эллипсоид и параметры ГСК-2011 в ГИС программах

Тема в разделе "ГИС", создана пользователем X-Y-H, 5 июн 2019.

  1. Slinger

    Форумчанин

    Регистрация:
    13 окт 2010
    Сообщения:
    1.085
    Симпатии:
    453
    Гаврилов Максим (Oiz), Допустим. То есть раньше координаты СК42 были схожи с СК95 (и наоборот). А теперь они еще и схожи с ГСК2011. И какой в этом толк, создавать копию СК42. Тем более сейчас выполнил пару тестовых вычислений и оказалось, что разница между датумом СК42 и датумом ГСК2011 (в одной и той же проекции) ничтожна


    stout, А как думаете Вы ?
     
    #41
  2. zvezdochiot

    Форумчанин

    Регистрация:
    27 июн 2014
    Сообщения:
    1.962
    Симпатии:
    714
    Адрес:
    г. Москва
    "Ты мне не юли, ты пальцем покажи!"

    Выложите их (датумы) уже наконец!
     
    #42
  3. Гаврилов Максим (Oiz)

    Форумчанин

    Регистрация:
    19 авг 2009
    Сообщения:
    590
    Симпатии:
    202
    Адрес:
    Новосибирск
    Эллипсойд то другой. У меня плановые координаты примерно на 40-50 м по x - y изменились.
    --- Сообщения объединены, 26 фев 2021, Оригинальное время сообщения: 26 фев 2021 ---
    За отсчетную поверхность в ГСК-2011 принят общеземной эллипсоид с большой полуосью равной, 6378136,5 м, и сжатием равным 2564151
     
    #43
  4. StudentX

    Форумчанин

    Регистрация:
    1 июн 2019
    Сообщения:
    1.470
    Симпатии:
    573
    Адрес:
    Пермь
    Оффтоп
    ?!::blink.gif:: Это что за параметр такой? Полярное сжатие земных эллипсоидов порядка 1/300. А 2564151 - это эллипсоид, сжатый до состояния плоскости меридиана.
     
    #44
    zvezdochiot нравится это.
  5. zvezdochiot

    Форумчанин

    Регистрация:
    27 июн 2014
    Сообщения:
    1.962
    Симпатии:
    714
    Адрес:
    г. Москва
    Оффтоп
    Это "хвост" после запятой. Но написано вырвиглазно!
     
    #45
    StudentX и Qvinto нравится это.
  6. Slinger

    Форумчанин

    Регистрация:
    13 окт 2010
    Сообщения:
    1.085
    Симпатии:
    453
    zvezdochiot, Для ArcGIS параметры трансформаций есть тут https://blogs.esri-cis.ru/2018/02/14/geo-transformations/

    Для MapInfo свежие параметры выставлены здесь http://mapinfo.ru/articles/gauss-kruger-gsk-2011 Эти же параметры подтверждают, что проекция осталось той же. Поменялся только эллипсоид и датум

    Для QGIS руки пока не дошли собрать

    Для Global Mapper почти собрал


    Гаврилов Максим (Oiz), Да, я теперь тоже убедился в смещении. Потестил в MapInfo. Но разница сильно заметна на восточных регионах, то есть начиная хотя бы с 10-й зоны. (в моем примере были объекты из 14-й зоны и разница составила примерно метров 30) На территории Москвы (7-я зона) смещение и разница между СК-42 и ГСК-2011 от силы сантиметров 10 у меня получилась
     
    #46
    GeoCod нравится это.
  7. X-Y-H

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

    Регистрация:
    18 май 2007
    Сообщения:
    18.862
    Симпатии:
    5.647
    Адрес:
    Россия
    по цифрам да, по факту на земле до 100 м
     
    #47
  8. zvezdochiot

    Форумчанин

    Регистрация:
    27 июн 2014
    Сообщения:
    1.962
    Симпатии:
    714
    Адрес:
    г. Москва
    Оффтоп
    Ну а иначе зачем тогда этот "переход" нужен был бы? Работали бы "по-старинке".
     
    #48
    Qvinto нравится это.
  9. uriygr

    Форумчанин

    Регистрация:
    14 апр 2008
    Сообщения:
    1.078
    Симпатии:
    1.390
    Адрес:
    Россия, Самара
    Оффтоп
    Это теодолит, светодальномер, нивелир?
     
    #49
  10. zvezdochiot

    Форумчанин

    Регистрация:
    27 июн 2014
    Сообщения:
    1.962
    Симпатии:
    714
    Адрес:
    г. Москва
    Оффтоп
    Светодальномер здесь лишний! )))
     
    #50
  11. geo-prog

    Форумчанин

    Регистрация:
    17 сен 2014
    Сообщения:
    170
    Симпатии:
    64
    Должно быть (знаки разворота перепутали):
    +proj=tmerc +lat_0=0 +lon_0=39 +k=1 +x_0=7500000 +y_0=0 +a=6378136.5 +rf=298.2564151 +towgs84=0.013,-0.092,-0.03,-0.001738,0.003559,-0.004263,0.0074 +units=m +no_defs
     
    #51
    MARK SHEYDER нравится это.
  12. Гаврилов Максим (Oiz)

    Форумчанин

    Регистрация:
    19 авг 2009
    Сообщения:
    590
    Симпатии:
    202
    Адрес:
    Новосибирск
    Да вы все верно написали. У себя исправил, брал данные как из ТБС, но как обсудили выше в этой теме в qgis используется другое преобразование. Поэтому знаки разворота осей противоположные.
     
    #52
  13. zvezdochiot

    Форумчанин

    Регистрация:
    27 июн 2014
    Сообщения:
    1.962
    Симпатии:
    714
    Адрес:
    г. Москва
    Э-э-эм. Что то я запутался, "знаки разворота осей" или "знаки разворота ВОКРУГ осей"?
     
    #53
  14. X-Y-H

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

    Регистрация:
    18 май 2007
    Сообщения:
    18.862
    Симпатии:
    5.647
    Адрес:
    Россия
    Ну для QGIS и некоторых других программ актуальнее WKT
    WKT
    Код:
    BOUNDCRS[
        SOURCECRS[
            PROJCRS["unknown",
                BASEGEOGCRS["unknown",
                    DATUM["unknown",
                        ELLIPSOID["GSK-2011",6378136.5,298.2564151,
                            LENGTHUNIT["metre",1,
                                ID["EPSG",9001]]]],
                    PRIMEM["Greenwich",0,
                        ANGLEUNIT["degree",0.0174532925199433],
                        ID["EPSG",8901]]],
                CONVERSION["unknown",
                    METHOD["Transverse Mercator",
                        ID["EPSG",9807]],
                    PARAMETER["Latitude of natural origin",0,
                        ANGLEUNIT["degree",0.0174532925199433],
                        ID["EPSG",8801]],
                    PARAMETER["Longitude of natural origin",51,
                        ANGLEUNIT["degree",0.0174532925199433],
                        ID["EPSG",8802]],
                    PARAMETER["Scale factor at natural origin",1,
                        SCALEUNIT["unity",1],
                        ID["EPSG",8805]],
                    PARAMETER["False easting",9500000,
                        LENGTHUNIT["metre",1],
                        ID["EPSG",8806]],
                    PARAMETER["False northing",0,
                        LENGTHUNIT["metre",1],
                        ID["EPSG",8807]]],
                CS[Cartesian,2],
                    AXIS["(E)",east,
                        ORDER[1],
                        LENGTHUNIT["metre",1,
                            ID["EPSG",9001]]],
                    AXIS["(N)",north,
                        ORDER[2],
                        LENGTHUNIT["metre",1,
                            ID["EPSG",9001]]]]],
        TARGETCRS[
            GEOGCRS["WGS 84",
                DATUM["World Geodetic System 1984",
                    ELLIPSOID["WGS 84",6378137,298.257223563,
                        LENGTHUNIT["metre",1]]],
                PRIMEM["Greenwich",0,
                    ANGLEUNIT["degree",0.0174532925199433]],
                CS[ellipsoidal,2],
                    AXIS["latitude",north,
                        ORDER[1],
                        ANGLEUNIT["degree",0.0174532925199433]],
                    AXIS["longitude",east,
                        ORDER[2],
                        ANGLEUNIT["degree",0.0174532925199433]],
                ID["EPSG",4326]]],
        ABRIDGEDTRANSFORMATION["Transformation from unknown to WGS84",
            METHOD["Position Vector transformation (geog2D domain)",
                ID["EPSG",9606]],
            PARAMETER["X-axis translation",0.013,
                ID["EPSG",8605]],
            PARAMETER["Y-axis translation",-0.092,
                ID["EPSG",8606]],
            PARAMETER["Z-axis translation",-0.03,
                ID["EPSG",8607]],
            PARAMETER["X-axis rotation",-0.001738,
                ID["EPSG",8608]],
            PARAMETER["Y-axis rotation",0.003559,
                ID["EPSG",8609]],
            PARAMETER["Z-axis rotation",-0.004263,
                ID["EPSG",8610]],
            PARAMETER["Scale difference",1.0000000074,
                ID["EPSG",8611]]]]
     
    #54
  15. zvezdochiot

    Форумчанин

    Регистрация:
    27 июн 2014
    Сообщения:
    1.962
    Симпатии:
    714
    Адрес:
    г. Москва
    Так @cr2 вроде шаблоны уже выкладывал? Или @ErnieBoyd , я уже не помню.
     
    #55
  16. Гаврилов Максим (Oiz)

    Форумчанин

    Регистрация:
    19 авг 2009
    Сообщения:
    590
    Симпатии:
    202
    Адрес:
    Новосибирск
    Да что же такое, конечно термины надо использовать верно. Знаки разворота вокруг координатных осей.
     
    #56
  17. ErnieBoyd

    Форумчанин

    Регистрация:
    10 июн 2014
    Сообщения:
    215
    Симпатии:
    107
    В базе EPSG есть датум "Geodezicheskaya Sistema Koordinat 2011" и всё, что с ним связано, так что не парьтесь о разворотах.
    Пользовательская проекция "GSK-2011 / Gauss-Kruger zone 7":
    Код:
    PROJCRS["GSK-2011 / Gauss-Kruger zone 7",
        BASEGEOGCRS["GSK-2011",
            DATUM["Geodezicheskaya Sistema Koordinat 2011",
                ELLIPSOID["GSK-2011",6378136.5,298.2564151,
                    LENGTHUNIT["metre",1]]],
            PRIMEM["Greenwich",0,
                ANGLEUNIT["degree",0.0174532925199433]],
            ID["EPSG",7683]],
        CONVERSION["6-degree Gauss-Kruger zone 7",
            METHOD["Transverse Mercator",
                ID["EPSG",9807]],
            PARAMETER["Latitude of natural origin",0,
                ANGLEUNIT["degree",0.0174532925199433],
                ID["EPSG",8801]],
            PARAMETER["Longitude of natural origin",39,
                ANGLEUNIT["degree",0.0174532925199433],
                ID["EPSG",8802]],
            PARAMETER["Scale factor at natural origin",1,
                SCALEUNIT["unity",1],
                ID["EPSG",8805]],
            PARAMETER["False easting",7500000,
                LENGTHUNIT["metre",1],
                ID["EPSG",8806]],
            PARAMETER["False northing",0,
                LENGTHUNIT["metre",1],
                ID["EPSG",8807]]],
        CS[Cartesian,2],
            AXIS["northing (X)",north,
                ORDER[1],
                LENGTHUNIT["metre",1]],
            AXIS["easting (Y)",east,
                ORDER[2],
                LENGTHUNIT["metre",1]],
        USAGE[
            SCOPE["Topographic mapping (medium scale)."],
            AREA["Russian Federation - onshore 36°E to 42°E"],
            BBOX[41.43,36,69.23,42]]]
     
    #57
    Последнее редактирование: 2 мар 2021
  18. geo-prog

    Форумчанин

    Регистрация:
    17 сен 2014
    Сообщения:
    170
    Симпатии:
    64
    Вот вам контрольный пример посчитанный в QGIS и photomod geokalсulator (ГОСТ 32453-2017) с моими параметрами Контрольный пример.jpg
     
    #58
  19. ErnieBoyd

    Форумчанин

    Регистрация:
    10 июн 2014
    Сообщения:
    215
    Симпатии:
    107
    Беру свои слова обратно. Система координат есть, но к ней, похоже, нет трансформации. Поэтому надо париться создавать системы координат с пользовательской трансформацией.

    Кстати, в WKT знаки компонент вращения могут задаваться по методу "Coordinate Frame rotation", как в ГОСТе. Пример - пользовательская проекция "GSK-2011 / Gauss-Kruger zone 7":
    Код:
    BOUNDCRS[
        SOURCECRS[
            PROJCRS["unknown",
                BASEGEOGCRS["unknown",
                    DATUM["Unknown based on GSK-2011 ellipsoid",
                        ELLIPSOID["GSK-2011",6378136.5,298.2564151,
                            LENGTHUNIT["metre",1,
                                ID["EPSG",9001]]]],
                    PRIMEM["Greenwich",0,
                        ANGLEUNIT["degree",0.0174532925199433],
                        ID["EPSG",8901]]],
                CONVERSION["unknown",
                    METHOD["Transverse Mercator",
                        ID["EPSG",9807]],
                    PARAMETER["Latitude of natural origin",0,
                        ANGLEUNIT["degree",0.0174532925199433],
                        ID["EPSG",8801]],
                    PARAMETER["Longitude of natural origin",39,
                        ANGLEUNIT["degree",0.0174532925199433],
                        ID["EPSG",8802]],
                    PARAMETER["Scale factor at natural origin",1,
                        SCALEUNIT["unity",1],
                        ID["EPSG",8805]],
                    PARAMETER["False easting",7500000,
                        LENGTHUNIT["metre",1],
                        ID["EPSG",8806]],
                    PARAMETER["False northing",0,
                        LENGTHUNIT["metre",1],
                        ID["EPSG",8807]]],
                CS[Cartesian,2],
                    AXIS["(E)",east,
                        ORDER[1],
                        LENGTHUNIT["metre",1,
                            ID["EPSG",9001]]],
                    AXIS["(N)",north,
                        ORDER[2],
                        LENGTHUNIT["metre",1,
                            ID["EPSG",9001]]]]],
        TARGETCRS[
            GEOGCRS["WGS 84",
                DATUM["World Geodetic System 1984",
                    ELLIPSOID["WGS 84",6378137,298.257223563,
                        LENGTHUNIT["metre",1]]],
                PRIMEM["Greenwich",0,
                    ANGLEUNIT["degree",0.0174532925199433]],
                CS[ellipsoidal,2],
                    AXIS["geodetic latitude (Lat)",north,
                        ORDER[1],
                        ANGLEUNIT["degree",0.0174532925199433]],
                    AXIS["geodetic longitude (Lon)",east,
                        ORDER[2],
                        ANGLEUNIT["degree",0.0174532925199433]],
                ID["EPSG",4326]]],
        ABRIDGEDTRANSFORMATION["Transformation from unknown to WGS84",
            METHOD["Coordinate Frame rotation (geog2D domain)",
                ID["EPSG",9607]],
            PARAMETER["X-axis translation",0.013,
                ID["EPSG",8605]],
            PARAMETER["Y-axis translation",-0.092,
                ID["EPSG",8606]],
            PARAMETER["Z-axis translation",-0.03,
                ID["EPSG",8607]],
            PARAMETER["X-axis rotation",0.001738,
                ID["EPSG",8608]],
            PARAMETER["Y-axis rotation",-0.003559,
                ID["EPSG",8609]],
            PARAMETER["Z-axis rotation",0.004263,
                ID["EPSG",8610]],
            PARAMETER["Scale difference",1.0000000074,
                ID["EPSG",8611]]]]
    --- Сообщения объединены, 7 мар 2021, Оригинальное время сообщения: 7 мар 2021 ---
    "Coordinate Frame rotation" - это вращение координатных осей, а "Position Vector transformation" - вращение пространства вогруг осей.

    Способ Гельмерта - это математически точное конформное преобразование, использующее матрицу вращения. Формула Бурша-Вольфа - упрощённое представление способа Гельмерта, при котором компоненты матрицы вращения раскладываются в ряды Тейлора по степеням углов вращения с отбрасыванием компонент второго и более высоких порядков.
     
    #59
  20. X-Y-H

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

    Регистрация:
    18 май 2007
    Сообщения:
    18.862
    Симпатии:
    5.647
    Адрес:
    Россия
    Рисунка про это нет? Своим надо объяснить вечно путают.
     
    #60
    zvezdochiot нравится это.

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

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