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

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

  1. Slinger

    Slinger Форумчанин

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


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

    zvezdochiot Форумчанин

    "Ты мне не юли, ты пальцем покажи!"

    Выложите их (датумы) уже наконец!
     
  3. Эллипсойд то другой. У меня плановые координаты примерно на 40-50 м по x - y изменились.
    --- Сообщения объединены, 26 фев 2021, Оригинальное время сообщения: 26 фев 2021 ---
    За отсчетную поверхность в ГСК-2011 принят общеземной эллипсоид с большой полуосью равной, 6378136,5 м, и сжатием равным 2564151
     
  4. Deleted member 122005

    Deleted member 122005 Только чтение

    Оффтоп
    ?!::blink.gif:: Это что за параметр такой? Полярное сжатие земных эллипсоидов порядка 1/300. А 2564151 - это эллипсоид, сжатый до состояния плоскости меридиана.
     
    zvezdochiot нравится это.
  5. zvezdochiot

    zvezdochiot Форумчанин

    Оффтоп
    Это "хвост" после запятой. Но написано вырвиглазно!
     
    Deleted member 122005 и Qvinto нравится это.
  6. Slinger

    Slinger Форумчанин

    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 у меня получилась
     
    GeoCod нравится это.
  7. X-Y-H

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

    по цифрам да, по факту на земле до 100 м
     
  8. zvezdochiot

    zvezdochiot Форумчанин

    Оффтоп
    Ну а иначе зачем тогда этот "переход" нужен был бы? Работали бы "по-старинке".
     
    Qvinto нравится это.
  9. uriygr

    uriygr Форумчанин

    Оффтоп
    Это теодолит, светодальномер, нивелир?
     
  10. zvezdochiot

    zvezdochiot Форумчанин

    Оффтоп
    Светодальномер здесь лишний! )))
     
  11. Deleted member 51242

    Deleted member 51242 Только чтение

    Должно быть (знаки разворота перепутали):
    +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
     
    zvezdochiot и MARK SHEYDER нравится это.
  12. Да вы все верно написали. У себя исправил, брал данные как из ТБС, но как обсудили выше в этой теме в qgis используется другое преобразование. Поэтому знаки разворота осей противоположные.
     
  13. zvezdochiot

    zvezdochiot Форумчанин

    Э-э-эм. Что то я запутался, "знаки разворота осей" или "знаки разворота ВОКРУГ осей"?
     
  14. X-Y-H

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

    Ну для 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]]]]
     
  15. zvezdochiot

    zvezdochiot Форумчанин

    Так @cr2 вроде шаблоны уже выкладывал? Или @ErnieBoyd , я уже не помню.
     
  16. Да что же такое, конечно термины надо использовать верно. Знаки разворота вокруг координатных осей.
     
  17. ErnieBoyd

    ErnieBoyd Форумчанин

    В базе 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]]]
     
    Последнее редактирование: 2 мар 2021
  18. Deleted member 51242

    Deleted member 51242 Только чтение

    Вот вам контрольный пример посчитанный в QGIS и photomod geokalсulator (ГОСТ 32453-2017) с моими параметрами Контрольный пример.jpg
     
  19. ErnieBoyd

    ErnieBoyd Форумчанин

    Беру свои слова обратно. Система координат есть, но к ней, похоже, нет трансформации. Поэтому надо париться создавать системы координат с пользовательской трансформацией.

    Кстати, в 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" - вращение пространства вогруг осей.

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

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

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