Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP7688907B2 - Sugar content measurement method, sugar content measurement device, and sugar content measurement program - Google Patents
[go: Go Back, main page]

JP7688907B2 - Sugar content measurement method, sugar content measurement device, and sugar content measurement program - Google Patents

Sugar content measurement method, sugar content measurement device, and sugar content measurement program Download PDF

Info

Publication number
JP7688907B2
JP7688907B2 JP2021204633A JP2021204633A JP7688907B2 JP 7688907 B2 JP7688907 B2 JP 7688907B2 JP 2021204633 A JP2021204633 A JP 2021204633A JP 2021204633 A JP2021204633 A JP 2021204633A JP 7688907 B2 JP7688907 B2 JP 7688907B2
Authority
JP
Japan
Prior art keywords
sugar content
light
light source
vegetables
measurement
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2021204633A
Other languages
Japanese (ja)
Other versions
JP2022095598A (en
Inventor
秀和 伊藤
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.)
National Agriculture and Food Research Organization
Original Assignee
National Agriculture and Food Research Organization
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 National Agriculture and Food Research Organization filed Critical National Agriculture and Food Research Organization
Publication of JP2022095598A publication Critical patent/JP2022095598A/en
Application granted granted Critical
Publication of JP7688907B2 publication Critical patent/JP7688907B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigating Or Analysing Materials By Optical Means (AREA)

Description

本発明は、糖度計測方法、糖度計測装置、及び糖度計測プログラムに関する。 The present invention relates to a sugar content measurement method, a sugar content measurement device, and a sugar content measurement program.

例えば、特許文献1の非破壊糖度測定装置では、860nmから960nmの範囲にある3種類の波長の光を出射する単一若しくは複数の光源を使用すると共に、上記光吸収を検出する検出器の配置位置が、光源から出射され上記青果物へ入射する入射光の青果物表面における照射領域の中心点と青果物の中心とを結ぶ直線の延長線上の位置を除き、かつ入射光の青果物表面における照射領域と検出器が受光する上記青果物からの出射光の青果物表面における検出領域とが重ならない位置に設定されている。 For example, the non-destructive sugar content measuring device in Patent Document 1 uses a single or multiple light sources that emit light of three wavelengths in the range of 860 nm to 960 nm, and the detector that detects the light absorption is positioned so that it is located at a position other than the position on the extension line of the straight line connecting the center point of the irradiated area on the surface of the fruit or vegetable of the incident light emitted from the light source and incident on the fruit or vegetable and the center of the fruit or vegetable, and so that the irradiated area on the surface of the incident light does not overlap with the detection area on the surface of the fruit or vegetable of the emitted light received by the detector from the fruit or vegetable.

しかしながら、特許文献1では、860~890nm、900~920nm、920を越えて~960nmの3種類の光源と説明変数を使用する。しかし、各波長域は広いため各波長域で有効な波長をさらに検討する必要があり、かつ、3波長域を適用した場合に糖度の非破壊計測精度が低いという問題がある。 However, in Patent Document 1, three types of light sources and explanatory variables are used: 860-890 nm, 900-920 nm, and over 920 to 960 nm. However, because each wavelength range is broad, it is necessary to further consider the effective wavelengths in each wavelength range, and there is also the problem that the accuracy of non-destructive measurement of sugar content is low when three wavelength ranges are applied.

特許第3346142号公報Patent No. 3346142 特開2000-221134号公報JP 2000-221134 A

岩元、日本食品工業学会誌、27、pp464-471(1980)Iwamoto, Journal of the Japanese Society for Food Technology, 27, pp. 464-471 (1980) S.Morimoto、et al.、Near Infrared Spectroscopy:Proceedings of the 10th International Conference pp155-159(2002)S. Morimoto et al. , Near Infrared Spectroscopy: Proceedings of the 10th International Conference pp155-159 (2002) 伊藤、森本、照明学会誌、98、pp581-584(2014)Ito, Morimoto, Journal of the Illuminating Engineering Society of Japan, 98, pp581-584 (2014)

本発明は、上記に鑑みてなされたもので、青果物の糖度を計測する場合に、高精度に計測することが可能な糖度計測方法、糖度計測装置、及び糖度計測プログラムを提供することを目的とする。 The present invention has been made in consideration of the above, and aims to provide a sugar content measurement method, sugar content measurement device, and sugar content measurement program that are capable of measuring the sugar content of fruits and vegetables with high accuracy.

上述した課題を解決し、目的を達成するために、本発明は、光源からの光源光を青果物に照射し、その拡散反射を含む反射光を受光して当該青果物の糖度を計測する糖度計測方法であって、1つの光源又は複数の光源から照射される近赤外線波長域の光を青果物に照射し、その拡散反射を含む反射光から波長の異なる4つの、又はこれら4つの波長を含む波長域の分光吸光スペクトルの光のシグナル値を取得する取得工程と、取得した波長の異なる4つの吸光度(光のシグナル値)を説明変数として使用して、多変量解析を行うことで、糖度を推定するための推定モデルを作成する推定モデル作成工程と、を含むことを特徴とする。 In order to solve the above problems and achieve the object, the present invention provides a sugar content measurement method that irradiates fruit or vegetables with light from a light source and receives the reflected light, including diffuse reflection, to measure the sugar content of the fruit or vegetables, and is characterized by including an acquisition step of irradiating the fruit or vegetables with light in the near-infrared wavelength range from one light source or multiple light sources and acquiring light signal values of a spectral absorbance spectrum of four different wavelengths or a wavelength range including these four wavelengths from the reflected light, including diffuse reflection, and an estimation model creation step of performing multivariate analysis using the acquired absorbances (light signal values) of the four different wavelengths as explanatory variables to create an estimation model for estimating sugar content.

また、本発明の一態様によれば、前記異なる4つの波長は、856±2nm、876±2nm又は884±2nm、900~918nm、926±2nmであることにしてもよい。 Furthermore, according to one aspect of the present invention, the four different wavelengths may be 856±2 nm, 876±2 nm or 884±2 nm, 900 to 918 nm, and 926±2 nm.

また、本発明の一態様によれば、前記推定モデル作成工程では、前記多変量解析として重回帰分析を行って、前記推定モデルとして重回帰式を作成することにしてもよい。 According to one aspect of the present invention, in the estimation model creation step, a multiple regression analysis may be performed as the multivariate analysis to create a multiple regression equation as the estimation model.

また、本発明の一態様によれば、前記青果物は、当該青果物の赤道部、花痕部、加工品、及び破砕物を含むことにしてもよい。 According to one aspect of the present invention, the fruits and vegetables may include the equator portion, the flower scar portion, processed products, and crushed products of the fruits and vegetables.

また、本発明の一態様によれば、前記青果物は、西洋ナシ、カキ、マンゴー、イチゴ、パプリカ、ミカン(Tangerineを含む)、不知火、トマト、サクランボ、モモ、ナシ、リンゴ、スモモ、及びメロンを含むことにしてもよい。 According to one aspect of the present invention, the fruits and vegetables may include pears, persimmons, mangoes, strawberries, paprika, mandarins (including tangerines), shiranui, tomatoes, cherries, peaches, pears, apples, plums, and melons.

また、本発明の一態様によれば、前記光源は、リング状の光源であることにしてもよい。 According to one aspect of the present invention, the light source may be a ring-shaped light source.

また、本発明の一態様によれば、前記取得工程では、前記リング状の光源から照射される光に対する青果物の拡散反射を含む反射光を、当該リング状の光源の略中心で検出することにしてもよい。 According to one aspect of the present invention, in the acquisition process, reflected light, including diffuse reflection from fruits and vegetables, of the light irradiated from the ring-shaped light source may be detected at approximately the center of the ring-shaped light source.

また、本発明の一態様によれば、前記1つの光源は、ハロゲンランプであることにしてもよい。 According to one aspect of the present invention, the one light source may be a halogen lamp.

また、本発明の一態様によれば、前記複数の光源は、発光波長の異なる複数のLEDであることにしてもよい。 According to one aspect of the present invention, the multiple light sources may be multiple LEDs with different emission wavelengths.

また、本発明の一態様によれば、前記複数のLEDは、計測対象の青果物側に平面に対して所定角度(但し、所定角度は10度より大きい)を有して配置されることにしてもよい。 According to one aspect of the present invention, the LEDs may be arranged at a predetermined angle (however, the predetermined angle is greater than 10 degrees) with respect to a plane on the side of the fruit or vegetable being measured.

また、上述した課題を解決し、目的を達成するために、本発明は、光源からの光源光を青果物に照射し、その拡散反射を含む反射光を受光して当該青果物の糖度を計測する糖度計測装置であって、1つの光源又は複数の光源から照射される近赤外線波長域の光を青果物に照射し、その拡散反射を含む反射光から波長の異なる4つの、又はこれら4つの波長を含む波長域の分光吸光スペクトルの光のシグナル値を取得する分光検出手段と、取得した波長の異なる4つの吸光度(光のシグナル値)を説明変数として使用して、多変量解析を行うことで、糖度を推定するための推定モデルを作成する推定モデル作成手段と、を備えたことを特徴とする。 In order to solve the above-mentioned problems and achieve the object, the present invention provides a sugar content measuring device that irradiates fruit or vegetables with light from a light source and receives the reflected light, including diffuse reflection, to measure the sugar content of the fruit or vegetables, and is characterized by comprising: a spectroscopic detection means that irradiates the fruit or vegetables with light in the near-infrared wavelength range from one light source or multiple light sources and acquires light signal values of a spectral absorbance spectrum of four different wavelengths or a wavelength range including these four wavelengths from the reflected light, including diffuse reflection; and an estimation model creation means that creates an estimation model for estimating sugar content by performing multivariate analysis using the acquired absorbances (light signal values) of the four different wavelengths as explanatory variables.

また、本発明の一態様によれば、前記異なる4つの波長は、856±2nm、876±2nm又は884±2nm、900~918nm、926±2nmであることにしてもよい。 Furthermore, according to one aspect of the present invention, the four different wavelengths may be 856±2 nm, 876±2 nm or 884±2 nm, 900 to 918 nm, and 926±2 nm.

また、本発明の一態様によれば、前記推定モデル作成手段は、前記多変量解析として重回帰分析を行って、前記推定モデルとして重回帰式を作成することにしてもよい。 According to one aspect of the present invention, the estimation model creation means may perform a multiple regression analysis as the multivariate analysis and create a multiple regression equation as the estimation model.

また、本発明の一態様によれば、前記青果物は、当該青果物の赤道部、花痕部、加工品、及び破砕物を含むことにしてもよい。 According to one aspect of the present invention, the fruits and vegetables may include the equator portion, the flower scar portion, the processed product, and the crushed product of the fruits and vegetables.

また、本発明の一態様によれば、前記青果物は、西洋ナシ、カキ、マンゴー、イチゴ、パプリカ、ミカン(Tangerineを含む)、不知火、トマト、サクランボ、モモ、ナシ、リンゴ、スモモ、及びメロンを含むことにしてもよい。 According to one aspect of the present invention, the fruits and vegetables may include pears, persimmons, mangoes, strawberries, paprika, mandarins (including tangerines), shiranui, tomatoes, cherries, peaches, pears, apples, plums, and melons.

また、本発明の一態様によれば、前記光源は、リング状の光源であることにしてもよい。 According to one aspect of the present invention, the light source may be a ring-shaped light source.

また、本発明の一態様によれば、前記分光検出手段は、前記リング状の光源から照射される光に対する青果物の拡散反射を含む反射光を、当該リング状の光源の略中心で検出することにしてもよい。 According to one aspect of the present invention, the spectroscopic detection means may detect reflected light, including diffuse reflection from fruits and vegetables, of the light irradiated from the ring-shaped light source at approximately the center of the ring-shaped light source.

また、本発明の一態様によれば、前記1つの光源は、ハロゲンランプであることにしてもよい。 According to one aspect of the present invention, the one light source may be a halogen lamp.

また、本発明の一態様によれば、前記複数の光源は、発光波長の異なる複数のLEDであることにしてもよい。 According to one aspect of the present invention, the multiple light sources may be multiple LEDs with different emission wavelengths.

また、本発明の一態様によれば、前記複数のLEDは、計測対象の青果物側に平面に対して所定角度(但し、所定角度は10度より大きい)を有して配置されることにしてもよい。 According to one aspect of the present invention, the LEDs may be arranged at a predetermined angle (however, the predetermined angle is greater than 10 degrees) with respect to a plane on the side of the fruit or vegetable being measured.

また、上述した課題を解決し、目的を達成するために、本発明は、青果物の糖度を計測するための糖度計測プログラムであって、1つの光源又は複数の光源から照射される近赤外線波長域の光に対する、青果物からの拡散反射を含む反射光から波長の異なる4つの、又はこれら4つの波長を含む波長域の分光吸光スペクトルの光のシグナル値を取得する取得工程と、取得した波長の異なる4つの吸光度(光のシグナル値)を説明変数として使用して、多変量解析を行うことで、糖度を推定するための推定モデルを作成する推定モデル作成工程と、コンピュータに実行させるための糖度計測プログラムであることを特徴とする。 In order to solve the above-mentioned problems and achieve the object, the present invention is a sugar content measurement program for measuring the sugar content of fruits and vegetables, which includes an acquisition step of acquiring light signal values of a spectral absorbance spectrum of four different wavelengths or a wavelength range including these four wavelengths from reflected light, including diffuse reflection, from fruits and vegetables in response to light in the near-infrared wavelength range irradiated from one light source or multiple light sources, and an estimation model creation step of creating an estimation model for estimating sugar content by performing multivariate analysis using the acquired four absorbances (light signal values) of different wavelengths as explanatory variables, and a sugar content measurement program for causing a computer to execute the process.

この発明によれば、青果物の糖度を計測する場合に、高精度に計測することが可能になるという効果を奏する。 This invention has the effect of making it possible to measure the sugar content of fruits and vegetables with high accuracy.

図1は、本発明の概略を説明するための説明図である。FIG. 1 is an explanatory diagram for explaining an outline of the present invention. 図2は、本実施の形態に係る糖度計測方法を説明するためのフローチャートである。FIG. 2 is a flowchart for explaining the sugar content measuring method according to the present embodiment. 図3は、本実施の形態に係る糖度計測装置の外観構成例を示す図である。FIG. 3 is a diagram showing an example of the external configuration of a sugar content measuring device according to this embodiment. 図4は、本実施の形態に係る糖度計測装置の構成の一例を示すブロック図である。FIG. 4 is a block diagram showing an example of the configuration of a sugar content measuring device according to this embodiment. 図5は、高精度実用機の要部を説明するための図である。FIG. 5 is a diagram for explaining the main parts of a high-precision practical machine. 図6-Aは、LED試作機の要部を説明するための図であり、LED試作機を上から見た概略の平面図である。FIG. 6-A is a diagram for explaining the main parts of the LED prototype, and is a schematic plan view of the LED prototype as seen from above. 図6-Bは、LED試作機の要部を説明するための図であり、LED試作機の概略の断面図である。FIG. 6-B is a diagram for explaining the main parts of the LED prototype, and is a schematic cross-sectional view of the LED prototype. 図6-Cは、LED試作機の要部を説明するための図であり、LED試作機の試料台に試料をセットして、LEDから照射された光の経路を説明するための模式図である。FIG. 6-C is a diagram for explaining the main parts of the LED prototype, and is a schematic diagram for explaining the path of light irradiated from the LED when a sample is set on the sample stage of the LED prototype. 図6-Dは、LED試作機の要部を説明するための図であり、LED試作機の概略の斜視図である。FIG. 6-D is a diagram for explaining the main parts of the LED prototype, and is a schematic perspective view of the LED prototype. 図7は、モモの高精度実用機での計測結果を示す図である。FIG. 7 is a diagram showing the measurement results of peaches using a high-precision practical machine. 図8は、モモのLED試作機での計測結果を示す図である。FIG. 8 is a diagram showing the measurement results of the peach LED prototype. 図9は、ナシの高精度実用機での計測結果を示す図である。FIG. 9 is a diagram showing the results of measurement of pears using a high-precision practical machine. 図10は、ナシのLED試作機での計測結果を示す図である。FIG. 10 is a diagram showing the measurement results of the pear LED prototype. 図11は、リンゴの高精度実用機での計測結果を示す図である。FIG. 11 is a diagram showing the measurement results of an apple using a high-precision practical machine. 図12は、リンゴのLED試作機での計測結果を示す図である。FIG. 12 is a diagram showing the measurement results of the LED prototype for apples. 図13は、西洋ナシの高精度実用機での計測結果を示す図である。FIG. 13 is a diagram showing the measurement results of a pear using a high-precision practical machine. 図14は、カキの高精度実用機での計測結果を示す図である。FIG. 14 shows the results of measurements of oysters taken using a high-precision practical machine. 図15は、スモモの高精度実用機での計測結果を示す図である。FIG. 15 shows the measurement results of plums using a high-precision practical machine. 図16は、スモモのLED試作機での計測結果を示す図である。FIG. 16 is a diagram showing the measurement results of plums using the LED prototype. 図17は、マンゴーの高精度実用機での計測結果を示す図である。FIG. 17 shows the measurement results of mango using a high-precision practical machine. 図18は、トマト赤道部の高精度実用機での計測結果を示す図である。FIG. 18 is a diagram showing the measurement results of the equator of a tomato using a high-precision practical device. 図19は、トマト花痕部の高精度実用機での計測結果を示す図である。FIG. 19 shows the results of measurement of the flower scar of a tomato using a high-precision practical machine. 図20は、トマト花痕部のLED試作機での計測結果を示す図である。FIG. 20 shows the measurement results of the tomato flower scar using the LED prototype. 図21は、イチゴの高精度実用機での計測結果を示す図である。FIG. 21 shows the results of measuring strawberries using a high-precision practical machine. 図22は、パプリカの高精度実用機での計測結果を示す図である。FIG. 22 is a diagram showing the measurement results of paprika using a high-precision practical machine. 図23は、メロンの高精度実用機での計測結果を示す図である。FIG. 23 is a diagram showing the measurement results of a melon using a high-precision practical machine. 図24は、メロンのLED試作機での計測結果を示す図である。FIG. 24 is a diagram showing the measurement results of a melon using an LED prototype. 図25は、ミカン花痕部の高精度実用機での計測結果を示す図である。FIG. 25 shows the results of measurement of the tangerine blossom scar using a high-precision practical machine. 図26は、ミカン赤道部のLED試作機での計測結果を示す図である。FIG. 26 shows the measurement results of the LED prototype at the equator of a mandarin orange. 図27は、不知火の高精度実用機での計測結果を示す図である。FIG. 27 shows the measurement results of the Shiranui high-precision practical instrument. 図28は、サクランボのLED試作機での計測結果を示す図である。FIG. 28 is a diagram showing the measurement results of cherries using the LED prototype. 図29は、トマト破砕物及びトマト加工品の高精度実用機での計測結果を示す図である。FIG. 29 shows the results of measurement of crushed tomatoes and processed tomato products using a high-precision practical machine.

以下に、本発明に係る糖度計測方法、糖度計測装置、及び糖度計測プログラムの好適な実施の形態の例を、図1~図29を参照して詳細に説明する。なお、この実施の形態によりこの発明が限定されるものではない。 Below, examples of preferred embodiments of the sugar content measurement method, sugar content measurement device, and sugar content measurement program according to the present invention will be described in detail with reference to Figs. 1 to 29. Note that the present invention is not limited to these embodiments.

[本発明の概略]
まず、図1を参照して、本発明の概略を説明する。図1は、本発明の概略を説明するための説明図である。
[Outline of the present invention]
First, an outline of the present invention will be described with reference to Fig. 1. Fig. 1 is an explanatory diagram for explaining an outline of the present invention.

上述したように、特許文献1では、860~890nm、900~920nm、920を越えて~960nmの3種類の光源と説明変数(吸光度)を使用しているが、波長域は広いため各波長域で有効な波長をさらに検討する必要があり、かつ、これらの3波長域を適用した場合に糖度の非破壊計測精度が低いという問題がある。 As mentioned above, Patent Document 1 uses three types of light sources, 860-890 nm, 900-920 nm, and over 920 to 960 nm, and explanatory variables (absorbance). However, because the wavelength ranges are wide, it is necessary to further consider the effective wavelengths in each wavelength range, and there is a problem that the accuracy of non-destructive measurement of sugar content is low when these three wavelength ranges are applied.

そこで、本発明では、3種類の波長域に加えて、更に、第4の波長域の使用した組み合わせを検討した。 Therefore, in this invention, in addition to the three wavelength ranges, we also considered a combination that uses a fourth wavelength range.

本発明では、説明変数の組合せとして、876±2nmまたは884±2nm、900~918nm、926±2nmの吸光度(光のシグナル値)に、第4の説明変数として856±2nmの吸光度(光のシグナル値)を組み合わせることで、高精度に糖度を非破壊計測できる点を見い出した。 In the present invention, it has been discovered that sugar content can be measured non-destructively with high accuracy by combining the absorbance (light signal value) at 876±2 nm or 884±2 nm, 900-918 nm, and 926±2 nm as a combination of explanatory variables with the absorbance (light signal value) at 856±2 nm as a fourth explanatory variable.

さらに、本発明では、リング状光源を使用し、リング状光源から青果物にリング状の光を照射し、青果物からの拡散反射を含む反射光をリング状光源の略中心で検出することで、より高精度に糖度を非破壊計測できることを見い出した。 Furthermore, in this invention, it has been discovered that sugar content can be measured non-destructively with higher accuracy by using a ring-shaped light source, irradiating a ring-shaped light from the ring-shaped light source onto the fruit or vegetable, and detecting the reflected light, including diffuse reflection, from the fruit or vegetable at approximately the center of the ring-shaped light source.

他方、新たな光源として普及しているLED(Light Emitting Diode)は、小型で長寿命、低消費電力であることが特徴である。従来、糖度の高精度な非破壊計測が可能な近赤外分光光度計(拡散反射光計測モードによる非接触計測(800-1000nm):(株)クボタ製K-BA100R(高精度実用機))では、光源としてリング状のハロゲンランプを用いて、リング状の光を試料に照射し、その拡散反射を含む反射光を光源の中心で検出し光ファイバーを用いて本体に伝達する。 On the other hand, LEDs (Light Emitting Diodes), which are becoming popular as a new light source, are characterized by their small size, long life, and low power consumption. Conventionally, near-infrared spectrophotometers capable of high-precision non-destructive measurement of sugar content (non-contact measurement using diffuse reflected light measurement mode (800-1000 nm): Kubota Corporation K-BA100R (high-precision practical model)) use a ring-shaped halogen lamp as a light source, irradiate the sample with a ring-shaped light, and detect the reflected light, including the diffuse reflection, at the center of the light source and transmit it to the main body using optical fibers.

しかしながら、LEDの普及によりハロゲンランプは品薄となる可能性がある。また、光ファイバーを使って導光する場合は、機器の体積もその分増加する。一方、MEMS(Micro Electro Mechanical Systems)技術により超小型かつ低コストな分光器が製造されるようになった。 However, with the widespread use of LEDs, halogen lamps may become scarce. Also, if optical fibers are used to guide the light, the volume of the device increases accordingly. On the other hand, MEMS (Micro Electro Mechanical Systems) technology has made it possible to manufacture ultra-compact, low-cost spectroscopes.

そこで、LEDリング光源の中央に超小型分光器を組み込んだ光ファイバーを用いない、分光器の真上を試料台とするコンパクトな近赤外分光光度計(LED試作機)を試作した。 Therefore, we developed a compact near-infrared spectrophotometer (LED prototype) that does not use optical fiber and has a super-compact spectrometer built into the center of an LED ring light source, with the sample stage placed directly above the spectrometer.

従来は近赤外分光法を用いて果実糖度を非破壊計測するために、非破壊計測された光吸収スペクトルに習慣的に前処理が実施されてきた。その前処理方法として主に、2次微分が多用されてきたが、前処理法は他に中心化や1次微分など多くの方法があり、微分の条件のみならず前処理法の組み合わせも検討しなければならない。また、小さな情報を再現性良く計測することにより高精度な非破壊計測が可能となっている(非特許文献1参照)が、光吸収スペクトルの前処理(微分、平滑化)はノイズを削除する効果があるものの、同時にシグナルの一部も削除してしまう。しかし、分光光度計の高性能化によりスペクトル前処理を適用する必要性は低くなってきている。 Conventionally, in order to measure fruit sugar content nondestructively using near-infrared spectroscopy, preprocessing has been customarily performed on the nondestructively measured light absorption spectrum. Second-order differentiation has been the most widely used preprocessing method, but there are many other preprocessing methods, such as centering and first-order differentiation, and it is necessary to consider not only the differentiation conditions but also the combination of preprocessing methods. In addition, high-precision nondestructive measurement is possible by measuring small information with good reproducibility (see Non-Patent Document 1), but while preprocessing of light absorption spectra (differentiation, smoothing) has the effect of removing noise, it also removes part of the signal at the same time. However, the need to apply spectral preprocessing has decreased due to the improvement in the performance of spectrophotometers.

非破壊計測用検量線(「推定モデル」や「モデル式」ともいう)の研究開発においては、PLS回帰分析のように多数の説明変数を用いる方法が主流となったがモデル式の完成度が見えにくく、重回帰分析のように少ない説明変数でモデル式を開発した場合はモデル式の構造も説明しやすく実用性が高い。 In the research and development of calibration curves for nondestructive measurement (also called "estimated models" or "model formulas"), methods that use a large number of explanatory variables, such as PLS regression analysis, have become mainstream, but it is difficult to see the completeness of the model formula. However, when a model formula is developed with a small number of explanatory variables, such as multiple regression analysis, the structure of the model formula is easier to explain and is highly practical.

また、近赤外光吸収スペクトルは温度の影響を受けることが知られており、実用上は試料温度の異なる果実を非破壊計測してモデル式を開発することにより温度の影響を受けにくくすることが行われている(非特許文献2参照)。 In addition, it is known that near-infrared light absorption spectra are affected by temperature, and in practice, models are developed by non-destructively measuring fruits with different sample temperatures to reduce the effect of temperature (see Non-Patent Document 2).

そこで、前述の高精度実用機とLED試作機を用いて、温度の異なる青果物を非破壊計測して取得した光吸収スペクトルを微分や平滑化する前処理無しに青果物の糖度を非破壊推定する説明変数の組合せを新たに見出した。 Therefore, using the aforementioned high-precision practical device and LED prototype device, we discovered a new combination of explanatory variables that can non-destructively estimate the sugar content of fruits and vegetables without pre-processing such as differentiating or smoothing the light absorption spectra obtained by non-destructively measuring fruits and vegetables at different temperatures.

未検討のスイカ等を除き、青果物の品目が異なっても同様の説明変数を採用可能であり、統一したハードウエア・ソフトウエア設計が可能となる。また、計測した近赤外光吸収スペクトルのシグナルを削ること無く有効に利用可能となる。 Except for watermelon, which was not examined, the same explanatory variables can be used for different types of fruit and vegetables, making it possible to design unified hardware and software. In addition, the signal of the measured near-infrared light absorption spectrum can be effectively used without being reduced.

高精度実用機は、小果実(およそ10g未満)では集光性が悪化するので治具を開発して集光性を改善する必要があった(非特許文献3参照)。加えて、メロンのような大きな果実や果皮の厚いミカンでは計測時間がおよそ300ms(ミリ秒)程度と長くなることがある。一方、LED試作機では、LEDの設置角度を平面に対して、10度以上、望ましくは32度(例えば、特許文献2では、計測用光が試料を長距離にわたって透過させるために10度に設定している)に設置することで、高精度実用機に比して、短い計測時間(5~30ms)でサクランボ、スモモ、モモ、ナシ、リンゴ、トマト、メロン、ミカン等の非破壊計測が可能となる。 The light-gathering ability of the high-precision practical machine deteriorates with small fruits (less than about 10 g), so it was necessary to develop a jig to improve the light-gathering ability (see Non-Patent Document 3). In addition, the measurement time can be as long as about 300 ms (milliseconds) for large fruits such as melons and mandarin oranges with thick skin. On the other hand, with the LED prototype machine, the installation angle of the LED is set at 10 degrees or more, preferably 32 degrees, with respect to the plane (for example, in Patent Document 2, the angle is set to 10 degrees so that the measurement light passes through the sample over a long distance), making it possible to perform non-destructive measurement of cherries, plums, peaches, pears, apples, tomatoes, melons, mandarin oranges, etc. in a shorter measurement time (5 to 30 ms) than with the high-precision practical machine.

本発明では、計測対象の青果物は、果実及び野菜を含み、また、当該青果物の赤道部、花痕部、加工品、及び破砕物を含むものとする。計測対象の青果物として、西洋ナシ、カキ、マンゴー、イチゴ、パプリカ、ミカン、不知火、トマト、サクランボ、モモ、ナシ、リンゴ、スモモ、及びメロンを計測する場合を一例に挙げて説明するが、本発明はこれらの青果物以外にも適用可能である。 In the present invention, the fruits and vegetables to be measured include the equator, flower scar, processed products, and crushed parts of the fruits and vegetables. The following examples are given of measuring pears, persimmons, mangoes, strawberries, peppers, mandarin oranges, shiranui, tomatoes, cherries, peaches, pears, apples, plums, and melons, but the present invention can be applied to other fruits and vegetables as well.

図1において、本発明は、対象の青果物の876±2nmまたは884±2nm、900~918nm、926±2nmの吸光度(光のシグナル値)の3つの説明変数に加えて、第4の説明変数である856±2nmの吸光度(光のシグナル値)から糖度を推定するための推定モデルを作成する推定モデル作成工程(S1)と、測定対象の青果物の876±2nmまたは884±2nm、900~918nm、926±2nm、856±2nmの光のシグナル値を測定し、測定した876±2nmまたは884±2nm、900~918nm、926±2nm、856±2nmの吸光度(光のシグナル値)を推定モデルに適用して糖度を推定する計測工程(S2)と、を備える。 In FIG. 1, the present invention includes an estimation model creation process (S1) for creating an estimation model for estimating sugar content from three explanatory variables, namely, absorbance (light signal value) at 876±2 nm or 884±2 nm, 900-918 nm, and 926±2 nm of the target fruit or vegetable, as well as a fourth explanatory variable, absorbance (light signal value) at 856±2 nm, and a measurement process (S2) for measuring the light signal values at 876±2 nm or 884±2 nm, 900-918 nm, 926±2 nm, and 856±2 nm of the target fruit or vegetable, and applying the measured absorbance (light signal value) at 876±2 nm or 884±2 nm, 900-918 nm, 926±2 nm, and 856±2 nm to the estimation model to estimate sugar content.

推定モデル作成工程(S1)では、対象の青果物・青果物加工品の876±2nmまたは884±2nm、900~918nm、926±2nm、856±2nmの光のシグナル値を測定して、吸光度(光のシグナル値)を多変量解析して青果物中の糖度を推定するための推定モデルを作成する。吸光度(光のシグナル値)に対しては必要によりデータ前処理を行う。多変量解析では、既知の青果物に対して876±2nmまたは884±2nm、900~918nm、926±2nm、856±2nmの光のシグナル値を測定し、吸光度(光のシグナル値)に対して多変量解析を行うことで推定モデルを作成する。 In the estimation model creation step (S1), the light signal values of 876±2 nm or 884±2 nm, 900-918 nm, 926±2 nm, and 856±2 nm of the target fresh produce or fresh produce processed product are measured, and the absorbance (light signal value) is subjected to multivariate analysis to create an estimation model for estimating the sugar content in the fresh produce. Data preprocessing is performed on the absorbance (light signal value) as necessary. In the multivariate analysis, light signal values of 876±2 nm or 884±2 nm, 900-918 nm, 926±2 nm, and 856±2 nm are measured for known fresh produce, and a multivariate analysis is performed on the absorbance (light signal value) to create an estimation model.

計測工程(S2)では、糖度を計測したい青果物について876±2nmまたは884±2nm、900~918nm、926±2nm、856±2nmの吸光度(光のシグナル値)に対しては、必要によりデータ前処理を行う。取得した吸光度(光のシグナル値)を推定モデルに適用して、糖度を推定する。このように、一旦、推定モデルを作成すると、測定対象の青果物の876±2nmまたは884±2nm、900~918nm、926±2nm、856±2nmの吸光度(光のシグナル値)を取得するだけで、その糖度を高精度に推定(計測)することが可能となる。 In the measurement step (S2), data preprocessing is performed as necessary on the absorbance (light signal value) at 876±2 nm, 884±2 nm, 900-918 nm, 926±2 nm, and 856±2 nm for the fruits and vegetables whose sugar content is to be measured. The acquired absorbance (light signal value) is applied to an estimation model to estimate the sugar content. In this way, once an estimation model is created, it is possible to estimate (measure) the sugar content with high accuracy simply by acquiring the absorbance (light signal value) at 876±2 nm, 884±2 nm, 900-918 nm, 926±2 nm, and 856±2 nm for the fruits and vegetables to be measured.

[糖度計測方法]
図2を参照し、本実施の形態に係る糖度計測方法について説明する。図2は、本実施の形態に係る糖度計測方法を説明するためのフローチャートである。
[Sugar content measurement method]
The sugar content measuring method according to the present embodiment will be described with reference to Fig. 2. Fig. 2 is a flow chart for explaining the sugar content measuring method according to the present embodiment.

図2に示すように、本実施の形態に係る糖度計測方法は、対象の青果物の876±2nmまたは884±2nm、900~918nm、926±2nm、856±2nmの吸光度(光のシグナル値)を多変量解析して、糖度を推定するための推定モデルを作成する推定モデル作成工程(S1)と、測定対象の青果物の876±2nmまたは884±2nm、900~918nm、926±2nm、856±2nmの光のシグナル値を測定し、吸光度(光のシグナル値)を推定モデルに適用して糖度を推定する計測工程(S2)とに大別される。 As shown in FIG. 2, the sugar content measurement method according to this embodiment is roughly divided into an estimation model creation process (S1) in which the absorbance (light signal value) of the target fruit or vegetable at 876±2 nm or 884±2 nm, 900-918 nm, 926±2 nm, and 856±2 nm is subjected to multivariate analysis to create an estimation model for estimating sugar content, and a measurement process (S2) in which the light signal values of the target fruit or vegetable at 876±2 nm or 884±2 nm, 900-918 nm, 926±2 nm, and 856±2 nm are measured, and the absorbance (light signal value) is applied to the estimation model to estimate sugar content.

多変量解析には、回帰分析等がある。回帰分析は、特定の数値を推定する回帰式(検量線)を作成するものであり、糖度を推定するのに使用することができる。回帰分析には、PLS回帰分析や重回帰分析などがあり、糖度の数値を推定するのに使用することができる。本実施の形態の多変量解析では、重回帰分析を使用した場合について説明するが、本発明はこれに限られず、判別分析、主成分分析、クラスター分析等を使用してもよい。なお、多変量解析の代わりに機械学習により、又は、多変量解析に機械学習を組み込んでモデルを作成してもよい。 Multivariate analysis includes regression analysis. Regression analysis creates a regression equation (calibration curve) that estimates a specific value, and can be used to estimate sugar content. Regression analysis includes PLS regression analysis and multiple regression analysis, and can be used to estimate sugar content values. In the multivariate analysis of this embodiment, a case where multiple regression analysis is used will be described, but the present invention is not limited to this, and discriminant analysis, principal component analysis, cluster analysis, etc. may also be used. Note that a model may be created by machine learning instead of multivariate analysis, or by incorporating machine learning into multivariate analysis.

測定対象の青果物は、果実及び野菜を含み、また、測定部位については、当該青果物の赤道部、花痕部、加工品、及び破砕物を含むものとする。計測対象の青果物は、例えば、西洋ナシ、カキ、マンゴー、イチゴ、パプリカ、ミカン、不知火、トマト、サクランボ、モモ、ナシ、リンゴ、スモモ、及びメロン等である。 The fruits and vegetables to be measured include fruits and vegetables, and the measurement areas include the equator, flower scar, processed products, and crushed products of the fruits and vegetables. Examples of the fruits and vegetables to be measured include pears, persimmons, mangoes, strawberries, peppers, mandarin oranges, shiranui, tomatoes, cherries, peaches, pears, apples, plums, and melons.

推定モデル作成工程では、まず、推定モデルを作成するために、測定対象の青果物を準備する(ステップS11)。測定対象の青果物を計測部の試料台にセットして、光源から青果物に近赤外光(ハロゲン光又は複数のLED光)を照射し、その拡散反射を含む反射光を検出して、近赤外光波長域の4つの異なる波長(例えば、876±2nmまたは884±2nm、900~918nm、926±2nm、856±2nm)の、又はこれら4つの異なる波長を含む波長域の光のシグナル値を分光検出装置で取得する(ステップS13)。ここで、ここで、「拡散反射光を含む反射光」は、拡散反射光のみの場合と、拡散反射光+拡散反射光以外の反射光(例えば、全反射光)の場合を含む。また、光のシグナル値は、反射光を光電変換した電気信号値(元データ)である。この光のシグナル値は、所定の変換式で吸光度に変換することができる。以下の前処理や多変量解析等の処理では、光のシグナル値をそのまま使用してもよく、また、吸光度に変換して使用してもよい。反射光は、試料表面で反射する光及び試料内部で拡散反射する光の両方を含むものである。光源は、リング状の光源としてもよい。リング状の光源から照射される光に対する青果物の拡散反射を含む反射光を、当該リング状の光源の略中心で検出することにしてもよい。リング状の光源は、ハロゲンランプで構成してもよい。また、リング状の光源は、発光波長の異なる複数のLEDで構成してもよい。当該複数のLEDは、計測対象の青果物側に平面に対して所定角度(10度より大きく、望ましくは32度)を有して配置されることにしてもよい。 In the estimation model creation process, first, the fruit or vegetable to be measured is prepared in order to create the estimation model (step S11). The fruit or vegetable to be measured is set on the sample stage of the measurement unit, and near-infrared light (halogen light or multiple LED lights) is irradiated from the light source onto the fruit or vegetable, and the reflected light including the diffuse reflection is detected, and the signal value of the light of four different wavelengths in the near-infrared wavelength range (e.g., 876±2 nm or 884±2 nm, 900-918 nm, 926±2 nm, 856±2 nm) or a wavelength range including these four different wavelengths is obtained by the spectroscopic detection device (step S13). Here, "reflected light including diffuse reflected light" includes the case of only diffuse reflected light and the case of diffuse reflected light + reflected light other than diffuse reflected light (e.g., total reflected light). The signal value of the light is an electrical signal value (original data) obtained by photoelectric conversion of the reflected light. This signal value of the light can be converted into absorbance using a predetermined conversion formula. In the following pre-processing, multivariate analysis, and other processing, the light signal value may be used as is, or may be converted to absorbance before use. Reflected light includes both light reflected from the sample surface and light diffusely reflected inside the sample. The light source may be a ring-shaped light source. Reflected light, including diffuse reflection from the fruit or vegetable relative to the light irradiated from the ring-shaped light source, may be detected approximately at the center of the ring-shaped light source. The ring-shaped light source may be a halogen lamp. The ring-shaped light source may also be composed of multiple LEDs with different emission wavelengths. The multiple LEDs may be arranged on the side of the fruit or vegetable to be measured at a predetermined angle (greater than 10 degrees, preferably 32 degrees) with respect to the plane.

つぎに、取得した吸光度(光のシグナル値)に対して、必要によりデータ前処理を実行する(ステップS14)。データ前処理では、例えば、中心化、標準化、規格化、及びベースライン補正等の1つ又は組み合わせて信号処理演算を行う。 Next, data preprocessing is performed on the acquired absorbance (light signal value) as necessary (step S14). In the data preprocessing, signal processing calculations are performed using, for example, one or a combination of centering, standardization, normalization, and baseline correction.

他方、測定対象の青果物について糖度の実測値を取得する(ステップS12)。ステップS15では、データ無処理又はデータ前処理が行われた近赤外光波長域の吸光度(光のシグナル値)を多変量解析して青果物中の糖度を推定するための推定モデルをそれぞれ作成する。 On the other hand, the actual sugar content of the fruits and vegetables being measured is obtained (step S12). In step S15, a multivariate analysis is performed on the absorbance (light signal value) in the near-infrared wavelength range for unprocessed or preprocessed data to create an estimation model for estimating the sugar content in each of the fruits and vegetables.

多変量解析として回帰分析を使用する場合は、例えば、分光吸光スペクトルから糖度の数値を推定するための回帰式(検量線)を推定モデルとして作成する。回帰分析としては、例えば、重回帰分析やPLS回帰分析を使用することができる。 When regression analysis is used as the multivariate analysis, for example, a regression equation (calibration curve) for estimating the sugar content value from the spectroscopic absorption spectrum is created as an estimation model. For example, multiple regression analysis or PLS regression analysis can be used as the regression analysis.

つぎに、計測工程では、糖度を推定したい青果物(測定対象物)を準備する(ステップS21)。測定対象の青果物を計測部の試料台にセットし、光源から青果物に近赤外光波長域の光(ハロゲン光又は波長域の異なる複数のLED光)を照射し、その拡散反射光を含む反射光を検出して、近赤外光波長域の4つの異なる波長(例えば、876±2nmまたは884±2nm、900~918nm、926±2nm、856±2nm)の、又はこれら4つの異なる波長を含む波長域の光のシグナル値を分光検出装置で測定して取得する(ステップS22)。ここで、「拡散反射光を含む反射光」は、拡散反射光のみの場合と、拡散反射光+拡散反射光以外の反射光(例えば、全反射光)の場合を含む。 Next, in the measurement process, the fruit or vegetable (measurement target) whose sugar content is to be estimated is prepared (step S21). The fruit or vegetable to be measured is set on the sample stage of the measurement unit, and light in the near-infrared wavelength range (halogen light or multiple LED lights with different wavelength ranges) is irradiated from the light source onto the fruit or vegetable. The reflected light including the diffuse reflected light is detected, and the signal values of the light at four different wavelengths in the near-infrared wavelength range (e.g., 876±2 nm or 884±2 nm, 900-918 nm, 926±2 nm, 856±2 nm) or in a wavelength range including these four different wavelengths are measured and obtained by a spectroscopic detection device (step S22). Here, "reflected light including diffuse reflected light" includes cases where there is only diffuse reflected light, and cases where there is diffuse reflected light plus reflected light other than diffuse reflected light (e.g., total reflected light).

つぎに、4つの異なる波長の吸光度(光のシグナル値)に対して、必要によりデータ前処理を実行する(ステップS23)。データ前処理は、S14と同様である。データ前処理が行われた4つの異なる波長の、又はこれら4つの異なる波長を含む波長域の吸光度(光のシグナル値)を、ステップS15で作成した糖度を推定するための推定モデルに適用して糖度を推定する(ステップS24)。 Next, data preprocessing is performed on the absorbance (light signal value) of the four different wavelengths as necessary (step S23). The data preprocessing is the same as in S14. The absorbance (light signal value) of the four different wavelengths or the wavelength range including these four different wavelengths that have been subjected to data preprocessing is applied to the estimation model for estimating sugar content created in step S15 to estimate the sugar content (step S24).

[糖度計測装置]
次に、本発明の糖度計測装置の構成について図3及び図4を参照し実施形態を例に挙げて説明する。なお、本実施の形態に係る糖度計測装置は、前述の糖度計測方法に好適に使用できるものであるが、本実施の形態に係る糖度計測方法に用いる装置はこれに限定されるものではない。
[Sugar content measuring device]
Next, the configuration of the sugar content measuring device of the present invention will be described by taking an embodiment as an example with reference to Figures 3 and 4. Note that the sugar content measuring device according to this embodiment can be suitably used in the sugar content measuring method described above, but the device used in the sugar content measuring method according to this embodiment is not limited to this.

ここで、図3は、本実施の形態に係る糖度装置の外観構成例を示す模式図である。図4は、本実施の形態に係る糖度計測装置の構成の一例を示すブロック図であり、該構成のうち本発明に関係する部分のみを概念的に示している。 Here, FIG. 3 is a schematic diagram showing an example of the external configuration of a sugar content device according to this embodiment. FIG. 4 is a block diagram showing an example of the configuration of a sugar content measuring device according to this embodiment, conceptually showing only the parts of the configuration that are related to the present invention.

図3及び図4に示すように、糖度計測装置1は、光源12を備えた計測部11、分光検出装置10と、データ処理装置20と、を備えている。 As shown in Figures 3 and 4, the sugar content measuring device 1 includes a measuring unit 11 having a light source 12, a spectroscopic detection device 10, and a data processing device 20.

データ処理装置20は、分光検出装置10で取得した近赤外線領域の波長の異なる4つの、又はこれら4つの異なる波長を含む波長域の分光吸光スペクトルから、測定対象の青果物の糖度を計測する装置である。データ処理装置20は、メモリ21、制御部23、計算処理部24を備えており、測定者はキーボード・マウス22により、データ処理装置20に測定条件等を入力する。分光検出装置10は、計測部11に接続されている。光源12は、計測部11に設けられている。光源12は、計測部11の試料台にセットされる測定対象の青果物に所定の波長(例えば、近赤外光波長域を含む波長)の光を照射する装置である。光源12は、リング状の光源としてもよい。リング状の光源から照射される光に対する青果物の拡散反射を含む反射光を、当該リング状の光源の略中心で検出することにしてもよい。リング状の光源は、ハロゲンランプで構成してもよい。また、リング状の光源は、発光波長の異なる複数のLEDで構成してもよい。当該複数のLEDは、計測対象の青果物側に平面に対して所定角度(10度より大きく、望ましくは32度)を有して配置されることにしてもよい。 The data processing device 20 is a device that measures the sugar content of the fruit or vegetable to be measured from the spectral absorption spectrum of four different wavelengths in the near-infrared region or a wavelength range including these four different wavelengths acquired by the spectroscopic detection device 10. The data processing device 20 is equipped with a memory 21, a control unit 23, and a calculation processing unit 24, and the measurer inputs the measurement conditions, etc. into the data processing device 20 using the keyboard/mouse 22. The spectroscopic detection device 10 is connected to the measurement unit 11. The light source 12 is provided in the measurement unit 11. The light source 12 is a device that irradiates light of a predetermined wavelength (for example, a wavelength including a near-infrared wavelength range) to the fruit or vegetable to be measured that is set on the sample stage of the measurement unit 11. The light source 12 may be a ring-shaped light source. The reflected light, including the diffuse reflection of the fruit or vegetable against the light irradiated from the ring-shaped light source, may be detected at approximately the center of the ring-shaped light source. The ring-shaped light source may be composed of a halogen lamp. The ring-shaped light source may also be composed of multiple LEDs with different emission wavelengths. The multiple LEDs may be positioned at a predetermined angle (greater than 10 degrees, preferably 32 degrees) with respect to the plane on the side of the fruit or vegetable being measured.

分光検出装置10は、計測部11に接続されている。分光検出装置10は、光源12から照射される光に対する測定対象の青果物・青果物加工品の反射光(直接反射光と拡散反射光を含む)を受光して近赤外光波長域の4つの異なる波長の吸光度(876±2nmまたは884±2nm、900~918nm、926±2nm、856±2nmの、又はこれら4つの異なる波長を含む波長域の吸光度(光のシグナル値))を取得し、データ処理装置20に送信する装置である。分光はグレーティング、フィルタ方式、前分光、後分光方式等の各種方式を使用することができる。また、光センサはイメージセンサ(ポリクロメーター)、フォトダイオード(モノクロメーター)等を使用することができる。 The spectroscopic detection device 10 is connected to a measurement unit 11. The spectroscopic detection device 10 receives the light (including direct reflected light and diffuse reflected light) reflected by the fruit or vegetable product to be measured in response to the light irradiated from the light source 12, obtains the absorbance (light signal value) of four different wavelengths in the near-infrared light wavelength range (876±2 nm or 884±2 nm, 900-918 nm, 926±2 nm, 856±2 nm, or a wavelength range including these four different wavelengths), and transmits it to a data processing device 20. Spectroscopic analysis can be performed using various methods such as grating, filter method, pre-spectroscopic method, and post-spectroscopic method. In addition, the optical sensor can be an image sensor (polychromator), photodiode (monochromator), etc.

計測部11は、光源12を備えており、試料台にセットされる測定対象に対して、光源12から入力される光源光を照射し、測定対象の青果物からの拡散反射を含む反射光を分光検出装置10に送る。 The measurement unit 11 is equipped with a light source 12, and irradiates the measurement object set on the sample stage with light source light input from the light source 12, and sends reflected light, including diffuse reflection, from the measurement object fruit or vegetable to the spectroscopic detection device 10.

データ処理装置20は、例えば、パーソナルコンピュータ等で構成することができ、メモリ21、制御部23、及び計算処理部24を備えており、キーボード・マウス22、I/Oポート(例えば、USBポート等)26、及びディスプレイ30等が接続されている。 The data processing device 20 can be configured, for example, as a personal computer, and is equipped with a memory 21, a control unit 23, and a calculation processing unit 24, and is connected to a keyboard/mouse 22, an I/O port (e.g., a USB port, etc.) 26, a display 30, etc.

メモリ21は、分光検出装置10からデータ処理装置20へ転送され、計算処理部24の分光吸光スペクトル取得部24-1により取得された分光吸光スペクトルや、計算処理部24の推定モデル作成部24-2で推定モデルを作成する際に使用する青果物の糖度の実測値や作成した推定モデル27等を格納する。推定モデル27を作成する際に使用する青果物の糖度の実測値は、キーボード・マウス22やI/Oポート26から入力することができる。 The memory 21 stores the spectroscopic absorption spectrum transferred from the spectroscopic detection device 10 to the data processing device 20 and acquired by the spectroscopic absorption spectrum acquisition unit 24-1 of the calculation processing unit 24, the actual measured sugar content of fruits and vegetables used when creating an estimation model in the estimation model creation unit 24-2 of the calculation processing unit 24, and the created estimation model 27. The actual measured sugar content of fruits and vegetables used when creating the estimation model 27 can be input from the keyboard/mouse 22 or the I/O port 26.

制御部23は、オペレータのキーボード・マウス22の操作に応じて、光源12に対する光の照射のON/OFFの指示や分光検出装置10のスペクトル検出の開始指示等の制御を行うことができ、また、計算処理部24に処理を行うよう命令することができる。 The control unit 23 can perform control operations such as instructing the light source 12 to turn on/off the light irradiation and instructing the spectroscopic detection device 10 to start spectral detection in response to the operator's operation of the keyboard/mouse 22, and can also command the calculation processing unit 24 to perform processing.

分光検出装置10から転送された青果物の4つの異なる波長の、又はこれら4つの異なる波長を含む波長域の分光吸光スペクトルは、データ処理装置20のメモリ21に格納される。測定者がキーボード・マウス22を通じて、計算処理部24に対して処理を行うよう命令すると、まず、計算処理部24の分光吸光スペクトル取得部24-1が、メモリ21に格納された分光吸光スペクトルから、必要によりデータ前処理により解析に必要な近赤外光波長域の分光吸光スペクトルを抽出する。 The spectral absorption spectrum of the four different wavelengths of the fruit or vegetable transferred from the spectroscopic detection device 10, or the wavelength range including these four different wavelengths, is stored in the memory 21 of the data processing device 20. When the measurer issues a command to the calculation processing unit 24 to perform processing via the keyboard/mouse 22, the spectral absorption spectrum acquisition unit 24-1 of the calculation processing unit 24 first extracts the spectral absorption spectrum in the near-infrared wavelength range required for analysis from the spectral absorption spectrum stored in the memory 21, by data preprocessing as necessary.

計算処理部24の推定モデル作成部24-2は、抽出された近赤外光波長域の分光吸光スペクトルに対して多変量解析を行い、推定モデルを作成する。本実施形態において、計算処理部24の推定モデル作成部24-2は、例えば、多変量解析としてPLS回帰分析や重回帰分析で作成した回帰式(検量線)を推定モデルとして作成してもよい。 The estimation model creation unit 24-2 of the calculation processing unit 24 performs multivariate analysis on the extracted spectroscopic absorption spectrum in the near-infrared light wavelength range to create an estimation model. In this embodiment, the estimation model creation unit 24-2 of the calculation processing unit 24 may create, for example, a regression equation (calibration curve) created by PLS regression analysis or multiple regression analysis as the multivariate analysis as the estimation model.

計算処理部24の推定モデル作成部24-2は、多変量解析の結果(推定モデル27)を、メモリ21に格納してもよく、ディスプレイ30上に出力してもよく、また、プリンタ(図示せず)を介して印刷してもよい。 The estimation model creation unit 24-2 of the calculation processing unit 24 may store the results of the multivariate analysis (estimation model 27) in the memory 21, output them on the display 30, or print them out via a printer (not shown).

計算処理部24の推定部24-3は、メモリ21に格納された推定モデル27に抽出された近赤外光波長域の4つの異なる波長の、又はこれら4つの異なる波長を含む波長域の分光吸光スペクトルを適用して、青果物の糖度を推定する。なお、メモリ21には、複数の種類の青果物の推定モデル27を格納しておくことで、1台の計測装置で、複数の種類の青果物の糖度を計測することが可能となる。なお、分光検出装置10に、データ処理装置20の機能を搭載して、分光検出装置10でデータ処理装置20の処理を実行して推定モデルの作成等を行ってもよい。 The estimation unit 24-3 of the calculation processing unit 24 estimates the sugar content of fruits and vegetables by applying the spectral absorption spectrum of the four different wavelengths in the near-infrared light wavelength range, or of a wavelength range including these four different wavelengths, extracted to the estimation model 27 stored in the memory 21. Note that by storing estimation models 27 for multiple types of fruits and vegetables in the memory 21, it becomes possible to measure the sugar content of multiple types of fruits and vegetables with a single measuring device. Note that the spectroscopic detection device 10 may be equipped with the functions of the data processing device 20, and the spectroscopic detection device 10 may execute the processing of the data processing device 20 to create the estimation model, etc.

[実施例]
図6~図29を参照して実施例を説明する。本実施例において、試料となる青果物として、西洋ナシ、カキ、マンゴー、イチゴ、パプリカ、ミカン、不知火、トマト、サクランボ、モモ、ナシ、リンゴ、スモモ、及びメロンを使用した。分光吸光スペクトル取得部24-1の処理により青果物の近赤外光波長域の4つの異なる波長の、又はこれら4つの異なる波長を含む波長域の分光吸光スペクトルを取得した。そして、推定モデル作成部24-2の処理により、青果物の近赤外光波長域の4つの異なる吸光度(シグナル値)を説明変数として重回帰分析を行うことで、青果物の糖度を推定するための重回帰式をそれぞれ作成した。また、重回帰式の評価を行った。さらに、推定部24-3により、作成した重回帰式を使用して、青果物の近赤外光波長域の4つの異なる吸光度(シグナル値)から当該青果物の糖度を推定した。
[Example]
An embodiment will be described with reference to FIG. 6 to FIG. 29. In this embodiment, as the fruits and vegetables to be samples, pear, persimmon, mango, strawberry, paprika, mandarin orange, Shiranui, tomato, cherry, peach, pear, apple, plum, and melon were used. The spectroscopic absorption spectrum acquisition unit 24-1 acquired the spectroscopic absorption spectrum of four different wavelengths in the near-infrared wavelength range of the fruits and vegetables, or the wavelength range including these four different wavelengths. Then, the estimation model creation unit 24-2 performed multiple regression analysis using the four different absorbances (signal values) in the near-infrared wavelength range of the fruits and vegetables as explanatory variables to create multiple regression equations for estimating the sugar content of the fruits and vegetables. In addition, the multiple regression equations were evaluated. Furthermore, the estimation unit 24-3 used the created multiple regression equations to estimate the sugar content of the fruits and vegetables from the four different absorbances (signal values) in the near-infrared wavelength range of the fruits and vegetables.

(1.糖度計測装置)
糖度計測装置1の一部又は全部について、拡散反射光非接触計測方式の近赤外分光光度計である、(1)高精度実用機と、(2)LED試作機を使用した。
(1. Sugar content measuring device)
For part or all of the sugar content measuring device 1, (1) a high-precision practical device and (2) an LED prototype device, which are near-infrared spectrophotometers using a diffuse reflection light non-contact measurement method, were used.

(1)高精度実用機
高精度実用機として、(株)クボタ製K-BA100Rを使用した。図5は、高精度実用機の要部の斜視図を示しており、図5(A)は計測部の概略の斜視図、図5(B)はカバーを示している。高精度実用機の計測部は、図5(A)に示すように、本体部41と、リング状の光源12(本体のハロゲンランプから光ファイバーを経由して試料に照射)と、試料をセットするクッションで構成される試料台42、光源12から照射される近赤外線領域の光に対する試料の拡散反射を含む反射光を受光して、測定光を本体に導出する受光ファイバー(光ファイバー)43と、等を備えている。また、高精度実用機の計測部には、図5(B)に示すカバー44が設けられており、計測時には取り外して使用する。
(1) High-precision practical instrument As the high-precision practical instrument, K-BA100R manufactured by Kubota Corporation was used. Figure 5 shows a perspective view of the main part of the high-precision practical instrument, with Figure 5(A) being a schematic perspective view of the measurement part and Figure 5(B) showing the cover. As shown in Figure 5(A), the measurement part of the high-precision practical instrument includes a main body part 41, a ring-shaped light source 12 (illuminating the sample from the halogen lamp of the main body via an optical fiber), a sample stage 42 made of a cushion on which the sample is set, and a light-receiving fiber (optical fiber) 43 that receives reflected light including diffuse reflection of the sample against the light in the near-infrared region irradiated from the light source 12 and outputs the measurement light to the main body. In addition, the measurement part of the high-precision practical instrument is provided with a cover 44 shown in Figure 5(B), which is removed when used for measurement.

トマト破砕物およびトマト加工品は、セル(例えば、コニカミノルタ製CR-A50)を使用して、セルにトマト破砕物およびトマト加工品を充填して、セルを試料台42にセットして計測した。受光したシグナルから吸光度を計算し、例えば、2nm間隔のデータになるよう線形補間した。 The crushed tomatoes and processed tomato products were measured using a cell (e.g., Konica Minolta CR-A50) by filling the cell with the crushed tomatoes and processed tomato products and setting the cell on the sample stage 42. The absorbance was calculated from the received signal and linearly interpolated to obtain data at intervals of, for example, 2 nm.

(2)LED試作機
LED試作機を試作した。図6は、LED試作機の要部を示す図であり、図6-Aは、
LED試作機を上から見た概略の平面図、図6-Bは、LED試作機の概略の断面図、図6-Cは、LED試作機の試料台に試料(例えば、トマト)をセットして、LEDから照射された光の経路を説明するための模式図、図6-Dは、LED試作機の概略の斜視図を示している。
(2) LED Prototype An LED prototype was produced. Figure 6 shows the main parts of the LED prototype. Figure 6-A shows the main parts of the LED prototype.
FIG. 6-A is a schematic plan view of the LED prototype seen from above, FIG. 6-B is a schematic cross-sectional view of the LED prototype, FIG. 6-C is a schematic diagram for explaining the path of light irradiated from the LED when a sample (e.g., a tomato) is set on the sample stage of the LED prototype, and FIG. 6-D is a schematic perspective view of the LED prototype.

図6-A~図6-Dにおいて、LED試作機は、リング状の本体部51と、本体部51の上部に設けられた複数のLEDで構成されるリング状の光源12と、本体部51の内側に設けられ、測定対象の試料がセットされる試料台52と、本体部51の光入射口53に入射される試料の反射光(拡散反射光)を受光する超小型分光器(分光検出装置)54と、を備えている。 In Figures 6-A to 6-D, the LED prototype is equipped with a ring-shaped main body 51, a ring-shaped light source 12 consisting of multiple LEDs provided on the upper part of the main body 51, a sample stage 52 provided inside the main body 51 on which the sample to be measured is set, and an ultra-compact spectroscope (spectroscopic detection device) 54 that receives the reflected light (diffuse reflected light) from the sample that is incident on the light entrance 53 of the main body 51.

光源12は、3つのLED(中心波長:850、870、940nm)をリング状に交互に90度に1個配置した構成となっており、直径は例えば36mmとした。この3つのLED(中心波長:850、870、940nm)により、約800nm~約1000nmの波長範囲で山なりの連続光を得ることができ、上述の4つの異なる分光吸光スペクトルの光のシグナル値を取得することができる。なお、3つのLEDの中心波長の例はこれに限られるものではない。各LEDの設置角度は平面(XY面)に対して、試料台52側にα(例えば、α=32度)度とした。 The light source 12 is configured with three LEDs (center wavelengths: 850, 870, 940 nm) arranged alternately at 90 degrees in a ring shape, with a diameter of, for example, 36 mm. These three LEDs (center wavelengths: 850, 870, 940 nm) can produce continuous light with a peak in the wavelength range of approximately 800 nm to approximately 1000 nm, and can obtain the signal values of light with the four different spectral absorption spectra described above. Note that the central wavelengths of the three LEDs are not limited to these examples. The installation angle of each LED is set to α (for example, α = 32 degrees) degrees toward the sample stage 52 with respect to the plane (XY plane).

試料台52は、ゴム板(例えば、厚さ3mm・直径24mm)の中央に直径8.5mmの穴を空け、穴の中心が光入射口53の中央になるように構成した。メロンでは、この試料台52の上にゴムワッシャ(例えば、内径8mm、外径18mm、厚さ1.5mm)を設置して計測してもよい。 The sample stage 52 is configured by drilling a hole of 8.5 mm diameter in the center of a rubber plate (e.g., 3 mm thick and 24 mm diameter) so that the center of the hole is the center of the light entrance 53. For melons, measurements may be taken by placing a rubber washer (e.g., inner diameter 8 mm, outer diameter 18 mm, thickness 1.5 mm) on the sample stage 52.

超小型分光器54として、浜松ホトニクス製C11708MA(640-1050nm、波長分解能最大20nm、画素数256)を用い、光入射口53から入射される測定光を検出する。 The ultra-compact spectrometer 54 is a Hamamatsu Photonics C11708MA (640-1050 nm, maximum wavelength resolution of 20 nm, number of pixels of 256) and detects the measurement light incident from the light entrance port 53.

図6-Cに示すように、上記構成のLED試作機では、試料台52の上に測定対象物(青果物)をセットする。光源12の各LEDは、水平面(XY面)に対して、本体部51から所定角度α(例えば、α=32度)で光を照射して、測定対象物からの反射光(拡散反射光)を、光入射口53を介して超小型分光器54で検出する。 As shown in Figure 6-C, in the LED prototype configured as above, the measurement object (fruit or vegetable) is placed on the sample stage 52. Each LED of the light source 12 irradiates light from the main body 51 at a predetermined angle α (for example, α = 32 degrees) with respect to the horizontal plane (XY plane), and the reflected light (diffuse reflected light) from the measurement object is detected by the ultra-compact spectrometer 54 via the light entrance port 53.

スペクトルの計測条件は、ゲインLow、平均回数3または1、計測回数1とし、レファレンス板はLabshere製SRS-99-020を使用した。ダークとレファレンスは5ms計測、試料は5ms(スモモ、サクランボ、トマト)では集光性が悪い時は、10~30msとした。また、青果物を非破壊計測して得られた光吸収スペクトルは補間なし(前処理なし)とした。 The spectral measurement conditions were low gain, 3 or 1 average count, 1 measurement count, and a Labshere SRS-99-020 reference plate. Dark and reference were measured for 5 ms, while samples were measured for 10-30 ms when light collection was poor (plums, cherries, tomatoes). Additionally, no interpolation (no pre-processing) was performed on the optical absorption spectra obtained by non-destructive measurement of fruit and vegetables.

(2.計測方法)
(2-1.試料、試料温度、計測部位)
試料温度を2段階(低温、室温)に変えて、室温で、高精度実用機またはLED試作機を用いて果実赤道部(モモ、ナシ、リンゴ、西洋ナシ(果実下部)、カキ、スモモ、マンゴー、ミカン、不知火、トマト、イチゴ、パプリカ)または花痕部(ミニトマトを含むトマト、メロン、ミカン、サクランボ)を試料台の光入射口に向かい合わせて計測した。トマトは品温を3段階に変えて非破壊計測する試験も実施した。
(2. Measurement Method)
(2-1. Sample, sample temperature, measurement location)
The sample temperature was changed to two stages (low temperature, room temperature) and at room temperature, the fruit equator (peach, pear, apple, European pear (lower part of fruit), persimmon, plum, mango, mandarin orange, Shiranui, tomato, strawberry, paprika) or flower scar (tomato including cherry tomato, melon, mandarin orange, cherry) was measured facing the light entrance of the sample stage using a high-precision practical machine or an LED prototype machine. A test was also conducted for non-destructive measurement of tomatoes by changing the product temperature to three stages.

モモ、ナシ、リンゴ、西洋ナシ、カキ、スモモ、マンゴー、メロン、パプリカ、ミカン(赤道部計測)は非破壊計測部位の、トマト、イチゴ、ミカン(花痕部計測)、サクランボは果実全体の、不知火は果実を縦に1/2に切断し、各糖度を破壊測定し目的変数とした。 For peaches, pears, apples, pears, persimmons, plums, mangoes, melons, peppers, and mandarins (equatorial measurements) non-destructive measurements were made; for tomatoes, strawberries, mandarins (flower scar measurements) and cherries, the whole fruit was measured; and for shiranui, the fruit was cut in half lengthwise, and the sugar content of each was measured destructively and used as the objective variable.

(2-2.近赤外線スペクトルの計測)
西洋ナシ、カキ、マンゴー、イチゴ、パプリカ、ミカン花痕部、不知火、トマト赤道部(およそ10g以上の果実)、トマト破砕物やトマト加工品(トマトジュース、トマトピューレ、トマトケチャップ)は高精度実用機のみで、ミカン赤道部、サクランボ、ミニトマトを含むトマトはLED試作機のみで非破壊計測した。モモ、ナシ、リンゴ、スモモ、メロンは高精度実用機及びLED試作機を用いて非破壊計測した。
(2-2. Measurement of near-infrared spectrum)
Pears, persimmons, mangoes, strawberries, paprika, mandarin orange blossoms, shiranui, tomato equator (fruits weighing 10g or more), crushed tomatoes, and tomato products (tomato juice, tomato puree, tomato ketchup) were measured non-destructively using only the high-precision practical machine, while mandarin orange equator, cherries, and tomatoes including cherry tomatoes were measured non-destructively using only the LED prototype machine. Peaches, pears, apples, plums, and melons were measured non-destructively using both the high-precision practical machine and the LED prototype machine.

(2-3.非破壊計測用検量線の開発および評価)
検量線を作成するための試料(糖度既知)を「非破壊計測法開発用試料」と称し、作成した検量線を評価するための試料(糖度未知)を「非破壊計測法評価用試料」と称する。非破壊計測法開発用試料により得られたデータをMicrosoft Excelを用いて重回帰分析し非破壊計測用検量線(重回帰式:推定モデル)を開発し、得られた検量線(重回帰式)を未知試料(非破壊計測法評価用試料)に適用して評価した。評価指標としてRMSE(Root Mean Squared Error)を計算した。
(2-3. Development and evaluation of calibration curves for non-destructive measurements)
The sample (with known sugar content) for creating the calibration curve is called the "nondestructive measurement method development sample", and the sample (with unknown sugar content) for evaluating the created calibration curve is called the "nondestructive measurement method evaluation sample". The data obtained from the nondestructive measurement method development sample was subjected to multiple regression analysis using Microsoft Excel to develop a nondestructive measurement calibration curve (multiple regression equation: estimation model), and the obtained calibration curve (multiple regression equation) was applied to an unknown sample (nondestructive measurement method evaluation sample) for evaluation. RMSE (Root Mean Squared Error) was calculated as an evaluation index.

RMSEの計算式(非破壊計測法評価用試料の非破壊計測精度)は次式(1)のようになる。 The formula for calculating RMSE (nondestructive measurement accuracy of samples used for evaluating nondestructive measurement methods) is given by the following formula (1):

Figure 0007688907000001
Figure 0007688907000001

但し、Y:非破壊計測法による推定糖度
X:実際の糖度
n:非破壊計測法評価用試料の試料数
Where Y is the estimated sugar content by the non-destructive measurement method, X is the actual sugar content, and n is the number of samples for evaluation of the non-destructive measurement method.

(3.計測結果)
(3-1.モモ)
白桃および黄桃の計測では、非破壊計測法開発用試料(高精度実用機n=208、LED試作機n=142)、非破壊計測法評価用試料(高精度実用機n=144、LED試作機n=167)とした。
(3. Measurement Results)
(3-1. Momo)
For measurements of white peaches and yellow peaches, samples were used for developing non-destructive measurement methods (high-precision practical machine n = 208, LED prototype machine n = 142) and for evaluating non-destructive measurement methods (high-precision practical machine n = 144, LED prototype machine n = 167).

図7は、モモの高精度実用機での計測結果を示している。図7(A)は、説明変数(各波長の吸光度)を、2つ(876、902nm)、3つ(876、902、926nm)、4つ(856、876、902、926nm)とした場合の重回帰分析結果を示している。図7(B)は、開発時の糖度の実測値と非破壊計測値(推定値)の相関(重回帰分析)を示す図である。図7(C)は、評価時の糖度の実測値と非破壊計測値(推定値)の相関を示す図である。図7(B)及び(C)において、横軸は、糖度の実測値[Brix%]、縦軸は、糖度の非破壊計測値(推定値)[Brix%]を示している。図7(D)は、4つの説明変数(各波長の吸光度)を異なる波長の組み合わせとした場合の重回帰分析結果及び評価結果を示している。 Figure 7 shows the measurement results of peaches using a high-precision practical device. Figure 7 (A) shows the results of multiple regression analysis when the explanatory variables (absorbance at each wavelength) are two (876, 902 nm), three (876, 902, 926 nm), and four (856, 876, 902, 926 nm). Figure 7 (B) shows the correlation (multiple regression analysis) between the actual sugar content value at the time of development and the non-destructive measurement value (estimated value). Figure 7 (C) shows the correlation between the actual sugar content value at the time of evaluation and the non-destructive measurement value (estimated value). In Figures 7 (B) and (C), the horizontal axis shows the actual sugar content value [Brix %], and the vertical axis shows the non-destructive measurement value (estimated value) of sugar content [Brix %]. Figure 7 (D) shows the results of multiple regression analysis and evaluation when the four explanatory variables (absorbance at each wavelength) are combinations of different wavelengths.

図7(A)に示すように、説明変数を4つとした場合は、2つや3つの場合に比して高い相関係数となった。 As shown in Figure 7 (A), when there were four explanatory variables, the correlation coefficient was higher than when there were two or three explanatory variables.

図7(B)に示すように、非破壊計測法開発用試料において856、876、902、926nmの吸光度を説明変数に採用した時に、切片10.9、係数625.6、-1298.7、886.5、-216.4、相関係数0.96となった。図7(C)に示す非破壊計測法評価用試料においてもRMSEは0.72と良好であった。図7(D)に示すように、4つの説明変数(各波長の吸光度)を異なる波長の組み合わせとした場合でも、相関係数及びRMSEともに良好な結果を示した。 As shown in Figure 7 (B), when absorbance at 856, 876, 902, and 926 nm was used as explanatory variables for the sample for nondestructive measurement method development, the intercept was 10.9, the coefficients were 625.6, -1298.7, 886.5, and -216.4, and the correlation coefficient was 0.96. The RMSE for the sample for nondestructive measurement method evaluation shown in Figure 7 (C) was also good at 0.72. As shown in Figure 7 (D), even when the four explanatory variables (absorbance at each wavelength) were used in combination with different wavelengths, both the correlation coefficient and RMSE showed good results.

図8は、モモのLED試作機での計測結果を示している。図8(A)は、説明変数(各波長の吸光度)を、2つ(876、909nm)、3つ(876、909、926nm)、4つ(857、876、909、926nm)とした場合の重回帰分析結果を示している。図8(B)は、開発時の糖度の実測値と非破壊計測値(推定値)の相関(重回帰分析)を示す図である。図8(C)は、評価時の糖度の実測値と非破壊計測値(推定値)の相関を示す図である。図8(B)及び(C)において、横軸は、糖度の実測値[Brix%]、縦軸は、糖度の非破壊計測値(推定値)[Brix%]を示している。図8(D)は、4つの説明変数(各波長の吸光度)を異なる波長の組み合わせとした場合の重回帰分析結果及び評価結果を示している。 Figure 8 shows the measurement results of the LED prototype for peaches. Figure 8 (A) shows the results of multiple regression analysis when the explanatory variables (absorbance at each wavelength) are two (876, 909 nm), three (876, 909, 926 nm), and four (857, 876, 909, 926 nm). Figure 8 (B) shows the correlation (multiple regression analysis) between the actual sugar content value at the time of development and the non-destructive measurement value (estimated value). Figure 8 (C) shows the correlation between the actual sugar content value at the time of evaluation and the non-destructive measurement value (estimated value). In Figures 8 (B) and (C), the horizontal axis shows the actual sugar content value [Brix %], and the vertical axis shows the non-destructive measurement value (estimated value) of sugar content [Brix %]. Figure 8 (D) shows the results of multiple regression analysis and evaluation when the four explanatory variables (absorbance at each wavelength) are combinations of different wavelengths.

図8(A)に示すように、説明変数を4つとした場合は、2つや3つの場合に比して高い相関係数となった。 As shown in Figure 8 (A), when there were four explanatory variables, the correlation coefficient was higher than when there were two or three explanatory variables.

図8(B)に示すように、非破壊計測法開発用試料において856.53(857)、876.35(876)、909.33(909)、926.42(926)nmの吸光度を説明変数として採用した時に、切片14.8、係数519.3、-1049、810.6、-279.8、相関係数0.94となった。図8(C)に示す非破壊計測法評価用試料においてもRMSEは0.80と良好であった。開発した複数の検量線から最終的に一つを選択する場合は、高精度実用機と同様に相関係数やRMSEのみならず、適用試料やバイアスなどに留意してX:Y=1:1のラインを決定する必要がある。 As shown in Figure 8 (B), when the absorbance at 856.53 (857), 876.35 (876), 909.33 (909), and 926.42 (926) nm was used as explanatory variables for the nondestructive measurement method development sample, the intercept was 14.8, the coefficients were 519.3, -1049, 810.6, and -279.8, and the correlation coefficient was 0.94. The RMSE for the nondestructive measurement method evaluation sample shown in Figure 8 (C) was also good at 0.80. When ultimately selecting one of the multiple calibration curves developed, it is necessary to determine the X:Y = 1:1 line by taking into consideration not only the correlation coefficient and RMSE, but also the applicable sample and bias, just as with the high-precision practical device.

図8(D)に示すように、4つの説明変数(各波長の吸光度)を異なる波長の組み合わせとした場合でも、相関係数及びRMSEともに良好な結果を示した。 As shown in Figure 8 (D), even when the four explanatory variables (absorbance at each wavelength) were combined with different wavelengths, good results were obtained for both the correlation coefficient and RMSE.

(3-2.ナシ)
赤ナシの計測では、非破壊計測法開発用試料(n=144)、赤ナシおよび青ナシを非破壊計測法評価用試料(n=96)とした。
(3-2. None)
For the measurements of red pears, samples (n=144) were used for developing nondestructive measurement methods, and red and green pears were used for evaluating nondestructive measurement methods (n=96).

図9は、ナシの高精度実用機での計測結果を示している。図9(A)は、説明変数(各波長の吸光度)を、2つ(876、910nm)、3つ(876、910、926nm)、4つ(856、876、910、926nm)とした場合の重回帰分析結果を示している。図9(B)は、開発時の糖度の実測値と非破壊計測値(推定値)の相関(重回帰分析)を示す図である。図9(C)は、評価時の糖度の実測値と非破壊計測値(推定値)の相関を示す図である。図9(B)及び(C)において、横軸は、糖度の実測値[Brix%]、縦軸は、糖度の非破壊計測値(推定値)[Brix%]を示している。 Figure 9 shows the results of measurements of pears using a high-precision practical device. Figure 9 (A) shows the results of multiple regression analysis when the explanatory variables (absorbance at each wavelength) are two (876, 910 nm), three (876, 910, 926 nm), and four (856, 876, 910, 926 nm). Figure 9 (B) shows the correlation (multiple regression analysis) between the actual sugar content measured during development and the non-destructive measurement value (estimated value). Figure 9 (C) shows the correlation between the actual sugar content measured during evaluation and the non-destructive measurement value (estimated value). In Figures 9 (B) and (C), the horizontal axis shows the actual sugar content value [Brix %], and the vertical axis shows the non-destructive measurement value (estimated value) of sugar content [Brix %].

図9(A)に示すように、説明変数を4つとした場合は、2つや3つの場合に比して高い相関係数となった。 As shown in Figure 9 (A), when there were four explanatory variables, the correlation coefficient was higher than when there were two or three explanatory variables.

図9(B)に示すように、非破壊計測法開発用試料において856、876、910、926nmの吸光度を説明変数に採用した時に、切片7.1、係数502.3、-1015.4、792.8、-288.6、相関係数0.96となった。図9(C)に示すように、非破壊計測法評価用試料においてもRMSEは0.58と良好であった。 As shown in Figure 9 (B), when the absorbance at 856, 876, 910, and 926 nm was used as an explanatory variable for the sample for nondestructive measurement method development, the intercept was 7.1, the coefficients were 502.3, -1015.4, 792.8, and -288.6, and the correlation coefficient was 0.96. As shown in Figure 9 (C), the RMSE for the sample for nondestructive measurement method evaluation was also good at 0.58.

図10は、ナシのLED試作機での計測結果を示している。図10(A)は、説明変数(各波長の吸光度)を、2つ(876、909nm)、3つ(876、909、926nm)、4つ(857、876、909、926nm)とした場合の重回帰分析結果を示している。図10(B)は、開発時の糖度の実測値と非破壊計測値(推定値)の相関(重回帰分析)を示す図である。図10(C)は、評価時の糖度の実測値と非破壊計測値(推定値)の相関を示す図である。図10(B)及び(C)において、横軸は、糖度の実測値[Brix%]、縦軸は、糖度の非破壊計測値(推定値)[Brix%]を示している。 Figure 10 shows the measurement results of the LED prototype for pears. Figure 10 (A) shows the results of multiple regression analysis when the explanatory variables (absorbance at each wavelength) are two (876, 909 nm), three (876, 909, 926 nm), and four (857, 876, 909, 926 nm). Figure 10 (B) shows the correlation (multiple regression analysis) between the actual sugar content measured during development and the non-destructive measurement value (estimated value). Figure 10 (C) shows the correlation between the actual sugar content measured during evaluation and the non-destructive measurement value (estimated value). In Figures 10 (B) and (C), the horizontal axis shows the actual sugar content value [Brix %], and the vertical axis shows the non-destructive measurement value (estimated value) of sugar content [Brix %].

図10(A)に示すように、説明変数を4つとした場合は、2つや3つの場合に比して高い相関係数となった。 As shown in Figure 10 (A), when there were four explanatory variables, the correlation coefficient was higher than when there were two or three explanatory variables.

図10(B)に示すように、非破壊計測法開発用試料において856.53(857)、876.35(876)、909.33(909)、926.42(926)nmの吸光度を説明変数として採用した時に、切片12.9、係数406.6、-877.7、725.7、-253.4、相関係数0.86となった。図10(C)に示すように、非破壊計測法評価用試料においてもRMSEは0.87と良好であった。 As shown in Figure 10 (B), when the absorbances at 856.53 (857), 876.35 (876), 909.33 (909), and 926.42 (926) nm were used as explanatory variables for the samples used in developing nondestructive measurement methods, the intercept was 12.9, the coefficients were 406.6, -877.7, 725.7, and -253.4, and the correlation coefficient was 0.86. As shown in Figure 10 (C), the RMSE for the samples used in evaluating nondestructive measurement methods was also good at 0.87.

(3-3.リンゴ)
赤・黄・青リンゴの計測では、非破壊計測法開発用試料(n=134)、非破壊計測法評価用試料(高精度実用機n=164、LED試作機n=132)とした。
(3-3. Apple)
For measurements of red, yellow, and green apples, samples were used for developing non-destructive measurement methods (n = 134) and for evaluating non-destructive measurement methods (high-precision practical machines, n = 164, LED prototype machines, n = 132).

図11は、リンゴの高精度実用機での計測結果を示している。図11(A)は、説明変数(各波長の吸光度)を、2つ(876、902nm)、3つ(876、902、926nm)、4つ(856、876、902、926nm)とした場合の重回帰分析結果を示している。図11(B)は、開発時の糖度の実測値と非破壊計測値(推定値)の相関(重回帰分析)を示す図である。図11(C)は、評価時の糖度の実測値と非破壊計測値(推定値)の相関を示す図である。図11(B)及び(C)において、横軸は、糖度の実測値[Brix%]、縦軸は、糖度の非破壊計測値(推定値)[Brix%]を示している。 Figure 11 shows the measurement results of apples using a high-precision practical device. Figure 11 (A) shows the results of multiple regression analysis when the explanatory variables (absorbance at each wavelength) are two (876, 902 nm), three (876, 902, 926 nm), and four (856, 876, 902, 926 nm). Figure 11 (B) shows the correlation (multiple regression analysis) between the actual sugar content value at the time of development and the non-destructive measurement value (estimated value). Figure 11 (C) shows the correlation between the actual sugar content value at the time of evaluation and the non-destructive measurement value (estimated value). In Figures 11 (B) and (C), the horizontal axis shows the actual sugar content value [Brix %], and the vertical axis shows the non-destructive measurement value (estimated value) of sugar content [Brix %].

図11(A)に示すように、説明変数を4つとした場合は、2つや3つの場合に比して高い相関係数となった。 As shown in Figure 11 (A), when there were four explanatory variables, the correlation coefficient was higher than when there were two or three explanatory variables.

図11(B)に示すように、非破壊計測法開発用試料において856、876、902、926nmの吸光度を説明変数に採用した時に、切片14.0、係数575.7、-1154.4、780.1、-201.1、n=134、相関係数0.96となった。図11(C)に示すように、非破壊計測法評価用試料においてもRMSEは0.48と良好であった。910nm(相関係数0.97)よりも902nm(相関係数0.96)の方がわずかに相関係数は下がるが、モモと同様に902nmの方が非破壊計測法評価用試料のスロープ(バイアス)を改善した。 As shown in Figure 11 (B), when the absorbance at 856, 876, 902, and 926 nm was used as an explanatory variable for the samples for nondestructive measurement method development, the intercept was 14.0, the coefficients were 575.7, -1154.4, 780.1, -201.1, n = 134, and the correlation coefficient was 0.96. As shown in Figure 11 (C), the RMSE for the samples for nondestructive measurement method evaluation was also good at 0.48. The correlation coefficient was slightly lower at 902 nm (correlation coefficient 0.96) than at 910 nm (correlation coefficient 0.97), but as with peaches, 902 nm improved the slope (bias) of the samples for nondestructive measurement method evaluation.

図12は、リンゴのLED試作機での計測結果を示している。図12(A)は、説明変数(各波長の吸光度)を、2つ(876、909nm)、3つ(876、909、926nm)、4つ(857、876、909、926nm)とした場合の重回帰分析結果を示している。図12(B)は、開発時の糖度の実測値と非破壊計測値(推定値)の相関(重回帰分析)を示す図である。図12(C)は、評価時の糖度の実測値と非破壊計測値(推定値)の相関を示す図である。図12(B)及び(C)において、横軸は、糖度の実測値[Brix%]、縦軸は、糖度の非破壊計測値(推定値)[Brix%]を示している。 Figure 12 shows the measurement results of apples using the LED prototype. Figure 12 (A) shows the results of multiple regression analysis when the explanatory variables (absorbance at each wavelength) are two (876, 909 nm), three (876, 909, 926 nm), and four (857, 876, 909, 926 nm). Figure 12 (B) shows the correlation (multiple regression analysis) between the actual sugar content value at the time of development and the non-destructive measurement value (estimated value). Figure 12 (C) shows the correlation between the actual sugar content value at the time of evaluation and the non-destructive measurement value (estimated value). In Figures 12 (B) and (C), the horizontal axis shows the actual sugar content value [Brix %], and the vertical axis shows the non-destructive measurement value (estimated value) of sugar content [Brix %].

図12(A)に示すように、説明変数を4つとした場合は、2つや3つの場合に比して高い相関係数となった。 As shown in Figure 12 (A), when there were four explanatory variables, the correlation coefficient was higher than when there were two or three explanatory variables.

図12(B)に示すように、非破壊計測法開発用試料において856.53(857)、876.35(876)、909.33(909)、926.42(926)nmの吸光度を説明変数として採用した時に、切片14.5、係数378.1、-771.8、609、-214.3、相関係数0.92となった。図12(C)に示すように、非破壊計測法評価用試料においてもRMSEは0.88と良好であった。非破壊計測法評価用試料でRMSEがやや大きくなった理由の一つは、近赤外光吸収スペクトル計測時の平均回数を3回から1回に減らしたためと考えられる。 As shown in Figure 12 (B), when the absorbance at 856.53 (857), 876.35 (876), 909.33 (909), and 926.42 (926) nm was used as the explanatory variables for the nondestructive measurement method development sample, the intercept was 14.5, the coefficients were 378.1, -771.8, 609, and -214.3, and the correlation coefficient was 0.92. As shown in Figure 12 (C), the RMSE for the nondestructive measurement method evaluation sample was also good at 0.88. One of the reasons why the RMSE was slightly larger for the nondestructive measurement method evaluation sample is thought to be because the number of averages when measuring the near-infrared light absorption spectrum was reduced from three to one.

(3-4.西洋ナシ)
西洋ナシの計測では、非破壊計測法開発用試料はn=120、非破壊計測法評価用試料はn=100とした。
(3-4. Pear)
For the measurements of pears, n = 120 samples were used for developing the nondestructive measurement method, and n = 100 samples were used for evaluating the nondestructive measurement method.

図13は、西洋ナシの高精度実用機での計測結果を示している。図13(A)は、説明変数(各波長の吸光度)を、2つ(874、904nm)、3つ(874、904、926nm)、4つ(856、874、904、926nm)とした場合の重回帰分析結果を示している。図13(B)は、開発時の糖度の実測値と非破壊計測値(推定値)の相関(重回帰分析)を示す図である。図13(C)は、評価時の糖度の実測値と非破壊計測値(推定値)の相関を示す図である。図13(B)及び(C)において、横軸は、糖度の実測値[Brix%]、縦軸は、糖度の非破壊計測値(推定値)[Brix%]を示している。 Figure 13 shows the measurement results of pears using a high-precision practical device. Figure 13 (A) shows the results of multiple regression analysis when the explanatory variables (absorbance at each wavelength) are two (874, 904 nm), three (874, 904, 926 nm), and four (856, 874, 904, 926 nm). Figure 13 (B) shows the correlation (multiple regression analysis) between the actual sugar content value at the time of development and the non-destructive measurement value (estimated value). Figure 13 (C) shows the correlation between the actual sugar content value at the time of evaluation and the non-destructive measurement value (estimated value). In Figures 13 (B) and (C), the horizontal axis shows the actual sugar content value [Brix %], and the vertical axis shows the non-destructive measurement value (estimated value) of sugar content [Brix %].

図13(A)に示すように、説明変数を4つとした場合は、2つや3つの場合に比して高い相関係数となった。 As shown in Figure 13 (A), when there were four explanatory variables, the correlation coefficient was higher than when there were two or three explanatory variables.

図13(B)に示すように、非破壊計測法開発用試料において856、874、904、926nmの吸光度を説明変数に採用した時に、切片14.1、係数609.2、-1086.6、681.2、-204.8、相関係数0.94となった。図13(C)に示すように、品種「オーロラ」は糖度を確認する際に汁液が得にくい果実であり、非破壊計測値(推定値)との間に誤差が大きい試料が認められたが、非破壊計測法評価用試料においてもRMSEは0.74と良好であった。 As shown in Figure 13 (B), when absorbance at 856, 874, 904, and 926 nm was used as an explanatory variable in the samples for developing the nondestructive measurement method, the intercept was 14.1, the coefficients were 609.2, -1086.6, 681.2, and -204.8, and the correlation coefficient was 0.94. As shown in Figure 13 (C), the Aurora variety is a fruit from which it is difficult to obtain juice when checking the sugar content, and samples with large errors between the nondestructive measurement value (estimated value) were found, but even in the samples for evaluating the nondestructive measurement method, the RMSE was good at 0.74.

(3-5.カキ)
カキの計測では、非破壊計測法開発用試料はn=76、非破壊計測法評価用試料はn=178とした。
(3-5. Oysters)
For the measurements of oysters, n = 76 samples were used for developing non-destructive measurement methods, and n = 178 samples were used for evaluating non-destructive measurement methods.

図14は、カキの高精度実用機での計測結果を示している。図14(A)は、説明変数(各波長の吸光度)を、2つ(886、918nm)、3つ(886、918、928nm)、(856、886、918nm)、4つ(856、886、918、928nm)とした場合の重回帰分析結果を示している。図14(B)及び(D)は、開発時の糖度の実測値と非破壊計測値(推定値)の相関(重回帰分析)を示す図である。図14(C)及び(E)は、評価時の糖度の実測値と非破壊計測値(推定値)の相関を示す図である。図14(B)、(C)、(D)、及び(E)において、横軸は、糖度の実測値[Brix%]、縦軸は、糖度の非破壊計測値(推定値)[Brix%]を示している。 Figure 14 shows the measurement results of persimmons using a high-precision practical device. Figure 14 (A) shows the results of multiple regression analysis when the explanatory variables (absorbance at each wavelength) are two (886, 918 nm), three (886, 918, 928 nm), (856, 886, 918 nm), and four (856, 886, 918, 928 nm). Figures 14 (B) and (D) show the correlation (multiple regression analysis) between the actual sugar content value and the non-destructive measurement value (estimated value) during development. Figures 14 (C) and (E) show the correlation between the actual sugar content value and the non-destructive measurement value (estimated value) during evaluation. In Figures 14 (B), (C), (D), and (E), the horizontal axis shows the actual sugar content value [Brix %], and the vertical axis shows the non-destructive measurement value (estimated value) of sugar content [Brix %].

図14(A)に示すように、説明変数を4つとした場合は、2つや3つの場合に比して高い相関係数となった。 As shown in Figure 14 (A), when there were four explanatory variables, the correlation coefficient was higher than when there were two or three explanatory variables.

図14(B)に示すように、非破壊計測法開発用試料において、856、886、918、928nmの吸光度を説明変数に採用した時に、切片15.4、係数243.3、-689.6、885、-439.5、相関係数0.97となった。図14(C)に示すように、硬い果実で比較的大きな誤差を発生することがあるが、渋柿にも適用でき、非破壊計測法評価用試料もRMSEは1.02と良好であった。なお、このように光吸収スペクトルが測定可能な場合にスペクトル前処理が有効な例を以下に示す。スペクトルを2次微分処理後に856、876、906、926nmの2次微分値を説明変数として採用すると、開発時の相関係数は0.98(図14(D)参照)、評価時のRMSEは0.86(図14(E)参照)へと非破壊計測値を改善できた。 As shown in Figure 14 (B), when the absorbance at 856, 886, 918, and 928 nm was used as an explanatory variable in the sample for developing the nondestructive measurement method, the intercept was 15.4, the coefficient was 243.3, -689.6, 885, and -439.5, and the correlation coefficient was 0.97. As shown in Figure 14 (C), a relatively large error may occur with hard fruits, but it can also be applied to astringent persimmons, and the RMSE of the sample for evaluating the nondestructive measurement method was good at 1.02. An example of how spectrum preprocessing is effective when the optical absorption spectrum can be measured is shown below. When the second derivative values at 856, 876, 906, and 926 nm after the second derivative processing of the spectrum were used as explanatory variables, the correlation coefficient during development was 0.98 (see Figure 14 (D)) and the RMSE during evaluation was 0.86 (see Figure 14 (E)), improving the nondestructive measurement values.

(3-6.スモモ)
スモモの計測では、非破壊計測法開発用試料は高精度実用機n=84、LED試作機n=72、非破壊計測法評価用試料はn=100とした。
(3-6. Plum)
For the measurements of plums, the samples for developing the nondestructive measurement method were a high-precision practical machine (n = 84), an LED prototype machine (n = 72), and the samples for evaluating the nondestructive measurement method (n = 100).

図15は、スモモの高精度実用機での計測結果を示している。図15(A)は、説明変数(各波長の吸光度)を、2つ(876、900nm)、3つ(876、900、926nm)、4つ(856、876、900、926nm)、(856、886、900、926nm)とした場合の重回帰分析結果を示している。図15(B)は、開発時の糖度の実測値と非破壊計測値(推定値)の相関(重回帰分析)を示す図である。図15(C)は、評価時の糖度の実測値と非破壊計測値(推定値)の相関を示す図である。図15(B)及び(C)において、横軸は、糖度の実測値[Brix%]、縦軸は、糖度の非破壊計測値(推定値)[Brix%]を示している。 Figure 15 shows the measurement results of plums using a high-precision practical device. Figure 15 (A) shows the results of multiple regression analysis when the explanatory variables (absorbance at each wavelength) are two (876, 900 nm), three (876, 900, 926 nm), four (856, 876, 900, 926 nm), and (856, 886, 900, 926 nm). Figure 15 (B) shows the correlation (multiple regression analysis) between the actual sugar content value at the time of development and the non-destructive measurement value (estimated value). Figure 15 (C) shows the correlation between the actual sugar content value at the time of evaluation and the non-destructive measurement value (estimated value). In Figures 15 (B) and (C), the horizontal axis shows the actual sugar content value [Brix %], and the vertical axis shows the non-destructive measurement value (estimated value) of sugar content [Brix %].

図15(A)に示すように、説明変数を4つとした場合は、2つや3つの場合に比して高い相関係数となった。 As shown in Figure 15 (A), when there were four explanatory variables, the correlation coefficient was higher than when there were two or three explanatory variables.

図15(B)に示すように、非破壊計測法開発用試料において、856、876、900、926nmの吸光度を説明変数に採用した時に、切片14.1、係数1522.3、-3098.5、1998.3、-421.9、相関係数0.97と良好となった。図15(C)に示すように、非破壊計測法評価用試料においてもRMSEは0.66と良好であった。 As shown in Figure 15 (B), when the absorbance at 856, 876, 900, and 926 nm was used as an explanatory variable for the sample for developing the nondestructive measurement method, the intercept was 14.1, the coefficients were 1522.3, -3098.5, 1998.3, and -421.9, and the correlation coefficient was 0.97, which were good results. As shown in Figure 15 (C), the RMSE for the sample for evaluating the nondestructive measurement method was also good at 0.66.

図16は、スモモのLED試作機での計測結果を示している。図16(A)は、説明変数(各波長の吸光度)を、2つ(886、911nm)、3つ(886、911、928 nm)、4つ(855、886、911、928nm)とした場合の重回帰分析結果を示している。図16(B)及び(D)は、開発時の糖度の実測値と非破壊計測値(推定値)の相関(重回帰分析)を示す図である。図16(C)及び(E)は、評価時の糖度の実測値と非破壊計測値(推定値)の相関を示す図である。図16(B)、(C)、(D)、及び(E)において、横軸は、糖度の実測値[Brix%]、縦軸は、糖度の非破壊計測値(推定値)[Brix%]を示している。 Figure 16 shows the measurement results of plums using the LED prototype. Figure 16 (A) shows the results of multiple regression analysis when the explanatory variables (absorbance at each wavelength) are two (886, 911 nm), three (886, 911, 928 nm), and four (855, 886, 911, 928 nm). Figures 16 (B) and (D) show the correlation (multiple regression analysis) between the actual sugar content measured during development and the non-destructive measurement value (estimated value). Figures 16 (C) and (E) show the correlation between the actual sugar content measured during evaluation and the non-destructive measurement value (estimated value). In Figures 16 (B), (C), (D), and (E), the horizontal axis shows the actual sugar content value [Brix %], and the vertical axis shows the non-destructive measurement value (estimated value) of sugar content [Brix %].

図16(A)に示すように、説明変数を4つとした場合は、2つや3つの場合に比して高い相関係数となった。 As shown in Figure 16 (A), when there were four explanatory variables, the correlation coefficient was higher than when there were two or three explanatory variables.

図16(B)に示すように、非破壊計測法開発用試料において854.53(855)、886.15(886)、911.25(911)、928.31(928)nmの吸光度を説明変数に採用した時に、切片9.4、係数339.4、-1462.5、1505.6、-393、相関係数0.91と良好となった。図16(C)に示すように、非破壊計測法評価用試料においてもRMSEは0.97と良好であった。高精度実用機と比較して相関係数が低くなった理由の一つは、近赤外光吸収スペクトル計測時の平均回数を1回に設定したためと考えられる。なお、このように光吸収スペクトルが測定可能な場合にスペクトル前処理が有効な例を以下に示す。非破壊計測法評価用試料では、非破壊計測値がやや高い傾向があるが(バイアス=0.395)、スペクトルをSNV(Standard Normal Variate)処理後に同じ説明変数を採用すると開発時の相関係数は0.90と若干低下するが(図16(D)参照)、評価時のRMSEは0.96、バイアスは0.281(図16(E)参照)へと非破壊計測値が高くなる傾向を改善できた。 As shown in Figure 16 (B), when the absorbance at 854.53 (855), 886.15 (886), 911.25 (911), and 928.31 (928) nm was used as the explanatory variable for the sample for developing the nondestructive measurement method, the intercept was 9.4, the coefficient was 339.4, -1462.5, 1505.6, and -393, and the correlation coefficient was 0.91, which were good results. As shown in Figure 16 (C), the RMSE was also good at 0.97 for the sample for evaluating the nondestructive measurement method. One of the reasons for the lower correlation coefficient compared to the high-precision practical machine is thought to be that the averaging number of times when measuring the near-infrared light absorption spectrum was set to one. An example in which spectrum preprocessing is effective when the light absorption spectrum can be measured like this is shown below. For samples used for evaluating nondestructive measurement methods, the nondestructive measurement values tended to be slightly higher (bias = 0.395), but when the same explanatory variables were used after the spectrum was processed with SNV (Standard Normal Variate), the correlation coefficient during development dropped slightly to 0.90 (see Figure 16 (D)), but the RMSE during evaluation was 0.96 and the bias was 0.281 (see Figure 16 (E)), improving the tendency for nondestructive measurement values to be higher.

(3-7.マンゴー)
マンゴーの計測では、非破壊計測法開発用試料はn=52、非破壊計測法評価用試料はn=34とした。
(3-7. Mango)
For mango measurements, n = 52 samples were used for developing the nondestructive measurement method, and n = 34 samples were used for evaluating the nondestructive measurement method.

図17は、マンゴーの高精度実用機での計測結果を示している。図17(A)は、説明変数(各波長の吸光度)を、2つ(876、904nm)、3つ(876、904、924nm)、4つ(858、876、904、924nm)とした場合の重回帰分析結果を示している。図17(B)は、開発時の糖度の実測値と非破壊計測値(推定値)の相関(重回帰分析)を示す図である。図17(C)は、評価時の糖度の実測値と非破壊計測値(推定値)の相関を示す図である。図17(B)及び(C)において、横軸は、糖度の実測値[Brix%]、縦軸は、糖度の非破壊計測値(推定値)[Brix%]を示している。 Figure 17 shows the results of measurements of mangoes using a high-precision practical device. Figure 17 (A) shows the results of multiple regression analysis when the explanatory variables (absorbance at each wavelength) are two (876, 904 nm), three (876, 904, 924 nm), and four (858, 876, 904, 924 nm). Figure 17 (B) shows the correlation (multiple regression analysis) between the actual sugar content measured during development and the non-destructive measurement value (estimated value). Figure 17 (C) shows the correlation between the actual sugar content measured during evaluation and the non-destructive measurement value (estimated value). In Figures 17 (B) and (C), the horizontal axis shows the actual sugar content value [Brix %], and the vertical axis shows the non-destructive measurement value (estimated value) of sugar content [Brix %].

図17(A)に示すように、説明変数を4つとした場合は、2つや3つの場合に比して高い相関係数となった。 As shown in Figure 17 (A), when there were four explanatory variables, the correlation coefficient was higher than when there were two or three explanatory variables.

図17(B)に示すように、非破壊計測法開発用試料において858、876、904、924nmの吸光度を説明変数に採用した時に、切片13.3、係数458、-893.8、639.5、-208、相関係数0.95であった。また、図17(C)に示すように、非破壊計測法評価用試料においてもRMSEは0.98と良好であった。 As shown in Figure 17 (B), when the absorbance at 858, 876, 904, and 924 nm was used as an explanatory variable for the sample for nondestructive measurement method development, the intercept was 13.3, the coefficients were 458, -893.8, 639.5, and -208, and the correlation coefficient was 0.95. Also, as shown in Figure 17 (C), the RMSE for the sample for nondestructive measurement method evaluation was good at 0.98.

(3-8.トマト)
(1)赤道部計測
トマトの赤道部計測では、試料温度は3段階、約10g以上のトマトを非破壊計測法開発用試料(n=94)、非破壊計測法評価用試料(n=177)とした。
(3-8. Tomato)
(1) Equatorial Measurement For measuring the equatorial part of tomatoes, the sample temperature was set at three levels, and tomatoes weighing approximately 10 g or more were used as samples for developing the nondestructive measurement method (n=94) and for evaluating the nondestructive measurement method (n=177).

図18は、トマト赤道部の高精度実用機での計測結果を示している。図18(A)は、説明変数(各波長の吸光度)を、2つ(876、900nm)、3つ(876、900、926nm)、4つ(856、876、900、926nm)とした場合の重回帰分析結果を示している。図18(B)は、開発時の糖度の実測値と非破壊計測値(推定値)の相関(重回帰分析)を示す図である。図18(C)は、評価時の糖度の実測値と非破壊計測値(推定値)の相関を示す図である。図18(B)及び(C)において、横軸は、糖度の実測値[Brix%]、縦軸は、糖度の非破壊計測値(推定値)[Brix%]を示している。 Figure 18 shows the measurement results of the tomato equator using a high-precision practical device. Figure 18 (A) shows the results of multiple regression analysis when the explanatory variables (absorbance at each wavelength) are two (876, 900 nm), three (876, 900, 926 nm), and four (856, 876, 900, 926 nm). Figure 18 (B) shows the correlation (multiple regression analysis) between the actual sugar content value at the time of development and the non-destructive measurement value (estimated value). Figure 18 (C) shows the correlation between the actual sugar content value at the time of evaluation and the non-destructive measurement value (estimated value). In Figures 18 (B) and (C), the horizontal axis shows the actual sugar content value [Brix %], and the vertical axis shows the non-destructive measurement value (estimated value) of sugar content [Brix %].

図18(A)に示すように、説明変数を4つとした場合は、2つや3つの場合に比して高い相関係数となった。 As shown in Figure 18 (A), when there were four explanatory variables, the correlation coefficient was higher than when there were two or three explanatory variables.

図18(B)に示すように、非破壊計測法開発用試料において856、876、900、926nmの吸光度を説明変数に採用した時に、切片8.5、係数1525.3、-3181.3、2079.5、-419.4、相関係数0.97であった。図18(C)に示すように、非破壊計測法評価用試料においてもRMSEは0.61と良好であった。 As shown in Figure 18 (B), when absorbance at 856, 876, 900, and 926 nm was used as explanatory variables for the samples for nondestructive measurement method development, the intercept was 8.5, the coefficients were 1525.3, -3181.3, 2079.5, and -419.4, and the correlation coefficient was 0.97. As shown in Figure 18 (C), the RMSE for the samples for nondestructive measurement method evaluation was also good at 0.61.

(2)花痕部計測
トマトの花痕部計測では、非破壊計測開発用試料は高精度実用機で試料温度3段階、約10g以上のトマト(n=72)、LED試作機で試料温度2段階、10g未満のトマトも含めてn=73とした。
(2) Measurement of the flower scar For measuring the flower scar of tomatoes, samples for non-destructive measurement development were prepared using a high-precision practical device with three sample temperature levels and tomatoes weighing approximately 10 g or more (n = 72), and an LED prototype device with two sample temperature levels and including tomatoes weighing less than 10 g (n = 73).

図19は、トマト花痕部の高精度実用機での計測結果を示している。図19(A)は、説明変数(各波長の吸光度)を、2つ(876、902nm)、3つ(876、902、926nm)、(856、876、902nm)、4つ(856、876、902、926nm)とした場合の重回帰分析結果を示している。図19(B)は、開発時の糖度の実測値と非破壊計測値(推定値)の相関(重回帰分析)を示す図である。図19(B)において、横軸は、糖度の実測値[Brix%]、縦軸は、糖度の非破壊計測値(推定値)[Brix%]を示している。 Figure 19 shows the results of measurements of the tomato flower scar using a high-precision practical device. Figure 19 (A) shows the results of multiple regression analysis when the explanatory variables (absorbance at each wavelength) are two (876, 902 nm), three (876, 902, 926 nm), (856, 876, 902 nm), and four (856, 876, 902, 926 nm). Figure 19 (B) shows the correlation (multiple regression analysis) between the actual sugar content value during development and the non-destructive measurement value (estimated value). In Figure 19 (B), the horizontal axis shows the actual sugar content value [Brix %], and the vertical axis shows the non-destructive measurement value (estimated value) of sugar content [Brix %].

図19(A)に示すように、説明変数を4つとした場合は、2つや3つの場合に比して高い相関係数となった。 As shown in Figure 19 (A), when there were four explanatory variables, the correlation coefficient was higher than when there were two or three explanatory variables.

図19(B)に示すように、非破壊計測法開発用試料において856、876、902、926nmの吸光度を説明変数に採用した時に、切片9.1、係数1456.6、-2917、1881.4、-422.3、相関係数0.95と良好であった。 As shown in Figure 19 (B), when the absorbance at 856, 876, 902, and 926 nm was used as an explanatory variable for the sample for developing a nondestructive measurement method, the intercept was 9.1, the coefficients were 1456.6, -2917, 1881.4, and -422.3, and the correlation coefficient was 0.95, which was good.

図20は、トマト花痕部のLED試作機での計測結果を示している。図20(A)は、説明変数(各波長の吸光度)を、2つ(882、907nm)、3つ(882、907、925nm)、4つ(855、882、907、925nm)、(855、878、907、925nm)とした場合の重回帰分析結果を示している。図20(B)は、開発時の糖度の実測値と非破壊計測値(推定値)の相関(重回帰分析)を示す図である。図20(B)において、横軸は、糖度の実測値[Brix%]、縦軸は、糖度の非破壊計測値(推定値)[Brix%]を示している。 Figure 20 shows the measurement results of the LED prototype on the tomato flower scar. Figure 20 (A) shows the results of multiple regression analysis when the explanatory variables (absorbance at each wavelength) are two (882, 907 nm), three (882, 907, 925 nm), four (855, 882, 907, 925 nm), and (855, 878, 907, 925 nm). Figure 20 (B) shows the correlation (multiple regression analysis) between the actual sugar content value and the non-destructive measurement value (estimated value) during development. In Figure 20 (B), the horizontal axis shows the actual sugar content value [Brix %], and the vertical axis shows the non-destructive measurement value (estimated value) of sugar content [Brix %].

図20(A)に示すように、説明変数を4つとした場合は、2つや3つの場合に比して高い相関係数となった。 As shown in Figure 20 (A), when there were four explanatory variables, the correlation coefficient was higher than when there were two or three explanatory variables.

図20(B)に示すように、非破壊計測法開発用試料において854.53(855)、882.24(882)、907.42(907)、924.54(925)nmの吸光度を説明変数に採用した時に、切片16.5、係数644.4、-2101.4、2149.7、-675.1、相関係数0.90と良好であった。高精度実用機と異なり10g未満のトマトも同じ計測時間(5ms)で計測可能である。糖度が低い場合は誤差が大きくなることがある。 As shown in Figure 20 (B), when the absorbance at 854.53 (855), 882.24 (882), 907.42 (907), and 924.54 (925) nm was used as the explanatory variable for the sample for developing the non-destructive measurement method, the intercept was 16.5, the coefficients were 644.4, -2101.4, 2149.7, and -675.1, and the correlation coefficient was 0.90, which was good. Unlike the high-precision practical machine, tomatoes weighing less than 10 g can also be measured in the same measurement time (5 ms). When the sugar content is low, the error may be large.

(3-9.イチゴ)
イチゴの計測では、非破壊計測法開発用試料はn=110、非破壊計測法評価用試料はn=150とした。
(3-9. Strawberries)
For the measurements of strawberries, n = 110 samples were used for developing the non-destructive measurement method, and n = 150 samples were used for evaluating the non-destructive measurement method.

図21は、イチゴの高精度実用機での計測結果を示している。図21(A)は、説明変数(各波長の吸光度)を、2つ(886、902nm)、3つ(886、902、926nm)、4つ(856、886、902、926nm)、(856、876、902、926nm)とした場合の重回帰分析結果を示している。図21(B)は、開発時の糖度の実測値と非破壊計測値(推定値)の相関(重回帰分析)を示す図である。図21(C)は、評価時の糖度の実測値と非破壊計測値(推定値)の相関を示す図である。図21(B)及び(C)において、横軸は、糖度の実測値[Brix%]、縦軸は、糖度の非破壊計測値(推定値)[Brix%]を示している。 Figure 21 shows the results of measurements of strawberries using a high-precision practical device. Figure 21 (A) shows the results of multiple regression analysis when the explanatory variables (absorbance at each wavelength) are two (886, 902 nm), three (886, 902, 926 nm), four (856, 886, 902, 926 nm), and (856, 876, 902, 926 nm). Figure 21 (B) shows the correlation (multiple regression analysis) between the actual sugar content measured during development and the non-destructive measurement value (estimated value). Figure 21 (C) shows the correlation between the actual sugar content measured during evaluation and the non-destructive measurement value (estimated value). In Figures 21 (B) and (C), the horizontal axis shows the actual sugar content value [Brix %], and the vertical axis shows the non-destructive measurement value (estimated value) of sugar content [Brix %].

図21(A)に示すように、説明変数を4つとした場合は、2つや3つの場合に比して高い相関係数となった。 As shown in Figure 21 (A), when there were four explanatory variables, the correlation coefficient was higher than when there were two or three explanatory variables.

図21(B)に示すように、非破壊計測法開発用試料において856、886、902、926nmの吸光度を説明変数に採用した時に、切片8.5、係数1525.3、-3181.3、2079.5、-419.4、相関係数0.94であった。図21(C)に示すように、非破壊計測法評価用試料においてもRMSEは0.49と良好であった。 As shown in Figure 21 (B), when absorbance at 856, 886, 902, and 926 nm was used as explanatory variables for samples used in developing nondestructive measurement methods, the intercept was 8.5, the coefficients were 1525.3, -3181.3, 2079.5, and -419.4, and the correlation coefficient was 0.94. As shown in Figure 21 (C), the RMSE for samples used in evaluating nondestructive measurement methods was also good at 0.49.

(3-10.パプリカ)
パプリカの計測では、非破壊計測法開発用試料はn=112、非破壊計測法評価用試料はn=174とした。
(3-10. Paprika)
For the measurements of paprika, n = 112 samples were used for developing the nondestructive measurement method, and n = 174 samples were used for evaluating the nondestructive measurement method.

図22は、パプリカの高精度実用機での計測結果を示している。図22(A)は、説明変数(各波長の吸光度)を、2つ(876、902nm)、3つ(876、902、926nm)、4つ(856、876、902、926nm)とした場合の重回帰分析結果を示している。図22(B)は、開発時の糖度の実測値と非破壊計測値(推定値)の相関(重回帰分析)を示す図である。図22(C)は、評価時の糖度の実測値と非破壊計測値(推定値)の相関を示す図である。図22(B)及び(C)において、横軸は、糖度の実測値[Brix%]、縦軸は、糖度の非破壊計測値(推定値)[Brix%]を示している。 Figure 22 shows the measurement results of paprika using a high-precision practical device. Figure 22 (A) shows the results of multiple regression analysis when the explanatory variables (absorbance at each wavelength) are two (876, 902 nm), three (876, 902, 926 nm), and four (856, 876, 902, 926 nm). Figure 22 (B) shows the correlation (multiple regression analysis) between the actual sugar content value at the time of development and the non-destructive measurement value (estimated value). Figure 22 (C) shows the correlation between the actual sugar content value at the time of evaluation and the non-destructive measurement value (estimated value). In Figures 22 (B) and (C), the horizontal axis shows the actual sugar content value [Brix %], and the vertical axis shows the non-destructive measurement value (estimated value) of sugar content [Brix %].

図22(A)に示すように、説明変数を4つとした場合は、2つや3つの場合に比して高い相関係数となった。 As shown in Figure 22 (A), when there were four explanatory variables, the correlation coefficient was higher than when there were two or three explanatory variables.

図22(B)に示すように、非破壊計測法開発用試料において856、876、902、926nmの吸光度を説明変数に採用した時に、切片7.9、係数1797.4、-3605.4、2311.3、-504.2、相関係数0.88であった。図22(C)に示すように、非破壊計測法評価用試料においてもRMSEは0.59と良好であった。 As shown in Figure 22 (B), when absorbance at 856, 876, 902, and 926 nm was used as an explanatory variable for the sample for developing the nondestructive measurement method, the intercept was 7.9, the coefficients were 1797.4, -3605.4, 2311.3, and -504.2, and the correlation coefficient was 0.88. As shown in Figure 22 (C), the RMSE for the sample for evaluating the nondestructive measurement method was also good at 0.59.

(3-11.メロン)
メロンの計測では、非破壊計測法開発用試料は高精度実用機n=37、LED試作機n=42、非破壊計測法評価用試料は高精度実用機n=52とした。
(3-11. Melon)
For the melon measurements, the samples for developing the nondestructive measurement method were the high-precision practical machine n = 37 and the LED prototype machine n = 42, and the samples for evaluating the nondestructive measurement method were the high-precision practical machine n = 52.

図23は、メロンの高精度実用機での計測結果を示している。図23(A)は、説明変数(各波長の吸光度)を、2つ(876、902nm)、3つ(876、902、926nm)、4つ(856、876、902、926nm)、(856、884、902、926nm)とした場合の重回帰分析結果を示している。図23(B)は、開発時の糖度の実測値と非破壊計測値(推定値)の相関(重回帰分析)を示す図である。図23(C)は、評価時の糖度の実測値と非破壊計測値(推定値)の相関を示す図である。図23(B)及び(C)において、横軸は、糖度の実測値[Brix%]、縦軸は、糖度の非破壊計測値(推定値)[Brix%]を示している。 Figure 23 shows the measurement results of melons using a high-precision practical device. Figure 23 (A) shows the results of multiple regression analysis when the explanatory variables (absorbance at each wavelength) are two (876, 902 nm), three (876, 902, 926 nm), four (856, 876, 902, 926 nm), and (856, 884, 902, 926 nm). Figure 23 (B) shows the correlation (multiple regression analysis) between the actual sugar content value at the time of development and the non-destructive measurement value (estimated value). Figure 23 (C) shows the correlation between the actual sugar content value at the time of evaluation and the non-destructive measurement value (estimated value). In Figures 23 (B) and (C), the horizontal axis shows the actual sugar content value [Brix %], and the vertical axis shows the non-destructive measurement value (estimated value) of sugar content [Brix %].

図23(A)に示すように、説明変数を4つとした場合は、2つや3つの場合に比して高い相関係数となった。 As shown in Figure 23 (A), when there were four explanatory variables, the correlation coefficient was higher than when there were two or three explanatory variables.

図23(B)に示すように、非破壊計測法開発用試料において856、876、902、926nmの吸光度を説明変数に採用した時に、切片9.8、係数575.8、-1223.7、850.5、-205.5、相関係数0.83であった。図23(C)に示すように、非破壊計測法評価用試料においてもRMSEは1.21と良好であった。 As shown in Figure 23 (B), when the absorbance at 856, 876, 902, and 926 nm was used as an explanatory variable for the sample for nondestructive measurement method development, the intercept was 9.8, the coefficients were 575.8, -1223.7, 850.5, and -205.5, and the correlation coefficient was 0.83. As shown in Figure 23 (C), the RMSE for the sample for nondestructive measurement method evaluation was also good at 1.21.

図24は、メロンのLED試作機での計測結果を示している。図24(A)は、説明変数(各波長の吸光度)を、2つ(886、904nm)、3つ(886、904、926nm)、4つ(857、886、904、926nm)、(857、874、904、926nm)とした場合の重回帰分析結果を示している。図24(B)は、開発時の糖度の実測値と非破壊計測値(推定値)の相関(重回帰分析)を示す図である。図24(B)において、横軸は、糖度の実測値[Brix%]、縦軸は、糖度の非破壊計測値(推定値)[Brix%]を示している。 Figure 24 shows the measurement results of melons using the LED prototype. Figure 24 (A) shows the results of multiple regression analysis when the explanatory variables (absorbance at each wavelength) are two (886, 904 nm), three (886, 904, 926 nm), four (857, 886, 904, 926 nm), and (857, 874, 904, 926 nm). Figure 24 (B) shows the correlation (multiple regression analysis) between the actual sugar content values and non-destructive measurement values (estimated values) during development. In Figure 24 (B), the horizontal axis shows the actual sugar content values [Brix %], and the vertical axis shows the non-destructive measurement values (estimated values) of sugar content [Brix %].

図24(A)に示すように、説明変数を4つとした場合は、2つや3つの場合に比して高い相関係数となった。 As shown in Figure 24 (A), when there were four explanatory variables, the correlation coefficient was higher than when there were two or three explanatory variables.

計測部位である花痕部と試料台の間に隙間が発生するメロンが存在したため、ゴムワッシャを試料台52に設置して計測した結果、図24(B)に示すように、856.53(857)、886.15(886)、903.58(904)、926.43(926)nmの吸光度を説明変数に採用した場合、切片13.1、係数234.4、-937.6、851.9、-143.5、相関係数0.84と良好な結果を得た。 Since there were melons where a gap occurred between the measurement site, the flower scar, and the sample stage, a rubber washer was placed on the sample stage 52 and measurements were taken. As a result, as shown in Figure 24 (B), when the absorbance at 856.53 (857), 886.15 (886), 903.58 (904), and 926.43 (926) nm was used as the explanatory variables, good results were obtained with an intercept of 13.1, coefficients of 234.4, -937.6, 851.9, and -143.5, and a correlation coefficient of 0.84.

(3-12.ミカン)
(1)花痕部計測
ミカン花痕部の計測では、非破壊計測法開発用試料はn=95とした。
(3-12. Tangerine)
(1) Measurement of the flower scar When measuring the mandarin flower scar, n = 95 samples were used to develop a non-destructive measurement method.

図25は、ミカン花痕部の高精度実用機での計測結果を示している。図25(A)は、説明変数(各波長の吸光度)を、2つ(874、906nm)、3つ(874、906、926nm)、4つ(856、874、906、926nm)とした場合の重回帰分析結果を示している。図25(B)は、開発時の糖度の実測値と非破壊計測値(推定値)の相関(重回帰分析)を示す図である。図25(B)において、横軸は、糖度の実測値[Brix%]、縦軸は、糖度の非破壊計測値(推定値)[Brix%]を示している。 Figure 25 shows the results of measurements of the tangerine blossom scar using a high-precision practical device. Figure 25(A) shows the results of multiple regression analysis when the explanatory variables (absorbance at each wavelength) are two (874, 906 nm), three (874, 906, 926 nm), and four (856, 874, 906, 926 nm). Figure 25(B) shows the correlation (multiple regression analysis) between the actual sugar content value during development and the non-destructive measurement value (estimated value). In Figure 25(B), the horizontal axis shows the actual sugar content value [Brix %], and the vertical axis shows the non-destructive measurement value (estimated value) of sugar content [Brix %].

図25(A)に示すように、説明変数を4つとした場合は、2つや3つの場合に比して高い相関係数となった。 As shown in Figure 25 (A), when there were four explanatory variables, the correlation coefficient was higher than when there were two or three explanatory variables.

図25(B)に示すように、非破壊計測法開発用試料において856、874、906、926nmの吸光度を説明変数に採用した時に、切片10.2、係数848.4、-1570.8、1086.8、-368.5、相関係数0.96と良好であった。 As shown in Figure 25 (B), when the absorbance at 856, 874, 906, and 926 nm was used as an explanatory variable for the sample for developing the nondestructive measurement method, the intercept was 10.2, the coefficients were 848.4, -1570.8, 1086.8, and -368.5, and the correlation coefficient was 0.96, which was good.

(2)赤道部計測
ミカン赤道部の計測では、非破壊計測法開発用試料はn=72とした。
(2) Equatorial measurement For the measurement of the equatorial part of mandarin oranges, the number of samples for developing a non-destructive measurement method was n=72.

図26は、ミカン赤道部のLED試作機での計測結果を示している。図26(A)は、説明変数(各波長の吸光度)を、2つ(878、907nm)、3つ(878、907、928nm)、4つ(857、878、907、928nm)とした場合の重回帰分析結果を示している。図26(B)は、開発時の糖度の実測値と非破壊計測値(推定値)の相関(重回帰分析)を示す図である。図26(B)において、横軸は、糖度の実測値[Brix%]、縦軸は、糖度の非破壊計測値(推定値)[Brix%]を示している。 Figure 26 shows the measurement results of the LED prototype at the equator of a mandarin orange. Figure 26 (A) shows the results of multiple regression analysis when the explanatory variables (absorbance at each wavelength) are two (878, 907 nm), three (878, 907, 928 nm), and four (857, 878, 907, 928 nm). Figure 26 (B) shows the correlation (multiple regression analysis) between the actual sugar content value during development and the non-destructive measurement value (estimated value). In Figure 26 (B), the horizontal axis shows the actual sugar content value [Brix %], and the vertical axis shows the non-destructive measurement value (estimated value) of sugar content [Brix %].

図26(A)に示すように、説明変数を4つとした場合は、2つや3つの場合に比して高い相関係数となった。 As shown in Figure 26 (A), when there were four explanatory variables, the correlation coefficient was higher than when there were two or three explanatory variables.

図26(B)に示すように、非破壊計測法開発用試料において、856.53(857)、878.32(878)、907.42(907)、928.31(928)nmの吸光度を説明変数に採用した時に、切片10.8、係数587.5、-1586.5、1384.6、-387.4、相関係数0.93と良好であった。 As shown in Figure 26 (B), when the absorbances at 856.53 (857), 878.32 (878), 907.42 (907), and 928.31 (928) nm were used as explanatory variables in the samples for developing nondestructive measurement methods, the intercept was 10.8, the coefficients were 587.5, -1586.5, 1384.6, and -387.4, and the correlation coefficient was 0.93, which was good.

(3-13.不知火)
不知火の計測では、非破壊計測法開発用試料はn=36とした。
(3-13. Shiranui)
In the measurement of Shiranui, the number of samples for developing non-destructive measurement methods was n=36.

図27は、不知火の高精度実用機での計測結果を示している。図27(A)は、説明変数(各波長の吸光度)を、2つ(876、908nm)、3つ(876、908、926nm)、(856、876、908nm)、4つ(856、876、908、926nm)とした場合の重回帰分析結果を示している。図27(B)は、開発時の糖度の実測値と非破壊計測値(推定値)の相関(重回帰分析)を示す図である。図27(B)において、横軸は、糖度の実測値[Brix%]、縦軸は、糖度の非破壊計測値(推定値)[Brix%]を示している。 Figure 27 shows the measurement results using Shiranui's high-precision practical device. Figure 27(A) shows the results of multiple regression analysis when the explanatory variables (absorbance at each wavelength) are two (876, 908 nm), three (876, 908, 926 nm), (856, 876, 908 nm), and four (856, 876, 908, 926 nm). Figure 27(B) shows the correlation (multiple regression analysis) between the actual sugar content value and the non-destructive measurement value (estimated value) during development. In Figure 27(B), the horizontal axis shows the actual sugar content value [Brix %], and the vertical axis shows the non-destructive measurement value (estimated value) of sugar content [Brix %].

図27(A)に示すように、説明変数を4つとした場合、2つや3つの場合に比して高い相関係数となった。 As shown in Figure 27 (A), when there were four explanatory variables, the correlation coefficient was higher than when there were two or three explanatory variables.

図27(B)に示すように、非破壊計測法開発用試料において856、876、908、926nmの吸光度を説明変数に採用した時に、切片5.8、係数183.5、-418.5、339、-102.7、相関係数0.82と良好であった。 As shown in Figure 27 (B), when the absorbance at 856, 876, 908, and 926 nm was used as an explanatory variable for the sample for developing a nondestructive measurement method, the intercept was 5.8, the coefficients were 183.5, -418.5, 339, and -102.7, and the correlation coefficient was 0.82, which was good.

(3-14.サクランボ)
サクランボの計測では、非破壊計測法開発用試料はn=50、非破壊計測法評価用試料はn=60とした。
(3-14. Cherries)
For the measurements of cherries, n = 50 samples were used for developing the nondestructive measurement method, and n = 60 samples were used for evaluating the nondestructive measurement method.

図28は、サクランボのLED試作機での計測結果を示している。図28(A)は、説明変数(各波長の吸光度)を、2つ(884、909nm)、3つ(884,909、926nm)、4つ(857、884、909、926nm)、(857、874、909、928nm)とした場合の重回帰分析結果を示している。図28(B)は、開発時の糖度の実測値と非破壊計測値(推定値)の相関(重回帰分析)を示す図である。図28(C)は、評価時の糖度の実測値と非破壊計測値(推定値)の相関を示す図である。図28(B)及び(C)において、横軸は、糖度の実測値[Brix%]、縦軸は、糖度の非破壊計測値(推定値)[Brix%]を示している。 Figure 28 shows the measurement results of cherries using the LED prototype. Figure 28 (A) shows the results of multiple regression analysis when the explanatory variables (absorbance at each wavelength) are two (884, 909 nm), three (884, 909, 926 nm), four (857, 884, 909, 926 nm), and (857, 874, 909, 928 nm). Figure 28 (B) shows the correlation (multiple regression analysis) between the actual sugar content measured during development and the non-destructive measurement value (estimated value). Figure 28 (C) shows the correlation between the actual sugar content measured during evaluation and the non-destructive measurement value (estimated value). In Figures 28 (B) and (C), the horizontal axis shows the actual sugar content value [Brix %], and the vertical axis shows the non-destructive measurement value (estimated value) of sugar content [Brix %].

図28(A)に示すように、説明変数を4つとした場合、2つや3つの場合に比して高い相関係数となった。 As shown in Figure 28 (A), when there were four explanatory variables, the correlation coefficient was higher than when there were two or three explanatory variables.

図28(B)に示すように、非破壊計測法開発用試料において856.53(857)、884.19(884)、909.33(909)、926.43(926)nmの吸光度を説明変数に採用した時に、切片17.1、係数925.7、-3121.7、3003.1、-819.5、相関係数0.96であった。図28(C)に示すように、非破壊計測法評価用試料においてもRMSEは1.46と良好であった。 As shown in Figure 28 (B), when the absorbances at 856.53 (857), 884.19 (884), 909.33 (909), and 926.43 (926) nm were used as explanatory variables for the samples used in developing nondestructive measurement methods, the intercept was 17.1, the coefficients were 925.7, -3121.7, 3003.1, -819.5, and the correlation coefficient was 0.96. As shown in Figure 28 (C), the RMSE for the samples used in evaluating nondestructive measurement methods was also good at 1.46.

(3-15.トマト破砕物及びトマト加工品)
トマト破砕物及びトマト加工品の計測では、計測法開発用試料はn=13、計測法評価用試料はn=31とした。
(3-15. Crushed tomatoes and processed tomato products)
For the measurements of crushed tomatoes and processed tomato products, the number of samples for measurement method development was n=13 and the number of samples for measurement method evaluation was n=31.

図29は、トマト破砕物およびトマト加工品の高精度実用機での計測結果を示している。図29(A)は、説明変数(各波長の吸光度)を、2つ(876、902nm)、3つ(876、902、926nm)、4つ(856、876、902、926nm)とした場合の重回帰分析結果を示している。図29(B)は、開発時の糖度の実測値と計測値(推定値)の相関(重回帰分析)を示す図である。図29(C)は、評価時の糖度の実測値と計測値(推定値)の相関を示す図である。図29(B)及び(C)において、横軸は、糖度の実測値[Brix%]、縦軸は、糖度の計測値(推定値)[Brix%]を示している。 Figure 29 shows the measurement results of crushed tomatoes and processed tomato products using a high-precision practical machine. Figure 29 (A) shows the results of multiple regression analysis when the explanatory variables (absorbance at each wavelength) are two (876, 902 nm), three (876, 902, 926 nm), and four (856, 876, 902, 926 nm). Figure 29 (B) shows the correlation (multiple regression analysis) between the actual sugar content value and the measured value (estimated value) during development. Figure 29 (C) shows the correlation between the actual sugar content value and the measured value (estimated value) during evaluation. In Figures 29 (B) and (C), the horizontal axis shows the actual sugar content value [Brix %], and the vertical axis shows the measured sugar content value (estimated value) [Brix %].

図29(A)に示すように、説明変数を4つとした場合、2つや3つの場合に比して高い相関係数となった。 As shown in Figure 29 (A), when there were four explanatory variables, the correlation coefficient was higher than when there were two or three explanatory variables.

図29(B)に示すように、計測法開発用試料において856、876、902、926nmの吸光度を説明変数に採用した時に、切片10.6、係数4681.3、-10265.9、7346.7、-1730.2、n=13、相関係数0.994であった。図29(C)に示すように、計測法評価用試料においてもRMSEは0.98と良好であった。 As shown in Figure 29 (B), when the absorbance at 856, 876, 902, and 926 nm was used as an explanatory variable in the measurement method development sample, the intercept was 10.6, the coefficients were 4681.3, -10265.9, 7346.7, and -1730.2, n = 13, and the correlation coefficient was 0.994. As shown in Figure 29 (C), the RMSE was also good at 0.98 in the measurement method evaluation sample.

(4.まとめ)
876±2nmまたは884±2nm、900~918nm、926±2nmの3つの吸光度を重回帰式の説明変数として採用した場合、相関係数が0.34~0.96と高い値も含むが低い値が多かった。加えて、これら3波長を説明変数として採用した場合に、カキ、トマト花痕部計測(高精度実用機)では2波長と同程度の相関係数であった。そこで、第4の説明変数として856±2nmの吸光度を採用した場合、相関係数は0.82~0.994と改善した。
(4. Summary)
When the three absorbances of 876±2 nm or 884±2 nm, 900-918 nm, and 926±2 nm were used as explanatory variables in a multiple regression equation, the correlation coefficients were 0.34-0.96, including high values, but mostly low values. In addition, when these three wavelengths were used as explanatory variables, the correlation coefficients were similar to those of the two wavelengths in persimmon and tomato flower scar measurements (high-precision practical machine). Therefore, when the absorbance of 856±2 nm was used as the fourth explanatory variable, the correlation coefficient improved to 0.82-0.994.

ダークやリファレンス値が一定の場合は、分光された光のシグナル値を吸光度等に変換せずに直接説明変数として利用できる。波長を絞り込んだこれら4つの説明変数を採用して果実糖度を非破壊推定する重回帰式を開発した場合、高精度実用機及びLED試作機共に一貫してそれらの係数の符号は波長の低い方から順にプラス、マイナス、プラス、マイナスとなる。特に炭水化物(糖)の吸収帯である900~918nmの吸光度の係数はプラスとなり吸光度が高いほど糖の含有量が高いことを意味する。加えて、トマト果実破砕物およびトマト加工品においても同様に推定精度の高い結果が得られたため、これらの説明変数は普遍性が高い。 When the dark and reference values are constant, the signal values of the dispersed light can be used directly as explanatory variables without being converted to absorbance, etc. When a multiple regression equation for non-destructively estimating fruit sugar content is developed using these four explanatory variables with narrowed-down wavelengths, the signs of the coefficients are consistently plus, minus, plus, minus, from lowest wavelength, for both the high-precision practical machine and the LED prototype machine. In particular, the coefficient for absorbance in the 900-918 nm absorption band for carbohydrates (sugar) is positive, meaning that the higher the absorbance, the higher the sugar content. In addition, similarly high estimation accuracy results were obtained for crushed tomato fruit and processed tomato products, so these explanatory variables are highly universal.

スモモ、サクランボ、10g未満の果実を含むトマト、イチゴなど小果実やカキ、メロンを対象とする場合、876±2nmよりも884±2nmの方が非破壊計測精度を改善可能、または両波長で同程度の非破壊計測精度であった。 When measuring small fruits such as plums, cherries, tomatoes and strawberries containing fruits weighing less than 10g, as well as persimmons and melons, 884±2 nm was found to provide better non-destructive measurement accuracy than 876±2 nm, or both wavelengths provided similar non-destructive measurement accuracy.

上述のように、高精度実用機、LED試作機共に一貫して同様の結果が得られるが、多少違いも認められる。結果に示したように、供試試料、分光器や平均回数等により採用する説明変数が多少ずれる可能性があり、高精度実用機では902nmを採用した場合にLED試作機では909nmが採用された事例がモモ、リンゴで認められ、スモモでは高精度実用機で900nm、LED試作機で911nmが採用された。また、カキでは918nmが採用された。このように、糖の吸収帯では900~918nmと広い範囲をカバーする必要があった。一方、他の波長(856、876、884、926nm)では結果に示したように±2nm程度ずれる可能性がある。 As mentioned above, both the high-precision practical machine and the LED prototype machine consistently give similar results, but some differences are also observed. As shown in the results, there is a possibility that the explanatory variables used may vary slightly depending on the test sample, spectrometer, number of averages, etc., and there were cases for peaches and apples where 902 nm was used with the high-precision practical machine but 909 nm was used with the LED prototype machine, and for plums, 900 nm was used with the high-precision practical machine and 911 nm was used with the LED prototype machine. Also, 918 nm was used for persimmons. In this way, it was necessary to cover a wide range of 900 to 918 nm for the sugar absorption band. On the other hand, there is a possibility of deviation of about ±2 nm for other wavelengths (856, 876, 884, 926 nm), as shown in the results.

高精度実用機を用いて小さい果実(およそ10g未満)やメロン、ミカンを非破壊計測する場合に集光性が悪くなり計測時間がおよそ300~350ms程度と長くなる。一方、LED試作機では、前述のようにLEDの設置角度を平面に対して32度とすることで、集光性が向上し、短い計測時間(5~30ms)でモモ、ナシ、リンゴのみならず、サクランボ、スモモ、ミニトマトを含むトマト、メロン、ミカンの計測が可能である。 When using a high-precision practical device to non-destructively measure small fruits (approximately less than 10g), melons, or mandarins, the light gathering ability is poor and the measurement time is long, at approximately 300 to 350 ms. On the other hand, with the LED prototype, the LED is installed at an angle of 32 degrees to the plane as mentioned above, improving the light gathering ability and making it possible to measure not only peaches, pears, and apples, but also tomatoes, melons, and mandarins, including cherries, plums, and cherry tomatoes, in a short measurement time (5 to 30 ms).

加えて、LED試作機では、サクランボやミニトマト等の小果実で高精度実用機とは異なり専用の治具を用いずに短い時間で計測可能である。メロンではその形状により試料台と果実花痕部との間に隙間が発生する場合があったため、試料台にゴムワッシャを設置して外乱光の影響が少なくなるように統一して非破壊計測する必要がある。ミカンの花痕部においては、メロン花痕部よりもへこみの程度が甚だしいものがある。結果に示したように、ミカンの花痕部計測においては高精度実用機の方が適していたが、LED試作機を用いるミカンの花痕部計測においても、このような方法で試料台と試料の間の隙間を減らすことにより、精度の高い非破壊計測が可能である。 In addition, the LED prototype can measure small fruits such as cherries and cherry tomatoes in a short time without using special tools, unlike the high-precision practical machine. In melons, gaps can sometimes occur between the sample stage and the blossom scar of the fruit due to their shape, so it is necessary to install rubber washers on the sample stage and perform standardized non-destructive measurements to reduce the effects of ambient light. The blossom scar of mandarin oranges has more severe dents than the melon blossom scar. As the results show, the high-precision practical machine is more suitable for measuring the blossom scar of mandarin oranges, but by reducing the gap between the sample stage and the sample in this way, highly accurate non-destructive measurements are also possible when measuring the blossom scar of mandarin oranges using the LED prototype.

本発明によれば、876±2nmまたは884±2nm、900~918nm、926±2nmの吸光度に第4の説明変数として856±2nmの吸光度(ダークやリファレンス値が一定の場合は、分光されたシグナル値を直接説明変数として利用できる)を組み合わせることで、高精度に糖度を計測することが可能となる。 According to the present invention, by combining absorbance at 876±2 nm, 884±2 nm, 900-918 nm, and 926±2 nm with absorbance at 856±2 nm as a fourth explanatory variable (when the dark or reference value is constant, the spectroscopic signal value can be used directly as an explanatory variable), it is possible to measure sugar content with high accuracy.

本発明は、果実糖度の非破壊計測だけでなく、トマト破砕物やトマト加工品(トマトジュース、トマトピューレ、トマトケチャップ)の糖度の推定にも利用でき、適用試料はこれらの試料に限定されない。 The present invention can be used not only for non-destructive measurement of fruit sugar content, but also for estimating the sugar content of crushed tomatoes and tomato products (tomato juice, tomato puree, tomato ketchup), and the applicable samples are not limited to these.

[他の実施の形態]
さて、これまで本発明の実施の形態について説明したが、本発明は、上述した実施の形態以外にも、特許請求の範囲に記載した技術的思想の範囲内において種々の異なる実施の形態にて実施されてよいものである。
[Other embodiments]
While the embodiments of the present invention have been described above, the present invention may be embodied in various different embodiments other than those described above within the scope of the technical ideas set forth in the claims.

また、実施の形態において説明した各処理のうち、自動的に行われるものとして説明した処理の全部又は一部を手動的に行うこともでき、或いは、手動的に行われるものとして説明した処理の全部又は一部を公知の方法で自動的に行うこともできる。 Furthermore, among the processes described in the embodiments, all or part of the processes described as being performed automatically can be performed manually, or all or part of the processes described as being performed manually can be performed automatically using a known method.

このほか、上記文献中や図面中で示した処理手順、制御手順、具体的名称、各処理の登録データや検索条件等のパラメータを含む情報、画面例、データベース構成については、特記する場合を除いて任意に変更することができる。 In addition, the processing procedures, control procedures, specific names, registered data for each process, information including search conditions and other parameters, screen examples, and database configurations shown in the above documents and drawings may be changed as desired unless otherwise specified.

また、糖度計測装置1に関して、図示の各構成要素は機能概念的なものであり、必ずしも物理的に図示の如く構成されていることを要しない。 In addition, with regard to the sugar content measuring device 1, each component shown in the figure is a functional concept, and does not necessarily have to be physically configured as shown in the figure.

例えば、糖度計測装置1の各装置が備える処理機能、特に計算処理部24にて行われる各処理機能については、その全部又は任意の一部を、CPU(Central Processing Unit)及び当該CPUにて解釈実行されるプログラムにて実現してもよく、また、ワイヤードロジックによるハードウェアとして実現してもよい。尚、プログラムは、後述する記録媒体に記録されており、必要に応じて糖度計測装置1に機械的に読み取られる。すなわち、ROM又はHDなどのメモリ21などは、OS(Operating System)として協働してCPUに命令を与え、各種処理を行うためのコンピュータプログラムが記録されている。このコンピュータプログラムは、RAMにロードされることによって実行され、CPUと協働して制御部を構成する。 For example, all or any part of the processing functions of each device of the sugar content measuring device 1, particularly the processing functions performed by the calculation processing unit 24, may be realized by a CPU (Central Processing Unit) and a program interpreted and executed by the CPU, or may be realized as hardware using wired logic. The program is recorded on a recording medium described below, and is mechanically read by the sugar content measuring device 1 as necessary. That is, memory 21 such as a ROM or HD records a computer program that works together as an OS (Operating System) to give instructions to the CPU and perform various processes. This computer program is executed by being loaded into RAM, and works together with the CPU to form a control unit.

また、このコンピュータプログラムは、糖度計測装置1に対して任意のネットワークを介して接続されたアプリケーションプログラムサーバに記憶されていてもよく、必要に応じてその全部又は一部をダウンロードすることも可能である。 This computer program may also be stored in an application program server connected to the sugar content measuring device 1 via any network, and all or part of it may be downloaded as needed.

また、本発明に係るプログラムを、コンピュータ読み取り可能な記録媒体に格納してもよく、また、プログラム製品として構成することもできる。ここで、この「記録媒体」とは、メモリーカード、USBメモリ、SDカード、フレキシブルディスク、光磁気ディスク、ROM、EPROM、EEPROM、CD-ROM、MO、DVD、及び、Blu-ray(登録商標) Disc等の任意の「可搬用の物理媒体」を含むものとする。 The program according to the present invention may be stored on a computer-readable recording medium, or may be configured as a program product. Here, the term "recording medium" includes any "portable physical medium" such as a memory card, USB memory, SD card, flexible disk, magneto-optical disk, ROM, EPROM, EEPROM, CD-ROM, MO, DVD, and Blu-ray (registered trademark) Disc.

また、「プログラム」とは、任意の言語や記述方法にて記述されたデータ処理方法であり、ソースコードやバイナリコード等の形式を問わない。なお、「プログラム」は必ずしも単一的に構成されるものに限られず、複数のモジュールやライブラリとして分散構成されるものや、OS(Operating System)に代表される別個のプログラムと協働してその機能を達成するものをも含む。なお、実施の形態に示した各装置において記録媒体を読み取るための具体的な構成、読み取り手順、あるいは、読み取り後のインストール手順等については、周知の構成や手順を用いることができる。 A "program" is a data processing method written in any language or description method, and may be in any format, such as source code or binary code. Note that a "program" is not necessarily limited to a single configuration, but also includes a distributed configuration consisting of multiple modules or libraries, or a program that works in conjunction with a separate program, such as an OS (Operating System), to achieve its function. Note that the specific configurations and reading procedures for reading a recording medium in each device shown in the embodiments, as well as installation procedures after reading, can use well-known configurations and procedures.

メモリ21に格納される各種のデータベース等は、RAM、ROM等のメモリ装置、ハードディスク等の固定ディスク装置、フレキシブルディスク、光ディスク等のストレージ手段であり、各種処理やウェブサイト提供に用いる各種のプログラムやテーブルやデータベースやウェブページ用ファイル等を格納する。 The various databases stored in memory 21 are storage devices such as memory devices such as RAM and ROM, fixed disk devices such as hard disks, flexible disks, optical disks, etc., and store various programs, tables, databases, web page files, etc. used for various processes and website provision.

また、糖度計測装置1は、既知のパーソナルコンピュータ、ワークステーション等の情報処理装置として構成してもよく、また、該情報処理装置に任意の周辺装置を接続して構成してもよい。また、糖度計測装置1は、該情報処理装置に本発明の方法を実現させるソフトウェア(プログラム、データ等を含む)を実装することにより実現してもよい。 The sugar content measuring device 1 may be configured as a known information processing device such as a personal computer or a workstation, or may be configured by connecting any peripheral device to the information processing device. The sugar content measuring device 1 may also be realized by implementing software (including programs, data, etc.) that realizes the method of the present invention in the information processing device.

更に、装置の分散・統合の具体的形態は図示するものに限られず、その全部又は一部を、各種の付加等に応じて、又は、機能負荷に応じて、任意の単位で機能的又は物理的に分散・統合して構成することができる。すなわち、上述した実施形態を任意に組み合わせて実施してもよく、実施形態を選択的に実施してもよい。 Furthermore, the specific form of distribution and integration of the devices is not limited to that shown in the figures, and all or part of them can be functionally or physically distributed and integrated in any unit depending on various additions, etc., or depending on the functional load. In other words, the above-mentioned embodiments may be implemented in any combination, or the embodiments may be implemented selectively.

1 糖度計測装置
10 分光検出装置
11 計測部
12 光源
20 データ処理装置
21 メモリ
22 キーボード・マウス
23 制御部
24 計算処理部
24-1 分光吸光スペクトル取得部
24-2 推定モデル作成部
24-3 推定部
26 I/Oポート
27 推定モデル
30 ディスプレイ
41 本体部
42 試料台
43 受光ファイバー
44 カバー
51 本体部
52 試料台
53 光入射口
54 超小型分光器
REFERENCE SIGNS LIST 1 Sugar content measuring device 10 Spectroscopic detection device 11 Measuring unit 12 Light source 20 Data processing device 21 Memory 22 Keyboard/mouse 23 Control unit 24 Calculation processing unit 24-1 Spectroscopic absorption spectrum acquisition unit 24-2 Estimation model creation unit 24-3 Estimation unit 26 I/O port 27 Estimation model 30 Display 41 Main unit 42 Sample stage 43 Light receiving fiber 44 Cover 51 Main unit 52 Sample stage 53 Light entrance 54 Ultra-compact spectrometer

Claims (21)

光源からの光源光を青果物に照射し、その拡散反射を含む反射光を受光して当該青果物の糖度を計測する糖度計測方法であって、
1つの光源又は複数の光源から照射される近赤外線波長域の光を青果物に照射し、その拡散反射を含む反射光から波長の異なる4つのみの分光吸光スペクトルの光のシグナル値を取得する取得工程と、
取得した波長の異なる4つのみの吸光度(光のシグナル値)を説明変数として使用して、多変量解析を行うことで、糖度を推定するための推定モデルを作成する推定モデル作成工程と、
を含むことを特徴とする糖度計測方法。
A method for measuring sugar content of fruits or vegetables by irradiating the fruits or vegetables with light from a light source and receiving reflected light including diffuse reflection, comprising:
an acquisition step of irradiating fruit or vegetables with light in the near-infrared wavelength range from one light source or multiple light sources, and acquiring signal values of light having only four different wavelengths of spectral absorption spectra from the reflected light, including diffuse reflection;
A process of creating an estimation model for estimating sugar content by performing multivariate analysis using only the four absorbance values (light signal values) obtained at different wavelengths as explanatory variables;
A method for measuring sugar content comprising the steps of:
前記異なる4つの波長は、856±2nm、876±2nm又は884±2nm、900~918nm、926±2nmであることを特徴とする請求項1に記載の糖度計測方法。 The sugar content measurement method according to claim 1, characterized in that the four different wavelengths are 856±2 nm, 876±2 nm or 884±2 nm, 900-918 nm, and 926±2 nm. 前記推定モデル作成工程では、前記多変量解析として重回帰分析を行って、前記推定モデルとして重回帰式を作成することを特徴とする請求項1又は2に記載の糖度計測方法。 The sugar content measurement method according to claim 1 or 2, characterized in that in the estimation model creation step, a multiple regression analysis is performed as the multivariate analysis to create a multiple regression equation as the estimation model. 前記青果物は、当該青果物の赤道部、花痕部、加工品、及び破砕物を含むことを特徴とする請求項1~3のいずれか1つに記載の糖度計測方法。 The method for measuring sugar content according to any one of claims 1 to 3, characterized in that the fruits and vegetables include the equator, flower scar, processed products, and crushed products of the fruits and vegetables. 前記青果物は、西洋ナシ、カキ、マンゴー、イチゴ、パプリカ、ミカン(Tangerineを含む)、不知火、トマト、サクランボ、モモ、ナシ、リンゴ、スモモ、及びメロンを含むことを特徴とする請求項1~請求項4のいずれか1つに記載の糖度計測方法。 The method for measuring sugar content according to any one of claims 1 to 4, characterized in that the fruits and vegetables include pears, persimmons, mangoes, strawberries, paprikas, mandarin oranges (including tangerines), shiranui, tomatoes, cherries, peaches, pears, apples, plums, and melons. 前記光源は、リング状の光源であることを特徴とする請求項1~5のいずれか1つに記載の糖度計測方法。 The sugar content measurement method according to any one of claims 1 to 5, characterized in that the light source is a ring-shaped light source. 前記取得工程では、前記リング状の光源から照射される光に対する青果物の拡散反射を含む反射光を、当該リング状の光源の略中心で検出することを特徴とする請求項6に記載の糖度計測方法。 The sugar content measurement method according to claim 6, characterized in that in the acquisition process, reflected light, including diffuse reflection from fruits and vegetables, of the light irradiated from the ring-shaped light source is detected at approximately the center of the ring-shaped light source. 前記1つの光源は、ハロゲンランプであることを特徴とする請求項6又は7に記載の糖度計測方法。 The method for measuring sugar content according to claim 6 or 7, characterized in that the one light source is a halogen lamp. 前記複数の光源は、発光波長の異なる複数のLEDであることを特徴とする請求項6又は7に記載の糖度計測方法。 The method for measuring sugar content according to claim 6 or 7, characterized in that the multiple light sources are multiple LEDs with different emission wavelengths. 前記複数のLEDは、計測対象の青果物側に平面に対して所定角度(但し、所定角度は10度より大きい)を有して配置されることを特徴とする請求項9に記載の糖度計測方法。 The method for measuring sugar content according to claim 9, characterized in that the LEDs are arranged at a predetermined angle (however, the predetermined angle is greater than 10 degrees) with respect to a plane on the side of the fruit or vegetable to be measured. 光源からの光源光を青果物に照射し、その拡散反射を含む反射光を受光して当該青果物の糖度を計測する糖度計測装置であって、
1つの光源又は複数の光源から照射される近赤外線波長域の光を青果物に照射し、その拡散反射を含む反射光から波長の異なる4つのみの分光吸光スペクトルの光のシグナル値を取得する分光検出手段と、
取得した波長の異なる4つのみの吸光度(光のシグナル値)を説明変数として使用して、多変量解析を行うことで、糖度含有量を推定するための推定モデルを作成する推定モデル作成手段と、
を備えたことを特徴とする糖度計測装置。
A sugar content measuring device that irradiates a fruit or vegetable with light from a light source and receives reflected light including diffuse reflection to measure the sugar content of the fruit or vegetable,
a spectroscopic detection means for irradiating fruit or vegetables with light in the near-infrared wavelength range from one light source or a plurality of light sources, and acquiring signal values of light having only four different wavelengths of spectral absorption spectra from the reflected light including diffuse reflection;
an estimation model creation means for creating an estimation model for estimating sugar content by performing multivariate analysis using only the four absorbances (light signal values) obtained at different wavelengths as explanatory variables;
A sugar content measuring device comprising:
前記異なる4つの波長は、856±2nm、876±2nm又は884±2nm、900~918nm、926±2nmであることを特徴とする請求項11に記載の糖度計測装置。 The sugar content measuring device according to claim 11, characterized in that the four different wavelengths are 856±2 nm, 876±2 nm or 884±2 nm, 900 to 918 nm, and 926±2 nm. 前記推定モデル作成手段は、前記多変量解析として重回帰分析を行って、前記推定モデルとして重回帰式を作成することを特徴とする請求項11又は12に記載の糖度計測装置。 The sugar content measuring device according to claim 11 or 12, characterized in that the estimation model creation means performs multiple regression analysis as the multivariate analysis and creates a multiple regression equation as the estimation model. 前記青果物は、当該青果物の赤道部、花痕部、加工品、及び破砕物を含むことを特徴とする請求項11~13のいずれか1つに記載の糖度計測装置。 The sugar content measuring device according to any one of claims 11 to 13, characterized in that the fruit or vegetable includes the equator, flower scar, processed products, and crushed products of the fruit or vegetable. 前記青果物は、西洋ナシ、カキ、マンゴー、イチゴ、パプリカ、ミカン(Tangerineを含む)、不知火、トマト、サクランボ、モモ、ナシ、リンゴ、スモモ、及びメロンを含むことを特徴とする請求項11~請求項14のいずれか1つに記載の糖度計測装置。 The sugar content measuring device according to any one of claims 11 to 14, characterized in that the fruits and vegetables include pears, persimmons, mangoes, strawberries, paprika, mandarin oranges (including tangerines), shiranui, tomatoes, cherries, peaches, pears, apples, plums, and melons. 前記光源は、リング状の光源であることを特徴とする請求項11~15のいずれか1つに記載の糖度計測装置。 The sugar content measuring device according to any one of claims 11 to 15, characterized in that the light source is a ring-shaped light source. 前記分光検出手段は、前記リング状の光源から照射される光に対する青果物の拡散反射を含む反射光を、当該リング状の光源の略中心で検出することを特徴とする請求項16に記載の糖度計測装置。 The sugar content measuring device according to claim 16, characterized in that the spectroscopic detection means detects reflected light, including diffuse reflection from fruits and vegetables, of the light irradiated from the ring-shaped light source at approximately the center of the ring-shaped light source. 前記1つの光源は、ハロゲンランプであることを特徴とする請求項16又は17に記載の糖度計測装置。 The sugar content measuring device according to claim 16 or 17, characterized in that the one light source is a halogen lamp. 前記複数の光源は、発光波長の異なる複数のLEDであることを特徴とする請求項16又は17に記載の糖度計測装置。 The sugar content measuring device according to claim 16 or 17, characterized in that the multiple light sources are multiple LEDs with different emission wavelengths. 前記複数のLEDは、計測対象の青果物側に平面に対して所定角度(但し、所定角度は10度より大きい)を有して配置されることを特徴とする請求項19に記載の糖度計測装置。 The sugar content measuring device according to claim 19, characterized in that the LEDs are arranged at a predetermined angle (however, the predetermined angle is greater than 10 degrees) with respect to a plane on the side of the fruit or vegetable to be measured. 青果物の糖度を計測するための糖度計測プログラムであって、
1つの光源又は複数の光源から照射される近赤外線波長域の光に対する、青果物からの拡散反射を含む反射光から波長の異なる4つのみの分光吸光スペクトルの光のシグナル値を取得する取得工程と、
取得した波長の異なる4つのみの吸光度(光のシグナル値)を説明変数として使用して、多変量解析を行うことで、糖度を推定するための推定モデルを作成する推定モデル作成工程と、
コンピュータに実行させるための糖度計測プログラム。
A sugar content measurement program for measuring the sugar content of fruits and vegetables,
an acquisition step of acquiring signal values of light having only four different wavelengths of spectral absorption spectra from reflected light, including diffuse reflection, from fruits and vegetables in response to light in the near-infrared wavelength range irradiated from one light source or a plurality of light sources;
A process of creating an estimation model for estimating sugar content by performing multivariate analysis using only the four absorbance values (light signal values) obtained at different wavelengths as explanatory variables;
A sugar content measurement program to be executed by a computer.
JP2021204633A 2020-12-16 2021-12-16 Sugar content measurement method, sugar content measurement device, and sugar content measurement program Active JP7688907B2 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2020208786 2020-12-16
JP2020208786 2020-12-16

Publications (2)

Publication Number Publication Date
JP2022095598A JP2022095598A (en) 2022-06-28
JP7688907B2 true JP7688907B2 (en) 2025-06-05

Family

ID=82163073

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2021204633A Active JP7688907B2 (en) 2020-12-16 2021-12-16 Sugar content measurement method, sugar content measurement device, and sugar content measurement program

Country Status (1)

Country Link
JP (1) JP7688907B2 (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116148210A (en) * 2023-02-10 2023-05-23 北京市农林科学院信息技术研究中心 Fruit quality detection method and device
JP2025087961A (en) * 2023-11-30 2025-06-11 国立研究開発法人農業・食品産業技術総合研究機構 Non-destructive measurement device fixture
JP2025099256A (en) * 2023-12-21 2025-07-03 株式会社クボタ Fruit and vegetable quality measurement system, calibration curve switching method, and computer program

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006226775A (en) 2005-02-16 2006-08-31 Toyohashi Univ Of Technology Fruit taste component evaluation method and evaluation apparatus
JP2007300105A (en) 2006-05-03 2007-11-15 Samsung Electronics Co Ltd Substrate inspection apparatus and its illumination unit
JP2008209726A (en) 2007-02-27 2008-09-11 Olympus Corp Illuminating device
JP2015108508A (en) 2012-03-14 2015-06-11 千代田電子工業株式会社 Non-destructive measurement apparatus
WO2018047366A1 (en) 2016-09-06 2018-03-15 株式会社アタゴ Nondestructive measuring device
JP2019211456A (en) 2018-06-07 2019-12-12 学校法人 東洋大学 Quality inspection device for fruits, vegetables, and fruit juice

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3346142B2 (en) * 1994-12-28 2002-11-18 住友金属鉱山株式会社 Non-destructive sugar content measuring device

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006226775A (en) 2005-02-16 2006-08-31 Toyohashi Univ Of Technology Fruit taste component evaluation method and evaluation apparatus
JP2007300105A (en) 2006-05-03 2007-11-15 Samsung Electronics Co Ltd Substrate inspection apparatus and its illumination unit
JP2008209726A (en) 2007-02-27 2008-09-11 Olympus Corp Illuminating device
JP2015108508A (en) 2012-03-14 2015-06-11 千代田電子工業株式会社 Non-destructive measurement apparatus
WO2018047366A1 (en) 2016-09-06 2018-03-15 株式会社アタゴ Nondestructive measuring device
JP2019211456A (en) 2018-06-07 2019-12-12 学校法人 東洋大学 Quality inspection device for fruits, vegetables, and fruit juice

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
伊藤秀和、森本進、堀江秀樹,近赤外分光法によるメロン糖度の非破壊計測法の開発,第48回自動制御連合講演会講演論文集,日本,自動制御連合講演会,2005年11月26日,p.981-985,https://www.jstage.jst.go.jp/article/jacc/48/0/48_0_282/_pdf/-char/ja

Also Published As

Publication number Publication date
JP2022095598A (en) 2022-06-28

Similar Documents

Publication Publication Date Title
Fan et al. Non-destructive evaluation of soluble solids content of apples using a developed portable Vis/NIR device
Cortés et al. Monitoring strategies for quality control of agricultural products using visible and near-infrared spectroscopy: A review
Cayuela Vis/NIR soluble solids prediction in intact oranges (Citrus sinensis L.) cv. Valencia Late by reflectance
JP7688907B2 (en) Sugar content measurement method, sugar content measurement device, and sugar content measurement program
Wang et al. Vis/NIR optical biosensors applications for fruit monitoring
Liu et al. Use of FT-NIR spectrometry in non-invasive measurements of internal quality of ‘Fuji’apples
Xie et al. Quantification of glucose, fructose and sucrose in bayberry juice by NIR and PLS
Magwaza et al. Evaluation of Fourier transform-NIR spectroscopy for integrated external and internal quality assessment of Valencia oranges
Yu et al. Nondestructive determination of SSC in Korla fragrant pear using a portable near-infrared spectroscopy system
Ignat et al. Forecast of apple internal quality indices at harvest and during storage by VIS-NIR spectroscopy
US20020011567A1 (en) Apparatus and method and techniques for measuring and correlating characteristics of fruit with visible/near infra-red spectrum
Nturambirwe et al. Non-destructive measurement of internal quality of apple fruit by a contactless NIR spectrometer with genetic algorithm model optimization
Luo et al. Robustness improvement of NIR-based determination of soluble solids in apple fruit by local calibration
WO2020196690A1 (en) Product inspection method and product inspection device
Zude Non-destructive prediction of banana fruit quality using VIS/NIR spectroscopy
Huishan et al. Application Fourier transform near infrared spectrometer in rapid estimation of soluble solids content of intact citrus fruits
JP7169643B2 (en) Nitrate ion concentration nondestructive measuring method, nitrate ion concentration nondestructive measuring device, and nitrate ion concentration nondestructive measuring program
Han et al. Evaluation of the optical layout and sample size on online detection of apple watercore and SSC using Vis/NIR spectroscopy
Xu et al. Nondestructive detection of internal flavor in ‘Shatian’pomelo fruit based on visible/near infrared spectroscopy
Nturambirwe et al. Detecting bruise damage and level of severity in apples using a contactless nir spectrometer
JP2025061815A (en) Functional component measuring method, functional component measuring device, and functional component measuring program
Sun et al. Non-destructive detection of blackheart and soluble solids content of intact pear by online NIR spectroscopy: X. Sun et al.
PUANGSOMBUT et al. Evaluation of internal quality of fresh‐cut pomelo using vis/nir transmittance
Wang et al. Nondestructive Quality Detection of Characteristic Fruits Based on Vis/NIR Spectroscopy: Principles, Systems, and Applications.
Sharabiani et al. Non-destructive assessment of quality parameters in Javadi cv. Peach fruits using Vis/NIR spectroscopy and multiple regression analysis

Legal Events

Date Code Title Description
AA64 Notification of invalidation of claim of internal priority (with term)

Free format text: JAPANESE INTERMEDIATE CODE: A241764

Effective date: 20220201

RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20220217

RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7422

Effective date: 20220222

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20240419

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20250226

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20250304

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20250404

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20250519

R150 Certificate of patent or registration of utility model

Ref document number: 7688907

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150