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

Самой распространенной ошибкой перевода в "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])
и свои собственные для максимальной точности.

 

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

 

Давайте рассмотрим для примера 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 в сервисах ПКК
для перевода данных ЕГРН в сферический меркатор.

 

Вы ссылаетесь на ресурс 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 ---

Вот не сказал бы, глядя на ПКК и сравнивая объекты на ней с координатами из выписок.

 

Согласен. На официальном сайте 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. В других регионах возможно это не так,
приношу извинения за чрезмерное обобщение.

 

Сами параметры epsg.io выдаёт правильно. Посмотрите ту же EPSG:5044. Косяки только в экспортных примерах.

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

 

Здравствуйте
подскажите я правильно сделал 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"]]

 

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

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

 

Спасибо, я этот prj делал для Фотоскана, или в фотоскан можно загружать систему координат в других форматах?

 

Извиняюсь, за TOWGS84 в Фотоскане не скажу.

 

Добрый день!
Подскажите, пожалуйста, верно ли указаны параметры местной системы координат в формате 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. Прикрепляю скрины экрана.
Подскажите, пожалуйста, в чем проблема.

 

Различаются знаки параметров вращения, поменяйте .

 

Спасибо! Можно мне подробнее, как первоклашке =))

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

Какие должны быть значения?

--- Сообщения объединены, 16 ноя 2021 ---

Вот эти значкния изменила, ничего не изменилось.

 

-0,35,
-0,82
попробуйте так

 

Спасибо!

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

НЕ ПОМОГЛО =(

 

смотрите в личке

 

Поменяла на эти:

+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

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

 

И не будут совпадать, так как на публичке задано смещение, по каждому из регионов свое.