Форумы GIS-Lab.info

Уважаемые ГИСовцы на этот раз хотелось бы попробовать создать локальную систему координат в формате 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

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

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

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

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

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