CN102937593A - Ceramic radome crack automatic detection method - Google Patents
Ceramic radome crack automatic detection method Download PDFInfo
- 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
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
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:
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:
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:
And
In formula, B is structural element,
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
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:
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:
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:
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:
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:
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
In formula, B is structural element,
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
Individual pixel will, according to the connected mode between pixel, be calculated respectively distance when measuring.
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)
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)
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 |
-
2012
- 2012-10-20 CN CN201210402918XA patent/CN102937593A/en active Pending
Patent Citations (1)
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)
Title |
---|
于彬: "基于数字图像处理的陶瓷瓶裂纹检测研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》, no. 05, 15 November 2007 (2007-11-15) * |
刘刚等: "《MATLAB数字图像处理》", 31 July 2007, article "形态学图像处理", pages: 270 * |
Cited By (11)
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 |