(57)【特許請求の範囲】
【請求項１】
セルロース系繊維やタンパク質系繊維など同系に分類される同系異種繊維を鑑別する繊維鑑別方法であって、
鑑別しようとする２種類（２グループ）以上の同系異種繊維に対して、それぞれ繊維の種類が既知の単一繊維を比較繊維として複数準備し、各比較繊維に対して赤外線又は近赤外線を照射してそれぞれの吸収スペクトルを求め、
これらの吸収スペクトルから得られたスペクトルデータＸを用いて、式（１）によりグループ間は離れ且つ各グループは纏まる軸ｗを求める判別分析を行い、得られたスコアプロットから判別モデルを作成しておき、
ＳＢｗ＝λ（ＳＷ＋ζＩ）ｗ  ・・・（１）
（ここで、ＳＢ・ＳＷは、分散共分散行列又は変動行列であって、分散共分散行列を使用する場合には、ＳＢはグループ間分散共分散行列、ＳＷはグループ内分散共分散行列であり、変動行列を使用する場合には、ＳＢはグループ間変動行列、ＳＷはグループ内変動行列である。一方、ζは正則化係数、Ｉは単位行列である。）
次に、繊維の種類が未知の繊維を被検繊維とし、前記判別モデルの作成と同様にして当該被検繊維のスペクトルデータＹから求めたスコアを前記判別モデルに当て嵌めて、前記被検繊維がいずれのグループに属するかを照合して、繊維の種類を鑑別する繊維鑑別方法。
【請求項２】
前記判別モデルの作成段階において、各比較繊維のスペクトルデータＸから求めたスコアプロットの等確率楕円を作成し、
前記被検繊維の鑑別段階において、当該被検繊維のスペクトルデータＹから求めたスコアを前記等確率楕円に当て嵌めて、前記被検繊維がいずれのグループに属するかを照合することを特徴とする請求項１に記載の繊維鑑別方法。
【請求項３】
前記判別モデルの作成段階において、各比較繊維のスペクトルデータＸから求めたスコアからｎ次元（ｎは１以上の整数）の正規分布の推定を行い、
前記被検繊維の鑑別段階において、当該被検繊維のスペクトルデータＹから求めたスコアから前記ｎ次元の正規分布に対する確率密度を算出することにより、前記被検繊維がいずれのグループに属するかを照合することを特徴とする請求項１に記載の繊維鑑別方法。
【請求項４】
前記比較繊維及び被検繊維に対して、近赤外線を除く波数４０００ｃｍ−１〜６００ｃｍ−１の範囲内の赤外線を照射して前記吸収スペクトルを求め、
これらの吸収スペクトルから所定の波数域における吸収スペクトルを抽出して前記スペクトルデータＸ及びＹを求めることを特徴とする請求項１〜３のいずれか１つに記載の繊維鑑別方法。
【請求項５】
前記判別モデルの作成にあたり、
第１回判別分析において各比較繊維の元のスペクトルデータをデータＸ１として求めた第１軸ｗ１とデータＸ１とから、式（２）によりデータＸ１の第１軸座標（ｗ１軸座標）の値ｔ１を算出し、
ｔ１＝ｗ１ＴＸ１  ・・・（２）
（ここで、ｗ１Ｔはｗ１の転置ベクトル）
式（３）により前記データＸ１からｗ１軸方向の情報を引き抜いたデータＸ２を求め、
Ｘ２＝Ｘ１−ｗ１ｔ１  ・・・（３）
次に、第２回判別分析においてスペクトルデータをデータＸ２として求めた第２軸ｗ２とデータＸ２とから、式（４）によりデータＸ２の第２軸座標（ｗ２軸座標）の値ｔ２を算出して、
ｔ２＝ｗ２ＴＸ２  ・・・（４）
（ここで、ｗ２Ｔはｗ２の転置ベクトル）
直交分解のなされた判別モデルを作成することを特徴とする請求項１〜４のいずれか１つに記載の繊維鑑別方法。
【請求項６】
上記式（３）及び式（４）の操作を複数回繰り返すことにより３次元以上の互いに直交する軸座標を有する判別モデルを作成することを特徴とする請求項５に記載の繊維鑑別方法。
【請求項７】
セルロース系繊維において、下記の各組み合わせに係る２種類（２グループ）の比較繊維、
（１）天然繊維、対、再生繊維、
（２）綿、対、麻類、
（３）亜麻、対、苧麻、
（４）ビスコース系レーヨン、対、銅アンモニアレーヨン又は溶剤紡糸セルロース繊維、
（５）銅アンモニアレーヨン、対、溶剤紡糸セルロース繊維、
のスペクトルデータＸを判別分析して得られた各判別モデルを使用して、前記被検繊維のスペクトルデータＹを前記各判別モデルに当て嵌めて、前記被検繊維がいずれのグループに属するかを照合して、繊維の種類を鑑別することを特徴とする請求項１〜６のいずれか１つに記載の繊維鑑別方法。
【請求項８】
セルロース系繊維において、天然繊維と再生繊維との鑑別には主に波数１３００〜８５０ｃｍ−１の範囲内若しくはその近傍の範囲内を含む１組又は２組以上のスペクトルデータを使用して鑑別することを特徴とする請求項１〜７のいずれか１つに記載の繊維鑑別方法。
【請求項９】
セルロース系繊維において、綿と麻類との鑑別には主に波数１６００〜８００ｃｍ−１の範囲内若しくはその近傍の範囲内を含む１組又は２組以上のスペクトルデータを使用して鑑別することを特徴とする請求項１〜７のいずれか１つに記載の繊維鑑別方法。
【請求項１０】
セルロース系繊維において、再生繊維どうしの鑑別には主に波数１６００〜９００ｃｍ−１の範囲内若しくはその近傍の範囲内を含む１組又は２組以上のスペクトルデータを使用して鑑別することを特徴とする請求項１〜７のいずれか１つに記載の繊維鑑別方法。
【請求項１１】
セルロース系繊維において、銅アンモニアレーヨンと溶剤紡糸セルロース繊維との鑑別には主に波数１５００〜８００ｃｍ−１の範囲内若しくはその近傍の範囲内を含む１組又は２組以上のスペクトルデータを使用して鑑別することを特徴とする請求項１〜７のいずれか１つに記載の繊維鑑別方法。
【請求項１２】
セルロース系繊維において、亜麻と苧麻との鑑別には主に波数１７００〜９００ｃｍ−１の範囲内若しくはその近傍の範囲内を含む１組又は２組以上のスペクトルデータを使用して鑑別することを特徴とする請求項１〜７のいずれか１つに記載の繊維鑑別方法。
【請求項１３】
前記式（１）における正則化係数ζの値を、１〜０の範囲内として鑑別することを特徴とする請求項１〜１２のいずれか１つに記載の繊維鑑別方法。
【請求項１４】
セルロース系繊維において、比較繊維及び被検繊維に対してアルカリ性物質による前処理を施してから吸収スペクトルを求めることを特徴とする請求項１〜１３のいずれか１つに記載の繊維鑑別方法。
【請求項１５】
前記比較繊維及び前記被検繊維の吸収スペクトルを求める方法は、ＡＴＲ法（全反射測定法）であることを特徴とする請求項１〜１４のいずれか１つに記載の繊維鑑別方法。
【請求項１６】
前記セルロース系繊維に分類される同系異種繊維としては、綿、亜麻、苧麻、黄麻、大麻、ビスコースレーヨン、ハイウェットモジュラスレーヨン、ポリノジックレーヨン、銅アンモニアレーヨン、及び、溶剤紡糸セルロース繊維が含まれることを特徴とする請求項１〜１５のいずれか１つに記載の繊維鑑別方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、繊維製品或いは織編物などに使用されている繊維の種類を鑑別する繊維鑑別方法に関するものである。特に、セルロース系繊維やタンパク質系繊維など同系に分類される同系異種繊維を鑑別する繊維鑑別方法に関するものである。
【背景技術】
【０００２】
市場には多くの繊維製品が広い用途に流通している。また、繊維製品の生産地と消費地がグローバルに展開される今日においては、繊維製品の輸出入の際に取引の安全や信頼を確保するために、各国の繊維関係の検査機関で繊維鑑別が行われている。【０００３】
これらの検査機関、例えば、日本の検査機関においては、ＪＩＳ Ｌ １０３０‐１（繊維製品の混用率試験方法‐第１部：繊維鑑別）、及び、ＪＩＳ Ｌ １０３０‐２（繊維製品の混用率試験方法‐第２部：繊維混用率）に基づいて鑑別を行っている。【０００４】
例えば、ＪＩＳ Ｌ １０３０‐１（繊維製品の混用率試験方法‐第１部：繊維鑑別）における鑑別方法には、燃焼試験、繊維中の塩素の確認試験、繊維中の窒素の確認試験、顕微鏡試験、よう素‐よう化カリウム溶液による着色試験、キサントプロテイン反応試験、赤外吸収スペクトルの測定試験などがある。【０００５】
これらの試験法はそれぞれ有効なものであり、これらを組み合わせることにより多くの繊維が鑑別できる。しかし、化学的組成が同じ繊維（以下「同系異種繊維」という）、例えば、綿、各種麻、各種レーヨン、銅アンモニアレーヨン、溶剤紡糸セルロース繊維などのセルロース系繊維や、カシミヤ、ウール、ヤク、モヘア、アンゴラ、アルパカ、ビキューナ、キャメル、リャマなどの獣毛繊維は、上記各試験法のうち化学的試験法では同一繊維或いは類似繊維と鑑別される。従って、これらの化学的試験法では同系異種繊維を明確に区別することはできない。【０００６】
また、上記試験法には赤外吸収スペクトルの測定試験も含まれる。しかし、この測定試験は、専ら化学的組成が異なる繊維（以下「異系繊維」という）を鑑別するものであり、同系異種繊維を明確に区別することはできないとされている。【０００７】
そこで、これらのセルロース系繊維やタンパク質系繊維などの同系異種繊維の鑑別には、主にその外観的特徴の違いを指標とする顕微鏡試験が有効であり、広く行われている。この顕微鏡試験による繊維鑑別を行うには、検査員が光学顕微鏡を用いて目視により検査対象の繊維を標準写真見本と対比させて行っている。【０００８】
従って、これらの方法においては、各検査機関の検査員の経験とノウハウの違いによる鑑別結果のバラツキが生じるという問題があった。更に、高価なタンパク質系繊維、特に獣毛繊維などには手の込んだ偽装が行われていることがあり、上記方法のみでは正確な鑑別が行えないという問題があった。【０００９】
更に、セルロース系繊維の中でもテンセル（登録商標）、リヨセル（登録商標）などの溶剤紡糸セルロース繊維（以下、代表して「リヨセル」ともいう）の外観的特徴は、銅アンモニアレーヨン（以下「キュプラ」ともいう）とほぼ同じ円形断面をしている。また、タンパク質系繊維の中でも特に高級品とされるカシミヤを用いた繊維製品には、カシミヤと見分けが付きにくいヤクの毛を混合し、或いは、ウールのスケールを除去（「脱スケール」という）して混合するなど手の込んだ偽装が行われている。このような場合には、顕微鏡試験に経験を積んだ検査員でも正確な判断が困難である。【００１０】
これに対して、下記特許文献１においては、セルロース系の同系異種繊維である繊維素繊維（溶剤紡糸セルロース繊維に同じ）とキュプラの鑑別方法が提案されている。この鑑別方法は、両繊維が６１％以上の硫酸に浸漬した時の溶解状態を顕微鏡下で観察してセルロース系の同系異種繊維を鑑別するというものである。【００１１】
また、下記特許文献２においては、獣毛系の同系異種繊維に対する繊維鑑別方法および繊維鑑別装置が提案されている。この繊維鑑別方法は、一般の赤外線より波長の長い電磁波を利用したテラヘルツ分光法を用いて獣毛繊維の細胞構造（一次構造）や細胞の集合形式（高次構造）を解析して獣毛系の同系異種繊維を鑑別するというものである。【００１２】
更に、下記非特許文献１においては、本発明の発明者により、一般の赤外線より波長の短い近赤外分光法を用いた布地材質の判別が提案されている。この判別方法は、近赤外分光法にフーリエ変換によるスペクトルのゆらぎ解析の手法を応用したものである。【００１３】
また、下記非特許文献２においては、本発明の発明者により、近赤外分光法による混紡繊維布地の繊維鑑別と混用率測定が提案されている。この鑑別方法は、近赤外分光法とケモメトリックス（化学を示すｃｈｅｍｉｓｔｒｙと計量学を示すｍｅｔｒｉｃｓからなる造語）を用いた解析手法を応用したものである。【００１４】
下記特許文献１の鑑別方法は、繊維の溶解状態を顕微鏡下で観察するというものであり、この場合にも検査員の経験とノウハウの違いによる鑑別結果のバラツキが生じるという問題があった。また、繊維に染色や樹脂加工などが施されている場合には、溶解状態が変化して正確な鑑別が行えないという問題があった。【００１５】
また、下記特許文献２の繊維鑑別方法においては、鑑別対象である繊維をテラヘルツ電磁波の波長より十分に小さいサイズ（１０μｍ程度）にしなければ、入射テラヘルツ電磁波が散乱して検出器に入射するテラヘルツ電磁波の強度が減衰する。その為、粉砕に伴う温度上昇を防ぎながら凍結粉砕する方法が要求される。これらの操作は煩雑であり、また、粉砕により繊維の高次構造が破壊され、情報量が減少するという問題があった。【００１６】
一方、下記非特許文献１の判別方法においては、判別の可能性は示唆されるものの、未だ正確な鑑別をするには至っていない。また、下記非特許文献２の鑑別方法においては、主に異系繊維である綿‐ポリエステル混紡布地の混用率を求める可能性を示唆するものであり、同系異種繊維の鑑別や混用率の鑑別を提案するものではない。【００１７】
このような問題に対処して、本発明者らは、先に下記特許文献３の繊維鑑別方法を提案した。この繊維鑑別方法においては、対比する２種類の同系異種繊維（比較繊維）に対して波数４０００ｃｍ−１〜６００ｃｍ−１の範囲内の赤外線を照射して求めた吸収スペクトルデータを主成分分析してデータベースを蓄積する。次に、繊維の種類が未知の繊維（被検繊維）の吸収スペクトルを求め、データベースのデータ群と照合して、被検繊維の種類及び混用率を鑑別するというものである。
【先行技術文献】
【特許文献】
【００１８】
【特許文献１】特開平１０−３３２６８４号公報
【特許文献２】特開２０１１−２０３１３８号公報
【特許文献３】ＷＯ２０１６／０２７７９２号公報
【非特許文献】
【００１９】
【非特許文献１】分光研究、第５３巻、第４号、Ｐ.２４９−２５６（２００４）
【非特許文献２】分光研究、第５８巻、第６号、Ｐ.２６８−２７４（２００９）
【発明の概要】
【発明が解決しようとする課題】
【００２０】
上記特許文献３の繊維鑑別方法は、鑑別操作が比較的簡単で客観性を有し、検査員の経験やノウハウに頼ることなく同系異種繊維の鑑別をすることのできる良好な方法である。しかし、主成分分析による鑑別では、鑑別する繊維の組み合わせによって鑑別精度にバラツキが生じることがあった。特に、セルロース系繊維の鑑別において、レーヨンと（キュプラ及びリヨセル）との鑑別や、キュプラとリヨセルとの鑑別で高度な鑑別精度（例えば、正答率９０％以上）を保証することができない場合があった。【００２１】
そこで、本発明は、上記問題に対処して、鑑別操作が比較的簡単で客観性を有し、検査員の経験やノウハウに頼ることなく同系異種繊維の鑑別をすることができ、且つ、鑑別する繊維の組み合わせによって生じる鑑別精度のバラツキを極力小さくして、高度な鑑別精度を実現することのできる繊維鑑別方法を提供することを目的とする。
【課題を解決するための手段】
【００２２】
上記課題の解決にあたり、本発明者らは、鋭意研究の結果、吸収スペクトルの情報量が豊富な赤外分光法を採用し、得られた吸収スペクトルを判別分析するにあたり正則化項を採用することにより同系異種繊維を高度な鑑別精度で鑑別できることを見出し、本発明の完成に至った。【００２３】
即ち、本発明に係る繊維鑑別方法は、請求項１の記載によると、
セルロース系繊維やタンパク質系繊維など同系に分類される同系異種繊維を鑑別する繊維鑑別方法であって、
鑑別しようとする２種類（２グループ）以上の同系異種繊維に対して、それぞれ繊維の種類が既知の単一繊維を比較繊維として複数準備し、各比較繊維に対して赤外線又は近赤外線を照射してそれぞれの吸収スペクトルを求め、
これらの吸収スペクトルから得られたスペクトルデータＸを用いて、式（１）によりグループ間は離れ且つ各グループは纏まる軸ｗを求める判別分析を行い、得られたスコアプロットから判別モデルを作成しておき、
ＳＢｗ＝λ（ＳＷ＋ζＩ）ｗ  ・・・（１）
（ここで、ＳＢ・ＳＷは、分散共分散行列又は変動行列であって、分散共分散行列を使用する場合には、ＳＢはグループ間分散共分散行列、ＳＷはグループ内分散共分散行列であり、変動行列を使用する場合には、ＳＢはグループ間変動行列、ＳＷはグループ内変動行列である。一方、ζは正則化係数、Ｉは単位行列である。）
次に、繊維の種類が未知の繊維を被検繊維とし、前記判別モデルの作成と同様にして当該被検繊維のスペクトルデータＹから求めたスコアを前記判別モデルに当て嵌めて、前記被検繊維がいずれのグループに属するかを照合して、繊維の種類を鑑別する。【００２４】
また、本発明は、請求項２の記載によると、請求項１に記載の繊維鑑別方法であって、
前記判別モデルの作成段階において、各比較繊維のスペクトルデータＸから求めたスコアプロットの等確率楕円を作成し、
前記被検繊維の鑑別段階において、当該被検繊維のスペクトルデータＹから求めたスコアを前記等確率楕円に当て嵌めて、前記被検繊維がいずれのグループに属するかを照合することを特徴とする。【００２５】
また、本発明は、請求項３の記載によると、請求項１に記載の繊維鑑別方法であって、
前記判別モデルの作成段階において、各比較繊維のスペクトルデータＸから求めたスコアからｎ次元（ｎは１以上の整数）の正規分布の推定を行い、
前記被検繊維の鑑別段階において、当該被検繊維のスペクトルデータＹから求めたスコアから前記ｎ次元の正規分布に対する確率密度を算出することにより、前記被検繊維がいずれのグループに属するかを照合することを特徴とする。【００２６】
また、本発明は、請求項４の記載によると、請求項１〜３のいずれか１つに記載の繊維鑑別方法であって、
前記比較繊維及び被検繊維に対して、近赤外線を除く波数４０００ｃｍ−１〜６００ｃｍ−１の範囲内の赤外線を照射して前記吸収スペクトルを求め、
これらの吸収スペクトルから所定の波数域における吸収スペクトルを抽出して前記スペクトルデータＸ及びＹを求めることを特徴とする。【００２７】
また、本発明は、請求項５の記載によると、請求項１〜４のいずれか１つに記載の繊維鑑別方法であって、
前記判別モデルの作成にあたり、
第１回判別分析において各比較繊維の元のスペクトルデータをデータＸ１として求めた第１軸ｗ１とデータＸ１とから、式（２）によりデータＸ１の第１軸座標（ｗ１軸座標）の値ｔ１を算出し、
ｔ１＝ｗ１ＴＸ１  ・・・（２）
（ここで、ｗ１Ｔはｗ１の転置ベクトル）
式（３）により前記データＸ１からｗ１軸方向の情報を引き抜いたデータＸ２を求め、
Ｘ２＝Ｘ１−ｗ１ｔ１  ・・・（３）
次に、第２回判別分析においてスペクトルデータをデータＸ２として求めた第２軸ｗ２とデータＸ２とから、式（４）によりデータＸ２の第２軸座標（ｗ２軸座標）の値ｔ２を算出して、
ｔ２＝ｗ２ＴＸ２  ・・・（４）
（ここで、ｗ２Ｔはｗ２の転置ベクトル）
直交分解のなされた判別モデルを作成することを特徴とする。【００２８】
また、本発明は、請求項６の記載によると、請求項５に記載の繊維鑑別方法であって、
上記式（３）及び式（４）の操作を複数回繰り返すことにより３次元以上の互いに直交する軸座標を有する判別モデルを作成することを特徴とする。【００２９】
また、本発明は、請求項７の記載によると、請求項１〜６のいずれか１つに記載の繊維鑑別方法であって、
セルロース系繊維において、下記の各組み合わせに係る２種類（２グループ）の比較繊維、
（１）天然繊維、対、再生繊維、
（２）綿、対、麻類、
（３）亜麻、対、苧麻、
（４）ビスコース系レーヨン、対、銅アンモニアレーヨン又は溶剤紡糸セルロース繊維、
（５）銅アンモニアレーヨン、対、溶剤紡糸セルロース繊維、
のスペクトルデータＸを判別分析して得られた各判別モデルを使用して、前記被検繊維のスペクトルデータＹを前記各判別モデルに当て嵌めて、前記被検繊維がいずれのグループに属するかを照合して、繊維の種類を鑑別することを特徴とする。【００３０】
また、本発明は、請求項８の記載によると、請求項１〜７のいずれか１つに記載の繊維鑑別方法であって、
セルロース系繊維において、天然繊維と再生繊維との鑑別には主に波数１３００〜８５０ｃｍ−１の範囲内若しくはその近傍の範囲内を含む１組又は２組以上のスペクトルデータを使用して鑑別することを特徴とする。【００３１】
また、本発明は、請求項９の記載によると、請求項１〜７のいずれか１つに記載の繊維鑑別方法であって、
セルロース系繊維において、綿と麻類との鑑別には主に波数１６００〜８００ｃｍ−１の範囲内若しくはその近傍の範囲内を含む１組又は２組以上のスペクトルデータを使用して鑑別することを特徴とする。【００３２】
また、本発明は、請求項１０の記載によると、請求項１〜７のいずれか１つに記載の繊維鑑別方法であって、
セルロース系繊維において、再生繊維どうしの鑑別には主に波数１６００〜９００ｃｍ−１の範囲内若しくはその近傍の範囲内を含む１組又は２組以上のスペクトルデータを使用して鑑別することを特徴とする。【００３３】
また、本発明は、請求項１１の記載によると、請求項１〜７のいずれか１つに記載の繊維鑑別方法であって、
セルロース系繊維において、銅アンモニアレーヨンと溶剤紡糸セルロース繊維との鑑別には主に波数１５００〜８００ｃｍ−１の範囲内若しくはその近傍の範囲内を含む１組又は２組以上のスペクトルデータを使用して鑑別することを特徴とする。【００３４】
また、本発明は、請求項１２の記載によると、請求項１〜７のいずれか１つに記載の繊維鑑別方法であって、
セルロース系繊維において、亜麻と苧麻との鑑別には主に波数１７００〜９００ｃｍ−１の範囲内若しくはその近傍の範囲内を含む１組又は２組以上のスペクトルデータを使用して鑑別することを特徴とする。【００３５】
また、本発明は、請求項１３の記載によると、請求項１〜１２のいずれか１つに記載の繊維鑑別方法であって、
前記式（１）における正則化係数ζの値を、１〜０の範囲内として鑑別することを特徴とする。【００３６】
また、本発明は、請求項１４の記載によると、請求項１〜１３のいずれか１つに記載の繊維鑑別方法であって、
セルロース系繊維において、比較繊維及び被検繊維に対してアルカリ性物質による前処理を施してから吸収スペクトルを求めることを特徴とする。【００３７】
また、本発明は、請求項１５の記載によると、請求項１〜１４のいずれか１つに記載の繊維鑑別方法であって、
前記比較繊維及び前記被検繊維の吸収スペクトルを求める方法は、ＡＴＲ法（全反射測定法）であることを特徴とする。【００３８】
また、本発明は、請求項１６の記載によると、請求項１〜１５のいずれか１つに記載の繊維鑑別方法であって、
前記セルロース系繊維に分類される同系異種繊維としては、綿、亜麻、苧麻、黄麻、大麻、ビスコースレーヨン、ハイウェットモジュラスレーヨン、ポリノジックレーヨン、銅アンモニアレーヨン、及び、溶剤紡糸セルロース繊維が含まれることを特徴とする。
【発明の効果】
【００３９】
上記構成によれば、本発明は、セルロース系繊維やタンパク質系繊維など同系に分類される同系異種繊維を鑑別する繊維鑑別方法である。化学的組成が同じ繊維であって起源を異にする繊維どうしを赤外吸収スペクトル（近赤外吸収スペクトルを含む）により鑑別することができる。上記構成によれば、まず、鑑別しようとする２種類（２グループ）以上の同系異種繊維に対して、それぞれ繊維の種類が既知の単一繊維を比較繊維として複数準備する。次に、これらの比較繊維に対して赤外線又は近赤外線を照射してそれぞれの吸収スペクトルを求める。次に、これらの吸収スペクトルから得られたスペクトルデータＸを基にして、上記式（１）により軸ｗを求める判別分析を行い、得られたスコアプロットから判別モデルを作成しておく。この軸ｗにより、グループ間は離れ且つ各グループは纏まることができるので各グループの繊維を明確に区別することができる。【００４０】
次に、繊維の種類が未知の繊維を被検繊維とし、判別モデルの作成と同様にして当該被検繊維のスペクトルデータＹから求めたスコアを判別モデルに当て嵌めて、被検繊維がいずれのグループに属するかを照合する。このことにより、被検繊維の種類を客観的に鑑別することができる。このように、比較繊維から得られた複数のスペクトルデータＸを判別分析し、得られた判別モデルを利用することにより、より正確な繊維鑑別を比較的簡単、且つ、客観的に行うことができる。【００４１】
従って、上記構成によれば、客観的な繊維鑑別をすることができ、被検繊維の種類を正確に鑑別することができる。これらの操作は比較的簡単であり、また、機器分析であることから、検査員の経験とノウハウの違いによる鑑別結果のバラツキが生じるということがない。【００４２】
また、上記構成によれば、判別モデルの作成段階において、各比較繊維のスペクトルデータＸから求めたスコアプロットの等確率楕円を作成するようにしてもよい。次に、被検繊維の鑑別段階において、当該被検繊維のスペクトルデータＹから求めたスコアを作成した等確率楕円に当て嵌めて、被検繊維がいずれのグループに属するかを照合する。このことにより、上記効果をより一層具体的に発揮することができる。【００４３】
また、上記構成によれば、判別モデルの作成段階において、各比較繊維のスペクトルデータＸから求めたスコアからｎ次元（ｎは１以上の整数）の正規分布の推定するようにしてもよい。次に、被検繊維の鑑別段階において、当該被検繊維のスペクトルデータＹから求めたスコアから推定したｎ次元の正規分布に対する確率密度を算出して、被検繊維がいずれのグループに属するかを照合する。このことにより、上記効果をより一層具体的に発揮することができる。【００４４】
また、上記構成によれば、比較繊維及び被検繊維に対して、近赤外線を除く波数４０００ｃｍ−１〜６００ｃｍ−１の範囲内の赤外線を照射するようにしてもよい。この範囲の赤外線で得られる吸収スペクトルからは豊富な情報が得られる。このことにより、上記効果をより一層具体的に発揮することができる。【００４５】
また、上記構成によれば、２回の判別分析を行って判別モデルを作成するようにしてもよい。具体的には、第１回判別分析において各比較繊維の元のスペクトルデータをデータＸ１として上記式（１）により第１軸ｗ１を求める。次に、この第１軸ｗ１とデータＸ１とから、上記式（２）によりデータＸ１の第１軸座標（ｗ１軸座標）の値ｔ１を算出する。更に、上記式（３）によりデータＸ１からｗ１軸方向の情報を引き抜いたデータＸ２を求める。【００４６】
次に、第２回判別分析においてスペクトルデータをデータＸ２として上記式（１）により第２軸ｗ２を求める。次に、この第２軸ｗ２とデータＸ２とから、式（４）によりデータＸ２の第２軸座標（ｗ２軸座標）の値ｔ２を算出する。このようにして、直交分解のなされた判別モデルを作成することができる。このことにより、上記効果をより一層具体的に発揮することができる。【００４７】
また、上記構成によれば、上記式（３）及び式（４）の操作を複数回繰り返すようにしてもよい。よって、３次元以上の互いに直交する軸座標を有する判別モデルを作成することができる。このことにより、上記効果をより一層具体的に発揮することができる。【００４８】
また、上記構成によれば、比較繊維及び被検繊維の判別分析を行う際に、鑑別対象とする同系異種繊維の組み合わせにより抽出する所定の波数域のスペクトルデータを変化させるようにしてもよい。このように、判別分析に用いる波数域を選択することで鑑別の精度を更に向上させることができる。【００４９】
また、上記構成によれば、上記式（１）における正則化係数ζの値を、１〜０の範囲内の範囲内として判別分析するようにしてもよい。このように、正則化係数ζの値を選択することで鑑別の精度を更に向上させることができる。【００５０】
また、上記構成によれば、比較繊維及び被検繊維に対してアルカリ性物質による前処理を施してから吸収スペクトルを求めるようにしてもよい。このように、前処理を施すことにより鑑別の精度を更に向上させることができる。【００５１】
よって、本発明によれば鑑別操作が比較的簡単で客観性を有し、検査員の経験やノウハウに頼ることなく同系異種繊維の鑑別をすることができ、且つ、鑑別する繊維の組み合わせによって生じる鑑別精度のバラツキを極力小さくして、高度な鑑別精度を実現することのできる繊維鑑別方法を提供することができる。
【図面の簡単な説明】
【００５２】
【図１】各種セルロース系繊維の吸収スペクトル（平均スペクトル）を示す図である。
【図２】図１の各スペクトルを１次微分したスペクトル（微分スペクトル）を示す図である。
【図３】第１実施形態において被検繊維を鑑別する解析手順を表す鑑別フロー図である。
【図４】図３の鑑別フロー図の一部を抽出した部分フロー図である。
【図５】実施例１で得られた「天然繊維」と「再生繊維」の判別モデル（ＦＤＯＤプロット）である。
【図６】実施例１で得られた「綿」と「麻類」の判別モデル（ＦＤＯＤプロット）である。
【図７】実施例１で得られた「リネン」と「ラミー」の判別モデル（ＦＤＯＤプロット）である。
【図８】実施例１で得られた「レーヨン」と「キュプラ及びリヨセル」の判別モデル（ＦＤＯＤプロット）である。
【図９】実施例１で得られた「キュプラ」と「リヨセル」の判別モデル（ＦＤＯＤプロット）である。
【図１０】第２実施形態において被検繊維を鑑別するイメージを表す概念図である。
【図１１】実施例２で得られた「天然繊維」と「再生繊維」のＦＤＯＤプロットである。
【図１２】実施例２で得られた「綿」と「麻類」のＦＤＯＤプロットである。
【図１３】実施例２で得られた「リネン」と「ラミー」のＦＤＯＤプロットである。
【図１４】実施例２で得られた「レーヨン」と「キュプラ及びリヨセル」のＦＤＯＤプロットである。
【図１５】実施例２で得られた「キュプラ」と「リヨセル」のＦＤＯＤプロットである。
【発明を実施するための形態】
【００５３】
本発明においては、化学的組成が同じ繊維を「同系異種繊維」として定義する。この同系異種繊維としては、以下の述べるセルロース系繊維及びタンパク質系繊維などが挙げられる。本発明において、セルロース系繊維とは、天然セルロース繊維や再生セルロース繊維などのセルロースからなる全ての繊維をいう。天然セルロース繊維としては、綿、及び、亜麻（リネン）、苧麻（ラミー）、黄麻（ジュート）、大麻（ヘンプ）などの麻類が挙げられる。再生セルロース繊維としては、ビスコースレーヨン、ハイウェットモジュラスレーヨン（「ＨＷＭレーヨン」ともいう）、ポリノジックレーヨン、銅アンモニアレーヨン（キュプラ）、溶剤紡糸セルロース繊維（テンセル及びリヨセル）などが挙げられる。以下、これらの再生セルロース繊維のうちビスコース法により再生される繊維であるビスコースレーヨン、ハイウェットモジュラスレーヨン（ＨＷＭレーヨン）、及び、ポリノジックレーヨンを「ビスコース系レーヨン」ともいう。
【００５４】
本発明において、タンパク系質繊維とは、羊毛や絹などの動物性タンパク質からなる動物繊維だけでなく、植物性タンパク質からなる繊維、及び、これらのタンパク質を含有する全ての繊維をいう。ここで、動物繊維には、羊毛以外の獣毛繊維（後述する）や、絹以外に昆虫やクモなどが作る繊維も含まれる。特に、クモ糸繊維は、工業的な生産が検討されており、例えば、クモ糸遺伝子をカイコのゲノムＤＮＡに組み込み、絹と同様にしてカイコにクモ糸のタンパク質を含有した糸を吐かせる研究がある（信州大学、中垣雅雄教授）。また、微生物培養で製造したクモ糸のタンパク質を紡糸してクモ糸を製造する企業が出現している（スパイバー株式会社）。
【００５５】
これらのタンパク系質繊維の中で、特に重要と考えられるのが獣毛繊維である。獣毛繊維とは、動物より得られる天然ケラチン質を主成分とする毛繊維の全てを含むものであり、ウール（羊の羊毛）、カシミヤ（カシミヤ山羊の毛）、ヤク（牛の一種ヤクの毛）、モヘア（アンゴラ山羊の毛）、アンゴラ（アンゴラ兎の毛）、アルパカ（小型こぶなしラクダのアルパカの毛）、ビキューナ（小型こぶなしラクダのビクーニャの毛）、キャメル（ラクダの毛）、リャマ（小型こぶなしラクダのリャマの毛）、フォックス（キツネの毛）、ミンク（イタチの一種ミンクの毛）、チンチラ（ネズミの一種チンチラの毛）、ラビット（ウサギの毛）などが挙げられる。
【００５６】
これらの同系異種繊維は、化学的組成が同じであり鑑別が容易ではない。特に、操作が簡単で客観的な鑑別法である赤外吸収スペクトルなどの光学的測定法は、化学的組成が異なる異系繊維の鑑別に効果を発揮するが、同系異種繊維の鑑別は難しいとされてきた。これに対して、本発明は同系異種繊維であるセルロース系繊維どうし、或いは、タンパク質系繊維どうしの鑑別などにおいて効果が発揮される。
【００５７】
例えば、セルロース系繊維においては、天然セルロース繊維と再生セルロース繊維との識別が可能となる。また、同じ天然セルロース繊維である綿と麻類との識別が可能となる。また、同じ麻類である亜麻（リネン）と苧麻（ラミー）との識別などが可能となる。更に、同じ再生セルロース繊維であるビスコース系レーヨンと「キュプラ及びリヨセル」との識別が可能となる。また、同じ再生セルロース繊維であり繊維断面形状が円形であり顕微鏡観察では判定し辛い、キュプラとリヨセルとの識別が可能となる。更に、タンパク質系繊維においては、高級なカシミヤとこれに酷似するヤク、脱スケールしたウールとの識別が可能となる。
【００５８】
ここで、本発明において判別分析（後述する）によって繊維鑑別をする際に、繊維の種類が既知の繊維（以下「比較繊維」という）、及び、鑑別対象である繊維（以下「被検繊維」という）から赤外吸収スペクトル（近赤外スペクトルを含む）を得る方法について説明する。この方法は、一般に赤外分光法（以下「ＩＲ分光法」という）又は近赤外分光法（以下「NＩＲ分光法」という）といわれ、測定対象の物質に赤外線（近赤外線）を照射し、透過光或いは反射光を分光することでスペクトルを得て対象物の特性を知る方法である。特に、ＩＲ分光法は、対象物の分子構造や状態を知るために使用され、化学的組成が異なる有機物の分析には極めて一般的な方法である。
【００５９】
一般に分光法は、ラジオ波からガンマ線まで広く電磁波の放出或いは吸収を測定する方法であって、その中でＩＲ分光法は、吸収スペクトルから得られる情報量が多く、産業において研究部門だけでなく製造部門、品質管理部門でも広く普及している。一方、NＩＲ分光法は、ＩＲ分光法に比べ得られる情報量は少ないが、近年の多変量解析によるケモメトリックスの発展により、各産業において鑑別への使用が提案されている。本発明においては、情報解析法として独自の工夫を加えた判別分析（後述する）を使用してＩＲ分光法及びNＩＲ分光法のいずれの分光法にも対応することができる。なお、本発明においても情報量の多さと鑑別精度とを考慮すると、ＩＲ分光法を使用することが好ましい。
【００６０】
一般に赤外線は、近赤外線、中赤外線、及び遠赤外線として区別されるが、更に細かく区別する場合もあり、その波長範囲の定義が明確でない。そこで、本発明においては、ＩＲ分光法で使用する波長範囲として、近赤外線を除く２５００ｎｍ〜１７０００ｎｍ（２５００ｎｍ〜２００００ｎｍとしてもよい）の範囲内とする。また、ＩＲ分光法では、波長よりも波数（本発明においては「ＷＮ」ともいう）を使用することが多い。本発明においては、ＩＲ分光法で使用する波数範囲として、近赤外線を除く４０００ｃｍ−１〜６００ｃｍ−１の範囲内とする。一方、本発明においては、NＩＲ分光法で使用する波長範囲として、８００ｎｍ〜２５００ｎｍの範囲内とする。また、NＩＲ分光法で使用する波数範囲として、１２５００ｃｍ−１〜４０００ｃｍ−１の範囲内とする。
【００６１】
なお、ＩＲ分光法とNＩＲ分光法とで使用する波数域の境界は、４０００ｃｍ−１とされている。従って、４０００ｃｍ−１自体が近赤外線であるのか、或いは、中赤外線であるのかで疑義が生じる場合がある。そこで、本発明においては、「近赤外線を除く４０００ｃｍ−１〜６００ｃｍ−１の範囲内」という意味を厳格に解釈する必要があるときには、「波数（ＷＮ）が４０００ｃｍ−１＞ＷＮ≧６００ｃｍ−１の範囲内」とすることができる。
【００６２】
上述のように、NＩＲ分光法は、繊維の鑑別への利用が提案されている（上記特許文献１及び２）。しかし、その情報量の少なさから同系異種繊維の実用的な鑑別には至っていない。そこで、本発明は、繊維から得られる情報量が多いにもかかわらず、これまで化学的組成が異なる異系繊維の鑑別のみに利用され、同系異種繊維の鑑別に利用されていない一般のＩＲ分光法を利用できるものである。また、本発明に使用する赤外分光光度計は、多くの産業で使用されていることから廉価であり、鑑別のための設備投資やメンテナンスの費用が抑えられる。なお、本発明は、独自の工夫を加えた判別分析（後述する）を使用することにより、情報量の少ないNＩＲ分光法を利用できるものである。
【００６３】
以下、ＩＲ分光法によって本発明を詳細に説明する。本発明においてＩＲ分光法で鑑別するには、一般に使用されているフーリエ変換赤外分光光度計（以下「ＦＴ／ＩＲ分光光度計」という）を使用することができる。比較繊維及び被検繊維の吸収スペクトルを得るには、繊維を微粉砕して臭化カリウム（ＫＢｒ）粉末と共に錠剤を形成して、その透過光を測定するＫＢｒ錠剤法などの透過法で測定してもよい。或いは、織編物を破壊することなくそのまま反射光を測定できるＡＴＲ法（全反射測定法）や、ＦＴ／ＩＲ分光光度計の普及と共に広く使われるようになった拡散反射法などの反射法で測定するようにしてもよい。また、得られた吸収スペクトルは、必要により通常の方法により大気補正、ベースライン補正、平滑化補正、潜り込み深さ補正などを行うことが好ましい。
【００６４】
なお、ＫＢｒ錠剤法において繊維を粉砕する場合には、ハサミでの切断、ミルによる粉砕、或いは、凍結粉砕などを利用することができる。また、粉砕後の繊維の長さ及び径は任意でよいが、微粉砕することが好ましい。また、ＫＢｒ錠剤法とＡＴＲ法のいずれを採用する場合であっても、繊維の吸湿状態を制御することが好ましい。特に、吸湿性の高いセルロース系繊維やタンパク質系繊維においては、同系異種繊維間で吸湿性が異なり吸収スペクトルが繊維中の自由水の影響を受けることがあるからである。
【００６５】
ここで、セルロース系繊維を例として、これに含まれる綿、麻類（リネン及びラミー）、ビスコース系レーヨン（以下、単に「レーヨン」という）、キュプラ、リヨセルの各繊維の吸収スペクトルを比較する。図１は、各種セルロース系繊維の吸収スペクトル（平均スペクトル）を示す図である。図１の各吸収スペクトル（１〜６）は、ＦＴ／ＩＲ分光光度計のＡＴＲ法で波数４０００ｃｍ−１〜６００ｃｍ−１の範囲内で測定した。なお、図１においては、標準サンプルとして結晶セルロース（６）であるアビセル（登録商標、旭化成工業）を標準サンプルとして各繊維（１〜５）と比較した。これらの吸収スペクトルは、標準偏差で規格化した平均スペクトルである。
【００６６】
図１において、３６００ｃｍ−１〜３１００ｃｍ−１にかけて、再生繊維であるレーヨン（３）、キュプラ（４）、リヨセル（５）の吸収スペクトルは、天然繊維である綿（１）、麻類（２）に比べ、広いピーク幅を示している。また、天然繊維である綿（１）、麻類（２）には、３３３０ｃｍ−１に鋭いピークが現れるが、再生繊維であるレーヨン（３）、キュプラ（４）、リヨセル（５）には現れない。更に、１２００ｃｍ−１〜８００ｃｍ−１にかけて現れるメインピークにおいて、天然繊維である綿（１）、麻類（２）には頂点に２つのピークが現れるが、再生繊維であるレーヨン（３）、キュプラ（４）、リヨセル（５）には１つのピークしか現れない。
【００６７】
このように、波数４０００ｃｍ−１〜６００ｃｍ−１の範囲内の吸収スペクトルを詳細に観察することにより、天然繊維と再生繊維とを鑑別できる可能性がある。しかし、これらの吸収スペクトルを比較するだけでは、セルロース系繊維において、天然繊維と再生繊維とを明確に鑑別することは容易ではない。また、同じ天然繊維である綿と麻類との鑑別、同じ麻類であるリネンとラミーとの鑑別、再生繊維どうしの鑑別、及び形状の類似したキュプラとリヨセルとの鑑別などは難しい。
【００６８】
本発明において、化学的組成が同じセルロース系繊維である天然繊維と再生繊維とがＩＲ分光法で得られた吸収スペクトルから鑑別できる理由は、これらの繊維を構成するセルロースの結晶状態が異なることによるものと思われる。一般にセルロース系繊維の結晶構造については、天然繊維がＩ型結晶であり、再生繊維がII型結晶を構成している。また、セルロース分子の配向状態について、天然繊維は分子の還元末端を同一方向に向けるパラレル構造を示し、再生繊維は分子の還元末端が交互に位置を変えるアンチパラレル構造を示している。
【００６９】
これらの結晶状態の違いにより、天然繊維と再生繊維の吸収スペクトルにおいては吸収位置のシフトや吸収強度の変化が生じるものと思われる。更に、本発明者らは、同じセルロース系天然繊維である綿と麻類において、或いは、同じセルロース系再生繊維であるレーヨン、キュプラ、リヨセルにおいても、結晶化度、結晶の大きさ、結晶の分布が異なることから、ＩＲ分光法の吸収スペクトルに何らかの情報が含まれているものと予測した。
【００７０】
以下、本発明に係る繊維鑑別方法について、各実施形態により詳細に説明する。なお、本発明は、下記の各実施形態にのみ限定されるものではない。本発明においては、各種類の比較繊維の吸収スペクトルから得られたスペクトルデータＸを判別分析し、得られた判別モデルとして準備する。この工程を各実施形態においては「判別モデル準備工程」という。次に、同様にして被検繊維の吸収スペクトルから得られたスペクトルデータＹを判別モデルと照合する。この照合において、被検繊維のスペクトルデータＹと比較繊維のスペクトルデータＸから得られた判別モデルとの一致性を指標として、被検繊維の種類を鑑別する。この工程を各実施形態においては「鑑別工程」という。
【００７１】
《第１実施形態》
本第１実施形態においては、複数種類の単一繊維に対して、各単一繊維間の鑑別を行うものである。例えば、セルロース系繊維においては、天然繊維と再生繊維との鑑別、綿と麻類との鑑別、麻類におけるリネンとラミーとの鑑別、再生繊維どうしの鑑別、及び、キュプラとリヨセルとの鑑別などについて説明する。なお、本第１実施形態においては、判別モデルのスコアの等確率楕円に対する被検繊維のスコアの適合度合から鑑別を行うものである。
【００７２】
（１）判別モデル準備工程
本第１実施形態においては、鑑別しようとする繊維の組み合わせ、例えば、綿と麻類との鑑別を行う場合には、これらの繊維を比較繊維として吸収スペクトルを求める。なお、これらの吸収スペクトルに対しては、所定の方法で各種補正を行うようにしてもよい。これらの補正としては、例えば、波数の変化による赤外光の潜り込み深さの補正、或いは、乗算的散乱補正（ＭＳＣ）などがある。
【００７３】
ここで、綿の吸収スペクトルを求める場合には、綿からなる複数の試料を１つのグループとしてそれぞれの吸収スペクトルを求める。また、麻類の吸収スペクトルを求める場合にも、麻類からなる複数の試料を１つのグループとしてそれぞれの吸収スペクトルを求める。ここで、麻類としてリネン、ラミー、ジュート、ヘンプなどを１つの麻類グループとして吸収スペクトルを求めるようにしてもよい。また、例えば、リネンとラミーとを別グループとして吸収スペクトルを求め、それぞれのグループを綿のグループと鑑別するようにしてもよい。
【００７４】
また、セルロース系繊維製品の場合には、織編物となった段階で染色性や物性を向上させるためにアルカリ性物質による繊維加工が施されることがある。この繊維加工としては、一般にアルカリ水溶液、又は、液体アンモニアなどによる浸漬処理が行われる。特に、綿に対する水酸化ナトリウム水溶液による浸漬処理を「マーセライズ加工」といい、広く行われている。なお、鑑別対象であるセルロース系繊維にアルカリ性物質による繊維加工が施されているか否かにより、赤外吸収スペクトルが若干変化することが考えられる。従って、アルカリ性物質による繊維加工の有無が混在した試料群による鑑別を行うと、鑑別精度が低下することも考えられる。
【００７５】
また、市場にはアルカリ性物質による繊維加工が施されたセルロース系繊維製品とアルカリ性物質による繊維加工が施されていないセルロース系繊維製品が流通する。そこで、本発明者らは、比較繊維及び被検繊維に対して所定濃度のアルカリ性物質による前処理をしてから赤外吸収スペクトルを求めることにより、鑑別精度が向上することを見出した。すなわち、アルカリ性物質による繊維加工の有無が混在した比較繊維及び被検繊維に対してアルカリ性物質による前処理を行うことにより、これらがアルカリ性物質による繊維加工が施されたセルロース系繊維として統一され、吸収スペクトルが近似して鑑別精度が向上するものと考える。
【００７６】
このアルカリ性物質による前処理の条件は、特に限定するものではなく、使用するアルカリ性物質の種類や濃度、処理する繊維の種類などにより適宜選定すればよい。例えば、セルロース系繊維が綿である場合には、通常のマーセライズ加工と同様にアルカリ性物質として８重量％〜２４重量％濃度、好ましくは１０重量％〜２０重量％濃度、更に好ましくは１５重量％〜２０重量％濃度の水酸化ナトリウム水溶液を使用し、室温にて浸漬処理することにより鑑別精度が向上する。    
【００７７】
なお、判別モデルを準備する際の比較繊維、及び、鑑別対象である被検繊維に対して上記のアルカリ性物質による前処理を施した赤外吸収スペクトルと施していない赤外吸収スペクトルの両方を求め、これらを併用して鑑別することにより鑑別精度は更に向上する。
【００７８】
更に、各繊維には、織編物の区別、糸の太さの違い、染色の有無、繊維加工剤の有無などがある。従って、これらの情報が吸収スペクトルに与える影響と、解析に使用する波数域（後述する）との関係を考慮して各グループの吸収スペクトルを求めるようにする。
【００７９】
本第１実施形態においては、比較繊維の吸収スペクトルを求める際にＦＴ／ＩＲ分光光度計を使用する。ＦＴ／ＩＲ分光光度計は、多くの産業で使用されていることから廉価であり、鑑別のための設備投資やメンテナンスの費用が抑えられるからである。また、本第１実施形態においては、測定にＡＴＲ法を使用する。鑑別しようとする繊維は、多くの場合、既に繊維製品として織編物の状態にあり、この織編物の状態のまま非破壊で鑑別できるという利点があるからである。なお、ＡＴＲ法に使用するプリズムの種類は、特に限定するものではないが、繊維試料の場合、一般にＺｎＳｅプリズムを使用することが好ましい。
【００８０】
本第１実施形態においては、近赤外線を除く波数４０００ｃｍ−１〜６００ｃｍ−１の範囲内の赤外線を照射してそれぞれの吸収スペクトルを求める。ここで、「近赤外線を除く」とは、吸収スペクトルを測定する際に「近赤外線部分の吸収スペクトルを測定しない」という意味に解するものではない。「鑑別に使用する吸収スペクトルの範囲が波数４０００ｃｍ−１〜６００ｃｍ−１の範囲内のものである」と解するものである。
【００８１】
また、本発明においては、波数４０００ｃｍ−１〜６００ｃｍ−１の範囲内の全ての領域の吸収スペクトルで解析を行うようにしてもよい。或いは、特定の波数域を選択して解析を行うようにしてもよい。本第１実施形態においては、所定の波数域におけるスペクトルデータを解析に使用する。解析に使用する波数域を限定することにより、鑑別に必要な赤外吸収が強調されると共にノイズを排除して鑑別精度が向上するからである。
【００８２】
ここで、所定の波数域におけるスペクトルデータの抽出には、予め測定領域を決めておき、その範囲内でのみ吸収スペクトルを求めるようにしてもよい。しかし、一般には使用する分光光度計の全測定領域で得られた吸収スペクトルから解析に必要とする所定の波数域のスペクトルデータを抽出することが好ましい。
【００８３】
なお、ここでいう所定の波数域は、鑑別する繊維の種類と組み合わせにより適宜選択することが好ましい。本発明者らは、セルロース系繊維において、天然繊維と再生繊維との鑑別、綿と麻類との鑑別、リネンとラミーとの鑑別、再生繊維どうしの鑑別、及び、キュプラとリヨセルとの鑑別などにおいて、それぞれに適切な波数域が存在することを確認した。なお、鑑別に使用される所定の波数域については、１つの波数域だけで解析するようにしてもよく、或いは、２又はそれ以上の複数の波数域を組み合わせて解析するようにしてもよい。
【００８４】
ここで、比較繊維のスペクトルデータを解析するに当たり、まず、スペクトルデータの処理をすることが好ましい。スペクトルデータの処理には、例えば、１次微分、２次微分、或いはそれ以上の高次微分処理などを行うことが好ましい。これらの微分処理により、スペクトルに埋もれたピークの先鋭化やバックグランドの影響の消去など、スペクトル情報の増強を行うことができるからである。微分処理で使用する方法と次元は、特に限定するものではなく、鑑別する繊維の種類と組み合わせにより適宜選択することが好ましい。一例として、図１の各スペクトル（１〜５）をSavitzky-Golay法により１次微分した各スペクトル（１〜５）を図２に示す。
【００８５】
次に、得られた微分スペクトルから解析に有効な所定の波数域のスペクトルデータを抽出する。スペクトルデータを抽出する方法は特に限定するものではない。鑑別する繊維の種類と組み合わせに対して、解析に必要と考えられる特定の官能基による波数域を選定するようにしてもよい。また、解析に対してノイズを有すると考えられる波数域のスペクトルデータを積極的に排除するようにしてもよい。例えば、セルロース系繊維の場合には、波数２７５０〜１８５０ｃｍ−１の範囲内のスペクトルデータを解析に使用しないようにしてもよい。この波数域からは解析に有効な情報が少なく、逆に、ＣＯ２の吸収などがノイズとして現れることが考えられるからである。
【００８６】
例えば、セルロース系繊維の鑑別においては、波数３５００ｃｍ−１〜３０００ｃｍ−１の範囲内（Ｏ-Ｈ）のスペクトルデータ、及び、波数１２００ｃｍ−１〜１０００ｃｍ−１の範囲内（Ｃ-ＯＨ、Ｃ-Ｏ-Ｃ、Ｃ-Ｃ）のスペクトルデータなどが重要と考えられるので、これらを組み合わせて解析するようにしてもよい。但し、吸湿性の高いセルロース系繊維の鑑別においては、波数３５００ｃｍ−１〜３０００ｃｍ−１の範囲内のスペクトルデータが吸湿状態の影響を受けやすく、これらを排除して鑑別することも考えられる。
【００８７】
本発明者らは、赤外吸収の大きな波数域の組み合わせを用いて判別分析を行い、特に有効と考えられる特定の波数域を抽出することができた。例えば、天然繊維と再生繊維との鑑別において、波数１３００〜８５０ｃｍ−１の範囲内（Ｃ-ＯＨ、Ｃ-Ｏ-Ｃ、Ｃ-Ｃ）のスペクトルデータを抽出した。また、綿と麻類との鑑別において、波数１６００〜８００ｃｍ−１の範囲内（Ｏ-Ｈ、Ｃ-ＯＨ、Ｃ-Ｏ-Ｃ、Ｃ-Ｃ）のスペクトルデータを抽出した。また、リネンとラミーとの鑑別において、波数１６００〜９００ｃｍ−１の範囲内（Ｏ-Ｈ、Ｃ-ＯＨ、Ｃ-Ｏ-Ｃ、Ｃ-Ｃ）のスペクトルデータを抽出した。また、再生繊維どうしの鑑別すなわちレーヨンと「キュプラ及びリヨセル」との鑑別において、波数１５００〜８００ｃｍ−１の範囲内（Ｏ-Ｈ、Ｃ-ＯＨ、Ｃ-Ｏ-Ｃ、Ｃ-Ｃ）のスペクトルデータを抽出した。また、キュプラとリヨセルとの鑑別において、波数１７００〜９００ｃｍ−１の範囲内（Ｏ-Ｈ、Ｃ-ＯＨ、Ｃ-Ｏ-Ｃ、Ｃ-Ｃ）のスペクトルデータを抽出した。
【００８８】
次に、抽出したスペクトルデータを判別分析して判別モデルを作成する。ここで、判別分析とは、事前に与えられているデータが異なるグループに分かれる場合、新しいデータが得られた際に、どちらのグループに入るのかを判別するための基準（判別関数）を得るための
1495603678060_0
の手法をいう。判別分析として使用される手法は、特に限定するものではなく、どのような手法を採用するようにしてもよい。例えば、フィッシャーの判別分析であるＦＤＡ（Fisher’s linear discriminant analysis）などが一般的である。これに対して、本発明者らは、判別分析の手法に独自の工夫を加えたＦＤＯＤ（Fisher’s linear discriminant orthogonal decomposition）を使用した。以下にＦＤＯＤによる判別分析を具体的に説明する。
【００８９】
ＦＤＯＤにおいては、２つのグループを分離する軸ｗ（グループ間は離れ且つ各グループは纏まる軸ｗ）を求めるにあたり、直交分解により互いに直交する２つの軸ｗ１、ｗ２を求めることを特徴とする。また、本第１実施形態においては、ＦＤＯＤにおいて正則化を行うことを特徴とする。具体的には、下記式（１）において、正則化項ζＩを導入する。
【００９０】
ＳＢｗ＝λ（ＳＷ＋ζＩ）ｗ  ・・・（１）
上記式（１）において、ＳＢ・ＳＷは、分散共分散行列又は変動行列である。ＳＢ・ＳＷに分散共分散行列を使用する場合には、ＳＢはグループ間分散共分散行列、ＳＷはグループ内分散共分散行列を表している。一方、ＳＢ・ＳＷに変動行列を使用する場合には、ＳＢはグループ間変動行列、ＳＷはグループ内変動行列を表している。また、ζは正則化係数、Ｉは単位行列を表している。なお、本第１実施形態においては、前記式（１）における正則化係数ζの値を、１〜０の範囲内とすることが好ましく、また、１×１０−３〜１×１０−８の範囲内とすることがより好ましい。
【００９１】
一般に、判別分析において正則化を行うことにより、モデルサンプル（この場合は比較繊維）の判別精度が低下するといわれている。このことは、一方で過学習を防止するという意味も考えられる。過学習を起こすと、モデルサンプル（この場合は比較繊維）に過度に適合した判別モデルができてしまい、テストサンプル（この場合は被検繊維）への適合性が悪くなり汎化性能の悪い状態となる。そこで、本発明者らは適切な正則化係数ζを選択することにより、高い判別精度と良好な汎化性能を両立できることを見出した。
【００９２】
本第１実施形態においては、上述の軸の直交分解を行うために２回の判別分析を行う。まず、各比較繊維の吸収スペクトルから上述の特定波数域のスペクトルデータＸを抽出し、これをデータＸ１とする。このデータＸ１を用いて上記式（１）により第１回判別分析を行って第１軸ｗ１を求める。このとき、予め選定した正則化係数ζ１を使用する。なお、判別分析に使用する解析ソフトについては、特に限定するものではない。本第１実施形態においては、発明者自らが構成したプログラムを使用して解析した。
【００９３】
次に、第１回判別分析で得られた第１軸ｗ１とデータＸ１とから、下記式（２）によりデータＸ１の第１軸座標（ｗ１軸座標）の値ｔ１（射影ｔ１ともいう）を算出する。
【００９４】
ｔ１＝ｗ１ＴＸ１  ・・・（２）
ここで、ｗ１Ｔはｗ１の転置ベクトルを表している。
【００９５】
次に、第２回判別分析に使用するデータを求めるためにデフレーションという操作を行う。具体的には、下記式（３）により先のデータＸ１からｗ１軸方向の情報を引き抜いたものをデータＸ２（残差行列ともいう）とする。
【００９６】
Ｘ２＝Ｘ１−ｗ１ｔ１  ・・・（３）
次に、このデータＸ２を用いて上記式（１）により第２回判別分析を行って第２軸ｗ２を求める。このとき、予め選定した正則化係数ζ２を使用する。この時使用する正則化係数ζ２は、第１回判別分析で使用した正則化係数ζ１と異なる値であってもよい。このようにして求めた第２軸ｗ２とデータＸ２とから、式（４）によりデータＸ２の第２軸座標（ｗ２軸座標）の値ｔ２（射影ｔ２ともいう）を算出する。
【００９７】
ｔ２＝ｗ２ＴＸ２  ・・・（４）
ここで、ｗ２Ｔはｗ２の転置ベクトルを表している。
【００９８】
このようにして求めた互いに直交する２つの軸ｗ１、ｗ２を用いて、直交分解のなされた判別モデルを作成することができる。この判別モデルにおいては、軸ｗ１、ｗ２を直交軸として２つのグループが明確に区別される。
【００９９】
各グループに属する各比較繊維の射影ｔ１及び射影ｔ２に対応する値をその比較繊維のスコアという。２つの軸ｗ１、ｗ２に対して、各比較繊維のスコアをプロットしたものをスコアプロット（ＦＤＯＤプロット）という。本第１実施形態においては、判別モデルとして得られた２つのスコアプロットについて、それぞれ等確率楕円を作成する。なお、これらの判別モデル及び等確率楕円については、下記の実施例１において詳述する。
【０１００】
また、本第1実施形態においては、上述のように互いに直交する２つの軸ｗ１、ｗ２を用いて、直交分解のなされた判別モデルを作成するものであるが、上記式（３）及び式（４）の操作を繰り返すことにより、互いに直交する３軸以上の高次元の判別モデルを作成するようにしてもよい。
【０１０１】
（２）鑑別工程
本第１実施形態の鑑別工程においては、まず、鑑別しようとする被検繊維の吸収スペクトルを求める。吸収スペクトルを求める方法、吸収スペクトルに各種補正をする方法、及び、得られた吸収スペクトルに微分処理などの処理を行う方法は、上述の比較繊維に対する方法と同様である。次に、求めた微分スペクトルから比較繊維と同じ波数域のスペクトルデータＹを抽出し、比較繊維と同様にして被検繊維のスコア（射影ｔ１及び射影ｔ２に対応）を算出する。次に、得られた被検繊維のスコアを判別モデルの等確率楕円に当て嵌めて、被検繊維がいずれのグループに属する繊維であるかを鑑別する。
【０１０２】
なお、本鑑別工程においては、被検繊維の種類は未知であるが、比較的簡単な顕微鏡法などで被検繊維がセルロース系繊維であることが判明している。しかし、被検繊維がセルロース系繊維のうちのいずれの繊維であるかが不明である。そこで、本第１実施形態においては、次のような手順で繊維鑑別することが好ましい。
【０１０３】
図３は、本第１実施形態において被検繊維を鑑別する解析手順を表す鑑別フロー図である。図３において、まず、被検繊維が天然繊維であるか再生繊維であるかを鑑別する。このときには、天然繊維と再生繊維とから得られた判別モデルを使用する。ここで、被検繊維が天然繊維であると鑑別された場合には、続いて被検繊維が綿であるか麻類であるかを鑑別する。このときには、綿と麻類とから得られた判別モデルを使用する。更に、被検繊維が麻類であると鑑別された場合には、続いて被検繊維がどのような麻であるかを鑑別する。なお、図３においては、麻類として最も一般的なリネンであるかラミーであるかを鑑別する。このときには、リネンとラミーとから得られた判別モデルを使用する。
【０１０４】
同様に、図３において、被検繊維が再生繊維であると鑑別された場合には、続いて被検繊維がレーヨンであるか「キュプラ又はリヨセル」であるかを鑑別する。このときには、レーヨンと「キュプラ及びリヨセル」とから得られた判別モデルを使用する。ここで、被検繊維が「キュプラ又はリヨセル」であると鑑別された場合には、続いて被検繊維がキュプラであるかリヨセルであるかを鑑別する。このときには、キュプラとリヨセルとから得られた判別モデルを使用する。
【０１０５】
また、図４は、図３の鑑別フロー図の一部を抽出した部分フロー図である。図４において、被検繊維が天然繊維であるか再生繊維であるかを鑑別した際に、天然繊維と再生繊維とから得られた判別モデルでは対応できないものがある。この場合、図４の判別不可の被検繊維は、天然繊維と再生繊維の混合繊維であるか、或いは、セルロース系繊維以外の繊維が混合されている可能性がある。
【０１０６】
このようにして、図３の鑑別フロー図に従って鑑別を行うことにより、被検繊維の種類を客観的に鑑別することができる。次に、本第１実施形態の鑑別方法を実施例１により具体的に説明する。
【実施例１】
【０１０７】
本実施例１は、上記第１実施形態に係る各単一繊維間の鑑別を行うものであり、複数の被検繊維に対して鑑別フロー図（図３参照）に従って鑑別を行うものである。なお、各被検繊維は、いずれも顕微鏡法などの予備鑑定においてセルロース系繊維であることが判明している。【０１０８】
（１）判別モデル準備工程
本実施例１においては、まず、繊維の種類が既知の単一繊維として、綿２７点、リネン２５点、ラミー２１点の天然繊維７３点を準備した。また、繊維の種類が既知の単一繊維として、レーヨン４８点、キュプラ３０点、リヨセル３１点の再生繊維１０９点を準備した。これら合計１８２点の単一繊維を本実施例１の比較繊維とした。なお、綿及び麻類（リネン及びラミー）７３点、並びに、再生繊維の一部４０点（レーヨン１４点、キュプラ１３点、リヨセル１３点）に関しては、室温において水酸化ナトリウム水溶液（１７重量％）による前処理を行った。【０１０９】
次に、これら合計１８２点の織編物について、その吸収スペクトルを得た。吸収スペクトルの測定は、赤外分光光度計ＦＴ／ＩＲ‐４７００（日本分光株式会社）を使用し、ＺｎＳｅプリズムによるＡＴＲ法で、波数４０００ｃｍ−１〜６００ｃｍ−１の吸収スペクトルを測定した。次に、乗算的散乱補正（ＭＳＣ）後の各吸収スペクトルをSavitzky-Golay法により１次微分して微分スペクトルを得た。【０１１０】
次に、これらの吸収スペクトル（微分スペクトル）をＦＤＯＤによる判別モデル作成用グループと等確率楕円作成用グループに分類した。吸収スペクトルを判別モデル作成用グループと等確率楕円作成用グループに分類する理由については後述する。【０１１１】
［判別モデル作成用グループ］
・第１グループ（Ａ１）：３０点の天然繊維（綿１５点、リネン７点及びラミー８点）
・第２グループ（Ａ２）：３０点の再生繊維（前処理レーヨン１０点、前処理キュプラ１０点及び前処理リヨセル１０点）
・第３グループ（Ｂ１）：１５点の綿
・第４グループ（Ｂ２）：１５点の麻類（リネン７点及びラミー８点）
・第５グループ（Ｃ１）：１５点のリネン
・第６グループ（Ｃ２）：１５点のラミー
・第７グループ（Ｄ１）：３５点の未処理レーヨン   
・第８グループ（Ｄ２）：３５点のキュプラ及びリヨセル（未処理キュプラ１８点及び未処理リヨセル１７点）
・第９グループ（Ｅ１）：１５点の未処理キュプラ
・第１０グループ（Ｅ２）：１５点の未処理リヨセル
このような上記１０グループを得た。【０１１２】
［等確率楕円作成用グループ］
・第１１グループ（Ａ３）：１０点の天然繊維（綿５点、リネン２点及びラミー３点）
・第１２グループ（Ａ４）：１０点の再生繊維（前処理レーヨン４点、前処理キュプラ３点及び前処理リヨセル３点）
・第１３グループ（Ｂ３）：１０点の綿
・第１４グループ（Ｂ４）：１０点の麻類（リネン７点及びラミー３点）
・第１５グループ（Ｃ３）：５点のリネン
・第１６グループ（Ｃ４）：５点のラミー
・第１７グループ（Ｄ３）：１０点の未処理レーヨン 
・第１８グループ（Ｄ４）：１０点のキュプラ及びリヨセル（未処理キュプラ５点及び未処理リヨセル５点）
・第１９グループ（Ｅ３）：１０点の未処理キュプラ
・第２０グループ（Ｅ４）：１０点の未処理リヨセル
このような上記１０グループを得た。【０１１３】
次に、判別モデル作成用グループを用いて判別モデルを作成した。上記第１グループ（Ａ１）〜第１０グループ（Ｅ２）のうち２つずつ組み合わせて、５組の組み合わせとした。次に、これらの微分スペクトルから各組み合わせを特徴付ける所定の波数域を抽出した。更に、各組み合わせに対する判別分析用の２つの正則化係数ζ１、ζ２を特定した。本実施例１において、各判別モデルＡ〜Ｅと抽出した波数域及び２つの正則化係数ζ１、ζ２の値を下記に示す。【０１１４】
［判別モデルＡ］：天然繊維と再生繊維
・第１グループ（Ａ１）と第２グループ（Ａ２）：波数１２００〜８５０ｃｍ−１
・第１回判別分析用の正則化係数：ζ１＝１×１０−５
・第２回判別分析用の正則化係数：ζ２＝１×１０−５
［判別モデルＢ］：綿と麻類
・第３グループ（Ｂ１）と第４グループ（Ｂ２）：波数１４００〜９００ｃｍ−１
・第１回判別分析用の正則化係数：ζ１＝１×１０−７
・第２回判別分析用の正則化係数：ζ２＝１×１０−７
［判別モデルＣ］：リネンとラミー
・第５グループ（Ｃ１）と第６グループ（Ｃ２）：波数１４００〜９００ｃｍ−１
・第１回判別分析用の正則化係数：ζ１＝１×１０−７
・第２回判別分析用の正則化係数：ζ２＝１×１０−７
［判別モデルＤ］：レーヨンと（キュプラ及びリヨセル）
・第７グループ（Ｄ１）と第８グループ（Ｄ２）：波数１４００〜９００ｃｍ−１
・第１回判別分析用の正則化係数：ζ１＝１×１０−６
・第２回判別分析用の正則化係数：ζ２＝１×１０−６
［判別モデルＥ］：キュプラとリヨセル
・第９グループ（Ｅ１）と第１０グループ（Ｅ２）：波数１４００〜９００ｃｍ−１ 
・第１回判別分析用の正則化係数：ζ１＝１×１０−６
・第２回判別分析用の正則化係数：ζ２＝１×１０−６
このようにして抽出した波数域における各グループのスペクトルデータＸを第１回目の判別分析に使用するデータＸ１とした。このデータＸ１を使用して５組の組み合わせに対する２回の判別分析を行い、それぞれの比較繊維のスコア（射影ｔ１及び射影ｔ２に対応）を算出した。各グループの比較繊維のスコアを互いに直交する２つの軸ｗ１、ｗ２にプロットして判別モデル（仮モデル）Ａ〜Ｅを得た。【０１１５】
次に、判別モデル（仮モデル）Ａ〜Ｅの各グループに対して、等確率楕円作成用グループを用いて等確率楕円を作成した。本実施例１において、判別モデル作成用の比較繊維とは別に等確率楕円作成用の比較繊維を使用する。その理由は、次のように考えられる。すなわち、判別分析において判別モデルを作成した場合には、グループ間は離れ且つ各グループは纏まる軸ｗが求められる。従って、この狭い領域のデータのみで等確率楕円を作成した場合には、被検繊維のデータ上の僅かの違いで等確率楕円への当てはまりが悪くなる場合が生じる。【０１１６】
これを解消するには、２つの方法が考えられる。１つは、各グループの比較繊維の数（サンプル数）を非常に多くして判別モデルを作成し、これらから作成した等確率楕円からの当てはまりを良くする方法である。いま１つは、限られた数の比較繊維で判別モデルを作成し、判別モデル作成用の比較繊維とは別に等確率楕円作成用の比較繊維を使用して当てはまりの悪い情報を考慮した等確率楕円を作成する方法である。本実施例１においては、比較繊維の数が限られることから、後者の方法を採用した。【０１１７】
本実施例１においては、９５％信頼水準の等確率楕円を作成した。具体的には、上記第１１グループ（Ａ３）〜第２０グループ（Ｅ４）の等確率楕円用の各比較繊維に対して、上記と同じ手法によりスコア（射影ｔ１及び射影ｔ２に対応）を算出した。これらのスコアをそれぞれのグループがプロットされた各判別モデルＡ〜Ｅにプロットした。次に、発明者自らが構成したプログラムを使用して９５％信頼水準の等確率楕円を作成した。【０１１８】
図５〜図９は、本実施例１で得られた各判別モデル（完成モデル）Ａ〜Ｅ（ＦＤＯＤプロット）である。また、各判別モデルの記載されたグループを囲む楕円は、各グループに対する９５％信頼水準の等確率楕円である。図５は、判別モデルＡ：「天然繊維」と「再生繊維」のＦＤＯＤプロットであり、第１判別軸（ｗ１軸）を横軸とし第２判別軸（ｗ２軸）を縦軸として、「天然繊維」と「再生繊維」の２つのグループが明確に層別された。図６は、判別モデルＢ：「綿」と「麻類」の散布図であり、第１判別軸を横軸とし第２判別軸を縦軸として、「綿」と「麻類」の２つのグループが明確に層別された。図７は、判別モデルＣ：「リネン」と「ラミー」の散布図であり、第１判別軸を横軸とし第２判別軸を縦軸として、「リネン」と「ラミー」の２つのグループが明確に層別された。図８は、判別モデルＤ：「レーヨン」と「キュプラ及びリヨセル」の散布図であり、第１判別軸を横軸とし第２判別軸を縦軸として、「レーヨン」と「キュプラ及びリヨセル」の２つのグループが明確に層別された。図９は、判別モデルＥ：「キュプラ」と「リヨセル」の散布図であり、第１判別軸を横軸とし第２判別軸を縦軸として、「キュプラ」と「リヨセル」の２つのグループが明確に層別された。このようにして求めた解析データ群を各組み合わせに対する判別モデルとして本実施例１のデータベースとして蓄積した。【０１１９】
（２）鑑別工程
本実施例１の鑑別においては、セルロース系繊維の単一繊維からなる５つの被検繊維Ｚ１〜Ｚ５を準備した。まず、判別モデルの作成と同様にして、鑑別しようとする被検繊維Ｚ１〜Ｚ５の吸収スペクトルを求め、微分スペクトルを得た。【０１２０】
まず、被検繊維Ｚ１〜Ｚ５が「天然繊維」と「再生繊維」のいずれのグループに属するかを鑑別した。具体的には、求めた微分スペクトルから判別モデルＡと同様の波数域（波数１２００〜８５０ｃｍ−１）のスペクトルデータＹを抽出した。このようにして得られた被検繊維Ｚ１〜Ｚ５の各スペクトルデータＹから上記と同じ手法によりスコア（射影ｔ１及び射影ｔ２に対応）を算出した。これら被検繊維Ｚ１〜Ｚ５のスコアを図５の「天然繊維」と「再生繊維」の判別モデルＡにプロットした（図５のＺ１〜Ｚ５）。図５において、本実施例１の被検繊維Ｚ１〜Ｚ５のうち、被検繊維Ｚ１とＺ２は「天然繊維」の第１グループ（Ａ１）に属する繊維であることが分かる。一方、被検繊維Ｚ３〜Ｚ５は「再生繊維」の第２グループ（Ａ２）に属する繊維であることが分かる。【０１２１】
次に、「天然繊維」と鑑別された被検繊維Ｚ１、Ｚ２が「綿」と「麻類」のいずれのグループに属するかを鑑別した。具体的には、求めた微分スペクトルから判別モデルＢと同様の波数域（波数１４００〜９００ｃｍ−１）のスペクトルデータＹを抽出した。このようにして得られた被検繊維Ｚ１、Ｚ２の各スペクトルデータＹから上記と同じ手法によりスコア（射影ｔ１及び射影ｔ２に対応）を算出した。これら被検繊維Ｚ１、Ｚ２のスコアを図６の「綿」と「麻類」の判別モデルＢにプロットした（図６のＺ１、Ｚ２）。図６において、本実施例１の被検繊維Ｚ１、Ｚ２のうち、被検繊維Ｚ１は「綿」の第３グループ（Ｂ１）に属する繊維であることが分かる。一方、被検繊維Ｚ２は「麻類」の第４グループ（Ｂ２）に属する繊維であることが分かる。【０１２２】
次に、「麻類」と鑑別された被検繊維Ｚ２が「リネン」と「ラミー」のいずれのグループに属するかを鑑別した。具体的には、求めた微分スペクトルから判別モデルＣと同様の波数域（波数１４００〜９００ｃｍ−１）のスペクトルデータＹを抽出した。このようにして得られた被検繊維Ｚ２の各スペクトルデータＹから上記と同じ手法によりスコア（射影ｔ１及び射影ｔ２に対応）を算出した。この被検繊維Ｚ２のスコアを図７の「リネン」と「ラミー」の判別モデルＣにプロットした（図７のＺ２）。図７において、本実施例１の被検繊維Ｚ２は「リネン」の第５グループ（Ｃ１）に属する繊維であることが分かる。【０１２３】
一方、「再生繊維」と鑑別された被検繊維Ｚ３〜Ｚ５が「レーヨン」と「キュプラ及びリヨセル」のいずれのグループに属するかを鑑別した。具体的には、求めた微分スペクトルから判別モデルＤと同様の波数域（波数１４００〜９００ｃｍ−１）のスペクトルデータＹを抽出した。このようにして得られた被検繊維Ｚ３〜Ｚ５の各スペクトルデータＹから上記と同じ手法によりスコア（射影ｔ１及び射影ｔ２に対応）を算出した。これら被検繊維Ｚ３〜Ｚ５のスコアを図８の「レーヨン」と「キュプラ及びリヨセル」の判別モデルＤにプロットした（図８のＺ３〜Ｚ５）。図８において、本実施例１の被検繊維Ｚ３〜Ｚ５のうち、被検繊維Ｚ３は「レーヨン」の第７グループ（Ｄ１）に属する繊維であることが分かる。一方、被検繊維Ｚ４、Ｚ５は「キュプラ及びリヨセル」の第８グループ（Ｄ２）に属する繊維であることが分かる。【０１２４】
次に、「キュプラ及びリヨセル」と鑑別された被検繊維Ｚ４、Ｚ５が「キュプラ」と「リヨセル」のいずれのグループに属するかを鑑別した。具体的には、求めた微分スペクトルから判別モデルＥと同様の波数域（波数１４００〜９００ｃｍ−１）のスペクトルデータＹを抽出た。このようにして得られた被検繊維Ｚ４、Ｚ５のスペクトルデータＹから上記と同じ手法によりスコア（射影ｔ１及び射影ｔ２に対応）を算出した。これら被検繊維Ｚ４、Ｚ５のスコアを図９の「キュプラ」と「リヨセル」の判別モデルＥにプロットした（図９のＺ４、Ｚ５）。図９において、本実施例１の被検繊維Ｚ４、Ｚ５のうち、被検繊維Ｚ４は「キュプラ」の第９グループ（Ｅ１）に属する繊維であることが分かる。一方、被検繊維Ｚ５、は「リヨセル」の第１０グループ（Ｅ２）に属する繊維であることが分かる。【０１２５】
以上説明したように、本第１実施形態においては、単一繊維からなる被検繊維に対して、繊維の種類を容易かつ正確に鑑別することができた。よって、本発明においては、鑑別操作が比較的簡単で客観性を有し、検査員の経験やノウハウに頼ることなく同系異種繊維の鑑別をすることができ、且つ、鑑別する繊維の組み合わせによって生じる鑑別精度のバラツキを極力小さくして、高度な鑑別精度を実現することのできる繊維鑑別方法を提供することができる。【０１２６】
《第２実施形態》
本第２実施形態においては、上記第１実施例と同様に、複数種類の単一繊維に対して、各単一繊維間の鑑別を行うものである。例えば、セルロース系繊維においては、天然繊維と再生繊維との鑑別、綿と麻類との鑑別、麻類におけるリネンとラミーとの鑑別、再生繊維どうしの鑑別、及び、キュプラとリヨセルとの鑑別などについて説明する。なお、本第２実施形態においては、判別モデルのスコア群の正規分布に対する被検繊維のスコアの確率密度から鑑別を行うものである。【０１２７】
（１）判別モデル準備工程
本第２実施形態における判別モデル準備工程の各操作は、基本的には上記第１実施形態と同様である。まず、比較繊維の各グループの吸収スペクトルを求める。吸収スペクトルの補正についても、上記第１実施形態と同様である。また、アルカリ性物質による前処理をする方法についても、上記第１実施形態と同様である。【０１２８】
本第２実施形態においては、上記第１実施例と同様に、比較繊維の吸収スペクトルを求める際にＦＴ／ＩＲ分光光度計でＡＴＲ法を使用する。また、上記第１実施形態と同様に、近赤外線を除く波数４０００ｃｍ−１〜６００ｃｍ−１の範囲内の赤外線を照射してそれぞれの吸収スペクトルを求める。また、スペクトルデータの微分処理などに関しても、上記第１実施形態と同様である。【０１２９】
次に、得られた微分スペクトルから解析に有効な所定の波数域のスペクトルデータを抽出する。本第２実施形態においては、比較繊維の各組み合わせに対して抽出したスペクトルデータの波数域は、基本的には上記第１実施形態と同様である。但し、本第２実施形態においては、更に狭い波数域に分割して複数の波数域を抽出したものもある。なお、使用した波数域については、下記の実施例２において詳述する。【０１３０】
次に、抽出したスペクトルデータを判別分析して判別モデルを作成する。本第２実施形態においては、上記第１実施例と同様に、本発明者らが独自の工夫を加えた判別分析の手法であるＦＤＯＤ（Fisher’s linear discriminant orthogonal decomposition）を使用した。また、本第２実施形態においても、２回の判別分析を行う。【０１３１】
まず、各比較繊維の吸収スペクトルから特定波数域のスペクトルデータＸを抽出し、これをデータＸ１とする。このデータＸ１を用いて上記第１実施形態と同じ式（１）により第１回判別分析を行って第１軸ｗ１を求める。このとき、予め選定した正則化係数ζ１を使用する。なお、本第２実施形態においても、発明者自らが構成したプログラムを使用して解析した。【０１３２】
次に、第１回判別分析で得られた第１軸ｗ１とデータＸ１とから、上記第１実施形態と同じ式（２）によりデータＸ１の第１軸座標（ｗ１軸座標）の値ｔ１（射影ｔ１ともいう）を算出する。次に、上記第１実施形態と同じ式（３）により先のデータＸ１からｗ１軸方向の情報を引き抜いたものをデータＸ２（残差行列ともいう）とする。次に、このデータＸ２を用いて上記第１実施形態と同じ式（１）により第２回判別分析を行って第２軸ｗ２を求める。このとき、予め選定した正則化係数ζ２を使用する。このようにして求めた第２軸ｗ２とデータＸ２とから、上記第１実施形態と同じ式（４）によりデータＸ２の第２軸座標（ｗ２軸座標）の値ｔ２（射影ｔ２ともいう）を算出する。【０１３３】
このようにして求めた互いに直交する２つの軸ｗ１、ｗ２を用いて、直交分解のなされた判別モデルを作成することができる。この判別モデルにおいては、軸ｗ１、ｗ２を直交軸として２つのグループが明確に区別される。各グループに属する各比較繊維の射影ｔ１及び射影ｔ２に対応する値をその比較繊維のスコアという。本第２実施形態においては、判別モデルの２つのグループについて、それぞれスペクトルデータから求めたスコア群の２次元の正規分布の推定を行う。具体的には、判別モデルの２つのグループ（グループａ、グループｂ）について、それぞれ各スコア群の２次元の正規分布の確率密度関数Ｆa(ｘ)、Ｆb(ｘ)を求める。なお、これらの判別モデル及び確率密度関数Ｆa(ｘ)、Ｆb(ｘ)については、下記の実施例２において詳述する。【０１３４】
また、本第２実施形態においては、上述のように互いに直交する２つの軸ｗ１、ｗ２を用いて、直交分解のなされた判別モデルを作成するものであるが、上記式（３）及び式（４）の操作を繰り返すことにより、互いに直交する３軸以上の高次元の判別モデルを作成するようにしてもよい。なお、２次元ではなくｎ次元の判別モデルを作成した場合には、各スコア群からｎ次元（ｎは１以上の整数）の正規分布の推定を行い、ｎ次元の確率密度関数を求めるようにする。【０１３５】
（２）鑑別工程
本第２実施形態の鑑別工程においては、まず、鑑別しようとする被検繊維Ｚの吸収スペクトルを求める。吸収スペクトルを求める方法、吸収スペクトルに各種補正をする方法、及び、得られた吸収スペクトルに微分処理などの処理を行う方法は、上述の比較繊維に対する方法と同様である。次に、求めた微分スペクトルから比較繊維と同じ波数域のスペクトルデータＹを抽出し、比較繊維と同様にして被検繊維Ｚのスコア（射影ｔ１及び射影ｔ２に対応）を算出する。次に、得られた被検繊維Ｚのスコアを用いて、判別モデルの２つのグループ（グループａ、グループｂ）それぞれの確率密度関数Ｆa(ｘ)、Ｆb(ｘ)から各グループに対する被検繊維Ｚの確率密度ｄz,a、ｄz,bを算出する。【０１３６】
次に、得られた被検繊維の確率密度ｄz,a、ｄz,bから、下記式（５）により被検繊維の一方のグループ（例えば、グループａ）に対する確率密度比Ｒz,aを算出する。【０１３７】
Ｒz,a＝ｄz,a／（ｄz,a＋ｄz,b）  ・・・（５）
本第２実施形態においては、このようにして得られた確率密度比Ｒz,aによって、被検繊維Ｚがいずれのグループ（グループａ又はグループｂ）に属する繊維であるかを判別する。【０１３８】
まず、被検繊維Ｚのグループａに対する確率密度比Ｒz,a（又は、グループｂに対する確率密度比Ｒz,b）による任意の判別基準を設定する。例えば、グループａに対する確率密度比Ｒz,aによる判別基準を、
Ｒz,a≧０．９  ・・・グループａに属する
Ｒz,a≦０．１  ・・・グループｂに属する
０．１＜Ｒz,a＜０．９  ・・・判別不能
とする。【０１３９】
この場合、Ｒz,a≧０．９であれば、被検繊維Ｚのグループａに属する確率は９０％以上であって、グループｂに属する確率は１０％未満と判別できる。逆に、Ｒz,a≦０．１であれば、被検繊維Ｚのグループｂに属する確率は９０％より大きく、グループａに属する確率は１０％以下と判別できる。一方、０．１＜Ｒz,a＜０．９であれば、判別不能と判断する。なお、判別基準は、９０％に限るものではなく、任意の基準であってもよい。一般に、８５％〜９５％程度であることが好ましい。【０１４０】
図１０は、本第２実施形態において被検繊維を鑑別するイメージを表す概念図である。図１０において、グループａの分布とグループｂの分布に対して、確率密度比Ｒz,a＝０．９の２本の判別ラインＬ１、Ｌ２と、確率密度比Ｒz,a＝０．１の２本の判別ラインＬ３、Ｌ４とが存在する。グループａに近い判別ラインＬ１のグループａの側の領域Ｇ１は、グループａと判別される。また、グループｂに近い２本の判別ラインＬ３と判別ラインＬ４とで囲まれた領域Ｇ３は、グループｂと判別される。これに対して、判別ラインＬ１と判別ラインＬ３とで囲まれた領域Ｇ２は、０．１＜Ｒz,a＜０．９の領域であって判別不能と判別される。同様に、判別ラインＬ２と判別ラインＬ４とで囲まれた領域Ｇ４も、０．１＜Ｒz,a＜０．９の領域であって判別不能と判別される。このような被検繊維は、２種以上の混合繊維であるか、或いは、セルロース系繊維以外の繊維である可能性がある。【０１４１】
一方、グループａの分布とグループｂの分布から大きく離れた領域Ｇ５が存在する。この領域Ｇ５においては、確率密度ｄz,a、ｄz,bの値がいずれも小さくなり、確率密度比Ｒz,aの値が０．９以上又は０．１以下となることがある。このような被検繊維は、２種以上の混合繊維であるか、或いは、セルロース系繊維以外の繊維である可能性がある。この場合、領域Ｇ５はグループｂに相対的に近いにもかかわらず、確率密度比Ｒz,aの値からはグループａと判別されることになる。このような誤判別を避けるために、本第２実施形態においては、確率密度比Ｒz,aの値のみから判断するのではなく、判別モデルのスコア群のＦＤＯＤプロット及び必要により等確率楕円により確認することが好ましい。
【実施例２】
【０１４２】
次に、本第２実施形態の鑑別方法を実施例２により具体的に説明する。本実施例２は、上記実施例１と同様に各単一繊維間の鑑別を行うものであり、複数の被検繊維に対して鑑別フロー図（図３参照）に従って鑑別を行うものである。なお、各被検繊維は、いずれも顕微鏡法などの予備鑑定においてセルロース系繊維であることが判明している。【０１４３】
（１）判別モデル準備工程
本実施例２においては、まず、繊維の種類が既知の単一繊維として、綿２５点、リネン２５点、ラミー２５点の天然繊維７５点を準備した。また、繊維の種類が既知の単一繊維として、レーヨン４５点、キュプラ２５点、リヨセル２５点の再生繊維９５点を準備した。これら合計１７０点の単一繊維を本実施例２の比較繊維とした。なお、綿及び麻類（リネン及びラミー）７５点、並びに、再生繊維の一部４０点（レーヨン１４点、キュプラ１３点、リヨセル１３点）に関しては、室温において水酸化ナトリウム水溶液（１７重量％）による前処理を行った。【０１４４】
次に、これら合計１７０点の織編物について、その吸収スペクトルを得た。吸収スペクトルの測定は、上記実施例１と同様に赤外分光光度計ＦＴ／ＩＲ‐４７００（日本分光株式会社）を使用し、ＺｎＳｅプリズムによるＡＴＲ法で、波数４０００ｃｍ−１〜６００ｃｍ−１の吸収スペクトルを測定した。次に、乗算的散乱補正（ＭＳＣ）後の各吸収スペクトルをSavitzky-Golay法により１次微分して微分スペクトルを得た。得られた吸収スペクトル（微分スペクトル）をＦＤＯＤによる判別モデル作成用グループとして分類した。【０１４５】
［判別モデル作成用グループ］
・第１グループ（Ａ３）：４０点の天然繊維（綿２０点、リネン１０点及びラミー１０点）
・第２グループ（Ａ４）：４０点の再生繊維（前処理レーヨン１４点、前処理キュプラ１３点及び前処理リヨセル１３点）
・第３グループ（Ｂ３）：２５点の綿
・第４グループ（Ｂ４）：２５点の麻類（リネン１２点及びラミー１３点）
・第５グループ（Ｃ３）：２５点のリネン
・第６グループ（Ｃ４）：２５点のラミー
・第７グループ（Ｄ３）：４５点の未処理レーヨン   
・第８グループ（Ｄ４）：４５点のキュプラ及びリヨセル（未処理キュプラ２４点及び未処理リヨセル２１点）
・第９グループ（Ｅ３）：２５点の未処理キュプラ
・第１０グループ（Ｅ４）：２５点の未処理リヨセル
このような上記１０グループを得た。【０１４６】
次に、判別モデル作成用グループを用いて判別モデルを作成した。上記第１グループ（Ａ３）〜第１０グループ（Ｅ４）のうち２つずつ組み合わせて、５組の組み合わせとした。次に、これらの微分スペクトルから各組み合わせを特徴付ける所定の波数域を抽出した。更に、各組み合わせに対する判別分析用の２つの正則化係数ζ１、ζ２を特定した。本実施例２において、各判別モデルＡ〜Ｅと抽出した波数域及び２つの正則化係数ζ１、ζ２の値を下記に示す。但し、本実施例２においては、抽出した波数域が上記実施例１とは異なり、狭く区分した波数域を複数抽出したものもある。【０１４７】
［判別モデルＡ］：天然繊維と再生繊維
・第１グループ（Ａ３）と第２グループ（Ａ４）：波数１３００〜９００ｃｍ−１
・第１回判別分析用の正則化係数：ζ１＝１×１０−６
・第２回判別分析用の正則化係数：ζ２＝１×１０−３
［判別モデルＢ］：綿と麻類
・第３グループ（Ｂ３）と第４グループ（Ｂ４）：波数１６００〜１４００ｃｍ−１、波数１２００〜１１００ｃｍ−１、波数１０００〜８００ｃｍ−１
・第１回判別分析用の正則化係数：ζ１＝１×１０−６
・第２回判別分析用の正則化係数：ζ２＝１×１０−４
［判別モデルＣ］：リネンとラミー
・第５グループ（Ｃ３）と第６グループ（Ｃ４）：波数１６００〜１５００ｃｍ−１、波数１４００〜１１００ｃｍ−１、波数１０００〜９００ｃｍ−１
・第１回判別分析用の正則化係数：ζ１＝１×１０−６
・第２回判別分析用の正則化係数：ζ２＝１×１０−６
［判別モデルＤ］：レーヨンと（キュプラ及びリヨセル）
・第７グループ（Ｄ３）と第８グループ（Ｄ４）：波数１５００〜１４００ｃｍ−１、波数１３００〜１０００ｃｍ−１、波数９００〜８００ｃｍ−１
・第１回判別分析用の正則化係数：ζ１＝１×１０−４
・第２回判別分析用の正則化係数：ζ２＝１×１０−７
［判別モデルＥ］：キュプラとリヨセル
・第９グループ（Ｅ３）と第１０グループ（Ｅ４）：波数１７００〜１６００ｃｍ−１、波数１５００〜１３００ｃｍ−１、波数１１００〜１０００ｃｍ−１ 
・第１回判別分析用の正則化係数：ζ１＝１×１０−６
・第２回判別分析用の正則化係数：ζ２＝１×１０−５
このようにして抽出した波数域における各グループのスペクトルデータＸを第１回目の判別分析に使用するデータＸ１とした。このデータＸ１を使用して５組の組み合わせに対する２回の判別分析を行い、それぞれの比較繊維のスコア（射影ｔ１及び射影ｔ２に対応）を算出した。【０１４８】
本実施例２においては、判別する２つのグループに含まれる比較繊維の２つのスコア群が判別モデルに対応する。例えば、判別モデルＡにおいては、第１グループ（Ａ３）と第２グループ（Ａ４）の各比較繊維のスコア群の組み合わせに対応する。このようにして、判別モデル（判別する各比較繊維のスコア群）Ａ〜Ｅを得た。従って、本実施例２においては、各グループの比較繊維のスコア群を互いに直交する２つの軸ｗ１、ｗ２にプロットした判別モデル（ＦＤＯＤプロット）Ａ〜Ｅを作成することを必須の要件とするものではない。但し、判別したグループに属すること、又はいずれのグループにも属さないこと（判別不可）を確認するためにＦＤＯＤプロット及び必要により等確率楕円を作成しておくことが好ましい。【０１４９】
次に、判別モデル（判別する各比較繊維のスコア群）Ａ〜Ｅの各グループ（Ａ３〜Ｅ４）について、それぞれスコア群の２次元の正規分布の推定を行った。例えば、判別モデルＡにおいては、２つのグループＡ３、Ａ４の各スコア群について、２次元の正規分布の確率密度関数ＦA3(ｘ)、ＦA4(ｘ)を求めた。確率密度関数を求めるにあたっては、通常の方法を使用すればよい。同様にして、判別モデルＢにおいて確率密度関数ＦB3(ｘ)、ＦB4(ｘ)、判別モデルＣにおいて確率密度関数ＦC3(ｘ)、ＦC4(ｘ)、判別モデルＤにおいて確率密度関数ＦD3(ｘ)、ＦD4(ｘ)、判別モデルＥにおいて確率密度関数ＦE3(ｘ)、ＦE4(ｘ)を求めた。【０１５０】
（２）鑑別工程
本実施例２の鑑別においては、上記実施例１と同じ、セルロース系繊維の単一繊維からなる５つの被検繊維Ｚ６〜Ｚ１０を準備した。まず、判別モデルの作成と同様にして、鑑別しようとする被検繊維Ｚ６〜Ｚ１０の吸収スペクトルを求め、微分スペクトルを得た。【０１５１】
まず、被検繊維Ｚ６〜Ｚ１０が「天然繊維」と「再生繊維」のいずれのグループに属するかを鑑別した。具体的には、求めた微分スペクトルから判別モデルＡと同様の波数域（波数１３００〜９００ｃｍ−１）のスペクトルデータＹを抽出した。このようにして得られた被検繊維Ｚ６〜Ｚ１０の各スペクトルデータＹから上記と同じ手法によりスコア（射影ｔ１及び射影ｔ２に対応）を算出した。【０１５２】
得られた被検繊維Ｚ６のスコアを用いて、確率密度関数ＦA3(ｘ)、ＦA4(ｘ)から各グループＡ３、Ａ４に対する被検繊維Ｚ６の確率密度ｄZ6,A3、ｄZ6,A4を算出した。同様にして、被検繊維Ｚ７の確率密度ｄZ7,A3、ｄZ7,A4、被検繊維Ｚ８の確率密度ｄZ8,A3、ｄZ8,A4、被検繊維Ｚ９の確率密度ｄZ9,A3、ｄZ9,A4、被検繊維Ｚ１０の確率密度ｄZ10,A3、ｄZ10,A4を算出した。算出した各確率密度の値を表１に示す。【０１５３】
得られた被検繊維Ｚ６の確率密度ｄZ6,A3、ｄZ6,A4から、下記式（５−１）により被検繊維Ｚ６のグループＡ３に対する確率密度比ＲZ6,A3を算出した。【０１５４】
ＲZ6,A3＝ｄZ6,A3／（ｄZ6,A3＋ｄZ6,A4）  ・・・（５−１）
同様にして、被検繊維Ｚ７のグループＡ３に対する確率密度比ＲZ7,A3、被検繊維Ｚ８のグループＡ３に対する確率密度比ＲZ8,A3、被検繊維Ｚ９のグループＡ３に対する確率密度比ＲZ9,A3、被検繊維Ｚ１０のグループＡ３に対する確率密度比ＲZ10,A3を算出した。算出した各確率密度比の値を表１に示す。【０１５５】
【表１】
0006771576.tif 000002
【０１５６】
表１において、グループＡ３に対する確率密度比の判別基準を次のように定めた場合、
ＲZ,A3≧０．９  ・・・グループＡ３に属する
ＲZ,A3≦０．１  ・・・グループＡ４に属する
０．１＜ＲZ,A3＜０．９  ・・・判別不能
本実施例２の被検繊維Ｚ６〜Ｚ１０のうち、被検繊維Ｚ６とＺ７は「天然繊維」の第１グループ（Ａ３）に属する繊維であることが分かる。一方、被検繊維Ｚ８〜Ｚ１０は「再生繊維」の第２グループ（Ａ４）に属する繊維であることが分かる。【０１５７】
なお、本実施例２においては、上述のように、判別結果を確認するためにＦＤＯＤプロット及びこれに基づいた９９．７％信頼水準の等確率楕円を作成しておいた。図１１は、本実施例２で得られた「天然繊維」と「再生繊維」のＦＤＯＤプロットである。図１１において、第１判別軸（ｗ１軸）を横軸とし第２判別軸（ｗ２軸）を縦軸として、「天然繊維」と「再生繊維」の２つのグループＡ３、Ａ４が明確に層別された。また、被検繊維Ｚ６〜Ｚ１０の各プロットが２つのグループＡ３、Ａ４のいずれか一方の等確率楕円の中に納まっており、確率密度比による判別結果と一致した。このことから、本実施例２の確率密度比による判別結果は、正確なものであり誤判別ではないことが確認された。【０１５８】
次に、「天然繊維」と鑑別された被検繊維Ｚ６、Ｚ７が「綿」と「麻類」のいずれのグループに属するかを鑑別した。具体的には、求めた微分スペクトルから判別モデルＢと同様の波数域（波数１６００〜１４００ｃｍ−１、波数１２００〜１１００ｃｍ−１、波数１０００〜８００ｃｍ−１）のスペクトルデータＹを抽出した。このようにして得られた被検繊維Ｚ６、Ｚ７の各スペクトルデータＹから上記と同じ手法によりスコア（射影ｔ１及び射影ｔ２に対応）を算出した。【０１５９】
得られた被検繊維Ｚ６のスコアを用いて、確率密度関数ＦB3(ｘ)、ＦB4(ｘ)から各グループＢ３、Ｂ４に対する被検繊維Ｚ６の確率密度ｄZ6,B3、ｄZ6,B4を算出した。同様にして、被検繊維Ｚ７の確率密度ｄZ7,B3、ｄZ7,B4を算出した。算出した各確率密度の値を表２に示す。次に、得られた被検繊維Ｚ６の確率密度ｄZ6,B3、ｄZ6,B4から、下記式（５−２）により被検繊維Ｚ６のグループＢ３に対する確率密度比ＲZ6,B3を算出した。【０１６０】
ＲZ6,B3＝ｄZ6,B3／（ｄZ6,B3＋ｄZ6,B4）  ・・・（５−２）
同様にして、被検繊維Ｚ７のグループＢ３に対する確率密度比ＲZ7,B3を算出した。算出した各確率密度比の値を表２に示す。【０１６１】
【表２】
0006771576.tif 000003
【０１６２】
表２において、グループＢ３に対する確率密度比の判別基準を次のように定めた場合、
ＲZ,B3≧０．９  ・・・グループＢ３に属する
ＲZ,B3≦０．１  ・・・グループＢ４に属する
０．１＜ＲZ,B3＜０．９  ・・・判別不能
本実施例２の被検繊維Ｚ６、Ｚ７のうち、被検繊維Ｚ６は「綿」の第３グループ（Ｂ３）に属する繊維であることが分かる。一方、被検繊維Ｚ７は「麻類」の第４グループ（Ｂ４）に属する繊維であることが分かる。【０１６３】
図１２は、本実施例２で得られた「綿」と「麻類」のＦＤＯＤプロットである。図１２において、第１判別軸（ｗ１軸）を横軸とし第２判別軸（ｗ２軸）を縦軸として、「綿」と「麻類」の２つのグループＢ３、Ｂ４が明確に層別された。また、被検繊維Ｚ６、Ｚ７の各プロットが２つのグループＢ３、Ｂ４のいずれか一方の等確率楕円の中に納まっており、確率密度比による判別結果と一致した。このことから、本実施例２の確率密度比による判別結果は、正確なものであり誤判別ではないことが確認された。【０１６４】
次に、「麻類」と鑑別された被検繊維Ｚ７が「リネン」と「ラミー」のいずれのグループに属するかを鑑別した。具体的には、求めた微分スペクトルから判別モデルＣと同様の波数域（波数１６００〜１５００ｃｍ−１、波数１４００〜１１００ｃｍ−１、波数１０００〜９００ｃｍ−１）のスペクトルデータＹを抽出した。このようにして得られた被検繊維Ｚ７の各スペクトルデータＹから上記と同じ手法によりスコア（射影ｔ１及び射影ｔ２に対応）を算出した。【０１６５】
得られた被検繊維Ｚ７のスコアを用いて、確率密度関数ＦC3(ｘ)、ＦC4(ｘ)から各グループＣ３、Ｃ４に対する被検繊維Ｚ７の確率密度ｄZ7,C3、ｄZ7,C4を算出した。算出した各確率密度の値を表３に示す。次に、得られた被検繊維Ｚ７の確率密度ｄZ7,C3、ｄZ7,C4から、下記式（５−３）により被検繊維Ｚ７のグループＣ３に対する確率密度比ＲZ7,C3を算出した。【０１６６】
ＲZ7,C3＝ｄZ7,C3／（ｄZ7,C3＋ｄZ7,C4）  ・・・（５−３）
算出した確率密度比の値を表３に示す。【０１６７】
【表３】
0006771576.tif 000004
【０１６８】
表３において、グループＣ３に対する確率密度比の判別基準を次のように定めた場合、
ＲZ,C3≧０．９  ・・・グループＣ３に属する
ＲZ,C3≦０．１  ・・・グループＣ４に属する
０．１＜ＲZ,C3＜０．９  ・・・判別不能
本実施例２の被検繊維Ｚ７は「ラミー」の第６グループ（Ｃ４）に属する繊維であることが分かる。【０１６９】
図１３は、本実施例２で得られた「リネン」と「ラミー」のＦＤＯＤプロットである。図１３において、第１判別軸（ｗ１軸）を横軸とし第２判別軸（ｗ２軸）を縦軸として、「リネン」と「ラミー」の２つのグループＣ３、Ｃ４が明確に層別された。また、被検繊維Ｚ７のプロットがグループＣ４の等確率楕円の中に納まっており、確率密度比による判別結果と一致した。このことから、本実施例２の確率密度比による判別結果は、正確なものであり誤判別ではないことが確認された。【０１７０】
一方、「再生繊維」と鑑別された被検繊維Ｚ８〜Ｚ１０が「レーヨン」と「キュプラ及びリヨセル」のいずれのグループに属するかを鑑別した。具体的には、求めた微分スペクトルから判別モデルＤと同様の波数域（波数１５００〜１４００ｃｍ−１、波数１３００〜１０００ｃｍ−１、波数９００〜８００ｃｍ−１）のスペクトルデータＹを抽出した。このようにして得られた被検繊維Ｚ８〜Ｚ１０の各スペクトルデータＹから上記と同じ手法によりスコア（射影ｔ１及び射影ｔ２に対応）を算出した。【０１７１】
得られた被検繊維Ｚ８のスコアを用いて、確率密度関数ＦD3(ｘ)、ＦD4(ｘ)から各グループＤ３、Ｄ４に対する被検繊維Ｚ８の確率密度ｄZ8,D3、ｄZ8,D4を算出した。同様にして、被検繊維Ｚ９の確率密度ｄZ9,D3、ｄZ9,D4、被検繊維Ｚ１０の確率密度ｄZ10,D3、ｄZ10,D4を算出した。算出した各確率密度の値を表４に示す。次に、得られた被検繊維Ｚ８の確率密度ｄZ8,D3、ｄZ8,D4から、下記式（５−４）により被検繊維Ｚ８のグループＤ３に対する確率密度比ＲZ8,D3を算出した。【０１７２】
ＲZ8,D3＝ｄZ8,D3／（ｄZ8,D3＋ｄZ8,D4）  ・・・（５−４）
同様にして、被検繊維Ｚ９のグループＤ３に対する確率密度比ＲZ9,D3、被検繊維Ｚ１０のグループＤ３に対する確率密度比ＲZ10,D3を算出した。算出した各確率密度比の値を表４に示す。【０１７３】
【表４】
0006771576.tif 000005
【０１７４】
表４において、グループＤ３に対する確率密度比の判別基準を次のように定めた場合、
ＲZ,D3≧０．９  ・・・グループＤ３に属する
ＲZ,D3≦０．１  ・・・グループＤ４に属する
０．１＜ＲZ,D3＜０．９  ・・・判別不能
本実施例２の被検繊維Ｚ８〜Ｚ１０のうち、被検繊維Ｚ８は「レーヨン」の第７グループ（Ｄ３）に属する繊維であることが分かる。一方、被検繊維Ｚ９、Ｚ１０は「キュプラ及びリヨセル」の第８グループ（Ｄ４）に属する繊維であることが分かる。【０１７５】
図１４は、本実施例２で得られた「レーヨン」と「キュプラ及びリヨセル」のＦＤＯＤプロットである。図１４において、第１判別軸（ｗ１軸）を横軸とし第２判別軸（ｗ２軸）を縦軸として、「レーヨン」と「キュプラ及びリヨセル」の２つのグループＤ３、Ｄ４が明確に層別された。また、被検繊維Ｚ８〜Ｚ１０の各プロットが２つのグループＢ３、Ｂ４のいずれか一方の等確率楕円の中に納まっており、確率密度比による判別結果と一致した。このことから、本実施例２の確率密度比による判別結果は、正確なものであり誤判別ではないことが確認された。【０１７６】
次に、「キュプラ及びリヨセル」と鑑別された被検繊維Ｚ９、Ｚ１０が「キュプラ」と「リヨセル」のいずれのグループに属するかを鑑別した。具体的には、求めた微分スペクトルから判別モデルＥと同様の波数域（波数１７００〜１６００ｃｍ−１、波数１５００〜１３００ｃｍ−１、波数１１００〜１０００ｃｍ−１）のスペクトルデータＹを抽出た。このようにして得られた被検繊維Ｚ９、Ｚ１０のスペクトルデータＹから上記と同じ手法によりスコア（射影ｔ１及び射影ｔ２に対応）を算出した。【０１７７】
得られた被検繊維Ｚ９のスコアを用いて、確率密度関数ＦE3(ｘ)、ＦE4(ｘ)から各グループＥ３、Ｅ４に対する被検繊維Ｚ９の確率密度ｄZ9,E3、ｄZ9,E4を算出した。同様にして、被検繊維Ｚ１０の確率密度ｄZ10,E3、ｄZ10,E4を算出した。算出した各確率密度の値を表５に示す。次に、得られた被検繊維Ｚ９の確率密度ｄZ9,E3、ｄZ9,E4から、下記式（５−５）により被検繊維Ｚ９のグループＥ３に対する確率密度比ＲZ9,E3を算出した。【０１７８】
ＲZ9,E3＝ｄZ9,E3／（ｄZ9,E3＋ｄZ9,E4）  ・・・（５−５）
同様にして、被検繊維Ｚ１０のグループＥ３に対する確率密度比ＲZ10,E3を算出した。算出した各確率密度比の値を表５に示す。【０１７９】
【表５】
0006771576.tif 000006
【０１８０】
表５において、グループＥ３に対する確率密度比の判別基準を次のように定めた場合、
ＲZ,E3≧０．９  ・・・グループＥ３に属する
ＲZ,E3≦０．１  ・・・グループＥ４に属する
０．１＜ＲZ,E3＜０．９  ・・・判別不能
本実施例２の被検繊維Ｚ９、Ｚ１０のうち、被検繊維Ｚ９は「キュプラ」の第９グループ（Ｅ３）に属する繊維であることが分かる。一方、被検繊維Ｚ１０、は「リヨセル」の第１０グループ（Ｅ４）に属する繊維であることが分かる。【０１８１】
図１５は、本実施例２で得られた「キュプラ」と「リヨセル」のＦＤＯＤプロットである。図１５において、第１判別軸（ｗ１軸）を横軸とし第２判別軸（ｗ２軸）を縦軸として、「キュプラ」と「リヨセル」の２つのグループＥ３、Ｅ４が明確に層別された。また、被検繊維Ｚ９、Ｚ１０の各プロットが２つのグループＥ３、Ｅ４のいずれか一方の等確率楕円の中に納まっており、確率密度比による判別結果と一致した。このことから、本実施例２の確率密度比による判別結果は、正確なものであり誤判別ではないことが確認された。【０１８２】
以上説明したように、本第２実施形態においては、単一繊維からなる被検繊維に対して、繊維の種類を容易かつ正確に鑑別することができた。よって、本発明においては、鑑別操作が比較的簡単で客観性を有し、検査員の経験やノウハウに頼ることなく同系異種繊維の鑑別をすることができ、且つ、鑑別する繊維の組み合わせによって生じる鑑別精度のバラツキを極力小さくして、高度な鑑別精度を実現することのできる繊維鑑別方法を提供することができる。【０１８３】
なお、本発明の実施にあたり、上記各実施形態及び各実施例に限らず次のような種々の変形例が挙げられる。
（１）上記各実施形態においては、化学的組成が同じ繊維どうしの鑑別例としてセルロース系繊維で説明するものであるが、これに限るものではなく、タンパク質系繊維としての絹や各種獣毛繊維どうしの鑑別や、その他の化学的組成が同じ繊維どうしの鑑別を行うようにしてもよい。
（２）上記各実施形態においては、２つのグループ間の鑑別を判別分析するものであるが、これに限るものではなく、３つ或いはそれ以上のグループ間の鑑別を同時に判別分析するようにしてもよい。
（３）上記各実施形態においては、ＩＲ分光分析により得られる吸収スペクトルにより鑑別を行うものであるが、これに限るものではなく、ＩＲ分光分析により得られる透過スペクトルにより鑑別を行うようにしてもよい。
（４）上記各実施形態においては、波数４０００〜６００ｃｍ−１の範囲内の吸収スペクトルを測定し、これから解析に使用する所定の波数域におけるスペクトルデータを抽出するものであるが、これに限るものではなく、解析に使用する１つ又は２つ以上の所定の波数域における吸収スペクトルのみを測定するようにしてもよい。
（５）上記各実施形態においては、ＩＲ分光分析により得られる吸収スペクトルにより鑑別を行うものであるが、これに限るものではなく、ＮＩＲ分光分析により得られる吸収スペクトルにより鑑別を行うようにしてもよい。このＮＩＲ分光分析により吸収スペクトルを得る場合には、ＩＲ分光分析と同様にＦＴ／ＩＲ分光光度計を使用するようにしてもよく、或いは、ＦＴ／ＮＩＲ分光光度計、分散型ＮＩＲ分光光度計、分散型ＵＶ−ＶＩＳ−ＮＩＲ分光光度計などを使用するようにしてもよい。また、ＮＩＲ分光分析においてこれらの分光光度計を使用して吸収スペクトルを測定する場合には、拡散反射法などを使用することが一般的である。
（６）上記各実施形態においては、ＡＴＲ法により吸収スペクトルを求めるものであるが、これに限るものではなく、その他の方法、例えば、比較繊維或いは被検繊維を微粉砕してからＫＢｒ錠剤法などで吸収スペクトルを求めるようにしてもよい。
（７）上記各実施例においては、赤外吸収の大きな波数域の組み合わせを用いて判別分析を行い、特に有効と考えられる特定の波数域を抽出するものであるが、これに限るものではなく、解析ソフトによる波数域の抽出などを利用して波数域を抽出するようにしてもよい。
（８）上記各実施例においては、互いに直交する第１軸と第２軸との組み合わせによる２次元の分析を使用したが、これらに限るものではなく、互いに直交する３つ以上の軸を用いて３次元以上で分析するようにしてもよい。
（９）上記実施例１においては、天然繊維と再生繊維との鑑別に波数１２００〜８５０ｃｍ−１の範囲内のスペクトルデータのみを使用し、上記実施例２においては、波数１３００〜９００ｃｍ−１のみを使用するものであるが、これに限るものではなく、例えば、波数３５００〜３０００ｃｍ−１の範囲内若しくはその近傍の範囲内のスペクトルデータを組み合わせて使用するようにしてもよい。
（１０）上記実施例１においては、綿と麻類との鑑別に波数１４００〜９００ｃｍ−１の範囲内のスペクトルデータのみを使用し、上記実施例２においては、波数１６００〜１４００ｃｍ−１、波数１２００〜１１００ｃｍ−１、波数１０００〜８００ｃｍ−１の組み合わせを使用するものであるが、これに限るものではなく、例えば、波数３５００〜３０００ｃｍ−１の範囲内若しくはその近傍の範囲内のスペクトルデータを組み合わせて使用するようにしてもよい。
（１１）上記実施例１においては、リネンとラミーとの鑑別に波数１４００〜９００ｃｍ−１の範囲内のスペクトルデータのみを使用し、上記実施例２においては、波数１６００〜１５００ｃｍ−１、波数１４００〜１１００ｃｍ−１、波数１０００〜９００ｃｍ−１の組み合わせを使用するものであるが、これに限るものではなく、例えば、波数３５００〜３０００ｃｍ−１の範囲内若しくはその近傍の範囲内のスペクトルデータを組み合わせて使用するようにしてもよい。
（１２）上記実施例１においては、再生繊維どうしの鑑別に波数１４００〜９００ｃｍ−１の範囲内のスペクトルデータのみを使用し、上記実施例２においては、波数１５００〜１４００ｃｍ−１、波数１３００〜１０００ｃｍ−１、波数９００〜８００ｃｍ−１の組み合わせを使用するものであるが、これに限るものではなく、例えば、波数３５００〜３０００ｃｍ−１の範囲内若しくはその近傍の範囲内のスペクトルデータを組み合わせて使用するようにしてもよい。
（１３）上記実施例１においては、キュプラとリヨセルとの鑑別に波数１４００〜９００ｃｍ−１の範囲内のスペクトルデータのみを使用し、上記実施例２においては、波数１７００〜１６００ｃｍ−１、波数１５００〜１３００ｃｍ−１、波数１１００〜１０００ｃｍ−１の組み合わせを使用するものであるが、これに限るものではなく、例えば、波数３５００〜３０００ｃｍ−１の範囲内若しくはその近傍の範囲内のスペクトルデータを組み合わせて使用するようにしてもよい。
【産業上の利用可能性】
【０１８４】
市場には多くの繊維製品が広い用途に流通している。また、繊維製品の生産地と消費地がグローバルに展開される今日においては、繊維製品の輸出入の際に取引の安全や信頼を確保するために、輸出入の際に迅速、且つ、正確な鑑別方法が望まれている。特に、化学的組成が同じセルロース系繊維どうしの過誤混入や、化学的組成が同じカシミヤなどの高級獣毛繊維と安価な他の獣毛繊維との正確な鑑別が望まれている。【０１８５】
本発明は、このような市場の要求に対して的確な鑑別手段を提供するものであり、また、従来法のように検査員の経験やノウハウに頼ることがない。特に、化学的組成が同じ繊維どうしの種類を客観的に鑑別できること、また、過誤混入や偽装がなされていても正確な鑑別結果が得られるということは、これまでにない画期的な鑑別手段となる。【０１８６】
このことから、本発明は、市場の安定や国際間の公正取引に有効な鑑別手段を提供するものであり、単に従来法であるＪＩＳ Ｌ １０３０‐１（繊維製品の混用率試験方法‐第１部：繊維鑑別）、及び、ＪＩＳ Ｌ １０３０‐２（繊維製品の混用率試験方法‐第２部：繊維混用率）を補完する鑑別手段に留まらず、国際標準として利用可能な鑑別手段を提供することができる。
【符号の説明】
【０１８７】
１…綿、２…麻類、３…レーヨン、４…キュプラ、５…リヨセル、６…結晶セルロース、
Ａ〜Ｅ…判別モデル、
Ａ１、Ａ３…天然繊維のグループ、Ａ２、Ａ４…再生繊維のグループ、
Ｂ１、Ｂ３…綿のグループ、Ｂ２、Ｂ４…麻類のグループ、
Ｃ１、Ｃ３…リネンのグループ、Ｃ２、Ｃ４…ラミーのグループ、
Ｄ１、Ｄ３…レーヨンのグループ、Ｄ２、Ｄ４…キュプラ及びリヨセルのグループ、
Ｅ１、Ｅ３…キュプラのグループ、Ｅ２、Ｅ４…リヨセルのグループ、
Ｇ１〜Ｇ５…領域、L１〜L４…判別ライン、
Ｘ、Y…スペクトルデータ、Ｘ１、Ｘ２…データ、
Ｚ１〜Ｚ１０…被検繊維。
【図１】
0006771576.tif 000007
【図２】
0006771576.tif 000008
【図３】
0006771576.tif 000009
【図４】
0006771576.tif 000010
【図５】
0006771576.tif 000011
【図６】
0006771576.tif 000012
【図７】
0006771576.tif 000013
【図８】
0006771576.tif 000014
【図９】
0006771576.tif 000015
【図１０】
0006771576.tif 000016
【図１１】
0006771576.tif 000017
【図１２】
0006771576.tif 000018
【図１３】
0006771576.tif 000019
【図１４】
0006771576.tif 000020
【図１５】
0006771576.tif 000021