JP6822691B2

JP6822691B2 - 種卵検査システム及び種卵検査プログラム

Info

Publication number: JP6822691B2
Application number: JP2018553805A
Authority: JP
Inventors: 伸一藤谷
Original assignee: NABEL CO.,LTD.
Current assignee: NABEL CO.,LTD.
Priority date: 2016-11-30
Filing date: 2017-11-22
Publication date: 2021-01-27
Anticipated expiration: 2037-11-22
Also published as: CN109862783A; CN109862783B; WO2018101139A1; JPWO2018101139A1

Description

本発明は、種卵の情報を取得して、その種卵を選別する種卵検査システムに関するものである。

食用卵の生産用に養鶏場で飼養される鶏の系統を採卵鶏（レイヤー）と呼ぶ。食用卵の生産者は、採卵用の実用鶏の親となる種鶏群を飼養する種鶏会社から実用鶏の雛を購入して、その雛を大雛まで育成して食用卵の生産に使用している。

種鶏会社において、実用鶏の親の種鶏群が飼養されている農場を種鶏農場と言う。この種鶏農場で産まれた種卵は、図１に示すように、一旦貯卵場所に集められ、その後に予備加温と呼ばれる工程を経て、孵卵場のセッターに入卵されて孵卵が開始される。セッターに入卵してから１８日目又は１９日目に孵卵場のハッチャーに移卵されて、入卵から約２１日目に雛が孵化する。

この孵卵工程では、種卵をセッターからハッチャーに移卵する際、又はその他の孵卵途中において、生存卵と非生存卵（内部の胚が死亡した種卵または未受精卵）とに分ける作業が行われている。

生存卵と非生存卵とを高速かつ非侵襲的に選別するものとして、特許文献１に示す種卵検査装置が考えられている。この種卵検査装置は、入卵から所定日数経過した時点（例えば、１８日目）で種卵に光を照射し、その時点で得られた透過光の光強度の時系列データ中の変動成分の有無と、透過光の光強度が適切な範囲に入るように個別の卵毎にＬＥＤ光量を制御したときの光源光量の制御量を用いて生存卵と非生存卵とを判定するものである。
このような移卵時（１８日目）に検査するのと同様な方法により、中間時（例えば１１日目や１４日目）に検査して、生存卵と非生存卵とを判定する場合もある。これは、検査時期を１８日目よりも早めることにより、（１）移卵作業時点での腐敗卵の発生の可能性を減らす、（２）そのロット中の未受精卵や中止卵の数量を知り、おおよその雛の歩留まりを予測する、などの利用価値がある。

特許第４８５８８６３号公報

このように、従来の種卵検査装置は、生存卵と非生存卵とを判別したいという主目的のもと、移卵時または中間時の一時点での情報に基づいて判別するものであったが、本発明では、生存卵と非生存卵との判別以外にも、様々な観点から高速かつ非侵襲的に種卵を選別することを課題とする。

具体的に言えば、単一の日での測定だけではわからない変化率の把握や、計測データの卵座に対する依存性の解消、また、複数の証拠の利用による性別判定精度の向上を図ることができる種卵検査装置を提供する。また、データを測定する日とトレイ間で卵を移動させる選別日とを分けて、種卵に特有な不安定な時期を避けて卵の仕分けが可能となる種卵検査装置を提供する。

すなわち本発明に係る種卵検査システムは、種卵の状態を示す状態情報を計測する計測部と、前記計測部により得られた複数日の状態情報を個別の種卵毎に関連づけて記録する記録部と、前記記録部に記録された前記複数日の状態情報の個別の種卵毎に関連づけられた記録に基づいて前記種卵を選別する種卵選別部とを備えることを特徴とする。

ここで、「状態情報」とは、種卵の状態を示すものであればどのようなものであってもよく、例えば、胚や、血管又は血液の状態を示すもの、卵重、卵の長径並びに短径や、卵の体積など卵の大きさを示すものや卵殻色などが挙げられる。

また、「複数日」とは、種卵が孵化するまでの２時点以上を指し、孵卵器に入った後（入卵後）のみならず、それらのうち少なくとも１日は孵卵器に入る前（入卵前）でもかまわない。

このような種卵検査システムであれば、異なる複数日における状態情報を用いて様々な観点から種卵を選別することができる。例えば、異なる複数日における状態情報から胚の成長速度などの時間的変化を把握することができる。この時間的変化は、雛の性別、孵化時間、雛の増体性などの属性、胚の発育の中止、中止時期などに関係を持つ。

以下、具体的に説明すれば、本願発明者は、胚が生存して発育が進んでいる種卵の内部では、胚の成長並びに血管の生成及び伸長が時々刻々と進行していることに着目し、種卵の内部の胚、血管や血液の状態の時間的変化を利用することで、生存卵と非生存卵との判別以外にも、様々な観点から種卵を分離できるのではないかと考えた。

例えば、本願発明者は、種卵の内部の血管の成長度合い（成長速度）にオスの胚とメスの胚とで有意な差があるのではないかという仮説を立てた。
そして、本願発明者は、血管の成長度合いが血管の成長に伴う血液の量から見積もれること、つまり、血管の成長度合いが血液中のヘモグロビンの吸収のピーク波長の一つである５７８ｎｍにおける透過率から見積もれることに着目して、波長５７８ｎｍを有する光を種卵に照射して、波長５７８ｎｍにおける透過率の変化を計測した。

その結果、孵卵開始（入卵）から７日目の時点の波長５７８ｎｍにおける透過率において、確かにオスの胚とメスの胚とで有意な差異が存在することを見出した。なお、詳細については後述する。
つまり、本願発明者は、血管や血液の形成の進行の仕方に性差があり、血管の成長度合いからオスの胚とメスの胚が判別できるという新しい知見を得た。

また、本願発明者は、種卵の胚の成長度合い（成長速度）にオスの胚とメスの胚とで有意な差があるのではないかという仮説を立てた。
そして、本願発明者は、胚の成長度合いが血液中のヘモグロビンや水分の吸収を受けにくい波長８７０ｎｍにおける不透明度から見積もれることに着目して、波長８７０ｎｍを有する光を種卵に照射して、波長８７０ｎｍにおける不透明度の変化を計測した。

その結果、孵卵開始から７日目、８日目の時点の波長８７０ｎｍにおける不透明度において、確かにオスの胚とメスの胚とで有意な差異が存在することを見出した。なお、詳細については後述する。
つまり、本願発明者は、胚の形成の進行の仕方に性差があり、胚の成長度合いからオスの胚とメスの胚が判別できるという新しい知見を得た。

上記の新しい知見に基づいてなされた本発明によれば、複数日における種卵の内部の胚、血管や血液の状態を示す状態情報から得られる種卵の内部の胚の成長情報や血管の成長情報を用いて種卵を選別することができる。

また、複数日の状態情報を記録部に記録させているので、異なる日に計測された状態情報の管理を容易にすることができる。

さらに、孵卵初期の種卵は、胚が小さく外部からの振動等の刺激に鋭敏であり、この時期の卵を仕分けるために卵を別々のトレイに移し変える作業を行うと孵化率の低下を招く恐れがある。種卵を移し変える作業を行うときに種卵の選別を行う場合には、前記移し変える作業は孵化率に影響を与えない時期まで待つことが望ましい。本発明では、複数日の状態情報を記録部に記録させているので、孵化率に影響を与えない時期まで待った後に、種卵の選別及び移し変える作業を行うことができる。

具体的に種卵選別部は、前記複数日の状態情報により求まる前記種卵の内部の胚の成長状態に基づいて前記種卵を選別するものであることが望ましい。
このように胚の成長状態に基づいて種卵を選別することにより、雛の性別、孵化時間、雛の増体性などの属性、胚の発育の中止、中止時期の予測を行うことができる。

採卵鶏の場合は、オスの雛の経済的価値は乏しい。採卵鶏のオスの雛は、採卵の目的に使用できないのは勿論のこと、食肉用としても増体性が肉用鶏に比較して劣るからである。したがって、孵化したオスの雛は、雌雄鑑別の後に廃棄されており、オスの孵卵コストなど経済的な損失のみならず、生命倫理上の観点からも問題視されている。
この状況を考えれば、仮に孵卵中の種卵であっても、ワクチン卵（ワクチン製造に使用される種卵）に転用可能な日限である孵卵開始から９日目以前に孵化する雛の性別が高い確率で予測できれば、オスの孵化する確率の高い種卵をワクチン卵に転用することが可能であり、経済的な損失を少なくし、新たな収益を確保できるとともに、生命倫理上の問題も緩和できる。
このため、前記種卵選別部は、前記複数日の状態情報により予想される前記種卵から孵化する雛の性別に基づいて前記種卵を選別するものであることが望ましい。
この構成であれば、オスの雛が孵化する種卵の割合が高い種卵の集まりを選抜して、ワクチン卵として活用又は販売することができる。さらに、限られた孵卵場の収容能力をメスの雛が孵化する種卵に優先的に割り当てることができる。その上、オスの雛を殺処分する生命倫理上の問題を緩和することができる。

また、本発明に係る種卵検査プログラムは、種卵の状態を示す状態情報を計測する計測部により得られた前記複数日の状態情報を個別の種卵毎に関連づけて記録する記録部と、前記記録部に記録された前記複数日の状態情報の個別の種卵毎に関連づけられた記録に基づいて前記種卵を選別する種卵選別部と、としての機能をコンピュータに備えさせることを特徴とする。

このように構成した本発明によれば、様々な観点から高速かつ非侵襲的に種卵を選別することができる。

１００・・・種卵検査システム
２００・・・セッタートレイ
２・・・光照射部
３・・・光検出部
４・・・記録部
５・・・種卵選別部

一般的な孵卵場における雛の生産工程を示すフロー図である。検証実験に用いた種卵の非破壊検査装置の構成を模式的に示す図である。孵卵０日目から８日目までの相対透過率スペクトルを示す図である。特定波長毎の相対透過率の時間変化を示す図である。種卵の雌雄別における孵卵７日目の５７８ｎｍの相対透過率を示す図である。種卵の雌雄別における初期卵重と孵卵７日目の５７８ｎｍの相対透過率との関係を示す図である。種卵の雌雄別における孵卵０日目の５７８ｎｍの相対透過率に対する孵卵７日目の５７８ｎｍの相対透過率の比を示す図である。種卵の雌雄別における初期卵重と孵卵０日目の５７８ｎｍの相対透過率に対する孵卵７日目の５７８ｎｍの相対透過率の比との関係を示す図である。種卵の雌雄及び未受精卵における孵卵６日目の５７８ｎｍの相対透過率を示す図である。種卵の雌雄及び未受精卵における孵卵０日目の５７８ｎｍの相対透過率に対する孵卵６日目の相対透過率の比を示す図である。初期卵重と孵卵０日目の相対透過率に対する孵卵７日目の相対透過率の比とを説明変数とし、性別（オス＝１、メス＝０）を目的変数としたロジスティック回帰分布の判別値の分布を示す図である。本実施形態の種卵の非破壊検査装置の構成（ヘッド退避位置Ｎ）を模式的に示す図である。同実施形態の種卵の非破壊検査装置の構成（ヘッド計測位置Ｍ）を模式的に示す図である。本実施形態の非破壊検査装置に用いるセッタートレイの平面図、及びセッタートレイの卵座に種卵を載置した状態の平面図である。同実施形態の複数の種卵を同時に計測している状態を示す図である。種卵の雌雄別の孵卵０日目の５７４ｎｍの相対透過率に対する孵卵７日目の５７４ｎｍの相対透過率の比を示す図である。種卵の雌雄別の初期卵重と孵卵０日目の５７４ｎｍの相対透過率に対する孵卵７日目の５７４ｎｍの相対透過率の比との関係を示す図である。雌雄判別及び未受精卵の閾値を示す模式図である。ＲＦＩＤを設けたセッタートレイを用いた種卵検査システムを示す模式図である。変形実施形態の種卵の非破壊検査装置を組み込んだ種卵仕分けシステムを模式的に示す図である。変形実施形態の雛の生産工程を示すフロー図である。種卵の雌雄別における初期卵重とＤａｙ３に対するＤａｙ８の不透明度の比との関係を示す図である。種卵の雌雄別におけるＤａｙ３に対するＤａｙ８の不透明度の比を示す図である。種卵の雌雄別におけるＤａｙ３に対するＤａｙ８の不透明度の比とＤａｙ０に対するＤａｙ７の５７４ｎｍの相対透過率の比との関係を示す図である。Ｄａｙ３に対するＤａｙ８の不透明度の比と、Ｄａｙ０に対するＤａｙ７の相対透過率の比とを説明変数とし、性別（オス＝１、メス＝０）を目的変数としたロジスティック回帰分布の判別値の分布を示す図である。Ｄａｙ３に対するＤａｙＮ（Ｎ＝３〜１８）の不透明度の比の変化を示す成長曲線を示す図である。

＜検証実験＞
まず、本発明に至った仮説の１つに対する検証実験について説明する。
採卵鶏（ジュリア・ライト）の種卵を９０個入手した。なお、入手時に全ての種卵の卵重（以下、初期卵重という。）を測定した。また、本実験に用いた種卵は予備加温をしていない。

孵卵器に入卵する直前（以下、Ｄａｙ０という。）から４８時間以降２４時間毎にハロゲンランプ光源の光を種卵の側方から照射し、当該種卵からの透過光を分光器により分光して透過光の分光データを測定した。なお、測定装置を図２に示す。
これらのデータは、孵卵器に入れてから何日目のデータかを識別する目的で、Ｎ日目のデータならＤａｙＮのデータとして参照する。

測定環境の校正の目的で、測定前に長さ３０ｍｍ、直径４５ｍｍの円筒状の合成樹脂（ポリテトラフルオロエチレン）ブロックを種卵の代わりにブロックの底面が投光側に来るように測定台に置いて分光データを測定した。

５００ｎｍから９００ｎｍの波長範囲において、１ｎｍ刻みで、種卵の分光データを合成樹脂ブロックの分光データで割って、合成樹脂ブロックをリファレンスにした相対透過率のスペクトルを各卵毎に求めた。

相対透過率Ｔ（λ）
＝種卵の波長λの分光データ／合成樹脂ブロックの波長λの分光データ

図３に、種卵のＤａｙ０からＤａｙ８までの相対透過率スペクトルを示す。この図３に示すように、内部で胚が成長している種卵は、相対透過率スペクトルの波形が毎日変化する。

この変化を５７８ｎｍ、６２３ｎｍ、７５０ｎｍ、８１０ｎｍ、８７０ｎｍの波長毎に分けて示したのが図４である。この図４に示すように、全ての波長に対して、Ｄａｙ４からＤａｙ５にかけての相対透過率の変化が著しいが、特に両日の５７８ｎｍでの相対透過率の比で計算される減少率が他の波長に比較して著しく大きい。

この理由は、次の通りである。
Ｄａｙ３から卵黄の胚盤の周りに血管の生成が見られる。この血管が日々伸びることによって、血液の成分であるヘモグロビンの量が増える。そうすると、ヘモグロビンによる光の吸収のピーク波長の付近において、光の吸収が大きくなる。ヘモグロビンの吸収のピーク波長は、可視域では５７８ｎｍ、５４０ｎｍ、４１０ｎｍ付近であることが知られており５７８ｎｍは、その一つであるからである。

したがって、５７８ｎｍの相対透過率の時間変化は、血管の生成及び伸長に対応している。

次に、この血管の成長にオスの胚とメスの胚とで違いがあるかを調べた。
９０個の種卵の孵卵を継続して、雛を孵化させた。９０個の種卵から７６羽の雛が孵化した。雌雄鑑別の結果、３５羽がメスで４１羽がオスであった。なお、５個が未受精卵であった。残りは、胚が孵化前に死んだ発育中止卵であった。

図５は、雛が孵化した７６個の種卵において、Ｄａｙ７における５７８ｎｍの相対透過率を雛の性別に分けて図示したものである。図５中の白抜きのひし形のシンボルは、雛の雌雄鑑別の結果オスと鑑別された種卵を示すものであり、黒塗りの丸のシンボルはメスと鑑別された種卵を示すものである。以下、他の散布図における雌雄鑑別結果における２種のシンボルの使い分けはこれに準じている。
この図５から、メスの５７８ｎｍの相対透過率がオスのものより有意に小さいことがわかる（ｔ検定ｐ値＜０．００２）。
この図５において例えばＤａｙ７における５７８ｎｍの相対透過率の値が、０．００５より大きい種卵を抽出すれば、２６個が抽出される。これは全体７６個中の３４%である。２６個中１９個がオスであり、オスの比率は、７３％である。すなわち、特定時点（この場合はＤａｙ７）の５７８ｎｍの相対透過率で血管や血液の生成の程度を推定し、その値が定められた閾値より大きいものを抽出すれば、オスの孵化する確率の高い卵が抽出される。

図６は、横軸に初期卵重をとり、縦軸にＤａｙ７における５７８ｎｍの相対透過率をとった散布図であり、この図６ではオスの分布とメスの分布の分離がより顕著である。

以上は、Ｄａｙ７における５７８ｎｍの相対透過率のみを用いたが、５７８ｎｍの相対透過率は、血管以外にも種卵の大きさなどの影響を受けている。そこで、この影響を除く目的で、Ｄａｙ０における５７８ｎｍの相対透過率に対する比率を求めた。以下、Ｄａｙ０における５７８ｎｍの相対透過率に対するＤａｙ７における５７８ｎｍの相対透過率の比を、単にＤａｙ０に対するＤａｙ７の相対透過率の比という。

図７は、Ｄａｙ０に対するＤａｙ７の相対透過率の比（減少率）を雛の性別に分けて図示したものである。この図７から、メスの相対透過率の比がオスのものより有意に小さいことがわかる。
この図７において例えばＤａｙ０に対するＤａｙ７の相対透過率の比が、０．０２２より大きい種卵を抽出すれば、２９個が抽出される。これは全体７６個中の３８%である。２９個中２１個がオスであり、オスの比率は、７２％である。すなわち、特定時点（この場合はＤａｙ７）の５７８ｎｍの相対透過率とそれより前の時点（この場合はＤａｙ０）の５７８ｎｍの相対透過率の比を用いて血管や血液の生成の程度を推定し、その値が定められた閾値より大きいものを抽出すれば、オスの孵化する確率の高い卵が抽出される。先のＤａｙ７のみの場合と比較して、オス率が変化せず、抽出率が向上している。

図８は、横軸に初期卵重をとり、縦軸にＤａｙ０に対するＤａｙ７の相対透過率の比をとった散布図であり、この図８ではオスの分布とメスの分布の分離がより顕著である。

図９は、Ｄａｙ６における５７８ｎｍの相対透過率である。また、図１０は、Ｄａｙ０における５７８ｎｍの相対透過率に対するＤａｙ６における５７８ｎｍの相対透過率の比である。これらの図９及び図１０から生存卵と未受精卵の差異が顕著である。未受精卵であれば、例えば、図９、図１０で示したＤａｙ６やその前後（後述する第１時点付近）では、種卵内の血管及び／又は血液の形成が受精卵と異なり成されないからである。なお、生存卵と未受精卵の差異は、最も早い段階で、かつ、差異が明確となっている時点がＤａｙ６であるが、Ｄａｙ７でもＤａｙ６に準じた結果が得られる。

図８の散布図の横軸の変数Ｘ１（＝初期卵重）と、縦軸の変数Ｘ２（＝Ｄａｙ０に対するＤａｙ７の相対透過率の比）との２つの変数を説明変数とし、性別（オス＝１、メス＝０）を目的変数として、ロジスティック回帰分析で判別式を求めた。

ロジスティック回帰分析におけるオスメス判別値Ｙの分布を図１１に示す。
ロジスティック回帰分析におけるオスメス判別値Ｙは、Ｙの値が１に近いほどオスである確率が高く、逆に０に近いほどメスである確率が高い。
したがって、Ｙ≧０．５をオスと判定し、Ｙ＜０．５をメスとした場合、種卵は２つのグループ、Ｍ判定グループとＦ判定グループに２分される。すなわちＭ判定グループは、オスの雛が孵化する確率が高い種卵が含まれるグループであり、Ｆ判定グループは、メスの雛が孵化する確率が高い種卵が含まれるグループである。この分類を２分法と言う。

２分法による判別結果を以下の表１に示す。
表１の羽根鑑別オス、羽根鑑別メスは、雛が孵化したのちに羽根鑑別によりオスとメスが鑑別された雛の羽数である。Ｍ判定、Ｆ判定はそれぞれ、Ｍ判定グループとＦ判定グループに判定された種卵の数である。抽出数は、各グループに分類された種卵の個数であり、抽出率は全体に対する各グループが占める割合である。オス率は各グループ中でのオスが孵化した割合である。

Ｙ≧０．６をオスと判定し、Ｙ＜０．４をメスと判定し、それ以外を不明と判定した場合、種卵は、３つのグループ、Ｍ判定グループ、Ｇ判定グループ、Ｆ判定グループに３分される。ここでＭ判定グループは、オスの雛が孵化する確率が高い種卵が含まれるグループであり、Ｆ判定グループは、メスの雛が孵化する確率が高い種卵が含まれるグループであり、Ｇ判定グループは、オスとメスの確率が拮抗するグループである。この分類を３分法と言う。

３分法による判別結果を以下の表２に示す。
表２の羽根鑑別オス、羽根鑑別メスは、雛が孵化したのちに羽根鑑別によりオスとメスが鑑別された雛の羽数である。Ｍ判定、Ｇ判定、Ｆ判定はそれぞれ、Ｍ判定グループ、Ｇ判定グループ、Ｆ判定グループに判定された種卵の数である。抽出数は、各グループに分類された種卵の個数であり、抽出率は全体に対する各グループが占める割合である。オス率は各グループ中でのオスが孵化した割合である。

以上より、Ｄａｙ０に対するＤａｙ７の相対透過率の比の性差、及び、これらと初期卵重との組み合わせにより、種卵から孵化する雛の性別を予測し、３分法で４２％の卵をＭ判定グループとして抽出して、そこから約８４％の卵からオスが孵化することが分かった。上述したＤａｙ０に対するＤａｙ７の相対透過率の比のみの場合と比較して、抽出率もオス率も向上している。

以上、図８の散布図の横軸の変数Ｘ１と縦軸の変数Ｘ２を説明変数とし、性別を目的変数とするロジスティック回帰分析による雌雄の判別結果を示した。同様に図６の散布図の横軸の変数Ｘ１と縦軸の変数Ｘ２（Ｄａｙ７の５７８ｎｍの相対透過率）を説明変数として、性別を目的変数とするロジスティック回帰分析により雌雄判別を行っても、図８に基づく場合と同様の結果が得られる。
このようにＤａｙ０に対するＤａｙ７の５７８ｎｍの相対透過率の比、又は、Ｄａｙ７の５７８ｎｍの相対透過率に加えて、卵重を性別判定の変数に加えることにより、判別の性能の向上を図ることができる。
ここでは、卵重を判定に追加する変数としたが、卵重以外の卵の長径又は短径や、卵の体積など卵の大きさを示す変数であれば良い。

＜第１実施形態＞
以下に、本発明に係る種卵検査システムの第１実施形態について、図面を参照して説明する。

本実施形態の種卵検査システム１００は、種卵の孵卵段階において、種卵内の血管及び／又は血液の形成の程度を推定することにより、種卵から孵化する雛の性別を非破壊で選別するものである。

具体的に種卵検査システム１００は、図１２及び図１３に示すように、セッタートレイ２００に載置された複数の種卵を一挙に検査可能に構成されており、図示しない搬送機構により搬送されたセッタートレイ２００の下方から種卵に向かって光を照射する光照射部２と、セッタートレイ２００の上方に設けられて種卵を透過した光の強度を検出する光検出部３とを備えている。これら光照射部２及び光検出部３は、種卵の内部の胚、血管や血液の状態を示す状態情報をそれぞれ計測する計測部となる。つまり、種卵の内部の胚、血管や血液の状態を示す状態情報は、種卵を透過した光の強度又はそれを用いて求まる値（例えば透過率）である。なお、図１２は、セッタートレイ２００の搬送状態を示しており、図１３は、複数の種卵の計測状態を示している。

まず、セッタートレイ２００について説明すると、このセッタートレイ２００は、図１４に示すように、種卵が載置される例えば正６角形の卵座２０１を同一平面上に複数有するものである。なお、図１４には、計４２個の卵座が６行×７列で設けられており、４２個の種卵を載置可能に構成されたものを示している。

また、各卵座２０１は、その底面が下方に開口するとともに、種卵を保持する１又は複数の突起部２０２を有している。そして各卵座２０１は、図１５に示すように、上下方向において突起部２０２のほかに光を遮るものが無いように構成されている。なお、孵卵場で用いられるセッタートレイ２００には、ここに例示するもの以外に種々の形状のものが存在するが、各卵座２０１には卵を保持する突起部２０２のほかに光を遮るものは存在しないなど後述の光計測に必要な要件は、共通に満たされており、本実施形態は、このセッタートレイ２００の形状に限定されるものではない。

そして、このセッタートレイ２００は、図示しない搬送機構により所定のトレイ搬送方向に沿って搬送され（図１２参照）、所定の検出位置において一旦停止されて、光照射部２により光が照射されるとともに、種卵を透過した光が光検出部３により検出される（図１３参照）。

光照射部２は、血管や血液に吸収される波長を有する光を照射するものである。血管や血液に吸収される波長は、具体的にはヘモグロビンやミオグロビンに吸収される波長である。具体的には、検出位置にあるセッタートレイ２００に載置された複数の種卵に対応して設けられた複数の発光ダイオード（ＬＥＤ）２１である。複数のＬＥＤ２１は、５７８ｎｍ付近に発光中心波長を有するＬＥＤであり、本実施形態では、５７４ｎｍに発光中心波長を有するものである。なお、光照射部２としては、５７８ｎｍ付近に発光中心波長を有するレーザであっても良い。

光検出部３は、各ＬＥＤ２１に正対するように設けられた複数のフォトダイオード（ＰＤ）３１である。各ＰＤ３１は、個々に独立した黒色の遮光性と柔軟性を備えた素材でできた吸盤３２内に収められ、各吸盤３２とともにヘッド３３に固定されている。このヘッド３３は、図示しない昇降機構により吸盤３２が種卵に密着する計測位置Ｍ（図１３、図１５参照）と、当該計測位置Ｍから上方に離間してセッタートレイ２００が搬送される退避位置Ｎ（図１２参照）との間を移動する。

本実施形態の光照射部２は５７４ｎｍに発光中心波長をもつＬＥＤであるが、この５７４ｎｍであっても、上述した検証実験における５７８ｎｍの相対透過率とほぼ同様の結果が得られる。その結果を図１６及び図１７に示す。なおこの場合において、セッタートレイ２００上に合成樹脂製の同一形状の模擬卵を載置して、種卵の計測前にその透過光の強度を計測することにより、５７４ｎｍの相対透過率を求めることができる。

図１６は、Ｄａｙ０に対するＤａｙ７の相対透過率の比（減少率）を雛の性別に分けて図示したものである。この図１６から、メスの相対透過率の比がオスのものより有意に小さいことがわかる。

図１７は、横軸に初期卵重をとり、縦軸にＤａｙ０に対するＤａｙ７の相対透過率の比をとった散布図であり、この図１７ではオスの分布とメスの分布の分離がより顕著である。

そして、本実施形態の種卵検査システム１００は、図１２及び図１３に示すように、計測部である光照射部２及び光検出部３により得られた複数日の状態情報を個別の種卵毎に関連づけて記録する記録部４と、記録部４に記録された複数日の状態情報の個別の種卵毎に関連づけられた記録に基づいて種卵を選別する種卵選別部５とを備えている。

同一セッタートレイ２００上の個別の種卵は、卵座２０１の位置によって特定される。具体的には、セッタートレイ２００内の列の番号と行の番号を指定することで卵座２０１の位置が特定され、その上の種卵が特定される。
複数のセッタートレイ２００が使用される場合は、セッタートレイ２００を識別するトレイＩＤをバーコードなどの手段で付与することによりトレイを識別できる。
したがって、記録日などの記録の順序情報、トレイＩＤ、卵座の列番号、行番号に対応づけて状態情報を記録することで、複数日の状態情報を個別の種卵毎に関連づけて記録することができる。
また、記録の順序情報、トレイＩＤ、卵座の列番号、行番号をもとに記録を検索することで、同一の種卵に対する複数日にわたる状態情報を関連付けて参照できる。

なお、この記録部４及び種卵選別部５は、ＣＰＵ、内部メモリ、入出力インターフェース、ＡＤ変換部等の専用乃至汎用のコンピュータにより構成されている。そして、内部メモリに格納された種卵検査プログラムにしたがってＣＰＵやその他の周辺機器が協働することによって、記録部４及び種卵選別部５としての機能が発揮される。また、記録部４及び種卵選別部５は、物理的に一体のコンピュータにより構成されたものであっても良いし、それぞれ物理的に別体をなすコンピュータにより構成されたものであっても良い。

以下、各部４、５について説明する。
記録部４は、光検出部３のＰＤ３１により得られた電圧値、又は、ＰＤ３１により得られた電圧値と事前に求めておいた合成樹脂ブロックなどの模擬卵における電圧値との比である相対透過率の値を個別の種卵毎に関連づけて記録するものである。なお、模擬卵における電圧値は、種卵と同様に、セッタートレイ２００上に模擬卵を載置して取得する。

また、記録部４は、セッタートレイ２００の識別子（例えばセッタートレイ２００に設けられたバーコード情報など。図示していない。）とともに当該セッタートレイ２００に載置された種卵の位置情報（卵座位置）及びその種卵を透過した光の強度信号（電圧値）を記録する。

本実施形態の記録部４には、複数日における状態情報として、以下の（１）及び（２）が記録される。なお、記録部４には、入卵前に測られた初期卵重も種卵の状態情報として併せて記録される。
（１）孵卵開始から第１設定日（本実施形態では７日目、以下、同じ）におけるＰＤ３１で検出した電圧値、又はそれにより得られた相対透過率
（２）第１設定日以前の第２設定日（本実施形態では孵卵器に入卵する直前（孵卵開始前）である０日目）におけるＰＤ３１で検出した電圧値、又はそれにより得られた相対透過率

前記第１設定日及び第２設定日における電圧値の取得は、同じセッタートレイ２００を検出位置に搬送して、複数のＬＥＤ２１により当該セッタートレイ２００上の複数の種卵に光を照射して、それぞれの種卵を透過した光を複数のＰＤ３１で検出することにより行われる。

前記第１設定日は、孵卵７日目に限られず、その前後（例えば孵卵５日目〜８日目のいずれか）であってもよい。性別判別に最適な第１設定日は、鶏種や入卵前の予備加温の条件で変化し、別途、本発明の検証実験と同様の方法で決定される。
また、前記第２設定日は、孵卵０日目に限られないが好ましくは、孵卵０日目よりも前あるいは直後であればよい。性別判別に最適な第２設定日は、鶏種や入卵前の予備加温の条件で変化し、別途、本発明の検証実験と同様の方法で決定される。

種卵選別部５は、記録部４に記録された複数日の状態情報（ＰＤ３１の電圧値又は相対透過率）の個別の種卵毎に関連づけられた記録に基づいて、種卵内の血管及び／又は血液の形成の程度を推定する形成程度推定部５１と、形成程度の推定値から種卵から孵化する雛の性別を判別する性別判別部５２とを有する。

形成程度推定部５１は、第１設定日の状態情報である電圧値と第２設定日の状態情報である電圧値との比を形成程度の推定値とする。なお、両者の比を取ることにより、測定機器の機差及びトレイの種類などの影響に加えて、種卵のサイズや色などの種卵の属性の影響を低減することができ、また、模擬卵を用いた測定の手間を省くことができる。

性別判別部５２は、形成程度推定部５１により得られた推定値を取得して、当該推定値に基づいて、種卵から孵化する雛の性別を判別するものである。また、本実施形態の性別判別部５２は、推定値に基づいて、未受精卵も判別できるように構成されている。言い換えれば、この種卵検査システム１００は、第１設定日の種卵を透過した光の強度に基づいて、未受精卵か否かを判別する未受精判別部を備えている。

ここで、性別判別部５２は、閾値の設定により、Ｍ判定グループ及びＦ判定グループに選別する２分法、又は、Ｍ判定グループ、Ｇ判定グループ及びＦ判定グループに選別する３分法の何れかとして構成される。

Ｍ判定グループは、オスの雛が孵化する確率が高い種卵が含まれるグループである。
Ｆ判定グループは、メスの雛が孵化する確率が高い種卵が含まれるグループである。
３分法にした場合、Ｇ判定グループは、Ｍ判定グループ、Ｆ判定グループの何れにも属さない種卵が含まれるグループである。なお、Ｇ判定グループは、Ｍ判定グループとＦ判定グループとの中間のグループであり、雌雄いずれの雛が孵化する確率が高いとは言えないグループである。

つまり、性別判別部５２が、Ｍ判定グループ及びＦ判定グループに選別する２分法の場合には、Ｍ判定グループとＦ判定グループとに選別するためのＭ／Ｆ閾値が設定される。なお、このＭ／Ｆ閾値は、予め性別判別部に入力されている。

一方、性別判別部５２が、Ｍ判定グループ、Ｇ判定グループ及びＦ判定グループに選別する３分法の場合には、図１８に示すように、Ｍ判定グループとＧ判定グループとに選別するためのＭ／Ｇ閾値、及び、Ｇ判定グループとＦ判定グループとに選別するためのＧ／Ｆ閾値が設定される。なおこれらのＭ／Ｇ閾値及びＧ／Ｆ閾値は、予め性別判別部５２に入力されている。
また、性別判別部５２は、上記の２分法及び３分法のいずれにおいても、未受精卵を選別するための未受精卵閾値が設定されている。この未受精卵閾値は、Ｍ／Ｆ閾値及びＭ／Ｇ閾値よりも大きい値としてある。

以上により種卵選別部５において形成程度推定部５１により得られた推定値を用いて性別判別部５２により雌雄判別を行うことができる。

＜第１実施形態の効果＞
このように構成した本実施形態の種卵検査システム１００によれば、異なる複数日（Ｄａｙ０及びＤａｙ７。ただし、これに限られない。）における状態情報から、胚の成長速度などの時間的変化を把握することができる。この時間的変化は、種卵内の血管及び／又は血液の形成の程度を表しており、オスの胚とメスの胚とで異なることから、種卵から孵化する雛の性別を判別することができ、種卵を選別することができる。

すなわち、早い段階で、おおよその雛の歩留まりを予測できる。具体的には、そのロット中からオスの多く含まれるグループを選別して、ワクチン卵へと転用したり、販売可能な雛の予測数量が受注数より多い場合は転売先を探したり、少ない場合は同業他社から雛を調達して不足分を補うなどの行動をとることができる。

なお、本実施形態では、推定値に基づいて、未受精卵も判別できるように構成されているので、早い段階で、未受精卵や中止卵を除くことにより、移卵作業時点での腐敗卵の発生の可能性を減らせたり、早い段階で、雛の歩留まり予測や未受精卵率を把握して、雛の過不足が生じないような最適な種卵の調達数量の決定や、種鶏農場の親鶏群の管理に役立てたりすることができる。

また、複数日の状態情報を記録部４に記録させているので、異なる日に計測された状態情報の管理を容易にすることができる。さらに、複数日の状態情報を記録部４に記録させているので、孵化率に影響を与えない時期まで待った後に、種卵の選別及び移し変える作業を行うことができる。

＜第２実施形態＞
次に、本発明に係る種卵検査システムの第２実施形態について、図面を参照して説明する。

本実施形態の種卵検査システム１００は、記録部４の構成及び記録部４にデータを書き込む構成が前記実施形態とは異なる。
具体的に本実施形態の記録部４は、図１９に示すように、セッタートレイ２００に設けられたＲＦＩＤ等の非接触ＩＣタグである。

図１９に示す例では、種卵検査システム１００は、光検出部３により得られた光の強度（電圧値）を非接触ＩＣタグ４に書き込む光強度書込部６と、非接触ＩＣタグ４に記録された光強度を読み取る光強度読取部７と、性別判別部５２により得られた判別結果を非接触ＩＣタグ４に書き込む判別結果書込部８とを備えている。なお、光強度書込部６、光強度読取部７及び判別結果書込部８は、非接触ＩＣタグ４に対してデータの書き込み及び読み取りを行うリーダライタにより構成されている。この場合の種卵の非破壊検査の手順は以下である。

入卵前のＤａｙ０においてセッタートレイ２００を検出位置に搬送して、複数のＬＥＤ２１により当該セッタートレイ２００上の複数の種卵に光を照射して、それぞれの種卵を透過した光を複数のＰＤ３１で検出する。種卵検査システム１００は、そのＰＤ３１の電圧値を取得する。このＤａｙ０の電圧値は、種卵検査システム１００の光強度書込部６により、そのセッタートレイ２００に設けられた非接触ＩＣタグ４に、種卵の位置情報（卵座位置）とともに書き込まれる。なお、この非接触ＩＣタグ４には、初期卵重も併せて記録される。

孵卵開始から７日目（Ｄａｙ７）においてセッタートレイ２００を再び検出位置に搬送して、複数のＬＥＤ２１により当該セッタートレイ２００上の複数の種卵に光を照射して、それぞれの種卵を透過した光を複数のＰＤ３１で検出する。種卵検査システム１００は、そのＰＤ３１の電圧値を取得する。このＤａｙ７の電圧値は、種卵検査システム１００の光強度書込部６により、そのセッタートレイ２００に設けられた非接触ＩＣタグ４に、種卵の位置情報（卵座位置）とともに書き込まれる。

種卵検査システム１００の光強度読取部７は、非接触ＩＣタグ４からＤａｙ０の電圧値及びＤａｙ７の電圧値を取得する。そして、形成程度推定部５１は、Ｄａｙ０の電圧値及びＤａｙ７の電圧値を用いて推定値を算出し、性別判別部５２は、その推定値と所定の閾値とを比較して、各種卵毎に孵化する雛の性別を判別する。この判別結果は、種卵検査システム１００の判別結果書込部８により、そのセッタートレイ２００の非接触ＩＣタグ４に、種卵の位置情報（卵座位置）とともに書き込まれる。

その後、孵卵開始から例えば９日目（第１設定日以降の第３設定日の一例、９日目には限られない。）に、セッタートレイ２００は、仕分け装置４００に搬送される。この仕分け装置４００に設けられた判別結果読取部（具体的にはリーダライタ）４０１により、非接触ＩＣタグ４に記録された判別結果が読み取られて、仕分け装置４００により複数の種卵は、少なくともＭ判定グループ及びＦ判定グループに仕分けされる。なお、図１９では、Ｇ判定グループ及び未受精卵（廃棄）にも仕分ける場合を示している。

＜第２実施形態の効果＞
本実施形態の種卵検査システム１００は、セッタートレイ２００に設けられた非接触ＩＣタグ４に、光検出部３により得られた光の強度及び性別判別部５２により得られた判別結果を書き込むように構成されているので、セッタートレイ２００毎の種卵のデータ管理が容易となる。また、セッタートレイ２００に設けられた非接触ＩＣタグ４に性別判別部５２により得られた判別結果を書き込み、後日、その判別結果を用いて種卵の仕分けを行うことができるので、雌雄判別の時期が孵卵初期の場合であっても、種卵の移し変える作業は孵化率に影響を与えない時期まで待つことができる。

なお、本発明は前記第１、第２実施形態に限られるものではない。
＜種卵仕分けシステムの変形例（図２０）＞
図２０に種卵検査システム１００を組み込んだ種卵仕分けシステムを示す。このシステムでは、光照射部２及び光検出部３と、種卵選別部５とが物理的に分離されて構成されている。この例では、種卵仕分けシステムは、光検出部３により得られた光の強度を非接触ＩＣタグ４に書き込む光強度書込部６を備えている。また、非破壊検査装置１００とは別の装置（例えば仕分け装置４００）は、非接触ＩＣタグ４に記録された光強度を読み取る光強度読取部７と、形成程度推定部５１及び性別判別部５２とを備えている。なお、光強度書込部６及び光強度読取部７は、同一又は別々にリーダライタにより構成されている。この場合の種卵の非破壊検査の手順は以下である。

入卵前のＤａｙ０においてセッタートレイ２００を検出位置に搬送して、複数のＬＥＤ２１により当該セッタートレイ２００上の複数の種卵に光を照射して、それぞれの種卵を透過した光を複数のＰＤ３１で検出する。種卵検査システム１００は、そのＰＤ３１の電圧値を取得する。このＤａｙ０の電圧値は、種卵検査システム１００の光強度書込部６により、そのセッタートレイ２００の非接触ＩＣタグ４に、種卵の位置情報（卵座位置）とともに書き込まれる。なお、この非接触ＩＣタグ４には、初期卵重も併せて記録される。

孵卵開始から７日目（Ｄａｙ７）においてセッタートレイ２００を再び検出位置に搬送して、複数のＬＥＤ２１により当該セッタートレイ２００上の複数の種卵に光を照射して、それぞれの種卵を透過した光を複数のＰＤ３１で検出する。種卵検査システム１００は、そのＰＤ３１の電圧値を取得する。このＤａｙ７の電圧値は、種卵検査システム１００の光強度書込部６により、そのセッタートレイ２００の非接触ＩＣタグ４に、種卵の位置情報（卵座位置）とともに書き込まれる。

その後、孵卵開始から例えば９日目に、セッタートレイ２００は、仕分け装置４００に搬送される。この仕分け装置４００に設けられた光強度読取部７により、非接触ＩＣタグ４に記録されたＤａｙ０及びＤａｙ７の電圧値が読み取られる。そして、仕分け装置４００の形成程度推定部４は、Ｄａｙ０の電圧値及びＤａｙ７の電圧値を用いて推定値を算出し、性別判別部５は、その推定値と所定の閾値とを比較して、各種卵毎に孵化する雛の性別を判別する。この判定結果により、仕分け装置４００は、複数の種卵を、少なくともＭ判定グループ及びＦ判定グループに仕分ける。なお、図２０では、Ｇ判定グループ及び未受精卵（廃棄）にも仕分ける場合を示している。
この変形例や前述した第２実施形態のように、記録部としてＲＦＩＤ等の非接触ＩＣタグを用いる場合には、非接触ＩＣタグに全ての情報を記録させるものに限らず、例えば、非接触ＩＣタグにトレイＩＤや識別コードなど一部の情報のみを記憶させる一方、この非接触ＩＣタグと別の場所に第２の記録部を設けて前記一部の情報に関連した個別情報を記録するようにしてもよい。

＜サイズ情報をさらに利用して雛の性別を予測する変形例＞
種卵の透過率は種卵の卵重や長径、短径などのサイズにも依存する。このサイズの影響を考慮して、形成程度推定部５１は、変数Ｘ１（＝初期卵重）と、変数Ｘ２（＝Ｄａｙ０に対するＤａｙ７の相対透過率の比）との２つの変数を説明変数とし、性別（オス＝１、メス＝０）の予測値Ｙを求めるロジスティック回帰分析を行うものであり、性別判別部５２は、この予測値Ｙを用いて、種卵から孵化する雛の性別を判別するものであってもよい。この場合、性別判別部５２は、例えばＹ≧０．５をオス、Ｙ＜０．５をメスとした２分法による判別を行ってもよいし、例えばＹ≧０．６をオス、Ｙ＜０．４をメスとした３分法による判別を行ってもよい。

このように、形成程度推定部５１は、前記第１設定日以前の第２設定日又は孵卵開始前の前記種卵のサイズ情報をさらに用いるものであってもよく、上述した変数Ｘ１は、初期卵重以外に卵の長径や短径、又は卵の体積などが考えられ、そのサイズ情報の計測時期も種々変更可能である。

この変形例では、入卵前に測定するのが適切な種卵のサイズ情報と、Ｄａｙ０とＤａｙ７に測定するのが適切な他の情報とを組み合わせて用いている。これは、入卵前または入卵後の最適な日にちに測定された複数の証拠の利用により、性別やその他の判定精度の向上を図ることができ、複数日の状態情報の個別の種卵毎に関連づけられた記録に基づいて種卵を選別するという種卵検査システムを有効に活用した例と言える。同様に、図６で示したような初期卵重とＤａｙ７（単一の測定日）における５７８ｎｍの相対透過率との関係のようなものでも、この種卵検査システムの利点を活かすことができるのはもちろんである。

なお、変数Ｘ２についても、Ｄａｙ０に対するＤａｙ７の相対透過率の比に限られず、Ｄａｙ０とＤａｙ７の相対透過率の差であってもよい。Ｄａｙ０やＤａｙ７は、第１設定日や第２設定日の一例であることは言うまでもない。さらに、予測値Ｙを求める際には、変数Ｘ１と変数Ｘ２との関係性を求めることができるものであればロジスティック回帰分析に限られず種々変更可能である。

＜複数の波長の光を利用する変形例＞
前記実施形態の光照射部２は、血管や血液に吸収される波長を有する光（５７４ｎｍ）を射出するＬＥＤを用いたものであったが、このＬＥＤに加えて、孵卵後期の胚の心拍測定に適した波長を有する光（例えば８７０ｎｍ）を射出するＬＥＤを備えるものであってもよい。この場合、孵卵初期においては５７４ｎｍのＬＥＤを用いて種卵の雌雄判別及び未受精卵の判別を行うとともに、孵卵後期においては８７０ｎｍのＬＥＤを用いて非生存卵（発育中止卵および未受精卵）の判別を行うことができる。このような構成であれば、図２１に示す雛の生産工程において、入卵作業時に５７４ｎｍのＰＤの電圧値を測定し、中間検卵作業時に５７４ｎｍのＰＤの電圧値を測定して、未受精卵を選別して取り除くことができる。また、その後の性別仕分け作業によりオスの種卵を取り除いてワクチン卵に回すことができる。メスの種卵はそのまま孵卵が継続されて、移卵時検卵作業時に８７０ｎｍのＰＤの電圧値を測定して非生存卵を選別して取り除くことができ、移卵作業によりメスの生存卵をハッチャーに移卵して孵化させることができる。

前記実施形態の形成程度推定部５１は、前記第１設定日と第２設定日又は孵卵開始前の光の強度に基づいて、前記種卵内の血管及び／又は血液の形成の程度を推定する際に、各設定日の光の強度の比を用いていたが、各設定日間の光の強度の差を用いるようにするなど種々変更可能である。

＜胚の成長状態から雛の雌雄を予測する変形例＞
前記実施形態の種卵選別部は、血管や血液の状態情報から予想される種卵から孵化する雛の性別に基づいて種卵を選別するものであったが、複数日の状態情報により求まる種卵の内部の胚の成長状態（例えば成長速度のレベル）に基づいて、種卵を選別するものであっても良い。
このとき、記録部４には、複数日における状態情報として、以下の（１）及び（２）が記録される。
（１）孵卵開始から第１設定日（例えば１４日目）におけるＰＤ３１で検出した電圧値、又はそれにより得られた相対透過率
（２）第１設定日以後の第２設定日（例えば１８日目）におけるＰＤ３１で検出した電圧値、又はそれにより得られた相対透過率

前記実施形態では、孵卵開始から７日目（Ｄａｙ７）の相対透過率を用いて雌雄判別をするものであったが、図１２に示す検査装置において、セッタートレイ上の種卵を測定する際に個々の卵に対して得られるＬＥＤ電流（測定時にＬＥＤに流す電流値）と、ＰＤ受光電圧（ＰＤ受光電圧の時系列データの平均値）とを用いて、次式で計算される不透明度を用いて雌雄判別しても良い。なお、この場合のＬＥＤは、血液中のヘモグロビンや水分の吸収を受けにくい近赤外光を照射するものであり、具体的には、８７０ｎｍに発光中心波長を有するものである。

不透明度＝ＬＥＤ電流／ＰＤ受光電圧
＝ＬＥＤ電流／（卵の透過率×ＬＥＤ電流×ＬＥＤ発光効率×ＰＤ感度）
＝Ｋ×（１／卵の透過率）
ただし、Ｋ＝１／（ＬＥＤ発光効率×ＰＤ感度）

なお、同一卵座での異なる測定日の不透明度の比を求めることにより、係数Ｋを消去することができる。

ＤａｙＭに対するＤａｙＮの不透明度の比
＝ＤａｙＮの不透明度／ＤａｙＭの不透明度
＝ＤａｙＭの卵の透過率／ＤａｙＮの卵の透過率

この式を用いることにより、測定に用いた卵座に依存しないので、卵固有の属性である透過率の比として、卵間で大小を論じることができる。

以下の実施例においては、Ｄａｙ３から不透明度の測定を開始したので、Ｄａｙ３に対するＤａｙＮの不透明度の比を用いる。

この実施例では、上述した検証実験と同様に、採卵鶏（ジュリア・ライト）の種卵を９０個入手して、入手時に全ての種卵の卵重（以下、初期卵重という。）を測定した。

Ｄａｙ３からＤａｙ１８まで２４時間毎に、発光中心波長８７０ｎｍのＬＥＤの光を種卵に照射して、当該種卵からの透過光をＰＤにより受光した。なお、この検査装置は、卵毎にＰＤの受光電圧が所定の範囲になるように、ＬＥＤに流す電流値を変化させた。

図２２は、雛が孵化した７６個の種卵において、横軸に初期卵重をとり、縦軸にＤａｙ３に対するＤａｙ８の不透明度の比を取った散布図であり、図２３は、雛が孵化した７６個の種卵において、Ｄａｙ３に対するＤａｙ８の不透明度の比を雛の性別に分けて図示したものである。この図２２及び図２３から、メスのＤａｙ３に対するＤａｙ８の不透明度の比がオスのものより有意に大きいことがわかる（ｔ検定ｐ値＜０．０３７）。

この図２３において例えばＤａｙ３に対するＤａｙ８の不透明度の比が、２．２５より小さい種卵を抽出すれば、３７個が抽出される。これは全体７６個中の４９%である。３７個中２３個がオスであり、オスの比率は、６２％である。このように、Ｄａｙ３に対するＤａｙ８の不透明度の比によっても性別判別することができる。

図２４は、横軸にＤａｙ３に対するＤａｙ８の不透明度の比をとり、縦軸にＤａｙ０に対するＤａｙ７の５７４ｎｍの相対透過率の比をとった散布図である。
図２４の散布図の横軸の変数Ｘ１（＝Ｄａｙ３に対するＤａｙ８の不透明度の比）と、縦軸の変数Ｘ２（＝Ｄａｙ０に対するＤａｙ７の相対透過率の比）との２つの変数を説明変数とし、性別（オス＝１、メス＝０）を目的変数として、ロジスティック回帰分析で判別式を求めた。

ロジスティック回帰分析におけるオスメス判別値Ｙの分布を図２５に示す。なお、ロジスティック回帰分析におけるオスメス判別値Ｙは、Ｙの値が１に近いほどオスである確率が高く、逆に０に近いほどメスである確率が高い。

Ｙ≧０．５５をオスと判定し、Ｙ＜０．５５をメスとした場合、種卵は２つのグループ、Ｍ判定グループとＦ判定グループに２分される（２分法）。

この２分法による判別結果を以下の表３に示す。

Ｙ≧０．６２をオスと判定し、Ｙ＜０．４２をメスと判定し、それ以外を不明と判定した場合、種卵は、３つのグループ、Ｍ判定グループ、Ｇ判定グループ、Ｆ判定グループに３分される（３分法）。

この３分法による判別結果を以下の表４に示す。

以上より、Ｄａｙ０に対するＤａｙ７の相対透過率の比の性差、及び、これらとＤａｙ３に対するＤａｙ８の不透明度の比の性差との組み合わせにより、種卵から孵化する雛の性別を予測し、３分法で３３％の卵をＭ判定グループとして抽出して、そこから約８０％の卵からオスが孵化することが分かった。上述したＤａｙ３に対するＤａｙ８の不透明度の比の性差のみの場合と比較して、抽出率もオス率も向上している。

なお、ここでは、胚の成長度合いを推定するために不透明度を用いていたが、「胚の状態情報」は他に、代謝の強さ（一例として心拍の強さ）などであってもよい。

＜胚の成長状態から後期中止胚を予測し選別する変形例＞
前記実施形態の種卵選別部は、種卵から孵化する雛の性別に基づいて種卵を選別するものであったが、複数日の状態情報により求まる種卵の内部の胚の成長状態（例えば成長速度のレベル）に基づいて、雄雌以外の観点から種卵を選別するものであっても良い。
以下、異なる日の不透明度の比の雌雄判別以外の応用例について説明する。
孵卵開始１８日目（Ｄａｙ１８）の移卵時の生死判別の検査において、心臓の動きや胚の活動が観測されて生存卵と判断されても、その全ての生存卵が孵化するとは限らない。

通常、移卵後に呼吸が肺呼吸に変わると、卵内の雛は内側から卵殻を破り小さな穴を開ける。この穴を啄痕（たくこん）という。多くの場合、啄痕が出ると数時間以内に雛が孵化するが、中には孵化しない場合がある。
これらの移卵時に生存していたにも関わらず孵化しなかった卵を「死籠り卵」という。死籠りになるか否かは、Ｄａｙ１８のデータのみでは判別しづらい。

一方、Ｄａｙ１４に対するＤａｙ１８の不透明度の比を見ることにより、孵卵後期の胚の急成長期の成長速度の違いを見ることができる。そして、本願発明者は、Ｄａｙ１４に対するＤａｙ１８の不透明度の比が小さい場合に「死籠り卵」になる可能性が高いことを見出した。

図２６に、Ｄａｙ３に対するＤａｙＮ（Ｎ＝３〜１８）の不透明度の比の変化を示す成長曲線を示す。この図２６から分かるように、Ｄａｙ１４に対するＤａｙ１８の不透明度の比が小さい場合は、「死籠り卵」になる可能性が高い。また、Ｄａｙ１４に対するＤａｙ１８の不透明度の比が小さい場合は、「後期中止卵」の可能性が高い。なお、後期中止卵は、バイタルサインの存在でも識別可能であるが、不透明度の比を用いることで、判別精度の改善を図ることができる。

前記第１設定日及び第２設定日における電圧値の取得は、前記実施形態と同様に、同じセッタートレイ２００を検出位置に搬送して、複数のＬＥＤ２１により当該セッタートレイ２００上の複数の種卵に光を照射して、それぞれの種卵を透過した光を複数のＰＤ３１で検出することにより行われる。

前記第１設定日は、孵卵１４日目に限られず、その前後（例えば孵卵１２日目〜１５日目のいずれか）であってもよい。また、第２設定日は、孵卵１８日目に限られず、その前後（例えば孵卵１６日目〜１９日目のいずれか）であってもよい。胚の成長速度の判別に最適な第１設定日及び第２設定日は、鶏種や入卵前の予備加温の条件で変化し、別途、本発明の検証実験と同様の方法で決定される。

このように種卵選別部が成長速度のレベルに応じて種卵を選別することにより、種卵の孵化時間を予測して孵化時間ごとに種卵を仕分けたり、雛の増体性等の属性を予測して属性ごとに種卵を仕分けたりできる。また、種卵を属性ごとに仕分けることにより、種卵ごとに最適な孵卵条件を設定することもできる。

＜その他の変形例＞
例えば、光照射部２として、ハロゲンランプなどの広帯域の波長の光源を用いてもよい。この場合、光検出部３としては、多波長の分光スペクトルを求める分光器を用いた構成とする。この構成の場合、光の強度とは、分光強度スペクトルを意味する。

このとき記録部４には、複数日における状態情報として、以下の（１）及び（２）が記録される。
（１）孵卵開始から第１設定日（本実施形態では７日目、上述した実施形態に準ずる。）での分光強度スペクトルの５７８ｎｍの値、又は、その近傍の値の平均値、或いは、孵卵開始から第１設定日での種卵に対する分光強度スペクトルと事前に求めておいた合成樹脂ブロックなどの模擬卵の分光強度スペクトルとの各波長の比である透過率スペクトルにおける５７８ｎｍの値、又は、その近傍の値の平均値
（２）第１設定日以前の第２設定日（本実施形態では孵卵器に入卵する直前（孵卵開始前）である０日目、上述した実施形態に準ずる。）での分光強度スペクトルの５７８ｎｍの値、又は、その近傍の値の平均値、或いは、孵卵開始から第１設定日での種卵に対する分光強度スペクトルと事前に求めておいた合成樹脂ブロックなどの模擬卵の分光強度スペクトルとの各波長の比である透過率スペクトルにおける５７８ｎｍの値、又は、その近傍の値の平均値

前記第１設定日及び第２設定日における分光強度スペクトルの取得は、同じセッタートレイ２００を検出位置に搬送して、光照射部２により当該セッタートレイ２００上の複数の種卵に光を照射して、それぞれの種卵を透過した光を光検出部３で検出することにより行われる。

また、種卵検査システム１００は、形成程度推定部５１により得られた形成程度を非接触ＩＣタグ４に書き込む形成程度推定値書込部を有しており、当該形成程度推定値書込部により非接触ＩＣタグ４に推定値を書き込むように構成してもよい。この場合、非接触ＩＣタグ４に書き込まれた推定値は、形成程度推定値読取部が設けられた種卵検査システム１００又は仕分け装置などにより読み取られて、性別判別部５２による性別判別に利用される。さらに、種卵検査システム１００は、判別結果書込部及び判別結果読取部をさらに備えてもよい。

前記各実施形態では、２つの設定日又は３つの設定日それぞれで取得した状態情報から種卵を選別するものであったが、４つ以上の設定日それぞれで取得した状態情報から種卵を選別するものであっても良い。このように４つ以上の状態情報を比較することによって、雛の性別、孵化時間、雛の増体性などの属性、胚の発育の中止、中止時期などを精度よく予測することができる。

上述した実施形態では、セッタートレイ上に静止した状態の種卵に光を照射して検査を行うようにしていたが、これに限られず、搬送途中の種卵に光を照射して次々と検査するようにしてもよいのはもちろんである。

また、状態情報を得るための計測部は、透過光を用いて計測するものには限られない。さらに、複数日の状態情報を得るために、複数種類の計測部を用いてもよいのはもちろんである。
その他、本発明は前記実施形態に限られず、その趣旨を逸脱しない範囲で種々の変形が可能であるのは言うまでもない。

本発明によれば、様々な観点から高速かつ非侵襲的に種卵を選別する種卵検査システムを提供することができる。

Claims

種卵の状態を示す状態情報を計測する計測部と、
前記計測部により得られた複数日の状態情報を個別の種卵毎に関連づけて記録する記録部と、
前記記録部に記録された前記複数日の状態情報の個別の種卵毎に関連づけられた記録から、前記種卵の互いに異なる設定日の状態情報同士の比較に基づいて、前記種卵を選別する種卵選別部とを備える種卵検査システム。
前記種卵選別部は、前記種卵の互いに異なる設定日の状態情報同士の比較により求まる前記種卵の内部の胚の成長状態に基づいて前記種卵を選別するものである、請求項１記載の種卵検査システム。
前記種卵選別部は、前記種卵の互いに異なる設定日の状態情報同士の比較により予想される前記種卵から孵化する雛の性別に基づいて前記種卵を選別するものである、請求項１記載の種卵検査システム。
種卵の状態を示す状態情報を計測する計測部により得られた複数日の状態情報を個別の種卵毎に関連づけて記録する記録部と、
前記記録部に記録された前記複数日の状態情報の個別の種卵毎に関連づけられた記録から、前記種卵の互いに異なる設定日の状態情報同士の比較に基づいて、前記種卵を選別する種卵選別部と、としての機能をコンピュータに備えさせることを特徴とする種卵検査プログラム。