Method of and apparatus for establishing highlight and shadow densities

Abstract: Highlight density (X H ) (or shadow density) is established on the basis of a cumulative histogram (H(X)) of density in an original: First, reference densities (X HR , X HL and X H1 - X Hn ) associated with prescribed reference statistic values (Y HR , Y HL and Y H1 - Y Hn ), respectively, and threshold values (W H1 - W Hn ) are computed. Then, the reference density (X HR ) is compared with a prescribed threshold density (X HS ), and a difference (ΔX Hi ) concerning the reference density (X Hi ) is compared with the threshold value (W Hi ) in a procedure, to thereby select one of procedures for establishing the highlight density.

Description: BACKGROUND OF THE INVENTION Field of the Invention The present invention relates to a method of and an apparatus for establishing highlight and shadow densities used for gradation control of image data obtained through reading an original. Description of the Prior Art In an image processing apparatus such as a process scanner or the like, an image data obtained through reading an original is subjected to gradation control, whereby a recorded image having a desired gradation expression is obtained. In the gradation control, a highlight density and a shadow density are needed to be previously set in the image processing apparatus in order to define both ends of that density range of the original which is faithfully reproduced with the recorded image. The highlight and shadow densities are respectively measured with a densitometer at a highlight area and a shadow area which are selected by an operator as the brightest area and the darkest area in the original. The measured highlight and shadow densities are set in the image processing apparatus to be used for the gradation control. Because the highlight and shadow densities are manually measured and set in the image processing apparatus by an operator so far, its operation takes much time and skill. Furthermore, a densitometer is necessary for the measurement besides the image processing apparatus. If the highlight area and the shadow area are not large enough for an aperture of the densitometer to be set in those areas, which aperture is usually a circle of a few millimeters in diameter, the highlight and shadow densities cannot be set at appropriate values. In order to overcome those problems, there has been proposed some techniques for establishing highlight and shadow densities properly and easily, see in this respect also EP-A-0255949. Japanese Patent Laying Open Gazette 1984-15939 According to the technique disclosed in this gazette, a highlight point and a shadow point where the highlight and shadow densities are to be measured are first specified in an original by an operator, respectively. Then, brightnesses of the respective highlight and shadow points are automatically read with a photoelectric converter device in a graphic arts scanner, to thereby be converted to the highlight and shadow densities in the process scanner. In order to prevent measurement error in the highlight and shadow densities because of dust or the like on the original, the brightnesses of the highlight and shadow points are determined as average values of brightnesses in the vicinity of the highlight and shadow points, respectively. The highlight and shadow densities are automatically determined on the basis of the average values. As a result, measurement error due to dust or the like is reduced and the highlight and shadow densities are obtained with considerable accuracy. Although appropriate highlight and shadow densities are established by this technique, manual operation is still necessary for specifying the highlight and shadow points. Japanese Patent Laying Open Gazette 1981-87044 Japanese Patent Laying Open Gazette 1986-20042 The applicant of the present invention is also the applicant of these gazettes. According to the technique disclosed in these gazettes, sizes of candidate areas for the highlight and shadow points are measured first. If the candidate areas have enough size for measuring accurate densities, they are selected as the highlight and shadow points, respectively. Then, densities of the candidate areas are determined to be highlight and shadow densities, respectively. However, an apparatus for only measuring the sizes of the candidate areas is necessary in the above techniques. Further, it is difficult to set highlight and shadow densities appropriately for such an original as has no candidate areas having enough sizes for measuring accurate densities. As described above, techniques in the art are insufficient for establishing proper highlight and shadow densities easily. There is another issue concerning highlight and shadow densities for an original where average density is shifted from a normal value. Highlight and shadow densities for such an abnormal original need to be set shifted from a normal value according to magnitude of the shift of the average density in order to reproduce a relatively normal image. If an area of low density does not exist in an original or if the area has only an extremely small size, highlight density should be set at a lower value than usual in order for a recorded image to be reproduced normally. On the other hand, if an area of high density does not exist or if the area has only an extremely small size, shadow density should be set at a higher value than usual in order for a recorded image to be reproduced normally. When an original image is abnormal where average density is shifted from a normal value, operation for establishing proper highlight and shadow densities depends on a skilled operator. Therefore, a method of and an apparatus for easily establishing proper highlight and shadow densities for an abnormal image has been an issue for years. SUMMARY OF THE INVENTION The invention is as set out in the claims. Accordingly, an object of the present invention is to provide a method of and an apparatus for easily establishing proper highlight and shadow densities according to density characteristics of an original. These and other objects, features, aspects and advantages of the present invention will become more apparent from the following detailed description of the present invention when taken in conjunction with the accompanying drawings. BRIEF DESCRIPTION OF THE DRAWINGS Figs. 1A and 1B are block diagrams showing the structure of a graphic arts scanner applied to preferred embodiments of the present invention. Figs. 2A and 2B are graphs showing a density histogram and a cumulative histogram, Figs. 3A and 3B are graphs showing cumulative histograms in a first embodiment of the present invention, Fig. 4 is a flow chart showing the procedure according to the first embodiment, Fig. 5 is a graph showing a gradation correction curve, Figs. 6A and 6B are enlarged views showing an end part of a cumulative histogram, Fig. 7 is a flow chart showing the procedure according to a second embodiment of the present invention, Figs. 8A through 8D are schematic views showing the procedure for establishing a highlight density according to the second embodiment, Figs. 9 and 11 are graphs showing a cumulative histogram in the automatic setting of highlight and shadow densities, Figs. 10 and 12 are flow charts showing the procedures according to a third embodiment of the present invention, Fig. 13 is a graph showing a gradation correction curve, Figs. 14A and 14B are flow charts showing another procedure according to the third embodiment, and Fig. 15 is a block diagram showing the structure of another apparatus applied to the preferred embodiments of the present invention. RELATIONSHIP BETWEEN SYMBOLS In the following description of preferred embodiments various, symbols expressing mathematical or phisical quantities are used in order to clarify quantitive character of the preferred embodiment. However, the symbols in the preferred embodiment are different from those used in the section of "SUMMARY OF INVENTION" and in the appended claims. This is because the symbols in the preferred embodiments are so provided as to express the quantities under embodied conditions, while those in the claims are provided for representing global concept. For reference, the relationship or correspondence between these symbols are listed in Table 1. DESCRIPTION OF THE PREFERRED EMBODIMENTS Figs. 1A and 1B are schematic block diagrams showing a graphic arts scanner in which an apparatus according to a preferred embodiment of the present invention is assembled. The graphic arts scanner shown in Fig. 1 comprises a transparent original-receiving glass plate 2 provided in an upper opening of a frame box 1, so that an original 3 is placed on the glass plate 2, being downwardly directed. The original-receiving glass plate 2 is provided in its upper end surface with a white reference plate 4 having reference optical density. The white reference plate 4 is employed for shading correction. Illumination light 6 from a light source 5, which is formed by a halogen lamp or the like, is reflected by the surface of the original 3 to be a reflected light 7 having image information. The reflected light 7 enters a first mirror 9 included in an optical system 8 to be reflected. The light 7 reflected by the first mirror 9 is successively reflected by second and third mirrors 10 and 11, to form an image on a light receiving surface of a charge-coupled device (CCD) line image sensor 13 serving as a photoelectric converter through an image forming lens 12. The CCD line image sensor 13 comprises a plurality of CCD elements one-dimensionally arrayed in a direction perpendicular to the plane of the drawing. The direction perpendicular to the plane of the drawing corresponds to a main scanning direction. The incident light focused on the light receiving surface of the CCD line image sensor 13 is photoelectrically converted into an image signal V CCD for each pixel by the CCD line image sensor 13. The image signal V CCD is digitized by an A-D converter 14, to be supplied to a shading correction circuit 15 for each pixel. The shading correction circuit 15 is adapted to correct nonuniformity in illuminaiton on the surface of the original 3, nonuniformity in image-forming function of the image-forming optical system 8 and nonuniformity in sensitivity of the respective CCD elements forming the CCD line sensor 13. The shading correction circuit 15 generates an image signal V N after the shading correction. The image signal V N is selectively supplied to a histogram counter circuit 17 or a gradation correction circuit 18 having a look-up table RAM 18a through a switching circuit 16. The histogram counter circuit 17 is employed to obtain a cumulative histogram as statistics expressing density characteristics of the original 3. As is described later in detail, highlight and shadow densities are determined on the basis of the result obtained in the histogram counter circuit 17, and a gradation correction table is generated using the highlight and shadow densities. The gradation correction table is stored in the RAM 18a. An image signal V₁ generated in the gradation correction circuit 18 is supplied to an image processing circuit 19, to be subjected to processing such as contour enhancement, that is, sharpness enhancement and magnification change in the image processing circuit 19. An image signal V₂ outputted from the image processing circuit 19 is supplied to a halftone signal generator 20. A halftone dot signal V dot generated in the halftone signal generator 20 serves as a modulation control signal for an acoustic optical modulator 24. The acoustic optical modulator 24 receives laser beam 22 from a laser source 21 through a mirror 23. The acoustic optical modulator 24 modulated the laser beam 22 on the basis of the halftone dot signal V dot , to provide an exposure beam 25. The exposure beam 25 is horizontally oscillated by vibration of a galvano mirror 26, to be emitted on a recording photosensitive material 28 through an image forming optical system 27 formed by an fϑ lens and the like. The vibration of the galvano mirror 26 is synchronous with output timing of the CCD elements in the CCD line sensor 13, whereby optical scanning in the main scanning direction α is achieved. The light source 5 and the first mirror 9 are fixed to a mechanism (not shown) transferring relatively to the original 3 in a direction A, whereby the original 3 is scanned in the direction A. In synchronization with this, the photosensitive material 28 is moved in the downward direction (-β), to thereby achieve reading subscanning in the direction A and recording subscanning in a direction β synchronously. This apparatus is also provided with a microcomputer 30 for controlling the respective circuits and processing the respective data, especially processing for establishing highlight and shadow densities, and processing for generating the gradation correction table as described later. The microcomputer 30 has a CPU 31 and a memory 32, and is connected to the A-D converter 14, the shading correction circuit 15, the switching circuit 16, the histogram counter circuit 17, the gradation correction circuit 18, the image processing circuit 19 and the halftone signal generator 20 through an I/O port 33. A keyboard panel 34 for inputting various data as hereinafter described is also connected to the I/O port 33. In a first stage for establishing highlight and shadow densities according to preferred embodiments of the present invention, a cumulative histogram about an object area in the original 3 is obtained through prescanning the original 3 as described below. The object area is a part of the original 3 as an object in reading and reproducing. First, the switching circuit 16 is switched to the histogram counter circuit 17 according to a switching signal SW₁ supplied from the microcomputer 30. Then, the object area in the original 3 is scanned and image signal V CCD is sequentially generated in the CCD line sensor 13 accordingly. The image data V CCD is A-D converted in the A-D converter 14 and then supplied to the shading correction circuit 15 to be subjected to shading correction. The shading correction is performed on the basis of reference data obtained by reading the white reference plate 4. The image signal V N , which is obtained as a result of the shading correction, is supplied to the histogram counter circuit 17 through the switching circuit 16. When the prescanning of the original 3 is completed, the histogram counter circuit 17 holds a density histogram shown in Fig. 2A for example. The abscissa axis of the graph is density X. The ordinate axis of the graph is frequency of occurence y of pixels having the density X. The frequency of occurence y is obtained by dividing sum of pixels having the density X by the total number of pixels. A density interval ΔX of the graph is a unit of density width used for constructing the density histogram. The histogram counter circuit 17 further generates a cumulative histogram h(X) shown in Fig. 2B. The ordinate axis of the graph is cumulative frequency of occurence Y concerning the density X. In the Figs. 2A and 2B and the other figures hereinafter described, symbols X min and X max are defined as follows: X min : Density at which the cumulative density histogram h(X) steps up from zero percent. X min is hereinafter referred to as "minimum density". X max : Density at which the cumulative density histogram h(X) reaches a hundred percent. X max is hereinafter referred to as "maximum density". Fig. 3A is a graph showing an example of the cumulative histogram to be processed in a first embodiment of the present invention. Fig. 3B is an enlarged view showing a part of Fig. 3A near the origin. A cumulative histogram h₀(X) is associated with an ordinary original which is appropriately exposed. A cumulative histogram h₁(X) is associated with an extraordinarily-dark (or "low-key") original such as a discolored original and an original underexposed in taking picture of a subject. The cumulative histograms h₀(X) and h₁(X) are simplified graphs showing examples of the histogram h(X) of Fig. 2B. A minimum value Y min and a maximum value Y max of cumulative frequency of occurence is zero and a hundred percent, respectively, for example. There are some items considered in a first embodiment of the present invention. First, since the minimum density X min is generally easy to fluctuate according to dust, defects or the like on an original, it is not good to adapt the minimum density X min for the highlight density in consideration of reproducibility of the gradation correction processing. Accordingly, the highlight density is usually set at a value a little higher than the minimum density X min . Second, if density character of an original is abnormal, value of the highlight density is determined as follows: The cumulative histogram h₁(X) of an extraordinarily-dark original is shifted to the right, that is, higher density side from the histogram h₀(X) of a normal original, and its minimum density X 1min is much larger than a minimum density X 0min of the normal original. For the extraordinarily-dark original, the highlight density is set at as small value as the minimum density X 1min in order for a brighter part to be emphasized in a reproduced image, whereby the reproduced image can be produced comparatively beautiful. The treatment described above is achievable only when it is recognizable whether an original to be treated is normal or not. However, it can be judged whether density character of the original is normal or not on the basis of a density value associated with a prescribed cumulative frequency in its histogram. Taking the above items into consideration, a reference density (RD) associated with a prescribed value of reference cumulative frequency (RCF) of occurrence is first obtained in a procedure of the first embodiment, and then the density character of the original is recognized through comparing the reference density with a prescribed threshold value. According to the density character, one of procedures prepared in advance for establishing the highlight density is selected as appropriate one to the original ,and the highlight density is obtained through the selected procedure. Fig. 4 is a flow chart showing the procedure according to a method of the first embodiment. At a step S1, a reference density X HR for highlight density ("highlight RD") and a minimum RD X HL associated with prescribed values of a reference cumulative frequency Y HR for highlight density ("highlight RCF") and a minimum RCF Y HL , respectively, are obtained on a cumulative histogram h(X). In Figs 3A and 3B, X HR =X 0HR and X HL =X 0HL in the histogram h₀(X), and X HR =X 1HR and X HL =X 1HL in the histogram h₁(X). The highlight RD X HR is used as an index showing how many pixels of comparatively low densities exist in the original. The minimum RD X HL is used as an index showing a density where the cumulative histogram h(X) substantially steps up. Accordingly, the minimum RCF Y HL is set as low as practicable and the highlight RCF Y HR is set higher than that. Detailed procedure for obtaining the highlight RD X HR and the minimum RD X HL is described later. At a step S2, the highlight RD X HR is compared with a prescribed threshold density X HS for highlight density, whereby it is judged whether the original is normal one or abnormal one to be subjected to exceptional processing. If the highlight RD X HR is lower than the threshold density X HS , the original is processed as a normal one at a step S3. For example, the following relation holds for the the cumulative histogram h₀(X): X 0HR < X HS In this case, the original is processed as a normal one at the step S3. The processing for a normal original may be performed according to a second embodiment of the present invention described later or some other method. On the other hand, the following relation holds for the cumulative histogram h₁(X): X HS < H 1HR Then, a step S4 is executed as exceptional processing for an abnormal original. At the step S4, a highlight density X 1H (X H ) is set at the minimum RD X 1HL obtained at the step S1. Because the highlight density X 1H is set close to the minimum density X 1min , pixels in a density range down to the highlight density X 1H is faithfully reproduced, that is, comparatively bright part of the original can be faithfully reproduced to have given gradation with halftone dots. Accordingly, a reproduced image can be produced beautifully. If a color film is used as an original, for example, the highlight RCF Y HR may be preferably set at 0.1 % - 2.0 %, especially 0.2 % - 0.8 %. The threshold density X HS may be preferably set at about 0.35. The minimum RCF Y HL may be preferably set at 0.01 % - 0.5 %, especially 0.02 % - 0.1 %. The values of Y HR and Y HL are empirically determined as cumulative frequencies close to a minimum value which can invariably specify one typical density of an original according to density character of the original without influence of dust, defects and the like on the original. Consequently, the values of Y HR and Y HL are not limited in the above stated ranges. Further, the values of Y HR and Y HL can be equal to each other. The value of X HS is also determined empirically and not limited in the above stated range. A shadow density for an extraordinarily-dark original such as an underexposed original, which has a histogram like h₁(X) shown in Fig. 3A, can be set according to a second embodiment or some other method for a normal original. It is because, as shown in Fig. 3A, a maximum density X 1max of the extraordinarily-dark original is not so much different from a maximum density X 0max of the normal original. Incidentally, the highlight density X H may be computed by a prescribed function f H of the densities X HR and X HL as follows: X H = f H (X HR , X HL ) For example, the highlight density X H may be computed by the following linear equation with a weighting coefficient k : X H = kX HR + (1-k) X HL where 0 < k < 1 According to the equation (2), the highlight density X H is determined at a value between the minimum RDX HL and the highlight RD X HR . As to an extraordinarily-bright (or "high-key") original such as an overexposed original and a faded original, a shadow density can be appropriately determined through similar processing to the above processing at higher density range of its cumulative histogram. A method according to a second embodiment of the present invention is explained along with Fig. 5 through 8. Fig. 5 is a graph showing a gradation correction curve F(X) employed for gradation control. The abscissa axis of the graph is density X, and the ordinate axis is halftone area rate Q in a recorded image. The gradation correction curve F(X) is determined as a curve drawn through coordinate points P H and P S which are specified by the values of the highlight density X H , the shadow density X S , a halftone area rate Q H for the highlight density and a halftone area rate Q S for the shadow density. As is known in the art, the halftone area rates Q H and Q L are specified so that halftone dots having area rates between those values Q H and Q L can be invariably reproducible. The values of Q H and Q L are set at 5% and 95%, respectively, for example. Accordingly, by using the gradation correction curve F(X), image areas of an original having densities between the highlight density X H and the shadow density X S are invariably reproduced with halftone dots having area rates between Q H and Q S . Further, image areas in a density range between the highlight density X H and a density X a associated with a minimum halftone area rate Q min (0% for example) can be reproduced with halftone dots expressing some gradation, although reliability and reproducibility of the halftone dots may be worse. This is true with image areas in a density range between the shadow density X S and a density X b associated with a maximum halftone area rate Q max (100% for example). Those density ranges have the width D H and D S as shown in Fig. 5, respectively. It can be seen that if the minimum density X min of Fig. 2B is equal to the density X a of Fig. 5 and the maximum density X max of Fig. 2B is equal to the density X b of Fig. 5, density range of an original is reproduced most widely and effectively. Hence, the main purpose of the second embodement is to equalize the minimum density X min and the density X a , and the maximum density X max and the density X b . However, the densities X a and X b are not directly set for the gradation correction curve F(X), but the highlight density X H and the shadow density X S are set for it. Accordingly, for the above purpose, the highlight and shadow densities X H and X S need to be set appropriately according to shapes of both ends of the cumulative histogram h(X) so that the densities X a and X b of the gradation correction curve F(X) becomes equal to the minimum density X min and the maxmum density X max of the cumulative histogram h(X), respectively. According to a method of the second embodiment, one of processing procedures prepared in advance is selected according to shape of end portion of the cumulative histogram h(X). The the cumulative histogram h(X) is processed by the selected procedure so that appropriate highlight and shadow densities are determined. Fig. 6A is an enlarged view showing a highlight part of the cumulative histogram h(X), where the histogram h(X) is approximated as a straight line. Now, two reference points P₀ (X H0 , Y H0 ) and P i (X Hi , Y Hi ) are selected on the histogram h(X), where the density X Hi is higher than the density X H0 . A density X im is defined as a density at which a straight line through the reference points P₀ and P i meets with the X axis. The density X im is equal to the munimum density X min because the histogram h(X) is also a straight line in the graph. Now a difference D i between the densities X Hi and X im is defined as: D i = X Hi - X im A difference ΔX Hi between the densities X Hi and X H0 is also defined as: ΔX Hi = X Hi - X H0 When the difference D i is equal to the width D H of the gradation correction curve F(X) shown in Fig. 5, the above stated purpose is achieved through setting the highlight density X H at the density X Hi . In this case, the relation between ΔX Hi and D H is derived as follows: First, the following equation holds because the histogram h(X) is approximated as a straight line: Y Hi - Y min Y H0 - Y min = D i D i - ΔX Hi In the case stated above, the following relation holds: D H = D i Solving for ΔX Hi using the equations (5) and (6): ΔX Hi = D H ( 1 - Y H0 - Y min Y Hi - Y min ) Because Y min is equal to 0% in usual, the following equation holds in place of the equation (7): ΔX Hi = D H (1 - Y H0 Y Hi ) The above purpose is achieved through setting the highlight density X H at the density X Hi which meets the equations (4) and (7) or (8). However the equation (7) or (8) does not necessarily hold strictly. For example, the righthand side of the equation (7) may be defined as W Hi : w Hi = D H ( 1 - Y H0 - Y min Y Hi - Y min ) If the value ΔX Hi obtained by the equation (4) agrees with the value W Hi obtained by the equation (9) within a prescribed error, the above purpose can be well achieved. Fig. 7 is a flow chart showing the procedure according to the second embodiment. Fig. 6B is an enlarged view showing a highlight part (low density part) of the cumulative histogram h(X). The point of the processing described below is that, as shown in Figs. 8A through 8D, a temporary value X₁, X₂ or X₃ of the highlight density X H is successively shifted to the lower density side of the cumulative histogram h(X) so that the highlight density X H is determined in order for the left side range having the width D H of the gradation correction curve F(X) to include the highlitht part of the cumulative histogram h(X) at which the histogram h(X) starts to increase. When the density X₂ is selected as the highlight density X H , as shown in Fig. 8C, a minimum density X 1m of the histogram h(X) corresponds to the density X a of the gradation correction curve F(X) so that the gradation corrective curve F(X) does not include useless part at its highlight end. The density X₂ does not satisfy the above point as shown in Fig. 8B and the density X₃ includes useless part between the densities X a and X 1m as shown in Fig. 8C. Detailed procedure of the processing for establishing the highlight density X H , as described above, is shown in Fig. 7. At a step S10, reference densities (RD) X H0 - X Hn corresponding to reference cumulative frequencies (RCF) Y H0 - Y Hn , respectively, are computed on the cumulative histogram h(X), where the RCF Y H0 - Y Hn have respective predetermined values at the lower density range and the subscript n is an integer. Then at a step S11, threshold values W H1 - W Hn are computed for respective RD's X H1 - X Hn by the equation (9). Fig. 6B shows the case in which n = 4 and five RCF's Y H0 - Y H4 are preset. Reference points P₀ - P₄ on the histogram h(X) have respective coordinate values of (X H0 , Y H0 ) - (X H4 , Y H4 ). The values X H0 and Y H0 , which are the smallest of all, work as basic values in the procedure shown in Fig. 7. Therefore, the values X H0 and Y H0 are referred to as "basic reference density (basic RD)" and "basic reference cumulative frequency (basic RCF)", respectively. The following relation among the RCF's Y H1 - Y H4 is determined for convenience with respect to the processing procedure described below: Y H4 < Y H3 < Y H2 < Y H1 In Fig. 6B, a density X im , where the subscript i is an integer, is defined as a density at which a straight line through points P i and P₀ intersects the X axis. If the end portion of the histogram h(x) is nearly a straight line, a difference between the density X im and the minimum density X min is small. At steps S12 and S13, the value of the subscript i is set at "1". Then a judgement whether the highlight density X H can be set at the value of the RD X H1 or not is substantially performed at a step S14. In this step, a difference ΔX H1 between the densities X H1 and X H0 is compared with the threshold value W H1 which is determined in advance in the step S11. If the difference ΔX H1 is smaller than or equal to the threshold value W H1 , the highlight density X H is set at the value of the density X H1 at a step S15, and the processing is ended. This is because a difference D₁ between the density X H1 and the minimum density X 1m is smaller than or equal to the difference D H , and because the initial purpose is achieved so that density range from the density X im to the density X H1 is reproduced with certain gradation. If the difference ΔX H1 is larger than the threshold value W H1 , the hilight X H is not set at the density X H1 and the processing returns from the step S16 to the step S13 where similar judgement on the next RD X H2 is performed. Because the difference ΔX H1 is larger than the threshold value W H1 , the next RD X H2 has to be smaller than the first RD X H1 . Consequently, the RCF Y H1 - Y H4 are predetermined according to the equation (10). In other words, the highlight density X H needs to be set as high as possible on the cumulative histogram h(X) in order that the density range having the width D H shown in Fig. 5 is utilized as much as possible in gradation control. Accordingly, the judgement is sequentially executed on the densities X H1 through X H4 in order of decreasing value. Then, the highlight density X H is set at a density X Hi when its difference ΔX Hi is smaller than or equal to its threshold value W Hi . If the subscript i is equal to a maximum number n (n = 4 in Fig. 6B), the value of the basic RD X H0 is selected as the highlight density X H at a step S17, and the processing is ended. Therefore, the basic RD X H0 is a minimum value to be as the highlight density X H . Accordingly, the smaller the RCF Y H0 is predetermined, the wider becomes a density range in which the value of the highlight density is selected. The RCF Y H0 is predetermined as a minimum cumulative frequency by which a proper density is invariably specified in the highlight portion of the histogram h(X) according to density characteristics of an original without reference to dust, defects and the like on the original. Preferably, the value of the basic RCF Y H0 may be set at 0.1 % - 0.5% and the RCF Y H1 - Y Hn are at 0.1 % - 2%. The threshold value W Hi is not always calculated by the equation (9), but it may be given by the following equation having a semi-empirical value α i : W Hi = D H ( 1 - Y H0 - Y min Y Hi - Y min ) + α i The value α i is given for each RD X i , where the subscript i is an integer, and the value α i is predetermined under consideration of the curvature at the end portion of the cumulative histogram h(X). In short, if the difference ΔX Hi agrees with the threshold value W Hi given by the equation (9) within the prescribed error α i or the former is smaller than the latter, the highlight density X H can be set at the value of the density X Hi . Although one of the RD's X H0 - X Hn is selected as the highlight density X H in the above mentioned procedure, other value calculated on the basis of the RD's X H0 - X Hn may be selected instead. For example, a straight line may be drawn through two reference points P₀ (X H0 , Y H0 ) and P₁ (X H1 , Y H1 ) as shown in Fig. 6B to intersect the X axis at the density X 1m . Then the highlight density X H may be set at the value of a density which is higher than the density X 1m by the difference D H . Thus, a proper value for the highlight density X H can be easily set according to the purpose described above. By a combination of the procedures of the first and second embodiments, appropriate highlight and shadow densities according to density characteristics of an original can be automatically established. Fig. 9 is a graph showing a cumulative histogram H(X) on which several values used in automatic setting of a highlight density are shown. Fig. 10 is a flow chart showing the procedure according to a third embodiment of the present invention. At a step S20, a highlight RD X HR , a minimum RD X HL and RD's X H1 - X Hn associated with a highlight RCF Y HR , a minimum RCF Y HL and RCF's Y H1 - Y Hn , respectively, are computed. This step includes the step S1 shown in Fig. 4 and the step S10 shown in Fig. 7. The highlight RCF Y HR of Fig. 9 works not on y as the basic RCF Y H0 in the second embodiment but also as the highlight RCF in the first embodiment. Therefore, the highlight RD X HR is the same with the basic RD X H0 in Fig. 9. Relations among the RCF's and among the RD's are as follows: Y HL < Y HR = Y H0 < Y H3 < Y H2 < Y H1 X HL ≦ X HR = X H0 ≦ X H3 ≦ X H2 ≦ X H1 Then, threshold values W H1 - W Hn are calculated by the equation (9) at a step S21. At a step 22, the highlight RD X HR is compared with a predetermined threshold density X HS . As described in the first embodiment, the highlight RCF Y HR is predetermined preferably at a range of 0.1 % - 2.0 %, and a threshold density X HS is at about 0.35. If the highlight RD X HR is judged to be higher than the threshold density X HS , the original is subjected to a exceptional processing at a step S23. Namely, the minimum RD X HL is selected as the highlight density X H at the step S23 and the overall processing is ended. The step S23 corresponds to the step S4 in Fig. 4. Since the highlight RD X HR is lower than the threshold density X HS in Fig. 9, the original is judged to have normal density characteristics at the highlight portion of the histogram H(X), and steps S24 - S29 are executed. The steps S24 - S29 as a whole corresponds to the step S3 of Fig. 4 as well as the steps S12 - S17 of Fig. 7. A parameter i for subscript in Fig. 10 is initialized to be "1" at the steps S24 and S25. A density difference ΔX Hi (=X Hi - X HR ) is compared with a threshold value W Hi at the step S26 which is computed in advance at the step S21. If the difference ΔX Hi is smaller than or equal to the threshold value W Hi , the step S27 is executed that the RD X Hi is selected as the highlight density X H . On the other hand, if the difference ΔX Hi is greater than the threshold value W Hi , the next RD X Hi+1 is subjected to the processing on and after the step S25 again. If the RD X Hi is the last RD X Hn at the step S28, where n = 3 in Fig. 9, the step S29 is executed that the highlight RD X HR is selected as the highlight density X H . Preferably, the minimum RCF Y HL is predetermined at a range of 0.01 % - 0.5%, the RCF's Y H1 - Y Hn are at a range of 0.1 % - 2%, and the threshold density X HS is at about 0.35, respectively, as described before. As has been described above, appropriate highlight density X H can be easily set in consideration of density characteristics of an original according to the third embodiment including methods of the first and second embodiments. Fig. 11 is a graph showing the cumulative histogram H(X) with several values used in automatic setting of a shadow density. Fig. 11 corresponds to Fig. 9 for highlight density setting. Fig. 12 is a flow chart showing the similar procedure to that of Fig. 10, applied to shadow density setting. The steps S21 - S29 of Fig. 10 correspond to steps S210 - S290 of Fig. 12, respectively. The values shown in Figs. 11 and 12 correspond to those shown in Figs. 9 and 10, respectively, as follows: Reference Cumulative Frequency Y HR → Y SR Reference Density X HR → X SR Threshold Density X HS → X SS Reference Cumulative Frequency Y Hi → Y Si (i=1-3) Reference Density X Hi → X Si (i=1-3) Minimum Reference Cumulative Frequency Y HL → Maximum Reference Cumulative Frequency Y SL Minimum Reference DensiTy X HL → Maximum Reference Density X SL Minimum Cumulative Frequency Y min → Maximum Cumulative Frequency Y max Density Width D H → D S Threshold Value W Hi → W Si The relations (12) and (13) are replaced with the following relations (14) and (15): Y SL > Y SR = Y S0 > Y S3 > Y S2 > Y S1 X SL ≧ X SR = X S0 ≧ X S3 ≧ X S2 ≧ X S1 Preferably, the shadow RCF Y SR is predetermined at a range of 99.5 % - 99.9%, the shadow RD Y Si is at a range of 98 % - 99.9 %, the maximum RCF Y SL is at a range of 99.5 % - 99.9% and the threshold density X SS is at about 2.0, respectively. The procedure for establishing shadow density shown in Figs. 11 and 12 is almost the same with that for highlight density shown in Fig. 9 and 10, and therefor detail explanation will be omitted. A gradation correction table is generated on the basis of the highlight and shadow densities obtained according to the above procedures. For example, data expressing relation between density X and halftone area rate Q, which relation is shown as a gradation correction curve G(X) in Fig. 13, is stored in the RAM 18a of the gradation correction circuit 18 shown in Fig. 1 in form of a look-up table. The gradation correction curve G(X) is drawn through points P H and P S which are defined by a combination of the highlight density X H and a halftone area rate Q H therefor, and that of the shadow density X S and a halftone area rate Q S therefor, respectively. The halftone area rates Q H and Q S are predetermined in advance. After a gradation correction table expressing the gradation correction curve G(X) is stored in the RAM 18a, the switching circuit 16 is switched to the gradation correction circuit 18. Then the original 3 is scanned and read again, whereby image data V N thus obtained is given to the gradation correction circuit 18. The image data V N is corrected in the circuit 18 according to the gradation correction curve G(X), to thereby be changed into a corrected image data V₁. Recorded image is finally produced as described before on the basis of the corrected image data V₁. In the first embodiment, the reference density X HR is compared with a single value of the threshold density X HS . However, the reference density X HR may be compared with plurality of threshold densities in order to select one of plurality of procedures preset associated with the plurality of threshold densities, respectively. Although a shape of the cumulative histogram is approximated as a straight line at its end portion in the first embodiment and the value of the highlight density (or the shadow density) is determined on the basis of a difference between two of reference densities, the shape of the histogram may be approximated to a line expressed by a polynominal or the like. In such a case, the highlight density X H may be calculated by a function of the values of X H0 - X Hn and Y H0 - Y Hn as follows: X H = X H ( X H0 , ..., X Hn , Y H0 , ..., Y Hn ) The highlight density X H may be also set on the basis of a differencial coefficient at the end portion of the histogram h(X). Further, statistics r(X) expressed by the following equation may be utilized in place of the cumulative histogram h(X): Since a change rate of staistic value due to a density change on the statistics r(X) is less than on the histogram h(X), the highlight and shadow densities are established more steadily without reference to noise component on the original if the statistics r(X) is utilized instead of the histogram h(X). In general, the highlight and shadow densities may be set according to a shape of an end portion of density statistics of an original. The method of the second embodiment can be applied to another procedure where different processings are performed on a normal original, an underexposed original and an overexposed original, respectively. Fig. 14 is a flow chart showing this procedure, and it is partly transformed from Fig. 10. Namely, a step S30(S31 - S34) is inserted between the step S21 and S22 in Fig. 10. The cumulative histogram h(X) of Fig. 9 is also subjected to the procedure of Fig. 14. At a step S31 in Fig. 14, a threshold value W HL associated with a minimum RD X HL is calculated by the following equation: W HL = D H ( 1 - Y H0 - Y min Y HL - Y min ) Then, the absolute value of the threshold W HL is compared with that of the difference ΔX HL (= X HL - X H0 ) at a step S32. The reason why the absolute values are compared in the step S32 is that both the values of ΔX HL and W HL are negative because relations X HL < X H0 and Y HL < Y H0 hold. If the absolute value of ΔX HL is larger than that of W HL , the original is processed as a normal one by the steps S22 - S29. The value of ΔX HL expresses a gradient of the histogram H(X) at a density range from X HL to X H0 and the step S32 means that the original is processed as a normal one if the gradient is greater than a certain value. On the other hand, if the gradient is smaller than the certain value, the highlight end portion of the histogram H(X) is gentle to have a wide density range thereabout where the cumulative frequency is fairly low. Hence, the original is subjected to an exceptional processing for an abnormal original. Namely, the minimum RD X HL is first compared with a predetermined threshold density X HLS at a step S33. If the minimum RD X HL is higher than the threshold density X HLS , the step S23 is executed that the minimum RD X HL is selected as the highlight density X H . This is similar to the step S22 and S23 in Fig. 10. If the minimum X HL is lower than or equal to the threshold density X HLS , the highlight density X H is set at the average value X LR of the minimum RD X HL and the highlight RD X HR at a step S34. It is considered in this procedure that the highlight density X H would not be set excessively low. The value X LR may be calculated by the following equation at the step S34: X LR = k X HL + (1 - k) X HR where 0 ≦ k ≦ 1. Although a combination of the methods of the first and second embodiment is used in the third embodiment in order to establishing the highlight and shadow densities automatically, those methods can be independently employable and produce respective effects separately. However, if the combination of the methods is used as in the third embodiment, appropriate highlight and shadow densities can be automatically established, whether density character of an original is normal or not. Further, other combinations of the methods of the first and second embodiments are available which are different from that of the third embodiments. The apparatus used for the present invention is not restricted to the graphic arts scanner shown in Fig. 1, but the present invention may be applicable to other apparatuses such as a drum type image scanner, facsimile terminal equipment, a copying machine, an image transmission apparatus and the like. The image data to be processed by the present invention is not restricted to what is obtained through scanning and reading by a graphic arts scanner as described in the embodiments, but various types of image data such as those obtained by so-called High Definition Television and other apparatuses referred to as New Media, can be also processed by the present invention. Fig. 15 is a block diagram showing the structure of one of such apparatuses. A magnetic tape 41 stores image data obtained by the other apparatus. The magnetic tape 41 is controlled by the magnetic tape controller 42 to transmit the image data V N to the switching circuit 16 or to the I/O port 33 of the microcomputer 30. That is, the image processing apparatus of Fig. 15 is provided with the magnetic tape 41 and the magnetic tape controller 42 in place of the internal structure of the frame box 1, the A-D converter 14 and the shading correction circuit 15 of the graphic arts scanner shown in Fig. 1, and the other structure of the apparatus is the same with that shown in Fig. 1. The magnetic tape 41 and the controller 42 may be also replaceable by a magnetic disc, an optical disc, other memory means with controller and an apparatus for generating image data such as High Definition Television. According to the present invention, appropriate highlight and shadow densities can be easily established in consideration of density characteristics of an original because the processing procedure is selected on the basis of a certain density of or a shape of density statistics of the original. Further, the present invention does not need manual operation for measuring densities, thereby reducing time and labor of an operator remarkably and providing a recorded image and a reproduced image having invariable quality.

Claims: A method of establishing a gradation control density X c employed as a highlight density or a shadow density on a gradation correction curve F(X) which is used in gradation correction of image data expressing density of each pixel in an image, said method comprising the steps of: (a) on the basis of said image data, generating a statistical curve Y = S(X) expressing a cumulative density distribution on said image, where said statistical curve Y = S(X) is defined on a X - Y coordinate plane having a X-axis representing density levels and a Y-axis representing a cumulative number of pixels, (b) finding a line approximate to an end portion of said statistical curve Y = S(X), (c) presetting a first difference D between said gradation control density X c and a terminal density X a of said gradation correction curve F(X), (d) finding a terminal value Y m on said Y-axis at which said statistical curve Y = S(X) is terminated, (e) computing a terminal density X im of said line at which said line reaches said terminal value Y m , and (f) establishing said gradation control density X c so that said first difference D agrees substantially with a second difference D i between said gradation control density X c and said terminal density X im of said line. A method in accordance with claim 1, wherein said line is a straight line drawn through a first point P₀ and a second point P i , where said first point P₀ is specified by a combination of a prescribed first reference value Y 2R on said Y-axis and a first reference density X 2R computed by converting said first reference value Y 2R through said statistical curve Y = S(X) and said second point P i is specified by a combination of a prescribed second reference value Y i on said Y-axis and a second reference density X i computed by converting said second reference value Y i through said statistical curve Y = S(X), while said second reference value Y i is preset farther from said terminal value Y m than said first reference value Y 2R is. A method in accordance with claim 1 or 2, characterized in that the method comprises the steps of: (1) on the basis of said image data, generating a statistical curve Y = S(X) expressing a cumulative density distribution on said image, where said statistical curve Y = S(X) is defined on a X - Y coordinate plane having a X-axis representing density levels and a Y-axis representing a cumulative number of pixels, (2) providing a threshold density X TH on said X-axis, (3) finding a terminal value Y m on said Y-axis at which said statistical curve Y = S(X) is terminated, (4) providing a first reference value Y 1R on said Y-axis, (5) finding a terminal density X m on said X-axis at which said statistical curve Y = S(X) reaches said terminal value Y m , (6) converting said first reference value Y 1R through said statistical curve Y = S(X), to thereby obtain a first reference density X 1R , (7) comparing said first reference density X 1R with said threshold density X TH , to thereby designate said statistical curve Y = S(X) as a first type corresponding to an extraordinarily-bright image or an extraordinarily-dark image when said first reference density X 1R is farther from said terminal density X m than said threshold density X TH is, or as a second type corresponding to an ordinary image other than said extraordinarily-bright image and said extraordinarily-dark image when said first reference density X 1R is closer to said terminal density X m than said threshold density X TH is, (8) selecting a first procedure for setting said gradation control density X c at a terminal reference density X TR predetermined between said terminal density X m and said first reference density X 1R when said statistical curve Y = S(X) is said first type, or selecting a second procedure other than said first procedure for setting said gradation control density X c when said statistical curve Y = S(X) is said second type, and (9) establishing said gradation control density X c through said first procedure or said second procedure selected at said step (8). An apparatus for establishing a gradation control density X c employed as a highlight density or a shadow density on a gradation correction curve F(X) which is used in gradation correction of image data expressing density of each pixel in an image, said apparatus comprising: (1) a first processor means for generating a statistical curve Y = S(X) expressing a cumulative density distribution on said image on the basis of said image data, where said statistical curve Y = S(X) is defined on a X - Y coordinate plane having a X-axis representing density levels and a Y-axis representing a cumulative number of pixels, (2) a second processor means for finding a line approximate to an end portion of said statistical curve Y = S(X), (3) a third processor means for computing a terminal density X im , at which said line is terminated, and (4) a fourth processor means for establishing said gradation control density X c so that a first difference D between said gradation control density X c and a terminal density X a of said gradation correction curve F(X) agrees substantially with a second difference D i between said gradation control density X c and said terminal density X im of said line. An apparatus in accordance with claim 4, characterized in that said apparatus comprises: (1) a first processor means for generating a statistical curve Y = S(X) expressing a cumulative density distribution on said image on the basis of said image data, where said statistical curve Y = S(X) is defined on a X - Y coordinate plane having a X-axis representing density levels and a Y-axis representing a cumulative number of pixels, (2) a second processor means for finding a terminal density X m on said X-axis at which said statistical curve Y = S(X) is terminated, (3) a third processor means for computing a first reference density X 1R by converting a prescribed first reference value Y 1R on said Y-axis through said statistical curve Y = S(X), (4) a fourth processor means for comparing said first reference density X 1R with a prescribed threshold density X TH , to thereby designate said statistical curve Y = S(X) as a first type corresponding to an extraordinarily-bright image or an extraordinarily-dark image when said first reference density X 1R is farther from said terminal density X m than said threshold density X TH is, or as a second type corresponding to an ordinary image other than said extraordinarily-bright image and said extraordinarily-dark image when said first reference density X 1R is closer to said terminal density X m than said threshold density X TH is, (5) a fifth processor for establishing said gradation control density X c through a first procedure for setting said gradation control density X c at a terminal reference density X TR predetermined between said terminal density X m and said first reference density X 1R when said statistical curve Y = S(X) is said first type, or through a second procedure other than said first procedure for setting said gradation control density X c when said statistical curve Y = S(X) is said second type.

Inventors: NOMURA AKIHIRO DAINIPPON SCREE, NOMURA, AKIHIRO DAINIPPON SCREEN MFG. CO., LTD.
Applicants: DAINIPPON SCREEN MFG, DAINIPPON SCREEN MFG. CO., LTD.
Application year: 1989
Main classification: H04N1
IPC codes: G06T5, G03F3, H04N1