(57)【要約】【課題】家畜の行動を特定することができる。【解決手段】家畜の行動を特定する行動特定装置であって、前記家畜に装着された加速度センサが測定した加速度データを記憶する記憶手段と、前記記憶手段に記憶されている加速度データのうち、所定の時間の間における１以上の加速度データを取得する取得手段と、前記取得手段により取得された１以上の加速度データの標準偏差と最大値との組に基づき、前記家畜の行動を特定する特定手段と、を有する。【選択図】図１
2020036627.tif 000002
【特許請求の範囲】
【請求項１】
家畜の行動を特定する行動特定装置であって、
前記家畜に装着された加速度センサが測定した加速度データを記憶する記憶手段と、
前記記憶手段に記憶されている加速度データのうち、所定の時間の間における１以上の加速度データを取得する取得手段と、
前記取得手段により取得された１以上の加速度データの標準偏差と最大値との組に基づき、前記家畜の行動を特定する特定手段と、
を有する行動特定装置。
【請求項２】
前記記憶手段は、
更に、前記家畜に装着された電波センサが測定した電波強度データを記憶し、
前記特定手段は、
前記特定された家畜の行動が所定の行動である場合、前記電波強度データに基づき、前記家畜の詳細な行動を特定する、請求項１に記載の行動特定装置。
【請求項３】
前記行動には、前記家畜が立っていることを示す起立行動と、前記家畜が横臥していることを示す横臥行動とが含まれ、
前記記憶手段は、
更に、前記家畜に装着された気圧センサが測定した気圧データを記憶し、
前記特定手段は、
前記気圧データに基づき、前記家畜の前記起立行動又は前記横臥行動を特定する、請求項１又は２に記載の行動特定装置。
【請求項４】
前記特定手段は、
前記加速度データの変化パターン又は前記気圧データの変化パターンに基づいて、前記家畜の前記起立行動又は前記横臥行動を特定する、請求項３に記載の行動特定装置。
【請求項５】
前記家畜は、牛である、請求項１乃至４の何れか一項に記載の行動特定装置。
【請求項６】
家畜の行動を特定する行動特定装置が、
前記家畜に装着された加速度センサが測定した加速度データを記憶部に記憶させる記憶手順と、
前記記憶部に記憶されている加速度データのうち、所定の時間の間における１以上の加速度データを取得する取得手順と、
前記取得手順で取得された１以上の加速度データの標準偏差と最大値との組に基づき、前記家畜の行動を特定する特定手順と、
を実行する行動特定方法。
【請求項７】
コンピュータを、請求項１乃至５の何れか一項に記載の行動特定装置における各手段として機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、行動特定装置、行動特定方法、及びプログラムに関する。
【背景技術】
【０００２】
乳用牛の搾乳期間は、乳用牛を交配（人工授精）させて子牛を生ませた後（分娩後）であり、この期間は泌乳期と呼ばれている。泌乳期は約３０５日程度であることが知られている。また、次の出産に備えて、次回の分娩予定日の約６０日前で搾乳を止めるのが一般的である。この６０日の期間は乾乳期と呼ばれている。乳用牛は、分娩、泌乳期、乾乳期を毎年繰り返している。
【０００３】
乳用牛を飼育する酪農家等は、乳用牛から継続的に搾乳するため、毎年乾乳期に乳用牛を交配（人工授精）させる必要がある。乳用牛の交配（人工授精）にあたっては、乳用牛の発情を発見する必要がある。乳用牛の発情を発見するには、乳用牛の外部的兆候を確認したり、発情時に多く見られる行動を確認したりすることが行われている。また、台帳を用いた繁殖記録による発情時期の管理等も行われている。
【先行技術文献】
【非特許文献】
【０００４】
【非特許文献１】阿部 亮著、「農学基礎セミナー  家畜飼育の基礎 」、新版、社団法人  農産漁村文化協会、2008 年 4 月, p.109, p.122-124.
【発明の概要】
【発明が解決しようとする課題】
【０００５】
ここで、例えば、飼養頭数が数千頭以上にも及ぶ大規模な酪農場では、１頭１頭の乳用牛の外部的兆候や行動を確認することは酪農家の大きな負担となっていた。これに対して、例えば、飼養している乳用牛の行動を特定することで、発情時に多く見られる行動の確認を容易することができれば、酪農家の負担を軽減することができる。
【０００６】
また、乳用牛の行動を特定することで、発情の発見のみならず、病気やケガ等による異常行動の発見にも繋がり、乳用牛の健康管理にも資することができる。
【０００７】
本発明の一実施形態は、上記の点に鑑みてなされたもので、家畜の行動を特定することを目的とする。
【課題を解決するための手段】
【０００８】
上記課題を解決するため、家畜の行動を特定する行動特定装置であって、前記家畜に装着された加速度センサが測定した加速度データを記憶する記憶手段と、前記記憶手段に記憶されている加速度データのうち、所定の時間の間における１以上の加速度データを取得する取得手段と、前記取得手段により取得された１以上の加速度データの標準偏差と最大値との組に基づき、前記家畜の行動を特定する特定手段と、を有する。
【発明の効果】
【０００９】
家畜の行動を特定することができる。
【図面の簡単な説明】
【００１０】
【図１】第一の実施形態に係る行動特定システムの全体構成の一例を示す図である。
【図２】牛の行動特定の一例を説明する図である。
【図３】第一の実施形態に係る測定データ記憶部に記憶されている測定データの一例を示す図である。
【図４】行動特定モデルの一例を示す図である。
【図５】第一の実施形態に係る行動特定処理部の機能構成の一例を示す図である。
【図６】第一の実施形態に係る行動特定処理の一例を示すフローチャートである。
【図７】第二の実施形態に係る行動特定システムの全体構成の一例を示す図である。
【図８】気圧分析により起立又は横臥を特定する場合の一例を示す図である。
【図９】第二の実施形態に係る測定データ記憶部に記憶されている測定データの一例を示す図である。
【図１０】第二の実施形態に係る行動特定処理部の機能構成の一例を示す図である。
【図１１】第二の実施形態に係る行動（起立又は横臥）特定処理の一例を示すフローチャートである。
【図１２】第三の実施形態に係る測定データ及び基準気圧データの取得範囲の一例を説明する図である。
【発明を実施するための形態】
【００１１】
以下、本発明の実施の形態について、図面を参照しながら説明する。以降では、家畜の一例として牛の行動を特定する場合について説明する。ただし、家畜は、牛に限られない。
【００１２】
［第一の実施形態］
＜全体構成＞
まず、本実施形態に係る行動特定システム１の全体構成について、図１を参照しながら説明する。図１は、第一の実施形態に係る行動特定システム１の全体構成の一例を示す図である。
【００１３】
図１に示すように、本実施形態に係る行動特定システム１は、牛の行動を特定する行動特定装置１０と、牛に装着された１以上のタグ２０と、所定の電波を周囲に発信する1以上の漏洩同軸アンテナ３０とが含まれる。なお、タグ２０は、牛の首部分に固定して装着されることが好ましい。
【００１４】
タグ２０は、牛に装着される機器である。１頭の牛に対して１つのタグ２０が装着されている。タグ２０には、当該タグ２０を装着した牛の加速度（Ｘ軸、Ｙ軸、及びＺ軸の３軸の加速度）を測定する加速度センサと、漏洩同軸アンテナ３０から受信する電波の受信強度（ＲＳＳＩ：Received Signal Strength Indicator）を測定する電波センサとが含まれる。
【００１５】
タグ２０は、所定の時間毎（例えば２秒毎）に、加速度センサにより測定した加速度センサ値と、電波センサにより測定したＲＳＳＩ値とを含む測定データを行動特定装置１０に送信する。行動特定装置１０に送信された測定データは、後述する測定データ記憶部２００に蓄積（記憶）される。
【００１６】
漏洩同軸アンテナ３０は、牛舎内の餌場や水場等の所定の場所に設置されるアンテナである。漏洩同軸アンテナ３０は、漏洩同軸ケーブルにより周囲に所定の電波を発信する。なお、以降では、餌場に設置されている漏洩同軸アンテナ３０と、水場に設置されている漏洩同軸アンテナ３０とを区別する場合には、それぞれ「漏洩同軸アンテナ３０１」と「漏洩同軸アンテナ３０２」と表す。
【００１７】
行動特定装置１０は、牛の行動を特定する１以上のコンピュータである。行動特定装置１０は、行動特定処理部１００と、測定データ記憶部２００と、行動特定モデル３００とを有する。
【００１８】
行動特定処理部１００は、測定データ記憶部２００に記憶されている測定データと、行動特定モデル３００とに基づいて、牛の行動を特定する。行動特定処理部１００は、行動特定装置１０にインストールされた１以上のプログラムが、ＣＰＵ（Central Processing Unit）等に実行させる処理により実現される。
【００１９】
測定データ記憶部２００は、タグ２０から受信した測定データを記憶する。測定データ記憶部２００は、例えばＨＤＤ（Hard Disk Drive）やＳＳＤ（Solid State Drive）等の補助記憶装置等を用いて実現可能である。
【００２０】
行動特定モデル３００は、測定データに含まれる加速度センサ値から行動を特定するためのモデルである。行動特定モデル３００は、例えば、家畜の種類や特定したい行動の種類に応じて、ＳＶＭ（Support Vector Machine）等の機械学習の手法により予め作成される。行動特定モデル３００は、例えばＨＤＤやＳＳＤ等の補助記憶装置等に記憶されている。
【００２１】
なお、図１に示す行動特定システム１の構成は一例であって、他の構成であっても良い。例えば、行動特定装置１０は、複数台のコンピュータで構成されていても良い。また、例えば、行動特定処理部１００が有する機能の一部を、行動特定装置１０とネットワークを介して接続される装置（クラウドサーバ等）が有していても良い。更に、例えば、牛には、タグ２０の代わりに、加速度センサと、電波センサとがそれぞれ装着されていても良い。
【００２２】
＜牛の行動の特定＞
ここで、本実施形態に係る行動特定システム１により特定される牛の行動と、特定方法の概要とについて、図２を参照しながら説明する。図２は、牛の行動特定の一例を説明する図である。
【００２３】
まず、本実施形態に係る行動特定システム１が特定する牛の行動には、動作が小さい順に、「起立」、「横臥」、「反芻」、及び「活動」の４つがあるものとする。また、このうち「活動」は、更に、「採食」、「歩行」、及び「飲水」の３つに詳細化されるものとする。したがって、本実施形態では、「起立」、「横臥」、「反芻」、「採食」、「歩行」、及び「飲水」の６つの行動のうち、いずれの行動を牛が行っているかを特定するものとする。
【００２４】
「起立」とは、牛が立っている状態のことである。「横臥」とは、牛が横たわっている状態のことである。「反芻」とは、牛が反芻（一度飲み込んだ食物を口に戻して咀嚼し直す動作）を行っている状態のことである。「活動」とは、「採食」、「歩行」、又は「飲水」のいずれかを牛が行っている状態のことである。「採食」とは、牛が食物を食べている状態のことである。「歩行」とは、牛が歩いている状態のことである。「飲水」とは、牛が水を飲んでいる状態のことである。
【００２５】
「起立」、「横臥」、「反芻」、及び「活動」のいずれの行動を牛が行っているかは、動作強度分析により特定される。動作強度分析では、測定データに含まれる加速度センサ値と、行動特定モデル３００とに基づいて、牛の行動が、「起立」、「横臥」、「反芻」、及び「活動」のいずれであるかを特定する。より具体的には、動作強度分析では、所定の時間（例えば１０分）の間における加速度センサ値から所定の指標値（所定の時間の間における加速度センサ値のＬ２ノルムの振れ幅の最大値及び加速度センサ値のＬ２ノルムの標準偏差）を算出する。なお、Ｌ２ノルムの振れ幅の最大値とは、所定の時間の間における加速度センサ値のＬ２ノルムの平均値と、当該所定の時間の間における各速度センサ値のＬ２ノルムの値との差の絶対値が最大となる値である。
【００２６】
そして、「起立」、「横臥」、「反芻」、及び「活動」の４つの領域に分類される行動特定モデル３００において、算出した指標値を示す点がこれら４つの領域のいずれに含まれるかにより行動を特定する。
【００２７】
また、動作強度分析により「活動」と特定された場合に、「採食」、「歩行」、及び「飲水」のいずれの行動を牛が行っているかは、位置分析により特定される。位置分析では、測定データに含まれるＲＳＳＩ値と、所定の閾値とに基づいて、牛の行動が、「採食」、「歩行」、及び「飲水」のいずれの行動であるかを特定する。より具体的には、位置分析では、餌場に設置された漏洩同軸アンテナ３０１から受信した電波のＲＳＳＩ値が、所定の閾値を超えている場合（すなわち、牛が餌場にいる場合）、「採食」と特定する。同様に、水場に設置された漏洩同軸アンテナ３０２から受信した電波のＲＳＳＩ値が、所定の閾値を超えている場合（すなわち、牛が水場にいる場合）、「飲水」と特定する。一方で、これらのいずれでも無い場合（すなわち、餌場にも水場にもいない場合）、「歩行」と特定する。
【００２８】
なお、例えば、ＧＰＳ（Global Positioning System）受信機等がタグ２０に含まれている場合は、位置分析は、ＧＰＳ受信機から取得した位置情報を用いて行われても良い。この場合、タグ２０から取得される位置情報から、牛が餌場にいるのか、水場にいるのか、又は餌場にも水場にもいないのかを特定すれば良い。
【００２９】
以上のように、本実施形態に係る行動特定システム１は、牛に装着されたタグ２０から受信した加速度センサ値から、「起立」、「横臥」、「反芻」、及び「活動」の４つの行動を特定する。また、本実施形態に係る行動特定システム１は、牛の行動が「活動」であると特定された場合に、牛に装着されたタグ２０から受信したＲＳＳＩ値から、更に、「採食」、「歩行」、及び「飲水」の３つの行動に詳細化する。これにより、牛の行動が、「起立」、「横臥」、「反芻」、「採食」、「歩行」、及び「飲水」のいずれであるかを特定することができる。
【００３０】
したがって、これら特定した牛の行動を、例えば、行動特定装置１０の表示装置（ディスプレイ等）に表示することで、発情時に多く見られる行動の確認を容易に行うことができるようになる。また、病気やケガ等による異常行動の早期発見にも繋がり、牛の健康管理も容易に行うことができるようになる。
【００３１】
＜測定データ記憶部２００に記憶されている測定データ＞
ここで、本実施形態に係る測定データ記憶部２００に記憶されている測定データについて、図３を参照しながら説明する。図３は、第一の実施形態に係る測定データ記憶部２００に記憶されている測定データの一例を示す図である。なお、行動特定装置１０は、タグ２０から測定データを受信した場合、行動特定処理部１００により、受信した測定データを測定データ記憶部２００に記憶（蓄積）させれば良い。
【００３２】
図３に示すように、測定データ記憶部２００には、タグを識別するタグＩＤ毎に、１以上の測定データが記憶されている。なお、１頭の牛に対して１つのタグ２０が装着されていることから、タグＩＤは、牛を識別する情報（牛の個体識別情報）であっても良い。
【００３３】
測定データには、日時と、加速度センサ値と、ＲＳＳＩ値とが含まれる。日時は、例えば、タグ２０が測定データを送信した日時である。なお、日時は、行動特定装置１０が測定データを受信した日時であっても良い。
【００３４】
加速度センサ値は、タグ２０に含まれる加速度センサにより測定された加速度の値である。加速度センサ値には、Ｘ軸方向の加速度成分を示すＸ成分と、Ｙ軸方向の加速度成分を示すＹ成分と、Ｚ軸方向の加速度成分を示すＺ成分とが含まれる。例えば、日時「ｔ１」の測定データには、加速度センサ値のＸ成分「Ｘ１」と、Ｙ成分「Ｙ１」と、Ｚ成分「Ｚ１」とが含まれる。
【００３５】
ＲＳＳＩ値は、タグ２０に含まれる電波センサにより測定されたＲＳＳＩの値である。ＲＳＳＩ値には、例えば、餌場に設置されている漏洩同軸アンテナ３０１から受信した電波のＲＳＳＩ値と、水場に設置されている漏洩同軸アンテナ３０２から受信した電波のＲＳＳＩ値とが含まれる。例えば、日時「ｔ１」の測定データには、餌場に設置されている漏洩同軸アンテナ３０１から受信したＲＳＳＩ値「ｒ１１」と、水場に設置されている漏洩同軸アンテナ３０２から受信したＲＳＳＩ値「ｒ１２」とが含まれる。
【００３６】
このように、本実施形態に係る測定データ記憶部２００に記憶されている測定データには、タグＩＤ毎に、日時と、加速度センサ値と、ＲＳＳＩ値とが含まれる測定データが蓄積（記憶）されている。
【００３７】
＜行動特定モデル３００＞
ここで、測定データに含まれる加速度センサ値から行動を特定するための行動特定モデル３００について、図４を参照しながら説明する。図４は、行動特定モデル３００の一例を示す図である。
【００３８】
図４に示すように、行動特定モデル３００は、所定の時間（例えば１０分）の間における加速度センサ値のＬ２ノルムの振れ幅の最大値を横軸、当該Ｌ２ノルムの標準偏差を縦軸とした場合に、これら振れ幅の最大値及び標準偏差の値と、行動との関係を示す関係グラフとして表される。
【００３９】
例えば、ある１０分間における加速度センサ値のＬ２ノルムの振れ幅の最大値と標準偏差とが領域Ｄ１に含まれる場合、当該１０分間における牛の行動は「起立」と特定される。また、例えば、ある１０分間における加速度センサ値のＬ２ノルムの振れ幅の最大値と標準偏差とが領域Ｄ２に含まれる場合、当該１０分間における牛の行動は「横臥」と特定される。
【００４０】
同様に、例えば、ある１０分間における加速度センサ値のＬ２ノルムの振れ幅の最大値と標準偏差とが領域Ｄ３に含まれる場合、当該１０分間における牛の行動は「反芻」と特定される。なお、「反芻」と特定される領域Ｄ３は、更に、起立した状態で反芻を行う「反芻（起立）」と、横臥した状態で反芻を行う「反芻（横臥）」とに分けられても良い。
【００４１】
また、同様に、例えば、ある１０分間における加速度センサ値のＬ２ノルムの振れ幅の最大値と標準偏差とが領域Ｄ４に含まれる場合、当該１０分間における牛の行動は「活動」と特定される。
【００４２】
なお、図４に示す例では、行動特定モデル３００の各領域Ｄ１〜Ｄ４が互いに重なっていない場合を示しているが、各領域Ｄ１〜Ｄ４のうちの２以上の領域が互いに重なっている部分が存在しても良い。
【００４３】
このように、行動特定モデル３００は、所定の時間（例えば１０分）の間における加速度センサ値のＬ２ノルムの最大値及び標準偏差と、牛の行動との関係を表す領域が定義されたモデルである。このような行動特定モデル３００は、ＳＶＭ等の機械学習の手法により予め作成される。なお、ＳＶＭは一例であって、例えば、ニューラルネットワーク等の種々の機械学習の手法により作成されても良い。また、例えば、牛が同一の動作を行った場合であっても、加速度センサの種類等によって、具体的な加速度センサ値は異なる。ただし、各行動をそれぞれ示す領域Ｄ１〜Ｄ４の位置関係は、略一定に定まっている。
【００４４】
なお、以降では、Ｌ２ノルムの振れ幅の最大値を、単に、「Ｌ２ノルムの最大値」又は「最大値」とも表す。
【００４５】
＜行動特定処理部１００の機能構成＞
次に、本実施形態に係る行動特定処理部１００の機能構成について、図５を参照しながら説明する。図５は、第一の実施形態に係る行動特定処理部１００の機能構成の一例を示す図である。
【００４６】
図５に示すように、行動特定処理部１００は、取得部１０１と、前処理部１０２と、指標値算出部１０３と、行動特定部１０４とを有する。
【００４７】
取得部１０１は、測定データ記憶部２００から測定データを取得する。このとき、取得部１０１は、例えば、タグＩＤ毎に、所定の時間（例えば１０分）の間の測定データを測定データ記憶部２００から取得する。
【００４８】
前処理部１０２は、取得部１０１により取得された測定データに対して前処理を行う。前処理とは、例えば、測定データの欠損補完（リサンプリング）処理やノイズ除去処理等である。前処理部１０２は、欠損補完後の測定データに含まれる加速度センサ値のＬ２ノルムを算出して、算出したＬ２ノルムを用いてノイズ除去処理を行う。
【００４９】
指標値算出部１０３は、前処理部１０２による前処理後のＬ２ノルムから指標値を算出する。例えば、指標値算出部１０３は、指標値として、前処理後のＬ２ノルムから最大値及び標準偏差を算出する。
【００５０】
行動特定部１０４は、指標値算出部１０３により算出された最大値及び標準偏差と、行動特定モデル３００とから４つの行動（「起立」、「横臥」、「反芻」、及び「活動」）を特定する。また、行動特定部１０４は、更に、上記で「活動」が特定された場合に、測定データに含まれるＲＳＳＩ値から詳細な行動（「採食」、「歩行」、又は「飲水」）を特定する。なお、行動特定部１０４は、当該最大値及び標準偏差と、行動特定モデル３００と同等のデータを生成するプログラム等の処理結果とから上記４つの行動を特定しても良い。
【００５１】
なお、行動特定部１０４により特定された行動を示す情報は、例えば、ＨＤＤやＳＳＤ等の補助記憶装置等で実現されるＤＢ（データベース）やファイル等に格納されても良いし、ディスプレイ等の表示装置に表示されても良い。また、行動特定装置１０に接続される他の装置（例えば、ＰＣやスマートフォン、タブレット端末等）に出力されても良い。
【００５２】
＜行動特定処理＞
次に、本実施形態に係る行動特定処理について、図６を参照しながら説明する。図６は、第一の実施形態に係る行動特定処理の一例を示すフローチャートである。なお、図６に示す行動特定処理は、例えば、予め設定された日時に実行されても良いし、所定の時間（例えば２４時間）毎に繰り返し実行されても良い。
【００５３】
まず、取得部１０１は、タグＩＤ毎に、所定の時間（例えば１０分）の間の測定データを測定データ記憶部２００から取得する（ステップＳ１１）。このように、取得部１０１は、牛（タグＩＤ）毎に、所定の時間単位（例えば１０分単位）の測定データを測定データ記憶部２００から取得する。なお、このような所定の時間は、１０分に限られず、例えば行動特定装置１０のユーザが任意の時間に設定することができる。
【００５４】
以降では、タグＩＤ「Ｔａｇ１」の１０分間の測定データ１，測定データ２，・・・，測定データＮ１が取得部１０１により取得された場合について説明を続ける。
【００５５】
次に、前処理部１０２は、取得部１０１により取得された測定データに対して前処理を行う（ステップＳ１２）。すなわち、前処理部１０２は、測定データの欠損補完（リサンプリング）処理やノイズ除去処理等を行う。なお、欠損補完処理は、収集データの精度を上げるため、データが取得できなかった場合の欠損を補完する処理である。また、ノイズ除去処理は、牛の瞬間的な動作（例えば、瞬間的に身体を震わせる動作や瞬間的に大きく身体をびくつかせる動作等）を示すデータを除去する処理である。
【００５６】
次に、指標値算出部１０３は、前処理部１０２による前処理後のＬ２ノルムから指標値（Ｌ２ノルムの最大値及び標準偏差）を算出する（ステップＳ１３）。すなわち、前処理後のＬ２ノルムａ１，ａ２，・・・ａＮ３の標準偏差σ及び最大値ｍを算出する。なお、最大値ｍは、Ｌ２ノルムａ１，ａ２，・・・ａＮ３の平均値と、各Ｌ２ノルムａｉ（ｉ＝１，・・・，Ｎ３）との差の最大値である。
【００５７】
次に、行動特定部１０４は、指標値算出部１０３により算出された標準偏差σ及び最大値ｍと、行動特定モデル３００とから４つの行動（「起立」、「横臥」、「反芻」、及び「活動」）を特定する（ステップＳ１４）。すなわち、行動特定部１０４は、指標値算出部１０３により算出された標準偏差σ及び最大値ｍが、行動特定モデル３００上の領域Ｄ１〜Ｄ４のいずれの領域に含まれるかを特定することで、行動を特定する。
【００５８】
次に、行動特定部１０４は、上記のステップＳ１４で特定した行動が「活動」であるか否かを判定する（ステップＳ１５）。
【００５９】
ステップＳ１５において、特定した行動が「活動」でないと判定された場合（すなわち、特定した行動が「起立」、「横臥」、又は「反芻」である場合）、行動特定処理部１００は、処理を終了させる。
【００６０】
一方、ステップＳ１５において、特定した行動が「活動」であると判定された場合、行動特定部１０４は、上記のステップＳ１２の欠損補完後の測定データ１，測定データ２，・・・，測定データＮ２に含まれるＲＳＳＩ値から詳細な行動（「採食」、「歩行」、又は「飲水」）を特定する（ステップＳ１６）。
【００６１】
すなわち、例えば、測定データ１，測定データ２，・・・，測定データＮ２にそれぞれ含まれるＲＳＳＩ値のうち、漏洩同軸アンテナ３０１や漏洩同軸アンテナ３０２から受信した電波のＲＳＳＩ値が所定の閾値を超えている場合、行動特定部１０４は、詳細な行動を「採食」や「飲水」と特定する。一方で、漏洩同軸アンテナ３０１から受信した電波のＲＳＳＩ値と、漏洩同軸アンテナ３０２から受信した電波のＲＳＳＩ値とがいずれも所定の閾値を超えていない場合、行動特定部１０４は、詳細な行動を「歩行」と特定する。
【００６２】
以上により、本実施形態に係る行動特定システム１は、牛が装着するタグ２０から受信した測定データから当該牛の行動を特定することができる。このとき、本実施形態に係る行動特定システム１では、測定データに含まれる加速度センサ値から算出された指標値と、機械学習の手法により予め作成された行動特定モデル３００とを用いることで、牛の行動を特定する。特に、指標値として、Ｌ２ノルムの標準偏差と最大値とを用いることで、牛に特有の行動である「反芻」を高い精度で特定することができる。言い換えれば、牛に特有の行動である「反芻」を、「起立」や「横臥」等の行動と高い精度で区別することができるようになる。
【００６３】
このように特定された牛の行動を、例えば行動特定装置１０の表示装置（ディスプレイ等）に表示して、酪農家等に提供することで、牛の発情時に多く見られる行動の確認を容易に行うことができるようになる。また、病気やケガ等による異常行動の早期発見にも繋がり、牛の健康管理も容易に行うことができるようになる。
【００６４】
［第二の実施形態］
次に、第二の実施形態について説明する。第一の実施形態では、加速度センサ値を用いて牛の４つの行動（「起立」、「横臥」、「反芻」、及び「活動」）を特定したが、これらのうち「起立」と「横臥」とは区別が難しい場合がある。すなわち、例えば、牛が「起立」している場合でも小刻みに身体を動かしている場合がある一方で、牛が「横臥」している場合でも大きく身体を動かしている場合がある。このような場合には、上記の動作強度分析では、「起立」と「横臥」との区別が難しい場合がある（言い換えれば、「起立と「横臥」とが誤って特定される場合がある。）。
【００６５】
同様に、例えば、牛が「反芻」している場合でも、起立した状態で反芻を行っている場合もあれば、横臥した状態で反芻を行っている場合もある。このような場合に、上記の動作強度分析では、同様に、「起立」と「横臥」との区別が難しい場合もある。
【００６６】
そこで、第二の実施形態では、気圧センサを用いることで、「起立」と「横臥」とをより高い精度で特定する場合について説明する。
【００６７】
なお、第二の実施形態では、主に、第一の実施形態との相違点について説明し、第一の実施形態と同様の機能を有する箇所や同様の処理を行う箇所については、適宜、その説明を省略するものとする。
【００６８】
＜全体構成＞
まず、本実施形態に係る行動特定システム１の全体構成について、図７を参照しながら説明する。図７は、第二の実施形態に係る行動特定システム１の全体構成の一例を示す図である。
【００６９】
図７に示すように、本実施形態に係る行動特定システム１は、基準となる気圧を測定する基準気圧センサ４０が含まれる。
【００７０】
基準気圧センサ４０は、牛舎内の所定の位置（例えば、牛舎内の地面上）に設置され、基準となる気圧を測定する。基準気圧センサ４０は、所定の時間毎（例えば２秒毎）に、測定した気圧を示す基準気圧センサ値を含む基準気圧データを行動特定装置１０に送信する。行動特定装置１０に送信された基準気圧データは、後述する基準気圧データ記憶部４００に蓄積（記憶）される。
【００７１】
また、本実施形態に係るタグ２０には、更に、気圧を測定する気圧センサが含まれる。タグ２０は、所定の時間毎（例えば２秒毎）に、更に、気圧センサ値を含む測定データを行動特定装置１０に送信する。なお、タグ２０には気圧センサが必ずしも含まれている必要はなく、例えば、気圧センサがタグ２０とは別体で牛に装着されていても良い。
【００７２】
更に、本実施形態に係る行動特定装置１０は、基準気圧データ記憶部４００を有する。基準気圧データ記憶部４００は、基準気圧センサ４０から受信した基準気圧データを記憶する。基準気圧データ記憶部４００は、例えばＨＤＤやＳＳＤ等の補助記憶装置等を用いて実現可能である。
【００７３】
本実施形態に係る行動特定処理部１００は、更に、測定データ記憶部２００に記憶されている測定データに含まれる気圧センサ値と、基準気圧データ記憶部４００に記憶されている基準気圧データに含まれる基準気圧センサ値とに基づいて、牛の行動（「起立」及び「横臥」）を特定する。
【００７４】
＜気圧センサによる「起立」及び「横臥」の特定＞
ここで、本実施形態に係る行動特定システムにより特定される牛の行動のうち、気圧センサによる「起立」及び「横臥」の特定方法の概要について、図８を参照しながら説明する。図８は、気圧分析により起立又は横臥を特定する場合の一例を示す図である。
【００７５】
「起立」及び「横臥」のいずれの行動を牛が行っているかは、第一の実施形態で説明した動作強度分析に加えて、気圧分析でも特定される。気圧分析では、測定データに含まれる気圧センサ値と、基準気圧データに含まれる基準気圧センサ値とに基づいて、牛の行動が、「起立」又は「横臥」のいずれであるかを特定する。より具体的には、気圧分析では、所定の時間（例えば２４時間）の間における各気圧センサ値と各基準センサ値との差分を算出する。そして、算出した差分（この差分を「差分気圧センサ値」と表す。）が、所定の閾値を超えている場合、「横臥」と特定する。一方で、差分気圧センサ値が、所定の閾値を超えていない場合、「起立」と特定する。基準気圧センサ４０により測定された差分気圧センサ値を用いることで、例えば、天候（低気圧の接近等）に依らずに、牛に装着されたタグ２０に含まれる気圧センサの測定値が相対的に上がったのか又は下がったのかを判定することができるようになる。
【００７６】
このように気圧センサにより測定された気圧センサ値を用いて、差分気圧センサ値が閾値を超えた場合に「横臥」、閾値を超えていない場合に「起立」と特定することで、牛の行動のうち、「起立」と「横臥」とを高い精度で特定することができる。
【００７７】
＜測定データ記憶部２００に記憶されている測定データ＞
ここで、本実施形態に係る測定データ記憶部２００に記憶されている測定データについて、図９を参照しながら説明する。図９は、第二の実施形態に係る測定データ記憶部２００に記憶されている測定データの一例を示す図である。なお、行動特定装置１０は、タグ２０から測定データを受信した場合、行動特定処理部１００により、受信した測定データを測定データ記憶部２００に記憶（蓄積）させれば良い。
【００７８】
図９に示すように、測定データ記憶部２００には、タグを識別するタグＩＤ毎に、１以上の測定データが記憶されている。本実施形態に係る測定データには、更に、気圧センサ値が含まれる。気圧センサ値は、タグ２０に含まれる気圧センサにより測定された気圧の値である。
【００７９】
このように、本実施形態に係る測定データ記憶部２００に記憶されている測定データには、タグＩＤ毎に、日時と、加速度センサ値と、ＲＳＳＩ値と、気圧センサ値とが含まれる測定データが蓄積（記憶）されている。
【００８０】
＜行動特定処理部１００の機能構成＞
次に、本実施形態に係る行動特定処理部１００の機能構成について、図１０を参照しながら説明する。図１０は、第二の実施形態に係る行動特定処理部１００の機能構成の一例を示す図である。
【００８１】
図１０に示すように、本実施形態に係る行動特定処理部１００は、更に、差分値算出部１０５を有する。また、本実施形態に係る取得部１０１及び行動特定部１０４は、第一の実施形態と機能が異なる。
【００８２】
本実施形態に係る取得部１０１は、位置分析による行動特定が行われる場合、タグＩＤ毎に、所定の時間（例えば２４時間）の間の測定データを測定データ記憶部２００から取得すると共に、所定の時間の間の基準気圧データを基準気圧データ記憶部４００から取得する。
【００８３】
差分値算出部１０５は、取得部１０１により取得された測定データに含まれる気圧センサ値と、基準気圧データに含まれる基準気圧センサ値との差分を示す差分気圧センサ値を算出する。
【００８４】
なお、前処理部１０２は、取得部１０１により取得された測定データに対して前処理を行っても良い。この場合、前処理部１０２は、測定データの欠損補完処理やノイズ除去処理（すなわち、気圧センサ値のノイズ除去）等を行えば良い。そして、この場合、差分値算出部１０５は、前処理部１０２による前処理後の測定データに含まれる気圧センサ値と、基準気圧データに含まれる基準気圧センサ値とから差分気圧センサ値を算出すれば良い。
【００８５】
行動特定部１０４は、差分値算出部１０５により算出された差分気圧センサ値が所定の閾値を超えているか否かを判定することで、牛の行動が「起立」又は「横臥」のいずれであるかを特定する。なお、所定の閾値としては、例えば、取得部１０１により取得された測定データに含まれる気圧センサ値から算出される差分気圧センサ値の最大値と最小値との中間値とすれば良い。これは、例えば、２４時間の間に、牛は、少なくとも１度は「起立」と「横臥」とを行うと考えられるためである。
【００８６】
なお、例えば、タグ２０が牛の首輪等に装着されている場合、牛の「採食」行動や「飲水」行動により気圧センサ値が上がることがある。これは、牛が起立した状態から首を下げて採食や飲水を行うためである。牛の首輪等にタグ２０が装着されている場合、「採食」又は「飲水」行動を行った場合におけるタグ２０の地面からの高さは、牛が「横臥」した場合におけるタグ２０の地面からとほぼ同じである場合が多い。したがって、取得部１０１は、動作強度分析及び位置分析により「採食」又は「飲水」行動が特定された前後の時間帯を含む測定データを取得して、取得した測定データに含まれる気圧センサ値から算出される差分気圧センサ値の最大値と最小値との中間値を閾値とすることもできる。
【００８７】
＜行動（起立又は横臥）特定処理＞
次に、本実施形態に係る行動特定処理について、図１１を参照しながら説明する。図１１は、第二の実施形態に係る行動（起立又は横臥）特定処理の一例を示すフローチャートである。なお、図１１に示す行動特定処理は、例えば、予め設定された日時に実行されても良いし、所定の時間（例えば２４時間）毎に繰り返し実行されても良い。
【００８８】
まず、取得部１０１は、タグＩＤ毎に、所定の時間（例えば２４時間）の間の測定データを測定データ記憶部２００から取得すると共に、所定の時間の間の基準気圧データを基準気圧データ記憶部４００から取得する（ステップＳ２１）。
【００８９】
以降では、タグＩＤ「Ｔａｇ１」の２４時間の間の測定データ１，測定データ２，・・・，測定データＮ４と、２４時間の間の基準気圧データ１，基準気圧データ２，・・・，基準気圧データＮ４が取得部１０１により取得された場合について説明を続ける。なお、簡単のため、測定データｉ及び基準気圧データｉ（ｉ＝１，２，・・・，Ｎ４）は、同一日時におけるデータであるものとして説明する。
【００９０】
次に、差分値算出部１０５は、取得部１０１により取得された測定データに含まれる気圧センサ値と、基準気圧データに含まれる基準気圧センサ値との差分を示す差分気圧センサ値を算出する（ステップＳ２２）。
【００９１】
すなわち、差分値算出部１０５は、測定データ１に含まれる気圧センサ値と、基準気圧データ１に含まれる基準気圧センサ値との差分気圧センサ値１を算出する。また、差分値算出部１０５は、測定データ２に含まれる気圧センサ値と、基準気圧データ２に含まれる基準気圧センサ値との差分気圧センサ値２を算出する。以降も同様に、差分値算出部１０５は、測定データＮ４に含まれる気圧センサ値と、基準気圧データＮ４に含まれる基準気圧センサ値との差分気圧センサ値Ｎ４を算出する。
【００９２】
次に、行動特定部１０４は、差分値算出部１０５により算出された差分気圧センサ値が所定の閾値を超えているか否かを判定することで、牛の行動が「起立」又は「横臥」のいずれであるかを特定する（ステップＳ２３）。すなわち、行動特定部１０４は、差分気圧センサ値１，差分気圧センサ値２，・・・，差分気圧センサ値Ｎ４が所定の閾値を超えているか否かを判定することで、牛の行動が「起立」又は「横臥」のいずれであるかを特定する。
【００９３】
例えば、ある差分気圧センサ値が所定の閾値を超えている場合、行動特定部１０４は、当該差分気圧センサ値が算出された日時における牛の行動を「横臥」と特定する。一方で、ある差分気圧センサ値が所定の閾値を超えていない場合、行動特定部１０４は、当該差分気圧センサ値が算出された日時おける牛の行動を「起立」と特定する。
【００９４】
行動特定部１０４は、動作強度分析の結果を考慮して、牛の行動が「起立」又は「横臥」のいずれであるかを考慮しても良い。より具体的には、動作強度分析により「反芻」又は「活動」と行動が特定された日時においては、気圧分析は行わないようにしても良い。一方で、動作強度分析により「起立」又は「横臥」と行動が特定された場合、当該特定された日時を含む所定の範囲の測定データ及び基準気圧データを取得して、気圧分析を行うようにすれば良い。
【００９５】
すなわち、動作強度分析により「反芻」又は「活動」が特定された場合には、気圧分析を行わない一方で、動作強度分析により「起立」又は「横臥」が特定された場合には、気圧分析を行うようにしても良い。これにより、動作強度分析により「起立」又は「横臥」が特定された場合に、より高い精度で「起立」又は「横臥」を特定するための気圧分析が行われる。このように、気圧分析は、動作強度分析の結果に応じて行われるようにしても良い。
【００９６】
以上により、本実施形態に係る行動特定システム１は、気圧センサを用いた気圧分析により、牛の「起立」及び「横臥」行動を特定することができる。しかも、本実施形態に係る行動特定システム１は、気圧分析単独で、又は動作強度分析と併用することで、動作強度分析単独による場合と比較して、高い精度で「起立」及び「横臥」行動を特定することができる。
【００９７】
［第三の実施形態］
次に、第三の実施形態について説明する。第三の実施形態では、牛舎がフリーバーン式である場合について説明する。
【００９８】
牛舎には、一般に、フリーストール式牛舎と、フリーバーン式牛舎とが知られている。フリーバーン式牛舎では、フリーストール式牛舎と異なり、牛舎内に丘が存在する場合が多い。このため、フリーバーン式牛舎では、気圧分析による行動特定が誤る場合がある。
【００９９】
例えば、牛が丘の下の平面上で「起立」している場合における気圧センサ値と、丘の上で「横臥」している場合における気圧センサ値とが同程度であったとする。この場合、丘の上での「横臥」が、「起立」と誤って特定される場合がある。
【０１００】
そこで、第三の実施形態では、測定データ及び基準気圧データの取得範囲を工夫することで、「起立」と「横臥」とをフリーバーン式牛舎においても高い精度で特定する場合について説明する。
【０１０１】
なお、第三の実施形態では、主に、第二の実施形態との相違点について説明し、第二の実施形態と同様の機能を有する箇所や同様の処理を行う箇所については、適宜、その説明を省略するものとする。
【０１０２】
＜フリーバーン牛舎における測定データの取得＞
牛舎がフリーバーン式である場合には、図１１のステップＳ２１で、取得部１０１により取得される測定データと基準気圧データとの範囲を工夫する。これについて、図１２を参照しながら説明する。図１２は、フリーバーン式牛舎における測定データ及び基準気圧データの取得範囲の一例を説明する図である。
【０１０３】
図１２に示すように、時間Ｔ１〜Ｔ２の間は牛が平面（丘の下）上におり、時間Ｔ２〜Ｔ３の間に丘を登って、時間Ｔ３〜Ｔ４の間は牛が丘上にいるものとする。この場合、平面上における牛の行動特定には、牛が平面上にいる時間Ｔ１〜Ｔ２の間の測定データ及び基準気圧データを用いて、気圧分析を行う。一方で、丘上にいる牛の行動特定には、牛が丘上にいる時間Ｔ３〜Ｔ４の間の測定データ及び基準気圧データを用いて、気圧分析を行う。
【０１０４】
平面上にいるか、丘を登っている最中か、又は丘上にいるかは、加速度センサ値と気圧センサ値とを用いて特定することができる。すなわち、例えば、ある時間の間、気圧センサ値が殆ど変化しない一方で、加速度センサ値が変化する場合、牛は、平面上にいると特定することができる。また、例えば、ある時間の間、気圧センサ値と加速度センサ値とが共に変化する場合、牛は、丘を登っている最中と特定することができる。更に、例えば、牛が丘を登っている最中であると特定された後、牛が平面上にいると特定された場合、牛は、丘上にいると特定することができる。なお、丘を下っている最中については、丘を登っている最中と同様に特定することができる。
【０１０５】
なお、牛が丘を登っている最中（又は丘を下っている最中）に「横臥」する場合も有り得る。また、牛が丘を登っている最中（又は丘を下っている最中）に「横臥」した後、その後、「起立」する場合も有り得る。この場合、例えば、丘を登っている最中（又は丘を下っている最中）における加速度センサ値の変化パターン又は気圧センサ値の変化パターンから特定すれば良い。
【０１０６】
以上により、本実施形態に係る行動特定システム１は、フリーバーン式牛舎においても気圧センサを用いた気圧分析により、牛の「起立」及び「横臥」行動を特定することができる。
【０１０７】
本発明は、具体的に開示された上記の実施形態に限定されるものではなく、特許請求の範囲から逸脱することなく、種々の変形や変更が可能である。
【符号の説明】
【０１０８】
１        行動特定システム
１０      行動特定装置
２０      タグ
３０      漏洩同軸アンテナ
１００    行動特定処理部
１０１    取得部
１０２    前処理部
１０３    指標値算出部
１０４    行動特定部
２００    測定データ記憶部
３００    行動特定モデル
【図１】
2020036627.tif 000003
【図２】
2020036627.tif 000004
【図３】
2020036627.tif 000005
【図４】
2020036627.tif 000006
【図５】
2020036627.tif 000007
【図６】
2020036627.tif 000008
【図７】
2020036627.tif 000009
【図８】
2020036627.tif 000010
【図９】
2020036627.tif 000011
【図１０】
2020036627.tif 000012
【図１１】
2020036627.tif 000013
【図１２】
2020036627.tif 000014