JP6171732B2 - Crop discrimination device, crop discrimination program and crop discrimination method - Google Patents
Crop discrimination device, crop discrimination program and crop discrimination method Download PDFInfo
- Publication number
- JP6171732B2 JP6171732B2 JP2013175908A JP2013175908A JP6171732B2 JP 6171732 B2 JP6171732 B2 JP 6171732B2 JP 2013175908 A JP2013175908 A JP 2013175908A JP 2013175908 A JP2013175908 A JP 2013175908A JP 6171732 B2 JP6171732 B2 JP 6171732B2
- Authority
- JP
- Japan
- Prior art keywords
- intensity
- component
- specific wavelength
- crop
- wavelength band
- 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
- 238000012850 discrimination method Methods 0.000 title description 3
- 238000000034 method Methods 0.000 claims description 61
- 238000004364 calculation method Methods 0.000 claims description 46
- 238000003860 storage Methods 0.000 claims description 24
- 230000007423 decrease Effects 0.000 claims description 13
- 238000003384 imaging method Methods 0.000 claims description 5
- 241000219310 Beta vulgaris subsp. vulgaris Species 0.000 description 89
- 235000002595 Solanum tuberosum Nutrition 0.000 description 89
- 244000061456 Solanum tuberosum Species 0.000 description 89
- 235000021536 Sugar beet Nutrition 0.000 description 89
- 241000209140 Triticum Species 0.000 description 88
- 235000021307 Triticum Nutrition 0.000 description 88
- 238000001228 spectrum Methods 0.000 description 38
- 238000002474 experimental method Methods 0.000 description 33
- 238000012545 processing Methods 0.000 description 22
- 238000010586 diagram Methods 0.000 description 19
- 238000004519 manufacturing process Methods 0.000 description 14
- 230000003595 spectral effect Effects 0.000 description 6
- 239000000284 extract Substances 0.000 description 5
- 241000196324 Embryophyta Species 0.000 description 4
- 238000007796 conventional method Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 239000004973 liquid crystal related substance Substances 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 235000012015 potatoes Nutrition 0.000 description 2
- 238000009826 distribution Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
Images
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Description
本発明は、作物判別装置、作物判別プログラム及び作物判別方法に関する。 The present invention relates to a crop discrimination device, a crop discrimination program, and a crop discrimination method.
圃場に植えられた作物を判別する従来技術が存在する。例えば、人工衛星のスペクトルカメラによって得られた判別対象の圃場における太陽光の反射光に含まれる特定の波長帯域の成分の強度と、作物毎に設定された基準となる特定の波長帯域の成分の強度とを比較して、圃場に植えられた作物を判別する従来技術がある。 There is a conventional technique for discriminating crops planted in a field. For example, the intensity of a specific wavelength band component included in reflected sunlight of a discrimination target field obtained by a spectrum camera of an artificial satellite and a specific wavelength band component that is a reference set for each crop There is a conventional technique for comparing the strength and discriminating the crops planted in the field.
しかしながら、以下に説明するように、上述した従来技術では、精度良く作物を判別することができないという問題がある。 However, as described below, the above-described conventional technique has a problem that crops cannot be accurately identified.
すなわち、従来技術では、スペクトルカメラによって、所定の範囲の領域(例えば、10m×10mの領域)ごとに、太陽光の反射光に含まれる複数の波長帯域の成分の強度が含まれるピクセルを取得する。そして、従来技術では、判別対象の圃場における複数のピクセルのそれぞれに含まれる特定の波長帯域の成分の強度の平均値を、判別対象の圃場における特定の波長帯域の成分の強度として算出する。そして、従来技術では、算出した強度と、作物ごとに設定された基準となる特定の波長帯域の成分の強度とを比較して、判別対象の圃場に植えられた作物を判別する。 That is, in the related art, a pixel including intensities of components in a plurality of wavelength bands included in reflected sunlight is acquired for each predetermined range of region (for example, a region of 10 m × 10 m) by a spectrum camera. . In the related art, the average value of the components of the specific wavelength band included in each of the plurality of pixels in the discrimination target field is calculated as the intensity of the component of the specific wavelength band in the discrimination target field. In the related art, the calculated intensity is compared with the intensity of a component in a specific wavelength band as a reference set for each crop, and the crop planted in the field to be determined is determined.
しかしながら、判別対象の圃場のピクセルの一部には、圃場以外の領域が少なくとも一部に含まれてしまう場合がある。この場合には、判別対象の圃場以外の領域を少なくとも一部に含むピクセルに含まれる特定の波長帯域の成分の強度は、判別対象の圃場以外の領域を含まない場合における特定の波長帯域の成分の強度と大きく異なる場合がある。このような場合であっても、従来技術では、一部のピクセルに判別対象の圃場以外の領域が含まれるような複数のピクセルのそれぞれに含まれる特定の波長帯域の成分の強度の平均値を、判別対象の圃場における特定の波長帯域の成分の強度として算出する。この結果、算出された特定の波長帯域の成分の強度は、判別対象の圃場以外の領域を含まない場合における特定の波長帯域の成分の強度と大きく異なる場合がある。したがって、上述した従来技術による作物の判別では、精度良く作物を判別することができない場合がある。 However, there are cases where at least a part of the region other than the field is included in a part of the pixels of the field to be determined. In this case, the intensity of the component in the specific wavelength band included in the pixel including at least a part of the region other than the field to be discriminated is the component of the specific wavelength band in the case where the region other than the field to be discriminated is not included. There are cases where it differs greatly from the strength. Even in such a case, in the conventional technology, the average value of the intensities of the components in the specific wavelength band included in each of the plurality of pixels in which some pixels include regions other than the field to be discriminated is obtained. The intensity of the component in the specific wavelength band in the field to be discriminated is calculated. As a result, the calculated intensity of the component in the specific wavelength band may be significantly different from the intensity of the component in the specific wavelength band when the region other than the field to be discriminated is not included. Therefore, there is a case where the crop cannot be accurately discriminated by the above-described conventional crop discrimination.
1つの側面では、本発明は、作物の判別の精度の低下を抑制することを目的とする。 In one aspect, an object of the present invention is to suppress a decrease in the accuracy of crop discrimination.
1態様では、作物判別装置は、記憶部と、生成部と、特定部と、算出部と、判別部とを有する。記憶部は、ある地域をカメラによって撮像することにより得られた、一部のピクセルが該ある地域の少なくとも一部の領域を含む複数のピクセルのそれぞれにおける複数の特定の波長帯域の成分の強度を記憶する。生成部は、記憶部に記憶された記憶内容に基づいて、特定の波長帯域ごとに、特定の波長帯域の成分の強度と、該強度に対応するピクセルの個数との関係を示す情報を生成する。特定部は、生成部により生成された情報に基づいて、特定の波長帯域ごとに、ピクセルの個数が最大である強度を特定する。算出部は、生成部により生成された情報に基づいて、特定の波長帯域ごとに、次の処理を行う。すなわち、算出部は、特定部により特定された強度を含む、最大のピクセルの個数からピクセルの個数が減少する範囲内の特定の波長帯域の成分の強度に対応するピクセルにおける特定の波長帯域の成分の強度の平均値を算出する。判別部は、算出部により特定の波長帯域ごとに算出された平均値に基づいて、ある地域に所定の作物が植えられているか否かを判別する。 In one aspect, the crop discrimination device includes a storage unit, a generation unit, a specifying unit, a calculation unit, and a discrimination unit. The storage unit obtains the intensities of components in a plurality of specific wavelength bands in each of a plurality of pixels, in which some pixels include at least a part of the region, obtained by imaging a region with a camera. Remember. The generation unit generates, for each specific wavelength band, information indicating the relationship between the intensity of the component in the specific wavelength band and the number of pixels corresponding to the intensity, based on the storage content stored in the storage unit. . The specifying unit specifies the intensity with the maximum number of pixels for each specific wavelength band based on the information generated by the generating unit. The calculation unit performs the following processing for each specific wavelength band based on the information generated by the generation unit. That is, the calculation unit includes a component in a specific wavelength band in a pixel corresponding to the intensity of a component in a specific wavelength band within a range in which the number of pixels decreases from the maximum number of pixels, including the intensity specified by the specifying unit. The average value of the intensity is calculated. The determining unit determines whether or not a predetermined crop is planted in a certain area based on the average value calculated for each specific wavelength band by the calculating unit.
作物の判別の精度の低下を抑制することができる。 A decrease in the accuracy of crop discrimination can be suppressed.
以下に、本願の開示する作物判別装置、作物判別プログラム及び作物判別方法の実施例を図面に基づいて詳細に説明する。なお、実施例は開示の技術を限定するものではない。 Hereinafter, embodiments of a crop discrimination device, a crop discrimination program, and a crop discrimination method disclosed in the present application will be described in detail with reference to the drawings. The embodiments do not limit the disclosed technology.
[作物判別装置10の機能構成の一例]
図1は、実施例に係る作物判別装置10の機能構成の一例を示す図である。図1の例に示すように、作物判別装置10は、入力部11、出力部12、記憶部13及び制御部14を有する。本実施例に係る作物判別装置10は、判別の対象の圃場に植えられている作物が、小麦、馬鈴薯、テンサイの3種類の作物のいずれであるかを判別する。
[Example of functional configuration of crop discrimination device 10]
FIG. 1 is a diagram illustrating an example of a functional configuration of a
入力部11は、各種情報を制御部14に入力する。例えば、入力部11は、作物判別装置10のユーザから後述の作物判別処理を実行するための指示を受け付けて、受け付けた指示を制御部14に入力する。入力部11のデバイスの一例としては、マウスやキーボードなどが挙げられる。
The input unit 11 inputs various information to the
出力部12は、各種の情報を出力する。例えば、出力部12は、後述の出力制御部14eの制御により、どの圃場にどの種類の作物が植えられているかを示す情報を表示する。この情報については後述する。出力部12のデバイスの一例としては、液晶ディスプレイなどが挙げられる。
The
記憶部13は、各種情報を記憶する。例えば、記憶部13は、基準強度データ13a、圃場別スペクトルデータ13b、判別用圃場別スペクトルデータ13cを記憶する。
The
基準強度データ13aは、判別の対象の圃場の後述の特定の波長帯域の成分の強度と比較されるデータである。図2は、基準強度データ13aの一例を示す図である。図2に示す基準強度データ13aには、「小麦」について後述する「G」の範囲の最小値「68.3」及び最大値「92.4」と、後述する「R」の範囲の最小値「61.3」及び最大値「82.9」とが登録されている。また、図2に示す基準強度データ13aには、「小麦」について後述する「NIR」の範囲の最小値「85.0」及び最大値「115.0」と、後述する「SWIR」の範囲の最小値「56.7」及び最大値「76.7」とが登録されている。
The
ここで、「小麦」についての「G」は、小麦が植えられた複数の圃場における太陽光の反射光に含まれる緑成分の波長帯域の強度の平均値である。例えば、「小麦」についての「G」は、波長が550nmの緑成分の強度の平均値である。なお、緑成分の波長帯域は、500nm〜590nmである。また、「小麦」についての「R」は、小麦が植えられた複数の圃場における太陽光の反射光に含まれる赤成分の波長帯域の強度の平均値である。例えば、「小麦」についての「R」は、波長が640nmの赤成分の強度の平均値である。なお、赤成分の波長帯域は、610nm〜680nmである。また、「小麦」についての「NIR」は、小麦が植えられた複数の圃場における太陽光の反射光に含まれる近赤外成分の波長帯域の強度の平均値である。例えば、「小麦」についての「NIR」は、波長が840nmの近赤外成分の強度の平均値である。なお、近赤外成分の波長帯域は、780nm〜890nmである。また、「小麦」についての「SWIR」は、小麦が植えられた複数の圃場における太陽光の反射光に含まれる短波長赤外成分の波長帯域の強度の平均値である。例えば、「小麦」についての「SWIR」は、例えば、波長が1660nmの近赤外成分の強度の平均値である。なお、短波長赤外成分の波長帯域は、1580nm〜1750nmである。 Here, “G” for “wheat” is an average value of the intensity of the wavelength band of the green component contained in the reflected light of sunlight in a plurality of fields where wheat is planted. For example, “G” for “wheat” is an average value of the intensity of the green component having a wavelength of 550 nm. The wavelength band of the green component is 500 nm to 590 nm. Further, “R” for “wheat” is an average value of the intensity of the wavelength band of the red component contained in the reflected light of sunlight in a plurality of fields where wheat is planted. For example, “R” for “wheat” is an average value of the intensity of the red component having a wavelength of 640 nm. The wavelength band of the red component is 610 nm to 680 nm. In addition, “NIR” for “wheat” is an average value of the intensity in the wavelength band of the near-infrared component included in the reflected sunlight of a plurality of fields in which wheat is planted. For example, “NIR” for “wheat” is an average value of the intensity of near-infrared components having a wavelength of 840 nm. The wavelength band of the near infrared component is 780 nm to 890 nm. In addition, “SWIR” for “wheat” is an average value of the intensity of the wavelength band of the short-wavelength infrared component included in the reflected light of sunlight in a plurality of fields in which wheat is planted. For example, “SWIR” for “wheat” is, for example, an average value of the intensity of near-infrared components having a wavelength of 1660 nm. The wavelength band of the short wavelength infrared component is 1580 nm to 1750 nm.
そして、「小麦」についての「G」の範囲の最小値は、「小麦」についての「G」を0.85倍した値である。また、「小麦」についての「G」の範囲の最大値は、「小麦」についての「G」を1.15倍した値である。また、「小麦」についての「R」の範囲の最小値は、「小麦」についての「R」を0.85倍した値である。また、「小麦」についての「R」の範囲の最大値は、「小麦」についての「R」を1.15倍した値である。また、「小麦」についての「NIR」の範囲の最小値は、「小麦」についての「NIR」を0.85倍した値である。また、「小麦」についての「NIR」の範囲の最大値は、「小麦」についての「NIR」を1.15倍した値である。また、「小麦」についての「SWIR」の範囲の最小値は、「小麦」についての「SWIR」を0.85倍した値である。また、「小麦」についての「SWIR」の範囲の最大値は、「小麦」についての「SWIR」を1.15倍した値である。 And the minimum value of the range of “G” for “wheat” is a value obtained by multiplying “G” for “wheat” by 0.85. In addition, the maximum value of the “G” range for “wheat” is a value obtained by multiplying “G” for “wheat” by 1.15. Further, the minimum value in the range of “R” for “wheat” is a value obtained by multiplying “R” for “wheat” by 0.85. The maximum value of the “R” range for “wheat” is a value obtained by multiplying “R” for “wheat” by 1.15. Further, the minimum value of the “NIR” range for “wheat” is a value obtained by multiplying “NIR” for “wheat” by 0.85. In addition, the maximum value of the “NIR” range for “wheat” is a value obtained by multiplying “NIR” for “wheat” by 1.15. Further, the minimum value in the range of “SWIR” for “wheat” is a value obtained by multiplying “SWIR” for “wheat” by 0.85. In addition, the maximum value of the “SWIR” range for “wheat” is a value obtained by multiplying “SWIR” for “wheat” by 1.15.
例えば、図2の例において、ある圃場における太陽光の反射光に含まれる特定の成分が、以下に説明する4つの条件を満たすときには、この圃場に、小麦が植えられていると考えられる。1つ目の条件は、ある圃場における太陽光の反射光に含まれる緑成分の強度が68.3以上、かつ、92.4以下であるという条件である。2つ目の条件は、ある圃場における太陽光の反射光に含まれる赤成分の強度が61.3以上、かつ、82.9以下であるという条件である。3つ目の条件は、ある圃場における太陽光の反射光に含まれる近赤外成分の強度が85.0以上、かつ、115.0以下であるという条件である。4つ目の条件は、ある圃場における太陽光の反射光に含まれる短波長赤外成分の強度が56.7以上、かつ、76.7以下であるという条件である。このようなある圃場に小麦が植えられているか判定するための4つの条件を以下の説明では、「小麦についての4つの条件」と表記する場合がある。 For example, in the example of FIG. 2, when a specific component included in the reflected light of sunlight in a certain field satisfies the following four conditions, it is considered that wheat is planted in this field. The first condition is that the intensity of the green component contained in the reflected sunlight of a certain field is 68.3 or more and 92.4 or less. The second condition is that the intensity of the red component contained in the reflected sunlight of a certain field is 61.3 or more and 82.9 or less. The third condition is that the intensity of the near-infrared component included in the reflected sunlight of a certain field is 85.0 or more and 115.0 or less. The fourth condition is a condition that the intensity of the short-wavelength infrared component contained in the reflected sunlight of a certain field is 56.7 or more and 76.7 or less. In the following description, the four conditions for determining whether wheat is planted in a certain field may be referred to as “four conditions for wheat”.
また、図2に示す基準強度データ13aには、「馬鈴薯」について後述する「G」の範囲の最小値「53.8」及び最大値「72.8」と、後述する「R」の範囲の最小値「42.0」及び最大値「56.9」とが登録されている。また、図2に示す基準強度データ13aには、「馬鈴薯」について後述する「NIR」の範囲の最小値「85.0」及び最大値「115.0」と、後述する「SWIR」の範囲の最小値「66.5」及び最大値「90.0」とが登録されている。
In addition, the
「馬鈴薯」についての「G」は、馬鈴薯が植えられた複数の圃場における太陽光の反射光に含まれる緑成分の波長帯域の強度の平均値である。例えば、「馬鈴薯」についての「G」は、波長が550nmの緑成分の強度の平均値である。また、「馬鈴薯」についての「R」は、馬鈴薯が植えられた複数の圃場における太陽光の反射光に含まれる赤成分の波長帯域の強度の平均値である。例えば、「馬鈴薯」についての「R」は、波長が640nmの赤成分の強度の平均値である。また、「馬鈴薯」についての「NIR」は、馬鈴薯が植えられた複数の圃場における太陽光の反射光に含まれる近赤外成分の波長帯域の強度の平均値である。例えば、「馬鈴薯」についての「NIR」は、波長が840nmの近赤外成分の強度の平均値である。また、「馬鈴薯」についての「SWIR」は、馬鈴薯が植えられた複数の圃場における太陽光の反射光に含まれる短波長赤外成分の波長帯域の強度の平均値である。例えば、「馬鈴薯」についての「SWIR」は、例えば、波長が1660nmの近赤外成分の強度の平均値である。 “G” for “potato” is the average value of the intensity of the wavelength band of the green component contained in the reflected sunlight of the plurality of fields in which the potato is planted. For example, “G” for “potato” is the average value of the intensity of the green component having a wavelength of 550 nm. In addition, “R” for “potato” is an average value of the intensity of the wavelength band of the red component contained in the reflected sunlight of the plurality of fields in which the potato is planted. For example, “R” for “potato” is the average value of the intensity of the red component having a wavelength of 640 nm. “NIR” for “potato” is the average value of the intensity of the wavelength band of the near-infrared component included in the reflected light of sunlight in a plurality of fields in which the potato is planted. For example, “NIR” for “potato” is the average value of the intensity of the near infrared component having a wavelength of 840 nm. In addition, “SWIR” for “potato” is an average value of the intensity of the wavelength band of the short wavelength infrared component contained in the reflected light of sunlight in a plurality of fields in which the potato is planted. For example, “SWIR” for “potato” is, for example, the average value of the intensity of the near-infrared component having a wavelength of 1660 nm.
そして、「馬鈴薯」についての「G」の範囲の最小値は、「馬鈴薯」についての「G」を0.85倍した値である。また、「馬鈴薯」についての「G」の範囲の最大値は、「馬鈴薯」についての「G」を1.15倍した値である。また、「馬鈴薯」についての「R」の範囲の最小値は、「馬鈴薯」についての「R」を0.85倍した値である。また、「馬鈴薯」についての「R」の範囲の最大値は、「馬鈴薯」についての「R」を1.15倍した値である。また、「馬鈴薯」についての「NIR」の範囲の最小値は、「馬鈴薯」についての「NIR」を0.85倍した値である。また、「馬鈴薯」についての「NIR」の範囲の最大値は、「馬鈴薯」についての「NIR」を1.15倍した値である。また、「馬鈴薯」についての「SWIR」の範囲の最小値は、「馬鈴薯」についての「SWIR」を0.85倍した値である。また、「馬鈴薯」についての「SWIR」の範囲の最大値は、「馬鈴薯」についての「SWIR」を1.15倍した値である。 And the minimum value of the range of “G” for “potato” is a value obtained by multiplying “G” for “potato” by 0.85. Further, the maximum value of the range of “G” for “potato” is a value obtained by multiplying “G” for “potato” by 1.15. The minimum value of the range of “R” for “potato” is a value obtained by multiplying “R” for “potato” by 0.85. The maximum value of the range of “R” for “potato” is a value obtained by multiplying “R” for “potato” by 1.15. Further, the minimum value of the “NIR” range for “potato” is a value obtained by multiplying “NIR” for “potato” by 0.85. In addition, the maximum value of the “NIR” range for “potato” is a value obtained by multiplying “NIR” for “potato” by 1.15. The minimum value of the “SWIR” range for “potato” is a value obtained by multiplying “SWIR” for “potato” by 0.85. The maximum value of the “SWIR” range for “potato” is a value obtained by multiplying “SWIR” for “potato” by 1.15.
例えば、図2の例において、ある圃場における太陽光の反射光に含まれる特定の成分が、以下に説明する4つの条件を満たすときには、この圃場に、馬鈴薯が植えられていると考えられる。1つ目の条件は、ある圃場における太陽光の反射光に含まれる緑成分の強度が53.8以上、かつ、72.8以下であるという条件である。2つ目の条件は、ある圃場における太陽光の反射光に含まれる赤成分の強度が42.0以上、かつ、56.9以下であるという条件である。3つ目の条件は、ある圃場における太陽光の反射光に含まれる近赤外成分の強度が85.0以上、かつ、115.0以下であるという条件である。4つ目の条件は、ある圃場における太陽光の反射光に含まれる短波長赤外成分の強度が66.5以上、かつ、90.0以下であるという条件である。このようなある圃場に馬鈴薯が植えられているか判定するための4つの条件を以下の説明では、「馬鈴薯についての4つの条件」と表記する場合がある。 For example, in the example of FIG. 2, when a specific component included in reflected sunlight of a certain field satisfies four conditions described below, it is considered that a potato is planted in this field. The first condition is that the intensity of the green component contained in the reflected sunlight of a certain field is 53.8 or more and 72.8 or less. The second condition is that the intensity of the red component contained in the reflected sunlight of a certain field is 42.0 or more and 56.9 or less. The third condition is that the intensity of the near-infrared component included in the reflected sunlight of a certain field is 85.0 or more and 115.0 or less. The fourth condition is a condition that the intensity of the short wavelength infrared component contained in the reflected sunlight of a certain field is 66.5 or more and 90.0 or less. In the following description, the four conditions for determining whether a potato is planted in a certain field may be referred to as “four conditions for potato”.
また、図2に示す基準強度データ13aには、「テンサイ」について後述する「G」の範囲の最小値「61.1」及び最大値「82.7」と、後述する「R」の範囲の最小値「49.9」及び最大値「67.5」とが登録されている。また、図2に示す基準強度データ13aには、「テンサイ」について後述する「NIR」の範囲の最小値「85.0」及び最大値「115.0」と、後述する「SWIR」の範囲の最小値「56.8」及び最大値「76.9」とが登録されている。
Further, the
「テンサイ」についての「G」は、テンサイが植えられた複数の圃場における太陽光の反射光に含まれる緑成分の波長帯域の強度の平均値である。例えば、「テンサイ」についての「G」は、波長が550nmの緑成分の強度の平均値である。また、「テンサイ」についての「R」は、テンサイが植えられた複数の圃場における太陽光の反射光に含まれる赤成分の波長帯域の強度の平均値である。例えば、「テンサイ」についての「R」は、波長が640nmの赤成分の強度の平均値である。また、「テンサイ」についての「NIR」は、テンサイが植えられた複数の圃場における太陽光の反射光に含まれる近赤外成分の波長帯域の強度の平均値である。例えば、「テンサイ」についての「NIR」は、波長が840nmの近赤外成分の強度の平均値である。また、「テンサイ」についての「SWIR」は、テンサイが植えられた複数の圃場における太陽光の反射光に含まれる短波長赤外成分の波長帯域の強度の平均値である。例えば、「テンサイ」についての「SWIR」は、例えば、波長が1660nmの近赤外成分の強度の平均値である。 “G” for “sugar beet” is the average value of the intensity of the wavelength band of the green component contained in the reflected light of sunlight in a plurality of fields where sugar beet is planted. For example, “G” for “sugar beet” is the average value of the intensities of the green component having a wavelength of 550 nm. Also, “R” for “sugar beet” is the average value of the intensity of the wavelength band of the red component contained in the reflected light of sunlight in a plurality of fields where the sugar beet is planted. For example, “R” for “sugar beet” is the average value of the intensity of the red component having a wavelength of 640 nm. Further, “NIR” for “sugar beet” is an average value of the intensity of the wavelength band of the near infrared component included in the reflected light of sunlight in a plurality of fields where sugar beet is planted. For example, “NIR” for “sugar beet” is the average value of the intensity of the near infrared component having a wavelength of 840 nm. In addition, “SWIR” for “sugar beet” is an average value of the intensity of the short wavelength infrared component wavelength band included in the reflected light of sunlight in a plurality of fields where sugar beet is planted. For example, “SWIR” for “sugar beet” is, for example, the average value of near-infrared component intensities having a wavelength of 1660 nm.
そして、「テンサイ」についての「G」の範囲の最小値は、「テンサイ」についての「G」を0.85倍した値である。また、「テンサイ」についての「G」の範囲の最大値は、「テンサイ」についての「G」を1.15倍した値である。また、「テンサイ」についての「R」の範囲の最小値は、「テンサイ」についての「R」を0.85倍した値である。また、「テンサイ」についての「R」の範囲の最大値は、「テンサイ」についての「R」を1.15倍した値である。また、「テンサイ」についての「NIR」の範囲の最小値は、「テンサイ」についての「NIR」を0.85倍した値である。また、「テンサイ」についての「NIR」の範囲の最大値は、「テンサイ」についての「NIR」を1.15倍した値である。また、「テンサイ」についての「SWIR」の範囲の最小値は、「テンサイ」についての「SWIR」を0.85倍した値である。また、「テンサイ」についての「SWIR」の範囲の最大値は、「テンサイ」についての「SWIR」を1.15倍した値である。 The minimum value of the “G” range for “sugar beet” is a value obtained by multiplying “G” for “sugar beet” by 0.85. The maximum value of the “G” range for “sugar beet” is a value obtained by multiplying “G” for “sugar beet” by 1.15. Further, the minimum value of the “R” range for “sugar beet” is a value obtained by multiplying “R” for “sugar beet” by 0.85. In addition, the maximum value of the range of “R” for “sugar beet” is a value obtained by multiplying “R” for “sugar beet” by 1.15. Further, the minimum value of the “NIR” range for “sugar beet” is a value obtained by multiplying “NIR” for “sugar beet” by 0.85. In addition, the maximum value of the “NIR” range for “sugar beet” is a value obtained by multiplying “NIR” for “sugar beet” by 1.15. Further, the minimum value in the range of “SWIR” for “sugar beet” is a value obtained by multiplying “SWIR” for “sugar beet” by 0.85. The maximum value of the “SWIR” range for “sugar beet” is a value obtained by multiplying “SWIR” for “sugar beet” by 1.15.
例えば、図2の例において、ある圃場における太陽光の反射光に含まれる特定の成分が、以下に説明する4つの条件を満たすときには、この圃場に、テンサイが植えられていると考えられる。1つ目の条件は、ある圃場における太陽光の反射光に含まれる緑成分の強度が61.1以上、かつ、82.7以下であるという条件である。2つ目の条件は、ある圃場における太陽光の反射光に含まれる赤成分の強度が49.9以上、かつ、67.5以下であるという条件である。3つ目の条件は、ある圃場における太陽光の反射光に含まれる近赤外成分の強度が85.0以上、かつ、115.0以下であるという条件である。4つ目の条件は、ある圃場における太陽光の反射光に含まれる短波長赤外成分の強度が56.8以上、かつ、76.9以下であるという条件である。このようなある圃場にテンサイが植えられているか判定するための4つの条件を以下の説明では、「テンサイについての4つの条件」と表記する場合がある。 For example, in the example of FIG. 2, when a specific component included in reflected sunlight of a certain field satisfies the following four conditions, it is considered that sugar beet is planted in the field. The first condition is that the intensity of the green component contained in the reflected sunlight of a certain field is 61.1 or more and 82.7 or less. The second condition is that the intensity of the red component contained in the reflected sunlight of a certain field is 49.9 or more and 67.5 or less. The third condition is that the intensity of the near-infrared component included in the reflected sunlight of a certain field is 85.0 or more and 115.0 or less. The fourth condition is a condition that the intensity of the short wavelength infrared component included in the reflected sunlight of a certain field is 56.8 or more and 76.9 or less. In the following description, the four conditions for determining whether sugar beets are planted in a certain field may be referred to as “four conditions for sugar beet”.
圃場別スペクトルデータ13bには、圃場毎に、人工衛星のスペクトルカメラによって、圃場が撮像された結果得られた複数のピクセルのそれぞれにおける太陽光の反射光に含まれる緑成分、赤成分、近赤外成分及び短波長赤外成分の各強度が登録されている。一例としては、圃場別スペクトルデータ13bには、太陽光の反射光に含まれる波長が550nmの緑成分、波長が640nmの赤成分、波長が840nmの近赤外成分、波長が1660nmの短波長赤外成分の各強度が登録されている。
The field-
図3は、圃場別スペクトルデータ13bの一例を示す図である。図3の例に示す圃場別スペクトルデータ13bは、「圃場ID」、「ピクセルID」、「G」、「R」、「NIR」及び「SWIR」の各項目を有する。圃場別スペクトルデータ13bの1レコードには、1つの圃場の1つのピクセルにおける太陽光の反射光に含まれる緑成分、赤成分、近赤外成分及び短波長赤外成分の各強度が予め登録されている。
FIG. 3 is a diagram showing an example of the field-
「圃場ID」の項目には、圃場を識別するためのID(Identification)が予め登録されている。「ピクセルID」の項目には、「圃場ID」の項目に登録された圃場の領域を少なくとも一部に含むピクセルを識別するためのIDが予め登録されている。「G」の項目には、「ピクセルID」の項目に登録されたIDが示すピクセルにおける太陽光の反射光に含まれる緑成分の強度が予め登録されている。「R」の項目には、「ピクセルID」の項目に登録されたIDが示すピクセルにおける太陽光の反射光に含まれる赤成分の強度が予め登録されている。「NIR」の項目には、「ピクセルID」の項目に登録されたIDが示すピクセルにおける太陽光の反射光に含まれる近赤外成分の強度が予め登録されている。「SWIR」の項目には、「ピクセルID」の項目に登録されたIDが示すピクセルにおける太陽光の反射光に含まれる短波長赤外成分の強度が予め登録されている。 In the item “Field ID”, an ID (Identification) for identifying the field is registered in advance. In the “pixel ID” item, an ID for identifying a pixel that includes at least a part of the field area registered in the “field ID” item is registered in advance. In the item “G”, the intensity of the green component included in the reflected sunlight of the pixel indicated by the ID registered in the item “pixel ID” is registered in advance. In the item “R”, the intensity of the red component included in the reflected sunlight of the pixel indicated by the ID registered in the item “pixel ID” is registered in advance. In the “NIR” item, the intensity of the near-infrared component included in the reflected sunlight of the pixel indicated by the ID registered in the “pixel ID” item is registered in advance. In the item “SWIR”, the intensity of the short wavelength infrared component included in the reflected sunlight of the pixel indicated by the ID registered in the item “pixel ID” is registered in advance.
例えば、図3の例に示す圃場別スペクトルデータ13bは、ID「AA」が示す圃場の領域を少なくとも一部に含むピクセルが、ID「AAA」が示すピクセル、ID「AAB」が示すピクセル、及び、ID「AAC」が示すピクセルであることを示す。また、圃場別スペクトルデータ13bは、ID「AAA」が示すピクセルにおける太陽光の反射光に含まれる緑成分、赤成分、近赤外成分、短波長赤外成分のそれぞれの強度が、「90」、「81」、「103」、「61」である場合を示す。また、圃場別スペクトルデータ13bは、ID「AAB」が示すピクセルにおける太陽光の反射光に含まれる緑成分、赤成分、近赤外成分、短波長赤外成分のそれぞれの強度が、「87」、「71」、「100」、「71」である場合を示す。また、圃場別スペクトルデータ13bは、ID「AAC」が示すピクセルにおける太陽光の反射光に含まれる緑成分、赤成分、近赤外成分、短波長赤外成分のそれぞれの強度が、「85」、「72」、「101」、「66」である場合を示す。
For example, the field-specific
また、図3の例に示す圃場別スペクトルデータ13bは、ID「BB」が示す圃場の領域を少なくとも一部に含むピクセルが、ID「BAA」が示すピクセル、ID「BAB」が示すピクセル、及び、ID「BAC」が示すピクセルであることを示す。また、圃場別スペクトルデータ13bは、ID「BAA」が示すピクセルにおける太陽光の反射光に含まれる緑成分、赤成分、近赤外成分、短波長赤外成分のそれぞれの強度が、「71」、「50」、「101」、「57」である場合を示す。また、圃場別スペクトルデータ13bは、ID「BAB」が示すピクセルにおける太陽光の反射光に含まれる緑成分、赤成分、近赤外成分、短波長赤外成分のそれぞれの強度が、「70」、「51」、「100」、「58」である場合を示す。また、圃場別スペクトルデータ13bは、ID「BAC」が示すピクセルにおける太陽光の反射光に含まれる緑成分、赤成分、近赤外成分、短波長赤外成分のそれぞれの強度が、「71」、「52」、「100」、「57」である場合を示す。
Further, the field-
判別用圃場別スペクトルデータ13cについては後述する。 The field-specific spectrum data 13c for determination will be described later.
記憶部13は、例えば、フラッシュメモリなどの半導体メモリ素子、または、ハードディスク、光ディスクなどの記憶装置である。
The
制御部14は、各種の処理手順を規定したプログラムや制御データを格納するための内部メモリを有し、これらによって種々の処理を実行する。図1に示すように、制御部14は、生成部14aと、特定部14bと、算出部14cと、判別部14dと、出力制御部14eとを有する。
The
生成部14aは、圃場別スペクトルデータ13bに基づいて、特定の波長帯域ごとに、特定の波長帯域の成分の強度と、強度に対応するピクセルの個数との関係を示す情報であるヒストグラムを生成する。
The generation unit 14a generates, for each specific wavelength band, a histogram that is information indicating the relationship between the intensity of a component in a specific wavelength band and the number of pixels corresponding to the intensity, based on the field-specific
生成部14aの一態様について説明する。例えば、生成部14aは、入力部11から作物判別処理を実行する指示が入力されると、圃場別スペクトルデータ13bの「圃場ID」の項目に登録された圃場を識別するためのIDを全て取得する。なお、生成部14aは、重複するIDについては、1つのIDとして取得する。そして、生成部14aは、取得したIDのうち、未選択のIDを1つずつ選択する処理を行う。なお、生成部14aにより選択されたIDが示す圃場は、植えられている作物が判別される対象の圃場(判別対象の圃場)である。そして、生成部14aは、未選択のIDを1つ選択するたびに、以下で説明する処理を行う。
One aspect of the generation unit 14a will be described. For example, when an instruction to execute the crop discrimination process is input from the input unit 11, the generation unit 14a acquires all IDs for identifying the fields registered in the item “Field ID” of the field-
すなわち、生成部14aは、圃場別スペクトルデータ13bの全てのレコードの中から、選択したIDが「圃場ID」の項目に登録された全てのレコードを特定する。例えば、図3の例に示す圃場別スペクトルデータ13bが記憶部13に記憶され、生成部14aがID「AA」を選択した場合には、生成部14aは、ID「AA」が「圃場ID」の項目に登録された3つのレコードを特定する。
In other words, the generation unit 14a identifies all the records in which the selected ID is registered in the item “Field ID” from among all the records of the field-
そして、生成部14aは、特定した全てのレコードのそれぞれの「G」の項目に登録された緑成分の強度を全て取得する。例えば、図3の例に示す圃場別スペクトルデータ13bが記憶部13に記憶され、生成部14aが、ID「AA」が「圃場ID」の項目に登録された3つのレコードを特定した場合には、生成部14aは、次の処理を行う。すなわち、生成部14aは、緑成分の3つの強度「90」、「87」、「85」を取得する。
And the production | generation part 14a acquires all the intensity | strengths of the green component registered into each item of "G" of all the specified records. For example, when the field-
そして、生成部14aは、取得した全ての緑成分の強度のうち、同一の緑成分の強度の数を取得した緑成分の強度ごとに計数する。ここで、生成部14aは、例えば、「N」よりも大きく、かつ、「N+10」以下の緑成分の強度については、同一の緑成分の強度「N+5」として扱う。例えば、生成部14aは、例えば、「120」よりも大きく、かつ、「130」以下の緑成分の強度については、同一の緑成分の強度「125」として扱う。なお、計数した緑成分の強度の数は、同一の緑成分の強度を有するピクセルの個数である。そして、生成部14aは、横軸を緑成分の強度、縦軸をピクセルの数とするヒストグラムを生成する。すなわち、生成部14aは、選択したIDが示す圃場における緑成分の強度についてのヒストグラムを生成する。 And the production | generation part 14a counts the number of the intensity | strengths of the same green component among the intensity | strengths of all the acquired green components for every intensity | strength of the acquired green component. Here, for example, the intensity of the green component greater than “N” and equal to or less than “N + 10” is handled as the intensity “N + 5” of the same green component. For example, the generation unit 14a treats the intensity of the green component that is greater than “120” and equal to or less than “130” as the intensity “125” of the same green component. The counted number of green component intensities is the number of pixels having the same green component intensity. Then, the generation unit 14a generates a histogram with the horizontal axis representing the green component intensity and the vertical axis representing the number of pixels. That is, the generation unit 14a generates a histogram regarding the intensity of the green component in the field indicated by the selected ID.
また、生成部14aは、特定した全てのレコードのそれぞれの「R」の項目に登録された赤成分の強度を全て取得する。そして、生成部14aは、取得した全ての赤成分の強度のうち、同一の赤成分の強度の数を、取得した赤成分の強度ごとに計数する。ここで、生成部14aは、例えば、「N」よりも大きく、かつ、「N+10」以下の赤成分の強度については、同一の赤成分の強度「N+5」として扱う。なお、計数した赤成分の強度の数は、同一の赤成分の強度を有するピクセルの個数である。そして、生成部14aは、横軸を赤成分の強度、縦軸をピクセルの個数とするヒストグラムを生成する。すなわち、生成部14aは、選択したIDが示す圃場における赤成分の強度についてのヒストグラムを生成する。 Further, the generation unit 14a acquires all the intensities of the red components registered in the respective “R” items of all the specified records. And the production | generation part 14a counts the number of the intensity | strength of the same red component among the intensity | strengths of all the acquired red components for every intensity | strength of the acquired red component. Here, for example, the intensity of the red component greater than “N” and equal to or less than “N + 10” is handled as the intensity “N + 5” of the same red component. The counted number of red component intensities is the number of pixels having the same red component intensity. Then, the generation unit 14a generates a histogram with the horizontal axis representing the red component intensity and the vertical axis representing the number of pixels. That is, the generation unit 14a generates a histogram about the intensity of the red component in the field indicated by the selected ID.
また、生成部14aは、特定した全てのレコードのそれぞれの「NIR」の項目に登録された近赤外成分の強度を全て取得する。そして、生成部14aは、取得した全ての近赤外成分の強度のうち、同一の近赤外成分の強度の数を、取得した近赤外成分の強度ごとに計数する。ここで、生成部14aは、例えば、「N」よりも大きく、かつ、「N+10」以下の近赤外成分の強度については、同一の近赤外成分の強度「N+5」として扱う。なお、計数した近赤外成分の強度の数は、同一の近赤外成分の強度を有するピクセルの個数である。そして、生成部14aは、横軸を近赤外成分の強度、縦軸をピクセルの個数とするヒストグラムを生成する。すなわち、生成部14aは、選択したIDが示す圃場における近赤外成分の強度についてのヒストグラムを生成する。 Further, the generation unit 14a acquires all the intensities of the near-infrared components registered in the “NIR” items of all the specified records. And the production | generation part 14a counts the number of the intensity | strength of the same near-infrared component among the intensity | strengths of all the acquired near-infrared components for every intensity | strength of the acquired near-infrared component. Here, for example, the intensity of the near-infrared component greater than “N” and equal to or less than “N + 10” is handled as the intensity “N + 5” of the same near-infrared component. Note that the counted number of near-infrared component intensities is the number of pixels having the same near-infrared component intensity. Then, the generation unit 14a generates a histogram with the horizontal axis indicating the intensity of the near infrared component and the vertical axis indicating the number of pixels. That is, the generation unit 14a generates a histogram about the intensity of the near infrared component in the field indicated by the selected ID.
また、生成部14aは、特定した全てのレコードのそれぞれの「SWIR」の項目に登録された短波長赤外成分の強度を全て取得する。そして、生成部14aは、取得した全ての短波長赤外成分の強度のうち、同一の短波長赤外成分の強度の数を、取得した短波長赤外成分の強度ごとに計数する。ここで、生成部14aは、例えば、「N」よりも大きく、かつ、「N+10」以下の短波長赤外成分の強度については、同一の短波長成分の強度「N+5」として扱う。なお、計数した短波長赤外成分の強度の数は、同一の短波長赤外成分の強度を有するピクセルの個数である。そして、生成部14aは、横軸を短波長赤外成分の強度、縦軸をピクセルの個数とするヒストグラムを生成する。すなわち、生成部14aは、選択したIDが示す圃場における短波長赤外成分の強度についてのヒストグラムを生成する。 Further, the generation unit 14a acquires all the intensities of the short-wavelength infrared components registered in the “SWIR” items of all the specified records. And the production | generation part 14a counts the number of the intensity | strength of the same short wavelength infrared component among the intensity | strengths of all the acquired short wavelength infrared components for every intensity | strength of the acquired short wavelength infrared component. Here, for example, the intensity of the short wavelength infrared component that is larger than “N” and equal to or less than “N + 10” is handled as the intensity “N + 5” of the same short wavelength component. The counted number of short wavelength infrared component intensities is the number of pixels having the same short wavelength infrared component intensity. Then, the generation unit 14a generates a histogram with the horizontal axis representing the intensity of the short wavelength infrared component and the vertical axis representing the number of pixels. That is, the production | generation part 14a produces | generates the histogram about the intensity | strength of the short wavelength infrared component in the field which selected ID shows.
図4は、実施例に係る生成部14aにより生成されたヒストグラムの一例を示す図である。図4の例は、横軸を近赤外成分の強度、縦軸をピクセルの個数とするヒストグラム20を示す。図4に示すヒストグラム20では、ピクセルの個数についてのピークが2つある。図4の例に示すヒストグラム20では、1つ目のピークにおける近赤外成分の強度は、「65」であり、2つ目のピークにおける近赤外成分の強度は、「125」である。ここで、2つのピークのうち、ピクセルの個数が少ない方のピーク(近赤外成分の強度「65」に対応するピーク)は、次のような理由により発生することが考えられる。すなわち、判別対象の圃場のピクセルの一部に、判別対象の圃場以外の領域が少なくとも一部に含まれてしまった場合に、発生すると考えられる。例えば、図5や図6に示すように、判別対象の圃場30のピクセル35の一部に、判別対象の圃場以外のあぜ道40や他の圃場41の領域が少なくとも一部に含まれてしまった場合に、発生すると考えられる。本実施例では、判別対象の圃場のピクセルの一部に、判別対象の圃場以外の領域が少なくとも一部に含まれてしまった場合であっても、精度良く作物の種別を判別できるように、後述する各処理を行う。なお、上述した図5及び図6は、判別対象の圃場のピクセルの一部に、判別対象の圃場以外の領域が少なくとも一部に含まれてしまった場合の一例を説明するための図である。
FIG. 4 is a diagram illustrating an example of a histogram generated by the generation unit 14a according to the embodiment. The example of FIG. 4 shows a
生成部14aは、上述した処理を、IDを選択するたびに行う。これにより、圃場別スペクトルデータ13bにIDが登録された全ての圃場における、緑成分、赤成分、近赤外成分、短波長赤外成分のそれぞれについてのヒストグラムが生成される。
The generation unit 14a performs the above-described process every time an ID is selected. Thereby, the histogram about each of a green component, a red component, a near infrared component, and a short wavelength infrared component is produced | generated in all the fields in which ID was registered into the
特定部14bは、生成部14aにより生成されたヒストグラムに基づいて、特定の波長帯域ごとに、ピクセルの個数が最大である強度を特定する。 The specifying unit 14b specifies the intensity at which the number of pixels is the maximum for each specific wavelength band based on the histogram generated by the generating unit 14a.
特定部14bの一態様について説明する。例えば、特定部14bは、生成部14aにより選択されたIDが示す圃場における緑成分、赤成分、近赤外成分及び短波長赤外成分についての4つのヒストグラムが生成されるたびに、ヒストグラムごとに、次の処理を行う。すなわち、特定部14bは、ヒストグラムにおいて、ピクセルの個数についてピークとなる特定の波長帯域の成分の強度を抽出する。ここで、ヒストグラムにおいてピークが複数ある場合には、特定部14bは、複数のピークを抽出する。なお、ピークとは、例えば、ヒストグラムが示す、特定の波長帯域の成分の強度と、特定の波長帯域の成分の強度に対応するピクセルの個数との関係をグラフで表した場合に、極大となる点を指す。例えば、図4の例に示すヒストグラム20が生成された場合には、特定部14bは、ピクセルの個数についてピークとなる2つの近赤外成分の強度「65」及び「125」を抽出する。
One aspect of the specifying unit 14b will be described. For example, the specifying unit 14b generates a histogram for each histogram every time four histograms are generated for the green component, the red component, the near infrared component, and the short wavelength infrared component in the field indicated by the ID selected by the generating unit 14a. The following processing is performed. That is, the specifying unit 14b extracts the intensity of a component in a specific wavelength band that peaks for the number of pixels in the histogram. Here, when there are a plurality of peaks in the histogram, the specifying unit 14b extracts a plurality of peaks. Note that the peak is a maximum when, for example, a graph shows the relationship between the intensity of a component in a specific wavelength band and the number of pixels corresponding to the intensity of a component in a specific wavelength band indicated by a histogram. Point to a point. For example, when the
そして、特定部14bは、抽出した特定の波長帯域の成分の強度が複数であるか否かを判定する。抽出した特定の波長帯域の成分の強度が複数であると判定した場合には、特定部14bは、次の処理を行う。すなわち、特定部14bは、ヒストグラムが示す特定の波長帯域の成分の強度とピクセルの個数との関係に基づいて、抽出した複数の特定の波長帯域の成分の強度の中から、対応するピクセルの個数が最大である特定の波長帯域の成分の強度を特定する。例えば、図4の例に示すヒストグラム20が生成されて、ピクセルの個数についてピークとなる2つの近赤外成分の強度「65」及び「125」を抽出した場合には、特定部14bは、次の処理を行う。すなわち、特定部14bは、2つの近赤外成分の強度「65」及び「125」の中から、ピクセルの個数が最大(60個)である近赤外成分の強度「125」を特定する。
Then, the specifying unit 14b determines whether or not the intensity of the extracted component in the specific wavelength band is plural. If it is determined that the extracted component of the specific wavelength band has a plurality of intensities, the specifying unit 14b performs the following process. That is, the specifying unit 14b determines the number of corresponding pixels from the extracted component intensities in the specific wavelength band based on the relationship between the intensity of the component in the specific wavelength band indicated by the histogram and the number of pixels. The intensity of the component in the specific wavelength band where is the maximum is specified. For example, when the
算出部14cは、ヒストグラムに基づいて、特定の波長帯域ごとに、次の処理を行う。すなわち、算出部14cは、特定部14bにより特定された強度を含む、最大のピクセルの個数からピクセルの個数が減少する範囲内の特定の波長帯域の成分の強度に対応するピクセルにおける特定の波長帯域の成分の強度の平均値を算出する。 The calculation unit 14c performs the following processing for each specific wavelength band based on the histogram. That is, the calculation unit 14c includes the specific wavelength band in the pixel corresponding to the intensity of the component in the specific wavelength band within the range in which the number of pixels decreases from the maximum number of pixels, including the intensity specified by the specifying unit 14b. The average value of the intensities of the components is calculated.
算出部14cの一態様について説明する。例えば、算出部14cは、特定部14bにより特定の波長帯域の成分の強度が特定されると、次の処理を行う。すなわち、算出部14cは、ヒストグラムが示す特定の波長帯域の成分の強度とピクセルの個数との関係において、特定部14bにより特定された強度に対応するピクセルの個数からピクセルの個数が減少する範囲内の特定の波長帯域の成分の強度を特定する。図7は、実施例に係る算出部14cが実行する処理の一例について説明するための図である。例えば、先の図4の例に示すヒストグラム20において、特定部14bにより近赤外成分の強度「125」が特定された場合には、算出部14cは、次の処理を行う。すなわち、算出部14cは、図7の例に示すように、強度「125」に対応するピクセルの個数「60」からピクセルの個数が減少する範囲内の近赤外成分の強度「115」、「105」、「125」、「135」、「145」、「155」を特定する。このように、算出部14cは、特定部14bにより特定された強度から、横軸方向に順に、各強度に対応するピクセルの個数を参照し、特定部14bにより特定された強度に対応するピクセルの個数から、ピクセルの個数が減少傾向にある強度を全て特定する。この場合に、算出部14cは、特定部14bにより特定された強度も特定する。
One aspect of the calculation unit 14c will be described. For example, when the specifying unit 14b specifies the intensity of a component in a specific wavelength band, the calculating unit 14c performs the following process. That is, the calculation unit 14c has a range in which the number of pixels decreases from the number of pixels corresponding to the intensity specified by the specifying unit 14b in the relationship between the intensity of the component in the specific wavelength band indicated by the histogram and the number of pixels. The intensity of the component in the specific wavelength band is specified. FIG. 7 is a diagram for explaining an example of a process executed by the calculation unit 14c according to the embodiment. For example, in the
そして、算出部14cは、特定した全ての強度のそれぞれに対応するピクセルにおける特定の波長帯域の成分の強度の平均値を算出する。例えば、図7の例に示すように、強度「115」、「105」、「125」、「135」、「145」、「155」を特定した場合には、算出部14cは、次の処理を行う。すなわち、算出部14cは、図7の例に示すヒストグラム20から、強度「115」に対応するピクセルの個数「5」を取得し、強度「105」に対応するピクセルの個数「2」を取得し、強度「125」に対応するピクセルの個数「60」を取得する。また、算出部14cは、図7の例に示すヒストグラム20から、強度「135」に対応するピクセルの個数「50」を取得し、強度「145」に対応するピクセルの個数「10」を取得し、強度「155」に対応するピクセルの個数「5」を取得する。そして、算出部14cは、強度「115」、「105」、「125」、「135」、「145」、「155」に対応するピクセルの強度の平均値を算出する。すなわち、算出部14cは、「224」(115×5+105×2+125×60+135×50+145×10+155×5)/(5+2+60+50+10+5)を算出する。
Then, the calculating unit 14c calculates the average value of the intensities of the components in the specific wavelength band in the pixels corresponding to all the specified intensities. For example, as shown in the example of FIG. 7, when the intensities “115”, “105”, “125”, “135”, “145”, and “155” are specified, the calculation unit 14 c I do. That is, the calculation unit 14c obtains the number of pixels “5” corresponding to the intensity “115” and the number of pixels “2” corresponding to the intensity “105” from the
そして、平均値を算出すると、算出部14cは、生成部14aにより選択されたIDと、算出した平均値とを対応付けて記憶部13に記憶された判別用圃場別スペクトルデータ13cに登録する。
When the average value is calculated, the calculation unit 14 c registers the ID selected by the generation unit 14 a and the calculated average value in the discrimination-specific field-specific spectrum data 13 c stored in the
図8は、判別用圃場別スペクトルデータ13cのデータ構造の一例を示す図である。図8の例に示す判別用圃場別スペクトルデータ13cは、「圃場ID」、「G」、「R」、「NIR」、「SWIR」の項目を有する。算出部14cは、平均値を算出すると、生成部14aにより選択されたIDが「圃場ID」の項目に登録されたレコードがあるか判別用圃場別スペクトルデータ13cの全レコードを検索する。検索の結果、生成部14aにより選択されたIDが「圃場ID」の項目に登録されたレコードが得られた場合には、算出部14cは、次のような処理を行う。例えば、算出部14cは、算出した平均値が、緑成分の強度の平均値である場合には、得られたレコードの「G」の項目に、算出した平均値を登録する。また、算出部14cは、算出した平均値が、赤成分の強度の平均値である場合には、得られたレコードの「R」の項目に、算出した平均値を登録する。また、算出部14cは、算出した平均値が、近赤外成分の強度の平均値である場合には、得られたレコードの「NIR」の項目に、算出した平均値を登録する。また、算出部14cは、算出した平均値が、短波長赤外成分の強度の平均値である場合には、得られたレコードの「SWIR」の項目に、算出した平均値を登録する。 FIG. 8 is a diagram illustrating an example of a data structure of the spectrum data 13c for each discrimination field. The field-specific spectral data 13c for discrimination shown in the example of FIG. 8 includes items of “field ID”, “G”, “R”, “NIR”, and “SWIR”. When calculating the average value, the calculation unit 14c searches all records of the field-specific spectrum data 13c for determination whether there is a record in which the ID selected by the generation unit 14a is registered in the item “field ID”. As a result of the search, when a record in which the ID selected by the generation unit 14a is registered in the item “field ID” is obtained, the calculation unit 14c performs the following process. For example, when the calculated average value is the average value of the intensity of the green component, the calculation unit 14c registers the calculated average value in the item “G” of the obtained record. Further, when the calculated average value is the average value of the red component intensity, the calculation unit 14c registers the calculated average value in the item “R” of the obtained record. In addition, when the calculated average value is the average value of the near-infrared component intensity, the calculation unit 14c registers the calculated average value in the item “NIR” of the obtained record. Further, when the calculated average value is the average value of the intensity of the short wavelength infrared component, the calculation unit 14c registers the calculated average value in the item “SWIR” of the obtained record.
一方、検索の結果、生成部14aにより選択されたIDが「圃場ID」の項目に登録されたレコードが得られなかった場合には、算出部14cは、次のような処理を行う。例えば、算出部14cは、判別用圃場別スペクトルデータ13cに新たなレコードを1つ追加する。そして、算出部14cは、生成部14aにより選択されたIDを、追加したレコードの「圃場ID」の項目に登録する。そして、算出部14cは、算出した平均値が、緑成分の強度の平均値である場合には、追加したレコードの「G」の項目に、算出した平均値を登録する。また、算出部14cは、算出した平均値が、赤成分の強度の平均値である場合には、追加したレコードの「R」の項目に、算出した平均値を登録する。また、算出部14cは、算出した平均値が、近赤外成分の強度の平均値である場合には、追加したレコードの「NIR」の項目に、算出した平均値を登録する。また、算出部14cは、算出した平均値が、短波長赤外成分の強度の平均値である場合には、追加したレコードの「SWIR」の項目に、算出した平均値を登録する。 On the other hand, as a result of the search, when the record registered in the item “field ID” as the ID selected by the generation unit 14a is not obtained, the calculation unit 14c performs the following process. For example, the calculation unit 14c adds one new record to the discrimination-specific field-specific spectrum data 13c. Then, the calculation unit 14c registers the ID selected by the generation unit 14a in the item “Field ID” of the added record. Then, when the calculated average value is the average value of the intensity of the green component, the calculation unit 14c registers the calculated average value in the item “G” of the added record. In addition, when the calculated average value is the average value of the intensity of the red component, the calculation unit 14c registers the calculated average value in the item “R” of the added record. Further, when the calculated average value is the average value of the intensity of the near-infrared component, the calculation unit 14c registers the calculated average value in the item “NIR” of the added record. In addition, when the calculated average value is the average value of the intensity of the short wavelength infrared component, the calculation unit 14c registers the calculated average value in the item “SWIR” of the added record.
算出部14cは、上述したような方法で、判別の対象となる全ての圃場について、緑成分、赤成分、近赤外成分、短波長赤外成分のそれぞれの強度の平均値を判別用圃場別スペクトルデータ13cに登録する。このようにして、算出部14cは、作物の判別を精度良く行うことが可能なデータを算出する。 The calculation unit 14c calculates the average value of the intensity of each of the green component, the red component, the near infrared component, and the short wavelength infrared component for each field for discrimination by the method described above. Register in the spectrum data 13c. In this way, the calculation unit 14c calculates data that enables accurate crop discrimination.
また、算出部14cは、抽出された特定の波長帯域の成分の強度が複数でないと特定部14bにより判定された場合には、次の処理を行う。すなわち、算出部14cは、ヒストグラムが示す全ての強度のそれぞれに対応するピクセルにおける特定の波長帯域の成分の強度の平均値を、上述した平均値の算出方法と同様の方法で算出する。そして、平均値を算出すると、算出部14cは、生成部14aにより選択されたIDと、算出した平均値とを対応付けて上述した判別用圃場別スペクトルデータ13cに登録する。 In addition, the calculating unit 14c performs the following process when the specifying unit 14b determines that the extracted component of the specific wavelength band has multiple intensities. That is, the calculation unit 14c calculates the average value of the intensities of the components in the specific wavelength band in the pixels corresponding to all the intensities indicated by the histogram by the same method as the above-described average value calculation method. When the average value is calculated, the calculation unit 14c associates the ID selected by the generation unit 14a with the calculated average value and registers the ID in the above-described determination-specific field-specific spectrum data 13c.
判別部14dは、基準強度データ13aが示す所定の作物が植えられた地域における太陽光の反射光に含まれる特定の3種類以上の波長帯域の成分の強度に、算出部14cにより算出された3種類以上の波長帯域の成分の強度が合致するか否かを判定する。そして、判別部14dは、合致すると判定した場合には、判定対象の圃場に、所定の作物が植えられていると判別する。
The discriminating unit 14d is calculated by the calculating unit 14c based on the intensities of the components of the specific three or more wavelength bands included in the reflected light of sunlight in the area where the predetermined crop indicated by the
判別部14dの一態様について説明する。例えば、判別対象の圃場について、緑成分、赤成分、近赤外成分、短波長赤外成分のそれぞれの強度の平均値が算出部14cにより判別用圃場別スペクトルデータ13cに登録されると、判別部14dは、次の処理を行う。すなわち、判別部14dは、生成部14aにより選択されたIDが「圃場ID」の項目に登録されたレコードの「G」、「R」、「NIR」、「SWIR」の各項目に登録された強度の平均値を判別用圃場別スペクトルデータ13cから取得する。そして、判別部14dは、取得した4つの強度の平均値が、上述した「小麦についての4つの条件」を満たすか否かを判定する。4つの強度の平均値が、「小麦についての4つの条件」を満たすと判定した場合には、判別部14dは、生成部14aにより選択されたIDと、作物の種類「小麦」とを対応付けて、制御部14の内部メモリに格納する。上述したような方法で、判別部14dは、判別対象の圃場に植えられている作物が小麦であると判別する。
One aspect of the determination unit 14d will be described. For example, when the average value of the intensity of each of the green component, red component, near infrared component, and short wavelength infrared component is registered in the discriminating field-specific spectrum data 13c by the calculation unit 14c, the discrimination is performed. The unit 14d performs the following process. That is, the determination unit 14d is registered in the items “G”, “R”, “NIR”, and “SWIR” of the record in which the ID selected by the generation unit 14a is registered in the item “Field ID”. The average value of intensity is acquired from the field-specific spectral data 13c for discrimination. Then, the determination unit 14d determines whether or not the acquired average value of the four intensities satisfies the “four conditions for wheat” described above. When it is determined that the average value of the four intensities satisfies “four conditions for wheat”, the determination unit 14 d associates the ID selected by the generation unit 14 a with the crop type “wheat”. And stored in the internal memory of the
一方、取得した4つの強度の平均値が、上述した「小麦についての4つの条件」を満たさないと判定した場合には、判別部14dは、次の処理を行う。すなわち、判別部14dは、取得した4つの強度の平均値が、上述した「馬鈴薯についての4つの条件」を満たすか否かを判定する。4つの強度の平均値が、「馬鈴薯についての4つの条件」を満たすと判定した場合には、判別部14dは、生成部14aにより選択されたIDと、作物の種類「馬鈴薯」とを対応付けて、制御部14の内部メモリに格納する。上述したような方法で、判別部14dは、判別対象の圃場に植えられている作物が馬鈴薯であると判別する。
On the other hand, when it is determined that the acquired average value of the four intensities does not satisfy the above-mentioned “four conditions for wheat”, the determination unit 14d performs the following process. That is, the determination unit 14d determines whether or not the acquired average value of the four intensities satisfies the above-mentioned “four conditions for potato”. When it is determined that the average value of the four intensities satisfies “four conditions for potato”, the determination unit 14d associates the ID selected by the generation unit 14a with the crop type “potato”. And stored in the internal memory of the
一方、取得した4つの強度の平均値が、上述した「馬鈴薯についての4つの条件」を満たさないと判定した場合には、判別部14dは、次の処理を行う。すなわち、判別部14dは、取得した4つの強度の平均値が、上述した「テンサイについての4つの条件」を満たすか否かを判定する。4つの強度の平均値が、「テンサイについての4つの条件」を満たすと判定した場合には、判別部14dは、生成部14aにより選択されたIDと、作物の種類「テンサイ」とを対応付けて、制御部14の内部メモリに格納する。上述したような方法で、判別部14dは、判別対象の圃場に植えられている作物がテンサイであると判別する。
On the other hand, when it is determined that the acquired average value of the four intensities does not satisfy the above-mentioned “four conditions for potato”, the determination unit 14d performs the following process. That is, the determination unit 14d determines whether or not the acquired average value of the four intensities satisfies the “four conditions for sugar beet” described above. When it is determined that the average value of the four intensities satisfies “four conditions for sugar beet”, the determination unit 14 d associates the ID selected by the generation unit 14 a with the crop type “sugar beet”. And stored in the internal memory of the
一方、取得した4つの強度の平均値が、上述した「テンサイについての4つの条件」を満たさないと判定した場合には、判別部14dは、次の処理を行う。すなわち、判別部14dは、生成部14aにより選択されたIDと、メッセージ「小麦、馬鈴薯、テンサイのいずれの作物も植えられていません。」とを対応付けて、制御部14の内部メモリに格納する。
On the other hand, when it is determined that the acquired average value of the four intensities does not satisfy the “four conditions for sugar beet” described above, the determination unit 14d performs the following process. That is, the determination unit 14d associates the ID selected by the generation unit 14a with the message “None of wheat, potato, or sugar beet is planted” and stores it in the internal memory of the
上述したような方法で、判別部14dは、圃場に植えられている作物の種類を判定する。 With the method described above, the determination unit 14d determines the type of crop planted in the field.
出力制御部14eは、基準強度データ13aが示す所定の作物が植えられた地域における太陽光の反射光に含まれる特定の3種類以上の波長帯域の成分の強度に、算出された3種類以上の波長帯域の成分の強度が合致すると判定された場合に、次の処理を行う。すなわち、出力制御部14eは、判別対象の地域に所定の作物が植えられていることを示す情報を出力するように出力部12を制御する。
The output control unit 14e calculates three or more types of calculated intensities of components of three or more types of wavelength bands included in the reflected light of sunlight in the area where the predetermined crop indicated by the
出力制御部14eの一態様について説明する。例えば、出力制御部14eは、判別部14dにより制御部14の内部メモリにIDに対応付けられて格納された各種の情報を取得する。例えば、出力制御部14eは、内部メモリに圃場のIDと作物の種類とが対応付けられて格納されている場合には、圃場のIDと作物の種類とを取得する。また、出力制御部14eは、内部メモリに、IDに対応付けられて、メッセージ「小麦、馬鈴薯、テンサイのいずれの作物も植えられていません。」が格納されている場合には、圃場のIDとメッセージとを取得する。そして、出力制御部14eは、取得した圃場のIDと作物の種類とを対応付けて表示するように出力部12を制御する。これに加えて、出力制御部14eは、取得した圃場のIDとメッセージ「小麦、馬鈴薯、テンサイのいずれの作物も植えられていません。」とを対応付けて表示するように出力部12を制御する。
One aspect of the output control unit 14e will be described. For example, the output control unit 14e acquires various types of information stored in association with IDs in the internal memory of the
図9は、出力制御部14eによる制御によって出力部12の画面に表示される内容の一例を示す図である。図9の例に示す画面は、ID「AA」が示す圃場に小麦が植えられていることを示す。また、図9の例に示す画面は、ID「BB」が示す圃場にテンサイが植えられていることを示す。また、図9の例に示す画面は、ID「CC」が示す圃場に馬鈴薯が植えられていることを示す。また、図9の例に示す画面は、ID「EE」が示す圃場に、小麦、馬鈴薯、テンサイのいずれの作物も植えられていないことを示す。
FIG. 9 is a diagram illustrating an example of contents displayed on the screen of the
制御部14は、ASIC(Application Specific Integrated Circuit)、FPGA(Field Programmable Gate Array)、CPU(Central Processing Unit)、MPU(Micro Processing Unit)などの回路である。
The
[処理の流れ]
次に、本実施例に係る作物判別装置10が実行する処理の流れについて説明する。図10A及び図10Bは、実施例に係る作物判別処理の手順を示すフローチャートである。実施例に係る作物判別処理は、例えば、入力部11から、作物判別処理を実行するための指示が制御部14に入力された場合に、制御部14により実行される。
[Process flow]
Next, the flow of processing executed by the
図10Aに示すように、生成部14aは、圃場別スペクトルデータ13bの「圃場ID」の項目に登録された圃場を識別するためのIDを全て取得する(S101)。なお、生成部14aは、重複するIDについては、1つのIDとして取得する。そして、生成部14aは、取得したIDの中に、未選択のIDがあるか否かを判定する(S102)。未選択のIDがあると判定した場合(S102;Yes)には、生成部14aは、未選択のIDを1つ選択する(S103)。
As illustrated in FIG. 10A, the generation unit 14a acquires all IDs for identifying the fields registered in the item “Field ID” of the field-
そして、生成部14aは、圃場別スペクトルデータ13bの全てのレコードの中から、選択したIDが「圃場ID」の項目に登録された全てのレコードを特定する(S104)。そして、生成部14aは、特定した全てのレコードのそれぞれの「G」の項目に登録された緑成分の強度を全て取得する(S105)。そして、生成部14aは、取得した全ての緑成分の強度のうち、同一の緑成分の強度の数を取得した緑成分の強度ごとに計数する(S106)。ここで、生成部14aは、例えば、「N」よりも大きく、かつ、「N+10」以下の緑成分の強度については、同一の緑成分の強度「N+5」として扱う。そして、生成部14aは、横軸を緑成分の強度、縦軸をピクセルの数とするヒストグラムを生成する。すなわち、生成部14aは、選択したIDが示す圃場における緑成分の強度についてのヒストグラムを生成する(S107)。
And the production | generation part 14a specifies all the records in which the selected ID was registered into the item of "field ID" from all the records of the
そして、生成部14aは、特定した全てのレコードのそれぞれの「R」の項目に登録された赤成分の強度を全て取得する(S108)。そして、生成部14aは、取得した全ての赤成分の強度のうち、同一の赤成分の強度の数を、取得した赤成分の強度ごとに計数する(S109)。ここで、生成部14aは、例えば、「N」よりも大きく、かつ、「N+10」以下の赤成分の強度については、同一の赤成分の強度「N+5」として扱う。そして、生成部14aは、横軸を赤成分の強度、縦軸をピクセルの個数とするヒストグラムを生成する。すなわち、生成部14aは、選択したIDが示す圃場における赤成分の強度についてのヒストグラムを生成する(S110)。 Then, the generation unit 14a acquires all the intensities of the red components registered in the “R” items of all the specified records (S108). And the production | generation part 14a counts the number of the intensity | strengths of the same red component among the intensity | strengths of all the acquired red components for every intensity | strength of the acquired red component (S109). Here, for example, the intensity of the red component greater than “N” and equal to or less than “N + 10” is handled as the intensity “N + 5” of the same red component. Then, the generation unit 14a generates a histogram with the horizontal axis representing the red component intensity and the vertical axis representing the number of pixels. That is, the generation unit 14a generates a histogram about the intensity of the red component in the field indicated by the selected ID (S110).
そして、生成部14aは、特定した全てのレコードのそれぞれの「NIR」の項目に登録された近赤外成分の強度を全て取得する(S111)。そして、生成部14aは、取得した全ての近赤外成分の強度のうち、同一の近赤外成分の強度の数を、取得した近赤外成分の強度ごとに計数する(S112)。ここで、生成部14aは、例えば、「N」よりも大きく、かつ、「N+10」以下の近赤外成分の強度については、同一の近赤外成分の強度「N+5」として扱う。そして、生成部14aは、横軸を近赤外成分の強度、縦軸をピクセルの個数とするヒストグラムを生成する。すなわち、生成部14aは、選択したIDが示す圃場における近赤外成分の強度についてのヒストグラムを生成する(S113)。 Then, the generation unit 14a acquires all the intensities of the near infrared components registered in the “NIR” items of all the specified records (S111). And the production | generation part 14a counts the number of the intensity | strength of the same near-infrared component among the intensity | strengths of all the acquired near-infrared components for every intensity | strength of the acquired near-infrared component (S112). Here, for example, the intensity of the near-infrared component greater than “N” and equal to or less than “N + 10” is handled as the intensity “N + 5” of the same near-infrared component. Then, the generation unit 14a generates a histogram with the horizontal axis indicating the intensity of the near infrared component and the vertical axis indicating the number of pixels. That is, the generation unit 14a generates a histogram about the intensity of the near infrared component in the field indicated by the selected ID (S113).
そして、生成部14aは、特定した全てのレコードのそれぞれの「SWIR」の項目に登録された短波長赤外成分の強度を全て取得する(S114)。そして、生成部14aは、取得した全ての短波長赤外成分の強度のうち、同一の短波長赤外成分の強度の数を、取得した短波長赤外成分の強度ごとに計数する(S115)。ここで、生成部14aは、例えば、「N」よりも大きく、かつ、「N+10」以下の短波長赤外成分の強度については、同一の短波長成分の強度「N+5」として扱う。そして、生成部14aは、横軸を短波長赤外成分の強度、縦軸をピクセルの個数とするヒストグラムを生成する。すなわち、生成部14aは、選択したIDが示す圃場における短波長赤外成分の強度についてのヒストグラムを生成する(S116)。 Then, the generation unit 14a acquires all the intensities of the short wavelength infrared components registered in the “SWIR” items of all the specified records (S114). And the production | generation part 14a counts the number of the intensity | strength of the same short wavelength infrared component among the intensity | strengths of all the acquired short wavelength infrared components for every intensity | strength of the acquired short wavelength infrared component (S115). . Here, for example, the intensity of the short wavelength infrared component that is larger than “N” and equal to or less than “N + 10” is handled as the intensity “N + 5” of the same short wavelength component. Then, the generation unit 14a generates a histogram with the horizontal axis representing the intensity of the short wavelength infrared component and the vertical axis representing the number of pixels. That is, the generation unit 14a generates a histogram about the intensity of the short wavelength infrared component in the field indicated by the selected ID (S116).
そして、特定部14bは、生成部14aにより選択されたIDが示す圃場における、緑成分、赤成分、近赤外成分、短波長赤外成分の各成分についての4つのヒストグラムのうち、未選択のヒストグラムがあるか否かを判定する(S117)。未選択のヒストグラムがあると判定した場合(S117;Yes)には、特定部14bは、未選択のヒストグラムを1つ選択する(S118)。 And the specific | specification part 14b is unselected among four histograms about each component of a green component, a red component, a near infrared component, and a short wavelength infrared component in the field which ID selected by the production | generation part 14a shows. It is determined whether there is a histogram (S117). If it is determined that there is an unselected histogram (S117; Yes), the specifying unit 14b selects one unselected histogram (S118).
そして、特定部14bは、選択したヒストグラムにおいて、ピクセルの個数についてピークとなる特定の波長帯域の成分の強度を抽出する(S119)。ここで、ヒストグラムにおいてピークが複数ある場合には、特定部14bは、複数のピークを抽出する。 Then, the specifying unit 14b extracts the intensity of a component in a specific wavelength band that peaks for the number of pixels in the selected histogram (S119). Here, when there are a plurality of peaks in the histogram, the specifying unit 14b extracts a plurality of peaks.
そして、特定部14bは、抽出した特定の波長帯域の成分の強度が複数であるか否かを判定する(S120)。抽出した特定の波長帯域の成分の強度が複数であると判定した場合(S120;Yes)には、特定部14bは、次の処理を行う。すなわち、特定部14bは、選択したヒストグラムが示す強度とピクセルの個数との関係に基づいて、抽出した複数の特定の波長帯域の成分の強度の中から、対応するピクセルの個数が最大である特定の波長帯域の成分の強度を特定する(S121)。 And specific part 14b judges whether the intensity of the component of the extracted specific wavelength band is plurality (S120). When it determines with the intensity | strength of the component of the extracted specific wavelength band being plurality (S120; Yes), the specific | specification part 14b performs the following process. That is, the specifying unit 14b specifies that the number of corresponding pixels is the maximum from the extracted intensities of the components of the specific wavelength bands based on the relationship between the intensity indicated by the selected histogram and the number of pixels. The intensity of the component in the wavelength band is specified (S121).
そして、算出部14cは、選択されたヒストグラムが示す強度とピクセルの個数との関係において、特定部14bにより特定された強度に対応するピクセルの個数からピクセルの個数が減少する範囲内の特定の波長帯域の成分の強度を特定する(S122)。 The calculating unit 14c then selects a specific wavelength within a range in which the number of pixels decreases from the number of pixels corresponding to the intensity specified by the specifying unit 14b in the relationship between the intensity indicated by the selected histogram and the number of pixels. The intensity of the band component is specified (S122).
そして、算出部14cは、特定した全ての強度のそれぞれに対応するピクセルにおける特定の波長帯域の成分の強度の平均値を算出する(S123)。そして、算出部14cは、生成部14aにより選択されたIDと、算出した平均値とを対応付けて記憶部13に記憶された判別用圃場別スペクトルデータ13cに登録し(S124)、S117に戻る。 Then, the calculating unit 14c calculates the average value of the intensities of the components in the specific wavelength band in the pixels corresponding to all the specified intensities (S123). Then, the calculation unit 14c registers the ID selected by the generation unit 14a and the calculated average value in association with the discriminating field-specific spectrum data 13c stored in the storage unit 13 (S124), and returns to S117. .
一方、抽出された特定の波長帯域の成分の強度が複数でないと特定部14bにより判定された場合(S120;No)には、次の処理を行う。すなわち、算出部14cは、ヒストグラムが示す全ての強度のそれぞれに対応するピクセルにおける特定の波長帯域の成分の強度の平均値を、上述した平均値の算出方法と同様の方法で算出する(S125)。そして、算出部14cは、生成部14aにより選択されたIDと、算出した平均値とを対応付けて判別用圃場別スペクトルデータ13cに登録し(S126)、S117に戻る。 On the other hand, when it is determined by the specifying unit 14b that the extracted component of the specific wavelength band has a plurality of intensities (S120; No), the following processing is performed. That is, the calculating unit 14c calculates the average value of the components of the specific wavelength band in the pixels corresponding to all the intensities indicated by the histogram by the same method as the above-described average value calculating method (S125). . Then, the calculation unit 14c associates the ID selected by the generation unit 14a with the calculated average value and registers it in the discrimination-specific field-specific spectrum data 13c (S126), and returns to S117.
また、未選択のヒストグラムがないと特定部14bにより判定された場合(S117;No)には、図10Bに示すように、判別部14dは、次の処理を行う。すなわち、判別部14dは、生成部14aにより選択されたIDが「圃場ID」の項目に登録されたレコードの「G」、「R」、「NIR」、「SWIR」の各項目に登録された強度の平均値を判別用圃場別スペクトルデータ13cから取得する(S127)。そして、判別部14dは、取得した4つの強度の平均値が、「小麦についての4つの条件」を満たすか否かを判定する(S128)。4つの強度の平均値が、「小麦についての4つの条件」を満たすと判定した場合(S128;Yes)には、判別部14dは、次の処理を行う。すなわち、判別部14dは、生成部14aにより選択されたIDと、作物の種類「小麦」とを対応付けて、制御部14の内部メモリに格納し(S129)、S102に戻る。上述したような方法で、判別部14dは、判別対象の圃場に植えられている作物が小麦であると判別する。 In addition, when the specifying unit 14b determines that there is no unselected histogram (S117; No), as illustrated in FIG. 10B, the determination unit 14d performs the following process. That is, the determination unit 14d is registered in the items “G”, “R”, “NIR”, and “SWIR” of the record in which the ID selected by the generation unit 14a is registered in the item “Field ID”. The average value of intensity is acquired from the field-specific spectral data 13c for discrimination (S127). Then, the determination unit 14d determines whether or not the acquired average value of the four intensities satisfies “four conditions for wheat” (S128). When it is determined that the average value of the four intensities satisfies “four conditions for wheat” (S128; Yes), the determination unit 14d performs the following process. That is, the determination unit 14d associates the ID selected by the generation unit 14a with the crop type “wheat” and stores them in the internal memory of the control unit 14 (S129), and returns to S102. By the method as described above, the determination unit 14d determines that the crop planted in the determination target field is wheat.
一方、取得した4つの強度の平均値が、「小麦についての4つの条件」を満たさないと判定した場合(S128;No)には、判別部14dは、次の処理を行う。すなわち、判別部14dは、取得した4つの強度の平均値が、「馬鈴薯についての4つの条件」を満たすか否かを判定する(S130)。4つの強度の平均値が、「馬鈴薯についての4つの条件」を満たすと判定した場合(S130;Yes)には、判別部14dは、次の処理を行う。すなわち、判別部14dは、生成部14aにより選択されたIDと、作物の種類「馬鈴薯」とを対応付けて、制御部14の内部メモリに格納し(S131)、S102に戻る。上述したような方法で、判別部14dは、判別対象の圃場に植えられている作物が馬鈴薯であると判別する。 On the other hand, when it is determined that the acquired average value of the four intensities does not satisfy the “four conditions for wheat” (S128; No), the determination unit 14d performs the following process. That is, the determination unit 14d determines whether or not the acquired average value of the four intensities satisfies “four conditions for potato” (S130). When it is determined that the average value of the four intensities satisfies “four conditions for potato” (S130; Yes), the determination unit 14d performs the following process. That is, the determination unit 14d associates the ID selected by the generation unit 14a with the crop type “potato” and stores it in the internal memory of the control unit 14 (S131), and returns to S102. By the method as described above, the determination unit 14d determines that the crop planted in the determination target field is a potato.
一方、取得した4つの強度の平均値が、上述した「馬鈴薯についての4つの条件」を満たさないと判定した場合(S130;No)には、判別部14dは、次の処理を行う。すなわち、判別部14dは、取得した4つの強度の平均値が、「テンサイについての4つの条件」を満たすか否かを判定する(S132)。4つの強度の平均値が、「テンサイについての4つの条件」を満たすと判定した場合(S132;Yes)には、判別部14dは、次の処理を行う。すなわち、判別部14dは、生成部14aにより選択されたIDと、作物の種類「テンサイ」とを対応付けて、制御部14の内部メモリに格納し(S133)、S102に戻る。上述したような方法で、判別部14dは、判別対象の圃場に植えられている作物がテンサイであると判別する。 On the other hand, if it is determined that the acquired average value of the four intensities does not satisfy the above-mentioned “four conditions for potato” (S130; No), the determination unit 14d performs the following process. That is, the determination unit 14d determines whether or not the acquired average value of the four strengths satisfies “four conditions for sugar beet” (S132). When it is determined that the average value of the four intensities satisfies “four conditions for sugar beet” (S132; Yes), the determination unit 14d performs the following process. That is, the determination unit 14d associates the ID selected by the generation unit 14a with the crop type “sugar beet” and stores them in the internal memory of the control unit 14 (S133), and returns to S102. With the method described above, the determination unit 14d determines that the crop planted in the field to be determined is a sugar beet.
一方、取得した4つの強度の平均値が、「テンサイについての4つの条件」を満たさないと判定した場合(S132;No)には、判別部14dは、次の処理を行う。すなわち、判別部14dは、生成部14aにより選択されたIDと、メッセージ「小麦、馬鈴薯、テンサイのいずれの作物も植えられていません。」とを対応付けて、制御部14の内部メモリに格納し(S134)、S102に戻る。
On the other hand, when it is determined that the acquired average value of the four intensities does not satisfy “four conditions for sugar beet” (S132; No), the determination unit 14d performs the following process. That is, the determination unit 14d associates the ID selected by the generation unit 14a with the message “None of wheat, potato, or sugar beet is planted” and stores it in the internal memory of the
また、未選択のIDがないと生成部14aにより判定された場合(S102;No)には、出力制御部14eは、判別部14dにより制御部14の内部メモリにIDに対応付けられて格納された各種の情報を取得する。例えば、出力制御部14eは、内部メモリに圃場のIDと作物の種類とが対応付けられて格納されている場合には、圃場のIDと作物の種類とを取得する。また、出力制御部14eは、内部メモリに、IDに対応付けられて、メッセージ「小麦、馬鈴薯、テンサイのいずれの作物も植えられていません。」が格納されている場合には、圃場のIDとメッセージとを取得する(S135)。そして、出力制御部14eは、取得した圃場のIDと作物の種類とを対応付けて表示するように出力部12を制御する。これに加えて、出力制御部14eは、取得した圃場のIDとメッセージ「小麦、馬鈴薯、テンサイのいずれの作物も植えられていません。」とを対応付けて表示するように出力部12を制御し(S136)、作物判別処理を終了する。
If the generation unit 14a determines that there is no unselected ID (S102; No), the output control unit 14e is stored in the internal memory of the
上述してきたように、本実施例に係る作物判別装置10は、圃場別スペクトルデータ13bに基づいて、特定の波長帯域ごとに、特定の波長帯域の成分の強度と、強度に対応するピクセルの個数との関係を示す情報であるヒストグラムを生成する。そして、作物判別装置10は、生成されたヒストグラムに基づいて、特定の波長帯域ごとに、ピクセルの個数が最大である強度を特定する。そして、作物判別装置10は、ヒストグラムに基づいて、特定の波長帯域ごとに、次の処理を行う。すなわち、作物判別装置10は、特定された強度を含む、最大のピクセルの個数からピクセルの個数が減少する範囲内の特定の波長帯域の成分の強度に対応するピクセルにおける特定の波長帯域の成分の強度の平均値を算出する。そして、作物判別装置10は、特定の波長帯域ごとに算出された平均値に基づいて、判別対象の圃場に所定の作物(小麦、馬鈴薯、テンサイ)が植えられているか否かを判別する。このように、作物判別装置10は、判別対象の圃場のピクセルの一部に、圃場以外の領域が少なくとも一部に含まれてしまう場合であっても、次のような処理を行う。すなわち、作物判別装置10は、判別対象の圃場を少なくとも一部に含むピクセルのうち、判別対象の圃場の領域のみを含むピクセルにおける強度の平均値を算出する。そして、作物判別装置10は、算出した平均値に基づいて、判別対象の圃場に所定の作物が植えられているか否かを判別する。このように、作物判別装置10によれば、精度良く作物の判別を行うことが可能な、判別対象の圃場の領域のみを含むピクセルにおける強度の平均値に基づいて、作物の判別を行うので、精度良く作物を判別することができる。したがって、作物判別装置10によれば、作物の判別の精度の低下を抑制することができる。
As described above, the
また、作物判別装置10は、判別対象の圃場のピクセルの一部に、圃場以外の領域が少なくとも一部に含まれてしまう場合であっても、次のような処理を行う。すなわち、作物判別装置10は、判別対象の圃場のピクセルのうち外周部分のピクセルを全て除くことなく、圃場以外の領域の一部が少なくとも含まれるピクセルのみを除いた上で、強度の平均値を算出する。それゆえ、作物判別装置10によれば、外周部分のピクセルを全て除いて強度の平均値を算出した場合と比較して、強度の平均値を算出する際の強度の数が多いため、判定の精度が良い強度の平均値を算出することができる。
In addition, the
ここで、実施例1に係る作物判別装置10が、100個の圃場に対して上述した作物判別処理による作物の判別を行う実験を行ったところ、100個の圃場全てに対して、植えられている作物を正しく判別することができた。一方、判別対象の圃場のピクセルの一部に判別対象の圃場以外の領域が含まれてしまう場合であっても、判別対象の圃場の全ピクセルにおける強度の平均値を算出する従来技術では、100個の圃場のうち、25個の圃場に植えられた作物を正しく判別できなかった。すなわち、従来技術による作物の判別の精度は、75%の精度である。
Here, when the
また、「小麦についての4つの条件」において、1番目の条件、2番目の条件及び4番目の条件の3つの条件を満たしているか否かを判定することにより、小麦が植えられているか否かを判定する実験を行った。また、「馬鈴薯についての4つの条件」において、1番目の条件、2番目の条件及び4番目の条件の3つの条件を満たしているか否かを判定することにより、馬鈴薯が植えられているか否かを判定する実験を行った。また、「テンサイについての4つの条件」において、1番目の条件、2番目の条件及び4番目の条件の3つの条件を満たしているか否かを判定することにより、テンサイが植えられているか否かを判定する実験を行った。これらの実験の結果、93%の圃場に対して、植えられている作物を正しく判別することができ、精度良く作物を判別できた。そこで、作物判別装置10は、「小麦についての4つの条件」、「馬鈴薯についての4つの条件」及び「テンサイについての4つの条件」のそれぞれの1番目、2番目及び4番目の条件の3つの条件を判定して植えられている作物を判別するようにしてもよい。
Whether or not wheat is planted by determining whether or not three conditions of the first condition, the second condition, and the fourth condition are satisfied in the “four conditions for wheat” An experiment was conducted to determine the above. Whether or not potatoes are planted by determining whether or not the three conditions of the first condition, the second condition, and the fourth condition are satisfied in the “four conditions for potatoes” An experiment was conducted to determine the above. Whether or not sugar beet is planted by determining whether or not the three conditions of the first condition, the second condition, and the fourth condition are satisfied in the “four conditions for sugar beet”. An experiment was conducted to determine the above. As a result of these experiments, it was possible to correctly discriminate the planted crops with respect to 93% of the fields and to discriminate the crops with high accuracy. Therefore, the
また、「小麦についての4つの条件」において、1番目の条件、2番目の条件及び3番目の条件の3つの条件を満たしているか否かを判定することにより、小麦が植えられているか否かを判定する実験を行った。また、「馬鈴薯についての4つの条件」において、1番目の条件、2番目の条件及び3番目の条件の3つの条件を満たしているか否かを判定することにより、馬鈴薯が植えられているか否かを判定する実験を行った。また、「テンサイについての4つの条件」において、1番目の条件、2番目の条件及び3番目の条件の3つの条件を満たしているか否かを判定することにより、テンサイが植えられているか否かを判定する実験を行った。これらの実験の結果、97%の圃場に対して、植えられている作物を正しく判別することができ、精度良く作物を判別できた。そこで、作物判別装置10は、「小麦についての4つの条件」、「馬鈴薯についての4つの条件」及び「テンサイについての4つの条件」のそれぞれの1番目、2番目及び3番目の条件の3つの条件を判定して植えられている作物を判別するようにしてもよい。
Whether or not wheat is planted by determining whether or not the three conditions of the first condition, the second condition, and the third condition are satisfied in the “four conditions for wheat” An experiment was conducted to determine the above. Whether or not a potato is planted by determining whether or not the three conditions of the first condition, the second condition, and the third condition are satisfied in the “four conditions for potato” An experiment was conducted to determine the above. Whether or not sugar beet is planted by determining whether or not the three conditions of the first condition, the second condition, and the third condition are satisfied in the “four conditions for sugar beet” An experiment was conducted to determine the above. As a result of these experiments, it was possible to correctly discriminate the planted crops with respect to 97% of the fields and to discriminate the crops with high accuracy. Therefore, the
また、「小麦についての4つの条件」において、1番目の条件、3番目の条件及び4番目の条件の3つの条件を満たしているか否かを判定することにより、小麦が植えられているか否かを判定する実験を行った。また、「馬鈴薯についての4つの条件」において、1番目の条件、3番目の条件及び4番目の条件の3つの条件を満たしているか否かを判定することにより、馬鈴薯が植えられているか否かを判定する実験を行った。また、「テンサイについての4つの条件」において、1番目の条件、3番目の条件及び4番目の条件の3つの条件を満たしているか否かを判定することにより、テンサイが植えられているか否かを判定する実験を行った。これらの実験の結果、95%の圃場に対して、植えられている作物を正しく判別することができ、精度良く作物を判別できた。そこで、作物判別装置10は、「小麦についての4つの条件」、「馬鈴薯についての4つの条件」及び「テンサイについての4つの条件」のそれぞれの1番目、3番目及び4番目の条件の3つの条件を判定して植えられている作物を判別するようにしてもよい。
Whether or not wheat is planted by determining whether the three conditions of the first condition, the third condition, and the fourth condition are satisfied in the “four conditions for wheat” An experiment was conducted to determine the above. Whether or not a potato is planted by determining whether or not three conditions of the first condition, the third condition, and the fourth condition are satisfied in the “four conditions for potato” An experiment was conducted to determine the above. Whether or not sugar beet is planted by determining whether or not the three conditions of the first condition, the third condition, and the fourth condition are satisfied in the “four conditions for sugar beet”. An experiment was conducted to determine the above. As a result of these experiments, it was possible to correctly discriminate the planted crops with respect to 95% of the fields and to discriminate the crops with high accuracy. Therefore, the
また、「小麦についての4つの条件」において、2番目の条件、3番目の条件及び4番目の条件の3つの条件を満たしているか否かを判定することにより、小麦が植えられているか否かを判定する実験を行った。また、「馬鈴薯についての4つの条件」において、2番目の条件、3番目の条件及び4番目の条件の3つの条件を満たしているか否かを判定することにより、馬鈴薯が植えられているか否かを判定する実験を行った。また、「テンサイについての4つの条件」において、2番目の条件、3番目の条件及び4番目の条件の3つの条件を満たしているか否かを判定することにより、テンサイが植えられているか否かを判定する実験を行った。これらの実験の結果、96%の圃場に対して、植えられている作物を正しく判別することができ、精度良く作物を判別できた。そこで、作物判別装置10は、「小麦についての4つの条件」、「馬鈴薯についての4つの条件」及び「テンサイについての4つの条件」のそれぞれの2番目、3番目及び4番目の条件の3つの条件を判定して植えられている作物を判別するようにしてもよい。
Whether or not wheat is planted by determining whether or not three conditions of the second condition, the third condition, and the fourth condition are satisfied in the “four conditions for wheat” An experiment was conducted to determine the above. Whether or not a potato is planted by determining whether or not three conditions of the second condition, the third condition, and the fourth condition are satisfied in the “four conditions for potato” An experiment was conducted to determine the above. Whether or not sugar beet is planted by determining whether or not the three conditions of the second condition, the third condition, and the fourth condition are satisfied in the “four conditions for sugar beet” An experiment was conducted to determine the above. As a result of these experiments, it was possible to correctly discriminate the planted crops with respect to 96% of the fields and to discriminate the crops with high accuracy. Therefore, the
また、「小麦についての4つの条件」、「馬鈴薯についての4つの条件」及び「テンサイについての4つの条件」のそれぞれの1番目及び2番目の条件の2つの条件を判定して、植えられている作物を判別することができるか否かの実験を行った。しかしながら、この実験では、24%の正解率しか得られなかった。 In addition, the two conditions of the first and second conditions of “four conditions for wheat”, “four conditions for potato” and “four conditions for sugar beet”, respectively, are judged and planted. An experiment was conducted to determine whether or not crops could be identified. However, in this experiment, only a correct answer rate of 24% was obtained.
また、「小麦についての4つの条件」、「馬鈴薯についての4つの条件」及び「テンサイについての4つの条件」のそれぞれの1番目及び3番目の条件の2つの条件を判定して、植えられている作物を判別することができるか否かの実験を行った。しかしながら、この実験では、25%の正解率しか得られなかった。 Also, the two conditions of the first and third conditions of “four conditions for wheat”, “four conditions for potato” and “four conditions for sugar beet”, respectively, are judged and planted. An experiment was conducted to determine whether or not crops could be identified. However, in this experiment, only a correct answer rate of 25% was obtained.
また、「小麦についての4つの条件」、「馬鈴薯についての4つの条件」及び「テンサイについての4つの条件」のそれぞれの1番目及び4番目の条件の2つの条件を判定して、植えられている作物を判別することができるか否かの実験を行った。しかしながら、この実験では、20%の正解率しか得られなかった。 In addition, the two conditions of the first and fourth conditions of “four conditions for wheat”, “four conditions for potato” and “four conditions for sugar beet”, respectively, are determined and planted. An experiment was conducted to determine whether or not crops could be identified. However, in this experiment, only a correct answer rate of 20% was obtained.
また、「小麦についての4つの条件」、「馬鈴薯についての4つの条件」及び「テンサイについての4つの条件」のそれぞれの2番目及び3番目の条件の2つの条件を判定して、植えられている作物を判別することができるか否かの実験を行った。しかしながら、この実験では、16%の正解率しか得られなかった。 In addition, the two conditions of the second and third conditions of “four conditions for wheat”, “four conditions for potato” and “four conditions for sugar beet”, respectively, are judged and planted. An experiment was conducted to determine whether or not crops could be identified. However, in this experiment, only a correct answer rate of 16% was obtained.
また、「小麦についての4つの条件」、「馬鈴薯についての4つの条件」及び「テンサイについての4つの条件」のそれぞれの2番目及び4番目の条件の2つの条件を判定して、植えられている作物を判定することができるか否かの実験を行った。しかしながら、この実験では、14%の正解率しか得られなかった。 Also, the two conditions of the second and fourth conditions of “four conditions for wheat”, “four conditions for potato” and “four conditions for sugar beet”, respectively, are judged and planted. An experiment was conducted to determine whether or not crops can be judged. However, in this experiment, only a correct answer rate of 14% was obtained.
また、「小麦についての4つの条件」、「馬鈴薯についての4つの条件」及び「テンサイについての4つの条件」のそれぞれの3番目及び4番目の条件の2つの条件を判定して、植えられている作物を判定することができるか否かの実験を行った。しかしながら、この実験では、12%の正解率しか得られなかった。 In addition, the three conditions of the “four conditions for wheat”, “four conditions for potato” and “four conditions for sugar beet”, and the third condition, respectively, are judged and planted. An experiment was conducted to determine whether or not crops can be judged. However, in this experiment, only a correct answer rate of 12% was obtained.
また、「小麦」についての「G」、「R」、「NIR」、「SWIR」のそれぞれの範囲の最小値を、「小麦」についての「G」、「R」、「NIR」、「SWIR」のそれぞれを0.9倍した値とした。また、「小麦」についての「G」、「R」、「NIR」、「SWIR」のそれぞれの範囲の最大値を、「小麦」についての「G」、「R」、「NIR」、「SWIR」のそれぞれを1.1倍した値とした。また、「馬鈴薯」についての「G」、「R」、「NIR」、「SWIR」のそれぞれの範囲の最小値を、「馬鈴薯」についての「G」、「R」、「NIR」、「SWIR」のそれぞれを0.9倍した値とした。また、「馬鈴薯」についての「G」、「R」、「NIR」、「SWIR」のそれぞれの範囲の最大値を、「馬鈴薯」についての「G」、「R」、「NIR」、「SWIR」のそれぞれを1.1倍した値とした。また、「テンサイ」についての「G」、「R」、「NIR」、「SWIR」のそれぞれの範囲の最小値を、「テンサイ」についての「G」、「R」、「NIR」、「SWIR」のそれぞれを0.9倍した値とした。また、「テンサイ」についての「G」、「R」、「NIR」、「SWIR」のそれぞれの範囲の最大値を、「テンサイ」についての「G」、「R」、「NIR」、「SWIR」のそれぞれを1.1倍した値とした。そして、上述した実施例と同様の方法で、植えられている作物の判定を行う実験を行ったところ、97%の正解率が得られた。そこで、作物判定装置10は、「小麦」、「馬鈴薯」、「テンサイ」のそれぞれについての「G」、「R」、「NIR」、「SWIR」の範囲を上述した範囲とすることができる。
Further, the minimum value of each range of “G”, “R”, “NIR”, and “SWIR” for “wheat” is set to “G”, “R”, “NIR”, “SWIR” for “wheat”. Each of the values was multiplied by 0.9. In addition, the maximum value of each range of “G”, “R”, “NIR”, “SWIR” for “wheat” is set to “G”, “R”, “NIR”, “SWIR” for “wheat”. Each value of "" was multiplied by 1.1. Further, the minimum value of each range of “G”, “R”, “NIR”, “SWIR” for “potato” is set to “G”, “R”, “NIR”, “SWIR” for “potato”. Each of the values was multiplied by 0.9. In addition, the maximum value of each range of “G”, “R”, “NIR”, “SWIR” for “potato” is set to “G”, “R”, “NIR”, “SWIR” for “potato”. Each value of "" was multiplied by 1.1. Also, the minimum value of each range of “G”, “R”, “NIR”, “SWIR” for “sugar beet” is set to “G”, “R”, “NIR”, “SWIR” for “sugar beet”. Each of the values was multiplied by 0.9. In addition, the maximum value of each range of “G”, “R”, “NIR”, “SWIR” for “sugar beet” is set to “G”, “R”, “NIR”, “SWIR” for “sugar beet”. Each value of "" was multiplied by 1.1. And when the experiment which determines the planted plant by the method similar to the Example mentioned above was conducted, the correct answer rate of 97% was obtained. Therefore, the
また、「小麦」についての「G」、「R」、「NIR」、「SWIR」のそれぞれの範囲の最小値を、「小麦」についての「G」、「R」、「NIR」、「SWIR」のそれぞれを0.7倍した値とした。また、「小麦」についての「G」、「R」、「NIR」、「SWIR」のそれぞれの範囲の最大値を、「小麦」についての「G」、「R」、「NIR」、「SWIR」のそれぞれを1.3倍した値とした。また、「馬鈴薯」についての「G」、「R」、「NIR」、「SWIR」のそれぞれの範囲の最小値を、「馬鈴薯」についての「G」、「R」、「NIR」、「SWIR」のそれぞれを0.7倍した値とした。また、「馬鈴薯」についての「G」、「R」、「NIR」、「SWIR」のそれぞれの範囲の最大値を、「馬鈴薯」についての「G」、「R」、「NIR」、「SWIR」のそれぞれを1.3倍した値とした。また、「テンサイ」についての「G」、「R」、「NIR」、「SWIR」のそれぞれの範囲の最小値を、「テンサイ」についての「G」、「R」、「NIR」、「SWIR」のそれぞれを0.7倍した値とした。また、「テンサイ」についての「G」、「R」、「NIR」、「SWIR」のそれぞれの範囲の最大値を、「テンサイ」についての「G」、「R」、「NIR」、「SWIR」のそれぞれを1.3倍した値とした。そして、上述した実施例と同様の方法で、植えられている作物の判定を行う実験を行ったところ、95%の正解率が得られた。そこで、作物判定装置10は、「小麦」、「馬鈴薯」、「テンサイ」のそれぞれについての「G」、「R」、「NIR」、「SWIR」の範囲を上述した範囲とすることができる。
Further, the minimum value of each range of “G”, “R”, “NIR”, and “SWIR” for “wheat” is set to “G”, “R”, “NIR”, “SWIR” for “wheat”. Each value of “is a value multiplied by 0.7. In addition, the maximum value of each range of “G”, “R”, “NIR”, “SWIR” for “wheat” is set to “G”, “R”, “NIR”, “SWIR” for “wheat”. Each of the values was multiplied by 1.3. Further, the minimum value of each range of “G”, “R”, “NIR”, “SWIR” for “potato” is set to “G”, “R”, “NIR”, “SWIR” for “potato”. Each value of “is a value multiplied by 0.7. In addition, the maximum value of each range of “G”, “R”, “NIR”, “SWIR” for “potato” is set to “G”, “R”, “NIR”, “SWIR” for “potato”. Each of the values was multiplied by 1.3. Also, the minimum value of each range of “G”, “R”, “NIR”, “SWIR” for “sugar beet” is set to “G”, “R”, “NIR”, “SWIR” for “sugar beet”. Each value of “is a value multiplied by 0.7. In addition, the maximum value of each range of “G”, “R”, “NIR”, “SWIR” for “sugar beet” is set to “G”, “R”, “NIR”, “SWIR” for “sugar beet”. Each of the values was multiplied by 1.3. And when the experiment which judges the planted plant by the method similar to the Example mentioned above was conducted, the correct answer rate of 95% was obtained. Therefore, the
また、「小麦」についての「G」、「R」、「NIR」、「SWIR」のそれぞれの範囲の最小値を、「小麦」についての「G」、「R」、「NIR」、「SWIR」のそれぞれを0.95倍した値とした。また、「小麦」についての「G」、「R」、「NIR」、「SWIR」のそれぞれの範囲の最大値を、「小麦」についての「G」、「R」、「NIR」、「SWIR」のそれぞれを1.05倍した値とした。また、「馬鈴薯」についての「G」、「R」、「NIR」、「SWIR」のそれぞれの範囲の最小値を、「馬鈴薯」についての「G」、「R」、「NIR」、「SWIR」のそれぞれを0.95倍した値とした。また、「馬鈴薯」についての「G」、「R」、「NIR」、「SWIR」のそれぞれの範囲の最大値を、「馬鈴薯」についての「G」、「R」、「NIR」、「SWIR」のそれぞれを1.05倍した値とした。また、「テンサイ」についての「G」、「R」、「NIR」、「SWIR」のそれぞれの範囲の最小値を、「テンサイ」についての「G」、「R」、「NIR」、「SWIR」のそれぞれを0.95倍した値とした。また、「テンサイ」についての「G」、「R」、「NIR」、「SWIR」のそれぞれの範囲の最大値を、「テンサイ」についての「G」、「R」、「NIR」、「SWIR」のそれぞれを1.05倍した値とした。そして、上述した実施例と同様の方法で、植えられている作物の判定を行う実験を行ったところ、55%の正解率しか得られなかった。 Further, the minimum value of each range of “G”, “R”, “NIR”, and “SWIR” for “wheat” is set to “G”, “R”, “NIR”, “SWIR” for “wheat”. Each of the above values was multiplied by 0.95. In addition, the maximum value of each range of “G”, “R”, “NIR”, “SWIR” for “wheat” is set to “G”, “R”, “NIR”, “SWIR” for “wheat”. Each of the values was multiplied by 1.05. Further, the minimum value of each range of “G”, “R”, “NIR”, “SWIR” for “potato” is set to “G”, “R”, “NIR”, “SWIR” for “potato”. Each of the above values was multiplied by 0.95. In addition, the maximum value of each range of “G”, “R”, “NIR”, “SWIR” for “potato” is set to “G”, “R”, “NIR”, “SWIR” for “potato”. Each of the values was multiplied by 1.05. Also, the minimum value of each range of “G”, “R”, “NIR”, “SWIR” for “sugar beet” is set to “G”, “R”, “NIR”, “SWIR” for “sugar beet”. Each of the above values was multiplied by 0.95. In addition, the maximum value of each range of “G”, “R”, “NIR”, “SWIR” for “sugar beet” is set to “G”, “R”, “NIR”, “SWIR” for “sugar beet”. Each of the values was multiplied by 1.05. And when the experiment which judges the planted plant by the method similar to the Example mentioned above was conducted, only the correct answer rate of 55% was obtained.
また、「小麦」についての「G」、「R」、「NIR」、「SWIR」のそれぞれの範囲の最小値を、「小麦」についての「G」、「R」、「NIR」、「SWIR」のそれぞれを0.65倍した値とした。また、「小麦」についての「G」、「R」、「NIR」、「SWIR」のそれぞれの範囲の最大値を、「小麦」についての「G」、「R」、「NIR」、「SWIR」のそれぞれを1.35倍した値とした。また、「馬鈴薯」についての「G」、「R」、「NIR」、「SWIR」のそれぞれの範囲の最小値を、「馬鈴薯」についての「G」、「R」、「NIR」、「SWIR」のそれぞれを0.65倍した値とした。また、「馬鈴薯」についての「G」、「R」、「NIR」、「SWIR」のそれぞれの範囲の最大値を、「馬鈴薯」についての「G」、「R」、「NIR」、「SWIR」のそれぞれを1.35倍した値とした。また、「テンサイ」についての「G」、「R」、「NIR」、「SWIR」のそれぞれの範囲の最小値を、「テンサイ」についての「G」、「R」、「NIR」、「SWIR」のそれぞれを0.65倍した値とした。また、「テンサイ」についての「G」、「R」、「NIR」、「SWIR」のそれぞれの範囲の最大値を、「テンサイ」についての「G」、「R」、「NIR」、「SWIR」のそれぞれを1.35倍した値とした。そして、上述した実施例と同様の方法で、植えられている作物の判定を行う実験を行ったところ、52%の正解率しか得られなかった。 Further, the minimum value of each range of “G”, “R”, “NIR”, and “SWIR” for “wheat” is set to “G”, “R”, “NIR”, “SWIR” for “wheat”. Each value of “is a value multiplied by 0.65. In addition, the maximum value of each range of “G”, “R”, “NIR”, “SWIR” for “wheat” is set to “G”, “R”, “NIR”, “SWIR” for “wheat”. Each of the values was multiplied by 1.35. Further, the minimum value of each range of “G”, “R”, “NIR”, “SWIR” for “potato” is set to “G”, “R”, “NIR”, “SWIR” for “potato”. Each value of “is a value multiplied by 0.65. In addition, the maximum value of each range of “G”, “R”, “NIR”, “SWIR” for “potato” is set to “G”, “R”, “NIR”, “SWIR” for “potato”. Each of the values was multiplied by 1.35. Also, the minimum value of each range of “G”, “R”, “NIR”, “SWIR” for “sugar beet” is set to “G”, “R”, “NIR”, “SWIR” for “sugar beet”. Each value of “is a value multiplied by 0.65. In addition, the maximum value of each range of “G”, “R”, “NIR”, “SWIR” for “sugar beet” is set to “G”, “R”, “NIR”, “SWIR” for “sugar beet”. Each of the values was multiplied by 1.35. And when the experiment which determines the planted plant by the method similar to the Example mentioned above was conducted, only the correct answer rate of 52% was obtained.
さて、これまで開示の装置に関する実施例について説明したが、本発明は上述した実施例以外にも、種々の異なる形態にて実施されてよいものである。 Although the embodiments related to the disclosed apparatus have been described above, the present invention may be implemented in various different forms other than the above-described embodiments.
また、各実施例において説明した各処理のうち、自動的に行われるものとして説明した処理の全部または一部を手動的に行うこともできる。また、各実施例において説明した各処理のうち、手動的に行われるものとして説明した処理の全部または一部を公知の方法で自動的に行うこともできる。 In addition, among the processes described in the embodiments, all or a part of the processes described as being automatically performed can be manually performed. In addition, among the processes described in the embodiments, all or a part of the processes described as being performed manually can be automatically performed by a known method.
また、各種の負荷や使用状況などに応じて、各実施例において説明した各処理の各ステップでの処理を任意に細かくわけたり、あるいはまとめたりすることができる。また、ステップを省略することもできる。 In addition, the processing at each step of each processing described in each embodiment can be arbitrarily finely divided or combined according to various loads and usage conditions. Also, the steps can be omitted.
また、各種の負荷や使用状況などに応じて、各実施例において説明した各処理の各ステップでの処理の順番を変更できる。 Further, the order of processing at each step of each processing described in each embodiment can be changed according to various loads and usage conditions.
また、図示した各装置の各構成要素は機能概念的なものであり、必ずしも物理的に図示の如く構成されていることを要しない。すなわち、各装置の分散・統合の具体的状態は図示のものに限られず、その全部または一部を、各種の負荷や使用状況などに応じて、任意の単位で機能的または物理的に分散・統合して構成することができる。 Further, each component of each illustrated apparatus is functionally conceptual, and does not necessarily need to be physically configured as illustrated. In other words, the specific state of distribution / integration of each device is not limited to the one shown in the figure, and all or a part thereof may be functionally or physically distributed or arbitrarily distributed in arbitrary units according to various loads or usage conditions. Can be integrated and configured.
[作物判別プログラム]
また、上記の各実施例で説明した作物判別装置10の各種の処理は、あらかじめ用意されたプログラムをパーソナルコンピュータやワークステーションなどのコンピュータシステムで実行することによって実現することもできる。そこで、以下では、図11を用いて、上記の実施例で説明した作物判別装置10と同様の機能を有する作物判別プログラムを実行するコンピュータの一例を説明する。図11は、作物判別プログラムを実行するコンピュータを示す図である。
[Crop identification program]
Moreover, the various processes of the
図11に示すように、コンピュータ300は、CPU310、ROM320、HDD(Hard Disk Drive)330、RAM340、入力装置350、出力装置360を有する。これら各機器310〜360は、バス370を介して接続されている。
As illustrated in FIG. 11, the
入力装置350は、キーボードやマウスなどである。入力装置350は、入力部11に対応する。出力装置360は、液晶ディスプレイなどである。出力装置360は、出力部12に対応する。
The
ROM320には、OS(Operating System)などの基本プログラムが記憶されている。また、HDD330には、上記の各実施例で示す生成部14a、特定部14b、算出部14c、判別部14d、出力制御部14eと同様の機能を発揮する作物判別プログラム330aが予め記憶される。また、HDD330には、記憶部13に記憶された各種のデータが設けられる。
The
そして、CPU310が、作物判別プログラム330aをHDD330から読み出して実行する。
Then, the
そして、CPU310は、HDD330から各種のデータを読み出してRAM340に格納する。さらに、CPU310は、RAM340に格納された各種のデータを用いて、作物判別プログラム330aを実行する。なお、RAM340に格納されるデータは、常に全てのデータがRAM340に格納されなくともよい。処理に用いられるデータがRAM340に格納されれば良い。
Then, the
なお、上記した作物判別プログラム330aについては、必ずしも最初からHDD330に記憶させなくともよい。
Note that the
例えば、コンピュータ300に挿入されるフレキシブルディスク(FD)、CD−ROM、DVDディスク、光磁気ディスク、ICカードなどの「可搬用の物理媒体」にプログラムを記憶させておく。そして、コンピュータ300がこれらからプログラムを読み出して実行するようにしてもよい。
For example, the program is stored in a “portable physical medium” such as a flexible disk (FD), a CD-ROM, a DVD disk, a magneto-optical disk, or an IC card inserted into the
さらには、公衆回線、インターネット、LAN、WANなどを介してコンピュータ300に接続される「他のコンピュータ(またはサーバ)」などにプログラムを記憶させておく。そして、コンピュータ300がこれらからプログラムを読み出して実行するようにしてもよい。
Furthermore, the program is stored in “another computer (or server)” connected to the
10 作物判別装置
13 記憶部
13a 基準強度データ
13b 圃場別スペクトルデータ
13c 判別用圃場別スペクトルデータ
14 制御部
14a 生成部
14b 特定部
14c 算出部
14d 判別部
14e 出力制御部
DESCRIPTION OF
Claims (6)
前記記憶部に記憶された記憶内容に基づいて、前記特定の波長帯域ごとに、前記特定の波長帯域の成分の強度と、該強度に対応するピクセルの個数との関係を示す情報を生成する生成部と、
前記生成部により生成された情報に基づいて、前記特定の波長帯域ごとに、ピクセルの個数が最大である強度を特定する特定部と、
前記生成部により生成された情報に基づいて、前記特定の波長帯域ごとに、前記特定部により特定された強度を含む、前記最大のピクセルの個数からピクセルの個数が減少する範囲内の前記特定の波長帯域の成分の強度に対応するピクセルにおける前記特定の波長帯域の成分の強度の平均値を算出する算出部と、
前記算出部により前記特定の波長帯域ごとに算出された平均値に基づいて、前記ある地域に所定の作物が植えられているか否かを判別する判別部と、
を有することを特徴とする作物判別装置。 Some regions were obtained by imaging by the camera, a storage unit for storing the intensity of the components of a plurality of specific wavelength band in each of the plurality of pixels including at least a part of the area of the region where the certain,
Generation that generates information indicating the relationship between the intensity of the component in the specific wavelength band and the number of pixels corresponding to the intensity for each specific wavelength band based on the storage content stored in the storage unit And
Based on the information generated by the generating unit, for each of the specific wavelength bands, a specifying unit that specifies the intensity with the maximum number of pixels,
Based on the information generated by the generation unit, the specific wavelength within a range in which the number of pixels decreases from the maximum number of pixels, including the intensity specified by the specification unit, for each specific wavelength band. A calculation unit that calculates an average value of the intensity of the component of the specific wavelength band in the pixel corresponding to the intensity of the component of the wavelength band;
Based on the average value calculated for each of the specific wavelength band by the calculation unit, a determination unit for determining whether the region to a Jo Tokoro crop is planted in the,
A crop discriminating apparatus comprising:
前記判別部は、前記算出部により前記特定の波長帯域ごとに算出された平均値と、前記記憶部に記憶された所定の作物が植えられた地域における太陽光の反射光に含まれる前記複数の特定の波長帯域の成分の強度とが合致する場合に、前記ある地域に前記所定の作物が植えられていると判別する
ことを特徴とする請求項1に記載の作物判別装置。 The storage unit stores intensities of components of the plurality of specific wavelength bands included in reflected sunlight of an area where the predetermined crop is planted,
The determination unit includes the average value calculated for each of the specific wavelength bands by the calculation unit, and the plurality of sunlight included in reflected sunlight of an area where the predetermined crop stored in the storage unit is planted. The crop discrimination device according to claim 1, wherein when the intensity of a component in a specific wavelength band matches, it is determined that the predetermined crop is planted in the certain area.
ことを特徴とする請求項1または請求項2に記載の作物判別装置。 The said storage part memorize | stores the intensity | strength of the component of three or more types of wavelength bands in each of the some pixel containing the at least one part area | region of the said certain area. The Claim 1 or Claim 2 characterized by the above-mentioned. Crop discrimination device.
ことを特徴とする請求項1〜請求項3のいずれか1つに記載の作物判別装置。 The storage unit stores the intensity of the green component, the intensity of the red component, the intensity of the near infrared component, and the intensity of the short wavelength infrared component in each of a plurality of pixels including at least a part of the certain area. Do
Crop discriminating device according to any one of claims 1 to 3, characterized in that.
ある地域をカメラによって撮像することにより得られた、該ある地域の少なくとも一部の領域を含む複数のピクセルのそれぞれにおける複数の特定の波長帯域の成分の強度を記憶する記憶部に記憶された記憶内容に基づいて、前記特定の波長帯域ごとに、前記特定の波長帯域の成分の強度と、該強度に対応するピクセルの個数との関係を示す情報を生成し、
生成された情報に基づいて、前記特定の波長帯域ごとに、ピクセルの個数が最大である強度を特定し、
生成された情報に基づいて、前記特定の波長帯域ごとに、特定された強度を含む、前記最大のピクセルの個数からピクセルの個数が減少する範囲内の前記特定の波長帯域の成分の強度に対応するピクセルにおける前記特定の波長帯域の成分の強度の平均値を算出し、
前記特定の波長帯域ごとに算出された平均値に基づいて、前記ある地域に所定の作物が植えられているか否かを判別する
処理を実行させることを特徴とする作物判別プログラム。 On the computer,
Obtained by imaging an area by a camera, a plurality of memory stored in the storage unit for storing the intensity of the components of a plurality of specific wavelength bands in each pixel including at least a partial area of the one region Based on the content, for each specific wavelength band, generating information indicating the relationship between the intensity of the component of the specific wavelength band and the number of pixels corresponding to the intensity,
Based on the generated information, for each of the specific wavelength bands, identify the intensity with the maximum number of pixels,
Corresponding to the intensity of the component of the specific wavelength band within the range in which the number of pixels decreases from the maximum number of pixels, including the specified intensity, for each specific wavelength band based on the generated information Calculating an average value of the components of the specific wavelength band in the pixel to be
The specific based on the average value calculated for each wavelength band, the certain crop judgment program for causing to execute a process for determining whether the region to a Jo Tokoro crop is planted.
ある地域をカメラによって撮像することにより得られた、該ある地域の少なくとも一部の領域を含む複数のピクセルのそれぞれにおける複数の特定の波長帯域の成分の強度を記憶する記憶部に記憶された記憶内容に基づいて、前記特定の波長帯域ごとに、前記特定の波長帯域の成分の強度と、該強度に対応するピクセルの個数との関係を示す情報を生成し、
生成された情報に基づいて、前記特定の波長帯域ごとに、ピクセルの個数が最大である強度を特定し、
生成された情報に基づいて、前記特定の波長帯域ごとに、特定された強度を含む、前記最大のピクセルの個数からピクセルの個数が減少する範囲内の前記特定の波長帯域の成分の強度に対応するピクセルにおける前記特定の波長帯域の成分の強度の平均値を算出し、
前記特定の波長帯域ごとに算出された平均値に基づいて、前記ある地域に所定の作物が植えられているか否かを判別する
処理を実行することを特徴とする作物判別方法。 Computer
Obtained by imaging an area by a camera, a plurality of memory stored in the storage unit for storing the intensity of the components of a plurality of specific wavelength bands in each pixel including at least a partial area of the one region Based on the content, for each specific wavelength band, generating information indicating the relationship between the intensity of the component of the specific wavelength band and the number of pixels corresponding to the intensity,
Based on the generated information, for each of the specific wavelength bands, identify the intensity with the maximum number of pixels,
Corresponding to the intensity of the component of the specific wavelength band within the range in which the number of pixels decreases from the maximum number of pixels, including the specified intensity, for each specific wavelength band based on the generated information Calculating an average value of the components of the specific wavelength band in the pixel to be
The specific based on the average value calculated for each wavelength band, the certain crop discriminating method characterized by performing a process for determining whether the region to a Jo Tokoro crop is planted.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013175908A JP6171732B2 (en) | 2013-08-27 | 2013-08-27 | Crop discrimination device, crop discrimination program and crop discrimination method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013175908A JP6171732B2 (en) | 2013-08-27 | 2013-08-27 | Crop discrimination device, crop discrimination program and crop discrimination method |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2015043707A JP2015043707A (en) | 2015-03-12 |
JP6171732B2 true JP6171732B2 (en) | 2017-08-02 |
Family
ID=52669832
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2013175908A Active JP6171732B2 (en) | 2013-08-27 | 2013-08-27 | Crop discrimination device, crop discrimination program and crop discrimination method |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP6171732B2 (en) |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003006612A (en) * | 2001-06-20 | 2003-01-10 | Ntt Data Corp | Device and method for predicting harvest |
JP2012196167A (en) * | 2011-03-20 | 2012-10-18 | Fujitsu Ltd | Plant species identification method |
JP5673823B2 (en) * | 2011-06-29 | 2015-02-18 | 富士通株式会社 | Plant species identification device, method and program |
JP5886054B2 (en) * | 2012-01-16 | 2016-03-16 | 株式会社日立製作所 | Image analysis system and method |
-
2013
- 2013-08-27 JP JP2013175908A patent/JP6171732B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
JP2015043707A (en) | 2015-03-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhao et al. | Immature green citrus detection based on colour feature and sum of absolute transformed difference (SATD) using colour images in the citrus grove | |
EP3159790B1 (en) | Program generation device, program generation method, and program | |
JP4869843B2 (en) | Cell image processing apparatus and cell image processing method | |
EP3664019A1 (en) | Information processing device, information processing program, and information processing method | |
US20170186201A1 (en) | Image processing apparatus, control method and record medium | |
US11645743B2 (en) | Method, medium, and system for detecting potato virus in a crop image | |
US9996554B2 (en) | Search apparatus, search method, and storage medium | |
Sabrol et al. | Intensity based feature extraction for tomato plant disease recognition by classification using decision tree | |
JP6299300B2 (en) | Image processing apparatus and image processing method | |
US8175335B2 (en) | Content adaptive detection of images with stand-out object | |
JP2009086901A (en) | Age estimation system and age estimation method | |
EP3156897B1 (en) | Program generation device, program generation method and program | |
Bankar et al. | Plant disease detection techniques using canny edge detection & color histogram in image processing | |
EP3316184B1 (en) | Program generating device, program generating method, and generating program | |
Müller et al. | Multi-disease detection in retinal imaging based on ensembling heterogeneous deep learning models | |
JP2013090752A (en) | Spot classification method, spot classification device and spot classification program | |
JP2015088040A (en) | Authentication device, authentication method, and program | |
JP6586782B2 (en) | Information processing apparatus, information processing program, and information processing system | |
US11036974B2 (en) | Image processing apparatus, image processing method, and storage medium | |
WO2016158768A1 (en) | Clustering device and machine learning device | |
US20220222820A1 (en) | Image processing apparatus, image processing method, and program | |
US20170236035A1 (en) | System, apparatus and method for configuration of industrial vision control modules | |
JP2017053699A (en) | Wavelength decision method of near-infrared image capturing used in substance determination, and substance determination method using near-infrared image | |
JP6171732B2 (en) | Crop discrimination device, crop discrimination program and crop discrimination method | |
Martínez et al. | Improving unsupervised saliency detection by migrating from RGB to multispectral images |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20160510 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20170313 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20170328 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20170523 |
|
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: 20170606 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20170619 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6171732 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |