ITMO20060093A1 - Metodo per l'individuazione di aree di disuniformita' su una superficie - Google Patents

Metodo per l'individuazione di aree di disuniformita' su una superficie Download PDF

Info

Publication number
ITMO20060093A1
ITMO20060093A1 IT000093A ITMO20060093A ITMO20060093A1 IT MO20060093 A1 ITMO20060093 A1 IT MO20060093A1 IT 000093 A IT000093 A IT 000093A IT MO20060093 A ITMO20060093 A IT MO20060093A IT MO20060093 A1 ITMO20060093 A1 IT MO20060093A1
Authority
IT
Italy
Prior art keywords
dimensional
value
areas
elements
domain
Prior art date
Application number
IT000093A
Other languages
English (en)
Inventor
Davide Lampugnani
Franco Stefani
Original Assignee
System Spa
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 System Spa filed Critical System Spa
Priority to IT000093A priority Critical patent/ITMO20060093A1/it
Priority to DE602007001533T priority patent/DE602007001533D1/de
Priority to ES07705596T priority patent/ES2329194T3/es
Priority to CN2007800071842A priority patent/CN101395635B/zh
Priority to PCT/IB2007/000362 priority patent/WO2007107816A1/en
Priority to AT07705596T priority patent/ATE436059T1/de
Priority to EP07705596A priority patent/EP1997071B1/en
Publication of ITMO20060093A1 publication Critical patent/ITMO20060093A1/it

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0004Industrial image inspection
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/84Systems specially adapted for particular applications
    • G01N21/88Investigating the presence of flaws or contamination
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/90Determination of colour characteristics
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30108Industrial image inspection
    • G06T2207/30124Fabrics; Textile; Paper

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Biochemistry (AREA)
  • Immunology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Quality & Reliability (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Image Analysis (AREA)
  • Image Processing (AREA)
  • Investigating Materials By The Use Of Optical Means Adapted For Particular Applications (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
  • Inspection Of Paper Currency And Valuable Securities (AREA)

Description

DESCRIZIONE
annessa a domanda di brevetto per INVENZIONE INDUSTRIALE avente per titolo: METODO PER L'INDIVIDUAZIONE DI AREE DI DISUNIFORMITÀ SU UNA SUPERFICIE.
Oggetto della presente invenzione è un metodo per rindividuazione di aree di disuniformità su una superficie.
In particolare il metodo si rivolge all'individuazione di disuniformità cromatiche, intese come difetti o imperfezioni, presenti nel motivo ornamentale riprodotto su una superficie. Preferibilmente, ma non esclusivamente, tale superficie è la superficie di una piastrella ceramica. Generalmente i difetti in oggetto assumono la forma di macchie dal contorno e dalla colorazione più o meno regolari ed intense.
Sono noti metodi per l'individuazione di difetti sulla superficie decorata di una piastrella ceramica che consistono sommariamente nell'acquisizione e nella digitalizzazione di un'immagine della superficie decorata della piastrella, nella misurazione, effettuata su tale immagine della piastrella, del valore di alcuni parametri che sono indicativi di determinate caratteristiche della superficie della piastrella, nel confronto di tali valori misurati con valori di soglia che sono indicativi di una deviazione tollerabile (o intollerabile) del corrispondente parametro da un valore corrispondente all'assenza di difetti.
Le immagini delle piastrelle vengono acquisite mediante una serie di telecamere almeno una delle quali è puntata perpendicolarmente alla superficie recante il decoro sul quale si vogliono individuare eventuali difetti.
La misurazione dei parametri indicativi ed il loro confronto con i valori di soglia vengono eseguiti in base ad algoritmi matematici tra loro differenti a seconda del parametro indicativo analizzato. Parametri indicativi diversi possono essere costituiti ad esempio da un valore medio di grigio, da un valore medio di contrasto o da una percentuale di somiglianza della superficie rispetto ad una immagine di riferimento prestabilita.
Quando la differenza tra il valore di un parametro indicativo ed il corrispondente valore di soglia supera un predeterminato livello di tolleranza, viene segnalata la potenziale presenza di un difetto. L’affidabilità di un metodo di individuazione è dato dalla capacità di indicare un superamento dei livelli di tolleranza predeterminati solo quando viene rilevato un reale difetto, ovvero quando il superamento dei livelli di tolleranza è relativamente netto.
Alcuni esempi di algoritmi matematici noti utilizzati per la misurazione dei parametri indicativi ed il loro confronto con i valori di soglia sono l'algoritmo “Contrasto”, l'algoritmo “Lum” e l'algoritmo “Macchie”.
L'algoritmo Contrasto si basa su una convoluzione deH’immagine ottenuta applicando diversi filtri. L’immagine risultante è un’immagine in cui tutto ciò che presenta un alto gradiente, ovvero un elevato contrasto rispetto allo sfondo, viene evidenziato con un maggiore livello di grigio mentre tutto il resto viene attenuato. Questo tipo di algoritmo è particolarmente efficace per la ricerca di difetti in decori realizzati su superfici lisce in tinta unita. In tali condizioni è possibile infatti effettuare un confronto relativamente accurato tra le misurazioni fatte sulla superficie analizzata ed un insieme di valori standard corrispondenti ad una decorazione priva di difetti.
L’algoritmo Lum si basa sulla misurazione della luminanza globale della superficie decorata. In sostanza l'algoritmo si basa sulla ricerca di pixel più scuri o più chiari rispetto alla media delle piastrelle senza difetti. Tali pixel sono individuabili mediante il confronto con livelli di tolleranza superiore ed inferiore della luminanza globale. L'algoritmo è efficacemente utilizzabile anche in presenza di decorazioni riprodotte su uno sfondo non a tinta unita. Al contrario l'algoritmo non è efficace in presenza di decorazioni casuali della superficie.
L’algoritmo di Macchie si basa su un'analisi locale di luminosità. A tal fine l'algoritmo effettua una divisione dell’ immagine in tanti tasselli, ciascuno dei quali viene identificato da un valore binario in base al valore medio di grigio più alto trovato tra i tasselli più vicini. Questo algoritmo permette di individuare agglomerati di pixel aventi valori di grigio che si discostano rispetto allo sfondo della piastrella. Va riscontrato che in presenza di piastrelle con decorazioni casuali l'accuratezza dei risultati fomiti dall'algoritmo è relativamente bassa a causa della difficoltà nella corretta impostazione dei valori di soglia corrispondenti al contrasto tra gli agglomerati di pixel individuati e lo sfondo della piastrella.
I metodi descritti, in presenza di superimi lucide di colore variabile o recanti decorazioni casuali, non consentono di individuare in maniera affidabile i difetti, soprattutto se essi si presentano in forma di macchie di colore molto simile al colore dalla superficie e con bordi sfumati. In tal caso i riconoscimenti erronei di presunti difetti aumentano notevolmente all'aumentare della severità dei valori di soglia impostati. Spesso è necessario accettare valori di compromesso non soddisfacenti che tra l'altro possono essere definiti solo empiricamente, con grande dispendio di tempo, attraverso un'analisi per campionamento di superimi decorate e modificando tutti i possibili parametri a disposizione.
Scopo della presente invenzione è proporre un metodo per l'individuazione di aree di disuniformità su una superficie che consenta di superare gli inconvenienti dei metodi di tipo noto. Un vantaggio del metodo è che non è necessaria alcuna impostazione empirica di parametri significativi e/o valori di soglia.
Un altro vantaggio del metodo è di essere efficace su qualunque tipo di superficie e per qualunque tipo di decorazione. Ulteriori caratteristiche e vantaggi dell’invenzione meglio appariranno dalla descrizione dettagliata che segue di un esempio di attuazione preferito, ma non esclusivo, del metodo in oggetto. Il metodo comprende le fasi principali che verranno di seguito elencate e descritte.
Una fase del metodo consiste nella definizione delle principali caratteristiche geometriche e dimensionali delle aree di disuniformità da individuare. In questa fase si tratta in sostanza di predefinire la forma dei difetti che possono interessare la superficie decorata da esaminare. La forma di tali difetti può, con ottima approssimazione, essere assunta come circolare o ellittica, in quanto in generale essa deriva dalla formazione di gocce di colorante che si depositano sulla superficie decorata. Le caratteristiche geometriche e dimensionali delle aree di disuniformità comprendono quindi il diametro e/o i semiassi delle aree stesse.
Il metodo comprende inoltre una fase di definizione di un modello matematico il quale è rappresentato da una finizione predeterminata che stabilisce una corrispondenza o correlazione tra elementi di un dominio o sottodominio bidimensionale, costituito da un'area avente dimensioni e caratteristiche geometriche (di forma) simili a quelle delle aree di disuniformità da individuare, e corrispondenti elementi di un condominio o sottocodominio bidimensionale ciascuno dei quali esprime un valore di intensità di tono di grigio.
Tale funzione predeterminata che stabilisce una corrispondenza o correlazione tra elementi di un dominio o sottodominio bidimensionale, costituito da un'area avente dimensioni e caratteristiche geometriche (forma) simili a quelle delle aree di disuniformità da individuare, e corrispondenti elementi di un condominio bidimensionale ciascuno dei quali esprime un valore di intensità di tono di grigio, è una funzione bidimensionale di Gauss che può essere espressa nel modo seguente:
in cui δ indica la deviazione standard (ovvero la distanza rispetto alla media in cui la funzione presenta il flesso) ed Xo ed y0rappresentano le coordinate del centro della goccia. Tutte le coordinate possono essere espresse in pixel con riferimento airimmagine della superficie che viene acquisita in una fase ulteriore del metodo.
Le coordinate del centro della goccia possono essere considerate (0,0) ed in tal caso l'equazione (1.1) diviene:
Il parametro Κ=1/δ(2π)<-1>nella (1.2) esprime il valore di picco della funzione stessa. Il parametro K può assumere un valore compreso tra 1 e 255. Si è trovato che il valore K=255 rappresenta il valore ottimale per la funzione in quanto, per tale valore, il metodo è sostanzialmente invariante al valore di grigio del difetto.
Quando si assume che l'area di disuniformità, o difetto, ha forma circolare, si applica la formula (1.2).
Se invece si assume che la forma del difetto è ellittica, allora, unitamente alla (1.2), valgono le seguenti condizioni di sistema:
dove D=2a è l'asse maggiore dell'ellisse e 2b è l'asse minore dell'ellisse.
Il metodo comprende altresì una fase di esplorazione della superficie da esaminare per ricavarne un'immagine. L'esplorazione della superficie può essere effettuata mediante una telecamera, a colori o in bianco e nero.
L'immagine ricavata viene confrontata, per zone o "finestre" di dimensioni prestabilite, con il condomino o sottocodominio rappresentante il modello matematico e vengono rilevate le differenze tra i due oggetti del confronto.
Per ogni zona o "finestra" si esegue il calcolo di un indice numerico complessivo r che esprime il valore di correlazione con il modello matematico. Tale calcolo può essere effettuato in forma di convoluzione del tipo seguente:
La (1.4) rappresenta una convoluzione non normalizzata, in cui Mi sono gli N pixel del modello ed li sono gli N pixel dell'immagine presenti nella finestra di confronto.
Poiché la correlazione viene effettuata per ogni pixel dell' immagine di ricerca, le posizioni della finestra dove il risultato acquista valori più alti sono quelle in cui è presente il modello di difetto da trovare. Per ottenere stabilità nelle misurazioni dell’algoritmo, cioè per ottenere invarianza rispetto al valore di grigio dell’ immagine e rispetto al modello, occorre modificare il calcolo di r normalizzandolo:
In questo caso il risultato r raggiunge il massimo, pari a 1, quando il modello ed il difetto sono uguali, diventa nullo quando non sono correlati.
Ciascun indice numerico complessivo viene quindi confrontato con un valore di soglia di accettabilità prestabilito. Se un indice numerico supera il valore di soglia di accettabilità prestabilito si può dedurre l'informazione che la corrispondente zona o finestra di confronto presenta un difetto.
Il metodo descritto consegue importanti vantaggi.
Innanzitutto poiché, si basa su un modello matematico, non è necessario effettuare operazioni empiriche di impostazione del modello, come invece accade per i metodi di tipo noto in cui è necessario impostare il modello sulla base di un certo numero di difetti effettivamente presenti su superfici campione.
Il metodo è sostanzialmente indipendente dal tipo di superficie e dal tipo di decorazione su essa riprodotta. Quando non sono presenti difetti che corrispondono alla caratteristiche di forma preimpostate l'indice numerico complessivo r è molto basso, mentre diventa significativo in caso contrario. In questo modo è possibile separare nettamente i difetti effettivi dai difetti erroneamente riconosciuti. Inoltre, impostando un campo di variazione corretto delle dimensioni dei difetti, è possibile individuare contemporaneamente diversi diffetti, incrementando notevolmente l'efficacia del metodo rispetto ai metodi di tipo noto.
Il metodo funziona correttamente sia in presenza di immagini a colori, sia in presenza di immagini in bianco e nero.
Ulteriore importante vantaggio è dato dalla semplicità di impostazione dei parametri, in quanto è sufficiente impostare il diametro minimo e massimo, e/o l'ovalità minima e massima, dei difetti per ottenere un'elevata efficacia nel riconoscimento dei difetti stessi.

Claims (5)

  1. RIVENDICAZIONI 1) Metodo per l'individuazione di aree di disuniformità su una superfìcie, caratterizzato dal fatto che comprende le seguenti fasi principali: - definizione delle principali caratteristiche geometriche e dimensionali delle aree di disuniformità da individuare; - definizione di un modello matematico il quale è rappresentato da una funzione predeterminata che stabilisce una corrispondenza o correlazione tra elementi di un dominio o sottodominio bidimensionale, costituito da un'area avente dimensioni e caratteristiche geometriche (di forma) simili a quelle delle aree di disuniformità da individuare, e corrispondenti elementi di un condominio o sottocodominio bidimensionale ciascuno dei quali esprime un valore di intensità di tono di grigio; - esplorazione della detta superficie per ricavarne un'immagine; - confronto, per zone o "finestre" di dimensioni prestabilite, della detta immagine con il condomino o sottocodominio rappresentante il detto modello matematico e rilevamento delle differenze; - calcolo, per ogni detta zona o "finestra", di un indice numerico complessivo esprimente il valore di correlazione con detti elementi di un condominio o sottocodominio bidimensionale ciascuno dei quali esprime un valore di intensità di tono di grigio; - confronto tra ciascun detto indice numerico complessivo con un valore di soglia di accettabilità prestabilito; - segnalazione dell'eventuale superamento del valore di soglia di accettabilità prestabilito.
  2. 2) Metodo secondo la rivendicazione 1, caratterizzato dal fatto che detta funzione predeterminata che stabilisce una corrispondenza o correlazione tra elementi di un dominio o sottodominio bidimensionale, costituito da un'area avente dimensioni e caratteristiche geometriche (forma) simili a quelle delle aree di disuniformità da individuare, e corrispondenti elementi di un condominio bidimensionale ciascuno dei quali esprime un valore di intensità di tono di grigio, è una funzione bidimensionale di Gauss:
  3. 3) Metodo secondo la rivendicazione 1, caratterizzato dal fatto che detta funzione predeterminata che stabilisce una corrispondenza o correlazione tra elementi di un dominio o sottodominio bidimensionale, costituito da un'area avente dimensioni e caratteristiche geometriche (forma) simili a quelle delle aree di disuniformità da individuare, e corrispondenti elementi di un condominio bidimensionale ciascuno dei quali esprime un valore di intensità di tono di grigio, è espressa dalle condizioni di sistema:
  4. 4) Metodo secondo la rivendicazione 1, in cui detto calcolo, per ogni detta zona o "finestra", di un indice numerico complessivo esprimente il valore di correlazione con detti elementi di un condominio o sottocodominio bidimensionale ciascuno dei quali esprime un valore di intensità di tono di grigio, è effettuato in base alla funzione normalizzata:
  5. 5) Metodo secondo le rivendicazioni 2 o 3, in cui il parametro K assume un valore tra 1 e 255.
IT000093A 2006-03-21 2006-03-21 Metodo per l'individuazione di aree di disuniformita' su una superficie ITMO20060093A1 (it)

Priority Applications (7)

Application Number Priority Date Filing Date Title
IT000093A ITMO20060093A1 (it) 2006-03-21 2006-03-21 Metodo per l'individuazione di aree di disuniformita' su una superficie
DE602007001533T DE602007001533D1 (de) 2006-03-21 2007-02-06 Verfahren zum identifizieren ungleichförmiger bereiche auf einer oberfläche
ES07705596T ES2329194T3 (es) 2006-03-21 2007-02-06 Metodo para identificar areas no uniformes de una superficie.
CN2007800071842A CN101395635B (zh) 2006-03-21 2007-02-06 用于识别表面上的不均匀区域的方法
PCT/IB2007/000362 WO2007107816A1 (en) 2006-03-21 2007-02-06 A method for identifying non-uniform areas on a surface
AT07705596T ATE436059T1 (de) 2006-03-21 2007-02-06 Verfahren zum identifizieren ungleichförmiger bereiche auf einer oberfläche
EP07705596A EP1997071B1 (en) 2006-03-21 2007-02-06 A method for identifying non-uniform areas on a surface

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
IT000093A ITMO20060093A1 (it) 2006-03-21 2006-03-21 Metodo per l'individuazione di aree di disuniformita' su una superficie

Publications (1)

Publication Number Publication Date
ITMO20060093A1 true ITMO20060093A1 (it) 2007-09-22

Family

ID=38016793

Family Applications (1)

Application Number Title Priority Date Filing Date
IT000093A ITMO20060093A1 (it) 2006-03-21 2006-03-21 Metodo per l'individuazione di aree di disuniformita' su una superficie

Country Status (7)

Country Link
EP (1) EP1997071B1 (it)
CN (1) CN101395635B (it)
AT (1) ATE436059T1 (it)
DE (1) DE602007001533D1 (it)
ES (1) ES2329194T3 (it)
IT (1) ITMO20060093A1 (it)
WO (1) WO2007107816A1 (it)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI489420B (zh) * 2012-02-29 2015-06-21 Chroma Ate Inc 表面花紋之檢測方法

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE3639636C2 (de) * 1986-11-20 1996-04-18 Robert Prof Dr Ing Massen Automatische Inspektion von Textilbahnen
US5481619A (en) * 1992-07-07 1996-01-02 Schwartz; Nira Inspection method using templates images, unique histogram analysis, and multiple gray levels correlated to addresses of template images
DE19638065A1 (de) * 1996-09-18 1998-03-19 Massen Machine Vision Systems Automatische Qualitätskontrolle von Fliesen
DE19708582A1 (de) * 1997-03-03 1998-09-10 Bauer Ernst & Sohn Gmbh Co Kg Qualitätskontrolle für Kunststeine
US20010048765A1 (en) * 2000-02-29 2001-12-06 Steven Yi Color characterization for inspection of a product having nonuniform color characteristics
WO2003031956A1 (en) * 2001-10-05 2003-04-17 Millennium Venture Holdings Ltd System and method for classifying workpieces according to tonal variations
CN1306456C (zh) * 2002-12-13 2007-03-21 佳能株式会社 图像处理方法和装置
US7751643B2 (en) * 2004-08-12 2010-07-06 Semiconductor Insights Inc. Method and apparatus for removing uneven brightness in an image

Also Published As

Publication number Publication date
EP1997071B1 (en) 2009-07-08
ES2329194T3 (es) 2009-11-23
ATE436059T1 (de) 2009-07-15
CN101395635B (zh) 2011-04-13
WO2007107816A1 (en) 2007-09-27
DE602007001533D1 (de) 2009-08-20
EP1997071A1 (en) 2008-12-03
CN101395635A (zh) 2009-03-25

Similar Documents

Publication Publication Date Title
CN107705288B (zh) 伪目标运动强干扰下的危险气体泄漏红外视频检测方法
KR101730787B1 (ko) 결함 검사 장치 및 생산 시스템
CN104504388B (zh) 一种路面裂缝识别和特征提取算法及系统
US8121400B2 (en) Method of comparing similarity of 3D visual objects
CN103886597B (zh) 一种基于边缘检测与拟合曲线聚类的圆检测方法
CN109214380B (zh) 车牌倾斜校正方法
KR20150021049A (ko) 검사대상 물체의 표면 검사 방법 및 장치
CN108550160B (zh) 基于光强模板的非均匀光条特征区域提取方法
ITMI20081296A1 (it) &#34;procedimento per filtrare immagini sar dal rumore di speckle e relativo dispositivo&#34;
KR102073229B1 (ko) 표면 결함 검출 장치 및 표면 결함 검출 방법
CN110530278B (zh) 利用多线结构光测量间隙面差的方法
CN117095009B (zh) 一种基于图像处理的pvc装饰板缺陷检测方法
CN106643555A (zh) 基于结构光三维测量系统的连接件识别方法
CN104063707B (zh) 基于人类视觉多尺度感知特性的彩色图像聚类分割方法
CN107036560B (zh) 光学玻璃精密磨削加工的表面粗糙度检测方法
CN106815602B (zh) 一种基于多层次特征描述的跑道fod图像探测方法及装置
CN107610174A (zh) 一种鲁棒的基于深度信息的平面检测方法及系统
ITMO20060093A1 (it) Metodo per l&#39;individuazione di aree di disuniformita&#39; su una superficie
CN108109138B (zh) 一种针对类镜面物体的高光区域自适应匀光的方法
JP5858524B2 (ja) 化粧材の塗装色検出方法
CN116051992A (zh) 一种裂缝识别方法
CN113870131A (zh) 多光谱激光雷达彩色点云的高光去除方法
JP7187830B2 (ja) 画像処理プログラム、画像処理装置、及び画像処理方法
CN109544547A (zh) 一种卡口式汽车灯泡灯头顶部外观缺陷检测方法
Chen et al. Shadow detection based on rgb color model