
- •Федеральное агентство морского и речного транспорта
- •Определение искажений в заданной точке карты
- •2. Графоаналитический способ определения m, n и
- •3. Определение элементов эллипса искажений
- •4. Построение эллипса искажений
- •5. Анализ искажений в заданной точке карты
- •5.1. Анализ искажения длин
- •5.2. Анализ искажения площадей
- •5.3. Анализ искажения направлений
- •5.4. Анализ искажения углов
- •Литература
2. Графоаналитический способ определения m, n и
Если для анализируемой карты не известны уравнения проекции, то значения m, n и можно определить графоаналитическим способом.
В
ариант
1.
В данной точке карты с географическими
координатами (,
)
строится угол
как острый угол между касательными к
изображению меридиана и параллели,
который затем измеряется транспортиром
(рис.1).
Р
ис.1
Для определения m и n достаточно измерить отрезки изображения меридиана l' и параллели k' между ближайшими линиями картографической сетки 1 и 2 (для l'), 1 и 2 (для k'). Измерения можно выполнять при помощи курвиметра или способом Шокальского с использованием циркуля-измерителя.
Значения m и n вычисляются по формулам:
m = = ;
(2)
n = = .
где l = ; k = соответственно длины отрезков дуг меридиана параллели в натуре (на сфероиде);
N0 знаменатель главного масштаба карты ( 0 = 1: N0);
L длина дуги меридиана на земном сфероиде между параллелями с широтами 1 и 2;
K длина дуги параллели на земном сфероиде между меридианами с долготами 1 и 2.
Значение L может быть определено как разность длин отрезков меридиана на земном сфероиде от экватора до параллелей с широтами 1 и 2:
L = X2 – X1 , (3)
где
X1 = ;
X2 = .
Здесь M , как и в выражениях (1), радиус кривизны меридиана в точке с заданной широтой. Значения X1 и X2 могут быть выбраны из "Картографических таблиц" [ 7 ] по аргументам 1 и 2. Значение
K = (r arc1') (2 - 1)', (4)
где (r arc1') - значение длины одной минуты параллели с широтой, выбирается из "Картографических таблиц";
(2 - 1 )' разность долгот, выраженная в угловых минутах.
Пример. На карте с 0 = 1 : 20 000 000 в точке с координатами = 60° с.ш., = 30° в.д. измерены величины:
l' = 107,7 мм ,
k' = 60,6 мм ,
= -72,8° (sin = 0,95528).
Из "Картографических таблиц" выбраны значения:
X1 = 5 540 944 м ,
X2 = 7 769 116 м ,
r arc1' = 930 015 мм/мин ,
с учётом которых по формулам (2)(4) получаем:
L = 2 228 172 м = 2 228,17 106 мм ;
K = 930 015 мм/мин1200'= 1 116,02 106 мм ;
m = = = 0,9667 ;
n = = = 1,0860 .
Вариант 2. В исследуемой точке карты с географическими координатами (, ) задаётся прямоугольная система координат. В этой системе координат путём измерения определяются (с учётом знака) приращения координат соответственно l'x, l'y для изменения широты от значения 1 до значения 2, и k'x, k'y для изменения долготы со значения 1 до значения 2 (рис. 2).
Значения m и n, которые вычисляются по формулам
m = ;
(5)
n = ,
Рис. 2
где
mx = ; nx = ;
(6)
my = ; ny =
являются проекциями относительных масштабов по меридианам и параллелям m и n на оси X и Y заданной системы прямоугольных координат. Значения l и k определяются так же как в варианте 1.
Значение угла в данном случае вычисляется по формуле
= arctg() , (7)
где в числителе p = mxny - mynx относительный масштаб площади.
Пример. На карте с 0 = 1 : 20 000 000 в точке с координатами = 60° с.ш., = 30° в.д. измерены величины:
l'x = +95,9 мм , k'x = + 9,8 мм ,
l'y = -48,8 мм , k'y = +60,0 мм .
Из "Картографических таблиц" выбраны значения:
X1 = 5 540 944 м ,
Y1 = 7 769 116 м ,
r=60° arc1' = 930 015 мм/мин ,
с учётом которых по формулам (3) и (4) получаем:
L = 2 228 172 м = 2228,17 106 мм ;
K = 930 015 мм/мин1200'= 1116,02 106 мм ;
Откуда
l = = 111,40 мм ;
k = = 55,80 мм ;
По формулам (5)(7) вычисляем:
mx = = +0,8609; nx = = +0,1756;
my = = 0,4381; ny = = +1,0753;
m = = 0,9659 ;
n = = 1,0895 ;
p = 0,86091,0753 (0,4381)0,1756 = 1,0026 ;
= arctg() = 72° 18,4'.
Рассматриваемый пример соответствует приведённому выше примеру аналитического определения m, n и . Результаты графоаналитического способа хорошо согласуются с аналитическим решением.