CN105136295A - Solution method and apparatus for spectrum non-uniformity in same image of AOTF - Google Patents

Solution method and apparatus for spectrum non-uniformity in same image of AOTF Download PDF

Info

Publication number
CN105136295A
CN105136295A CN201510593719.5A CN201510593719A CN105136295A CN 105136295 A CN105136295 A CN 105136295A CN 201510593719 A CN201510593719 A CN 201510593719A CN 105136295 A CN105136295 A CN 105136295A
Authority
CN
China
Prior art keywords
aotf
spectrum
diffraction
ccd
driving frequency
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
CN201510593719.5A
Other languages
Chinese (zh)
Other versions
CN105136295B (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.)
North University of China
Original Assignee
North University of China
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 North University of China filed Critical North University of China
Priority to CN201510593719.5A priority Critical patent/CN105136295B/en
Publication of CN105136295A publication Critical patent/CN105136295A/en
Application granted granted Critical
Publication of CN105136295B publication Critical patent/CN105136295B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Spectrometry And Color Measurement (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

The invention relates to the technical field of AOTF spectrum imaging, and more specifically relates to a solution method and an apparatus for spectrum non-uniformity in the same image of an AOTF. The solution method is mainly used for solving spectrum non-uniformity of the same image in the imaging of the AOTF and is a method for realizing spectrum correction via subsequent data processing, and the central wavelength of AOTF diffraction corresponding to a pixel xi in the x direction of any CCD with any driving frequency f is obtained via a formula shown in the description. The solution method and the apparatus are mainly applied to spectrum imaging of the AOTF, compared with the conventional AOTF imaging spectrum, the average error of the spectrum in the same image can reduce one order of magnitude, and the spectrum measuring precision of the AOTF is obviously improved.

Description

The solution that in a kind of AOTF same width figure, spectrum is uneven and device
Technical field
The present invention relates to AOTF spectral imaging technology field, more specifically, relate to solution and device that in a kind of AOTF same width figure, spectrum is uneven, the spectrum being mainly used to solve same width figure in AOTF imaging is uneven, is a kind of method being realized middle-low alloy steels by follow-up data process.
Background technology
Acousto-optic tunable filter (Acousto-optictunablefilter, AOTF) has that volume is little, tuned speed is fast, spectral range is wide, diffraction efficiency advantages of higher.Imaging spectrometer based on AOTF is applied in remote sensing, environmental monitoring, biomedicine and food inspection.
In AOTF, spectral measurement precision plays critical effect to light spectrum image-forming quality.In practical study, for AOTF imaging spectral, when light impinges perpendicularly in AOTF, a certain fixing wavelength that its a certain driving frequency is corresponding.But the light of space diverse location target can not all vertical incidence, this just causes the incident angle entering AOTF light different, thus cause its diffraction centre wavelength of the light of different incidence angles also to change, make CCD diverse location detect the light of different diffraction wavelength.And in existing survey AOTF imaging spectral, the detection of CCD diverse location obtains: during vertical incidence, the light of wavelength corresponding to fixed drive frequency.Certain error is there is in this with regard to making the measurement of theoretical spectral measurement and actual spectrum.Pursuing today of high precision spectral measurement, thisly can not to be ignored by theoretical and the actual error brought.So propose the method that can improve AOTF spectral measurement precision by theoretical and experimental study, the method can reduce measuring error to greatest extent, for later high precision spectral measurement is laid a good foundation.
Summary of the invention
In order to overcome the problem that in prior art, in AOTF imaging, same width figure spectrum is uneven, the invention provides and a kind ofly adopt quadratic fit to revise the method for spectrum, thus improve the measuring accuracy of AOTF spectrum.
In order to solve the problems of the technologies described above, the technical solution used in the present invention is:
The solution that in a kind of AOTF same width figure, spectrum is uneven, in the spectral range and angular range of AOTF, under any driving frequency f, diffraction central wavelength lambda (f in incident angle α and this situation, α) can parabolic relation be approximately, be specially:
λ(f,α)=Aα 20(f)(1)
In above formula, AOTF diffraction light centre wavelength when λ (f, α) is driving frequency is f, incident angle is α; λ 0f () is driving frequency diffraction light centre wavelength when being f, vertical incidence AOTF; Described λ 0f (), for carrying when AOTF dispatches from the factory, also can adopt polychromatic light vertical incidence AOTF and be recorded by high resolution spectrometer; Described A adopts following steps to obtain:
AOTF is placed on high-accuracy angle measurement unit and makes the n of AOTF in a spectral range driving frequency f respectively 1, f 2, f 3... f nlower work, for each driving frequency f nunder, parallel polychromatic light by AOTF diffraction, and is obtained the diffraction wavelength λ (f under different incidence angles α by high resolution spectrometer n, α), to α and λ (f n, α) and carry out least square fitting, and require that quafric curve must be crossed [α=0, λ 0(f n)] point, the A obtained after matching n, A nfor A in (1) formula, can approximate representation be:
A n=Bλ 0(f n)+C(2)
B and C in above-mentioned (2) formula can by the n obtained above (λ 0(f 1), A 1), (λ 0(f 2), A 2), (λ 0(f 3), A 3) ... (λ 0(f n), A n) adopt least square fitting to obtain, in conjunction with above-mentioned (1) formula and (2) Shi Ke get:
λ(f,α)=[Bλ 0(f)+C]α 20(f)(3)
Because CCD is placed in the focus of imaging len, so the picture dot x of CCD iwith the pass of the focal length F of incident angle α and imaging len be:
α = a r c t a n ( i l F ) - - - ( 4 )
In above formula, l is the length of side of CCD picture dot; I Wei Yi CCD center is 0 at diffraction direction picture dot x isequence number;
In λ (f, α), λ is the function about f and α, λ (f, x i) in λ be relevant f and picture dot x ithe function of sequence number i, therefore can obtain the picture dot x of any driving frequency f, arbitrarily CCD by above-mentioned (4) formula (3) formula of bringing into ithe centre wavelength of corresponding A OTF diffraction:
λ ( f , x i ) = [ Bλ 0 ( f ) + C ] arctan 2 ( i l F ) + λ 0 ( f ) - - - ( 5 ) .
Described high-accuracy angle measurement unit is spectroscope, described spectroscope angle-measurement accuracy 1'=0.0167 °.
The solving device that in a kind of AOTF same width figure, spectrum is uneven, comprise preposition optical system, the first polaroid, diaphragm, AOTF, the second polaroid, imaging len, CCD and computing machine, described preposition optical system is made up of the first lens, the second lens and the 3rd lens, and described first lens, the second lens, the 3rd lens, the first polaroid, diaphragm, AOTF, the second polaroid, imaging len, CCD and computing machine are arranged in order.
The polarization direction of described first polaroid and the second polaroid is orthogonal.
The beneficial effect that compared with prior art the present invention has is:
The present invention adopts quadratic fit to revise the problem that in AOTF imaging, same width figure spectrum is uneven, and compared with existing AOTF imaging spectral, the average error of same width figure spectrum can reduce an order of magnitude, significantly improves the spectral measurement precision of AOTF.
Accompanying drawing explanation
Below by accompanying drawing, the specific embodiment of the present invention is described in further detail.
Fig. 1 is the basic index path of apparatus of the present invention;
Fig. 2 is AOTF angle of diffraction spectral distribution graph of the present invention;
Fig. 3 is CCD polishing wax distribution schematic diagram;
Fig. 4 is driving frequency when being 116MHz vertical incidence, reality and revised numerical value schematic diagram;
Fig. 5 is driving frequency when being 120MHz vertical incidence, reality and revised numerical value schematic diagram;
Fig. 6 is driving frequency when being 136MHz vertical incidence, reality and revised numerical value schematic diagram;
Fig. 7 is driving frequency when being 140MHz vertical incidence, reality and revised numerical value schematic diagram.
In figure: 1 be diaphragm for preposition optical system, 2 is the first polaroid, 3,4 be imaging len for AOTF, 5 is the second polaroid, 6,7 be CCD, 8 be the 3rd lens for computing machine, 9 be the first lens, 10 is the second lens, 11.
Embodiment
The invention will be further described by reference to the accompanying drawings for embodiment below.
As shown in Figure 1, the solving device that in a kind of AOTF same width figure, spectrum is uneven, comprise preposition optical system 1, first polaroid 2, diaphragm 3, AOTF4, the second polaroid 5, imaging len 6, CCD7 and computing machine 8, preposition optical system 1 is made up of the first lens 9, second lens 10 and the 3rd lens 11, measured target light is successively after preposition optical system 1, first polaroid 2 and diaphragm 3 at a distance, enter AOTF4 diffraction to filter, by the second polaroid 5, be imaged on CCD7 through imaging len 6, finally by computing machine 8 by middle-low alloy steels.
Preposition optical system 1 is compressed to by distant place measured target angle of light within the scope of the field angle of AOTF4 permission, and a kind of parallel beam is become another kind of parallel beam, sees Fig. 1.Because measured target distance AOTF imaging spectrum system is distant, add the restriction of diaphragm 3, diverse location target can be similar to and see as: under different field angle, a branch of very thin parallel beam.The namely measured target of distant place diverse location, the incident angle entering AOTF imaging spectrum system is different, and then the incident angle entering AOTF4 after preposition optical system 1 is different, is imaged on the different picture dot of CCD7, as shown in Figure 1 through imaging len 6; Again due to when fixed drive frequency, AOTF4 diffraction wavelength is relevant with the angle of incident AOTF4, and wavelength just can be revised in therefore different with CCD7 according to the focal length of imaging len 6 picture dot positions, as shown in Figure 2.
The polarization direction of described first polaroid 2 and the second polaroid 5 is completely orthogonal, mainly in order to eliminate in AOTF4 0 grade and-1 order diffraction light to the impact of AOTF light spectrum image-forming.
As shown in Figure 2, suppose that AOTF diffraction direction is x direction, because AOTF4 only has the centre wavelength of diffraction direction relevant with incident angle, the practical center wavelength of the different picture dots therefore on CCD7 as shown in Figure 3.
Due to impacts such as crystal pro cessing precision, installation accuracies, there is some difference for the AOTF that actual production goes out and theoretical meeting, therefore for the AOTF4 of reality, under each driving frequency, the relation entering AOTF4 diffraction centre wavelength and incident angle can adopt the spectrometer of High-precision angle measurement mechanism and high spectral resolution to record, due under will recording each driving frequency, diffraction centre wavelength corresponding to different incidence angles, quantities is very large, therefore adopts the mode of matching.Specific as follows: actual AOTF, due in certain spectral range and angular range, can parabolic relation be approximately for diffraction wavelength λ (f, α) in incident angle α under any driving frequency f and this situation, and about incident angle α=0 ° both sides symmetry:
λ(f,α)=Aα 20(f)(1)
In above formula, AOTF4 diffraction light centre wavelength when λ (f, α) is driving frequency is f, incident angle is α; λ 0f () is driving frequency diffraction light centre wavelength when being f, vertical incidence AOTF; λ 0f () is also the wavelength that existing AOTF imaging spectral generally adopts, do not consider the impact of incident angle α.
Wherein, in (1) formula, λ 0f () producer when buying AOTF can provide, or also can adopt polychromatic light vertical incidence AOTF and be recorded by high resolution spectrometer.
Wherein, (1) in formula, A can adopt high resolution spectrometer and High-precision angle measurement mechanism measurement data least square fitting to obtain, concrete grammar is: be placed on high-accuracy angle measurement unit AOTF4 (as spectroscope, angle-measurement accuracy 1'=0.0167 °), make the n of AOTF in a spectral range driving frequency f respectively 1, f 2, f 3... f nlower work, for each driving frequency f nunder, parallel polychromatic light by AOTF diffraction, and is obtained the diffraction wavelength λ (f under different incidence angles α by high resolution spectrometer n, α), to α and λ (f n, α) and carry out least square fitting, and require that quafric curve must be crossed [α=0, λ 0(f n)] point, obtain the A in (1) formula after matching n, see the following form 1:
Table 1 different diffraction light wavelength lambda 0a value corresponding to (f)
Diffraction wavelength λ 0(f)/nm 617.41 600.88 585.62 571.10 557.59 545.10 533.37
A 0.1890 0.1774 0.1634 0.1522 0.1481 0.1319 0.1198
Wherein, due to the A under different driving frequency ndifference, but can approximate representation be:
A n=Bλ 0(f n)+C(2)
(2) B and C in formula can by the n obtained above (λ 0(f 1), A 1), (λ 0(f 2), A 2), (λ 0(f 3), A 3) ... (λ 0(f n), A n) adopt least square fitting to obtain, after matching, (2) are: A=0.000802 λ 0(f)-0.3048.
Wherein, in conjunction with (1) formula and (2) Shi Ke get:
λ(f,α)=[0.000802λ 0(f)-0.3048]α 20(f)(3)
Wherein, vertical incidence value (i.e. current AOTF generally adopt numerical value), actual measured value and modified value driving frequency be the result of 116MHz, 120MHz, 136MHz and 140MHz respectively as shown in Figure 4, Figure 5, Figure 6 and Figure 7.
Wherein, after matching correction with not have to revise after the relative error of (i.e. current AOTF generally adopt value) contrast as shown in table 2, revised relative error compare revise before to reduce an order of magnitude.
Relative error magnitudes before and after revising under table 2 different diffraction optical wavelength
Wherein, composition graphs 2 and Fig. 3 known, the picture dot x in the x direction of CCD7 iwith the pass of the focal length F of incident angle α and imaging len 6 be:
α = a r c t a n ( i l F ) - - - ( 4 )
L is the length of side of CCD picture dot, and i Wei Yi CCD center is 0 at diffraction direction (i.e. x direction) picture dot x isequence number.
By (4) formula, (3) formula of bringing into just can obtain the picture dot x in any driving frequency f, arbitrarily the x direction of CCD ithe centre wavelength of corresponding A OTF diffraction:
λ ( f , x i ) = [ 0.000802 λ 0 ( f ) - 0.3048 ] arctan 2 ( i l F ) + λ 0 ( f ) - - - ( 5 ) .

Claims (4)

1. the solution that in AOTF same width figure, spectrum is uneven, is characterized in that:
In the spectral range and angular range of AOTF, under any driving frequency f, in incident angle α and this situation, diffraction central wavelength lambda (f, α) can be approximately parabolic relation, is specially:
λ(f,α)=Aα 20(f)(1)
In above formula, AOTF diffraction light centre wavelength when λ (f, α) is driving frequency is f, incident angle is α; λ 0f () is driving frequency diffraction light centre wavelength when being f, vertical incidence AOTF; Described λ 0f (), for carrying when AOTF dispatches from the factory, also can adopt polychromatic light vertical incidence AOTF and be recorded by high resolution spectrometer; Described A adopts following steps to obtain:
AOTF is placed on high-accuracy angle measurement unit and makes the n of AOTF in a spectral range driving frequency f respectively 1, f 2, f 3... f nlower work, for each driving frequency f nunder, parallel polychromatic light by AOTF diffraction, and is obtained the diffraction wavelength λ (f under different incidence angles α by high resolution spectrometer n, α), to α and λ (f n, α) and carry out least square fitting, and require that quafric curve must be crossed [α=0, λ 0(f n)] point, the A obtained after matching n, A nfor A in (1) formula, can approximate representation be:
A n=Bλ 0(f n)+C(2)
B and C in above-mentioned (2) formula can by the n obtained above (λ 0(f 1), A 1), (λ 0(f 2), A 2), (λ 0(f 3), A 3) ... (λ 0(f n), A n) adopt least square fitting to obtain, in conjunction with above-mentioned (1) formula and (2) Shi Ke get:
λ(f,α)=[Bλ 0(f)+C]α 20(f)(3)
Because CCD is placed in the focus of imaging len, so the picture dot x of CCD iwith the pass of the focal length F of incident angle α and imaging len be:
α = a r c t a n ( i l F ) - - - ( 4 )
In above formula, l is the length of side of CCD picture dot; I Wei Yi CCD center is 0 at diffraction direction picture dot x isequence number;
In λ (f, α), λ is the function about f and α, λ (f, x i) in λ be relevant f and picture dot x ithe function of sequence number i, therefore can obtain the picture dot x of any driving frequency f, arbitrarily CCD by above-mentioned (4) formula (3) formula of bringing into ithe centre wavelength of corresponding A OTF diffraction:
λ ( f , x i ) = [ Bλ 0 ( f ) + C ] arctan 2 ( i l F ) + λ 0 ( f ) - - - ( 5 ) .
2. the solution that in a kind of AOTF according to claim 1 same width figure, spectrum is uneven, is characterized in that: described high-accuracy angle measurement unit is spectroscope, described spectroscope angle-measurement accuracy 1'=0.0167 °.
3. the solving device that in an AOTF same width figure, spectrum is uneven, it is characterized in that: comprise preposition optical system (1), first polaroid (2), diaphragm (3), AOTF (4), second polaroid (5), imaging len (6), CCD (7) and computing machine (8), described preposition optical system (1) is by the first lens (9), second lens (10) and the 3rd lens (11) composition, described first lens (9), second lens (10), 3rd lens (11), first polaroid (2), diaphragm (3), AOTF (4), second polaroid (5), imaging len (6), CCD (7) and computing machine (8) are arranged in order.
4. the solving device that in a kind of AOTF according to claim 3 same width figure, spectrum is uneven, is characterized in that: the polarization direction of described first polaroid (2) and the second polaroid (5) is orthogonal.
CN201510593719.5A 2015-09-17 2015-09-17 The uneven solution of spectrum and device in a kind of same width figure of AOTF Expired - Fee Related CN105136295B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510593719.5A CN105136295B (en) 2015-09-17 2015-09-17 The uneven solution of spectrum and device in a kind of same width figure of AOTF

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510593719.5A CN105136295B (en) 2015-09-17 2015-09-17 The uneven solution of spectrum and device in a kind of same width figure of AOTF

Publications (2)

Publication Number Publication Date
CN105136295A true CN105136295A (en) 2015-12-09
CN105136295B CN105136295B (en) 2017-03-29

Family

ID=54721730

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510593719.5A Expired - Fee Related CN105136295B (en) 2015-09-17 2015-09-17 The uneven solution of spectrum and device in a kind of same width figure of AOTF

Country Status (1)

Country Link
CN (1) CN105136295B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105717076A (en) * 2016-05-04 2016-06-29 中北大学 Spectrum SPR imaging sensing system based on acousto-optic light filtration
CN109115339A (en) * 2018-05-16 2019-01-01 中北大学 One kind being based on AOTF and intensity modulated high speed high spectrum full polarization imaging device and method

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5828451A (en) * 1997-09-30 1998-10-27 Northrop Grumman Corporation Spectral imaging system and method employing an acousto-optic tunable filter for wavelength selection with increased field of view brightness
CN201041514Y (en) * 2007-06-14 2008-03-26 武汉大学 Portable imaging spectrum instrument
CN101561388A (en) * 2008-04-18 2009-10-21 中国科学院上海技术物理研究所 Acousto-optic tunable filter imaging spectrometer
CN201497574U (en) * 2009-04-24 2010-06-02 中国科学院西安光学精密机械研究所 Programmable polarization hyperspectral imager
CN102645278A (en) * 2012-05-07 2012-08-22 北京航空航天大学 Radio frequency drive control method for hyper-spectral imager of acousto-optic tunable filter

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5828451A (en) * 1997-09-30 1998-10-27 Northrop Grumman Corporation Spectral imaging system and method employing an acousto-optic tunable filter for wavelength selection with increased field of view brightness
CN201041514Y (en) * 2007-06-14 2008-03-26 武汉大学 Portable imaging spectrum instrument
CN101561388A (en) * 2008-04-18 2009-10-21 中国科学院上海技术物理研究所 Acousto-optic tunable filter imaging spectrometer
CN201497574U (en) * 2009-04-24 2010-06-02 中国科学院西安光学精密机械研究所 Programmable polarization hyperspectral imager
CN102645278A (en) * 2012-05-07 2012-08-22 北京航空航天大学 Radio frequency drive control method for hyper-spectral imager of acousto-optic tunable filter

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105717076A (en) * 2016-05-04 2016-06-29 中北大学 Spectrum SPR imaging sensing system based on acousto-optic light filtration
CN109115339A (en) * 2018-05-16 2019-01-01 中北大学 One kind being based on AOTF and intensity modulated high speed high spectrum full polarization imaging device and method
CN109115339B (en) * 2018-05-16 2021-04-20 中北大学 High-speed hyperspectral full-polarization imaging device and method based on AOTF and intensity modulation

Also Published As

Publication number Publication date
CN105136295B (en) 2017-03-29

Similar Documents

Publication Publication Date Title
CN103363901B (en) A kind of scaling method towards coaxial alignment microassembly system
US9574967B2 (en) Wavefront measurement method, shape measurement method, optical element manufacturing method, optical apparatus manufacturing method, program, and wavefront measurement apparatus
CN101334267A (en) Digital image feeler vector coordinate transform calibration and error correction method and its device
CN107462166B (en) Long stroke, high precision displacement measurement method based on diffraction grating
CN108510546B (en) Camera calibration method suitable for map and structure information synchronous detection system
CN109780992A (en) Interferometer measuration system error calibrating method based on the processing of optical flat stripe pattern
CN104457600A (en) Testing device of optical fiber collimator array
CN104567752A (en) Stray light elimination dual-optical-path optical centering instrument
CN105066910A (en) Electro-optic crystal Z axis deviation angle measurement device and measurement method
CN106932173A (en) High-precision heavy-caliber optical grating five degree of freedom splices the measuring method of precision
CN109579779A (en) High-precision high-frequency rings anti-interference big working distance autocollimation and method
CN105352915A (en) Refractive index two-dimensional distribution dynamic measurement method
Huo et al. A flexible calibration method based on small planar target for defocused cameras
CN104931138A (en) Method of using prism to increase AOTF spectrum imaging quality and apparatus thereof
CN105136295A (en) Solution method and apparatus for spectrum non-uniformity in same image of AOTF
CN105450912A (en) Scanning method-based real-time field butting device and method for surface array CCD detectors
CN109579776A (en) Anti-interference big working distance autocollimation and method in high precision
CN111220088A (en) Measurement system and method
CN103438803B (en) Computer vision technique accurately measures the method for Rectangular Parts size across visual field
CN206146834U (en) V V -prism refractometer based on auto -collimation and CCD vision technique
US11536560B2 (en) Optical apparatus, optical system, and method for measuring an amount of strain of an object
CN106247992B (en) A kind of high-precision, wide scope and big working distance autocollimation and method
CN108362222A (en) Based on the multidirectional novel point-diffraction interference measuring system of nonzero digit for tilting carrier frequency
CN108955880B (en) Coaxiality calibration method for high-precision ultraviolet double-grating spectrometer
CN102269582A (en) Spatial three-dimensional angle measurement apparatus

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170329

Termination date: 20170917

CF01 Termination of patent right due to non-payment of annual fee