finding hidden components

Upload: lucas-gallindo

Post on 07-Jul-2018

217 views

Category:

Documents


0 download

TRANSCRIPT

  • 8/18/2019 Finding Hidden Components

    1/21

      

      

    N SF- CBM SL

     e c t ur e S eri  e s

    L  e c t  ur e1 : 

    Fi  n d i  n gHi   d  d  en C om p on en t  s 

    B y

    D a t  aAn al   y s i   s 

    D a vi   d D on oh  o

     S  t   a t  i   s  t  i   c s D e p ar  t  m en t  

     S  t   anf   or  d  Uni   v er  s i   t   y

      

      

    A g en

     d  a

    Hi   d  d  en C om p on en t   s 

    Pr i  n ci   p al   C om p on en t   s 

    B e y on d Pr i  n ci   p al   C om p on en t   s 

     S  t   a t  i   s  t  i   c al  An al   y s i   s  of  I  m a g e s 

    L  o ok Ah  e a d 

    N S F- CBM S T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

     3 

       

       

    I  .AHi   d 

     d  en C om p on en t  s i  nL i   gh  t 

      (   a  )  N e w t  on

      (   b   )  Pr i   s m

    N S F- CBM S 

    T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

       

       

    Fr omN e w t  on

    ’   s N o t  e b  o ok 

  • 8/18/2019 Finding Hidden Components

    2/21

      

      

    I  .BHi   d 

     d  en C om p on en t  s i  nH e a t 

    F i   g ur  e2 : F  o ur i   er 

      

      

    Di  ff  u s i   on

     of  H e a t 

    I  ni   t  i   al  T em p er  a t   ur  ePr  ofi l   e: 

      f    (   θ    , t  = 0   )    ,f   or  θ   ∈  [    0   ,2  π  )  

    Pr  o

    fi l   e a t  l   a t   er  t  i  m e: 

      f    (   θ    , t    )    ,f   or  θ  

     ∈  [    0   ,2  π  )  

     S h  o

     w s 

    l   e s  s  v ar i   a t  i   on

    m or  e s m o o t  h n e s  s 

     e v en t   u al   uni  f   or mi   t   y

    N S F- CBM S T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

     7 

       

       

     0 

     3 

     5 

     6 

    - 0 . 5    0 

     0 . 5 1 

    1 . 5 

    f(θ,t)

     θ  

     S  t   a g e s  of  D i  f  f   u s i   on

    N S F- CBM S 

    T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

     8 

       

       

    F o uri   er’   s 

    I  n s i   gh  t 

    Ex p an d i  n S i  n u s  oi   d  al   C om p on en t   s 

      f    (   θ    , t  = 0   )  =      

     n

     c  n ex p  {   i   n θ    }  

    Pr  o

    fi l   em a d  e of  ‘   w a v e s ’  .

    Pr  o

    fi l   e a t  l   a t   er  t  i  m e: 

      f    (   θ    , t    )  =      

     n

     λ  n  (   t  

      )  ·  c  n ex p  {   i   n θ    }  

    D am pi  n gf   a c t   or : 

     λ  n  (   t    )  = ex p  {  − n2  t    }  

    D e cr  e a s  e s r  a pi   d l   y wi   t  h  w a v el   en g t  h  ,

     t  i  m e

  • 8/18/2019 Finding Hidden Components

    3/21

      

      

     0 

     5 

    - 0 . 5    0 

     0 . 5 1 

    1 . 5 

    f(θ,0)

     θ  

    I  ni   t  i   al  D  en s i   t   y  t  = 0 

     0 

     5 

     0 1 2 3 4 5 6 7 8 9 1  0 

    Component

     θ  

     C  om p on en t   s  t  = 0 

     0 

     5 

     0 1 2 3 4 5 6 7 8 9 1  0 

    Component

     θ  

     C  om p on en t   s  t  =. 0 4  6  9 

     0 

     5 

    - 0 . 5    0 

     0 . 5 1 

    1 . 5 

    f(θ,t=.)

     θ  

    D  en s i   t   y 

     a t   t  =. 0 4  6  9 

      

      

    I  n t  er pr e

     t  a t i   on

    F  o ur i   er An al   y s i   s i   s  ak i  n d  of  Pr i   s m

    Hi   d  d  en C om p on en t   s :  S i  n u s  oi   d  a

    l   C om p on en t   s 

    P o s  t  F  a c t   o:  S i  n u s  oi   d  s F  un d  am e

    n t   al  

    D omi  n an t   t  r  en d  ,1  9  t  h  c en t   ur  yM a t  h  em a t  i   c al  Ph  y s i   c s 

    N S F- CBM S T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

    1 1 

       

       

    M o d  ern Vi   e w p oi  n t –I  

      f    (  ·   , 0   )  

    F o u ri  e r

    −→

      (   c  n  (   0   )    )  

    Di  ff  u s i   on

      ↓  

      ↓  

    × λ  n  (   t    )  

      f    (  ·   , t    )  

    F o u ri  e r−1

    ←−

      (   c  n  (   t    )    )  

    N S F- CBM S 

    T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

    1 2 

       

       

    M o d  ern Vi   e w p oi  n t –I  I  

    L i  n

     e ar Tr  an s f   or m a t  i   on

      (  T  c   )   n=       n

     m

    T  n  , m c  m

    Di   a

     g on al  Tr  an s f   or m a t  i   on

      (  T  c   )   n=T  n  , n c  n

    Di   a

     g on al  –n oi  n t   er  a c t  i   on b  e t   w e en c om p on en t   s .

  • 8/18/2019 Finding Hidden Components

    4/21

      

      

    M o

     d  ern Vi   e w p oi  n t –I  I  I  

      f    (  ·   ,

     0   )  F o u ri  e r

    −→

      (   c  n  (   0   )    )  

    Di  ff  u s i   on

      ↓  

      ↓  

    Di    a

      g on al   

      f    (  ·   ,

     t    )  F o u ri  e r−1

    ←−

      (   c  n  (   t    )    )  

      

      

    I  I  .Pri  n ci   p al   C om p on en t  s 

      (   a  )  K ar h  un en

      (   b   )  L  ò  e v e

    N S F- CBM S T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

    1  5 

       

       

      (   a  )  

    P e ar  s  on

      (   b   )  H o t  el  l  i  n g

    N S F- CBM S 

    T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

    1  6 

       

       

    K arh  un en-L  ò  e v eTr an s f   orm

    X  (   t    )   , t   ∈T  ar  an d  om pr  o c e s  s i  nL 2   (  T   )  

    Γ   (   t    , s  )  i   t   s  c o v ar i   an c e: 

     C o v  (  X  (   t    )    ,X  (   s  )    )  =E  (  X  (   t  

      )  −EX  (   t    )    )    (  X  (   s  )  −EX  (   s  )    )  

     Un d  er  g en er  al   c on d i   t  i   on s  ,∃  or t  h  o b  a si    s  φ 1  ,  φ 2  ,... ,

          φ  n  ,  φ  m

        = δ   n  , m

    Ei   g en v e c t   or  s 

             Γ   (   t    , s  )    φ  n  (   s  )   d  s= λ  n×  φ  n  (   t    )  

    Ei   g en v al   u e s 

     λ 1  ≥ λ 2 

     ≥ λ  3 ...

  • 8/18/2019 Finding Hidden Components

    5/21

      

      

    M o d  ern Vi   e w p oi  n t 

      f  

    KLT

    −→

      (   c  n  )  

     C o v ar i   an

     c e  ↓  

      ↓  

    × λ  n

    Γ   f  

    KLT−1

    ←−

      (   d  n  )  

      

      

    KLT R e pr e

     s  en t  a t i   on

    X  (   t    )   , t   ∈T  ar  an d  om pr  o c e s  s i  nL 2   (  T   )  

    X n=    X  ,  φ  n     ,i   t   s  n- t  h  c o or  d i  n a t   ei  nKL  b  a s i   s 

     V ar i   an c e s 

     V  a rX n

    = λ  n

     C o v ar i   an c e s 

     C o v  (  X n  ,X m

      )  = λ  n× δ   n  , m

     R e pr  e s  en t   a t  i   on

    X  (   t    )  =      

     n

    X n·   φ  n  (   t    )  

    h  er  e t  h  eX n

     ar  e or  t  h  o g on al  .

    N S F- CBM S T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

    1  9 

       

       

     S  om eK e yA p pl  i   c a t i   on s 

    •N oi   s  e R em o v al  – Wi   en er 

    •D a t   a C om pr  e s  s i   on

    – S h  ann on

    N S F- CBM S 

    T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

    2  0 

       

       

    N oi   s  e R em o v al  

    Y   (   t    )   o b  s  er  v e d  d  a t   a

    X  (   t    )   d  e s i  r  e d  s i   gn al   t   or  e c o v er 

    Y   (   t    )  =X  (   t  

      )   +N  (   t    )  

    N  (   t    )   G a u s  s i   ann oi   s  e pr  o c e s  s  , wi   t  h 

     S  am eKL  b  a s i   s  a s X:   (    φ  n  )  

    Ei   g en v al   u e s   (    µ n  )  E

      {  X  (   t    )    |   Y   }  =      

     n

     w nY  n  φ  n  (   t    )  

     wh  er  e

     w n=

     λ  n

     λ  n

     +  µ n  .

  • 8/18/2019 Finding Hidden Components

    6/21

      

      

    M o d  ern Vi   e w p oi  n t 

    KLT

    −→

      (  Y  n  )  

    Ex p e c t   a t  i   on

      ↓  

      ↓  

    × w n

    E  {  X  (   t    )    |   Y   }  

    KLT−1

    ←−

      (  ˆ X

     n  )  

      

      

    D a t  a C om

     pr e s  s i   on

    N S F- CBM S T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

    2  3 

       

       

    K e yEx am pl   e

    B  (   t    )  Br  o wni   anBr i   d  g e

    .

    • G a u s  s i   an

    •B  (   0   )  =B  (  1   )  = 0 

    •Γ   (   t    , s  )  = m i   n  (   t    , s  )  

    − t   s.

    • Cl   o s  el   yr  el   a t   e d  t   o

    Br  o wni   anM o t  i   on

    N S F- CBM S 

    T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

    2 4 

       

       

     0 

     0 .1 

     0 .2 

     0 . 3 

     0 .4 

     0 . 5 

     0 . 6 

     0 .7 

     0 . 8 

     0 . 9 

    -2  0 

    -1  5 

    -1  0  - 5    0 5 

    1  0 

    1  5 

    2  0 

    2  5 

     3  0 

     t  

    B(t)

    B r  owni   anB r i   d  g e

  • 8/18/2019 Finding Hidden Components

    7/21

      

      

    KL T of  Br o wni   an

    •  φ  n  (   t    )  =      √2  s i  n  (   π n t    )  

    •  φ  n  (   0   )  =  φ  n  (  1   )  = 0 

    • λ  n= C o n s t    /   n2 

    • R e pr  e s  en t   a t  i   on: 

    B  (   t    )  =      

     n

    Z  n  φ  n  (   t    )  

    •Pr i  n ci   p al   S  c or  e s Z  n

    i  n d  e pN  (   0   , λ  n  )  .

    “  S i  n u s  oi   d  s  ar  e t  h  ePr i  n ci   p al   C om p on en t   s  of  Br  o wni   an

    Br i   d  g e.” 

      

      

    I  I  I  .B e y on d Pri  n ci   p al   C om p on en t  s 

    I  n d  e p en d  en t   C om p on en t   s 

     S  p ar  s  e C om p on en t   s 

    N S F- CBM S T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

    2  7 

       

       

    I  I  I  .AI  n d  e p en d  en t  C om p on en t 

     s 

    X  (   t    )  =      

     n

      ξ   n  ψ  n  (   t    )  

    •  ξ   n

    I  n d  e p en d  en t  

    •  ψ  n

    B a s i   s F  un c t  i   on

     s 

     R em ar k  s 

    •N e e d n o t   exi   s  t  

    •N e e d n o t   b  e or  t  h  o

     g on al  

    •N o pr  o c e d  ur  ei  n g

     en er  al  

    N S F- CBM S 

    T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

    2  8 

       

       

    I  n d  e p en d  en t  C om p on en t  s  of   G a u s  s i   anPr o c e s  s 

    F  or 

     G a u s  s i   anPr  o c e s  s  ,I   C=P C.

    Pr  o

     of  : 

    I  f    (  X n  )   un c or r  el   a t   e d  & j   oi  n t  l   y G a u s  s i   an , t  h  e  (  X n  )   ar  e

    i  n d  e p en d  en t  .

    F  or  a G a u s  s i   anPr  o c e s  s  , t  h  e pr i  n ci   p al   s  c or  e s  ar  e un c or r  el   a t   e d 

     an d  j   oi  n t  l   y G a u s  s i   an.

  • 8/18/2019 Finding Hidden Components

    8/21

      

      

     R am p: 

    AN on G a u s  s i   anPr o c e s 

     s 

     S  c al   ar  R V ω∼ U   (   0   ,1   )  

    .

     S  t   o ch  a s  t  i   cPr  o c e s  s   (  Y.M e y er   )  

    X

      (   t    , ω  )  =                t  

     t   < ω

     t  −1 

     t   ≥ ω

    Pr  o p er  t  i   e s : 

    •EX  (   t    )  = 0 

    •X  (  1   )  =X  (   0   )  = 0 

      

      

     0 

     0 .1 

     0 .2 

     0 . 3 

     0 .4 

     0 . 5 

     0 . 6 

     0 .7 

     0 . 8 

     0 . 9 

    - 0 .4 

    - 0 .2    0 

     0 .2 

     0 .4 

     0 . 6 

     0 . 8 1 

     t  

    Ramp(t)

    R  am

     p

    N S F- CBM S T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

     3 1 

       

       

     C

     o v ari   an c e of   R am p

    •Γ   (   t    , s  )  = C o v  (  X  (   t  

      )    ,X  (   s  )    )  = m i   n  (   t    , s  )  − t   s

    • S  e c on d - Or  d  er E q u

    i   v al   en t   t   oBr  o wni   anBr i   d  g e!   

    N S F- CBM S 

    T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

     3 2 

       

       

     0 

     0 .1 

     0 .2 

     0 . 3 

     0 .4 

     0 . 5 

     0 . 6 

     0 .7 

     0 . 8 

     0 . 9 

    -2  0 

    -1  0    0 

    1  0 

    2  0 

     3  0 

    4  0 

     5  0 

     t  

    B(t)

    B r  owni   an

    B r i   d  g e

     0 

     0 .1 

     0 .2 

     0 . 3 

     0 .4 

     0 . 5 

     0 . 6 

     0 .7 

     0 . 8 

     0 . 9 

    - 0 .4 

    - 0 .2    0 

     0 .2 

     0 .4 

     0 . 6 

     0 . 8 1 

     t  

    Ramp(t)

    R  am

     p

  • 8/18/2019 Finding Hidden Components

    9/21

      

      

     C on t r a s  t 

    •Pr i  n ci   p al   C om p on

     en t   s  ar  e s  e c on d - or  d  er i  n d  e p en d  en t   , b  u t  

     p er h  a p s n o t  i  n d  e p

     en d  en t  .

    •I  n d  e p en d  en t   C om p on en t   s  ar  ei  n d  e p en d  en t   , b  u t   p er h  a p s  c ann o t  

     b  ef   o un d  or  d  on o

     t   exi   s  t  .

      

      

    I  I  I  .BA p pr ox

    i  m a t  eI   CA

    A t   t  

     em p t   t   o a p pr  oxi  m a t   eHi   gh  er - Or 

     d  er I  n d  e p en d  en c e

    F  or 

    m C u m u l    a n t  T e n s  o r  of   or  d  er  >

    A t   t  

     em p t   t   oDi   a g on al  i  z eT en s  or 

    N S F- CBM S T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

     3  5 

       

       

    F o ur t h 

    - or d  er J  oi  n t  C um ul   an t  s 

    K ur  t   o s i   s : K4   (  X  )  =E  (  X−  µ  )  4 − 3  V  a r  (  X  )  2 

     C u m4   (  X1   ,...  ,X4   )  

    =

    E  (  X1 X2 X 3 X4   )  

    E  (  X1   )  · E  (  X2 X 3 X4   )  −E  (  X2 

      )  · E  (  X1 X 3 X4   )  

    E  (  X 3   )  · E  (  X1 X2 X4   )  −E  (  X4 

      )  · E  (  X1 X2 X 3   )  

     +

    E  (  X1 X2   )  · E  (  X 3 X4   )  

     +

    E  (  X1 X 3   )  · E  (  X2 X4   )   +E  (  X1 X4   )  · E  (  X2 X 3   )  

     +

    2   (  EX1   )    (  EX2   )    (  EX 3   )    (  EX4   )  

    N S F- CBM S 

    T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

     3  6 

       

       

    I  n d  e p en d  en c e an d  C um ul   an t T en s  or s 

    L  e t  

    X n=    X  ,  φ  n     c o effi ci   en t  i  n b  a s i   s   (    φ  n  )  .

     C um ul   an t  T en s  or K4   ,  φ 

    4 - w a y ar r  a y

    K4   ,  φ   (   i  1   , i  2   , i   3   , i  4   )  = C u m

    4   (  X i  1  ,X i  2  ,X i   3  ,X i  4  )  

    I    f     (  X n  )   a r  e i    n d   e  p e n d   e n t   r  a n d   o m

     v a r  i    a b   l    e s   ,K4   ,  φ 

     i    s  h    y  p e r  d   i    a  g o n a l   : 

    K4   ,  φ   (   i  1   , i  2   , i   3   , i  4   )  =                0 

    n o t   al  l  i  n d i   c e s  s  am e

     i  1 = i  2 = i   3 = i  4 

  • 8/18/2019 Finding Hidden Components

    10/21

      

      

    N u

    m eri   c al  Ex p eri  m en t 

    F i  ni   t   e R am pPr  o c e s  s : 

    •T =  {  1   /   3 2   ,...  , 3 2   /   3 2   }  

    •Ω=  {  1   /   3 2   ,...  , 3 2   /   3 2   }  

    •D a t   a b  a s  eX  (   i    ,  j    )  =X  (   ω i    , t    j   )  .

    –N

    = 3 2 r  e al  i  z a t  i   on s X  (   i    ,·   )  

    –N

    = 3 2  v ar i   a b l   e s X  (  ·   , t    )  

    • U s  e J  e an-F r  an c oi   s 

     C ar  d  o s  o“  J ADE” Al   g or i   t  h m

    •A t   t   em p t  Di   a g on al  i  z a t  i   on of  Em pi  r i   c al   C um ul   an t  T en s  or 

    •M o v e s  9  0  % of  K ur  t  i   cEn er  g y onDi   a g on al  .

      

      

     R e s  ul   t  s f  r om J ADE

     an al   y s i   s  of   R am p

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

    N S F- CBM S T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

     3  9 

       

       

    D a u b  e ch i   e s N e arl   y S  ymm e t ri   c W a v el   e t  s 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

     0 

    2  0 

    -1    0 1 

    N S F- CBM S 

    T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

    4  0 

       

       

     C om p ari   s  on s 

     0 

    2  0 

    -1 - 0 . 5    0 

     0 . 5 1 

     J  a d  el   e t  1  6 

     0 

    2  0 

    -1 - 0 . 5    0 

     0 . 5 1 

    W av  el   e t  1 2 

     0 

    2  0 

    -1 - 0 . 5    0 

     0 . 5 1 

     J  a d  el   e t   5 

     0 

    2  0 

    -1 - 0 . 5    0 

     0 . 5 1 

    W av  el   e t  2 1 

     0 

    2  0 

    -1 

    - 0 . 5    0 

     0 . 5 1 

     J  a d  el   e t   9 

     0 

    2  0 

    -1 

    - 0 . 5    0 

     0 . 5 1 

    W av  el   e t  2  3 

     0 

    2  0 

    -1 - 0 . 5    0 

     0 . 5 1 

     J  a d  el   e t  1  3 

     0 

    2  0 

    -1 - 0 . 5    0 

     0 . 5 1 

    W av  el   e t  2  5 

  • 8/18/2019 Finding Hidden Components

    11/21

      

      

    I  n t  er pr e t  a t i   on s 

    Hi   d  d  en C om p on en t   s  d 

    i   s  c o v er  e d  b  y J ADE

    •D y a d i   c S  c al   e s 

    •N e ar -Tr  an s l   a t   e s 

    •A p pr  ox.D a u b  e ch i   e s N e ar l   y S  ymm e t  r i   c

    “  Q u a s i  - W a v el   e t   s  ar  e t  h  e al  m o s  t  -i  n d  e p en d  en t   c om p on e

    n t   s  of   R am p” 

      

      

    I  I  I  .D S  p ar s  e C om p on en t  s An al   y s i   s 

    X  (   t    )  =      

     n

    X n  φ  n  (   t    )  

      (  X n  )   a s  p ar  s  e v e c t   or : f   e wl   ar  g e

     en t  r i   e s .

      φ  n a b  a s i   s  or  d i   c t  i   on ar  y– t   o b  ef   o un d 

     G o al  – ch  o o s  e  (    φ  n  )   t   o o p t  i  mi  z e

     s  p ar  s i   t   y

     R em ar k  s : 

    N e e d n o t   exi   s  t  

    N o g en er  al   d  efi ni   t  i   on

    N o g en er  al   al   g or i   t  h m

    N S F- CBM S T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

    4  3 

       

       

     C om p a

    ri   s  on of   S  CA wi   t h P CA

    P CA

    •A

      (   m  )    (  X  ,  (    φ  n  )    )  fi r s

     t  m- t   er m a p pr  oxi  m a t  i   on ,

    A  (   m  )  X

    =

     m       n

    =1   c  n  φ  n

    •F i  n d B a s i   s   (    φ  n  )   o p t  i  mi  zi  n g

     ar  gmi  n

      (    φ  n  )  E    X−A

      (   m  )    (  X  ,  (    φ  n  )    )      2 2 

    N S F- CBM S 

    T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

    4 4 

       

       

     C om p ari   s  on of   S  CA

     wi   t h P CA  (   c on t   )  

     S  CA  (   On e a p pr  o a ch   )  

    B  (   m  )    (  X  ,  (    φ  n  )    )   b  e s t  m- t   er m a p pr  oxi  m a t  i   on ,

        X−B

      (   m  )    (  X  ,  (    φ  n  )    )      2 = a

    r  gmi  n

      (   n  j   )      X−

     m        j 

    =1   c  n  j   φ  n  j     2 

    F i  n d B a s i   s  o p t  i  mi  zi  n g

     ar  gmi  n

      (    φ  n  )  E    X−

    B  (   m  )    (  X  ,  (    φ  n  )    )      2 2 

  • 8/18/2019 Finding Hidden Components

    12/21

      

      

    Ex am pl   e of  A p pr oxi  m a t  e S  CA:  R am p

     W a v el   e t   R e pr  e s  en t   a t  i   o

    nX  (   t    )  =      

     n

     α  j   , k   ψ   j   , k   (   t    )  

    •  (    ψ   j   , k   )   Or  t  h  on or m al   W a v el   e t   s  b  a s i   s  , 3  v ani   s h i  n gM

     om en t   s .

    •  (   α  j   , k   )   s  p ar  s  e v e c t   or : f   e wl   ar  g e en t  r i   e s .

    •A t  l   e v el    j   a t  m o s  t  

     Cn onz er  o s  , ∀  k 

    • Wi   t  h  5 4  c o effi ci   en

     t   s  c anr  e pr  e s  en t   t   o a c c ur  a c y1  0 − 6 

      

      

     S  p ar s i   t  y C o

    m p ari   s  on

     W a

     v el   e t  B a s i   s   (   q u a s i  - S  CA  )  

     Wi   t  h  5 4  c o effi ci   en t   s  c anr  e pr  e s  en t   t   o a c c ur  a c y1  0 − 6 

        X−B

      (   m  )    (  X  , wa velet s   )      2  ≤ C1  ex p  (  − C2  m  )  

    F  o u

    r i   er B a s i   s   (  P CA  )  

    N e e d 1  9  0 2  c o effi ci   en t   s f   or  a c c ur  a c y1  0 − 6 .

        X−A

      (   m  )    (  X  , s i  n u s  oi  d s   )      2  ≤ C  /        √ m

        X−B

      (   m  )    (  X  , s i  n u s  oi  d s   )      2  ≤ C  /        √ m

    N S F- CBM S T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

    4  7 

       

       

     0 

    2  0 

    4  0 

     6  0 

     8  0 

    1  0  0 

    1 2  0 

    1 4  0 

    1  6  0 

    1  8 

     0 

    2  0  0 

    - 8  -7  - 6  - 5  -4  - 3  -2  -1    0 

    n

    Log(cn)

     C  om pr  e s  s i   onN  um b  er  s f   or D  S T  an d W av  el   e t  A n al   y  s i   s  of  R  am p

     S i  n eT r  an s f   or m

    W av  el   e t  T r  an s f   or m

    N S F- CBM S 

    T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

    4  8 

       

       

     0 

     0 .2 

     0 .4 

     0 . 6 

     0 . 8 

    - 0 .4 

    - 0 .2    0 

     0 .2 

     0 .4 

     0 . 6 

     0 . 8 1 

    R  am p

     t  

    Ramp(t)

     0 

     0 .2 

     0 .4 

     0 . 6 

     0 . 8 

    -1  0  - 8  - 6  -4  -2 

    W av  el   e t  T r  an s f   or m of  R  am p

     0 

    1  0  0 

    2  0  0 

     3  0  0 

    - 0 . 5    0   0 . 5 

     S i  n eT r  an s f   or m of  R  am p

     0 

    1  0  0 

    2  0  0 

     3  0  0 

    1  0 - 6  1  0 -4  1  0 -2  1  0  0 

     Or  d  er  e d T r  an s f   or m C  o ef  f  i   c i   en t   s 

    k - t  h l   ar  g e s  t   c  o ef  f  i   c i   en t  

    Coefficient Amplitude

     S i  n eT r  an s f   or m

    W av  el   e t  T r  an s f   or m

  • 8/18/2019 Finding Hidden Components

    13/21

      

      

     C on t r a s  t  of  I   CA  /  S  CA  / P CA

    I   CA

     S  CA

    N e w C on c e p

     t  Hi   gh - Or  d  er 

    B e s  t  -mT er 

    m

    I  n d  e p en d  en c e

    A p pr  oxi  m a t  i   on

    P CA C on c e p t   S  e c on d - Or  d  er 

    F i  r  s  t  -mT er m

    I  n d  e p en d  en c e

    A p pr  oxi  m a t  i   on

      

      

     C om p ari   s  on o

    f  I   CA  /  S  CA

    A s  p e c t  

    I   CA

     S  CA

    P o p ul   ar i   t   y

    I  mm en s  e

     R ar i  fi  e d 

    M o d  e

     S  t   a t  i   s  t  i   c al  

    A p pr  oxi  m a t  i   onTh  e or  y

    H ar m oni   cAn al   y s i   s 

    Pl   a u s i   b i  l  i   t   y

    N on e

     S  uffi ci   en t  

    Al   g or i   t  h m s 

    M an yH e ur i   s  t  i   c

     s F  e wH e ur i   s  t  i   c s 

    N S F- CBM S T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

     5 1 

       

       

     S i  mi  l   ar O u t  c om ef  r omI   CA & S  CA: 

     R am p

    • W a v el   e t   s  gi   v e v er  yh i   gh r  a t   e of   m- t   er mA p pr  oxi  m

     a t  i   on.

    • W a v el   e t   s  al  m o s  t   d 

    i   a g on al  i  z e t  h  e c um ul   an t   t   en s  or .

     C anfi n d  q u al  i   t   a t  i   v el   y

     s  am e o u t   c om e b  y ei   t  h  er  pr i  n ci   p

    l   e.

    N S F- CBM S 

    T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

     5 2 

       

       

    Al  m o s  t -Di   a g on al  i   t  y o

    f   C um ul   an t T en s  or s 

        K4   ,  φ     AD  (    p  )  = s  u p

     i  1

      (    i  1  , i  

    2  , i   3

      |   K4   ,  φ   (   i  1   , i  2   , i   3   , i  4   )    |     p  )  1   /   p

    I  n t   er  e s  t  i  n g c a s  e ;   p≈ 0  ;  s  a y s  e a ch  s l  i   c e s  p ar s e.

    Th  e or em.F  o r  c i    r  c u l    a r  R a m  p  p r  o c e s  s  o n  [    0   ,2  π  )    ,

    D a u b   e c h   i    e s  n e a r  l     y

     s   y m

     m e t   r  i    c w a v e l    e t   s  w i    t   h   ≥ 3  v a n i    s  h  

     i    n  g m o m e n t   s 

        K4   ,  ψ    AD  (    p  )   < ∞

    Al  m

     o s  t  - d i   a g on al  i  n s  en s  e of   g o o d  d  e c a y a w a yf  r  omh  y p er - d i   a g on al  .

  • 8/18/2019 Finding Hidden Components

    14/21

      

      

    I  n t  er pr e t  a t i   on

    • S  t   u d  y1 − d  o b  j   e c t   s  wi   t  h  s i  n g ul   ar i   t  i   e s 

    •Th  e y ar  en on g a u s  s i   an

    •P CA w o ul   d  gi   v e s 

    i  n s  u oi   d  al   c om p on en t   s 

    •I   CA & S  CA gi   v e w a v el   e t   c om p on en t   s 

     Wh  a t  i  f   w e s  t   u d  y2 − d  o b  j   e c t   s –I  m a  g e s?  

      

      

    I   V. Vi   s  u al  N e ur o s  ci   en c e

    I  nf   or m a t  i   onF l   o wi  n t  h  e Vi   s  u al   S  y s 

     t   em

    I  m a g e

     R e t  i  n a

    ...

     V1 →

    ...

     V1 

     C om p on en t   of   Vi   s  u al   C or  t   ex: 

    L i  n e ar  R e c e p t  i   v eF i   el   d  s –H u b  el   an d  Wi   e s  el  .

    Ne u ral O utp ut n=    I  ma ge  ,  φ  n    

    Pr  o p er  t  i   e s  of    (    φ  n  )  

    –L  o c al  i  z e d 

    – Or i   en t   e d 

    – W a v el   e t    /   G a b  or M o d  el   s 

    N S F- CBM S T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

     5  5 

       

       

     C ari   c a t  ur e of   V1 

    N S F- CBM S 

    T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

     5  6 

       

       

    I   V.AT w o

     S l   o g an s 

      (   a  )  B ar l   o w

      (   b   )  Fi   el   d 

  • 8/18/2019 Finding Hidden Components

    15/21

      

      

    B arl   o w

    • Vi   s  u al   S  y s  t   em d  o e s  ar  em ar k  a b l   e j   o b  of  i  m a g e c om

     pr  e s  s i   on.

    •1  0  7  b i   t   s   /   s  e c a t  i  n p u t  

    •2  0 −4  0  b i   t   s   /   s  e c a t   en d  of   pr  o c e s  s i  n g

     Vi   s  u al   s  y s  t   em t   ak  e s h i   gh l   yr  e d  un d  an t   d  a t   a ,r  em o v e s  t  h  e

    r  e d  un d  an c y , pr  o d  u c e s 

    n e ar l   yi  n d  e p en d  en t   d  a t   a.

      

      

    Fi   el   d 

     S  t  r  u c t   ur  e of  n a t   ur  al  i  m a g e s  d r i   v e s  s  t  r  u c t   ur  e of   vi   s  u al   s  y s  t   em

     R e d  un d  an c y of  n a t   ur  al  i  m a g e s 

     c o ul   d  t   e a ch  vi   s  u al   s  y s  t   emh  o w

     t   o b  e s  p ar  s  e

     S  t   a

     t  i   s  t  i   c al   an al   y s i   s  of  i  m a g e s  an d n e ur  o d  e v el   o pm en t  ?    (   C a t   s   )  

    N S F- CBM S T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

     5  9 

       

       

     S  p ar s  e C om p on en t  s  of  N a t  ur al  I  m

     a g e s 

      (   a  )   Ol   s h  a u s  en

      (   b   )  Fi   el   d 

    N S F- CBM S 

    T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

     6  0 

       

       

     Ol   s h  a u s  en-Fi   el   d Ex p eri  m en t 

    A p p e ar  e d i  nN a t   u r  e ,1  9  9  6 

     C ol  l   e c t  i   on of  n a t   ur  al  i  m a g e s 

    Ex t  r  a c t  1  6 ×1  6  p a t   ch  e s .

    D a t   a b  a s  e of   p a t   ch  e s 

     S  p ar  s  e C om p on en t   s An al   y s i   s 

    D e pi   c t  r  e s  ul   t   s 

  • 8/18/2019 Finding Hidden Components

    16/21

      

      

    D a t  a

     b  a s  e of  I  m a g eP a t  ch  e s 

    • C ol  l   e c t  i   on of  n a t   u

    r  al  i  m a g e s 

    •Ex t  r  a c t  1  6 ×1  6  p a t   ch  e s .

    •P a t   ch P   (    p  )  =  (  P   (    p  )  

      (   s  , t    )    ,1  ≤ s  , t   ≤1  6 

    •X  (    p  ,·   )  = en t  r i   e s  o

    f  P   (    p  )  i  nr  o w-m a j   or  or  d  er .

    • R e s  ul   t  : X

     anN

     b 

     y2  5  6 m a t  r i  x

      

      

    Em pi  ri   c al   S  p ar s  e C o

    m p on en t  s An al   y s i   s 

     W an t  X  (    p  ,  j    )  ≈     θ    j   (    p  )    φ   j 

    F i  n d  b  o t  h  c o effi ci   en t   s  θ    j   (    p  )   an

     d  b  a s i   s   φ   j .

     W an t   θ    j   (    p  )   s  p ar  s  e.

    F  or m al   O b  j   e c t  i   v e

    mi  n

      (    φ   j   )  

    mi  n

      (   θ   j   (    p  )    )    

      p

        X  (    p  ,·   )  −

          X  j   (    p  )    φ   j     2  + λ  S   (   θ    (    p  )    )  

     S  p ar  s i   t   ym e a s  ur  e

     S   (   θ    )  =      

      j 

    l   o g  (  1  + θ  2   j   )  

    N S F- CBM S T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

     6  3 

       

       N S F- CBM S 

    T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

     6 4 

       

       

    I  n t  er pr e

     t  a t i   on

    M ul   t  i  - Or i   en t   e d 

    M ul   t  i   s  c al   e

    B an d  p a s  s 

    I  ni   t  i   al  l   y , w a v el   e t   s .L  a t   er n o t  i  n t  i   t  l   e.

     Cl   e

     ar l   y  (   t   om e  )  n o t   w a v el   e t   s .

  • 8/18/2019 Finding Hidden Components

    17/21

      

      

     Wh  a t Ar eTh  e s  eTh i  n g s ? 

    •N o t   W a v el   e t   s – t   o

     om an y d i  r  e c t  i   on s 

    •N o t   G a b  or – t   o om an y s  c al   e s 

    • S i  mi  l   ar i   t  i   e s  t   o R e c en t   CHAi   d  e a s 

    – Ri   d  g el   e t   s 

    – C ur  v el   e t   s 

    –B e aml   e t   s 

      

      

    N S F- CBM S T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

     6  7 

       

       

    I  n d  e p en d  en t  C om p on en t  s  of  N a t  ur al  

    I  m a g e s 

      (   a  )  B el  l  

      (   b   )   S  e j  n o wk  s i  

    N S F- CBM S 

    T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

     6  8 

       

       

    I  n d  e p en d  en t  C om p on en t  s An al   y s i   s 

    I   CA-2  0  0  0  C onf   er  en c e

    NI  P S  C onf   er  en c e s 

    H un d r  e d  s  of   p a p er  s 

    L  e wi   ck i  - S  e j  n o wk  s i  

    H y v ar i  n en

     v an d  er  S  ch  a af   an d  v anH a t   er  en

  • 8/18/2019 Finding Hidden Components

    18/21

      

      

    I  n d  e p en d  en t  C om p on en t  s  of  N a t  ur al  

    M o vi   e s 

    F i   g ur  e 8 : H an s  v anH a t   er  en

      

      

    A R ol   ef   orM a t h ? 

    P CAD a t   aAn al   y s i   s fi n d  s  S i  n u s 

     oi   d  s 

    I   CA  /   S  CAD a t   aAn al   y s i   s  of  1 - d 

    Di   s  c on t  i  n ui   t  i   e s fi n d  s  W a v el   e t   s 

    I   CA  /   S  CAD a t   aAn al   y s i   s  of  2 - d 

    Di   s  c on t  i  n ui   t  i   e s fi n d  s ?  ?  ?  

    I   CA  /   S  CAD a t   aAn al   y s i   s  of   3 - d 

    Di   s  c on t  i  n ui   t  i   e s fi n d  s ?  ?  ?  

    L i   t   t  l   ePr  o s  p e c t   of   s  u b  s  t   an t  i   al  l   y

    l   ar  g er  d  a t   a an al   y s i   s 

    N S F- CBM S T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

     7 1 

       

       

    L  e c t  ur e s  t  o C om e

    •An al   o g y b  e t   w e en S  t   a t  i   s  t  i   c al  An al   y s i   s  of  I  m a g eD a

     t   a b  a s  e s  &

    H ar m oni   cAn al   y s i   s  of  I  m a g eM o d  el   s .

    • C omm onF  a c t   or :  S  p ar  s  e R e pr  e s  en t   a t  i   on

    •P oi  n t   S i  n g ul   ar i   t  i   e s  v s .L i  n e S i  n g ul   ar i   t  i   e s 

    •L i  n e S i  n g ul   ar i   t  i   e s 

     v s . C ur  v e S i  n g ul   ar i   t  i   e s 

    •H ar m oni   cAn al   y s i   s  of  L i  n e & C ur  v e S i  n g ul   ar i   t  i   e s 

    • Wh  a t  L i   e s B e y on d 

    N S F- CBM S 

    T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

     7 2 

       

       

    L  e c t  ur e

    Pl   an

    M

     on.

    T u e s .

     W e d .

    Th  ur  s .

    F r i  .

    D a t   aAn al  .

     Ri   d  g el   e t   s 

     Ri   d  g el   e t   s  C ur  v el   e t   s 

    B e aml   e t   s 

    M

     a t  h An al  .

     Ri   d  g el   e t   s 

     C ur  v el   e

     t   s  C ur  v el   e t   s 

     C om b i  n e d 

  • 8/18/2019 Finding Hidden Components

    19/21

      

      

     Wh 

     a t ’   s Th  e C onn e c t i   on? 

      (   a  )   Ri   d  g el   e t  s 

      (   b   )  I  m a g el   e t  s 

      

      

     Wh  a t ’   s  a

     t  S  t  ak  e

    A p pl  i   c a t  i   on s 

    –D a t   a C om pr  e s  s i   on

    –N oi   s  e R em o v al  

    –D e t   e c t  i   on , Cl   a s  s i  fi  c a t  i   on

    N e wA p pr  oxi  m a t  i   onTh  e or  y

    N e w C onn e c t  i   on s  wi   t  h M a t  h  em

     a t  i   c s 

    N e w R e s  e ar  ch i  nH ar m oni   cAn al   y s i   s 

    N S F- CBM S T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

     7  5 

       

       

     C ol  l   a b  or a t  or s   (  1   )  

    F i   g ur  e 9 : I   ai  n J  oh n s  t   on e

    N S F- CBM S 

    T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

     7  6 

       

       

     C ol  l   a b  or a

     t  or s   (  2   )  

      (   a  )  E. C an d `  e s 

      (   b   )  X.H u o

  • 8/18/2019 Finding Hidden Components

    20/21

      

      

     C

     ol  l   a b  or a t  or s   (   3  a b   )  

      (   a  )  A v er  b  u ch 

      (   b   )   C oi  f  m an

      

      

     C ol  l   a b  or a t 

     or s   (   3  c d   )  

      (   a  )  I   s r  a el  i  

      (   b   )  

     W al   d ´  en

    N S F- CBM S T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

     7  9 

       

       

     C ol  l   a b  or a t  or s   (  4   )  

      (   a  )   G.Fl   e s i   a

      (   b   )  E.Ar i   a s 

      (   c  )   O.L  e vi  

    N S F- CBM S 

    T al  k 1 : Fi  n d i  n  gHi   d  d  en C om  p on en t  s b   yD a t  aAn al    y si   s

     8  0 

       

       

     C ol  l   a b  or a

     t  or s   (   5   )  

      (   a  )  H.H el  - Or 

      (   b   )  Y.H el  - Or 

  • 8/18/2019 Finding Hidden Components

    21/21

      

      

    Th  ank  s ...

      (   a  )   S  t  o d  d  en

      (   b   )  D un c an