Ортотрансформирование снимков

МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ УКРАИНЫ
Харьковская национальная академия городского хозяйства.
МЕТОДИЧЕСКИЕ УКАЗАНИЯ
К ВЫПОЛНЕНИЮ КУРСОВОЙ РАБОТЫ ПО
ФОТОГРАММЕТРИИ И ДИСТАНЦИОННОМУ
ЗОНДИРОВАНИЮ
(для студентов 4 курса дневной формы обучения
спец. 7.070908
«Геоинформационные системы и технологии»).
Харьков – ХНАГХ-2006
1
Методические указания к выполнению курсовой работы по
фотограмметрии и дистанционному зондированию (для студентов 4
курса дневной формы обучения спец. 7.070908 «геоинформационные
системы и технологии»). Сост. Лисицин В.Э. Харьков: ХНАГХ, 2006.
– 90 с.
Составитель В.Э. Лисицин
Рецензент: профессор, канд. техн. наук В.Д. Шипулин
Рекомендовано кафедрой геоинформационных систем и геодезии,
протокол N4 от 02.12.2005 г.
2
Содержание
ЗАДАНИЕ К КУРСОВОЙ РАБОТЕ………………………………..…....4
ГЕОМЕТРИЧЕСКОЕ ТРАНСФОРМИРОВАНИЕ ОДИНОЧНОГО
СНИМКА…………………………………………………………………..6
Полиномиальные преобразования……………………...………...……9
Трансформирование………………………………………………....10
Расчет потребного числа опорных точек при трансформировании
снимка полиномами различных порядков………………………...13
Недостатки нелинейного трансформирования……………………14
Оценка среднеквадратической ошибки трансформирования……16
Контрольные точки (Check Points)……………………………...…18
ПРАКТИЧЕСКОЕ ВЫПОЛНЕНИЕ КООРДИНАТНОЙ ПРИВЯЗКИ
И ГЕОМЕТРИЧЕСКОГО ТРАНСФОРМИРОВАНИЯ СНИМКОВ….18
ОБРАБОТКА ФОТОГРАММЕТРИЧЕСКИХ БЛОКОВ ПРИ
ПОМОЩИ МОДУЛЯ ORTHOBASE………………………………….27
Общие понятия…….......……………………………………………....28
Параметры, регулирующие процесс автоматического
генерирования связующих точек в OrthoBase…………………....28
Основные схемы аэрокосмического фотографирования………...31
Выбор геометрической модели съемочной аппаратуры…….…..32
Системы координат используемые в OrthoBASE…………….….33
Элементы внутреннего ориентирования……………………….…36
Элементы внешнего ориентирования………………………….….39
ВЫПОЛНЕНИЕ ОРТОТРАНСФОРМИРОВАНИЯ И ПРОВЕРКА
ТОЧНОСТИ СТЕРЕОМОДЕЛИ ДЛЯ СНИМКОВ ПОЛУЧЕНЫХ
ПРИ ПОМОЩИ КАДРОВОЙ КАМЕРЫ…………………………….…42
Подготовка файла блока………………...………………………….…43
Ввод опорных и контрольных точек……....……………………….....54
Выполнение автоматического сбора связующих точек……....…....69
Выполнение аэротриангуляции……………………………..……......70
Проверка точности стереомодели……………..................................73
Ортотрансформирование снимков…..………………………………79
СПИСОК ЛИТЕРАТУРЫ……………………………..…………………82
ПРИЛОЖЕНИЕ 1. ВАРИАНТЫ ЗАДАНИЙ ДЛЯ ВЫПОЛНЕНИЯ
ГЕОМЕТРИЧЕСКОГО ТРАНСФОРМИРОВАНИЯ СНИМКА………83
ПРИЛОЖЕНИЕ 2. ВАРИАНТЫ ЗАДАНИЙ ДЛЯ ВЫПОЛНЕНИЯ
ОРТОТРАНСФОРМИРОВАНИЯ СНИМКА…………………………..83
ПРИЛОЖЕНИЕ 3. НАЗВАНИЯ КОНТРОЛЬНЫХ И ОПОРНЫХ
ТОЧЕК ДЛЯ ОРТОТРАНСФОРМИРОВАНИЯ / НАЗВАНИЕ
ФРАГМЕНТА СНИМКА……………………………….………………..86
3
Задание к курсовой работе
Курсовая работа состоит из двух частей. В первой части,
используя инструментарий OrthoBase, обучающиеся выполняют
геометрическое трансформирование предложенного снимка. Отчет по
первой части включает следующие пункты.
1. Информация по исходному (подлежащему трансформированию)
снимку.
2. Краткая характеристика исходного снимка (типы рельефа,
характерные области, присутствующие на снимке). Схематическое
изображение нетрансформированного снимка. Указание на схеме
опорных и контрольных точек. Схематическое (или при помощи
Coverage либо Shape - слоя) изображение характерных областей
рельефа на исходном снимке.
3. Информация по опорному снимку.
4. Информация по трансформированному снимку.
Имя
трансформированного снимка, имя файла модели, имена INPUT и
REFERENCE-файлов.
5. Выбор и обоснование на основе характеристик исходного снимка
модели трансформирования.
6. Система полиномов выбранной модели трансформирования (с
подстановкой реальных коэффициентов) в виде:
xисх1  A0  A1* xt1  A2 * y t1  ...
yисх1  B0  B1* xt1  B2 * yt1  ...
.............................................
xисхN  A0  A1* xtN  A2 * y  ...
(1)
tN
y исхN  B0  B1* xtN  B2 * y tN  ...
Примечание. Для проверки качества трансформирования в ERDAS
используется обратное преобразование, поэтому,
x ti , y –
ti
географические координаты на опорном снимке а
x исхi , yисхi -
пиксельные координаты на исходном снимке.
7. Обоснование выбора посадочных мест для контрольных и
опорных точек, используемых при трансформировании снимка и
4
контроле результатов (по каждой точке). Например - опорная точка
№1 выбирается на углу строения и т.д.
8. Результирующая таблица, содержащая данные автоматического
расчета и точки.
9. Трансформационная матрица.
10. Ручной расчет невязок для трех контрольных точек Xисхi, Yисхi.
11. Выводы.
а) оценка качества трансформирования;
б) анализ качества трансформирования характерных областей
рельефа, присутствующих на снимке;
в) возможности улучшения качества трансформирования;
Во второй части работы Вам предлагается три снимка. При
помощи модуля OrthoBase определяется модель съемочной камеры,
составляется съемочный блок, выполняется измерение опорных и
контрольных точек, а также автоматический сбор связующих точек.
Выполняется триангуляция, затем при помощи модуля StereoAnalyst
проверяется точность созданной модели. И, наконец, выполняется
ортотрансформирование снимков. Отчет по второй части работы
включает следующее.
1. Принципиальная схема фотографирования (кадровое или щелевое).
2. Категория, модель и основные свойства модели для аппаратуры
съемки. Элементы внутреннего ориентирования - F, x0, y0,таблица
координат для координатных меток. Схема матрицы координатных
меток, система координат координатных меток. Схема соотношения
координатных меток и параметров внутреннего ориентирования.
3. Схема и свойства блока изображений (сколько файлов, их общая
характеристика - , количество каналов (зон) отображаемые объекты)
4. Сводная таблица начальных значений параметров внешнего
ориентирования для каждого снимка.
5. Сводная таблица опорных и контрольных точек - номер пары
снимков, тип точки, координаты пиксельные, координаты
географические. Схематическое расположение контрольных точек.
6. Параметры автоматической генерации связующих точек.
7. Отчет аэротриангуляции. Комментарии к отдельным блокам отчета
и их анализ (обоснование значений параметров внешнего
ориентирования, невязки опорных и контрольных точек).
8. Сводная таблица окончательных значений параметров внешнего
ориентирования для каждого снимка.
9. Проверка качества блока в стереомодели с использованием двух
контрольных точек.
5
10. Параметры ортотрансформирования (имя внешнего DEM-файла,
размер ячейки выходного файла, метод передискретизации - Nearest
Neighbor и т.д.).
11. Анализ качества ортотрансформирования, выполненный при
помощи инструмента "шторка".
12. Выводы.
Геометрическое трансформирование
одиночного снимка
Авиационные и спутниковые снимки несут информацию о части
земной поверхности. Даже если рельеф данной области земной
поверхности относительно равномерный и искажения, вызванные
рельефом на снимках минимальные, приходиться принимать во
внимание искажения вызванные учетом кривизны Земли,
погрешности, вносимые в процессе сканирования и искажения,
вносимые съемочной аппаратурой. Ниже даны сведения о процессе
геометрического
трансформирования
(rectifiсation)
исходных
изображений, минимизирующем вышеуказанные искажения.
Подчеркнем, что при геометрическом трансформировании мы не
учитываем ошибки снимка, вызванные рельефом местности. Для этого
существует
другая
разновидность
трансформирования
ортотрансформирование
(orthorectification)
–
разновидность
трансформирования корректирующего на снимке ошибки вызванные
рельефом местности. Такая коррекция выполняется только при
наличии цифровой модели местности (Digital Elevation Model – DEM).
Для
областей
с
относительно
равномерным
рельефом
ортотрансформирование не требуется, но при обработке снимков,
например горных районов или городских зданий, существенно
повышает точность обработки. В некоторых случаях, если DEM на
данную территорию отсутствует, допускается задание высоты в виде
константы.
Геометрическое трансформирование
следующих случаях:
снимков
необходимо
в
- выполнение по снимкам точных измерений расстояний и площадей;
- сравнение изображений с существенно различающимися масштабами
съемки;
6
- создание точных карт на основе фотографического материала;
- наложение растровых и векторных данных;
- создание мозаики изображений (склейка областей из разных
снимков);
Геометрическое трансформирование снимков не требуется, если
отсутствует искажение изображения. Например, если изображение
получено сканированием бумажной карты, имеющей известную
картографическую
проекцию,
то
обычно,
дальнейшее
трансформирование не требуется (если только при сканировании не
было привнесено некоторое смещение или разворот).
Введем понятие опорных точек (Ground Control Points - GCP).
Это пиксели изображения, для которых известны наземные
координаты в единицах проекции карты. Любой снимок в формате
ERDAS IMAGINE может иметь один набор GCP, ассоциированных с
ним. Этот набор точек, если он существует, храниться внутри imgфайла, наряду с другой информацией. Фактически, каждая GCP
описывается двумя парами X,Y-координат:
- исходные координаты – обычные файловые координаты внутри
снимка подлежащего трансформированию;
- опорные (референц-)координаты – координаты в единицах проекции
карты, в которой зарегестрирован снимок;
В случае существования встроенного набора опорных точек, их
можно просматривать в Визуализаторе при помощи инструментария
GCP Tool. При этом в качестве имени опорной точки используется ее
идентификатор - ID.
Набор GCP можно получить в результате полевых измерений, при
помощи GPS, при помощи ортотрансформированного изображения
либо снять с бумажной карты
Однако исходный снимок, подвергаемый трансформированию, как
правило, не привязан к какой-либо картографической проекции. То
есть, мы работаем в файловой системе (начало координат в левом
верхнем углу окна, единицы измерения – пиксели). Где же взять
референц-координаты? Выходную систему координат удобно задавать
при помощи вспомогательного изображения, с известной
картографической проекцией и датумом. Такое изображение будем
называть референц-снимком (или опорным снимком), а его систему
координат – референц-системой.
7
8
Иллюстрация процесса геометрического трансформирования снимка.
Теперь процесс геометрического трансформирования будет содержать
следующие шаги:
- нанесение опорных точек на исходный и опорный снимки;
- вычисление параметров трансформирования (об этом ниже);
-создание выходного файла, пиксели которого подвергнуться
передискретизации, чтобы вписаться в сетку новой координатной
системы;
Для точного трансформирования, следует аккуратно наносить
опорные
точки.
По
ним
будут
интерполироваться
трансформированные координаты остальных точек изображения. GCP
должны быть равномерно распределены по всему снимку, а не
сконцентрированы в какой-то одной области.
В качестве GCP на снимке удобно выбирать перекрестки дорог,
углы зданий, взлетных полос и стадионов. Не следует задавать GCP на
кромках озер, краях зон растительности и прочих объектах,
подверженных изменениям с течением времени.
В
качестве
математического
аппарата
трансформирования
снимков
используют
преобразования.
геометрического
полиномиальные
Полиномиальные преобразования
Полиномиальные уравнения используются для того, чтобы
выполнить преобразование координат снимка из системы исходного
файла в систему координат трансформированного снимка.
По
координатам опорных точек строится математическая модель
трансформирования снимка. В зависимости от искажений снимка,
числа опорных точек и их взаимного расположения для построения
уравнений трансформации могут применяться полиномы различной
степени (1-й, 2-й, а иногда и выше).
При
полиномиальном
преобразовании
над
исходным
изображением могут быть выполнены операции следующего типа:
- тождественное преобразование;
- смещение;
- масштабирование;
- зеркальное отображение;
- поворот;
9
- сдвиг;
- любые комбинации перечисленных операций
Трансформирование
Линейное трансформирование является простейшим образцом
полиномиальных преобразований. При этом применяются полиномы
1-го порядка. Линейное трансформирование может применяться при
перенесении исходного, необработанного снимка в плоскую
картографическую проекцию или при преобразовании из одной
плоской
картографической
проекции
в
другую
плоскую
картографическую
проекцию
или
при
трансформировании
относительно малых площадей снимка.
По набору GCP рассчитывается матрица трансформирования,
содержащая коэффициенты полиномов. Каждая опорная точка вносит
свой вклад в расчетные величины полиномиальных коэффициентов.
Матрица
линейного
трансформирования
содержит
шесть
коэффициентов:
a0
b0
a1
b1
a2
b2
На базе этих коэффициентов построим, например, систему
полиномов первого порядка, используемую для трансформирования
по трем опорным точкам изображения, задаваемым координатами
(xi,yi):
xt1  a0  a1* x1  a2 * y1
yt1  b0  b1* x1  b2 * y1
xt 2  a0  a1* x2  a2 * y2
yt 2  b0  b1* x2  b2 * y2
(2)
xt 3  a0  a1* x3  a2 * y3
yt 3  b0  b1* x3  b2 * y3
10
где
x ti , yti – координаты трансформированного изображения;
xi , y i – исходные координаты;
При линейном трансформировании в ERDAS IMAGINE можно
задавать параметры для операций масштабирования, смещения
поворота и зеркального отображения.
На рисунке, приведенном выше, показаны результаты следующих
операций: 1 – исходное изображение; 2 – масштабирование по X; 3 –
масштабирование по Y; 4 – сдвиг по X; 5 – сдвиг по Y; 6 – поворот;
Нелинейное трансформирование выполняется при помощи
полиномов второго порядка и выше. При этом могут быть
скорректированы и нелинейные искажения. Трансформирование при
помощи полиномов 2 – го порядка может быть использовано для
преобразования данных в координатной системе широта/долгота в
плоскую картографическую проекцию, для преобразования данных
покрывающих большие
площади (учет кривизны Земли),
преобразования искаженных данных (например, по причине дисторсии
линзы съемочной камеры). Полиномы 3 – го порядка и выше
используются для трансформации сильно искаженных снимков.
11
Пример системы уравнений с применением полиномов третьего
порядка для трансформирования по одной опорной точке
изображения, задаваемой координатами (x,y):
xt  5  4 x  6 y  10 x
2
 5xy 
y
2
2
 3 x  7 x y  11 x y  4 y
3
2
3
(3)
yt  13  12 x  4 y  x
2
2
2
 21 xy  11 y  x  2 x y  5 x y  12
3
2
y
3
Ниже представлены некоторые примеры нелинейных трансформаций.
1 – исходное изображение; 2,3,4.5 – выходные изображения;
При
трансформировании
при
помощи
полиномов,
трансформационная матрица содержит для каждой из координат (X и
Y) коэффициенты, количество которых определяется по формуле:
t 1
N  2 * i
(4)
i 1
где N – количество коэффициентов трансформационной матрицы;
t – порядок полинома;
12
Расчет потребного числа опорных точек при
трансформировании снимка полиномами
различных порядков
Чем выше порядок полиномов трансформирования, тем более
сложные искажения могут быть исправлены. Однако с повышением
порядка полиномов требуется большее число опорных точек.
Например, при трансформировании первого порядка (линейном),
полиномы для X и Y описывают уравнение плоскости. А плоскость
однозначно определяется тремя точками. Значит, необходимы, по
крайней
мере,
три
опорных
точки.
Аналогично,
при
трансформировании второго порядка, полиномы для X и Y описывают
уравнение параболоида. А его однозначно определяют шесть точек. И
требуются не менее шести GCP.
К тем же выводам можно прийти, рассмотрев число неизвестных
элементов в системе полиномиальных уравнений. Например,
проанализируем рассмотренные выше полиномы третьего порядка.
2
2
3
2
3
2
xt  5  4 x  6 y  10 x  5xy  y  3 x  7 x y  11x y  4 y
2
2
3
yt  13  12 x  4 y  x2  21 xy  11 y  x3  2 x2 y  5x y  12 y
Здесь x и y – картографические координаты текущей опорной
точки. Они берутся с референц-изображения и значит известны. А
коэффициенты трансформирования a 0  a 9 , b0  b9 (всего 20
коэффициентов) неизвестные величины. Для решения такой
системы потребуется 20 уравнений (по два на каждую опорную точку).
Следовательно, при трансформировании с помощью полиномов
третьего порядка минимально потребное число опорных точек - 10.
Вообще, минимальное количество опорных точек, необходимое
для выполнения трансформирования определяется как:
N min 
(( t  1) * (t  2))
2
где Nmin – минимальное количество опорных точек;
t – порядок полинома;
13
(5)
Недостатки нелинейного трансформирования
Рассмотрим одномерный случай – изменение координаты X на
исходном и трансформированном снимках. По осям графика
откладываются значения Xc (исходные координаты X) и Xr (опорные
координаты X).
На правом графике показана следующая ситуация. Требуемое
трансформирование невозможно выполнить при помощи полиномов
первой степени. Какие бы коэффициенты не подбирались, полином
первой степени будет задавать уравнение прямой линии и опорная
точка номер 2 (ось Xc) не будет трансформирована правильно.
Требуется применение полиномов более высокого порядка, что и
отображено на рисунке, приведенном ниже.
14
Теперь, полином второй степени отлично справляется с ситуацией
до тех пор, пока не будет добавлена дополнительная четвертая точка.
Повторив рассуждения, приведенные выше, придем к
необходимости использования полиномов третьей степени, что и
отображено на двух рисунках, приведенных ниже.
Из
последних
двух
рисунков
видны
недостатки
трансформирования при помощи полиномов более высокого, чем
первый порядков. Может получиться так, что на трансформированном
изображении изменится порядок следования пикселей. Данное явление
иллюстрируется на двух рисунках, приведенных ниже.
Видно, что порядок точек 3 и 4 на правом (трансформированном)
изображении нарушен.
Такие явления на изображении называются артефактами.
Следует вдумчиво подходить к подбору степени полиномов в процессе
геометрического трансформирования снимков.
15
Оценка среднеквадратической ошибки
трансформирования
В нашем случае среднеквадратическая ошибка (СКО) это
эвклидово расстояние между заданным местоположением опорной
точки на исходном снимке и действительным местоположением той
же точки после геометрического трансформирования. СКО для данной
точки рассчитывается по формуле:
RMSError 
( x t  x i)  ( y t  y i)
2
2
(6)
где RMSError - СКО точки;
xi , y – исходные координаты;
i
xt , y t
– действительные, трансформированные координаты;
СКО выражается в единицах исходной координатной системы.
Например, если опорный снимок записан в файловой системе
координат (данные о картографической проекции отсутствуют), то
СКО будет выражаться в пикселях. Например, СКО = 2 означает, что
точка на референц-изображении отстоит на расстоянии 2 пикселя от
точки на трансформированном снимке.
СКО по определенному направлению (по оси X или Y) называется
невязкой (residual). Невязки вычисляются для каждой точки.
Например невязка по X данной точки это расстояние по оси X между
исходной и трансформированной X-координатами точки. СКО каждой
точки может быть выражено через невязки:
2
2
R i  XR i  YR i
где R i – СКО текущей точки с индексом i;
XRi , YRi – X – невязка и Y – невязка точки с индексом i;
16
(7)
Общее СКО для всех точек можно выразить через невязки
отдельных точек:
Rx 
1 n
 XR i 2
n i 1
T  XR i 2  YR i 2
Ry 
(8)
1 n
 YR i 2
ni 1
(9)
(10)
где R x , R y – общее СКО по X и Y;
T - общее СКО;
n – общее количество опорных точек;
XRi , YRi – X – невязка и Y – невязка точки с индексом i;
Вклад(Contribution) точки с индексом i в СКО:
Ei 
Ri
T
(11)
где Ei - вклад точки с индексом i в СКО;
T - общее СКО;
Ri - общее СКО точки с индексом i;
После того, как вы получили в первом приближении СКО, у Вас есть
четыре возможности.
1. Отбросить все точки, где значение СКО самое большое и
провести новое геометрическое трансформирование с оставшимися
точками. Однако если удаляемые точки единственные в некотором
регионе, то после их удаления трансформирование региона может
привести к еще большим ошибкам.
2. Назначить более высокий допуск на СКО.
3. Усложнить модель трансформирования, подвергая изображение
более сложным преобразованиям.
4. Выбрать и оставить только те точки, которые вы считаете
полностью «надежными».
17
Контрольные точки (Check Points)
Мы определились с минимальным количеством опорных точек,
необходимым для трансформирования при помощи полиномов
степени t. А нужно ли задавать GCP сверх минимального количества?
К чему это приведет? Поскольку все опорные точки принимаются в
расчет при построении полиномов, может случиться так, что
некоторые «лишние» точки плохо впишутся в создаваемую модель.
Это происходит потому, что процесс вычисления коэффициентов
полиномов по своей сути является интерполяционным, а не
аппроксимирующим. То есть дополнительные опорные точки могут не
вписаться в полином степени t, понадобятся полиномы более высоких
порядков.
Тогда, для того чтобы оценить матрицу трансформирования,
следует превратить некоторые опорные точки (GCP) из опорных
(Contгol) в контрольные (Check). Их отличие в том, что контрольные
точки
не
используются
при
вычислении
матрицы
трансформирования. Следовательно, вы можете использовать
контрольные
точки
для
независимой
оценки
точности
трансформирования, просмотрев ошибки контрольных точек (RMS
еггог) и общую среднеквадратическую ошибку по ним (Total).
Практическое выполнение
координатной привязки и
геометрического
трансформирования снимков
1. Из главного меню ERDAS IMAGINE выберите пункт Main |
DataPrep Image | Geometric Correction (Главное | Подготовка данных |
Геометрическая коррекция изображений). В появившемся диалоге
выберите опцию From Image File (Из файла) и выберите имя файла
tmAtlanta.img. Нажмите ОК. Вы выбрали нетрансформированный
космический снимок Landsat ТМ.
18
2. В списке моделей трансформирования выберите Polynomial и
нажмите OK.
3. В верхней части экрана появилась панель с инструментами
геометрической коррекции (Geo Correction Tools), а в центре - диалог
для выбора свойств полиномиального преобразования (Polynomial
Model Propeгties). Используйте полином первого порядка (Polynomial
Order = 1).
19
4. Откройте вкладку Projection (Проекция). Так как проекция для
данного изображения не определена, на этой вкладке нет информации
о какой-либо проекции. Нажмите кнопку Set Projection from GCP Tool
(Установить проекцию с помощью инструмента работы с опорными
точками). Открывшийся диалог GCP Tool Reference Setup
предоставляет возможность выбора источника получения опорных
координат. В данном упражнении выполняется привязка снимка к
снимку, поэтому выберите как источник опоры Image Layer (New
Viewer) (Изображение в новом Визуализаторе) и нажмите ОК.
5. Выберите
снимок
panAtlanta.img,
как
референцизображение(reference image), т.е. изображение, к которому вы будете
привязывать снимок tmAtlanta.img, и нажмите ОК. Нажмите ОК в
появившемся диалоге Reference Мар Information.
20
6. ERDAS IMAGINE организует на экране Визуализатор с
опорным изображением и дополнительные Визуализаторы с
увеличенными изображениями основных Визуализаторов, а также
GCP Тооl (Инструмент работы с опорными точками).
7. Нажмите Apply (Применить) и Close (Закрыть) в диалоге
Polynomial
Model
Propeгties
(Свойства
полиномиального
преобразования). Если вы планируете использовать те же параметры
модели преобразования при привязке снимков в будущем, то перед
закрытием диалога сохраните (Save) параметры в файл.
8. Если вы не изменяли установки по умолчанию для редактора
опорных точек (GCP Tool preferences), то при открытии инструментов
работы с опорными точками будет выбрана кнопка
Режим
автоматического редактирования
(Toggle Fully Automatic GCP
Editing Mode ). Этот режим удобно использовать для предварительной
оценки (прогноза) местоположения опорных точек и для уточнения
соответствия опорных точек на исходном (Souгce) и опорном
(Refeгenced) снимке.
9. Перемещайте связующую рамку в обоих Визуализаторах
таким образом, чтобы она покрыла один и тот же опознаваемый на
21
обоих снимках объект (например, хорошо опознаваемый перекресток
дорог). Теперь этот объект отобразился в окнах Визуализаторов
подробного вида. На панели инструментов редактора опорных точек
нажмите кнопку
Кеер Current Tool (Продолжать работать с
выбранным инструментом). Затем используйте кнопку
Create GCP
(Создать опорную точку). Теперь щелчком мыши разместите
соответствующие
точки
в
обоих
увеличительных
окнах
Визуализаторов. Как только это сделано, можно перейти к
опознаванию следующего объекта на снимках и использовать для
этого кнопку
Select GCP (она автоматически перетащит
связующую рамку к указанному в Визуализаторе объекту и отобразит
увеличенный фрагмент изображения в увеличительном окне
Визуализатора). После ввода З-х опорных точек, при нажатой кнопке
Toggle Fully Automatic GCP Editing Mode, ERDAS IMAGINE после
ввода опорной точки в одном окне Визуализатора автоматически
определит приблизительное (вероятное) местоположение этой точки
во втором Визуализаторе.
10. Если белые маркеры точек плохо видны на изображениях, вы
можете изменить их цвет. Для этого выделите в таблице одну или
несколько записей, содержащих опорные точки, щелкните в поле
таблицы цвет(Coloг) и выберите новый цвет из контекстного меню.
11. По аналогии с первой, наберите еще 8 опорных точек.
Постарайтесь расположить их равномерно по всему снимку. Как и в
любой другой таблице CellAггay, в Редакторе опорных точек
возможно удаление и редактирование опорных точек. Для удаления
точки выберите соответствующую строку в таблице, а затем,
используйте функцию контекстного меню таблицы Delete Selection
(удалить выбранное).
12. Для того чтобы оценить матрицу трансформирования,
следует превратить некоторые опорные точки (GCP) из опорных
22
(Contгol) в контрольные (Check). Их отличие в том, что контрольные
точки не используются при вычислении матрицы трансформирования.
Следовательно, вы можете использовать контрольные точки для
независимой оценки точности трансформирования, просмотрев
ошибки контрольных точек (RMS еггог) и среднеквадратическую
ошибку по всем ним (Total). Выберите любые четыре точки в таблице
CellAггay (выбранные строки выделяются желтым цветом). В меню
редактора опорных точек выберите Edit | Set Point Туре | Check
(Редактировать | Установить тип точки | Контрольная). Содержимое
столбца Туре (тип точки) изменится.
13. В меню редактора опорных точек выберите Edit | Point
Matching. Откроется диалог задания параметров установки опорных
точек. В данном диалоге установите переключатель Discard Unmatched
Point и введите пороговое значение 0.8 в поле Correlation Threshold.
14. Нажмите кнопку Вычисления ошибки для контрольных точек
(Compute Error For Check Points). Ошибки опорных точек (Contгol
Point Еггог) будут заменены на ошибки контрольных точек (Check
Point Еггог). Если не установлен режим Set Automatic Transformation
Calculation при помощи кнопки
,
то для того, чтобы снова
посмотреть ошибки опорных точек, щелкните кнопку
23
Solve
Geometric Model with Control Points (Рассчитать модель по опорным
точкам) - значение Total RMS изменится.
15. Если среднеквадратическая ошибка слишком велика, вы
можете удалить опорную точку с наибольшим вкладом (поле Contrib в
таблице). Для этого выбирают точку (строку в таблице) и используют
функцию Delete Selection (Удалить выделенное) из контекстного
меню, доступного через нажатие правой клавиши на первом столбце
Point #. Вы можете удалять и создавать новые опорные точки до тех
пор, пока среднеквадратическая ошибка не станет приемлемой.
Трансформирование
полиномом
первого
порядка
требует
использования как минимум З-х опорных точек.
Во время редактирования опорных точек следует иметь в виду,
что большая ошибка точки (RMS еrrоr) указывает лишь на то, что
заданное преобразование не может точно совместить исходную
точку с опорной, что чаще всего указывает на ошибку ввода. Но если
вы абсолютно уверены, что положение точки в обоих Визуализаторах
указано правильно, то это положение должно остаться неизменным,
несмотря на большую величину ошибки в этой точке.
16. Когда вы достигнете желаемой точности опорных точек,
сохраните их. Для этого выберите в меню File | Save Input (Файл |
Сохранить исходные точки), затем - File | Save Reference (Файл I
Сохранить точки на опорном изображении). Можно ответить Yes,
чтобы сохранить опорные точки в структуру файла .img.
17. На панели инструментов геокоррекции (Geo Соггесtiоп
Tools)
щелкните
(Передискретизация)
кнопку
.
Display
Появится
трансформирования.
24
Resample
окно
Image
задания
Dialog
параметров
18. Задайте имя выходного файла (xrectify.img) в каталоге, где у
Вас есть право на запись. Укажите метод билинейной интерполяции
(Bilinear Interpolation) в строке Resample Method (Метод
передискретизации). Включите флажок Ignore Zeros in Stats
(Игнорировать нули при расчете статистики). Нажмите ОК, чтобы
начать трансформирование.
19. Когда процесс трансформирования завершится, нажмите ОК
в окне со строкой состояния процесса. Закройте диалог геокоррекции
(GеоСоггесtiоп dialog Bох) кнопкой Exit (Выход). Если вы вносили
какие-то изменения в редакторе опорных точек после последнего
сохранения, вам будет предложено их сохранить. Нажмите Yes (Да)
для сохранения текущей геометрической модели и назовите ее
xrectify.gms. Закройте все окна, за исключением главного окна ERDAS
IMAGINE.
25
20. Откройте новый Визуализатор. Поместите в него
трансформированный снимок xrectify.img и опорное изображение
panAtlanta.img. Используйте инструмент Шторка(Swipe) для
визуальной проверки качества трансформирования.
26
Обработка фотограмметрических
блоков при помощи модуля
OrthoBASE
В противовес обработке одиночных снимков, OrthoBASE
использует подход называемый блочной триангуляцией – то есть
процесс установления математической связи между снимками,
содержащимися в проекте, моделью камеры/сенсора и поверхностью
земли. Данные, полученные в результате блочной триангуляции,
необходимы
в
качестве
исходной
информации
для
ортотрансформирования, создания ЦММ и стереопары. Таким
образом, OrthoBASE позволяет обрабатывать сотни фотографий и
снимков с использованием небольшого количества опорных точек,
исключая, тем самым, проблемы смещения при создании мозаики
снимков. Это приводит к уменьшению времени, затрат и усилий при
получении географической информации большой точности.
Классическая триангуляция, использует метод сбора точек,
известный как расширение опорной сети. Расширение опорной сети
включает в себя измерение вручную по фотоснимкам наземных точек,
которые находятся в области перекрытия снимков. Для сбора опорных
точек на больших площадях фотограмметрические методы более
эффективны, чем традиционные полевые методы наземной съемки, так
как затраты при использовании наземной съемки очень велики.
27
Общие понятия
Параметры, регулирующие процесс
автоматического генерирования связующих
точек в OrthoBase
Наземные точки, используемые в процессе триангуляции,
называются связующими точками (tie points). В OrthoBASE, в
понятие связуюшие точки вкладывается следующий смысл. Это
точки, наземные координаты которых неизвестны, но положение этих
точек визуально распознается в области перекрытия стереопары.
Для повышения эффективности работы с большими и
сверхбольшими блоками изображений в OrthoBASE используется
режим автоматического сбора связующих точек. Для работы в
режиме автоматического сбора связующие точки должны быть хорошо
распознаваемы на всех снимках. В идеале, они должны иметь хорошо
выраженный контраст по двум направлениям, например, угол здания или
пересечение дорог. Связующие точки должны быть также равномерно
распределены по областям перекрытия блока. Пользователю OrthoBASE
в первом приближении необходимо задать количество связующих точек
на блок. Обычно для блочной триангуляции достаточно девяти
связующих точек на каждом снимке. В процессе автоматического сбора
связующих точек решаются следующие задачи:
- автоматическая компоновка блока - основываясь на анализе
начальных данных, IMAGINE OrthoBASE автоматически определяет
соотношение между соседними снимками в блоке;
- выделение особенных точек в качестве кандидатов на роль
связующих точек;
- перенос точки (переколка) - автоматически подбираются и
идентифицируются особенные точки, имеющиеся в областях
продольного и поперечного перекрытии снимков;
- отбор связующих точек - из всего числа намеченных точек
автоматически отбирается то количество точек, которое изначально
задавалось пользователем;
Выбор и измерение связующих точек весьма трудоемкий и
дорогостоящий процесс. Поэтому в последние годы одним из
главных направлений исследований и развития фотограмметрии
28
стала автоматическая триангуляция, в которой основным
вопросом является накопление связующих точек.
Связующие точки на снимках могут быть получены различными
способами:
- ручная «посадка» точки на нескольких или всех снимках, где
присутствует объект, определяющий положение данной связующей
точки;
- импорт некоторого количества связующих точек из внешнего
источника (например, другого блока);
- автоматическая генерация в пределах заданного количества
связующих точек;
Возникают некоторые вопросы. Что, например, случится, если
после импорта или ручного ввода точек мы запустим процесс
автогенерирования? Как будет при этом осуществляться перенос
(«переколка») точек вдоль маршрута и между маршрутами?
Еще один режим работы – обработка сверхбольших блоков.
Например, блок, содержащий 40 маршрутов и 80 снимков в каждом
маршруте, в итоге будет содержать 3200 снимков. Запустив
автоматическую генерацию связующих точекдля всех снимков сразу,
получим массив из 30000 – 40000 точек. Дальнейший анализ такого
блока будет весьма затруднителен.
Кроме этого, Вы можете осуществлять ручное редактирование,
удаляя «плохие» связующие точки с некоторых снимков.
Для регулирования параметров процесса автогенерирования и
переноса связующих точек в OrthoBase существует возможность
выбора вариантов переноса:
-
-
-
No Transfer, New Points Only – не выполнять перенос вдоль
маршрута и между маршрутами для уже созданных или
импортированных точек. При этом разрешено генерировать и
переносить новые точки. Эта опция установлена по умолчанию.
Transfer Only, No New Points – в процессе автогенерирования,
уже созданные или импортированные точки «перекалываются»
вдоль маршрута и между маршрутами. Но новые точки не
генерируются.
Transfer With New Points –
существующие или
импортированные точки переносятся вдоль маршрута и между
маршрутами. Кроме того, генерируются новые точки.
29
В IMAGINE OrthoBASE был разработан общий метод для
автоматического накопления связующих точек. Он создан для работы
с различными цифровыми снимками, такими как аэрофотографии,
спутниковые снимки, снимки, полученные цифровой камерой.
При
автоматическом
поиске
связующих
точек
на
соответствующих
снимках
OrthoBASE
выполняет
процесс
согласования по площадям - также еще называемый согласованием
по знакам (signal). Этот способ определяет соответствие между двумя
областями снимков по подобию величин яркости. Дадим некоторые
понятия процесса согласования областей снимков по площадям.
Эталонное окно - это исходное окно на первом снимке, которое имеет
постоянное местоположение. Его стороны обычно составляют квадрат
(например, 3х3, 5х5 и т.д.).
Окнами поиска являются окна кандидатов на втором снимке, которые
связаны с эталонным окном. В процессе корреляции проверяется
множество различных окон поиска до тех пор, пока не установится
положение, наилучшим образом согласованное с эталонным окном.
Размер окна поиска (Search size) –. По умолчанию задается размером
21*21 пиксель. Чем больше размер окна поиска, тем больше время
поиска и тем выше вероятность неправильного определения
соответствующих точек. Чем меньше размер окна поиска, тем меньше
точек-кандидатов будет найдено.
Пересекающаяся корреляция – процесс вычисления критерия
сходства между эталонным окном и окном поиска, при котором
подсчитывается коэффициент значения яркости между эталонным
окном и окном поиска:
Параметры распределения точек на снимках при их
автоматическом генерировании (Distribution) – параметры
распределения точек на снимках при их автоматическом
генерировании. Параметр распределения может принимать два
значения.
1. Default Distribution – распределение по умолчанию, при
котором точки распределяются случайным образом на всех снимках
блока. При этом будет использоваться матрица из 5*5 шаблонов
30
области генерации (patterns). То есть связующие точки будут
генерироваться в областях снимка соответствующих 10%, 30%, 50%,
70% и 90% в направлении координат X и Y.
2. Defined Pattern – предопределенный шаблон. При этом надо
дополнительно задать опции Starting Column, Starting Line, Column
Increment (Начальная строка, начальный столбец, шаг приращения
для шаблона). При помощи этой опции вы можете вручную
регулировать расположение областей в пределах которых будут
генерироваться связующие точки.
Основные схемы аэрокосмического
фотографирования
В аэрофотоаппаратах (АФА) реализуются три основные схемы
фотографирования: кадровое, щелевое и панорамное. Применяемые
при этом аэрофотоаппараты называют, соответственно, кадровыми,
щелевыми и панорамными.
При кадровом фотографировании (рис. а) участок земной
поверхности 1 одновременно проецируется при помощи аэрофотообъектива 2 на плоскость 3, с которой совмещается
светочувствительный слой аэрофотоматериала. Фотографическое
изображение местности получается в виде отдельных кадров
аэрофотоснимков или космических снимков, оптическая ось
съемочной камеры перпендикулярна к плоскости 3. Так как при
кадровом фотографировании осуществляется одновременное проектирование участка земной поверхности на весь кадр изображения,
получаемый в этом случае аэрофотоснимок является центральной
31
проекцией местности, т. е. земная поверхность представляется суммой
отдельных изображений, полученных по закону центральной
проекции. Размеры кадра изображения в плоскости 3 ограничены прикладной рамкой.
При щелевом аэрофотографировании (рис. б) полоса земной
поверхности непрерывно последовательно проектируется при помощи
поступательно перемещающегося аэрофотообъектива 2 на непрерывно
движущийся в том же направлении аэрофотоматериал 3 через щель 4.
Скорости перемещения объектива и аэрофотоматериала согласованы.
При щелевом аэрофотографировании изображения земная поверхность
представляется не в виде отдельных кадров, а сплошной лентой —
непрерывной последовательностью изображений узких полос земной
поверхности 7; закону центральной проекции соответствует только
изображение узкой полосы местности 5, получаемое в плоскости щели
4 в некоторый момент времени.
При панорамном аэрофотографировании (рис. в) участок земной
поверхности 1 проектируется при помощи аэрофотообъектива 2 на
аэрофотоматериал 3. При этом обеспечивается, практически
мгновенное, последовательное проектирование участка местности на
аэрофотоматериал,
что
достигается,
например,
вращением
аэрофотообъектива или другими способами.
В результате
панорамного аэрофотографирования получают отдельные полосыпанорамы, изображающие земную поверхность от горизонта до
горизонта.
Выбор геометрической модели съемочной
аппаратуры
Работа по созданию нового блока изображений в модуле
OrthoBASE начинается с выбора категории геометрической модели
датчика съемочной аппаратуры. Категории включают:
- Camera – съемочные камеры различных типов. Как правило,
для этой категории необходимо определить параметры внутреннего
ориентирования;
- Rational Functions – модели, описываемые при помощи
коэффициентов рациональных полиномов;
Polynomial-based Pushbroom – полиномиальные модели
сканирующих разверток. Сбор данных развертки производится вдоль
полосы сканирования. Для каждой полоса сканирования задаются
32
центр проекции и угол вращения. Эти данные задаются при помощи
полиномов;
- Orbital Pushbroom - орбитальные сканирующие развертки –
модели использующие информацию орбитальных эфемерид и
учитывающие кривизну Земли;
- Mixed Sensor датчик смешанного типа. Данная категория не
является отдельной геометрической моделью, а просто использует
геометрическую модель или информацию о модели датчика,
связанного с определенным снимком.
На следующем шаге выбирается непосредственно геометрическая
модель. Опишем лишь один тип геометрических моделей,
предлагаемых в OrthoBASE и используемый в данной курсовой
работе.
Frame Camera – Кадровая камера (категория Camera). Этот тип
камеры ассоциируется с аэрофотокамерами, имеющими фокусное
расстояние 152 мм. Размер кадрирующей рамки – 9x9 дюймов. В
плоскости прикладной рамки нанесены координатные метки, которые
фиксируются на каждом снимке и определяют его координатную
систему. Для различных аэрофотокамер могут задаваться матрицы,
содержащие 4 (2x2), 8 или 289 (17x17) координатных меток. Ниже, на
рисунке приведены примеры некоторых схем расположения
координатных меток на кадрирующей рамке.
Системы координат используемые в OrthoBASE
Файловая система координат.
Файловая система координат, для
цифрового снимка,
определяется оконной (пиксельной) системой координат Windows, то
33
есть, имеет начало в левом верхнем углу окна, содержащего
изображения. Ось x направлена вправо, а ось у - перпендикулярно
вниз. Единицы измерения - пиксели.
Система координат снимка.
Система координат снимка или плоская система координат снимка
определяется как двумерная система координат на плоскости. Начало
системы находится в центре снимка и обычно определяется главной
точкой или точкой пересечения линий проведенных через угловые
(или центральные) координатные метки. Координаты снимка
используются для определения местоположения на плоскости снимка.
Единицами измерения в системе координат снимка являются
миллиметры или микроны. Далее, такая система координат
определяется координатами снимка (х, у).
Соотношение между файловой системой координат и плоской
системой координат снимка представлено на рисунке приведенном
ниже.
Пространственная система координат снимка.
Пространственная система координат снимка
определяется
аналогично плоской системе координат снимка, за исключением того,
34
что добавлена третья ось (z). Центр такой системы координат
определяется центром проекции S. Оси x и у параллельны
соответствующим осям плоской системы координат снимка. Ось z это оптическая ось, следовательно, значение z в любой точке
пространственных координат снимка обычно эквивалентно фокусному
расстоянию f. Пространственные координаты снимка используются
для описания положения снимка внутри камеры. Единицы измерения миллиметры или микроны. Далее, такая система координат
определяется координатами (x, y, z).
Наземная система координат.
Наземная система координат обычно определяется как трехмерная
система координат, которая использует известные картографические
проекции. Наземные координаты (X,Y,Z) обычно выражаются в
метрах. Величина Z определяет превышение над уровнем моря для
данного исходного горизонта. Соотношение между пространственной
системой координат снимка и наземной системой координат
представлено на рисунке приведенном ниже.
35
Элементы внутреннего ориентирования
Внутреннее ориентирование определяет внутреннюю геометрию
камеры или сенсора на момент съемки. Процесс внутреннего
ориентирования заключается в нахождении элементов, определяющих
положение снимка. Внутреннее ориентирование главным образом
используется для преобразования файловой или других систем
координат снимка в систему пространственных координат снимка.
Внутренняя геометрия камеры определяется следующими элементами:
-
координатами главной точки;
фокусным расстоянием;
координатными метками;
дисторсией объектива;
Ниже показано соотношение между
ориентирования и геометрией снимка.
параметрами
внутреннего
Главная точка математически определяется как пересечение
перпендикуляра, проходящего через центр проекции, и плоскости
снимка. Расстояние от главной точки до центра проекции называется
фокусным расстоянием.
36
Плоскость снимка обычно определяется как фокальная плоскость.
Например, для многих широкоформатных аэрофотокамер фокусное
расстояние равно примерно 152 мм. Еще в качестве примера, - для
некоторых цифровых камер оно составляет 28 мм. Основной задачей
перед проведением фотограмметрических работ является точное
определение или калибровка фокусного расстояния метрической
камеры, которое проводится в лабораторных условиях.
Один из этапов внутреннего ориентирования включает в себя
определение положения главной точки для каждого снимка в проекте.
Это положение, как уже указывалось, определяется как точка
пересечения линий, проведенных через угловые (или центральные)
координатные метки.
Таким образом, на снимке измеряются положения координатных
меток, которые затем сравниваются с калибровочными координатами
каждой метки. Так как в данном процессе для каждого снимка система
пространственных координат еще не определена, измеренные
координатные метки связываются с файловой (пиксельной) системой
координат. Файловая система координат имеет координаты х
(столбцы) и у (строки). Ниже показана связь между файловой
системой координат и системой пространственных координат.
Используя двумерное аффинное преобразование, определяют
связь между файловой системой координат и системой
пространственных координат снимка. С помощью двумерного
37
аффинного преобразования могут быть определены коэффициенты,
требующиеся для перевода файловых координат в координаты снимка:
x  a1  a2 * X  a3 * Y
(12)
y  b1  b2 * X  b3 *Y
где x , у - координаты калиброванных координатных меток в
системе координат снимка; X и Y - файловые координаты измеренных
координатных меток. Все эти координаты используются для
определения шести коэффициентов аффинного преобразования.
Полученные коэффициенты могут быть использованы для
преобразования каждой строки (y) и каждого столбца (x) файловых
координат в координаты изображения.
38
Элементы внешнего ориентирования
Внешнее ориентирование определяет положение и угловую
ориентацию снимка. Переменные, характеризующие положение и
ориентацию
снимка,
называются
элементами
внешнего
ориентирования. Элементы внешнего ориентирования определяют
параметры снимка на момент съемки. Рисунок, приведенный ниже,
показывает элементы внешнего ориентирования.
Линейными элементами внешнего ориентирования являются Хо,
Уо и Zo. Они задают положение центра проекции (О) относительно
39
наземной системы координат (X, Y, и Z). Zo обычно определяет
высоту камеры над уровнем моря, которая задается датумом.
Угловые элементы внешнего ориентирования характеризуют
соотношение между наземной системой координат X, Y, и Z) и
системой пространственных координат изображения (х, у и z). Для
определения углового ориентирования обычно используются три угла
поворота омега (v), фи (f), и каппа (k). Угол v определяет поворот
вокруг оси фотографирования х, f определяет поворот вокруг оси
фотографирования у, а k определяет поворот вокруг оси
фотографирования z. Эти повороты являются положительными, когда
они производятся против часовой стрелки, если смотреть из
положительного конца соответствующей оси.
Еще один важный параметр, задаваемый в OrthoBase при работе с
элементами внешнего ориентирования это средняя высота полета.
Средняя высота полета – заданное значение используется как
начальное приближение в процессе автогенерирования связующих
точек. Среднюю высоту полета можно получить тремя способами.
1.
2.
Воспользоваться данными, предоставляемыми поставщиком
снимков, либо взять значение высоты из полетного плана.
Вычислить по формуле:
H = f * m;
где
(13)
f – фокусное расстояние объектива;
m – масштаб снимка;
40
Действительно, согласно рисунку, приведенному выше:
L
H
l 
m*l l
tg   
  H  m* f
l 
H
f
L  m*l

tg 
Например: f = 153.124 мм;
m: 1мм = 40 м;
H =153.124*40 = 6124.96 м;
3.
Получить как результат вычитания средней высоты рельефа
местности из среднего значения Z параметра внешнего
ориентирования
(средней
высоты
центра
проекции
относительно наземной системы координат).
Например:
Z =(Z1+Z2+Z3) / 3 = (1851.948+1853.9934+1853.561) / 3=1853.167;
P = 605м;
H = Z – P = 1248.167 м;
41
Выполнение
ортотрансформирования и проверка
точности стереомодели для снимков
полученных при помощи кадровой
камеры
На рисунке, показанном ниже, представлено геометрическое
соотношение между тремя снимками и их опорными точками, а также
двумя контрольными точками. Опорные точки выбираются в тех
областях, которые легко опознаются на снимке. Исходные снимки не
привязаны в какой-либо системе координат, связанной с опорными
точками. Первоначально, как правило, они привязаны в файловой
системе координат, которая может быть определена как произвольная
координатная сетка.
42
В процессе триангуляции опорные точки используются для
установления математического соотношения между геометрическими
параметрами камеры, фотоснимком и земной поверхностью в трех
измерениях. Полученное математическое соотношение используется
для
проведения
ортотрансформирования.
В
процессе
ортотрансформирования координатная сетка исходного изображения
пересчитывается в картографическую систему координат опорных
точек.
Подготовка файла блока
1. Запустите ERDAS IMAGINE.
2. Закройте Визуализатор, который автоматически открывается
при старте BRDAS IMAGINE.
3. Выберите из главного меню ERDAS IMAGINE пункт Main |
OrthoBASE. Откроется основное окно проекта модуля OrthoBASE.
4. В Основном меню OrthoBASE выберите пункт File | New.
Откроется диалог создания нового файла блока. Когда вы используете
IMAGINE OrthoBASE в вашей работе, вы создаете файлы блоков. Файлы
43
блоков имеют расширение .blk. Один файл может содержать один снимок,
полосу смежных снимков или несколько полос снимков. Файл .blk
является двоичным файлом. В нем содержится информация о блоке, которая
включает в себя информацию о расположении снимков, параметры
камеры, измерения координатных меток, опорных и связующих точек и
тому подобное.
5. Укажите имя создаваемого файла как frame_tour.blk в папке,
где у Вас есть разрешение на запись. Откроется диалог настройки модели
(Model Setup).
6. В
списке
Geometric
Model
Category
(Категория
геометрической модели) выберите опцию Camera. Из списка моделей
камеры (Geometric Model) выберите Кадровую камеру
(Frame
Camera). Нажмите ОК. Откроется диалог настройки свойств блока
(Block Property Setup).
44
7. Нажмите кнопку Set на панели Horizontal Reference
Coordinate System. Откроется диалог выбора проекции Projection
Chooser.
8.
Перейдите на вкладку стандартных проекций Standard и в
списке Categories выберите опцию US State Plane - NAD27 – Old
USGS(DO154)Zone Numbers. Затем, из списка с прокруткой Projection,
выберите COLORADO CENTRAL (3476). Нажмите ОК, для закрытия
диалога выбора проекций.
9. В поле Horizontal Units диалога настройки свойств блока
выберите значение Meters. OrthoBASE использует горизонтальный
датум, заданный Вами на предыдущем шаге, в качестве базового
датума. На его основе, вертикальная координатная система
устанавливается автоматически и в большинстве случаев ее менять не
надо.
10. Нажмите кнопку Next и перейдите в диалог установки
параметров съемки (Frame Specific Information).
45
11. Из всех параметров диалога установки параметров съемки,
измените только поле Average Flying Height (Средняя полетная
высота в метрах). Задайте в этом поле значение 7000. и нажмите ОК.
12. Диалоги настройки будут закрыты и мы снова увидим
основное окно проекта модуля OrthoBASE. Теперь необходимо
добавить в проект набор снимков, составляющих съемочный блок. Для
этого в древовидном списке проекта выделите опцию Images.
13. В меню OrthoBASE выберите пункт Edit | Add Frame
(Добавить снимок). В диалоге выбора файлов переместитесь в
поддиректорию Examples\Orthobase\Frame демонстрационных данных
и при помощи клавиши SHIFT выберите файлы col90p1.img,
col91p1.img, col92p1.img.
Всякий раз, когда Вы
добавляете новые снимки в блок, они
отображаются в виде табличного списка в главном диалоге IMAGINE
OrthoBASE.
К каждому снимку относятся несколько колонок. Колонка номер
строки Row # позволяет вам выбрать конкретный снимок для
использования в IMAGINE OrthoBASE. Например, вы можете захотеть
посчитать пирамидные слои только для одного снимка. Колонка
идентификатор снимка Image ID порядковый номер каждого снимка
блока. Вы можете изменить этот номер, если захотите. Колонка >
позволяет вам определить текущий активный снимок. Колонка названия
снимка Image Name указывает директорию и имя файла каждого снимка.
46
Когда полный путь к снимку точно установлен, соответствующая колонка
Online становится зеленой.
В колонке Active отображается X, обозначая, какой из снимков будет
использоваться для проведения автоматического сбора связующих точек,
триангуляции и ортотрансформирования в IMAGINE OrthoBASE.
Последние пять колонок показывают статус выполнения процесса,
зеленый цвет означает, что обработка выполнена и является точной;
красный цвет обозначает, что обработка не выполнена.
Колонка Руг показывает наличие пирамидных слоев. Колонка Int
указывает - были ли измерены координатные метки. Колонка Ext
показывает - вычислены ли элементы внутреннего ориентирования.
Колонка DTM указывает, создавалась ли цифровая модель местности (в
виде DEM или TIN).
Колонка Ortho указывает - были ли
ортотрансформированы снимки. Колонка Online показывает, что путь к
файлу снимка точно установлен.
14. Затем необходимо построить пирамидные слои для
каждого снимка в блоке. Пирамидные слои предназначены для
оптимизации отображения снимка и определения связующих
точек. Выберите в меню редактирования опцию Edit | Compute
Pyramid Layers. Откроется диалог расчета пирамидных слоев
(Compute Pyramid Layers).
15. Установите в диалоговом окне расчета пирамидных слоев
переключатель в положение All Images Without Pyramids и
нажмите ОК. Этим, Вы запустите процесс для тех снимков, для
которых пирамидные слои еще отсутствуют. Когда процесс будет
завершен, соответствующая колонка Руг. каждого снимка
закрасится в зеленый цвет.
16. Теперь зададим модель съемочной камеры. Когда вы
определяете модель камеры, которой были получены снимки, вы
вводите информацию о положении координатных меток, а также
задаете параметры внешнего ориентирования.
47
17. Выберите Edit | Frame Editor в главном меню OrthoBASE
или щелкните по иконке
свойств снимка (Frame Properties).
Откроется диалог редактирования снимка (Camera Frame Editor),
отображающий информацию о текущем снимке в списке таблицы
IMAGINE OrthoBASE, который выделен знаком >.
18. Сделайте текущим снимок col90p1.img. Далее, в пунктах 19 –
32 для этого снимка мы построим модель камеры, и зададим
параметры внутреннего и внешнего ориентирования. Для задания
названия датчика (поле Sensor Name), нажмите кнопку New Camera.
Откроется диалог ввода информации о камере (Camera
Information).
19. На вкладке General установите следующие параметры:
- в поле Camera Name введите Zeiss RMK A 15/23 в качестве
названия камеры;
- в поле Focal Length введите величину фокуса равную 153.124;
48
- в поле Principal Point X0 введите -0.002 в качестве координаты
X главной точки;
- в поле Principal Point Y0 введите 0.002 в качестве координаты
Y главной точки;
20. Перейдите на вкладку Fiducials. Эта вкладка задает положения
координатных меток.
21. Введите число координатных меток равное 8 в числовое поле
Number of Fiducials, затем нажмите Enter. В таблице появятся
дополнительные координатные метки, которые перечислены колонке
Row#. Введите в таблицу следующую информацию:
Row#
1
Film X (mm)
-103.947
2
103.945
103.924
3
-103.937
103.927
4
103.958
-103.952
5
-112.996
-0.005
6
112.990
-0.015
7
0.003
113.001
8
0.026
-112.971
49
Film Y (mm)
-103.952
Примечание: в данном задании не определяется радиальная дисторсия
объектива. Поэтому, параметры вкладки Radial Lens Distorsion не
устанавливаем.
22. Нажмите ОК для закрытия диалога Camera Information.
Вернитесь в диалог редактирования снимка (Camera Frame Editor).
Перейдите на вкладку Interior Orientation. Отобразятся инструменты
для задания параметров внутреннего ориентирования снимка.
23. Убедитесь, что установлен режим ориентирования меток,
соответствующий кнопке
(панель Fiducial Orientation). Нажав
кнопку вызова Визуализатора,
перейдите в режим просмотра.
Поверх диалога Frame Editor откроется главное окно (Main View), обзорное
окно (Over View), которое показывает входной снимок и окно детального
просмотра (Detail View), показывающее часть снимка, на которой
находится связанный курсор. Любой из трех Визуализаторов может быть
использован для измерения координатных меток, однако обычно легче
сажать метки в окне детального просмотра.
50
24.Примерная область расположения первой координатной метки
определяется связанным курсором в главном окне и отображается в окне
детального просмотра.
25. Нажмите в центр связанного курсора в главном окне и
перетащите его так, чтобы крестообразная координатная метка
оказалась в центре. Теперь координатная метка окажется в центре окна
детального просмотра. Щелкните по иконке размещения метки
(Place Image Fiducial) на панели Interior Orientation для того, чтобы
поместить изображение координатной метки. На рисунке, приведенном
51
ниже, показана форма связанного курсора и координатная метка,
посаженная точно в центр перекрестия.
26. Для измерения первой координатной метки нажмите на ее
изображение в центре окна детального просмотра. Координатная метка
измерена, и ее координаты переданы в соответствующие колонки Image X
и Image Y. Если нажата кнопка включения режима автоматического
перемещения
, изображение автоматически перемещается в область
приблизительного расположения следующей координатной метки.
27. Для посадки остальных меток повторите шаги 25 и 26.
Примечание: Координатные метки с 5 no 8 будут отображаться на
фотоснимке в виде точек, a не в виде перекрестий. Они являются
дополнительными координатными метками.
28. Когда все восемь координатных меток будут измерены,
связанный курсор вернется к первой координатной метке.
Среднеквадратическая ошибка определения положения координатных
меток (выводится кнопкой Solve в панели Interior Orientation диалога
редактирования снимка) обычно должна быть меньше 0 33 пикселей. Если
ошибка больше. 0.3J пикселей, повторите измерения некоторых
координатных меток. Зачастую первые несколько измерений не
так точны, как последующие. Для проведения повторного
измерения, выделите строку с номером точки Point #, которую вы
хотите перемерить, нажмите инструментальную кнопку выделения
52
объектов
, затем переместите координатную метку в окне
детального просмотра.
29. Отожмите кнопку вызова Визуализатора
, чтобы выйти
из режима просмотра.
30. Для задания элементов внешнего ориентирования
перейдите на закладку Exterior Information в диалоге Frame
Editor.
31. В числовое поле Value введите следующую информацию.
Эти значения соответствуют снимку со190р1.
Xo
666700.000
Yo
115900.000
Zo
8800.000
Omega
0.0000
Phi
Kappa
0.0000 90.0000
32. Установите
переключатель
Set
Value
и
в
соответствующем списке, выберите для всех параметров внешнего
ориентирования тип Initial.
33. Построение модели для снимка col90p1.img закончено.
Щелкните в диалоге Frame Editor кнопку Next для перехода к
следующему снимку - col91p1.img.
53
34. Повторите шаги 22 – 32 для изображения col91p1.img. При
выполнении 31 шага, введите следующую информацию:
Xo
666700.000
Yo
119400.000
Zo
8800.000
Omega
0.0000
Phi
0.0000
Kappa
90.0000
35. Щелкните в диалоге Frame Editor кнопку Next для
перехода к следующему снимку - col92p1.img. Повторите шаги 22
– 32 для изображения col92p1.img. При выполнении 31 шага, введите
следующую информацию:
Xo
666800.000
Yo
122900.000
Zo
Omega
8800.000 0.0000
Phi
0.0000
Kappa
90.0000
36. Ввод данных для всех трех снимков закончен. Нажмите ОК в
диалоге Frame Editor. В главном меню OrthoBASE выберите пункт
File | Save, для сохранения промежуточных результатов построения
модели.
Примечание: другим вариантом представленного здесь процесса является
обработка каждого снимка целиком, (то есть, задание всех сенсоров,
измерение всех координатных меток, затем, ввод всех элементов
внешнего ориентировании).
Ввод опорных и контрольных точек
37. Теперь, когда вы измерили координатные метки и задали
элементы внешнего ориентирования для каждого снимка, составляющего
блок, вы готовы к использованию инструмента для измерения точек (Point
Measurement Tool) для определения положения наземных опорных
(Ground Control Point - GCP), связующих (Tie) и контрольных (Check)
точек по снимку. При работе с IMAGINE OrthoBASE наряду с
другими
источниками
информации,
предназначенными для
определения расположения опорных точек, вы можете использовать
обзорный фотоснимок. В этом упражнении моделируется именно такой
тип работы, использующий снимок, набросок и детальный обзор для
каждой опорной точки.
54
Обычно опорные точки располагаются в местах пересечения дорог,
на углах зданий или межевых знаках. Вы должны избегать размещения
опорных точек в местах, которые подвержены изменениям, таких как
граница леса или береговая линия. Вы можете использовать
топографическую карту масштаба 1:24000, результаты государственной
съемки или полевые наброски, которые помогут вам определить
местоположение соответствующих точек.
В основном, процесс сбора опорных точек состоит из
изучения топографических карт для определения положения
характерных точек, таких как пересечение дорог. Определив их
местоположение, вы используете IMAGINE OrthoBASE для записи
этой информации и ввода опорных координат. Затем вы проверяете по
обзорной фотографии, где расположены точки. Сбор опорных точек
может быть трудоемким процессом, однако от его точности зависит
точность триангуляции и ортотрансформирования снимков в вашем
блоке.
Когда вы определите все опорные точки, общие для двух или более
снимков в блоке, вы можете приступить к выполнению триангуляции.
Минимальное количество опорных точек на снимок - три
вертикальные и две горизонтальные. Если вы используете функцию
автоопределения связующих точек, вам понадобится меньшее
количество опорных точек.
Триангуляция связывает снимки друг с другом, таким образом они
могут быть ортотрансформированы вместе. Точно заданные опорные
точки делают проведение триангуляции и ортотрансформирования
успешным.
38. Выберите в главном меню IMAGINE OrthoBASE пункт Edit |
Point Measurement. Откроется диалог измерения точек (Point Measurement),
состоящий из набора инструментов, таблиц с опорными и файловыми
координатами и шести Визуализаторов, составляющих две группы,
предназначенные для отображения пары снимков в IMAGINE
OrthoBASE. В данном случае будут отображены первые два файла блока
со190р1 и со191р1. Задавать пары изображений возможно, выбирая
названия правого и левого снимков из списков Left View и Right View.
55
ВНИМАНИЕ! В данном упражнении приведен один из возможных
вариантов набора опорных и контрольных точек. При выполнении
курсовой работы, выберите вариант соответствующий Вашему
заданию. Для всех вариантов нумерация точек сохраняется. В
качестве первоначального приближения при нахождении
местоположения точек Вы можете использовать справочные
схемы данного варианта (приводятся ниже). Для окончательного
выбора местоположения точек используйте схемы размещения
точек, приведенные для каждого варианта.
39. Измерение точки ID 1002. Нажмите на кнопку Add в правом
верхнем углу панели инструментов Point Measurement Tool для ввода
первой точки. Это добавит новую строку в таблицу координат в нижней
части панели инструментов.
40. Позиционируйтесь мышью на ячейку Point ID и исправьте
идентификатор на ID 1002.
56
41. Позиционируйтесь мышью на ячейку Type и выберите в
соответствующем списке тип Full (этот тип указывает, что заданы все
три координаты опорной точки X, Y и Z).
42. Позиционируйтесь мышью на ячейку Usage и выберите в
соответствующем списке тип Control (этот тип указывает, что точка
будет использоваться в качестве опорной).
43. Следующая справочная схема покажет вам, где расположены
точки на первом снимке блока со190р1. Вы будете обращаться к этому
рисунку при определении точек. Местоположение точек выделено
квадратами и помечено соответствующим образом.
44. Вы будете использовать опорный снимок местности, а
также набросок и детальный вид опорной точки 1002 для того,
чтобы определить ее местоположение
57
45. Следующий рисунок поможет Вам определить примерное
местоположение точки. Также помещен детальный вид точки ID
1002.
46. На основании вспомогательных данных найдем положение
точки ID 1002 на левом снимке пары (со190р1). Для этого
используем три левых Визуализатора. В главном окне перемещайте
связанный курсор до тех пор, пока опорная точка не будет видна
в окне детального просмотра. Измените размер связанного курсора, если
необходимо
изменить масштаб в окне детального просмотра. Нажмите на иконку
создания точки
(Create Point) на панели инструментов и посадите
точку на левый снимок. Файловые координаты на левом снимке
созданной точки появятся в связующей таблице.
47. Аналогичным образом найдем положение точки ID 1002 на
правом снимке пары (со191р1). Для этого используем три правых
Визуализатора. Файловые координаты на правом снимке посаженой
точки появятся в связующей таблице. Для контроля можете использовать
следующие значения файловых координат на левом и правом снимках:
Image
col90p1
col91p1
X File
952.625
165.875
Y File
819.625
846.625
48. Введем опорные координаты для точки ID 1002:
- нажмите на колонку X Reference, соответствующую точке ID 1002;
- напечатайте 665228.955 в колонке X Reference;
- нажмите на колонку Y Reference, соответствующую точке ID 1002;
- напечатайте 115012.472 в колонке Y Reference;
58
- нажмите на колонку Z Reference, соответствующую точке ID 1002;
- напечатайте 1947.672 в колонке Z Reference, затем нажмите Enter;
- когда вы закончите позиционирование опорной точки ID 1002,
нажмите Save на панели инструментов;
49. Измерение точки ID 1003. Нажмите на кнопку Add в правом
верхнем углу панели инструментов Point Measurement Tool для ввода
точки. Следующий рисунок покажет вам, где расположены точки на
втором снимке блока со191р1.
50. Позиционируйтесь мышью на ячейку Point ID и исправьте
идентификатор на ID 1003. Позиционируйтесь мышью на ячейку Type
59
и выберите в соответствующем списке тип Full. Позиционируйтесь
мышью на ячейку Usage и выберите в соответствующем списке тип
Control.
51. Следующий рисунок поможет Вам определить примерное
местоположение точки. Также помещен детальный вид точки ID
1003.
52. На основании вспомогательных данных найдем положение
точки ID 1003 на левом снимке пары (со190р1). Для этого используем
три левых Визуализатора. В главном окне перемещайте связанный
курсор до тех пор, пока опорная точка не будет видна
в окне детального просмотра. Измените размер связанного курсора,
если необходимо изменить масштаб в окне детального просмотра.
Нажмите на иконку создания точки Create Point
на панели
инструментов и посадите точку на левый снимок. Файловые
координаты на левом снимке посаженой точки появятся в связующей
таблице.
53. Аналогичным образом найдем положение точки ID 1003 на
правом снимке пары (со191р1). Для этого используем три правых
Визуализатора. Файловые координаты созданной на правом снимке
точки, появятся в связующей таблице.
54. Опорная точка ID 1003 может быть позиционирована на всех
трех снимках, поэтому сменим пару снимков. В соответствующем
списке Right View выберите файл правого снимка со192р1. Найдем
положение точки ID 1003 на правом снимке пары (со192р1) и посадим
на нем данную опорную точку. Следующий рисунок покажет вам, где
расположены точки на третьем снимке блока со192р1.
60
Для контроля можете использовать следующие значения файловых
координат на трех снимках:
Image
col90p1
col91p1
col92p1
X File
1857.875
1064.875
286.875
Y File
639.125
646.375
639.125
55. Введем опорные координаты для точки ID 1003:
- нажмите на колонку X Reference, соответствующую точке ID 1003;
- напечатайте 664456.22 в колонке X Reference;
- нажмите на колонку Y Reference, соответствующую точке ID 1003;
- напечатайте 119052.15 в колонке Y Reference;
- нажмите на колонку Z Reference, соответствующую точке ID 1003;
61
- напечатайте 1988.820 в колонке Z Reference, затем нажмите Enter;
- когда вы закончите позиционирование опорной точки ID 1003,
нажмите Save на панели инструментов;
Местоположение точки ID 1003 определено на всех трех снимках.
56. IMAGINE OrthoBASE обеспечивает накоторые автоматические
функции для более быстрого выделения опорных точек. Одной из таких
функций является автоматическое слежение (х, у). Нажмите иконку
автоматического слежения (х, у)
.
Используя функцию
автоматического слежения (х, у), IMAGINE OrthoBASE
аппроксимирует положение точки на втором снимке, основываясь
на положении точки на первом снимке.
57. Измерение точки ID 1004. Точка ID 1004 расположена на
снимках со191р1 и со192р1. Последняя точка, которую вы
определили, находилась на снимке со192р1, который отображен в
правом Визуализаторе. Прежде чем начать, вы должны
настроить Left View на отображение снимка со191р1. Нажмите на
спадающее меню Left View и выберите со191 р1. Снимок со191р1
отобразится в трех левых Визуализаторах; снимок со192р1 уже
отображен в трех правых Визуализаторах.
58. Нажмите на кнопку Add в правом верхнем углу панели
инструментов Point Measurement Tool для ввода точки ID 1004.
Позиционируйтесь мышью на ячейку Point ID и исправьте
идентификатор на ID 1004. Позиционируйтесь мышью на ячейку Type
и выберите в соответствующем списке тип Full. Позиционируйтесь
мышью на ячейку Usage и выберите в соответствующем списке тип
Control.
59. Следующий рисунок поможет Вам определить примерное
местоположение точки. Также помещен детальный вид точки ID
1004.
62
60. Найдем положение точки ID 1004 на левом (со191р1) и
правом (со192р1) снимках пары. Обратите внимание на то, что так
как Вы включили функцию автоматического слежения (х, у),
IMAGINE OrthoBASE установит снимки в Визуализаторах так, что
отобразится примерное расположение точки на втором снимке
со192р1. Для контроля можете использовать следующие значения
файловых координат на левом и правом снимках:
Image
col91p1
col92p1
X File
1839.52
1050.60
Y File
1457.43
1465.23
61. Введем опорные координаты для точки ID 1004:
- нажмите на колонку X Reference, соответствующую точке ID 1004;
- напечатайте 668150.61 в колонке X Reference;
- нажмите на колонку Y Reference, соответствующую точке ID 1004;
- напечатайте 122404.68 в колонке Y Reference;
- нажмите на колонку Z Reference, соответствующую точке ID 1004;
- напечатайте 1972.056 в колонке Z Reference, затем нажмите Enter;
- когда вы закончите позиционирование опорной точки ID 1004,
нажмите Save на панели инструментов;
62. Следующая опорная точка ID 1005 расположена на
всех трех снимках блока. Обычно такой случай считается
наилучшим для выделения опорной точки на снимках блока (то есть,
определение точки на снимке со190р1, затем на со191р1, затем на со192р1).
Вы должны использовать ниспадающее меню Left View и Right View для
того, чтобы вернуть порядок расположения снимков. Нажмите на
ниспадающее меню Left View и выберите со190р1. Нажмите на
ниспадающее меню Right View и выберите со191р1.
63. Нажмите на кнопку Add в правом верхнем углу панели
инструментов Point Measurement Tool для ввода точки ID 1005.
Позиционируйтесь мышью на ячейку Point ID и исправьте
идентификатор на ID 1005. Позиционируйтесь мышью на ячейку Type
и выберите в соответствующем списке тип Full. Позиционируйтесь
мышью на ячейку Usage и выберите в соответствующем списке тип
Control.
63
64. Следующий рисунок поможет Вам определить примерное
местоположение точки. Также помещен детальный вид точки ID 1005.
65. Найдем положение точки ID 1005 на левом (со190р1) и правом
(со191р1) снимках пары. Так как Вы включили функцию
автоматического слежения (х, у), IMAGINE OrthoBASE установит
снимки в Визуализаторах так, что отобразится примерное
расположение точки на втором снимке со192р1. Используйте
соответствующий список Right View для того, чтобы определить точку
ID 1005 на снимке col92p1. Для контроля можете использовать
следующие значения файловых координат на левом и правом снимках:
Image
col90p1
col91p1
col92p1
X File
1769.450
1007.250
224.670
Y File
1508.430
1518.170
1510.670
66. Введем опорные координаты для точки ID 1005:
- нажмите на колонку X Reference, соответствующую точке ID 1005;
- напечатайте 668338.22 в колонке X Reference;
- нажмите на колонку Y Reference, соответствующую точке ID 1005;
- напечатайте 118685.9 в колонке Y Reference;
- нажмите на колонку Z Reference, соответствующую точке ID 1005;
- напечатайте 1886.712 в колонке Z Reference, затем нажмите Enter;
- когда вы закончите позиционирование опорной точки ID 1005,
нажмите Save на панели инструментов;
67. Следующая опорная точка ID 1006 расположена на
всех трех снимках блока. Вы должны использовать ниспадающее
64
меню Left View и Right View для того, чтобы вернуть порядок
расположения снимков. Нажмите на ниспадающее меню Right View
и выберите со191р1, затем на ниспадающее меню Left View и
выберите со190р1.
68. Нажмите на кнопку Add в правом верхнем углу панели
инструментов Point Measurement Tool для ввода точки ID 1006.
Позиционируйтесь мышью на ячейку Point ID и исправьте
идентификатор на ID 1006. Позиционируйтесь мышью на ячейку Type
и выберите в соответствующем списке тип Full. Позиционируйтесь
мышью на ячейку Usage и выберите в соответствующем списке тип
Control.
69. Следующий рисунок поможет Вам определить примерное
местоположение точки. Также помещен детальный вид точки ID
1006.
70. Найдем положение точки ID 1006 на левом (со190р1) и
правом (со191р1) снимках пары. Так как Вы включили функцию
автоматического слежения (х, у), IMAGINE OrthoBASE установит
снимки в Визуализаторах так, что отобразится примерное
расположение точки на втором снимке со192р1. Используйте
соответствующий список Right View для того, чтобы определить
точку ID 1005 на снимке col92p1. Для контроля можете
использовать следующие значения файловых координат на левом и
правом снимках:
Image
col90p1
col91p1
col92p1
X File
1787.875
1023.625
215.125
Y File
2079.625
2091.390
2083.790
71. Введем опорные координаты для точки ID 1006:
65
- нажмите на колонку X Reference, соответствующую точке ID 1006;
- напечатайте 670841.48 в колонке X Reference;
- нажмите на колонку Y Reference, соответствующую точке ID 1006;
- напечатайте 118696.89 в колонке Y Reference;
- нажмите на колонку Z Reference, соответствующую точке ID 1006;
- напечатайте 2014.0 в колонке Z Reference, затем нажмите Enter;
- когда вы закончите позиционирование опорной точки ID 1006,
нажмите Save на панели инструментов;
Теперь вы введете информацию о двух контрольных точках в
таблицы опорных и файловых координат. Информация о контрольных
точках вводится так же, как и в случае с опорными точками,
единственным различием является указание типа точки Check
(контрольная) в колонке Usage. Контрольные точки дополняют опорные
и используются для определения точности триангуляции. Контрольные
точки фактически не являются необходимыми условием для выполнения
триангуляции.
72. Ввод контрольной точки ID 2001. Нажмите на ниспадающее
меню Right View и выберите со191р1. Нажмите на кнопку Add в
правом верхнем углу панели инструментов Point Measurement Tool для
ввода точки ID 2001. Позиционируйтесь мышью на ячейку Point ID и
исправьте идентификатор на ID 2001. Позиционируйтесь мышью на
ячейку Type и выберите в соответствующем списке тип Full.
Позиционируйтесь мышью на ячейку Usage и выберите в
соответствующем списке тип Check.
73. Следующий рисунок поможет Вам определить примерное
местоположение точки. Также помещен детальный вид точки ID
2001.
66
74. На основании вспомогательных данных найдем положение
точки ID 2001 на левом снимке пары (со190р1). Для контроля можете
использовать следующие значения файловых координат на левом и
правом снимках:
Image
col90p1
col91p1
X File
915.02
160.90
Y File
2095.71
2127.84
- нажмите на колонку X Reference, соответствующую точке ID 2001;
- напечатайте 670970.45 в колонке X Reference;
- нажмите на колонку Y Reference, соответствующую точке ID 2001;
- напечатайте 114815.23 в колонке Y Reference;
- нажмите на колонку Z Reference, соответствующую точке ID 2001;
- напечатайте 1891.888 в колонке Z Reference, затем нажмите Enter;
- когда вы закончите позиционирование опорной точки ID 2001,
нажмите Save на панели инструментов;
75. Ввод контрольной точки ID 2002. Нажмите на ниспадающее
меню Right View и выберите со192р1, а затем на ниспадающее
меню Left View и выберите со191р1. Нажмите на кнопку Add в
правом верхнем углу панели инструментов Point Measurement Tool для
ввода точки ID 2002. Позиционируйтесь мышью на ячейку Point ID и
исправьте идентификатор на ID 2002. Позиционируйтесь мышью на
ячейку Type и выберите в соответствующем списке тип Full.
Позиционируйтесь мышью на ячейку Usage и выберите в
соответствующем списке тип Check.
76. Следующий рисунок поможет Вам определить примерное
местоположение точки. Также помещен детальный вид точки ID
2002.
77. На основании вспомогательных данных найдем положение
точки ID 2002 на левом снимке пары (со191р1). Для контроля можете
67
использовать следующие значения файловых координат на левом и
правом снимках:
Image
X File
Y File
col91p1
2032.030
2186.530
col92p1
1227.375
2199.125
- нажмите на колонку X Reference, соответствующую точке ID 2002;
- напечатайте 671408.73 в колонке X Reference;
- нажмите на колонку Y Reference, соответствующую точке ID 2002;
- напечатайте 123166.52 в колонке Y Reference;
- нажмите на колонку Z Reference, соответствующую точке ID 2002;
- напечатайте 1983.762 в колонке Z Reference, затем нажмите Enter;
- когда вы закончите позиционирование опорной точки ID 2002,
нажмите Save на панели инструментов;
По окончанию ввода точек, диалог измерения точек будет выглядеть
следующим образом:
68
Выполнение автоматического сбора
связующих точек
В процессе автоматического сбора связующих точек измеряются
положения опорных точек в координатах снимка, представленных на двух
или более снимках с перекрытием.
78. На панели инструментов измерения точек, нажмите
иконку
Automatic Tie Point Collection Properties, которая
сделает доступным диалог свойств автоматического сбора
связующих точек.
79. На вкладке General оставьте все параметры без изменений.
Перейдите на вкладку Distribution.
80. В поле Intended Number of Points Per Image задайте примерное
количество точек на снимок равное 15. Сбросьте переключатель Keep
All Points. Нажмите RUN в диалоге свойств автоматического создания
связующих точек, чтобы запустить процесс. Фактическое количество
точек, созданных автоматически, будет больше или меньше 15. Это
количество частично зависит от величины перекрытия между
снимками. В данном случае вы должны получить 24 сгенерированных
связующих точки.
ВНИМАНИЕ! В данном упражнении приведен один из возможных
вариантов задания количества автоматически генерируемых
связующих точек. При выполнении курсовой работы, выберите
вариант соответствующий Вашему заданию.
81. Программа
IMAGINE
OrthoBASE
запустит
процесс
автоматического создания связующих точек, выполнение этого процесса
будет отображено на индикаторе состояния и покажет, когда процесс будет
69
завершен. Тип связующих точек указан в колонке Туре как None, а
применение - в колонке Usage как Tie. Точки добавляются в таблицу
опорных координат с соответствующими значениями X File и Y File в
таблице файловых координат.
Вы должны проверить несколько связующих точек, чтобы убедиться
в их точности. Если точка не обладает желательной точностью, вы можете
внести поправки с помощью инструментов выделения или удаления точки,
выделяя строку в таблице или нажимая на кнопку Delete.
82. Для проверки точности связующих точек, в таблице опорных
координат нажмите на полосу прокрутки и выберите какую-нибудь точку,
например ID 2022. Нажмите на колонку > чтобы посмотреть точку ID
2022. Точка отобразится в Визуализаторах. Можно увидеть, что она
является вполне точной.
83. Проверьте несколько других точек. При необходимости,
выполните коррекцию.
84. Нажмите кнопку Close для того, чтобы закрыть диалог измерения
точек. Вы вернетесь в главный диалог IMAGINE OrthoBASE.
Выполнение аэротриангуляции
Теперь, когда найдены все опорные, контрольные и связующие точки,
IMAGINE OrthoBASE имеет всю необходимую информацию для
проведения аэротриангуяции. Триангуляция является этапом процесса
установления математического соотношения между снимками,
составляющими блок.
85. Из меню Edit главного диалога IMAGINE OrthoBASE выберите
опцию свойств триангуляции (Triangulation Properties). Откроется
диалог аэротриангуляции.
70
86. Перейдите на вкладку Point. В поле Type выберите из
соответствующего списка значение Same Weighted Values. Нажмите
кнопку Run для запуска аэротриангуляции. Создастся и откроется диалог
отчета триангуляции.
87. Вы можете захотеть обратиться к отчету для более
детального знакомства с результатами триангуляции. Вы можете
сохранить отчет в виде текстового файла для того, чтобы
обращаться к нему в будущем. В диалоге отчета триангуляции
(Triangulaiion Summary), нажмите кнопку Report. Откроется отчет
триангуляции.
71
88. Прокрутитe отчет вниз до тex пор, пока не увидите раздел
отчета The Output of Self-calibrating Bundle Block Adjustment, в
котором представлены выходные данные самокалибрующегося
уравнивания блока по методу связок.
Обратите внимание на стандартную ошибку (standard error) итерации
номер 3. В отчете стандартная ошибка равна 0.2568. Это стандартное
отклонение единицы веса, которое показывает общее качество этой
итерации 3. Обратите внимание на строчку элементов внешнего
ориентирования exterior orientation parameters. Здесь представлены
элементы внешнего ориентирования для каждого снимка в блоке.
Обратите внимание на невязки опорных точек (residuals of the control
point) и невязки контрольных точек (residuals of the check points).
Невязки опорных точек X, Y и 2 определяют точность качества и как
среднюю невязку и как невязку геометрического среднего. Невязки X, Y, и
Z контрольных точек служат в качестве независимой проверки при
определении среднего качества решения.
Снова прокрутите вниз до тех пор, пока не увидите раздел невязок
точек снимка residuals of image point. Это невязки координат снимка или
фотографии. Их точность может быть гораздо меньше. Как правило, эти
точки имеют большую невязку.
89. Из меню File выберите Save As (сохранить как). Откроется
диалог сохранения отчета. Выберите директорию, в которую хотите
записать файл. В текстовом
поле File name пиелите имя файла
frame_report.txt, затем нажмите Enter.
90. В диалоге отчета триангуляции нажмите кнопку Update для того,
чтобы обновить значения элементов внешнего ориентирования. Эта
операция удалит значения элементов внешнего ориентирования, которые
вы ввели в процессе измерения координатных меток, и запишет вместо
них значения, вычисленные IMAGINE OrthoBASE на основе опорных и
связующих точек в снимках, составляющих файл блока.
Нажмите кнопку Close для того, чтобы закрыть диалог отчета
триангуляции.
В диалоге аэротриангуляции (Aerial Triangulation) нажмите кнопку
Accept для того, чтобы принять параметры триангуляции.
Нажмите ОК в диалоге аэротриангуляции, чтобы закрыть его.
91. Главный диалог IMAGINE OrthoBASE отражает выполнение этапа по
установлению элементов внешнего ориентирования.
72
Проверка точности стереомодели
Выполнив триангуляцию и сохранив файл блока, Вы имеете
возможность оценить точность созданной модели в стереорежиме. Для
этого открыть blk-файл в модуле Stereo Analyst. Инструментарий
позиционирования Stereo Analyst используется для того, чтобы задать
трехмерные координаты контрольных точек. Эти координаты, затем
следует сравнить с ожидаемым положением контрольных точек в
стереовизуализаторе. Если положение контрольной точки правильное,
плавающий курсор должен находиться в той же позиции.
92. Сохраните файл проекта выбрав из меню File | Save.
93. Из главного меню ERDAS выберите Main | Stereo Analyst, чтобы
начать работу в Стерео-Аналитике.
94. Из главного меню Stereo Analyst выберите пункт File | Open |
Imagine OrthoBase Block File. Укажите в качестве имени blk-файла Ваш
файл проекта, в котором уже выполнена триангуляция и нажмите OK. Ваш
проект откроется в стереовизуализаторе.
73
95. Выберите стереопару col90p1_col91p1, на которой задается
первая контрольная точка ID 2001. Для этого нажмите инструментальную
кнопку выбора стереопары
(Choose Another Stereo Pair) в панели Stereo
Analyst. Откроется диалоговое окно, позволяющее задать снимки,
включаемые в состав стереопары.
В открывшемся диалоге выберите первую пару из предлагаемого
списка, нажмите Apply, а затем Close.
74
96. В панели кнопок быстрого доступа Stereo Analyst выберите
инструментальную кнопку ручного ввода координат
(Show Tool For
Manualy Ajusting Image Position Information). Откроется диалоговое окно
ручного задания координат.
97. Обязательно убедитесь, что переключатель Enable Update
сброшен. В поле Zoom введите масштаб отображения 4.0. Установите
координаты первой контрольной точки (ID 2001) в диалоговом окне
ручного задания координат. В полях, задающих координаты, введите
следующие исходные значения:
X: 670970.45
Y: 114815.23
Z: 1891.888
75
ВНИМАНИЕ! В данном упражнении приведен один из возможных
вариантов задания координат контрольных точек. При
выполнении данного этапа курсовой работы, выберите вариант
расположения контрольных точек и справочных рисунков,
соответствующий Вашему заданию.
Большое перекрестие курсора
должно указывать на область,
соответствующую положению первой контрольной точки. Для справки
можно воспользоваться рисунком, приведенным ниже.
98. Если положение перекрестия курсора визуально соответствует
координатам точки ID 2001, перейдите к пункту 104. В противном
случае, выполните следующие действия.
99. Проверим точность стереомодели по направлениям X и Y.
Запишите исходные координаты X и Y. Установите переключатель
Enable Update. После этого, перетяните изображения стереопары (нажав и
удерживая левую клавишу мыши) так, чтобы положение контрольной
точки ID 2001 совпало с перекрестием курсора.
100. Снова сбросьте переключатель Enable Update. Запишите новые
координаты X и Y. Разница между соответствующими исходными и
новыми координатами позволяет оценить точность стереомодели в
направлении X (Y) для текущей контрольной точки.
101. Проверим точность стереомодели по высоте Z. Поместите
плавающий курсор в предполагаемое расположение контрольной точки
ID 2001, сверяясь со справочной схемой приведенной выше. Запишите
исходную высоту Z.
102. На правом и левом визуализаторах в режиме моно плавающий
курсор должен находиться над одним и тем же объектом. Если это не так,
требуется регулирование продольного параллакса. Для этого нажмите
клавишу X на клавиатуре компьютера, затем левую клавишу мыши и,
удерживая их, перемещайте мышь в продольном направлении.
76
77
Работа c Визуализатором модуля Stereo Analyst
103. Когда на правом и левом моновизуализаторах положение
курсора станет идентичным, запишите новое значение Z. Разница между
соответствующими исходными и новыми координатами Z позволяет
оценить точность стереомодели по высоте для текущей контрольной
точки.
104. Выберите стереопару col91p1_col92p1, на которой задается
вторая контрольная точка ID 2002. Для этого нажмите инструментальную
кнопку выбора стереопары
(Choose Another Stereo Pair) в панели Stereo
Analyst. В открывшемся диалоге выберите первую пару из предлагаемого
списка, нажмите Apply, а затем Close.
105. Убедитесь, что переключатель Enable Update сброшен. В поле
Zoom введите масштаб отображения 4.0.
Установите координаты второй контрольной точки (ID 2002) в диалоговом
окне ручного задания координат. В полях, задающих координаты, введите
следующие исходные значения:
X: 671408.73
Y: 123166.52
Z: 1983.762
ВНИМАНИЕ! Выберите вариант расположения контрольных
точек и справочных рисунков, соответствующий Вашему заданию.
Большое перекрестие курсора
должно указывать на область,
соответствующую положению первой контрольной точки. Для справки
можно воспользоваться рисунком приведенным ниже.
106. Если положение перекрестия курсора визуально соответствует
координатам точки ID 2002, перейдите к пункту 107. В противном
случае, выполните пункты 99 -103 для второй контрольной точки.
По окончании проверки точности стереомодели закройте Stereo Analyst.
78
Ортотрансформирование снимков
Следующим шагом является создание ортотрансформированных
снимков из Вашего блока. Ортотрансформированные снимки имеют
меньшее смещение за рельеф и меньшие геометрические ошибки, чем
неортотрансформированные снимки, что делает их более точными.
Объекты, отображенные на ортотрансформированных снимках, имеют
точные координаты X, Y, и Z.
107. В главном диалоге IMAGINE OrthoBASE, нажмите иконку
ортотрансформирования
(Ortho
Resampling).
Откроется
диалог
ортотрансформирования. Перейдите на вкладку General.
108. В соответствующем списке DTM Source выберите опцию
DEM. DEM-файл будет использоваться для обеспечения информации
о высотах для ортотрансформируемых снимков. В поле Vertical Units
установите единицы измерений Meters. В соответствующем списке
DEM File Name выберите опцию Find DEM и откройте файл
colspr_dem.img из подкаталога /examples/orthobase/frame как DEMфайл.
79
109.
Перейдите на вкладку Advanced. В поле Resample
Method соответствующего списка, выберите Nearest Neighbor.
Установите переключатель Ignore Value, и игнорируемое значение 0. По умолчанию, OrthoBASE предполагает, что Вы хотите
генерировать ортоизображение только для первого снимка в блоке.
Используя кнопку Add Multiple Вы можете добавить все остальные
снимки в Ваш блок ортотрансформирования.. Нажмите кнопку Add
Multiple. Откроется диалог Add Multiple Outputs.
110. Убедитесь, что префикс ortho находится в поле File name.
Файлы, содержащие ортотрансформированные снимки, будут иметь в
своем названии этот префикс. Установите переключатель Use Current
Cell Sizes и нажмите ОК.
111. Откроется диалог состояния, показывающий выполнение
процесса ортотрансформирования. Когда диалог состояния покажет,
что
выполнено
100%.
нажмите
кнопку
ОК,
чтобы
закрыть его. Сохраните проект съемочного блока.
112. Теперь, когда снимки были ортотрансформированы, вы
можете просмотреть их совместно в Визуализаторе. Для этого в
ERDAS IMAGINE запустите Viewer.
113. Выберите директорию, в которой вы сохранили
ортофансформпроманные сннмки. Выберите при помощи клавиши
SHIFT
снимки
orthocol90p1.img,
orthocol91p1.img,
orthocol92p1.img. На закладке Raster Options сбросьте переключатель
Clear Display и установите переключатели Fit to Frame и Background
Transparent. Нажмпте ОК.
80
114. Используйте утилиты Swipe и Flicker, для просмотра областей
перекрытия снимков.
81
Список литературы:
1. Erdas Field Guide – GIS & Mapping, LLC Atlanta, Georgia, Seventh
Edition, 2003 – 698 p.
2. Erdas Imagine Tour Guides – GIS & Mapping, LLC Atlanta, Georgia,
2003 – 670 p.
3. Leica Photogrammetry Suite OrthoBASE & OrthoBASE Pro User’s
Guide – GIS & Mapping, LLC Atlanta, Georgia, 2003 – 489 p.
4. Leica Photogrammetry Suite Stereo Analyst User’s Guide – GIS &
Mapping, LLC Atlanta, Georgia, 2003 – 267 p.
5. Сборник упражнений по работе с ERDAS IMAGINE. – Москва:
Data+, 2004. -108 стр.
6. Савиных В.П., Кучко А.С., Стеценко А.Ф. Аэрокосмическая
фотосъемка. – Москва: Картогеоцентр – Геоиздат, 1997. -378 стр.
7. Матиясевич Л.Н. Введение в космическую фотографию. – Москва:
Недра, 1989.-148стр.
8. Лобанов А.Н., Буров М.И., Краснопевцев Б.В. Фотограмметрия. –
Москва: Недра, 1987. – 309 стр.
9. Быков Р.Е., Фрайер Р., Иванов К.В., Манцветов А.А. Цифровое
преобразование изображений. – Москва: Горячая линия-Телеком, 2003.
– 227 стр.
82
ПРИЛОЖЕНИЕ 1. ВАРИАНТЫ ЗАДАНИЙ ДЛЯ ВЫПОЛНЕНИЯ ГЕОМЕТРИЧЕСКОГО ТРАНСФОРМИРОВАНИЯ СНИМКА.
ВАРИАНТ
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
Название
img-файла
p1_2
p1_5
p1_3
p1_1
p2_4
p2_5
p3_5
p3_4
p2_2
p1_4
p3_3
p3_2
p3_1
p2_3
p2_1
ПРИЛОЖЕНИЕ 2. ВАРИАНТЫ ЗАДАНИЙ ДЛЯ ВЫПОЛНЕНИЯ ОРТОТРАНСФОРМИРОВАНИЯ СНИМКА.
Название фрагмента снимка и
координаты контрольных
точек.
15
4
18
Точка
p2002,
пара
снимков
1,2
3
Точка
p2001,
пара
снимков
0,1
12
Точка
p1006,
пара
снимков
1,2
2
Точка
p1005,
пара
снимков
1,2
15
Точка
p1004,
пара
снимков
1,2
1
Точка
p1003,
пара
снимков
0,1
Число
связующих
точек
Точка
p1002,
пара
снимков
0,1
Вариант
Название фрагмента снимка и координаты опорных точек.
Pt 2.
665534.3453
115137.0433
1937.9880
Pt 3
665621.8938
115375.3598
1941.6121
Pt 4
665304.2205
115560.2221
1964.0570
Pt 2
665003.7491
119471.1389
1964.6235
Pt 3
664328.1797
119233.4255
1996.6287
Pt 5
664251.3688
119037.1018
1990.6149
Pt2
668150.61
122404.68
1972.056
Pt2
668150.61
122404.68
1972.056
Pt2
668150.61
122404.68
1972.056
Pt 2
668402.5778
118510.9363
1880.9208
Pt 3
667917.5955
117994.2351
1884.6889
Pt 4
668104.2699
118889.4491
1895.8103
Pt 2
670743.252
119046.3593
1997.8586
Pt 4
670990.8787
118491.4566
2018.6760
Pt 3
670633.8227
118490.7914
2034.3120
Pt 2
671312.5644
115774.3454
1907.5248
Pt 4
670166.7720
116125.4597
1882.8912
Pt 3
670354.3834
115550.5878
1886.4154
Pt 2
671504.2822
123474.3582
1984.6181
Pt 3
671038.2103
123390.7810
1959.3059
Pt 2
671504.2822
123474.3582
1984.6181
Pt 5
665142.1149
115308.4675
1962.0192
Pt 4
664041.1887
119252.8084
2000.3467
Pt2
668150.61
122404.68
1972.056
Pt 2
668402.5778
118510.9363
1880.9208
Pt 4
670990.8787
118491.4566
2018.6760
Pt 2
671312.5644
115774.3454
1907.5248
Pt 2
671504.2822
123474.3582
1984.6181
83
5
14
6
12
7
16
8
18
9
14
10
14
11
13
12
17
Pt 2.
665534.3453
115137.0433
1937.9880
Pt 3
664328.1797
119233.4255
1996.6287
Pt2
668150.61
122404.68
1972.056
Pt 3
667917.5955
117994.2351
1884.6889
Pt 3
670633.8227
118490.7914
2034.3120
Pt 4
670166.7720
116125.4597
1882.8912
Pt 3
671038.2103
123390.7810
1959.3059
Pt 4
665304.2205
115560.2221
1964.0570
Pt 3
665621.8938
115375.3598
1941.6121
Pt 3
664328.1797
119233.4255
1996.6287
Pt 5
664251.3688
119037.1018
1990.6149
Pt2
668150.61
122404.68
1972.056
Pt2
668150.61
122404.68
1972.056
Pt 4
668104.2699
118889.4491
1895.8103
Pt 2
668402.5778
118510.9363
1880.9208
Pt 4
670990.8787
118491.4566
2018.6760
Pt 3
670633.8227
118490.7914
2034.3120
Pt 3
670354.3834
115550.5878
1886.4154
Pt 2
671312.5644
115774.3454
1907.5248
Pt 2
671504.2822
123474.3582
1984.6181
Pt 3
671038.2103
123390.7810
1959.3059
Pt 2.
665534.3453
115137.0433
1937.9880
Pt 2
665003.7491
119471.1389
1964.6235
Pt2
668150.61
122404.68
1972.056
Pt 4
668104.2699
118889.4491
1895.8103
Pt 2
670743.2524
119046.3593
1997.8586
Pt 3
670354.3834
115550.5878
1886.4154
Pt 2
671504.2822
123474.3582
1984.6181
Pt 5
665142.1149
115308.4675
1962.0192
Pt 3
665621.8938
115375.3598
1941.6121
Pt 2.
665534.3453
115137.0433
1937.9880
Pt 4
665304.2205
115560.2221
1964.0570
Pt 4
664041.1887
119252.8084
2000.3467
Pt 2
665003.7491
119471.1389
1964.6235
Pt 4
664041.1887
119252.8084
2000.3467
Pt 5
664251.3688
119037.1018
1990.6149
Pt2
668150.61
122404.68
1972.056
Pt2
668150.61
122404.68
1972.056
Pt2
668150.61
122404.68
1972.056
Pt2
668150.61
122404.68
1972.056
Pt 3
667917.5955
117994.2351
1884.6889
Pt 2
668402.5778
118510.9363
1880.9208
Pt 2
668402.5778
118510.9363
1880.9208
Pt 3
667917.5955
117994.2351
1884.6889
Pt 3
670633.8227
118490.7914
2034.3120
Pt 2
670743.2524
119046.3593
1997.8586
Pt 3
670633.8227
118490.7914
2034.3120
Pt 4
670990.8787
118491.4566
2018.6760
Pt 4
670166.7720
116125.4597
1882.8912
Pt 3
670354.3834
115550.5878
1886.4154
Pt 2
671312.5644
115774.3454
1907.5248
Pt 4
670166.7720
116125.4597
1882.8912
Pt 3
671038.2103
123390.7810
1959.3059
Pt 3
671038.2103
123390.7810
1959.3059
Pt 3
671038.2103
123390.7810
1959.3059
Pt 2
671504.2822
123474.3582
1984.6181
84
13
14
15
13
Pt 5
665142.1149
115308.4675
1962.0192
Pt 3
664328.1797
119233.4255
1996.6287
Pt2
668150.61
122404.68
1972.056
Pt 3
667917.5955
117994.2351
1884.6889
Pt 2
670743.2524
119046.3593
1997.8586
Pt 2
671312.5644
115774.3454
1907.5248
Pt 2
671504.2822
123474.3582
1984.6181
17
Pt 4
665304.2205
115560.2221
1964.0570
Pt 2
665003.7491
119471.1389
1964.6235
Pt2
668150.61
122404.68
1972.056
Pt 4
668104.2699
118889.4491
1895.8103
Pt 3
670633.8227
118490.7914
2034.3120
Pt 4
670166.7720
116125.4597
1882.8912
Pt 3
671038.2103
123390.7810
1959.3059
16
Pt 3
665621.8938
115375.3598
1941.6121
Pt 5
664251.3688
119037.1018
1990.6149
Pt2
668150.61
122404.68
1972.056
Pt 4
668104.2699
118889.4491
1895.8103
Pt 2
670743.2524
119046.3593
1997.8586
Pt 3
670354.3834
115550.5878
1886.4154
Pt 2
671504.2822
123474.3582
1984.6181
85
ПРИЛОЖЕНИЕ 3. НАЗВАНИЯ КОНТРОЛЬНЫХ И ОПОРНЫХ ТОЧЕК ДЛЯ ОРТОТРАНСФОРМИРОВАНИЯ / НАЗВАНИЕ
ФРАГМЕНТА СНИМКА.
P1002 / Pt2
P1002 / Pt 4
P1002 / Pt 3
P1002 / Pt 5
86
P1003 / Pt 2
P1003 / Pt 3
P1003 / Pt 4
P1003 / Pt 5
87
P1004 / Pt2
P1005 / Pt2
P1005 / Pt3
P1005 / Pt4
88
P1006 / Pt2
P1006 / Pt3
P1006 / Pt4
P2001 / Pt2
89
P2001 / Pt3
P2002 / Pt2
P2001 / Pt4
P2002 / Pt3
90
Учебное издание
Методические указания к выполнению курсовой работы по
фотограмметрии и дистанционному зондированию (для студентов 4
курса дневной формы обучения спец. 7.070908 «Геоинформационные
системы и технологии»).
Составитель : Лисицин Владимир Эдуардович
Ответственный за выпуск: Л.Г. Запара
Печатается в авторской редакции
План 2006, поз.396
Подп. к печ.
Печать на ризографе.
Зак.№
Формат 60х84х1/16
Уч.-изд. л. 2,5
Цена договорная
Бумага офисная
Тираж 70 экз.
61002, Харков, ул.Революции,12
Сектор оперативной полиграфии ИВЦ ХНАГХ