JP3548152B2 - 植物又は植物群の3次元構造測定方法 - Google Patents
植物又は植物群の3次元構造測定方法 Download PDFInfo
- Publication number
- JP3548152B2 JP3548152B2 JP2001349857A JP2001349857A JP3548152B2 JP 3548152 B2 JP3548152 B2 JP 3548152B2 JP 2001349857 A JP2001349857 A JP 2001349857A JP 2001349857 A JP2001349857 A JP 2001349857A JP 3548152 B2 JP3548152 B2 JP 3548152B2
- Authority
- JP
- Japan
- Prior art keywords
- light
- plant
- measurement
- measuring
- reflected
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Description
【発明の属する技術分野】
本発明は、植物又は植物群の3次元構造を測定する方法に関し、詳しくは測定対象である植物又は植物群を破壊することなく非接触で測定するための技術に関する。
【0002】
【従来の技術】
木(植物の一種)や森(植物群の一種)の活力(健康度)を評価する一つの指標として、従来から葉の状態(葉の数、生え方等)や枝の状態(枝の数、枝の伸び方等)が用いられる。このため、従来から木(又は森)を構成する葉や枝の状態(すなわち、木や森の構造)を測定するための種々の技術が開発されている。それら従来の技術の一つとして、本発明者等は、レーザ光切断法を用いた木や森の構造を測定する技術を開発している(Agricultural and Forest Meteorology 91(1998) 149−160)。
この技術では、測定対象物にスリット光を照射する光源と、スリット光の光面に対して光軸が斜交するような位置関係で配置される撮像器とを用いる。具体的には、光源から測定対象物(木や森等)にスリット光を照射し、測定対象物(木や森等)によって反射された光(反射光の輝度)を撮像器で撮像する。そして、撮像器で撮像された画像(スリット光によって形成される光断面像)から、測定対象物(木や森等)の3次元構造を求める。すなわち、光源の位置と撮像器の位置と撮像された画像上の光の位置から、測定対象物(木や森等)の3次元位置を算出する。このようなスリット光の照射、反射光の取得、取得した画像のデータ解析を、スリット光の照射方向を変えながら複数の光切断面に対して行うことで、測定対象物の3次元構造を測定する。
【0003】
【発明が解決しようとする課題】
ところで、木(又は森)は、多数の葉、多数の枝、多数の幹等により構成され、これらの各構成要素(葉、枝・幹)は、それぞれ光の反射特性が異なる性質を有している。したがって、上述した従来の技術を発展させ、測定された光の輝度(光の強度,エネルギ)と測定対象物までの距離から反射特性を算出し、その算出された反射特性から測定対象物が葉なのか枝・幹なのかを判定することが考えられている。
しかしながら、木や森は多数の葉、多数の枝・幹等が何重にも重なって構成される。したがって、測定装置(光源,撮像器)と測定対象物(木や森)との距離は、数メートルから20メートル以上の範囲に及ぶこととなる。このため、測定装置から近くにある反射率の高い対象物の場合と、測定装置から遠くにある反射率の低い対象物の場合では、撮像器で測定される光の輝度が大きく異なることとなる。
このため、現在の撮像器の解像度では、木や森を広範囲に渡って測定(葉・枝・幹等の判定)を行うことは不可能であって、測定対象物の一部しか測定することできなかった。
【0004】
本発明はこのような事情に鑑みてなされたものであって、その目的は、簡易な方法によって、測定対象物(植物又は植物群)を広範囲に渡って計測することができる技術を実現する。
【0005】
【課題を解決するための手段、作用及び効果】
上記課題を解決するため本発明に係る方法は、測定対象物に光を照射する光源と、光源から照射される光に対し光軸が斜交するような位置関係に配置される撮像器とを有する光学式測定装置を用い、光源から植物又は植物群に光を照射し、植物又は植物群から反射される光の輝度を撮像器で測定することで植物又は植物群の3次元構造を測定する方法であり、下記の位置決め工程、第1測定工程、第2測定工程及びデータ処理工程を有する。
位置決め工程は、植物又は植物群に対する光の照射方向が所定の方向となるように光源を位置決めする。第1測定工程は、位置決め工程で位置決めされた光源から光を照射し、第1の測定時間で植物又は植物群から反射される反射光の輝度を測定する。第2測定工程は、位置決め工程で位置決めされた光源から光を照射し第2の測定時間で植物又は植物群から反射される反射光の輝度を測定する。データ処理工程は、第1の測定時間と第2の測定時間のうち短いほうの測定時間で測定された反射光の輝度から測定装置の近くにある植物又は植物群の構造を求めると共に、第1の測定時間と第2の測定時間のうち長いほうの測定時間で測定された反射光の輝度から測定装置の遠くにある植物又は植物群の構造を求める。そして、位置決め工程を繰り返すことで光の照射方向を変えながら、位置決めされた各照射方向について第1測定工程、第2測定工程及びデータ処理工程を行うことで、植物又は植物群の3次元構造を測定する。
【0006】
上記の方法では、光源から照射される光の方向が決まると、その照射方向に対してそれぞれ、第1の測定時間で測定対象物(植物又は植物群)から反射される反射光の輝度と、第2の測定時間で測定対象物(植物又は植物群)から反射される反射光の輝度が測定される。そして、第1の測定時間と第2の測定時間のうち短いほうの測定時間の測定結果から測定装置の近くにある測定対象物(植物又は植物群)の構造を求め、第1の測定時間と第2の測定時間のうち長いほうの測定時間の測定結果から測定装置の遠くにある測定対象物(植物又は植物群)の構造を求める。
すなわち、短い測定時間(撮像器の露出時間が短い)では、撮像器から近い距離にある植物(特に、反射率の高い部分)から反射された光の輝度だけが測定され、これによって測定装置から近い距離にある植物(特に、反射率が高い部分)の構造を求めることができる。一方、長い測定時間(撮像器の露出時間が長い)では、撮像器から近い距離にある植物と遠い距離にある植物の両者から反射された光の輝度が測定される。ここで、長い測定時間とすると(すなわち、撮像器の露出時間を長くすると)、近い距離にある植物(特に、反射率の高い部分)から反射された光の輝度は撮像器の解像度を越えて飽和することとなる(葉と枝・幹の判別が不能となる)。しかしながら、近い距離にある植物(植物群)の構造については、短い測定時間による測定結果から求めることができるため、長い測定時間の測定結果からは遠い距離にある植物(特に、反射率の低い部分)の構造について求めることができれば良い。
したがって、上記の方法では測定時間を変えることによって、測定装置と測定対象物(木や森等)との距離が広範囲に渡っても測定対象物の構造を求めることができる。
【0007】
ここで、上記「第1の測定時間」と「第2の測定時間」は、測定装置から測定対象物までの距離(測定範囲)や撮像器の解像度に応じて適宜設定することができる。例えば、一桁程度異なる測定時間とすることで広範囲に渡って測定対象物を測定することができる。また、測定対象物との距離(測定範囲)によっては、第3の測定時間を設定して計測を行うようにしても良い。すなわち、測定回数は、測定対象物、撮像器の解像度等に応じて種々に変更することができる。
また、「光源」としては、光源からの光の照射角、反射光の撮像器への入射角から測定対象物の3次元位置が特定できるものであればどのようなものでも良く、例えば、スリット状(直線状)の光を照射するスリット光源や点状の光を照射する光源等を用いることができる。
【0008】
前記光源は少なくとも第1の波長の光と第2の波長の光が照射可能となっており、前記の第1測定工程と第2測定工程では、第1の波長と第2の波長の光をそれぞれ照射し、それぞれの場合について反射される反射光の輝度を測定することが好ましい。
このような構成によると、第1の波長の光を照射した場合と第2の波長の光を照射した場合における反射光の輝度をそれぞれ測定することで、測定対象物(植物又は植物群)の構造をより正確に測定することができる。
なお、第1の波長の光として600〜700nmの波長帯域の光(赤色光)を用い、第2の波長の光として700〜900nmの波長帯域の光(近赤外)を用いることが好ましい。植物(詳しくは、植物を構成する葉)の反射特性は、赤色と近赤外の光のときに大きな違いがあり、この二つの波長帯域の光を利用することで植物の構造をより正確に測定することができる。
【0009】
前記第1測定工程で測定される反射光の輝度は、光源から光を照射しない状態で第1の測定時間に植物又は植物群から反射される反射光の輝度によって補正され、前記第2測定工程で測定される反射光の輝度は、光源から光を照射しない状態で第2の測定時間に植物又は植物群から反射される反射光の輝度によって補正されていることが好ましい。
このような構成によると、測定装置の光源以外の外光(例えば、月光、屋外広告灯等)の影響を取除くことができ、より正確に測定対象物(植物又は植物群)の構造を測定することができる。
また、本願は上述した植物又は植物群の3次元構造測定方法を好適に実施するための測定装置を提供する。すなわち、本発明に係る測定装置は、植物又は植物群の3次元構造を測定する装置であり、測定対象物に光を照射する光源と、光源から照射される光に対し光軸が斜交するような位置関係に配置される撮像器と、植物又は植物群に対する光の照射方向が所定の方向となるように光源を位置決めする位置決め手段と、位置決め手段で位置決めされた光源から光を照射させ、第1の測定時間で植物又は植物群から反射される反射光の輝度を撮像器に測定する第1測定手段と、位置決め手段で位置決めされた光源から光を照射させ、第2の測定時間で植物又は植物群から反射される反射光の輝度を撮像器に測定する第2測定手段と、第1の測定時間と第2の測定時間のうち短いほうの測定時間で測定された反射光の輝度から測定装置の近くにある植物又は植物群の構造を求めると共に、第1の測定時間と第2の測定時間のうち長いほうの測定時間で測定された反射光の輝度から測定装置の遠くにある植物又は植物群の構造を求めるデータ処理手段と、前記位置決め手段を駆動して光の照射方向を変えながら、位置決めされた各照射方向について第1測定手段、第2測定手段及びデータ処理手段を実行させる手段と、を有する。
この測定装置によると、同一の光照射方向について測定時間を変えて植物又は植物群からの反射光の輝度を測定するため、測定装置と測定対象物(木や森等)との距離が広範囲に渡っても測定対象物の構造を求めることができる。
【0010】
【発明の実施の形態】
本発明は、下記に示す形態によって好適に実施することができる。
(形態1) 第1の波長の光は685nmの波長の光(赤色レーザ光)であり、第2の波長の光は830nm(近赤外レーザ光)である。木や森等の一構成要素である葉は、赤色レーザ光に対する反射率が低く、近赤外レーザ光に対する反射率が高い。一方、木や森等の他の構成要素である枝・幹等は、赤色レーザ光に対する反射率が高く、近赤外レーザ光に対する反射率も高い。したがって、赤外レーザ光と近赤外レーザ光を用いて測定することで、測定対象物が葉であるのか、それとも枝・幹であるのかの判別を正確に行うことができる。
(形態2) 測定装置の光源からはスリット状の光が照射される。スリット状の光を照射することで、点状の光を照射する場合に比較して、光の照射・反射光の測定の回数を少なくすることができる。
【0011】
【実施例】
本発明の測定方法に係る具体的な実施例について図面を参照して説明する。図1に本発明の測定方法を実施するための測定装置の構成を模式的に示している。
図1に示すように、本実施例に係る測定装置は、プラットホーム30上の測定ヘッド(図示省略)に取付けられるレーザユニット10及びCCDカメラ20と、CCDカメラ20で撮影された画像データを記憶すると共に、レーザユニット10,CCDカメラ20の動作を制御する制御装置24とで構成される。
レーザユニット10は、測定対象である植物群に対しスリット光を照射するレーザラインマーカ12(スリットレーザ)を備える。レーザラインマーカ12は、波長685nmの赤色レーザ光と波長830nmの近赤外レーザ光のいずれかを選択的に出射可能となっており、制御装置24の指示に応じて一方の波長のレーザ光を出射する。また、レーザラインマーカ12は、ポテンションメータ14及びステッピングモータ16を介してプラットホーム30の測定ヘッドに取付けられており、ステッピングモータ16が回転することでレーザラインマーカ12から照射されるレーザ光の照射方向が変わるようになっている。ステッピングモータ16の回転角度はポテンションメータ14によって検出され、ポテンションメータ14で検出された回転角度は制御装置24に出力されるようになっている。制御装置24は、ポテンションメータ14から出力される検出信号に基づいてステッピングモータ16を駆動することで、レーザラインマーカ12から照射されるレーザ光の照射角度が所望の角度となるように制御する。
CCDカメラ20は、測定対象物から反射された光の像(反射輝度)を撮像する装置であって、プラットホーム30の測定ヘッドに対して撮影方向を調整可能に取付けられている。CCDカメラ20とレーザユニット10間の距離は、測定前に予め所定の間隔となるように調整され、また、CCDカメラ20の撮影方向は、CCDカメラ20の光軸とレーザユニット10から照射されるスリット光の光面とが斜めに交わるよう(斜交するよう)に調整されている。CCDカメラ20には制御装置24が接続されており、CCDカメラ20で撮像された画像データは制御装置24に記憶される。
制御装置24は、汎用のパーソナルコンピュータであり、レーザユニット10の動作を制御することでレーザ光の照射方向を調整し、同時に、CCDカメラ20で撮影された画像データを制御装置内のハードディスクに格納する処理と、制御装置24のハードディスクに格納された画像データを処理することで植物群の3次元構造を算出する処理を行う。
【0012】
次に、上述した測定装置を用いた測定方法を説明する。なお、本発明に係る測定方法では、レーザユニット10以外の光源による光の影響を少なくするため、測定対象物の測定は夜間に行われる。
まず、図2〜図4を用いて、レーザユニット10とCCDカメラ20の配置からCCDカメラ20による画像の撮像、撮像した画像データから測定対象物の3次元位置の算出までの手順を説明する。ここで、図2はレーザユニットから照射されたスリットレーザ光により測定対象物に描かれる光切断像をCCDカメラで撮像する状態を模式的に示す図であり、図3はCCDカメラで撮像された画像から測定対象物の3次元的位置を算出する際の原理を説明するための図であり、図4はレーザユニットから照射されるスリットレーザの照射角を変えながら測定対象物の全体を測定する際の手順を説明するための図である。
上述した測定装置によって測定対象物(植物群)を測定するには、まず、図2に示すように、測定対象物Bに対してスリット光を照射するレーザユニット10を配置する。次いで、レーザユニット10から測定対象物Bに対して照射されるスリットレーザ光の光面とCCDカメラ20の光軸(撮像方向)が斜交するようにCCDカメラ20を配置する。なお、CCDカメラ20は、測定対象物Bの全体が撮像できる位置に配置され、測定中にその光軸(撮像方向)が変更されることはない。また、既に説明したようにレーザユニット10とCCDカメラ20は、プラットホーム30上に所定の間隔を設けて取付けられているため、レーザユニット10とCCDカメラ20の間隔(位置)が測定中に変更されることはない。
【0013】
上述のようにしてレーザユニット10とCCDカメラ20の位置が位置決めされると、レーザユニット10からスリット光Lが照射される。レーザユニット10から照射されたスリット光Lの光面に測定対象物Bがあると、測定対象物Bから光が反射され、これによって測定対象物Bの上に光切断像Tが描かれる。この光切断像T(測定対象物Bにより反射された光)がCCDカメラ20で撮像され、撮像された画像データは制御装置24のハードディスクに記憶される。このようなスリット光Lの照射と、測定対象物B上に描かれた光切断像Tの撮像とを、スリット光Lの照射方向を変えながら行うことで、測定対象物Bの全体について光切断像Tを撮像する。すなわち、図4に示すように、ステッピングモータ16を駆動することでレーザラインマーカ12から照射されるレーザ光の照射方向をL1→L2→L3→L4→L5・・と変えながら、各照射方向(レーザユニット10の各走査角)における光切断像TをCCDカメラ20で撮像し、撮像した画像を制御装置24のハードディスクに記憶する。
【0014】
レーザユニット10の各照射方向(各走査角)における光切断像Tが撮像されると、撮像された各画像から測定対象物Bの3次元の位置座標が算出される。3次元の位置座標を算出する原理を図3に基づいて説明する。図3には、A点にレーザユニット10(光源)が、C点にCCDカメラ20(撮像器)が配置され、測定対象物Bの座標(x,y,z)を算出する場合を示している。
図3に示すように、レーザユニット10(光源)から照射されるスリット光の照射角θはポテンションメータ14によって検出されて既知であり、また、レーザユニット10とCCDカメラ20はプラットホーム30上に取付けられて両者の距離dは既知である。さらに、CCDカメラ20で撮像した画像内における座標値(画像内のピクセルの座標値)から、CCDカメラ20から見た対象物Bの方向〔すなわち、反射光のCCDカメラ20への入射角(図3に示すαとβで規定される)〕を求めることができる。このように3つの角度θ,α,βと、レーザユニット10とCCDカメラ20間の距離dが分かると、図3に示す幾何学的関係から以下に示す関係式が成立する。ここで、測定対象物Bの位置座標は(x,y,z)とする。
【0015】
【式1】
【0016】
したがって、上記の連立方程式を解くことによって、測定対象物Bの位置(x,y,z)が算出される。このような計算を撮像した画像内の全てのピクセルについて行うことで画像内の各点の3次元の位置座標が算出され、その点の反射輝度から各位置に何があるのか(測定対象物の3次元構造)を特定することができる。
【0017】
上述したレーザユニット10のレーザ照射角を所定の角度に位置決めするレーザユニット10の位置決め動作と、レーザユニット10から照射されたレーザ光による光切断像のCCDカメラ20による撮像と、撮像された画像の制御装置24のハードディスクへの記憶は、制御装置24がレーザユニット10及びCCDカメラ20を制御することにより自動的に行われる。以下、測定対象物測定時における制御装置24の処理フローを図5,6を用いて説明する。ここで、図5は制御装置24における画像撮影処理全体の手順を示すフローチャートであり、図6は撮影した画像のデータ取込処理の手順を示すフローチャートである。
図5に示すように、制御装置24は、まずレーザユニット10の照射方向が所定の角度となるようレーザユニット10を位置決めする(S01)。具体的には、ステッピングモータ16を駆動することでレーザラインマーカ12を回転させる。この際、ポテンションメータ14でステッピングモータ16の回転角を検出し、検出されたステッピングモータ16(レーザラインマーカ12)の回転角から、レーザラインマーカ12が所定の角度となる位置でステッピングモータ16の駆動を停止する。なお、レーザユニット10の位置決め(ただし、最初の走査線上への位置決めのみ)はレーザ光を照射しながら行われ、これにより測定者は意図した測定範囲にレーザ光が照射されるか否かを確認することができる。
次に、CCDカメラ20による測定対象物の撮影時間(露出時間)を決定し、ステップS03で行われるCCDカメラ20の撮影のためにレーザユニット10から照射されるレーザをOFFする(S02)。すなわち、本実施例では、測定対象物を撮影する撮影時間(CCDカメラ20の露出時間)を予め設定された複数の条件(本実施例では2条件)に変えて、それぞれの撮影時間で測定対象物を撮影する。このため、CCDカメラ20の露出時間を予め設定された複数の露出時間の中の一つに決定する。また、ステップS03では測定対象物にレーザが照射されていない状態で画像を撮影するため、ステップS02ではレーザユニット10をOFFする。
【0018】
次いで、制御装置24はCCDカメラ20にデータ取込指令を出力することでデータ取込処理を行う(S03)。このデータ取込処理については図6を用いて説明する。図6に示すように、データ取込処理では、まず制御装置24はCCDカメラ20に対してデータ取込指令を出力する(S11)。この制御装置24からCCDカメラ20に出力されるデータ取込指令には、図5のステップS02で決定された露出時間が含まれている。
データ取込指令がCCDカメラ20に出力されると、CCDカメラ20は、測定対象物の撮影を開始し、撮影開始からの経過時間がデータ取込指令で指示された露出時間となったか否かを判断する(S12)。
そして、指示された露出時間が経過していない場合〔ステップS12でNOの場合〕にはCCDカメラ20による撮影を続行する。一方、指示された露出時間が経過している場合〔ステップS12でYESの場合〕には撮影(データ取込)を終了する(S13)。したがって、このステップS03のデータ取込処理によって、レーザユニット10からレーザ光が照射されていない状態でステップS02で設定された露出時間だけ測定対象物が撮影され、撮影された画像データが取込まれることとなる。なお、CCDカメラ20により撮影された画像データは制御装置24に転送され、制御装置24のハードディスク上の所定のアドレスに格納される。
【0019】
上述のステップS03のデータ取込処理が終了すると、次に、制御装置24はレーザユニット10から照射されるレーザ光の波長を設定する(S04)。すなわち、本実施例では、レーザユニット10から波長685nmの赤色レーザ光と、波長830nmの近赤外レーザ光を交互に照射し、それぞれのレーザ光照射時において測定対象物の反射光の輝度をCCDカメラ20で撮像する。したがって、このステップS04では、どちらの波長のレーザ光を測定対象物に照射するかを設定する。
ステップS04でレーザユニット10から照射するレーザ光の波長が決まると、制御装置24はレーザユニット10にステップS04で決定された波長のレーザ光を測定対象物に照射するよう指令を出力する(S05)。これにより、レーザユニット10からは、ステップS04で決定された波長のレーザ光が照射される。
レーザ光が照射が開始され、照射されるレーザ光の出力が安定すると、次に、前述したデータ取込処理(図6参照)を行う(S06)。したがって、このステップS06の処理では、ステップS04で決定された波長のレーザ光がレーザユニット10から測定対象物に照射された状態で、ステップS02で決定された露出時間だけ測定対象物を撮像する。
【0020】
データ取込処理が終了すると、レーザユニット10から照射されるレーザをオフし(S07)、次いで、全ての波長についてステップS06のデータ取込処理を行ったか否かを判断する(S08)。すなわち、本実施例では、赤色レーザと近赤外レーザの2種類のレーザ光を測定対象物に照射し、それぞれの場合についてCCDカメラ20で画像を撮像する。したがって、2種類の波長のレーザ光を照射したそれぞれの場合について測定対象物を撮像しているか否かを判定する。
全ての波長について画像を撮像していない場合〔ステップS08でNOの場合〕には、ステップS04に戻ってステップS04からの処理を繰返す。これによって、レーザユニット10から照射されるレーザ光の波長が変更されて、測定対象物の画像が撮像される。
一方、全ての波長について画像を撮影している場合〔ステップS08でYESの場合〕にはステップS09に進んで、撮像しなければならない全ての露出時間で測定対象物を撮像したか否かが判断される(S09)。すなわち、本実施例では、測定対象物を予め設定された複数の露出時間のそれぞれの条件で撮影する。したがって、ステップS09では、全ての露出時間(本実施例では2条件)で測定対象物を撮影したか否かを判断する。
【0021】
全ての露出時間で測定対象物を撮影していない場合〔ステップS09でNOの場合〕には、ステップS02に戻ってステップS02からの処理を繰返す。したがって、ステップS02で露出時間が変更され、変更された露出時間についてステップS03からの処理が行われる。一方、全ての露出時間で測定対象物を撮影している場合〔ステップS09でYESの場合〕にはステップS10に進んで、予め設定された全てのレーザ照射角(ステッピングモータ16の走査角)についてステップS01〜ステップS09までの処理が終了しているか否かが判断される(S10)。
全てのレーザ照射角について処理が終了していない場合〔ステップS10でNOの場合〕は、ステップS01に戻ってステップS01からの処理が繰返される。したがって、レーザユニット10のレーザ照射角が変更されて、ステップS02からの処理が繰返される。一方、全てのレーザ照射角について処理が終了している場合〔ステップS10でYESの場合〕には画像撮影処理を終了する。
【0022】
上述の説明から明らかなように、本実施例に係る方法では、1つのスリット光(レーザユニット10の1つの照射方向)について、各露出時間毎に、(1)レーザ光が照射されていない状態、(2)赤色レーザ光が照射されている状態、(3)近赤外レーザ光が照射されている状態の3条件で測定対象物が撮影される。なお、本実施例では、露出時間は2条件であるので、レーザユニット10の1つの照射方向について撮影条件を変えて計6枚の画像が撮影され、次に説明する画像解析処理が行われる。
【0023】
制御装置24で行われる画像解析処理について図7に示すフローチャートを参照して説明する。
図7に示すように画像解析処理では、まず、解析対象となる画像を選択するための条件として、まず、撮影した画像の各レーザ照射角(ステッピングモータ16の走査角)の中から一つのレーザ照射角を選択する(S14)。すなわち、図4に示すようにレーザユニット10でL1〜L5・・とレーザ照射角を変えて測定対象物を撮影しているので、これらのレーザ照射角L1〜L5・・の一つを選択する。
レーザ照射角が選択されると、そのレーザ照射角で撮影された画像内の各点の座標(2次元座標)を3次元の位置座標に変換する(S15)。すなわち、既に説明した式1の連立方程式を解くことによって、画像内の各点の座標を3次元の位置座標に変換する(図3参照)。
【0024】
次に、撮影条件の一つであるCCDカメラ20の露出時間を選択する(S16)。具体的には、撮影した2つの露出時間の中の一つの露出時間を選択する。ステップS14とステップS16により解析対象となる画像のレーザ照射角と露出時間が決定されると、まずは、当該レーザ照射角及び当該露出時間で撮影した各画像(具体的には、レーザ非照射時,赤色レーザ照射時,近赤外レーザ照射時の3画像)の中からレーザ非照射時の画像データを制御装置24のハードディスクから読み込む(S17)。
次に、照射されたレーザ光の波長の一つを選択する(S16)。すなわち、赤色レーザ照射時か近赤外レーザ照射時のいずれかを選択する。
レーザ光の波長が選択されると、その波長のレーザ光を照射した時の画像データ(ステップS14で選択されたレーザ照射角でレーザが照射され、ステップS16で選択された露出時間で撮影した画像のデータ)を制御装置24のハードディスクより読み込む(S19)。
次に、ステップS19で読み込んだ画像データ(ステップS18で選択された波長のレーザ光照射時の画像データ)の各点の反射輝度を、ステップS16で読み込んだ画像データ(レーザ非照射時の画像データ)の対応する点の反射輝度によって修正する(S20)。すなわち、レーザ照射時の画像データの反射輝度からレーザ非照射時の画像データの反射輝度を減算することによって、レーザ以外の光源(例えば、月光、屋外広告灯等の外光)による影響を除去する。
各画像データの反射輝度が修正されると、次に、その修正された各点の反射輝度をステップS15で算出された3次元の位置座標とを関連付けて記憶する(S21)。これによって、各位置座標(x,y,z)における反射輝度とが関連付けられて記憶されることなる。
【0025】
ステップS22では、ステップS14で選択されたレーザ照射角で、かつ、ステップS16で選択された露出時間で撮影された全ての波長の画像データについて、上記ステップS21の処理が行われたか否かを判断する(S22)。全ての波長(本実施例では、赤色レーザ光と近赤外レーザ光の2種類)についてステップS21の処理が終了していない場合〔ステップS22でNOの場合〕には、ステップS18に戻ってステップS18からの処理を繰り返す。これによって、ステップS14で選択されたレーザ照射角で、かつ、ステップS16で選択された露出時間で撮影された赤色レーザ光照射時の画像データと近赤外レーザ光照射時の画像データについて、3次元の位置座標とその位置座標における反射輝度とが関連付けられて記憶される。
全波長についてステップS21の処理が終了している場合〔ステップS22でYESの場合〕には、ステップS14で選択されたレーザ照射角について,全ての露出時間の画像データ(すなわち、短露出時間と長露出時間)についてデータ処理が終了しているか否かを判断する(S23)。全ての露出時間の画像データについてデータ処理が終了していない場合〔ステップS23でNOの場合〕には、ステップS16に戻ってステップS16からの処理を繰り返す。これによって、ステップS14で選択されたレーザ照射角で撮影された全ての画像データについてステップS16からの処理が実行される。
全ての露出時間についてデータ処理が終了している場合〔ステップS22でYESの場合〕には、次に、上述したデータ処理が全てのレーザ照射角の画像データについて終了しているか否かかを判断する(S24)。全てのレーザ照射角の画像データについてデータ処理が終了していない場合〔ステップS24でNOの場合〕には、ステップS14に戻ってステップS14からの処理を繰り返す。これによって、撮影した全画像データの全てについてステップS21の処理が実行されることとなる。
【0026】
全ての画像データについてデータ処理が終了している場合〔ステップS24でYESの場合〕には、ステップS21で格納された位置座標(x,y,z)と、その位置座標の各反射輝度(赤色レーザ光の反射輝度,近赤外レーザ光の反射輝度)から、その位置に葉があるのか枝・幹があるのかを判別する(S25)。すなわち、葉がある場合は、赤色レーザ光の反射輝度は低く、近赤外レーザ光の反射輝度は高くなる。一方、枝・幹がある場合は、赤色レーザ光と近赤外レーザ光の反射輝度はともに高くなる。また、葉・枝・幹のいずれもが無い場合は、赤色レーザ光と近赤外レーザ光の反射輝度はいずれも低くなる。したがって、これらの特徴から、各位置に葉があるのか枝・幹があるのかを判別する。これによって、測定対象物(植物又は植物群)の3次元構造(葉と枝・幹の判別)が測定される。
なお、衛星画像解析の分野で提案されているNDVI(正規化植生指標)を用いて、求めた反射輝度から植物(植物群)の分光反射特性を評価することもできる。具体的には、露出時間を短くして撮影した画像データの反射輝度をその露出時間で除算することで正規化し、また、露出時間を長くして撮影した画像データの反射輝度をその露出時間で除算することで正規化する。そして、この正規化した反射輝度でNDVIを求め、測定した植物(植物群)を総合的に評価する。NDVIは、次に示す式2であらわされる。
【0027】
【式2】
【0028】
式2中、Lrefはレーザ光非照射時の正規化された反射輝度を示し、Lredは赤色レーザ光照射時の正規化された反射輝度を示し、Lirは近赤外レーザ光照射時の正規化された反射輝度を示している。
上述した葉と枝・幹の相違による反射特性の相違から明らかなように、NDVIの値が大きくなると、その位置に葉が存在すると判定することができる。また、NDVIの値が小さくなると、その位置には葉が存在せず、枝・幹又は植物が存在しないと判定することができる。したがって、NDVIの値によって植物(又は植物群)の葉の状態が判定することができ、これによって植物の活性度(健康度)を判定することができる。
【0029】
以上本実施例に係る測定装置の構造及び動作を説明してきたが、次に、上述の測定装置を用いて森林(植物群)を測定した例について説明する。まず、本実施例に係る測定装置を用いた森林測定の応用例について図8を参照して説明する。図8は、測定装置を用いた森林の測定の応用例を模式的に示した図である。
図8(a)に示す例は、森林内にレーザユニット10とCCDカメラ20を地面上に、地面と水平方向を撮影するように配置し、森林の低層部分を測定した例である。この測定によって森林の低層部の植物群の状態(樹木の植生,低層木の活性度)を判定することができる。図8(b)に示す例は、森林内にレーザユニット10とCCDカメラ20を地面から上方に向って撮影するように配置し、森林の上層部分を測定した例である。この測定によって森林内における高さ方向の葉の状態を判定することができる。図8(c)に示す例は、森林の上方にレーザユニット10とCCDカメラ20をセットし、森林の林冠部を測定した例である。この測定によって林冠部の葉の状態を判定することができる。
図9及び図10は、実際に測定した実験結果を示している。すなわち、森林の林床(地面)に測定装置を真上に向けて設置した場合(図8(b)の場合)の測定結果を3次元的に横から示した図である。図中左側に落葉樹(コナラ)が、右側に常緑樹(クスノキ)が示されており、植生指標(NDVI)が高いほど濃い色で示されている。各図の測定条件は、それぞれ、図9は落葉期に露出時間0.03秒(単一露出時間)で測定した例であり、図10は落葉期に露出時間0.063秒(単一露出時間)で測定した例である。
図9,図10から明らかなように、露出時間が短いと測定装置から距離が近い植物群のみが観測され、一方露出時間が長くなると測定装置から距離が離れた植物群までが観測されている。ただし、露出時間が長くなると、測定装置から距離が近い植物群の植生指標が全て飽和し、葉と枝・幹の区別が不能となっている(一方、図9では測定装置から近い植物群の濃淡が表れている)。したがって、露出条件を変えて複数の露出時間で植物群を測定することで、植物群の3次元構造を広範囲に渡って測定できることが確認された。
【0030】
上述の説明から明らかなように、本実施例に係る測定方法では、CCDカメラ20の露出時間を変えて複数条件で撮影することで、短露出時間の測定によって測定装置から近い距離の植物(又は植物群)の3次元構造を測定することができ、長露出時間の測定によって測定装置から遠い距離の植物(又は植物群)の3次元構造を測定することができる。
また、赤色レーザ光と近赤外レーザ光を照射して得られた輝度を、レーザ非照射時に得られた輝度で修正するため、レーザ以外の光源(月光、屋外広告灯等)の影響を除去することができる。
【0031】
以上、本発明の具体例を詳細に説明したが、これらは例示にすぎず、特許請求の範囲を限定するものではない。特許請求の範囲に記載の技術には、以上に例示した具体例を様々に変形、変更したものが含まれる。
また、本明細書または図面に説明した技術要素は、単独であるいは各種の組み合わせによって技術的有用性を発揮するものであり、出願時請求項記載の組み合わせに限定されるものではない。また、本明細書または図面に例示した技術は複数の目的を同時に達成するものであり、そのうちの一つの目的を達成すること自体で技術的有用性を持つものである。
【図面の簡単な説明】
【図1】本発明の測定方法を実施するための測定装置の構成を模式的に示す図。
【図2】レーザユニットから照射されたスリットレーザ光による測定対象物の光切断像をCCDカメラで撮像する状態を模式的に示す図。
【図3】CCDカメラで撮像された画像から測定対象物の3次元的位置を算出する際の原理を説明するための図。
【図4】レーザユニットから照射されるレーザ光の照射角を変えながら測定対象物の全体を測定する際の手順を説明するための図。
【図5】画像撮影処理の手順を示すフローチャート。
【図6】データ取込処理の手順を示すフローチャート。
【図7】画像解析処理の手順を示すフローチャート。
【図8】本実施例に係る測定装置による森林測定の応用例を模式的に示す図。
【図9】測定結果の一例を示す図。
【図10】測定結果の一例を示す図。
【符号の説明】
10・・レーザユニット
12・・レーザラインマーカ
14・・ポテンションメータ
16・・ステッピングモータ
20・・CCDカメラ
24・・制御装置(パーソナルコンピュータ)
30・・プラットホーム
Claims (4)
- 測定対象物に光を照射する光源と、光源から照射される光に対し光軸が斜交するような位置関係に配置される撮像器とを有する光学式測定装置を用い、光源から植物又は植物群に光を照射し、植物又は植物群から反射される光の輝度を撮像器で測定することで植物又は植物群の3次元構造を測定する方法であり、
植物又は植物群に対する光の照射方向が所定の方向となるように光源を位置決めする位置決め工程と、
その位置決め工程で位置決めされた光源から光を照射し、第1の測定時間で植物又は植物群から反射される反射光の輝度を測定する第1測定工程と、
前記位置決め工程で位置決めされた光源から光を照射し、第2の測定時間で植物又は植物群から反射される反射光の輝度を測定する第2測定工程と、
第1の測定時間と第2の測定時間のうち短いほうの測定時間で測定された反射光の輝度から測定装置の近くにある植物又は植物群の構造を求めると共に、第1の測定時間と第2の測定時間のうち長いほうの測定時間で測定された反射光の輝度から測定装置の遠くにある植物又は植物群の構造を求めるデータ処理工程と、を有し、
前記位置決め工程を繰り返すことで光の照射方向を変えながら、位置決めされた各照射方向について第1測定工程、第2測定工程及びデータ処理工程を行うことで、植物又は植物群の3次元構造を測定することを特徴とする植物又は植物群の3次元構造測定方法。 - 前記光源は少なくとも第1の波長の光と第2の波長の光が照射可能となっており、前記の第1測定工程と第2測定工程では、第1の波長と第2の波長の光をそれぞれ照射し、それぞれの場合について反射される反射光の輝度を測定することを特徴とする請求項1に記載の植物又は植物群の3次元構造測定方法。
- 前記第1測定工程で測定される反射光の輝度は、光源から光を照射しない状態で第1の測定時間に植物又は植物群から反射される反射光の輝度によって補正され、前記第2測定工程で測定される反射光の輝度は、光源から光を照射しない状態で第2の測定時間に植物又は植物群から反射される反射光の輝度によって補正されていることを特徴とする請求項1又は2に記載の植物又は植物群の3次元構造測定方法。
- 植物又は植物群の3次元構造を測定する装置であり、
測定対象物に光を照射する光源と、
光源から照射される光に対し光軸が斜交するような位置関係に配置される撮像器と、
植物又は植物群に対する光の照射方向が所定の方向となるように光源を位置決めする位置決め手段と、
位置決め手段で位置決めされた光源から光を照射させ、第1の測定時間で植物又は植物群から反射される反射光の輝度を撮像器に測定する第1測定手段と、
位置決め手段で位置決めされた光源から光を照射させ、第2の測定時間で植物又は植物群から反射される反射光の輝度を撮像器に測定する第2測定手段と、
第1の測定時間と第2の測定時間のうち短いほうの測定時間で測定された反射光の輝度から測定装置の近くにある植物又は植物群の構造を求めると共に、第1の測定時間と第2の測定時間のうち長いほうの測定時間で測定された反射光の輝度から測定装置の遠くにある植物又は植物群の構造を求めるデータ処理手段と、
前記位置決め手段を駆動して光の照射方向を変えながら、位置決めされた各照射方向について第1測定手段、第2測定手段及びデータ処理手段を実行させる手段と、
を有する植物又は植物群の3次元構造測定装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2001349857A JP3548152B2 (ja) | 2001-11-15 | 2001-11-15 | 植物又は植物群の3次元構造測定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2001349857A JP3548152B2 (ja) | 2001-11-15 | 2001-11-15 | 植物又は植物群の3次元構造測定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2003143961A JP2003143961A (ja) | 2003-05-20 |
JP3548152B2 true JP3548152B2 (ja) | 2004-07-28 |
Family
ID=19162464
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2001349857A Expired - Fee Related JP3548152B2 (ja) | 2001-11-15 | 2001-11-15 | 植物又は植物群の3次元構造測定方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3548152B2 (ja) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3728511B2 (ja) * | 2003-07-07 | 2005-12-21 | 国立大学法人名古屋大学 | 移動式3次元構造測定装置 |
JP4599590B2 (ja) * | 2005-05-10 | 2010-12-15 | 独立行政法人農業・食品産業技術総合研究機構 | 植物の生育度測定装置 |
JP2007171033A (ja) * | 2005-12-22 | 2007-07-05 | Univ Nagoya | 葉面積指数の間接測定方法および間接測定システム |
JP5180772B2 (ja) * | 2008-10-03 | 2013-04-10 | 本田技研工業株式会社 | 穂の構造の解析方法および装置 |
US9207072B2 (en) * | 2010-12-02 | 2015-12-08 | Nec Corporation | Leaf area index measurement system, device, method, and program |
JP2012189535A (ja) * | 2011-03-14 | 2012-10-04 | Ihi Corp | 植生検出装置及び植生検出方法 |
KR101319011B1 (ko) * | 2011-12-20 | 2013-10-17 | 충북대학교 산학협력단 | 레이저 스캐너를 이용한 폭 및 직경 측정장치 및 측정방법 |
CN104884937B (zh) * | 2012-12-26 | 2019-01-08 | 日本电气株式会社 | 图像测量方法、系统、设备 |
JP6489604B2 (ja) * | 2014-09-04 | 2019-03-27 | 国立研究開発法人農業・食品産業技術総合研究機構 | 植物群落構造解析システム |
JP6338184B2 (ja) * | 2014-09-04 | 2018-06-06 | 国立研究開発法人農業・食品産業技術総合研究機構 | 植物栽培用の環境制御システム |
JP6895905B2 (ja) * | 2018-01-18 | 2021-06-30 | ヤンマーパワーテクノロジー株式会社 | 植物長尺部検出システム |
CN109374856B (zh) * | 2018-09-25 | 2020-04-07 | 大连理工大学 | 观测透明土模型内部三维空间变形的试验装置及使用方法 |
JP7290415B2 (ja) * | 2018-12-06 | 2023-06-13 | 三星電子株式会社 | 三次元計測装置および加熱調理器 |
CN112540060A (zh) * | 2020-12-03 | 2021-03-23 | 吉林农业大学 | 一种归一化植被指数获取装置及方法 |
JP7527005B2 (ja) | 2020-12-23 | 2024-08-02 | 国立研究開発法人農業・食品産業技術総合研究機構 | 植物体撮像装置 |
-
2001
- 2001-11-15 JP JP2001349857A patent/JP3548152B2/ja not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP2003143961A (ja) | 2003-05-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3548152B2 (ja) | 植物又は植物群の3次元構造測定方法 | |
CN105286785B (zh) | 多光谱医学成像装置及其方法 | |
CN106546568B (zh) | 一种获取植物三维叶绿素荧光图像信息的方法及装置 | |
EP3938802A1 (en) | Detector for identifying at least one material property | |
Haala et al. | Combination of terrestrial laser scanning with high resolution panoramic images for investigations in forest applications and tree species recognition | |
Ma et al. | Determining woody-to-total area ratio using terrestrial laser scanning (TLS) | |
US20110134114A1 (en) | Depth-based gain control | |
JP6454489B2 (ja) | 観察システム | |
JP5133626B2 (ja) | 表面反射特性測定装置 | |
US6907138B1 (en) | Method and device for analyzing surface structure | |
JP2009511897A (ja) | 対象物再構成方法およびシステム | |
US20190124914A1 (en) | Improvements in or relating to optical remote sensing systems for aerial and aquatic fauna, and use thereof | |
JP2003009664A (ja) | 作物生育量測定装置、作物生育量測定方法、作物生育量測定プログラム及びその作物生育量測定プログラムを記録したコンピュータ読取可能な記録媒体 | |
AU2021204034B2 (en) | Information processing device, information processing method and program | |
JP4507806B2 (ja) | 指紋画像撮像装置 | |
CA2315357A1 (fr) | Dispositif d'analyse non destructive de plantes et vehicule comportant un tel dispositif embarque | |
JP2004309240A (ja) | 3次元形状測定装置 | |
CN108154126A (zh) | 虹膜成像系统及方法 | |
JP2006285763A (ja) | 被写体についての陰影のない画像を生成する方法および装置、並びにそれに用いる白色板 | |
US20040120553A1 (en) | Device for contactless optical acquisition of biometric characteristics of at least one body part | |
JP2006106979A (ja) | 指紋画像撮像装置 | |
KR101341189B1 (ko) | 표면각도 측정 광센서 및 이를 이용한 감지 방법 | |
JP3638569B2 (ja) | 3次元計測装置 | |
ES2966765T3 (es) | Procedimiento, sistema y programas de ordenador para el recuento automático del número de insectos en una trampa | |
JP6542173B2 (ja) | 角度測定装置およびその作動方法並びに角度測定プログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A712 Effective date: 20031031 |
|
RD03 | Notification of appointment of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7423 Effective date: 20031210 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20040121 |
|
RD02 | Notification of acceptance of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7422 Effective date: 20040121 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20040121 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20040330 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20040401 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20040415 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313115 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100423 Year of fee payment: 6 |
|
LAPS | Cancellation because of no payment of annual fees |