(57)【特許請求の範囲】
【請求項１】
導管の流れ方向とほぼ直交する方向に変化させた異なる照明角度で撮影した２つの画像それぞれに対して微分演算を行う段階と、前記微分演算により得られる両微分データの減算を行う段階と、得られた減算データを絶対値化する段階と、絶対値化したデータに対して前記照明角度の変化方向についてスムージング処理する第１のスムージング段階と、
前記第１のスムージング段階により得られるスムージングデータの全画素を分母とした全体ヒストグラムを算出する段階と、
前記スムージングデータの特定の位置ごとに所定の画素数を分母とした部分ヒストグラムを算出する段階と、
前記全体ヒストグラムと前記部分ヒストグラムの関係に基づいて、前記スムージングデータ上の各画素の値を修正する画素修正段階と、
前記画素修正が行われたスムージングデータに対して、前記スムージング処理の処理方向と直交する方向にスムージング処理を行う第２のスムージング段階を有することを特徴とする導管画像作成方法。
【請求項２】
導管の流れ方向とほぼ直交する方向に変化させた異なる照明角度で撮影した２つの画像を入力する画像入力手段と、
入力された２つの画像それぞれに対して微分演算を行った後、微分演算された両微分データの減算を行い、得られた減算データを絶対値化する導管部強調手段と、
前記導管部強調手段により得られた導管部強調データに対して前記照明角度の変化方向についてスムージング処理する第１のスムージング処理手段と、
前記第１のスムージング処理手段により得られたスムージングデータの全画素を分母とした全体ヒストグラムを算出する全体ヒストグラム算出手段と、
前記スムージングデータの特定の位置ごとに所定の画素数を分母とした部分ヒストグラムを算出する部分ヒストグラム算出手段と、
前記全体ヒストグラム算出手段により算出された全体ヒストグラムと、前記部分ヒストグラム算出手段により算出された部分ヒストグラムの関係に基づいて、前記スムージングデータ上の各画素の値を修正する画素修正手段と、
前記修正後のスムージングデータに対して、前記第１のスムージング処理手段におけるスムージング処理と直交する方向にスムージング処理を行う第２のスムージング処理手段と、
を有することを特徴とする導管画像作成システム。
【請求項３】
コンピュータに、導管の流れ方向とほぼ直交する方向に変化させた異なる照明角度で撮影した２つの画像それぞれに対して微分演算を行う段階、前記微分演算により得られる両微分データの減算を行う段階、得られる減算データを絶対値化する段階、絶対値化したデータに対して前記照明角度の変化方向についてスムージング処理する第１のスムージング段階、前記第１のスムージング段階により得られるスムージングデータの全画素を分母とした全体ヒストグラムを算出する段階、前記スムージングデータの特定の位置ごとに所定の画素数を分母とした部分ヒストグラムを算出する段階、前記全体ヒストグラムと前記部分ヒストグラムの関係に基づいて、前記スムージングデータ上の各画素の値を修正する画素修正段階、前記画素修正が行われたスムージングデータに対して、前記スムージング処理の処理方向と直交する方向にスムージング処理を行う第２のスムージング段階、を実行させるためのプログラム。
【発明の詳細な説明】
【０００１】
【産業上の利用分野】
本発明は、木材を撮影することにより得られる木目画像から導管の部分の像のみを抽出するための技術に関する。
【０００２】
【従来の技術】
木目柄を印刷した化粧紙に対してエンボス加工を施す場合には当該木目の導管部分が凹むように加工することが望ましいことが知られている。そのため、従来より、砥の粉やトナー等の目止め剤を撮影する木材の導管に充填して撮影し、撮影して得た画像をスキャナで読み取ってデジタル化し、そのデジタル化した画像を所定の閾値によって２値化することによって導管部分のみの画像を得ている。撮影に先立って木材の導管に砥の粉等の目止め剤を充填するのは、撮影した画像において導管の部分が特に際立った明るさ、あるいは暗さになるようにして、２値化によって導管の部分のみを抽出できるようにするためである。
【０００３】
この従来からの手法では、木材の導管部分に砥の粉等の目止め剤を充填する作業が作業者にとって大きな負担となるばかりでなく、導管部分に目止め剤を充填した木材はその後何らかの用途に再利用ができないので、資源の有効利用の観点から望ましいものとはいえないものである。
【０００４】
また、撮影する木材の導管部分に目止め剤を充填して撮影したとしても、デジタル化して２値化したときに導管部分のみだけではなく、その他の部分も抽出されてしまうことがある。このような場合には、適宜な手法によって、導管部分以外の部分を削除したり、あるいは導管の形状を修正したりする、いわゆるレタッチ作業を行う必要があるが、導管の模様は非常に細かなものであるのでレタッチ作業は非常に面倒で時間が掛かるものであり、レタッチ作業を行う作業者の負担が大きいものであった。
【０００５】
そこで、本出願人は、特開平１０−３０７９０６号公報において、照明条件を変えた２枚の撮影像からコンピュータにより導管部分を強調してから２値化を行うことにより導管部分の画像を得る手法を提案した。
【０００６】
【発明が解決しようとする課題】
しかしながら、上記公報において提案した手法では、入力装置の解像度が低い場合、微細な導管をもつ木材に対して適用した場合に、充分に導管を捉えきることができず、処理結果にノイズ成分を多く含んでしまうという問題が生じる。
【０００７】
上記のような点に鑑み、本発明は、導管の細かさに対して入力装置が低解像度であるような場合であっても、ノイズの少ない導管の画像を得ることが可能な導管画像作成方法およびシステムを提供することを課題とする。
【０００８】
【課題を解決するための手段】
上記課題を解決するため、本発明では、異なる照明角度で撮影した２つの画像それぞれに対して微分演算を行う段階、前記微分演算により得られる両微分データの減算を行う段階、得られた減算データを絶対値化する段階、絶対値化したデータに対して前記照明角度の変化方向についてスムージング処理する第１のスムージング段階、前記スムージング処理により得られるスムージングデータに対して、前記スムージング処理の処理方向と直交する方向にスムージング処理を行う第２のスムージング段階、を実行するようにしたことを特徴とする。
【０００９】
本発明によれば、異なる照明角度で撮影した２枚の画像に対して、微分演算、減算処理、導管と直交する方向にスムージング処理を施して導管部を強調したスムージングデータを得た後、得られたスムージングデータに対して上記スムージングと直交（導管と平行）する方向にスムージング処理を行うようにしたので、導管の細かさに対して入力装置が低解像度であるような場合であっても、ノイズの少ない導管の画像を得ることが可能となる。
【００１０】
【発明の実施の形態】
以下、図面を参照して本発明の実施形態について詳細に説明する。
（撮影像の特性）
以下、本発明に係る導管画像作成方法について説明するが、それに先立って、通常通りに木目を正面からの照明により撮影した画像において、木目の各部がどのような明るさとなるかについて検討すると次のようである。
【００１１】
木目には導管の他に、年輪や放射組織と呼ばれる組織がある。そこで、図１（ａ）に示すような木材の断面を考える。図１（ａ）において、１は導管を示し、２の斜線を施した部分は年輪を示している。導管１は孔であるので断面をとると凹んでいるが、年輪２は平坦である。
【００１２】
そして、この木材に表面の明るさが均一になるように正面（図の矢印の方向）から照明を当て、図の矢印の方向から撮影したとすると、当該部分の画像の輝度のプロファイルは図１（ｂ）に示すようになる。即ち、導管１の部分は中央部は若干明るくなるが全体としては暗い像となり、年輪２の部分も全体的に暗い像となる。従って、このような画像に対して、導管１の部分のみを抽出しようとしてある閾値を設定して２値化し、暗い部分、即ち輝度の低い部分を抽出しても、導管１だけではなく、年輪２の部分も抽出されてしまうことになる。上述した、デジタル化して２値化したときに導管部分のみだけではなく、その他の部分も抽出されてしまうことがあるというのはこのようなことを指している。
【００１３】
ところで、図１（ｂ）に示すように導管１の部分が全体として暗い像になるのは、導管１が溝状になっているために陰影を生じているためであり、年輪２の部分が全体として暗い像となるのは、年輪２の部分の反射率がその周辺より低いからであると考えられる。
【００１４】
つまり、導管の部分と年輪の部分とは像としては共に暗い像となるが、その原因が異なっていると考えられるのである。そこで、この原因の相違を利用して導管の像と年輪の像とを分離できれば、導管の像のみを抽出できることになる。
【００１５】
そこで、従来のように木目の撮影に際して木材の正面から照明を当てるのではなく、正面から角度を有した方向から照明を当てたときに画像の輝度のプロファイルがどのようになるか考えてみる。いま、図２（ａ）の▲１▼で示す矢印の方向（以下、この方向を左方向または左側と称す）から照明を当てて撮影したとする。この画像（これを左画像と称す）の輝度のプロファイルは図２（ｂ）に示すようであり、年輪２の像は図１（ｂ）に示すと同様に全体的に暗い像となっているが、導管１の像の右側の部分は左側からの照明によって輝度が高くなっている。
【００１６】
同様にして、図２（ａ）の▲２▼で示す矢印の方向（以下、この方向を右方向または右側と称す）から照明を当てて撮影したとすると、この画像（これを右画像と称す）の輝度のプロファイルは図２（ｃ）に示すようであり、年輪２の像は図１（ｂ）に示すと同様に全体的に暗い像となっているが、導管１の像の左側の部分は右側からの照明によって輝度が高くなっている。
【００１７】
このように互いに反対方向から照明を当てて撮影した二つの画像では年輪の像の輝度のプロファイルは同様であるが、導管の像の輝度のプロファイルは異なっているのである。なお、撮影方向は従来と同様であり、通常は正面である。また、左方向、右方向は導管１の長手方向と略直交する方向にとるものとする。なお、図２（ａ）は図１（ａ）と全く同じものとなっている。
【００１８】
このことから、導管の像の輝度のプロファイルは導管の溝状の構造に基づいているのに対して、年輪の像の輝度のプロファイルは周辺の組織との反射率の相違に基づいているという仮説が正しいことが分かる。
【００１９】
このことを利用すれば、輝度のプロファイルは同様であるが、その輝度のプロファイルが生じる原因となる構造が相違する部分、具体的にはその輝度のプロファイルが表面の凹凸に起因する部分の像と、周辺部分と反射率が異なることに起因する部分の像とを分離することが可能である。そして、このことによって、木目柄の画像において導管の像と年輪の像を分離できる可能性があることが分かる。
【００２０】
（本発明に係る方法の具体的処理）
ここから、本発明に係る導管画像作成方法の具体的処理について、図６のフローチャートを用いて説明する。まず、上述したように左画像と右画像を得る（ステップＳ１）。ここで、左画像のあるラインの輝度のプロファイルは図３（ａ）のようであるとする。また、同一ラインにおける右画像の輝度のプロファイルが図４（ａ）のようであるとする。図３（ａ）、図４（ａ）において３は導管の像の輝度のプロファイルを示し、４は年輪の像の輝度のプロファイルを示している。なお、図３、図４はあくまでも輝度のプロファイルを模式的に示すものである。
【００２１】
次に、左画像と右画像の画像データに対して微分演算を施す（ステップＳ２）。画像データに微分演算を施す手法としては、例えば水平方向一次微分オペレータ等の微分オペレータを作用させても良いことは当然であるが、単純に隣接する画素間で画素値の減算をとる演算を行っても良い。図３（ａ）の画像データに微分演算を施すと図３（ｂ）に示すようになり、図４（ａ）の画像データに微分演算を施すと図４（ｂ）に示すようになる。
【００２２】
このようにして、左画像の画像データに対して微分演算を施したデータ（以下、これを左微分データと称す）と、右画像の画像データに対して微分演算を施したデータ（以下、これを右微分データと称す）が得られるが、次に、位置を合わせて左微分データと右微分データとの間で減算を行う（ステップＳ３）。以下、この減算をとって得られたデータを左右減算データと称することにする。
【００２３】
図３（ｂ）、図４（ｂ）から明らかなように、左微分データと右微分データとでは年輪の像に対応する微分データは略同一であるのに対して、導管の像に対応する微分データは異なっているので、この減算の演算によって、年輪の像に対応する微分データはほぼ除去され、導管の像に対応する微分データのみが残り、この場合の減算データは図５（ａ）に示すようになる。
【００２４】
次に、減算データを絶対値化する（ステップＳ４）。これによって、図５（ｂ）に示すデータが得られる。このデータは導管部が強調された導管部強調データとなるが、この状態で２値化すると図５（ｂ）の例では、導管の部分が３つに分割された状態で抽出されることになる。これを避けるために、導管部強調データにぼかし処理を施してスムージングする（ステップＳ５）。ここで、具体的なスムージングの手法としては、適宜のものを適用することができるが、本実施形態では、導管部強調データ上の隣接する両画素の値を利用して平均化する。具体的には、ある画素のスムージング後の画素値Ｃは、その画素のスムージング処理前の画素値をＣS、左側の画素の値をＣL、右側の画素の値ＣRとすると、以下の〔数式１〕により算出される。
【００２５】
〔数式１〕
Ｃ＝ （ＣL ＋ ２ＣS ＋ ＣR ）／４
【００２６】
すなわち、その画素の元の値に２倍の重み付けをして両隣の画素との平均を求めることになる。これによって導管の像に対応する位置には図５（ｃ）の実線で示すように一つの大きなピークが得られる。
【００２７】
なお、上述したように、左画像の画像データと右画像の画像データにそれぞれ微分演算を施して左微分データと右微分データを得、その左微分データと右微分データとの間で減算するのは次の理由による。
【００２８】
例えば、図３（ａ）に示す左画像の画像データと、図４（ａ）の右画像の画像データとの間で減算を位置を合わせて行えば、年輪の像に対応するデータを除去し、導管の像に対応するデータのみを抽出できる可能性はあるが、そのためには左画像と右画像の直流分、即ち平均輝度が同じでなければならない。
【００２９】
しかし、左画像を撮影するときの木材表面の平均輝度と、右画像を撮影するときの木材の平均輝度とは一般には同一ではない。更に、照明が近い方が遠い方より輝度が高いので、左画像では左側程輝度が高く、右画像では右側程輝度が高い。このように、左画像と右画像では直流分が異なり、しかも一つの画像の中でも位置により直流分が異なるので、単純に左画像の画像データと右画像の画像データとの減算を行っただけでは年輪の像の部分には直流分の差に対応する値が現れてくるので、２値化する際の閾値を低くすると年輪の像も抽出されてしまうことになる。閾値を高くすれば年輪の像の部分を除くことはできるが、その場合には抽出される導管の像は実際の像より細くなってしまう。
【００３０】
そこで、上述したように左右の画像の画像データに微分演算を施して直流分を除去した後に減算を行うようにしているのであり、これによれば、閾値を小さくできるので、抽出される導管の像は実際の像と略同じになる。上記ステップＳ１〜ステップＳ５の処理により撮影画像から年輪部分を除去できることになる。
【００３１】
上記ステップＳ１〜ステップＳ５までは、特開平１０−３０７９０６号公報において提案したものとほぼ同様の手順で行われる。これにより導管部分を強調し、導管と直交する方向にスムージングされたデータが得られることになる。ただ、この時点で得られるスムージングデータは、導管と直交する方向のみにスムージング処理されたものであるため、このスムージングデータを２値化して２値画像として出力すると、入力装置の解像度が低い場合には、導管の流れ方向において、抽出される導管の像が途切れてしまう現象が生じる。例えば、上記ステップＳ５までの処理において、図７（ａ）に示すようなスムージングデータが得られたとする。図７（ａ）において、各格子は画素を示しており、濃い網掛けがしてある画素ほど値が大きいことを示している。なお、このスムージングデータは、上述のように左右両微分データの減算により得られたものであるので、厳密にいえば、画像ではないが、入力された左右の画像に対応する画素を有し、各画素にスカラー値が与えられているため、演算処理上は、画像と同等に扱うことができる。
【００３２】
図７（ａ）からわかるように、スムージングデータにおいて、導管と直交する方向（図面左右方向）については、スムージング処理が行われているため、中央部の値が大きく、両側の値が小さくなっているが、導管の流れ方向（図面上下方向）については、導管に相当する部分であっても両端より値が小さくなっている部分がある。このスムージングデータを、ニ値化すると、図７（ｂ）に示すような導管画像が得られる。図７（ｂ）に示すように、上記手法で得た導管画像は、導管の一部が途切れてしまった画像となってしまう。そこで、本発明では、スムージング処理された画像に対して以下の処理が行われる。
【００３３】
まず、上記ステップＳ５の処理により得られたスムージングデータ（図７（ａ）に相当）に対して、上記ステップＳ５で行ったのと直交する方向に、スムージング処理を行う（ステップＳ６）。すなわち、上記ステップＳ５では導管の流れ方向と垂直な方向に対して行ったので、ここでは、導管の流れ方向に対して行うことになる。具体的な処理としては、ある画素の画素値Ｃは、その画素のスムージング処理前の画素値をＣS、上側の画素の値をＣU、右側の画素の値ＣDとすると、以下の〔数式２〕により算出される。
【００３４】
〔数式２〕
Ｃ＝ （ ＣU ＋ ２ＣS ＋ ＣD ）／４
【００３５】
上記〔数式２〕は、ステップＳ５において利用した〔数式１〕と同様の式であり、スムージングデータ上の隣接する両画素の値を利用して平均化を行っていることを示している。〔数式１〕では、左右の画素の値を利用したのに対して、〔数式２〕においては、上下の画素の値を利用している点が異なっている。なお、スムージングの手法としては、〔数式２〕に示したもの以外の適宜のものを適用することもできる。このようにして得られたスムージングデータに対して、あらかじめ設定された閾値により２値化を行うことにより導管部分の像を抽出した２値画像が得られる（ステップＳ７）。
【００３６】
本発明によるステップＳ６の処理を、図７（ａ）に示したスムージングデータに対して行うと、図７（ｃ）に示すようなスムージングデータが得られることになる。図７（ｃ）では、図７（ａ）に示すものに比べて中央の値が高くなっていることがわかる。このようなスムージングデータに対して２値化を行うことにより、図７（ｄ）に示すような導管画像が得られることになる。図７（ｂ）と比べると導管が途中で途切れるようなことがなくなることがわかる。以上のようにして、入力装置の解像度が導管の細かさに比べて低い場合であっても、鮮明な導管画像を得ることが可能となる。
【００３７】
（応用形態）
続いて、本発明の応用形態について以下に説明する。以下に説明する処理は、撮影時の光学系の収差による影響の除去を行うために効果的な処理であり、本発明と組み合わせることにより、より鮮明な導管画像を得ることが可能となるものである。
【００３８】
図８に本発明の導管画像作成方法の応用例のフローチャートを示す。図８において、ステップＳ１〜ステップＳ５までの処理は、図６に示した処理と同一である。この応用形態においては、ステップＳ５のスムージング処理が行われた後、スムージングデータの画像領域全体の画素を基に全体ヒストグラムを算出する（ステップＳ１１）。この全体ヒストグラムとは、スムージング処理後の画像（スムージングデータ）の全画素を母集団とするヒストグラムである。ここで、本実施形態で求める全体ヒストグラムの一例を図９（ａ）に示す。図９（ａ）において、横軸は画素値であり、縦軸は各値をとる画素数である。
【００３９】
次に、画像内の特定位置ごとにおける部分ヒストグラムを算出する（ステップＳ１２）。本実施形態では、画像全体をいくつかのブロックに分割し、その各ブロックを特定位置として、各ブロックに存在する画素を母集団とする部分ヒストグラムを算出する。例えば、図１０に示したように画像をＢ１〜Ｂ２０の２０個のブロックに分割した場合、２０個の異なる部分ヒストグラムが求められることになる。
【００４０】
続いて、全体ヒストグラム、部分ヒストグラムを用いて各画素の画素値を修正する（ステップＳ１３）。この画素値の修正は、画素Ｐの存在するブロックの部分ヒストグラムにおける相対的な位置を求め、全体ヒストグラムにおける同位置の画素値を画素Ｐの画素値とすることにより行う。両ヒストグラムを利用した画素値の修正は以下のような考え方に基づいて行われる。
【００４１】
ヒストグラムを求めることにより、その母集団における画素値の平均、および画素値の分散を算出することが可能となる。すなわち、全体ヒストグラムにおいては、全体平均Ａｇ、全体分散Ｄｇが求まり、部分ヒストグラムにおいては、ブロック平均Ａｐ、ブロック分散Ｄｐが求まる。このとき、画素Ｐの画素値の部分ヒストグラムにおける相対位置は、画素Ｐの画素値とブロック平均Ａｐの差Ｌｐと、ブロック分散Ｄｐの関係で表現することができる。ここで、あるブロックにおいて算出した部分ヒストグラムを図９（ｂ）に示す。図９（ｂ）からわかるように、この平均との差Ｌｐと、ブロック分散Ｄｐの関係を用いると、画素Ｐの画素値の部分ヒストグラム上における相対位置に対応する全体ヒストグラムにおける位置を求めることができる。全体ヒストグラムにおける位置と全体平均Ａｇとの差をＬｇとすると、両者の関係は、以下の〔数式３〕で表現される。
【００４２】
〔数式３〕
Ｌｇ／Ｌｐ＝Ｄｇ／Ｄｐ
【００４３】
画素Ｐの修正前の画素値をＶ、修正後の画素値をＶ´とすると、上記〔数式３〕は以下の〔数式４〕のように変形することができる。
【００４４】
〔数式４〕
（Ｖ−Ａｇ）／（Ｖ´−Ａｐ）＝ Ｄｇ／Ｄｐ
【００４５】
上記〔数式４〕はさらに変形することができ、結局、修正後の画素値Ｖ´は、以下の〔数式５〕により求められることになる。
【００４６】
〔数式５〕
Ｖ´＝ （Ｄｐ／Ｄｇ）×（Ｖ−Ａｇ）＋ Ａｐ
【００４７】
上記〔数式５〕を利用して画像上の全画素について、画素値の修正を行うことにより、光学系の収差による影響を除去したスムージングデータが得られることになる。このようにして得られたスムージングデータに対して上記ステップＳ６、ステップＳ７の処理を行うことにより、光学系の収差による影響を除去されると共に、入力装置が低解像度であっても鮮明な導管画像を得ることができる。
【００４８】
（システム構成）
次に、上記導管画像作成方法を実行するためのシステム構成について説明する。図１１は、本発明による導管画像作成システムの一実施形態を示す構成図である。図１１において、１１は画像入力手段、１２は導管部強調手段、１３は第１のスムージング処理手段、１４は第２のスムージング処理手段、１５は２値化手段、１６は画像出力手段、１７は全体ヒストグラム算出手段、１８は部分ヒストグラム算出手段、１９は画素修正手段である。画像入力手段１１は、異なる照明条件において撮影された２つの画像を入力する機能を有する。画像入力手段１１としては、撮影画像がデジタルデータとして既に用意されている場合は、そのデジタルデータを記録した電子記録媒体を読取り可能な読取り装置等で実現できる。また、デジタルカメラ自体を画像入力手段１１として利用し、撮影したデジタル画像データをそのまま入力するようにすることもできる。本発明においては、この画像入力手段１１が低解像度のものであっても鮮明な導管画像が得られることは上述の通りである。
【００４９】
導管部強調手段１２は、図６および図８のステップＳ２〜ステップＳ４の処理を実行する機能を有しており、入力された左画像、右画像それぞれに対して微分演算を行い、得られた両微分データの減算処理を行った後、得られた減算データの絶対値化を行うものである。第１のスムージング処理手段１３は、図６および図８のステップＳ５の処理を実行する機能を有しており、導管の流れと垂直な方向（上記の例では左右方向）にぼかし処理を行うものである。第２のスムージング処理手段１４は、図６および図８のステップＳ６の処理を実行する機能を有しており、導管の流れ方向（上記の例では上下方向）にぼかし処理を行うものである。２値化手段１５は、図６および図８のステップＳ７の処理を実行する機能を有しており、第２のスムージング処理手段１４によりぼかし処理された多値画像を所定の閾値で２値化するものである。画像出力手段１６は、作成された導管画像を出力する機能を有するものであり、ハードディスク等の外部記憶装置、携帯可能な電子記録媒体への記録装置、さらにネットワークを介して遠隔地へ送信するための通信装置等が適用可能である。全体ヒストグラム算出手段１７は、図８のステップＳ１１の処理を実行する機能を有しており、導管と直交する方向にスムージング処理された画像の全画素を母集団としてその画素値のヒストグラムを算出するものである。部分ヒストグラム算出手段１８は、図８のステップＳ１２の処理を実行する機能を有しており、スムージング処理された画像の一部の画素を母集団としてその画素値のヒストグラムを算出するものである。画素修正手段１９は、図８のステップＳ１３の処理を実行する機能を有しており、全体ヒストグラムおよび部分ヒストグラムを利用して、スムージング処理手段１３によりスムージング処理された画像の各画素の値を修正するものである。
【００５０】
図１１に示した導管部強調手段１２、第１のスムージング処理手段１３、第２のスムージング処理手段１４、２値化手段１５、全体ヒストグラム算出手段１７、部分ヒストグラム算出手段１８、画素修正手段１９は、コンピュータのＣＰＵ、メモリ等で実現され、各機能はコンピュータに搭載された専用のソフトウェアプログラムにより実現される。
【００５１】
以上、本発明の好適な実施形態について説明したが、本発明は上記実施形態に限定されず、種々の変形が可能である。例えば、上記実施形態では、ステップＳ１２において特定位置ごとにおけるヒストグラムとして、画像全体をいくつかのブロックに分割し、その各ブロックごとにヒストグラムを算出するようにしたが、各画素ごとに、当該画素の周辺の所定の範囲に属する全画素を母集団としたヒストグラムを算出するようにしても良い。
【００５２】
【発明の効果】
以上、説明したように本発明によれば、異なる照明角度で撮影した２つの画像それぞれに対しての微分演算、微分演算により得られる両微分データの減算、得られた減算データの絶対値化、絶対値化したデータに対して照明角度の変化方向についてのスムージング処理を順次行い、このスムージング処理により得られるスムージングデータに対して、前記スムージング処理の処理方向と直交する方向にスムージング処理を行うようにしたので、導管の細かさに対して入力装置が低解像度であるような場合であっても、ノイズの少ない導管の画像を得ることが可能となるという効果を奏する。
さらに、最初のスムージングにより得られるスムージングデータの全画素を分母とした全体ヒストグラムを算出すると共に、スムージングデータの特定の位置ごとに所定の画素数を分母とした部分ヒストグラムを算出し、全体ヒストグラムと部分ヒストグラムの関係に基づいて、スムージングデータ上の各画素の値を修正する処理を行うことにより、レンズ等の光学系の収差に起因する不均一性を除去し、撮影像の全体に渡って同程度で強調された導管の画像を得ることが可能となるという効果も奏する。
【図面の簡単な説明】
【図１】木材に表面の明るさが均一になるように照明を当てて撮影した画像の導管の部分と年輪の部分の輝度のプロファイルの例を示す図である。
【図２】左画像と右画像における輝度のプロファイルの相違を説明するための図である。
【図３】左画像における輝度のプロファイルと、それを微分演算した左微分データの例を示す図である。
【図４】右画像における輝度のプロファイルと、それを微分演算した右微分データの例を示す図である。
【図５】左微分データと右微分データを得た後の処理を説明するための図である。
【図６】本発明に係る導管画像作成方法の概要を示すフローチャートである。
【図７】ステップＳ６によるスムージング処理の効果を説明するための図である。
【図８】本発明に係る導管画像作成方法の応用例を示すフローチャートである。
【図９】全体ヒストグラムおよび部分ヒストグラムの一例を示す図である。
【図１０】画像を複数のブロックに分割した状態を示す図である。
【図１１】本発明に係る導管画像作成システムの構成図である。
【符号の説明】
１・・・導管
２・・・年輪
３・・・導管の像の輝度のプロファイル
４・・・年輪の像の輝度のプロファイル
１１・・・画像入力手段
１２・・・導管部強調手段
１３・・・第１のスムージング処理手段
１４・・・第２のスムージング処理手段
１５・・・２値化手段
１６・・・画像出力手段
１７・・・全体ヒストグラム算出手段
１８・・・部分ヒストグラム算出手段
１９・・・画素修正手段
Ａｇ・・・全体平均
Ａｐ・・・部分平均
Ｄｇ・・・全体分散
Ｄｐ・・・部分分散
【図１】
0004139601.tif 000002
【図２】
0004139601.tif 000003
【図３】
0004139601.tif 000004
【図４】
0004139601.tif 000005
【図５】
0004139601.tif 000006
【図６】
0004139601.tif 000007
【図７】
0004139601.tif 000008
【図８】
0004139601.tif 000009
【図９】
0004139601.tif 000010
【図１０】
0004139601.tif 000011
【図１１】
0004139601.tif 000012