▶ 学校法人慶應義塾の特許一覧 ▶ 株式会社チノーの特許一覧
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公開特許公報(A)
(11)【公開番号】P2023128967
(43)【公開日】2023-09-14
(54)【発明の名称】バナナ追熟予測装置およびバナナ追熟予測方法
(51)【国際特許分類】
A23N 15/06 20060101AFI20230907BHJP
G01N 21/27 20060101ALI20230907BHJP
【ＦＩ】
A23N15/06 B
G01N21/27 Z
【審査請求】未請求
【請求項の数】2
【出願形態】ＯＬ
(21)【出願番号】P 2022033675
(22)【出願日】2022-03-04
(71)【出願人】
【識別番号】598121341
【氏名又は名称】慶應義塾
(71)【出願人】
【識別番号】000133526
【氏名又は名称】株式会社チノー
(74)【代理人】
【識別番号】100067323
【弁理士】
【氏名又は名称】西村 教光
(74)【代理人】
【識別番号】100124268
【弁理士】
【氏名又は名称】鈴木 典行
(72)【発明者】
【氏名】山田 秀
(72)【発明者】
【氏名】上田 真由
(72)【発明者】
【氏名】石橋 政三
【テーマコード（参考）】
2G059
4B061
【Ｆターム（参考）】
2G059AA05
2G059BB11
2G059DD01
2G059DD16
2G059EE13
2G059FF04
2G059KK04
2G059MM03
2G059MM05
4B061AA01
4B061AB03
4B061AB07
4B061BA17
4B061BB20
(57)【要約】
【課題】追熟室内のバナナの追熟状態を予測する。
【解決手段】追熟室２内のバナナ３のカラー画像から色彩を判定する色彩判定部１２ａと、追熟スケジュールと追熟室２内の温度に基づいて温度累積値を算出する温度累積値算出部１１ａと、バナナ３の個体スコアの平均値を算出する平均値算出部１２ｂと、個体スコアの平均値が閾値を超えたか否かを判定する平均値判定部１２ｃと、個体スコアの平均値が閾値を超えたときに、追熟室２内の温度累積値を説明変数データ、個体スコアの平均値を目的変数データとして出力するデータ出力部１３と、ロジスティック回帰モデル：１／（１＋ｅｘｐ（－ａ（ｘ－ｂ）））において、説明変数データと目的変数データをもとに、最尤推定法によりａ，ｂを算出する追熟モデル生成部１４と、算出したａ，ｂをロジスティック回帰モデルに当てはめてバナナ３の追熟状態を予測する追熟状態予測部１５と、を備える。
【選択図】図１
【解決手段】追熟室２内のバナナ３のカラー画像から色彩を判定する色彩判定部１２ａと、追熟スケジュールと追熟室２内の温度に基づいて温度累積値を算出する温度累積値算出部１１ａと、バナナ３の個体スコアの平均値を算出する平均値算出部１２ｂと、個体スコアの平均値が閾値を超えたか否かを判定する平均値判定部１２ｃと、個体スコアの平均値が閾値を超えたときに、追熟室２内の温度累積値を説明変数データ、個体スコアの平均値を目的変数データとして出力するデータ出力部１３と、ロジスティック回帰モデル：１／（１＋ｅｘｐ（－ａ（ｘ－ｂ）））において、説明変数データと目的変数データをもとに、最尤推定法によりａ，ｂを算出する追熟モデル生成部１４と、算出したａ，ｂをロジスティック回帰モデルに当てはめてバナナ３の追熟状態を予測する追熟状態予測部１５と、を備える。
【選択図】図１
【特許請求の範囲】
【請求項1】
追熟室内で追熟されるバナナの追熟状態を予測するバナナ追熟予測装置であって、
追熟スケジュールの設定、前記追熟室内の温度累積値の初期化、バナナの色彩の正規化を含む初期設定を行う設定部と、
前記追熟室内のバナナに照明をあてて撮影するバナナ撮影用カメラと、
前記追熟室内の温度を計測する温度計と、
前記バナナ撮影用カメラにより撮影したカラー画像から前記追熟室内のバナナの色彩を判定する色彩判定部と、
前記追熟スケジュールと前記温度計により計測した前記追熟室内の温度に基づいて前記追熟室内の温度累積値を算出する温度累積値算出部と、
前記色彩判定部の判定結果に基づいて前記追熟室内のバナナの個体スコアの平均値を算出する平均値算出部と、
前記平均値算出部により算出した前記個体スコアの平均値が閾値を超えたか否かを判定する平均値判定部と、
前記個体スコアの平均値が閾値を超えたと前記平均値判定部が判定したときに、前記追熟室内の温度累積値を説明変数データ、前記個体スコアの平均値を目的変数データとして出力するデータ出力部と、
前記データ出力部が出力する前記説明変数データと前記目的変数データをもとに、最尤推定法によりロジスティック回帰モデル：ｐ＝１／（１＋ｅｘｐ（－ａ（ｘ－ｂ）））（但し、ｐ：個体スコアの平均値、ｘ：追熟室内の温度累積値、ａ：増加率、ｂ：変曲点）のａ，ｂを算出する追熟モデル生成部と、
前記個体スコアの平均値または前記追熟室内の温度累積値を設定し、前記追熟モデル生成部にて算出したａ，ｂを前記ロジスティック回帰モデルに当てはめてバナナの追熟状態を予測する追熟状態予測部と、を備えたことを特徴とするバナナ追熟予測装置。
追熟スケジュールの設定、前記追熟室内の温度累積値の初期化、バナナの色彩の正規化を含む初期設定を行う設定部と、
前記追熟室内のバナナに照明をあてて撮影するバナナ撮影用カメラと、
前記追熟室内の温度を計測する温度計と、
前記バナナ撮影用カメラにより撮影したカラー画像から前記追熟室内のバナナの色彩を判定する色彩判定部と、
前記追熟スケジュールと前記温度計により計測した前記追熟室内の温度に基づいて前記追熟室内の温度累積値を算出する温度累積値算出部と、
前記色彩判定部の判定結果に基づいて前記追熟室内のバナナの個体スコアの平均値を算出する平均値算出部と、
前記平均値算出部により算出した前記個体スコアの平均値が閾値を超えたか否かを判定する平均値判定部と、
前記個体スコアの平均値が閾値を超えたと前記平均値判定部が判定したときに、前記追熟室内の温度累積値を説明変数データ、前記個体スコアの平均値を目的変数データとして出力するデータ出力部と、
前記データ出力部が出力する前記説明変数データと前記目的変数データをもとに、最尤推定法によりロジスティック回帰モデル：ｐ＝１／（１＋ｅｘｐ（－ａ（ｘ－ｂ）））（但し、ｐ：個体スコアの平均値、ｘ：追熟室内の温度累積値、ａ：増加率、ｂ：変曲点）のａ，ｂを算出する追熟モデル生成部と、
前記個体スコアの平均値または前記追熟室内の温度累積値を設定し、前記追熟モデル生成部にて算出したａ，ｂを前記ロジスティック回帰モデルに当てはめてバナナの追熟状態を予測する追熟状態予測部と、を備えたことを特徴とするバナナ追熟予測装置。
【請求項2】
追熟室内で追熟されるバナナの追熟状態を予測するバナナ追熟予測方法であって、
追熟スケジュールの設定、前記追熟室内の温度累積値の初期化、バナナの色彩の正規化を含む初期設定を行うステップと、
前記追熟室内の温度を読み込むステップと、
前記追熟室内のバナナに照明をあてて撮影するステップと、
前記撮影したカラー画像から前記追熟室内のバナナの色彩を判定するステップと、
前記追熟スケジュールと前記追熟室内の温度に基づいて前記追熟室内の温度累積値を算出するステップと、
前記追熟室内のバナナの色彩の判定結果に基づいて当該バナナの個体スコアの平均値を算出するステップと、
前記個体スコアの平均値が閾値を超えたか否かを判定するステップと、
前記個体スコアの平均値が前記閾値を超えたと判定したときに、前記追熟室内の温度累積値を説明変数データ、前記個体スコアの平均値を目的変数データとして出力するステップと、
前記説明変数データと前記目的変数データをもとに、最尤推定法によりロジスティック回帰モデル：ｐ＝１／（１＋ｅｘｐ（－ａ（ｘ－ｂ）））（但し、ｐ：個体スコアの平均値、ｘ：追熟室内の温度累積値、ａ：増加率、ｂ：変曲点）のａ，ｂを算出するステップと、
前記個体スコアの平均値または前記追熟室内の温度累積値を設定し、前記算出したａ，ｂを前記ロジスティック回帰モデルに当てはめてバナナの追熟状態を予測するステップと、を含むことを特徴とするバナナ追熟予測方法。
追熟スケジュールの設定、前記追熟室内の温度累積値の初期化、バナナの色彩の正規化を含む初期設定を行うステップと、
前記追熟室内の温度を読み込むステップと、
前記追熟室内のバナナに照明をあてて撮影するステップと、
前記撮影したカラー画像から前記追熟室内のバナナの色彩を判定するステップと、
前記追熟スケジュールと前記追熟室内の温度に基づいて前記追熟室内の温度累積値を算出するステップと、
前記追熟室内のバナナの色彩の判定結果に基づいて当該バナナの個体スコアの平均値を算出するステップと、
前記個体スコアの平均値が閾値を超えたか否かを判定するステップと、
前記個体スコアの平均値が前記閾値を超えたと判定したときに、前記追熟室内の温度累積値を説明変数データ、前記個体スコアの平均値を目的変数データとして出力するステップと、
前記説明変数データと前記目的変数データをもとに、最尤推定法によりロジスティック回帰モデル：ｐ＝１／（１＋ｅｘｐ（－ａ（ｘ－ｂ）））（但し、ｐ：個体スコアの平均値、ｘ：追熟室内の温度累積値、ａ：増加率、ｂ：変曲点）のａ，ｂを算出するステップと、
前記個体スコアの平均値または前記追熟室内の温度累積値を設定し、前記算出したａ，ｂを前記ロジスティック回帰モデルに当てはめてバナナの追熟状態を予測するステップと、を含むことを特徴とするバナナ追熟予測方法。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、追熟室内に保管されるバナナの追熟状態を予測するバナナ追熟予測装置およびバナナ追熟予測方法に関するものである。
【背景技術】
【0002】
バナナは、収穫されると、追熟室内に保管され、１週間～１０日かけて糖度を上げた後に出荷される。
【0003】
さらに説明すると、産地で収穫されたバナナは、下記特許文献１の従来の技術の欄にも記載されているように、そのままの状態で市場に出荷されるのではなく、収穫時に未熟成の青いバナナを店頭市場に出す前に追熟室内に保管し、追熟室内の温度を制御することで追熟状態を調整し、適度に追熟されたときに出荷される。
【0004】
ところで、下記特許文献１では、追熟室内におけるバナナの追熟状態を知るため、換気可能手段を有し、多数のバナナ果実が収納される熟成加工室内に、バナナの果皮に接触する測定端部を配置し、加工室外の加工室に近接した位置に、クロロフィル又はカロチノイドに吸収される波長を有する単色光を生成する単色光生成手段および光電変換手段を備えた光度計を配置し、測定端部と単色光生成手段および光電変換手段との間を光伝達手段で接続し、加工室から離れた場所にある管理室内に配置した表示手段に光度計で得られた測定情報に基づくバナナ果皮色相状態を表示するように構成している。
【先行技術文献】
【特許文献】
【0005】
【特許文献1】特開平７－１２８２２６号公報
【発明の概要】
【発明が解決しようとする課題】
【0006】
しかしながら、上述した特許文献１に開示される従来の方法では、バナナの果皮色を直接測定して表示することにより、バナナの果皮色から追熟状態を外部から知ることができるが、バナナの追熟状態を予測することはできなかった。
【0007】
そこで、本発明は上記問題点に鑑みてなされたものであって、追熟室内のバナナの追熟状態を予測することができるバナナ追熟予測装置およびバナナ追熟予測方法を提供することを目的としている。
【課題を解決するための手段】
【0008】
上記目的を達成するため、本発明の請求項１に記載されたバナナ追熟予測装置は、追熟室内で追熟されるバナナの追熟状態を予測するバナナ追熟予測装置であって、
追熟スケジュールの設定、前記追熟室内の温度累積値の初期化、バナナの色彩の正規化を含む初期設定を行う設定部と、
前記追熟室内のバナナに照明をあてて撮影するバナナ撮影用カメラと、
前記追熟室内の温度を計測する温度計と、
前記バナナ撮影用カメラにより撮影したカラー画像から前記追熟室内のバナナの色彩を判定する色彩判定部と、
前記追熟スケジュールと前記温度計により計測した前記追熟室内の温度に基づいて前記追熟室内の温度累積値を算出する温度累積値算出部と、
前記色彩判定部の判定結果に基づいて前記追熟室内のバナナの個体スコアの平均値を算出する平均値算出部と、
前記平均値算出部により算出した前記個体スコアの平均値が閾値を超えたか否かを判定する平均値判定部と、
前記個体スコアの平均値が閾値を超えたと前記平均値判定部が判定したときに、前記追熟室内の温度累積値を説明変数データ、前記個体スコアの平均値を目的変数データとして出力するデータ出力部と、
前記データ出力部が出力する前記説明変数データと前記目的変数データをもとに、最尤推定法によりロジスティック回帰モデル：ｐ＝１／（１＋ｅｘｐ（－ａ（ｘ－ｂ）））（但し、ｐ：個体スコアの平均値、ｘ：追熟室内の温度累積値、ａ：増加率、ｂ：変曲点）のａ，ｂを算出する追熟モデル生成部と、
前記個体スコアの平均値または前記追熟室内の温度累積値を設定し、前記追熟モデル生成部にて算出したａ，ｂを前記ロジスティック回帰モデルに当てはめてバナナの追熟状態を予測する追熟状態予測部と、を備えたことを特徴とする。
追熟スケジュールの設定、前記追熟室内の温度累積値の初期化、バナナの色彩の正規化を含む初期設定を行う設定部と、
前記追熟室内のバナナに照明をあてて撮影するバナナ撮影用カメラと、
前記追熟室内の温度を計測する温度計と、
前記バナナ撮影用カメラにより撮影したカラー画像から前記追熟室内のバナナの色彩を判定する色彩判定部と、
前記追熟スケジュールと前記温度計により計測した前記追熟室内の温度に基づいて前記追熟室内の温度累積値を算出する温度累積値算出部と、
前記色彩判定部の判定結果に基づいて前記追熟室内のバナナの個体スコアの平均値を算出する平均値算出部と、
前記平均値算出部により算出した前記個体スコアの平均値が閾値を超えたか否かを判定する平均値判定部と、
前記個体スコアの平均値が閾値を超えたと前記平均値判定部が判定したときに、前記追熟室内の温度累積値を説明変数データ、前記個体スコアの平均値を目的変数データとして出力するデータ出力部と、
前記データ出力部が出力する前記説明変数データと前記目的変数データをもとに、最尤推定法によりロジスティック回帰モデル：ｐ＝１／（１＋ｅｘｐ（－ａ（ｘ－ｂ）））（但し、ｐ：個体スコアの平均値、ｘ：追熟室内の温度累積値、ａ：増加率、ｂ：変曲点）のａ，ｂを算出する追熟モデル生成部と、
前記個体スコアの平均値または前記追熟室内の温度累積値を設定し、前記追熟モデル生成部にて算出したａ，ｂを前記ロジスティック回帰モデルに当てはめてバナナの追熟状態を予測する追熟状態予測部と、を備えたことを特徴とする。
【0009】
本発明の請求項２に記載されたバナナ追熟予測方法は、追熟室内で追熟されるバナナの追熟状態を予測するバナナ追熟予測方法であって、
追熟スケジュールの設定、前記追熟室内の温度累積値の初期化、バナナの色彩の正規化を含む初期設定を行うステップと、
前記追熟室内の温度を読み込むステップと、
前記追熟室内のバナナに照明をあてて撮影するステップと、
前記撮影したカラー画像から前記追熟室内のバナナの色彩を判定するステップと、
前記追熟スケジュールと前記追熟室内の温度に基づいて前記追熟室内の温度累積値を算出するステップと、
前記追熟室内のバナナの色彩の判定結果に基づいて当該バナナの個体スコアの平均値を算出するステップと、
前記個体スコアの平均値が閾値を超えたか否かを判定するステップと、
前記個体スコアの平均値が前記閾値を超えたと判定したときに、前記追熟室内の温度累積値を説明変数データ、前記個体スコアの平均値を目的変数データとして出力するステップと、
前記説明変数データと前記目的変数データをもとに、最尤推定法によりロジスティック回帰モデル：ｐ＝１／（１＋ｅｘｐ（－ａ（ｘ－ｂ）））（但し、ｐ：個体スコアの平均値、ｘ：追熟室内の温度累積値、ａ：増加率、ｂ：変曲点）のａ，ｂを算出するステップと、
前記個体スコアの平均値または前記追熟室内の温度累積値を設定し、前記算出したａ，ｂを前記ロジスティック回帰モデルに当てはめてバナナの追熟状態を予測するステップと、を含むことを特徴とする。
追熟スケジュールの設定、前記追熟室内の温度累積値の初期化、バナナの色彩の正規化を含む初期設定を行うステップと、
前記追熟室内の温度を読み込むステップと、
前記追熟室内のバナナに照明をあてて撮影するステップと、
前記撮影したカラー画像から前記追熟室内のバナナの色彩を判定するステップと、
前記追熟スケジュールと前記追熟室内の温度に基づいて前記追熟室内の温度累積値を算出するステップと、
前記追熟室内のバナナの色彩の判定結果に基づいて当該バナナの個体スコアの平均値を算出するステップと、
前記個体スコアの平均値が閾値を超えたか否かを判定するステップと、
前記個体スコアの平均値が前記閾値を超えたと判定したときに、前記追熟室内の温度累積値を説明変数データ、前記個体スコアの平均値を目的変数データとして出力するステップと、
前記説明変数データと前記目的変数データをもとに、最尤推定法によりロジスティック回帰モデル：ｐ＝１／（１＋ｅｘｐ（－ａ（ｘ－ｂ）））（但し、ｐ：個体スコアの平均値、ｘ：追熟室内の温度累積値、ａ：増加率、ｂ：変曲点）のａ，ｂを算出するステップと、
前記個体スコアの平均値または前記追熟室内の温度累積値を設定し、前記算出したａ，ｂを前記ロジスティック回帰モデルに当てはめてバナナの追熟状態を予測するステップと、を含むことを特徴とする。
【発明の効果】
【0010】
本発明によれば、追熟期間における前半のバナナの追熟状態のデータを用いて後半のバナナの追熟状態を予測することができる。これにより、バナナの追熟状態の予測を元にして、出荷先相手が所望する追熟状態のバナナの出荷を行うことが可能となる。
【図面の簡単な説明】
【0011】
【図1】本発明に係るバナナ追熟予測装置の概略構成を示すブロック図である。
【図2】本発明に係るバナナ追熟予測方法の手順を示すフローチャートである。
【図3】追熟期間全てのデータを用いたときの追熟室内の温度累積値と追熟室内のバナナの個体スコアの平均値の関係の具体例を示す図である。
【発明を実施するための形態】
【0012】
以下、本発明を実施するための形態について、添付した図面を参照しながら詳細に説明する。
【0013】
図１に示すように、本実施の形態のバナナ追熟予測装置１は、収穫後に追熟室２内に保管されるバナナ３の追熟状態を予測するもので、追熟室２内の環境を計測する各種の測定計４、追熟室２内に保管されるバナナ３に照明を当てる照明器５、照明器５で照らされた追熟室２内のバナナ３を撮影するバナナ撮影用カメラ６、設定部７、中央演算処理装置８、表示部９を備えて概略構成される。
【0014】
測定計４は、追熟室２内環境の物理量を所定時間（例えば１時間）ごとに計測するもので、温度計４ａ、湿度計４ｂ、ＣＯ2 計４ｃ、エチレン濃度計４ｄから構成される。
【0015】
温度計４ａは、追熟室２内のバナナ３近傍に設けられ、バナナ３を追熟する上で最も重要な物理量である追熟室２内の温度を計測する。湿度計４ｂは、追熟室２内の湿度を計測する。ＣＯ2 計４ｃは、追熟室２内のＣＯ2 濃度を計測する。エチレン濃度計４ｄは、追熟室２内のエチレン濃度を計測する。これら所定時間（例えば１時間）ごとに計測された情報（追熟室２内の温度、湿度、ＣＯ2 濃度、エチレン濃度）は中央演算処理装置８に逐次送られる。
【0016】
設定部７は、初期設定を含め、追熟室２内に保管されたバナナ３を追熟するために必要な各種設定を行う。初期設定としては、追熟スケジュールの設定、追熟室２内の温度累積値の初期化、バナナ３の色彩の正規化がある。
【0017】
追熟スケジュールは、追熟期間（例えば７日～１０日）、エチレンガス注入スケジュール（例えば５００ｐｐｍ以上の濃度で最初の２４時間）、温度スケジュール（例えば２０℃（２４時間）→１８℃（２４時間）→以後１６℃、２２℃（２４時間）→２０℃（２４時間）→以後１８℃など）からなる。
【0018】
追熟室２内の温度累積値とは、バナナ３に与える熱量、すなわち、温度（追熟室２内の温度）×時間（温度スケジュールの時間）に相当するものである。
【0019】
バナナ３の個体スコアは、バナナ３の色彩を正規化することで決定される。具体的には、追熟に伴って変化するバナナ３の色彩を複数段階のカラーコードに分け、カラーコードを０．０～１．０に数値化する。カラーコードの例としては、ドール（登録商標）のバナナカラーチャートが挙げられる。本実施の形態では、バナナ３の色彩をオールグリーン（カラーコード１）、ライトグリーン（カラーコード２）、ハーフグリーン（カラーコード３）、ハーフイエロー（カラーコード４）、グリーンチップ（カラーコード５）、フルイエロー（カラーコード６）の６段階に分け、カラーコード１→０．０、カラーコード２→０．２、カラーコード３→０．４、カラーコード４→０．６、カラーコード５→０．８、カラーコード６→１．０に数値化している。つまり、個体スコアが０では、バナナが未成熟であり、個体スコアが１に近づくほど、バナナの追熟が進んでいることを示す。
【0020】
中央演算処理装置８は、追熟室２内のバナナ３の追熟状態を予測する際に各部を統括制御するもので、追熟室環境計測部１１、追熟状態計測・制御部１２、データ出力部１３、追熟モデル生成部１４、追熟状態予測部１５を含んで構成される。
【0021】
追熟室環境計測部１１は、計測計４にて計測した情報（追熟室２内の温度、湿度、ＣＯ2 濃度、エチレン濃度）を取得して追熟室２内環境の物理量を具体的数値として計測しており、温度累積値算出部１１ａを備える。
【0022】
温度累積値算出部１１ａは、初期設定された追熟スケジュールの温度スケジュールと計測計４の温度計４ａにより計測した追熟室２内の温度に基づいて追熟室２内の温度累積値（追熟室２内の温度×温度スケジュールの時間）を算出する。
【0023】
追熟状態計測・制御部１２は、照明器５の点灯／消灯制御、バナナ撮影用カメラ６により撮影したバナナ３のカラー画像の取得を行っており、色彩判定部１２ａ、平均値算出部１２ｂ、平均値判定部１２ｃを備える。本実施の形態では、照明器５は白色光を発光する照明を用いた。
【0024】
色彩判定部１２ａは、照明器５により照らされたバナナ３をバナナ撮影用カメラ６により撮影したバナナ３のカラー画像から追熟室２内のバナナ３の色彩をバナナ３ごとに判定する。具体的には、追熟室２内に置いた、バナナ３とカラーコードチャートを同時に撮影し、撮影したバナナ３のカラー画像とカラーコードチャートのカラー画像を比較して、追熟室２内のバナナ３の色彩がカラーコードチャートのオールグリーン、ライトグリーン、ハーフグリーン、ハーフイエロー、グリーンチップ、フルイエローの何れに該当するかをバナナ３ごとに判定する。本実施の形態では、追熟室２内でバナナ３とカラーコードチャートを同時に撮影することで、照明器５による発色の影響をなくしている。
【0025】
平均値算出部１２ｂは、色彩判定部１２ａの判定結果に基づいて追熟室２内のバナナ３の個体スコアの平均値（０．０～１．０）を算出する。
【0026】
平均値判定部１２ｃは、平均値算出部１２ｂにより算出したバナナ３の個体スコアの平均値が予め設定される閾値（例えば０．５）を超えたか否かを判定する。
【0027】
データ出力部１３は、バナナ３の個体スコアの平均値が閾値を超えたと平均値判定部１２ｃが判定したときに、追熟室２内の温度累積値を説明変数データ、バナナ３の個体スコアの平均値を目的変数データとしてデータを出力する。
【0028】
追熟モデル生成部１４は、ロジスティック回帰モデル：ｐ＝１／（１＋ｅｘｐ（－ａ（ｘ－ｂ）））（但し、ｐ：バナナ３の個体スコアの平均値、ｘ：追熟室２内の温度累積値、ａ：増加率、ｂ：変曲点）において、データ出力部１３が出力する説明変数データと目的変数データをもとに、最尤推定法によりａ（増加率），ｂ（変曲点）を算出する。
【0029】
追熟状態予測部１５は、ロジスティック回帰モデルに対し、追熟モデル生成部１４にて算出したａ（増加率），ｂ（変曲点）を当てはめてバナナ３の追熟状態を予測する。具体的には、追熟モデル生成部１４にて算出したａ（増加率），ｂ（変曲点）をロジスティック回帰モデルに当てはめ、バナナ３の個体スコアの平均値ｐを設定し、設定したバナナ３の個体スコアの平均値ｐが得られるまでに要する追熟室２内の温度累積値ｘを予測する。または、追熟モデル生成部１４にて算出したａ（増加率），ｂ（変曲点）をロジスティック回帰モデルに当てはめ、追熟室２内の温度累積値ｘを設定し、設定した追熟室２内の温度累積値ｘで得られるバナナ３の個体スコアの平均値ｐを予測する。
【0030】
表示部９は、例えば液晶表示器などの各種表示器で構成され、設定画面の表示、追熟スケジュールの表示、計測計４が計測した情報（追熟室２内の温度、湿度、ＣＯ2 濃度、エチレン濃度）の表示、予測結果の表示などを行う。
【0031】
次に、上記のように構成されるバナナ追熟予測装置１を用いたバナナ３の追熟状態の予測方法について図２のフローチャートを参照しながら説明する。
【0032】
収穫した所定数のバナナ３を追熟室２に保管し、初期設定を行う（ＳＴ１）。初期設定では、追熟スケジュールの設定、追熟室２内の温度累積値の初期化、バナナ３の色彩の正規化を行う。
【0033】
初期設定を終え、追熟スケジュールに従ってバナナ３の追熟が開始されると、計測計４の温度計４ａにより追熟室２内の温度を所定時間（例えば１時間）ごとに読み込む（ＳＴ２）。
【0034】
続いて、追熟状態計測・制御部１２により照明器５を点灯してバナナ３を照らし、バナナ撮影用カメラ６によりバナナ３を所定時間（例えば６時間）ごとに撮影する。撮影したカラー画像は追熟状態計測・制御部１２の色彩判定部１２ａに逐次送られる。そして、色彩判定部１２ａは、バナナ撮影用カメラ６により撮影したカラー画像から追熟室２内のバナナ３の色彩をバナナ３ごとに判定する（ＳＴ３）。
【0035】
また、追熟室環境計測部１１の温度累積値算出部１１ａは、追熟室２内の温度累積値（追熟室２内の温度×温度スケジュールの時間）を算出する（ＳＴ４）。
【0036】
そして、追熟状態計測・制御部１２の平均値算出部１２ｂは、色彩判定部１２ａの判定結果に基づいて追熟室２内のバナナ３の個体スコアの平均値を算出する（ＳＴ５）。
【0037】
続いて、追熟状態計測・制御部１２の平均値判定部１２ｃは、平均値算出部１２ｂにより算出したバナナ３の個体スコアの平均値が閾値（０．５）以上か否かを判別する（ＳＴ６）。
【0038】
そして、バナナ３の個体スコアの平均値が閾値（０．５）以上であると判別すると（ＳＴ６－Ｙｅｓ）、データ出力部１３が追熟室２内の温度累積値を説明変数データ、バナナ３の個体スコアの平均値を目的変数データとしてデータを出力する（ＳＴ７）。
【0039】
なお、バナナ３の個体スコアの平均値が閾値（０．５）以上でないと判別すると（ＳＴ６－Ｎｏ）、ＳＴ２に戻り、ＳＴ２～ＳＴ５の処理を繰り返し実行する。
【0040】
続いて、追熟モデル生成部１４は、ロジスティック回帰モデル：ｐ＝１／（１＋ｅｘｐ（－ａ（ｘ－ｂ）））において、データ出力部１３が出力する説明変数データと目的変数データをもとに、最尤推定法によりａ（増加率）とｂ（変曲点）を算出する（ＳＴ８）。
【0041】
なお、最尤推定法は、与えられたデータ（追熟室２内の温度累積値による説明変数データ、バナナ３の個体スコアの平均値による目的変数データ）からそれが従う確率分布の母数（ａ（増加率）、ｂ（変曲点））を点推定する周知の方法である。
【0042】
そして、追熟状態予測部１５は、ロジスティック回帰モデルにａ（増加率），ｂ（変曲点）を当てはめ、追熟状態を予測する（ＳＴ９）。この追熟状態の予測では、ロジスティック回帰モデルにａ，ｂを当てはめ、バナナ３の個体スコアの平均値ｐまたは追熟室２内の温度累積値ｘを決めることにより、バナナ３の個体スコアｐの平均値が得られるまでに要する追熟室２内の温度累積値ｘを予測したり、追熟室２内の温度累積値ｘで得られるバナナ３の個体スコアの平均値ｐを予測することができる。すなわち、追熟期間の途中までのデータを用いた処理により、所望とするカラーコードのバナナ３に追熟させるまでに要する時間や、追熟室２内の温度累積値によって得られるバナナ３のカラーコードを予測することができる。
【0043】
［追熟状態の予測結果の具体例］
追熟期間：７日間、エチレンガスの注入スケジュール：所定濃度で最初の２４時間のみ、温度スケジュール：２２℃（２４時間）→２０℃（２４時間）→１８℃（５日間）として追熟スケジュールを設定し、６時間おきに追熟室２内のバナナ３を撮影して追熟状態の調査を行った。
追熟期間：７日間、エチレンガスの注入スケジュール：所定濃度で最初の２４時間のみ、温度スケジュール：２２℃（２４時間）→２０℃（２４時間）→１８℃（５日間）として追熟スケジュールを設定し、６時間おきに追熟室２内のバナナ３を撮影して追熟状態の調査を行った。
【0044】
そして、追熟室２内の温度累積値（追熟室２内の温度×温度スケジュールの時間）と追熟室２内のバナナ３の個体スコアの平均値の関係を調べたところ、追熟期間全て（７日間）のデータを用いた場合は図３に示す結果が得られ、図２のフローチャートに基づくバナナ追熟予測方法により追熟期間における前半（４日間）のデータを用いて後半（３日間）を予測した場合は図４に示す結果が得られた。
【0045】
両者を比較すると、ロジスティック回帰モデルのａ（増加率）に関しては、図３の追熟期間全てのデータを用いた場合：ａ［１０-3］＝１．８２５に対し、図４の追熟期間の前半のデータを用いて予測した場合：ａ［１０-3］＝１．５０２となり、両者で多少のずれはあるものの特徴を捉えた値の予測結果が得られた。また、ロジスティック回帰モデルのｂ（変曲点）に関しては、図３の追熟期間全てのデータを用いた場合：ｂ［１０3 ］＝０．９０９に対し、図４の追熟期間の前半のデータを用いて予測した場合：ｂ［１０3 ］＝０．９６１となり、両者でほぼ一致するという結果が得られた。
【0046】
このように、上述した実施の形態によれば、追熟期間における前半のバナナの追熟状態のデータを用いて後半のバナナの追熟状態を予測することができる。これにより、バナナの追熟状態の予測を元にして、出荷先相手が所望する追熟状態のバナナの出荷を行うことが可能となる。
【0047】
以上、本発明に係るバナナ追熟予測装置およびバナナ追熟予測方法の最良の形態について説明したが、この形態による記述及び図面により本発明が限定されることはない。すなわち、この形態に基づいて当業者等によりなされる他の形態、実施例及び運用技術などはすべて本発明の範疇に含まれることは勿論である。
【符号の説明】
【0048】
１ バナナ追熟予測装置
２ 追熟室
３ バナナ
４ 測定計
４ａ 温度計
４ｂ 湿度計
４ｃ ＣＯ2 計
４ｄ エチレン濃度計
５ 照明器
６ バナナ撮影用カメラ
７ 設定部
８ 中央演算処理部
９ 表示部
１１ 追熟室環境計測部
１１ａ 温度累積値算出部
１２ 追熟状態計測・制御部
１２ａ 色彩判定部
１２ｂ 平均値算出部
１２ｃ 平均値判定部
１３ データ出力部
１４ 追熟モデル生成部
１５ 追熟状態予測部
２ 追熟室
３ バナナ
４ 測定計
４ａ 温度計
４ｂ 湿度計
４ｃ ＣＯ2 計
４ｄ エチレン濃度計
５ 照明器
６ バナナ撮影用カメラ
７ 設定部
８ 中央演算処理部
９ 表示部
１１ 追熟室環境計測部
１１ａ 温度累積値算出部
１２ 追熟状態計測・制御部
１２ａ 色彩判定部
１２ｂ 平均値算出部
１２ｃ 平均値判定部
１３ データ出力部
１４ 追熟モデル生成部
１５ 追熟状態予測部
【図1】
【図2】
【図3】
【図4】