Параметры искажения для проекций карты
areascale = distortcalc(lat,long)
areascale = distortcalc(mstruct,lat,long)
[areascale,angdef,maxscale,minscale,merscale,parscale]
= distortcalc(...)
areascale = distortcalc(lat,long) вычисляет искажение области для текущей проекции карты в заданном географическом положении. Шкала области 1 не указывает ни на какое искажение шкалы. Широта и долгота могут быть скалярами, векторами или матрицами в угловых модулях заданной проекции карты.
areascale = distortcalc(mstruct,lat,long) использует проекцию, заданную в структуре карты mstruct.
[areascale,angdef,maxscale,minscale,merscale,parscale]
= distortcalc(...) вычисляет шкалу области, максимальную угловую деформацию прямых углов (в угловых модулях заданной проекции), конкретные максимальные и минимальные искажения шкалы в любом направлении и конкретная шкала вдоль меридиана и параллели. Можно также вызвать distortcalc с меньшим количеством выходных аргументов, в показанном порядке.
Сопоставьте проекции, неизбежно вводят искажения в формах и размерах объектов, когда они преобразовываются от 3D сферических координат до двумерных Декартовых координат. Сумма и тип искажения варьируются между проекциями по проекции, и с выбором параметров проекции, такими как стандартные параллели. Эта функция позволяет количественную оценку параметров искажения.
На экватор Меркаторская проекция свободна и от области и от углового искажения:
axesm mercator
[areascale,angdef] = distortcalc(0,0)
areascale =
1.0000
angdef =
8.5377e-007В 60 градусах на север объекты показывают в 400% их истинной области. Проекция конформна, таким образом, угловое искажение является все еще нулевым.
[areascale,angdef] = distortcalc(60,0)
areascale =
4.0000
angdef =
4.9720e-004Эта функция использует метод конечной разности. Географические координаты встревожены немного в различных направлениях и спроектированы. Небольшое количество ошибки введено численным расчетом производных и изменением параметров искажения карты.