Локальная система координат в формате proj4

[Эльдар]

Уважаемые ГИСовцы на этот раз хотелось бы попробовать создать локальную систему координат в формате PROJ4.
есть параметры перехода.
lat_0=0 lon_0=46.8 x_0=350000.39 y_0=-4542621.65 и еще есть угол поворота локальной СК Angle="45.89"
по этим параметрам получена следующая проекция
+proj=tmerc +lat_0=0 +lon_0=46.8 +k=1 +x_0=350000.39 +y_0=-4542621.65 +ellps=krass +towgs84=24.83,-130.97,-81.74,0,0,0.13,-0.22 +units=m +no_defs
Подскажите пожалуйста как и куда прописать этот параметр Angle="45.89"
Проекция будет использоваться в QGIS

[Игорь Белов]

Это для Вас написано.
Кратко суть: для tmerc нет параметра вращения, что плохо.
Территория под локальной СК ограничена, что хорошо, так как позволяет использовать другую конформную проекцию без значительных искажений.
Выбор проекций невелик: omerc. Параметры придётся посчитать переносом с экватора на центр территории.

[Эльдар]

Полезная статья, попробую разобраться

[Игорь Белов]

Давайте вернёмся к исходной теме. Вам нужно решить свою задачу, а не обсуждать статью.

Начать нужно с того, чем Вы располагаете. В данном случае одной строки параметров PROJ.4 недостаточно. В идеале необходимо знать координаты центра локальной проекции, а также иметь надёжные пары координат в обеих системах. Если координаты центра неизвестны, нужно хотя бы иметь представление, какой географический объект обеспечивается локальной координатной системой.

Итак, какими данными Вы располагаете?

[Эльдар]

Есть следующие параметры.
координаты центра B 42,985472 L 47.500106
угол поворота 3,253402778
и 2 точки с известными координатами.

1
x=-7010.47577 B42.918967
y=6771.15086 L47.578003

2
x=-3900.46336 B42.947343
y=5949.95408 L47.570166

[Игорь Белов]

Очень хорошо, вот строка параметров:

Код: Выделить всё

+proj=omerc +lat_0=42.985472 +lonc=47.500106 +alpha=3.253402778 +gamma=0 +k_0=1 +x_0=0 +y_0=0 +ellps=krass

Пересчитав географические координаты с этими параметрами, получаем прямоугольные:

Код: Выделить всё

1. 6769.126735 -7012.483760
2. 5948.777730 -3902.204889

Отличия от исходных прямоугольных координат порядка двух метров.

Если усомниться в надёжности значения угла разворота и параметров k_0, x_0, y_0, можно получить нулевые невязки при

Код: Выделить всё

+alpha=3.26921513 +k_0=0.9999872043 +x_0=0.176 +y_0=0.050

Косвенное подтверждение справедливости такого подхода — получающиеся малые значения параметров x_0, y_0. Однако исходных данных для такой математики маловато, нужно больше точек, расположенных в разных направлениях от центра проекции.

[Эльдар]

как я понимаю надо определить a0 и a1 как тут, а потом перенести значения к x_0, y_0

[Игорь Белов]

a₀ → x_0
b₀ → y_0
m → k_0
θ → alpha

[Эльдар]

благодарен за помощь))))