CN114354607A - Luminosity three-dimensional flaw detection method based on spiral phase contrast filtering algorithm - Google Patents

Luminosity three-dimensional flaw detection method based on spiral phase contrast filtering algorithm Download PDF

Info

Publication number
CN114354607A
CN114354607A CN202111174769.1A CN202111174769A CN114354607A CN 114354607 A CN114354607 A CN 114354607A CN 202111174769 A CN202111174769 A CN 202111174769A CN 114354607 A CN114354607 A CN 114354607A
Authority
CN
China
Prior art keywords
phase contrast
light source
detection method
method based
camera
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.)
Granted
Application number
CN202111174769.1A
Other languages
Chinese (zh)
Other versions
CN114354607B (en
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.)
Seven Ocean Metrology Shenzhen Co ltd
Original Assignee
Seven Ocean Metrology Shenzhen Co ltd
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 Seven Ocean Metrology Shenzhen Co ltd filed Critical Seven Ocean Metrology Shenzhen Co ltd
Priority to CN202111174769.1A priority Critical patent/CN114354607B/en
Publication of CN114354607A publication Critical patent/CN114354607A/en
Application granted granted Critical
Publication of CN114354607B publication Critical patent/CN114354607B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30Computing systems specially adapted for manufacturing

Landscapes

  • Length Measuring Devices By Optical Means (AREA)
  • Investigating Materials By The Use Of Optical Means Adapted For Particular Applications (AREA)

Abstract

The invention discloses a luminosity three-dimensional flaw detection method based on a spiral phase contrast filtering algorithm, which comprises a data acquisition part and a data processing part, wherein the data acquisition part is used for acquiring luminosity three-dimensional flaw data; the data acquisition part comprises: obliquely emitting light emitted by the light source to the surface of a Lambert reflecting object to form diffuse reflection; the resolution required by object defects is calculated, and the proper multiplying power and the camera size of the required objective lens are calculated; shooting and storing a reflected light intensity image by the objective lens and the camera group; the light source is respectively placed on a plurality of directions of the object to obliquely irradiate the object, and a plurality of reflected light intensity graphs corresponding to the plurality of directions are acquired. The invention comprises the following steps: 1. spiral match filtering is added into the luminosity stereo algorithm, the flaw position edge is enhanced, the surface background of the object is removed, background interference is avoided, and the image contrast is enhanced. 2. Compared with the existing threshold segmentation background and flaw method, the method has more obvious background removing effect. 3. And gamma compression is introduced to process the image, so that the image is smoother and the human visual characteristics are compensated.

Description

Luminosity three-dimensional flaw detection method based on spiral phase contrast filtering algorithm
Technical Field
The invention relates to the fields of spiral phase contrast filtering algorithm in the field of phase imaging and luminosity three-dimensional algorithm in machine vision, in particular to a luminosity three-dimensional flaw detection method based on the spiral phase contrast filtering algorithm.
Background
In the last century, in order to perform label-free imaging on an object, zernike developed a phase contrast microscopic imaging technology, Gabor proposed holography, researchers developed a spiral phase contrast filtering imaging method combining phase contrast microscopy and holography, and an image generated by the label-free technology contains refractive index information of the thickness and the surface structure of the object, enhances and highlights the edge of the surface structure of the object, removes and balances the surface background, so that the edge of the surface structure has higher contrast with the background, and simultaneously avoids background interference caused by a strong background.
In recent years, it has become increasingly important to capture the "look" of an object from an image. By "appearance" is generally meant a model that is able to predict the image of an object under all possible fields of view and lighting conditions. The appearance of an object is adequately sampled because it does vary arbitrarily in shape and reflectivity, requiring images combined from various views and light sources, which is impractical in most cases. Fortunately, objects in the real world often exhibit regularity, which can be exploited to greatly reduce the number of images required. Therefore, effective (or near effective) constraints are chosen and strong enough to be used in practical systems, which is crucial for appearance capture. Currently, there are stereovision inspection works that have shown that under known lighting conditions, camera viewing angle determination is feasible to recover a definite appearance model from the changing image. This is an important special case appearance capture because it relies only on photometric cues to avoid solving the correspondence problem. This is also important because such explicit appearance models have proven to be useful for visual tasks such as flaw detection. However, in the field of defect detection, the image processed by the photometric stereo method needs background elimination and balanced smoothing, only the defect position is highlighted, and the background and the defect position form the sharpest contrast. The general approach to achieve this is to set a threshold to separate the background from the defect, but this approach sometimes causes misjudgment due to the close pixel values of the background and the defect, and the background and the defect cannot be completely separated. Therefore, a photometric stereo flaw detection method based on a spiral phase contrast filtering algorithm is urgently needed to solve the problem existing in the prior art.
In order to solve the technical problems, a new technical scheme is especially provided.
Disclosure of Invention
The invention aims to provide a luminosity three-dimensional flaw detection method based on a spiral phase contrast filtering algorithm, so as to solve the problems in the background technology.
In order to achieve the purpose, the invention provides the following technical scheme: a luminosity three-dimensional flaw detection method based on a spiral phase contrast filtering algorithm comprises a data acquisition part and a data processing part;
the data acquisition part comprises the following steps:
step a, obliquely emitting light emitted by a light source to the surface of a Lambertian reflecting object;
c, forming diffuse reflection on the surface of the object by the incident light;
step c, counting the required resolution of object flaws, and calculating the required lens magnification and the number of camera pixels according to the required resolution;
d, placing the lens and the camera combination right above the object to form orthogonal projection and collect reflected light carrying the surface structure information of the object;
step e, connecting the camera with a computer, shooting through camera software on the computer and storing the intensity map of the reflected light;
and f, respectively placing the light sources on a plurality of directions of the object to obliquely irradiate the object, repeating the steps a to e for a plurality of times, and acquiring a plurality of reflected light intensity maps corresponding to the plurality of directions.
Preferably, the data processing section includes the steps of:
step 1, calibrating the direction of a light source: first orthogonal projection scenePlacing a smooth ball in the system, shooting a reflected light intensity map under the conditions of 4 different light source directions, placing the 4 intensity maps into MATLAB to calculate and fit the circle boundary of the smooth ball, and positioning the coordinate (xc) of the circle centeri,yci) I is an index of 4 images;
and 2, respectively positioning the highlight positions of the surfaces of the spheres on the 4 graphs, wherein the coordinates of the highlight positions are expressed as (hxc)i,hyci) And the direction of the light source is reflected by the highlight position of the surface;
and 3, calculating a normal vector of the surface of the smooth sphere according to the following formula:
n:nx=hxci-xci,ny=hyci-yci
Figure BDA0003294935210000031
step 4, estimating the illumination intensity: an illumination parameter cost function is defined by the following formula:
Figure BDA0003294935210000032
where N is the number of pixels per intensity map, j is the pixel index on the intensity map, I is the intensity value (pixel value) on each intensity map, ρ is the reflectivity of the object,
Figure BDA0003294935210000033
is the unit vector of the ith light source direction, λiIs the ith intensity of illumination, λi=||ei||;
Step 5, performing object surface structure reflectivity rho obtained by photometric stereojSum normal vector nj:ρj=||bj||,
Figure BDA0003294935210000034
Step 6, judging the defect type and the characteristics of the position of the defect type in the image according to the two-dimensional image, and measuring the third channel data n (3) of the defect type according to the obtained normal amount;
step 7, writing a spiral phase plate function
Figure BDA0003294935210000035
Step 8, writing a 4f system, placing n (3) as an object at the position of an object plane of the 4f system, performing spiral phase contrast imaging, and solving a complex function U at the image plane:
Figure BDA0003294935210000036
f represents a Fourier transform;
step 9, taking the amplitude A of the image plane complex function U:
Figure BDA0003294935210000037
step 10, carrying out gamma compression on the obtained amplitude A to obtain a final object flaw image Iout:Iout=BAγ. B is a constant, gamma is a correction parameter, and a fraction less than 1 is taken and is called the encoding gamma value.
Preferably, the formula of the direction of the light source is: l is 2nzn-R, wherein R represents a camera direction of (0,0,1) and coordinates of (0,0, -1) if viewed from the camera perspective;
preferably, the minimization formula of the cost function is:
Figure BDA0003294935210000041
Figure BDA0003294935210000042
bj=ρjnjrepresents an unknown scene attribute at j pixels;
preferably, the first and second electrodes are formed of a metal,
Figure BDA0003294935210000043
is a non-linear least squares problem that can be solved by the Levenberg-Marquardt algorithm.
Preferably, for all eiHas Ij=LTbjThen b can be obtained according to the least square methodj=(LLT)-1LIj
Preferably, the plurality of orientations of the object in step f refers to selecting a plurality of angles in the 360 ° direction of the object, for example, four angles of 0 °, 90 °, 180 °, and 270 °.
Preferably, said plurality of steps f refers to 3 or more, preferably 4.
Preferably, the step f is to place the light source on the front, back, left, and right of the object respectively to irradiate the object obliquely, repeat the above steps a to e four times, and acquire 4 reflected light intensity maps corresponding to the 4 orientations.
Compared with the prior art, the invention has the beneficial effects that:
1. spiral match filtering is added into the luminosity stereo algorithm, the flaw position edge is enhanced, the surface background of the object is removed, background interference is avoided, and the image contrast is enhanced.
2. Compared with the existing threshold segmentation background and flaw method, the method has more obvious background removing effect.
3. And gamma compression is introduced to process the image, so that the image is smoother and the human visual characteristics are compensated.
Drawings
FIG. 1 is a schematic diagram of a photometric stereo imaging optical set providing illumination with a light source tilted at a 0 position above an object from a camera perspective.
FIG. 2 is a schematic diagram of a photometric stereo imaging optical set providing illumination with a light source tilted at a 90 position above an object from a camera perspective.
FIG. 3 is a schematic diagram of a photometric stereo imaging optical set providing illumination with a light source tilted 180 above an object from a camera perspective.
FIG. 4 is a schematic diagram of a photometric stereo imaging optical set providing illumination of a light source tilted at a 270 position above an object from a camera perspective.
Fig. 5 is a diagram illustrating a lambertian reflection characteristic.
FIG. 6 is a schematic view of the light source direction calibration.
Fig. 7 is a schematic diagram of a spiral phase plate.
FIG. 8 is a camera shot intensity chart of the light source irradiating the steel material of the automobile door at the 0 degree position.
FIG. 9 is a camera shot intensity chart of the light source irradiating the steel material of the automobile door at the 90 degree position.
FIG. 10 is a camera shot intensity chart of the light source irradiating the steel material of the automobile door at the 180 degree position.
FIG. 11 is a photograph of the camera intensity of the light source irradiating the steel material of the door of the automobile at 270 deg.
Fig. 12 is an image of the 4 shot intensity maps of fig. 8 to 11 after being processed by the photometric stereo algorithm.
Fig. 13 is the image of fig. 12 after being processed by the spiral phase contrast filtering algorithm.
FIG. 14 is the gamma compressed image of FIG. 13.
Fig. 15 is a photograph intensity chart of the light source illuminating the hoop camera at the 0 ° position.
Fig. 16 is a photograph intensity diagram of a light source illuminating a hoop camera at a 90 ° position.
Fig. 17 is a diagram of the shooting intensity of the light source illuminating the hoop camera at the 180 ° position.
Fig. 18 is a photograph of the intensity of the light source illuminating the hoop camera at the 270 ° position.
Fig. 19 is an image of the 4 shot intensity maps of fig. 15 to 18 after being processed by a photometric stereo algorithm.
Fig. 20 is the image of fig. 19 after being processed by the spiral phase contrast filtering algorithm.
FIG. 21 is the gamma compressed image of FIG. 20.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Referring to the attached drawings of the specification, the invention provides a technical scheme that: a luminosity three-dimensional flaw detection method based on a spiral phase contrast filtering algorithm comprises the following steps:
firstly, a data acquisition part.
1. The light emitted by the light source is obliquely incident on the surface of the lambertian reflecting object.
2. The incident light forms diffuse reflection on the surface of the object.
3. And counting the required resolution of the object defects, and calculating the required lens magnification and the number of camera pixels according to the required resolution.
4. The lens and the camera are combined and placed right above the object to form orthogonal projection to collect reflected light carrying the surface structure information of the object.
5. The camera is connected with a computer, and the intensity map of the reflected light is shot and stored through camera software on the computer.
6. The light source is respectively placed in 4 directions of the front, the back, the left and the right of the object to obliquely irradiate the object, the steps 1 to 5 are repeated four times, and 4 reflected light intensity graphs corresponding to the 4 directions are acquired (if necessary, a plurality of reflected light intensity graphs can be acquired, but at least 3 reflected light intensity graphs are acquired and are uniformly distributed on the 360-degree direction of the object).
And secondly, a data processing part.
8. Calibrating the light source direction: firstly, placing a smooth ball in an orthogonal projection scene, shooting a reflected light intensity map under the conditions of 4 different light source directions, putting 4 intensity maps into MATLAB to calculate and fit a circle boundary of the smooth ball, and positioning a circle center coordinate (xc)i,yci) And i is an index of 4 images.
9. The surface highlights of the sphere on the 4 graphs are located separately, and the coordinates of the highlights are represented as (hxc)i,hyci) And the direction of the light source is reflected by the highlight position of the surface.
10. Calculating the normal vector n of the smooth sphere surfacex=hxci-xci,ny=hyci-yci
Figure BDA0003294935210000061
11. Light source direction L2 nzn-R, R indicates that the camera direction is (0,0,1), and the coordinates are (0,0, -1) when viewed from the camera perspective.
12. Estimating the illumination intensity: an illumination parameter cost function is defined,
Figure BDA0003294935210000062
n is the number of pixels per intensity map, j is the pixel index on the intensity map, I is the intensity value (pixel value) on each intensity map, p is the reflectivity of the object,
Figure BDA0003294935210000071
is the unit vector of the ith light source direction, λiIs the ith intensity of illumination, λi=||ei||。
Minimization of cost function
Figure BDA0003294935210000072
Figure BDA0003294935210000073
bj=ρjnjRepresenting the unknown scene properties at j pixels.
14.
Figure BDA0003294935210000074
This is a non-linear least squares problem that can be solved by the Levenberg-Marquardt algorithm.
15. For all eiHas Ij=LTbjThen b can be obtained according to the least square methodj=(LLT)-1LIj
16. Object surface structure reflectivity rho for executing photometric stereo solutionjSum normal vector nj:ρj=||bj||,
Figure BDA0003294935210000075
17. Since only two-dimensional images are needed to determine the defect type and its position in the image, the obtained normal amount is taken as the third channel data n (3).
18. Writing a spiral phase plate function
Figure BDA0003294935210000076
19. Writing a 4f system, placing n (3) as an object at the position of an object plane of the 4f system, performing spiral phase contrast imaging, and solving a complex function U at an image plane:
Figure BDA0003294935210000077
f denotes fourier transform.
20. Taking the amplitude A of the image plane complex function U:
Figure BDA0003294935210000078
21. carrying out gamma compression on the obtained amplitude A to obtain a final object flaw image Iout:Iout=BAγ. B is a constant, gamma is a correction parameter, and a fraction less than 1 is taken and is called the encoding gamma value.
The image forming apparatus of the present invention includes: 4 light sources with consistent light intensity, object to be measured, objective lens, industrial camera set and computer
Wherein the objective lens and the industrial camera are connected in combination, and the sensor of the industrial camera is located at the focal plane position of the objective lens. The computer is connected with the industrial camera through a data line and is used for storing the 4 reflection intensity images shot by the camera into the computer and carrying out algorithm processing on the collected images on the computer.
As shown in fig. 1, the method comprises the following steps:
s1, from the camera view angle, 4 light sources with consistent light intensity are respectively placed in four directions of 0 degree, 90 degrees, 180 degrees and 270 degrees of an object to be detected for illumination. The 4 light sources are sequentially turned on to irradiate the object to be detected, and only one light source is turned on at a time.
And S2.4, the light emitted by the light sources is sequentially incident on the surface of the object to be measured to form diffuse reflection.
And S3, the objective lens and the camera group are positioned right above the object to be detected and used for sequentially collecting reflected light transmitted from the object, and 4 signals are sequentially formed on the camera.
And S4, sequentially transmitting 4 signals on the camera to a computer through a data line to form 4 images.
In use, the invention: 1. spiral match filtering is added into the luminosity stereo algorithm, the flaw position edge is enhanced, the surface background of the object is removed, background interference is avoided, and the image contrast is enhanced. 2. Compared with the existing threshold segmentation background and flaw method, the method has more obvious background removing effect. 3. And gamma compression is introduced to process the image, so that the image is smoother and the human visual characteristics are compensated.
Although embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that changes, modifications, substitutions and alterations can be made in these embodiments without departing from the principles and spirit of the invention, the scope of which is defined in the appended claims and their equivalents.

Claims (8)

1. A luminosity three-dimensional flaw detection method based on a spiral phase contrast filter algorithm is characterized by comprising a data acquisition part and a data processing part;
the data acquisition part comprises the following steps:
step a, obliquely emitting light emitted by a light source to the surface of a Lambertian reflecting object;
b, forming diffuse reflection on the surface of the object by the incident light;
step c, counting the required resolution of object flaws, and calculating the required lens magnification and the number of camera pixels according to the required resolution;
d, placing the lens and the camera combination right above the object to form orthogonal projection and collect reflected light carrying the surface structure information of the object;
step e, connecting the camera with a computer, shooting through camera software on the computer and storing the intensity map of the reflected light;
and f, respectively placing the light sources on a plurality of directions of the object to obliquely irradiate the object, repeating the steps a to e for a plurality of times, and acquiring a plurality of reflected light intensity maps corresponding to the plurality of directions.
The data processing section includes the steps of:
step 1, calibrating the direction of a light source: firstly, placing a smooth ball in an orthogonal projection scene, shooting a reflected light intensity map under the conditions of 4 different light source directions, putting 4 intensity maps into MATLAB to calculate and fit a circle boundary of the smooth ball, and positioning a circle center coordinate (xc)t,yct) I is an index of 4 images;
and 2, respectively positioning the highlight positions of the surfaces of the spheres on the 4 graphs, wherein the coordinates of the highlight positions are expressed as (hxc)i,hyci) And the direction of the light source is reflected by the highlight position of the surface;
and 3, calculating a normal vector of the surface of the smooth sphere according to the following formula:
Figure RE-FDA0003531682000000011
step 4, estimating the illumination intensity: an illumination parameter cost function is defined by the following formula:
Figure RE-FDA0003531682000000012
where N is the number of pixels per intensity map, j is the pixel index on the intensity map, I is the intensity value (pixel value) on each intensity map, ρ is the reflectivity of the object,
Figure RE-FDA0003531682000000021
Figure RE-FDA0003531682000000022
is the unit vector of the ith light source direction, λiIs the ith intensity of illumination, λi=||ei||,
Step 5, performing object surface structure reflectivity rho obtained by photometric stereojSum normal vector nj:ρj=||bj||,
Figure RE-FDA0003531682000000023
Step 6, judging the defect type and the characteristics of the position of the defect type in the image according to the two-dimensional image, and measuring the third channel data n (3) of the defect type according to the obtained normal amount;
step 7, writing a spiral phase plate function
Figure RE-FDA0003531682000000024
Step 8, writing a 4f system, placing n (3) as an object at the position of an object plane of the 4f system, performing spiral phase contrast imaging, and solving a complex function U at the image plane:
Figure RE-FDA0003531682000000025
f represents a Fourier transform;
step 9, taking the amplitude A of the image plane complex function U:
Figure RE-FDA0003531682000000026
step 10, carrying out gamma compression on the obtained amplitude A to obtain a final object flaw image Iout:IoutBA γ. B is a constant, gamma is a correction parameter, and a fraction less than 1 is taken and is called the encoding gamma value.
2. The photometric stereo flaw detection method based on the spiral phase contrast filtering algorithm according to claim 1, characterized in that: the direction formula of the light source is as follows: l is 2nzn-R, where R denotes the camera direction as (0,0,1) and coordinates as (0,0, -1) if viewed from the camera perspective.
3. The photometric stereo flaw detection method based on the spiral phase contrast filtering algorithm according to claim 1, characterized in that: the minimization formula of the cost function is:
Figure RE-FDA0003531682000000027
Figure RE-FDA0003531682000000028
bj=ρjnjrepresenting the unknown scene properties at j pixels.
4. The photometric stereo flaw detection method based on the spiral phase contrast filtering algorithm according to claim 3, characterized in that:
Figure RE-FDA0003531682000000031
is a non-linear least squares problem that can be solved by the Levenberg-Marquardt algorithm.
5. The photometric stereo flaw detection method based on the spiral phase contrast filtering algorithm according to claim 4, characterized in that: for all eiHas Ij=LTbjThen b can be obtained according to the least square methodj=(LLT)-1LIj
6. The photometric stereo flaw detection method based on the spiral phase contrast filtering algorithm according to claim 1, characterized in that: the multiple orientations of the object in step f refer to selecting multiple angles in the 360-degree direction of the object, and take four angles of 0 degrees, 90 degrees, 180 degrees and 270 degrees as examples.
7. The photometric stereo defect detection method based on the helical phase contrast filtering algorithm according to claim 6, characterized in that: the plurality of steps f refers to 3 or more, preferably 4.
8. The photometric stereo flaw detection method based on the helical phase contrast filtering algorithm according to claim 7, characterized in that: and step f specifically, placing the light source on the front, back, left and right of the object respectively to obliquely irradiate the object in 4 directions, repeating the steps a to e four times, and acquiring 4 corresponding reflected light intensity maps in 4 directions.
CN202111174769.1A 2021-10-09 2021-10-09 Luminosity three-dimensional flaw detection method based on spiral phase contrast filtering algorithm Active CN114354607B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111174769.1A CN114354607B (en) 2021-10-09 2021-10-09 Luminosity three-dimensional flaw detection method based on spiral phase contrast filtering algorithm

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111174769.1A CN114354607B (en) 2021-10-09 2021-10-09 Luminosity three-dimensional flaw detection method based on spiral phase contrast filtering algorithm

Publications (2)

Publication Number Publication Date
CN114354607A true CN114354607A (en) 2022-04-15
CN114354607B CN114354607B (en) 2023-10-31

Family

ID=81096345

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111174769.1A Active CN114354607B (en) 2021-10-09 2021-10-09 Luminosity three-dimensional flaw detection method based on spiral phase contrast filtering algorithm

Country Status (1)

Country Link
CN (1) CN114354607B (en)

Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010025652A (en) * 2008-07-17 2010-02-04 Nippon Steel Corp Surface flaw inspection device
CN101872491A (en) * 2010-05-21 2010-10-27 清华大学 Free view angle relighting method and system based on photometric stereo
US20110292179A1 (en) * 2008-10-10 2011-12-01 Carlos Hernandez Imaging system and method
CN103363924A (en) * 2013-07-15 2013-10-23 中国科学院空间科学与应用研究中心 Compressing three-dimension calculation ghost imaging system and method
US20140376003A1 (en) * 2012-01-05 2014-12-25 Helmee Imaging Oy Arrangement for optical measurements and related method
US20150346115A1 (en) * 2014-05-30 2015-12-03 Eric J. Seibel 3d optical metrology of internal surfaces
JP2015232487A (en) * 2014-06-09 2015-12-24 株式会社キーエンス Image inspection device, image inspection method, image inspection program and computer readable recording medium, and apparatus having image inspection program recorded therein
JP2016186421A (en) * 2015-03-27 2016-10-27 株式会社Ihi Image processing method
CN109596637A (en) * 2017-09-30 2019-04-09 湖南海擎智能科技有限责任公司 Object bumps flaw visible detection method and detection device
CN110455815A (en) * 2019-09-05 2019-11-15 西安多维机器视觉检测技术有限公司 A kind of method and system of electronic component open defect detection
CN111896550A (en) * 2020-03-31 2020-11-06 广西师范大学 Surface defect detection device and method
CN112734756A (en) * 2021-03-29 2021-04-30 惠州高视科技有限公司 Detection method and system based on photometric stereo vision
CN112734824A (en) * 2021-01-26 2021-04-30 中国科学院空天信息创新研究院 Three-dimensional reconstruction method based on generalized luminosity stereo model
CN112986284A (en) * 2021-02-20 2021-06-18 苏州天准科技股份有限公司 Combined type light source structure and optical detection device
CN113252567A (en) * 2021-06-08 2021-08-13 菲特(天津)检测技术有限公司 Method, system, medium and terminal for rapidly detecting multiple defects on surface of aluminum valve plate

Patent Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010025652A (en) * 2008-07-17 2010-02-04 Nippon Steel Corp Surface flaw inspection device
US20110292179A1 (en) * 2008-10-10 2011-12-01 Carlos Hernandez Imaging system and method
CN101872491A (en) * 2010-05-21 2010-10-27 清华大学 Free view angle relighting method and system based on photometric stereo
US20140376003A1 (en) * 2012-01-05 2014-12-25 Helmee Imaging Oy Arrangement for optical measurements and related method
CN103363924A (en) * 2013-07-15 2013-10-23 中国科学院空间科学与应用研究中心 Compressing three-dimension calculation ghost imaging system and method
US20150346115A1 (en) * 2014-05-30 2015-12-03 Eric J. Seibel 3d optical metrology of internal surfaces
JP2015232487A (en) * 2014-06-09 2015-12-24 株式会社キーエンス Image inspection device, image inspection method, image inspection program and computer readable recording medium, and apparatus having image inspection program recorded therein
JP2016186421A (en) * 2015-03-27 2016-10-27 株式会社Ihi Image processing method
CN109596637A (en) * 2017-09-30 2019-04-09 湖南海擎智能科技有限责任公司 Object bumps flaw visible detection method and detection device
CN110455815A (en) * 2019-09-05 2019-11-15 西安多维机器视觉检测技术有限公司 A kind of method and system of electronic component open defect detection
CN111896550A (en) * 2020-03-31 2020-11-06 广西师范大学 Surface defect detection device and method
CN112734824A (en) * 2021-01-26 2021-04-30 中国科学院空天信息创新研究院 Three-dimensional reconstruction method based on generalized luminosity stereo model
CN112986284A (en) * 2021-02-20 2021-06-18 苏州天准科技股份有限公司 Combined type light source structure and optical detection device
CN112734756A (en) * 2021-03-29 2021-04-30 惠州高视科技有限公司 Detection method and system based on photometric stereo vision
CN113252567A (en) * 2021-06-08 2021-08-13 菲特(天津)检测技术有限公司 Method, system, medium and terminal for rapidly detecting multiple defects on surface of aluminum valve plate

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
BANG-HWAN KIM等: "Shape from Shading and Photometric Stereo Using Surface Approximation by Legendre Polynomials", COMPUTER VISION AND IMAGE UNDERSTANDING, vol. 66, no. 03 *
GEORG RADOW等: "Optimisation of Classic Photometric Stereo by Non-convex Variational Minimisation", JOURNAL OF MATHEMATICAL IMAGING AND VISION, vol. 61 *
LONG MA等: "A fast LED calibration method under near field lighting based on photometric stereo", OPTICS AND LASERS IN ENGINEERING, vol. 147 *
程耀瑜等: "基于光度立体技术的地板疵点检测方法", 科学技术与工程, vol. 17, no. 30 *

Also Published As

Publication number Publication date
CN114354607B (en) 2023-10-31

Similar Documents

Publication Publication Date Title
CN107607040B (en) Three-dimensional scanning measurement device and method suitable for strong reflection surface
JP4324238B2 (en) Image processing system, method and apparatus
JP6319329B2 (en) Surface attribute estimation using plenoptic camera
US9562857B2 (en) Specular object scanner for measuring reflectance properties of objects
EP3382645B1 (en) Method for generation of a 3d model based on structure from motion and photometric stereo of 2d sparse images
CN108846837B (en) Object surface defect detection method and device
CN113205593A (en) High-light-reflection surface structure light field three-dimensional reconstruction method based on point cloud self-adaptive restoration
JP2016186421A (en) Image processing method
CN113252568A (en) Lens surface defect detection method, system, product and terminal based on machine vision
JP2024507089A (en) Image correspondence analysis device and its analysis method
Ciortan et al. A practical reflectance transformation imaging pipeline for surface characterization in cultural heritage
KR20230096057A (en) Defect Layering Detection Method and System Based on Light Field Camera and Detection Production Line
CN113252567A (en) Method, system, medium and terminal for rapidly detecting multiple defects on surface of aluminum valve plate
CN111208138A (en) Intelligent wood recognition device
KR20060065800A (en) Apparatus and method for 3d reconstruction of objects with depth discontinuities using helmholtz stereopsis
CN114354607B (en) Luminosity three-dimensional flaw detection method based on spiral phase contrast filtering algorithm
CN116380918A (en) Defect detection method, device and equipment
WO2019238583A1 (en) Deflectometric techniques
CN115272258A (en) Metal cylindrical surface defect detection method, system and medium based on machine vision
CN114486939A (en) Lens scratch detection system and method
CN112782176A (en) Product appearance detection method and device
TWI826988B (en) System and method for three-dimensional image evaluation
CN113592995B (en) Multi-reflection light separation method based on parallel single-pixel imaging
CN116482109A (en) Surface defect detection method and device, storage medium and electronic equipment
CN114674244B (en) Coaxial normal incidence speckle deflection measurement method and device

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