Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

MI6UGRU

.pdf
Скачиваний:
16
Добавлен:
13.02.2015
Размер:
14.71 Mб
Скачать

Приложение 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

707

Приложение 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

708

Приложение 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

 

 

 

 

710

Приложение 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

 

 

 

 

711

Приложение 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

712

Приложение 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. В результате можно использовать одну из двух последовательностей чисел:

713

Приложение 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

714

Приложение 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”? Представим, например, что имеется база данных о шоссейных дорогах. Каждому шоссе соответствует единственное название, но шоссе

716

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]