無線位置システムによる放飼鶏の軌跡データ取得とその空間分析
岡部篤行 佐藤俊明 岡部佳世 今村栄二 Siripun MORATHOP Charan JAILANGKA Surachai RATANASERMPONG 林良博 秋篠宮文仁
Acquisition of spatio-temporal data of free-range chicken using a Wireless Fidelity (WiFi) positioning system and spatial analysis of the obtained data
Atsuyuki OKABE, Toshiaki SATOH, Kayo OKABE, Eiji IMAMURA, Siripun MORATHOP, Charan JAILANKA, Surachai RATANASERMPONG, Yoshihiro HAYASHI and Fumihito, AKISHINONOMIYA
Abstract: The objective of this study is to acquire spatio-temporal data of free-range chicken kept in an experiment field in Chiang Rai, Thailand, using a Wireless Fidelity (WiFi) positioning system First, the paper describes the functions of the components of the system, the arrangement of the system components in the experiment field, and the procedure of the experiment. Second, the paper obtains the probability density function of observed locations for a fixed point. Third, based on this distribution, the daily life areas of chicken are estimated with the kernel density method, and findings obtained from this estimation are reported.
Keywords: 無線位置システム (WiFi positioning system), 放飼鶏(free-range chicken), 軌跡 追跡(tracking), 精度(accuracy), カーネル密度法 (kernel density method)
1. はじめに
動物生態学で一番の労力と時間を要する作業は,
動物の空間行動軌跡データをある程度の日にちに亘 って取得することである (Turchin, 1998).動物生態 学においては,主に二つの方法が頻繁に使われてい る.一つがGPSを使う方法であり(例えば,報告書
『GPSで動物を追跡する方法』,2001を参照),もう 一 つ がテ レメ ト リー を使 う 方法 であ る (例 えば Millsoaugh and Marzluff, 2001を参照).どちらの方法 も,それを利用する環境条件によって,長所と短所 がある.
GPSについては,重さが一つの問題点となる.研 究当初,知りうる限りで一番軽いものが65グラムで あったが,Ando and Ogasawara (1970)によれば,動物 に装着して動物の行動の自由を奪わない装着物の重
さは,個体の体重の約 5%ということである.従っ て,GPS を装着できる動物は,個体の体重が 1kg 以上のものに限られる.さらにGPSは,動物が木の 下や藪の中などに入ると,その位置データを取得す るのが著しく困難になる.
一方,テレメトリーであるが,研究当初,知りう る限りでは0.5グラム程度のものも販売されており,
かなり体重の少ない動物にも装着することが可能で ある.しかしテレメトリーは,位置データを取得す るのに常に調査者がアンテナを操作していなければ ならないという制約がある.この制約は,昼夜に亘 って多くの動物個体を同時に調査する場合,大きな 障害になる.
そこで当研究では,これらGPSとテレメトリーの 短所を克服できそうなWireless Fidelity(WiFi)位置シ ステムを実験的に放飼鶏群に適用し,実際に短所を 克服できるかどうかを調べた.
岡部 篤行:113-8656 東京都文京区本郷7-3-1
東京大学空間情報科学研究センター 同大大学院工学系研究科 電話:03-5841-6225 e-mail:[email protected]
論文は,第1にWiFi位置システムの適用方法を述 べ,第2にそのシステムで得られるデータの精度を 示す.第3に,その精度の下でのニワトリ空間行動 を把握する方法を検討し,第4にその方法を適用し てニワトリ空間行動の分析を行う.最後に,当研究 をまとめながら WiFi 位置システムの長所と短所を 述べる.
2.実験の方法
2.1. WiFi位置システムの構成と機能
当研究ではAeroScout社が開発したWiFi位置情報 システムを用いた.このシステムは,タグ,アクテ ィベータ,位置受信機,Power over Ethernet (PoE) ハ ブ,WiFiアクセスポイント,管理エンジンの6つの 機器で構成され,主な機器の機能は次の通りである.
タグは動物(今回の実験ではニワトリ)に装着さ れるもので,位置点データ(x-y座標)を無線で送 信する.タグの重さは35グラムで,大きさは62mm
×40mm×17mmである(図1).
位置受信機は,タグから送信されてきた位置デー タを受信し,その受信信号を管理エンジンに送る機 能を持つ.受信機のアンテナには指向性の違いで,
360度領域,135度領域,60度領域の3種類のアン テナがある.
管理エンジンには WiFi アクセスポイント経由で 送られてくる位置データを処理し,タグの位置を表 示する機能を持つ.
位置データは,1 メートルグリッド上の離散点で
与えられるので,解像度は1メートル以下とはなり えず(実際の解像度は後述),測定の時間間隔は,最
短で0.125秒であるが,当研究では1秒に設定した.
2.2. 実験場所とWiFi位置システムの配置
実験を行った場所は,タイの北部チェンラーイ県 にあるチェンラーイ畜産研究技術移転センターの一 画で,広さは約200メートル四方の開墾地である(図 2).この地区に8軒のコンクリートブロック造りの 平屋があり(図2),南の2軒(図2のH1, H5)は居住 家屋で,残りの6軒は空屋となっている.そのうち の2軒(図2のH2, H6)をそれぞれ実験機器の準備 室,管理エンジン室とした.
実験地は,道路を挟んで東側は開墾地で裸地と草 地が広がり,周辺は藪で囲まれている(図3).西側 は,丈の低い木とバナナがまばらに植えられた畑地 である.人の出入りは,2 軒の住人だけであり,日 中は全員が農作業に出るため,人影は極めて少ない.
ニワトリは自由に動き回ることができる地区となっ ている.
図 2 実験場における施設,機器の配置
100m
0 50
人家(H: House)
鶏小屋(CH: Chicken House)
餌場 道路
360°受信機 135°受信機 60°受信機 N
H8 H7
H6 H5
H4 H3 H2
H1 CH1
CH2 CH3
H5の横の木 東の木 H5床下
図3 実験場の東側から望んだ景観 図 1 ニワトリの背中に装着したタグ
タグ
WiFi位置システムは,13個の受信機,2個のPoE ハブ,1個のアクセスポイント,1台の管理エンジン で構成し,受信機は図2のように配置した.
2.3. ニワトリ
被験動物はタイ語で Kai Chon と呼ばれるシャモ 系の18羽のニワトリである.これらのニワトリは,
チェンラーイ畜産研究技術移転センターで飼われて いたニワトリであるが,2005年10月15日から実験 場の図2に示された3つの鶏小屋で放し飼いにせず 飼育した.それぞれの鶏小屋には雄1羽と雌5羽が 飼われ,月齢とひな連れのあり・なしは表1の通り である.
これら18羽のニワトリの背中に図1で示すように タグを装着した.装着には,ニワトリがテープをつ い て ほ つ れ な い よ う に 縫 い 目 が 一 切 な い Bally
Ribbon Mills 社のテフロンテープを使用した.この
テープを含むタグの重さは40グラムとなった.装着 直後,ニワトリは違和感を感じるようであったが,1 時間もするとタグの装着を気にしていないような振 る舞いに見えた.
2.4. 実験の経緯
実験は2005年11月2日から9日にかけて行った.
11月2日にWiFi位置システムの設置作業を行い,
ニワトリにタグを装着し,鶏小屋に戻した.
3 日朝に鶏小屋の扉を開放し,放し飼い状態とし た.3日から4日の午後4時30分までは,WiFi位置 システムの調整を行った.この間,データの取得は 行ったが,アンテナの本数,位置などを変えている ので,データを比較するには,それなりの調整が必 要となる.
4日午後4時30分から9日の午前9時15分まで は WiFi 位置システムを同一状態に固定してデータ の取得をすることができた.ただし,完全なデータ が取得できているかのチェックとハードウェアトラ ブルにより万が一にも取得済みのデータを損なわな いようにするため,毎日,ニワトリが寝床につきお
わっている午後6時頃から約30分間,システムを止 めて,静的なバックアップをとることとした.
9日午前9時15分以降,タグの取り外しとWiFi 位置システムの撤収を行った.
餌は,5日まで午前と午後に1回ずつ,鶏小屋の 中と,図2の▲の位置に置いた.6日から9日まで は,餌を与えなかった.ただし6日には前日の餌が 残っていた.
3.固定点の位置精度
タグの位置精度を解析するために,3 軒の鶏小屋 の屋根にニワトリに装着したタグと同じものをそれ
ニワトリID 飼われていた
小屋番号 雄雌 月齢 ひな連れ
CH1̲C1 CH1 雄 12 なし
CH1̲H1 CH1 雌 6〜7 なし
CH1̲H2 CH1 雌 6〜7 なし
CH1̲H3 CH1 雌 6〜7 なし
CH1̲H4 CH1 雌 6〜7 なし
CH1̲H5 CH1 雌 6〜7 なし
CH2̲C1 CH2 雄 12 なし
CH2̲H̲M CH2 雌 6〜7 あり
CH2̲H2̲Y CH2 雌 2〜3 なし
CH2̲H3̲Y CH2 雌 2〜3 なし
CH2̲H4 CH2 雌 6〜7 なし
CH2̲H5 CH2 雌 6〜7 なし
CH3̲C1 CH3 雄 12 なし
CH3̲H1 CH3 雌 6〜7 なし
CH3̲H2 CH3 雌 6〜7 なし
CH3̲H3 CH3 雌 6〜7 なし
CH3̲H4̲M CH3 雌 6〜7 あり
CH3̲H5 CH3 雌 6〜7 なし
表 1 放飼鶏の所属小屋と特性
*ニワトリIDの命名規則は,最初の3文字が「小屋番号」
を,次の1文字のC(COCK)が「雄鶏」を,H(HEN)が
「雌鶏」を表し,次の一文字は小屋内でのユニーク番号を 示す.また一部に付いている M(MOTHER)および Y
(YOUNG)は,それぞれ「雛連れあり」,「若鶏」を示す.
H5 H1 H6 H2 H7 H3
図4 鶏小屋の固定点の位置点データ分布(分 布点の中にある小さい□は真の位置で,△は 平均位置を示す.)
ぞれ1個ずつ設置した.
図4は,11月7日の3固定点のシステムで取得し た位置点の分布である.この図から分かるように,
固定点でもシステムに起因するランダム要因によっ てばらつきがあり,タグの位置点は確率的である.5 日から8日のデータを同様に調べた結果,次の特性 があることが分かった.
i. 真の位置(図 4 の□)と分布の平均位置(図 4 の△)が一致していないことから,位置点データ にはバイアスがある.
ii. バイアスの方向は,場所によって異なる.
iii. 位置点の分布は真の位置を中心として等方性が ない.非等方性の度合いは3地点によって異なる.
iv. 3地点の分布形は,地点ごとに異なる.同一地点 では,日々似てはいるものの,完全に同じという わけではない.
4. カーネル関数の推定
以上の位置点分布の特性から,取得したニワトリ の位置データは,確率的に扱うのが適当であると言 える(Okabe et al., 2006).具体的には,データを基 にカーネル密度法で位置を推定するのが一つの方法 である.特性iiiからすると,カーネル関数は,非等 方性の関数形が考えられるが,3 地点での分布形か ら実験地の任意の地点のカーネル関数を推定するの は困難である.そういう状況の下では,等方性のカ ーネル関数を使うのがセカンドベストと言えよう.
カーネル関数としては多用されている bi-weight 関数を採用した(Silverman, 1986).そのバンド幅は,
3地点の真の位置を合わせ3地点の分布を集約した 分布より95%を含む半径から求めた.その結果,バ ンド幅は5.63メートルとなった.
5.ニワトリの一日行動圏による分析
図 5 に,カーネル密度法で推定した 3 羽の雄鶏 CH1_C1, CH2_C1, CH3_C1と雌鶏CH_H1の午前6 時から午後4時までの位置点密度分布を示す.グレ ーで塗られている領域は,密度が 0.1以上の領域で
ある.この領域を,「一日行動圏」と呼ぶことにする.
図5と同様な図を全てのニワトリについて求めた.
その結果,それぞれのニワトリの一日行動圏と位置 点密度分布は相違が視覚で明確に判断できることが 判明した.そこで,ニワトリの一日行動圏・位置点 密度分布が似ているものをグループ化し,そのグル ープを表2に示した.この表をもとに主な結果をま とめると次のようになる.
i. 3羽の雄鶏の一日行動圏は,ある程度重なるもの の,異なった形状になっている(図5a,b,c).
ii. それぞれの雄鶏の一日行動圏とほぼ同じ雌鶏が いる(雄鶏に「従う」雌鶏がいる)(CH1_C1 と CH1_H1〜H5,CH2_H4;CH2_C1とCH2_H5; CH3_C1とCH3_H1〜H3,H5).
iii. 若鶏は,雄鶏と異なった一日行動圏を持っている.
iv. ひな連れ雌鶏は,雄鶏と異なった一日行動圏を持 っている.ひな連れ雌鶏CH_H3_Mは,6〜7日 は若鶏と似た一日行動圏をし,8日には異なる一 日行動圏となった.
v. 日中のグループが必ずしも同じところで寝ると はかぎらない.CH2_H4,H5は,夜CH2_C1と異 なる所で寝ている.また日中の異なるグループが 同じ所に寝ることがある.CH3_H4_Mは夜G4グ
(a) (b)
(c) (d)
図 5 雄 鶏 (a) CH1_C1 (b) CH2_C1, (c) CH3_C3,雌鶏(d) CH1_H1の11月7日の一日 行動圏と位置点密度分布
表2 日中におけるグループと夜間の寝場所
ニワトリID 2日まで飼わ
れていた小屋 5日昼 5日夜 6日昼 6日夜 7日昼 7日夜 8日昼 8日夜
CH1̲C1 CH1 G 1 CH 1 G 1 CH 1 G 1 CH 1 G 1 CH 1
CH1̲H1 CH1 G 1 CH 1 G 1 CH 1 G 1 CH 1 屋 根 G 1 CH 1 屋 根
CH1̲H2 CH1 G 1 CH 1 G 1 CH 1 G 1 CH 1 G 1 CH 1
CH1̲H3 CH1 G 1 CH 1 G 1 CH 1 G 1 CH 1 G 1 CH 1
CH1̲H4 CH1 G 1 CH 1 G 1 CH 1 G 1 CH 1 G 1 CH 1
CH1̲H5 CH1 G 1 CH 1 G 1 CH 1 G 1 CH 1 G 1 CH 1
CH2̲C1 CH2 G 2 H 5 床 下 G 2 H 5 床 下 G 2 H 5 床 下 G 2 H 5 床 下
CH2̲H1̲M CH2 G 3 CH 2 G 3 CH 2 G 3 CH 2 G 6 CH 2
CH2̲H2̲Y CH2 G 3 CH 2 G 3 CH 2 G 3 CH 2 G 3 CH 2
CH2̲H3̲Y CH2 G 3 CH 2 G 3 CH 2 G 3 CH 2 G 3 CH 2
CH2̲H4 CH2 G 1 H 5 横 の 木 G 1 H 5 横 の 木 G 1 H 5 横 の 木 G 1 H 5 横 の 木
CH2̲H5 CH2 G 2 H 5 横 の 木 G 2 H 5 横 の 木 G 2 H 5 横 の 木 G 2 東 の 木
CH3̲C1 CH3 G 4 CH 3 G 4 CH 3 G 4 CH 3 G 4 CH 3
CH3̲H1 CH3 G 4 CH 3 G 4 CH 3 G 4 CH 3 G 4 CH 3
CH3̲H2 CH3 G 4 CH 3 G 4 CH 3 G 4 CH 3 G 4 CH 3
CH3̲H3 CH3 G 4 CH 3 G 4 CH 3 G 4 CH 3 G 4 CH 3
CH3̲H4̲M CH3 G 5 CH 3 G 5 CH 3 G 5 CH 3 G 5 CH 3
CH3̲H5 CH3 G 4 CH 3 G 4 CH 3 G 4 CH 3 G 4 CH 3
*表内の「G」は日中の行動グループを示す.
ループのニワトリと共にCH3で寝ている.
vi. ニワトリは実験開始日まで飼われていた鶏小屋 で寝るとはかぎらない.観察よりCH2_C1はH5 の床下(図2の+)で,CH2_H4はH5の横の木 の上(図2の☆)で,CH2_H5はH5の横の木の 上と東の木の上(図2の☆および△)で寝ており,
CH1_H1 は鶏小屋 CH1 から「閉め出されて」そ
の小屋の屋根の上で寝ている.
vii. 餌を断つと一日行動圏は広がるが,グループによ
って程度の差がある.特に若鶏の一日行動圏は拡 大率が大きい(図6).
6.分析方法の検討
一日行動圏は,視覚的な比較でニワトリ空間行動 の特色をかなり明らかにすることができるが,この 視覚的判断をより客観的にするには,いくつかの方 法が考えられる.例えば,2 羽の一日行動圏の形状
を重ね,そのズレで2羽の一日行動圏の近さを計量 化するのも 1方法である.また,2羽の密度分布の 違いを,Kullback-Leibler統計量で計量化することも できよう(但し,密度ゼロの扱いを工夫する必要が ある).このような方法で先に述べた結果をより計量 的に述べることが可能となるが,定性的な結果は変 わらないように思われる.
一日行動圏は,ニワトリの位置の1日に亘る集約 データであるということに注意する必要がある.一 日行動圏が重なるということは,同時間に近い位置 にいるということを意味しない.ということは,一 日行動圏から同じグループに属するとされたニワト リが,時系列的に見ると異なったグループで行動し ているということはあり得る.今後,時系列分析を 加える必要がある.
7.おわりに
WiFi位置システムの適用には,アンテナの位置の 誤認や,測量の座標系の読み間違い,アンテナの再 設置作業などがあり,最終形態で定常的データを得 られるようになるまでにはほぼ 3日を要した(第2 節).しかし一旦,システムを設置すると,位置デー タは人手を借りずに18羽の行動軌跡を昼夜連続(バ ックアップ時を除く)して1秒間隔で取得すること ができた.この労力削減は,準備の労力を上回るも
図6 若鶏(CH2_H3_Y)の一日行動圏
(a)5日,(b)8日
(a) (b)
のがある.
当論文の成果を要約すると,次のようになる.
WiFi 位置システムから得られる位置データは確 率的である.そこで当論文では,定点のデータから 得られる確率密度関数の特性(第3節)を吟味し(第 4節),それからカーネル関数を設定し,カーネル密 度法で位置点密度を推定する方法を採用した.
この位置点密度から導出される一日行動圏を使っ て,ニワトリの行動領域の分析を行った結果,第 5 節に述べたニワトリの空間行動の特色を得ることが できた.
第6節に述べたように,当論文で提案した方法は,
ニワトリの1日の行動位置を時間軸で集約したデー タに基づく分析で,時系列的軌跡のデータに基づい た分析ではない.結果の解釈においては,その限界 をわきまえる必要がある.
最後に現時点での WiFi 位置システムの短所を述 べておこう.図2に示したように,今回,位置を測 定できる範囲は約200メートル四方である.この範 囲をさらに広げることは,ある程度までは可能であ るが,所詮,受信機が設置されている範囲でしか位 置を測定できない.その範囲を越えて行動する動物 には適用をすることが出来ない.また,受信機の設 置は人手で行うので,人が入れる場所でなければな らず,人が容易に行けない場所を動き回る動物には 適用がむずかしい.
WiFi 位置システムの野生動物行動把握への適用 にあたっては,以上の長所と短所を踏まえて行うの が肝要と思われる.
謝辞
当 研 究 は Human-Chicken Multi-relationship Research (HCMR) Project (代表:秋篠宮文仁)の一環と して行われたものである.当研究の一部は,科学研 究費補助金(基盤A)「鶏の形態嗜好に関する日本と タ イ の 多 面 的 比 較 感 性 モ デ ル 」( 代 表 : 原 田 昭,
No.17200015)と家禽資源研究会の補助を受けており,
感謝する次第である.
当研究は多くの方々の共同作業によって行われた ものである.HCMR日本メンバーの赤木攻,原田昭,
川嶋舟,高松淳,Chiang Rai Livestock Research and Technology Transfer Center の C. Jailangka, N.
Mahakanta, S. Trimanee, V. Chintanawat,NECネッツ エスアイの中川貴之,松下和弘,長野一博,石渡祥 嗣,飴本幸司,NESIC THAIの小田浩次, 福安智亮, B.
Kumkhean, J. Krutmoung, PASCO の 佐 藤 圭, K.
Khanophet,Mongkon Surveying & Consulting の M.
Krachangtoy, V. Boonlanasri, C. Kanobdee, M. Viennok, DISTDA のT. Suepa,またI. Rittaporn, 以上の諸氏 に様々なご指導,ご支援を得た.ここに感謝の意を 表する.
参考文献
Ando, S. and Ogasawara, A. (1970) Telemetry of the behavior of Pale thrush(Turdus pallidus). Japanese Journal of Ecology , 20, 137-144.
Millspaugh, J. and Marzluff , J. (2001) Radio Tracking and Animal Populations, New York: Academic Press.
Atsuyuki Okabe, Toshiaki Satoh, Kayo Okabe, Takayuki Nakagawa, Eiji Imamura, Kazuhiro Matsushita, Kazuhiro Nagano, Yoshitsugu Ishiwatari, Koji Amemoto, Yoshihiro Hayashi and Akishinonomiya Fumihito (2006) Applicability of a Wireless Fidelity Positioning System to Tracking Free-range Chickens (Gallus gallus domesticus) and Helmeted Guineafowl (Numida meleagris galeata) Journal of the Yamashina Institute for Ornithology (to appear).
Tracking Animals with GPS, An International Conference held at the Macaulay Land Use Institute, Abadeen, 12-13 March 2001.
Turchin, P. (1998) Quantitative Analysis of Movement, Sunderland: Sinauer Associates.
Silverman, B.W. (1986). Density Estimation for Statistics and Data Analysis, London: Chapman and Hall.