CN113340856B - Noise removal algorithm of optical sensor and visibility calculation method thereof - Google Patents

Noise removal algorithm of optical sensor and visibility calculation method thereof Download PDF

Info

Publication number
CN113340856B
CN113340856B CN202110620068.XA CN202110620068A CN113340856B CN 113340856 B CN113340856 B CN 113340856B CN 202110620068 A CN202110620068 A CN 202110620068A CN 113340856 B CN113340856 B CN 113340856B
Authority
CN
China
Prior art keywords
light
infrared
red
photons
red light
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
Application number
CN202110620068.XA
Other languages
Chinese (zh)
Other versions
CN113340856A (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.)
Wuhan Zhiteng Technology Co ltd
Original Assignee
Wuhan Zhiteng Technology Co ltd
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 Wuhan Zhiteng Technology Co ltd filed Critical Wuhan Zhiteng Technology Co ltd
Priority to CN202110620068.XA priority Critical patent/CN113340856B/en
Publication of CN113340856A publication Critical patent/CN113340856A/en
Application granted granted Critical
Publication of CN113340856B publication Critical patent/CN113340856B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/59Transmissivity
    • G01N21/61Non-dispersive gas analysers

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Photometry And Measurement Of Optical Pulse Characteristics (AREA)

Abstract

The invention relates to a noise removal algorithm of an optical sensor, which is realized by the following steps: step one, firstly, calculating the actual number of received photons
Figure DDA0003099463080000011
Step two, utilizing a mathematical statistical method to compare the number of received photons C in step one i And the output light intensity value I i Linear regression is carried out, B is obtained by calculation i =C r0 (ii) a Step three, calculating the true value of the number of received photons of the receiving device of the photoelectric sensor
Figure DDA0003099463080000012
According to the noise removal algorithm for the photoelectric sensor, the true value of the number of received photons is used for calculating the possible detection device more accurately when the numerical value of the photoelectric sensor is calculated, so that the detection device can remove the numerical value of noise generated by changing the ambient temperature when the numerical value of the photoelectric sensor is calculated. The noise removal algorithm of the present invention can be applied to various detection devices using a photosensor.

Description

Noise removal algorithm of optical sensor and visibility calculation method thereof
Technical Field
The invention relates to the field of optical detection, in particular to a noise removal algorithm of an optical sensor and a computing method for computing visibility by adopting the noise removal algorithm.
Background
As is well known, different types of sensors are used for collecting measurement data in some detection devices, wherein the principle of a photoelectric sensor is that light emitted by an emitting device of the photoelectric sensor is attenuated by air, is received by a receiving device of the photoelectric sensor and converted into an electric signal, and then is received by a detection circuit such as an input circuit and an operational amplifier and finally output as photon number, and the obtained photon number is used for calculating the result of the detection device. When the photoelectric sensor is used in a detection device, the noise generated by the ambient temperature can bring certain influence on detection data, and cause inaccurate detection. Especially, if this optical sensor is used for visibility to detect time, can make the detection of visibility less accurate when the data that optical sensor obtained is inaccurate to bring the potential safety hazard for driving.
Therefore, the present inventors have made extensive studies on the above problems, and as a result, the present invention has been made.
Disclosure of Invention
The invention aims to provide a noise removal algorithm of an optical sensor, which eliminates the influence of noise generated by environment temperature change on the detection of the optical sensor.
Another object of the present invention is to provide a visibility calculation method using an optical sensor that can accurately detect visibility.
To achieve the above object, the solution of the present invention is as follows: a noise removal algorithm for a light sensor is implemented by the following steps:
step one, inputting n (n) to a transmitting device of a photoelectric sensor>3) The linear current value makes the emitting device of the photoelectric sensor output different output light intensity I i (I =1,2,3 …, n) at n different intensity values I i Next, the receiving device of the photoelectric sensor actually receives different current values A i (i =1,2,3 …, n) and results in a different number of received photons, C i (i =1,2,3 …, n), dark current due to ambient temperature change is defined as a T Then, the receiving current value of the receiving device of the photoelectric sensor
Figure BDA0003099463060000021
(i =1,2,3 …, n), wherein
Figure BDA0003099463060000022
Is the true value of the receiving current obtained when the light of the emitting device is irradiated to the light receiving device under the current input light intensity, and the actual number of the received photons of the receiving device
Figure BDA0003099463060000023
(i=1,2,3 …, n), wherein
Figure BDA0003099463060000024
Is based on the true value of the received current
Figure BDA0003099463060000025
Generating a true value of the number of received photons, C r0 Is caused by dark current A T The number of generated photons is the noise value generated by the change of the environmental temperature;
step two, utilizing a mathematical statistical method to compare the number of received photons C in step one i And the output light intensity value I i Linear regression is performed, and C is set i =K i I i +B i The slope of the fitted line is
Figure BDA0003099463060000026
Figure BDA0003099463060000027
And step one
Figure BDA0003099463060000028
I.e. B can be calculated i =C r0;
Step three, when I i The value of (a) is 0,
Figure BDA0003099463060000029
the number of received photons C of the receiving means of the photosensor i =B i =C r0 In this case, since the dark current is constant in a short time and the number of photons due to noise is constant, let B i =0, namely the true value of the number of received photons of the receiving device of the photoelectric sensor
Figure BDA00030994630600000210
A visibility calculation method using an optical sensor is realized by the following steps:
the method comprises the following steps: the detection system is installed, a detection device, a convex lens, a red light filter and a reflector are sequentially arranged at a detection position, the detection device, the convex lens, the red light filter and the reflector are sequentially arranged on the same straight line at intervals, the convex lens, the red light filter and the reflector are vertically arranged, the detection device is provided with two light emitters and two light detectors, one light emitter is a red light emitter capable of emitting red light, the other light emitter is an infrared light emitter capable of emitting infrared light, one light detector is a red light detector capable of receiving red light, the other light detector is an infrared light detector capable of receiving infrared light, the emitting direction of the red light emitter, the emitting direction of the infrared light emitter, the receiving direction of the red light detector and the receiving direction of the infrared light detector all face towards the convex lens, the distance between the detection device and the convex lens is 3-10cm, the distance between the convex lens and the red light filter is 15-25cm, the visible light filter and the reflector are arranged close to each other, and the reflector is an infrared light reflector for reflecting infrared light back to the infrared light in the original way;
step two, the noise removal algorithm is specifically realized by the following steps:
step a) inputting n (n) to the infrared emitter and the red light emitter respectively>3) A linearly varying current value to make the infrared emitter and the red emitter output different output light intensities I i (I =1,2,3 …, n) at n different intensity values I i Next, the infrared light detector receives different current values A1 i (i =1,2,3 …, n), and different numbers of infrared light receiving photons C1 are obtained i (i =1,2,3 …, n), the red light detector receives different current values A2 i (i =1,2,3 …, n) and a different number of red-light-receiving photons C2 is obtained i (i =1,2,3 …, n), dark current due to environmental temperature change, where a is assumed to be dark current T Then, the received current value of the infrared light detector
Figure BDA0003099463060000031
(i =1,2,3 …, n) wherein
Figure BDA0003099463060000032
Is received when the light of the infrared light emitter irradiates the infrared light detector under the current input light intensityTrue current value, number of infrared receiving photons of infrared light detector
Figure BDA0003099463060000033
(i =1,2,3 …, n), wherein
Figure BDA0003099463060000034
Is based on the true value of the received current
Figure BDA0003099463060000035
True value of the number of infrared photons generated, C1 r0 Is infrared light in dark current A T The number of infrared photons generated at the lower part is the noise value generated by the change of the environmental temperature; received current value of red light detector
Figure BDA0003099463060000036
(i =1,2,3 …, n), wherein
Figure BDA0003099463060000037
Is the true value of the receiving current obtained when the light of the red light emitter irradiates the red light detector under the current input light intensity, and the number of the red light receiving photons of the red light detector
Figure BDA0003099463060000038
(i =1,2,3 …, n) wherein
Figure BDA0003099463060000039
Is based on the true value of the received current
Figure BDA00030994630600000310
True value of the number of red-light-receiving photons, C2, is generated r0 Is red light at dark current A T The number of red photons generated;
step b), respectively counting the number of infrared light receiving photons C1 in the step a) by using a mathematical statistical method i And the output light intensity value I i And the number of red light-receiving photons C2 i And the output light intensity value I i Linear regression is performed, and C1 is set i =K1 i I i +B1 i ,C2 i =K2 i I i +B2 i The slope of the fitted line of infrared light is
Figure BDA0003099463060000041
In step a)
Figure BDA0003099463060000042
B1 can be calculated i =C1 r0 (ii) a Similarly, the slope of the fitted line of red light is
Figure BDA0003099463060000043
In step a)
Figure BDA0003099463060000044
B2 can be calculated i =C2 r0
Step c) when I i The value of (a) is 0,
Figure BDA0003099463060000045
the number of infrared receiving photons of the infrared light detector C1 i =B1 i =C1 r0 The number of red light receiving photons C2 of the red light detector i =B2 i =C2 r0 Since the dark current is constant in a short time and the number of photons due to noise is constant, B1 is set to i =0, namely the true value of the infrared light receiving photon number of the infrared light detector
Figure BDA0003099463060000046
In the same way, let B2 i =0, the true value of the red light receiving photon number of the red light detector is obtained
Figure BDA0003099463060000047
Step three, eliminating the infrared back scattering light, because the infrared light receiving photon number true value of the infrared light detector
Figure BDA00030994630600000414
The number of infrared backward radiation photons is C1S i And the number of infrared transmission photons is CT i Then, then
Figure BDA0003099463060000049
Figure BDA00030994630600000410
And in step c)
Figure BDA00030994630600000411
Namely CT i =K1 i I i -C1S i Furthermore, the red light receives the photon number true value due to the arrangement of the red light filter
Figure BDA00030994630600000412
Let the number of red backscattered photons be C2S i I.e. by
Figure BDA00030994630600000413
Meanwhile, setting the distance between the convex lens and the red light filter as d, and the backscattering energy equation p (d) of the light detector at the distance d is as follows:
Figure BDA0003099463060000051
in the formula, p 0 Is the emission power of the photodetector, c is the speed of light, τ is the pulse width, A d For the effective receiving area, Y (d) is the optical characteristic of the photodetector, β (d) is the backscattering coefficient, T d Is the transmission factor; according to the transmission factor of light in the aerosol particles
Figure BDA0003099463060000052
Wherein I is the transmitted light intensity, I 0 Is the incident light intensity, this incident light intensity I 0 The output light intensity I of the step two i
The conversion of p (d) at this distance d into the number of photons C takes the form
Figure BDA0003099463060000053
Eta is the quantity of the light detectorThe sub-efficiency, λ is the wavelength of light, h is the Planck constant, and t is the time between the emission of the light emitter and the reception of the light detector; the number of red backscattered photons C2S at the distance d can be known i (d) And the number of infrared backscattered photons C1S i (d) In a ratio of
Figure BDA0003099463060000054
Figure BDA0003099463060000055
In the formula, λ R Wavelength of red light, λ IR Is the wavelength of infrared light, beta R ) d) is the red backscattering coefficient, beta IR (d) Is the infrared light backscattering coefficient, T rR Is the red light transmission factor, T rIR For the infrared light transmission factor, the infrared backscattered photon number C1S can be obtained by using the formula of the ratio of red backscattered photon number to infrared backscattered photon number i By using
Figure BDA0003099463060000056
The infrared transmission photon number CT can be obtained by the formula i
Step four, according to the infrared transmission photon number CT obtained in step three i The transmission factor of infrared light can be obtained
Figure BDA0003099463060000057
Wherein C is 0 At a light intensity value I 0 The number of received photons of the infrared light detector when the lower light propagation distance is 0; then the Beer-Bouguer-Lambert theorem T rIR =e -σ2d Obtaining extinction coefficient sigma of infrared light, where e is the base of natural number logarithm and 2d is optical path, and substituting the obtained extinction coefficient sigma into visibility calculation formula
Figure BDA0003099463060000058
The visibility can be obtained.
By adopting the technical scheme, the noise removal algorithm of the optical sensor provided by the invention can be used for calculating the value of noise generated due to the change of the environmental temperature in the measured data of the optical sensor to obtain the true test value of the optical sensor, eliminating the noise generated due to the environmental temperature, ensuring that the test data of the optical sensor is calculated by the detection device more accurately, and improving the detection accuracy of the detection device.
According to the visibility calculation method adopting the optical sensor, the number of backward radiation photons of infrared light and the noise value generated by the change of the ambient temperature can be considered during the visibility calculation, so that the visibility calculation accuracy is greatly improved.
Drawings
FIG. 1 is a schematic structural diagram of the present invention.
Detailed Description
In order to further explain the technical solution of the present invention, the present invention is explained in detail by the following specific examples.
A noise removal algorithm for a light sensor is implemented by the following steps:
step one, inputting n (n) to a transmitting device of a photoelectric sensor>3) The linear current value makes the emitting device of the photoelectric sensor output different output light intensity I i (I =1,2,3 …, n) at n different intensity values I i Next, the receiving device of the photoelectric sensor actually receives different current values A i (i =1,2,3 …, n) and results in a different number of received photons, C i (i =1,2,3 …, n), dark current due to environmental temperature change, where a is assumed to be dark current T Then, the receiving current value of the receiving device of the photoelectric sensor
Figure BDA0003099463060000061
(i =1,2,3 …, n), wherein
Figure BDA0003099463060000062
Is the true value of the receiving current obtained when the light of the emitting device is irradiated to the light receiving device under the current input light intensity, and the actual number of the received photons of the receiving device
Figure BDA0003099463060000063
(i =1,2,3 …, n), wherein
Figure BDA0003099463060000064
Is based on the true value of the received current
Figure BDA0003099463060000065
Generating a true value of the number of received photons, C r0 Is caused by dark current A T The number of generated photons is the noise value generated by the change of the environmental temperature;
step two, utilizing a mathematical statistical method to compare the number of the received photons C in the step one i And the output light intensity value I i Linear regression is performed, and C is set i =K i I i +B i The slope of the fitted line is
Figure BDA0003099463060000066
Figure BDA0003099463060000071
And step one
Figure BDA0003099463060000072
I.e. B can be calculated i =C r0; Implementation of this step reduces the effect of white noise;
step three, when I i The value of (a) is 0,
Figure BDA0003099463060000073
the number of received photons C of the receiving means of the photosensor i =B i =C r0 In this case, since the dark current is constant in a short time and the number of photons due to noise is constant, let B i =0, namely the true value of the number of received photons of the receiving device of the photoelectric sensor
Figure BDA0003099463060000074
By adopting the steps, the true value of the number of received photons of the photoelectric sensor can be obtained, so that the detection device can accurately calculate the possible detection device by utilizing the true value of the number of received photons when calculating the numerical value of the photoelectric sensor, and the detection device can remove the numerical value of noise caused by the change of the ambient temperature when calculating the numerical value of the photoelectric sensor. The noise removal algorithm of the present invention can be applied to various detection devices using a photosensor.
The invention relates to a noise removal algorithm of an optical sensor, comprising the following steps: the detection system is installed, as shown in fig. 1, a detection device 1, a convex lens 2, a red light filter 3 and a reflector 4 are arranged at a detection position, the detection device 1, the convex lens 2, the red light filter 3 and the reflector 4 are sequentially arranged on the same straight line at intervals, the convex lens 2, the red light filter 3 and the reflector 4 are vertically arranged, the detection device 1 is provided with two light emitters and two light detectors, one light emitter is a red light emitter (not shown) capable of emitting light 100, the other light emitter is an infrared emitter (not shown) capable of emitting infrared light 200, one light detector is a red light detector for receiving red light, the other light detector is an infrared light detector for receiving infrared light, the red light emitter, the infrared emitter, the red light detector and the infrared light detector are all known, the emitting direction of the red light emitter, the emitting direction of the infrared emitter, the receiving direction of the red light detector and the receiving direction of the infrared light detector all face the convex lens 2, the distance between the detection device 1 and the convex lens 2 is 3-10cm, preferably, the emitting end of the infrared emitter, the receiving direction of the red light detector and the receiving direction of the infrared light detector are all facing the convex lens 2, the reflector 4, preferably, the distance between the detecting device is a distance between the infrared light filter 4 and the reflector 4, the reflector 4 is preferably, the distance between the reflector 4 is preferably, the reflector 4 is set for detecting device is set for detecting the reflector 4 and the reflector 4, and the reflector 4;
step two, the noise removal algorithm is specifically realized by the following steps:
step a) inputting n (n) to the infrared emitter and the red light emitter respectively>3) A linearly varying current value to make the infrared emitter and the red emitter output different output light intensities I i (I =1,2,3 …, n) at n different intensity values I i Next, the infrared light detector receives different current values A1 i (i =1,2,3 …, n), and different numbers of infrared light receiving photons C1 are obtained i (i =1,2,3 …, n), the red light detector receives different current values A2 i (i =1,2,3 …, n) and results in a different number of red light-receiving photons C2 i (i =1,2,3 …, n), dark current due to ambient temperature change is defined as a T Then, the received current value of the infrared light detector
Figure BDA0003099463060000081
(i =1,2,3 …, n), wherein
Figure BDA0003099463060000082
The receiving current true value is obtained when the light of the infrared light emitter irradiates the infrared light detector under the current input light intensity, and the infrared light receiving photon number of the infrared light detector
Figure BDA0003099463060000083
(i =1,2,3 …, n), wherein
Figure BDA0003099463060000084
Is based on the true value of the received current
Figure BDA0003099463060000085
True value of the number of infrared photons generated, C1 r0 Is infrared light in dark current A T The number of infrared photons generated at the lower part is the noise value generated by the change of the environmental temperature; received current value of red light detector
Figure BDA0003099463060000086
(i =1,2,3 …, n), wherein
Figure BDA0003099463060000087
Is the true value of the receiving current obtained when the light of the red light emitter irradiates the red light detector under the current input light intensity, and the number of the red light receiving photons of the red light detector
Figure BDA0003099463060000088
(i =1,2,3 …, n), wherein
Figure BDA0003099463060000089
Is based on the true value of the received current
Figure BDA00030994630600000810
True value of the number of red-light-receiving photons, C2, is generated r0 Is red light in dark current A T The number of red photons generated;
step b), respectively counting the number of infrared light receiving photons C1 in the step a) by using a mathematical statistical method i And the output light intensity value I i And the number of red light-receiving photons C2 i And the output light intensity value I i Linear regression is performed, and C1 is set i =K1 i I i +B1 i ,C2 i =K2 i I i +B2 i The slope of the fitted line of infrared light is
Figure BDA0003099463060000091
In step a)
Figure BDA0003099463060000092
B1 can be calculated i =C1 r0; Similarly, the slope of the fitted line of red light is
Figure BDA0003099463060000093
In step a)
Figure BDA0003099463060000094
B2 can be calculated i =C2 r0
Step c) when I i The value of (a) is 0,
Figure BDA0003099463060000095
the number of infrared receiving photons of the infrared light detector C1 i =B1 i =C1 r0 Number of red light receiving photons C2 of red light detector i =B2 i =C2 r0 In this case, since the dark current is constant in a short time and the number of photons due to noise is constant, B1 is set i =0, namely the true value of the infrared light receiving photon number of the infrared light detector
Figure BDA0003099463060000096
In the same way, let B2 i =0, the true value of the red light receiving photon number of the red light detector is obtained
Figure BDA0003099463060000097
Step three, eliminating the infrared back scattering light, because the infrared light receiving photon number true value of the infrared light detector
Figure BDA0003099463060000098
The number of infrared backward radiation photons is C1S i The number of infrared transmission photons is CT i Then, then
Figure BDA0003099463060000099
Figure BDA00030994630600000910
And in step c)
Figure BDA00030994630600000911
Namely CT i =K1 i I i -C1S i Furthermore, the red light receives the photon number true value due to the arrangement of the red light filter
Figure BDA00030994630600000912
Let the number of red backscattered photons be C2S i I.e. by
Figure BDA00030994630600000913
Meanwhile, setting the distance between the convex lens and the red light filter as d, and the backscattering energy equation p (d) of the light detector at the distance d is as follows:
Figure BDA00030994630600000914
in the formula, p 0 Is the emission power of the photodetector, c is the speed of light, τ is the pulse width, A d For the effective receiving area, Y (d) is the optical characteristics of the photodetector, β (d) is the backscattering coefficient, T d Is the transmission factor; according to the transmission factor of light in the aerosol particles
Figure BDA0003099463060000101
Wherein I is the transmitted light intensity, I 0 Is the incident light intensity, I 0 The output light intensity I of the step two i
The conversion of p (d) at this distance d into the number of photons C is in the form of
Figure BDA0003099463060000102
Eta is the quantum efficiency of the optical detector, lambda is the wavelength of light, h is the Planck constant, and t is the time between the emission of the optical transmitter and the reception of the optical detector; the number of red backscattered photons C2S at the distance d can be known i (d) And the number of infrared backscattered photons C1S i (d) In a ratio of
Figure BDA0003099463060000103
Figure BDA0003099463060000104
In the formula of lambda R Wavelength of red light, λ IR Is the wavelength of infrared light, beta R (d) Is the red light backscattering coefficient, beta IR (d) Is the infrared light backscattering coefficient, T rR Is the red light transmission factor, T rIR For the infrared light transmission factor, the infrared backscattered photon number C1S can be obtained by using the formula of the ratio of red backscattered photon number to infrared backscattered photon number i By using
Figure BDA0003099463060000105
The infrared transmission photon number CT can be obtained by the formula i
Step four, according to the infrared transmitted photon number CT obtained in step three i The transmission factor of infrared light can be obtained
Figure BDA0003099463060000106
Wherein C 0 At a light intensity value I 0 The number of received photons of the infrared light detector when the lower light propagation distance is 0; then the Beer-Bouguer-Lambert theorem T rIR =e -σ2d Obtaining extinction coefficient sigma of infrared light, where e is the base of natural number logarithm and 2d is optical path, and substituting the obtained extinction coefficient sigma into visibility calculation formula
Figure BDA0003099463060000107
The visibility can be obtained.
The visibility calculation method adopting the optical sensor comprises the steps that when the visibility calculation method is applied, a detection device emits red light and infrared light, the red light and the infrared light pass through a convex lens 2 to a red light filter 3, at the moment, the red light and the infrared light are influenced by particles in the air in the process of transmitting to the red light filter, a small amount of backward scattering of the red light and the infrared light occurs, most of the red light is filtered by the red light filter 3 and cannot be reflected to an infrared light detector, the number of photons received by the red light detector is the number of photons of red light backward radiation, the infrared light passes through the red light filter 3 to a reflector 4 and is reflected to the detection device by a reflector 4 in an original way, the number of photons received by the infrared light detector at the detection device is the sum of the number of infrared transmission photons and the number of infrared light backward radiation photons, the number of backward radiation photons can be calculated when the visibility is calculated, the problem that backward radiation is not considered in the traditional visibility calculation is solved, the calculation accuracy of the visibility calculation is greatly improved, and the noise value generated by the ambient temperature is further removed when the visibility calculation is more accurate.
In the visibility calculation method, the light reflecting piece in the first step is a cylinder, one surface of the light reflecting piece facing the red light filter is adhered with a 3MM light reflecting film, preferably the cylinder is a prism, infrared light at different angles can be reflected back by utilizing the matching of the 3MM light reflecting film and the prism, and the visibility calculation accuracy is further ensured.
The above examples and drawings are not intended to limit the scope of the present invention, and any suitable changes or modifications thereof by one of ordinary skill in the art should be considered as not departing from the scope of the present invention.

Claims (1)

1. A visibility calculation method using an optical sensor is characterized by comprising the following steps:
the method comprises the following steps: the detection system is installed, a detection device, a convex lens, a red light filter and a light reflecting piece are sequentially arranged at a detection position, the detection device, the convex lens, the red light filter and the light reflecting piece are sequentially arranged on the same straight line at intervals, the convex lens, the red light filter and the light reflecting piece are vertically arranged, the detection device is provided with two light emitters and two light detectors, one light emitter is a red light emitter capable of emitting red light, the other light emitter is an infrared emitter capable of emitting infrared light, one light detector is a red light detector capable of receiving red light, the other light detector is an infrared light detector capable of receiving infrared light, the emitting direction of the red light emitter, the emitting direction of the infrared emitter, the receiving direction of the red light detector and the receiving direction of the infrared light detector all face towards the convex lens, the distance between the detection device and the convex lens is 3-10cm, the distance between the convex lens and the red light filter is 15-25cm, the red light filter and the light reflecting piece are arranged close to each other, and the light reflecting piece is an infrared light reflecting piece capable of reflecting infrared light back to the infrared light in the original way;
step two, the noise removal algorithm is specifically realized by the following steps:
step a) inputting n current values which change linearly to the infrared emitter and the red light emitter respectively, wherein n>3, making the infrared emitter and the red emitter output different output light intensity I i I =1,2,3 …, n, at n different intensity values I i Next, the infrared light detector receives different current values A1 i I =1,2,3 …, n, and obtains different numbers of infrared light receiving photons C1 i I =1,2,3 …, n, the red light detector receives different current values A2 i I =1,2,3 …, n, and yields different numbers of red light receiving photons C2 i I =1,2,3 …, n, dark current due to ambient temperature change, where a is the dark current T Then, the first step is executed,received current value of infrared light detector
Figure FDA0003897114830000011
Wherein
Figure FDA0003897114830000012
Is the true value of the receiving current obtained when the light of the infrared light emitter irradiates the infrared light detector under the current input light intensity, and the number of the infrared light receiving photons of the infrared light detector
Figure FDA0003897114830000013
Figure FDA0003897114830000021
Wherein
Figure FDA0003897114830000022
Is based on the true value of the received current
Figure FDA0003897114830000023
The generated infrared photon number true value, C1 r0 Is infrared light in dark current A T The number of infrared photons generated is the noise value generated by the change of the environmental temperature; received current value of red light detector
Figure FDA0003897114830000024
Figure FDA0003897114830000025
Wherein
Figure FDA0003897114830000026
Is the true value of the receiving current obtained when the light of the red light emitter irradiates the red light detector under the current input light intensity, and the number of the red light receiving photons of the red light detector
Figure FDA0003897114830000027
Wherein
Figure FDA0003897114830000028
Is based on receiving current true value
Figure FDA0003897114830000029
True value of the number of red-light-receiving photons, C2, is generated r0 Is red light at dark current A T The number of red photons generated;
step b), respectively comparing the infrared light receiving photon number C1 in the step a) by using a mathematical statistical method i And the output light intensity value I i And the number of red light-receiving photons C2 i And the output light intensity value I i Linear regression is performed, and C1 is set i =K1 i I i +B1 i ,C2 i =K2 i I i +B2 i The slope of the fitted line of infrared light is
Figure FDA00038971148300000210
In step a)
Figure FDA00038971148300000211
Figure FDA00038971148300000212
B1 can be calculated i =C1 r0 (ii) a Similarly, the slope of the fitted line of red light is
Figure FDA00038971148300000213
Figure FDA00038971148300000214
In step a)
Figure FDA00038971148300000215
B2 can be calculated i =C2 r0
Step c) when I i The value of (a) is 0,
Figure FDA00038971148300000216
the infrared light receiving photon number C1 of the infrared light detector i =B1 i =C1 r0 The number of red light receiving photons C2 of the red light detector i =B2 i =C2 r0 In this case, since the dark current is constant in a short time and the number of photons due to noise is constant, B1 is set i =0, namely the true value of the infrared receiving photon number of the infrared detector is obtained
Figure FDA00038971148300000219
In the same way, let B2 i =0, the true value of the red light receiving photon number of the red light detector is obtained
Figure FDA00038971148300000217
Step three, eliminating the infrared back scattering light, because the infrared light receiving photon number true value of the infrared light detector
Figure FDA00038971148300000218
The number of infrared backward radiation photons is C1S i The number of infrared transmission photons is CT i Then, then
Figure FDA0003897114830000031
And in step c) of
Figure FDA0003897114830000032
Namely CT i =K1 i I i -C1S i Furthermore, the red light receives the photon number true value due to the arrangement of the red light filter
Figure FDA0003897114830000033
Let the number of red backscattered photons be C2S i I.e. by
Figure FDA0003897114830000034
Meanwhile, setting the distance between the convex lens and the red light filter as d, and the backscattering energy equation p (d) of the light detector at the distance d is as follows:
Figure FDA0003897114830000035
in the formula, p 0 Is the emission power of the photodetector, c is the speed of light, τ is the pulse width, A d For the effective receiving area, Y (d) is the optical characteristics of the photodetector, β (d) is the backscattering coefficient, T d Is the transmission factor; according to the transmission factor of light in the aerosol particles
Figure FDA0003897114830000036
Wherein I is the transmitted light intensity, I 0 Is the incident light intensity, this incident light intensity I 0 Output light intensity I of step two i
The conversion of p (d) at this distance d into the number of photons C is in the form of
Figure FDA0003897114830000037
Eta is the quantum efficiency of the optical detector, lambda is the wavelength of light, h is the Planck constant, and t is the time between the emission of the optical transmitter and the reception of the optical detector; the number of red backscattered photons C2S at the distance d can be known i (d) And the number of infrared backscattered photons C1S i (d) In a ratio of
Figure FDA0003897114830000038
Figure FDA0003897114830000039
In the formula of lambda R Wavelength of red light, λ IR Is the wavelength of infrared light, beta R (d) Is the red light backscattering coefficient, beta IR (d) Is the infrared light backscattering coefficient, T rR Is the red light transmission factor, T rIR For the infrared light transmission factor, the infrared backscattered photon number C1S can be obtained by using the formula of the ratio of the red backscattered photon number to the infrared backscattered photon number i By using
Figure FDA00038971148300000310
The infrared transmission photon number CT can be obtained by the formula i
Step four, according to the infrared transmitted photon number CT obtained in step three i The transmission factor of infrared light can be obtained
Figure FDA00038971148300000311
Wherein C is 0 At a light intensity value I 0 The number of received photons of the infrared light detector when the lower light propagation distance is 0; then the Beer-Bouguer-Lambert theorem T rIR =e -σ2d Obtaining extinction coefficient sigma of infrared light, where e is the base of natural number logarithm and 2d is optical path, substituting the obtained extinction coefficient sigma into visibility calculation formula
Figure FDA0003897114830000041
The visibility can be obtained.
CN202110620068.XA 2021-06-03 2021-06-03 Noise removal algorithm of optical sensor and visibility calculation method thereof Active CN113340856B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110620068.XA CN113340856B (en) 2021-06-03 2021-06-03 Noise removal algorithm of optical sensor and visibility calculation method thereof

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110620068.XA CN113340856B (en) 2021-06-03 2021-06-03 Noise removal algorithm of optical sensor and visibility calculation method thereof

Publications (2)

Publication Number Publication Date
CN113340856A CN113340856A (en) 2021-09-03
CN113340856B true CN113340856B (en) 2022-12-06

Family

ID=77473491

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110620068.XA Active CN113340856B (en) 2021-06-03 2021-06-03 Noise removal algorithm of optical sensor and visibility calculation method thereof

Country Status (1)

Country Link
CN (1) CN113340856B (en)

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100561190C (en) * 2006-09-15 2009-11-18 深圳迈瑞生物医疗电子股份有限公司 Dark current detection device and absorbency detection method
CN201503396U (en) * 2009-08-31 2010-06-09 四川科奥达技术有限公司 Atmospheric transmission visibility meter
CN107359174B (en) * 2017-07-11 2023-07-25 展谱光电科技(上海)有限公司 Multispectral image pickup device
CN109297590B (en) * 2018-09-14 2021-01-01 北京遥感设备研究所 Method for calculating background radiation noise of infrared detector
CN114235173B (en) * 2021-11-17 2024-04-09 自然资源部第二海洋研究所 Photon counting satellite-borne ocean laser radar detection simulation method

Also Published As

Publication number Publication date
CN113340856A (en) 2021-09-03

Similar Documents

Publication Publication Date Title
CN108205867B (en) Early fire smoke detection method with interference particle recognition capability
CN109085599B (en) Distance measuring device with SPAD array and distance walk compensation
CN109324318B (en) Distance meter for measuring distance to target object
CN108627812A (en) A kind of laser radar atmospheric visibility measurement method and device
US20190049355A1 (en) Method and Sensor System for Detecting Particles
RU177419U1 (en) Lidar for remote measurement of temperature and humidity with minimal dead zone sounding
Yang et al. The design of rapid turbidity measurement system based on single photon detection techniques
CN107356914B (en) Calibration system for satellite-borne laser radar detector
CN101629917B (en) Method and device for measuring effective atomic number of substance
CN113340856B (en) Noise removal algorithm of optical sensor and visibility calculation method thereof
CN113138398A (en) Aerosol extinction coefficient inversion method based on unmanned aerial vehicle atmospheric laser radar
Kuldashov et al. Infrared sensor for remote monitoring of moisture content in raw cotton
CN110987736B (en) Aerosol particle spectrum and concentration measuring device and method
CN112903547B (en) High-concentration cloud and mist particle concentration measuring device based on double light sources
CN208314186U (en) A kind of laser radar atmospheric visibility measuring device
CN112268871B (en) Method for simultaneously measuring concentration of various polluted gases in atmosphere
RU169314U1 (en) Lidar for remote measurement of temperature and humidity
EP3789793B1 (en) An optical proximity sensor and corresponding method of operation
US20220155442A1 (en) Light detection device, lidar device including the same, and method of measuring distance
CN113805159A (en) Failure detection method, device, equipment and storage medium of signal receiving assembly
CN209656554U (en) Optical thickness spectromonitor for satellite calibration
CN109900655A (en) Optical thickness spectromonitor for satellite calibration
CN209416914U (en) Forward scattering visibility meter linearity detection device
CN101694526A (en) Reflection type optical fiber nuclear radiation sensor
CN104931997B (en) The flat field image measurement method of X-ray energy

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
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Wang Hua

Inventor after: Quan Wei

Inventor after: Tian Jinbao

Inventor after: Wei Zhichao

Inventor after: Man Yongxing

Inventor before: Wang Hua

Inventor before: Quan Wei

Inventor before: Tian Jinbao

Inventor before: Gai Zhichao

Inventor before: Man Yongxing

GR01 Patent grant
GR01 Patent grant