(57)【要約】【課題】定性的な本物感を構成している特徴を明らかにし、定性的な本物感を、物理特徴量を用いて目視と合うように定量的に評価する。【解決手段】対象物の表面を撮像した画像に基づいて前記対象物の本物感を評価する評価装置において、前記対象物の撮像画像内の色度分布に基づく統計量を算出する統計量算出手段と、前記対象物の撮像画像内の表面性状の変動量に基づく空間周波数特性を、前記対象物の表面性状に応じて少なくとも２つ以上の周波数帯に分けて積分し、それぞれの積分値を得る分割積分値取得手段と、本物感に対する主観評価点を用い、前記統計量および前記積分値に基づいて前記対象物の本物感を示す評価値を算出する評価値算出手段と、を備える。【選択図】図６
2020153759.tif 000002
【特許請求の範囲】
【請求項１】
対象物の表面を撮像した画像に基づいて前記対象物の本物感を評価する評価装置において、
前記対象物の撮像画像内の色度分布に基づく統計量を算出する統計量算出手段と、
前記対象物の撮像画像内の表面性状の変動量に基づく空間周波数特性を、前記対象物の表面性状に応じて少なくとも２つ以上の周波数帯に分けて積分し、それぞれの積分値を得る分割積分値取得手段と、
本物感に対する主観評価点を用い、前記統計量および前記積分値に基づいて前記対象物の本物感を示す評価値を算出する評価値算出手段と、
を備えることを特徴とする評価装置。
【請求項２】
前記対象物の撮像画像内の色度分布を取得する色度分布取得手段と、
前記対象物の撮像画像内の表面性状を取得する表面性状取得手段と、
を備え、
前記色度分布取得手段と前記表面性状取得手段とは、少なくとも２角度以上の計測角度条件から取得された前記対象物の撮像画像を取得する、
ことを特徴とする請求項１に記載の評価装置。
【請求項３】
前記２角度以上の計測角度条件は、正反射条件に近いハイライト角条件と正反射条件から遠いシェード角条件とを含む、
ことを特徴とする請求項２に記載の評価装置。
【請求項４】
前記表面性状の変動量に基づいて空間周波数特性を算出し、視覚の空間周波数特性で重み付けをする空間周波数特性算出手段を備え、
前記空間周波数特性算出手段は、前記対象物の観察距離に応じて前記視覚の空間周波数特性を変化させ、前記表面性状の変動量の空間周波数特性に対して重み付けする、
ことを特徴とする請求項１ないし３の何れか一項に記載の評価装置。
【請求項５】
前記分割積分値取得手段は、前記少なくとも２つ以上に分ける周波数帯として０．１ｃｙｃｌｅ／ｍｍより波長が長い周波数帯と、０．１ｃｙｃｌｅ／ｍｍより波長が短い周波数帯とを含む、
ことを特徴とする請求項１ないし４の何れか一項に記載の評価装置。
【請求項６】
前記対象物の撮像画像内の色度分布および前記対象物の撮像画像内の表面性状は、Ｌ*ａ*ｂ*表色系におけるＬ*，ａ*，ｂ*値を用いる、
ことを特徴とする請求項１ないし５の何れか一項に記載の評価装置。
【請求項７】
対象物の表面を撮像した画像に基づいて前記対象物の本物感を評価する評価装置における本物感評価方法であって、
前記対象物の撮像画像内の色度分布に基づく統計量を算出する統計量算出工程と、
前記対象物の撮像画像内の表面性状の変動量に基づく空間周波数特性を、前記対象物の表面性状に応じて少なくとも２つ以上の周波数帯に分けて積分し、それぞれの積分値を得る分割積分値取得工程と、
本物感に対する主観評価点を用い、前記統計量および前記積分値に基づいて前記対象物の本物感を示す評価値を算出する評価値算出工程と、
を含むことを特徴とする本物感評価方法。
【請求項８】
対象物の表面を撮像した画像に基づいて前記対象物の本物感を評価する評価装置を制御するコンピュータを、
前記対象物の撮像画像内の色度分布に基づく統計量を算出する統計量算出手段と、
前記対象物の撮像画像内の表面性状の変動量に基づく空間周波数特性を、前記対象物の表面性状に応じて少なくとも２つ以上の周波数帯に分けて積分し、それぞれの積分値を得る分割積分値取得手段と、
本物感に対する主観評価点を用い、前記統計量および前記積分値に基づいて前記対象物の本物感を示す評価値を算出する評価値算出手段と、
として機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、評価装置、本物感評価方法およびプログラムに関する。
【背景技術】
【０００２】
従来、天然皮革の質感を模倣した人工または合成皮革（以下、合皮という）は、安価な衣類、靴等に用いられている。加えて、合皮は、近年の天然皮革価格の高騰などにより、従来は天然皮革が用いられてきた製品にまで採用されるようになっている。そのため、合皮については、より天然皮革に質感を近づけることが求められている。
【０００３】
特許文献１には、合皮の質感を高めることを目的として、スエード調の合皮表面の樹脂塊によって生じる「イラツキ」という質感を、画像特徴量に基づいて評価する技術が開示されている。
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、合皮の質感については、あくまでも人間が天然皮革を見た際に感じる定性的な「らしさ」で判断しているが、人間が感じる天然皮革らしさを構成している特徴についてはまだ明らかになっていない。人間が天然皮革を見た際に感じる定性的な「らしさ」を、「本物感」と呼ぶ。
【０００５】
本発明は、上記に鑑みてなされたものであって、定性的な本物感を構成している特徴を明らかにし、定性的な本物感を、物理特徴量を用いて目視と合うように定量的に評価することを目的とする。
【課題を解決するための手段】
【０００６】
上述した課題を解決し、目的を達成するために、本発明は、対象物の表面を撮像した画像に基づいて前記対象物の本物感を評価する評価装置において、前記対象物の撮像画像内の色度分布に基づく統計量を算出する統計量算出手段と、前記対象物の撮像画像内の表面性状の変動量に基づく空間周波数特性を、前記対象物の表面性状に応じて少なくとも２つ以上の周波数帯に分けて積分し、それぞれの積分値を得る分割積分値取得手段と、本物感に対する主観評価点を用い、前記統計量および前記積分値に基づいて前記対象物の本物感を示す評価値を算出する評価値算出手段と、を備えることを特徴とする。
【発明の効果】
【０００７】
本発明によれば、計測した物理特徴量から、定性的であった本物感を定量的に評価することができる、という効果を奏する。
【図面の簡単な説明】
【０００８】
【図１】図１は、実施の形態に係る評価システムの構成例を示す概略図である。
【図２】図２は、制御部のハードウェア構成を示すブロック図である。
【図３】図３は、制御部の機能を示す機能ブロック図である。
【図４】図４は、主観評価実験のための装置例を示す図である。
【図５】図５は、主観評価試験の結果を示す図である。
【図６】図６は、制御部による本物感評価処理の流れを概略的に示すフローチャートである。
【図７】図７は、本物感と絶対値項の関係性を示す図である。
【図８】図８は、本物感評価量（評価値）と予め行った主観評価点との相関図である。
【発明を実施するための形態】
【０００９】
以下に添付図面を参照して、評価装置、本物感評価方法およびプログラムの実施の形態を詳細に説明する。
【００１０】
図１は、実施の形態に係る評価システム１００の構成例を示す概略図である。図１に示すように、評価システム１００は、第１の画像取得装置１０と、第２の画像取得装置２０と、制御部５と、を備えている。
【００１１】
第１の画像取得装置１０は、サンプルの表面性状に関する特徴量を得る画像を取得する。第１の画像取得装置１０は、投影部の一例である光源１と、サンプル台３と、光源１のサンプル台３に対する照射角を可変とする角度変更部２と、撮像部の一例であるカメラ４と、を備えている。
【００１２】
光源１としては、分光分布が太陽光に近いキセノン光源が用いられている。光源１は、サンプル台３の法線Ｌ方向に対して、角度変更部２により移動されて所望の角度方向に適宜位置づけられる。
【００１３】
サンプル台３には、表面性状を評価するサンプルＯが取り付けられている。カメラ４は、サンプル台３の法線Ｌ方向に配置されている。
【００１４】
カメラ４は、サンプル台３に取り付けられたサンプルＯを撮影して２次元画像を取得する。本実施形態では、カメラ４として一般的なデジタルカメラやＣＣＤ等であるＲＧＢカメラを用いている。
【００１５】
一方、第２の画像取得装置２０は、サンプルの色度分布に関する特徴量を得る画像を取得する。第２の画像取得装置２０は、投影部の一例である光源１１と、サンプル台１３と、光源１１のサンプル台１３に対する照射角を可変とする角度変更部１２と、撮像部の一例であるカメラ１４と、を備えている。
【００１６】
光源１１としては、分光分布が太陽光に近いキセノン光源が用いられている。光源１１は、サンプル台３の法線Ｌ方向に対して、角度変更部２により移動されて所望の角度方向に適宜位置づけられる。
【００１７】
サンプル台１３には、色度分布を評価するサンプルＯが取り付けられている。カメラ１４は、サンプル台１３の法線Ｌ方向に対して４５°傾けた方向に配置されている。
【００１８】
カメラ１４は、サンプル台１３に取り付けられたサンプルＯを撮影して２次元画像を取得する。本実施形態では、カメラ１４としてハイパースペクトルカメラ（分光カメラ）を用いている。ハイパースペクトルカメラとは、光の波長ごとの輝度を計測できるカメラであり、各画素で分光スペクトルを得ることができる。本実施形態の分光カメラは、波長４００ｎｍ〜７００ｎｍの範囲を１０ｎｍ毎に取得することができる。
【００１９】
なお、カメラ１４がハイパースペクトルカメラに限定されることはなく、デジタルカメラやＣＣＤ等のＲＧＢカメラを用いることもできる。ただし、スペクトルではなくではなくＲＧＢの３ｃｈのみの情報からＬ*ａ*ｂ*に変換するため、測色精度は悪くなる。
【００２０】
ここで、図２は制御部５のハードウェア構成を示すブロック図である。図２に示すように、制御部５は、ＣＰＵ５ａ（Central  Processing  Unit）５ａや、ＲＯＭ（Read  Only  Memory）５ｂ、ＲＡＭ（Random  Access  Memory）５ｃ、ＨＤＤ５ｄ等を備えている。制御部５は、ＲＯＭ５ｂやＨＤＤ５ｄに予め記憶されている本物感評価プログラムに従って、ＲＡＭ５ｃをワークメモリとして用いて、光源１、角度変更部２、カメラ４等、第１の画像取得装置１０の各部を駆動制御する。また、制御部５は、ＲＯＭ５ｂやＨＤＤ５ｄに予め記憶されている本物感評価プログラムに従って、ＲＡＭ５ｃをワークメモリとして用いて、光源１１、角度変更部１２、カメラ１４等、第２の画像取得装置２０の各部を駆動制御する。制御部５として、例えば、パーソナルコンピュータ（ディスクトップ、ノートパソコン）を用いることができる。
【００２１】
本実施形態の制御部５で実行される本物感評価プログラムは、インストール可能な形式又は実行可能な形式のファイルでＣＤ−ＲＯＭ、フレキシブルディスク（ＦＤ）、ＣＤ−Ｒ、ＤＶＤ（Digital  Versatile  Disk）等のコンピュータで読み取り可能な記録媒体に記録して提供するように構成してもよい。
【００２２】
さらに、本実施形態の制御部５で実行される本物感評価プログラムを、インターネット等のネットワークに接続されたコンピュータ上に格納し、ネットワーク経由でダウンロードさせることにより提供するように構成してもよい。また、本実施形態の制御部５で実行される本物感評価プログラムをインターネット等のネットワーク経由で提供または配布するように構成してもよい。
【００２３】
制御部５は、カメラ４やカメラ１４で撮影した２次元画像に基づいて、サンプルＯの表面性状および色度分布を評価する。すなわち、制御部５は、評価装置としても機能する。
【００２４】
制御部５の評価装置としての機能について、図３の機能ブロック図を用いて説明する。この図３に示すように、制御部５は、ＣＰＵ５ａが本物感評価プログラムに従って動作することにより、統計量算出手段５１、分割積分値取得手段５２、評価値算出手段５３、色度分布取得手段５４、表面性状取得手段５５、空間周波数特性算出手段５６、記憶部５７として機能する。
【００２５】
色度分布取得手段５４は、対象物となるサンプルＯの撮像画像内の色度分布を取得する。また、表面性状取得手段５５対象物となるサンプルＯの撮像画像内の表面性状を取得する。色度分布取得手段５４と表面性状取得手段５５とは、少なくとも２角度以上の計測角度条件から取得された対象物となるサンプルＯの撮像画像を取得する。
【００２６】
統計量算出手段５１は、対象物となるサンプルＯの撮像画像内の色度分布に基づく統計量を算出する。
【００２７】
空間周波数特性算出手段５６は、対象物となるサンプルＯの撮像画像内の表面性状の変動量に基づいて空間周波数特性を算出し、視覚の空間周波数特性で重み付けをする。空間周波数特性算出手段５６は、対象物となるサンプルＯの観察距離に応じて視覚の空間周波数特性を変化させ、表面性状の変動量の空間周波数特性に対して重み付けする。
【００２８】
分割積分値取得手段５２は、対象物となるサンプルＯの撮像画像内の表面性状の変動量に基づく空間周波数特性を、対象物となるサンプルＯの表面性状に応じて少なくとも２つ以上の周波数帯に分けて積分し、それぞれの積分値を得る。
【００２９】
評価値算出手段５３は、本物感に対する主観評価点を用い、対象物となるサンプルＯの撮像画像内の色度分布に基づく統計量および積分値に基づいて対象物となるサンプルＯの本物感を示す評価値を算出する。
【００３０】
記憶部５７には、予め取得された官能評価情報が記憶される。
【００３１】
本発明の評価システム１００で本物感評価を行うにあたり、サンプルＯの物理的量計測のほかに、尺度化された官能評価点（主観評価点）が必要となる。そこで、予め本物感の主観評価実験を行い、実験によって得られた官能評価点（主観評価点）を、官能評価情報として記憶部５７に記憶する。
【００３２】
以下、主観評価実験について説明する。主観評価実験には、一対比較法と呼ばれる手法を用いた。一対比較法は、被験者に対し、評価対象の物体を１対ずつランダムで提示し、例えば左の物体に対して、右の物体は本物感があるかどうか被検者に点数を付けて貰う手法である。
【００３３】
このような一対比較法であるシェッフェ法による主観評価実験を行うため、天然皮革と合皮の計１０種類のサンプルを作成し、用意した。下記に示す表１に、今回用いたサンプル一覧を示す。
【００３４】
【表１】
2020153759.tif 000003
【００３５】
今回用いたサンプルの色はすべて最も流通量の多い黒、サンプルの大きさはＡ４サイズ（縦：約３００ｍｍ，横：約２１０ｍｍ）とした。図４は、主観評価実験のための装置例を示す図である。図４に示すように、皮革は観察する角度によって質感が変化するため、冶具５０にサンプルＯを湾曲させて設置した。観察距離は、３００ｍｍとした。被験者は、提示された湾曲された状態の一対の左側サンプルＡ及び右側サンプルＢを比較し、左側サンプルＡに対する右側サンプルＢの表面から視覚的に認知される本物感に基づいて、下記のように点数を付与した。
【００３６】
ＡよりもＢの方が非常に強く本物感を感じたとき      ＋３点
ＡよりもＢの方が強く本物感を感じたとき            ＋２点
ＡよりもＢの方が本物感を感じたとき                ＋１点
Ａ＝Ｂのとき                                        ０点
ＢよりもＡの方が本物感を感じたとき                −１点
ＢよりもＡの方が強く本物感を感じたとき            −２点
ＢよりもＡの方が非常に強く本物感を感じたとき      −３点
【００３７】
なお、本実験では７段階で点数付けをさせたが、例えば９段階、１１段階等、そのほかの段階数でもよい。段階数を増やすほど、より微小な差を検出できるが、被験者への負担が大きくなる。
【００３８】
被験者は、３００ｍｍ離れた位置から評価を行った。本実験では、１０種類の観察条件があるため全部で１０Ｃ２＝４５対を評価させた。被験者として画像評価技術者計８名に対して実験を実施した。各サンプル物体に対する評価点は、上記評価実験の結果に数量化ＩＩＩ類の解析手法を適用し算出した。
【００３９】
評価結果を図５に示す。図５の横軸は、サンプル物体に付した番号（通番）であり、サンプル物体を主観評価点順に並べている。縦軸は、各サンプル物体の主観評価点を示しており、値が小さいほど人工的に感じられることを意味している。つまり、サンプル物体「１」が最も本物に、サンプル物体「４」が最も人工的に感じられるといった評価結果になった。
【００４０】
以下、上記構成の評価システム１００により実行される本物感評価方法の動作について説明する。
【００４１】
まず、第１の画像取得装置１０で取得する画像について説明する。
【００４２】
評価対象のサンプルＯは、サンプル台３にセットされる。第１の画像取得装置１０は、光源１を点灯し、サンプル台３に取り付けられたサンプルＯに対して光を照射し、サンプル台３の法線方向Ｌに設置したカメラ４で計測を行う。皮革は、観察する角度によって質感が変化することから、光源１は、角度変更部２によって移動され、カメラ４に対してハイライト角とシェード角に傾けて照明した。ここで、ハイライトとはカメラ４（サンプルＯの法線方向Ｌ）に対して１０〜２５°傾けた範囲、シェードとはカメラ４（サンプルＯの法線方向Ｌ）に対して４５〜６０°傾けた範囲を指す。このように正反射条件に近いハイライト角条件と正反射条件から遠いシェード角条件とを計測角度条件とすることにより、角度変化による質感変化がより大きい２角度の特徴量を用いて評価することができる。本実施形態では、１５°と６０°での計測値を用いた（ジオメトリ：１５°／０°、６０°／０°）。これにより、観察角度による質感の変化を考慮することができる。第１の画像取得装置１０で取得した画像は、計測サイズ：３０００×３０００ｐｉｘｅｌ（撮像範囲：約７５ｍｍ×７５ｍｍ）、解像度：約１０００ｄｐｉである。
【００４３】
次に、第２の画像取得装置２０で取得する画像について説明する。
【００４４】
評価対象のサンプルＯは、サンプル台１３にセットされる。第２の画像取得装置２０は、光源１１を点灯し、サンプル台１３に取り付けられたサンプルＯに対して光を照射し、サンプル台１３の法線方向Ｌから４５°傾けて設置したカメラ１４で計測を行う。皮革は、観察角度によって色が変化することから、光源１１は、角度変更部２によって移動され、カメラ１４と正反射角および正反射角に対して法線方向Ｌにハイライト角、シェード角と傾けて照明した。ここで、ハイライトとはサンプルＯの法線方向Ｌに対して２０〜３５°傾けた範囲、シェードとはサンプルＯの法線方向Ｌに対して０〜−１５°傾けた範囲を指す。このように正反射条件に近いハイライト角条件と正反射条件から遠いシェード角条件とを計測角度条件とすることにより、角度変化による質感変化がより大きい２角度の特徴量を用いて評価することができる。本実施形態では、ハイライトは３０°、シェードは−１５°での計測値を用いた（ジオメトリ：４５°／４５°、３０°／４５°、−１５°／４５°）。これにより、観察角度による質感の変化を考慮することができる。第２の画像取得装置２０で取得した画像は、計測サイズ：６００×６００ｐｉｘｅｌ（撮像範囲：約１５ｍｍ×１５ｍｍ）、解像度：約１０００ｄｐｉである。
【００４５】
続いて、制御部５により実行される本物感評価処理の流れについて説明する。
【００４６】
ここで、図６は制御部５による本物感評価処理の流れを概略的に示すフローチャートである。図６に示すように、制御部５は、第１の画像取得装置１０および第２の画像取得装置２０を用いて計測画像（撮像画像）を取得する（ステップＳ１）。
【００４７】
次に、制御部５は、計測画像（撮像画像）をＬ*ａ*ｂ*表色系におけるＬ*ａ*ｂ*画像に変換し、色度分布を取得する（ステップＳ２）。より詳細には、制御部５は、Ｌ*ａ*ｂ*を算出するために、下記に示す式（１）、式（２）を用いてＲＧＢ画像または分光画像をＸＹＺ画像に変換する。
【００４８】
【数１】
2020153759.tif 000004
【００４９】
【数２】
2020153759.tif 000005
【００５０】
ここで、Ｓ（λ）は照明の分光分布、ｘ（λ）、ｙ（λ）、ｚ（λ）は等色関数、Ｒ（λ）は分光反射率、ｋは係数である。ｋは、下記に示す式（３）から算出される。なお、分光反射率は、反射率が１００％で完全拡散放射に近い反射特性を持つ白色基準板の撮影データを用いて、サンプルの撮影データを規格化し、反射率化した。
【００５１】
【数３】
2020153759.tif 000006
【００５２】
そして、制御部５は、下記に示す式（４）を用いてＸＹＺ画像からＬ*ａ*ｂ*画像に変換する。
【００５３】
【数４】
2020153759.tif 000007
【００５４】
ここで、Ｘｎ、Ｙｎ、Ｚｎは国際照明委員会（ＣＩＥ）で定められた完全拡散反射面での３刺激値であり、Ｘｎ＝９６．４２、Ｙｎ＝１００、Ｚｎ＝８２．４９である。なお、本実施形態では、Ｓ（λ）はＤ６５の分光分布とした。このように人間の感覚に合う等色系に変換することで、目視評価との相関を向上させることができる。
【００５５】
図６に戻り、制御部５（統計量算出手段５１）は、第２の画像取得装置２０で得られたＬ*ａ*ｂ*画像のうち、Ｌ*画像を用いて尖度と歪度を算出する（ステップＳ３）。尖度と歪度は統計量の一種であり、母集団が正規分布からどれだけずれているかを表す指標である。
【００５６】
尖度は、分布が正規分布からどれだけ尖っているかを表す統計量で、山の尖り度と裾の広がり度を示す。正規分布より尖った分布のときには正の値を、正規分布より扁平な分布のときには負の値となる（正規分布の場合には０）。
【００５７】
一方、歪度は分布が正規分布からどれだけ歪んでいるかを表す統計量であり、左右対称性を示す指標である。右裾が長いもしくは左に偏った分布のときには正の値を、左裾が長いもしくは右に偏った分布のときには負の値をとる（左右対称の分布の場合には０）。
【００５８】
次に、制御部５（空間周波数特性算出手段５６）は、第１の画像取得装置１０で得られたＬ*ａ*ｂ*画像のうち、Ｌ*画像を用いて空間周波数特性を算出する（ステップＳ４）。より詳細には、制御部５（空間周波数特性算出手段５６）は、Ｌ*画像の各画素における値からＬ*画像平均値を差し引き、Ｌ*偏差画像を取得する。
【００５９】
Ｌ*の偏差ΔＬ*は、全体的に反射する金属のようなサンプルでは平均値が上昇するために一様に低い値をとり、平均値に影響の少ない程度の小さな部位が局所的に特に強く反射する、所謂きらきらと光って見える状態のときには、光って見える画素の部分において大きい値をとる。すなわち、このように局所的に反射が強くなる画素が、シボやムラとして人の目には感知される。制御部５（空間周波数特性算出手段５６）は、取得したＬ*偏差画像を全ての画素における偏差Δについてフーリエ変換する。制御部５（空間周波数特性算出手段５６）は、フーリエ変換後、振幅を各周波数で平均化して、１次元化された空間周波数特性を算出する。
【００６０】
次いで、制御部５（空間周波数特性算出手段５６）は、ステップＳ４で得られた空間周波数特性に対して、「観察距離に対する視覚の空間周波数特性（ＶＴＦ）」で重み付けを行う（ステップＳ５）。本実施形態では、視覚の空間周波数特性として、下記に示す式（５）に示す関数を用いている。なお、視覚の空間周波数特性としては、既に知られた他の視覚特性についての関数を利用しても良い。
【００６１】
【数５】
2020153759.tif 000008
【００６２】
ここでυの単位は、ｃｙｃｌｅ／ｄｅｇｒｅｅである。ｃｙｃｌｅ／ｄｅｇｒｅｅは、式（５）によって求めることができることから、サンプルの観察距離に応じて計算した式（５）を１次元化された空間周波数特性に掛け合わせることで、視覚特性を考慮したシボやムラの評価を行うことができる。これにより、想定する観察距離による人間の眼の解像度に応じて、空間周波数特性を算出することができる。
【００６３】
次いで、制御部５（分割積分値取得手段５２）は、ステップＳ５で取得した重み付け空間周波数特性を、皮革の表面性状に応じていくつかの周波数帯に分割する（ステップＳ６）。本実施形態では、使用する皮革サンプルを元に３つの周波数帯に分割した。１つ目は、０〜０．１ｃｙｃｌｅ／ｍｍであり、これは皮革のムラに相当する。２つ目は、０．１〜１ｃｙｃｌｅ／ｍｍであり、これは粗いシボに相当する。最後に、３つ目は１〜４ｃｙｃｌｅ／ｍｍであり、これは細かなシボに相当する。このようにすることで、分ける周波数帯の範囲をある程度規定できるとともに、それぞれの表面性状の特徴に対応した積分値を算出することができる。
【００６４】
次いで、制御部５（分割積分値取得手段５２）は、ステップＳ６にて分割した周波数帯に対して、積分値取得手段として、それぞれ積分を行う（ステップＳ７）。
【００６５】
最後に、制御部５（評価値算出手段５３）は、主観評価点と、ステップＳ２、ステップＳ３、ステップＳ６それぞれで算出した特徴量を用いて、以下に示す式（６）にて本物感評価値を算出する（ステップＳ８）。評価式の算出には、ステップワイズ法を用いた回帰分析を使用した。しかし、ステップワイズ法では、決定係数が最良になるモデルを出力するため、評価式の意味、つまり、説明変数選択のロジックに関しては考慮されていない。そこで、ステップワイズ法の出力に対して、統計的有意を保持しつつ、意味が通る評価式になるよう調整した。
【００６６】
【数６】
2020153759.tif 000009
【００６７】
下記の表２に示すように、ｘ１、ｘ２、ｘ３、ｘ４、ｘ５、ｘ６はそれぞれ説明変数であり、４５°／４５°Ｌ*、６０°／０°空間周波数特性（０．１〜１ｃｙｃｌｅ／ｍｍ）、−１５°／４５°尖度、１５°／０°空間周波数特性（０〜０．１ｃｙｃｌｅ／ｍｍ）、６０°／０°空間周波数特性（０〜０．１ｃｙｃｌｅ／ｍｍ）、６０°／０°空間周波数特性（１〜４ｃｙｃｌｅ／ｍｍ）である。
【００６８】
【表２】
2020153759.tif 000010
【００６９】
また、下記の表３に示すように、ｐ１、ｐ２、ｐ３、ｐ４、ｐ５、ｐ６、ｐ７、ｐ８、ｐ９は係数である。係数は、それぞれ主観評価点に基づいて決定されるパラメータである。そのため、主観評価実験の評価段階を例えば７段階から９段階に変えて主観評価点を算出した場合、パラメータは変化する。本実験では、下記の表３に示すように、ｐ１＝０．３、ｐ２＝−０．０２、ｐ３＝５．３７、ｐ４＝−１．０６、ｐ５＝−０．０６、ｐ６＝０．６９、ｐ７＝１３．２１、ｐ８＝−０．５、ｐ９＝−０．５（ｐ＜０．０５）である。
【００７０】
【表３】
2020153759.tif 000011
なお、式（６）のパラメータは、説明変数を標準化して決定してもよい。標準化を行うことにより、単位や平均値などが異なるデータ同士を単純に比較できるようになる。
【００７１】
式（６）から本物感は、シェード条件で撮影した時の０．１〜１ｃｙｃｌｅ／ｍｍの周波数帯、つまり、サンプル表面上の粗いシボに相当する特徴量と正反射条件での明度値とを足し合わせた値の絶対値と、ハイライト条件とシェード条件で撮影した時の０〜０．１ｃｙｃｌｅ／ｍｍの周波数帯、つまり、サンプル表面上のムラに相当する特徴量がそれぞれ大きいほど、値は大きくなることがわかる。
【００７２】
一方、シェード条件で撮影した時の１〜４ｃｙｃｌｅ／ｍｍの周波数帯、つまり、サンプル表面上の細かなシボに相当する特徴量が大きいと値は小さくなることがわかる。シェードの粗いシボと正反射条件での明度の項が絶対値となっているのは、主観評価結果から、絶対値項の正負の値が大きいほど、本物感の点数が大きくなることが分かったためである。
【００７３】
図７に本物感と絶対値項の関係性を示す。図７の横軸は、サンプル物体に付した番号（通番）であり、サンプル物体を主観評価点順に並べている。縦軸は、各サンプル物体の主観評価点（本物感）および絶対値項を示している。絶対値項の正負の値が大きいほど、本物感の点数が大きくなっている。
【００７４】
上述のような本物感評価方法によって算出された本物感評価量（評価値）と予め行った主観評価点との相関図を図８に示す。寄与率Ｒ２は０．９２となり、本実施形態で得られた評価量が、人間の主観と極めて近いことがわかる。以上のことから、本発明により、目視と相関の高い本物感評価方法を開発できた。
【００７５】
このように本実施の形態によれば、計測した物理特徴量から、定性的であった本物感（天然皮革らしさ）を定量的に評価することができる。要するに、従来評価法とは異なり、あらかじめ取得した本物感に対する主観評価点に対して、皮革サンプルに対して多角度条件計測で得た物理特徴量から、本物感の評価式を構築することで、今まで定性的であった本物感を定量的に評価することができる。
【００７６】
以上、本発明の実施形態を図面により詳述してきたが、上記各実施形態は本発明の例示にしか過ぎないものであり、本発明は上記各実施形態の構成にのみ限定されるものではない。本発明の要旨を逸脱しない範囲の設計の変更等があっても、本発明に含まれる。
【符号の説明】
【００７７】
５      評価装置
５１    統計量算出手段
５２    分割積分値取得手段
５３    評価値算出手段
５４    色度分布取得手段
５５    表面性状取得手段
５６    空間周波数特性算出手段
【先行技術文献】
【特許文献】
【００７８】
【特許文献１】特開２０１７−０１５４９６号公報
【図１】
2020153759.tif 000012
【図２】
2020153759.tif 000013
【図３】
2020153759.tif 000014
【図４】
2020153759.tif 000015
【図５】
2020153759.tif 000016
【図６】
2020153759.tif 000017
【図７】
2020153759.tif 000018
【図８】
2020153759.tif 000019