CN106908148B - A kind of infrared interference figure method for reconstructing - Google Patents
A kind of infrared interference figure method for reconstructing Download PDFInfo
- Publication number
- CN106908148B CN106908148B CN201710128061.XA CN201710128061A CN106908148B CN 106908148 B CN106908148 B CN 106908148B CN 201710128061 A CN201710128061 A CN 201710128061A CN 106908148 B CN106908148 B CN 106908148B
- Authority
- CN
- China
- Prior art keywords
- interference signal
- sampling
- infrared
- interpolation
- laser
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000005070 sampling Methods 0.000 claims abstract description 60
- 238000012952 Resampling Methods 0.000 claims abstract description 8
- 230000003287 optical effect Effects 0.000 claims description 7
- 238000001228 spectrum Methods 0.000 claims description 4
- 238000000605 extraction Methods 0.000 abstract description 2
- 230000001360 synchronised effect Effects 0.000 abstract 1
- 230000002123 temporal effect Effects 0.000 abstract 1
- 239000011521 glass Substances 0.000 description 4
- 108010074864 Factor XI Proteins 0.000 description 2
- 238000001157 Fourier transform infrared spectrum Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J3/45—Interferometric spectrometry
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Spectrometry And Color Measurement (AREA)
Abstract
The invention discloses a kind of infrared interference figure method for reconstructing, includes the laser interference signal and infrared interference signal that Michelson's interferometer light splitting generates, use identical sample rate while over-sampling laser interference signal and infrared interference signal;According to the needs of actual conditions, identical interpolation factor is set, interpolation operation is carried out to the two paths of signals after over-sampling;Laser zero crossing information is obtained in laser interference signal after interpolation, and resampling is then carried out to the infrared interference signal after also passing through interpolation according to this zero crossing information, finally obtains the infrared interference figure at the aplanatism difference interval after rebuilding.The present invention is synchronous with interpolation operation to the sampling of two paths of signals, ensure that temporal consistency, simplifies the extraction of laser interference signal zero-crossing information;The complexity of data processing is reduced using method disclosed by the invention.
Description
Technical field
The invention belongs to infrared signal acquisition and process fields, and in particular to a kind of infrared interference figure method for reconstructing.
Background technology
In FTIR spectrum instrument, the light that light source is sent out is divided into two bundles by the beam splitter in Michelson's interferometer, wherein one
Beam is transmitted to interferometer index glass through beam splitter reflection to interferometer horizontal glass, another beam through beam splitter.Then respectively along respective original
It light path reflected back into beam splitter and is divided again by beam splitter, finally generates interference signal.Since the sampling location of interference signal is to dry
The movement velocity of index glass is very sensitive in interferometer, and the minor change of movement velocity can all cause interference pattern that severe deviations occur, from
And influence to restore the precision of spectrum.Therefore, usually using in Michelson's interferometer index glass move caused by optical path difference come etc.
Interval sampling interference signal.It is used as in this process usually using an one-wavelength laser and carrys out measurement of optical path difference with reference to light source.Mesh
Before, main Model of Interferogram Sampling method has laser triggering method and " Brault " method of sampling in FTIR spectrum instrument.Laser triggering
For method due to the limitation of nyquist sampling theorem, sample frequency is relatively low and to the more demanding of instrument and hardware.And
The method of sampling that " Brault " is proposed uses nonsynchronous counter records laser and infrared signal, and nonsynchronous sampling is two
Delay is increased in the signal of road, causes to need to be corrected this delay in subsequent processing, and used hardware in the method, is increased
The complexity for having added apparatus structure reduces the adaptedness to ambient enviroment.Therefore it needs to study a kind of good environmental adaptability,
The higher interference pattern method for reconstructing of sample frequency.
Invention content
The object of the present invention is to provide one kind capable of improving instrument signal to noise ratio, reduces the infrared interference of data processing complexity
Figure method for reconstructing.
Brault in 1996 has been put forward for the first time using nonsynchronous separate counters recording laser and infrared interference signal, this
Method sampling process is fairly simple, and data volume to be treated is bigger, and this method can not only improve instrument signal to noise ratio, also
Sampling error can be reduced.But which employs nonsynchronous counters, it cannot be guaranteed that laser and infrared interference that sampling obtains
The consistency of signal in time, and the hardware used is higher to environmental requirement, apparatus structure is complex, and needs to calculate
The zero crossing information of laser interference signal, calculating complex and processing time is longer, needs research one kind that can simplify interference
Signal sampling and processing, ensure the interference pattern method for reconstructing for restoring spectral accuracy.
To achieve the above object, the present invention provides following technical schemes:
A kind of infrared interference figure method for reconstructing, it is characterised in that:Include the laser that Michelson's interferometer light splitting generates
Interference signal and infrared interference signal, wherein laser interference signal are reference signal for generating zero crossing information;Using identical
Over-sampling laser interference signal and infrared interference signal, over-sampling are used to improve instrument signal to noise ratio sample rate simultaneously;According to reality
The needs of situation are arranged identical interpolation factor and carry out interpolation operation to the two paths of signals after over-sampling, and interpolation is for keeping
Additional data point is generated on the basis of signal waveform, ensures the accuracy for restoring spectrum;Laser interference signal after interpolation
Then middle acquisition laser zero crossing information carries out weight according to this zero crossing information to the infrared interference signal after also passing through interpolation
Sampling finally obtains the infrared interference figure at the aplanatism difference interval after rebuilding.
The infrared interference figure method for reconstructing, it is characterised in that specific method includes the following steps:
(1) according to nyquist sampling theorem, sample frequency is at least 2 times of interference signal maximum frequency to be measured, uses
Identical sample rate carries out over-sampling to laser interference signal and infrared interference signal simultaneously, i.e.,
fs≥2fm (1)
The two-way interference signal obtained at this time is the function of time, and dt is the time interval of sampling, ti=0, dt, 2dt,
3dt ... is the sampling instant corresponding to each laser interference signal sampling point and infrared interference signal sampling point, then infrared interference
Signal is identical in the sampling instant of ith sample point as laser interference signal in the sampling instant of ith sample point, i.e.,
tii=tli(2);
(2) and then according to actual needs suitable interpolation factor two paths of signals is selected to be carried out at the same time interpolation operation, using slotting
The laser interference signal and infrared interference signal that value function obtains over-sampling generate additional additional data point into row interpolation;
(3) the laser interference signal data point after over-sampling and interpolation and infrared interference signal data point data amount compared with
Greatly, consider and calculate time and computational accuracy, the value nearest apart from zero is indexed directly in laser interference signal data point
Corresponding location information, it is final to obtain then according to the infrared interference signal after this location information resampling over-sampling and interpolation
To the infrared interference figure at aplanatism difference interval, i.e., infrared interference figure at this time is the function of optical path difference.
A kind of infrared interference figure method for reconstructing, it is characterised in that:The laser interference signal is by wavelength
632.8nm, the He-Ne laser light sources that modulating frequency is 5KHz generate, and infrared interference signal is generated by globar infrared light supply.
A kind of infrared interference figure method for reconstructing, it is characterised in that:The sample rate chooses interference signal to be measured most
5~10 times of big frequency.
A kind of infrared interference figure method for reconstructing, it is characterised in that:The sample rate chooses 60KS/s.
A kind of infrared interference figure method for reconstructing, it is characterised in that:The interpolation factor is set as 20.
A kind of infrared interference figure method for reconstructing, it is characterised in that:The interpolating function is using Fourier's interpolation, line
Property interpolation or cubic spline functions.
A kind of infrared interference figure method for reconstructing, it is characterised in that:The resampling infrared interference figure is according to sharp
Light zero crossing location index one-to-one infrared interference signal strength therewith.
The technical effects of the invention are that:
By the present invention in that synchronizing over-sampling laser and infrared interference signal with same sample rate constant duration, phase is set
Same interpolation factor, is carried out at the same time interpolation to two paths of signals, ensure that the time consistency of two paths of signals, simplify reference laser
The extraction of zero crossing reduces the complexity of interference pattern data processing.
Description of the drawings
Fig. 1 is the infrared interference figure method for reconstructing schematic diagram of the present invention.
Fig. 2 is the infrared interference figure reconstruction process schematic diagram of the present invention.
Fig. 3 is the infrared interference figure at the aplanatism difference interval after rebuilding.
Fig. 4 is the infrared interference figure partial enlarged view at the aplanatism difference interval after rebuilding.
Specific implementation mode
As shown in Figure 1, a kind of infrared interference figure method for reconstructing, specifically includes:What Michelson's interferometer light splitting generated swashs
Optical interference signals 1 and infrared interference signal 2, wherein laser interference signal 1 are reference signal for generating zero crossing information 7;Make
With identical sample rate 10, over-sampling 3,4 laser interference signals 1 and infrared interference signal 2, over-sampling are believed for improving instrument simultaneously
It makes an uproar ratio;According to the needs of actual conditions, identical interpolation factor 11 is arranged to grasp the two paths of signals after over-sampling into row interpolation 5,6
Make, interpolation ensures the accuracy for restoring spectrum for generating additional data point on the basis of keeping signal waveform;In interpolation
Laser zero crossing information 7 is obtained in laser interference signal 1 afterwards, then according to this zero crossing information 7 to also passing through interpolation after
Infrared interference signal 2 carry out resampling 8, finally obtain rebuild after aplanatism difference interval infrared interference Fig. 9.
The He-Ne laser light sources that laser interference signal by wavelength is 632.8nm, modulating frequency is 5KHz generate, infrared dry
Signal is related to be generated by globar infrared light supply.
Embodiment
Such as Fig. 2-4, according to nyquist sampling theorem, sample frequency is at least 2 times of interference signal maximum frequency to be measured,
I.e.
fs≥2fm (1)
Sample frequency generally chooses 5~10 times of interference signal maximum frequency to be measured.Use 60KS/s's in the present embodiment
Sample rate 10 carries out over-sampling to laser interference signal 1 and infrared interference signal 2 simultaneously, and the two-way interference signal obtained at this time is
The function of time, dt are the time interval of sampling, ti=0, dt, 2dt, 3dt ... it is each 12 He of laser interference signal sampling point
Sampling instant corresponding to infrared interference signal sampling point 13, then infrared interference signal ith sample point sampling instant with
Laser interference signal ith sample point sampling instant, i.e.,
tii=tli (2)
Then it is needed to select suitable interpolation factor according to instrument, it is 20 pairs of two-way that interpolation factor 11 is arranged in the present embodiment
Signal is carried out at the same time interpolation (5 and 6) operation, can use Fourier's interpolation, linear interpolation, cubic spline interpolation equal interpolation order
The laser interference signal 12 and infrared interference signal 13 obtained to over-sampling generates additional additional data point into row interpolation.By
Oversampling technique has been used in front, therefore can suitably reduce interpolation factor to improve the arithmetic speed of computer.
15 data volume of laser interference signal data point 14 and infrared interference signal data point after over-sampling and interpolation
It is larger, consider and calculate time and computational accuracy, can be indexed apart from zero directly in laser interference signal data point 14
Location information 16 corresponding to 18 nearest values indexes | Ili|minCorresponding sampling instant tlzeroi;Then according to this position
Infrared interference signal 6 after information resampling over-sampling and interpolation, according to the sampling of laser and infrared interference signal during this
Moment corresponds, and therefore, need to only use the zero passage point moment t of laser interference signallzeroiIn resampling infrared interference signal with this
Moment corresponding infrared interference signal light intensity;Finally obtain the infrared interference figure at aplanatism difference interval, i.e. at this time infrared
Interference pattern is the function of optical path difference.
Claims (7)
1. a kind of infrared interference figure method for reconstructing, it is characterised in that:Include that the laser that generates of Michelson's interferometer light splitting is dry
Signal and infrared interference signal are related to, wherein laser interference signal is reference signal for generating zero crossing information;It is adopted using identical
Over-sampling laser interference signal and infrared interference signal, over-sampling are used to improve instrument signal to noise ratio sample rate simultaneously;According to practical feelings
The needs of condition are arranged identical interpolation factor and carry out interpolation operation to the two paths of signals after over-sampling, and interpolation is used to keep believing
Additional data point is generated on the basis of number waveform, ensures the accuracy for restoring spectrum;In laser interference signal after interpolation
Laser zero crossing information is obtained, then the infrared interference signal after also passing through interpolation is adopted again according to this zero crossing information
Sample finally obtains the infrared interference figure at the aplanatism difference interval after rebuilding;
Specific method includes the following steps:
(1)According to nyquist sampling theorem, sample frequency is at least 2 times of interference signal maximum frequency to be measured, and use is identical
Sample rate over-sampling is carried out to laser interference signal and infrared interference signal simultaneously, i.e.,
(1)
The two-way interference signal obtained at this time is the function of time,For the time interval of sampling, t i =0,dt,2dt, 3dt,…
For the sampling instant corresponding to each laser interference signal sampling point and infrared interference signal sampling point, then infrared interference signal exists
The sampling instant of ith sample point is identical in the sampling instant of ith sample point as laser interference signal, i.e.,
(2);
(2)Then it selects suitable interpolation factor two paths of signals to be carried out at the same time interpolation operation according to actual needs, uses interpolation letter
Several laser interference signals obtained to over-sampling and infrared interference signal generate additional additional data point into row interpolation;
(3)Laser interference signal data point and infrared interference signal data point data amount after over-sampling and interpolation is larger,
Consider and calculate time and computational accuracy, it is right that the value institute nearest apart from zero is indexed directly in laser interference signal data point
The location information answered finally obtains then according to the infrared interference signal after this location information resampling over-sampling and interpolation
The infrared interference figure at optical path difference interval, i.e., infrared interference figure at this time are the functions of optical path difference.
2. a kind of infrared interference figure method for reconstructing according to claim 1, it is characterised in that:The laser interference signal by
The He-Ne laser light sources that wavelength is 632.8nm, modulating frequency is 5KHz generate, and infrared interference signal is by globar infrared light supply
It generates.
3. a kind of infrared interference figure method for reconstructing according to claim 2, it is characterised in that:The sample rate is chosen to be measured
5 ~ 10 times of interference signal maximum frequency.
4. a kind of infrared interference figure method for reconstructing according to claim 3, it is characterised in that:The sample rate is chosen
60KS/s。
5. a kind of infrared interference figure method for reconstructing according to claim 1, it is characterised in that:The interpolation factor is set as
20。
6. a kind of infrared interference figure method for reconstructing according to claim 1, it is characterised in that:The interpolating function uses Fu
In leaf interpolation, linear interpolation or cubic spline functions.
7. a kind of infrared interference figure method for reconstructing according to claim 1, it is characterised in that:The resampling infrared interference
Figure is according to laser zero crossing location index one-to-one infrared interference signal strength therewith.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710128061.XA CN106908148B (en) | 2017-03-06 | 2017-03-06 | A kind of infrared interference figure method for reconstructing |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710128061.XA CN106908148B (en) | 2017-03-06 | 2017-03-06 | A kind of infrared interference figure method for reconstructing |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106908148A CN106908148A (en) | 2017-06-30 |
CN106908148B true CN106908148B (en) | 2018-08-24 |
Family
ID=59186084
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710128061.XA Active CN106908148B (en) | 2017-03-06 | 2017-03-06 | A kind of infrared interference figure method for reconstructing |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106908148B (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107990983B (en) * | 2017-11-21 | 2020-05-12 | 北京空间机电研究所 | Digital undersampling method for narrow-spectrum interference signal |
CN108593111B (en) * | 2018-05-31 | 2019-09-24 | 北京航空航天大学 | The moving imaging simulation method and device of Space-Time Modulation inteference imaging spectrometer |
CN109186767B (en) * | 2018-09-25 | 2021-06-01 | 天津大学 | Method for sampling interferogram in software-triggered Fourier transform spectrometer |
CN116593417B (en) * | 2023-03-30 | 2023-10-31 | 无锡迅杰光远科技有限公司 | Aplanatic sampling method, device and medium of Fourier transform infrared spectrometer |
CN116662789B (en) * | 2023-07-31 | 2023-10-27 | 杭州泽天春来科技有限公司 | Signal processing method and system of Fourier infrared spectrometer and readable storage medium |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1575524A (en) * | 2001-08-23 | 2005-02-02 | 华盛顿州大学 | Image acquisition with depth enhancement |
-
2017
- 2017-03-06 CN CN201710128061.XA patent/CN106908148B/en active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1575524A (en) * | 2001-08-23 | 2005-02-02 | 华盛顿州大学 | Image acquisition with depth enhancement |
Non-Patent Citations (4)
Title |
---|
"软"过零检测在红外光谱测量系统中的应用研究;任利兵等;《中国光学学会2010年光学大会论文集2010年》;20101231;1-8 * |
《两种红外干涉图采集及光谱复原方法的对比研究》;李妍等;《光学学报》;20150930;第35卷(第9期);正文全文 * |
光纤傅里叶光谱仪干涉图均匀抽样方法;李保生等;《红外与毫米波学报》;20070630;第26卷(第3期);1-8 * |
基于光栅信号控制干涉图数据采样方法研究;任兴;《科技资讯》;20161231;1-8 * |
Also Published As
Publication number | Publication date |
---|---|
CN106908148A (en) | 2017-06-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106908148B (en) | A kind of infrared interference figure method for reconstructing | |
Vallisneri | Synthetic LISA: Simulating time delay interferometry in a model LISA | |
Neff et al. | VLA Observations of the Nearby Merger NGC 4038/4039: H II Regions and Supernova Remnants in the “Antennae” | |
CN105030201B (en) | Frequency sweep OCT digit phases antidote and system based on static region information | |
CN101996418B (en) | Flame sampling device with temperature information and simulation method | |
Castellanos et al. | Algorithmic and methodological developments towards full waveform inversion in 3D elastic media | |
JPS63282665A (en) | Signal analyser | |
CN105716536A (en) | Three-dimensional digital speckle pattern interferometry synchronous measurement method and device | |
CN108563422A (en) | Randomizer and random number generation method | |
Felipe et al. | Inversions of synthetic umbral flashes: effects of scanning time on the inferred atmospheres | |
Yu et al. | Off-the-grid vertical seismic profile data regularization by a compressive sensing method | |
Van Der Jeught et al. | Optimized loss function in deep learning profilometry for improved prediction performance | |
CN110865385A (en) | Coherent superposition state source super-resolution quantum ranging system | |
Kulkarni et al. | Measurement of in-plane and out-of-plane displacements and strains using digital holographic moiré | |
Juusola et al. | Visual acuity for moving objects in first-and second-order neurons of the fly compound eye | |
CN103791845B (en) | The measuring method of optical lateral thin tail sheep and device based on laser high-order transverse mode | |
CN103823241A (en) | Method for calculating refraction residual static correcting values in offset domains | |
CN106840039A (en) | A kind of quick three-dimensional measuring method based on composite grating projection | |
CN205195039U (en) | Spectrum synthesis system based on multiplexing body holographic grating | |
CN108871595B (en) | Super time resolution shock velocity calculation method | |
CN102184555A (en) | Color clustering method for central color fringes of De Bruijn color structure light coding image | |
WO2016174660A1 (en) | Intensity waveform reconstruction from intensity autocorrelation | |
Zavareh et al. | A novel continuous time ternary encoding based SS-OCT calibration | |
CN113011107B (en) | One-dimensional optical fiber sensing signal phase recovery method based on deep convolutional neural network | |
Barajas et al. | Towards an on-chip signal processing solution for the online calibration of SS-OCT systems |
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 |