Abstract:
The present invention relates to an interpolation method for enlarging a digital image or predicting a moving vector of a compressed image system as a sub-pixel unit when the image digitized through a CCD (Charge Coupled Device) camera ect. has a low resolution in a video phone or video conference or general digital video system, particularly the present invention can be adapted to a post processor of a compressed digital image in order to improve the image quality, and can be used for finding a moving vector of a moving picture compressed type, accordingly the present invention is capable of improving the image quality.

Description:
CROSS-REFERENCE TO RELATED APPLICATIONS 
     NOTICE: More than one reissue application has been filed for the reissue of U.S. Pat. No. 6,803,954, including reissue application Ser. No. 12/424,927; U.S. patent application Ser. No. 12/024,408, which is a divisional reissue of U.S. Pat. No. 6,803,954; and U.S. patent application Ser. No. 11/546,484, which is a reissue of U.S. Pat. No. 6,803,954. Accordingly, this application is a continuation application of U.S. patent application Ser. No. 11/546,484, filed Oct. 12, 2006, now U.S. Pat. No. Re. 42,045, which is a reissue of U.S. application Ser. No. 09/692,156, filed Oct. 20, 2000, now U.S. Pat. No. 6,803,954, which claims the benefit of a foreign priority application filed in Korea as Serial No. 1999/45805 on Oct. 21, 1999. The disclosures of the prior applications are considered part of (and are incorporated by reference in) the disclosure of this application. 
    
    
     BACKGROUND OF THE INVENTION 
     1. Field of the Invention 
     The present invention relates to an interpolation method adapted to enlargement of a low resolution image when the image digitized through a CCD (Charged-Coupled Device) has the low resolution, in particular to a filtering control method for improving the image quality of a bi-linear interpolated image which is capable of restoring a requested interpolated high resolution image from a low resolution image by finding a coefficient of a two-dimensional filter on the basis of a regularization image restoration method. 
     2. Description of the Prior Art 
     In the conventional technology, a still picture or a moving picture has or transmits a low resolution image because it can not physically satisfy a sensor having the low resolution or a nyquist value. 
     In addition, a compressed moving picture has or transmits the low resolution image due to its bit value problem. 
     For example, when the compressed moving picture having the low bit value is transmitted to a receiver and the receiver enlarges the transmitted moving picture, the resolution of the transmitted moving picture lowers due to a degradation phenomenon ect. 
     Accordingly, a method for getting a high resolution image from a low resolution image is required. 
     In the meantime, the method for getting the high resolution image from the low resolution image is largely divided into an image expansion type method and an image enhancement type method. 
     First, the image expansion type method converts the size of the low resolution image into a requested size. The bi-linear interpolation method, a zero order expansion method, and a cubic spline method are comprised in the image expansion type method. 
     However, as described above, the image expansion type method has an image visibility lowering problem because when the image is interpolation-restored by the above-mentioned method such as the bi-linear interpolation method, zero order hold expansion method, cubic spline method, the outlines of the image is over-blurred. 
     Meanwhile, the image enhancement type method comprises many methods, but the image enhancement type method causes a computational complexity, accordingly the method is not suited to a real-time processing due to the its computational complexity. 
     In addition, when the image enhancement type method is used for getting the high resolution image from the low resolution image, setting of each parameter is not adaptable. 
     For example, there is a POCS (Projection Onto Convex Set) method for increasing the resolution of an image. In the POCS method, in use of time region information, it is assumed as correlation between the images is uniformly same, but actually the correlation between the images is not uniform. 
     In addition, there is a mapping method for mapping a non-uniform sample of the low resolution image into a uniform sample of the high resolution image by using moving information and segmentation information of the image. However, the mapping method has the computational complexity problem, accordingly the mapping method is not suited to the real-time image data processing of the image processing system. 
     SUMMARY OF THE INVENTION 
     The object of the present invention is to provide a filtering control method for improving the image quality of a bi-linear interpolated image which is capable of improving the image quality of the interpolated image by using an interpolation method considering a real-time processing, a computational complexity and an efficiency when the digital video system seeks the interpolated image from the low resolution image. 
     The other object of the present invention is to provide the filtering control method for improving the image quality of the bi-linear interpolated image which is capable of finding a two-dimensional filter coefficient for getting the interpolated image from the low resolution image on the basis of a regularization image restoration method. 
     The other object of the present invention is to provide the filtering control method for improving the image quality of the bi-linear interpolated image which can approximate and find a PSF (Point Spread Function) for the bi-linear interpolated image from a modeling of the degraded image in the frequency region. 
     The other object of the present invention is to provide the filtering control method for improving the image quality of the bi-linear interpolated image which is capable of performing a real-time adaptive processing by finding a filter coefficient from the bi-linear interpolated image and approximated PSF. 
     In the present invention, in order to find a filter coefficient for finding the interpolated image from the low resolution image on the basis of the regularization image restoration method, when H is the PSF (Point Spread Function), f is a requested high resolution image, Z is the low resolution image, g is the high resolution image gotten from the bi-linear interpolation method, an added function M (f)=∥g−Hf∥ 2 +α∥Cf ∥ 2  for finding the PSF(H) from an equation g=Bz=Hf+n (B, H are bi-linear interpolated filters, n is a noise component generated by the assumed H) is defined. 
     The filtering control method for improving the image quality of the bi-linear interpolated image can be implemented by finding the PSF(H) from the added function M(f) by using an equation 
     
       
         
           
             
               H 
               ⁡ 
               
                 ( 
                 
                   k 
                   , 
                   l 
                 
                 ) 
               
             
             = 
             
               
                 
                   G 
                   ⁡ 
                   
                     ( 
                     
                       k 
                       , 
                       l 
                     
                     ) 
                   
                 
                 
                   F 
                   ⁡ 
                   
                     ( 
                     
                       k 
                       , 
                       l 
                     
                     ) 
                   
                 
               
               . 
             
           
         
       
     
     The filtering control method for improving the image quality of the bi-linear interpolated image can be implemented by finding a PSF(P) of a f=Pg function by using an equation 
               P   ⁡     (     k   ,   l     )       =       H   *     (     k   ,   l     )           H   *     (     k   ,   l     )     ⁢     H   ⁡     (     k   ,   l     )         +     C   *     (     k   ,   l     )     ⁢     C   ⁡     (     k   ,   l     )                   
after finding the PSF(H).
 
     The filtering control method for improving the image quality of the bi-linear interpolated image can restore the requested high resolution image(f) by finding an added filter coefficient Q of the PSF(P) and interpolation filter B from the equation f=Pg=PBz=Qz. 
    
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         FIG. 1  illustrates an image sample for getting a twice enlarged high resolution image according to the embodiment of the present invention. 
         FIG. 2  illustrates an interpolation filter coefficient for getting the twice enlarged image according to the embodiment of the present invention. 
         FIG. 3  illustrates an image sample for getting a three times enlarged high resolution image according to the other embodiment of the present invention. 
         FIG. 4  illustrates the interpolation filter coefficient for getting the three times enlarged image according to the other embodiment of the present invention. 
         FIG. 5  illustrates an image sample for getting a six times enlarged high resolution image according to the another embodiment of the present invention. 
         FIG. 6  illustrates the interpolation filter coefficient for getting the six times enlarged image according to the another embodiment of the present invention. 
     
    
    
     DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
       FIG. 1  illustrates an image sample for getting a twice enlarged high resolution image according to the embodiment of the present invention. 
     As depicted in  FIG. 1 , a˜i illustrate low resolution pixels, A˜D illustrate high resolution pixels. In addition, pixels depicted as ‘x’ illustrate pixels interpolated as twice by a twice interpolation filter coefficient. 
       FIG. 2  illustrates the interpolation filter coefficient for getting a twice enlarged image according to the embodiment of the present invention. In other words, the interpolation filter coefficient for interpolating the twice enlarged image of  FIG. 1  is depicted in  FIG. 2 . 
     As depicted in  FIG. 2 , the high resolution image is gotten from the low resolution pixels a˜i (3×3 pixels) inside of a circle of  FIG. 1  by using the interpolation filter coefficient. 
       FIG. 3  illustrates an image sample for getting a three times enlarged high resolution image according to the other embodiment of the present invention. 
     As depicted in  FIG. 3 , a˜p illustrate the low resolution pixels, A˜I illustrate the high resolution pixels using the filter according to the present invention. 
       FIG. 4  illustrates the interpolation filter coefficient for getting the three times enlarged image according to the other embodiment of the present invention. 
     As depicted in  FIG. 4 , three times enlarged pixels which are newly generated illustrated as triangles in  FIG. 3  are gotten from the low resolution pixels a˜p (4×4 pixels) by using the interpolation filter coefficient of  FIG. 4 . 
       FIG. 5  illustrates the image sample for getting a six times enlarged high resolution image according to the another embodiment of the present invention. In other words, it illustrates the image sample for getting the six times enlarged high resolution image from the twice and three times interpolation filter coefficients by using the bi-linear interpolation method. 
     As depicted in  FIG. 5 , pixels illustrated as a ‘X’ can be gotten by using the twice interpolation filter of  FIG. 2 , and pixels illustrated as a triangle can be gotten by using the three times interpolation filter coefficient of  FIG. 4 . 
     In addition, pixels illustrated as a quadrilateral can be gotten from the pixels generated by the twice and three times interpolation filter coefficients by using the bi-linear interpolation method. 
       FIG. 6  illustrates the interpolation filter coefficient for getting the six times enlarged image according to the another embodiment of the present invention. In other words, the interpolation filter coefficient for getting the six times enlarged image of  FIG. 5  is depicted in  FIG. 6 . 
     Meanwhile, as depicted in  FIG. 2 ,  FIG. 4  and  FIG. 6 , the value found by using the interpolation filter coefficient of the present invention has an integer value. 
     In addition, a 9 bit shift is performed to the value calculated by the interpolation filter coefficient, accordingly there is no need to perform a floating point operation processing. 
     The twice, three times, six times interpolated images are depicted in  FIG. 1˜FIG .  6 , however the present invention is not limited by that, it can be adapted freely to a certain interpolation value. 
     Hereinafter, the filtering control method for improving the image quality of the bi-linear interpolated image will be described in more detail. 
     First, a spatially invariant PSF (Point Spread Function) for finding the interpolation filter coefficient according to the each interpolation value can be easily analyzed and approximated in the frequency region, accordingly the spatially invariant PSF (Point Spread Function) is considered from the bi-linear interpolated image. 
     After that, when it is assumed as the low resolution image is z, high resolution image gotten by the bi-linear interpolation method is g, high resolution image to be restored is f, the relation between the each image can be described as below.
 
g=Bz=Hf+n  [Equation 1]
 
     Herein, the B, H, n are the bi-linear interpolation filters, H is the spatially invariant PSF defining the relation between the original high resolution image and high resolution image gotten by the interpolation method, and the n is a noise component generated by the assumed H. 
     Herein, when the noise component is neglected and a direct inverse is used in order to find the PSF(H), the PSF(H) can be described as below equation 2 in the frequency region. 
     
       
         
           
             
               
                 
                   
                     H 
                     ⁡ 
                     
                       ( 
                       
                         k 
                         , 
                         l 
                       
                       ) 
                     
                   
                   = 
                   
                     
                       
                         G 
                         ⁡ 
                         
                           ( 
                           
                             k 
                             , 
                             l 
                           
                           ) 
                         
                       
                       
                         F 
                         ⁡ 
                         
                           ( 
                           
                             k 
                             , 
                             l 
                           
                           ) 
                         
                       
                     
                     . 
                   
                 
               
               
                 
                   [ 
                   
                     Equation 
                     ⁢ 
                     
                         
                     
                     ⁢ 
                     2 
                   
                   ] 
                 
               
             
           
         
       
     
     Herein, the H(k,l) is the component in the k,l frequency region of the PSF(H), the G (k,l) is the component in the k,l frequency region of the bi-linear interpolated image. In addition, the F (k,l) is the component in the k,l frequency region of the high resolution image. 
     Meanwhile, the high resolution image f to be restored is unknown, the PSF(H) can be gotten from the bi-linear interpolated high resolution image through a statistical processing after performing an under-sample processing of various images as various value. 
     Herein, the high resolution image is gotten by using the PSF(H) found from the direct inverse. In other words, there is a system stabilization problem because the high resolution image gotten from the PSF(H) by using the direct inverse is overshoot in the region where the k,l have ‘0’ value (in general, high frequency region) in the frequency region, accordingly the regularization image restoration for improving the system stabilization is used to solve the problem. 
     The regularization image restoration method is used for restoring the image or finding a certain PSF, an added function M(f) for finding the PSF(H) by using the regularization image restoration method can be described as below equation 3.
 
M(f)=∥g−Hf∥ 2 +α∥Cf∥ 2   [Equation 3]
 
     Herein, the first term of the right side of Equation 3 illustrates the credibility of the bi-linear interpolated image, the second term of the right side illustrates increase of the stability of the system by providing the mitigation to the restored image. 
     In addition, the ∥.∥ means a norm, the α is a regularization parameter for determining the credibility and mitigation of the original image. In addition, the C is the two-dimensional high frequency filter for determining the mitigation of the original image, in the present invention a two-dimensional Gaussian filter is used as the C. 
     When a gradient operator is adapted to Equation 3 in order to get the high resolution image, it can be described as below equation 4.
 
□ f M(f)=−2H T (g−Hf)+2αC T Cf=0  [Equation 4]
 
     Herein, the T means a transpose of a matrix. 
     Meanwhile, conventionally a repetition method is used in order to get the high resolution image and regularization parameter, but it is not suited to the moving picture processing because the method causes lots of computational complexity. 
     Accordingly, in the present invention, the regularization parameter α is fixed as ‘1’, and the high resolution image f can be found as below equation 5. 
     
       
         
           
             
               
                 
                   f 
                   = 
                   
                     
                       
                         
                           H 
                           T 
                         
                         ⁢ 
                         g 
                       
                       
                         ( 
                         
                           
                             
                               H 
                               T 
                             
                             ⁢ 
                             H 
                           
                           + 
                           
                             
                               C 
                               T 
                             
                             ⁢ 
                             C 
                           
                         
                         ) 
                       
                     
                     = 
                     Pg 
                   
                 
               
               
                 
                   [ 
                   
                     Equation 
                     ⁢ 
                     
                         
                     
                     ⁢ 
                     5 
                   
                   ] 
                 
               
             
           
         
       
     
     When the PSF(P) is found by Equation 5, PSF(P)=H/(H T H+C T C) requires the lots of computational complexity for calculating an inverse matrix, however the PSF(P) in Equation 5 is a block-circulant matrix, accordingly it can be easily calculated in the frequency region. 
     Accordingly, the PSF(P) can be found finally as below Equation 6. 
     
       
         
           
             
               
                 
                   
                     P 
                     ⁡ 
                     
                       ( 
                       
                         k 
                         , 
                         l 
                       
                       ) 
                     
                   
                   = 
                   
                     
                       H 
                       * 
                       
                         ( 
                         
                           k 
                           , 
                           l 
                         
                         ) 
                       
                     
                     
                       
                         H 
                         * 
                         
                           ( 
                           
                             k 
                             , 
                             l 
                           
                           ) 
                         
                         ⁢ 
                         
                           H 
                           ⁡ 
                           
                             ( 
                             
                               k 
                               , 
                               l 
                             
                             ) 
                           
                         
                       
                       + 
                       
                         C 
                         * 
                         
                           ( 
                           
                             k 
                             , 
                             l 
                           
                           ) 
                         
                         ⁢ 
                         
                           C 
                           ⁡ 
                           
                             ( 
                             
                               k 
                               , 
                               l 
                             
                             ) 
                           
                         
                       
                     
                   
                 
               
               
                 
                   [ 
                   
                     Equation 
                     ⁢ 
                     
                         
                     
                     ⁢ 
                     6 
                   
                   ] 
                 
               
             
           
         
       
     
     Herein, the ‘*’ the means a complex-conjugate. 
     The PSF(P) can be found by using an IFT (Inverse Fourier Transform) from Equation 6. 
     The requested high resolution image f can be found as below Equation 7 by using the found PSF(P) and Equation 1.
 
f=Pg=PBz=Qz  [Equation 7]
 
     The PSF(P) is the spatially invariant function, the bi-linear interpolation filter B can be easily found by the conventional technology, accordingly the added filter coefficient Q of the PSF(P) and bi-linear interpolation filter B can be found. 
     Herein, in order to reduce the computational complexity, the number of a kernel of the PSF(P) is set in accordance with the up-sampling value. 
     When the up-sampling value is 2 in the present invention, the number of the kernel is limited as 3, when the up-sampling value is 3, the number of the kernel is limited as 4. 
     When the up-sampling value is 2, it can be used in an application segment for enlarging the size of the image as twice at a post processor of the compressed digital image and in finding of a sub-pixel moving vector in a H.263 moving picture compressed method. 
     In addition, when the up-sampling value is 3, it can be used in using of a ⅓ unit moving vector in a H.26L moving picture compressed method. 
     Herein, the H.263 and H.26L are moving picture compressed standards presented in the ITU-T (International Telecommunications Union-Telecommunication). 
     As described above, the present invention can be used for improving the image quality at the post processor of the compressed digital image by using the interpolation method for getting the interpolated high resolution image from the low resolution image when the resolution of the digital image lowers due to the low resolution image sensor. 
     In addition, the interpolation method of the present invention can improve the image quality by finding the moving vector of the moving picture compressed type.