- Картографические проекции. Урок 6
- Виды искажений при использовании картографических проекций
- Типы проекций по характеру искажений
- Равноугольные картографические проекции
- Равновеликие картографические проекции
- Произвольные картографические проекции
- Равнопромежуточные картографические проекции
- Классификация географических проекций по геометрической фигуре, являющейся вспомогательной поверхностью
- Ликбез по картографическим проекциям с картинками
- Легкое введение
- Проблема
- Варианты решения
- Немного экзотики и специальных случаев
Картографические проекции. Урок 6
Картографические проекции сегодня – это математические способы изображения всего земного эллипсоида или его части на плоскости, систематическое преобразование широт и долгот с поверхности сферы на плоскость.
Для создания географических карт выполняют две последовательных операции:
- Перенесение (проектирование) поверхности геоида с его сложным рельефом на поверхность эллипсоида вращения или шара.
- Дальнейшее его проектирование на плоскость (преобразование системы географических координат в декартовы) с использованием одной из картографических проекций.
При этом картографы пытаются добиться как можно меньшего количества искажений. Сделать мелкомасштабную карту совсем без искажений невозможно. На крупномасштабных (топографических) картах искажения почти отсутствуют. В зависимости от назначения карты одни погрешности допустимы, другие нет. Поэтому и существуют разные типы проекций, предназначенные для сохранения некоторых свойств сферы за счёт других её свойств.
Проекция на шар — глобус. Автор: UBC Library Digitization Centre
Виды искажений при использовании картографических проекций
Разложить на плоскости эллипс или шар очень трудно, для того, чтобы убедиться в этом, можно попробовать это сделать на практике. Сложить кусочки апельсиновой кожуры так, чтобы между ними не было пустых мест и попробовать получить непрерывную ровную плоскость. Корка соберётся в складки, она не уложится без промежутков.
При любом способе разложения шара на плоскость присутствует один или несколько типов искажения:
- базовое – искажение расстояний (длин линий), от него зависит степень других видов деформаций. Признак: между соседними параллелями отрезки меридианов неодинаковы по длине;
- площадей. При таком искажении между соседними параллелями форма и величина (а значит и площадь) ячеек неодинакова;
- углов – углы между определённым направлением на местности и на карте не совпадают. Узнать его можно по тому, что углы между параллелями и меридианами не являются прямыми;
- форм. При одинаковой площади форма клеток, находящихся на одной широте, разная.
При этом типы искажений взаимозависимы, при уменьшении одного из показателей увеличивается другой. В зависимости от назначения карты, на ней присутствуют места с нулевым искажением, с удалением от него количество искажений увеличивается. Поэтому на карте есть три вида масштаба:
- основной (тот, что подписан), действующий на линии нулевого искажения,
- частные (определяются при помощи эллипса искажений), их может быть бесконечно много;
- средний (совокупность частных масштабов отрезка).
При выборе типа картографической проекции сначала строят изоколы – изолинии, соединяющие точки с одинаковым искажением.
Изоколы искажения углов
Источник: https://ds04.infourok.ru/uploads/ex/0617/00148bfe-04623ef1/hello_html_329bd6b7.jpg
Типы проекций по характеру искажений
Для разных целей нужны карты с отсутствием тех или иных видов искажений. При помощи разных проекций можно сделать так чтобы на них отсутствовали погрешности либо углов, либо длин, либо площадей. Чем больше искажаются углы, тем меньше искажаются площади и наоборот. По характеру искажений все картографические проекции делят на:
- равноугольные (конформные);
- равновеликие (эквивалентные);
- произвольные:
- равнопромежуточные (эквидистантные).
Равноугольные картографические проекции
На картах, построенных по этому типу, нет искажений направлений и углов. Направления на местности совпадают с таковыми на карте, прямые линии на местности остаются прямыми на карте. Они используются для прокладки точных маршрутов и применяется на навигационных и топографических картах.
Зато на них сильно изменены площади объектов Земли и линейный масштаб карты зависит от положения на ней данной точки. Типичный пример равноугольной проекции – цилиндрическая проекция Герхарда Меркатора (Герарда Кремера), созданная ещё в 1569 г и используемая в морской навигации до сих пор. Примером использования Проекции Меркатора является равноугольная проекция Гаусса-Крюгера.
В этой проекции создаются отдельные океанологические, климатические и геофизические карты.
Проекция Меркатора.
Файл доступен по лицензии: Creative Commons Attribution-Share Alike 3.0 Unported
Равновеликие картографические проекции
Это проекции для построении карт, на которых нет искажения площадей (масштаб площадей имеет везде одну и ту же величину), зато сильно растёт погрешность форм и углов (материки и океаны в высоких широтах сплющиваются). Картами, построенными в равновеликих проекциях, удобно пользоваться для расчета площадей, например типов почв, посадок кукурузы, облесенности материков, загрязнения океана или радиоактивного загрязнения суши и др.
Их применяют для составления климатических, почвенных, геофизических, геологических, зоогеографических, геоботанических, экономических, исторических, этнографических, административных карт.
Пример равновеликой проекции. Автор: CC BY-SA 3.0
Произвольные картографические проекции
Углы и площади здесь искажаются, но значительно меньше, чем в предыдущих двух проекциях. Поэтому они наиболее используемы. Произвольные картографические проекции не относятся ни к равновеликим, ни к равноугольным.
Произвольная проекция Робинсона. Автор: CC BY-SA 3.0
Равнопромежуточные картографические проекции
Это тип произвольных картографических проекций. В них масштаб длин одного из главных направлений остаётся неизменным. Пример: прямая азимутальная проекция. Равнопромежуточные проекции используют для создания общегеографических, физических, тектонических, политических и др. видов карт.
Характер искажения всегда входит в общее название проекции (равновеликая азимутальная, равноугольная коническая, равновеликая цилиндрическая и т.д.).
Интересно,
что д ревнейшей картографической проекцией является гномическая проекция, применённая на картах звёздного неба Фалесом Милетским ещё в Древней Греции.
Равнопромежуточная коническая проекция. Автор: CC BY-SA 3.0
Классификация географических проекций по геометрической фигуре, являющейся вспомогательной поверхностью
На плоскость эллипсоид проектируют при помощи геометрических фигур, а поверхности, на которые он проектируется, могут быть секущими (разрезающей) фигуру или касательными (соприкасается, но не разрезает глобус) к ней. При этом на полученной карте касательные и секущие линии (стандартные) представлены неискажёнными.
Проекции также бывают по-разному ориентированы.
- Нормальными называют проекции, в которых оси вспомогательной поверхности совмещаются с осью земного эллипсоида или шара, а спроектированная поверхность размещается касательно к полюсу.
- Поперечными – ось располагают под прямым углом к оси Земли.
- Наклонными– под любым другим (непрямым) углом к оси Земли.
Нормальная, поперечная и косая ориентации площади проекции. Автор: Rylem
Поверхности, которые могут быть развёрнуты на плоскость или лист без растяжений, разрыва или усадки, называются разрабатываемыми поверхностями. Ими являются цилиндр, конус и плоскость. Поэтому по вспомогательной поверхности проекции делятся на:
- цилиндрические – вспомогательная поверхность – боковая цилиндра, касательная к эллипсоиду или секущая эллипсоида. Меридианы изображаются равностоящими параллельными прямыми, а параллели – прямыми, перпендикулярными меридианам. Пример – нормальная равноугольная цилиндрическая проекция Меркатора.
Цилиндрическая, коническая, азимутальная типы проекций.
Автор: Rylem — собственная работа, CC BY-SA 4.0
- псевдоцилиндрические – центральный меридиан на них представлен в виде отрезка прямой, другие меридианы длиннее центрального и изогнуты наружу. Параллели псевдоцилиндрических проекций – прямые линии.
Псевдоцилиндрическая проекция
- конические – боковая поверхность секущей или касательной конуса. Конической называется любая проекция, в которой меридианы представлены прямыми линиями, выходящими из одного центра и равноудаляющимися к периферии, а параллели – дуги, центрированные на вершине. При построении картографы чаще выбирают 2 основные параллели, которые могут быть секущими или касательными. Искажения масштаба и формы на них низкие. К северу и югу от стандартных параллелей расстояния растягиваются, а между стандартными параллелями расстояния сжимаются. Может использоваться и одна стандартная параллель, тогда с удалением от неё расстояния растягиваются. Применяется для территорий, вытянутых вдоль параллелей, например, все карты России построены в конических проекциях.
- псевдоконические – проекции, где центральный меридиан – прямая, остальные меридианы кривые линии, а параллели – прямые, промежутки между которыми уменьшаются к полюсам.
- азимутальные – вспомогательной поверхностью служит секущая или касательная плоскость. Параллели на них – полные окружности. Меридианы – их радиусы. По меридианам такая проекция является равнопромежуточной и сохраняет вдоль них главный масштаб. Именно разновидностью азимутальной проекции является первая известная на Земле гномическая проекция.
Типы азимутальных проекций
- поликонические – боковые вспомогательные поверхности нескольких касательных конусов, каждая из которых затем разворачивается на плоскость. Экватор и средний меридиан – перпендикулярные прямые, параллели – дуги, выпуклостью направленные к экватору, меридианы – кривые малой кривизны, направленные выпуклой стороной от центрального меридиана.
- условные – те, что ни входят ни в один из выше перечисленных классов. Параллели и меридианы на них являются кривыми очень разного вида.
Полное название проекций может быть следующим: косая азимутальная равновеликая, нормальная равноугольная цилиндрическая, произвольная поликоническая и т.д.
Источник
Ликбез по картографическим проекциям с картинками
Визуализация данных самого разного рода, имеющих некое географическое распределение, в последнее время получает все большее и большее распространение. Тут, на Хабре, статьи с картами встречаются чуть ли не каждую неделю. Карты в статьях очень разные, но роднит их одно: как правило, в них используются всего две картографические проекции, при том — не самые удачные из существующих. Мне бы хотелось дать несколько наглядных примеров проекций, которые выглядят более эстетично и лучше приспособлены для разных видов визуализации. В этой статье будут рассмотрены общемировые проекции и проекции большей части Земли, так как визуализация чего-либо на карте мира, пожалуй, является наиболее распространенной из подобных задач.
Легкое введение
Поскольку статья ориентирована на вопросы визуализации данных, я не буду касаться глубоко теории проекций (датумов, конформности, равноугольности и тому подобного), кроме общих принципов их построения. Также, я буду говорить тут о «проекциях», формально подразумевая «систему координат», coordinate reference system, потому что для карт таких масштабов не имеет смысла отдельно рассматривать проекцию и датум. Математики здесь тоже практически не будет, кроме простой геометрии. Желающие ознакомиться с математическими принципами, могут это сделать по статьям на Wolfram MathWorld. Так что изучающим программирование в области геоинформационных систем или их опытным пользователям, эта статья, возможно, будет не очень полезна.
Перед началом, объясню пару вещей. Все примеры будут даваться с использованием набора данных государственных границ с вот этого сайта и набора данных Blue Marble Next Generation с сайта NASA. Последний включает в себя синтезированные безоблачные снимки земной поверхности за каждый из двенадцати месяцев 2004-го года, что позволит внести некоторое разнообразие в иллюстрации.
Я очень люблю открытый софт, но использовать GDAL в данном случае мне показалось неэффективно — некоторых не очень ходовых, но полезных проекций в его реализации на данный момент либо нет, либо я плохо смотрел исходники, а потому иллюстрации я готовил в коммерческой программе GlobalMapper, которой пользуюсь уже много лет, и которая славится поддержкой внушительного списка систем координат.
Названия проекций и некоторые термины я буду давать и англоязычные, потому что если кому-то захочется поискать материалы по этой теме, русскоязычных источников в сети найдется несколько меньше (объем статей в Википедии на русском меньше в несколько раз). Для большинства проекций я постараюсь дать не только названия, но и коды EPSG и/или WKID, а также название проекции в библиотеке PROJ.4, широко используемой в открытом софте (например, в пакете R) для поддержки систем координат.
Некоторые проекции, возможно, окажутся кому-то знакомыми по картинке с xkcd, но все из них тут рассмотрены не будут.
Проблема
Начнем с того, что же это за самые распространенные проекции, и что с ними не так.
Первая проекция — так называемая «Географическая», она же – Geographic projection, Latitude/Longitude, Plate carrée EPSG:4326 WKID:54001 PROJ.4:longlat. Строго говоря, она даже не совсем является проекцией, потому что получается путем интерпретации полярных угловых координат, как линейных прямоугольных, без всяких вычислений. Эту проекцию используют, потому что она способна отобразить всю поверхность Земли целиком и потому, что она самая простая математически, а данные очень часто распространяются не спроецированными, то есть именно в географических координатах (градусах широты и долготы).
Что же получается? Получается прямоугольник, где точки полюсов обращены в линии (верхнюю и нижнюю границы). Чем дальше от экватора, тем сильнее любой объект на карте оказывается сплюснут по вертикали и растянут по горизонтали. Как я уже сказал, это худо-бедно годится для отображения глобальных наборов данных, но полярные территории (Канада, Норвегия, Швеция, север России, Финляндия, Гренландия, Антарктида, Исландия) оказываются искажены. Проекции, которые позволяют избежать этого, существуют, и о них пойдет речь дальше. Единственная причина использовать эту проекцию — ее предельная простота программной реализации — нужно просто отобразить систему координат от -180º до 180º по X и от -90º до 90º по Y на плоскость, считая угловые единицы линейными.
Другая весьма популярная проекция — «проекция Меркатора», Mercator projection PROJ.4:merc. Она также используется для визуализации данных, покрывающих весь мир, но ее популярность продиктована не только простотой — ее варианты являются стандартом де-факто для глобальных картографических сервисов, таких как Google Maps, Bing Maps, Here. С ней глубоко связаны картографические библиотеки OpenLayers, Leaflet, API упомянутых выше сервисов. В варианте Google и OpenStreetMap она носит название Web Mercator и имеет код EPSG/WKID:3857, иногда на нее также ссылаются, как на EPSG:900913. Принцип ее построения не сильно сложнее Географической – это проекция на цилиндр, чья ось совпадает с географической осью Земли, проецирование происходит линиями, выходящими из центра планеты, от чего ошибка растяжения приполярных областей по горизонтали оказывается скомпенсирована пропорциональным растяжением по вертикали. Проблема с этим только в том, что карта получится слишком большой по вертикали, если попытаться отобразить и север Гренландии. Потому обычно отбрасывают 16° полярных областей (в равной пропорции или больше — с юга).
На чей-то взгляд выглядит чуть лучше, чем Географическая, но одну проблему мы уже упомянули, а вторая — чем ближе объект к полюсам, тем он кажется больше, хотя его форма уже не так искажена. Потому, если предмет визуализации — плотность маркеров на единицу территории или расстояния, такой способ отображения будет вводить в заблуждение. При грамотном выборе способа визуализации, конечно, это можно скомпенсировать, а для каких-то случаев это вообще не проблема: например, если величина какого-то показателя в целой стране соотнесена с цветом этой страны на карте, эффект растяжения площадей не сказывается. Эта проекция сохраняет только форму объектов, потому очертания континентов и стран выглядят довольно узнаваемо. И, как я уже сказал, она — ваш первый и самый простой вариант при создании интерактивных веб-карт.
Варианты решения
Что же делать с глобальными данными, если нам по какой-то причине понадобилась проекция, лучше сохраняющая такие свойства объектов, как форма, площадь, расстояния и углы? Законы геометрии не дают нам сохранить все эти свойства сразу, развернув круглую поверхность Земли на плоскость. Однако, для визуализации данных более всего важна эстетика и восприятие, а не сохранение свойств, как для навигационных или измерительных задач. Потому становится возможным подобрать такую проекцию, искажения в которой были бы равномерно распределены по свойствам. И таких проекций существует довольно много. Существуют три самых известных, обладающих сходными свойствами: «Тройная проекция Винкеля» Winkel Tripel WKID:54042 PROJ.4:wintri, «проекция Робинсона» Robinson projection WKID:54030 PROJ.4:robin, «проекция Каврайского» (Kavrayskiy projection). Первая и последняя имеют визуально минимальные искажения, а неспециалисту, не видя градусной сетки, вообще весьма сложно различить их, потому я приведу иллюстрацию для Winkel Tripel, как той, которая лично мне нравится больше всего.
Вот так описание этой проекции выглядит в формате ESRI WKT:
PROJCS[«Robinson»,
GEOGCS[«GCS_WGS_1984»,
DATUM[«D_WGS84»,
SPHEROID[«WGS84»,6378137,298.257223563]
],
PRIMEM[«Greenwich»,0],
UNIT[«Degree»,0.017453292519943295]
],
PROJECTION[«Robinson»],
PARAMETER[«central_meridian»,0],
PARAMETER[«false_easting»,0],
PARAMETER[«false_northing»,0],
UNIT[«Meter»,1]
]
Как легко видеть, хотя искажение контуров и некоторое увеличение площади стран к полюсам здесь также наблюдаются, но это нельзя даже сравнивать с растяжением Географической проекции и пропорциональным увеличением проекции Меркатора.
Тут стоит сделать небольшое отступление и обратить внимание на то, что вид этой проекции по умолчанию страдает одним недостатком, который касается и других общемировых проекций. Дело в том, что если за центральный меридиан — линию, соединяющую северный и южный полюс через центр карты (longitude of origin) — принять нулевой меридиан, то карта будет разрезана по 180-му. Но при этом треть Чукотки окажется на левом краю карты, а две трети — на правом. Чтобы сделать карту красивее, разрез должен проходить где-то в районе 169-го западного меридиана восточнее острова Ратманова, для чего за центральный должен быть принят 11-й. Вот иллюстрация того, что получается:
А вот измененное для этого случая описание в ESRI WKT:
PROJCS[«Robinson»,
GEOGCS[«GCS_WGS_1984»,
DATUM[«D_WGS84»,
SPHEROID[«WGS84»,6378137,298.257223563]
],
PRIMEM[«Greenwich»,0],
UNIT[«Degree»,0.017453292519943295]
],
PROJECTION[«Robinson»],
PARAMETER[«central_meridian»,11],
PARAMETER[«false_easting»,0],
PARAMETER[«false_northing»,0],
UNIT[«Meter»,1]
]
В формате определения системы координат для PROJ.4 долгота центра проекции задается параметром +lon_0=.
11-й меридиан — «магическое» число: практически все мировые проекции, имеющие равномерный масштаб вдоль экватора, могут быть разрезаны по Берингову проливу, если за центральный принять именно его, а не нулевой.
Замечу, что задумываясь о выборе проекции, стоит принимать во внимание все существующие реальные требования к визуализации. Например, если данные касаются климата, то может иметь смысл либо нанести на карту линии широты, либо использовать проекцию, где они горизонтальны, а не загибаются к краям карты (то есть, отказаться от Тройной Винкеля в пользу, например, Робинсона). В данном случае, это позволит легче и точнее оценить относительную близость разных мест к полюсам и экватору. Еще один весомый плюс проекции Робинсона — то, что она поддерживается множеством софта, в том числе открытого, тогда как про некоторые другие этого сказать нельзя.
Иногда, когда требуется максимально сохранить какое-то свойство, например — соотношение площадей объектов (стран) — эстетическая сторона страдает. Но поскольку это все же может для чего-то понадобиться, я приведу один пример такой проекции — «проекцию Моллвейде», Mollweide projection WKID:54009 PROJ.4:moll.
Как видно, она довольно сильно напоминает проекцию Робинсона, но с той разницей, что полюса все же стянуты в точки, от чего форма приполярных областей выглядит сильно искаженной. Но пропорции площадей стран, как и требовалось, сохраняются куда лучше.
Самым молодым конкурентом этих проекций является проекция Natural Earth PROJ.4:natearth — она представляет из себя гибрид проекций Каврайского и Робинсона, а ее параметры были подобраны группой американских, швейцарских и словенских специалистов в 2007 году, тогда как возраст большинства картографических проекций — не менее полувека.
Для перепроецирования данных в нее существует некоторое количество инструментов, которые были написаны специально для этого, но ее поддержка еще далека от повсеместной.
Немного экзотики и специальных случаев
Конечно, все многообразие проекций на этом не заканчивается. Их изобретено немало. Некоторые просто выглядят странно (скажем, проекция Бонне изображает Землю в виде фигуры, напоминающей разрезанное яблоко или стилизованное сердце), некоторые — предназначены для особых ситуаций. Например, готов поспорить, что очень многие видели на картинках карту мира, которая похожа на корку мандарина, которую сняли и расплющили. Это, наверняка, была «Разрывная гомолосинусоидальная проекция Гуда» Interrupted Goode Homolosine projection WKID:54052.
Вид ее вполне достоин названия. Ее назначение — отображать размер объектов (и в некоторой степени — форму) близко к естественным пропорциям. Ее главная проблема, кроме названия и странного вида, состоит в том, что путем подбора центрального меридиана невозможно добиться того, чтобы ни один крупный кусок суши не был разрезан. Обязательно пострадает что-то из списка: Гренландия, Исландия, Чукотка, Аляска. Лично на мой взгляд, проще привести отдельно изображения стран, чем использовать такую карту, если вы не хотите стилизовать свою работу под середину XX века.
Существуют проекции, которые по своей природе никак не отнести к общемировым, но мне бы хотелось рассмотреть их здесь, потому что они способны показать земной шар, то есть как-бы вид планеты из космоса. Одна из них — «Проекция вертикальной ближней перспективы» Vertical Near-Side Perspective projection WKID:54049. Ее особое свойство — показывать земную поверхность в такой перспективе, как она выглядит с определенной высоты. Высота над эллипсоидом (идеализированной фигурой, моделирующей Землю) задается для этой проекции в явном виде.
На иллюстрации эта проекция имеет широту и долготу центра, равные широте и долготе Москвы, а высоту — 5000000 метров. Чем больше это расстояние, тем сильнее изображение Земли становится похоже на ее изображение в проекции, которую мы рассмотрим последней.
Проекция, которая показывает вид на Землю в параллельной перспективе, то есть как-бы с бесконечного расстояния, называется «Ортографическая проекция» Orthographic projection WKID:43041 PROJ.4:ortho. В каком-то смысле, она знакома всем, кто когда-либо пользовался Google Earth. Я говорю, что в каком-то смысле, потому что «направление взгляда» в этой проекции всегда перпендикулярно поверхности Земли, тогда как в Google Earth его можно наклонять как угодно.
Для нее, как и для предыдущей проекции, можно задать центральные широту и долготу, чтобы ориентировать Землю желаемым образом. Например, можно показать полушарие с центром в какой-то точке, о которой идет речь — скажем, иллюстрируя транспортные потоки континентального масштаба, исходящие от одного предприятия. Сделав две карты с противоположными значениями координат, можно получить карту всего мира (правда, на краях искажения будут очень велики). Генерация последовательности карт с плавным изменением центральной точки даст кадры для анимации вращающейся планеты без всякой трехмерной графики.
Если статья окажется интересной, постараюсь написать продолжение о проекциях, используемых для отображения отдельных стран или регионов, ориентированную, как и эта статья, на базовые свойства этих проекций для задачи визуализации данных, инфографики и тому подобного.
Источник