WO2015011991A1 - 閾値データ設定装置、方法及びプログラム、並びに画像形成システム - Google Patents

閾値データ設定装置、方法及びプログラム、並びに画像形成システム Download PDF

Info

Publication number
WO2015011991A1
WO2015011991A1 PCT/JP2014/064694 JP2014064694W WO2015011991A1 WO 2015011991 A1 WO2015011991 A1 WO 2015011991A1 JP 2014064694 W JP2014064694 W JP 2014064694W WO 2015011991 A1 WO2015011991 A1 WO 2015011991A1
Authority
WO
WIPO (PCT)
Prior art keywords
image signal
threshold
halftone dot
response characteristic
threshold data
Prior art date
Application number
PCT/JP2014/064694
Other languages
English (en)
French (fr)
Inventor
修平 堀田
Original Assignee
富士フイルム株式会社
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 富士フイルム株式会社 filed Critical 富士フイルム株式会社
Priority to EP14829027.3A priority Critical patent/EP3026882B1/en
Publication of WO2015011991A1 publication Critical patent/WO2015011991A1/ja
Priority to US14/990,679 priority patent/US9407788B2/en

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N1/00—Scanning, transmission or reproduction of documents or the like, e.g. facsimile transmission; Details thereof
    • H04N1/40—Picture signal circuits
    • H04N1/405—Halftoning, i.e. converting the picture signal of a continuous-tone original into a corresponding signal showing only two levels
    • H04N1/4055—Halftoning, i.e. converting the picture signal of a continuous-tone original into a corresponding signal showing only two levels producing a clustered dots or a size modulated halftone pattern
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B41—PRINTING; LINING MACHINES; TYPEWRITERS; STAMPS
    • B41C—PROCESSES FOR THE MANUFACTURE OR REPRODUCTION OF PRINTING SURFACES
    • B41C1/00—Forme preparation
    • B41C1/02—Engraving; Heads therefor
    • B41C1/04—Engraving; Heads therefor using heads controlled by an electric information signal
    • B41C1/05—Heat-generating engraving heads, e.g. laser beam, electron beam
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06K—GRAPHICAL DATA READING; PRESENTATION OF DATA; RECORD CARRIERS; HANDLING RECORD CARRIERS
    • G06K15/00—Arrangements for producing a permanent visual presentation of the output data, e.g. computer output printers
    • G06K15/02—Arrangements for producing a permanent visual presentation of the output data, e.g. computer output printers using printers
    • G06K15/18—Conditioning data for presenting it to the physical printing elements
    • G06K15/1867—Post-processing of the composed and rasterized print image
    • G06K15/1872—Image enhancement
    • G06K15/1876—Decreasing spatial resolution; Dithering
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06K—GRAPHICAL DATA READING; PRESENTATION OF DATA; RECORD CARRIERS; HANDLING RECORD CARRIERS
    • G06K15/00—Arrangements for producing a permanent visual presentation of the output data, e.g. computer output printers
    • G06K15/02—Arrangements for producing a permanent visual presentation of the output data, e.g. computer output printers using printers
    • G06K15/18—Conditioning data for presenting it to the physical printing elements
    • G06K15/1867—Post-processing of the composed and rasterized print image
    • G06K15/1872—Image enhancement
    • G06K15/1881—Halftoning
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B41—PRINTING; LINING MACHINES; TYPEWRITERS; STAMPS
    • B41C—PROCESSES FOR THE MANUFACTURE OR REPRODUCTION OF PRINTING SURFACES
    • B41C1/00—Forme preparation
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B41—PRINTING; LINING MACHINES; TYPEWRITERS; STAMPS
    • B41M—PRINTING, DUPLICATING, MARKING, OR COPYING PROCESSES; COLOUR PRINTING
    • B41M1/00—Inking and printing with a printer's forme
    • B41M1/02—Letterpress printing, e.g. book printing
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B41—PRINTING; LINING MACHINES; TYPEWRITERS; STAMPS
    • B41M—PRINTING, DUPLICATING, MARKING, OR COPYING PROCESSES; COLOUR PRINTING
    • B41M1/00—Inking and printing with a printer's forme
    • B41M1/02—Letterpress printing, e.g. book printing
    • B41M1/04—Flexographic printing
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B41—PRINTING; LINING MACHINES; TYPEWRITERS; STAMPS
    • B41P—INDEXING SCHEME RELATING TO PRINTING, LINING MACHINES, TYPEWRITERS, AND TO STAMPS
    • B41P2200/00—Printing processes
    • B41P2200/10—Relief printing
    • B41P2200/12—Flexographic printing

Definitions

  • the present invention relates to a threshold data setting device, method and program for setting threshold data for generating a binary image signal representing the shape and arrangement of each halftone dot constituting an image based on a continuous tone image signal, and image formation About the system.
  • CTP Computer-to-Plate
  • a screening process including a systematic dither method is performed on a continuous tone image signal.
  • various methods for determining threshold data considering characteristics of the whole or a part of the printing system have been proposed.
  • Patent Document 1 there is a method for determining a threshold array in which moiré or the like is unlikely to occur by sequentially performing (1) pseudo-generation of dots and (2) evaluation of a specific frequency component of a density pattern including the dots. Proposed.
  • types of printing devices or printing materials coloring materials, printing media, printing plates, etc.
  • types of printing apparatuses not only offset printing machines but also gravure printing machines, flexographic printing machines, digital printing machines (on-demand printing machines), and the like exist.
  • types of printing media include resins such as vinyl chloride and PET (polyethylene terephthalate), tarpaulins, and metal sheets.
  • Patent Document 1 since it is premised on creating threshold data suitable for each combination of a printing apparatus and a printing material, a lot of work man-hours are required for the preparation and creation processing, There was plenty of room for improvement.
  • An object of the present invention is to provide a threshold data setting apparatus, method and program, and an image forming system.
  • a threshold data setting device is a device for setting threshold data for creating a binary image signal representing the shape and arrangement of each halftone dot constituting an image based on a continuous tone image signal, Obtain halftone dot information related to the shape and arrangement of each halftone dot associated with the continuous tone image signal, A threshold candidate acquisition unit that acquires a plurality of threshold data matching the halftone dot information as each threshold candidate, and a halftone dot shape generated in the image forming process from the creation of the binary image signal to the formation of the image.
  • the change in the halftone dot shape is simulated by acting the response characteristic acquisition unit that acquires the change as a response characteristic in the spatial frequency domain and the response characteristic acquired by the response characteristic acquisition unit.
  • a simulation image creating unit that creates a simulated image signal, and an evaluation value obtained by executing a predetermined evaluation process on the simulated image signal created by the simulated image creating unit as a quantitative value indicating the degree of image quality
  • the halftone dot shape can be easily and accurately reproduced on the simulated image. Since a plurality of threshold data acquired as each threshold candidate is compared and / or evaluated based on a quantitative value indicating the degree of image quality, even if the number of types of printing apparatuses and printing materials is enormous, Threshold data suitable for the combination can be determined easily and easily.
  • the simulated image creating unit creates the simulated image signal including at least one flat mesh region, and the quantitative value calculating unit relates to the power spectrum of the simulated image signal in the flat mesh region. It is preferable to calculate the quantitative value based on the main component excluding the component corresponding to the point.
  • the simulated image creating unit creates the simulated image signal including two or more flat mesh regions having different gradation levels
  • the quantitative value calculating unit is a statistical value for the main component in each of the flat mesh regions. It is preferable to calculate the quantitative value based on the above.
  • a response characteristic estimating unit that estimates the response characteristic according to an event that changes the halftone dot shape in the image forming process, and the response characteristic acquiring unit is estimated by the response characteristic estimating unit; It is preferable to acquire the response characteristic.
  • the response characteristic estimation unit estimates the response characteristic when the color material is transferred to the printing material via the printing plate.
  • the response characteristic estimation unit estimates the response characteristic when the printing plate is formed by engraving the main surface of the plate material.
  • a binarization processing unit that converts the continuous tone image signal into the binary image signal by setting the threshold data determined by the threshold data determination unit.
  • An image forming system includes the threshold data setting device and an image forming device that forms the image based on the binary image signal converted by the threshold data setting device.
  • a threshold data setting method is a method for setting threshold data for creating a binary image signal representing the shape and arrangement of each halftone dot constituting an image based on a continuous tone image signal.
  • Obtaining halftone dot information relating to the shape and arrangement of each halftone dot associated with a continuous tone image signal, obtaining a plurality of threshold data matching the halftone dot information as each threshold candidate, and the binary image By obtaining the change in the halftone dot shape that occurs in the image formation process from the creation of the signal to the formation of the image as a response characteristic in the spatial frequency domain, and by operating the acquired response characteristic, A step of creating a simulated image signal in which the change in the halftone dot shape is simulated, and an evaluation value obtained by executing a predetermined evaluation process on the created simulated image signal, By applying the obtained threshold candidates and comparing and / or evaluating based on the calculated quantitative values, threshold data to be set from the threshold candidates is calculated. And causing the computer to execute the determining step.
  • a threshold data setting program is a program for setting threshold data for creating a binary image signal representing the shape and arrangement of each halftone dot constituting an image based on a continuous tone image signal.
  • To create a simulated image signal in which the change in the halftone dot shape is simulated, and an evaluation obtained by executing a predetermined evaluation process on the created simulated image signal Is calculated as a quantitative value indicating the degree of image quality, and each of the obtained threshold candidates is applied, and is set from among the threshold candidates by comparison and / or evaluation based on the calculated quantitative values.
  • causing the computer to execute the step of
  • the change in halftone dot shape that occurs in the image forming process from the creation of the binary image signal to the formation of the image is represented by the spatial frequency.
  • Threshold data suitable for the combination can be determined easily and easily.
  • FIG. 1 is an overall configuration diagram of an image forming system according to this embodiment. It is a schematic explanatory drawing regarding operation
  • FIG. 12A is a graph illustrating an example of a power spectrum in a simulated image signal.
  • FIG. 12B is a graph showing the peak characteristics of the power spectrum in each threshold candidate. It is a table which shows the quantitative result with respect to each threshold value candidate. It is a whole block diagram of a response characteristic estimation system.
  • FIG. 15A is a schematic front view of a sample image halftone visualized.
  • FIG. 15B is a front view of the sample letterpress formed based on the sample image halftone dot of FIG. 15A.
  • FIG. 15C is a front view of a sample printed matter formed based on the sample letterpress of FIG. 15B.
  • FIG. 16A is a first explanatory diagram relating to an estimation model of the height shape of a halftone protrusion.
  • FIG. 16B is a graph showing the Fourier transform characteristics of the height shape shown in FIG. 16A.
  • FIG. 17A is a second explanatory diagram relating to an estimation model of the height shape of the halftone protrusions.
  • FIG. 17B is a graph showing the Fourier transform characteristics of the height shape shown in FIG. 17A.
  • a threshold data setting method according to the present invention will be described in detail with reference to the accompanying drawings by giving preferred embodiments in relation to a threshold data setting device, a program, and an image forming system for executing the threshold data setting method.
  • forming an image may be referred to as “printing” or “printing”.
  • FIG. 1 is an overall configuration diagram of an image forming system 10 according to this embodiment.
  • the image forming system 10 forms a relief plate 14 (printing plate) by engraving a plate-like plate material 12, and transfers ink (coloring material) to a printing medium 16 through the relief plate 14 to produce a printed matter.
  • This is a system that produces 18.
  • This printed matter 18 is formed by forming at least one color halftone dot on the printing medium 16.
  • the image forming system 10 includes an image processing device 20 as a threshold data setting device, a relief forming device 22, a flexographic printing machine 24, a threshold database (hereinafter referred to as threshold DB 26), and a response characteristic database (hereinafter referred to as response characteristic DB 28).
  • an image processing device 20 as a threshold data setting device
  • a relief forming device 22 a flexographic printing machine 24, a threshold database (hereinafter referred to as threshold DB 26), and a response characteristic database (hereinafter referred to as response characteristic DB 28).
  • the relief forming apparatus 22 and the flexographic printing machine 24 constitute an image forming apparatus that forms an image on the printing medium 16.
  • the image processing device 20 creates depth data representing the three-dimensional shape on the main surface of the relief plate 14 based on the input printing data 30 and sends the data to the relief plate forming device 22 side. More specifically, the image processing apparatus 20 includes a RIP (Raster image processor) processing unit 32, a binarization processing unit 34, a threshold data setting unit 36, and a depth data creation unit 38.
  • RIP Raster image processor
  • the RIP processing unit 32 converts the continuous tone image signal 40 expressed in the page description language (PDL) into a raster format image signal (hereinafter also referred to as a raster image signal 42) by performing a rasterization process.
  • PDL page description language
  • the binarization processing unit 34 converts the raster image signal 42 into a binary image signal (hereinafter also referred to as a binary image signal 46) by performing a screening process.
  • This binary image signal 46 represents the shape and arrangement of each halftone dot constituting the image.
  • a systematic dither method using threshold data 44 is employed as a screening processing method.
  • FIG. 2 is a schematic explanatory diagram regarding the operation of the binarization processing unit 34 shown in FIG. This figure shows the concept of binarization using threshold data 44 (threshold matrix).
  • threshold data 44 threshold matrix
  • each matrix element of the threshold data 44 is associated with each address of the raster image signal 42.
  • the binarization processing unit 34 compares the pixel value G of the pixel of interest with the threshold value T of the corresponding matrix element. If G> T is satisfied, “1 (on)” is assigned, and otherwise (G ⁇ T), “0 (off)” is assigned. In this way, the raster image signal 42 with the number of gradation levels “26” is converted into a binary image signal 46 with the number of gradation levels “2”.
  • 13 “on” pixel clusters (AM halftone dots) with hatching are formed for each square area having 25 pixels as a unit.
  • the threshold data setting unit 36 acquires threshold data 44 suitable for the printing data 30 from the threshold DB 26 and supplies it to the binarization processing unit 34 side prior to the screening process described above.
  • the printing data 30 includes various conditions relating to the formation of the printed matter 18 (hereinafter, printing conditions 48).
  • the printing condition 48 includes information on the shape and arrangement of the halftone dots (hereinafter, halftone dot information 50) and information on the apparatus or printing material regarding the image forming system 10 (hereinafter, apparatus / printing material information 52). .
  • FIG. 3 is an explanatory diagram of each parameter related to the halftone dot information 50.
  • the “number of screen lines” is the density of halftone dots per unit length, and LPI (Line Per Inch) is generally used as a unit.
  • the “screen angle” is an angle formed by the arrangement direction of the halftone dots and the vertical axis (or horizontal axis), and a degree is generally used as a unit.
  • the depth data creation unit 38 creates depth data corresponding to the two-dimensional distribution of the engraving amount of the plate material 12 based on the binary image signal 46 acquired from the binarization processing unit 34. To do.
  • the value or range of the depth data can be arbitrarily defined according to the relief forming apparatus 22. Further, instead of the engraving amount, the plate material 12 may be defined as the height from the reference surface.
  • the relief plate forming apparatus 22 can form the relief plate 14 having a plurality of dot projections (not shown) by engraving the main surface of the plate material 12 based on the depth data supplied from the image processing device 20.
  • a CTP (Computer To Plate) drawing apparatus In the case of flexographic printing, as the plate material 12, a highly elastic material including a rubber sheet, a photocurable resin sheet, and the like is used.
  • FIG. 4 is a schematic configuration diagram of a laser engraving machine 60 constituting the relief forming apparatus 22 of FIG.
  • the laser engraving machine 60 includes an exposure head 62, a focus position changing mechanism 64, and an intermittent feeding mechanism 66 in the arrow Y direction (sub-scanning direction).
  • the focus position changing mechanism 64 has a motor 70 and a ball screw 72 that move the exposure head 62 back and forth with respect to the drum 68 to which the printing plate 12 is attached.
  • the focus position changing mechanism 64 moves the focus position under the control of the motor 70.
  • the intermittent feed mechanism 66 includes a ball screw 76 and a sub-scanning motor 78 that rotates the ball screw 76.
  • the intermittent feed mechanism 66 moves the stage 74 on which the exposure head 62 is mounted in the arrow Y direction under the control of the sub-scanning motor 78. As a result, the exposure head 62 is intermittently fed in the direction of the axis 80 of the drum 68.
  • the plate material 12 is chucked by the chuck member 82 on the drum 68.
  • the position of the chuck member 82 is set in a region where exposure by the exposure head 62 is not performed.
  • the main surface of the plate material 12 is laser engraved by irradiating the plate material 12 on the drum 68 with the laser beam 84 from the exposure head 62 while rotating the drum 68 around the axis 80.
  • laser engraving for the next line is performed by performing intermittent feeding in the arrow Y direction.
  • the position of the exposure scan is controlled by sequentially repeating the feeding of the plate 12 in the arrow X direction (main scanning direction) by the rotation of the drum 68 and the intermittent feeding of the exposure head 62 in the arrow Y direction (sub scanning direction).
  • depth data ie, engraving amount
  • the relief forming apparatus 22 forms the relief 14 for printing for each color plate.
  • the letterpress 14 is mounted on a flexographic printing machine 24.
  • the flexographic printing machine 24 is a printing device capable of forming a printed material 18 by transferring ink to the printing plate 16 of the printing medium 16 via the relief plate 14 formed by the relief plate forming device 22.
  • a printing medium having a large surface irregularity shape such as corrugated paper, film, or cloth, can be used as the printing medium 16.
  • FIG. 5 is a schematic configuration diagram of the flexographic printing machine 24 shown in FIG.
  • the flexographic printing press 24 pressurizes the relief plate 14, the plate cylinder 88 to which the relief plate 14 is attached via the cushion tape 86, the anilox roller 92 to which ink is supplied by the doctor chamber 90, and the printing medium 16.
  • the impression cylinder 94 is configured.
  • the ink supplied from the doctor chamber 90 is transferred in the order of the anilox roller 92 and the letterpress plate 14 (the top surface of the dot projection).
  • the printing medium 16 is conveyed in a state of being sandwiched between the plate cylinder 88 and the impression cylinder 94. That is, the ink transferred to the relief plate 14 is further transferred to the printing medium 16, whereby an image composed of a large number of halftone dots is formed on the printing medium 16.
  • FIG. 6 is an electrical block diagram of the image processing apparatus 20 shown in FIG.
  • the image processing apparatus 20 is a computer including a control unit 100, a communication I / F 102, an input unit 104, a display unit 106, an output I / F 108, and a memory 110 (storage medium).
  • the communication I / F 102 is an interface (I / F) that transmits and receives electrical signals from an external device.
  • the image processing apparatus 20 can acquire various data including the printing data 30 from the external apparatus.
  • the input unit 104 includes various input devices such as a mouse, a trackball, a keyboard, or a touch panel.
  • the display unit 106 displays an image corresponding to the display control signal via a control circuit (not shown).
  • a user interface is realized by combining an input function by the input unit 104 and a display function by the display unit 106.
  • the output I / F 108 is an interface (I / F) that transmits depth data as an electrical signal to the relief forming apparatus 22.
  • the memory 110 stores a program and data necessary for the control unit 100 to control each component.
  • printing data 30, a plurality of types of threshold data 44, and response characteristic data 112 described later are stored.
  • the memory 110 may be a non-transitory and computer-readable storage medium.
  • the computer-readable storage medium is a portable medium such as a magneto-optical disk, ROM, CD-ROM, or flash memory, or a storage device such as a hard disk built in the computer system.
  • this storage medium may hold a program dynamically in a short time or may hold a program for a certain period of time.
  • the control unit 100 is configured by a processor such as a CPU (Central Processing Unit).
  • the control unit 100 reads out and executes the program stored in the memory 110, whereby the RIP processing unit 32, the binarization processing unit 34, the threshold data setting unit 36, the depth data creation unit 38 (FIG. 1), and the response characteristics.
  • Each function of the estimation unit 114 (FIG. 14) can be realized.
  • the threshold data setting unit 36 is a threshold candidate acquisition unit 116 that acquires threshold data 44 (hereinafter referred to as threshold candidates), which is a candidate for setting, and a spatial frequency response characteristic (hereinafter simply referred to as “response characteristic”) that represents a change in halftone dot shape.
  • threshold candidates threshold data 44
  • response characteristic a spatial frequency response characteristic
  • a response characteristic acquisition unit 118 that acquires a threshold value
  • a threshold value candidate evaluation unit 120 that quantitatively evaluates the reproducibility of the halftone dot shape of each threshold value candidate.
  • the threshold candidate evaluation unit 120 includes a binary halftone dot creation unit 122, a simulated image creation unit 124, a quantitative value calculation unit 126, and a threshold data determination unit 128. The function of each part will be described later.
  • the image processing apparatus 20 as the halftone data setting apparatus according to this embodiment is configured as described above. Next, the operation of the image processing apparatus 20 (mainly the threshold data setting unit 36) shown in FIGS. 1 and 6 will be described in detail with reference mainly to the flowchart of FIG.
  • step S ⁇ b> 1 the threshold data setting unit 36 acquires the printing condition 48 associated with the continuous tone image signal 40.
  • the printing condition 48 includes halftone dot information 50 and device / printing material information 52.
  • Examples of the halftone dot information 50 include the screen line number, screen angle, halftone dot target shape (circular, elliptical, rectangular, and combinations thereof), or screen characteristics (AM screen or FM screen).
  • the device / printing material information 52 includes, for example, the device configuration of the image forming system 10, the printing medium 16, the plate material 12, the color material, or the component of the device (for example, the exposure head 62 in FIG. 4). .
  • the threshold candidate acquisition unit 116 acquires threshold data 44 (hereinafter also referred to as a threshold candidate) that partially matches the printing condition 48 obtained in step S1.
  • the threshold DB 26 executes a database search process in response to a request from the image processing apparatus 20 and supplies threshold data 44 matching the halftone dot information 50 to the image processing apparatus 20 side.
  • the threshold candidate acquisition unit 116 reads and acquires all threshold candidates temporarily stored in the memory 110.
  • step S3 the response characteristic acquisition unit 118 acquires response characteristic data 112 that partially matches the printing condition 48 obtained in step S1.
  • the response characteristic DB 28 executes a database search process in response to a request from the image processing apparatus 20 and supplies response characteristic data 112 matching the apparatus / printing material information 52 to the image processing apparatus 20 side. . Thereafter, the response characteristic acquisition unit 118 reads and acquires the response characteristic data 112 temporarily stored in the memory 110.
  • the response characteristic data 112 is data indicating a change in halftone dot shape that occurs in the image forming process from the creation of the binary image signal 46 (FIG. 2) to the formation of the image. Note that it depends on.
  • the response characteristic acquisition unit 118 uses the “engraving response characteristic” corresponding to the optical characteristic of the laser engraving machine 60 (FIG. 4) as the response characteristic data 112, and the transition of the flexographic printing machine 24 (FIG. 5). “Transition response characteristics” corresponding to the characteristics are acquired.
  • FIG. 8A and 8B are graphs showing an example of engraving response characteristics.
  • FIG. 8A shows a point spread function ⁇ hereinafter, PSF (Point (Spread Function) ⁇
  • FIG. 8B corresponds to a cross-sectional view of the Fourier transform function of the PSF shown in FIG. 8A.
  • the X axis and the Y axis in this graph correspond to the arrow X direction and the arrow Y direction shown in FIG.
  • PSF is a function having a substantially conical shape with the origin O as a peak.
  • PSF is a delta function in an ideal system, but in an actual system, a point image is expanded by various optical factors.
  • the extent of expansion in the X-axis direction depends on the effective beam diameter in the arrow X direction of the laser beam 84 (FIG. 4), the rotational speed of the drum 68 (FIG. 4), the rise and / or fall characteristics of a laser light source (not shown), and the like.
  • the extent of expansion in the Y-axis direction depends on the effective beam diameter of the laser beam 84 in the arrow Y direction, the exposure multiplicity of the laser beam 84, and the like.
  • the MTF (Modulation Transfer Function) shown in FIG. 8B corresponds to the sinusoidal response characteristics for each spatial frequency. As can be seen from the figure, the MTF gradually decreases as the spatial frequency increases. Then, the “sub-scan” MTF indicated by the solid line tends to be higher than the “main scan” indicated by the broken line.
  • the threshold candidate evaluation unit 120 creates a test image signal used for quantitative evaluation of each threshold candidate.
  • An arbitrary raster image signal 42 may be used as the test image signal.
  • an image signal including at least one flat halftone area (halftone solid area) is created.
  • the threshold candidate evaluation unit 120 designates one threshold candidate that has not yet been evaluated among the plurality of threshold candidates (three threshold candidates A to C) acquired in step S2.
  • the designated threshold candidate is referred to as a “designated threshold candidate”.
  • step S6 the binary halftone dot creating unit 122 performs a screening process using the designated threshold candidate on each test image signal generated in step S4, thereby obtaining a binary halftone dot signal for each gradation level.
  • FIG. 9 is an explanatory diagram showing the result of binarization processing performed on the test image signal.
  • five types of flat halftone images having different gray levels are represented from test image signals including five types of flat halftone regions having gradation levels of 10%, 30%, 50%, 70%, and 90%. Create a binary dot signal.
  • step S7 the simulated image creating unit 124 creates a simulated image signal in which the halftone dot shape on the printed matter 18 is simulated from the binary halftone dot signal created in step S6.
  • FIG. 10 is a schematic explanatory diagram regarding the operation of the simulated image creation unit 124 shown in FIG. This figure schematically shows the processing steps from the creation of a simulated image signal from a binary halftone signal.
  • the simulated image creating unit 124 converts the binary halftone dot signal indicating the on or off state of each pixel into a multi-value halftone dot signal corresponding to the height distribution of the relief plate 14. Thereafter, the simulated image creating unit 124 executes a process for applying a response characteristic in the spatial frequency domain to the obtained multi-value halftone dot signal.
  • a mask process, a filter process, or a combination thereof can be applied to this action process. This is because the calculation results in Euclidean space or Fourier space are substantially equivalent due to the nature of the mapping of Fourier transform.
  • mask processing means image processing for modulating a spatial frequency component (spectral intensity) of an image by executing a mask operation (convolution) on the Euclidean space.
  • Filter processing means image processing for modulating a spatial frequency component (spectral intensity) of an image signal by executing a filter operation (multiplication) in Fourier space.
  • the simulated image creating unit 124 sequentially performs a first mask process using a mask showing engraving response characteristics and a second mask process using a mask showing transfer response characteristics. Create a simulated image signal. Further, according to the processing flow in the lower part of the figure, the simulated image creating unit 124 sequentially performs the first filter process using the filter showing the engraving response characteristic and the second filter process using the filter showing the transfer response characteristic. To create a simulated image signal.
  • FFT Fast Fourier Transform
  • IFFT Inverse Fast Fourier Transform
  • step S8 the quantitative value calculation unit 126 calculates an evaluation value obtained by executing a predetermined evaluation process on the simulated image signal created in step S7 as a quantitative value indicating the degree of image quality.
  • the power spectrum of the simulated image signal is used to calculate the quantitative value. For example, in order to increase the degree of correlation with the visibility, a preprocessing for multiplying the power spectrum by a human standard visual response characteristic may be performed.
  • FIG. 11 is a graph of the Dooley-Shaw function with an observation distance of 300 mm.
  • This function is a kind of VTF (Visual Transfer Function), and is a typical function imitating human standard visual response characteristics. Specifically, this corresponds to the square value of the contrast ratio characteristic of luminance.
  • the horizontal axis of the graph is the spatial frequency (unit: Cycle / mm), and the vertical axis is the VTF value (unit is dimensionless).
  • the function shape of the VTF is not limited to this, and various characteristics derived from a mathematical model, experimental data, or the like may be applied. Further, the observation distance is not limited to 300 mm, but may be variously changed according to the observation mode of the image or the evaluation standard.
  • the pixel value of the simulated image signal is set to an amount highly correlated (preferably linear) with the amount of light that reflects or transmits the image. You may convert it. For example, RGB values, tristimulus values (XYZ), light reflectance in the case of a reflective document, or light transmittance in the case of a transmissive document may be used.
  • FIG. 12A is a graph showing an example of a power spectrum in a simulated image signal.
  • the horizontal axis of the graph is the spatial frequency (unit: Cycle / mm), and the vertical axis of the graph is the spectrum intensity (unit is arbitrary; hereinafter, also simply referred to as “component”).
  • the power spectrum shown in this figure has four peaks at spatial frequencies of about 1.9, about 2.4, about 3.6, and about 7.4 Cycle / mm.
  • the quantitative value calculation unit 126 removes a component corresponding to each halftone dot (here, a component around about 7.4 Cycle / mm), and then the maximum value of the spectral intensity (corresponding to the main component; hereinafter referred to as the maximum intensity). ).
  • the quantitative value calculation unit 126 extracts “60”, which is the peak maximum value at about 2.4 Cycle / mm, as the maximum intensity. It is evaluated that the larger the maximum intensity value is, the lower the image quality is, and the smaller the value is, the higher the image quality is.
  • This quantitative value is a value that correlates quantitatively or qualitatively with the image quality of the printed matter 18. This is because components other than the component corresponding to each halftone dot contribute as image noise including color unevenness, moire, graininess, and the like.
  • the graph in this example corresponds to the power spectrum of a simulated image signal created from a test image signal having a gradation level of 70% using the threshold value candidate B.
  • the quantitative value calculation unit 126 extracts the maximum intensity for each simulated image signal for each gradation level, including 70%.
  • the quantitative value calculation unit 126 calculates the quantitative value of the designated threshold value candidate based on the maximum intensity for each gradation level obtained as described above.
  • the quantitative value calculation unit 126 calculates an average value and / or maximum value over all gradation levels as a “quantitative value”. Note that the calculation formula for the quantitative value is not limited to this, and various values may be calculated based on the statistical value for the maximum intensity in each flat network region or a combination thereof (for example, the weighted addition value of the average value and the maximum value). Good.
  • the evaluation process for the simulated image signal is not limited to the above-described method, and includes known methods including granularity ⁇ RMS (Root (Mean Square) ⁇ , sharpness (MTF, etc.), or gradation (gradation smoothness, etc.). Needless to say, various image evaluation methods can be applied.
  • RMS Root (Mean Square) ⁇
  • MTF sharpness
  • gradation gradation smoothness, etc.
  • step S9 the threshold candidate evaluation unit 120 determines whether or not the designation for all threshold candidates is completed. If it is determined that it has not been completed (step S9: NO), the process returns to step S5, and steps S5 to S8 are sequentially repeated. On the other hand, if it is determined that all have been completed (step S9: YES), the process proceeds to the next step (S10).
  • step S10 the threshold data determination unit 128 determines one threshold data 44 by comparing and / or evaluating threshold candidates.
  • FIG. 12B is a graph showing the maximum intensity for each gradation level in each of the threshold candidates A to C.
  • FIG. 13 is a table showing calculation results of quantitative values for the respective threshold candidates A to C.
  • the threshold candidate evaluation unit 120 determines the threshold candidate A having the smallest “average value” or “maximum value” as the threshold data 44 based on the evaluation result shown in FIG.
  • step S11 the threshold data setting unit 36 performs a setting for associating the threshold data 44 determined in step S10 with the printing data 30. Thereby, when creating the depth data from the printing data 30, the threshold data 44 suitable for the image forming system 10 can be automatically set. In this way, the operation of the image processing apparatus 20 (mainly the threshold data setting unit 36) is completed.
  • FIG. 14 is an overall configuration diagram of the response characteristic estimation system 150.
  • the response characteristic estimation system 150 sequentially forms a sample letterpress 154 and a sample print 156 based on the sample halftone dot signal 152, and uses the measurement results obtained by the first measuring instrument 158 and the second measuring instrument 160 to This is a system for estimating response characteristics in whole or in part.
  • the sample halftone dot signal 152 is a binary image signal representing the shape and arrangement of each halftone dot constituting an arbitrary image.
  • the sample halftone dot signal 152 may be the same signal as the binary halftone dot signal (FIG. 9) in order to facilitate estimation processing described later.
  • the depth data creation unit 38 creates depth data corresponding to the two-dimensional distribution of the engraving amount of the plate material 12 based on the input sample halftone signal 152. Then, the letterpress forming apparatus 22 engraves the main surface of the plate material 12 based on the depth data supplied from the depth data creation unit 38, thereby providing a sample letterpress having a plurality of dot projections (not shown). 154 is formed. Then, the flexographic printing machine 24 forms a sample printed matter 156 by transferring ink to the plate surface of the printing medium 16 through the sample relief plate 154 formed by the relief plate forming apparatus 22.
  • the first measuring device 158 is a device that optically reads the uneven shape of the sample relief 154, and can measure the dot area ratio on the sample relief 154.
  • the second measuring device 160 is a device that optically reads the color characteristics of the sample printed matter 156, and can measure the dot area ratio on the sample printed matter 156.
  • the response characteristic estimation unit 114 estimates a response characteristic according to an event that changes the halftone dot shape in the image forming process. One or two or more events depending on the recording method of the image forming apparatus may be considered as an event that changes the halftone dot shape.
  • the response characteristic estimation unit 114 performs “engraving response characteristics” when forming the relief plate 14 by engraving the main surface of the plate material 12, and when transferring ink to the printing medium 16 via the relief plate 14. Estimate the “metastasis response characteristics” respectively.
  • the response characteristic estimation unit 114 is measured by the engraving characteristic estimation unit 162 that estimates the engraving response characteristic based on the halftone dot area ratio on the sample letterpress 154 measured by the first measuring unit 158 and the second measuring unit 160.
  • a transfer characteristic estimation unit 164 that estimates a transfer response characteristic based on the dot area ratio on the sample print 156;
  • FIG. 15A is a schematic front view of the sample halftone dot signal 152 shown in FIG.
  • the filled regions in FIG. 15A and FIGS. 15B and 15C described later indicate portions that contribute to the determination of each halftone dot shape.
  • each halftone dot is circular and the diameter is D (known value).
  • FIG. 15B is a schematic front view of the sample relief 154 formed based on the sample halftone dot signal 152 of FIG. 15A.
  • the halftone dot area ratio on the sample relief 154 is smaller than the design value of the halftone dot area ratio.
  • the average diameter corresponding to the measured dot area ratio is (Dd1).
  • FIG. 15C is a schematic front view of a specimen print 156 formed based on the specimen relief 154 of FIG. 15B.
  • the plate material 12 When a highly elastic material is used as the plate material 12, the dot gain tends to increase when the ink is transferred.
  • the halftone dot area ratio on the sample printed matter 156 is larger than the halftone dot area ratio on the sample letterpress 154.
  • the average diameter corresponding to the measured dot area ratio is (D ⁇ d2) (where d2 ⁇ d1).
  • FIG. 16A is a first explanatory diagram relating to an estimated model of the height shape of a halftone dot projection.
  • the horizontal axis of the graph is the position in the X-axis (or Y-axis) direction, and the vertical axis of the graph is the position in the Z-axis (height) direction.
  • the solid line in this figure indicates the ideal height shape of each halftone protrusion in the sample halftone dot signal 152 (hereinafter, ideal halftone dot shape 170).
  • the broken line in this figure indicates the virtual height shape (hereinafter, virtual dot shape 172) of each dot projection in the sample relief 154.
  • This estimation model is constructed on the premise that the height of the halftone dot projections is always constant, and mainly simulates the phenomenon of dot gain.
  • FIG. 16B is a graph showing the Fourier transform characteristics of the height shape shown in FIG. 16A.
  • the horizontal axis of the graph is the spatial frequency
  • the vertical axis of the graph is the Fourier transform value.
  • the solid line in this figure indicates the Fourier transform value (hereinafter, ideal frequency component 170f) for the ideal halftone dot shape 170.
  • the broken line in this figure indicates the Fourier transform value (hereinafter, virtual frequency component 172f) for the virtual halftone dot shape 172.
  • the engraving characteristic estimation unit 162 estimates the ratio of the virtual frequency component 172f to the ideal frequency component 170f, that is, the value of (virtual frequency component 172f / ideal frequency component 170f) for each spatial frequency as the engraving response characteristic.
  • FIG. 17A is a second explanatory diagram relating to an estimated model of the height shape of the halftone dot projections. Since the definition of the graph is the same as that in FIG. 16A, the description thereof is omitted.
  • the solid line in this figure indicates the ideal halftone dot shape 170, and the broken line in this figure indicates the virtual halftone dot shape 174.
  • This estimation model is built on the assumption that the volume (area) of the dot projection is always preserved regardless of the diameter of the top. As can be understood from this figure, the inclination of the side surface becomes smaller as the diameter of the top portion becomes smaller.
  • FIG. 17B is a graph showing the Fourier transform characteristics of the height shape shown in FIG. 17A.
  • the definition of the graph is the same as that in FIG.
  • the solid line in this figure indicates the ideal frequency component 170f, and the broken line in this figure indicates the virtual frequency component 174f.
  • the engraving characteristic estimation unit 162 estimates the ratio of the virtual frequency component 174f to the ideal frequency component 170f, that is, the value of (virtual frequency component 174f / ideal frequency component 170f) for each spatial frequency as the engraving response characteristic.
  • the transfer characteristic estimation unit 164 estimates the transfer response characteristic based on the response characteristic multiplication rule (cascade model). First, the transfer characteristic estimation unit 164 estimates response characteristics (hereinafter referred to as overall response characteristics) in the entire image forming system 10 using the same or different estimation models as in FIGS. 16A to 17B. Then, the transfer characteristic estimation unit 164 estimates the transfer response characteristic by dividing the engraving response characteristic from the overall response characteristic.
  • response characteristic multiplication rule cascade model
  • the response characteristic estimation unit 114 supplies the estimated engraving response characteristic and the transfer response characteristic to the response characteristic DB 28 side.
  • the response characteristic DB 28 associates the acquired engraving response characteristic with the types of the plate material 12 and the relief forming apparatus 22.
  • the response characteristic DB 28 associates the acquired transfer response characteristic with the types of the plate material 12 and the flexographic printing machine 24.
  • the response characteristic data 112 is stored and managed in the response characteristic DB 28 in a state associated with the device / printing material information 52.
  • the image processing apparatus 20 as a threshold data setting apparatus is an apparatus for setting threshold data 44 for creating a binary image signal 46 representing the shape and arrangement of each halftone dot constituting an image based on the continuous tone image signal 40. It is.
  • threshold value acquisition for acquiring halftone dot information 50 relating to the shape and arrangement of each halftone dot associated with the continuous tone image signal 40 and acquiring a plurality of threshold data 44 matching the halftone dot information 50 as each threshold candidate.
  • a response characteristic acquisition unit 118 that acquires, as response characteristic data 112 in the spatial frequency domain, a change in halftone dot shape that occurs in the image formation process from the creation of the binary image signal 46 to the formation of the image.
  • a simulated image creating unit 124 that creates a simulated image signal in which the change in the halftone dot shape is simulated is executed, and a predetermined evaluation process is executed on the simulated image signal.
  • a threshold data determination unit 128 that determines the threshold data 44 to be set from the threshold candidates.
  • the change in the halftone dot shape that occurs in the image forming process from the creation of the binary image signal 46 to the formation of the image is acquired as the response characteristic data 112 in the spatial frequency domain, and the response characteristic data 112
  • the halftone dot shape can be easily and accurately reproduced on the simulated image.
  • a plurality of threshold data 44 acquired as each threshold candidate is compared and / or evaluated based on a quantitative value indicating the degree of image quality, even if the number of types of printing apparatuses and printing materials is enormous, The threshold data 44 suitable for these combinations can be determined easily and easily.
  • flexographic printing has been described as an example.
  • the image forming apparatus is not limited to this printing method, and can be applied to offset printing, gravure printing, a digital printing machine (on-demand printing machine), and the like.
  • Quantitative value calculation unit 128 ... Threshold data determination unit 150 ... Response characteristic estimation system 152 . Sample halftone dot signal 154 ... Sample Letterpress 156 ... Printed specimen 158 ... First measuring device 160 ... Second measuring device 170 ... Ideal halftone dot shape 172, 174 ... Virtual halftone dot shape

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Signal Processing (AREA)
  • Multimedia (AREA)
  • Manufacturing & Machinery (AREA)
  • Optics & Photonics (AREA)
  • Plasma & Fusion (AREA)
  • Manufacture Or Reproduction Of Printing Formes (AREA)
  • Facsimile Image Signal Circuits (AREA)
  • Color, Gradation (AREA)
  • Image Processing (AREA)

Abstract

印刷装置及び印刷材の種類数が膨大な場合であっても、これらの組み合わせに適した閾値データを簡便且つ容易に決定可能な閾値データ設定装置、方法及びプログラム、並びに画像形成システムを提供する。 ２値画像信号４６の作成後から画像を形成するまでの画像形成過程にて生じる網点形状の変化を、空間周波数領域での応答特性データ１１２として取得する。応答特性データ１１２を作用することにより、網点形状の変化が模擬的に再現された模擬画像信号を作成する。模擬画像信号に対して所定の評価処理を実行することによりそれぞれ算出された定量値に基づく比較及び/又は評価によって、各閾値候補の中から設定対象である閾値データ４４を決定する。

Description

閾値データ設定装置、方法及びプログラム、並びに画像形成システム
　この発明は、画像を構成する各網点の形状及び配置を表す２値画像信号を、連続調画像信号に基づき作成するための閾値データを設定する閾値データ設定装置、方法及びプログラム、並びに画像形成システムに関する。
　印刷分野において、中間成果物を生成せずに電子データにより直接刷版を形成するＣＴＰ（Computer to Plate）が普及している。通常、刷版を形成する過程にて、連続調である画像信号に対して、組織的ディザ法を含むスクリーニング処理が実行される。従来から、高画質の印刷物を得るために、印刷システムの全体又は一部の特性を考慮した閾値データを決定する手法が種々提案されている。
　特許文献１では、（１）ドットの擬似発生、及び（２）該ドットを含む濃度パターンの特定周波数成分の評価、を順次行うことにより、モアレ等の発生が起こり難い閾値配列を決定する方法が提案されている。
特開２００２－３６８９９５号公報
　近時、使用可能な印刷装置又は印刷材（色材、印刷媒体、刷版等）の種類が多種多様になっている。例えば、印刷装置の種類として、オフセット印刷機のみならず、グラビア印刷機、フレキソ印刷機、デジタル印刷機（オンデマンド印刷機）等が存在する。また、印刷媒体の種類として、合成紙、厚紙、若しくはアルミ蒸着紙等の紙類のみならず、塩化ビニル・ＰＥＴ（ポリエチレンテレフタラート）等の樹脂、ターポリン、又は金属シート等が存在する。
　ところで、一般的には、使用される印刷装置及び印刷材の組み合わせに応じて、印刷特性が変化することが知られている。例えば、オフセット印刷及びフレキソ印刷では、版材の種類（特に、弾性特性）が異なるため、ドットゲインが変化する程度が異なる傾向がみられる。また、刷版が不要なデジタル印刷機の場合、当該刷版の形成に起因する画質変化を何ら考慮する必要はない。
　しかし、特許文献１に記載の方法によれば、印刷装置及び印刷材の組み合わせ毎に適した閾値データを作成することを前提としているため、その準備及び作成処理に多くの作業工数を必要とし、改良の余地が十分にあった。
　本発明は、上記した課題を解決するためになされたものであり、印刷装置及び印刷材の種類数が膨大な場合であっても、これらの組み合わせに適した閾値データを簡便且つ容易に決定可能な閾値データ設定装置、方法及びプログラム、並びに画像形成システムを提供することを目的とする。
　本発明に係る閾値データ設定装置は、画像を構成する各網点の形状及び配置を表す２値画像信号を、連続調画像信号に基づき作成するための閾値データを設定する装置であって、上記連続調画像信号に紐付けられた上記各網点の形状及び配置に関する網点情報を取得し、
該網点情報に合致する複数の閾値データを各閾値候補として取得する閾値候補取得部と、上記２値画像信号の作成後から上記画像を形成するまでの画像形成過程にて生じる網点形状の変化を、空間周波数領域での応答特性として取得する応答特性取得部と、上記応答特性取得部により取得された上記応答特性を作用することにより、上記網点形状の変化が模擬的に再現された模擬画像信号を作成する模擬画像作成部と、上記模擬画像作成部により作成された上記模擬画像信号に対して所定の評価処理を実行して得られる評価値を、画質の程度を示す定量値として算出する定量値算出部と、上記閾値候補取得部により取得された上記各閾値候補を適用し、上記定量値算出部によりそれぞれ算出された上記定量値に基づく比較及び/又は評価によって、上記各閾値候補の中から設定対象である閾値データを決定する閾値データ決定部とを備える。
　このように、２値画像信号の作成後から画像を形成するまでの画像形成過程にて生じる網点形状の変化を、空間周波数領域での応答特性として取得し、当該応答特性を作用することにより、網点形状の変化が模擬的に再現された模擬画像信号を作成するので、網点形状を模擬画像上に簡便に且つ精度よく再現可能である。そして、各閾値候補として取得された複数の閾値データを、画質の程度を示す定量値に基づき比較及び/又は評価するので、印刷装置及び印刷材の種類数が膨大な場合であっても、これらの組み合わせに適した閾値データを簡便且つ容易に決定できる。
　また、上記模擬画像作成部は、平網領域を少なくとも１つ含む上記模擬画像信号を作成し、上記定量値算出部は、上記模擬画像信号の上記平網領域でのパワースペクトルに関し、上記各網点に対応する成分を除いた主成分に基づいて上記定量値を算出することが好ましい。
　また、上記模擬画像作成部は、階調レベルが異なる平網領域を２つ以上含む上記模擬画像信号を作成し、上記定量値算出部は、各上記平網領域における上記主成分についての統計値に基づいて上記定量値を算出することが好ましい。
　また、上記画像形成過程にて上記網点形状に変化を及ぼす事象に応じた上記応答特性を推定する応答特性推定部を更に備え、上記応答特性取得部は、上記応答特性推定部により推定された上記応答特性を取得することが好ましい。
　また、上記応答特性推定部は、刷版を介して被印刷材に色材を転移させる際の上記応答特性を推定することが好ましい。
　また、上記応答特性推定部は、版材の主面を彫刻することにより上記刷版を形成する際の上記応答特性を推定することが好ましい。
　また、上記閾値データ決定部により決定された上記閾値データを設定することにより、上記連続調画像信号を上記２値画像信号に変換する２値化処理部を更に備えることが好ましい。
　本発明に係る画像形成システムは、上記の閾値データ設定装置と、上記閾値データ設定装置により変換された上記２値画像信号に基づいて上記画像を形成する画像形成装置とを備える。
　本発明に係る閾値データ設定方法は、画像を構成する各網点の形状及び配置を表す２値画像信号を、連続調画像信号に基づき作成するための閾値データを設定する方法であって、上記連続調画像信号に紐付けられた上記各網点の形状及び配置に関する網点情報を取得し、該網点情報に合致する複数の閾値データを各閾値候補として取得するステップと、上記２値画像信号の作成後から上記画像を形成するまでの画像形成過程にて生じる網点形状の変化を、空間周波数領域での応答特性として取得するステップと、取得された上記応答特性を作用することにより、上記網点形状の変化が模擬的に再現された模擬画像信号を作成するステップと、作成された上記模擬画像信号に対して所定の評価処理を実行して得られる評価値を、画質の程度を示す定量値として算出するステップと、取得された上記各閾値候補を適用し、それぞれ算出された上記定量値に基づく比較及び/又は評価によって、上記各閾値候補の中から設定対象である閾値データを決定するステップとをコンピュータに実行させる。
　本発明に係る閾値データ設定プログラムは、画像を構成する各網点の形状及び配置を表す２値画像信号を、連続調画像信号に基づき作成するための閾値データを設定するためのプログラムであって、上記連続調画像信号に紐付けられた上記各網点の形状及び配置に関する網点情報を取得し、該網点情報に合致する複数の閾値データを各閾値候補として取得するステップと、上記２値画像信号の作成後から上記画像を形成するまでの画像形成過程にて生じる網点形状の変化を、空間周波数領域での応答特性として取得するステップと、取得された上記応答特性を作用することにより、上記網点形状の変化が模擬的に再現された模擬画像信号を作成するステップと、作成された上記模擬画像信号に対して所定の評価処理を実行して得られる評価値を、画質の程度を示す定量値として算出するステップと、取得された上記各閾値候補を適用し、それぞれ算出された上記定量値に基づく比較及び/又は評価によって、上記各閾値候補の中から設定対象である閾値データを決定するステップとをコンピュータに実行させる。
　本発明に係る閾値データ設定装置、方法及びプログラム、並びに画像形成システムによれば、２値画像信号の作成後から画像を形成するまでの画像形成過程にて生じる網点形状の変化を、空間周波数領域での応答特性として取得し、当該応答特性を作用することにより、網点形状の変化が模擬的に再現された模擬画像信号を作成するので、網点形状を模擬画像上に簡便に且つ精度よく再現可能である。そして、各閾値候補として取得された複数の閾値データを、画質の程度を示す定量値に基づき比較及び/又は評価するので、印刷装置及び印刷材の種類数が膨大な場合であっても、これらの組み合わせに適した閾値データを簡便且つ容易に決定できる。
この実施形態に係る画像形成システムの全体構成図である。 図１に示す２値化処理部の動作に関する概略説明図である。 網点情報に関する各パラメータの説明図である。 図１に示す凸版形成装置を構成するレーザ彫刻機の概略構成図である。 図１に示すフレキソ印刷機の概略構成図である。 図１に示す画像処理装置の電気的なブロック図である。 図１及び図６に示す画像処理装置の動作説明に供されるフローチャートである。 図８Ａ及び図８Ｂは、彫刻応答特性の一例を示すグラフである。 各テスト画像信号に対して２値化処理を施した結果を示す説明図である。 図６に示す模擬画像作成部の動作に関する概略説明図である。 観察距離が３００ｍｍであるドゥーリー・ショー関数のグラフである。 図１２Ａは、模擬画像信号におけるパワースペクトルの一例を示すグラフである。図１２Ｂは、各閾値候補におけるパワースペクトルのピーク特性を示すグラフである。 各閾値候補に対する定量結果を示すテーブルである。 応答特性推定システムの全体構成図である。 図１５Ａは、標本画像網点を可視化した概略正面図である。図１５Ｂは、図１５Ａの標本画像網点に基づき形成された標本凸版の正面図である。図１５Ｃは、図１５Ｂの標本凸版に基づき形成された標本印刷物の正面図である。 図１６Ａは、網点突起部の高さ形状の推定モデルに関する第１の説明図である。図１６Ｂは、図１６Ａに示す高さ形状のフーリエ変換特性を示すグラフである。 図１７Ａは、網点突起部の高さ形状の推定モデルに関する第２の説明図である。図１７Ｂは、図１７Ａに示す高さ形状のフーリエ変換特性を示すグラフである。
　以下、本発明に係る閾値データ設定方法について、それを実施する閾値データ設定装置、プログラム及び画像形成システムとの関係において好適な実施形態を挙げ、添付の図面を参照しながら詳細に説明する。本明細書において、画像を形成することを「印刷」又は「印字」という場合がある。
［画像形成システム１０の全体構成］
　図１は、この実施形態に係る画像形成システム１０の全体構成図である。この画像形成システム１０は、平板状の版材１２を彫刻することにより凸版１４（刷版）を形成し、該凸版１４を介して被印刷体１６にインキ（色材）を転移させることにより印刷物１８を生産するシステムである。この印刷物１８は、被印刷体１６上に少なくとも１色の網点が形成されてなる。
　画像形成システム１０は、閾値データ設定装置としての画像処理装置２０と、凸版形成装置２２と、フレキソ印刷機２４と、閾値データベース（以下、閾値ＤＢ２６）と、応答特性データベース（以下、応答特性ＤＢ２８）を基本的に備えている。この実施形態にて、凸版形成装置２２及びフレキソ印刷機２４は、被印刷体１６上に画像を形成する画像形成装置を構成する。
　画像処理装置２０は、入力された印刷用データ３０に基づいて、凸版１４の主面上の三次元形状を表す深さデータを作成し、当該データを凸版形成装置２２側に送出する。画像処理装置２０は、より詳細には、ＲＩＰ（Raster Image Processor）処理部３２、２値化処理部３４、閾値データ設定部３６及び深さデータ作成部３８を備えている。
　ＲＩＰ処理部３２は、ページ記述言語（ＰＤＬ）で表現された連続調画像信号４０に対してラスタライズ処理を施すことにより、ラスタ形式の画像信号（以下、ラスタ画像信号４２ともいう）に変換する。
　２値化処理部３４は、ラスタ画像信号４２に対してスクリーニング処理を施すことにより、２値の画像信号（以下、２値画像信号４６ともいう）に変換する。この２値画像信号４６は、画像を構成する各網点の形状及び配置を表している。この実施形態では、スクリーニング処理手法として、閾値データ４４による組織的ディザ法を採用する。
　図２は、図１に示す２値化処理部３４の動作に関する概略説明図である。本図は、閾値データ４４（閾値マトリクス）を用いた２値化の概念を示す。先ず、ラスタ画像信号４２の各アドレスに閾値データ４４の各行列要素を対応付ける。そして、２値化処理部３４は、着目する画素での画素値Ｇと、対応する行列要素での閾値Ｔの大小関係をそれぞれ比較する。Ｇ＞Ｔを満たす場合には「１（オン）」を割り当て、それ以外（Ｇ≦Ｔ）の場合には「０（オフ）」を割り当てる。このようにして、階調レベルの数が「２６」であるラスタ画像信号４２は、階調レベルの数が「２」である２値画像信号４６に変換される。本図例では、２５個の画素を単位とする正方領域毎に、ハッチングが付された１３個の「オン」画素のクラスタ（ＡＭ網点）が形成されている。
　図１に戻って、閾値データ設定部３６は、上記したスクリーニング処理に先立ち、閾値ＤＢ２６から印刷用データ３０に適した閾値データ４４を取得し、２値化処理部３４側に供給する。この印刷用データ３０には、上記した連続調画像信号４０の他、印刷物１８の形成に関する各種条件（以下、印刷条件４８）が含まれている。そして、印刷条件４８は、網点の形状及び配置に関する情報（以下、網点情報５０）と、画像形成システム１０に関する装置又は印刷材に関する情報（以下、装置・印刷材情報５２）で構成される。
　図３は、網点情報５０に関する各パラメータの説明図である。「スクリーン線数」は単位長さ当たりの網点の密度であり、単位としてＬＰＩ（Line Per Inch）が一般に用いられる。「スクリーン角度」は網点の配列方向と垂直軸（又は水平軸）とのなす角であり、単位として度（degree）が一般に用いられる。
　図１に戻って、深さデータ作成部３８は、２値化処理部３４から取得した２値画像信号４６に基づいて、版材１２の彫刻量の２次元分布に相当する深さデータを作成する。深さデータの値又は範囲は、凸版形成装置２２に応じて任意に定義することができる。また、彫刻量に代わって、版材１２の基準面からの高さとして定義してもよい。
［凸版形成装置２２の概略構成］
　凸版形成装置２２は、画像処理装置２０から供給された深さデータに基づき、版材１２の主面を彫刻することにより、複数の網点突起部（図示しない）を備える凸版１４を形成可能なＣＴＰ（Computer To Plate）描画装置である。フレキソ印刷の場合は、版材１２として、ゴムシート、光硬化性樹脂シート等を含む高弾性材料が用いられる。
　図４は、図１の凸版形成装置２２を構成するレーザ彫刻機６０の概略構成図である。レーザ彫刻機６０は、露光ヘッド６２と、ピント位置変更機構６４と、矢印Ｙ方向（副走査方向）への間欠送り機構６６とを備える。
　ピント位置変更機構６４は、版材１２が取り付けられたドラム６８に対して、露光ヘッド６２を前後移動させるモータ７０及びボールネジ７２を有する。ピント位置変更機構６４は、モータ７０の制御によりピント位置を移動させる。
　間欠送り機構６６は、ボールネジ７６及び該ボールネジ７６を回転させる副走査モータ７８を有する。間欠送り機構６６は、副走査モータ７８の制御により、露光ヘッド６２が搭載されたステージ７４を矢印Ｙ方向に移動させる。これにより、露光ヘッド６２は、ドラム６８の軸線８０の方向に間欠送りされる。
　版材１２は、ドラム６８上でチャック部材８２によりチャックされる。チャック部材８２の位置は、露光ヘッド６２による露光が行われない領域に設定されている。この場合、軸線８０を中心としてドラム６８を回転させながら、該ドラム６８上の版材１２に対して、露光ヘッド６２からレーザビーム８４を照射することにより、版材１２の主面がレーザ彫刻される。そして、ドラム６８の回転により、露光ヘッド６２の前方をチャック部材８２が通過するときに、矢印Ｙ方向に間欠送りを行うことにより、次のライン分のレーザ彫刻が行われる。
　ドラム６８の回転による版材１２の矢印Ｘ方向（主走査方向）への送りと、露光ヘッド６２の矢印Ｙ方向（副走査方向）の間欠送りを逐次繰り返すことにより、露光走査の位置が制御される。この制御と同期して、画像処理装置２０（図１）から供給された位置毎の深さデータ（すなわち、彫刻量）に基づいて、レーザビーム８４のオン・オフ又は発光強度が制御される。この結果、版材１２の主面（版面）に所望の形状を有するレリーフ（網点突起部を含む）が形成される。
　このようにして、凸版形成装置２２は、印刷に供される凸版１４を色版毎に形成する。この凸版１４は、フレキソ印刷機２４に装着される。
［フレキソ印刷機２４の概略構成］
　フレキソ印刷機２４は、凸版形成装置２２により形成された凸版１４を介し、被印刷体１６の版面にインキを転移することにより、印刷物１８を形成可能な印刷装置である。フレキソ印刷の場合は、被印刷体１６として、印刷用途の枚葉紙、ロール紙のみならず、表面の凹凸形状が大きい印刷媒体、例えば、段ボール紙、フイルム、又は布等を用いることができる。
　図５は、図１に示すフレキソ印刷機２４の概略構成図である。フレキソ印刷機２４は、上記した凸版１４と、クッションテープ８６を介して凸版１４が取り付けられる版胴８８と、ドクターチャンバ９０によりインキが供給されるアニロックスローラ９２と、被印刷体１６を加圧するための圧胴９４とから構成される。
　ドクターチャンバ９０から供給されたインキは、アニロックスローラ９２、凸版１４（網点突起部の頂面）の順に転移される。これに併せて、被印刷体１６は、版胴８８と圧胴９４との間に挟持された状態下に搬送される。すなわち、凸版１４に転移されたインキは被印刷体１６に更に転移されることにより、多数の網点からなる画像が被印刷体１６上に形成される。
［画像処理装置２０の電気的なブロック図］
　図６は、図１に示す画像処理装置２０の電気的なブロック図である。画像処理装置２０は、制御部１００と、通信Ｉ／Ｆ１０２と、入力部１０４と、表示部１０６と、出力Ｉ／Ｆ１０８と、メモリ１１０（記憶媒体）とを備えるコンピュータである。
　通信Ｉ／Ｆ１０２は、外部装置からの電気信号を送受信するインターフェース（Ｉ／Ｆ）である。これにより、画像処理装置２０は、印刷用データ３０を含む各種データを外部装置から取得可能である。
　入力部１０４は、マウス、トラックボール、キーボード、又はタッチパネル等の種々の入力デバイスで構成される。表示部１０６は、図示しない制御回路を介して、表示制御信号に応じた画像を表示する。入力部１０４による入力機能及び表示部１０６による表示機能を組み合わせることにより、ユーザ・インターフェースを実現する。
　出力Ｉ／Ｆ１０８は、電気信号としての深さデータを凸版形成装置２２に向けて送信するインターフェース（Ｉ／Ｆ）である。
　メモリ１１０は、制御部１００が各構成要素を制御するのに必要なプログラム及びデータ等を記憶している。本図例では、印刷用データ３０、複数種類の閾値データ４４及び後述する応答特性データ１１２がそれぞれ格納されている。
　メモリ１１０は、非一過性であり、且つ、コンピュータ読み取り可能な記憶媒体で構成されてもよい。ここで、コンピュータ読み取り可能な記憶媒体は、光磁気ディスク、ＲＯＭ、ＣＤ－ＲＯＭ、フラッシュメモリ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置である。また、この記憶媒体は、短時間に且つ動的にプログラムを保持するものであっても、一定時間プログラムを保持するものであってもよい。
　制御部１００は、ＣＰＵ（Central Processing Unit）等のプロセッサによって構成されている。制御部１００は、メモリ１１０に格納されたプログラムを読み出し実行することにより、ＲＩＰ処理部３２、２値化処理部３４、閾値データ設定部３６、深さデータ作成部３８（図１）及び応答特性推定部１１４（図１４）の各機能を実現可能である。
　閾値データ設定部３６は、設定対象の候補である閾値データ４４（以下、閾値候補）を取得する閾値候補取得部１１６、網点の形状変化を表現する空間周波数応答特性（以下、単に「応答特性」ともいう）を取得する応答特性取得部１１８、及び、各閾値候補における網点形状の再現性を定量評価する閾値候補評価部１２０を備える。
　閾値候補評価部１２０は、２値網点作成部１２２、模擬画像作成部１２４、定量値算出部１２６及び閾値データ決定部１２８を備える。なお、各部の機能については後述する。
［画像処理装置２０の動作］
　この実施形態に係る網点データ設定装置としての画像処理装置２０は、以上のように構成される。続いて、図１及び図６に示す画像処理装置２０（主に、閾値データ設定部３６）の動作について、図７のフローチャートを主に参照しながら詳細に説明する。
　ステップＳ１において、閾値データ設定部３６は、連続調画像信号４０に紐付けられた印刷条件４８を取得する。図１に示すように、印刷条件４８には、網点情報５０と、装置・印刷材情報５２が含まれる。
　網点情報５０として、例えば、スクリーン線数、スクリーン角度、網点の目標形状（円形、楕円形、矩形、これらの組合せ）、又はスクリーン特性（ＡＭスクリーン又はＦＭスクリーン）が挙げられる。装置・印刷材情報５２として、例えば、画像形成システム１０の装置構成、被印刷体１６、版材１２、色材、又は装置の構成部品（例えば、図４の露光ヘッド６２）の種類が挙げられる。
　ステップＳ２において、閾値候補取得部１１６は、ステップＳ１にて得た印刷条件４８と部分的に一致する閾値データ４４（以下、閾値候補ともいう）を取得する。この取得に先立ち、閾値ＤＢ２６は、画像処理装置２０からの要求に応じてデータベースの検索処理を実行し、網点情報５０に合致する閾値データ４４を画像処理装置２０側に供給する。その後、閾値候補取得部１１６は、メモリ１１０に一時的に記憶された閾値候補をすべて読み出し取得する。ここでは、３つの閾値候補Ａ、Ｂ、及びＣをそれぞれ取得したとする。
　ステップＳ３において、応答特性取得部１１８は、ステップＳ１にて得た印刷条件４８と部分的に一致する応答特性データ１１２を取得する。この取得に先立ち、応答特性ＤＢ２８は、画像処理装置２０からの要求に応じてデータベースの検索処理を実行し、装置・印刷材情報５２に合致する応答特性データ１１２を画像処理装置２０側に供給する。その後、応答特性取得部１１８は、メモリ１１０に一時的に記憶された応答特性データ１１２を読み出し取得する。
　なお、応答特性データ１１２は、２値画像信号４６（図２）の作成後から画像を形成するまでの画像形成過程にて生じる網点形状の変化を示すデータであり、画像形成システム１０の構成に依存する点に留意する。図１例では、応答特性取得部１１８は、応答特性データ１１２として、レーザ彫刻機６０（図４）の光学的特性に相当する「彫刻応答特性」、及びフレキソ印刷機２４（図５）の転移特性に相当する「転移応答特性」をそれぞれ取得する。
　図８Ａ及び図８Ｂは、彫刻応答特性の一例を示すグラフである。具体的には、図８Ａは点像分布関数｛以下、ＰＳＦ（Point Spread Function）｝であり、図８Ｂは図８Ａに示すＰＳＦのフーリエ変換関数の断面図に相当する。また、本グラフにおけるＸ軸、Ｙ軸は、図４に示す矢印Ｘ方向、矢印Ｙ方向にそれぞれ対応する。
　図８Ａに示すように、ＰＳＦは、原点Ｏをピークとする概略円錐形状の関数である。ＰＳＦは、理想系ではデルタ関数であるが、実際系では種々の光学的因子により点像が拡がる。Ｘ軸方向の拡がり程度は、レーザビーム８４（図４）の矢印Ｘ方向の有効ビーム径、ドラム６８（同図）の回転速度、レーザ光源（不図示）の立ち上がり及び／又は立ち下がり特性等に依存する。また、Ｙ軸方向の拡がり程度は、レーザビーム８４の矢印Ｙ方向の有効ビーム径、レーザビーム８４の露光多重度等に依存する。
　図８Ｂに示すＭＴＦ（Modulation Transfer Function）は、空間周波数毎の正弦波応答特性に相当する。本図から理解されるように、空間周波数が高くなるにつれてＭＴＦが次第に小さくなっている。そして、実線で示す「副走査」のＭＴＦは、破線で示す「主走査」よりも高くなる傾向を有する。
　ステップＳ４において、閾値候補評価部１２０は、各閾値候補の定量評価に供されるテスト画像信号を作成する。テスト画像信号として任意のラスタ画像信号４２を用いてもよいが、ここでは、少なくとも１つの平網領域（ハーフトーンべた領域）を含む画像信号を作成しておく。
　ステップＳ５において、閾値候補評価部１２０は、ステップＳ２により取得された複数の閾値候補（３つの閾値候補Ａ～Ｃ）のうち、未だ評価されていない１つの閾値候補を指定する。以下、指定された閾値候補のことを「指定閾値候補」と称する。
　ステップＳ６において、２値網点作成部１２２は、ステップＳ４で生成された各テスト画像信号に対して、指定閾値候補を用いたスクリーニング処理を施すことにより、階調レベル毎の２値網点信号を作成する。
　図９は、テスト画像信号に対して２値化処理を施した結果を示す説明図である。本図例では、階調レベルが１０％、３０％、５０％、７０％、及び９０％である５種類の平網領域を含むテスト画像信号から、濃淡が異なる５種類の平網画像を表す２値網点信号を作成する。
　ステップＳ７において、模擬画像作成部１２４は、ステップＳ６で作成された２値網点信号から、印刷物１８上での網点形状が模擬的に再現された模擬画像信号を作成する。
　図１０は、図６に示す模擬画像作成部１２４の動作に関する概略説明図である。この図は、２値網点信号から模擬画像信号を作成するまでの処理過程を模式的に示している。
　先ず、模擬画像作成部１２４は、各画素のオン又はオフ状態を示す２値網点信号を、凸版１４の高さ分布に相当する多値網点信号に変換する。その後、模擬画像作成部１２４は、得られた多値網点信号に対して空間周波数領域での応答特性を作用する処理を実行する。この作用処理には、マスク処理、フィルタ処理、又はこれらの組合せを適用できる。なぜならば、フーリエ変換の写像の性質上、ユークリッド空間上又はフーリエ空間上での演算結果は実質的に等価だからである。
　ここで、「マスク処理」は、ユークリッド空間上においてマスク演算（畳み込み）を実行することにより、画像の空間周波数成分（スペクトル強度）を変調する画像処理を意味する。また、「フィルタ処理」は、フーリエ空間上においてフィルタ演算（乗算）を実行することにより、画像信号の空間周波数成分（スペクトル強度）を変調する画像処理を意味する。
　本図上段に示す処理流れによれば、模擬画像作成部１２４は、彫刻応答特性を示すマスクを用いた第１マスク処理、転移応答特性を示すマスクを用いた第２マスク処理を順次施すことにより模擬画像信号を作成する。また、本図下段の処理流れによれば、模擬画像作成部１２４は、彫刻応答特性を示すフィルタを用いた第１フィルタ処理、転移応答特性を示すフィルタを用いた第２フィルタ処理を順次施すことにより模擬画像信号を作成する。ここで、フィルタ処理の実行前にＦＦＴ（Fast Fourier Transform）を、フィルタ処理の実行後にＩＦＦＴ（Inverse Fast Fourier Transform）をそれぞれ施す点に留意する。
　ステップＳ８において、定量値算出部１２６は、ステップＳ７で作成された模擬画像信号に対して所定の評価処理を実行して得られる評価値を、画質の程度を示す定量値として算出する。この実施形態では、定量値を算出するために、模擬画像信号のパワースペクトルを利用する。例えば、視認性との相関度を高めるため、このパワースペクトルに対して人間の標準視覚応答特性を乗算する前処理を行ってもよい。
　図１１は、観察距離が３００ｍｍであるドゥーリー・ショー（Dooley-Shaw）関数のグラフである。この関数は、ＶＴＦ（Visual Transfer Function）の一種であり、人間の標準視覚応答特性を模した代表的な関数である。具体的には、輝度のコントラスト比特性の２乗値に相当する。グラフの横軸は空間周波数（単位：Ｃｙｃｌｅ／ｍｍ）であり、縦軸はＶＴＦの値（単位は無次元）である。
　なお、ＶＴＦの関数形状はこれに限られず、数理モデルや実験データ等から導出された種々の特性を適用してもよい。また、観察距離は３００ｍｍ相当のみならず、画像の観察態様又は評価基準等に応じて種々変更してもよい。
　また、上記の前処理による視覚的な補正効果を適切に反映させるため、模擬画像信号の画素値を、画像を反射又は透過する光量に対して相関性が高い（好ましくは線形的な）量に変換しておいてもよい。その一例として、ＲＧＢ値、三刺激値（ＸＹＺ）、反射原稿の場合は光反射率、又は透過原稿の場合は光透過率等を用いてもよい。
　図１２Ａは、模擬画像信号におけるパワースペクトルの一例を示すグラフである。グラフの横軸は空間周波数（単位：Ｃｙｃｌｅ／ｍｍ）であり、グラフの縦軸はスペクトル強度（単位は任意；以下、単に「成分」ともいう）である。本図に示すパワースペクトルは、空間周波数が約１．９、約２．４、約３．６、約７．４Ｃｙｃｌｅ／ｍｍの位置に４つのピークを有する。定量値算出部１２６は、各網点に対応する成分（ここでは、約７．４Ｃｙｃｌｅ／ｍｍ周辺の成分）を除いた上で、スペクトル強度の最大値（主成分に相当；以下、最大強度という）を抽出する。この場合、定量値算出部１２６は、約２．４Ｃｙｃｌｅ／ｍｍにおけるピーク最大値である「６０」を、最大強度として抽出しておく。この最大強度の値が大きいほど画質の程度が低く、値が小さいほど画質の程度が高いと評価する。この定量値は、印刷物１８の画質に対して定量的又は定性的に相関する値である。なぜならば、各網点に対応する成分以外の成分は、色ムラ、モアレ、又は粒状感等を含む画像ノイズとして寄与するからである。
　なお、本図例のグラフは、閾値候補Ｂを用いて、階調レベルが７０％であるテスト画像信号から作成された模擬画像信号のパワースペクトルに相当する。定量値算出部１２６は、７０％の場合も含め、階調レベル毎の模擬画像信号に対して最大強度をそれぞれ抽出する。
　その後、定量値算出部１２６は、上記のように得られた階調レベル毎の最大強度に基づいて、指定閾値候補の定量値を算出する。ここでは、定量値算出部１２６は、全階調レベルにわたる平均値及び／又は最大値を「定量値」として算出する。なお、定量値の算出式はこれに限られず、各平網領域における最大強度についての統計値、又はこれらの組み合わせ（例えば、平均値と最大値の重み付け加算値）に基づいて種々算出してもよい。また、模擬画像信号に対する評価処理は上記した手法に限られず、粒状性｛ＲＭＳ（Root Mean Square）等｝、鮮鋭性（ＭＴＦ等）、又は階調性（階調滑らかさ等）を含む公知の画像評価手法を種々適用できることは言うまでもない。
　ステップＳ９において、閾値候補評価部１２０は、すべての閾値候補に対する指定が完了したか否かを判定する。まだ完了していないと判定された場合（ステップＳ９：ＮＯ）、ステップＳ５に戻って、以下ステップＳ５～Ｓ８を順次繰り返す。一方、すべて完了したと判定された場合（ステップＳ９：ＹＥＳ）、次のステップ（Ｓ１０）に進む。
　ステップＳ１０において、閾値データ決定部１２８は、閾値候補の比較及び/又は評価を行うことにより、１つの閾値データ４４を決定する。
　図１２Ｂは、各閾値候補Ａ～Ｃにおける階調レベル毎の最大強度を示すグラフである。図１３は、各閾値候補Ａ～Ｃに対する定量値の算出結果を示すテーブルである。閾値候補評価部１２０は、図１３に示す評価結果を踏まえ、「平均値」或いは「最大値」が最も小さい閾値候補Ａを、閾値データ４４として決定する。
　ステップＳ１１において、閾値データ設定部３６は、ステップＳ１０で決定された閾値データ４４を印刷用データ３０に対応付ける設定を行う。これにより、印刷用データ３０から深さデータを作成する際に、画像形成システム１０に適した閾値データ４４を自動的に設定することができる。このようにして、画像処理装置２０（主に閾値データ設定部３６）の動作は完了する。
［応答特性推定システム１５０の構成及び動作］
　続いて、応答特性推定部１１４（図６）を組み込んだ応答特性推定システム１５０の構成及び動作について、図１４～図１７Ｂを参照しながら説明する。
　図１４は、応答特性推定システム１５０の全体構成図である。この応答特性推定システム１５０は、標本網点信号１５２に基づき標本凸版１５４及び標本印刷物１５６を順次形成し、第１測定器１５８、第２測定器１６０による測定結果を用いて、画像形成システム１０の全体又は一部における応答特性を推定するシステムである。
　標本網点信号１５２は、任意の画像を構成する各網点の形状及び配置を表す２値の画像信号である。ここでは、後述する推定処理を容易にするため、標本網点信号１５２は、２値網点信号（図９）と同一の信号であってもよい。
　深さデータ作成部３８は、入力された標本網点信号１５２に基づいて、版材１２の彫刻量の２次元分布に相当する深さデータを作成する。そして、凸版形成装置２２は、深さデータ作成部３８から供給された深さデータに基づき、版材１２の主面を彫刻することにより、複数の網点突起部（図示しない）を備える標本凸版１５４を形成する。そして、フレキソ印刷機２４は、凸版形成装置２２により形成された標本凸版１５４を介し、被印刷体１６の版面にインキを転移することにより、標本印刷物１５６を形成する。
　第１測定器１５８は、標本凸版１５４の凹凸形状を光学的に読み取る装置であり、標本凸版１５４上での網点面積率を測定可能である。また、第２測定器１６０は、標本印刷物１５６の色特性を光学的に読み取る装置であり、標本印刷物１５６上での網点面積率を測定可能である。
　応答特性推定部１１４は、画像形成過程にて網点形状に変化を及ぼす事象に応じた応答特性を推定する。網点形状に変化を及ぼす事象として、画像形成装置の記録方式に応じた１つ又は２つ以上の事象を考慮してもよい。ここでは、応答特性推定部１１４は、版材１２の主面を彫刻することにより凸版１４を形成する際の「彫刻応答特性」、凸版１４を介して被印刷体１６にインキを転移させる際の「転移応答特性」をそれぞれ推定する。
　応答特性推定部１１４は、第１測定器１５８により測定された標本凸版１５４上での網点面積率に基づき彫刻応答特性を推定する彫刻特性推定部１６２と、第２測定器１６０により測定された標本印刷物１５６上での網点面積率に基づき転移応答特性を推定する転移特性推定部１６４とを備える。
　以下、彫刻特性推定部１６２による推定方法について、図１５Ａ～図１６Ｂを参照しながら詳細に説明する。
　図１５Ａは、図１４に示す標本網点信号１５２を可視化した概略正面図である。ここで、図１５Ａ及び後述する図１５Ｂ、図１５Ｃにおける塗り潰し領域は、各網点形状の決定に寄与する部位を示す。本図例では、各網点が円形状であり、直径がＤ（既知の値）である場合を想定する。
　図１５Ｂは、図１５Ａの標本網点信号１５２に基づき形成された標本凸版１５４の概略正面図である。図８Ａ及び図８Ｂに示す応答特性の作用下、標本凸版１５４上での網点面積率は、網点面積率の設計値よりも小さくなる。その結果、測定された網点面積率に相当する平均直径が（Ｄ－ｄ１）であったとする。
　図１５Ｃは、図１５Ｂの標本凸版１５４に基づき形成された標本印刷物１５６の概略正面図である。版材１２として高弾性材料が用いられる場合、インキを転移する際にドットゲインが増加する傾向がある。このとき、標本印刷物１５６上での網点面積率は、標本凸版１５４上での網点面積率よりも大きくなる。その結果、測定された網点面積率に相当する平均直径が（Ｄ－ｄ２）（ただし、ｄ２＜ｄ１）であったとする。
　ここで、測定された網点面積率から網点突起部の高さ形状を推定し、該高さ形状の変化に基づいて彫刻応答特性を推定する手法を適用できる。以下、具体的な方法について、図１６Ａ～図１７Ｂを参照しながら説明する。
　図１６Ａは、網点突起部の高さ形状の推定モデルに関する第１の説明図である。グラフの横軸はＸ軸（又はＹ軸）方向の位置であり、グラフの縦軸はＺ軸（高さ）方向の位置である。本図の実線は、標本網点信号１５２における各網点突起部の理想的な高さ形状（以下、理想網点形状１７０）を示している。本図の破線は、標本凸版１５４における各網点突起部の仮想的な高さ形状（以下、仮想網点形状１７２）を示している。この推定モデルは、網点突起部の高さが常に一定である前提の下に構築されており、主にドットゲインの現象を模擬している。
　図１６Ｂは、図１６Ａに示す高さ形状のフーリエ変換特性を示すグラフである。グラフの横軸は空間周波数であり、グラフの縦軸はフーリエ変換値である。本図の実線は、理想網点形状１７０に対するフーリエ変換値（以下、理想周波数成分１７０ｆ）を示している。本図の破線は、仮想網点形状１７２に対するフーリエ変換値（以下、仮想周波数成分１７２ｆ）を示している。
　この場合、彫刻特性推定部１６２は、理想周波数成分１７０ｆに対する仮想周波数成分
１７２ｆの比、すなわち、空間周波数毎の（仮想周波数成分１７２ｆ／理想周波数成分１７０ｆ）の値を、彫刻応答特性として推定する。
　図１７Ａは、網点突起部の高さ形状の推定モデルに関する第２の説明図である。グラフの定義は図１６Ａと同様であるため、その説明を省略する。本図の実線は理想網点形状１７０を示すと共に、本図の破線は仮想網点形状１７４を示している。
　この推定モデルは、頂部の直径にかかわらず、網点突起部の体積（面積）が常に保存される前提の下に構築されている。本図から理解されるように、頂部の直径が小さくなるにつれて側面の傾斜が小さくなっている。
　図１７Ｂは、図１７Ａに示す高さ形状のフーリエ変換特性を示すグラフである。グラフの定義は図１６Ｂと同様であるため、その説明を省略する。本図の実線は理想周波数成分１７０ｆを示すと共に、本図の破線は仮想周波数成分１７４ｆを示している。
　この場合、彫刻特性推定部１６２は、理想周波数成分１７０ｆに対する仮想周波数成分１７４ｆの比、すなわち、空間周波数毎の（仮想周波数成分１７４ｆ／理想周波数成分１７０ｆ）の値を、彫刻応答特性として推定する。
　図１４に戻って、転移特性推定部１６４は、応答特性の乗算則（カスケードモデル）に基づいて転移応答特性を推定する。先ず、転移特性推定部１６４は、図１６Ａ～図１７Ｂと同一の又は異なる推定モデルを用いて、画像形成システム１０全体における応答特性（以下、全体応答特性）を推定する。そして、転移特性推定部１６４は、全体応答特性から彫刻応答特性を除算することにより、転移応答特性を推定する。
　最後に、応答特性推定部１１４は、それぞれ推定された彫刻応答特性及び転移応答特性を応答特性ＤＢ２８側に供給する。そして、応答特性ＤＢ２８は、取得した彫刻応答特性を、版材１２及び凸版形成装置２２の種類に対応付ける。応答特性ＤＢ２８は、取得した転移応答特性を、版材１２及びフレキソ印刷機２４の種類に対応付ける。このようにして、応答特性データ１１２は、装置・印刷材情報５２に紐付けられた状態下、応答特性ＤＢ２８内にて保存・管理される。
［本発明の効果］
　閾値データ設定装置としての画像処理装置２０は、画像を構成する各網点の形状及び配置を表す２値画像信号４６を、連続調画像信号４０に基づき作成するための閾値データ４４を設定する装置である。
　そして、連続調画像信号４０に紐付けられた各網点の形状及び配置に関する網点情報５０を取得し、網点情報５０に合致する複数の閾値データ４４を各閾値候補として取得する閾値候補取得部１１６と、２値画像信号４６の作成後から画像を形成するまでの画像形成過程にて生じる網点形状の変化を、空間周波数領域での応答特性データ１１２として取得する応答特性取得部１１８と、応答特性データ１１２を作用することにより、網点形状の変化が模擬的に再現された模擬画像信号を作成する模擬画像作成部１２４と、模擬画像信号に対して所定の評価処理を実行して得られる評価値を、画質の程度を示す定量値として算出する定量値算出部１２６と、各閾値候補を適用しそれぞれ算出された定量値に基づく比較及び/又は評価によって、各閾値候補の中から設定対象である閾値データ４４を決定する閾値データ決定部１２８とを備える。
　このように、２値画像信号４６の作成後から画像を形成するまでの画像形成過程にて生じる網点形状の変化を、空間周波数領域での応答特性データ１１２として取得し、当該応答特性データ１１２を作用することにより、網点形状の変化が模擬的に再現された模擬画像信号を作成するので、網点形状を模擬画像上に簡便に且つ精度よく再現可能である。そして、各閾値候補として取得された複数の閾値データ４４を、画質の程度を示す定量値に基づき比較及び/又は評価するので、印刷装置及び印刷材の種類数が膨大な場合であっても、これらの組み合わせに適した閾値データ４４を簡便且つ容易に決定できる。
［補足］
　なお、この発明は、上述した実施形態に限定されるものではなく、この発明の主旨を逸脱しない範囲で自由に変更できることは勿論である。
　この実施形態ではフレキソ印刷を例に説明したが、画像形成装置はこの印刷方式に限られず、オフセット印刷、グラビア印刷、デジタル印刷機（オンデマンド印刷機）等にも適用できることは言うまでもない。
１０…画像形成システム　　　　　　　　　　１４…凸版
１６…被印刷体　　　　　　　　　　　　　　１８…印刷物
２０…画像処理装置　　　　　　　　　　　　２２…凸版形成装置
２４…フレキソ印刷機　　　　　　　　　　　３０…印刷用データ
３４…２値化処理部　　　　　　　　　　　　３６…閾値データ設定部
４０…連続調画像信号　　　　　　　　　　　４４…閾値データ
５０…網点情報　　　　　　　　　　　　　　５２…装置・印刷材情報
６０…レーザ彫刻機　　　　　　　　　　　　１００…制御部
１１２…応答特性データ　　　　　　　　　　１１４…応答特性推定部
１１６…閾値候補取得部　　　　　　　　　　１１８…応答特性取得部
１２０…閾値候補評価部　　　　　　　　　　１２２…２値網点作成部
１２４…模擬画像作成部　　　　　　　　　　１２６…定量値算出部
１２８…閾値データ決定部　　　　　　　　　１５０…応答特性推定システム
１５２…標本網点信号　　　　　　　　　　　１５４…標本凸版
１５６…標本印刷物　　　　　　　　　　　　１５８…第１測定器
１６０…第２測定器　　　　　　　　　　　　１７０…理想網点形状
１７２、１７４…仮想網点形状

Claims (10)

  1. 　画像を構成する各網点の形状及び配置を表す２値画像信号を、連続調画像信号に基づき作成するための閾値データを設定する装置であって、
    　前記連続調画像信号に紐付けられた前記各網点の形状及び配置に関する網点情報を取得し、該網点情報に合致する複数の閾値データを各閾値候補として取得する閾値候補取得部と、
    　前記２値画像信号の作成後から前記画像を形成するまでの画像形成過程にて生じる網点形状の変化を、空間周波数領域での応答特性として取得する応答特性取得部と、
    　前記応答特性取得部により取得された前記応答特性を作用することにより、前記網点形状の変化が模擬的に再現された模擬画像信号を作成する模擬画像作成部と、
    　前記模擬画像作成部により作成された前記模擬画像信号に対して所定の評価処理を実行して得られる評価値を、画質の程度を示す定量値として算出する定量値算出部と、
    　前記閾値候補取得部により取得された前記各閾値候補を適用し、前記定量値算出部にてそれぞれ算出された前記定量値に基づく比較及び/又は評価によって、前記各閾値候補の中から設定対象である閾値データを決定する閾値データ決定部と
    　を備える閾値データ設定装置。
  2. 　請求項１記載の閾値データ設定装置において、
    　前記模擬画像作成部は、平網領域を少なくとも１つ含む前記模擬画像信号を作成し、
    　前記定量値算出部は、前記模擬画像信号の前記平網領域でのパワースペクトルに関し、前記各網点に対応する成分を除いた主成分に基づいて前記定量値を算出する
    　閾値データ設定装置。
  3. 　請求項２記載の閾値データ設定装置において、
    　前記模擬画像作成部は、階調レベルが異なる平網領域を２つ以上含む前記模擬画像信号を作成し、
    　前記定量値算出部は、各前記平網領域における前記主成分についての統計値に基づいて前記定量値を算出する
    　閾値データ設定装置。
  4. 　請求項１～３のいずれか１項に記載の閾値データ設定装置において、
    　前記画像形成過程にて前記網点形状に変化を及ぼす事象に応じた前記応答特性を推定する応答特性推定部を更に備え、
    　前記応答特性取得部は、前記応答特性推定部により推定された前記応答特性を取得する
    　閾値データ設定装置。
  5. 　請求項４記載の閾値データ設定装置において、
    　前記応答特性推定部は、刷版を介して被印刷材に色材を転移させる際の前記応答特性を推定する閾値データ設定装置。
  6. 　請求項５記載の閾値データ設定装置において、
    　前記応答特性推定部は、版材の主面を彫刻することにより前記刷版を形成する際の前記応答特性を推定する閾値データ設定装置。
  7. 　請求項１～６のいずれか１項に記載の閾値データ設定装置において、
    　前記閾値データ決定部により決定された前記閾値データを設定することにより、前記連続調画像信号を前記２値画像信号に変換する２値化処理部を更に備える閾値データ設定装置。
  8. 　請求項７記載の閾値データ設定装置と、
    　前記閾値データ設定装置により変換された前記２値画像信号に基づいて前記画像を形成する画像形成装置と
    　を備える画像形成システム。
  9. 　画像を構成する各網点の形状及び配置を表す２値画像信号を、連続調画像信号に基づき作成するための閾値データを設定する方法であって、
    　前記連続調画像信号に紐付けられた前記各網点の形状及び配置に関する網点情報を取得し、該網点情報に合致する複数の閾値データを各閾値候補として取得するステップと、
    　前記２値画像信号の作成後から前記画像を形成するまでの画像形成過程にて生じる網点形状の変化を、空間周波数領域での応答特性として取得するステップと、
    　取得された前記応答特性を作用することにより、前記網点形状の変化が模擬的に再現された模擬画像信号を作成するステップと、
    　作成された前記模擬画像信号に対して所定の評価処理を実行して得られる評価値を、画質の程度を示す定量値として算出するステップと、
    　取得された前記各閾値候補を適用し、それぞれ算出された前記定量値に基づく比較及び/又は評価によって、前記各閾値候補の中から設定対象である閾値データを決定するステップと
    　をコンピュータに実行させることを特徴とする閾値データ設定方法。
  10. 　画像を構成する各網点の形状及び配置を表す２値画像信号を、連続調画像信号に基づき作成するための閾値データを設定するためのプログラムであって、
    　前記連続調画像信号に紐付けられた前記各網点の形状及び配置に関する網点情報を取得し、該網点情報に合致する複数の閾値データを各閾値候補として取得するステップと、
    　前記２値画像信号の作成後から前記画像を形成するまでの画像形成過程にて生じる網点形状の変化を、空間周波数領域での応答特性として取得するステップと、
    　取得された前記応答特性を作用することにより、前記網点形状の変化が模擬的に再現された模擬画像信号を作成するステップと、
    　作成された前記模擬画像信号に対して所定の評価処理を実行して得られる評価値を、画質の程度を示す定量値として算出するステップと、
    　取得された前記各閾値候補を適用し、それぞれ算出された前記定量値に基づく比較及び/又は評価によって、前記各閾値候補の中から設定対象である閾値データを決定するステップと
    　をコンピュータに実行させることを特徴とする閾値データ設定プログラム。
PCT/JP2014/064694 2013-07-23 2014-06-03 閾値データ設定装置、方法及びプログラム、並びに画像形成システム WO2015011991A1 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
EP14829027.3A EP3026882B1 (en) 2013-07-23 2014-06-03 Threshold value data setting device, method and program, and image forming system
US14/990,679 US9407788B2 (en) 2013-07-23 2016-01-07 Threshold value data setting device, method and program, and image forming system

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2013152904A JP5977716B2 (ja) 2013-07-23 2013-07-23 閾値データ設定装置、方法及びプログラム、並びに画像形成システム
JP2013-152904 2013-07-23

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US14/990,679 Continuation US9407788B2 (en) 2013-07-23 2016-01-07 Threshold value data setting device, method and program, and image forming system

Publications (1)

Publication Number Publication Date
WO2015011991A1 true WO2015011991A1 (ja) 2015-01-29

Family

ID=52393045

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2014/064694 WO2015011991A1 (ja) 2013-07-23 2014-06-03 閾値データ設定装置、方法及びプログラム、並びに画像形成システム

Country Status (4)

Country Link
US (1) US9407788B2 (ja)
EP (1) EP3026882B1 (ja)
JP (1) JP5977716B2 (ja)
WO (1) WO2015011991A1 (ja)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6472336B2 (ja) 2014-06-18 2019-02-20 キヤノン株式会社 画像処理装置、画像処理方法およびプログラム
US20180108185A1 (en) * 2015-06-17 2018-04-19 International Boxes S.R.L. Method and computer system for virtually displaying images on carton boxes
EP3438755B1 (en) 2017-08-04 2022-06-22 Canon Kabushiki Kaisha Image forming apparatus
EP3438757B1 (en) 2017-08-04 2021-03-24 Canon Kabushiki Kaisha Image forming apparatus
EP3438756B1 (en) * 2017-08-04 2021-09-08 Canon Kabushiki Kaisha Image forming apparatus
CN111241663A (zh) * 2020-01-02 2020-06-05 北京航天测控技术有限公司 发动机数据评估方法、装置、计算机设备和存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002368995A (ja) 2001-06-12 2002-12-20 Fuji Photo Film Co Ltd 階調画像作成用閾値配列決定方法
JP2004064689A (ja) * 2002-07-31 2004-02-26 Canon Inc 画像評価装置及び画像評価方法
JP2011139364A (ja) * 2009-12-28 2011-07-14 Canon Inc 画像評価装置、画像評価方法、及びコンピュータプログラム
JP2011239193A (ja) * 2010-05-11 2011-11-24 Konica Minolta Holdings Inc 閾値マトリクス生成方法、閾値マトリクス生成装置、閾値マトリクス、量子化装置及び画像形成装置
JP2012054761A (ja) * 2010-09-01 2012-03-15 Canon Inc 画像処理装置および画像処理方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5557709A (en) * 1994-07-01 1996-09-17 Seiko Epson Corporation Method and apparatus for dither array generation to reduce artifacts in halftoned images
US5760920A (en) * 1994-09-19 1998-06-02 Lin; Qian System and method for generating calibrated dither matrices
US8520277B2 (en) 2010-08-03 2013-08-27 Canon Kabushiki Kaisha Image processing apparatus, image processing method, and storage medium
US8467099B2 (en) * 2011-05-16 2013-06-18 Xerox Corporation Directional halftone method with enhanced edge and detail rendition
JP5856593B2 (ja) * 2013-08-30 2016-02-10 富士フイルム株式会社 色分解装置、色分解方法および色分解プログラム
JP5856594B2 (ja) * 2013-08-30 2016-02-10 富士フイルム株式会社 色分解装置、色分解方法および色分解プログラム

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002368995A (ja) 2001-06-12 2002-12-20 Fuji Photo Film Co Ltd 階調画像作成用閾値配列決定方法
JP2004064689A (ja) * 2002-07-31 2004-02-26 Canon Inc 画像評価装置及び画像評価方法
JP2011139364A (ja) * 2009-12-28 2011-07-14 Canon Inc 画像評価装置、画像評価方法、及びコンピュータプログラム
JP2011239193A (ja) * 2010-05-11 2011-11-24 Konica Minolta Holdings Inc 閾値マトリクス生成方法、閾値マトリクス生成装置、閾値マトリクス、量子化装置及び画像形成装置
JP2012054761A (ja) * 2010-09-01 2012-03-15 Canon Inc 画像処理装置および画像処理方法

Also Published As

Publication number Publication date
US9407788B2 (en) 2016-08-02
US20160117576A1 (en) 2016-04-28
EP3026882A4 (en) 2016-06-08
EP3026882A1 (en) 2016-06-01
JP2015023551A (ja) 2015-02-02
EP3026882B1 (en) 2018-08-15
JP5977716B2 (ja) 2016-08-24

Similar Documents

Publication Publication Date Title
US9407788B2 (en) Threshold value data setting device, method and program, and image forming system
JP7471268B2 (ja) フレキソ印刷の色調反応をさらに滑らかにする方法
US6775029B1 (en) Method for efficient calibration of printing devices
US20030217662A1 (en) Dot gain calibration method and apparatus
JP6541417B2 (ja) 画像処理装置、画像形成装置、画像処理方法及びプログラム
US7522312B2 (en) Threshold matrix, a method of generating the same, and a method of reproducing color image
US20090034006A1 (en) Stochastic halftone images based on screening parameters
JP5845797B2 (ja) 画像処理装置、画像処理方法およびプログラム
JP2006042332A (ja) 階調依存型マルチ周波数ハーフトーンスクリーニングのためのシステム及び方法
WO2009020522A1 (en) Method for generating stochastic dither matrix
US20100046035A1 (en) Information processing apparatus and method
JP2007049443A5 (ja)
US20150036946A1 (en) Metrics to identify image smoothness
Hladnik et al. Assessment of paper surface topography and print mottling by texture analysis
JPH0785272A (ja) 周波数変調ハーフトーン画像および作成方法
JP2007021905A (ja) 印刷支援装置、方法及びプログラム
Namedanian et al. Physical and optical dot gain: characterization and relation to dot shape and paper properties
JP2005252893A (ja) 閾値マトリクス
JP6537323B2 (ja) 画質調整情報を作成する装置、方法、及びプログラム
JPH10191045A (ja) 印画可能な媒体用の目視管理ストリツプ
JP2016049768A (ja) 印刷システム及び管理装置並びに管理方法
JP6310821B2 (ja) 印刷システム及び管理装置並びに管理方法
Kuznetsov et al. Halftones
Dharavath et al. Analysis of print attributes of amplitude modulated (am) vs. frequency modulated (fm) screening of multicolor offset printing
US11446923B2 (en) Systems and methods for improved flexographic printing with reduced mottling and increased highlight stability

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 14829027

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 2014829027

Country of ref document: EP

NENP Non-entry into the national phase

Ref country code: DE