Трансформация описания систем координат из формата MapInfo в WKT и PROJ.4

Эта страница является черновиком статьи.

Примечание: Если хотите прочитать эту статью вы должны знать, что такое система координат, датум, проекция, эллипсоид.

Введение

В настоящее время самыми распространёнными форматами записи описания систем координат являются WKT и proj

WKT - Well-known text текстовый формат представления геометрии и описания СК, также существует WKB (well-known binary) – бинарный аналог этого описания. Данный формат впервые был официально представлен в документе OGC 99-049 (OGC это Open GIS Consortium до 2004г и Open Geospatial Consortium после), в настоящее время актуальная редакция - OGC 06-103r3 и проект новой редакции OGC 12-063r2. Данный формат является частью технологии Simple Feature Access – целью данной технологии является предоставления общего интерфейса для обмена пространёнными данными между различными программами и сервисами.

Также компания ESRI использует несколько иную реализацию формата WKT, для описания СК – эта реализация обычно называется ESRI WKT и используется, в частности, в ArcGIS.

Вторым, по распространённости, форматом является, так называемый proj-формат – это формат описания систем координат программы proj. В настоящее время библиотека proj часто используется в различных программах и сервисах, вследствие этого этот формат получил широкое распространение de facto. Структура этого формата описана в документации к proj.

Описание СК в программе MapInfo находится в файле mapinfo.prj, структура этого формата изложина в документации MapInfo.

Примеры описания СК (WGS84.LL):

OGC WKT

     GEOGCS["WGS 84",
   DATUM["WGS_1984",
       SPHEROID["WGS 84",6378137,298.257223563,
           AUTHORITY["EPSG","7030"]],
       AUTHORITY["EPSG","6326"]],
   PRIMEM["Greenwich",0,
       AUTHORITY["EPSG","8901"]],
   UNIT["degree",0.01745329251994328,
       AUTHORITY["EPSG","9122"]],
   AUTHORITY["EPSG","4326"]]

ESRI WKT

     GEOGCS["GCS_WGS_1984",
   DATUM["D_WGS_1984",
       SPHEROID["WGS_1984",6378137,298.257223563]],
   PRIMEM["Greenwich",0],
   UNIT["Degree",0.017453292519943295]]

proj

     +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs

MapInfo.prj

     "Долгота / Широта (WGS 84)", 1, 104