JP4482398B2 - Automatic mortality determination method and mortality automatic determination device by thermograph - Google Patents
Automatic mortality determination method and mortality automatic determination device by thermograph Download PDFInfo
- Publication number
- JP4482398B2 JP4482398B2 JP2004236134A JP2004236134A JP4482398B2 JP 4482398 B2 JP4482398 B2 JP 4482398B2 JP 2004236134 A JP2004236134 A JP 2004236134A JP 2004236134 A JP2004236134 A JP 2004236134A JP 4482398 B2 JP4482398 B2 JP 4482398B2
- Authority
- JP
- Japan
- Prior art keywords
- surface area
- mortality
- area
- calculated
- image
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Radiation Pyrometers (AREA)
- Image Processing (AREA)
Description
本発明は、サーモグラフによる死亡率自動判定方法および死亡率自動判定装置に関するものである。 The present invention relates to a mortality rate automatic determination method and a mortality rate automatic determination device using a thermograph.
大規模な養鶏場は多数（例えば数万羽）の肉用鶏（ブロイラ）を複数の飼育ケージに分けてそれぞれ自動給餌・自動給水などの自動化システムで飼育している。 Large-scale poultry farms divide a large number (for example, tens of thousands) of meat chickens (broilers) into a plurality of breeding cages and raise them using automated systems such as automatic feeding and automatic feeding.
ブロイラは動作が鈍く、給餌や給水時にその場所に歩行して餌を食べたり、水を飲んだりし、それ以外はほとんど横たわっている。そのため、管理者は見回り時などに一羽一羽、生存を確認して死亡している数を把握することは困難であり、強いて生存を確認しようとすると飼育ケージの中に入って一羽一羽、例えば間近で呼吸による微動作をしていない鶏を確認したり、時としてつついて動かないことを確認したりし、その生存を確認して死亡した数を把握する必要があった。 Broilers are sluggish and walk to the place when they are feeding or drinking water to eat or drink water, and otherwise lie down. For this reason, it is difficult for an administrator to check the number of deaths by checking the survival of each individual at the time of patrol and so on. It was necessary to check the number of deaths by checking their wings, for example, chickens that were not moving finely due to respiration, and that they sometimes confirmed that they did not move.
また、赤外線カメラで大きな室内の多数の観客を撮影してその画像中の所定閾値以上の温度の画素を人体を表す画素と判定してその総画素から総面積を求めて基準画素面積で除し、室内の総人数を計測する技術がある（特許文献１）。
上述した前者では、管理者が飼育ケージ内の多数のブロイラの生死を確認して死亡数を算出する必要があり、数百羽から数千羽も一箇所の飼育ケージ内で飼育されているブロイラの死亡数を算出することは極めて大変であり、このため、飼料要求率の低下や、腐臭の発生や、出荷時の計測により初めてその死亡率が予想した死亡率よりも大幅に多くて危険な状態（例えば鳥インフルエンザなどによる多量死の状態）をリアルタイムに監視して検出できないという問題があった。 In the former described above, it is necessary for the administrator to check the number of broilers in the breeding cage to calculate the number of deaths, and hundreds to thousands of broilers are housed in one breeding cage. Therefore, it is extremely difficult to calculate the number of deaths, and for this reason, the mortality rate is significantly higher than the expected death rate for the first time due to the decrease in feed demand rate, the occurrence of odors, and measurement at the time of shipment. There has been a problem that the state (for example, the state of mass death due to avian influenza) cannot be monitored and detected in real time.
また、上述した後者では、室内の観客数を赤外線カメラによって撮影した画像中から所定閾値温度以上の画素の総数をもとに顧客人数を計測しているため、ブロイラのように、日数ともにそのサイズが成長する（大きくなる）場合には適用できないという問題があった。また、ブロイラの日数に対応した基準となる面積を予め算出して保持してこれをもとに総人数（総ブロイラ数）を算出したのでは、現実の成長度合いが異なる場合には誤差が大きくなって実用に供し得ないという問題もあった。 Moreover, in the latter mentioned above, since the number of customers is measured based on the total number of pixels having a predetermined threshold temperature or more from the image taken by the infrared camera, the size of both the number of days as in a broiler. There is a problem that it cannot be applied when the growth (growth). In addition, if the total number of people (total number of broilers) is calculated based on the area that serves as a reference corresponding to the number of days of the broiler and is stored in advance, there will be a large error if the actual degree of growth differs. There was also a problem that it could not be put to practical use.
本発明は、これらの問題を解決するため、飼育ケージ内で飼育中の多数のブロイラを日毎などにサーモグラフで撮影した画像から所定温度以上の総面積を求めて鶏数（当初は飼育時の総数、それ以外は死亡数を除いた現生存数）に標準表面積を乗算して理論総表面積を算出して両者の比である画像上における誤差率を算出して記録すると共に、誤差率が所定値よりも小さくあるいは指定日数経過後に、日毎などに撮影した画像から求めた総面積が現日齢の標準表面積、直前に求めた総鶏数、死亡危険率から求めた現在の許容表面積よりも小さいときに危険アラームを発信するようにしている。 In order to solve these problems, the present invention obtains the total area above a predetermined temperature from images obtained by thermographing a large number of broilers being raised in a breeding cage every day, etc. Multiply the total number, otherwise the number of survivors excluding the number of deaths) by multiplying the standard surface area to calculate the theoretical total surface area and calculate and record the error rate on the image, which is the ratio of both, and the error rate is predetermined The total area obtained from images taken every day after the specified number of days has elapsed is smaller than the current allowable surface area obtained from the standard surface area of the current day, the total number of chickens obtained immediately before, and the risk of death I sometimes try to send out a danger alarm.
従って、日齢により成長して大きくなる度合いや画像上で重なってブロイラが撮影されることなどによる補正を、現状に合わせて正確に自動補正し、ブロイラの死亡率をリアルタイムに提示すると共に死亡率が危険率以上となったときに高精度に自動警報を発することが可能となる。 Therefore, corrections based on the degree of growth and growth due to age and the fact that the broiler is photographed over the image are automatically and accurately corrected according to the current situation, and the mortality rate of the broiler is presented in real time. It becomes possible to issue an automatic alarm with high accuracy when the value becomes higher than the danger rate.
本願発明は、日齢により成長して大きくなる度合いや画像上で重なってブロイラが撮影されることなどによる補正を、現状に合わせて正確に自動補正し、ブロイラの死亡率をリアルタイムに提示すると共に死亡率が危険率以上となったときに高精度に自動警報を発することが可能となる。 The present invention automatically corrects the degree to which the broiler grows and grows with age and the image of the broiler being photographed over the image, and automatically corrects it according to the current situation, and presents the mortality rate of the broiler in real time. An automatic alarm can be issued with high accuracy when the mortality rate exceeds the risk rate.
本発明は、飼育ケージ内で飼育中の多数のブロイラを日毎などにサーモグラフで撮影した画像から所定温度以上の総面積を求めて鶏数（当初は飼育時の総数、それ以外は死亡数を除いた現生存数）に標準表面積を乗算して理論総表面積を算出して両者の比である画像上における誤差率を算出して記録すると共に、誤差率が所定値よりも小さくあるいは指定日数経過後に、日毎などに撮影した画像から求めた総面積が現日齢の標準表面積、直前に求めた総鶏数、死亡危険率から求めた現在の許容表面積よりも小さいときに危険アラームを発信し、ブロイラの死亡率をリアルタイムに提示すると共に死亡率が危険率以上となったときに高精度に自動警報を発することを実現した。 The present invention obtains the total area above a predetermined temperature from images obtained by thermographing a large number of broilers kept in a breeding cage on a daily basis and the number of chickens (initially the total number at breeding, otherwise the number of deaths). The total number of survivors) is multiplied by the standard surface area to calculate the theoretical total surface area, and the error rate on the image, which is the ratio between the two, is calculated and recorded, and the error rate is less than the specified value or the specified number of days has elapsed. Later, when the total area calculated from images taken every day etc. is smaller than the current surface area calculated from the standard surface area of the current age, the total number of chickens calculated immediately before, the death risk rate, In addition to presenting the mortality rate of the broiler in real time, it realized the automatic warning with high accuracy when the mortality rate exceeded the risk rate.
図１は、本発明のシステム構成図を示す。
図１において、処理装置１は、プログラムに従い各種処理を実行するものであって、ここでは、サーモグラフィ読み込み処理手段１１、誤差率算出手段１２、理論表面積算出処理手段１３、安全領域判定処理手段１４、危険アラーム発信処理手段１５、標準表面積テーブル１６、誤差率テーブル１７、および死亡率テーブル１８などから構成されるものである。
FIG. 1 shows a system configuration diagram of the present invention.
In FIG. 1, a
サーモグラフィ読み込み処理手段１１は、サーモグラフィ装置２を制御し、鶏舎内で飼育されているブロイラの画像を撮影してその画像を読み込むものである。
The thermography reading processing means 11 controls the
誤差率算出処理手段１２は、撮影した画像から算出したブロイラの高温部の画素に対応する総実測面積と、ブロイラの総数に対する理論総表面積とから誤差率を算出するものである（図２から図４、および図５の（ｃ）の誤差率テーブル１７参照）。 The error rate calculation processing means 12 calculates an error rate from the total measured area corresponding to the pixels in the high temperature portion of the broiler calculated from the captured image and the theoretical total surface area with respect to the total number of broilers (FIG. 2 to FIG. 2). 4 and the error rate table 17 of FIG.
理論表面積算出処理手段１３は、ブロイラの総数に対応する理論表面積を算出するものである（図２から図４、および図５の（ａ）の標準表面積テーブル１６参照）。 The theoretical surface area calculation processing means 13 calculates a theoretical surface area corresponding to the total number of broilers (see the standard surface area table 16 in FIGS. 2 to 4 and FIG. 5A).
安全領域判定処理手段１４は、理論表面積および実際に撮影した画像から算出した総実測面積などをもとに安全領域か否かを判定するものである（図４参照）。 The safe area determination processing means 14 determines whether or not the area is a safe area based on the theoretical surface area and the total actually measured area calculated from the actually captured image (see FIG. 4).
危険アラーム発信処理手段１５は、安全領域判定処理手段１４により安全領域でない（危険領域）と判定されたときに、死亡率が危険領域（ブロイラの死亡率が大）にある旨のアラームを発信するものである（図４参照）。 The danger alarm transmission processing means 15 transmits an alarm indicating that the mortality rate is in the dangerous area (the mortality rate of the broiler is large) when the safety area determination processing means 14 determines that the safety area is not a safety area (risk area). (See FIG. 4).
標準表面積テーブル１６は、ブロイラの日齢毎の標準表面積を測定して登録したものである（図５の（ａ）参照）。 The standard surface area table 16 is obtained by measuring and registering the standard surface area for each age of the broiler (see FIG. 5A).
誤差率テーブル１７は、総実測面積と、理論総表面積とから算出して登録したものである（図２から図４、および図５の（ｃ）の誤差率テーブル１７参照）。 The error rate table 17 is calculated and registered from the total measured area and the theoretical total surface area (see the error rate table 17 in FIGS. 2 to 4 and FIG. 5C).
死亡率テーブル１８は、日齢毎の死亡率を登録したものである（図５の（ｂ）参照）。
サーモグラフィ装置２は、図示外の鶏舎内で飼育されているブロイラを、赤外線で撮影しその温度分布を持つ画像を生成する公知の装置である。
The mortality rate table 18 registers the mortality rate for each age (see FIG. 5B).
The
端末３は、処理装置１に接続（直接に接続、あるいはネットワークを介して接続）し、管理者４が操作して各種指示、データを入力したり、結果を表示したりなどするものである。
The
管理者４は、端末３を操作して処理装置１に各種処理を実行させるものである（図２から図４参照）。
The
次に、図２から図４のフローチャートの順番に従い、図１の構成の動作を順次詳細に説明する。 Next, the operation of the configuration of FIG. 1 will be described sequentially in detail according to the order of the flowcharts of FIGS.
図２から図４は、本発明の動作説明フローチャートを示す。
図２において、Ｓ１は、搬入日に入雛数（Ｈ０）および入雛時の日齢を入力する。これは、管理者４が端末３から、鶏舎にブロイラの雛を搬入した搬入日の雛数Ｈ０および雛の日齢を入力する。
2 to 4 are flowcharts for explaining the operation of the present invention.
In FIG. 2, S1 inputs the number of entering chicks (H0) and the age of entering dolls on the date of delivery. The
Ｓ２は、サーモグラフィ装置２よりデータを収集し、高温部の面積（総実測面積Ｓ０＝ΣＳｉ）を算出する。これは、Ｓ１で雛を鶏舎に入れて飼育を開始した時点で、サーモグラフィ装置２で当該鶏舎内の雛の画像を撮影し、撮影した画像上で、所定温度以上の高温部の画素の全てに対応する面積（雛の実測された総実測面積Ｓ０）を算出する。
In S2, data is collected from the
Ｓ３は、搬入日の日齢の標準表面積（Ｍ０）を標準表面積テーブル１６より読み込む。これは、Ｓ１で入力された搬入日の雛の日齢をもとに、後述する図５の（ａ）の標準表面積テーブル１６から当該日齢の標準表面積（Ｍ０）を取り込む。
In
Ｓ４は、入雛総数の理論総表面積（Ｐ０）を算出する。
Ｐ０＝Ｈ０×Ｍ０（総雛数×標準表面積）
これは、Ｓ３で取り込んだ搬入日の日齢の標準表面積Ｍ０に、搬入した雛の総数Ｍ０を乗算し、理論総表面積Ｐ０を算出する。
In S4, the theoretical total surface area (P0) of the total number of nestlings is calculated.
P0 = H0 × M0 (total number of chicks × standard surface area)
In this case, the theoretical total surface area P0 is calculated by multiplying the standard surface area M0 of the day of the loading date taken in S3 by the total number M0 of the loaded chicks.
Ｓ５は、搬入初日の誤差率（Ｒ０）を算出し、誤差率テーブル１７に書き込む。ここで、誤差率Ｒ０＝（総実測面積Ｓ０）／（理論総表面積Ｐ０）で当該誤差率Ｒ０を算出する（誤差率Ｒ０は、式が表すように、実際に撮影した画像上で高温部の総実測面積Ｓ０を標準的な雛の面積に雛の総数を乗算した雛の重なりなどの無い理想的な状態の理論総面積Ｐ０で除算した値（誤差率）である）。
In
Ｓ６は、死亡雛ありか判別する。これは、管理者が鶏舎を監視などして死亡雛があるか判別する。ＹＥＳの場合には、Ｓ７で死亡雛（鶏）数Ｋｎを入力し、現生存数（Ｌｎ）を算出する（Ｌ０＝Ｈ０−Ｋｎ）。そして、Ｓ８に進む。一方、Ｓ６のＮＯの場合には、そのままＳ８に進む。 In S6, it is determined whether there are dead chicks. In this case, the manager monitors the poultry house to determine whether there are dead chicks. In the case of YES, the number of dead chicks (chicken) Kn is input in S7, and the current survival number (Ln) is calculated (L0 = H0-Kn). Then, the process proceeds to S8. On the other hand, if NO in S6, the process proceeds directly to S8.
Ｓ８は、サーモグラフィ装置２よりデータを収集し、高温部の総実測面積Ｓｎを算出する（Ｓｎ＝ΣＳｉ）。これは、鶏舎のブロイラをサーモグラフィ装置２で撮影してその画像上で、高温部の総実測面積Ｓｎを算出（画像上で、所定温度以上の全ての画素に対応するブロイラに相当する総実測面積Ｓｎを算出）する。
In S8, data is collected from the
Ｓ９は、搬入ｎ日の日齢の標準表面積（Ｍｎ）を標準表面積テーブル１６から読み込む。
In
図３において、Ｓ１０は、現鶏の理論表面積の合計（理論総表面積Ｐｎ）を算出する（Ｐｎ＝Ｌｎ×Ｍｎ）。これは、現生存鶏数Ｌｎ（死亡鶏が発生したときはＳ７で当該死亡鶏数を減算した残りの現生存鶏数）に、Ｓ９で図５の（ａ）の標準表面積テーブル１６から取り出した現鶏の日齢の理論表面積Ｍｎを乗算して算出した全ての現鶏の理論総表面積Ｐｎを算出する。 In FIG. 3, S10 calculates the total theoretical surface area (theoretical total surface area Pn) of the current chicken (Pn = Ln × Mn). This is taken out from the standard surface area table 16 of FIG. 5 (a) in S9 to the current number of live chickens Ln (the number of remaining live chickens obtained by subtracting the number of dead chickens in S7 when dead chickens are generated). The theoretical total surface area Pn of all current chickens calculated by multiplying the theoretical surface area Mn of the age of the current chicken is calculated.
Ｓ１１は、搬入ｎ日の誤差率（Ｒｎ）を算出し、誤差率テーブル１７に書き込む（Ｒｎ＝Ｓｎ／Ｐｎ）。これは、Ｓ８で算出した総実測表面積Ｓｎを、Ｓ１０で算出した理論総表面積Ｐｎで除算し、誤差率Ｒｎを算出し、図５の（ｃ）の誤差率テーブル１７の該当する日齢の欄に書き込む。 In S11, the error rate (Rn) of the carry-in n day is calculated and written in the error rate table 17 (Rn = Sn / Pn). This is obtained by dividing the total measured surface area Sn calculated in S8 by the theoretical total surface area Pn calculated in S10 to calculate the error rate Rn, and the corresponding age column of the error rate table 17 in FIG. Write to.
Ｓ１２は、ｎ日＜ｍ日か判別する。ここで、ｍは誤差率が安定するまでの日数や鶏舎の詳細な確認可能な日までの日数をパラメータとして設定（入力など）する。また、ｎは雛を鶏舎に搬入した日から現在までの経過日数である。ＹＥＳの場合には、ｎ日（搬入してからの経過日数）が、誤差率が安定するまでの日数（ｍ日で、管理者が設定できるパラメータ）よりも小さく、未だ、誤差率が変動して画像上の総実測面積で正確に雛の数を推測するのに十分でないと判明したので、図２のＳ６に戻り、日毎に誤差率を算出して図５の（ｃ）の誤差率テーブル１７の該当日齢の欄に書き込むことを繰り返す。一方、ＮＯの場合には、ｎ日（搬入してからの経過日数）が、誤差率が安定するまでの日数（ｍ日で、管理者が設定できるパラメータ）を経過し、誤差率の変動が所定閾値以下となり、画像上の総実測面積で正確に雛の数を推測するのに十分であると判明したので、Ｓ１３に進む。 S12 determines whether n days <m days. Here, m is set (input, etc.) as parameters such as the number of days until the error rate is stabilized and the number of days until the day when the poultry house can be confirmed in detail. N is the number of days that have passed since the day when the chicks were carried into the poultry house. In the case of YES, n days (the number of days since the delivery) is smaller than the number of days until the error rate stabilizes (a parameter that can be set by the administrator in m days), and the error rate still fluctuates. Thus, since it has been found that the total actual area on the image is not sufficient to accurately estimate the number of chicks, the process returns to S6 in FIG. 2 to calculate the error rate for each day, and the error rate table in FIG. Repeat writing in the column for 17 days of age. On the other hand, in the case of NO, n days (the number of days since the delivery) has passed the number of days until the error rate is stabilized (a parameter that can be set by the administrator in m days), and the error rate fluctuates. Since it has become less than the predetermined threshold and is found to be sufficient to accurately estimate the number of chicks with the total measured area on the image, the process proceeds to S13.
Ｓ１３は、死亡鶏の確認可能か判別する。ＹＥＳの場合には、Ｓ１４に進む。ＮＯの場合には、図４のＳ２１に進む。 S13 discriminate | determines whether the dead chicken can be confirmed. If YES, the process proceeds to S14. In the case of NO, the process proceeds to S21 in FIG.
Ｓ１４は、Ｓ１３のＹＥＳで死亡鶏の確認可能と判明したので、死亡鶏ありか判別する。ＹＥＳの場合には、Ｓ１５で死亡雛（鶏）数Ｋｘを入力し、現生存数（Ｌｘ）を算出する（Ｌｘ＝Ｌｎ−Ｋｘ、Ｌｎは既述したＳ７で算出したなどの現生存数である）。そして、Ｓ１６に進む。一方、Ｓ１４のＮＯの場合には、そのままＳ１６に進む。 In S14, it is determined that the dead chicken can be confirmed by YES in S13, so it is determined whether there is a dead chicken. In the case of YES, the number Kx of dead chicks (chicken) is input in S15, and the current survival number (Lx) is calculated (Lx = Ln−Kx, Ln is the current survival number such as that calculated in S7 described above) is there). Then, the process proceeds to S16. On the other hand, in the case of NO in S14, the process proceeds to S16 as it is.
Ｓ１６は、サーモグラフィ装置２よりデータを収集し、高温部の総実測面積Ｓｘを算出する（Ｓｘ＝ΣＳｉ）。これは、鶏舎のブロイラをサーモグラフィ装置２で撮影してその画像上で、高温部の総実測面積Ｓｘを算出（画像上で、所定温度以上の全ての画素に対応するブロイラに相当する総実測面積Ｓｘを算出）する。
In S16, data is collected from the
Ｓ１７は、現在（ｘ）の日齢の標準表面積（Ｍｘ）を標準表面積テーブル１６から読み込む。
In
Ｓ１８は、現鶏の理論表面積の合計（理論総表面積Ｐｘ）を算出する（Ｐｘ＝Ｌｘ×Ｍｘ）。これは、現生存鶏数Ｌｘ（死亡鶏が発生したときはＳ１５で当該死亡鶏数を減算した残りの現生存鶏数）に、Ｓ１７で図５の（ａ）の標準表面積テーブル１６から取り出した現鶏の日齢の理論表面積Ｍｘを乗算して算出した全ての現鶏の理論総表面積Ｐｘを算出する。 In S18, the total theoretical surface area (theoretical total surface area Px) of the current chicken is calculated (Px = Lx × Mx). This is taken out from the standard surface area table 16 of FIG. 5 (a) in S17 to the current number of live chickens Lx (the number of remaining live chickens obtained by subtracting the number of dead chickens in S15 when dead chickens are generated). The theoretical total surface area Px of all current chickens calculated by multiplying the theoretical surface area Mx of the age of the current chicken is calculated.
図４のＳ１９は、現在の誤差率（Ｒｘ）を算出し、誤差率テーブル１７に書き込む（Ｒｘ＝Ｓｘ／Ｐｘ）。これは、Ｓ１６で算出した総実測表面積Ｓｘを、Ｓ１８で算出した理論総表面積Ｐｘで除算し、誤差率Ｒｘを算出し、図５の（ｃ）の誤差率テーブル１７の該当する日齢の欄に書き込む。 In S19 of FIG. 4, the current error rate (Rx) is calculated and written in the error rate table 17 (Rx = Sx / Px). This is obtained by dividing the total measured surface area Sx calculated in S16 by the theoretical total surface area Px calculated in S18 to calculate the error rate Rx, and the corresponding age column of the error rate table 17 in FIG. Write to.
Ｓ２０は、誤差率関数を作成（算出）する（Ｒｘ＝ｆ（ｘ））。ここで、誤差率関数ｆ（ｘ）は、図５の（ｃ）の誤差率テーブル１７に書き込まれた日齢（ｘ日）をＸ軸、そのときの誤差率ＲｘをＹ軸としてプロットしこれらを曲線で結んだ（あるいは近似した）誤差率関数を作成する。 In S20, an error rate function is created (calculated) (Rx = f (x)). Here, the error rate function f (x) is plotted with the age (x days) written in the error rate table 17 of FIG. 5C as the X axis and the error rate Rx at that time as the Y axis. Create an error rate function that connects (or approximates) with a curve.
Ｓ２１は、サーモグラフィ装置２よりデータを収集し、高温部の総実測面積ＳＧを算出する（ＳＧ＝ΣＳｉ）。これは、鶏舎のブロイラをサーモグラフィ装置２で撮影してその画像上で、高温部の総実測面積ＳＧを算出（画像上で、所定温度以上の全ての画素に対応するブロイラに相当する総実測面積ＳＧを算出）する。
In S21, data is collected from the
Ｓ２２は、現在（Ｇ）の日齢の標準表面積（ＭＧ）を標準表面積テーブル１６から読み込む。 In S <b> 22, the standard surface area (MG) of the current age (G) is read from the standard surface area table 16.
Ｓ２３は、現鶏の理論表面積の合計（理論総表面積ＰＧ）を算出する（ＰＧ＝Ｌｘ×ＭＧ）。これは、現生存鶏数Ｌｘ（死亡鶏が発生したときはＳ１５などで当該死亡鶏数を減算した残りの現生存鶏数）に、Ｓ２２で図５の（ａ）の標準表面積テーブル１６から取り出した現鶏の日齢の理論表面積ＭＧを乗算して算出した全ての現鶏の理論総表面積ＰＧを算出する。 In S23, the total theoretical surface area (theoretical total surface area PG) of the current chicken is calculated (PG = Lx × MG). This is taken out from the standard surface area table 16 of FIG. 5 (a) in S22 to the current number of live chickens Lx (the number of remaining live chickens obtained by subtracting the number of dead chickens in S15 when a dead chicken is generated). The theoretical total surface area PG of all the current chickens calculated by multiplying the theoretical surface area MG of the age of the current chicken is calculated.
Ｓ２４は、現日齢での差分の危険死亡率（ＤＲ）より安全限界の理論表面積（Ｆｓ）を算出する（Ｆｓ＝ＰＧ×（１−ＤＲ）。ここで、安全限界の理論表面積Ｆｓは、Ｓ２３で算出した現鶏の理論総表面積ＰＧに、１から危険死亡率ＤＲを減算してその減算した値を乗算して算出する。死亡率は、ブロイラの標準死亡率（図５の（ｂ）の標準死亡率テーブル１８）を使用するが、社会情勢の変化（特定の疾病が流行中）によりパラメータ（α）を設定可能にする。ＤＲ＝死亡率×α（α＜＝１）とする。 S24 calculates the theoretical surface area (Fs) of the safety limit from the difference risk mortality rate (DR) at the current day age (Fs = PG × (1−DR), where the theoretical surface area Fs of the safety limit is The theoretical total surface area PG of the current chicken calculated in S23 is calculated by subtracting the risk mortality DR from 1 and multiplying the result by subtracting the risk mortality DR. The standard mortality table 18) is used, but the parameter (α) can be set by changing the social situation (a specific disease is prevalent): DR = mortality × α (α <= 1).
Ｓ２５は、理論安全限界面積（Ｆｓ）を補正（Ｒ（Ｇ））し、現在の許容表面積（Ｆ（ｋ））を算出する（Ｆ（ｋ）＝ＰＧ×Ｒ（Ｇ））。ここで、理論安全限界面積Ｆｓは、Ｓ２３で算出した現鶏の理論総表面積ＰＧに、補正係数Ｒ（Ｇ）を乗算して算出する。補正係数Ｒ（Ｇ）は、実験、経験などにより理論的なな安全限界面積として妥当となるように予め設定する。 In S25, the theoretical safety limit area (Fs) is corrected (R (G)), and the current allowable surface area (F (k)) is calculated (F (k) = PG × R (G)). Here, the theoretical safety limit area Fs is calculated by multiplying the theoretical total surface area PG of the current chicken calculated in S23 by the correction coefficient R (G). The correction coefficient R (G) is set in advance so as to be appropriate as a theoretical safety limit area based on experiments and experiences.
Ｓ２６は、Ｆ（ｋ）＜ＳＧか判別する。これは、Ｓ２５で算出した現鶏の理論安全限界面積Ｆｓが、Ｓ２１で算出した現鶏の画像上で算出した総実測面積ＳＧよりも小さくて危険死亡率以下の鶏が死亡している（鶏が多量死亡でない）か判別する。ＹＥＳの場合には、現鶏の死亡率が安全限界以下と判明したので、図３のＳ１３に戻り繰り返す。一方，ＮＯの場合には、現鶏の死亡率が安全限界以上と判明したので、Ｓ２７で危険アラーム（現鶏が多量死した旨のアラーム）を発信する。 S26 determines whether F (k) <SG. This is because the theoretical safety limit area Fs of the current chicken calculated in S25 is smaller than the total actually measured area SG calculated on the image of the current chicken calculated in S21, and the chickens having a risk mortality or less have died (chicken Is not a large death). In the case of YES, since the mortality rate of the current chicken is found to be below the safety limit, the process returns to S13 in FIG. 3 and is repeated. On the other hand, in the case of NO, since the mortality rate of the current chicken is found to be above the safety limit, a danger alarm (alarm indicating that the current chicken has died in large quantities) is transmitted in S27.
以上によって、雛を鶏舎に搬入した日から日毎（日齢毎）にサーモグラフィ装置２で撮影した画像上で現鶏の総実測面積Ｓを算出および現日齢の理論総表面積Ｐを算出してその誤差率Ｒを求めて誤差率テーブル１７に書き込み、誤差率が安定した日以降（あるいは指定された日以降）に測定を開始し、日齢毎に測定した現鶏の総実測面積Ｓとその理論総表面積Ｐ、更に、死亡率などをもとに当該日齢毎に安全限界を予測してそれを越える現鶏の総実測面積Ｓとなったときに多量死亡と判定（Ｓ２６のＮＯ）して、危険アラームを自動発信することが可能となる。
As described above, the total actual measurement area S of the current chicken is calculated and the theoretical total surface area P of the current day is calculated on the image taken by the
図５は、本発明のテーブル例を示す。
図５の（ａ）は、標準表面積テーブルの例を示す。図示の標準面積テーブル１６は、ブロイラ（肉鶏）の１羽の標準表面積を日齢毎に登録したものであって、ここでは、図示の下記の情報を対応づけて登録したものである。
FIG. 5 shows an example table of the present invention.
FIG. 5A shows an example of a standard surface area table. The illustrated standard area table 16 is a table in which the standard surface area of one broiler (meat chicken) is registered for each age, and here, the following information shown in the figure is associated and registered.
・日齢：
・面積：
ここで、日齢は鶏（ブロイラ）が生まれてからの日数である。面積は、ブロイラの１羽当りの画像として撮影したときのその標準的な面積である。
-Age:
·area:
Here, the age is the number of days since the chicken (broiler) was born. The area is a standard area taken as an image per broiler.
以上のように、鶏の日齢に対応づけて１羽当りの標準的な面積を予め登録しておくことにより、日々成長して大きくなる鶏について日齢毎に鶏数の全体の理論総面積を正確に算出などすることが可能となる。 As described above, by registering the standard area per bird in advance in association with the age of the chicken, the theoretical total area of the total number of chickens for each age for chickens that grow and grow daily. Can be calculated accurately.
図５の（ｂ）は、死亡率テーブルの例を示す。図示の死亡率テーブル１８は、ブロイラ（肉鶏）の日齢に従う死亡率を登録したものであって、ここでは、図示の下記の情報を対応づけて登録したものである。 FIG. 5B shows an example of a mortality table. The illustrated mortality table 18 registers the mortality according to the age of the broiler (meat chicken), and here, the following information shown in FIG.
・日齢：
・死亡率（％）：
ここで、日齢は鶏（ブロイラ）が生まれてからの日数である。死亡率は日齢に従い死亡する鶏の割合（％）を予め実測して標準的な値を登録したものである。
-Age:
-Mortality rate (%):
Here, the age is the number of days since the chicken (broiler) was born. The mortality rate is obtained by registering a standard value by actually measuring the proportion (%) of chickens that die according to age.
以上のように、鶏の日齢に対応づけて死亡率を予め登録しておくことにより、日々成長して大きくなると共に鶏舎で飼育中に死亡する鶏の死亡率を考慮して日齢毎の現生存率を正確に算出などすることが可能となる。 As described above, by registering the mortality rate in advance in association with the age of the chicken, the mortality rate of the chicken that grows and grows daily while taking into account the mortality rate of the chickens that die during breeding in the poultry house It becomes possible to accurately calculate the current survival rate.
図５の（ｃ）は、誤差率テーブルの例を示す。図示の誤差率テーブル１７は、鶏舎で飼育中のブロイラ（肉鶏）群をサーモグラフィ装置２で撮影した画像から求めた総実測面積と、理論総表面積との誤差率を日齢毎に測定して書き込んだものであって、ここでは、図示の下記の情報を対応づけて登録したものである。
FIG. 5C shows an example of the error rate table. The illustrated error rate table 17 measures the error rate between the total measured area obtained from the image obtained by photographing the broiler (meat chicken) group being raised in the poultry house with the
・搬入日：
・誤差率（％）：
ここで、搬入日は鶏（ブロイラ）を鶏舎に搬入して飼育を開始した日である。搬入日に、搬入日の雛の日齢を加算することで、当該鶏の日齢を算出できる。誤差率は、既述したように、総実測面積と理論総面積との誤差率である。
·Installation date:
·Error rate(%):
Here, the carry-in date is the day when the chicken (broiler) is carried into the poultry house and breeding is started. The age of the chicken can be calculated by adding the age of the chicks on the date of delivery. As described above, the error rate is an error rate between the total measured area and the theoretical total area.
以上のように、鶏の搬入日（日齢）に対応づけて誤差率を書き込んでチェックすることにより、総実測面積と理論総面積との誤差が所定範囲に収まり安定した場合や安定すると予測される日齢を指定し、当該安定した日以降あるいは指定された日以降に、総実測面積と理論総面積をもとに日々成長して大きくなる鶏の多量死の判定の正確性を担保することが可能となる。 As described above, it is predicted that the error between the total measured area and the theoretical total area falls within a predetermined range and is stable or stable by writing and checking the error rate in association with the day of chicken delivery (age). To ensure the accuracy of the determination of the mass death of chickens that grow and grow daily based on the total measured area and the theoretical total area after the stable day or after the specified date Is possible.
（付記１）
飼育中の多数の鶏をサーモグラフで撮影した画像から鶏の死亡率を自動判定する死亡率自動判定方法において、
サーモグラフで多数の鶏を撮影して画像を取得するステップと、
前記取得した画像をもとに所定温度以上の画素に対応する鶏の総面積を算出するステップと、
テーブルから取り込んだ鶏の現日齢の標準表面積、直前に算出された総鶏数、および危険死亡率から現日齢の許容表面積を算出するステップと、
前記算出した総面積が前記算出した許容表面積よりも小さいときに死亡率が大きいと判定するステップと
を有するサーモグラフによる死亡率自動判定方法。
(Appendix 1)
In the automatic mortality determination method that automatically determines the mortality of chickens from images taken by thermographs of many chickens being raised,
Taking a large number of chickens with a thermograph and acquiring images;
Calculating the total area of the chicken corresponding to pixels above a predetermined temperature based on the acquired image;
Calculating an acceptable surface area for the current age from the standard surface area for the current age of chickens taken from the table, the total number of chickens calculated immediately before, and the risk mortality rate;
A method for automatically determining a mortality rate by a thermograph, comprising: determining that the mortality rate is large when the calculated total area is smaller than the calculated allowable surface area.
（付記２）
前記標準表面積および危険死亡率は、鶏が成長する日齢により変化させることを特徴とする付記１記載のサーモグラフによる死亡率自動判定方法。
(Appendix 2)
The method for automatically determining a mortality rate by a thermograph according to
（付記３）
前記判定は、前記算出した総面積と、前記取り込んだ現日齢の標準表面積に総鶏数を乗算した面積との比である誤差率が所定範囲内に収まったとき以降に開始あるいは指定された日齢以降に開始することを特徴とする付記１あるいは付記２記載のサーモグラフによる死亡率自動判定方法。
(Appendix 3)
The determination is started or specified after an error rate that is a ratio of the calculated total area and the area obtained by multiplying the standard surface area of the captured current age by the total number of chickens falls within a predetermined range. The method for automatically determining a mortality rate according to the thermograph according to
（付記４）
飼育中の多数の鶏をサーモグラフで撮影した画像から鶏の死亡率を自動判定する死亡率自動判定装置において、
サーモグラフで多数の鶏を撮影して画像を取得する手段と、
前記取得した画像をもとに所定温度以上の画素に対応する鶏の総面積を算出する手段と、
テーブルから取り込んだ鶏の現日齢の標準表面積、直前に算出された総鶏数、および危険死亡率から現日齢の許容表面積を算出する手段と、
前記算出した総面積が前記算出した許容表面積よりも小さいときに死亡率が大きいと判定する手段と
を備えたことを特徴とするサーモグラフによる死亡率自動判定装置。
(Appendix 4)
In the automatic mortality determination device that automatically determines the mortality rate of chickens from images taken with a thermograph of a large number of chickens being raised,
Means for taking images of a large number of chickens with a thermograph and acquiring images;
Means for calculating the total area of the chicken corresponding to pixels of a predetermined temperature or higher based on the acquired image;
Means for calculating the allowable surface area of the current day from the standard surface area of the current day of the chicken taken from the table, the total number of chickens calculated immediately before, and the risk mortality rate;
An automatic mortality determination apparatus using a thermograph, comprising: means for determining that the mortality rate is large when the calculated total area is smaller than the calculated allowable surface area.
本発明は、日齢により成長して大きくなる度合いや画像上で重なってブロイラが撮影されることなどによる補正を、現状に合わせて正確に自動補正し、ブロイラの死亡率をリアルタイムに提示すると共に死亡率が危険率以上となったときに高精度に自動警報を発することが可能となる。 The present invention automatically corrects the degree of growth and enlargement with age and the fact that the broiler is photographed over the image and accurately corrects it according to the current situation, and presents the mortality rate of the broiler in real time. An automatic alarm can be issued with high accuracy when the mortality rate exceeds the risk rate.
１：処理装置
１１：サーモグラフィ読み込み処理手段
１２：誤差率算出処理手段
１３：理論表面積算出処理手段
１４：安全領域判定処理手段
１５：危険アラーム発信処理手段
１６：標準表面積テーブル
１７：誤差率テーブル
１８：死亡率テーブル
２：サーモグラフィ装置
３：端末
４：管理者
1: Processing device 11: Thermography reading processing means 12: Error rate calculation processing means 13: Theoretical surface area calculation processing means 14: Safety area determination processing means 15: Danger alarm transmission processing means 16: Standard surface area table 17: Error rate table 18: Mortality table 2: Thermography device 3: Terminal 4: Administrator
Claims (3)
サーモグラフで複数の鶏を撮影した画像を取得するステップと、
前記取得した画像をもとに所定温度以上の画素に対応する面積を算出するステップと、
鶏の日齢と標準面積とを対応付けて格納する標準表面積テーブルから取り込んだ、前記画像の撮影時の日齢に対応する標準表面積と、前記画像の撮影時の鶏数と、日齢と対応付けられた死亡率から算出した、前記画像の撮影時の日齢に対応する危険死亡率とから、前記画像の撮影時の日齢の許容表面積（Ｆｋ）を算出するステップと、
前記算出した面積が前記算出した許容表面積よりも小さいときに鶏の死亡率が大きいと判定するステップと
を有する死亡率自動判定方法。 In the automatic mortality determination method for determining the mortality rate of chickens executed by a computer,
Acquiring a picture of multiple chickens taken with a thermograph;
Calculating an area corresponding to a pixel having a predetermined temperature or higher based on the acquired image;
Taken from a standard surface area table for storing in association with each other ages and the standard area of the chicken, and the standard surface area corresponding to the ages at the time of photographing of the image, and chicken number at the time of photographing of the image, corresponding to the day-old was calculated from the Tagged mortality, and calculating from the danger of death rate corresponding to the ages at the time of photographing of the image, the allowable surface area of the ages at the time of photographing of the image (Fk),
And a step of determining that the mortality rate of the chicken is large when the calculated area is smaller than the calculated allowable surface area.
サーモグラフで複数の鶏を撮影した画像を取得する手段と、
前記取得した画像をもとに所定温度以上の画素に対応する面積を算出する手段と、
鶏の日齢と標準面積とを対応付けて格納する標準表面積テーブルと、
前記標準表面積テーブルから取り込んだ、前記画像の撮影時の日齢に対応する標準表面積と、前記画像の撮影時の鶏数と、日齢と対応付けられた死亡率から算出した、前記画像の撮影時の日齢に対応する危険死亡率とから、前記画像の撮影時の日齢の許容表面積を算出する手段と、
前記算出した面積が前記算出した許容表面積よりも小さいときに鶏の死亡率が大きいと判定する手段と
を有することを特徴とする死亡率自動判定装置。 In the automatic mortality determination device that determines the mortality rate of chickens,
Means for acquiring images of a plurality of chickens taken with a thermograph;
Means for calculating an area corresponding to a pixel having a predetermined temperature or higher based on the acquired image;
A standard surface area table for storing the age of chicken and the standard area in association with each other;
It is taken from the standard surface area table, and the standard surface area corresponding to the ages at the time of photographing of the image, was calculated from the chicken numbers and, associated with the day-old mortality at the time of photographing of the image, capturing of the image and a risk mortality corresponding to age day time, means for calculating an allowable surface area of the ages at the time of photographing of the image,
Mortality automatic determination apparatus characterized by the area that the calculated and a determining means and chicken mortality is greater when smaller than the allowable surface area the calculated.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2004236134A JP4482398B2 (en) | 2004-08-13 | 2004-08-13 | Automatic mortality determination method and mortality automatic determination device by thermograph |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2004236134A JP4482398B2 (en) | 2004-08-13 | 2004-08-13 | Automatic mortality determination method and mortality automatic determination device by thermograph |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2006050989A JP2006050989A (en) | 2006-02-23 |
| JP4482398B2 true JP4482398B2 (en) | 2010-06-16 |
Family
ID=36028850
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2004236134A Active JP4482398B2 (en) | 2004-08-13 | 2004-08-13 | Automatic mortality determination method and mortality automatic determination device by thermograph |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4482398B2 (en) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9961883B2 (en) | 2014-06-16 | 2018-05-08 | Allan Schaefer | Rapid and automatic determination of metabolic efficiency in livestock |
| US9955672B2 (en) | 2014-06-16 | 2018-05-01 | Nigel Cook | Infrared thermography and behaviour information for identification of biologically important states in animals |
| CN104621004B (en) * | 2014-12-22 | 2017-01-25 | 中国农业大学 | Living thing and environment monitoring intelligent robot |
| JP6183734B1 (en) * | 2016-12-15 | 2017-08-23 | 山田 泰弘 | Method and structure for preventing flight of harmful birds such as starlings and sparrows |
| WO2019146416A1 (en) | 2018-01-26 | 2019-08-01 | パナソニックＩｐマネジメント株式会社 | Poultry raising system, poultry raising method, and program |
| JP6712660B2 (en) * | 2018-05-24 | 2020-06-24 | 大豊産業株式会社 | Dead chicken detection system, live chicken death program, and live chicken death device |
| WO2020110616A1 (en) * | 2018-11-29 | 2020-06-04 | パナソニックＩｐマネジメント株式会社 | Poultry raising system, poultry raising method, and program |
| CN113273178A (en) | 2019-01-30 | 2021-08-17 | 松下知识产权经营株式会社 | Animal house monitoring method and animal house monitoring system |
| WO2022070128A1 (en) * | 2020-09-30 | 2022-04-07 | Pal Singh Mahender | System and method for automatic mortality detection and management |
| WO2023195063A1 (en) * | 2022-04-05 | 2023-10-12 | 株式会社ハイテム | Dead fowl detection method |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0198562U (en) * | 1987-12-17 | 1989-06-30 | ||
| JPH05268853A (en) * | 1992-03-27 | 1993-10-19 | Koshin Seiki:Kk | Monitoring device for barn |
| JPH11219437A (en) * | 1998-01-30 | 1999-08-10 | Yamatake Corp | Number of people counting system |
| JPH11276011A (en) * | 1998-03-31 | 1999-10-12 | Kubo Systems:Kk | Carcass detection in henhouse and device therefor |
| JP2002153162A (en) * | 2000-11-22 | 2002-05-28 | Fuji Electric Co Ltd | Birth sign detector |
-
2004
- 2004-08-13 JP JP2004236134A patent/JP4482398B2/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0198562U (en) * | 1987-12-17 | 1989-06-30 | ||
| JPH05268853A (en) * | 1992-03-27 | 1993-10-19 | Koshin Seiki:Kk | Monitoring device for barn |
| JPH11219437A (en) * | 1998-01-30 | 1999-08-10 | Yamatake Corp | Number of people counting system |
| JPH11276011A (en) * | 1998-03-31 | 1999-10-12 | Kubo Systems:Kk | Carcass detection in henhouse and device therefor |
| JP2002153162A (en) * | 2000-11-22 | 2002-05-28 | Fuji Electric Co Ltd | Birth sign detector |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2006050989A (en) | 2006-02-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Fernández et al. | Real-time monitoring of broiler flock's welfare status using camera-based technology | |
| Maselyne et al. | Measuring the drinking behaviour of individual pigs housed in group using radio frequency identification (RFID) | |
| JP4482398B2 (en) | Automatic mortality determination method and mortality automatic determination device by thermograph | |
| Oczak et al. | Classification of nest-building behaviour in non-crated farrowing sows on the basis of accelerometer data | |
| CN113287140A (en) | Insect breeding control device and method | |
| Van der Zande et al. | Individual detection and tracking of group housed pigs in their home pen using computer vision | |
| Mei et al. | Automated radial 8-arm maze: A voluntary and stress-free behavior test to assess spatial learning and memory in mice | |
| Adriaens et al. | Milk losses and dynamics during perturbations in dairy cows differ with parity and lactation stage | |
| De Montis et al. | Analysis of poultry eating and drinking behavior by software eYeNamic | |
| Romanini et al. | Monitoring the hatch time of individual chicken embryos | |
| KR20220137137A (en) | Livestock information management system, livestock information management server, livestock information management method, and livestock information management program stored in a storage medium | |
| Papadakis et al. | Sub-second analysis of fish behavior using a novel computer-vision system | |
| Moura et al. | Noise analysis to evaluate chick thermal comfort | |
| Hartmann et al. | Detection of the time of foaling by accelerometer technique in horses (Equus caballus)—a pilot study | |
| JPWO2020158307A1 (en) | Barn monitoring method and barn monitoring system | |
| Vafidis et al. | Supplemental food alters nest defence and incubation behaviour of an open‐nesting wetland songbird | |
| Larsen et al. | Pen fouling in finisher pigs: Changes in the lying pattern and pen temperature prior to fouling | |
| Schomburg et al. | Elevated platforms with integrated weighing beams allow automatic monitoring of usage and activity in broiler chickens | |
| JP7129675B2 (en) | Poultry farming system, poultry farming method and program | |
| CN114375158A (en) | Sow estrus judgment device, sow estrus judgment method, and sow estrus judgment program | |
| Galeana et al. | Mother-young spatial association and its relation with proximity to a fence separating ewes and lambs during enforced weaning in hair sheep (Ovis aries) | |
| US9700257B1 (en) | Method of measuring gut transit time in animals | |
| CN116963595A (en) | Pig raising support device, pig raising support method, and pig raising support program | |
| JP2023033103A (en) | Information processing system, information processing method and program | |
| WO2022181132A1 (en) | Body weight estimation system and body weight estimation method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20070720 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20090225 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20090825 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20091019 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20100316 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20100319 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130326 Year of fee payment: 3 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 4482398 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140326 Year of fee payment: 4 |