Patent Publication Number: US-2005141778-A1

Title: Image processing method, image processing apparatus and image processing program

Description:
FIELD OF THE INVENTION  
      The present invention relates to an image processing method, image processing apparatus and image processing program.  
     BACKGROUND OF THE INVENTION  
      Conventionally, in a portrait photographing using a silver halide film in a portrait studio, a correction processing is applied to a print or a negative film by a manual working so that a puckering or spotting of a face portion or a neck portion of a person is lightened or erased.  
      Further, recently, a development of a digital still camera (hereinafter, called DSC) is conspicuous, and a single-lens reflex type DSC having a resolving power not smaller than 10 7  pixels is put on the market, and also in the portrait photographing in the portrait studio, by using a high resolving power DSC, the photographing is conducted.  
      Accordingly, in Patent Document 1, an image processing method by which a removal of a puckering or spotting portion can be effectively applied to the face of a person included in an image made by DSC, is proposed.  
      (Patent Document 1) Tokkai 2003-209683  
      However, in the method of Patent Document 1, an operator indicates and inspects a correction of red-eye or a removal of a puckering or spotting after an extraction of a face area included in the image. Therefore, the productivity is low, and depending on the operator&#39;s inspection, there is a possibility that the sharpness of the whole image is deteriorated, or an outline structure of the face of a person is changed, and the image becomes an image having a sense of incompatibility.  
     SUMMARY OF THE INVENTION  
      A problem of the present invention is to make it possible to remove easily the puckering or spotting without the sense of incompatibility when the puckering or spotting is removed from a face or neck of the portrait image, without deteriorating the sharpness feeling of the whole image, or without changing the outline structure of a face of a person.  
      To solve the above problem, the present invention provides an image processing method, apparatus and program in which an image signal expressing a color image is obtained, the pixels forming the skin-colored area are extracted from the pixels included in the image signal, and when the pixels has the spatial frequency within a specific range and the signal intensity whose variation amount is not larger than specific threshold value among the extracted pixels, a variation amount of the signal intensity of the pixels are reduced. 
    
    
     BRIEF DESCRIPTION OF THE DRAWINGS  
       FIG. 1  is a view showing the wavelet function used in the multiresolution transform according to the present invention.  
       FIG. 2  is a view showing a flow of a filter processing of 1-level wavelet transform.  
       FIG. 3  is a system block diagram showing a flow of a filter processing of 1-level wavelet transform.  
       FIG. 4  is a typical view showing a process in which a signal is decomposed by 3-level wavelet transform.  
       FIG. 5  is a system block diagram for re-structuring a signal before decomposition by the filter processing of the inverse wavelet transform.  
      FIGS.  6 ( a )- 6 ( e ) are views showing a waveform of an input signal S 0  and a waveform of corrected higher frequency band component W i -γ i  of each level obtained by the wavelet transform.  
       FIG. 7  is a system block diagram showing a flow of the filter processing of a 1-level dyadic wavelet transform.  
       FIG. 8  is a view showing a flow of the filter processing of a 1-level inverse dyadic wavelet transform.  
       FIG. 9  is a system block diagram showing a flow of a processing from the dyadic wavelet transform to a time point when the signal S 0 ′ on which the image processing is applied, is obtained.  
       FIG. 10  is a perspective view showing an outline structure of an image recording apparatus.  
       FIG. 11  is a block diagram showing an internal structure of the image recording apparatus.  
       FIG. 12  is a block diagram showing a functional structure of an image processing section of  FIG. 11 .  
       FIG. 13  is a system block diagram according to the internal processing of an image adjustment processing section in Example 1.  
       FIG. 14  is a view showing an image evaluation result when a plurality of image processing in which image processing conditions in Example 1 are different, are applied.  
       FIG. 15  is a system block diagram according to the internal processing of the image adjustment processing section in Example 2. 
    
    
     DETAILED DESCRIPTION OF THE INVENTION  
      Preferred embodiments of the present invention will be described below.  
      To solve the above problem, an embodiment written in item 1 is an image processing method including an obtaining step of obtaining an image signal representing a color image, an extracting step of extracting pixels which form a skin colored area from pixels which are included in the obtained image signal, a signal intensity adjusting step of reducing a variation amount of a signal intensity of a pixel whose spatial frequency is in a predetermined range and whose signal intensity variation amount is a predetermined threshold or below among the extracted pixels.  
      Further, as an embodiment written in item 2, in the embodiment written in item 1, it is preferable that the extracting step extracts pixels which form a skin colored area by applying a dyadic wavelet transform to the obtained image signal and applying a simple area expansion to a lower frequency band component of the transformed image signal.  
      As an embodiment written in item 3, in the embodiment written in item 1 or item 2, it is preferable that the spatial frequency within the predetermined range is obtained by a higher frequency band component calculated by a multiresolution transform.  
      As an embodiment written in item 4, in the embodiment written in item 3, it is preferable that the signal intensity adjusting step determines the predetermined threshold based on a statistic value calculated by the higher frequency band component of an image signal intensity.  
      As an embodiment written in item 5, in the embodiment written in item 4, it is preferable that the statistic value is one of a standard deviation, an average, a median and a mode.  
      As an embodiment written in item 6, in the embodiment written in any one of items 1-5, it is preferable that when a pixel among the extracted pixels has a spatial frequency of 1.5-3.0 lines/mm and a signal intensity variation is in a range between 0-20% of a maximum signal intensity, the signal intensity adjusting step reduces a variation amount of a signal intensity of the pixel by 0.5-0.75 times.  
      As an embodiment written in item 7, in the embodiment written in any one of items 1-6, it is preferable that the image processing method further comprise: a color information converting step of converting the obtained image signal into a brightness signal and a color-difference signal, wherein the signal intensity regulating step is applied to the brightness signal and/or the color-difference signal of the extracted pixels which form a skin colored area.  
      As an embodiment written in item 8, in the embodiment written in item 7, it is preferable that the image processing method further comprise: a contrast controlling step of controlling a contrast of an image signal of the extracted pixels which form a skin colored area and whose image signal is reduced in the signal intensity adjusting step.  
      As an embodiment written in item 9, in the embodiment written in item 7 or item 8, it is preferable that the image processing method further comprise: a noise adding step of adding a noise signal to the image signal of the extracted pixels which form a skin colored area and whose image signal is reduced by the signal intensity regulating step.  
      To solve the above problem, an embodiment written in item 10 is an image processing method of item 1 further including a color information converting step of converting the obtained image signal into a brightness signal and a color-difference signal; and a wavelet transform step of decomposing the converted image signal into image signals with different frequency band components by applying at least a second level dyadic wavelet transform to the converted image signal, wherein when an image signal intensity in at least a higher frequency band component at level 2 of the extracted pixels among the decomposed frequency band components is predetermined threshold value or below, the image intensity adjusting step reduces the signal intensity.  
      In other words, the embodiment written in item 10 is an image processing method including an obtaining step of obtaining an image signal representing a color image; an extracting step of extracting pixels which form a skin colored area from pixels which are included in the obtained image signal; a color information converting step of converting the obtained image signal into a brightness signal and a color-difference signal; a wavelet transform step of decomposing the converted image signal into image signals with different frequency band components by applying at least a second level dyadic wavelet transform to the converted image signal; and a signal intensity adjusting step of reducing a variation amount of a signal intensity of a pixel whose spatial frequency is in a predetermined range and whose signal intensity variation is a predetermined threshold or below among the transformed pixels.  
      As an embodiment written in item 11, in the embodiment written in items 10, it is preferable that the wavelet transform step applies the dyadic wavelet transform to the brightness signal and/or the color-difference signal.  
      As an embodiment written in item 12, in the embodiment written in item 11, it is preferable that the image processing method further comprise a contrast controlling step of controlling a contrast of a image signal in the lower frequency band component of the extracted pixels which form a skin colored area among the image signals with different frequency band components decomposed by at least a first level of dyadic wavelet transform in the wavelet transform step.  
      As an embodiment written in item 13, in the embodiment written in item 11 or item 12, it is preferable that the image processing method further comprise: a noise adding step of adding a noise signal to the image signal of the extracted pixels which form a skin colored area and whose image signal is reduced by the signal intensity regulating step.  
      To solve the above problem, an embodiment written in item 14 is an image processing method of item 1 further comprise: a first transform step and a second transform step, wherein the first transform step decomposes the obtained image signal into image signals with different frequency band components by applying at least the first level multiresolution transform which is a method of reducing an image size, the extracting step extracts pixels which form skin colored area from pixels included in the image signal with a lower frequency band component among image signals with different frequency band components decomposed by the multiresoluting transform, the second transform step decomposes the image signal with the lower frequency band component into image signals with different frequency band components using at least a first level dyadic wavelet transform, and when an intensity of an image signal in a higher frequency band component of the extracted pixels which form a skin colored area among image signals in different frequency band components decomposed by the dyadic wavelet transform, the signal intensity adjusting step reduces the signal intensity.  
      In other words, the embodiment written in item 14 is an image processing method including:  
      an obtaining step of obtaining an image signal representing a color image; a first transform step of decomposing the obtained image signal into image signals with different frequency band components by applying at least the first level multiresolution transform which is a method of reducing an image size; an extracting step of extracting pixels which form a skin colored area from pixels included in the image signal in a lower frequency band component among image signals with different frequency band components decomposed by the multiresoluting transform; a second transform step of decomposing the image signal in the lower frequency band component into image signals with different frequency band components using at least a first level dyadic wavelet transform; a signal intensity adjusting, step of reducing the signal intensity when an intensity of an image signal in a higher frequency band component of the extracted pixels which form a skin colored area among image signals in different frequency band components decomposed by the dyadic wavelet transform.  
      As an embodiment written in item 15, in the embodiment written in item 14, it is preferable that the image processing method further comprise: a color information converting step of converting the image signal with a lower frequency band component among the decomposed frequency band components into a brightness signal and a color-difference signal, wherein the second transform step applies at least a first level dyadic wavelet transform to the brightness signal and/or the color-difference signal.  
      As an embodiment written in item 16, in the embodiment written in item 15, it is preferable that the image processing method further comprises: a contrast controlling step of controlling a contrast of a image signal of the extracted pixels which form a skin colored area.  
      As an embodiment written in item 17, in the embodiment written in item 15 or item 16, it is preferable that the image processing method further comprise: a noise adding step of adding a noise signal to the image signal of the extracted pixels which form a skin colored area and whose image signal is reduced by the signal intensity regulating step.  
      To solve the above problem, an embodiment written in item 18 is an image processing apparatus including an obtaining an obtaining section of obtaining an image signal representing a color image, an extracting section of extracting pixels which form a skin colored area from pixels which are included in the obtained image signal; a signal intensity adjusting section of reducing a variation amount of a signal intensity of a pixel whose spatial frequency is in a predetermined range and whose signal intensity variation amount is a predetermined threshold or below among the extracted pixels.  
      Further, as an embodiment written in item 19, in the embodiment written in item 18, it is preferable that the extracting section extracts pixels which form a skin colored area by applying a dyadic wavelet transform to the obtained image signal and applying a simple area expansion to a lower frequency band component of the transformed image signal.  
      As an embodiment written in item 20, in the embodiment written in item 18 or item 19, it is preferable that the spatial frequency within the predetermined range is obtained by a higher frequency band component calculated by a multiresolution transform.  
      As an embodiment written in item 21, in the embodiment written in item 20, it is preferable that the signal intensity adjusting section determines the predetermined threshold based on a statistic value calculated by the higher frequency band component of an image signal intensity.  
      As an embodiment written in item 22, in the embodiment written in item 21, it is preferable that the statistic value is one of a standard deviation, an average, a median and a mode.  
      As an embodiment written in item 23, in the embodiment written in any one of items 18-22, it is preferable that when a pixel among the extracted pixels has a spatial frequency of 1.5-3.0 lines/mm and a signal intensity variation is in a range between 0-20% of a maximum signal intensity, the signal intensity adjusting section reduces a variation amount of a signal intensity of the pixel by 0.5-0.75 times.  
      As an embodiment written in item 24, in the embodiment written in any one of items 18-23, it is preferable that the image processing apparatus further comprise: a color information converting section of converting the obtained image signal into a brightness signal and a color-difference signal, wherein the signal intensity regulating section is applied to the brightness signal and/or the color-difference signal of the extracted pixels which form a skin colored area.  
      As an embodiment written in item 25, in the embodiment written in item 24, it is preferable that the image processing apparatus further comprise: a contrast controlling section of controlling a contrast of an image signal of the extracted pixels which form a skin colored area and whose image signal is reduced in the signal intensity adjusting section.  
      As an embodiment written in item 26, in the embodiment written in item 24 or item 25, it is preferable that the image processing apparatus further comprise: a noise adding section of adding a noise signal to the image signal of the extracted pixels which form a skin colored area and whose image signal is reduced by the signal intensity regulating section.  
      To solve the above problem, an embodiment written in item 27 is an image processing apparatus of item 18 further including a color information converting section of converting the obtained image signal into a brightness signal and a color-difference signal; and a wavelet transform section of decomposing the converted image signal into image signals with different frequency band components by applying at least a second level dyadic wavelet transform to the converted image signal, wherein when an image signal intensity in at least a higher frequency band component at level 2 of the extracted pixels among the decomposed frequency band components is predetermined threshold value or below, the image intensity adjusting section reduces the signal intensity.  
      In other words, the embodiment written in item 27 is an image processing apparatus including an obtaining section for obtaining an image signal representing a color image; an extracting section for extracting pixels which form a skin colored area from pixels which are included in the obtained image signal; a color information converting section for converting the obtained image signal into a brightness signal and a color-difference signal; a wavelet transform section for decomposing the converted image signal into image signals with different frequency band components by applying at least a second level dyadic wavelet transform to the converted image signal; and a signal intensity adjusting section for reducing a variation amount of a signal intensity of a pixel whose spatial frequency is in a predetermined range and whose signal intensity variation is a predetermined threshold or below among the transformed pixels.  
      As an embodiment written in item 28, in the embodiment written in items 27, it is preferable that the wavelet transform section applies the dyadic wavelet transform to the brightness signal and/or the color-difference signal.  
      As an embodiment written in item 29, in the embodiment written in item 28, it is preferable that the image processing apparatus further comprise a contrast controlling section of controlling a contrast of a image signal in the lower frequency band component of the extracted pixels which form a skin colored area among the image signals with different frequency band components decomposed by at least a first level of dyadic wavelet transform in the wavelet transform section.  
      As an embodiment written in item 30, in the embodiment written in item 28 or item 29, it is preferable that the image processing apparatus further comprise: a noise adding section of adding a noise signal to the image signal of the extracted pixels which form a skin colored area and whose image signal is reduced by the signal intensity regulating section.  
      To solve the above problem, an embodiment written in item 31 is an image processing apparatus of item 18 further comprise: a first transform section and a second transform section, wherein the first transform section decomposes the obtained image signal into image signals with different frequency band components by applying at least the first level multiresolution transform which is a apparatus of reducing an image size, the extracting section extracts pixels which form skin colored area from pixels included with the image signal in a lower frequency band component among image signals with different frequency band components decomposed by the multiresoluting transform, the second transform section decomposes the image signal with the lower frequency band component into image signals with different frequency band components using at least a first level dyadic wavelet transform, and when an intensity of an image signal with a higher frequency band component of the extracted pixels which form a skin colored area among image signals in different frequency band components decomposed by the dyadic wavelet transform, the signal intensity adjusting section reduces the signal intensity.  
      In other words, the embodiment written in item 31 is an image processing apparatus including: an obtaining section for of obtaining an image signal representing a color image; a first transform section for decomposing the obtained image signal into image signals with different frequency band components by applying at least the first level multiresolution transform which is a method of reducing an image size; an extracting section for extracting pixels which form a skin colored area from pixels included in the image signal in a lower frequency band component among image signals with different frequency band components decomposed by the multiresoluting transform; a second transform section for decomposing the image signal in the lower frequency band component into image signals with different frequency band components using at least a first level dyadic wavelet transform; a signal intensity adjusting section for reducing the signal intensity when an intensity of an image signal in a higher frequency band component of the extracted pixels which form a skin colored area among image signals in different frequency band components decomposed by the dyadic wavelet transform.  
      As an embodiment written in item 32, in the embodiment written in item 31, it is preferable that the image processing apparatus further comprise: a color information converting section of converting the image signal with a lower frequency band component among the decomposed frequency band components into a brightness signal and a color-difference signal, wherein the second transform section applies at least a first level dyadic wavelet transform to the brightness signal and/or the color-difference signal.  
      As an embodiment written in item 33, in the embodiment written in item 32, it is preferable that the image processing apparatus further comprise: a contrast controlling section of controlling a contrast of a image signal of the extracted pixels which form a skin colored area.  
      As an embodiment written in item 34, in the embodiment written in item 32 or item 33, it is preferable that the image processing apparatus further comprise: a noise adding section of adding a noise signal to the image signal of the extracted pixels which form a skin colored area and whose image signal is reduced by the signal intensity regulating section.  
      To solve the above problem, an embodiment written in item 35 is an image processing program for use in a computer executing an image processing, including an obtaining an obtaining step of obtaining an image signal representing a color image, an extracting step of extracting pixels which form a skin colored area from pixels which are included in the obtained image signal; a signal intensity adjusting step of reducing a variation amount of a signal intensity of a pixel whose spatial frequency is in a predetermined range and whose signal intensity variation amount is a predetermined threshold or below among the extracted pixels.  
      Further, as an embodiment written in item 36, in the embodiment written in item 35, it is preferable that the extracting step extracts pixels which form a skin colored area by applying a dyadic wavelet transform to the obtained image signal and applying a simple area expansion to a lower frequency band component of the transformed image signal.  
      As an embodiment written in item 37, in the embodiment written in item 35 or item 36, it is preferable that the spatial frequency within the predetermined range is obtained by a higher frequency band component calculated by a multiresolution transform.  
      As an embodiment written in item 38, in the embodiment written in item 37, it is preferable that the signal intensity adjusting step determines the predetermined threshold based on a statistic value calculated by the higher frequency band component of an image signal intensity.  
      As an embodiment written in item 39, in the embodiment written in item 38, it is preferable that the statistic value is one of a standard deviation, an average, a median and a mode.  
      As an embodiment written in item 40, in the embodiment written in any one of items 35-39, it is preferable that when a pixel among the extracted pixels has a spatial frequency of 1.5-3.0 lines/mm and a signal intensity variation is in a range between 0-20% of a maximum signal intensity, the signal intensity adjusting step reduces a variation amount of a signal intensity of the pixel by 0.5-0.75 times.  
      As an embodiment written in item 41, in the embodiment written in any one of items 35-40, it is preferable that the image processing program further comprise: a color information converting step of converting the obtained image signal into a brightness signal and a color-difference signal, wherein the signal intensity regulating step is applied to the brightness signal and/or the color-difference signal of the extracted pixels which form a skin colored area.  
      As an embodiment written in item 42, in the embodiment written in item 41, it is preferable that the image processing program further comprise: a contrast controlling step of controlling a contrast of an image signal of the extracted pixels which form a skin colored area and whose image signal is reduced in the signal intensity adjusting step.  
      As an embodiment written in item 43, in the embodiment written in item 41 or item 42, it is preferable that the image processing program further comprise: a noise adding step of adding a noise signal to the image signal of the extracted pixels which form a skin colored area and whose image signal is reduced by the signal intensity regulating step.  
      To solve the above problem, an embodiment written in item 44 is an image processing program of item 35 further including a color information converting step of converting the obtained image signal into a brightness signal and a color-difference signal; and a wavelet transform step of decomposing the converted image signal into image signals with different frequency band components by applying at least a second level dyadic wavelet transform to the converted image signal, wherein when an image signal intensity in at least a higher frequency band component at level 2 of the extracted pixels among the decomposed frequency band components is predetermined threshold value or below, the image intensity adjusting step reduces the signal intensity.  
      In other words, an embodiment written in item 44 is an image processing program including an obtaining step of obtaining an image signal representing a color image; an extracting step of extracting pixels which form a skin colored area from pixels which are included in the obtained image signal; a color information converting step of converting the obtained image signal into a brightness signal and a color-difference signal; a wavelet transform step of decomposing the converted image signal into image signals with different frequency band components by applying at least a second level dyadic wavelet transform to the converted image signal; and a signal intensity adjusting step of reducing a variation amount of a signal intensity of a pixel whose spatial frequency is in a predetermined range and whose signal intensity variation is a predetermined threshold or below among the transformed pixels.  
      As an embodiment written in item 45, in the embodiment written in items 44, it is preferable that the wavelet transform step applies the dyadic wavelet transform to the brightness signal and/or the color-difference signal.  
      As an embodiment written in item 46, in the embodiment written in item 45, it is preferable that the image processing program further comprise a contrast controlling step of controlling a contrast of a image signal in the lower frequency band component of the extracted pixels which form a skin colored area among the image signals with different frequency band components decomposed by at least a first level of dyadic wavelet transform in the wavelet transform step.  
      As an embodiment written in item 47, in the embodiment written in item 45 or item 46, it is preferable that the image processing program further comprise: a noise adding step of adding a noise signal to the image signal of the extracted pixels which form a skin colored area and whose image signal is reduced by the signal intensity regulating step.  
      To solve the above problem, an embodiment written in item 48 is an image processing program of item 35 further comprise: a first transform step and a second transform step, wherein the first transform step decomposes the obtained image signal into image signals with different frequency band components by applying at least the first level multiresolution transform which is a program of reducing an image size, the extracting step extracts pixels which form skin colored area from pixels included in the image signal with a lower frequency band component among image signals with different frequency band components decomposed by the multiresoluting transform, the second transform step decomposes the image signal in the lower frequency band component into image signals with different frequency band components using at least a first level dyadic wavelet transform, and the signal intensity adjusting step reduces the signal intensity when an intensity of an image signal with a higher frequency band component of the extracted pixels which form a skin colored area among image signals in different frequency band components decomposed by the dyadic wavelet transform.  
      In other words, the embodiment written in item 48 is an image processing program including: an obtaining step of obtaining an image signal representing a color image; a first transform step of decomposing the obtained image signal into image signals with different frequency band components by applying at least the first level multiresolution transform which is a method of reducing an image size; an extracting step of extracting pixels which form a skin colored area from pixels included in the image signal in a lower frequency band component among image signals with different frequency band components decomposed by the multiresoluting transform; a second transform step of decomposing the image signal in the lower frequency band component into image signals with different frequency band components using at least a first level dyadic wavelet transform; a signal intensity adjusting step of reducing the signal intensity when an intensity of an image signal in a higher frequency band component of the extracted pixels which form a skin colored area among image signals in different frequency band components decomposed by the dyadic wavelet transform.  
      As an embodiment written in item 49, in the embodiment written in item 48, it is preferable that the image processing program further comprise: a color information converting step of converting the image signal with a lower frequency band component among the decomposed frequency band components into a brightness signal and a color-difference signal, wherein the second transform step applies at least a first level dyadic wavelet transform to the brightness signal and/or the color-difference signal.  
      As an embodiment written in item 50, in the embodiment written in item 49, it is preferable that the image processing program further comprise: a contrast controlling step of controlling a contrast of a image signal of the extracted pixels which form a skin colored area.  
      As an embodiment written in item 51, in the embodiment written in item 49 or item 50, it is preferable that the image processing program further comprise: a noise adding step of adding a noise signal to the image signal of the extracted pixels which form a skin colored area and whose image signal is reduced by the signal intensity regulating step.  
      Herein, terms written in items will be additionally described.  
      The “skin-colored area” means an area of the skin of the face or neck of a person included in an image signal. Hereupon, as a method to extract the skin-colored area, an image area obtained from the image signal with the lower frequency band components obtained when the dyadic wavelet transform is applied on the image signal, and from the pixels, when a simple range extension is applied, may also be defined as a skin-colored area, or centering around a most suitable eye of the object included in the image signal, the logarithmic polar coordinate transform is applied, and from the image signal after the processing, when a template matching or simple range extension is applied, a range judged as a most suitable face may also be defined as a skin-colored area, or the publicly-known extraction method may also be used.  
      “Spatial frequency” means a spatial frequency when the image signal is outputted on a printing paper, hard copy, and display device, and “variation amount of the signal intensity” means the difference between a signal intensity of a certain pixel, and the signal intensity of the pixel which is regulated by the spatial frequency.  
      A phrase “the image signal is converted into the brightness signal and the color difference signal” means, for example, a intensity signal for 3 colors of RGB representing the image signal of an object of the image processing, is converted into a intensity signal for the publicly known YIQ base or YUV base, or converted into a intensity signal for XYZ base of the CIE 1931 color specification system or L*a*b base or L*u*V* base which is advised by CIE 1976 according to a standard such as sRGB (standard RGB) or NTSC (National TV Standards Committee). Further, it may also means a conversion such that an average of RGB values of the signal is defined to the brightness signal, and 2-axis perpendicular to the brightness signal is defined to a color difference signal.  
      “Image size” means the number of pixels when the image including a color photographic film is photo-electrically read by a CCD sensor and converted into the image signal.  
      Further, “multiresolution transform” is remarked a general name of methods represented by the wavelet transform, perfect reconstruction filter banks, Laplacian pyramid. “multiresolution transform” is a transform for obtaining a multiresolution signal having a plurality of signal components with different frequency bands from an input signal, such that the input signal is decomposed to the lower frequency band component signal and the higher frequency band component signal by 1-time transform operation and the same transform operations are applied to the obtained lower frequency band component signal. When the inverse multiresolution transform is applied to the obtained multiresolution signal without any processing as it is, the original signal is re-structured.  
      Herein, as a representative example of the multiresolution transform, an outline of the wavelet transform will be described. The wavelet transform is a transform for decomposing an input signal f(x) into the total sum of wavelet functions shown in the following expression (3), such that the wavelet transform function (f, φ a,b ) for the input signal f(x) is found as in the following expression (2) by using the wavelet function (the following expression (1)) which is oscillated in the finite range as shown in  FIG. 1 .  
             (Math-1)                               Ψ     a   ,   b       ⁡     (   x   )       =     Ψ   (       x   -   b     a     )             (   1   )               (Math-2)                             〈     f   ,     Ψ     a   ,   b         〉     ≡       1   a     ⁢           ⁢     ∫         f   ⁡     (   x   )       ·     Ψ   (       x   -   b     a     )       ⁢     ⅆ   x                   (   2   )               (Math-3)                               f   ⁡     (   x   )       =       ∑     a   ,   b       ⁢       〈     f   ,     Ψ   a     ⁢     ,   b       〉     ·     Ψ   a           ⁢     ,   b     ⁢     (   x   )             (   3   )             
 
      In the above expressions (1)-(3), a expresses a scale of the wavelet function and b shows a position of the wavelet function. As illustrated in  FIG. 1 , the larger a value of the scale a is, the smaller the frequency of the wavelet function φ a,b  (X) is. According to a value of a position b, a position at which the the wavelet function φ a,b  (x) is oscillated, is shifted. Accordingly, the above expression (3) means that the input signal f(x) is decomposed to the total sum of the wavelet function φ a,b  (x) having a various scales and positions.  
      As “the multiresolution transform which is a method of reducing the image size” among above-described types of the wavelet transforms, the orthogonal wavelet transform, biorthogonal wavelet transform are well known. The outline of the calculation of the orthogonal and biorthogonal wavelet transform will be described below.  
             (Math-4)                               ψ     i   ,   j       ⁡     (   x   )       =       2     -   i       ⁢     ψ   (       x   -     j   ·     2   i           2   i       )               (   4   )             
 
      Where, i is a natural number.  
      When comparing the expression (4) to the expression (1), it can be seen that a value of the scale a is discretely defined by i-th power of 2, and further, the minimum movement unit of the position b is discretely defined by 2 i  in the orthogonal wavelet transform and the biorthogonal wavelet transform. This value of i is called a level.  
      When the level i is limited up to a finite upper limit N, the input signal f(x) is expressed as the following expressions (5)-(7).  
             (Math-5)                             f   ⁡     (   x   )       ≡     S   0             (   5   )                       ⁢     =         ∑   j     ⁢       (       S   0     ,     Ψ     1   ,   j         )     ·       Ψ     1   ,   j       ⁡     (   x   )           +       ∑   j     ⁢       (       S   0     ,     ϕ     1   ,   j         )     ·       ϕ     1   ,   j       ⁡     (   x   )                                               ⁢     ≡         ∑   j     ⁢         W   1     ⁡     (   j   )       ·       Ψ     1   ,   j       ⁡     (   x   )           +       ∑   j     ⁢         S   1     ⁡     (   j   )       ·       ϕ     1   ,   j       ⁡     (   x   )                                       (Math-6)                             S     i   -   1       =         ∑   j     ⁢       〈       S     i   -     1   ,         ⁢     Ψ       i   ,     ⁢   j         〉     ·       Ψ     i   ,   j       ⁡     (   x   )           +       ∑   j     ⁢       〈       S       i   -   1     ,       ⁢     ϕ     i   ,   j         〉     ·       ϕ     i   ,   j       ⁡     (   x   )                     (   6   )                       ⁢     ≡         ∑   j     ⁢         W   i     ⁡     (   j   )       ·       Ψ     i   ,   j       ⁡     (   x   )           +       ∑   j     ⁢         S   i     ⁡     (   j   )       ·       ϕ     i   ,   j       ⁡     (   x   )                                       (Math-7)                             f   ⁡     (   x   )       ≡         S   0     ⁢           ⁢     =     ⁢     ≡     ⁢           ⁢       ∑     i   =   1     N     ⁢       ∑   j     ⁢         W   i     ⁡     (   j   )       ·       ψ     i   ,   j       ⁡     (   x   )               +       ∑   j     ⁢         S   N     ⁡     (   j   )       ·       ϕ     i   ,   j       ⁡     (   x   )                     (   7   )             
 
      In the second term of the expression (5), the lower frequency band components of residuals which can not be expressed by the total sum of the wavelet function Φ l,j (x) of level 1, are expressed by the total sum of the scaling function Φ l,j (x) of level 1. An appropriate scaling function is used corresponding to the wavelet function. By applying the wavelet transform of level 1 shown in the expression (5), the input signal f(x)=S 0  is signal-decomposed to the higher frequency band components W 1  and the lower frequency band components S 1  of level 1.  
      Because the minimum shift unit of the wavelet function Φ i,j (x) is 2 i , the signal amounts of the higher frequency band components W 1  and the lower frequency band components S 1  are respectively {fraction (1/2)} to the signal amount of the input signal S 0 , and the total sum of the signal amounts of the higher frequency band components W 1  and the lower frequency band components S 1  equals to the signal amount of the input signal S 0 . The lower frequency band components S 1  of level 1 are decomposed to the higher frequency band components W 2  and the lower frequency band components S 2  by the expression (6), and hereinafter, by repeating the transform up to the level N in the same manner, the input signal S 0 , as shown in expression (7), is decomposed to the total sum of the higher frequency band components of level 1−N, and the sum of the lower frequency band components of level N.  
      Herein, it is well known that the wavelet transform shown in the expression (6), can be calculated by the filter processing as shown in  FIG. 2 . In  FIG. 2 , an LPF shows a lower pass filter, and a HPF shows a high pass filter. Filer coefficients of the low pass filter LPF and the high pass filter are appropriately determined corresponding to the wavelet function. In  FIG. 2 , “2θ” shows the down sampling which thins out the signals every other.  
      As shown in  FIG. 2 , when the input signal S n-1  is processed by the low pass filter LPF and the high pass filter HPF, and the signals are thinned out every other, the input signal S n-1  can be decomposed to the higher frequency band components W n  and the lower frequency band components S n .  
      The first level wavelet transform in the two dimensional signal as the image signal, is calculated by the filter processing as shown in  FIG. 3 . In  FIG. 3 , LPFx, HPFx and 2↓x show the processing in the x-direction, and LPFy, HPFy and 2↓y show the processing in the y-direction. Initially, the input signal S n-1  is filter processed by the low pass filter LPFx and high pass filter HPFx, in the x-direction, and the down-sampling is applied in the x-direction. Hereby, the input signal S n-1  is decomposed to the lower frequency band components SX n , and the higher frequency band components WX n . On each of the lower frequency band components SX n  and the higher frequency band components WX n , the filter processing by the low pass filter LPFy and high pass filter HPFy in the y-direction is applied, and the down-sampling is applied in the y-direction.  
      By this first level wavelet transform, the lower frequency band components S n-1  is decomposed to 3 higher frequency band components Wh n , Wv n , Wd n , and 1 lower frequency band components S n . Because each of signal amounts of Wh n , Wv n , Wd n  and S n , which are generated by the decomposition is {fraction (1/2)} in the length and width, as compared to S n-1  before the decomposition, the total sum of the signal amounts of 4 components after the decomposition is equal to the signal of s n-1  before decomposition.  
      A process in which the input signal S 0  is signal-decomposed by the third level wavelet transform, is typically shown in  FIG. 4 . It can be seen that, as the number of level is larger, the image signal is thinned out by the down-sampling, and the decomposed image becomes small as shown in  FIG. 4 .  
      Further, as shown in  FIG. 5 , it is well known that the signal S n-1  before the decomposition can be perfectly re-structured by applying the inverse wavelet transform calculated by the filter processing to Wh n , Wv n , Wd n  and S n  generated by the decomposition. In  FIG. 5 , LPF′ shows the low pass filter for the inverse transform, and HPF′ shows the high pass filter for the inverse transform. Further, “2↑” shows the up-sampling processing by which zero is inserted every other into the signal. Further, the LPF′x, HPF′x, 2↑x show a flow of the processing in the x-direction, and the LPF′y, HPF′y, 2↑y show a flow of the processing in the y-direction.  
      As shown in  FIG. 5 , SX n  is obtained by adding the signal obtained by processing the up-sampling processing and filter-processing by the low pass filter LPF′y in the y-direction on S n , and the signal obtained by applying the up-sampling processing and the filter-processing by the high pass filter HPF′y in the y-direction on Wh n . In the same manner as this, WX n  is generated from Wv n  and Wd n .  
      Further, the signal S n-1  before the decomposition can be re-structured by adding the signal obtained by applying the up-sampling processing and the filter processing by the low pass filter LPF′x in the x-direction on SX n , and the signal obtained by applying the up-sampling processing and the filter processing by the high pass filter HPF′x in the x-direction on WX n .  
      When the orthogonal wavelet transform has been applied, as a filter used for the inverse wavelet transform, a filter with the same coefficient as the coefficient used in the orthogonal wavelet transform is used. When the biorthogonal wavelet transform has been applied, a filter with the different coefficients from the coefficients used in the biorthogonal wavelet transform is used at the time of the inverse transform.  
      Next, an outline of a dyadic wavelet transform will be described. The wavelet function used in the dyadic wavelet transform is defined as the following expression (8).  
             (Math-8)                               ψ     i   ,   j       ⁡     (   x   )       =       2     -   i       ⁢     ψ   ⁡     (       x   -   j       2   i       )                 (   8   )             
 
      Where, i is a natural number.  
      As described above, the minimum movement unit of the position at level i is discretely defined in the wavelet functions including the orthogonal wavelet transform and biorthogonal wavelet transform. However, in the dyadic wavelet transform, the minimum shift unit of the position is constant regardless of level i. By this difference, the dyadic wavelet transform has the following features.  
      As the first feature, each of signal amounts of the higher frequency band components W i  and the lower frequency band components S i , generated by the first level dyadic wavelet transform shown in the following expression (9) is the same as the signal S i-1  before the transform.  
             (Math-9)                             S     i   -   1       =         ∑   j     ⁢       〈       S       i   -   1     ,       ⁢     ψ     i   ,   j         〉     ·       ψ     i   ,   j       ⁡     (   x   )           +       ∑   j     ⁢       〈       S       i   -   1     ,       ⁢     ϕ     i   ,   j         〉     ·       ϕ     i   ,   j       ⁡     (   x   )                     (   9   )                       ⁢     ≡         ∑   j     ⁢         W   i     ⁡     (   j   )       ·       ψ     i   ,   j       ⁡     (   x   )           +       ∑   j     ⁢         S   i     ⁡     (   j   )       ·       ϕ     i   ,   j       ⁡     (   x   )                                     
 
      Like this, the dyadic wavelet transform is different from the biorthogonal and orthogonal wavelet transform, and the image size after transform is not decreased from the original image.  
      As the second feature, the following relational expression (10) is satisfied for the scaling function φ i,j (x) and the wavelet function φ i,j (x).  
             (Math-10)                               ψ     i   ,   j       ⁡     (   x   )       =       δ     δ   ⁢           ⁢   x       ⁢       ϕ     i   ,   j       ⁡     (   x   )                 (   10   )             
 
      Accordingly, the higher frequency band components W i  generated by the dyadic wavelet transform is expressed by the first order differentiation (gradient) of the lower frequency band components S i .  
      As the third feature, when Wi-γi (hereinafter, called the corrected higher frequency band components) is calculated by multiplying the higher frequency band components to the coefficient γ i  (hereinafter, called correction coefficient) determined corresponding to the level i of the wavelet transform, corresponding to the singularity of a signal change of the input signal, the relationship between the levels of the signal density of the corrected higher frequency band components W i -γ i  follows predetermined formula.  
      In FIGS.  6 ( a )- 6 ( e ), a waveform of the input signal So and a waveform of the corrected higher frequency band components of each level obtained by the wavelet transform are shown.  FIG. 6 ( a ) shows the input signal S 0 ,  FIG. 6 ( b ) shows the corrected higher frequency band components W 1 -γ 1  obtained by the level-i dyadic wavelet transform,  FIG. 6 ( c ) shows the corrected higher frequency band components W 2 -γ 2  obtained by the level-2 dyadic wavelet transform,  FIG. 6 ( d ) shows the corrected higher frequency band components W 3 -γ 3  obtained by the level-3 dyadic wavelet transform, and  FIG. 6 ( e ) shows the corrected higher frequency band components W 4 -γ 4  obtained by the level-4 dyadic wavelet transform.  
      When a change of the signal intensity in each level is observed, the corrected higher frequency band components W i -γ i  is corresponding to a smooth (differentiable) signal variation shown in “1” or “4” in  FIG. 6 ( a ). The signal intensity of the corrected higher frequency band components W i -γ i , is increased as a value of i (level) is increased, as shown in FIGS.  6 ( b ) through  6 ( e ).  
      In the input signal S 0 , the corrected higher frequency band components W i -γ i  is corresponding to a step-like signal change shown in “2” of  FIG. 6 ( a ). The signal intensity of the corrected higher frequency band components W i -γ i  becomes constant irrespective of a value of i. In the input signal S 0 , the corrected higher frequency band components W i -γ i  is corresponding to a δ-functional signal change shown in “3”. The signal intensity of the corrected higher frequency band components W i -γ i  is reduced as a number of level i is increased, as shown in FIGS.  6 ( b ) through  6 ( e ).  
      As the fourth feature, a method of the first level dyadic wavelet transform in the second dimensional signal as the image signal is different from the above-described orthogonal wavelet transform or the biorthonal wavelet transform, and applied by a method shown in  FIG. 7 .  
      As shown in  FIG. 7 , the lower frequency band components S n  is obtained by processing the input signal S n-1  by the low pass filter LPFx in x direction and the low pass filter LPFy in y direction using the first level dyadic wavelet transform. Further, by processing the input signal S n-1  is processed using the high pass filter in x direction, the higher frequency band components W xn  are obtained. Further, by processing the input signal S n-1  using the high pass filter HPFy in y direction, the another higher frequency band components W yn  are obtained.  
      In this manner, the input signal S n-1  is decomposed to 2 higher frequency band components W xn , W yn , and one lower frequency components S n  by the first level dyadic wavelet transform. The 2 higher frequency band components W xn , W yn  correspond to x component and y component of the change vector V n  in 2 dimension of the lower frequency band components S n . The magnitude M n  and the deflection angle An of the change vector V n  are given in the following expressions (11) and (12).  
             (Math-11)                             M   n     =         W   ⁢           ⁢     x   n   2       +     W   ⁢           ⁢     y   n   2                   (   11   )             
 
 (Math-12) 
 
 An=arg ( Wx   n   +iWy   n )  (12) 
 
      Further, when the inverse dyadic wavelet transform shown in  FIG. 8  is applied to the 2 higher frequency band components Wx n  and Wy n  and one lower frequency band components S n  obtained by the dyadic wavelet transform, the signal S n-1  before the transform can be re-structured. That is, the signal S n-1  before the dyadic wavelet transform can be obtained by adding a signal obtained by processing S n  by the low pass filter LPFx in x-direction and the-low pass filter LPFy in y-direction, a signal obtained by processing Wx n  by the high pass filter HPF′x in x-direction and the low pass filter LPF′y in y-direction, and a signal obtained by processing Wy n  by the low pass filter LPF′x in x-direction and the high pass filter HPF′y in y-direction.  
      Next, according to the block diagram in  FIG. 9 , a method of obtaining the output signal S 0 ′ from the input signal S 0 , such that after the n-th level (n is a natural number) dyadic wavelet transform is applied to the input signal S 0 , and any one of image processing (in  FIG. 9 , “editing” is written) is applied to the obtained higher frequency components, lower frequency components, the n-th level inverse dyadic wavelet transform is applied, and the output signal S 0 ′ is obtained.  
      By the level-1 dyadic wavelet transform to the input signal S 0 , the input signal S 0  is decomposed to the signal of 2 higher frequency band components Wx 1 , Wy 1  and the lower frequency components S 1 . The lower frequency band components S 1  obtained by the level-1 dyadic wavelet transform is further decomposed to 2 higher frequency band components Wx 2 , Wy 2  and the lower frequency band components S 2  using the wavelet transform of level 2. When such a decomposition operation is repeated to the level n, the input signal S 0  is decomposed to a plurality of the higher frequency band components Wx 1 , Wx 2 , . . . , Wx n , Wy 1 , Wy 2 , . . . , Wy n  and one lower frequency band components S n .  
      After the image processing (editing) is applied to the higher frequency ban components Wx 1 , Wx 2 , . . . , Wx n , Wy 1 , Wy 2 , . . . , Wy n  and one lower frequency band components Sn, which are obtained in this manner, the higher frequency band components Wx 1 ′, Wx 2 ′, . . . , WX n ′, Wy 1 ′, Wy 2 ′, . . . , Wy n ′ and one lower frequency band components S n ′ are obtained.  
      Then, the inverse dyadic wavelet transform is applied to these higher frequency band components Wx 1 ′, Wx 2 ′, . . . , Wx n ′, Wy 1 ′, Wy 2 ′, . . . , Wy n ′ and one lower frequency band components S n ′. That is, the lower frequency components S n-1 ′ of the level 1 that has been processed by the image processing (editing), is structured from the 2 higher frequency band components Wxn′, Wyn′ in the level n and the lower frequency band components Sn′ in the level n, that have been processed by the image processing (editing). After such an operation is repeated, the lower frequency components S 1 ′ of the level 1 after the image processing is structured from the 2 higher frequency band components Wx 2 ′, Wy 2 ′ in the level 2 after the image processing, and the lower frequency band components S2′. From this the lower frequency band components S 1 ′ and 2 higher frequency band components Wx 1 ′, Wy 1 ′ in the level 1 after the image processing, the image signal S 0 ′ is structured.  
      Hereupon, a filter coefficient of each filter used in  FIG. 9  is adequately determined corresponding to the dyadic wavelet transform. Further, in the dyadic wavelet transform, the filter coefficient used for each level is different. In the filter coefficients used in the level n, coefficients in which zeros of 2 n-1  pieces are inserted between each of coefficients of filters of level 1 are used.  
      Further, in  FIG. 9 , a flow of the image processing (editing) applied to the image of the higher frequency band components obtained by the dyadic wavelet transform and the lower frequency band components of the final level is shown, however, the image processing may be applied to the image of the lower frequency band components synthesized after the inverse dyadic wavelet transform. Furthermore, the image processing may also be applied to the image of the lower frequency band components on the halfway of the dyadic wavelet transform.  
      As “the multiresolution transform which is a method of reducing the image size” in embodiments written in items 14, 31, 48, it is preferable that the biorthogonal and orthogonal wavelet transform is used. Further, in the biorthogonal and orthogonal wavelet transform, because the image size of the lower frequency band component image after the transform can be made ¼ of the original image, it is preferable also from the view point of the processing load.  
      A phrase of “to reduce the signal intensity” in embodiments written in items 10, 14, 27, 31, 44, 48 means to process the absolute value of the signal intensity of pixels so that it is decreased. Further, the pixels whose “signal intensity variation is a predetermined threshold or below” may be selected, for example, by a threshold value determined according to a standard deviation σ of the signal intensity of the higher frequency band components, a threshold value determined by using the mean value, median, or mode of the signal intensity, or a threshold value obtained based on a comparison of the image signal of pixels corresponding to the image signal of the corrected higher frequency band components of the P-th level (P is a natural number) obtained by the dyadic wavelet transform, to the image signal of the corrected higher frequency band components of the (P+1)th level or the (P−1)th level.  
      Further, “the first” level and “the second” level written in each of items mean the number of times at which the multiresolution transform or the dyadic wavelet transform is applied.  
      According to embodiments written in items 1, 18, 35, the image signal expressing the color image is obtained, pixels forming the skin-colored area are extracted from pixels included in the image signal, and the image processing to reduce the variation amount of the signal intensity is applied to the pixels which have the spatial frequency within the specific range and have the variation amount of the signal intensity with a specific threshold value or below, among the extracted pixels. Because the sharpness feeling of the whole image is not deteriorated, further, the outline structure of the face of the person is not changed, the spotting or puckering can be removed without a sense of incompatibility.  
      Hereupon, as a threshold value of the range of the spatial frequency and a variation amount of the signal intensity, it is preferable that the threshold value is appropriately changed corresponding to the age, sex, gradation of the skin, an amount of spotting or puckering, of the person recorded in the image signal. Further, also for an adjusting amount to reduce the variation amount of the signal intensity, in the same manner, it is preferable that it is appropriately changed corresponding to the age, sex, gradation of the skin, an amount of spotting or puckering, of the person recorded in the image signal. Further, because the human puckering is many in the lateral direction (a parallel direction to the eyes and mouth), the puckering in the lateral direction may also be adjusted so that it can be intensively removed.  
      According to embodiments written in items 7, 24, 41, by converting the image signal into the brightness signal and the color difference signal, and applying the image processing reducing the variation amount of the signal intensity to the brightness signal and/or color difference signal of pixels forming the skin-colored area, the sharpness feeling of the whole image is not deteriorated, and the outline structure of the face of the person is not changed. Therefore, the spotting or puckering can be removed without a sense of the incompatibility. Hereupon, it is preferable that the image processing is applied to the brightness signal for the puckering or chapping of the skin of a person, and it is preferable that the image processing is applied to the color difference signal for the spotting or freckle.  
      According to embodiments written in items 8, 25, 42, by adjusting the contrast for the image signal of pixels forming the skin-colored area after the variation amount of the signal intensity in the signal intensity adjusting process is reduced, the image without a sense of incompatibility can be presented even when the result of the image processing of the present invention is perceived as if the contrast is changed.  
      According to embodiments written in items 9, 26, 43, by adding the noise signal to the image signal of pixels forming the skin-colored area after the variation amount of the signal intensity is reduced, the image without a sense of incompatibility can be presented even when, a result of the image processing of the present invention is perceived as if the contrast is changed.  
      According to embodiments written in items 10, 27, 44, the image signal expressing the color image is obtained, pixels forming the skin-colored area is extracted from the pixels included in the image signal, and the image signal is converted into the brightness signal and the color difference signal, the image processing for reducing the signal intensity is applied in case that the signal intensity of the image signal of at least level-2 higher frequency band components of pixels forming the skin-colored area is not larger than a specific threshold value among the frequency band components obtained by applying at least the second level dyadic wavelet transform is applied to the converted image signal and decomposing the image signal into the image signal with different frequency band components. Therefore, because the sharpness feeling of the whole image is not deteriorated and the outline structure of the face of the person is not changed, the spotting and the puckering can be removed without any sense of incompatibility.  
      Hereupon, it is preferable that the threshold value is appropriately changed corresponding to the age, sex, gradation of the skin, and an amount of the spotting or puckering, of the person recorded in the image signal as the threshold value of the signal intensity. Further, also for an adjusting amount to reduce the signal intensity, in the same manner, it is preferable that it is appropriately changed corresponding to the age, sex, gradation of the skin, and an amount of the spotting or puckering, of the person recorded in the image signal. Further, because the human puckering is many in the lateral direction (a parallel direction to the eyes and mouth), the puckering in the lateral direction may also be adjusted so that it can be intensively removed.  
      According to embodiments written in items 11, 28, 45, the dyadic wavelet transform is applied to the brightness signal and/or color difference signal to decompose the brightness signal and/or color difference signal into the image signal with the various frequency band components, and when the signal intensity of the image signal of at least level-2 higher frequency band components of pixels forming the skin-colored area is not larger than a specific threshold value, the image processing to reduce the signal intensity is applied. Therefore the sharpness feeling of the whole image is not deteriorated, further, the outline structure of the face of the person is not changed, the spotting and the puckering can be removed without any sense of incompatibility. Hereupon, for the puckering or freckle of the skin of a person, it is preferable that the image processing is applied to the brightness signal, and for the spotting or chapping, the image processing is applied on color difference signal.  
      According to embodiments written in items 12,29,46, the contrast is adjusted for the image signal of the lower frequency band components of pixels forming the skin-colored area among the image signals with the various frequency band components decomposed by the first level dyadic wavelet transform. Therefore, the image not having any sense of incompatibility can be presented even when a result of the image processing of the present invention is perceived as if the contrast is changed.  
      According to embodiments written in items 13, 30, 47, the noise signal is added to the image signal of pixels forming the skin-colored area after the variation amount of the signal intensity is reduced. The image not having any sense of incompatibility can be presented even when a result of the image processing of the present invention is perceived as if the contrast is changed.  
      According to embodiments written in items 14, 31, 49, the image signal expressing the color image is obtained, pixels forming the skin-colored area are extracted from the pixels included in the image signal of the lower frequency band components, among the various frequency band components obtained by applying the multiresolution transform at least at the first level which is a method of reducing an image size to the obtained image signal, and the image processing to reduce the signal intensity is applied when the signal intensity of the image signal of the higher frequency band components of pixels forming the skin-colored area is not larger than a specific threshold value among image signals with the various frequency band components obtained by applying at least first level dyadic wavelet transform to the image signal with the lower frequency band components. Therefore, the productivity of the processing can be increased because the image processing is applied after the image size is decreased, and further, the outline structure of the face of the person is not changed, the spotting and the puckering can be removed without any sense of incompatibility because the sharpness feeling of the whole image is not deteriorated.  
      Hereupon, it is preferable that the threshold value is appropriately changed corresponding to the age, sex, gradation of the skin, and an amount of the spotting or puckering, of the person recorded in the image signal as the threshold value of the signal intensity. Further, also for an adjusting amount to reduce the signal intensity, in the same manner, it is preferable that it is appropriately changed corresponding to the age, sex, gradation of the skin, and an amount of the spotting or puckering, of the person recorded in the image signal. Further, because the human puckering is many in the lateral direction (a parallel direction to the eyes and mouth), the puckering in the lateral direction may also be adjusted so that it can be intensively removed.  
      According to embodiments written in items 15, 32, 49, the image signal of the lower frequency band components obtained by the multiresolution transform is converted into the brightness signal and the color difference signal, and on the brightness signal and/or color difference signal, in the frequency band components obtained when at least the first level dyadic wavelet transform is applied, in the case where the signal intensity of the image signal of the higher frequency band components is not larger than a specific threshold value, when the image processing to reduce the signal intensity is applied, because the sharpness feeling of the whole image is not deteriorated, further, the outline structure of the face of the person is not changed, the spotting and the puckering can be removed without any sense of incompatibility. Hereupon, for the puckering or freckle of the skin of a person, it is preferable that the image processing is applied to the brightness signal, and for the spotting or chapping, the image processing is applied to the color difference signal.  
      According to embodiments written in items 16, 33, 50, the contrast is adjusted for the image signal of pixels forming the skin-colored area. Even when a result of the image processing of the present invention is perceived as if the contrast is changed, the image not having any sense of incompatibility can be presented.  
      According to embodiments written in items 17, 34, 51, the noise signal is added to the image signal of pixels forming the skin-colored area after the variation amount of the signal intensity is reduced. Even when a result of the image processing of the present invention is perceived as if the contrast is changed, the image not having any sense of incompatibility can be presented.  
     DETAILED DESCRIPTION OF THE INVENTION  
      Referring to the drawings, the preferred embodiments to carry out the present invention will be detailed below. The external structure of the image processing apparatus  1 :  
      Initially, referring to  FIG. 10 , the external structure of an image processing apparatus will be described.  
      In the image processing apparatus  1 , as shown in  FIG. 10 , a magazine loading section  3  is provided. On one side surface of a casing  2 , the magazine loading section  3  for loading a photosensitive material is provided. Inside of the casing  2 , an exposure processing section  4  for exposing the photosensitive material, and a print making section  5  for development processing the exposed photosensitive material and drying, and for making a print, are provided. On the other side surface of the casing  2 , a tray  6  for discharging the print made in the print making section  5  is provided.  
      Further, on an upper portion of the casing  2 , a CRT (Cathode Ray Tube)  8  as a display device, a film scanner section  9  which is a device for reading-in a transmission document, a reflection document input device  10 , and an operation section  11  are provided. Further, in the casing  2 , an image reading-in section  14  which can read the image information recorded in each kind of digital recording medium, and an image writing section  15  which can write an image signal in each kind of digital recording medium, are provided. Further, inside the casing  2 , a control section  7  for overall-controlling each section of them is provided.  
      In the image reading-in section  14 , an adapter for a PC card  14   a , an adapter for a floppy disk (registered trade mark, same hereinafter)  14   b  are provided, and a PC card  13   a  or a floppy disk  13   b  can insert into them. For example, the PC card  13   a  has a memory in which the information of a plurality of frame images image-picked up by a digital camera is recorded. In the floppy disk  13   b , for example, the information of a plurality of frame images image-picked up by the digital camera is recorded.  
      In the image writing section  15 , an adapter  15   a  for the floppy disk, an adapter for MO (Magneto-Optical)  15   b , an adapter  15   c  for an optical disk are provided, and respectively, a floppy disk  16   a , MO  16   b , and optical disk  16   c  can insert into them. As the optical disk  16   c , there are a CD-R (CD Recordable), CD-RW (CD Rewritable), DVD±R, DVD±RW, and Blu-ray Disc.  
      Hereupon, in  FIG. 10 , a structure in which the operation section  11 , CRT  8 , film scanner section  9 , reflection document input device  10 , image reading section  14  are integrally provided in the casing  2 , is applied, however, any one of them may also be provided as a separated body.  
      Hereupon, in the image processing apparatus  1  shown in  FIG. 10 , an apparatus which exposes the photosensitive material, develops, and makes a print, is illustrated, however, a print making method is not limited to this, for example, a system such as an inkjet system, electronic photographing system, heat sensitive system, sublimation system, may also be used.  
      Inside Structure of the Image Processing Apparatus  1 :  
      Next, referring to  FIG. 11 , an inside structure of the image processing apparatus  1  will be described. The image processing apparatus  1  is provided with, as shown in  FIG. 11 , the control section  7 , exposure processing section  4 , print making section  5 , film scanner section  9 , reflection document input device  10 , image reading section  14 , communicating means (input)  32 , image writing section  15 , data accumulation means  71 , operating section  11 , CRT  8 , communicating means (output)  33 .  
      The control section  7  is structured by a microcomputer, and each kind of control programs such as an image processing program stored in a ROM (Read Only Memory), (not shown), and by a cooperation with a CPU (Central Processing Unit), (not shown), an operation of each section structuring the image processing apparatus  1 , is overall controlled.  
      The control section  7  has the image processing section  70 , and according to the input signal (instruction information) from the operating section  11 , to the image data obtained by the film scanner section  9  or the reflection document input device  10 , image data read from the image reading section  14 , and image data inputted from the external equipment through the communication means  32 , forms the image information for exposure, and outputs to the exposure processing section  4 . Further, the image processing section  70  applies the transform processing corresponding to the output mode on the image processed image data, and outputs it. As the output destination of the image processing section  70 , there are the CRT  8 , image writing section  15 , communication means (output)  33 , and data accumulation means  71 .  
      The exposure processing section  4  exposes the image on the photosensitive material, and outputs it to the print making section  5 . The making section  5  developing-processes the exposed photosensitive material, dries, and makes the prints P 1 , P 2 , P 3 . The print P 1  is a print such as service size, hi-vision size, panorama size, the print P 2  is a print of A4-size print, and the print P 3  is a print of name card size.  
      The film scanner section  9  reads-in the information of the frame image from each kind of transmission-type documents such as the developed negative film or reversal film, image picked-up by an analog camera. The reflection document input device  10  reads-in each kind of images formed in the print P (photographic print).  
      The image reading section  14  has an image transfer means  30 , and reads-out the frame image information recorded in the PC card  13   a  or floppy disk  13   b , and transfers it to the control section  7 . The image transfer means  30  has an adapter  14   a  for the PC card, and an adapter  14   b  for the floppy disk. The image reading section  14  reads out the information of the frame image recorded in the PC card  13   a  inserted into the adapter  14   a  for the PC card, or the floppy disk  13   b  inserted into an adapter  14   b  for the floppy disk, and transfers it to the control section  7  by using the image transfer means  30 . As the adapter  14   a  for the PC card, for example, a PC card reader or a PC card slot is used.  
      The communication means (input)  32  receives an image signal expressing the picked-up image or a print command signal from another computer in the facility in which the image processing apparatus  1  is installed, or a remote computer through the internet.  
      The image writing section  15  is provided with an adapter  15   a  for the floppy disk, adapter  15   b  for MO, and adapter  15   c  for the optical disk, as the image conveying section  31 . The image writing section  15 , according to a writing signal inputted from the control section  7 , writes each kind of data in the floppy disk  16   a  inserted into the adapter  15   a  for the floppy disk, MO  16   b  inserted into the adapter  15   b  for MO, and the optical disk  16   c  inserted into the adapter  15   c  for the optical disk.  
      The data accumulate means  71  accumulates the image information and an order information (the information in which, from which frame image, which number of prints are made, or the information of a print size) corresponding to it.  
      The operating section  11  has an information input means  12 . The information input means  12  is structured by, for example, a touch panel, and outputs the pressing signal of the information input means  12  as the input signal to the control section  7 . Hereupon, the operating section  11  may also be structured by providing a key-board or mouse. The CRT  8  displays the image information according to a display control signal inputted from the control section  7 .  
      The communication means (output)  33  sends the image signal expressing the photographed image, and the order information accompanied to it, to another computer in the facility in which the image processing device  1  is installed, or a remote computer through internet. Structure of the image processing section  70 :  
      Next, referring to  FIG. 12 , a structure of the image processing section  70  will be described.  
      The image processing section  70  is, as shown in  FIG. 12 , provided with a film scan data processing section  701 , reflection document scan data processing section  702 , image data format decode processing section  703 , image data adjustment processing section  704 , CRT proper processing section  705 , printer proper processing section  706 , printer proper processing section  707 , and image data format making section  708 .  
      The film scan data processing section  701  applies a correction operation proper to the film scanner section  9 , a negative positive reversal in the case of a negative document, a gray balance adjustment, a contrast adjustment, on the image information inputted from the film scanner section  9 , and outputs it to the image adjustment processing section  704 . Further, the film scan data processing section  701  also outputs the information relating to a film size, kind of negative film and positive film, ISO (International Organization for Standardization) sensitivity optically or magnetically recorded in the film, name of manufacturer, main object, the information relating to the photographing condition (for example, a content of the written information of APS (Advanced-Photo System)) in addition to it to the image adjustment processing section  704 .  
      The reflection document scan data processing section  702  applies the correction operation proper to the reflection document input device  10 , negative, positive reversal in the case of the negative document, gray balance adjustment, contrast adjustment on the image information inputted from the reflection document input device  10 , and outputs it to the image adjustment processing section  704 .  
      The image data format decode processing section  703  applies the reversion of the compression sign, transform of the expression method of the color data, according to the data format of the image data inputted from the image transfer means  30  or communication means (input)  32 , and outputs it to the image adjustment processing section  704 .  
      The image adjustment processing section  704  applies each kind of image processing on the image signal expressing the color image inputted from the film scanner section  9 , reflection document input device  10 , image transfer means  30 , communication means (input)  32 , according to a command of the operating section  11  or control section  7 . Specifically, the image adjustment processing section  704  extracts the skin-colored area from the image signal when the skin-colored area of a person is included in the image signal.  
      Herein, the skin-colored area means a range of the skin such as a face or-neck of the person included in the image signal. Herein, the skin-colored area can be determined by an image area extracted by applying the simple range extension to the image signal in the lower frequency band components obtained by the dyadic wavelet transform on the image signal, or by an image area decided as the most suitable face by applying a logarithmic polar coordinate transform to around the most suitable eye of the object included in the image signal and applying the template matching or the simple area extension to the transformed image signal.  
      Further, the image adjustment processing section  704  converts the 3-color intensity signal of RGB of the image signal expressing the color image into the brightness signal and the color difference signal. The conversion into the brightness signal and the color difference signal may also be a method of converging an intensity signal for 3 colors of RGB representing the image signal of an object of the image processing into a intensity signal for the publicly known YIQ base or YUV base, or an intensity signal for XYZ base of the CIE 1931 color specification system or L*a*b base or L*u*V* base which is advised by CIE 1976 according to a standard such as sRGB (standard RGB) or NTSC (National TV Standards Committee). Further, it may also be a conversion such that an average of RGB values of the signal is defined to the brightness signal, and 2-axis perpendicular to the brightness signal is defined to a color difference signal.  
      Further, the image adjustment processing section  704  applies processing to reduce a variation amount of the signal intensity on pixels in which the spatial frequency is within a specific range, and the variation amount of the signal intensity is not larger than a specific threshold value. This reducing processing corresponds to the processing to remove the puckering or spotting included in the image signal of the higher frequency band components. Hereupon, even when the variation amount of the spatial frequency and the signal intensity is out of a specific range, a slightly weak reducing processing may also be applied to the variation amount of the signal intensity for the pixels in the vicinity of a specific range. Further, on the pixels not corresponding to even the pixels within the specific range and in the vicinity of the specific range, inversely, the processing to emphases the variation amount of the signal intensity may also be applied.  
      Herein, variation amounts of the spatial frequency and the signal intensity may be measured, by pasting the spatial frequency of a sinusoidal wave and a plurality of image signals on the image signal before the processing, using a soled retouch software, applying the image processing to the image signal, and measuring a change of an amplitude value after the processing to the amplitude value before processing.  
      Further, the image adjustment processing section  704  applies the adjustment of the contrast or a processing to add the noise signal.  
      Further, the image adjustment processing section  704  applies at least the second level dyadic wavelet transform to the image signal and decomposes the image signal into the frequency band components which are different from each other. Alternately, at least the first level multiresolution transform (for example, biorthogonal wavelet transform) which is a method of reducing the image size is applied to the image signal, and the image signal is decomposed into the image signal with various frequency band components, at least first level dyadic wavelet transform is further applied to the image signal of the lower frequency band components among decomposed components, the image signal of the lower frequency band components is decomposed into the image signal with the more various frequency band components. In this case, the image adjustment processing section  704  determines, the number of levels at which the the multiresolution transform which is a method of reducing the image size is decreased is switched to the dyadic wavelet transform according to the reading resolution of the image signal. According to the determined level number, the multiresolution transform and the dyadic wavelet transform are applied. Further, the inverse dyadic wavelet transform and the inverse multiresolution transform are applied to the image signal which has been applied the image processing.  
      The image adjustment processing section  704  outputs the image signal after the processing to the CRT proper processing section  705 , printer proper processing section  706 , printer proper processing section  707 , image data format making processing section  708 , and data accumulation means  71 .  
      The CRT proper processing section  705  applies a processing such as a pixel number change, color matching, on the image signal after the image processing inputted from the image adjustment processing section  704 , and outputs it to CRT  8  in accompany with each kind of display information.  
      The printer proper processing section  706  applies the correction processing proper to the printer, color matching, pixel number change, on the image signal after the image processing inputted from the image adjustment processing section  704 , and outputs it to the exposure processing section  4 .  
      When an external printer  34  such as an inkjet printer is connected to the image processing apparatus  1  of the present invention, the printer proper processing section  707  is provided for each of connected printers. This printer proper processing section  707  applies the correction processing proper to the printer, color matching, pixel amount change, on the image after the-image processi ng inputted from the image adjustment processing section  704 .  
      The image data format making processing section  708  transforms the image adjustment processing section  704  into each kind of general use image formats represented by JPEG (Joint Photographic Experts Group), TIFF (Tagged Image File Format), Exif (Exchangeable Image File Format) on the image signal after the image processing inputted, and outputs it to the image conveying section  31 , or communication means (output)  33 .  
      Hereupon, a division like as the film scan data processing section  701 , reflection document scan data processing section  702 , image data format decode processing section  703 , image adjustment processing section  704 , CRT proper processing section  705 , printer proper processing sections  706  and  707 , image data format making processing section  708 , is a division provided for a help of an understand of the functions of the image processing section  70  of the present embodiment, and it is not necessarily to be realized as a physically independent device. For example, it may also be realized as a division of the kind of a single software processing. Further, the image processing apparatus  1  in the present embodiment is not limited to the above-described content, but it can be applied to various embodiments such as a digital photo-printer, printer driver, plug-in of each kind of image processing software.  
      Next, as a specific execution method of the processing to remove the puckering or spotting, carried out in the image adjustment processing section  704  in  FIG. 12 , an example in which dyadic wavelet transform which is one of the wavelet transform, is used, will be described by dividing into Example 1 and Example 2.  
     EXAMPLE 1  
      Referring to  FIG. 13 , operations in the present example will be described below. Operations described in the present example are applied by the image adjustment processing section  704 . The present example 1 shows an example for removing the puckering or spotting is removed from the portrait, in which pixels forming the skin-colored area are extracted from pixels included in the image signal expressing the obtained color image, and the image signal expressing the color image is converted into the brightness signal and the color difference signal, and on the converted brightness signal, the third level dyadic wavelet transform is applied, and in the image signal of the frequency band components which are different from each other, obtained by level-2 and level-3 dyadic wavelet transforms, when the signal density of the image signal of the higher frequency band components of prior obtained pixels forming the skin-colored area is not larger than a specific threshold level, an example in which, by applying the processing to reduce the signal intensity.  
      Initially, pixels forming the skin-colored area are extracted from pixels included in the image signal, and the image signal is decomposed to the brightness signal S 0  and the color difference signal (not shown). On the brightness signal S 0 , the level-1 dyadic wavelet transform is applied through a filter group F 1 , and the brightness signal S 0  is decomposed to image signals W x1 , W y1  of the higher frequency band components and the image signal S 1  of the lower frequency band components. In the same manner, on the image signal S 1  of the lower frequency band components, the level-2 dyadic wavelet transform is applied through a filter group F 2 , and the image signal S 1  is decomposed to image signals W x2 , W y2  of the higher frequency band components and the image signal S 2  of the lower frequency band components, further, on the image signal S 2  of the lower frequency band components, the level-3 dyadic wavelet transform is applied through a filter group F 3 , and the image signal S 2  is decomposed to image signals W x3 , W y3  of the higher frequency band components and the image signal S 3  of the lower frequency band components.  
      In this case, from image signals W x2 , W y2 , W x3 , W y3  in the higher frequency band components decomposed by the level-2 and level-3 dyadic wavelet transforms, the standard deviations σ of the absolute value of the signal intensity of the image signals of pixels forming the previously extracted skin-colored area are respectively calculated, and the threshold values which are references of the removing processing of the puckering or spotting are determined. Then, when image signals of pixels forming the skin colored area included in each of image signals W x2 , W y2 , W x3 , W y3 , of the higher frequency band components have the signal intensity not larger than threshold values, the processing to reduce the signal intensity of pixels is applied (refer to F 4  in the drawing), the attenuation processed image signals W x2 ′, W y2 ′ W x3 ′, W y3 ′ are respectively generated.  
      Next, W x3 ′, W y3 ′ and the image signal S 3  of the lower frequency band components are level-3 inverse wavelet transformed through a filter group FS, and S 2 ′ is generated, next, W x2 ′, W y2 ′ and S 2 ′ are level-2 inverse wavelet transformed through a filter group F 6 , and S 1 ′ is generated, and next, W x1 , W y1  and S 1 ′ are level-1 inverse wavelet transformed through a filter group F 7 , and S 0 ′ is generated. Herein, S 0 ′ is a brightness signal in which the puckering or spotting is removed from the brightness signal S 0 . After that, the brightness signal S 0 ′ and the color difference signal (not shown) are converted into RGB signal, and the color image signal in which the puckering or spotting is removed, can be obtained.  
      Herein, filter groups F 1 -F 3  respectively apply level-1, 2, 3 dyadic wavelet transforms to the input signal, and filter groups F 5 -F 7  respectively apply level-3, 2, 1 inverse dyadic wavelet transforms to the input signal.  
      In each of filters of filter groups F 1 -F 3  and F 5 -F 7 , D_HPFkl, D_HPF′kl (k=1, 2, . . . n (n is a natural number); 1=x, y) are high pass filters, D_LPFkl, D_LPF′kl (k=1, 2, . . . n (n is a natural number); 1=x, y) are low pass filters. An example of the high pass filter and low pass filter is shown in Table 1.  
                               TABLE 1                       c   D_HPFl   D_LPFl   D_HPF′l   D_LPF′l                                                    −3           0.0078125   0.0078125       −2           0.054685   0.046875       −1       0.125   0.171875   0.1171875       0   −2.0   0.375   −0.171875   0.65625       1   2.0   0.375   −0.054685   0.1171875       2       0.125   −0.0078125   0.046875       3               0.0078125                  
 
 In Table 1, a filter coefficient of c=0 is a filter coefficient for pixels processing at this time, filter coefficient of c=−1 is a filter coefficient for one-preceding pixels to pixels processing at this time, and a filter coefficient of c=+1 is a filter coefficient for one-succeeding pixels to pixels processing at this time. 
 
      In the dyadic wavelet transform, the filter coefficient is different for each level. For the filter coefficient of level-i, a coefficient in which 2 i-1  zeros are inserted between respective coefficients of filters of level-1, is used.  
      Further, an example of the correction coefficient γ i  determined corresponding to the level-i of the dyadic wavelet transform is shown in Table 2.  
                           TABLE 2                                   i   γ                                                    1   0.66666667           2   0.89285714           3   0.97087379           4   0.99009901           5   1                      
 
      Hereupon, even when, the change of contrast is not applied to the color image signal from which the puckering and spotting obtained by the above processing are removed, there is a case where it is perceived as if the contrast is changed depending on conditions such as largeness and brightness of the face of a person included in the image signal. Therefore, the processing to change the contrast is applied to the image signal of the lower frequency band components obtained by the level-1 dyadic wavelet transform of pixels included in the skin-colored area. Alternatively, the processing to add the faint noise signal to the image signal is applied to the image signal of pixels included in the skin-colored area of the processed color image signal. When such a processing is applied, the contrast without visually a sense of incompatibility can be maintained. Hereupon, it is described herein that either one of the change of contrast or the addition of noise signal is applied, however, both processing may also be applied, or both processing may also not be applied.  
      In  FIG. 14 , an example of the image evaluation result in the present example 1 is shown.  FIG. 14  shows after the image processing of the present example 1 is applied to the image signal in which the object is recorded, after the face of 20-ages woman which is an object, is photographed in 2656 pixels*3992 pixels (approximate 10,600×10 3  pixels) by using Phase One H10 camera pack (made by Phase One co.) for Mamiya RZ67 Pro2 (made by Mamiya Co.) camera, the image evaluation result when it is outputted on the 2L-sized silver halide printing paper in the output resolving power of about 300 dpi.  
      In the evaluation result of  FIG. 14, 5  kinds of image processing (experiment 1-experiment 5) with various image processing conditions are different are applied, and the mean value of 5-grade evaluation by 10 subjects is made as the evaluation result. Herein, by indexes of a feeling of the puckering and a sense of incompatibility, a grade of removal of the puckering or spotting of the image, a sense of incompatibility for the image are respectively shown. Five kinds of image processing are attained by changing the level of the wavelet to reduce the signal intensity, the suppression amount to reduce the signal intensity, and the threshold value of the signal intensity to be reduced.  
      Herein, the image processing condition shows a range of the spatial frequency of the image signal, a range of a variation amount of the signal intensity to the maximum signal variation amount, a value of suppression amount to multiply a changing mount of the signal intensity to reduce the signal intensity. The spatial frequency corresponds to the frequency of the image signal of the higher frequency band components of each level obtained by the dyadic wavelet transform, and the variation amount range corresponds to the signal intensity of the image signal of the higher frequency band components obtained at each level. Herein, the maximum signal variation amount means the maximum value of the variation amount of the signal intensity (signal value), for example, in the case of 8-bit system, because the signal value of the image signal can take the range of 0-255, the maximum signal variation amount is 255.  
      According to the evaluation result shown in  FIG. 14 , the evaluation result when the image processing is applied by the image processing conditions of the experiment 1 and the experiment 3, is higher than the evaluation result by the other experiments 2, 4, 5. It indicates that it is preferable that the suppression processing to reduce the variation amount of the signal intensity to 0.5 times-0.75 times is applied, to the pixels in which the variation amount of the signal intensity is within a range of 0-20% of the maximum signal variation amount.  
      As described above, according to the image processing apparatus  1  of Example 1, when pixels forming the skin-colored area is extracted from pixels included in the obtained image signal, and the obtained image signal is decomposed to the brightness signal and the color difference signal, and the level-1 through level-3 dyadic wavelet transforms are applied to the decomposed brightness signal and color difference signal, and to pixels of the skin-colored area in which the spatial frequency of the image signal of the higher frequency band components obtained by the level-2 and level-3 dyadic wavelet transform is 1.5-3.0 lines/mm, and the variation amount of the signal intensity is within the range of 0-20% of the maximum signal variation amount, the processing to reduce the variation amount of the signal intensity is applied, because a sense of sharpness of the whole image is not deteriorated, further, the outline structure of the face of a person is not changed, the spotting or puckering can be removed without a sense of incompatibility.  
      Hereupon, in the present Example 1, although the processing to reduce the variation amount of the signal intensity is applied, it is not limited to this. A slightly weak suppression processing may also be applied to the variation amount of the signal intensity, for the pixels which are in the vicinity of a specific range, even when variation amounts of the spatial frequency and the signal intensity are out of a specific range. Further, inversely, the processing to emphasize the variation amount of the signal intensity may also be applied to pixels which does not correspond to also ones within a specific range and in the vicinity of a specific range.  
      Further, in the present Example 1, although the level-1 through level 3 dyadic wavelet transforms and the processing to reduce the signal intensity is applied to only the brightness signal, it is not limited to this. It may also be applied to only the color difference signal, and may also be applied to both the brightness signal and the color difference signal.  
      Further, it is preferable that the number of levels for applying the dyadic wavelet transform is appropriately changed correspondingly to the kind of objects included in the image signal, number of pixels of the image signal, and output size of the output resolving power.  
     EXAMPLE 2  
      Next, referring to  FIG. 15 , operations in Example 2 will be described. Operations described in the present Example 2 are applied by the image adjustment processing section  704 . The present Example 2 shows the first level biothogonal wavelet transform is applied to the image signal expressing the obtained color image, and in image signals of 2 different wavelength band components obtained by the 1-level biothogonal wavelet transform, from pixels included in the image signal of the lower frequency band components, pixels forming the skin-colored area are extracted, and the image signal of the lower frequency band components is converted into the brightness signal and the color difference signal, and the second level dyadic wavelet transform is applied to the converted brightness signal, and in the image signal of the frequency band components which are different from each other, obtained by the level-2 and level-3 dyadic wavelet transform, when the signal intensity of the image signal of the higher frequency band components of pixels forming the skin-colored area previously obtained, is not larger than a specific threshold value, by applying the processing to reduce the signal intensity, an example in which the puckering or spotting is removed from the portrait.  
      Initially, the level-1 biorthogonal wavelet transform is applied through the filter groups F 11   a  and F 11   b  to the image signal So. Herein, the image signal So is decomposed to the image signal W x1  of higher frequency band components and the image signal S x1  of lower frequency band components through the filter group F 11   a , and further, W x1  is decomposed to the image signals W d1  and W v1  of higher frequency band components through the filter group F 11   b , and S x1  is decomposed to the image signal W h1  of higher frequency band components and the image signal S 1  of lower frequency band components. Next, pixels forming the skin-colored area are extracted from pixels included in the decomposed image signal S 1  of the lower frequency components, and the image signal S 1  of the lower frequency components is converted into the brightness signal and the color difference signal. In  FIG. 15 , the converted brightness signal is re-defined as S 1 .  
      Next, on the brightness signal S 1 , the level-2 dyadic wavelet transform is applied through the filter group F 12 , and the brightness signal S 1  is decomposed to the image signals W x2 , W y2  of the higher frequency band components, and the image signal S 2  of the lower frequency components, further, on the image signal S 2  of the lower frequency components, the level-3 dyadic wavelet transform is applied through the filter group F 13 , and the image signal S 2  is decomposed to image signals W x3 , W y3  of the higher frequency band components and the image signal S 3  of the lower frequency band components.  
      In this case, from image signals Wx 2 , Wy 2 , Wx 3 , Wy 3  of the higher frequency band components decomposed by the level-2 and level-3 dyadic wavelet transforms, the standard deviations σ of the absolute values of the signal intensity of the image signal of pixels forming the skin colored area previously extracted, are respectively calculated, and the threshold value which is the reference of the removal processing of the puckering and spotting, is determined. Then, when the image signal of pixels forming the skin colored area included in each of image signals W x2 , W y2 , W x3 , W y3  of the higher frequency band components has the signal density not larger than the threshold value, the processing to reduce the signal intensity of pixels is applied (refer to the sign F 14  in the drawing), and the attenuation processed image signals W x2 ′, W y2 ′, W x3 ′, W y3 ′ are respectively generated.  
      Next, on W x3 ′, W y3 ′ and the image signal S 3  of the lower frequency band components, the level-3 inverse dyadic wavelet transform is applied through the filter group F 15 , and S 2 ′ is generated, next, on W x2 ′, W y2 ′ and the image signal S 2 ′, the level-2 inverse dyadic wavelet transform is applied through the filter group F 16 , and S 1 ′ is generated, and the brightness signal S 1 ′ and color difference signal (not shown) are converted into RGB signal, and the image signal S 1 ′ (re-defined) of the lower frequency band components in which the puckering and spotting have been removed, can be obtained.  
      Next, on W d1 , W v1 , W h1  and S 1 ′, the level-1 inverse biorthogonal wavelet transform is applied through the filter group F 17   a  and the filter group F 17   b . Herein, W d1  and W v1  are generated in W x1 , W h1  and S 1 ′ are generated in S x1 ′, respectively through the filter F 17   a , and further, W x1  and S x1 ′ are generated in the image signal S 0 ′ expressing the color image from which the puckering and spotting are removed, through the filter group  17   b.    
      Herein, the filter group F 11  applies the level-1 biorthogonal wavelet transform on the input signal, and the filter groups F 12  and F 13  respectively apply the level-2, 3 dyadic wavelet transform on the input signal. Further, filter groups F 15  and F 16  respectively apply the level-3, 2 inverse dyadic wavelet transform on the input signal, and the filter group F 17  applies the level-1 inverse biorthogonal wavelet transform.  
      In each of filters of the filter groups F 11  and F 17 , O_HPF 1 , O_LPFl (1=x, y) are respectively the high pass filter and low pass filter for biorthogonal wavelet transform, and O_HPF′ 1 , O_LPF′l are respectively the high pass filter and low pass filter for inverse biorthogonal wavelet transform. An example of the high pass filter and low pass filter is shown in Table 3.  
                               TABLE 3                       c   O_HPF   O_LPF   O_HPF′   O_LPF′                                                    −4       0.037829   −0.037829           −3   −0.064539   −0.023849   −0.023849   −0.064539       −2   0.04069   −0.110624   0.110624   −0.04069       −1   0.418092   0.377403   0.377403   0.418092       0   −0.788485   0.852699   −0.852699   0.788485       1   0.418092   0.377403   0.377403   0.418092       2   0.04069   −0.110624   0.110624   −0.04069       3   −0.064539   −0.023849   −0.023849   −0.064539       4       0.037829   −0.037829                  
 
 In Table 3, a filter coefficient of c=0 is a filter coefficient for pixels processing at this time, filter coefficient of c=−1 is a filter coefficient for one-pixels preceding pixels to pixels processing at this time, and a filter coefficient of c=+1 is a filter coefficient for one-pixels succeeding pixels to pixels processing at this time. Hereupon, for the simplification of the description, the filter coefficients D_HPFkl, D_HPF′kl, D_LPFkl, D_LPF′kl, (k=1, 2, . . . , n (n is a natural number); 1=x, y) for the dyadic wavelet transform and the inverse dyadic wavelet transform are expressed by the same name as Example 1, and an example of the filters is shown in Table 1. Further, an example of the correction coefficient γ i  is shown in Table 2. 
 
      Herein, there is a case where the processed color image signal is perceived as if the contrast is changed, depending on conditions such as the largeness of the face of a person included in the image signal, and the brightness even a case where the change of the contrast is not applied. Therefore, the processing to change the contrast is applied to the image signal of pixels included in the skin-colored area of the extracted lower frequency band components. Alternatively, the processing to add the faint noise signal to it, is applied to the image signal of pixels included in the skin-colored area of the processed color image signal. By applying such a processing, the contrast visually not having a sense of incompatibility can be maintained. Hereupon, it is described herein that either one of the change of contrast or the addition of noise signal is applied, however, both processing may also be applied, or both processing may also not be applied.  
      As described above, according to the image processing apparatus  1  of Example 2, the biorthogonal wavelet transform which is the multiresolution transform which is a method decreasing the image size is applied to the obtained image signal, and the skin-colored area is extracted from the image signal of the lower frequency band components obtained by the transform, and the image signal of the lower frequency band components is decomposed to the brightness signal and the color difference signal. Further, on the image signal of pixels corresponding to the skin-colored area of the image signal of the higher frequency band components obtained when the dyadic wavelet transforms are applied to the converted brightness signal, when the signal intensity is not larger than a specific threshold value, the processing to reduce the signal intensity and to suppress, is applied. When processed in this manner, the image size is decreased and the productivity of the processing can be increased, and because a sense of sharpness of the whole image is not deteriorated, further, the outline structure of the face of a person is not changed, the spotting and puckering can be removed without a sense of incompatibility.  
      Hereupon, in the present Example 2, although the level-2, 3 dyadic wavelet transforms and the processing to reduce the signal intensity are applied to only the brightness signal, it is not limited to this. It may also be applied to only the color difference signal, or may also be applied to both the brightness signal and the color difference signal.  
      Hereupon, the content of description in the detailed description to carry out the present invention, can be appropriately modified and changed within a scope which does not depart from the true spirit of the present invention.