Этот пример показывает, как выполнить те же вычисления проекции, которые сделаны в рамках команд отображения 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