TWI398632B - 粒子統計方法及裝置 - Google Patents
粒子統計方法及裝置 Download PDFInfo
- 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
- shape factor
- image
- equal
- Prior art date
Links
Landscapes
- Image Analysis (AREA)
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項所述之粒子統計方法,其中該粒子的數量的計算是根據該至少兩個面積與粒子數量和粒子半徑或邊長的關係式求得。
- 如申請專利範圍第1項所述之粒子統計方法,其中若縮/放前粒子的總面積等於S0;縮/放後粒子的總面積等於S1;縮/放量為N,粒子的數量為M,則:
- 如申請專利範圍第1項所述之粒子統計方法,其中若縮/放前粒子的總面積等於S0;放大後粒子的總面積等於S1;縮小後粒子的總面積等於S2;縮/放量均為N,粒子的數量為M,則:
- 如申請專利範圍第1項所述之粒子統計方法,其中獲取粒子圖像後還包括將粒子圖像作二值化處理的步驟。
- 一種粒子統計裝置,包括:圖像獲取單元,用於獲取待統計的粒子圖像,其改良在於:該粒子統計裝置還包括:縮/放單元,用於對粒子圖像進行縮/放處理;面積計算單元,用於計算縮/放前後粒子圖像中粒子的總面積;粒子統計單元,用於根據該面積計算單元算得的至少兩個粒子總面積以及對應的粒子圖像縮/放量計算出粒子圖像中粒子的數量。
- 如申請專利範圍第6項所述之粒子統計裝置,其中該粒子的數量的計算是根據該至少兩個面積與粒子數量和粒子半徑或邊長的關係式求得。
- 如申請專利範圍第6項所述之粒子統計裝置,其中若縮/放前粒子的總面積等於S0;縮/放後粒子的總面積等於S1;縮/放量為N,粒子的數量為M,則:
- 如申請專利範圍第6項所述之粒子統計裝置,其中若縮/放前粒子的總面積等於S0;放大後粒子的總面積等於S1;縮小後粒子的總面積等於S2;縮/放量均為N,粒子的數量為M,則:
- 如申請專利範圍第6項所述之粒子統計裝置,其中該粒子統計裝置還包括圖像處理單元,用於對粒子圖像進行二值化處理。
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)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS62239044A (ja) * | 1986-04-11 | 1987-10-19 | Jeol Ltd | 画像中における粒子像分析方法 |
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 |
-
2008
- 2008-03-14 TW TW97109212A patent/TWI398632B/zh not_active IP Right Cessation
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS62239044A (ja) * | 1986-04-11 | 1987-10-19 | Jeol Ltd | 画像中における粒子像分析方法 |
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 |
---|---|---|
CN109544533B (zh) | 一种基于深度学习的金属板缺陷检测和度量方法 | |
JP6121425B2 (ja) | ラスタ画像のエッジ検出によるベルト摩耗の測定 | |
JP4425347B1 (ja) | 液体中の粒子のサイズの検出方法および装置 | |
CN104713885A (zh) | 一种用于pcb板在线检测的结构光辅助双目测量方法 | |
CN103673934A (zh) | 一种基于网格投影的pcb板平整度检测方法 | |
CN115797359A (zh) | 基于电路板上锡膏的检测方法、设备和存储介质 | |
CN112258507A (zh) | 互联网数据中心的目标对象检测方法、装置和电子设备 | |
CN103247548A (zh) | 一种晶圆缺陷检测装置及方法 | |
TWI398632B (zh) | 粒子統計方法及裝置 | |
JP3904796B2 (ja) | 異物または欠陥検査装置、および、異物または欠陥検査方法 | |
JP6478281B2 (ja) | 粗大粒子の測定方法および粗大粒子の測定装置 | |
CN112066888A (zh) | 输送机链条的磨损检测***、方法及装置 | |
US10062155B2 (en) | Apparatus and method for detecting defect of image having periodic pattern | |
CN109000560B (zh) | 基于三维相机检测包裹尺寸的方法、装置以及设备 | |
CN107734324B (zh) | 一种闪光灯照度均匀性的测量方法、***及终端设备 | |
JP7039899B2 (ja) | 識別装置、識別方法およびプログラム | |
TWI493177B (zh) | 一種檢測具週期性結構光學薄膜的瑕疵檢測方法及其檢測裝置 | |
CN115511718A (zh) | Pcb图像校正方法、装置、终端设备和存储介质 | |
RU2519005C1 (ru) | Способ допускового контроля печатных плат | |
WO2016188176A1 (zh) | 一种测量方法、装置及存储介质 | |
CN103837548B (zh) | 用于检测具有周期图案的图像的缺陷的装置和方法 | |
TW202400993A (zh) | Pcb缺陷檢測方法、電路板計量測量方法及製作方法 | |
JP2004294076A (ja) | Mtf測定装置及び測定方法 | |
JPH0794971B2 (ja) | 断面形状検知方法 | |
CN101520401B (zh) | 粒子统计方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
MM4A | Annulment or lapse of patent due to non-payment of fees |