(57)【要約】【課題】木材の年輪画像等から取得した濃度波形信号等のように波形の特徴点の値または特徴点間距離が不規則な波形信号であっても、波形の特徴点（ピーク点またはトラフ点）を迅速かつ高精度に検出することができる計測方法および装置を提供する。【解決手段】所定のマザーウェーブレットと複数段階のスケールレベルを用いて所定区間内の波形信号をウェーブレット変換し、このウェーブレット変換により生成された各スケールレベル毎のウェーブレット変換信号に対し、該スケールレベルに応じた区間幅毎に区間二乗平均を算出し、前記所定区間内の任意の点において前記区間二乗平均の算出値が最大となるウェーブレット変換信号のスケールレベルをその点における支配レベルとし、この支配レベルに応じた区間幅毎に前記波形信号が最大値または最小値となる点を波形信号の特徴点として検出する。【選択図】  図２
2004164346.tif 000002
【特許請求の範囲】
【請求項１】
波形上の特徴点の値または特徴点間距離が不規則な波形信号の特徴点を検出する方法であって、所定のマザーウェーブレットと複数段階のスケールレベルを用いて所定区間内の波形信号をウェーブレット変換し、このウェーブレット変換により生成された各スケールレベル毎のウェーブレット変換信号に対し、該スケールレベルに応じた区間幅毎に区間二乗平均を算出し、前記所定区間内の任意の点において前記区間二乗平均の算出値が最大となるウェーブレット変換信号のスケールレベルをその点における支配レベルとし、この支配レベルに応じた区間幅毎に前記波形信号が最大値または最小値となる点を波形信号の特徴点として検出することを特徴とする波形信号の特徴点計測方法。
【請求項２】
請求項１記載の波形信号の特徴点計測方法において、前記ウェーブレット変換に次式（１）、すなわち、
【数１】
2004164346.tif 000019
ここで、ｆ（ｘ）は波形信号、ψ（ｘ）はマザーウェーブレット、ｂｊはスケールパラメータであり、ｂは定数（ただし、ｂ＞１）、ｊは０または負の整数からなるスケールレベル、ｋはトランスレートパラメータである
を用いることを特徴とする波形信号の特徴点計測方法。
【請求項３】
請求項２記載の波形信号の特徴点計測方法において、前記マザーウェーブレットが次式（２）、すなわち、
【数２】
2004164346.tif 000020
で定義されるフレンチハット型の関数であることを特徴とする波形信号の特徴点計測方法。
【請求項４】
請求項２記載の波形信号の特徴点計測方法において、前記マザーウェーブレットが次式（３）、すなわち、
【数３】
2004164346.tif 000021
で定義されるメキシカンハット型の関数であることを特徴とする波形信号の特徴点計測方法。
【請求項５】
請求項２記載の波形信号の特徴点計測方法において、前記
区間二乗平均の演算に、次式（４）、すなわち、
【数４】
2004164346.tif 000022
ここで、ｊは式（１）で用いたスケールレベル、ｋはトランスレートパラメータ、ｐｊはスケールレベルｊが低レベルになるほど大きくなるようにスケールレベルｊに応じて設定される定数である
を用いることを特徴とする波形信号の特徴点計測方法。
【請求項６】
請求項５記載の波形信号の特徴点計測方法において、前記
区間二乗平均の演算式（４）におけるｐｊが次式（５）、すなわち、
【数５】
2004164346.tif 000023
ここで、ａは式（１）で用いたマザーウエーブレットψ（ｘ）のサポートによって決定される定数、ｂは式（１）で用いた定数、ｊは式（１）で用いたスケールレベルである
で定義されるものであることを特徴とする波形信号の特徴点計測方法。
【請求項７】
請求項２または５記載の波形信号の特徴点計測方法において、前記式（１）におけるｂの値を２とすることを特徴とする波形信号の特徴点計測方法。
【請求項８】
請求項１記載の波形信号の特徴点計測方法において、前記波形信号は、木材年輪画像等の対象画像から該画像上に設定した計測線に沿って取得した画素上の輝度または濃度情報信号であることを特徴とする波形信号の特徴点計測方法。
【請求項９】
請求項１記載の波形信号の特徴点計測方法において、前記波形信号は、木材年輪画像等の対象画像から該画像上に設定した計測線に沿って取得した画素上の輝度または濃度情報信号をさらに微分処理したものであることを特徴とする波形信号の特徴点計測方法。
【請求項１０】
請求項９記載の波形信号の特徴点計測方法において、前記微分処理は複数画素分の間隔を隔てた複数の画素同士の差分演算処理であることを特徴とする波形信号の特徴点計測方法。
【請求項１１】
請求項８記載の波形信号の特徴点計測方法において、前記波形信号が濃度情報信号であるとき、該濃度情報信号をｆ（ｘ）とし、前記支配レベルをｊｄとし、該支配レベルに応じた定数をｑｊｄとし、前記計測線上の任意の点をｘとしたとき、ｆ（ｘ）の値が区間［ｘ−ｑｊｄ，ｘ＋ｑｊｄ］におけるｆ（ｘ）の最大値ｆｍａｘ（ｘ）と等しければ、その点ｘを、年輪層内の最大濃度箇所を示す特徴点と決定することを特徴とする波形信号の特徴点計測方法。
【請求項１２】
請求項９記載の波形信号の特徴点計測方法において、前記波形信号が濃度情報信号を微分した微分信号であるとき、該微分信号をｆ（ｘ）とし、前記支配レベルをｊｄとし、該支配レベルに応じた定数をｑｊｄとし、前記計測線上の任意の点をｘとしたとき、ｆ（ｘ）の値が区間［ｘ−ｑｊｄ，ｘ＋ｑｊｄ］におけるｆ（ｘ）の最小値ｆｍｉｎ（ｘ）と等しければ、その点ｘを、年輪層内の晩材終端を示す特徴点と決定することを特徴とする波形信号の特徴点計測方法。
【請求項１３】
請求項８または９記載の波形信号の特徴点計測方法において、前記計測線が主計測線と、該主計測線の両側に各々等間隔を隔てて並列に設定される複数の副計測線とからなり、該主計測線および副計測線の始端から各々略同一距離の点において波形信号の特徴点が検出されたとき、少なくとも該特徴点の数が主計測線および副計測線の本数に対し過半数であることを条件の１つとして、その特徴点を主計測線上の特徴点として決定することを特徴とする波形信号の特徴点計測方法。
【請求項１４】
請求項１３記載の波形信号の特徴点計測方法において、前記主計測線の両側に各々２本の副計測線を設定し、該主計測線および副計測線の始端から各々略同一距離の点において波形信号の特徴点が検出されたとき、少なくとも主計測線上と該主計測線に近い位置にある２本の副計測線のうちのいずれか一方上に特徴点があり且つ他のいずれかの副計測線上に特徴点がある場合、または、主計測線に近い位置にある２本の副計測線上と他のいずれかの副計測線上に特徴点がある場合に、その特徴点を主計測線上の特徴点として決定することを特徴とする波形信号の特徴点計測方法。
【請求項１５】
請求項８または９記載の波形信号の特徴点計測方法において、対象画像から該画像上に設定した計測線に沿って取得した画素上の輝度または濃度情報に対し周辺画素の情報を用いて平滑化することを特徴とする波形信号の特徴点計測方法。
【請求項１６】
波形上の特徴点の値または特徴点間距離が不規則な波形信号の特徴点を検出する装置であって、所定のマザーウェーブレットと複数段階のスケールレベルを用いて所定区間内の波形信号をウェーブレット変換するウェーブレット変換手段と、このウェーブレット変換手段により生成された各スケールレベル毎のウェーブレット変換信号に対し、該スケールレベルに応じた区間幅毎に区間二乗平均を算出する区間二乗平均演算手段と、前記所定区間内の任意の点において前記区間二乗平均の算出値が最大となるウェーブレット変換信号のスケールレベルをその点における支配レベルとして決定する支配レベル決定演算手段と、この支配レベルに応じた区間幅毎に前記波形信号が最大値または最小値となる点を波形信号の特徴点として検出する特徴点検索手段と、を有することを特徴とする波形信号の特徴点計測装置。
【請求項１７】
請求項１５記載の波形信号の特徴点計測装置において、検出された波形信号の特徴点に基づき特徴点間距離を算出する特徴点間距離手段をさらに有することを特徴とする波形信号の特徴点計測装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、波形上の特徴点の値または特徴点間距離が不規則な波形信号の特徴点を精度良く検出する方法および装置に関し、更に詳しくは、例えば木材年輪数、年輪幅等の計測に好適に用いることができる波形信号の特徴点計測方法および装置に関する。
【０００２】
【従来の技術】
年輪年代学では、年代確定の基準となる標準的な年輪幅変動のデータベースと照合することにより、各年輪層が何年に形成されたものかを一年単位で確定することができる。このデータベースは、奈良文化財研究所による精力的な調査研究によって、現在のところ日本では檜材で紀元前９１２年、杉材で紀元前１３１３年まで遡ることができる。ちなみに、年輪年代学の先進地域であるドイツでは、約一万年にわたる標準データベースが作成されている。年輪年代学では、主に以下のような事例が取り扱われている。
【０００３】
（ａ）木材の伐採年の推定
（ｂ）木造文化財（建築物、仏像、美術工芸品等）の制作年代や修理経緯の推定、真贋鑑定等
（ｃ）過去の長期にわたる気候変動の研究、地球環境温暖化の研究等
年輪年代学に使用される各年輪幅の時系列データは、誤検出（年輪以外のものを年輪として認識する誤り）も、非検出（年輪を認識しない誤り）も共に０という究極の検出性能が必要とされる。
【０００４】
このため、年輪幅の計測は、計測顕微鏡を用いた専用システムによって目視で行われており、その作業には熟練した技術と膨大な時間（３００層クラスの試料で約１時間）を要していた。また、システム構成が大掛かりなものになるという問題もかかえていた。
【０００５】
計測作業を自動化するためには、デジタルカメラやスキャナ等の画像取得装置を用いて取得した年輪画像から、パソコンレベルのコンピュータを用いて各年輪幅の計測を行う方法が考えられており、これまでにもいくつかの試みがなされているが、前述の検出性能やシステム規模、価格などの点で問題があり、研究手段として一般的に広く用いられるには至っていないのが現状である。
【０００６】
とりわけ檜材については、文化財等に広く利用されていて年輪年代学上重要な樹種であるにもかかわらず、年輪幅が狭小であること、杉材に比べて年輪が不明瞭なことなどの条件から、自動計測の実用化が極めて困難であった。
【０００７】
本特許案件に類似した公知の技術の主要なもの以下に示す。
【０００８】
（１）「Ｗｉｎ  ＤＥＮＤＲＯ」  カナダ  Ｒｅｇｅｎｔ  Ｉｎｓｔｒｕｍｅｎｔｓ社  １９８８年
（ｈｔｔｐ：／／ｗｗｗ．ｒｅｇｅｎｔ．ｑｃ．ｃａ．ｐｒｏｄｕｃｔｓ／ｄｅｎｄｒｏ／ＤＥＮＤＲＯ．ｈｔｍｌ参照）
ケベック大学のＤｒ．  Ｒｅｊｅａｎ  Ｇａｇｎｏｎ  と  Ｄｒ．  Ｈｕｂｅｒｔ  Ｍｏｒｉｎ  が設計し、Ｒｅｇｅｎｔ  Ｉｎｓｔｒｕｍｅｎｔｓ社が商品化した年輪年代学用ソフトである。このソフトは、年輪画像の輝度変化情報をもとに、年輪計測や木材組織の解析を行うことができる。このソフトのアルゴリズムの詳細は不明であるが、カタログ文面から察するところ、ウェーブレット処理および複数計測線情報の統合技術は用いていないと思われる。
【０００９】
（２）「画像処理システムを用いた年輪幅計測」  野田真人  １９９０年  樹木年輪研究会にて発表
この発表では、杉年輪の計測は可能であったが、檜年輪は計測不可能と結論づけている。この発表に係る計測方法においては、ウェーブレット処理と複数計測線情報の統合技術は用いられていない。
【００１０】
（３）特開平１１−２３２４２７号公報
年輪数測定のために画像輝度情報を用いる点が記載されているが、同技術は上記（２）によりすでに公知である。この公報記載の公知技術においても、ウェーブレット処理と複数計測線情報の統合技術は何ら用いられていない。
【００１１】
【発明が解決しようとする課題】
年輪画像から画像上の計測線に沿って画素情報を取得することにより、輝度変化情報の波形信号または濃度変化情報の波形信号を得ることができる。濃度波形の極大箇所（輝度波形の場合は極小箇所）は各年輪層の各年輪層のうちの最も濃い部分（晩材最高濃度部）に対応する。したがって、濃度信号波形の極大箇所または輝度信号波形の極小箇所を認識していくことにより、各年輪層を識別することができる。
【００１２】
一方、濃度波形信号をさらに微分処理すると、暗部から明部へと移行する点（微分波形の極小箇所）が各年輪層の終端点（晩材終端）に対応することなる。したがって、計測線上でこの微分波形信号の極小箇所間距離を計測することにより、年輪幅をより正確に計測することができる。
【００１３】
年輪画像から画像上の計測線に沿って画素情報を取得することにより、輝度変化情報の波形信号または濃度変化情報の波形信号等を得ることは、上記した公知文献等により公知である。
【００１４】
しかしながら、木材の年輪幅は一般に不規則であり、最大年輪幅が最小年輪幅の数倍から約１００倍程度まで異なることも少なくない、したがって、年輪画像から得た波形信号に対する特徴点（極大箇所であるピーク点または極小箇所であるトラフ点）の検出精度を小さい特徴点間距離のレベルに合わせると、特徴点間距離の大きい領域では年輪に無関係なノイズを拾いやすくなる。一方、波形信号に対する特徴点（ピーク点またはトラフ点）の検出精度を大きい特徴点間距離のレベルに合わせると、特徴点間距離の小さい領域では特徴点の検出漏れが生じやすくなる。したがって、最小年輪幅と最大年輪幅との差が大きい計測対象物に対し年輪数や年輪幅を正確に計測することが困難となっている。
【００１５】
さらに、年輪画像においては、濃度のレベルが一定ではないため、年輪画像から得られる上記輝度波形信号、濃度波形信号、微分波形信号等はいずれも計測区間全体にわたりうねったような変化特性を有する場合が多い。このため、一定の閾値を用いて波形信号における特徴点（極大箇所であるピーク点または極小箇所であるトラフ点）を検出しようとすると、検出漏れが生じる原因となり、年輪数や年輪幅を正確に計測できない原因となる。
【００１６】
したがって、年輪年代学にとって最も基礎的なデータとなる各年輪幅の時系列データを迅速かつ高精度に得ることができる計測方法および装置が要望されている。
【００１７】
本発明は、上記の要望に応え得る波形信号の特徴点計測方法および装置を提供することを目的とするものであり、波形の特徴点の値または特徴点間距離が不規則な波形信号であっても、波形の特徴点を迅速かつ高精度に検出することができる計測方法および装置を提供しようとするものである。
【００１８】
【課題を解決するための手段】
上記課題を解決するために、本発明は、第１発明として、波形上の特徴点の値または特徴点間距離が不規則な波形信号の特徴点を検出する方法であって、所定のマザーウェーブレットと複数段階のスケールレベルを用いて所定区間内の波形信号をウェーブレット変換し、このウェーブレット変換により生成された各スケールレベル毎のウェーブレット変換信号に対し、該スケールレベルに応じた区間幅毎に区間二乗平均を算出し、前記所定区間内の任意の点において前記区間二乗平均の算出値が最大となるウェーブレット変換信号のスケールレベルをその点における支配レベルとし、この支配レベルに応じた区間幅毎に前記波形信号が最大値または最小値となる点を波形信号の特徴点として検出することを特徴とする波形信号の特徴点計測方法を提供する。
【００１９】
さらに、第２発明として、上記第１発明の構成を有する計測方法において、前記ウェーブレット変換に次式（６）、すなわち、
【００２０】
【数６】
2004164346.tif 000003
ここで、ｆ（ｘ）は波形信号、ψ（ｘ）はマザーウェーブレット、ｂｊはスケールパラメータであり、ｂは定数（ただし、ｂ＞１）、ｊは０または負の整数からなるスケールレベル、ｋはトランスレートパラメータである
を用いた波形信号の特徴点計測方法を提供する。
【００２１】
さらに、第３発明として、上記第２発明の構成を有する計測方法において、前記マザーウェーブレットが次式（７）、すなわち、
【００２２】
【数７】
2004164346.tif 000004
で定義されるフレンチハット型の関数を用いることを特徴とする波形信号の特徴点計測方法を提供する。
【００２３】
さらに、第４発明として、上記第２発明の構成を有する計測方法において、前記マザーウェーブレットが次式（８）、すなわち、
【００２４】
【数８】
2004164346.tif 000005
で定義されるメキシカンハット型の関数であることを特徴とする波形信号の特徴点計測方法を提供する。
【００２５】
さらに、第５発明として、上記第２発明の構成を有する計測方法において、前記区間二乗平均の演算に、次式（９）、すなわち、
【００２６】
【数９】
2004164346.tif 000006
ここで、ｊは式（６）で用いたスケールレベル、ｋはトランスレートパラメータ、ｐｊはスケールレベルｊが低レベルになるほど大きくなるようにスケールレベルｊに応じて設定される定数である
を用いることを特徴とする波形信号の特徴点計測方法を提供する。
【００２７】
さらに、第６発明として、上記第５発明の構成を有する計測方法において、前記区間二乗平均の演算式（９）におけるｐｊが次式（１０）、すなわち、
【００２８】
【数１０】
2004164346.tif 000007
ここで、ａは式（６）で用いたマザーウエーブレットψ（ｘ）のサポートによって決定される定数、ｂは式（６）で用いた定数、ｊは式（６）で用いたスケールレベルである
で定義されるものであることを特徴とする波形信号の特徴点計測方法を提供する。
【００２９】
さらに、第７発明として、上記第２発明または第５発明の構成を有する計測方法において、前記式（６）におけるｂの値を２とすることを特徴とする波形信号の特徴点計測方法を提供する。
【００３０】
さらに、第８発明として、上記第１発明の構成を有する計測方法において、前記波形信号は、木材年輪画像等の対象画像から該画像上に設定した計測線に沿って取得した画素上の輝度または濃度情報信号であることを特徴とする波形信号の特徴点計測方法を提供する。
【００３１】
さらに、第９発明として、上記第１発明の構成を有する計測方法において、前記波形信号は、木材年輪画像等の対象画像から該画像上に設定した計測線に沿って取得した画素上の輝度または濃度情報信号をさらに微分処理したものであることを特徴とする波形信号の特徴点計測方法を提供する。
【００３２】
さらに、第１０発明として、上記第９発明の構成を有する計測方法において、前記微分処理は複数画素分の間隔を隔てた複数の画素同士の差分演算処理であることを特徴とする波形信号の特徴点計測方法を提供する。
【００３３】
さらに、第１１発明として、上記第８発明の構成を有する計測方法において、前記波形信号が濃度情報信号であるとき、該濃度情報信号をｆ（ｘ）とし、前記支配レベルをｊｄとし、該支配レベルに応じた定数をｑｊｄとし、前記計測線上の任意の点をｘとしたとき、ｆ（ｘ）の値が区間［ｘ−ｑｊｄ，ｘ＋ｑｊｄ］におけるｆ（ｘ）の最大値ｆｍａｘ（ｘ）と等しければ、その点ｘを、年輪層内の最大濃度箇所を示す特徴点と決定することを特徴とする波形信号の特徴点計測方法を提供する。
【００３４】
さらに、第１２発明として、上記第９発明の構成を有する計測方法において、前記波形信号が濃度情報信号を微分した微分信号であるとき、該微分信号をｆ（ｘ）とし、前記支配レベルをｊｄとし、該支配レベルに応じた定数をｑｊｄとし、前記計測線上の任意の点をｘとしたとき、ｆ（ｘ）の値が区間［ｘ−ｑｊｄ，ｘ＋ｑｊｄ］におけるｆ（ｘ）の最小値ｆｍｉｎ（ｘ）と等しければ、その点ｘを、年輪層内の晩材終端を示す特徴点と決定することを特徴とする波形信号の特徴点計測方法を提供する。
【００３５】
さらに、第１３発明として、上記第８発明または第９発明の構成を有する計測方法において、前記計測線が主計測線と、該主計測線の両側に各々等間隔を隔てて並列に設定される複数の副計測線とからなり、該主計測線および副計測線の始端から各々同一距離の点において波形信号の特徴点が検出されたとき、少なくとも該特徴点の数が主計測線および副計測線の本数に対し過半数であることを条件の１つとして、その特徴点を主計測線上の特徴点として決定することを特徴とする波形信号の特徴点計測方法を提供する。
【００３６】
さらに、第１４発明として、上記第１３発明の構成を有する計測方法において、前記主計測線の両側に各々２本の副計測線を設定し、該主計測線および副計測線の始端から各々略同一距離の点において波形信号の特徴点が検出されたとき、少なくとも主計測線上と該主計測線に近い位置にある２本の副計測線のうちのいずれか一方上に特徴点があり且つ他のいずれかの副計測線上に特徴点がある場合、または、主計測線に近い位置にある２本の副計測線上と他のいずれかの副計測線上に特徴点がある場合に、その特徴点を主計測線上の特徴点として決定することを特徴とする波形信号の特徴点計測方法を提供する。
【００３７】
さらに、第１５発明として、上記第８発明または第９発明の構成を有する計測方法において、対象画像から該画像上に設定した計測線に沿って取得した画素上の輝度または濃度情報に対し周辺画素の情報を用いて平滑化することを特徴とする波形信号の特徴点計測方法を提供する。
【００３８】
さらに、第１６発明として、波形上の特徴点の値または特徴点間距離が不規則な波形信号の特徴点を検出する装置であって、所定のマザーウェーブレットと複数段階のスケールレベルを用いて所定区間内の波形信号をウェーブレット変換するウェーブレット変換手段と、このウェーブレット変換手段により生成された各スケールレベル毎のウェーブレット変換信号に対し、該スケールレベルに応じた区間幅毎に区間二乗平均を算出する区間二乗平均演算手段と、前記所定区間内の任意の点において前記区間二乗平均の算出値が最大となるウェーブレット変換信号のスケールレベルをその点における支配レベルとして決定する支配レベル決定演算手段と、この支配レベルに応じた区間幅毎に前記波形信号が最大値または最小値となる点を波形信号の特徴点として検出する特徴点検索手段と、を有することを特徴とする波形信号の特徴点計測装置を提供する。
【００３９】
さらに、第１７発明として、上記第１５発明の構成を有する計測装置において、検出された波形信号の特徴点に基づき特徴点間距離を算出する特徴点間距離手段をさらに有することを特徴とする波形信号の特徴点計測装置を提供する。
【００４０】
【発明の実施の形態】
以下、図面を参照して本発明の好ましい実施形態を説明する。
【００４１】
図１および図２は、本発明に係る波形信号の特徴点計測方法を年輪計測に適用した場合の一実施形態を示すフローチャートである。図１および図２を参照すると、年輪計測の手順は次のプロセスからなる。
【００４２】
１．年輪画像の取得（ステップＳ１）
２．計測箇所の指定（ステップＳ２）
３．画素情報の取得（ステップＳ３）
４．濃度変換（ステップＳ４、Ｓ５）
５．平滑化処理（ステップＳ６、Ｓ７）
６．ピーク信号またはエッジ信号の取得（ステップＳ８〜Ｓ１２）
７．マザーウェーブレットの指定（ステップＳ１３）
８．ウェーブレット変換（ステップＳ１４）
９．区間二乗平均の演算（ステップＳ１５、Ｓ１６）
１０．支配レベルの決定（ステップＳ１７）
１１．区間内最大値（エッジ信号の場合は区間内最小値）の決定（ステップＳ１８〜Ｓ２０）
１２．主計測線と副計測線の相互参照（ステップＳ２１）
１３．副計測線間の相互参照（ステップＳ２２）
１４．年輪箇所の決定（ステップＳ２３）
１５．年輪幅の計測および結果出力（ステップＳ２４、Ｓ２５）
また、本発明の一実施形態に係る波形信号の特徴点計測装置は、上記ステップＳ２からステップＳ２５までのプロセスを実行することができるプログラムであり、基本システムがインストールされたパソコンなどのコンピュータ装置にインストールすることにより、実行可能となる。
【００４３】
この実施形態におけるプログラムは、計測箇所の指定手段、画素情報取得手段、濃度変換手段、平滑化処理手段、ピーク信号またはエッジ信号取得手段、マザーウェーブレット指定手段、ウェーブレット変換手段、区間二乗平均の演算手段、支配レベル決定手段、区間内最大値（エッジ信号の場合は区間内最小値）の決定手段、主計測線と副計測線の相互参照手段、副計測線間の相互参照手段、年輪箇所の決定手段、年輪幅の計測手段および結果出力手段を有する。
【００４４】
次に、上記各プロセスについて詳述する。
【００４５】
１．年輪画像の取得（ステップＳ１）
デジタルカメラやスキャナ等の撮像装置を用いて年輪画像を取得する。例えば小型の移動可能な試料の場合にはスキャナを使用し、大型の試料には高解像度タイプのデジタルカメラを使用する。この際、最も注意すべき点は、所定の画像解像度を得ることである。檜材では、最も年輪間隔の狭小な部分で１層あたり約０．１ｍｍという箇所を生ずることがある。このような部分を識別するためには、ナイキストの標本化定理により、１画素あたりの大きさをその半分の０．０５ｍｍ以下にすべきである。
【００４６】
さらに安全と計測誤差を考慮すると、１画素あたりの大きさ約０．０２ｍｍに相当する１２００ｄｐｉの解像度で画像取得を行うことが好ましい。大面積の試料でこの解像度を計測面全域にわたって保持するためには、一回あたりの撮影を小領域に限定し、画像を継ぎ合わせるという手法を用いることが好ましい。
【００４７】
本発明で扱う年輪画像としては、上記のようなデジタル画像取得装置を用いて試料表面から直接取得した画像のほかに、銀塩写真（通常のフィルムを使用する写真）等をデジタル変換した画像も同様に扱うことができる。また、Ｘ線写真、Ｘ線ＣＴ画像（Ｘ線断層写真）、ＭＲＩ画像のような試料内部の年輪画像も扱うことができる。
【００４８】
２．計測箇所の指定（ステップＳ２）
年輪幅計測は、図３に示すように、少なくとも１本の計測線１に沿って行われる。木口（横断面）円盤上では、通常、計測線は樹木中心（図３中左方向）から樹皮方向（図３中右方向）に向かって設定される。また、柾目（縦断面）では、木理と直交するように樹木中心方向から樹皮方向に向かって設定される。本発明においては、好ましくは、認識性能の向上を図るため、１本の主計測線１に対し、それと平行に適切な間隔を置いて主計測線１の両側に各々複数の副計測線２、３、４、５を設定する（図３参照）。
【００４９】
副計測線の数は、木口で２本、柾目面で４本が適当であることが我々の実験によって確かめられている。ただし、試料による個体差などもあり、必ずしもこの限りではない。実際のプログラム上では、副計測線を設定する間隔と本数、および年輪画像上で主計測線１の開始点１ａと終端点１ｂの２点を指定することにより、計測箇所を指定可能とすることが好ましい。
【００５０】
３．画素情報の取得（ステップＳ３）
上記したように、主計測線１および副計測線２〜５に沿って各画像の画素情報（通常のパソコンの場合には、ＢＧＲ３チャネルの８ビットデジタル信号）を取得する。画像取得に関しては、計測線上の画素情報のみでなく、後述の平滑化処理のために計測線周辺画素（具体的な周辺画素数は、平滑化の程度によって決定される）の情報も同時に取得することが好ましい。
【００５１】
なお、本発明方法による年輪認識では色情報を使用しないため、ＢＧＲ信号の平均信号や混合信号、あるいは特定チャネルの信号のみ（たとえばＧ信号のみ）の輝度に基づく信号を使用することができる。通常は、Ｂ：Ｇ：Ｒ＝１：１：１の割合で平均混合した信号を使用すればよい。
【００５２】
４．濃度変換（ステップＳ５）
上記方法で取得された画素情報は、各画素の輝度に起因しているため、値が大きくなるほど明るくなるという特性を有する。さらに、パソコンシステム上のＲＧＢ信号には、通常、ＣＲＴモニタの特性を補正するためにガンマ補正という非線形変換が施されている。一方、年輪の濃淡変化は年輪の細胞密度に起因しているため、密度あるいは密度に関連している画像濃度をもとに濃淡を記述するのが、年輪年代学の慣例である。
【００５３】
濃度記述では、値が大きくなるほど暗くなる。本発明による実施形態では、こうした年輪年代学の慣例に従い、ＲＧＢ平均信号を濃度信号に変換している（図４参照）。変換方法としては、画像中の濃度既知部分（通常はステップタブレットチャート）から導出される変換式を用いることができる。変換式の一例をあげると、
【００５４】
【数１１】
2004164346.tif 000008
また、非線形の変換ではなく、下記のような線形変換で代用してもよい。
【００５５】
【数１２】
2004164346.tif 000009
なお、２値化等の画像処理を施す場合や、従来の方法では年輪識別のための閾値設定の必要性から、年輪信号波形の振幅特性をどのように記述するかが重要であった。
【００５６】
しかし、本発明では後述するとおり、年輪信号波形の値そのものとは無関係に年輪識別を行うことができるため、ＲＧＢ、輝度あるいは濃度いずれの信号を用いても、年輪識別性能に大差はない。したがって、処理ステップや時間短縮を図るためには、このプロセスは省略可能である。ただし、このプロセスを省略すると年輪信号波形は反転することに留意する必要がある。
【００５７】
５．平滑化処理（ステップＳ６、Ｓ７）
上記プロセスによって取得した主計測線１および副計測線２〜５上の各画素値に対し、周辺画素の情報を用いて平滑化を行う。平滑化を施すことによって、試料自体や測定機器に起因するノイズによる誤認識を軽減することができる（図５参照）。
【００５８】
平滑化方法には、移動平均処理または、移動メジアン処理を用いることができる。周辺画素数としては、例えば１２００ｄｐｉ解像度の場合、我々の実験では木口面で、計測線直交方向に５画素、計測線方向に３画素程度、柾目面で計測線直交方向に１５画素、計測線方向に５画素程度が適切であったが、これらの画素数に限定されない。なお、平滑化処理を行わない場合には、平滑化処理範囲を計測線直交方向および計測線方向に各々１画素として設定する。
【００５９】
６．ピーク信号またはエッジ信号の取得（ステップＳ８〜Ｓ１２）
上記プロセスによって取得した信号は、波形の極大箇所が各年輪層の最も濃い部分（晩材最高濃度部）を表している。この濃度情報信号を、以下、ピーク信号と称することにする。このピーク信号の特徴点であるピーク点（波形の極大箇所）を認識していくことにより、各年輪層を識別することができる。
【００６０】
一方、このピーク信号の微分信号（実際の演算は差分信号）における特徴点であるトラフ点（波形の極小箇所）は、年輪画像の暗部（高濃度部）から明部（低濃度部）へと移行する点に対応し、各年輪層の終端点（晩材終端）に対応する。この微分信号を、以下、エッジ信号と称することにする。
【００６１】
通常の微分信号は、隣接画素どうしの差分演算によって得られるが、隣接差分演算ではノイズの影響を受けやすくなる。そこで、本発明の好ましい実施形態では、複数画素の間隔を隔てた複数の画素同士の差分演算を行うことにより、ノイズの軽減を図っている（図６中央の波形説明参照）。
【００６２】
ピーク信号（図６上部の波形）とエッジ信号（図６下部の波形）を比較すると、波形ノイズの影響などにより、各年輪の識別にはピーク信号の方が検出性能は良好であるが、年輪幅を正確に計測するという観点では晩材終端に対応しているエッジ信号の方が精度は高い。そこで、本発明の好ましい実施形態では、このような特徴を活かして、ピーク信号による認識、エッジ信号による認識、両者を併用した認識をメニュー上で選択可能にしている。
【００６３】
後述の「８．ウェーブレット変換」の箇所で述べるウェーブレット変換（畳み込み積分演算）では、ある１点の画素に対して周辺画素を含む演算を行う。したがって、計測線両端部（図３の開始点１ａおよび終端点１ｂ）では、画素不足による演算不能を防止するため、信号自身を折り返してダミーの信号（図７の波形の両端の点線部分参照）を作成しておくことが望ましい。ダミー信号として必要な長さは、後述の「７．マザーウェーブレットの指定」で指定するマザーウェーブレットのサポートと「８．ウェーブレット変換」で指定するレベル数に依存するが、最も長い端部所理の長さは最低位のレベルｊの時であるから、これに合わせてダミー信号を延長しておけば十分である。
【００６４】
７．マザーウェーブレットの指定（ステップＳ１３）
マザーウェーブレットとして用いることができる関数ψ（ｘ）は、
（ａ）端部が０または０に収束する関数であること
（ｂ）全区間での総和（積分値）が０であること
の２つの条件を満たす関数である。
【００６５】
また次の条件は、必ずしも満たす必要はないが、可能であれば満たしていることが望ましい。
【００６６】
（ｃ）サポート（関数の値が０ではない区間）がコンパクトに纏まっていること
マザーウェーブレットの具体例として、次のフレンチハット型の関数（式（１３））とメキシカンハット型の関数（式（１４））をあげることができる（図８（ａ）、（ｂ）参照）。
【００６７】
【数１３】
2004164346.tif 000010
【００６８】
【数１４】
2004164346.tif 000011
本発明の好ましい実施形態においては、生成が容易でサポートもコンパクトにまとまっているため、式（１３）のフレンチハット型関数をデフォルトとして用いている。もちろん、指定すればメキシカンハット型関数やその他のマザーウェーブレットも使用可能である。
【００６９】
８．ウェーブレット変換（ステップＳ１４）
１次元の画像信号ｆ（ｘ）に対するウェーブレット変換を次式（１５）で定義する。
【００７０】
【数１５】
2004164346.tif 000012
式（１４）において、ｆ（ｘ）は計測線上のピーク信号またはエッジ信号である。また、上記「６．ピーク信号またはエッジ信号の取得」で述べたように、信号両端部は、信号自身を折り返すことによって適切に処理されている。ψ（ｘ）は「７．マザーウェーブレットの指定」で述べたマザーウェーブレットである。２ｊはスケールパラメータで、レベル（スケールレベル）ｊは０または負の整数である。このスケールパラメータによって、マザーウエーブレットの幅および高さを同一面積に保持したまま２のべき乗則にしたがって変化させることができる（図９参照）。ｋはトランスレートパラメータを表している。なお、式（１５）ではスケーリングパラメータを２ｊとしているが、これに限られず、ｂｊ（ただし、ｂ＞１）としておいて次式（１６）のように表すこともできる。
【００７１】
【数１６】
2004164346.tif 000013
式（１６）では、１＜ｂ＜２のとき、スケーリングパラメータのステップ間隔は式（１４）の時よりも細かくなり、２＜ｂのとき、スケーリングパラメータのステップ間隔は式（１５）の時よりも粗くなる。
【００７２】
我々の実験で確認されているように、通常は式（１５）で表される２のべき乗ステップで良好な検出性能を得ている。
【００７３】
式（１５）や式（１６）の意味するところは、大きさ（スケーリングパラメータによって規定される）と位置（トランスレートパラメータによって規定される）を変化させたマザーウエーブレットψ（ｘ）と画像から形成される信号ｆ（ｘ）の畳み込み積分演算である。この演算によって生成されるｄｊ（ｘ）は、ｆ（ｘ）の局所的な周期とマザーウェーブレットψ（ｘ）のサポート（ψが０でない区間）が概ね一致したときに激しく波立つという特徴がある（図１０参照）。
【００７４】
図９や図１０では、説明および理解を容易にするためにレベルｊを３段階にしか変化させていないが、理論的にはレベルｊを８段階に変化させると、
【００７５】
【数１７】
2004164346.tif 000014
となり、１００倍以上の年輪幅に対応することができる。我々の実験でも、８段階に設定したレベルｊによって、年輪幅０．１ｍｍから１ｃｍを超える箇所まで対応できることが確認された。したがって、レベルｊの段階数は、計測対象物の特性に応じて任意に設定することができる。
【００７６】
９．区間二乗平均の演算（ステップＳ１６）
「８．ウェーブレット変換」で生成されたｄｊ（ｘ）に対して、次式（１８）に示すように、区間二乗平均を演算する（図１１参照）。
【００７７】
【数１８】
2004164346.tif 000015
ただし、ａはマザーウエーブレットψ（ｘ）のサポートによって決定される定数である。
【００７８】
式（１８）では、二乗平均を求める区間は、ｘ−２−ｊａからｘ＋２−ｊａまでになっているが、これに限定されず、次式（１９）のように定義してもよい。
【００７９】
【数１９】
2004164346.tif 000016
すなわち、ｐｊはスケールレベルｊが低レベルになるほど大きくなるようにスケールレベルｊに応じて設定されるので、区間二乗平均の演算は、レベルｊの高い高周波のｄｊ（ｘ）に対しては狭い区間［ｘ−ｐｊ，  ｘ＋ｐｊ］で、レベルｊの低い低周波のｄｊ（ｘ）に対しては広い区間［ｘ−ｐｊ，  ｘ＋ｐｊ］でそれぞれ求めればよい。式（１８）に示すｘ−２−ｊａからｘ＋２−ｊａまでというのは、あくまでも標準となる積分区間である。積分区間［ｘ−ｐｊ，  ｘ＋ｐｊ］をこれよりも狭くすると、局所的な濃度変化に対して鋭敏になりやすいがノイズの影響も受けやすくなる。一方、積分区間［ｘ−ｐｊ，  ｘ＋ｐｊ］をこれよりも広くすると、局所的な濃度変化に対しては鋭敏ではなくなるものの、ノイズに対して強い安定志向になる傾向がある。したがって、積分区間［ｘ−ｐｊ，  ｘ＋ｐｊ］を上記の標準的な状態から加減してやることによって、検出の感度調節を図ることもできる。
【００８０】
１０．支配レベルの決定（ステップＳ１７）
ある点ｘにおいて、式（１８）または式（１９）で定義されたｇｊ（ｘ）が最大となるレベルｊを点ｘにおける支配レベルｊｄとする。これにより、計測線上のすべての点が、ある一つの支配レベルを持つことになる（図１２参照）。支配レベルｊｄは、年輪幅が広い部分では低くなる傾向があるのに対し、年輪幅が狭い部分や晩材周辺の信号が著しく変化する部分では高くなる傾向がある。
【００８１】
１１．区間内最大値（エッジ信号の場合は区間内最小値）の決定（ステップＳ１８〜Ｓ２０）
図１３に示すように、
（ａ）ｆ（ｘ）がピーク信号の場合には、ｆ（ｘ）の値が区間［ｘ−２−ｊｄａ，  ｘ＋２−ｊｄａ］における  ｆ（ｘ）の最大値ｆｍａｘ（ｘ）と等しければ、点ｘを年輪層内の最大濃度箇所と決定する。
【００８２】
（ｂ）ｆ（ｘ）がエッジ信号の場合には、ｆ（ｘ）の値が区間［ｘ−２−ｊｄａ，  ｘ＋２−ｊｄａ］におけるｆ（ｘ）の最小値ｆｍｉｎ（ｘ）と等しければ、点ｘ年輪層内の晩材終端と決定する。
【００８３】
これを全範囲にわたって行うことで、計測線上のすべての年輪を検出することができる。
【００８４】
なお、上記（ａ）、（ｂ）は上記したスケールパラメータｂｊにおけるｂが２の場合であるが、下記のように表すこともできる。
【００８５】
（ｃ）ｆ（ｘ）がピーク信号の場合には、ｆ（ｘ）の値が区間［ｘ−ｂ−ｊｄ，  ｘ＋ｂ−ｊｄ］におけるｆ（ｘ）の最大値ｆｍａｘ（ｘ）と等しければ、点ｘを年輪層内の最大濃度箇所と決定する。
【００８６】
（ｄ）ｆ（ｘ）がエッジ信号の場合には、ｆ（ｘ）の値が区間［ｘ−ｂ−ｊｄ，  ｘ＋ｂ−ｊｄ］におけるｆ（ｘ）の最小値ｆｍｉｎ（ｘ）と等しければ、点ｘを年輪層内の晩材終端と決定する。
【００８７】
これをさらに一般化すると、支配レベルｊｄに応じた定数ｑｊｄを用いて下記のように表すことができる。
【００８８】
（ｅ）ｆ（ｘ）がピーク信号の場合には、ｆ（ｘ）の値が区間［ｘ−ｑｊｄ，  ｘ＋ｑｊｄ］におけるｆ（ｘ）の最大値ｆｍａｘ（ｘ）と等しければ、点ｘを年輪層内の最大濃度箇所と決定する。
【００８９】
（ｆ）ｆ（ｘ）がエッジ信号の場合には、ｆ（ｘ）の値が区間［ｘ−ｑｊｄ，  ｘ＋ｑｊｄ］におけるｆ（ｘ）の最小値ｆｍｉｎ（ｘ）と等しければ、点ｘを年輪層内の晩材終端と決定する。
【００９０】
１２．主計測線と副計測線の相互参照（ステップＳ２１）
上記「２．計測箇所の指定」で定めた主計側線１上の年輪が不明瞭であった場合には、上記「３．画素情報の取得」から「１１．区間内最大値（エッジ信号の場合は区間内最小値）の決定」の手順で各層の最大濃度箇所や晩材終端を検出しようとしても検出されないことがある。また、年輪層の部分的な不明瞭は、ランダムに発生する場合が多い。この状態を非検出と呼ぶことにする。
【００９１】
また、主計側線１上に木材組織の不均一な箇所があったり、割れや研磨痕などの本来の年輪以外濃度変化（以下ノイズと称する）があったりした場合には、最大濃度箇所や晩材終端以外の箇所を検出してしまうことがある。この状態を誤検出と呼ぶことにする。
【００９２】
通常、年輪は木口では同心円状に、柾目では平行状に生じるのに対し、ノイズの発生パターンは不確定的である。これらの違いを利用して、主計測線上に非検出や誤検出を生じても、正しい検出結果を出力できるように、「２．計測箇所の指定」で設定した副計測線上の検出結果（「３．画素情報の取得」から「１１．区間内最大値（エッジ信号の場合は区間内最小値）の決定」のプロセスをそれぞれの副計測線２〜５上に対して行ったもの）と相互参照しながら、年輪として確定してもよいか否かを検証する。その具体的な手順を以下に記す。
【００９３】
実際のプログラム上では相互参照を効率よく進めるため、あらかじめ、すべての主および副計測線上の検出点の位置と計測線番号を確定候補点としてメモリーに記憶しておく。
【００９４】
木口面では、通常年輪は同心円状であるが、微小な幅の中で考えた場合には、ほぼ平行とみなすことができる。また、柾目では平行状である。ところが、いずれも計測線方向に直交する完全な平行状態ではないため、主計測線上の検出位置と副計測線上の検出位置とで若干のズレを生じることがある。本来同一年輪層であるべき検出点が、主または副計測線設定位置の違いによって生じる誤差を吸収するために、すなわち、各計測線の始端から略同一距離の範囲内にある検出点か否かを決定するために、同一とみなしうる幅ｈを設定する。
【００９５】
次に、複数存在する主および計測線上のうち、何箇所以上に検出点を生じた場合に、その検出点を正しいものとして確定してよいかを定義する。我々の実験では、国産檜材の場合、木口では主および副計測線合計３本中２箇所以上に検出点を生じたもの、柾目では主および副計測線合計５本中３箇所以上に検出点を生じたものとした場合、最も良好な結果が得られている。より一般的には、同一とみなしうる幅±ｈ以内で主および副計測線合計ｎ本中ｍ箇所以上に検出点のあるものと、定義することができる。
【００９６】
上記規則に従って、主計測線上の全検出点に対し、幅±ｈ以内で主および副計測線合計ｎ本中ｍ以上に検出点のあるものを、正しい年輪（各層内最大濃度箇所または各層内晩材終端）として確定していく（図１４参照）。これにより、
（１）主計測線上の検出点が正しい年輪であった場合には、相互参照により正しいものとして確定される。例えば、図１５の例では、検出点ａ１〜ａ３においては主計測線上に検出点があり且つ５本の計測線中３つ以上の検出点があるので、検出点が確定する。
【００９７】
（２）主計測線上の検出点が誤検出であった場合には、相互参照により排除される。例えば、図１５の例では、検出点ａ５は５本の計測線中３箇所未満であるため、確定しないで次の「１３．副計測線間の相互参照」を行う。
【００９８】
（３）主計測線上に非検出があった場合（図１５における検出点ａ４参照）には、この過程では考慮されないので、次の「１３．副計測線間の相互参照」を行う。
【００９９】
「１３．副計測線間の相互参照」においては、図１６に示すように、検出点ａ４においては３箇所以上に検出点があるので、優先順位１を検索中に確定する。しかし、検出点ａ５においては５本の計測線中３箇所未満であるので確定しない。
【０１００】
図１７に示す段階においては、検索点ａ１〜ａ４は確定済みであり、一方、検出点ａ５は検出点として確定せずに相互参照が終了する。
【０１０１】
なお、次の「１３．副計測線間の相互参照」を効率よく行うために「１２．主計測線と副計測線の相互参照」の過程で正しい年輪として確定されたものは、確定候補点から順次消去する。
【０１０２】
１３．副計測線間の相互参照（ステップＳ２２）
複数存在する副計測線２〜５に、あらかじめ優先順位をつける。通常は、主計測線１に近いほど、高い優先順位とする。このとき、同順位のないようにしておく。
【０１０３】
まず、確定候補点として残っている検出点のなかで最も優先順位の高い副計測線上の全検出点を、「１２．主計測線と副計測線の相互参照」と同様の規則に従って正しい年輪として確定していく（図１６における検出点ａ４参照）。なお、この過程で正しい年輪として確定されたものは、確定候補点から順次消去する。
【０１０４】
以下、優先順位の高いものから順に同様のことを行う。これらの一連のプロセスは、確定候補点がなくなった時点、または最も優先順位の低い副計測線に対して行った時点で終了する（図１７参照）。なお、実際には、ｎ本中ｍ以上の規則の場合には、最低順位を含む順位下位のｍ−１本に対しては無意味であるから、この順位で打切りにしてよい。
【０１０５】
上記した用に、「１２．主計測線と副計測線の相互参照」のプロセスで主計測線上に非検出があった場合にも、「１３．副計測線間の相互参照」を行うことで正しい年輪として確定される（図１５、図１６における検出点ａ４参照）。
【０１０６】
１４．年輪箇所の決定（ステップＳ２３）
「１２．主計測線と副計測線の相互参照」と「１３．副計測線間の相互参照」の年輪確定点を統合して、年輪箇所（各層内最大濃度箇所または各層内晩材終端）とする（図１４参照）。なお、「１３．副計測線間の相互参照」の過程で確定された副計測線上の確定点は、統合の際に主計測線に射影しておく。
【０１０７】
１５．年輪幅の計測および結果の出力（ステップＳ２４、Ｓ２５）
すべての年輪確定点の情報は主計測線１上に集められているため、主計測線上の各年輪確定点の座標値から各層の画素数をカウントしていく。「１．年輪画像の取得」で設定した画像解像度によって画素１個あたりの大きさが決められるので、
【０１０８】
【数２０】
2004164346.tif 000017
によって各層の年輪幅を求めることができる。この結果を、出力することにより、一連の全プロセスが完了する。
【０１０９】
ところで、年輪年代学に有用なデータは、年輪幅だけではない。層内最大濃度（すなわち層内最大密度）、層内最小濃度（すなわち層内最小密度）、早材晩材比なども、重要なデータとなる。
【０１１０】
上記したプロセスによる波形信号の特徴点計測方法の効果を図１８および図１９に示す。図１８は、波形信号に上記濃度信号を用いた場合の最大濃度点検出結果を従来方法による検出結果と比較し表にして示したものである。また、図１９は、波形信号に上記微分信号を用いた場合の最大晩材終端点検出結果を従来方法による検出結果と比較し表にして示したものである。図１８および図１９において、本発明による新しい方法は、ウェーブレット変換のみ（複数の計測線の相互参照は省略）の場合と、複数の計測線の相互参照のみ（ウェーブレット変換は無し）の場合と、ウェーブレット変換および複数の計測線の相互参照を併用した場合とに分けている。
【０１１１】
また、図１８および図１９に示す検出率、誤検出率、非検出率は、図２０に示す、Ｓｓ、Ｓｎ、Ｎｓ、Ｎｎの関係から、次式（２１）、（２２）、（２３）のように定義したものである。
【０１１２】
【数２１】
2004164346.tif 000018
図１８および図１９からわかるように、本発明の計測方法を用いれば、いずれの場合も従来の方法に比べて検出率の増加、誤検出率および非検出率の低下が確認され、検出性能が全般的に向上することが実証された。
【０１１３】
以上、年輪幅を計測する場合について詳述したが、「１４．年輪箇所の決定」までを正確に行うことができれば、ピーク信号やエッジ信号からこれら年輪幅以外の特徴量も容易に求めることができる。
【０１１４】
また、本発明による計測方法および計測装置は、木材の年輪計測のみに用途を限定されず、例えば個人認証のための指紋、音紋、網膜等の計測、電子部品における配線パターンの計測、脳波信号等の生体信号の計測等にも適用することができる。
【０１１５】
【発明の効果】
以上の説明から明らかなように、本発明による波形信号の特徴点計測方法および計測装置によれば、特徴点間距離が不規則な波形信号であっても、ピーク波形信号の特徴点を正確に検出することができるので、木材の年輪等の検出性能を著しく向上することができる。
【０１１６】
すなわち、計測法にウェーブレット変換を導入する前は、ｆ（ｘ）上にある閾値を定めて閾値以上か閾値以下かで区別する単純二値化で年輪検出を行っていたため、ノイズや不鮮明な年輪に対して非常に弱く、また閾値が個体差に左右されやすい欠点を持っていたが、本発明によるウェーブレット変換を用いた計測法を導入することにより、その場その場その場の年輪の細かさ（あるいは粗さ）や年輪の鮮明さ（あるいは不鮮明さ）に自動的に適応しながら年輪検出を行えるようになり、これが検出性能の向上に寄与するのである。
【０１１７】
また、この検出性能は、閾値設定不要なため、個体差に左右されにくい頑健な面も持ちあわせており、この点においてもきわめて画期的かつ効果的である。
【図面の簡単な説明】
【図１】本発明を年輪計測に適用した場合の一実施形態に係る波形信号の特徴点計測方法のプロセス前半部分を示すフローチャートである。
【図２】図１のフローチャートに続くプロセス後半部分を示すフローチャートである。
【図３】主計測線及び副計測線の設定例を模式的に示す説明図である。
【図４】計測線に沿って取得したＲＧＢ信号から濃度信号へ変換する様子を示す説明図であり、（ａ）はＲＢＧ信号の波形説明図、（ｂ）はＲＢＧ信号から作成した濃度信号の波形説明図、（ｃ）はＲＧＢ信号から濃度信号へ変換に用いられるＲＧＢ−濃度変換曲線の説明図である。
【図５】（ａ）は濃度信号の変化を示す波形説明図であり、（ｂ）は平均化処理した濃度信号を示す波形説明図である。
【図６】（ａ）は濃度信号の波形説明図であり、（ｂ）は濃度信号の微分処理により作成された微分信号の波形説明図、（ｃ）は微分処理（差分処理）方法を説明する説明図である。
【図７】濃度信号の両端部に延長ダミー信号を付加した例を示す波形説明図である。
【図８】（ａ）はウェーブレット変換に用いられるフレンチハット型のマザーウェーブレットの代表例を示す説明図であり、（ｂ）はメキシカンハット型マザーウェーブレットの代表例を示す説明図である。
【図９】フレンチハット型のマザーウェーブレット波形がレベル値に応じて段階的に変化する様子を示す説明図である。
【図１０】濃度信号を３段階のレベル値を用いてウェーブレット変換することにより得られる信号波形を示す説明図である。
【図１１】図１０のウェーブレット変換により得られる信号波形をレベル値に応じて区間二乗平均処理した結果を示す説明図である。
【図１２】図１１の区間二乗平均処理結果から支配レベルを決定演算した結果を示す説明図である。
【図１３】図１２の支配レベル決定演算結果を用いて濃度信号波形上の区間内最大値および微分信号波形上の区間内最小値を検索し決定するプロセスを模式的に示す説明図であり、（ａ）は支配レベル説明図、（ｂ）は濃度信号の波形説明図、（ｃ）は微分信号の波形説明図である。
【図１４】（ａ）は濃度信号からの最大濃度点（特徴点）検出結果を示す説明図であり、（ｂ）は微分信号からの晩材終端点（特徴点）検出結果を示す説明図である。
【図１５】主計測線および副計測線の相互参照による検出点（特徴点）確定プロセスを模式的に示す説明図である。
【図１６】図１５のプロセスに続けて副計測線の相互参照による検出点（特徴点）確定プロセスを模式的に示す説明図である。
【図１７】図１６のプロセスに続く相互参照による検出点（特徴点）確定プロセスの最終結果を模式的に示す説明図である。
【図１８】濃度信号を用いた場合の本発明方法による最大濃度点（特徴点）検出結果を従来方法による検出結果と比較して示す説明図である。
【図１９】濃度信号を用いた場合の本発明方法による晩材終端点（特徴点）検出結果を従来方法による検出結果と比較して示す説明図である。
【図２０】検出率と誤検出率と非検出率の定義方法を説明するための説明図である。
【符号の説明】
１    主計測線
２〜５    副計測線
【図１】
2004164346.tif 000024
【図２】
2004164346.tif 000025
【図３】
2004164346.tif 000026
【図４】
2004164346.tif 000027
【図５】
2004164346.tif 000028
【図６】
2004164346.tif 000029
【図７】
2004164346.tif 000030
【図８】
2004164346.tif 000031
【図９】
2004164346.tif 000032
【図１０】
2004164346.tif 000033
【図１１】
2004164346.tif 000034
【図１２】
2004164346.tif 000035
【図１３】
2004164346.tif 000036
【図１４】
2004164346.tif 000037
【図１５】
2004164346.tif 000038
【図１６】
2004164346.tif 000039
【図１７】
2004164346.tif 000040
【図１８】
2004164346.tif 000041
【図１９】
2004164346.tif 000042
【図２０】
2004164346.tif 000043