CN103489186A - Spatial position matching method of dynamic interferometer child interferograms - Google Patents

Spatial position matching method of dynamic interferometer child interferograms Download PDF

Info

Publication number
CN103489186A
CN103489186A CN201310422931.6A CN201310422931A CN103489186A CN 103489186 A CN103489186 A CN 103489186A CN 201310422931 A CN201310422931 A CN 201310422931A CN 103489186 A CN103489186 A CN 103489186A
Authority
CN
China
Prior art keywords
interferogram
subimage
sub
child
neutron
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201310422931.6A
Other languages
Chinese (zh)
Inventor
陈磊
郑东晖
周雅娴
陈宇琪
周舒
李金鹏
李建欣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Science and Technology
Original Assignee
Nanjing University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN201310422931.6A priority Critical patent/CN103489186A/en
Publication of CN103489186A publication Critical patent/CN103489186A/en
Pending legal-status Critical Current

Links

Images

Abstract

The invention discloses a spatial position matching method of dynamic interferometer child interferograms. According to the method, a phase correlation method is used for determining spatial position matching errors of the child interferograms of a dynamic interferometer. The method comprises the following steps that spatial phase shifting interferograms of the dynamic interferometer are collected, the child interferograms of the spatial phase shifting interferograms are divided into child images with the same size, each child image comprises a whole child interferogram at a corresponding position; one of the child images is used as a standard, the phase correlation method is used for determining spatial position matching errors of the interferograms in other child images and the child interferogram in the standard child image; according to the spatial position matching errors, the phase position formula is corrected, and spatial position matching of the child interferograms is completed. The method improves the automation level of dynamic interferometer spatial position matching, lowers influence on spatial position matching caused by light split ratio, image gradation and other factors, and has strong applicability.

Description

The locus matching process of the sub-interferogram of dynamic interferometer
Technical field
The invention belongs to interference of light metering field, particularly the locus matching process of the sub-interferogram of dynamic interferometer.
Background technology
Thereby dynamic interferometer is widely used among the on-line measurement of optics processing site due to its good resistance to shock.Yet no matter present stage is to adopt the many CCD of single CCD to gather the dynamic interferometer of interferogram, although light splitting means different (prismatic decomposition, one-dimensional grating light splitting, two-dimensional quadrature grating beam splitting) but all have the locus matching problem of phase-shift interference.The locus method for registering adopted at present generally can be divided into following two kinds:
A kind of is the space-artifact position registration.This scheme often adds certain sharp keen target (as needle point) in the light path before light splitting, due to each image device after light splitting, all to this target imaging, the picture of this target is arranged in sub-interferogram.Will be in four width interferograms take out separately, artificially identify the most advanced and sophisticated of target picture and read this most advanced and sophisticated position coordinates in sub-interferogram, because these coordinate points are same points on space actually, utilize these coordinates to carry out sub-interferogram locus coupling.The method does not need to carry out the consideration on algorithm, but higher to sharp keen target imaging quality requirements, and the identification of the sharp keen target in four width interferograms is comprised to human factor, and automaticity is not high.
Another kind is to adopt the method for image correlation to carry out spatial registration.This scheme generally first adopts certain means to make the same image of sub-hot spot imaging, be that sub-hot spot has identical gray level image to distribute, adopt the image correlation parameter to mate evaluating as locus, by the related operation of image, realize the locus coupling between each subgraph.The method depends on gradation of image, and reliability easily is subject to the interference of image quality, light splitting ratio error factors.
Summary of the invention
The object of the present invention is to provide the locus matching process of the sub-interferogram of dynamic interferometer that a kind of precision is high, reliability is strong, convenient, as to realize efficiently dynamic interferometer neutron interferogram locus coupling.
The technical scheme that realizes the object of the invention is: the locus matching process of the sub-interferogram of a kind of dynamic interferometer comprises the following steps:
Step 1, gather the spatial Phase-shifting Method interferogram of dynamic interferometer, and each the sub-interferogram in the partition space phase-shift interference be the subimage that size is identical, and each sub-picture pack contains the whole sub-interferogram of correspondence position;
Step 2, the wherein width subimage of usining, as benchmark, adopts phase correlation method to determine the locus matching error of other subimage neutron interferogram and benchmark subimage neutron interferogram;
Step 3, according to the locus matching error, revise the locus coupling that the phase place formula completes each sub-interferogram.
Compared with prior art, its remarkable advantage is in the present invention: (1) has largely improved the automaticity of dynamic interferometer locus coupling, eliminates the impact that human factor causes the locus coupling; (2) reduce the impact that the factors such as light splitting proportional error, gradation of image cause the locus coupling, improved reliability; (3) method is quick and easy, and not additional any ancillary hardware, is applicable to the dynamic interferometer that great majority adopt the light splitting scheme.
Below in conjunction with accompanying drawing, the present invention is described in further detail.
The accompanying drawing explanation
Fig. 1 is the hardware experiments device schematic diagram of the locus matching process of the sub-interferogram of dynamic interferometer of the present invention.
Fig. 2 is the spatial Phase-shifting Method interferogram that dynamic interferometer of the present invention records, and comprises A, B, C, tetra-sub-interferograms of D.
Fig. 3 is the phase-shift interference of Spatial Mismatch, and wherein (a), (b), for two width subimages after cutting apart, (c) are the coordinate system of image (a), (b).
Fig. 4 is the dynamic interferometer movable phase interfere map space matching result based on grating beam splitting in embodiment 1.
Embodiment
Below in conjunction with drawings and the specific embodiments, the present invention is described in further detail.
Dynamic interferometer can collect several width simultaneously and have the fixedly sub-interferogram of amount of phase shift, thereby organizes sub-interferogram and can resolve the face shape that PHASE DISTRIBUTION to be measured is determined measured piece from this.Adopt the dynamic interferometer of light splitting scheme generally to collect the mutually sub-interferogram of four amplitude shifts, be respectively I 1(x, y), I 2(x, y), I 3(x, y), I 4(x, y), the phase shift phasor magnitude is pi/2, if not Existential Space location matches error, tested phase place can be solved by following formula:
If yet Existential Space location matches error is calculated phase place by (1) formula and will be introduced measuring error.
In conjunction with Fig. 1~4, the locus matching process of the sub-interferogram of dynamic interferometer of the present invention, utilize the locus of the sub-interferogram of phase place related algorithm implementation space phase shift to mate, and comprises the following steps:
Step 1, gather the spatial Phase-shifting Method interferogram of dynamic interferometer, and each the sub-interferogram in the partition space phase-shift interference be the subimage that size is identical, and each sub-picture pack contains the whole sub-interferogram of correspondence position, is specially:
In conjunction with Fig. 1, place standard spherical mirror in the dynamic interferometer test arm, make the axial distance of standard spherical mirror and dynamic interferometer equal the radius-of-curvature of standard spherical mirror, at this moment can collect the spatial Phase-shifting Method interferogram by CCD, the inclination of adjusting axial defocusing and standard spherical mirror makes the closed fringe number of each sub-interferogram all be less than 3, and closed center all is positioned at the field of view center of corresponding sub-interferogram; Gather spatial Phase-shifting Method interferogram now, wherein comprise some sub-interferograms, if take Fig. 2 as example, wherein comprise A, B, C, D tetra-width interferograms, as shown in Figure 3, the sub-interferogram of each in the partition space phase-shift interference is the subimage that size is identical, each sub-picture pack is containing the whole sub-interferogram of correspondence position, Fig. 3 (a), 3 (b) are two width subimages after cutting apart, and in figure, the white edge scope comprises corresponding locus, and does not mate the locus of known two width interferograms.
Step 2, the wherein width subimage of usining, as benchmark, adopts phase correlation method to determine the locus matching error of other subimage neutron interferogram and benchmark subimage neutron interferogram, and detailed process is as follows:
(1) choose wherein a width subimage as benchmark subimage f 0(x, y), by benchmark subimage f 0(x, y) and another width subimage f 1(x, y) carries out Fourier transform under as Fig. 3 (c) at the same coordinate system:
F 0 ( u , v ) = FFT { f 0 ( x , y ) } F 1 ( u , v ) = FFT { f 1 ( x , y ) } - - - ( 2 )
(2) determine F 0(u, v) and F 1the phase place correlation matrix of (u, v), and the phase place correlation matrix is obtained to Spatial Mismatch function C (x, y) do inverse Fourier transform:
C ( x , y ) = FFT - 1 { F 0 ( u , v ) * F 1 ( u , v ) | F 0 ( u , v ) | | F 1 ( u , v ) | } - - - ( 3 )
In formula *f 1(u, v) is F 1the conjugation of (u, v);
(3) find coordinate (Δ x, Δ y) corresponding to Spatial Mismatch function C (x, y) maximal value, be subimage f 1(x, y) neutron interferogram and benchmark subimage f 0the locus matching error of (x, y) neutron interferogram:
(Δx,Δy)=Findcoordinate{max{C(x,y)}} (4)
In formula, Findcoordinate means to find coordinate computation;
(4), according to the method for (1)~(3), determine successively the locus matching error of other subimage neutron interferogram and benchmark subimage neutron interferogram.
Step 3, according to the locus matching error, revise the locus coupling that the phase place formula completes each sub-interferogram.Take that to collect the mutually sub-interferogram of four amplitude shifts be example, establish subimage f 1(x, y), f 2(x, y), f 3(x, y) neutron interferogram and benchmark subimage f 0the locus matching error of (x, y) neutron interferogram is respectively (Δ x 1, Δ y 1), (Δ x 2, Δ y 2), (Δ x 3, Δ y 3), the tested phase place of revising
Figure BDA0000382743230000041
can be solved by following formula:
Figure BDA0000382743230000042
Below in conjunction with specific embodiment, the present invention is described in further detail.
Embodiment 1
The dynamic interferometer of utilization based on grating beam splitting realized the locus coupling of sub-interferogram, and this interferometer utilizes a two-dimensional quadrature grating to generate the spatial Phase-shifting Method interferogram of 2 * 2 arrays.
Step 1, used the proving installation shown in Fig. 1, and wherein dynamic interferometer is the dynamic interferometer based on grating beam splitting, gathers the spatial Phase-shifting Method interferogram, as shown in Figure 4 interferogram is divided into to four number of sub images, number consecutively 0,1,2,3.
Step 2, using subimage 0 as benchmark, adopt phase correlation method to determine that subimage 1,2,3 neutron interferograms are respectively (26 with respect to the locus matching error of benchmark subimage 0 neutron interferogram,-20), (12 ,-17), (17,19), Fig. 4 has shown the space matching result of interferogram 0 with interferogram 1,2,3.
Step 3, according to the locus matching error, revise the locus coupling that the phase place formula completes each sub-interferogram.Phase place formula to be measured is revised:
Figure BDA0000382743230000043
In sum, the locus matching process of the sub-interferogram of dynamic interferometer of the present invention, utilize the locus coupling of the sub-interferogram of phase place related algorithm implementation space phase shift, largely improve the automaticity of dynamic interferometer locus coupling, eliminated the impact that human factor causes the locus coupling; Reduce the impact that the factors such as light splitting proportional error, gradation of image cause the locus coupling, improved reliability; The method is quick and easy, and not additional any ancillary hardware, is applicable to the dynamic interferometer that great majority adopt the light splitting scheme.

Claims (3)

1. the locus matching process of the sub-interferogram of dynamic interferometer, is characterized in that, comprises the following steps:
Step 1, gather the spatial Phase-shifting Method interferogram of dynamic interferometer, and each the sub-interferogram in the partition space phase-shift interference be the subimage that size is identical, and each sub-picture pack contains the whole sub-interferogram of correspondence position;
Step 2, the wherein width subimage of usining, as benchmark, adopts phase correlation method to determine the locus matching error of other subimage neutron interferogram and benchmark subimage neutron interferogram;
Step 3, according to the locus matching error, revise the locus coupling that the phase place formula completes each sub-interferogram.
2. the locus matching process of the sub-interferogram of dynamic interferometer according to claim 1, is characterized in that, the closed fringe number of described each the sub-interferogram of step 1 all is less than 3, and closed center all is positioned at the field of view center of corresponding sub-interferogram.
3. the locus matching process of the sub-interferogram of dynamic interferometer according to claim 1, it is characterized in that, the described employing phase correlation method of step 2 is determined the locus matching error of other subimage neutron interferogram and benchmark subimage neutron interferogram, and detailed process is as follows:
(1) choose wherein a width subimage as benchmark subimage f 0(x, y), by benchmark subimage f 0(x, y) and another width subimage f 1(x, y) carries out Fourier transform under the same coordinate system:
F 0(u,v)=FFT{f 0(x,y)}
F 1(u,v)=FFT{f 1(x,y)}
(2) determine F 0(u, v) and F 1the phase place correlation matrix of (u, v), and the phase place correlation matrix is obtained to Spatial Mismatch function C (x, y) do inverse Fourier transform:
C ( x , y ) = FFT - 1 { F 0 ( u , v ) * F 1 ( u , v ) | F 0 ( u , v ) | | F 1 ( u , v ) | }
In formula *f 1(u, v) is F 1the conjugation of (u, v);
(3) find coordinate (Δ x, Δ y) corresponding to Spatial Mismatch function C (x, y) maximal value, be subimage f 1(x, y) neutron interferogram and benchmark subimage f 0the locus matching error of (x, y) neutron interferogram:
(Δx,Δy)=Findcoordinate{max{C(x,y)}}
In formula, Findcoordinate means to find coordinate computation;
(4), according to the method for (1)~(3), determine successively the locus matching error of other subimage neutron interferogram and benchmark subimage neutron interferogram.
CN201310422931.6A 2013-09-16 2013-09-16 Spatial position matching method of dynamic interferometer child interferograms Pending CN103489186A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310422931.6A CN103489186A (en) 2013-09-16 2013-09-16 Spatial position matching method of dynamic interferometer child interferograms

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310422931.6A CN103489186A (en) 2013-09-16 2013-09-16 Spatial position matching method of dynamic interferometer child interferograms

Publications (1)

Publication Number Publication Date
CN103489186A true CN103489186A (en) 2014-01-01

Family

ID=49829384

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310422931.6A Pending CN103489186A (en) 2013-09-16 2013-09-16 Spatial position matching method of dynamic interferometer child interferograms

Country Status (1)

Country Link
CN (1) CN103489186A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107121205A (en) * 2017-05-05 2017-09-01 南京理工大学 A kind of spot light dislocation type Mach-Zehnder interferometers measurement apparatus and method

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5534704A (en) * 1993-04-21 1996-07-09 Sharp Kabushiki Kaisha Optical image correlator and system for performing parallel correlation
CN102507020A (en) * 2011-11-01 2012-06-20 南京理工大学 Microlens array-based synchronized phase-shifting interference test method and test device

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5534704A (en) * 1993-04-21 1996-07-09 Sharp Kabushiki Kaisha Optical image correlator and system for performing parallel correlation
CN102507020A (en) * 2011-11-01 2012-06-20 南京理工大学 Microlens array-based synchronized phase-shifting interference test method and test device

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
BO LI等: "Spatial mismatch calibration using circular carriercircular technique in the simultaneous phase shifting interferometry", 《APPLIED OPTICS》, vol. 51, no. 8, 31 December 2012 (2012-12-31), pages 1037 - 1044, XP001574463, DOI: doi:10.1364/AO.51.001037 *
孙辉: "相位相关技术实现离焦模糊图像运动估计", 《液晶与显示》, vol. 27, no. 2, 30 April 2012 (2012-04-30), pages 223 - 228 *
徐晨: "动态干涉测试技术与应用研究", 《中国博士学位论文全文数据库》, no. 01, 15 January 2011 (2011-01-15) *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107121205A (en) * 2017-05-05 2017-09-01 南京理工大学 A kind of spot light dislocation type Mach-Zehnder interferometers measurement apparatus and method

Similar Documents

Publication Publication Date Title
CN103713287B (en) A kind of height reconstruction method based on relatively prime many baselines and device
Jackson et al. LBCS: the LOFAR long-baseline calibrator survey
CN106526593B (en) Sub-pixel-level corner reflector automatic positioning method based on the tight imaging model of SAR
CN103592108A (en) CCD chip modulation transfer function test device and method
CN105842694A (en) FFBP SAR imaging-based autofocus method
CN108535097A (en) A kind of method of triaxial test sample cylindrical distortion measurement of full field
CN105606217A (en) Image-spectrum-polarization-state integrated obtaining apparatus and method
CN102662171A (en) Synthetic aperture radar (SAR) tomography three-dimensional imaging method
US10551177B2 (en) Apparatus and method for measuring 3D form or deformation of an object surface using a grid pattern and reference plane
CN104730520A (en) Circumference SAR back projection self-focusing method based on subaperture synthesis
CN102788979B (en) GPU (graphic processing unit) implementing method based on backward projection InSAR (interfering synthetic aperture radar) imaging registration
CN104748835A (en) Interference-mount separating type nonlinear error correcting method and device for laser interference vibration tester
GB2554796A (en) Testing 3D imaging systems
CN106249236A (en) A kind of spaceborne InSAR long-short baselines image associating method for registering
CN107783086B (en) Method for diagnosing distorted position of antenna array aperture amplitude phase field
Simon et al. CFHTLenS: higher order galaxy–mass correlations probed by galaxy–galaxy–galaxy lensing
CN105487065A (en) Time sequence satellite borne radar data processing method and device
CN104422916A (en) Two-point Dixon technique
KR100921368B1 (en) Enhanced sound source localization system and method by using a movable microphone array
CN102096920B (en) Target image-based sub-pixel registering method
CN108508468A (en) Medium-long baselines single epoch Ambiguity Solution Methods, system, equipment, storage medium
CN113607346A (en) Three-dimensional acoustic imaging real-time signal processing device based on embedded processor
CN103886600A (en) Continuous phase plate machining error identification method
CN104237876B (en) Metal pylon is to the determination of distance method and device of calibration support
CN103489186A (en) Spatial position matching method of dynamic interferometer child interferograms

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
C53 Correction of patent of invention or patent application
CB03 Change of inventor or designer information

Inventor after: Chen Lei

Inventor after: Zheng Donghui

Inventor after: Zhou Yaxian

Inventor after: Zhou Shu

Inventor after: Chen Yuqi

Inventor after: Li Jinpeng

Inventor after: Guo Renhui

Inventor after: Li Jianxin

Inventor before: Chen Lei

Inventor before: Zheng Donghui

Inventor before: Zhou Yaxian

Inventor before: Chen Yuqi

Inventor before: Zhou Shu

Inventor before: Li Jinpeng

Inventor before: Li Jianxin

COR Change of bibliographic data

Free format text: CORRECT: INVENTOR; FROM: CHEN LEI ZHENG DONGHUI ZHOU YAXIAN CHEN YUQI ZHOU SHU LI JINPENG LI JIANXIN TO: CHEN LEI ZHENG DONGHUI ZHOU YAXIAN ZHOU SHU CHEN YUQI LI JINPENG GUO RENHUI LI JIANXIN

SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20140101