Визуализация пространственной ошибки с помощью индикаторов ткани
Стандартным методом визуализации искажений, вносимых проекцией карты, является отображение небольших кругов с регулярными интервалами по всему земному шару. После проекции небольшие круги появляются как эллипсы различных размеров, удлинений и ориентаций. Размеры и формы эллипсов отражают проекционные искажения. Конформные выступы имеют круговые эллипсы, а равные по площади выступы имеют эллипсы одинаковой площади. Этот метод был изобретён Николя Тиссо в XIX веке, и эллипсы в его честь называются указателями Тиссо. Мера является тензорной функцией местоположения, которая изменяется в зависимости от места и отражает тот факт, что, если карта не является конформной, шкала карты отличается в каждом направлении в местоположении.
Визуализируйте проекционные искажения с помощью индикаторов Tissot
В этом примере показано, как добавить индикаторы Tissot к отображению карты.
Установите Синусоидальную проекцию в перекосном аспекте, построив график графика.
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