CN110146027B - Method for measuring thickness of celadon glaze layer of SD-OCT image - Google Patents
Method for measuring thickness of celadon glaze layer of SD-OCT image Download PDFInfo
- Publication number
- CN110146027B CN110146027B CN201910491840.5A CN201910491840A CN110146027B CN 110146027 B CN110146027 B CN 110146027B CN 201910491840 A CN201910491840 A CN 201910491840A CN 110146027 B CN110146027 B CN 110146027B
- Authority
- CN
- China
- Prior art keywords
- glaze layer
- celadon
- pixel
- image
- celadon glaze
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000003708 edge detection Methods 0.000 claims abstract description 13
- 238000001914 filtration Methods 0.000 claims abstract description 5
- 238000012545 processing Methods 0.000 claims description 10
- 230000003287 optical effect Effects 0.000 claims description 6
- 230000010339 dilation Effects 0.000 claims description 4
- 230000003628 erosive effect Effects 0.000 claims description 4
- 238000012217 deletion Methods 0.000 claims description 3
- 230000037430 deletion Effects 0.000 claims description 3
- 238000012014 optical coherence tomography Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 13
- 229910052573 porcelain Inorganic materials 0.000 description 7
- 230000000694 effects Effects 0.000 description 6
- 239000000919 ceramic Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 238000000053 physical method Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/02—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness
- G01B11/06—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness for measuring thickness ; e.g. of sheet material
- G01B11/0616—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness for measuring thickness ; e.g. of sheet material of coating
- G01B11/0625—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness for measuring thickness ; e.g. of sheet material of coating with measurement of absorption or reflection
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
- G01B21/02—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring length, width, or thickness
- G01B21/04—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring length, width, or thickness by measuring coordinates of points
- G01B21/042—Calibration or calibration artifacts
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/41—Refractivity; Phase-affecting properties, e.g. optical path length
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
The invention discloses a method for measuring the thickness of a celadon glaze layer of an SD-OCT image. Establishing a database of refractive indexes of different types of celadon glaze layers by measuring the thickness of the sample celadon glaze layer; collecting an SD-OCT image of the celadon glaze layer; positioning the upper boundary of the glaze layer by filtering and binaryzation of the SD-OCT image of the celadon glaze layer; designing structural elements by edge detection, performing closed operation on the image subjected to edge detection by using the structural elements, and extracting the lower boundary of the celadon glaze layer; and the thickness of the celadon glaze layer is calculated by multiplying the difference value of the upper boundary pixel and the lower boundary pixel of the glaze layer by the physical depth of each pixel. The method realizes the nondestructive real-time measurement of the thickness of the celadon glaze layer, the measurement precision reaches micron level, the accuracy is high, the thickness of the glaze layer of various celadons can be measured according to the established celadon glaze layer refraction rate database, the method has strong adaptability, and the measurement efficiency is improved.
Description
Technical Field
The invention belongs to the field of automatic measurement of the thickness of a celadon glaze layer, relates to a processing method of an SD-OCT image, and particularly relates to a method for measuring the thickness of the celadon glaze layer of the SD-OCT image.
Background
China has a long history of ceramics, and celadon is one of the main types of ceramics, and is widely concerned and researched by people at home and abroad. The thickness of the glaze layer and the uniformity of the whole glaze layer of the porcelain are one of main factors for judging the quality of the celadon, and the traditional method for measuring the thickness of the celadon glaze layer is used for destroying the porcelain to carry out physical measurement, has certain precision, but cannot measure the thickness of the glaze layer of all parts of the porcelain and is a method for destructive measurement, so that the method for realizing nondestructive measurement has great significance.
The method of nondestructive measurement is not invented so far, the spectral domain optical coherence tomography (SD-OCT) shows the internal structure morphology and distribution of the porcelain by measuring the optical reflection scattering property of the substance, the SD-OCT image of the porcelain is used for the structural research, classification and qualitative identification of the porcelain at present, and reports show that the image can clearly show the internal structure of the porcelain. The SD-OCT imaging technology is used for measuring the thickness of the celadon glaze layer, and the measurement which is lossless, real-time and not limited by a measuring position has important research value.
Disclosure of Invention
Aiming at the problems in the background art, the invention aims to provide a method for measuring the thickness of the celadon glaze layer of an SD-OCT image, which can detect the thickness of the celadon glaze layer in the SD-OCT image in a real-time and nondestructive manner, and lays a technical foundation for the detection of the glazing uniformity of the celadon by matching with a method for manufacturing the celadon.
As shown in fig. 1, the technical scheme adopted by the invention comprises the following frame steps:
1) establishing a database of refractive indexes of different types of celadon glaze layers;
2) collecting an SD-OCT image of the celadon glaze layer;
3) extracting the upper boundary of the celadon glaze layer;
4) extracting the lower boundary of the celadon glaze layer;
5) and calculating the thickness of the celadon glaze layer.
The step 1) is specifically as follows:
1.1) collecting celadon samples of different types and with known glaze layer thickness;
1.2) acquiring an SD-OCT image of the celadon glaze layer by using an OQ LabScope system, taking a 'Depth Per Pixel' parameter in the OQ LabScope system as the physical resolution Pr of each Pixel in the SD-OCT image, measuring the thickness of the celadon glaze layer by using a system self-provided caliper tool, adjusting the 'Depth Per Pixel' parameter, and further adjusting the physical resolution Pr until the thickness of the celadon glaze layer measured by the system self-provided caliper tool is the same as the actual thickness of the celadon glaze layer;
1.3) recording the physical resolution Pr of each pixel of the SD-OCT image of the celadon glaze layer, and calculating the refractive index n of each celadon glaze layer according to the optical resolution Or and the physical resolution Pr as follows:
the optical resolution Or is obtained by processing of an optical acquisition system.
1.4) sequentially scanning all types of celadon glaze layers, calculating the refractive indexes of all types of celadon glaze layers according to the step 1.3) and establishing a celadon glaze layer refractive index database.
The step 2) is specifically as follows: and (3) acquiring an SD-OCT (secure digital-optical coherence tomography) image of the celadon with the glaze layer thickness to be measured by using an OQ LabScope system, and performing median filtering on the image by using a 3 multiplied by 5 size template.
The step 3) is specifically as follows:
3.1) carrying out binarization processing on the filtered image;
3.2) searching each row of pixel points in the binarized image from top to bottom to present a pixel point with a first gray value of 1, recording the point as a pixel point to be fitted, and recording the number of rows Top (i) of the pixel point to be fitted in each row, wherein i represents the number of rows of the pixel point to be fitted;
3.3) recording the median of the number of rows top (i) where the to-be-fitted pixel points of all columns are located as m, and starting to check and process the value of the number of rows top (i) where the to-be-fitted pixel points of each column are located from i ═ 1: if the absolute value of the difference between the number of rows Top (1) where the pixel point to be fitted in the first row is located and the median m is larger than 30, updating and replacing the number of rows Top (1) with the median between the number of rows Top (2) where the pixel point to be fitted in the second row is located and the number of rows Top (I) where the pixel point to be fitted in the last row is located, wherein I represents the number of columns in the last row; if the absolute value of the difference between the number of rows Top (2) where the pixel points to be fitted in the second row are located and the median m is greater than 30, updating the value of the number of rows Top (2) into the number of rows Top (1); if i is greater than 2 and the absolute value of the difference between the row number Top (i) of the pixel point to be fitted in the ith column and the row number Top (i-1) of the pixel point to be fitted in the i-1 column is greater than 30, updating the value of the row number Top (i) by using a Lagrange interpolation method, and calculating as follows:
Top(i)=2×Top(i-1)-Top(i-2)
wherein i represents the row ordinal number of the pixel point to be fitted;
3.4) after updating the number of rows Top (i) where the pixel points to be fitted of all columns are located, forming an upper boundary pixel point by all the updated pixel points to be fitted, and calculating to obtain the mean value Top of the number of rows Top (i) where the pixel points to be fitted of all columns are locatedAverage。
The step 4) is specifically as follows:
4.1) carrying out edge detection on the filtered image by using a canny operator, and taking the generated image as an edge detection image Fb;
4.2) defining a circular structural element se with the radius of five pixels, and performing a closing operation on the edge detection image Fb by using the structural element se to obtain an image Fc, wherein the closing operation is calculated as follows:
Fc=Fb·se
wherein,symbolRepresenting dilation operations on images, symbolsRepresenting an erosion operation on the image; the effect of the dilation operation is to "grow" or "coarsen" objects in the binary image, and the effect of the erosion operation is to "shrink" or "refine" objects in the binary image, both morphological operations of the image.
4.3) taking connected pixel points as a target in the image Fc, wherein the connection means that eight neighborhoods of the Sudoku are connected around the pixel points; then counting the number of pixel points of the target, arranging all targets from small to large according to the number of the pixel points, taking the number of the pixel points positioned in the median as a point threshold, deleting the targets with the number of the pixel points smaller than the point threshold, and then dividing the targets into two types: the first-class target is a target with the number of lines occupied by the contained pixel points larger than the number of columns, and is reserved; the other type of target is a target that the number of lines occupied by the contained pixel points is less than the number of columns, and deletion is carried out;
4.4) counting the number of rows of the most-lower pixel points at the row positions in all the targets remained in the image Fc, and solving the average value Bot of the row positions of all the pixel pointsAverage;
4.5) in the SD-OCT image of the celadon glaze layer, the D pixel distance between each column of the upper glaze layer boundary and the lower glaze layer boundary is calculated as follows:
D=BotAverage-TopAverage
wherein, TopAverageRepresenting the mean value of the number of rows top (i) where the pixel points to be fitted of all the columns after the updating is finished;
4.6) translating the upper boundary pixel point obtained in the step 3) downwards by a distance of D pixels according to a row to serve as the position of a lower boundary pixel point of the celadon glaze layer, and recording the position as bot (i).
The step 5) is specifically as follows:
5.1) obtaining the refractive index n corresponding to the celadon according to the type x of the measured celadonxThe physical resolution Pr of each pixel of the collected celadon SD-OCT image is obtained by the processing of the step 1)x;
5.2) converting the distance between the upper boundary and the lower boundary of the glaze layer on the celadon glaze layer SD-OCT image obtained in the step 4) into the actual glaze layer thickness T, and calculating as follows:
T=D×Prx
wherein D represents the distance between the upper and lower boundaries of the glaze layer on the celadon glaze layer SD-OCT image.
The method comprises the steps of calculating the refractive indexes of different kinds of celadons by measuring the thickness of a sample celadon glaze layer, establishing a database of the refractive indexes of different celadon glaze layers, filtering and binarizing an SD-OCT image of the celadon glaze layer, positioning the upper boundary of the glaze layer, carrying out edge detection on the filtered image, designing a structural element, carrying out closed operation on the image after edge detection by using the structural element, marking all targets in the image after closed operation, deleting the targets which do not meet the requirements, counting the positions of all pixels of the targets which meet the requirements in the image, calculating the longitudinal mean value of the minimum longitudinal position pixel in all the targets, taking the longitudinal mean value as the mean value of the longitudinal positions of the lower boundary pixels of the glaze layer, calculating the mean value of the longitudinal positions of all the upper boundary pixels of the glaze layer, and multiplying the difference value of the longitudinal position mean value of the lower boundary pixels of the glaze layer by the physical depth of each pixel to obtain the thickness of the glaze layer And (4) degree. The key points are the calculated refractive index of the celadon glaze layer, a software algorithm and the logic of the whole process.
The method realizes the nondestructive real-time measurement of the thickness of the celadon glaze layer, the measurement precision reaches micron level, the accuracy is high, the thickness of the glaze layer of various celadons can be measured according to the established celadon glaze layer refraction rate database, the method has strong adaptability, and the measurement efficiency is improved.
The invention has the beneficial effects that:
the method can be used for nondestructively measuring the thickness of the glaze layer at any position of the celadon glaze layer in real time, and the problem that the thickness of the glaze layer can be detected only by destroying the celadon is solved.
The method accurately positions the upper boundary and the lower boundary of the celadon glaze layer in the image, ensures the accuracy and the robustness of measurement, establishes the refraction rate databases of different celadon glaze layers, can measure the thickness of various types of celadon glaze layers, and improves the measurement universality.
Drawings
FIG. 1 is a flow chart of the method of the present invention.
FIG. 2 is a sample drawing of a celadon sample with a glaze layer to be measured, wherein (a) - (d) correspond to the samples respectively.
Fig. 3 is an SD-OCT image of the collected celadon glaze layer, and (a) - (d) correspond to (a) - (d) in fig. 2 in sequence.
FIG. 4 shows the SD-OCT image after noise reduction according to the embodiment, wherein (a) - (d) correspond to (a) - (d) in FIG. 3 in sequence.
Fig. 5 is a graph showing the effect of extracting the upper boundary of the celadon glaze layer according to the embodiment, wherein (a) - (d) correspond to (a) - (d) in fig. 2 in sequence.
Fig. 6 is a drawing showing the effect of extraction of the lower boundary of the celadon glaze layer according to the embodiment, wherein (a) - (d) correspond to (a) - (d) in fig. 2 in sequence.
Detailed Description
The present invention will be described in further detail below with reference to the accompanying drawings and examples. The specific embodiments described herein are merely illustrative of the invention and do not delimit the invention.
The specific embodiment of the invention is as follows:
1. a method for measuring the thickness of a celadon glaze layer of an SD-OCT image is characterized by comprising the following steps
1) Establishing a database of refractive indexes of different types of celadon glaze layers;
1.1) collecting celadon samples of different types and with known glaze layer thickness;
1.2) an OQ LabScope system produced by Lumedica company is used for collecting an SD-OCT image of the celadon glaze layer, a parameter of 'Depth Per Pixel' in the OQ LabScope system is used as the physical resolution Pr of each Pixel in the SD-OCT image, a system self-provided caliper tool is used for measuring the thickness of the celadon glaze layer to adjust the parameter of 'Depth Per Pixel', and then the physical resolution Pr is adjusted until the thickness of the celadon glaze layer measured by the system self-provided caliper tool is the same as the actual thickness of the celadon glaze layer.
1.3) recording the physical resolution Pr of each pixel of the SD-OCT image of the celadon glaze layer, and calculating the refractive index n of each celadon glaze layer according to the optical resolution Or and the physical resolution Pr.
1.4) sequentially scanning all types of celadon glaze layers, calculating the refractive indexes of all types of celadon glaze layers according to the step 1.3) and establishing a celadon glaze layer refractive index database.
2) And collecting an SD-OCT image of the celadon glaze layer. An OQ LabScope system produced by Lumedica is used for collecting SD-OCT images of celadon with the thickness of a glaze layer to be measured, and a 3 x 5 size template is used for carrying out median filtering on the images.
Fig. 2 is a sample drawing of a celadon sample with a glaze layer to be measured, wherein (a) - (d) correspond to the sample, fig. 3 is a collected SD-OCT image of the celadon glaze layer, wherein (a) - (d) correspond to (a) - (d) in fig. 2 in sequence, fig. 4 is a SD-OCT image after noise reduction, and wherein (a) - (d) correspond to (a) - (d) in fig. 3 in sequence.
3) Extracting the upper boundary of the celadon glaze layer;
3.1) carrying out binarization processing on the filtered image, wherein the threshold value of the binarization of the image is set to be 0.1;
3.2) searching each row of pixel points in the binarized image from top to bottom to present a pixel point with a first gray value of 1, recording the point as a pixel point to be fitted, and recording the number of rows Top (i) of the pixel point to be fitted in each row, wherein i represents the number of rows of the pixel point to be fitted;
3.3) recording the median of the number of rows top (i) where the to-be-fitted pixel points of all columns are located as m, and starting to check and process the value of the number of rows top (i) where the to-be-fitted pixel points of each column are located from i ═ 1:
if the absolute value of the difference between the number of rows Top (1) where the pixel point to be fitted in the first row is located and the median m is larger than 30, updating and replacing the number of rows Top (1) with the median between the number of rows Top (2) where the pixel point to be fitted in the second row is located and the number of rows Top (I) where the pixel point to be fitted in the last row is located, wherein I represents the number of columns in the last row; if not, not processing;
if the absolute value of the difference between the number of rows Top (2) where the pixel points to be fitted in the second row are located and the median m is greater than 30, updating the value of the number of rows Top (2) into the number of rows Top (1); if not, not processing;
and if i is greater than 2 and the absolute value of the difference between the row number Top (i) of the pixel point to be fitted in the ith column and the row number Top (i-1) of the pixel point to be fitted in the (i-1) th column is greater than 30, updating the value of the row number Top (i) by using a Lagrange interpolation method. If not, not processing;
3.4) after updating the number of rows Top (i) where the pixel points to be fitted of all columns are located, forming an upper boundary pixel point by all the updated pixel points to be fitted, and calculating to obtain the mean value Top of the number of rows Top (i) where the pixel points to be fitted of all columns are locatedAverage。
In this embodiment, fig. 5 is a graph showing the fitting effect of the upper boundary of the celadon glaze layer, where (a) - (d) correspond to (a) - (d) in fig. 2, respectively.
4) Extracting the lower boundary of the celadon glaze layer;
4.1) carrying out edge detection on the filtered image by using a canny operator, and taking the generated image as an edge detection image Fb;
4.2) defining a circular structure element se with the radius of five pixel lengths, wherein the pixel point of a circular area in the circular structure element se is 1, the surrounding pixel points are 0, and performing a closing operation on the edge detection image Fb by using the structure element se to obtain an image Fc, wherein the calculation is as follows:
Fc=Fb·se
wherein,symbolRepresenting dilation operations on images, symbolsRepresenting an erosion operation on the image;
4.3) taking connected pixel points as a target in the image Fc, wherein the connection means that eight neighborhoods of the Sudoku are connected around the pixel points; then counting the number of pixel points of the target, arranging all targets from small to large according to the number of the pixel points, taking the number of the pixel points positioned in the median as a point threshold, deleting the targets with the number of the pixel points smaller than the point threshold, and then dividing the targets into two types:
the first-class target is a target with the number of lines occupied by the contained pixel points larger than the number of columns, and is reserved;
the other type of target is a target that the number of lines occupied by the contained pixel points is less than the number of columns, and deletion is carried out;
4.4) counting the number of rows of the most-lower pixel points at the row positions in all the targets remained in the image Fc, and solving the average value Bot of the row positions of all the pixel pointsAverage;
4.5) in the SD-OCT image of the celadon glaze layer, the distance of D pixels between each column of the upper boundary and the lower boundary of the glaze layer is calculated.
4.6) translating the upper boundary pixel point obtained in the step 3) downwards by a distance of D pixels according to a row to serve as the position of a lower boundary pixel point of the celadon glaze layer, and recording the position as bot (i).
In this embodiment, fig. 6 is a fitting effect diagram of the lower boundary of the celadon glaze layer, where (a) - (d) correspond to (a) - (d) in fig. 2, respectively.
5) And calculating the thickness of the celadon glaze layer.
5.1) obtaining the refractive index n corresponding to the celadon according to the type x of the measured celadonxObtaining the physical resolution Pr of each pixel of the collected celadon SD-OCT image by the processing of the step 1.3) of the step 1)x;
5.2) converting the distance between the upper boundary and the lower boundary of the glaze layer on the celadon glaze layer SD-OCT image obtained in the step 4) into the actual glaze layer thickness T.
Claims (4)
1. A method for measuring the thickness of a celadon glaze layer of an SD-OCT image is characterized by comprising the following steps
1) Establishing a database of refractive indexes of different types of celadon glaze layers;
the step 1) is specifically as follows:
1.1) collecting celadon samples of different types and with known glaze layer thickness;
1.2) acquiring an SD-OCT image of the celadon glaze layer by using an OQ LabScope system, taking a 'Depth Per Pixel' parameter in the OQ LabScope system as the physical resolution Pr of each Pixel in the SD-OCT image, measuring the thickness of the celadon glaze layer by using a system self-provided caliper tool, adjusting the 'Depth Per Pixel' parameter, and further adjusting the physical resolution Pr until the thickness of the celadon glaze layer measured by the system self-provided caliper tool is the same as the actual thickness of the celadon glaze layer;
1.3) recording the physical resolution Pr of each pixel of the SD-OCT image of the celadon glaze layer, and calculating the refractive index n of each celadon glaze layer according to the optical resolution Or and the physical resolution Pr as follows:
1.4) sequentially scanning all types of celadon glaze layers, calculating the refractive indexes of all types of celadon glaze layers according to the step 1.3) and establishing a celadon glaze layer refractive index database;
2) collecting an SD-OCT image of the celadon glaze layer;
3) extracting the upper boundary of the celadon glaze layer;
the step 3) is specifically as follows:
3.1) carrying out binarization processing on the filtered image;
3.2) searching each row of pixel points in the binarized image from top to bottom to present a pixel point with a first gray value of 1, recording the point as a pixel point to be fitted, and recording the number of rows Top (i) of the pixel point to be fitted in each row, wherein i represents the number of rows of the pixel point to be fitted;
3.3) recording the median of the number of rows top (i) where the to-be-fitted pixel points of all columns are located as m, and starting to check and process the value of the number of rows top (i) where the to-be-fitted pixel points of each column are located from i ═ 1: if the absolute value of the difference between the number of rows Top (1) where the pixel point to be fitted in the first row is located and the median m is larger than 30, updating and replacing the number of rows Top (1) with the median between the number of rows Top (2) where the pixel point to be fitted in the second row is located and the number of rows Top (I) where the pixel point to be fitted in the last row is located, wherein I represents the number of columns in the last row; if the absolute value of the difference between the number of rows Top (2) where the pixel points to be fitted in the second row are located and the median m is greater than 30, updating the value of the number of rows Top (2) into the number of rows Top (1); if i is greater than 2 and the absolute value of the difference between the row number Top (i) of the pixel point to be fitted in the ith column and the row number Top (i-1) of the pixel point to be fitted in the i-1 column is greater than 30, updating the value of the row number Top (i) by using a Lagrange interpolation method, and calculating as follows:
Top(i)=2×Top(i-1)-Top(i-2)
wherein i represents the row ordinal number of the pixel point to be fitted;
3.4) after updating the number of rows Top (i) where the pixel points to be fitted of all columns are located, forming an upper boundary pixel point by all the updated pixel points to be fitted, and calculating to obtain the mean value Top of the number of rows Top (i) where the pixel points to be fitted of all columns are locatedAverage;
4) Extracting the lower boundary of the celadon glaze layer;
5) and calculating the thickness of the celadon glaze layer.
2. The method for measuring the thickness of the celadon glaze layer of the SD-OCT image as claimed in claim 1, wherein: the step 2) is specifically as follows: and (3) acquiring an SD-OCT (secure digital-optical coherence tomography) image of the celadon with the glaze layer thickness to be measured by using an OQ LabScope system, and performing median filtering on the image by using a 3 multiplied by 5 size template.
3. The method for measuring the thickness of the celadon glaze layer of the SD-OCT image as claimed in claim 1, wherein: the step 4) is specifically as follows:
4.1) carrying out edge detection on the filtered image by using a canny operator, and taking the generated image as an edge detection image Fb;
4.2) defining a circular structural element se with the radius of five pixels, and performing a closing operation on the edge detection image Fb by using the structural element se to obtain an image Fc, wherein the closing operation is calculated as follows:
Fc=Fb·se
wherein,symbolRepresenting dilation operations on imagesMake, symbolRepresenting an erosion operation on the image;
4.3) taking connected pixel points as a target in the image Fc, wherein the connection means that eight neighborhoods of the Sudoku are connected around the pixel points; then counting the number of pixel points of the target, arranging all targets from small to large according to the number of the pixel points, taking the number of the pixel points positioned in the median as a point threshold, deleting the targets with the number of the pixel points smaller than the point threshold, and then dividing the targets into two types: the first-class target is a target with the number of lines occupied by the contained pixel points larger than the number of columns, and is reserved; the other type of target is a target that the number of lines occupied by the contained pixel points is less than the number of columns, and deletion is carried out;
4.4) counting the number of rows of the most-lower pixel points at the row positions in all the targets remained in the image Fc, and solving the average value Bot of the row positions of all the pixel pointsAverage;
4.5) in the SD-OCT image of the celadon glaze layer, the D pixel distance between each column of the upper glaze layer boundary and the lower glaze layer boundary is calculated as follows:
D=BotAverage-TopAverage
wherein, TopAverageRepresenting the mean value of the number of rows top (i) where the pixel points to be fitted of all the columns after the updating is finished;
4.6) translating the upper boundary pixel point obtained in the step 3) downwards by a distance of D pixels according to a row to serve as the position of a lower boundary pixel point of the celadon glaze layer, and recording the position as bot (i).
4. The method for measuring the thickness of the celadon glaze layer of the SD-OCT image as claimed in claim 1, wherein: the step 5) is specifically as follows:
5.1) obtaining the refractive index n corresponding to the celadon according to the type x of the measured celadonxThe physical resolution Pr of each pixel of the collected celadon SD-OCT image is obtained by the processing of the step 1)x;
5.2) converting the distance between the upper boundary and the lower boundary of the glaze layer on the celadon glaze layer SD-OCT image obtained in the step 4) into the actual glaze layer thickness T, and calculating as follows:
T=D×Prx
wherein D represents the distance between the upper and lower boundaries of the glaze layer on the celadon glaze layer SD-OCT image.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910491840.5A CN110146027B (en) | 2019-06-06 | 2019-06-06 | Method for measuring thickness of celadon glaze layer of SD-OCT image |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910491840.5A CN110146027B (en) | 2019-06-06 | 2019-06-06 | Method for measuring thickness of celadon glaze layer of SD-OCT image |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110146027A CN110146027A (en) | 2019-08-20 |
CN110146027B true CN110146027B (en) | 2021-02-02 |
Family
ID=67590752
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910491840.5A Active CN110146027B (en) | 2019-06-06 | 2019-06-06 | Method for measuring thickness of celadon glaze layer of SD-OCT image |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110146027B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117419650B (en) * | 2023-12-18 | 2024-02-27 | 湖南西欧新材料有限公司 | Alumina ceramic surface glaze layer thickness measuring method based on visual analysis |
CN117968545A (en) * | 2024-03-28 | 2024-05-03 | 钛玛科(北京)工业科技有限公司 | Thickness measuring method and system for foam cotton with patterns |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103047938B (en) * | 2013-01-05 | 2015-08-12 | 山西省电力公司大同供电分公司 | Electric power line ice-covering thickness detection method and pick-up unit |
US9615735B2 (en) * | 2014-01-31 | 2017-04-11 | University Of Rochester | Measurement of the lipid and aqueous layers of a tear film |
US20180315185A1 (en) * | 2015-10-28 | 2018-11-01 | Loma Linda University | System and method to analyze various retinal layers |
CN106023158B (en) * | 2016-05-10 | 2018-09-18 | 浙江科技学院 | The fresh water pipless pearl pearly layer defect identification method of SD-OCT images |
CN106332713B (en) * | 2016-08-16 | 2019-06-11 | 浙江科技学院 | A kind of loquat early stage bruise discrimination method of SD-OCT image |
-
2019
- 2019-06-06 CN CN201910491840.5A patent/CN110146027B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN110146027A (en) | 2019-08-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yiyang | The design of glass crack detection system based on image preprocessing technology | |
CN109804232B (en) | Asphalt pavement crack development degree detection method based on infrared thermography analysis | |
CN109658398B (en) | Part surface defect identification and evaluation method based on three-dimensional measurement point cloud | |
CN108764229B (en) | Water gauge image automatic identification method based on computer vision technology | |
CN107527067B (en) | Railway roadbed disease intelligent identification method based on ground penetrating radar | |
Wang et al. | Comparison analysis on present image-based crack detection methods in concrete structures | |
WO2018040302A1 (en) | Method and system for measuring width of cut tobacco piece or cut stem, and server having said system | |
CN110146027B (en) | Method for measuring thickness of celadon glaze layer of SD-OCT image | |
CN105574161B (en) | A kind of brand logo key element recognition methods, device and system | |
CN112926563B (en) | Fault diagnosis system for steel coil spray printing mark | |
CN114549446A (en) | Cylinder sleeve defect mark detection method based on deep learning | |
CN111882573B (en) | Cultivated land block extraction method and system based on high-resolution image data | |
CN111968079B (en) | Three-dimensional pavement crack extraction method based on local extremum of section and segmentation sparsity | |
CN115731257A (en) | Leaf form information extraction method based on image | |
CN106324708B (en) | Digitizing solution, the device of rainfall record drawing | |
CN105913426B (en) | A kind of shallow lake purse seine area extracting method based on ZY-3 images | |
CN107993243B (en) | Wheat tillering number automatic detection method based on RGB image | |
CN116883402B (en) | Corrugated paper production quality real-time detection method based on computer vision | |
CN111539951B (en) | Visual detection method for outline size of ceramic grinding wheel head | |
CN113076959A (en) | Concrete structure surface microcrack feature extraction method based on convolutional neural network | |
CN110298857B (en) | Automatic celadon glaze layer thickness measuring method based on SD-OCT image | |
CN106446913A (en) | Pointer position extracting method of multi-dial pointer type instrument | |
Di Angelo et al. | A review on computer-based methods for archeological pottery classification and reconstruction | |
CN113409302B (en) | OCT image-based corn kernel early mildew identification method | |
CN109916923A (en) | A kind of customization plate automatic defect detection method based on machine vision |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |