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

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


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

 

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

Выложите их (датумы) уже наконец!

 

Эллипсойд то другой. У меня плановые координаты примерно на 40-50 м по x - y изменились.

--- Сообщения объединены, 26 фев 2021, Оригинальное время сообщения: 26 фев 2021 ---

За отсчетную поверхность в ГСК-2011 принят общеземной эллипсоид с большой полуосью равной, 6378136,5 м, и сжатием равным 2564151

 

Оффтоп

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

 

Оффтоп

Это "хвост" после запятой. Но написано вырвиглазно!

 

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 у меня получилась

 

Оффтоп

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

 

Оффтоп

Это теодолит, светодальномер, нивелир?

 

Оффтоп

Светодальномер здесь лишний! )))

 

Должно быть (знаки разворота перепутали):
+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

 

Да вы все верно написали. У себя исправил, брал данные как из ТБС, но как обсудили выше в этой теме в qgis используется другое преобразование. Поэтому знаки разворота осей противоположные.

 

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

 

Так @cr2 вроде шаблоны уже выкладывал? Или @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]]]

 

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

 

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

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

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