Як створити точну Tissot Indicatrix?


31

Тіссо Індикатриса є корисним методом для передачі на першому погляді видів спотворень даної проекція схильна (на малюнку нижче, кожен з червоних кіл займає ту ж площу). Мені сказали, що у популярних методів для генерування ТІ є свої проблеми, які часом бувають жахливо неточними.

У чому полягає проблема з популярними методами та який найбільш правильний спосіб генерування TI, доступного для вашого середнього GIS-чувака (ette)?

Меркатор і глобуси з тисотами


2
велике запитання. Я також хотів би це знати.
Джордж Сільва

Відповіді:


23

Будь-яке програмне забезпечення, яке може точно проектувати координати, може обчислити точні показники Tissot .

Хорошим джерелом для формул є Snyder, John, Map Projections - Робочий посібник , головним чином, на стор. 20-26. (Я не відтворюю їх тут, оскільки на цьому веб-сайті немає відповідних інструментів для передачі математичних формул.) Вони потребують усіх чотирьох перших похідних проектованих координат (x, y) щодо сферичних координат (lat, lon) = (фі, лямбда):

dx / d(phi), dx / d(lambda);
dy / d(phi), dy / d(lambda).

Все інше про ТІ обчислюється з точки зору цих (використовуючи деякі арифметичні та тригонометричні функції: косинус, головний зворотний синус та головний зворотний дотичний). Для обчислень потрібен опис форми землі. Для найбільшої точності використовуйте еліпсоїдальну дату із напіввеликою віссю a та ексцентриситетом e. (Вони будуть відомі програмному забезпеченню.)

У книзі Снайдера є інструкції, як обчислити все, крім цих похідних. Зробіть це чисельно. Я мав відмінні результати, використовуючи оцінки центральних кінцевих різниць першого порядку на відстані h = 10 ^ (- 5,2) радіанів (як правило, близько 50 метрів): це хороший компроміс між спробою наблизитись до нескінченно малої кількості і втратити занадто велику точність від округлення з плаваючою точкою (якщо припустити подвійну точність), оскільки допущена помилка пропорційна (10 ^ (- 5.2)) ^ 2 = 10 ^ (- 10.4), а 10 ^ (- 5.2) дорівнює 10 ^ 10.4 кратній точності подвійної точності IEEE 10 ^ (- 15,6), і це все-таки набагато більше, ніж типова точність проекцій, які зазвичай проходять від 10 ^ (- 10) до приблизно 10 ^ (- 14).

Отже, як ви обчислюєте кінцеві оцінки різниці? Ця частина напрочуд проста. Щоб отримати dx / d (phi) у точці (phi, лямбда), попросіть ваш ГІС спроектувати бали

(phi - h/2, lambda) --> (x0,y0),
(phi + h/2, lambda) --> (x1,y1).

Використовуйте кошториси

dx / d(phi) = (x1 - x0)/h,
dy / d(phi) = (y1 - y0)/h.

Аналогічно спроектуйте бали

(phi, lambda - h/2) --> (x2,y2),
(phi, lambda + h/2) --> (x3,y3)

і використовувати оцінки

dx / d(lambda) = (x3 - x2)/h,
dy / d(lambda) = (y3 - y2)/h.

Це займає чотири проекції та невелику арифметику. (Ви можете зменшити його до трьох, використовуючи не центральні відмінності, але точність трохи знижується. Розумно прагнути до високої точності, не дозволяючи h отримати занадто малим, якщо ви не впевнені, що ваш ГІС використовує оцінку (міліметр) точність у формулах прогнозування.)

З цих похідних, поряд із формулами Снайдера (звертаючи увагу на модифікації, описані в 4-19 та 4-21), можна отримати довжину осей індикатрії Tissot at (phi, лямбда) та її орієнтацію. На картах світового масштабу TI буде настільки малим, щоб бути непомітним, тому остаточне, що потрібно зробити, - це вирішити, скільки ви хочете змінити масштаб кожного TI. Я визначаю коефіцієнт масштабу, з’ясовуючи, наскільки великою буде карта, знаходя розміри типових TI на карті та масштабуючи так, що ці ТІ будуть приблизно на 6% ширшими від карти. Все-таки вдалий початок; Я дозволяю користувачеві регулювати розмір TI звідти. Звичайно, ви будете масштабувати всі TI на однакову кількість, щоб їх можна порівняти, і кожен буде розміщений навколо власного центру (який отримується за п’ятою проекцією, (phi, лямбда) -> (x, y) ).

Приємним доповненням до еліптичного зображення ТІ є показ напрямків локального меридіана та паралелі: тоді, з першого погляду, можна оцінити конвергенцію сітки . Я також показую стандартне коло (не представляє спотворень), концентричне для кожного ТІ, оскільки воно покращує здатність читача оцінювати кількість спотворень, представлених кожним еліпсом.

alt текст

Слід зазначити, що ця проекція Mollweide - крайня TI біля південного полюса. Це все ще ідеальний еліпс і точно описує спотворення карти там.


2
whuber найкращий: P.
Джордж Сільва

1
Я помітив, що ESRI опублікував статтю про створення індикатриси tissot з буферами, чи правильний буферний метод, оскільки індикатрикс і буфер «зовсім не однакові»? blogs.esri.com/Support/blogs/mappingcenter/archive/2011/03/28/…
SaultDon

1
@Sault Буфери є заміною правильних TI. Загалом вони працюватимуть досить добре, за умови, що вони не простягнуться більше ніж на кілька градусів, за винятком випадків, коли викривлення стають крайніми, де або вони взагалі не працюватимуть (бо буфер розливається по області нескінченного спотворення, як показано на ця карта ESRI ) або вони призводять до нееліптичних фігур. Невелика зміна зробить такий підхід набагато кращим: обчисліть крихітні буфери, такі як буфери на 50 м, і розгорніть кожен рівномірно (на карті) навколо його центру, щоб зробити його видимим.
whuber

1
Як ви обчислили орієнтацію еліпсів?
Джейсон Дейвіс

3
@ Джейсон Покажчик Тіссота - це не що інше, як мальовниче зображення похідної проекції, прикладеної до кола в дотичній площині в точці. Найпростіший спосіб зробити це, то, щоб обчислити цю похідну (це 2 на 2 матриці ), параметризацію кола, і застосовувати похідну для параметрезованих вершин. Вуаля, у тебе еліпс. (Саме так намальовані еліпси на моїй фігурі.) Його осі можна обчислити як власні вектори A'A ; це простий розрахунок, що включає нічого складнішого, ніж квадратичну формулу.
whuber
Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.