(57)【特許請求の範囲】
【請求項１】
加熱した殻付卵の卵黄粘度判定方法であって、
殻付卵に光を照射する工程と、殻付卵の透過光を受光する工程と、前記透過光のスペクトルを取得する工程と、前記スペクトルに基づき卵黄粘度を判定する工程とを備え、
前記殻付卵に光を照射する工程は、殻付卵の膨らみ部分を通過するように光を照射し、
前記スペクトルに基づき卵黄粘度を判定する工程は、前記スペクトルを多変量解析する工程を含み、該多変量解析がＰＬＳ回帰分析であって、前記スペクトルの可視領域を含む波長を用い、該波長が６００ｎｍから９００ｎｍの領域を含む波長であり、波長の吸光度を抽出したデータと、見かけの卵黄粘度ηを常用対数変換したデータｌｏｇ（η）とに基づくＰＬＳ回帰分析である、
卵黄粘度判定方法。
【請求項２】
前記スペクトルに基づき卵黄粘度を判定する工程は、
前記スペクトルの二次微分値を算出する工程を含む、
請求項１に記載の卵黄粘度判定方法。
【請求項３】
前記スペクトルに基づき卵黄粘度を判定する工程は、
前記スペクトルを平滑化処理及びＳＮＶ処理する工程を含む、
請求項２に記載の卵黄粘度判定方法。
【請求項４】
前記殻付卵に光を照射する工程と、前記殻付卵の透過光を受光する工程と、前記透過光のスペクトルを取得する工程とは、
殻付卵の長軸を中心に殻付卵を回転させ、殻付卵の角度を変えて複数回行う、
請求項１乃至請求項３のいずれか１項に記載の卵黄粘度判定方法。
【請求項５】
加熱した殻付卵の卵黄粘度判定装置であって、
殻付卵に光をあてる投光部と、前記投光部と対向して配置され殻付卵の透過光を受ける受光部と、前記透過光のスペクトルを取得して該スペクトルに基づいて卵黄粘度を判定する判定部とを備え、
前記投光部は、光が殻付卵の膨らみ部分を通過するように配置され、
前記判定部は、スペクトルを多変量解析する手段を有していて、スペクトルの可視領域を含む波長に基づいて卵黄粘度を判定し、該多変量解析がＰＬＳ回帰分析であって、前記スペクトルの可視領域を含む波長を用い、該波長が６００ｎｍから９００ｎｍの領域を含む波長であり、波長の吸光度を抽出したデータと、見かけの卵黄粘度ηを常用対数変換したデータｌｏｇ（η）とに基づくＰＬＳ回帰分析である、
卵黄粘度判定装置。
【請求項６】
前記判定部は、
前記スペクトルの二次微分値を算出する手段を含む、
請求項５に記載の卵黄粘度判定装置。
【請求項７】
前記判定部は、
前記スペクトルを平滑化処理及びＳＮＶ処理する手段を含む、
請求項５又は請求項６に記載の卵黄粘度判定装置。
【請求項８】
殻付卵を長軸を中心に回転させる回転機構と、
殻付卵を搬送する搬送機構と、
を更に備える、
請求項５乃至請求項７のいずれか１項に記載の卵黄粘度判定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、殻付卵の判定方法及び判定装置に関し、特に、殻付卵の卵黄粘度を非破壊により判定する判定方法及び判定装置に関する。
【背景技術】
【０００２】
温泉卵などの加熱調理された殻付卵は、その卵黄粘度の違いにより味や食感が大きく異なる。そのため、事業者等には各々の用途に適した卵黄粘度を有する殻付卵を製造する需要がある。しかしながら、加熱調理された殻付卵の卵黄粘度は、卵自体の大きさや卵黄の大きさにバラつきがあるため、同一条件で加熱調理しても製造された卵の卵黄粘度にバラつきが生じる。
【０００３】
従来、殻付卵の卵黄粘度は、殻付卵を割卵して測定していた。すなわち、製造品から１０００個中１個など、任意抽出した殻付卵を割卵し卵黄膜を破膜させ、卵黄が一定時間に流れる距離を測定することにより卵黄粘度を測定していた。
【０００４】
一方、卵黄粘度とは異なるが、殻付卵を非破壊によって血卵や破卵を検出する方法及び装置が開示されている（例えば、特許文献１乃至特許文献３参照）。特許文献１の血卵の非破壊検出方法は、透過光のスペクトルにより血卵を検出するものである。特許文献２の血卵や破卵を検出する検査方法及び検査装置は、殻付卵の透過光を撮像し、画像解析を用いて殻付卵の血卵や破卵を検出するものである。特許文献３の非破壊検卵装置も、透過光を撮像し、画像解析によって検卵するものである。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００６−１６２４５４号公報
【特許文献２】特開平１１−０５６１５９号公報
【特許文献３】特開２００３−０６５９６１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
従来の卵黄粘度の測定は破壊を伴うものであるため、商品が損失する、時間がかかる、全ての製造品の粘度のバラつきについては言及できない、という課題がある。
【０００７】
また、特許文献１乃至特許文献３の方法又は装置は、血卵や破卵など異物の検出は可能であるが、卵黄粘度を判定することができない。すなわち、これらの方法又は装置は、正常卵と異常卵の○×を検出するものであり、殻付卵の卵黄粘度等の卵の状態を判定することはできない。
【０００８】
本発明は、上記課題に鑑みなされたものであって、非破壊による殻付卵の卵黄粘度判定方法及び判定装置を提供することを目的とする。
【課題を解決するための手段】
【０００９】
本発明の一の態様に係る卵黄粘度判定方法は、殻付卵の卵黄粘度判定方法であって、殻付卵に光を照射する工程と、殻付卵の透過光を受光する工程と、前記透過光のスペクトルを取得する工程と、前記スペクトルに基づき卵黄粘度を判定する工程と、を備える。この構成によれば、殻付卵の卵黄粘度を非破壊により判定することができる。
【００１０】
好ましくは、この卵黄粘度判定方法は、前記スペクトルに基づき卵黄粘度を判定する工程が、前記スペクトルを多変量解析する工程を含む。
【００１１】
また、この卵黄粘度判定方法は、前記多変量解析が回帰分析である。回帰分析には、多重線形回帰分析、主成分回帰分析、ＰＬＳ回帰分析などがある。
【００１２】
好ましくは、この卵黄粘度判定方法は、前記スペクトルに基づき卵黄粘度を判定する工程が、前記スペクトルの可視領域を含む波長を用いる。可視領域の波長に、卵黄粘度による影響が現れやすいためである。ここで、可視領域を含む波長とは、可視領域の波長の一部又は全部を含む波長をいう。
【００１３】
好ましくは、この卵黄粘度判定方法は、前記スペクトルの可視領域を含む波長が、５８８ｎｍから１０８４ｎｍの一部または全部の領域を含む波長である。この領域のスペクトルに基づくデータに、卵黄粘度の影響が現れやすいためである。
【００１４】
より好ましくは、この卵黄粘度判定方法は、前記スペクトルの可視領域を含む波長が、６００ｎｍ−９００ｎｍの一部または全部の領域を含む波長である。この領域のスペクトルに基づくデータに、卵黄粘度の影響がより現れやすいためである。
【００１５】
また、この卵黄粘度判定方法は、前記スペクトルに基づき卵黄粘度を判定する工程が、前記スペクトルの二次微分値を算出する工程を含んでもよい。スペクトルの二次微分値を用いても、卵黄粘度による違いが現れるためである。
【００１６】
また、この卵黄粘度判定方法は、前記スペクトルに基づき卵黄粘度を判定する工程が、前記スペクトルを平滑化処理及びＳＮＶ処理する工程を含んでもよい。スペクトルを平滑化処理及びＳＮＶ処理した値を用いても、卵黄粘度による違いが現れるためである。
【００１７】
また、この卵黄粘度判定方法は、前記殻付卵に光を照射する工程と、前記殻付卵の透過光を受光する工程と、前記透過光のスペクトルを取得する工程と、が殻付卵を回転させ、殻付卵の角度を変えて複数回行ってもよい。殻付卵の角度を変えて複数のスペクトルを取得し、これらの平均スペクトルを用いる。
【００１８】
好ましくは、この卵黄粘度判定方法は、前記スペクトルに基づき卵黄粘度を判定する工程が、あらかじめ実測した殻付卵サンプルの実測卵黄粘度データと、予測卵黄粘度データと、の相関に基づいて、判定対象の殻付卵の卵黄粘度を判定する。ここで、実測した卵黄粘度とは、割卵した卵黄を粘度計等により実測する。実測卵黄粘度データには、卵黄粘度の実測値に基づいて算出されたデータが含まれる。
【００１９】
好ましくは、この卵黄粘度判定方法は、前記実測卵黄粘度データが、実測した卵黄粘度を対数変換したデータであり、前記予測卵黄粘度データは、予測される卵黄粘度を対数変換したデータである。対数変換には例えば常用対数変換を用いる。対数変換することにより、数百から数百万（１０2から１０6）ｍＰａ・ｓまで変化する卵黄粘度の範囲を小さくして、誤差の範囲を小さく、計算をしやすくするためである。これにより、数百ｍＰａ・ｓ程度の低粘度から百万ｍＰａ・ｓの高粘度までを推定することができる。
【００２０】
本発明の一の態様に係る卵黄粘度判定装置は、殻付卵の卵黄粘度の判定装置であって、殻付卵に光をあてる投光部と、前記投光部と対向して配置され、殻付卵の透過光を受ける受光部と、前記透過光のスペクトルを取得して、該スペクトルを用いて卵黄粘度を判定する判定部と、を備える。この構成によれば、殻付卵の卵黄粘度を非破壊により判定することができる。
【００２１】
好ましくは、この卵黄粘度判定装置は、前記判定部が、前記判定部がスペクトルを多変量解析する手段を有する。
【００２２】
好ましくは、この卵黄粘度判定装置は、前記多変量解析が回帰分析である。
【００２３】
好ましくは、この卵黄粘度判定装置は、スペクトルの可視領域を含む波長に基づいて卵黄粘度を判定する。
【００２４】
好ましくは、この卵黄粘度判定装置は、前記スペクトルの可視領域を含む波長が、５８８ｎｍから１０８４ｎｍの一部または全部の領域を含む波長である。
【００２５】
より好ましくは、この卵黄粘度判定装置は、前記スペクトルの可視領域を含む波長が、６００ｎｍから９００ｎｍの一部または全部の領域を含む波長である。
【００２６】
また、この卵黄粘度判定装置は、前記判定部が、前記スペクトルの二次微分値を算出する手段を含んでもよい。
【００２７】
また、この卵黄粘度判定装置は、前記判定部が、前記スペクトルを平滑化処理及びＳＮＶ処理する手段を含んでもよい。
【００２８】
好ましくは、この卵黄粘度判定装置は、殻付卵を回転させる回転機構と、殻付卵を搬送する搬送機構と、を更に備える。
【００２９】
好ましくは、この卵黄粘度判定装置は、前記判定部が、記憶部を有し、前記記憶部には、あらかじめ実測した数種の粘度の卵黄を具備する殻付卵サンプルの実測卵黄粘度データと、殻付卵サンプルのスペクトルに基づく予測卵黄粘度データと、の相関データが保存されている。
【００３０】
好ましくは、この卵黄粘度判定装置は、前記実測卵黄粘度データが、実測した卵黄粘度を対数変換したデータであり、前記予測卵黄粘度データが、スペクトルに基づく卵黄粘度を対数変換したデータである。
【発明の効果】
【００３１】
本発明の卵黄粘度方法及び卵黄粘度装置によれば、殻付卵の卵黄粘度を判定することができる。
【図面の簡単な説明】
【００３２】
【図１】本発明の一実施形態に係る判定装置の概略構成を示す図である。
【図２】本発明の一実施形態に係る判定装置の投光部及び受光部を示す図である。
【図３】本発明の一実施形態に係る判定装置の搬送部を示す図である。
【図４Ａ】本発明の一実施形態に係る判定方法に係るフロー１を示す図である。
【図４Ｂ】本発明の一実施形態に係る判定方法に係るフロー２を示す図である。
【図５】殻付卵サンプルを割卵した状態を示す図である。
【図６】殻付卵サンプルのスペクトルを示す図である。
【図７Ａ】殻付卵サンプルのスペクトルの二次微分値と卵黄粘度の実測値（対数変換後）との相関を示す図である。
【図７Ｂ】殻付卵サンプルの卵黄粘度の実測値（対数変換後）と予測値（対数変換後）を示す図である。
【図８Ａ】変形例の殻付卵サンプルを割卵した状態を示す図である。
【図８Ｂ】変形例の殻付卵サンプルのスペクトルを示す図である。
【図８Ｃ】変形例の殻付卵サンプルの実測値と予測値の相関を示す図である。
【図９】ＰＬＳ回帰分析結果を示す図である。
【図１０】ＰＬＳ回帰分析における回帰ベクトルを示す図である。
【図１１】マーテンスの不確かさ検定前後の回帰ベクトルを示す図である。
【図１２】スペクトル測定器、セルを示す図である。
【図１３】卵黄と卵白をそれぞれ加熱した場合の各波長とスペクトルとの関係を示す図である。
【図１４】平滑化及びＳＮＶ処理を用いた殻付卵サンプルの卵黄粘度の実測値と予測値を示す図である。
【図１５】平滑化処理とＳＮＶ処理の数式である。
【図１６】殻付卵サンプルのビタミンＥ含有量の実測値と予測値を示す図である。
【図１７】従来の落球式粘度計を示す図である。
【発明を実施するための形態】
【００３３】
以下、本発明に係る一実施形態を図面に基づき説明するが、本発明は下記実施形態に限定されるものではない。
【００３４】
＜１．判定装置の構成＞
図１は、本実施形態の判定装置の概略構成を示す図である。図２は、同判定装置の投光部及び受光部を示す図である。図３は、判定装置の搬送部を示す図である。図１に示すように、本実施形態の判定装置１０は、投光器２０と受光器３０と判定部４０とを主要な構成として備える。
【００３５】
投光器２０は、判定対象の殻付卵Ｘに光を照射するためのものである。投光器２０は、投光部２１と光源２２とを有する。光源２２は、種々の波長を発するものであって、少なくとも可視領域の光を発するものであれば特に限定されない。なお、本実施形態では、一例として５８８ｎｍ−１０９２ｎｍのものを用いている。
【００３６】
また、投光器２０は、光ファイバなどを用いた細い光を照射するものが望ましい。広い光を用いた場合の、殻の表面に反射した反射光の影響を抑えるためである。また、投光器２０は、投光部２１と光源２２とが一体のものでもよく、投光部２１と光源２２とが光ファイバ等で接続されたものでもよい。
【００３７】
受光器３０は、殻付卵Ｘを透過した透過光を受光するためのものである。受光器３０は、受光部３１と変換部３２とを有する。受光部３１は、投光部２１と対向配置されていて、殻付卵Ｘの透過光を受光する（図２参照）。変換部３２は、受光部３１が受光した光を、図示しない光電変換素子を介して、各波長ごとにデジタル信号に変換するものである。これにより、殻付卵Ｘの透過光のスペクトルが得られる。
【００３８】
判定部４０は、透過光のスペクトルを演算処理し、該スペクトルから殻付卵の卵黄粘度を数値化する。卵黄粘度の数値化については後述する。判定部４０は、演算部４１と記憶部４２と表示部４３とを有し、受光器３０と接続されている。なお、判定部４０は、判定装置専用のものでもよいし、汎用のコンピュータでもよい。また、表示部４３は必須の構成ではなく、例えば結果が印刷される構成などとしてもよい。
【００３９】
本実施形態の判定装置１０は、回転機構５０と搬送機構６０とを更に備える。回転機構５０は、殻付卵を回転させる機構である。本実施形態の回転機構５０は、一例として回転軸５１と一対のアーム５２とを具備する機構である。アーム５２の先端には、ゴム等からなる弾性部材５３が取り付けられていて、殻付卵を両側から保持する。アーム５２は、開閉可能に回転軸５１に接続されている。なお、本実施形態のアーム５２は一対に構成されているが、アーム５２の数は３でも４でもよい。
【００４０】
回転軸５１は、図示しないモータに接続されている。アーム５２が殻付卵を保持した状態で回転軸５１が所定角度回転して、投光部２１及び受光部３１に対する殻付卵の角度を変える。なお、本実施形態の判定装置１０においては、殻付卵を回転させる構成であるが、これに限られず、投光部２１及び受光部３１を回転させる構成としてもよい。
【００４１】
搬送機構６０は、殻付卵を搬送する機構である。本実施形態の搬送機構６０は、一例としてベルトコンベアであるが、これに限られない。図３に示すように、搬送機構６０により殻付卵が搬送され、投光部２１と受光部３１の間で停止して、殻付卵の判定が行われる。判定が終わった殻付卵は搬送され、次の殻付卵が投光部２１と受光部３１の間まで搬送される。これが順に繰り返され、複数の殻付卵が判定される。なお、複数の殻付卵の全量判定が望ましいが、所定数に１つの割合で殻付卵を判定することも可能である。
【００４２】
＜２．判定方法＞
次に、本発明に係る一実施形態である卵黄粘度の判定方法を図面に基づき説明する。一例として、判定装置１０を使用する。図４Ａ・Ｂは、本実施形態の判定方法に係るフローチャートである。
【００４３】
［判定前準備］
（１）判定前準備として、加熱レベルが６段階に異なる殻付卵サンプル１−７を用意する。図５は殻付卵サンプル１−７を割卵した図である。殻付卵サンプル１は生卵であり、殻付卵サンプル２−７は順に加熱レベルを高くした殻付卵である。殻付卵サンプル１が最も卵黄粘度が低く、殻付卵サンプル７が最も卵黄粘度が高い。なお、ここでいう加熱レベルとは、加熱時間の長短をいう。
【００４４】
（２−１）殻付卵サンプル１−７について、判定装置１０を用いて各スペクトルを収集する。具体的には、殻付卵サンプル１を判定装置１０にセットして、搬送機構６０により投光器２０と受光器３０の間まで搬送し、停止する。
【００４５】
（２−２）投光器２０の投光部２１より殻付卵サンプル１に光を照射し、殻付卵サンプル１を通過した透過光を受光器３０の受光部３１が受光する。
【００４６】
（２−３）受光部３１が受光した透過光を、光電変換素子を介して、変換部３２が各波長ごとにデジタル信号に変換し、殻付卵サンプル１の透過光のスペクトルを得る。そして、スペクトルが判定部４０の記憶部４２に記憶される。
【００４７】
（２−４）続き、回転機構５０が殻付卵サンプル１を略１２０度回転させ、上記工程（２−２）及び（２−３）を行い、スペクトルが記憶部４２に記憶される。更に、再度回転機構５０が殻付卵サンプル１を略１２０度回転させ、スペクトルが記憶部４２に記憶される。
【００４８】
（２−５）殻付卵サンプル２−７についても、同様にスペクトルの取得を行う。以上の工程により、殻付卵サンプル１−７のスペクトルが得られる（図６参照）。
【００４９】
（３−１）次に、殻付卵サンプル１−７をそれぞれ割卵して、落球式粘度計１００を用いて卵黄粘度を算出する。具体的には、殻付卵サンプル１の卵黄を落球式粘度計１００の測定管１０１に入れて、本体１０２に取り付ける（図１７参照）。そして、上部に錘が取り付けられたニードル１０３が、上方のスタート位置から一定距離落下するまでの時間ｔを計測する。
【００５０】
（３−２）この時間ｔを用いて数式１により粘度を算出する。以下、これを見かけの卵黄粘度ηと称し、必要に応じて見かけの粘度ηを常用対数変換したｌｏｇ（η）を用いる。また、殻付卵サンプル２−７についても、同様に見かけの卵黄粘度ηを算出する。なお、下記式において、ｐｓ＝４７．１４４４ｇ／ｃｍ3、ｐｆ＝１．０００ｇ／ｃｍ3を用いる。
【００５１】
【数１】
0006753651.tif 000002
【００５２】
（４−１）次に、得られた殻付卵サンプル１−７のスペクトルと見かけの卵黄粘度ηとに基づき、多変量解析を用いてモデル化する。また、本実施形態においては、多変量解析の一例としてＰＬＳ回帰分析を用いる。なお、ＰＬＳ回帰分析にかえて、多重線形回帰分析、主成分回帰分析を用いてもよい。
【００５３】
（４−２）殻付卵サンプル１−７のスペクトルは、上記工程により判定部４０の記憶部４２に保存されているので、見かけの卵黄粘度ηを入力して記憶部４２に保存する。
【００５４】
（４−３）演算部４１は、スペクトルと見かけの卵黄粘度ηとを用いてＰＬＳ回帰分析を行う。具体的には、演算部４１は、スペクトルを二次微分して６００ｎｍ―７２０ｎｍにおける波長の吸光度を抽出したデータと、見かけの卵黄粘度ηを常用対数変換したデータｌｏｇ（η）とに基づき、ＰＬＳ回帰分析を行う。以下、必要に応じて、見かけの卵黄粘度ηを対数変換したデータｌｏｇ（η）を、見かけの卵黄粘度ｌｏｇ（η）や、ｌｏｇ（η）実測値と称する場合がある。なお、説明変量をスペクトルの二次微分値とし、目的変量を見かけの卵黄粘度ｌｏｇ（η）とする。
【００５５】
（４−４）これにより、見かけの卵黄粘度ｌｏｇ（η）とスペクトルの二次微分値との相関が得られる（図７Ａ）。この結果に基づき、見かけの卵黄粘度ｌｏｇ（η）と予測の卵黄粘度ｌｏｇ（η）との相関を示す回帰直線を得る。なお、本実施形態の回帰直線は、見かけの卵黄粘度ｌｏｇ（η）と予測の卵黄粘度ｌｏｇ（η）とが略１：１に対応する（図７Ｂ参照）。
【００５６】
・上記において、殻付卵サンプルを７個としているが、これに限られず、サンプル数を少なくしてもよいし、サンプル数を多くしてもよい。なお、図８は、殻付卵サンプルを６個にした場合の結果である。図８Ａは殻付卵サンプルを割卵した状態を示す図である。図８Ｂは、これら６個の殻付卵サンプルのスペクトルを示す図である。図８Ｃは、これら６個の殻付卵サンプルの実測値と予測値の相関を示す図である。サンプル数を多くすればより精度の高い卵黄粘度の判定が可能となる。
【００５７】
・上記（２−４）において、殻付卵サンプルのスペクトルを１２０度毎に３回測定しているが、１回や２回でもよく、４回以上としてもよい。
【００５８】
・上記（４−１）乃至（４−４）において、ＰＬＳ回帰分析を判定装置１０の判定部４０が行っているが、これに限られず、別途用意したコンピュータ等により行ってもよい。また、判定装置１０に測定結果やＰＬＳ回帰分析の結果等を表示するディスプレイや、結果を印刷する機構を設けてもよい。
【００５９】
・上記（４−３）において、スペクトルにおける６００ｎｍ−７２０ｎｍの波長領域の吸光度を抽出しているが、これに限られず、例えば５８８ｎｍ−１０８４ｎｍ、５９８ｎｍ−９４６ｎｍ、６１０ｎｍ−９４６ｎｍ、６００ｎｍ−９００ｎｍ、６１０ｎｍ−７２０ｎｍ、５９８ｎｍ−７０６ｎｍなどとしてもよい。すなわち、可視領域の波長を含むものであればよい（図９、図１０、図１１参照）。
【００６０】
・上記（４−３）において、スペクトルを二次微分したデータを用いているが、これにかえて、例えばスペクトルをそのまま用いることや、スペクトルを平滑化及びＳＮＶ処理したデータを用いることもできる。図１３に示すように、スペクトルを平滑化及びＳＮＶ処理したデータによっても、卵黄粘度が予測可能である。なお、平滑化処理及びＳＮＶ処理は、例えば図１４の式を用いる。
【００６１】
［殻付卵の判定］
（１）調査対象の殻付卵Ｘを判定装置１０にセットして、搬送機構６０により投光器２０と受光器３０の間まで搬送し、停止する。
【００６２】
（２）投光器２０の投光部２１より、近赤外線を殻付卵Ｘに照射し、殻付卵Ｘの透過光を受光器３０の受光部３１が受光する。
【００６３】
（３）受光部３１が受光した透過光を、光電変換素子を介して、変換部３２が各波長ごとにデジタル信号に変換し、殻付卵Ｘの透過光のスペクトルを得る。
【００６４】
（４）透過光のスペクトルが判定部４０の記憶部４２に記憶され、演算部４１がスペクトルにおける６００ｎｍ−７２０ｎｍの二次微分値に基づいて、殻付卵Ｘの予測の卵黄粘度ｌｏｇ（η）を算出する。なお、必要に応じて、殻付卵Ｘの予測の卵黄粘度ｌｏｇ（η）を、ｌｏｇ（η）予測値と称する場合がある。
【００６５】
（５）以上の工程により得られた殻付卵Ｘの予測の卵黄粘度ｌｏｇ（η）を殻付卵サンプル１−７の卵黄粘度と照らし合わせれば、殻付卵Ｘの卵黄の状態が判定できる。例えば、卵黄粘度ｌｏｇ（η）＝３．０であれば殻付卵サンプル３の状態に近似する。
【００６６】
［検証］
上記工程により、殻付卵の卵黄粘度を判定できることがわかったが、念のため発明者らは、卵黄粘度に好適な波長の範囲と、卵黄粘度の卵白のゲル化による影響について検証した。なお、本実施形態における好適な波長の範囲の検証には、一例としてマーテンスの不確かさ検定手法を用いた。
【００６７】
（１）まず、卵黄粘度に好適な波長の範囲の検証について、本実施形態におけるＰＬＳ回帰分析で得られた回帰ベクトル（図１１の上グラフ参照）から、重要と思われる波長を選定し、選定した波長のみを用いて再度ＰＬＳ回帰分析を行った。これにより得られた回帰ベクトルを図１１に示す（図１１の下グラフ参照）。この結果から、６００ｎｍ−９００ｎｍにおける波長領域のスペクトルが、卵黄粘度による影響を反映していることがわかった。すなわち、卵黄粘度の判定においては、９００ｎｍを超える波長については、除外することが好ましい。
【００６８】
（２−１）次に、卵白のゲル化の影響について検証した。図１２は、スペクトル測定器、セル、セルに入れた卵黄及び卵白の状態を示す図である。図１２に示すように、生卵を割卵し、卵黄と卵白に分けてそれぞれをセルに入れる。各セルをスペクトル測定器に入れ、卵黄と卵白を各セル内で加熱しつつ、加熱時間を０分から２分間隔で３０分間、透過スペクトルを採取する。
【００６９】
（２−２）これにより得られた卵黄と卵白それぞれのｔ分後−０分後のスペクトルの二次微分値の差スペクトルと、スペクトルの各波長との相関を得る（図１３参照）。図１３に示すように、卵白の加熱程度による変動が大きい波長領域が、上記（１）のマーテンスの不確かさ検定で得られた卵黄粘度で除外された波長領域（９００ｎｍを超える波長（９００−９７０ｎｍ付近）と略一致することがわかった。
【００７０】
（２−３）この結果から、６００ｎｍ−９００ｎｍ付近の波長領域は卵黄粘度による影響が支配的であり、９００ｎｍ−９７０ｎｍ付近や９００ｎｍ−１０００ｎｍ付近の波長領域は卵白のゲル化による影響が支配的であることがわかった。すなわち、卵黄と卵白それぞれが支配的である波長領域が異なることから、上記した卵黄粘度の判定において、卵白のゲル化の影響を受けないことが検証によりわかった。一方、卵白のゲル化についても、上記した卵黄粘度の判定と同様の方法を用いて、選定するスペクトルの波長領域を変更することで、卵黄粘度の影響を受けずに判定することが可能である。
【００７１】
［３．本実施形態の効果］
本実施形態の判定装置及び判定方法によれば、殻付卵の卵黄粘度を非破壊により判定することができる。また、多変量回帰分析を利用することにより、卵白のゲル化の程度に影響されることなく卵黄粘度を推定可能であり、併せて卵黄粘度を常用対数変換したデータを用いることで、幅広い卵黄粘度の推定に対応可能である。また、本実施形態では主としてスペクトルを二次微分したデータを用いているが、スペクトルをそのまま用いることもでき、または平滑化及びＳＮＶ処理したデータを用いることもできる（図９参照）。
【００７２】
＜その他の実施形態＞
以上のとおり、図面を参照しながら本発明の好適な実施形態を説明したが、本発明の趣旨を逸脱しない範囲内で、種々の追加、変更又は削除が可能である。したがって、そのようなものも本発明の範囲内に含まれる。
【００７３】
また、本発明の判定方法及び判定装置は、卵黄粘度以外にも適用可能である。実測、予測、判定する対象を卵黄粘度にかえて、卵白のゲル化度やビタミンＥ含有量とすることが可能である。卵白のゲル化度については、上述した通り卵黄粘度と同様の方法又は装置により、選定するスペクトルの波長領域を例えば９００ｎｍ−９７０ｎｍ付近や９００ｎｍ−１０００ｎｍ付近とすることで判定が可能である。
【００７４】
また、図１６に示すように、ビタミンＥの含有量についても予測が可能であることがわかる。図１６は、本実施形態の卵黄粘度と同様の方法又は装置により、ビタミンＥ含有量の実測値と予測値の相関を現した図である。なお、ビタミンＥの含有量の実測については、既知の方法により測定している。
【符号の説明】
【００７５】
１０  判定装置
２０  投光器
２１  投光部
２２  光源
３０  受光器
３１  受光部
３２  変換部
４０  判定部
４１  演算部
４２  記憶部
４３  表示部
５０  回転機構
５１  回転軸
５２  アーム
５３  弾性部材
６０  搬送機構
１００  落球式粘度計
１０１  測定管
１０２  本体
１０３  ニードル
Ｘ  殻付卵
【図１】
0006753651.tif 000003
【図２】
0006753651.tif 000004
【図３】
0006753651.tif 000005
【図４Ａ】
0006753651.tif 000006
【図４Ｂ】
0006753651.tif 000007
【図５】
0006753651.tif 000008
【図６】
0006753651.tif 000009
【図７Ａ】
0006753651.tif 000010
【図７Ｂ】
0006753651.tif 000011
【図８Ａ】
0006753651.tif 000012
【図８Ｂ】
0006753651.tif 000013
【図８Ｃ】
0006753651.tif 000014
【図９】
0006753651.tif 000015
【図１０】
0006753651.tif 000016
【図１１】
0006753651.tif 000017
【図１２】
0006753651.tif 000018
【図１３】
0006753651.tif 000019
【図１４】
0006753651.tif 000020
【図１５】
0006753651.tif 000021
【図１６】
0006753651.tif 000022
【図１７】
0006753651.tif 000023