Стандартный метод визуализации искажений, введенных проекцией карты, состоит в том, чтобы отобразить маленькие круги равномерно по всему миру. После проекции маленькие круги появляются как замещающие знаки различных размеров, удлинений и ориентаций. Размеры и формы замещающих знаков отражают искажения проекции. Конформные проекции имеют круговые замещающие знаки, в то время как проекции равной области имеют замещающие знаки той же области. Этот метод был изобретен Николя Тиссо в 19-м веке, и замещающие знаки называются Тиссо indicatrices в его честь. Мерой является функция тензора местоположения, которое варьируется с места на место и отражает то, что, если карта не конформна, шкала карты отличается в каждом направлении в местоположении.
В этом примере показано, как добавить Тиссо indicatrices в отображение карты.
Настройте Синусоидальную проекцию в скошенном аспекте, строя координатную сетку.
figure axesm sinusoid gridm on framem on setm(gca,'Origin',[20 30 45])
Загрузите набор данных побережья и постройте его как зеленые закрашенные фигуры.
load coastlines patchm(coastlat,coastlon,'g')
Постройте значение по умолчанию схема Тиссо. Заметьте, что круги значительно варьируются по форме. Это указывает, что Синусоидальная проекция не конформна. Несмотря на искажения, однако круги все покрытие равняется суммам области на карте, потому что проекция имеет свойство равной области. Значение по умолчанию схемы Тиссо чертятся с синими незаполненными кругами с 100 точками, расположило 30 градусов с интервалами независимо в области обоих направлений. Круговой радиус по умолчанию является 1/10 текущего радиуса ссылочного эллипсоида (по умолчанию, что радиус равняется 1).
tissot
Очистите схему Тиссо, вращайте проекцию к полярному аспекту и постройте новую схему Тиссо с помощью кругов, по которым следуют 20 градусов независимо, вдвое менее большой, чем прежде, чертивший с 20 точками и подошедшим к концу красным. В результате обратите внимание, что круги нарисованы быстрее, потому что меньше точек вычисляется для каждого. Также обратите внимание, что искажения являются все еще самыми маленькими близко к источнику карты и все еще самыми большими около системы координат карты.
clmo tissot setm(gca,'Origin',[90 0 45]) tissot([20 20 .05 20],'Color','r')
distortcalc | mdistort | tissot
