CN107741280A - Light beam forward direction turbulence transfer drift angle and reverse transfer angle of arrival coefficient correlation measuring method - Google Patents
Light beam forward direction turbulence transfer drift angle and reverse transfer angle of arrival coefficient correlation measuring method Download PDFInfo
- Publication number
- CN107741280A CN107741280A CN201711174454.0A CN201711174454A CN107741280A CN 107741280 A CN107741280 A CN 107741280A CN 201711174454 A CN201711174454 A CN 201711174454A CN 107741280 A CN107741280 A CN 107741280A
- Authority
- CN
- China
- Prior art keywords
- laser
- array
- light spot
- spot image
- ccd
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000012546 transfer Methods 0.000 title claims abstract description 27
- 230000003287 optical effect Effects 0.000 claims abstract description 49
- 238000004088 simulation Methods 0.000 claims abstract description 49
- 238000005259 measurement Methods 0.000 claims abstract description 12
- 239000011159 matrix material Substances 0.000 claims description 38
- 239000011521 glass Substances 0.000 claims description 27
- 239000002184 metal Substances 0.000 claims description 15
- 238000001514 detection method Methods 0.000 claims description 3
- 239000000523 sample Substances 0.000 claims description 3
- 230000005855 radiation Effects 0.000 claims description 2
- 238000009738 saturating Methods 0.000 claims description 2
- 238000007619 statistical method Methods 0.000 abstract 1
- 230000005540 biological transmission Effects 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 6
- 230000003044 adaptive effect Effects 0.000 description 5
- 238000011160 research Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000004020 conductor Substances 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000012360 testing method Methods 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
- G01J9/00—Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/26—Measuring arrangements characterised by the use of optical techniques for measuring angles or tapers; for testing the alignment of axes
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
- H04B10/50—Transmitters
- H04B10/501—Structural aspects
- H04B10/503—Laser transmitters
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Optics & Photonics (AREA)
- Electromagnetism (AREA)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Length Measuring Devices By Optical Means (AREA)
- Image Analysis (AREA)
Abstract
The present invention discloses the coefficient correlation measuring method of a kind of light beam forward direction turbulence transfer drift angle and reverse transfer angle of arrival, it is characterised in that：Using convection type atmospheric turbulence simulation device, simulated atmosphere turbulent flow, beam Propagation pass through the atmospheric turbulance of lab simulation indoors, produce beam drift and angle of arrival change.Measure the light angle of arrival along two beam laser of identical optical axis reverse transfer and beam drift angle respectively first, and measurement result is saved in computer documents；Then, according to the measurement result being stored in computer documents, the coefficient correlation of light beam forward direction turbulence transfer drift angle and reverse transfer angle of arrival is drawn using statistical analysis technique.
Description
Technical field
The present invention relates to a kind of light beam forward direction turbulence transfer drift angle and reverse transfer angle of arrival coefficient correlation measuring method,
Belong to laser atmospheric turbulence transmission technique field.
Background technology
Beam drift phenomenon can occur when being transmitted in atmospheric turbulance for laser beam, i.e., the beam cross-section light in receiving plane
There is random movement in spot centroid position.Assuming that sometime, the beam cross-section facula mass center position in receiving plane is deviateed
The distance of beam cross-section facula mass center position during without Turbulent Flow Effects is rc, beam emissions plane to beam reception plane away from
From for L, then the beam drift angle at the moment is β=arctan (rc/L).R under normal circumstancescMuch smaller than L, therefore β ≈ rc/L.Swash
Light beam is received after atmospheric turbulance transmits by a finite aperture receiver, and focuses on a ccd detector, CCD detection
With the time randomized jitter can occur for the barycenter of the light spot image of device output, and which reflects the change of the aperture averaging angle of arrival of light wave.
Theoretical research shows that beam drift is mainly the contribution for being derived from the turbulent eddy close to transmitting terminal, and light arrival angle fluctuation is mainly
Contribution from the turbulent eddy close to receiving terminal., can be in an optical transmitter and receiver in two-way laser atmospheric turbulance Transmission system
It is middle to measure the light angle of arrival received, and carry out deflection angle using transmitting light beam of the quick titling mirror to the optical transmitter and receiver accordingly and mend in advance
Repay, so as to reach the purpose of compensation beam drift.This is actually a kind of adaptive light beam lift-off technology.This adaptive light beam
The correlation that the performance of lift-off technology is depended between the light angle of arrival received and the drift angle for launching light beam.If launch light
The drift angle of beam and the light arrival angle dependence received are very low, then it is pre- to carry out deflection angle to transmitting light beam using quick titling mirror
The effect of compensation can be very poor.Therefore, the transmitting beam drift angle of an optical transmitter and receiver and the coefficient correlation for receiving light angle of arrival are measured
It is significant to the adaptive light beam lift-off technology of research and analysis.As shown in figure 1, optical transmitter and receiver A is a branch of to optical transmitter and receiver B transmittings
Laser, optical transmitter and receiver B launch beam of laser to optical transmitter and receiver A, are referred to as positive transmission from optical transmitter and receiver A to optical transmitter and receiver B transmission, from light end
Machine B to optical transmitter and receiver A transmission are referred to as reverse transfer, beam drift angle of the laser that optical transmitter and receiver A launches after atmospheric turbulance transmits
Referred to as light beam forward direction turbulence transfer drift angle, arrival of the laser that optical transmitter and receiver B launches after atmospheric turbulance transmits at optical transmitter and receiver A
Angle is referred to as the reverse turbulence transfer angle of arrival of light beam.It is rapid for the ease of control air in Study of Laser atmospheric turbulance transmission characteristic
Stream condition, atmospheric turbulence simulation device simulated atmosphere turbulent flow indoors can be utilized.The present invention provides a kind of light beam forward direction turbulent flow
Drift angle and the coefficient correlation measuring method of reverse transfer angle of arrival are transmitted, measuring beam passes through lab simulation atmospheric turbulance is positive
The statistical correlation coefficient of the angle of arrival of drift angle and light beam after lab simulation atmospheric turbulance reverse transfer after defeated.
The content of the invention
It is an object of the invention to provide a kind of light beam forward direction turbulence transfer drift angle and reverse transfer angle of arrival phase relation
Number measuring method, experiment number is provided for adaptive light beam lift-off technology of the research and analysis based on transmitting beam deflection angle precompensation
According to support.
The technical proposal of the invention is realized in this way：Light beam forward direction turbulence transfer drift angle and reverse transfer angle of arrival phase
Close coefficient measuring method, it is characterised in that：The laser A001 that first laser device is sent enters after the first lens and the second lens
It is mapped on the first optical splitter, from the laser A001 of the first optical splitter transmissive from the first of convection type atmospheric turbulence simulation device
Glass window, which enters in the atmospheric turbulance of simulation, to be transmitted, and is transmitted from the second glass window of convection type atmospheric turbulence simulation device
Go out, the laser A001 transmitted away from the second glass window is incided on the second optical splitter, by swashing for the second spectrophotometric reflection
Light A001 is incided in rectangle diffuse reflecting screen；The hot spot in rectangle diffuse reflecting screen is continuously shot with digital camera, digital camera connects
The continuous obtained light spot image sequence that shoots is transmitted to second computer system；The laser A002 that second laser is sent is by the
Incided after three lens and the 4th lens on the second optical splitter, it is big from convection type from the laser A002 of the second optical splitter transmissive
Second glass window of gas turbulent flow simulation device, which enters in the atmospheric turbulance of simulation, to be transmitted, and is simulated and filled from convection type atmospheric turbulance
The first glass window put is transmitted away, and the laser A002 transmitted away from the first glass window is incided on the first optical splitter,
Incided by the laser A002 of the first spectrophotometric reflection on the 5th lens, the laser A002 after the 5th lens focus is incided
On ccd detector, the light spot image that ccd detector detects is transmitted to first computer system；Ccd detector is saturating to the 5th
The distance of mirror is the focal length of the 5th lens；Convection type atmospheric turbulence simulation device for closing square metal bucket, normal work
Atmospheric turbulance is produced when making in metal bucket.
1) determine in the case of the atmospheric turbulance in the absence of simulation, the barycenter for the hot spot that ccd detector detects, specifically
Method is as follows：
Step Step101：Convection type atmospheric turbulence simulation device is set not work, in convection type atmospheric turbulence simulation device
Air in metal bucket is tranquil air；Make second laser normal work and launch laser A002；Make ccd detector and first
Computer system normal work；
Step Step102：The frame light spot image B001 that ccd detector is detected is sent to first computer system
In, and light spot image B001 pixel grey scale matrix is stored in the two-dimensional array I in the internal memory of first computer systemccdIn,
Array IccdLine number equal to light spot image B001 pixel grey scale matrix line number, array IccdColumns be equal to light spot image
B001 pixel grey scale matrix column number；
Step Step103：CalculateWherein Iccd(i1,j1)
Represent array IccdI-th1Row, jth1The value of column element, Nr1Represent array IccdLine number, Nc1Represent array IccdColumns；
irz1And jcz1Line position corresponding to being the barycenter for the hot spot that ccd detector detects is put and column position.
2) adjust and demarcate digital camera, specific method is as follows：
Step Step201：Convection type atmospheric turbulence simulation device is set not work, in convection type atmospheric turbulence simulation device
Air in metal bucket is tranquil air；Make first laser device normal work and launch laser A001；Make digital camera and second
Computer system normal work；The position of rectangle diffuse reflecting screen is adjusted, makes to incide laser A001's in rectangle diffuse reflecting screen
Hot spot is located at the core of rectangle diffuse reflecting screen；
Step Step202：First laser device is set not work, with the whole rectangle diffuse reflecting screen of white light；
Step Step203：The mark point of four black is added in rectangle diffuse reflecting screen, first mark point is P1, the
Two mark points are P2, the 3rd mark point is P3, the 4th mark point is P4, mark point P1、P2、P3、P4It is unrestrained anti-positioned at rectangle
Penetrate the core of screen and spaced apart, P1And P2Line perpendicular to P3And P4Line, P1And P2Line parallel to
One rectangular edges EV1, P of rectangle diffuse reflecting screen3And P4Line parallel to another rectangular edges EH1 of rectangle diffuse reflecting screen, square
Shape side EV1 is perpendicular to rectangular edges EH1；Digital camera is adjusted, makes connection P1And P2Line segment in the image that digital camera photographs
Same one-row pixels are extended only through, make connection P3And P4Line segment same row pixel is extended only through in the image that digital camera photographs；
Step Step204：Calculate connection P1And P2The number of pixels that is passed through in the image that digital camera photographs of line segment
NP12, calculate connection P3And P4The number of pixels N that is passed through in the image that digital camera photographs of line segmentP34；It is unrestrained anti-to measure rectangle
The connection P penetrated on screen1And P2Line segment length L12, measure the connection P in rectangle diffuse reflecting screen3And P4Line segment length L34；
Step Step205：Mark point P is removed from rectangle diffuse reflecting screen1、P2、P3、P4；
3) determine in the case of the atmospheric turbulance in the absence of simulation, in the rectangle diffuse reflecting screen that digital camera photographs
The barycenter of hot spot, specific method are as follows：
Step Step301：Convection type atmospheric turbulence simulation device is set not work, in convection type atmospheric turbulence simulation device
Air in metal bucket is tranquil air；Make first laser device normal work and launch laser A001；Make digital camera and second
Computer system normal work；
Step Step302：The frame light spot image B002 that digital camera is photographed is sent in second computer system,
And light spot image B002 pixel grey scale matrix is stored in the two-dimensional array I in the internal memory of second computer systemcamIn, number
Group IcamLine number equal to light spot image B002 pixel grey scale matrix line number, array IcamColumns be equal to light spot image B002
Pixel grey scale matrix column number；
Step Step303：CalculateWherein Icam
(i2,j2) represent array IcamI-th2Row, jth2The value of column element, Nr2Represent array IcamLine number, Nc2Represent array Icam's
Columns；irz2And jcz2Line position corresponding to being the barycenter of the hot spot for the light spot image B002 that digital camera photographs is put and column position.
4) while measurement is incided the angle of arrival of the laser A002 on ccd detector and incided in rectangle diffuse reflecting screen
Laser A001 beam drift angle, specific method are as follows：
Step Step401：Make convection type atmospheric turbulence simulation device normal work, make first laser device normal work concurrent
Penetrate laser A001；Make second laser normal work and launch laser A002；Make ccd detector and first computer system normal
Work；Make digital camera and second computer system worked well；A new disk file is created in first computer system
FILE1, a new disk file FILE2 is created in second computer system；
Step Step402：The every frame light spot image B001 detected to ccd detector, perform following operation：
Step Step402-1：Light spot image B001 is sent in first computer system, and light spot image B001's
Pixel grey scale matrix is stored in the two-dimensional array I in the internal memory of first computer systemccdIn, array IccdLine number be equal to hot spot
The line number of image B001 pixel grey scale matrix, array IccdColumns be equal to light spot image B001 pixel grey scale matrix column
Number；
Step Step402-2：CalculateWherein Iccd(i1,
j1) represent array IccdI-th1Row, jth1The value of column element, Nr1Represent array IccdLine number, Nc1Represent array IccdRow
Number；
Step Step402-3：Orderβ1=arctan
(dw/f5) value be added to disk file FILE1 end, β1Represent light angle of arrival, f5The focal length of the 5th lens is represented,
Δh1For the pixels tall of ccd detector, Δw1For the pixel wide of ccd detector；For light spot image B001 operation knot
Beam；
Step Step403：The every frame light spot image B002 photographed to digital camera, perform following operation：
Step Step403-1：Light spot image B002 is sent in second computer system, and light spot image B002's
Pixel grey scale matrix is stored in the two-dimensional array I in the internal memory of second computer systemcamIn, array IcamLine number be equal to hot spot
The line number of image B002 pixel grey scale matrix, array IcamColumns be equal to light spot image B002 pixel grey scale matrix column
Number；
Step Step403-2：CalculateWherein Icam
(i2,j2) represent array IcamI-th2Row, jth2The value of column element, Nr2Represent array IcamLine number, Nc2Represent array Icam's
Columns；
Step Step403-3：OrderΔh2=L34/
NP34, Δw2=L12/NP12, β2=arctan (Dw/ F) value be added to disk file FILE2 end, β2Represent light beam drift
Move angle, F=da+db+dc, daRepresent from the first glass window to the optical axis distance of the second glass window, dbRepresent from the second glass
Window is to the optical axis distance of the second optical splitter, dcRepresent from the second optical splitter to the optical axis distance of rectangle diffuse reflecting screen；For light
Spot image B002 operation terminates.
5) after the measurement operation execution of the 4) part terminates, file FILE2 is copied in first computer system, meter
Count the angle of arrival for the laser A002 being mapped on ccd detector and incide the light beam drift of the laser A001 in rectangle diffuse reflecting screen
The coefficient correlation at angle is moved, specific method is as follows：
Step Step501：All smooth angle of arrival β being stored in file FILE11Value and be stored in file FILE2
In all beam drift angle betas2Value be loaded into the internal memory of first computer system；In the internal memory of first computer system
Create an one-dimension array AaoaTo store all smooth angle of arrival β loaded from file FILE11Value, create a dimension
Group AbwTo store all beam drift angle betas loaded from file FILE22Value；M-th be stored in file FILE1
β1Value be stored in AaoaIn [m], Aaoa[m] represents array AaoaM-th of element, m=1,2,3 ..., Num, Num represent preserve
β in file FILE11Number；M-th of β being stored in file FILE22Value be stored in AbwIn [m], Abw[m] table
Registration group AbwM-th of element, m=1,2,3 ..., Num, Num represents to be stored in β in file FILE22Number；
Step Step502：Incide the angle of arrival of the laser A002 on ccd detector and incide in rectangle diffuse reflecting screen
The laser A001 coefficient correlation at beam drift angle be：
Wherein
The measurement operation of the 4) part continuously performs TmTerminate after minute, i.e. ccd detector continuous probe TmMinute, numeral
Camera is continuously shot TmMinute；It is utilized respectively first computer system and second computer system, step Step402 and step
Step403 starts simultaneously at execution, and ccd detector detection two continuous frames light spot image B001 time interval is clapped equal to digital camera
Take the photograph two continuous frames light spot image B002 time interval；The β being stored in file FILE11Number be equal to be stored in file
β in FILE22Number.
The positive effect of the present invention is can be arrived by measuring acquisition light beam forward direction turbulence transfer drift angle with reverse transfer
Up to the coefficient correlation at angle, provided for adaptive light beam lift-off technology of the research and analysis based on transmitting beam deflection angle precompensation real
Test data.By the turbulent parameters that can be simulated with manual control convection type atmospheric turbulence simulation device, therefore can be easily
Obtain the light beam forward direction turbulence transfer drift angle under the conditions of different turbulent parameters and the coefficient correlation of reverse transfer angle of arrival measures
Value.
Brief description of the drawings
Fig. 1 is positive transmission and reverse transfer schematic diagram.
Fig. 2 is measuring system structural representation.
Fig. 3 is four density bullet point schematic diagrames in rectangle diffuse reflecting screen.
Fig. 4 is the picture element matrix of ccd detector and the width and high-level schematic of pixel.
Embodiment
In order that the feature and advantage of this method are more clearly understood, this method is made into one with reference to specific embodiment
The description of step.In the present embodiment, first laser device 101 and second laser 112 select continuous conductor laser, Tm=10
Minute, the time interval that ccd detector 106 detects two continuous frames light spot image B001 is 1 millisecond, and the shooting of digital camera 113 connects
Continuous two frame light spot image B002 time interval is 1 millisecond, Nr1=800, Nc1=800, Nr2=1000, Nc2=1000.
As shown in Fig. 2 the laser A001 that first laser device 101 is sent enters after the first lens 102 and the second lens 103
It is mapped on the first optical splitter 104, from the laser A001 of the transmissive of the first optical splitter 104 from convection type atmospheric turbulence simulation device
108 the first glass window 116, which enters in the atmospheric turbulance of simulation, to be transmitted, and from convection type atmospheric turbulence simulation device 108
Second glass window 117 transmits away, and the laser A001 transmitted away from the second glass window 117 incides the second optical splitter
On 109, incided by the laser A001 that the second optical splitter 109 reflects in rectangle diffuse reflecting screen 115；It is continuous with digital camera 113
The hot spot in rectangle diffuse reflecting screen 115 is shot, the light spot image sequence that digital camera 113 is continuously shot to obtain is transmitted to second
Computer system 114；The laser A002 that second laser 112 is sent is incided after the 3rd lens 111 and the 4th lens 110
On second optical splitter 109, from the laser A002 of the transmissive of the second optical splitter 109 from convection type atmospheric turbulence simulation device 108
The second glass window 117 enter in the atmospheric turbulance of simulation and transmit, and from the first of convection type atmospheric turbulence simulation device 108
Glass window 116 is transmitted away, and the laser A002 transmitted away from the first glass window 116 is incided on the first optical splitter 104,
Incided by the laser A002 that the first optical splitter (104) reflects on the 5th lens 105, the laser after the focusing of the 5th lens 105
A002 is incided on ccd detector 106, and the light spot image that ccd detector 106 detects is transmitted to first computer system
107；The distance of ccd detector 106 to the 5th lens 105 is the focal length of the 5th lens 105；Convection type atmospheric turbulance is simulated
Device 108 is the square metal bucket of closing, and atmospheric turbulance is produced in metal bucket during normal work.
1) determine in the case of the atmospheric turbulance in the absence of simulation, the barycenter for the hot spot that ccd detector 106 detects,
Specific method is as follows：
Step Step101：Convection type atmospheric turbulence simulation device 108 is not worked, simulate and fill in convection type atmospheric turbulance
It is tranquil air to put the air in 108 metal bucket；Make the normal work of second laser 112 and launch laser A002；Visit CCD
Survey device 106 and the normal work of first computer system 107；
Step Step102：The frame light spot image B001 that ccd detector 106 is detected is sent to the first department of computer science
In system 107, and light spot image B001 pixel grey scale matrix is stored in the two-dimemsional number in the internal memory of first computer system 107
Group IccdIn, array IccdLine number equal to light spot image B001 pixel grey scale matrix line number, array IccdColumns be equal to light
Spot image B001 pixel grey scale matrix column number；
Step Step103：CalculateWherein Iccd(i1,
j1) represent array IccdI-th1Row, jth1The value of column element, Nr1Represent array IccdLine number, Nc1Represent array IccdRow
Number；irz1And jcz1Line position corresponding to being the barycenter for the hot spot that ccd detector 106 detects is put and column position.
2) digital camera 113 is adjusted and demarcates, specific method is as follows：
Step Step201：Convection type atmospheric turbulence simulation device 108 is not worked, simulate and fill in convection type atmospheric turbulance
It is tranquil air to put the air in 108 metal bucket；Make the normal work of first laser device 101 and launch laser A001；Make numeral
Camera 113 and the normal work of second computer system 114；The position of rectangle diffuse reflecting screen 115 is adjusted, makes to incide rectangle unrestrained anti-
The laser A001 penetrated on screen 115 hot spot is located at the core of rectangle diffuse reflecting screen 115；
Step Step202：First laser device 101 is set not work, with the whole rectangle diffuse reflecting screen 115 of white light；
Step Step203：The mark point of four black is added in rectangle diffuse reflecting screen 115, first mark point is P1,
Second mark point is P2, the 3rd mark point is P3, the 4th mark point is P4, mark point P1、P2、P3、P4Overflow positioned at rectangle
The core of radiation shield 115 simultaneously spaced apart, P1And P2Line perpendicular to P3And P4Line, P1And P2Line put down
Row is in rectangular edges an EV1, P of rectangle diffuse reflecting screen 1153And P4Line parallel to rectangle diffuse reflecting screen 115 another square
Shape side EH1, rectangular edges EV1 is perpendicular to rectangular edges EH1；Digital camera 113 is adjusted, makes connection P1And P2Line segment in digital camera
Same one-row pixels are extended only through in 113 images photographed, make connection P3And P4The image that is photographed in digital camera 113 of line segment
In extend only through same row pixel；
Step Step204：Calculate connection P1And P2The pixel that is passed through in the image that digital camera 113 photographs of line segment
Number NP12, calculate connection P3And P4The number of pixels N that is passed through in the image that digital camera 113 photographs of line segmentP34；Measurement
Connection P in rectangle diffuse reflecting screen 1151And P2Line segment length L12, measure the connection P in rectangle diffuse reflecting screen 1153And P4
Line segment length L34；
Step Step205：Mark point P is removed from rectangle diffuse reflecting screen 1151、P2、P3、P4；
3) determine in the case of the atmospheric turbulance in the absence of simulation, the rectangle diffuse reflecting screen that digital camera 113 photographs
The barycenter of hot spot on 115, specific method are as follows：
Step Step301：Convection type atmospheric turbulence simulation device 108 is not worked, simulate and fill in convection type atmospheric turbulance
It is tranquil air to put the air in 108 metal bucket；Make the normal work of first laser device 101 and launch laser A001；Make numeral
Camera 113 and the normal work of second computer system 114；
Step Step302：The frame light spot image B002 that digital camera 113 is photographed is sent to second computer system
In 114, and light spot image B002 pixel grey scale matrix is stored in the two-dimensional array in the internal memory of second computer system 114
IcamIn, array IcamLine number equal to light spot image B002 pixel grey scale matrix line number, array IcamColumns be equal to hot spot
Image B002 pixel grey scale matrix column number；
Step Step303：CalculateWherein Icam
(i2,j2) represent array IcamI-th2Row, jth2The value of column element, Nr2Represent array IcamLine number, Nc2Represent array Icam's
Columns；irz2And jcz2Position is put and arranged to line position corresponding to being the barycenter of the hot spot for the light spot image B002 that digital camera 113 photographs
Put.
4) while measurement incides the angle of arrival of the laser A002 on ccd detector 106 and incides rectangle diffuse reflecting screen
(115) the beam drift angle of the laser A001 on, specific method are as follows：
Step Step401：Make the normal work of convection type atmospheric turbulence simulation device 108, make 101 normal work of first laser device
Make and launch laser A001；Make the normal work of second laser 112 and launch laser A002；Count ccd detector 106 and first
The normal work of calculation machine system 107；Make digital camera 113 and the normal work of second computer system 114；In first computer system
A new disk file FILE1 is created in 107, a new disk file FILE2 is created in second computer system 114；
Step Step402：The every frame light spot image B001 detected to ccd detector 106, perform following operation：
Step Step402-1：Light spot image B001 is sent in first computer system 107, and light spot image
B001 pixel grey scale matrix is stored in the two-dimensional array I in the internal memory of first computer system 107ccdIn, array IccdRow
The line number of pixel grey scale matrix of the number equal to light spot image B001, array IccdColumns equal to light spot image B001 pixel ash
Spend matrix column number；
Step Step402-2：CalculateWherein Iccd(i1,
j1) represent array IccdI-th1Row, jth1The value of column element, Nr1Represent array IccdLine number, Nc1Represent array IccdRow
Number；
Step Step402-3：Orderβ1=arctan
(dw/f5) value be added to disk file FILE1 end, β1Represent light angle of arrival, f5Represent the focal length length of the 5th lens 105
Degree, Δh1For the pixels tall of ccd detector 106, Δw1For the pixel wide of ccd detector 106；For light spot image B001
Operation terminate；
Step Step403：The every frame light spot image B002 photographed to digital camera 113, perform following operation：
Step Step403-1：Light spot image B002 is sent in second computer system 114, and light spot image
B002 pixel grey scale matrix is stored in the two-dimensional array I in the internal memory of second computer system 114camIn, array IcamRow
The line number of pixel grey scale matrix of the number equal to light spot image B002, array IcamColumns equal to light spot image B002 pixel ash
Spend matrix column number；
Step Step403-2：CalculateWherein Icam
(i2,j2) represent array IcamI-th2Row, jth2The value of column element, Nr2Represent array IcamLine number, Nc2Represent array Icam's
Columns；
Step Step403-3：OrderΔh2=L34/
NP34, Δw2=L12/NP12, β2=arctan (Dw/ F) value be added to disk file FILE2 end, β2Represent light beam drift
Move angle, F=da+db+dc, daRepresent the optical axis distance from first the 116 to the second glass window of glass window 117, dbRepresent from the
The optical axis distance of two the 117 to the second optical splitters of glass window 109, dcRepresent from the second optical splitter 109 to rectangle diffuse reflecting screen
(115) optical axis distance；Operation for light spot image B002 terminates.
5) after the measurement operation execution of the 4) part terminates, file FILE2 is copied to first computer system 107
In, calculate the laser A002 incided on ccd detector 106 angle of arrival and incide the laser in rectangle diffuse reflecting screen 115
The coefficient correlation at A001 beam drift angle, specific method are as follows：
Step Step501：All smooth angle of arrival β being stored in file FILE11Value and be stored in file FILE2
In all beam drift angle betas2Value be loaded into the internal memory of first computer system 107；In first computer system 107
An one-dimension array A is created in internal memoryaoaTo store all smooth angle of arrival β loaded from file FILE11Value, create one
One-dimension array AbwTo store all beam drift angle betas loaded from file FILE22Value；Being stored in file FILE1
M-th of β1Value be stored in AaoaIn [m], Aaoa[m] represents array AaoaM-th of element, m=1,2,3 ..., Num, Num tables
Show the β being stored in file FILE11Number；M-th of β being stored in file FILE22Value be stored in AbwIn [m],
Abw[m] represents array AbwM-th of element, m=1,2,3 ..., Num, Num represents to be stored in β in file FILE22
Number；
Step Step502：Incide the angle of arrival of the laser A002 on ccd detector 106 and incide rectangle diffusing reflection
The coefficient correlation at beam drift angle of laser A001 on screen 115 is：
Wherein
The measurement operation of the 4) part continuously performs TmTerminate after minute, i.e. the continuous probe T of ccd detector 106mMinute, number
Word camera 113 is continuously shot TmMinute；It is utilized respectively first computer system 107 and second computer system 114, step
Step402 and step Step403 starts simultaneously at execution, and ccd detector 106 was detected between two continuous frames light spot image B001 time
Every the time interval that two continuous frames light spot image B002 is shot equal to digital camera 113；The β being stored in file FILE11
Number is equal to the β being stored in file FILE22Number.
As shown in Fig. 2 dotted line represents optical axis, d in figureaRepresent from PATo PBLine segment length, dbRepresent from PBTo PCLine
Segment length, dcRepresent from PCTo PDLine segment length.As shown in figure 3, P1、P2、P3、P4It is to be located at rectangle diffuse reflecting screen core
Four points；As shown in figure 4, the picture element matrix of ccd detector (106) is arranged in rows, the pixel of ccd detector (106)
Height and width mark in figure.
Claims (1)
1. light beam forward direction turbulence transfer drift angle and reverse transfer angle of arrival coefficient correlation measuring method, it is characterised in that：First
The laser A001 that laser is sent is incided on the first optical splitter after the first lens and the second lens, saturating from the first optical splitter
Passed during past laser A001 is penetrated from the first glass window of convection type atmospheric turbulence simulation device into the atmospheric turbulance of simulation
It is defeated, and transmitted away from the second glass window of convection type atmospheric turbulence simulation device, transmitted away from the second glass window
Laser A001 is incided on the second optical splitter, is incided by the laser A001 of the second spectrophotometric reflection in rectangle diffuse reflecting screen；With
Digital camera is continuously shot the hot spot in rectangle diffuse reflecting screen, and the light spot image sequence that digital camera is continuously shot to obtain is transmitted
Give second computer system；The laser A002 that second laser is sent incides second point after the 3rd lens and the 4th lens
On light device, enter from the laser A002 of the second optical splitter transmissive from the second glass window of convection type atmospheric turbulence simulation device
Enter and transmitted in the atmospheric turbulance of simulation, and transmitted away from the first glass window of convection type atmospheric turbulence simulation device, from
The laser A002 that one glass window transmits away is incided on the first optical splitter, incident by the laser A002 of the first spectrophotometric reflection
Onto the 5th lens, the laser A002 after the 5th lens focus is incided on ccd detector, the light that ccd detector detects
Spot image is transmitted to first computer system；The distance of ccd detector to the 5th lens is the focal length of the 5th lens；It is right
Streaming atmospheric turbulence simulation device is the square metal bucket of closing, and atmospheric turbulance is produced in metal bucket during normal work；
1) determine in the case of the atmospheric turbulance in the absence of simulation, the barycenter for the hot spot that ccd detector detects, specific method
It is as follows：
Step Step101：Convection type atmospheric turbulence simulation device is set not work, in the metal bucket of convection type atmospheric turbulence simulation device
Interior air is tranquil air；Make second laser normal work and launch laser A002；Calculate ccd detector and first
Machine system worked well；
Step Step102：The frame light spot image B001 that ccd detector is detected is sent in first computer system, and
Two-dimensional array I light spot image B001 pixel grey scale matrix being stored in the internal memory of first computer systemccdIn, array
IccdLine number equal to light spot image B001 pixel grey scale matrix line number, array IccdColumns equal to light spot image B001
Pixel grey scale matrix column number；
Step Step103：CalculateWherein Iccd(i1,j1) represent
Array IccdI-th1Row, jth1The value of column element, Nr1Represent array IccdLine number, Nc1Represent array IccdColumns；irz1With
jcz1Line position corresponding to being the barycenter for the hot spot that ccd detector detects is put and column position；
2) adjust and demarcate digital camera, specific method is as follows：
Step Step201：Convection type atmospheric turbulence simulation device is set not work, in the metal bucket of convection type atmospheric turbulence simulation device
Interior air is tranquil air；Make first laser device normal work and launch laser A001；Calculate digital camera and second
Machine system worked well；The position of rectangle diffuse reflecting screen is adjusted, makes to incide the hot spot of the laser A001 in rectangle diffuse reflecting screen
Positioned at the core of rectangle diffuse reflecting screen；
Step Step202：First laser device is set not work, with the whole rectangle diffuse reflecting screen of white light；
Step Step203：The mark point of four black is added in rectangle diffuse reflecting screen, first mark point is P1, second mark
Note point is P2, the 3rd mark point is P3, the 4th mark point is P4, mark point P1、P2、P3、P4Positioned at rectangle diffuse reflecting screen
Core simultaneously spaced apart, P1And P2Line perpendicular to P3And P4Line, P1And P2Line overflow parallel to rectangle
One rectangular edges EV1, P of radiation shield3And P4Line parallel to another rectangular edges EH1 of rectangle diffuse reflecting screen, rectangular edges
EV1 is perpendicular to rectangular edges EH1；Digital camera is adjusted, makes connection P1And P2Line segment only worn in the image that digital camera photographs
Same one-row pixels are crossed, make connection P3And P4Line segment same row pixel is extended only through in the image that digital camera photographs；
Step Step204：Calculate connection P1And P2The number of pixels N that is passed through in the image that digital camera photographs of line segmentP12,
Calculate connection P3And P4The number of pixels N that is passed through in the image that digital camera photographs of line segmentP34；Measure rectangle diffuse reflecting screen
On connection P1And P2Line segment length L12, measure the connection P in rectangle diffuse reflecting screen3And P4Line segment length L34；
Step Step205：Mark point P is removed from rectangle diffuse reflecting screen1、P2、P3、P4；
3) determine in the case of the atmospheric turbulance in the absence of simulation, the hot spot in rectangle diffuse reflecting screen that digital camera photographs
Barycenter, specific method is as follows：
Step Step301：Convection type atmospheric turbulence simulation device is set not work, in the metal bucket of convection type atmospheric turbulence simulation device
Interior air is tranquil air；Make first laser device normal work and launch laser A001；Calculate digital camera and second
Machine system worked well；
Step Step302：The frame light spot image B002 that digital camera is photographed is sent in second computer system, and handle
Light spot image B002 pixel grey scale matrix is stored in the two-dimensional array I in the internal memory of second computer systemcamIn, array Icam
Line number equal to light spot image B002 pixel grey scale matrix line number, array IcamColumns be equal to light spot image B002 picture
The columns of plain gray matrix；
Step Step303：CalculateWherein Icam(i2,j2)
Represent array IcamI-th2Row, jth2The value of column element, Nr2Represent array IcamLine number, Nc2Represent array IcamColumns；
irz2And jcz2Line position corresponding to being the barycenter of the hot spot for the light spot image B002 that digital camera photographs is put and column position；
4) while measurement incides the angle of arrival of the laser A002 on ccd detector and incides the laser in rectangle diffuse reflecting screen
A001 beam drift angle, specific method are as follows：
Step Step401：Make convection type atmospheric turbulence simulation device normal work, make first laser device normal work and launch to swash
Light A001；Make second laser normal work and launch laser A002；Make ccd detector and the normal work of first computer system
Make；Make digital camera and second computer system worked well；A new disk file is created in first computer system
FILE1, a new disk file FILE2 is created in second computer system；
Step Step402：The every frame light spot image B001 detected to ccd detector, perform following operation：
Step Step402-1：Light spot image B001 is sent in first computer system, and light spot image B001 pixel
Gray matrix is stored in the two-dimensional array I in the internal memory of first computer systemccdIn, array IccdLine number be equal to light spot image
The line number of B001 pixel grey scale matrix, array IccdColumns be equal to light spot image B001 pixel grey scale matrix column number；
Step Step402-2：CalculateWherein Iccd(i1,j1)
Represent array IccdI-th1Row, jth1The value of column element, Nr1Represent array IccdLine number, Nc1Represent array IccdColumns；
Step Step402-3：Orderβ1=arctan (dw/f5)
Value be added to disk file FILE1 end, β1Represent light angle of arrival, f5Represent the focal length of the 5th lens, Δh1For
The pixels tall of ccd detector, Δw1For the pixel wide of ccd detector；Operation for light spot image B001 terminates；
Step Step403：The every frame light spot image B002 photographed to digital camera, perform following operation：
Step Step403-1：Light spot image B002 is sent in second computer system, and light spot image B002 pixel
Gray matrix is stored in the two-dimensional array I in the internal memory of second computer systemcamIn, array IcamLine number be equal to light spot image
The line number of B002 pixel grey scale matrix, array IcamColumns be equal to light spot image B002 pixel grey scale matrix column number；
Step Step403-2：CalculateWherein Icam(i2,
j2) represent array IcamI-th2Row, jth2The value of column element, Nr2Represent array IcamLine number, Nc2Represent array IcamRow
Number；
Step Step403-3：OrderΔh2=L34/NP34,
Δw2=L12/NP12, β2=arctan (Dw/ F) value be added to disk file FILE2 end, β2Beam drift angle is represented,
F=da+db+dc, daRepresent from the first glass window to the optical axis distance of the second glass window, dbRepresent from the second glass window to
The optical axis distance of second optical splitter, dcRepresent from the second optical splitter to the optical axis distance of rectangle diffuse reflecting screen；For light spot image
B002 operation terminates；
5) after the measurement operation execution of the 4) part terminates, file FILE2 is copied in first computer system, calculate into
The angle of arrival for the laser A002 being mapped on ccd detector and the beam drift angle for inciding the laser A001 in rectangle diffuse reflecting screen
Coefficient correlation, specific method is as follows：
Step Step501：All smooth angle of arrival β being stored in file FILE11Value and the institute that is stored in file FILE2
There is beam drift angle beta2Value be loaded into the internal memory of first computer system；One is created in the internal memory of first computer system
Individual one-dimension array AaoaTo store all smooth angle of arrival β loaded from file FILE11Value, create an one-dimension array AbwCome
Store all beam drift angle betas loaded from file FILE22Value；M-th of β being stored in file FILE11Value protect
A be presentaoaIn [m], Aaoa[m] represents array AaoaM-th of element, m=1,2,3 ..., Num, Num represent be stored in file
β in FILE11Number；M-th of β being stored in file FILE22Value be stored in AbwIn [m], Abw[m] represents array
AbwM-th of element, m=1,2,3 ..., Num, Num represents to be stored in β in file FILE22Number；
Step Step502：Incide the angle of arrival of the laser A002 on ccd detector and incide swashing in rectangle diffuse reflecting screen
The coefficient correlation at light A001 beam drift angle is：
Wherein
The measurement operation of the 4) part continuously performs TmTerminate after minute, i.e. ccd detector continuous probe TmMinute, digital camera
It is continuously shot TmMinute；It is utilized respectively first computer system and second computer system, step Step402 and step
Step403 starts simultaneously at execution, and ccd detector detection two continuous frames light spot image B001 time interval is clapped equal to digital camera
Take the photograph two continuous frames light spot image B002 time interval；The β being stored in file FILE11Number be equal to be stored in file
β in FILE22Number.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711174454.0A CN107741280B (en) | 2017-11-22 | 2017-11-22 | Light beam forward direction turbulence transfer drift angle and reverse transfer angle of arrival related coefficient measurement method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711174454.0A CN107741280B (en) | 2017-11-22 | 2017-11-22 | Light beam forward direction turbulence transfer drift angle and reverse transfer angle of arrival related coefficient measurement method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107741280A true CN107741280A (en) | 2018-02-27 |
| CN107741280B CN107741280B (en) | 2019-05-17 |
Family
ID=61239160
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711174454.0A Active CN107741280B (en) | 2017-11-22 | 2017-11-22 | Light beam forward direction turbulence transfer drift angle and reverse transfer angle of arrival related coefficient measurement method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107741280B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105046724A (en) * | 2015-07-21 | 2015-11-11 | 长春理工大学 | Dual-image-point centroid extraction method in atmospheric coherent parameter test based on differential image motion |
| CN105610493A (en) * | 2015-12-21 | 2016-05-25 | 西安空间无线电技术研究所 | Atmosphere turbulence simulation system and method based on inverse self-adaptation technology |
| CN106059678A (en) * | 2015-04-13 | 2016-10-26 | 波音公司 | System for bidirectional free-space laser communication of gigabit Ethernet telemetry data |
-
2017
- 2017-11-22 CN CN201711174454.0A patent/CN107741280B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106059678A (en) * | 2015-04-13 | 2016-10-26 | 波音公司 | System for bidirectional free-space laser communication of gigabit Ethernet telemetry data |
| CN105046724A (en) * | 2015-07-21 | 2015-11-11 | 长春理工大学 | Dual-image-point centroid extraction method in atmospheric coherent parameter test based on differential image motion |
| CN105610493A (en) * | 2015-12-21 | 2016-05-25 | 西安空间无线电技术研究所 | Atmosphere turbulence simulation system and method based on inverse self-adaptation technology |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107741280B (en) | 2019-05-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106679940B (en) | A kind of high-precision laser angle of divergence parameter calibration device | |
| CN104567738B (en) | Parallelism of optical axis accurate measuring systems and method | |
| CN105865755B (en) | A kind of display device measuring device and measuring method of simulation human eyes structure | |
| CN105910712B (en) | Five-channel self-adaption two-dimensional temperature field measuring apparatus and its measurement method | |
| CN109632104A (en) | A kind of infrared array radiation source calibrating installation and calibration method | |
| CN102853916B (en) | Method and system for conducting remote infrared temperature measurement on coal pile surfaces | |
| CN104848801B (en) | A kind of line structured light vision sensor calibration method based on parallel bicylindrical target | |
| CN103994719A (en) | High-precision three-dimensional imaging device based on Geiger APD arrays and using method thereof | |
| CN109883391B (en) | Monocular distance measurement method based on digital imaging of microlens array | |
| CN107422392A (en) | It is a kind of to be positioned and tracing system and method around angle based on single photon detection | |
| CN101520343A (en) | Assembling and aligning device and method for thermal infrared spectrum imaging system | |
| CN102564611A (en) | High-power laser wave front measuring instrument and wave front measuring method | |
| Settles | Smartphone schlieren and shadowgraph imaging | |
| CN101706951B (en) | Method, device and system for objectively evaluating pneumatic optical image quality based on feature fusion | |
| Li et al. | Correction model for microlens array assembly error in light field camera | |
| CN203011419U (en) | An optical axis parallelism digital calibration instrument for multiple optical sensors | |
| CN105928456A (en) | Optical detection device, system and method | |
| CN103778637B (en) | A kind of in-orbit 90 ° of radiation calibration processing methods based on histogram balance | |
| CN107741280A (en) | Light beam forward direction turbulence transfer drift angle and reverse transfer angle of arrival coefficient correlation measuring method | |
| CN108007387B (en) | Surface shape measurement device and method based on Structured Illumination | |
| CN104267002B (en) | Improve transmission-type visibility meter and the measuring method thereof of accuracy of measurement | |
| Growcott et al. | Development and assessment of a digital stereo photogrammetric system to measure cetaceans at sea | |
| CN102998093A (en) | Device and method for quickly automatically measuring focal distance | |
| CN207783024U (en) | The EMCCD blind elements that collimator is homogenized based on light test system | |
| Haussmann et al. | Streak detection of space debris by a passive optical sensor |
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 |