- •ПРЕДИСЛОВИЕ
- •Глава 9. МЕТОДЫ ЦИФРОВОЙ ФОТОГРАММЕТРИИ
- •9.1. Понятие о цифровом изображении
- •9.2. Характеристики цифрового изображения
- •9.4. Источники цифровых изображений
- •9.5. Стереоскопические наблюдения и измерения цифровых изображений
- •9.7.1. Внутреннее ориентирование снимков
- •9.8.1. Способы представления цифровой модели рельефа
- •9.9. Ортотрансформирование снимков
- •9.10.4. Создание цифровых трансформированных изображений
- •9.10.5. Создание цифровых фотопланов
- •Глава 10. МЕТОДЫ ИНЕРЦИАЛЬНОЙ И СПУТНИКОВОЙ НАВИГАЦИИ
- •10.2. Инерциальные навигационные системы
- •10.2.1. Общие принципы инерциальной навигации
- •10.2.4. Обработка инерциальных данных
- •10.3. Спутниковые навигационные системы
- •10.3.1. Действующие и разрабатываемые СНС
- •10.3.2. Основные компоненты СНС
- •10.3.3. Навигационные сигналы GPS, ГЛОНАСС и Galileo
- •10.4. Интеграция инерциальных и спутниковых систем
- •10.4.1. Достоинства и недостатки навигационных систем
- •10.4.2. Фильтр Калмана
- •10.4.3. Элементы модели интеграции ИНС и СНС
- •Глава 11. МЕТОД АЭРОГЕОДЕЗИЧЕСКИХ РАБОТ НА ОСНОВЕ ВОЗДУШНОЙ ЛАЗЕРНОЙ ЛОКАЦИИ И ЦИФРОВОЙ АЭРОФОТОСЪЁМКИ
- •11.2.1. Установка и наладка оборудования на борту летательного аппарата
- •11.2.2. Геодезическое обеспечение аэросъемочных работ
- •11.2.3. Производство измерений на борту летательного аппарата
- •11.2.7. Тематическая обработка
- •11.2.8. Обработка цифровых фотоснимков
- •11.3. Программный комплекс ALTEXIS
- •11.4. Основные возможности воздушных сканеров ALTM
- •11.5. Инструментальные средства лазерной локации
- •11.5.2. Методы выполнения развертки
- •Глава 12. СИСТЕМЫ НАЗЕМНОГО МОБИЛЬНОГО ЛАЗЕРНОГО СКАНИРОВАНИЯ
- •12.1. Особенности и преимущества наземных мобильных систем
- •12.3. Процесс съемки и получаемый результат
- •12.4. Испытание системы StreetMapper
- •Глава 13. ГЕОИНФОРМАЦИОННОЕ ОБЕСПЕЧЕНИЕ ТЕРРИТОРИИ ГОРОДА
- •14.1. Обзор информации о БПЛА
- •Библиографический список
9.10.4. Создание цифровых трансформированных изображений
Рассмотрим процесс цифрового трансформирования исходного снимка на примере его преобразования в цифровое трансформированное изображение, представляющее собой ортогональную проекцию местности на горизонтальную плоскость (ортофотоснимок).
Принципиальная схема этого |
|
|
Z |
Y |
|
|
||||
процесса представлена на рис. 9.23. |
|
|
|
|
|
|||||
|
z |
|
y |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
В |
результате |
цифрового |
|
S |
|
|
|
|
|
|
|
|
|
|
x |
X |
|
||||
трансформирования |
исходный |
|
|
|
|
|
||||
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
||||
снимок преобразуется в цифровое |
|
|
|
|
|
|
y xp |
|||
изображение местности, |
представ- |
|
|
|
|
op |
|
|
||
|
|
|
|
|
|
|
||||
ляющее собой ортогональную про- |
|
|
|
|
a |
|
x |
|||
|
|
|
|
|
|
|||||
екцию местности на горизонталь- |
yp |
|
|
|
|
|
|
|||
ную плоскость. |
|
|
|
|
|
|
X, Y, Z |
|
|
|
Исходными материалами при |
|
|
|
|
И |
|||||
цифровом |
|
трансформировании |
Pij |
|
|
|
||||
|
|
|
|
|
|
|
||||
снимков служат: |
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|||
• цифровое изображение ис- |
Y |
|
|
|
|
|
|
|||
ходного фотоснимка; |
|
Z |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|||
• цифровая модель рельефа |
|
|
|
|
|
|
X |
|||
|
|
|
|
X, Y |
|
|||||
|
|
|
|
Д |
|
|||||
(ЦМР); |
|
|
|
|
Y |
|
3 |
|
|
1 |
• |
|
элементов внут- |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|||
реннего и внешнего ориентирова- |
|
|
|
|
|
A |
|
|||
|
|
4 |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
н я сн мков; |
|
|
O |
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
X |
|||
|
|
|
А |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
||
• значен е параметров внут- |
|
|
|
|
|
|
|
|||
реннего ор ент рован я снимка в |
Рис. 9. 23. Схема трансформирования |
|||||||||
с стеме |
коорд нат |
ц фрового |
|
|
|
|
|
|
|
|
|
|
б |
|
|
|
|
|
|
|
|
значение |
|
|
|
|
|
|
|
|
||
С |
|
|
|
|
|
|
|
|
|
зображен я.
В больш нстве случаев при трансформировании снимков используется ц фровая модель местности в виде матрицы высот, представляющей собой регулярную сетку квадратов на местности, стороны которых параллельны осям X и Y системы координат объекта 0XYZ. Координаты и высоты узлов сетки квадратов определены в системе координат объекта.
Для формирования ЦМР в виде матрицы высот в большинстве случаев используют цифровые модели рельефа, созданные в результате стереофотограмметрической обработки снимков или по уже существующим топографическим картам.
47
Наиболее распространённым методом построения цифровых |
|||||||
моделей рельефа местности является в настоящее время метод триан- |
|||||||
гуляции Делоне, в котором рельеф местности представлен в виде про- |
|||||||
странственной сети треугольников, координаты и высоты вершин |
|||||||
которых определены в системе координат объекта. Рельеф местности |
|||||||
в пределах треугольника в этом виде ЦМР аппроксимируется плоско- |
|||||||
стью, проведённой через его вершины. |
|
|
|
||||
При формировании ЦМР этим методом по высотным пикетам |
|||||||
треугольники генерируются под условием, чтобы в окружность, про- |
|||||||
веденную через вершины треугольников, не должны попадать верши- |
|||||||
ны других треугольников. |
|
|
|
|
|
||
По цифровой модели рельефа в виде триангуляции Делоне мож- |
|||||||
но сформировать ЦМР в виде матрицы высот. |
|
|
|||||
На рис. 9.24 представлены фрагменты ЦМР местности в виде |
|||||||
|
|
|
матрицы высот и в виде триангуляции Де- |
||||
|
|
|
лоне. |
|
|
И |
|
|
|
|
|
|
|
||
k |
|
|
Для определения высоты узла i |
мат- |
|||
|
|
рицы высот, по координатам xi и yi |
этого |
||||
|
|
|
|||||
|
i |
m |
узла в l системе координат объекта находят |
||||
|
|
вершины треугольника триангуляции Де- |
|||||
|
|
|
|||||
|
|
|
лоне, в котором находится узел i. |
|
|||
|
|
|
|
|
Д |
|
|
l |
|
|
Значение высоты узла i определяют |
||||
|
|
по формуле |
|
|
|
||
|
|
|
|
|
|
||
Рис. 9.24. Фрагмент ЦМР |
|
|
|
|
|
||
|
|
|
|
Zi |
= A + BXi + CYi . |
(9.34) |
|
|
|
|
А |
|
|
||
Выражен е (9.34) представляет собой уравнение плоскости, |
|||||||
проведенной через верш ны l, k, m треугольника, внутри которого |
|||||||
наход тся узел i. |
|
|
|
|
|
||
|
|
б |
|
|
|
|
|
Коэфф ц енты уравнения (9.34) A, B и C получают в результате |
|||||||
решен я с стемы з трёх уравнений |
|
|
|
|
|||
иA + BX + CY − Z = 0, |
|
(9.35) |
|||||
составленных по значениям координат X, Y и высот Z каждой из вер- |
|||||||
шин l, k, m треугольника. |
|
|
|
|
|
||
С |
|
|
|
|
|
|
|
|
|
|
48 |
|
|
|
|
Цифровое трансформирование снимка выполняется следующим образом.
Сначала формируется прямоугольная матрица цифрового
трансформированного изображения, строки и столбцы которой параллельны осям X и Y системы координат объекта, а координаты одного из углов матрицы заданы в этой же системе координатИ. Размер
элементов (пикселей) матрицы обычно выбирают приблизительно равным величине ×m, в которой:
• − размер пикселя цифрового изображения исходного снимка;
• m − знаменатель среднего масштаба снимка.
Значения координат начала системы координат создаваемой
матрицы выбирают кратными величине элементов матрицы.
|
|
|
|
|
|
ля формирования |
||
|
|
|
|
|
|
цифрового трансформи- |
||
|
|
|
|
А |
|
|
||
|
|
|
|
|
|
рованного |
изображения |
|
|
|
|
|
|
|
каждому элементу циф- |
||
|
|
|
|
|
|
рового изображения a ij |
||
|
|
|
|
|
|
необходимо |
присвоить |
|
|
|
|
ображен |
|
оптическую |
плотность |
||
|
|
|
|
изображения |
соответ- |
|||
|
|
|
|
|
Дствующего участка объ- |
|||
|
Определение |
|
|
екта на исходном циф- |
||||
|
|
|
ровом снимке. Эта опе- |
|||||
|
Рис. 9.25. Определение значения высоты точки i |
рация выполняется сле- |
||||||
|
|
|
|
|
|
дующим |
образом. По |
|
|
|
|
|
|
|
|||
|
значен ям |
ндексов i |
j элементов матрицы a ij определяются коор- |
|||||
|
д наты X, |
Y центра соответствующего пикселя цифрового трансфор- |
||||||
С |
|
|
|
|
|
|
||
|
м рованного з |
я в системе координат объекта. |
|
По коорд натам Xi, Yi точки объекта, соответствующей центру п кселя, по ц фровой модели рельефа определяется высота этой точки Zi.
значения Zi по ЦМР в виде матрицы высот выполняется методом билинейного интерполирования (рис. 9.25).
На рис. 9.25 X = Xi – X1; Y= Yi – Y1, где X1 и Y1 – координаты узла i цифровой модели рельефа.
Высота точки Zi вычисляется по формуле
Zi = Z K + ((ZM − Z K ) D) X i , |
(9.36) |
49
в которой
Z |
|
= Z |
|
|
+ |
Z2 − Z1 |
|
Y |
|
||
|
|
|
|
D |
|
||||||
|
K |
|
1 |
|
|
|
i |
||||
|
|
|
|
|
|
|
Z4 − Z3 |
|
|
|
|
Z |
|
= Z |
|
+ |
|
Y |
|||||
M |
3 |
|
|
||||||||
|
|
|
|
|
D |
|
|
i |
|||
|
|
|
|
|
|
|
|
|
|
По координатам Xi, Yi, Zi и значениям элементов внутреннего и внешнего ориентирования снимка вычисляют координаты х,у соответствующей точки на исходном цифровом снимке в системе координат снимка Sхуz.
Вычисления производят по формулам: |
|
|
|
|||||||
|
|
x = x0 |
− f |
x * |
|
|
|
|||
|
|
|
|
|
|
, |
|
(9.37) |
||
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|||
|
|
|
|
z * |
|
|||||
|
|
y = y0 |
− f |
|
y * |
|
|
И |
||
|
|
|
|
|
|
|
||||
|
|
|
|
|
z * |
|
||||
в которых |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
x * |
X − X S |
|
|||||||
|
|
|
T |
|
|
|
|
|
|
|
|
y * = |
A |
|
Y − YS |
. |
|
||||
|
z * |
|
Z − Z |
|
|
|
||||
|
|
|
|
|
|
|
|
S |
|
|
|
|
|
|
|
Д |
|
||||
По |
х,у и значениям параметров внутреннего ориен- |
сн мка в с стеме коорд нат цифрового изображения осхсус.
тирования цифрового изо ражения определяют координаты точки |
|
|
А |
б |
|
координатам |
|
В случае спользован |
я аффинных преобразований при выпол- |
нен внутреннего ор ент |
рования определение координат выпол- |
няют по формуле: |
|
|
|
|
|
|
|
|
|
|
|
|
|
С |
x |
|
|
x − a |
|
A |
A |
x − a |
|
||||
|
|
c |
= P−1 |
|
0 |
|
= |
1 |
2 |
|
0 |
|
|
|
|
y |
c |
|
|
y − b |
|
|
B |
B |
|
y − b |
|
|
|
|
0 |
1 |
2 |
0 |
Затем по координатам хс и ус вычисляют пиксельные координаты точки:
xp y p
50
По значениям пиксельных координат xp, yp точки цифрового изображения снимка, которая является проекцией центра пикселя матрицы цифрового трансформированного изображения, находят
ближайшие к этой точке четыре пикселя цифрового изображения |
|||
снимка. Затем методом билинейной интерполяции, изложенным вы- |
|||
|
|
|
И |
ше, по формулам (9.20) определяют значение оптической плотности |
|||
Di или цвета соответствующего пикселя матрицы цифрового транс- |
|||
формированного изображения. |
|
|
|
При этом значение величин xp, yp в формулах (6.20) опреде- |
|||
ляют по формулам: |
Д |
||
|
xp = xpi |
||
|
− xpk |
|
|
|
yp = ypi |
|
|
|
− ypk |
|
|
|
А |
|
|
Таким же образом определяют оптические плотности или цвет |
|||
всех остальных пикселей цифрового трансформированного изобра- |
|||
жения. |
|
|
|
Помимо метода билинейной интерполяции для формирования |
|||
|
б |
|
|
цифрового трансформированного изображения |
применяют метод |
||
«ближайшего соседа», в котором по пиксельным координатам xp,yp |
|||
находят пиксель цифрового изо ражения снимка, |
на который проек- |
тируется точка, соответствующая центру пикселя цифрового транс- |
||
изображения |
|
|
формированного изо ражения; значение его оптической плотности |
||
ли цвета пр сва вается п кселю цифрового трансформированного |
||
зображен я. |
|
|
Метод « |
л жайшего соседа» позволяет сократить время форми- |
|
рован я ц фрового трансформированного изображения по сравнению |
||
С |
|
|
с методом линейной |
нтерполяции, однако изобразительные свой- |
|
ства форм руемого ц фрового изображения при этом ухудшаются. |
||
Если превышен |
точек на участке местности, изображенной на |
|
снимке, незначительны, то при создании цифрового трансформиро- |
||
ванного |
поступают упрощенно. Значения высот точек |
местности, соответствующие центрам пикселей трансформированного изображения, принимаются равными для всех смежных участков.
В этом случае нет необходимости в создании цифровой модели рельефа местности, так как трансформированное цифровое изображение представляет собой центральную проекцию исходного снимка на горизонтальную плоскость, расположенную на высоте Z, равной
51
среднему значению высоты участка местности. Такой метод трансформирования допустим в случае, если ошибки в положении точек на трансформированном изображении, вызываемые рельефом местности, не превышают допустимых значений.
ошибки в положении точек на трансформированном изображении не будут превышать установленного допуска Rmax, можно определить
Величины максимально допустимых значений превышений точек местности hmax относительно средней плоскостиИ, при которых
по формуле
|
f |
Д |
|
h max = |
r |
Rmax , |
(9.38) |
где f – фокусное расстояние съёмочной камеры; r – расстояние на исходном снимке от главной точки до точки на снимке.
Как следует из формулы (9.38), величина ошибки Rmax прямо пропорциональна значению r. Поэтому при определении hmax измеряется значение r до наиболее удаленной от главной точки снимка точки, участ-
вующей в формировании трансформированного изображения.
Следует отметить, что формулу (9.38) используют только в случае, если трансформирование выполняется по снимкам, углы наклона которых не превышают 3 – 5 .
Аналогичным о разом можно определить величину допустимой ошибки hmax рельефа определения высот точек местности, соответствующих центрам пикселей трансформированного изображения, по
ц фровой модели. |
|
А |
|
|||
|
|
|
|
|||
|
|
h max = |
f |
R max. |
(9.39) |
|
|
|
|
||||
|
|
|
|
r |
|
|
В случае |
бтрансформирование снимков выполняется с це- |
|||||
лью создан я |
ли обновлен я карт и планов, значение Rmax |
выбира- |
||||
если |
|
|
|
|
||
С |
|
|
|
|
|
|
ется равным вел ч не 0,2 мм на карте или плане:
где М – знаменатель масштаба создаваемой карты.
При создании цифровых трансформированных изображений местности в проекции карты плановые координаты узлов цифровой модели рельефа определяют в системе координат карты. В России то-
52