Patent Publication Number: US-5841909-A

Title: Method of generating an orthogonal basis function set in an image processing system

Description:
This application is a continuation of application Ser. No. 08/266,765, filed Jun. 28, 1994, now abandoned. 
    
    
     BACKGROUND OF THE INVENTION 
     1. Field of the Invention 
     The present invention relates generally to improvements in a digital image processing system and more specifically to a method of coding arbitrarily shaped image segments wherein calculation errors can markedly be reduced. The present invention is also applicable to a decoder for reproducing an original image signal while holding the same features of the present invention. 
     2. Description of the Prior Art 
     Digital image coding has found effective applications ranging from medical imagery, remote sensing up to image communication such as video conference and picture telephone. 
     Prior to image coding it is necessary to divide each of incoming digital frame images into a plurality of segments. Typical image segmentation is to partition the frame image into rectangular blocks (viz., block-segmentation), after which each block is coded separately. 
     The block-segmentation has several advantages of no transmission of shape or contour information due to fixed block-size and low computational complexity. However, the independent treatment of each block often results in a visibility of the block boundaries in the reproduced Image, called blocking-effect, due to discontinuities in the gray values between adjacent blocks. 
     In recent years, in order to overcome the problems inherent in the block-segmentation, a region-oriented segmentation has been proposed. One example of which is described in a paper (PRIOR PAPER 1) by Michael GILGE, et al., entitled &#34;CODING OF ARBITRARILY SHAPED IMAGE SEGMENTS BASED ON A GENERALIZED ORTHOGONAL TRANSFORM&#34; published 1989 by Elsevier Science Publishers B.V. 
     Before turning to the present invention it is deemed preferable to briefly discuss an image signal encoder with reference to FIG. 1. 
     The present invention is directly concerned with an orthogonal basis function generator 14 of FIG. 1, and as such, the other blocks 10, 12, 16, and 18 will not be discussed in detail merely for the sake of simplifying the disclosure. 
     In FIG. 1, a digital image signal S1 is successively applied to the arbitrarily shaped image segmentation circuit 10 on a frame basis. 
     The circuit 10 stores one frame image data in a suitable memory (not shown) and then carries out image segmentation utilizing techniques such as measurement space guided clustering or single linkage region growing. These techniques are discussed in detail in a paper (PRIOR PAPER 2) by Robert M. HARALICK, et al. entitled &#34;Survey: Image Segmentation Techniques&#34;, computer vision, graphics, and image processing 29, 100-132 (1985), published by Academic Press, Inc. 
     The above mentioned clustering technique defines a partition in measurement space. Then each pixel is assigned the label of the cell in the measurement space partition to which it belongs. The image segments are defined as the connected components of the pixels having the same label. The image segmentation circuit 10 outputs a shape signal S2 and a segment content signal S3. 
     The contour coding circuit 12 is supplied with the shape signal S2 and codes a contour thereof utilizing a chain coding scheme (for example) described in a paper by Herbert Freeman, entitled &#34;Computer Processing of Line-Drawing Images&#34;, Computing Surveys, Vol. 6, No. 1, March 1974 (PRIOR PAPER 3). The circuit 12 outputs a contour code which is applied to the multiplexer 18. 
     The orthogonal basis function generator 14 receives the shape signal S2 and outputs a basis function set which is applied to the transform coefficient generator 16. A set of basis functions each of which is orthogonal with respect to the shape of the segment to be coded can be obtained using orthogonalization schemes. The generator 14 will be discussed in detail later and thus, further descriptions thereof are not given here. 
     The transform coefficient generator 16 receives the content signal S3 which represents one segmented data, and thereafter, outputs transform coefficients after orthogonally coding the content signal S3 using the basis function set applied from the generator 14. 
     One of the features of the coder shown in FIG. 1 is that the transmission of the basis functions for each segment is no longer required as opposed to the block segmentation as mentioned in PRIOR PAPER 1. 
     A known orthogonal basis function set generator 14 will be discussed. 
     As is well known in the art, a Gram-Schmidt process is to generate an orthonormal basis function from a given basis function. That is, an orthogonal basis function Y 
     
         Y={b.sub.1, b.sub.2, . . . , b.sub.i, . . . , b.sub.n }    (1) 
    
     is generated from a given basis function X 
     
         X={a.sub.1, a.sub.2, . . . , a.sub.i, . . . , a.sub.n }    (2) 
    
     using the following relationship: ##EQU1## Normalizing the orthogonal basis function given in expression (1) yields an orthonormal basis function. 
     According to the Gram-Schmidt process, the second term vector (depicted by /a i ) of expression (3) ##EQU2## is a shadow of the vector a i  on a hyperplane spanned by the vector b j  (j=1, 2, . . . , i-1) as shown in FIG. 2 wherein /a i  is represented by a i  with a hat &#34;&#34;. Therefore, a vector e i  which indicates the difference between the vectors /a i  and a i , is a component which is orthogonal with respect to the above mentioned hyperplane spanned by the vector b j . The vector e i  is then normalized. 
     It is assumed that the image segmentation circuit 10 (FIG. 1) outputs a reverse L-shaped block (segment) as shown in FIG. 3, which is applied to the basis function set generator 14. As shown, the reverse L-shaped block, depicted by RLB, includes five pixels. 
     Firstly, the generator 14 determines nine (9) orthogonal basis functions (FIG. 4) with respect to a rectangle region, formed of 3×3 pixels, which includes the reverse L-shaped block RLB (FIG. 3). The nine orthogonal basis functions are generated using DCT (Discrete Cosine Transform) functions by way of example. In FIG. 4, bold arrows indicate a so-called zigzag scan of the nine pixels within the rectangle of FIG. 3, while the parenthesized numeral at the top of the left of each of the nine blocks indicates the order of the scan. 
     FIG. 5 is a table showing nine (9) sets of basis functions, each set of which is obtained by extracting the reverse L-shaped block RLB (FIG. 3) from the corresponding function set of FIG. 4 and then normalizing same. Merely for the convenience of description, each of the basis function sets shown in FIG. 5 is called a basis function candidate in the instant disclosure. 
     Although the number of the basis function candidate sets shown in FIG. 5 is nine, five basis function candidate sets are sufficient in that the number of the pixels within the reverse L-shaped block RLB is five in this particular case. Therefore, the five sets of the basis function candidates  1!- 5! along the zigzag scan are utilized. Thus, the following matrix can be obtained wherein the elements of the first column: 
     0.4472, 0.5000, 0.5000, 0.3536, and 0.5774 
     correspond to the leftmost vectors of the uppermost lines of the basis function candidates  1!- 5! of FIG. 5, respectively. Similarly, the second column&#39;s elements correspond to the middle vectors of the first lines of the basis candidate sets  1!- 5!, and so on. ##EQU3## 
     It is understood that the determinant of expression (5) equals zero (0) while the rank thereof is 4. This implies that the basis function candidates as discussed above are not linearly independent. 
     FIG. 6 is a table showing five basis functions which are obtained from the basis function candidates of FIG. 5 using equation (3) by way of 64-bit double precision floating point arithmetic operations. In FIG. 6, each of the numbers parenthesized (viz.,  1!- 5!) indicates that the basis function with the numeral is obtained from the basis function candidate with the same number of FIG. 5. 
     An inner product matrix obtained from the basis functions of FIG. 6 is shown in FIG. 7. The matrix of FIG. 7 indicates that the No. 5 basis function is not orthogonal with respect to the other basis functions Nos. 1-4. 
     The vector e i  is relatively large as shown in FIG. 2 and thus, any undesirable result raised from the fact that the vector e i  is small is not noticeable in the above case. However, if the orthogonal component (viz., the vector e i ) is small as shown in FIG. 8, calculation errors are liable to occur. That is to say, if a vector a i  to be orthogonalized has a component which is perpendicular to the hyperplane spanned by the following vector Y&#39; which has already been orthogonalized 
     
         Y&#39;={b.sub.1, b.sub.2, . . . , b.sub.i-1 }                  (7) 
    
     then vector a i  can be approximated by the vector /a i  of equation (4) and thus, the orthogonal component vector e i  becomes small. As a result, when the vector b i  is normalized using 
     
         b.sub.i =e.sub.i /|e.sub.i |             (8) 
    
     there exists the possibility that a basis function, which is not orthogonal with respect to the hyperplane spanned by Y&#39;, nay be generated. In other words, calculation errors occur. 
     SUMMARY OF THE INVENTION 
     It is therefore an object of the present invention to provide a method of generating an orthogonal basis function set wherein calculation errors can be extensively reduced. 
     In brief, the above object is achieved by a method of generating a set of orthogonal basis functions in an image processing system. A contour signal of an arbitrarily shaped image segment is generated. Thereafter, a set of orthonormal basis functions is determined with respect to pixels within a rectangle which includes said arbitrarily shaped image segment. Further, a set of orthogonal basis function candidates is determined using a set of the orthonormal basis functions. The orthogonal basis function candidates are numbered along a pixel scan. 
    
    
     BRIEF DESCRIPTION OF THE DRAWINGS 
     The features and advantages of the present invention will become more clearly appreciated from the following description taken in conjunction with the accompanying drawings in which like elements are denoted by like reference numerals and in which: 
     FIG. 1 is a block diagram schematically showing an encoder of a digital image processing system to which the present invention is applicable, this figure having been described in the opening paragraphs of the instant disclosure; 
     FIG. 2 is a sketch for describing the operations of one block of FIG. 1; 
     FIG. 3 shows an arbitrarily shaped segment (hatched portion) for discussing the operation of the arrangement of FIG. 1; 
     FIG. 4 is a table showing nine orthogonal basis functions with respect to a rectangle which includes the arbitrarily shaped segment of FIG. 3, having been referred to in the opening paragraphs of the instant disclosure; 
     FIG. 5 is a table showing nine orthogonal basis function candidates with respect to the arbitrarily shaped segment of FIG. 3, having been referred to in the opening paragraphs of the instant disclosure; 
     FIG. 6 is a table showing five orthogonal basis functions with respect to the arbitrarily shaped segment of FIG. 3, this figure having been referred to in the opening paragraphs of the instant disclosure; 
     FIG. 7 is an inner product matrix obtained from the orthogonal basis functions of FIG. 6; 
     FIG. 8 is a sketch for describing the operation of one block of FIG. 1; 
     FIG. 9 is a flow chart which depicts steps which characterize a first embodiment of the present invention; 
     FIG. 10 is a table showing five orthogonal basis functions obtained with the first embodiment; 
     FIG. 11 shows an inner product matrix obtained from the orthogonal basis functions of FIG. 10; 
     FIG. 12 is a flow chart which depicts steps which characterize a second embodiment of the present invention; 
     FIG. 13 is a table showing five orthogonal basis functions obtained with the second embodiment; 
     FIG. 14 shows an inner product matrix obtained from the orthogonal basis functions of FIG. 13; 
     FIG. 15 is a flow chart which depicts steps which characterize a third embodiment of the present invention; 
     FIG. 16 is a table showing five orthogonal basis functions obtained with the third embodiment; and 
     FIG. 17 shows an inner product matrix obtained from the orthogonal basis functions of FIG. 13. 
    
    
     DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
     FIG. 9 is a flow chart which depicts the operation which characterizes the first embodiment of the present invention. 
     As mentioned above, the present invention is directly concerned with the operation at the orthogonal basis function set generator 14 of FIG. 1. 
     Merely for the convenience of description, it is assumed, as in the opening paragraphs of the instant disclosure, that an arbitrarily shaped segment (viz., reverse L-shaped block) of FIG. 3 has been applied to the generator 14 from the image segmentation circuit 10 of FIG. 1. 
     In FIG. 9, step 100 of the routine illustrated, comprises the operation wherein a set of the nine orthonormal basis functions shown in FIG. 4 is determined in connection with the rectangle which includes five pixels as shown in FIG. 3. Following this, at step 101 a set of the nine orthogonal basis function candidates is determined based on the orthonormal basis functions which are defined by the region of the reverse L-shaped block (FIG. 5). At step 102 the basis function candidates thus obtained are numbered along the zigzag-like scan as indicated by  1!- 9! in FIG. 5. At step 103, the first basis function candidate (depicted by  1! in FIG. 5) is extracted and defined as the first orthogonal basis function, and the routine flows to step 104 wherein the next basis function candidate (depicted by  2! at this time) is extracted. 
     At step 105, determined is an orthogonal component which is perpendicular to a hyperplane spanned by the basis function which has been determined. This operation is implemented using the known Gram-Schmidt orthogonalizing process. Subsequently, at step 106 it is determined if the absolute value of the orthogonal component obtained at step 105 exceeds a threshold value. At step 106, a basis function candidate, whose orthogonality is small (viz., the orthogonal component e i  does not exceed the threshold value), can effectively be canceled. Further, a basis function candidate, which Is not linearly independent (viz., the orthogonal component is zero), can be omitted. 
     If the answer is not affirmative (viz., NO) at step 106, the routine goes back to step 104 wherein the next orthogonal basis function candidate  3! is extracted and the above mentioned operation is repeated. Contrarily, if the answer is affirmative at step 106, at step 107 the orthogonal component obtained at step 105 is adopted as the next orthogonal basis function at step 105 (if the routine goes to this step), Following this, it is determined if the number of the orthogonal basis functions obtained is equal to the number of the pixels in the arbitrarily shaped segment (viz., the reverse L-shaped block). In the event that the answer is affirmative at step 108, the program is terminated. Otherwise, the routine returns to step 104. 
     At step 106, if the i-th basis function candidate is depicted as a i  and if the basis functions b k  (k=1, 2, . . . ,j (j&lt;i)) have been determined, then the orthogonal component e i  is given by ##EQU4## It the vector e i  does not exceed the predetermine value, it is canceled at step 106 as mentioned above. 
     The threshold value used at step 106 is determined considering what value is appropriate in order to achieve effective reduction of calculation errors without adversely affecting the operation results. 
     FIG. 10 is a table showing five orthogonal basis functions obtained, using the first embodiment, with respect to the reverse L-shaped segment shown in FIG. 3. In FIG. 10, each of the numerals parenthesized (viz.,  1!- 4! and  6!) at the top portions of the orthogonal basis functions, indicates the number when the zigzag-like scan is implemented as shown in FIGS. 4 and 5. It is understood that the fifth basis function candidate is not selected at step 106 in this case. Instead, the sixth basis function candidate is selected. 
     FIG. 11 shows an inner product matrix obtained from the orthogonal basis functions shown in FIG. 10. It is understood that the inner product matrix of FIG. 11 is a unit matrix. 
     FIG. 12 is a flow chart which depicts the operation which characterizes the second embodiment of the present invention. 
     As in the first embodiment, an arbitrarily shaped segment (viz., reverse L-shaped block) of FIG. 3 has been applied to the generator 14 from the image segmentation circuit 10 of FIG. 1. 
     In FIG. 12, steps 110 and 111 are same as steps 100 and 101 of FIG. 9. That is to say, at step 110, a set of the nine orthonormal basis functions shown in FIG. 4 is determined in connection with the rectangle which includes five pixels shown in FIG. 3. Following this, at step 111, a set of the nine orthogonal basis function candidates is determined based on the orthonormal basis functions which are defined by the region of the reverse L-shaped block (FIG. 5). 
     At step 112, one of the basis function candidates is extracted and will be used as the first orthogonal basis function at step 113. At step 113, an orthogonal component is determined which is perpendicular to a hyperplane spanned by the basis function, which has been determined, with respect to each of the remaining basis function candidates. The operation at step 113 is implemented using the known Gram-Schmidt orthogonalizing process. 
     At step 114, the maximum orthogonal component among those obtained at step 113 is selected and also used as the next orthogonal basis function at step 113 (if the routine goes back to step 113). 
     At step 115 it is determined if the number of the orthogonal basis functions obtained is equal to the number of the pixels in the arbitrarily shaped segment (viz., the reverse L-shaped block). If the answer is affirmative at step 115, the program is terminated. Otherwise (viz., the answer is NO), the routine returns to step 113 wherein the orthogonal component of each of all the remaining basis function candidates is determined. That is, the number of the basis function candidates is reduced by one every time the one looped operations at steps 113-115 is completed. 
     FIG. 13 is a table showing five orthogonal basis functions obtained with respect to the reverse L-shaped segment shown in FIG. 3, using the second embodiment. As shown in FIG. 13, the orthogonal basis function candidates  1!,  9!,  5!,  2!, and  8! are extracted in this order and then utilized to obtain the five orthogonal basis functions. 
     FIG. 14 shows an inner product matrix obtained from the orthogonal basis functions shown in FIG. 13. It is understood that the inner product matrix of FIG. 11 is a unit matrix. 
     Finally, the third embodiment of the present invention will be discussed with reference to a flow chart shown in FIG. 15. 
     The flow chart of FIG. 15 includes steps 120 and 121 in addition to steps 110-115 of FIG. 12. Thus, steps 110-115 will not be described for the sake of simplifying the disclosure. Step 120 corresponds to step 102 of FIG. 9 (first embodiment). That is, at step 120 the basis function candidates obtained at step 111 are numbered along the zigzag-like scan as indicated by  1!- 9! in FIG. 5. Further, at step 121, the orthogonal basis functions are rearranged in the order of the numbered basis function candidates. 
     FIG. 16 is a table showing five orthogonal basis functions obtained with respect to the reverse L-shaped segment, using the third embodiment. The orthogonal basis function candidates correspond to those shown in FIG. 13 but arranged in the order of  1!,  2!,  5!,  8!, and  9!. 
     FIG. 17 shows an inner product matrix obtained from the orthogonal basis functions shown in FIG. 16. It is understood that the inner product matrix of FIG. 11 is a unit matrix. 
     It will be understood that the above disclosure is representative of three possible embodiments of the present invention and that the concept on which the invention is based is not specifically limited thereto.