Приложение H: Создание собственной координатной системы
зависимости от проекции, которая используется, значение может быть от 20000 до 32760. Коды не могут использоваться более одного раза в файле mapinfow.prj. Для более подробной информации о файлах GeoTIFF и кодах PCS обратитесь к web - странице:
http://www.remotesensing.org/geotiff/geotiff.html
Вопрос: ”Как я могу конвертировать систему координат в метрах в координатную систему в футах?”
Ответ: Вам надо открыть файл Mapinfow.prj и изменить в нем эту же проекцию, переведя ее в другие единицы измерения, а также переустановить Восточное и Северное смещения.
В примере, описываемом ниже, первая строка это исходная проекция, выражаемая в метрах. Вторая строка - это измененная первая строка, где единицы измерения поменялись с ”7” (метры) на ”8” (футы), что означает использование единиц измерения в футах, Восточное Смещение (2000000) и Северное (500000) теперь надо разделить на
.3048, при переходе от метров к футам.
Пример:
”California, Zone I (1983)\p26941”, 3, 74, 7, -122, 39.3333333333, 40, 41.6666666667, 2000000, 500000
”California, Zone I FT (1983)\p26941”, 3, 74, 8 , -122, 39.3333333333, 40, 41.6666666667, 6561679.7, 164041.99
Вопрос: ”Я выбрал проекцию Широта/Долгота (NAD 83) и теперь, когда хочу поменять ее в диалоге "Выбор проекции" на Широта/Долгота (GRS 80), почему моя проекция не меняется?”
Ответ: Регион GRS 80 точно соответствует региону NAD 83. MapInfo использует числовые параметры, а не имена, для определения координатной системы, указанной в диалоге. Поскольку GRS 80 имеет теже числовые параметры, как и NAD 83, а GRS 80 находится в списке проекций выше, то MapInfo выбирает GRS 80 вместо NAD 83.
Для более подробной информации о проекциях
Рекомендуем Вам книгу Картография Л.С.Гараевской, Недра, 1971.
Картография. Выпуск 1. Москва, Прогресс, 1983 275 с.
Справочник геодезиста. Под ред. В.Д.Большакова и Г.П.Левчука. Москва, Недра, 1985, 454 с.
Адреса и телефоны:
American Congress on Surveying and Mapping 5410 Grosvenor Lane, Suite 100
Приложение H: Создание собственной координатной системы
Bethesda, MD 20814?2212 301-493-0200
Earth Science Information Center
U.S. Geological Survey
507 National Center
Reston, VA 22092
703-860-6045 or 1?800-USA-MAPS
Для подробной информации о геодезических регионах
http://www.utexas.edu/depts/grg/gcraft/notes/datum/datum.html
Для подробной информации о системах координат
http://www.utexas.edu/depts/grg/gcraft/notes/coordsys/coordsys.html
Для подробной информации о проекциях карт
http://www.utexas.edu/depts/grg/gcraft/notes/mapproj/mapproj.html
Приложение I: Создание собственного
референц-эллипсоида
Референц-эллипсоид – это математическая модель земной поверхности. Так как земная поверхность неоднородна, существует много моделей ее описания, приспособленных к разным областям мира. Эти модели позволяют достичь максимального приближения к реальности на некоторой выбранной местности.
Каждая координатная система использует некоторый локальный референц-эллипсоид, т.е. предназначенный для наилучшего моделирования части земного шара. Если две координатные системы используют разные референц-эллипсоиды, то MapInfo каждый раз проводит преобразование моделей при смене координатных систем. В MapInfo используется метод Бурса-Вольфа для преобразования референц-эллипсоидов, который гарантирует точность в 10 метров. (Если MapInfo проводит преобразование систем координат без смены референц-эллипсоида, то преобразования Бурса-Вольфа не проводится и ошибка не превышает 0,1 м.)
Большинство координатных систем используют один из референц-эллипсоидов MapInfo, перечисленных в Приложении F. Если Вы желаете использовать референцэллипсоид, не присутствующий в списке, то Вы можете создать собственный. Для этого нужно использовать следующие параметры:
Эллипсоид или сфероид. Это поверхность объемной фигуры, полученной вращением эллипса вокруг своей меньшей оси. Эллипсоид описывается двумя математическими параметрами: длиной, в метрах, его главной оси (обозначенной буквой a) и степенью сжатия (обозначенной буквой f). MapInfo поддерживает более 40 предопределенных эллипсоидов, представленных в таблице ниже.
Три параметра, описывающие сдвиг эллипсоида вдоль его осей. Обычно они обозначаются как dX, dY и dZ. Они также обозначаются как DX, DYи DZ или u, v и w.
Три параметра, описывающие углы поворота эллипсоида вокруг каждой из его осей в секундах. Они обозначаются либо как EX, EY и EZ, либо как eX, eY и eZ, либо как e, y и w.
Коэффициент масшабной коррекции, используемый для приближения эллипсоида к реальной земной поверхности. Он измеряется в частях на миллион (parts per million или ppm) и обозначается либо буквой m, либо буквой k.
Долгота начального меридиана, отсчитываемая к востоку от Гринвича. Этот начальный меридиан становится нулевым и его долгота принимается равной 0°. Большинство эллипсоидов используют Гринвичский меридиан в качестве начального, то есть начальный меридиан равен нулю. Другие эллипсоиды, однако, используют ненулевые начальные меридианы. Например, эллипсоид NTF использует как начальный Парижский меридиан, проходящий на 2.33722917 градуса к востоку от Гринвича. Все координатные системы, основанные на эллипсоиде NTF, отсчитывают долготы от Парижского меридиана, а не от Гринвича.
Собственный эллипсоид можно создать для любой координатной системы. В
Приложении F описано, как создавать координатную систему. Чтобы создать
Приложение I: Создание собственного референц-эллипсоида
координатную систему с собственным эллипсоидом, нужно использовать номер 9999, после которого указываются параметры в следующем порядке:
9999, НомерЭллипсоида, dX, dY, dZ, EX, EY, EZ, m, НачальныйМеридиан
Некоторые референц-эллипсоиды используют только эллипсоид и параметры сдвига (dX, dY, dZ), без параметров поворота, коэффициента коррекции и начального меридиана. В этом случае нужно использовать число 999 вместо 9999:
999, НомерЭллипсоида, dX, dY, dZ
Номер эллипсоида извлекается из следующей таблицы. Свой эллипсоид Вы задать не можете. Если Вы хотите, чтобы Ваш эллипсоид можно было использовать в следующих версиях MapInfo, известите об этом службу технической поддержки MapInfo.
Эллипсоид # |
Старое имя |
Новое имя/Определение NIMA |
11 |
Everest (1830) |
Everest (India 1830) |
|
|
|
17 |
Everest (Kertau) |
Everest (W. Malaysia and Singapore 1948) |
|
|
|
39 |
Everest (Timbalai) |
Everest (Brunei and East Malaysia (Sabah and |
|
|
Sarawak)) |
|
|
|
48 |
Everest(West Malaysia) |
Everest (West Malaysia 1969) |
|
|
|
40 |
Everest (Kalianpur) |
Everest (India 1956) |
|
|
|
Дополнения и новые эллипсоиды описаны в следующей таблице.
Номер |
Эллипсоид |
a |
1/f |
|
|
|
|
9 |
Airy 1930 |
6377563.396 |
299.3249646 |
13 |
Airy 1930 (modified for Ireland 1965) |
6377340.189 |
299.3249646 |
51 |
ATS7 77 |
6378135.0 |
298.257 |
2 |
Australian |
6378160.0 |
298.25 |
10 |
Bessel 1841 |
6377397.155 |
299.1528128 |
35 |
Bessel 1841 (modified for NGO 1948) |
6377492.0176 |
299.15281 |
14 |
Bessel 1841 (modified for Schwarzeck) |
6377483.865 |
299.1528128 |
36 |
Clarke 1858 |
6378293.639 |
294.26068 |
7 |
Clarke 1866 |
6378206.4 |
294.9786982 |
8 |
Clarke 1866 (modified for Michigan) |
6378450.047484481 |
294.9786982 |
6 |
Clarke 1880 |
6378249.145 |
293.465 |
|
|
|
|
15 |
Clarke 1880 (modified for Arc 1950) |
6378249.145326 |
293.4663076 |
|
|
|
|
30 |
Clarke 1880 (modified for IGN) |
6378249.2 |
293.4660213 |
|
|
|
|
37 |
Clarke 1880 (modified for Jamaica) |
6378249.136 |
293.46631 |
|
|
|
|
Приложение I: Создание собственного референц-эллипсоида
16 |
Clarke 1880 (modified for Merchich) |
6378249.2 |
293.46598 |
|
|
|
|
38 |
Clarke 1880 (modified for Palestine) |
6378300.79 |
293.46623 |
|
|
|
|
39 |
Everest (Brunei and East Malaysia (Sabah |
6377298.556 |
300.8017 |
|
and Sarawak)) |
|
|
|
|
|
|
11 |
Everest (India 1830) |
6377276.345 |
300.8017 |
|
|
|
|
40 |
Everest (India 1956) |
6377301.243 |
300.80174 |
|
|
|
|
50 |
Everest (Pakistan) |
6377309.613 |
300.8017 |
|
|
|
|
17 |
Everest (W. Malaysia and Singapore |
6377304.063 |
300.8017 |
|
1948) |
|
|
|
|
|
|
48 |
Everest (West Malaysia 1969) |
6377295.664 |
300.8017 |
|
|
|
|
18 |
Fischer 1960 |
6378166.0 |
298.3 |
|
|
|
|
19 |
Fischer 1960 (modified for South Asia) |
6378155.0 |
298.3 |
|
|
|
|
20 |
Fischer 1968 |
6378150.0 |
298.3 |
|
|
|
|
21 |
GRS 67 |
6378160.0 |
298.247167427 |
|
|
|
|
0 |
GRS 80 |
6378137.0 |
298.257222101 |
|
|
|
|
5 |
Hayford |
6378388.0 |
297.0 |
|
|
|
|
22 |
Helmert 1906 |
6378200.0 |
298.3 |
|
|
|
|
23 |
Hough |
6378270.0 |
297.0 |
|
|
|
|
31 |
IAG 75 |
6378140.0 |
298.257222 |
|
|
|
|
41 |
Indonesian |
6378160.0 |
298.247 |
|
|
|
|
4 |
International 1924 |
6378388.0 |
297.0 |
|
|
|
|
49 |
Irish (WOFO) |
6377542.178 |
299.325 |
|
|
|
|
3 |
Krassovsky |
6378245.0 |
298.3 |
|
|
|
|
32 |
MERIT 83 |
6378137.0 |
298.257 |
|
|
|
|
33 |
New International 1967 |
6378157.5 |
298.25 |
|
|
|
|
43 |
NWL 10D |
6378135.0 |
298.26 |
|
|
|
|
42 |
NWL 9D |
6378145.0 |
298.25 |
|
|
|
|
44 |
OSU86F |
6378136.2 |
298.25722 |
|
|
|
|
45 |
OSU91A |
6378136.3 |
298.25722 |
|
|
|
|
46 |
Plessis 1817 |
6376523.0 |
308.64 |
|
|
|
|
24 |
South American |
6378160.0 |
298.25 |
|
|
|
|
12 |
Sphere |
6370997.0 |
0.0 |
|
|
|
|
47 |
Struve 1860 |
6378297.0 |
294.73 |
|
|
|
|
34 |
Walbeck |
6376896.0 |
302.78 |
|
|
|
|
Приложение I: Создание собственного референц-эллипсоида
25 |
War Office |
6378300.583 |
296.0 |
|
|
|
|
26 |
WGS 60 |
6378165.0 |
298.3 |
|
|
|
|
27 |
WGS 66 |
6378145.0 |
298.25 |
|
|
|
|
1 |
WGS 72 |
6378135.0 |
298.26 |
|
|
|
|
28 |
WGS 84 |
6378137.0 |
298.257223563 |
|
|
|
|
Параметры сдвига и поворота, определяющие положение эллипсоида в пространстве, отсчитываются от референц-эллипсоида WGS 84. Важно следить за тем, чтобы эти параметры имели правильные знаки (положительные или отрицательные). Обычно в технической документации локальный референц-эллипсоид описывается как набор параметров преобразования в референц-эллипсоид WGS 84. (Фактически это означает, что эти параметры вычисляются как разница значений для локального РЭ и WGS 84). В технической документации также применяется описание обратного преобразования – из WGS 84 в локальный референц-эллипсоида. Поэтому важно не перепутать знаки (абсолютные значения при этом не меняются).
Также важно не перепутать порядок перечисления параметров. В некоторых документах параметры поворота перечисляются, начиная с EZ, т.е. в порядке EZ, EY, EX. Вы должны перечислять их в определении собственного референц-эллипсоида в обратном порядке. Легче всего запутаться с греческими буквами. Если в документе параметры перечислены как e, y, w, то Вы должны применять их в обратном порядке.
Вот пример описания локального референц-эллипсоида (назовем его LD-1) на языке, который принят в технической литературе:
Эллипсоид LD-1: Международный a 6378388.0 м
f1 / 297.0
Преобразование из WGS 84 в LD-1 DX 93.5 м
DY 103.5 ì
DZ 123.3 ì eX –0.25'' eY 0.11'' eZ 0.07ppm
m–2.1 ppm
Этот референц-эллипсоид использует международный эллипсоид, который в таблице имеет номер 4. Так как остальные параметры описывают преобразование из WGS 84, то все знаки надо поменять на обратные. Так как не указан нулевой меридиан, то значит используется Гринвичский. Вот как будет выглядеть соответствующее описание нашего референц-эллипсоида в MapInfo:
9999, 4, -93.5, -103.5, -123.3, 0.25, -0.11, -0.07, 2.1, 0
Приложение I: Создание собственного референц-эллипсоида
Эти числа можно поместить в строку описания координатной системы в файл
MAPINFOW.PRJ, как показано в следующих примерах: (Внимание: описание вводится
единой строкой).
”Долгота / Широта (LD-1)”, 1, 9999, 4, -93.5, -103.5, -123.3, 0.25, - 0.11, -0.07, 2.1, 0
”UTM Zone 30 (LD-1)”, 8, 9999, 4, -93.5, -103.5, -123.3, 0.25, -0.11, - 0.07, 2.1, 0, 7, -3, 0, 0.9996, 500000, 0
Рассмотрим другой пример определения локального референц-эллипсоида. Назовем его LD-2:
Эллипсоид LD-2: Красовский a 6378245.0 m
f 1 / 298.3 Преобразование из LD-2 в WGS 84
u +24 ì v –123 ì
w–94 ì
w+0.13''
y+0.25''
e –0.02''
m+1.1 · 10-6
Этот референц-эллипсоид использует эллипсоид Красовского, который в нашей таблице имеет номер 3. Обращать знаки в этом случае не нужно, так как параметры описывают преобразование в WGS 84. Однако, требуется перечислить в обратном порядке параметры поворота, так как параметр w идет первым. Вот как должен выглядеть результат:
9999, 3, 24, -123, -94, -0.02, 0.25, 0.13, 1.1, 0
В третьем примере, обозначенном как LD-3, описан только эллипсоид и сдвиги:
Эллипсоид LD-3: Кларк 1880 a 6378249.145 м
f1 / 293.465
Значения WGS 84 минус значения LD-3 dX –7 м
dY 36 ì dZ 225 ì
Этот референц-эллипсоид использует эллипсоид Кларка 1880, который в нашей таблице имеет номер 6. Обращать знаки и следить за порядком параметров поворота не нужно. Нужно только использовать число 999 вместо 9999. В результате можно использовать одну из двух последовательностей чисел:
Приложение I: Создание собственного референц-эллипсоида
999, 6, -7, 36, 225
9999, 6, -7, 36, 225, 0, 0, 0, 0, 0
Эти последовательности можно использовать в описании координатных систем:
”Долгота / Широта (LD-3)”, 1, 999, 6, -7, 36, 225
”UTM Zone 30 (LD-3)”, 8, 999, 6, -7, 36, 225, 7, -3, 0, 0.9996, 500000, 0
Преобразование референц-эллипсоида
Когда преобразуются координаты из одного референц-эллипсоида в другой, MapInfo использует методы Молоденского (3-параметра) и Бурса-Вольфа (7-параметров). Суть этих методов в преобразовании координат из одного референц-эллипсоида в другой.
После введения референц-эллипсоида NAD 83, NOAA разработала программу NADCON, по первым буквам North American Datum CONversion. Это весьма узкоспециализированная программа, преобразующая координаты только из NAD 27 в NAD 83 и обратно. Для этой специализированной задачи используется высокоточный алгоритм, чем общий метод Молоденского; NADCON дает точность около 0.1 метра, а метод Молоденского - точность только 10-30 метров. Большинство американских правительственных учереждений используют NADCON для конвертации между NAD 27 и NAD 83.
Начиная с MapInfo 4.1.2, алгоритм NADCON используется для преобразования координат между NAD 27 и NAD 83, если эти координаты внутри области, покрываемой NADCON (США, Пуэрто Рико и Виргинские острова). Если координаты вне этих областей, или используют референц-эллипсоид отличный от NAD 27 или NAD 83, то MapInfo использует методы преобразования Молоденского и Бурса-Вольфа.
Метод преобразования NADCON работает медленее, чем метод Молоденского. Если Вы хотите отключить конвертацию NADCON, добавьте в registry “NADCON”. Путь должен выглядеть так:
HKEY_LOCAL_MACHINE\Software\MapInfo\MapInfo\Common\NADCON
Если это значение равно 0, то будет использоваться метод Молоденского вместо метода NADCON.
При использовании 16-разрядной версии MapInfo, добавьте строчку в раздел [MAPINFO COMMON] в файле MAPINFO.INI, чтобы отключить преобразование NADCON:
NADCON=0
Приложение J: Формат обмена данными
MapInfo
В этом приложении описан формат обмена данными MapInfo. В частности, здесь описаны:
∙Заголовок файла формата MIF.
∙Секция данных в файлах формата MIF.
∙Коды типов линий, штриховки, символов и шрифтов в файлах формата MIF.
∙Файл формата MID.
∙Пример файла.
Этот универсальный формат позволяет сопоставлять различным графическим элементам настраиваемые данные. MIF-файл является текстовым (ASCII), поэтому файлы в таком формате можно редактировать, достаточно легко генерировать, причем они будут работать в любых системах, где работает MapInfo. Возможно, оптимальным способом изучения MapInfo Interchange Format (MIF) является просмотр примера, приведенного в конце этого приложения, по ходу чтения описания формата. Вы можете также создать свои примеры, экспортировав какие-нибудь файлы в формат MIF, а затем просматривать их в текстовом редакторе.
Данные MapInfo хранятся в двух файлах — графическая информация содержится в файлах с расширением .MIF, а текстовая в файлах с расширением .MID. Каждая строка текстовых данных отделяется от следующей строки либо символом возврата каретки (Carriage Return), либо символом возврата каретки вместе с новой строкой (Line Feed), либо только символом новой строки. MIF-файл состоит из двух частей: заголовка и секции данных. Заголовок файла содержит информацию о том, как создавать таблицы MapInfo, а секция данных содержит определения графических объектов.
Заголовок файла формата MIF
Заголовок файла формата MIF имеет следующий вид (в квадратные скобки взята информация, которая может отсутствовать).
VERSION n
Charset ”ИмяНабораСимволов” [ DELIMITER "<c>" ]
[ UNIQUE n,n.. ] [ INDEX n,n.. ] [ COORDSYS...] [ TRANSFORM...] COLUMNS n
Приложение J: Формат обмена данными MapInfo
<èìÿ> <òèï> <èìÿ> <òèï>
ž
ž
DATA
Version (Версия)
Для правильной обработки MIF-файла MapInfo необходимо знать, в какой версии программы он был создан. Номер версии задается предложении VERSION: VERSION 1, VERSION 2, VERSION 300 или VERSION 450. В третьей версии формата MIF-файла впервые стало возможным описание объектов, состоящих из нескольких секций (таких, как ломаные). В версии 450 повысился лимит максимально допустимого количества узлов для регионов и полилиний: теперь это 1048572 узлов для региона из одного полигона или для полилинии. Предел понижается до 7 узлов для каждых двух добавочных полигонов. Если объект с более, чем 32К узлов сохранен и таблица считывается в версии MapInfo, ранее, чем 4.5, то такой объект будет невидим. Объекты в таблице, не превышающие 32K будут видимы. Версия 4.5 содержит определение ширины линий в пунктах.
Charset (Набор символов)
В предложении Charset указывается, набор символов какой кодовой страницы используется в таблице. Например, “WindowsLatin1” означает, что используется стандартная для США и Западной Европы кодовая страница Windows (так называемая “ANSIкодировка” или “кодовая страница 1251”); “WindowsCyrillic” обозначает стандартный для России набор символов Windows (кодовая страница 1252); если указать слово “Neutral”, то MapInfo не будет выполнять преобразование символов в текущую кодовую страницу Windows. Если Вы не знаете точно, как обозначается используемый Вами набор символов, то выполните экспорт таблицы и изучите MIF– файл с помощью текстового редактора. Название используемого Вами набора символов будет указано в предложении Charset.
Delimiter (Разделитель)
Здесь в двойных кавычках указывается символ, использующийся в качестве разделителя, например:
DELIMITER ";".
Стандартное значение разделителя — символ табуляции; при использовании стандартного разделителя строка DELIMITER может быть опущена.
Unique (Уникальная колонка)
Здесь задается число. Это число указывает на столбец в базе данных: 3 указывает на третий столбец, 7 — на седьмой и так далее. Что происходит со столбцами, перечисленными в списке UNIQUE”? Представим, например, что имеется база данных о шоссейных дорогах. Каждому шоссе соответствует единственное название, но шоссе