(57)【特許請求の範囲】
【請求項１】
毛髪に直線偏光を照射し、毛髪を通過した光を、照射光の偏光面に平行な透過軸を有する偏光子をかけて照射光の光軸上でＣＣＤカメラを用いて検出することによって、毛髪からの散乱光が除かれた透過光の強度に依存した画像を得るとともに、かかる散乱光が除かれた透過光の強度を測定し、毛髪内の特定の物質の存在を画像として得ることを特徴とする毛髪内部の評価方法。
【請求項２】
照射光が、波長２００ｎｍ〜２０００ｎｍの範囲を含む光である請求項１に記載の毛髪内部の評価方法。
【請求項３】
毛髪を互いの透過軸が平行である２枚の偏光子で挟み、偏光子に垂直方向から光を照射する請求項１又は請求項２記載の毛髪内部の評価方法。
【請求項４】
毛髪内部に含まれる物質の吸光係数が異なる複数の波長によって、透過光の強度を測定し、演算を行うことにより、その物質のみを強調した画像を得る請求項１ないし請求項３の何れかの請求項記載の毛髪内部の評価方法。
【請求項５】
吸光係数が異なる複数の波長が、吸光係数が異なる２種の波長である請求項４記載の毛髪内部の評価方法。
【請求項６】
時間を追って特定の波長の透過光の強度を測定し、毛髪内部に含まれる特定物質の量の時間変化を画像でモニターする請求項１ないし請求項５の何れかの請求項記載の毛髪内部の評価方法。
【請求項７】
毛髪を回転させながら透過光の強度を測定し、毛髪内部の３次元構造を得る請求項１ないし請求項６の何れかの請求項記載の毛髪内部の評価方法。
【請求項８】
毛髪を回転させながら測定して得た透過光の複数の強度データからなる連立方程式を、ムーアーペンローズ（Ｍｏｏｒｅ−Ｐｅｎｒｏｓｅ）型一般逆行列を用いて解き、毛髪内の各微小部分の光透過係数を得る請求項７記載の毛髪内部の評価方法。
【請求項９】
毛髪内の各微小部分の光透過係数の値から、毛髪の特定断面上の微小部分の光透過係数を抽出することによって、毛髪を破壊することなくその断面の画像を得る請求項７又は請求項８記載の毛髪内部の評価方法。
【発明の詳細な説明】
【技術分野】  本発明は、対象物体に光を照射し、それを通過した光を検出することによって対象物体内部を評価する方法に関する。詳細には、特定の方法で散乱光の影響を除いた透過光を抽出し、対象物体内部に含まれる特定の物質に着目した情報を得ることができる方法に関し、更に詳細には、対象物体を回転させながら透過光の強度を測定し、対象物体を破壊することなしに、その内部の３次元構造を知ることができる方法に関する。
【背景技術】  従来より、物体内部の情報を得る手段としては、物体を薄く切断し透過型電子顕微鏡（ＴＥＭ）を用いて観察したり、切断面を光学顕微鏡や走査型電子顕微鏡（ＳＥＭ）を用いて観察する方法が用いられてきた。しかし、これらの方法は対象物体を破壊することが必要であり、対象物体を元に戻せないだけでなく、時間とともに対象物体の変化を追うこともできなかった。
これを解決するものとして、Ｘ線断層撮影法（Ｘ線ＣＴ）を用いる方法が知られているが、やはり対象物体を侵襲するものであり、装置としても大きくならざるを得なかった。
一方、微小部分の対象物体内部の情報を非破壊で、すなわちリアルタイムで得たいという要望が高まっており、特に対象物体が毛髪である場合にその要望が顕著である。
すなわち、近年ヘアカラーに対する関心が高まり、多彩な色を備えたカラーリング剤や、カラーリングの際に毛髪が受けるダメージを少なくしたカラーリング剤や、毛髪内部を修復するヘアトリートメント剤等の開発の重要性が高まっている。そのため、実際に毛髪内部で起こっているカラーリング剤の脱色・染色や、ヘアートリートメントによる補修プロセスを明らかにすることは、製品を開発する上で重要な技術である。
そして通常、カラーリング剤やヘアートリートメント剤が毛髪に与える影響は時間依存性があるため、その影響をリアルタイムに測定する必要があった。そこで、非破壊、非侵襲でこれを実現するものとして、通常の光を光源として用いたものが原理的には考えられるが、物体内部で散乱する光と、直線的に透過する光とを効率よく分離する方法が確立されていないために、分解能が低く問題点のあるものであった。
従って、物体内部で散乱する光を除き、直線的に透過する光の強度を測定して画像化する方法の開発が望まれていた。さらにはこれにより、光を利用して微小部分の情報を得たり、特定物質を強調した画像情報を得たり、光ＣＴによる対象物体内部の３次元構造を得る方法の開発が望まれていた。
【発明の開示】
本発明者は、上記課題を解決するために、鋭意研究を行ったところ、対象物体に直線偏光を照射し、対象物体を通過した光を、透過軸が照射光の偏光と平行であるような偏光子をかけて測定したり、光の波長を変化させて対象物体を通過した光の強度を測定することによって、散乱光と直線透過光を効率よく分離できることを知った。
そして、散乱光の影響が除かれた透過光の強度を、波長を変えて測定したり、時間を追って測定したり、対象物体を回転させて測定したりすることによって、それぞれ特定物質に着目した情報が得られたり、リアルタイムの情報が得られたり、３次元のすなわち対象物体内部の情報を非破壊で得られることを見出し本発明に至った。
すなわち、本発明は、対象物体に直線偏光を照射し、対象物体を通過した光を、透過軸が、照射光の偏光と平行であるような偏光子をかけて照射光の光軸上でＣＣＤカメラを用いて検出することによって、対象物体からの散乱光が除かれた透過光の強度に依存した画像を得る対象物体内部の評価方法を提供するものである。
また、対象物体を回転させながら散乱光の影響が除かれた透過光の強度をＣＣＤカメラで測定し、対象物体内部の３次元構造を得る対象物体内部の評価方法を提供するものである。
【図面の簡単な説明】  図１は、本発明の評価方法に用いられる装置の概略図である。
図２は、対象物体（試料）セット用ホルダーを上から見た図面である。
図３は、図２のＡ−Ａ’断面図である。
図４は、対象物体中の特定物質３種の吸光係数波長依存性の例を示す図である。
図５は、対象物体を回転させながら測定し、３次元情報を得るアルゴリズム説明図である。
図６は、対象物体を回転させながら測定し、３次元情報を得るアルゴリズム説明図である。
図７は、実施例１における、各ブリーチ時間毎の撮像結果を示す図である。
図８は、実施例１における、ブリーチ時間＝０分のＲＧＢの輝度分布（透過光強度）を示す図である。
図９は、実施例１における、各ブリーチ時間毎のメラニン色素の分布を示す図である。
図１０は、実施例２における、各ブリーチ時間毎の、毛髪内部のメラニン色素のうち、脱色反応によって吸光係数が変化したものの分布を示す図である。
図１１は、実施例３で使用された毛髪回転装置を示す図である。
図１２は、実施例３の説明で使用される３次元座標系とボクセル分割を示す図である。
なお、図面中の符号は以下の通りである。
１ａ・・・・・・・・照射光側偏光板
１ｂ・・・・・・・・ＣＣＤカメラ側偏光板
１ｄ、１ｅ・・・・・互いに透過軸を平行にした偏光子
２・・・・・・・・・光源
３・・・・・・・・・対象物体
４・・・・・・・・・透過光用ＣＣＤカメラ
１０・・・・・・・・・対象物体セット用ホルダー
１１・・・・・・・・・プレパラート
１３ａ、１３ｂ・・・・固定用ガラス
１４・・・・・・・・・空隙
２０・・・・・・・・・毛髪回転装置
２１・・・・・・・・・毛髪
２２・・・・・・・・・回転部
２３・・・・・・・・・回転部受け
【発明を実施するための最良の形態】  本発明は、対象物体に直線偏光を照射し、対象物体を通過した光を、透過軸が、照射光の偏光と平行であるような偏光子をかけて照射光の光軸上で検出することによって、対象物体からの散乱光を除くことによって、散乱光の影響が除かれた透過光の強度を測定するものである。散乱光は波長依存的に緩やかに変化するが、散乱光の影響が除かれた透過光の波長依存性は時に大きく変化するので、散乱光の影響が除かれた透過光の強度を測定することで、対象物体内の特定の物質の存在を画像として強調して得ることもできる。
図１において詳述すると、光源（２）から出た光は偏光子（１ａ）を通過して直線偏光となり、対象物体（３）に照射される。照射光は、対象物体内部で散乱、吸収を受けて、後部の光軸上に設けた偏光子（１ｂ）を通りＣＣＤカメラ（４）で検出される。このとき、偏光子（１ａ）と偏光子（１ｂ）を、透過軸が平行になるように設置すると、ＣＣＤカメラ（４）には対象物体内部での散乱光が除かれた透過光が主に入射することになり、散乱光の影響が除かれた透過光による画像を得ることができる。
さらに、図２及び図３に示したように、互いの透過軸が平行な２枚の偏光子（１ｄ）、（１ｅ）で対象物体を挟み、そこに光照射する方法も対象物体が毛髪等の小さいものの場合に特に好ましい。
また、散乱光の強度は一般に波長依存性が小さく、散乱光以外の透過光は物質により、その波長依存性が大きいことを利用して、波長を変化させて対象物体を通過した光強度を測定し、さらに散乱光を演算で除くことで散乱光の影響が除かれた透過光のデータとすることもできる。波長を変化させて対象物体を通過した光強度を測定する方法としては、対象物体に照射する光の波長を変化させる方法と白色光を照射して検出する方の光の波長を変化させる方法があり、何れも好適に用いられる。ここで、ＣＣＤカメラを用いる場合には、赤、緑、青の３チャンネルで一度に光強度を測定できるので後者の方法が好ましい。
照射光としてＸ線を用いることはすでに知られており、本発明はＸ線より長波長の光を用いることを特徴とする。照射光の波長としては、２００ｎｍ〜２０００ｎｍであることが好ましい。特に好ましくは、可視光領域の波長である。
本発明では、対象物体を通過した光をＣＣＤカメラで検出し、主に散乱光の影響が除かれた透過光に依存した結果を２次元イメージとして表示するが、その際ＣＣＤカメラの前に対物レンズを設置することが、小さい対象物体の透過率を測定できるために好ましい。対物レンズの倍率には、特に限定はないが５倍〜８００倍とすることが好ましい。
本発明においては、さらに、対象物体内部に含まれる物質の吸光係数が異なる２種類の波長、すなわち図４において波長λ１とλ２の光でそれぞれ照射するか、又は白色光を照射して、波長λ１とλ２の光のみ検知して、λ１とλ２それぞれの透過光の強度を測定し、その値を演算することにより、その物質（ａ）のみを強調した情報を得ることができる。このとき、物質（ｂ）や物質（ｃ）のように波長λ１とλ２の光の吸光係数に差が小さい物質が対象物体内部に含まれていても、２種の波長での透過光の強度差が小さいため、物質（ｂ）や物質（ｃ）の存在を強調させない画像を得ることができる。
さらに、本発明は非破壊で測定できるため、時間を追って特定の波長の透過光の強度を測定し、対象物体内部に含まれる一の物質の量の時間変化をモニターすることができる。このとき、時間変化を求めようとする物質の吸光度の差が最も大きい波長を選ぶことが好ましい。
本発明においては、対象物体を回転させながら透過光の強度を測定し、対象物体内部の３次元構造に関する情報を得ることもできる。
すなわち、対象物体の１断面をｎ×ｎの有限要素に分割し、それぞれの要素（ｉ，ｊ）の散乱定数をｆｄ（ｉ，ｊ）、吸収定数をｆａ（ｉ，ｊ）とすると、要素（ｉ，ｊ）をｊ方向に平行に透過する前の光強度Ｉｉｎ（ｉ，ｊ）と透過後の光強度Ｉｏｕｔ（ｉ，ｊ）の関係は、式（１）で示される（図５参照）。
0004575296.tif 000002
ここで、すべての要素において、式（２）が成り立つ。
0004575296.tif 000003
従って、対象物体に入射する光強度を、Ｉｉｎ、ｉに対応する透過光強度を、Ｉｏｕｔ（ｉ）とすると、式（１）は式（３）で表すことができる。
0004575296.tif 000004
さらに、各要素の透過定数をｆｔ（ｉ，ｊ）＝１−ｆｄ（ｉ，ｊ）ｆａ（ｉ，ｊ）とすると、
0004575296.tif 000005
この計算をすべてのｉについて行う。また、入射角度を変化させる場合は、入射光線の幅を変化させずに回転させ、１本１本の光線が交差するすべての要素について式（４）を計算する（図６参照）。その結果、式（５）のような連立方程式が算出される。
0004575296.tif 000006
0004575296.tif 000007
（ここで、Ａは、（ｊ・ｍ）×ｎ２の行列（ｍは、入射角のパターン数）、ＸおよびＢは、ｎ２次の縦ベクトルである）
ｊ・ｍがｎ２より大きくなるように条件を設定し、Ｍｏｏｒｅ−Ｐｅｎｒｏｓｅ型の一般逆行列を算出して、解Ｘを得る。
このような一連の作業を異なる波長に対して行い、それぞれのＸの値（Ｘλ＝ａ’Ｘλ＝ｂ）を算出する。このとき、吸収の波長依存性に比較して、散乱の波長依存性が充分小さいλを選択すると、（Ｘλ＝ａ’Ｘλ＝ｂ）から、各要素における、吸収に依存したパラメーターを算出し、対象物体内部の３次元吸収係数分布を可視化できる。
このようにして得られた対象物体内の各微小部分（各要素）の光透過係数の値から、対象物体の特定断面上の光透過係数を抽出することによって、対象物体を破壊することなくその断面の画像を得ることもできる。
本発明の対象物体内部の評価方法は、散乱光の影響が除かれた透過光の強度を測定できるため、分解能が高い画像情報を得ることができる。また、散乱光の影響が除かれた透過光を波長を変えて測定することによって、特定物質に着目した情報が得られたり、さらには非破壊であるため時間を追って測定することでリアルタイムの情報が得られる。また、対象物体を回転させて測定することによって、対象物体内部の３次元情報を得ることができる。
そして、本発明の対象物体内部の評価方法は特に微小部分の測定に適しており、毛髪内部に与えるカラーリング剤やヘアートリートメント剤の影響をリアルタイムで知ることができる。
【実施例】  以下に実施例を挙げて本発明を更に詳しく説明するが、本発明はこれら実施例に何ら制約されるものではない。
［実施例１］
＜毛髪の脱色過程の２次元動的解析＞
図２及び図３に示す対象物体セット用ホルダー（１０）を用い、対象物体として未処理の黒髪を用いた。対象物体セット用ホルダー（１０）は光学顕微鏡観察に汎用されるプレパラート（１１）の上に、固定用ガラス（１３ａ）、（１３ｂ）を用いて約１ｍｍの隙間を空けて固定された２枚の偏光子（１ｄ）、（１ｅ）から構成されており。２枚の偏光子間には、約２ｍｍ（幅）×約２５ｍｍ（奥行き）×約１ｍｍ（高さ）の空隙（１４）ができている。ここで、２枚の偏光子は互いに透過軸を平行とした。
対象物体セット用ホルダーに毛髪をセットした。撮像には、光学顕微鏡（オリンパス光学工業株式会社製、ＢＸ５１）を取り付けたＣＣＤカメラ（オリンパス光学工業株式会社製、ＤＰ１２）を用いた。対物レンズの倍率は１０倍とした。空隙（１４）にブリーチ溶液（４％Ｈ２Ｏ２，４％ＮＨ３水溶液）を満たし、セット直後、３分後、５分後、１０分後、２０分後に毛髪の同一箇所での撮像を行った。
得られたカラー画像から、毛髪の写っていない領域を取り除いた後（図７）、画像中の各チャネル（ＲＧＢ）の８ビット輝度分布（図８）を取得し、それぞれＲｔ（ｘ，ｙ）、Ｇｔ（ｘ，ｙ）、Ｂｔ（ｘ，ｙ）とした（（ｘ，ｙ）は画像中のピクセルの位置、ｔはブリーチ時間を示す）。その後、以下に示すアルゴリズムを用いて毛髪内部のメラニン色素の分布を示す画像（図９）を得た。
（アルゴリズム）
ｔ＝０の時、全ての（ｘ，ｙ）についてＤ０（ｘ，ｙ）を算出し（式（７））、その最大値をＤｍａｘとする。
0004575296.tif 000008
次に、全てのｔ、（ｘ，ｙ）についてＤｔを算出した後、８ビットに規格化したＭｔ（ｘ，ｙ）を得る（式（８））。
0004575296.tif 000009
黒色毛髪が時間と共に、脱色していく様子が明確に評価できた。この方法を用いると、カラーリング剤やヘアートリートメント剤が毛髪に与える影響をリアルタイムで知ることができる。
［実施例２］
＜毛髪の脱色過程の２次元動的解析＞
試料及び対象物体セット用ホルダーは実施例１と同様のものを用いた。毛髪を対象物体セット用ホルダーにセットしたのち、実施例１と同様の条件で撮像を行った。
得られたカラー画像から、毛髪の写っていない領域を取り除いた後、画像中の各チャネル（ＲＧＢ）の８ビット輝度分布を取得し、それぞれＲｔ（ｘ，ｙ）、Ｇｔ（ｘ，ｙ）、Ｂｔ（ｘ，ｙ）とした（（ｘ，ｙ）は画像中のピクセルの位置、ｔはブリーチ時間を示す）。その後、以下に示すアルゴリズムを用いて、毛髪内部のメラニン色素のうち、脱色反応によって吸光特性が変化したものの分布を示す画像（図１０）を得た。
（アルゴリズム）
ｔ＝２０の時、全ての（ｘ，ｙ）についてＤｔ（ｘ，ｙ）を算出し（式（９））、その最大値をＤｍａｘとする。Ｄｔ（ｘ，ｙ）の値が０より小さい時は、Ｄｔ（ｘ，ｙ）の値は０とする。
0004575296.tif 000010
次に、全てのｔ、（ｘ，ｙ）についてＤｔを算出した後、８ビットに規格化したＭｔ（ｘ，ｙ）を得る（式（１０））。Ｍｔ（ｘ，ｙ）が０よりも小さい時は０とする。
0004575296.tif 000011
毛髪が時間と共に、変化していく様子が明確に評価できた。この方法を用いると、カラーリング剤やヘアートリートメント剤が毛髪に与える影響をリアルタイムで知ることができる。
［実施例３］
＜毛髪中のメラニン色素の３次元分布の可視化＞
図２及び図３に示す対象物体セット用ホルダー（１０）を用い、光強度測定装置も実施例１と同様のものを用いた。対象物体は毛髪である。また、試料毛髪を回転させるために図１１に示す毛髪回転装置（２０）を用いた。試料毛髪（２１）を対象物体セット用ホルダー（１０）にセットしたのち、空隙（１４）に水を満たし、毛髪の両端を回転部（２２）の中心に固定した後、毛髪回転装置（２０）を用いて毛髪を回転させながら撮像を行った。回転装置は回転部受（２３）の部分を顕微鏡のステージ上に固定した。
得られたカラー画像から、毛髪の写っていない領域を取り除いた後、画像中の各チャネル（ＲＧＢ）の８ビット輝度分布を取得し、それぞれＲ（ｘ，ｙ，θ）、Ｇ（ｘ，ｙ，θ）、Ｂ（ｘ，ｙ，θ）とする（（ｘ，ｙ）は画像中のピクセルの位置、θは毛髪の回転角度（０≦θ≦２π）を示す）。その後、以下に示すアルゴリズムを用いて毛髪内部のメラニン色素の３次元分布を示す画像を得ることができた。
（アルゴリズム）
全ての（ｘ，ｙ）についてＤ（ｘ，ｙ）を算出し（式（１１））、その最大値をＤｍａｘとする。
0004575296.tif 000012
次に、Ｄを８ビットに規格化したＭ（ｘ，ｙ）を得る（式（１２））。
0004575296.tif 000013
その後、毛髪を含む３次元領域を適当な解像度、例えば断面５０×５０ピクセル、長さ２００ピクセルの解像度を持つボクセルで分割し、そのボクセルの集合をＱとする。Ｑ内のボクセルの位置を（ｘ，ｙ，ｚ）、Ｑ内の各ボクセルの光透過係数をＴ（ｘ，ｙ，ｚ）で表す（座標系は図１２参照）。各ボクセル内は均質であると仮定する。このとき、Ｍ（ｘ，ｙ，θ）の対数は、画像中の点（ｘ，ｙ）に結像した光が、試料内を透過する際に通過したＱ内の全ての点におけるＴ（ｘ，ｙ，ｚ）の積分値とみなせる。例えばθ＝０の時は、式（１３）のように表せる。
0004575296.tif 000014
同様に全てのｘ，ｙ，θにおいて計算し、得られたＴの連立方程式を適当な手法（例えば、ムーアーペンローズ（Ｍｏｏｒｅ−Ｐｅｎｒｏｓｅ）の一般逆行列を求める手法）を用いて解くことで、Ｑ内のＴの値を推定する。また、ｙ＝ｙ０とすればＱを平面ｙ＝ｙ０で切断した時の断面のＴを得ることができる。
【産業上の利用可能性】  本発明によれば、物体内部で散乱する光を除き、直線的に透過する光の強度を測定して画像化できるので、光を利用して微小部分の情報を得たり、特定物質を強調した画像情報を得たり、対象物体を回転させながら透過光の強度を測定し、対象物体を破壊することなしに、その内部の３次元構造を知ることができる。
また、非破壊であるため対象物体の変化をリアルタイムで知ることができる。さらに、対象物体を毛髪にすることで、実際に毛髪内部で起こっているカラーリング剤の脱色・染色状況や、ヘアートリートメントによる補修プロセスを明らかにでき、かかる製品を開発する上で広く利用できる。
【図１】
0004575296.tif 000015
【図２】
0004575296.tif 000016
【図３】
0004575296.tif 000017
【図４】
0004575296.tif 000018
【図５】
0004575296.tif 000019
【図６】
0004575296.tif 000020
【図７】
0004575296.tif 000021
【図８】
0004575296.tif 000022
【図９】
0004575296.tif 000023
【図１０】
0004575296.tif 000024
【図１１】
0004575296.tif 000025
【図１２】
0004575296.tif 000026