EP1642087B1 - Method and apparatus for automatically detecting and mapping, particularly for burnt areas without vegetation - Google Patents
Method and apparatus for automatically detecting and mapping, particularly for burnt areas without vegetation Download PDFInfo
- Publication number
- EP1642087B1 EP1642087B1 EP04745201A EP04745201A EP1642087B1 EP 1642087 B1 EP1642087 B1 EP 1642087B1 EP 04745201 A EP04745201 A EP 04745201A EP 04745201 A EP04745201 A EP 04745201A EP 1642087 B1 EP1642087 B1 EP 1642087B1
- Authority
- EP
- European Patent Office
- Prior art keywords
- image
- images
- nir
- vegetation
- reflectance
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 67
- 238000013507 mapping Methods 0.000 title claims abstract description 14
- 230000003595 spectral effect Effects 0.000 claims abstract description 39
- 239000002028 Biomass Substances 0.000 claims description 18
- 230000007423 decrease Effects 0.000 claims description 17
- 238000012545 processing Methods 0.000 claims description 12
- 238000001514 detection method Methods 0.000 claims description 11
- 230000000694 effects Effects 0.000 claims description 8
- 230000001932 seasonal effect Effects 0.000 claims description 8
- 239000002689 soil Substances 0.000 claims description 7
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 6
- 239000011707 mineral Substances 0.000 claims description 6
- 230000003287 optical effect Effects 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 3
- 230000004044 response Effects 0.000 abstract description 20
- 238000001228 spectrum Methods 0.000 abstract description 6
- 238000012731 temporal analysis Methods 0.000 abstract description 6
- ZPUCINDJVBIVPJ-LJISPDSOSA-N ***e Chemical compound O([C@H]1C[C@@H]2CC[C@@H](N2C)[C@H]1C(=O)OC)C(=O)C1=CC=CC=C1 ZPUCINDJVBIVPJ-LJISPDSOSA-N 0.000 description 9
- 238000004458 analytical method Methods 0.000 description 6
- 238000002834 transmittance Methods 0.000 description 6
- 239000011435 rock Substances 0.000 description 4
- 230000004075 alteration Effects 0.000 description 3
- 238000005286 illumination Methods 0.000 description 3
- 240000007817 Olea europaea Species 0.000 description 2
- 239000005427 atmospheric aerosol Substances 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 229930002875 chlorophyll Natural products 0.000 description 2
- 235000019804 chlorophyll Nutrition 0.000 description 2
- ATNHDLDRLWWWCB-AENOIHSZSA-M chlorophyll a Chemical compound C1([C@@H](C(=O)OC)C(=O)C2=C3C)=C2N2C3=CC(C(CC)=C3C)=[N+]4C3=CC3=C(C=C)C(C)=C5N3[Mg-2]42[N+]2=C1[C@@H](CCC(=O)OC\C=C(/C)CCC[C@H](C)CCC[C@H](C)CCCC(C)C)[C@H](C)C2=C5 ATNHDLDRLWWWCB-AENOIHSZSA-M 0.000 description 2
- 239000002131 composite material Substances 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 239000000443 aerosol Substances 0.000 description 1
- 239000003054 catalyst Substances 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000000513 principal component analysis Methods 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000005437 stratosphere Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 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 1
- 238000012546 transfer Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C11/00—Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
- G01C11/04—Interpretation of pictures
- G01C11/06—Interpretation of pictures by comparison of two or more pictures of the same area
-
- 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
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J3/30—Measuring the intensity of spectral lines directly on the spectrum itself
- G01J3/32—Investigating bands of a spectrum in sequence by a single detector
-
- 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
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J3/30—Measuring the intensity of spectral lines directly on the spectrum itself
- G01J3/36—Investigating two or more bands of a spectrum by separate detectors
Definitions
- the present invention refers to an automatic method of detecting and mapping, particularly for burnt areas without vegetation, enabling, in a precise, simple, versatile, reliable, and efficient way, the discrimination of areas, in images acquired from satellite, where fires occurred on the basis of a comparison of images overcoming possible problems of dishomogeneity arising in the case when the conditions of acquisition of the images to be compared are different.
- the automatic mapping method is based on the multi-temporal analysis of multi spectral data, or on the changes in the spectral response of the Earth features in combined near infrared, or NIR (Near Infra-Red), and short wave infrared, or SWIR (Short Wave Infra-Red), or medium infrared, or MIR (Medium Infra-Red), windows of the electromagnetic spectrum.
- NIR Near Infra-Red
- SWIR Short Wave Infra-Red
- MIR Medium Infra-Red
- the present invention further refers to an apparatus performing such a method, which is simple, fast, and inexpensive.
- a spaceborne hot temperature event detection arrangement is described in document EP-A-1 048 928.
- Multi spectral radiometers are generally constituted by one or more passive optical sensors that register the radiance emitted by the Earth in various windows of the electromagnetic spectrum.
- the principal windows of the electromagnetic spectrum commonly used for land and vegetation analysis are the Visible (with wavelength from 0.4 ⁇ m to 0.7 ⁇ m, hereinafter referred to as VIS), the Near Infra-Red (with wavelength from 0.7 ⁇ m to 1.3 ⁇ m, hereinafter referred to as NIR), the Short Wavelength Infra-Red (with wavelength from 1.3 ⁇ m to 2.5 ⁇ m, hereinafter referred to as SWIR), the Medium Infra-Red (with wavelength from 2.5 ⁇ m to 4.5 ⁇ m, hereinafter referred to as MIR), and the Thermal Infra-Red (with wavelength up to 14.5 ⁇ m, hereinafter referred to as TIR).
- VIS Visible
- NIR Near Infra-Red
- SWIR Short Wavelength Infra-Red
- MIR Medium Infra-Red
- TIR Thermal Infra-Red
- Radiometers operating at these wavelengths are installed on airborne platforms - such as microlights, airplanes, balloons and dirigibles, operating between the ground free surface and the stratosphere - and spaceborne platforms, including all type of orbiting and sub-orbital spacecrafts at altitudes up to about 36.000 Km above the Earth's ground free surface.
- All the afore-said systems look at the Earth orthogonally to the direction of the platform's motion with an angle ranging from 0 to 35 degrees with respect to Nadir - that is, perpendicularly to the Earth.
- the aforementioned angle is called "Off Nadir".
- resolution cells or grids of the Earth surface can be generated with a spatial resolution (that is, the spacing of the elementary area, or pixel, within which it is possible to distinguish two objects) ranging from few centimetres to several kilometres.
- multi spectral radiometers are coherent imaging systems. This means that they allow measuring radiated energy of each Earth surface cell at different wavelengths, as well as measuring the difference in radiated energy between two independent cells in each wavelength.
- converting radiance in reflectance allows measuring the time variation in energy at each wavelength and in each Earth surface cell, and measuring the time variation of energy difference between two independent cells at each wavelength.
- Lp is the path radiance (that is the additional radiance due to the path in the atmosphere and hence to the diffusion or scattering phenomenon)
- E d is the downward irradiance
- E 0 is exo-atmospheric solar irradiance
- ⁇ z is the solar zenith angle.
- the atmospheric transmittance T v from the target toward the sensor depends on the off-nadir angle.
- the atmospheric transmittance T z in the illumination direction is a function of ⁇ z , and of the optical thickness to be accounted for in Rayleigh scattering and atmospheric aerosol effects.
- Radiance Lp can be evaluated by using the so-called Dark Object Subtraction or DOS method, which, as shown by P.S. Chavez in "An improved Dark-Object Subtraction Technique for Atmospheric Scattering Correction of Multispectral Data", Remote Sensing of the Environment 24, pp. 459-479, 1988, relies upon the hypothesis that some ground types - as some types of rocks and soils - behave as almost perfect "black bodies": in such cases, the associated reflectance measured at the satellite is almost entirely due to scattered path reflectance.
- Atmospheric radiative transfer codes ( Lowtran, Modtran, Code-6S, etc.) make use of standard atmosphere models, or of available measures of vertical atmosphere profiles, for determining the downward irradiance E d .
- the accuracy of atmospheric corrections for Rayleigh scattering and aerosol effects strictly depend upon the availability of data on atmosphere physical and chemical parameters at the precise time of the remote sensing image acquisition
- the fundamental key is to distinguish reflectance changes originated by illumination and atmospheric conditions, by sensor drift or by sensor differences, from reflectance changes due to decrease in biomass content, to chlorophyll absorption, to water content and to density in the vegetation cover.
- a method for separating changes due to intrinsic reflectance value of vegetation from changes due to other factors influencing the at-satellite reflectance is that of identifying "pseudo-invariant” features to be used as "reference targets” in cross-rectification of reflectance in different images. Identification of "pseudo-invariant" features requires accurate analysis of the whole of remote sensing data involved in the multi temporal analysis, following a linear radiometric transformation and aimed to implement relative reflectance retrieval. Accuracy of radiometric rectification strongly depends on the selection of pseudo-invariant features or of reference targets.
- mapping burnt areas starting from multi spectral remote sensing data were carried out on single date image analysis, on the basis of so-called Principal Components, band ratio or vegetation index, unsupervised or supervised classification clustering.
- use of a single remote sensing image for mapping burnt and de-vegetated areas does not allow these ones to be distinguished from urban areas, water surfaces, sparsely vegetated areas and rock outcrops.
- Multi-temporal methods based on the evaluation of changes of the spectral response between NIR and SWIR channels seem advantageous, since fire is an agent of land cover alteration, and therefore change detection methods reduce the likelihood of confusion with non vegetated or low vegetated static land cover types, such as urban areas, and rock outcrops
- NDVI index Normalized Difference Vegetation Index
- said one or more detection frequency channels may be at least two infrared channels, preferably a near infrared or NIR (Near Infra Red) channel, and a short wave infrared or SWIR (Short Wave Infra Red) channel, and/or a medium infrared or MIR (Medium Infra Red) channel.
- NIR Near Infra Red
- SWIR Short Wave Infra Red
- MIR Medium Infra Red
- Permanent Reflectors PRs may be determined as those resolution cells p xy which, in both images A i and A i +1 of the pair, have values of reflectance ⁇ in each one of said one or more channels which correspond to densely built urban areas, preferably presenting no or little vegetation with respect to the size of the resolution cells.
- said one or more first quantities depending on a variation, considered in steps C, D, F, and G, may comprise at least one quantity identically equal to the corresponding variation.
- said one or more first quantities depending on a variation, considered in steps C, D, F, and G comprise at least one percentage of the corresponding variation.
- said one or more second quantities depending on at least one reflectance ⁇ may comprise at least one quantity identically equal to the corresponding reflectance ⁇ of a channel.
- said one or more second quantities depending on at least one reflectance ⁇ may comprise at least one ratio between two reflectances ⁇ of two channels.
- the acquired data may comprise the radiance L of each one of the N images A i on each one of said one or more detection frequency channels
- the method before executing step B, may transforms the radiance L into the reflectance ⁇ for each one of said one or more channels, preferably taking account of the exo-atmospheric solar irradiance E 0 and of the effects of the sun zenith angle ⁇ z , and comprising the estimation of the path radiance L p , more preferably through the DOS method.
- the method takes vegetation seasonal differences into account.
- the method may take vegetation seasonal differences into account through the use of one or more corrective factors, preferably depending on the solar elevation angles ⁇ SEL_i and ⁇ SEL_ ( i+ 1) and/or on their difference ( ⁇ SEL_i - ⁇ SEL _( i +1) ) at the acquisition instants of the two images A i and A i +1 .
- the method executes steps L, M, and N only if the following checking step gives a positive outcome:
- step K checking that the two images A i and A i +1 of the pair have been acquired in periods of the year in which vegetation is in a comparable blooming condition, while in the case step K gives a negative outcome the following step is executed:
- Such radiometers used I the method according to the invention have at least two Infra-red channels, among which a NIR channel and a SWIR or MIR channel.
- the method discriminates: the set BA of resolution cells or pixels p xy of the image A i +1 , identified by coordinates ( xy ), corresponding to zones which have been subject to fire during the time interval between the acquisition of image A i and the acquisition of image A i +1 ; the set DB of resolution cells p xy of the image A i +1 corresponding to zones which have been subject to a decrease in biomass content or to an increase in soil or mineral cover during the time interval between the acquisitions of the two successive images A i and A i +1 ; and the set DV of resolution cells p xy of the image A i +1 corresponding to zones which have been subject to a loss of vegetation during the time interval between the acquisitions of the two successive images A i and A i +1 .
- the method identifies the three sets BA, DB, and DV in relation to pseudo invariant features or Permanent Reflectors, also indicated hereinafter as "PRs", by performing the steps illustrated in the following with reference to Figure 1.
- radiance L NIR of the NIR channel and radiance L SWIR of the SWIR channel are transformed, respectively, into reflectance p NIR and into reflectance ⁇ SWIR (and/or in reflectance ⁇ MIR ), respectively, preferably taking account of the exo-atmospheric solar irradiance E 0 and of the effects of the sun zenith angle ⁇ z , and estimating the path radiance L p , more preferably through the DOS method, still more preferably by employing equation [1].
- the resolution cells p xy of the N images A i are assumed as PRs.
- Such difference gives an estimation of the maximum difference of the cumulative contribution in atmospheric transmittance T v from the target towards the sensor, in atmospheric transmittance T z in the illumination direction, in optical thickness t R for Rayleigh scattering, and atmospheric aerosol effects E A .
- the percentage P xy of variation, from the i -th image A i to the time successive ( i +1)-th image A i +1 , of reflectances ⁇ specifically: the percentage P xy_NIR _ i of variation of the reflectance ⁇ NIR of pixel p xy in NIR channel, and the percentage P xy_NIR / SWIR_i (and/or P xy_NIR / MIR_i ) of variation of the ratio ⁇ NIR / ⁇ SWIR (and/or of the ratio ⁇ NIR / ⁇ MIR ) of the reflectances ⁇ NIR and ⁇ SWIR (and/or ⁇ MIR ) of pixel p xy in the NIR/SWIR channel ratio (and/or in the NIR/MIR channel ratio).
- step G it is determined a first set S 1 of resolution cells p xy of the image A i +1 corresponding to zones which have been subject to a decrease in biomass content or an increase in soil or mineral cover during the time interval between the acquisitions of the two successive images A i and A i+ 1 .
- the first set S 1 comprises the resolution cells p xy of the image A i +1 for which the percentages calculated in the sixth step F are lower than the difference ( M PR_Pi - ⁇ Pi ) calculated in the preceding step D in the corresponding channels of the response of the Permanent Reflectors PRs, calculated in the third step C.
- the first set S 1 comprises the pixels p xy for which the following inequalities are satisfied: P x ⁇ y - ⁇ N ⁇ I ⁇ R - ⁇ i ⁇ M P ⁇ R - ⁇ P - ⁇ N ⁇ I ⁇ R - ⁇ i - ⁇ P ⁇ R - ⁇ P - ⁇ N ⁇ I ⁇ R - ⁇ i , and P xy - NIR / SWIR - ⁇ i ⁇ M PR - ⁇ P - ⁇ NIR / SWIR - ⁇ i - ⁇ P ⁇ R - ⁇ P ⁇ ⁇ NIR / SWIR - ⁇ i and / or P xy - ⁇ NIR / MIR - ⁇ i ⁇ M PR - ⁇ P - ⁇ NIR / MIR - ⁇ i - ⁇ P ⁇ R - ⁇ P ⁇ NIR / SWIR -
- step H it is determined a second set S 2 of resolution cells p xy of the image A i +1 corresponding to zones which have been subject to a loss of vegetation during the time interval between the acquisitions of the two successive images A i and A i +1 .
- the second set S 2 comprises the pixels p xy for which the following inequalities are satisfied: ⁇ xy - NIR - ⁇ i + 1 ⁇ M PR - ⁇ ⁇ NIR - ⁇ i + 1 ⁇ xy - NIR - ⁇ SWIR - ⁇ i + 1 ⁇ M PR - ⁇ ⁇ NIR - ⁇ ⁇ SWIR - ⁇ i + 1 and / or ⁇ xy - ⁇ NIR - ⁇ MIR - ⁇ i + 1 ⁇ M PR - ⁇ ⁇ NIR - ⁇ ⁇ MIR - ⁇ i + 1 .
- step I it is determined a third set S 3 of resolution cells p xy of the image A i +1 corresponding to zones which have been subject to a damage during the time interval between the acquisitions of the two successive images A i and A i +1 .
- the set DB of the resolution cells p xy of the image A i +1 corresponding to zones which have been subject to a decrease in biomass content or to an increase in soil or mineral cover during the time interval between the acquisitions of the two successive images A i and A i +1 .
- a check is carried out for ascertaining whether the solar elevation angles ⁇ SEL_i and ⁇ SEL_ ( i+ 1) at the acquisition instants of the two images A i and A i +1 are both higher or both lower than a threshold value ⁇ .
- Such check assures that the two images A i and A i +1 have been acquired in periods of the year in which vegetation is in a comparable blooming condition, in order to grant reliability of mapping.
- Such check may be easily carried out by checking the following inequality: ⁇ SEL - ⁇ i - ⁇ * ⁇ SEL - ⁇ i + 1 - ⁇ > 0
- the method may comprise the use of corrective factors in calculations for determining the first set S 1 and/or the second set S 2 , corrective factors preferably depending on the solar elevation angles ⁇ SEL_i and ⁇ SEL_ ( i +1) , and possibly on their difference ( ⁇ SEL_i - ⁇ SEL_ ( i +1) ), at the acquisition instants of the two images A i and A i +1 .
- step K gives a positive outcome
- the three steps L, M and N are performed; otherwise, if the check of step K gives a negative outcome, step P is performed.
- step L it is determined, within the third set S 3 , a fourth set S 4 ⁇ S 3 of resolution cells p xy corresponding to zones wherein values of reflectance ⁇ xy of the preceding image A i identify them as cells corresponding to areas provided with "burnable" vegetation, such as forestry, agro-forestry, and semi-natural areas, and zones of permanent crops (preferably no zone is subject to human changes).
- step M it is determined the set BA of the resolution cells p xy of the image A i +1 corresponding to zones which have been subject to a fire during the time interval between the acquisition of the image A i and the acquisition of the image A i +1 .
- step N it is determined the set DV of the resolution cells p xy of the image A i +1 corresponding to zones which have been subject to a loss of vegetation during the time interval between the acquisitions of the two successive images A i and A i +1 .
- step P is performed, with which it is determined only the set DV of the resolution cells p xy of the image A i +1 corresponding to zones which have been subject to a loss of vegetation during the time interval between the acquisitions of the two successive images A i and A i +1 .
- the method comprises a further final step wherein the pixels p xy belonging to the set BA which are singly isolated from the other pixels of the set BA (or which belong to groups of pixels p xy of the set BA creating a zone having a reduced surface extension) are re-assigned to the set to which the surrounding pixels belong.
- This is done in order to reduce wrong assignments due to occurrence of limit values of reflectance or of the related statistical parameters used in calculations performed by the method and which may introduce evaluation ambiguities.
- this final step of the method may re-assign these pixels to the set BA. Similar processings may be carried out for the sets DB and DV.
- the method according to the invention enables estimating biomass content and vegetation cover from one date to another and ascertaining presence of areas which have been subject to loss of vegetation or fires during the same period.
- FIG 2 shows the so-called spectral reflectance signature ⁇ xy of 12 zones pertaining to different vegetation classes for two acquisition dates and the spectral response of the same zone after having burned.
- Data have been collected through satellite platform Landsat-5 TM.
- Landsat 5 TM channels number 1-5, and 7 are indicated as 1-6 in X coordinate, while the corresponding reflectance percentage is reported in Y coordinate.
- the wavelength interval of channels 1-5, and 7 are respectively 0.45-0.52 (Visible, or VIS, Blue), 0.53-0.61 (VIS-Green), 0.63-0.69 (VIS-Red), 0.78-0.90 (NIR), 1.55-1.75 (SWIR) and 2.09-2.35 (SWIR)
- the burned areas are obtained by a visual comparison of the September 12, 1997, and the August 30, 1998, Landsat 5 TM images, made by an operator, by using various RGB combination in order to visually enhance the vegetation alteration.
- Such areas which are sites of fire events reported in National Forest Guards and/or National Fire Brigades databases, are selected as unequivocally identified shapes of fire-damaged surfaces.
- Vegetation classes are indicated as land use Corine classes (according to the nomenclature of the known CORINE European project of coordinated information on the land use in the European context: COoRdinated INformation on the European Environment).
- the vegetation classes presenting a minimum of three distinct burned areas with a surface greater than 3 hectares are considered as fire-damaged areas.
- Resolution cells p xy situated in the inner central part of areas identified as fire-damaged have been selected and their associated reflectance percentages in channels 1-5, and 7 have been adopted as characteristic spectral signature of burnt area for the relevant vegetation class.
- spectral signatures of vegetated resolution cells on the two Landsat 5 images acquired on two different dates are different. This means that a quantitative evaluation of the damages occurred from one date to another by using multi-temporal analysis of remote sensing data can be performed only if the vegetation variation is calculated in relationship to a reference. Intrinsic reflectance spectral signature of this reference target homogeneoc!sly distributed over the whole remote sensing image must be independent from the acquisition date of the remote sensing data. Such reference targets are just the Permanents Reflectors PRs. The analysis of the vegetation variation from one date to another must be calculated by considering the distance from the vegetation spectral variation to the PRs spectral variation.
- the detection change from one date to another is decomposed in two components: one as representative of the variations due to difference of sensor sensitivity and radiometric range, to difference of wavelength centre value, to difference of up-welling and down-welling transmittance along the ground sensor path between the two image acquisition dates, and the second component as representative of the effective variation of the vegetation state between the two remote sensed image acquisition dates.
- the reflectance variation from vegetation to burned area is strictly related to the vegetation type and biomass content originally present, whereas the spectral signature of fire-damaged areas, on the post-fire image, appears slightly depending on the initial vegetation cover.
- the NIR reflectance decrease associated to fire damage is particularly significant for vegetal land classes characterised by the presence of trees and/or shrubs, such as: vineyards, fruit trees and olive groves (respectively identified as classes 221, 222, and 223 in the CORINE classification) or agro-forestry areas, transitional woodland and shrubs (respectively identified as classes 244, 323 and 324).
- the poorly defined class 242, including all complex cultivation patterns, is also characterised by a notable NIR reflectance decrease.
- Figure 3 shows the spectral reflectance signature in some Corine vegetation classes and the response from densely urbanised areas, always obtained in channels 1-5, and 7 of the satellite platform Landsat 5 TM radiometer. Said urban areas homogeneously distributed over the whole sensed image are considered as PRs.
- channels 1-5 and 7 of Landsat-5 TM are numbered as 1 to 6 in X coordinate.
- NIR reflectance percentage of PRs and reflectance percentage ratio in NIR/SWIR channels of PRs are lower than the corresponding values of the vegetation classes taken under consideration.
- Figure 4 shows some spectral reflectance signatures of the fire-damaged areas corresponding to the preceding Corine classes, considered in Figure 3, and the spectral response of densely urbanised areas PR areas, still obtained in channels 1-5, 7 (indicated with numbers 1 to 6 in X coordinate) through the satellite platform Landsat 5 TM.
- Figure 4 shows that fire-damaged areas have the same general form of spectra independently from the original vegetation cover class. Moreover, both NIR reflectance percentage of PRs and reflectance percentage ratio in NIR/SWIR channels of PRs are higher than the corresponding values of all fire-damaged areas, independently from the original vegetation cover class.
- Figure 5 represents a false colour image of Mount Etna acquired by means of the satellite platform SPOT-4 in channels NIRSWIR-RED on 7th June 2001.
- SPOT-4 borne multi spectral radiometer is characterised by four channels in the following wavelength intervals: 0.50-0.59 (VIS-Green), 0.61-0.68 (VIS-Red), 0.78-0.89 (NIR) and 1.58-1.75 (SWIR). Pixels size of SPOT images is equal to 20 meters* 20 meters.
- Figure 6 represents a false colour composite image of Mount Etna acquired by means of the satellite platform Landsat-7 in channels SWIR-NIR-RED on 29th July 2001.
- Landsat-7 borne radiometer pixels originally having size equal to 30 meters*30 meters, have been resized to pixels of size equal to 20 meters*20 meters with re-sampling method that does not alter the original land radiative spectrum.
- the two remote sensed images are acquired by two different multi-spectral radiometers, i.e. SPOT-4 (radiometer HRVIR-1) and Landsat-7 (radiometer ETM+), which are furthermore differen in radiometric response distribution, different in channel central wavelength, different in radiometric range, and different in pixel size.
- SPOT-4 radioometer HRVIR-1
- Landsat-7 radiometric ETM+
- the two remote sensed images are acquired in different seasons since the pre-fire image has been acquired in late spring at the beginning of June, when the vegetation is characterised by a high level of biomass content and ground coverage, whereas the post-fire image has been acquired in summer at the end of July, when the vegetation is characterized by high hydric stress. Due to seasonal changes, all vegetation Corine classes are characterised by a substantial decrease in both the biomass content and the vegetation cover.
- Figure 7 represents the image comprising resolution cells, of size equal to 20 meters*20 meters, classified through variation analysis as damaged vegetation cells, which are characterised by a decrease in percentage of biomass content (decrease in percentage of NIR reflectance value) and in vegetation cover (decrease in percentage of NIR/SWIR ratio) higher than the PRs spectral variation calculated between the two remote sensed images of 7th June 2001, acquired by SPOT-4, and 29th July 2001, acquired by Landsat 7 ETM+.
- PRs are dense urban areas, the spectral variation of which is assumed as representative of differences in sensor sensitivity and radiometric range, in wavelength centre value, or in response histogram in the wavelength interval of NIR and SWIR channels, in up-welling and down-welling transmittance along the ground sensor path between the two remote sensing images of 7th June 2001, acquired by SPOT-4, and 29th July 2001, acquired by Landsat 7 ETM+.
- the North-West part of the image is mainly occupied by evergreen forests and it can be noted that for pixels classified as damaged resolution cells from 7th June to 29th July, the decrease in biomass content and vegetation cover is principally due to seasonal changes occurring between spring and summer vegetation states.
- Figure 8 shows all the de-vegetated resolution cells present in the Landsat 7 ETM+ image acquired on 29th July 2001.
- Said resolution cells are characterised by biomass content (given by the NIR reflectance value) and vegetation cover (given by NIR/SWIR ratio) lower than the one of PRs (dense urban areas) on final remote sensing data Landsat 7 ETM+ of 29th July.
- PRs dense urban areas
- Figure 9 shows all the resolution cells wherein a decrease in biomass content and vegetation cover occurs from 7th June 2001 to 29th July 2001 and which are classified as de-vegetated resolution cells on 29th July 2001.
- the pixel clusters constituted by more than 10 pixels selected by both multi-temporal variation ( Figure 7) and post fire analysis ( Figure 8) are classified as fire damaged areas.
- Figure 10 represents the false colour composite image of 29th July 2001 acquired by Landsat-7 ETM + in channels SWIR-NIR-Blue overlaying the layer GIS of fire damaged surfaces larger than 1 hectare, defined in Figure 9.
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Image Processing (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Mobile Radio Communication Systems (AREA)
- Stereo-Broadcasting Methods (AREA)
Abstract
Description
- The present invention refers to an automatic method of detecting and mapping, particularly for burnt areas without vegetation, enabling, in a precise, simple, versatile, reliable, and efficient way, the discrimination of areas, in images acquired from satellite, where fires occurred on the basis of a comparison of images overcoming possible problems of dishomogeneity arising in the case when the conditions of acquisition of the images to be compared are different. More in particular, the automatic mapping method is based on the multi-temporal analysis of multi spectral data, or on the changes in the spectral response of the Earth features in combined near infrared, or NIR (Near Infra-Red), and short wave infrared, or SWIR (Short Wave Infra-Red), or medium infrared, or MIR (Medium Infra-Red), windows of the electromagnetic spectrum.
- The present invention further refers to an apparatus performing such a method, which is simple, fast, and inexpensive.
- A spaceborne hot temperature event detection arrangement is described in document EP-A-1 048 928.
- Multi spectral radiometers are generally constituted by one or more passive optical sensors that register the radiance emitted by the Earth in various windows of the electromagnetic spectrum.
- The principal windows of the electromagnetic spectrum commonly used for land and vegetation analysis are the Visible (with wavelength from 0.4 µm to 0.7 µm, hereinafter referred to as VIS), the Near Infra-Red (with wavelength from 0.7 µm to 1.3 µm, hereinafter referred to as NIR), the Short Wavelength Infra-Red (with wavelength from 1.3 µm to 2.5 µm, hereinafter referred to as SWIR), the Medium Infra-Red (with wavelength from 2.5 µm to 4.5 µm, hereinafter referred to as MIR), and the Thermal Infra-Red (with wavelength up to 14.5 µm, hereinafter referred to as TIR).
- Radiometers operating at these wavelengths are installed on airborne platforms - such as microlights, airplanes, balloons and dirigibles, operating between the ground free surface and the stratosphere - and spaceborne platforms, including all type of orbiting and sub-orbital spacecrafts at altitudes up to about 36.000 Km above the Earth's ground free surface.
- All the afore-said systems look at the Earth orthogonally to the direction of the platform's motion with an angle ranging from 0 to 35 degrees with respect to Nadir - that is, perpendicularly to the Earth. The aforementioned angle is called "Off Nadir".
- With such system, resolution cells or grids of the Earth surface can be generated with a spatial resolution (that is, the spacing of the elementary area, or pixel, within which it is possible to distinguish two objects) ranging from few centimetres to several kilometres.
- The main characteristic of multi spectral radiometers is that they are coherent imaging systems. This means that they allow measuring radiated energy of each Earth surface cell at different wavelengths, as well as measuring the difference in radiated energy between two independent cells in each wavelength.
- In theory, converting radiance in reflectance allows measuring the time variation in energy at each wavelength and in each Earth surface cell, and measuring the time variation of energy difference between two independent cells at each wavelength.
-
- In equation [1], Lp is the path radiance (that is the additional radiance due to the path in the atmosphere and hence to the diffusion or scattering phenomenon), Ed is the downward irradiance, E0 is exo-atmospheric solar irradiance, and θz is the solar zenith angle. The atmospheric transmittance Tv from the target toward the sensor depends on the off-nadir angle. The atmospheric transmittance Tz in the illumination direction is a function of θz, and of the optical thickness to be accounted for in Rayleigh scattering and atmospheric aerosol effects. Finally, d is the radius vector in astronomic units, given by
- Radiance Lp can be evaluated by using the so-called Dark Object Subtraction or DOS method, which, as shown by P.S. Chavez in "An improved Dark-Object Subtraction Technique for Atmospheric Scattering Correction of Multispectral Data", Remote Sensing of the Environment 24, pp. 459-479, 1988, relies upon the hypothesis that some ground types - as some types of rocks and soils - behave as almost perfect "black bodies": in such cases, the associated reflectance measured at the satellite is almost entirely due to scattered path reflectance. Atmospheric radiative transfer codes (Lowtran, Modtran, Code-6S, etc.) make use of standard atmosphere models, or of available measures of vertical atmosphere profiles, for determining the downward irradiance Ed. Finally, the accuracy of atmospheric corrections for Rayleigh scattering and aerosol effects strictly depend upon the availability of data on atmosphere physical and chemical parameters at the precise time of the remote sensing image acquisition
- In case of multi-temporal analysis, the fundamental key is to distinguish reflectance changes originated by illumination and atmospheric conditions, by sensor drift or by sensor differences, from reflectance changes due to decrease in biomass content, to chlorophyll absorption, to water content and to density in the vegetation cover. A method for separating changes due to intrinsic reflectance value of vegetation from changes due to other factors influencing the at-satellite reflectance, is that of identifying "pseudo-invariant" features to be used as "reference targets" in cross-rectification of reflectance in different images. Identification of "pseudo-invariant" features requires accurate analysis of the whole of remote sensing data involved in the multi temporal analysis, following a linear radiometric transformation and aimed to implement relative reflectance retrieval. Accuracy of radiometric rectification strongly depends on the selection of pseudo-invariant features or of reference targets.
- A large number of algorithms and processings intended for mapping burnt areas starting from multi spectral remote sensing data were carried out on single date image analysis, on the basis of so-called Principal Components, band ratio or vegetation index, unsupervised or supervised classification clustering. However, use of a single remote sensing image for mapping burnt and de-vegetated areas does not allow these ones to be distinguished from urban areas, water surfaces, sparsely vegetated areas and rock outcrops.
- Multi-temporal methods based on the evaluation of changes of the spectral response between NIR and SWIR channels seem advantageous, since fire is an agent of land cover alteration, and therefore change detection methods reduce the likelihood of confusion with non vegetated or low vegetated static land cover types, such as urban areas, and rock outcrops
- Multi temporal analysis were performed through Principal Component analysis followed by maximum likelihood classification, band ratio or vegetation index as NDVI index (Normalized Difference Vegetation Index) variations. The reflectance variation in areas passing from pre-fire situation to post-fire situation is strictly related to the vegetation type and the originally present biomass content. If only changes in reflectance response were considered, confusion may be then made between variations due to fire damage and variations due to hydric stress or seasonal vegetation changes.
- It is therefore an object of the present invention to provide an automatic method capable to significantly improve the unsteadiness, and hence the poor reliability, of present processing techniques, allowing detecting, measuring, and mapping burned and de-vegetated areas in a precise, simple, versatile, reliable, and efficient way.
- It is specific subject matter of this invention an automatic method of detecting and mapping, particularly for burnt areas without vegetation, in N images Ai, with N ≥ 2, for i = 1, 2,..., N, obtained from data acquired along a certain time period on the same geographical area by means of one or more spectral passive radiometers on one or more detection frequency channels, each image Ai comprising a same plurality of resolution cells or pixels pxy, the reflectance p of each one of said one or more channels being available for each image, characterised in that it comprises the following steps:
- B. for each pair of successive images Ai and A i+1, for i = 1, 2,..., N-1, assuming as Permanent Reflectors PRs the resolution cells pxy which have optical and electromagnetic properties constant in both images;
- C. for each pair of successive images Ai and A i+1, for i = 1, 2,..., N-1, calculating the mean MPR_P_i and the standard deviation σPR_P_i of one or more first quantities depending on the variation, from the i-th image Ai to the time successive (i+1)-th image A i+1, of one or more second quantities depending on at least one reflectance ρ of the Permanent Reflectors PRs in said one or more channels;
- D. calculating the difference (MPR_Pi - σPi ) between the mean MPR_P_i and the standard deviation σPR_P_i, calculated in step C, for each one of said one or more first quantities;
- E. calculating the mean values MPR_ρ of said one or more second quantities, depending on at least one reflectance ρ of the Permanent Reflectors PRs in said one or more channels, in each one of the second images A i+1 of the pairs of successive images (Ai and A i+1), for i =1,2,..., N-1;
- F. for each resolution cell pxy, calculating said one or more first quantities, already considered in step C for the Permanent Reflectors PRs, depending on the variation, from the i-th image Ai to the time successive (i+1)-th image A i+1, of said one or more second quantities depending on at least one reflectance ρ of the pixels pxy in said one or more channels;
- G. determining a first set S1 of resolution cells pxy of the image A i+1, comprising the resolution cells pxy of the image A i+1 for which said one or more first quantities calculated in step F are lower, in each one of said one or more channels, than the difference (MPR_Pi - σPi ) calculated in step D for each corresponding quantity of said one or more first quantities related to the Permanent Reflectors PRs;
- H. determining a second set S2 of resolution cells pxy of the image A i+1 comprising the resolution cells pxy of the image A i+1 (for i = 1, 2,..., N-1) for which said one or more second quantities, depending on at least one reflectance ρ of the resolution cell pxy in said one or more channels in each one of the second images A i+1 of the pairs of successive images, for i = 1, 2,..., N-1, are lower than all the mean values MPR_ρ of the corresponding one or more second quantities of the Permanent Reflectors PRs calculated in step E;
- I. determining a third set S3 of resolution cells pxy of the image A i+1 equal to the set resulting from the intersection of the first set S1 with the second set S2:
- J. determining a set DB of the resolution cells pxy of the image A i+1 classified as corresponding to zones which have been subject to a decrease in biomass content or to an increase in soil or mineral cover during the time interval between the acquisitions of the two successive images Ai and A i+1, the set DB being equal to the set of the pixels pxy belonging to the first set S1 which do not belong to the third set S3 :
- L. determining, within the third set S3, a fourth set S4 ⊆ S3 of resolution cells pxy corresponding to areas provided with "burnable" vegetation in the first images Ai of the pairs of successive images, for i = 1, 2,..., N-1;
- M. determining a set BA of the resolution cells pxy of the image A i+1 classified as corresponding to zones which have been subject to a fire during the time interval between the acquisition of the image Ai and the acquisition of the image A i+1, the set BA being equal to the fourth set S4:
- N. determining a set DV of the resolution cells pxy of the image A i+1 classified as corresponding to zones which have been subject to a loss of vegetation during the time interval between the acquisitions of the two successive images Ai and A i+1, the set DV being equal to the set of pixels pxy belonging to the second set S2 which do not belong to the fourth set S4 :
- Always according to the invention, said one or more detection frequency channels may be at least two infrared channels, preferably a near infrared or NIR (Near Infra Red) channel, and a short wave infrared or SWIR (Short Wave Infra Red) channel, and/or a medium infrared or MIR (Medium Infra Red) channel.
- Still according to the invention, the method may preliminarily process the data acquired by said one or more spectral passive radiometers for obtaining said N images Ai, with N ≥ 2, for i = 1, 2,..., N, of the same geographical area having the same plurality of resolution cells or pixels pxy and the same size thereof.
- Furthermore according to the invention, in step B Permanent Reflectors PRs may be determined as those resolution cells pxy which, in both images Ai and A i+1 of the pair, have values of reflectance ρ in each one of said one or more channels which correspond to densely built urban areas, preferably presenting no or little vegetation with respect to the size of the resolution cells.
- Always according to the invention, said one or more first quantities depending on a variation, considered in steps C, D, F, and G, may comprise at least one quantity identically equal to the corresponding variation.
- Preferably according to the invention, said one or more first quantities depending on a variation, considered in steps C, D, F, and G, comprise at least one percentage of the corresponding variation.
- Still according to the invention, said one or more second quantities depending on at least one reflectance ρ, considered in steps C, E, F, and H, may comprise at least one quantity identically equal to the corresponding reflectance ρ of a channel.
- Furthermore according to the invention, said one or more second quantities depending on at least one reflectance ρ, considered in steps C, E, F, and H, may comprise at least one ratio between two reflectances ρ of two channels.
- Always according to the invention, the acquired data may comprise the radiance L of each one of the N images Ai on each one of said one or more detection frequency channels, and the method, before executing step B, may transforms the radiance L into the reflectance ρ for each one of said one or more channels, preferably taking account of the exo-atmospheric solar irradiance E0 and of the effects of the sun zenith angle θz, and comprising the estimation of the path radiance Lp, more preferably through the DOS method.
- Preferably according to the invention, for determining the first set S1 in step G and/or the second set S2 in step H, the method takes vegetation seasonal differences into account.
- Still according to the invention, the method may take vegetation seasonal differences into account through the use of one or more corrective factors, preferably depending on the solar elevation angles θ SEL_i and θ SEL_(i+1) and/or on their difference (θ SEL_i - θ SEL_(i+1)) at the acquisition instants of the two images Ai and A i+1.
- Preferably according to the invention, the method executes steps L, M, and N only if the following checking step gives a positive outcome:
- K. checking that the two images Ai and A i+1 of the pair have been acquired in periods of the year in which vegetation is in a comparable blooming condition,
while in the case step K gives a negative outcome the following step is executed: -
- Furthermore according to the invention, step K may check whether the solar elevation angles θSEL_i and θ SEL_(i+1) at the acquisition instants of the two images Ai and A i+1 are both higher or both lower than a threshold value Θ, preferably equal to 35° (Θ= 35°).
- It is further subject matter of this invention an Electronic apparatus for processing data acquired along a certain time period on the same geographical area by means of one or more spectral passive radiometers on one or more detection frequency channels, for supplying a mapping, in particular for burnt areas without vegetation, of N images Ai, with N ≥ 2, for i = 1, 2,..., N, each image Ai comprising a same plurality of resolution cells or pixels pxy , the apparatus comprising:
- a data processing unit;
- at least one memory unit connected to said data processing unit; and
- at least one input/output, or I/O, interface connected to said data processing unit,
- the electronic apparatus being characterised in that it executes the previously illustrated automatic method of detecting and mapping, particularly for burnt areas without vegetation.
- The present invention will be now described, by way of illustration and not by way of limitation, according to its preferred embodiments, by particularly referring to the Figures of the enclosed drawings, in which:
- Figure 1 schematically shows a flow diagram of a preferred embodiment of the method according to the invention;
- Figure 2 shows the spectral reflectance signature ρxy of 12 zones belonging to different land cover vegetation classes for two acquisition dates and the spectral response of the same zone after fire;
- Figure 3 shows the spectral reflectance signature of some land cover classes of vegetation and the response of dense urban areas;
- Figure 4 shows the spectral reflectance signatures of the areas corresponding to Figure 3 after having been damaged by fire and the spectral response of dense urban PR areas;
- Figure 5 shows a first image of Mount Etna acquired by means of the multi spectral radiometer SPOT-4;
- Figure 6 shows a second image of Mount Etna acquired by means of the multi spectral radiometer Landsat-7 ETM+;
- Figure 7 shows a third image of Mount Etna obtained by means of the method of Figure 1 applied to the two images of Figures 5 and 6;
- Figure 8 shows a fourth image of Mount Etna obtained by means of the method of Figure 1 applied to the two images of Figures 5 and 6;
- Figure 9 shows a fifth image of Mount Etna obtained by means of the method of Figure 1 applied to the two images of Figures 5 and 6; and
- Figure 10 shows a sixth image of Mount Etna obtained by means of the method of Figure 1 applied to the two images of Figures 5 and 6.
- The method according to the invention performs a multi spectral measuring of burned and de-vegetated zones on the basis of N images Ai, with N ≥ 2, and for i = 1, 2,..., N, acquired over a certain time span by means of spectral passive radiometers. Such radiometers used I the method according to the invention have at least two Infra-red channels, among which a NIR channel and a SWIR or MIR channel.
process for multi spectral measuring of burned and de-vegetated zones which, with N ≥ 2 images acquired over an assigned time span, with multi spectral passive radiometers presenting at least two Infrared channels - one NIR and one SWIR or MIR - identifies for every resolution cell, the fire damaged ones in relationship to a pseudo invariant features or Permanent Reflectors PRs, characterized in that said fire damaged cells are identified through the following steps - In particular, for each pair of successive images Ai and A i+1, for i = 1, 2,..., N-1, the method discriminates:
the set BA of resolution cells or pixels pxy of the image A i+1, identified by coordinates (xy), corresponding to zones which have been subject to fire during the time interval between the acquisition of image Ai and the acquisition of image A i+1;
the set DB of resolution cells pxy of the image A i+1 corresponding to zones which have been subject to a decrease in biomass content or to an increase in soil or mineral cover during the time interval between the acquisitions of the two successive images Ai and A i+1; and
the set DV of resolution cells pxy of the image A i+1 corresponding to zones which have been subject to a loss of vegetation during the time interval between the acquisitions of the two successive images Ai and A i+1. - The method identifies the three sets BA, DB, and DV in relation to pseudo invariant features or Permanent Reflectors, also indicated hereinafter as "PRs", by performing the steps illustrated in the following with reference to Figure 1.
- In the first step A, for each cell pxy of each one of the N images Ai, radiance LNIR of the NIR channel and radiance LSWIR of the SWIR channel (and/or radiance LMIR of the MIR channel) are transformed, respectively, into reflectance pNIR and into reflectance ρSWIR (and/or in reflectance ρMIR ), respectively, preferably taking account of the exo-atmospheric solar irradiance E0 and of the effects of the sun zenith angle θz, and estimating the path radiance Lp, more preferably through the DOS method, still more preferably by employing equation [1].
- In the second step B, the resolution cells pxy of the N images Ai, the values of reflectance ρ of which identify them as cells corresponding to densely built urban areas, presenting no or little vegetation, are assumed as PRs.
- In the third step C, for each pair of successive images Ai and A i+1, for i = 1, 2,..., N-1, some statistical parameters of percentage P of variation, from i-th image Ai to the time successive (i+1)-th image A i+1, of reflectances ρ of the response of the Permanent Reflectors PRs are calculated. In particular, the following are calculated:
the mean MPR_P_NIR_i of percentage P of variation of the reflectance ρNIR in NIR channel of the response of the Permanent Reflectors PRs, for all the pairs of successive images Ai and A i+1, for i = 1,2,..., N-1 ;
the standard deviation σPR_P_NIR_i of percentage P of variation of the reflectance ρNIR in NIR channel of the response of the Permanent Reflectors PRs, for all the pairs of successive images Ai and A i+1, for i = 1, 2,..., N-1 ;
the mean M PR_P_NIR/SWIR_i (and/or the mean M PR_P_NIR/MIR_i ) of percentage P of variation of the ratio ρNIR / ρSWIR (and/or of the ratio ρNIR / ρMIR ) of the reflectances ρNIR and ρSWIR (and/or ρMIR ) in NIR/SWIR channel ratio (and/or in NIR/MIR channel ratio), for all the pairs of successive images Ai and A i+1, for i = 1, 2,..., N-1 ; and
the standard deviation σPR_P_NIR/SWIR_i (and/or the standard deviation σ PR_P_NIR/SWIR_i ) of percentage P of variation of the ratio ρNIR / ρSWIR (and/or of the ratio ρNIR / ρMIR ) of the reflectances ρNIR and ρWEIR (and/or ρ MIR ) in NIR/SWIR channel ratio (and/or in NIR/MIR channel ratio), for all the pairs of successive images Ai and A i+1, for i = 1, 2,..., N-1. - In the fourth step D, it is calculated the difference
between the values of the means MPR_Pi and the corresponding standard deviation σPi for the NIR channel and the NIR/SWIR channel ratio (and/or the NIR/MIR channel ratio), calculated in step C for each pair of successive images Ai and A i+1, for i = 1, 2,...,N-1. Such difference gives an estimation of the maximum difference of the cumulative contribution in atmospheric transmittance Tv from the target towards the sensor, in atmospheric transmittance Tz in the illumination direction, in optical thickness tR for Rayleigh scattering, and atmospheric aerosol effects EA. - In the fifth step E, the mean values MPR_ρ of reflectances ρ of the response of the Permanent Reflectors PRs in each one of the second images A i+1 of the pairs of successive images (Ai and A i+1), for i = 1, 2,..., N-1 are calculated. In particular, the following are calculated:
the mean M PR_ρNIR_(i+1) of the reflectance ρNIR in NIR channel of the response of the Permanent Reflectors PRs in all the images A i+1, for i = 1, 2,..., N-1; and
the mean M PR_ρNIR_ρSWIR_(i+1) (and/or M PR_ρNIR_ρMIR_(i+1) of the ratio ρNIR / ρSWIR (and/or of the ratio ρNIR / ρ MIR) of the reflectances ρNIR and ρSWIR (and/or ρMIR ) in NIR/SWIR channel ratio (and/or in NIR/MIR channel ratio) of the response of the Permanent Reflectors PRs in all the images A i+1, for i=1,2,...,N-1. - In the sixth step F, for each resolution cell pxy, it is calculated the percentage Pxy of variation, from the i-th image Ai to the time successive (i+1)-th image A i+1, of reflectances ρ, specifically: the percentage P xy_NIR_i of variation of the reflectance ρNIR of pixel pxy in NIR channel, and the percentage P xy_NIR/SWIR_i (and/or P xy_NIR/MIR_i ) of variation of the ratio ρNIR / ρSWIR (and/or of the ratio ρNIR / ρMIR ) of the reflectances ρNIR and ρSWIR (and/or ρMIR ) of pixel pxy in the NIR/SWIR channel ratio (and/or in the NIR/MIR channel ratio).
- In the seventh step G, it is determined a first set S1 of resolution cells pxy of the image A i+1 corresponding to zones which have been subject to a decrease in biomass content or an increase in soil or mineral cover during the time interval between the acquisitions of the two successive images Ai and A i+1 .
- In particular, the first set S1 comprises the resolution cells pxy of the image A i+1 for which the percentages calculated in the sixth step F are lower than the difference (MPR_Pi- σ Pi ) calculated in the preceding step D in the corresponding channels of the response of the Permanent Reflectors PRs, calculated in the third step C. That is, the first set S1 comprises the pixels pxy for which the following inequalities are satisfied:
- In the eighth step H, it is determined a second set S2 of resolution cells pxy of the image A i+1 corresponding to zones which have been subject to a loss of vegetation during the time interval between the acquisitions of the two successive images Ai and A i+1.
- In particular, the second set S2 comprises the resolution cells pxy of the image A i+1 (for i = 1, 2,..., N-1) for which the reflectances ρxy in the various channels are lower than the mean value MPR_ρ of reflectances p of the response of the Permanent Reflectors PRs in the corresponding channels of the corresponding images A i+1, calculated in the fifth step E. That is, the second set S2 comprises the pixels pxy for which the following inequalities are satisfied:
- In the subsequent ninth step I, it is determined a third set S3 of resolution cells pxy of the image A i+1 corresponding to zones which have been subject to a damage during the time interval between the acquisitions of the two successive images Ai and A i+1.
-
- In the tenth step J, it is determined the set DB of the resolution cells pxy of the image A i+1 corresponding to zones which have been subject to a decrease in biomass content or to an increase in soil or mineral cover during the time interval between the acquisitions of the two successive images Ai and A i+1. In particular, the set DB is equal to the set of the pixels pxy belonging to the first set S1 which do not belong to the third set S3, (i.e. which do not belong to the second set S2 ):
- In the eleventh step K, a check is carried out for ascertaining whether the solar elevation angles θSEL_i and θ SEL_(i+1) at the acquisition instants of the two images Ai and A i+1 are both higher or both lower than a threshold value Θ.
- Such check assures that the two images Ai and A i+1 have been acquired in periods of the year in which vegetation is in a comparable blooming condition, in order to grant reliability of mapping. By way of example, in Italy it may be deemed that the two images Ai and A i+1 have been both acquired in spring or summer if θSEL_i and θ SEL_(i+1) are both higher than Θ= 35°, or they have been both acquired in autumn or winter if θSEL_i and θ SEL_(i+1) are both lower than Θ = 35°. Such check may be easily carried out by checking the following inequality:
- Alternatively, in order to consider the vegetation seasonal differences, the method may comprise the use of corrective factors in calculations for determining the first set S1 and/or the second set S2, corrective factors preferably depending on the solar elevation angles θSEL_i and θ SEL_(i+1), and possibly on their difference (θSEL_i - θ SEL_(i+1)), at the acquisition instants of the two images Ai and A i+1.
- In the case when the check of step K gives a positive outcome, the three steps L, M and N are performed; otherwise, if the check of step K gives a negative outcome, step P is performed.
- In step L it is determined, within the third set S3, a fourth set S4 ⊆ S3 of resolution cells pxy corresponding to zones wherein values of reflectance ρxy of the preceding image Ai identify them as cells corresponding to areas provided with "burnable" vegetation, such as forestry, agro-forestry, and semi-natural areas, and zones of permanent crops (preferably no zone is subject to human changes).
- In step M, it is determined the set BA of the resolution cells pxy of the image A i+1 corresponding to zones which have been subject to a fire during the time interval between the acquisition of the image Ai and the acquisition of the image A i+1. In particular, the set BA is equal to the fourth set S4
- In step N, it is determined the set DV of the resolution cells pxy of the image A i+1 corresponding to zones which have been subject to a loss of vegetation during the time interval between the acquisitions of the two successive images Ai and A i+1. In particular, the set DV is equal to the set of pixels pxy belonging to the second set S2 which do not belong to the fourth set S4 :
- In the case when the check of step K gives a negative outcome, it is not possible to determine with sufficient reliability the set BA, and step P is performed, with which it is determined only the set DV of the resolution cells pxy of the image A i+1 corresponding to zones which have been subject to a loss of vegetation during the time interval between the acquisitions of the two successive images Ai and A i+1. In particular, the set DV is equal to the second set S2 :
- Preferably, the method comprises a further final step wherein the pixels pxy belonging to the set BA which are singly isolated from the other pixels of the set BA (or which belong to groups of pixels pxy of the set BA creating a zone having a reduced surface extension) are re-assigned to the set to which the surrounding pixels belong. This is done in order to reduce wrong assignments due to occurrence of limit values of reflectance or of the related statistical parameters used in calculations performed by the method and which may introduce evaluation ambiguities. Similarly, in the case when pixels not assigned to the set BA are present within wide areas of pixels of the set BA, this final step of the method may re-assign these pixels to the set BA. Similar processings may be carried out for the sets DB and DV.
- It is evident that the method according to the invention enables estimating biomass content and vegetation cover from one date to another and ascertaining presence of areas which have been subject to loss of vegetation or fires during the same period.
- Characteristics and advantages of the method according to the invention are illustrated in the following with reference to Figures 2-10, which show some applications thereof.
- Figure 2 shows the so-called spectral reflectance signature ρxy of 12 zones pertaining to different vegetation classes for two acquisition dates and the spectral response of the same zone after having burned. Data have been collected through satellite platform Landsat-5 TM. In Figure 2,
Landsat 5 TM channels number 1-5, and 7 are indicated as 1-6 in X coordinate, while the corresponding reflectance percentage is reported in Y coordinate. The wavelength interval of channels 1-5, and 7 are respectively 0.45-0.52 (Visible, or VIS, Blue), 0.53-0.61 (VIS-Green), 0.63-0.69 (VIS-Red), 0.78-0.90 (NIR), 1.55-1.75 (SWIR) and 2.09-2.35 (SWIR) - The burned areas are obtained by a visual comparison of the September 12, 1997, and the August 30, 1998,
Landsat 5 TM images, made by an operator, by using various RGB combination in order to visually enhance the vegetation alteration. Such areas, which are sites of fire events reported in National Forest Guards and/or National Fire Brigades databases, are selected as unequivocally identified shapes of fire-damaged surfaces. Vegetation classes are indicated as land use Corine classes (according to the nomenclature of the known CORINE European project of coordinated information on the land use in the European context: COoRdinated INformation on the european Environment). In particular, the vegetation classes presenting a minimum of three distinct burned areas with a surface greater than 3 hectares are considered as fire-damaged areas. Resolution cells pxy situated in the inner central part of areas identified as fire-damaged have been selected and their associated reflectance percentages in channels 1-5, and 7 have been adopted as characteristic spectral signature of burnt area for the relevant vegetation class. - For each land cover Corine class, spectral signatures of vegetated resolution cells on the two
Landsat 5 images acquired on two different dates are different. This means that a quantitative evaluation of the damages occurred from one date to another by using multi-temporal analysis of remote sensing data can be performed only if the vegetation variation is calculated in relationship to a reference. Intrinsic reflectance spectral signature of this reference target homogeneoc!sly distributed over the whole remote sensing image must be independent from the acquisition date of the remote sensing data. Such reference targets are just the Permanents Reflectors PRs. The analysis of the vegetation variation from one date to another must be calculated by considering the distance from the vegetation spectral variation to the PRs spectral variation. Thus, for each pixel the detection change from one date to another is decomposed in two components: one as representative of the variations due to difference of sensor sensitivity and radiometric range, to difference of wavelength centre value, to difference of up-welling and down-welling transmittance along the ground sensor path between the two image acquisition dates, and the second component as representative of the effective variation of the vegetation state between the two remote sensed image acquisition dates. - The reflectance variation from vegetation to burned area is strictly related to the vegetation type and biomass content originally present, whereas the spectral signature of fire-damaged areas, on the post-fire image, appears slightly depending on the initial vegetation cover.
- In the NIR window, reflectance is directly proportional to the amount of biomass. Since green vegetation is very reflective in the NIR channel, the fire process typically gives rise to a reflectance decrease.
- The NIR reflectance decrease associated to fire damage is particularly significant for vegetal land classes characterised by the presence of trees and/or shrubs, such as: vineyards, fruit trees and olive groves (respectively identified as
classes classes class 242, including all complex cultivation patterns, is also characterised by a notable NIR reflectance decrease. - For all vegetation classes that underwent this analysis, NIR reflectance sharply decreases with respect to SWIR reflectance after fire damage
- Figure 3 shows the spectral reflectance signature in some Corine vegetation classes and the response from densely urbanised areas, always obtained in channels 1-5, and 7 of the
satellite platform Landsat 5 TM radiometer. Said urban areas homogeneously distributed over the whole sensed image are considered as PRs. As in Figure 2, also in Figure 3 channels 1-5 and 7 of Landsat-5 TM are numbered as 1 to 6 in X coordinate. - It can be noted in Figure 3 that, except for the case of
Corine class 333, corresponding to sparsely vegetated areas, and ofCorine class 223, corresponding to olive groves, NIR reflectance percentage of PRs and reflectance percentage ratio in NIR/SWIR channels of PRs are lower than the corresponding values of the vegetation classes taken under consideration. - Figure 4. shows some spectral reflectance signatures of the fire-damaged areas corresponding to the preceding Corine classes, considered in Figure 3, and the spectral response of densely urbanised areas PR areas, still obtained in channels 1-5, 7 (indicated with
numbers 1 to 6 in X coordinate) through thesatellite platform Landsat 5 TM. - Figure 4 shows that fire-damaged areas have the same general form of spectra independently from the original vegetation cover class. Moreover, both NIR reflectance percentage of PRs and reflectance percentage ratio in NIR/SWIR channels of PRs are higher than the corresponding values of all fire-damaged areas, independently from the original vegetation cover class.
- Figure 5 represents a false colour image of Mount Etna acquired by means of the satellite platform SPOT-4 in channels NIRSWIR-RED on 7th June 2001. SPOT-4 borne multi spectral radiometer is characterised by four channels in the following wavelength intervals: 0.50-0.59 (VIS-Green), 0.61-0.68 (VIS-Red), 0.78-0.89 (NIR) and 1.58-1.75 (SWIR). Pixels size of SPOT images is equal to 20 meters* 20 meters.
- Figure 6 represents a false colour composite image of Mount Etna acquired by means of the satellite platform Landsat-7 in channels SWIR-NIR-RED on 29th July 2001. Landsat-7 borne radiometer pixels, originally having size equal to 30 meters*30 meters, have been resized to pixels of size equal to 20 meters*20 meters with re-sampling method that does not alter the original land radiative spectrum.
- The subsequent Figures 7-10 show the results of an application of the method according to the invention in critical conditions with respect to better processing procedures for the method, according to the following reasons.
- First of all, the two remote sensed images are acquired by two different multi-spectral radiometers, i.e. SPOT-4 (radiometer HRVIR-1) and Landsat-7 (radiometer ETM+), which are furthermore differen in radiometric response distribution, different in channel central wavelength, different in radiometric range, and different in pixel size.
- Moreover, the two remote sensed images are acquired in different seasons since the pre-fire image has been acquired in late spring at the beginning of June, when the vegetation is characterised by a high level of biomass content and ground coverage, whereas the post-fire image has been acquired in summer at the end of July, when the vegetation is characterized by high hydric stress. Due to seasonal changes, all vegetation Corine classes are characterised by a substantial decrease in both the biomass content and the vegetation cover.
- Finally, on the year 2001 Etna volcano was the site of eruptive activity during summer and particularly at the end of July. 2001 eruptive activity of Etna was characterised by high level of ash fall, that is a strong catalyst for vegetation cover alteration, since it contributes to the increase of mineral cover on the ground and to the local drop of chlorophyll transformation of biomass.
- Figure 7 represents the image comprising resolution cells, of size equal to 20 meters*20 meters, classified through variation analysis as damaged vegetation cells, which are characterised by a decrease in percentage of biomass content (decrease in percentage of NIR reflectance value) and in vegetation cover (decrease in percentage of NIR/SWIR ratio) higher than the PRs spectral variation calculated between the two remote sensed images of 7th June 2001, acquired by SPOT-4, and 29th July 2001, acquired by
Landsat 7 ETM+. In the case under consideration, PRs are dense urban areas, the spectral variation of which is assumed as representative of differences in sensor sensitivity and radiometric range, in wavelength centre value, or in response histogram in the wavelength interval of NIR and SWIR channels, in up-welling and down-welling transmittance along the ground sensor path between the two remote sensing images of 7th June 2001, acquired by SPOT-4, and 29th July 2001, acquired byLandsat 7 ETM+. - The North-West part of the image is mainly occupied by evergreen forests and it can be noted that for pixels classified as damaged resolution cells from 7th June to 29th July, the decrease in biomass content and vegetation cover is principally due to seasonal changes occurring between spring and summer vegetation states.
- Differently, for pixels classified as damaged resolution cells from 7th June to 29th July in the Central East part of the image, the decrease in biomass content and vegetation cover is principally due to ash fall occurred during the eruption period of Etna in 2001.
- Figure 8 shows all the de-vegetated resolution cells present in the
Landsat 7 ETM+ image acquired on 29th July 2001. Said resolution cells are characterised by biomass content (given by the NIR reflectance value) and vegetation cover (given by NIR/SWIR ratio) lower than the one of PRs (dense urban areas) on final remotesensing data Landsat 7 ETM+ of 29th July. On this image, cold lava flows, eroded soils, rock areas and few urban sites are classified as de-vegetated resolution cells. - Figure 9 shows all the resolution cells wherein a decrease in biomass content and vegetation cover occurs from 7th June 2001 to 29th July 2001 and which are classified as de-vegetated resolution cells on 29th July 2001. In particular, in Figure 9 the pixel clusters constituted by more than 10 pixels selected by both multi-temporal variation (Figure 7) and post fire analysis (Figure 8) are classified as fire damaged areas.
- Figure 10 represents the false colour composite image of 29th July 2001 acquired by Landsat-7 ETM + in channels SWIR-NIR-Blue overlaying the layer GIS of fire damaged surfaces larger than 1 hectare, defined in Figure 9.
- The preferred embodiments have been above described and some modifications of this invention have been suggested, but it should be understood that those skilled in the art can make variations and changes, without so departing from the related scope of protection, as defined by the following claims.
Claims (18)
- Automatic method of detecting and mapping, particularly for burnt areas without vegetation, in N images Ai, with N ≥ 2, for i = 1, 2,..., N, obtained from data acquired along a certain time period on the same geographical area by means of one or more spectral passive radiometers on one or more detection frequency channels, each image Ai comprising a same plurality of resolution cells or pixels pxy, the reflectance p of each one of said one or more channels being available for each image, characterised in that it comprises the following steps:B. for each pair of successive images Ai and A i+1, for i = 1, 2,..., N-1, assuming as Permanent Reflectors PRs the resolution cells pxy which have optical and electromagnetic properties constant in both images;C. for each pair of successive images Ai and A i+1, for i = 1, 2,..., N-1, calculating the mean MPR_P_i and the standard deviation σPR_P_i of one or more first quantities (P) depending on the variation, from the i-th image Ai to the time successive (i+1)th image A i+1, of one or more second quantities (ρ NIR , ρNIR / ρSWIR, ρNIR / ρMIR ) depending on at least one reflectance ρ of the Permanent Reflectors PRs in said one or more channels;D. calculating the difference (MPR_Pi - σPi ) between the mean MPR_P_i and the standard deviation σPR_P_i, calculated in step C, for each one of said one or more first quantities (P);E. calculating the mean values MPR_ρ of said one or more second quantities (ρNIR , ρNIR /ρSWIR, ρNIR / ρMIR ), depending on at least one reflectance ρ of the Permanent Reflectors PRs in said one or more channels, in each one of the second images A i+1 of the pairs of successive images (Ai and A i+1), for i = 1, 2,...,N-1 ;F. for each resolution cell pxy, calculating said one or more first quantities (Pxy ), already considered in step C for the Permanent Reflectors PRs, depending on the variation, from the i-th image Ai to the time successive (i+1)-th image A i+1, of said one or more second quantities (ρNIR , ρNIR / ρSWIR, ρNIR / ρMIR) depending on at least one reflectance p of the pixels pxy in said one or more channels;G. determining a first set S1 of resolution cells pxy of the image A i+1, comprising the resolution cells pxy of the image A i+1 for which said one or more first quantities (Pxy ) calculated in step F are lower, in each one of said one or more channels, than the difference (MPR_Pi - σPi ) calculated in step D for each corresponding quantity (P) of said one or more first quantities (P) related to the Permanent Reflectors PRs;H. determining a second set S2 of resolution cells pxy of the image A i+1 comprising the resolution cells pxy of the image A i+1 (for i = 1, 2,..., N-1) for which said one or more second quantities (ρNIR, ρNIR / ρSWIR, ρNIR /ρMIR), depending on at least one reflectance ρ of the resolution cell pxy in said one or more channels in each one of the second images A i+1 of the pairs of successive images, for i = 1, 2,..., N-1, are lower than all the mean values MPR_p of the corresponding one or more second quantities (ρNIR, ρNIR / ρSWIR, ρNIR /ρMIR) of the Permanent Reflectors PRs calculated in step E;I. determining a third set S3 of resolution cells pxy of the image A i+1 equal to the set resulting from the intersection of the first set S1 with the second set S 2 :J. determining a set DB of the resolution cells pxy of the image A i+1 classified as corresponding to zones which have been subject to a decrease in biomass content or to an increase in soil or mineral cover during the time interval between the acquisitions of the two successive images Ai and A i+1, the set DB being equal to the set of the pixels pxy belonging to the first set S1 which do not belong to the third set S3 :L. determining, within the third set S3, a fourth set S4 ⊆ S3 of resolution cells pxy corresponding to areas provided with "burnable" vegetation in the first images Ai of the pairs of successive images, for i = 1, 2,..., N-1;M. determining a set BA of the resolution cells pxy of the image A i+1 classified as corresponding to zones which have been subject to a fire during the time interval between the acquisition of the image Ai and the acquisition of the image A i+1, the set BA being equal to the fourth set S4:N. determining a set DV of the resolution cells pxy of the image A i+1 classified as corresponding to zones which have been subject to a loss of vegetation during the time interval between the acquisitions of the two successive images Ai and A i+1, the set DV being equal to the set of pixels pxy belonging to the second set S2 which do not belong to the fourth set S4 :
- Method according to claim 1, characterised in that said one or more detection frequency channels are at least two infrared channels.
- Method according to claim 2, characterised in that said at least two infrared channels are a near infrared or NIR (Near Infra Red) channel, and a short wave infrared or SWIR (Short Wave Infra Red) channel, and/or a medium infrared or MIR (Medium Infra Red) channel.
- Method according to any one of the preceding claims, characterised in that the method preliminarily processes the data acquired by said one or more spectral passive radiometers for obtaining said N images Ai, with N ≥ 2, for i = 1, 2,..., N , of the same geographical area having the same plurality of resolution cells or pixels pxy and the same size thereof.
- Method according to any one of the preceding claims, characterised in that in step B Permanent Reflectors PRs are determined as those resolution cells Pxy which, in both images Ai and A i+1 of the pair, have values of reflectance ρ in each one of said one or more channels which correspond to densely built urban areas, preferably presenting no or little vegetation with respect to the size of the resolution cells.
- Method according to any one of the preceding claims, characterised in that said one or more first quantities (P) depending on a variation, considered in steps C, D, F, and G, comprise at least one quantity identically equal to the corresponding variation.
- Method according to any one of the preceding claims, characterised in that said one or more first quantities (P) depending on a variation, considered in steps C, D, F, and G, comprise at least one percentage (P) of the corresponding variation.
- Method according to any one of the preceding claims, characterised in that said one or more second quantities (ρNIR, ρNIR / ρSWIR, ρNIR / ρMIR) depending on at least one reflectance ρ, considered in steps C, E, F, and H, comprise at least one quantity identically equal to the corresponding reflectance ρ (ρNIR ) of a channel.
- Method according to any one of the preceding claims, characterised in that said one or more second quantities (ρNIR, ρNIR / ρSWIR, ρNIR /ρMIR) depending on at least one reflectance ρ, considered in steps C, E, F, and H, comprise at least one ratio between two reflectances p (ρNIR / ρSWIR, ρNIR / ρMIR) of two channels.
- Method according to any one of the preceding claims, characterised in that the acquired data comprise the radiance L of each one of the N images Ai on each one of said one or more detection frequency channels, and in that the method, before executing step B, transforms ([1]) the radiance L into the reflectance ρ for each one of said one or more channels.
- Method according to claim 10, characterised in that the transformation of the radiance L into the reflectance ρ takes account of the exo-atmospheric solar irradiance E0 and of the effects of the sun zenith angle θz, and it comprises the estimation of the path radiance Lp, preferably through the DOS method.
- Method according to any one of the preceding claims, characterised in that, for determining the first set S1 in step G and/or the second set S2 in step H, the method takes vegetation seasonal differences into account.
- Method according to claim 12, characterised in that the method takes vegetation seasonal differences into account through the use of one or more corrective factors.
- Method according to claim 13, characterised in that said one or more corrective factors depends on the solar elevation angles θSEL_i and θ SEL_(i+1) and/or on their difference (θSEL_i - θ SEL_(i+1)) at the acquisition instants of the two images Ai and A i+1.
- Method according to any one of claims 12 to 14, characterised in that it executes steps L, M, and N only if the following checking step gives a positive outcome:K. checking that the two images Ai and A i+1 of the pair have been acquired in periods of the year in which vegetation is in a comparable blooming condition,
while in the case step K gives a negative outcome the following step is executed: - Method according to claim 15, characterised in that step K checks whether the solar elevation angles θSEL_i and θ SEL_(i+1) at the acquisition instants of the two images Ai and A i+1 are both higher or both lower than a threshold value Θ.
- Method according to claim 16, characterised in that said threshold value Θ is equal to 35° (Θ= 35°).
- Electronic apparatus for processing data acquired along a certain time period on the same geographical area by means of one or more spectral passive radiometers on one or more detection frequency channels, for supplying a mapping, in particular for burnt areas without vegetation, of N images Ai, with N ≥ 2, for i = 1, 2,..., N , each image Ai comprising a same plurality of resolution cells or pixels pxy, the apparatus comprising:a data processing unit;at least one memory unit connected to said data processing unit; andat least one input/output, or I/O, interface connected to said data processing unit,the electronic apparatus being characterised in that it executes the automatic method of detecting and mapping, particularly for burnt areas without vegetation, according to any one of claims 1 to 17.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
IT000336A ITRM20030336A1 (en) | 2003-07-09 | 2003-07-09 | AUTOMATIC DETECTION AND MAPPING METHOD, |
PCT/IT2004/000376 WO2005005926A1 (en) | 2003-07-09 | 2004-07-07 | Method and apparatus for automatically detecting and mapping, particularly for burnt areas without vegetation |
Publications (2)
Publication Number | Publication Date |
---|---|
EP1642087A1 EP1642087A1 (en) | 2006-04-05 |
EP1642087B1 true EP1642087B1 (en) | 2006-12-27 |
Family
ID=29765907
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
EP04745201A Active EP1642087B1 (en) | 2003-07-09 | 2004-07-07 | Method and apparatus for automatically detecting and mapping, particularly for burnt areas without vegetation |
Country Status (5)
Country | Link |
---|---|
EP (1) | EP1642087B1 (en) |
AT (1) | ATE349677T1 (en) |
DE (1) | DE602004003944D1 (en) |
IT (1) | ITRM20030336A1 (en) |
WO (1) | WO2005005926A1 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107368617A (en) * | 2016-05-13 | 2017-11-21 | 南京理工大学 | Air-ground detection infrared imaging system operating distance computational methods based on Lowtran7 air softwares |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102074093B (en) * | 2010-12-29 | 2012-05-23 | 神华集团有限责任公司 | Satellite remote sensing-based mine fire monitoring method |
CN103700092B (en) * | 2013-12-04 | 2016-06-29 | 中国科学院遥感与数字地球研究所 | A kind of forest brulee extraction method based on sequential remote sensing image |
CN104318707B (en) * | 2014-11-04 | 2016-12-07 | 河南城建学院 | Forest fire detection platform based on satellite remote sensing images |
CN106840409A (en) * | 2017-01-23 | 2017-06-13 | 北京师范大学 | A kind of forest fire point detection method based on MODIS |
CN109829872B (en) * | 2019-01-22 | 2022-11-22 | 河南农业大学 | Multi-temporal multi-source remote sensing image fusion method for remote sensing of inland water body |
CN111008565B (en) * | 2019-11-04 | 2023-05-02 | 同济大学 | Fire trace detection method based on short wave infrared and thermal infrared data feature fusion |
CN112116242B (en) * | 2020-09-17 | 2022-08-16 | 福州福大经纬信息科技有限公司 | Bare soil change identification method combining multiple remote sensing indexes |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5878356A (en) * | 1995-06-14 | 1999-03-02 | Agrometrics, Inc. | Aircraft based infrared mapping system for earth based resources |
US5832187A (en) * | 1995-11-03 | 1998-11-03 | Lemelson Medical, Education & Research Foundation, L.P. | Fire detection systems and methods |
IL117521A0 (en) * | 1996-03-17 | 1996-10-31 | Israel Aircraft Ind Ltd Malat | A fire imaging system and method |
EP1048928A1 (en) * | 1999-04-27 | 2000-11-02 | Deutsches Zentrum für Luft- und Raumfahrt e.V. | Spaceborne hot temperature event (HTE) detection arrangement and multi-satellite HTE detection system |
-
2003
- 2003-07-09 IT IT000336A patent/ITRM20030336A1/en unknown
-
2004
- 2004-07-07 AT AT04745201T patent/ATE349677T1/en not_active IP Right Cessation
- 2004-07-07 EP EP04745201A patent/EP1642087B1/en active Active
- 2004-07-07 WO PCT/IT2004/000376 patent/WO2005005926A1/en active IP Right Grant
- 2004-07-07 DE DE602004003944T patent/DE602004003944D1/en active Active
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107368617A (en) * | 2016-05-13 | 2017-11-21 | 南京理工大学 | Air-ground detection infrared imaging system operating distance computational methods based on Lowtran7 air softwares |
CN107368617B (en) * | 2016-05-13 | 2020-08-11 | 南京理工大学 | Ground-air detection infrared imaging system action distance calculation method based on Lowtran7 atmospheric software |
Also Published As
Publication number | Publication date |
---|---|
WO2005005926A1 (en) | 2005-01-20 |
EP1642087A1 (en) | 2006-04-05 |
ITRM20030336A1 (en) | 2005-01-10 |
ITRM20030336A0 (en) | 2003-07-09 |
DE602004003944D1 (en) | 2007-02-08 |
ATE349677T1 (en) | 2007-01-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Giglio et al. | Active fire detection and characterization with the advanced spaceborne thermal emission and reflection radiometer (ASTER) | |
Wooster et al. | Sentinel-3 SLSTR active fire detection and FRP product: Pre-launch algorithm development and performance evaluation using MODIS and ASTER datasets | |
Xu et al. | New GOES imager algorithms for cloud and active fire detection and fire radiative power assessment across North, South and Central America | |
Garcia-Haro et al. | Monitoring fire-affected areas using Thematic Mapper data | |
Li et al. | A technique for detecting burn scars using MODIS data | |
Lyapustin et al. | Discrimination of biomass burning smoke and clouds in MAIAC algorithm | |
US8395119B2 (en) | Airborne/spaceborne oil spill determining system | |
Hutchinson et al. | Threshold functions for automated cloud analyses of global meteorological satellite imagery | |
Matheson et al. | Evaluating the effects of spatial resolution on hyperspectral fire detection and temperature retrieval | |
Spencer et al. | Exploring aerosols near clouds with high‐spatial‐resolution aircraft remote sensing during SEAC4RS | |
US8515716B2 (en) | Remote material identification process performance prediction tool | |
Zarei et al. | Validation of non-linear split window algorithm for land surface temperature estimation using Sentinel-3 satellite imagery: Case study; Tehran Province, Iran | |
EP1642087B1 (en) | Method and apparatus for automatically detecting and mapping, particularly for burnt areas without vegetation | |
Li et al. | Multi-information collaborative cloud identification algorithm in Gaofen-5 Directional Polarimetric Camera imagery | |
Lachérade et al. | ICARE: A physically-based model to correct atmospheric and geometric effects from high spatial and spectral remote sensing images over 3D urban areas | |
Egorov et al. | Generation and comprehensive validation of 30 m conterminous United States Landsat percent tree cover and forest cover loss annual products | |
Libonati et al. | Validation of the burned area ‘(V, W)’Modis algorithm in Brazil | |
Bostater Jr | Imaging derivative spectroscopy for vegetation dysfunction assessments | |
Xu et al. | Detecting forest fire omission error based on data fusion at subpixel scale | |
Team et al. | Discriminating clear-sky from cloud with modis algorithm theoretical basis document (mod35) | |
Buongiorno et al. | Recovery of spectral emissivity from thermal infrared multispectral scanner imagery acquired over a mountainous terrain: A case study from Mount Etna Sicily | |
Laneve et al. | SIGRI project: Products validation results | |
Ismaya et al. | New automated cloud and cloud-shadow detection using Landsat imagery | |
Gosteva et al. | Search of Changes in the Temperature of Urban Environment with Use of Satellite Data on the Example of the Krasnoyarsk | |
Paris | Remote sensing applications for freshwater systems |
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 |
|
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 |
|
17P | Request for examination filed |
Effective date: 20060109 |
|
GRAP | Despatch of communication of intention to grant a patent |
Free format text: ORIGINAL CODE: EPIDOSNIGR1 |
|
DAX | Request for extension of the european patent (deleted) | ||
GRAS | Grant fee paid |
Free format text: ORIGINAL CODE: EPIDOSNIGR3 |
|
GRAA | (expected) grant |
Free format text: ORIGINAL CODE: 0009210 |
|
AK | Designated contracting states |
Kind code of ref document: B1 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 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: BE Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20061227 Ref country code: DK Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20061227 Ref country code: FI Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20061227 Ref country code: CH Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20061227 Ref country code: CZ Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20061227 Ref country code: SI Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20061227 Ref country code: NL Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20061227 Ref country code: RO Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20061227 Ref country code: SK Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20061227 Ref country code: PL Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20061227 Ref country code: AT Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20061227 Ref country code: LI Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20061227 |
|
REG | Reference to a national code |
Ref country code: GB Ref legal event code: FG4D |
|
REG | Reference to a national code |
Ref country code: IE Ref legal event code: FG4D |
|
REF | Corresponds to: |
Ref document number: 602004003944 Country of ref document: DE Date of ref document: 20070208 Kind code of ref document: P |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: SE Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20070327 Ref country code: BG Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20070327 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: DE Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20070328 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: ES Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20070407 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: PT Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20070528 |
|
NLV1 | Nl: lapsed or annulled due to failure to fulfill the requirements of art. 29p and 29m of the patents act | ||
REG | Reference to a national code |
Ref country code: CH Ref legal event code: PL |
|
ET | Fr: translation filed | ||
PLBE | No opposition filed within time limit |
Free format text: ORIGINAL CODE: 0009261 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT |
|
26N | No opposition filed |
Effective date: 20070928 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: MC Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20070731 Ref country code: GR Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20070328 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: IE Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20070709 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: EE Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20061227 |
|
GBPC | Gb: european patent ceased through non-payment of renewal fee |
Effective date: 20080707 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: GB Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20080707 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: LU Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20070707 Ref country code: CY Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20061227 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: TR Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20061227 Ref country code: HU Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20070628 |
|
REG | Reference to a national code |
Ref country code: FR Ref legal event code: PLFP Year of fee payment: 13 |
|
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: FR Payment date: 20160613 Year of fee payment: 13 |
|
REG | Reference to a national code |
Ref country code: FR Ref legal event code: ST Effective date: 20180330 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: FR Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20170731 |
|
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: IT Payment date: 20230704 Year of fee payment: 20 |