Эта статья посвящаеться вопросам касающихся привязки подробных топографических карт. Постараемся узнать суть привязки топографических карт, для чего это нужно и как ее выполнить в ArcGis, Erdas Imagine, Envi. Для чего нужны подробные топографические карты с привязкой? Какую полезную информацию могут они нести? Интересные вопросы, на которые мы сейчас ответим. Итак, для того чтобы можно было использовать производную информацию с топографической карты или саму топокарту в уже готовых или новых проэктах ГИС, нужно чтобы топокарты были привязаны(имели ситему координат). Нами была составленна подробная инструкция по привязке топографических карт в 3-х програмных продуктах
ArcGis, Erdas Imagine, Envi.
С топокарты можно получить следующую информацию:
Растровая топографическая карта района, области, интересующей територии;
Вектор изолиний, с которого в дальнейшем можно получить ЦМР;
Вектор дорог, населенных пунктов, лесных насаждений, водных объектов и т.д.
Процес привязки топографических карт основан на том, чтобы назначить определенным точкам на карте(опорным точкам) известные координаты.
На каждой топокарте есть уже пять точек с известными координатоми(4 по краям и 1 посредине), которые в первую очередь будут использоваться в процессе привязки. Остальные координаты точек можно вычислить самому, применив при этом небольшие знание в матиматике и понимание топокарты. Например: если координаты левого угла листа топокарты 22°30´, а правого угла 23°00´, то координаты точки расположенной ровно посредине будут 22°45´.
Первым делом открываем новый Viewer, и загружаем в него топографическую карту, которая нуждается в привязке.
Во вьювере запускаем модель географической привязки командой Raster|GeometricCorrection …, где нужно выбрать тип трансформирования – полиноминальный и нажать ОК.
После этого появятся два дополнительных окна : GeoCorrectionTools(Инструменты гео коррекции) и PolynomialModelProperties (Свойства полиноминальной модели).
Geo CorrectionTools
- PolynomialModelProperties - вызов окна Polynomial Model Properties
-GCP-Tool- вызов окна GCP-Tool
- ResampleIcon - выбор метода трансформирования
- CalibrateImage - калибровка изображения.
В окне PolynominalModelProperties на вкладке Parameters выбираем порядок полинома равным 2, на вкладке Projection нажимаем Add/ChangeProjection…, для того чтобы выбрать систему координат в которой будет производится привязка карты.
Задаем систему координат с следующими параметрами
Projection Type– Geographic (Lat/Lon)
Spheroid Name – Krasovskiy
Datum Name –Pulkovo 1942
Если вы планируете привязывать несколько топографических карт, то желательно будет сохранить свойства полиноминальной модели, чтобы каждый раз их не указывать.
После сохранения или указания системы координат жмем Close.
Сразу после этого появляется окно GCPToolReferenceSetup, в котором нужно сделать выбор откуда вы желаете проводить чтение опорных точек привязки. В нашем случае, нужно выбрать пункт KeyboardOnly и нажать Ok .После этого перед вами откроется окно Reference Map Projection, где вы можете отредактировать установленную проекцию. Нажимаете ОК.
Автоматически будут открыты следующие окна и панели:
Viewer с загруженной топокартой ;
Magnifier (Окно с увеличенным масштабом отображения карты и связано с основным окном);
панель GCP-Tool, с помощью которой будет осуществляться ввод опорных точек.
В начале процесс привязки, кнопки Resamle и Calibrate будут недоступны, так как для выполнения этих процессов требуется некоторое количество опорных точек.
Панель GCP-Tool представляет собой таблицу со следующими основными полями:
Point ID - идентификационный номер точки
Color - цвет точки
X, Y input - координаты точки в системе файловых координат входного изображения
X, Y Reference – координаты точки в выходной системе координат , в данном случае в географической системе координат.
Type - тип точки (контрольная - control или проверочная - check)
XResidual, YResidual – ошибки привязки по оси X и Y
RMS Error - вклад в общую RMS-погрешность
Сам процесс привязки подробных топографических карт не вызывает никаких трудностей.
В первую очередь на панели инструментов GCP-Tool необходимо отключить кнопки и , таким образом в процессе набора точек не будут использоваться режимы предсказание положения точек и автоматический пересчет RMS-погрешности, после изменения положения какой-либо точки.
Ну а дальше, увеличиваем во Viewer тот фрагмент топокарты, в которой желаем поставить опорную точку. Желательно работать с теми участками топокарты, на которых можно узнать координаты опорной точки (на топокарте такими точками являются углы карты, пересечения параллелей и меридианов ). Затем нажать кнопку (создание GCP-точки) и поставить точку во вьюере на найденном пересечении. После каждого создания GCP-точки в столбцах X input, Y input будут проставлены координаты точки в пиксельной системе координат. Теперь в столбцах X, Y Reference нужно ввести координаты точки через пробел.
Например:
X Reference= 25 30
Y Reference = 50 20
В процессе ввода и редактирования точек привязки необходимо руководствоваться следующими принципами:
-Опорные точки должны быть размещены на равноудаленных растояних, и должна также присутствовать как минимум одна точка по средине карты
-Количество опорных точек выбирается с тех предположений, что бы обеспечить минимальную RMS-ошибку.
Для полиномиального преобразования 2-й степени минимальное количество точек равно 6. Но для получения хорошего результата число точек должно быть больше, например 21 точка. Но если несколько точек вносят существенную ошибку привязки (большое СКО), то их нужно либо перепривязать либо удалить. Одним словом нужно добиться приемлемой величины СКО.
После того, как все точки будут набраны, необходимо проверить правильность проставленных точек привязки и введенных значений координат. Для этого нужно рассчитать RMS-ошибку нажав кнопку на панели инструментов GCP-Tool. После этого в верхней части панели GCP-Tool появится рассчитанная ошибка.
Среднеквадратическая ошибка метода полиноминального трансформирования для нормальной привязки не должна превышать значение 0,55-0,7.
Последняя стадия привязки топокарт- это определение метода трансформирования (Resamplemethod). Для этого на панели Geo Correction Tool нажимаете Resample. Перед вами откроется окно Resample
Здесь необходимо ввести имя выходного файла и выбрать один с методов трансформирования изображения Resample Method. Также, чтобы в процессе расчета статистики, нулевые значения пикселов были проигнорированы нужно поставить галочку рядом с полем Ignore Zero in Stats.
Дальше жмем ОК и ждем выполнения процесса привязки.
Добавьте растр в таблицу содержания, не обращая внимая на сообщение, что у изображения отсутствует информация о географической привязке
Увеличьте растр так, чтобы точно поставить точку
С помощью инструмента Add Control Points создайте необходимое количество точек привязки (для более точного указания координат используйте инструмент View Link Table )
Разместите курсор в углу карты и щелкните левой кнопкой мыши. Цвет перекрестья сменится, и не двигая мышку, нажмите на правую кнопку. Появится контекстное меню> Input DMS of Lon and Lat
Введите координаты точек привязки в поля Latitude=50040’ Longitude=28000’
ОК
После ввода координат первой точки растр может пропасть с экрана. В этом случае нажмите на кнопку Полный экстент
Посмотреть данные о точке
Привязка топографической карты:
Выбрать Georeferencing
Выбрать Rectify
Выбрать метод привязки –Resample Type – Nearest Neighbor
Выбрать каталог и задать имя привязанной топографической карте