теории зрения с использованием оптико-электронных или програм- |
|||||||
мных блоков, называемых корреляторами. В ее основе лежит по- |
|||||||
нятие образа – произвольной по форме и размерам области снимка |
|||||||
|
|
|
вместе со всей имею- |
||||
|
|
|
щейся |
информацией. С |
|||
Образ R |
|
Образ R |
математической |
точки |
|||
|
|
Зона |
зрения |
|
образ |
пред- |
|
|
|
поиска |
ставляет |
собой |
много- |
||
|
|
|
мерный вектор R как со- |
||||
Левый снимок |
Правый снимок |
вокупность |
элементов |
||||
|
изображения (пикселей), |
||||||
|
|
|
|||||
|
|
|
каждый из которых ха- |
||||
Рис. 9.10. Автоматическая идентификация точек |
рактеризуется своим по- |
||||||
|
|
|
ложением и оптической |
||||
плотностью ij согласно формуле (9.6). |
|
|
|
|
|
||
Опознавание точки левого снимка на правом сводится к опреде- |
|||||||
|
|
|
|
И |
|||
лению некоторого образа R на левом снимке и поиску на правом |
|||||||
снимке такого образа R , чтобы расхождение между ними было ми- |
|||||||
нимальным: |
|
R – R =Дmin. |
|
|
|||
|
|
|
(9.7) |
||||
Практически для автоматического опознавания точки необхо- |
|||||||
димо: |
|
|
|
|
|
|
|
1. Вы рать на левом снимке о раз R |
в виде области размером |
||||||
|
|
А |
|
|
|
точка |
|
n n п кселей, в центре которой размещена опознаваемая |
|||||||
(р с. 9.10), |
определ ть его характеристики, на основе которых бу- |
||||||
дет выполняться проверка условия (9.7). |
|
|
|
|
|
||
2. Определбть зону по ска размером m m пикселей (m n) веро- |
|||||||
ятного расположен я |
скомой точки на правом снимке (см. рис. 9.10) с |
||||||
коорд натами центра xп xл − bсн; yп yл. |
|
|
|
|
|
||
3. Последовательным перемещением области размером n n пик- |
|||||||
селейграницахв зоны поиска размером m m создать на правом сним- |
|||||||
ке серию образов R и определить характеристики каждого из них с |
|||||||
целью проверки условия (9.7). |
|
|
|
|
|
||
4. Сопоставить характеристики каждого образа R с характеристи- |
|||||||
кой-эталоном вектора R. Искомая точка на правом снимке будет распо- |
|||||||
С |
|
|
|
|
|
|
|
ложена в центре образа R , для которого выполняется условие (9.7). |
|||||||
21
Установление степени соответствия векторов R и R представляет основную трудность и может быть выполнено различными путями. Например, один из них основан на расчете для образа-эталона R и каждого образа R правого снимка ковариационной матрицы
c11 |
c12 |
c13 |
|
|
|
|
|
|
|
|
|
C = c21 |
c22 |
c23 |
|
, |
(9.8) |
|
c32 |
c33 |
|
|
|
c31 |
|
|
|
||
где |
|
|
|
|
|
c11 = k (xi − x0 )2 , |
c12 = c21 = k (xi |
− x0 )( yi − y0 ), |
x0 = xi |
/ N |
|||||||||||||
c22 = k ( yi − y0 ) |
2 |
, |
|
|
|
|
|
|
|
)( |
|
− |
|
|
y0 = yi |
|
|
|
c = c |
31 |
= k |
|
(x |
i |
− x |
i |
0 |
), |
/ N |
||||||
|
|
|
13 |
|
|
0 |
|
|
|
|
|
|
|||||
c33 = k ( i − 0 )2 , |
c23 = c32 = k ( yi |
|
|
|
|
|
|
|
0 = i |
||||||||
− y0 )( i |
− 0 ), |
/ N |
|||||||||||||||
|
|
|
k =1/(N −1), |
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
И |
|
где N – число элементов в образах R и R ; |
xi, yi, i – координаты эле- |
||||||||||||||||
ментов образа и их оптические плотности. |
|
|
|
|
|
|
|
||||||||||
Если C0 – ковариационная матрица |
(9.8) |
|
образа R на |
левом |
|||||||||||||
|
|
|
|
|
|
|
|
|
Д |
|
|||||||
снимке, то условие (9.7) удет выполненным для вектора R , для ко- |
|||||||||||||||||
торого разность матриц V = C0–Ci минимальна.
Сгде 0ихи 0 – средние оптические плотности элементов зон (фрагментов), соответствующ образам левого R и правого R снимков.
Более широко применяется метод идентификации точек сним- |
||||||||||||||||
ков, основанный на расчете коэффициентов корреляции между опти- |
||||||||||||||||
ческ ми плотностями элементов о раза R и каждого из образов R на |
||||||||||||||||
|
|
|
А |
|
|
|
||||||||||
правом сн мке с спользованием формулы |
|
|
|
|
||||||||||||
|
|
|
|
( |
|
− |
|
)( |
− ) |
|
|
|
|
|||
r = |
|
i |
|
|
|
0 |
i |
0 |
|
|
, |
(9.9) |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||
б2 |
( − |
) |
2 |
|
|
|||||||||||
|
|
|
( |
i |
− |
) |
|
|
|
|
||||||
|
|
|
|
|
|
0 |
|
|
|
|
i |
0 |
|
|
|
|
Коэффициенты, подсчитанные для сочетаний образа R со всеми образами R , образуют корреляционную матрицу: максимальное значение ее элемента соответствует наиболее тесной связи оптических плотностей сравниваемых участков левого и правого снимков и, следовательно, выполнению условия (9.7). Так что искомая точка лежит в центре образа R с максимальным коэффициентом корреляции.
22
Некоторые способы предполагают отыскание искомой точки по максимуму корреляционной функции, составленной на основе анализа элементов корреляционной матрицы с частными коэффициентами корреляции (9.9), соответствующей искомому образу R .
Размер области снимка, отождествляемой с образом R, обычно составляет 20 20 пикселей; при малом числе контуров местности он увеличивается до 40 40 пикселей. Программы обработки обычно запрашивают полуразмер корреляционной матрицы (образа R). Размер области поиска на правом снимке выбирается таким, чтобы он был в два раза больше ожидаемого смещения точки из-за влияния рельефа местности. Так, при высоте фотографирования H = 2000 м, максимальном превышении рельефа над средней плоскостью h = 50 м и расстоянии от центра снимка до угла рабочей площади r = 100 мм (формат кадра 23 23 см) будем иметь h = 2,5 мм, а размер области
Рассмотренные варианты установления степениИсоответствия образов R и R далеко не единственные: в специальной литературе име-
поиска – 5 5 мм.
ются упоминания о методеАиерархическойДрелаксации (методе пирамид HRC), методе вертикальной линии (ULL), методе динамического программирования и т.п. Некоторые из них основаны на анализе тех или иных признаков, характеризующих выделенные по группе пикселей элементы: их форму, взаимное расположение, ориентировку и др. Все методы идентификации точек о ъединяют два обстоятельства:
• отождествление выполняется на основе более или менее строгого анал за элементов изо ражений – их оптических плотностей
геометр ческого положен я;
|
• корреляторов, |
о еспечивающих 100% гарантию качества ав- |
|||
томат ческой дент ф кац и точек, не существует. |
|||||
|
|
обстоятельствабпредопределяют целесообразность стерео- |
|||
скоп ческого контроля положения измерительной марки после рабо- |
|||||
|
|
|
|
|
ты коррелятора. В большинстве случаев |
|
Эти |
опытный наблюдатель ее обязательно |
|||
|
чуть-чуть подправит по высоте. |
||||
|
|
|
|
|
|
|
|
|
|
|
Искажения точек, вызванные вли- |
|
|
|
|
|
|
|
Рис. 9.11. Геометрическое |
янием угла наклона снимка и рельефа |
|||
|
местности, вызывают смещения отдель- |
||||
|
положение элементов |
|
|||
|
образа R на правом снимке |
ных элементов растра, и это не может не |
|||
С |
|
сказываться на качестве идентификации. |
|||
|
искажается рельефом |
|
|||
23
Так, если на левом снимке образ R всегда представлен квадратом, то на правом снимке (рис. 9.11) из-за влияния угла наклона и рельефа местности соответствующие элементы образа R получают смещения и в общем случае образуют многоугольник. Поскольку анализируемый образ R геометрически не всегда подобен образу R, то вероятность полного совпадения искомой точки уменьшается. Еще более сложная ситуация возникает в случае, если изображения левого и правого снимков развернуты одно относительно другого.
Однако из-за малости элементов изображения эти искажения не столь значительны, а их влияние можно компенсировать геометрической коррекцией положения пика корреляции и соответствующей ему точки правого снимка. Одно из таких решений основано на том, что смещения элементов образа являются следствием неровностей рельефа в пределах соответствующего участка местности и интерпретиру-
Основная задача фотограмметрииД– определение формы, размеров и пространственного положения объектов местности – при использовании цифровых изо ражений решается путем определения
ются как параллактические смещения. |
И |
|
|
9.7. Фотограмметрическая обработка |
|
цифровых снимков |
|
координат этих о ъектов или их элементов с применением рассмот- |
|||
ренных ранее аналитических спосо ов. Важнейшей их особенностью |
|||
является о ъед нен е выч слительной обработки (внутреннее, вза- |
|||
|
|
|
А |
мное, внешнее ор ент рование и т. п.) со сбором нужной для этого |
|||
нформац |
– |
змерен ем координат необходимых точек и автома- |
|
т ческой х |
дент ф кац ей на смежных снимках. Это обстоятель- |
||
ство в сочетан |
с макс мальной автоматизацией технологических |
||
|
|
б |
|
процессов делают ц фровую фотограмметрическую обработку высо- |
|||
коэффект вной |
технолог чной, не требующей высокой квалифика- |
||
ции исполнителей. |
|
||
рассмотрены некоторые элементы вычислительной обра- |
|||
Ниже |
|
||
С |
|
|
|
ботки цифровых снимков, аналогичные процессам аналитической обработки данных в главе 5.
24
Внутреннее ориентирование цифровых изображений (снимков) выполняется с целью установления соответствия между координатными
системами растра oPxPyP и снимка oxy (рис. 9.12). |
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
y |
|
|
|
|
|
|
|
|
|
|
|
|
Математическая |
|
|
|
модель |
внутреннего |
|||||||||||||||||||||
op |
|
|
|
|
|
|
|
|
|
|
|
|
xp |
|
|
|
ориентирования |
|
цифрового |
изображения |
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
аналогична применяемой при аналитическом |
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
построении сетей фототриангуляции и пред- |
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
o |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
полагает определение параметров ортого- |
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
нального, аффинного или проективного пре- |
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||
yp |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
образования по избыточному числу измере- |
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||
Рис. 9.12. Координатные |
|
|
|
ний методом наименьших квадратов. |
|
||||||||||||||||||||||||||||||||||||||||||
|
|
|
системы снимка oxy |
|
|
|
|
|
|
|
|
Так, |
наиболее часто применяемые аф- |
||||||||||||||||||||||||||||||||||
|
|
|
|
и растра opxpyp |
|
|
|
|
|
финные преобразования основаны на исполь- |
|||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||
зовании следующих формул связи исходных xPyP и преобразованных |
|||||||||||||||||||||||||||||||||||||||||||||||
x,y координат: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
И |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
= |
|
a |
0 |
|
+ |
a |
|
a |
|
|
|
xp |
|
= |
a |
0 |
|
+ P |
xp |
|
, |
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
(9.10) |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
y |
|
|
|
b0 |
|
|
b1 |
|
b2 |
|
yp |
|
|
b0 |
|
|
|
|
|
|
yp |
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Д |
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
А |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
где ai и bi (i=0,1,2) – коэффициенты аффинного преобразования, |
|||||||||||||||||||||||||||||||||||||||||||||||
определяемые из решения системы уравнений поправок, составлен- |
|||||||||||||||||||||||||||||||||||||||||||||||
ных по результатам измерения координат оптических меток снимка. |
|||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
Л нейные коорд наты xP ,yP точек снимка в системе oPxPyP свя- |
||||||||||||||||||||||||||||||||||||||||||
заны с геометр ческ |
|
|
разрешением растра и номерами столбцов iX |
||||||||||||||||||||||||||||||||||||||||||||
|
|
|
строк iY |
|
|
|
|
зав |
|
с |
мостями: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
бx = i , |
|
|
i |
X |
= x |
p |
/ , |
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
p |
|
|
|
X |
|
|
|
|
|
и |
|
|
|
|
|
|
|
|
|
|
|
|
|
(9.11) |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
yp |
= iY |
|
|
|
|
|
iY |
= yp |
/ |
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
простымиДля обратного преобразования координат точек из системы oxy |
||||||||||||||||||||||||||||||||||||||||||
в систему растра oPxPyP используются формулы |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||
С |
|
|
|
xp |
|
= |
|
A |
A |
|
|
|
x − a |
0 |
|
= P−1 |
|
|
x − a |
0 |
|
, |
|
(9.12) |
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||
|
|
|
y |
|
|
|
|
|
1 |
2 |
|
|
|
|
|
|
|
y |
|
|
|
|
|||||||||||||||||||||||||
|
|
|
p |
|
|
|
|
B |
B |
|
|
|
y − b |
|
|
|
|
|
|
|
|
|
− b |
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
1 |
2 |
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
где P и P-1 – прямая и обратная матрицы аффинного преобразования.
25