JP2023121504A - 牛の評価装置及び牛の評価方法 - Google Patents

Abstract

【課題】多様な飼育現場に適用することが可能な牛の評価装置及び牛の評価方法を提供すること。【解決手段】牛の評価装置１００は、立位にある牛を略真後ろから撮像した二次元の撮像画像５０を取得し、撮像画像５０から牛の輪郭線６０を抽出し、輪郭線６０の上部の左右両側に存在する２つの山部６１ａ、６１ｂのうちの一方の山部６１ａの輪郭線６０上の略中央の点である頂点６２ａから他方の山部６１ｂ側に伸びる第１曲線部６３ａ上の点４０３ａと第１曲線部６３ａとは反対側に伸びる第２曲線部６４ａ上の点４０４ａとを決定し、点４０３ａと点４０４ａとを結ぶ線分４１０ａの中点４０５ａを通り線分４１０ａに垂直な直線４２４ａが山部６１ａと交差する点４０６ａを特定し、点４０３ａと点４０６ａを結ぶ線分４１２ａと、点４０４ａと点４０６ａを結ぶ線分４１４ａと、がなす角度ω１を算出し、算出した角度ω１から牛の肉付きを評価する。【選択図】図６

Description

本発明は、牛の評価装置及び牛の評価方法に関する。

牛の飼養管理におけるエネルギーバランスを評価するために、牛の体脂肪の付着状態（肉付き状態）を評価することが行われている。体脂肪の付着状態を示す指標の１つとしてボディコンディションスコアがある。通常、ボディコンディションスコアの判定は獣医師や酪農の専門家によって行われるが、人が判定をするために判定基準があいまいとなり、判定にばらつきが生じてしまう。一方、動物の三次元画像に基づいて身体の状態を示すスコアを評価する方法（例えば特許文献１）や、動物を撮像した二次元画像に基づいて身体の状態を示すスコアを評価する方法（例えば特許文献２）が知られている。

特表２０１２－５１０２７８号公報 欧州特許第２０２７７７０号明細書

特許文献１に記載の評価方法は、動物の三次元画像を用いて身体の状態を示すスコアを算出するため、大掛かりな設備が設置された特定の場所での撮像が必要であり、多様な飼育現場に適用することは難しい。特許文献２に記載の評価方法では、動物の二次元画像を用いるが、動物を上方から撮像した二次元画像を用いるため、動物を上方から撮像できる設備が設置された特定の場所での撮像が必要であり、特許文献１と同様に、多様な飼育現場に適用することは難しい。

本発明は上記の課題に鑑みなされたものであり、多様な飼育現場に適用することが可能な牛の評価装置及び牛の評価方法を提供することを目的とする。

本発明の牛の評価装置は、立位にある牛を略真後ろ又は略正面から撮像した二次元の撮像画像を取得する画像取得部と、前記撮像画像から、前記牛の輪郭線を抽出する抽出部と、前記輪郭線の上部の左右両側に存在する２つの山部のうちの一方の山部の前記輪郭線上の略中央の点から他方の山部側に伸びる第１曲線部上の第１の点と前記第１曲線部とは反対側に伸びる第２曲線部上の第２の点とを決定する決定部と、前記第１の点と前記第２の点とを結ぶ第１線分の長さ、前記第１線分の略中点を通り前記第１線分に略垂直な直線が前記一方の山部と交差する第３の点と前記第１線分の略中点とを結ぶ第２線分の長さ、前記第２線分の長さに対する前記第１線分の長さの比、及び前記第１の点と前記第３の点を結ぶ第３線分と前記第２の点と前記第３の点を結ぶ第４線分とがなす角度、の少なくとも１つを算出する算出部と、前記算出部が算出した値から、前記牛の肉付きを評価する評価部と、を備える。

本発明によれば、多様な飼育現場に適用することが可能な牛の評価装置及び牛の評価方法を得ることができる。

図１は、第１の実施形態に係る牛の評価装置の構成を示すブロック図である。 図２は、第１の実施形態における制御部の機能構成を示すブロック図である。 図３は、画像取得部が取得する画像の一例を示す図である。 図４は、抽出部が実行する処理を説明するための図である。 図５は、決定部が実行する処理を説明するための図である。 図６（ａ）及び図６（ｂ）は、第１の実施形態における算出部が実行する処理を説明するための図である。 図７（ａ）から図７（ｃ）は、肉付きの異なる牛の撮像画像である。 図８（ａ）から図８（ｃ）は、牛の肉付きと、算出部が算出する角度ω１と、の関係を示す図である。 図９は、ボディコンディションスコアと、算出部が算出した角度ω１、角度ω２、及び角度ω１と角度ω２の平均値と、の相関を調査した実験結果である。 図１０は、実施形態に係る牛の評価装置のハードウェア構成を示すブロック図である。 図１１は、第１の実施形態における牛の肉付きの評価方法の一例を示すフローチャートである。 図１２（ａ）から図１２（ｃ）は、決定部による他の決定方法を示す図である。 図１３（ａ）及び図１３（ｂ）は、第２の実施形態における算出部が実行する処理を説明するための図である。 図１４（ａ）から図１４（ｃ）は、牛の肉付きと、算出部が算出する長さＸ１ａ、Ｘ２ａと、の関係を示す図である。 図１５は、ボディコンディションスコアと、算出部が算出した長さＸ１ａ、長さＸ１ｂ、及び長さＸ１ａと長さＸ１ｂの平均値と、の相関を調査した実験結果である。 図１６は、ボディコンディションスコアと、算出部が算出した長さＸ２ａ、長さＸ２ｂ、及び長さＸ２ａと長さＸ２ｂの平均値と、の相関を調査した実験結果である。 図１７は、ボディコンディションスコアと、算出部が算出した長さＸ２ａに対する長さＸ１の長さ比（Ｘ１ａ／Ｘ２ａ）、長さＸ２ｂに対する長さＸ１ｂの長さ比（Ｘ１ｂ／Ｘ２ｂ）、及び長さ比（Ｘ１ａ／Ｘ２ａ）と長さ比（Ｘ１ｂ／Ｘ２ｂ）の平均値と、の相関を調査した実験結果である。 図１８（ａ）は、立位にある牛を略真後ろから撮像した撮像画像、図１８（ｂ）は、図１８（ａ）の牛を略正面から撮像した撮像画像である。 図１９（ａ）及び図１９（ｂ）は、図１８（ａ）及び図１８（ｂ）の撮像画像から牛の輪郭線を抽出した抽出画像である。

《第１の実施形態》
以下、第１の実施形態に係る牛の評価装置１００について、図１～図１２に基づいて説明する。図１には、第１の実施形態に係る牛の評価装置１００の構成がブロック図にて示されている。図２には、制御部２２の機能構成がブロック図にて示されている。本第１の実施形態では、評価装置１００は、作業者が携帯して利用する携帯型情報機器であるとする。携帯型情報機器として、例えばスマートフォン、タブレット型パソコン、ＰＤＡ（Personal Digital Assistant）、スマートグラス等が挙げられる。評価装置１００は、例えば、電話機能やメール機能、インターネット等に接続するための通信機能、及びプログラムを実行するためのデータ処理機能等を有する。

評価装置１００は、図１に示すように、表示部１２と、操作部１４と、通信部１６と、撮像部１８と、記憶部２０と、制御部２２と、を備える。

表示部１２は、画像や、各種情報、及びタッチ操作ボタン等の操作入力用画像等を表示する表示デバイスである。

操作部１４は、タッチパネルやスイッチ等のデバイスである。タッチパネルは、作業者が触れたことに応じて情報入力を受け付け、受け付けた情報を制御部２２に送信する。タッチパネルは、例えば表示部１２に組み込まれている。したがって、タッチパネルは、作業者が表示部１２の表面をタッチすることに応じて、種々の情報入力を受け付ける。スイッチは、作業者から評価装置１００に対する操作を受け付ける操作部材であり、受け付けた情報を制御部２２に送信する。

通信部１６は、他の機器と近距離無線通信（例えばＮＦＣ（Near Field Communication））を行ったり、ネットワークに接続された他の機器と無線通信（例えば携帯電話回線や無線ＬＡＮ（Local Area Network）等を用いた通信）を行ったりする通信インターフェースである。

撮像部１８は、静止画や動画等の画像を撮像するデバイスである。撮像部１８は、評価装置１００を携帯する作業者が操作部１４を操作することに応じて撮像する。作業者は、牛の肉付きを評価する場合に、撮像部１８によって立位にある牛を略真後ろから撮像する。略真後ろから撮像するとは、牛の左右の腰角周りの輪郭が写る範囲で撮像することである。したがって、作業者は、例えば、牛の後方数ｍ（例えば１ｍ～２ｍ）の位置から、牛の左右の腰角周りの輪郭が写るように撮像する。

記憶部２０は、牛の撮像画像や牛の輪郭線を抽出した抽出画像、牛の肉付きを評価した評価結果等、各種情報を記憶する。

制御部２２は、牛の評価装置１００の各部を制御する処理部である。図２に示すように、制御部２２は、画像取得部３０、抽出部３２、決定部３８、算出部３４、及び評価部３６を有する。

図３には、画像取得部３０が取得する画像の一例が示されている。図３に示すように、画像取得部３０は、例えば、立位にある牛を略真後ろから撮像した撮像画像５０を取得する。取得する牛の撮像画像５０は、撮像部１８が撮像した画像であってもよいし、他の機器が撮像した画像を通信部１６により受信した画像であってもよい。また、取得する牛の撮像画像５０は、記憶部２０に記憶されていたものの中から作業者が操作部１４を用いて選択したものであってもよい。画像取得部３０は、取得した牛の撮像画像５０を表示部１２に表示してもよい。これにより、作業者は、これから肉付きを評価する牛を表示部１２上で視覚的に認識することができる。

図４は、抽出部３２が実行する処理を説明するための図である。図４に示すように、抽出部３２は、画像取得部３０が取得した牛の撮像画像５０から、牛の輪郭をなぞった輪郭線６０を抽出する。輪郭線６０の抽出は、一般的に知られた画像処理技術を用いることができる。抽出部３２は、抽出した牛の輪郭線６０のみを表した抽出画像５１を表示部１２に表示してもよい。これにより、作業者は、これから肉付きを評価する牛の輪郭を表示部１２上で視覚的に認識することができる。抽出画像５１の外形は撮像画像５０の外形と同じであり、画像の枠に対する牛の位置は抽出画像５１と撮像画像５０とで同じである。

輪郭線６０の上部の左右両側には山部６１ａ、６１ｂが存在する。山部６１ａ、６１ｂは、牛の腰角に対応する部分である。山部６１ａは、山部６１ａの輪郭線６０上の略中央の点を頂点６２ａとし、頂点６２ａから輪郭線６０の上部中央側（牛の背骨側）に伸びる第１曲線部６３ａと、頂点６２ａから輪郭線６０の側部下側（牛の腹側）に伸びる第２曲線部６４ａと、を有する。同様に、山部６１ｂは、山部６１ｂの輪郭線６０上の略中央の点を頂点６２ｂとし、頂点６２ｂから輪郭線６０の上部中央側（牛の背骨側）に伸びる第１曲線部６３ｂと、頂点６２ｂから輪郭線６０の側部下側（牛の腹側）に伸びる第２曲線部６４ｂと、を有する。なお、第１曲線部６３ａ、６３ｂは、山部６１ａ、６１ｂの輪郭線６０上の略中央の点（頂点６２ａ、６２ｂ）から他方の山部側（牛の背骨側）に伸びる曲線部であるといえる。また、第２曲線部６４ａ、６４ｂは、山部６１ａ、６１ｂの輪郭線６０上の略中央の点（頂点６２ａ、６２ｂ）から、第１曲線部６３ａ、６３ｂとは反対側に伸びる曲線部であるといえる。

山部６１ａの頂点６２ａと山部６１ｂの頂点６２ｂは、以下の方法により特定することができる。例えば、抽出画像５１内で上から左下に向かって傾いた傾斜直線５５を抽出画像５１の左上端点５２から輪郭線６０に向かって近づけて行き、傾斜直線５５が輪郭線６０に最初に接した点を左側に存在する山部６１ａの頂点６２ａと特定することができる。同様に、抽出画像５１内で上から右下に向かって傾いた傾斜直線５６を抽出画像５１の右上端点５４から輪郭線６０に向かって近づけて行き、傾斜直線５６が輪郭線６０に最初に接した点を右側に存在する山部６１ｂの頂点６２ｂと特定することができる。傾斜直線５５、５６は、抽出画像５１内の水平線から例えば３０°～６０°傾いた直線とすることが好ましく、４０°～５０°傾いた直線とすることがより好ましく、４５°傾いた直線とすることが更に好ましい。

山部６１ａの頂点６２ａと山部６１ｂの頂点６２ｂは、上記の方法によって特定される場合に限られず、その他の方法によって特定されてもよい。例えば、抽出画像５１内で上から左下に向かって第１の所定角度で傾いた傾斜直線５５ａを左上端点５２から輪郭線６０に向かって近づけて最初に接する第１の点を特定し、第２の所定角度で傾いた傾斜直線５５ｂを左上端点５２から輪郭線６０に向かって近づけて最初に接する第２の点を特定し、輪郭線６０上において第１の点と第２の点の略中央に位置する点を山部６１ａの頂点６２ａと特定してもよい。同様に、抽出画像５１内で上から右下に向かって第３の所定角度で傾いた傾斜直線５６ａを右上端点５４から輪郭線６０に向かって近づけて最初に接する第３の点を特定し、第４の所定角度で傾いた傾斜直線５６ｂを右上端点５４から輪郭線６０に向かって近づけて最初に接する第４の点を特定し、輪郭線６０上において第３の点と第４の点の略中央に位置する点を山部６１ｂの頂点６２ｂと特定してもよい。第１の所定角度および第３の所定角度は抽出画像５１内の水平線から１５°～３０°とすることが好ましく、第２の所定角度および第４の所定角度は抽出画像５１内の水平線から６０°～７５°とすることが好ましい。

また、例えば、抽出画像５１をＸ－Ｙ座標系に変換し、抽出画像５１の各画素のＸ、Ｙ座標から２つの山部６１ａ、６１ｂの頂点６２ａ、６２ｂを特定してもよい。例えば、輪郭線６０のなかで抽出画像５１の左上端点５２からの距離が最も短い点を山部６１ａの頂点６２ａとし、輪郭線６０のなかで抽出画像５１の右上端点５４からの距離が最も短い点を山部６１ｂの頂点６２ｂとすることができる。

また、上記いずれかの方法によって山部６１ａまたは山部６１ｂのいずれか一方の頂点を特定し、特定した頂点を通る抽出画像５１内の水平線が他方の山部の輪郭線６０に交わる点を他方の山部の頂点と特定してもよい。なお、上記の方法以外の方法によって山部６１ａの頂点６２ａおよび山部６１ｂの頂点６２ｂを特定してもよい。なお、山部６１ａ、６１ｂの範囲については、様々な方法で特定することができる。例えば、山部６１ａの範囲は、山部６１ａの頂点６２ａを中心とした所定半径の円に含まれる輪郭線６０の範囲とすることができる。所定半径は、例えば、山部６１ａの頂点６２ａと輪郭線６０のうち最も上側に位置する最上点９１との間の距離に所定の値（例えば０．５など）を掛けた値とすることができるが、その他の方法で得られる値であってもよい。また、山部６１ａの範囲は、山部６１ａの頂点６２ａから両側に所定長さの範囲としてもよい。所定長さは、例えば、上記所定半径の長さとしてもよいし、輪郭線６０全体の長さに所定の値を掛けた長さとしてもよい。山部６１ｂについても同様である。また、山部６１ａ、６１ｂの範囲については、機械学習等により特定することとしてもよい。

図５は、決定部３８が実行する処理を説明するための図である。図５は、輪郭線６０の上部を拡大した図であり、図の明瞭化のために模式的に図示している（以下の同様な図においても同じ）。図５に示すように、決定部３８は、輪郭線６０の上部の左右両側に存在する山部６１ａ、６１ｂのうちの一方の山部６１ａの第１曲線部６３ａ上の点４０３ａと第２曲線部６４ａ上の点４０４ａを決定する。本第１の実施形態では、決定部３８は、抽出画像５１（撮像画像５０）内で左右方向（水平方向）に伸びる水平線４２０を山部６１ａの上方から山部６１ａに向かって近づけて最初に山部６１ａに接する点を第１曲線部６３ａ上の点４０３ａと決定する。また、決定部３８は、抽出画像５１（撮像画像５０）内で上下方向（垂直方向）に伸びる垂直線４２２ａを山部６１ａの左側方から山部６１ａに向かって近づけて最初に山部６１ａに接する点を第２曲線部６４ａ上の点４０４ａと決定する。このように、本第１の実施形態では、決定部３８は、第１曲線部６３ａのうち最も上側に張り出した点を点４０３ａと決定し、第２曲線部６４ａのうち最も横側に張り出した点を点４０４ａと決定する。

決定部３８は、山部６１ａの第１曲線部６３ａ上の点４０３ａと第２曲線部６４ａ上の点４０４ａの決定の代わりに又は決定に加えて、山部６１ｂの第１曲線部６３ｂ上の点４０３ｂと第２曲線部６４ｂ上の点４０４ｂを決定してもよい。決定部３８は、水平線４２０を山部６１ｂの上方から山部６１ｂに向かって近づけて最初に山部６１ｂに接する点を第１曲線部６３ｂ上の点４０３ｂと決定し、抽出画像５１内で上下方向（垂直方向）に伸びる垂直線４２２ｂを山部６１ｂの右側方から山部６１ｂに向かって近づけて最初に山部６１ｂに接する点を第２曲線部６４ｂ上の点４０４ｂと決定することができる。

図６（ａ）及び図６（ｂ）は、第１の実施形態における算出部３４が実行する処理を説明するための図である。図６（ａ）に示すように、算出部３４は、点４０３ａと点４０４ａを結ぶ線分４１０ａの中点４０５ａを通り線分４１０ａに垂直な直線４２４ａが山部６１ａに交差する点４０６ａを特定する。そして、算出部３４は、点４０３ａと点４０６ａを結ぶ線分４１２ａと、点４０４ａと点４０６ａを結ぶ線分４１４ａと、がなす角度ω１を算出する。算出部３４が算出する角度ω１は、山部６１ａの傾斜の度合いを示す指標値の１つである。

評価部３６は、算出部３４が算出した角度ω１から、撮像画像５０に写された牛の肉付きを評価する。肉付きとして、例えばボディコンディションスコアを評価する。

また、図６（ｂ）に示すように、決定部３８が山部６１ｂの第１曲線部６３ｂ上の点４０３ｂと第２曲線部６４ｂ上の点４０４ｂを決定した場合では、算出部３４は、点４０３ｂと点４０４ｂを結ぶ線分４１０ｂの中点４０５ｂを通り線分４１０ｂに垂直な直線４２４ｂが山部６１ｂに交差する点４０６ｂを特定する。そして、算出部３４は、点４０３ｂと点４０６ｂを結ぶ線分４１２ｂと、点４０４ｂと点４０６ｂを結ぶ線分４１４ｂと、がなす角度ω２を算出する。算出部３４が算出する角度ω２は、山部６１ｂの傾斜の度合いを示す指標値の１つである。この場合、評価部３６は、算出部３４が算出した角度ω２から、撮像画像５０に写された牛の肉付きを評価する。

また、決定部３８が山部６１ａの第１曲線部６３ａ上の点４０３ａと第２曲線部６４ａ上の点４０４ａに加えて、山部６１ｂの第１曲線部６３ｂ上の点４０３ｂと第２曲線部６４ｂ上の点４０４ｂを決定した場合では、算出部３４は、角度ω１と角度ω２の両方を算出して角度ω１と角度ω２の平均値を算出してもよい。この場合、評価部３６は、算出部３４が算出した角度ω１と角度ω２の平均値から、撮像画像５０に写された牛の肉付きを評価してもよい。

ここで、算出部３４が算出する角度ω１、ω２と、ボディコンディションスコアと、の関係について説明する。なお、ボディコンディションスコアは、体脂肪の付着状態を示す指標であり、例えば１．００～５．００の点数を付ける場合があり、値が大きいほど肉付きが良いことを意味する。

図７（ａ）から図７（ｃ）は、肉付きの異なる牛の撮像画像５０である。図７（ａ）は、肉付きが悪い牛の撮像画像５０、図７（ｂ）は、図７（ａ）の牛よりも肉付きが良い牛の撮像画像５０、図７（ｃ）は、肉付きが更に良い牛の撮像画像５０である。図７（ａ）の牛は、専門家によるボディコンディションスコアの評価は２．２５であった。図７（ｂ）の牛は３．００であり、図７（ｃ）の牛は４．００であった。図７（ａ）から図７（ｃ）に示すように、肉付きが異なる（すなわちボディコンディションスコアが異なる）ことで、腰角周りの肉付きが変化し、立位にある牛を略真後ろから見たときの腰角の形状が異なることが確認できる。

図８（ａ）から図８（ｃ）は、牛の肉付きと、算出部３４が算出する角度ω１と、の関係を示す図である。図８（ａ）、図８（ｂ）、図８（ｃ）の順に肉付きの良い牛の輪郭線６０を示している。図８（ａ）から図８（ｃ）のように、肉付きが良いほど、腰角に対応する部分である山部６１ａが丸みを帯びた形状となる。このため、角度ω１は、肉付きが良いほど小さな値となることが分かった。なお、図示は省略するが、牛の肉付きと、算出部３４が算出する角度ω２と、の関係についても同じである。これらのことから、角度ω１、角度ω２、又は角度ω１と角度ω２の平均値を求め、これらを用いて牛の肉付きを評価できると考えられる。

［実験］
獣医師等の専門家が評価したボディコンディションスコアと、算出部３４が算出した角度ω１、角度ω２、及び角度ω１と角度ω２の平均値と、の相関を調査する実験を行った。実験は、乳牛（ホルスタイン種）に対して行った。

図９には、ボディコンディションスコアと、算出部３４が算出した角度ω１、角度ω２、及び角度ω１と角度ω２の平均値と、の相関を調査した実験結果が示されている。図９の横軸（ｘ軸）は専門家が評価したボディコンディションスコアであり、縦軸は角度である。図９中の複数の黒丸は角度ω１を用いたときの乳牛１頭１頭の結果を示していて、点線の直線はそのときの近似直線（回帰直線）を示している。複数の黒四角は角度ω２を用いたときの乳牛１頭１頭の結果を示していて、一点鎖線の直線はそのときの近似直線（回帰直線）を示している。複数の黒三角は角度ω１と角度ω２の平均値を用いたときの乳牛１頭１頭の結果を示していて、破線の直線はそのときの近似直線（回帰直線）を示している。図９に示すように、角度ω１、角度ω２、及び角度ω１と角度ω２の平均値のいずれを用いた場合であっても、角度が小さくなるほどボディコンディションスコアは大きくなっている。角度ω１を用いた点線の近似直線はｙ＝－７．４７２４ｘ＋１６２．８５と求まり、決定係数Ｒ^２は０．６４７８であった。角度ω２を用いた一点鎖線の近似直線はｙ＝－１０．５９４ｘ＋１７１．７６と求まり、決定係数Ｒ^２は０．７０４８であった。角度ω１と角度ω２の平均値を用いた破線の近似直線はｙ＝－９．０３３１ｘ＋１６７．３１と求まり、決定係数Ｒ^２は０．７３９７であった。このことから、ボディコンディションスコアと、算出部３４が算出する角度ω１、角度ω２、及び角度ω１と角度ω２の平均値と、には強い相関関係があることが確認された。

このように、ボディコンディションスコアと、算出部３４が算出する角度ω１、角度ω２、及び角度ω１と角度ω２の平均値と、には強い相関関係があることから、評価部３６が角度ω１、角度ω２、又は角度ω１と角度ω２の平均値から牛の肉付きを評価することで、牛の肉付きを精度良く評価できることが分かる。評価部３６は、例えば、記憶部２０に記憶された、角度とボディコンディションスコアとが関連付けられたスコア情報を用いてボディコンディションスコアを評価する。

表１には、記憶部２０に記憶されたスコア情報の一例が示されている。表１に示すように、スコア情報においては、角度とボディコンディションスコアとが関連付けられている。例えば、角度がａ°未満の場合のボディコンディションスコアは５．００、ａ°以上ｂ°未満の場合のボディコンディションスコアは４．５０等のように、角度とボディコンディションスコアとが関連付けられている。なお、表１では、ボディコンディションスコアが０．５刻みの場合を例に示したが、その他の場合でもよい。例えば、０．２５刻みの場合でもよいし、２．００～４．００の間は０．２５刻みで、それ以外は０．５０刻みの場合でもよい。表１のようなスコア情報として、角度ω１、角度ω２、及び角度ω１と角度ω２の平均値とで別々のスコア情報を作成してもよいし、角度ω１、角度ω２、又は角度ω１と角度ω２の平均値のいずれかスコア情報を作成して共用してもよい。スコア情報を共用する場合、角度ω１と角度ω２の平均値で作成したスコア情報を共用することで、角度ω１と角度ω２からボディコンディションスコアを評価する場合でも精度良く評価できる。

なお、記憶部２０に記憶されるスコア情報は、角度とボディコンディションスコアとが関連付けられた情報であれば、表１に示したスコア情報以外であってもよい。例えば、図９に示す近似直線のような一次関数（ｙ＝ａｘ＋ｂ）がスコア情報として記憶部２０に記憶されていてもよい。この場合、ｙに角度を代入することで得られたｘの値に対して四捨五入や、切り上げ、切り捨て等の端数処理を行ってボディコンディションスコアを評価してもよい。

図１０は、第１の実施形態に係る牛の評価装置１００のハードウェア構成を示すブロック図である。図１０に示すように、牛の評価装置１００は、記憶装置７０、メモリ７１、プロセッサ７２、通信インターフェース７３、入力装置７４、表示装置７５、及び撮像装置７６を有する。これらの各部はバス７７により相互に接続される。

記憶装置７０は、フラッシュメモリ等の不揮発性の半導体メモリであって、牛の評価プログラム７８を記憶する。

メモリ７１は、ＤＲＡＭ（Dynamic Random Access Memory）等のようなデータを一時的に記憶するハードウェアである。図１における記憶部２０は、記憶装置７０とメモリ７１によって実現される。

プロセッサ７２は、牛の評価装置１００の各部を制御するＣＰＵ（Central Processing Unit）等のハードウェアである。プロセッサ７２は、メモリ７１と協働して牛の評価プログラム７８を実行する。このように、メモリ７１とプロセッサ７２とが協働して牛の評価プログラム７８を実行することにより、図２における制御部２２の各機能が実現される。

通信インターフェース７３は、牛の評価装置１００を他の機器と近距離無線通信又はネットワークを介した無線通信をさせるためのハードウェアである。通信インターフェース７３により図１における通信部１６が実現される。

入力装置７４は、作業者が牛の評価装置１００を操作するためのタッチパネルやスイッチ等の入力デバイスである。入力装置７４により図１における操作部１４が実現される。

表示装置７５は、図１における表示部１２を実現するためのハードウェアであり、牛の撮像画像５０や抽出画像５１、牛の肉付きの評価結果等を表示するための液晶ディスプレイ等の表示デバイスである。

撮像装置７６は、図１における撮像部１８を実現するためのハードウェアであり、レンズや撮像素子等を備える。

次に、第１の実施形態における牛の肉付きの評価方法の一例について、図１１のフローチャートに沿って説明する。図１１の処理は、評価装置１００の制御部２２によって実行される。なお、記憶部２０には、一例として表１のようなスコア情報が格納されているとする。

図１１の処理では、まず、ステップＳ１０において、制御部２２の画像取得部３０は、牛の評価プログラム７８が起動されるまで待機する。この場合、作業者の操作（例えば操作部１４への入力操作又はマイク（不図示）への音声入力等）に応じて制御部２２が牛の評価プログラム７８を起動した段階で、ステップＳ１２に移行する。

ステップＳ１２に移行すると、画像取得部３０は、肉付きの評価を行う牛が撮像された撮像画像５０を取得する（図３参照）。牛の撮像画像５０は、立位にある牛が略真後ろから撮像された画像である。画像取得部３０は、例えば、作業者が撮像部１８を用いて撮像した牛の撮像画像を取得する。なお、画像取得部３０は、通信部１６を介して他の機器から受信した牛の撮像画像、又は記憶部２０に記憶された複数の牛の撮像画像の中から作業者が操作部１４を操作することで選択した牛の撮像画像を取得してもよい。なお、画像取得部３０は、牛の評価プログラム７８が起動した後に、作業者が撮像部１８によって牛を撮像できるよう、カメラ機能を立ち上げてもよい。作業者が撮像部１８によって牛を撮像した場合、撮像した牛の画像と、この撮像画像を選択するか否かを問う選択ボタンと、を表示部１２に表示し、作業者が選択しないことを選んだ場合、カメラ機能を再度立ち上げるようにしてもよい。

制御部２２の抽出部３２は、ステップＳ１４において、ステップＳ１２で取得した牛の撮像画像５０から牛の輪郭をなぞった輪郭線６０を抽出する（図４参照）。これにより、牛の輪郭線６０が示された抽出画像５１が生成される。輪郭線６０を抽出するにあたり、抽出部３２は、輪郭線６０を良好な精度で抽出するために、撮像画像５０内の牛以外の背景を除去する処理を行ってもよい。

次いで、ステップＳ１６では、制御部２２の決定部３８は、本第１の実施形態では、山部６１ａの第１曲線部６３ａ上の点４０３ａと第２曲線部６４ａ上の点４０４ａを決定する。決定部３８は、水平線４２０を山部６１ａの上方から山部６１ａに向かって近づけて最初に山部６１ａに接する点を第１曲線部６３ａ上の点４０３ａと決定し、垂直線４２２ａを山部６１ａの左側方から山部６１ａに向かって近づけて最初に山部６１ａに接する点を第２曲線部６４ａ上の点４０４ａと決定する（図５参照）。

次いで、ステップＳ１８では、制御部２２の算出部３４は、まず、点４０３ａと点４０４ａを結ぶ線分４１０ａの中点４０５ａを通り線分４１０ａに垂直な直線４２４ａが山部６１ａに交差する点４０６ａを特定する。そして、算出部３４は、点４０３ａと点４０６ａを結ぶ線分４１２ａと、点４０４ａと点４０６ａを結ぶ線分４１４ａと、がなす角度ω１を算出する（図６（ａ）参照）。

次いで、ステップＳ２０では、制御部２２の評価部３６は、ステップＳ１８で算出した角度ω１と、記憶部２０に記憶された表１のようなスコア情報とから、撮像画像５０に写された牛の肉付き（ボディコンディションスコア）を評価する。

次いで、ステップＳ２２では、評価部３６は、ステップＳ２０で評価したボディコンディションスコアを表示部１２に表示すると共に、記憶部２０に記憶する。例えば、評価部３６は、牛の撮像画像５０と共にボディコンディションスコアを表示部１２に表示してもよい。また、評価部３６は、作業者が操作部１４を操作して入力した牛の個体識別番号等の個体識別情報と、評価したボディコンディションスコアと、を関連付けて記憶部２０に記憶してもよい。なお、評価部３６は、作業者が操作部１４を操作することで入力した個体識別情報を取得する場合に限られず、牛の体に取り付けられた識別票の画像を撮像部１８が撮像した撮像画像から個体識別情報を取得する場合や、ＲＦＩＤ（Radio Frequency Identification）技術を用いて、牛の体に取り付けられたＲＦタグから通信部１６が個体識別情報を受信することで取得する場合等でもよい。また、カメラや通信装置を牛舎内外の通路に設置し、通路を通過する牛をセンサで検出できるようにしておき、通路を通過する牛をカメラが撮像した撮像画像を通信部１６が取得し、その撮像画像に写された識別票の画像から個体識別情報を取得してもよいし、通路を通過する牛に取り付けられたＲＦタグから通信装置が個体識別情報を受信し、通信部１６が通信装置から個体識別情報を受け取ることで取得してもよい。また、評価部３６は、ボディコンディションスコアを評価した牛の個体識別番号と、当該牛を撮像した撮像画像５０と、撮像日と、当該牛のボディコンディションスコアと、を関連付けた１つの情報として記憶部２０に記憶してもよい。

以上説明したように、本第１の実施形態によれば、画像取得部３０は、立位にある牛を略真後ろから撮像した二次元の撮像画像５０を取得する（図３参照）。抽出部３２は、撮像画像５０から、牛の輪郭線６０を抽出する（図４参照）。決定部３８は、輪郭線６０の上部の左右両側に存在する２つの山部６１ａ、６１ｂのうちの山部６１ａの第１曲線部６３ａ上の点４０３ａ（第１の点）と、第２曲線部６４ａ上の点４０４ａ（第２の点）と、を決定する（図５参照）。算出部３４は、点４０３ａと点４０４ａとを結ぶ線分４１０ａ（第１線分）の中点４０５ａを通り線分４１０ａに垂直な直線４２４ａが山部６１ａに交差する点４０６ａ（第３の点）を特定する。そして、算出部３４は、点４０３ａと点４０６ａを結ぶ線分４１２ａ（第３線分）と、点４０４ａと点４０６ａを結ぶ線分４１４ａ（第４線分）と、がなす角度ω１を算出する（図６（ａ）参照）。評価部３６は、算出部３４が算出した角度ω１から、撮像画像５０に写された牛の肉付きを評価する。このように、本第１の実施形態では、立位にある牛を略後ろから撮像した二次元の撮像画像５０を用いて牛の肉付きを評価するため、牛を三次元に撮像するための特別な設備や牛を上方から撮像するための特別な設備を必要とせず、多様な飼育現場に適用することができる。また、２次元画像である撮像画像５０に牛が大きく写された場合でも、小さく写された場合でも、角度ω１は変化しないことから、角度ω１を用いて牛の肉付きを評価することで、正規化等の処理を行わなくても牛の肉付きを良好に評価することができる。更に、人が牛の肉付きを評価する場合は、判定基準があいまいとなって判定にばらつきが起こるが、本第１の実施形態では、制御部２２が牛の肉付きの判定を行うため、判定のばらつきを抑えることができる。更に、立位にある牛を略真後ろから撮像した二次元の撮像画像５０を用いて牛の肉付きを評価することから、スマートフォン等の携帯性に優れて汎用性の高い評価装置１００を用いることができる。

また、本第１の実施形態では、決定部３８は、図５のように、山部６１ａの第１曲線部６３ａのうち最も上側に張り出している点を、第１曲線部６３ａ上の点４０３ａと決定する。また、山部６１ａの第２曲線部６４ａのうち最も外側に張り出している点を、第２曲線部６４ａ上の点４０４ａと決定する。このような点４０３ａ、４０４ａを用いて算出した角度ω１は、図９のように、ボディコンディションスコアとの相関性が高い。このため、角度ω１を用いて牛の肉付きを評価することで、牛の肉付きを精度良く評価することができる。

また、本第１の実施形態では、決定部３８は、図５のように、水平線４２０を山部６１ａの上方から山部６１ａに向かって近づけて最初に接する点を第１曲線部６３ａ上の点４０３ａと決定する。また、決定部３８は、垂直線４２２ａを山部６１ａの側方から山部６１ａに向かって近づけて最初に接する点を、第２曲線部６４ａ上の点４０４ａと決定する。これにより、第１曲線部６３ａのうち最も上側に張り出している点を点４０３ａとすること、及び、第２曲線部６４ａのうち最も外側に張り出している点を点４０４ａとすることを簡易な処理により行うことができる。

また、決定部３８は、上述したように、輪郭線６０の上部の左右両側に存在する２つの山部６１ａ、６１ｂのうちの山部６１ｂの第１曲線部６３ｂ上の点４０３ｂ（第１の点）と、第２曲線部６４ｂ上の点４０４ｂ（第２の点）と、を決定してもよい（図５参照）。算出部３４は、点４０３ｂと点４０４ｂとを結ぶ線分４１０ｂ（第１線分）の中点４０５ｂを通り線分４１０ｂに垂直な直線４２４ｂが山部６１ｂに交差する点４０６ｂ（第３の点）を特定し、点４０３ｂと点４０６ｂを結ぶ線分４１２ｂ（第３線分）と、点４０４ｂと点４０６ｂを結ぶ線分４１４ｂ（第４線分）と、がなす角度ω２を算出してもよい（図６（ｂ）参照）。そして、評価部３６は、算出部３４が算出した角度ω２から、撮像画像５０に写された牛の肉付きを評価してもよい。この場合でも、角度ω１から牛の肉付きを評価する場合と同じ効果が得られる。

また、決定部３８は、上述したように、輪郭線６０の上部の左右両側に存在する２つの山部６１ａ、６１ｂの両方について、第１曲線部６３ａ、６３ｂ上の点４０３ａ、４０３ｂ（第１の点）と、第２曲線部６４ａ、６４ｂ上の点４０４ａ、４０４ｂ（第２の点）と、を決定してもよい（図５参照）。算出部３４は、山部６１ａにおいて線分４１２ａと線分４１４ａとがなす角度ω１と、山部６１ｂにおいて線分４１２ｂと線分４１４ｂとがなす角度ω２と、の平均値を算出してもよい。評価部３６は、算出部３４が算出した角度ω１と角度ω２の平均値から、撮像画像５０に写された牛の肉付きを評価してもよい。角度ω１と角度ω２の平均値は、図９で説明したように、ボディコンディションスコアとの相関を示す決定係数Ｒ^２が高い。このため、牛の肉付きを精度良く評価することができる。また、角度ω１と角度ω２の平均値を用いることで、牛を略真後ろから撮像するときの撮影角度の影響を受け難くなる。

なお、上記第１の実施形態において、評価部３６は、角度ω１、角度ω２、及び角度ω１と角度ω２の平均値の少なくとも２つからボディコンディションスコアを求め、求めたボディコンディションスコアの平均値を評価してもよい。

なお、牛の肉付きを評価するにあたり、角度ω１を用いるか、角度ω２を用いるか、角度ω１と角度ω２の平均値を用いるかは、ユーザーが適宜選択することができる。

［決定部３８による他の決定方法］
図１２（ａ）から図１２（ｃ）は、決定部３８による他の決定方法を示す図である。図１２（ａ）に示すように、山部６１ａの第２曲線部６４ａの一部が垂直線４２２ａに対し平行となっていて、垂直線４２２ａが第２曲線部６４ａに最初に接する箇所が点ではなく範囲４２６である場合がある。この場合、決定部３８は、範囲４２６の中央の点を点４０４ａと決定する。点４０３ａ、４０３ｂ、４０４ｂについても同様である。

図１２（ｂ）に示すように、山部６１ａの第１曲線部６３ａが頂点６２ａから最上点９１に向かって凹みを形成せずに上がり続ける場合がある。この場合、決定部３８は、第１曲線部６３ａのうちの最も牛の背骨側に位置する点である輪郭線６０の頂点６２ａ側の変曲点を点４０３ａと決定する。点４０３ｂについても同様である。

図１２（ｃ）に示すように、山部６１ａの第２曲線部６４ａが頂点６２ａから下側に向かって凹みを形成せずに下がり続ける場合がある。この場合、決定部３８は、第２曲線部６４ａのうちの最も腹側に位置する点である腰角と腹部の境目の点を点４０４ａと決定する。点４０４ｂについても同様である。

図１２（ａ）から図１２（ｃ）の場合でも、決定部３８は、第１曲線部６３ａのうち最も上側に張り出している点を点４０３ａと決定し、第２曲線部６４ａのうち最も横側に張り出している点を点４０４ａと決定している。したがって、上記第１の実施形態と同様に、牛の肉付きを精度良く評価することができる。

《第２の実施形態》
第２の実施形態に係る牛の評価装置の機能構成およびハードウェア構成は、第１の実施形態の図１および図１０と同じである。第２の実施形態における制御部２２の機能構成は第１の実施形態の図２と同じであるが、算出部３４が実行する処理が第１の実施形態と異なっている。以下においては、算出部３４が実行する処理を中心に説明する。

図１３（ａ）及び図１３（ｂ）は、第２の実施形態における算出部３４が実行する処理を説明するための図である。図１３（ａ）に示すように、決定部３８が山部６１ａの第１曲線部６３ａ上の点４０３ａと第２曲線部６４ａ上の点４０４ａを決定した場合、算出部３４は、点４０３ａと点４０４ａとを結ぶ線分４１０ａの長さＸ１ａ、中点４０５ａと点４０６ａとを結ぶ線分４１６ａの長さＸ２ａ、又は長さＸ２ａに対する長さＸ１ａの長さ比（Ｘ１ａ／Ｘ２ａ）を算出する。算出部３４が算出する長さＸ１ａ、長さＸ２ａ、及び長さ比（Ｘ１ａ／Ｘ２ａ）は、山部６１ａの傾斜の度合いを示す指標値の１つである。

評価部３６は、算出部３４が算出した長さＸ１ａ、長さＸ２ａ、又は長さ比（Ｘ１ａ／Ｘ２ａ）から、撮像画像５０に写された牛の肉付きを評価する。

また、図１３（ｂ）に示すように、決定部３８が山部６１ｂの第１曲線部６３ｂ上の点４０３ｂと第２曲線部６４ｂ上の点４０４ｂを決定した場合では、算出部３４は、点４０３ｂと点４０４ｂとを結ぶ線分４１０ｂの長さＸ１ｂ、中点４０５ｂと点４０６ｂとを結ぶ線分４１６ｂの長さＸ２ｂ、又は長さＸ２ｂに対する長さＸ１ｂの長さ比（Ｘ１ｂ／Ｘ２ｂ）を算出する。算出部３４が算出する長さＸ１ｂ、長さＸ２ｂ、及び長さ比（Ｘ１ｂ／Ｘ２ｂ）は、山部６１ｂの傾斜の度合いを示す指標値の１つである。

評価部３６は、算出部３４が算出した長さＸ１ｂ、長さＸ２ｂ、又は長さ比（Ｘ１ｂ／Ｘ２ｂ）から、撮像画像５０に写された牛の肉付きを評価する。

また、決定部３８が山部６１ａの第１曲線部６３ａ上の点４０３ａと第２曲線部６４ａ上の点４０４ａ、及び、山部６１ｂの第１曲線部６３ｂ上の点４０３ｂと第２曲線部６４ｂ上の点４０４ｂ、の全てを決定した場合では、算出部３４は、長さＸ１ａと長さＸ１ｂの平均値、長さＸ２ａと長さＸ２ｂの平均値、又は長さ比（Ｘ１ａ／Ｘ２ａ）と長さ比（Ｘ１ｂ／Ｘ２ｂ）の平均値を算出してもよい。この場合、評価部３６は、算出部３４が算出した長さＸ１ａと長さＸ１ｂの平均値、長さＸ２ａと長さＸ２ｂの平均値、又は長さ比（Ｘ１ａ／Ｘ２ａ）と長さ比（Ｘ１ｂ／Ｘ２ｂ）の平均値から、撮像画像５０に写された牛の肉付きを評価してもよい。

図１４（ａ）から図１４（ｃ）は、牛の肉付きと、算出部３４が算出する長さＸ１ａ、Ｘ２ａと、の関係を示す図である。図１４（ａ）、図１４（ｂ）、図１４（ｃ）の順に肉付きの良い牛の輪郭線６０を示している。図１４（ａ）から図１４（ｃ）のように、点４０３ａと点４０４ａを結ぶ線分４１０ａの長さＸ１ａ、及び、中点４０５ａと点４０６ａを結ぶ線分４１６ａの長さＸ２ａは、肉付きが良いほど大きな値となることが分かった。また、長さＸ２ａに対する長さＸ１ａの長さ比（Ｘ１ａ／Ｘ２ａ）は、肉付きが良いほど小さな値となることが分かった。これらのことから、長さＸ１ａ、長さＸ２ａ、又は長さ比（Ｘ１ａ／Ｘ２ａ）を求めることで、牛の肉付きを評価できると考えられる。なお、図示は省略するが、牛の肉付きと、算出部３４が算出する長さＸ１ｂ、長さＸ２ｂ、長さ比（Ｘ１ｂ／Ｘ２ｂ）と、の関係についても同じである。

［実験］
獣医師等の専門家が評価したボディコンディションスコアと、算出部３４が算出した値と、の相関を調査する実験を行った。実験は、乳牛（ホルスタイン種）に対して行った。

図１５には、ボディコンディションスコアと、算出部３４が算出した長さＸ１ａ、長さＸ１ｂ、及び長さＸ１ａと長さＸ１ｂの平均値と、の相関を調査した実験結果が示されている。図１５の横軸（ｘ軸）は専門家が評価したボディコンディションスコアであり、縦軸は長さである。図１５中の複数の黒丸は長さＸ１ａを用いたときの乳牛１頭１頭の結果を示していて、点線の直線はそのときの近似直線（回帰直線）を示している。複数の黒四角は長さＸ１ｂを用いたときの乳牛１頭１頭の結果を示していて、一点鎖線の直線はそのときの近似直線（回帰直線）を示している。複数の黒三角は長さＸ１ａと長さＸ１ｂの平均値を用いたときの乳牛１頭１頭の結果を示していて、破線の直線はそのときの近似直線（回帰直線）を示している。図１５に示すように、長さＸ１ａ、長さＸ１ｂ、及び長さＸ１ａと長さＸ１ｂの平均値のいずれを用いた場合であっても、長さが長くなるほどボディコンディションスコアは大きくなっている。長さＸ１ａを用いた点線の近似直線はｙ＝５．８３１４ｘ＋３．８１３５と求まり、決定係数Ｒ^２は０．６０３７であった。長さＸ１ｂを用いた一点鎖線の近似直線はｙ＝５．９７７１ｘ＋２．７０２８と求まり、決定係数Ｒ^２は０．６０８４であった。長さＸ１ａと長さＸ１ｂの平均値を用いた破線の近似直線はｙ＝５．９０４２ｘ＋３．２５８２と求まり、決定係数Ｒ^２は０．６５３７であった。このことから、ボディコンディションスコアと、算出部３４が算出する長さＸ１ａ、長さＸ１ｂ、及び長さＸ１ａと長さＸ１ｂの平均値と、には強い相関関係があることが確認された。

このように、ボディコンディションスコアと、算出部３４が算出する長さＸ１ａ、長さＸ１ｂ、及び長さＸ１ａと長さＸ１ｂの平均値と、には強い相関関係があることから、評価部３６は長さＸ１ａ、長さＸ１ｂ、又は長さＸ１ａと長さＸ１ｂの平均値を用いて牛の肉付きを評価する。例えば、表１と同じような、長さＸ１ａ、長さＸ１ｂ、又は長さＸ１ａと長さＸ１ｂの平均値と、ボディコンディションスコアと、が関連付けられたスコア情報を用いて、牛の肉付きを評価する。これにより、牛の肉付きを精度良く評価することができる。

図１６には、ボディコンディションスコアと、算出部３４が算出した長さＸ２ａ、長さＸ２ｂ、及び長さＸ２ａと長さＸ２ｂの平均値と、の相関を調査した実験結果が示されている。図１６の横軸（ｘ軸）は専門家が評価したボディコンディションスコアであり、縦軸は長さである。図１６中の複数の黒丸は長さＸ２ａを用いたときの乳牛１頭１頭の結果を示していて、点線の直線はそのときの近似直線（回帰直線）を示している。複数の黒四角は長さＸ２ｂを用いたときの乳牛１頭１頭の結果を示していて、一点鎖線の直線はそのときの近似直線（回帰直線）を示している。複数の黒三角は長さＸ２ａと長さＸ２ｂの平均値を用いたときの乳牛１頭１頭の結果を示していて、破線の直線はそのときの近似直線（回帰直線）を示している。図１６に示すように、長さＸ２ａ、長さＸ２ｂ、及び長さＸ２ａと長さＸ２ｂの平均値のいずれを用いた場合であっても、長さが長くなるほどボディコンディションスコアは大きくなっている。長さＸ２ａを用いた点線の近似直線はｙ＝２．０３８８ｘ－２．２０４２と求まり、決定係数Ｒ^２は０．７７０７であった。長さＸ２ｂを用いた一点鎖線の近似直線はｙ＝２．５９３９ｘ－３．９３９２と求まり、決定係数Ｒ^２は０．７２９７であった。長さＸ２ａと長さＸ２ｂの平均値を用いた破線の近似直線はｙ＝２．３１６４ｘ－３．０７１７と求まり、決定係数Ｒ^２は０．８０９８であった。このことから、ボディコンディションスコアと、算出部３４が算出する長さＸ２ａ、長さＸ２ｂ、及び長さＸ２ａと長さＸ２ｂの平均値と、には強い相関関係があることが確認された。

このように、ボディコンディションスコアと、算出部３４が算出する長さＸ２ａ、長さＸ２ｂ、及び長さＸ２ａと長さＸ２ｂの平均値と、には強い相関関係があることから、評価部３６は長さＸ２ａ、長さＸ２ｂ、又は長さＸ２ａと長さＸ２ｂの平均値を用いて牛の肉付きを評価する。例えば、表１と同じような、長さＸ２ａ、長さＸ２ｂ、又は長さＸ２ａと長さＸ２ｂの平均値と、ボディコンディションスコアと、が関連付けられたスコア情報を用いて、牛の肉付きを評価する。これにより、牛の肉付きを精度良く評価することができる。

図１７には、ボディコンディションスコアと、算出部３４が算出した長さＸ２ａに対する長さＸ１ａの長さ比（Ｘ１ａ／Ｘ２ａ）、長さＸ２ｂに対する長さＸ１ｂの長さ比（Ｘ１ｂ／Ｘ２ｂ）、及び長さ比（Ｘ１ａ／Ｘ２ａ）と長さ比（Ｘ１ｂ／Ｘ２ｂ）の平均値と、の相関を調査した実験結果が示されている。図１７の横軸（ｘ軸）は専門家が評価したボディコンディションスコアであり、縦軸は長さ比である。図１７中の複数の黒丸は長さ比（Ｘ１ａ／Ｘ２ａ）を用いたときの乳牛１頭１頭の結果を示していて、点線の直線はそのときの近似直線（回帰直線）を示している。複数の黒四角は長さ比（Ｘ１ｂ／Ｘ２ｂ）を用いたときの乳牛１頭１頭の結果を示していて、一点鎖線の直線はそのときの近似直線（回帰直線）を示している。複数の黒三角は長さ比の平均値（（（Ｘ１ａ／Ｘ２ａ）＋（Ｘ１ｂ／Ｘ２ｂ））／２）を用いたときの乳牛１頭１頭の結果を示していて、破線の直線はそのときの近似直線（回帰直線）を示している。図１７に示すように、長さ比（Ｘ１ａ／Ｘ２ａ）、長さ比（Ｘ１ｂ／Ｘ２ｂ）、及び長さ比（Ｘ１ａ／Ｘ２ａ）と長さ比（Ｘ１ｂ／Ｘ２ｂ）の平均値のいずれを用いた場合であっても、長さ比が小さくなるほどボディコンディションスコアは大きくなっている。長さ比（Ｘ１ａ／Ｘ２ａ）を用いた点線の近似直線はｙ＝－１．０９６７ｘ＋９．０１６４と求まり、決定係数Ｒ^２は０．５６７６であった。長さ比（Ｘ１ｂ／Ｘ２ｂ）を用いた一点鎖線の近似直線はｙ＝－１．５３３８ｘ＋１０．４４１と求まり、決定係数Ｒ^２は０．５９８５であった。長さ比（Ｘ１ａ／Ｘ２ａ）と長さ比（Ｘ１ｂ／Ｘ２ｂ）の平均値を用いた破線の近似直線はｙ＝－１．３１５２ｘ＋９．７２８６と求まり、決定係数Ｒ^２は０．６３５３であった。このことから、ボディコンディションスコアと、算出部３４が算出する長さ比（Ｘ１ａ／Ｘ２ａ）、長さ比（Ｘ１ｂ／Ｘ２ｂ）、及び長さ比（Ｘ１ａ／Ｘ２ａ）と長さ比（Ｘ１ｂ／Ｘ２ｂ）の平均値と、には強い相関関係があることが確認された。

このように、ボディコンディションスコアと、算出部３４が算出する長さ比（Ｘ１ａ／Ｘ２ａ）、長さ比（Ｘ１ｂ／Ｘ２ｂ）、及び長さ比（Ｘ１ａ／Ｘ２ａ）と長さ比（Ｘ１ｂ／Ｘ２ｂ）の平均値と、には強い相関関係があることから、評価部３６は長さ比（Ｘ１ａ／Ｘ２ａ）、長さ比（Ｘ１ｂ／Ｘ２ｂ）、又は長さ比（Ｘ１ａ／Ｘ２ａ）と長さ比（Ｘ１ｂ／Ｘ２ｂ）の平均値を用いて牛の肉付きを評価する。例えば、表１と同じような、長さ比（Ｘ１ａ／Ｘ２ａ）、長さ比（Ｘ１ｂ／Ｘ２ｂ）、又は長さ比（Ｘ１ａ／Ｘ２ａ）と長さ比（Ｘ１ｂ／Ｘ２ｂ）の平均値と、ボディコンディションスコアと、が関連付けられたスコア情報を用いて、牛の肉付きを評価する。これにより、牛の肉付きを精度良く評価することができる。

以上説明したように、本第２の実施形態によれば、算出部３４は、点４０３ａと点４０４ａとを結ぶ線分４１０ａの長さＸ１ａ、中点４０５ａと点４０６ａとを結ぶ線分４１６ａ（第２線分）の長さＸ２ａ、又は長さＸ２ａに対する長さＸ１ａの長さ比（Ｘ１ａ／Ｘ２ａ）を算出する（図１３（ａ）参照）。評価部３６は、算出部３４が算出した長さＸ１ａ、長さＸ２ａ、又は長さ比（Ｘ１ａ／Ｘ２ａ）から、撮像画像５０に写された牛の肉付きを評価する。この場合でも、上記第１の実施形態と同じく、立位にある牛を略後ろから撮像した二次元の撮像画像５０を用いて牛の肉付きを評価するため、多様な飼育現場に適用することができる。また、長さ比（Ｘ１ａ／Ｘ２ａ）を用いて牛の肉付きを評価する場合は、上記第１の実施形態と同じく、正規化等の処理を行わなくても牛の肉付きを良好に評価することができる。更に、制御部２２が牛の肉付きの判定を行うため、判定のばらつきを抑えることができる。更に、立位にある牛を略真後ろから撮像した二次元の撮像画像５０を用いて牛の肉付きを評価することから、スマートフォン等の携帯性に優れて汎用性の高い評価装置を用いることができる。

また、決定部３８が山部６１ｂの第１曲線部６３ｂ上の点４０３ｂ（第１の点）と、第２曲線部６４ｂ上の点４０４ｂ（第２の点）とを決定した場合では、算出部３４は、点４０３ｂと点４０４ｂとを結ぶ線分４１０ｂ（第１線分）の長さＸ１ｂ、中点４０５ｂと点４０６ｂとを結ぶ線分４１６ｂ（第２線分）の長さＸ２ｂ、又は長さＸ２ｂに対する長さＸ１ｂの長さ比（Ｘ１ｂ／Ｘ２ｂ）を算出する（図１３（ｂ）参照）。評価部３６は、算出部３４が算出した長さＸ１ｂ、長さＸ２ｂ、又は長さ比（Ｘ１ｂ／Ｘ２ｂ）から、撮像画像５０に写された牛の肉付きを評価する。この場合でも、長さＸ１ａ、長さＸ２ａ、又は長さ比（Ｘ１ａ／Ｘ２ａ）から牛の肉付きを評価する場合と同じ効果が得られる。

また、決定部３８が山部６１ａ、６１ｂの両方について第１曲線部６３ａ、６３ｂ上の点４０３ａ、４０３ｂ（第１の点）と、第２曲線部６４ａ、６４ｂ上の点４０４ａ、４０４ｂ（第２の点）と、を決定した場合、算出部３４は、長さＸ１ａと長さＸ１ｂの平均値、長さＸ２ａと長さＸ２ｂの平均値、又は長さ比（Ｘ１ａ／Ｘ２ａ）と長さ比（Ｘ１ｂ／Ｘ２ｂ）の平均値を算出してもよい。評価部３６は、算出部３４が算出した長さＸ１ａと長さＸ１ｂの平均値、長さＸ２ａと長さＸ２ｂの平均値、又は長さ比（Ｘ１ａ／Ｘ２ａ）と長さ比（Ｘ１ｂ／Ｘ２ｂ）の平均値から、撮像画像５０に写された牛の肉付きを評価してもよい。長さＸ１ａと長さＸ１ｂの平均値、長さＸ２ａと長さＸ２ｂの平均値、及び長さ比（Ｘ１ａ／Ｘ２ａ）と長さ比（Ｘ１ｂ／Ｘ２ｂ）の平均値は、図１５から図１７で説明したように、ボディコンディションスコアとの相関を示す決定係数Ｒ^２が高い。このため、牛の肉付きを精度良く評価することができる。また、長さＸ１ａと長さＸ１ｂの平均値、長さＸ２ａと長さＸ２ｂの平均値、及び長さ比（Ｘ１ａ／Ｘ２ａ）と長さ比（Ｘ１ｂ／Ｘ２ｂ）の平均値を用いることで、牛を略真後ろから撮像するときの撮影角度の影響を受け難くなる。

なお、長さＸ１ａ、長さＸ１ｂ、長さＸ２ａ、長さＸ２ｂ、長さＸ１ａと長さＸ１ｂの平均値、及び長さＸ２ａと長さＸ２ｂの平均値を用いて牛の肉付きを評価する場合、牛を略真後ろから撮影するときに、長さの基準となる基準物を一緒に撮影したり、同一距離及び同一倍率で撮影したりすることが好ましい。

なお、牛の肉付きを評価するにあたり、長さＸ１ａ、長さＸ１ｂ、長さＸ２ａ、長さＸ２ｂ、長さ比（Ｘ１ａ／Ｘ２ａ）、長さ比（Ｘ１ｂ／Ｘ２ｂ）、長さＸ１ａと長さＸ１ｂの平均値、長さＸ２ａと長さＸ２ｂの平均値、又は長さ比（Ｘ１ａ／Ｘ２ａ）と長さ比（Ｘ１ｂ／Ｘ２ｂ）の平均値、のどれを用いるかは、ユーザーが適宜選択することができる。

なお、上記第１の実施形態及び第２の実施形態において、算出部３４は、山部６１ａに関して、角度ω１、長さＸ１ａ、長さＸ２ａ、及び長さ比（Ｘ１ａ／Ｘ２ａ）の２つ以上を算出してもよい。そして、評価部３６は、角度ω１、長さＸ１ａ、長さＸ２ａ、及び長さ比（Ｘ１ａ／Ｘ２ａ）の２つ以上の値それぞれから得られたボディコンディションスコアの最大値、最小値、又は平均値によって牛の肉付きを評価してもよい。同様に、算出部３４は、山部６１ｂに関して、角度ω２、長さＸ１ｂ、長さＸ２ｂ、及び長さ比（Ｘ１ｂ／Ｘ２ｂ）の２つ以上を算出してもよい。そして、評価部３６は、角度ω２、長さＸ１ｂ、長さＸ２ｂ、及び長さ比（Ｘ１ｂ／Ｘ２ｂ）の２つ以上の値それぞれから得られたボディコンディションスコアの最大値、最小値、又は平均値によって牛の肉付きを評価してもよい。また、評価部３６は、角度ω１と角度ω２の平均値、長さＸ１ａと長さＸ１ｂの平均値、長さＸ２ａと長さＸ２ｂの平均値、及び長さ比（Ｘ１ａ／Ｘ２ａ）と長さ比（Ｘ１ｂ／Ｘ２ｂ）の平均値の２つ以上の値それぞれから得られたボディコンディションスコアの最大値、最小値、又は平均値によって牛の肉付きを評価してもよい。

なお、上記第１の実施形態及び第２の実施形態において、線分４１０ａ、４１０ｂの中点４０５ａ、４０５ｂを用いる場合を例に示したが、線分４１０ａ、４１０ｂの中点から線分４１０ａ、４１０ｂの長さの１／５０以下の範囲に位置する略中点を用いる場合でもよい。直線４２４ａ、４２４ｂ及び線分４１６ａ、４１６ｂは、線分４１０ａ、４１０ｂの垂線である場合を例に示したが、線分４１０ａ、４１０ｂの垂線から±５°の範囲内で傾いた略垂線の場合でもよい。

なお、上記第１の実施形態及び第２の実施形態において、決定部３８は、抽出画像５１をＸ－Ｙ座標系に変換し、抽出画像５１の各画素のＸ、Ｙ座標から第１曲線部６３ａ、６３ｂのうち最も上側に張り出した点を点４０３ａ、４０３ｂと決定し、第２曲線部６４ａ、６４ｂのうち最も側方に張り出した点を点４０４ａ、４０４ｂと決定してもよい。

［水平線４２０及び垂直線４２２ａ、４２２ｂの許容角度］
上記第１の実施形態及び第２の実施形態では、抽出画像５１内の水平線４２０を用いて第１曲線部６３ａ、６３ｂ上の点４０３ａ、４０３ｂを決定する場合を示した。抽出画像５１内の垂直線４２２ａ、４２２ｂを用いて第２曲線部６４ａ、６４ｂ上の点４０４ａ、４０４ｂを決定する場合を示した。ここでは、水平線４２０の代わりに少し傾斜をした略水平線を用いた場合、及び、垂直線４２２ａ、４２２ｂの代わりに少し傾斜をした略垂直線を用いた場合、傾斜角度はどの程度まで許容されるかを調査した。

［実験］
抽出画像５１内の水平線からの傾斜角度を－２０°、－１０°、０°、１０°、２０°とした略水平線を用いて第１曲線部６３ａ、６３ｂ上の点４０３ａ、４０３ｂを決定した。略水平線が左下から右上に伸びている場合の傾斜角度を正とし、右下から左上に伸びている場合の傾斜角度を負とした。また、抽出画像５１内の水平線からの傾斜角度が８０°、９０°、１００°とした略垂直線を用いて第２曲線部６４ａ、６４ｂ上の点４０４ａ、４０４ｂを決定した。傾斜角度が８０°の場合は、略垂直線が左下から右上に伸び、傾斜角度が１００°の場合は、略垂直線が右下から左上に伸びている。

決定した点を用いて、上記第１の実施形態のように、角度ω１と角度ω２の平均値とボディコンディションスコアとの間の決定係数Ｒ^２を求めた。また、決定した点を用いて、上記第２の実施形態のように、長さＸ１ａと長さＸ１ｂの平均値とボディコンディションスコアとの間の決定係数Ｒ^２、長さＸ２ａと長さＸ２ｂの平均値とボディコンディションスコアとの間の決定係数Ｒ^２、長さＸ２ａに対する長さＸ１ａの長さ比（Ｘ１ａ／Ｘ２ａ）と長さＸ２ｂに対する長さＸ１ｂの長さ比（Ｘ１ｂ／Ｘ２ｂ）の平均値とボディコンディションスコアとの間の決定係数Ｒ^２を求めた。実験結果を表２に示す。

表２に示すように、角度ω１と角度ω２の平均値を用いた場合では、傾斜角度が－１０°～２０°の略水平線を用いかつ傾斜角度が８０°～９０°の略垂直線を用いた場合に大きな決定係数Ｒ^２が得られた。長さＸ１ａと長さＸ１ｂの平均値を用いた場合では、傾斜角度が－２０°～２０°の略水平線を用いかつ傾斜角度が９０°の略垂直線を用いた場合に大きな決定係数Ｒ^２が得られた。長さＸ２ａと長さＸ２ｂの平均値を用いた場合では、傾斜角度が－１０°～２０°の略水平線を用いかつ傾斜角度が８０°～９０°の略垂直線を用いた場合に大きな決定係数Ｒ^２が得られた。長さ比の平均値を用いた場合では、傾斜角度が－１０°～２０°の略水平線を用いかつ傾斜角度が９０°の略垂直線を用いた場合に大きな決定係数Ｒ^２が得られた。これらの結果から、抽出画像５１内の水平方向（左右方向）からの傾斜角度が－１０°～２０°の略水平線を用いて第１曲線部６３ａ、６３ｂ上の点４０３ａ、４０３ｂを決定する。抽出画像５１内の水平方向（左右方向）からの傾斜角度が８０°～９０°の略垂直線を用いて第２曲線部６４ａ、６４ｂ上の点４０４ａ、４０４ｂを決定する。これにより、牛の肉付きを精度良く評価することができる。牛の肉付きを精度良く評価する点から、略水平線の傾きは－５°～１０°の場合がより好ましく、－５°～５°の場合が更に好ましい、略垂直線の傾きは８５°～９０°の場合がより好ましい。

なお、上記第１の実施形態及び第２の実施形態では、画像取得部３０は、立位にある牛を略真後ろから撮像した撮像画像を取得する場合を例に示したが、立位にある牛を略正面から撮像した撮像画像を取得する場合でもよい。図１８（ａ）は、立位にある牛を略真後ろから撮像した撮像画像５０、図１８（ｂ）は、図１８（ａ）の牛を略正面から撮像した撮像画像５０ａである。図１８（ａ）に示すように、立位にある牛を略真後ろから撮像した場合に限られず、図１８（ｂ）に示すように、立位にある牛を略正面から撮像した場合であっても、牛の左右の腰角周りの輪郭が写るように牛を撮像することができる。

図１９（ａ）及び図１９（ｂ）は、図１８（ａ）及び図１８（ｂ）の撮像画像５０、５０ａから牛の輪郭線６０を抽出した抽出画像５１である。図１９（ａ）に示すように、立位にある牛を略真後ろから撮像した場合であっても、図１９（ｂ）に示すように、立位にある牛を略正面から撮像した場合であっても、輪郭線６０の上部の左右両側に腰角に対応する山部６１ａ、６１ｂを存在させることができる。よって、立位にある牛を略正面から撮像した撮像画像５０ａを用いても、略真後ろから撮像した撮像画像５０を用いる場合と同じ方法によって、牛の肉付きを評価することができる。また、立位にある牛を略正面から撮像した二次元の撮像画像５０ａを用いる場合でも、略真後ろから撮像した撮像画像５０を用いる場合と同様に、牛を上方から撮像するための特別な設備を必要とせず、多様な飼育現場に適用することができるとともに、携帯性に優れて汎用性の高いスマートフォン等の携帯型情報機器を評価装置とすることができる。

なお、上記第１の実施形態及び第２の実施形態では、牛の肉付きを評価する評価装置が、作業者が携帯して利用される、撮像部１８を備える携帯型情報機器である場合を例に示したが、その他の場合でもよい。例えば、ネットワークに接続された他の機器と通信可能なサーバ等の据置型情報機器の場合でもよい。この場合、作業者は例えばスマートフォン等の携帯型情報機器を用いて撮像した牛の撮像画像を、ネットワークを介してサーバ等の据置型情報機器に送信する。据置型情報機器では、評価装置１００と同様の処理を行って牛の肉付きを評価する。そして、据置型情報機器は、牛の肉付きの評価結果を、ネットワークを介して作業者が携帯する携帯型情報機器に送信する。これにより、スマートフォンに掛かる負荷を低減できる。また、牛舎内に固定又は移動可能にカメラを設置し、スマートフォン等の携帯型情報機器及びサーバ等の据置型情報機器は、牛舎内に設置されたカメラが撮像した立位にある牛の略真後ろ又は略正面の撮像画像を取得してもよい。

なお、上記第１の実施形態及び第２の実施形態では、牛の肉付きとしてボディコンディションスコアを評価する場合を例に示したが、牛の肉付きを評価するその他の指標を評価する場合でもよい。また、各実験は、乳牛に対して行ったが、同様の骨格を有する肉牛も同様の結果になると考えられる。したがって、上記各実施形態は乳牛だけでなく肉牛にも適用できる。

以上、本願発明の実施形態について詳述したが、本願発明はかかる特定の実施形態に限定されるものではなく、特許請求の範囲に記載された本願発明の要旨の範囲内において、種々の変形・変更が可能である。

１２ 表示部
１４ 操作部
１６ 通信部
１８ 撮像部
２０ 記憶部
２２ 制御部
３０ 画像取得部
３２ 抽出部
３４ 算出部
３６ 評価部
３８ 決定部
５０、５０ａ 撮像画像
５１ 抽出画像
５２ 左上端点
５４ 右上端点
５５、５５ａ、５５ｂ、５６、５６ａ、５６ｂ 傾斜直線
６０ 輪郭線
６１ａ、６１ｂ 山部
６２ａ、６２ｂ 頂点
６３ａ、６３ｂ 第１曲線部
６４ａ、６４ｂ 第２曲線部
９１ 最上点
１００ 評価装置
４０３ａ、４０３ｂ 点
４０４ａ、４０４ｂ 点
４０５ａ、４０５ｂ 中点
４０６ａ、４０６ｂ 点
４１０ａ、４１０ｂ 線分
４１２ａ、４１２ｂ 線分
４１４ａ、４１４ｂ 線分
４１６ａ、４１６ｂ 線分
４２０ 水平線
４２２ａ、４２２ｂ 垂直線
４２４ａ、４２４ｂ 直線
４２６ 範囲

図８（ａ）から図８（ｃ）は、牛の肉付きと、算出部３４が算出する角度ω１と、の関係を示す図である。図８（ａ）は、痩せている牛の輪郭線６０を示し、図８（ｂ）は、図８（ａ）の牛よりも太っている（肉付きの良い）牛の輪郭線６０を示し、図８（ｃ）は、図８（ｂ）の牛よりも太っている（肉付きの良い）牛の輪郭線６０を示している。図８（ａ）から図８（ｃ）のように、肉付きが良いほど、腰角に対応する部分である山部６１ａが丸みを帯びた形状となる。このため、角度ω１は、肉付きが良いほど小さな値となることが分かった。なお、図示は省略するが、牛の肉付きと、算出部３４が算出する角度ω２と、の関係についても同じである。これらのことから、角度ω１、角度ω２、又は角度ω１と角度ω２の平均値を求め、これらを用いて牛の肉付きを評価できると考えられる。

図１４（ａ）から図１４（ｃ）は、牛の肉付きと、算出部３４が算出する長さＸ１ａ、Ｘ２ａと、の関係を示す図である。図１４（ａ）は、痩せている牛の輪郭線６０を示し、図１４（ｂ）は、図１４（ａ）の牛よりも太っている（肉付きの良い）牛の輪郭線６０を示し、図１４（ｃ）は、図１４（ｂ）の牛よりも太っている（肉付きの良い）牛の輪郭線６０を示している。図１４（ａ）から図１４（ｃ）のように、点４０３ａと点４０４ａを結ぶ線分４１０ａの長さＸ１ａ、及び、中点４０５ａと点４０６ａを結ぶ線分４１６ａの長さＸ２ａは、肉付きが良いほど大きな値となることが分かった。また、長さＸ２ａに対する長さＸ１ａの長さ比（Ｘ１ａ／Ｘ２ａ）は、肉付きが良いほど小さな値となることが分かった。これらのことから、長さＸ１ａ、長さＸ２ａ、又は長さ比（Ｘ１ａ／Ｘ２ａ）を求めることで、牛の肉付きを評価できると考えられる。なお、図示は省略するが、牛の肉付きと、算出部３４が算出する長さＸ１ｂ、長さＸ２ｂ、長さ比（Ｘ１ｂ／Ｘ２ｂ）と、の関係についても同じである。

Claims (8)

1. 立位にある牛を略真後ろ又は略正面から撮像した二次元の撮像画像を取得する画像取得部と、
   前記撮像画像から、前記牛の輪郭線を抽出する抽出部と、
   前記輪郭線の上部の左右両側に存在する２つの山部のうちの一方の山部の前記輪郭線上の略中央の点から他方の山部側に伸びる第１曲線部上の第１の点と前記第１曲線部とは反対側に伸びる第２曲線部上の第２の点とを決定する決定部と、
   前記第１の点と前記第２の点とを結ぶ第１線分の長さ、前記第１線分の略中点を通り前記第１線分に略垂直な直線が前記一方の山部と交差する第３の点と前記第１線分の略中点とを結ぶ第２線分の長さ、前記第２線分の長さに対する前記第１線分の長さの比、及び前記第１の点と前記第３の点を結ぶ第３線分と前記第２の点と前記第３の点を結ぶ第４線分とがなす角度、の少なくとも１つを算出する算出部と、
   前記算出部が算出した値から、前記牛の肉付きを評価する評価部と、を備える牛の評価装置。
2. 前記決定部は、前記第１曲線部のうち最も上側に張り出した点を前記第１の点と決定し、前記第２曲線部のうち最も横側に張り出した点を前記第２の点と決定する、請求項１に記載の牛の評価装置。
3. 前記決定部は、前記撮像画像内で略左右方向に伸びる第１直線を前記一方の山部の上方から前記一方の山部に向かって近づけて最初に接する点を前記第１の点と決定し、前記撮像画像内で略上下方向に伸びる第２直線を前記一方の山部の側方から前記一方の山部に向かって近づけて最初に接する点を前記第２の点と決定する、請求項１または２に記載の牛の評価装置。
4. 前記決定部は、前記２つの山部のうちの他方の山部についての、前記第１曲線部上の前記第１の点と前記第２曲線部上の前記第２の点とを決定し、
   前記算出部は、他方の山部についての、前記第１線分の長さ、前記第２線分の長さ、前記第２線分の長さに対する前記第１線分の長さの比、及び前記第３線分と前記第４線分とがなす角度、の少なくとも１つを算出するとともに、前記２つの山部それぞれの前記第１線分の長さの平均値、前記２つの山部それぞれの前記第２線分の長さの平均値、前記２つの山部それぞれの前記第２線分の長さに対する前記第１線分の長さの比の平均値、及び前記２つの山部それぞれの前記第３線分と前記第４線分とがなす角度の平均値、の少なくとも１つを算出する、請求項１から３のいずれか一項に記載の牛の評価装置。
5. 前記２つの山部は前記牛の腰角に対応する部分である、請求項１から４のいずれか一項に記載の牛の評価装置。
6. 前記評価部は、前記牛の肉付きの評価としてボディコンディションスコアを評価する、請求項１から５のいずれか一項に記載の牛の評価装置。
7. 前記牛を撮像する撮像部を備える、請求項１から６のいずれか一項に記載の牛の評価装置。
8. 立位にある牛を略真後ろ又は略正面から撮像した二次元の撮像画像から前記牛の輪郭線を抽出し、
   抽出した前記輪郭線の上部の左右両側に存在する２つの山部のうちの一方の山部の前記輪郭線上の略中央の点から他方の山部側に伸びる第１曲線部上の第１の点と前記第１曲線部とは反対側に伸びる第２曲線部上の第２の点とを決定し、
   決定した前記第１の点と前記第２の点とを結ぶ第１線分の長さ、前記第１線分の略中点を通り前記第１線分に略垂直な直線が前記一方の山部と交差する第３の点と前記第１線分の略中点とを結ぶ第２線分の長さ、前記第２線分の長さに対する前記第１線分の長さの比、及び前記第１の点と前記第３の点を結ぶ第３線分と前記第２の点と前記第３の点を結ぶ第４線分とがなす角度、の少なくとも１つを算出し、
   算出した値から、前記牛の肉付きを評価する、
   処理をコンピュータが実行することを特徴とする牛の評価方法。
   

Images