Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
МУ - Анализ искажений (редакция 2012 года) - А4....doc
Скачиваний:
5
Добавлен:
23.04.2019
Размер:
1.76 Mб
Скачать

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 = X2X1 , (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,86091,0753  (0,4381)0,1756 = 1,0026 ;

 = arctg() = 72° 18,4'.

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