(57)【特許請求の範囲】
【請求項１】
植物の葉片の成長を測定するための方法であって、
ａ）撮影システムを較正し、撮影システムの較正パラメータを決定するステップと、
ｂ）前記較正された撮影システムを用いて、１枚の葉片の少なくとも１対のステレオ画像を撮像するステップと、
ｃ）撮影された各画像を処理するステップであって、
ｉ）しきい値セグメント化によって前記葉片をセグメント化することと、
ｉｉ）前記セグメント化された画像に対して形態学的侵食ステップを行うことと、
ｉｉｉ）前記形態学的侵食を施された画像上で、保管容器の縁部を除去するためにエッジ補正を行い、前記葉片の縁部を認識することと
を含むステップと、
ｄ）前記エッジ補正を施された画像から、前記葉片の投影される葉面積値Ａｐを算出するステップと、
ｅ）３Ｄ再構成を行うステップであって、前記撮影された１対の画像から、ステレオアルゴリズムによって視差画像Ｄを生成し、事前に決定された較正パラメータに基づいて、視差画像から３Ｄ表面モデルを計算するステップと、
ｆ）前記取得された３Ｄ表面モデルを平滑化するステップと、
ｇ）前記平滑化された表面モデルから、前記葉片の３Ｄ面積値Ａ３Ｄを算出するステップと、
ｈ）異なる時点で、前記ステップｂ）による撮像を繰り返すステップと、
ｉ）前記ステップｃ）〜ｇ）に従って異なる時点に撮影された撮像の解析によって、前記算出された面積値の時系列を生成するステップと、
ｊ）前記算出された面積値の前記時系列から、前記葉片の成長率を算出するステップと
を含む方法。
【請求項２】
撮影システムの較正時に、少なくとも３０枚の市松模様の画像が様々な向きから撮影され、カメラパラメータ、すなわち焦点距離、光軸と像面の交点、ならびに半径方向および接線方向歪みが決定されることを特徴とする請求項１に記載の方法。
【請求項３】
各撮影位置での撮像時に、２枚の画像が撮影され、基線長ｂが、左右の画像内で前記葉片を見ることができるような大きさに選択されることを特徴とする請求項１または２に記載の方法。
【請求項４】
しきい値セグメント化時に、最小および最大しきい値の使用によってバイナリマスク（Ｍｔ）が生成され、前記最小および最大しきい値が、生成されるバイナリマスクを段階的に変化させて、視覚的に認識できる葉片の境界と比較することによって決定されることを特徴とする請求項１〜３のいずれか一項に記載の方法。
【請求項５】
エッジ補正時に、葉片の縁部を認識するために、葉片の重心Ｃｉからの距離ｒの増加に伴う画素数または画素面積の増加
【数１】
0005396484.tif 000010
が使用されることを特徴とする請求項１〜４のいずれか一項に記載の方法。
【請求項６】
ある画像での各画像点に関するステレオアルゴリズムによる視差画像Ｄの生成時に、別の画像内での対応する点のずれが求められることを特徴とする請求項１〜５のいずれか一項に記載の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、葉片の成長を測定するための方法およびそれに適した装置に関する。
【背景技術】
【０００２】
植物の成長は、重要な栽培特性の１つであり、代謝過程および成長に対する作用物質の影響の研究において非常に興味深いものである。
【０００３】
植物をそのまま（ｇａｎｚ）使用する場合には、広い面積が必要となる。さらに、葉をそのまま測定するのに適しているのは、特定の成長形態（とりわけロゼット植物；非特許文献１）のみである。植物が年を取るにつれて日陰となる部分が増加するため、そのような方法は、非常に若い段階の植物に使用が限定される。
【０００４】
従来技術によれば、植物そのままではなく、打ち抜かれた葉片も使用されている。打ち抜かれた葉片は、適切な栄養液中でさらに数日間成長する（非特許文献２）。それにより、場所を取らずに、多くの植物の成長の測定を高いスループットで行うことができる。さらに、葉片は、作用物質を容易に使用できるようにする。
【０００５】
葉の成長を光学的に測定するための従来の方法は、通常、植物をそのまま使用する。個々の植物の表面の光学的決定は、２次元で、すなわち投影された葉片を測定することによって（非特許文献３；非特許文献１）、または３次元で、すなわちレーザスキャンまたはステレオカメラによって得られる表面モデルによって（非特許文献４）行われる。さらに、固定台の上に載せた葉片の面積を測定するための市販のシステムが存在する（例えば、ＬｅｍｎａＴｅｃ  ＧｍｂＨ）。そのようなシステムは、例えば、草食動物による食い荒らし被害、または細菌病によって引き起こされる被害の定量化に使用される。
【０００６】
固定台は成長を大きく妨げるので、成長実験のためには、葉片が溶液中で自由に浮動しなければならない。実験が進むうちに流体水位が低下すると、投影面積が減少し、実際の面積増加分がそこに重なる。したがって、２Ｄでの従来の光学的測定法は、実際の成長率を求めることができない。高い成長率の葉片の場合、流体水位の低下は、実際の成長率の過小評価をもたらすだけであるが、（水位低下によって見掛け上は負の成長が測定されるようなオーダーの値を有する）低い成長率の葉片では、有意な測定を全く行うことができない。多くの植物が、１日の間に顕著な成長過程を示すので、実際の成長率の測定が非常に重要である。
【０００７】
葉片が斜めに置かれているとき、または湾曲しているとき（これは例えば植物ホルモンの作用下で生じる）、投影面積がより小さくなる。したがって、非ステレオ法を用いた成長の測定も可能でない。
【先行技術文献】
【非特許文献】
【０００８】
【非特許文献１】Ｗａｌｔｅｒ  Ａ．，  Ｓｃｈａｒｒ  Ｈ．，  Ｇｉｌｍｅｒ  Ｆ．，  Ｚｉｅｒｅｒ  Ｒ．，  Ｎａｇｅｌ  Ｋ．  Ａ．，  Ｅｒｎｓｔ  Ｍ．，  Ｗｉｅｓｅ  Ａ．，  Ｖｉｒｎｉｃｈ  Ｏ．，  Ｃｈｒｉｓｔ  Ｍ．  Ｍ．，  Ｕｈｌｉｇ  Ｂ．，  Ｊｕｅｎｇｅｒ  Ｓ．  ＆  Ｓｃｈｕｒｒ  Ｕ．  （２００７）  Ｄｙｎａｍｉｃｓ  ｏｆ  ｓｅｅｄｌｉｎｇ  ｇｒｏｗｔｈ  ａｃｃｌｉｍａｔｉｏｎ  ｔｏｗａｒｄｓ  ａｌｔｅｒｅｄ  ｌｉｇｈｔ  ｃｏｎｄｉｔｉｏｎｓ  ｃａｎ  ｂｅ  ｑｕａｎｔｉｆｉｅｄ  ｖｉａ  ＧＲＯＷＳＣＲＥＥＮ：  ａ  ｓｅｔｕｐ  ａｎｄ  ｐｒｏｃｅｄｕｒｅ  ｄｅｓｉｇｎｅｄ  ｆｏｒ  ｒａｐｉｄ  ｏｐｔｉｃａｌ  ｐｈｅｎｏｔｙｐｉｎｇ  ｏｆ  ｄｉｆｆｅｒｅｎｔ  ｐｌａｎｔ  ｓｐｅｃｉｅｓ．  Ｎｅｗ  Ｐｈｙｔｏｌ．  １７４，  ４４７−４５５．
【非特許文献２】Ｗａｌｔｅｒ  Ａ．  （２０００）  Ｒａｅｕｍｌｉｃｈｅ  ｕｎｄ  ｚｅｉｔｌｉｃｈｅ  Ｗａｃｈｓｔｕｍｓｍｕｓｔｅｒ  ｉｎ  Ｗｕｒｚｅｌｎ  ｕｎｄ  Ｂｌａｅｔｔｅｒｎ  ｄｉｋｏｔｙｌｅｒ  Ｐｆｌａｎｚｅｎ．  Ｄｉｓｓｅｒｔａｔｉｏｎ，  Ｕｎｉｖ．  Ｈｅｉｄｅｌｂｅｒｇ．
【非特許文献３】Ｇｒａｎｉｅｒ  Ｃ，  Ａｇｕｉｒｒｅｚａｂａｌ  Ｌ．，  Ｃｈｅｎｕ  Ｋ．，  Ｃｏｏｋｓｏｎ  Ｓ．  Ｊ．，  Ｄａｕｚａｔ  Ｍ．，  Ｈａｍａｒｄ  Ｐ．，  Ｔｈｉｏｕｘ  Ｊ．，  Ｒｏｌｌａｎｄ  Ｇ．，  Ｂｏｕｃｈｉｅｒ−Ｃｏｍｂａｕｄ  Ｓ．，  Ｌｅｂａｕｄｙ  Ａ．，  Ｍｕｌｌｅｒ  Ｂ．，  Ｓｉｍｏｎｎｅａｕ  Ｔ．  ＆  Ｔａｒｄｉｅｕ  Ｆ．  （２００６）  ＰＨＥＮＯＰＳＩＳ，  ａｎ  ａｕｔｏｍａｔｅｄ  ｐｌａｔｆｏｒｍ  ｆｏｒ  ｒｅｐｒｏｄｕｃｉｂｌｅ  ｐｈｅｎｏｔｙｐｉｎｇ  ｏｆ  ｐｌａｎｔ  ｒｅｓｐｏｎｓｅｓ  ｔｏ  ｓｏｉｌ  ｗａｔｅｒ  ｄｅｆｉｃｉｔ  ｉｎ  Ａｒａｂｉｄｏｐｓｉｓ  ｔｈａｌｉａｎａ  ｐｅｒｍｉｔｔｅｄ  ｔｈｅ  ｉｄｅｎｔｉｆｉｃａｔｉｏｎ  ｏｆ  ａｎ  ａｃｃｅｓｓｉｏｎ  ｗｉｔｈ  ｌｏｗ  ｓｅｎｓｉｔｉｖｉｔｙ  ｔｏ  ｓｏｉｌ  ｗａｔｅｒ  ｄｅｆｉｃｉｔ．  Ｎｅｗ  Ｐｈｙｔｏｌ．  １６９，  ６２３−６３５．
【非特許文献４】Ｋａｍｉｎｕｍａ  Ｅ．，  Ｈｅｉｄａ  Ｎ．，  Ｔｓｕｍｏｔｏ  Ｙ．，  Ｙａｍａｍｏｔｏ  Ｎ．，  Ｇｏｔｏ  Ｎ．，  Ｏｋａｍｏｔｏ  Ｎ．，  Ｋｏｎａｇａｙａ  Ａ．，  Ｍａｔｓｕｉ  Ｍ．，  ａｎｄ  Ｔｏｙｏｄａ  Ｔ．  （２００４）  Ａｕｔｏｍａｔｉｃ  ｑｕａｎｔｉｆｉｃａｔｉｏｎ  ｏｆ  ｍｏｒｐｈｏｌｏｇｉｃａｌ  ｔｒａｉｔｓ  ｖｉａ  ｔｈｒｅｅ−ｄｉｍｅｎｓｉｏｎａｌ  ｍｅａｓｕｒｅｍｅｎｔ  ｏｆ  Ａｒａｂｉｄｏｐｓｉｓ．  Ｐｌａｎｔ  Ｊ．  ３８，  ３５８−６５．
【非特許文献５】Ｂｏｕｇｕｅｔ  ２００５
【非特許文献６】Ｆａｕｇｅｒａｓ  ｅｔ  ａｌ．  １９９６
【非特許文献７】Ｓｃｈａｒｒ  ｅｔ  ａｌ．，  ２００７
【非特許文献８】Ｓｃｈａｒｒ  Ｈ．  （２００４）．  Ｏｐｔｉｍａｌ  ｆｉｌｔｅｒｓ  ｆｏｒ  ｅｘｔｅｎｄｅｄ  ｏｐｔｉｃａｌ  ｆｌｏｗ．  Ｉｎ  Ｉｎｔｅｒｎａｔｉｏｎａｌ  Ｗｏｒｋｓｈｏｐ  ｏｎ  Ｃｏｍｐｌｅｘ  Ｍｏｔｉｏｎ  （ＩＷＣＭ）  １４−２９．
【非特許文献９】Ｊａｅｈｎｅ  Ｂ．  （２００２）  Ｄｉｇｉｔａｌ  Ｉｍａｇｅ  Ｐｒｏｃｅｓｓｉｎｇ  （３ｒｄ  ｅｄ．）：  Ｃｏｎｃｅｐｔｓ，  ａｌｇｏｒｉｔｈｍｓ，  ａｎｄ  ｓｃｉｅｎｔｉｆｉｃ  ａｐｐｌｉｃａｔｉｏｎｓ．  Ｓｐｒｉｎｇｅｒ−Ｖｅｒｌａｇ，  Ｌｏｎｄｏｎ，  ＵＫ
【非特許文献１０】Ｂｉｓｋｕｐ  Ｂ．，  Ｓｃｈａｒｒ  Ｈ．，  Ｓｃｈｕｒｒ  Ｕ．  ＆  Ｒａｓｃｈｅｒ  Ｕ．  （２００７）  Ａ  Ｓｔｅｒｅｏ  ｉｍａｇｉｎｇ  Ｓｙｓｔｅｍ  ｆｏｒ  ｍｅａｓｕｒｉｎｇ  ｓｔｒｕｃｔｕｒａｌ  ｐａｒａｍｅｔｅｒｓ  ｏｆ  ｐｌａｎｔ  ｃａｎｏｐｉｅｓ．  Ｐｌａｎｔ，  Ｃｅｌｌ  ａｎｄ  Ｅｎｖｉｒｏｎ．  ３０，  １２９９−１３０８．
【非特許文献１１】Ｋｎｕｔｓｏｎ  ＆  Ｗｅｓｔｉｎ，  １９９３
【非特許文献１２】Ｊｅｎｎｅｓｓ  ＪＳ  （２００４）  Ｃａｌｃｕｌａｔｉｎｇ  ｌａｎｄｓｃａｐｅ  ｓｕｒｆａｃｅ  ａｒｅａ  ｆｒｏｍ  ｄｉｇｉｔａｌ  ｅｌｅｖａｔｉｏｎ  ｍｏｄｅｉｓ．  Ｗｉｌｄｌｉｆｅ  Ｓｏｃｉｅｔｙ  Ｂｕｌｌｅｔｉｎ  ３２，  ８２９−８３９
【非特許文献１３】Ｒｕｓｓ  Ｊ．  Ｃ．  （２００２）  Ｔｈｅ  ｉｍａｇｅ  ｐｒｏｃｅｓｓｉｎｇ  ｈａｎｄｂｏｏｋ．  ＣＲＣ  Ｐｒｅｓｓ，  Ｉｎｃ．，  Ｂｏｃａ  Ｒａｔｏｎ，  ＦＬ
【非特許文献１４】Ｅｌｌｉｓ  ＣＡ．  （１９９９）  Ｗｏｒｋｆｌｏｗ  ｔｅｃｈｎｏｌｏｇｙ．  Ｉｎ：  Ｂｅａｕｄｏｕｉｎ−Ｌａｆｏｎ，  Ｍ．  （Ｅｄ．），  Ｃｏｍｐｕｔｅｒ  Ｓｕｐｐｏｒｔｅｄ  Ｃｏ−ｏｐｅｒａｔｉｖｅ  Ｗｏｒｋ  Ｔｒｅｎｄｓ  ｉｎ  Ｓｏｆｔｗａｒｅ  ７，  ２９−５４．  Ｊｏｈｎ  Ｗｉｌｅｙ  ＆  Ｓｏｎｓ  Ｌｔｄ．，  Ｃｈｉｃｈｅｓｔｅｒ，  ＵＫ
【非特許文献１５】Ｚｈａｎｇ  Ｚ．  （１９９９）  Ｆｌｅｘｉｂｌｅ  ｃａｍｅｒａ  ｃａｌｉｂｒａｔｉｏｎ  ｂｙ  ｖｉｅｗｉｎｇ  ａ  ｐｌａｎｅ  ｆｒｏｍ  ｕｎｋｎｏｗｎ  ｏｒｉｅｎｔａｔｉｏｎｓ．  Ｉｎ：  Ｉｎｔ．  Ｊ．  Ｃｏｍｐｕｔ．  Ｖｉｓｉｏｎ  ６６６−６７３．
【非特許文献１６】Ｚｈａｎｇ  Ｚ．  （２０００）  Ａ  ｆｌｅｘｉｂｌｅ  ｎｅｗ  ｔｅｃｈｎｉｑｕｅ  ｆｏｒ  ｃａｍｅｒａ  ｃａｌｉｂｒａｔｉｏｎ．  ＩＥＥＥ  Ｔ．  Ｐａｔｔｅｒｎ  Ａｎａｌ．  ２２：  １３３０−１３３４．
【非特許文献１７】Ｄｉｊｋｓｔｒａ  Ｅ．Ｗ．  （１９５９）  Ａ  ｎｏｔｅ  ｏｎ  ｔｗｏ  ｐｒｏｂｌｅｍｓ  ｉｎ  ｃｏｎｎｅｘｉｏｎ  ｗｉｔｈ  ｇｒａｐｈｓ．  Ｎｕｍｅｒ．  Ｍａｔｈ．  １：  ２６９−２７１
【非特許文献１８】Ｆａｕｇｅｒａｓ  Ｏ．，  Ｈｏｔｚ  Ｂ．，  Ｍａｔｈｉｅｕ  Ｈ．，  Ｖｉｅｖｉｌｌｅ  Ｔ．，  Ｚｈａｎｇ  Ｚ．，  Ｆｕａ  Ｐ．，  Ｔｈｅｒｏｎ  Ｅ．，  Ｍｏｌｌ  Ｌ．，  Ｂｅｒｒｙ  Ｇ．，  Ｖｕｉｌｌｅｍｉｎ  Ｊ．，  Ｂｅｒｔｉｎ  Ｐ．，  Ｐｒｏｙ  Ｃ．  （１９９３）  Ｒｅａｌ  ｔｉｍｅ  ｃｏｒｒｅｌａｔｉｏｎ−ｂａｓｅｄ  Ｓｔｅｒｅｏ：  ａｌｇｏｒｉｔｈｍ  ｉｍｐｌｅｍｅｎｔａｔｉｏｎｓ  ａｎｄ  ａｐｐｌｉｃａｔｉｏｎｓ．  Ｔｅｃｈｎｉｃａｌ  Ｒｅｐｏｒｔ  ２０１３，  ＩＮＲＩＡ  Ｓｏｐｈｉａ  Ａｎｔｉｐｏｌｉｓ
【非特許文献１９】Ｈａｒｔｌｅｙ  ａｎｄ  Ｚｉｓｓｅｒｍａｎ，  ２００４
【非特許文献２０】Ｔｒｕｃｃｏ  ａｎｄ  Ｖｅｒｒｉ，  １９９８
【非特許文献２１】Ｋｌｅｔｔｅ  Ｒ．，  Ｚａｍｐｅｒｏｎｉ  Ｐ．，  （１９９５）  Ｈａｎｄｂｕｃｈ  ｄｅｒ  Ｏｐｅｒａｔｏｒｅｎ  ｆｕｅｒ  ｄｉｅ  Ｂｉｌｄｖｅｒａｒｂｅｉｔｕｎｇ．  Ｂｒａｕｎｓｃｈｗｅｉｇ−Ｖｉｅｗｅｇ，  ２．Ａｕｆｌ．
【非特許文献２２】Ｗａｌｔｅｒ  Ａ，  Ｓｃｈｕｒｒ  Ｕ  （１９９９）  Ｔｈｅ  ｍｏｄｕｌａｒ  ｃｈａｒａｃｔｅｒ  ｏｆ  ｇｒｏｗｔｈ  ｉｎ  Ｎｉｃｏｔｉａｎａ  ｔａｂａｃｕｍ  ｐｌａｎｔｓ  ｕｎｄｅｒ  ｓｔｅａｄｙ−ｓｔａｔｅ  ｎｕｔｒｉｔｉｏｎ．  Ｊ．  Ｅｘｐ．  Ｂｏｔ．  ５０：  １  １６９−１１７７
【非特許文献２３】Ｖａｎ  Ｒｏｓｓｕｍ，  ２００６
【発明の概要】
【発明が解決しようとする課題】
【０００９】
したがって、本発明の課題は、成長（すなわち葉片の面積成長）を正確に測定し、１日の間の成長のダイナミクスを解明することができるようにする方法を提供することである。さらに、本発明の目的は、多様な植物種の標準化された成長測定を行うことができるようにすることである。
【課題を解決するための手段】
【００１０】
この課題は、本発明によれば、請求項１の前提部を基に、請求項１の特徴部に記載の特徴によって解決される。さらに、この課題は、本発明によれば、請求項７の前提部を基に、請求項７の特徴部に記載の特徴によって解決される。
【００１１】
本発明による方法および装置を用いて、高いスループットで、例えば約５００枚の葉片という標本量で成長測定を行うことができ、その際、成長および約１〜１．５時間の間隔での１日の間の成長のダイナミクスを把握または解明することができることが示された。本発明による方法および装置は、時間に関して高い分解能での葉片成長の分析を可能にする。本方法および装置は、調節および較正後、さらに手を加える必要はなく、人手をかけずに自動で２〜５日から約２週間にわたって動作させることができる。
【００１２】
本発明による方法および装置により、葉片に基づいて、植物の（光学的検出に関する）実際の成長率を実現することができる。これは、これまで、従来技術で知られている光学的測定法では可能でなかった。
【００１３】
本発明の方法は、葉片の成長を測定するための方法であって、以下の方法ステップを特徴とする方法に関する。
【００１４】
１．撮影システムの較正：
カメラの較正は、例えば非特許文献５に従って標準的な方法で行うことができる。しかし、従来技術により知られている別の方法も可能である。カメラ位置の相対位置が事前に与えられているので、別途のステレオ較正は必要ない。後の３Ｄ再構成のために、較正パラメータが分かっている必要がある。このステップで決定されるパラメータは、焦点距離、中心点（＝光軸と像面の交点）、半径方向および接線方向歪みである。
【００１５】
２．撮像：
選択するステレオ法に応じて、１枚の葉片につき、ａ）相関ベースの方法（非特許文献６）では２枚の画像を撮影することができ、ｂ）微分フィルタによってオプティカルフローを求める方法（非特許文献７）では３つ以上の画像を撮影することができる。本方法において画像数をより多くすると、系統的誤差を大幅に減少させることができるようになる（非特許文献８）。上記のカメラを第３のスライドテーブルによって摺動させることによって、十分に小さいステレオ基線が実現される。第３のスライドテーブルの代わりに、２台のカメラを使用することもできる。カメラとしては、例えば産業用カメラまたはステレオカメラを使用することができる。
【００１６】
３．画像データの処理：
画像データを解析するために、栄養液中で自由に浮動する葉片と背景とを区別しなければならず、または葉片の画像データをセグメント化しなければならない。赤外線フィルタを使用するので、セグメント化のためにグレースケール画像しか使用することができない。
【００１７】
まず、ａ）しきい値セグメント化（例えば、非特許文献９）を行う。従来技術で知られているしきい値セグメント化のための別の方法も適している。この場合、画素が順次に選択されて、背景か葉片かに従って選別される。
【００１８】
その後、はじめに、ｂ）複数回の形態学的侵食ステップによって縁部領域（例えば容器の壁、流体表面での反射）を除去し、次いで、残っている物体を（画素単位で）最小限の大きさで選択することによって、葉片と考えられる部分を識別する。
【００１９】
次いで、ｃ）葉片が入っている容器の縁部を除外するエッジ補正を行う。（図９参照）。
【００２０】
４．３Ｄ再構成：
ステレオアルゴリズム（例えば、非特許文献１０）によって、視差マップを作成する。
【００２１】
ステレオカメラ、ステレオビジョン、深さ画像解析に関連して視差があり、あるいはまた２つの異なる像面上での投影図中で同一の物体が占める位置のオフセットのずれがある。その際、像面に属する焦点は、基線長（Ｂａｓｉｓ）ｂによって空間的に互いに分離される。次いで、２つのレンズの焦点距離ｆを用いて、以下の関係が得られる：ｒ＝ｂｆ／Ｄ。ここで、Ｄは視差を表す。したがって、ステレオ画像における視差の測定により、物体に対する距離ｒを算出することができる。
【００２２】
視差画像の後処理
外れ値、すなわち周辺の点から大きくかけ離れた個々の３Ｄ点を、メジアンフィルタ（例えば、非特許文献９）によって除去する。例えばコントラストが低すぎるために視差評価が可能でない領域は、正規化された畳み込み（非特許文献１１）によって閉じる。次いで、事前に決定された較正パラメータに基づいて、視差マップから３Ｄ表面モデルを計算する（例えば、非特許文献１２）。
【００２３】
５．表面の平滑化および測定：
雑音による面積測定の過大評価を妨げるために、前に得られた表面モデルを、Ｚ画像のガウスフィルタ処理によって平滑化する（例えば、非特許文献９）。
【００２４】
６．成長率：
相対成長率を、面積値の時系列から求める。測定した成長率は、さらなる解析のためにデータバンク内に用意しておく。
【００２５】
本発明の核心は、正確な面積測定を可能にする葉片の３Ｄ表面モデルが得られるステレオ撮像の使用である。
【００２６】
葉片を、例えば水平な面上でマイクロタイタープレートまたはペトリ皿に載せる。少なくとも１台のカメラと、ＬＥＤアレイと、スライドテーブルとからなる撮影システムによって、葉片の光学的検出を行う。Ｘ−Ｙ−スライドテーブルによって、カメラシステムをそれぞれの葉片の上に位置決めすることができる。カメラシステムは、例えば、撮影対象の植物を均質に照明する複数の赤外ＬＥＤ照明群からなる。照明の位置は静止したままであり、一方、様々なカメラの位置は３基のスライドテーブルによって移動される。植物への生理学的な作用を避けるために、撮像用の照明が赤外範囲内で行われることが好ましい。それに対応して、外部照明に依存せずにできるだけ一定の照明強度を保証するために、カメラに赤外線フィルタが装備される。光合成で利用できる波長を放出する蛍光灯または別の光源によって、規定の明暗周期で光合成活性照明を行うことができる。規定の温度および湿度を保証するために、この方法は人工気候室内で実施することができる。
【００２７】
本発明による方法は、例えば約５００枚の葉片の同時測定を可能にする。様々な作用物質を用いた実験から、除草剤作用の下で成長率がより低い場合に、本発明による方法によって正確な（すなわち投影効果によって歪曲されていない）成長率が求められることが示されている。それに対して、投影面積を使用する場合には、除草剤作用の下では、見掛け上は負の成長率が生じた。
【００２８】
例えば、ベルトコンベアによって、周期的に測定されるマイクロタイタープレートの数を増加させることができる。
【００２９】
例えば、マイクロタイタープレート用の保管装置およびロボットアームによって、周期的に測定されるマイクロタイタープレートの数を増加させることができる。
【００３０】
本発明による装置は、大型で高価な人工気候室ではなく、例えば空調管理された筐体内に置くことができる。それにより、本発明による方法および装置は、例えば空調管理されていない実験室で、ほぼ自立式に使用可能である。
【００３１】
さらに、本発明は、上述した方法に適した装置に関する。
【００３２】
この装置は、赤外線フィルタを有する少なくとも１台のカメラと、赤外線照明と、カメラおよび赤外線照明用のＸ／Ｙ平面内でのスライドテーブルと、葉片用の保管容器と、電子解析および制御ユニットとを備える。
【００３３】
有利な発展形態は、従属請求項に記載される。
【００３４】
図面は、本発明による方法および装置の例示的実施形態、ならびに成長測定の結果を示す。
【図面の簡単な説明】
【００３５】
【図１】撮影システムを備える概略構造を示す図である。Ｘ、Ｙ、およびＸ２＝スライドテーブル；Ｘスライドテーブルは構造に取り付けられている。矢印＝スライドテーブルの移動方向；Ｌ＝赤外線ランプ；Ｃ＝カメラ。
【図２】葉片が入っているマイクロタイタープレートの上に配置された、４つのＬＥＤアレイによって取り囲まれたカメラを示す図である。
【図３】赤外線照明で撮影された３９５枚の個別画像からなる概要画像を示す図である。
【図４】溶液上で浮動するタバコ（Ｎｉｃｏｔｉａｎａ  Ｔａｂａｃｕｍ）の葉片を示す図である。表面がわずかに湾曲している。
【図５】（子葉を含めて、植物の根元から数えた）葉番号を付した葉と、中央脈に沿った穿孔位置を示す図である（Ｂ：付け根；Ｍ：中央；Ｔ：先端）。
【図６】複数の個別画像からまとめられた概要画像の一部と、そこに重ねて記された個々の葉片に関する情報（植物識別番号、処理、葉番号、穿孔位置）を示す図である。
【図７】様々な流体液位での一定の大きさの葉片の面積測定を示す図である。ｘ：流体液位（ｍｍ）、ｙ：（水位０ｍｍでの面積に対する）相対面積。水位は、０ｍｍ（マイクロタイタープレートの凹部の底部）〜１０ｍｍの間で変化する。Ａｐ（白四角）およびＡ３Ｄ（黒丸）に関する平均値±標準誤差；ｎ＝２４
【図８】葉片を入れたマイクロタイタープレートを示す図である。左半分：栄養液を用いた処理。右半分：除草剤グリホサートを用いた処理（実験の開始から７２時間後に画像を撮影）。
【図９】グレースケール画像のしきい値セグメント化によって生成されるバイナリマスクＭｔを示す図である。（白：前景、黒：背景；灰色：葉片。
【数１】
0005396484.tif 000002
の計算によって認識される）。最終的なセグメント化マスクＭｓは、認識された葉片の重心を中心としてｒの円内にあり、Ｍｔにおいて値１を有していたＭｔからの画素を含む。
【図１０】重心からの距離の増加に伴う面積成長を示す図である。ｘ：
【数２】
0005396484.tif 000003
（画素）、ｙ：ｒ（画素）。×印：測定された増大；点線：ｒの増加に伴う黒丸の理論上の増大２πｒ。
【図１１】背景（植物でない部分）が黒く表された、後処理された視差画像を示す図である。スカラ：視差範囲［６１０：６２０］画素。
【図１２】再構成された葉片の３Ｄ図である。
【図１３】本発明による方法の分散図である。ギガビットネットワークが、ＭＡＳＴＥＲ（撮像、データバンクシステム、およびワークフローエンジンを備える）、ＣＬＩＥＮＴコンピュータ、データサーバ、バックアップシステム、およびＬｉｎｕｘクラスターを接続する。長方形のボックス：スクリーニングプログラム。
【図１４】ＵＭＬ（統一モデリング言語（Ｕｎｉｆｉｅｄ  Ｍｏｄｅｌｉｎｇ  Ｌａｎｇｕａｇｅ））標記での、主要なワークフローステップのアクティビティ図である。ワークフロー規則は、データバンクテーブルに規則を記憶されている。２種類の異なるワークフロー規則がある。（１）ジョブのステップと状態の特定の組合せに関して発動される状態ベースの規則と、（２）新規のジョブを生成しなければならないかどうか決定するためにプログラムコードが実行されるコードベースの規則（例えば、撮像工程の終了後の概要画像の生成）である。
【図１５】スケーラビリティ、すなわち入力量が増加したときのリソース要件、ここでは特に時間要件に関するプログラムの挙動を示す図である。ｘ：使用されるＣＰＵの数；ｙ：％実行速度（順次処理に対して正規化）。白四角：ＥＭＰＴＹジョブ；白丸：ＣＰＵ＿ＢＯＵＮＤジョブ；黒四角：ＩＯ＿ＢＯＵＮＤジョブ；黒丸：ＭＩＸＥＤジョブ。誤差バーが、標準誤差（ｎ＝３）を示す。
【発明を実施するための形態】
【００３６】
以下、本発明による方法および本発明による装置を詳細に説明する例示的実施形態を提示する。
【００３７】
・例示的実施形態
・本発明による方法の精度
投影される葉面積Ａｐ（投影葉面積）を求める精度のテストとして、タバコ（Ｎｉｃｏｔｉａｎａ  ｔａｂａｃｕｍ）の葉片２４枚に穿孔し、すぐにフラットベッドスキャナ（Ｅｐｓｏｎ  Ｐｅｒｆｅｃｔｉｏｎ  ２４８０  Ｐｈｏｔｏ，Ｅｐｓｏｎ（日本国長野県諏訪市））を用いて、４００ｄｐｉの解像度で白い背景の前でスキャンした。スキャナによって生成されたＲＧＢ画像を、ＨＳＶ（色相、彩度、明度値）色空間に変換した。次いで、チャネルＨ、Ｓ、Ｖでのしきい値（最小／最大）を使用することによって、葉標本を葉と背景にセグメント化した。この方法は、フルカラー情報を利用できるので、葉緑素の識別に最適である（非特許文献１３；非特許文献１；非特許文献１０）。このようにして求められたＡｐのばらつき係数（ＣＶ）は３．８％であり、すなわち、コルク穿孔器によって穿孔された葉片の面積は３．８％のばらつきがあった。それに対し、本発明による方法によって測定された２４枚の葉片の面積ＡｐのＣＶは３．２％であり、ＣＶ  Ａ３Ｄ（本発明による３次元測定法によって再構成された葉面積の面積値）は３．８％であった。一般に、ＣＶ（Ａ３Ｄ）はＣＶ（ＡＰ）よりも大きい値を取るはずである。なぜなら、コルク穿孔器が、葉片をほぼ完全な円形の周縁で抜き取るが、表面は、隆起していたり波打っていたりする可能性があるからである。この精度決定から、本発明による方法により、色情報がないにも関わらず、従来の画像処理法よりも良くはないにせよ同等に良好なセグメント化結果を得られることが分かった。
【００３８】
面積決定の精度に対するカメラと対象物の離隔距離の影響を調べるために、様々な流体液位で浮動する葉片の画像を撮影した（図７）。投影された葉面積Ａｐは、明らかに離隔距離に依存していた。すなわち、カメラ距離を１０ｍｍ縮めると、葉面積は１４％増大した。逆に、Ａ３Ｄはほぼ一定であり、最大のばらつきが２％であった。これは、測定期間が非常に長く、蒸発により流体水位が著しく低下するときに正確な成長測定を行うために、ここで提示する本発明による３Ｄ再構成法が必要であることを明らかに示す。
【００３９】
成長分析
成長測定の事例研究のために、計４５８枚のタバコ（Ｎｉｃｏｔｉａｎａ  ｔａｂａｃｕｍ）の葉片を使用した。可能な最速の撮像では、各葉片を１．５時間で一回撮影した（１枚の葉片につき１対の画像）。これは、カメラの位置決めを含めて、１枚の葉片につき１１．８秒の平均撮影時間に相当する。７２時間後にこの実験を終了した。合計で約２２０００枚の画像を撮影した。
【００４０】
画像処理システム
ハードウェア構造（図１および図２）は、２ＤスクリーニングシステムＧＲＯＷＳＣＲＥＥＮ（非特許文献１）の設計に基づく。本発明によるシステムは、５８ｃｍ×５８ｃｍの矩形の測定領域を検出する。Ｐｏｉｎｔ  Ｇｒｅｙ  Ｓｃｏｒｐｉｏｎ  ＳＣＯＲ−２０ＳＯＭカメラ（Ｐｏｉｎｔ  Ｇｒｅｙ  Ｒｅｓｅａｒｃｈ，Ｖａｎｃｏｕｖｅｒ，ＢＣ，カナダ）を用いて画像を撮影し（２メガ画素；最大解像度１６００×１２００画素）、ＦｉｒｅＷｉｒｅ４００ケーブルによって制御計算機に伝送する（ＩＥＥＥ  １３９４ａ；最大転送率：４００Ｍビット）。赤外線（ＩＲ）フィルタ（Ｍｏｄｅｌ  ２７  ０９３；Ｓｃｈｎｅｉｄｅｒ−Ｋｒｅｕｚｎａｃｈ，Ｂａｄ−Ｋｒｅｕｚｎａｃｈ，Ｇｅｒｍａｎｙ;８００ｎｍよりも長い波長で透過性）を有するＣ−Ｍｏｕｎｔ−Ｏｂｊｅｋｔｉｖ（１５ｍｍ；Ｒｏｄｅｎｓｔｏｃｋ，Ｍｕｅｎｃｈｅｎ,Ｇｅｒｍａｎｙ）によって画像を撮影する。カメラは、垂直下方向に向ける。撮像用の照明として、それぞれ１８個のＬＥＤ（λｍａｘ＝８００ｎｍ；Ｃｏｎｒａｄ  Ｅｌｅｃｔｏｒｎｉｃ，Ｈｉｒｓｃｈａｕ，Ｇｅｒｍａｎｙ）を有する４つのＩＲ−ＬＥＤアレイを使用する。ＩＲ照明は、生理学的に活性でなく、かつより高いコントラストで植物標本が見えるようにするという追加の利点を有する。
【００４１】
照明強度の温度依存ドリフトを防ぐために、実験全体を通じてＬＥＤ照明を点けておく。使用したＬＥＤアレイは、２０分超経過してからようやく、安定強度の１％に近づいた。画像雑音を減少させるために、１つの撮影位置につき５枚の画像を撮影し、その平均値を取る。画像は、１２８０×９６０画素の解像度で撮影し、ＴＩＦＦ形式で、損失のないＰａｃｋＢｉｔｓ圧縮で記憶した。
【００４２】
２基の高精度スライドテーブル（「Ｘ」および「Ｙ」；Ｐｉｃｏ−Ｍａｘｉ，Ｔｙｐｅ  ＦＭＤ−ＬＰＴ８０．５５０．１２０５−ＳＭ，Ｌａｓｅｒ  ２０００  ＧｍｂＨ，Ｍｕｅｎｃｈｅｎ，Ｇｅｒｍａｎｙ；マイクロステップモータ：ＭＤｒｉｖｅ  ２３  Ｐｌｕｓ；Ｉｎｔｅｌｌｉｇｅｎｔ  Ｍｏｔｉｏｎ  Ｓｙｓｔｅｍｓ，Ｍａｒｌｂｏｒｏｕｇｈ，ＵＳＡ）によって、カメラとＩＲ照明ユニットを水平面内で移動させる。第３のスライドテーブル（「Ｘ２」）が、１対または１列のステレオ画像の撮影時にカメラをＸ方向に摺動させるように働き、一方、ＩＲ照明ユニットの位置は変化しない。これは、照明が一定であることがステレオ撮影の前提条件であるので必須である。撮影および位置決めシステム全体が、Ｘ−９５プロファイル要素（Ｌｉｎｏｓ  Ｐｈｏｔｏｎｉｃｓ，Ｇｏｅｔｔｉｎｇｅｎ，Ｇｅｒｍａｎｙ）からなる安定な骨組みに取り付けられる。
【００４３】
調べる葉片は、２４個の凹部を有するマイクロタイタープレートに保管する（図８；Ｎｕｎｃ，Ｒｏｓｋｉｌｄｅ，デンマーク）。凹部は、均等に配列され、それぞれ体積が２．５ｍｌである。実験期間にわたって、プレートの蓋を外して、妨げのないガス交換を可能にし、曇りを防止する。リフトテーブル（Ｓｗｉｓｓ  ＢＯＹ  １１０，Ｒｕｄｏｌｆ  Ｇｒａｕｅｒ  ＡＧ，Ｄｅｇｅｒｓｈｅｉｍ，スイス）の上に検査対象物を載せ、このリフトテーブルによって撮像に適した距離を調節することができる。水準器によって、位置決めシステムも検査対象物も水平向きであることが保証される。コントラストを向上させる（ひいては、それによりセグメント化を容易にする）ために、台を黒色のプラスチックで覆う。
【００４４】
コンピュータシステム
図１３は、本発明による方法のハードウェアコンポーネントと個別の計算機システム上で実行されるソフトウェアコンポーネントとを示す概略図である。撮像は、パーソナルコンピュータ（ＰＣ）上で行われる（「ＭＡＳＴＥＲ」と表されている；Ｉｎｔｅｌ（登録商標）Ｐｅｎｔｉｕｍ（登録商標）４  ＤｕａｌＣｏｒｅ，２  ＧＢ  ＤＤＲ４００  ＲＡＭ（Ｉｎｔｅｌ，Ｓａｎｔａ  Ｃｌａｒａ，ＵＳＡ），ＳｕＳＥ  Ｌｉｎｕｘ  １０．２；ＳｕＳＥ，２００６）。このＰＣ上で、マルチユーザデータバンクが実行される（下記参照）。スクリーニングシステムは、自動画像解析のために、複数のさらなるＰＣ（クラスターノード；１  ＧＢ  ＤＤＲ  ４００  ＲＡＭ，ＳｕＳＥ  Ｌｉｎｕｘ  １０．２  （Ｎｏｖｅｌｌ，Ｐｒｏｖｏ，ＵＳＡ））を利用できる。大量のデータを処理するために、使用するすべての計算機がギガビットネットワークを介して接続される。
【００４５】
データ管理
得られたデータは、ＭｙＳＱＬデータバンク（ＭｙＳＱＬ  ＡＢ，Ｕｐｐｓａｌａ，Ｓｗｅｄｅｎ）およびデータサーバ上で管理する。スクリーニングシステムの基本データおよび動的データをデータバンクに記憶し、ネットワークを介して呼び出すことができる。データバンクは、コアアプリケーションとクラスターノードの間の通信の主インターフェースを構成する。
【００４６】
画像データおよび他の大容量データがデータサーバ上に保存される。このサーバには、ＭＡＳＴＥＲおよびすべてのノードが、ＣＩＦＳ（共通インターネットファイルシステム（Ｃｏｍｍｏｎ  ｉｎｔｅｒｎｅｔ  ｆｉｌｅ  ｓｙｓｔｅｍ））を介してアクセスできる。このサーバ上のデータは、そのパス名に基づいて参照される。データバンクは、大容量データに関して、メタ情報（データタイプ、作成日、作成者）を記憶する。
【００４７】
並列化およびワークフロー概念
撮影された画像の処理を、ワークフローパラダイム（非特許文献１４）に従って実施した。ここで、タスク全体を、一連の別個のワークフローステップに分ける。このようにすると、１．）処理ステップを別のワークフローで再利用することができ、２．）解析を容易に並列化することができ、３．）例えば構成パラメータの変更に従って、個別のステップを簡単に再び行うことができる。図１４は、葉片の３Ｄ面積を計算する個別のステップを示す。
【００４８】
１．）撮影システムの較正
ａ）カメラ較正
実験を行う前に、確立されている方法によってカメラを較正した（非特許文献１５；非特許文献１６）。様々な向きから、少なくとも３０枚の市松模様の画像を撮影した。Ｊ．Ｙ．Ｂｏｕｇｕｅｔ（Ｉｎｔｅｌ，Ｓａｎｔａ  Ｃｌａｒａ，ＣＡ，ＵＳＡ）のカメラ較正ツールボックスのＯｐｅｎＣＶ実装によって、カメラの内因性パラメータ（焦点距離、中心点、半径方向および接線方向歪み）を評価した。
【００４９】
ｂ）ステレオ較正
ステレオシステムのパラメータ（左のカメラシステムに対する右のカメラシステムの回転および並進）を以下のように決定した。ステレオ基線は、使用される高精度スライドテーブルにより既知である。しかし、カメラのＣＣＤセンサは、取付けにより、摺動方向に対してわずかに回転していることがあるので、修正が必要なことがある。回転している場合、これは摺動距離に線形に依存する。これは、ステレオ対応付け試行の質を低下させることがある。なぜなら、使用されるアルゴリズムは、（Ｘ方向での）カメラセンサの列に沿ったカメラ摺動を行うことに基づいており、したがって対応付け試行も１列でしか予想することができないからである。したがって、修正のために、パノラマステッチソフトウェアＨｕｇｉｎを用いて少なくとも４０点の対応付けを求め、その後、回転を評価する。次いで、画像列と摺動方向が再び一致するように右の画像を回転させる。このプロセスは、カメラの取付け後に１回だけ実施すればよい。あるいは、例えば上記のＪ．Ｙ．Ｂｏｕｇｕｅｔ（Ｉｎｔｅｌ，Ｓａｎｔａ  Ｃｌａｒａ，ＣＡ，ＵＳＡ）によるカメラ較正ツールボックスを用いてステレオ較正を行うこともできる。
【００５０】
２．）撮像
位置決め時間を最短にするために、（Ｄｉｊｋｓｔｒａのアルゴリズムに従って；非特許文献１７）最適な順序で、予め定義された撮像位置（図３）を取る。各撮影位置で、基線長（カメラ中心間の距離）ｂ＝２０ｍｍで２枚の画像を撮影する。それによって得られるより大きな視差が、より良い深さ解像度を可能にするので、ｂをできるだけ大きく選択することが望ましい。他方で、ｂは、視野および葉片の大きさによって制限される。基線長ｂ＝２０ｍｍが選択される場合、実験期間にわたって左右の画像内で葉片を見ることができることが保証された。
【００５１】
３．）葉片のセグメント化／認識
撮像時に白色光ではなく赤外線照明を使用するという生物学的な必要性が、セグメント化を難しくする。カラー画像は、簡単に、かつ正確にセグメント化することができる（非特許文献１３；非特許文献１）。それに対し、ＩＲ照明の下では、限定されたスペクトルにより、グレー値画像しか撮影することができない。暗い背景の使用は有利であるが、グレー値だけに基づいてマイクロタイタープレートの凹部壁と葉片を区別することはできない。したがって、以下のセグメント化プロセスを使用した。
【００５２】
・撮影された画像Ｉａｃｑ内の強度のばらつき（不均質な照明による）を、基準背景画像によって較正する（非特許文献９）。
【００５３】
・しきい値セグメント化：最小および最大しきい値を使用することによって、撮影された画像のバイナリマスク（Ｍｔ）を生成する。これらの値は、葉片の照明強度および輝度に応じて、生成されたマスクを段階的に変化させて、視覚的に認識できる葉片の境界と比較することによって求める。
【００５４】
・形態学的侵食：マスクの縁部を形態学的に侵食する（非特許文献９）。それにより、物体輪郭の不整をなくす。次いで、残った画素から、葉片と考えられる部分それぞれの重心Ｃｉを計算する。
【００５５】
・エッジ補正：葉片の縁部を認識するために、ｃｉからの距離ｒの増加に伴う画素数（画素面積）の増加
【数３】
0005396484.tif 000004
を使用する。
【数４】
0005396484.tif 000005
が葉片の理論上の増大２πｒの０．６倍未満に低下するときの半径ｒｍａｘを、セグメント化の縁部基準として使用する（図９および図１０）。このようにして、最終的なセグメント化マスクＭｓが得られる。１対の画像の両画像からの葉片と考えられる部分を、実際の予想される視差に基づいて一致させる。
【００５６】
所定の位置で事前に葉片が検出されなかった場合、一意の識別子、すなわち植物識別番号を有する新たな葉片をデータバンクに登録する。
【００５７】
投影された葉面積Ａｐは、所定の位置での１対の画像のうち第１の画像のＭｓにおける非ゼロ画素の和として計算される。
【００５８】
データ量を減少させるために、さらなる処理ステップに関して、葉片がある画像部分のみを使用する。
【００５９】
４．）ステレオ対応付け試行および３Ｄ再構成
前の処理ステップからの葉片の１対の画像を使用して、２カメラステレオによって深さマップを計算する。
【００６０】
・画像を調整する、すなわち真直ぐにする。その際、ずれ（視差）がｘ方向のみで生じるように画像の向きを真直ぐにするか、または修正する。それにより、カメラセンサのスキャンラインがカメラの移動方向と一致する（「ステレオ較正」参照）。
【００６１】
・相関ベースのステレオアルゴリズム（非特許文献１８；非特許文献１０）が、調整された１対の画像に対して使用される。アルゴリズムは、画像内の各画像点ごとに、別の画像内での対応する点のずれを求めることを試行する（＝対応付け試行）。例えば非特許文献１８からの相関関数ｃ２を使用して、矩形の画像範囲を比較する。使用するステレオ幾何形状および事前に実施された調整に基づいて、水平方向での摺動のみを行うことができるようになる。したがって、対応付け試行は、２Ｄ問題ではなく１Ｄ問題のみを提示する（非特許文献１９；非特許文献２０）。対応付け試行の結果が、視差画像Ｄである。
【００６２】
・後処理：Ｄでの外れ値をメジアンフィルタ（半径：２５画素）によって取り除く（非特許文献２１）。外れ値とは、暗に、周辺に比べて非常に小さいか、または非常に大きい画素である。周辺にあるすべての画素をそれらの値に従ってソートし、リストの中央の値を取る。そのようにして順序付けを行う。
【００６３】
処理ステップの最後に面積値を計算すべきであるので、再構成された葉片において、存在することがあり得る穴は閉じなければならない。これは、正規化された畳み込みにより、視差評価が可能でなかった点の補間によって行う（非特許文献１１）。後処理により、平滑化された視差画像Ｄｇｅｇが得られる。
【００６４】
・ステレオシステムの内因性および外因性パラメータに基づいて、Ｄｇｅｇ内のすべての点の３Ｄ位置を三角測量で測定する。このようにして、ユークリッド座標系Ｘ、Ｙ、およびＺが定められた座標格子Ｗが得られる。
【００６５】
５．）表面モデルの平滑化
再較正された３Ｄ点群の深さＺを、円板状のメジアンフィルタ（非特許文献９）（半径：１５画素）によって平滑化し、表面の小さい起伏を取り除く。
【００６６】
・３Ｄモデルから副標本を抽出し（ｕｎｔｅｒａｂｇｅｔａｓｔｅｔ）（Ｘ方向とＹ方向で各７点）、葉片の３Ｄ面積Ａ３Ｄを、表面上の個々の点の既知の隣接関係に基づいて求める。
【数５】
0005396484.tif 000006
ここで、
ｈ＝画像の高さ
ｗ＝画像の幅
ＡΔ＝三角形面積
Ｗ＝座標格子
Ｗｘｙ＝格子位置
三角形面積ＡΔは３つのベクトルによって定義される。
【００６７】
６．）成長率の算出
式
【数６】
0005396484.tif 000007
の指数的な成長率を仮定すると、葉組織の相対成長率（ｒｅｌａｔｉｖｅ  ｇｒｏｗｔｈ  ｒａｔｅｓ；ＲＧＲ）は、
【数７】
0005396484.tif 000008
によって定義され、ここで、Ａ（ｔ２）およびＡ（ｔ１）は、それぞれ２つの連続する時点ｔ２およびｔ１に関する投影されたか、または３Ｄ葉面積であり、ｄ−１＝１日当たりである（例えば、非特許文献２２参照）
【００６８】
ソフトウェアアーキテクチャ
本発明による方法は、主にオープンソースソフトウェアを使用する。クライアントソフトウェア（プログラムＣＬＩＥＮＴ）は、Ｃ＋＋／Ｑｔ（Ｔｒｏｌｌｔｅｃｈ，Ｏｓｌｏ，Ｎｏｒｗｅｇｅｎ）で実装される。ＣＬＩＥＮＴプログラムの主なタスクは、１．）新規の測定キャンペーンをパラメータ化すること（測定時間間隔など）；２．）カメラおよびスライドテーブルを調節すること；３．）作業範囲全体の概要画像を生成し、それを用いて、次いで撮影すべき葉片を選択することができること；４．）実行されているシステムを監視すること（解析やシステム負荷などの状態）。プログラムＣＬＩＥＮＴは、ＬｉｎｕｘおよびＷｉｎ３２オペレーティングシステム上で実行される。
【００６９】
性能特性
分散システムの性能特性を測定するために、様々な人工的な解析ジョブの平均実行時間を分析した（表１）。様々な数のクラスターノードで測定を繰り返して、スケーラビリティを求めた。性能テスト用の分散システムは、ＭｙＳＱＬデータバンクおよび撮影ソフトウェアを備える６４ビット２．６ＧＨｚ  ＤｕａｌＣｏｒｅコンピュータと、解析用の５台の３２ビットＩｎｔｅｌ（登録商標）Ｐｅｎｔｉｕｍ（登録商標）４  ＤｕａｌＣｏｒｅコンピュータ（クラスターノード）とからなる。最大構成では、各クラスターノードで、並列して２つのコンシューマプロセスを行った。
【００７０】
【表１】
0005396484.tif 000009
【００７１】
図１５は、様々な人工的なジョブを実行する際のスケーラビリティを示す。ＥＭＰＴＹとＣＰＵ＿ＢＯＵＮＤジョブは、類似した特性を有する。それらは、１〜１０台の範囲の使用されるＣＰＵにわたってほぼ線形にスケールする。ＩＯ＿ＢＯＵＮＤジョブは、それらに比べてスケーリングが良好でない。使用されるＣＰＵが１０台の場合に、使用されるＣＰＵが１台の場合に比べてスループットが２００％しか上昇しない。典型的な解析ジョブは、むしろＣＰＵ＿ＢＯＵＮＤジョブのように機能する。なぜなら、データサーバからの、およびデータサーバへのデータトラフィックのみが別のノードの実行速度に影響を及ぼし、ローカルデータトラフィックは影響を及ぼさないからである。
【００７２】
サーバプログラム
サーバプログラムは、プログラミング言語Ｐｙｔｈｏｎ（非特許文献２３）で実装されている。性能が重視される部分は、Ｃ＋＋で実装されている。ＣおよびＣ＋＋ライブラリ用のＰｙｔｈｏｎ接続は、ＳＷＩＧ  Ｗｒａｐｐｅｒ  Ｇｅｎｅｒａｔｏｒによって生成した。
【００７３】
プログラムＷＯＲＫＦＬＯＷは、タイムアウトに達しているジョブ、または実行したところエラーが返された（状態ＥＲＲＯＲの）ジョブについて、データバンクテーブルｊｏｂを定期的に検査し、それらのジョブを新たに実行させる。さらに、ワークフロー規則を使用し、それにより、場合によっては新規のジョブを生成する。
【００７４】
プログラムＣＯＮＳＵＭＥＲは、各クラスターノード上で実行され、そこで、状態ＮＥＷでジョブを実行する。マルチコアマシン上で、ＣＯＮＳＵＭＥＲを、存在するプロセッサコア１つにつき一回開始させる。
【００７５】
プログラムＷＡＴＣＨＤＯＧは、データバンクおよびデータシステムの問題（機能不良、メモリ領域の不足）を定期的に検査し、場合によっては、電子メールを介してシステム管理者に通知する。
【図１】
0005396484.tif 000011
【図２】
0005396484.tif 000012
【図３】
0005396484.tif 000013
【図４】
0005396484.tif 000014
【図５】
0005396484.tif 000015
【図６】
0005396484.tif 000016
【図７】
0005396484.tif 000017
【図８】
0005396484.tif 000018
【図９】
0005396484.tif 000019
【図１０】
0005396484.tif 000020
【図１１】
0005396484.tif 000021
【図１２】
0005396484.tif 000022
【図１３】
0005396484.tif 000023
【図１４】
0005396484.tif 000024
【図１５】
0005396484.tif 000025