CN102937593A - Ceramic radome crack automatic detection method - Google Patents

Ceramic radome crack automatic detection method Download PDF

Info

Publication number
CN102937593A
CN102937593A CN201210402918XA CN201210402918A CN102937593A CN 102937593 A CN102937593 A CN 102937593A CN 201210402918X A CN201210402918X A CN 201210402918XA CN 201210402918 A CN201210402918 A CN 201210402918A CN 102937593 A CN102937593 A CN 102937593A
Authority
CN
China
Prior art keywords
image
crackle
pixel
gray
ceramic radome
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.)
Pending
Application number
CN201210402918XA
Other languages
Chinese (zh)
Inventor
赵玉刚
李业富
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shandong University of Technology
Original Assignee
Shandong University of Technology
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 Shandong University of Technology filed Critical Shandong University of Technology
Priority to CN201210402918XA priority Critical patent/CN102937593A/en
Publication of CN102937593A publication Critical patent/CN102937593A/en
Pending legal-status Critical Current

Links

Abstract

The invention relates to a ceramic radome crack automatic detection method. The method is characterized in that the detection is carried out with a transmission image processing method. The method comprises the specific steps that: (1) a radome transmission image is obtained; (2) the image is subjected to gray scale and contrast enhancement processing, such that a crack area is underlined; (3) a binarization treatment is carried out; (4) a median filtering noise removing treatment is carried out; (5) a crack framework abstraction treatment is carried out; and (6) a crack length statistic treatment is carried out. With the treatment process, ceramic radome tiny cracks can be detected. According to the invention, with the radome crack transmission image processing method, radome crack automatic detection is realized.

Description

Ceramic radome crackle automatic testing method
Technical field
The present invention relates to a kind of ceramic radome crackle automatic testing method, belong to technical field of nondestructive testing.
Background technology
Thin-walled ceramic radome production is by often there being crack defect, but in current thin-walled ceramic radome product detects, detection to crackle is just undertaken by the method for manual observation, and it is large that this measuring method detects error, can't guarantee the detection quality of product.
Summary of the invention
The purpose of this invention is to provide a kind ofly can overcome above-mentioned defect, detect accurately, ceramic radome crackle automatic testing method that efficiency is high.Its technical scheme is:
A kind of ceramic radome crackle automatic testing method, is characterized in that adopting the transmission image disposal route to detect, and concrete steps are: 1) the antenna house transmission image obtains; 2) gradation of image contrast enhancement processing, outstanding slit region; 3) binary conversion treatment; 4) medium filtering is removed noise processed; 5) the key extraction process of crackle; 6) each Crack length of statistical treatment, obtain the length value of each Crack.
Described ceramic radome crackle automatic testing method, in step 1), arrange light source in ceramic radome inside, at the outside ccd video camera that adopts of ceramic radome, automatically carries out image acquisition, obtains the gray scale transmission image of ceramic radome different parts.
Described ceramic radome crackle automatic testing method, step 2) in, for crack defect district and non-defect area contrast are increased, outstanding crack defect district, need to carry out contrast enhancement processing, and contrast enhancement processing is the gray-scale value according to each pixel in certain rule pointwise modification image, thereby change gradation of image dynamic range and contrast, if original image is f (x, y), the image after processing is g (x, y), contrast increases and can be expressed as:
g(x,y)=T[f(x,y)]
In formula, T means the variation of the image corresponding pixel points grey scale mapping after original image and processing.
Described ceramic radome crackle automatic testing method, in step 3), detection for the ease of crack defect, crackle gray level and non-defect area gray level strictly need to be distinguished, binary conversion treatment is exactly that the gray level image that contrast is strengthened arranges a gray threshold S, the gray-scale value of each pixel on image and gray threshold are compared, and image is converted into the bianry image that only contains two gray levels the most at last:
g ( x , y ) = 255 f ( x , y ) &GreaterEqual; S 0 f ( x , y ) < S
Wherein f (x, y) is original image, and g (x, y) is bianry image, and gray level " 255 " is non-defect area, and gray level " 0 " is the crack defect district, and gray threshold S adopts the iteration threshold method to try to achieve.
Described ceramic radome crackle automatic testing method, in step 4), in order to remove the noise in binary image, need to be removed noise processed, the medium filtering Denoising disposal is the nonlinear smoothing method that suppresses noise, effective impulse noise mitigation, good keep the edge information simultaneously, for given n number { a 1, a 2..., a n, by their ordered arrangement by size, when n is odd number, the number that is positioned at centre position is called the intermediate value of this n number, and when n is even number, the mean value that is positioned at centre position two numbers is called the intermediate value of this n number, is designated as med[a 1, a 2..., a n], to set the pixel neighborhood of a point in medium filtering, to this neighborhood, be all rank-ordered pixels in moving window, the gray-scale value by its gray scale intermediate value as processed pixel, array [x (i, j)] M * NThrough window, be A nMedium filtering after, the output of pixel (i, j) is designated as:
y ( i , j ) = med A n ( i , j ) [ x ( i , j ) ]
In formula, A n(i, j) means the neighborhood of point (i, j), and it contains n pixel.
Described ceramic radome crackle automatic testing method, in step 5), the key extraction of crackle is to carry out corrosion treatment on the Width of crackle, until crack width only has a pixel unit, it has kept the topological property of crackle, and the process that the set A skeleton is turned to S (A) can be expressed as:
S ( A ) = &cup; k = 0 K S k ( A )
And
In formula, B is structural element,
Figure BDA00002283387500025
Continuous k the corrosion of expression to A.K is the number of times that A is etched into the front last iteration of empty set.
Described ceramic radome crackle automatic testing method, in step 6), for every Crack crack length statistics wherein, be that the crackle that the backbone has been changed is traveled through since an end, after the whole crackle of traversal, the length of crackle has just come out, wherein the statistics of length comprises two kinds of situations, it between pixel, is the upper and lower, left and right four direction of mode arranged side by side, distance between this pixel arranged side by side is 1 pixel, vergence direction has upper left, lower-left, upper right, bottom right four direction, and the distance between vergence direction is like this
Figure BDA00002283387500026
Individual pixel will, according to the connected mode between pixel, be calculated respectively distance when measuring.
Compared with prior art, its advantage is in the present invention: adopted image to process automatic testing method the crackle of ceramic radome is detected, measurement result is more accurate, objective.
Embodiment
Inside at ceramic radome arranges light source, then, at the outside transmission image that adopts ccd video camera automatically to gather the antenna house different parts of antenna house, the gray level image gathered is processed according to following steps:
Step 1): for the crack defect district in the gray level image that makes to collect and the increase of non-defect area contrast and outstanding crack defect district, need to carry out contrast enhancement processing, adopt piecewise linearity to change and is processed, transformation for mula is:
g ( x , y ) = c a &times; f ( x , y ) [ 0 , a ) d - c b - a &times; ( f ( x , y ) - a ) + c [ a , b ] 255 - d 255 - b &times; ( f ( x , y ) - b ) + d ( b , 255 ]
Wherein f (x, y) means the gray level image collected, and g (x, y) means the image that contrast strengthens, and waypoint is elected as according to the grey value profile situation of the image obtained: a=70, b=170, c=30, d=110.
Step 2): step 1) is processed to the gray level image obtained a gray threshold is set, the gray-scale value of all pixels on image and gray threshold are compared, and when gray-scale value is more than or equal to gray threshold, this pixel grey scale becomes " 255 ", otherwise gray scale becomes " 0 ", process formula as follows:
g ( x , y ) = 255 f ( x , y ) &GreaterEqual; S 0 f ( x , y ) < S
Image before f (x, y) means to process, the bianry image after g (x, y) means to process, S means the gray threshold of setting, and wherein the selection of gray threshold S adopts process of iteration to obtain, and implementation step is as follows:
(1) obtain maximum gray scale and the minimum gradation value of image, be designated as respectively ZMax and ZMin, making initial threshold is S K=(Zmax+Zmin)/2;
(2) according to threshold value S KImage is divided into to the 1st family and the 2nd family, obtains respectively both average gray value Z1 and Z2.
(3) obtain new threshold value S K+1=(Z1+Z2)/2.
(4) if specify a minimal value ε, have | S K+1-S K|<ε, if the value of approaching meets the demands substantially, gained is threshold value, S K+1Be last iteration result, otherwise make S K=S K+1, re-execute top computation process, until meet the error requirements condition, wherein minimal value ε is taken as 10.
Step 3): to step 2) binary image obtained is removed noise processed, adopts median filter to carry out 4 filtering to image and removes noise.
Step 4): after step 3) is removed noise, in order accurately to obtain the length of crackle, first will carry out key extraction process, crackle is carried out to the key image that 5 corrosion treatments have obtained crackle afterwards.
Step 5): after the crackle backboneization, crackle is carried out to the length statistics, employing is from an end of crackle, traveling through successively the backbone of whole crackle, thereby obtain the length of crackle, is wherein that distance between pixel has two kinds of situations in length statistics, it between pixel, is the upper and lower, left and right four direction of mode arranged side by side, distance between this pixel arranged side by side is 1 pixel, and vergence direction has upper left, lower-left, upper right, bottom right four direction, and the distance between vergence direction is like this Individual pixel will be calculated respectively distance according to the connected mode between pixel when measuring, and the length of the crackle of measurement is 25.32.
Evidence: crackle physical length value is 24.83, and measuring length is 25.32, and its relative error is calculated as: &delta; = | 25.32 - 24.83 | 24.83 &times; 100 % = 1.97 % , Error is controlled in 5%.

Claims (7)

1. a ceramic radome crackle automatic testing method, is characterized in that adopting the transmission image disposal route to detect, and concrete steps are: 1) the antenna house transmission image obtains; 2) gradation of image contrast enhancement processing, outstanding slit region; 3) binary conversion treatment; 4) medium filtering is removed noise processed; 5) the key extraction process of crackle; 6) each Crack length of statistical treatment, obtain the length value of each Crack.
2. ceramic radome crackle automatic testing method as claimed in claim 1, it is characterized in that: in step 1), in ceramic radome inside, light source is set, automatically carry out image acquisition at the outside ccd video camera that adopts of ceramic radome, obtain the gray scale transmission image of the whole different parts of ceramic radome.
3. ceramic radome crackle automatic testing method as claimed in claim 1, it is characterized in that: step 2) in, the gray-scale value of each pixel in the pointwise modification image, thus gradation of image dynamic range and contrast changed, if image is f (x, y), the image after processing is g (x, y), contrast increases and can be expressed as g (x, y)=T[f (x, y)], wherein T means the variation of input picture and output image corresponding pixel points grey scale mapping.
4. ceramic radome crackle automatic testing method as claimed in claim 1, it is characterized in that: in step 3), gray threshold S of image setting that gray scale is strengthened, gray threshold S adopts the iteration threshold method to try to achieve, then allow gray-scale value and the gray threshold S of each pixel on image compare, image is converted into the bianry image that only contains two gray levels the most at last:
g ( x , y ) = 255 f ( x , y ) &GreaterEqual; S 0 f ( x , y ) < S
Wherein f (x, y) is original image, and g (x, y) is bianry image, and gray level " 255 " is non-defect area, and gray level " 0 " is the crack defect district.
5. ceramic radome crackle automatic testing method as claimed in claim 1, is characterized in that: in step 4), for given n number { a 1, a 2..., a n, by their ordered arrangement by size, when n is odd number, the number that is positioned at centre position is called the intermediate value of this n number, and when n is even number, the mean value that is positioned at centre position two numbers is called the intermediate value of this n number, is designated as med[a 1, a 2..., a n], to set the pixel neighborhood of a point in medium filtering, to this neighborhood, be all rank-ordered pixels in moving window, the gray-scale value by its gray scale intermediate value as processed pixel, array [x (i, j)] M * NThrough window, be A nMedium filtering after, the output of pixel (i, j) is designated as:
Figure FDA00002283387400012
In formula, A n(i, j) means the neighborhood of point (i, j), and it contains n pixel.
6. ceramic radome crackle automatic testing method as claimed in claim 1, it is characterized in that: in step 5), the key extraction of crackle is to carry out corrosion treatment on the Width of crackle, until crack width only has a pixel unit, it has kept the topological property of crackle, and the process that the set A skeleton is turned to S (A) can be expressed as: And
Figure FDA00002283387400022
In formula, B is structural element,
Figure FDA00002283387400024
Continuous k the corrosion of expression to A, k is the number of times that A is etched into the front last iteration of empty set.
7. ceramic radome crackle automatic testing method as claimed in claim 1, it is characterized in that: in step 6), for every Crack length statistics wherein, be that the crackle that the backbone has been changed is traveled through since an end, after traversal whole piece crackle, the length of crackle has just come out, wherein the statistics of length comprises two kinds of situations, it between pixel, is the upper and lower, left and right four direction of mode arranged side by side, distance between this pixel arranged side by side is 1 pixel, vergence direction has upper left, lower-left, upper right, bottom right four direction, and the distance between vergence direction is like this
Figure FDA00002283387400025
Individual pixel will, according to the connected mode between pixel, be calculated respectively distance when measuring.
CN201210402918XA 2012-10-20 2012-10-20 Ceramic radome crack automatic detection method Pending CN102937593A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210402918XA CN102937593A (en) 2012-10-20 2012-10-20 Ceramic radome crack automatic detection method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210402918XA CN102937593A (en) 2012-10-20 2012-10-20 Ceramic radome crack automatic detection method

Publications (1)

Publication Number Publication Date
CN102937593A true CN102937593A (en) 2013-02-20

Family

ID=47696504

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210402918XA Pending CN102937593A (en) 2012-10-20 2012-10-20 Ceramic radome crack automatic detection method

Country Status (1)

Country Link
CN (1) CN102937593A (en)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105335744A (en) * 2015-11-10 2016-02-17 佛山科学技术学院 One-dimensional code region location based on image backbone extraction strip distribution features
CN106124742A (en) * 2016-07-11 2016-11-16 无锡市华东电力设备有限公司 A kind of pigment and the method for inspection
CN106226157A (en) * 2016-08-31 2016-12-14 孙金更 Concrete structure member crevices automatic detection device and method
CN107389695A (en) * 2015-04-20 2017-11-24 赵媛媛 Detection method based on visible light communication technology porcelain crack detecting device
CN107993223A (en) * 2017-11-27 2018-05-04 歌尔股份有限公司 Scratch detection method, apparatus and electronic equipment
CN107991307A (en) * 2017-08-28 2018-05-04 中国人民解放军总后勤部油料研究所 A kind of soft material face crack automatic detection device and method
CN108020548A (en) * 2017-08-03 2018-05-11 刘素兰 Sinking degree measuring system
CN110009606A (en) * 2019-03-22 2019-07-12 北京航空航天大学 A kind of crack propagation dynamic monitoring method and device based on image recognition
CN113866700A (en) * 2021-10-11 2021-12-31 上海霍莱沃电子系统技术股份有限公司 Device and method for calibrating mechanical precision of antenna array surface test based on laser range finder

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11326225A (en) * 1998-05-11 1999-11-26 Hitachi Ltd Method and apparatus for inspecting tubular member

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11326225A (en) * 1998-05-11 1999-11-26 Hitachi Ltd Method and apparatus for inspecting tubular member

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
于彬: "基于数字图像处理的陶瓷瓶裂纹检测研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》, no. 05, 15 November 2007 (2007-11-15) *
刘刚等: "《MATLAB数字图像处理》", 31 July 2007, article "形态学图像处理", pages: 270 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107389695A (en) * 2015-04-20 2017-11-24 赵媛媛 Detection method based on visible light communication technology porcelain crack detecting device
CN105335744A (en) * 2015-11-10 2016-02-17 佛山科学技术学院 One-dimensional code region location based on image backbone extraction strip distribution features
CN105335744B (en) * 2015-11-10 2018-09-21 佛山科学技术学院 A kind of one-dimension code zone location extracting band distribution characteristics based on image backbone
CN106124742A (en) * 2016-07-11 2016-11-16 无锡市华东电力设备有限公司 A kind of pigment and the method for inspection
CN106226157A (en) * 2016-08-31 2016-12-14 孙金更 Concrete structure member crevices automatic detection device and method
CN108020548A (en) * 2017-08-03 2018-05-11 刘素兰 Sinking degree measuring system
CN107991307A (en) * 2017-08-28 2018-05-04 中国人民解放军总后勤部油料研究所 A kind of soft material face crack automatic detection device and method
CN107991307B (en) * 2017-08-28 2020-04-28 中国人民解放军总后勤部油料研究所 Automatic detection device and method for surface cracks of soft material
CN107993223A (en) * 2017-11-27 2018-05-04 歌尔股份有限公司 Scratch detection method, apparatus and electronic equipment
CN110009606A (en) * 2019-03-22 2019-07-12 北京航空航天大学 A kind of crack propagation dynamic monitoring method and device based on image recognition
CN113866700A (en) * 2021-10-11 2021-12-31 上海霍莱沃电子系统技术股份有限公司 Device and method for calibrating mechanical precision of antenna array surface test based on laser range finder

Similar Documents

Publication Publication Date Title
CN102937593A (en) Ceramic radome crack automatic detection method
CN111429403B (en) Automobile gear finished product defect detection method based on machine vision
CN111080582B (en) Method for detecting defects of inner and outer surfaces of workpiece
CN105719259A (en) Pavement crack image detection method
CN108333488B (en) Arc detection method based on fusion of ultraviolet, infrared and optical images
CN111274939B (en) Automatic extraction method for road pavement pothole damage based on monocular camera
CN106780514B (en) The calculation method of the area heavy rain Ji Lao depth of accumulated water based on monitor video image
CN105069778B (en) Based on the industrial products detection method of surface flaw that target signature notable figure builds
CN105608429A (en) Differential excitation-based robust lane line detection method
CN104239870A (en) Curve arc segmentation based ellipse detection method
CN115311507B (en) Building board classification method based on data processing
CN103778611A (en) Switch weighting vector median filter method utilizing edge detection
CN101527043A (en) Video picture segmentation method based on moving target outline information
CN114332081B (en) Textile surface abnormity determination method based on image processing
CN103116747B (en) Automatically the method and system of corn stem leaf image is identified
CN103455985B (en) Crack on road Enhancement Method based on Hessian structural analysis
CN104680488A (en) Method for performing satellite-borne SAR (Synthetic Aperture Radar) image ortho-rectification based on DSM (Digital Surface Model) of median filtering
CN104848917A (en) Image measuring method for analyzing water level in glass water tank
CN102509265B (en) Digital image denoising method based on gray value difference and local energy
CN103996177A (en) Snow noise removing algorithm free of reference detection
CN104598906B (en) Vehicle outline detection method and its device
CN108898080B (en) Ridge line neighborhood evaluation model-based crack connection method
CN103455986B (en) Random noise point detecting method based on fractional order differential gradient
CN113076959A (en) Concrete structure surface microcrack feature extraction method based on convolutional neural network
CN105405112A (en) Multispectral satellite image distance deviation index defogging method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20130220