2.22 ГСК 2011 - настройка на примере Global mapper и Qgis
1. Global mapper
При попытке открыть файл в неизвестной системе координат, Global Mapper выводит стандартное окно выбора проекции. Давайте заведем параметры для ГСК-2011.
Нажимаем Add datum, а в следующем окне нажимаем Add ellipsoid:
Нажимаем ОК, на базе созданного эллипсоида создаем новый датум со следующими параметрами:
Параметры:
Эллипсоид: GSK-2011
Параметры:
Большая полуось=6378136.5 (м)
Малая полуось=6356751.758 (м)
Геометрическое сжатие: 298,2564151 (м) - данный параметр в рамках статьи не используется
Датум (имя):GSK-2011_Gauss-Kruger
7 параметров, метод - coordinate frame rotation
dX= 0.013
dY= -0.092
dZ= -0.03
Rotation (arc-ceconds)
X= 0.001738
Y= -0.003559
Z= 0.004263
Scale(ppm)= 1.0000000074
Далее используем либо проекцию Гаусса-Крюгера (6-ти градусные зоны)
Наша зона - № 11
Либо Transverse Mercator со следующими дополнительными параметрами:
Scale factor = 1
Central meridian = 63.000000
Origin latitude = 0
False easting(m) = 11500000
False northing(m) = 0
Rotation angle = 0
2.Qgis
Все в некотором роде проще:
В окно (4) вводим следующий текст: (форматирование похерилось, но на работоспособность это не повлияло)
BOUNDCRS[
SOURCECRS[
PROJCRS["GSK-2011 / Gauss-Kruger zone 11",
BASEGEOGCRS["Unknown datum based upon the GSK-2011 ellipsoid",
DATUM["Not specified (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["6-degree Gauss-Kruger zone 11",
METHOD["Transverse Mercator",
ID["EPSG",9807]],
PARAMETER["Latitude of natural origin",0,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8801]],
PARAMETER["Longitude of natural origin",63,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8802]],
PARAMETER["Scale factor at natural origin",1,
SCALEUNIT["unity",1],
ID["EPSG",8805]],
PARAMETER["False easting",11500000,
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]]]]
Проверяем:
Если все ок, можем прямо в этом окошке рассчитать образец координат (5).