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

Введение в геоинформаци онные системы

.doc
Скачиваний:
26
Добавлен:
18.03.2015
Размер:
2.6 Mб
Скачать

Запросы и выборки

<<< предыдущая глава | оглавление | следующая глава >>>

Теоретическая часть

Язык запросов - мощный инструмент по работе с условиями и критериями. Использование запросов позволяет отбирать записи в таблице и объекты на карте по определенным условиям. Это может быть одно условие/критерий или целое выражение, содержащее большое число условий.

Запрос - выражение на специальном языке запросов, включающее логические операторы и название слоев\полей из которых осуществляется выборка.

Выборка - процесс и результат выполнения запроса. Набор объектов\записей в таблице ему соотвествующий.

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

Атрибутивные запросы

Работа осуществляется с атрибутивной таблицей. Выбор записей в таблице позволяет работать с заданным подмножеством данных. Если таблица принадлежит теме, то при выборе записей в таблице окажутся выделенными и объекты в слое, которые соответствуют выбранным записям.

Атрибутивный запрос представляет собой некоторое логическое условие-команду, которое формируется также с применением команд стандартной логики (И, НЕ, ИЛИ).

Пример простейшего запроса:

На человеческом языке: Выбрать все квартала с возрастом преобладающей породы более 80 лет

На языке запросов: [age] > 80

Как видно, в атрибутивной таблице слоя «кварталы» должно содержаться поле age, в котором должны быть проставлены значения возраста для каждого квартала, если значение определенной строки в поле age будет отсутствовать, то запрос проигнорирует эту строку.

Логических условий в одном запросе может быть несколько. Например:

На человеческом языке: Выбрать все квартала с возрастом преобладающей породы между 80-ю и 120 годами

На языке запросов: ([age] > 80) AND ([age] < 120)

Как видно, в запросе было использовано логическое условие И (AND), то есть выбираемое значение должно удовлетворять сразу двум условиям, быть больше 80 лет и быть меньше 120 лет.

В качестве данных в запросе может быть использованы данные из нескольких полей. Например:

На человеческом языке: Выбрать все квартала с возрастом преобладающей породы более 80 лет и бонитетом равным единице

На языке запросов: ([age] > 80) AND ([bonitet] = 1)

Как видно, в атрибутивной таблице слоя «кварталы» также должно содержаться и поле bonitet, в котором должны быть проставлены значения бонитета для каждого квартала.

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

Пространственные запросы

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

Для выбора в исходном слое могут быть использованы следующие запросы:

Выбрать объекты

  • пересекающиеся с границей объекта другого слоя (are crossed by the outline of)

  • пересекающиеся с объектом другого слоя (intersect), отличие от предыдущего метода заключается в том , что при данном типе выборки будут также выбраны объекты касающиеся границ объекта другого слоя;

  • находящиеся на расстоянии Х от объекта другого слоя (Are within a distance of);

  • чей центр находится внутри объекта другого слоя - только для полигональных слоев (Have their center in);

  • находящиеся полностью внутри объекта другой темы (Are completely within), объект не будет выбран если его границы касаются границ объекта другой темы;

  • находящиеся полностью внутри объекта другой темы (Are contained by), объект будет выбран даже если его границы касаются границ объекта другой темы;

  • полностью включающие объект ( ы ) другой темы (Completely contain), объект не будет выбран если его границы касаются границ объекта другой темы;

  • полностью включающие объект ( ы ) другой темы (Contain), объект будет выбран даже если его границы касаются границ объекта другой темы;

  • имеющие общий фрагмент линии с объектом другой темы (Share a line segment with);

  • имеющие общий фрагмент линии с объектом другой темы , но не пересекающие его (Touch the boundary of)

  • равные объекту другой темы (Are identical to)

Выборка может осуществляться как для всех объектов одной и другой темы, так и для выделенных объектов в одной и другой теме. Тем, в которых будут выделяться объекты может быть несколько. Некоторые правила не работают для объектов определенного типа, например правило «выбрать объекты имеющие общий фрагмент линии с объектом другой темы» не будет работать с точечным слоем.

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

Практическая часть

Скачать учебные материалы для этой главы.

В практической части этой главы Вы научитесь:

  • Работать с атрибутивными запросами;

  • Работать с пространственными запросами;

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

Упражнение 1. Атрибутивные запросы.

  1. Откройте проект chapt16.mxd. В нем содержится тема представляющая собой фрагмент лесоустройства на территорию Сукпайского лесничества - les.shp. Элементарной единицей этого слоя является лесной выдел.

  2. Перед тем как производить атрибутивные запросы, необходимо изучить на базе какой информации мы можем их производить. Откроем атрибутивную таблицу слоя, щелкнув правой кнопкой мыши на слое и выбрав Open Attribute Table. В атрибутивной таблице слоя les.shp содержится следующая информация:

    Название поля

    Тип поля

    Описание поля

    Kvartal

    Числовой

    Номер квартала

    Vydel

    Числовой

    Номер выдела

    Area

    Числовой

    Площадь выдела

    Main

    Текстовый

    Тип преобладающей породы

    Age

    Числовой

    Возраст преобладающей породы

  3. Выберем все выдела площадь больше 50 га. Для этого, выберем из меню Selection пункт Select by Attributes... Layer: les (выборка производится из слоя les) Method: Create new selection (создается новая выборка, если до этого было что-то выбрано, то выделение его исчезнет) Два раза щелкнем на поле Area, его название в кавычках добавится в поле запроса внизу, выберем знак "больше" (>), и впишем значение 50. Целиком наш запрос должен приобрести следующий вид: "Area" > 50

  4. После ввода всех параметров нажмите Ок.

  5. Обратите внимание как изменился вид. В лесоустроительной карте выбрались некоторые выдела. Переключитесь в таблицу, чтобы определить сколько записей (т.е. выделов выделилось).

  6. Допустим нам обходимо далее выбрать из всех выбранных выделов только выдела с елью в качестве преобладающей породы. Опять выбираем из меню Selection пункт Select by Attributes... Layer: les (выборка производится из слоя les) Method: Select from current selection (выбрать из существующей выборки, перед выполнением такого запроса, должна существовать некая предыдущая выборка) Два раза щелкнем на поле Main, его название в кавычках добавится в поле запроса внизу, выберем знак "равно" (=) и нажмем на кнопку Get Unique Values (Получить уникальные значения). Эта кнопка , для выделенного поля осуществит поиск всех уникальных значений и выведет их в списке справа. Из этого списка можно будет потом выбрать нужное нам значение и щелкнув на нем 2 раза, вставить в наш запрос. Проделаем это для значения "Ель". Целиком наш запрос должен приобрести следующий вид: "Main" = 'Ель' Нажмем Ок.

  7. Как видно в окне фрейма данных часть выделов, которые были выделены перестали быть выделенными, это те выдела, которые не соответствуют нашему последнему запросу. Открыв таблицу можно определить, сколько осталось выделенных выделов.

  8. Мы произвели выделение нужных нам выделов в 2 этапа, а можем произвести и в один. Для этого, снимем выделение с выделенных объектов. Selection\Clear Selected Features

  9. Опять выбираем из меню Selection пункт Select by Attributes... Layer: les (выборка производится из слоя les) Method: Create new selection (создается новая выборка, если до этого было что-то выбрано, то выделение его исчезнет) Далее либо выбирая соответствующие поля и значения, либо просто вписав запрос в поле ввода запроса получаем следующее выражение: "Area" > 50 AND "Main" = 'Ель' Нажимаем Ок.

  10. Мы получили в один шаг результат, который до этого получили за два.

  11. Создадим из полученной выборки виртуальный слой. Для этого правой кнопкой мыши щелкнем на слой содержащий выбранные объекты, в нашем случае это слой les и выберем Selection\Create Layer from Selected Features

  12. Во фрейм данных будет добален новый слой с именем les selection, в котором будут содержаться только выделенные объекты. Полученный слой является виртуальным и существует только в данном проекте, если мы закроем наш проект без сохранения, то наш результат пропадет.

  13. Для того, чтобы сохранить результат в виде отдельного слоя, существующего в виде соответствующего набора файлов на диске, проделаем следующую операцию. Правой кнопкой мыши щелкнем на слой содержащий выбранные объекты, в нашем случае это слой les и выберем Data\Export Data... В появившемся окне в пункте Export (Экспорт) должно стоять Selected Features (выбранные объекты) что означает, выделенные объекты из слоя les будут преобразованы в новый слой. В поле Output shapefile or feature class (выходной слой) можно задать имя и путь к новому слою.

Упражнение 2. Пространственные запросы, определение количества и номенклатуры необходимых для исследуемой территории топокарт масштаба 1:200000.

  1. Создайте новый проект

  2. Загрузите слой oblasts.shp

  3. Загрузите слой grid2km.shp. Этот слой является стандартной разграфкой топографических карт выпускаемых Военно-Топографическим Управлением Генерального Штаба РФ. Если вы работаете с картами этого производителя, то с помощью такой разграфки вы легко можете определить количество и номенклатуру нужных вам топографических карт. Такие разграфки можно получить здесь. В атрибутивной таблице таких слоев для каждого топографического листа обычно значится его стандартное название, согласно номенклатуре.

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

  5. Создайте атрибутивную выборку, которую вы научились делать в предыдущем упражении, выберите из слоя oblasts.shp Челябинскую область. Selection\Select by Attributes...

  6. После того, как необходимый объект был выделен, увеличьте его с помощью инструмента Zoom to Selected Selection\Zoom to Selected

  7. Произведем пространственную выборку. Выберем из меню Selection пункт Select by Location... (выбрать по местоположению)

  8. Зададим условия выборки на человеческом языке звучащей как: Выбрать из темы grid2km (разграфка карт 2км масштаба) все объекты пересекающиеся с выделенными объектами темы oblasts (области). Для этого выберем: I want to: select features from (Я хочу выбрать из) the following layers: grid2km (из следующих слоев: поставим галочку напротив слоя grid2km (разграфка карт 2км масштаба)) that: intersect (которые: пересекают) the features in this layer: oblasts (объекты в этом слое: oblasts (области)) Переключатель use selected features (использовать выделенные объекты) под слоем oblasts должен быть включен (установлена галочка), иначе выделены будут объекты пересекающиеся с любыми объектами слоя областей, а не только выделенной нами областью. В результате окно настроек должно выглядеть так:

  9. После того, как все условия введены, нажмите Apply (применить).

  10. Снимите выделение для слоя областей, для этого, выделите слой oblasts и выберите из меню Selection\Clear Selected Features (Снять выделение)

  11. Откройте таблицу, вы увидите, что в ней выделились объекты являющиеся топографическими листами на нужную вам территорию. По атрибутивной таблице вы можете определить сколько вам необходимо стандартных топографических листов (т.е. сколько выделилось объектов) и что это за листы (по полю Name)

Упражнение 3. Визуальная выборка, скрывание объектов по атрибутивной информации

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

Работа с табличными данными

<<< предыдущая глава | оглавление | следующая глава >>>

Теоретическая часть

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

  • Целочисленные (integer)

  • Десятичные (float)

  • Текстовые (text)

  • Логические (boolean)

  • Дата (Date)

Редактирование данных

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

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

Установление связи между таблицами

Связи между таблица бывают двух видов:

1. Соединение

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

2. Связывание

Также как и соединение, связывание (Link) устанавливает связь одна-ко-многим между таблицей назначения и таблицей источника. Одна запись в таблице назначения связывается с одной или несколькими записями в таблице источнике. При связывании таблиц добавления данных из одной таблицы к другой не происходит (в отличие от операции соединения). Между ними только устанавливается связь, поэтому выбирая запись в одной из таблиц Вы автоматически выбираете запись в связанной с ней таблице.

Пространственное соединение таблиц

Если в качестве общего поля используется поле Shape, соединение происходит по принципу пространственной связи между объектами в соединяемых темах. Эта операция называется пространственным соединением. Пространственное соединение подобно соединению по атрибутам, то есть в результате соединения атрибуты из одного слоя присоединяются к атрибутам другого связанным по записям; однако основано оно на пространственной связи между объектами в двух темах. Для каждого объекта, представленного в результирующей таблице, осуществляется поиск наличия пространственных связей с любым объектом, представленным в исходной таблице. Если такие связи существуют, то запись объекта из исходной таблицы присоединяется в результирующую таблицу. Соединяться могут слои разных типов, соединение может происходить по следующим принципам:

 

точка

линия

полигон

точка

близость

близость

внутри

линия

близость

часть от

внутри

полигон

-

-

внутри

Когда пространственное соединение основано на связи типа "близость" (один из двух связываемых слоев должен содержать точки и ни один из слоев не должен быть полигональным), связываются объекты находящиеся на кратчайшем расстоянии друг от друга. Расстояние вычисляется в Виде в картографических единицах.

Когда пространственное соединение основано на связи типа "часть от" (оба слоя должны быть линейными), линейные объекты, представленные в результирующей таблице связываются с линейными объектами, представленными в исходной таблице, если линия в исходной таблице является частью линии другого слоя. Например, если первая тема содержит реку Ока, а вторая - все реки бассейна Волги, то все объекты в первой теме рассматриваются как "часть от" объектов второй темы, потому что река Ока - приток Волги, и объекты, представляющие Волгу, присутствуют и во второй теме.

Работа с таблицей в процессе создания макета карты В макете карты таблица может использоваться как графический элемент. Она может быть добавлена в любое место карты, масштабирована, связана с исходным документом Word или Excel, связь позволяет открывать и редактировать таблицу не покидая режима верстки карты.

Практическая часть

Скачать учебные материалы для этой главы (0.1 Мб).

В практической части этой главы Вы научитесь:

  • Работать с атрибутивной таблицей слоя - открывать ее, выделять объекты;

  • Добавлять и удалять новые записи в таблице;

  • Создавать новые поля в таблице;

  • Записывать значения в созданное поле.

Упражнение 1. Работа с табличными данными

  1. Откройте проект ex12.mxd из папки chapt17. В этом проекте находится только один слоя - линейный слой дорог.

  2. Открйте атрибутивную таблицу слоя (нажмите левой кнопкой мыши на названию слоя в таблице содержания вида и выберите Open Attribute Table) В атрибутивной таблице выделите несколько строчек щелкая мышкой по самому правому столбцу таблицы. Удерживая клавишу crtl выделите одновременно несколько строк. Обратите внимание, что в виде все объекты соответствующие выделенным в таблице строкам также выделились

  3. Начните редкатирование слоя дорог (Editor/Start Editing). Создайте (нарисуйте) новый объект в виде. Обратите внимание, что в таблице добавилась новая пустая строка соответсвующая созданному только что объекту. В поле TYPE укажите тип дороги - грунтовая Закончите редактирование слоя без сохранения изменений (Editor/Stop editing).

  4. Откройте атрибутивную таблицу слоя road-l. Нажмите на кнопку Options и выберите Add Field.. В окне создания нового поля наберите название слоя Subclass и задайте тип поля Text с длиной поля 100 Нажмите ок. В таблице добавилось новое поле Subclass. Левой кнопкой мыши щелкните по названию поля Subclass и выберите Delete Field

  5. Запустите ArcCatolog. В окне дерева каталога найдите папку chapt17/vector. В окне Contents щелкните левой кнопкой мыши по названию слоя road-l и выберите из списка Properties В окне Shapefile Properties в закладке Fields просмотрите какие поля, содержаться в атрибутивной таблице слоя и какой тип (Data Type) имеет каждое поле. В одной из пустых строк столбца Field Name попробуйте написать название нового столбца - как Вы увидите, это сделать невозможно, так как в данный момент слой загружен в открытый проект

  6. Перейдите в окно ArcMap и удалите слой road-l из таблицы содержания вида. Затем опять переключитесь в окно ArcCatalog, откройте окно Shapefile Properties слоя road-l и создайте новое поле Subclass в столбце Field Name. Тип поля (Data Type) выбирите Text, длину поля задайте равной 100. Нажмите Apply и загрузите тему в проект ex12.mxd. Откройте атрибутивную таблицу слоя. Как Вы видете в таблице повился новый столбец Subclass.

  7. Начните редактирование слоя road-l. Нажмите кнопку Options и выбирите Select By Attributes.. В окне построения запроса введите выражение "TYPE" = 'проселочная'. Таким образом, все дороги этого типа будут выбраны Нажмите Apply

  8. Левой кнопкой мыши щелкните по названию столбца Subclass в атрибутивной таблице. В раскрывшемся списке выберите Calculate Values... В окне Field Calculator в поле ввода значения "Subсlass =" наберите слово "ремонт" (обязательно в двойных кавычках) и нажмите ок Как Вы видите всем дорогам типа "проселочная" в поле Subclass было присвоено значение "ремонт". Закончите редактирование слоя без сохранения изменений.

  9. Довольно часто может возникать необходимость присоединения таблицы одного слоя к таблице другого слоя. Для того, чтобы осуществить такое присоединение, нажмите левой кнопкой мыши по названию слоя в таблице содержания вида, выберите в списке Joins and Relates и затем Join.. В открывшемся окне Join Data укажите то, по какому поля слоя будет проводится присоединение, укажите слой, таблица которого будет присоединена и поле в этой таблице.

  10. Атрибутивные таблицы можно добавлять в макет карты. Для этого откройте атрибутивную таблицу слоя road-l, нажмите кнопку Options и выберите Add Table To Layout. Автоматически таблица будет добавлена в центр макета карты и вы будете переключены в режим работы с макетом карты.

Работа с растровыми данными

<<< предыдущая глава | оглавление | следующая глава >>>

Теоретическая часть

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

Процесс привязки растровых данных состоит из нескольких этапов:

Этап 1 Расстановка точек привязки - т.е. нахождение и точек, координаты которых известных и ввод для них этих координат. Источниками координат могут быть углы и точки пересечения координатной сетки (на топографических картах и любых других картах с координатной сеткой), другие, уже привязанные, изображения, данные, полученные с помощью GPS, любые другие географически привязанные данные.

Этап 2 Проверка точности и правильности расстановки точек - расчет среднеквадратичной ошибки.

Этап 3 выбор метода пересчета значений элементов изображения при трансформации, (resampling, передискретизация) - метод ближайшего соседа, кубической свертки, билинейной интерполяции. Для отсканированных материалов обычно выбирается метод ближайшего соседа.

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

Этап 5 Выбор размера ячейки результирующего изображения (т.е. выбор размера пиксела)

Этап 6 Осуществление трансформации. Используя параметры указанные выше, программа осуществляет пересчет нового местоположения для каждого пиксела\ исходного растра и расчитывает для новых пикселов их новые значения. Процесс может занимать длительное время.

Практическая часть

Скачать учебные материалы для этой главы.

В практической части этой главы Вы научитесь:

  • Привязывать топографические карты по координатам;

  • Привязывать другие растровые данные к уже привязанным материалам.

Упражнение 1. Привязка топографической карты с известными координатами

  1. Запустите ArcMap и начните новый проект.

  2. Добавьте в проект тему chapt18\topo.tif На вопрос-сообщение "Would you like to create pyramid for..." (создать пирамидные слои) ответьте No (нет).

  3. Запустите из окна ArcMap ArcCatalog Tools\ArcCatalog

  4. Найдите в дереве папок каталог chapt18 и щелкнув на нем правой кнопкой мыши выберите New - Shapefile (Новый шейп-файл)

  5. Введите в поле Name (имя) - points (точки), Feature type (тип объектов) должен быть Point (точка).

  6. В разделе Spatial Reference (Пространственная привязка), нажав кнопку Edit (редактировать) и затем кнопку Select (выбрать) выберите географическую систему координат Pulkovo 1942. Нам известно, что наша топографическая карта, как и другие подобные сделана именно в этой системе координат, поэтому мы ее и выбираем. Нажав Ок мы получим новый, пустой слой.

  7. Добавим созданную тему во фрейм данных ArcMap File\Add Data... Выберите тему chapt18\points.shp

  8. Начнем ее редактирование. На панели редактирования нажмем на кнопку Editor (редактор) и выберем Start Editing (начать редактирование). Если панели Editor (редактор) в вашем ArcMap нет, включите ее. View\Toolbars\Editor

  9. С помощью инструмента Zoom In (Увеличительное стекло) увеличим фрагмент карты с известными координатами, например верхний левый угол.

  10. Выберем инструмент редактирования и щелкнув правой кнопкой мыши, выберем из выпадающего меню Absolute X,Y...

  11. В появившемся окне, введите координаты нужной вам точки и нажмите клавишу Enter. Значения широты и долготы должны вводиться в формате десятичных градусов, но не градусы-минуты-секунды. Если на вашей карте значения координат градусы-минуты-секунды, то вам необходимо перевести их в десятичные градусы. После ввода, точка исчезнет из вида, не обращайте на это внимание.

  12. Вернемся к общему виду щелкнув правой кнопкой мыши на слое карты (имя слоя) и выбрав Zoom To Layer (Увеличить до полного охвата слоя).

  13. Повторив шаги 9-12, введем таким же образом три остальные угловые точки с привязываемой топографической карты. Можно сказать, что в этом случае источником координат является сама карта.

  14. После окончания ввода координат, если вам известна система координат растровых данных которые вы привязываете, а в нашем случае она известна, задайте фрейму данных спроектированную систему координат: Predefined - Projected Coordinate Systems - Gauss Kruger - Pulkovo 1942 - Pulkovo 1942 GK Zone 15

  15. Щелкните правой кнопкой мыши на слое точек (points.shp) и выберите Zoom To Layer (Увеличить до полного охвата слоя). Слой точек, который мы создали отобразится на экране.

  16. Откройте панель Georeferencing (Географическая привязка) View\Toolbars\Georeferencing

  17. На появившейся панели из меню Georeferencing (Географическая привязка) выберите пункт Fit to Display (Вписать в экран), эта функция примерно впишет привязываемый растр в уже показываемый слой точек (см. пункт 15)

  18. Далее наша задача будет состоять в том, чтобы более точно "подтянуть" растр к точечному слою. Для этого будем использовать инструмент Add Control Points (добавить линию привязки) находящийся на панели инструментов Georeferencing (Географическая привязка). Выберем этот инструмент.

  19. С помощью инструмента Zoom In (Увеличительное стекло) увеличим фрагмент карты с известными координатами, например верхний левый угол и одновременно наиболее близкую к нему точку, представляющую реальное положение этого угла. Затем, щелкнув на соответствующем пересечении линий растра протянем к точке линию привязки. Мы увидим как сместился растр.

  20. Вернемся к полному охвату точечной темы щелкнув правой кнопкой мыши на слое карты (имя слоя) и выбрав Zoom To Layer (Увеличить до полного охвата слоя).

  21. Проделаем операции 19-20 еще 3 раза, для остальных углов карты.

  22. После того, как карта "встала" на место, то есть привязалась к нужным координатам (=точкам), можно оценить математическую точность привязки, выбрав инструмент View Link Table (посмотреть таблицу линий привязки) . В этой таблице, помимо списка исходных и конечных точек и относительных ошибок каждой из них, можно также посмотреть общую среднеквадратичную ошибку (Total RMS Error) в единицах карты. В нашем случае это должна быть ошибка порядка 50-70 метров.

  23. Непосредственно процесс привязки (расстановки контрольных точек) завершен, осталось только сохранить измененную карту под новым именем. Georeferencing\Rectify... Cell Size: 45 (размер ячейки выходного растра в метрах) Resample Type: Nearest Neighbor (тип интерполяции значений пикселов - метод ближайшего соседа) Output Raster: (имя новой привязанной карты - введите новое имя и поместите карту в папку chapt18)

  24. Процесс займет некоторое время, после создания нового, привязанного, растрового слоя, он будет добавлен во фрейм данных.

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