JP 4229325 B2 2009.2.25(57)【 特 許 請 求 の 範 囲 】【請求項１】
任意の物理量の空間における分布に対応する濃淡画像を複数の小区画に分割した各区画が有する濃淡を表す濃淡データを各区画の位置を表す位置データと関連付けて含む濃淡画
像データ群の中の１つの区画の濃淡データ値を当該１つの区画の周囲を囲む全ての周囲区画の各濃淡データ値と個別に比較演算するステップと、前記１つの区画のデータ値が前記周囲区画の全てのデータ値より大きいか又は等しい場
合、又は前記１つの区画のデータ値が前記周囲区画の全てのデータ値より小さいか又は等しい場合、前記１つの区画のデータ値をピーク値と決定し、且つ前記１つの区画の位置をピーク位置と決定するピーク決定ステップと、10ピーク値として決定されたデータ値を有する複数の区画同士が互いに隣接又は隣り合う
場合任意の選択手法でもって前記複数の区画のうちの１つの区画を選択するステップとを備えるピーク検出画像処理方法。
【請求項２】
前記ピーク決定ステップが、前記１つの区画のデータ値が前記周囲区画の全てのデータ値より大きいか又は等しい場合に前記１つの区画のデータ値を極大のピーク値と決定し、且つ前記１つの区画のデータ値が前記周囲区画の全てのデータ値より小さいか又は等しい場合に前記１つの区画のデータ値を極小のピーク値と決定するステップを含む請求項１記載のピーク検出画像処理方法。
【請求項３】20

(2)JP 4229325 B2 2009.2.25

前記比較演算するステップ及び前記ピーク決定ステップを全ての区画に対して実行する
請求項１又は２記載のピーク検出画像処理方法。
【請求項４】
前記濃淡画像が２次元の矩形画像であり、前記各区画は、前記矩形画像の縦及び横方向それぞれを少なくとも３つに分割することにより形成されており、前記比較演算するステップは、前記１つの区画が前記矩形画像の外周に接していない区画である場合は前記１つの区画のデータ値を前記１つの区画を囲む８個の区画のデータ値と比較演算し、前記１つの区画が前記矩形画像の縦及び横方向のうちの一方の方向のみの外周に接している区画である場合は前記１つの区画のデータ値を前記１つの区画を囲む５10個の区画のデータ値と比較演算し、前記１つの区画が前記矩形画像の角にある区画である
場合前記１つの区画のデータ値を前記１つの区画を囲む３個の区画のデータ値と比較演算する
請求項１から３のいずれか一項に記載のピーク検出画像処理方法。
【請求項５】
前記ピーク決定ステップにより決定されたピーク値の数を計数するステップを更に備える請求項１から４のいずれか一項に記載のピーク検出画像処理方法。
【請求項６】
前記濃淡画像が、撮像管又はＣＣＤ等の光画像検出器により光学的に取得され、前記濃淡画像データ群の濃淡データ及び位置データが、前記の取得された濃淡画像から20ディジタル値として生成される請求項１から５のいずれか一項に記載のピーク検出画像処理方法。
【請求項７】
前記の決定されたピーク値及びピーク位置を記録するステップを更に備える請求項１から６のいずれか一項に記載のピーク検出画像処理方法。
【請求項８】
コンピュータに、請求項１から７のいずれか一項に記載のピーク検出画像処理方法を実行させるピーク検出画像処理プログラム。
【請求項９】
任意の物理量の空間における分布に対応する濃淡画像を複数の小区画に分割した各区画30が有する濃淡を表す濃淡データを各区画の位置を表す位置データと関連付けて含む濃淡画
像データ群の中の１つの区画の濃淡データ値を当該１つの区画の周囲を囲む全ての周囲区画の各濃淡データ値と個別に比較演算する手段と、
前記１つの区画のデータ値が前記周囲区画の全てのデータ値より大きいか又は等しい場
合、又は前記１つの区画のデータ値が前記周囲区画の全てのデータ値より小さいか又は等しい場合、前記１つの区画のデータ値をピーク値と決定し、且つ前記１つの区画の位置をピーク位置と決定するピーク決定手段と、ピーク値として決定されたデータ値を有する複数の区画同士が互いに隣接又は隣り合う
場合任意の選択手法でもって前記複数の区画のうちの１つの区画を選択する手段とを備えるピーク検出画像処理装置。

40【請求項１０】前記ピーク決定手段が、前記１つの区画のデータ値が前記周囲区画の全てのデータ値より大きいか又は等しい場合に前記１つの区画のデータ値を極大のピーク値と決定し、且つ前記１つの区画のデータ値が前記周囲区画の全てのデータ値より小さいか又は等しい場合に前記１つの区画のデータ値を極小のピーク値と決定する請求項９記載のピーク検出画像処理装置。
【発明の詳細な説明】【技術分野】【０００１】
本発明は、任意の物理量の空間における分布に対応する濃淡画像を複数の小区画に分割

50

(3)

JP 4229325 B2 2009.2.25

した各区画が有する濃淡を表す濃淡データを各区画の位置を表す位置データと関連付けて備える濃淡画像データ群の中で濃淡データのピーク値、その位置及びピークの数を検出するピーク検出画像処理方法に関する。
【０００２】
本発明は、工業品として生産される粉または細かな粒子が多数集まった粉体、粒子状物質等の画像計測で得られる画像から、画像処理・解析により（１）粉、粒子の数量、または画像中の粒子の輝度から（２）粉、粒子の大きさ、粒子径、及び（３）粒子の位置を検出する際に利用できる。または、空気清浄産業等の環境科学研究分野で、大気中または液中に浮遊する微粒子等を観測した斑点濃淡画像から微粒子の数、粒径及び画像中での位置を求める際に本発明を利用することができる。 さらに半導体シリコンウェハーのように10平板上の凹凸検査やそれらに付着する微小な粉塵等のゴミの検出にも使用できる。
【背景技術】【０００３】
従来の手法では、本発明が対象とするような斑点が多数集まった濃淡画像中に多数ある斑点のピーク値と数を正確に検出あるいは取得するために、一般的には数々の特殊かつ煩雑なフィルタ処理を多用している。その一例を挙げると、まず画像を移動平均法やメディアンフィルタなどのローパス特性のフィルタで処理することによって画像を平滑化して、エンベロープ画像を作成する。次に、元の画像との差分をとって、スポットまたは点状の部分だけを抜き出す処理を経る。そして、この差分画像に対してしきい値処理で画像を２値化処理した後、細線化処理を一つの画素まで削り込むために複数回繰り返し施す。最後20に、その一つの画素となった画素を数えることでピークを有する画素の数とピーク値を得る手法が必要となる。このとき、斑点ごとにピーク値を取得するとき、差分画像の斑点の領域の最高値（ピーク値）を残すような処理が必要となる。
【０００４】
本発明に関連する従来技術を記載したものとしては下記のものが挙げられる。
【特許文献１】特開２００２−２６２０９３号公報【特許文献２】特開２００３−２５９１２９号公報【非特許文献１】Ｋ．Ｒ．キャストルマ著、斎藤恒雄訳「ディジタル画像処理」科学技術
出 版 、 2002年 3月 刊 行【非特許文献２】Ａ. Ｒｏｓｅｎｆｅｌｄ著、長尾真監訳「ディジタル画像処理」近代科30

学 社 、 1978年 12月 刊 行【 非 特 許 文 献 ３ 】 河 田 聡 ／ 南 茂 夫 著 「 科 学 計 測 の た め の 画 像 デ ー タ 処 理 」 CQ出 版 社 、 1994
年 ４ 月 30日 刊 行【発明の開示】【発明が解決しようとする課題】【０００５】
従来の画像処理法で画像中に多数ある斑点のピーク値と数を求める場合、上述のようにフィルタ、２値化、細線化等の複数の画像処理法を繰り返し多用し、適用しなければならない。このため画像処理の手順、プロセスが煩雑になり、かつ膨大な画像処理時間が必要になるという問題が発生する。さらに、数々の画像処理により生のピーク値が保存・維持40されず、変化する可能性が発生する。
【０００６】
本発明は、任意の物理量の空間における分布に対応する濃淡画像中に多数ある斑点のピーク値及びその位置、並びに所望の場合ピーク値の数を求める場合の画像処理において、複数の画像処理法の適用手順、処理プロセスの煩雑さとそれに必要とされる膨大な画像処理時間の二つの大きな問題を解決して、簡単で且つ短時間で処理できるようにすることを課題とする。
【課題を解決するための手段】【０００７】
上記課題を解決するため、本発明のピーク検出画像処理方法は、任意の物理量の空間に50

(4)JP 4229325 B2 2009.2.25

おける分布に対応する濃淡画像を複数の小区画に分割した各区画が有する濃淡を表す濃淡データを各区画の位置を表す位置データと関連付けて含む濃淡画像データ群の中の１つの区画の濃淡データ値を当該１つの区画の周囲を囲む全ての周囲区画の各濃淡データ値と個別に比較演算するステップと、前記１つの区画のデータ値が前記周囲区画の全てのデータ値より大きいか又は等しい場合、又は前記１つの区画のデータ値が前記周囲区画の全てのデータ値より小さいか又は等しい場合、あるいはこれら両方の場合、前記１つの区画のデータ値をピーク値と決定し、且つ前記１つの区画の位置をピーク位置と決定するピーク決定ステップとを備えることを特徴とする。
【０００８】
前記比較演算するステップ及び前記ピーク決定ステップが全ての区画に対して実行され10ることが好ましい。
本発明の一局面においては、前記濃淡画像が２次元の矩形画像であり、前記区画は、前記矩形画像の縦及び横方向それぞれを少なくとも３つに分割することにより形成されており、前記比較演算するステップは、前記１つの区画が前記矩形画像の外周に接していない区画である場合は前記１つの区画のデータ値を前記１つの区画を囲む８個の区画のデータ値と比較演算し、前記１つの区画が前記矩形画像の縦及び横方向のうちの一方の方向のみの外周に接している区画である場合は前記１つの区画のデータ値を前記１つの区画を囲む５個の区画のデータ値と比較演算し、前記１つの区画が前記矩形画像の角にある区画である場合前記１つの区画のデータ値を前記１つの区画を囲む３個の区画のデータ値と比較演算することが好ましい。

20【０００９】
本発明のピーク検出画像処理方法は更に、ピーク値として決定されたデータ値を有する
複数の区画同士が互いに隣接又は隣り合う場合任意の選択手法でもって前記複数の区画のうちの１つの区画を選択するステップを備えることが好ましい。
【００１０】
本発明のピーク検出画像処理方法は更に、前記ピーク決定ステップにより決定されたピーク値の数を計数するステップを備えることが好ましい。
本発明の別の局面においては、前記濃淡画像が、撮像管又はＣＣＤ等の光画像検出器により光学的に取得され、前記濃淡画像データ群の濃淡データ及び位置データが、前記の取得された濃淡画像からディジタル値として生成されることが好ましい。

30【００１１】
本発明のピーク検出画像処理方法は更に、前記の決定されたピーク値及びピーク位置を記録するステップを備えることが好ましい。
上記課題は、コンピュータに、上記の種々の態様のピーク検出画像処理方法を実行させるピーク検出画像処理プログラムにより解決される。
【００１２】
上記課題を解決するため、本発明のピーク検出画像処理装置は、任意の物理量の空間における分布に対応する濃淡画像を複数の小区画に分割した各区画が有する濃淡を表す濃淡データを各区画の位置を表す位置データと関連付けて含む濃淡画像データ群の中の１つの区画の濃淡データ値を当該１つの区画の周囲を囲む全ての周囲区画の各濃淡データ値と個40別に比較演算する手段と、前記１つの区画のデータ値が前記周囲区画の全てのデータ値より大きいか又は等しい場合、又は前記１つの区画のデータ値が前記周囲区画の全てのデータ値より小さいか又は等しい場合、あるいはこれら両方の場合、前記１つの区画のデータ値をピーク値と決定し、且つ前記１つの区画の位置をピーク位置と決定するピーク決定手段とを備えることを特徴とする。
【発明の効果】【００１３】
本発明は、上記の構成を備えることにより、任意の物理量の空間における分布に対応する濃淡画像における濃淡データのピーク値及びその位置を簡単に且つ短時間で検出することができる。

50

(5)

JP 4229325 B2 2009.2.25【００１４】
また、本発明は、上記の構成を備えることにより、従来画像中に多数あるピーク値の検出や解析に不可欠となっていた多種多様な画像処理、煩雑なプロセスが不必要となり、且つ多大な処理計算時間が大幅に短縮できる。さらに、従来では不可能であった画像中のピーク値のリアルタイム検出、解析が可能となる。
【００１５】
例えば、濃淡画像に斑点が多数ある画像である場合、その画像中の各斑点のピーク値とピークを有する画素の数、及び画面上のピークの位置を同時に求めることができる。
また、本発明の一局面では、同じピーク値を有する画素が隣接し、画面上連続的に連なりあう領域を形成する場合、その中の一つの画素を選択する任意の方法または規準を付加10することによってピーク値を一つに絞り込むことができ、ピーク値を有する斑点の数を正確に取得できる。
【発明を実施するための最良の形態】【００１６】
本発明の好適な実施形態について図面を参照して以下に説明する。
図１は、本発明のピーク検出画像処理装置の好適な実施形態の構成を概略示したブロック図である。ピーク検出画像処理装置は、例えば、撮像管又はＣＣＤ（電荷結合素子）カメラ等の光画像検出器で光学的に取得された画像をディジタル形式で受け取る画像入力部
１０、画像入力部１０で取得されたディジタル画像、及び当該ディジタル画像の中の濃淡のピークを検出するピーク検出画像処理プログラムを格納し、また処理中のワークエリア20を含む外部記憶装置１２、画像処理のため外部記憶装置１２から呼び出されたピーク検出
画像処理プログラム及び処理すべき画像のデータ、あるいは処理中に発生するデータを記憶するメイン・メモリ１４、メイン・メモリ１４のピーク検出画像処理プログラムに従ってピーク検出処理を実行する中央処理装置（ＣＰＵ）１６、ピーク検出処理を実行させるために必要な制御情報等を入力するためのキーボード等の入力装置１８、及びピーク検出処理の中間結果、最終結果、入力された制御情報等を表示するためのディスプレイ２０を備える。これら各構成要素は、通常のバスを介して相互接続されている。
【００１７】
次に、本発明のピーク検出画像処理方法及びプログラムの好適な実施形態を図面を参照して以下に説明する。

30

本発明の処理対象である濃淡画像は、例えば、図４を参照して後述するように室内の任意の空間に浮遊する微粒子にレーザー光を照射し、その計測される微粒子からの散乱光を
画 素 数 324,064個 の Ｃ Ｃ Ｄ カ メ ラ で 撮 影 す る こ と に よ り 得 ら れ た 濃 淡 画 像 で あ る 。 図 ４ の濃淡画像中の各画素の濃淡は、散乱光の強度の強弱に対応し、各画素の濃淡を表すディジタル値は散乱光の強弱を表している。なお、本発明の処理対象である濃淡画像は、このようなＣＣＤカメラ等で撮影して得られた画像に限られず、任意の手段を用いて、２次元あるいは３次元に分布している任意の対象を撮像することにより得られ、その結果としての画像が濃淡で表現されるものであれば、いずれものでよい。また、濃淡画像は、ＣＣＤの場合のように直接ディジタル形式であってもよく、また通常のカメラの場合のようにアナログ形式であってもよい。後者の場合は、濃淡画像のアナログ・データをディジタル・デ40ータに変換する任意の手段を外部又は本発明のピーク検出画像処理装置に設け、その手段によりアナログ値をディジタル値に変換すればよい。また、画像入力部１０は、他の任意のシステムから直接又は通信回線等の通信媒体を介して送られてくる画像データを受け取るようにしてもよい。
【００１８】
本発明によるピーク検出画像処理方法を画像中にスポットまたは点状の画像が多数集まった斑点状の画像データとして取得された図２に示されるようなｎ行ｍ列のサンプル行列データに適用する。なお、ｎ及び／又はｍは２以上の整数であり、３以上であることが好ましい。図３にピーク検出画像処理フローを示す。ピーク検出画像処理の手順は、図３のピーク検出画像処理フローに示すように処理がスタートすると同時に、ステップＳ１にお50

(6)JP 4229325 B2 2009.2.25

いて、処理の対象となる画像行列データを選択する。次に、ステップＳ２において、まだピ ー ク 検 出 画 像 処 理 に 選 択 さ れ て い な い 、 即 ち ま だ 比 較 さ れ て い な い 任 意 の 一 つ の i行 j列の 画 素 デ ー タ 値 ｘ (i,j)を 選 択 す る 。 次 い で 、 ス テ ッ プ Ｓ ３ に お い て 、 選 択 さ れ た 画 素 デー タ 値 ｘ (i,j)を そ の 周 囲 に あ る ８ つ の 画 素 デ ー タ 値 ｘ (i‑1,j‑1)、 ｘ (i‑1,j)、 ｘ (i‑1,j+1)、 ｘ (i,j‑1)、 ｘ (i,j+1)、 ｘ (i+1,j‑1)、 ｘ (i+1,j)、 ｘ (i+1,j+1)と 各 々 比 較 し 、 画 素デ ー タ 値 ｘ (i,j)が 周 囲 の 画 素 デ ー タ 値 と 以 下 の 関 係 式 （ １ ） を 満 た す と き 、 す な わ ち 周
囲 の 画 素 デ ー タ 値 各 々 と 比 較 し て こ の 画 素 デ ー タ 値 ｘ (i,j)が 全 て の 周 囲 の 画 素 デ ー タ よ
り大きいかまたは周囲の画素データの中で最も大きいものと等しいとき、その画素データ値 ｘ (i,j)を ピ ー ク 値 と し て 検 出 し て 、 ス テ ッ プ Ｓ ４ に 進 み 、 こ れ を ピ ー ク 値 を 有 す る 画
素 デ ー タ 値 ｘ (i,j)と し て 外 部 記 憶 装 置 １ ２ に 記 録 、 保 存 す る 。

10【００１９】【数１】

20【００２０】
一 方 、 ス テ ッ プ Ｓ ３ に お い て 、 画 素 デ ー タ 値 ｘ (i,j)が 、 周 り の 画 素 デ ー タ 値 の ど れ かよりも小さい場合にはこの画素データ値をピーク値を有する画素データとしては記録、保存しない。ピーク検出画像処理は、次いで、ステップＳ５に進み、ｎ行ｍ列の画像行列デ30ータの全画像データ（ｎ×ｍ個）を全部処理したか判定する。未処理の画像データがある場合には、ステップＳ２に戻り、前述の処理を繰り返す。ステップＳ５において、全画像データが処理されたと判定された場合、ピーク検出画像処理は終了する。
【００２１】
上記実施形態のピーク検出画像処理方法は、簡単に述べると、画像行列データの全画素データについて任意の一つの画素のデータ値をその画素の周囲にあるか、または隣接する画素データ値一つ一つと比較演算することを全画素のデータに施すことにより、周囲にある 画 素 の デ ー タ 値 よ り も 高 い 値 の デ ー タ 、 ピ ー ク 値 (極 大 値 )を 有 す る 画 素 を 見 つ け 出 す 処理法である。
【００２２】

40

前述のような処理を行列データの全画素データに施すことにより、画像中にあるスポットまたは点の画素のピーク値を求めることが可能となり、同時にその数及び画素の位置も求められる。
【００２３】
図 ２ に 示 さ れ る ｎ 行 ｍ 列 の 画 像 行 列 デ ー タ の 角 に 位 置 す る 画 素 デ ー タ で あ る ｘ (1,1)、ｘ (n,1)、 ｘ (1,m)、 ｘ (n,m)の 場 合 に は 、 各 画 像 デ ー タ は そ れ ら の 画 素 を 囲 む 三 つ の 画 素デ ー タ と 比 較 さ れ る 。 例 え ば 、 ｘ (1,1)の 場 合 、 ｘ (2,1)、 ｘ (1,2) ｘ (2,2)と そ れ ぞ れ 比較 処 理 す る 。 さ ら に 、 画 像 行 列 デ ー タ の 端 に 位 置 す る １ 列 の 2行 １ 列 x(2,1)か ら n‑1行 １ 列x(n‑1,1)、 ｍ 列 の 2行 ｍ 列 x(2,m)か ら n‑1行 ｍ 列 x(n‑1,m)、 １ 行 の 1行 2列 x(1,2)か ら 1行 m‑1列 x(1,m‑1)、 及 び n行 の n行 ２ 列 x(n,2)か ら n行 m‑1列 ｘ (n,m‑1)の 画 素 デ ー タ の 場 合 に は 、50(7)JP 4229325 B2 2009.2.25

それらの各画像データをそれらを囲む五つの画素データと比較処理する。
【００２４】
仮 に 、 画 像 処 理 の 結 果 、 図 ２ で ｘ (i,j)、 ｘ (i,j+1)、 ｘ (i+1,j)、 ｘ (i+1,j+1)が 同 じ 値を有し、かつこれら隣接する四つの画素全てが同一のスポット又は斑点のピーク値を持つものであると記録、保存された場合には、任意の手法によりこの中の一つをピーク値として記録、保存することによりピークの数を正確に求めることができる。例えば、このような 場 合 に 行 と 列 の 数 値 が 大 き い 画 素 を 選 択 す る よ う な 手 法 を 付 加 す れ ば 、 ｘ (i,j)、 ｘ (i,j+1)、 ｘ (i+1,j)、 ｘ (i+1,j+1)の う ち ｘ (i+1,j+1)が ピ ー ク 値 を 有 す る 画 素 と し て 外 部 記
憶装置１２に記録、保存され、ピークを有する画素の数を正確に取得することが可能となる。

10【００２５】
次に、図３に示すピーク検出画像処理フローを、室内の任意の空間に浮遊する微粒子にレーザー光を照射し、その散乱光を撮影した画像に適用した場合の測定例又は実験データについて説明する。図４は、室内の任意の空間に浮遊する微粒子にレーザー光を照射し、そ の 計 測 さ れ る 微 粒 子 か ら の 散 乱 光 を 画 素 数 324,064個 の Ｃ Ｃ Ｄ カ メ ラ で 撮 影 し た 画 像 を示す。画像中に多数あるスポットまたは点状の斑点は、気中の微粒子一つ一つからの散乱光 で あ る 。 こ の 画 像 デ ー タ は 、 494行 656列 の 行 列 デ ー タ で あ る 。 そ の 一 部 分 を 図 ５ に 示 す。散乱光の強度、画面では白黒の濃淡で表され、最も黒く見える所は散乱光の強度が最も弱く、最も白く見える所は散乱光の強度が最も強い。図５は、最小強度であるＣＣＤの暗電流のみの状態を０とし、最大強度を２1 2 ＝４０９６として、図４に示す画像のそれぞれ20の画素の濃淡を相対値で表したものである。
【００２６】
その図５の行列データを三次元グラフに示すと図６のようになる。図６に多数見える山形状または尖塔状のものの一つは、微粒子一つからの散乱光を示し、これが図４の斑点一つを表す。この山の数は、微粒子の数を表し、その高さ、ピーク値は微粒子の大きさを示す。
【００２７】
図４の画像データである図５の行列データに本発明のピーク検出画像処理法を適用したものが図７の行列データである。図７の行列データ中、便宜のためピーク値以外の画素の行列データは０として表示し、ピーク値を有する画素にあたる行列データにのみピーク値30を記した。これにより、ピーク値を有する画素の数とそのピーク値、及びその画素の画面上での位置を求めることができる。
【００２８】
また、希望する場合には、この画像処理後の行列データをまとめてヒストグラム（度数分布グラフ）にした図８やピーク値の範囲ごとにその範囲に含まれるピークの数を表した図９に示す表を容易に得ることが可能となる。
【００２９】
本発明のピーク検出画像処理装置は、図１に示される構成に限られず、通常のパーソナル・コンピュータ、マイクロプロセッサ、ワークステーション、汎用コンピュータ等の任意の形式のコンピュータを用いて、前述のピーク検出画像処理プログラムを実行させ得る

40。
【００３０】
本発明は、前述の実施形態のように、ピーク値が極大値である場合に限られず、比較演算において、周囲の画素のデータ値の最も小さいデータ値と等しいかまたは全てのデータ値 よ り も 小 さ い 場 合 は そ の 画 素 を 極 小 値 を 有 す る デ ー タ 値 ｘ (i,j)と し て 記 録 ・ 保 存 す る
操作を繰り返し行うことによって、ピーク値を有する画素の検出と同様に画像中のマイナス側に窪んだ輝度分布を有する斑点の極小値（マイナスのピーク値）の数、ピーク値及び画面上での位置を取得することも可能である。
【００３１】
上記の実施形態においては、濃淡画像として図４に示されるように、その外形が矩形で50

(8)JP 4229325 B2 2009.2.25

ある場合について説明したが、本発明は、濃淡画像の形状は、矩形に限られず、円形等いずれの２次元の形状であってよく、その濃淡画像を小区画に分割した区画（上記の実施形態における画素に対応）の位置が特定されていればよい。また、区画は、上記の実施形態の矩形に限定されず、区画の大きさが一様であれば任意の形状でよい。
【００３２】
更に、上記の実施形態においては、濃淡画像は、図４に示されるように２次元の画像であるが、本発明は、濃淡画像が３次元の構成の場合にも適用し得る。その場合は、区画としては、例えば、小さい立体構造の濃淡画像を小さく分割した立方体等にし、その立方体の濃淡に対応するデータ値のピーク値を上記実施形態と同様にして当該立方体を立体的に囲む隣りの立方体のデータ値と比較することにより検出し、その立方体の位置を特定すれ10ばよい。
【図面の簡単な説明】【００３３】【図１】図１は、本発明のピーク検出画像処理装置の好適な実施形態の構成を概略示したブロック図である。
【図２】図２は、画像中にスポットまたは点状の画像が多数集まった斑点状の画像データとして取得されたｎ行ｍ列のサンプル行列データを示す。
【図３】図３は、本発明の好適な実施形態によるピーク検出画像処理フローを示す。
【図４】図４は、室内の任意の空間に浮遊する微粒子にレーザー光を照射し、その計測され る 微 粒 子 か ら の 散 乱 光 を 画 素 数 324,064個 の Ｃ Ｃ Ｄ カ メ ラ で 撮 影 し た 画 像 を 示 す 。

20【図５】図５は、図２に示す画像の任意の部分を表したサンプル行列データを示す。
【図６】図６は、図４あるいは図５の散乱光画像の任意の部分の３次元分布図を示す。
【図７】図７は、図４の画像データである図５の行列データに本発明のピーク検出画像処理法を適用して得られたた行列データを示す。
【図８】図８は、図７に示す画像処理後の結果をピーク値ごとにその数、度数を表したヒストグラムを示す。
【図９】図９は、図７に示すデータをピーク値の範囲に対して、その範囲に含まれるピークの数を表した図表を示す。
【符号の説明】【００３４】
１０

画像入力部

１２

外部記憶装置１４

メイン・メモリ

１６

中央処理装置（ＣＰＵ）

１８

入力装置

２０

ディスプレイ

30

(9)
【図１】

【図３】

【図２】

【図５】

【図７】

【図８】

JP 4229325 B2 2009.2.25

(10)
【図９】

【図４】

JP 4229325 B2 2009.2.25

(11)
【図６】

JP 4229325 B2 2009.2.25

(12)JP 4229325 B2 2009.2.25

フロントページの続き(74)代理人 100093713
弁理士 神田 藤博(74)代理人 100091063
弁理士 田中 英夫(74)代理人 100102727
弁理士 細川 伸哉(74)代理人 100117813
弁理士 深澤 憲広(74)代理人 100123548
弁理士 平山 晃二(72)発明者 大図 章
茨城県那珂郡東海村白方字白根２番地の４ 日本原子力研究所 東海研究所内審査官 八島 剛(56)参考文献 特開平０５−２７３１１０（ＪＰ，Ａ）
特開２０００−１５５０８７（ＪＰ，Ａ）
特開平０８−０５４３３３（ＪＰ，Ａ）
特開平０６−１６０５０７（ＪＰ，Ａ）(58)調査した分野(Int.Cl.，ＤＢ名)
Ｇ０１Ｂ１１／００−１１／３０
Ｇ０１Ｎ１５／００−１５／１４
Ｇ０１Ｎ２１／００−２１／９５８
Ｇ０６Ｔ１／００−９／４０