tymkul (729084), страница 3
Текст из файла (страница 3)
.
P’(N, L) = [ 1 - z / x2 + y2 ) k4+z2 ] [ 2 x2 / (x2 + y2) - 1 ] . ( 53 )
Для удобства вывода выражения для P’(N, L) в сферической системе координат, воспользуемся переводом компонент в другую систему координат:
X = sin cos ;
Y = sin cos ; ( 54 )
Z = cos .
.
( n* rн ) = 1 / 1 + k4 tg2
Поскольку угол для конуса является постоянным и равным отношению радиуса основания к высоте, то справедливо выражение:
tg = a / c. ( 55 )
Подставив формулу ( 55 ) в выражения ( 54 ), получим :
.
( n* rн ) =1 / 1 + k2, ( 56 )
( nxy i ) = cos. ( 57 )
Тогда
.
P’(N, L) = [ 1 - 1 / 1 + k2 ] cos. ( 58 )
Таким образом, формула ( 58 ) является оптико-математической моделью поляризационного тепловизионного изображения конуса. При этом угол связан с номером строки L и с номером элемента строки N кадра зависимостью:
= arctg[( L - L0 ) / ( N - N0 )]. ( 59 )
Эта формула получена на основе геометрии перевода объёмного изображения на плоский кадр и логических рассуждений.
2.5. Анализ результатов исследования поляризационных
тепловизионных изображений объектов простой формой.
Практической целью моделирования поляризационных тепловизионных изображений объектов является распознование их формы внутри контура. Если проанализировать полученные модели изображений эллипсоидов с различными значениями коэффициента сжатия k, то можно заметить по поверхности сферы равномерное распределение степени поляризации Р’ от 0 до 1 вдоль горизонтальной линии от центра к краю и от 0 до -1 вдоль вертикальной линии от центра к краю.
По мере вытягивания эллипсоида ( к >1 ) область небольших по модулю значений степени поляризации | P’ | < 0.09 снижается , при этом область значений 1< | P’ | < 0.09 расширяется. При сжатии эллипсоида наблюдается обратная картина. Так для диска почти по всей поверхности значения P’ близки к нулю и только область, близкая к краю, занята значениями | P’ |, близкими к 1.
Поляризационные тепловизионные изображения конуса также дают возможность интерпретации его формы внутри конуса. Распределение степени поляризации в модели диска, полученное по формулам для сильно сжатого конуса аналогично распределению в модели, полученной по формулам для сильно сжатого эллипсоида. Однако модель самого конуса имеет очевидное отличие от объектов в виде эллипсоида по распределению степени поляризации. Здесь наблюдаются одинаковые значения степени поляризации вдоль выбранного диаметра. Причём, чем более вытянут конус, тем больше в модели изображения области с | P’| близкими к 1 и наоборот.
Таким образом, приведённый анализ поляризации тепловизионных изображений объектов показал, что имеется существенная зависимость формы объектов внутри их контура от значений степени поляризации P’ по наблюдаемым участкам поверхности объектов.
2.6. Модифицированный метод моделирования
поляризационных тепловизионных изображений.
В приведённых выше математических выкладках для вывода основных формул моделирования поляризационных тепловизионных изображений использовались поляризационные свойства собственного излучения объекта. Эти свойства обычной тепловизионной обработкой выделить невозможно, поэтому необходимы дополнительные технические средства в качестве анализатора поляризационного излучения. Таким анализатором может служить поляризационный фильтр, азимут поляризации которого будет изменяться от 00 до 3600 . Формировать оптико-математическую модель изображения тепловизионной системы с поляризационным фильтром можно модифицированным методом моделирования не основе вектор-параметра Стокса и влияния на излучение от объекта поляризационного фильтра. Причём исходным выражением для видеосигнала будем считать:
( N, L ) = ( 1/ cosN,L)dS(N,L)SW(,T,y,z)a(d( 60 );
Вектор-параметр Стокса, описанный в разделе 2.1 формулой ( 4 ), в нормированном виде выглядит следующим образом:
1
Uj (N, L) = U0 | P cos 2 t | ( 61 );
| P sin 2 t |
0
где U0 - суммарный видеосигнал при азимутах поляризации излучения t=00 и t=900. U0 = U0 + U90;
P - степень поляризации излучения;
t - азимут поляризации излучения.
Вектор-параметр Стокса для яркости излучения объекта в таком случае будет следующим:
1
Uj (N, L) = [ W(,T,y,z) / ] | P cos 2 t | , ( 62 )
| P sin 2 t |
0
В свою очередь, влияние поляризационного фильтра на излучение от объекта описывается матрицей Мюллера:
1 cos 2sin20
ij = п | cos 2 cos2 2 sin2 cos20 | , ( 63 )
| sin 2 cos 2 sin 2 sin220 |
00 0 0
где п - энергетический коэффициент пропускания фильтра;
- азимут поляризации фильтра, отсчитываемой относительно плоскости референции.
Тогда при положении поляризационного фильтра с азимутами
= 00 и = 450, матрицы ij будут иметь вид:
1 1 0 0
ij(0) = п | 1 1 0 0 | ; ( 64 )
| 0 0 00 |
00 0 0
1 0 1 0
ij(45) = п | 0 0 0 0 | ; ( 65 )
| 1 0 10 |
00 0 0
Вектор-параметр Стокса для энергетической яркости излучения, прошедшего произвольный поляризационный фильтр, можно записать:
4
Li(,T,P) = ij Lj(,T,P). ( 66 )
j =1
Сигнал на выходе приёмника излучения запишется в виде:
2
U1 = c( Li(,T,P) d, ( 67 )
1
где c( = cos dS S 0 a
Тогда, вектор-параметры Стокса для яркости излучения, прошедшего поляризационный фильтр при азимутах поляризации =00 и =450, будут следующие:
1 + P cos 2 t
Li(0) = п W(,T,y,z) / ] | 1 + P cos 2 t | , ( 68 )
| 0 |
0
1 + P sin 2 t
Li(45) = п W(,T,y,z) / ] | 0 | , ( 69 )
| 1 + P sin 2 t |
0
Как известно, первая строка вектор-параметра Стокса характеризует энергетические характеристики излучения, поэтому выражение для сигналов приёмника при двух положениях поляризационного фильтра можно записать в виде:
U1 =п(1+Pcos2t)[(1/ cosdS ]Sa(W(,T,y,z)d
( 70 ).
U2 =п(1+Psin2t)[(1/ cosdS ]Sa(W(,T,y,z)d
Если обозначить одинаковые множители U1 и U2 в виде:
B( T )=п[(1/ cosdS ]Sa(W(,T,y,z)d
то формулы ( 70 ) примут вид:
U1 =B( T )( 1 + P cos2t )
( 71 )
U2 =B( T )( 1 + P sin2t ).
Упростим формулы ( 71 ), пронормировав их B( T ):
U1н = 1 + P cos2t ;
( 72 )
U2н = 1 + P sin2t .
Если эти выражения записать соответственно обозначениям, принятым в разделах 2.1 и 2.2 для каждого элемента разложения кадра, то получим:
U1н(N, L) = 1 + P(N, L) cos2t ;
( 73 )
U2н(N, L) = 1 + P(N, L) sin2t .
Как уже отмечалось, для формирования модели изображения очень важной является задача преобразования объёмного объекта в плоский кадр. В разделах 2.3 и 2.4 это преобразование проведено через координаты и сферической системы координат, связанные с элементами разложения кадра по строкам L и элементами в строках N на основе геометрических и логических связей.
Для модифицированной модели можно использовать в качестве системы преобразования объёмного объекта в плоский кадр декартову систему координат.
Пусть тепловизионная система с поляризационной насадкой строит кадр с изображением объекта, которому соответствует в плоскости объекта поле зрения с размерами по вертикали - ( zk -zн), а по горизонтали - ( yk - yн ), где zk, zн, yk, yн - конечные и начальные координаты поля зрения в системе координат объекта
( рис.6 ). Причём, расположение по вертикали вдоль оси OZ и по горизонтали вдоль оси OX, что соответствует геометрии наблюдения объекта по рисунку 6. Здесь показан общий принцип системы построения кадра. Более конкретное преобразование объекта в кадр будет представлено для каждого объекта конкретно.
2.7. Оптико-математическая модель поляризационных
изображений с учётом эллиптичности поляризации
теплового излучения.
Все выводы, приведённые выше, представлены для частично линейно-поляризованного излучения. В случае эллиптично-поляризованного излучения окончательные формулы будут иметь иной вид.
Если обозначить эллиптичность через, то для линейно-поляризованного излучения =0, а для эллиптично-поляризованного имеет значения, которые отличны от нуля. Поэтому вектор-параметр Стокса для эллиптично-поляризованного излучения в нормированном виде представляются как:
1
U(N, L) = U0 | P(N, L) cos2tcos2 | . ( 74 )
| P(N, L) sin2tcos2 |
P(N, L) sin2
После вывода, аналогичного случаю линейно-поляризованного излучения можно получить выражения для нормированных видеосигналов U1н и U2н для эллиптично-поляризованного излучения:
U1н = 1 + P cos2cos2t ;
U2н = 1 + P cos2sin2t . ( 75 )
Таким образом, на основании формул ( 75 ) и ( 16 ) - ( 19 ) можно формировать оптико-математические модели поляризационных тепловизионных изображений объектов с учётом эллиптичности поляризации излучения.
2.8. Модифицированная формула моделирования
изображения диска, сферы и эллипсоида.
Модифицированная модель изображения предполагает иное, чем в разделе 2.3 преобразование объёмного тела, наблюдаемого тепловизионной системой, в плоский кадр с изображением этого объекта, поэтому для пояснения процесса преобразования воспользуемся рис.6. Из данного рисунка видно, что ( z0, y0 ) - это координаты центра объекта и кадра, R - радиус самой сферы, а rt - радиус-вектор плоского кадра, связывающий координаты вертикали Z и координаты горизонтали Y с центром элемента разложения кадра. Из геометрических связей можно определить rt :
.
rt = ( y-y0)2 + ( z-z0)2 . ( 76 )
Элементу разложения кадра dS’ с координатами ( y, z ) будет соответствовать элементарная площадка поверхности сферы с координатой Х. Поскольку уравнение сферы в декартовой системе координат, с геометрией наблюдения, аналогичной рис. 7, имеет вид:
x2 ( y-y0)2 ( z-z0)2
f( x, y, z ) = ---- + --------- + --------- = 1, ( 77 )
R2 R2 R2
то координата Х элементарной площадки dS определяется по формуле:
.
x = R2-( y-y0)2 + ( z-z0)2 . ( 78 )















