USRE42716E1 - Method for filtering an image - Google Patents

Abstract

The present invention relates to a method for recovering a compressed image for an image processing technique and an apparatus therefor. In the present invention, a cost function is defined in consideration with a directional characteristic of the pixels which will be recovered and a plurality of pixels of the recovering pixels. In addition, a regularization parameter variable having a certain weight is obtained from the cost function, and the regularization parameter variable is approximated using the compressed pixel for thereby obtaining a recovering pixel. The regularization parameter variable has a weight of a reliability with respect to the original image and a weight of a smoothing degree of the original image.In the method, a pixel of an image is filtered using a filtering methodology that adjusts a degree of filtering based on a difference value. The difference value may be based on the pixel being filtered and a neighboring pixel. As an example, the neighboring pixel may be a pixel adjacent to the pixel being filtered.

Description

DIVISIONAL REISSUE APPLICATIONS

Notice: More than one reissue application has been filed for the reissue of U.S. Pat. No. 6,535,642. The reissue applications are application Ser. Nos. 11/081,073; 11/081,075 (the subject application); and 11/338,905.

BACKGROUND OF THE INVENTION

1. Field of the Invention

The present invention relates to an image processing technique, and in particular to a method for recovering a compressed video signal and an apparatus therefor.

2. Description of the Prior Art

The image compression technique of MPEG, MPEG2, H261, H263, etc. is implemented by a Hybrid MC DCT (Motion Compensation Discrete Cosine Transform) technique. This hybrid MC DCT is classified into an encoding process and decoding processes. In the encoding process, the original image is divided into a plurality of blocks for compressing the information of a spacious region, and a two-dimensional DCT is performed with respect to each block, and a redundancy is decreased in the image or between the images using a correlation on a time axis between the images for decreasing the information of the time region. In addition, in the decoding process, the reverse sequence of the decoding process is performed. In order to implement the MCDCT technique, an encoder and decoder are required.

FIG. 1 is a block diagram illustrating a conventional image encoder. As shown therein, an input video signal is subtracted by a subtractor 1 with a motion compensated video signal from a video memory 9 and is inputted via a first switching unit 2 and a DCT unit 3. The DCT unit 3 processes the inputted video signal based on a DCT, and a quantization unit 4 quantizes a DCT-processed video signal and outputs a DCT coefficient q. This coefficient is reversely quantized by a reverse quantizing unit 6 and is processed based on a reverse DCT by a reverse DCT unit 7 for thereby recovering the original video signal. The thusly recovered video signal is summed by a summing unit 8 with a video signal recovered in the earlier process via a second switching unit 10 and is inputted into the video memory. A controller 5 controls the first and second switching units 2 and 10 and transmits an intra/inter information (p=mtype; flag for INTRA/INTER), a transmission information (q; flag for transmitted or not), and an quantizing information (qz=Qp; quantizer indication) to a decoder (not shown in FIG. 1). The video memory 9 outputs a motion vector information (v=MV; motion vector) to the decoder. The DCT unit 3 outputs a DCT coefficient q to the decoder.

While the video signal is being coded, the information may be lost during the quantizing process. Therefore, the video signals reconstructed by the decoder may cause blocking artifacts and ring effects. The block artifacts occur when quantizing a low frequency DCT coefficient, and the ring effects occur due to the information loss of the original video in the quantizing process for a high frequency DCT coefficient.

Namely, in the case of the coding technique using a DCT in a system which is capable of coding a still picture or a motion picture, the entire image is divided into a plurality of small images (for example, 8×8 blocks), and then a transforming operation is performed with respect to the divided blocks, and the original image is processed based on a DCT, and an important information of the original image based on a result of the conversion is included in the low frequency component. As the component becomes high frequency, the important information is decreased. The low frequency component includes an information related to the neighboring block. The DCT transform is performed without considering a correlation between the blocks. Namely, the low frequency components are quantized by the blocks, so that a continuity between the neighboring blocks is lost. This phenomenon is called as the referred to as blocking artifacts.

In addition, when quantizing the coefficients obtained when performing from the DCT operation, as the interval of the quantizing operation is increased, the components to be coded is decreased. Therefore, the number of bits which will be processed is decreased, so that a distortion occurs in the reconstructed original image. This phenomenon is called as the referred to as ring effects. The ring effects, which occur when increasing the intervals of the quantizing operations, are increased at a contour line of an object among in the images.

As a technique for removing the above-described block artifacts and ring effects, a low pass filter technique and a regularization recovering technique are generally used.

The low pass filter sets a filter tap or a filter coefficient based on or by selecting (filter mask) a plurality of pixels near a certain pixel and obtaining an average of the pixels. The recovered images are over smoothed in accordance with the kinds of images, and a compression ratio.

In the regularization recovering method, the block artifacts are adaptively processed in accordance with the statistical characteristic of the images. Namely, a non-uniform information is all computed at all direction boundary areas and in the interior of the block. However, since the computed values have a matrix form, it is impossible to implement a real time computation due to a the large amount of computation. In addition, with an exception except for the amount of non-uniformity, since an average is comprehensively adapted based on a result of the computation of the non-uniform information, in the block having a large amount of non-uniformity, the degree of the non-uniformity is decreased. On the contrary, the degree of the non-uniformity may be increased. Therefore, it is hard to say whether it is well adaptive adapted well to the system.

The above-described two techniques have advantages and disadvantages in view of a complexity and performance increase of the system. Namely, the low pass filter technique has requires less computation amount compared to the regularization recovering technique and has, but has a small capacity for adaptively processing the images, so that the information is lost at an edge portion. The regularization recovering method has an excellent performance and, but requires a large amount of computation when computing regularization parameters.

SUMMARY OF THE INVENTION

Accordingly, in the present invention, it is possible to removing a block artifact and ring effect which occur in a decoded video signal.

In addition, it is possible to define a cost function having a directional feature by the unit of pixels during a decoding operation and obtain a regularization parameter based on the cost function.

To achieve the above objects, there is provided a method for recovering a compressed motion picture according to an embodiment of the invention, comprising the steps of defining a cost function having a smoothing degree of an image and a reliability with respect to an original image in consideration of the directional characteristics of the pixels which will be recovered and a plurality of pixels near the recovering pixels, obtaining a regularization parameter variable having a weight value of a reliability with respect to an original image based on the cost function, and approximating the regularization parameter variable using the compressed pixel and obtaining a recovering pixel.

These and other objects of the present application will become more readily apparent from the detailed description given hereinafter. However, it should be understood that the detailed description and specific examples, while indicating preferred embodiments of the invention, are given by way of illustration only, since various changes and modifications within the spirit and scope of the invention will become apparent to those skilled in the art from this detailed description.

The present invention relates to filtering an image. In one embodiment, a pixel of an image is filtered using a filtering methodology that adjusts a degree of filtering based on a difference value. Here, the difference value may be based on the pixel being filtered and a neighboring pixel. For example, the neighboring pixel may be a pixel adjacent to the pixel being filtered.

In one embodiment, the method includes determining the difference value. For example, the difference between the pixel being filtered and a neighboring pixel may be determined as the difference value.

In another embodiment, the pixel is filtered based on a quantization parameter used in processing a portion of an image including the pixel. For example, the portion of an image including the pixel may be a macroblock.

In yet another embodiment, the filtering methodology includes determining at least one boundary value based on a quantization parameter of a portion of the image including the pixel. Again, as an example, the portion of the image including the pixel may be a macroblock.

BRIEF DESCRIPTION OF THE DRAWINGS

The present invention will become more fully understood from the detailed description given hereinbelow and the accompanying drawings which are given by way of illustration only, and thus as not limitative of the present invention and wherein:

FIG. 1 is a block diagram illustrating a conventional video encoder;

FIG. 2 is a block diagram illustrating an apparatus for recovering a compressed motion picture according to an embodiment of the present invention;

FIG. 3 is a view illustrating pixels and a position information of the pixels for explaining the operation of an embodiment of the present invention;

FIG. 4 is a flow chart of a method for recovering a compressed motion picture according to an embodiment of the present invention; and

FIG. 5 is a flow chart of a method for recovering a compressed motion picture according to another embodiment of the present invention.

DETAILED DESCRIPTION OF THE PREFERRED EXAMPLE EMBODIMENTS

FIG. 2 is a block diagram illustrating an apparatus for recovering a compressed motion picture according to an embodiment of the present invention. As shown therein, a hardware decoder 201 receives an intra/inter information (p=mtype), a transmission information (t), a quantizing information (qz=Qp), a DCT coefficient q, and a motion vector information (v=MV; motion vector) from an a hardware encoder as shown in FIG. 1 and decodes the thusly received information. The hardware encoder and decoder 201 are connected by a communication channel or network. A hardware block removing filter 202 receives a video signal (Y,U,V), a quantizing variable (qz=Qp), a macro block type (mtype), and a motion vector (v=MV) from the decoder 201 and performs an image compressing process according to the present invention for thereby outputting a recovered video signal.

FIG. 3 illustrates pixels and the position of the pixels for explaining the operation of the present invention. As shown therein, assuming the original pixels f(i,j) at the center portion as a reference, f(i,j−1) represents a pixel near the left side, and f(i,j+1) represents a pixel near the right side, and f(i−1,j) represents a pixel near the upper side, and f(i+1,j) represents a pixel near the lower side. Here, i, j represent a position information of each pixel.

A first embodiment of the present invention will be explained with reference to the accompanying drawings.

In the first embodiment of the present invention, a cost function having a directional feature by the unit of pixels is defined, and a regularization parameter is obtained based on the cost function. A recoverable pixel is obtained using a value which is actually adapted to the regularization parameter and is processed based on a DCT and a projection. Then a resultant data is processed based on a reverse DCT for thereby recovering an image similar to the original image. The above-described operation will be explained in detail.

Definition of Cost Function

When the original image f is compressed and transmitted, the image g which is reconstructed by the decoder 201 may be expressed as follows.
g=f+n   (1)
where, g, f, and n have a size of MM×1 rearranged in a scanning sequence, and n represents a quantizing difference.

In order to process the original image f by the unit of pixels, the original pixels f(i,j) having a certain position information(i,j) is adapted. The recovered pixel g(i,j) may be expressed using the original pixel(i,j) and a quantizing difference n(i,j) with respect to the original pixel(i,j).
g(i,j)=f(i,j)+n(i,j)   (2)

As seen Equation 2, a smoothing which represents a non-uniformity degree with respect to the original pixel f(i,j) and the neighboring pixels of the original pixel f(i,j) and a cost function including includes a reliability with respect to the original pixel f(i,j) and the compressed pixel g(i,j).

First, in order to consider the directional features of four pixels f(i,j+1), f(i+1,j), f(i,j−1), and f(i−1,j) with respect to the original pixel f(i,j), the cost functions of M_HL(f(i,j)), M_HR(f(i,j)), M_VT(f(i,j)), M_VD(f(i,j)), M_T(f(i,j)) are defined with respect to the neighboring pixels. In order to set a time based region relationship of the original pixel f(i,j), the cost function M_T(f(i,j)) is defined. Next, the cost functions of M_HL(f(i,j)), M_HR(f(i,j)), M_VT(f(i,j)), M_VD(f(i,j)), M_T(f(i,j)) with respect to the neighboring pixels and the cost function M_T(f(i,j)) of the time region are summed, so that it is possible to obtain the cost function M(f(i,j)) with respect to the original pixel f(i,j) may be obtained based on as shown in Equation (3).
M(f(i,j))=M_HL(f(i,j))+M_HR(f(i,j))+M_VT(f(i,j))+M_VD(f(i,j))+M_T(f(i,j))   (3)
where M_HLrepresents a cost function having a relationship between the pixel f(i,j) and the left side neighboring pixel f(i,j−1), M_HR(f(i,j)) represents a cost function having a relationship between the pixel f(i,j) and the right side neighboring pixel f(i,j+1), M_VT(f(i,j)) represents a cost function having a relationship between the pixel f(i,j) and the upper side neighboring pixel f(i−1,j), M_VD(f(i,j)) represents a cost function having a relationship between the pixel f(i,j) and the lower side neighboring pixel f(i+1,j), and M_T(f(i,j)) represents a cost function having a relationship of the time region.

The cost function having a smoothing degree and reliability may be expressed as the following equation 4.
M_HL(f(i,j))=[f(i,j)−f(i,j−1)]²+α_HL[g(i,j)−f(i,j)]²
M_HR(f(i,j))=[f(i,j)−f(i,j+1)]²+α_HR[g(i,j)−f(i,j)]²
M_VT(f(i,j))=[f(i,j)−f(i−1,j)]²+α_VT[g(i,j)−f(i,j)]²
M_VD(f(i,j))=[f(i,j)−f(i+1,j)]²+α_VD[g(i,j)−f(i,j)]²
MT(f(i,j))=[f(i,j)−f_MC(i,j)]²+α_T[g(i,j)−f(i,j)]²   (4)

As seen in Equation 4, the first term of the right side of each cost function represents a smoothing degree with respect to the original pixel and the neighboring pixel, and the second term of the right side represents a reliability with respect to the original pixel and the recovered pixel.

The first term of the right side of the cost function M_HL(f(i,j)) represents a square value of the difference between the original pixel f(i,j) and the left side neighboring pixel f(i,j−1) and represents a uniformity degree, namely, a smoothed degree of the original pixel f(i,j) and the left side neighboring pixel f(i,j−1) based on the error component between the original pixel f(i,j) and the left side neighboring pixel f(i,j−1). In addition, the second term of the right side of the cost function M_HL(f(i,j)) represents a square value of the difference between the original pixel f(i,j) and the compressed pixel g(i,j) and represents a value for comparing whether a certain difference exists between the compressed pixel g(i,j) and the original pixel f(i,j) based on a difference component between the original pixel f(i,j) and the compressed pixel g(i,j) and represents a reliability of the original pixel f(i,j) and the compressed pixel g(i,j).

In addition, the first term of the right side of M_HR(f(i,j)) represents a smoothing degree of the original pixel f(i,j) and the right side neighboring pixel f(i,j+1), and the second term of the right side represents a reliability of the original pixel f(i,j) and the compressed pixel g(i,j). The first term of the right side of the cost function M_VT(f(i,j)) represents a smoothing degree of the original pixel f(i,j) and the upper side neighboring pixel f(i−1,j), and the second term of the right side represents a reliability of the original pixel, and the compressed pixel g(i,j). The first term of the right side of the cost function M_VT(f(i,j)) represents a smoothing degree of the original pixel f(i,j) and the lower side neighboring pixel f(i+1,j), and the second term of the right side represents a reliability of the original pixel f(i,j) and the compressed pixel g(i,j). M_T(f(i,j)) represents a cost function for setting a relationship of the time region.

The values of α_HL, α_HR, α_VT, α_VD α_T of the second term of the right side represents a regularization parameter and a ratio of a smoothing degree and reliability. These values represent a difference component. In addition, these values represent a weight value with respect to the reliability. As these values are increased, the reliability is enhanced. Since the smoothing degree and the reliability are opposite to each other, the ratio of the smoothing degree and reliability is determined when the regularization parameter is determined. Each regularization parameter may be expressed as the following Equation 5.

$\begin{array}{cc}{\alpha }_{\mathrm{HL}}=\frac{{\left[f\left(i,j\right)-f\left(i,j-1\right)\right]}^2}{{\left[g\left(i,j\right)-f\left(i,j\right)\right]}^2},{\alpha }_{\mathrm{HR}}=\frac{{\left[f\left(i,j\right)-f\left(i,j+1\right)\right]}^2}{{\left[g\left(i,j\right)-f\left(i,j\right)\right]}^2}{\alpha }_{\mathrm{VT}}=\frac{{\left[f\left(i,j\right)-f\left(i-1,j\right)\right]}^2}{{\left[g\left(i,j\right)-f\left(i,j\right)\right]}^2},{\alpha }_{\mathrm{VD}}=\frac{{\left[f\left(i,j\right)-f\left(i+1,j\right)\right]}^2}{{\left[g\left(i,j\right)-f\left(i,j\right)\right]}^2}{\alpha }_T=\frac{{\left[f\left(i,j\right)-f_{\mathrm{MC}}\left(i,j\right)\right]}^2}{{\left[g\left(i,j\right)-f\left(i,j\right)\right]}^2}& \left(5\right)\end{array}$

In the above Equation 5, the denominators of the above-equations represents a difference between the original pixel and the compressed pixel, and the numerator represents a difference between the original pixel and the neighboring pixel.

Computation of Recovering Pixels Based on Cost Function

It is needed to obtain the recovering pixels which is the original pixels. However, the cost function includes a square with respect to the original pixel. Therefore, the cost function is partially differentiated with respect to the original pixel, so that it is possible to obtain the original pixels based on the differentiated values. The cost function M(f(i,j)) may be differentiated based on Equation 3.

$\begin{array}{cc}\partial M\frac{\left(f\left(i,j\right)\right)}{\partial f\left(i,j\right)}=\frac{\partial M_{\mathrm{HL}}\left(f\left(i,j\right)\right)}{\partial f\left(i,j\right)}+\frac{\partial M_{\mathrm{HR}}\left(f\left(i,j\right)\right)}{\partial f\left(i,j\right)}+\frac{\partial M_{\mathrm{VT}}\left(f\left(i,j\right)\right)}{\partial f\left(i,j\right)}+\frac{\partial M_{\mathrm{VD}}\left(f\left(i,j\right)\right)}{\partial f\left(i,j\right)}+\frac{\partial M_T\left(f\left(i,j\right)\right)}{\partial f\left(i,j\right)}=0& \left(6\right)\end{array}$

Each term of the right side of the cost function with respect to the neighboring pixels is as follows.

$\begin{array}{cc}\frac{\partial M_{\mathrm{HL}}\left(f\left(i,j\right)\right)}{\partial f\left(i,j\right)}=2\left[f\left(i,j\right)-f\left(i,j-1\right)\right]-2{\alpha }_{\mathrm{HL}}\left[g\left(i,j\right)-f\left(i,j\right)\right]\frac{\partial M_{\mathrm{HR}}\left(f\left(i,j\right)\right)}{\partial f\left(i,j\right)}=2\left[f\left(i,j\right)-f\left(i,j+1\right)\right]-2{\alpha }_{\mathrm{HR}}\left[g\left(i,j\right)-f\left(i,j\right)\right]\frac{\partial M_{\mathrm{VT}}\left(f\left(i,j\right)\right)}{\partial f\left(i,j\right)}=2\left[f\left(i,j\right)-f\left(i-1,j\right)\right]-2{\alpha }_{\mathrm{VT}}\left[g\left(i,j\right)-f\left(i,j\right)\right]\frac{\partial M_{\mathrm{VD}}\left(f\left(i,j\right)\right)}{\partial f\left(i,j\right)}=2\left[f\left(i,j\right)-f\left(i+1,j\right)\right]-2{\alpha }_{\mathrm{VD}}\left[g\left(i,j\right)-f\left(i,j\right)\right]\frac{\partial M_T\left(f\left(i,j\right)\right)}{\partial f\left(i,j\right)}=2\left[f\left(i,j\right)-f_{\mathrm{MC}}\left(i,j\right)\right]-2{\alpha }_T\left[g\left(i,j\right)-f\left(i,j\right)\right]& \left(7\right)\end{array}$

The values of Equation 7 are substituted for Equation 6, and the pixels which will be finally recovered are in the following Equation 8.

$\begin{array}{cc}f\left(i,j\right)=\frac{\begin{array}{c}f\left(i,j-1\right)+f\left(i,j+1\right)+f\left(i-1,j\right)+\\ f\left(i+1,j\right)+f_{\mathrm{MC}}\left(i,j\right)+{\alpha }_{\mathrm{TOT}}g\left(i,j\right)\end{array}}{5+{\alpha }_{\mathrm{TOT}}}{\alpha }_{\mathrm{TOT}}={\alpha }_{\mathrm{HL}}+{\alpha }_{\mathrm{HR}}+{\alpha }_{\mathrm{VT}}+{\alpha }_{\mathrm{VD}}+{\alpha }_T& \left(8\right)\end{array}$

The pixels expressed by Equation 8 are the pixels included in the inter macro block. However, the pixels of the macro block coded into the intra macro type based on Equation 6 is

$\frac{\partial M_T\left(f\left(i,j\right)\right)}{\partial f\left(i,j\right)}=0$
because there is not a motion information on tile time axis. Therefore, the pixels included in the intra macro block may be expressed in the following Equation 9.

$\begin{array}{cc}f\left(i,j\right)=\frac{\begin{array}{c}f\left(i,j-1\right)+f\left(i,j+1\right)+f\left(i-1,j\right)+\\ f\left(i+1,j\right)+{\alpha }_{\mathrm{TOT}}g\left(i,j\right)\end{array}}{4+{\alpha }_{\mathrm{TOT}}}{\alpha }_{\mathrm{TOT}}={\alpha }_{\mathrm{HL}}+{\alpha }_{\mathrm{HR}}+{\alpha }_{\mathrm{VT}}+{\alpha }_{\mathrm{VD}}& \left(9\right)\end{array}$

Therefore, the pixels included in the inter macro block are obtained based on a filter strength of Equation 8, and the pixels included in the intra macro block are obtained based on a filter strength of Equation 9. Whether the pixels of the macro block are coded in the intra macro type or in the inter macro type are determined by the intra inter information (p=mtype).

As seen in Equations 8 and 9, the recovering pixels include a regularization parameter α, and each regularization parameter variable is approximated as follows.

Approximation of Regularization Parameter Variable

As seen in Equation 5, each regularization parameter variable includes an original pixel, a neighboring pixel, and a recovering pixel (compressed pixel). In addition, since the original pixel f(i,j) and four neighboring pixels f(i,j−1), f(i,j+1), f(i−1,j), f(i+1,j) are the original pixels, these values do not exist in the decoder. Therefore, the pixels f(i,j), f(i,j−1), f(i,j+1), f(i−1,j), f(i+1,j) may not be used for an actual computation. Therefore, in order to actually use the pixels f(i,j), f(i,j−1), f(i,j+1), f(i−1,j), f(i+1,j), the compressed pixels g(i,j), g(i,j−1), g(i,j+1), g(i−1,j), g(i+1,j) must be approximated. To implement the above-described approximation, the following three cases are assumed.

First, the quantizing maximum difference of the macro block unit is a quantizing variable (Qp).

Second, a quantizing difference of each DCT coefficient is uniformly allocated to each pixel of a corresponding macro block,

Third, the non-uniform values between two pixels of the original image are statistically similar to the non-uniform values between two pixels of the compressed image.

As seen in the following Equation 10, each regularization variable is approximated based on the above-described three cases.

$\begin{array}{cc}{\alpha }_{\mathrm{HL}}=\frac{{\left[f\left(i,j\right)-f\left(i,j-1\right)\right]}^2}{{\left[g\left(i,j\right)-f\left(i,j\right)\right]}^2}\approx \frac{{\left[g\left(i,j\right)-g\left(i,j-1\right)\right]}^2}{Q_{\mathrm{pl}}^2}{\alpha }_{\mathrm{HR}}=\frac{{\left[f\left(i,j\right)-f\left(i,j+1\right)\right]}^2}{{\left[g\left(i,j\right)-f\left(i,j\right)\right]}^2}\approx \frac{{\left[g\left(i,j\right)-g\left(i,j+1\right)\right]}^2}{Q_{\mathrm{pl}}^2}{\alpha }_{\mathrm{VT}}=\frac{{\left[f\left(i,j\right)-f\left(i-1,j\right)\right]}^2}{{\left[g\left(i,j\right)-f\left(i,j\right)\right]}^2}\approx \frac{{\left[g\left(i,j\right)-g\left(i-1,j\right)\right]}^2}{Q_{\mathrm{pl}}^2}{\alpha }_{\mathrm{VD}}=\frac{{\left[f\left(i,j\right)-f\left(i+1,j\right)\right]}^2}{{\left[g\left(i,j\right)-f\left(i,j\right)\right]}^2}\approx \frac{{\left[g\left(i,j\right)-g\left(i+1,j\right)\right]}^2}{Q_{\mathrm{pl}}^2}{\alpha }_T=\frac{{\left[f\left(i,j\right)-f_{\mathrm{MC}}\left(i,j\right)\right]}^2}{{\left[g\left(i,j\right)-f\left(i,j\right)\right]}^2}\approx \frac{{\left[g\left(i,j\right)-f_{\mathrm{MC}}\left(i,j\right)\right]}^2}{Q_{\mathrm{pl}}^2}& \left(10\right)\end{array}$

• • where 1 represents the l-th macro block, and Qpl represents a quantizing variable of the l-th macro block. As seen in Equation 10, the difference between the original pixel which is the denominator component of each regularization parameter variable and the compressed pixel is approximated based on the quantizing maximum difference, and the difference between the original pixel which is the numerator component and the compressed pixel is approximated based on the difference with respect to the difference value between the compressed pixel and the neighboring pixel.

The thusly approximated regularization parameter variable is substituted for Equation 8 or 9 for thereby obtaining a result value f(i,j).

FIG. 4 is a flow chart illustrating a method for recovering a compressed motion picture according to the present invention.

As shown therein, in Step ST1, whether the processing pixels are referred to the pixels of the intra macro block or the pixels of the inter macro block is judged. As a result of the judgement, in Steps ST2 and ST3, the regularization parameter variable is obtained. Namely, if the processing pixels are referred to the pixels of the intra macro block, in Step ST2, the regularization parameter variables α_HL, α_HR, α_VT, α_VDar are obtained based on Equation 9. In addition, if the processing pixels are referred to the pixels of the inter macro block, the regularization parameter variables α_HL, α_HR, α_VT, α_VD, α_T are obtained in Step ST3. In addition, the pixel f(i,j) is obtained in Step ST4 based on the obtained regularization parameter variable. At this time, if the processing pixels are referred to the pixels of the inter macro block, and the pixels are obtained based on Equation 8, and if the processing pixels are referred to the pixels of the inter macro block, the pixels are obtained based on Equation 9.

Recovering the Images Using a Projection Technique

In Step ST5, a DCT is performed with respect to the pixel f(ij) f(i,j), and then a quantizing process is performed therefor. Here, the DCT coefficient of the pixel f(i,j) may be expressed as F(u,v).

The value G(u,v) which is DCT-processed with respect to the compressed image g(i,j) may be expressed in the DCT region based on the following Equation 11.
G=QBf   (11)

• • where B represents a DCT process, and Q represents a quantizing process.

The DCT coefficient of the original image and the DCT coefficient of the compressed image have the following interrelationship as seen in Equation 12.
G(u,v)−Qpl≦F(u,v)≦G(u,v)+Qpl   (12)

• • where G(u,v) represents a (u,u)-th value of the two-dimensional DCT coefficient of the compressed image, F(u,v) represents a (u,v)-th value of the two-dimensional DCT coefficient of the original image, Qpl represents the quantizing maximum difference of the l-th macro block, and each DCT coefficient value represents a subset for setting the range of the DCT coefficient of the recovered images. Therefore, the recovered images must be projected based on the subset of Equation 12, and this process is performed in Step ST6 as seen in the following Equation 13.
    P(F(u,v))=G(u,v)−Qpl if F(u,v)<G(u,v)−Qpl
    P(F(u,v))=G(u,v)+Qpl if F(u,v)>G(u,v)+Qpl
    P(F(u,v))=F(u,v) otherwise   (13)

The Equation 13 will be explained in detail.

If F(u,v) is smaller than G(u,v)−Qpl, the projected recovering image P(F(u,v) is mapped based on G(u,v)−Qpl, and if F(u,v) is larger than G(u,v)−Qpl, the projected recovering image P(F(u,v)) is mapped based on G(u,v)+Qpl, otherwise P(F(u,v)) is directly mapped based on the projected recovering image F(u,v).

The mapped image P(F(u,v)) is reversely DCT-processed in the spacious region in Step ST7, and the finally recovered image may be expressed by the following Equation 14.
f′=B^TPBf=B^TPBK(g)   (14)

• • where K(g) represents a computation of the recovering pixels of Equation 8 or 9, BK(g) represents a block DCT coefficient, PBK(g) represents a projected block DCT coefficient, and B^TPBK(g) represents that the projected block DCT coefficient is recovered in the spacious region. The recovered image is stored in the image memory and is outputted.

In the present invention, it is possible to eliminate a block artifact and ring effect based on an non-uniform degree and reliability of the recovered image using a plurality of information from the decoder.

Repetition Technique

If the block artifact and ring effect are not fully eliminated from the recovered pixels, he the above-described processes may be repeatedly performed. As the process for eliminating the block artifact and ring effect is repeatedly performed, the block artifact and ring effect of the recovering image is more further eliminated. In this case, a blurring phenomenon occurs in the edge region of the image. Therefore, the number of repetition must be determined based on the block artifact and ring effect and the blurring phenomenon which is opposite thereto.

The recovering image f_k+1(i,j) is as follows based on Equations 15 and 16 by repeating the above-described process by k-times.

$\begin{array}{cc}{f}_{k+1}\left(i,j\right)=\frac{\begin{array}{c}{f}_k\left(i,j-1\right)+f_k\left(i,j+1\right)+f_k\left(i-1,j\right)+\\ f_k\left(i+1,j\right)+f_{\mathrm{MC}}\left(i,j\right)+{\alpha }_{k,\mathrm{TOT}}g\left(i,j\right)\end{array}}{5+{\alpha }_{k,\mathrm{TOT}}}{\alpha }_{k,\mathrm{TOT}}={\alpha }_{k,\mathrm{HL}}+{\alpha }_{k,\mathrm{HR}}+{\alpha }_{k,\mathrm{VT}}+{\alpha }_{k,\mathrm{VD}}+{\alpha }_{k,T}& \left(15\right)\\ f_{k+1}\left(i,j\right)=\frac{\begin{array}{c}{f}_k\left(i,j-1\right)+f_k\left(i,j+1\right)+f_k\left(i-1,j\right)+\\ f_k\left(i+1,j\right)+{\alpha }_{k,\mathrm{TOT}}g\left(i,j\right)\end{array}}{4+{\alpha }_{k,\mathrm{TOT}}}{\alpha }_{k,\mathrm{TOT}}={\alpha }_{k,\mathrm{HL}}+{\alpha }_{k,\mathrm{HR}}+{\alpha }_{k,\mathrm{VT}}+{\alpha }_{k,\mathrm{VD}}& \left(16\right)\end{array}$

The image expressed based on Equation 15 is related to the pixels included in the inter macro block, and the image expressed based on Equation 16 is related to the pixels included in the intra macro block.

The regularization parameter variables are obtained based on Equations (15) 15 and 16, and the DCT is performed with respect thereto, and the projection technique is adapted. Then, and then the reverse DCT is performed therefor, so that the final recovering image is obtained based on Equation 17.
f_k+1′=B^TPB{circumflex over (f)}=B^TPBK(g_k+1)   (17)

Namely, the block artifacts and ring effects are eliminated from the recovered images by an adaptive decoding operation, so that a real time process is implemented in the digital video apparatus. In particular, it is possible to enhance the resolution in the compression images which require a low bit ration ratio or high speed process.

Next, another embodiment of the present invention will be explained. This embodiment of the present invention is basically directed to decreasing the computation amount and time compared to the earlier embodiment of the present invention. The operation thereof is performed by the recovering apparatus, as shown in FIG. 2, of the compression motion picture according to the present invention. First, the cost function may be defined as seen in Equation 18.
M(f(i,j)=M_L(f(i,j))+M_R(f(i,j))+M_U(f(i,j))+M_D(f(i,j))   (18)

• • where M_L represents a cost function having an interrelationship between the pixel f(i,j) and the left side neighboring pixel f(i,j−1), M_R(f(i,j)) represents a cost function having an interrelationship between the pixel f(i,j) and the right side neighboring pixel f(i,j+1), M_U(f(i,j)) represents a cost function having an interrelationship between the pixel f(i,j) and the upper side neighboring pixel f(i−1,j), and M_D(f(i,j)) is a cost function having an interrelationship between the pixel f(i,j) and the lower side neighboring pixel f(i+1,j).

Next, the cost functions including a smoothing degree and reliability are defined. The regularization parameter variable is included in only the portion (the second term of the right side in Equation 4) of the reliability with respect to the original pixel and recovered pixel. Differently from this construction, in another embodiment of the present invention, the regularization parameter variable is included in the portion which represents a reliability of the original pixel and recovered pixel as well as is included in the portion which represents the smoothing degree with respect to the original pixel and the neighboring pixel. In addition, the smoothing degree and the reliability of the pixel are opposite each other inn in their meaning. Each cost function may be expressed based on Equation 19 as follows:. Equation 19.
M_L(f(i,j))=α_L(f(i,j))[f(i,j)−f(i,j−1)]²+(1−α_L(f(i,j)))[g(i,j)−f(i,j)]²
M_R(f(i,j))=α_R(f(i,j))[f(i,j)−f(i,j+1)]²+(1−α_R(f(i,j)))[g(i,j)−f(i,j)]²
M_U(f(i,j))=α_U(f(i,j))[i f(i,j)−f(i−1i,j)]²+(1−α_U(f(i,j)))[g(i,j)−f(i,j)]² M_U(f(i,j))=α_U(f(i,j))[f(i,j)−f(i−1,j)]²+(1−α_U(f(i,j)))[g(i,j)−f(i,j)]²
M_D(f(i,j))=α_D(f(i,j))[f(i,j)−f(i+1,j)]²+(1−α_D(f(i,j)))[g(i,j)−f(i,j)]²  (19)

As seen in Equation 19, the first term of the right side represents a smoothing degree with respect to the original pixel and the neighboring pixel, and the second term of the right side represents a reliability with respect to the original pixel and the recovered pixel. Here, α_L, α_R, α_U, α_D represent a regularization parameter variable with respect to each cost function and represent a ratio of a smoothing degree and reliability as a difference component. For example, α_L represents a weight value with respect to the smoothing degree, and 1−α_L represents a weight value with respect to the reliability. Therefore, as the regularization parameter variable is increased, the smoothing degree is increased, and the reliability is decreased. Since the regularization includes the right side first term and the left side term of the cost function, it is possible to implement more stable smoothing degree and reliability compared to the earlier embodiment of the present invention.

Next, in order to obtain the recovering pixel, the cost function is partially differentiated with respect to the original pixel. The thusly differentiated value is obtained by the following Equation 20.

$\begin{array}{cc}\partial M\frac{\left(f\left(i,j\right)\right)}{\partial f\left(i,j\right)}=\frac{\partial M_L\left(f\left(i,j\right)\right)}{\partial f\left(i,j\right)}+\frac{\partial M_R\left(f\left(i,j\right)\right)}{\partial f\left(i,j\right)}+\frac{\partial M_R\left(f\left(i,j\right)\right)}{\partial f\left(i,j\right)}+\frac{\partial M_U\left(f\left(i,j\right)\right)}{\partial f\left(i,j\right)}+\frac{\partial M_D\left(f\left(i,j\right)\right)}{\partial f\left(i,j\right)}=0& \left(20\right)\end{array}$

The terms of the right side of Equation 20 are as follows:

$\begin{array}{cc}\frac{\partial M_L\left(f\left(i,j\right)\right)}{\partial f\left(i,j\right)}=2{\alpha }_L\left(f\left(i,j\right)\right)\left[f\left(i,j\right)-f\left(i,j-1\right)\right]-2\left(1-{\alpha }_L\left(f\left(i,j\right)\right)\right)\left[g\left(i,j\right)-f\left(i,j\right)\right]\frac{\partial M_R\left(f\left(i,j\right)\right)}{\partial f\left(i,j\right)}=2{\alpha }_R\left(f\left(i,j\right)\right)\left[f\left(i,j\right)-f\left(i,j+1\right)\right]-2\left(1-{\alpha }_R\left(f\left(i,j\right)\right)\right)\left[g\left(i,j\right)-f\left(i,j\right)\right]\frac{\partial M_U\left(f\left(i,j\right)\right)}{\partial f\left(i,j\right)}=2{\alpha }_U\left(f\left(i,j\right)\right)\left[f\left(i,j\right)-f\left(i-1,j\right)\right]-2\left(1-{\alpha }_U\left(f\left(i,j\right)\right)\right)\left[g\left(i,j\right)-f\left(i,j\right)\right]\frac{\partial M_D\left(f\left(i,j\right)\right)}{\partial f\left(i,j\right)}=2{\alpha }_D\left(f\left(i,j\right)\right)\left[f\left(i,j\right)-f\left(i+1,j\right)\right]-2\left(1-{\alpha }_D\left(f\left(i,j\right)\right)\right)\left[g\left(i,j\right)-f\left(i,j\right)\right]& \left(21\right)\end{array}$

When the values expressed based on Equation 21 are substituted for Equation 20, the finally recovered pixels are obtained based on the following Equation 22.

$\begin{array}{cc}f\left(i,j\right)=\frac{\begin{array}{c}{\alpha }_{L}f\left(i,j-1\right)+{\alpha }_{R}f\left(i,j+1\right)+{\alpha }_{U}f\left(i-1,j\right)+\\ {\alpha }_{D}f\left(i+1,j\right)+\left(4-{\alpha }_{\mathrm{TOT}}\right)g\left(i,j\right)\end{array}}{4}{\alpha }_{\mathrm{TOT}}={\alpha }_L+{\alpha }_R+{\alpha }_U+{\alpha }_D& \left(22\right)\end{array}$

In addition, in the macro block type (mtype), the bit value which is defined as COD is included. This COD includes an information of the macro block. If COD value is ‘0’, it means the coded macro block, and if COD value is ‘1’, it means the non-coded macro block (not coded). Namely, it is possible to Recognize recognize whether the pixels of the current macro block are the same as the pixels of the previously transmitted macro block. If COD value is ‘0’, it means that the macro block of the previous compressed image is different from the macro block of the current image, and if COD value is ‘1’, it means that the macro block of the previous image is the same as the macro block of the current image. Therefore, if COD value is ‘0’, the value is recovered based on Equation 22, and if COD value is ‘1’, as seen in Equation 23, the recovered pixel value fp(i,j) is substituted for the current pixel value with respect to the macro block of the previous image.
f(i,j)=fp(i,j)   (23)

Next, as seen in Equation 22, the recovering pixel includes a regularization parameter variable α, and each regularization parameter variable is obtained as follows.

The regularization parameter variable is obtained based on Equation 19. Namely, since the smoothing degree and reliability are opposite to each other, the regularization parameter variable may be arranged according to Equation 24 as follows based on a ratio of the smoothing degree and the reliability. Equation 24 may be expressed as follows.

$\begin{array}{cc}\frac{1-{\alpha }_L}{{\alpha }_L}=\frac{{\left[f\left(i,j\right)-f\left(i,j-1\right)\right]}^2}{{\left[f\left(i,j\right)-g\left(i,j\right)\right]}^2}\frac{1-{\alpha }_R}{{\alpha }_R}=\frac{{\left[f\left(i,j\right)-f\left(i,j+1\right)\right]}^2}{{\left[f\left(i,j\right)-g\left(i,j\right)\right]}^2}\frac{1-{\alpha }_U}{{\alpha }_U}=\frac{{\left[f\left(i,j\right)-f\left(i-1,j\right)\right]}^2}{{\left[f\left(i,j\right)-g\left(i,j\right)\right]}^2}\frac{1-{\alpha }_D}{{\alpha }_D}=\frac{{\left[f\left(i,j\right)-f\left(i+1,j\right)\right]}^2}{{\left[f\left(i,j\right)-g\left(i,j\right)\right]}^2}& \left(24\right)\end{array}$

In order to obtain the regularization parameter variable expressed as Equation 24, the pixels f(i,j), f(ij−1), f(i,j+1), f(i−1,j), f(i+1,j) must be approximated based on the compressed pixels g(i,j), g(i,j−1), g(i,j+1), g(i−1,j), g(i+1,j) which may be actually used. For implementing the above-described operation, the following three cases are assumed.

First, a quantization difference of each pixel is a function of a quantization variable Qp which is set by the unit of macro blocks.

Second, since the block artifacts generating at a block boundary has a certain non-uniformity degree which is larger than the ring effect occurring in the interior of the block, the difference with respect to the pixels positioned at the block boundary is more largely reflected compared to the pixels positioned in the interior of the block. Namely, a weight value is provided to the difference based on the position of the pixels.

Equation 24 is approximated to Equation 25 based on the above-described two assumptions.

$\begin{array}{cc}\frac{1-{\alpha }_L}{{\alpha }_L}=\frac{{\left[f\left(i,j\right)-f\left(i,j-1\right)\right]}^2}{{\left[f\left(i,j\right)-g\left(i,j\right)\right]}^2}\approx \frac{{\left[g\left(i,j\right)-g\left(i,j-1\right)\right]}^2}{\Phi \left(Q_p\right)}\frac{1-{\alpha }_R}{{\alpha }_R}=\frac{{\left[f\left(i,j\right)-f\left(i,j+1\right)\right]}^2}{{\left[f\left(i,j\right)-g\left(i,j\right)\right]}^2}\approx \frac{{\left[g\left(i,j\right)-g\left(i,j+1\right)\right]}^2}{\Phi \left(Q_p\right)}\frac{1-{\alpha }_U}{{\alpha }_U}=\frac{{\left[f\left(i,j\right)-f\left(i-1,j\right)\right]}^2}{{\left[f\left(i,j\right)-g\left(i,j\right)\right]}^2}\approx \frac{{\left[g\left(i,j\right)-g\left(i-1,j\right)\right]}^2}{\Phi \left(Q_p\right)}\frac{1-{\alpha }_D}{{\alpha }_D}=\frac{{\left[f\left(i,j\right)-f\left(i+1,j\right)\right]}^2}{{\left[f\left(i,j\right)-g\left(i,j\right)\right]}^2}\approx \frac{{\left[g\left(i,j\right)-g\left(i+1,j\right)\right]}^2}{\Phi \left(Q_p\right)}& \left(25\right)\end{array}$

• • where Φ(Qp) is a function of the quantizing variable Qp and is different based on the position of a pixel.
    Therefore, with consideration of the position of each pixel in the function Φ(Qp), Φ(Qp) may be expressed as K_LQP² with respect to α_L, and Φ(Qp) is expressed as K_RQp² with respect to α_R, and Φ(Qp) is expressed as K_UQp², with respect to α_U, and Φ(Qp) is expressed as K_DQp² with respect to α_D. Here, constants K_L, K_R, K_U, K_D are weight values and are different based on whether the neighboring pixel is positioned at the block boundary or in the interior of the block. With consideration to the position of each pixel, type regularization parameter variable is approximated based on the following Equation 26.

$\begin{array}{cc}{\alpha }_L=\frac{K_L{Q}_P^2}{{\left[g\left(i,j\right)-g\left(i,j-1\right)\right]}^2+K_L{\mathrm{<figure-callout\; id="2"\; label="Q\; P"\; filenames="USRE042716-20110920-D00001.png"\; state="\{\{state\}\}">Q</figure-callout>}}_P^{}}{\alpha }_R=\frac{K_R{Q}_P^2}{{\left[g\left(i,j\right)-g\left(i,j+1\right)\right]}^2+K_R{\mathrm{<figure-callout\; id="2"\; label="Q\; P"\; filenames="USRE042716-20110920-D00001.png"\; state="\{\{state\}\}">Q</figure-callout>}}_P^{}}{\alpha }_U=\frac{K_U{Q}_P^2}{{\left[g\left(i,j\right)-g\left(i-1,j\right)\right]}^2+K_U{\mathrm{<figure-callout\; id="2"\; label="Q\; P"\; filenames="USRE042716-20110920-D00001.png"\; state="\{\{state\}\}">Q</figure-callout>}}_P^{}}{\alpha }_D=\frac{K_D{Q}_P^2}{{\left[g\left(i,j\right)-g\left(i+1,j\right)\right]}^2+K_D{Q}_P^2}& \left(26\right)\end{array}$

Assuming that one block is formed of 8×8 number of pixels, namely, assuming that I and j of f(i,j) is 8, respectively, the weight values K_L, K_R, K_U, K_D may be expressed as follows.

• • K_L={9, if j mod 8=0; 1, otherwise}
  • K_R={9, if j mod 8=7; 1, otherwise}
  • K_U={9, if i mod 8=0; 1, otherwise}
  • K_D={9, if i mod 8=7; 1, otherwise}

For example, in the Equation related to K_L, if the residual is 0 when dividing j by 8, K_L is 9, and otherwise, K_L is 1.

When the approximated regularization parameter values are substituted for Equation 22, it is possible to obtain a resultant value f(i,j).

FIG. 5 is a flow chart illustrating a method for recovering a compressed image for an image processing system according to another embodiment of the present invention.

In Step ST10, it is judged whether the pixels of the current macro block are the same as the pixels of the previously transmitted macro block based on the COD value. If they are same, in Step ST11, the recovering pixel values are substituted for the pixel values which are previously recovered based on Equation 23. If they are not the same, in Step ST12, the regularization parameter variables α_L, α_R, α_U, α_D are obtained based on Equation 26, and the recovering pixel f(i,j) is obtained based on Equation 22 in Step ST13.

As described above, in the present invention, a certain weight is provided to the regularization parameter variable, which will be approximated, based on the position of the pixels in consideration with the reliability and smoothing degree as well as the regularization parameter variables, so that it is possible to obtain a value which is near the actual pixel value. Therefore, in the present invention, it is not needed to perform a projection method and a repetition method. In addition, in the present invention, the computation amount and time are significantly decreased.

The invention being thus described, it will be obvious that the same may be varied in many ways. Such variations are not to be regarded as a departure from the spirit and scope of the invention, and all such modifications as would be obvious to one skilled in the art are intended to be included within the scope of the following claims invention.

Claims (44)

1. A method for recovering a compressed motion picture, comprising the steps of:

defining a cost function having a smoothing degree of an image and a reliability with respect to an original image in consideration of the directional characteristics of the pixels which will be recovered and a plurality of pixels near the pixels which will be recovered;

obtaining a regularization parameter variable having a weight value of the reliability with respect to the original image based on a cost function; and

approximating the regularization parameter variable using the compressed pixel and obtaining a pixel which will be recovered,

wherein said regularization parameter variable is a weight value with respect to reliability and is determined based on a difference between the original pixel and the compressed pixel and a difference value between the original pixel and the neighboring pixel.

2. The method of claim 1, wherein said cost function includes another cost function for setting an interrelationship of a time region with respect to the recovering pixel when the pixel which will be recovered is in an inter macro block.

3. The method of claim 1, wherein said cost function includes another cost function which is defined based on a smoothing degree which is obtained by computing a difference between the recovering pixel and the neighboring pixel, a reliability of the original image obtained by computing a difference between the original image and the compressed image, and an interrelationship of a time region of the pixels of the block having a motion information.

4. The method of claim 1, wherein said plurality of neighboring pixels are the pixels which are neighboring in the upper, lower, left and right side directions of the recovering pixels.

5. The method of claim 1, wherein said difference value between the original pixel and the compressed pixel is approximated based on a quantizing maximum difference, and a difference value between the original pixel and the neighboring pixel is approximated based on a difference value between the compressed pixel and the neighboring compressed pixel.

6. The method of claim 1, after the step for obtaining the recovering pixel, further comprising a step for performing a DCT with respect to the recovering pixels, projecting the recovering pixels in accordance with pixel value which will be processed, and performing a reverse DCT with respect to the projected images, and in said projecting step, a recovering image is projected at a subset for setting a range of DCT coefficients of a compressed image, and a maximum quantizing difference of the macro block is included in the subset.

7. The method of claim 1, wherein in said step for approximating the regularization parameter variable, a quantizing maximum difference of a macro block unit is a quantizing variable, a quantizing difference is uniformly allocated to each pixel in a corresponding macro block, and the non-uniform values between two pixels of the original image are statistically similar to the non-uniform values between two pixels of the compressed image.

8. The method claim 1, wherein said regularization parameter variable includes a weight value of a smoothing degree of the original image based on the cost function.

9. The method of claim 8, wherein when the pixels of the current macro block are the same as the pixels of the previously transmitted macro block, the recovered pixel values are substituted for the current pixel values with respect to the macro block of the previous image.

10. The method of claim 8, wherein in said step for approximating the regularization parameter variable, a quantizing difference of each pixel is set based on a function of a quantizing variable set by the unit of a macro block, and a weight value is added to the pixel based on the pixel position.

11. In a method for recovering a compressed motion image for processing an original pixel f(i,j) based on a DCT by the unit of macro blocks of a M×M size, quantizing the DCT-processed coefficient, transmitting together with motion vector information, reversely quantizing and reversely DCT-processing the compressed pixel g(i,j) and recovering an image similar to the original image, a method for recovering a compressed motion picture, comprising the steps of:

defining a cost function M(i,j) having a smoothing degree of an image and a reliability with respect to an original image as a pixel unit in consideration of a directional characteristic between the pixels which will be recovered and the pixels neighboring the pixels which will be recovered;

adaptively searching a regularization parameter variable having a weight of a reliability with respect to the original image from the cost function M(i,j); and

obtaining a projected pixel P(F(u,v)) using a projection method for mapping the pixels which will be recovered in accordance with a range value of the pixels which will be recovered,

wherein said regularization parameter variable is a weight value with respect to reliability and is determined based on a difference between the original pixel and the compressed pixel and a difference value between the original pixel and the neighboring pixel.

12. The method of claim 11, wherein said cost function M(i,j) is formed of a cost function M_HL(f(i,j)) which represents a smoothing degree and a reliability with respect to an original pixel f(i,j) and a left side neighboring pixel f(i,j−1), a cost function M_HR(f(i,j)) which represents a smoothing degree and a reliability with respect to the original pixel f(ij)and a right side neighboring pixel f(i,j+1), a cost function M_VT(f(i,j)) which represents a smoothing degree and a reliability with reflect to the original pixel f(i,j) and an upper side neighboring pixel f(i−1,j), a cost function M_VD(f(i,j)) which represents a smoothing degree and a reliability with respect to the original pixel f(i,j) and a lower side neighboring pixel f(i+1,j), and a cost function M_T(f(i,j)) for setting an interrelationship of a time region with respect to the original pixel.

13. The method of claim 12, wherein each cost function is obtained according to the following equations:

M_HL(f(i,j))=[f(i,j)−f(i,j−1)]²+α_HL[g(i,j)−f(i,j)]²

M_HR(f(i,j))=[f(i,j)−f(i,j+1)]²+α_HR[g(i,j)−f(i,j)]²

M_VT(f(i,j))=[f(i,j)−f(i−1,j)]²+α_VT[g(i,j)−f(i,j)]²

M_VD(f(i,j))=[f(i,j)−f(i+1,j)]²+α_VD[g(i,j)−f(i,j)]²

M_T(f(i,j))=[f(i,j)−f_MC(i,j)]²+α_T[g(i,j)−f(i,j)]²

where f_MC(i,j) represents a motion compensated pixel, α_{HL, αHR}, α_VT, α_VD and α_T represent a regulation parameter variable with respect to each cost function.

14. The method of claim 13, wherein the pixel f(i,j) which will be recovered is obtained based on the following equation when the pixel is included in an inter macro block,

$f\left(i,j\right)=\frac{\begin{array}{c}f\left(i,j-1\right)+f\left(i,j+1\right)+f\left(i-1,j\right)+\\ f\left(i+1,j\right)+f_{\mathrm{MC}}\left(i,j\right)+{\alpha }_{\mathrm{TOT}}g\left(i,j\right)\end{array}}{5+{\alpha }_{\mathrm{TOT}}}$

where, α_TOT=α_HL+α_HR+α_VT+α_VD+α_T, and the pixel f(i,j) which will be recovered is obtained based on the following equation when the pixel is included in an intra macro block,

$f\left(i,j\right)=\frac{\begin{array}{c}f\left(i,j-1\right)+f\left(i,j+1\right)+f\left(i-1,j\right)+\\ f\left(i+1,j\right)+{\alpha }_{\mathrm{TOT}}g\left(i,j\right)\end{array}}{4+{\alpha }_{\mathrm{TOT}}}$

where α_TOT=α_HL+α_HR+α_VT+α_VD.

15. The method of claim 13, wherein said regularization parameter variables α_HL, α_HR, α_VT, α_VD, α_T are obtained by approximations as follows:

${\alpha }_{\mathrm{HL}}=\frac{{\left[g\left(i,j\right)-g\left(i,j-1\right)\right]}^2}{Q_{\mathrm{pl}}^2},{\alpha }_{\mathrm{HR}}=\frac{{\left[g\left(i,j\right)-g\left(i,j+1\right)\right]}^2}{Q_{\mathrm{pl}}^2}$ ${\alpha }_{\mathrm{VT}}=\frac{{\left[g\left(i,j\right)-g\left(i-1,j\right)\right]}^2}{Q_{\mathrm{pl}}^2},{\alpha }_{\mathrm{VD}}=\frac{{\left[g\left(i,j\right)-g\left(i+1,j\right)\right]}^2}{Q_{\mathrm{pl}}^2}$ ${\alpha }_T=\frac{{\left[g\left(i,j\right)-f_{\mathrm{MC}}\left(i,j\right)\right]}^2}{Q_{\mathrm{pl}}^2},$

where Q_pl represents a quantizing variable of the l-th macro block.

16. The method of claim 11, wherein in said step for obtaining the projected pixel P(F(u,v)), when (u,v)-th value F(u,v) of two-dimensional DCT coefficient of the original image is smaller than G(u,v)−Qpl, the projected pixel P(F(u,v)) is mapped to G(u,v)−Qpl, and when the value F(u,v) is larger than G(u, v)+Qpl, the projected pixel P(F(u, v)) is mapped to G(u, v)+Qpl, otherwise, the projected pixel P(F(u,v)) is mapped to F(u,v), where G(u,v) represents (u,v)th value of the two-dimensional DCT coefficient of the compression image, and Qpl represents a quantizing maximum difference of the l-th macro block.

17. The method of claim 11, further comprising the following steps which are repeatedly performed by k-times:

defining a cost function M(i,j) having a smoothing degree of an image and a reliability with respect to the original image by the unit of pixels in consideration with a directional characteristic between the pixels which will be recovered and the pixels neighboring the pixels which will be recovered;

adaptively searching a regularization parameter variable having a weight value of a reliability with respect to the original image from the cost function M(i,j); and

obtaining a projected pixel P(F(u,v) using a projection method for mapping the recovering pixel in accordance with a range value of the pixel which will be recovered, for thereby finally obtaining a recovering image.

18. In a method for recovering a compressed motion image for processing an original pixel f(i,j) based on a DCT by the unit of macro blocks of a M×M size, quantizing the DCT-processed coefficient, transmitting together with motion vector information, reversely quantizing and reversely DCT-processing the compressed pixel g(i,j) and recovering an image similar to the original image, a method for recovering a compressed motion picture, comprising the steps of:

defining a cost function M(i,j) having a smoothing degree of an image and a reliability with respect to an original image as a pixel unit in consideration of a directional characteristic between the pixels which will be recovered and the pixels neighboring the pixels which will be recovered;

adaptively searching a regularization parameter variable having a weight of a reliability with respect to the original image from the cost function M(i,j); and

obtaining a finally recovered image of a spatial region by obtaining a block DCT coefficient based on a block DCT and obtaining a projected pixel P(F(u,v)) by a projection method for mapping the pixels which will be recovered in a range value of the pixel for processing the block DCT coefficient, and performing a reverse DCT,

wherein said regularization parameter variable is a weight value with respect to reliability and is determined based on a difference between the original pixel and the compressed pixel and a difference value between the original pixel and the neighboring pixel.

19. An apparatus for recovering a compressed motion picture, comprising:

an image decoding unit for outputting an information with respect to an image which will be recovered such as a decoded image, a quantized variable, a macro block type, and a motion type by decoding a coded image signal; and

a block process eliminating filter for defining a cost function based on a smoothing degree of an image and a reliability with respect to an original pixel in consideration of a directional characteristic between the neighboring pixel and the pixel which will be processed based on the pixels which will be recovered using an information with respect to the image which will be recovered inputted from the image decoding unit, adaptively searching a regularization parameter variable which provides a weight of a reliability with respect to the original image for each cost function, and recovering an original pixel using a projection method for mapping the pixels which will be recovered in accordance with a range value of the pixels which will be processed,

wherein said regularization parameter variable is a weight value with respect to reliability and is determined based on a difference between the original pixel and the compressed pixel and a difference value between the original pixel and the neighboring pixel.

20. The apparatus of claim 19, further comprising:

a DCT unit for performing a DCT with respect to an image recovered by the block process eliminating filter;

a vector projection unit for projecting a pixel which will be recovered in accordance with a pixel value after the DCT process is performed; and

an IDCT unit for performing a reverse DCT with respect to the image projected by the vector projection unit.

21. In a method for recovering a compressed motion image for processing an original pixel f(i,j) based on a DCT by the unit of macro blocks of a M×M size, quantizing the DCT-processed coefficient, transmitting together with motion vector information, reversely quantizing and reversely DCT-processing the compressed pixel g(i,j) and recovering an image similar to the original image, a method for recovering a compressed motion picture, comprising the steps of:

defining a cost function M(i,j) having a smoothing degree of an image and a reliability with respect to an original image as a pixel unit in consideration with a directional characteristic between the pixels which will be recovered and the pixels neighboring the pixels which will be recovered; and

adaptively searching a regularization parameter variable having a weight of a reliability with respect to the original image from the cost function M(i,j) and a weight value of a smoothing degree of the original image,

wherein said regularization parameter variable is a weight value with respect to reliability and is determined based on a difference between the original pixel and the compressed pixel and a difference value between the original pixel and the neighboring pixel.

22. The method of claim 21, wherein said cost function is obtained based on the following equations:

M_L(f(i,j))=α_L(f(i,j))[f(i,j)−f(i,j−1)]²+(1−α_L(f(i,j)))[g(i,j)−f(i,j)]²

M_R(f(i,j))=α_R(f(i,j))[f(i,j)−f(i,j+1)]²+(1−α_R(f(i,j)))[g(i,j)−f(i,j)]²

M_U(f(i,j))=α_U(f(i,j))[i f(i,j)−f(i−1i,j)]²+(1−α_U(f(i,j)))[g(i,j)−f(i,j)]²

M_D(f(i,j))=α_D(f(i,j))[f(i,j)−f(i−1,j)]²+(1−α_D(f(i,j)))[g(i,j)−f(i,j)]²

where α_L, α_R, α_U, α_D are regularization parameter variables with respect to each cost function.

23. The method of claim 22, wherein when the pixel of the current macro block is the same as the pixel of the previously transmitted macro block, in said pixel f(i,j) which will be recovered, the pixel value which is previously recovered with respect to the macro block of the previous image is substituted for the current pixel value, and otherwise the following Equation is obtained:

$f\left(i,j\right)=\frac{\begin{array}{c}{\alpha }_{L}f\left(i,j-1\right)+{\alpha }_{R}f\left(i,j+1\right)+{\alpha }_{U}f\left(i-1,j\right)+\\ {\alpha }_{D}f\left(i+1,j\right)+\left(4-{\alpha }_{\mathrm{TOT}}\right)g\left(i,j\right)\end{array}}{4}$

where α_TOT=α_L+α_R+α_U+α_D.

24. The method of claim 22, wherein said regularization parameter variables α_L, α_R, α_U, α_D are approximated as follows:

${\alpha }_L=\frac{K_L{Q}_P^2}{{\left[g\left(i,j\right)-g\left(i,j-1\right)\right]}^2+K_L{Q}_P^2}$ ${\alpha }_R=\frac{K_R{Q}_P^2}{{\left[g\left(i,j\right)-g\left(i,j+1\right)\right]}^2+K_R{Q}_P^2}$ ${\alpha }_U=\frac{K_U{Q}_P^2}{{\left[g\left(i,j\right)-g\left(i-1,j\right)\right]}^2+K_U{Q}_P^2}$ ${\alpha }_D=\frac{K_D{Q}_P^2}{{\left[g\left(i,j\right)-g\left(i+1,j\right)\right]}^2+K_D{Q}_P^2}$

where K_LQp², K_RQp², K_UQp², K_DQp² are functions of the quantizing variable Qp, and constants K_L, K_R, K_U, K_D are weight values with respect to the regularization parameter variables α_L, α_R, α_U, α_D, and have different values based on whether the neighboring pixel is positioned at the block boundary or in the interior of the block.

25. The method of claim 24, wherein the weight values K_L, K_R, K_U, K_D are expressed as follows, assuming that i and j of the pixel f(i,j) are 8, respectively,

K_L={9, if j mod 8=0; 1, otherwise}

K_R={9, if j mod 8=7; 1, otherwise}

K_U={9, if i mod 8=0; 1, otherwise}

K_D={9, if i mod 8=7; 1, otherwise}.

26. An apparatus for recovering a compressed motion picture, comprising:

an image decoding unit for outputting an information with respect to an image which will be recovered, a quantized variable, a macro block type, and a motion type by decoding a coded image signal; and

a block process eliminating filter for defining a cost function based on a smoothing degree of an image and a reliability with respect to an original pixel in consideration of a directional characteristic between a neighboring pixel and the pixel which will be processed based on the pixels which will be recovered using an information with respect to the image which will be recovered inputted from the image decoding unit, and adaptively searching a regularization parameter variable which has a weight of a reliability with respect to the original image from each cost function and a weight of a smoothing degree of the original image for thereby recovering an original pixel,

wherein said regularization parameter variable is a weight value with respect to reliability and is determined based on a difference between the original pixel and the compressed pixel and a difference value between the original pixel and the neighboring pixel.

27. A method of filtering an image, comprising:

determining, with a filter apparatus, a filter strength based on whether a block including a pixel to be filtered is intra-coded; and

filtering, with the filter apparatus, the pixel to be filtered using a filtering methodology that adjusts a degree of filtering based on the determined filter strength and a difference value, the difference value being based on the pixel to be filtered and a neighboring pixel.

28. The method of claim 27, wherein the neighboring pixel is a pixel adjacent to the pixel to be filtered.

29. The method of claim 27, further comprising:

determining the difference value.

30. The method of claim 29, wherein the determining step determines the difference between the pixel to be filtered and the neighboring pixel as the difference value.

31. The method of claim 27, wherein the filtering step filters the pixel to be filtered based on a quantization parameter used in processing a portion of an image including the pixel to be filtered.

32. The method of claim 31, wherein the portion of an image including the pixel to be filtered is a macroblock.

33. The method of claim 27, wherein the filtering methodology includes determining at least one boundary value based on a quantization parameter of a portion of the image including the pixel to be filtered.

34. The method of claim 33, wherein the portion of the image including the pixel to be filtered is a macroblock.

35. The method of claim 33, wherein the filtering methodology includes determining more than one boundary value based on the quantization parameter.

36. The method of claim 35, wherein the portion of the image including the pixel to be filtered is a macroblock.

37. A method of filtering an image, comprising:

determining, with a filter apparatus, a filter strength based on whether a block including a pixel to be filtered is intra-coded;

obtaining, with the filter apparatus, a difference between values of two pixels in a portion of an image; and

filtering, with the filter apparatus, the pixel to be filtered using a filtering methodology that adjusts a degree of filtering based on the difference and the determined filter strength.

38. The method of claim 37, wherein the two pixels are the pixel to be filtered and a neighboring pixel near the pixel to be filtered.

39. The method of claim 37, wherein the portion of the image is an image block.

40. The method of claim 39, wherein the image block is a macro block.

41. The method of claim 37, wherein the degree of filtering is an amount of filtering.

42. The method of claim 27, wherein the degree of filtering is an amount of filtering.

43. The method of claim 27, further comprising:

outputting a filtered image including the filtered pixel.

44. The method of claim 37, further comprising:

outputting a filtered image including the filtered pixel.

Images