/ *   
 
   *   P R O J E C T :   N y A R T o o l k i t C S ( E x t e n s i o n ) 
 
   *   - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
 
   *   T h e   N y A R T o o l k i t C S   i s   J a v a   e d i t i o n   A R T o o l K i t   c l a s s   l i b r a r y . 
 
   *   C o p y r i g h t   ( C ) 2 0 0 8 - 2 0 0 9   R y o   I i z u k a 
 
   * 
 
   *   T h i s   p r o g r a m   i s   f r e e   s o f t w a r e :   y o u   c a n   r e d i s t r i b u t e   i t   a n d / o r   m o d i f y 
 
   *   i t   u n d e r   t h e   t e r m s   o f   t h e   G N U   G e n e r a l   P u b l i c   L i c e n s e   a s   p u b l i s h e d   b y 
 
   *   t h e   F r e e   S o f t w a r e   F o u n d a t i o n ,   e i t h e r   v e r s i o n   3   o f   t h e   L i c e n s e ,   o r 
 
   *   ( a t   y o u r   o p t i o n )   a n y   l a t e r   v e r s i o n . 
 
   *   
 
   *   T h i s   p r o g r a m   i s   d i s t r i b u t e d   i n   t h e   h o p e   t h a t   i t   w i l l   b e   u s e f u l , 
 
   *   b u t   W I T H O U T   A N Y   W A R R A N T Y ;   w i t h o u t   e v e n   t h e   i m p l i e d   w a r r a n t y   o f 
 
   *   M E R C H A N T A B I L I T Y   o r   F I T N E S S   F O R   A   P A R T I C U L A R   P U R P O S E .     S e e   t h e 
 
   *   G N U   G e n e r a l   P u b l i c   L i c e n s e   f o r   m o r e   d e t a i l s . 
 
   * 
 
   *   Y o u   s h o u l d   h a v e   r e c e i v e d   a   c o p y   o f   t h e   G N U   G e n e r a l   P u b l i c   L i c e n s e 
 
   *   a l o n g   w i t h   t h i s   p r o g r a m .     I f   n o t ,   s e e   < h t t p : / / w w w . g n u . o r g / l i c e n s e s / > . 
 
   *   
 
   *   F o r   f u r t h e r   i n f o r m a t i o n   p l e a s e   c o n t a c t . 
 
   * 	 h t t p : / / n y a t l a . j p / n y a t o o l k i t / 
 
   * 	 < a i r m a i l ( a t ) e b o n y . p l a l a . o r . j p >   o r   < n y a t l a ( a t ) n y a t l a . j p > 
 
   *   
 
   * / 
 
 n a m e s p a c e   j p . n y a t l a . n y a r t o o l k i t . c s . c o r e 
 
 { 
 
 
 
 
 
         / * * 
 
           *   S0n0000o00$R%Rl0(uD0f0weE\$P0Bl00_j0cOW0~0Y00
 
           *   < p > m e m o : ;u }peL02 0 4 8 ^ 2 k0n_0j0D04XTo00f i x e d i n t ( 2 4 - 8 ) g0{g0M0~0Y00< / p > 
 
           * / 
 
         p u b l i c   c l a s s   N y A R H i s t o g r a m A n a l y z e r _ D i s c r i m i n a n t T h r e s h o l d   :   I N y A R H i s t o g r a m A n a l y z e r _ T h r e s h o l d 
 
         { 
 
                 p r i v a t e   d o u b l e   _ s c o r e ; 
 
                 / * * 
 
                   *   S0n0peo00$R%Rl0(uD0f0weE\$P01 PBl0~0Y00weE\$Pn0{Vo00i _ h i s t o g r a m _pen0{Vh0TX0g0Y00
 
                   *   peo00t h i s n0000000feW0~0Y00
 
                   * / 
 
                 p u b l i c   i n t   g e t T h r e s h o l d ( N y A R H i s t o g r a m   i _ h i s t o g r a m ) 
 
                 { 
 
                         i n t [ ]   h i s t   =   i _ h i s t o g r a m . d a t a ; 
 
                         i n t   n   =   i _ h i s t o g r a m . l e n g t h ; 
 
                         i n t   d a ,   s a ,   d b ,   s b ,   d t ,   p t ,   s t ; 
 
                         i n t   i ; 
 
                         i n t   t h   =   0 ; 
 
                         / / _g0OF0
 
                         d t   =   p t   =   0 ; 
 
                         f o r   ( i   =   0 ;   i   <   n ;   i + + ) 
 
                         { 
 
                                 i n t   h   =   h i s t [ i ] ; 
 
                                 d t   + =   h   *   i ; 
 
                                 p t   + =   h   *   i   *   i ; / / ckSn0Bfk0OF0[pe
 
                         } 
 
                         s t   =   i _ h i s t o g r a m . t o t a l _ o f _ d a t a ; 
 
                         / / L o w tP( 0 < = i < = n - 2 ) 
 
                         d a   =   d t ; 
 
                         s a   =   s t ; 
 
                         / / H i g h tP( i = n - 1 ) 
 
                         d b   =   s b   =   0 ; 
 
 
 
                         d o u b l e   m a x   =   - 1 ; 
 
                         d o u b l e   m a x _ m t   =   0 ; 
 
                         / / T000000n0R▦^0{Y00( 1 < = i < = n - 1 n0{Vg0UO) 
 
                         f o r   ( i   =   n   -   1 ;   i   >   0 ;   i - - ) 
 
                         { 
 
                                 / / !kn00000000{
 
                                 i n t   h i s t _ c o u n t   =   h i s t [ i ] ; 
 
                                 i n t   h i s t _ v a l   =   h i s t _ c o u n t   *   i ; 
 
                                 d a   - =   h i s t _ v a l ; 
 
                                 s a   - =   h i s t _ c o u n t ; 
 
                                 d b   + =   h i s t _ v a l ; 
 
                                 s b   + =   h i s t _ c o u n t ; 
 
 
 
                                 / / 000Rce0{
 
                                 d o u b l e   d v   =   ( s a   +   s b ) ; 
 
                                 d o u b l e   m t   =   ( d o u b l e ) ( d a   +   d b )   /   d v ; 
 
                                 d o u b l e   m a   =   ( s a   ! =   0   ?   ( ( d o u b l e ) d a   /   ( d o u b l e ) s a )   :   0 )   -   m t ; 
 
                                 d o u b l e   m b   =   ( s b   ! =   0   ?   ( ( d o u b l e ) d b   /   ( d o u b l e ) s b )   :   0 )   -   m t ; 
 
                                 d o u b l e   k a i   =   ( ( d o u b l e ) ( s a   *   ( m a   *   m a )   +   s b   *   ( m b   *   m b ) ) )   /   d v ; 
 
                                 i f   ( m a x   <   k a i ) 
 
                                 { 
 
                                         m a x _ m t   =   m t ; 
 
                                         m a x   =   k a i ; 
 
                                         t h   =   i ; 
 
                                 } 
 
                                 / / S y s t e m . o u t . p r i n t l n ( k a i ) ; 
 
                         } 
 
                         / / m a x _ m t 0CQk0ckS
 
                         t h i s . _ s c o r e   =   m a x   /   ( ( d o u b l e ) ( p t   +   m a x _ m t   *   m a x _ m t   *   s t   -   2   *   m a x _ m t   *   d t )   /   s t ) ; / / 1 2 9 , 0 . 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 
 
                         r e t u r n   t h ; 
 
                 } 
 
                 / * * 
 
                   *    g_k0[LW0_0{ @ l i n k   # g e t T h r e s h o l d } n00000ԏtSW0~0Y00
 
                   *   @ r e t u r n 
 
                   *   000$P0{Vo00 & l t ; = n & l t ; = 1 . 0 n0g0Y00
 
                   * / 
 
                 p u b l i c   d o u b l e   g e t L a s t S c o r e ( ) 
 
                 { 
 
                         r e t u r n   t h i s . _ s c o r e ; 
 
                 } 
 
                 / * * 
 
                   *   0000(upe
 
                   *   @ p a r a m   a r g s 
 
                   *   m a i n pe_pe
 
                   * / 
 
                 p u b l i c   s t a t i c   v o i d   m a i n ( s t r i n g [ ]   a r g s ) 
 
                 { 
 
                         N y A R H i s t o g r a m   d a t a   =   n e w   N y A R H i s t o g r a m ( 2 5 6 ) ; 
 
                         f o r   ( i n t   i   =   0 ;   i   <   2 5 6 ;   i + + ) 
 
                         { 
 
                                 d a t a . d a t a [ i ]   =   1 2 8   -   i   >   0   ?   1 2 8   -   i   :   i   -   1 2 8 ; 
 
                         } 
 
                         d a t a . t o t a l _ o f _ d a t a   =   d a t a . g e t T o t a l ( 0 ,   2 5 5 ) ; 
 
                         N y A R H i s t o g r a m A n a l y z e r _ D i s c r i m i n a n t T h r e s h o l d   a n   =   n e w   N y A R H i s t o g r a m A n a l y z e r _ D i s c r i m i n a n t T h r e s h o l d ( ) ; 
 
                         i n t   t h   =   a n . g e t T h r e s h o l d ( d a t a ) ; 
 
                         / / S y s t e m . o u t . p r i n t ( t h ) ; 
 
                         r e t u r n ; 
 
                 } 
 
         } 
 
 } 
 
 