CN112729168A - F-P etalon interference imaging quality evaluation method for micro angle measurement - Google Patents

F-P etalon interference imaging quality evaluation method for micro angle measurement Download PDF

Info

Publication number
CN112729168A
CN112729168A CN202011533267.9A CN202011533267A CN112729168A CN 112729168 A CN112729168 A CN 112729168A CN 202011533267 A CN202011533267 A CN 202011533267A CN 112729168 A CN112729168 A CN 112729168A
Authority
CN
China
Prior art keywords
image
axis
fwhmy
fwhmx
ring
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
CN202011533267.9A
Other languages
Chinese (zh)
Other versions
CN112729168B (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.)
China Jiliang University
Original Assignee
China Jiliang University
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 China Jiliang University filed Critical China Jiliang University
Priority to CN202011533267.9A priority Critical patent/CN112729168B/en
Publication of CN112729168A publication Critical patent/CN112729168A/en
Application granted granted Critical
Publication of CN112729168B publication Critical patent/CN112729168B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/26Measuring arrangements characterised by the use of optical techniques for measuring angles or tapers; for testing the alignment of axes
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00Measuring instruments characterised by the use of optical techniques
    • G01B9/02Interferometers
    • G01B9/02055Reduction or prevention of errors; Testing; Calibration
    • G01B9/02062Active error reduction, i.e. varying with time

Abstract

The invention discloses an F-P etalon interference imaging quality evaluation method for measuring a tiny angle, which greatly improves the imaging quality by measuring and evaluating the light intensity distribution uniformity evaluation, the image definition evaluation, the image ovalization evaluation and the optimal selection evaluation between imaging rings of images of an interference imaging image of the F-P etalon.

Description

F-P etalon interference imaging quality evaluation method for micro angle measurement
Technical Field
The invention relates to the field of micro-angle measurement methods, in particular to an F-P etalon interference imaging quality evaluation method for micro-angle measurement.
Background
The micro angle measurement technology is widely applied to the occasions of precision mechanical component installation, precision and ultra-precision machining, aiming and positioning and the like, and has extremely important significance and effect in the fields of machinery, aerospace, military and the like. The F-P etalon interference method is one of the technical means for realizing high-precision measurement of the tiny angle at present, and the method utilizes the interference imaging principle of the F-P etalon, uses an objective lens to project the deflection angle of incident light onto an image plane, and uses an area array device to receive the circle center displacement of a concentric interference ring so as to realize the measurement of the deflection angle.
For the F-P etalon interferometry to measure the tiny angle, the previous measurement research mostly focuses on the improvement of the measurement principle, the improvement of the processing algorithm of the interference image and the like, and the research on the imaging quality aspects such as the imaging light intensity distribution uniformity, the imaging definition degree, the ovalization of the imaging ring, the optimal selection between the imaging rings and the like of the interference image is lacked, so that the tracing to the tiny angle measured by the F-P etalon interferometry is not facilitated, and the accuracy of the tiny angle measurement is influenced.
In order to solve the existing problems, a novel F-P etalon interference imaging quality evaluation method for measuring a tiny angle is needed.
Disclosure of Invention
Aiming at the defects of the prior art and a processing method, the invention provides an F-P etalon interference imaging quality evaluation method for micro-angle measurement, solves the imaging quality problems of the uniformity of imaging light intensity distribution, the imaging definition degree, the imaging ring ovalization degree, the optimal selection between imaging rings and the like of an F-P etalon interference image in micro-angle measurement, makes up for the defect that an interference image quality evaluation method is lacked in the process of measuring micro-angles by an F-P etalon interference method, and has the advantages of low cost, high efficiency and high cost performance.
In order to solve the problems, the invention is realized by the following technical scheme:
an F-P etalon interference imaging quality evaluation method for micro-angle measurement comprises the following steps:
the method comprises the following steps: and acquiring a concentric interference circular ring image. Receiving a concentric interference circular ring image before a deflection angle is generated by an area array device;
step two: and judging the light intensity uniformity of the concentric interference circular ring image. Taking the acquired concentric interference circular ring image as central axes along the horizontal direction and the vertical direction, measuring the maximum light intensity of all circular rings on the two central axes, establishing a rectangular coordinate system, taking the horizontal central axis as an x axis, the vertical central axis as a y axis, and the intersection point of the two central axes as an origin, and respectively calculating the positive direction, the negative direction and the light intensity of the x axis,The average value of the maximum light intensity value of each circular ring in the positive direction of the y axis and the negative direction of the y axis is respectively
Figure BDA0002852373500000021
By comparison
Figure BDA0002852373500000022
Figure BDA0002852373500000023
The proximity between the four values is to a threshold L2Comparing, and judging whether the uniformity of the image light intensity is qualified or not, as follows;
Figure BDA0002852373500000024
Figure BDA0002852373500000025
computing
Figure BDA0002852373500000026
Arithmetic mean of
Figure BDA0002852373500000027
Respectively calculate
Figure BDA0002852373500000028
Figure BDA0002852373500000029
And
Figure BDA00028523735000000210
so that the absolute value of the largest difference is equal to
Figure BDA00028523735000000211
Is divided to obtain L1I.e. the light intensity distribution uniformity parameter of the image, when L1>L2Then, the light intensity distribution of the image is judgedAnd if not, the uniformity is qualified.
Step three: and judging the definition of the concentric interference circular ring image. The full width at half maximum FWHM of the interference fringe, i.e. the phase difference corresponding to when the light intensity falls to half of the maximum. And moving the area array device near the ideal focal plane, receiving the concentric interference ring images of the area array device at different positions, and calculating the FWHM values of the rings with the assigned serial numbers in all the concentric interference ring images, wherein the FWHM values are respectively FWHMx +, FWHMx-, FWHMy + and FWHMy-. Comparing the FWHM arithmetic mean values of all the images, wherein the image with the minimum FWHM arithmetic mean value is the closest to an ideal focal plane, so as to obtain the clearest image, and the method comprises the following steps;
FWHMn=(FWHMx+n)+(FWHMx-n)+(FWHMy+n)+(FWHMy-n)(n=1,2,3…)
FWHM min=min(FWHM1,FWHM2,FWHM3,FWHM4…FWHMn)
the n serial numbers are used for marking different focusing positions, and the FWHMn is an arithmetic mean value of FWHM values of the serial number rings passing through the vicinity of two central axes in the concentric interference ring image at the same focusing position. FWHMmin is the minimum FWHMn value of concentric interference circular ring images at different focusing positions, namely the definition parameter of the images, and at the moment, the corresponding n is the serial number of the optimal focusing position, so that the clearest images can be obtained.
Step four: and judging the ovalization degree of the concentric interference circular ring image. And calculating the FWHM values of the designated sequence number rings passing through the vicinity of the two central axes in the concentric interference ring image, namely FWHMx +, FWHMx-, FWHMy + and FWHMy-. By comparing the closeness between four values FWHMx +, FWHMx-, FWHMy +, FWHMy-with a threshold T2And comparing to judge whether the ovalization of the ring is qualified or not, and the method comprises the following steps:
Figure BDA0002852373500000031
Figure BDA0002852373500000032
calculating the arithmetic mean of FWHMx +, FWHMx-, FWHMy +, FWHMy-)
Figure BDA0002852373500000033
FWHMx +, FWHMx-, FWHMy +, FWHMy-and
Figure BDA0002852373500000041
so that the absolute value of the largest difference is equal to
Figure BDA0002852373500000042
Is divided to obtain T1I.e. the ovalization parameter of the image, when T1>T2And judging that the ovalization degree of the image is unqualified, otherwise, judging that the ovalization degree is qualified.
Step five: and judging the optimal selection among the imaging rings of the concentric interference ring image. Respectively establishing an x-axis rectangular coordinate system and a y-axis rectangular coordinate system in the horizontal direction and the vertical direction of the diameter of the approximate concentric interference circular ring, rotating the coordinate systems by 45 degrees anticlockwise to obtain the x '-axis rectangular coordinate system and the y' -axis rectangular coordinate system, and adopting a method for interpolating, subdividing and smoothing the array virtual pixels, which is disclosed in the patent application specification with the application date of 201710374595.0. And respectively establishing N parallel lines on two sides of the approximate diameter of the circular ring of the coordinate axis of the x 'axis and the y' axis, wherein N is a positive integer. For the same ring, after the parallel line intersects with each ring, there are 8N +4 photoelectric signal small line segments, and the peak position coordinates and the standard deviation S of the photoelectric signals on each small line segment of all the rings in the concentric interference ring image are obtained by adopting the method for obtaining the peak coordinate values of the photoelectric signals in the patent application specification with the application number of 201510217472.7. Calculating the arithmetic mean value of the standard deviation of the peak position coordinates of 8N +4 small segments of the photoelectric signals in each circle
Figure BDA0002852373500000043
Comparing all rings
Figure BDA0002852373500000044
Value, select minimum
Figure BDA0002852373500000045
The circle where the value is located is used as the optimal circle for subsequently calculating the radius of the circle center, and the method comprises the following steps:
Figure BDA0002852373500000046
(m-1, 2,3 … 8N +4, N is a natural number, k-1, 2,3 …)
Figure BDA0002852373500000047
And m is the peak position coordinate and the standard deviation number on one circular ring. k is the serial number of the circular ring,
Figure BDA0002852373500000048
and calculating the arithmetic mean value of all the peak position coordinate standard deviations on the circular ring where the designated circular ring serial number is located in the clearest concentric interference circular ring image.
Figure BDA0002852373500000049
For each ring in the image of the most clearly concentric interfering rings
Figure BDA00028523735000000410
And the minimum value is the optimal selected parameter between the imaging rings of the image, and the corresponding k is the optimal ring serial number at the moment and is used for calculating the diameter of the subsequent circle center.
The invention has the following beneficial effects:
(1) the F-P etalon interference imaging quality evaluation method for measuring the tiny angle overcomes the defect that an interference image quality evaluation method is lacked in the process of measuring the tiny angle by the F-P etalon interference method.
(2) The evaluation process is standardized, the influence of human factors is reduced, the accuracy of image quality judgment is high, and a foundation is laid for the improvement of the accuracy of the measurement of the tiny angle by the F-P etalon interference method in the future.
(3) The evaluation process is easy to program, the image quality can be conveniently and rapidly judged, and the efficiency is high.
(4) The method has strong universality, can carry out measurement under different conditions by replacing a light source, an F-P etalon, an imaging objective lens and the like, and provides a certain foundation for the deep research of the technology for measuring the tiny angle by the F-P etalon interferometry.
Drawings
FIG. 1 is a flow chart of the F-P etalon interference imaging quality evaluation method in the invention.
FIG. 2 is a diagram of an apparatus used in an embodiment of the present invention.
FIG. 3 is a detailed flowchart of a specific process for recording evaluation parameters according to the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings.
As shown in fig. 2, the apparatus used in the embodiment of the present invention includes a point light source 1, an interference filter 2, a converging lens 3, an F-P etalon 4, a mirror 5, an imaging objective lens 6, and an area array device 8, the point light source 1 being a mercury lamp; the interference filter 2 is a monochromatic filter with the wavelength of 546 nm; the F-P etalon 5 is a 2mm spaced air gap F-P etalon; the flatness of the reflector 5 is better than lambda/20; the imaging objective 6 is an industrial fixed-focus lens with a focal length of 80 mm; the area array size of the area array device 9 is 22.16mm multiplied by 15.22mm, the area array pixel number is 7904 multiplied by 5432, and the average pixel pitch of the area array is w ≈ 2.8 μm. The reflector 5 deflects, the area array device 9 obtains two concentric interference ring images with changed circle center positions, the circle center displacement δ and the focal length f of the imaging objective lens 6, and the deflection angle α of the reflector 5 is δ/f. Experiments show that within the rotation angle measurement range 600 ″, the angle is not determined to be greater than 0.132 ″. The components can be stably measured in real time and have high precision, and the device is only used for explaining the specific technical scheme of the invention.
Specifically, the F-P etalon interference imaging quality evaluation method for measuring a minute angle, as shown in fig. 1 and 3, includes the following steps:
the method comprises the following steps: acquiring a concentric interference circular ring image through an area array device 9;
step two: for the obtained concentric stemsRegarding a ring image, taking the horizontal direction and the vertical direction as central axes, measuring the maximum light intensity of all rings on the two central axes, establishing a rectangular coordinate system, taking the horizontal central axis as an x axis, the vertical central axis as a y axis, taking the intersection point of the two central axes as an origin, respectively calculating the average value of the maximum light intensity value of each ring in the positive direction of the x axis, the negative direction of the x axis, the positive direction of the y axis and the negative direction of the y axis, respectively
Figure BDA0002852373500000061
By comparison
Figure BDA0002852373500000062
The proximity between the four values is to a threshold L2Comparing the image intensity values with 0.200 to judge whether the uniformity of the image intensity is qualified or not, as follows;
Figure BDA0002852373500000063
Figure BDA0002852373500000064
computing
Figure BDA0002852373500000065
Arithmetic mean of
Figure BDA0002852373500000066
Respectively calculate
Figure BDA0002852373500000067
Figure BDA0002852373500000068
And
Figure BDA0002852373500000069
so that the absolute value of the largest difference is equal to
Figure BDA00028523735000000610
Is divided to obtain L10.038 is the light intensity distribution uniformity parameter of the image, since L1<L2Judging that the uniformity of the light intensity distribution of the image is qualified;
step three: and judging the definition of the concentric interference circular ring image. The full width at half maximum FWHM of the interference fringe, i.e. the phase difference corresponding to when the light intensity falls to half of the maximum. And moving the area array device near the ideal focal plane, receiving the concentric interference ring images of the area array device at different positions, and calculating the FWHM values of the rings with the assigned serial numbers in all the concentric interference ring images, wherein the FWHM values are respectively FWHMx +, FWHMx-, FWHMy + and FWHMy-. By comparing the FWHM arithmetic mean size of all the images, the image with the smallest FWHM arithmetic mean is the closest to the ideal focal plane, thereby obtaining the clearest image, and FWHMx +, FWHMx-, FWHMy +, FWHMy-, and FWHMy-of the specified 10 th circle are calculated for 11 focus positions as shown in table 1 below
Figure BDA0002852373500000071
TABLE 1 FWHM values at different focus positions
w
Figure BDA0002852373500000072
FWHMmin=min(FWHM1,FWHM2,FWHM3,FWHM4…FWHM10)=5.681
FWHMmin 5.681 is the concentric interference ring image minimum for these 11 focus positions
Figure BDA0002852373500000073
The value is the definition parameter of the image, and at this time, the corresponding n is 6 and is the serial number of the optimal focusing position, so that the clearest image can be obtained;
step four: and judging the ovalization degree of the concentric interference circular ring image. Calculating the FWHM values of the designated sequence number rings passing through the vicinity of the two central axes in the concentric interference ring image, wherein the FWHM values are respectively FWHMx +, (zero mean square root and zero mean square root),FWHMx-, FWHMy +, FWHMy-. By comparing the closeness between four values FWHMx +, FWHMx-, FWHMy +, FWHMy-with a threshold T2And (5) comparing the two values to be 0.200, and judging whether the ovalization of the circular ring is qualified or not, wherein the ovalization is as follows:
Figure BDA0002852373500000081
Figure BDA0002852373500000082
calculating the arithmetic mean of FWHMx +, FWHMx-, FWHMy +, FWHMy-)
Figure BDA0002852373500000083
FWHMx +, FWHMx-, FWHMy +, FWHMy-and
Figure BDA0002852373500000084
so that the absolute value of the largest difference is equal to
Figure BDA0002852373500000085
Is divided to obtain T1Is 0.065, i.e. the ovalization parameter of the image, due to T1<T2And judging that the ovalization degree of the image is qualified.
Step five: and judging the optimal selection among the imaging rings of the concentric interference ring image. Respectively establishing an x-axis rectangular coordinate system and a y-axis rectangular coordinate system in the horizontal direction and the vertical direction of the diameter of the approximate concentric interference circular ring, rotating the coordinate systems by 45 degrees anticlockwise to obtain the x '-axis rectangular coordinate system and the y' -axis rectangular coordinate system, and adopting a method for interpolating, subdividing and smoothing the array virtual pixels, which is disclosed in the patent application specification with the application date of 201710374595.0. And respectively establishing N parallel lines on two sides of the approximate diameter of the circular ring of the coordinate axis of the x 'axis and the y' axis, wherein N is 1. For the same ring, after the parallel line intersects with each ring, there are 8N + 4-12 photoelectric signal small line segments, and the method for solving the peak coordinate value of the photoelectric signal is adopted in the patent application specification with application number 201510217472.7,and solving the peak position coordinates and the standard deviation S of the peak position coordinates of the photoelectric signals on each small line segment of all the rings in the concentric interference ring image. Calculating the standard deviation arithmetic mean value of the peak position coordinates of 8N + 4-12 photoelectric signal small line segments in each circle
Figure BDA0002852373500000091
Comparing all rings
Figure BDA0002852373500000092
Value, select minimum
Figure BDA0002852373500000093
The circle where the value is located is used as the optimal circle for subsequently calculating the radius of the circle center, and the method comprises the following steps:
Figure BDA0002852373500000094
Figure BDA0002852373500000095
m is the peak position coordinate on a circle and the number of standard deviations 12. k is the serial number of the ring, 3 rings are provided in total,
Figure BDA0002852373500000096
and calculating the arithmetic mean value of all the peak position coordinate standard deviations on the circular ring where the designated circular ring serial number is located in the clearest concentric interference circular ring image.
Figure BDA0002852373500000097
For 3 rings in the clearest concentric interference ring image
Figure BDA0002852373500000098
The optimal selection parameter between the imaging rings of the image is obtained, and the corresponding k is 2 at the moment and is the optimal ring serial number for calculating the subsequent circle center displacement.
Finally, it should be noted that: the above-mentioned embodiments are only specific embodiments of the present invention, which are used for illustrating the technical solutions of the present invention and not for limiting the same, and the protection scope of the present invention is not limited thereto, although the present invention is described in detail with reference to the foregoing embodiments, those skilled in the art should understand that: any person skilled in the art can modify or easily conceive the technical solutions described in the foregoing embodiments or equivalent substitutes for some technical features within the technical scope of the present disclosure; such modifications, changes or substitutions do not depart from the spirit and scope of the embodiments of the present invention, and they should be construed as being included therein. Therefore, the protection scope of the present invention shall be subject to the protection scope of the appended claims.

Claims (1)

1. An F-P etalon interference imaging quality evaluation method for micro-angle measurement comprises the following steps:
the method comprises the following steps: and acquiring a concentric interference circular ring image. Receiving a concentric interference circular ring image before a deflection angle is generated by an area array device;
step two: and judging the light intensity uniformity of the concentric interference circular ring image. Taking the obtained concentric interference ring image as central axes along the horizontal direction and the vertical direction, measuring the maximum light intensity of all rings on the two central axes, establishing a rectangular coordinate system, taking the horizontal central axis as an x axis, the vertical central axis as a y axis, and the intersection point of the two central axes as an origin, respectively calculating the average value of the maximum light intensity value of each ring in the positive direction of the x axis, the negative direction of the x axis, the positive direction of the y axis and the negative direction of the y axis, respectively
Figure FDA0002852373490000019
By comparison
Figure FDA00028523734900000110
The proximity between the four values is to a threshold L2Comparing, and judging whether the uniformity of the image light intensity is qualified or not, as follows;
Figure FDA0002852373490000011
Figure FDA0002852373490000012
computing
Figure FDA0002852373490000013
Arithmetic mean of
Figure FDA0002852373490000014
Respectively calculate
Figure FDA0002852373490000015
Figure FDA0002852373490000016
And
Figure FDA0002852373490000017
so that the absolute value of the largest difference is equal to
Figure FDA0002852373490000018
Is divided to obtain L1I.e. the light intensity distribution uniformity parameter of the image, when L1>L2And judging that the uniformity of the light intensity distribution of the image is unqualified, otherwise, judging that the uniformity is qualified.
Step three: and judging the definition of the concentric interference circular ring image. The full width at half maximum FWHM of the interference fringe, i.e. the phase difference corresponding to when the light intensity falls to half of the maximum. And moving the area array device near the ideal focal plane, receiving the concentric interference ring images of the area array device at different positions, and calculating the FWHM values of the rings with the assigned serial numbers in all the concentric interference ring images, wherein the FWHM values are respectively FWHMx +, FWHMx-, FWHMy + and FWHMy-. Comparing the FWHM arithmetic mean values of all the images, wherein the image with the minimum FWHM arithmetic mean value is the closest to an ideal focal plane, so as to obtain the clearest image, and the method comprises the following steps;
FWHMn=(FWHMx+n)+(FWHMx-n)+(FWHMy+n)+(FWHMy-n)(n=1,2,3…)
FWHMmin=min(FWHM1,FWHM2,FWHM3,FWHM4…FWHMn)
the n serial numbers are used for marking different focusing positions, and the FWHMn is an arithmetic mean value of FWHM values of the serial number rings passing through the vicinity of two central axes in the concentric interference ring image at the same focusing position. FWHMmin is the minimum FWHMn value of concentric interference circular ring images at different focusing positions, namely the definition parameter of the images, and at the moment, the corresponding n is the serial number of the optimal focusing position, so that the clearest images can be obtained.
Step four: and judging the ovalization degree of the concentric interference circular ring image. And calculating the FWHM values of the designated sequence number rings passing through the vicinity of the two central axes in the concentric interference ring image, namely FWHMx +, FWHMx-, FWHMy + and FWHMy-. By comparing the closeness between four values FWHMx +, FWHMx-, FWHMy +, FWHMy-with a threshold T2And comparing to judge whether the ovalization of the ring is qualified or not, and the method comprises the following steps:
Figure FDA0002852373490000021
Figure FDA0002852373490000022
calculating the arithmetic mean of FWHMx +, FWHMx-, FWHMy +, FWHMy-)
Figure FDA0002852373490000031
FWHMx +, FWHMx-, FWHMy +, FWHMy-and
Figure FDA0002852373490000032
so that the absolute value of the largest difference is equal to
Figure FDA0002852373490000033
Is divided to obtain T1I.e. the ovalization parameter of the image, when T1>T2And judging that the ovalization degree of the image is unqualified, otherwise, judging that the ovalization degree is qualified.
Step five: and judging the optimal selection among the imaging rings of the concentric interference ring image. Respectively establishing an x-axis rectangular coordinate system and a y-axis rectangular coordinate system in the horizontal direction and the vertical direction of the diameter of the approximate concentric interference circular ring, rotating the coordinate systems by 45 degrees anticlockwise to obtain the x '-axis rectangular coordinate system and the y' -axis rectangular coordinate system, and adopting a method for interpolating, subdividing and smoothing the array virtual pixels, which is disclosed in the patent application specification with the application date of 201710374595.0. And respectively establishing N parallel lines on two sides of the approximate diameter of the circular ring of the coordinate axis of the x 'axis and the y' axis, wherein N is a positive integer. For the same ring, after the parallel line intersects with each ring, there are 8N +4 photoelectric signal small line segments, and the peak position coordinates and the standard deviation S of the photoelectric signals on each small line segment of all the rings in the concentric interference ring image are obtained by adopting the method for obtaining the peak coordinate values of the photoelectric signals in the patent application specification with the application number of 201510217472.7. Calculating the arithmetic mean value of the standard deviation of the peak position coordinates of 8N +4 small segments of the photoelectric signals in each circle
Figure FDA00028523734900000310
Comparing all rings
Figure FDA0002852373490000038
Value, select minimum
Figure FDA0002852373490000039
The circle where the value is located is used as the optimal circle for subsequently calculating the radius of the circle center, and the method comprises the following steps:
Figure FDA0002852373490000034
Figure FDA0002852373490000035
and m is the peak position coordinate and the standard deviation number on one circular ring. k is the serial number of the circular ring,
Figure FDA0002852373490000036
and calculating the arithmetic mean value of all the peak position coordinate standard deviations on the circular ring where the designated circular ring serial number is located in the clearest concentric interference circular ring image.
Figure FDA0002852373490000041
For each ring in the image of the most clearly concentric interfering rings
Figure FDA0002852373490000042
And the minimum value is the optimal selected parameter between the imaging rings of the image, and the corresponding k is the optimal ring serial number at the moment and is used for calculating the subsequent circle center displacement.
CN202011533267.9A 2020-12-23 2020-12-23 F-P etalon interference imaging quality evaluation method for micro angle measurement Active CN112729168B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011533267.9A CN112729168B (en) 2020-12-23 2020-12-23 F-P etalon interference imaging quality evaluation method for micro angle measurement

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011533267.9A CN112729168B (en) 2020-12-23 2020-12-23 F-P etalon interference imaging quality evaluation method for micro angle measurement

Publications (2)

Publication Number Publication Date
CN112729168A true CN112729168A (en) 2021-04-30
CN112729168B CN112729168B (en) 2022-11-22

Family

ID=75604147

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011533267.9A Active CN112729168B (en) 2020-12-23 2020-12-23 F-P etalon interference imaging quality evaluation method for micro angle measurement

Country Status (1)

Country Link
CN (1) CN112729168B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070127245A1 (en) * 2005-12-06 2007-06-07 3M Innovative Properties Company Illumination system incorporating collimated light source
WO2014121181A1 (en) * 2013-02-01 2014-08-07 The General Hospital Corporation Synthetic apertures for depth-of-focus tomography imaging
WO2015040618A1 (en) * 2013-09-17 2015-03-26 Ramot At Tel-Aviv University Ltd. A system and a method for quantitative sample imaging using off-axis interferometry with extended field of view or faster frame rate
TWI596446B (en) * 2016-04-08 2017-08-21 Luminous uniformity correction equipment and its correction method
CN113177935A (en) * 2021-05-21 2021-07-27 陕西利丰恒信生物科技发展有限公司 Near-infrared light intensity uniformity detection method and device and computer equipment

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070127245A1 (en) * 2005-12-06 2007-06-07 3M Innovative Properties Company Illumination system incorporating collimated light source
WO2014121181A1 (en) * 2013-02-01 2014-08-07 The General Hospital Corporation Synthetic apertures for depth-of-focus tomography imaging
WO2015040618A1 (en) * 2013-09-17 2015-03-26 Ramot At Tel-Aviv University Ltd. A system and a method for quantitative sample imaging using off-axis interferometry with extended field of view or faster frame rate
TWI596446B (en) * 2016-04-08 2017-08-21 Luminous uniformity correction equipment and its correction method
CN113177935A (en) * 2021-05-21 2021-07-27 陕西利丰恒信生物科技发展有限公司 Near-infrared light intensity uniformity detection method and device and computer equipment

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
杨淼等: "无参考水下静态图像质量评价与仿真分析", 《系统仿真学报》 *
黎明等: "利用光强信息的结构光图像轮廓提取修正方法", 《光电工程》 *

Also Published As

Publication number Publication date
CN112729168B (en) 2022-11-22

Similar Documents

Publication Publication Date Title
US10412365B2 (en) Calibration arrangement for structured light system using a tele-centric lens
Palousek et al. Effect of matte coating on 3D optical measurement accuracy
CN110411346B (en) Method for quickly positioning surface micro-defects of aspheric fused quartz element
JP2008275612A (en) Device equipped with high resolution measurement structure on substrate for manufacturing semiconductor, and use of aperture in measuring device
CN104655051A (en) High-speed structured light three-dimensional surface shape vertical measurement method
CN110702009A (en) Three-dimensional measurement system based on reverse Hartmann computer-aided method
CN109238659A (en) A kind of focal length of lens measurement technology and device based on experiment ray tracing principle
CN110702026A (en) Flatness three-dimensional shape detection device based on complex beam angle adaptive optics and processing method thereof
CN107063644B (en) Finite object distance distortion measuring method and system
CN112066913A (en) Method for measuring spherical curvature radius
CN113251953B (en) Mirror included angle measuring device and method based on stereo deflection technology
Ma et al. Non-diffracting beam based probe technology for measuring coordinates of hidden parts
CN112729168B (en) F-P etalon interference imaging quality evaluation method for micro angle measurement
CN111366079B (en) Calibration method for rapidly determining geometric positions of components in deflection measurement system
CN110428471B (en) Accurate self-positioning method for optical free-form surface sub-aperture deflection measurement
TW201641914A (en) Full-range image detecting system and method thereof
CN105547653A (en) Double-telecentric-lens multiplying power measurement method based on More interference
KR101833245B1 (en) METHOD AND SYSTEM FOR CALIBRATION OF SURFACE MEASUREMENT DEVICE USING Phase measuring deflectometry
CN112781502B (en) Method for verifying and adjusting grating parallelism based on Moire fringe equation
CN114659455A (en) Diffraction stripe recognition processing method for measuring diameter of filament
CN115164775A (en) Large-caliber convex aspheric reflector surface shape detection device and detection method thereof
KR20150029424A (en) Appapatus for three-dimensional shape measurment and method the same
Zhao et al. Relationship between the line structured light vision calibration accuracy and hardware parameters
Zhao et al. Study on the Technologies of Close Range Photogrammetry and Applications in the Manufacture of Aviation
CN107941147B (en) Non-contact online measurement method for three-dimensional coordinates of large-scale system

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