Документация

Координаты проекта без осей карты

В этом примере показано, как выполнить те же расчеты проекции, которые сделаны в рамках команд отображения Mapping Toolbox путем вызова defaultm и mfwdtran функции.

Создайте пустую структуру проекции карты для Синусоидальной проекции, с помощью defaultm функция. Функция возвращает mstruct .

mstruct = defaultm('sinusoid');

Установите пределы карты для mstruct . Чтобы заполнить поля структуры проекции карты и гарантировать эффекты настроек свойства правильно реализованы, вызывают defaultm во второй раз. Обратите внимание на то, что долгота источника сосредоточена между пределами долготы.

mstruct.maplonlimit = [-150 -30];
mstruct.geoid = referenceEllipsoid('grs80','kilometers');
mstruct = defaultm(mstruct);
mstruct.origin

ans = 1×3

     0   -90     0

Обрежьте побережье к пределам карты.

load coastlines
[latt,lont] = maptriml(coastlat,coastlon, ...
     mstruct.maplatlimit,mstruct.maplonlimit);

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

[x,y] = mfwdtran(mstruct,latt,lont);
figure
plot(x,y)
axis equal