Реферат: Методика моделирования тепловизионных изображений
--PAGE_BREAK--rн определяется как скалярное произведение этих векторов, получим:P(y) = [ 1-(n*rн) ] ×a . ( 13 )
Подставив это выражение в формулу (12) получим:
e÷÷ (y) 1+ [ 1 — (n*rн)]×a
--------- = ------------------------- . ( 14 )
eûë (y) 1 — [ 1 — (n*rн)]×a
Тогда, с учётом соотношения (12), из формулы (11) получим основное уравнение, выражающее зависимость между степенью поляризации P’(N, L) и формой объекта через функцию распределения нормали n для каждого элемента поверхности объекта:
1+ [ 1 — (n*rн)]×a
------------------------ [(n*j)2 — (n*k)2] +[(eûë*j)2 — (eûë*k)2]
1-[ 1 — (n*rн)]×a
P’(N, L) = ---------------------------------------------------------------------- . ( 15 )
1+ [ 1 — (n*rн)]×a
— [(n*j)2 — (n*k)2] +[(eûë*j)2 + (eûë*k)2]
1-[1 — (n*rн)]×a
С помощью этой формулы можно определить степень поляризации всех элементов наблюдаемой тепловизором части поверхности объекта любой формы. Для этого нужно знать направление нормали n для каждого элемента поверхности в зависимости от его положения в декартовой системе координат. Оно определяется как оператор Гамильтона ( набла-оператор ) от функции f(x,y,z) = 0, описывающий форму объекта:
[( df/dx ) ×i + ( df/dy ) ×j + ( df/dz )×k ]
n= ---------------------------------------------------- . ( 16 )
[( df/dx )2+ ( df/dy )2 + ( df/dz )2] 1/2
Единичный вектор наблюдения rн определяется как разница векторов l и R по формуле:
rн = ( l — R ) / |( l — R )|, ( 17 )
где l — вектор, определяющий положение декартовой системы координат по отношению к точке наблюдения H;
R — радиус-вектор элемента dS поверхности объекта, определяющий его положение в декартовой системе координат x, y, z с единичными ортами i, j, k.
Радиус-вектор задаётся R формулой :
R = x ×i+ y ×j + z ×k . ( 18 )
Если направление наблюдения центра декартовой системы координат выбрано вдоль оси х, то есть направление вектора l и оси х совпадают, то вектор l выразится в виде:
l = l×i , ( 19 )
где l — расстояние от центра декартовой системы координат О до точки наблюдения Н;
i — единичный орт оси ОХ .
В этом случае выражение (17) примет вид:
rн = [( l-x)i + y×j +z×k] / [( l-x)2+ y2+ z2]1/2. ( 20 )
Вектор перпендикулярной составляющей коэффициента излучения eûëперпендикулярен плоскости, определяемой векторами n и rн ( плоскости наблюдения ), и находится как векторное произведение этих векторов по формуле:
eûë = [ n* rн ] /
|
[ n* rн ]
|. ( 21 )
Таким образом, определив степень поляризации P’ от всех элементов видимой части объекта, можно построить оптико-математическую модель поляризационных тепловизионных изображений объектов любой формы.
2.1. Теория моделирования поляризационных тепловизионных
изображений на основе степени и азимута поляризации
теплового изображения.
Для описания этого метода воспользуемся рис. 3.
Допустим, что азимут поляризации излучения элемента dS поверхности объекта составляет угол t с поверхностью референции.
Для определения степени поляризации P’ необходимо найти величины видеосигналов U0и U90 поляризационных тепловизионных изображений элементов dS поверхности объекта при азимутах поляризатора t=00и t=900. Выразим U0и U90 через параллельную и перпендикулярную составляющие коэффициента излучения элемента dS и азимут t поляризации этого элемента, который представляет собой угол между плоскостью поляризации ( ось ОА ) и плоскостью референции ( ось OY ). В общем случае, когда азимут t поляризации излучения элемента dS не совпадает с азимутом поляризатора, обе компоненты коэффициента излучения дают вклады в величины видеосигналов U0и U90 следующим образом:
U0(N, L) = Umax×cos2 t + Umin×sin2 t = A(N, L)× ( e÷÷ ×cos2 t + eûë×sin2 t) ; ( 22 )
U90(N, L) = Umax×sin2 t + Umin×cos2 t = A(N, L)× ( e÷÷ ×sin2 t + eûë×cos2t) ; ( 23 )
где Umax= A(N, L)× e÷÷, Umin= A(N, L)× eûë.
Согласно формуле (6) найдем степень поляризации P’(N, L) излучения элемента dS объекта в виде:
P’(N, L) = [e÷÷ — eûë ] / [e÷÷ + eûë]×cos(2×t) = P ×cos(2×t) , ( 24 )
где P = [e÷÷ — eûë ] / [e÷÷ + eûë ] — распределение степени поляризации излучения элементов dS объекта.
Так как cosy = ( n* rн ), то с учётом формулы (12) имеем:
P’(N, L) = [ 1- ( n* rн ) ]× а×cos(2×t); ( 25 )
В связи с тем, что вдоль оси ОА расположен вектор nyz, являющийся проекцией вектораn на плоскость xyz, то справедливо выражение:
cos t = ( nyz*j ) , ( 26 )
тогда, приняв во внимание тождество
cos(2×t) = 2×cos2t — 1,
выражение (25) для расчёта степени поляризации всех элементов поверхности объекта примет вид:
P’(N, L) = а×[ 1- ( n* rн ) ]× [ 2×( nyz*j )2 -1 ]. ( 27 )
Таким образом, формулы (15) и (27) с учётом формул (16) — (21) являются оптико-математической моделью поляризационных тепловизионных изображений излучающих объектов [5,6]. В тех случаях, когда необходимо моделировать поляризационные тепловизионные изображения по распределению степени поляризации, можно воспользоваться выражением:
P(N, L) = а×[ 1- ( n* rн ) ]. ( 28 )
2.3. Формулы для моделирования изображения
диска, сферы и эллипсоида.
Для подтверждения теории моделирования поляризационных тепловизионных изображений рассмотрим объекты в виде сферы, эллипсоида и диска. Как уже отмечалось раньше, традиционный тепловизионный метод при наблюдении этих объектов сверху даёт одинаковое изображение как по контуру, так и внутри контура, несмотря на явное различие формы этих объектов внутри контура изображения видимой части их поверхности. Для подробного вывода остановимся на сфере, как наиболее наглядном и симметричном объекта ( рис. 4).
Уравнение сферы в декартовых координатах имеет вид:
f(x,y,z) =x2+ y2+ z2 — R2= 0. ( 29 )
Тогда n = (x ×i+ y ×j + z ×k) /R — вектор нормали сферы,
где R = (x2+ y2+ z2)1/2 — радиус сферы.
Вектор наблюдения rн можно определить из формулы (17):
rн= [( l-x)×i— y×j — z×k] / [R2+ l2+ 2×l ×x]1/2. ( 30 )
Тогда по правилам векторного умножения:
e = [ n* rн] = ( ny ×rнz — nz ×rнy)×i +( nz ×rнx — nx ×rнz)×j +( nx ×rнy — ny ×rнx)×k ;
в нормированном виде:
_____________
e
ûë= ( lz×i — ly×j ) / (R × Ö R2+ l2 — 2 ×l ×x ), ( 32 )
Теперь определим все остальные недостающие выражения для формулы (15):
_____________
( n* rн ) = (x ×l -R2)/ (R × Ö R2+ l2 — 2 ×l ×x ), ( 33 )
( n* j)2 = y2 / R2 ; ( 34 )
( n* k)2 = z2 / R2 ; ( 35 )
(
eûë * j)2 = l2×z2/ (R2× ( R2+ l2 — 2 ×l ×x ); ( 36 )
(
e÷÷* k)2 = l2×z2/ (R2× ( R2+ l2 — 2 ×l ×x ); ( 37 )
После подстановки формул (30) — (37) в выражение (15), получим:
l ×x — R2
2 — ---------------------------------
R2× ( R2+ l2 — 2 ×l ×x )1/2 æ y2 — z2 ö é l2 ×z2 — l2 ×y2 ù
----------------------------------------- × ï---------
ê
+ ï--------------------------- ç
l ×x — R2 è R2 ø ëR2×( R2+ l2 — 2 ×l ×x )û
---------------------------------
R2× ( R2+ l2 — 2 ×l ×x )1/2
P’ (N, L) = ---------------------------------------------------------------------------------------------- .
l ×x — R2
2 — ---------------------------------
R2× ( R2+ l2 — 2 ×l ×x )1/2 æ y2+ z2 ö é l2 ×z2 + l2 ×y2 ù
----------------------------------------- × ï продолжение
--PAGE_BREAK-----------
ê
— ï--------------------------- ç
l ×x — R2 è R2 ø ëR2×( R2+ l2 — 2 ×l ×x)û
---------------------------------
R2× ( R2+ l2 — 2 ×l ×x )1/2
После упрощения это выражение принимает вид:
P’(N, L) = [( y2 — z2 ) / ( y2 + z2 )] ×( 1 — x/R ). ( 38 )
Это есть степень поляризации теплового изображения сферы в декартовых координатах.
Перейдем к сферическим координатам:
X = R×sinq ×cosj;
Y = R×sinq ×cosj;
Z = R ×cosq.
Тогда выражение (38) принимает вид:
sin2q ×sin2j— cos2q
P’(N, L) = — ( 1 — sinq×cosj) . ( 39 )
sin2q×sin2j +cos2q
Это и есть степень поляризации теплового изображения сферы в сферических координатах.
Аналогично можно получить формулы для эллипсоида. Для этого необходимо начать вывод с функции:
f(x,y,z) =x2 / b2+ y2 / a2+ z2 / c2 — 1= 0. ( 40 )
С учётом обозначения K = b/a — коэффициента сжатия эллипсоида ( b — большая полуось эллипсоида, a — малая ), получим формулу для степени поляризации в декартовых координатах:
________________
P’(N, L) = [( y2 — z2) / ( y2 + z2)] ×[ 1 — ( x / Öx2 + k2×y2 + k2×z2)] . ( 41 )
C учётом сферических координат для эллипсоида:
X = b×sinq ×cosj;
Y = a×sinq ×cosj;
Z = a ×cosq.
степень поляризации принимает вид:
sin2q ×sin2j— cos2q é sinq×cosj ù
P’(N, L) = — × ê 1- ------------------------------------------------------ ç(42)
sin2q×sin2j +cos2q ë Ö sin2q×cos2j +k2×(sin2q×sin2j +cos 2q) û
Что касается диска, то для него используется формула ( 42 ), с учётом, что коэффициент сжатия k := 0.1, т.е. эллипсоид сжатый до состояния диска, когда большая полуось составляет всего лишь 10-ю часть от малой полуоси; для сферы формула ( 42 ) справедлива при k = 1. Таким образом, для получения модели поляризационного тепловизионного изображения диска, сферы и эллипсоида можно пользоваться формулой ( 42 ) с использованием различных значений k. При этом необходима связь углов qи jс номерами строк L и номерами элементов в строках N тепловизионного кадра. На основе геометрии наблюдения и логических рассуждений были получены следующие связи:
q= L× p /L0 ; ( 43 )
j= ( N × p/ N0) — p/2 ; ( 44 )
где L0 — число всех строк в кадре;
N0 — число элементов в каждой строке.
2.4. Формула моделирования изображений конуса.
Вывод формулы моделирования изображений конуса аналогичен выводу формулы для тел типа эллипсоида, но для разнообразия расположим конус по другой оси координат — вдоль оси OZ ( рис. 5).
В декартовой системе координат уравнение конуса имеет вид:
f(x,y,z) = x2 / a2+ y2 / a2 — z2 / c2 = 0. ( 45 )
где а — радиус основания конуса;
с — высота конуса.
Вектор нормали n в соответствии с формулой (16), имеет вид:
[(-2×z/c2)×k+ (2×x/a2)×i+ (2×y/a2)×j]
n = ------------------------------------------------- . ( 46 )
Ö (2×x / a2 )2+(2×y / a2 )2+(2×z / c2 )2
В свою очередь вектор наблюдения для конуса данного расположения в декартовой системе координат имеет вид:
.
rн = — xн ×i — yн ×j — ( l — zн )×k / Ö x2н + y2н + ( 1 — z2н) , ( 47 )
Если конус наблюдается из бесконечности, то упрощение в формулах можно произвести в процессе вывода, а не в окончательном виде, как в случае эллипсоида. Так, при l стремящемся в бесконечность, rн = -k.
Тогда произведение ( n* rн) принимает вид:
.
( n* rн )= (2×z/c2) / Ö(2×x / a2 )2+(2×y / a2 )2+(2×z / c2 )2 ( 48 )
принимая во внимание то, что коэффициент сжатия конуса k = c / a, тогда
.
( n* rн )= z / Ö(x2 + y2 )×k4+z2 . ( 49 )
Если применить способ формирования изображения на основе степени и азимута поляризации, то необходимо для конечной формулы пользоваться формулой ( 27 ), которая для случая наблюдения объекта вдоль оси OZ примет вид:
P(N, L) = [ 1 — ( n×rн )] × [ 2 ×( nxy×i )2 — 1 ]. ( 50 )
в этом случае
.
nxy = (x×i+ y×j) / Öx2 + y2 ; ( 51 )
.
( nxy×i ) = x / Öx2 + y2 ; ( 52 )
Соединив формулы ( 49 ) — ( 51 ), получим степень поляризации в виде:
.
P’(N, L) = [ 1 — z / Ö(x2 + y2 )×k4+z2 ] × [ 2 ×x2 / (x2 + y2) — 1 ] . ( 53 )
Для удобства вывода выражения для P’(N, L) в сферической системе координат, воспользуемся переводом компонент в другую систему координат:
X = sinq ×cosj;
Y = sinq ×cosj; ( 54 )
Z = cosq.
.
( n* rн )= 1 / Ö1 + k4 tg2q .
Поскольку угол q для конуса является постоянным и равным отношению радиуса основания к высоте, то справедливо выражение:
tgq= a / c. ( 55 )
Подставив формулу ( 55 ) в выражения ( 54 ), получим :
.
( n* rн )=1 / Ö1 + k2, ( 56 )
( nxy×i ) = cosj. ( 57 )
Тогда
.
P’(N, L) = [ 1 — 1 / Ö1 + k2 ]×cosj. ( 58 )
Таким образом, формула ( 58 ) является оптико-математической моделью поляризационного тепловизионного изображения конуса. При этом угол jсвязан с номером строки L и с номером элемента строки N кадра зависимостью:
j= 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. Формировать оптико-математическую модель изображения тепловизионной системы с поляризационным фильтром можно модифицированным методом моделирования не основе вектор-параметра Стокса и влияния на излучение от объекта поляризационного фильтра. Причём исходным выражением для видеосигнала будем считать:
l2
( N, L ) = ( 1/ p)×w ×cosy(N,L)×dS(N,L)×òSl×W(l,T,y,z)×t(l)×ta(l)×dl ( 60 );
l1
Вектор-параметр Стокса, описанный в разделе 2.1 формулой ( 4 ), в нормированном виде выглядит следующим образом:
é 1 ù
Uj (N, L) = U0× | P × cos 2 × t | продолжение
--PAGE_BREAK--( 61 );
| P × sin 2 × t |
ë 0 û
где U0 — суммарный видеосигнал при азимутах поляризации излучения t=00и t=900. U0= U0+ U90;
P — степень поляризации излучения;
t - азимут поляризации излучения.
Вектор-параметр Стокса для яркости излучения объекта в таком случае будет следующим:
é 1 ù
Uj (N, L) = [ W(l,T,y,z) / p ] × | P × cos 2 × t | , ( 62 )
| P × sin 2 × t |
ë 0 û
В свою очередь, влияние поляризационного фильтра на излучение от объекта описывается матрицей Мюллера:
é 1 cos 2× d sin2× d 0 ù
tij= tп× |cos 2× d cos2 2× d sin2× d× cos2× d 0 | , ( 63 )
|sin 2× d cos 2× d× sin 2× d sin22× d 0 |
ë 0 0 0 û
где tп— энергетический коэффициент пропускания фильтра;
d— азимут поляризации фильтра, отсчитываемой относительно плоскости референции.
Тогда при положении поляризационного фильтра с азимутами
d= 00 и d= 450, матрицы tijбудут иметь вид:
é 1 1 0 0 ù
tij(0)= tп× | 1 1 0 0 | ; ( 64 )
| 0 0 0 0 |
ë 0 0 0 û
é 1 0 1 0 ù
tij(45)=tп× | 0 0 0 0 | ; ( 65 )
| 1 0 1 0 |
ë 0 0 0 û
Вектор-параметр Стокса для энергетической яркости излучения, прошедшего произвольный поляризационный фильтр, можно записать:
4
Li(l,T,P) = S tij ×Lj(l,T,P). ( 66 )
j =1
Сигнал на выходе приёмника излучения запишется в виде:
l2
U1 = òc(l) ×Li(l,T,P) ×dl, ( 67 )
l1
где c(l) = w ×cosy×dS ×Sl× t(l) × ta(l) .
Тогда, вектор-параметры Стокса для яркости излучения, прошедшего поляризационный фильтр при азимутах поляризации d=00и d=450, будут следующие:
é1 + P × cos 2 × tù
Li(0) = tп×W(l,T,y,z) / p ] × |1 + P × cos 2 × t | , ( 68 )
| 0 |
ë 0 û
é1 + P × sin 2 × tù
Li(45) = tп×W(l,T,y,z) / p ] ×| 0 | , ( 69 )
| 1 + P × sin 2 × t |
ë 0 û
Как известно, первая строка вектор-параметра Стокса характеризует энергетические характеристики излучения, поэтому выражение для сигналов приёмника при двух положениях поляризационного фильтра можно записать в виде:
l2
U1 =tп×(1+P×cos2×t)×[(1/ p)×w ×cosy×dS ]×òSl×t(l)×ta(l)×W(l,T,y,z)×dl
l1
( 70 ).
l2
U2 =tп×(1+P×sin2×t)×[(1/ p)×w ×cosy×dS ]×òSl×t(l)×ta(l)×W(l,T,y,z)×dl
l1
Если обозначить одинаковые множители U1 и U2 в виде:
l2
B( T )=tп×[(1/ p)×w ×cosy×dS ]×òSl×t(l)×ta(l)×W(l,T,y,z)×dl
l1
то формулы ( 70 ) примут вид:
U1 =B( T )× ( 1 + P × cos2×t )
продолжение
--PAGE_BREAK--( 71 )
U2 =B( T )× ( 1 + P × sin2×t ).
Упростим формулы ( 71 ), пронормировав их B( T ):
U1н = 1 + P × cos2×t ;
( 72 )
U2н = 1 + P × sin2×t .
Если эти выражения записать соответственно обозначениям, принятым в разделах 2.1 и 2.2 для каждого элемента разложения кадра, то получим:
U1н(N, L) = 1 + P(N, L) × cos2×t ;
( 73 )
U2н(N, L) = 1 + P(N, L) × sin2×t .
Как уже отмечалось, для формирования модели изображения очень важной является задача преобразования объёмного объекта в плоский кадр. В разделах 2.3 и 2.4 это преобразование проведено через координаты q и jсферической системы координат, связанные с элементами разложения кадра по строкам L и элементами в строках N на основе геометрических и логических связей.
Для модифицированной модели можно использовать в качестве системы преобразования объёмного объекта в плоский кадр декартову систему координат.
Пусть тепловизионная система с поляризационной насадкой строит кадр с изображением объекта, которому соответствует в плоскости объекта поле зрения с размерами по вертикали — ( zk -zн), а по горизонтали — ( yk — yн ), где zk, zн, yk, yн — конечные и начальные координаты поля зрения в системе координат объекта
( рис.6 ). Причём, расположение по вертикали вдоль оси OZ и по горизонтали вдоль оси OX, что соответствует геометрии наблюдения объекта по рисунку 6. Здесь показан общий принцип системы построения кадра. Более конкретное преобразование объекта в кадр будет представлено для каждого объекта конкретно.
2.7. Оптико-математическая модель поляризационных
изображений с учётом эллиптичности поляризации
теплового излучения.
Все выводы, приведённые выше, представлены для частично линейно-поляризованного излучения. В случае эллиптично-поляризованного излучения окончательные формулы будут иметь иной вид.
Если обозначить эллиптичность черезg, то для линейно-поляризованного излучения g=0, а для эллиптично-поляризованного g имеет значения, которые отличны от нуля. Поэтому вектор-параметр Стокса для эллиптично-поляризованного излучения в нормированном виде представляются как:
é 1 ù
U(N, L) = U0 | P(N, L) × cos2×t× cos2×g | . ( 74 )
| P(N, L) × sin2×t× cos2×g |
ë P(N, L) × sin2×g û
После вывода, аналогичного случаю линейно-поляризованного излучения можно получить выражения для нормированных видеосигналов U1н и U2н для эллиптично-поляризованного излучения:
U1н = 1 + P × cos2×g × cos2×t ;
U2н = 1 + P × cos2×g × sin2×t . ( 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 )
Подстановкой формулы ( 76 ) в выражение ( 78 ) преобразуем выражение для Х:
.
x =ÖR2 — rt2 . ( 79 )
Далее для определения степени и степени поляризации воспользуемся формулами ( 16 ) — ( 19 ) применительно к конкретному объекту:
df df df 2 ×x 2 ×(y-y0) 2 ×(z-z0)
n= -----×i+ — ×j+ — ×k= — ×i+ — ×j+ -----------×k; ( 80 )
dx dy dz R2 R2 R2
2 ×(y-y0) 2 ×(z-z0)
nyz= — ×j+ -----------×k; ( 81 )
R2 R2
é (y-y0) ù y-y0
t = arccos | ------------------- | = arccos ------------ ; ( 82 )
ë Ö(y-y0) + (z-z0) û rt
( n* rн ) .
cosy= — = x / Öx2 + rt2 . ( 83 )
|( n* rн ) |
По формуле ( 12 ) можно найти Р:
ì x ü
Р = a ( 1- cosy) = a × | 1 - ---------- | . ( 84 )
î Öx2 + rt2þ
Для получения оптико-математической модели достаточно подставить формулы ( 80 ) — ( 84 ) в формулы для видеосигналов ( 73 ) ( или ( 74 ) для случая эллиптично-поляризованного излучения ).
Вернёмся теперь к формуле ( 82 ) для азимута поляризации излучения. Как видно из этой формулы, t зависит только от y и z, а от координаты х зависимости нет. Поскольку в данной работе рассматриваются объекты, различающиеся по форме именно вдоль оси Х, а в плоскости осей Y и Z ( т.е. в кадре ) имеющие одинаковый контур, то можно сделать вывод, что значение азимута поляризации t для всех рассматриваемых здесь объектов ( конус, эллипсоид, сфера ) будет одинаковым.
Для полной ясности необходимо установить распределение азимута поляризации по поверхности этих фигур. По формуле ( 82 ) рассмотрим некоторые конкретные случаи. Например, при z=z0и y>y0, t=0 при z=z0и y< y0, t = p;при y=y0и z> z0, t= -p/2; при y=y0, z<z0, t= p/2.
Если попробовать свести эти результаты к схематичному распределению азимута поляризации излучения внутри контура с учётом того, что в случаях, не указанных в примере, азимут поляризации принимает промежуточные положения, то получается рисунок 7.
Чтобы сформировать оптико-математическую модель для эллипсоида, воспользуемся рисунком 8 и уравнением эллипсоида в декартовой системе координат:
x2 ( y-y0)2 ( z-z0)2
f( x, y, z ) = ---- + --------- + --------- = 1, ( 85 )
a2 b2 c2
При моделировании для упрощения примем:
b = c = R ; ( 86 )
a = k ×R, ( 87 )
где k — коэффициент сжатия.
Тогда уравнение ( 85 ) примет вид:
x2 ( y-y0)2 ( z-z0)2
f( x, y, z ) = -------- + --------- + --------- = 1, ( 88 )
k2 ×R2 R2 R2
Уравнение для координаты х, исходя из выражения ( 88 ), будет следующим:
.
x = k × ÖR2 + rt2 . ( 89 )
Выражение для азимута поляризации в случае объекта типа эллипсоида, будет таким же, как для случая со сферой ( 88 ), так как азимут поляризации не зависит от координаты х:
t = arccos [(y — y0) / rt ]. ( 90 )
Степень поляризации для каждого элемента разложения кадра с координатами ( у, z ) можно определить аналогично сфере из формул ( 16 ) — ( 19 ) и ( 25 ) — ( 27 ):
.
cosy= x / Öx2 + k4×( y-y0 )2+ k4×( z-z0 )2= x / Öx2 + k4×rt2 . ( 91 )
Степень поляризации, соответственно, равняется
.
P = a ×( 1 — x / Öx2 + k4×rt2 ) . ( 92 )
Далее, по выражения ( 73 ) ( или ( 75 ), в случае эллиптичной поляризации ) можно получить модели изображений эллипсоида при азимутах фильтра d = 00и d= 450соответственно. Причём, при к = 1 формулы для эллипсоида становятся аналогичными для сферы. Если в формулы ( 73 ) или ( 75 ) подставить к = 0.1, то это будет модель изображения диска. Во всех остальных случаях можно получить модели изображений эллипсоида с различными коэффициентами сжатия.
2.9. Модифицированная формула моделирования
изображения конуса.
Рассмотрим, согласно рис. 9, уравнение конуса в декартовых координатах:
f(x, y, z) = — ( h- x )2 / h2 + ( y — y0 )2 / R2 + ( z — z0)2 / R2 = 1, ( 93 )
где R — радиус основания конуса;
h — высота конуса.
Уравнение для координаты х в случае конуса будет иметь вид:
x = h × ( 1 — rt / R) . ( 94 )
Значение степени поляризации определим аналогичным образом. Для этого найдём вектора n и rt :
n = — 2 ×( h — x ) × i/ h2 + 2 ×( y — y0) × j/ R2 + 2 ×( z — z0) ×k / R2, rн = i. ( 95 )
Тогда
продолжение
--PAGE_BREAK--
еще рефераты
Еще работы по производству
Реферат по производству
Электропривод грузового подъёмника
3 Сентября 2013
Реферат по производству
Температурный расчет сварки
3 Сентября 2013
Реферат по производству
Теория развития вселенной
3 Сентября 2013
Реферат по производству
Классическая физика самоорганизующиеся системы и микромир
3 Сентября 2013