книги2 / 424
.pdf. 3.10. @ # # + " " " # ##:
– # + " #" # /;– # + " " "
( N + );
– (N + ) " + #" /
( ( +-(()7 -3 ; ( 5
# & $* ! !!& ! 4 ! 4 , & # ! # *. & ! !! ! ! - # ! ;, # ! ! !!&, & # ! - , « ! » ( , , U , !. .) #! ! & # &$ # ! $. # !! # ! !!& & & 4 ! & 4 , ! 4 ! 4 &$ !!& MATLAB (4 & ' ! «mat»), & # ! - # ' & #! % 4 &$. < #! - # &$ # % ( ! &$ ! ) ! % & # , # * *+ 4 ! &: 320×240, 640×480, 128×128, 512×512. =-! & ! # ! # $ # % 4 4 4 ! JPAG ( ! JPEG),
& ! |
4 4 ! ( # % ) ! - |
% & |
! , U $! & - |
! !!& ! + * # !! . K ! # # * $ 4 & ! !!
4 '-, *+ #V ! . =- $ # % + % ! ( & # # &$ $,$ *+ V ) ! + * ! - , & ! % & # (! % -
71
# ! ! , * ! - & $ V ). & # & #- * $ ! !!& % ! *# ! 4 PCM=A # #-(~ 6 ? ). ! % & - * & #! *+$ , ' # ! * & &$&. U ! ! & # % -! % 30…50 ? &' , 4 ! ! -. & 4 # - ! !! ! ! # USB ! * FireWire.
-4 3 3 +-(()7 -3 ; ( 5
; &! ! # &' - ' -' !, . . #& ! « » 4 .!& # # * * -# ! & $ ! #- % !, & ! ' ! -' !.
, ! # * U-! &$ $, &$ # # - &! ! ! 4 4 . < ! ! ! - &$ ! , % *+ # #, & % 5 4 ( . 3.11). # - # ! ! ! , % &! # .
|
|
|
. 3.11. 0 # ##/ " ( # ###" #" # ):
– # ## #" ;
– #" # ;
– #" 8 (0) 8 (?0)
. ! « ! # » # ! $ ! !! * U44 - ! . < # * ! # :
72
1)! # $ # % & # % , #& * « ! # *+ !»;
2)$! (3D) ! # , !$ , # &$ U ! - , – [5].
< ! ! "-# % T(i, j, ) , # -
&$ ! ( h ).!$ $ ! , - 4 & * $% , . . m h . R ! ! # # % ! *+ :
Cnorm |
(i, j, ) |
T(i, j, ) Tamb (i, |
j) |
, |
||
Tnorm (i, j, norm ) |
Tamb (i, j) |
|||||
|
|
|
Cnorm (i, j, ) – ! # & ; Tamb (i, j) – ! - % *+ &, + ! ! % #! ; Tnorm (i, j, norm ) – # *$ # % . , ! # &4 ' ! % ! # &- &! ! ! , U ! * #& !& «!- & » ' !&, ! &$ # & ! -. A ! # *+ # % & ! # !,$& & $ 4 + #- , ! ! # & # % % & -! & 4 ! &! ! 4 -! ! & ! , # m . ! , ! ! -! # *+ & * ! !! . =$ - ! ! # # % # . 3.12.
! ! # ! !!,! ! ! # , . 3.13.$ ! !! #& - ! & ! ! ! , # &$ ! + & #-! ( . 3.13, ). ! & !$ ! ! ! m . !- # ! !! ( h ). U ! !- # & , h , &! ( !. %# % ), ! ! # ! !! ,
m , #& 4 & * $, #% ! $, ! &! !.
73
. 3.12. % # # + "
|
|
. 3.13. % # # ##:
– # ## ;
– # / # ##
?-. K #" 8 # 8 (pulse phase thermography), % & K !, K - * , ! # R 4 #! -! & ! Tij ( ) [5, 18].
"! 4 # ! 4 + / # $ ! ó # ! 4 &$ 4 &$ #$. "# , + & $- * ! U * 4 !& ( -! ! &! ! ! !- ! ) $ # ! , & *, U ( ! ) # # &-! ( 4 # ). ; & & ! & * % &! # !, # #-# 4 . k ! % -&$ , ! + *+$ ! % $! # 4 . C ! && &$ ! ! !-& ! ' . ! # !, # & *
74
# &$ &$ 4 . U # , -$ # ! % # &- &! ! , # ! &$ , & 4 & & # &! - ! . U ! 4 * 4 ! * % 4 #& ! % # % *+ *+ .
R ! # !! 4 4 $ ! % **+ ! % !. "# , : 1) 4 & $ -
# * |
&! ! &! ! ; 2) |
|
R ! &$ ! &$ ! # & !- |
||
& ' |
! + , ! & & . = ! , |
|
% |
$ 4 * #$ , |
|
&' # |
# , ! # ! &! !- |
! ( ! ! ) ! % & &$ 4 # . ; , #! % & , # & $ ! +- * , &! ! , !&! -! # , ' ! ! . .
; R -# , # ! , &$! &$ # % # * « ! - », ! ! ó! # * .
! # !, R -# # * ! $ "-# % . ;# % 4 #& (8 ##/) !- & (# ##/) R - !, #- & # R # & 4-! &$ &$, U$ . ! - !&! ! * #! % + &$ - &$ # ! , & # &, ! , + -* ! 4 * ! # 4 #.
< ! ! + * R - #! ó ! !!. " # ! #-
|
R |
Fs , |
4 # : |
|
||
|
|
|
|
1 |
N |
|
|
|
Fs |
|
" fn e2 j (n 1) (s 1) / N , |
||
|
|
|
||||
|
|
|
|
N n 1 |
|
|
fn – $ 4 , % + N ( N - |
||||||
! !! |
); s – |
; j – ! ! ; |
||||
n – ! . |
|
|
|
|
||
|
T # % . 3.14, . K M - |
|||||
! (4 # ) 7 # R |
|
|||||
& . 3.14, , . |
|
|
||||
|
|
|
75 |
|
a |
|
. 3.14. " + 8 # " #" / # ( ) " # / !-: # ( ) 8 ( )
$ , 4 ! 4 #& - # R % 20 , ! + 1 -. U44 % (aliasing), # * 4 !- * &, *+ ! 1 10 ( 4 ! + 11 20 * # &! % ! &$ 10- ).
% , $ ! , - ! , ! ! * ! - ! 4 , ! 4 # ! 4 4 !-! , # # # &, $!& % % +$ 4 . = -&, ! ! % ! ! !! ! % & ! , & & & R - , $ - !& % $&$ 4 ( & U. 3.14, ). U ! ! 4 # ! 4 & ! -
! , + ! # $ # % .
+ ! ! ! + ! 4 # ! 4 & ! & ! ! , -! &! 4 ! & , - # * $ ! .
76
|
|
. 3.15. # !- " :
– 8 / 3 + ;
– # #" / + /
4 : . <= (Principal Component Analysis) – ! - , ' ! *+ ! ! &$ &$ ( . 3.16) [19]. $ &$ ! % % , ! , 4 ! ! 4 !& ( # &, # , . .). U !&$ ( ! ! ! ) -! ' !&$ , #! &$ -! , &! ' * .
|
. 1 |
. 2 |
. 3 |
… |
. N |
|
|
|
|
|
|
|
x |
x |
x |
… |
x |
|
|||||
|
|
|
|
|
|
- |
|
|
|
|
|
|
x |
x |
x |
… |
x |
|
|||||
|
|
|
|
|
|
. 3.16. # # / / ( 8 + + )
+ ! !& ! % ! ! #! &$ U -#! N ( ! ! #! -M · N). & & #! # & - ! (. . ! # & &$), ! #! % ! ' #! &$ ( U), < U &$ , & * $& & .
! # & ! % 4 ! * ! # &$ $ $ ($ # 4-). G , U ! ! ! %-! ! .
& & #& * &! ! ! . ? ! * ! - . ! * !& ( # ! &$ ! ), . . % &-
V V T E , |
} – ! , V – ! , - |
+ # , & |
* &! ! ! . ; |
|
77 |
&$ ! !, $$, % & , -* ! ! ! &$ &$.
A , $ &$ ! + * ! & & ( #! M ·N), & &$ !*+ : 1) & « *+ » ! S
( #! * U · U): S (A A |
)(A A |
)T , A – ! , + |
mean |
mean |
mean |
# #! % N # $ # % ! & &; |
||
# # ; 2) ! & S & * - |
& & # . % ! ! # * !- & S & . & &* &! ! ! !& . J ! ' # , ! ' -$&$ &$ *+ .
< & & &$ # - , ! ! U, #! M · N ! % & ! '#! M · , * & , *+ !- ! &! &! # ! ( * ! % -! – & 20 % &$ - ). & &$, *+ $ &$ & & , ! ! % ! A0 UT A .
, 4 : -: : * 4 D . A ! ! !! #! * nr · nc · nt,nr – ! !! , nc – , nt – - ! !!, # ! <= #! % $:
1.# ! &$ ! : U ! # ! nt -
!&$ ( #! $&$ &$ ! ! !! nt),% # &$ & nr · nc #! . -
$, # ! Amean , U !! 4 . < #! ! &U ! nt · nt. ! <= ! ' #- ! * ! &$, % *+$ # $&$! &$ ! % ! . U !$ ! P & 4 & #! ! !!! &!, . . nr·nc U ! , ! &$ ( ! !!
#! $&$ &$) ! % & ! ' .
2.# &$ ! : U ! # ! nr · nc ! &$, % # &$ nt #! . -- $, # ! U, U ! -
78
# % . < #! ! & U ! ne · ne, ne = nr · nc. - ! <= ! ' #! * ! &$, % *+$ # ! & . -$ ! H ! ' &$ 4-#! * ! #! & ! !!& ( nr nc),% * #! ! !!&, % % -&$ # . = & # #&$ ! % & !. A nr = 240, nc = 320, nt = 150, ! % $ # * -! !!, # ! &$ ! #! -! & 1508150 U ! , # -&$ ! #! ! & 76 800876 800 U ! , & ! % ! &$ &$ V !$ ! ! * .
"! #! % # U ! #$ 4 - ! . U ! & ! & & -, ! % !
&! ! 4 . |
U ! - |
& |
& ! &, & , |
! # !, # U ! # % |
|
4 . |
|
% m-4 P= MATLAB, - |
# & &$ ! $! ! &$ &$ & ! &.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [ newData, pcs, vals] = Pca( data, cType) %function [ newData, pcs, vals] = Pca( data, cType) %data – ! &$ #! * m * n %m – ! &$, n – *
%cType – ! & & &$ ! %newData – $ &$ &$ & ! & %(4 ! data)
%pcs – & ! & ( - &)
%vals – & |
# |
% & |
# , . . |
% ! & U !
%[nr, nc] = size( data);
%mn = mean( data, 2);
%data = data – repmat( mn, 1, nc);
% & ! * * ! if( cType == 'cov')
cMat = cov( data'); else
if( cType == 'cor')
79
cMat = corrcoef( data'); else
error('cType must be: cor or cov'); end
end;
[eVect, eVals] = eig( cMat); eVals = diag( eVals);
[temp, inds] = sort( -1*eVals); eVect = eVect( :, inds);
pcs = eVect; newData = pcs' * data; vals = eVals( inds);
! # $ $! -! !! ! * ! ! % #-m-4 «SeqForTimePca» «TimePcaForSeq»:
function Res = SeqForTimePca( Seq)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%function Res = SeqForTimePca( Seq)
% # 3-! [nr, nc, nt] % [nt, nr*nc]
[nr, nc, nt] = size( Seq);
Res = reshape( Seq, nr*nc, nt); Res = Res';
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function Res = TimePcaForSeq( M, nr, nc) %function Res = TimePcaForSeq( M, nr, nc)
% # 2-! [nt, nr*nc] % [nr, nc, nt]
nt = size( M, 1);
M = M';
Res = reshape( M, nr, nc, nt);
. 3.17 & # & ! ! # - &$ ! . ! % ,! # &$ ! ' - ' /' ! * $&! ! !! ! .
. 3.17. / 1, 2, 3 4-* / #" /
! ! &' &$ & # &- * !& &$, , ! ! , # , > , P *, < ! . ., %& [7, 19, 20]. "$ &$ # ! .
80