▶ 日本ハム株式会社の特許一覧 ▶ 株式会社エヌ・ティ・ティ・データの特許一覧
特許7162749雌豚用発情判定装置、雌豚の発情判定方法および雌豚の発情判定プログラム
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2022-10-20
(45)【発行日】2022-10-28
(54)【発明の名称】雌豚用発情判定装置、雌豚の発情判定方法および雌豚の発情判定プログラム
(51)【国際特許分類】
A01K 29/00 20060101AFI20221021BHJP
【ＦＩ】
A01K29/00 D
【請求項の数】
12
(21)【出願番号】P 2021543770
(86)(22)【出願日】2020-09-01
(86)【国際出願番号】 JP2020033004
(87)【国際公開番号】W WO2021045033
(87)【国際公開日】2021-03-11
【審査請求日】2021-12-15
(31)【優先権主張番号】P 2019162720
(32)【優先日】2019-09-06
(33)【優先権主張国・地域又は機関】JP
(73)【特許権者】
【識別番号】000229519
【氏名又は名称】日本ハム株式会社
(73)【特許権者】
【識別番号】000102728
【氏名又は名称】株式会社エヌ・ティ・ティ・データ
(74)【代理人】
【識別番号】100079108
【弁理士】
【氏名又は名称】稲葉 良幸
(74)【代理人】
【識別番号】100109346
【弁理士】
【氏名又は名称】大貫 敏史
(74)【代理人】
【識別番号】100117189
【弁理士】
【氏名又は名称】江口 昭彦
(74)【代理人】
【識別番号】100134120
【弁理士】
【氏名又は名称】内藤 和彦
(72)【発明者】
【氏名】助川 慎
(72)【発明者】
【氏名】内田 大介
(72)【発明者】
【氏名】奥田 雅貴
(72)【発明者】
【氏名】吉田 樹
(72)【発明者】
【氏名】森田 尚樹
(72)【発明者】
【氏名】大城 祐介
【審査官】竹中 靖典
(56)【参考文献】
【文献】特開２００７－３７４５４（ＪＰ，Ａ）
【文献】特表２００９－５０１０３５（ＪＰ，Ａ）
【文献】特開２００５－２１１０８（ＪＰ，Ａ）
【文献】国際公開第２０１２／０１４４１９（ＷＯ，Ａ１）
(58)【調査した分野】(Int.Cl.，ＤＢ名)
Ａ０１Ｋ ２９／００
(57)【特許請求の範囲】
【請求項1】
ストールで飼育されている雌豚の単位時間あたりの起臥の頻度を計測する計測部と、設定された一定期間において前記計測部で繰り返し計測された複数の前記頻度に基づいて前記雌豚の発情を判定する判定部と
を備える雌豚用発情判定装置。
【請求項2】
前記計測部は、前記雌豚の起臥の状態として臥位、着座および起立のいずれであるかを認定し、認定したそれぞれの状態に対して予め設定されたスコアを累積加算した累積値を前記頻度とする請求項１に記載の雌豚用発情判定装置。【請求項3】
前記計測部は、臥位、着座および起立のそれぞれの状態が継続された時間に基づいて前記スコアを補正する請求項２に記載の雌豚用発情判定装置。【請求項4】
前記一定期間は、３日以上７日以下の期間であり、前記単位時間は、前記一定期間の間毎日設定される２時間以上６時間以下の時間である請求項１から３のいずれか１項に記載の雌豚用発情判定装置。【請求項5】
前記単位時間は、前記雌豚が餌を与えられて１２時間経過後から次に餌を与えられる前までの期間内に設定される請求項１から４のいずれか１項に記載の雌豚用発情判定装置。【請求項6】
前記単位時間は、前記ストールが設置された環境が暗い状態から明るい状態へ変化する時間帯を含むように設定される請求項１から５のいずれか１項に記載の雌豚用発情判定装置。【請求項7】
前記単位時間は、前記ストールが設置された周辺環境に人が居ない期間内に設定される請求項１から６のいずれか１項に記載の雌豚用発情判定装置。【請求項8】
前記判定部は、前記一定期間において繰り返し計測された複数の前記頻度を用いて前記雌豚の発情を識別する、予め機械学習された識別器を用いて前記雌豚の発情を判定する請求項１から７のいずれか１項に記載の雌豚用発情判定装置。【請求項9】
前記識別器は、サポートベクターマシンである請求項８に記載の雌豚用発情判定装置。【請求項10】
前記判定部は、前記一定期間における前記頻度が一旦低下し再度上昇する変化を示した場合に発情と判定する請求項１から９のいずれか１項に記載の雌豚用発情判定装置。【請求項11】
設定された一定期間において、ストールで飼育されている雌豚の単位時間あたりの起臥の頻度を繰り返して計測する計測ステップと、前記計測ステップで繰り返し計測された複数の前記頻度に基づいて前記雌豚の発情を判定する判定ステップと
を有する雌豚の発情判定方法。
【請求項12】
設定された一定期間において、ストールで飼育されている雌豚の単位時間あたりの起臥の頻度を繰り返して計測する計測ステップと、前記計測ステップで繰り返し計測された複数の前記頻度に基づいて前記雌豚の発情を判定する判定ステップと
をコンピュータに実行させる雌豚の発情判定プログラム。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、雌豚用発情判定装置、雌豚の発情判定方法および雌豚の発情判定プログラムに関する。
【背景技術】
【0002】
家畜を飼育する施設にセンサを設置して家畜の異常行動を検出するシステムが知られている（例えば、特許文献１参照）。
【先行技術文献】
【特許文献】
【0003】
【文献】特開２０１９－２４４８２号公報
【発明の概要】
【発明が解決しようとする課題】
【0004】
養豚において雌豚の発情の見極めは、年間の産子数を増やす観点や繁殖サイクルを適正な間隔に維持する観点等から重要である。しかし、従来技術では、対象雌豚の異常行動を検出することができても、発情といった特定の状態を検出することはできない。経験や勘に頼らず、雌豚の発情を精度良く見極めたいという要望が、飼育者から多く挙げられている。
【0005】
本発明は、このような問題を解決するためになされたものであり、雌豚の発情を精度良く見極める技術を提供するものである。
【課題を解決するための手段】
【0006】
本発明の第１の態様における雌豚用発情判定装置は、ストールで飼育されている雌豚の単位時間あたりの起臥の頻度を計測する計測部と、設定された一定期間において計測部で繰り返し計測された複数の頻度に基づいて雌豚の発情を判定する判定部とを備える。
【0007】
本発明の第２の態様における雌豚の発情判定方法は、設定された一定期間において、ストールで飼育されている雌豚の単位時間あたりの起臥の頻度を繰り返して計測する計測ステップと、計測ステップで繰り返し計測された複数の頻度に基づいて雌豚の発情を判定する判定ステップとを有する。
【0008】
本発明の第３の態様における雌豚の発情判定プログラムは、設定された一定期間において、ストールで飼育されている雌豚の単位時間あたりの起臥の頻度を繰り返して計測する計測ステップと、計測ステップで繰り返し計測された複数の頻度に基づいて雌豚の発情を判定する判定ステップとをコンピュータに実行させる。
【発明の効果】
【0009】
本発明により、観察者の経験や勘に頼らずとも、雌豚の発情を精度良く見極めることができる。
【図面の簡単な説明】
【0010】
【図1】本実施形態に係る判定装置を採用した養豚環境の全体像を示す図である。
【図2】判定装置のハードウェア構成を示す図である。
【図3】ストール内の雌豚の状態を説明する図である。
【図4】雌豚の起臥の頻度を累積スコアに換算する手法を説明する図である。
【図5】発情していない雌豚の累積スコアの推移を示す図である。
【図6】発情している雌豚の累積スコアの推移を示す図である。
【図7】起臥頻度の計測時間帯を説明する図である。
【図8】計測時間帯の違いによる累積スコアの推移の違いを説明する図である。
【図9】発情識別器を用いた発情判別の手法を説明する図である。
【図10】累積スコアを計測する計測ステップの処理について説明する図である。
【図11】対象となる雌豚の発情を判定する判定ステップの処理について説明する図である。
【図12】他の発情判定手法を説明する図である。
【発明を実施するための形態】
【0011】
以下、発明の実施の形態を通じて本発明を説明するが、特許請求の範囲に係る発明を以下の実施形態に限定するものではない。また、実施形態で説明する構成の全てが課題を解決するための手段として必須であるとは限らない。
【0012】
図１は、本実施形態に係る判定装置２００を採用した養豚環境の全体像を示す図である。養豚場では、観察対象となる雌豚１０１がストール１０２に収容されている。一つのストール１０２には一頭の雌豚１０１が収容されている。一つのストール１０２の大きさは、収容されている性成熟期の雌豚１０１が自力では旋回できない大きさであり、例えば幅８０ｃｍ、奥行き２００ｃｍ程度である。したがって、雌豚１０１の動作は、起き上がったり臥したりする動作に制限される。観察対象となる雌豚１０１が複数存在する場合には、それぞれが収容された複数のストール１０２が並設される。なお、ストール１０２の大きさは、観察対象となる雌豚の品種や個体差、飼育環境に応じて設定され得る。
【0013】
カメラユニット１１０は、観察対象となる雌豚１０１の全身を俯瞰して撮像できる撮像センサを備えており、撮像センサで撮像した画像を画像データに変換し、インターネット９００を介してサーバ２１０へ送信する。観察対象となる雌豚１０１が複数存在する場合には、それぞれのストール１０２に対してカメラユニット１１０を一つずつ設置しても良いし、複数のストール１０２の群のそれぞれに対してカメラユニットを設置しても良い。複数のストール１０２の群のそれぞれに対してカメラユニットを設置する場合には、群を構成するストール１０２のそれぞれに収容された雌豚１０１を同時に観察できるように、カメラユニット１１０の画角を調整する。
【0014】
管理施設には、観察対象となる雌豚１０１の発情を判定する判定装置２００が設置されている。判定装置２００は、サーバ２１０と、サーバ２１０に接続された表示モニタ２２０等によって構成され、サーバ２１０は、インターネット９００と接続されている。サーバ２１０は、インターネット９００を介してカメラユニット１１０から送られてくる画像データを受け取り、当該画像データから雌豚１０１の起臥の頻度を計測し、計測された頻度に基づいて雌豚１０１が発情しているか否かを判定する。サーバ２１０は、判定結果を表示モニタ２２０へ表示する。養豚場で作業する作業員から作業員端末１２０を介して判定結果を求められた場合には、インターネット９００を介して、判定結果を作業員端末１２０の表示部へ表示する。作業員端末１２０は、例えば、タブレット端末やスマートフォンである。
【0015】
なお、カメラユニット１１０と判定装置２００を接続する回線は、インターネット９００に限らず、イントラネット等であっても良い。養豚場内に管理施設が設けられるような場合には、近距離無線通信が採用されても良い。
【0016】
図２は、判定装置２００のハードウェア構成を示す図である。判定装置２００は、上述のように主にサーバ２１０と表示モニタ２２０によって構成される。表示モニタ２２０は、例えば液晶パネルを備え、演算部２３０が生成する映像信号を視認可能な映像に変換して表示する。サーバ２１０は、主に、演算部２３０、画像処理部２４０、データ蓄積部２５０、メモリ２６０および通信ユニット２７０を備える。
【0017】
演算部２３０は、例えばＣＰＵであり、メモリ２６０から読み込んだ各種プログラムを実行することにより、判定装置２００の全体を制御したり諸々の演算処理を実行したりする。例えば、計測部２３１としての処理を実行する場合には、画像処理部２４０と協働して、対象となる雌豚１０１の単位時間あたりの起臥の頻度を計測する。判定部２３２としての処理を実行する場合には、計測部２３１の計測結果を用いて当該雌豚１０１の発情を判定し、表示モニタ２２０や作業員端末１２０へ出力する。具体的な処理については、後に詳述する。
【0018】
画像処理部２４０は、例えば画像処理用のＡＳＩＣであり、カメラユニット１１０から受け取った画像データから対象雌豚の画像領域を切り出した姿勢判定画像を生成するなどの画像処理を実行する。データ蓄積部２５０は、例えばＨＤＤ（Hard Disc Drive）であり、観察対処となる雌豚１０１の識別情報と、識別情報に紐付けられた累積スコアとを蓄積する。累積スコアは、予め定められた単位時間の間に観察された対象雌豚１０１の起臥状態を数値化して累積加算したものである。累積スコアは、対象雌豚１０１がストール１０２に収容されている間は日ごとにカウントされ、日付情報と共にその累積スコアがデータ蓄積部２５０に記録される。
【0019】
メモリ２６０は、例えばＳＳＤ（Solid State Drive）であり、判定装置２００を制御するための制御プログラムや対象雌豚１０１の発情判定を行う発情判定プログラムの他にも、様々なパラメータ値、関数、ルックアップテーブル等を記憶している。特に、学習済みモデルである起臥識別器２６１と発情識別器２６２を格納している。起臥識別器２６１は、入力された姿勢判定画像に写る雌豚の起臥状態を識別する。発情識別器２６２は、入力された連日の累積スコアから対象となる雌豚が発情しているか否かを識別する。詳細については後述する。
【0020】
通信ユニット２７０は、例えば有線ＬＡＮユニットである。演算部２３０は、通信ユニット２７０を介して、インターネット９００に接続されたカメラユニット１１０へ画像データを要求し、これに呼応してカメラユニット１１０から送られてくる画像データを受信する。また、判定部２３２は、通信ユニット２７０を介して受信した作業員端末１２０からの要求に応じて、発情判定結果を当該通信ユニット２７０へ送信する。
【0021】
図３は、ストール１０２内の雌豚１０１の状態を説明する図である。本実施形態においては、観察対象となる雌豚１０１の起臥の状態を、臥位、着座および起立のいずれかの状態として認定する。
【0022】
図３（Ａ）は、雌豚１０１が起立状態である様子を示す。起立状態は、前肢および後肢ともに起き上がった姿勢の状態である。本実施形態においては、起立状態が観察された場合に、スコアα＝１．０を与える。図３（Ｂ）は、雌豚１０１が着座状態である様子を示す。着座状態は、前肢および後肢の一方が起き上がり他方が曲げられた姿勢の状態である。図は、雌豚１０１が後肢を曲げて臀部を接地させている様子を示している。本実施形態においては、着座状態が観察された場合に、スコアα＝０．５を与える。図３（Ｃ）は、雌豚１０１が臥位状態である様子を示す。臥位状態は、前肢および後肢を共に曲げたり側方へ投げ出したりして、脚に荷重が掛からない姿勢の状態である。図は、雌豚１０１が全脚を側方へ投げ出して側腹を接地させている様子を示している。本実施形態においては、臥位状態が観察された場合に、スコアα＝０とする。
【0023】
カメラユニット１１０は、このように変化する雌豚１０１の様子を、計測部２３１からの制御に従って周期的に撮像する。計測部２３１は、カメラユニット１１０から送られてくる画像データを取得し、画像処理部２４０へ引き渡す。画像処理部２４０は、画像データから観察対象である雌豚１０１の領域を切り出し、予め設定された画像処理を施して姿勢判定画像を生成する。予め設定された画像処理は、画像サイズの調整や特定色（例えば対象雌豚の肌の色）周りの輪郭強調等である。計測部２３１は、メモリ２６０から起臥識別器２６１を読み出して、画像処理部２４０が生成した姿勢判定画像を入力する。起臥識別器２６１は、姿勢判定画像に写る雌豚の起臥状態が、起立状態、着座状態、臥位状態のいずれかを認定結果として出力する。計測部２３１は、認定結果に応じてスコアαを確定させる。
【0024】
起臥識別器２６１は、機械学習によって学習された学習済みモデルである。起臥識別器２６１は、事前に学習装置によって作成される。具体的には、学習装置は、姿勢判定画像とその正解（起立状態、着座状態、臥位状態のいずれか）のセットである教師データを多数与えられ、機械学習の一つである教師あり学習を実行する。ここでは、教師あり学習として、画像認識に好適なＣＮＮ（畳み込みニューラルネットワーク）を用いる。教師あり学習によって学習を終えた起臥識別器２６１は、学習装置からメモリ２６０へ移され、上述の利用に供される。
【0025】
単位時間Ｔ0として例えば２時間を設定し、この２時間の間に対象雌豚１０１が起臥状態を何回変化させたかを評価する。すなわち、起臥の頻度を評価する。この評価の評価値として、上述のように、対象雌豚１０１が１回起立すればスコアα＝１．０を与え、１回着座すればスコアα＝０．５を与える。臥位している回数はスコアα＝０なので、評価値としてはカウントしない。このようにして単位時間Ｔ0の間に累積した累積スコアは、計測日における対象雌豚１０１の起臥の頻度を表す評価値になり得る。
【0026】
本実施形態においては、評価値の精度をさらに高めるために、起立、着座および臥位のそれぞれの状態が継続された時間も考慮して評価値を計算する。つまり、臥位、着座および起立のそれぞれの状態が継続された時間に基づいて累積されるスコアを補正する手法を採用する。図４は、本実施形態における、対象雌豚１０１の起臥の頻度を累積スコアに換算する手法を説明する図である。横軸は経過時間を表し、縦軸は起立（α＝１．０）、着座（α＝０．５）および臥位（α＝０）に対応するスコアαを表す。
【0027】
図は、時刻ｔｓに対象となる雌豚１０１の観察を開始し、単位時間Ｔ0が経過する時刻ｔｅまで観察を続けた観察結果の一例である。カメラユニット１１０は、例えば１フレーム／秒で雌豚１０１を撮像し、その画像データを判定装置２００へ送信する。計測部２３１は、受信するそれぞれの画像データを用いて雌豚１０１が起立、着座、臥位のいずれの状態であるかを認定し、スコアαを確定させる。すなわち、１秒ごとに雌豚１０１の状態に応じたスコアαを確定させる。そして、時刻ｔｓからの累積スコアに今回のスコアαを加算して累積スコアを更新する。この処理を単位時間Ｔ0が経過する時刻ｔｅまで続ける。
【0028】
計測部２３１は、このように累積スコアをカウントすることにより、起立、着座および臥位のそれぞれの状態が継続された時間を考慮した起臥の頻度を表す評価値を計算することができる。例えば、図において時刻ｔｓから時刻ｔ1までは起立状態（α＝１．０）が継続しているが、時刻ｔｓから時刻ｔ1までの時間が１０００秒だとすると、１秒ごとに１．０が加算されるので、１．０×１０００＝１０００が累積スコアに加算される。
【0029】
このようにして計算された累積スコアは、雌豚１０１の活動強度を表していると言える。すなわち、累積スコアが大きければ雌豚１０１は活発に活動していたと言え、累積スコアが小さければおとなしくしていたと言える。計測部２３１が一定期間に亘って毎日決められた時間帯に対象である雌豚１０１の累積スコアを計算し、その累積スコアを日ごとにデータ蓄積部２５０へ蓄積すれば、その間における雌豚１０１の活動強度の変化を知ることができる。
【0030】
雌豚の年間の産子数を増やす観点や繁殖サイクルを適正な間隔に維持する観点から、ストールに収容された雌豚の発情を精度良く見極めることは重要である。これまでは、雌豚が発情したか否かの判断は、熟練作業員の経験や勘に大きく依存していた。したがって、熟練作業員の経験や勘に頼らず、雌豚の発情を精度良く見極めたいという要望が、多くの飼育者から挙げられている。このような背景のもと、本願発明者は、研究を重ねた結果、雌豚の活動強度の変化と発情に相関があることを見出した。同時に、上述のような活動強度の評価手法を開発した。活動強度の変化と発情の関係を、上述の累積スコアを用いて説明する。
【0031】
図５は、発情していない雌豚の累積スコアの推移を示す図である。横軸は観察経過日数を表し、縦軸はそれぞれの観察日における累積スコアを表している。ここでは、４頭の雌豚（豚Ａ、豚Ｂ、豚Ｃ、豚Ｄ）の観察結果をプロットしている。
【0032】
それぞれの雌豚には個性があり、もともと活発に活動する雌豚もいれば、おとなしい雌豚もいる。そこで、それぞれの雌豚にとって、他の観察日と比較して相対的に活発に活動していると考えられる累積スコアを示した観察日をそれぞれ１日目としてプロットした。そして、それに続く３日分の累積スコアをプロットして、合計４日分の累積スコアの変化を表した。
【0033】
図示するように、これら４頭の雌豚の累積スコアの推移は、右肩下がりか横ばいである。これら４頭の雌豚のいずれにも発情は認められなかった。ここでは４頭の例を示したが、発情が認められなかった他の雌豚の累積スコアの推移も、およそ同様の結果であった。
【0034】
図６は、発情している雌豚の累積スコアの推移を示す図である。図５と同様に、横軸は観察経過日数を表し、縦軸はそれぞれの観察日における累積スコアを表している。ここでも、４頭の雌豚（豚Ｅ、豚Ｆ、豚Ｇ、豚Ｈ）の観察結果をプロットしている。
【0035】
図６の場合も図５の場合と同様に、それぞれの雌豚にとって、他の観察日と比較して相対的に活発に活動していると考えられる累積スコアを示した観察日をそれぞれ１日目としてプロットした。そして、それに続く３日分の累積スコアをプロットして、合計４日分の累積スコアの変化を表した。
【0036】
図示するように、これら４頭の雌豚の累積スコアの推移は、いずれも、一旦下がってから回復するＶ字状を示している。これら４頭の雌豚のいずれにも発情が認められた。ここでは４頭の例を示したが、発情が認められた他の雌豚の累積スコアの推移も、およそ同様の結果であった。なお、本願発明者は、１０５頭の雌豚から得られた１８１個の累積スコアの推移データを用いて図５および図６の検証を行った。
【0037】
以上の検証結果から、本願発明者は、一定期間に亘る活動強度の変化がＶ字状となった場合にその雌豚は発情している可能性が高い、という知見を得た。すなわち、上述のように観察対象となる雌豚の累積スコアを毎日計測して蓄積しておき、発情の有無を知りたくなった時点で一定期間分の累積スコアを抽出してその変化を確認すれば、その雌豚が発情しているか否かの判定を行うことができる。
【0038】
計測部２３１が累積スコアを繰り返し計測する一定期間は、Ｖ字状を初めて認識し得る３日以上が必要であり、また、２日目、３日目、４日目が累積スコアの底値となる場合を考慮して、最大で７日以下程度とすることが好ましい。また、雌豚の起臥頻度を計測する上記の単位時間Ｔ0を、どれくらいの長さにして、一日のうちのいつ設定するかについても、判定の精度を高める上で重要である。まず、単位時間Ｔ0を一日のうちのいつ設定するかについて説明する。なお、以下の説明においては、単位時間Ｔ0を、起臥頻度を計測する時間として計測時間Ｔ0と称する場合がある。
【0039】
図７は、起臥頻度の計測時間帯を説明する図である。原則的には、観察対象となる雌豚が、外部からの影響ではなく自発的に起きたり臥したりする活動を計測すべきであるので、自然環境のリズムに即した静かな時間帯に計測時間Ｔ0を設定することが好ましいと言える。図７（Ａ）は、雌豚に餌を与える給餌時間と計測時間Ｔ0の好ましい関係を示す図である。給餌時間は、発情の有無に関わらず雌豚が活発に活動する時間帯である。また、採餌後もしばらくの間は活発に活動する。したがって、計測時間Ｔ0は、雌豚が餌を与えられて１２時間経過後から次に餌を与えられる前までの期間内に設定されることが好ましい。図の例では、８時から９時の間が給餌時間に設定されており、９時から１２時間後の２１時から２４時の３時間が計測時間Ｔ0として設定されている。
【0040】
図７（Ｂ）は、養豚場付近の夜明け時間帯と計測時間Ｔ0の好ましい関係を示す図である。ここでは養豚場には窓が設けられており、夜明けと共に場内が暗い状態から明るい状態へ変化することを想定している。雌豚は、一般的に、暗い状態では活動が停滞し、明るくなってくると徐々に活動を開始する。したがって、計測時間Ｔ0は、ストールが設置された環境が暗い状態から明るい状態へ変化する夜明け時間帯を含むように設定されることが好ましい。図の例では、５時から５時半が夜明け時間帯である場合に、この時間帯を含む５時から７時までの２時間が計測時間Ｔ0として設定されている。なお、養豚場に窓が存在せず、屋内の明るさが照明によって制御されているような場合には、計測時間Ｔ0は、ストールが設置された環境が照明によって暗い状態から明るい状態へ変化する時間帯を含むように設定されると良い。
【0041】
図７（Ｃ）は、養豚場内で作業員が作業をする作業時間と計測時間Ｔ0の好ましい関係を示す図である。作業時間は、作業員が養豚場内の掃除をしたり飼育されている雌豚の健康状態をチェックしながら周回したりする時間であり、発情の有無に関わらず雌豚が活発に活動する時間帯である。したがって、計測時間Ｔ0は、ストールが設置された周辺環境に人が居ない期間内に設定されることが好ましい。図の例では、８時から１０時と１６時から１８時が作業時間に設定されており、４時から８時の４時間が計測時間Ｔ0として設定されている。
【0042】
計測時間Ｔ0として、図７（Ａ）の例では３時間、図７（Ｂ）の例では２時間、図７（Ｃ）の例では４時間を設定したが、具体的にどれくらいの時間に設定するかについては、給餌時間など考慮すべき要因との関係で決定すると良い。いずれの場合にしても、計測時間Ｔ0は、２時間以上６時間以下の時間であることが好ましいことが実験を通じて見出された。
【0043】
次に、計測時間帯の違いが累積スコアの推移にどのように現れるかについて説明する。図８は、計測時間帯の違いによる累積スコアの推移の違いを説明する図である。横軸は観察経過日数を表し、縦軸はそれぞれの観察日における累積スコアを表している。ここでは、発情を示す特定の雌豚に対して、毎日５時から７時までの間に計測した累積スコアをプロットして結んだ推移と、毎日１６時から１８時までの間に計測した累積スコアをプロットして結んだ推移とを表す。
【0044】
５時から７時までは、図７（Ａ）から（Ｃ）のいずれにおいても好ましい時間帯であり、１６時から１８時までは、給餌の影響が残り、作業員の作業時間と被る時間帯である。１６時から１８時までの間に計測した累積スコアの推移も確かにＶ字状を成すが、５時から７時までの間に計測した累積スコアの推移に比べて変化が緩やかである。後述する判定手法においては、Ｖ字状の変化が顕著に現れるほど精度良く発情を判定できるので、計測時間Ｔ0を、１６時から１８時までの時間帯に設定するよりは５時から７時までの時間帯に設定する方が好ましいと言える。
【0045】
計測部２３１が日々決められた時間帯に対象雌豚１０１の起臥状態を観察し、累積スコアを計算してデータ蓄積部２５０へ蓄積すれば、判定部２３２は、利用者の要求に従って、対象雌豚１０１の発情の有無を判定することができる。本実施形態において判定部２３２は、発情識別器２６２を利用して発情の有無を判定する。図９は、発情識別器２６２を用いた発情判別の手法を説明する図である。
【0046】
判定部２３２は、メモリ２６０から発情識別器２６２を読み出して、データ蓄積部２５０に蓄積された一定期間分の累積スコアを入力する。ここでは、一定期間を４日間として、１日目の累積スコアｘ1、２日目の累積スコアｘ2、３日目の累積スコアｘ3、４日目の累積スコアｘ4の４次元ベクトル（ｘ1，ｘ2，ｘ3，ｘ4）を発情識別器２６２へ入力する。
【0047】
本実施形態の発情識別器２６２はサポートベクターマシン（ＳＶＭ）であり、事前に学習装置で学習された学習済みモデルである。具体的には、一定期間分の累積スコアが「発情あり」または「発情なし」の正解ラベルと共に入力ベクトルとして与えられ、「発情あり」の空間と「発情なし」の空間を区分する識別関数が学習によって確定されている。このような学習を経て作成された発情識別器２６２は、学習装置からメモリ２６０へ移され、判定の利用に供される。
【0048】
判定部２３２が４次元ベクトル（ｘ1，ｘ2，ｘ3，ｘ4）を発情識別器２６２へ入力すると、発情識別器２６２は、「発情あり」を表す「１」または「発情なし」を表す「－１」を出力する。判定部２３２は、発情識別器２６２の出力である判定結果を表示モニタ２２０や作業員端末１２０へ出力する。例えば、表示モニタ２２０には、「管理番号○○／発情あり」のように、対象雌豚１０１の管理番号と共に発情の有無が表示される。なお、ここでは、４日分の累積スコアを４次元ベクトルとして入力したが、累積スコアに加え他の特徴量も入力ベクトルに加えたＳＶＭを作成することもできる。例えば、離乳日や気温を特徴量として入力するＳＶＭを作成することができる。また、学習済みモデルは、「発情あり」と「発情なし」の二分判定の出力を行うＳＶＭに限らず、例えば「発情確率８０％」のように多段階判定の出力を行うものであっても良い。
【0049】
以上説明したように、本実施形態における判定装置２００は、大きく分けて２つの処理ステップ、すなわち、累積スコアを計測する計測ステップと、対象となる雌豚の発情を判定する判定ステップとを実行する。そこで、それぞれの処理の流れを整理する。
【0050】
図１０は、累積スコアを計測する計測ステップの処理について説明するフロー図である。フローは、観察対象となる雌豚１０１がストール１０２へ収容され、システムオペレータにより継続的な観察の開始指示がなされた時点から開始する。
【0051】
計測部２３１は、ステップＳ１０１で、現在時刻を確認し、予め設定されている計測開始時刻になったか否かを判断する。計測開始時刻になったと判断したらステップＳ１０２へ進み、そうでないと判断したらステップＳ１０７へ進む。ステップＳ１０２へ進んだ場合、計測部２３１は、撮像を指示する指示信号をカメラユニット１１０へ送信する。カメラユニット１１０は、当該指示信号を受信すると雌豚１０１の撮像を実行し、生成した画像データを計測部２３１へ送信する。
【0052】
計測部２３１は、カメラユニット１１０から画像データを受信すると、ステップＳ１０３へ進み、当該画像データを画像処理部２４０へ引き渡して姿勢判定画像を生成させ、生成された姿勢判定画像を起臥識別器２６１へ入力する。そして、起臥識別器２６１が起立状態、着座状態、臥位状態のいずれかを認定結果として出力したら、当該認定結果に応じてスコアαを確定させる。
【0053】
計測部２３１は、ステップＳ１０４で、それまでの累積スコアに今回確定したスコアαを加算して累積スコアを更新する。ステップＳ１０５へ進み、計測部２３１は、現在時刻を確認し、計測開始時刻から予め設定された単位時間Ｔ0が経過したか否かを判断する。位時間Ｔ0が経過したと判断したらステップＳ１０６へ進み、そうでないと判断したらステップＳ１０２へ戻って雌豚１０１の観察を継続する。なお、ステップＳ１０２へ戻って再びカメラユニット１１０へ撮像を指示する指示信号を送信する場合には、予め設定された撮像周期（例えば１秒）となるように、タイミングを調整する。
【0054】
計測部２３１は、ステップＳ１０６へ進むと、累積スコアを確定させて、雌豚１０１の識別情報および日付情報と共に、データ蓄積部２５０へ記録する。演算部２３０は、ステップＳ１０７で、一連の計測ステップ処理の終了指示が発生しているか否かを確認する。終了指示は、システムオペレータによるメニュー操作や、制御プログラムによる終了判定により発生する。終了指示が発生していなければステップＳ１０１へ戻って一連の処理を継続し、終了指示が発生していれば計測ステップの処理を終了する。
【0055】
図１１は、対象となる雌豚の発情を判定する判定ステップの処理について説明する図である。判定ステップは、システムオペレータや作業員からの要求に応じて開始される。システムオペレータや作業員は、判定装置２００や作業員端末１２０を操作して、発情の有無を知りたい対象となる雌豚を指定する。フローは、判定対象となる雌豚が確定した時点から開始する。
【0056】
判定部２３２は、ステップＳ２０１で、指定された雌豚の直近における一定期間分の累積スコアをデータ蓄積部２５０から読み出す。具体的には、一定期間が例えば４日と設定されている場合には、直近の４日分の累積スコアを読み出す。ステップＳ２０２へ進み、判定部２３２は、発情識別器２６２をメモリ２６０から読み出す。そして、ステップＳ２０３で、データ蓄積部２５０から読み出した一定期間分の累積スコアを、メモリ２６０から読み出した発情識別器２６２へ入力し、判定計算を実行する。判定計算の結果「発情あり」または「発情なし」の出力を得たら、ステップＳ２０４で、当該判定結果を表示モニタ２２０や作業員端末１２０へ出力し、一連の処理を終了する。
【0057】
なお、ここではシステムオペレータや作業員からの要求に応じて判定ステップの処理を実行する例を説明したが、予め定められた条件を満たす場合に自動的に判定ステップの処理を実行しても良い。例えば、判定部２３２は、毎日決められた時間に観察対象のすべての雌豚に対して判定ステップの処理を実行し、「発情あり」と判定された雌豚グループと、「発情なし」と判定された雌豚グループを一覧にして表示モニタ２２０や作業員端末１２０へ出力しても良い。あるいは、日ごとの累積スコアがデータ蓄積部２５０へ記録されるタイミングで自動的に判定ステップの処理を実行し、「発情あり」と判定された場合に限ってシステムオペレータや作業員へ通知するようにしても良い。
【0058】
以上説明した本実施形態においてはＳＶＭである発情識別器２６２を用いて発情を判定したが、発情識別器２６２はＳＶＭに限らない。機械学習によって生成する識別器としては、他にもロジスティック回帰やランダムフォレストの手法も採用し得る。さらには、識別器を用いず、解析的に発情を判定する手法を採用しても良い。
【0059】
図１２は、解析的に発情を判定する発情判定手法を説明する図である。横軸は観察経過日数を表し、縦軸はそれぞれの観察日における累積スコアを表している。ここでは、発情を示す特定の雌豚の１日目から５日目までの累積スコアをプロットして結んだ推移を表す。上述のように、発情している雌豚は一定期間に亘る活動強度の変化がＶ字状となるので、発情しているか否かの判定は、累積スコアの推移がＶ字状になっているか否かを分析すれば良い。
【0060】
そこで、一定期間（図の場合は５日間）の中間日（図の場合は２日目から４日目）に、累積スコアの最小値が存在するか否かを確認する。最小値が中間日に存在しない場合は、「発情なし」と判定する。最小値が中間日に存在する場合は、その最小値が計測された観察日よりも前を前期とし、後を後期とする。図の場合は、３日目の累積スコアが最小値であるので、１日目と２日目が前期であり、４日目と５日目が後期である。
【0061】
そして、前期の累積スコアの最大値である前期最大値と後期の累積スコアの最大値である後期最大値を決定する。図の場合は、１日目の累積スコアが前期最大値であり、５日目の累積スコアが後期最大値である。この前期最大値→最小値→後期最大値の変化が、「発情あり」と判定される程度のＶ字状であるか否かを分析する。具体的には、前期最大値から最小値への低下率Ｔαが、予め設定された基準低下率Ｔα0よりも大きく、最小値から後期最大値への上昇率Ｔβが、予め設定された基準上昇率Ｔβ0よりも大きい場合に、「発情あり」と判定する。そうでない場合には、「発情なし」と判定する。なお、基準低下率Ｔα0および基準上昇率Ｔβ0は、結果が既知の確定データを用いて設定される。
【0062】
以上いくつかの判定手法を説明したが、判定装置においていずれの判定手法を採用するかは、観察対象となる雌豚の頭数や品種、判定装置の規模、コスト、要求される判定精度等を踏まえて決定され得る。また、判定手法に限らず、起臥状態の計測手法も上記の手法に限らない。例えば、カメラユニット１１０の代わりに、ストール１０２に収容された雌豚１０１の頭部付近と臀部付近の高さを計測する距離センサを天井に設置し、計測部２３１は、当該距離センサの出力から起臥状態を認定しても良い。
【0063】
また、本実施形態においては、起立状態と臥位状態に加えて着座状態も起臥の頻度を計測する対象とし、着座状態に対してスコアα＝０．５を与えた。しかし、スコアの数値は雌豚の品種や体重等を考慮して０．５以外の値に調整しても良く、また、着座状態も更に姿勢に応じて細分化して、それぞれに異なる値のスコアを与えても良い。逆に、脚が短い品種のような場合や、それほどの判定精度が必要でない場合には、着座状態を省いて起立状態と臥位状態の２つの状態のみで起臥の頻度を計測しても良い。
【符号の説明】
【0064】
１０１ 雌豚、１０２ ストール、１１０ カメラユニット、１２０ 作業員端末、２００ 判定装置、２１０ サーバ、２２０ 表示モニタ、２３０ 演算部、２３１ 計測部、２３２ 判定部、２４０ 画像処理部、２５０ データ蓄積部、２６０ メモリ、２６１ 起臥識別器、２６２ 発情識別器、２７０ 通信ユニット、９００ インターネット
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】