Вы знаете, что кратчайший путь между двумя географическими точками - это большой круг. Моряки и авиаторы заинтересованы в минимизации пройденного расстояния и, следовательно, прошедшего времени. Вы также знаете, что румбная линия - это путь постоянного курса, естественное средство передвижения. В общем, чтобы пройти по великому круговому пути, нужно было бы непрерывно менять курс. Это нецелесообразно. Однако можно аппроксимировать большую траекторию окружности с помощью отрезков римбованной линии так, чтобы добавленное расстояние было меньшим, а количество курсов менялось минимальным.
Удивительно, но для близкой аппроксимации расстояния пути большой окружности требуется очень мало трапа.
Рассмотрим путешествие из Норфолка, штат Виргиния (37 ° N, 76 ° W), на мыс Сент - Винсент, Португалия (37 ° N, 9 ° W), один из наиболее интенсивных маршрутов в Атлантике. Протяжённость трассы с вершиной на восток составляет 3 213 морских миль, в то время как оптимальное большое расстояние от круга составляет 3 141 морскую милю.
Хотя длина пути линии румба лишь немногим более 2%, это дополнительные 72 мили по ходу путешествия. Для 12-узлового танкера это приводит к 6-часовой задержке, а в отгрузке время - деньги. Если для аппроксимации большого круга используются только три отрезка линии румба, общее расстояние поездки составляет 3 147 морских миль. Наш танкер потерпел бы лишь получасовую задержку по сравнению с непрерывным курсом румба. Вот код для вычисления трех типов треков между Норфолком и Сент-Винсентом:
figure('color','w'); ha = axesm('mapproj','mercator',... 'maplatlim',[25 55],'maplonlim',[-80 0]); axis off, gridm on, framem on; setm(ha,'MLineLocation',15,'PLineLocation',15); mlabel on, plabel on; load coastlines; hg = geoshow(coastlat,coastlon,'displaytype','line','color','b'); % Define point locs for Norfolk, VA and St. Vincent, Portugal norfolk = [37,-76]; stvincent = [37, -9]; geoshow(norfolk(1),norfolk(2),'DisplayType','point',... 'markeredgecolor','k','markerfacecolor','k','marker','o') text(-0.61,0.66,'Norfolk','HorizontalAlignment','left') geoshow(stvincent(1),stvincent(2),'DisplayType','point',... 'markeredgecolor','k','markerfacecolor','k','marker','o') text(0.54,0.66,'St. Vincent','HorizontalAlignment','right') % Compute and draw 100 points for great circle gcpts = track2('gc',norfolk(1),norfolk(2),... stvincent(1),stvincent(2)); geoshow(gcpts(:,1),gcpts(:,2),'DisplayType','line',... 'color','red','linestyle','--') text(-0.02,0.85,'Great circle: 3,141 nm (optimal)',... 'color','r','HorizontalAlignment','center') % Compute and draw 100 points for rhumb line rhpts = track2('rh',norfolk(1),norfolk(2),... stvincent(1),stvincent(2)); geoshow(rhpts(:,1),rhpts(:,2),'DisplayType','line',... 'color',[.7 .1 0],'linestyle','-.') text(-0.02,0.66,'Direct course: 3,213 nm',... 'color',[.7 .1 0],'HorizontalAlignment','center') % Compute and draw path along three waypoints [latpts,lonpts] = gcwaypts(norfolk(1),norfolk(2),... stvincent(1),stvincent(2),3); geoshow(latpts,lonpts,'DisplayType','line',... 'color',[.4 .2 0],'linestyle','-') text(-0.02,0.75,'3-leg approximation: 3,149 nm',... 'color',[.4 .2 0],'HorizontalAlignment','center')
Функция Toolbox™ мэппинга gcwaypts вычисляет ППМ в формате навигационной дорожки для аппроксимации большого круга с отрезками линии румба. Он использует следующий синтаксис:
[latpts,lonpts] = gcwaypts(lat1,lon1,lat2,lon2,numlegs)
Все входы для этой функции являются скалярами a (начальная и конечная позиция). numlegs введите требуемое количество ветвей равной длины, которое по умолчанию равно 10. Выходами являются векторы столбцов, представляющие ППМ в формате навигационной дорожки ([heading distance]). Размер каждого из этих векторов равен [(numlegs+1) 1]. Вот пункты для этого примера:
[latpts,lonpts] = gcwaypts(norfolk(1),norfolk(2),... stvincent(1),stvincent(2),3) % Compute 3 waypoints
latpts = 37.0000 41.5076 41.5076 37.0000 lonpts = -76.0000 -54.1777 -30.8223 -9.0000
Эти точки представляют ППМ вдоль большой окружности, между которыми аппроксимирующая траектория следует по осенним линиям. Четыре очка нужны для трех ног, потому что финальная точка на мысе Сент-Винсент должна быть включена.
Теперь мы можем вычислить расстояние в морских милях (нм) вдоль каждой дорожки и через ППМ:
drh = distance('rh',norfolk,stvincent); % Get rhumb line dist (deg) dgc = distance('gc',norfolk,stvincent); % Get gt. circle dist (deg) % Compute headings and distances for the waypoint legs [course distnm] = legs(latpts,lonpts,'rh');
distrhnm = deg2nm(drh) % Nautical mi along rhumb line distgcnm = deg2nm(dgc) % Nautical mi along great circle distlegsnm = sum(distnm) % Total dist along the 3 legs rhgcdiff = distrhnm - distgcnm % Excess rhumb line distance trgcdiff = distlegsnm - distgcnm % Excess distance along legs
distrhnm =
3.2127e+003
distgcnm =
3.1407e+003
distlegsnm =
3.1490e+003
rhgcdiff =
71.9980
trgcdiff =
8.3446