MATLAB.Exponenta
–Û·Ë͇ Matlab&Toolboxes

Математика\Mapping Toolbox

В оглавление \ К предыдущему разделу \ К следующему разделу

Основные понятия сферической географии. Виды проекций

Земля в первом приближении напоминает сферу. Такая модель может применяться только для приближенных географических расчетов. На самом деле Земля является слегка приплюснутой сферой – эллипсоидом. Такое модельное представление является более точным.

Героид – это эмпирическое приближение к Земле без учета особенностей картографического рельефа. Для построения героида используется следующий набор команд.

clear;load geroid;load coast
figure;axesm robinson
meshm(geoid,geoidlegend)
colorbar('horiz');
plotm(lat,long,'k')

Как было сказано выше, при расчетах повышенной точности Земля представляется в виде эллипсоида и описывается трема основными параметрами: ось вращения, малая полярная ось и большая экваториальная ось (экваториальный радиус). Вектором эллипсоида в Mapping Toolbox называется вектор, содержащий значение экваториального радиуса и отклонение от него полярного радиуса.

Рассмотрим пример вычисления такого радиуса:

almanac('earth','ellipsoid','kilometers')
ans =
  1.0e+003 *
    6.3781    0.0001

Этот пример относится к мировой геодезической системе 1972 г. Такие же параметры можно найти и для других геодезических систем. Они немного отличаются.

almanac('earth','sphere','kilometers')
ans =
        6371           0

Второй параметр в этом случае рамен нулю, поскольку экваториальный и полярный радиусы для сферы одинаковы.

Основные виды проекций

Земной шар – это трехмерный объект. Однако карты, как правило, объекты плоские. Возникает вопрос – как наиболее оптимально с точки зрения точности спроецировать Земной шар на поверхность. Пакет Mapping Toolbox включает более 60 наиболее широко распространенных проекций.

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

В качестве примера рассмотрим построение карты с синусоидальной проекцией.

close all;clear all; axesm sinusoid
grid on; framem on;
setm(gca,'Origin',[20 30 45]);
load coast; patchm(lat,long,'g');tissot

 

Последняя команда tissot наносит на карту множество окружностей, преобразуемых в соответствии с видом проекции. Это позволяет оценить степень геометрических искажений в каждой области карты. Рассмотрим еще несколько видов различных проекций, которые доступны в Mapping Toolbox.

Цилиндрические проекции

Проекция Balthasart

Проекция Balthasart была представлена в 1935 году и является специальной формой равноплощадных цилиндрических проекций.

landareas = shaperead('landareas.shp','UseGeoCoords',true);
axesm ('balthsrt', 'Frame', 'on', 'Grid', 'on');
geoshow(landareas,'FaceColor',[1 1 .5],'EdgeColor',[.6 .6 .6]);
tissot;

 

Проекция Behrmann

Эта проекция названа именем Вальтера Бермана и была представлена в 1910 году.

landareas = shaperead('landareas.shp','UseGeoCoords',true);
axesm ('behrmann', 'Frame', 'on', 'Grid', 'on');
geoshow(landareas,'FaceColor',[1 1 .5],'EdgeColor',[.6 .6 .6]);
tissot;

 

Проекция Bolshoi Sovietskii Atlas Mira Projection

Эта проекция была впервые представлена в 1937 году.

landareas = shaperead('landareas.shp','UseGeoCoords',true);
axesm ('bsam', 'Frame', 'on', 'Grid', 'on');
geoshow(landareas,'FaceColor',[1 1 .5],'EdgeColor',[.6 .6 .6]);
tissot;

 

Проекция Braun

Эта проекция была впервые представлена в 1867 году.

landareas = shaperead('landareas.shp','UseGeoCoords',true);
axesm ('braun', 'Frame', 'on', 'Grid', 'on');
geoshow(landareas,'FaceColor',[1 1 .5],'EdgeColor',[.6 .6 .6]);
tissot;

 

Проекция Cassini

Эта проекция разработана Cesar Francois Cassini de Thury (1714-84).

landareas = shaperead('landareas.shp','UseGeoCoords',true);
axesm ('cassini', 'Frame', 'on', 'Grid', 'on');
geoshow(landareas,'FaceColor',[1 1 .5],'EdgeColor',[.6 .6 .6]);
tissot;

 

Проекция Cassini (стандарт)

Эта проекция разработана Cesar Francois Cassini de Thury (1714-84).

landareas = shaperead('landareas.shp','UseGeoCoords',true);
axesm ('cassinistd', 'Frame', 'on', 'Grid', 'on');
geoshow(landareas,'FaceColor',[1 1 .5],'EdgeColor',[.6 .6 .6]);
tissot;

 

Центральная цилиндрическая проекция

Эта проекция доступна только на сфере. Данные, которые представлены на широтах больше чем 750, являются сильно деформированными.

landareas = shaperead('landareas.shp','UseGeoCoords',true);
axesm ('ccylin', 'Frame', 'on', 'Grid', 'on');
geoshow(landareas,'FaceColor',[1 1 .5],'EdgeColor',[.6 .6 .6]);
tissot;

 


В оглавление \ К предыдущему разделу \ К следующему разделу

 


Поиск по сайту:

Система Orphus

Яндекс.Метрика