(57)【要約】【課題】観察条件の不安定さや写像の非線形性という条件があったとしても、見本の色に合わせて物品に着色を施すことを簡単に行うことができる着色材料の配合量の算出方法を提供すること。【解決手段】見本物体を測色した複数の第１の分光数値データと、見本物体の色を目標として着色した際の複数の着色材料の配合量のデータである第１の配合数値データとのセットに基づいて分光数値データと配合数値データの近似関数を最適化計算によって作成し、新たな見本物体と同等といえる色の目的物体を着色する際に、この近似関数に新たな見本物体を測色した第１の分光数値データを適用して目的物体を着色するための着色材料の配合量のデータを正確かつ簡単に算出する。【選択図】図４
2018191073.tif 000002
【特許請求の範囲】
【請求項１】
所定の照明環境で人が目視によって観察したときに見本物体と同じ色に見えるように複数の着色材料を配合して目的物体を色付けする際の前記着色材料の配合量を算出するためのコンピュータで実行される算出方法であって、
ある１つの前記見本物体を測色することで得られる第１の分光数値データと、前記目的物体を前記ある１つの前記見本物体と同じ色に見えるように前記複数の着色材料を配合して色付けした際の前記複数の着色材料の配合量のデータ（以下、第１の配合数値データ）とを含む数値データを１セットとして、複数の前記見本物体と対応する複数の前記目的物体について第１の分光数値データと第１の配合数値データの複数のセットを取得し、得られた第１の分光数値データと第１の配合数値データの複数のセットに基づいて分光数値データと配合数値データの近似関数を最適化計算によって作成し、得られた前記近似関数に新たな前記見本物体を測色することで得られる第１の分光数値データを適用することで新たな前記見本物体と同じ色に見える新たな前記目的物体のための前記着色材料の配合量のデータを算出するようにしたことを特徴とする着色材料の配合量の算出方法。
【請求項２】
前記目的物体に色付けをし、その色付けをした前記目的物体を測色することで得られる第２の分光数値データと、前記目的物体を色付けした際の第２の配合数値データとを含む数値データを１セットとして、このセットを更に前記近似関数を最適化する計算に用いるようにしたことを特徴とする請求項１に記載の着色材料の配合量の算出方法。
【請求項３】
第１の色付け処理で、ある１つの前記見本物体と同じ色に見えるように色付けしたにも関わらず前記目的物体の色が前記ある１つの前記見本物体と同じ色には見えないとき、第２の色付け処理として、前記ある１つの前記見本物体と同じ色に見えるように追加的に色付けしその追加的な色付けをした前記目的物体を測定することで第３の分光数値データを得、
前記近似関数に前記第３の分光数値データを適用することで追加的に色付けした前記見本物体と同じ色に見える新たな前記目的物体のための前記着色材料の配合量のデータ（以下、第３の配合数値データ）量を算出し、
前記ある１つの前記見本物体を測色することで得られる前記第１の分光数値データと前記第３の配合数値データとを含む数値データをセットとして、このセットを更に前記近似関数を最適化する計算に用いるようにしたことを特徴とする請求項１又は２に記載の着色材料の配合量の算出方法。
【請求項４】
前記近似関数に前記ある１つの前記見本物体を測色することで得られる前記第１の分光数値データを適用することで新たな前記見本物体と同じ色に見える新たな前記目的物体のための前記着色材料の配合量のデータ（以下、第４の配合数値データ）量を算出し、
前記第３の分光数値データと前記第４の配合数値データをセットとして、前記近似関数を最適化する計算に用いるようにしたことを特徴とする請求項３に記載の着色材料の配合量の算出方法。
【請求項５】
第１の色付け処理で、ある１つの前記見本物体と同じ色に見えるように色付けしたにも関わらず前記目的物体の色が前記ある１つの前記見本物体と同じ色には見えないとき、第２の色付け処理として、前記ある１つの前記見本物体と同じ色に見えるように追加的に色付けしその追加的な色付けをした前記目的物体を測色することで第３の分光数値データを得、
前記近似関数に前記ある１つの前記見本物体を測色することで得られる前記第１の分光数値データを適用することで新たな前記見本物体と同じ色に見える新たな前記目的物体のための前記着色材料の配合量のデータ（以下、第４の配合数値データ）量を算出し、
前記第３の分光数値データと前記第４の配合数値データをセットとして、前記近似関数を最適化する計算に用いるようにしたことを特徴とする請求項１又は２に記載の着色材料の配合量の算出方法。
【請求項６】
前記第３の分光数値データは分光透過率データ又は分光反射率データであることを特徴とする請求項３〜５のいずれかに記載の着色材料の配合量の算出方法。
【請求項７】
前記第１又は第２の分光数値データは分光透過率データ又は分光反射率データであることを特徴とする請求項１〜６のいずれかに記載の着色材料の配合量の算出方法。
【請求項８】
数値データとして分光数値データ以外の前記見本物体の数値データを用いることを特徴とする請求項１〜７のいずれかに記載の着色材料の配合量の算出方法。
【請求項９】
前記近似関数は、機械学習により得られる非線形的な計算手段を含むことを特徴とする請求項１〜８のいずれかに記載の着色材料の配合量の算出方法。
【請求項１０】
前記近似関数としてニューラルネットワークを用いることを特徴とする請求項１〜９のいずれかに記載の着色材料の配合量の算出方法。
【請求項１１】
前記見本物体及び前記目的物体は透明体であり、染料を染みこませて色付けすることを特徴とする請求項１〜１０のいずれかに記載の着色材料の配合量の算出方法。
【請求項１２】
前記見本物体及び前記目的物体は眼鏡用レンズであることを特徴とする請求項１１に記載の着色材料の配合量の算出方法。
【請求項１３】
前記見本物体及び前記目的物体は非透明体であり、表面に色素を付着させて色付けすることを特徴とする請求項１〜１０のいずれかに記載の着色材料の配合量の算出方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、複数の着色材料を使用して見本の色に合わせて物品に着色を施して人が見本と同じ色であると認識できる再現色を調整するための着色材料の配合量の算出方法に関するものである。
【背景技術】
【０００２】
例えば、多くのカラーバリエーションを持つ物品（工業製品）では、その物品の製造コストが高ければ在庫を持つことは効率が悪いため、注文に応じてその都度彩色を施すことが合理的である。彩色を施す場合には注文に応じた見本を参照し、複数の異なる着色材料を混ぜ合わせたりあるいは重ね合わせることで物品に彩色を施してその見本と同じ色となるように（つまり同一性のあるように）再現色を調整することとなる。
しかし、実際には見本と同一性のある再現色を調整することは難しい。その理由としては、例えば、見本の素材と実際に製造する物品との素材の違いがあること、新しいカラーの種類ができたり、染料の種類が変わったりすること、見本がどのような染料や素材を用いているのかまったく不明であること、等が挙げられる。
そのため、見本の色に合わせて物品に着色を施すことを簡単に行う技術が求められていた。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１６−２２５９４０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
そのような技術の一例として特許文献１を挙げる。特許文献１は色推定システムや色推定方法に関する技術である。特許文献１では、濃度階調出現率テーブルをもとに原色インキＣＭＹの量を決定している。網点を構成する各濃度階調の濃度階調領域の面積の比率）から分光反射率Ｒ'（ｓ，λ）を算出する。そして、算出された算出分光反射率Ｒ'（ｓ，λ）の各々と、測定分光反射率Ｒｓ（λ）との平均二乗誤差ＲＭＳＥを、指令網点面積率毎に所定の波長範囲において求める。そして、濃度階調出現率算出部は、算出分光反射率Ｒ'（ｓ，λ）と測定分光反射率Ｒｓ（λ）との平均二乗誤差（ＲＭＳＥ）が最も小さくなる濃度階調領域の各々の出現率を求めるものである。
しかしながら、分光反射率パターンを似せても、観察条件の不安定さがあるため、必ずしも同じ色に見えるとは限らないものである。「観察条件の不安定さ」とは、例えば人間が対象物に相対して観察するため、眼と対象物と照明の位置関係が毎回変動すること。見本物体と目的物体の配置（どちらを右側に置いて見るか）。光源となるランプの劣化に伴う分光強度の経時的な変化。色の見え方（感じ方）が人によって微妙に異なること等である。
【０００５】
また、観察条件の不安定だけではなく、そもそも色を写像する場合には非線形性が大きいことも挙げられる。
例えば、眼鏡レンズを例に取ると、眼鏡レンズにおける規格化された見本カラーレンズは、例えば数十種類から、百種類以上もある。それらのカラーを適当な色空間、例えばＣＩＥ１９７６ (Ｌ＊、ａ＊、ｂ＊) 色空間に配置することで、各カラーのＬ＊ａ＊ｂ＊座標からＲＹＢ設定量を得るための近似関数を作ることができる。例えば、見本カラー注文の見本レンズの分光透過率を測定し、その結果をもとにＬ＊ａ＊ｂ＊座標を求め、さらに近似関数によってＲＹＢ設定量を算出する方法が考えられる。
しかし、実際にその方法では、見本レンズと加工レンズの見た目の色を一致させることはできない。色空間の領域によってＬ＊ａ＊ｂ＊座標からＲＹＢ設定量を写像する近似関数が大きく異なるためである（非線形性が大きいためである）。この場合に、十分に有用な近似関数を得るためには、極めて膨大な数のデータが必要になる。さらに上記した観察条件の不安定性があるので、バラツキを含むデータをもとにするので、非常に多くのデータが必要となる。
これらのような課題から、観察条件の不安定さや写像の非線形性という条件があったとしても、見本の色に合わせて物品に着色を施すことを簡単に行うことができる技術が求められていた。
【課題を解決するための手段】
【０００６】
上記課題を解決するために手段１では、所定の照明環境で人が目視によって観察したときに見本物体と同じ色に見えるように複数の着色材料を配合して目的物体を色付けする際の前記着色材料の配合量を算出するためのコンピュータで実行される算出方法であって、ある１つの前記見本物体を測色することで得られる第１の分光数値データと、前記目的物体を前記ある１つの前記見本物体と同じ色に見えるように前記複数の着色材料を配合して色付けした際の前記複数の着色材料の配合量のデータ（以下、第１の配合数値データ）とを含む数値データを１セットとして、複数の前記見本物体と対応する複数の前記目的物体について第１の分光数値データと第１の配合数値データの複数のセットを取得し、得られた第１の分光数値データと第１の配合数値データの複数のセットに基づいて分光数値データと配合数値データの近似関数を最適化計算によって作成し、得られた前記近似関数に新たな前記見本物体を測色することで得られる第１の分光数値データを適用することで新たな前記見本物体と同じ色に見える新たな前記目的物体のための前記着色材料の配合量のデータを算出するようにした。
【０００７】
このような構成においては、見本物体を測色した複数の第１の分光数値データと見本物体の色を目標として着色した際の第１の配合数値データのセットに基づいて分光数値データと配合数値データの近似関数を最適化する計算によって作成することによって、新たな見本物体と同等といえる色の目的物体を着色する際に、この近似関数に新たな見本物体を測色した第１の分光数値データを適用することで、そのための着色材料の配合量のデータを正確かつ簡単に算出することができる。
このような手法であれば、物体を目視する際の照明光源が多種類であっても、あるいは反射率や透過率が角度によって変わる影響で色の見え方が変わったとしても、着色材料の配合量のデータは妥当なものとなる。
手段１における「第１の分光数値データと第１の配合数値データの複数のセットに基づいて分光数値データと配合数値データの近似関数を最適化計算によって作成」することをわかりやすく模式化すると図４における左上Ａからのデータセットによる最適化に相当する。
ここに、「近似関数」とは、本発明では着色材料の配合量を算出するために必要な精度を持った関数であって、新たなデータ又はデータセットを適用して最適化する計算をすることで精度を上げることができる。近似関数は例えば機械学習、特に深層学習により得られる非線形的な計算手段を含んでいる。
「第１の分光数値データと第１の配合数値データ」のセットは、複数セットが必要である。この発明は近似関数を最適化するための計算を行うために複数の見本物体の色を前提としているためである。見本物体は多いほど計算上精度が増す。
実際に目的物体をある１つの見本物体と同じ色に見えるように着色する場合においては、過去に同じ見本を着色した経緯があれば過去の配合量のデータを使用することができるが、まったく過去にデータがない場合や、材質も使用した着色材料の種類も不明な場合等では、目的物体を見本物体と同じ色に見えるように着色することはなかなか困難である。最適化計算の初期段階では精度を上げるために、そのような場合には見本物体と同じ色に見えるものを得るまで（つまり、第１の配合数値データを取得できるまで）何度も着色を繰り返すことも必要となる。見本物体と同じ色に見えない目的物体を得た場合には、手段２〜５において最適化の計算に使用することが可能である。
ここに「見本物体」とは製品の発注者と受注者が同一物（またはその複製）として参照するサンプルをいう。「見本物体の色」が必ずしも実際の製品に着色されているとは限らない。単に色見本である場合もある。見本物体の色は毎回異なり、過去に使用した見本物体であってもよい。
また「目的物体」は、基本的に「見本物体」と同じ色になるように着色したものであり、透明であっても非透明であってもよい。また、手段２のように「見本物体」と同じ色にならない場合であっても近似関数を最適化する計算に使用することができる。
また「複数の着色材料を配合して目的物体を色付け」する際の前提条件として複数の着色材料のそれぞれの配合量がデータとして把握できる必要がある。つまり、複数の着色材料は数値化できなければならない。複数の着色材料は複数であればいくつであってもよい。見本物品と同じ着色材料でもよく、異なる着色材料であってもよい。
発明に使用される「コンピュータ」はＣＰＵ（中央処理装置）及びＲＯＭやＲＡＭ等の周辺装置によって構成される。ＣＰＵは各種プログラムやオペレータの操作に従って入力装置（例えば、キーボードや、マウス）からのデータ入力に基づいて計算を行う。
【０００８】
また、手段２では、前記目的物体に色付けをし、その色付けをした前記目的物体を測色することで得られる第２の分光数値データと、前記目的物体を色付けした際の第２の配合数値データとを含む数値データを１セットとして、このセットを更に前記近似関数を最適化する計算に用いるようにした。
手段１では「見本物体を測色した複数の第１の分光数値データと見本物体の色を目標として着色した際の第１の配合数値データのセット」を複数セット取得することで分光数値データと配合数値データの近似関数を最適化計算によって算出するようにしていたが、手段２では、「色付けした目的物体を測色した複数の第２の分光数値データとその色付けした目的物体を着色した際の第２の配合数値データのセット」を手段１の近似関数の最適化計算のために手段１のセットに追加して加えるというものである。このようなセットのデータを加えることで、データ数が増えて近似関数を最適化する計算に貢献し、近似関数の精度が増すこととなる。「色付けした目的物体」は「見本物体」と同じ色であってもよく、違う色であってもよい。
手段２における最適化計算をわかりやすく模式化すると図４における左下Ｂからのデータセットによる最適化に相当する。
【０００９】
また、手段３では、第１の色付け処理で、ある１つの前記見本物体と同じ色に見えるように色付けしたにも関わらず前記目的物体の色が前記ある１つの前記見本物体と同じ色には見えないとき、第２の色付け処理として、前記ある１つの前記見本物体と同じ色に見えるように追加的に色付けしその追加的な色付けをした前記目的物体を測定することで第３の分光数値データを得、前記近似関数に前記第３の分光数値データを適用することで追加的に色付けした前記見本物体と同じ色に見える新たな前記目的物体のための前記着色材料の配合量のデータ（以下、第３の配合数値データ）量を算出し、前記ある１つの前記見本物体を測色することで得られる前記第１の分光数値データと前記第３の配合数値データとを含む数値データをセットとして、このセットを更に前記近似関数を最適化する計算に用いるようにした。
第１の色付け処理と第２の色付け処理では着色材料の配合量が別個独立であって、このように２つの処理を行うことで、追加的に色付けした目的物体を測色して第３の分光数値データを得ることはできるが、第１の色付け処理と第２の色付け処理の処理方法が異なる場合では色付けに使用したトータルの配合数値データを得ることはできない。
しかし、このように追加的に着色した目的物体を測色して得た第３の分光数値データを近似関数に入力した場合に出力される配合量と、見本物体を測色して得た第１の分光数値データを近似関数に入力した場合に出力される配合量は近似しているものである。そのため、この段階で第３の分光数値データを現段階で最適化して得ている近似関数に適用することで第３の配合数値データを算出し、第１の分光数値データとこの第３の配合数値データをセットとして、更に前記近似関数を最適化する計算に用いるようにして近似関数の精度を上げることができる。
手段３における最適化計算をわかりやすく模式化すると図４における右下Ｃからのデータセットによる最適化に相当する。
手段３は、例えば見本物体と同じ色に見えるように着色してみたところ、見本物体と同じ色に見えない目的物体を得てしまった場合に活用できる手法である。
ここに、第１の色付け処理とは例えば後述する手段１１ではドライ染色が相当し、ウェット染色が第２の色付け処理に相当する。
【００１０】
また、手段４では、前記近似関数に前記ある１つの前記見本物体を測色することで得られる前記第１の分光数値データを適用することで新たな前記見本物体と同じ色に見える新たな前記目的物体のための前記着色材料の配合量のデータ（以下、第４の配合数値データ）量を算出し、前記第３の分光数値データと前記第４の配合数値データをセットとして、前記近似関数を最適化する計算に用いるようにした。
また、手段５では、第１の色付け処理で、ある１つの前記見本物体と同じ色に見えるように色付けしたにも関わらず前記目的物体の色が前記ある１つの前記見本物体と同じ色には見えないとき、第２の色付け処理として、前記ある１つの前記見本物体と同じ色に見えるように追加的に色付けしその追加的な色付けをした前記目的物体を測色することで第３の分光数値データを得、前記近似関数に前記ある１つの前記見本物体を測色することで得られる前記第１の分光数値データを適用することで新たな前記見本物体と同じ色に見える新たな前記目的物体のための前記着色材料の配合量のデータ（以下、第４の配合数値データ）量を算出し、前記第３の分光数値データと前記第４の配合数値データをセットとして、前記近似関数を最適化する計算に用いるようにした。
追加的に着色した目的物体を測色して得た第３の分光数値データを近似関数に入力した場合に出力される配合量と、見本物体を測色して得た第１の分光数値データを近似関数に入力した場合に出力される配合量は近似しているものである。そのため手段４及び手段５において、第１の分光数値データを現段階で最適化して得ている近似関数に適用することで、第４の配合数値データを算出し、第３の分光数値データとこの第４の配合数値データをセットとして、更に前記近似関数を最適化する計算に用いるようにして近似関数の精度を上げることができる。
手段４又は５における最適化計算をわかりやすく模式化すると図４における右上Ｄからのデータセットによる最適化に相当する。
【００１１】
また、手段６では、前記第３の分光数値データは分光透過率データ又は分光反射率データであるようにした。
また、手段７では、前記第１又は第２の分光数値データは分光透過率データ又は分光反射率データであるようにした。
第１〜第３の分光数値データとして分光透過率データ又は分光反射率データを使用することは正確な分光数値データの取得手段として望ましい。
また、手段８では、数値データとして分光数値データ以外の前記見本物体の数値データを用いるようにした。
上記第１〜第３の分光数値データ以外の前記見本物体の数値データを用いることで近似関数の精度を上げることができるからである。例えば、基材の種類などのカテゴリデータやレンズであればレンズの度数等の数値データである。
また、手段９では、前記近似関数は、機械学習により得られる非線形的な計算手段を含むようにした。
機械学習の推定アルゴリズムは、線形回帰、ボルツマンマシン、ニューラルネットワーク、サポートベクターマシン、ベイジアンネットワークを用いた統計的推定、強化学習、深層学習の内の１つまたは複数を含む。
機械学習とは、学習用データセットを使ってコンピュータプログラムを訓練した後に、未知の例について正確に判断できるようにすることである。本発明においてはコンピュータに分光数値データ（加えて、基材の種類などのカテゴリデータやレンズの度数等の）その他の数値データ）を近似関数に入力したときに出力される値が配合数値データに近づくよう機械学習をさせるようにすることがよい。
機械学習により得られる非線形写像を用いて近似関数を更新していくことが可能となる。機械学習では第１の分光数値データと第１の配合数値データのデータセット又は／及び第２の分光数値データと第２の配合数値データのデータセットによる教師あり学習によって近似関数を収束させることができ、更に手段３〜５のように第１の分光数値データと第３の配合数値データのデータセット又は／及び第３の分光数値データと前記第４の配合数値データのデータセット（これらは完全な教師あり学習ではなく、一部のみ教師ありの半教師あり学習となる）を組み合わせることによって、さらに好ましい推定値に収束させることができる。
また、手段１０では、近似関数としてニューラルネットワークを用いるようにした。
ニューラルネットワークとは、脳機能に見られるいくつかの特性を計算機上のシミュレーションによって表現することを目指した数学モデルである。てニューラルネットワークは機械学習の一種である。ニューラルネットワークとしては３層以上のニューロンを重ねるディープニューラルネットワーク（以下、ＤＮＮ）がよりよい。ニューラルネットワークを最適化する計算には誤差逆伝播法を用いることがよい。以下の実施の形態ではＤＮＮを一例とした。
【００１２】
また、手段１１では、前記見本物体及び前記目的物体は透明体であり、染料を染みこませて色付けするようにした。
また、手段１２では、前記見本物体及び前記目的物体は眼鏡用レンズであるようにした。
特に眼鏡用レンズでは、次のような課題があるため、本発明を適用することが非常に適している。
眼鏡レンズを色付けして納品する「カラー注文」には、規格カラー注文と見本カラー注文がある。以後これらに対応する染色加工のことをそれぞれ規格染色・見本染色と呼ぶ。各種の規格カラーを表す見本レンズを眼鏡店とレンズメーカーが所持し、各カラーを表す名前または番号を用いて受発注を行い、製造時はメーカー側が、納品後は眼鏡店側がそれぞれの所持する見本レンズと色を比較する。これに対して見本染色用の見本レンズは、毎回異なるものが眼鏡店側から提示される。
規格染色のドライ染色条件を数値的に設定することは比較的容易である。ドライ染色後のレンズの色を所定の照明環境で見本と比較観察して、同じ色に見えるようにＲＹＢの各設定量を定めれば良い。これは設定量を変えながら何回も加工を繰り返すことで、カラーごとの条件を決定することができる。ただし、レンズ基材の種類は屈折率別などで複数種類を扱うので、それぞれの基材に対応する条件を定めなくてはならない。だから「見本レンズがどこかのメーカーで加工されたのだから、それを作った条件で加工すれば良い」ということではない。また、レンズの厚さや色の付きやすさは度数や厚さによって少しずつ異なるので、それらに対応した条件を決定することは難しい。そのため、標準的なレンズ（たとえば度数の無いレンズ）で定めた条件をすべての度数のレンズに適用するなどして、妥協している。見本染色に関してこのようなことを行うと、受注するごとに何枚もレンズを無駄にすることになる。
一方、ウェット染色は再現性に劣るので、加工条件を数値で確定することは難しい。また、染色液を一定温度で維持する設備のコストがかかる。加工枚数が少なければ、１枚あたりコストが高くなる。加工枚数が多ければ、同じカラーのレンズを多く加工するのであればドライ染色のほうが有利になる。しかし、規格染色において加工したレンズの色の仕上がり具合が悪い場合などは、さらにウェット染色を追加で行って色を整える。これを調色という。
このような状況から、眼鏡店から提示される毎回カラーが異なる見本レンズに関して、そのドライ染色条件を簡単に決定する方法が望まれていた。しかし、上記のように観察条件の不安定さがあることや、写像の非線形性が大きいこと等からドライ染色条件を簡単に決定することは困難であった。更に、レンズでは入力データとしてＬ＊ａ＊ｂ＊座標または分光透過率データだけでなく、レンズの度数や厚さを含める必要がある。
本発明を眼鏡用レンズに適用すれば、基材・度数・厚さ・ハーフカラー（ハーフカラーについては後述する）などの多様な条件に対応し、近似関数に求められる非線形性にも対応し、観察条件の不安定性に由来するデータのバラツキを統計的に処理することで、見本レンズを測色したデータに基づいて最適なドライ染色条件を算出することができる。
また、この発明を眼鏡用レンズに適用すれば、規格カラーの条件を決定することもできる。カラーの規格が一斉に変わるときは、新しい見本レンズのセットに切り替えるが、その新しい見本レンズと同じ色に見えるようなレンズを加工するための近似関数は、見本染色の加工を行ううちに高い精度レベルで最適化されるからである。
ここに、「ドライ染色」とは、所定の用紙に染料をインックジェットでプリントし、用紙をレンズに近接させた状態で保持して電気炉に入れて熱することで、染料をレンズに付着させて、基材中に侵入させる色付け方法である。「ウェット染色」とは溶媒に染料を溶かしこんだ液体に透明な目的物体を浸漬して染料を染みこませる方法である。
【００１３】
また、手段１３では、前記見本物体及び前記目的物体は非透明体であり、表面に色素を付着させて色付けするようにした。
これは、例えば、物体表面に印刷するようなケースである。物体表面が白いことがよいが、色素を付着させることができれば必ずしも白くなくともよい。
【発明の効果】
【００１４】
本発明では、観察条件の不安定さや写像の非線形性があっても、正確かつ簡単に見本の色に合わせて物品に着色を施すことができる。
【図面の簡単な説明】
【００１５】
【図１】横軸を波長、縦軸を透過率×１００％とした分光透過率のグラフ。
【図２】実施の形態におけるＤＮＮの学習手法における入力層から複数の隠れ層を経て出力されるまでのネットワーク状態を説明する模式図。
【図３】他の実施の形態におけるＤＮＮの学習手法における入力層から複数の隠れ層を経て出力されるまでのネットワーク状態を説明する模式図。
【図４】本発明における近似関数の最適化の概要を説明する模式図。
【図５】近似関数の計算に使用する入力値データの一例のグラフ。
【図６】近似関数の計算に使用する入力値データの一例のグラフ。
【図７】近似関数の計算に使用する入力値データの一例のグラフ。
【図８】近似関数の計算に使用する入力値データの一例のグラフ。
【発明を実施するための形態】
【００１６】
以下、図面に基づいて本発明の具体的な実施の形態の説明をする。
本実施の形態では、眼鏡店を通じてユーザーが選択した見本レンズをレンズメーカーに送付し、レンズメーカーにてこの見本レンズと同じ色になるようにドライ染色する場合のドライ染色用のＲ（赤）Ｙ（黄）Ｂ（青）の三色の配合量を算出するものとする。レンズ情報としてはその他、基材の種類や処方等が提供されるものとする。近似関数としてＤＮＮを用いる。
１．着色材料の配合量の算出について
ＤＮＮでは入力値を設定し、それに重みを乗じてノード（ニューロン）の値を算出し、最終段階でＲＹＢの配合量を出力値として得る。出力値が適切になるように近似関数を最適化する点については後述する。まず、ここで使用する入力値の種類について説明する。入力値の一般形を次の数１のように表示する。  この一般形では入力の順序を表すｋを下に添えて表す。本実施の形態では４種類４１個の入力値を設定するが、これらは適宜変更可能である。
【００１７】
【数１】
2018191073.tif 000003
【００１８】
１）見本レンズの分光透過率
例えば測定装置である分光透過率計を使用して分光透過率を測定することができる。図１は、横軸を波長、縦軸を透過率×１００％としたカラーレンズの分光透過率の一例である。この図１ではａ．見本レンズを測定した測色データ、ｂ．見本レンズと同じ色に見えるようにドライ染色した新たなレンズを測定した測色データ、ｃ．ドライ染色したが見本レンズと同じ色に見えなかったためウェット調色した新たなレンズを測定した測色データの３種類のレンズ特性が表示されたグラフである。本発明では、入力値としてはａ．の見本レンズの分光透過率を使用するわけである。そして、出力値としてｂ．のようなレンズを作製するためのＲＹＢの配合量を正確に計算で求めるというものである。本実施の形態では３９０〜７３０ｎｍの帯域において１０ｎｍ間隔で、３５個の波長の分光数値データを取得したａ．の見本レンズの３５個の波長の分光数値データを取得する。この間隔や取得する波長帯域や取得個数は一例であって変更可能である。
【００１９】
２）見本レンズの度数
眼鏡店を通じて預かった見本レンズの度数をレンズメーターにて測定する。
３）加工するレンズの度数
ユーザーの処方に基づいてレンズメーカーで加工するレンズの度数。この値は眼鏡店からの発注情報にて指定される。
４）基材の種類を表すカテゴリ値
基材の種類は眼鏡店からの発注情報にて指定される。加工するレンズの種類はコートのバリエーションを加えると数が多くなるが、本実施の形態では４種類で例示する。定数を４個設定し、１番目の基材を｛１，０，０，０｝、２番目の基材を｛０，１，０，０｝・・・として表わす。例えば下記式のように１番目の基材のデータは表される。
【００２０】
【数２】
2018191073.tif 000004
【００２１】
次に、入力値の基準化を図ることについて説明する。基準化をすることによって近似関数を最適化する計算を効率良く実行できることがわかっているためである。以下では、分光透過率の場合を例にとって説明する。基準化は各波長の透過率について、全データの平均値と標準偏差を求める。そして、各データから平均値を引いて標準偏差で割るという作業となる。数式で表すと下記数１〜数３となる。
【００２２】
【数３】
2018191073.tif 000005
【００２３】
【数４】
2018191073.tif 000006
【００２４】
【数５】
2018191073.tif 000007
【００２５】
ｙの添え字１は、波長３９０ｎｍ（１番目の波長）を基準化した値であることを表している。尚、ｙにはデータの順番を表すｐを添えない。１つのデータを基準化して入力値を作ったら、その入力値を近似関数に入力して出力を得るようにしているためである。見本レンズの度数と加工するレンズの度数を基準化する手順は同様なので省略する。また、本実施の形態では基材を表すカテゴリ値は、基準化しないで０と１のまま入力値ｙ３８〜ｙ４１とした。
【００２６】
次に、近似関数の計算プロセスについて説明する。以下の計算は複雑であるためコンピュータで実行される。伝達関数としてＲｅＬＵ（rectified linear unit整流された線形素子）を使用する。本実施の形態では隠れ層は４層で構成されている。以下の数６〜数１７の計算を模式化したネットワーク図が図２である。
まず、４１個の入力層の各ノードから第１層のノード２０個のうち、第１〜１９番目の１９個のノードに対して４１×１９＝７７９個の重みを設定する。入力ノードの値に重みの値を乗じて合計した値を第１層の各ノードの値とする。ただしマイナスの値は０にする。第１層の２０番目のノードの値は常に１とする。
各層の重みの一般形は下記数６で示す。重みの値ｗに対して上の層を表すｉを上に、下の層の順序を表すｋと上の層の順序を表すｊを下に添える。下記数７において関数ｆはＲｅＬＵを表す。
重みはこの段階までの近似関数（つまり、ここではＤＮＮ）の最適化の計算において更新された最新の重みが適用される。
【００２７】
【数６】
2018191073.tif 000008
【００２８】
【数７】
2018191073.tif 000009
【００２９】
次に第１層の出力の一般形を表す式を示す。ｊ＝１〜１９に関する式が数８で、ｊ＝２０に関する式が数９である。
第１層のｊ番目の出力値を表す式ではｉ＝１、入力層のノード数Ｎｉ−１＝４１とする。
【００３０】
【数８】
2018191073.tif 000010
【００３１】
【数９】
2018191073.tif 000011
【００３２】
次に、第１層のノード２０個から第２層のノード１５個のうち、第１〜１４番目の１４個のノードに対して２０×１４＝２８０個の重みを設定する。ここでも、入力ノードの値に重みの値を乗じて合計した値を第２層の各ノードの値とする。ただしマイナスの値は０にする。第２層の１５番目のノードの値は常に１とする。
ここで、数１０に示すように第２層より上の第ｉ層の出力の一般形は、入力ｙをｚにおきかえて表す。第ｉ−１層（下の層）のノード数をＮｉ−１で表わす。
【００３３】
【数１０】
2018191073.tif 000012
【００３４】
第２層のｊ番目の出力値を表す式では、ｉ＝２、Ｎｉ−１＝２０となる。ｊ＝１〜１４に関する式が数１１で、ｊ＝１５に関する式が数１２である。
【００３５】
【数１１】
2018191073.tif 000013
【００３６】
【数１２】
2018191073.tif 000014
【００３７】
次に、第２層のノード１５個から第３層のノード１０個のうち、第１〜９番目の９個のノードに対して１５×９＝１３５個の重みを設定する。ここでも、入力ノードの値に重みの値を乗じて合計した値を第３層の各ノードの値とする。ただしマイナスの値は０にする。第３層の１０番目のノードの値は常に１とする。第３層のｊ番目の出力値を表す式では、ｉ＝３、Ｎｉ−１＝１５となる。ｊ＝１〜９に関する式が数１３で、ｊ＝１０に関する式が数１４である。
【００３８】
【数１３】
2018191073.tif 000015
【００３９】
【数１４】
2018191073.tif 000016
【００４０】
次に、第３層のノード１０個から第４層のノード７個のうち、第１〜６番目の６個のノードに対して１０×６＝６０個の重みを設定する。ここでも、入力ノードの値に重みの値を乗じて合計した値を第４層の各ノードの値とする。ただしマイナスの値は０にする。第４層の７番目のノードの値は常に１とする。第４層のｊ番目の出力値を表す式では、ｉ＝４、Ｎｉ−１＝１０とする。ｊ＝１〜６に関する式が数１５で、ｊ＝７に関する式が数１６である。
【００４１】
【数１５】
2018191073.tif 000017
【００４２】
【数１６】
2018191073.tif 000018
【００４３】
最後に、第４層のノード７個から第５層（出力層）の３個のノードに対して７×３＝２１個の重みを設定する。ここでも、入力ノードの値に重みの値を乗じて合計した値を第５層の各ノードの値とし、そのまま出力値とする。出力層においては、マイナスの値もそのまま出力する。したがって、ＲｅＬＵを表す関数ｆは用いない。第５層のｊ番目の出力値を表す式では、ｉ＝５、Ｎｉ−１＝７とする。ｊ＝１〜３に関する式が数１７である。
３個の出力値の指数から１を引けば、数１８のようにＲＹＢの設定量（配合量）を得る。出力値の指数をとって１を引く理由は段落００４７で後述する。
この計算では重みの数は７７９（第１層）＋２８０（第２層）＋１３５（第３層）＋６０（第４層）＋２１（第５層）＝１２７５個であった。
【００４４】
【数１７】
2018191073.tif 000019
【００４５】
【数１８】
2018191073.tif 000020
【００４６】
２．近似関数の更新について
次に、上記のようなＤＮＮの計算を行う際に用いる重みを更新して計算精度を向上させる手法について説明する。出力された値がＲＹＢの参照値に近似するように、すべての重みの値を最適化するようにした。本実施の形態ではバックプロパゲーション（誤差逆伝播法）を用いて最適化を行った。バックプロパゲーションでは出力値と参照値との差の二乗和が最小になるような最適化計算を行う。下記数１９で最小化すべき残差２乗和の一般式を示す。ｎはデータの総数。係数１／２は計算の都合で付してある。ｏｐｊはｐ番目のＲＹＢ配合量の出力値ｊ（但し、ｊ＝１〜３）であり、下記数２０のように上記「１．着色材料の配合量の算出について」の出力と等価である。ｔｐｊは同じくｐ番目のＲＹＢ配合量の参照値ｊである。
【００４７】
【数１９】
2018191073.tif 000021
【００４８】
【数２０】
2018191073.tif 000022
【００４９】
出力値とはここでは上記「１．着色材料の配合量の算出について」で直近のある見本レンズＡの分光数値データ及びその他の数値データを入力値としてＤＮＮによって算出したＲＹＢの配合量である。
一方、参照値は近似関数の出力すべき目標値であって、以下の１）〜４）のいずれかに属するものである。
１）見本レンズＡと同じ色に見えるようにドライ染色したレンズに使用したＲＹＢの配合量
２）任意のある色にドライ染色した際のそのレンズのドライ染色に使用したＲＹＢの配合量
３）見本レンズＡを目標にドライ染色したが同じ色に見えなかったのでウェット調色して見本レンズＡと同じ色に見えるように調色したレンズについて分光数値データ及びその他の数値データを入力値として上記「１．着色材料の配合量の算出について」でＤＮＮによって算出したＲＹＢの配合量
４）見本レンズＡの分光数値データ及びその他の数値データを入力値として上記「１．着色材料の配合量の算出について」でＤＮＮによって算出したＲＹＢの配合量
図５〜図８に、本実施の形態で扱う具体的な入力値データの一例を示す。図５〜図７はＤＮＮに必要な項目のみを表示するようにして、計算に無関係なデータは割愛した。また、このようなデータは近似関数の更新については多ければ多いほど精度が上がる。図５は上記の参照値１）のケースでの入力値データの例である。また、図６は上記の参照値２）のケースでの入力値データの例である。これらは実際に測色しているためＲＹＢの配合量をデータとして使用している。図７は上記の参照値３）の図８は４）のケースでの入力値データの例である。図７及び図８ではＲＹＢの配合量は計算で得られた出力値を使用しているため、データ情報としては使用していない。
【００５０】
尚、参照値は下記式２１のようにＲＹＢ設定量（配合量）に１を加えた値の対数で表される。
ＲＹＢ設定量（配合量）に１を加えるのは、設定量を０とすることもあるので、対数関数の引数が０になることが無いようにするための、便宜上の処理である。また、１を加えた後でその対数をとる理由は次の通りである。まず１つは、ランベルト・ベールの法則から、ある物質の透過率は  ｅｘｐ（−αｘ）で表わされるためである。ここでαは吸収係数でｘは経路長である。レンズの染色においては染色の程度を表す量がｘに相当するので、対数をとることによって設定量と透過率の関係を線形に近くできる。その結果、近似関数をより精度よく作ることができる。もう１つの理由は、設定量の数値が大きいデータと小さいデータに関して、誤差の比率を平等にして最適化するためである。もし対数をとらなければ、数１９の式で与えられるＥへの寄与は、設定量の数値が大きいデータのほうが大きくなってしまうため、設定量の小さいデータに関して誤差の割合が大きくなり、うまく最適化することができない。
以上の理由から参照値はＲＹＢ設定量に１を加えた値の対数とするので、近似関数の３個の出力値の指数から１を引いてＲＹＢの設定量を得る。
【００５１】
【数２１】
2018191073.tif 000023
【００５２】
近似関数を最適化するためには、全データに関する残差２乗和を最小化するが、近似関数の更新は個別のデータについて逐次的に行う。すなわち、第ｐ番目のデータに関して、下記数２２のＥｐを最小化する。その更新手順を全データに関して繰り返し行うことによって、近似関数を最適化する。その手順においては、重要なデータ（類似のデータが少ない稀なデータ）の使用頻度を多くしたり、確率的にランダムにデータを選択したりする工夫をしてもよい。
【００５３】
【数２２】
2018191073.tif 000024
【００５４】
まず、第４層のノード７個から第５層（出力層）の３個のノードへの２１個の重みを更新する計算方法を示す。数２３は重みを表す一般式である（但し、ｊ＝１〜３、ｋ＝１〜７）。尚、以下の計算は複雑であるためコンピュータで実行される。
【００５５】
【数２３】
2018191073.tif 000025
【００５６】
Ｅｐを特定の重み（第４層ｋ番目のノードから第５層ｊ番目のノードへの重み）で微分したときの値（微分値）を求める
【００５７】
【数２４】
2018191073.tif 000026
【００５８】
ここで数２４の右辺の左側は、数２５のようにＥｐを実際に偏微分することで、「近似関数のｊ番目の出力」と「出力するべき目標値」との差であることがわかる。「近似関数の出力」は第ｐ番目の測色データを近似関数に入力して得られる。「出力するべき目標値」は第ｐ番目の数値データに１を加えて対数をとった値である。
【００５９】
【数２５】
2018191073.tif 000027
【００６０】
また、数２４の右辺の右側は、数２６のように「第５層ｊ番目のノードの出力値の変化量」の「第４層ｋ番目のノードから第５層ｊ番目のノードへの重みの変化量」に対する比である。それは、第４層ｋ番目のノードの値に相当する。
【００６１】
【数２６】
2018191073.tif 000028
【００６２】
こうして、Ｅｐを特定の重み（第４層ｋ番目のノードから第５層ｊ番目のノードへの重み）で微分したときの値を求めることができる。その値に学習係数ηをかけて、重みから引くことにより更新する（数２７）。学習係数ηは近似関数を構成するすべての重みについて共通に適用する値であり、典型的には１０の−３乗から−８乗といった小さい値を用いる。近似関数を最適化する学習の初期においては学習係数ηを大きな値とし、最適化が進んだら小さくするといった工夫をしてもよい。
【００６３】
【数２７】
2018191073.tif 000029
【００６４】
次に、第３層のノード１０個から第４層の７個のノードへの７０個の重みを更新する計算方法を示す。下記数２８は重みを表す一般式である（但し、ｊ＝１〜７、ｋ＝１〜１０）。
【００６５】
【数２８】
2018191073.tif 000030
【００６６】
数２９のようにＥｐを特定の重み（第３層ｋ番目のノードから第４層ｊ番目のノードへの重み）で微分したときの値（微分値）を求める。
【００６７】
【数２９】
2018191073.tif 000031
【００６８】
数２９の右辺には項が３つある。以下では第１項の値を求める計算についてのみ説明するが、第２項と第３項の計算も同様に行うことができる。
右辺第１項の左側は「第５層１番目のノードの値」と「第１番目のノードが出力するべき目標値」との差である。その値を得るための式を数３０として示す。
【００６９】
【数３０】
2018191073.tif 000032
【００７０】
右辺第１項の右側は、「第５層の第１番目のノードの出力値の変化量」の「第３層ｋ番目のノードから第４層ｊ番目のノードへの重みの変化量」に対する比である。それは、次の式で与えられる。
【００７１】
【数３１】
2018191073.tif 000033
【００７２】
数３１において、右辺の左側は、「第５層１番目のノードの値の変化量」の「第４層ｊ番目のノードの値の変化量」に対する比である。その値はＲｅＬＵ関数の出力なので、数３２のように「第４層ｊ番目のノードの値」が負のときは０であり、正のときは「第４層ｊ番目のノードから第５層１番目のノードへの重み」に相当し、０のときは微分値を定義できないが便宜上どちらかを選択して値を決める。
【００７３】
【数３２】
2018191073.tif 000034
【００７４】
また、数３１において、右辺の右側は、「第４層ｊ番目のノードの出力値の変化量」の「第３層ｋ番目のノードから第４層ｊ番目のノードへの重みの変化量」に対する比である。そしてそれは、下記数３３のように第３層ｋ番目のノードの値に相当する。
【００７５】
【数３３】
2018191073.tif 000035
【００７６】
こうして、第３層ｋ番目のノードから第４層ｊ番目のノードへの重みでＥｐを微分したときの値を求めることができる。その値に学習係数ηをかけて、重みから引くことにより更新する（数３４）。
【００７７】
【数３４】
2018191073.tif 000036
【００７８】
以下、計算方法が複雑になっていくが、第２層のノードから第３層のへの重み、第１層のノードから第２層のへの重みを更新する式を同様の計算で得ることができる。但し、入力ノードから第１層のへの重みを更新する際は、入力ノードの値にＲｅＬＵ関数を適用しないことを考慮する。
【００７９】
上記実施の形態は本発明の原理およびその概念を例示するための具体的な実施の形態として記載したにすぎない。つまり、本発明は上記の実施の形態に限定されるものではない。本発明は、例えば次のように変更した態様で具体化することも可能である。
・上記実施の形態では物体の例として透明な眼鏡レンズを挙げたが、非透明な物体の表面に着色するようにしてもよい。図３はそのような場合のＤＮＮによるネットワーク図である。ここでは上記実施の形態と同様に入力値を決めて、出力値としての４つの色Ｃ（シアン）Ｍ（マゼンタ）Ｙ（イエロー）Ｋ（ブラック）のプリンターで使用する着色剤の配合量を算出するというものである。ここでは、実施の形態に倣って見本について波長ごとの３５個の分光反射率（分光反射率は分光反射率測定装置によって測定）と、紙やインクの種類等のカテゴリ値（変数）を数値データとして入力値とし、重みを適用してＤＮＮを行って出力値としてＣ、Ｍ、Ｙ、Ｋの各色の配合量を算出する。また、上記実施の形態に倣ってＤＮＮの重みを更新する。また、上記実施の形態ではドライ染色（第１の色付け処理）に対して追加的な調色（第２の色付け処理）としてウェット調色を採用したが、印刷においての追加的な調色では、例えばプリント後に追加印刷で調色するようにする。
尚、印刷装置を制御する条件を表す量的変数を入力に追加してもよい。
・上記実施の形態で使用した近似関数では、重みの値を合計で１２９５個とした。これは若干多いため適宜重みを減らすようにしてもよい。しかし、パラメタの最適化を無理なく行うには、参照データは少なくとも５００〜１０００個は必要である。少ないデータを用いて近似関数を作ると過学習を生じるが、その場合は参照データとの違いが大きい特殊な入力値（参照データに含まれないような分光透過率パターン）から、異常な値が出力されてしまうことがあるためである。この問題に対応するためには、使用可能なデータが少ないうちは重みの数が少ないような近似関数を作るとよい。すなわち、５層よりも少ない層数にしたり、あるいは測定波長を２０ｎｍ間隔にして１８個とする方法が考えられる。初期は小さなモデルで受注・加工・測定を行い、データが蓄積されたらモデルを変えて計算の精度を高めればよい。
・上記実施の形態ではＤＮＮは５層で計算したが、５層よりも多層にしてもよい。
・上記実施の形態の入力値の条件は適宜変更可能である。例えば、レンズ基材のカテゴリ数を増やすこと、カテゴリをレンズ基材に限定せずに、コート別に設定することなども自由である。
・実施の形態のレンズの染色におけるドライ加工は赤（Ｒ）黄（Ｙ）青（Ｂ）の３色で行う例を示したが、黒（Ｋ）緑（Ｇ）シアン（Ｃ）マゼンタ（Ｍ）などを加えて４色以上で色付けを行ってもよい。
【００８０】
図６〜図８の入力値データを近似関数の入力として使用する際、見本レンズの度数を常に０としてもよい。それは、これら入力値データは、染色加工したレンズの分光透過率を測定したものであり、近似関数の出力は「結果的にレンズの色がそのように見える」条件のＲＹＢ設定量なので、入力・出力ともに見本レンズの度数とは無関係だからという考え方に基づく。
また、図６〜図８の入力値データを近似関数の入力として使用する際、見本レンズの度数を加工するレンズの度数と同じ値としても良い。それは、前回のドライ染色によってできたレンズと同じ色に見える同度数のレンズを新たな加工において製造するとしたら、前回のドライ染色によってできたレンズを見本レンズとするので、前回の図６の入力値データをそのまま図５の入力値データとして用いればよく、前回のドライ染色加工にて用いたＲＹＢ設定量を参照値とするべきという考え方に基づく。
・実施の形態ではレンズ全面を同じ色に染める全面加工の例を示したが、部分的な加工に適用するようにしてもよい。部分的な加工とは例えばハーフ加工がある。これは一般にレンズ上部を濃いカラーに、レンズ下部を薄いカラーに染める方法である。色の濃さが上部から下部にかけて徐々に変わる（グラデーション）効果がある。これをドライ染色で行うには、用紙の領域別に色の濃さを変えてプリントする。徐々に変える方法や、３〜５段階程度に分けて数値を設定する方法により加工できる。ウェット染色で行うには、レンズを上部のみを染色液につけた状態で、上下に揺動する方法により加工できる。
ハーフ加工に本発明を適用するには、レンズの上部と下部の所定の位置でそれぞれ分光透過率測定を行い、分光透過率から加工条件を算出する近似関数の最適化を、上部と下部それぞれで行うとよい。ある見本レンズについて所定位置２ヶ所を分光透過率測定し、２組のＲＹＢ設定量を算出する。そして上中下の３領域にプリントするのであれば、２組の値を上部と下部に適用し、中部は中間的な値とする方法が考えられる。中間的な値とは、たとえば平均値であり、あるいは２つの値の対数の平均値のさらに指数を用いる方法が考えられる。
本願発明は上述した実施の形態に記載の構成に限定されない。上述した各実施の形態や変形例の構成要素は任意に選択して組み合わせて構成するとよい。また各実施の形態や変形例の任意の構成要素と、発明を解決するための手段に記載の任意の構成要素または発明を解決するための手段に記載の任意の構成要素を具体化した構成要素とは任意に組み合わせて構成するとよい。これらについても本願の補正または分割出願等において権利取得する意思を有する。
【図１】
2018191073.tif 000037
【図２】
2018191073.tif 000038
【図３】
2018191073.tif 000039
【図４】
2018191073.tif 000040
【図５】
2018191073.tif 000041
【図６】
2018191073.tif 000042
【図７】
2018191073.tif 000043
【図８】
2018191073.tif 000044