(57)【要約】  本技術は、生体の内部の状態を判断することができるようにする情報処理装置、情報処理方法、プログラム、及び、センシング装置に関する。情報処理装置は、測定対象の生体から得られる過渡応答の測定データに対し、生体が有する機能を示す過渡モデルに応じた計算アルゴリズムを適用して、過渡モデルに関するパラメータのうち、未知のパラメータを算出することで、生体の内部の状態を判断することができるようにする。本技術は、例えば、植物の光合成の量子収率を算出する情報処理装置に適用することができる。
2018056102.tif 000002
【特許請求の範囲】
【請求項１】
測定対象の生体から得られる過渡応答の測定データに対し、前記生体が有する機能を示す過渡モデルに応じた計算アルゴリズムを適用して、前記過渡モデルに関するパラメータのうち、未知のパラメータを算出する算出部を備える
情報処理装置。
【請求項２】
前記生体は、植物であり、
前記過渡モデルは、前記植物の光合成での光化学反応の電子伝達系を表現したモデルであり、
前記計算アルゴリズムは、前記過渡モデルを用い、前記植物のクロロフィル蛍光の過渡応答から光合成の量子収率を算出するためのアルゴリズムであり、
前記測定データは、前記植物のクロロフィル蛍光の過渡応答の測定データであり、
前記算出部は、前記測定データを入力とし、前記計算アルゴリズムを適用することで、未知のパラメータとして、前記植物の光合成の量子収率を算出する
請求項１に記載の情報処理装置。
【請求項３】
前記植物の光合成の量子収率は、光合成の電子伝達系の最大量子収率（Fv/Fm）及び光合成の電子伝達系の量子収率（ΦPSII）のうち少なくとも一方である
請求項２に記載の情報処理装置。
【請求項４】
前記過渡モデルは、１又は複数用意され、
前記計算アルゴリズムは、前記過渡モデルに応じて１又は複数用意される
請求項１に記載の情報処理装置。
【請求項５】
前記過渡モデルと前記計算アルゴリズムを記憶する記憶部をさらに備える
請求項１に記載の情報処理装置。
【請求項６】
前記植物に対する光を発光する発光部を、前記測定データが得られるように制御する発光制御部をさらに備える
請求項２に記載の情報処理装置。
【請求項７】
前記発光部と、
前記植物のクロロフィル蛍光をセンシングするセンサと
をさらに備える請求項６に記載の情報処理装置。
【請求項８】
レーザ光を発光する前記発光部と、
前記レーザ光を走査するための走査部と、
前記植物のクロロフィル蛍光をセンシングするセンサと
をさらに備える請求項６に記載の情報処理装置。
【請求項９】
前記発光部は、前記植物の育成用の光を発光する外部の発光装置であり、
前記発光装置と接続されるインターフェースと、
前記植物のクロロフィル蛍光をセンシングするセンサと
をさらに備え、
前記発光制御部は、前記インターフェースを介して、前記発光装置を、前記測定データが得られるように制御する
請求項６に記載の情報処理装置。
【請求項１０】
前記植物のクロロフィル蛍光をセンシングする第１のセンサと、
前記第１のセンサによるセンシング時に、環境光をセンシングする第２のセンサと
をさらに備え、
前記算出部は、前記環境光の変化時に得られた前記測定データに基づいて、前記植物の光合成の量子収率を算出する
請求項２に記載の情報処理装置。
【請求項１１】
前記環境光は、太陽光である
請求項１０に記載の情報処理装置。
【請求項１２】
前記測定データを提供する外部の処理装置と接続されるインターフェースをさらに備え、
前記算出部は、前記インターフェースを介して、前記処理装置から提供される前記測定データに基づいて、前記植物の光合成の量子収率を算出する
請求項２に記載の情報処理装置。
【請求項１３】
前記処理装置は、前記測定データを計測する計測装置、又は前記測定データを記録する記録装置である
請求項１２に記載の情報処理装置。
【請求項１４】
情報処理装置の情報処理方法において、
前記情報処理装置が、
測定対象の生体から得られる過渡応答の測定データに対し、前記生体が有する機能を示す過渡モデルに応じた計算アルゴリズムを適用して、前記過渡モデルに関するパラメータのうち、未知のパラメータを算出する
ステップを含む情報処理方法。
【請求項１５】
コンピュータを、
測定対象の生体から得られる過渡応答の測定データに対し、前記生体が有する機能を示す過渡モデルに応じた計算アルゴリズムを適用して、前記過渡モデルに関するパラメータのうち、未知のパラメータを算出する算出部を備える情報処理装置として
機能させるためのプログラム。
【請求項１６】
測定対象の生体をセンシングするセンサと、
前記センサによるセンシングによって前記生体から得られる過渡応答の測定データに対し、前記生体が有する機能を示す過渡モデルに応じた計算アルゴリズムを適用して、前記過渡モデルに関するパラメータのうち、未知のパラメータを算出する算出部と
を備えるセンシング装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本技術は、情報処理装置、情報処理方法、プログラム、及び、センシング装置に関し、特に、生体の内部の状態を判断することができるようにした情報処理装置、情報処理方法、プログラム、及び、センシング装置に関する。
【背景技術】
【０００２】
植物の成長は、日照や温度などの環境の状態に影響を受けるが、同じ環境においても、植物の種類や環境への順化状態によって異なる挙動を示すため、植物を良好に育成するには、環境の状態を把握するだけでなく、植物の生育を診断し、その状態に応じて育成環境を整える必要がある。【０００３】
この種の植物の生育の診断方法としては、例えば、特許文献１に開示された技術が知られている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開平11-332375号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
ところで、一般的に植物の生育などの生体の内部の状態を、迅速かつ正確に判断することは難しく、そのような判断をできるようにすることが望まれていた。【０００６】
本技術はこのような状況に鑑みてなされたものであり、生体の内部の状態を判断することができるようにするものである。
【課題を解決するための手段】
【０００７】
本技術の一側面の情報処理装置は、測定対象の生体から得られる過渡応答の測定データに対し、前記生体が有する機能を示す過渡モデルに応じた計算アルゴリズムを適用して、前記過渡モデルに関するパラメータのうち、未知のパラメータを算出する算出部を備える情報処理装置である。【０００８】
本技術の一側面の情報処理方法又はプログラムは、上述した本技術の一側面の情報処理装置に対応する情報処理方法又はプログラムである。【０００９】
本技術の一側面の情報処理装置、情報処理方法、及び、プログラムにおいては、測定対象の生体から得られる過渡応答の測定データに対し、前記生体が有する機能を示す過渡モデルに応じた計算アルゴリズムが適用され、前記過渡モデルに関するパラメータのうち、未知のパラメータが算出される。【００１０】
本技術の一側面のセンシング装置は、測定対象の生体をセンシングするセンサと、前記センサによるセンシングによって前記生体から得られる過渡応答の測定データに対し、前記生体が有する機能を示す過渡モデルに応じた計算アルゴリズムを適用して、前記過渡モデルに関するパラメータのうち、未知のパラメータを算出する算出部とを備えるセンシング装置である。【００１１】
本技術の一側面のセンシング装置においては、測定対象の生体がセンシングされ、センシングによって前記生体から得られる過渡応答の測定データに対し、前記生体が有する機能を示す過渡モデルに応じた計算アルゴリズムが適用され、前記過渡モデルに関するパラメータのうち、未知のパラメータが算出される。【００１２】
なお、情報処理装置及びセンシング装置は、独立した装置であってもよいし、１つの装置を構成している内部ブロックであってもよい。
【発明の効果】
【００１３】
本技術の一側面によれば、生体の内部の状態を判断することができる。【００１４】
なお、ここに記載された効果は必ずしも限定されるものではなく、本開示中に記載されたいずれかの効果であってもよい。
【図面の簡単な説明】
【００１５】
【図１】植物の光合成の仕組みを示す図である。
【図２】植物のクロロフィル蛍光を説明するための図である。
【図３】ショートレスポンスとロングレスポンスを説明するための図である。
【図４】光合成の量子収率を説明するための図である。
【図５】パルス変調蛍光測定を用いた光合成計測装置の構成を示す図である。
【図６】植物に照射される光の発光強度とクロロフィル蛍光の強度との関係を示す図である。
【図７】光合成の量子収率を算出するための計算アルゴリズムを電気的なモデルで表した図である。
【図８】パルス変調蛍光測定に関する補足を説明するための図である。
【図９】本技術の過渡モデルと光合成量子収率計算アルゴリズムを用いたパラメータ推定を概念的に表した図である。
【図１０】第１の実施の形態の光合成計測装置の構成例を示す図である。
【図１１】植物に照射される光の発光強度とクロロフィル蛍光の強度との関係を示す図である。
【図１２】光合成計測処理を説明するフローチャートである。
【図１３】過渡モデルの構成を表した図である。
【図１４】過渡モデルのシミュレーションの結果の例を示す図である。
【図１５】時定数の計算アルゴリズムを説明する図である。
【図１６】Fv/Fmを算出する場合の計算アルゴリズムの適用までの流れを説明する図である。
【図１７】ΦPSIIを算出する場合の計算アルゴリズムの適用までの流れを説明する図である。
【図１８】現状の方式と本技術の方式の測定時に用いられる光の違いを説明する図である。
【図１９】第２の実施の形態の光合成計測装置の構成例を示す図である。
【図２０】第３の実施の形態の光合成計測装置の構成例を示す図である。
【図２１】第４の実施の形態の光合成計測装置の構成例を示す図である。
【図２２】第５の実施の形態の情報処理装置の構成例を示す図である。
【図２３】植物の測定を行う測定装置の具体例を示す図である。
【図２４】コンピュータの構成例を示す図である。
【発明を実施するための形態】
【００１６】
以下、図面を参照しながら本技術の実施の形態について説明する。なお、説明は以下の順序で行うものとする。
【００１７】
１．本技術の概要
２．第１の実施の形態：基本構成
３．第２の実施の形態：レーザ光を用いた構成
４．第３の実施の形態：外部の発光装置を用いた構成
５．第４の実施の形態：環境光の強度変化を用いた構成
６．第５の実施の形態：外部からの測定データを用いた構成
７．変形例
８．コンピュータの構成
【００１８】
＜１．本技術の概要＞
【００１９】
植物の光合成活動及びそれに伴う植物体の成長は、日照や温度、飽差、二酸化炭素の濃度（CO2濃度）、土壌水分、土壌中の肥料成分といった環境の状態に影響を受けるとともに、同じ環境においても、植物の種類や環境への順化状態によって異なる挙動を示す。従って、植物の良好な育成のためには、環境の状態を把握するだけでなく、植物の種類や育成ステージごとに、時間的変化を伴う植物の内部状態をモニタ（監視）し、その状態に応じて育成環境を制御することが必要となる。
【００２０】
しかしながら、一般的に、植物の内部状態を、迅速にかつ正確にモニタすることは、極めて難しい。例えば、日照は十分であっても、温度が低すぎて光合成が十分に行えない状態にあることは、過去の経験から判断するか、あるいはその状態が一定のまま数日が経過して成長が十分でなかった結果が観測できない限りは、その判別が難しい。また、環境条件が一定ではない野外環境等においては、植物の育成に影響を与えている環境条件の特定や定量化、あるいは状態に応じた対処を行うことは、困難である。
【００２１】
本技術は、このような現状を鑑みて、植物の内部の状態を判断（診断）することができるようにするものであるが、ここでは、まず、本技術を説明する上で前提となる光合成の仕組みや現状の技術などについて説明してから、本技術の具体的な内容について説明する。
【００２２】
（光合成の仕組み）
図１は、植物の光合成の仕組みを説明するための図である。
【００２３】
図１に示すように、光合成は、大別すると、光化学反応（Photochemical System）と炭素還元反応（Calvin-Benson Cycle）から構成される。光化学反応は、光のエネルギーで水（H2O）を分解してエネルギーを得る過程である。炭素還元反応は、前段の光化学反応で得られたエネルギーを用いて、二酸化炭素（CO2）を分解し、糖（CH2O）を生成する過程である。
【００２４】
より具体的には、光化学反応は、受け取る光の波長帯域の違いにより、光化学系I（PSI）と、光化学系II（PSII）に分けられる。図１に示すように、光化学反応は、主に、光を集める集光アンテナ（LHCI，LHCII）と、光のエネルギーで電子を励起させる反応中心（RCI，RCII：図１ではP680，P700と記述）と、電子伝達系に沿って、電子の伝搬や一時貯蔵を行う中間プール（QA，QB，PQ，A0等）により構成される。
【００２５】
そして、この電子伝達系での電子のエネルギー勾配を示すＺ機構（Z-scheme）の酸化還元反応により、ニコチンアミドアデニンジヌクレオチドリン酸（NADPH：Nicotinamide Adenine Dinucleotide Phosphate）が合成される。ニコチンアミドアデニンジヌクレオチドリン酸は、生体内のどこにでも存在し、酸化型（NADP+）と還元型（NADPH）の２つの状態があり、電子や水素を運ぶ役割を有している。ただし、酸化型（NADP+）では電子を受け取りやすく、還元型（NADPH）では電子を放出しやすい。
【００２６】
また、光化学反応では、酸素発生複合体（OEC：Oxygen Evolving Complex）によって、水（H2O）が分解され、酸素（CO2）が発生する。このときに得られるプロトン（H+）のエネルギーによって、ATP合成酵素（ATP synthase）が、アデノシン二リン酸（ADP：Adenosine Diphosphate）から、アデノシン三リン酸（ATP：Adenosine Triphosphate）を合成している。
【００２７】
なお、光化学反応には、処理できない余剰のエネルギーを熱（Heat）として発散するために、余剰のエネルギーの大きさに従い、化学組成を変えながら適切な量のエネルギーを熱として放出することができる、キサントフィルサイクル（Xanthophyll Cycle）が存在している。
【００２８】
一方で、炭素還元反応では、光化学反応で作られたニコチンアミドアデニンジヌクレオチドリン酸（NADPH）とアデノシン三リン酸（ATP）を使って、二酸化炭素（CO2）が分解され、糖（CH2O）が作られる。
【００２９】
（クロロフィル蛍光）
クロロフィル蛍光（Chlorophyll Fluorescence）は、植物の光合成に伴い、植物から発せられる蛍光であり、光により電子が励起した反応中心（RC：Reaction Center）から一定時間内にエネルギーが抜き取られないと、高等植物では、680〜750nm前後の波長の蛍光としてエネルギーが放出される現象である。なお、放射されるエネルギーは、入射光のエネルギーに対し、0.5〜3％とされ、植物の光合成の状態に応じて変動するものである。図２には、植物のクロロフィル蛍光を模式的に表している。
【００３０】
ここで、電子伝達系の中間プールが酸化状態で、電子を受け入れることが可能な状態である場合を、反応中心（RC）がオープン（Open）であるといい、光により励起した電子は、電子伝達系を通じて流れ出し、励起状態が解消される。これを、光化学消光（PhotoChemicalQuenching）と呼ぶ。一般的に、光合成が効率よく進行している場合、光化学消光の度合いが大きいため、クロロフィル蛍光の強度は小さくなる。
【００３１】
一方、電子伝達系の中間プールが還元状態で、電子を受け入れることができない場合を、反応中心（RC）がクローズ（Close）であるという。すべての反応中心（RC）がクローズしたとき、光化学消光がなくなり（電子伝達系への電子伝達が止まり）、クロロフィル蛍光の強度は最大となる。
【００３２】
光より励起した電子が、クロロフィル蛍光（Fluorescence）又は熱（Heat）として発散されるとき、これを、非化学消光（Non- PhotoChemicalQuenching）と呼ぶ。植物に光化学系で処理しきれない強い光を与えた場合、一時的にクロロフィル蛍光の強度が上昇するが、時間が経過するにつれて、炭素還元反応の活性化や、キサントフィルサイクルの化学組成が変化し、熱を効率よく放出できるようになると、クロロフィル蛍光の強度は減少する。
【００３３】
（電子伝達系の反応速度）
光化学反応において、光の入力に対する電子伝達系の反応速度は、ピコ秒〜数十ミリ秒と高速な反応を示し、これを、ショートレスポンス（Short Response）と呼ぶ。図３には、光化学反応における電子伝達系の反応速度が、ピコ秒〜数十ミリ秒となることが記述されている。
【００３４】
一方で、炭素還元反応と、キサントフィルサイクルは、活性化に数秒以上かかり、光の強さなどに応じて数分かけて定常状態に遷移する。このことを、ロングレスポンス（Long Response）と呼ぶ。図３には、炭素還元反応とキサントフィルサイクルの活性化に、数秒〜数分かかることが記述されている。なお、図３に示したタイミングパラメータの値は一例である。
【００３５】
（光合成の量子収率）
光合成の量子収率としては、例えば、Fv/Fmと、ΦPSIIがあり、これらの値により、現在の環境条件において、植物がどの程度の効率で光合成を行っているかを、リアルタイムに把握すること可能となる。
【００３６】
Fv/Fmは、光合成の光化学系II（PSII）の電子伝達系の最大量子収率であり、炭素還元反応やキサントフィルサイクルに律速することのない暗順応状態における光化学系IIの量子収率（入力光子数に対する出力電子数）を表している。
【００３７】
ΦPSIIは、光合成の光化学系II（PSII）の電子伝達系の量子収率であり、一定の強度の光が当たって光合成が行われているときの量子収率を表している。
【００３８】
例えば、図４に示すように、入力光の強度I（PPFD[μmol/m2/s]）、植物（群落）の光の吸収率（一般的に0.84がよく用いられる）、入力光の光化学系I（PSI）と光化学系II（PSII）との分配率（一般的に定数0.5がよく用いられる）、光合成の光化学系IIの電子伝達系の量子収率（ΦPSII）を用いると（仮に0.7が計測されたものとする）、電子伝達速度（ETR：Electron Transport Rate）を算出することができる。
【００３９】
例えば、PPFD値が1000[μmol/m2/s]となる場合の電子伝達速度（ETR）は、ETR ＝ 1000×0.84×0.5×0.7 ＝ 294[μmol/m2/s]となる。つまり、ETR ＝ I×0.84×0.5×ΦPSII[μmol/m2/s]となる。電子伝達速度（ETR）は、現在の環境条件や植物の状態において、植物がどのような速度で光を電子に変えているかを、リアルタイムで定量的に表現するパラメータである。
【００４０】
ここで、PPFD値は、光合成光量子密度（PPFD：Photosynthetic Photon Flux Density）の値であって、葉緑素（クロロフィル）の吸収波長である400nmから700nmの波長での光量子が、単位時間、単位面積あたりに入射する個数で表したものである。
【００４１】
さらに、電子伝達系からニコチンアミドアデニンジヌクレオチドリン酸（NADPH）の効率を0.5、炭素還元反応の効率を0.5とした場合、GPP(Gross Primary Production)は、GPP ＝ ETR×0.5×0.5 ＝ 73.5[μmol/m2/s]となる。
【００４２】
なお、炭素還元反応の効率は、温度やCO2濃度により変動するが、炭素還元反応の効率が落ち、受け入れられるエネルギーの上限に制限が生じている場合、電子伝達速度（ETR）が低下することにつながるため、ETRとGPPとの相関が高いことが知られている。
【００４３】
（パルス変調蛍光測定）
ところで、環境や植物の状態に応じた光合成の状態をリアルタイムに計測するための手法として、パルス変調蛍光測定と称される手法がある。このパルス変調蛍光測定の手法を用いた光合成計測器は既に実用化されている。パルス変調蛍光測定は、太陽光の強度を超えるような非常に強い人工光（以下、飽和光という）を植物に照射することで、当該植物の光合成の状態を測定するものである。
【００４４】
図５は、パルス変調蛍光測定を用いた光合成計測装置９０の構成を示す図である。
【００４５】
光合成計測装置９０は、飽和光を用いて、植物の光合成の量子収率の計測を行う。光合成計測装置９０は、中央制御部９００、発光部９０１、発光制御部９０２、測定光生成部９０３、合成部９０４、光学フィルタ９０５、光センサ９０６、測定光分離電気的フィルタ９０７、サンプリング部９０８、光合成量子収率演算部９０９、及び結果表示部９１０から構成される。
【００４６】
発光制御部９０２は、中央制御部９００からの制御に従い、測定光生成部９０３により生成されるパルス変調された測定光を、合成部９０４により合成することで、発光部９０１から、所定の強度の光が発光されるようにする。
【００４７】
発光部９０１からは、LED(Light Emitting Diode)によって、３種類の強度の光が発光される。３種類の強度の光には、測定光からなる第１の光（以下、「測定光」と記述する）と、定常光と測定光からなる第２の光（以下、「定常光＋測定光」と記述する）と、飽和光と測定光からなる第３の光（以下、「飽和光＋測定光」と記述する）とが含まれる。
【００４８】
ここで、測定光は、反応中心（RC）が、オープン（Open）の状態のクロロフィル蛍光（F0等）を測定するための非常に弱い光である。飽和光は、反応中心（RC）を完全にクローズ（Close）させた状態のクロロフィル蛍光（Fm，Fm'）を測定するための非常に強い光である。この飽和光は、0.5〜1秒前後の継続時間を持った光となる。
【００４９】
定常光は、植物１に光合成を定常的に行わせた状態のクロロフィル蛍光（Fs）を測定するための光である。この定常光は、ロングレスポンスを観測するため、植物１の光合成が安定するまで数分間継続して発光させる光となる。
【００５０】
光学フィルタ９０５は、発光部９０１からの光や外乱光を除去し、植物１からのクロロフィル蛍光のみを取り出すためのフィルタである。光センサ９０６は、光学フィルタ９０５を透過したクロロフィル蛍光の強度を、電気的な信号に変換する。
【００５１】
測定光分離電気的フィルタ９０７は、光センサ９０６により変換された電気信号のDC成分を除去して、測定光に対するクロロフィル蛍光の強度を取り出すためのフィルタである。サンプリング部９０８は、中央制御部９００からの制御に従い、測定光分離電気的フィルタ９０７によるフィルタ処理で得られる電気信号を、アナログ信号からデジタル信号に変換する。
【００５２】
光合成量子収率演算部９０９は、中央制御部９００からの制御に従い、サンプリング部９０８によりサンプリングされたデジタル信号から、光合成の量子収率を算出する。結果表示部９１０は、中央制御部９００からの制御に従い、光合成量子収率演算部９０９により算出された光合成の量子収率の演算結果を表示する。ここで、光合成の量子収率としては、Fv/Fmと、ΦPSIIが表示される。
【００５３】
光合成計測装置９０は、以上のように構成される。
【００５４】
（発光強度とクロロフィル蛍光強度との関係）
ここで、図５の光合成計測装置９０においては、植物１に対して照射する光の種類を変えながら、植物１のクロロフィル蛍光のサンプリングが行われるが、植物１に照射される光の発光強度と、クロロフィル蛍光の強度との関係を表すと、図６に示すような関係となる。
【００５５】
図６において、図６のＡは、発光部９０１から発光される光の強度を示し、図６のＢは、発光部９０１から発光される光を制御することで引き起こされる、植物１のクロロフィル蛍光の強度を示している。また、図６においては、横方向が時間を表し、時間の方向は、図中の左側から右側に向かう方向とされる。
【００５６】
図６のＡにおいては、発光部９０１から発光される光の種類（強度）に応じて、測定光のみが発光される区間と、飽和光と測定光（飽和光＋測定光）が発光される区間と、定常光と測定光（定常光＋測定光）が発光される区間がある。
【００５７】
ここで、光合成の量子収率として、Fv/Fmの計測を行う場合、当該計測の開始直後の測定光のみを発光している区間で、測定光によるクロロフィル蛍光（F0）が得られ、飽和光＋測定光を発光している区間で、飽和光によるクロロフィル蛍光（Fm）が得られる。
【００５８】
そして、このような発光シーケンスで得られるクロロフィル蛍光（F0，Fm）を用い、下記の式（１）を演算することで、光合成の量子収率として、Fv/Fmが求められる。
【００５９】
Fv/Fm ＝ (Fm - F0) / Fm       ・・・（１）
【００６０】
また、光合成の量子収率として、ΦPSIIの計測を行う場合、植物１の光合成が定常状態時の定常光＋測定光を発光している区間で、定常光によるクロロフィル蛍光（Fs）が得られ、光合成が定常状態時の飽和光＋測定光を発光している区間で、飽和光によるクロロフィル蛍光（Fm'）が得られる。
【００６１】
そして、このような発光シーケンスで得られるクロロフィル蛍光（Fs，Fm'）を用い、下記の式（２）を演算することで、光合成の量子収率として、ΦPSIIが求められる。
【００６２】
ΦPSII ＝ (Fm' - Fs) / Fs       ・・・（２）
【００６３】
このように、図５の光合成計測装置９０では、短時間（例えば、2〜3秒）で電子伝達系を詰まらせるために、非常に強い光である飽和光を与えている。なお、計測の最後には、定常光の照射が停止され、測定光によるクロロフィル蛍光（F0'）が得られる。
【００６４】
（電気的なモデル）
図７は、植物１の光合成の量子収率（Fv/Fm）を算出するための計算アルゴリズムを電気的なモデルで表した図である。
【００６５】
図７においては、入力光の大きさをIと表現し、クロロフィル蛍光、熱、光化学反応としてエネルギーが吸収される度合いを、それぞれ、GF，GH，GPと表現した場合に、電気回路的には、Iが電流に相当し、GF，GH，GPが並列に接続されたコンダクタンスに相当することを表している。
【００６６】
特に、反応中心（RC：Reaction Center）がオープン（Open）の状態であるとき、クロロフィル蛍光（F0）は、電気的な等価回路として考えると、図７のＡに示した等価回路で表すことができるので、下記の式（３）のように表すことができる。
【００６７】
【数１】
2018056102.tif 000003
・・・（３）
【００６８】
一方で、反応中心（RC）がクローズ（Close）の状態であるとき、クロロフィル蛍光（Fm）は、電気的な等価回路として考えると、図７のＢに示した等価回路で表すことができるので、下記の式（４）のように表すことができる。
【００６９】
【数２】
2018056102.tif 000004
・・・（４）
【００７０】
そして、上述の式（３）と、式（４）から、光合成の量子収率（Fv/Fm）は、下記の式（５）のように表すことができる。
【００７１】
【数３】
2018056102.tif 000005
・・・（５）
【００７２】
ただし、式（５）において、右辺は、入力光（I）のうち、反応中心（RC）に分割される比率、すなわち、収率を示している。
【００７３】
（パルス変調蛍光測定に関する補足）
図８は、図５の光合成計測装置９０により行われるパルス変調蛍光測定に関する補足を説明するための図である。
【００７４】
光合成計測装置９０においては、クロロフィル蛍光（F0，Fm）を測定するために、入力光（I）として、測定光と、飽和光＋測定光が用いられる。ここで、図８のＡに示すように、クロロフィル蛍光（Fm）の測定では、交流の測定光に加えて、飽和光が印加されるため、クロロフィル蛍光の応答も、飽和光に対する応答（DC成分）が重ね合わされたかたちで、測定される。
【００７５】
しかしながら、図８のＢに示すように、測定光分離電気的フィルタ９０７の働きにより、飽和光に対するクロロフィル蛍光の応答（DC成分）は、除去されるため、飽和光は、反応中心（RC）をクローズの状態にさせることで、上述の図７のＢに示した電気的なモデルで表現可能な状態にする効果にのみ働いて、サンプリング部９０８による測定光及びクロロフィル蛍光のサンプリング結果には、影響しないことになる。よって、光合成の量子収率（Fv/Fm）は、上述したアルゴリズムにより演算することができる。
【００７６】
以上のようにして、図５の光合成計測装置９０では、パルス変調蛍光測定が行われる。しかしながら、このパルス変調蛍光測定では、太陽光の強度を超えるような非常に強い人工光である飽和光を用いるため、現実的には、限られた領域での測定は可能であるが、例えば、圃場の植物を対象とした場合のような、飽和光を当てることが困難となる広範囲の領域での測定には向いておらず、現時点では実用化には至っていない。
【００７７】
そのため、飽和光を用いずに、光合成の量子収率を求めて、植物の内部の状態を判断（診断）できるようにすることが望まれる。そこで、本技術では、植物の光合成を示す過渡モデルを１又は複数用意し、対象の過渡モデルに応じた計算アルゴリズムを、測定対象の植物から得られるクロロフィル蛍光の過渡応答の測定データに適用することで、当該過渡モデルに関するパラメータのうち、未知のパラメータである光合成の量子収率が推定されるようにする。これにより、飽和光を用いずに、光合成の量子収率を算出することが可能となって、植物の内部の状態を判断（診断）することができる。
【００７８】
（過渡モデルと計算アルゴリズムによるパラメータ推定）
図９は、本技術の過渡モデルと光合成量子収率計算アルゴリズムを用いたパラメータ推定を概念的に表した図である。
【００７９】
図９において、入力生成Ｃ１は、測定対象植物Ｃ３に対し、例えば測定光や育成光などの入力光Ｃ２が照射されることを表している。また、入力生成Ｃ１は、過渡モデルＣ５が生成され、あらかじめ用意されることを表している。
【００８０】
クロロフィル蛍光出力Ｃ４は、入力光Ｃ２が照射された測定対象植物Ｃ３をサンプリングすることで得られる、測定対象植物Ｃ３のクロロフィル蛍光の過渡応答の測定データである。このクロロフィル蛍光出力Ｃ４が、光合成量子収率計算アルゴリズムＣ６の入力データとなる。
【００８１】
過渡モデルＣ５は、植物の光合成での光化学反応の電子伝達系を表現したモデルである。過渡モデルＣ５としては、例えば、パドルモデル（Puddle Model）、レイクモデル（Lake Model）、又はコネクテッドユニットモデル（Connected Unit Model）等の植物生理学上のモデルを用いることができる。なお、過渡モデルＣ５は、１又は複数用意することができる。
【００８２】
ここで、過渡モデルＣ５のモデリングの手法としては、グレーボックスモデルが用いられる。グレーボックスモデルとは、ホワイトボックスモデル（第一原理モデル）とブラックボックスモデル（システム同定）の中間に位置づけられるモデルであって、システムの構造の一部が不明なモデルである。
【００８３】
例えば、グレーボックスモデルでは、対象の第一原理モデルを基本として、未知のパラメータを、入出力データ（入力データ）から推定する。このグレーボックスに関する技術としては、例えば、特開2011-137627号に開示されている技術がある。この文献には、建造物における熱伝達を、グレーボックスモデルを用いてモデル化する手法が開示されている。
【００８４】
光合成量子収率計算アルゴリズムＣ６は、過渡モデルＣ５を用い、クロロフィル蛍光出力Ｃ４から、光合成の量子収率を算出するためのアルゴリズムである。光合成量子収率計算アルゴリズムＣ６は、例えば、過渡モデルＣ５の種類や、算出すべき光合成の量子収率（Fv/Fm，ΦPSII）などに応じて、１又は複数用意することができる。
【００８５】
ここで、過渡モデルＣ５は、グレーボックスモデルに相当し、複数のパラメータのうち、未知のパラメータとして、測定対象植物Ｃ３の光合成の量子収率を含んでいる。そのため、クロロフィル蛍光出力Ｃ４（測定対象植物Ｃ３のクロロフィル蛍光の過渡応答の測定データ）を入力データとして、この過渡モデルＣ５に応じた光合成量子収率計算アルゴリズムＣ６を適用することで、未知のパラメータとしての測定対象植物Ｃ３の光合成の量子収率を算出することが可能となる。
【００８６】
光合成量子収率演算結果Ｃ７は、光合成量子収率計算アルゴリズムＣ６の演算結果として、測定対象植物Ｃ３の光合成の量子収率（Fv/Fm，ΦPSII）が算出されることを表している。
【００８７】
以上のように、本技術では、クロロフィル蛍光出力Ｃ４に対し、過渡モデルＣ５に応じた光合成量子収率計算アルゴリズムＣ６を適用することで、過渡モデルＣ５に関するパラメータのうち、未知のパラメータとして、測定対象植物Ｃ３の光合成の量子収率を推定（算出）することになる。以下、このような過渡モデルＣ５と光合成量子収率計算アルゴリズムＣ６を用いたパラメータ推定を実現するための具体的な実施の形態として、第１の実施の形態乃至第５の実施の形態について説明する。
【００８８】
＜２．第１の実施の形態＞
【００８９】
まず、図１０乃至図１８を参照して、第１の実施の形態について説明する。
【００９０】
（光合成計測装置の構成例）
図１０は、第１の実施の形態の光合成計測装置１０の構成例を示す図である。
【００９１】
光合成計測装置１０は、測定対象の植物１から得られるクロロフィル蛍光の過渡応答の測定データに対し、植物の光合成を示す過渡モデルに応じた計算アルゴリズムを適用して、当該過渡モデルに関するパラメータのうち、未知のパラメータとして、光合成の量子収率を推定（算出）する。
【００９２】
図１０において、光合成計測装置１０は、中央制御部１００、発光部１０１、発光制御部１０２、光学フィルタ１０３、光センサ１０４、サンプリング部１０５、過渡モデル保持部１０６、計算アルゴリズム保持部１０７、セレクタ１０８、セレクタ１０９、光合成量子収率演算部１１０、及び結果表示部１１１から構成される。
【００９３】
中央制御部１００は、例えば、CPU(Central Processing Unit)やFPGA(Field Programmable Gate Array)等の回路から構成される。中央制御部１００は、光合成計測装置１０の各部の動作を制御する。例えば、中央制御部１００は、発光制御部１０２による発光制御、サンプリング部１０５によるサンプリング、セレクタ１０８による過渡モデルの選択、セレクタ１０９による計算アルゴリズムの選択、及び結果表示部１１１による結果の表示を制御する。
【００９４】
発光部１０１は、例えば、LED(Light Emitting Diode)による光を発光（出射）する発光装置である。発光部１０１は、発光制御部１０２からの制御に従い、例えば３種類の強度の光を発光する。発光部１０１から発光された光は、植物１に照射される。発光制御部１０２は、中央制御部１００からの制御に従い、発光部１０１から発光される光の発光強度と発光シーケンスを制御する。
【００９５】
光学フィルタ１０３は、環境光とクロロフィル蛍光を分離するためのフィルタである。例えば、人工的な励起光を植物１に与える場合には、その励起光（可視光）と、クロロフィル蛍光（680〜750nm前後の波長）とを分離する波長のバンドパスフィルタ又はローパスフィルタを使用する。また、太陽光の下で測定を行う場合には、太陽光が環境光となるので、太陽暗線（例えば、酸素（O2）、吸収帯687nm）の蛍光を選択的に透過させる、波長の狭いバンドパスフィルタを使用してもよい。
【００９６】
光センサ１０４は、クロロフィル蛍光の強度を測定（センシング）するために、クロロフィル蛍光の波長に感度を有する１次元又は２次元のセンサである。光センサ１０４は、光学フィルタ１０３を透過した光を、センシング素子により検出し、その結果得られるクロロフィル蛍光の測定信号（強度値）を、サンプリング部１０５に供給する。
【００９７】
ここで、光センサ１０４としては、１次元のセンシングを目的とする場合には、例えば、フォトダイオードが搭載されたセンシング素子を用いることができる。また、２次元のセンシングを目的とする場合には、光センサ１０４として、例えば、CMOS(Complementary Metal-Oxide Semiconductor)やCCD(Charge Coupled Device)など、複数の画素が繰り返しパターンで２次元配列されたセンシング素子から構成されるイメージセンサを用いることができる。
【００９８】
なお、センシングとは、測定対象の植物１を測定することを意味する。また、センシングには、測定対象の植物１を撮像するという意味も含まれる。
【００９９】
サンプリング部１０５は、中央制御部１００からの制御に従い、光センサ１０４からのクロロフィル蛍光の測定信号（強度値）に対するサンプリングを行うことで、クロロフィル蛍光の測定信号を、時間経過の順に、アナログ信号からデジタル信号に変換する。当該サンプリングの結果得られたクロロフィル蛍光の測定データは、光合成量子収率演算部１１０に供給される。
【０１００】
なお、このサンプリングで得られる測定データは、１次元又は２次元の時系列のクロロフィル蛍光の測定データとされる。また、測定データは、数値データのほか、例えば、画像データなど、各種の形式のデータとすることができる。
【０１０１】
過渡モデル保持部１０６は、１又は複数の過渡モデルを保持する。セレクタ１０８は、中央制御部１００からの制御に従い、過渡モデル保持部１０６に保持されている過渡モデルの中から、対象の過渡モデルを選択し、光合成量子収率演算部１１０に供給する。
【０１０２】
計算アルゴリズム保持部１０７は、１又は複数の光合成量子収率計算アルゴリズムを保持する。セレクタ１０９は、中央制御部１００からの制御に従い、計算アルゴリズム保持部１０７に保持されている光合成量子収率計算アルゴリズムの中から、対象の光合成量子収率計算アルゴリズムを選択し、光合成量子収率演算部１１０に供給する。
【０１０３】
なお、過渡モデル保持部１０６及び計算アルゴリズム保持部１０７は、半導体メモリ（例えばRAM(Random Access Memory)等）などから構成される記憶部１３０として構成される。
【０１０４】
光合成量子収率演算部１１０には、サンプリング部１０５からのクロロフィル蛍光の測定データと、セレクタ１０８からの過渡モデルと、セレクタ１０９からの光合成量子収率計算アルゴリズムが供給される。ただし、サンプリング部１０５からの測定データには、植物１に対して照射される光を変化させることで得られる、植物１のクロロフィル蛍光の過渡応答の測定データが含まれる。
【０１０５】
光合成量子収率演算部１１０は、植物１のクロロフィル蛍光の過渡応答の測定データに対し、植物の光合成を示す過渡モデルに応じた光合成量子収率計算アルゴリズムを適用することで、当該過渡モデルに関するパラメータのうち、未知のパラメータである光合成の量子収率を算出する。光合成量子収率演算部１１０は、算出した光合成の量子収率を、結果表示部１１１に供給する。
【０１０６】
この光合成の量子収率としては、例えば、光合成の電子伝達系の最大量子収率（Fv/Fm）及び光合成の電子伝達系の量子収率（ΦPSII）のうち少なくとも一方が算出される。
【０１０７】
結果表示部１１１は、例えば、LCD(Liquid Crystal Display)やOELD(Organic Electroluminescence Display)などのディスプレイにより構成される。結果表示部１１１は、中央制御部１００からの制御に従い、光合成量子収率演算部１１０から供給される、光合成の量子収率（Fv/Fm，ΦPSII）の演算結果に関するデータ（例えば、数値データや画像データ）を表示する。
【０１０８】
なお、図１０において、サンプリング部１０５、セレクタ１０８、セレクタ１０９、及び、光合成量子収率演算部１１０の機能の全部又は一部は、例えば、CPUとして構成される中央制御部１００が実行するプログラム（ソフトウェア）により実現されるようにすることができる。
【０１０９】
また、図１０においては、結果表示部１１１と記憶部１３０は、光合成計測装置１０の内部に設けられるとして説明したが、表示装置や記憶装置として、光合成計測装置１０の外部に設けられるようにしてもよい。この場合、光合成量子収率演算部１１０は、ネットワークを介して外部の記憶装置から、過渡モデルや計算アルゴリズムを取得することになる。また、中央制御部１００は、光合成量子収率演算部１１０により算出された光合成の量子収率の演算結果のデータを、外部の表示装置に表示させたり、あるいは外部の記憶装置に記憶させたりすることができる。
【０１１０】
光合成計測装置１０は、以上のように構成される。
【０１１１】
（発光強度とクロロフィル蛍光強度との関係）
ここで、図１０の光合成計測装置１０においては、植物１に対して照射する光の種類を変えながら、植物１のクロロフィル蛍光のサンプリングが行われるが、植物１に照射される光の発光強度と、クロロフィル蛍光の強度との関係を表すと、図１１に示すような関係となる。
【０１１２】
ただし、図１１においては、縦方向が強度のレベルを表し、図中の上側に行くほど、そのレベルが大きくなることを表している。また、図１１においては、横方向が時間を表し、時間の方向は、図中の左側から右側に向かう方向とされる。
【０１１３】
図１１において、図１１のＡは、発光部１０１から発光される光の強度を示し、図１１のＢは、発光部１０１から発光される光を制御することにより引き起こされる、植物１のクロロフィル蛍光の強度を示している。
【０１１４】
図１１のＡにおいては、発光部１０１から発光される光の種類（強度）に応じて３つの区間がある。第１区間は、微弱な光である測定光が発光される区間である。この測定光としては、例えば、10[μmol/m2/s]前後の光が発光される。また、第２区間は、第１育成光が発光される区間であり、第３区間は、第１育成光よりも強度が強い光である第２育成光が発光される区間である。育成光は、植物の育成の際に使用される強度の光であり、例えば、最低で、100[μmol/m2/s]程度とされる。ただし、第２育成光のほうが、第１育成光よりも強度が強くなる。
【０１１５】
ここで、光合成の量子収率として、Fv/Fmの計測を行う場合、一般的に、植物１を、暗状態に適切な時間（例えば30分以上）おいて、暗順化させた後に、計測を開始する。そして、当該計測が開始されると、第１区間で測定光を照射し、その後に続く第２区間で、第１育成光を照射することで、光の強度を変化させて、クロロフィル蛍光の測定データを取得する。
【０１１６】
すなわち、図１１のＡに示すように、第１区間で測定光が照射され、第１区間の次の第２区間で第１育成光が照射された場合、図１１のＢに示すように、クロロフィル蛍光の強度は、微弱な光である測定光を照射する第１区間では、レベルがほとんど変化しないが、第２区間で第１育成光が照射された直後に、そのレベルが大きくなっている。
【０１１７】
例えば、光合成計測装置１０においては、サンプリング部１０５によって、第１区間（又は第１区間と第２区間の境界）における時刻t0と、第２区間における、第１育成光が照射された直後のクロロフィル蛍光の強度のレベルが上昇したとき（時刻t1，時刻t2，時刻t3）に、クロロフィル蛍光のサンプリングが行われるようにする。なお、時刻t1，時刻t2，時刻t3は、時刻t0を起点として、最低250マイクロ秒で、かつ、最高500ミリ秒の条件を満たすようにする。
【０１１８】
このサンプリングによって、植物１に対して照射される光を変化（測定光から第１育成光に変化）させることで得られる、f0の値と、f(t)（t ＝ t1，t2，t3）の値を、植物１のクロロフィル蛍光の過渡応答の測定データとして取得することができる。そして、光合成量子収率演算部１１０は、このクロロフィル蛍光の過渡応答の測定データ（クロロフィル蛍光の強弱の値）に基づいて、Fv/Fmを算出することができる。つまり、ここでは、光合成の暗順応状態の量子収率（Fv/Fm）が求められる。
【０１１９】
また、光合成の量子収率として、ΦPSIIの計測を行う場合、一般的に、植物１を、明状態に適切な時間（例えば5分以上）おいておく必要がある。そのため、第２区間の間は、第１育成光を照射し続けて、その後、第３区間において、第２育成光を照射することで、光の強度を変化させて、クロロフィル蛍光の測定データを取得する。
【０１２０】
すなわち、図１１のＡに示すように、第２区間で第１育成光が照射され続け、第２区間の次の第３区間で第２育成光が照射された場合、図１１のＢに示すように、クロロフィル蛍光の強度は、第２区間ではレベルがほとんど変化しないが、第３区間で第２育成光が照射された直後に、そのレベルが大きくなっている。
【０１２１】
例えば、光合成計測装置１０においては、サンプリング部１０５によって、第２区間（又は第２区間と第３区間の境界）における時刻ts0と、第３区間における、第２育成光が照射された直後のクロロフィル蛍光の強度のレベルが上昇したとき（時刻ts1，時刻ts2，時刻ts3）に、クロロフィル蛍光のサンプリングが行われるようにする。なお、時刻ts1，時刻ts2，時刻ts3は、時刻ts0を起点として、最低250マイクロ秒で、かつ、最高500ミリ秒の条件を満たすようにする。
【０１２２】
このサンプリングによって、植物１に対して照射される光を変化（第１育成光から第２育成光に変化）させることで得られる、fsの値と、f(t)（t ＝ ts1，ts2，ts3）の値を、植物１のクロロフィル蛍光の過渡応答の測定データとして取得することができる。そして、光合成量子収率演算部１１０は、植物１のクロロフィル蛍光の過渡応答の測定データ（クロロフィル蛍光の強弱の値）に基づいて、ΦPSIIを算出することができる。つまり、ここでは、一定の強度の光が当たって光合成が行われているときの量子収率（ΦPSII）が求められる。
【０１２３】
なお、図１１の発光シーケンスにおいては、光の発光強度を大きく（強く）することで、植物１のクロロフィル蛍光の過渡応答の測定データが取得されるようにしているが、要は、光の発光強度に変化が起きればよく、例えば、光の発光強度を小さく（弱く）することでも、過渡応答の測定データは取得できる。
【０１２４】
（光合成計測処理）
次に、図１２のフローチャートを参照して、図１０の中央制御部１００が、光合成計測装置１０の各部を制御することにより実現される光合成計測処理の流れについて説明する。
【０１２５】
なお、図１２の光合成計測処理を行うに際し、光合成計測装置１０による計測可能な位置に、高等植物である植物１が存在するものとする。また、発光部１０１による発光のタイミングや、サンプリング部１０５によるサンプリングのタイミングなどについては、適宜、上述した図１１を参照しながら説明する。
【０１２６】
ステップＳ１０１において、中央制御部１００は、各部を制御して、植物１に対し、測定光を照射して、図１１の時刻t0のタイミングでのクロロフィル蛍光のサンプリングが行われるようにする。
【０１２７】
具体的には、中央制御部１００は、発光制御部１０２を制御して、植物１に対し、発光部１０１からの測定光が照射されるようにする。これにより、サンプリング部１０５には、光センサ１０４から、測定光が照射された植物１からのクロロフィル蛍光の測定信号が入力される。
【０１２８】
また、中央制御部１００は、サンプリング部１０５を制御して、光センサ１０４からのクロロフィル蛍光の測定信号のサンプリングが行われるようにする。これにより、図１１の時刻t0において、測定光が照射された植物１のクロロフィル蛍光の測定信号が、デジタル信号に変換され、その結果得られる測定データが、光合成量子収率演算部１１０に出力される。
【０１２９】
ステップＳ１０２において、中央制御部１００は、各部を制御して、植物１に対し、第１育成光を照射して、図１１の時刻t1，時刻t2，時刻t3のタイミングでのクロロフィル蛍光のサンプリングが行われるようにする。
【０１３０】
具体的には、中央制御部１００は、発光制御部１０２を制御して、植物１に対し、発光部１０１からの第１育成光が照射されるようにする。これにより、サンプリング部１０５には、光センサ１０４から、第１育成光が照射された植物１からのクロロフィル蛍光の測定信号が入力される。
【０１３１】
また、中央制御部１００は、サンプリング部１０５を制御して、光センサ１０４からのクロロフィル蛍光の測定信号のサンプリングが行われるようにする。これにより、図１１の時刻t1において、第１育成光が照射された植物１のクロロフィル蛍光の測定信号が、デジタル信号に変換され、その結果得られる測定データが、光合成量子収率演算部１１０に出力される。図１１の時刻t2，時刻t3のタイミングにおいても、時刻t1のタイミングと同様に、サンプリングが行われる。
【０１３２】
このようにして、光合成の量子収率（Fv/Fm）を、過渡モデルに対応した光合成量子収率計算アルゴリズムにより算出する場合に、入力データとして必要となる、図１１のf0の値と、f(t)（t ＝ t1，t2，t3）の値が、植物１のクロロフィル蛍光の過渡応答の測定データとして得られる。
【０１３３】
ステップＳ１０３において、中央制御部１００は、発光制御部１０２を制御して、植物１に対し、発光部１０１からの第１育成光の照射が継続されるようにする。すなわち、ステップＳ１０２の処理で、植物１に対し、第１育成光の照射が開始されたが、ステップＳ１０３の処理でも、植物１に対する第１育成光の照射が継続される。
【０１３４】
ステップＳ１０４において、中央制御部１００は、サンプリング部１０５を制御して、光センサ１０４からのクロロフィル蛍光の測定信号のサンプリングが行われるようにする。これにより、図１１の時刻ts0において、第１育成光が照射された植物１のクロロフィル蛍光の測定信号が、デジタル信号に変換され、その結果得られる測定データが、光合成量子収率演算部１１０に出力される。
【０１３５】
ステップＳ１０５において、中央制御部１００は、各部を制御して、植物１に対し、第２育成光を照射して、図１１の時刻ts1，時刻ts2，時刻ts3のタイミングでのクロロフィル蛍光のサンプリングが行われるようにする。
【０１３６】
具体的には、中央制御部１００は、発光制御部１０２を制御して、植物１に対し、発光部１０１からの第２育成光が照射されるようにする。これにより、サンプリング部１０５には、光センサ１０４から、第２育成光が照射された植物１からのクロロフィル蛍光の測定信号が入力される。
【０１３７】
また、中央制御部１００は、サンプリング部１０５を制御して、光センサ１０４からのクロロフィル蛍光の測定信号のサンプリングを行う。これにより、図１１の時刻ts1において、第２育成光が照射された植物１のクロロフィル蛍光の測定信号が、デジタル信号に変換され、その結果得られる測定データが、光合成量子収率演算部１１０に出力される。また、図１１の時刻ts2，時刻ts3のタイミングにおいても、時刻ts1のタイミングと同様に、サンプリングが行われる。
【０１３８】
このようにして、光合成の量子収率（ΦPSII）を、過渡モデルに対応した光合成量子収率計算アルゴリズムにより算出する場合に、入力データとして必要となる、図１１のfsの値と、f(t)（t ＝ ts1，ts2，ts3）の値が、植物１のクロロフィル蛍光の過渡応答の測定データとして得られる。
【０１３９】
ステップＳ１０６において、中央制御部１００は、セレクタ１０８を制御して、過渡モデル保持部１０６に保持されている１又は複数の過渡モデルの中から、光合成の量子収率の計測に用いられる過渡モデルが選択されるようにする。ここでは、例えば、パドルモデル、レイクモデル、及びコネクテッドユニットモデルの過渡モデルが保持されている場合に、パドルモデル（の過渡モデル）を選択することができる。
【０１４０】
なお、過渡モデルの詳細については、図１３及び図１４を参照して後述する。
【０１４１】
ステップＳ１０７において、中央制御部１００は、セレクタ１０９を制御して、計算アルゴリズム保持部１０７に保持されている１又は複数の光合成量子収率計算アルゴリズムの中から、ステップＳ１０６の処理で選択された過渡モデルに対応する光合成量子収率計算アルゴリズムが選択されるようにする。ここでは、例えば、過渡モデルとしてパドルモデルが選択された場合に、パドルモデル用の光合成量子収率計算アルゴリズムを選択することができる。
【０１４２】
ステップＳ１０８において、中央制御部１００は、光合成量子収率演算部１１０を制御して、ステップＳ１０１乃至Ｓ１０５の処理で得られた過渡応答の測定データに対し、ステップＳ１０７の処理で得られた光合成量子収率計算アルゴリズムを適用することで、光合成の量子収率が算出されるようにする。
【０１４３】
ここでは、例えば、ステップＳ１０１乃至Ｓ１０２の処理で得られた植物１のクロロフィル蛍光の過渡応答の測定データ（例えば、図１１のf0の値と、f(t)（t ＝ t1，t2，t3）の値）を入力データとして、パドルモデル用の光合成量子収率計算アルゴリズムに従い、処理することで、光合成の量子収率（Fv/Fm）を算出することができる。
【０１４４】
また、例えば、ステップＳ１０３乃至Ｓ１０５の処理で得られた植物１のクロロフィル蛍光の過渡応答の測定データ（例えば、図１１のfsの値と、f(t)（t ＝ ts1，ts2，ts3）の値）を入力データとして、パドルモデル用の光合成量子収率計算アルゴリズムに従い、処理することで、光合成の量子収率（ΦPSII）を算出することができる。
【０１４５】
なお、クロロフィル応答の過渡応答の測定データの測定から、光合成量子収率計算アルゴリズムの適用までの流れについては、図１５乃至図１７を参照して後述する。
【０１４６】
ステップＳ１０９において、中央制御部１００は、結果表示部１１１を制御して、ステップＳ１０８の処理の演算結果が、結果表示部１１１に表示されるようにする。これにより、結果表示部１１１には、光合成の量子収率（Fv/Fm，ΦPSII）の演算結果に関するデータが表示される。
【０１４７】
以上、光合成計測処理の流れについて説明した。
【０１４８】
なお、図１２の光合成計測処理では、光合成の量子収率として、Fv/FmとΦPSIIの両方を算出する場合を説明したが、光合成の量子収率としては、Fv/Fm及びΦPSIIの少なくとも一方が算出されればよく、ステップＳ１０８の処理で算出される光合成の量子収率の内容に応じて、ステップＳ１０１乃至Ｓ１０５の処理のサンプリングの内容も異なる。
【０１４９】
（過渡モデルの構成）
次に、図１３及び図１４を参照して、過渡モデルの構成について説明する。
【０１５０】
（反応中心（RC）のモデル）
過渡モデルとしては、例えば、パドルモデルやレイクモデル、コネクテッドユニットモデル等の植物生理学上のモデルを用いることができる。図１３のＡには、パドルモデルに対応した反応中心（RC）のモデルの構成の例が示されている。
【０１５１】
図１３のＡに示すように、パドルモデルは、入力光の大きさ（PPFD値）をIで表した場合に、時間的な要素を考慮すれば、反応中心（RC）が、入力された光子を優先的に取得し（光子を電子に変換し）、反応中心（RC）が取得できなかった光子が、クロロフィル蛍光（Fluorescence）と熱（Heat）側に流れ込むような等価回路により表すことができる。なお、反応中心（RC）は、複数のモジュールから構成される（QA，QB，PQ等）。
【０１５２】
すなわち、図１３のＡの等価回路において、反応中心（RC）は、キャパシタとして表され、優先的に入力された光子に相当する電流（i(t)）を、電荷（q(t)）として蓄積する。また、反応中心（RC）側に、光子に相当する電流（i(t)）が流れ込むことで、反応中心（RC）のキャパシタの最大電荷保持量を超えた場合には、光子に相当する電流（I-i(t)）が、蛍光（Fluorescence）と熱（Heat）側に流れ込むことになる。
【０１５３】
ここで、図１３のＡに示した反応中心（RC）の非線形な挙動を、確率演算を用いて記述する。すなわち、反応中心（RC）への吸収は、入力とオープン率との乗算により確率的に決まるとして、その特性を下記の式（６）のように定義する。
【０１５４】
【数４】
2018056102.tif 000006
・・・（６）
【０１５５】
ただし、式（６）において、右辺の第１項である（I×PM）は、入力を表し、右辺の第２項である(QM-q(t))/QMは、オープン率を表している。また、PMは、オープン時（q(t)＝0）の収率を表し、QMは、反応中心（RC）の最大電荷保持量を表している。
【０１５６】
ここで、式（６）について、i(t)を、q(t)で表現すると、下記の式（７）のように表すことができる。
【０１５７】
【数５】
2018056102.tif 000007
・・・（７）
【０１５８】
そして、式（７）を解くと、下記の式（８）のように表すことができる。
【０１５９】
【数６】
2018056102.tif 000008
・・・（８）
【０１６０】
なお、上述した説明では、図１３のＡに示した等価回路のような、反応中心（RC）が単体の場合のパドルモデルを一例に説明したが、反応中心（RC）が複数存在する場合も想定される。ここでは、反応中心（RC）が複数存在する場合（一定の葉面積当たりのモデルを考える場合）、下記の式（９）のように表すことができる。
【０１６１】
【数７】
2018056102.tif 000009
・・・（９）
【０１６２】
ただし、式（９）において、Nは、反応中心（RC）の個数を表している。また、パドルモデルでは、反応中心（RC）の間のエネルギーを共有しないモデルであるから、式（８）と式（９）とは、数式形状に違いがないため、以下の説明では、反応中心（RC）が複数存在する場合も、式（８）で代用するものとする。
【０１６３】
ここで、図１４には、本技術の発明者による過渡モデルのシミュレーションの結果の例が示されている。図１４では、論理最高性能を持つ光合成が行われた場合におけるモデルごとのクロロフィル応答の波形を表しており、横軸は、時間（単位：秒）とされ、縦軸は、蛍光の強度を表している。
【０１６４】
図１４において、６つのクロロフィル応答の波形のうち、「0」で表された波形は、パドルモデルに対応している。パドルモデルの場合には、反応中心（RC）の間でのエネルギー伝達がなく、クロロフィル応答の波形は、自然対数の形状となる。
【０１６５】
また、「0.2」，「0.4」，「0.6」，「0.8」で表された４つの波形は、コネクテッドユニットモデルに対応している。コネクテッドユニットモデルの場合には、反応中心（RC）の間でのエネルギー伝達が一定の確率で起こり、クロロフィル応答の波形は、自然対数の形状とシグモイド関数の形状が合わさった形状となる。
【０１６６】
さらに、「1」で表された波形は、レイクモデルに対応している。レイクモデルの場合には、反応中心（RC）の間のエネルギー伝達が、100％の確率で起こり、クロロフィル応答の波形は、シグモイド関数の形状となる。
【０１６７】
（クロロフィル蛍光のモデル）
図１３に戻り、図１３のＢには、パドルモデルに対応したクロロフィル蛍光のモデルの構成の例が示されている。
【０１６８】
図１３のＢに示すように、パドルモデルでは、反応中心（RC）に伝達できなかったエネルギーが、クロロフィル蛍光（Fluorescence）と熱（Heat）となって放出されるとして、モデル化することができる。すなわち、図１３のＢの等価回路では、クロロフィル蛍光と熱の側に入力される入力光の大きさ（PPFD値）をI-i(t)と表現し、クロロフィル蛍光、熱としてエネルギーが吸収される度合いを、それぞれ、GF，GHと表現した場合に、電気回路的には、I-i(t)が電流に相当し、GF，GHが並列に接続されたコンダクタンスに相当することを表している。
【０１６９】
ここで、図１３のＢに示したクロロフィル蛍光の応答（F(t)）は、下記の式（１０）により表すことができる。
【０１７０】
【数８】
2018056102.tif 000010
・・・（１０）
【０１７１】
ここで、先に述べたクロロフィル蛍光の測定（図７等）におけるF0，Fmを、このモデルで表現すると、下記の式（１１）と式（１２）で表すことができる。
【０１７２】
【数９】
2018056102.tif 000011
・・・（１１）
【０１７３】
【数１０】
2018056102.tif 000012
・・・（１２）
【０１７４】
ここで、上述した式（１）より、Fv/Fm ＝ (Fm - F0) / Fmであるから、これを、式（１０）と式（１１）に代入すると、下記の式（１３）が求められる。
【０１７５】
Fv/Fm ＝ PM       ・・・（１３）
【０１７６】
すなわち、この式（１３）から、このクロロフィル蛍光のモデルにおいて、飽和光を用い、Fv/Fmを求めることと、何らか別の手法で、PMを求めることは等価であると言える。
【０１７７】
以上のように、通常の電気回路では、上述したような、反応中心（RC）、クロロフィル蛍光、熱という３つの出力先が存在する場合には、エネルギーは、並列回路として３つの出力素子のパラメータに互いに影響されながら分割される。しかしながら、光合成においては、「励起されたエネルギーの抜き取り時間」という要素があり、反応中心（RC）のほうが、クロロフィル蛍光と熱よりも「抜き取り時間」が早いため、この時間的な優先度によって、反応中心（RC）に、優先的にエネルギーが分配される。そのため、上述のパドルモデルでは、反応中心（RC）で受け取れなかった残りのエネルギーが、独立回路としてクロロフィル蛍光と熱に分配されるものとしている。
【０１７８】
なお、上述した説明では、過渡モデルとして、パドルモデルを用いた場合を例示したが、レイクモデル又はコネクテッドユニットモデルを用いるようにしてもよい。また、過渡モデルは、パドルモデルやレイクモデル、コネクテッドユニットモデルに限らず、例えば、植物の葉が大きなダメージを受けたときのモデルであるエクセスモデル（Excess Model）など、植物生理学上の他のモデルを用いることができる。さらに、これらのモデルは、一般的には非線形であるため、本技術の実施の形態で記載しているような解析解を用いたパラメータ同定アルゴリズムではなく、回帰分析を使ったパラメータ同定アルゴリズムを用いてもよい。
【０１７９】
（光合成量子収率計算アルゴリズムの詳細：解析解を用いた手法例）
次に、図１５乃至図１７を参照して、光合成量子収率計算アルゴリズムの詳細について説明する。ここでは、まず、光合成の量子収率として、Fv/Fmを算出する場合の計算アルゴリズムを先に説明し、その後、ΦPSIIを算出する場合についても説明する。ただし、ここでは、上述したパドルモデルに対応した過渡モデルが用いられる場合における、光合成量子収率計算アルゴリズムを説明するものとする。
【０１８０】
（１）Fv/Fmの計算アルゴリズム
【０１８１】
（Fv/Fmの計算アルゴリズム）
rFvを、F(t)とF0を用い、下記の式（１４）で表すように定義する。
【０１８２】
【数１１】
2018056102.tif 000013
・・・（１４）
【０１８３】
ここで、上述した式（１０）に対し、上述した式（１１）と、式（８）を代入して、整理すると、式（１４）は、下記の式（１５）のように表すことができる。
【０１８４】
【数１２】
2018056102.tif 000014
・・・（１５）
【０１８５】
そして、式（１５）を整理すると、下記の式（１６）のように表すことができる。
【０１８６】
【数１３】
2018056102.tif 000015
・・・（１６）
【０１８７】
すなわち、Fv/Fmの計算アルゴリズムでは、「rFv」及び「I×PM/QM」を計測し、式（１６）で表される計算アルゴリズム（パドルモデル用の光合成量子収率計算アルゴリズム）に従い、PMを算出することになる。つまり、上述したように、飽和光を用いてFv/Fmを求めることと、何らか別の手法で、PMを求めることは等価であるのだから、ここでは、このパドルモデル用の光合成量子収率計算アルゴリズムに従い、PMを算出することで、光合成の量子収率としてのFv/Fmを求めている。
【０１８８】
（時定数計算アルゴリズム）
ここで、上述した式（１６）における時定数（k：I×PM/QM）を算出するためのアルゴリズムについて説明する。
【０１８９】
一般的に、R,C,L等価回路の過渡解は、自然対数を用いて表現することができる。ここでは、仮に、解を、下記の式（１７）で表した場合を想定する。
【０１９０】
【数１４】
2018056102.tif 000016
・・・（１７）
【０１９１】
そして、この式（１７）により解が表される場合に、式（１７）に対し、１階微分と２階微分を適用し、その結果得られる値の比をとることで、時定数（k）を求めることができる。すなわち、ここでの１階微分と２階微分は、下記の式（１８）と式（１９）により表される。
【０１９２】
【数１５】
2018056102.tif 000017
・・・（１８）
【０１９３】
【数１６】
2018056102.tif 000018
・・・（１９）
【０１９４】
また、式（１８）で表される１階微分と、式（１９）で表される２階微分との比によって、下記の式（２０）で表すように、時定数（k）が求められる。
【０１９５】
【数１７】
2018056102.tif 000019
・・・（２０）
【０１９６】
ここで、式（２０）の値は、３点の測定データから求めることができる。例えば、図１５に示すように、測定データがf(t)で表される場合に、時刻t1，時刻t2，時刻t3における測定データが、a，b，cであったとき、時刻t1〜t2のグラフの傾きと、時刻t2〜t3のグラフの傾きは、次のように表される。
【０１９７】
すなわち、時刻t1から時刻t2までの時間、及び、時刻t2から時刻t3までの時間が、Δtで表される場合、時刻t1〜t2のグラフの傾きは、b-a/Δt＝mで表され、時刻t2〜t3のグラフの傾きは、c-b/Δt＝nで表される。このとき、式（１９）で表される２階微分は、n-m/Δt＝lと表すことができるので、上述した式（２０）は、下記の式（２１）のように表すことができる。
【０１９８】
【数１８】
2018056102.tif 000020
・・・（２１）
【０１９９】
このようにして、３点の測定データから、時定数（k）を求めることができる。
【０２００】
（計算アルゴリズムの適用までの流れ）
次に、図１６を参照して、Fv/Fmを算出する場合における、測定データの測定から、計算アルゴリズムの適用までの流れについて説明する。
【０２０１】
光合成計測装置１０（図１０）では、図１６のステップＳ１５１乃至Ｓ１５４の処理を実行することで、植物１のクロロフィル蛍光の過渡応答の測定データの測定から、パドルモデル用の光合成量子収率計算アルゴリズムの適用までの処理を行う。なお、ここでは、光合成の量子収率として、Fv/Fmを算出するので、ステップＳ１５１に示した発光シーケンスは、上述した図１１の発光シーケンスにおける第１区間乃至第２区間の一部の区間に対応している。
【０２０２】
まず、サンプリング部１０５が、時刻t0，時刻t1，時刻t2，時刻t3のタイミングでのサンプリングを行うことで、f0の値と、f(t)（t ＝ t1，t2，t3）の値が、植物１のクロロフィル蛍光の過渡応答の測定データとして得られる（Ｓ１５１）。
【０２０３】
次に、光合成量子収率演算部１１０は、上述した図１５に示した手法を用い、サンプリング部１０５からの過渡応答の測定データ（３点の測定データ）を、式（２１）に適用することで、時定数（k）を算出することができる（Ｓ１５２）。
【０２０４】
次に、光合成量子収率演算部１１０は、サンプリング部１０５からの過渡応答の測定データを、式（１４）に適用することで、rFvを算出することができる（Ｓ１５３）。
【０２０５】
そして、光合成量子収率演算部１１０は、ステップＳ１５２の処理で得られた時定数（k）と、ステップＳ１５３の処理で得られたrFvを、パドルモデル用の光合成量子収率計算アルゴリズムとしての下記の式（２２）に適用することで、PMを算出することができる（Ｓ１５４）。このようにして求められるPMが、光合成の量子収率としてのFv/Fmに相当することは、先に述べた通りある。
【０２０６】
【数１９】
2018056102.tif 000021
・・・（２２）
【０２０７】
このようにして、光合成の量子収率としてのFv/Fmを求めることができる。なお、この例においては、f0の値と、f(t)（t ＝ t1，t2，t3）の値をサンプリングしているが、原理的には、時刻t0におけるf0と、時刻t1におけるf(t1)との値を共用することも可能である。
【０２０８】
（２）ΦPSIIの計算アルゴリズム
【０２０９】
光合成の量子収率として、ΦPSIIを算出する場合には、ΦPSIIの計算アルゴリズムが用いられる。また、ΦPSIIを算出する場合にも、時定数（ｋ）を求める必要があるが、上述した図１５に示した手法により、時定数（ｋ）を算出することができる。
【０２１０】
（計算アルゴリズムの適用までの流れ）
ここで、図１７を参照して、ΦPSIIを算出する場合における、測定データの測定から、計算アルゴリズム（パドルモデル用の光合成量子収率計算アルゴリズム）の適用までの流れについて説明する。
【０２１１】
光合成計測装置１０（図１０）では、図１７のステップＳ１７１乃至Ｓ１７４の処理を実行することで、植物１のクロロフィル蛍光の過渡応答の測定データの測定から、パドルモデル用の光合成量子収率計算アルゴリズムの適用までの処理を行う。なお、ここでは、光合成の量子収率として、ΦPSIIを算出するので、ステップＳ１７１に示した発光シーケンスは、上述した図１１の発光シーケンスにおける第２区間乃至第３区間の一部の区間に対応している。
【０２１２】
まず、サンプリング部１０５が、時刻ts0，時刻ts1，時刻ts2，時刻ts3のタイミングでのサンプリングを行うことで、fsの値と、f(t)（t ＝ ts1，ts2，ts3）の値が、植物１のクロロフィル蛍光の過渡応答の測定データとして得られる（Ｓ１７１）。
【０２１３】
次に、光合成量子収率演算部１１０は、上述した図１５に示した手法を用い、サンプリング部１０５からの過渡応答の測定データ（３点の測定データ）に基づいて、時定数（k）を算出することができる（Ｓ１７２）。
【０２１４】
次に、光合成量子収率演算部１１０は、サンプリング部１０５からの過渡応答の測定データを基づいて、rFvを算出することができる（Ｓ１７３）。
【０２１５】
そして、光合成量子収率演算部１１０は、ステップＳ１７２の処理で得られた時定数（k）と、ステップＳ１７３の処理で得られたrFvを、パドルモデル用の光合成量子収率計算アルゴリズムとしての下記の式（２４）に適用することで、ΦPSIIを算出することができる（Ｓ１７４）。
【０２１６】
【数２０】
2018056102.tif 000022
・・・（２４）
【０２１７】
このようにして、光合成の量子収率としてのΦPSIIを求めることができる。
【０２１８】
以上のように、第１の実施の形態では、植物１に対し、測定光と育成光（定常光）を照射することで、クロロフィル蛍光の過渡応答の測定データを取得し、当該測定データを入力データとして、パドルモデル等の過渡モデルに応じた光合成量子収率計算アルゴリズムに従い、処理することで、飽和光がなくても、光合成の量子収率を算出することができる。
【０２１９】
ここで、現状の技術（パルス変調蛍光測定）を採用している光合成計測装置９０（図５）と、本技術を採用している光合成計測装置１０（図１０）での測定に使用される光を比較すれば、図１８に示すようになる。すなわち、図１８においては、現状と本技術の方式ごとに、計測対象としての光合成の量子収率（Fv/Fm，ΦPSII）を計測する際の測定に使用される光が、強度別に表されている。
【０２２０】
現状の方式では、上述した図６の発光シーケンスに示したように、Fv/Fmを計測する場合には、測定光と飽和光の２種類の光が用いられる。ここで用いられる測定光の強度は、10[μmol/m2/s]前後とされる。また、飽和光の強度は、2000[μmol/m2/s]以上とされる。飽和光は、飽和パルス光とも称される。なお、現状の方式で、Fv/Fmを計測する場合、定常光は不要となる。
【０２２１】
また、現状の方式で、ΦPSIIを計測する場合には、定常光と飽和光の２種類の光が用いられる。ここで用いられる定常光の強度は、植物１の育成に使う強度であって、最低で100[μmol/m2/s]程度とされる。また、ここでも、飽和光の強度は、2000[μmol/m2/s]以上とされる。なお、現状の方式で、ΦPSIIを計測する場合、測定光は不要となる。
【０２２２】
一方で、本技術の方式では、上述した図１１の発光シーケンスに示したように、Fv/Fmを計測する場合には、測定光と育成光（第１育成光）の２種類の光が用いられる。ここで用いられる測定光の強度は、10[μmol/m2/s]前後とされる。また、第１育成光の強度は、植物１の育成に使う強度であって、最低で100[μmol/m2/s]程度とされる。第１育成光は、第１定常光とも称される。なお、本技術の方式で、Fv/Fmを計測する場合、飽和光は不要となる。
【０２２３】
また、本技術の方式で、ΦPSIIを計測する場合には、２種類の育成光（第１育成光、第２育成光）が用いられる。ここで用いられる、第１育成光と第２育成光の強度は、植物１の育成に使う強度であって、最低で100[μmol/m2/s]程度とされる。ただし、第２育成光の強度は、第１育成光の強度よりも大きくなる（第１育成光 < 第２育成光）。第２育成光は、第２定常光とも称される。なお、本技術の方式で、ΦPSIIを計測する場合、測定光と飽和光は不要となる。
【０２２４】
このように、現状の方式においては、Fv/Fmの計測とΦPSIIの計測でともに飽和光が必要になる。一方で、本技術の方式においては、Fv/Fmの計測とΦPSIIの計測でともに飽和光が不要になる。そのため、本技術を採用している光合成計測装置１０（図１０）では、現状の技術（パルス変調蛍光測定）を採用している光合成計測装置９０（図５）と比べて、測定に必要な発光装置（発光部）の単位面積当たりの出力が小さくて済むため、より広範囲の測定が可能となる。
【０２２５】
＜３．第２の実施の形態＞
【０２２６】
（光合成計測装置の構成例）
図１９は、第２の実施の形態の光合成計測装置２０の構成例を示す図である。
【０２２７】
図１９において、光合成計測装置２０は、中央制御部２００、レーザ発光部２０１、発光制御部２０２、光学フィルタ２０３、光センサ２０４、サンプリング部２０５、過渡モデル保持部２０６、計算アルゴリズム保持部２０７、セレクタ２０８、セレクタ２０９、光合成量子収率演算部２１０、結果表示部２１１、アクチュエータ２１２、及び画像合成部２１３から構成される。ただし、過渡モデル保持部２０６と、計算アルゴリズム保持部２０７は、半導体メモリ等の記憶部２３０として構成される。
【０２２８】
図１９の光合成計測装置２０は、図１０の光合成計測装置１０と比べて、発光部１０１の代わりに、レーザ発光部２０１が設けられ、アクチュエータ２１２と、画像合成部２１３が新たに設けられる点が異なっている。なお、光合成計測装置２０において、それ以外の構成については、光合成計測装置１０の構成と、基本的に同様とされるため、その説明は適宜省略するものとする。
【０２２９】
レーザ発光部２０１は、発光制御部２０２からの制御に従い、例えば３種類（測定光、第１育成光、及び第２育成光）の強度のレーザ光を発光（出射）する。アクチュエータ２１２は、レーザ発光部２０１から発光されるレーザ光によって、測定対象の植物１の所定の領域を走査する機能を有している。
【０２３０】
レーザ光が照射された植物１のクロロフィル蛍光は、光学フィルタ２０３を介して、２次元のセンサとして構成される光センサ２０４により検出される。そして、サンプリング部２０５によって、光センサ２０４からのクロロフィル蛍光の測定信号（強度値）のサンプリングが行われる。
【０２３１】
光合成量子収率演算部２１０は、植物１のクロロフィル蛍光の過渡応答の測定データに対し、植物の光合成を示す過渡モデルに応じた光合成量子収率計算アルゴリズムを適用することで、当該過渡モデルに関するパラメータのうち、未知のパラメータである光合成の量子収率（Fv/Fm，ΦPSII）を算出する。ただし、光合成量子収率演算部２１０では、レーザ光による走査ごとに、光合成の量子収率が算出され、画像合成部２１３に供給される。
【０２３２】
画像合成部２１３は、中央制御部２００からの制御に従い、レーザ光による走査ごとに得られる光合成の量子収率を合成して、１枚の画像を生成する。画像合成部２１３により生成された画像は、結果表示部２１１に表示される。
【０２３３】
光合成計測装置２０は、以上のように構成される。
【０２３４】
（光合成計測処理）
図１９の光合成計測装置２０により実行される光合成計測処理は、上述した光合成計測処理（図１２）と基本的に同様であるが、植物１に対して照射される光が、レーザ光である点が異なっている。すなわち、図１２のステップＳ１０１乃至Ｓ１０５の処理でサンプリングを行う際に、発光部１０１（図１０）からのLED等による光の代わりに、レーザ発光部２０１（図１９）からのレーザ光が用いられる。
【０２３５】
以上のように、第２の実施の形態では、植物１に対し、レーザ光による走査を行うことで、クロロフィル蛍光の過渡応答の測定データを取得し、当該測定データを入力データとして、パドルモデル等の過渡モデルに応じた光合成量子収率計算アルゴリズムに従い、処理することで、飽和光がなくても、光合成の量子収率を算出することができる。
【０２３６】
ここで、レーザ光を用いて、光合成の量子収率を計測する技術としては、特表2006-504956号に開示されている技術がある。この文献には、レーザのビームを動かしてドットごとにカメラでクロロフィル蛍光をスキャンし、光合成の量子収率の２次元計測を行う装置が開示されている。
【０２３７】
具体的には、この文献に開示されている技術では、植物材料を、Ffast計測と称される、約1Hz〜約10kHz（すなわち、ドットあたり1秒〜0.1ミリ秒）のレーザでスキャンし、その後に、Fslow計測と称される、約0.01〜約1Hz（すなわち、ドットあたり100秒〜1秒）のレーザでスキャンする。そして、光合成活性の量子効率（光合成のイメージング量子効率（IQP））を、下記の式（２５）により演算している。
【０２３８】
IQP ＝ (Fslow - Ffast) / Fslow       ・・・（２５）
【０２３９】
この文献に開示されている技術においては、1ドットの計測に1秒以上を要するため、解像度を上げようとした場合に計測時時間を要する。例えば、100m2（10m×10m）の圃場の計測を、縦横それぞれ1000ポイントを測定した場合、278時間を要するが、本技術を使用した場合、光センサ１０４として２次元配列されたセンシング素子を用いれば、一度に１画面分を計測できるため、最大で500ミリ秒以内に計測を終えることができる。
【０２４０】
また、この文献によれば、上記の式（２５）の計算式で、光合成活性の量子効率を求めるために、Fv/Fm＝(Fm-F0)/Fmという計算式が定義されているため、Fmの測定のために飽和光の強度を持つレーザを使用する必要があると考えられるが、本技術を使用した場合、レーザの強度を弱くできるため、レーザの安全規格の観点から優位性があると考えられる。
【０２４１】
一方で、第２の実施の形態の光合成計測装置３０（図１９）においても、植物１に照射するための光としてレーザ光を用いてはいるが、上記の文献に開示されている技術と比べて、飽和光を用いずに測定が可能であることから、レーザ発光部２０１の単位面積当たりの出力が小さくて済むため、野外での測定の際に安全性が向上するほか、測定にかかる時間を短縮することが可能となる。
【０２４２】
＜４．第３の実施の形態＞
【０２４３】
（光合成計測装置の構成例）
図２０は、第３の実施の形態の光合成計測装置３０の構成例を示す図である。
【０２４４】
図２０において、光合成計測装置３０は、中央制御部３００、発光制御部３０２、光学フィルタ３０３、光センサ３０４、サンプリング部３０５、過渡モデル保持部３０６、計算アルゴリズム保持部３０７、セレクタ３０８、セレクタ３０９、光合成量子収率演算部３１０、結果表示部３１１、及びインターフェース３１２から構成される。ただし、過渡モデル保持部３０６と、計算アルゴリズム保持部３０７は、半導体メモリ等の記憶部３３０として構成される。
【０２４５】
図２０の光合成計測装置３０は、図１０の光合成計測装置１０と比べて、発光部１０１の代わりに、外部の装置と接続するためのインターフェース３１２が設けられる点が異なっている。なお、光合成計測装置３０において、それ以外の構成については、光合成計測装置１０の構成と、基本的に同様とされるため、その説明は適宜省略するものとする。
【０２４６】
インターフェース３１２は、入出力インターフェース回路などにより構成され、外部の植物育成用発光装置３１と接続可能とされる。植物育成用発光装置３１は、例えばLEDによる光など、植物の育成のための光を発光する装置である。つまり、環境によっては、植物の日照が不足し、その生育に影響がでる場合があるので、植物育成用発光装置３１によって、一定時間LED等による光を照射することで、光合成に必要な波長の光が効率的に照射され、植物の成長を促進することができる。
【０２４７】
発光制御部３０２は、インターフェース３１２を介して、植物育成用発光装置３１から発光される光を制御する。植物育成用発光装置３１は、発光制御部３０２からの制御に従い、例えば３種類（測定光、第１育成光、及び第２育成光）の強度の光を発光（出射）する。
【０２４８】
植物育成用発光装置３１からの光が照射された植物のクロロフィル蛍光は、光学フィルタ３０３を介して、光センサ３０４により検出される。そして、サンプリング部３０５によって、光センサ３０４からのクロロフィル蛍光の測定信号（強度値）のサンプリングが行われる。
【０２４９】
光合成量子収率演算部３１０は、植物１のクロロフィル蛍光の過渡応答の測定データに対し、植物の光合成を示す過渡モデルに応じた光合成量子収率計算アルゴリズムを適用することで、当該過渡モデルに関するパラメータのうち、未知のパラメータである光合成の量子収率（Fv/Fm，ΦPSII）を算出する。光合成量子収率演算部３１０による演算結果は、結果表示部３１１に表示される。
【０２５０】
光合成計測装置３０は、以上のように構成される。
【０２５１】
（光合成計測処理）
図２０の光合成計測装置３０により実行される光合成計測処理は、上述した光合成計測処理（図１２）と基本的に同様であるが、植物１に対して照射される光が、外部の植物育成用発光装置３１により照射される点が異なっている。すなわち、図１２のステップＳ１０１乃至Ｓ１０５の処理でサンプリングを行う際に、発光部１０１（図１０）からのLED等による光の代わりに、植物育成用発光装置３１（図２０）からのLED等による光が用いられる。
【０２５２】
以上のように、第３の実施の形態では、植物１に対し、外部の植物育成用発光装置３１からの光を照射することで、クロロフィル蛍光の過渡応答の測定データを取得し、当該測定データを入力データとして、パドルモデル等の過渡モデルに応じた光合成量子収率計算アルゴリズムに従い、処理することで、飽和光がなくても、光合成の量子収率を算出することができる。
【０２５３】
また、第３の実施の形態では、植物１に照射する光が、外部の植物育成用発光装置３１から発光されるようにしているため、光合成計測装置３０側では、光合成計測を行うための特別な光源を用意することなく、光合成の量子収率を求めることが可能となる。
【０２５４】
＜５．第４の実施の形態＞
【０２５５】
（光合成計測装置の構成例）
図２１は、第４の実施の形態の光合成計測装置４０の構成例を示す図である。
【０２５６】
図２１において、光合成計測装置４０は、中央制御部４００、光学フィルタ４０３、光センサ４０４、サンプリング部４０５、過渡モデル保持部４０６、計算アルゴリズム保持部４０７、セレクタ４０８、セレクタ４０９、光合成量子収率演算部４１０、結果表示部４１１、環境光センサ４１２、サンプリング部４１３、及び最適変化検出部４１４から構成される。ただし、過渡モデル保持部４０６と、計算アルゴリズム保持部４０７は、半導体メモリ等の記憶部４３０として構成される。
【０２５７】
図２１の光合成計測装置４０は、図１０の光合成計測装置１０と比べて、発光部１０１及び発光制御部１０２の代わりに、環境光センサ４１２、サンプリング部４１３、及び最適変化検出部４１４が設けられる点が異なっている。なお、光合成計測装置４０において、それ以外の構成については、光合成計測装置１０の構成と、基本的に同様とされるため、その説明は適宜省略するものとする。
【０２５８】
環境光センサ４１２は、太陽光等の環境光の強度を測定するためのセンサである。環境光センサ４１２は、環境光の強度に応じた測定信号を、サンプリング部４１３に供給する。
【０２５９】
サンプリング部４１３は、環境光センサ４１２からの環境光の測定信号（強度値）に対するサンプリングを行うことで、環境光の測定信号を、時間経過の順に、アナログ信号からデジタル信号に変換する。当該サンプリングの結果得られた環境光の測定データは、最適変化検出部４１４に供給される。
【０２６０】
最適変化検出部４１４は、サンプリング部４１３からの環境光の測定データに基づいて、環境光の強度の変化のうち、植物１の光合成計測を行うのに必要な急激な変化を検出し、その検出結果を、中央制御部４００に供給する。
【０２６１】
すなわち、例えば、光合成の量子収率として、Fv/Fmを算出する場合には、図１１に示した発光シーケンスにおける、第１区間の測定光から、第２区間の第１育成光への発光強度の変化に相当する環境光の強度変化が検出される。また、例えば、光合成の量子収率として、ΦPSIIを算出する場合には、図１１に示した発光シーケンスにおける、第２区間の第１育成光から、第３区間の第２育成光への発光強度の変化に相当する環境光の強度変化が検出される。
【０２６２】
太陽光等の環境光が照射された植物のクロロフィル蛍光は、光学フィルタ４０３を介して、光センサ４０４により検出される。サンプリング部４０５は、中央制御部４００からの制御に従い、光センサ４０４からのクロロフィル蛍光の測定信号（強度値）のサンプリングを行う。ここで、中央制御部４００は、最適変化検出部４１４からの検出結果に基づいて、環境光の強度の変化に応じたサンプリングが行われるようにする。
【０２６３】
すなわち、例えば、光合成の量子収率として、Fv/Fmを算出する場合には、図１１に示した発光シーケンスにおける、発光強度の変化に相当する環境光の強度変化を検出することで、図１１の時刻t0，時刻t1，時刻t2，時刻t3のタイミングに相当する所定のタイミングで、クロロフィル蛍光のサンプリングが行われるようにする。これにより、図１１のf0の値とf(t)（t ＝ t1，t2，t3）の値に相当する所定の値が、植物１のクロロフィル蛍光の過渡応答の測定データとして得られる。
【０２６４】
また、例えば、光合成の量子収率として、ΦPSIIを算出する場合には、図１１に示した発光シーケンスにおける、発光強度の変化に相当する環境光の強度変化を検出することで、図１１の時刻ts0，時刻ts1，時刻ts2，時刻ts3のタイミングに相当する所定のタイミングで、クロロフィル蛍光のサンプリングが行われるようにする。これにより、図１１のfsの値とf(t)（t ＝ ts1，ts2，ts3）の値に相当する所定の値が、植物１のクロロフィル蛍光の過渡応答の測定データとして得られる。
【０２６５】
光合成量子収率演算部４１０は、植物１のクロロフィル蛍光の過渡応答の測定データに対し、植物の光合成を示す過渡モデルに応じた光合成量子収率計算アルゴリズムを適用することで、当該過渡モデルに関するパラメータのうち、未知のパラメータである光合成の量子収率（Fv/Fm，ΦPSII）を算出する。光合成量子収率演算部４１０による演算結果は、結果表示部４１１に表示される。
【０２６６】
以上、光合成計測装置４０は、以上のように構成される。
【０２６７】
（光合成計測処理）
図２１の光合成計測装置４０により実行される光合成計測処理は、上述した光合成計測処理（図１２）と基本的に同様であるが、植物１に対して照射される光が、太陽光等の環境光のみである点と、植物１のクロロフィル蛍光のサンプリングのタイミングが、環境光の強度変化のタイミングに応じて決定される点が異なっている。
【０２６８】
すなわち、図１２のステップＳ１０１乃至Ｓ１０５の処理でサンプリングを行う際に、発光部１０１（図１０）からのLED等による光の代わりに、太陽光等の環境光が用いられる。また、時刻t0〜t3や、時刻ts0〜ts3などのクロロフィル蛍光のサンプリングのタイミングが、発光部１０１（図１０）から発光される光の発光強度の変化ではなく、太陽光等の環境光の強度変化のタイミングに応じて決定される。
【０２６９】
以上のように、第４の実施の形態では、植物１に対し、太陽光等の環境光のみが照射されるようにして、当該環境光の強度変化に応じたタイミングでサンプリングされたクロロフィル蛍光の過渡応答の測定データを取得し、当該測定データを入力データとして、パドルモデル等の過渡モデルに応じた光合成量子収率計算アルゴリズムに従い、処理することで、光合成の量子収率を算出することができる。
【０２７０】
すなわち、上述した光合成計測装置１０（図１０）等では、植物１に対し照射される光の強度やタイミングを制御することで、植物１に照射される光が瞬間的に強まるようにし、そのタイミングでサンプリングを行っていた。一方で、第４の実施の形態の光合成計測装置４０では、太陽光等の環境光は、その強度変化を制御することはできないので、強度変化のタイミングを検出して、その検出結果に応じたタイミングで、サンプリングが行われるようにしている。
【０２７１】
例えば、環境光が太陽光である場合には、雲の通過などによって引き起こされる強度の変化を検出すればよい。換言すれば、第４の実施の形態の光合成計測装置４０では、太陽光等の環境光の変化を、環境光センサ４１２により継続して測定し、たまたま都合のよい環境光の強度変化が起こったときを選択して、光合成の計測が行われることになる。
【０２７２】
このように、第４の実施の形態では、飽和光だけでなく、測定光と定常光（育成光）をも用いずに、太陽光等の環境光の強度変化を用いて、光合成の量子収率を算出している。すなわち、太陽光等の環境光の強度変化などを利用する場合には、図１１に示した発光シーケンスは、必須とはならない。
【０２７３】
なお、第４の実施の形態においては、環境光の強度の変化を測定する手法として、環境光センサ４１２により環境光の強度を測定する手法を説明したが、例えば、植物１の葉面照度を測定する手法など、他の手法によって、環境光の強度の変化が測定されるようにしてもよい。
【０２７４】
＜６．第５の実施の形態＞
【０２７５】
（情報処理装置の構成例）
図２２は、第５の実施の形態の情報処理装置５０の構成例を示す図である。
【０２７６】
図２２において、情報処理装置５０は、中央制御部５００、過渡モデル保持部５０６、計算アルゴリズム保持部５０７、セレクタ５０８、セレクタ５０９、光合成量子収率演算部５１０、及び結果表示部５１１から構成される。ただし、過渡モデル保持部５０６と、計算アルゴリズム保持部５０７は、半導体メモリ等の記憶部５３０として構成される。
【０２７７】
図２２の情報処理装置５０は、図１０の光合成計測装置１０と比べて、発光部１０１乃至サンプリング部１０５の代わりに、インターフェース５１２が設けられる点が異なっている。なお、情報処理装置５０において、それ以外の構成については、光合成計測装置１０の構成と、基本的に同様とされるため、その説明は適宜省略するものとする。
【０２７８】
インターフェース５１２は、入出力インターフェース回路などにより構成され、外部ストレージ５１、又はネットワーク５３を介して外部クロロフィル蛍光計測装置５２と接続可能とされる。
【０２７９】
外部ストレージ５１は、例えば、ハードディスクや半導体メモリなどの大容量の記録装置である。外部ストレージ５１は、植物１のクロロフィル蛍光の過渡応答の測定データのファイル（以下、過渡応答測定ファイルという）を蓄積している。外部ストレージ５１は、情報処理装置５０からの要求などに応じて、蓄積している過渡応答測定ファイルを、情報処理装置５０に提供する。情報処理装置５０においては、外部ストレージ５１からの過渡応答測定ファイルが、インターフェース５１２介して、光合成量子収率演算部５１０に供給される。
【０２８０】
外部クロロフィル蛍光計測装置５２は、植物１のクロロフィル蛍光を計測（測定）する装置である。例えば、外部クロロフィル蛍光計測装置５２は、図１０に示した中央制御部１００、発光部１０１、発光制御部１０２、光学フィルタ１０３、光センサ１０４、及びサンプリング部１０５に相当する機能を有し、測定光や育成光が照射された植物１のクロロフィル蛍光のサンプリングを行うことが可能である。
【０２８１】
外部クロロフィル蛍光計測装置５２は、情報処理装置５０からの要求などに応じて、植物１のクロロフィル蛍光の過渡応答の測定データを、ネットワーク５３を介して、情報処理装置５０に提供する。情報処理装置５０においては、外部クロロフィル蛍光計測装置５２からの過渡応答の測定データが、インターフェース５１２介して、光合成量子収率演算部５１０に供給される。
【０２８２】
なお、ここでは、過渡応答の測定データの代わりに、当該測定データのファイル（過渡応答測定ファイル）が提供されるようにしてもよい。また、例えば、インターネット等のネットワーク５３に接続されたサーバ等の情報処理装置が、過渡応答の測定データや過渡応答測定ファイルを、ネットワーク５３を介して提供するようにしてもよい。
【０２８３】
情報処理装置５０において、光合成量子収率演算部５１０には、外部ストレージ５１又は外部クロロフィル蛍光計測装置５２から、過渡応答測定ファイル又は過渡応答の測定データが提供される。ただし、過渡応答測定ファイルには、過渡応答の測定データが含まれる。
【０２８４】
光合成量子収率演算部５１０は、外部から提供された過渡応答の測定データに対し、植物の光合成を示す過渡モデルに応じた光合成量子収率計算アルゴリズムを適用することで、当該過渡モデルに関するパラメータのうち、未知のパラメータである光合成の量子収率（Fv/Fm，ΦPSII）を算出する。光合成量子収率演算部５１０による演算結果は、結果表示部５１１に表示される。
【０２８５】
情報処理装置５０は、以上のように構成される。
【０２８６】
（情報処理）
図２２の情報処理装置５０により実行される情報処理は、上述した光合成計測処理（図１２）と比べると、前段のステップＳ１０１乃至Ｓ１０５の処理が実行されずに、後段のステップＳ１０６乃至Ｓ１０９の処理のみが実行される点が異なっている。
【０２８７】
すなわち、前段のステップＳ１０１乃至Ｓ１０５の処理は、外部クロロフィル蛍光計測装置５２により実行されるか、あるいは、前段のステップＳ１０１乃至Ｓ１０５の処理で得られた過渡応答測定ファイルが、外部ストレージ５１に格納されている。そして、情報処理装置５０では、外部ストレージ５１又は外部クロロフィル蛍光計測装置５２から得られた過渡応答の測定データを用いて、後段のステップＳ１０６乃至Ｓ１０９の処理が実行されることになる。
【０２８８】
以上のように、第５の実施の形態では、外部から提供された過渡応答の測定データを取得し、当該測定データを入力データとして、パドルモデル等の過渡モデルに応じた光合成量子収率計算アルゴリズムに従い、処理することで、光合成の量子収率を算出することができる。
【０２８９】
＜７．変形例＞
【０２９０】
（植物の測定装置の具体例）
図２３には、上述した第１の実施の形態乃至第５の実施の形態における、測定対象の植物１の測定を行う測定装置の具体例として、定点観測を行う定点測定装置６０Ａと、移動観測を行う移動測定装置６０Ｂを例示している。
【０２９１】
図２３のＡに示した定点測定装置６０Ａは、固定脚６１Ａによって、測定対象の植物１（例えば圃場の植物など）を測定（センシング）することが可能な位置に固定され、そこで測定された植物１のクロロフィル蛍光の過渡応答の測定データを、例えば無線通信等を利用して、情報処理装置５０（図２２）に送信する。情報処理装置５０は、定点測定装置６０Ａから送信されてくる過渡応答の測定データを処理することで、定点測定装置６０Ａにより定点観測された植物１の光合成の量子収率を求めることができる。
【０２９２】
図２３のＢに示した移動測定装置６０Ｂは、例えば無人航空機（UAV：Unmanned Aerial Vehicle）であって、プロペラ状の回転翼６１Ｂが回転することで飛行し、上空から、測定対象の植物１（例えば、圃場の植物など）を測定（空撮）する。移動測定装置６０Ｂは、そこで測定された植物１のクロロフィル蛍光の過渡応答の測定データを、例えば無線通信等を利用して、例えば、情報処理装置５０（図２２）に送信する。情報処理装置５０は、移動測定装置６０Ｂから送信されてくる過渡応答の測定データを処理することで、移動測定装置６０Ｂにより移動観測された植物１の光合成の量子収率を求めることができる。
【０２９３】
なお、移動測定装置６０Ｂは、無線操縦のほか、例えば、飛行ルートを座標データとしてあらかじめ記憶しておくことで、GPS(Global Positioning System)などの位置情報を用いて自律飛行するようにしてもよい。また、図２３のＢでは、移動測定装置６０Ｂが、回転翼６１Ｂを有する回転翼機であるとして説明したが、移動測定装置６０Ｂは、固定翼機であってもよい。
【０２９４】
また、上述した説明では、定点測定装置６０Ａ又は移動測定装置６０Ｂにより測定された過渡応答の測定データが、情報処理装置５０（図２２）に送信され、情報処理装置５０によって、植物１の光合成の量子収率を求めるとして説明したが、定点測定装置６０Ａ又は移動測定装置６０Ｂが、光合成計測装置１０（図１０）や情報処理装置５０（図２２）などと同様の機能を有することで、植物１の光合成の量子収率を求めるようにしてもよい。
【０２９５】
（植物以外の測定対象）
上述した説明では、測定対象を植物とした場合を説明したが、植物に限らず、生体全般を測定対象とすることができる。すなわち、本技術は、上述した図９に示したように、測定対象から得られる過渡応答の測定データを入力として、過渡モデルに応じた計算アルゴリズムを適用することで、当該過渡モデルに関するパラメータのうち、未知のパラメータを推定（算出）することを特徴としており、この原理は、生体全般に適用することができる。
【０２９６】
したがって、本技術によれば、測定対象の生体から得られる過渡応答の測定データを入力として、対象の生体が有する機能を示す過渡モデルに応じた計算アルゴリズムを適用することで、当該過渡モデルに関するパラメータのうち、未知のパラメータを算出することができる。
【０２９７】
＜８．コンピュータの構成＞
【０２９８】
上述した一連の処理（例えば、図１２に示した光合成計測処理の少なくとも一部の処理）は、ハードウェアにより実行することもできるし、ソフトウェアにより実行することもできる。一連の処理をソフトウェアにより実行する場合には、そのソフトウェアを構成するプログラムが、コンピュータにインストールされる。図２４は、上述した一連の処理をプログラムにより実行するコンピュータのハードウェアの構成例を示す図である。
【０２９９】
コンピュータ１０００において、CPU(Central Processing Unit)１００１、ROM(Read Only Memory)１００２、RAM(Random Access Memory)１００３は、バス１００４により相互に接続されている。バス１００４には、さらに、入出力インターフェース１００５が接続されている。入出力インターフェース１００５には、入力部１００６、出力部１００７、記録部１００８、通信部１００９、及び、ドライブ１０１０が接続されている。
【０３００】
入力部１００６は、キーボード、マウス、マイクロフォンなどよりなる。出力部１００７は、ディスプレイ、スピーカなどよりなる。記録部１００８は、ハードディスクや不揮発性のメモリなどよりなる。通信部１００９は、ネットワークインターフェースなどよりなる。ドライブ１０１０は、磁気ディスク、光ディスク、光磁気ディスク、又は半導体メモリなどのリムーバブル記憶媒体１０１１を駆動する。
【０３０１】
以上のように構成されるコンピュータ１０００では、CPU１００１が、ROM１００２や記録部１００８に記録されているプログラムを、入出力インターフェース１００５及びバス１００４を介して、RAM１００３にロードして実行することにより、上述した一連の処理が行われる。
【０３０２】
コンピュータ１０００(CPU１００１)が実行するプログラムは、例えば、パッケージメディア等としてのリムーバブル記憶媒体１０１１に記録して提供することができる。また、プログラムは、ローカルエリアネットワーク、インターネット、デジタル衛星放送といった、有線又は無線の伝送媒体を介して提供することができる。
【０３０３】
コンピュータ１０００では、プログラムは、リムーバブル記憶媒体１０１１をドライブ１０１０に装着することにより、入出力インターフェース１００５を介して、記録部１００８にインストールすることができる。また、プログラムは、有線又は無線の伝送媒体を介して、通信部１００９で受信し、記録部１００８にインストールすることができる。その他、プログラムは、ROM１００２や記録部１００８に、あらかじめインストールしておくことができる。
【０３０４】
ここで、本明細書において、コンピュータがプログラムに従って行う処理は、必ずしもフローチャートとして記載された順序に沿って時系列に行われる必要はない。すなわち、コンピュータがプログラムに従って行う処理は、並列的あるいは個別に実行される処理(例えば、並列処理あるいはオブジェクトによる処理)も含む。また、プログラムは、１のコンピュータ(プロセッサ)により処理されるものであってもよいし、複数のコンピュータによって分散処理されるものであってもよい。
【０３０５】
なお、本技術の実施の形態は、上述した実施の形態に限定されるものではなく、本技術の要旨を逸脱しない範囲において種々の変更が可能である。例えば、上述した複数の実施の形態の全て又は一部を組み合わせた形態を採用することができる。
【０３０６】
なお、本技術は、以下のような構成をとることができる。
【０３０７】
（１）
測定対象の生体から得られる過渡応答の測定データに対し、前記生体が有する機能を示す過渡モデルに応じた計算アルゴリズムを適用して、前記過渡モデルに関するパラメータのうち、未知のパラメータを算出する算出部を備える
情報処理装置。
（２）
前記生体は、植物であり、
前記過渡モデルは、前記植物の光合成での光化学反応の電子伝達系を表現したモデルであり、
前記計算アルゴリズムは、前記過渡モデルを用い、前記植物のクロロフィル蛍光の過渡応答から光合成の量子収率を算出するためのアルゴリズムであり、
前記測定データは、前記植物のクロロフィル蛍光の過渡応答の測定データであり、
前記算出部は、前記測定データを入力とし、前記計算アルゴリズムを適用することで、未知のパラメータとして、前記植物の光合成の量子収率を算出する
（１）に記載の情報処理装置。
（３）
前記植物の光合成の量子収率は、光合成の電子伝達系の最大量子収率（Fv/Fm）及び光合成の電子伝達系の量子収率（ΦPSII）のうち少なくとも一方である
（２）に記載の情報処理装置。
（４）
前記過渡モデルは、１又は複数用意され、
前記計算アルゴリズムは、前記過渡モデルに応じて１又は複数用意される
（１）乃至（３）のいずれかに記載の情報処理装置。
（５）
前記過渡モデルと前記計算アルゴリズムを記憶する記憶部をさらに備える
（１）乃至（４）のいずれかに記載の情報処理装置。
（６）
前記植物に対する光を発光する発光部を、前記測定データが得られるように制御する発光制御部をさらに備える
（２）乃至（５）のいずれかに記載の情報処理装置。
（７）
前記発光部と、
前記植物のクロロフィル蛍光をセンシングするセンサと
をさらに備える（６）に記載の情報処理装置。
（８）
レーザ光を発光する前記発光部と、
前記レーザ光を走査するための走査部と、
前記植物のクロロフィル蛍光をセンシングするセンサと
をさらに備える（６）に記載の情報処理装置。
（９）
前記発光部は、前記植物の育成用の光を発光する外部の発光装置であり、
前記発光装置と接続されるインターフェースと、
前記植物のクロロフィル蛍光をセンシングするセンサと
をさらに備え、
前記発光制御部は、前記インターフェースを介して、前記発光装置を、前記測定データが得られるように制御する
（６）に記載の情報処理装置。
（１０）
前記植物のクロロフィル蛍光をセンシングする第１のセンサと、
前記第１のセンサによるセンシング時に、環境光をセンシングする第２のセンサと
をさらに備え、
前記算出部は、前記環境光の変化時に得られた前記測定データに基づいて、前記植物の光合成の量子収率を算出する
（２）に記載の情報処理装置。
（１１）
前記環境光は、太陽光である
（１０）に記載の情報処理装置。
（１２）
前記測定データを提供する外部の処理装置と接続されるインターフェースをさらに備え、
前記算出部は、前記インターフェースを介して、前記処理装置から提供される前記測定データに基づいて、前記植物の光合成の量子収率を算出する
（２）に記載の情報処理装置。
（１３）
前記処理装置は、前記測定データを計測する計測装置、又は前記測定データを記録する記録装置である
（１２）に記載の情報処理装置。
（１４）
情報処理装置の情報処理方法において、
前記情報処理装置が、
測定対象の生体から得られる過渡応答の測定データに対し、前記生体が有する機能を示す過渡モデルに応じた計算アルゴリズムを適用して、前記過渡モデルに関するパラメータのうち、未知のパラメータを算出する
ステップを含む情報処理方法。
（１５）
コンピュータを、
測定対象の生体から得られる過渡応答の測定データに対し、前記生体が有する機能を示す過渡モデルに応じた計算アルゴリズムを適用して、前記過渡モデルに関するパラメータのうち、未知のパラメータを算出する算出部を備える情報処理装置として
機能させるためのプログラム。
（１６）
測定対象の生体をセンシングするセンサと、
前記センサによるセンシングによって前記生体から得られる過渡応答の測定データに対し、前記生体が有する機能を示す過渡モデルに応じた計算アルゴリズムを適用して、前記過渡モデルに関するパラメータのうち、未知のパラメータを算出する算出部と
を備えるセンシング装置。
【符号の説明】
【０３０８】
１  植物，  １０，２０，３０，４０  光合成計測装置，  ３１  植物育成用発光装置，  ５０  情報処理装置，  ５１  外部ストレージ，  ５２  外部クロロフィル蛍光計測装置，  ５３  ネットワーク，  １００，２００，３００，４００，５００  中央制御部，  １０１  発光部，  １０２，２０２，３０２  発光制御部，  １０３，２０３，３０３，４０３  光学フィルタ，  １０４，２０４，３０４，４０４  光センサ，  １０５，２０５，３０５，４０５  サンプリング部，  １０６，２０６，３０６，４０６，５０６  過渡モデル保持部，  １０７，２０７，３０７，４０７，５０７  計算アルゴリズム保持部，  １０８，２０８，３０８，４０８，５０８  セレクタ，  １０９，２０９，３０９，４０９，５０９  セレクタ，  １１０，２１０，３１０，４１０，５１０  光合成量子収率演算部，  １１１，２１１，３１１，４１１，５１１  結果表示部，  １３０，２３０，３３０，４３０，５３０  記憶部，  ２０１  レーザ発光部，  ２１２  アクチュエータ，  ２１３  画像合成部，  ３１２  インターフェース，  ４１２  環境光センサ，  ４１３  サンプリング部，  ４１４  最適変化検出部，  ５１２  インターフェース，  １０００  コンピュータ，  １００１  CPU
【図１】
2018056102.tif 000023
【図２】
2018056102.tif 000024
【図３】
2018056102.tif 000025
【図４】
2018056102.tif 000026
【図５】
2018056102.tif 000027
【図６】
2018056102.tif 000028
【図７】
2018056102.tif 000029
【図８】
2018056102.tif 000030
【図９】
2018056102.tif 000031
【図１０】
2018056102.tif 000032
【図１１】
2018056102.tif 000033
【図１２】
2018056102.tif 000034
【図１３】
2018056102.tif 000035
【図１４】
2018056102.tif 000036
【図１５】
2018056102.tif 000037
【図１６】
2018056102.tif 000038
【図１７】
2018056102.tif 000039
【図１８】
2018056102.tif 000040
【図１９】
2018056102.tif 000041
【図２０】
2018056102.tif 000042
【図２１】
2018056102.tif 000043
【図２２】
2018056102.tif 000044
【図２３】
2018056102.tif 000045
【図２４】
2018056102.tif 000046
【国際調査報告】
2018056102.tif 000047
2018056102.tif 000048
2018056102.tif 000049
2018056102.tif 000050