Laser safety with regard to the human eye is a well-known topic. Everybody working with laser sources has to follow the long-established occupational safety rules to prevent people from eye damage by accidental irradiation. These rules comprise, for example, the use of laser safety eyewear and the calculation of the maximum permissible exposure (MPE) and its corresponding hazard distance, the nominal ocular hazard distance. At exposure levels below the MPE, glare effects may occur if the laser wavelengths are in the visible spectral range. The physical effects of laser dazzling on the human eye are described by a quite new concept, which defines the maximum dazzle exposure (MDE) and the corresponding nominal ocular dazzle distance (NODD). Triggered by the MDE/NODD concept, we investigated whether similar laser safety calculations could be performed for electro-optical imaging systems. In this publication, we will review our approach for laser safety calculations for such systems. We have succeeded in finding closed-form equations, allowing calculations of exposure limits to prevent electro-optical imaging systems from damage and/or dazzle. Furthermore, we found some interesting effects related to the corresponding hazard distances, which are also discussed. |
1.IntroductionAnyone who works professionally with lasers knows the long-established occupational safety rules to prevent eye damage by accidental laser irradiation. These include, among others, the calculation of the maximum permissible exposure (MPE) and its corresponding hazard distance, the nominal ocular hazard distance (NOHD)1,2 as well as the choice of proper laser safety eyewear.3,4 Quantities equivalent to MPE and NOHD but related to the reversible effect of laser eye dazzle were recently introduced by Williamson and McLin:5–7 The maximum dazzle exposure (MDE) and nominal ocular dazzle distance (NODD). In the last years, we have attempted to define such quantities for electro-optical imaging systems and to derive equations to calculate those quantities, see Fig. 1. These are:
In this publication, we review our work on laser safety calculations for electro-optical imaging systems, which was presented in detail in several journal articles and conference proceedings.8–12 This includes the derivation of equations to estimate exposure limits (ELs; Sec. 4) and dazzle spot sizes (Sec. 5), approximations to estimate hazard distances (Sec. 6), suggestions for values of nonstandard parameters (Sec. 7), and the validation of this work by modeling and field trials (Sec. 8). We will not go into details here but mainly summarize the results of our work. A brief note on terminology: We use the term electro-optical imaging system for a combination of a camera lens and a camera. The most important component of the camera for further discussion is the image sensor. Sometimes, we denote an electro-optical imaging system simply as a sensor. 2.Fundamentals2.1.ObjectivesWhen we started our work on this topic some years ago, we initially defined the objectives we wanted to reach. The main goal was to develop a tool for laser safety calculations for electro-optical imaging systems that can be applied also by users, who are not experts in this field, even though this may be associated with some limitation of accuracy. The objectives defined were: Objective 1. Equivalent to laser safety calculations for the human eye, the values of and shall be related to the position of the entrance aperture of the camera lens. Objective 2. All equations to calculate the laser safety quantities shall be given as closed-form expressions containing only commonly known operations and functions. Objective 3. All equations to calculate the laser safety quantities should contain, as far as practical, only standard parameters of the involved devices (laser, camera lens, image sensor/camera) and the underlying scenario (e.g., distance, atmospheric extinction). The formulated objectives have the following background: Objective 1 allows the user to position a power meter at the typically easily accessible place in front of the camera lens to compare calculated ELs with the incident laser irradiance. Objective 2 ensures that users who do not have relevant experience in this field can still perform laser safety calculations for sensors. In principle, everybody should be able to perform laser safety calculations using a sheet of paper and a pocket calculator. Therefore, we want to avoid numerical calculations that can only be performed with the help of a computer. Finally, Objective 3 shall enable the user to perform such calculations for a wide range of electro-optical imaging systems, i.e., for different combinations of camera lenses and image sensors, without the necessity to measure unknown parameters beforehand. 2.2.ApproachTo derive damage-related laser safety quantities for electro-optical imaging systems analogous to those of the human eye (EL MPE and hazard distance NOHD), we had to start from the laser-induced damage threshold (LIDT) of the image sensor, which is located at the focal plane of the camera lens. Such damage thresholds for imaging sensors are usually not known and have to be determined by appropriate measurements, e.g., see the work of Becker et al.,13,14 Théberge et al.,15 Burgess et al.,16 Westgate and James,17 and Schwarz et al.18–20 In the next step, we had to transfer the image sensor’s damage threshold to the corresponding value at the position of the camera lens’ entrance aperture to achieve Objective 1. However, this required finding out how the irradiance distribution of the threatening laser beam is related to the irradiance distribution in the focal plane of the camera lens. This issue was the biggest step in our efforts to define the laser safety quantities for sensors since there is not, as with the human eye, only one single kind of camera lens but many different ones. The irradiance distribution at the focal plane depends on the design and quality of the camera lens, comprising its scattering characteristics and image distortions. The scattering properties of camera lenses are usually only very rarely known and must otherwise be determined in time-consuming measurements with a dedicated setup; see, e.g., Ref. 10. The same considerations are also valid for dazzle-related laser safety quantities. For example, the laser dazzle threshold can be defined as the irradiance, where the pixels of the imaging sensor start to saturate. Such saturation thresholds can be calculated easily from the specifications of the image sensor.9 However, to describe the extent of laser dazzle, i.e., the size of the dazzle spot on the image sensor, we again need the quantitative irradiance distribution at the focal plane to compare it with the saturation threshold. In summary: the basis for achieving our goal was to set up a theoretical model (based on closed-form equations to achieve Objective 2) that quantitatively describes the irradiance distribution at the focal plane of a camera lens in the case of laser irradiation. For our application, the main focus of the theoretical model was that it is capable of describing the irradiance distribution in an appropriate manner for our laser safety calculations. There was no intent that the theoretical model describes the irradiance distribution perfectly accurate, so it may also be used for other purposes (e.g., stray light analysis in optics design). Therefore, we concentrated on the estimation of the peak irradiance (which is crucial for laser damage) and the estimation of the large-scale irradiance distribution (which is crucial for the laser dazzle effect). This means that we included diffraction effects in the theoretical model regarding the peak irradiance and scatter/stray light effects for the large-scale irradiance distribution. The importance of including stray light for the calculation of the focal plane irradiance distribution was shown, for example, by Benoist and Schleijpen,21 who modeled the size of laser dazzle spots for charge-coupled device (CCD) cameras. However, we did not include aberrations in our theoretical model since aberrations have typically only an influence on the spatial irradiance distribution near the center of the laser spot at the image sensor as long as the imaging system is focused. The inclusion of aberrations is only required for an accurate description of very small laser dazzle spots in the order of some pixels in diameter, which is rather irrelevant for practical use. Regarding laser dazzle, our interest laid on the description of dazzle spots that affect a considerable amount of the sensor’s field of view (FOV; ). Also, specific optical effects, such as diffraction spikes, have not been considered since they have an influence on the sensor image only in very limited areas and thus have little impact on the overall perceptibility of the scene. 2.3.ScenarioFor our laser safety calculations, we assumed the scenario shown in Fig. 2. A laser emits a Gaussian beam characterized by the output power, ; the laser wavelength, ; the output beam diameter, ; and the full angle divergence, . The laser radiation illuminates an electro-optical imaging system consisting of a camera lens and a camera with an image sensor. The laser beam diameter at the camera lens is denoted by , which refers (in the case of Gaussian beams) to the two opposing positions at the radial irradiance profile, where the irradiance dropped to of the peak irradiance. The camera lens is described by the focal length and the diameter of the entrance pupil . We would like to point out that the location and size of the entrance pupil and the rear principle plane as well as the beam paths in Fig. 2 are drawn for illustrative purposes only. Further parameters of the camera lens are the -number , the number of optical elements , the transmittance , and three scatter parameters , (or alternatively ), and , which will be described in more detail in Sec. 3.1. The ratio of the beam diameter to the diameter of the camera lens’ entrance pupil is called the truncation factor and, as we will see later, has a determining influence on the distribution of the laser light in the focal plane of the camera lens where the laser beam diameter in a distance to the laser source can be calculated byUsing the truncation factor , the fraction of the laser power that enters the lens can be calculated by The image sensor is characterized by several parameters: the number of pixel columns and rows , , pixel size , total quantum efficiency , saturation capacity , and the integration time . We assume that the sensor system is well focused, i.e., the imaging sensor is placed at (or very near to) the focal plane of the camera lens. The onset of damage and dazzle is described by irradiance values (LIDT) and (saturation irradiance). The attenuation of laser power by the atmosphere may be included by substituting with , where is the atmospheric extinction coefficient. 2.4.ParametersTable 1 summarizes all parameters that we use for our laser safety calculations. Table 1Parameters used for the laser safety calculations.
According to Objective 3, our equations shall be based only on standard parameters that are typically specified by the manufacturers of lasers, camera lenses, or image sensors. However, besides such standard parameters, some of the aforementioned parameters are usually not specified or known:
In Sec. 7, we state values/equations for these nonstandard parameters that may be used, if measured values for a specific electro-optical imaging system are not available. 3.Estimation of the Focal Plane Irradiance DistributionIn our theoretical model, the incident power contributes to the focal plane irradiance distribution by two components: (a) the scatter/stray light component and (b) the diffraction component : where describes the fraction of the incident laser power that is diffracted and is the radial coordinate. In Eq. (4), as well as in all subsequent equations, the dependency on the radial coordinate (in the focal plane) can be replaced by the dependency on the viewing angle using the relationship3.1.Stray Light ComponentTo estimate the contribution of the stray light to the focal plane irradiance distribution , we relied on the work of Peterson,22 who published an analytical approach for this task using the three-parameter Harvey scatter model as a bidirectional scattering distribution function (BSDF). This type of BSDF describes the angular distribution of stray light caused by scatter from smooth surfaces of optical elements using three parameters, , , and . Other sources of stray light such as multiple reflections or scatter from the lens housing are not considered. For a detailed explanation of the BSDF and the meaning of the scatter parameters, we refer the reader to other publications, e.g., Ref. 23. For our theoretical model, we applied some simplifications (see Ref. 9 for more details) to the work of Peterson22 to keep equations manageable for typical camera lenses with five or more optical elements. As a result, we calculate the stray light component by where is defined byis the number of scattering surfaces of the camera lens, which we assume is twice the number of optical elements . The modified truncation factor considers that the beam diameter within the camera lens cannot increase, if a laser beam, which is already much larger than the entrance aperture, would expand further, e.g., with increasing distance to the laser source. The ratio of scattered power to the incident power (for a single scattering surface) is called the total integrated scatter (TIS) and can be calculated by23 Please note: Eq. (8) uses an alternative scatter parameter that has its origins in the two-parameter Harvey scatter model (see, e.g., Ref. 21). This scatter parameter is linked to scatter parameter by the following equation: The quantity is used to calculate the fraction of the incident power that is diffracted; see Eq. (4) The scatter parameters are wavelength dependent; the wavelength dependency is described by the equation24 Reference 24 stated that these wavelength scaling laws are valid only “over a limited range of ,” but limits are not given. Wein,25 for example, stated that the “wavelength scaling approximately predicts the scattering in the visible but not in the far infrared.” 3.2.Diffraction ComponentFor the diffraction component , we assumed in our theoretical model a Gaussian beam and used for our calculations the work of Urey26 regarding the diffraction pattern of a truncated Gaussian beam. This diffraction pattern can be imagined to look like something between the Airy diffraction pattern and a pure Gaussian distribution. The shape depends on the value of the truncation factor and consists of a central lobe that can be approximated by a Gaussian distribution and diffraction rings of lower power similar to the Airy diffraction pattern. The central lobe is approximated by where the peak irradiance of the diffraction pattern is given by and the spatial extent of the central lobe is calculated by26Here, is a spot size constant, which also depends on the truncation factor , see Ref. 26 for details For a typical camera lens designed for the visible spectral range, the diffraction spot size is usually in the order of some micrometers. Since the pixel size of common CCD or complementary metal-oxide-semiconductor (CMOS) image sensors is of the same order, the central lobe is usually not resolved by an image sensor. Thus, for our laser safety calculations, we relied on the peak irradiance only. Outside the central lobe, the wings of the diffraction pattern are described by the mean of the diffraction ring irradiance, which is given by Since Eq. (16) describes the mean of the diffraction ring irradiance, there is no oscillating term in the formula. Here, the diffraction ring irradiance is a monotonically decreasing function with a dependence. 3.3.Focal Plane Irradiance DistributionSome examples of radial irradiance profiles as calculated with Eq. (4) are plotted in Fig. 3. The three graphs show the normalized irradiance as a function of the radial coordinate for different values of the truncation factor . For the calculations, the parameters listed in Table 2 were used. Table 2Parameters used for the example calculations of Figs. 3 and 5.
4.Exposure LimitsAfter setting up the theoretical model for the focal plane irradiance distribution, the EL and could be derived. For the derivation of these two quantities, we did not employ the exact equation for but used only the respective essential components. This means that we used the peak irradiance only for the derivation of the since the onset of laser damage depends on the peak irradiance. By contrast, the derivation of the is based on the irradiance of the diffraction rings and the stray light irradiance since the spatial extent of the dazzle spot essentially depends only on these components; see Fig. 3. Only for the treatment of tiny dazzle spots, the irradiance of the central lobe would be needed. As mentioned before, we considered this case as not relevant for practical use and focused on larger dazzle spots that comprise considerable amounts of the sensors FOV (>10%). 4.1.MPE for a SensorThe can be calculated by equalizing the focal plane peak irradiance of Eq. (13) and the image sensor’s damage threshold . This leads to Resolving for the maximum permissible laser power and calculating the corresponding laser peak irradiance by dividing by result in the searched quantity From Eq. (18), we see that the depends on the truncation factor, which means that the also depends on the distance of the laser to the imaging system. This has a crucial impact on the calculation of the corresponding hazard distance, as we will explain in more detail in Sec. 6. The lowest value of the (worst case) occurs for Equation (19) is applicable when the laser source is far away from the electro-optical imaging system and the laser beam overfills the aperture of the camera lens. In this case, no distance dependence exists. 4.2.MDE for a SensorFor the human eye, the MDE is not just a single value but is stated for specific dazzle levels.6,7 The dazzle levels range from very low to low, medium, and high, which corresponds to angular dazzle fields of 2, 10, 20, and 40 deg, respectively. For electro-optical imaging systems, an equivalent definition of such default values for the dazzle field is not useful since the system’s FOV changes with the focal length of the camera lens. For an imaging system, we defined the dazzle level as the fraction of the system’s FOV that is dazzled. This means, e.g., that for an incident irradiance of , and , a tenth of the FOV, half of the FOV, and the full FOV are dazzled, respectively. Figure 4 illustrates that approach. The fraction shall be understood as the diameter of the dazzle spot divided by the size of the longer side of the image sensor. Using this definition, the angular radius of the dazzle spot is where the FOV of the imaging system can be calculated byWe can find the for a specific dazzle level by equating the focal plane irradiance distribution and the image sensor’s saturation threshold . As previously addressed, we approximate the focal plane irradiance distribution for this case through the diffraction ring irradiance and the stray light component . Thus, we obtain the following equation: Equivalent to the , the searched quantity can be derived by resolving for the maximum permissible laser power and calculating the corresponding laser irradiance by dividing by As in the case of the , the equation simplifies for the case of , which gives the minimum value of the Equations (24) and (25) should not be used for vanishing values of , such as or very close to zero, see Ref. 9 for more details. The value would correspond to the onset of laser dazzle. This onset can be estimated using Eq. (18) or (19) but replacing the focal plane damage threshold by the saturation threshold 5.Dazzle Spot SizeThe size of a dazzle spot can be calculated by solving Eq. (22) for Unfortunately, there is no closed-form expression for this solution. This equation can be solved numerically using a computer. However, an approximate analytic solution can be given by solving the equations for the diffracted irradiance and the scattered irradiance individually and choosing the maximum value as the dazzle radius. The equation to solve for the irradiance of the diffraction rings is and the equation to solve for the scattered irradiance isSolving Eqs. (29) and (30) for and, respectively, we obtain and finally6.Hazard DistancesUnfortunately, the calculation of hazard distances for electro-optical imaging systems is different compared with those for the human eye. As we can see from Eqs. (18) and (24), both the and the depend on the truncation factor . This implicates a dependence of these ELs on the distance between the laser source and imaging system since the beam diameter changes with distance and thus the truncation factor does. The conventional equations to calculate the human eye-related NOHD and NODD (neglecting atmospheric extinction) are Please note that the beam diameter and divergence in Eqs. (34) and (35) have to refer to the opposite positions of the laser beam profile at which the irradiance has dropped to of the peak irradiance. This is a common practice in laser safety, whereas, in this publication, we typically refer to the points unless otherwise stated. Nevertheless, these equations cannot be applied in our case, even if the correct quantities are used since Eqs. (34) and (35) were derived for distance-independent ELs. The distance-dependent values of and are plotted in Fig. 5 for a generic electro-optical imaging system and laser source; the parameters assumed for this calculation are listed in Table 2. For the graph in Fig. 5(a), the independent parameter is the distance of the electro-optical imaging system to the laser source, whereas, for the graph in Fig. 5(b), we used the truncation factor as the independent parameter. While the is a single curve (plotted in red), the is represented by a light blue band, which indicates the range of values for various dazzle levels . The curve for is highlighted by a blue line; the upper and lower border of the blue band correspond to and , respectively. In addition, the laser peak irradiance at the front face of the camera lens is shown as a green curve. Please note that the parameter in Table 2 is related to the focal plane, whereas the ordinate of the plot is related to the front face of the camera lens. In Fig. 5, we can recognize how the and vary with distance. Following the curves starting from large distance values, the ELs are quite constant, which correspond to the minimum values as given by Eqs. (19) and (25). For closer distances from to , the ELs increase strongly (for the given example) with decreasing distance. Finally, for distances below , the ELs stay constant again. A very important result concerns the hazard distance regarding the curve: Looking at the green curve, which indicates the laser peak irradiance, we can see that this curve intersects the curve twice. As long as the irradiance curve is below the curve, the imaging system is safe from damage. In the distance range where the irradiance curve is above the curve, the imaging system is not safe from damage. By the example of Fig. 5, we can recognize two intersection points , which cut the curve into three sections (see vertical dashed red lines): In the first section, the image sensor is safe from damage within the distance from the laser to the first intersection point. Between the two intersection points, the image sensor is not safe from damage, and in section three, beyond the second intersection point, the image sensor is again safe from damage. This means that the laser hazard zone for imaging systems may have a distance dependent upper as well as a lower limit. This is in contrast to laser hazard distances valid for the human eye. We acknowledge that we intentionally chose the laser parameters in such a way that this effect is eye-catching in the example of Fig. 5. If we choose a sufficiently higher value of laser power , the irradiance curve would be shifted so much upward that the irradiance curve would intersect the curve only once. In our example, this is the case for the curves in Fig. 5. For each dazzle level (0.1, 0.5, or 1.0), the irradiance curve will intersect a specific curve only once. For another choice of system parameters (mainly the scatter parameter ) that may be different. 6.1.Derivation of Hazard DistancesThe hazard distances NSeHD and NSeDD are defined by the fact that the laser irradiance at the system’s camera lens is equal to the and , respectively. Using EL as a synonym for either the or the , this means that with where is the truncation factor in the hazard distance . Resolving Eq. (36) for would result in the hazard distance NSeHD/NSeDD. Unfortunately, there is no solution in closed form, when we apply Eqs. (18) and (24) for the and the in Eq. (36), respectively.We recognized that the truncation factor is the determining parameter for the and . Thus, we chose in the first step the approach not to determine but the corresponding truncation factor . This led to somewhat easier equations for the derivation of the NSeHD and the NSeDD, although these equations were still not resolvable analytically but only numerically. Then, in a second step, the terms of the equations to be solved, which depended on the truncation factor , were substituted by approximations that allowed solutions for in closed form. We refer the reader to Ref. 8 for detailed explanations and will present here the results only. Finally, having an approximation for , the corresponding hazard distance (NSeHD or NSeDD) can then be calculated by 6.2.Nominal Sensor Hazard DistanceFor the NSeHD, we found approximate estimates both for the upper and lower values of 8 which can be used to calculate the upper and lower values of the laser hazard zone byThese approximations are well usable to determine hazard distances for sensor damage. In Ref. 8, approximate values of NSeHD calculated using Eqs. (39) and (40) () were compared with values computed by numerically solving Eq. (36) (). The calculations were performed for a total of 5734 different parameter sets. The relative error of the approximate NSeHD values was between [0, +0.5] in more than 90% of cases. The results are graphically displayed in Fig. 6: Fig. 6(a) as a histogram of the values (bin width 0.5) and Fig. 6(b) a plot of the relative error as a function of the numerically calculated value of the truncation factor . In both graphs, red and blue color correspond to the upper and lower values of the hazard distance, respectively.We found an overestimation of values rather than an underestimation, which is good in terms of laser safety. However, it should be noted that these investigations were based on a limited database and the results may not be generally valid. 6.3.Nominal Sensor Dazzle DistanceThe search for a closed-form expression to calculate the NSeDD was quite demanding. This is due to the fact that the equation that has to be solved to calculate the NSeDD is quite complex. As for the NSeHD, there may be an upper and lower limit of the laser hazard zone. Furthermore, for a specific dazzle level , there may be areas within the laser hazard zone where the dazzle level is lower than that specific value: . This means that the size of the dazzle spot in the camera image may decrease within the laser hazard zone (for the specific dazzle level ) as the imaging system approaches the laser source. For more details, see Ref. 8. Due to this complexity, we decided to state only an approximation for the upper limit of the laser hazard zone in case of laser dazzle. Even for this limited case, we could not find a generally applicable approximation in a closed form. Nevertheless, we can give approximations for the NSeDD for two special cases: (a) For extended laser beams, i.e., the laser beam diameter is much larger than the aperture of the camera lens (). (b) For beam diameters that are smaller than the lens aperture () in combination with larger dazzle levels. However, to anticipate, the accuracy of these approximations is not as good as in the case of the NSeHD. The NSeDD in the case of extended laser beams () can be approximately calculated by For laser beams with a diameter smaller than the aperture of the camera lens (), the NSeDD can be approximated by For Eq. (43), the following constraints have to be fulfilled: Both constraints imply that the dazzle level should be sufficiently large, particularly that the influence of the scatter component outweighs that of the diffraction component regarding the focal plane irradiance distribution. Finally, the NSeDD can be calculated using Eq. (42) or (43) by Equivalently to the NSeHD, we also investigated the accuracy of these approximations. The NSeDD was first estimated using the approximation of Eqs. (42) and (43). Since both equations give an approximation for the upper limit of the laser hazard zone, we chose the maximum of both values for further processing (). Second, the NSeDD was computed by numerically solving Eq. (36) (). Subsequently, the relative error of was calculated as before for the NSeHD These calculations were performed for 2,292,924 different parameter sets. In Fig. 7, the results are plotted for the NSeDD. In these graphs, the red color corresponds to the approximation for calculated using Eq. (42). The blue color corresponds to the approximation for calculated using Eq. (43). In the case of the NSeDD, we see a similar distribution in the histogram but not as good as for the NSeHD. Larger relative errors occur for the approximation for . We can recognize in the graph on the right-hand side that these cases mainly occur for values of the truncation factor but not exclusively. Furthermore, we can see that the relative error is negative for a considerable amount of the calculations. A negative relative error corresponds to an underestimation of the hazard distance. We conclude that Eqs. (42) and (43) have to be used with care for the estimation of hazard distances for sensor dazzle. 7.Nonstandard Parameters for Laser Safety Calculations7.1.Damage Thresholds of Image SensorsInformation on LIDTs of image sensors for continuous wave (cw) laser radiation in the visible spectral range is rare. Publications related to measured cw LIDT of CCD and CMOS cameras are, for example, Becker et al.,13,14 Théberge et al.,15 Burgess et al.,16 Westgate and James17 and Schwarz et al.18–20 Here, we refer to a publication by Schwarz et al.18 and summarize these threshold values in Table 3. Please note that Schwarz et al. measured these values for specific image sensors (CCD sensor: Sony ICX098, CMOS sensor: Aptina MT9V024). Laser damage thresholds for other types of image sensors may vary, but the order of magnitude () should be similar. Table 3One-on-one LIDT measured for specific image sensors (CCD: Sony ICX098, CMOS: Aptina MT9V024). Data reproduced with permission from Ref. 18.
7.2.Saturation Thresholds of Image SensorsThe saturation threshold of an image sensor may be calculated using its technical specifications. At a specific irradiance (i.e., the saturation irradiance), the number of photoelectrons generated by the incident photons within the camera’s exposure time will equal the saturation capacity of a pixel. Using this relation, the saturation irradiance of a pixel of size can then be estimated by Putting just the pixel saturation irradiance from Eq. (48) into the equations [Eqs. (24)–(27)] would imply that the image sensor is illuminated by the dazzle laser only. In a real situation, the imaging system typically observes a scene, which means that the charge generated by a pixel is determined by the ambient light of the scene and the laser light. To meet this condition, we assume that the operator or the camera’s automatic exposure control will set the exposure time to a level, such that the image sensor’s mean pixel signal equals roughly half of the maximum pixel signal. Thus, applying a factor of 0.5 to Eq. (48), we obtain a first estimate for the saturation irradiance If the saturation level of a camera image is known in more detail, e.g., when using a specific test chart in a laboratory environment, the saturation threshold may be estimated more precisely by where a saturation level means not saturated and means completely saturated.7.3.Scatter ParametersIn our simple theoretical model in Sec. 3, the scatter component contributing to the focal plane irradiance distribution depends on the scatter parameters and of the camera lens used. Unfortunately, these scatter parameters are usually not known for commercial off-the-shelf (COTS) camera lenses. Thus, we performed a series of measurements where we recorded images of the focal plane irradiance distribution for a selection of seven typical camera lenses with focal lengths ranging from 25 to 100 mm. We then derived radial irradiance profiles from the image data and subsequently fitted our Eq. (4) to these irradiance profiles, where we used the quantities and as fit parameters. The outcome of this work was an individual set of scatter parameters for each camera lens under test; see Table 4. Since these sets describe the corresponding camera lens as a whole (and not the scattering properties of the single scattering surfaces), we denote them as integrated scatter parameters. Details can be found in Ref. 10. Table 4Integrated scatter parameters for various commercial off-the-shelf camera lenses. Data reproduced with permission from Ref. 10.
We noticed that the integrated scatter parameters for the different camera lenses are of the same order of magnitude. Therefore, we additionally derived a generic set of integrated scatter parameters based on a statistical analysis of the individual sets. This generic set may be applied together with our equations to predict the laser dazzle of an electro-optical imaging system when a camera lens with unknown scatter parameters is used: 8.ValidationThe approximate equations for and NSeHD/NSeDD presented above only apply under the condition that the underlying theory of Sec. 3 is appropriate. To some extent, the equations were validated by our work published in Ref. 10, where we measured the integrated scatter parameters for various camera lenses; see Sec. 7.3. In this publication, we showed that the radial irradiance profile caused by incident laser irradiation can be described by our theoretical model. Nevertheless, it was our endeavor to test their usability following different approaches. The first test (Sec. 8.1) was focused on the usability of the integrated scatter parameters as measured according to Ref. 10. Our theoretical model contains a couple of simplifications and assumptions, which may cast doubt on its suitability. Thus, we compared calculations of the focal plane irradiance distribution based on our theoretical model with the outcome of the optics design software FRED;27 for a two-element lens and a lens of Double-Gauss type. The second test (Sec. 8.2) was to compare the size of laser dazzle spots in camera images acquired during a field trial with the theoretical predictions of the dazzle spot size according to Eq. (33). 8.1.Validation Using the Optical Design Software FREDThe focal plane irradiance distribution was simulated for two different optical systems using the FRED optical design software from Photon Engineering: The advantage of using the two-element laser focusing optics was its availability and the disclosed optical design. This allowed us to measure (a) the integrated scatter parameters of the optics as well as (b) its modeling using the FRED software. We then compared the FRED simulation results (based on the measured integrated scatter parameters) with the outcome of our theoretical model to show that the theoretical model is reliable. The FRED simulation results shall also serve to validate that our theoretical model, in conjunction with our generic set of integrated scatter parameters (see Sec. 7.3), is able to describe even quite complex multielement camera lenses as well. Hence, we simulated the irradiance distribution at the focal plane of the generic Double-Gauss camera lens using our generic set of scatter parameters; see Eq. (51). The design data of the Double-Gauss camera lens were taken from a standard textbook on optical design.28 The computation of the focal plane irradiance distribution using the FRED software was implemented by means of a multistep simulation. One simulation was focused on the diffraction component and the second one on the stray light component. The output of both simulation runs was then combined and subsequently compared with calculations based on Eq. (4); for details, see Ref. 11. Some results are presented in Fig. 8. The scatter plots of Fig. 8 show the radial irradiance distribution as computed by the FRED software for the aforementioned optical systems and different values of the -number . Each data point corresponds to a simulated sensor pixel. The green data points correspond to the results of the simulation run focused on the scatter component, while the blue data points represent the result of the simulation run focused on the diffraction component. Their combination is plotted by red data points. Furthermore, the output of the theoretical model is plotted as a black solid curve. For the LINOS 033486 laser focusing lens [Figs. 8(a) and 8(b)], we can see that the theoretical model predicts quite well the envelope of the FRED-simulation data. Of course, obtaining such a good result for this focusing lens is not really surprising since it has a rather simple optical design, and the integrated scatter parameters have been experimentally determined and served as an input for the FRED simulation. Regarding the Double-Gauss camera lens [Figs. 8(c) and 8(d)], one may recognize for radial coordinates that there is a larger deviation between the results of the model and the simulations, which we attribute to aberrations that are not covered by our theoretical model. For radial coordinates larger than , we can see from the scatter plots that the theoretical model fits reasonably to the simulation results. For small -numbers, see Fig. 8(c), there seems to be an overestimation of the irradiance, while for large -numbers, see Fig. 8(d), an underestimation occurs. Nevertheless, we got an adequate agreement between our simplified theoretical model and the output of the corresponding FRED simulation. 8.2.Validation Using Data from Outdoor MeasurementsDuring a field trial, an electro-optical imaging system at a distance of 660 m was dazzled on a slant path by cw laser radiation.29 Camera images were acquired for various laser wavelengths and camera exposure times. In the subsequent data analysis, dazzle spot sizes were assessed using overexposed pixel counting (OPC) and then compared to the theoretical predictions of the dazzle spot size according to Eq. (33). As dazzler, we used the Toptica iChrome MLE multiwavelength laser device. This laser source comprises four different lasers (wavelengths 488, 515, 561, and 640 nm) with output powers ranging from 40 to 100 mW. The data of the laser device are listed in Table 5. The actually used laser output power was measured during the field trial, whereas the output diameter and full angle divergence were measured in the laboratory using an measurement system (Thorlabs M2MS-BC106VIS/M). The laser beam diameter at the imaging system was calculated from these data using Eq. (2). For the atmospheric attenuation of the laser radiation, we assumed an extinction coefficient of since the weather conditions were good. Table 5Device parameters used for the evaluation of the measurement data gained during the field trial.
The imaging system consisted of a Schneider-Kreuznach Xenoplan 2.8/50 COTS camera lens30 mounted to an Allied Vision Mako G-223B NIR camera.31 The specifications of both devices are also listed in Table 5. The transmittance of the camera lens was not known. Thus, we assumed a broadband antireflection coating with a reflection below per surface, resulting in . Furthermore, we intentionally did not measure the scatter parameters of this camera lens but simply used the generic set of integrated scatter parameters as stated in Eq. (51) to investigate their applicability. The data of the camera were taken from the datasheet, see Ref. 31. From the corresponding graph of the data sheet, we estimated a quantum efficiency of for the laser wavelengths used. In summary, our set of measured, specified, and estimated values represents a quite suitable mix for testing the usability of our approach. In Fig. 9, some examples of camera images are shown for the laser wavelengths of 515 and 640 nm. The images were taken using three different exposure times of , 1 ms, and 10 ms. The data was processed according to the following way:
The three calculated dazzle spot sizes are drawn into the camera images of Fig. 9 as colored circles. Red circles indicate the dazzle spot size as calculated by the OPC method. The dazzle spot sizes according to our theoretical model are drawn as green and blue circles for a saturation level of 0.5 and the image-specific saturation level, respectively. Looking at Fig. 9, we can see that our approach for laser safety calculations can indeed give a reasonable estimate for the dazzle spot sizes. In the case of the underexposed images ( and ), the use of ISSL leads to a very good agreement between the OPC result and the calculated dazzle spot size using the theoretical model. Using the saturation level slightly overestimates the dazzle effect but not to an extraordinary degree. In the case of the slightly overexposed images (), the theoretical prediction of the dazzle spot size for saturation level 0.5 fits quite well to the dazzled area of the camera image, whereas there is an overestimation using the image-specific saturation level. At first glance, one might say that the calculated dazzle spot size using the ISSL fits quite well with the dazzle spot size computed using the OPC method. However, one has to keep in mind that the camera image was slightly overexposed, which means that the OPC method also counts pixels overexposed by ambient light and not only by laser light. Thus, the dazzle spot size computed with the OPC method overestimates the laser effect. The results for all laser wavelengths used at the field trial and for all values of camera exposure time are presented in the graphs of Fig. 10. For each laser wavelength, there is a graph that contains the results of the calculations according to the theoretical model and the OPC evaluation applied to the image data. The theoretical results comprise both the use of saturation level and the ISSL. The graphs also confirm the above statements for the other laser wavelengths of 488 and 561 nm. In summary, there is a good agreement between theory and experiment, considering that not all parameters involved were exactly specified. Since the underlying theory for the prediction of dazzle spot sizes is the same as for the estimation of the ELs and , we assume that the equations for estimating the ELs can also be considered validated. 9.ConclusionsIn this publication, we presented an approach for laser safety calculations for electro-optical imaging systems. It includes a set of equations to estimate ELs and hazard distances regarding laser-induced damage and dazzle. Furthermore, dazzle spot sizes can be predicted. All equations are in closed form, which allows their use, in principle, with a pocket calculator, a pen, and a sheet of paper. A computer is not necessary to numerically solve equations, even though we have to admit that some of the equations are somewhat longer and a computer may still be helpful. Furthermore, we tried to include in our equations only standard parameters of the involved systems (laser, camera lens, image sensor) that are typically specified by the manufacturer. However, some nonstandard parameters have to be known, which are: the LIDT and saturation threshold of the image sensor as well as the (integrated) scatter parameters of the camera lens. Values for LIDTs of image sensors can be found in the literature; saturation thresholds can be calculated from the image sensor’s specifications. Regarding scatter parameters for commercial off-the-shelf camera lenses, we performed a series of experiments to measure them. From these measurements, we derived a generic set of scatter parameters that may be used if the exact values for specific camera lenses are not available. To validate our theory, i.e., the derived closed-form equations as well as the generic set of scatter parameters, we performed several experimental tests, field trials, and optics modeling. In our opinion, our studies indicate the suitability of our equations for the desired application: laser safety calculations for electro-optical imaging systems. DisclosuresThe authors have no relevant financial interests in the study and no other potential conflicts of interest. Code and Data AvailabilityThe data that support the findings of this article are not publicly available due to classification. They can be requested from the author at gunnar.ritt@iosb.fraunhofer.de. ReferencesAmerican National Standard,
“Safe use of lasers outdoors,”
(2015). Google Scholar
Bundesanstalt für Arbeitsschutz und Arbeitsmedizin (BAuA),
“TROS Laserstrahlung Teil 2: Messungen und Berechnungen von Expositionen gegenüber Laserstrahlung,”
(2018). Google Scholar
American National Standard,
“Safe use of lasers,”
(2022). Google Scholar
European Standard,
“Personal eye-protection equipment—filters and eye-protectors against laser radiation (laser eye-protectors),”
(2017). Google Scholar
C. A. Williamson and L. N. McLin,
“Nominal ocular dazzle distance (NODD),”
Appl. Opt., 54 1564
–1572 https://doi.org/10.1364/AO.54.001564 APOPAI 0003-6935
(2015).
Google Scholar
C. A. Williamson and L. N. McLin,
“Determination of a laser eye dazzle safety framework,”
J. Laser Appl., 30 032010 https://doi.org/10.2351/1.5029384 JLAPEN 1042-346X
(2018).
Google Scholar
C. A. Williamson and L. N. McLin,
“Laser eye dazzle safety framework,”
in Int. Laser Saf. Conf.,
164
–165
(2017). https://doi.org/10.2351/1.5056871 Google Scholar
G. Ritt,
“Laser safety—what is the laser hazard distance for an electro-optical imaging system?,”
Sensors, 23 7033 https://doi.org/10.3390/s23167033 SNSRES 0746-9462
(2023).
Google Scholar
G. Ritt,
“Laser safety calculations for imaging sensors,”
Sensors, 19 3765 https://doi.org/10.3390/s19173765 SNSRES 0746-9462
(2019).
Google Scholar
G. Ritt, B. Schwarz and B. Eberle,
“Estimation of lens stray light with regard to the incapacitation of imaging sensors,”
Sensors, 20 6308 https://doi.org/10.3390/s20216308 SNSRES 0746-9462
(2020).
Google Scholar
B. Schwarz, G. Ritt and B. Eberle,
“Estimation of lens stray light with regard to the incapacitation of imaging sensors—part 2: validation,”
Sensors, 22 9447 https://doi.org/10.3390/s22239447 SNSRES 0746-9462
(2022).
Google Scholar
G. Ritt et al.,
“An approach for laser safety calculations for electro-optical imaging systems,”
Proc. SPIE, 12738 127380F https://doi.org/10.1117/12.2680178 PSISDG 0277-786X
(2023).
Google Scholar
M. F. Becker et al.,
“Laser-induced damage to silicon CCD imaging sensors,”
Proc. SPIE, 1105 68
–77 https://doi.org/10.1117/12.960613 PSISDG 0277-786X
(1989).
Google Scholar
M. F. Becker et al.,
“Laser-induced functional damage to silicon CCD sensor arrays,”
Proc. SPIE, 1624 67
–79 https://doi.org/10.1117/12.60092 PSISDG 0277-786X
(1992).
Google Scholar
F. Théberge et al.,
“Damage thresholds of silicon-based cameras for in-band and out-of-band laser expositions,”
Appl. Opt., 61 2473
–2482 https://doi.org/10.1364/AO.450317 APOPAI 0003-6935
(2022).
Google Scholar
C. Burgess et al.,
“Modelled and experimental laser-induced sensor damage thresholds to continuous wave infrared sources,”
Proc. SPIE, 10797 107970T https://doi.org/10.1117/12.2326401 PSISDG 0277-786X
(2018).
Google Scholar
C. Westgate and D. James,
“Visible-band nanosecond pulsed laser damage thresholds of silicon 2D imaging arrays,”
Sensors, 22 2526 https://doi.org/10.3390/s22072526 SNSRES 0746-9462
(2022).
Google Scholar
B. Schwarz et al.,
“Laser-induced damage threshold of camera sensors and micro-optoelectromechanical systems,”
Opt. Eng., 56 034108 https://doi.org/10.1117/1.OE.56.3.034108
(2017).
Google Scholar
B. Schwarz et al.,
“Further investigation on laser-induced damage thresholds of camera sensors and micro-optomechanical systems,”
Proc. SPIE, 11161 111610A https://doi.org/10.1117/12.2533142 PSISDG 0277-786X
(2019).
Google Scholar
B. Schwarz, G. Ritt and B. Eberle,
“Impact of threshold assessment methods in laser-induced damage measurements using the examples of CCD, CMOS, and DMD,”
Appl. Opt., 60 F39
–F49 https://doi.org/10.1364/AO.423791 APOPAI 0003-6935
(2021).
Google Scholar
K. W. Benoist and R. (H.) M. A. Schleijpen,
“Modelling of the over-exposed pixel area of CCD cameras caused by laser dazzling,”
Proc. SPIE, 9251 92510H https://doi.org/10.1117/12.2066305 PSISDG 0277-786X
(2014).
Google Scholar
G. L. Peterson,
“Analytic expressions for in-field scattered light distributions,”
Proc. SPIE, 5178 184
–193 https://doi.org/10.1117/12.509120 PSISDG 0277-786X
(2004).
Google Scholar
E. C. Fest, Stray Light Analysis and Control, SPIE Press, Bellingham, Washington, USA
(2013). Google Scholar
, “Scattering in ASAP,”
http://www.breault.com/knowledge-base/scattering-asap
().
Google Scholar
S. J. Wein,
“Small-angle scatter measurement,”
University of Arizona, Tucson, AZ, USA,
(1989).
Google Scholar
H. Urey,
“Spot size, depth-of-focus, and diffraction ring intensity formulas for truncated Gaussian beams,”
Appl. Opt., 43 620
–625 https://doi.org/10.1364/AO.43.000620 APOPAI 0003-6935
(2004).
Google Scholar
M. Laikin, Lens Design, 4th ed.CRC Press, Boca Raton, FL, USA
(2007). Google Scholar
B. Eberle et al.,
“NATO SET-249 joint measurement campaign on laser dazzle effects in airborne scenarios,”
Proc. SPIE, 11161 111610C https://doi.org/10.1117/12.2533744 PSISDG 0277-786X
(2019).
Google Scholar
, “Datasheet Schneider-Kreuznach Tourmaline 2.8/50 (formerly Schneider-Kreuznach Xenoplan 2.8/50,”
https://schneiderkreuznach.com/application/files/6216/4241/8065/TOURMALINE_28_50_C_1001976_datasheet.pdf
().
Google Scholar
, “Datasheet allied vision Mako G-223B NIR,”
https://cdn.alliedvision.com/fileadmin/pdf/de/Mako_G-223B_NIR_DataSheet_de.pdf
().
Google Scholar
BiographyGunnar Ritt is a research associate at Fraunhofer IOSB, Ettlingen, Germany. He received his diploma and PhD degrees in physics from the University of Tübingen, Germany, in 1999 and 2007, respectively. His main research focus is on laser protection. Bastian Schwarz has been a research scientist at Fraunhofer IOSB, Ettlingen, Germany, since 2013 and has been the head of the optical countermeasure group since 2022. He graduated in physics from the University of Freiburg in 2012. Since 2013, his research areas include laser protection and laser damage performance. Bernd Eberle has been a senior scientist at Fraunhofer IOSB in Ettlingen, Germany, where he headed the optical countermeasure and laser protection group. He received his diploma and PhD degrees in physics from the University of Konstanz in 1983 and 1987, respectively. His research activities included laser technology, laser spectroscopy, nonlinear optics, femtosecond optics, optical countermeasures, including protection against laser radiation, and imaging laser sensors. Michael Henrichsen is a research associate at Fraunhofer IOSB, Ettlingen, Germany. He has been the head of the optronic analysis systems and laser effects group since 2022. He received his master of science degree in physics from the University of Constance, Germany, in 2012. His main research focus is on laser effects and infrared imaging. |