С января 2021 года в РФ обязательна к использованию система координат ГСК-2011 (ГОСТ 32453-2017). При этом, в большинстве софта с поддержкой систем координат данная система координат или отсутствует вообще, или определена с ошибками. Частый пример такой ошибки - использование нулевых параметров трансформации для пересчета на WGS-1984, основываясь на том, что отличия датумов WGS-1984 и ГСК-2011 незначительны. Как, например, сделано в некоторых сборках QGIS. Столкнувшись с этой проблемой при работе в QGIS и Petrel, был найден способ добавления туда корректных параметров ГСК-2011 в соответствии с ГОСТ.
Для QGIS данный вопрос очень хорошо рассмотрен на форуме сайта geodesist.ru geodesist.ru/...6/.
Но можно сделать чуть проще и добавить в дефолтный профиль эллипсоид GSK-2011 / GOST. Строка WKT для добавления кастомной проекции можно найти ниже.
А вот для Petrel пришлось написать небольшую инструкцию, которую также можно найти ниже.
Полезные ссылки:
ГОСТ 32453-2017
Приказ Росреестра от 12.11.2021 N П/0516
GSK2011 GOST WKT:
BOUNDCRS[
SOURCECRS[
GEOGCRS["GSK2011 GOST",
DATUM["GSK2011 GOST 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]],
CS[ellipsoidal,2],
AXIS["longitude",east,
ORDER[1],
ANGLEUNIT["degree",0.0174532925199433,
ID["EPSG",9122]]],
AXIS["latitude",north,
ORDER[2],
ANGLEUNIT["degree",0.0174532925199433,
ID["EPSG",9122]]]]],
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 GSK2011 GOST 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]]]]
Как добавить ГСК-2011 в Petrel
------------------------------------------
Шаг 1. добавить эллипсоид ГСК-2011:
Открыть Petrel Coordinate System Manager.
Вкладка Coordinate Reference Systems and Transforms.
Раздел Coordinate Reference Systems. Кликнуть зеленую кнопку с плюсом (Create new coordinate reference system).
В открывшемся окне вставить этот текст:
GEOGCS["GCS_GSK_2011",DATUM["D_GSK_2011",SPHEROID["GSK_2011",6378136.5,298.2564151]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433],AUTHORITY["CUSTOM",8901]]
------------------------------------------
Шаг 2. добавить трансформацию из ГСК-2011 в WGS-1984 в соответствии с ГОСТ 32453-2017:
Petrel Coordinate System Manager. Вкладка Coordinate Reference Systems and Transforms.
Раздел Transforms. Кликнуть зеленую кнопку с плюсом (Create new transform).
В открывшемся окне вставить этот текст:
GEOGTRAN["GSK_2011_To_WGS_1984_GOST",GEOGCS["GCS_GSK_2011",DATUM["D_GSK_2011",SPHEROID["GSK_2011",6378136.5,298.2564151]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],METHOD["Coordinate_Frame"],PARAMETER["X_Axis_Translation",0.013],PARAMETER["Y_Axis_Translation",-0.092],PARAMETER["Z_Axis_Translation",-0.03],PARAMETER["X_Axis_Rotation",0.001738],PARAMETER["Y_Axis_Rotation",-0.003559],PARAMETER["Z_Axis_Rotation",0.004263],PARAMETER["Scale_Difference",1.0000000074],AUTHORITY["CUSTOM",8611]]
------------------------------------------
Шаг 3. Создать проекцию Гаусса-Крюгера зона 8 на эллипсоиде ГСК-2011:
(в данном случае зона 8 используется как пример. можно создавать любую зону, меняя название и центральный меридиан Central_Meridian)
Coordinate System Manager. Вкладка Coordinate Reference Systems and Transforms.
Раздел Corrdinate Reference Systems. Кликнуть зеленую кнопку с плюсом (Create new coordinate reference system).
В открывшемся окне вставить этот текст:
PROJCS["GSK_2011_Gauss_Kruger_Zone_8N",GEOGCS["GCS_GSK_2011",DATUM["D_GSK_2011",SPHEROID["GSK_2011",6378136.5,298.2564151]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433],AUTHORITY["CUSTOM",8901]],PROJECTION["Gauss_Kruger"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",45.0],PARAMETER["Scale_Factor",1.0],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0],AUTHORITY["CUSTOM",100001]]
------------------------------------------
Шаг 4. Создать Conflation Policy для новой проекции, то есть связать ее с параметрами трансформации, указанными на шаге 2:
(в данном случае зона 8 используется как пример, вместо нее может быть любая другая)
Coordinate System Manager. Вкладка Geodetic Catalog. Кликнуть зеленую кнопку с плюсом (Create new conflation Policy).
Появится окно Create Conflation Policy. Указать в нем следующие данные:
Name: GSK_2011_Gauss_Kruger_Zone_8N_GOST
Authority: CUSTOM
Description: Gauss-Kruger zone 8 projection using GSK-2011 GCS with transformation from GOST 32453-2017
В разделе Select Coordinate Reference System сначала сбросить фильтр, кликнув слева кнопку
Click here to clear all filter criteria.
Потом в столбце Name под ячейкой Contains (это ячейка фильтра) вписать GSK.
В отфильтрованном списке выделить мышкой созданную на шаге 3 проекцию GSK_2011_Gauss_Kruger_Zone_8N
В разделе Select Transform убрать галочку
Show transforms whose area of use overlap with coordinate reference system (Esri only).
Сбросить фильтр, кликнув слева кнопку Clicke here to clear all filter criteria.
Потом в столбце Name под ячейкой Contains (это ячейка фильтра) вписать GSK.
В отфильтрованном списке выделить мышкой созданную на шаге 2 трансформацию GSK_2011_To_WGS_1984_GOST
-------------------------------------------
После этого, при загрузке в Petrel данных в 8-й зоне Гаусса-Крюгера ГСК-2011, нужно будет указать систему координат данных
GSK_2011_Gauss_Kruger_Zone_8N_GOST, выбрав ее из списка.
(8-я зона в данном случае как пример)
-------------------------------------------
ВАЖНО:
В данном случае зона 8 используется как пример. Указанные здесь названия проекции, центральный меридиан и номер зоны (шаги 3-4) актуальны только для 8-й зоны.
GSK_2011_Gauss_Kruger_Zone_8N_GOST - это только 8-я зона проекции Гаусса-Крюгера (центральный меридиан 45 в.д.).
Если данные В ДРУГОЙ ЗОНЕ, то нужно будет повторить шаги 3-4 для создания соответствующей проекции и conflation, но с указанием другого названия и центрального меридиана.