EP1673613A1 - Infrarotnachweisvorrichtung - Google Patents
InfrarotnachweisvorrichtungInfo
- Publication number
- EP1673613A1 EP1673613A1 EP04761373A EP04761373A EP1673613A1 EP 1673613 A1 EP1673613 A1 EP 1673613A1 EP 04761373 A EP04761373 A EP 04761373A EP 04761373 A EP04761373 A EP 04761373A EP 1673613 A1 EP1673613 A1 EP 1673613A1
- Authority
- EP
- European Patent Office
- Prior art keywords
- calibration
- infrared
- temperature difference
- field
- view
- 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.)
- Withdrawn
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 25
- 230000005855 radiation Effects 0.000 claims abstract description 43
- 230000002411 adverse Effects 0.000 claims abstract description 31
- RAHZWNYVWXNFOC-UHFFFAOYSA-N Sulphur dioxide Chemical compound O=S=O RAHZWNYVWXNFOC-UHFFFAOYSA-N 0.000 claims description 40
- 238000000034 method Methods 0.000 claims description 33
- 239000000428 dust Substances 0.000 claims description 26
- 238000005259 measurement Methods 0.000 claims description 21
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 20
- 238000012544 monitoring process Methods 0.000 claims description 13
- 230000000694 effects Effects 0.000 claims description 12
- 238000010521 absorption reaction Methods 0.000 claims description 9
- 235000010269 sulphur dioxide Nutrition 0.000 claims description 9
- 239000004291 sulphur dioxide Substances 0.000 claims description 9
- 238000012937 correction Methods 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000012546 transfer Methods 0.000 claims description 2
- 239000002956 ash Substances 0.000 description 78
- 238000009826 distribution Methods 0.000 description 11
- CCEKAJIANROZEO-UHFFFAOYSA-N sulfluramid Chemical group CCNS(=O)(=O)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)F CCEKAJIANROZEO-UHFFFAOYSA-N 0.000 description 11
- 206010037844 rash Diseases 0.000 description 8
- 239000002245 particle Substances 0.000 description 6
- 238000004880 explosion Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 239000007789 gas Substances 0.000 description 4
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 238000012935 Averaging Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 3
- 230000006399 behavior Effects 0.000 description 3
- 230000007423 decrease Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 230000001960 triggered effect Effects 0.000 description 3
- 206010015946 Eye irritation Diseases 0.000 description 2
- 238000009529 body temperature measurement Methods 0.000 description 2
- 239000003818 cinder Substances 0.000 description 2
- 239000003086 colorant Substances 0.000 description 2
- 230000001276 controlling effect Effects 0.000 description 2
- 238000001816 cooling Methods 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 229910052732 germanium Inorganic materials 0.000 description 2
- GNPVGFCGXDBREM-UHFFFAOYSA-N germanium atom Chemical compound [Ge] GNPVGFCGXDBREM-UHFFFAOYSA-N 0.000 description 2
- 238000010438 heat treatment Methods 0.000 description 2
- 230000007794 irritation Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 239000004065 semiconductor Substances 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 239000002341 toxic gas Substances 0.000 description 2
- 206010061218 Inflammation Diseases 0.000 description 1
- QAOWNCQODCNURD-UHFFFAOYSA-N Sulfuric acid Chemical compound OS(O)(=O)=O QAOWNCQODCNURD-UHFFFAOYSA-N 0.000 description 1
- 206010043521 Throat irritation Diseases 0.000 description 1
- 230000002547 anomalous effect Effects 0.000 description 1
- 208000006673 asthma Diseases 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 231100000013 eye irritation Toxicity 0.000 description 1
- 230000010006 flight Effects 0.000 description 1
- 230000004054 inflammatory process Effects 0.000 description 1
- 230000033001 locomotion Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 230000007257 malfunction Effects 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 210000002345 respiratory system Anatomy 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 238000009991 scouring Methods 0.000 description 1
- SBIBMFFZSBJNJF-UHFFFAOYSA-N selenium;zinc Chemical compound [Se]=[Zn] SBIBMFFZSBJNJF-UHFFFAOYSA-N 0.000 description 1
- 239000000377 silicon dioxide Substances 0.000 description 1
- 239000012798 spherical particle Substances 0.000 description 1
- 239000001117 sulphuric acid Substances 0.000 description 1
- 235000011149 sulphuric acid Nutrition 0.000 description 1
- 239000002887 superconductor Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G08—SIGNALLING
- G08B—SIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
- G08B21/00—Alarms responsive to a single specified undesired or abnormal condition and not otherwise provided for
- G08B21/02—Alarms for ensuring the safety of persons
- G08B21/12—Alarms for ensuring the safety of persons responsive to undesired emission of substances, e.g. pollution alarms
- G08B21/14—Toxic gas alarms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J5/60—Radiation pyrometry, e.g. infrared or optical thermometry using determination of colour temperature
-
- G—PHYSICS
- G08—SIGNALLING
- G08B—SIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
- G08B17/00—Fire alarms; Alarms responsive to explosion
- G08B17/12—Actuation by presence of radiation or particles, e.g. of infrared radiation or of ions
- G08B17/125—Actuation by presence of radiation or particles, e.g. of infrared radiation or of ions by using a video camera to detect fire or smoke
-
- G—PHYSICS
- G08—SIGNALLING
- G08B—SIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
- G08B21/00—Alarms responsive to a single specified undesired or abnormal condition and not otherwise provided for
- G08B21/02—Alarms for ensuring the safety of persons
- G08B21/10—Alarms for ensuring the safety of persons responsive to calamitous events, e.g. tornados or earthquakes
-
- G—PHYSICS
- G08—SIGNALLING
- G08B—SIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
- G08B21/00—Alarms responsive to a single specified undesired or abnormal condition and not otherwise provided for
- G08B21/02—Alarms for ensuring the safety of persons
- G08B21/12—Alarms for ensuring the safety of persons responsive to undesired emission of substances, e.g. pollution alarms
Definitions
- the present invention relates to an infrared detection apparatus for monitoring a field of view in order to detect an adverse atmospheric condition.
- the apparatus has particular application to detecting sulphur dioxide, volcanic ash and wind-blown dust.
- Volcanic ash is a hazard to jet aircraft, causing engines to stall when ingested, scouring windows and the leading edges of the wings and causing instrument malfunctions. Damage to aircraft can be counted in the millions of dollars. Most serious aircraft encounters with ash clouds have been at cruise altitudes, but there is also a hazard to aircraft at airports affected by volcanic ash. These airports are usually close to an active volcano (e.g. Anchorage and Kagoshima) but they can also be at some distance from the source of the eruption due to atmospheric transport that brings ash into the region.
- active volcano e.g. Anchorage and Kagoshima
- Wind-blown dust from desert regions or semiarid lands can be a hazard to aircraft, reduces visibility significantly and can cause eye and throat irritation to humans.
- Large parts of the habitable earth are prone to dust storms, including northern Africa, the Mediterranean islands, southern Italy, Spain and France, southwestern USA, central and southern Australia, western parts of South America, central China, Japan and south and north Korea and the central deserts of Asia.
- the wind-blown dust can also be transported long distances - dust from China has been detected in North America.
- the dust consists of nearly spherical particles of Si0 2 in concentrations that can limit visibility to a few 10 's of metres. Accordingly, wind-blown dust can be a significant hazard to aircraft, vehicles and the like.
- An infrared detection apparatus for detecting an adverse atmospheric condition comprising: a plurality of filters corresponding to different ones of a plurality of wavelengths and at least including filters which enable the adverse atmospheric condition to be detected; an infrared array, said infrared array producing signals representative of infrared radiation reaching said array; radiation control means for controlling the infrared radiation received by the infrared array, the radiation control means including means for changing the filters so that said infrared array can produce wavelength signals representative of infrared radiation from each of said wavelengths corresponding to the adverse atmospheric array, and means for enabling said infrared array to produce calibration signals for each wavelength signal; calibration means for performing a calibration of each wavelength signal to correct for radiation from the infrared detection apparatus on the basis of at least the corresponding calibration signal to thereby produce a calibrated wavelength signal representative of radiation from the region of sky; and output means for producing an output indicative of the presence of the adverse atmospheric condition in the region of sky based on the calibrated wavelength signals .
- the apparatus also includes correction means for correcting the calibrated signals for water vapour absorption and/or viewing angle effects
- the correction means comprises a pre-computed look-up table (LUT) .
- the LUT is based on off-line detailed radiative transfer calculations which account for the effects of water vapour absorption and for viewing geometry for each of the preferred filtered infrared wavelengths .
- the means for enabling said infrared array to produce calibration signals comprises a shutter having infrared emissivity which is high in each of said wavelengths .
- said radiation control means controls said shutter to shut prior to each wavelength measurement to thereby enable said array to produce a calibration measurement corresponding to the preceding wavelength measurement .
- the calibration means also calibrates on the basis of a pre-calibration.
- the calibration means alters the coefficient of the calibration equation on the basis of the calibration signal.
- calibration coefficients are calculated for each filter and for each pixel of the infrared array.
- said output means outputs temperature difference images derived from at least two wavelengths.
- said output means produces an alarm if the adverse atmospheric condition is present.
- said alarm is produced if the temperature differences exceeds a predetermined value.
- the algorithm for enabling the output means is dependent on the atmospheric condition. We have determined algorithms suitable for a number of atmospheric conditions.
- the algorithm is based on the temperatures T 8 .6# Tio.o * Tu, 0 and T12.0 at four wavelengths, 8.6 ⁇ m, 10.0 ⁇ m, 11.0 ⁇ m and 12.0 ⁇ m.
- said algorithm is based on a temperature difference between two wavelengths, Tn.o and T ⁇ 2 .n.
- Figure la is a block diagram that shows the apparatus of a preferred embodiment
- Figure lb is a schematic diagram which shows the camera portion of the apparatus of the preferred embodiment
- Figure 2 shows the filter functions of the apparatus
- Figure 3 shows the variation of elevation angle with temperature difference
- Figures 4a and 4b show variation of elevation angle with temperature difference for two different channel differences
- Figure 5 shows the apparatus of the invention viewing a sulphur dioxide plume
- Figure 6 shows an image produced of S0 2 plume using the apparatus of the present invention
- Figure 7 shows a view of a plume of S0 2 with an explosion
- Figure 8 shows an image of an explosion produced using an ash algorithm
- Figure 9 shows a different image for the sky when there is no ash or sulphur dioxide present
- Figure 10 is a temperature histogram for the image of Figure 9
- Figure 11 illustrates the Gaussian and thresholding technique for setting an ash alarm
- Figures 12a and 12b show a raw and a calibrated image respectively
- Figure 17 has ash and visible images from Matupit;
- Figure 18 shows the alarm histogram for Figure 17;
- Figure 19 is a view of an eruption from Rabaul Volcanological Officer;
- Figure 20 is a graph of an alarm time-series;
- Figure 21 has ash and visible images from Hamamas Hotel; and
- Figure 22 shows the alarm histogram for Figure 22.
- the preferred embodiment provides an infrared detection apparatus for monitoring infrared signals from a field of view in order to detect an adverse atmospheric condition such as volcanic ash, sulphur dioxide or wind-blown dust.
- FIG la is a block diagram of the infrared detection apparatus of a preferred embodiment.
- the infrared detection apparatus 30 comprises an infrared camera 110, processor 140 and satellite modem 160 for transmitting signals from the infrared detection apparatus 130 to a monitoring station 150.
- the camera portion 110 of infrared detection apparatus 130 is shown in Figure 1.
- the camera 110 has a filter wheel housing 2, a window 3 that is transmissive in the infrared (e.g. a Germanium window) , a shutter 4, a filter wheel cover 6 and a filter wheel.
- Infrared array housing 8 contains an infrared array 9 and a signal processing unit 10.
- Output from the signal processing unit 10 is on signal lines 11a which have an Ethernet interface to the computer 140 so that it can process the signals.
- Lines lib allow the shutter control signals to be received by the infrared detection apparatus 130 from the computer 140 and also for equipment temperature measurements passed to the computer 140.
- the temperature measurement corresponds to the temperature of the shutter 4 and are obtained by a contact thermometer (not shown) .
- the infrared detection apparatus 130 operates by monitoring infrared signals from a field of view (e.g. a region of sky) being monitored at up to five pre-defined wavelengths. While the wavelengths are monitored, the wavelengths which are used in detecting the adverse atmospheric condition depend on the adverse atmospheric condition which is being monitored.
- the central wavelength and wavelength intervals are given in Table 1. It will be appreciated that a band of wavelengths surround a central wavelength, but for convenience the terms "wavelength" or "channel” are used herein to refer to the central wavelength and surrounding band unless the context implies otherwise.
- the infrared radiation measured by the camera 110 is linearly proportional to the resistance change in the detector, which is recorded and logged by the signal processing unit 10.
- the infrared array is an uncooled microBolometer staring array of 320 x 240 elements sensitive to radiation in the 6-14 ⁇ wavelength interval is used to detect filtered radiation.
- the detection apparatus uses a filter wheel to filter radiation.
- the radiation from the field of view is focussed onto the array by means of focussing optics in the form of lens 5 and the field of view is a cone of up to 90 degrees.
- the infrared array will be an uncooled microBolometer array of dimensions at least 320 x 240 elements, but 640 x 480 elements is also possible. There is a trade-off between the number of elements, cost and maximum spatial resolution per pixel for a fixed optical arrangement.
- the microBolometer operates on the principle that a temperature change produced by radiation falling on the detector produces a linear resistance change in the material.
- the camera 110 can be used to obtain measurements at up to five separate wavelengths to be filtered and also can have a single broadband channel depending on what atmospheric condition is being monitored.
- the camera 110 uses a filter wheel mounted with circular interference filters. Table 1 provides the information for the selection of the filters .
- the filters are 50 mm diameter germanium/ZnSe multi-layer interference filters mounted on a rotating wheel and driven by a stepper-motor. (But, smaller or larger diameter filters may be used depending on the field of view required and the focusing power of the optics) .
- the filter wheel provides a radiation control means for controlling the infrared radiation received by the infrared array 9.
- the array 9 has a nominal noise temperature of no less than »50 mK in the broadband channel.
- frame averaging is employed.
- Table 1 shows the theoretical noise equivalent temperatures (NE ⁇ T's) expected for various frame averaging values in 5 narrow wavelength bands or channels and one broadband channel 20.
- N ⁇ T's noise equivalent temperatures
- Table 1 shows noise equivalent temperatures over 1, 4, 8, 16, 32 and 64 frames. It will be seen that for each doubling in the number of frames the noise equivalent temperature is reduced by approximately 30%.
- the raw data produced by the infrared array is calibrated (or corrected) .
- the apparatus of the preferred embodiment is calibrated so that the processing means can produce corrected radiance values which then can be used to produce scene temperatures which can subsequently be processed using an algorithm specific to the atmospheric condition in order to determine the presence of the adverse atmospheric conditions.
- the calibration process consists of a pre-calibration and a field calibration. The calibrations correct for radiation from the infrared detection apparatus.
- the field calibration is required in order to correct for changes in the radiation from the infrared detection apparatus during operation.
- the field calibration technique is a key feature which allows sufficiently accurate determination of the presence (or not) of the adverse atmospheric conditions.
- Figure 12a shows a typical scene where the image is constructed from the raw signals, without calibration.
- Figure 12b shows the same scene this time after calibration and conversion to temperature units. Aspects of the scene not visible in the raw data are now clearly noticeable in the calibrated data. For example, roof 200 is now visible.
- the camera 110 provides raw digital counts as output of detector array 9 that have also had some corrections applied. These counts can be related to the scene radiance through a linear calibration process, and scene then to temperature through use of the Planck function.
- Equation (3) To convert from the calibrated radiances to scene temperature, Equation (3) is inverted. This is a non-linear problem which requires a minimisation procedure. A series of look-up tables were generated that give radiances equivalent to a series of pre-specified temperatures. Once the measured radiance is known by combining (1) and (2) , the look-up table is searched and interpolated (if necessary) to determine the closest scene temperature. The procedure is accurate to 10 mK over the range of observable temperatures 220 K to 330 K.
- the procedure to pre-calibrate the camera 110 is:
- a separate set of calibration coefficients a ⁇ , b ⁇ is developed for each pixel within the 320 x 240 image. Towards the edges of the image the quality of the calibration degrades due to image distortions and non-uniformity of the blackbodies.
- the blackbody temperature is measured in one place on each blackbody and non-uniformity of the temperature field will occur to some degree .
- the pre-calibration procedure was repeated several times and average look-up tables were generated. Typically the pre-calibration is performed in the laboratory and not done in the field.
- a field calibration technique is used to alter the coefficients to account for the optics (particularly the lens) that may be heating up or cooling down and thus be at a different temperature to its value when calibrated in the laboratory. This causes an off-set in the measured signals.
- Our innovative field calibration procedure makes use of a single shutter measurement just before measurements are taken at each wavelength. The shutter fills or slightly overfills the field of view of the instrument and provides a uniform radiation source to the detector. The temperature of the side of the shutter facing the lens is continuously monitored using a contact temperature probe.
- the shutter side facing the lens is blackened so that its infrared emissivity is high (exceeding 0.98) and uniform across the region 6-14 ⁇ m.
- the calibration is performed by making a single measurement of the shutter, followed by a measurement of the scene and then application of the calibration equations and shutter measurement which accounts for the off-set generated by any change in temperature of the lens or other radiating surfaces in front of the detector.
- the raw data is also calibrated for background atmospheric conditions and viewing angle by the calibration means. That is, the temperature differences on a single channel will vary depending on the channel measured.
- Figure 4a shows the difference as determined from measurements made at Saipan.
- the variation with elevation angle mimics the theoretical behaviour.
- the same effect with elevation can be seen for 8.6-12 ⁇ m temperature differences (Fig. 4b), except that after 60 degrees the difference starts to increase rather than decrease. This is not seen in the modelling results, and more data are required to determine the cause of this effect.
- temperature difference can be corrected in order to correct the output images so that the output temperature difference images correctly reflect the presence of the adverse atmospheric condition being detected.
- Figure 5 shows a digital image of the camera viewing towards the Etna volcano in the background with a plume of S0 2 30 emitted from the crater.
- An S0 2 image produced by the apparatus within 30 minutes of the digital image is shown in Figure 6.
- the colour scale 22 is drawn to indicate the amount of S0 2 in the plume - from yellow, indicating low amount, to brown, indicating high amounts.
- the background to the images comprises light areas of blue and green indicating a colder background whereas the bottom of the image which is dark in colour represents the ground. These areas are labelled 33 and 34. respectively.
- the left vertical axis represents elevation in degrees and the horizontal axis represents Azimuth.
- the images are typically produced in colour as indicated above however it will be appreciated that appropriate grey scale images can also be produced. The colour images make it easier to discern the plumes from the background.
- An S0 2 index is based on a 4 -channel algorithm, while the temperature difference images utilize all 5 channels, the 4 narrow band channels and a wide band channel .
- the S0 2 index is derived by: (1) forming the temperature difference between a channel centred at 8.6 ⁇ m and a channel centred at lO.O ⁇ m, label this difference as dl l f (2) forming the temperature difference between a channel at 11.0 ⁇ m and a channel centred at 12.0 ⁇ m, label this difference as ⁇ i (3) adding temperature differences d ⁇ i and 8l 2 ⁇ label this d ⁇ 3 , subtracting a reference value that depends on the viewing elevation of the camera, and has a typical range of 1-3°C, label this ⁇ T 4 to get the final temperature difference (this is the S0 2 index) .
- the displayed temperature difference image is produced by scaling 5r 4 and overlaying this scaled image into the broadband image so that all pixels in the ⁇ 5r 4 temperature difference image with a scaled value in the range 1-32 are preserved and all pixels outside this range are replaced by image pixel from the broadband filter image.
- a suitable colour table is then attached to the image and a reference grid and scale are incorporated.
- the resulting image ( Figure 6) shows S0 2 plume 30 in yellow to brown colour, water vapour (in various degrees of amount) in grey colours, the background (colder) sky as blue and green 33 and topographic features (mountain, ground, trees etc., which are generally warmer than the plume) as dark grey to black 34.
- ⁇ T t is a temperature threshold that depends on the water vapour content of the atmosphere and on the viewing elevation angle of the camera.
- the nominal value for ⁇ T t is 0°C.
- Data (pixels) with values above the threshold are regarded as volcanic ash.
- Data (pixels) below the threshold are regarded as not volcanic ash.
- the apparatus was also able to capture discrete explosions from the Stromboli Volcano.
- An example of this is shown in Figure 7.
- the pyroclastic material is mostly volcanic hot rocks, cinders and ash and reveals itself as grey to black colours 41 when the S0 2 algorithm is used.
- An S0 2 plume 42 is clearly shown.
- the ash algorithm that is, by taking temperature differences using two channels, specifically ⁇ T 2 , the image shown in Figure 7 is obtained.
- the colour scale now shows positive temperature differences in shades of orange and red, and negative differences as blue to yellow.
- the algorithm identifies the hot rocks and cinders as positive differences 40 (high ash content) , and resuspended ash as slightly negative (similar to the material on the surface of the mountain slopes) .
- the sky has markedly negative differences .
- Desert dust has a high silica (Si0 2 ) content and when small particles (diameters less than 10 ⁇ m) are suspended in the atmosphere they disperse infrared radiation in a similar fashion to volcanic ash particles. Consequently, the algorithm used to identify ash in the atmosphere can also be used to identify wind-blown dust.
- Dust storms are a frequent and global phenomenon. Most of the dust is confined to the boundary layer - the part of the atmosphere closest to the surface and generally not extending more than 5 km upwards. Occasionally, large dust storms can be transported vast horizontal distances (may 1000 's of kilometres) and be lifted to heights greater than 5 km. Dust storms have been identified using passive infrared radiation from satellites.
- the dust algorithm differs from the ash algorithm in one important aspect. Since it is unlikely that wind-blown dust will contain any appreciable amounts of S0 2 gas (the reverse being true for ash) , a channel at 8.6 ⁇ m, can be used in conjunction with the 11 and 12 ⁇ m channels. The dust algorithm thus uses three channels rather than two.
- an automated algorithm can be developed in order to initiate an alarm either so the person can be alerted to the need to monitor the alarm or to send data from the apparatus 130 to a monitoring station 150.
- This alarm is based on analysis of the difference images produced in accordance with the algorithms.
- the images do not actually have to be produced in order for the alarm to be triggered, that is, the data can be processed to determine whether an alarm condition is met.
- An automated alarm algorithm has been developed minimise the requirement for operator intervention of the infrared image data and/or to trigger the transmission of an image to the users - i.e. to transmit the image data to monitoring station 150 when there is sufficient reason for a user to inspect the image.
- the algorithm or ⁇ alarm' is based on a histogramming technique that takes into account the viewing elevation and the amount of water vapour in the atmosphere.
- the histogramming technique accounts for these anomalies.
- Figure 9 was obtained at 20 degrees elevation and as the field of view of the infrared camera is roughly 24 degrees in the vertical direction, some parts of the image view land surfaces.
- the histogram in Figure 10 has prominent peaks at roughly -1 K 51 and -5 K 50 which correspond to clouds and clear skies, respectively. In this case the least negative peak has a tail that includes some positive pixels. In the corresponding image ( Figure 9) these pixels are due to viewing features that are low on the horizon and include ground targets. Such anomalies' are difficult to isolate in an automated manner and could give rise to false alarms if a straightforward pixel thresholding technique were employed.
- the scheme chosen to automatically determine whether an image has detected ash is a statistically based method. This is the method of choice because by the nature of the problem there is often going to be a distribution of pixels that can be flagged as ash, within an image that has many pixels that are definitely ash or definitely not ash. In addition, because of the likelihood that pixels will contain mixtures, a simple threshold and binary decision process would be inappropriate.
- the histogram shown in Figure 10 has two prominent peaks with a spread of pixels around these peaks. If the detection apparatus 130 viewed a target of constant temperature (e.g. a uniform cloud or the clear sky), then simply because of the fact that the camera has a wide field of view and there is water vapour absorption along the differing paths to the target, the resulting difference image would be non-uniform. In practice it is unlikely that the sky would present a uniform target and even less likely that a cloud would by perfectly uniform. The combination of these effects leads to a natural spread in the histogram of the temperature differences, with a central peak corresponding to the mode temperature difference. For a relatively uniform scene the peak would be high and the spread (or standard deviation of the distribution) would be low. We have selected a Gaussian distribution to model the distribution. The Gaussian distribution in the mathematical terms is,
- ⁇ T is the temperature difference
- ⁇ ⁇ ⁇ is the mean temperature difference
- ⁇ ⁇ is the standard deviation
- a set of Gaussian distribution is fitted to the frequency distribution data and the parameters, &o, ⁇ , ⁇ ,i, and ⁇ ⁇ ,i derived. The linear combination of these distributions is the model- fit to the data.
- a threshold Gaussian is set with a mean and standard deviation derived from modelling. This threshold Gaussian is compared with the n-Gaussian data- fit. The region between the pixels bounded by the threshold Gaussian mean value, and the overlap region between the two Gaussians (the threshold and the data- fit) is calculated. This area (or number of pixels) is subtracted from the number of pixels that exceed the threshold Gaussian mean value and lie within the data- fit Gaussian (see Fig. 11) .
- P 0r ⁇ is the number of overlap pixels for Gaussian i
- P ⁇ is the number of pixels that exceed the threshold mean
- a 0 ,i are the maxima for the Gaussian fits. The purpose of normalising by the maximum is to ensure that more weight is given to distributions that have well-defined and dominant peaks .
- results indicate that the detection apparatus can image ash plumes and clouds and clearly discern these from meteorological clouds. Results are best at closest proximity to the ash cloud, but good results were obtained at distances greater than 5 km from the active crater. The ash alarm algorithm was also tested in an autonomous manner overnight from a distance of ⁇ 8 km.
- Figure 13 shows the locations of the measurement sites (six in all) used to image the ash-rich eruptions from Tavurvur. They are listed in Table 2.
- Tavurvur has been active since a major eruption took place in September 1994 which threatened the town of Rabaul and destroyed the airport.
- a new airport, Tokua has been constructed and is located about 20 km southwest of Tavurvur. With the crater still active, flights in and out of Tokua only take place in daylight hours and not at all if the winds move the ash towards the airport runway.
- Figure 14 shows a digital photograph taken from the runway at Tokua. A plume 60 from Tavurvur is noticeable in the background.
- the Rabaul Volcanological Observatory operates on a hill overlooking the active crater and at about 8 km distance from it. Economic pressures in PNG have meant that only limited resources are available at RVO for operating geophysical equipment and power failures are also common. The main means of transport throughout PNG is by jet and light aircraft and the economy is highly dependent on air transportation. Thus there is an urgent need to monitor the volcanoes in New England (there are many) and throughout PNG.
- the apparatus of the preferred embodiment operates off batteries for up to 16 hours and can be deployed in relatively hostile environments, rapidly by a single user.
- the apparatus was deployed at a variety of locations and in a variety of viewing configurations .
- the measurements were made in an atmosphere with quite high water vapour amounts and at elevation angles varying from 10° to nearly vertical viewing.
- the atmosphere around the instrument was filled with ash particles, making the atmosphere appear grey and causing irritation to the eyes and lungs .
- Figure 15 shows a typical results obtained from Rababa
- the ash image (left-panel) correctly identifies the plume 71 and clouds of ash from Tavurvur. Grey to black coloured regions of the image are identified as having no ash. The mountainside 72 is also identified as ash — this is not surprising since the mountain is covered in ash particles .
- the automatic alarm algorithm was used on all images and the alarm generated for the image shown in Fig. 15 is shown in Fig. 16.
- the alarm is being generated because of the difference between the actual histogram 210 and threshold histogram 211. About 43% of the pixels are identified as ash in the image and a clear and unambiguous ash alarm has been triggered.
- Fig. 17 shows an ash image and corresponding photograph of a dispersing ash cloud 75 viewed from 40° elevation. In this case 73% of the pixels were identified as ash (see Fig. 18) .
- Ash images from RVO (E) were the most difficult to obtain because of the distance from the active vent ( ⁇ 8 km) and because the Observatory is perched on a hill; thus the camera could only view at relatively small elevation angles ( ⁇ 10° or less) .
- the combination of the greater distance and small elevation angles means that a considerable water vapour path is traversed (the water vapour absorption masks the positive temperature differences expected from the ash signal) .
- the frequency of eruption was so high that the air between the camera and the eruption column was often filled with fine ash particles. This has the effect of making the atmosphere appear grey and also causes large absorption in the infrared.
- the camera was operated continuously overnight at RVO.
- An example of the scene viewed by the apparatus from RVO is shown in Figure 19.
- Figure 20 shows a time-series of alarms detected by the apparatus from RVO.
- the series of triangles and circles 81 represent alternate 5-min sampled data and their similarity gives confidence in the results.
- the threshold for the alarm has been arbitrarily set to a value of 10%. In practice this will have to be amended to suit the viewing conditions.
- the plot suggests that there were continuous ash emissions during the night - in agreement with what was observed during the day. Notice that the highest alarm percentages never exceed 35% or so and this is simply a function of the viewing attitude of the instrument. If the instrument were sited closer to the volcano, then more of the ash would fill the field of view of the instrument and there would be a higher percentage for the alarm. Again this is a feature that need to be considered within the context of the viewing geometry and location of the instrument with respect to the activity.
Landscapes
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Emergency Management (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Toxicology (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Multimedia (AREA)
- Radiation Pyrometers (AREA)
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
AU2003905315A AU2003905315A0 (en) | 2003-09-29 | An Alarm System for Remote Sensing Equipment | |
AU2004900214A AU2004900214A0 (en) | 2004-01-16 | An infrared detection apparatus | |
PCT/AU2004/001338 WO2005031323A1 (en) | 2003-09-29 | 2004-09-29 | An infrared detection apparatus |
Publications (1)
Publication Number | Publication Date |
---|---|
EP1673613A1 true EP1673613A1 (de) | 2006-06-28 |
Family
ID=34394898
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
EP04761373A Withdrawn EP1673613A1 (de) | 2003-09-29 | 2004-09-29 | Infrarotnachweisvorrichtung |
Country Status (2)
Country | Link |
---|---|
EP (1) | EP1673613A1 (de) |
WO (2) | WO2005031321A1 (de) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106781430A (zh) * | 2016-11-15 | 2017-05-31 | 北京空间机电研究所 | 一种高灵敏度红外遥感器性能测试装置 |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2564215A1 (de) | 2010-04-29 | 2013-03-06 | Norsk Institutt For Luftforskning | System und verfahren zur erkennung von ungünstigen wetterbedingungen im vorfeld eines flugzeuges |
WO2012170093A2 (en) * | 2011-03-25 | 2012-12-13 | Exxonmobil Upstream Research Company | Autonomous detection of chemical plumes |
WO2015199911A1 (en) | 2014-06-23 | 2015-12-30 | Exxonmobil Upstream Research Company | Methods and systems for detecting a chemical species |
US9442011B2 (en) | 2014-06-23 | 2016-09-13 | Exxonmobil Upstream Research Company | Methods for calibrating a multiple detector system |
WO2015199912A1 (en) | 2014-06-23 | 2015-12-30 | Exxonmobil Upstream Research Company | Image quality enhancement of a differential image for a multiple detector system |
US9448134B2 (en) | 2014-06-23 | 2016-09-20 | Exxonmobil Upstream Research Company | Systems for detecting a chemical species and use thereof |
JP6750672B2 (ja) * | 2016-04-20 | 2020-09-02 | コニカミノルタ株式会社 | ガス観測方法 |
CN114018982B (zh) * | 2021-10-14 | 2023-11-07 | 国网江西省电力有限公司电力科学研究院 | 一种空预器积灰可视化监测方法 |
CN114235690B (zh) * | 2021-11-25 | 2023-11-21 | 中国人民解放军空军工程大学 | 航空器涂层的面红外发射率测量方法及装置 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US565470A (en) * | 1896-08-11 | Envelop | ||
JPS5435426A (en) * | 1977-08-24 | 1979-03-15 | Showa Yuka Kk | Apparatus for monitoring flame from flare stack |
US4725733A (en) * | 1983-07-18 | 1988-02-16 | The United States Of America As Represented By The Secretary Of The Navy | Apparatus and method for remotely detecting the presence of chemical warfare nerve agents in an air-released thermal cloud |
US4965572A (en) * | 1988-06-10 | 1990-10-23 | Turbulence Prediction Systems | Method for producing a warning of the existence of low-level wind shear and aircraftborne system for performing same |
SU1712837A1 (ru) * | 1989-10-20 | 1992-02-15 | Конотопский Электромеханический Завод "Красный Металлист" | Способ автоматического контрол запыленности шахтной атмосферы |
WO1991015739A1 (en) * | 1990-04-09 | 1991-10-17 | Commonwealth Scientific And Industrial Research Organisation | A detection system for use in an aircraft |
FR2696939B1 (fr) * | 1992-10-16 | 1995-01-06 | Bertin & Cie | Procédé et dispositif de détection automatique rapide de feux de forêt. |
US6023061A (en) * | 1995-12-04 | 2000-02-08 | Microcam Corporation | Miniature infrared camera |
DE29917646U1 (de) * | 1999-10-06 | 1999-12-23 | Hwg Telekommunikationssysteme | Multifunktionelle Fernüberwachungsvorrichtung |
-
2004
- 2004-09-29 EP EP04761373A patent/EP1673613A1/de not_active Withdrawn
- 2004-09-29 WO PCT/AU2004/001340 patent/WO2005031321A1/en active Application Filing
- 2004-09-29 WO PCT/AU2004/001338 patent/WO2005031323A1/en not_active Application Discontinuation
Non-Patent Citations (1)
Title |
---|
See references of WO2005031323A1 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106781430A (zh) * | 2016-11-15 | 2017-05-31 | 北京空间机电研究所 | 一种高灵敏度红外遥感器性能测试装置 |
Also Published As
Publication number | Publication date |
---|---|
WO2005031323A1 (en) | 2005-04-07 |
WO2005031321A1 (en) | 2005-04-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10440291B2 (en) | System and method for detecting adverse atmospheric conditions ahead of an aircraft | |
EP0525007B1 (de) | Detektionssystem zum gebrauch im flugzeug | |
Minnis et al. | Near-real time cloud retrievals from operational and research meteorological satellites | |
JP6680788B2 (ja) | 空の観察領域内の雲を特定及び監視するための検出装置及び方法 | |
EP1673613A1 (de) | Infrarotnachweisvorrichtung | |
AU2011100797A4 (en) | System and method for detecting adverse atmospheric conditions ahead of an aircraft | |
Viollier et al. | Radiometric and spectral characteristics of the ScaRaB-3 instrument on Megha-Tropiques: comparisons with ERBE, CERES, and GERB | |
AU2004276374A1 (en) | An infrared detection apparatus | |
Dandini et al. | Halo ratio from ground-based all-sky imaging | |
AU2004276376A1 (en) | Apparatus for remote monitoring of a field of view | |
Putra et al. | Visualization of volcanic ash distribution based on multispectral satellite imagery: A comparing method | |
Mace | Level 2 GEOPROF product process description and interface control document | |
Halthore et al. | Overview of atmospheric correction and radiometric calibration efforts during FIFE | |
AU2012359085B2 (en) | System and method for detecting adverse atmospheric conditions ahead of an aircraft | |
Thurairajah | Thermal infrared imaging of the atmosphere: The Infrared Cloud Imager | |
AU654666B2 (en) | A detection system for use in an aircraft | |
Ahn et al. | Effects of the tropospheric dust on the sea surface temperature derivation from the GMS-5 IR data | |
Vogelmann | RACORO data guide | |
Hall et al. | Characterization of aerosol-containing chemical simulant clouds using a sensitive, thermal infrared imaging spectrometer | |
Platt et al. | A case study of cirrus layers with variable 3.74 µm reflection properties in the first FIRE experiment, 2 November 1986 | |
Schade et al. | IR Sky Radiance Distributions In The Marine Boundary Layer | |
IES85926Y1 (en) | System and method for detecting adverse atmospheric conditions ahead of an aircraft | |
IE20110213U1 (en) | System and method for detecting adverse atmospheric conditions ahead of an aircraft | |
Noel | Utilization of IR imagery in tactical reconnaissance | |
Adhikari | Cloudy condition assessment within an AIRS pixel by combining MODIS and ARM ground-based lidar and radar measurements |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PUAI | Public reference made under article 153(3) epc to a published international application that has entered the european phase |
Free format text: ORIGINAL CODE: 0009012 |
|
17P | Request for examination filed |
Effective date: 20060427 |
|
AK | Designated contracting states |
Kind code of ref document: A1 Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IT LI LU MC NL PL PT RO SE SI SK TR |
|
DAX | Request for extension of the european patent (deleted) | ||
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: THE APPLICATION IS DEEMED TO BE WITHDRAWN |
|
18D | Application deemed to be withdrawn |
Effective date: 20080401 |