# Patent

## Title
Image processing method using foreground probability

## Abstract

An image processing method of separating an input image into a foreground image and a background image, the method including determining a pixel of the input image as a pixel of the foreground image if a foreground probability value of the pixel of the foreground image determined by using the Gaussian mixture model or the pixel determined to be included in a motion region is greater than a setting threshold.

## Background
CROSS-REFERENCE TO RELATED PATENT APPLICATION
This application claims priority from Korean Patent Application No. 10-2012-0025660, filed on Mar. 13, 2012, in the Korean Intellectual Property Office, the disclosure of which is incorporated herein in its entirety by reference.
BACKGROUND
1. Field
Methods consistent with exemplary embodiments relate to image processing, and more particularly, to receiving an image from an image capturing device such as a surveillance camera for capturing an inside of an elevator and separating the received image into a foreground image and a background image.
2. Description of the Related Art
Surveillance systems need to separate a received image into a foreground image and a background image to perform various surveillance functions. That is, host apparatuses receive an image from a surveillance camera, determine whether each pixel of the received image is a pixel of a foreground image or a pixel of a background image, and separate the received image into the foreground image and the background image.
To separate the received image into the foreground image and the background image, conventionally, a Gaussian model was used to determine a pixel of the foreground image or a pixel of the background image (see Korean Patent Registration No. 1038650).
Such an image separation method is not robust and has a low separation accuracy with respect to an environment including a reflector and a frequent change in illumination, such as an inside of an elevator.
SUMMARY
One or more exemplary embodiments provide an image processing method of separating a received image into a foreground image and a background image, that is robust and enhances separation accuracy with respect to an environment including a reflector and a frequent change in illumination, such as the inside of an elevator.
According to an aspect of an exemplary embodiment, there is provided an image processing method including: (a) determining whether a pixel of an input image is a pixel of a foreground image or a pixel of the background image by using a Gaussian mixture model; (b) determining whether the pixel of the input image is included in a motion region; (c) obtaining a foreground probability value of the pixel of the input image according to a foreground probability histogram with respect to a correlation between the input image and a reference image of the input image; and (d) determining the pixel of the input image as the pixel of the foreground image if the foreground probability value of the pixel of the foreground image determined by using the Gaussian mixture model or the pixel determined to be included in the motion region is greater than a setting threshold.
The input image may include a series of frame images of video, wherein operations (a) through (d) are repeatedly performed.
In operation (c), the reference image may be separated into a foreground and a background, and renewed as a previous image during a process of repeatedly performing operations (a) through (d).
In operation (c), the correlation between the input image and the reference image of the input image may be calculated with respect to the pixel of the input image by using a texture value of the pixel that is a difference between a mean gradation value of the pixel and neighboring pixels and gradation values of the neighboring pixels.
The input image may be an image of an inside of a closed space captured by an image obtaining unit.
The method may further include: (e), if the pixel determined as the pixel of the foreground image is not connected to an image of a floor of the closed space, setting the pixel determined as the pixel of the foreground image as a pixel of the background image.
In operation (c), the correlation between the input image and the reference image may be calculated with respect to the pixel of the input image by using a texture value of the pixel that is a difference between a mean gradation value of the pixel and neighboring pixels and gradation values of the neighboring pixels.
The input image may include a series of frame images of video, wherein operations (a) through (d) are repeatedly performed.
In operation (c), the reference image may be separated into a foreground and a background, and renewed as a previous image during a process of repeatedly performing operations (a) through (d).
According to an aspect of another exemplary embodiment, there is provided an image processing method of receiving an image from a surveillance camera for capturing an inside of a closed space having a door to an outside, determining whether each pixel of the received image is a pixel of a foreground image or a pixel of a background image, and separating the received image into the foreground image and the background image, the method including: (a) determining whether a pixel of the received image is a pixel of a foreground image or a pixel of a background image by using a Gaussian mixture model; (b) determining whether the pixel of the received image is included in a motion region; (c) obtaining a foreground probability value of the pixel of the received image according to a foreground probability histogram with respect to correlation between the received image and a reference image of the received image; (d) if the foreground probability value of the pixel of the foreground image determined by using the Gaussian mixture model or the pixel determined to be included in the motion region is greater than a setting threshold, obtaining a foreground-information binary number “1” as a first foreground determination value of the pixel of the received image, and, if the foreground probability value is not greater than the setting threshold, obtaining a background-information binary number “0” as the first foreground determination value; (e) determining whether the door of the closed space is open or closed; (f) if the door of the closed space is open, setting an auxiliary foreground determination value in the same way as the first foreground determination value, and, if the door of the elevator is closed, obtaining a final foreground determination value of a previous image and the auxiliary foreground determination value set based on the first foreground determination value; and (g) if a result obtained by multiplying the auxiliary foreground determination value and the foreground probability value is greater than a final threshold, obtaining a final foreground determination value as the foreground-information binary number “1”, and, if the result is not greater than the final threshold, obtaining the final foreground determination value as the background-information binary number “0”.
The method may further include: (h) if the pixel determined as the pixel of the foreground image is not connected to an image of a floor of the closed space, setting the pixel determined as the pixel of the foreground image as the pixel of the background image.
The received image may include a series of frame images of video, wherein operations (a) through (g) are repeatedly performed.
In operation (c), the reference image may be separated into a foreground and a background, and renewed as a previous image during a process of repeatedly performing operations (a) through (g).
Operation (e) may include: if a mean foreground probability value of pixels corresponding to top portions of the door of the closed space is greater than a setting mean threshold, determining the door of the closed space to be open.
In operation (c), the correlation between the received image and the reference image of the received image may be calculated with respect to the pixel of the received image by using a texture value of the pixel that is a difference between a mean gradation value of the pixel and neighboring pixels and gradation values of the neighboring pixels.



BRIEF DESCRIPTION OF THE DRAWINGS
The above and other aspects will become more apparent by describing in detail exemplary embodiments with reference to the attached drawings, in which:
FIG. 1 is a diagram for explaining a received image to which an image processing method is to be applied, according to an exemplary embodiment;
FIG. 2 is a flowchart of an image processing method, according to an exemplary embodiment;
FIG. 3 is a diagram for explaining an operation of calculating a foreground probability value of FIG. 2, according to an exemplary embodiment;
FIG. 4 shows an example of a histogram of a foreground probability with respect to a correlation between a reference image and a received image;
FIG. 5 shows an example of a histogram of a background probability with respect to a correlation between a reference image and a received image;
FIG. 6 is a diagram for explaining the image processing method of FIG. 2, according to an exemplary embodiment;
FIG. 7 is a flowchart of an image processing method, according to another exemplary embodiment;
FIG. 8 is a block flowchart for explaining the image processing method of FIG. 7, according to an exemplary embodiment;
FIG. 9 is a screen image block diagram for explaining a determination operation of FIG. 8, according to an exemplary embodiment;
FIG. 10 is a screen image for explaining some operations of the image processing method of FIG. 7, according to an exemplary embodiment; and
FIG. 11 is a screen image for explaining a post-processing operation of FIG. 8, according to an exemplary embodiment.



DETAILED DESCRIPTION OF EXEMPLARY EMBODIMENTS
The inventive concept will now be described more fully with reference to the accompanying drawings, in which exemplary embodiments of the invention are shown. The inventive concept may, however, be embodied in many different forms and should not be construed as being limited to the embodiments set forth herein; rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the concept of the inventive concept to those of ordinary skill in the art.
FIG. 1 is a diagram for explaining a received image XT to which an image processing method is to be applied, according to an exemplary embodiment.
Referring to FIG. 1, the received image XT includes images x(t), . . . , x(t−T) of a series of frames, for example, F1-F15, of video input during a setting period T. In this regard, t denotes time. All of the images x(t), . . . , x(t−T) of the frames F1-F15 may be used or a representative image thereof may be used.
FIG. 2 is a flowchart of an image processing method, according to an exemplary embodiment. FIG. 3 is a diagram for explaining an operation S23 of calculating a foreground probability value PFG(fi|xi) of FIG. 2. FIG. 4 shows an example of a histogram of a foreground probability with respect to a correlation fi between a reference image 31 and a received image 32 in FIG. 3 or 61 in FIG. 6. FIG. 5 shows an example of a histogram of a background probability with respect to the correlation f, between the reference image 31 and the received image 32 or 61. FIG. 6 is a diagram for explaining the image processing method of FIG. 2.
The same reference numerals of FIGS. 1 through 6 denote the same elements.
Referring to FIGS. 1 through 6, the image processing method of FIG. 2 receives an image 61 from a surveillance camera for capturing, for example, an inside of an elevator, determines whether each pixel of the received image 61 is a pixel of a foreground image or a pixel of a background image, and separates the received image 61 into the foreground image and the background image. The image processing method of FIG. 2 includes operations S21 through S26.
In operation S21, a Gaussian mixture model is used to determine whether each pixel of the received image 61 is the pixel of the foreground image or the pixel of the background image.
In more detail, in a case where an M number of image frames are input during the setting period T, a Gaussian mixture distribution of a pixel at a time t may be expressed according to Equation 1 below.




  
    
      
        
          p
          ⁡
          
            (
            
              
                
                  x
                  _
                
                |
                
                  X
                  T
                
              
              ,
              
                BG
                +
                FG
              
            
            )
          
        
        =
        
          
            ∑
            
              m
              =
              1
            
            M
          
          ⁢
          
            
          
          ⁢
          
            
              π
              m
            
            ⁢
            
              N
              ⁡
              
                (
                
                  
                    x
                    _
                  
                  ,
                  
                    
                      
                        μ
                        _
                      
                      ^
                    
                    m
                  
                  ,
                  
                    
                      
                        
                          σ
                          _
                        
                        ^
                      
                      m
                      2
                    
                    ⁢
                    I
                  
                
                )
              
            
          
        
      
    
    
      
        [
        
          Equation
          ⁢
          
            
          
          ⁢
          1
        
        ]
      
    
  




In Equation 1 above, p denotes a probability, x denotes gradation of a pixel, XT denotes the received image, BG denotes a pixel of a background image probability, FG denotes a pixel of a background image probability, πm denotes a weighted value of an mth Gaussian distribution, N denotes a normal distribution, μm denotes a mean value of the mth Gaussian distribution, σm2 denotes a variance value x(t)-μm of the mth Gaussian distribution, and I denotes a constant.
As is well known, in a case where gradation of a pixel is input at a next time t+1, a Gaussian distribution is recursively renewed according to Equations 2 through 4 below.

{circumflex over (π)}m←{circumflex over (π)}m+α(om(t)−{circumflex over (π)}m)  [Equation 2]

 {circumflex over (μ)}m← {circumflex over (μ)}m+om(t)(α/{circumflex over (π)}m){circumflex over ({right arrow over (o)})}m  [Equation 3]

{circumflex over (σ)}m2←{circumflex over (σ)}m2+om(t)(a/{circumflex over (π)}m)({right arrow over (δ)}mT{right arrow over (δ)}m−{circumflex over (σ)}m2)  [Equation 4]

In Equations 2 through 4 above, α denotes a setting constant. om has a value “1” in a case where the gradation of the pixel at the time t matches the Gaussian distribution, and has a value “0” in a case where the gradation of the pixel at the time t does not match the Gaussian distribution.
In this regard, a Gaussian distribution of gradation of the pixel corresponding to the background image has a large weighted value πm and a small variance value σm2, compared to a Gaussian distribution of gradation of the pixel corresponding to the foreground image.
Such characteristics may be used to obtain a background probability of the pixel by summing a B number of weighted values πm in a sequence of ascending weighted values πm. That is, in a case where a value of the background probability is greater than a setting threshold, the pixel is determined as a pixel of the background image, and, in a case where the value of the background probability is not greater than the setting threshold, the pixel is determined as a pixel of the foreground image according to Equation 5 below.




  
    
      
        
          p
          ⁡
          
            (
            
              
                
                  x
                  _
                
                |
                
                  X
                  T
                
              
              ,
              BG
            
            )
          
        
        ≅
        
          
            ∑
            
              m
              =
              1
            
            B
          
          ⁢
          
            
          
          ⁢
          
            
              
                π
                ^
              
              m
            
            ⁢
            
              N
              ⁡
              
                (
                
                  
                    x
                    |
                    
                      
                        
                          μ
                          _
                        
                        ^
                      
                      m
                    
                  
                  ,
                  
                    
                      σ
                      m
                      2
                    
                    ⁢
                    I
                  
                
                )
              
            
          
        
      
    
    
      
        [
        
          Equation
          ⁢
          
            
          
          ⁢
          5
        
        ]
      
    
  




In operation S22, it is determined whether the received image 61 is included in a motion region. That is, the motion region is separated from the received image 61. In addition to an algorithm for detecting the motion region that uses optical flow, various other algorithms are well known, and thus, detailed descriptions thereof are omitted here.
In operation S23, the foreground probability value PFG(fi|xi) of each of pixels of the received images 32 or 61 is calculated according to the foreground probability histogram (FIG. 4) regarding the correlation fi between the received images 32 or 61 and the reference image 31 of the received images 32 and 61.
In the present embodiment, the reference image 31 is separated into a foreground and a background, and is renewed (operation S25) to a previous image during a process of repeatedly performing operations S21 through S24 until an end signal is generated (operation S26). Thus, accuracy of the foreground probability value PFG(fi|xi) increases.
In operation S23, a texture value of a target pixel that is a difference between a mean gradation value ( mi, ūi) of the target pixel and eight neighboring pixels and gradation values mj, uj of the eight neighboring pixels is used to calculate the correlation fi between the received images 32 or 61 and the reference image 31 of the received images 32 or 61. Accordingly, the image processing state may be made further robust and separation accuracy may be enhanced.
The correlation fi between the received images 32 or 61 and the reference image 31 of the received images 32 or 61 includes a normalized cross correlation (NCC) function fi1 and a texture feature function fi2 of Equations 6 and 7, respectively, below.




  
    
      
        
          f
          i
          1
        
        =
        
          
            
              ∑
              
                j
                ∈
                
                  ω
                  1
                
              
            
            ⁢
            
              
            
            ⁢
            
              
                (
                
                  
                    u
                    j
                  
                  -
                  
                    
                      u
                      _
                    
                    i
                  
                
                )
              
              ⁢
              
                (
                
                  
                    m
                    j
                  
                  -
                  
                    
                      m
                      _
                    
                    i
                  
                
                )
              
            
          
          
            
              
                ∑
                
                  j
                  ∈
                  
                    ω
                    1
                  
                
              
              ⁢
              
                
              
              ⁢
              
                
                  
                    (
                    
                      
                        u
                        j
                      
                      -
                      
                        
                          u
                          _
                        
                        i
                      
                    
                    )
                  
                  2
                
                ⁢
                
                  
                    ∑
                    
                      j
                      ∈
                      
                        ω
                        1
                      
                    
                  
                  ⁢
                  
                    
                  
                  ⁢
                  
                    
                      (
                      
                        
                          m
                          j
                        
                        -
                        
                          
                            m
                            _
                          
                          i
                        
                      
                      )
                    
                    2
                  
                
              
            
          
        
      
    
    
      
        [
        
          Equation
          ⁢
          
            
          
          ⁢
          6
        
        ]
      
    
  
  
    
      
        
          f
          i
          2
        
        =
        
          
            
              
                ∑
                
                  j
                  ∈
                  
                    ω
                    1
                  
                
              
              ⁢
              
                
              
              ⁢
              
                
                  (
                  
                    
                      u
                      j
                    
                    -
                    
                      
                        u
                        _
                      
                      i
                    
                  
                  )
                
                2
              
            
          
          +
          
            
              
                ∑
                
                  j
                  ∈
                  
                    ω
                    1
                  
                
              
              ⁢
              
                
              
              ⁢
              
                
                  (
                  
                    
                      m
                      j
                    
                    -
                    
                      
                        m
                        _
                      
                      i
                    
                  
                  )
                
                2
              
            
          
        
      
    
    
      
        [
        
          Equation
          ⁢
          
            
          
          ⁢
          7
        
        ]
      
    
  




In Equations 6 and 7 above, i denotes the target pixel, j denotes neighboring pixels, ω denotes a window, mi denotes a mean gradation value of a window 31i of the reference image 31, ūi denotes a mean gradation value of a window 32i of the received image 32, mj denotes a neighboring gradation value of the window 31i of the reference image 31, and uj denotes a neighboring gradation value of the window 32i of the received image 32.
Therefore, the foreground probability value PFG(fi|xi) of the target pixel may be calculated according to the foreground probability histogram (FIG. 4) regarding correlation values obtained by using Equations 6 and 7 above.
In operation S24, each of the pixels of the received images 32 or 61 is determined as a pixel of the foreground image if the foreground probability value PFG(fi|xi) of the pixel of the foreground image determined by using the Gaussian mixture model (operation S21) or the pixel determined to be included in the motion region (operation S22) is greater than the setting threshold Tm. Therefore, if information regarding a foreground image of a binary number “0” or “1” determined by using the Gaussian mixture model (operation S21) is Mi,GMM, and information regarding a motion region of the binary number “0” or “1” is Di with respect to an ith pixel, a foreground determination value Mi, motion may be obtained according to Equation 3 below.




  
    
      
        
          M
          
            i
            ,
            motion
          
        
        =
        
          {
          
            
              
                1
              
              
                
                  
                    if
                    ⁢
                    
                      
                    
                    ⁢
                    
                      
                        (
                        
                          
                            M
                            
                              i
                              ,
                              GMM
                            
                          
                          ⋁
                          
                            D
                            i
                          
                        
                        )
                      
                      ·
                      
                        
                          p
                          FG
                        
                        ⁡
                        
                          (
                          
                            
                              f
                              i
                            
                            |
                            
                              x
                              i
                            
                          
                          )
                        
                      
                    
                  
                  ≥
                  
                    T
                    m
                  
                
              
            
            
              
                0
              
              
                otherwise
              
            
          
        
      
    
    
      
        [
        
          Equation
          ⁢
          
            
          
          ⁢
          8
        
        ]
      
    
  




In Equation 8 above, an OR operation is performed on the information Mi,GMM regarding the foreground image determined by using the Gaussian mixture model and the information Di regarding the motion region (operation S24a of FIG. 6). Then, multiplication is performed on a result of the OR operation and the foreground probability value PFG(fi|xi), and an operation of determining whether a pixel is a foreground pixel is performed on a result of the multiplication (operation S24b of FIG. 6).
Therefore, according to an embodiment described with reference to FIGS. 1 through 6, the image processing method may be robust and enhance separation accuracy with respect to an environment including a reflector and a frequent change in illumination such as the inside of an elevator.
In addition, in a case where Equation 8 above is applied to each pixel of images of the inside of the elevator, if a determined foreground image is not connected to an image of a floor of the elevator, the determined foreground image is set as the background image. That is, characteristics of the elevator may be used to correct a separation error due to the reflector. This will now be described in detail with reference to FIGS. 7 through 11.
FIG. 7 is a flowchart of an image processing method, according to another exemplary embodiment. FIG. 8 is a block flowchart for explaining the image processing method of FIG. 7. FIG. 9 is a screen image block diagram for explaining a determination operation S82 of FIG. 8. FIG. 10 is a screen image for explaining operations S705 through S707 of FIG. 7. FIG. 11 is a screen image for explaining a post-processing operation S83 of FIG. 8. The post-processing operation S83 of FIG. 8 includes operations S709 and S710 of FIG. 7.
The same reference numerals of FIGS. 7 through 11 denote the same elements.
Referring to FIGS. 7 through 11, the image processing method of FIG. 7 receives an image 81, 91, 101a, 102a or 111 from a surveillance camera for capturing an inside of an elevator, determines whether each pixel of the received image 81, 91, 101a, 102a or 111 is a pixel of a foreground image or a pixel of a background image, and separates the received image 81, 91, 101a, 102a or 111 into a foreground image and a background image, and includes operations S701 through S712.
In operation S701, a Gaussian mixture model is used to determine whether each pixel of the received images 81, 91, 101a, 102a or 111 is a pixel of the foreground image or a pixel of the background image. Operation S701 is the same as described with reference to S21 of FIG. 2.
In operation S702, it is determined whether each of the received images 81, 91, 101a, 102a, or 111 is included in a motion region. That is, the motion region is separated from the received image 81, 91, 101a, 102a or 111. In addition to an algorithm for detecting the motion region that uses optical flow, various other algorithms are well known, and thus, detailed descriptions thereof are omitted here.
In operation S703, the foreground probability value PFG(fi|xi) of each of pixels of the received image 81, 91, 101a, 102a or 111 is calculated according to the foreground probability histogram (FIG. 4) regarding the correlation fi between the received image 81, 91, 101a, 102a or 111 and the reference image 31 of the received images 81, 91, 101a, 102a or 111.
In the present embodiment, the reference image 31 is separated into a foreground and a background and is renewed (operation S711) to a previous image during a process of repeatedly performing operations S701 through S712 until an end signal is generated (operation S712). Thus, accuracy of the foreground probability value PFG(fi|xi) increases. Operation S703 is the same as described with reference to S23 of FIG. 2.
In operation S704, a first foreground determination value Mi(t), motion is obtained. That is, if the foreground probability value PFG(fi|xi) of the pixel of the foreground image determined by using the Gaussian mixture model (operation S21) or the pixel determined to be included in the motion region (operation S22) is greater than the setting threshold Tm, a foreground-information binary number “1” is obtained as a first foreground determination value Mi(t), motion of each of pixels of the received images 81, 91, 101a, 102a or 111, and, if the foreground probability value PFG(fi|xi) is not greater than the setting threshold Tm, a background-information binary number “0” is obtained as the first foreground determination value Mi(t), motion.
Operation S704 is the same as described with reference to S24 of FIG. 2. That is, although Mi, motion denotes the foreground determination value in Equation 8 above, Mi(t), motion denotes the foreground determination value at the time t that is a time at which a currently received image is processed, i.e., the first foreground determination value.
Operations S701 through S704 above correspond to a first determination operation (S91 of FIG. 9) of producing a preliminary separation result (93 of FIG. 9). Also, operations S705 through S708 below correspond to a determination operation (S92 of FIG. 9) of using the preliminary separation result (93 of FIG. 9) and information regarding a previous image (92 of FIG. 9).
Thus, the determination operation S82 of FIG. 8 includes the first determination operation (S91 of FIG. 9) and the determination operation (S92 of FIG. 9).
Operations S709 and S710 below correspond to the post-processing operation S83 of FIG. 8.
Operations S705 through S710 will now be described in detail.
In operation S705, it is determined whether a door of the elevator is open or closed.
In this regard, if a mean foreground probability value PFG—MEAN of pixels corresponding to top portions Aup1 and Aup2 of FIG. 10, of the door of the elevator is greater than a setting mean threshold, the door of the elevator is determined to be open. That is, images of the top portions Aup1 and Aup2 of the door of the elevator scarcely change due to presence of people, and thus, the mean foreground probability value PFG—MEAN is low when the door of the elevator is closed, and the mean foreground probability value PFG—MEAN is high when the door of the elevator is open, thereby accurately determining whether the door of the elevator is open or closed.
In operation S706, which is performed when the door of the elevator is open, an auxiliary foreground determination value Q(t) is set in the same way as a first foreground determination value Mi(t), motion.
In operation S707 performed when the door of the elevator is closed, a final foreground determination value Mi(t−1), final of the previous image and the auxiliary foreground determination value Q(t) set based on the first foreground determination value Mi(t), motion are obtained.
That is, operations S706 and S707 above may be performed according to Equation 9 below.

Q(t)=[(1−λ)Mi,motion(t)+λMi,final(t−1)]  [Equation 9]

In Equation 9 above, a parameter λ has a value “0” when the door of the elevator is open. Thus, in operation S706 performed when the door of the elevator is open, the auxiliary foreground determination value Q(t) is set in the same way as a first foreground determination value Mi(t), motion (see 101a and 101b of FIG. 10).
Also, in Equation 9 above, the parameter λ has a value λclose greater than “0” and smaller than “1” when the door of the elevator is closed. According to an experiment of the present embodiment, the parameter λ may have about “0.65” when the door of the elevator is closed. Thus, in operation S707 performed when the door of the elevator is closed, the final foreground determination value Mi(t−1), final of the previous image and the auxiliary foreground determination value Q(t) set based on the first foreground determination value Mi(t), motion are obtained (see 102a and 102b of FIG. 10).
Next, in operation S708, if a result obtained by multiplying the auxiliary foreground determination value Q(t) and the foreground probability value PFG(fi|xi) is greater than a final threshold Tb, a foreground-information binary number “1” is obtained as a final foreground determination value Mi(t), final, and, if the result is not greater than the final threshold Tb, a background-information binary number “0” is obtained as the final foreground determination value Mi(t), final. That is, operation S708 may be expressed according to Equation 10 below.




  
    
      
        
          M
          
            i
            ,
            final
          
          
            (
            t
            )
          
        
        =
        
          {
          
            
              
                1
              
              
                
                  
                    if
                    ⁢
                    
                      
                    
                    ⁢
                    
                      
                        Q
                        
                          (
                          t
                          )
                        
                      
                      ·
                      
                        
                          p
                          FG
                        
                        ⁡
                        
                          (
                          
                            
                              f
                              i
                            
                            |
                            
                              x
                              i
                            
                          
                          )
                        
                      
                    
                  
                  ≥
                  
                    T
                    b
                  
                
              
            
            
              
                0
              
              
                otherwise
              
            
          
        
      
    
    
      
        [
        
          Equation
          ⁢
          
            
          
          ⁢
          10
        
        ]
      
    
  




In summary, since a subject hardly moves when the door of the elevator is closed, accuracy of the first foreground determination value Mi(t), motion may deteriorate. In this case, the final foreground determination value Mi(t−1), final of the previous image and the auxiliary foreground determination value Q(t) set based on the first foreground determination value Mi(t), motion are used, thereby preventing a separation error from occurring when the door of the elevator is closed.
Operations S709 and S710 below correspond to the post-processing operation (S83 of FIG. 8). In the post-processing operation (S83 of FIG. 8), small foreground errors in drop shapes may be removed by using a well known morphological image process before operation S709 below is performed.
That is, an image 112 of FIG. 11 indicating a first separation result is obtained by performing Equation 10. Also, an image 113 of FIG. 11 indicates a result obtained by removing small foreground errors in drop shapes by using the well known morphological image process. In this regard, an image 114 indicating a final separation result may be obtained by performing operations S709 and S710 below.
In operation S709, it is determined whether a determined foreground image is connected to an image Ifl of FIG. 11 of a floor of the elevator.
If a determined foreground image is not connected to the image Ifl of FIG. 11 of the floor of the elevator, the determined foreground image is set as the background image (operation S710).
Accordingly, the image 114 indicating the final separation result is obtained. That is, the separation error due to a reflector may be corrected by using characteristics of the elevator.
As described above, the image processing methods according to the exemplary embodiments obtain a foreground probability value of each pixel of a received image according to a foreground probability histogram of a correlation between the received image and a reference image of the received image. Also, each pixel of the received image is determined as a pixel of a foreground image if the foreground probability value of a pixel of a foreground image determined by using a Gaussian mixture model or a pixel determined to be included in a motion region is greater than a setting threshold.
Therefore, the image processing methods may be robust and enhance separation accuracy with respect to an environment including a reflector and a frequent change in illumination such as the inside of an elevator.
Furthermore, the reference image is renewed to a previous image including foreground and background information, and thus, accuracy of the foreground probability value may increase.
Furthermore, a texture value of a target pixel that is a difference between a mean gradation value of the target pixel and neighboring pixels and gradation values of the neighboring pixels is used to calculate the correlation between the received image and the reference image of the received image with respect to each pixel of the received image. Thus, the image processing methods may be robust and further enhance separation accuracy.
In a case where a determined foreground image is not connected to an image of a floor of the elevator, the determined foreground image is set as a background image. That is, a separation error due to a reflector may be corrected by using characteristics of the elevator.
Meanwhile, when a door of the elevator is open, the auxiliary foreground determination value Q(t) is set in the same way as a first foreground determination value Mi(t), motion, and when the door of the elevator is closed, the final foreground determination value Mi(t−1), final of the previous image and the auxiliary foreground determination value Q(t) set based on the first foreground determination value Mi(t), motion are obtained. If a result obtained by multiplying the auxiliary foreground determination value Q(t) and the foreground probability value PFG(fi|xi) is greater than the final threshold Tb, the foreground-information binary number “1” is obtained as a final foreground determination value Mi(t), final, and, if the result is not greater than the final threshold Tb, the background-information binary number “0” is obtained as the final foreground determination value Mi(t), final.
In this regard, since a subject hardly moves when the door of the elevator is closed, accuracy of the first foreground determination value Mi(t), motion may deteriorate. In this case, the final foreground determination value Mi(t−1), final of the previous image and the auxiliary foreground determination value Q(t) set based on the first foreground determination value Mi(t), motion are used, thereby preventing a separation error from occurring when the door of the elevator is closed.
Meanwhile, if the mean foreground probability value PFG—MEAN of pixels corresponding to top portions of the door of the elevator is greater than the setting mean threshold, the door of the elevator is determined to be open. That is, images of the top portions of the door of the elevator scarcely change due to the presence of people, and thus the mean foreground probability value PFG—MEAN is low when the door of the elevator is closed, and the mean foreground probability value PFG—MEAN is high when the door of the elevator is open, thereby accurately determining whether the door of the elevator is open or closed.
While the inventive concept has been particularly shown and described with reference to exemplary embodiments thereof, it will be understood by those of ordinary skill in the art that various changes in form and details may be made therein without departing from the spirit and scope of the inventive concept as defined by the following claims.

## Claims

1. An image processing method comprising:
(a) determining whether a pixel of an input image is a pixel of a foreground image or a pixel of a background image by using a Gaussian mixture model;
(b) determining whether the pixel of the input image is included in a motion region;
(c) obtaining a foreground probability value of the pixel of the input image according to a foreground probability histogram with respect to a correlation between the input image and a reference image of the input image; and
(d) determining the pixel of the input image as the pixel of the foreground image if the foreground probability value of the pixel of the foreground image determined as the pixel of the foreground image at operation (a) or the foreground probability value of the pixel determined to be included in the motion region at operation (b) is greater than a setting threshold.

2. The method of claim 1, wherein the input image comprises a series of frame images of video, and
wherein operations (a) through (d) are repeatedly performed.

3. The method of claim 2, wherein, in operation (c), the reference image is separated into a foreground and a background, and is renewed as a previous image during a process of repeatedly performing operations (a) through (d).

4. The method of claim 1, wherein, in operation (c), the correlation between the input image and the reference image is calculated with respect to the pixel of the input image by using a texture value of the pixel that is a difference between a mean gradation value of the pixel and neighboring pixels and gradation values of the neighboring pixels.

5. The method of claim 1, wherein the input image is an image of an inside of a closed space captured by an image obtaining unit.

6. The method of claim 5, further comprising: (e), if the pixel determined as the pixel of the foreground image is not connected to an image of a floor of the closed space, setting the pixel determined as the pixel of the foreground image as a pixel of the background image.

7. The method of claim 5, wherein, in operation (c), the correlation between the input image and the reference image is calculated with respect to the pixel of the input image by using a texture value of the pixel that is a difference between a mean gradation value of the pixel and neighboring pixels and gradation values of the neighboring pixels.

8. The method of claim 5, wherein the input image comprises a series of frame images of video, and
wherein operations (a) through (d) are repeatedly performed.

9. The method of claim 8, wherein, in operation (c), the reference image is separated into a foreground and a background, and is renewed as a previous image during a process of repeatedly performing operations (a) through (d).

10. An image processing method of receiving an image from a surveillance camera for capturing an inside of a closed space having a door to an outside, determining whether each pixel of the received image is a pixel of a foreground image or a pixel of a background image, and separating the received image into the foreground image and the background image, the method comprising:
(a) determining whether a pixel of the received image is a pixel of a foreground image or a pixel of a background image by using a Gaussian mixture model;
(b) determining whether the pixel of the received image is included in a motion region;
(c) obtaining a foreground probability value of the pixel of the received image according to a foreground probability histogram with respect to correlation between the received image and a reference image of the received image;
(d) if the foreground probability value of the pixel of the foreground image determined as the pixel of the foreground image or the foreground probability value of the pixel determined to be included in the motion region is greater than a setting threshold, obtaining a foreground-information binary number “1” as a first foreground determination value of the pixel of the received image, and, if the foreground probability value is not greater than the setting threshold, obtaining a background-information binary number “0” as the first foreground determination value;
(e) determining whether the door of the closed space is open or closed;
(f) if the door of the closed space is open, setting an auxiliary foreground determination value in the same way as the first foreground determination value, and, if the door of the elevator is closed, obtaining a final foreground determination value of a previous image and the auxiliary foreground determination value set based on the first foreground determination value; and
(g) if a result obtained by multiplying the auxiliary foreground determination value and the foreground probability value is greater than a final threshold, obtaining a final foreground determination value as the foreground-information binary number “1”, and, if the result is not greater than the final threshold, obtaining the final foreground determination value as the background-information binary number “0”.

11. The method of claim 10, further comprising: (h) if the pixel determined as the pixel of the foreground image is not connected to an image of a floor of the closed space, setting the pixel determined as the pixel of the foreground image as the pixel of the background image.

12. The method of claim 10, wherein the received image comprises a series of frame images of video,
wherein operations (a) through (g) are repeatedly performed.

13. The method of claim 12, wherein, in operation (c), the reference image is separated into a foreground and a background, and is renewed as a previous image during a process of repeatedly performing operations (a) through (g).

14. The method of claim 10, wherein operation (e) comprises: if a mean foreground probability value of pixels corresponding to top portions of the door of the closed space is greater than a setting mean threshold, determining the door of the closed space to be open.

15. The method of claim 10, wherein, in operation (c), the correlation between the received image and the reference image of the received image is calculated with respect to the pixel of the received image by using a texture value of the pixel that is a difference between a mean gradation value of the pixel and neighboring pixels and gradation values of the neighboring pixels.

16. An image processing method comprising:
(a) determining whether a pixel of an input image is a pixel of a foreground image or a pixel of a background image;
(b) determining whether the pixel is included in a motion region of the input image;
(c) obtaining a foreground probability value of the pixel according to a foreground probability histogram with respect to a correlation between the input image and a reference image of the input image; and
(d) determining that the pixel is included in the foreground image if the foreground probability value of the pixel determined as the pixel of the foreground image at operation (a) or the pixel determined to be included in the motion region at operation (b) is greater than a setting threshold.

17. The method of claim 16, wherein operation (a) is performed using a Gaussian mixture model.

18. The method of claim 16, wherein, in operation (c), the correlation between the input image and the reference image is calculated with respect to the pixel of the input image by using a texture value of the pixel that is a difference between a mean gradation value of the pixel and neighboring pixels and gradation values of the neighboring pixels.

19. The method of claim 16, wherein the input image is an image of an inside of a closed space captured by an image obtaining unit, and wherein the method further comprises:
(e) if the pixel determined as the pixel of the foreground image is not connected to an image of a floor of the closed space, setting the pixel determined as the pixel of the foreground image as a pixel of the background image.

20. The method of claim 19, wherein, in operation (c), the correlation between the input image and the reference image is calculated with respect to the pixel of the input image by using a texture value of the pixel that is a difference between a mean gradation value of the pixel and neighboring pixels and gradation values of the neighboring pixels.