JP2020041841A - Change rate estimating apparatus, change rate estimating method, state estimating device, and state estimating method - Google Patents

Change rate estimating apparatus, change rate estimating method, state estimating device, and state estimating method Download PDF

Info

Publication number
JP2020041841A
JP2020041841A JP2018167813A JP2018167813A JP2020041841A JP 2020041841 A JP2020041841 A JP 2020041841A JP 2018167813 A JP2018167813 A JP 2018167813A JP 2018167813 A JP2018167813 A JP 2018167813A JP 2020041841 A JP2020041841 A JP 2020041841A
Authority
JP
Japan
Prior art keywords
rice
change rate
rate
target rice
image data
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2018167813A
Other languages
Japanese (ja)
Other versions
JP7085952B2 (en
Inventor
晴雄 冨田
Haruo Tomita
晴雄 冨田
沙希 中山
Saki NAKAYAMA
沙希 中山
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Osaka Gas Co Ltd
Original Assignee
Osaka Gas Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Osaka Gas Co Ltd filed Critical Osaka Gas Co Ltd
Priority to JP2018167813A priority Critical patent/JP7085952B2/en
Publication of JP2020041841A publication Critical patent/JP2020041841A/en
Application granted granted Critical
Publication of JP7085952B2 publication Critical patent/JP7085952B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Length Measuring Devices By Optical Means (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

To provide a change rate estimating apparatus with which it is possible to easily estimate a change rate of a shape in the thickness direction of rice without requiring a device of complex configuration.SOLUTION: The change rate estimating apparatus comprises: an imaging device 2 for imaging the state of rice immersed in water; an image data acquisition unit 11 for acquiring the two-dimensional image data of the rice imaged by the imaging device 2 as seen from the thickness direction of the rice; a crack rate calculation unit 16 for referring to the two-dimensional image data acquired by the image data acquisition unit 11 and calculating a vertical crack rate; and an estimation unit 19 for acquiring the vertical crack rate from the crack rate calculation unit 16, and estimating the change rate of shape of the rice in the thickness direction on the basis of the correlation of the vertical crack rate with a predetermined vertical crack rate and change rate of shape in the thickness direction of rice.SELECTED DRAWING: Figure 1

Description

本発明は、対象米に関する情報を推定する変化率推定装置及び変化率推定方法に関し、また、推定した情報を基に、対象米の比重、密度又は吸水率を算出して対象米の状態を推定する状態推定装置及び状態推定方法に関する。   The present invention relates to a change rate estimating apparatus and a change rate estimating method for estimating information on a target rice, and further calculates a specific gravity, density or water absorption of the target rice to estimate a state of the target rice based on the estimated information. And a state estimating method.

例えば、酒造において、精米工程や吸水工程、蒸工程といった原料米に対する加工工程は、その後の醸造工程を左右する重要なプロセスである。大吟醸酒などの製造をする際の吸水工程では、原料米が吸水し過ぎないよう浸漬時間を短縮して原料米の吸水を制限する、いわゆる限定吸水という作業を行うことがあり、この限定吸水では1%単位での吸水率の調整が求められることから、製造する酒類の品質向上を図る上では、原料米の吸水状態をリアルタイムで詳細に把握することが極めて重要となる。   For example, in sake brewing, a processing step for raw rice, such as a rice polishing step, a water absorption step, and a steaming step, is an important process that determines the subsequent brewing step. In the water absorption process when manufacturing Daiginjo sake etc., so-called limited water absorption is sometimes performed, which is to limit the water absorption of the raw rice by shortening the immersion time so that the raw rice does not absorb too much water. In this case, it is necessary to adjust the water absorption in units of 1%. Therefore, in order to improve the quality of liquor to be produced, it is extremely important to grasp the state of water absorption of raw rice in detail in real time.

原料米の吸水変化は、原料米の状態(品種や精米歩合、含水率など)や浸漬条件(時間や温度、濃度など)により大きく変化することが知られており、従来、このような原料米の吸水変化を評価する手法としては、浸漬前後の原料米の重量を計量機で計測し、浸漬前後の重量差から吸水量を計算する手法が一般的であった。しかしながら、このような手法では、吸水量の経時変化や吸水分布などを正確に把握することができない。そのため、実際の酒造現場では、目視により原料米の色や形状の変化から吸水分布を推測し、上述した吸水量の計算結果と合わせて吸水変化を推測しており、原料米の吸水状態を客観的にリアルタイムで評価可能な手法の開発が求められていた。   It is known that the change in water absorption of raw rice varies greatly depending on the state of the raw rice (variety, polishing rate, moisture content, etc.) and immersion conditions (time, temperature, concentration, etc.). As a method of evaluating the change in water absorption, a method of measuring the weight of raw rice before and after immersion with a weighing machine and calculating the water absorption from the difference in weight before and after immersion was generally used. However, with such a method, it is not possible to accurately grasp the change over time in water absorption and the distribution of water absorption. Therefore, at an actual brewing site, the water absorption distribution is estimated visually from the change in the color and shape of the raw rice, and the change in water absorption is estimated in combination with the above-described calculation result of the water absorption. The development of a method that can be evaluated in real time was required.

ところで、本願出願人は、画像観察により炊飯対象米の吸水量や吸水分布などを判定する手法として、特許文献1記載の含水状態判定方法を提案している。この含水状態判定方法によれば、水に浸漬した状態に配置された炊飯対象米を光学的に監視して、炊飯対象米の含水率に関する含水率関連情報を経時的に測定することができるため、含水率関連情報を基に炊飯対象米の吸水量や吸水分布を評価することが可能となる。具体的には、炊飯対象米の吸水に伴う当該炊飯対象米の大きさの変化を経時的に測定して、吸水前の炊飯対象米の大きさに対して、吸水後の炊飯対象米の大きさが変化(膨張)した比率を膨張率とすることができるので、この膨張率を基にして炊飯対象米の吸水量や吸水分布を評価することができる。   Incidentally, the applicant of the present application has proposed a water-containing state determination method described in Patent Literature 1 as a method of determining the water absorption and the water absorption distribution of rice to be cooked by image observation. According to this water content determination method, since the rice to be cooked placed in a state of being immersed in water is optically monitored, the water content related information on the water content of the rice to be cooked can be measured over time. In addition, it is possible to evaluate the water absorption and the water absorption distribution of the rice to be cooked based on the water content-related information. Specifically, a change in the size of the rice to be cooked due to water absorption of the rice to be cooked is measured over time, and the size of the rice to be cooked after water absorption is compared with the size of the rice to be cooked before water absorption. Since the ratio of the change (expansion) can be used as the expansion rate, the water absorption and the water absorption distribution of the rice to be cooked can be evaluated based on the expansion rate.

特開2016−105074号公報JP-A-2006-105074

上記特許文献1記載の含水状態判定方法によれば、炊飯対象米の吸水量や吸水分布(吸水状態)を評価することができる。しかしながら、米の吸水状態は米の品種によって大きく変わるものであり、炊飯対象米は、吸水によって主に鉛直方向に直交する方向に膨張するのに対し、酒造対象米は、吸水によって鉛直方向に直交する方向だけでなく、鉛直方向にも膨張する。したがって、上記特許文献1記載の含水状態判定方法を用いた場合、例えば、水に浸漬した状態に配置された炊飯対象米を鉛直方向上方から光学的に監視するだけでは、鉛直方向に直交する方向への膨張の度合いを知ることはできても、鉛直方向への膨張の度合いを知ることができない。そのため、上記特許文献1記載の含水状態判定方法を使用して、鉛直方向へ膨張する米(例えば、酒造対象米)の含水状態を評価した場合、含水率関連情報は鉛直方向への膨張が反映されないため、測定した含水率関連情報と実際の酒造対象米の含水率との間にズレが生じて、結果的に、酒造対象米の吸水状態を正確に評価することができない。   According to the method for determining a water-containing state described in Patent Document 1, it is possible to evaluate the water absorption amount and water absorption distribution (water absorption state) of the rice to be cooked. However, the water absorption state of rice varies greatly with rice varieties, and rice to be cooked expands mainly in the direction perpendicular to the vertical direction due to water absorption, while rice to be brewed is orthogonal to the vertical direction due to water absorption. Expand not only in the vertical direction but also in the vertical direction. Therefore, when the method for determining a water-containing state described in Patent Document 1 described above is used, for example, simply monitoring the rice to be cooked placed in a state of being immersed in water optically from above in the vertical direction, the direction orthogonal to the vertical direction. It is possible to know the degree of expansion in the vertical direction, but not to know the degree of expansion in the vertical direction. Therefore, when the water-containing state of the rice expanding in the vertical direction (eg, rice for sake brewing) is evaluated using the water-containing state determination method described in Patent Document 1, the water content-related information reflects the vertical expansion. As a result, a difference occurs between the measured moisture content-related information and the actual moisture content of the rice to be brewed, and as a result, the water absorption state of the rice to be brewed cannot be accurately evaluated.

つまり、米の品種等を問わずに米の吸水状態を精度よく評価する上では、評価対象である米の2次元方向における形状の変化率(2次元形状の変化率)だけを基に考えたのでは不十分であり、2次元方向と直交する方向における形状の変化率や、3次元方向における形状の変化率(3次元形状の変化率)を考慮することが肝要である。   In other words, in order to accurately evaluate the water absorption state of rice irrespective of rice varieties, etc., it was considered based only on the rate of change in the two-dimensional shape of the rice to be evaluated (change rate of the two-dimensional shape). Is not sufficient, and it is important to consider the change rate of the shape in the direction orthogonal to the two-dimensional direction and the change rate of the shape in the three-dimensional direction (change rate of the three-dimensional shape).

上記特許文献1記載の含水状態判定方法を用いて、鉛直方向と直交する方向からも米を光学的に監視するようにすれば、鉛直方向における形状の変化率や、3次元形状の変化率を得ることができると考えられるが、この場合、米を光学的に監視するための装置(対象米に光を照射する照明装置や対象米からの透過光を撮像するカメラ装置など)を各方向ごとに設ける必要があるため、装置構成が複雑となり、装置コストも嵩むという懸念がある。   If the rice is optically monitored from a direction perpendicular to the vertical direction using the water content determination method described in Patent Document 1, the rate of change of the shape in the vertical direction and the rate of change of the three-dimensional shape can be reduced. In this case, a device for optically monitoring the rice (such as a lighting device for irradiating the target rice with light or a camera device for capturing the transmitted light from the target rice) may be provided in each direction. Therefore, there is a concern that the device configuration becomes complicated and the device cost increases.

また、上記2次元方向と直交する方向における形状の変化率や3次元形状の変化率は、米の吸水状態の評価に利用できるだけでなく、例えば、上記2次元方向と直交する方向(米の厚み方向)における形状の変化率や3次元形状の変化率を基に米の比重を算出すれば、心白の状態を評価することも可能になるなど、米の2次元方向と直交する方向における形状の変化率や米の3次元形状の変化率を求める利点は多く、これらを容易に求めることができる手法の開発が求められているという実情もある。   In addition, the rate of change of the shape and the rate of change of the three-dimensional shape in the direction orthogonal to the two-dimensional direction can be used not only for evaluating the water absorption state of rice, but also for example, in the direction orthogonal to the two-dimensional direction (rice thickness). Calculating the specific gravity of rice based on the rate of change of the shape in the direction) or the rate of change of the three-dimensional shape, it is possible to evaluate the whiteness of the rice. There are many advantages of obtaining the rate of change of the rice and the rate of change of the three-dimensional shape of rice, and there is a situation in which the development of a method that can easily obtain these rates is required.

本発明は以上の実情に鑑みなされたものであり、複雑な構成の装置を要することなく、米の厚み方向における形状の変化率を容易に求める、或いは、米の2次元形状の変化率を補正することによって容易に3次元形状の変化率を求めることができる変化率推定装置及び変化率推定方法の提供、並びに、上記2次元方向と直交する方向における形状の変化率や3次元形状の変化率を基に、米の比重、密度又は吸水率を容易に推定することができる状態推定装置及び状態推定方法の提供を、その目的とする。   SUMMARY OF THE INVENTION The present invention has been made in view of the above circumstances, and can easily determine the rate of change in the shape of rice in the thickness direction or correct the rate of change of the two-dimensional shape of rice without the need for a device having a complicated configuration. Of a change rate estimating apparatus and a change rate estimating method that can easily obtain a change rate of a three-dimensional shape by performing the method, and a change rate of a shape and a change rate of a three-dimensional shape in a direction orthogonal to the two-dimensional direction. It is an object of the present invention to provide a state estimating apparatus and a state estimating method which can easily estimate the specific gravity, density or water absorption of rice based on the above.

上記目的を達成するための本発明に係る変化率推定装置の特徴構成は、
水に浸漬した状態の対象米を撮影する撮影手段と、
前記撮影手段で撮影した前記対象米の、当該対象米を厚み方向から視た2次元画像データを取得する画像データ取得手段と、
前記画像データ取得手段で取得された2次元画像データを参照して、当該2次元画像データを解析し、複数の対象米のうち、前記対象米の長軸方向に沿う縦断面が開いた縦割れが発生している対象米の割合である縦割れ率を算出する割れ率算出手段と、
前記縦割れ率を前記割れ率算出手段から取得し、前記縦割れ率、及び予め定められた前記対象米の縦割れ率と前記厚み方向における形状の変化率との相関関係を基に、前記対象米の、前記厚み方向における形状の変化率を推定する推定手段とを備える点にある。
また、上記目的を達成するための本発明に係る変化率推定方法の特徴構成は、
水に浸漬した状態の対象米を撮影する撮影工程と、
前記撮影した前記対象米の、当該対象米を厚み方向から視た2次元画像を基に、複数の対象米のうち、前記対象米の長軸方向に沿う縦断面が開いた縦割れが発生している対象米の割合である縦割れ率を算出する割れ率算出工程と、
前記算出した縦割れ率、及び予め定められた前記対象米の縦割れ率と前記厚み方向における形状の変化率との相関関係を基に、前記対象米の、前記厚み方向における形状の変化率を推定する推定工程とを実行する点にある。
The characteristic configuration of the change rate estimation device according to the present invention for achieving the above object is as follows.
Photographing means for photographing the target rice immersed in water,
Image data acquisition means for acquiring two-dimensional image data of the target rice photographed by the photographing means, as viewed from the thickness direction of the target rice;
The two-dimensional image data is analyzed with reference to the two-dimensional image data acquired by the image data acquiring means, and a longitudinal crack having an open longitudinal section along a long axis direction of the target rice among a plurality of target rice is opened. Crack rate calculating means for calculating a vertical crack rate which is a ratio of target rice in which cracks have occurred,
The vertical cracking rate is obtained from the cracking rate calculating means, and the vertical cracking rate, and the predetermined correlation between the vertical cracking rate of the target rice and the change rate of the shape in the thickness direction are determined based on the correlation. Estimating means for estimating the rate of change in the shape of rice in the thickness direction.
Further, the characteristic configuration of the change rate estimation method according to the present invention for achieving the above object,
A photographing process of photographing the target rice immersed in water,
On the basis of a two-dimensional image of the target rice that has been photographed and the target rice viewed from the thickness direction, a vertical crack having an open longitudinal section along the long axis direction of the target rice among a plurality of target rice occurs. Cracking ratio calculating step of calculating a vertical cracking ratio which is a ratio of target rice having
Based on the calculated vertical cracking rate, and a predetermined correlation between the vertical cracking rate of the target rice and the change rate of the shape in the thickness direction, the target rice, the change rate of the shape in the thickness direction is calculated. And an estimating step of estimating.

ここで、本願発明者は、米の吸水率や吸水分布を画像観察により判定する手法について研究する過程で、水に浸漬した状態の米の割れ方が米によって異なり、米の割れ方によって吸水率や吸水状態、言い換えれば、米の3次元形状の変化率に違いがあることを発見した。割れ方の判断方法としては、目視で判断する方法もあるが、刻一刻と変化する米の割れ方を客観的に数値として記録することは困難である。そこで、本願発明者らは鋭意研究を重ね、その結果、2次元画像データについて人工知能を利用した自動解析を行って米の割れ方を数値化することに成功し、更に、対象米のうち縦割れが発生している対象米の割合が、対象米の厚み方向における形状の変化率と相関関係があるという新たな知見を得た。   Here, the inventor of the present application, in the course of studying a method of determining the water absorption rate and water absorption distribution of rice by image observation, the manner of cracking of the rice immersed in water differs depending on the rice, and the rate of water absorption depends on the way the rice cracks. And the state of water absorption, in other words, the rate of change in the three-dimensional shape of rice was different. As a method of determining how to break, there is a method of visually checking, but it is difficult to objectively record the changing manner of rice that changes every moment as a numerical value. Therefore, the inventors of the present application have conducted intensive studies, and as a result, succeeded in numerically analyzing the cracking of rice by performing automatic analysis using two-dimensional image data on the two-dimensional image data. New knowledge was obtained that the ratio of cracked rice was correlated with the rate of shape change in the thickness direction of the rice.

本発明に係る変化率推定装置及び変化率推定方法は、かかる知見に基づいてなされたものであり、上記特徴構成においては、まず、水に浸漬した状態の対象米の、厚み方向から視た2次元画像を撮影手段によって撮影し、ついで、撮影手段が撮影した2次元画像データを画像データ取得手段が取得する。次に、割れ率算出手段は、2次元画像データを参照して当該2次元画像データを解析し、複数の対象米のうち、縦割れが発生している対象米の割合である縦割れ率を算出する。   The change rate estimating apparatus and the change rate estimating method according to the present invention are based on such knowledge, and in the above-described feature configuration, first, the target rice in a state of being immersed in water is viewed from the thickness direction. A two-dimensional image is photographed by the photographing means, and two-dimensional image data photographed by the photographing means is acquired by the image data acquiring means. Next, the cracking ratio calculating means analyzes the two-dimensional image data with reference to the two-dimensional image data, and calculates a vertical cracking ratio, which is a ratio of target rice having a vertical crack, among a plurality of target rices. calculate.

しかる後、推定手段は、縦割れ率を取得し、縦割れ率、及び、予め定められた対象米の縦割れ率と厚み方向における形状の変化率との相関関係を基にして、対象米の厚み方向における形状の変化率を推定する。   Thereafter, the estimating means obtains the vertical cracking rate, and based on the vertical cracking rate, and a predetermined correlation between the vertical cracking rate of the target rice and the shape change rate in the thickness direction, the target rice The rate of change of the shape in the thickness direction is estimated.

このように、上記特徴構成によれば、対象米の厚み方向における形状の変化率と相関関係を有する縦割れ率を算出し、縦割れ率、及び上記相関関係を基に、対象米の厚み方向における形状の変化率を推定することができるため、従来のように厚み方向と直交する方向から対象米を光学的に監視することなく、当該対象米の厚み方向における形状の変化率を容易に求めることができる。よって、求めた厚み方向における形状の変化率を基にして、例えば、対象米の吸水状態や対象米の心白の状態などを評価することができるようになる。   As described above, according to the characteristic configuration, the vertical crack rate having a correlation with the rate of change in the shape of the target rice in the thickness direction is calculated, and based on the vertical crack rate and the correlation, the thickness direction of the target rice is calculated. Since the rate of change of the shape of the target rice can be estimated, the rate of change of the shape of the target rice in the thickness direction can be easily obtained without optically monitoring the target rice from a direction orthogonal to the thickness direction as in the related art. be able to. Therefore, for example, the water absorption state of the target rice, the whiteness of the target rice, and the like can be evaluated based on the obtained change rate of the shape in the thickness direction.

上記目的を達成するための本発明に係る状態推定装置の特徴構成は、
上記変化率推定装置と、
前記厚み方向における形状の変化率を前記推定手段から取得し、前記取得した厚み方向における形状の変化率を基に、前記対象米の比重、密度又は吸水率を算出する状態算出手段とを備える点にある。
また、上記目的を達成するための本発明に係る状態推定方法の特徴構成は、
上記変化率推定方法を行った後に、前記変化率推定方法で推定した前記厚み方向における形状の変化率を基に、前記対象米の比重、密度又は吸水率を算出する状態算出工程を実行する点にある。
The feature configuration of the state estimation device according to the present invention for achieving the above object is as follows.
The change rate estimation device,
A point calculating unit that obtains the shape change rate in the thickness direction from the estimating unit and calculates the specific gravity, density, or water absorption of the target rice based on the obtained shape change rate in the thickness direction. It is in.
The feature configuration of the state estimation method according to the present invention for achieving the above object is as follows.
After performing the change rate estimation method, based on the change rate of the shape in the thickness direction estimated by the change rate estimation method, performing a state calculation step of calculating the specific gravity, density, or water absorption of the target rice. It is in.

上記特徴構成によれば、光学的に監視する装置を要することなく容易に求められる厚み方向と直交する方向における形状の変化率を基に、対象米の比重、密度又は吸水率を算出することができるため、比重、密度、吸水率という対象米の状態を従来よりも容易に推定することができる。   According to the above feature configuration, it is possible to calculate the specific gravity, density, or water absorption of the target rice based on the shape change rate in a direction orthogonal to the thickness direction, which is easily obtained without the need for an optical monitoring device. Therefore, the state of the target rice such as specific gravity, density and water absorption can be more easily estimated than before.

更に、本発明に係る状態推定装置の更なる特徴構成は、
前記画像データ取得手段で取得された2次元画像データを参照して、当該2次元画像データを解析し、前記厚み方向と直交する平面に沿った方向における前記対象米の面積を取得する面積情報取得手段を備え、
前記状態算出手段において、前記面積を前記面積情報取得手段から取得し、前記取得した厚み方向における形状の変化率、前記取得した面積、及び予め記憶された水に浸漬する前の前記対象米の厚みを基に算出される前記対象米の体積を基に、前記対象米の比重又は密度を算出する点にある。
また、本発明に係る状態推定方法の更なる特徴構成は、
前記撮影した前記対象米の、当該対象米を厚み方向から視た2次元画像を基に、前記厚み方向と直交する平面に沿った方向における前記対象米の面積を取得する面積情報取得工程を備え、
前記状態算出工程において、前記厚み方向における形状の変化率、前記取得した面積、及び水に浸漬する前の前記対象米の厚みを基に算出される前記対象米の体積を基に、前記対象米の比重又は密度を算出する点にある。
Further, a further characteristic configuration of the state estimation device according to the present invention is as follows.
Area information acquisition for analyzing the two-dimensional image data with reference to the two-dimensional image data acquired by the image data acquisition means and acquiring the area of the target rice in a direction along a plane orthogonal to the thickness direction With means,
In the state calculating means, the area is obtained from the area information obtaining means, and the obtained shape change rate in the thickness direction, the obtained area, and the thickness of the target rice before immersion in pre-stored water. The specific gravity or density of the target rice is calculated based on the volume of the target rice calculated based on the following formula.
Further, a further characteristic configuration of the state estimation method according to the present invention includes:
An area information acquisition step of acquiring an area of the target rice in a direction along a plane orthogonal to the thickness direction based on a two-dimensional image of the target rice that has been photographed and viewing the target rice from the thickness direction. ,
In the state calculation step, based on the shape change rate in the thickness direction, the acquired area, and the volume of the target rice calculated based on the thickness of the target rice before immersion in water, the target rice The point is to calculate the specific gravity or density.

上記特徴構成によれば、2次元画像データを基に厚み方向と直交する平面に沿った方向における対象米の面積を取得し、この取得した面積と、従来よりも容易に求められる厚み方向における形状の変化率とに加えて、水に浸漬する前の対象米の厚みを予め求めておき、これらを基に体積を算出して、この体積を基にして、対象米の比重又は密度を算出する。したがって、対象米の比重又は密度という状態を従来よりも容易に推定することができる。   According to the above feature configuration, the area of the target rice in a direction along a plane orthogonal to the thickness direction is obtained based on the two-dimensional image data, and the obtained area and the shape in the thickness direction that are more easily obtained than before are obtained. In addition to the rate of change, the thickness of the target rice before immersion in water is determined in advance, the volume is calculated based on these, and the specific gravity or density of the target rice is calculated based on this volume. . Therefore, the state of specific gravity or density of the target rice can be estimated more easily than before.

上記目的を達成するための本発明に係る変化率推定装置の特徴構成は、
水に浸漬した状態の対象米を撮影する撮影手段と、
前記撮影手段で撮影した前記対象米の、当該対象米を厚み方向から視た2次元画像データを取得する画像データ取得手段と、
前記画像データ取得手段で取得された2次元画像データを参照して、当該2次元画像データを解析し、前記厚み方向と直交する平面に沿った方向における前記対象米の2次元形状の変化率を算出する2次元形状変化率算出手段と、
前記2次元形状変化率算出手段から前記2次元形状の変化率を取得し、複数の対象米のうち、前記対象米の長軸方向に沿う縦断面が開いた縦割れが発生している対象米の割合である縦割れ率と前記厚み方向における形状の変化率との相関関係、及び前記取得した2次元形状の変化率を基に、前記2次元形状の変化率を補正して3次元形状の変化率として推定する推定手段とを備える点にある。
また、上記目的を達成するための本発明に変化率推定方法の特徴構成は、
水に浸漬した状態の対象米を撮影する撮影工程と、
前記撮影した前記対象米の、当該対象米を厚み方向から視た2次元画像を基に、前記厚み方向と直交する平面に沿った方向における前記対象米の2次元形状の変化率を算出する2次元形状変化率算出工程と、
複数の対象米のうち、前記対象米の長軸方向に沿う縦断面が開いた縦割れが発生している対象米の割合である縦割れ率と前記厚み方向における形状の変化率との相関関係、及び前記算出した2次元形状の変化率を基に、前記2次元形状の変化率を補正して3次元形状の変化率として推定する推定工程とを実行する点にある。
The characteristic configuration of the change rate estimation device according to the present invention for achieving the above object is as follows.
Photographing means for photographing the target rice immersed in water,
Image data acquisition means for acquiring two-dimensional image data of the target rice photographed by the photographing means, as viewed from the thickness direction of the target rice;
The two-dimensional image data is analyzed with reference to the two-dimensional image data acquired by the image data acquiring means, and the rate of change of the two-dimensional shape of the target rice in a direction along a plane orthogonal to the thickness direction is calculated. A two-dimensional shape change rate calculating means for calculating;
The rate of change of the two-dimensional shape is obtained from the two-dimensional shape change rate calculating means, and a target rice having a vertical crack having an open longitudinal section along a long axis direction of the target rice among a plurality of target rices The rate of change of the two-dimensional shape is corrected based on the correlation between the vertical cracking rate, which is the ratio of the shape, and the rate of change of the shape in the thickness direction, and the obtained rate of change of the two-dimensional shape. Estimating means for estimating the rate of change.
Further, the characteristic configuration of the change rate estimation method according to the present invention for achieving the above object,
A photographing process of photographing the target rice immersed in water,
Calculating a change rate of a two-dimensional shape of the target rice in a direction along a plane orthogonal to the thickness direction based on a two-dimensional image of the target rice taken in the thickness direction of the target rice; Dimensional shape change rate calculating step,
Among a plurality of target rice, a correlation between a vertical cracking rate which is a ratio of target rice in which a vertical section along a longitudinal direction along the long axis direction of the target rice is open and a rate of shape change in the thickness direction. And an estimation step of correcting the change rate of the two-dimensional shape based on the calculated change rate of the two-dimensional shape and estimating the change rate of the two-dimensional shape as the change rate of the three-dimensional shape.

上述したように、本願発明者らは、米の割れ方によって米の3次元形状の変化率に違いがあり、また、米のうち縦割れが発生している米の割合が、米の厚み方向における形状の変化率と相関関係があるという新たな知見を得た。   As described above, the inventors of the present invention differ in the rate of change in the three-dimensional shape of rice depending on how the rice is cracked, and the ratio of rice that has vertical cracks in the rice varies in the thickness direction of the rice. New knowledge that there is a correlation with the shape change rate in.

本発明に係る変化率推定装置及び変化率推定方法は、かかる知見に基づきなされたものであり、上記構成においては、まず、水に浸漬した状態の対象米の、厚み方向から視た2次元画像を撮影手段によって撮影し、ついで、撮影手段が撮影した2次元画像データを画像データ取得手段が取得する。次に、2次元形状変化率算出手段は、2次元画像データを参照して当該2次元画像データを解析し、厚み方向と直交する平面に沿った方向における対象米の2次元形状の変化率を算出する。   The change rate estimating apparatus and the change rate estimating method according to the present invention are based on such knowledge, and in the above configuration, first, a two-dimensional image of the target rice in a state of being immersed in water, viewed from the thickness direction. Is photographed by the photographing means, and the image data acquiring means acquires two-dimensional image data photographed by the photographing means. Next, the two-dimensional shape change rate calculating means analyzes the two-dimensional image data with reference to the two-dimensional image data, and calculates a two-dimensional shape change rate of the target rice in a direction along a plane orthogonal to the thickness direction. calculate.

しかる後、推定手段は、2次元形状の変化率を取得し、縦割れ率と厚み方向における形状の変化率との相関関係、及び取得した2次元形状の変化率を基にして、2次元形状変化率算出手段で算出された2次元形状の変化率を3次元形状の変化率として推定する。   Thereafter, the estimating means obtains the change rate of the two-dimensional shape, and calculates the two-dimensional shape based on the correlation between the vertical crack rate and the change rate of the shape in the thickness direction and the obtained change rate of the two-dimensional shape. The change rate of the two-dimensional shape calculated by the change rate calculating means is estimated as the change rate of the three-dimensional shape.

このように、上記特徴構成によれば、上記相関関係及び取得した2次元形状の変化率を基に、2次元形状の変化率を補正して3次元形状の変化率とすることができるため、対象米の品種等を問わず、当該対象米の3次元形状の変化率を容易に推定することができる。よって、推定した3次元形状の変化率を基にして、例えば、対象米の吸水状態や対象米の心白の状態などを評価することができるようになる。   As described above, according to the characteristic configuration, the change rate of the two-dimensional shape can be corrected to the change rate of the three-dimensional shape based on the correlation and the obtained change rate of the two-dimensional shape. Regardless of the variety of the target rice, the change rate of the three-dimensional shape of the target rice can be easily estimated. Therefore, based on the estimated change rate of the three-dimensional shape, for example, the water absorption state of the target rice, the whiteness of the target rice, and the like can be evaluated.

本発明に係る変化率推定装置の更なる特徴構成は前記推定手段において、前記相関関係に基づいて定まる縦割れ影響度を、前記取得した2次元形状の変化率に乗じて、前記2次元形状の変化率を前記3次元形状の変化率に補正する点にある。
また、本発明に係る変化率推定方法の更なる特徴構成は、前記推定工程において、前記相関関係に基づいて定まる縦割れ影響度を、前記算出した2次元形状の変化率に乗じて、前記2次元形状の変化率を前記3次元形状の変化率に補正する点にある。
A further characteristic configuration of the change rate estimating apparatus according to the present invention is that the estimating means multiplies the degree of influence of the vertical crack determined based on the correlation by the obtained change rate of the two-dimensional shape to obtain the two-dimensional shape. The point is that the rate of change is corrected to the rate of change of the three-dimensional shape.
Further, a further characteristic configuration of the change rate estimating method according to the present invention is that, in the estimating step, the degree of vertical crack influence determined based on the correlation is multiplied by the calculated change rate of the two-dimensional shape to obtain the two-dimensional shape. The point is that the change rate of the three-dimensional shape is corrected to the change rate of the three-dimensional shape.

ここで、「縦割れ影響度」とは、縦割れの存在が3次元形状の変化率に与える影響を示す指標であって、具体的には、予め複数の品種に関する測定結果から得られたデータベースを基に、縦割れ率と厚み方向の変化率との相関関係を示すグラフを作成し、このグラフの傾きを縦割れ影響度として予め決定されるものである。   Here, the “vertical crack influence degree” is an index indicating the influence of the presence of a vertical crack on the rate of change of the three-dimensional shape, and specifically, a database obtained from measurement results for a plurality of varieties in advance. A graph showing the correlation between the vertical cracking rate and the rate of change in the thickness direction is created based on the above, and the inclination of this graph is determined in advance as the vertical cracking influence degree.

上記特徴構成によれば、対象米の2次元形状の変化率に縦割れ影響度を乗じることによって、当該2次元形状の変化率を補正して対象米の3次元形状の変化率を求めることができる。よって、上記と同様に、対象米の3次元形状の変化率を容易に求めることができ、米の吸水状態や米の心白の状態などを評価することができるようになる。   According to the above feature configuration, it is possible to correct the change rate of the two-dimensional shape by multiplying the change rate of the two-dimensional shape of the target rice by the degree of vertical cracking to obtain the change rate of the three-dimensional shape of the target rice. it can. Therefore, similarly to the above, the change rate of the three-dimensional shape of the target rice can be easily obtained, and the water absorption state of the rice and the whiteness of the rice can be evaluated.

また、上記目的を達成するための本発明に係る状態推定装置の特徴構成は、上記いずれかの変化率推定装置と、
前記3次元形状の変化率を前記推定手段から取得し、前記取得した3次元形状の変化率を基に、前記対象米の比重、密度又は吸水率を算出する状態算出手段とを備える点にある。
また、上記目的を達成するための本発明に係る状態推定方法の特徴構成は、上記いずれかの変化率推定方法を行った後に、前記変化率推定方法で推定した前記3次元形状の変化率を基に、前記対象米の比重、密度又は吸水率を算出する状態算出工程を実行する点にある。
Further, the characteristic configuration of the state estimation device according to the present invention for achieving the above object, any one of the above change rate estimation device,
The three-dimensional shape change rate is obtained from the estimating means, and a state calculating means for calculating the specific gravity, density or water absorption of the target rice based on the obtained three-dimensional shape change rate. .
In addition, a characteristic configuration of the state estimation method according to the present invention for achieving the above object includes, after performing any one of the change rate estimation methods described above, calculating the change rate of the three-dimensional shape estimated by the change rate estimation method. Based on this, a state calculation step of calculating the specific gravity, density or water absorption of the target rice is performed.

上記特徴構成によれば、2次元形状の変化率を補正して求められる3次元形状の変化率、即ち、従来よりも容易に求められる3次元形状の変化率を基に、対象米の比重、密度又は吸水率を算出することができるため、比重、密度、吸水率という対象米の状態を従来よりも容易に推定することが可能となる。   According to the above-mentioned feature configuration, the specific gravity of the target rice is determined based on the change rate of the three-dimensional shape obtained by correcting the change rate of the two-dimensional shape, that is, the change rate of the three-dimensional shape easily obtained than before. Since the density or the water absorption can be calculated, the state of the target rice such as specific gravity, density, and water absorption can be more easily estimated than before.

本発明に係る状態推定装置の更なる特徴構成は、前記状態算出手段において、前記取得した3次元形状の変化率、並びに予め記憶された水に浸漬する前の前記対象米の厚み及び厚み方向と直交する平面に沿った方向における面積を基に算出される前記対象米の体積を基に、前記対象米の比重又は密度を算出する点にある。
また、本発明に係る状態推定方法の更なる特徴構成は、前記状態算出工程において、前記3次元形状の変化率、並びに予め記憶された水に浸漬する前の前記対象米の厚み及び厚み方向と直交する平面に沿った方向における面積を基に算出される前記対象米の体積を基に、前記対象米の比重又は密度を算出する点にある。
A further characteristic configuration of the state estimation device according to the present invention is that, in the state calculation means, the acquired rate of change of the three-dimensional shape, and the thickness and thickness direction of the target rice before immersion in water stored in advance. The point is that the specific gravity or density of the target rice is calculated based on the volume of the target rice calculated based on the area in the direction along the orthogonal plane.
Further, a further characteristic configuration of the state estimation method according to the present invention includes, in the state calculation step, a change rate of the three-dimensional shape, and a thickness and a thickness direction of the target rice before immersion in pre-stored water. The point is that the specific gravity or density of the target rice is calculated based on the volume of the target rice calculated based on the area in the direction along the orthogonal plane.

上記特徴構成によれば、従来よりも容易に求められる3次元形状の変化率に加えて、水に浸漬する前の対象米の厚み及び厚み方向と直交する平面に沿った方向における面積を予め求めておき、これらを基に体積を算出して、この体積を基にして、対象米の比重又は密度を算出する。したがって、対象米の比重、密度という状態を従来よりも容易に推定することができる。   According to the above-mentioned characteristic configuration, in addition to the change rate of the three-dimensional shape, which is more easily obtained than before, the thickness of the target rice before immersion in water and the area in a direction along a plane orthogonal to the thickness direction are obtained in advance. In advance, the volume is calculated based on these, and the specific gravity or density of the target rice is calculated based on the volume. Therefore, the state of specific gravity and density of the target rice can be more easily estimated than before.

尚、本願において「比重」とは、密度(g/cm)の逆数としての比重(cm/g)を意味する。 In the present application, “specific gravity” means specific gravity (cm 3 / g) as a reciprocal of density (g / cm 3 ).

第1実施形態に係る状態推定装置及び変化率推定装置の構成を示す図である。It is a figure showing composition of a state estimating device and a change rate estimating device concerning a 1st embodiment. 状態推定方法及び変化率推定方法の過程を説明するフローチャートである。It is a flowchart explaining the process of a state estimation method and a change rate estimation method. 画像処理の例を説明する図である。FIG. 4 is a diagram illustrating an example of image processing. 酒米の形態の例を示す図である。It is a figure showing an example of a form of sake rice. 酒米の2次元形状の変化率及び3次元形状の変化率と、吸水率との関係を示すグラフである。It is a graph which shows the relationship between the rate of change of the two-dimensional shape and the rate of change of the three-dimensional shape of sake rice, and water absorption. 縦割れ率と厚み変化率との関係を示すグラフである。It is a graph which shows the relationship between a vertical crack rate and a thickness change rate. 機械学習処理を説明するフローチャートである。It is a flowchart explaining a machine learning process. 学習回数と正解率との関係を示すグラフである。It is a graph which shows the relationship between the number of times of learning and the correct answer rate. 学習回数と正解率との関係を示すグラフである。It is a graph which shows the relationship between the number of times of learning and the correct answer rate. 学習回数と正解率との関係を示すグラフである。It is a graph which shows the relationship between the number of times of learning and the correct answer rate. 第2実施形態に係る状態推定装置及び変化率推定装置の構成を示す図である。It is a figure showing composition of a state estimating device and a change rate estimating device concerning a 2nd embodiment. 機械学習処理を説明するフローチャートである。It is a flowchart explaining a machine learning process. 妥当性判定処理を説明するフローチャートである。It is a flowchart explaining a validity determination process. 第3実施形態に係る状態推定装置及び変化率推定装置の構成を示す図である。It is a figure showing composition of a state estimating device and a change rate estimating device concerning a 3rd embodiment. 状態推定方法及び変化率推定方法の過程を説明するフローチャートである。It is a flowchart explaining the process of a state estimation method and a change rate estimation method.

[第1実施形態]
以下、図面を参照して本発明の一実施形態に係る状態推定装置1及び変化率推定装置20Aについて説明する。尚、本実施形態では、酒米の厚み方向における形状の変化率を推定する場合を例にとって説明する。
[First Embodiment]
Hereinafter, a state estimation device 1 and a change rate estimation device 20A according to an embodiment of the present invention will be described with reference to the drawings. In this embodiment, a case where the rate of change in the shape of sake rice in the thickness direction is estimated will be described as an example.

図1に示すように、本実施形態に係る状態推定装置1A(1)は、撮影装置(撮影手段)2と、画像データ取得部(画像データ取得手段)11と、画像処理部12と、面積情報取得部(面積情報取得手段)13と、割れ率算出部(割れ率算出手段)16と、推定部(推定手段)19と、状態算出部(状態算出手段)30とを備えるとともに、取り扱う情報を記憶する記憶装置35を備える。図1では、状態推定装置1が、撮影装置2、画像データ取得部11、画像処理部12、割れ率算出部16、推定部19及び記憶装置35を含む変化率推定装置20A(20)と、面積情報取得部13と、状態算出部30とを備えて構成され、また、画像データ取得部11と、画像処理部12と、面積情報取得部13と、割れ率算出部16と、推定部19と、状態算出部30と、記憶装置35とが、米分析装置10A(10)を構成する形態で図示している。   As shown in FIG. 1, a state estimation device 1A (1) according to the present embodiment includes an imaging device (imaging unit) 2, an image data obtaining unit (image data obtaining unit) 11, an image processing unit 12, an area An information acquisition unit (area information acquisition unit) 13, a crack ratio calculation unit (crack ratio calculation unit) 16, an estimation unit (estimation unit) 19, and a state calculation unit (state calculation unit) 30. Is provided. In FIG. 1, the state estimating device 1 includes a change rate estimating device 20A (20) including a photographing device 2, an image data acquiring unit 11, an image processing unit 12, a crack rate calculating unit 16, an estimating unit 19, and a storage device 35; It is configured to include an area information acquisition unit 13 and a state calculation unit 30. The image data acquisition unit 11, the image processing unit 12, the area information acquisition unit 13, the crack rate calculation unit 16, and the estimation unit 19 , The state calculator 30 and the storage device 35 are shown in a form that constitutes the rice analyzer 10A (10).

状態推定装置1の一部を構成する米分析装置10Aは、情報の演算処理機能や情報の入出力機能、情報の記憶機能などを備える1台又は複数台のコンピュータ装置などを用いて実現される。その場合、画像データ取得部11の機能と、画像処理部12の機能と、面積情報取得部13の機能と、割れ率算出部16(形態判定部17の機能、形態データ解析部18の機能)の機能と、推定部19の機能と、状態算出部30の機能をコンピュータ装置に実現させるプログラムを、そのコンピュータ装置にインストールしておけばよい。   The rice analyzing apparatus 10A constituting a part of the state estimating apparatus 1 is realized using one or a plurality of computer devices having an information processing function, an information input / output function, an information storage function, and the like. . In that case, the function of the image data acquisition unit 11, the function of the image processing unit 12, the function of the area information acquisition unit 13, and the crack rate calculation unit 16 (the function of the form determination unit 17, the function of the form data analysis unit 18) A program that causes a computer device to implement the functions of the above, the function of the estimation unit 19, and the function of the state calculation unit 30 may be installed in the computer device.

本実施形態の記憶装置35は、撮影装置2が撮影した画像データを記憶する画像データ記憶部35aと、面積情報取得部13で取得された面積情報を記憶する面積記憶部35bと、機械学習で利用する訓練データを記憶する訓練データ記憶部35cと、形態判定部17の判定結果を記憶する形態記憶部35dとを備える。ここで、面積記憶部35b及び形態記憶部35dは、それぞれ面積情報取得部13で取得された面積情報及び形態判定部17の判定結果を、水への浸漬を開始した後での一つの対象米の面積及び形態の変化履歴が分かる状態で記憶しても良い。また、本実施形態の記憶装置35には、水に浸漬する前の酒米の厚みが予め記憶されている。   The storage device 35 of the present embodiment includes an image data storage unit 35a that stores image data captured by the imaging device 2, an area storage unit 35b that stores the area information acquired by the area information acquisition unit 13, and machine learning. It includes a training data storage unit 35c that stores training data to be used, and a form storage unit 35d that stores the determination result of the form determination unit 17. Here, the area storage unit 35b and the form storage unit 35d store the area information acquired by the area information acquisition unit 13 and the determination result of the form determination unit 17 in one target rice after the start of immersion in water. May be stored in a state in which the change history of the area and the form can be understood. The thickness of the sake rice before immersion in water is stored in the storage device 35 of the present embodiment in advance.

撮影装置2は、酒米を水に浸漬した状態に配置する米保持部の鉛直方向上方に配置されており、酒米からの透過光を撮影して当該酒米の2次元画像を撮影することができる装置である。   The photographing device 2 is arranged vertically above a rice holding unit that arranges sake rice in a state of being immersed in water, and shoots a two-dimensional image of the sake rice by photographing transmitted light from the sake rice. It is a device that can do.

出力装置3は、画像情報や文字情報などを表示できる表示装置や、それら画像情報や文字情報などを紙などに印刷できる印刷装置などである。   The output device 3 is a display device that can display image information, character information, and the like, and a printing device that can print the image information, character information, and the like on paper.

図2は、本実施形態に係る変化率推定方法及び状態推定方法に含まれる、酒米の面積情報を取得する工程(面積情報取得工程)、酒米の割れ率を算出する工程(割れ率算出工程)、酒米の厚み方向における形状の変化を推定する工程(推定工程)、及び推定した厚み方向における形状の変化率を基に、酒米の比重を算出する工程(状態算出工程)を説明するフローチャートである。図2の工程#10が行われた後、即ち、水に浸漬した状態に配置された酒米を撮影装置2によって撮影する撮影工程が実行された後、工程#11において、画像データ取得部11は、水に浸漬した状態に配置された酒米を当該酒米の厚み方向から撮影装置2で撮影した2次元画像データを取得する。つまり、水に浸漬した状態に配置された酒米を、当該酒米の厚み方向から視た2次元画像データを取得する画像データ取得工程が実行される。例えば、撮影装置2は、水に浸漬した状態に配置された複数の酒米を撮影して1つの2次元画像データを得る。そして、その2次元画像データは米分析装置10Aに伝達され、画像データ取得部11が取得する。画像データ取得部11が取得した2次元画像データを記憶装置35の画像データ記憶部35aに記憶しておくこともできる。例えば、図3に示す画像データAは、画像データ取得部11が取得する2次元画像データである。   FIG. 2 shows a step of acquiring the area information of the sake rice (area information acquisition step) and a step of calculating the crack rate of the sake rice (crack rate calculation) included in the change rate estimation method and the state estimation method according to the present embodiment. Step), a step of estimating a change in the shape of the sake rice in the thickness direction (estimation step), and a step of calculating the specific gravity of the sake rice based on the estimated rate of change of the shape in the thickness direction (state calculation step). FIG. After the process # 10 of FIG. 2 is performed, that is, after the photographing process of photographing the sake rice arranged in the state of being immersed in water by the photographing device 2 is performed, in the process # 11, the image data acquiring unit 11 is performed. Acquires two-dimensional image data obtained by photographing the sake rice arranged in a state of being immersed in water with the photographing device 2 from the thickness direction of the sake rice. That is, an image data acquisition step of acquiring two-dimensional image data of the sake rice placed in a state of being immersed in water as viewed from the thickness direction of the sake rice is executed. For example, the imaging device 2 obtains one two-dimensional image data by imaging a plurality of sake rice arranged in a state of being immersed in water. Then, the two-dimensional image data is transmitted to the rice analyzing apparatus 10A, and the image data acquiring unit 11 acquires the two-dimensional image data. The two-dimensional image data acquired by the image data acquisition unit 11 can be stored in the image data storage unit 35a of the storage device 35. For example, image data A shown in FIG. 3 is two-dimensional image data acquired by the image data acquisition unit 11.

次に、図2の工程#12において、画像処理部12は、面積情報取得部13によって面積情報を取得する処理及び割れ率算出部16によって縦割れ率を算出する処理に先立って、画像データ取得部11が取得した2次元画像データにおいて、1つの酒米のみを含む領域の2次元酒米画像データを抽出する処理を含む画像処理を行って、その1つの酒米の処理済みの2次元画像データを生成する。つまり、画像データ取得工程で取得した2次元画像データにおいて、1つの酒米のみを含む領域の2次元酒米画像データを抽出する処理を含む画像処理を行って、その1つの酒米の処理済みの2次元画像データを生成する画像処理工程が実行される。図3は、画像処理の例を説明する図である。本実施形態では、画像処理部12が行う画像処理として、画像データ取得部11が取得した画像データにおいて、1つの酒米のみを含む領域の2次元酒米画像データを抽出する処理、その1つの酒米が所定の姿勢になるように酒米の2次元画像を回転させる処理、及び画像のコントラストを調整する処理を含む。   Next, in step # 12 of FIG. 2, the image processing unit 12 obtains image data prior to the process of obtaining area information by the area information obtaining unit 13 and the process of calculating the vertical cracking ratio by the cracking ratio calculation unit 16. In the two-dimensional image data acquired by the unit 11, image processing including a process of extracting two-dimensional sake rice image data of an area including only one sake rice is performed, and a processed two-dimensional image of the one sake rice is performed. Generate data. That is, in the two-dimensional image data acquired in the image data acquiring step, image processing including a process of extracting two-dimensional sake rice image data of a region including only one sake rice is performed, and the processing of the one sake rice is completed. An image processing step of generating the two-dimensional image data is performed. FIG. 3 is a diagram illustrating an example of image processing. In the present embodiment, as the image processing performed by the image processing unit 12, in the image data acquired by the image data acquiring unit 11, a process of extracting two-dimensional sake rice image data of an area including only one sake rice, The process includes a process of rotating the two-dimensional image of the sake rice so that the sake rice has a predetermined posture, and a process of adjusting the contrast of the image.

例えば、図3に示す2次元画像データBは、画像データ取得部11が取得した画像データAにおいて、1つの酒米の2次元全体画像を含むデータを抽出した画像である。図3に示す2次元画像データCは、2次元画像データBにおいて、1つの酒米以外の他の酒米の2次元画像を除去することで、1つの酒米のみを含む領域の2次元酒米画像データを抽出する処理を行った後の画像である。図3に示す2次元画像データDは、2次元画像データCに対して、1つの酒米が所定の姿勢になるように酒米の2次元画像を回転させる処理を行った後の画像である。この場合、酒米の長軸が鉛直方向に沿うような姿勢に酒米画像を回転させている。図3に示す2次元画像データEは、2次元画像データに対して、酒米の向きを調整する処理(例えば、上下反転処理や左右反転処理など)を行った後の画像である。図3に示す2次元画像データFは、2次元画像データEに対して、画像のコントラストを調整する処理を行った後の画像である。   For example, the two-dimensional image data B illustrated in FIG. 3 is an image obtained by extracting data including a two-dimensional whole image of one sake rice from the image data A acquired by the image data acquiring unit 11. The two-dimensional image data C shown in FIG. 3 is obtained by removing a two-dimensional image of sake rice other than one sake rice from the two-dimensional image data B to obtain a two-dimensional sake of an area including only one sake rice. It is an image after performing the process of extracting rice image data. The two-dimensional image data D illustrated in FIG. 3 is an image obtained by performing a process on the two-dimensional image data C to rotate the two-dimensional image of the sake rice so that one sake rice has a predetermined posture. . In this case, the sake rice image is rotated so that the long axis of the sake rice is along the vertical direction. The two-dimensional image data E illustrated in FIG. 3 is an image obtained by performing a process of adjusting the direction of sake rice (for example, a vertical inversion process or a horizontal inversion process) on the two-dimensional image data. The two-dimensional image data F shown in FIG. 3 is an image obtained by performing a process of adjusting the contrast of the image on the two-dimensional image data E.

ついで、面積情報取得部13は、画像データ取得部11が取得した2次元画像データを参照して、当該2次元画像データを解析し、酒米の厚み方向と直交する平面に沿った方向における酒米の面積を算出する。   Next, the area information acquisition unit 13 analyzes the two-dimensional image data with reference to the two-dimensional image data acquired by the image data acquisition unit 11, and outputs the sake in a direction along a plane orthogonal to the thickness direction of the sake rice. Calculate the area of rice.

本実施形態では、面積情報取得部13は、画像データ取得部11が取得し、画像処理部12が処理した2次元画像データを参照して、酒米の厚み方向と直交する平面に沿った方向における酒米の面積を、様々な画像認識技術を用いて自動で判定することができる。   In the present embodiment, the area information acquisition unit 13 refers to the two-dimensional image data acquired by the image data acquisition unit 11 and processed by the image processing unit 12, and refers to a direction along a plane orthogonal to the thickness direction of the sake rice. Can be automatically determined using various image recognition techniques.

上記のように、面積情報取得部13が酒米の面積を判定するために参照する入力データとしての2次元画像データは、画像処理部12が、画像データ取得部11が取得した画像データにおいて、1つの酒米のみを含む領域の2次元酒米画像データを抽出する処理、1つの酒米が所定の姿勢になるように酒米の2次元画像を回転させる処理、及び画像のコントラストを調整する処理を含む画像処理を行った後の処理済みの2次元画像データである。つまり、面積情報取得部13は、酒米の面積を判定するのに適した状態にした2次元画像データ(処理済みの2次元画像データ)から酒米の面積を判定できるので、その判定結果が適切である可能性が高まる。   As described above, the two-dimensional image data as the input data referred to by the area information acquisition unit 13 to determine the area of the sake rice is the image data acquired by the image data acquisition unit 11 by the image processing unit 12. A process of extracting two-dimensional sake rice image data of an area including only one sake rice, a process of rotating a two-dimensional image of sake rice so that one sake rice has a predetermined posture, and adjusting the image contrast. This is processed two-dimensional image data after performing image processing including processing. That is, the area information acquisition unit 13 can determine the area of the sake rice from the two-dimensional image data (processed two-dimensional image data) that is in a state suitable for determining the area of the sake rice, and the determination result is More likely to be appropriate.

図2の工程#13において、面積情報取得部13は、酒米の面積を判定する。つまり、画像処理工程で生成した処理済みの2次元画像データにおける1つの酒米の面積を判定する面積判定工程(面積情報取得工程)が実行される。本実施形態では、面積情報取得部13が生成した処理済みの2次元画像データ(図3の2次元画像データF)における1つの酒米の面積を判定する。尚、水に浸漬された状態にある酒米が複数ある場合、面積情報取得部13がそのうちのいくつの酒米の面積を判定するのかは適宜設定可能である。そして、面積情報取得部13による判定結果は、複数の酒米のそれぞれに付与された識別子と関連付けられて酒米ごとに面積記憶部35bに記憶される。面積記憶部35bには、面積情報取得部13で取得された面積情報が、水への浸漬を開始した後での1つの酒米の面積の変化履歴が分かる状態で記憶される。   In step # 13 of FIG. 2, the area information acquisition unit 13 determines the area of the sake rice. That is, an area determining step (area information obtaining step) for determining the area of one sake rice in the processed two-dimensional image data generated in the image processing step is executed. In the present embodiment, the area of one sake rice in the processed two-dimensional image data (two-dimensional image data F in FIG. 3) generated by the area information acquisition unit 13 is determined. When there are a plurality of sake rice immersed in water, how many of the sake rice the area information acquisition unit 13 determines the area of the sake rice can be appropriately set. Then, the determination result by the area information acquisition unit 13 is stored in the area storage unit 35b for each sake rice in association with the identifier assigned to each of the plurality of sake rice. The area storage unit 35b stores the area information acquired by the area information acquisition unit 13 in a state where the change history of the area of one sake rice after the start of immersion in water is known.

このように、本実施形態では、面積情報取得部13が、画像処理部12が生成した処理済みの2次元画像データにおける1つの酒米の面積を判定する。その結果、酒米の面積に関して、短時間で判定結果を得ることができる。   As described above, in the present embodiment, the area information obtaining unit 13 determines the area of one sake rice in the processed two-dimensional image data generated by the image processing unit 12. As a result, a determination result can be obtained in a short time for the area of sake rice.

尚、面積情報取得部13で取得された面積情報は、出力装置3から出力することもできる。   The area information acquired by the area information acquiring unit 13 can be output from the output device 3.

次に、割れ率算出部16は、画像データ取得部11が取得した2次元画像データを参照して、当該2次元画像データを解析し、複数の酒米のうち、縦割れが発生している酒米の割合である縦割れ率を算出する。尚、「縦割れ」は、酒米に生じた縦方向の亀裂の隙間が大きくなった状態である。   Next, the cracking rate calculating unit 16 analyzes the two-dimensional image data with reference to the two-dimensional image data acquired by the image data acquiring unit 11, and a vertical crack has occurred in a plurality of sake rice. Calculate the vertical cracking rate, which is the ratio of sake rice. The “vertical crack” is a state in which the gap of the vertical crack generated in the sake rice is large.

本実施形態では、割れ率算出部16は、形態判定部17と形態データ解析部18とを備える。例えば、形態判定部17は、酒米に割れが発生しているか否かを、様々な画像認識技術を用いて自動で判定することができる。後述する例では、形態判定部17は、CNN(Convolutional Neural Network)などのアルゴリズムを用いて酒米の形態を予測判定する場合について説明する。また、形態判定部17は、画像処理部12が生成した処理済みの2次元画像データが入力された場合に、酒米の形態に関する情報が出力されるように機械学習を行う場合を説明する。例えば、機械学習では、形態判定部17のアルゴリズムを構成するCNNの複数の重みパラメータなどが調整される。   In the present embodiment, the crack rate calculation unit 16 includes a shape determination unit 17 and a shape data analysis unit 18. For example, the morphological determination unit 17 can automatically determine whether cracked rice wine has been generated using various image recognition techniques. In the example described later, a case will be described where the morphological determination unit 17 predicts and determines the morphology of sake rice using an algorithm such as CNN (Convolutional Neural Network). Also, a case will be described where the morphological determination unit 17 performs machine learning such that when processed two-dimensional image data generated by the image processing unit 12 is input, information on the morphology of sake rice is output. For example, in machine learning, a plurality of weight parameters of the CNN constituting the algorithm of the form determination unit 17 are adjusted.

上記のように、形態判定部17が酒米の形態を判定するために参照する入力データとしての2次元画像データは、画像処理部12が、画像データ取得部11が取得した画像データにおいて、1つの酒米のみを含む領域の2次元酒米画像データを抽出する処理、1つの酒米が所定の姿勢になるように酒米の2次元画像を回転させる処理、及び画像のコントラストを調整する処理を含む画像処理を行った後の処理済みの2次元画像データである。つまり、形態判定部17は、酒米の形態を判定するのに適した状態にした2次元画像データ(処理済みの2次元画像データ)から酒米の形態を判定できるので、その判定結果が適切である可能性が高まる。   As described above, the two-dimensional image data as the input data that is referred to by the morphological determination unit 17 to determine the morphology of the sake rice is such that the image processing unit 12 uses one of the image data acquired by the image data acquisition unit 11. A process of extracting two-dimensional sake rice image data of an area including only one sake rice, a process of rotating a two-dimensional image of sake rice so that one sake rice has a predetermined posture, and a process of adjusting image contrast Is processed two-dimensional image data after image processing including. That is, the morphological determination unit 17 can determine the morphology of the sake rice from the two-dimensional image data (processed two-dimensional image data) in a state suitable for determining the morphology of the sake rice. Is more likely to be.

図2の工程#14において、形態判定部17は、割れがあるか否かを判定する。つまり、画像処理工程で生成した処理済みの2次元画像データにおける1つの酒米に割れがあるか否かを判定する形態判定工程(割れ率算出工程)が実行される。本実施形態では、形態判定部17は、蓄積された訓練データの機械学習の実行結果を用いて、画像処理部12が生成した処理済みの2次元画像データ(図3の2次元画像データF)における1つの酒米に割れがあるか否かを判定する。   In step # 14 of FIG. 2, the form determining unit 17 determines whether there is a crack. That is, a form determination step (crack rate calculation step) for determining whether or not one sake rice in the processed two-dimensional image data generated in the image processing step has a crack is executed. In the present embodiment, the form determination unit 17 uses the execution result of the machine learning of the accumulated training data to process the two-dimensional image data generated by the image processing unit 12 (two-dimensional image data F in FIG. 3). It is determined whether or not one sake rice has cracks.

本実施形態では、形態判定部17は、酒米が損傷していない形態か、酒米に亀裂が生じた形態であるか、酒米の長軸方向に沿う縦断面が開いた縦割れが生じた形態であるか、又は、酒米の短軸方向に沿う横断面が開いた横割れが生じた形態であるかの4種類の何れかを判定する。このうち、下記表1に示すように、酒米が損傷していない形態、及び、酒米に亀裂が生じた形態は、「割れが無い」という判定になり、酒米に縦割れがある形態、及び、酒米に横割れがある形態は、「割れが有る」という判定になる。尚、水に浸漬された状態にある酒米が複数ある場合、形態判定部17がそのうちのいくつの酒米の形態を判定するのかは適宜設定可能である。   In the present embodiment, the morphological determination unit 17 determines whether the sake rice is in an intact form, the sake rice has a crack, or a vertical crack having an open longitudinal section along the long axis direction of the sake rice occurs. One of the four types is determined, that is, the type of the brewed rice, or the type in which the horizontal cross section along the short axis direction of the sake rice has an open lateral crack. Among them, as shown in Table 1 below, the form in which the sake rice is not damaged and the form in which the sake rice has cracks are determined to be "no crack", and the form in which the sake rice has vertical cracks. , And the mode in which the sake rice has horizontal cracks is determined to be "cracked". When there are a plurality of sake rice immersed in water, it is possible to appropriately set how many of the sake rice the form determination unit 17 determines the form of the sake rice.

形態判定部17は、画像処理部12が生成した処理済みの2次元画像データが入力されると、酒米の形態に関する情報が出力されるようにアルゴリズムが構築されている。例えば、形態判定部17で実行されるCNNは、画像処理部12が生成した処理済みの2次元画像データが入力される入力層と、入力層に入力された処理済みの2次元画像データに対して畳み込み処理を行って特徴マップを得る畳み込み層と、畳み込み層から出力された特徴マップを縮小するプール層と、全ユニットを結合する全結合層と、出力層とを備えて階層型のネットワークを構成している。畳み込み層及びプール層の組み合わせは複数回繰り返し設けられ、全結合層も複数層設けられる。   The algorithm of the morphological determination unit 17 is constructed such that when the processed two-dimensional image data generated by the image processing unit 12 is input, information on the morphology of sake rice is output. For example, the CNN executed by the form determination unit 17 performs processing on the input layer to which the processed two-dimensional image data generated by the image processing unit 12 is input and the processed two-dimensional image data input to the input layer. A hierarchical network including a convolution layer that performs a convolution process to obtain a feature map, a pool layer that reduces a feature map output from the convolution layer, a fully connected layer that connects all units, and an output layer. Make up. The combination of the convolutional layer and the pool layer is repeatedly provided a plurality of times, and the total connection layer is also provided in a plurality of layers.

本実施形態では、形態判定部17の判定結果となる酒米の形態に関する情報として4種類の情報(損傷無し、亀裂、縦割れ、横割れ)を想定しているため、CNNの出力層のユニット数は4になる。出力層ではソフトマックス関数を尤度関数として用いることができる。ソフトマックス関数は、出力層の全ユニットの値を累積した値で各ユニットの値を除算して正規化するので、出力層の各ユニットの尤度は0から1の間の値をとる。そして、形態判定部17は、尤度が最大となるクラスを分類クラス、即ち、酒米の形態(判定結果)として決定する。そして、形態判定部17による判定結果は、複数の酒米のそれぞれに付与された識別子と関連付けられて酒米ごとに形態記憶部35dに記憶される。形態記憶部35dには、形態判定部17の判定結果が、水への浸漬を開始した後での1つの酒米の形態の変化履歴が分かる状態で記憶される。   In the present embodiment, four types of information (no damage, cracks, vertical cracks, and horizontal cracks) are assumed as the information on the morphology of the sake rice, which is the determination result of the morphological determination unit 17, and therefore, the unit of the output layer of the CNN. The number becomes four. In the output layer, a softmax function can be used as a likelihood function. The softmax function divides the value of each unit by the value obtained by accumulating the values of all units in the output layer and normalizes the value, so that the likelihood of each unit in the output layer takes a value between 0 and 1. Then, the form determining unit 17 determines the class having the maximum likelihood as the classification class, that is, the form (determination result) of the sake rice. The determination result by the form determining unit 17 is stored in the form storage unit 35d for each sake rice in association with the identifier given to each of the plurality of sake rice. The determination result of the configuration determination unit 17 is stored in the configuration storage unit 35d in a state where the change history of the configuration of one sake rice after the start of immersion in water is known.

図4は、酒米の状態の例を示す図である。具体的には、図4(a)は酒米が損傷していない形態を示す図であり、図4(b)は酒米において短軸方向に沿う亀裂が生じた形態である。また、図4(c)及び図4(d)は酒米に生じた亀裂の隙間が大きくなった状態である割れが生じた形態であり、割れの形態として、酒米の短軸方向に沿う横断面が開いた横割れが生じた形態(図4(c))と、酒米の長軸方向に沿う縦断面が開いた縦割れが生じた形態(図4(d))とがある。このように、形態判定部17は、割れの形態として、酒米の短軸方向に沿う横断面が開いた横割れが生じた形態であるか、又は、酒米の長軸方向に沿う縦断面が開いた縦割れが生じた形態であるかも判定する。   FIG. 4 is a diagram illustrating an example of a state of sake rice. Specifically, FIG. 4A is a view showing a form in which the sake rice is not damaged, and FIG. 4B is a form in which cracks are generated in the sake rice along the short axis direction. 4 (c) and 4 (d) show a form in which cracks are generated in a state in which the gaps of cracks generated in sake rice are large, and the form of the cracks is along the short axis direction of the sake rice. There is a form in which a horizontal crack with an open cross section is generated (FIG. 4C), and a form in which a vertical crack with an open vertical cross section along the long axis direction of sake rice is generated (FIG. 4D). As described above, the shape determining unit 17 may determine whether the crack is in a form in which a horizontal cross section along the short axis direction of the sake rice is open or a vertical cross section along the long axis direction of the sake rice. It is also determined whether or not is a form in which an open vertical crack has occurred.

以上のように、本実施形態では、人間が酒米の形態を判定するのではなく、形態判定部17が、蓄積された訓練データの機械学習の実行結果を用いて、画像処理部12が生成した処理済みの2次元画像データにおける1つの酒米の形態を判定する。その結果、酒米の形態に関して、短時間で、且つ、一貫性のある判定結果を得ることができる。   As described above, in the present embodiment, instead of the human judging the form of the sake rice, the form judging unit 17 generates the image processing unit 12 using the machine learning execution result of the accumulated training data. The form of one sake rice in the processed two-dimensional image data is determined. As a result, a consistent determination result can be obtained in a short time with respect to the form of sake rice.

図2の工程#15において、形態データ解析部18は、画像データ取得部11が取得した2次元画像データに含まれる複数の酒米のそれぞれについての形態判定部17の判定結果に基づいて縦割れ率を算出する。つまり、形態判定工程の判定結果に基づいて縦割れ率を算出する形態データ解析工程(割れ率算出工程)が実行される。形態記憶部35dには、複数の酒米のそれぞれについての形態判定部17の判定結果が、水への浸漬を開始した後での1つの酒米の形態の変化履歴が分かる状態で記憶されている。また、形態記憶部35dには、各酒米について、水への浸漬を開始した時点から現在(最新)までの少なくとも2以上の形態が変化履歴として記憶されている。そして、形態データ解析部18は、形態記憶部35dに記憶されている複数の酒米の現在(最新)の形態を読み出して、複数の酒米のうち、縦割れが発生している酒米の割合である縦割れ率を算出する。   In step # 15 of FIG. 2, the morphological data analysis unit 18 performs vertical cracking based on the determination result of the morphological determination unit 17 for each of the plurality of sake rice included in the two-dimensional image data acquired by the image data acquisition unit 11. Calculate the rate. That is, a form data analysis step (a crack rate calculation step) of calculating a vertical crack rate based on the determination result of the form determination step is executed. In the form storage unit 35d, the determination result of the form determination unit 17 for each of the plurality of sake rice is stored in a state where the change history of the form of one sake rice after the start of immersion in water is known. I have. Further, in the form storage unit 35d, at least two or more forms from the start of immersion in water to the present (latest) of each sake rice are stored as a change history. Then, the morphological data analysis unit 18 reads out the current (latest) morphology of the plurality of sake rice stored in the morphology storage unit 35d, and among the plurality of sake rice, the one of the sake rice in which the vertical crack occurs. Calculate the vertical crack ratio which is a ratio.

尚、形態判定部17の判定結果を、出力装置3から出力することもできる。例えば、複数の酒米のうち、割れが発生している酒米の割合(割れ率)などを随時算出して出力装置3からリアルタイムで出力することもできる。また、形態記憶部35dに、形態判定部17の判定結果が、水への浸漬を開始した後での1つの酒米の形態の変化履歴が分かる状態で記憶されている場合には、割れ率の時間的な変化などを算出して出力装置3から出力することもできる。   Note that the determination result of the form determination unit 17 can be output from the output device 3. For example, the ratio (crack rate) of cracked rice among a plurality of sake rice may be calculated as needed and output from the output device 3 in real time. If the determination result of the shape determination unit 17 is stored in the form storage unit 35d in a state in which the change history of the shape of one sake rice after the start of immersion in water is known, the cracking rate Can be calculated and output from the output device 3.

推定部19は、割れ率算出部16が算出した酒米の縦割れ率を取得し、この取得した縦割れ率、及び予め定められた酒米の縦割れ率と厚み方向における形状の変化率との相関関係を基に、酒米の厚み方向における形状の変化率を推定する。   The estimating unit 19 obtains the vertical cracking rate of the sake rice calculated by the cracking rate calculation unit 16 and obtains the obtained vertical cracking rate, and the predetermined vertical cracking rate of the sake rice and the change rate of the shape in the thickness direction. The change rate of the shape of the sake rice in the thickness direction is estimated based on the correlation.

本願発明者は、酒米の吸水率や吸水分布といった吸水状態を画像観察により判定する手法について研究する過程で、酒米のうち縦割れが発生している酒米の割合である縦割れ率が、酒米の厚み方向における形状の変化率と相関関係があるという新たな知見を得た。具体的に、本願発明者は、水に浸漬した状態に配置された複数の酒米について、2次元画像データを基に測定した酒米の厚み方向と直交する平面における2次元形状の変化率の平均値を算出するとともに、2次元形状の変化率を測定した時点で水分率計を用いて測定した複数の酒米の吸水率の平均値を算出することで、酒米の2次元形状の変化率と吸水率との関係をプロットした。また、水に浸漬した状態に配置された複数の酒米について、既知の手法を用いて測定した3次元形状の変化率の平均値を算出するとともに、3次元形状の変化率を測定した時点で水分率計を用いて測定した複数の酒米の吸水率の平均値を算出することで、酒米の3次元形状の変化率と吸水率との関係をプロットした(図5参照)。図5に示すように、2次元形状の変化率と吸水率とに関する近似線の決定係数が0.30であるのに対し、3次元形状の変化率と吸水率とに関する近似線の決定係数が0.91であり、このことから、2次元形状の変化率は3次元形状の変化率と比較して吸水率との相関が低いことが分かった。   The inventor of the present application has studied a method of determining a water absorption state such as a water absorption rate and a water absorption distribution of sake rice by image observation. , A new finding that there is a correlation with the rate of change of the shape of sake rice in the thickness direction. Specifically, the inventor of the present application has determined that the rate of change of the two-dimensional shape in a plane orthogonal to the thickness direction of the sake rice measured based on the two-dimensional image data for a plurality of sake rice arranged in a state of being immersed in water. By calculating the average value and calculating the average value of the water absorption rates of a plurality of sake rice measured using a moisture meter at the time when the change rate of the two-dimensional shape is measured, the change in the two-dimensional shape of the sake rice is calculated. The relationship between the percentage and the water absorption was plotted. In addition, for a plurality of sake rice arranged in a state of being immersed in water, the average value of the change rate of the three-dimensional shape measured using a known method is calculated, and at the time when the change rate of the three-dimensional shape is measured. The relationship between the change rate of the three-dimensional shape of the sake rice and the water absorption was plotted by calculating the average value of the water absorption of a plurality of sake rice measured using a moisture meter (see FIG. 5). As shown in FIG. 5, the coefficient of determination of the approximate line for the change rate of the two-dimensional shape and the water absorption is 0.30, whereas the coefficient of determination of the approximate line for the change rate of the three-dimensional shape and the water absorption is 0.91, which indicates that the change rate of the two-dimensional shape has a lower correlation with the water absorption rate than the change rate of the three-dimensional shape.

そして、本願発明者は、2次元形状の変化率及び3次元形状の変化率を基に、厚み方向における形状の変化率(厚み変化率)を算出し、この厚み方向における形状の変化率に対して、横割れ率及び縦割れ率をそれぞれプロットし(図6参照)、厚み方向における形状の変化率に対する横割れ率及び縦割れ率の相関を調べた。その結果、図6に示すように、縦割れ率と厚み方向における形状の変化率との間に高い相関関係があることが分かった。したがって、複数の品種について、上記と同様の手法により、縦割れ率と厚み方向における形状の変化率との相関関係に関するデータベースを予め作成して、適宜記憶装置35に記憶させておくことで、この相関関係を基に、縦割れ率から厚み方向における形状の変化率を推定することができる。   The inventor of the present application calculates the shape change rate in the thickness direction (thickness change rate) based on the two-dimensional shape change rate and the three-dimensional shape change rate, and calculates the shape change rate in the thickness direction. Then, the lateral cracking ratio and the vertical cracking ratio were plotted (see FIG. 6), and the correlation between the lateral cracking ratio and the vertical cracking ratio with respect to the shape change rate in the thickness direction was examined. As a result, as shown in FIG. 6, it was found that there was a high correlation between the vertical cracking rate and the shape change rate in the thickness direction. Therefore, for a plurality of varieties, a database relating to the correlation between the vertical cracking rate and the shape change rate in the thickness direction is created in advance by the same method as described above, and stored in the storage device 35 as appropriate. Based on the correlation, the change rate of the shape in the thickness direction can be estimated from the vertical crack rate.

そこで、図2の工程#16において、推定部19は、割れ率算出部16によって算出された縦割れ率、及び上述した縦割れ率と厚み方向における形状の変化率との相関関係とを基に、酒米の厚み方向における形状の変化率を推定する。このようにして、推定部19によって推定された厚み方向における形状の変化率が変化率推定装置20Aにより推定された厚み方向における形状の変化率に相当する。このようにして推定された厚み方向における形状の変化率を、記憶装置35に記憶すること、及び、出力装置3から出力することができる。また、厚み方向における形状の変化率の時間的な変動を監視して、厚み方向における形状の変化率が所定の値になった時点で出力装置3から出力することもできる。   Therefore, in step # 16 of FIG. 2, the estimating unit 19 determines the vertical cracking rate calculated by the cracking rate calculating unit 16 and the above-described correlation between the vertical cracking rate and the shape change rate in the thickness direction. Then, the change rate of the shape of the sake rice in the thickness direction is estimated. Thus, the change rate of the shape in the thickness direction estimated by the estimation unit 19 corresponds to the change rate of the shape in the thickness direction estimated by the change rate estimation device 20A. The rate of change of the shape in the thickness direction estimated in this way can be stored in the storage device 35 and output from the output device 3. Further, it is also possible to monitor the temporal change of the shape change rate in the thickness direction and to output from the output device 3 when the shape change rate in the thickness direction becomes a predetermined value.

状態算出部30は、推定部19で推定された厚み方向における形状の変化率(変化率推定装置20Aで推定された厚み方向における形状の変化率)を基に、酒米の比重、密度又は吸水率を算出する。尚、本実施形態では、酒米の比重を算出するものとする。   The state calculation unit 30 determines the specific gravity, density, or water absorption of the sake rice based on the shape change rate in the thickness direction estimated by the estimation unit 19 (the shape change rate in the thickness direction estimated by the change rate estimation device 20A). Calculate the rate. In the present embodiment, the specific gravity of sake rice is calculated.

具体的に、図2の工程#17において、状態算出部30は、面積情報取得部13から酒米の面積を取得するとともに、推定部19から酒米の厚み方向における形状の変化率を取得し、更に、記憶装置35に記憶されている、水に浸漬する前の酒米の厚みを取得して、これら取得した厚み方向における形状の変化率、面積、厚みを掛け合わせて酒米の体積を算出して、算出した体積を、別途計測して適宜入力される酒米の重量で割ることによって、酒米の比重を算出できる。算出された酒米の比重を、記憶装置35に記憶すること、及び出力装置3から出力することができる。   Specifically, in step # 17 of FIG. 2, the state calculation unit 30 acquires the area of the sake rice from the area information acquisition unit 13 and acquires the rate of change in the shape of the sake rice in the thickness direction from the estimation unit 19. Further, the thickness of the sake rice before immersion in water stored in the storage device 35 is obtained, and the volume of the sake rice is multiplied by the obtained rate of change in shape, area, and thickness in the thickness direction. The specific gravity of the sake rice can be calculated by dividing the calculated volume by the weight of the sake rice that is separately measured and input as appropriate. The calculated specific gravity of the sake rice can be stored in the storage device 35 and output from the output device 3.

図7は、変化率推定方法及び状態推定方法の一部である機械学習処理を説明するフローチャートである。工程#20において画像データ取得部11は、酒米の2次元画像データを取得する。つまり、水に浸漬した状態に配置された酒米を撮影した2次元画像データを取得する画像データ取得工程が実行される。そして、工程#21において画像処理部12は、画像処理を施して、処理済みの2次元画像データを生成する。その後、工程#22において形態判定部17は、画像処理部12が生成した処理済みの2次元画像データとその処理済みの2次元画像データに含まれる酒米の形態に関する情報との組み合わせで構成される訓練データを訓練データ記憶部35cに記憶し、工程#23においてその訓練データの機械学習を実行する。例えば、一つの処理済みの2次元画像データに対して、その処理済みの2次元画像データに含まれている一つの酒米の形態について、人間が判定した正解データがラベル付けされることで、一つの訓練データが作成される。訓練データ記憶部35cには、このようにして作成した複数の訓練データが記憶されている。   FIG. 7 is a flowchart illustrating a machine learning process that is part of the change rate estimation method and the state estimation method. In step # 20, the image data acquisition unit 11 acquires two-dimensional image data of sake rice. That is, an image data acquisition step of acquiring two-dimensional image data of the sake rice placed in the state of being immersed in water is executed. Then, in step # 21, the image processing unit 12 performs image processing to generate processed two-dimensional image data. Thereafter, in step # 22, the morphological determination unit 17 is configured by a combination of the processed two-dimensional image data generated by the image processing unit 12 and information on the morphology of the sake rice contained in the processed two-dimensional image data. The training data is stored in the training data storage unit 35c, and machine learning of the training data is executed in step # 23. For example, for one processed two-dimensional image data, the correct data determined by a human is labeled with respect to one form of sake rice included in the processed two-dimensional image data, One training data is created. The training data storage unit 35c stores a plurality of training data created in this way.

本実施形態では、形態判定部17で実行されるCNNの出力層ではソフトマックス関数が尤度関数として用いられる。ソフトマックス関数は、出力層の全ユニットの値を累積した値で各ユニットの値を除算して正規化するので、出力層の各ユニットの尤度は0〜1の間の値をとる。そして、形態判定部17は、出力層で生成される推定データ(「尤度」)と正解データ(「1」)との誤差に対して誤差逆伝播法を用いてその誤差を小さくするように、実行するCNNの複数の重みパラメータを修正して、画像処理部12が生成した処理済みの2次元画像データが入力されると酒米の形態に関する情報が出力されるアルゴリズムを機械学習の実行結果として得る。   In the present embodiment, a softmax function is used as a likelihood function in the output layer of the CNN executed by the form determining unit 17. Since the softmax function divides the value of each unit by the value obtained by accumulating the values of all units in the output layer and normalizes it, the likelihood of each unit in the output layer takes a value between 0 and 1. Then, the form determination unit 17 reduces the error between the estimated data (“likelihood”) generated in the output layer and the correct data (“1”) by using an error back propagation method. The algorithm that corrects a plurality of weight parameters of the CNN to be executed and outputs information on the morphology of the sake rice when the processed two-dimensional image data generated by the image processing unit 12 is input is executed. Get as.

このように、機械学習は、酒米の形態を判定するのに適した状態にした2次元画像データ(処理済みの2次元画像データ)を訓練データとして行われ、上述した酒米の形態の判定も、酒米の形態を判定するのに適した状態にした2次元画像データ(処理済みの2次元画像データ)に対して行われる。その結果、形態判定部17の判定結果が適切である可能性が高まる。   As described above, the machine learning is performed by using the two-dimensional image data (processed two-dimensional image data) in a state suitable for determining the morphology of the sake rice as the training data, and determining the morphology of the sake rice described above. This is also performed on the two-dimensional image data (processed two-dimensional image data) that is in a state suitable for determining the form of sake rice. As a result, the possibility that the determination result of the form determination unit 17 is appropriate is increased.

次に、本実施形態のような処理済みの2次元画像データを訓練データとして機械学習を行うことにより得られる効果を説明する。図8〜図10は、学習回数と正解率との関係を示すグラフである。具体的には、図8は、図3の2次元画像データBに正解データがラベル付けされた訓練データを用いて機械学習を行った場合の学習回数と正解率の関係を示すグラフである。図9は、図3の2次元画像データDを訓練データとして機械学習を行った場合の学習回数と正解率の関係を示すグラフである。図10は、図3の2次元画像データFを訓練データとして機械学習を行った場合の学習回数と正解率の関係を示すグラフである。   Next, effects obtained by performing machine learning using the processed two-dimensional image data as the training data as in the present embodiment will be described. 8 to 10 are graphs showing the relationship between the number of times of learning and the correct answer rate. Specifically, FIG. 8 is a graph showing the relationship between the number of times of learning and the correct answer rate when machine learning is performed using training data in which correct answer data is labeled on the two-dimensional image data B in FIG. FIG. 9 is a graph showing the relationship between the number of times of learning and the accuracy rate when machine learning is performed using the two-dimensional image data D of FIG. 3 as training data. FIG. 10 is a graph showing the relationship between the number of times of learning and the accuracy rate when machine learning is performed using the two-dimensional image data F of FIG. 3 as training data.

図8に示すように、2次元画像データBを訓練データとした場合、学習回数が2万回に達しても、正解率は80%未満である。この80%という正解率は、例えば、ある程度熟練した人間が2次元画像データから酒米の形態を判定する場合の正解率に対応する。それに対して、2次元画像データDに正解データがラベル付けされた訓練データを用いて機械学習を行った場合、学習回数が約12000回程度で正解率が80%を上回る。また、2次元画像データFに正解データがラベル付けされた訓練データを用いて機械学習を行った場合、学習回数が約5000回程度で正解率が80%を上回る。このように、画像処理を行った後の処理済みの2次元画像データ、即ち、酒米の形態を判定するのに適した状態にした2次元画像データを訓練データとすることで、学習回数が少なくても正解率が上昇するという効果を得ることができた。   As shown in FIG. 8, when the two-dimensional image data B is used as training data, the correct answer rate is less than 80% even if the number of times of learning reaches 20,000. The correct answer rate of 80% corresponds to, for example, the correct answer rate in the case where a somewhat skilled person determines the form of sake rice from two-dimensional image data. On the other hand, when machine learning is performed using training data in which the correct answer data is labeled on the two-dimensional image data D, the number of times of learning is about 12,000, and the correct answer rate exceeds 80%. When machine learning is performed using training data in which the correct data is labeled on the two-dimensional image data F, the number of times of learning is about 5,000 and the correct answer rate exceeds 80%. As described above, the processed two-dimensional image data after the image processing is performed, that is, the two-dimensional image data in a state suitable for determining the form of sake rice is used as training data, so that the number of times of learning is reduced. At least, the effect of increasing the correct answer rate was obtained.

以上のように、変化率推定装置20A及び変化率推定方法によれば、水に浸漬した状態に配置された複数の酒米の2次元画像データを参照して、複数の酒米のうち、縦割れが発生している酒米の割合である縦割れ率を算出した上で、この算出した縦割れ率、及び予め定められた酒米の縦割れ率と厚み方向における形状の変化率との相関関係を基にして、酒米の厚み方向における形状の変化率を推定することができる。したがって、複数の撮影装置2などを備えるような複雑な構成の装置を要することなく、容易に、酒米の厚み方向における形状の変化率を推定することができる。尚、変化率推定装置20A及び変化率推定方法を使用して、吸水によって主に鉛直方向に直交する方向に膨張する炊飯対象米について、厚み方向における形状の変化率を推定することも当然に可能である。   As described above, according to the change rate estimating apparatus 20A and the change rate estimating method, the two-dimensional image data of the plurality of sake rice arranged in the state of being immersed in water is referred to, and among the plurality of sake rice, After calculating the vertical cracking rate, which is the percentage of cracked rice, the calculated vertical cracking rate, and the correlation between the predetermined vertical cracking rate of sake rice and the rate of change in shape in the thickness direction. Based on the relationship, the change rate of the shape of the sake rice in the thickness direction can be estimated. Therefore, the change rate of the shape of the sake rice in the thickness direction can be easily estimated without the need for a device having a complicated configuration including a plurality of photographing devices 2 and the like. In addition, it is naturally possible to estimate the change rate of the shape in the thickness direction of the rice to be cooked, which expands mainly in the direction perpendicular to the vertical direction due to water absorption, using the change rate estimation device 20A and the change rate estimation method. It is.

また、状態推定装置1A及び状態推定方法によれば、上記推定された厚み方向における形状の変化率を基にして酒米の比重を算出することができ、酒米の比重を容易に推定することができる。酒米などを含む一般的な米の比重は、米中のデンプンの粗密に大きく関わっており、心白のようなデンプンが疎な部分が多いほど、比重は大きくなる。したがって、比重を推定することで、心白の状態(例えば、大きさ)を評価することが可能となる。そして、例えば、酒造に好適な酒米は、粒が大きく、心白も大きいこととされているが、心白が大きすぎると精米時の割れが増えるため、心白の大きさは一定の範囲内であることが好ましいと考えられるため、推定した比重を基に心白の状態を評価して、この評価結果を利用して酒造に好適な米の選別を行うこともできる。   Further, according to the state estimation device 1A and the state estimation method, it is possible to calculate the specific gravity of sake rice based on the estimated change rate of the shape in the thickness direction, and to easily estimate the specific gravity of sake rice. Can be. The specific gravity of general rice, including sake rice, is largely related to the density of the starch in the rice, and the more sparse the starch, such as heart white, the greater the specific gravity. Therefore, by estimating the specific gravity, it is possible to evaluate the state of white heart (for example, the size). And, for example, sake rice suitable for sake brewing is said to have a large grain and a large white heart, but if the white heart is too large, cracking during rice polishing will increase, so the size of the white heart is within a certain range Since it is considered preferable to be within the range, it is also possible to evaluate the whiteness of the heart based on the estimated specific gravity, and to use the evaluation result to select rice suitable for sake brewing.

[第2実施形態]
第2実施形態の状態推定装置1B(1)は機械学習処理の内容が上記実施形態と異なっている。以下に第2実施形態の状態推定装置1Bについて説明するが、上記実施形態と同様の構成については説明を省略する。
[Second embodiment]
The state estimating device 1B (1) of the second embodiment differs from the above embodiment in the content of the machine learning process. Hereinafter, a state estimation device 1B according to the second embodiment will be described, but the description of the same configuration as the above embodiment will be omitted.

図11は第2実施形態の状態推定装置1Bの構成を示す図である。図示するように、第2実施形態の状態推定装置1Bは、形態判定部17が判定した酒米の形態の妥当性の高低を判定する妥当性判定部36を更に備える。この妥当性判定部36は米分析装置10B(10)の割れ率算出部16に設けられている。そして、形態判定部17は、妥当性判定部36によって妥当性が高いと判定された酒米の形態に関する情報とその形態が判定された処理済みの2次元画像データとの組み合わせを既存の訓練データに追加して機械学習を実行する。   FIG. 11 is a diagram illustrating a configuration of a state estimation device 1B according to the second embodiment. As shown in the figure, the state estimation device 1B of the second embodiment further includes a validity determining unit 36 that determines the validity of the rice rice form determined by the form determining unit 17. This validity determination unit 36 is provided in the crack rate calculation unit 16 of the rice analyzer 10B (10). Then, the morphological determination unit 17 compares the combination of the information on the morphology of the sake rice determined to have high validity by the validity determination unit 36 and the processed two-dimensional image data whose morphology has been determined to the existing training data. And perform machine learning.

本実施形態でも、米分析装置10Bは、情報の演算処理機能及び情報の入出力機能及び情報の記憶機能などを備える1台又は複数台のコンピュータ装置などを用いて実現される。その場合、画像データ取得部11の機能と、画像処理部12の機能と、面積情報取得部13の機能と、割れ率算出部16(形態判定部17の機能、形態データ解析部18の機能)の機能と、推定部19の機能と、状態算出部30の機能をコンピュータ装置に実現させるプログラムを、そのコンピュータ装置にインストールしておけばよい。   Also in the present embodiment, the rice analyzer 10B is realized using one or a plurality of computer devices having an information processing function, an information input / output function, and an information storage function. In that case, the function of the image data acquisition unit 11, the function of the image processing unit 12, the function of the area information acquisition unit 13, and the crack rate calculation unit 16 (the function of the form determination unit 17, the function of the form data analysis unit 18) A program that causes a computer device to implement the functions of the above, the function of the estimation unit 19, and the function of the state calculation unit 30 may be installed in the computer device.

形態記憶部35dは、形態判定部17の判定結果を、水への浸漬を開始した後での一つの酒米の形態の変化履歴が分かる状態で記憶する。例えば、形態記憶部35dには、形態の判定対象とする複数の酒米のそれぞれについて、水への浸漬を開始した時点から現在(最新)に至るまでの、形態判定部17の判定結果を、その形態の変化履歴が分かる状態で記憶している。その結果、一つの酒米の形態がどのような変化履歴を示したのかを知ることができる。   The form storage unit 35d stores the determination result of the form determination unit 17 in a state where the change history of the form of one sake rice after the start of immersion in water is known. For example, in the form storage unit 35d, the determination result of the form determination unit 17 from the time when the immersion in water is started to the present (latest) for each of a plurality of sake rice to be determined for the form is The change history of the form is stored in a state where it can be understood. As a result, it is possible to know what change history the form of one sake rice shows.

図12は、第2実施形態の機械学習処理を説明するフローチャートである。この機械学習処理は、第1実施形態で説明した機械学習処理(図7)に加えて行われる。具体的に説明すると、工程#30において妥当性判定部36は、形態判定部17が判定した酒米の現在の形態の妥当性を判定する。つまり、形態判定部17が実行する形態判定工程で判定した酒米の形態の妥当性の高低を判定する妥当性判定工程(割れ率算出工程)が実行される。図13は、妥当性判定部36が行う妥当性判定処理を説明するフローチャートである。工程#40において妥当性判定部36は、形態判定部17が判定した酒米の現在の形態の確信度が所定値以上であるか否かを判定する。例えば、形態判定部17の出力層で用いられるソフトマックス関数により算出される最大の尤度を上記確信度とした場合、そのソフトマックス関数により算出される最大の尤度が所定値以上(例えば0.8以上など)であれば、妥当性判定部36は、酒米の現在の形態の確信度が所定値以上であると判定し、工程#41に移行する。   FIG. 12 is a flowchart illustrating a machine learning process according to the second embodiment. This machine learning process is performed in addition to the machine learning process (FIG. 7) described in the first embodiment. More specifically, in step # 30, the validity determining unit 36 determines the validity of the current form of the sake rice determined by the form determining unit 17. That is, a validity determination step (crack rate calculation step) of determining the validity of the form of the sake rice determined in the form determination step performed by the form determination unit 17 is performed. FIG. 13 is a flowchart illustrating a validity determination process performed by the validity determination unit 36. In step # 40, the validity determination unit 36 determines whether the certainty factor of the current form of the sake rice determined by the form determination unit 17 is equal to or greater than a predetermined value. For example, when the maximum likelihood calculated by the softmax function used in the output layer of the form determination unit 17 is set as the certainty, the maximum likelihood calculated by the softmax function is equal to or more than a predetermined value (for example, 0). If it is .8 or more, the validity determination unit 36 determines that the certainty factor of the current form of the sake rice is equal to or greater than a predetermined value, and proceeds to step # 41.

次に工程#41において妥当性判定部36は、酒米の形態の変化履歴は所定の基準に適合したものであるか否かを判定する。この所定の基準は以下の6パターンである。そして、妥当性判定部36は、酒米の形態の変化履歴が6パターンの何れかの基準に適合していれば、工程#42に移行して「妥当性は高い」と判定する。   Next, in step # 41, the validity determination unit 36 determines whether or not the change history of the form of the sake rice conforms to a predetermined standard. The predetermined criterion is the following six patterns. If the change history of the shape of the sake rice conforms to any of the six patterns, the validity determination unit 36 proceeds to step # 42 and determines that the validity is high.

具体的に説明すると、例えば一つの酒米の形態を継続的に見た場合、酒米に割れが生じた形態から、酒米が損傷していない形態に戻ることはないし、酒米に亀裂が生じた形態に戻ることもない。つまり、水への浸漬を開始した後での酒米の形態が形態判定部17によって適切に判定されていれば、その変化履歴は所定の変化基準に適合しているはずである。そこで本特徴構成では、妥当性判定部36は、形態判定部17が判定した酒米の形態の確信度が所定値以上であり、且つ、水への浸漬を開始した後での酒米の形態の変化履歴が所定の変化基準に適合している場合には、形態判定部17が判定した酒米の形態の妥当性が高いと判定する。その結果、適切な訓練データの蓄積と機械学習とを自動で行うことができる。   More specifically, for example, when one sake rice is continuously viewed, the sake rice does not return from the cracked form to the non-damaged form, and cracks are formed in the sake rice. There is no return to the form that occurred. That is, if the morphology of the sake rice after the start of immersion in water is appropriately determined by the morphology determination unit 17, the change history should conform to the predetermined change criterion. Therefore, in the present characteristic configuration, the validity determination unit 36 determines that the degree of certainty of the morphology of the sake rice determined by the morphology determination unit 17 is equal to or more than a predetermined value, and that the morphology of the sake rice after starting immersion in water If the change history conforms to a predetermined change criterion, it is determined that the form of the sake rice determined by the form determining unit 17 is highly appropriate. As a result, accumulation of appropriate training data and machine learning can be performed automatically.

それに対して、妥当性判定部36は、工程#40において確信度が所定値以上ではないと判定した場合、及び、工程#41において酒米の形態の変化履歴が所定の基準に適合したものでないと判定した場合には、工程#43に移行して「妥当性は低い」と判定する。   On the other hand, the validity determination unit 36 determines that the certainty factor is not equal to or more than the predetermined value in step # 40, and that the change history of the sake rice form does not conform to the predetermined standard in step # 41. When the determination is made, the process proceeds to step # 43, where it is determined that “the validity is low”.

以上のように、妥当性判定部36は、形態判定部17が判定した酒米の形態の確信度が所定値以上であり、且つ、水への浸漬を開始した後での酒米の形態の変化履歴が所定の変化基準に適合している場合には、形態判定部17が判定した酒米の形態の妥当性が高いと判定し、形態判定部17が判定した酒米の形態の確信度が所定値未満である場合、又は、水への浸漬を開始した後での酒米の形態の変化履歴が所定の変化基準に適合していない場合には、形態判定部17が判定した酒米の形態の妥当性が低いと判定する。   As described above, the validity determination unit 36 determines that the degree of certainty of the morphology of the sake rice determined by the morphology determination unit 17 is equal to or greater than the predetermined value, and that the morphology of the sake rice after the start of immersion in water. If the change history meets the predetermined change criterion, it is determined that the form of the sake rice determined by the form determination unit 17 is high in validity, and the certainty factor of the form of the sake rice determined by the form determination unit 17 Is less than a predetermined value, or when the change history of the morphology of the sake rice after the start of immersion in water does not meet the predetermined change criteria, the sake rice determined by the form determination unit 17 Is determined to be low in validity.

次に、工程#31において、妥当性が高いと判定された酒米の形態に関する情報と、処理済みの2次元画像データとの組み合わせを訓練データとする。そして、工程#32において、その訓練データを既存の訓練データに追加して機械学習を実行する。   Next, in step # 31, the combination of the information on the morphology of the sake rice determined to be highly relevant and the processed two-dimensional image data is used as training data. Then, in step # 32, the training data is added to the existing training data to execute machine learning.

このように、妥当性判定部36は、形態判定部17が判定した酒米の形態の妥当性の高低を判定し、形態判定部17は、妥当性判定部36によって妥当性が高いと判定された酒米の形態に関する情報とその形態が判定された処理済みの2次元画像データとの組み合わせを既存の訓練データに追加して機械学習を実行する。つまり、適切な訓練データの蓄積と機械学習とを自動で行うことができる。   As described above, the validity determination unit 36 determines the level of validity of the sake rice form determined by the form determination unit 17, and the form determination unit 17 determines that the validity is high by the validity determination unit 36. Machine learning is performed by adding a combination of the information on the shape of the sake rice and the processed two-dimensional image data in which the shape is determined to the existing training data. That is, accumulation of appropriate training data and machine learning can be automatically performed.

[第3実施形態]
第3実施形態の状態推定装置1C(1)は、推定部19において3次元形状の変化率を推定する点が上記実施形態と異なっている。以下に第3実施形態の状態推定装置1Cについて説明するが、上記実施形態と同様の構成については説明を省略する。
[Third embodiment]
The state estimating device 1C (1) of the third embodiment differs from the above embodiment in that the estimating unit 19 estimates the rate of change of the three-dimensional shape. Hereinafter, a state estimating device 1C according to the third embodiment will be described, but the description of the same configuration as the above embodiment will be omitted.

図14は第3実施形態の状態推定装置1Cの構成を示す図である。図示するように、第3実施形態の状態推定装置1Cは、米分析装置10Cに設けられた2次元形状変化率算出部(2次元形状変化率算出手段)40を備える。この2次元形状変化率算出部40は、画像データ取得部11が取得した2次元画像データを参照して、当該2次元画像データを解析し、酒米の厚み方向と直交する平面に沿った方向における酒米の2次元形状の変化率を算出する。   FIG. 14 is a diagram illustrating a configuration of a state estimation device 1C according to the third embodiment. As illustrated, the state estimation device 1C of the third embodiment includes a two-dimensional shape change rate calculation unit (two-dimensional shape change rate calculation means) 40 provided in the rice analysis device 10C. The two-dimensional shape change rate calculating unit 40 analyzes the two-dimensional image data with reference to the two-dimensional image data acquired by the image data acquiring unit 11, and calculates a direction along a plane orthogonal to the thickness direction of the sake rice. The change rate of the two-dimensional shape of the sake rice in is calculated.

したがって、本実施形態において、米分析装置10Cは、情報の演算処理機能及び情報の入出力機能及び情報の記憶機能などを備える1台又は複数台のコンピュータ装置などを用いて実現され、その場合、画像データ取得部11の機能と、画像処理部12の機能と、2次元形状変化率算出部40(面積判定部41の機能、面積データ解析部42の機能)の機能と、推定部19の機能と、状態算出部30の機能をコンピュータ装置に実現させるプログラムを、そのコンピュータ装置にインストールしておけばよい。   Therefore, in the present embodiment, the rice analyzer 10C is realized using one or more computer devices having an information processing function, an information input / output function, an information storage function, and the like. The function of the image data acquisition unit 11, the function of the image processing unit 12, the function of the two-dimensional shape change rate calculation unit 40 (the function of the area determination unit 41, the function of the area data analysis unit 42), and the function of the estimation unit 19 A program for causing the computer device to realize the function of the state calculation unit 30 may be installed in the computer device.

本実施形態において、2次元形状変化率算出部40は、面積判定部41と面積データ解析部42とを備えており、面積判定部41は、画像データ取得部11が取得し、画像処理部12が処理した2次元画像データを参照して、酒米の厚み方向と直交する平面に沿った方向における酒米の面積を、様々な画像認識技術を用いて自動で判定することができる。   In the present embodiment, the two-dimensional shape change rate calculation unit 40 includes an area determination unit 41 and an area data analysis unit 42. The area determination unit 41 acquires the image data With reference to the processed two-dimensional image data, the area of the sake rice in a direction along a plane orthogonal to the thickness direction of the sake rice can be automatically determined using various image recognition techniques.

上記のように、面積判定部41が酒米の面積を判定するために参照する入力データとしての2次元画像データは、画像処理部12が、画像データ取得部11が取得した画像データにおいて、1つの酒米のみを含む領域の2次元酒米画像データを抽出する処理、1つの酒米が所定の姿勢になるように酒米の2次元画像を回転させる処理、及び画像のコントラストを調整する処理を含む画像処理を行った後の処理済みの2次元画像データである。つまり、面積判定部41は、酒米の面積を判定するのに適した状態にした2次元画像データ(処理済みの2次元画像データ)から酒米の面積を判定できるので、その判定結果が適切である可能性が高まる。   As described above, the two-dimensional image data as the input data referred to by the area determination unit 41 to determine the area of the sake rice is one in the image data acquired by the image data acquisition unit 11 by the image processing unit 12. A process of extracting two-dimensional sake rice image data of an area including only one sake rice, a process of rotating a two-dimensional image of sake rice so that one sake rice has a predetermined posture, and a process of adjusting image contrast Is processed two-dimensional image data after image processing including. That is, since the area determination unit 41 can determine the area of the sake rice from the two-dimensional image data (processed two-dimensional image data) that is in a state suitable for determining the area of the sake rice, the determination result is appropriate. Is more likely to be.

また、本実施形態において、記憶装置35には、水に浸漬する前の酒米の厚み及び厚み方向と直交する平面に沿った面積が予め記憶されている。   In the present embodiment, the thickness of the sake rice before immersion in water and the area along a plane orthogonal to the thickness direction are stored in the storage device 35 in advance.

図15は、本実施形態に係る変化率推定方法及び状態推定方法に含まれる、酒米の2次元形状の変化率を算出する工程(2次元形状変化率算出工程)、2次元形状の変化率を3次元形状の変化率に補正して3次元形状の変化率を推定する工程(推定工程)、及び推定した3次元形状の変化率を基に、酒米の比重を算出する工程(状態算出工程)を説明するフローチャートである。図15の工程#50から工程#52は、図2の工程#10から工程#12と同様であるため、詳細については説明を省略する。   FIG. 15 shows a process of calculating the change rate of the two-dimensional shape of sake rice (two-dimensional shape change rate calculating step), which is included in the change rate estimation method and the state estimation method according to the present embodiment. Estimating the rate of change of the three-dimensional shape by correcting the three-dimensional shape to the rate of change of the three-dimensional shape (estimating step), and calculating the specific gravity of sake rice based on the estimated rate of change of the three-dimensional shape (state calculation) 3 is a flowchart illustrating a step). Steps # 50 to # 52 in FIG. 15 are the same as steps # 10 to # 12 in FIG.

図15の工程#53において、面積判定部41は、酒米の面積を判定する。つまり、画像処理工程で生成した処理済みの2次元画像データにおける1つの酒米の面積を判定する面積判定工程(2次元形状変化率算出工程)が実行される。本実施形態では、面積判定部41が生成した処理済みの2次元画像データ(図3の2次元画像データF)における1つの酒米の面積を判定する。尚、水に浸漬された状態にある酒米が複数ある場合、面積判定部41がそのうちのいくつの酒米の面積を判定するのかは適宜設定可能である。そして、面積判定部41による判定結果は、複数の酒米のそれぞれに付与された識別子と関連付けられて酒米ごとに面積記憶部35bに記憶される。面積記憶部35bには、面積判定部41の判定結果が、水への浸漬を開始した後での1つの酒米の面積の変化履歴が分かる状態で記憶される。   In step # 53 of FIG. 15, the area determining unit 41 determines the area of the sake rice. That is, an area determining step (a two-dimensional shape change rate calculating step) of determining an area of one sake rice in the processed two-dimensional image data generated in the image processing step is executed. In the present embodiment, the area of one sake rice in the processed two-dimensional image data (two-dimensional image data F in FIG. 3) generated by the area determination unit 41 is determined. When there are a plurality of sake rice immersed in water, how many of the sake rice the area determination unit 41 determines the area of the sake rice can be appropriately set. The determination result by the area determination unit 41 is stored in the area storage unit 35b for each sake rice in association with the identifier given to each of the plurality of sake rice. The determination result of the area determination unit 41 is stored in the area storage unit 35b in a state where the change history of the area of one sake rice after the start of immersion in water is known.

このように、本実施形態では、面積判定部41が、画像処理部12が生成した処理済みの2次元画像データにおける1つの酒米の面積を判定する。その結果、酒米の面積に関して、短時間で判定結果を得ることができる。   As described above, in the present embodiment, the area determination unit 41 determines the area of one sake rice in the processed two-dimensional image data generated by the image processing unit 12. As a result, a determination result can be obtained in a short time for the area of sake rice.

図15の工程#54において、面積データ解析部42は、画像データ取得部11が取得した2次元画像データに含まれる複数の酒米のそれぞれについての面積判定部41の判定結果に基づいて酒米の2次元形状の変化率を算出する。つまり、面積判定工程の判定結果に基づいて2次元形状の変化率を算出する面積データ解析工程(2次元形状変化率算出工程)が実行される。面積記憶部35bには、複数の酒米のそれぞれについての面積判定部41の判定結果(即ち、酒米ごとの面積)が、水への浸漬を開始した後での1つの酒米の面積の変化履歴が分かる状態で記憶されている。また、面積記憶部35bには、各酒米について、水への浸漬を開始した時点から現在(最新)までの少なくとも2以上の面積が変化履歴として記憶されている。そして、面積データ解析部42は、面積記憶部35bに記憶されている、水への浸漬を開始した時点の複数の酒米の面積及び現在(最新)の複数の酒米の面積を読み出して、各酒米に関するこれら2つの面積を基に、各酒米ごとの2次元形状の変化率を算出する。   In step # 54 of FIG. 15, the area data analysis unit 42 determines whether or not the sake rice is based on the determination result of the area determination unit 41 for each of the plurality of sake rice included in the two-dimensional image data acquired by the image data acquisition unit 11. Of the two-dimensional shape are calculated. That is, an area data analysis step (two-dimensional shape change rate calculation step) of calculating a change rate of the two-dimensional shape based on the determination result of the area determination step is executed. In the area storage unit 35b, the determination result of the area determination unit 41 for each of the plurality of sake rice (that is, the area for each sake rice) is stored as the area of one sake rice after the start of immersion in water. The change history is stored in a state where the change history can be understood. Further, in the area storage unit 35b, for each sake rice, at least two or more areas from the start of immersion in water to the present (latest) are stored as a change history. Then, the area data analysis unit 42 reads the area of the plurality of sake rice at the time of starting the immersion in water and the area of the current (latest) plurality of sake rice stored in the area storage unit 35b, Based on these two areas for each sake rice, the change rate of the two-dimensional shape of each sake rice is calculated.

尚、面積判定部41の判定結果は、出力装置3から出力することもでき、また、酒米の2次元形状の変化率を随時算出して出力装置3からリアルタイムで出力することもできる。更に、面積記憶部35bに、面積判定部41の判定結果が、水への浸漬を開始した後での1つの酒米の面積の変化履歴が分かる状態で記憶されている場合には、面積や2次元形状の変化率の時間的な変化など算出して出力装置3から出力することもできる。   The determination result of the area determination unit 41 can be output from the output device 3, or the change rate of the two-dimensional shape of the sake rice can be calculated at any time and output from the output device 3 in real time. Furthermore, when the determination result of the area determination unit 41 is stored in the area storage unit 35b in a state where the change history of the area of one sake rice after the start of immersion in water is known, the area and the area are determined. The change rate of the two-dimensional shape over time can be calculated and output from the output device 3.

本実施形態に係る推定部19は、2次元形状変化率算出部40が算出した酒米の2次元形状の変化率を取得し、2次元形状の変化率、及び予め定められた縦割れ率と厚み方向における形状の変化率との相関関係を基に、2次元形状の変化率を補正して3次元形状の変化率とする点が上記各実施形態とは異なっている。   The estimating unit 19 according to the present embodiment acquires the two-dimensional shape change rate of the sake rice calculated by the two-dimensional shape change rate calculating unit 40, and obtains the two-dimensional shape change rate and a predetermined vertical crack rate. The present embodiment is different from the above embodiments in that the change rate of the two-dimensional shape is corrected to the change rate of the three-dimensional shape based on the correlation with the change rate of the shape in the thickness direction.

上述したように、本願発明者は、2次元形状の変化率が3次元形状の変化率と比較して吸水率との相関が低いことを確認し、厚み方向における形状の変化率に対する横割れ率及び縦割れ率の相関を調べ、縦割れ率と厚み方向における形状の変化率との間に高い相関関係があることを確認した。そして、本願発明者は、縦割れ率が酒米の3次元形状の変化率に与える影響を示す指標たる縦割れ影響度を、上記相関関係に基づいて決定できることを見出した。本実施形態において、縦割れ影響度は、予め複数の品種に関する測定結果から得られたデータベースを基に、縦割れ率と厚み方向の形状の変化率との相関関係を示すグラフ(図6参照)を作成し、このグラフの傾きを縦割れ影響度として予め決定されるものである。   As described above, the inventor of the present application has confirmed that the rate of change of the two-dimensional shape has a lower correlation with the water absorption rate than the rate of change of the three-dimensional shape. The correlation between the vertical cracking rate and the vertical cracking rate was examined, and it was confirmed that there was a high correlation between the vertical cracking rate and the shape change rate in the thickness direction. And the inventor of the present application has found that the degree of vertical cracking, which is an index indicating the effect of the vertical cracking rate on the change rate of the three-dimensional shape of sake rice, can be determined based on the correlation. In the present embodiment, the degree of vertical cracking is a graph showing a correlation between the vertical cracking rate and the rate of change in shape in the thickness direction based on a database obtained in advance from measurement results for a plurality of types (see FIG. 6). And the inclination of this graph is determined in advance as the degree of vertical crack influence.

そこで、図15の#55において、推定部19は、2次元形状変化率算出部40によって算出された2次元形状の変化率に縦割れ影響度を乗じることで、2次元形状の変化率を補正して、3次元形状の変化率として推定する。このようにして、推定部19によって推定された3次元形状の変化率が変化率推定装置20Cにより推定された3次元形状の変化率に相当する。このようにして推定された3次元形状の変化率を、記憶装置35に記憶すること、及び、出力装置3から出力することができる。また、3次元形状の変化率の時間的な変動を監視して、3次元形状の変化率が所定の値になった時点で出力装置3から出力することもできる。   Therefore, in # 55 of FIG. 15, the estimation unit 19 corrects the change rate of the two-dimensional shape by multiplying the change rate of the two-dimensional shape calculated by the two-dimensional shape change rate calculation unit 40 by the degree of vertical crack influence. Then, the change rate of the three-dimensional shape is estimated. Thus, the change rate of the three-dimensional shape estimated by the estimating unit 19 corresponds to the change rate of the three-dimensional shape estimated by the change rate estimating device 20C. The change rate of the three-dimensional shape estimated in this manner can be stored in the storage device 35 and output from the output device 3. Further, it is also possible to monitor the temporal change of the change rate of the three-dimensional shape and output the change from the output device 3 when the change rate of the three-dimensional shape becomes a predetermined value.

本実施形態の状態算出部30は、推定部19で求められた3次元形状の変化率(変化率推定装置20Cで推定された3次元形状の変化率)を基に、酒米の比重、密度又は吸水率を算出する。尚、本実施形態では、酒米の比重を算出するものとする。   The state calculator 30 of the present embodiment calculates the specific gravity and density of sake rice based on the change rate of the three-dimensional shape obtained by the estimator 19 (the change rate of the three-dimensional shape estimated by the change rate estimating device 20C). Alternatively, the water absorption is calculated. In the present embodiment, the specific gravity of sake rice is calculated.

具体的に、図15の工程#56において、状態算出部30は、推定部19から酒米の3次元形状の変化率を取得し、更に、記憶装置35に予め記憶されている、水に浸漬する前の酒米の厚み及び厚み方向と直交する平面に沿った方向における面積を取得して、これら取得した3次元形状の変化率、並びに水に浸漬する前の厚み及び面積を掛け合わせて酒米の体積を算出して、算出した体積を、別途計測して適宜入力される酒米の重量で割ることによって、酒米の比重を算出できる。尚、算出された酒米の比重を、記憶装置35に記憶すること、及び出力装置3から出力することもできる。   Specifically, in step # 56 of FIG. 15, the state calculation unit 30 acquires the change rate of the three-dimensional shape of the sake rice from the estimation unit 19, and further immerses the rice in water stored in the storage device 35 in advance. The area of the sake rice in the direction along the plane perpendicular to the thickness direction and the thickness direction of the sake rice is obtained, and the rate of change of the obtained three-dimensional shape is multiplied by the thickness and area before immersion in water. By calculating the volume of rice and dividing the calculated volume by the weight of sake rice that is separately measured and appropriately input, the specific gravity of sake rice can be calculated. Note that the calculated specific gravity of the sake rice can be stored in the storage device 35 and output from the output device 3.

以上のように、変化率推定装置20C及び変化率推定方法によれば、水に浸漬した状態に配置された複数の酒米の2次元画像データを参照して、酒米の2次元形状の変化率を算出した上で、予め定められた酒米の縦割れ率と厚み方向における形状の変化率との相関関係により決定される縦割れ影響度を、2次元形状の変化率に乗じることにより、2次元形状の変化率を補正して3次元形状の変化率として推定することができる。したがって、吸水によって鉛直方向に直交する方向だけでなく、鉛直方向にも膨張する酒米について、複数の撮影装置2などを備えるような複雑な構成の装置を要することなく、容易に、2次元形状の変化率を補正して酒米の3次元形状の変化率を容易に推定することができる。尚、変化率推定装置20A及び変化率推定方法を使用して、吸水によって主に鉛直方向に直交する方向に膨張する炊飯対象米について、2次元形状の変化率を補正して3次元形状の変化率を求めることも当然に可能である。   As described above, according to the change rate estimating apparatus 20C and the change rate estimating method, the change in the two-dimensional shape of the sake rice is referred to by referring to the two-dimensional image data of the plurality of sake rice arranged in a state of being immersed in water. After calculating the rate, by multiplying the change rate of the two-dimensional shape by the vertical crack influence degree determined by the correlation between the predetermined vertical crack rate of sake rice and the change rate of the shape in the thickness direction, The change rate of the two-dimensional shape can be corrected and estimated as the change rate of the three-dimensional shape. Therefore, sake rice that expands not only in the direction perpendicular to the vertical direction but also in the vertical direction due to water absorption can be easily formed into a two-dimensional shape without the need for a device having a complicated configuration including a plurality of photographing devices 2 and the like. Is corrected, the change rate of the three-dimensional shape of the sake rice can be easily estimated. In addition, using the change rate estimating apparatus 20A and the change rate estimating method, the change rate of the two-dimensional shape is corrected by correcting the change rate of the two-dimensional shape for the rice to be cooked, which expands mainly in the direction perpendicular to the vertical direction due to water absorption. Obviously it is also possible to determine the rate.

また、状態推定装置1C及び状態推定方法によれば、上記推定された3次元形状の変化率を基にして酒米の比重を算出することができ、酒米の比重を容易に推定することができる。したがって、上述したように、比重を推定することで、心白の状態(例えば、大きさ)を評価することが可能となり、例えば、推定した比重を基に心白の状態を評価して、この評価結果を利用して酒造に好適な米の選別を行うこともできる。   Further, according to the state estimation device 1C and the state estimation method, the specific gravity of the sake rice can be calculated based on the estimated change rate of the three-dimensional shape, and the specific gravity of the sake rice can be easily estimated. it can. Therefore, as described above, by estimating the specific gravity, it is possible to evaluate the state of the white heart (for example, the size). For example, the state of the white heart is evaluated based on the estimated specific gravity. Utilizing the evaluation results, rice suitable for sake brewing can be selected.

[別実施形態]
[1]上記第1及び第2実施形態では、酒米について厚み方向における形状の変化率を推定し、上記第3実施形態では、酒米について3次元形状の変化率を推定するようにしているが、これに限られるものではなく、酒米以外の米、例えば、炊飯対象米などについて、厚み方向における形状の変化率や、3次元形状の変化率を推定することもできる。
[Another embodiment]
[1] In the first and second embodiments, the change rate of the shape of the sake rice in the thickness direction is estimated, and in the third embodiment, the change rate of the three-dimensional shape of the sake rice is estimated. However, the present invention is not limited to this, and it is also possible to estimate the rate of change in the shape in the thickness direction and the rate of change in the three-dimensional shape of rice other than sake rice, for example, rice to be cooked.

[2]上記各実施形態では、状態算出部30において、酒米の体積を算出して、算出した体積を別途計測して適宜入力される酒米の重量で割ることによって、酒米の比重を算出するようにしているが、酒米の重量を体積で割ることによって、酒米の密度を算出するようにしても良い。 [2] In each of the above embodiments, the specific gravity of the sake rice is calculated by calculating the volume of the sake rice in the state calculation unit 30 and separately measuring the calculated volume and dividing the calculated volume by the weight of the sake rice that is appropriately input. Although the calculation is performed, the density of the sake rice may be calculated by dividing the weight of the sake rice by the volume.

[3]上記第3実施形態において、状態算出部30は、推定した3次元形状の変化率と、予め記憶している酒米の3次元形状と酒米の吸水率との間の相関関係とに基づいて、酒米の含水率を算出することもできる。即ち、水に浸漬した状態に配置された複数の酒米について、予め既知の手法によって求めた3次元形状の変化率の平均値を算出し、3次元形状の変化率を求めるのに用いた2次元形状の変化率を測定した時点で水分率計を用いて測定した複数の酒米の吸水率の平均値を算出し、酒米の3次元形状の変化率と吸水率との関係をプロットしておく。そして、プロットした複数の点にとっての近似線が、状態算出部30が参照する、酒米の3次元形状の変化率と吸水率との間の相関関係として、記憶装置35に記憶しておくことで、推定した3次元形状の変化率と、上記相関関係とに基づいて、酒米の吸水率を算出することもできる。
尚、上記第1及び第2実施形態において、水に浸漬する前の酒米の厚み方向と直交する平面に沿った面積を予め記憶しておけば、状態算出部30が、予め記憶された面積と面積情報取得部13で取得した面積とを基に算出される2次元形状の変化率に、推定部19で推定された厚み方向における形状の変化率を乗じて3次元形状の変化率を算出するように構成することで、上記と同様に、酒米の吸水率を算出することもできる。
[3] In the third embodiment, the state calculator 30 calculates the correlation between the estimated change rate of the three-dimensional shape and the previously stored three-dimensional shape of the sake rice and the water absorption of the sake rice. , The moisture content of sake rice can be calculated. That is, for a plurality of sake rice arranged in a state of being immersed in water, the average value of the change rate of the three-dimensional shape obtained in advance by a known method was calculated and used to obtain the change rate of the three-dimensional shape. Calculate the average value of the water absorption rate of a plurality of sake rice measured using a moisture meter at the time when the change rate of the three-dimensional shape is measured, and plot the relationship between the change rate of the three-dimensional shape of the sake rice and the water absorption rate. Keep it. Then, the approximate lines for the plotted points are stored in the storage device 35 as the correlation between the change rate of the three-dimensional shape of the sake rice and the water absorption referred to by the state calculation unit 30. Then, the water absorption of the sake rice can also be calculated based on the estimated change rate of the three-dimensional shape and the correlation.
In the first and second embodiments, if the area along a plane orthogonal to the thickness direction of the sake rice before immersion in water is stored in advance, the state calculation unit 30 can store the stored area in advance. The rate of change of the three-dimensional shape is calculated by multiplying the rate of change of the two-dimensional shape calculated based on the area obtained by the area information obtaining unit 13 and the rate of change of the shape in the thickness direction estimated by the estimating unit 19. With such a configuration, the water absorption of the sake rice can be calculated in the same manner as described above.

[4]上記実施形態では、状態推定装置1A,1B,1C及び変化率推定装置20A,20B,20Cの構成について具体例を挙げて説明したが、その構成は適宜変更可能である。 [4] In the above embodiment, the configurations of the state estimation devices 1A, 1B, and 1C and the change rate estimation devices 20A, 20B, and 20C have been described with specific examples. However, the configurations can be appropriately changed.

[5]上記第1及び第2実施形態では、形態判定部17が、CNN(Convolutional Neural Network)を用いて酒米の形態を予測判定する例を説明したが、形態判定部17を構成するアルゴリズムは上述したCNNに限定されず、サポートベクトルマシン(SVM)、k近傍法(k-nearest neighbor algorithm)、判別関数などの他の様々なアルゴリズムを用いて構成することができる。
また、形態判定部17は、酒米に割れが発生しているか否かを、上記実施形態とは異なる様々な画像認識技術を用いて判定するように構成してもよい。或いは、形態判定部17は、人間により行われた、酒米に割れが発生しているか否かの判定結果を受け付けて、その結果を自身の判定結果とするように構成してもよい。
[5] In the first and second embodiments, an example has been described in which the morphological determining unit 17 predicts and determines the morphology of sake rice using a CNN (Convolutional Neural Network). Is not limited to the above-described CNN, and can be configured using various other algorithms such as a support vector machine (SVM), a k-nearest neighbor algorithm, and a discriminant function.
Further, the form determining unit 17 may be configured to determine whether or not cracked rice has been generated using various image recognition techniques different from the above-described embodiment. Alternatively, the form determining unit 17 may be configured to receive a result of determination made by a human as to whether or not cracked rice is generated, and use the result as its own determination result.

[6]上記実施形態(別実施形態を含む)で開示される構成は、矛盾が生じない限り、他の実施形態で開示される構成と組み合わせて適用することが可能であり、また、本明細書において開示された実施形態は例示であって、本発明の実施形態はこれに限定されず、本発明の目的を逸脱しない範囲内で適宜改変することが可能である。 [6] The configuration disclosed in the above embodiment (including another embodiment) can be applied in combination with the configuration disclosed in another embodiment unless there is a contradiction. The embodiment disclosed in the present specification is an exemplification, and the embodiment of the present invention is not limited thereto, and can be appropriately modified without departing from the object of the present invention.

本発明は、米の厚み方向における形状の変化率、又は、米の3次元形状の変化率を推定できる変化率推定装置に利用でき、また、推定した厚み方向における形状の変化率や3次元形状の変化率を基にして米の比重、密度及び吸水率を推定できる状態推定装置に利用できる。   INDUSTRIAL APPLICABILITY The present invention can be used for a change rate estimating device capable of estimating a change rate of a shape in a thickness direction of rice or a change rate of a three-dimensional shape of rice. It can be used for a state estimating device that can estimate the specific gravity, density and water absorption of rice based on the rate of change of the rice.

1(1A,1B,1C) 状態推定装置
2 撮影装置
3 出力装置
10(10A,10B,10C) 米分析装置
11 画像データ取得部
13 面積情報取得部
16 割れ率算出部
19 推定部
20(20A,20B,20C) 変化率推定装置
30 状態算出部
40 2次元形状変化率算出部
1 (1A, 1B, 1C) State estimation device 2 Imaging device 3 Output device 10 (10A, 10B, 10C) Rice analyzer 11 Image data acquisition unit 13 Area information acquisition unit 16 Crack rate calculation unit 19 Estimation unit 20 (20A, 20B, 20C) Change rate estimating device 30 State calculating unit 40 Two-dimensional shape change rate calculating unit

Claims (14)

水に浸漬した状態の対象米を撮影する撮影手段と、
前記撮影手段で撮影した前記対象米の、当該対象米を厚み方向から視た2次元画像データを取得する画像データ取得手段と、
前記画像データ取得手段で取得された2次元画像データを参照して、当該2次元画像データを解析し、複数の対象米のうち、前記対象米の長軸方向に沿う縦断面が開いた縦割れが発生している対象米の割合である縦割れ率を算出する割れ率算出手段と、
前記縦割れ率を前記割れ率算出手段から取得し、前記縦割れ率、及び予め定められた前記対象米の縦割れ率と前記厚み方向における形状の変化率との相関関係を基に、前記対象米の、前記厚み方向における形状の変化率を推定する推定手段とを備える変化率推定装置。
Photographing means for photographing the target rice immersed in water,
Image data acquisition means for acquiring two-dimensional image data of the target rice photographed by the photographing means, as viewed from the thickness direction of the target rice;
The two-dimensional image data is analyzed with reference to the two-dimensional image data acquired by the image data acquiring means, and a longitudinal crack having an open longitudinal section along a long axis direction of the target rice among a plurality of target rice is opened. Crack rate calculating means for calculating a vertical crack rate which is a ratio of target rice in which cracks have occurred,
The vertical cracking rate is obtained from the cracking rate calculating means, and the vertical cracking rate, and the predetermined correlation between the vertical cracking rate of the target rice and the change rate of the shape in the thickness direction are determined based on the correlation. A change rate estimating device comprising: estimating means for estimating a change rate of the shape of rice in the thickness direction.
請求項1に記載の変化率推定装置と、
前記厚み方向における形状の変化率を前記推定手段から取得し、前記取得した厚み方向における形状の変化率を基に、前記対象米の比重、密度又は吸水率を算出する状態算出手段とを備える状態推定装置。
A change rate estimating device according to claim 1,
A state calculating unit that obtains the shape change rate in the thickness direction from the estimating unit, and calculates a specific gravity, density, or water absorption of the target rice based on the obtained shape change rate in the thickness direction. Estimation device.
前記画像データ取得手段で取得された2次元画像データを参照して、当該2次元画像データを解析し、前記厚み方向と直交する平面に沿った方向における前記対象米の面積を取得する面積情報取得手段を備え、
前記状態算出手段において、前記面積を前記面積情報取得手段から取得し、前記取得した厚み方向における形状の変化率、前記取得した面積、及び予め記憶された水に浸漬する前の前記対象米の厚みを基に算出される前記対象米の体積を基に、前記対象米の比重又は密度を算出する請求項2に記載の状態推定装置。
Area information acquisition for analyzing the two-dimensional image data with reference to the two-dimensional image data acquired by the image data acquisition means and acquiring the area of the target rice in a direction along a plane orthogonal to the thickness direction With means,
In the state calculating means, the area is obtained from the area information obtaining means, and the obtained shape change rate in the thickness direction, the obtained area, and the thickness of the target rice before immersion in pre-stored water. The state estimating device according to claim 2, wherein the specific gravity or density of the target rice is calculated based on the volume of the target rice calculated based on the following.
水に浸漬した状態の対象米を撮影する撮影手段と、
前記撮影手段で撮影した前記対象米の、当該対象米を厚み方向から視た2次元画像データを取得する画像データ取得手段と、
前記画像データ取得手段で取得された2次元画像データを参照して、当該2次元画像データを解析し、前記厚み方向と直交する平面に沿った方向における前記対象米の2次元形状の変化率を算出する2次元形状変化率算出手段と、
前記2次元形状変化率算出手段から前記2次元形状の変化率を取得し、複数の対象米のうち、前記対象米の長軸方向に沿う縦断面が開いた縦割れが発生している対象米の割合である縦割れ率と前記厚み方向における形状の変化率との相関関係、及び前記取得した2次元形状の変化率を基に、前記2次元形状の変化率を補正して3次元形状の変化率として推定する推定手段とを備える変化率推定装置。
Photographing means for photographing the target rice immersed in water,
Image data acquisition means for acquiring two-dimensional image data of the target rice photographed by the photographing means, as viewed from the thickness direction of the target rice;
The two-dimensional image data is analyzed with reference to the two-dimensional image data acquired by the image data acquiring means, and the rate of change of the two-dimensional shape of the target rice in a direction along a plane orthogonal to the thickness direction is calculated. A two-dimensional shape change rate calculating means for calculating;
The rate of change of the two-dimensional shape is obtained from the two-dimensional shape change rate calculating means, and a target rice having a vertical crack having an open longitudinal section along a long axis direction of the target rice among a plurality of target rices The rate of change of the two-dimensional shape is corrected based on the correlation between the vertical cracking rate, which is the ratio of the shape, and the rate of change of the shape in the thickness direction, and the obtained rate of change of the two-dimensional shape. A change rate estimating device comprising: estimating means for estimating a change rate.
前記推定手段において、前記相関関係に基づいて定まる縦割れ影響度を、前記取得した2次元形状の変化率に乗じて、前記2次元形状の変化率を前記3次元形状の変化率に補正する請求項4に記載の変化率推定装置。   The estimation means corrects the change rate of the two-dimensional shape to the change rate of the three-dimensional shape by multiplying the degree of influence of a vertical crack determined based on the correlation by the obtained change rate of the two-dimensional shape. Item 5. The change rate estimation device according to item 4. 請求項4又は5に記載の変化率推定装置と、
前記3次元形状の変化率を前記推定手段から取得し、前記取得した3次元形状の変化率を基に、前記対象米の比重、密度又は吸水率を算出する状態算出手段とを備える状態推定装置。
A change rate estimating device according to claim 4 or 5,
A state estimating device comprising: a state calculator that obtains the change rate of the three-dimensional shape from the estimator; and calculates a specific gravity, density, or water absorption of the target rice based on the obtained change rate of the three-dimensional shape. .
前記状態算出手段において、前記取得した3次元形状の変化率、並びに予め記憶された水に浸漬する前の前記対象米の厚み及び厚み方向と直交する平面に沿った方向における面積を基に算出される前記対象米の体積を基に、前記対象米の比重又は密度を算出する請求項6に記載の状態推定装置。   The state calculation means calculates the change rate of the acquired three-dimensional shape, and the thickness of the target rice before immersion in water stored in advance and the area in a direction along a plane orthogonal to the thickness direction of the target rice. The state estimating device according to claim 6, wherein a specific gravity or a density of the target rice is calculated based on a volume of the target rice. 水に浸漬した状態の対象米を撮影する撮影工程と、
前記撮影した前記対象米の、当該対象米を厚み方向から視た2次元画像を基に、複数の対象米のうち、前記対象米の長軸方向に沿う縦断面が開いた縦割れが発生している対象米の割合である縦割れ率を算出する割れ率算出工程と、
前記算出した縦割れ率、及び予め定められた前記対象米の縦割れ率と前記厚み方向における形状の変化率との相関関係を基に、前記対象米の、前記厚み方向における形状の変化率を推定する推定工程とを実行する変化率推定方法。
A photographing process of photographing the target rice immersed in water,
On the basis of a two-dimensional image of the target rice that has been photographed and the target rice viewed from the thickness direction, a vertical crack having an open longitudinal section along the long axis direction of the target rice among a plurality of target rice occurs. Cracking ratio calculating step of calculating a vertical cracking ratio which is a ratio of target rice having
Based on the calculated vertical cracking rate, and a predetermined correlation between the vertical cracking rate of the target rice and the change rate of the shape in the thickness direction, the target rice, the change rate of the shape in the thickness direction is calculated. A rate-of-change estimating method for performing an estimating step of estimating.
請求項8に記載の変化率推定方法を行った後に、前記変化率推定方法で推定した前記厚み方向における形状の変化率を基に、前記対象米の比重、密度又は吸水率を算出する状態算出工程を実行する状態推定方法。   A state calculation for calculating a specific gravity, density, or water absorption of the target rice based on the change rate of the shape in the thickness direction estimated by the change rate estimation method after performing the change rate estimation method according to claim 8. State estimation method for executing the process. 前記撮影した前記対象米の、当該対象米を厚み方向から視た2次元画像を基に、前記厚み方向と直交する平面に沿った方向における前記対象米の面積を取得する面積情報取得工程を備え、
前記状態算出工程において、前記厚み方向における形状の変化率、前記取得した面積、及び水に浸漬する前の前記対象米の厚みを基に算出される前記対象米の体積を基に、前記対象米の比重又は密度を算出する請求項9に記載の状態推定方法。
An area information acquisition step of acquiring an area of the target rice in a direction along a plane orthogonal to the thickness direction based on a two-dimensional image of the target rice that has been photographed and viewing the target rice from the thickness direction. ,
In the state calculation step, based on the shape change rate in the thickness direction, the acquired area, and the volume of the target rice calculated based on the thickness of the target rice before immersion in water, the target rice The state estimating method according to claim 9, wherein the specific gravity or the density is calculated.
水に浸漬した状態の対象米を撮影する撮影工程と、
前記撮影した前記対象米の、当該対象米を厚み方向から視た2次元画像を基に、前記厚み方向と直交する平面に沿った方向における前記対象米の2次元形状の変化率を算出する2次元形状変化率算出工程と、
複数の対象米のうち、前記対象米の長軸方向に沿う縦断面が開いた縦割れが発生している対象米の割合である縦割れ率と前記厚み方向における形状の変化率との相関関係、及び前記算出した2次元形状の変化率を基に、前記2次元形状の変化率を補正して3次元形状の変化率として推定する推定工程とを実行する変化率推定方法。
A photographing process of photographing the target rice immersed in water,
Calculating a change rate of a two-dimensional shape of the target rice in a direction along a plane orthogonal to the thickness direction based on a two-dimensional image of the target rice taken in the thickness direction of the target rice; Dimensional shape change rate calculating step,
Among a plurality of target rice, a correlation between a vertical cracking rate which is a ratio of target rice in which a vertical section along a longitudinal direction along the long axis direction of the target rice is open and a rate of shape change in the thickness direction. And an estimation step of correcting the change rate of the two-dimensional shape based on the calculated change rate of the two-dimensional shape to estimate the change rate of the two-dimensional shape as the change rate of the three-dimensional shape.
前記推定工程において、前記相関関係に基づいて定まる縦割れ影響度を、前記算出した2次元形状の変化率に乗じて、前記2次元形状の変化率を前記3次元形状の変化率に補正する請求項11に記載の変化率推定方法。   In the estimating step, a vertical crack influence degree determined based on the correlation is multiplied by the calculated two-dimensional shape change rate, and the two-dimensional shape change rate is corrected to the three-dimensional shape change rate. Item 12. The change rate estimation method according to Item 11. 請求項11又は12に記載の変化率推定方法を行った後に、前記変化率推定方法で推定した前記3次元形状の変化率を基に、前記対象米の比重、密度又は吸水率を算出する状態算出工程を実行する状態推定方法。   A state in which after performing the change rate estimation method according to claim 11 or 12, the specific gravity, density, or water absorption of the target rice is calculated based on the change rate of the three-dimensional shape estimated by the change rate estimation method. A state estimating method for executing the calculating step. 前記状態算出工程において、前記3次元形状の変化率、並びに予め記憶された水に浸漬する前の前記対象米の厚み及び厚み方向と直交する平面に沿った方向における面積を基に算出される前記対象米の体積を基に、前記対象米の比重又は密度を算出する請求項13に記載の状態推定方法。
In the state calculation step, the change rate of the three-dimensional shape, and the thickness of the target rice before immersion in water stored in advance and the area in a direction along a plane orthogonal to the thickness direction of the target rice is calculated based on the area 14. The state estimation method according to claim 13, wherein the specific gravity or density of the target rice is calculated based on the volume of the target rice.
JP2018167813A 2018-09-07 2018-09-07 Change rate estimation device, change rate estimation method, state estimation device and state estimation method Active JP7085952B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2018167813A JP7085952B2 (en) 2018-09-07 2018-09-07 Change rate estimation device, change rate estimation method, state estimation device and state estimation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2018167813A JP7085952B2 (en) 2018-09-07 2018-09-07 Change rate estimation device, change rate estimation method, state estimation device and state estimation method

Publications (2)

Publication Number Publication Date
JP2020041841A true JP2020041841A (en) 2020-03-19
JP7085952B2 JP7085952B2 (en) 2022-06-17

Family

ID=69798039

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2018167813A Active JP7085952B2 (en) 2018-09-07 2018-09-07 Change rate estimation device, change rate estimation method, state estimation device and state estimation method

Country Status (1)

Country Link
JP (1) JP7085952B2 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021166357A1 (en) * 2020-02-18 2021-08-26 パナソニックIpマネジメント株式会社 Measurement method and measurement device

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005035786A (en) * 2003-06-22 2005-02-10 Pusai:Kk Posture control method and posture control device for shape-anisotropic particle
JP2006071552A (en) * 2004-09-03 2006-03-16 Yamamoto Co Ltd Apparatus for determining state of particulate object to be inspected
JP2014163865A (en) * 2013-02-27 2014-09-08 Iwate Industrial Research Center Sake rice inspection apparatus
JP2016105074A (en) * 2014-11-21 2016-06-09 大阪瓦斯株式会社 Moisture state determination method of rice to be cooked, immersion time determination method and rice cooking facility
JP2019197000A (en) * 2018-05-10 2019-11-14 大阪瓦斯株式会社 Sake rice analysis device
JP2019200158A (en) * 2018-05-17 2019-11-21 大阪瓦斯株式会社 Water content estimation device

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005035786A (en) * 2003-06-22 2005-02-10 Pusai:Kk Posture control method and posture control device for shape-anisotropic particle
JP2006071552A (en) * 2004-09-03 2006-03-16 Yamamoto Co Ltd Apparatus for determining state of particulate object to be inspected
JP2014163865A (en) * 2013-02-27 2014-09-08 Iwate Industrial Research Center Sake rice inspection apparatus
JP2016105074A (en) * 2014-11-21 2016-06-09 大阪瓦斯株式会社 Moisture state determination method of rice to be cooked, immersion time determination method and rice cooking facility
JP2019197000A (en) * 2018-05-10 2019-11-14 大阪瓦斯株式会社 Sake rice analysis device
JP2019200158A (en) * 2018-05-17 2019-11-21 大阪瓦斯株式会社 Water content estimation device

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"高品質の大吟醸酒造りを画像認識技術により支援 伏見酒造組合と大阪ガスが酒米評価の新手法を開発", 大阪ガス プレスリリース(WEB), JPN6022016782, 20 July 2018 (2018-07-20), ISSN: 0004766050 *
中山繁喜、高橋亨: "酒造用白米の浸漬割れ測定法", 岩手県工業技術センター研究報告, vol. 13号, JPN6022016781, May 2006 (2006-05-01), pages 45 - 48, ISSN: 0004766051 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021166357A1 (en) * 2020-02-18 2021-08-26 パナソニックIpマネジメント株式会社 Measurement method and measurement device

Also Published As

Publication number Publication date
JP7085952B2 (en) 2022-06-17

Similar Documents

Publication Publication Date Title
JP7257745B2 (en) Moisture content estimator
US20180303409A1 (en) Osteoporosis diagnostic support apparatus
JP6608138B2 (en) Image processing apparatus, image processing method, image processing program, and recording medium storing the program
JP6688962B2 (en) Judgment device, judgment method, and judgment program
WO2022105676A1 (en) Method and system for measuring wear of workpiece plane
CN112802071B (en) Three-dimensional reconstruction effect evaluation method and system
JP2015219603A5 (en)
CN112734858B (en) Binocular calibration precision online detection method and device
JP2017227606A (en) Defect detection device and defect detection method
CN108369737A (en) Using heuristic graph searching quickly and automatically to divide layered image
JP2020186957A (en) Vibration analysis system, vibration analysis method, and program
JP6278880B2 (en) Water level measuring device
JP6700790B2 (en) Method for determining water content of rice to be cooked, water content determination device, immersion time determination device, and rice cooking equipment
JP7085952B2 (en) Change rate estimation device, change rate estimation method, state estimation device and state estimation method
JP7023180B2 (en) Sake rice analyzer
CN111311598A (en) Multi-view-angle-based fastener surface defect detection method and device
JPWO2019172206A5 (en)
JP7263149B2 (en) Image processing device, image processing method, and program
CN117058106A (en) Method for measuring flatness and surface defects of flexible glass based on random forest
JP7275595B2 (en) Ozone degradation evaluation system for vulcanized rubber materials
JP2004109105A (en) Flaw type classification boundary setting method in surface defect detection, and defect detection method
JP6904062B2 (en) Information processing equipment, information processing methods and programs
EP2693397B1 (en) Method and apparatus for noise reduction in an imaging system
JP5102139B2 (en) Thin part detection method and thin part detection device
CN112700446A (en) Algorithm model training method and device for industrial quality inspection

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20210524

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20220408

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: 20220510

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20220607

R150 Certificate of patent or registration of utility model

Ref document number: 7085952

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150