(57)【要約】      （修正有）【課題】アルカリ溶液による米の溶解過程を、精度よく、且つ、定量的に評価可能なアルカリ崩壊試験方法を提供する。【解決手段】米をアルカリ溶液に浸漬する工程Ｓ１２、経時的に撮像することにより、複数の画像を取得する工程Ｓ１３、取得した複数の画像の各々について輝度分布を求める工程Ｓ１４、輝度分布の時間的変化に基づく特徴量のプロファイルを求める工程Ｓ１５、プロファイルに基づいて、輝度分布に基づいて求められる評価パラメータが基準値に達するまでに要する特性時間を求める工程Ｓ１７を含む、アルカリ崩壊試験方法。【選択図】図１
2020115122.tif 000002
【特許請求の範囲】
【請求項１】
米をアルカリ溶液に浸す工程と、
前記アルカリ溶液に浸された前記米を経時的に撮像することにより、複数の画像を取得する工程と、
前記複数の画像の各々について輝度分布を求める工程と、
前記輝度分布の時間的変化に基づく特徴量のプロファイルを求める工程と、
前記プロファイルに基づいて、前記輝度分布に基づいて求められる評価パラメータが基準値に達するまでに要する特性時間を求める工程と、
を備える、アルカリ崩壊試験方法。
【請求項２】
前記輝度分布は、前記画像をアナログ処理することにより前記画像に含まれるピクセル毎の輝度値を算出することによって求められる、請求項１に記載のアルカリ崩壊試験方法。
【請求項３】
前記プロファイルにおける前記評価パラメータのピーク値を求める工程を更に備え、
前記特性時間は、前記ピーク値を用いて正規化された前記プロファイルに基づいて求められる、請求項１又は２に記載のアルカリ崩壊試験方法。
【請求項４】
前記基準値は５０％である、請求項３に記載のアルカリ崩壊試験方法。
【請求項５】
前記評価パラメータは、前記輝度分布の輝度分散値である、請求項１から４のいずれか一項に記載のアルカリ崩壊試験方法。
【請求項６】
前記米は、一粒ずつ独立した容器内において前記アルカリ溶液に浸される、請求項１から５のいずれか一項に記載のアルカリ崩壊試験方法。
【請求項７】
同じ品種の前記米を複数用意し、前記米の各々について求めた前記特性時間を統計処理することにより前記品種の評価を行う、請求項１から６のいずれか一項に記載のアルカリ崩壊試験方法。
【請求項８】
異なる条件に対応する前記米を用意し、前記米の各々について求めた前記特性時間を比較することにより前記異なる条件間における相対的評価を行う、請求項１から７のいずれか一項に記載のアルカリ崩壊試験方法。
【請求項９】
前記アルカリ溶液は、３％以上のアルカリ濃度を有する、請求項１から８のいずれか一項に記載のアルカリ崩壊試験方法。
【請求項１０】
前記特性時間に基づいてＤＳＣを評価する工程を更に備える、請求項１から９のいずれか一項に記載のアルカリ崩壊試験方法。
【請求項１１】
前記プロファイルに基づいて前記米における浸漬割れの有無を判定する工程を更に備える、請求項１から１０のいずれか一項に記載のアルカリ崩壊試験方法。
【請求項１２】
前記浸漬割れを判定する工程では、前記特性時間が予め設定された閾値以下である場合に、前記米に前記浸漬割れが有ると判定する、請求項１１に記載のアルカリ崩壊試験方法。
【請求項１３】
前記浸漬割れを判定する工程では、前記プロファイルの時間微分関数のピーク値に到達する時間が閾値以下である場合に、前記米に前記浸漬割れが有ると判定する、請求項１１に記載のアルカリ崩壊試験方法。
【請求項１４】
前記浸漬割れを判定する工程では、前記プロファイルの時間微分関数のピーク値が閾値以上である場合に、前記米に前記浸漬割れが有ると判定する、請求項１１に記載のアルカリ崩壊試験方法。
【請求項１５】
複数の前記米について前記浸漬割れを判定する工程を実施することで、前記浸漬割れの発生率を求める工程を更に備える、請求項１１から１４のいずれか一項に記載のアルカリ崩壊試験方法。
【請求項１６】
前記特性時間を求める工程では、前記浸漬割れを判定する工程において前記浸漬割れが有ると判定された前記米を除く前記米について求められた前記プロファイルを用いて前記特性時間を求める、請求項１１から１５のいずれか一項に記載のアルカリ崩壊試験方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は、アルカリ溶液を用いて米等の溶解特性を評価するためのアルカリ崩壊試験方法に関する。
【背景技術】
【０００２】
米の品質評価手法の一つとして、アルカリ崩壊試験が知られている。アルカリ崩壊試験では、米のアルカリ溶液に対する抵抗性が、品種、産地、生産年等の様々な条件によって変化する特性を利用し、アルカリ溶液に浸漬された米の崩壊性に基づいて品質評価が行われる。典型的なアルカリ崩壊試験では、アルカリ溶液に対して数十粒の米を浸漬させ、数時間放置する。その過程において、米が溶解することでアルカリ溶液に白濁が生じる様子を目視観察することによって、崩壊性が評価される。このようなアルカリ崩壊試験によって評価される崩壊性は、食用米では炊飯特性や食味などの嗜好性に関する質的特性の指標とされ、酒用米では醸造条件と密接に関連することが知られている。
【０００３】
従来のアルカリ崩壊試験では、上述のように目視観察に基づいた評価であるため、作業者の経験的要素や主観的要素が少なからず含まれることから安定性に課題があった。このような課題を解決するために、例えば特許文献１では、異なる濃度のアルカリ溶液を用意し、各アルカリ溶液において米が溶解するまでに要した時間に基づいて崩壊性を評価することで、酒造適正の評価精度の向上を図っている。また非特許文献１では、アルカリ溶液中に浸漬された米を撮像して取得した画像をデジタル処理することにより、アルカリ溶液における白濁の広がりを定量的に評価することが提案されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１７−１６１３９９号公報
【非特許文献】
【０００５】
【非特許文献１】Ｊｏｒｎａｌ  ｏｆ  Ｃｅｒｅａｌ  Ｃｉｅｎｃｅ  ５２（２０１０）２２７−２３５、Ａｌｋａｌｉ  ｓｐｒｅｄｉｎｇ  ｖａｌｕｅ  ａｎｄ  Ｉｍａｇｅ  Ａｎａｌｙｓｉｓ，Ｍ．Ｍａｒｉｏｔｔｉ、Ｌ．Ｆｏｎｇａｒｏ，Ｆ．Ｃａｔｅｎａｃｃｉ
【発明の概要】
【発明が解決しようとする課題】
【０００６】
上記特許文献１では、様々な濃度のアルカリ溶液に米を浸漬した際に、米が溶解するまでに要する時間に基づいて評価を行っているが、米が溶解したか否かの判断は依然として目視観察によって行われている。そのため、従来のアルカリ崩壊試験の安定性に関する課題が残っている。
【０００７】
また上記非特許文献１では、アルカリ溶液に白濁が生じている範囲を、画像中のピクセル毎に輝度が閾値以上であるか否かを判定することでデジタル的に特定している。ここで閾値が高く設定された場合には、白濁が濃い領域は特定できるものの、白濁が薄い領域の特定は難しい。逆に、閾値が低く設定された場合には、白濁が薄い領域は特定できるものの、白濁が濃い領域の特定は難しい。このように非特許文献１の技術では、アルカリ溶液に米を浸漬した際に、白濁範囲が次第に広がり、濃淡が変化する一連の様子を有効に評価することが難しい。
【０００８】
本発明の少なくとも一実施形態は上述の事情に鑑みなされたものであり、アルカリ溶液による米の溶解過程を、精度よく、且つ、定量的に評価可能なアルカリ崩壊試験方法を提供することを目的とする。
【課題を解決するための手段】
【０００９】
（１）本発明の少なくとも一実施形態に係るアルカリ崩壊試験方法は上記課題を解決するために、
米をアルカリ溶液に浸す工程と、
前記アルカリ溶液に浸された前記米を経時的に撮像することにより、複数の画像を取得する工程と、
前記複数の画像の各々について輝度分布を求める工程と、
前記輝度分布の時間的変化に基づく特徴量のプロファイルを求める工程と、
前記プロファイルに基づいて、前記輝度分布に基づいて求められる評価パラメータが基準値に達するまでに要する特性時間を求める工程と、
を備える。
【００１０】
上記（１）の方法によれば、アルカリ溶液に浸漬された米が溶解する様子を撮像することで取得した画像を解析することにより、崩壊性の定量的な評価が可能となる。特に、各画像に含まれる輝度分布の時間変化に基づく特徴量のプロファイルに基づいて特性時間を求めることで、米が溶解するに従ってアルカリ溶液中で白濁が広がりながら濃淡が変化する時間的変化を適切に評価できる。
【００１１】
（２）幾つかの実施形態では上記（１）の方法において、
前記輝度分布は、前記画像をアナログ処理することにより前記画像に含まれるピクセル毎の輝度値を算出することによって求められる。
【００１２】
上記（２）の方法によれば、各画像の輝度分布をアナログ処理することで算出されたプロファイルに基づいた評価が可能となる。これにより、米が溶解するに従ってアルカリ溶液中で白濁が広がりながら濃淡が変化する時間的変化を適切に評価できる。
【００１３】
（３）幾つかの実施形態では上記（１）又は（２）の方法において、
前記プロファイルにおける前記評価パラメータのピーク値を求める工程を更に備え、
前記特性時間は、前記ピーク値を用いて正規化された前記プロファイルに基づいて求められる。
【００１４】
上記（３）の方法によれば、各プロファイルをピーク値で正規化することで、異なる条件（例えば、異なる品種、産地、生産年等）に対応する米同士を、定量的な評価パラメータのもとで、同等に評価できる。
【００１５】
（４）幾つかの実施形態では上記（３）の方法において、
前記基準値は５０％である。
【００１６】
上記（４）の方法によれば、正規化されたプロファイルに基づいて評価パラメータを求める際に用いられる基準値として、各条件（品種、産地、生産年等）の変化が最も反映されやすい５０％が採用される。これにより、ピーク値で正規化することで、異なる条件（例えば、異なる品種、産地、生産年等）に対応する米同士を効果的に評価できる。
【００１７】
（５）幾つかの実施形態では上記（１）から（４）のいずれか一方法において、
前記評価パラメータは、前記輝度分布の輝度分散値である。
【００１８】
上記（５）の方法によれば、評価パラメータとして輝度分布の輝度分散値を選定することで、好適な評価が可能となる。 
【００１９】
（６）幾つかの実施形態では上記（１）から（５）のいずれか一方法において、
前記米は、一粒ずつ独立した容器内において前記アルカリ溶液に浸される。
【００２０】
上記（６）の方法によれば、複数粒の米に対して本方法を実施する際には、独立した容器内に一粒ずつ米を浸すことで、粒間の相互作用の影響を排除した精度のよい評価が可能となる。 
【００２１】
（７）幾つかの実施形態では上記（１）から（６）のいずれか一方法において、
同じ品種の前記米を複数用意し、前記米の各々について求めた前記特性時間を統計処理することにより前記品種の評価を行う。
【００２２】
上記（７）の方法によれば、複数粒の米について本方法による評価を行う場合には、得られた特性時間を統計処理することで、より信頼性の高い評価が可能となる。 
【００２３】
（８）幾つかの実施形態では上記（１）から（７）のいずれか一方法において、
異なる条件に対応する前記米を用意し、前記米の各々について求めた前記特性時間を比較することにより前記異なる条件間における相対的評価を行う。
【００２４】
上記（８）の方法によれば、条件の異なる米に対してそれぞれ特性時間を求め、得られた各特性時間を比較する。これにより、異なる条件に対応する米同士を定量的に評価できる。
【００２５】
（９）幾つかの実施形態では上記（１）から（８）のいずれか一方法において、
前記アルカリ溶液は、３％以上のアルカリ濃度を有する。
【００２６】
上記（９）の方法によれば、３％以上（より好ましくは５％以上）のアルカリ濃度を有するアルカリ溶液を用いて試験を実施することで、従来のアルカリ崩壊試験に比べて短い試験時間で評価を行うことができる。上述のように本方法では評価パラメータによる定量的評価が可能であるため、このように濃度が高いアルカリ溶液を用いた場合でも、精度のよい評価が可能となる（従来のアルカリ崩壊試験では目視観察による評価であったため、濃度が高いアルカリ溶液を用いると評価精度が低下してしまう）。
【００２７】
（１０）幾つかの実施形態では上記（１）から（９）のいずれか一方法において、
前記特性時間に基づいてＤＳＣを推定する工程を更に備える。
【００２８】
本発明者の研究によれば、本方法で求められる特性時間はＤＳＣと相関を有することが見出された。そのため、本方法で特性時間を求め、予め把握された相関に基づいてＤＳＣを定量的に推定することが可能となる。
【００２９】
（１１）幾つかの実施形態では上記（１）から（１０）のいずれか一方法において、
前記プロファイルに基づいて前記米における浸漬割れの有無を判定する工程を更に備える。
【００３０】
上記（１１）の方法によれば、各画像に含まれる輝度分布の時間変化に基づく特徴量であるプロファイルに基づいて、アルカリ溶液における米の溶けやすさを定量的に評価することで、米の内部に存在するクラックのような潜在的要素が発展して生じる浸漬割れの有無を適切に判定することができる。
【００３１】
（１２）幾つかの実施形態では上記（１１）の方法において、
前記浸漬割れを判定する工程では、前記特性時間が予め設定された閾値以下である場合に、前記米に前記浸漬割れが有ると判定する。
【００３２】
上記（１２）の方法によれば、浸漬割れが生じる米は溶解の進行度が早いという特性に鑑み、特性時間が閾値以下であるか否かに基づいて浸漬割れの有無を適切に判定できる。
【００３３】
（１３）幾つかの実施形態では上記（１１）の方法において、
前記浸漬割れを判定する工程では、前記プロファイルの時間微分関数のピーク値に到達する時間が閾値以下である場合に、前記米に前記浸漬割れが有ると判定する。
【００３４】
上記（１３）の方法によれば、浸漬割れが生じる米はプロファイルの形状に特異な振る舞いを示すという特性に鑑み、プロファイルの時間微分関数のピーク値に到達する時間が閾値以下であるか否かに基づいて浸漬割れの有無を適切に判定できる。
【００３５】
（１４）幾つかの実施形態では上記（１１）の方法において、
前記浸漬割れを判定する工程では、前記プロファイルの時間微分関数のピーク値が閾値以上である場合に、前記米に前記浸漬割れが有ると判定する。
【００３６】
上記（１４）の方法によれば、浸漬割れが生じる米はプロファイルの形状に特異な振る舞いを示すという特性に鑑み、プロファイルの時間微分関数のピーク値が閾値以上であるか否かに基づいて浸漬割れの有無を適切に判定できる。
【００３７】
（１５）幾つかの実施形態では上記（１１）から（１４）のいずれか一方法において、
複数の前記米について前記浸漬割れを判定する工程を実施することで、前記浸漬割れの発生率を求める工程を更に備える。
【００３８】
上記（１５）の方法によれば、複数の米について浸漬割れの有無を判定することで、浸漬割れの発生率による米の品質評価を行うことができる。
【００３９】
（１６）幾つかの実施形態では上記（１１）から（１５）のいずれか一方法において、
前記特性時間を求める工程では、前記浸漬割れを判定する工程において前記浸漬割れが有ると判定された前記米を除く前記米について求められた前記プロファイルを用いて前記特性時間を求める。
【００４０】
上記（１６）の方法によれば、浸漬割れの影響を排除した米の品質評価が可能となる。
【発明の効果】
【００４１】
本発明の少なくとも一実施形態によれば、アルカリ溶液による米の溶解過程を、精度よく、且つ、定量的に評価可能なアルカリ崩壊試験方法を提供できる。
【図面の簡単な説明】
【００４２】
【図１】本発明の少なくとも一実施形態に係るアルカリ崩壊試験方法を工程毎に示すフローチャートである。
【図２】複数粒の米が互いに独立した容器内でそれぞれアルカリ溶液に浸漬される様子を示す模式図である。
【図３Ａ】図２の一つのセルでアルカリ溶液１２に浸漬される一粒の米の様子を撮像した画像の一例である。
【図３Ｂ】図２の一つのセルでアルカリ溶液１２に浸漬される一粒の米の様子を撮像した画像の一例である。
【図４】図１のステップＳ１５で求められるプロファイルの一例である。
【図５】図４のプロファイルの正規化プロファイルである。
【図６】幾つかの品種について得られた正規化プロファイルを比較して示す試験結果の一例である。
【図７】図６で求められた各品種の特性時間を比較した結果である。
【図８】図６で求められた特性時間とＢｒｉｘ値とを品種毎に比較した結果である。
【図９】互いに異なる３品種について求められた特性時間とＤＣＳ値（熱伝導性）とを品種毎に比較した結果である。
【図１０】ＫＯＨ濃度がそれぞれ２％、１０％のアルカリ溶液を用いた場合に得られた正規化プロファイルの一例である。
【図１１】アルカリ崩壊試験を利用した浸漬割れの判定方法を工程毎に示すフローチャートである。
【図１２】浸漬割れが生じる前後における撮像イメージである。
【図１３】図１２の各サンプルについて得られたプロファイルである。
【図１４】他のサンプルについて取得されたプロファイル（規格化済）に基づいて判定例である。
【図１５】図１４のプロファイルの時間微分関数のピーク値に到達する時間を用いた判定例である。
【図１６】図１４のプロファイルの時間微分関数のピーク値を用いた判定例である。
【発明を実施するための形態】
【００４３】
以下、添付図面に従って本発明の実施形態について説明する。ただし、この実施形態に記載されている構成部品の寸法、材質、形状、その相対的配置等は、特定的な記載がない限り本発明の範囲をこれに限定する趣旨ではなく、単なる説明例にすぎない。
【００４４】
本発明の少なくとも一実施形態に係るアルカリ崩壊試験方法は、アルカリ溶液に浸漬された米の崩壊性（溶解性）を評価するための試験である。図１は本発明の少なくとも一実施形態に係るアルカリ崩壊試験方法を工程毎に示すフローチャートである。
【００４５】
まず試験対象となる米を用意する（ステップＳ１１）。米は食用米であってもよく、醸造米であってもよく、予め品種、生産地、生産年、生産者等の関連する条件が特定される。また用意する粒数は一粒であってもよいが、本実施形態では、個体差による影響を排除するために複数粒を用意し、後述するように各粒に対する評価結果を統計的に処理することで、評価精度を向上させている。
【００４６】
続いてステップＳ１１で用意した米を、試験溶液であるアルカリ溶液に浸漬する（ステップＳ１２）。アルカリ溶液は所定の容器内に用意され、これに試験対象となる米を投入することで浸漬される。この容器は、後述するように内部でアルカリ溶液中に浸漬される米の様子を撮像するために適した形状を有することが好ましい。本実施形態では、撮像が行われる上方側に向けて開口を有するとともに、内部に投入された米の姿勢が安定するように平坦な底部を有する有底容器が用いられる。このような容器の一例としては、透明なプラスチック材料やガラス材料から形成されたシャーレや試験管がある。
【００４７】
試験に用いられるアルカリ溶液は、例えば強アルカリ成分が好ましく、例えば、ＫＯＨ、ＮａＯＨ、ＬｉＯＨ、ＲｂＯＨ、ＣｓＯＨ、Ｃａ（ＯＨ）２、Ｓｒ（ＯＨ）２等が用いられる。アルカリ溶液のアルカリ濃度は、例えば１．７重量％の範囲が好ましい。
【００４８】
ここでステップＳ１１において複数粒の米を用意した場合には、各粒が互いに独立した容器内でアルカリ溶液に浸漬されることが好ましい。図２は複数粒の米１が互いに独立した容器内でそれぞれアルカリ溶液（透明なため不図示）に浸漬される様子を示す模式図である。この例では、互いに独立した容器として、マトリクス状に隣接配置された複数のセル１０が用いられており、各セル１０にはアルカリ溶液１２が貯留される。そして、各セル１０には一粒ずつ米１が投入されることで、アルカリ溶液１２に浸漬される。仮に複数粒の米１を同一の容器内でアルカリ溶液１２に浸漬すると、粒間の相互作用によって溶解過程に影響を及ぼすおそれがある。その点、本実施形態では、各米１を独立した容器で浸漬することで、２粒間の相互作用の影響を排除し、精度の高い評価が可能となる。
【００４９】
続いてアルカリ溶液１２に浸された米１を撮像することにより、画像を取得する（ステップＳ１３）。ステップＳ１３における米の撮像は、経時的に繰り返し行われることにより、米１がアルカリ溶液１２中で溶解する様子が、複数の画像にわたって撮影される。すなわち、これら複数の画像には、時間の経過に伴ってアルカリ溶液１２中で溶解する米１の各瞬間における様子が記録される。
【００５０】
また図２のように各セル１０に米１を一粒ずつ配置した場合には、ステップＳ１３における撮像もセル１０毎に実施される。図３Ａ及び図３Ｂは図２の一つのセル１０でアルカリ溶液１２に浸漬される一粒の米の様子を撮像した画像の一例である。この例では、安定した撮像を行うために、撮像装置としてイメージスキャナを用いるとともに、黒紙を背景とした暗所で試験及び撮影を実施することで外乱光の影響を排除している。図３Ａはアルカリ溶液１２に浸漬された直後（すなわち溶解の初期段階）における米１を撮像した画像であり、溶解が進行していないことから、本来の米１の輪郭が明確に残存している様子が映し出されている。一方の図３Ｂはアルカリ溶液１２に浸漬されてからしばらく時間が経過した段階における米１を撮像した画像であり、米１が溶解することで生じた白濁範囲１４が周囲に広がっている様子が映し出されている。図３Ａ及び図３Ｂでは、代表的に２つの時間における画像を示しているが、このように異なる時間にわたって複数の画像を取得することで、一粒の米１が時間とともに溶解していく様子が記録される。
【００５１】
尚、アルカリ溶液１２に米１を浸漬しながら撮像を行う際には、安定的な測定環境が維持されることが好ましい。例えば、周囲温度は米１の溶けやすさに影響を与えるため（例えば、周囲温度が高いほど米１の溶解速度が速くなるため）、試験が実施されている最中は周囲温度が安定するように管理するとよい。ただし、周囲温度が変化する場合には、当該周囲温度の変化による影響を予め考慮しておき、試験結果に対して適宜補正処理を行うことで最終的な試験結果に反映するようにしてもよい。
【００５２】
ステップＳ１３では、このように各セル１０において米１が溶解する過程を経時的に撮像することにより、複数の画像が取得される。このように取得された複数の画像は、撮像時間、セル番号等によって個別に特定可能な態様で管理される。
【００５３】
続いてステップＳ１３で取得した各画像について輝度分布を求める（ステップＳ１４）。典型的には、画像は多数のピクセルから構成されており、当該画像をアナログ解析することにより各ピクセルにおける輝度値を特定することにより、輝度分布が求められる。例えば各ピクセルの輝度値が８ビットで表される場合には、各ピクセルにおける輝度値が０〜２５５の範囲で特定されることとなる。これにより、例えば、各ピクセルの座標データと、各ピクセルの輝度値とが関連付けられた輝度分布が求められる。
【００５４】
尚、ステップＳ１３で取得された画像がカラー画像である場合には、グレースケール変換処理やＹＵＶ変換による輝度情報への変換を行ってからアナログ解析を実施することで、各ピクセルの輝度値を特定してもよい。
【００５５】
続いてステップＳ１４で求めた輝度分布に関する時間的変化を示すプロファイルＰを求める（ステップＳ１５）。本実施形態では、異なる時刻に撮像された各画像の輝度分布から輝度分散値を求め、その時間変化をプロファイルＰとして求める。輝度分散値ＩＤｉｓｐは、例えば平均輝度ＩＡｖｇを用いて次式により求められる。
2020115122.tif 000003
2020115122.tif 000004
ここでＩｘｙは二次元平面上の座標値（ｘ、ｙ）における輝度であり、Ｘ、ＹはそれぞれＸ軸方向及びＹ軸方向の画素数である。
【００５６】
図４は図１のステップＳ１５で求められるプロファイルＰの一例である。図４では、互いに異なる２つの品種（品種Ａ、品種Ｂ）に属する米に関するプロファイルＰＡ、ＰＢが示されている。品種ＡのプロファイルＰＡに着目すると、時刻０で初期値から輝度分散値が増加し、時刻ｔｐＡでピーク値ＤｐＡを示した後、輝度分散値がゼロに向けて次第に収束する振る舞いを示している。このようなプロファイル形状は、初期段階では米１がアルカリ溶液１２を吸収することで膨張したり、米の内部が白化したりすることで輝度分散値が増加した後、溶解が進行するに従って白濁範囲が次第に周囲に広がり、最終的にアルカリ溶液１２に完全に溶解することで均一な透明状態や均一な輝度分布状態に収束する、という一連の振る舞いに対応している。
【００５７】
一方の品種ＢのプロファイルＰＢは、時刻０で初期値から輝度分散値が増加し、時刻ｔｐＢでピーク値ＤｐＢを示した後、輝度分散値がゼロに向けて次第に収束する振る舞いを示している。このように品種Ａ及び品種Ｂは、輝度分散値のピーク値が少なからず異なるものの全体的に類似のプロファイル形状を有する。このようにステップＳ１５で求められるプロファイルＰは、試験対象である米１に関する条件に依存する溶解速度を可視化したものに相当する。
【００５８】
尚、本実施形態では、このように米が溶解する過程において白濁範囲が変化していく振る舞いを輝度分散値という評価パラメータを用いて数値化しているが、当該評価パラメータとして、例えば、各画像における輝度値の最大値と最小値との差分など、他のパラメータを採用することも可能である。
【００５９】
続いてステップＳ１５で求めたプロファイルＰに対して正規化処理を実施する（ステップＳ１６）。図５は図４のプロファイルＰＡ、ＰＢの正規化プロファイルＰＡ’、ＰＢ’である。本実施形態では、プロファイルＰＡ、ＰＢをそれぞれピーク値ＤｐＡ、ＤｐＢで割り算する正規化処理が実施される。このような正規化処理で得られる正規化プロファイルＰＡ’、ＰＢ’では、それぞれのピーク値が「１」に正規化される。このように正規化処理を実施することで、異なる条件に属する米同士を同等に評価することが可能となる。
【００６０】
尚、本実施形態では異なる条件に属する米同士の相対的評価を行うために正規化処理を実施しているが、絶対的評価（他の条件に属する米との比較を伴わない、ある条件に属する米そのものに関する評価）を行う場合には、上述の正規化処理は省略してもよい。
【００６１】
続いて正規化プロファイルＰＡ’、ＰＢ’に基づいて、評価パラメータが基準値に達するまでに要する特性時間を求める（ステップＳ１７）。本実施形態では、評価パラメータである正規化された輝度分散値（正規化処理を省略した場合には、輝度分散値）について、基準値Ｄｒを予め設定し、当該基準値Ｄｒに正規化された輝度分散値が到達した時間を特性時間として求める。図５では、品種Ａに対応する正規化プロファイルＰＡ’が基準値Ｄｒに到達する特性時間ｔＡ、品種Ｂに対応する正規化プロファイルＰＢ’が基準値Ｄｒに到達する特性時間ｔＢが示されている。このように特性時間ｔＡ、ｔＢは、試験対象の米１によって依存するパラメータであり、特性時間ｔＡ、ｔＢに基づいて溶解特性を評価することができる。
【００６２】
本実施形態では、特性時間を求めるための基準値Ｄｒとして、正規化された輝度分散値「０．５（５０％）」が設定されている。図５に示す正規化プロファイルＰＡ’、ＰＢ’の形状に示されるように、正規化された輝度分散値０．５（５０％）の近傍では、変化率（正規化プロファイルＰＡ’、ＰＢ’の傾斜）が大きくなっている。そのため、基準値Ｄｒとして５０％を設定することで、米１の溶解特性を効果的に評価できる。
【００６３】
尚、ステップＳ１７における特性時間は、複数粒の米１を用意した場合には、各粒に対して同等の試験を実施し、それぞれで得られた複数の試験結果を統計的に処理することで求めてもよい。例えば、図２のように同一品種に属する米を複数粒用意した場合には、各粒から求めた特性時間を平均化したものや群の中央値などを、最終的な評価パラメータとして採用してもよい。このように複数の試験結果を統計的に処理したものを最終的な試験結果として取り扱うことで、より試験精度を高めることができる。
【００６４】
続いてステップＳ１７で求めた特性時間に基づいて、品質評価を行う（ステップＳ１８）。ここで図６は幾つかの品種について得られた正規化プロファイルを比較して示す試験結果の一例であり、図７は図６で求められた各品種の特性時間を比較した結果である。図６では、より多くの品種Ａ〜Ｆについて同様に正規化プロファイルが求められており、図７では、図６の各プロファイルから求められた特性時間が品種ごとに比較されている。
尚、図６及び図７では、各品種について同様の試験を２回ずつ実施しており、各回における試験結果を異なるシンボルで区別可能なように示している。
【００６５】
図７では、他の試験結果又は経験則に従って、溶解しやすい順に各品種が並べられている。一方で、図７に示すように、図６のプロファイルに基づいて求められた特性時間もまた左から順に短くなる傾向が示されている。これは、本方法による定量的評価（特性時間に基づく評価）は、他の試験結果や経験則による評価結果と概ね一致していることを意味し、両者の間には高い相関があることが確認された。これにより、本方法で取り扱われる特性時間が、アルカリ崩壊性を評価するためのパラメータとして有効であることが確認された。
【００６６】
ただし図７の品種Ａに関しては、例外的な振る舞いが見られる。つまり、従来手法では品種Ａは品種Ｄより溶けやすいと考えられていたが、本方法による評価では品種Ａは品種Ｄより溶けにくく、逆の結果が得られている。このような例外的振る舞いが生じている要因は明らかではないが、発明者の研究によれば、品種Ａは酒造分野において、いわゆる「後溶け」と称される品種であった。そのため、後溶けに該当する品種については、本方法の評価対象から除外する方が好ましい可能性が示唆されている（確定的なことを言うためには、更なる研究が必要と思われる）。
【００６７】
図８は図６で求められた特性時間とＢｒｉｘ値とを品種毎に比較した結果である。図８では、Ｂｒｉｘ値は公設試験場における測定結果を酒類総合研究所で補正された値を用いており、本方法で得られた特性時間との間に高い相関が確認された。これにより、特性時間とＢｒｉｘ値との相関を予め把握しておくことで、ステップＳ１６で求められた特性時間を用いることで、特性時間からＢｒｉｘ値を推定することができる。
【００６８】
また図９は互いに異なる３品種について求められた特性時間とＤＣＳ値（熱伝導性）とを品種毎に比較した結果である。ＤＣＳ値は主に測定試料（各品種の米）と基準物質との間での熱量を計測して、測定試料の発熱・吸熱の特性変化を観察することで得られ、他の試験結果によって得られた評価値を用いている。具体的には米の場合、米粉水溶液を用い、一定速度で加熱した際に糊状への状態変化が生じる温度に基づいて特定されたものを用いている。図９に示すように、本方法で得られた特性時間は、このように別試験で得られたＤＣＳ値との間に高い相関が確認された。これにより、特性時間とＤＣＳ値との相関を予め把握しておくことで、ステップＳ１６で求められた特性時間を用いることで、特性時間からＤＣＳ値を推定することができる。
【００６９】
尚、ここでは米の品質に関する指標のうちＢｒｉｘ値やＤＣＳ値について言及したが、ＲＶＡ（粘度）のような他の化学的分析値についても同様に相関を有する可能性が期待される。この場合も、特性時間と対象となるパラメータとの相関を解析することで、特性時間に基づいたパラメータの推定が可能となる。
【００７０】
尚、上述のアルカリ崩壊試験に用いられるアルカリ溶液１２として、３％以上（より好ましくは５％以上）のアルカリ濃度のアルカリ溶液を用いることで、加速試験として実施することも可能である。図１０はＫＯＨ濃度がそれぞれ２％、１０％のアルカリ溶液１２を用いた場合に得られた正規化プロファイルの一例である。図１０に示すように、アルカリ溶液１２におけるＫＯＨ濃度を２％から１０％にすることで、試験に要する時間を大幅に短縮することができる。従来の目視観察が必要な試験では、このような加速試験を適用しようとすると溶解ポイントの判定精度が著しく低下してしまい、評価精度が更に低下してしまう。その点、本方法は特性時間による定量的評価を行うため、このようにアルカリ溶液１２を高濃度化した加速試験を適用した場合であっても適切な評価が可能となる。
【００７１】
以上説明したように、試験対象である米が溶解する過程を撮像した画像を処理することによって、アルカリ崩壊性を定量的に評価することができる。このような画像処理は、イメージスキャナのような撮像装置と、画像解析用のソフトウェアがインストールされた演算装置（例えばＰＣなど）を用いて簡易的に行うことができる。これにより、従来は各年度に品種ごとの抜き取りでしか評価されていなかった米の品質が、例えば生産者、生産地等のより細かいロット単位で評価を実施することが可能となる。このような評価結果は、例えば、米の流通時に試験結果を添付データとすることで、消費者に対して客観的な指標として活用することも可能となる。
【００７２】
ここで試験対象である米には、乾燥時に割れが生じていない状態にあっても、吸水（例えば醸造や炊飯）によって内部に潜在的に有するクラックなどが進行することによって割れが後発的に生じる、いわゆる浸漬割れが発生することがある。このような浸漬割れは、米の品種や状態によって所定の割合で発生し、米の品質（食用米では嗜好性、酒用米では醸造条件）と密接に関連することが知られている。
【００７３】
前述のアルカリ崩壊試験方法は、このような浸漬割れによる影響を評価するために、プロファイルに基づいて浸漬割れの有無を判定する工程を更に備えてもよい。図１１はアルカリ崩壊試験を利用した浸漬割れの判定方法を工程毎に示すフローチャートである。
【００７４】
まず浸漬割れの判定対象となる米について、プロファイルを取得する（ステップＳ２１）。ステップＳ２１におけるプロファイルの取得は、図１のステップＳ１１〜Ｓ１６を参照して前述した手順と同様に行われる。簡潔に説明すると、米をアルカリ溶液に浸漬することでアルカリ崩壊試験を実施し、その様子を経時的に撮像して米が溶けてなくなっていく様子を画像解析することによりプロファイルが得られる。
【００７５】
ここで図１２は浸漬割れが生じる前後における撮像イメージであり、図１３は図１２の各サンプルのうち２つのサンプルについて得られたプロファイルを比較して示す測定結果である。尚、図１３では図５及び図６のようにプロファイルをピーク値で正規化した正規化プロファイルが示されている。
【００７６】
図１２では、アルカリ溶液に浸漬された６個のサンプルＡ１〜Ａ６の撮像イメージが示されている。これらのサンプルＡ１〜Ａ６のうちサンプルＡ３は、アルカリ溶液に浸漬しておくと、浸漬割れによって複数の破片に分割されて溶解する様子が示されている。図１３では、図１２において浸漬割れが生じていないサンプルＡ１と、浸漬割れが生じているサンプルＡ３とについて、正規化されたプロファイルＰ１、Ｐ３が比較して示されている。両プロファイルＰ１，Ｐ３の評価方法については詳しくは後述するが、例えば、浸漬割れが発生しているサンプルＡ３では米が複数の破片に分割されることで、アルカリ溶液に対する接触面積が増加し、サンプルＡ１に比べて溶解速度が速くなる傾向を示すプロファイルが得られている。
【００７７】
続いてステップＳ２１で取得したプロファイルに基づいて浸漬割れの有無を判定する（ステップＳ２２）。プロファイルは、各画像に含まれる輝度分布の時間変化に基づく特徴量であることから、前述のように、アルカリ溶液における米の溶けやすさを定量的に評価することができる。浸漬割れは、図１２を参照して前述したように、米が複数の破片に分割されることで、アルカリ溶液への接触面積が急激に変化し、アルカリ溶液中の溶解進行度が影響を受ける。このような影響はプロファイルに反映されることから、プロファイルを分析することで浸漬割れの有無を適切に判定することができる。
【００７８】
ステップＳ２２では、特性時間が予め設定された閾値以下である場合に、米に浸漬割れが有ると判定してもよい。特性時間は、図１のステップＳ１７を参照して前述したように、正規化されたプロファイルに基づいて、評価パラメータが基準値Ｄｒに達するまでに要する時間として求められる。図１３では基準値Ｄｒとして正規化された輝度分散値「０．５（５０％）」が設定されており、プロファイルＰ１に対応する特性時間ｔｓ１と、プロファイルＰ３に対応する特性時間ｔｓ３とが示されている。浸漬割れが生じていないサンプルＡ１に対応するプロファイルＰ１では特性時間ｔｓ１が閾値ｔｓｒより大きいが、浸漬割れが生じているサンプルＡ３に対応するプロファイルＰ３では特性時間ｔｓ３が閾値ｔｓｒ以下になっている。浸漬割れが生じる米はアルカリ溶液に接触する表面積が増えることで溶解の進行度が早くなることから、この態様では、特性時間が閾値ｔｓｒ以下であるか否かに基づいて浸漬割れの有無を適切に判定できる。
【００７９】
図１４は他のサンプルについて取得されたプロファイル（規格化済）に基づく浸漬割れの判定例である。図１４では、異なる４つのサンプルＡ７〜Ａ１０に対応するプロファイルＰ７〜Ｐ１０が示されており、図１３と同様に、基準値Ｄｒ（規化された輝度分散値「０．５（５０％）」）に達するまでの特性時間ｔｓ７〜ｔｓ１０が示されている。この判定例では、サンプルＡ７の特性時間ｔｓ７が閾値ｔｓｒ以下であるため、サンプルＡ７について浸漬割れがあり、残りのサンプルＡ８〜Ａ１０については浸漬割れが無いと判定される。
【００８０】
またステップＳ２２では、プロファイルの時間微分関数のピーク値に到達する時間が閾値以下である場合に、米に浸漬割れが有ると判定してもよい。図１５は図１４のプロファイルＰ７〜Ｐ１０の時間微分関数Ｆｄ７〜Ｆｄ１０のピーク値ＰＥ７〜ＰＥ１０に到達する時間ｔｐ７〜ｔｐ１０を用いた判定例である。時間微分関数Ｆｄ７〜Ｆｄ１０は、時間ｔｐ７〜ｔｐ１０においてそれぞれピーク値ＰＥ７〜ＰＥ１０に達している。閾値ｔｐｒは予め設定される基準値であり、時間ｔｐ７〜ｔｐ９が閾値ｔｐｒ以下となっていることから、これらに対応するサンプルＡ７〜Ａ９について浸漬割れがあると判定され、残りのサンプルＡ１０については浸漬割れが無いと判定される。
【００８１】
ステップＳ２２では、プロファイルの時間微分関数のピーク値が閾値以上である場合に、米に浸漬割れが有ると判定してもよい。図１６は図１４のプロファイルＰ７〜Ｐ１０の時間微分関数Ｆｄ７〜Ｆｄ１０のピーク値ＰＥ７〜ＰＥ１０を用いた判定例である。時間微分関数Ｆｄ７〜Ｆｄ１０は、それぞれピーク値ＰＥ７〜ＰＥ１０を有している。閾値ＰＥｒは予め設定される基準値であり、ピーク値ＰＥ７〜ＰＥ９が閾値ＰＥｒ以下となっていることから、これらに対応するサンプルＡ７〜Ａ９について浸漬割れがあると判定され、残りのサンプルＡ１０については浸漬割れが無いと判定される。
【００８２】
尚、ステップＳ２２では上述の各態様によって判定結果にばらつきがあるが、各態様は状況に応じて適宜使い分けてもよい。例えば、アルカリ溶液への溶解速度を重視して浸漬割れを判定したい場合には、図１４を参照して前述したように、特性時間が予め設定された閾値以下であるか否かによって判定してもよい。またプロファイルの形状を重視して浸漬割れを判定したい場合には、図１５や図１６を参照して前述したように、プロファイルの時間微分関数の特性（ピーク値やピーク値に到達するまでの時間）に基づいて判定してもよい。
【００８３】
また上記の各態様は組み合わせて用いられてもよい。この場合、例えば、少なくとも一つの態様によって浸漬割れがあると判定された場合に、浸漬割れがあるとの最終判定を行ってもよいし、全ての態様で浸漬割れがあると判定された場合に浸漬割れがあるとの最終判定を行ってもよい。
【００８４】
図１１に戻って、続いてステップＳ２２の判定結果を用いて、米の品質評価を行う（ステップＳ２３）。例えばステップＳ２３では、品質評価の一つの指標として、浸漬割れの発生率を求めてもよい。浸漬割れの発生率は、複数のサンプルに対して前述の浸漬割れの判定を実施し、浸漬割れが発生したサンプル数の割合を求めることで算出される。浸漬割れの発生率は、例えば、精米ロット毎に求められてもよい。
【００８５】
またステップＳ２３では、浸漬割れが有るサンプルを特定し、これら浸漬割れがあるサンプルを除いた残りのサンプルに対して品質評価を行ってもよい。これにより、浸漬割れの影響を排除した米の品質評価が可能になる。尚、この場合、残りのサンプルに対して行われる品質評価方法は、前述の各実施形態を適用可能である。
【００８６】
以上説明したように各画像に含まれる輝度分布の時間変化に基づく特徴量であるプロファイルに基づいて、アルカリ溶液における米の溶けやすさを定量的に評価することで、米の内部に存在するクラックのような潜在的要素が発展して生じる浸漬割れの有無を適切に判定することができる。浸漬割れは、米の品種や状態によって所定の割合で発生し、米の品質（食用米では嗜好性、酒用米では醸造条件）と密接に関連することから、これを考慮することで、より高度な品質評価が可能となる。
【産業上の利用可能性】
【００８７】
本発明の少なくとも一実施形態は、アルカリ溶液を用いて米等の溶解特性を評価するためのアルカリ崩壊試験方法に利用可能である。
【符号の説明】
【００８８】
１  米
１０  セル
１２  アルカリ溶液
１４  白濁範囲
【図１】
2020115122.tif 000005
【図２】
2020115122.tif 000006
【図３Ａ】
2020115122.tif 000007
【図３Ｂ】
2020115122.tif 000008
【図４】
2020115122.tif 000009
【図５】
2020115122.tif 000010
【図６】
2020115122.tif 000011
【図７】
2020115122.tif 000012
【図８】
2020115122.tif 000013
【図９】
2020115122.tif 000014
【図１０】
2020115122.tif 000015
【図１１】
2020115122.tif 000016
【図１２】
2020115122.tif 000017
【図１３】
2020115122.tif 000018
【図１４】
2020115122.tif 000019
【図１５】
2020115122.tif 000020
【図１６】
2020115122.tif 000021