(57)【要約】 【課題】牛のＢＣＳ（ボディコンディションスコア）評価に客観性を持たせると同時に、その自動化により効率的な酪農経営の支援を図ることができる牛のＢＣＳの評価装置、評価方法及び評価プログラムを提供する。【解決手段】牛のＢＣＳの評価装置であって、評価装置を画像入力部１と、画像入力部１で入力された画像８のエッジ抽出部２と、エッジからＲＯＩ（処理対象領域）９を抽出し、ＲＯＩ９内の輪郭線１０を抽出する輪郭線抽出部３と、輪郭線１０を細線化及び滑線化する細線化部４と、細線化及び滑線化された輪郭線１１から近似曲線を求める近似曲線生成部５と、近似曲線１２と輪郭線１０の関係から牛のＢＣＳを評価するために必要な特徴パラメータを求める特徴パラメータ抽出部６と、特徴パラメータからＢＣＳを求めるＢＣＳ評価部７とから構成されることを特徴とする。【選択図】図１
2019187277.tif 000002
【特許請求の範囲】
【請求項１】
牛の臀部後方から撮影された写真あるいはビデオフレーム、又は、牛の臀部上方から撮影された距離画像の少なくともいずれか一方で構成された入力画像を取り込む画像入力部と、
該取り込まれた入力画像が前記牛の臀部後方から撮影された写真あるいはビデオフレームであることに応じてエッジ抽出及びノイズ除去を行うエッジ抽出部と、
前記エッジ抽出部で抽出されたエッジ、又は、前記牛の臀部上方から撮影された距離画像のいずれか一方からＲＯＩ（処理対象領域）を抽出し、該ＲＯＩ内の輪郭線を抽出する輪郭線抽出部と、
該輪郭線化部で抽出された輪郭線を細線化及び滑線化する細線化部と、
該細線化部で得られた輪郭線から近似曲線を求める近似曲線生成部と、
該近似曲線生成部で生成された近似曲線と前記輪郭線との関係から牛のＢＣＳ（ボディコンディションスコア）を評価するために必要な特徴パラメータを抽出する特徴パラメータ抽出部と、
該抽出された特徴パラメータから前記ＢＣＳを求めるＢＣＳ評価部と、
を備えることを特徴とする牛のボディコンディションスコアの評価装置。
【請求項２】
牛の臀部後方から撮影された写真あるいはビデオフレーム、又は、牛の臀部上方から撮影された距離画像の少なくともいずれか一方で構成された入力画像を取り込むステップと、
該取り込まれた入力画像が前記牛の臀部後方から撮影された写真あるいはビデオフレームであることに応じてエッジ抽出及びノイズ除去を行うステップと、
前記抽出したエッジ、又は、前記牛の臀部上方から撮影された距離画像の少なくとも一方からＲＯＩ（処理対象領域）を抽出し、該抽出した前記ＲＯＩ内の輪郭線を抽出するステップと、
該抽出された輪郭線を細線化及び滑線化するステップと、
該細線化及び滑線化された前記輪郭線から近似曲線を生成するステップと、
該近似曲線と前記輪郭線の関係から牛のＢＣＳを評価するために必要な特徴パラメータを抽出するステップと、
該抽出した特徴パラメータから前記ＢＣＳを求めると共に評価するステップと、
を備えたことを特徴とする牛のボディコンディションスコアの評価方法。
【請求項３】
前記輪郭線を抽出するステップにおいて、
垂直線から一定の角度傾いた直線と前記輪郭線との接点を用いて、前記ＲＯＩの開始点と終了点を決定するステップを備えている
ことを特徴とする請求項２に記載の牛のボディコンディションスコアの評価方法。
【請求項４】
牛の臀部後方から撮影された写真あるいはビデオフレーム、又は、牛の臀部上方から撮影された距離画像の少なくともいずれか一方で構成された入力画像を取り込むステップ処理と、
該取り込まれた入力画像が前記牛の臀部後方から撮影された写真あるいはビデオフレームであることに応じてエッジ抽出及びノイズ除去を行うステップ処理と、
該抽出したエッジ、又は、前記牛の臀部上方から撮影された距離画像の少なくとも一方からＲＯＩ（処理対象領域）を抽出するステップ処理と、
該抽出した前記ＲＯＩ内の輪郭線を抽出するステップ処理と、
該抽出した前記輪郭線から近似値を求め近似曲線を生成するステップ処理と、
該近似曲線と前記輪郭線との関係から牛のＢＣＳコンディションスコアを評価するために必要な特徴パラメータを求めるステップ処理と、
該求めた特徴パラメータから前記ＢＣＳを求めると共に評価するステップ処理と、
をコンピュータに実行させることを特徴とする牛のボディコンディションスコアの評価プログラム。
【請求項５】
前記ＲＯＩ内の輪郭線を抽出するステップ処理において、
垂直線から一定の角度傾いた直線と前記輪郭線との接点を用いて、前記ＲＯＩの開始点と終了点を決定するステップ処理を備えている
ことを特徴とする請求項４に記載の牛のボディコンディションスコアの評価プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、家畜、とくに牛の臀部後方の画像あるいは距離画像を用いて、牛のボディコンディションスコアを評価するための牛のボディコンディションの評価装置、評価方法及び評価プログラムに関する。
【背景技術】
【０００２】
家畜を飼育する場合には、肉や脂肪の付き、肥り具合等の生育状態を判断することが必要となる。例えば、肥痩状態評価方法として、ＢＣＳ（Ｂｏｄｙ  Ｃｏｎｄｉｔｉｏｎ  Ｓｃｏｒｅ：ボディコンディションスコア）を求めて肥痩状態を判断する方法が知られている。牛のＢＣＳは、体脂肪の蓄積状態を示す重要な指標である。その値は給餌の量や効率、健康状態などを示すため、定期的にチェックすることが重要である。しかしながら、従来は、酪農家が目視や触診によって、一頭ずつ目視や触診に基づいてＢＣＳを評価しているので、一頭あたりの評価時間がかかり、評価する頭数が多ければ膨大な時間と労力を要し、また、評価者による個人差も生じている。
【０００３】
例えば、乳用牛は泌乳期に多くの栄養を使って乳を生産している。とくに、泌乳最盛期では、餌から得られる栄養だけでは消費した分に追いつけないため、痩せてきて、乾乳期に太ってくるという傾向がある。このような身体の状態の変化が激しい牛は、健康状態を損ねやすい傾向にある。また、乾乳期に太りすぎてしまった牛は、分娩時に事故を起こす可能性が高くなる。空胎期、妊娠末期、授乳期などに応じて、最適なＢＣＳがあり、安定した状態に保つことが、健康で生産性の高い牛に育てる秘訣である。
【０００４】
そのため、一日当たりの乳量に応じた飼料給与が大切になる。肥り過ぎは飼料の給与過多が原因である場合が多く、逆に、やせすぎは飼料不足による場合が多い。いすれの場合も繁殖成績が悪くなるケースが多くみられる。一般に妊娠末期の牛は過肥（肥りすぎ）の傾向があるが、肥り過ぎは難産になりやすく、分娩後には食欲低下や代謝病のケトージス症に罹りやすくなる。また、授乳期はやせ過ぎになりやすく、そのため泌乳量が減少してしまう。
【０００５】
とくに多頭飼養の経営では、ＢＣＳが自動的に評価できれば、時系列情報を用いた個体ごとのＢＣＳ制御が容易になる。また、システムコストが下げられれば、小規模な農家でも導入が可能になり、省力化に繋がる。牛を健康な（病気になりにくい）状態で飼うためにＩＴ技術を活用した酪農は欠かせないものである。
【０００６】
以上のように、ＢＣＳは極めて重要であるにも関わらず、その評価には多くの時間と手間を要するため、従来多くの酪農家で定期的な評価は行われてこなかったのが現状である。そこで、ＢＣＳ評価に客観性を持たせると同時にその自動化により、効率的な酪農経営の支援を行うシステムの構築が切望されてきた。そこで、ＢＣＳについて、牛を三次元計測した三次元画像に基づいて算出する方法や管理システム及び管理プログラムが種々提案されている。
【０００７】
例えば、特許文献１には、簡便に牛体のさまざまな特性値を推定することができる牛体診断システムが開示されている。このシステムは、診断対象の牛を撮像して、その牛の臀部の体表面の３次元画像データを生成する３次元スキャン装置と、３次元画像データから、牛の臀部の体表面における任意の２点間の測地線長を測定する測地線長測定部と、体重を少なくとも一つの目的変数とし、左右腰角間の測地線長を少なくとも一つの説明変数とし、牛の個体情報を固定効果とする一般化線形回帰モデルに基づいて、測地線長測定部によって測定された測地線長のうち少なくとも左右腰角間の測地線長から牛体特性値として少なくとも診断対象の牛の体重を推定する牛体特性値推定部とを備え、牛のさまざまな特性値を推定することができ、酪農の大規模化に対応して大量の家畜牛の健康状態を効率よく診断することができるとされている（請求項１、段落「０００１」、「００１７」等参照。）。
【０００８】
特許文献２には、牛の撮像画像上に入力された線画の形状を判別し、判別された線画の形状に基づいて、撮像画像の牛のスコアを決定する処理をコンピュータで実行し、牛の肥痩状態評価のバラツキを抑制することができる管理方法、管理プログラム、管理装置および管理システムが開示されている（請求項１、段落「００１０」等参照。）。
【０００９】
特許文献３には、動物（牛）の身体のスコア（ｂｏｄｙ  ｃｏｎｄｉｔｉｏｎ  ｓｃｏｒｅ（ＢＳＣ））を決定するための装置であって、動物の方へ向けられかつ動物の少なくとも一つの三次元画像を記録するために与えられた三次元カメラシステム；及び三次元カメラシステムによって記録された三次元画像から動物の一部分の三次元表面表現を形成し、三次元表面表現を統計的に分析し、三次元表面表現の統計的に分析された表面に基づいて動物の身体の状態のスコアを決定するために与えられ、かつ三次元カメラシステムに接続された画像処理装置を含み、統計的分析は、三次元表面表現の表面のパラメータのヒストグラムの統計的特性の計算を含む装置及び方法が開示されている（請求項１、段落「０００１」等参照。）。
【００１０】
特許文献４には、動物（牛）の歩行数をカウントしてカウントデータを出力するカウントセンサと、一日を複数個に区分した各基準時刻から過去にさかのぼった所定の基準時刻間のカウントデータが積算された各積算カウントデータを出力する積算カウンタ手段と、各積算カウントデータを過去の所定日数間における各積算カウントデータの平均値又はその平均値と標準偏差に正の係数を乗じた数値との和で除した相対比を演算する相対比演算手段と、相対比が予め設定した設定値と比較して大きいとき警報信号を出力する出力手段とを備え、動物の歩行数の変化が高精度で検知されるので、放牧された家畜類等の動物が発情期であることを高精度で発見し管理して、家畜を合理的に生産できるとともに、家畜飼養管理の省力化を行うことができる運動量管理装置、及びその運動量管理装置を用いた運動量管理システム、並びに運動量管理プログラムが開示されている（請求項１、段落「０００１」、「００４０」等参照。）。
【００１１】
特許文献５には、放牧地に放牧されている家畜（牛）に装着された歩数計測手段により計測された歩数の測定結果を取得し、その時事系列変化に基づいて放牧地の植生状態が放牧に適した状態であるか否かを判定する植生判定プログラムを搭載した植生判定装置が開示されており、この植生判定装置が判定した結果を可搬型の通信端末に送信することが記載されている（請求項１、段落「００３０」等参照。）。
【先行技術文献】
【特許文献】
【００１２】
【特許文献１】特開２０１６−５９３００号公報
【特許文献２】特開２０１５−１７３７３２号公報
【特許文献３】特許第５５１９６９３号公報
【特許文献４】特開平１１−１２８２１０号公報
【特許文献５】特許第５９６２８５２号公報
【発明の概要】
【発明が解決しようとする課題】
【００１３】
本発明もまた、上記のような従来技術の課題に鑑み、牛のＢＣＳ評価に客観性を持たせると同時に、その自動化により効率的な酪農経営の支援を図ることができる牛のボディコンディションスコアの評価装置、評価方法及び評価プログラムを提供することを目的とするものである。
【課題を解決するための手段】
【００１４】
そこで、本発明の請求項１に係る牛のボディコンディションスコアの評価装置は、牛の臀部後方から撮影された写真あるいはビデオフレーム、又は、牛の臀部上方から撮影された距離画像の少なくともいずれか一方で構成された入力画像を取り込む画像入力部と、該取り込まれた入力画像が前記牛の臀部後方から撮影された写真あるいはビデオフレームであることに応じてエッジ抽出及びノイズ除去を行うエッジ抽出部と、前記エッジ抽出部で抽出されたエッジ、又は、前記牛の臀部上方から撮影された距離画像のいずれか一方からＲＯＩ（処理対象領域）を抽出し、該ＲＯＩ内の輪郭線を抽出する輪郭線抽出部と、該輪郭線化部で抽出された輪郭線を細線化及び滑線化する細線化部と、該細線化部で得られた輪郭線から近似曲線を求める近似曲線生成部と、該近似曲線生成部で生成された近似曲線と前記輪郭線との関係から牛のＢＣＳ（ボディコンディションスコア）を評価するために必要な特徴パラメータを抽出する特徴パラメータ抽出部と、該抽出された特徴パラメータから前記ＢＣＳを求めるＢＣＳ評価部と、を備えることを特徴とする、
【００１５】
また、本発明の請求項２に係る牛のボディコンディションスコアの評価方法は、牛の臀部後方から撮影された写真あるいはビデオフレーム、又は、牛の臀部上方から撮影された距離画像の少なくともいずれか一方で構成された入力画像を取り込むステップと、該取り込まれた入力画像が前記牛の臀部後方から撮影された写真あるいはビデオフレームであることに応じてエッジ抽出及びノイズ除去を行うステップと、前記抽出したエッジ、又は、前記牛の臀部上方から撮影された距離画像の少なくとも一方からＲＯＩ（処理対象領域）を抽出し、該抽出した前記ＲＯＩ内の輪郭線を抽出するステップと、該抽出された輪郭線を細線化及び滑線化するステップと、該細線化及び滑線化された前記輪郭線から近似曲線を生成するステップと、該近似曲線と前記輪郭線の関係から牛のＢＣＳを評価するために必要な特徴パラメータを抽出するステップと、該抽出した特徴パラメータから前記ＢＣＳを求めると共に評価するステップと、を備えたことを特徴とする。
【００１６】
また、本発明の請求項３に係る牛のボディコンディションスコアの評価方法は、請求項２に係る牛のボディコンディションスコアの評価方法であって、前記輪郭線を抽出するステップにおいて、垂直線から一定の角度傾いた直線と前記輪郭線との接点を用いて、前記ＲＯＩの開始点と終了点を決定するステップを備えていることを特徴とする。
【００１７】
そして、本発明の請求項４に係る牛のボディコンディションスコアの評価プログラムは、牛の臀部後方から撮影された写真あるいはビデオフレーム、又は、牛の臀部上方から撮影された距離画像の少なくともいずれか一方で構成された入力画像を取り込むステップ処理と、該取り込まれた入力画像が前記牛の臀部後方から撮影された写真あるいはビデオフレームであることに応じてエッジ抽出及びノイズ除去を行うステップ処理と、該抽出したエッジ、又は、前記牛の臀部上方から撮影された距離画像の少なくとも一方からＲＯＩ（処理対象領域）を抽出するステップ処理と、該抽出した前記ＲＯＩ内の輪郭線を抽出するステップ処理と、該抽出した前記輪郭線から近似値を求め近似曲線を生成するステップ処理と、該近似曲線と前記輪郭線との関係から牛のボディコンディションスコアを評価するために必要な特徴パラメータを求めるステップ処理と、該求めた特徴パラメータから前記ＢＣＳを求めると共に評価するステップ処理と、をコンピュータに実行させることを特徴とする。
【００１８】
さらに、本発明の請求項５に係る牛のボディコンディションスコアの評価プログラムは、請求項４に係る牛のボディコンディションスコアの評価プログラムであって、前記ＲＯＩ内の輪郭線を抽出するステップ処理において、垂直線から一定の角度傾いた直線と前記輪郭線との接点を用いて、前記ＲＯＩの開始点と終了点を決定するステップ処理を備えていることを特徴とする。
【発明の効果】
【００１９】
本発明によれば、以下の優れた効果がある、
（１）測定者の違いによる個人差がなくなり、ＢＣＳの客観的評価が可能となることで、ＢＣＳと事故の相関関係などのテータ共有が可能となる。
（２）汎用品だけでシステムを構成でき、小型化・コストダウンが可能であるので、酪農経営の効率化に繋がり、大規模農家だけでなく、小・中規模農家にもシステムを導入することができる。
（３）自動化が促進されることで、時系列で牛の身体変化の様子が記録・観察でき、より正確で詳細なデータに基づく体調管理が可能となる。
（４）３Ｄカメラの導入により、夜間や塵埃に対して頑健なシステム構成が可能となり、また、ビデオカメラだけによる構成も可能であるので、システム導入の際の選択肢が多く、適宜必要に応じた構成が可能である。
【図面の簡単な説明】
【００２０】
【図１】本発明を実施するためのシステム構成及び処理プロセスの一例を示すブロック図兼フローチャートである。
【図２】（ａ）は入力画像を示す図、（ｂ）は抽出したエッジ画像データを示す図、（ｃ）は輪郭線を示す図、（ｄ）は輪郭線の細線化データを示す図、（ｅ）は生成した近似曲線のデータを示す図である。
【図３】データ取得のための牛の撮像環境を模式的に示す概要図である。
【図４】（ａ）は牛の臀部からみた輪郭線とＢＣＳとの相関を示す図、（ｂ）はＢＣＳの違いによる牛の臀部を比較した図である。
【図５】（ａ）〜（ｄ）は、近似曲線の両端点を求める方法を示す説明図、（ｅ）は線形回帰分析を行うための輪郭線と近似曲線の概念図である。
【図６】ＢＣＳの指標としてヒストグラムを用いる方法を示す概念図である。
【図７】決定係数、多項式近似値、ＢＣＳの相関を示すグラフと図である。
【発明を実施するための形態】
【００２１】
以下、実施例を挙げて本発明を具体的に説明するが、本発明はこれらの実施例に限定されるものではない。
【００２２】
本発明では、ＢＣＳを評価する３つの手法「（１）線形回帰分布を用いて評価する、（２）ヒストグラムを用いてその特徴から評価する、（３）決定係数から評価する」と、輪郭線のＲＯＩ（解析対象領域）を決める手法「（１）一定の緩い傾斜をもった直線を用いて決定する」を提供する。
【００２３】
図４は、ＢＣＳの具体例であるが、この図から分かるように牛の臀部を後方から見た輪郭線とＢＣＳには強い相関があり、本方式の根拠を示している。これらは２次元画像を元にした説明図であるが、３次元カメラからの距離画像にも適用できる。牛の臀部を上から見た場合は、ＢＣＳの違いは３Ｄカメラの距離画像に現れる。
【００２４】
本発明では、通常のビデオカメラあるいは３次元カメラで撮像した画像データを元に画像処理技術を用いて処理し、以下に挙げる４手法でＢＣＳを評価する方法を提供する。
（１）近似曲線の両端点を求める方法
（２）線形回帰分析を用いる方法
（３）ヒストグラムを用いる方法
（４）決定係数を用いる方法
【００２５】
本発明で用いた３次元カメラのデータは牛の背骨と垂直にスライスした断面（ｚ＝一定）、ｘ−ｙ空間での輪郭線データとして用いる。一番簡単な方法は、それぞれのスライス画像に応じて誤差を計算し、さらにｚ軸方向に累積した結果に適当な閾値を設定して、ＢＣＳを計算する。また、ＢＣＳを従属変数とする回帰分析により求めることもできる。回帰分析に用いる線形回帰式、非線形回帰式は、周知であるので、説明を省略する。
【実施例】
【００２６】
本発明では、牛の後方臀部周辺の脂肪蓄積及び骨の凹凸の激しさから、ＢＣＳ（肥り過ぎ、痩せ過ぎを判断する一般的な指標）を評価する方法を提供する。入力装置としては、普通のビデオカメラを用いる場合と３次元カメラを用いる場合を考える。ここでは、画像処理技術によりＢＣＳを評価するための特徴パラメータを抽出し、簡易かつ非侵襲型のＢＣＳ評価システムを提供する。図１及び図２は、それぞれ処理プロセスの例、図３は牛の撮影環境の概念を示している。
【００２７】
本発明に係る牛のボディコンディションスコアの評価装置は、図３に示すように、ビデオカメラ１３を使用して牛１５の後方臀部周辺を後方から撮像したビデオ画像又は３次元カメラ１４を使用して牛１５の後方臀部周辺を上方から撮像した距離画像を入力画像として取り込む。
評価装置の構成は、図１及び図２に示すように、入力画像（「牛の臀部後方から撮影された写真あるいはビデオフレーム」又は「牛の臀部上方から撮影された距離画像」）８を取り込む画像入力部１と、画像入力部１で入力された画像８が「牛の臀部後方から撮影された写真あるいはビデオフレーム」であることに応じてエッジ抽出及びノイズ除去を行うエッジ抽出部２と、エッジ抽出部２で抽出されたエッジ、又は、「牛の臀部上方から撮影された距離画像」の少なくもいずれか一方からＲＯＩ（処理対象領域；Ｒｅｇｉｏｎ  ｏｆ  ｉｎｔｅｒｅｓｔ）９を抽出し、ＲＯＩ９内の輪郭線１０を抽出する輪郭線抽出部３と、輪郭線化部３で得られた輪郭線１０を細線化及び滑線化する細線化部４と、細線化部４で得られた細線化及び滑線化された輪郭線１１から近似曲線を求める近似曲線生成部５と、近似曲線生成部５で生成された近似曲線１２と輪郭線１０の関係から牛のボディコンディションスコアを評価するために必要な特徴パラメータを抽出する特徴パラメータ抽出部６と、抽出された特徴パラメータからＢＣＳ（ボディコンディションスコア）を求めるＢＣＳ評価部７とから構成される。なお、図２（ｃ）、図２（ｅ）において符号１０ａで示す曲線は、説明を分かりやすくするために、輪郭線１０の凹凸を上下に拡大して示したものである。また、符号１〜７で示される評価装置の各構成部分は、周知のコンピュータ、演算装置、記憶装置、表示装置等を含む電子装置、電子機器、電子回路等で構成することができる。
【００２８】
すなわち、本発明の牛のボディコンディションスコアの評価方法は、牛の臀部後方から撮影された写真あるいはビデオフレーム又は牛の臀部上方から撮影された距離画像のいずれかで構成された入力画像を取り込むステップＳ１と、該取り込まれた入力画像が「牛の臀部後方から撮影された写真あるいはビデオフレーム」であることに応じてエッジ抽出及びノイズ除去を行うステップＳ２と、抽出したエッジ、又は、「牛の臀部後方から撮影された写真あるいはビデオフレーム」の少なくともいずれか一方からＲＯＩ（処理対象領域）を抽出し、抽出した前記ＲＯＩ内の輪郭線を抽出するステップＳ３と、抽出された輪郭線を細線化及び滑線化するステップＳ４と、細線化及び滑線化された前記輪郭線から近似値を求め近似曲線を生成するステップＳ５と、該近似曲線と前記輪郭線の関係から牛のボディコンディションスコアを評価するために必要な特徴パラメータを抽出するステップＳ６と、抽出した特徴パラメータからＢＣＳを求めると共に評価するステップＳ７とからなり、これらステップ処理プログラムをコンピュータに実行させる。
【００２９】
「撮像環境」
３次元カメラを用いる場合は、牛１５の後方臀部全体の距離画像が一度の撮影で得られるので、一般的には複数台の３次元カメラを用いる必要はなく１台で十分である。図３に示すように、牛の後方臀部を真上から撮影した距離画像を用いる。距離画像を入力として用いる場合には、エッジ抽出を不要とし、輪郭線１０の抽出から開始することができる。つまり、図１において、エッジ抽出部２、および、エッジ抽出及びノイズを抽出するステップＳ２は不要であり、画像入力部１で取り込まれた距離画像を、「エッジ抽出及びノイズを抽出」せずに、直接、輪郭線抽出部３に入力して、ＲＯＩ９を抽出し、ＲＯＩ９内の輪郭線１０を抽出することができる。なお、画像入力部１の出力部の後段に画像データ判別手段（不図示）を設け、画像入力部１から出力される画像データがビデオカメラ１３で撮像された写真あるいはビデオフレームの場合は出力データをエッジ抽出部２に入力し、画像入力部１から出力される画像データが３次元カメラ１４で撮像された距離画像の場合は出力データを輪郭線抽出部３に入力する構成にしてもよい。また、画像データを判別して取り込む手段を、エッジ抽出部２の入力部、輪郭線抽出部３の入力部に、それぞれ設け、画像入力部１から出力される画像データがビデオカメラ１３で撮像された写真あるいはビデオフレームの場合は、エッジ抽出部２で取り込み、画像入力部１から出力される画像データが３次元カメラ１４で撮像された距離画像の場合は、輪郭線抽出部３で取り込むように構成してもよい。
【００３０】
図４からも分かるように、ＢＣＳが低い個体ほど骨格が浮かび上がり、体表面の凹凸の差が大きくなる。つまり、この凹凸具合を数値化することでＢＣＳ評価のための特徴量とすることができる。なお、距離画像とは、撮像対象物までの距離を明るさに対応して表したもので、近くにあるものほど明るく、遠くのものは暗くなる。
【００３１】
「近似曲線の開始点と終了点の決め方」
解析の対象となる領域ＲＯＩ９、すなわち、近似すべき輪郭線の両端点（開始点と終了点）を決めるために、垂直線と輪郭線の接点を用いると、安定性の低下や誤差の増大につながる場合がある。そこで、本発明では、垂直からθ°傾いた直線と輪郭線との接点を採用することで、安定性及び精度の向上を図る方法を提供する。
【００３２】
垂直線と輪郭線との接点を用いた例として、図５（ａ）及び図５（ｂ）に示すが、図５（ａ）は比較的安定しているものの、解析の対象となる臀部より下の部分まで拾ってしまう。図５（ｂ）は不安定で、臀部のみでなく脚部まで拾ってしまうので、ＢＣＳ評価のための特徴パラメータの抽出に悪影響を与える。そこで、垂直線でなく一定の傾き（θ）を持った直線を使うことで、安定した端点が得られる。θには厳密さは要求されず、５°〜４５°の範囲で、牛種や個体差に応じて決定する。図５（ｃ）及び図５（ｄ）の例ではθ＝１０°としている。
【００３３】
「線形回帰分析によるＢＣＳ評価」
牛の臀部の輪郭の凹凸具合を数値化する方法として牛領域のみの距離画像からデータを一行ごとに取り出し、グラフ上にプロットする。その後プロットされたデータに最小二乗法によって近似される放物線を描き（高次の多項式を使ってもよい）、この放物線を用いて、牛の体表面のデータのＲＭＳＤを計算することで、牛体表面の凹凸具合を数値化することができる。体表面との差を累積することでＢＣＳを評価する。本手法の概念図を図５（ｅ）に、ＢＣＳ評価結果を表１に示す。複数行の距離画像から得られる累積誤差に重みづけ加算した値を用いることもできる。ＲＭＳＤ（標準偏差）の式は下記に示すとおり、周知のものである。近似曲線は輪郭をスムーズに近似できるものであれば何でもよく、制御点から構成される自由曲線等を用いてもよい。経験的には、「線形回帰分析によるＢＣＳ評価」には２次多項式が単純かつ効果的であり、高次の多項式を用いると輪郭の形状変化に過敏に追従しすぎて、逆に精度の低下を招く原因となる。１次多項式も２次多項式に比べ、精度は低下する。
【００３４】
【表１】
2019187277.tif 000003
【００３５】
表１は最小二乗法による線形回帰分析を行うことで、ＢＣＳ推定を行った結果であり、本発明による提案手法の有効性を確認できた。同様の手法は普通の写真を用いた場合にも適用できる。
【００３６】
具体的には、各点における近似曲線と輪郭線との差の絶対値を累積した値、すなわち、図５（ｅ）の実線と破線の間の面積を｜ＡＣ−ＴＦ｜で表す。図５（ｅ）の輪郭線に相当する部分ＴＣは、３Ｄカメラを用いる場合は、各スライスから得られるので、｜ＡＣ−ＴＣ｜の値を重み付け累積して用いることもできる。数１は、線形回帰モデルによるＢＣＳ評価式である。
【数１】
2019187277.tif 000004
【００３７】
「ヒストグラムを用いる方法」
図６は、ＢＣＳ評価のための指標とした例を示す。ヒストグラムには明らかな差がみられる。また、相関係数の値にもＢＣＳの大小が反映されているので、その評価のために用いることができる。ヒストグラムによるグラフ表示方法自体は周知のものである。
【００３８】
「決定係数を用いる方法」
ＢＣＳが小さい数値からだんだん大きくなり、４〜５に近づくと、蓄積した脂肪により体表面の凹凸が明瞭に見えなくなる。図７（ｂ）の輪郭線に示されるように体表面が凸体となり、上述の凸包を利用する手法では、近似曲線と輪郭線がほぼ同一となり、その差がなくなる。牛種や個体差による違いはあるが、上述したヒストグラムを用いる方法では、ＢＣＳが大きくなった時に、精度が低下する傾向がある。
【００３９】
そのため、決定係数を用いる方法を提供する。ｎ次多項式を用いて近似式を求めるまでは線形回帰分析を用いた方法と同じである。近似式が求まったら決定係数を求める。決定係数を用いると、４以上のＢＣＳ（４〜５）の評価が可能となる。輪郭線の形状によってはさらに高次の多項式を用いると精度が向上する場合がある。
【００４０】
図７（ｃ）のＣｒｏｓｓｂｒｅｄ  ＣｏｗｓのＢＣＳ＝２，３の場合は視覚的にも似通っており、決定係数Ｒ２もよく似た値となっている。
数２は決定係数Ｒ２の導出式であり、周知である。図７（ｄ）及び（ｅ）のグラフから分かるように、４次元多項式近似を用いるとＣｒｏｓｓｂｒｅｄ  ＣｏｗｓとＪｅｒｓｅｙ  Ｃｏｗｓ２種類共にＢＣＳ＝４と５は識別可能であるが、ＢＣＳ＝２と３は難しい場合がある。図７（ｅ）より、６次多項式近似を用いると、ＢＣＳ＝１〜３は識別可能であるが、ＢＣＳ４と５は識別が難しい場合がある。逆に４次元多項式近似を用いると、ＢＣＳ＝３〜５とそれ以下は識別可能であるが、ＢＣＳ＝１と２は識別が難しい場合がある。なお、図７（ｅ）において、２次近似式で求めた近似曲線と、３次近似式で求めた近似曲線とが重なって表示されている。
【数２】
2019187277.tif 000005
【００４１】
「ＢＣＳ評価方法の選択と閾値」
以上、前述したどの方法が有効か、ＢＣＳ評価にための抽出特徴の閾値をいくらに設定すればよいのか、牛種や個体差を考慮して、個別に実験的に決めることで良好な結果が得られる。一律に決めることもできるが、より高い精度を得るためには、経時変化を観察しながら学習的に決定する。このように学習的に決定するための手段、方法、プログラムを、それぞれ、上述した牛のボディコンディションスコアの評価装置、牛のボディコンディションスコアの評価方法、牛のボディコンディションスコアの評価プログラムに設けることができる。なお、この学習手段、学習方法、学習プログラムは、上述した多項近似式の次数決定にも用いることができる。
【００４２】
導入時に多少の手間が増えるが、一個体については一旦決めれば大きな変動はないと考えられる。また、評価パラメータを使い分けることで、より精度の高いＢＣＳ評価が可能となる。実用的には、絶対値でなく各個体の相対的な経時変化が重要である。図７（ｃ）及び（ｃ）のグラフから分かるように、個体差より牛種の違いが大きい。
【産業上の利用可能性】
【００４３】
本発明は、畜産分野において牛の栄養状態評価、健康維持管理、餌料給与量の適正化に用いることができる。ＢＣＳ評価の自動化により高齢酪農家の負担軽減と効率化、牛のより詳細な健康管理、餌料給与の効率化に繋がる。画像処理による身体パラメータの抽出と健康管理、飼料給与の効率化に繋がる。
画像処理による身体パラメータの抽出と健康管理を結び付ける手法は、他の動物の健康維持にも利用できる。例えば、イルカ等の海洋大型生物の保全が問題視されているが、専門家によると肥満度などの客観的なデータの経年変化が得られれば、健康管理に加えてストレス等の精神状態の管理にも大きく貢献する。
【符号の説明】
【００４４】
１    画像入力部
２    エッジ抽出部
３    輪郭線抽出部
４    細線化部
５    近似曲線生成部
６    特徴パラメータ抽出部
７    ＢＣＳ評価部
８    入力画像
９    ＲＯＩ（処理対象領域）
１０  輪郭線
１１  細線化及び滑線化された輪郭線
１２  近似曲線
１３  ビデオカメラ
１４  ３次元カメラ
１５  牛
【図１】
2019187277.tif 000006
【図２】
2019187277.tif 000007
【図３】
2019187277.tif 000008
【図４】
2019187277.tif 000009
【図５】
2019187277.tif 000010
【図６】
2019187277.tif 000011
【図７】
2019187277.tif 000012