Параметры искажения для проекций карты
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
Эта функция использует метод конечной разности. Географические координаты встревожены немного в различных направлениях и спроектированы. Небольшое количество ошибки введено численным расчетом производных и изменением параметров искажения карты.