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

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

Войти

Описание систем координат PROJ и WKT

Тема в разделе "Исходные данные", создана пользователем cr2, 10 ноя 2020.

  1. cr2

    cr2
    Форумчанин

    Регистрация:
    23 сен 2014
    Сообщения:
    781
    Симпатии:
    511
    Самой распространенной ошибкой перевода в "WGS84" является неучет смены знака
    углов поворота в параметрах Mapinfo и TOWGS84[] (он такой же как у +towgs84 из proj4).
    Таких ошибок полно и в базе EPSG, поэтому для QGIS надо всегда задавать свои
    пользовательские системы координат: ГОСТ для совместимости с Росреестром
    (+towgs84=23.57,-140.95,-79.8,0,0.35,0.79,-0.22 он же TOWGS84[23.57,-140.95,-79.8,0,0.35,0.79,-0.22])
    и свои собственные для максимальной точности.
     
    #1
  2. ErnieBoyd

    Форумчанин

    Регистрация:
    10 июн 2014
    Сообщения:
    271
    Симпатии:
    159
    • Вы сильно ошибаетесь насчёт ошибок (тем более их обилия) в базе EPSG. MapInfo и ГОСТ используют трансформацию EPSG:9607 "Coordinate frame rotation". Параметр TOWGS84 использовал трансформацию EPSG:9606 "Position vector transformastion". Они различаются знаками параметров вращения.
    • Параметр TOWGS84 больше не работает. Теперь в QGIS (как и в ArcGIS, Blue Marble Geographic Calculator и др.) системы координат связаны с наборами трансформаций, из которых надо выбирать подходящую. Разумеется, ГОСТ-овские среди них имеются. Чтобы не выбирать трансформацию для системы координат каждый раз, назначаете нужную по умолчанию.
    • И о какой такой совместимости Росреестра с ГОСТом Вы говорите?
     
    #2
    Последнее редактирование: 11 ноя 2020
  3. X-Y-H

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

    Регистрация:
    18 май 2007
    Сообщения:
    21.792
    Симпатии:
    7.071
    Адрес:
    Россия
    в qgis вроде не работает. в GM работает
     
    #3
  4. cr2

    cr2
    Форумчанин

    Регистрация:
    23 сен 2014
    Сообщения:
    781
    Симпатии:
    511
    Давайте рассмотрим для примера EPSG:3352
    в ее инкарнации
    GOST R 51794-2008, approved by GosStandard of Russia Decree #609 of December 18, 2008.
    http://epsg.io/3352-5044

    Я честно говоря и не заметил этого изменения во внутреннем-обменном формате представления параметров
    проекций в QGIS3.10. Как пользовался proj4 строками, так и продолжаю ими пользоваться

    Росреестр (точнее его субконтракторы) использует GOST R 51794-2008 в сервисах ПКК
    для перевода данных ЕГРН в сферический меркатор.
     
    #4
  5. ErnieBoyd

    Форумчанин

    Регистрация:
    10 июн 2014
    Сообщения:
    271
    Симпатии:
    159
    Вы ссылаетесь на ресурс epsg.io. Между тем его разработчик MapTiler Team не является официальным представителем EPSG. Не стоит доверять их экспортным примерам. Лучше посмотрите на представление в QGIS:
    Код:
    Pulkovo 1942 / CS63 zone C2
    WKT
    PROJCRS["Pulkovo 1942 / CS63 zone C2",
    …
        ID["EPSG",3352]] 
    Proj4
    +proj=tmerc +lat_0=0.1 +lon_0=27.95 +k=1 +x_0=2250000 +y_0=0 +ellps=krass +towgs84=23.57,-140.95,-79.8,0,0.35,0.79,-0.22 +units=m +no_defs
    …
    Как видите, в БД знаки правильные.
    --- Сообщения объединены, 12 ноя 2020, Оригинальное время сообщения: 12 ноя 2020 ---
    QGIS 3.10.x
    x > 6 ?

    Не существует данных, которые хранят описание системы координат в формате proj.
    Откройте слой, сохранённый в МСК, и посмотрите в свойствах систему координат. Она будет неопределённой.

    В стандартах OGC и ISO формат данных GeoPackage, формат представления систем координат WKT версия 2.
    --- Сообщения объединены, 12 ноя 2020 ---
    Вот не сказал бы, глядя на ПКК и сравнивая объекты на ней с координатами из выписок.
     
    #5
  6. cr2

    cr2
    Форумчанин

    Регистрация:
    23 сен 2014
    Сообщения:
    781
    Симпатии:
    511
    Согласен. На официальном сайте EPSG этих параметров вообще нет,
    а для скачивания базы надо регистрироваться.

    У меня QGIS 3.10.11-A и в нем отсутствуют 7 параметров для CS63/C2
    Код:
     
     Proj4
    +proj=tmerc +lat_0=0.1 +lon_0=27.95 +k=1 +x_0=2250000 +y_0=0 +ellps=krass +units=m +no_defs
    
    При этом proj.db содержит правильную информацию для EPSG:5044
    Код:
    INSERT INTO helmert_transformation_table 
    VALUES('EPSG','5044','Pulkovo 1942 to WGS 84 (20)','
    Derived through concatenation of Pulkovo 1942 to PZ-90.02 to WGS 84. Replaces Pulkovo 1942 to WGS 84 (17) (code 1267).',
    'EPSG','9607','EPSG','4284','EPSG','4326',3.0,23.570000000000000284,-140.94999999999998863,-79.799999999999997159,
    'EPSG','9001',0.0,-0.34999999999999997779,-0.79000000000000003552,
    'EPSG','9104',-0.22000000000000000111,
    'EPSG','9202',NULL,NULL,NULL,NULL,NULL,NULL,NULL,NULL,NULL,NULL,NULL,NULL,NULL,NULL,NULL,NULL,NULL,NULL,NULL,NULL,NULL,
    'GOST-Rus 2008',0);
    
    Надо отметить что исторически proj4 использовал EPSG:9606(в/для +towgs84=), поэтому нынешний переход в базе на 'EPSG','9607'
    требует очень внимательного подхода.
    В Санкт-Петербурге и Ленинградской области ПКК использует ГОСТ-2008. В других регионах возможно это не так,
    приношу извинения за чрезмерное обобщение.
     
    #6
  7. ErnieBoyd

    Форумчанин

    Регистрация:
    10 июн 2014
    Сообщения:
    271
    Симпатии:
    159
    Сами параметры epsg.io выдаёт правильно. Посмотрите ту же EPSG:5044. Косяки только в экспортных примерах.

    Об этом я и твержу! Pulkovo 1942 / CS63 zone C2 на датуме "Pulkovo 1942", ассоциированном с разными трансформациями. Программа даёт возможность выбрать россиянам одну трансформацию, эстонцам другую, латвийцам третью.


    epsg5044.png epsg3352.png
     
    #7
  8. Сантьяго

    Форумчанин

    Регистрация:
    9 янв 2008
    Сообщения:
    898
    Симпатии:
    354
    Адрес:
    Краснодар
    Здравствуйте
    подскажите я правильно сделал prj для СК 26 от СК1995г
    у меня сомнения еcть в кодах "EPSG"

    PROJCS["msk26_2",GEOGCS["Pulkovo 1995",DATUM["Pulkovo 1995",SPHEROID["Krassowsky 1940",6378245,298.3,AUTHORITY["EPSG","7024"]],TOWGS84[24.47,-130.89,-81.56,0,0,-0.13,0.22],AUTHORITY["EPSG","4200"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,AUTHORITY["EPSG","9102"]],AUTHORITY["EPSG","4200"]],PROJECTION["Transverse_Mercator",AUTHORITY["EPSG","9807"]],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",44.75],PARAMETER["scale_factor",1],PARAMETER["false_easting",2300000],PARAMETER["false_northing",-4512890.197998],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AUTHORITY["EPSG","msk26_2"]]
     
    #8
  9. X-Y-H

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

    Регистрация:
    18 май 2007
    Сообщения:
    21.792
    Симпатии:
    7.071
    Адрес:
    Россия
    Выкинь - пиши везде CUSTOM
     
    #9
  10. ErnieBoyd

    Форумчанин

    Регистрация:
    10 июн 2014
    Сообщения:
    271
    Симпатии:
    159
    Сантьяго, два момента. Во-первых, дескриптор AUTHORITY["EPSG","msk26_2"] в конце необходимо удалить.

    Во-вторых, дескриптор TOWGS84 более не поддерживается стандартом OGC. Современный QGIS его тоже игнорирует. Как и многие другие программы, он использует ассоциированные с датумами трансформации. С датумом "Pulkovo 1995" по умолчанию используется трансформация "Pulkovo 1995 to WGS 84 (2)" по ГОСТ Р 51794-2008, что логично.
     
    #10
  11. Сантьяго

    Форумчанин

    Регистрация:
    9 янв 2008
    Сообщения:
    898
    Симпатии:
    354
    Адрес:
    Краснодар
    Спасибо, я этот prj делал для Фотоскана, или в фотоскан можно загружать систему координат в других форматах?
     
    #11
  12. ErnieBoyd

    Форумчанин

    Регистрация:
    10 июн 2014
    Сообщения:
    271
    Симпатии:
    159
    Извиняюсь, за TOWGS84 в Фотоскане не скажу.
     
    #12
  13. annaoafanaseva

    Регистрация:
    16 ноя 2021
    Сообщения:
    17
    Симпатии:
    0
    Добрый день!
    Подскажите, пожалуйста, верно ли указаны параметры местной системы координат в формате WKT для QGIS:

    BOUNDCRS[
    SOURCECRS[
    PROJCRS["unknown",
    BASEGEOGCRS["unknown",
    DATUM["Unknown based on Krassovsky, 1942 ellipsoid",
    ELLIPSOID["Krassovsky, 1942",6378245,298.3,
    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",90,
    ANGLEUNIT["degree",0.0174532925199433],
    ID["EPSG",8802]],
    PARAMETER["Scale factor at natural origin",1,
    SCALEUNIT["unity",1],
    ID["EPSG",8805]],
    PARAMETER["False easting",107543.3,
    LENGTHUNIT["metre",1],
    ID["EPSG",8806]],
    PARAMETER["False northing",-5540944.5,
    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["Position Vector transformation (geog2D domain)",
    ID["EPSG",9606]],
    PARAMETER["X-axis translation",23.92,
    ID["EPSG",8605]],
    PARAMETER["Y-axis translation",-141.27,
    ID["EPSG",8606]],
    PARAMETER["Z-axis translation",-80.9,
    ID["EPSG",8607]],
    PARAMETER["X-axis rotation",0,
    ID["EPSG",8608]],
    PARAMETER["Y-axis rotation",0.35,
    ID["EPSG",8609]],
    PARAMETER["Z-axis rotation",0.82,
    ID["EPSG",8610]],
    PARAMETER["Scale difference",0.99999988,
    ID["EPSG",8611]]]]

    Возникли сомнения, так как увидела разницу (смещение) в границах участка с сайта росреестра и на тойже подложке в qgis. Прикрепляю скрины экрана.
    Подскажите, пожалуйста, в чем проблема.
     

    Вложения:

    #13
  14. Сантьяго

    Форумчанин

    Регистрация:
    9 янв 2008
    Сообщения:
    898
    Симпатии:
    354
    Адрес:
    Краснодар
    Различаются знаки параметров вращения, поменяйте .
     
    #14
    annaoafanaseva нравится это.
  15. annaoafanaseva

    Регистрация:
    16 ноя 2021
    Сообщения:
    17
    Симпатии:
    0
    Спасибо! Можно мне подробнее, как первоклашке =))
    --- Сообщения объединены, 16 ноя 2021, Оригинальное время сообщения: 16 ноя 2021 ---
    Какие должны быть значения?
    --- Сообщения объединены, 16 ноя 2021 ---
    Вот эти значкния изменила, ничего не изменилось.
     

    Вложения:

    #15
  16. Сантьяго

    Форумчанин

    Регистрация:
    9 янв 2008
    Сообщения:
    898
    Симпатии:
    354
    Адрес:
    Краснодар
    -0,35,
    -0,82
    попробуйте так
     
    #16
    annaoafanaseva нравится это.
  17. annaoafanaseva

    Регистрация:
    16 ноя 2021
    Сообщения:
    17
    Симпатии:
    0
    Спасибо!
    --- Сообщения объединены, 16 ноя 2021, Оригинальное время сообщения: 16 ноя 2021 ---
    НЕ ПОМОГЛО =(
     
    #17
  18. Сантьяго

    Форумчанин

    Регистрация:
    9 янв 2008
    Сообщения:
    898
    Симпатии:
    354
    Адрес:
    Краснодар
    смотрите в личке
     
    #18
  19. annaoafanaseva

    Регистрация:
    16 ноя 2021
    Сообщения:
    17
    Симпатии:
    0
    Поменяла на эти:

    +proj=tmerc +lat_0=0 +lon_0=90 +k=1 +x_0=107543.3 +y_0=-5540944.5 +ellps=krass +towgs84=23.57,-140.95,-79.8,0,0.35,0.79,-0.22 +units=m +no_defs

    Все равно не совпадают границы.
     
    #19
  20. cfifti

    Форумчанин

    Регистрация:
    2 май 2012
    Сообщения:
    113
    Симпатии:
    44
    И не будут совпадать, так как на публичке задано смещение, по каждому из регионов свое.
     
    #20

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

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