(57)【要約】【課題】植物工場などにおいて植物の生育（成長速度）にバラツキがある場合でも、収穫される植物の重量を所望の範囲内（規格内）に制御することのできる植物の生産方法を提供すること。【解決手段】少なくとも１つの植物からなる植物群の種子を播種し、発芽させて成長させることにより稚苗を得る、緑化工程と、稚苗を育苗室に移植して成育させることにより、定植苗を得る、育苗工程と、定植苗を栽培室に移植して成育させた後に収穫する本栽培工程と、を含む、植物の生産方法。植物群を構成する個々の植物について、育苗工程の複数の時点における重量のデータを収集し、収取された重量のデータから、緑化期間、育苗期間および本栽培期間の合計期間である総栽培期間と、植物群の重量と、の関係式を作成し、関係式に基づいて、植物群について、収穫される植物のうち所定範囲内の重量を有する植物の比率が最適化されるように、総栽培期間が調整される。【選択図】図２
2019041582.tif 000002
【特許請求の範囲】
【請求項１】
少なくとも１つの植物からなる植物群の種子を播種し、発芽させて成長させることにより稚苗を得る、緑化工程と、
前記稚苗を育苗室に移植して成育させることにより、定植苗を得る、育苗工程と、
前記定植苗を栽培室に移植して成育させた後に収穫する本栽培工程と、を含む、植物の生産方法であって、
前記植物群を構成する個々の植物について、前記育苗工程の複数の時点における重量のデータを収集し、
収集された前記重量のデータから、前記緑化工程を実施する期間である緑化期間、前記育苗工程を実施する期間である育苗期間、および、前記本栽培工程を実施する期間である本栽培期間の合計期間である、総栽培期間と、前記植物群の重量と、の関係式を作成し、
前記関係式に基づいて、前記植物群について、収穫される植物のうち所定範囲内の重量を有する植物の比率が最適化されるように、前記総栽培期間が調整される、植物の生産方法。
【請求項２】
前記関係式は、下記式（１）を前記重量のデータにフィッティングさせて、下記式（１）中の成長に関する係数Ｋ、Ｗ１を決定することにより得られる式である、請求項１に記載の植物の生産方法。
2019041582.tif 000012
〔式（１）において、
ｔは、前記総栽培期間（日）であり、
ｔ１は、前記緑化期間（日）であり、
Ｗ（ｔ）は、前記総栽培期間ｔ（日）経過時における前記植物群の重量の平均値であり、
Ｗｍａｘは、前記植物の想定される最大重量（ｇ）であり、
Ｗ１は、前記植物の前記稚苗の重量（ｇ）であり、
Ｋは、成長率である。〕
【請求項３】
前記重量のデータは、葉面積の測定値に基づく推定により得られる、請求項１または２に記載の植物の生産方法。
【請求項４】
前記植物群のうち袋売りされる植物の重量の最小値である閾値重量以上の植物の重量当たりの価格と、前記植物群のうち前記閾値重量未満の植物の重量当たりの価格と、の比率に基づいて、前記植物群の生産および販売によって得られる利益率が最大化するように、
前記植物群について、収穫される植物のうち所定範囲内の重量を有する植物の比率が最適化される、請求項１〜３のいずれか１項に記載の植物の生産方法。
【請求項５】
前記本栽培期間を調整することにより、前記総栽培期間を調整する、請求項１〜４のいずれか１項に記載の植物の生産方法。
【請求項６】
前記本栽培期間を調整することにより、前記総栽培期間を調整し、
下記式（２）に基づいて、前記利益率が最大化するように、前記本栽培期間が調整される、請求項４に記載の植物の生産方法。
2019041582.tif 000013
〔式（２）において、
Ｙ（ｔ，ｗ＊，λ）は、植物の１株当たりの利益率（／株）であり、
ｗ＊は、前記閾値重量（ｇ）であり、
ｒ＊は、重量がｗ＊以上の株（袋売りの株）の割合（／株）であり、
λは、袋売り商品の重量当たりの価格に対する重量売り商品の重量当たりの価格の割合であり、
<Ｗ>は、重量がｗ＊未満の株（重量売りの株）の平均重量（ｇ）であり、
Ｐは、袋売り商品を消費者が購入する割合である購買率であり、
αは、固定費（／株）であり、
βは、１日当りの１株ごとの前記緑化工程および前記育苗工程における光熱費（／（日・株））であり、
γは、１日当りの１株ごとの本栽培工程における光熱費（／（日・株））であり、
ｔ１は、前記緑化期間（日）であり、
ｔ２は、前記育苗期間（日）である。〕
【請求項７】
前記育苗期間を調整することにより、前記総栽培期間を調整する、請求項１〜４のいずれか１項に記載の植物の生産方法。
【請求項８】
前記育苗期間を調整することにより、前記総栽培期間を調整し、
下記式（３）に基づいて、前記利益率が最大化するように、前記本栽培期間が調整される、請求項４に記載の植物の生産方法。
2019041582.tif 000014
〔式（３）において、
Ｙ（ｔ，ｗ＊，λ）は、植物の１株当たりの利益率（／株）であり、
ｗ＊は、前記閾値重量（ｇ）であり、
ｒ＊は、重量がｗ＊以上の株（袋売りの株）の割合（／株）
λは、袋売り商品の重量当たりの価格に対する重量売り商品の重量当たりの価格の割合であり、
<Ｗ>は、重量がｗ＊未満の株（重量売りの株）の平均重量（ｇ）であり、
Ｐは、袋売り商品を消費者が購入する割合である購買率であり、
αは、固定費（／株）であり、
βは、１日当りの１株ごとの前記緑化工程および前記育苗工程における光熱費（／（日・株））であり、
γは、１日当りの１株ごとの本栽培工程における光熱費（／（日・株））であり、
ｔ３は、前記本栽培期間（日）である。〕
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、植物の生産方法に関する。
【背景技術】
【０００２】
植物工場（例えば、人工光型植物工場、ならびに太陽光利用型植物工場）は、農業の工業化ならびに企業経営を実現する一つの柱に位置付けられ、輸出産業としての期待も大きいため、植物工場の技術開発が産業的に注目を集めている。
【０００３】
しかし、植物工場での植物の生産では、植物の生育が不均一であり（成育バラツキがあり）、生産が安定しない（生育バラツキにより収益が安定しない）という問題があった。例えば、野菜の販売形態には、重量売り(規格外のもの)と袋売り(規格内のものを個別に袋売り)がある。植物の生育バラツキが大きいと、規格外のものが増えて、収益性の低い重量売りの比率が増えるため、植物工場の収益性は低下してしまう。
【０００４】
収穫量、収穫時期などを予測するシステムまたはソフトウエアは、既にいくつか提案されている（例えば、特許文献１（国際公開第２０１３／０８８５３９号）参照）。しかし、これらにおいては、集団平均または時間平均に基づく第一近似・準静的モデルを用いており、個体別の動的特性を十分に分析していない（個々の植物の重量測定などは行っていない）ため、生産の不安定性（突如、生育が不良になる現象）が十分に分析されていない。したがって、これらの従来のシステムでは、植物工場における生産不安定性などを十分に回避することはできないと考えられる。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】国際公開第２０１３／０８８５３９号
【発明の概要】
【発明が解決しようとする課題】
【０００６】
本発明の課題は、植物工場などにおいて植物の生育（成長速度）にバラツキがある場合でも、収穫される植物の重量を所望の範囲内（規格内）に制御することのできる植物の生産方法を提供することである。
【課題を解決するための手段】
【０００７】
［１］  少なくとも１つの植物からなる植物群の種子を播種し、発芽させて成長させることにより稚苗を得る、緑化工程と、
前記稚苗を育苗室に移植して成育させることにより、定植苗を得る、育苗工程と、
前記定植苗を栽培室に移植して成育させた後に収穫する本栽培工程と、を含む、植物の生産方法であって、
前記植物群を構成する個々の植物について、前記育苗工程の複数の時点における重量のデータを収集し、
収集された前記重量のデータから、前記緑化工程を実施する期間である緑化期間、前記育苗工程を実施する期間である育苗期間、および、前記本栽培工程を実施する期間である本栽培期間の合計期間である、総栽培期間と、前記植物群の重量と、の関係式を作成し、
前記関係式に基づいて、前記植物群について、収穫される植物のうち所定範囲内の重量を有する植物の比率が最適化されるように、前記総栽培期間が調整される、植物の生産方法。
【０００８】
［２］  前記関係式は、下記式（１）を前記重量のデータにフィッティングさせて、下記式（１）中の成長に関する係数Ｋ、Ｗ１を決定することにより得られる式である、［１］に記載の植物の生産方法。
【０００９】
【数１】
2019041582.tif 000003
【００１０】
〔式（１）において、
ｔは、前記総栽培期間（日）であり、
ｔ１は、前記緑化期間（日）であり、
Ｗ（ｔ）は、前記総栽培期間ｔ（日）経過時における前記植物群の重量の平均値であり、
Ｗｍａｘは、前記植物の想定される最大重量（ｇ）であり、
Ｗ１は、前記植物の前記稚苗の重量（ｇ）であり、
Ｋは、成長率である。〕
【００１１】
［３］  前記重量のデータは、葉面積の測定値に基づく推定により得られる、［１］または［２］に記載の植物の生産方法。
【００１２】
［４］  前記植物群のうち袋売りされる植物の重量の最小値である閾値重量以上の植物の重量当たりの価格と、前記植物群のうち前記閾値重量未満の植物の重量当たりの価格と、の比率に基づいて、前記植物群の生産および販売によって得られる利益率が最大化するように、
前記植物群について、収穫される植物のうち所定範囲内の重量を有する植物の比率が最適化される、［１］〜［３］のいずれかに記載の植物の生産方法。
【００１３】
［５］  前記本栽培期間を調整することにより、前記総栽培期間を調整する、［１］〜［４］のいずれかに記載の植物の生産方法。
【００１４】
［６］  前記本栽培期間を調整することにより、前記総栽培期間を調整し、
下記式（２）に基づいて、前記利益率が最大化するように、前記本栽培期間が調整される、［４］に記載の植物の生産方法。
【００１５】
【数２】
2019041582.tif 000004
【００１６】
〔式（２）において、
Ｙ（ｔ，ｗ＊，λ）は、植物の１株当たりの利益率（／株）であり、
ｗ＊は、前記閾値重量（ｇ）であり、
ｒ＊は、重量がｗ＊以上の株（袋売りの株）の割合（／株）
λは、袋売り商品の重量当たりの価格に対する重量売り商品の重量当たりの価格の割合であり、
<Ｗ>は、重量がｗ＊未満の株（重量売りの株）の平均重量（ｇ）であり、
Ｐは、袋売り商品を消費者が購入する割合である購買率であり、
αは、固定費（／株）であり、
βは、１日当りの１株ごとの前記緑化工程および前記育苗工程における光熱費（／（日・株））であり、
γは、１日当りの１株ごとの本栽培工程における光熱費（／（日・株））であり、
ｔ１は、前記緑化期間（日）であり、
ｔ２は、前記育苗期間（日）である。〕
【００１７】
［７］  前記育苗期間を調整することにより、前記総栽培期間を調整する、［１］〜［４］のいずれかに記載の植物の生産方法。
【００１８】
［８］  前記育苗期間を調整することにより、前記総栽培期間を調整し、
下記式（３）に基づいて、前記利益率が最大化するように、前記本栽培期間が調整される、［４］に記載の植物の生産方法。
【００１９】
【数３】
2019041582.tif 000005
【００２０】
〔式（３）において、
Ｙ（ｔ，ｗ＊，λ）は、植物の１株当たりの利益率（／株）であり、
ｗ＊は、前記閾値重量（ｇ）であり、
ｒ＊は、重量がｗ＊以上の株（袋売りの株）の割合（／株）
λは、袋売り商品の重量当たりの価格に対する重量売り商品の重量当たりの価格の割合であり、
<Ｗ>は、重量がｗ＊未満の株（重量売りの株）の平均重量（ｇ）であり、
Ｐは、袋売り商品を消費者が購入する割合である購買率であり、
αは、固定費（／株）であり、
βは、１日当りの１株ごとの前記緑化工程および前記育苗工程における光熱費（／（日・株））であり、
γは、１日当りの１株ごとの本栽培工程における光熱費（／（日・株））であり、
ｔ３は、前記本栽培期間（日）である。〕
【発明の効果】
【００２１】
本発明の植物の生産方法によれば、植物工場などにおいて植物の生育（成長速度）にバラツキがある場合でも、収穫される植物の重量を所望の範囲内（規格内）に制御することができる。
【図面の簡単な説明】
【００２２】
【図１】本発明の植物の生産方法の各工程を示すフロー図である。
【図２】本発明の植物の生産方法に用いられるシステムの一例を示す模式図である。
【図３】本発明の植物の生産方法に用いられるシステムの一例を説明するための模式図である。
【図４】植物の閾値重量ｗ＊と購買率Ｐとの関係を示すグラフである。
【図５】実施形態１における本栽培期間の調整を説明するためのグラフである。（ａ）は、主に育苗期間を説明するためのグラフであり、（ｂ）は、主に本栽培期間を説明するためのグラフである。
【図６】実施形態２における育苗期間の調整を説明するためのグラフである。（ａ）は、主に育苗期間を説明するためのグラフであり、（ｂ）は、主に本栽培期間を説明するためのグラフである。
【図７】（ａ）は、葉面積と重量との相関性を示すグラフである。（ｂ）は、画像処理による葉面積の測定を説明するための写真である。
【図８】育苗工程における葉面積の測定を説明するためのグラフである。（ａ）は、育苗期間２４０時間までのグラフであり、（ｂ）は、育苗期間３３６時間までのグラフである。
【図９】実施形態１における閾値重量ｗ＊および総栽培期間ｔと利益率Ｙとの関係を示すグラフである。
【図１０】従来の植物の生産方法を説明するための模式図である。
【図１１】実施形態１の植物の生産方法を説明するための模式図である。
【図１２】実施形態２における閾値重量ｗ＊および総栽培期間ｔと利益率Ｙとの関係を示すグラフである。
【図１３】従来の植物の生産方法を説明するための模式図である。
【図１４】実施形態２の植物の生産方法を説明するための模式図である。
【発明を実施するための形態】
【００２３】
（植物の生産方法）
図１を参照して、本発明の植物の生産方法は、
少なくとも１つの植物からなる植物群の種子を播種し、発芽させて成長させることにより稚苗を得る、緑化工程と、
稚苗を育苗室に移植して成育させることにより、定植苗を得る、育苗工程と、
定植苗を栽培室に移植して成育させた後に収穫する本栽培工程と、を含む。
【００２４】
なお、植物の生産効率を高める観点からは、稚苗および定植苗は、苗選別ロボットなどにより選別された苗が移植されることが好ましい。
【００２５】
本発明の栽培（生産）対象となる植物としては、特に限定されないが、例えば、食用の植物、薬用の植物、観賞用の植物などが挙げられる。植物は、好ましくは食用の植物である。食用の植物としては、例えば、レタス、小松菜、ホウレンソウ、サンチュ、水菜、春菊、ハーブ類（ルッコラ、バジル、シソ、ツボクサ、ドクダミ等）などの葉物野菜、ナス、キュウリ、トマト、ピーマンなどの実物野菜、イチゴ、ミカン、マンゴー、ブドウ、ナシなどの果物；コメ、コムギ、オオムギ、ライムギ、エンバク、トウモロコシ、モロコシ、アワ、ヒエ、キビ等の穀物が挙げられる。苗の葉面積から重量（生重量）を推定し易い点からは、葉物野菜が好ましい。観賞用の植物としては、例えば、バラ、カーネーション、洋ラン、ガーベラ、トルコキキョウ等の花卉類；ポトス、セローム、アジアンタム等の観葉植物が挙げられる。なお、植物は、遺伝子組換え植物であってもよい。
【００２６】
本発明者らは、育苗工程における苗の生育のバラツキが、植物工場の生産を不安定化する要因の１つであることを見出した。育苗工程では、種の個体差や緑化工程および育苗工程における僅かな栽培条件の違いによって、苗の生育のバラツキが生じ易いと考えられる。また、育苗工程では、幼葉から本葉の形成など、生育における質的な変化が伴い、生育のバラツキが大きくなり易い傾向がある。
【００２７】
このため、本発明の植物の生産方法においては、まず、植物群を構成する個々の植物について、育苗工程の複数の時点における重量（生重量）のデータ（例えば、重量と重量分布のデータ）を収集する。そして、収集された重量のデータから、本栽培工程での植物（植物群）の成育量（重量）を予測する。
【００２８】
なお、植物群とは、少なくとも１つの植物（好ましくは複数の植物）からなる群であれば特に限定されないが、例えば、植物工場内の所定の領域に同日に播種された少なくとも１つの植物からなる群などが挙げられる。
【００２９】
具体的には、総栽培期間と、植物群の重量と、の関係式を作成し、該関係式からある総栽培期間の経過後に収穫された植物の重量の予測値を算出することができる。ここで、総栽培期間とは、収集された重量のデータから、緑化工程を実施する期間である緑化期間と、育苗工程を実施する期間である育苗期間と、本栽培工程を実施する期間である本栽培期間と、の合計期間である。
【００３０】
関係式は、植物の重量の予測値を算出することが可能であれば特に限定されないが、好ましくは、下記式（１）を重量のデータにフィッティングさせて、下記式（１）中の成長に関する係数Ｋ、Ｗ１を決定することにより得られる式である。
【００３１】
【数４】
2019041582.tif 000006
【００３２】
上記式（１）における各記号は、以下のとおりである。
ｔ：  総栽培期間（日）〔播種から（収穫まで）の合計期間〕
ｔ１：  緑化期間（日）〔播種から育苗室への移植までの期間〕
Ｗ（ｔ）：  植物重量（ｇ）〔総栽培期間ｔ（日）経過時における植物の重量〕
Ｗｍａｘ：  最大重量（ｇ）〔植物種ごとに通常想定される重量の最大値〕
Ｗ１：  稚苗重量（ｇ）〔育苗室への移植時（ｔ１経過時）における植物（稚苗）の重量〕
Ｋ：  成長率
なお、植物群が複数の植物からなる場合、Ｗ（ｔ）は、ある植物群についての平均値である。
【００３３】
重量のデータは、例えば、個々の植物の葉面積が植物の重量（地上部の生重量）に対して相関性を有している場合は、植物の葉面積の測定値に基づく推定により得ることができる。
【００３４】
ただし、葉面積は二次元画像の解析によって求められるため、植物が成長して葉の数が増えると、葉の重なりが生じるため、二次元画像から正確な葉面積を求めることができなくなる。このため、植物の成長を、例えば、ロジスティクス曲線（上記式（１））で近似して、成長曲線（植物の重量の予測値を示す曲線）を求める必要がある。
【００３５】
植物の成長を求めるためには、式（１）の成長パラメータのＷ１とＫの両方が必要である。Ｗ１とＫの両方を決めるためには、図８（ａ）に示されるような時系列データが必要である。
【００３６】
図８（ａ）の葉面積Ｌｉ（ｍｍ２）の時系列データに式（１）をフィッティングさせて、係数ＫとＷ１を求め、図８（ｂ）に点線で示されるような近似曲線を得ることができる。なお、ここでは、式（１）におけるＷｍａｘは、これまでの実験結果に基づいて定数（１８０ｇ）とした。
【００３７】
なお、近似曲線は、非線形ＧＲＧ（非線形計画）を用いて求めることができる。非線形ＧＲＧについては、例えば、ＵＲＬ［http://www.orsj.or.jp/archive2/or57-04/or57_4_175.pdf#search=%27GRG%E9%9D%9E%E7%B7%9A%E5%BD%A2%27］のp.181などを参照できる。
【００３８】
なお、苗の重量（葉面積）の測定は、植物群の全ての株について行う必要はなく、植物群の一部の株について測定を行い、それらの測定値から植物群全体の測定値を推定してもよい。なお、植物工場内の生育バラツキ（分散）は非常に大きいため、例えば、少なくとも生産全体の１０％程度の苗について重量（葉面積）の測定を行うことが好ましい。例えば、日産５０００株の植物工場の場合、少なくとも５００株について測定を実施することが好ましい。
【００３９】
そして、本発明の植物の生産方法においては、さらに、該関係式に基づいて、植物群について、収穫される植物のうち所定範囲内の重量を有する植物の比率が最適化されるように、総栽培期間が調整される。
【００４０】
「収穫される植物のうち所定範囲内の重量を有する植物の比率が最適化される」とは、例えば、植物群のうち所定の閾値重量ｗ＊以上の植物の重量当たりの価格と、植物群のうち閾値重量未満の植物の重量当たりの価格と、の比率に基づいて、植物群の生産および販売によって得られる利益率が最大化するようにすることである。すなわち、できるかぎり短い総栽培期間で植物の生産に必要なコストを抑えつつ、収穫される植物群のうち所定の閾値重量ｗ＊以上の重量を有する植物の比率ができるだけ多くなるようにすることで、当該利益率を最大化することができる。
【００４１】
例えば、あらかじめ、植物の生育量（重量）に応じた異なる販売形態（袋売りと重量売り）を考慮した上で、販売形態ごとに収益性（利益率）を求めるダイアグラムを構築し、全ての販売形態の合計の収益率を求める計算式（例えば、上記式（２）および式（３））を作成しておく。なお、上記式（２）および式（３）においては、袋売りされる植物の重量の最小値である閾値重量ｗ＊（出荷閾値）が計算式の１つのパラメータとなっている。そして、このような計算式を用いて、全ての販売形態の合計の収益率が最大となる総栽培期間ｔ（最適栽培期間）等を求める。
【００４２】
このようにして求められた最適栽培期間に従って、植物群の総栽培期間を調整することで、植物群の利益率を最大化することができる。なお、例えば、後述する実施形態１のように、本栽培期間を調整することにより、総栽培期間を調整してもよく、実施形態２のように、育苗期間を調整することにより、総栽培期間を調整してもよい。
【００４３】
例えば、植物群が植物工場内の所定の領域に同日に播種された植物であり、植物工場内で複数の植物群が栽培されている場合、それらの植物群ごとに以上のプロセスによる総栽培期間の調整を継続的に実施することで、継続的に植物工場等の収益性を最大化することができる。また、植物工場等における各植物群ごとの収益性（利益率）の変動幅を最小限に抑えることができ、収益性を安定化することができる。
【００４４】
上記のとおり、本発明の植物の生産方法によれば、植物工場などにおいて植物の生育（成長速度）にバラツキがある場合でも、収穫される植物の重量を所望の範囲内（規格内）に制御することができ、生産安定性が向上する。したがって、販売予測が可能となり、販売計画を立てることができる（例えば、播種から３週間経過時に１箇月後の計画が立てられる等）という利点がある。このため、本発明の植物の生産方法は、契約販売等の場合に特に有利である。
【００４５】
なお、総栽培期間の調整において、利益率のダイアグラムに将来的な市場（需要状況）の情報等をフィードフォワードしておいてもよい。例えば、１週間後に袋物の価格が低下して、重量売りの価格が上がる等の情報を、あらかじめλ（袋売り商品の重量当たりの価格に対する重量売り商品の重量当たりの価格の割合）などに反映させておいてもよい。また、購買率Ｐなどの変動のあるパラメータについて、フィードバック制御を行うことで、より高精度に収益性を向上させることが可能となる。
【００４６】
（システム）
図２は、本発明の植物の生産方法に用いられるシステムの一例を示す模式図である。
【００４７】
図２および図３（ａ）を参照して、葉面積の測定には、例えば、超小型ＰＣとカメラから構成される３０台の撮像装置（ＰＩ：Plant  Imager）からなる撮像装置群（ＭＰＩ：Multi  Plant  Imager）が用いられる。ＭＰＩは、全体で５００株の植物の撮像が可能である。なお、本栽培工程において葉面積をモニターすると大面積で複雑な撮像装置が必要となってしまうため、育苗工程で葉面積をモニターする場合、撮像装置等の装置構成を簡略化できる利点もある。
【００４８】
ＭＰＩで得られた画像データは、例えば、無線ＬＡＮアクセスポイントおよびＨＵＢを介して、解析用のメインＰＣおよびデータ保管用のファイルサーバ（または、ＮＡＳ：ネットワークハードディスク）に送信される。なお、画像データの送信は、特に無線送信に限られず、有線送信で実施されてもよい。
【００４９】
そして、メインＰＣにて、ＭＰＩによって得られた個々の植物の画像に対して、個体認識および画像処理（葉面積計測）を行うことで、個々の植物の葉面積を測定することができる（図３（ｂ）参照）。次に、メインＰＣは、葉面積の測定値から、個々の植物の重量（地上部の生重量）の推定値を算出する。
【００５０】
図７を参照して、図７（ｂ）に示されるように画像の解析処理により測定された個々の植物の葉面積Ｌｉの測定値は、図７（ａ）に示されるように、植物の重量（地上部の生重量）の測定値に対して相関性（相関係数Ｒ＝０．８７）を有していることが分かる。したがって、図７（ａ）に示されるような検量線を用いて、葉面積の測定値から、個々の植物の重量の推定値を算出することができる。
【００５１】
さらに、メインＰＣは、上記のようにして得られた個々の植物の重量データ等に基づいて、植物群の成育速度（ある総栽培期間経過時の植物群の重量予測値）を求めるアルゴリズムを構築する。
【００５２】
メインＰＣには、植物群の生産および販売によって得られる利益率を求めるアルゴリズムと、その利益率が最大化するような最適収穫期間等を求める収益最大化アルゴリズム（図３（ｃ））が実装されている。その収益最大化アルゴリズムによって、最適収穫日が算出される。なお、メインＰＣには、例えば、そのような最適収穫日を工場長などにアラートする機能が実装されていてもよい。
【００５３】
収益最大化アルゴリズムは多数のパラメータを持ち、それは植物工場毎に異なる固有の数値となっている。例えば、作物品種、電力コスト、標準栽培期間などである。収益最大化アルゴリズムの構築には、これらの緻密な数値が重要であるが、電力コストなどは天候にも依存しており（例えば、猛暑日の冷房コスト）、一般には複雑である。植物工場運営会社とも連携した緻密な現地調査により、緻密なパラメータ数値の獲得を行うことが望ましい。
【００５４】
上記のシステムは、植物を栽培するための光を照射する光源と、光源の点灯を制御する光源制御手段を備えていてもよい。また、システムは、植物を収容する空間内の湿度や温度を管理する空調設備、植物栽培用の培地を供給する培養設備などを備えていてもよい。
【００５５】
光源はＬＥＤや蛍光灯などから構成され、植物の生育に必要な波長の光を照射する。光源は固定された波長の光を照射する光源や複数の波長の光を照射する光源であり得る。照射する光の波長を制御可能な光源が好ましく用いられる。植物の成長は照射される光の波長によって影響されることも知られており、光の波長は適宜選択される。
【００５６】
なお、植物工場における「夜」に相当する時間は、植物工場内は暗黒となっている。このため、「夜」に相当する時間に鮮明な植物の画像を得るためには、栽培用ＬＥＤ照明（光源）を点灯させる必要があるが、ＬＥＤ証明制御装置によるリレー制御などによって、栽培用ＬＥＤ証明の点灯が植物の生育に干渉しないように制御することが好ましい（図２参照）。
【００５７】
なお、植物の栽培培地としては、植物を栽培（本栽培）することのできる培地であれば特に限定されず、例えば、水耕栽培等に用いられる種々公知の培地を使用できる。栽培培地は、例えば、移植された植物を識別する識別情報を記録する記録媒体が設けられた栽培容器に収容される。ＭＰＩまたはメインＰＣは、このような記録媒体の識別情報を読み取ることで、植物の個体認識を行うことができる。
【００５８】
複数の栽培容器は、例えば、棚板が所定間隔を隔てて複数段設けられた栽培棚に載置される。上記のシステムは、さらに、複数の棚板間で栽培容器を搬送する搬送装置などを備えていてもよい。
【００５９】
以下、本発明の植物の生産方法（特に総栽培期間の調整）の具体例について、実施形態を挙げて説明する。
【００６０】
＜実施形態１＞
本実施形態においては、本栽培期間を調整することにより、総栽培期間が調整される。
【００６１】
本実施形態においては、本栽培期間の調整に用いられる利益率Ｙの式として、例えば、下記式（２）が用いられる。
【００６２】
【数５】
2019041582.tif 000007
【００６３】
上記式（２）における各記号は、以下のとおりである。
Ｙ（ｔ，ｗ＊，λ）：  植物の１株当たりの利益率（／株）
ｗ＊：  閾値重量（ｇ）〔袋詰めされる重さの最小値〕
ｒ＊：  重量がｗ＊以上の株（袋売りの株）の割合（／株）
λ：  袋売り商品の重量当たりの価格に対する重量売り商品（カット野菜等）の重量当たりの価格の割合
<Ｗ>：  重量がｗ＊未満の株（重量売りの株）の平均重量（ｇ）
Ｐ：  袋売り商品を消費者が購入する割合である購買率
α：  固定費（播種数一定の場合、光熱費以外はすべて含まれる）（／株）
β：  １日当りの１株ごとの緑化工程および育苗工程（播種から定植まで）における光熱費（／（日・株））
γ：  １日当りの１株ごとの本栽培工程（定植から収穫まで）における光熱費（／（日・株））
ｔ１は、前記緑化期間（日）〔播種から移植までの期間〕
ｔ２は、前記育苗期間（日）〔移植から定植までの期間〕
【００６４】
ここで、<Ｗ>は、例えば、下記式（４）を用いて算出することができる。下記式（４）において、σ（ｔ）は、Ｗ（ｔ）の標準偏差である。
【００６５】
【数６】
2019041582.tif 000008
【００６６】
また、ｒ＊は、例えば、下記式（５）を用いて算出することができる。
【００６７】
【数７】
2019041582.tif 000009
【００６８】
なお、Ｐは、例えば、下記式（６）を用いて算出することができる。式（６）において、ｗ’は、消費者が実際の商品（植物）を見て１００％購入する最小の重さの基準である。
【００６９】
【数８】
2019041582.tif 000010
【００７０】
上記式（６）の植物の閾値重量ｗ＊と購買率Ｐとの関係は、図４のグラフで示される。なお、図４では、式（６）のパラメータが、ａ＝０．０２、ｂ＝−０．５，ｗ’＝７５（ｇ）である。図４では、閾値重量ｗ＊をｗ’（図４では、ｗ’＝７５ｇ）以上に設定すれば、袋詰め商品の購買率Ｐが１００％になることが示されている。
【００７１】
本実施形態においては、図９に示されるような上記式（２）の利益率Ｙのグラフに基づいて、植物群の最適な総栽培期間（最適栽培期間）を決定することができる。
【００７２】
なお、図９は、上記式（１）で、Ｗｍａｘ＝１７８．２５（ｇ）、Ｗ１＝０．１５（ｇ）、Ｋ＝０．２０と設定し、上記式（２）で、α＝０．５５２、β＝０．０００８２、γ＝０．００５５８、λ＝０．３０、ｔ1＝６（日）、ｔ２＝１４（日）と設定されたときのグラフである。なお、Ｗ（ｔ）のｃ．ｖ．（変動係数）は０．２３４であった。そして、図９において、閾値重量ｗ＊が７５ｇのときの総栽培期間ｔの最適値（最適栽培期間）は、４７日であり、そのときの利益率Ｙは０．２５９（２５．９％）であった。
【００７３】
この最適栽培期間に基づいて、本栽培期間を調整する（本栽培工程の段階で総栽培期間を調整する）ことにより、植物群間で同等の成長を実現することが可能となる（図５（ａ）および（ｂ）参照）。すなわち、同じ育苗期間で本栽培を開始した（同日に栽培室への移植を行った）場合、（本栽培期間（定植から収穫までの期間）を調整することで、所望の重量を有する植物の収穫が可能となる（図５（ｂ）参照）。
【００７４】
図５（ａ）において、Ｗ０は、播種時重量（ｇ）〔ｔ＝０の時点における植物（種子）の重量〕である。また、Ｗ２は、定植苗重量（ｇ）〔栽培室への移植時（ｔ１＋ｔ２経過時：図５（ａ）では２０日経過時）における植物（定植苗）の重量〕である。ここで、ｔ２は、育苗期間（日）〔育苗室への移植から栽培室への移植（定植）までの期間〕である。
【００７５】
本実施形態においては、苗成育指標のバラツキに応じて、総栽培期間（本栽培期間）を調整することで、従来の方法（図１０参照）に比べて、植物群の利益率の平均値＜Ｙ＞が安定化され、生産安定化が実現される（図１１参照）。なお、図１０および図１１における横軸は、植物工場の稼働日数を示している。図１０および図１１は、例えば、１００日間定期的に一定量の植物群が播種された場合において、定期的に播種された各植物群がＭＰＩ計測終了時から何日後に収穫されたかを示している。このように、定期的に播種し栽培される植物群ごとに、総栽培期間の調整を継続的に実施することで、継続的に各植物群場等の収益性を最大化することができる。また、各植物群ごとの収益性（利益率）の変動幅を最小限に抑えることができ、収益性を安定化することができる。
【００７６】
＜実施形態２＞
実施形態２は、育苗期間を調整することにより、総栽培期間が調整される点で、実施形態１（本栽培期間を調整することにより、総栽培期間が調整される）とは異なる。それ以外の点は、実施形態１と基本的に同様であるため、重複する説明は繰り返さない。
【００７７】
本実施形態においては、栽培期間の調整に用いられる利益率Ｙの式として、例えば、下記式（３）が用いられる。
【００７８】
【数９】
2019041582.tif 000011
【００７９】
本実施形態においては、図１２に示されるような上記式（３）の利益率Ｙのグラフに基づいて、植物群の最適な総栽培期間（最適栽培期間）を決定することができる。
【００８０】
なお、図１２は、上記式（１）で、Ｗｍａｘ＝１７８．２５（ｇ）、Ｗ１＝０．１５（ｇ）、Ｋ＝０．２０と設定し、上記式（２）で、α＝０．５５２、β＝０．０００８２、γ＝０．００５５８、λ＝０．３０、ｔ３＝２０（日）と設定されたときのグラフである。なお、Ｗ（ｔ）のｃ．ｖ．（変動係数）は０．２３４であった。そして、図１２において、閾値重量ｗ＊が７５ｇのときの総栽培期間ｔの最適値（最適栽培期間）は、５０日であり、そのときの利益率は０．３０３（３０．３％）であった。
【００８１】
この最適栽培期間に基づいて、育苗期間を調整する（育苗工程の段階で総栽培期間を調整する）ことにより（図６（ａ）参照）、（本栽培期間の調整を必要とせずに）同じ本栽培期間で収穫しても、生育の異なる植物群間で同等の成長を実現することが可能となる（図６（ｂ）参照）。なお、植物群Ａ〜Ｃの播種から苗の定植までの日数は、植物群Ａで２０−４日、植物群Ｂで２０日、植物群Ｃで２０＋１１日であった。
【００８２】
なお、図６（ｂ）において、Ｗ３は、収穫時重量（ｇ）〔収穫時（ｔ＝４０日＋調整期間α日経過時］における植物の重量〕である。
【００８３】
本実施形態においては、苗成育指標のバラツキに応じて、総栽培期間（育苗期間）を調整することで、従来の方法（図１３参照）に比べて、植物群の利益率の平均値＜Ｙ＞が安定化され、生産安定化が実現される（図１４参照）。なお、図１３および図１４における横軸は、植物工場の稼働日数を示している。図１３および図１４は、例えば、１００日間定期的に一定量の植物群が播種された場合において、定期的に播種された各植物群がＭＰＩ計測終了時から何日後に収穫されたかを示している。このように、定期的に播種し栽培される植物群ごとに、総栽培期間の調整を継続的に実施することで、継続的に各植物群場等の収益性を最大化することができる。また、各植物群ごとの収益性（利益率）の変動幅を最小限に抑えることができ、収益性を安定化することができる。
【００８４】
今回開示された実施形態および実施例はすべての点で例示であって制限的なものではないと考えられるべきである。本発明の範囲は上記した説明ではなくて特許請求の範囲によって示され、特許請求の範囲と均等の意味および範囲内でのすべての変更が含まれることが意図される。
【図１】
2019041582.tif 000015
【図２】
2019041582.tif 000016
【図３】
2019041582.tif 000017
【図４】
2019041582.tif 000018
【図５】
2019041582.tif 000019
【図６】
2019041582.tif 000020
【図７】
2019041582.tif 000021
【図８】
2019041582.tif 000022
【図９】
2019041582.tif 000023
【図１０】
2019041582.tif 000024
【図１１】
2019041582.tif 000025
【図１２】
2019041582.tif 000026
【図１３】
2019041582.tif 000027
【図１４】
2019041582.tif 000028