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

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

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

Создайте пустую структуру проекции карты для Синусоидальной проекции, с помощью функции 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