JP3728184B2 - 画像処理システム - Google Patents
画像処理システム Download PDFInfo
- Publication number
- JP3728184B2 JP3728184B2 JP2000190747A JP2000190747A JP3728184B2 JP 3728184 B2 JP3728184 B2 JP 3728184B2 JP 2000190747 A JP2000190747 A JP 2000190747A JP 2000190747 A JP2000190747 A JP 2000190747A JP 3728184 B2 JP3728184 B2 JP 3728184B2
- Authority
- JP
- Japan
- Prior art keywords
- differential
- image
- grayscale image
- value
- processing system
- 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
- Image Processing (AREA)
- Image Analysis (AREA)
Description
【発明の属する技術分野】
本発明は、画像処理システムに関し、特に、画質に左右されることなく、被測定物のエッジ位置を正確かつ安定して検出することを可能にした画像処理システムに関する。
【0002】
【従来の技術】
CCDカメラなどで撮像された被測定物の濃淡画像に含まれるエッジ部の位置検出を行う場合、まず、濃淡画像を微分することにより微分画像を作成し、この微分画像のピーク位置(微分ピーク位置)を求め、この微分ピーク位置近傍の微分曲線に対して近似曲線を当てはめ、この近似曲線のピーク位置をエッジ位置とする方法が広く知られている。この方法は、微分ピーク位置近傍の微分値分布が急峻である場合には有効であるが、微分ピーク位置近傍の微分値分布がなだらかな場合や撮像された濃淡画像のピントが外れていた場合などでは、微分ピークが相対的に低くなり、なだらかな丘のような形状の微分曲線となる。この結果、微分ピーク位置近傍の微分値がほぼ等しくなり、撮像条件によっては微分ピーク位置が一定しなくなるため、安定してエッジ位置を検出することができないという問題があった。特に、コントラストが低い場合は、ピーク位置の濃度値と左右のすそ野領域の濃度値の差が小さく、近似曲線に近似する場合に、すそ野領域の濃度ばらつき(雑音)の影響が大きくなり、ピーク位置検出の安定性が低下する問題があった。
【0003】
そこで、特開平9−257422号公報に開示する方法では、まず、被測定対象の濃淡画像に含まれるエッジのおおよその位置(エッジ概略位置)を検出し、次に、濃淡画像を微分した微分画像においてエッジ概略位置の近傍から微分ピーク位置を探索し、続いて、微分画像のピーク値のレベルに対して所定の割合のレベルを閾値として設定し、このピーク値の周辺の微分値分布を構成する画素のうち前記閾値を超えるレベルの画素の数が所定値を下回ったときに当該微分値分布が急峻であると判別し、一方、前記閾値を超える画素の数が所定値以上であるとき、あるいは前記微分ピーク値が閾値を下回るときに当該微分値分布がなだらかであると判別する。そして、微分分布が急峻であると判別したときは、前記微分ピーク位置近傍に近似曲線を当てはめそのピーク位置をエッジ位置とし、また、なだらかであると判別したときは前記微分ピーク位置近傍の重心計算を行い、その重心位置をエッジ位置とし、エッジ検出を行うようにしている。これにより、微分ピーク位置近傍の微分値分布がなだらかな場合でも、エッジ位置を検出可能にするとされている。
【0004】
【発明が解決しようとする課題】
しかし、特開平9−257422号公報で開示されている方式では、微分ピーク値により2つの方法のエッジ位置検出方法を切り替えているため、微分ピーク値や微分ピーク位置近傍の微分値分布形状によっては、この2つの方法のどちらが選択されるかわからず、したがって、この2つの方法の遷移部分でのエッジ位置を安定して求めることができないという問題がある。
【0005】
本発明は、この従来技術の問題点を解決するためになされたものであり、エッジ部近傍の微分値分布が急峻である場合だけでなく、ピントが外れたりエッジ部近傍の微分値分布がなだらかな場合でも、撮影された濃淡画像に含まれるエッジ位置を、切り替え処理のない1つの画像処理シーケンスにより、正確かつ安定して検出することが可能な画像処理システムを提供することを目的とする。
【0006】
【課題を解決するための手段】
上記目的を達成するために、請求項1にかかる発明では、被測定物を撮影することにより被測定物の濃淡画像データを得るCCDカメラと、このCCDカメラにて得られた濃淡画像データを格納する濃淡画像メモリと、この濃淡画像メモリに格納された濃淡画像データに対して演算処理を行うCPUと、を有する画像処理システムにおいて、前記CPUでは以下の一連の演算処理、すなわち、前記濃淡画像メモリに格納された濃淡画像データに対して微分処理を施すことにより微分画像を作成し、この微分画像上を走査し、走査された画素の微分値の中から最も値の大きいものを抽出し、この最大の微分値を有する画素の位置をピーク位置とし、このピーク位置の近傍の微分画像に対して非線形演算を行い、この非線形演算により得られた分布曲線に対して近似曲線とのマッチング処理を施すことにより被測定物のエッジ位置を算出する、一連の演算処理が行われることを特徴とする画像処理システムを提供した。
【0007】
かかる構成によれば、非線形演算により得られた分布曲線はその元になった微分画像の分布曲線に比して、分布曲線のピーク領域がより大きく、一方、分布曲線のすそ野領域はより小さくなるので、すそ野領域の変動の影響が相対的に小さくなり、その結果、安定してエッジ位置を検出できるようになる。
【0008】
なお、前記非線形演算は、具体的には、微分画像の微分値の二乗を演算してもよいし(請求項2)、微分画像の微分値が大きいほど大きな定数を乗算するようにしてもよい(請求項3)。
【0009】
【発明の実施の形態】
以下、図面を参照して本発明の一実施形態について説明する。図2は、本発明の一実施形態に係る画像処理システムを示すブロック図である。この画像処理システムは、CCDカメラ1、コンピュータ本体2、マウス3、キーボード4、及びCRT画面8により構成されている。CCDカメラ1により撮像された濃淡画像データは、インターフェース(以下「I/F」と記す)5を介して濃淡画像メモリ6に格納される。濃淡画像メモリ6に格納された濃淡画像データは、表示制御部7を介してCRT画面8に表示される。一方、データ入力機器としてのマウス3及びキーボード4から入力される位置情報は、それぞれI/F9、I/F10を介してCPU11に入力される。CPU11では後述するエッジ検出方法にかかる一連の演算処理が行われる。なお、プログラムメモリ12にはCPU11にて演算される処理を規定するプログラムが格納されており、また、ワークメモリ12にはCPU11における演算過程で使用されるデータが一時的に格納される。
【0010】
図1は、前述の図2で示した画像処理システムにおいて行われるエッジ検出の演算フローである。また、図3〜6は、演算フローの各段階における具体的な画像処理データを示したものである。ここで、○、●は画素毎にサンプリングした画素値を表す。
【0011】
まず、図2のCCDカメラ1からI/F5を経由して濃淡画像メモリ6に濃淡画像が記録される(図1のステップ▲1▼)。このときの濃淡画像及び濃度値分布を図3に示す。なお、この濃度値分布において、横軸は濃淡画像の横方向のX座標、縦軸は濃度値をそれぞれ示す。
【0012】
次に、CPU11において、濃淡画像メモリ6に記録されている濃淡画像データに対して微分フィルタ処理を施すことにより微分画像を作成し(図1のステップ▲2▼)、これをワークメモリ12に一時的に記憶する。このときの微分画像及び微分値分布を図4に示す。なお、この微分値分布において、横軸は濃淡画像の横方向のX座標、縦軸は微分値をそれぞれ示す。
【0013】
さらに、この微分画像上をX軸方向に走査し、微分値のピーク位置X0 を探索する(図1のステップ▲3▼)。ここでの探索は、走査された画素の微分値の中から最も値の大きいものを抽出し、この最大の微分値を有する画素の位置をピーク位置X0 とする。
【0014】
次に、ステップ▲2▼で求められた図4の微分値分布で示される波形において、ステップ▲3▼で求められたピーク位置X0 の近傍のn個の画素に対して、非線形演算を行う(図1のステップ▲4▼)。ただし、このとき微分値の値が大きいほど、より演算結果が大きくなるような非線形演算処理を行う。このときの演算結果の分布を図5に示す。この図5の例では、微分後の各画素の値を二乗することでその効果を出すようにしている。図5において、横軸は濃淡画像の横方向のX座標、縦軸は微分値の二乗をそれぞれ示す。この結果、ピークを含むピーク近傍領域はその領域以外に比して波形の急峻度が増し、ピーク濃度値とすそ野領域のダイナミックレンジが拡大することになる。
【0015】
最後に、この結果すなわち図5に示した分布曲線に対して図6に示すような近似曲線とのマッチング処理を行うことによりピーク位置を計算すると、エッジの位置を正確に算出することができる(図1のステップ▲5▼)。例えば、近似曲線をガウス分布曲線とすると、算出すべきエッジ位置XMAX は式(1)から求められる。
【0016】
【数1】
【0017】
式(1)において、X0 はステップ▲3▼で求められたピーク位置、f0 はステップ▲4▼で求められたピーク位置X0 における微分値の二乗(以下「画素値」と記す)をそれぞれ示す。また、f1 及びf2 はそれぞれピーク位置X0 からプラス方向に1画素分及び2画素分離れた画素における画素値、f-1及びf-2はそれぞれピーク位置X0 からマイナス方向に1画素分及び2画素分離れた画素における画素値である。
【0018】
以下、上述のエッジ検出アルゴリズムにガウス近似を採用する場合を例にして、本実施形態の効果を述べる。ここで、ランダムな雑音αが各画素値(fi 、i=0、±1、±2)に重畳していると仮定すると、雑音による誤差の最大値は式(2)で表される。
【0019】
【数2】
【0020】
一方、本実施形態におけるステップ▲4▼での非線形演算終了後の雑音による誤差の最大値は式(3)で表される。
【0021】
【数3】
【0022】
次いで、式(2)から求められるNaと式(3)から求められるNbとを比較するわけであるが、「f0 >>α」(f0 がαに比して非常に大きい)という条件を考慮すると、式(2)及び式(3)はそれぞれ式(2)′及び式(3)′のように書き直すことができる。
【0023】
【数4】
【0024】
さらに、画像の性質上、「f0 >(f-1、f1 )>(f-2、f2 )」であるから、式(4)が成り立つ。
【0025】
【数5】
【0026】
この式(4)の関係を前記式(3)′に適用し、式(2)′から求められるNa′と式(3)′から求められるNb′とを比較すると、「Na′>Nb′」となり、結局「Na>Nb」となる。ゆえに、本実施形態における非線形演算終了後の雑音による誤差の最大値の方が小さくなるので、本実施形態を適用することにより雑音の影響度が小さくなり、その結果、エッジ測定の再現性すなわち安定性が向上することになる。
【0027】
以上、本発明の一実施形態について説明した。上記実施形態では、非線形演算の例として二乗演算を紹介したが、本発明の趣旨は画像の微分値のダイナミックレンジを拡大して雑音成分の影響を小さくすることであるので、非線形演算の方法は二乗演算に限定するものではない。例えば、式(5)に示すような演算テーブルを使用しても同じような効果が得られる。
【0028】
【数6】
【0029】
すなわち、式(5)に示す演算テーブルでは、各画素の濃度値の微分値fi が大きいほど定数Kの値を大きく設定するようにしている。
【0030】
【発明の効果】
本発明によれば、非線形演算により得られた分布曲線はその元になった微分画像の分布曲線に比して、分布曲線のピーク領域がより大きく、一方、分布曲線のすそ野領域はより小さくなったので、すそ野領域の変動の影響が相対的に小さくなり、その結果、エッジ部近傍の微分値分布が急峻である場合だけでなく、ピントが外れたりエッジ部近傍の微分値分布がなだらかな場合でも、撮影された濃淡画像に含まれるエッジ位置を、切り替え処理のない1つの画像処理シーケンスにより、正確かつ安定して検出することができるものとなった。
【図面の簡単な説明】
【図1】画像処理システムにおいて行われるエッジ検出の演算フローである。
【図2】本発明の一実施形態に係る画像処理システムを示すブロック図である。
【図3】濃淡画像メモリ6に記録された濃淡画像及び濃度値分布を示した図である。
【図4】微分フィルタ処理を施すことにより作成された微分画像及び微分値分布を示した図である。
【図5】非線形演算処理を施すことにより得られた分布を示した図である。
【図6】近似曲線とのマッチング処理を示した図である。
【符号の説明】
1 CCDカメラ
2 パソコン本体
3 マウス
4 キーボード
5 I/F
6 濃淡画像メモリ
7 表示制御部
8 CRT画面
9 I/F
10 I/F
11 CPU
12 プログラムメモリ
13 ワークメモリ
Claims (3)
- 被測定物を撮影することにより被測定物の濃淡画像データを得るCCDカメラと、該CCDカメラにて得られた前記濃淡画像データを格納する濃淡画像メモリと、該濃淡画像メモリに格納された前記濃淡画像データに対して演算処理を行うCPUと、を有する画像処理システムにおいて、
前記CPUでは以下の一連の演算処理、すなわち、
前記濃淡画像メモリに格納された前記濃淡画像データに対して微分処理を施すことにより微分画像を作成し、
該微分画像上を走査し、走査された画素の微分値の中から最も値の大きいものを抽出し、この最大の微分値を有する画素の位置をピーク位置とし、
該ピーク位置の近傍の微分画像に対して非線形演算を行い、
該非線形演算により得られた分布曲線に対して近似曲線とのマッチング処理を施すことにより被測定物のエッジ位置を算出する、一連の演算処理が行われることを特徴とする画像処理システム。 - 前記非線形演算は前記微分画像の微分値の二乗を演算することを特徴とする請求項1に記載の画像処理システム。
- 前記非線形演算は前記微分画像の微分値が大きいほど大きな定数を乗算することを特徴とする請求項1に記載の画像処理システム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2000190747A JP3728184B2 (ja) | 2000-06-26 | 2000-06-26 | 画像処理システム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2000190747A JP3728184B2 (ja) | 2000-06-26 | 2000-06-26 | 画像処理システム |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2002008047A JP2002008047A (ja) | 2002-01-11 |
JP3728184B2 true JP3728184B2 (ja) | 2005-12-21 |
Family
ID=18690162
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2000190747A Expired - Fee Related JP3728184B2 (ja) | 2000-06-26 | 2000-06-26 | 画像処理システム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3728184B2 (ja) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004038877A (ja) * | 2002-07-08 | 2004-02-05 | Yazaki Corp | 車両用周辺監視装置及び画像処理装置 |
JP5004120B2 (ja) * | 2007-03-20 | 2012-08-22 | 独立行政法人海上技術安全研究所 | 舶用波浪監視装置 |
-
2000
- 2000-06-26 JP JP2000190747A patent/JP3728184B2/ja not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP2002008047A (ja) | 2002-01-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4160258B2 (ja) | 勾配ベースの局部輪郭線検出のための新しい知覚的しきい値決定 | |
US7292365B2 (en) | Methods and systems for determining distribution mean level without histogram measurement | |
JP6448065B2 (ja) | 画像処理システム、画像処理方法及び画像処理プログラム | |
JP2005172559A (ja) | パネルの線欠陥検出方法及び装置 | |
JP3754408B2 (ja) | パターン評価装置、パターン評価方法およびプログラム | |
JP2005165387A (ja) | 画面のスジ欠陥検出方法及び装置並びに表示装置 | |
JP3659426B2 (ja) | エツジ検出方法及びエツジ検出装置 | |
JP3728184B2 (ja) | 画像処理システム | |
JPH06160067A (ja) | 回路パターンの寸法測定方法 | |
JPH05307609A (ja) | パターン特徴抽出方法 | |
US6693669B1 (en) | Method for reducing image blurring | |
JP2871601B2 (ja) | 文字列検出装置及び方法 | |
JP3127598B2 (ja) | 画像中の濃度変動構成画素抽出方法および濃度変動塊判定方法 | |
JP2004219291A (ja) | 画面の線欠陥検出方法及び装置 | |
JP3223385B2 (ja) | 濃淡画像のパターンマッチング装置 | |
JP2001209810A (ja) | エッジ検出方法 | |
JPH0624014B2 (ja) | 濃淡画像の処理方法 | |
JP2005071125A (ja) | 被写体検出装置、被写体検出方法、被写体データ選定プログラムおよび被写体位置検出プログラム | |
WO2024098188A1 (zh) | 图像修复模型的视觉解析方法和装置、电子设备 | |
JP2004127784A (ja) | 荷電粒子線装置 | |
JP2002183713A (ja) | 画像処理方法およびその装置 | |
JP2003346165A (ja) | ワーク境界の検出方法 | |
TWI450572B (zh) | 影像邊界掃描的電腦系統及方法 | |
JP2841866B2 (ja) | 画像処理装置および方法ならびにその応用装置 | |
JP2004062266A (ja) | 被撮像物のエッジ位置検出方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20041216 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20050104 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050131 |
|
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: 20050913 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20050930 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081007 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091007 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091007 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101007 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101007 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111007 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111007 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121007 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121007 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131007 Year of fee payment: 8 |
|
LAPS | Cancellation because of no payment of annual fees |