(57)【特許請求の範囲】
【請求項１】
検査員の官能試験により、各等級の基準標品を決定する工程と；
非破壊検査により、前記基準標品の輝度、水分、重さ、厚さ、を測定する工程と；前記工程で測定した該基準標品の輝度に基づき輝度指数を求め、前記工程で測定された重量を水分の値を用いて乾物重量に変換し、厚さを該乾物重量で除して粗度指数を求める工程と；
前記各基準標品の等級を等級指数に変換して目的変数とし、輝度指数、粗度指数、を説明変数として重回帰分析を行い、回帰係数を求めて等級付け用の簡易重回帰式を作成する工程と；
非破壊検査により、被検乾海苔の輝度指数及び粗度指数、を求める説明変数実測工程と；
前記説明変数実測工程で求めた各値と前記簡易重回帰式とに基づき、該被検乾海苔の等級指数を求める工程と；
を備えていることを特徴とする乾海苔の等級決定方法。
【請求項２】
検査員の官能試験により、各等級の基準標品を決定する工程と；
非破壊検査により、前記基準標品の輝度、色合い、タンパク質含有率、水分、重さ、厚さ、を測定する工程と；
前記工程で測定された該基準標品の輝度に基づき輝度指数を求め、色合いに基づき色合いの均一性指数を求め、又、前記工程で測定された重量を水分の値を用いて乾物重量に変換し、厚さを該乾物重量で除して粗度指数を求める工程と；
前記各基準標品の等級を等級指数に変換して目的変数とし、輝度指数、色合いの均一性指数、タンパク質含有率、粗度指数、を説明変数として重回帰分析を行い、回帰係数を求めて等級付け用の重回帰式を作成する工程と；
非破壊検査により、被検乾海苔の輝度指数、色合いの均一性指数、タンパク質含有率及び粗度指数、を求める説明変数実測工程と；
前記説明変数実測工程で求めた各値と前記重回帰式とに基づき、該被検乾海苔の等級指数を求める工程と；
を備えていることを特徴とする乾海苔の等級決定方法。
【請求項３】
求めた被検乾海苔の等級指数を量子化して等級を決定することを特徴とする請求項１、又は、２記載の乾海苔の等級決定方法。
【請求項４】
等級付け用の簡易重回帰式が、等級指数＝定数項＋粗度の回帰係数×粗度指数＋輝度の回帰係数×輝度指数であることを特徴とする請求項１記載の乾海苔の等級決定方法。
【請求項５】
等級付け用の重回帰式が、等級指数＝定数項＋粗度の回帰係数×粗度指数＋タンパク質含有率の回帰係数×タンパク質含有率＋輝度の回帰係数×輝度指数＋色合いの均一性の回帰係数×色合いの均一性指数であることを特徴とする請求項２記載の乾海苔の等級決定方法。
【請求項６】
簡易重回帰式が、被検乾海苔の産地と時期が近似している場合に用いられることを特徴とする請求項１記載の乾海苔の等級決定方法。
【請求項７】
全自動乾海苔製造装置を使用して製造された直後の被検乾海苔であって、該被検乾海苔の産地と時期が近似している場合には、水分の測定を省略し、経験値を用いることを特徴とする請求項１、又は、２記載の乾海苔の等級決定方法。
【請求項８】
非破壊検査が、一帖の表面１枚の被検乾海苔に対して行われることを特徴とする請求項１、又は、２項記載の乾海苔の等級決定方法。
【請求項９】
非破壊検査が、一帖の全部について行われ、それらの平均値を測定値とすることを特徴とする請求項１、又は、２記載の乾海苔の等級決定方法。
【請求項１０】
各等級の基準標品を採取した検査場で、同一の等級付け用の重回帰式又は簡易重回帰式を用いて、同日に被検乾海苔の等級付けを行うことを特徴とする請求項１、又は、２記載の乾海苔の等級決定方法。
【請求項１１】
各等級の基準標品を採取した検査場で、該基準標品を採取した日と異なる日に、以前使用した等級付け用の重回帰式又は簡易重回帰式を用いて、被検乾海苔の等級付けを行うことを特徴とする請求項１、又は、２記載の乾海苔の等級決定方法。
【請求項１２】
各等級の基準標品を採取した検査場と異なる検査場で、かつ、該基準標品を採取した日と異なる日に、前記基準標品の等級付け用の重回帰式又は簡易重回帰式を用いて、被検乾海苔の等級付けを行うことを特徴とする請求項１記載の乾海苔の等級決定方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、乾海苔の等級決定方法に関するものであり、更に述べると、海苔養殖業者の製造した乾海苔を流通に乗せる過程で実施される等級検査に於いて、検査員の視覚、触覚、に頼らずに、機械を使用し、非破壊で能率良く、しかも、客観的結果をもたらせる、乾海苔の等級決定方法に関するものである。
【０００２】
【従来の技術】
我が国では現在１漁期中に約１００億枚の乾海苔が生産され、販売されている。乾海苔の生産地は、東日本から九州に至る表日本の沿岸であり、消費地は全国である。従って、漁村において生産された乾海苔は、生産地に近い拠点に集荷され、卸売業者に買い取られ、ついで卸売業者から小売業者の手を経て消費者へ流通される。天産品である乾海苔の品質は多様性に富んだものであるから卸売業者は本来なら１枚づつ品質を見て値を付けなければならないが、短い漁期中に集荷拠点に集められる乾海苔の量は非常に多いので、その様なことは困難である。
【０００３】
そこで、生産者は同程度の品質を持つ製品（乾海苔）を集めて格付けし仕訳をして、荷口ごとに見本をつけて卸売業者の検品に供する。
卸売業者は、見本を見て値踏みをし入札によって製品を入手する。
この時の格付け作業が検査であって、もし検査の機能が麻痺すると、取引はきわめて能率の悪いものになり、現在のような大量の生産物を捌くことが不可能になる。
【０００４】
従来の品質検査は、習熟した乾海苔の検査員が行っている。この検査員は、検査対象の乾海苔（被検乾海苔）を１枚ずつ手にとり、視覚、触覚により色、艶、出来、などを総合的に検査した後、１０枚１帖を１０帖束ねた１束ごとに巻いた帯に等級、例えば、味推上１等、味推１等、推１等、上１等、１等〜７等、などの印を打つ。
【０００５】
【発明が解決しようとする課題】
しかし、従来例では次の様な問題がある。
（１）乾海苔は通常、検査場へ３６００枚入りの段ボール箱に納められて集荷される。検査場では、検査員はまず梱包を解き、検査台の上に並べて１帖づつ検査し、検印を押した後再び梱包する。
この検査に要する時間は１枚あたり約０.１秒程度であり、１箱当たり約６分かかるので、１日１人２５０箱程度が限界である。
【０００６】
ところが、乾海苔の１次加工は全自動機を使用して行われていて、小型のものでも１日に３０００枚程度の製品を作っているが、最盛期になると、３６００枚位の製品が作られる。
従って、組合単位で１０日に１回の集荷を行う検査場では１回の検査は１０００箱程度になり、２人の検査員で検査しても２昼夜を要することになる。
実際の検査は、１日で終わるように計画されているが、梱包解き、印打ち、梱包、記録等の作業は早朝から深夜におよび、その労力は多大である。
【０００７】
（２）長時間検査していると、検査員の眼の疲労が著しくなり、色、艶の判断がにぶり、常に一定の判断の基準に照らしたものとはなり難くなる。
そのため、しばしば判定結果に不平等を生じ、売り手、買い手の双方に不満の原因を作って公平な取り引きを阻害する。
【０００８】
（３）検査員の検査可能な範囲は、前述の様に通常１日１人２５０箱位であるので、これよりも検査対象の乾海苔が多くなると、１日で検査が完了する様に、乾海苔を検査台の上に流しておいて、気になるところだけ抜き取って確認する方法が採用される。
しかし、この方法では、不良品の混合を見逃し、信用を失う弊害を生じる。
【０００９】
（４）検査員が複数の場合には、検査員の能力の差による判定誤差が発生する恐れがある。但し、異なった検査場の結果は相互に独立して取り扱われるので、問題にはされてはいるが、容認されているのが現状である。もちろん、そのことも売り手買い手の双方に不満の原因を作っている。
【００１０】
（５）乾海苔の検査は、高度の品質鑑定能力が必要である。ところが、この能力を備えたベテランの検査員は老齢化する一方、その後継者がなかなか育たないのが現状である。そのため、必要な数の検査員を確保することが困難となると共に、検査員１人当たりの検査枚数が飛躍的に増加し、過重な労働を強いる結果となっている。
【００１１】
前記問題を解決する目的で、重量、水分、色、艶、タンパク質含有率、等個々の形質を非破壊で計測する機械が開発され、一部で使用されている。
それらの中には、単一形質の計測値の再現性に関しては十分に実用の域に達しているものも見受けられる。
【００１２】
しかしながら、乾海苔の等級は、形質の総合として決定されるものであるため、上記の計測器は、何れも検査の補助的な役割を担うに過ぎず、検査員の不安を若干軽減する程度の効果はあっても、労力の軽減、検査精度の向上には抜本的に貢献するものでは無い。
【００１３】
この発明は、上記事情に鑑み、乾海苔の等級決定の作業能率を向上させるとともに、客観的で公正な評価を行うことを目的とする。
【００１４】
この発明は、検査員の官能試験により、各等級の基準標品を決定する工程と；非破壊検査により、前記基準標品の輝度、水分、重さ、厚さ、を測定する工程と；前記工程で測定した該基準標品の輝度に基づき輝度指数を求め、前記工程で測定された重量を水分の値を用いて乾物重量に変換し、厚さを該乾物重量で除して粗度指数を求める工程と；前記各基準標品の等級を等級指数に変換して目的変数とし、輝度指数、粗度指数、を説明変数として重回帰分析を行い、回帰係数を求めて等級付け用の簡易重回帰式を作成する工程と；非破壊検査により、被検乾海苔の輝度指数及び粗度指数、を求める説明変数実測工程と；前記説明変数実測工程で求めた各値と前記簡易重回帰式とに基づき、該被検乾海苔の等級指数を求める工程と；を備えていることを特徴とする。
【００１５】
この発明は、検査員の官能試験により、各等級の基準標品を決定する工程と；非破壊検査により、前記基準標品の輝度、色合い、タンパク質含有率、水分、重さ、厚さ、を測定する工程と；前記工程で測定された該基準標品の輝度に基づき輝度指数を求め、色合いに基づき色合いの均一性指数を求め、又、前記工程で測定された重量を水分の値を用いて乾物重量に変換し、厚さを該乾物重量で除して粗度指数を求める工程と；前記各基準標品の等級を等級指数に変換して目的変数とし、輝度指数、色合いの均一性指数、タンパク質含有率、粗度指数、を説明変数として重回帰分析を行い、回帰係数を求めて等級付け用の重回帰式を作成する工程と；非破壊検査により、被検乾海苔の輝度指数、色合いの均一性指数、タンパク質含有率及び粗度指数、を求める説明変数実測工程と；前記説明変数実測工程で求めた各値と前記重回帰式とに基づき、該被検乾海苔の等級指数を求める工程と；を備えていることを特徴とする。
【００１６】
【発明の実施の形態】
ベテランの検査員は、乾海苔の格付け（等級付け）を行う場合、色（色合いの均一性）、艶（輝度）、目方（重量）、柔らかさ、焼き色、等を判断データとして、視覚、触覚により総合的に判断している。
本発明者は、これらの判断データに対応する乾海苔のデータを、被破壊検査法、即ち、検査対象の乾海苔（被検乾海苔）を破壊しないで検査する方法、により測定して該測定データを計数処理すれば、客観的に、乾海苔の等級付けができる、と考えた。そこで、等級付けするための理論（計数処理）を研究し、次のような実験を行ったところ、検査員の等級付けと同一レベルで、被検乾海苔の等級付けができることがわかった。
【００１７】
乾海苔の等級の相対性：
表１に示す１９等級、４８ロットの標品（乾海苔）を入手し、これらを材料にして、１ロットにつき１０枚の標品に対して非破壊検査法により輝度、色合いの均一性、タンパク質含有率、水分、重量及び厚さ、の６項目を測定した。  この表１は、「平成１１年度Ｍ漁協海苔品質評価  標品一覧」で、縦に等級、横に入札回及び入札日、が記載されている。
又、等級は、ロットに付けられる階級名であり、味推上一等から７等まであり、同じ等級であれば「推」の付くものはやや上位であり、「味」の付くのは更に上位となる。
なお、ロットとは、共販所における取引で同一品質を保証する複数の製品を一纏めにしたものをいう。
【００１８】
【表１】
0003765005.tif 000002
【００１９】
測定された重量を水分の値を用いて乾物重量に換算し、厚さを乾物重量で除して得られた値を粗度指数、と名付けて指標形質に加えた。
【００２０】
二次元配置分散分析；
表１において、第１回共販（入札日  平11.11.23)から第５回共販（入札日  平12.01.27）まで揃って出現する等級、即ち、上３等から４等までの４等級、について、輝度、色合い及びタンパク質含有率、の３項目について、それぞれ等級と入札回を因子として二次元配置分散分析によって等級間、入札回及び交互作用の有意差を検討した。その結果、何れの項目も危険率１％以下で帰無仮説は棄却された。
【００２１】
このことから海苔の等級は、至って相対的なものであって、産地が同じでも時期が異なれば同一等級が、個々の形式については同一品質を保証していないことがわかった。
【００２２】
乾海苔の等級の類縁関係；
等級が比較的多数揃っている表１の第１回共販の分について形質の意味を検討し、それらの測定値の総合効果と等級との関係を調べた。
検査では、表面観察のほとんどは１帖の表面１枚に限られているので、この表面の一枚目の乾海苔だけについて、画像解析相対評価総合による類縁関係、相関係数、群間平均法でクラスタ分析を行った。
【００２３】
その結果、
イ）味推１等を除く上３等以上のクラスタ、ロ）３等と上５等で構成されるクラスタ、ハ）上４等から５等までで構成されるクラスタ、
に分割された。特に、上４等から５等のクラスタは、上３等以上のクラスタとは類似性が低いことが示されたので、上級のものと下級のものでは色の構成が異なることが明らかになった。
【００２４】
輝度、色合い、タンパク質含有率及び粗度の４変数による類縁関係、標準ユークリッド、ワード法で、クラスタ分析を行った。
その結果、味推上２等がタンパク質含有率の低さ故に下等グループに位置ずけられたことを除くと、類似度７.０で推２等以上と上３等以下に分割され、類似度２で、５クラスタに整理された各クラスタの内訳は、検査等級にほぼ合致した。
【００２５】
乾海苔のうまみと外観；
表１の第１回共販分の一帖の表面１枚に関し、等級指数と水分を項目に加えて各指標主成分分析を行ったところ、２成分で累積寄与率が９５％に達した。主成分１の寄与率は７１％で、因子負荷量をみると、等級指数とタンパク質含有率が有意な負の値で関わり、他は正の値で関わっているが、水分の負荷重は他の項目の２〜３分の１、に過ぎない。
主成分２の寄与率は、２４％で、因子負荷量をみると、等級指数との関わりは小さく、タンパク質含有率が比較的大きな正の値を示した。
【００２６】
これから推測すると、主成分１は外観、主成分２はうまみの成分であると判断される。この主成分分析では、外観はうまみにほとんど関係がないが、うまみのある海苔は、タンパク質が多く、かつ、色、艶も一般に良いのが通則であること、検査では、うまみよりも外観が重視されていること、がうまく説明されている。
【００２７】
乾海苔の計算による等級決定；
等級に連続量を仮定した等級指数を求め、それを目的変数とする重回帰分析を行った。ここで用いられる重回帰分析の詳細な説明は省略するが、回帰係数を求め、下記等級付け用の重回帰式を作成した。
目的変数（等級指数）＝切片（定数項）＋回帰係数（α1）×説明変数（Ｘ1）＋回帰係数（α2）×説明変数（Ｘ2）＋………
この等級指数は、まず７等を２とし、１等上がるごとに１を加え、上には0、5、推には0、2、味には0、1を加えることにより求めた。
【００２８】
この重回帰式の説明変数には、主成分分析の結果からみて影響を無視できると考えられる水分を除き、粗度指数（Ｘ1）、タンパク質含有率（Ｘ2）、輝度指数（Ｘ3）、色合いの均一性指数（Ｘ4）、及び、輝度×色合い／100指数（Ｘ5）、の５項目を選び、それらを組み合わせて説明変数４の場合、３の場合、２の場合につき全部で１３ケースを計算し、その結果表２にまとめた。
【００２９】
【表２】
0003765005.tif 000003
【００３０】
これら１３ケースの計算では、全ての場合について分散分析の結果、危険率１％以下であり、相関関係が認められた。
この１３ケース中、表２のどの番号の式が最適であるか、を検討した。
その結果、残差平方和の小ささと重相関係数の大きさからは、４変数（粗度、タンパク質含有率、輝度、色合い）の場合（式１.１）が最適であり、又、自由度により調整された重相関係数は前記４変数でも２変数（粗度、輝度）の場合（式３.２）でもほとんど差が無く、ＡＩＧ及びＣpの値からは、２変数の方がむしろ適していることがわかった。
【００３１】
重回帰分析による等級決定の妥当性：
重回帰分析を使用する等級決定法の演繹性について検討した。
前記式1.1及び式3.2（基準標品の１帖の表面１枚だけに基づいて作成された重回帰式）に、第１回共販等級別各測定値平均（一帖の全部１０枚について測定し、それらの測定値を平均した値）を代入した結果、等級指数と予測値との相関係数は、前者が0.919、後者が0.890で大差はなかった。
【００３２】
一方、第１回入札等級別平均値を元にした重回帰分析では、重相関係数0.92が得られ、この式による予測値と先の式による予測値との相関は、0.9982であったことから、一枚目の標品測定値から求められる重回帰式は、同一入札回の等級を同じにする全標品に適用しても全く問題がないことがわかった。
【００３３】
検査時期の相違：
第５回共販の粗度タンパク質含有率の測定値を第１回共販と比較検討した。  粗度は、等級指数にほぼ反比例し、タンパク質含有率はほぼ正比例していることがわかった。
【００３４】
そこで、仮にこれらを直線回帰とみなして実測値と予測値の残差平方和を求め、等級数で除すると、第１回共販の粗度では、0.14、タンパク質含有率では、9.84、第５回共販の粗度では、2.30、タンパク質含有率では、12.85が保たれ、特に、第１回共販に対し第５回共販の粗度の変動が大きかった。
【００３５】
相関係数は、第１回共販の粗度、タンパク質含有率、とも比較的高い値を示したのに対し、第５回共販では、低い値であった。
これらの数値は、摘採回が進んで、原藻の葉質が低下し、高級贈答品としての販路が閉ざされると、製造が粗雑になり、葉質の良否が製品の仕上がりに反映しにくくなることを示しているものと考えられる。
【００３６】
結  論：
以上の実験の結果、本発明者は、次のことを確認することができた。即ち、
(1)乾海苔の等級検査を機械で置き換える際に必要な測定項目は、輝度、色合い、タンパク質含有率、水分、重量及び厚さ、の６項目であること。
(2)等級指数を目的変数とし、粗度、タンパク質含有率、輝度及び色合いを説明変数とする重回帰分析のモデル（重回帰式）が成り立つこと。
【００３７】
(3)等級が比較的揃っている共販回の標品を用いて重回帰分析を行った場合、回帰係数はその共販回に限らず適用可能であるが、予測計算結果は検査員の評価と若干そして時には大きな誤差が生じることがある。
この誤差の生じる理由は、検査員の感覚的な問題による場合と、製品の品質構成の内容が大きく変化する場合とがあるが、後者の場合は回帰係数をその共販回の資料から求めて計算し直す必要があること。
本発明は、前記知見に基づいてなされたものである。
【００３８】
【実施例】
この発明の実施例を図１〜図３により説明する。
被検乾海苔：
検査対象にした乾海苔（被検乾海苔）は、表３に示すように、全部で５５ロットあり、Ｍ漁協検査場に時期を違えて５回集荷され、検査員の視覚、触覚による試験、即ち、官能試験、により１９等級に仕分けられたものである。それぞれの入札日（５回）に、等級別に産地問屋による入札によって落札単価が付けられている。
【００３９】
【表３】
0003765005.tif 000004
【００４０】
等級基準標品の選択：
検査官は官能試験を行いながら、第１回入札（平成１１年１１月２３日）分の１４ロットの乾海苔の中から、各等級（味推上１等、味推１等、推１等、１等〜）の基準となる標品を１束（１０帖、１００枚）づつ抜き出し、等級基準品を決める。この１束のうちの１帖も、又、一帖の一番外側（表面）の１枚もほとんど品質が同じなので、同一等級とみなすことができる。
【００４１】
等級基準標品の判断データの測定：
各等級の基準標品の輝度、色合いの同一性、タンパク質含有率、水分、重量及び厚さ、の６項目を測定する。
この測定は、検査装置により標品を破壊しないで行う検査、即ち、非破壊検査法により行う。
【００４２】
輝度及び色合いの同一性は、ＣＣＤ等をセンサとして通常の光学的手法によってコンピュータに取り込まれた乾海苔の画像情報から計算によって求める。この測定手段として、例えば、先の出願（特願平９−１５２３７５号、特許第３１６８４０号）で使用した検査装置を、前記出願で述べた要領で使用した。なお、前記測定手段としてカラースキャナを用いても良い。
この検査装置は、標品１を載置する検査台２と、該標品１に光を照射する光源３と、該検査台２と対向する撮影装置５と、該撮影装置５にＡ／Ｄ変換器４を介して接続されたコンピュータ７と、を備えている。
【００４３】
光源３は、波長特性並びに照度が整えられている。通常の色を検出するためには、可視光領域を網羅し、可視光領域内において、各波長のエネルギー量が自然光のエネルギー量に比例していることが理想的である。
【００４４】
撮影装置５には、カラーＣＣＤ６を有するビデオカメラが設けられている。
このＣＣＤ６は、電荷結合素子（Charge  Coupled  Device）であり、乾海苔の反射光をレンズを通してＣＣＤ受光面に焦点を結ばせ、ドットサイズが海苔の断片よりも十分小さい、ＲＧＢ（赤、緑、青）の３種それぞれの光量に比例するアナログ信号電流として出力する。
【００４５】
この電流は、ＡＤ変換器４によりデジタル化された後、コンピュータ７のＲＡＭへ転送される。なお、計算速度を向上させるために、ドットの間隔が海苔の断片のサイズに相当する程度にデータの省略を行う。
【００４６】
ここで前記「海苔の断片のサイズ」について説明する。
海苔は乾海苔に加工される過程で、淡水で洗浄後ミンチと称する機械により細断される。裁断された断片のサイズは特定されていないが、通常長短辺ともに１ないし２ミリ程度である。その後の製法は製紙の方法とほぼ同じで、裁断された海苔を水に懸濁させ、均一に分散するように撹拌しながら、水切れのよい廉の上に流し込み、水を切って乾燥する。
【００４７】
乾海苔は、細断された断片が折り重なってできているので、乾海苔の中の断片の数は無数と言う表現によるのが適当なほどに極めて多い。
なお、ＣＣＤ６が取り出す情報の１ドット分が占める大きさは直径０.１ミリ以下である。
【００４８】
この撮影装置５では、ＣＣＤ６は検査対象の乾海苔１の反射光Ａを取り込み、該乾海苔１の色（色合いの均一性）及び艶（輝度）を数値化する。
この数値化は、前記RGB配列データをＹｕｖ変換式を用いて計算し、Ｙｕｖ色座標数値データに変換する。この数値化の詳細については、前記出願と同じなので、ここでは説明を省略する。
【００４９】
輝度指数：
前述のように、ＲＧＢ配列をＹｕｖ色座標計数値配列に変換し、Ｙの頻度分布を統計処理して輝度指数を求める。
被写体（乾海苔）を斜め方向（通常４５度）から照射して撮影すると、被写体の表面で乱反射が生じる時、大きな値となり、鏡面のように乱反射のない表面からは、反射光が入らないため０、又は、０に近い値となる。
【００５０】
乾海苔の表面は比較的なめらかで艶があるので、乱反射は少なく、Ｙの値は０に近いのが普通であり、６４以上になることはない。
そこで、Ｙ＝０を１００、Ｙ＝６４を０、になるような反比例配分した値を求め、その平均値を求めて前記輝度指数とする。
【００５１】
色合いの均一性指数：
ｕｖ座標の頻度分布を統計処理して色合いの均一性指数を求める。
計算されたｕ、ｖを色座標を水平面に、頻度を垂直面又は、鉛直面にした立体座標上にプロットすると、山脈状の図が得られる。この時、色合いが均一であればあるほど山の高さが高く、すそ野の広がりは狭くなる。
【００５２】
そこで、頻度が１０以上あるｕ、ｖ座標値の頻度合計数を調べ、全数に対する割合を求めて色合いの均一性指数とした。即ち、全ての頻度値が１０以上である場合を１００、全ての頻度値が１０未満である場合を０とする代表値を求めて色合いの均一性指数とする。
【００５３】
図２は、画像解析結果の例を示すもので、一枚の味推上一等の標品の輝度指数及び色合いの均一性指数を示している。
この図２の左側の図では、縦は頻度、横は輝度あり、又、右側の図では、衝立のように書かれている面にある横線が頻度で、平面の左下から右上へ向かう線がｕで、右下から左上へ向かう線がｖである。両者の交点から左、やや上に向かって並んで立っている棒状体が（ｕ、ｖ）値の頻度である。
【００５４】
なお、図２の下部に表示されてる、絶対評価値とは、その標品だけの測得値から計算される値で、輝度指数と色合いの均一性指数のことであり、また、相対評価値とは基準標品と比べた値で、被検乾海苔がどの基準標品と似ているのかを確かめるのに大変便利な指数である。
【００５５】
タンパク質含有率と水分は、海苔成分計（JTエンジニアリング社製MJ-２０００）を使用して測定した。重量は電子天秤１０、又、厚さはキャリパ１１で測定した。タンパク質含有量の値（タンパク質含有率）は、乾海苔の乾燥重量に対するタンパク質重量のパーセンテージ（％）であり、この値は連続量として得られるので、そのまま説明変数に使用する。
【００５６】
厚さを乾物重量で除して得られる値を粗度指数と定義した。
前記乾物重量は重量から水分の占める重量を差し引いて求める。
なお、全自動乾海苔製造装置を使用して製造された直後の製品では、経験的に知られている水分の値（経験値）を代入して得られる近似値で代用できる。
【００５７】
第１回入札分、一枚目非破壊検査結果；
第１回目入札（平11.11.23）における等級基準標品（１３ロット）の各測定データは、表４の通りであり、例えば、味推上１等の粗度指数は3.31、タンパク質含有率54.4、輝度指数83、色合いの均一性指数は76、水分は4.6、であった。
【００５８】
【表４】
0003765005.tif 000005
【００５９】
重回帰分析：
検査員の等級付けは、離散量であるから、連続量である等級指数に変換した上で目的変数とし、粗度指数、タンパク質含有率、輝度指数、色合いの均一性指数、を説明変数として重回帰分析を行い、回帰係数を求めて下記モデル式（重回帰式）を作った。
等級指数=21.656（定数項）＋1.221（粗度の回帰係数）×粗度指数＋0.064 （タンパク質含有率の回帰係数）×タンパク質含有率＋0.185（輝  度の回帰係数）×輝度指数＋0.035（色合いの均一性の回帰係数）  ×色合いの均一性指数
【００６０】
前記重回帰式に表４の各値を代入し演算したところ、表５の結果を得た。
【表５】
0003765005.tif 000006
【００６１】
表５では、縦に項目、横に係数等をとっているが、ｐ値は、各項目の分散から見た危険率で、小さいほど精度がよい。残差平方和は、実測値と重回帰式で求めた予測値との差の二乗を加算したものである。
【００６２】
重相関係数は実測値と予測値の相関の程度を示すもので、この値の二乗が一致率になる。調整済み係数は、重相関係数を自由度で調整して算出した係数で、求めた重回帰式を無限母集団のデータ群に適用したときの重相関係数の最低値と考えることができる。AICは赤堀の判定指数で、他の結果が同じ様な値であるときは、この値の小さい方を選ぶことになっている。
【００６３】
この様にして得られた等級指数を量子化して等級を決定する。
ここで、量子化とは、アナログデータをデジタルデータに変換することを意味する。等級は、味推上１等、味推１等、等と分けられ、それぞれの間は繋がっていないので離散量であるが、これに対し、等級指数は、8.8、8.3、などと書かれても連続量である説明変数を組み込んだ計算で求められる値なので、8.8、8.3の間には端数の付いた数が無数にあり、切れ目がないので連続量となる。この連続量、例えば、8.7772……、を四捨五入して8.8にするのが量子化である。
【００６４】
基準等級指数と等級との対応関係は、例えば、次の様にして決められる。
９から等級の数を引いた値、をその等級の等級指数とし、「上」には0.5、「推」には0.２、「味」には0.1、を更に加算する。例えば、味推上１等の等級指数は、  ９−１＋0.5＋0.2＋0.1＝８．８
となる。
【００６５】
次に、非破壊検査により、各帖の乾海苔の輝度、色合い、タンパク質含有率、水分、重さ、厚さ、を測定する。この時、１帖に纏められている１０枚の被検乾海苔全部について測定を行うのではなく、その一番外側（表面）の１枚に付いてだけ行う。
【００６６】
前記測定された被検乾海苔の輝度に基づき輝度指数を求め、色合いに基づき色合いの均一性指数を求め、又、重量、水分、厚さに基づき粗度指数を求める。
【００６７】
前記各被検乾海苔の輝度指数、色合いの均一性指数、タンパク質含有率及び粗度指数、の値を前記重回帰式に代入して演算し、各被検乾海苔の等級指数を求める。この求めた等級指数を量子化し、その値を前記基準等級指数と比較して等級を決定する。
【００６８】
例えば、求められた等級指数が8.7772である場合には、四捨五入して8.8に電子化されるが、この値は味推上１等の8.8に一致するので、該被検乾海苔の等級は、味推上１等、に等級付けされる。
【００６９】
前記重回帰式により求めた被検乾海苔の等級指数（予測値）と、検査員の官能試験により求めた等級指数（実測値）と、を比較したところ、図３に示す結果となった。図３において、横軸は等級指数実測値、縦軸は予測値、Ｒは重相関係数、Ｒ2は、重相関係数の説明率、を示す。
この図３から明らかなように、前記予測値は極めて高い精度で等級指数実測値を再現していることが確認できた。
【００７０】
この発明の第２実施例を図４により説明するが、この実施例は、プライベート等級に関するものである。即ち、
第１実施例では、被検乾海苔を測定する場合に、１帖の束のうち一番外側の１枚目だけを測定し、その値を重回帰式に代入して等級指数を求めたが、この第２実施例では、同日、同一検査場における、同じロット（この場合は第１回入札分）の等級別に１帖の乾海苔全部を計測して、等級別に求めた１０枚の乾海苔の測定値平均（１０枚平均値）を、前記重回帰式に代入して等級指数を計算した。
【００７１】
なお、第１回入札分非破壊検査における１０枚平均値は、表６に示す通りであるが、その値は、表４の測定値と僅かに異なっている。
【表６】
0003765005.tif 000007
【００７２】
この１０枚平均値を用いて求めた等級指数を、プライベート等級指数（予測値）と呼ぶことにするが、このプライベート等級指数を、等級指数実測値と比較したところ、図４に示す結果となった。
この図４から明らかなように、前記予測値は極めて高い精度で等級指数実測値を再現していることが確認できた。
なお、図４において、横軸は等級指数実測値、縦軸は予測値（プライベート等級指数）、Ｒは重相関係数、Ｒ2は、重相関係数の説明率、を示す。
【００７３】
本発明の第３実施例を図５、図６により説明するが、この第３実施例は、ローカル等級指数に関するものである。即ち、
このローカル等級指数は、表４の第１回入札分の等級指数を求めた、第１実施例の等級付け用重回帰式に、表７に示す第３回入札分非破壊検査（１０枚平均値）の測定値を代入して等級指数を演算し、等級を求めるものである。
この様に、同一検査場での検査員の等級評価を他日の標品に適用する処理を、ローカル処理と呼び、この処理で計算される等級をローカル等級と呼ぶ。
【００７４】
【表７】
0003765005.tif 000008
【００７５】
このローカル等級指数（予測値）と等級指数実測値とを比較したところ、図５に示す様に、上質の製品の等級指数が実測値よりも高い等級として算出された。これは、表３に示されている第３回入札分の落札価格を参照すると、検査員の感覚による誤差を修正した値になっていること、を示しているものと認められる。
【００７６】
図６は、前記第１回入札分の重回帰式を使用せずに、独自に第３回入札分の同じ製品群の重回帰式を新たに作り、表７の測定値を代入して求めた被検乾海苔の等級指数（予測値）、即ち、プライベート等級指数、と等級指数実測値とを比較したものである。
この図６から明らかな様に、前記予測値は、第１回入札分のプライベート等級指数（図３）と同様に、極めて高い精度を再現している。
【００７７】
この発明の実施例は、前記に限定されるものではなく、例えば、次のようにしても良い。
(１)グローバル等級；
例えば、平成1１年１1月２３日に、Ｍ漁協検査場で等級指数を求める際に利用した重回帰式を、平1１年1１月３０日に行われた、Ｓ漁協検査場での等級指数の演算に用いることも可能である。
この様に、検査日及び検査場が異なる場合、即ち、ある日にある検査場における検査員の等級評価を、他日の他の検査場の標品に適用する処理をグローバル処理と呼び、この処理で計算される等級をグローバル等級と呼ぶ。
【００７８】
(２)測定項目の省略：
1)等級付け用の簡易重回帰式
タンパク質含有率は、産地と時期の近似した製品では、輝度と高い相関関係を示すことから測定を省略できる。
この場合は、等級付け用の重回帰式は、下記の様な簡易重回帰式となり、説明変数は、輝度指数と粗度指数だけとなるので、計算が容易となる。
等級指数＝定数項＋粗度の回帰係数×粗度指数＋輝度の回帰係数×輝度指数
【００７９】
2)水分：
水分は、全自動乾海苔製造装置を使用して製造された直後の製品では産地と時期が近似している場合、ほぼ一定である。そのため、該水分の測定を省略し、経験的に知られている値（経験値）を適用することができる。
【００８０】
【発明の効果】
この発明は、乾海苔の等級付け用の重回帰式又は簡易重回帰式に、説明変数である輝度指数、粗度指数等を代入することにより、客観的に乾海苔の等級付けを行うことができる。そのため、乾海苔の等級付けを、公正に、かつ、能率良く行うことができる。
また、一帖の全部の乾海苔を測定することなく、その表面の一枚だけを測定し、その測定値に基づいて前記等級付けの式に代入しても、正確な等級付けができる。
【図面の簡単な説明】
【図１】本発明の第１実施例を示す図である。
【図２】輝度指数と色合いの均一性指数を示す画像解析結果を示す図である。
【図３】予測値と等級指数実測値との比較を示す図である。
【図４】本発明の第２実施例を示す図で、プライベート等級予測値と等級指数実測値との比較を示す図である。ある。
【図５】本発明の第３実施例を示す図で、ローカル等級予測値と等級指数実測値とを比較する図である。
【図６】本発明の第３実施例を示す図で、プライベート等級予測値と等級指数実測値との比較を示す図である。
【符号の説明】
１    標品
２    検査台
３    光源
５    撮影装置
７    コンピュータ
【図１】
0003765005.tif 000009
【図２】
0003765005.tif 000010
【図３】
0003765005.tif 000011
【図４】
0003765005.tif 000012
【図５】
0003765005.tif 000013
【図６】
0003765005.tif 000014