Конформное преобразование

Конформное преобразование на плоскости широко используется в геодезии при создании местных координатных систем на небольшие территории, ограниченные, как правило, размерами населённого пункта.

Введение

Следующие формулы связывают координаты точек x, y, заданные в местной системе координат (МСК), и координаты X, Y, заданные в государственной системе координат (ГСК):

${\begin{array}{lcl}X&=&a_{0}+m(x\cos \theta -y\sin \theta )\\Y&=&b_{0}+m(x\sin \theta +y\cos \theta )\end{array}}$

где a₀, b₀ — положение начала МСК в ГСК, m — масштабный множитель, θ — угол разворота.

Используемый в геодезии набор параметров называется «ключ». Он определяется следующими величинами: X₀, Y₀, x₀, y₀, m, θ. Первые четыре — обычно координаты одного из геодезических пунктов в двух системах.

Исходный материал для определения параметров — пары координат пунктов геодезической сети, полученные из независимого уравнивания одних и тех же измерений в МСК и в ГСК. В зависимости от класса пунктам (вернее, парам уравнений) назначаются веса p.

Алгоритм нахождения параметров

Для определения четырёх параметров принимается следующая математическая модель:

${\begin{array}{lcl}X'&=&a_{0}+a_{1}x-b_{1}y\\Y'&=&b_{0}+b_{1}x+a_{1}y\end{array}}$

Очевидно, конформное преобразование является частным случаем аффинного.

1. Вычисление взвешенных средних

${\bar {x}}={\frac {\sum \limits _{i=1}^{n}p_{i}x_{i}}{\sum \limits _{i=1}^{n}p_{i}}}\quad {\bar {y}}={\frac {\sum \limits _{i=1}^{n}p_{i}y_{i}}{\sum \limits _{i=1}^{n}p_{i}}}\quad {\bar {X}}={\frac {\sum \limits _{i=1}^{n}p_{i}X_{i}}{\sum \limits _{i=1}^{n}p_{i}}}\quad {\bar {Y}}={\frac {\sum \limits _{i=1}^{n}p_{i}Y_{i}}{\sum \limits _{i=1}^{n}p_{i}}}$

2. Перенос осей в центр масс

$\Delta x_{i}=x_{i}-{\bar {x}}\quad \Delta y_{i}=y_{i}-{\bar {y}}\quad \Delta X_{i}=X_{i}-{\bar {X}}\quad \Delta Y_{i}=Y_{i}-{\bar {Y}}$

3. Вычисление a₁ и b₁

$S_{1}=\sum \limits _{i=1}^{n}p_{i}\Delta X_{i}\Delta x_{i}$	Невозможно разобрать выражение (SVG (MathML можно включить с помощью плагина для браузера): Недопустимый ответ («Math extension cannot connect to Restbase.») от сервера «https://wikimedia.org/api/rest_v1/»:): {\displaystyle S_2 = \sum\limits_{i=1}^n p_i \Delta Y_i \Delta y_i}
Невозможно разобрать выражение (SVG (MathML можно включить с помощью плагина для браузера): Недопустимый ответ («Math extension cannot connect to Restbase.») от сервера «https://wikimedia.org/api/rest_v1/»:): {\displaystyle S_3 = \sum\limits_{i=1}^n p_i \Delta Y_i \Delta x_i}	Невозможно разобрать выражение (SVG (MathML можно включить с помощью плагина для браузера): Недопустимый ответ («Math extension cannot connect to Restbase.») от сервера «https://wikimedia.org/api/rest_v1/»:): {\displaystyle S_4 = \sum\limits_{i=1}^n p_i \Delta X_i \Delta y_i}
Невозможно разобрать выражение (SVG (MathML можно включить с помощью плагина для браузера): Недопустимый ответ («Math extension cannot connect to Restbase.») от сервера «https://wikimedia.org/api/rest_v1/»:): {\displaystyle S_5 = \sum\limits_{i=1}^n p_i \left( \Delta x_i^2 + \Delta y_i^2 \right)}
Невозможно разобрать выражение (SVG (MathML можно включить с помощью плагина для браузера): Недопустимый ответ («Math extension cannot connect to Restbase.») от сервера «https://wikimedia.org/api/rest_v1/»:): {\displaystyle a_1 = \frac{S_1 + S_2}{S_5}}	Невозможно разобрать выражение (SVG (MathML можно включить с помощью плагина для браузера): Недопустимый ответ («Math extension cannot connect to Restbase.») от сервера «https://wikimedia.org/api/rest_v1/»:): {\displaystyle b_1 = \frac{S_3 - S_4}{S_5}}