TWI398632B - 粒子統計方法及裝置 - Google Patents

粒子統計方法及裝置 Download PDF

Info

Publication number
TWI398632B
TWI398632B TW97109212A TW97109212A TWI398632B TW I398632 B TWI398632 B TW I398632B TW 97109212 A TW97109212 A TW 97109212A TW 97109212 A TW97109212 A TW 97109212A TW I398632 B TWI398632 B TW I398632B
Authority
TW
Taiwan
Prior art keywords
particle
particles
image
shape factor
equal
Prior art date
Application number
TW97109212A
Other languages
English (en)
Other versions
TW200938832A (en
Inventor
Pei Chong Tang
Original Assignee
Hon Hai Prec Ind Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hon Hai Prec Ind Co Ltd filed Critical Hon Hai Prec Ind Co Ltd
Priority to TW97109212A priority Critical patent/TWI398632B/zh
Publication of TW200938832A publication Critical patent/TW200938832A/zh
Application granted granted Critical
Publication of TWI398632B publication Critical patent/TWI398632B/zh

Links

Description

粒子統計方法及裝置

本發明涉及一種圖像分析技術,特別是一種統計圖像中粒子數量的粒子統計方法及裝置。

在工業自動化光學檢測領域,統計圖像中粒子數量是很常見的應用,例如,與液晶產品製造有關的一種電路技術ACF(Anisotropic Conductive Film,各向異性導電薄膜)就需要根據電路中導電粒子的數量來判斷產品是否品質合格。

一般而言,統計圖像中目標數目的處理演算法是比較耗時的,因為不知道粒子的半徑或者長寬,而圖像中粒子又通常會粘連在一起,電腦很難識別哪個才是單個粒子,也就很難得到粒子的半徑或者長寬值,為此,目前的統計方法都需要經過複雜度較高的區域標記等處理過程,不僅處理時間長,而且容易受到雜訊的干擾影響統計精度。

鑒於此,有必要提供一種快速且準確的粒子統計方法。

還有必要提供一種快速且準確的粒子統計裝置。

一種粒子統計方法,包括如下步驟:獲取粒子圖像;至少進行一次粒子圖像的縮/放處理;計算縮/放前後至少兩幅粒子圖像中粒子的總面積;根據所述至少兩個粒子總面積以及對應的粒子圖像縮/放量計算出粒子圖像中粒子的數量,以判斷被拍攝物質中粒子含量是否滿足要求。

一種粒子統計裝置,包括:圖像獲取單元,用於獲取待統計的粒子圖像;縮/放單元,用於對粒子圖像進行縮/放處理;面積計算單元,用於計算縮/放前後粒子圖像中粒子的總面積;粒子統計單元,用於根據所述面積計算單元算得的至少兩個粒子總面積以及對應的粒子圖像縮/放量計算出粒子圖像中粒子的數量,以判斷被拍攝物質中粒子含量是否滿足要求。

藉由上述粒子統計方法及裝置計算粒子數量時,無需識別圖像中的獨立的粒子,也不需要計算粒子的直徑或者長寬等參數,僅利用以此圖像的縮/放就能夠得準確的粒子數量,使得統計速度大大提高。

以下實施方式提出一種快速的粒子數量統計方法和裝置,基本原理是藉由對二值化的粒子圖像進行圖像縮/放獲得粒子面積總和的統計值,然後根據粒子形狀特徵以及對應的運算公式計算出粒子數量,以判斷被拍攝物質中粒子含量是否滿足要求。。此技術可以對圓形、矩形或其他任意形狀較統一的粒子進行數量統計計算。下面以圓形粒子為例,詳細推導粒子統計的具體方法。

請參閱圖1,其為一張利用AOI(Automatic Optic Inspection,自動光學檢測)設備獲得的原始粒子圖像示意圖,圖中顏色較深的黑點為需要統計的粒子。

請參閱圖2,其為將圖1所示的原始粒子圖像進行二值化處理得到的二值粒子圖像,圖中白點表示粒子,有些白點面積較大且形狀不規則表示該處為多個粒子的粘連。

以下以一個實例對粒子數量統計方法作詳細說明,請參閱圖3,其為將圖2所示的二值粒子圖像進行收縮處理後的收縮粒子圖像,即從整個圖像或每個粒子的上下左右邊緣向中間壓縮,使整個圖像或每個粒子縮小。圖3中的圖像內容與圖2一致,只是其中的白點(粒子)比圖2中對應的白點要小。

首先,藉由電腦得到圖2所示的二值粒子圖像中粒子總面積S0,圖3所示的收縮粒子圖像中粒子總面積S1。

假設圖像中有M個粒子,粒子直徑為D(如圖2所示),圖3為收縮N個單位的收縮粒子圖像,那麼圖3中粒子的直徑為D-N(如圖3所示)。設X=D-N/2,則:S0=M *(π /4)D2 =M *(π /4)(X+N/2)2 S1=M *(π /4)(D-N)2 =M *(π /4)(X-N/2)

上述兩個方程中只有兩個未知數,那麼就可以得出M的值,本實施方式的計算方法如下:

通常縮/放單位N都遠小於粒子直徑為D,也就遠小於X,則:

由公式(1)和(2)可得:X=N *(S0+S1)/(S0-S1) D=X+N/2=N *(S1+S2)/(S1-S2)+N/2

最後可得到粒子數量M為:

其中:k=4/π ,即為圓形粒子的形狀係數,其他形狀的粒子只需改變k值即可。因為S0,S1,N均為已知數,則可以很快地得出粒子數量M。為了計算方便,運算過程中,所述直徑為D,收縮量N的單位都為圖元。

請參閱圖4、5、6,其分別為圓形粒子,正方形粒子以及長方形粒子的收縮示意圖。利用上述同樣的計算原理可以得出:正方形粒子的形狀係數k=1,長方形粒子的形狀係數k等於其長寬比,此時D為長方形寬度。

請參閱圖7,其為藉由上述粒子統計方法所測得的實驗資料,從資料可以看出,粒子數量統計誤差值保持在3%以內,這種精度的統計計算已經足以應付工業應用中許多需要粗略估計粒子數量應用的需求,如對ACF技術中導電粒子數量的檢測。

因為上述計算方法中存在約等於,從而可能會增加一定誤差,以下提供另一種計算方法:S0=M *(π /4)D2 ………(3)S1=M *(π /4)(D-N)2

那麼:

由公式(3)和(4)可得:

上述演算法沒有對任何資料進行省略,理論上統計結果更為準確。

為了簡化演算法,以下再提供一種計算方法:即對原始二值粒子圖像再進行一次放大處理,放大量也為N個單位,電腦算出放大粒子圖像中粒子的總面積S2。

S2=M *(π /4)(D+N)2

那麼:S2-S1=M *(π /4)* 4ND

實際統計中,因為粒子形狀並非完全的規則圖像,所以還需要加入適當的實驗修正值。

請參閱圖8,其為基於上述各種運算原理的粒子統計方法流程圖,包括如下步驟:步驟S803,獲取粒子圖像。

步驟S805,將粒子圖像作二值化處理。

步驟S807,至少進行一次粒子圖像的縮/放處理。

步驟S809,計算縮/放前後至少兩幅粒子圖像中粒子的總面積。

步驟S811,根據所述至少兩個粒子總面積以及對應的粒子圖像縮/放量計算出粒子圖像中粒子的數量,以判斷被拍攝物質中粒子含量是否滿足要求。。因為粒子的數量等於其總面積除以單個粒子的面積,則可以列出兩個面積(已知數)與粒子數量(未知數)和粒子半徑或邊長(未知數)的關係式,兩個關係式兩個未知數即可求得到兩個未知數的值,即可以得到粒子的數量。多種求解運算方法如下:若,縮/放前粒子的總面積等於S0;縮/放後粒子的總面積等於S1;縮/放量為N(放大時N為正值,縮小時N為負值),粒子的數量為M,則:

其中k為形狀係數,圓形粒子的形狀係數k=4/π ;正方形粒子的形狀係數k=1:長方形粒子的形狀係數k等於其長寬比。

若,縮/放前粒子的總面積等於S0;放大後粒子的總面積等於S1;縮小後粒子的總面積等於S2;縮/放量均為N,粒子的數量為M,則:

其中k為形狀係數:圓形粒子的形狀係數k=4/π ,正方形粒子的形狀係數k=1,長方形粒子的形狀係數k等於其長寬比。

步驟S813,輸出粒子數量,如藉由顯示器顯示資料或傳輸資料給其他設備。

藉由上述粒子統計方法計算粒子數量時,無需識別圖像中的獨立的粒子,也不需要計算粒子的直徑或者長寬等參數,僅利用以此圖像的縮/放就能夠得準確的粒子數量,使得統計速度大大提高。

請參閱圖9,其為基於上述運算原理的粒子統計裝置10的功能模組圖,包括:圖像獲取單元102、圖像處理單元104、縮/放單元106、面積計算單元108、粒子統計單元110以及輸出單元112。

圖像獲取單元102用於獲取待統計的粒子圖像,如藉由AOI設備拍攝的導電粒子圖像。

圖像處理單元104用於對粒子圖像進行二值化處理,得到二值粒子圖像。

縮/放單元106用於對二值粒子圖像進行收縮/放大處理。

面積計算單元108用於計算縮/放前後粒子圖像中粒子的總面積。

粒子統計單元110用於根據所述面積計算單元算得的至少兩個粒子總面積以及對應的粒子圖像縮/放量計算出粒子圖像中粒子的數量,以判斷被拍攝物質中粒子含量是否滿足要求。所述粒子的數量的計算是根據所述至少兩個面積與粒子數量和粒子半徑或邊長的關係式求得。

若,縮/放前粒子的總面積等於S0;縮/放後粒子的總面積等於S1;縮/放量為N(放大時N為正值,縮小時N為負值),粒子的數量為M,則:

其中k為形狀係數,圓形粒子的形狀係數k=4/π ;正方形粒子的形狀係數k=1;長方形粒子的形狀係數k等於其長寬比。

若,縮/放前粒子的總面積等於S0;放大後粒子的總面積等於S1;縮小後粒子的總面積等於S2;縮/放量均為N,粒子的數量為M,則:

其中k為形狀係數:圓形粒子的形狀係數k=4/π ,正方形粒子的形狀係數k=1,長方形粒子的形狀係數k等於其長寬比。

輸出單元112用於輸出粒子數量。

藉由上述粒子統計裝置計算粒子數量時,無需識別圖像中的獨立的粒子,也不需要計算粒子的直徑或者長寬等參數,僅利用以此圖像的縮/放就能夠得準確的粒子數量,使得統計速度大大提高。

本技術領域的普通技術人員應當認識到,以上的實施方式僅是用來說明本發明,而並非用作為對本發明的限定,只要在本發明的實質精神範圍之內,對以上實施例所作的適當改變和變化都落在本發明要求保護的範圍之內。

粒子統計裝置...10

圖像獲取單元...102

面積計算單元...108

粒子統計單元...110

圖像處理單元...104

縮/放單元...106

粒子統計方法步驟流程...S803~S813

輸出單元...112

圖1為原始粒子圖像示意圖。

圖2為二值化後之粒子圖像示意圖。

圖3為縮/放後之二值粒子圖像示意圖。

圖4為圓形粒子縮/放示意圖。

圖5為正方形粒子縮/放示意圖。

圖6為長方形粒子縮/放示意圖。

圖7為粒子統計結果資料圖。

圖8為一較佳實施方式之粒子統計方法步驟流程圖。

圖9為一較佳實施方式之粒子統計裝置功能模組圖。

粒子統計方法步驟流程...S803~S813

Claims (10)

  1. 一種粒子統計方法,包括如下步驟:獲取粒子圖像;至少進行一次粒子圖像的縮/放處理;計算縮/放前後至少兩幅粒子圖像中粒子的總面積;根據該至少兩個粒子總面積以及對應的粒子圖像縮/放量計算出粒子圖像中粒子的數量,以判斷被拍攝物質中粒子含量是否滿足要求。
  2. 如申請專利範圍第1項所述之粒子統計方法,其中該粒子的數量的計算是根據該至少兩個面積與粒子數量和粒子半徑或邊長的關係式求得。
  3. 如申請專利範圍第1項所述之粒子統計方法,其中若縮/放前粒子的總面積等於S0;縮/放後粒子的總面積等於S1;縮/放量為N,粒子的數量為M,則: 其中k為形狀係數,圓形粒子的形狀係數k=4/π ;正方形粒子的形狀係數k=1;長方形粒子的形狀係數k等於其長寬比。
  4. 如申請專利範圍第1項所述之粒子統計方法,其中若縮/放前粒子的總面積等於S0;放大後粒子的總面積等於S1;縮小後粒子的總面積等於S2;縮/放量均為N,粒子的數量為M,則: 其中k為形狀係數:圓形粒子的形狀係數k=4/π ,正方形粒子的形狀係數k=1,長方形粒子的形狀係數k等於其長寬比。
  5. 如申請專利範圍第1項所述之粒子統計方法,其中獲取粒子圖像後還包括將粒子圖像作二值化處理的步驟。
  6. 一種粒子統計裝置,包括:圖像獲取單元,用於獲取待統計的粒子圖像,其改良在於:該粒子統計裝置還包括:縮/放單元,用於對粒子圖像進行縮/放處理;面積計算單元,用於計算縮/放前後粒子圖像中粒子的總面積;粒子統計單元,用於根據該面積計算單元算得的至少兩個粒子總面積以及對應的粒子圖像縮/放量計算出粒子圖像中粒子的數量。
  7. 如申請專利範圍第6項所述之粒子統計裝置,其中該粒子的數量的計算是根據該至少兩個面積與粒子數量和粒子半徑或邊長的關係式求得。
  8. 如申請專利範圍第6項所述之粒子統計裝置,其中若縮/放前粒子的總面積等於S0;縮/放後粒子的總面積等於S1;縮/放量為N,粒子的數量為M,則: 其中k為形狀係數,圓形粒子的形狀係數k=4/π ;正方形粒子的形狀係數k=1;長方形粒子的形狀係數k等於其長寬比。
  9. 如申請專利範圍第6項所述之粒子統計裝置,其中若縮/放前粒子的總面積等於S0;放大後粒子的總面積等於S1;縮小後粒子的總面積等於S2;縮/放量均為N,粒子的數量為M,則: 其中k為形狀係數:圓形粒子的形狀係數k=4/π ,正方形粒子的形狀係數k=1,長方形粒子的形狀係數k等於其長寬比。
  10. 如申請專利範圍第6項所述之粒子統計裝置,其中該粒子統計裝置還包括圖像處理單元,用於對粒子圖像進行二值化處理。
TW97109212A 2008-03-14 2008-03-14 粒子統計方法及裝置 TWI398632B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
TW97109212A TWI398632B (zh) 2008-03-14 2008-03-14 粒子統計方法及裝置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
TW97109212A TWI398632B (zh) 2008-03-14 2008-03-14 粒子統計方法及裝置

Publications (2)

Publication Number Publication Date
TW200938832A TW200938832A (en) 2009-09-16
TWI398632B true TWI398632B (zh) 2013-06-11

Family

ID=44867465

Family Applications (1)

Application Number Title Priority Date Filing Date
TW97109212A TWI398632B (zh) 2008-03-14 2008-03-14 粒子統計方法及裝置

Country Status (1)

Country Link
TW (1) TWI398632B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS62239044A (en) * 1986-04-11 1987-10-19 Jeol Ltd Analyzing method for particle image in image
US20050099626A1 (en) * 2003-11-10 2005-05-12 King Frederick D. Method and apparatus for particle measurement employing optical imaging
TW200517647A (en) * 2003-11-28 2005-06-01 Samsung Corning Prec Glass Co Method for measuring particles in glass substrate

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS62239044A (en) * 1986-04-11 1987-10-19 Jeol Ltd Analyzing method for particle image in image
US20050099626A1 (en) * 2003-11-10 2005-05-12 King Frederick D. Method and apparatus for particle measurement employing optical imaging
TW200517647A (en) * 2003-11-28 2005-06-01 Samsung Corning Prec Glass Co Method for measuring particles in glass substrate

Also Published As

Publication number Publication date
TW200938832A (en) 2009-09-16

Similar Documents

Publication Publication Date Title
CA2847825C (en) Measurement of belt wear through edge detection of a raster image
US6549647B1 (en) Inspection system with vibration resistant video capture
US8184848B2 (en) Liquid level detection method
EP1664749B1 (en) Apparatus and method for automated web inspection
Zhu et al. Machine vision-based concrete surface quality assessment
Chauhan et al. Detection of bare PCB defects by image subtraction method using machine vision
JP2007024737A (ja) 半導体の欠陥検査装置及びその方法
US20130128026A1 (en) Image processing device for defect inspection and image processing method for defect inspection
JP2742240B2 (ja) 構造物表面の検査における欠陥の検出方法
JP3948728B2 (ja) パターン検査装置
TWI557650B (zh) 用於識別材料中之缺陷之系統及方法
JP5178079B2 (ja) 欠陥検査方法およびその装置
JP2017515097A (ja) 射影画像を用いた自動インライン検査及び計測
CN107084666B (zh) 基于机器视觉的刹车片尺寸综合检测方法
JP2017125861A (ja) 基板の異物質検査方法
CN104539937A (zh) 一种摄像头模组污点检测的方法和装置
JP2010117185A (ja) 欠陥検査装置および欠陥検査方法
CN104713885A (zh) 一种用于pcb板在线检测的结构光辅助双目测量方法
CN1223826C (zh) 影像测量系统和方法
CN103500459A (zh) 一种基于规则化参数的印刷电路板缺陷检测方法
JP4909691B2 (ja) 電子部品の実装状態検査方法、電子部品の実装状態検査装置及び電子機器の製造方法
US20140002673A1 (en) Line pair based full field sharpness test
TW201520540A (zh) 用於機器視覺檢驗的檢驗設備、方法和電腦程式產品
US20120093372A1 (en) Distance measuring device and distance measuring method
JPH0926802A (ja) ニューラルネットワークの学習方法

Legal Events

Date Code Title Description
MM4A Annulment or lapse of patent due to non-payment of fees