JP2013191107A - Method for estimating harvesting period and program - Google Patents
Method for estimating harvesting period and program Download PDFInfo
- Publication number
- JP2013191107A JP2013191107A JP2012058011A JP2012058011A JP2013191107A JP 2013191107 A JP2013191107 A JP 2013191107A JP 2012058011 A JP2012058011 A JP 2012058011A JP 2012058011 A JP2012058011 A JP 2012058011A JP 2013191107 A JP2013191107 A JP 2013191107A
- Authority
- JP
- Japan
- Prior art keywords
- temperature
- day
- effective
- average
- date
- 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.)
- Granted
Links
- 238000003306 harvesting Methods 0.000 title claims abstract description 82
- 238000000034 method Methods 0.000 title claims abstract description 53
- 230000008569 process Effects 0.000 claims description 14
- 238000005259 measurement Methods 0.000 claims 2
- 240000007124 Brassica oleracea Species 0.000 description 11
- 235000003899 Brassica oleracea var acephala Nutrition 0.000 description 11
- 235000011301 Brassica oleracea var capitata Nutrition 0.000 description 11
- 235000001169 Brassica oleracea var oleracea Nutrition 0.000 description 11
- 230000010365 information processing Effects 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 230000007704 transition Effects 0.000 description 5
- 230000006870 function Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000009331 sowing Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
Description
本発明は、収穫時期予測方法及びプログラムに関する。 The present invention relates to a harvest time prediction method and program.
農作物を栽培する農家や、農作物の出荷先においては、農作物がいつ頃収穫可能であるかを予測することができると、その後の作業に関して効率的な計画の立案等を行うことができる。したがって、農作物の収穫時期を予測可能な技術が求められている。 If a farmer who grows crops or a shipping destination of crops can predict when the crops can be harvested, an efficient plan for the subsequent work can be made. Therefore, there is a need for a technique that can predict the crop harvest time.
農作物の収穫時期を知るためのパラメータとして、有効積算気温が知られている。有効積算気温とは、各日の日平均気温のうち、農作物の生育に寄与する温度分(以下、「有効気温」という。)の積算値である。有効積算気温が、当該農作物の収穫時期に対応する値に到達すると、当該品種の農作物は収穫可能であると判定される。 The effective accumulated temperature is known as a parameter for knowing the crop harvest time. The effective integrated temperature is an integrated value of a temperature component (hereinafter referred to as “effective temperature”) that contributes to the growth of crops among the daily average temperatures of each day. When the effective integrated temperature reaches a value corresponding to the harvest time of the crop, it is determined that the crop of the variety can be harvested.
斯かる有効積算気温を用いて農作物の収穫時期を予測することが考えられる。収穫時期を予測する場合、予測時以降の未来の日平均気温は、未知である。そこで、未来の各日の日平均気温に関しては、過去の同月同日の日平均気温を利用して有効積算気温を算出することが考えられる。 It is conceivable to predict the harvest time of crops using such effective integrated temperature. When the harvest time is predicted, the future daily average temperature after the prediction time is unknown. Thus, regarding the daily average temperature of each future day, it is conceivable to calculate the effective integrated temperature using the daily average temperature of the same month in the past.
しかしながら、発明者による実験の結果、このような方法では、収穫時期を精度良く予測できないことが確認された。 However, as a result of experiments by the inventor, it has been confirmed that the harvesting time cannot be accurately predicted by such a method.
そこで、一側面では、農作物の収穫時期の予測精度の向上の可能性を高めることを目的とする。 Therefore, in one aspect, the object is to increase the possibility of improving the prediction accuracy of the crop harvest time.
一つの案では、収穫時期予測方法は、過去の複数年の各日の平均気温のうち収穫時期の予測対象の農作物の生育限界温度を超える分の気温である有効気温に関して、同月同日ごとの平均値を記憶する第一の記憶部を用いて、前記農作物の栽培開始日から予測実施日前までの各日の平均気温の実測値に基づく有効積算気温に対して、前記平均値を積算することにより、前記農作物の収穫時期の予測値を算出する処理をコンピュータが実行する。 In one proposal, the harvest time prediction method uses the average temperature for the same day in the same month for the effective temperature, which is the temperature that exceeds the growth limit temperature of the crop for which the harvest time is predicted, of the average temperature for each day of the past multiple years. By integrating the average value with respect to the effective integrated temperature based on the measured value of the average temperature of each day from the cultivation start date of the crop before the prediction implementation date, using the first storage unit that stores the value The computer executes a process for calculating a predicted value of the crop harvest time.
一態様によれば、農作物の収穫時期の予測精度の向上の可能性を高めることができる。 According to one aspect, it is possible to increase the possibility of improving the prediction accuracy of the crop harvest time.
以下、図面に基づいて本発明の実施の形態を説明する。本実施の形態では、農作物の具体例として品種Jのキャベツを採用する。品種Jのキャベツを、以下「キャベツJ」という。 Hereinafter, embodiments of the present invention will be described with reference to the drawings. In the present embodiment, cabbage of variety J is adopted as a specific example of the crop. The cabbage of variety J is hereinafter referred to as “cabbage J”.
本実施の形態では、有効積算気温という考え方を用いて、キャベツJの収穫予測日が算出される。キャベツJについては、栽培開始日の一例である定植日以降の栽培期間における有効積算気温が1250℃・日を超えた日の翌日が収穫可能日と仮定する。また、栽培期間における有効積算気温が1400℃・日を超えた日の翌日が収穫限界日であると仮定する。収穫可能日とは、収穫に適した期間の開始日をいう。収穫限界日とは、収穫に適した期間の終了日をいう。また、有効積算気温とは、各日の生育に寄与する温度分の積算値をいう。例えば、キャベツJでは、日の平均気温が5℃を超えたときにその温度分の生長があるとする。ある日の平均気温が12℃の場合には12−5=7℃分が生長に寄与する気温として有効積算気温に加算される。5℃以下の平均気温は、有効積算気温に加算されない。作物の生育に寄与する温度分を、有効気温という。また、有効気温に対する閾値となる温度を、生育限界温度という。上記の例では、5℃が、生育限界温度となる。したがって、日平均気温から5℃が差し引かれた値が有効気温となる。なお、「℃・日」は、日平均気温での有効積算気温の単位である。 In the present embodiment, the predicted harvest date of cabbage J is calculated using the concept of effective integrated temperature. As for cabbage J, it is assumed that the day after the day when the effective integrated temperature in the cultivation period after the planting date, which is an example of the cultivation start date, exceeds 1250 ° C. · day is the harvestable day. Further, it is assumed that the next day after the effective integrated temperature during the cultivation period exceeds 1400 ° C. · day is the harvest limit date. The harvestable date refers to the start date of a period suitable for harvesting. The harvest limit date is the end date of a period suitable for harvesting. Moreover, effective integrated temperature means the integrated value of the temperature part which contributes to growth of each day. For example, in cabbage J, when the average daily temperature exceeds 5 ° C., it is assumed that there is growth for that temperature. When the average temperature on a certain day is 12 ° C., 12−5 = 7 ° C. is added to the effective integrated temperature as the temperature contributing to growth. Average temperatures below 5 ° C are not added to the effective integrated temperature. The temperature that contributes to the growth of crops is called effective temperature. Moreover, the temperature which becomes a threshold value with respect to effective temperature is called growth limit temperature. In the above example, 5 ° C. is the growth limit temperature. Therefore, the value obtained by subtracting 5 ° C. from the daily average temperature is the effective temperature. Note that “° C./day” is a unit of effective integrated temperature at the daily average temperature.
収穫前の或る日において、収穫可能日又は収穫限界日等の収穫時期を予測する場合、当日から、収穫可能日又は収穫限界日に対応した有効積算気温に到達するまでの各日の平均気温は未知である。そこで、未知の平均気温には、過去の年の実績値を代入して、未来の有効積算気温を算出することが考えられる。なお、以下の説明において、収穫可能日から収穫限界日までの期間を「収穫期間」という。 When predicting harvest time such as possible harvest date or harvest limit date on a certain day before harvest, the average temperature of each day from the current day until reaching the effective integrated temperature corresponding to the harvest possible date or harvest limit date Is unknown. Thus, it is conceivable to calculate the future effective integrated temperature by substituting the past average value for the unknown average temperature. In the following description, the period from the harvestable date to the harvest limit date is referred to as “harvest period”.
図1は、有効積算気温の算出例を示す図である。図1では、2011年9月12日に定植を行ったキャベツJについて、同年11月12日に収穫期間の予測が実施される例が示されている。 FIG. 1 is a diagram illustrating an example of calculating the effective integrated temperature. FIG. 1 shows an example in which the harvesting period is predicted on November 12, 2011 for cabbage J planted on September 12, 2011.
図1(A)は、キャベツJの栽培が行われる地域の一つである地域Aに関して、2011年の9月12日から11月11日までの各日の平均気温の実績値と、前年の2010年の11月12日以降の各日の平均気温の実績値とを示す。 FIG. 1 (A) shows the average temperature of each day from September 12th to November 11th, 2011, and the previous year's value for region A, one of the regions where cabbage J is grown. The actual value of the average temperature of each day after November 12, 2010 is shown.
(B)は、(A)に示される各日の平均気温から、生育限界温度である5℃を差し引いた結果である。したがって、(B)は、(A)に示される各日の有効気温を示す。 (B) is the result of subtracting 5 ° C., the growth limit temperature, from the average temperature of each day shown in (A). Therefore, (B) shows the effective temperature of each day shown in (A).
(C)は、(B)に示される有効気温の各日までの積算値を示す。すなわち、(C)は、各日までの有効積算気温の算出結果を示す。予測日である11月12日以降は、2010年の有効気温が積算されることにより、未来の有効積算気温が算出されている。 (C) shows the integrated value up to each day of the effective temperature shown in (B). That is, (C) shows the calculation result of the effective integrated temperature up to each day. After November 12, which is the predicted date, the effective temperature in the future is calculated by integrating the effective temperature in 2010.
但し、未知の平均気温に代入される実績値が、一年分では信頼性に欠ける可能性がある。気象条件は年ごとによって異なる可能性が有るからである。そこで、未来の有効積算気温を算出するためのサンプル数を増加させること考えられる。例えば、過去10年の各日の平均気温の実績値を用いて、未来の有効積算気温を算出することが考えられる。 However, the actual value assigned to the unknown average temperature may be unreliable for one year. This is because weather conditions may vary from year to year. Therefore, it is conceivable to increase the number of samples for calculating the future effective integrated temperature. For example, it is possible to calculate the future effective integrated temperature using the actual value of the average temperature of each day in the past 10 years.
図2は、過去10年の各日の平均気温を用いた有効積算気温の算出例を示す図である。図2においても、図1と同様に、2011年9月12日に定植を行ったキャベツJについて、同年11月12日に収穫期間の予測が実施される例が示されている。 FIG. 2 is a diagram illustrating a calculation example of the effective integrated temperature using the average temperature of each day in the past 10 years. FIG. 2 also shows an example in which the harvest period is predicted on November 12, 2011 for cabbage J planted on September 12, 2011, as in FIG.
図2において、表T1には、2011年11月11までの各日の平均気温の実績値が示されている。また、11月12以降に関しては、過去10年の各年における各日の平均気温の実績値が示されている。 In FIG. 2, the table T1 shows the actual value of the average temperature of each day until November 11, 2011. In addition, for November 12 and after, the actual temperature value of each day in each year of the past 10 years is shown.
表T2には、11月12以降に関しては、過去の各年における各日の平均気温の実績値を用いて有効積算気温を算出することにより、収穫期間が導出された例が示されている。すなわち、有効積算気温が、収穫可能日に対応する有効積算気温である1250℃・日を超える日の翌日と、収穫限界日に対応する有効積算気温である1400℃・日を超える日の翌日とが示されている。更に、収穫可能日と収穫限界日との平均、すなわち、収穫期間の中間の日が、収穫目安日として、最後の行に示されている。 Table T2 shows an example in which the harvest period is derived by calculating the effective integrated temperature using the actual value of the average temperature of each day in each past year for November 12 and thereafter. That is, the day after the day when the effective integrated temperature exceeds 1250 ° C · day, which is the effective integrated temperature corresponding to the harvestable day, and the day after the day when the effective integrated temperature exceeds 1400 ° C · day, which is the effective integrated temperature corresponding to the harvest limit date, It is shown. Furthermore, the average of the harvestable date and the harvest limit date, that is, the day in the middle of the harvest period, is shown in the last line as the harvest target date.
例えば、2010年の各日の平均気温の実績値が用いられた場合、収穫可能日は4月2日であり、収穫限界日は、4月20日であると予測されている。また、2008年の各日の平均気温の実績値が用いられた場合、収穫可能日は、3月18日であり、収穫限界日は、4月10であると予測されている。 For example, when the actual value of the average temperature of each day in 2010 is used, the harvestable date is predicted to be April 2 and the harvest limit date is predicted to be April 20. Moreover, when the actual value of the average temperature of each day in 2008 is used, it is predicted that the harvestable date is March 18 and the harvest limit date is April 10.
図2では、これらの予測日の10年分の平均が、最終的な予測日とされる。予測日の平均とは、各年の予測日の日付を単純に合計し、年数分(例えば、10)で除することにより得られる値である。予測日は、表T2の右端の列に示されている。すなわち、収穫可能日は、3月19日として予測され、収穫限界日は4月10日と予測されている。また、収穫目安日は、3月30日として予測されている。 In FIG. 2, the average of these prediction dates for 10 years is set as the final prediction date. The average predicted date is a value obtained by simply adding the dates of the predicted date of each year and dividing by the number of years (for example, 10). The predicted date is shown in the rightmost column of Table T2. That is, the harvestable date is predicted as March 19, and the harvest limit date is predicted as April 10. The estimated date for harvesting is predicted as March 30th.
図2に示されるように、過去の複数年の各年の予測日の平均を、当年の予測日とすることで、収穫期間の予測に関して、或る年に特有の気象条件による影響を小さくすることができる。 As shown in FIG. 2, by setting the average of the forecast days of each of the past years as the forecast date of the current year, the influence of the weather conditions peculiar to a certain year on the forecast of the harvest period is reduced. be able to.
しかしながら、図3の例では、収穫期間の予測が行われるたびに、未来分の有効積算気温に関して、複数年のデータを利用した計算が必要となる。 However, in the example of FIG. 3, every time the harvest period is predicted, a calculation using data for a plurality of years is required for the effective integrated temperature for the future.
そこで、本実施の形態では、更に工夫が行われる。具体的は、過去の複数年分の各日の有効気温について、各日に一つの代表値を算出し、当該代表値を用いて収穫期間の予測が行われる。 Therefore, further improvements are made in the present embodiment. Specifically, for each day of effective temperatures for a plurality of past years, one representative value is calculated for each day, and the harvest period is predicted using the representative value.
図3は、過去複数年分の各日の有効気温の代表値の算出例を示す図である。図3では、説明の便宜上、2001年から2003年までの過去3年分の実績値を利用する例が示されているが、図2と同様に過去10年分以上の実績値が利用されてもよい。 FIG. 3 is a diagram illustrating a calculation example of a representative value of effective temperatures for each day for the past plural years. For convenience of explanation, FIG. 3 shows an example using actual values for the past three years from 2001 to 2003, but actual values for the past 10 years or more are used as in FIG. Also good.
図3(A)では、過去の各年における各日の平均気温について、同月同日ごとに平均値が算出される。更に、平均値から生育限界温度が差し引かれた温が算出される。当該気温が、過去の複数年分の各日の有効気温の代表値とされる。すなわち、収穫期間の予測実施日以降に関しては、当該代表値の積算値が有効積算気温とされる。 In FIG. 3 (A), an average value is calculated for each day of the same month for the average temperature of each day in the past years. Furthermore, the temperature obtained by subtracting the growth limit temperature from the average value is calculated. The said temperature is made into the representative value of the effective temperature of each day for the past several years. In other words, the integrated value of the representative value is the effective integrated temperature for the harvest period and subsequent execution dates.
しかし、(A)の方法では、日平均気温が生育限界温度を超える年と、生育限界温度未満の年とが混在する日について、生育限界温度を超える分が、生育限界温度を下回る分によって相殺されてしまうという問題がある。 However, in the method of (A), for the day in which the daily average temperature exceeds the growth limit temperature and the year where the growth temperature is less than the growth limit temperature, the amount exceeding the growth limit temperature is offset by the amount below the growth limit temperature. There is a problem of being done.
例えば、12月20日の平均気温に注目すると、2001年及び2002年は、キャベツJの生育限界温度である5℃を超えているが、2003年は、5℃を1.9℃下回っている。これら3つの平均気温の平均値を単純に計算すると、2003年の平均気温によって、本来であれば有効積算気温には影響しないはずの−1.9℃が、有効積算気温に影響を与えてしまう。換言すれば、有効気温の最小値は0であるにも拘わらず、0未満の有効気温が考慮された結果が得られてしまう。 For example, focusing on the average temperature on December 20, 2001 and 2002 exceeded the cabbage J growth limit temperature of 5 ° C, but in 2003 it was 1.9 ° C below 5 ° C. . If the average value of these three average temperatures is simply calculated, the average temperature in 2003 would have an impact on the effective integrated temperature of -1.9 ° C, which would otherwise not affect the effective integrated temperature. . In other words, although the minimum value of the effective temperature is 0, a result in which an effective temperature less than 0 is considered is obtained.
そこで、本実施の形態では、図3(B)に示される方法が採用される。(B)では、まず、各年における各日について、有効気温が算出される。すなわち、各年の各日の平均気温から生育限界温度が差し引かれた0以上の値が算出される。その後、各年における各日の有効気温について、同月同日ごとに平均値が算出される。当該平均値が、過去の複数年分の各日の有効気温の代表値とされる。すなわち、収穫期間の予測実施日以降に関しては、当該代表値の積算値が有効積算気温とされる。 Therefore, in this embodiment, the method shown in FIG. In (B), first, the effective temperature is calculated for each day in each year. That is, a value of 0 or more is calculated by subtracting the growth limit temperature from the average temperature of each day of each year. Thereafter, an average value is calculated for each day of the same month for the effective temperature of each day in each year. The said average value is made into the representative value of the effective temperature of each day for the past several years. In other words, the integrated value of the representative value is the effective integrated temperature for the harvest period and subsequent execution dates.
(B)の方法によれば、2003年の12月19日〜21日のように、生育限界温度未満の平均気温に対する有効気温は、0となる。すなわち、2003年の12月19日〜21日における、生育限界温度未満の温度分は、有効積算気温に影響しなくなる。したがって、(A)における問題点を解消することができる。 According to the method of (B), the effective temperature with respect to the average temperature below the growth limit temperature is 0 as in December 19 to 21 of 2003. That is, the temperature less than the growth limit temperature on December 19 to 21 in 2003 does not affect the effective integrated temperature. Therefore, the problem in (A) can be solved.
(A)の方法と、(B)の方法との違いを具体的に説明する。図4は、生育限界温度未満の日平均気温を無視しない場合と無視する場合との収穫期間の予測結果の違いの例を示す図である。生育限界温度未満の日平均気温を無視しない場合とは、図3(A)の場合をいう。生育限界温度未満の日平均気温を無視する場合とは、図3(B)の場合をいう。 The difference between the method (A) and the method (B) will be specifically described. FIG. 4 is a diagram illustrating an example of the difference in the prediction result of the harvest period between when the daily average temperature below the growth limit temperature is not ignored and when it is ignored. The case where the daily average temperature below the growth limit temperature is not ignored refers to the case of FIG. The case of ignoring the daily average temperature below the growth limit temperature refers to the case of FIG.
図4のグラフG1は、2011年9月12日に定植を行ったキャベツJの有効積算気温の推移を示すグラフである。同年11月12日が収穫期間の予測実施日であるため、11月11日までは、同年の9月12日からの各日の平均気温の実績値に基づく有効気温が積算された結果が一点鎖線によって示されている。11月12日以降の有効積算気温については、図4(A)の方法による算出結果の推移と、(B)の方法による算出結果の推移との双方が示されている。すなわち、(A)の方法による算出結果の推移は、破線によって示され、(B)の方法による算出結果の推移は、実線によって示されている。 The graph G1 of FIG. 4 is a graph which shows transition of the effective integrated temperature of the cabbage J planted on September 12, 2011. Since November 12 of the same year is the estimated implementation date of the harvest period, up to November 11, there is one result of integrating the effective temperature based on the actual value of the average temperature of each day from September 12 of the same year. Indicated by the dashed line. Regarding the effective integrated temperature after November 12, both the transition of the calculation result by the method of FIG. 4A and the transition of the calculation result by the method of (B) are shown. That is, the transition of the calculation result by the method (A) is indicated by a broken line, and the transition of the calculation result by the method (B) is indicated by a solid line.
二つの算出結果の違いをより明確に把握するため、グラフG2には、グラフG1において楕円で囲まれた部分が拡大されて示されている。 In order to grasp the difference between the two calculation results more clearly, the graph G2 shows an enlarged portion surrounded by an ellipse in the graph G1.
グラフG2において、収穫可能日に対応する有効積算気温である1250℃・日を示す直線と各曲線とが交わる日の翌日が、(A)の方法又は(B)の方法による、収穫可能日の予測値となる。また、収穫限界日に対応する有効積算気温である1400℃・日を示す直線と各曲線とが交わる日の翌日が、(A)の方法又は(B)の方法による、収穫可能日の予測値となる。 In graph G2, the day after the day when the curve intersects with the straight line indicating 1250 ° C · day, which is the effective integrated temperature corresponding to the harvestable day, is the harvestable day according to the method (A) or (B). Predicted value. In addition, the day after the day when each curve intersects with the straight line indicating 1400 ° C · day, which is the effective integrated temperature corresponding to the harvest limit date, is the predicted value of the harvestable date according to the method (A) or the method (B). It becomes.
テーブルT3には、それぞれの方法に基づく予測結果が示されている。すなわち、テーブルT3の第一列目は、(A)の方法に基づく予測結果である。第二列目は、(B)の方法に基づく予測結果である。なお、第三列目は、予測結果の妥当性の判断基準として、図2の方法に基づく予測結果が示されている。 In the table T3, prediction results based on the respective methods are shown. That is, the first column of the table T3 is a prediction result based on the method (A). The second column shows the prediction result based on the method (B). In the third column, the prediction result based on the method of FIG. 2 is shown as a criterion for determining the validity of the prediction result.
テーブルT3より明らかなように、(A)の方法に基づく予測結果は、図2の方法に基づく予測結果に対して大きなずれを有している。一方、(B)の方法に基づく予測結果は、図2の方法に基づく予測結果とほぼ等しくなっている。なお、テーブルT3では、(B)の方法に基づく予測結果と、図2の方法に基づく予測結果とは完全に一致しているが、常にそのようになることを示すものではない。 As is clear from the table T3, the prediction result based on the method (A) has a large deviation from the prediction result based on the method of FIG. On the other hand, the prediction result based on the method (B) is substantially equal to the prediction result based on the method of FIG. In the table T3, the prediction result based on the method (B) and the prediction result based on the method of FIG. 2 are completely coincident with each other, but it does not always indicate that this is the case.
図4より明らかなように、図3(B)による方法は、図2の方法に比較して予測結果の妥当性を大きく損なうことなく、計算量を削減できる方法であるといえる。計算量を削減できるのは、図3(B)の方法によれば、一度各日に関して代表値が算出されれば、再び異なる日に収穫時期の予測が行われる場合であっても、改めて代表値を算出する必要は無いからである。すなわち、予測を行う度に、過去複数年分のデータを利用する必要は無いからである。 As is clear from FIG. 4, the method according to FIG. 3B can be said to be a method that can reduce the amount of calculation without significantly impairing the validity of the prediction result as compared with the method of FIG. The amount of calculation can be reduced according to the method of FIG. 3B, once the representative value is calculated for each day, even if the harvest time is predicted again on a different day, This is because it is not necessary to calculate a value. That is, it is not necessary to use data for a plurality of past years each time a prediction is made.
次に、図3(B)において説明した収穫期間の予測方法を実行するコンピュータについて説明する。 Next, a computer that executes the harvest period prediction method described with reference to FIG.
図5は、本発明の実施の形態における農作業支援システムの構成例を示す図である。図5に示される農作業支援システム1において、農作業支援装置10及と情報処理端末20とは、例えば、LAN(Local Area Network)又はインターネット等のネットワークN1を介して通信可能とされている。ネットワークN1の一部又は全部は、無線通信でもよい。
FIG. 5 is a diagram showing a configuration example of the farm work support system in the embodiment of the present invention. In the farm
農作業支援装置10は、上記した方法によって、農作物の収穫期間を予測するコンピュータである。情報処理端末20は、農作業支援装置10に対する入力情報の送信や、農作業支援装置10による処理結果の表示等を実行する。農作業支援装置10に対する入力情報は、例えば、農作業者等のユーザによって、情報処理端末20に入力される。情報処理端末20の一例として、PC、携帯電話、スマートフォン、タブレット型端末等が挙げられる。
The farm
図6は、本発明の実施の形態における農作業支援装置のハードウェア構成例を示す図である。図6の農作業支援装置10は、それぞれバスBで相互に接続されているドライブ装置100、補助記憶装置102、メモリ装置103、CPU104、及びインタフェース装置105等を有する。
FIG. 6 is a diagram illustrating a hardware configuration example of the farm work support apparatus according to the embodiment of the present invention. The farm
農作業支援装置10での処理を実現するプログラムは、記録媒体101によって提供される。プログラムを記録した記録媒体101がドライブ装置100にセットされると、プログラムが記録媒体101からドライブ装置100を介して補助記憶装置102にインストールされる。但し、プログラムのインストールは必ずしも記録媒体101より行う必要はなく、ネットワークを介して他のコンピュータよりダウンロードするようにしてもよい。補助記憶装置102は、インストールされたプログラムを格納すると共に、必要なファイルやデータ等を格納する。
A program that realizes processing in the farm
メモリ装置103は、プログラムの起動指示があった場合に、補助記憶装置102からプログラムを読み出して格納する。CPU104は、メモリ装置103に格納されたプログラムに従って農作業支援装置10に係る機能を実行する。インタフェース装置105は、ネットワークに接続するためのインタフェースとして用いられる。
The
なお、記録媒体101の一例としては、CD−ROM、DVDディスク、又はUSBメモリ等の可搬型の記録媒体が挙げられる。また、補助記憶装置102の一例としては、HDD(Hard Disk Drive)又はフラッシュメモリ等が挙げられる。記録媒体101及び補助記憶装置102のいずれについても、コンピュータ読み取り可能な記録媒体に相当する。
An example of the
なお、農作業支援装置10に表示装置や、キーボード及びマウス等の入力装置が接続されてもよい。
Note that a display device or an input device such as a keyboard and a mouse may be connected to the farm
図7は、本発明の実施の形態における農作業支援装置の機能構成例を示す図である。図7において、農作業支援装置10は、代表値算出条件受付部11、有効気温算出部12、代表値算出部13、予測条件受付部14、収穫時期予測部15、及び出力部16等を有する。これら各部は、農作業支援装置10にインストールされた一以上のプログラムが、CPU104に実行させる処理により実現される。農作業支援装置10は、また、過去気温履歴記憶部21、当年気温履歴記憶部22、及び代表値記憶部23等を利用する。これら各記憶部は、補助記憶装置102、又は農作業支援装置10にネットワークを介して接続される記憶装置等を用いて実現可能である。
FIG. 7 is a diagram illustrating a functional configuration example of the farm work support apparatus according to the embodiment of the present invention. In FIG. 7, the farm
過去気温履歴記憶部21は、過去の複数年の各年について、1月1日から12月31日までの各日の平均気温を示す気温データを記憶する。当年気温履歴記憶部22は、収穫時期の予測対象とされる栽培が行われている当年の1月1日又は栽培開始日以降の各日の平均気温の実績値を示す気温データを記憶する。当年気温履歴記憶部22の記憶内容は、例えば、毎日更新される。
The past temperature
代表値算出条件受付部11は、過去の各日の有効気温の代表値を算出するための前提条件となるパラメータの入力を受け付ける。例えば、生育限界温度や、使用する過去の気温データの範囲等が、当該パラメータの一例に該当する。
The representative value calculation
有効気温算出部12は、過去気温履歴記憶部21が記憶する、過去の複数年の各年各日の平均気温に関して、生育限界温度を超える分の気温である有効気温を算出する。代表値算出部13は、有効気温算出部12によって各年各日について算出された有効気温について、同月同日ごとに平均値を算出する。当該平均値は、過去の各日の有効気温の代表値として代表値記憶部23に記憶される。したがって、代表値記憶部23は、1月1日から12月31日までの各日の代表値を記憶する。
The effective
予測条件受付部14は、収穫可能日又は収穫限界日等の収穫時期を予測するための前提条件となるパラメータの入力を受け付ける。例えば、栽培開始日や目標となる有効積算気温等が、当該パラメータの一例に該当する。
The prediction
収穫時期予測部15は、予測条件受付部14が受け付けたパラメータ、当年気温履歴記憶部22が記憶する情報、及び代表値記憶部23が記憶する情報等に基づいて収穫時期を予測する。出力部16は、収穫時期の予測結果を出力する。
The harvest
以下、農作業支援装置10が実行する処理手順について説明する。図8は、過去複数年の各日の有効気温の代表値の算出処理の処理手順の一例を説明するためのフローチャートである。
Hereinafter, a processing procedure executed by the farm
ステップS101において、代表値算出条件受付部11は、生育限界温度の値を受け付ける。続いて、代表値算出条件受付部11は、代表値の算出に使用する過去の気温データの範囲を受け付ける(S102)。例えば、過去N年分のデータに基づいて、代表値を算出する場合、Nの値が受け付けられる。
In step S101, the representative value calculation
なお、ステップS101及びS102は、例えば、情報処理端末20においてユーザによって入力された値が、農作業支援装置10において受信されることにより実行される。または、生育限界温度及びNの値は、予め補助記憶装置102等に記憶されていてもよい。更に、生育限界温度の代わりに、農作物の品種が入力されてもよい。この場合、例えば、農作業支援装置10には、農作物の品種ごとに生育限界温度が記憶されていてもよい。そうすることにより、入力された品種に基づいて、生育限界温度を特定することができる。
Note that steps S <b> 101 and S <b> 102 are executed when, for example, a value input by the user at the
続いて、有効気温算出部12は、過去気温履歴記憶部21に記憶されている、過去N年の各年の各日の平均気温から、生育限界温度を差し引くことにより当該各年各日の有効気温を算出する(S103)。但し、生育限界温度未満の日については、有効気温として0が付与される。すなわち、ステップS103では、過去N年の各年の各日の平均気温について、生育限界温度を超える分の気温が算出される。例えば、図3(B)の2001年から2003年の各列に示される値は、ステップS103における算出結果の一例である。
Subsequently, the effective
続いて、代表値算出部13は、過去N年の各年各日に関して算出された有効気温について、同月同日ごとに平均値を算出する(S104)。その結果、各日について、有効気温の代表値が算出される。例えば、図3(B)の有効気温の平均の列に示される値は、ステップS104における算出結果の一例である。
Subsequently, the representative
続いて、代表値算出部13は、閏年のみに存在する2月29日に対する代表値を算出する(S105)。当該代表値は、例えば、ステップS104において算出された、2月28日に対する代表値と、3月1日に対する代表値との平均を求めることにより算出される。なお、ステップS105は、収穫期間の予測が行われる当年が閏年である場合に限って実行されてもよい。
Subsequently, the representative
続いて、代表値算出部13は、各日の有効気温の代表値を、代表値記憶部23に記憶する(S106)。
Subsequently, the representative
なお、図8の処理は、生育限界温度、及び過去の気温データの範囲が変わらない限り、一回実行されればよい。すなわち、後述される図9の処理が実行されるたびに実行される必要はない。 The process of FIG. 8 may be executed once as long as the growth limit temperature and the range of past temperature data are not changed. That is, it does not need to be executed every time the process of FIG. 9 described later is executed.
続いて、図8の処理によって生成された代表値を利用して、収穫時期を予測する処理について説明する。 Next, a process for predicting the harvest time using the representative value generated by the process of FIG. 8 will be described.
図9は、収穫時期の予測処理の処理手順の一例を説明するためのフローチャートである。 FIG. 9 is a flowchart for explaining an example of a processing procedure of harvest time prediction processing.
ステップS201において、予測条件受付部14は、定植日又は播種日等の栽培開始日を受け付ける。なお、栽培開始日は、農作物の品種ごとに有効積算気温の算出の起点となる日が選択されればよい。
In step S201, the prediction
続いて、予測条件受付部14は、目標有効積算気温TAを受け付ける(S202)。目標有効積算気温TAとは、収穫時期の予測対象となる品種について、収穫可能日に対応する有効積算気温若しくは収穫限界日に対応する有効積算気温、又はその双方をいう。
Subsequently, the prediction
なお、ステップS201及びS202は、例えば、情報処理端末20においてユーザによって入力された値が、農作業支援装置10において受信されることにより実行される。または、栽培開始日及び目標有効積算気温TAの値は、予め補助記憶装置102等に記憶されていてもよい。更に、目標有効積算気温TAの代わりに、農作物の品種が入力されてもよい。この場合、例えば、農作業支援装置10には、農作物の品種ごとに目標有効積算気温TAが記憶されていてもよい。そうすることにより、入力された品種に基づいて、目標有効積算気温TAを特定することができる。
In addition, step S201 and S202 are performed when the value input by the user in the
続いて、収穫時期予測部15は、当年気温履歴記憶部22に記憶されている、当年の各日の平均気温に基づいて、昨日までの有効積算気温を算出する(S203)。昨日とは、図9の処理が実行されている日の前の日をいう。なお、算出された有効積算気温は、変数TBに代入される。
Subsequently, the harvest
続いて、収穫時期予測部15は、変数Dateに対して、本日の日付を代入する(S204)。本日とは、図9の処理が実行されている日をいう。また、日付とは、月日をいい、年は含まない。
Subsequently, the harvest
続いて、収穫時期予測部15は、変数TBの値に代表値(Date)を加算した結果を、変数TBに代入する(S205)。代表値(Date)は、Dateが示す日付に対して代表値記憶部23に記憶されている代表値を示す。したがって、ステップS205では、これまでに算出された有効積算気温に対して、Dateが示す日付に対応する代表値が加算されることにより、有効積算気温が更新される。
Subsequently, the harvest
続いて、収穫時期予測部15は、変数TBの値が、目標有効積算気温TA以上であるか否かを判定する(S206)。変数TBの値が、目標有効積算気温TA未満である場合(S206でNo)、変数Dateの日付が1日進められ(S207)、ステップS205以降が繰り返される。
Subsequently, the harvest
変数TBの値が、目標有効積算気温TA以上となると(S206でYes)、出力部16は、その時点において変数Dateが示す日付の翌日を、目標有効積算気温TAの到達日として出力する(S207)。例えば、目標有効積算気温TAが、収穫可能日に対応する有効積算気温であれば、変数Dateが示す日付の翌日が、収穫可能日の予測値として出力される。また、目標有効積算気温TAが、収穫限界日に対応する有効積算気温であれば、変数Dateが示す日付の翌日が、収穫限界日の予測値として出力される。目標有効積算気温TAに、収穫可能日及び収穫限界日のそれぞれに対応する有効積算気温が含まれる場合、ステップS205以降は、それぞれごとに実行される。
When the value of the variable TB is equal to or higher than the target effective integrated temperature TA (Yes in S206), the
なお、ステップS208における出力形態は、所定のものに限定されない。例えば、出力対象の日付を示す情報が、補助記憶装置102等に記憶されてもよいし、情報処理端末20に送信されてもよい。
The output form in step S208 is not limited to a predetermined one. For example, information indicating an output target date may be stored in the
また、当年気温履歴記憶部22と、過去気温履歴記憶部21とは、明確に区別されていなくてもよい。すなわち、これらの記憶部は、一つの記憶部に統合されてもよい。特に、収穫時期の予測対象とされる農作物の栽培が、年を跨ぐ場合、当年気温履歴記憶部22と、過去気温履歴記憶部21とは統合されていた方が、都合が良いと考えられる。
Further, the current year temperature
上述したように、本実施の形態によれば、農作物の収穫時期の予測精度の向上の可能性を高めることができる。 As described above, according to the present embodiment, it is possible to increase the possibility of improving the prediction accuracy of the crop harvest time.
なお、本実施の形態において、代表値記憶部23は、第一の記憶部の一例である。過去気温履歴記憶部21は、第二の記憶部の一例である。
In the present embodiment, the representative
以上、本発明の実施例について詳述したが、本発明は斯かる特定の実施形態に限定されるものではなく、特許請求の範囲に記載された本発明の要旨の範囲内において、種々の変形・変更が可能である。 As mentioned above, although the Example of this invention was explained in full detail, this invention is not limited to such specific embodiment, In the range of the summary of this invention described in the claim, various deformation | transformation・ Change is possible.
1 農作業支援システム
10 農作業支援装置
11 代表値算出条件受付部
12 有効気温算出部
13 代表値算出部
14 予測条件受付部
15 収穫時期予測部
16 出力部
20 情報処理端末
21 過去気温履歴記憶部
22 当年気温履歴記憶部
23 代表値記憶部
100 ドライブ装置
101 記録媒体
102 補助記憶装置
103 メモリ装置
104 CPU
105 インタフェース装置
B バス
DESCRIPTION OF
105 Interface device B bus
Claims (4)
前記農作物の栽培開始日から予測実施日前までの各日の平均気温の実測値に基づく有効積算気温に対して、前記平均値を積算することにより、前記農作物の収穫時期の予測値を算出する処理をコンピュータが実行する収穫時期予測方法。 Using the first storage unit that stores the average value for the same day in the same month for the effective temperature, which is the temperature exceeding the growth limit temperature of the crop for which the harvest time is predicted, of the average temperature of each day of the past multiple years And
A process of calculating a predicted value of the harvest time of the crop by integrating the average value with respect to an effective integrated temperature based on an actual measurement value of the average temperature of each day from the cultivation start date to the prediction implementation date. The harvest time prediction method that the computer executes.
前記複数年の同月同日ごとに、算出された有効気温の平均値を算出し、
算出された平均値を前記第一の記憶部に記憶する処理を前記コンピュータが実行する請求項1記載の収穫時期予測方法。 With respect to the average temperature of each day in each of the past multiple years stored in the second storage unit, an effective temperature that is a temperature exceeding the growth limit temperature of the crop is calculated,
Calculate the average value of the effective temperatures calculated for the same day in the same month of the multiple years,
The harvest time prediction method according to claim 1, wherein the computer executes a process of storing the calculated average value in the first storage unit.
前記農作物の栽培開始日から予測実施日前までの各日の平均気温の実測値に基づく有効積算気温に対して、前記平均値を積算することにより、前記農作物の収穫時期の予測値を算出する処理をコンピュータに実行させるプログラム。 Using the first storage unit that stores the average value for the same day in the same month for the effective temperature, which is the temperature exceeding the growth limit temperature of the crop for which the harvest time is predicted, of the average temperature of each day of the past multiple years And
A process of calculating a predicted value of the harvest time of the crop by integrating the average value with respect to an effective integrated temperature based on an actual measurement value of the average temperature of each day from the cultivation start date to the prediction implementation date. A program that causes a computer to execute.
前記複数年の同月同日ごとに、算出された有効気温の平均値を算出し、
算出された平均値を前記第一の記憶部に記憶する処理を前記コンピュータに実行させる請求項3記載のプログラム。 With respect to the average temperature of each day in each of the past multiple years stored in the second storage unit, an effective temperature that is a temperature exceeding the growth limit temperature of the crop is calculated,
Calculate the average value of the effective temperatures calculated for the same day in the same month of the multiple years,
The program according to claim 3, wherein the computer executes a process of storing the calculated average value in the first storage unit.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012058011A JP5853780B2 (en) | 2012-03-14 | 2012-03-14 | Harvest time prediction method and program |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012058011A JP5853780B2 (en) | 2012-03-14 | 2012-03-14 | Harvest time prediction method and program |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2013191107A true JP2013191107A (en) | 2013-09-26 |
JP5853780B2 JP5853780B2 (en) | 2016-02-09 |
Family
ID=49391252
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2012058011A Active JP5853780B2 (en) | 2012-03-14 | 2012-03-14 | Harvest time prediction method and program |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5853780B2 (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105075692A (en) * | 2015-04-22 | 2015-11-25 | 华北水利水电大学 | Winter wheat annual ridge culture and furrow irrigation yield prediction and irrigation quality assessment method |
CN111428924A (en) * | 2020-03-20 | 2020-07-17 | 北京禾水源智慧农装科技股份有限公司 | Plant growth period prediction method and device, electronic equipment and storage medium |
JP2021051331A (en) * | 2018-09-21 | 2021-04-01 | 国立研究開発法人農業・食品産業技術総合研究機構 | Program for environmental assessment inside greenhouse, method for environmental assessment inside greenhouse and apparatus for environmental assessment inside greenhouse |
CN113128779A (en) * | 2021-04-29 | 2021-07-16 | 华东交通大学 | Decision support system for optimal harvest period of oranges |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004030449A (en) * | 2002-06-27 | 2004-01-29 | Nissan Chem Ind Ltd | System, method and program for managing non-tilled direct sowing and computer readable recording medium with the same program |
JP2007310463A (en) * | 2006-05-16 | 2007-11-29 | Hitachi Software Eng Co Ltd | Farm field management support method and system |
-
2012
- 2012-03-14 JP JP2012058011A patent/JP5853780B2/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004030449A (en) * | 2002-06-27 | 2004-01-29 | Nissan Chem Ind Ltd | System, method and program for managing non-tilled direct sowing and computer readable recording medium with the same program |
JP2007310463A (en) * | 2006-05-16 | 2007-11-29 | Hitachi Software Eng Co Ltd | Farm field management support method and system |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105075692A (en) * | 2015-04-22 | 2015-11-25 | 华北水利水电大学 | Winter wheat annual ridge culture and furrow irrigation yield prediction and irrigation quality assessment method |
JP2021051331A (en) * | 2018-09-21 | 2021-04-01 | 国立研究開発法人農業・食品産業技術総合研究機構 | Program for environmental assessment inside greenhouse, method for environmental assessment inside greenhouse and apparatus for environmental assessment inside greenhouse |
JP7377517B2 (en) | 2018-09-21 | 2023-11-10 | 国立研究開発法人農業・食品産業技術総合研究機構 | Greenhouse environment evaluation program, greenhouse environment evaluation method, and greenhouse environment evaluation device |
CN111428924A (en) * | 2020-03-20 | 2020-07-17 | 北京禾水源智慧农装科技股份有限公司 | Plant growth period prediction method and device, electronic equipment and storage medium |
CN113128779A (en) * | 2021-04-29 | 2021-07-16 | 华东交通大学 | Decision support system for optimal harvest period of oranges |
Also Published As
Publication number | Publication date |
---|---|
JP5853780B2 (en) | 2016-02-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5853780B2 (en) | Harvest time prediction method and program | |
Goldberg et al. | Tempo and mode in plant breeding system evolution | |
JP2020194767A (en) | Soc correction method and apparatus, battery management system and storage medium | |
KR101820567B1 (en) | Agricultural work support method and agricultural work support device | |
US20140298321A1 (en) | Installation control method and installation control apparatus | |
JP2015176339A (en) | Device, method and program for prediction of crop harvest time | |
JP2016161314A (en) | Assimilative data determination device, weather prediction device, assimilative data determination method and program | |
JP2008077561A (en) | Energy prediction method, energy prediction device and program | |
Su et al. | Reference gene selection for quantitative real-time PCR (qRT-PCR) expression analysis in Galium aparine L. | |
JP5772980B2 (en) | Agricultural work support method and agricultural work support apparatus | |
US11625635B2 (en) | Computer-readable storage medium storing control program, control method, and control device | |
JP5811188B2 (en) | Harvest support method and harvest support device | |
JP2020205073A (en) | Dataset normalization for predicting dataset attribute | |
JP4550758B2 (en) | Sensor data detection time correction device | |
KR102194871B1 (en) | Method of predicting crop yield and apparatus for embodying the same | |
Zhang et al. | Genetic progress battles climate variability: drivers of soybean yield gains in China from 2006 to 2020 | |
JP6222822B2 (en) | Deterioration function calculation device, deterioration rate estimation system, deterioration function calculation method, and program | |
JP6863144B2 (en) | Cultivation support programs, equipment, and methods | |
JP6203313B2 (en) | Feature selection device, feature selection method, and program | |
JP6795134B1 (en) | Power management device | |
JP6493904B2 (en) | Parameter selection method, parameter selection program, and parameter selection device | |
JP7089285B2 (en) | Crop growth period estimation device and growth period estimation method | |
WO2021106046A1 (en) | Maturation prediction system | |
US9779110B2 (en) | Placement device, placement method and recording medium, that stores placement program information processing device | |
KR20190128102A (en) | Method, system and non-transitory computer-readable recording medium for providing prediction about breeding |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20141112 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20150818 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20150819 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20151013 |
|
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: 20151110 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20151123 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5853780 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |