JP2010136022A - 印刷制御装置、印刷システムおよび印刷制御プログラム - Google Patents

印刷制御装置、印刷システムおよび印刷制御プログラム Download PDF

Info

Publication number
JP2010136022A
JP2010136022A JP2008309030A JP2008309030A JP2010136022A JP 2010136022 A JP2010136022 A JP 2010136022A JP 2008309030 A JP2008309030 A JP 2008309030A JP 2008309030 A JP2008309030 A JP 2008309030A JP 2010136022 A JP2010136022 A JP 2010136022A
Authority
JP
Japan
Prior art keywords
color material
amount set
color
material amount
printing
Prior art date
Legal status (The legal status 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 status listed.)
Granted
Application number
JP2008309030A
Other languages
English (en)
Other versions
JP5157856B2 (ja
Inventor
Takashi Ito
隆志 伊藤
Atsushi Hoshii
淳 星井
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Seiko Epson Corp
Original Assignee
Seiko Epson Corp
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 Seiko Epson Corp filed Critical Seiko Epson Corp
Priority to JP2008309030A priority Critical patent/JP5157856B2/ja
Priority to US12/567,378 priority patent/US20100134550A1/en
Publication of JP2010136022A publication Critical patent/JP2010136022A/ja
Application granted granted Critical
Publication of JP5157856B2 publication Critical patent/JP5157856B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B41PRINTING; LINING MACHINES; TYPEWRITERS; STAMPS
    • B41JTYPEWRITERS; SELECTIVE PRINTING MECHANISMS, i.e. MECHANISMS PRINTING OTHERWISE THAN FROM A FORME; CORRECTION OF TYPOGRAPHICAL ERRORS
    • B41J2/00Typewriters or selective printing mechanisms characterised by the printing or marking process for which they are designed
    • B41J2/005Typewriters or selective printing mechanisms characterised by the printing or marking process for which they are designed characterised by bringing liquid or particles selectively into contact with a printing material
    • B41J2/01Ink jet
    • B41J2/205Ink jet for printing a discrete number of tones
    • B41J2/2056Ink jet for printing a discrete number of tones by ink density change
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N1/00Scanning, transmission or reproduction of documents or the like, e.g. facsimile transmission; Details thereof
    • H04N1/46Colour picture communication systems
    • H04N1/54Conversion of colour picture signals to a plurality of signals some of which represent particular mixed colours, e.g. for textile printing
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N1/00Scanning, transmission or reproduction of documents or the like, e.g. facsimile transmission; Details thereof
    • H04N1/46Colour picture communication systems
    • H04N1/56Processing of colour picture signals
    • H04N1/60Colour correction or control
    • H04N1/603Colour correction or control controlled by characteristics of the picture signal generator or the picture reproducer
    • H04N1/6033Colour correction or control controlled by characteristics of the picture signal generator or the picture reproducer using test pattern analysis
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N1/00Scanning, transmission or reproduction of documents or the like, e.g. facsimile transmission; Details thereof
    • H04N1/46Colour picture communication systems
    • H04N1/56Processing of colour picture signals
    • H04N1/60Colour correction or control
    • H04N1/6083Colour correction or control controlled by factors external to the apparatus
    • H04N1/6088Colour correction or control controlled by factors external to the apparatus by viewing conditions, i.e. conditions at picture output

Landscapes

  • Engineering & Computer Science (AREA)
  • Multimedia (AREA)
  • Signal Processing (AREA)
  • Textile Engineering (AREA)
  • Facsimile Image Signal Circuits (AREA)
  • Color Image Communication Systems (AREA)

Abstract

【課題】精度の高い再現を効率よく実現する。
【解決手段】ルックアップテーブルには、ターゲット値を特定するインデックスとターゲット値に近似した色再現可能なインク量セットとの対応関係が規定されている。該ルックアップテーブルには、濃色インクへのインク量割当てを優先しつつ最適化されたインク量セットを求めた後、このインク量セットを初期値として淡色インクへのインク量割当てを優先しつつ最適化されたインク量セットを求め、このインク量セットインデックスとの対応関係が規定されている。
【選択図】図15

Description

本発明は、印刷制御装置、印刷システムおよび印刷制御プログラムに関し、少なくとも1つの色相に関して略同一の色相を有し濃度が互いに異なる濃色材と淡色材を含んで構成された複数の色材を印刷装置にて記録媒体に付着させて印刷を実行させるにあたり、前記色材の使用量の組合せである色材量セットを前記印刷装置に指定し、前記色材量セットに基づく印刷を実行させる印刷制御装置、印刷システムおよび印刷制御プログラムに関する。
分光的な再現性に注目した印刷方法が提案されている(特許文献1参照。)。この文献においては、ターゲット画像に分光的かつ測色的に一致するような印刷を行うために、プリンティングモデルを使用し、ターゲットの分光反射率(ターゲットスペクトル)にフィッティングするようにプリンタ色(CMYKOG)の組合せを最適化している。このようにすれば、当該プリンタ色(CMYKOG)に基づく印刷を行うことにより、分光的にターゲット画像が再現でき、結果として測色的にも再現性の高い印刷結果を得ることができる。
特表2005−508125号公報
プリンティングモデルによれば、実際に印刷を行うことなく印刷結果を予測することが可能であるが、プリンティングモデルの予測結果が実際の印刷結果とずれる場合も考えられる。例えば、プリンティングモデル自体の精度が悪い場合や、プリンティングモデルが高精度であってもプリンティングモデルの予測における最適化条件(初期条件や目的関数の設定の仕方)が悪い場合には、プリンティングモデルによって予測したとおりの再現結果が得られないという問題が生じていた。
本発明は、前記課題にかんがみてなされたもので、精度の高い色再現を効率よく実現可能な印刷制御装置、印刷システムおよび印刷制御プログラムを提供する。
前記課題を解決するために、本発明の印刷制御装置は、少なくとも1つの色相に関して略同一の色相を有し濃度が互いに異なる濃色材と淡色材を含んで構成された複数の色材を印刷装置にて記録媒体に付着させて印刷を実行させるにあたり、前記色材の使用量の組合せである色材量セットを前記印刷装置に指定し、前記色材量セットに基づく印刷を実行させる印刷制御装置であって、印刷手段を備える構成としてある。
前記印刷手段は、前記色材量セットとインデックスとの対応関係を規定したルックアップテーブルを参照することにより、指定されたインデックスに対応する前記色材量セットを前記印刷装置に指定して印刷させる。すなわち、前記ルックアップテーブルには、対象物の色を示す情報であるターゲット値を特定するインデックスが含まれている。このインデックスに対応付けられた色材量セットは、前記印刷装置にて前記記録媒体上に付着したときに前記ターゲット値に対する近似性が極大化する色材量セット(ターゲット色材量セット)である。
ここで、ターゲット色材量セットは、次のようにして算出される。まず、前記淡色材の使用量を抑制しつつ(逆に言えば、前記濃色材の使用量を優先的に増加することである。)前記近似性を極大化するように所定の予測モデルに基づいて色材量セット(第1色材量セット)を予測する。次に、この第1色材量セットを前記所定の予測モデルの初期値として、前記濃色材の使用量を抑制しつつ(逆に言えば、前記淡色材の使用量を優先的に増加することである。)前記近似性を極大化するように前記所定の予測モデルに基づいて色材量セット(第2色材量セット)を予測する。このようにして算出された第2色材量セットが、ターゲット色材量セットである。すなわち、ターゲット値に近似する色を前記記録媒体上に再現する色材量セットにおいて各色材量を決定するにあたり、濃色材を優先的に割り当てることにより前記記録媒体に対するトータルでの色材の付着量の増大を抑制した色材量セットを算出することができ、濃色材で表現できない細かい色の濃さについて淡色材を優先的に割り当てることにより精度の高い色再現が可能な色材量セットを算出できる。
なお、前記ターゲット値は、前記対象物の分光反射率や色彩値とすることができる。前記ターゲット値を分光反射率とすると、分光反射率の再現性の良好な印刷を前記印刷装置に実行させることができる。この場合、前記予測モデルは、任意の前記色材量セットで印刷を行った場合の分光反射率を予測する。さらに、前記ターゲット値は、前記ターゲットの複数光源下における色彩値とすることにより、複数光源下での色の再現性が良好な印刷を前記印刷装置に実行させることができる。この場合、前記予測モデルは、任意の前記色材量セットで印刷を行った場合の複数光源下における色彩値を予測する。また、前記印刷装置は少なくとも複数の前記色材を前記記録媒体に付着させることができればよく、インクジェットプリンタやレーザープリンタや昇華型プリンタ等の種々の印刷装置に本発明を適用することができる。
また、本発明の選択的な一態様として、前記ターゲット値は、次のようにして求められた補正ターゲット値であってもよい。補正ターゲット値は、前記印刷装置にて前記記録媒体上に前記ターゲット値を再現するための前記色材量セットを前記所定の予測モデルに基づいて予測し、この予測した色材量セットを前記印刷装置に指定して確認パッチを印刷させ、この確認パッチの色を示す情報である確認ターゲット値と前記対象物の測色値であるターゲット測定値との偏差に基づいて設定される値である。
このようにすることにより、前記偏差を解消するような予測結果を得ることができる。従って、プリンティングモデル自体の精度が悪い場合や、プリンティングモデルが高精度であってもプリンタの再現特性に経時変化がある場合や、個々のプリンタの再現特性にばらつきがある場合であっても、精度の高い予測結果を得ることができる。また、前記ターゲット値から前記偏差を単純に引くのではなく、例えば前記偏差の何割かを引くようにしてもよい。
また、本発明の選択的な一態様として、前記第2色材量セットを前記印刷装置に指定して再確認パッチを印刷させ、この再確認パッチの色を示す情報である再確認ターゲット値と前記ターゲット測定値との偏差に基づいて算出した再補正ターゲット値をターゲット値として前記第1色材量セットと前記第2色材量セットの再予測を行う構成としてもよい。
すなわち、予測された第2色材量セットに基づいて、前記印刷装置に実際に印刷を行わせて、その印刷結果を測色してこれを新たなターゲット値(再補正ターゲット値)とするのである。そしてこの新たなターゲット値を再現するための色材量セットを予測する。従って、予測された第2色材量セットの印刷結果に基づくフィードバックがかかり、より予測精度を向上することが期待できる。
また、本発明の選択的な一態様として、前記所定の予測モデルにおいては、色材量セットを予測する際に、前記印刷装置にて付着することが可能な最小の単位量よりも少量ずつ色材量を変化させながら色材量セットの前記近似性を評価しており、前記所定の予測モデルに基づいて予測される色材量セットは、前記近似性が極大化した色材量セットについて前記単位量を丸め幅とした端数処理を行ったものであってもよい。端数処理を行うと、丸め誤差が発生する。濃色材の色材量に起因して発生する丸め誤差は、濃色材の単位量よりも小さいので、再予測を行って濃色材の色材量を変更したとしても、同様の丸め誤差が発生する可能性が高い。そこで、前述したように、濃色材を優先的に割り当てた色材量セットについて端数処理を行ったものに対し、今度は、淡色材を優先的に割り当てる予測を行う。その結果、濃色材の丸め誤差分の色材量は、淡色材の色材量で補われて淡色材の色材量に転化される可能性が高くなる。よって、色再現性精度の高い色材量セットが予測できる。
また、本発明の選択的な一態様として、予測された前記第2色材量セットを初期値として前記第1色材量セットと前記第2色材量セットを予測する過程を複数回繰り返し、該繰り返しの過程において、前記第2色材量セットにおける前記濃色材の量が2回連続して一致したことを検出した場合は、それ以降の繰り返しにおいて前記濃色材の量を固定する構成としてもよい。
一般的に、最適化処理を複数回繰り返すことにより、最適化の精度が向上する。しかしながら、本発明の濃色材のように変化のステップが大きい場合は、繰り返し数を増やしてもあまり精度向上は見込めないことを考慮すると、濃色材については2回連続して同一の値に最適化された場合は、以降の最適化を行わないことにより処理時間の短縮を図る方のメリットが高い。これは、第2色材量セットが最適解近傍にある場合などに顕著であり、この場合に必要な処理は前記丸め誤差を補うように最適化する処理であるため、濃色材を変化させて予測を行うと、最適解からの誤差が増えるばかりで、メリットが全く無い。
また、本発明の選択的な一態様として、前記第1色材量セットと前記第2色材量セットの各々の予測において、使用量を抑制されたインクを除いた色材組合せによって、全色相方向への色変更が可能になっていてもてもよい。このように全色相方向への自由度を確保するための具体的な一態様として、前記複数の色材が、シアン(C)とマゼンタ(M)とイエロー(Y)とブラック(K)とライトシアン(lc)とライトマゼンタ(lm)とを含んで構成されており、前記第1色材量セットの予測においては、少なくともシアン(C)とマゼンタ(M)とブラック(K)の使用量を優先的に変化し、前記第2色材量セットの予測においては、少なくともライトシアン(lc)とライトマゼンタ(lm)とイエロー(Y)の使用量を優先的に変化する構成としてもよい。
前記構成において、CMYKは濃色材であり、lclmは淡色材である。ただし、このインクセット組合せでは、lcとlmだけではイエロー方向への色変更が出来ないため、この2色だけで予測モデルによる予測を行ったとしても、イエロー方向への誤差が解消されないため予測精度が低下する。そこで、淡色材を優先的に変化する場合には、Yをlclmと一緒に変化させながら予測を行うことにより、全色相方向への色変化を可能にするのである。無論、イエロー方向への色変化を発生することのできる淡色材(ライトイエロー等)を色材セットに備えていれば、Yを淡色材と一緒に変化させずに、CMYKで第1色材量セットの予測を行ってもよい。
さらに、本発明の技術的思想は、具体的な印刷制御装置にて具現化されるのみならず、その方法としても具現化することができる。すなわち、上述した印刷制御装置が行う各手段に対応する工程を有する方法としても本発明を特定することができる。むろん、上述した印刷制御装置がプログラムを読み込んで上述した各手段を実現する場合には、当該各手段に対応する機能を実行させるプログラムや当該プログラムを記録した各種記録媒体においても本発明の技術的思想が具現化できることは言うまでもない。なお、本発明の印刷制御装置は、単一の装置のみならず、複数の装置によって分散して存在可能であることはいうまでもない。例えば、印刷制御装置の状態を示す各手段が、パーソナルコンピュータ上で実行されるプリンタドライバと、プリンタの双方において分散することも可能である。また、プリンタ等の印刷装置に本発明の印刷制御装置の各手段を包含させることも可能である。
以下、下記の順序に従って本発明の実施形態を説明する。
1.印刷制御装置の構成:
2.印刷データ生成処理:
3.印刷制御処理:
3−1.1D−LUT作成処理:
3−2.印刷制御データ生成処理:
4.キャリブレーション処理:
5.分光プリンティングモデル:
6.変形例:
6−1:変形例1:
6−2:変形例2:
6−3:変形例3:
6−4:変形例4:
6−5:変形例5:
6−6:変形例6:
6−7:変形例7:
6−8:変形例8:
1.印刷制御装置の構成
図1は、本発明の一実施形態にかかる印刷制御装置のハードウェア構成を示している。同図において、印刷制御装置は主にコンピュータ10によって構成されており、コンピュータ10はCPU11とRAM12とROM13とハードディスクドライブ(HDD)14と汎用インターフェイス(GIF)15とビデオインターフェイス(VIF)16と入力インターフェイス(IIF)17とバス18とから構成されている。バス18は、コンピュータ10を構成する各要素11〜17の間でのデータ通信を実現するものであり、図示しないチップセット等によって通信が制御されている。HDD14には、オペレーティングシステム(OS)を含む各種プログラムを実行するためのプログラムデータ14aが記憶されており、当該プログラムデータ14aをRAM12に展開しながらCPU11が当該プログラムデータ14aに準じた演算を実行する。GIF 15は、例えばUSB規格に準じたインターフェイスを提供するものであり、外部のプリンタ20と分光反射率計30をコンピュータ10に接続させている。VIF 16はコンピュータ10を外部のディスプレイ40に接続し、ディスプレイ40に画像を表示するためのインターフェイスを提供する。IIF 17はコンピュータ10を外部のキーボード50aとマウス50bに接続し、キーボード50aとマウス50bからの入力信号をコンピュータ10が取得するためのインターフェイスを提供する。
図2は、コンピュータ10にて実行されるプログラムのソフトウェア構成を概略的なデータの流れとともに示している。同図において、コンピュータ10では、おもにOS P1と見本印刷アプリケーション(APL)P2と1D−LUT生成アプリケーション(LUG)P3aとプリンタドライバ(PDV)P3bと分光反射率計ドライバ(MDV)P4とディスプレイドライバ(DDV)P5が実行されている。OS P1は、各プログラムが使用可能なAPIのひとつとして画像機器インターフェイス(GDI) P1aとスプーラP1bを提供しており、APL P2の要求に応じてGDI P1aが呼び出され、さらにGDI P1aの要求に応じてPDV P3bやDDV P5が呼び出される。GDI P1aはコンピュータ10がプリンタ20やディスプレイ40等の画像出力装置における画像出力を制御するための汎用的な仕組みを提供し、一方のPDV P3bやDDV P5はプリンタ20やディスプレイ40の機種固有の処理等を提供する。また、スプーラP1bは、APL P2やPDV P3bやプリンタ20の間に介在し、ジョブのコントロール等を実行する。APL P2は、見本チャートSCを印刷するためのアプリケーションプログラムであり、RGBビットマップ形式の印刷データPDを生成し、GDI P1aに対して当該印刷データPDを出力する。また、印刷データPDを生成するにあたっては、MDV P4からターゲットの分光反射率データRDを取得する。MDV P4は、APL P2の要求に応じて分光反射率計30を制御し、当該制御によって得られた分光反射率データRDをAPL P2に出力する。
APL P2が生成した印刷データPDはGDI P1aやスプーラP1bを経由してPDV P3bに出力され、PDV P3bが印刷データPDに基づいてプリンタ20に出力可能な印刷制御データCDを生成する処理を実行する。PDV P3bが生成した印刷制御データCDはOS P1が提供するスプーラP1bを介してプリンタ20に出力され、プリンタ20が当該印刷制御データCDに基づく動作を行うことにより見本チャートSCを印刷用紙上に印刷させる。以上においては、全体の処理の流れ概略的に説明したが、以下、フローチャートを用いて各プログラムP1〜P4が実行する処理を詳細に説明する。
2.印刷データ生成処理
図3は、APL P2が実行する印刷データ生成処理の流れを示している。図2に示すようにAPL P2はUI部(UIM)P2aと測定制御部(MCM)P2bと印刷データ生成部(PDG)P2cとから構成されており、これらの各モジュールP2a,P2b,P2cが図3に示す各ステップを実行する。ステップS100においては、UIM P2aがGDI P1aおよびDDV P5を介して、見本チャートSCを印刷させる旨の印刷指示を受け付けるためのUI画面を表示させる。前記UI画面においては、見本チャートSCのテンプレートを示す表示が設けられている。
図4は、前記UI画面の一例を示している。同図において、前記テンプレートTPが表示されており、当該テンプレートTPにはカラーパッチをレイアウトするための12個の枠FL1〜FL12が設けられている。前記UI画面には各枠FL1〜FL12をマウス50bのクリックによって選択することが可能となっており、枠FL1〜FL12をクリックすると分光反射率測定を開始させるか否かを指示するための選択ウィンドウWが表示される。また、前記UI画面においては、見本チャートSCの印刷を実行させるか否かを指示するためのボタンBも設けられている。ステップS110においては、UIM P2aがマウス50bによる各枠FL1〜FL12のクリックが検出し、検出された場合にはステップS120にて分光反射率測定を開始させるか否かを指示するための選択ウィンドウWを表示させる。ステップS130においては選択ウィンドウWにおけるマウス50bのクリックを検出し、キャンセルがクリックされた場合にはステップS110に戻る。一方、分光反射率測定実行がクリックされた場合には、ステップS140においてMCM P2bがMDV P4を介して分光反射率計30にターゲットTGの分光反射率R(λ)としてターゲット分光反射率Rt(λ)の測定を実行させ、当該ターゲット分光反射率Rt(λ)を格納した分光反射率データRDを取得する。ターゲット分光反射率Rt(λ)は、本発明のターゲットの状態を示す状態値およびターゲット値に相当する。
ステップS140におけるターゲット分光反射率Rt(λ)の測定が完了すると、最も標準的な光源であるD65光源を照射したときのCIELAB色空間における色彩値(L***値)を算出する。そして、当該L***値を所定のRGBプロファイルを使用してRGB値に変換し、当該RGB値を表示用RGB値として取得する。なお、RGBプロファイルは絶対色空間としてのCIELAB色空間と本実施形態のRGB色空間との等色関係を規定したプロファイルであり、例えばICCプロファイルを使用することができる。
図5は、ステップS140において、分光反射率データRDから表示用RGB値を算出する様子を模式的に示している。ターゲットTGについてターゲット分光反射率Rt(λ)を測定した結果、図示するようターゲット分光反射率Rt(λ)の分布を示す分光反射率データRDが得られたとする。なお、ターゲットTGは、分光的な再現の目標とする物体表面を意味し、例えば他の印刷装置や塗装装置等で形成した人工的物体表面や、自然物の表面等が該当する。一方、D65光源は図示するような可視波長域において不均一な分光エネルギーP(λ)の分布を有しており、ターゲットTGにD65光源を照射したときの各波長の反射光の分光エネルギーは、ターゲット分光反射率Rt(λ)と分光エネルギーP(λ)を各波長について掛け合わせた値となる。さらに、反射光の分光エネルギーのスペクトルに対して人間の分光感度特性に応じた等色関数x(λ),y(λ),z(λ)をそれぞれ畳み込み積分し、係数kによって正規化することにより、3刺激値X,Y,Zを得る。以上を数式で表すと下記の(1)式となる。
Figure 2010136022
3刺激値X,Y,Zを所定の変換式によって変換することにより、ターゲットTGにD65光源を照射したときの色を示すL***値を得ることができ、さらにRGBプロファイルを使用することにより、表示用RGB値を得ることができる。ステップS145においては、テンプレートTPにおいてクリックされた枠FL1〜FL12を前記表示用RGB値で塗りつぶした表示に更新する。これにより、標準的な光源であるD65光源でのターゲットTGの色をUI画面にて感覚的に把握することが可能となる。ステップS145が完了すると、ステップS150において固有のインデックスを生成するとともに、当該インデックスと、ステップS110にてクリックされた枠FL1〜FL12の位置情報と表示用RGB値を分光反射率データRDに対応付けてRAM12に記憶する。ステップS150が完了すると、ステップS110に戻り、ステップS120〜S150を繰り返して実行する。これにより、他の枠FL1〜FL12を選択し、他の枠FL1〜FL12について他のターゲットTGのターゲット分光反射率Rt(λ)の測定を行うことができる。
本実施形態においては、それぞれ異なる12種類のターゲットTG1〜TG12が用意されており、ターゲットTG1〜TG12のそれぞれについてのターゲット分光反射率Rt(λ)が分光反射率測定データRDとして取得されるものとする。従って、ステップS150においては、各枠FL1〜FL12についての分光反射率測定データRDと固有のインデックスを対応付けたデータが順次RAMに記憶されていくこととなる。なお、インデックスはそれぞれの値が固有となるように生成されればよく、インクリメントによって生成してもよいし、重複しない乱数によって生成してもよい。
ステップS110において、各枠FL1〜FL12のクリックが検出されない場合には、ステップS160にて見本チャートSCの印刷を実行させる旨のボタンBのクリックを検出し、検出されない場合にはステップS110に戻る。一方、見本チャートSCの印刷を実行させる旨のボタンBのクリックが検出された場合には、ステップS170にてPDG P2cが印刷データPDを生成する。
図6は、印刷データPDの構成を模式的に示している。同図において、印刷データPDはドットマトリクス状に配列した多数の画素によって構成されており、各画素が4バイト(8ビット×4)の情報を有している。印刷データPDは、図4に示したテンプレートTPと同様の画像を示しており、テンプレートTPの各枠FL1〜FL12に対応する領域以外の画素は、テンプレートTPに対応する色のRGB値を有している。RGB各チャネルの階調値はそれぞれ8ビット(256階調)によって表現され、上述した4バイトのうち3バイトがRGB値を格納するために使用される。例えば、テンプレートTPの各枠FL1〜FL12以外の色が(R,G,B)=(128,128,128)の一様な中間グレーで表される場合、印刷データPDにおける各枠FL1〜FL12に対応する領域以外の画素は(R,G,B)=(128,128,128)の色情報を有することとなる。なお、残りの1バイトは使用されない。
一方、テンプレートTPの各枠FL1〜FL12に対応する画素も4バイトの情報を有しており、通常、RGB値が格納される3バイトを使用してインデックスを格納する。このインデックスは、ステップS150にて各枠FL1〜FL12ごとに生成した固有のものであり、PDG P2cはインデックスをRAM12から取得し、各枠FL1〜FL12に対応する画素に対応するインデックスを格納する。このようにインデックスをRGB値の代わりに格納した各枠FL1〜FL12に対応する画素については、残りの1バイトを使用してインデックスが格納された旨のフラグを立てる。これにより、各画素がRGB値を格納しているか、インデックスを格納しているかを判別することができる。本実施形態では、インデックスを格納するために3バイトを使用することができるため、3バイト以下の情報量で表現できるインデックスをステップS150で生成しておく必要がある。以上のようにしてビットマップ形式の印刷データPDが生成できると、ステップS180において、PDG P2cがインデックステーブルIDBを生成する。
図7は、インデックステーブルIDBの一例を示している。同図において、各枠FL1〜FL12に対応して生成された固有のインデックスのそれぞれに対して、測定によって得られたターゲット分光反射率Rt(λ)およびD65光源におけるL***値に対応する表示用RGB値が格納されている。インデックステーブルIDBの生成が完了すると、印刷データPDはGDI P1aやスプーラP1bを経由してPDV P3bに出力される。印刷データPDは、外形上、通常のRGBビットマップ形式と変わらないため、OS P1が提供するGDI P1aやスプーラP1bにおいても通常の印刷ジョブと同様に処理することができる。一方、インデックステーブルIDBは、直接、PDV P3bに出力される。なお、本実施形態においては、インデックステーブルIDBを新たに生成するようにしたが、既存のインデックステーブルIDBにインデックスとターゲット分光反射率Rt(λ)と表示用RGB値の新たな対応関係を追記するようにしてもよい。また、以上の印刷データ生成処理と後述する印刷制御処理は必ずしも同一の装置において連続して実行する必要はなく、印刷データ生成処理と印刷制御処理を例えばLANやインターネット等の通信回線によって接続された複数のコンピュータ上にて個別に実行してもよい。
3.印刷制御処理
図8は、LUG P3aとPDV P3bが実行する印刷制御処理の全体的な流れを示している。図8に示す1D−LUT生成処理(ステップS200)をLUG P3aが担当し、一方の印刷制御データ生成処理(ステップS300)をPDV P3bが担当する。1D−LUT生成処理は印刷制御データ生成処理に先行して行われてもよいし、1D−LUT生成処理と印刷制御データ生成処理を並行して行うようにしてもよい。
3−1.1D−LUT生成処理
図9は、1D−LUT生成処理の流れを示している。図2に示すようにLUG P3aは、インク量セット算出モジュール(ICM)P3a1と分光反射率予測モジュール(RPM)P3a2と評価値算出モジュール(ECM)P3a3とLUT出力モジュール(LOM) P3a4とから構成されている。ステップS210においては、ICM P3a1がインデックステーブルIDBを取得する。ステップS220においては、インデックステーブルIDBから一つのインデックスを選択し、当該インデックスに対応付けられている分光反射率データRDを取得する。ステップS230においては、ICM P3a1が前記分光反射率データRDが示すターゲット分光反射率Rt(λ)と同様の分光反射率R(λ)が再現可能なインク量セットを算出する処理を行う。その際に、上述したRPM P3a2とECM P3a3を使用する。
図10は、前記分光反射率データRDが示すターゲット分光反射率Rt(λ)と同様の分光反射率R(λ)が再現可能なインク量セットを算出する処理の流れを模式的に示している。RPM P3a2は、ICM P3a1からのインク量セットφの入力に応じて、当該インク量セットφに基づいてプリンタ20が所定の印刷用紙にインクを吐出させたときの分光反射率R(λ)を予測し、当該分光反射率R(λ)を予測分光反射率Rs(λ)としてECM P3a3に出力する。
ECM P3a3は、分光反射率データRDが示すターゲット分光反射率Rt(λ)と予測分光反射率Rs(λ)の差分D(λ)を各波長λについて算出し、各波長λごとに重みが課せられた重み関数w(λ)を当該差分D(λ)に乗算する。この値の二乗平均の平方根を評価値E(φ)として算出する。以上の計算を数式で表すと下記の(2)式のように表すことができる。
Figure 2010136022
前記の(2)式において、Nは波長λの有限区分数を意味する。前記の(2)式において、評価値E(φ)が小さければ小さいほど、各波長λにおけるターゲット分光反射率Rt(λ)と予測分光反射率Rs(λ)の差が少ないということができる。すなわち、評価値E(φ)が小さければ小さいほど、入力したインク量セットφによってプリンタ20が印刷したときに記録媒体上にて再現される分光反射率R(λ)と、対応するターゲットTGから得られたターゲット分光反射率Rt(λ)が近似しているということができる。さらに、前記の(1)式によれば、光源の変動に応じてインク量セットφによってプリンタ20が印刷したときの記録媒体と対応するターゲットTGが示す絶対的な色彩値は双方とも変動するものの、分光反射率R(λ)が近似すれば光源の変動に拘わらず相対的には同じ色に知覚されるということができる。従って、評価値(φ)が小さくなるインク量セットφによれば、あらゆる光源においてターゲットTGと同じ色に知覚される印刷結果を得ることができるということができる。
また、本実施形態において、重み関数w(λ)は下記の(3)式のものを使用する。
Figure 2010136022

前記の(3)式においては、等色関数x(λ),y(λ),z(λ)を加算することにより、重み関数w(λ)が定義されている。なお、前記の(3)式の右辺全体に所定の係数を乗算して、重み関数w(λ)の値の範囲を正規化してもよい。前記の(1)式によれば、等色関数x(λ),y(λ),z(λ)が大きい波長域ほど、色彩値(L***値)に大きく影響するということができる。従って、等色関数x(λ),y(λ),z(λ)を加算した重み関数w(λ)を使用すれば、色への影響が大きい波長域を重視した二乗誤差が評価可能な評価値E(φ)を得ることができる。例えば、人間の目に知覚されない近紫外波長域においてはw(λ)が0となり、当該波長域における差分D(λ)は評価値E(φ)の増大に寄与しないこととなる。
すなわち、必ずしも全可視波長域においてターゲット分光反射率Rt(λ)と予測分光反射率Rs(λ)との差が小さくなくても、人間の目に特に強く知覚される波長域においてターゲット分光反射率Rt(λ)と予測分光反射率Rs(λ)とが似ていれば小さい値の評価値E(φ)を得ることができ、人間の目に知覚に即した分光反射率R(λ)の近似性の指標として評価値E(φ)を使用することができる。算出された評価値E(φ)はICM P3a1に返される。すなわち、ICMP 3a1が任意のインク量セットφをRPM P3a2とECM P3a3に出力することにより、最終的に評価値E(φ)がICM P3a1に返される構成となっている。ICM P3a1は任意のインク量セットφに対応して評価値E(φ)を得ることを繰り返し実行することにより、目的関数としての評価値E(φ)が極小化するようなインク量セットφの最適解を算出する。この最適解を算出する手法としては、例えば勾配法といった非線形最適化手法を用いることができる。
図11は、ステップS230においてインク量セットφが最適化されていく様子を模式的に示している。同図において、インク量セットφが最適化されていくにつれて、インク量セットφで印刷を行った場合の予測分光反射率Rs(λ)がターゲット分光反射率Rt(λ)に近づいていく。また、重み関数w(λ)を使用することにより、等色関数x(λ),y(λ),z(λ)が大きい波長域ほど、予測分光反射率Rs(λ)のターゲット分光反射率Rt(λ)への拘束が強くなっており、予測分光反射率Rs(λ)がターゲット分光反射率Rt(λ)の差が小さくなっている。このように、等色関数x(λ),y(λ),z(λ)が大きく、視覚に大きく影響する波長域について優先的に予測分光反射率Rs(λ)をターゲットTGのターゲット分光反射率Rt(λ)に拘束することができるため、任意の光源を照射したときの見た目が近くなるようなインク量セットφを算出することができる。以上により、いずれの光源においてもターゲットTGと似たような見た目をプリンタ20にて再現させることが可能なインク量セットφを算出することができる。なお、最適化の終了条件は、インク量セットφ更新の繰り返し回数としてもよいし、評価値E(φ)の閾値としてもよい。
以上のようにして、ステップS230においてICM P3a1がターゲットTGと同様の分光反射率R(λ)が再現可能なインク量セットφを算出すると、ステップS240においてインデックステーブルIDBに記述されたインデックスのすべてがステップS220にて選択されたか否かを判定し、すべて選択されていない場合にはステップS220に戻り、次のインデックスを選択する。このようにすることにより、すべてのインデックスについてターゲットTGと同様の色が再現可能なインク量セットφを算出することができる。すなわち、印刷データ生成処理(図2)のステップS140において測色を行ったすべてのターゲットTG1〜TG12についてターゲットTG1〜TG12と同様の分光反射率R(λ)が再現可能なインク量セットφを算出することができる。ステップS240において、すべてのインデックスについて最適なインク量セットφが算出されたことが判定されると、ステップS250において、LOM P3a4が1D−LUTを生成し、当該1D−LUTをPDV P3bに出力する。
図12は、1D−LUTの一例を示している。同図において、各インデックスに対応して最適なインク量セットφが格納されている。すなわち、各ターゲットTG1〜TG12について、各ターゲットTG1〜TG12と似たような見た目をプリンタ20にて再現させることが可能なインク量セットφを記述した1D−LUTを用意することができる。1D−LUTをPDV P3bに出力すると、1D−LUT生成処理が完了し、次の印刷制御データ生成処理(ステップS300)を実行させる。
3−2.印刷制御データ生成処理
図13は、印刷制御データ生成処理の流れを示している。図2に示すようにPDV P3bは、モード判別モジュール(MIM)P3b1とインデックス分版モジュール(ISM)P3b2とRGB分版モジュール(CSM)P3b3とハーフトーンモジュール(HTM)P3b4とラスタ化モジュール(RTM)P3b5とから構成されている。ステップS310においては、モード判別モジュール(MIM)P3b1が印刷データPDを取得する。ステップS320において、MIM P3b1は印刷データPDから一つの画素を選択する。ステップS330において、MIM P3b1は当該選択した画素にインデックスが格納された旨のフラグが立っているか否かを判定する。当該フラグが立っていないと判定された場合には、ステップS340にてCSM P3b3が3D−LUTを参照して、当該画素についての色変換(分版)を実行する。
図14は、3D−LUTの一例を示している。同図において、3D−LUTはRGB値とインク量セットφ(d,d,d,d,dlc,dlm)との対応関係が色空間における複数の代表的な座標について記述されたテーブルであり、CSM P3b3は3D−LUTを参照して当該画素が有するRGB値に対応するインク量セットφを取得する。その際に、3D−LUTに直接記述されていないRGB値については補間演算を行うことにより、対応するインク量セットφを取得する。なお、3D−LUTの作成方法として、特開2006−82460号公報等を採用することができる。当該公報においては、特定光源における色の再現性や、再現色の階調性や、粒状性や、再現色の光源非依存性や、ガマットや、インクデューティが総合的に良好となる3D−LUTが作成される。
一方、ステップS330において、選択した画素にインデックスが格納された旨のフラグが立っていると判定された場合には、ステップS350にてISM P3b2が1D−LUTを参照して、当該画素についての色変換(分版)を実行する。すなわち、インデックスが格納された旨のフラグが立っている画素から、インデックスを取得するとともに、1D−LUTにて当該インデックスに対応付けられているインク量セットφを取得する。ステップS340またはステップS350のいずれかにおいて、当該画素についてのインク量セットφが取得できると、ステップS360においてすべての画素についてインク量セットφが取得できたか否かを判定する。ここでインク量セットφが未取得の画素が残っている場合には、ステップS320に戻り次の画素を選択する。
以上の処理を繰り返して実行することにより、すべての画素についてインク量セットφを取得することができる。すべての画素についてインク量セットφが取得できると、すべての画素がインク量セットφで表現された印刷データPDに変換されたこととなる。以上のように各画素について1D−LUTと3D−LUTのいずれを使用するかを判定することにより、インデックスが格納された枠F1〜F12に対応する画素については、各光源において各ターゲットTG1〜TG12に近い色が再現可能なインク量セットφを取得することができるとともに、RGB値が格納された画素については3D−LUTの作成指針(例えば、粒状性を重視する等。)に基づく色再現が可能なインク量セットφを取得させることができる。
ステップS370においては、各画素がインク量セットφで表現された印刷データPDをHTM P3b4が取得し、ハーフトーン処理を実行する。HTM P3b4はハーフトーン処理をするにあたっては公知のディザ法や誤差拡散法等を使用することができる。ハーフトーン処理が完了した印刷データPDにおいては、各画素が各インクを吐出させるか否を示す吐出信号を有している。ステップS380では、ハーフトーン処理が完了した印刷データPDをRTM P3b5が取得し、当該印刷データPDにおける吐出信号をプリンタ20が有する印刷ヘッドの各走査パスおよび各ノズルに割り振る処理を実行する。以上によりプリンタ20に出力可能な印刷制御データCDが生成でき、プリンタ20の制御に必要な信号を添付した印刷制御データCDをスプーラP1bおよびプリンタ20に出力する。これによりプリンタ20が印刷用紙上にインクを吐出して、見本チャートSCを形成する。
以上のようにして印刷用紙上に形成された見本チャートSCの枠FL1〜FL12に対応する領域において、各ターゲットTG1〜TG12のターゲット分光反射率Rt(λ)を再現することができる。すなわち、枠FL1〜FL12に対応する領域は、ターゲットTG1〜TG12の複数光源下での色に追従するようなインク量セットφで印刷されているため、各光源下においてターゲットTG1〜TG12と似たような色を再現することができる。例えば、見本チャートSCを室内で視認したときの各枠FL1〜FL12に対応する領域の色はターゲットTG1〜TG12を室内で視認したときの色を再現するし、見本チャートSCを室外で視認したときの各枠FL1〜FL12に対応する領域の色もターゲットTG1〜TG12を室外で視認したときの色を再現することができる。
なお、究極的には、ターゲットTG1〜TG12と完全に同じ分光反射率R(λ)を有する見本チャートSCを再現すれば、いかなる光源においてもターゲットTG1〜TG12と同様の色を再現することができる。しかしながら、プリンタ20が使用可能なインク(色材の種類)はCMYKlclmに限られているため、ターゲットTG1〜TG12と完全に同じ分光反射率R(λ)が再現可能なインク量セットφを求めることは実質的に不可能である。また、知覚色に影響しない波長域についてもターゲットTG1〜TG12と同様な分光反射率R(λ)が再現可能なインク量セットφを求めても、視覚的な再現精度の実現においては無駄となる。これに対して、本発明では、等色関数x(λ),y(λ),z(λ)に基づく重み付けを行った評価値E(φ)を利用してターゲット分光反射率Rt(λ)への近似性を評価しているため、視覚的に十分な精度が達成可能なインク量セットφを求めることができる。
一方、印刷用紙上に形成された見本チャートSCの枠FL1〜FL12に対応する領域においては、上述した3D−LUTに基づいたインク量セットφによって印刷がなされることとなる。そのため、当該領域についての印刷パフォーマンスは3D−LUTに基づくものとなる。上述したとおり本実施形態においては枠FL1〜FL12以外の領域は中間グレーの一様な画像を示すが、当該領域において3D−LUTが目標とする印刷パフォーマンスを満足させることができる。すなわち、再現色の階調性や、粒状性や、再現色の光源非依存性や、ガマットや、インクデューティが総合的に良好となる印刷を実現することができる。
4.キャリブレーション処理:
以上の処理によって印刷された見本チャートSCの枠FL1〜FL12においては、各ターゲットTG1〜TG12のターゲット分光反射率Rt(λ)を再現することができる。
しかしながら、場合によっては見本チャートSCの枠FL1〜FL12の実際の分光反射率R(λ)と、各ターゲットTG1〜TG12のターゲット分光反射率Rt(λ)との間に誤差が生じることもある。インク量セットφは、RPM P3a2が予測モデル(分光プリンティングモデル)を使用して予測したものであるため、分光プリンティングモデルを構築(分光反射率データベースRDBを作成)したプリンタと実際に印刷を行うプリンタ20の機体が異なったり、同一の機体であっても時期がずれたりすることによって誤差が生じることは避けられない。
そこで、本キャリブレーション処理においては、ターゲット分光反射率Rt(λ)の再現性をさらに向上すべく、見本チャートSCの枠FL1〜FL12が実際にターゲットTG1〜TG12と近いターゲット分光反射率Rt(λ)を再現しているか否かを確認する処理を行う。
図15,16は、キャリブレーション処理のフローチャートである。図2に示したようにLUG P3aは、キャリブレーション処理を行うためのモジュールとして、確認パッチ測定部(KPM)P3a5と補正ターゲット値取得部(MRA)P3a6とを備えている。
処理が開始されるとステップS400において、キャリブレーション処理の繰り返し回数を示すカウンタ値(n)を1にリセットする。
ステップS405においては、すでに印刷された見本チャートSCの枠FL1〜FL12について分光反射率R(λ)を測定する。ここでは、MDV P4がKPM P3a5の要求に応じて分光反射率計30を制御し、当該制御によって得られた分光反射率データRDをKPM P3a5が取得する。なお、このとき分光反射率R(λ)が測定される見本チャートSCの枠FL1〜FL12が本発明の確認パッチに相当する。また、確認パッチの測色により得られる分光反射率R(λ)を確認分光反射率Rc(λ)と記載する。上述した印刷制御処理によれば、理想的にはターゲットTG1〜TG12から測定したターゲット分光反射率Rt(λ)のそれぞれとステップS405にて測定した確認分光反射率Rc(λ)は同じとなる。しかしながら、上述したように誤差が生じることがあるため、ターゲット分光反射率Rt(λ)と確認分光反射率Rc(λ)は完全に同じとはならない。
図17は、ターゲットTG1(枠FL1)についてのターゲット分光反射率Rt(λ)と確認分光反射率Rc(λ)を対比して示している。同図に示すように、確認分光反射率Rc(λ)がターゲット分光反射率Rt(λ)におおむね追従できているものの、全体的に低反射率側に確認分光反射率Rc(λ)がシフトしている。例えば、プリンタ20が吐出する各インクのインク量が経時的に増加した場合には、全体的に低反射率側に確認分光反射率Rc(λ)がシフトすることとなる。
ステップS410においては、補正ターゲット値取得部(MRA)P3a6がターゲットTG1〜TG12(枠FL1〜FL12)を選択する。ステップS420においては、選択したターゲットTGについて、確認分光反射率Rc(λ)からターゲット分光反射率Rt(λ)を差し引くことにより、各波長の偏差ΔR(λ)を算出する。なお、ターゲット分光反射率Rt(λ)は、インデックステーブルIDBから得ることができる。
さらに、ステップS420においては、MRA P3a6がターゲット分光反射率Rt(λ)から前記偏差ΔR(λ)を引くことにより補正ターゲット分光反射率Rtm(λ)={Rt(λ)−ΔR(λ)}を算出する。
以上のようにして、補正ターゲット分光反射率Rtm(λ)が得られると、ステップS430aとステップS430bにおいて、ICM P3a1がRPM P3a2とECM P3a3を使用して、補正ターゲット分光反射率Rtm(λ)と同様の分光反射率R(λ)が再現可能なインク量セットを算出する処理を行う。ただし、前記ステップS230とは異なり、濃インクを優先的に使用して(淡インクの使用量を抑制して)インク量セットを算出する処理と、淡インクを優先的に使用して(濃インクの使用量を抑制して)インク量セットを算出する処理とを使い分けることにより、インク量セットφの最適解を算出する。そのために、インクセットをインク濃度に基づいて濃インクグループと淡インクグループとに分け、ステップS430aにおいては濃インクグループのインクが優先的に割り当てられるようにしてインク量セットを算出し、ステップS430bにおいては淡インクグループのインクが優先的に割り当てられるようにしてインク量セットを算出する。
ここで図18を参照して濃インクと淡インクの相違について説明する。同図は濃インクと淡インクの階調値変化に対する濃度変化を示したグラフである。濃インクと淡インクを比較すると、濃インクの濃度は淡インクの濃度に比べて相対的に高く、同じ量のインクを記録媒体へ付着したときに記録媒体に生じる色彩値の変動量は濃インクの方が大きい。例えば、各インクのインク濃度が256階調で表現されており、濃インクが淡インクの3倍の濃度であれば、濃インクの1階調分の濃度変化は、淡インクでは3階調分の濃度変化に相当する。具体例を挙げると、インクセットがCMYKlclmで構成されている場合は、CMYKインクが濃インクであり、lclmインクが淡インクである。無論、インクセットにライトイエローやライトブラック等を含む場合もあり、これらも淡インクになる。従って、インク打込量制限(単位面積あたりに付着可能なインク総量)を考慮すると、ある色を再現するためのインクセットを決定する際には、濃インクを優先的に使用した方がよいことになる。
しかし逆に言えば、淡インクは同じ量のインクを記録媒体へ付着したときに記録媒体に生じる濃度の変動量が相対的に小さいわけであるから、淡インクの方がきめ細かな濃度の違いを表現できる。濃インクでの1階調が淡インクでは3階調に相当するのであれば、淡インクは濃インクよりも1/3細かい濃度の違いを表せる。いわば、淡インクは濃インクに比べて濃度についての分解能が高いのである。以上の濃淡インクの特徴を鑑みると、ある色を再現するためのインクセットには、濃度についての分解能の低い濃インクを優先的に割り当てた後で、濃度についての分解能の高い淡インクで微調整を行うのが好ましいと言える。
なお、本実施形態においては、インクセットを濃インクグループと淡インクグループとにグループ分けしている。本実施形態の濃インクグループは、CインクとMインクとKインクで構成され、淡インクグループはYインクとlcインクとlmインクで構成される。Yインクは濃度の高いインクであるが、本実施形態においてはYインクを淡インクグループに入れてある。CMYKlclmのインクセットの場合、どんなに分割してもイエロー方向にはイエロー一個しか色がない。そこで、淡インクグループを構成するインクだけを変更しつつでインク量セット予測した場合にも、色相ズレが発生しにくくするためである。無論、インクセット中にly(ライトイエロー)等のように明度方向への変化量を発生しやすい淡インクがあれば、単純にインク濃度のみに着目して、CMYKで濃インクグループを構成し、lclmlyで淡インクグループを構成することも可能である。
さらに、本実施形態のインク量セットを算出する処理においては、プリンタ20にて印刷用紙上に吐出可能な最小インク量よりも少ない量ずつ変化させながら最適解を探索する。例えば、プリンタ20がインク吐出量を変化させることにより256階調を表現可能であれば、0.01階調ずつ変化させながら最適なインク量セットを探索する。このように探索ステップを小さくすることにより、最適解近傍での振動が抑制され、最適解へ到達しやすくしている。しかしながら、実際に設定される前に小数点以下のインク量は端数処理されるため、丸め誤差が発生する。後述する本実施形態においては、この丸め誤差の影響も極小化されることになる。なお、以下の説明においては、吐出可能な最小インク量を単位量とする数値を整数値、該単位量よりも小さい数値を小数値と記載することにする。
ステップS430aにおいては、濃インクセットを優先的に変化させる最適化処理を実行する。より具体的には、以下の目的関数を用いることにより、濃インクセットを優先的に変化させる最適化処理を実現する。
Figure 2010136022

すなわち、上記の(4)式で表された評価値E(φ)のターゲット分光反射率Rt(λ)を補正ターゲット分光反射率Rtm(λ)に置き換えた関数を目的関数とし、当該目的関数を極小化させるようなインク量セットφの最適解を算出する。この目的関数においては、インクグループ毎に重み付けなどを行っているわけではないが、最適解の近傍に至るまでは、濃インクと淡インクとを同じ量だけ変化させたときに目的関数をより大きく減少させる濃インクの方が優先的に利用される。無論、淡インクセットのインク量変化を妨げる項を目的関数に付け加えて(例えば、濃インクセットのインク量が変化すると減少する項を目的関数に付け加えたり、淡インクセットのインク量が変化すると増大する項を目的関数に付け加えたりする。)も構わない。インク量セットφの最適解が算出されると、このインク量セットφについて端数処理を行う。
端数処理では小数点以下の数値が丸められて、端数処理後のインク量セットは最適解からの丸め誤差が生じる。濃インクのインク量と淡インクのインク量とで同じ量の丸め誤差とが発生した場合、濃インクの方が再現される色に対する影響が大きい。例えば、濃インクで0.5階調が丸められた場合の濃度変動は、淡インクの1.5階調分の濃度変動にあたる。すなわち、濃インクで小数値として丸められた濃度変動は、淡インクの整数値として表現可能なことがある。よって、より最適解に近いインク量セットを得るために、この誤差分を淡インクのインク量に積極的に表現する。
そこで、ステップS430bにおいては、ステップS430aで算出された端数処理後のインク量セットを初期条件として淡インクセットを優先的に変化させる最適化処理を行う。無論、優先するのみならず、濃インクセットが完全に変化しないようにして淡インクセットのみが変化するようにしてもよい。より具体的には、以下の目的関数を用いることにより、淡インクセットを優先的に変化させる最適化処理を実現する。
Figure 2010136022
前記(5)式において、ΔC,ΔM,ΔKはそれぞれ最適化処理において目的関数が減少するか否かの判定時に、試験的に変化されたシアンインク,マゼンタインク、ブラックインクの量である。また、w,w,wは、それぞれシアン,マゼンタ,ブラックの重み関数である。前記(5)式の目的関数を用いることにより、CMKインクの試験的な変化に対しては目的関数が増大し、その変化を妨げることになる。すなわち目的関数の重み項に起因した目的関数の増分を打ち消す以上に目的関数の初項が減少しない限り、CMYインクを変化させる最適化が妨げられる。よって、濃インクセットよりも淡インクセットが優先的に変化される最適化処理が実現される。このようにして算出されたインク量セットφの最適解について端数処理を行う。この端数処理においては、主に淡インクセットにおいて丸め誤差が発生する。従って、端数処理の前後で発生する丸め誤差は、ステップS430で発生する丸め誤差よりも少ない。
ステップS440においては、LUT出力モジュール(LOM)P3a4が、1D−LUTにおいて対応するインデックスについてのインク量セットφを、当該最適化されたインク量セットφによって更新する。インク量セットφが更新できると、ステップS450において、すべてのターゲットTG1〜TG12(枠FL1〜FL12)が選択されたか否かを判定し、選択していない場合にはステップS420にて次のターゲットTG1〜TG12(枠FL1〜FL12)を選択する。これにより、すべてのターゲットTG1〜TG12について、インク量セットφを更新することができる。以上のように1D−LUTを更新することにより、以降において実行する印刷制御データ生成処理においては、更新されたインク量セットφに基づいて見本チャートSCの印刷を実行させることができる。
以上のようなキャリブレーション処理を実行することにより、より高精度の分光反射率R(λ)の再現を実現することができる。例えば、確認分光反射率Rc(λ)がターゲット分光反射率Rt(λ)よりも大きかった場合、確認分光反射率Rc(λ)とターゲット分光反射率Rt(λ)の偏差ΔR(λ)がもとのターゲット分光反射率Rt(λ)から差し引かれるため、補正ターゲット分光反射率Rtm(λ)はもとのターゲット分光反射率Rt(λ)よりも小さい値となる。従って、補正ターゲット分光反射率Rtm(λ)によって最適化されたインク量セットφによれば、再現される分光反射率R(λ)を偏差ΔR(λ)の大きさに応じて下方修正することができる。逆に、確認分光反射率Rc(λ)がターゲット分光反射率Rt(λ)よりも小さかった場合、補正ターゲット分光反射率Rtm(λ)がもとのターゲット分光反射率Rt(λ)よりも大きい値とされるため、偏差ΔR(λ)の大きさに応じて上方修正することができる。
さらに、本実施形態では、以上のキャリブレーション処理を繰り返して実行することにより、より高精度の分光反射率R(λ)の再現を実現する。ステップS460においては、キャリブレーション処理の繰り返し回数を示すカウンタnが3になったか否かを判定し、3になっていない場合には、カウンタ値nに1を加算して(ステップS470)、ステップS402に戻る。これにより、ステップS402における確認パッチの印刷が再度実行されることとなる。ここでは、初回のキャリブレーション処理によって更新されたインク量セットφに基づいて確認パッチの印刷を印刷させるため、前回よりもターゲット分光反射率Rt(λ)と確認分光反射率Rc(λ)との偏差ΔR(λ)の絶対値が減少していると予測される。ステップS420においては新たな確認分光反射率Rc(λ)について補正ターゲット分光反射率Rtm(λ)={Rt(λ)−ΔR(λ)}が設定され、ステップS430〜S440では減少した偏差ΔR(λ)をさらに打ち消すようなインク量セットφに更新することができる。このようなキャリブレーション処理の繰り返しは、カウンタ値nが3となるまで繰り返されるため、その間に偏差ΔR(λ)の絶対値を極めて小さくさせることができ、より高精度の分光反射率再現を実現することができる。
なお、以上においては、偏差ΔR(λ)そのものを、もとのターゲット分光反射率Rt(λ)から差し引くようにしたが、偏差ΔR(λ)の80%程度を差し引くようにしてもよい。むろん、繰り返し回数は3回に限られるものではない。以上のようなキャリブレーション処理は、同一機体のプリンタ20を長期間使用しなかったときや、異なる機体のプリンタに見本チャートSCの印刷をさせる場合に実行しておくのが望ましい。
以上、ステップS400〜S460のキャリブレーション処理においては、確認パッチの印刷と確認パッチの測色とを実行しつつ最適化処理を2度実行してインク量セットφを求めたが、さらに演算精度を向上するために、パッチ印刷とパッチ測色を伴わない最適化処理を複数回繰り返してもよい。そこで、ステップS480〜S530においては、直前に測色した結果(ステップS460が条件成立した時のループにおいてステップS405で測色した結果)を利用して、キャリブレーション精度を向上するための最適化処理を実行する。
ステップS460で条件成立するとステップS480でカウンタ値nに1を加算し、ステップS490に進む。
ステップS490においては、補正ターゲット値取得部(MRA)P3a6がターゲットTG1〜TG12(枠FL1〜FL12)を選択する。これはステップS420と同様である。
ステップS500においては、カウンタnがn−1のときに更新された1D−LUTのインク量セットφn−1と、カウンタnがn−2のときに更新された1D−LUTのインク量セットφn−2とを比較し、濃インクグループのインク量に変化があったか否かを判断する。変化があった場合は、条件成立としてステップS501に進み濃インクグループを優先した最適化と淡インクグループを優先した最適化とを順に実行し、変化が無かった場合は、条件不成立としてステップS505に進んで淡インクグループを優先した最適化だけを行う。分解能の低い濃インクグループのインク量セットが2回連続で同値に最適化されていれば、その値が最適解であると考えられるからである。なお、インク量セットφn−1とインク量セットφn−2とを比較するために、RAMには直近2回分の1D−LUTが一時的に記憶されている。
ステップS501においては、カウンタ値がn−1のときに作成された1D−LUTにおけるインク量セットを初期値としてステップS430aと同様の最適化処理を行い、ステップS502においては、ステップS430bと同様の最適化処理を行う。また。ステップS505においても、カウンタ値がn−1のときに作成された1D−LUTにおけるインク量セットを初期値としてステップS430bと同様の最適化処理を行う。
ステップS510においては、LUT出力モジュール(LOM)P3a4が、ステップS502もしくはS505において最適化されたインク量セットφで、1D−LUTにおいて対応するインデックスについてのインク量セットφを更新する
ステップS520においては、すべてのターゲットTG1〜TG12(枠FL1〜FL12)が選択されたか否かを判定し、選択していない場合にはステップS490にて次のターゲットTG1〜TG12(枠FL1〜FL12)を選択する。これにより、すべてのターゲットTG1〜TG12について、インク量セットφを更新することができる。
ステップS530においては、キャリブレーション処理の繰り返し回数を示すカウンタnがm(mは4以上の整数)になったか否かを判定し、mになっていない場合には、カウンタ値nに1を加算して(ステップS480)、ステップS490以降の処理を繰り返す。カウンタ値がmに達していれば、所定のループ回数の最適化が終了したと判断して、キャリブレーション処理を終了する。なお、このループ回数は、例えば、濃インクグループのインク量が変化しなくなってから所定回数行う、等のように設定しても構わない。
5.分光プリンティングモデル
図19は、本実施形態のプリンタ20の印刷方式を模式的に示している。同図において、プリンタ20は、CMYKlclmのインクごとに複数のノズル21a,21a・・・を備えた印刷ヘッド21を備えており、ノズル21a,21a・・・が吐出するCMYKlclmのインクごとのインク量を上述したインク量セットφ(dc,dm,dy,dk,dlc,dlm)によって指定された量とする制御が印刷制御データCDに基づいて行われる。各ノズル21a,21a・・・が吐出したインク滴は印刷用紙上において微細なドットとなり、多数のドットの集まりによってインク量セットφ(dc,dm,dy,dk,dlc,dlm)に応じたインク被覆率の印刷画像が印刷用紙上に形成されることとなる。
RPM P3a2が使用する予測モデル(分光プリンティングモデル)は、本実施形態のプリンタ20で使用され得る任意のインク量セットφ(dc,dm,dy,dk,dlc,dlm)で印刷を行った場合の分光反射率R(λ)を予測分光反射率Rs(λ)として予測するための予測モデルである。分光プリンティングモデルにおいては、インク量空間における複数の代表点について実際にカラーパッチを印刷し、その分光反射率R(λ)を分光反射率計によって測定することにより得られた分光反射率データベースRDBを用意する。そして、この分光反射率データベースRDBを使用したセル分割ユール・ニールセン分光ノイゲバウアモデル(Cellular Yule-Nielsen Spectral Neugebauer Model)による予測を行うことにより、正確に任意のインク量セットφ(dc,dm,dy,dk,dlc,dlm)で印刷を行った場合の分光反射率R(λ)を予測する。
図20は、分光反射率データベースRDBを示している。同図に示すように分光反射率データベースRDBはインク量空間(本実施形態では6次元であるが、図の簡略化のためCM面のみ図示。)における複数の格子点のインク量セットφ(dc,dm,dy,dk,dlc,dlm)について実際に印刷/測定をして得られた分光反射率R(λ)が記述されたルックアップテーブルとなっている。例えば、各インク量軸を分割する5グリッドの格子点を発生させる。ここでは513個もの格子点が発生し、膨大な量のカラーパッチの印刷/測定をすることが必要となるが、実際にはプリンタ20にて同時に搭載可能なインク数や同時に吐出可能なインクデューティの制限があるため、印刷/測定をする格子点の数は絞られることとなる。
さらに、一部の格子点のみ実際に印刷/測定をし、他の格子点については実際に印刷/測定を行った格子点の分光反射率R(λ)に基づいて分光反射率R(λ)を予測することにより、実際に印刷/測定を行うカラーパッチの個数を低減させてもよい。分光反射率データベースRDBは、プリンタ20が印刷可能な印刷用紙ごとに用意されている必要がある。厳密には、分光反射率R(λ)は印刷用紙上に形成されたインク膜(ドット)による分光透過率と印刷用紙の反射率によって決まるものであり、印刷用紙の表面物性(ドット形状が依存)や反射率の影響を大きく受けるからである。次に、分光反射率データベースRDBを使用したセル分割ユール・ニールセン分光ノイゲバウアモデルによる予測を説明する。
RPM P3a2は、ICM P3a1の要請に応じて分光反射率データベースRDBを使用したセル分割ユール・ニールセン分光ノイゲバウアモデルによる予測を実行する。この予測にあたっては、ICM P3a1から予測条件を取得し、この予測条件を設定する。具体的には、印刷用紙やインク量セットφを印刷条件として設定する。例えば、光沢紙を印刷用紙として予測を行う場合には、光沢紙にカラーパッチを印刷することにより作成した分光反射率データベースRDBが設定される。
分光反射率データベースRDBの設定ができると、ICM P3a1から入力されたインク量セットφ(dc,dm,dy,dk,dlc,dlm)を分光プリンティングモデルに適用する。セル分割ユール・ニールセン分光ノイゲバウアモデルは、よく知られた分光ノイゲバウアモデルとユール・ニールセンモデルとに基づいている。なお、以下の説明では、説明の簡略化のためCMYの3種類のインクを用いた場合のモデルについて説明するが、同様のモデルを本実施形態のCMYKlclmを含む任意のインクセットを用いたモデルに拡張することは容易である。また、セル分割ユール・ニールセン分光ノイゲバウアモデルについては、Color Res Appl 25, 4-19, 2000およびR Balasubramanian, Optimization of the spectral Neugebauer model for printer characterization, J. Electronic Imaging 8(2), 156-166 (1999)を参照。
図21は、分光ノイゲバウアモデルを示す図である。分光ノイゲバウアモデルでは、任意のインク量セット(dc,dm,dy)で印刷したときの印刷物の予測分光反射率Rs(λ)は、以下の(6)式で与えられる。
Figure 2010136022
ここで、aiはi番目の領域の面積率であり、Ri(λ)はi番目の領域の分光反射率である。添え字iは、インクの無い領域(w)と、シアンインクのみの領域(c)と、マゼンタインクのみの領域(m)と、イエローインクのみの領域(y)と、マゼンタインクとイエローインクが吐出される領域(r)と、イエローインクとシアンインクが吐出される領域(g)と、シアンインクとマゼンタインクが吐出される領域(b)と、CMYの3つのインクが吐出される領域(k)をそれぞれ意味している。また、fc,fm,fyは、CMY各インクを1種類のみ吐出したときにそのインクで覆われる面積の割合(「インク被覆率(Ink area coverage)」と呼ぶ)である。
インク被覆率fc,fm,fyは、図21(B)に示すマーレイ・デービスモデルで与えられる。マーレイ・デービスモデルでは、例えばシアンインクのインク被覆率fcは、シアンのインク量dcの非線形関数であり、例えば1次元ルックアップテーブルによってインク量dcをインク被覆率fcに換算することができる。インク被覆率fc,fm,fyがインク量dc,dm,dyの非線形関数となる理由は、単位面積に少量のインクが吐出された場合にはインクが十分に広がるが、多量のインクが吐出された場合にはインクが重なり合うためにインクで覆われる面積があまり増加しないためである。他の種類のMYインクについても同様である。
分光反射率に関するユール・ニールセンモデルを適用すると、前記(6)式は以下の(7a)式または(7b)式に書き換えられる。
Figure 2010136022

ここで、nは1以上の所定の係数であり、例えばn=10に設定することができる。前記の(7a)式および(7b)式は、ユール・ニールセン分光ノイゲバウアモデル(Yule-Nielsen Spectral Neugebauer Model)を表す式である。
本実施形態で採用するセル分割ユール・ニールセン分光ノイゲバウアモデル(Cellular Yule-Nielsen Spectral Neugebauer Model)は、上述したユール・ニールセン分光ノイゲバウアモデルのインク量空間を複数のセルに分割したものである。
図22(A)は、セル分割ユール・ニールセン分光ノイゲバウアモデルにおけるセル分割の例を示している。ここでは、説明の簡略化のために、CMインクのインク量dc,dmの2つの軸を含む2次元インク量空間でのセル分割を描いている。なお、インク被覆率fc,fmは上述したマーレイ・デービスモデルにてインク量dc,dmと一意の関係にあるため、インク被覆率fc,fmを示す軸と考えることもできる。白丸は、セル分割のグリッド点(「格子点」と呼ぶ)であり、2次元のインク量(被覆率)空間が9つのセルC1〜C9に分割されている。各格子点に対応するインク量セット(dc,dm)は、分光反射率データベースRDBに規定された格子点に対応するインク量セットとされている。すなわち、上述した分光反射率データベースRDBを参照することにより、各格子点の分光反射率R(λ)を得ることができる。従って、各格子点の分光反射率R(λ)00,R(λ)10,R(λ)20・・・R(λ)33は、分光反射率データベースRDBから取得することができる。
実際には、本実施形態ではセル分割もCMYKlclmの6次元インク量空間で行うとともに、各格子点の座標も6次元のインク量セットφ(dc,dm,dy,dk,dlc,dlm)によって表される。そして、各格子点のインク量セットφ(dc,dm,dy,dk,dlc,dlm)に対応する格子点の分光反射率R(λ)が分光反射率データベースRDB(例えば光沢紙のもの)から取得されることとなる。
図22(B)は、セル分割モデルにて使用するインク被覆率fcとインク量dcとの関係を示している。ここでは、1種類のインクのインク量の範囲0〜dcmaxも3つの区間に分割されており、各区間毎に0から1まで単調に増加する非線形の曲線によってセル分割モデルにて使用する仮想的なインク被覆率fcが求められる。他のインクについても同様にインク被覆率fm,fyが求められる。
図22(C)は、図22(A)の中央のセルC5内にある任意のインク量セット(dc,dm)にて印刷を行った場合の予測分光反射率Rs(λ)の算出方法を示している。インク量セット(dc,dm)にて印刷を行った場合の分光反射率Rs(λ)は、以下の(8)式で与えられる。
Figure 2010136022

ここで、(8)式におけるインク被覆率fc,fmは図22(B)のグラフで与えられる値である。また、セルC5を囲む4つの格子点に対応する分光反射率R(λ)11,(λ)12,(λ)21,(λ)22は分光反射率データベースRDBを参照することにより取得することができる。これにより、(8)式の右辺を構成するすべての値を確定することができ、その計算結果として任意のインク量セットφ(dc,dm)にて印刷を行った場合の予測分光反射率Rs(λ)を算出することができる。波長λを可視波長域にて順次シフトさせていくことにより、可視波長域における予測分光反射率Rs(λ)を得ることができる。インク量空間を複数のセルに分割すれば、分割しない場合に比べて予測分光反射率Rs(λ)をより精度良く算出することができる。以上のようにして、RPM P3a2がICM P3a1の要請に応じて予測分光反射率Rs(λ)を予測することができる。
6.変形例
6−1:変形例1
図23は、変形例においてECM P3a3が設定する重み関数w(λ)を模式的に示している。同図において、ターゲットTGから得られたターゲット分光反射率Rt(λ)が示されており、当該ターゲット分光反射率Rt(λ)と各等色関数x(λ),y(λ),z(λ)との相関係数cx,cy,czがECM P3a3によって算出されている。そして、下記の(9)式によって本変形にかかる重み関数w(λ)を算出する。
Figure 2010136022

前記の(9)式においては、ターゲットTGから得られたターゲット分光反射率Rt(λ)との相関が高い等色関数x(λ),y(λ),z(λ)ほど線形結合の際の重みが大きくなるようにされている。以上のようにして得られた重み関数w(λ)においては、ターゲットTGのターゲット分光反射率Rt(λ)が大きい波長域についての重みを強調することができる。従って、各光源下での反射光の分光エネルギーのスペクトルが強くなりがちな波長域を重視した評価値E(φ)を得ることができる。すなわち、特にターゲットTGのターゲット分光反射率Rt(λ)が大きい波長域については、ターゲットTGのターゲット分光反射率Rt(λ)と予測分光反射率Rs(λ)とのずれを許容しないようなインク量セットφの最適解を得ることができる。むろん、重み関数w(λ)は各等色関数x(λ),y(λ),z(λ)に由来しているため、人間の知覚に適合した評価値E(φ)を得ることができる。
6−2:変形例2
図24は、別の変形例においてECM P3a3が設定する重み関数w(λ)を模式的に示している。同図において、ターゲットTGから得られたターゲット分光反射率Rt(λ)をそのまま重み関数w(λ)として適用している。このようにすることによっても、特にターゲットTGのターゲット分光反射率Rt(λ)が大きい波長域についてターゲットTGの分光反射率R(λ)とターゲット分光反射率Rt(λ)とのずれを許容しないようなインク量セットφの最適解を得ることができる。
6−3:変形例3
図25は、別の変形例においてECM P3a3が設定する重み関数w(λ)を模式的に示している。同図において、5種類(標準昼光系のD50光源,D55光源,D65光源、白熱電球系のA光源、蛍光ランプ系のF11光源)の各光源の分光エネルギーP D50(λ),P D55(λ),P D65(λ),PA(λ),P F11(λ)が示されている。本変形例においては、下記の(10)式によって、これらの分光エネルギーP D50(λ),P D55(λ),P D65(λ),PA(λ),P F11(λ)を線形結合することにより重み関数w(λ)を算出する。
Figure 2010136022
前記の(10)式において、w1〜w5は各光源についての重みを設定する重み係数である。このように、光源の分光エネルギー分布P D50(λ),P D55(λ),P D65(λ),PA(λ),P F11(λ)に由来する重み関数w(λ)を設定することにより、各光源下での反射光の分光エネルギーのスペクトルが強くなりがちな波長域を重視した評価値E(φ)を得ることができる。また、重み係数w1〜w5を調整することも可能である。例えば全光源における色の再現性をバランスよく確保したい場合にはw1=w2=w3=w4=w5とすればよいし、人工光源における色の再現性を重視したい場合にはw1,w2,w3<w4,w5とすればよい。
6−4:変形例4
図26は、変形例においてディスプレイ40に表示されるUI画面を示している。同図において、UI画面において複数のターゲット分光反射率Rt(λ)のグラフが表示されている。このようなUI画面を表示させることにより、ユーザーがステップS140においてターゲットTGのターゲット分光反射率Rt(λ)を測定する代わりに、所望の波形のグラフをターゲットTGのターゲット分光反射率Rt(λ)として選択することができる。このようにすることにより、実際に分光反射率測定をしなくてもターゲット分光反射率Rt(λ)を設定することができる。むろん、グラフの波形をユーザーが直接編集できるようにしてもよい。例えば、新規の物体表面の開発を行う際に目標とするターゲット分光反射率Rt(λ)に編集しておけば、実際に物体表面を試作することなく、目標とするターゲット分光反射率Rt(λ)を有する見本チャートSCをプリンタ20によって印刷させることができる。
6−5:変形例5
図27は、変形例にかかる評価値(φ)を模式的に説明している。同図において、ターゲットTGのターゲット分光反射率Rt(λ)に対して前記5種類の光源を照射したときの色彩値(ターゲット色彩値)を上述した(1)式,図5によって算出する。一方、RPM P3a2が予測した予測分光反射率Rs(λ)に対して前記5種類の光源を照射したときの色彩値(予測色彩値)も上述した(1)式(Rt(λ)をRs(λ)に置き換えて使用),図5によって算出する。そして、各光源におけるターゲット色彩値と予測色彩値の色差ΔE(ΔE2000)をCIE DE2000の色差式に基づいて算出する。そして、各光源についての色差ΔEをΔE D50,ΔE D55,ΔE D65,ΔEA,ΔE F11とし、下記の(11)式によって評価値E(φ)を算出する。
Figure 2010136022

前記の(11)式において、w1〜w5は各光源の重みを設定する重み係数であり、上述した変形例3の重み係数w1〜w5とほぼ同様の性質を有する。ここでも全光源における色の再現性をバランスよく確保したい場合にはw1=w2=w3=w4=w5とすればよいし、人工光源における色の再現性を重視したい場合にはw1,w2,w3<w4,w5とすればよい。
本変形例においてキャリブレーションを行う場合、確認パッチとしての見本チャートSCを印刷し、その分光反射率R(λ)を確認分光反射率Rc(λ)として測定する。そして、ターゲットTGのターゲット分光反射率Rt(λ)に対して前記5種類の光源を照射したときのターゲット色彩値を上述した(1)式,図5によって算出するとともに、確認パッチに前記5種類の光源を照射したときの色彩値を、上述した(1)式(Rt(λ)をRc(λ)に置き換えて使用),図5によって算出する。なお、後者の色彩値を確認色彩値と表記する。そして、各光源について、確認色彩値からターゲット色彩値の偏差(CIELAB色空間の偏差ベクトル)を算出し、当該偏差をターゲット色彩値から引く(偏差ベクトルの逆向きのベクトルを加算する)ことにより、補正ターゲット色彩値を算出する。なお、ターゲットTGに実際に前記5種類の光源を照射しつつ測色を行うことにより、確認色彩値を直接得るようにしてもよい。
図28は、補正ターゲット色彩値を模式的に示している。同図において、一例としてD50光源におけるターゲット色彩値(L* t,a* t,b* t)と確認色彩値(L* c,a* c,b* c)が図示されており、これらの偏差ベクトルdf(ΔL*,Δa*,Δb*)に基づいて補正ターゲット色彩値(L* tm,a* tm,b* tm)が算出される様子がCIELAB空間にて図示されている。以上のようにして、補正ターゲット色彩値(L* tm,a* tm,b* tm)が算出できると、ステップS430において、上記の(11)式の評価値E(φ)を極小化させるインク量セットφを算出する。なお、上記の(11)式の評価値E(φ)において、もとのターゲット色彩値と予測色彩値との色差ΔEをΔE D50,ΔE D55,ΔE D65,ΔEA,ΔE F11として使用するのではなく、補正後の補正ターゲット色彩値と予測色彩値との色差ΔEをΔE D50,ΔE D55,ΔE D65,ΔEA,ΔE F11として使用する。なお、本変形例においては、状態値として色彩値のみを使用するため、必ずしも分光反射率R(λ)を取得する必要はない。従って、ターゲットTGについて最初から複数光源下での色彩値を測色等によって取得するようにしてもよい。
ところで、本変形例において、各光源についてのターゲット色彩値(L* t,a* t,b* t)と確認色彩値(L* c,a* c,b* c)が得られた時点で、これらの色差ΔE(ΔE2000)を各光源ごとに算出するようにしてもよい。なお、ここでの色差ΔEを各光源についてΔe D50,Δe D55,Δe D65,ΔeA,Δe F11と表記する。当該色差Δe D50,Δe D55,Δe D65,ΔeA,Δe F11によれば、見本チャートSCがどの程度精度よく再現されているかを色差ΔE2000によって把握することができる。また、光源について色差Δe D50,Δe D55,Δe D65,ΔeA,Δe F11を下記の(12)式のように平均した平均色差Δeによれば、複数光源の各ターゲットTGの再現精度を総合的に判断することができる。
Figure 2010136022
図29は、本変形例にかかるキャリブレーション処理の流れを示している。ここでは、見本チャートSCが印刷(ステップS300)されると、ステップS405において各確認パッチ(枠FL1〜FL12)の確認分光反射率Rc(λ)を測定する。そして、ステップS402においては、ターゲット色彩値(L* t,a* t,b* t)と確認色彩値(L* c,a* c,b* c)との平均色差Δeを各枠FL1〜FL12について算出する。そして、すべての枠FL1〜FL12についての平均色差Δeが所定の閾値Th(例えばΔE=1.0。)を超えているか否かをステップS404にて判定し、いずれかの枠FL1〜FL12の平均色差Δeが閾値Th超えている場合にはステップS410以降のキャリブレーション処理を実行する。キャリブレーション処理が完了すると再度ステップS300に戻り、更新された1D−LUTに基づいて見本チャートSCを再度印刷し、同様の処理を繰り返し実行する。このようにすることにより、平均色差Δeが閾値Thを満足するまで、キャリブレーション処理を繰り返すことができる。
図25に図示したように各光源はそれぞれ異なる分光エネルギーのスペクトルを有しているため、各光源の色差Δe D50,Δe D55,Δe D65,ΔeA,Δe F11が一様に大きくなったり小さくなったりするとは限らない。例えば、昼光系の色差Δe D50,Δe D55,Δe D65は大きくても、白熱電球系の色差ΔeAが小さい場合も考えられる。このような場合、色差ΔeAを小さく維持したまま、色差Δe D50,Δe D55,Δe D65を小さくするようなキャリブレーション処理を行うのが望ましい。そこで、本変形例においては、下記の(12)式の評価値E(φ)を使用して、ステップS430の最適化を行うものとする。
Figure 2010136022
前記の(12)式において、Δr(λ)はステップS230にて最適化されたインク量セットφによる予測分光反射率Rs(λ)と、ステップS430のキャリブレーションにて最適化されるインク量セットφによる予測分光反射率Rs(λ)との差分の絶対値を表している。また、w(λ)は各波長ごと重みを定義した重み関数を示している。
図30は、重み関数w(λ)の一例を示している。同図において、最も色差ΔeAが小さかったA光源の分光エネルギーのスペクトルとほぼ同じ傾向を示すように設定されている。なお、所定の値より分光エネルギーが小さい波長域については、重み関数w(λ)=0としている。このようにすることにより、A光源での色彩値に大きく寄与する波長域における分光反射率の変化に応じて評価値E(φ)を増大させることができる。すなわち、キャリブレーション処理における長波長領域での分光反射率の変化を制限し、結果的にA光源での色彩値をできるだけ変化しないようにしている。このようにすることにより、A光源で色差ΔeAを小さく維持しつつ、他の光源の色差Δe D50,Δe D55,Δe D65,Δe F11を減少させることができる。
6−6:変形例6
なお、上述した実施形態において選択されていない枠Fに対応する領域については、枠F以外の領域と同じ色で印刷を行うようにすればよい。むろん、選択されていない枠Fに対応する領域については、分光的な再現性を要求する必要がないため、枠F以外の領域と同様に3D−LUTを使用した色変換を行わせるようにすればよい。さらに、ターゲットTGが指定された枠Fに対応する領域以外において、模様や文字やマーク等を印刷するようにしてもよい。例えば、ターゲットTGが指定された枠Fの付近に、ターゲットTGがどのようなものであるかを示す文字が記載できるようにしてもよい。
6−7:変形例7
図31は、1D−LUT生成処理の変形例のフローチャートである。同図に示す1D−LUT生成処理は、ステップS230を除いて図9の1D−LUT生成処理と同様の処理である。本変形例においては、評価値を最小化する処理として、図15,16に示したキャリブレーション処理において行ったものと同様に、濃インクを優先的に使用して(淡インクの使用量を抑制して)インク量セットを算出する処理と、淡インクを優先的に使用して(濃インクの使用量を抑制して)インク量セットを算出する処理とを使い分けることにより、インク量セットφの最適解を算出するように構成してある。このように1D−LUTを作成する段階で、濃インク量と淡インク量の使い分けを最適化しておくことにより、キャリブレーション処理における演算量を少なく出来る。
6−8.変形例8
図32〜図33は、本発明の変形例にかかる印刷システムのソフトウェア構成を示している。図32に示すように、上述した実施形態のLUG P3aに相当する構成がPDV P3bの内部モジュール(1D−LUT作成部)として備えられていてもよい。また、図33に示すように、上述した実施形態のLUG P3aに相当する構成が他のコンピュータ110において実行されてもよい。この場合、コンピュータ10とコンピュータ110とが所定の通信インターフェイスCIFによって接続され、コンピュータ110のLUG P3aにて生成された1D―LUTが通信インターフェイスCIFを介してコンピュータ10に送信される。通信インターフェイスCIFは、インターネットを介在するものであってもよい。その場合、コンピュータ10はインターネット上のコンピュータ110から取得した1D−LUTを参照して色変換を行うことができる。さらに、プリンタ20が各ソフトウェア構成P1〜P5を実行するようにしてもよい。むろん、ソフトウェア構成P1〜P5と同等の処理を実行するハードウェアがプリンタ20に組み込まれる場合にも、本発明を実現することができる。
印刷制御装置のハードウェア構成を示すブロック図である。 印刷制御装置のソフトウェア構成を示すブロック図である。 印刷データ生成処理のフローチャートである。 UI画面の一例を示す図である。 分光反射率に基づいて色彩値を算出するための計算を説明する図である。 印刷データを示す図である。 インデックステーブルを示す図である。 印刷制御処理の全体的な流れを示すフローチャートである。 1D−LUT生成処理のフローチャートである。 インク量セットを最適化する処理の流れを示す模式図である。 インク量セットが最適化されていく様子を示す模式図である。 1D−LUTを示す図である。 印刷制御データ生成処理のフローチャートである。 3D−LUTを示す図である。 キャリブレーション処理のフローチャートである。 キャリブレーション処理のフローチャートである。 偏差を説明するグラフである。 所定色相において、各インクの単位量当りの色変化を示した概念図である。 プリンタの印刷方式を示す模式図である。 分光反射率データベースを示す図である。 分光ノイゲバウアモデルを示す図である。 セル分割ユール・ニールセン分光ノイゲバウアモデルを示す図である。 変形例にかかる重み関数を示す模式図である。 変形例にかかる重み関数を示す模式図である。 変形例にかかる重み関数を示す模式図である。 変形例にかかるUI画面を示す図である。 変形例にかかる評価値を示す模式図である。 変形例にかかる補正ターゲット色彩値を説明している。 変形例にかかるキャリブレーション処理のフローチャートである。 変形例にかかる重み関数を説明するグラフである。 変形例にかかる1D−LUT生成処理のフローチャートである。 変形例にかかる印刷システムのソフトウェア構成である。 変形例にかかる印刷システムのソフトウェア構成である。
符号の説明
10…コンピュータ、11…CPU、12…RAM、13…ROM、14…HDD、15…GIF、16…VIF、17…IIF、18…バス、P1…OS、P1a…GDI、P1b…スプーラ、P2…APL、P2a…UIM、P2b…MCM、P2c…PDG、P3…PDV、P3a1…ICM、P3a2…RPM、P3a3…ECM、P3a4…LOM、P3a5…KPM、P3a6…MRA、P4…MDV、P5…DDV。

Claims (9)

  1. 少なくとも1つの色相に関して略同一の色相を有し濃度が互いに異なる濃色材と淡色材を含んで構成された複数の色材を印刷装置にて記録媒体に付着させて印刷を実行させるにあたり、前記色材の使用量の組合せである色材量セットを前記印刷装置に指定し、前記色材量セットに基づく印刷を実行させる印刷制御装置であって、
    前記色材量セットとインデックスとの対応関係を規定したルックアップテーブルを参照することにより、指定されたインデックスに対応する前記色材量セットを前記印刷装置に指定して印刷させる印刷手段を備え、
    前記ルックアップテーブルにおいては、
    対象物の色を示す情報であるターゲット値を特定する前記インデックスと、前記印刷装置にて前記記録媒体上に付着したときに前記ターゲット値に対する近似性が極大化する色材量セットであるターゲット色材量セットと、の対応関係が規定されており、
    前記ターゲット色材量セットは、前記淡色材の使用量を抑制しつつ前記近似性を極大化するように所定の予測モデルに基づいて予測された第1色材量セットを、更に前記所定の予測モデルの初期値として前記濃色材の使用量を抑制しつつ前記近似性を極大化するように前記所定の予測モデルに基づいて予測された第2色材量セットであることを特徴とする印刷制御装置。
  2. 前記ターゲット値は、
    前記印刷装置にて前記記録媒体上に前記ターゲット値を再現するための前記色材量セットを前記所定の予測モデルに基づいて予測し、この予測した色材量セットを前記印刷装置に指定して確認パッチを印刷させ、この確認パッチの色を示す情報である確認ターゲット値と前記対象物の測色値であるターゲット測定値との偏差に基づいて設定された補正ターゲット値である請求項1に記載の印刷制御装置。
  3. 前記第2色材量セットを前記印刷装置に指定して再確認パッチを印刷させ、この再確認パッチの色を示す情報である再確認ターゲット値と前記ターゲット測定値との偏差に基づいて算出した再補正ターゲット値をターゲット値として前記第1色材量セットと前記第2色材量セットの再予測を行う請求項2に記載の印刷制御装置。
  4. 前記所定の予測モデルにおいては、色材量セットを予測する際に、前記印刷装置で付着させられる最小の単位量よりも少量ずつ色材量を変化させながら色材量セットの前記近似性を評価しており、
    前記所定の予測モデルに基づいて予測される色材量セットは、前記近似性が極大化した色材量セットについて前記単位量を丸め幅とした端数処理を行ったものである請求項1〜請求項3の何れか1項に記載の印刷制御装置。
  5. 予測された前記第2色材量セットを初期値として前記第1色材量セットと前記第2色材量セットを予測する過程を複数回繰り返し、
    該繰り返しの過程において、前記第2色材量セットにおける前記濃色材の使用量が2回連続して一致したことを検出した場合は、それ以降の繰り返しにおいて前記濃色材の使用量を固定する請求項1〜請求項4の何れか1項に記載の印刷制御装置。
  6. 前記第2色材量セットの予測において、使用量を抑制されたインクを除いた色材組合せによって全色相方向への色変更が可能である請求項1〜請求項5の何れか1項に記載の印刷制御装置。
  7. 前記複数の色材が、シアン(C)とマゼンタ(M)とイエロー(Y)とブラック(K)とライトシアン(lc)とライトマゼンタ(lm)とを含んで構成されており、
    前記第1色材量セットの予測においては、少なくともシアン(C)とマゼンタ(M)とブラック(K)の使用量を優先的に変化し、
    前記第2色材量セットの予測においては、少なくともライトシアン(lc)とライトマゼンタ(lm)とイエロー(Y)の使用量を優先的に変化する請求項1〜請求項6の何れか1項に記載の印刷制御装置。
  8. 少なくとも1つの色相に関して略同一の色相を有し濃度が互いに異なる濃色材と淡色材を含んで構成された複数の色材を印刷装置にて記録媒体に付着させて印刷を実行させるにあたり、前記色材の使用量の組合せである色材量セットを前記印刷装置に指定し、前記色材量セットに基づく印刷を実行させる印刷制御システムであって、
    前記印刷制御装置は、
    前記色材量セットとインデックスとの対応関係を規定したルックアップテーブルを参照することにより、指定されたインデックスに対応する前記色材量セットを前記印刷装置に指定して印刷させる印刷手段を備え、
    前記ルックアップテーブルにおいては、
    対象物の色を示す情報であるターゲット値を特定する前記インデックスと、前記印刷装置にて前記記録媒体上に付着したときに前記ターゲット値に対する近似性が極大化する色材量セットであるターゲット色材量セットと、の対応関係が規定されており、
    前記ターゲット色材量セットは、前記淡色材の使用量を抑制しつつ前記近似性を極大化するように所定の予測モデルに基づいて予測された第1色材量セットを、更に前記所定の予測モデルの初期値として前記濃色材の使用量を抑制しつつ前記近似性を極大化するように前記所定の予測モデルに基づいて予測された第2色材量セットであり、
    前記印刷装置は、
    前記色材量セットに基づいて印刷を実行する印刷実行手段を具備することを特徴とする印刷システム。
  9. 少なくとも1つの色相に関して略同一の色相を有し濃度が互いに異なる濃色材と淡色材を含んで構成された複数の色材を印刷装置にて記録媒体に付着させて印刷を実行させるにあたり、前記色材の使用量の組合せである色材量セットを前記印刷装置に指定し、前記色材量セットに基づく印刷を実行させる機能をコンピュータに実現させるための印刷制御プログラムであって、
    前記色材量セットとインデックスとの対応関係を規定したルックアップテーブルを参照することにより、指定されたインデックスに対応する前記色材量セットを前記印刷装置に指定して印刷させる印刷機能を備え、
    前記ルックアップテーブルにおいては、
    対象物の色を示す情報であるターゲット値を特定する前記インデックスと、前記印刷装置にて前記記録媒体上に付着したときに前記ターゲット値に対する近似性が極大化する色材量セットであるターゲット色材量セットと、の対応関係が規定されており、
    前記ターゲット色材量セットは、前記淡色材の使用量を抑制しつつ前記近似性を極大化するように所定の予測モデルに基づいて予測された第1色材量セットを、更に前記所定の予測モデルの初期値として前記濃色材の使用量を抑制しつつ前記近似性を極大化するように前記所定の予測モデルに基づいて予測された第2色材量セットであることを特徴とする印刷制御プログラム。
JP2008309030A 2008-12-03 2008-12-03 印刷制御装置、印刷システムおよび印刷制御プログラム Expired - Fee Related JP5157856B2 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2008309030A JP5157856B2 (ja) 2008-12-03 2008-12-03 印刷制御装置、印刷システムおよび印刷制御プログラム
US12/567,378 US20100134550A1 (en) 2008-12-03 2009-09-25 Printing control apparatus, printing system, and printing control program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2008309030A JP5157856B2 (ja) 2008-12-03 2008-12-03 印刷制御装置、印刷システムおよび印刷制御プログラム

Publications (2)

Publication Number Publication Date
JP2010136022A true JP2010136022A (ja) 2010-06-17
JP5157856B2 JP5157856B2 (ja) 2013-03-06

Family

ID=42222442

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2008309030A Expired - Fee Related JP5157856B2 (ja) 2008-12-03 2008-12-03 印刷制御装置、印刷システムおよび印刷制御プログラム

Country Status (2)

Country Link
US (1) US20100134550A1 (ja)
JP (1) JP5157856B2 (ja)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012186711A (ja) * 2011-03-07 2012-09-27 Fuji Xerox Co Ltd 色処理装置及びプログラム
US8711432B2 (en) 2011-06-14 2014-04-29 Seiko Epson Corporation Image processing device, printing apparatus, image processing method, and method of producing printing apparatus
JP2015091116A (ja) * 2013-11-07 2015-05-11 凸版印刷株式会社 色予測装置、色予測方法およびプログラム

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012138347A1 (en) * 2011-04-08 2012-10-11 Hewlett-Packard Development Company, L.P. Computing a spectrum of a sample
JP6108933B2 (ja) 2012-07-03 2017-04-05 キヤノン株式会社 キャリブレーション装置及びその制御方法
US9191552B2 (en) * 2012-11-30 2015-11-17 Kyocera Document Solutions Inc. Image processing apparatus that performs color conversion and image processing method
US9626607B2 (en) 2013-04-22 2017-04-18 Hewlett-Packard Development Company, L.P. Spectral print mapping
DE102014116201A1 (de) * 2014-11-06 2016-05-12 Krones Ag Vorrichtung und Verfahren zur Kontrolle von Direktdruckmaschinen
DE102017209860B4 (de) * 2017-06-12 2019-07-11 Henkel Ag & Co. Kgaa Verfahren und Vorrichtung zum Ermitteln eines Körperbereichzustands
JP6930238B2 (ja) * 2017-06-15 2021-09-01 富士フイルムビジネスイノベーション株式会社 色変換装置、画像形成装置およびプログラム
US10516809B2 (en) 2018-02-28 2019-12-24 Ricoh Company, Ltd. Optimizing number of grid points in multi-dimensional color lookup table
WO2019172919A1 (en) 2018-03-08 2019-09-12 Hewlett-Packard Development Company, L.P. Printer calibration
JP6996532B2 (ja) * 2019-07-19 2022-01-17 セイコーエプソン株式会社 情報処理装置、色変換プロファイル作成方法、色変換プロファイル作成プログラムおよび学習装置
WO2021061141A1 (en) * 2019-09-26 2021-04-01 Hewlett-Packard Development Company, L.P. Print coverage vectors
CN117549667B (zh) * 2024-01-11 2024-05-17 福建省海泓彩印有限公司 一种印刷生产线质量追踪管理系统及其使用方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH03194593A (ja) * 1989-12-25 1991-08-26 Fujitsu Ltd 指標表示と直値表示の同時表示制御方式
JP2001171182A (ja) * 1999-12-22 2001-06-26 Seiko Epson Corp 印刷制御装置、印刷装置、印刷方法、データ変換方法、および記録媒体
JP2007329899A (ja) * 2006-05-12 2007-12-20 Canon Inc 色処理装置及びその方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4192383B2 (ja) * 2000-02-24 2008-12-10 ソニー株式会社 デジタル放送受信システム、デジタル放送受信装置、並びにデジタル放送印刷装置
JP4040625B2 (ja) * 2002-09-12 2008-01-30 松下電器産業株式会社 画像処理装置、プリンタ装置、撮影装置、及びテレビ受像装置
JP4582310B2 (ja) * 2005-02-02 2010-11-17 セイコーエプソン株式会社 プロファイル作成方法、プロファイル作成装置、プロファイル作成プログラム、及び印刷制御装置
JP4595734B2 (ja) * 2005-08-03 2010-12-08 セイコーエプソン株式会社 プロファイル作成方法、プロファイル作成装置、プロファイル作成プログラム、印刷制御方法、印刷制御装置、及び印刷制御プログラム
JP4431897B2 (ja) * 2006-03-28 2010-03-17 セイコーエプソン株式会社 画像出力制御装置、及び画像出力制御方法
EP1973334A3 (en) * 2007-03-20 2010-01-27 Seiko Epson Corporation Printing control according to combinations of color materials
US8240795B2 (en) * 2007-07-18 2012-08-14 Canon Kabushiki Kaisha Printing method and printing apparatus

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH03194593A (ja) * 1989-12-25 1991-08-26 Fujitsu Ltd 指標表示と直値表示の同時表示制御方式
JP2001171182A (ja) * 1999-12-22 2001-06-26 Seiko Epson Corp 印刷制御装置、印刷装置、印刷方法、データ変換方法、および記録媒体
JP2007329899A (ja) * 2006-05-12 2007-12-20 Canon Inc 色処理装置及びその方法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012186711A (ja) * 2011-03-07 2012-09-27 Fuji Xerox Co Ltd 色処理装置及びプログラム
US8711432B2 (en) 2011-06-14 2014-04-29 Seiko Epson Corporation Image processing device, printing apparatus, image processing method, and method of producing printing apparatus
JP2015091116A (ja) * 2013-11-07 2015-05-11 凸版印刷株式会社 色予測装置、色予測方法およびプログラム

Also Published As

Publication number Publication date
JP5157856B2 (ja) 2013-03-06
US20100134550A1 (en) 2010-06-03

Similar Documents

Publication Publication Date Title
JP5157856B2 (ja) 印刷制御装置、印刷システムおよび印刷制御プログラム
US8199367B2 (en) Printing control device, printing system and printing control program
JP4561483B2 (ja) 指定色を再現するインク量の組み合わせの決定
JP4590424B2 (ja) 色処理装置及びその方法
JP2009230739A (ja) 印刷制御装置、印刷システムおよび印刷制御プログラム
US20090213392A1 (en) Printing control device, printing system and printing control program
JP2009169941A (ja) 印刷制御装置、印刷データ生成装置、印刷システムおよび印刷制御プログラム
JP2009171556A (ja) 印刷制御装置、印刷システムおよび印刷制御プログラム
JP2005178180A (ja) 色画像データ修正装置、色画像データ修正方法および色修正テーブル作成プログラム。
JP5182044B2 (ja) 印刷制御装置、印刷システムおよび印刷制御プログラム
JP2005178180A5 (ja)
US8243329B2 (en) Printing control device, print data generation device, printing system and printing control program
JP4946908B2 (ja) 印刷制御装置、印刷システムおよび印刷制御プログラム
JP5589520B2 (ja) 色変換プロファイル作成装置、色変換プロファイル作成方法、色変換プロファイル作成プログラムおよび印刷装置
JP2009177789A (ja) 印刷制御装置、印刷システムおよび印刷制御プログラム
JP5157873B2 (ja) 印刷制御装置、印刷システムおよび印刷制御プログラム
JP2010147783A (ja) 印刷結果予測方法および印刷制御装置
JP4924414B2 (ja) 印刷制御装置、印刷システムおよび印刷制御プログラム
JP5157872B2 (ja) 印刷制御装置、印刷システムおよび印刷制御プログラム
JP4985475B2 (ja) 印刷制御装置、印刷システムおよび印刷制御プログラム
JP2008072602A (ja) 印刷装置に対する色調整
US8488191B2 (en) Color mapping determination for an N-color marking device based upon image spatial noise defects
JP2008230050A (ja) 色材の組み合わせに応じた印刷制御
JP2005161550A (ja) 色変換テーブル補正装置、色変換テーブル補正方法、色変換テーブル補正プログラムおよび印刷制御装置
JP5050932B2 (ja) 印刷結果予測方法および印刷制御装置

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20111116

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20120824

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20120828

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20121016

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20121113

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20121126

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20151221

Year of fee payment: 3

S531 Written request for registration of change of domicile

Free format text: JAPANESE INTERMEDIATE CODE: R313531

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

LAPS Cancellation because of no payment of annual fees