Часто задаваемые вопросы по координатам, проекциям, системам координат: различия между версиями

Материал из GIS-Lab
Перейти к навигации Перейти к поиску
Строка 387: Строка 387:
|-
|-
|
|
'''Ответ: '''Перечисление параметров систем координат - неудобный процесс, который может занимать значительное время. Код EPSG - удобный общепринятый способ для обозначения системы координат со всеми параметрами. Гораздо удобнее просто назвать код системы координат, чем перечислять все ее параметры. EPSG - устаревшая аббревиатура обозначавшая ныне не существующую European Petroleum Survey Group, сейчас она переименована в OGP - Internaltional Association of Oil and Gas Producers. Тем не менее сама аббревиатура EPSG продолжает использоваться.  
'''Ответ: '''Перечисление параметров систем координат при обсуждении и переписке - неудобный процесс. Код EPSG - удобный общепринятый способ для обозначения системы координат сразу со всеми параметрами. Гораздо удобнее просто назвать код системы координат, чем перечислять все ее параметры. EPSG - устаревшая аббревиатура обозначавшая ныне не существующую European Petroleum Survey Group, переименованную сейчас в OGP - Internaltional Association of Oil and Gas Producers. Тем не менее сама аббревиатура EPSG продолжает использоваться.  


Например, код EPSG:26416 назначен проекции Гаусса-Крюгера зона 16, или, более точно: проекции Поперечного Меркатора с центральными меридианом 93 в.д., ложным смещением на восток 16500000 на эллипсоиде Крассовского и т.д.
Например, код EPSG:26416 назначен проекции Гаусса-Крюгера зона 16, или, более точно: проекции Поперечного Меркатора с центральными меридианом 93 в.д., ложным смещением на восток 16500000 на эллипсоиде Крассовского и т.д.
Полное описание этой системы координат (Proj.4) выглядит следующим образом:
   +proj=tmerc +lat_0=0 +lon_0=93 +k=1 +x_0=16500000 +y_0=0 +ellps=krass  
   +proj=tmerc +lat_0=0 +lon_0=93 +k=1 +x_0=16500000 +y_0=0 +ellps=krass  
   +towgs84=23.92,-141.27,-80.9,-0,0.35,0.82,-0.12 +units=m +no_defs
   +towgs84=23.92,-141.27,-80.9,-0,0.35,0.82,-0.12 +units=m +no_defs


Полную базу кодов можно найти [http://www.epsg.org здесь].
Полную базу кодов EPSG можно найти [http://www.epsg.org здесь].
|-
|-
| class="tr1" rowspan="2" align="center" valign="middle" |
| class="tr1" rowspan="2" align="center" valign="middle" |
|}
|}

Версия от 06:46, 6 декабря 2013

Эта страница опубликована в основном списке статей сайта
по адресу http://gis-lab.info/qa/proj-sk-faq.html


Дополнительные комментарии по часто задаваемым вопросам, некоторые из тем более подробно разъясняются в других статьях и разделах, которые обозначены дополнительными ссылками.



Вопросы

  1. Что такое координаты? >>>
  2. Какие бывают системы координат? >>>
  3. Если мои данные хранятся в WGS 1984 или Pulkovo 1942 значит ли это что они спроецированы? >>>
  4. Как определить спроецированы мои данные или нет? Как определить тип системы координат? >>>
  5. Мои данные в спроецированной системе координат, как определить в какой именно? >>>
  6. Что такое файл prj и для чего он нужен? >>>
  7. Как лучше хранить данные – спроецированными или неспроецированными? >>>
  8. Я хочу перевести данные из одной проекции в другую, как это сделать? >>>
  9. Я хочу привязать данные, как это сделать? >>>
  10. Может ли растровый слой хранится в географической системе координат? >>>
  11. Я привязал топографические карты (Генштаб СССР) и хочу проверить качество привязки. Как мне это сделать? >>>
  12. Как создать дополнительную систему координат в ArcGIS? >>>
  13. Чем отличаются следующие проекции в ArcGIS: Pulkovo_1942_GK_Zone_9N и Pulkovo_1942_GK_Zone_9? >>>
  14. Выбранная проекция искажает один из основных параметров (длину, угол, площадь) как провести измерения более точно? >>>
  15. Что такое система координат GCS_Assumed_Geographic_1/Datum: D_North_American_1927? >>>
  16. Я загрузил(а) данные, они показываются неправильно, все сплющенное, как исправить? >>>
  17. Я создал(а) с помощью программы X (получил/скачал) растр, но он не ложится туда куда мне нужно, как подготовить их к работе? >>>
  18. Мои геоданные не имеют явным образом обозначенной системы координат, что мне делать? >>>
  19. Как определить, какой зоне UTM, GK относятся мои данные? >>>
  20. В чем разница между Географической проекций и проекцией Меркатора >>>
  21. Как назначить слою систему координат/создать файл prj? >>>
  22. Я задаю слою разные системы координат, но он не меняется, почему? >>>
  23. Я посмотрел координаты слоя в проекте и использую их для вырезки нужного мне фрагмента, но ничего не происходит, почему? >>>
  24. Я измеряю размеры (длину, площадь) и они сильно отличаются от реальных, почему? >>>
  25. Я провожу расчеты связанные с расстояниями, площадями и результирующие значения очень сильно отличаются от ожидаемых, почему? >>>
  26. При расчете площади в проекции Меркатора результат расчета отличается от реальной площади в несколько раз, почему такая ошибка? >>>
  27. Как правильно сообщать результаты расчетов длин и площадей коллегам, чтобы они смогли при желании их воспроизвести? >>>
  28. Что такое код EPSG? Нужно ли мне его знать? >>>

Ответы

1


Вопрос: Что такое координаты?
Ответ: Координатами называют угловые или линейные величины, определяющие положение точки на какой-либо поверхности или в пространстве (Военная топография, 1958, с.125).

Географическими координатами называют угловые величины (долгота и широта), определяющие положение точек на эллипсоиде.

Спроецированными (плоскими, прямоугольными, декартовыми) координатами называют линейные величины (абсцисса и ордината, X и Y), определяющие относительное положение точки на плоскости.

2


Вопрос: Какие бывают системы координат?
Ответ: Существует много различных систем координат. Две из них являются наиболее распространенными: это географическая система координат (далее ГСК) и спроецированная (плоская, прямоугольная, декартова) система координат (ПСК).

В ГСК координаты хранятся в десятичных градусах и их иногда называют неспроецированными.

В ПСК единицами измерения могут быть метры, футы, километры и т.д. Данные находящиеся в этой системе координат часто называют спроецированными. В основе любой спроецированной системы координат лежит географическая, трехмерная система координат (ГСК), преобразованная в двухмерную с помощью проекции.

Таким образом, для описания ПСК, необходимо знать ГСК, на которой она основана и метод, которым координаты из ГСК переведены в ПСК (проекцию). Если проекция не известна, такие данные нельзя перевести в другую СК, иногда такую СК называют "локальной", "пиксельной" или "файловой" (например ПСК только что отсканированного изображения).

В настоящее время, наиболее распространенная географическая система координат в мире – система координат WGS 1984. В России чаще всего используется географическая система координат Pulkovo 1942. Подробнее

3

Вопрос: Если мои данные хранятся в WGS 1984 или Pulkovo 1942 значит ли это что они спроецированы?

Ответ: Не обязательно, это значит только то, что ваши данные основаны на данной географической системе координат. Данные могут быть как неспроецированными, так и спроецированными. Другими словами, знать то, в WGS84 ваши данные или Pulkovo 1942 - мало.

Примеры:

  • Данные в ГСК WGS 1984 - данные не спроецированы. Единицы измерения - градусы.
  • Данные в ПСК UTM, зона 34 (WGS 1984) - данные спроецированы и имеют проекцию UTM, зона 34, базирующуюся на ГСК WGS 1984, эллипсоид WGS 1984. Единицы измерения - метры.
  • Данные в ПСК ГК, зона 4 (Pulkovo 1942) - данные спроецированы и имеют проекцию ГК, зона 4, базирующуюся на ГСК Pulkovo 1942, эллипсоид Красовского. Единицы измерения - метры.

4

Вопрос: Как определить спроецированы мои данные или нет? Как определить тип системы координат?

Ответ: Самый простой способ определить спроецированы ли ваши данные, загрузить их в ПО по ГИС (QGIS, ArcView, ArcGis и т.д.), и обратить внимание на цифры, меняющиеся при перемещении курсора в области ваших данных, если Вы видите цифры:

  • сотни тысяч, миллионы - то ваши данные, скорее всего, находятся в спроецированной системе координат;
  • единицы, десятки, сотни (до +/-180 по горизонтали и до +/-90 по вертикали, но не больше) и координаты верхнего левого угла не равны 0,0 - то вероятнее всего ваши данные находятся в географической системе координат;
  • единицы, десятки, сотни, тысячи и координаты верхнего левого угла равны 0,0 или близки к ним - данные не привязаны и находятся в т.н. локальной (файловой) системе координат.

Этот способ не является универсальным, но может помочь в большинстве случаев. Этот способ позволяет определить только тип СК, он не позволяет уверенно определить проекцию и точную СК данных.

5

Вопрос: Мои данные в спроецированной системе координат, как определить в какой именно?

Ответ: Часто эта проблема является «тяжелым» случаем.

  1. Для векторных данных определить присутствует ли файл prj (см. следующий вопрос), для растровых данных посмотреть свойства растра с помощью gdalinfo (так же можно загрузить в Erdas Imagine или ArcCatalog и посмотреть свойства растра);
  2. Поискать описание системы координат и проекции в дополнительной информации (других файлах описания) или запросить в первоисточнике информации;
  3. Попробовать назначить данным основные распространенные спроецированные системы координат часто используемые для данной территории и наложить векторные данные в известной системе координат для проверки совпадают ли они.

6

Вопрос: Что такое файл prj и для чего он нужен?

Ответ: Если данные правильно документированы, то векторные слои должны содержать дополнительный файл с расширением prj (в покрытии Arcinfo этот файл находится внутри папки покрытия). Этот файл содержит исчерпывающее описание системы координат, в которой находятся ваши данные. Наличие этого файла позволяет ПО ГИС узнавать систему координат данных и при необходимости изменять ее (временно или постоянно).

Разные ГИС по разному работают с файлами описания проекций. В Arcview GIS он не создается при создании слоев, а только при использовании специального расширения (Projection Utility Wizard). В ArcGIS систему координат можно задать с помощью ArcCatalog.

Если данный файл отсутствует, то некоторое ПО, например ArcGIS присвоит данным некоторую СК по умолчанию. Отсутствие файла также сделает невозможным операции перевода данных из одной СК в другую в ArcGIS. В Arcview GIS это не проблема, если ваши данные находятся в географической СК и файл описания проекции отсутствует, программа будет изменять данные в зависимости от проекции вида.

Рекомендуется всегда иметь правильный prj-файл для каждого набора векторных данных, если их нет - создайте их.

7

Вопрос: Как лучше хранить данные – спроецированными или неспроецированными?

Ответ: Выбор типа системы координат зависит от нескольких переменных:

  1. Типа данных: вектор, растр или табличные данные;
  2. Идет ли речь о хранении, использовании или расчетах, использующих эти данные.

Обычно растровые данные хранятся спроецированными, векторные - неспроецированными в географической системе координат (долгота/широта). Однако для повышения скорости работы и простоты бывает полезным спроецировать все данные в определенную, одинаковую для всех данных, проекцию, базирующуюся на одной и той же СК. У этого подхода один большой недостаток, он плохо работает для данных покрывающих большую территорию, где часто приходится переводить их из одной проекции в другую.

В любом случае, какая бы не была система координат, она должна быть явным образом описана для всех слоёв (независимо от их типа) со всеми необходимыми подробностями.

8

Вопрос: Я хочу перевести данные из одной проекции в другую, как это сделать?

Ответ: Зависит от типа данных и доступного ПО:

Для векторных данных:

  • ArcView 3.x - использовать расширение Projection Utility Wizard (подробнее)
  • ArcGis 8.x, 9.x – инструмент Project (Data Management Tools\Projections and Transformations\Feature)
  • ArcInfo – команда project cover

Для растровых данных:

  • ArcView (только растры в формате GRID, при наличие Image Analysis, Spatial Analyst) – скрипт ProjectGridArcView.
  • ArcGis 8.x, 9.x /Arctoolbox – инструмент Project Raster (Data Management Tools\Projections and Transformations\Raster)
  • ArcGis 8.x, 9.x – загрузка в проект, изменение СК фрейма данных на необходимую, экспорт растра в новый с присвоением новой нужной СК
  • ArcInfo (только растры в формате GRID) – команда project grid
  • Erdas Imagine - DataPrep\Reproject images (подробнее)

9

'Вопрос: 'Я хочу привязать данные, как это сделать?

Ответ: Процесс имеет смысл, если вам неизвестна система координат исходных данных или она локальная (файловая), например только что отсканированный растр. Технологический процесс зависит от типа данных и доступного ПО:

Для векторных данных:

  • ArcView 3.x - использовать расширение ShapeWarp(подробнее)
  • Arcinfo Workstation - использовать команды ADJUST и TRANSFORM
  • ArcGis 8.x, 9.x – Spatial Adjustment Tool (подробнее)

Для растровых данных:

  • ArcView 3.x - при наличие Image Analysis, использовать инструмент Align, при наличие Spatial Analyst – использовать расширение ImageWarp (подробнее)
  • ArcGis 8.x, 9.x – Georeferencing Tool (подробнее)
  • ERDAS Imagine – geometric correction tool (Viewer\Raster\Geometric correction) (подробнее)

10

Вопрос: Может ли растровый слой хранится в географической системе координат?

Ответ: Да, технически такое возможно, однако проводить измерения по таким данным будет некорректно, так как площадь элемента растра (пиксел) будет зависеть от широты. Градусы - не линейные единицы изменения, а угловые.

11

Вопрос: Я привязал топографические карты (Генштаб СССР) и хочу проверить качество привязки. Как мне это сделать?

Ответ: Скачать километровые сетки для каждой конкретной зоны проекции Гаусса-Крюгера и проверить насколько точно она совпадает с метровой сеткой, нанесенной на топографическую карту (подробнее). Для проверки точности можно также использовать разграфки топокарт.

12

Вопрос: Как создать дополнительную систему координат в ArcGIS?

Ответ: Создать ее через свойства фрейма данных, закладка Coordinate Systems и переместить в папку c:\Program Files\ArcGIS\Coordinate Systems\MyProjections. Подробнее

13

Вопрос: Чем отличаются следующие проекции в ArcGIS: Pulkovo_1942_GK_Zone_9N и Pulkovo_1942_GK_Zone_9?

Ответ: Ложным смещением на восток (false easting), у 9N - 500000, а у 9 - 9500000. В практическом смысле это мало на что влияет, только то, как будут отображаться координаты по оси X.

Использование в качестве ложного смещения на восток значения [номер зоны]*106+500000 удобно например тем, что по координатам можно сразу определить с большой степенью вероятности проекцию и ее параметры.

14

Вопрос: Выбранная проекция искажает один из основных параметров (длину, угол, площадь), как провести измерения более точно?

Ответ: Один из способов получать результат в нужных единицах (например метрах), не используя при этом проекцию - проводить измерения прямо на сфере, это гарантирует заведомо более точный, а главное не зависящий от выбранной проекции результат. Например, измерения длин на сфере осуществляются с точностью около 0.5% (подробнее)

15

Вопрос: Что такое система координат GCS_Assumed_Geographic_1/Datum: D_North_American_1927 ?

Ответ: СК GCS_Assumed_Geographic_1/Datum: D_North_American_1927 появляется при открытии в ArcGIS любых данных в формате shape не имеющих явного описания СК (в виде файла *.prj). Данная СК указывает, что ArcGIS подразумевает, что данные находятся в этой географической СК, однако на самом деле они могут находится в любой другой географической СК.

Формат shape не требует наличия явного описания системы координат и/или проекции. Поэтому, получив из некоторого источника только файлы *.shp, *.dbf, *.shx, невозможно абсолютно точно сказать, в какой СК они находятся. Если СК данных известна, то это проблема решается созданием дополнительного файла *.prj, играющего роль такого описания проекции. Наличие этого файла с соответствующим описанием позволяет программному обеспечению (такому как ArcGIS или Arcview Projection Utility) "распознавать" СК и осуществлять перевод данных из одной СК в другую.

16

Я загрузил(а) данные, они показываются неправильно, все сплющенное, как исправить?

Ответ: Если данные привязаны, то показываются они, скорее всего, так, как и должны. Визуализация географических данных зависит от той системы координат, в которой они находятся исходно и той системы координат, в которой они отображаются. Эти две системы могут быть различными. Визуальная "привычность" картографической информации достигается подбором нужной системы (как правило, наиболее часто используемой для данной территории). Сами данные при этом остаются неизменными. См. также #22

Вот как выглядит п-ов Камчатка в географической системе координат:

Proj-sk-faq-01.gif

А вот те же самые данные, но спроецированные (проекция UTM, зона 54):

Proj-sk-faq-02.gif

Таким образом, для визуального изменения данных, измените проекцию в которой они отображаются.

17

Я создал с помощью программы X (получил/скачал) геоданные, как подготовить их к работе?

Ответ: Сначала необходимо выяснить, в какой системе координат находятся ваши данные. Эта информация может быть получена из описания программ, с помощью которых вы эти данные создали или в метаданных (дополнительной документации). Возможно удасться получить эту информацию от тех, кто поставил данные или, наконец, угадать систему координат.

После этого, если этого еще не сделано, необходимо явным образом "прописать" эту систему координат этому набору данных, чтобы она стала доступна не только Вам, но и программам, умеющим считыватья информацию о системе координат из заголовков или специальных дополнительных файлов.

18

Мои геоданные не имеют явным образом обозначенной системы координат, что мне делать?

Ответ: Если вы работаете в ГИС, которая не считывает информацию о системе координат из заголовков или специальных файлов, например Arcview, то можно ничего не делать, хотя для лучшей организации данных рекомендуется все-таки правильно подготовить данные. Если ваша ГИС требует явного описания системы координат данных см. 18.

19

Как определить какой зоне UTM или ГК относятся мои данные?

Ответ: Скачать разграфку зон UTM и наложить ее на свои данные, которые при этом должны находиться в географической системе координат. Посмотреть атрибутику зоны в которую попадают данные.

Можно также воспользоваться простой схемой (цифра указывает зону Гаусса-Крюгера, прибавьте к ней 30 и получите зону UTM).

Gk-grid01.gif

20

В чем разница между Географической проекций и проекцией Меркатора

Ответ: Во всем, кроме того, что они обе цилиндрические.

Географическая проекция (Plate Carrée, равнопромежуточная):

Proj-sk-faq-03.gif

Меркатор:

Proj-sk-faq-04.gif

21

Как назначить слою систему координат/создать файл prj?

Ответ: Назначить слою систему координат (создать файл prj) как правило можно в любой настольной ГИС или вручную. Для этого абсолютно необходимо точно знать, в какой СК и проекции находятся данные. Если вы этого не знаете, это необходимо сначала выяснить.

QGIS

Инструменты\Управление данными\Задать текущую проекцию
(Vector\Data Management Tools\Define current projection)

ArcCatalog

Назначение системы координат слою в ArcGIS

22

Я задаю слою разные системы координат, но его представление не меняется, почему?

Ответ: потому что у конкретного слоя физически только одна система координат и менять ее на другую (без изменения самих данных слоя) бесполезно. Следует различать родную систему координат данных (слоя) и систему координат представления слоя. Первая - та, в которой слой был создан и она описывается в его файле prj. Вторая - та, в которой слой показывается и она может быть изменена так, как вам необходимо.

Обычно, в ПО ГИС существует две возможности изменить СК представления слоя:

1. Инструмент, позволяющий перманентно пересчитать все данные из исходной СК в конечную. После этого нужно загрузить новый пересчитанный слой.

2. Свойства вида данных, где указывается новая СК для всех данных (может быть несколько слоёв). Этот способ позволяет не создавать копию данных в другой СК. Физически данные останутся в исходной СК, но показываться будут в нужной вам.

23

Я посмотрел координаты слоя в проекте и использую их для вырезки нужного мне фрагмента, но ничего не происходит, почему?

Ответ: потому что у слоя физически одна система координат, а в проекте он может быть показан в другой. Соответственно, вводить эти (производные) координаты бессмысленно. Используйте исходные координаты слоя, а не трансформированные (проецированные). Более подробно см. 22

24

Я измеряю размеры (длину, площадь) и они сильно отличаются от реальных, почему?

Ответ: как правило, измерения непосредственно по карте производятся в системе координат этой карты. Таким образом, если для карты выбрана проекция искажающая длины, площади или углы, то результаты измерений могут сильно отличаться от реальных.

Например, сегодня часто приходится сталкиваться с данными в проекции Меркатора (EPSG: 3857), эта проекция была создана для мореплавателей, для которых главное - постоянность углов, расстояния в этой проекции сильно искажаются. Если измерить длину любого объекта, она может в разы отличаться от реальной. Это не означает, что объект (данные) неправильный, просто его представление не подходит для задач измерений расстояний.

Проводите измерения в подходящих проекциях или на эллипсоиде, см. 14

25

Я провожу расчеты связанные с расстояниями, площадями и результирующие значения очень сильно отличаются от ожидаемых, почему?

Ответ: Так же как и измерения по карте в ответе выше, расчеты по объектам производятся в системе координат этих объектов. Таким образом, если единицы измерения данных - градусы, то построив буферную зону размером 500, вы получите зону радиусом вовсе не 500 метров, а 500 градусов.

Эта проблема может решаться двумя способами:

  1. Предварительным переводом данных в подходящую систему координат. Например, если нужно построить буферные зоны (т.е. главное - расстояния) переводите данные в одну из равнопромежуточных (equidistant) СК или, если масштаб достаточно детален, то можно использовать и соответствующую зону UTM, GK.
  2. Использованием ПО, которое умеет производить расчеты в СК вида данных, а не нативной СК данных. Таким образом, вам не придется делать лишнюю копию данных по совету выше.
  3. Использованием ПО, которое умеет производить расчеты на сфероиде (эллипсоиде) вообще без использования проекционных преобразований.

26

При расчете площади в проекции Меркатора результат расчета отличается от реальной площади в несколько раз, почему такая ошибка??

Ответ: Проекция Меркатора была изобретена для мореплавателей, которым было необходимо держать постоянное направление. Эта проекция плохо подходит для измерения площадей, поскольку в этой проекции масштаб равен единице только на экваторе. Измеренные в этой проекции расстояния и площади могут отличаться от ожидаемых в 2 и более раз.

В 25 можно найти ответ на вопрос "как же получить правильные площади".


27

Как правильно сообщать результаты расчетов длин и площадей коллегам, чтобы они смогли при желании их воспроизвести?

Ответ: При передаче таких данных коллегам, учтите, что проекции и эллипсоиды на которых вы проводили вычисления бывают разные, поэтому нужно обязательно сопровождать результаты расчетов пояснением в какой проекции или на каком эллипсоиде они производились. Иначе велик шанс, что ваш коллега при повторном расчете по вашим данным получит отличающиеся от ваших значения.

См. также: 25

28

Что такое код EPSG? Нужно ли мне его знать.

Ответ: Перечисление параметров систем координат при обсуждении и переписке - неудобный процесс. Код EPSG - удобный общепринятый способ для обозначения системы координат сразу со всеми параметрами. Гораздо удобнее просто назвать код системы координат, чем перечислять все ее параметры. EPSG - устаревшая аббревиатура обозначавшая ныне не существующую European Petroleum Survey Group, переименованную сейчас в OGP - Internaltional Association of Oil and Gas Producers. Тем не менее сама аббревиатура EPSG продолжает использоваться.

Например, код EPSG:26416 назначен проекции Гаусса-Крюгера зона 16, или, более точно: проекции Поперечного Меркатора с центральными меридианом 93 в.д., ложным смещением на восток 16500000 на эллипсоиде Крассовского и т.д.

Полное описание этой системы координат (Proj.4) выглядит следующим образом:

 +proj=tmerc +lat_0=0 +lon_0=93 +k=1 +x_0=16500000 +y_0=0 +ellps=krass 
 +towgs84=23.92,-141.27,-80.9,-0,0.35,0.82,-0.12 +units=m +no_defs

Полную базу кодов EPSG можно найти здесь.