CN115471743A - Method, device and equipment for determining open fire area in satellite image - Google Patents
Method, device and equipment for determining open fire area in satellite image Download PDFInfo
- Publication number
- CN115471743A CN115471743A CN202210986567.5A CN202210986567A CN115471743A CN 115471743 A CN115471743 A CN 115471743A CN 202210986567 A CN202210986567 A CN 202210986567A CN 115471743 A CN115471743 A CN 115471743A
- Authority
- CN
- China
- Prior art keywords
- infrared channel
- satellite image
- determining
- fire
- pixel
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 39
- 230000005855 radiation Effects 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 9
- 238000004590 computer program Methods 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 abstract description 9
- 230000001629 suppression Effects 0.000 abstract description 3
- 238000001514 detection method Methods 0.000 description 10
- 238000004364 calculation method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000000354 decomposition reaction Methods 0.000 description 3
- 230000006399 behavior Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000005215 recombination Methods 0.000 description 2
- 230000006798 recombination Effects 0.000 description 2
- 238000003915 air pollution Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/13—Satellite images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/188—Vegetation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Multimedia (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Astronomy & Astrophysics (AREA)
- Remote Sensing (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
The invention provides a method, a device and equipment for determining an open fire area in a satellite image. The method for determining the open fire area in the satellite image comprises the following steps: acquiring a satellite image shot by a preset satellite at a morning and evening moment, wherein the spatial resolution of a far infrared channel of the preset satellite is higher than that of a middle infrared channel; determining a target area where the fire point image element in the satellite image is located according to the intermediate infrared channel; and determining an open fire area in the target area according to the far infrared channel. The scheme of the invention can improve the position precision of the open fire area with strong fire behavior in the satellite remote sensing large-scale forest grassland fire at the morning and evening time, can provide more accurate positioning for fire monitoring, and provides powerful technical support for fire suppression and relief.
Description
Technical Field
The invention relates to the technical field of meteorological satellite observation, in particular to a method, a device and equipment for determining an open fire area in a satellite image.
Background
The forest and grassland fire is an important disaster for destroying forest and grassland resources, has the characteristics of strong burst, large influence range, serious loss and the like, and can cause serious air pollution.
The meteorological satellite has wide observation range, dense frequency and rich detection information, and is an important means for monitoring forest and grassland fires.
The existing meteorological satellite fire detection technology mainly uses an intermediate infrared channel which is sensitive to open fire information, but has the problems that the positioning precision is limited by the spatial resolution and the open fire information is easy to be exaggerated; when the fire point is large, the fire point is monitored only by using the intermediate infrared channel, and a plurality of pixels reach the fire point identification threshold. From the display effect of the fire condition multi-channel synthetic diagram and the given fire point positioning longitude and latitude information, a satellite remote sensing monitoring result which is obviously larger than an actual fire occurrence area appears, the fire field range is obviously overestimated, and the fire point monitoring precision is influenced.
Disclosure of Invention
The invention aims to solve the technical problem of providing a method, a device and equipment for determining an open fire area in a satellite image, which can improve the position precision of the open fire area with strong fire behavior in a satellite remote sensing large-range forest grassland fire at the morning and evening time, can provide more accurate positioning for fire monitoring, and provide powerful technical support for fire suppression and relief.
In order to solve the technical problems, the technical scheme of the invention is as follows:
a method for determining an open fire region in a satellite image comprises the following steps:
acquiring a satellite image shot by a preset satellite at morning and evening, wherein the spatial resolution of a far infrared channel of the preset satellite is higher than the spatial resolution of a middle infrared channel;
determining a target area where the fire point image element in the satellite image is located according to the intermediate infrared channel;
and determining an open fire area in the target area according to the far infrared channel.
Optionally, determining a target area where the fire pixel in the satellite image is located according to the mid-infrared channel includes:
determining the mixed pixel brightness temperature of the intermediate infrared channel, at least one pixel meeting a first preset condition as a fire point pixel in the satellite image; the first preset condition is determined according to the average value of the background brightness temperature of the middle infrared channel, the standard deviation of the background brightness temperature of the middle infrared channel, the brightness temperature difference of the middle infrared channel and the far infrared channel, the standard deviation of the background brightness temperature difference of the middle infrared channel and the far infrared channel and a first preset coefficient;
and determining a region formed by at least one fire point pixel in the satellite image as the target region.
Optionally, determining the mixed pixel brightness temperature of the mid-infrared channel and at least one pixel meeting a first preset condition as a fire pixel in the satellite image, including:
the first preset condition will be satisfied: t is 1 ≥(T 1BG +a×δT 1BG ) And T is 2 ≥(T 2BG +a×δT 2BG ) Is determined as a fire pixel in the satellite image;
wherein, T 1 Is the mixed pixel brightness temperature of the middle infrared channel;
T 1BG the average value of the background brightness temperature of the middle infrared channel is obtained;
δT 1BG standard deviation of background brightness temperature of the intermediate infrared channel;
T 2 the bright temperature difference of the middle infrared channel and the far infrared channel;
T 2BG the average value of the background bright temperature difference of the middle infrared channel and the far infrared channel is obtained;
δT 2BG the standard difference of the background bright temperature difference of the middle infrared channel and the far infrared channel is shown;
a is a first preset coefficient, and a is set as a first value for a region with solar radiation in the satellite image; for the area without sun irradiation in the satellite image, a is a second value, and the first value is larger than the second value.
Alternatively, if δ T 1BG When less than 2K, will delta T 1BG Set to 2K;
if delta T 1BG When greater than 3K, will delta T 1BG Set to 3K;
if delta T 2BG When less than 2K, will delta T 2BG Set to 2K;
if delta T 2BG When the temperature is more than 3K, the temperature will beδT 2BG Set to 3K;
the first value is 4 and the second value is 3.
Optionally, determining an open fire region in the target region according to the far infrared channel, including:
in the target area, the brightness temperature of the mixed pixel of the far infrared channel and at least one pixel meeting a second preset condition are determined as an open fire pixel in the satellite image; the second preset condition is determined according to the average value of the background brightness temperature of the far infrared channel, the standard deviation of the background brightness temperature of the far infrared channel and a second preset coefficient;
and determining the open fire area according to the area formed by at least one open fire point pixel.
Optionally, the second preset condition is to be satisfied: t is a unit of 3 ≥T 3BG +q×δT 3BG Is determined as the pixel of the open fire point in the satellite image;
wherein, T 3 The mixed pixel which is a far infrared channel is bright;
T 3BG the average value of the background brightness temperature of the far infrared channel;
δT 3BG standard deviation of background brightness temperature of a far infrared channel;
q is a second preset coefficient, and for the area with the sun irradiation in the satellite image, q is set as a third value; and for the area without sun irradiation in the satellite image, q is a fourth value, and the third value is larger than the fourth value.
Alternatively, when δ T 3BG When less than 1K, will δ T 3BG Is set to 1K;
when delta T 3BG When greater than 2K, will delta T 3BG Is set to be 2K;
the third value is 3 and the fourth value is 2.
The invention also provides a device for determining the open fire region in the satellite image, which comprises:
the system comprises an acquisition module, a processing module and a processing module, wherein the acquisition module is used for acquiring a satellite image at morning and evening time shot by a preset satellite, and the spatial resolution of a far infrared channel of the preset satellite is higher than the spatial resolution of a middle infrared channel;
the processing module is used for determining a target area where the fire point image element in the satellite image is located according to the intermediate infrared channel; and determining an open fire area in the target area according to the far infrared channel.
The present invention also provides a computing device comprising: a processor, a memory storing a computer program which, when executed by the processor, performs the method as described above.
The present invention also provides a computer readable storage medium storing instructions which, when executed on a computer, cause the computer to perform the method as described above.
The scheme of the invention at least comprises the following beneficial effects:
according to the scheme, the satellite image at the morning and evening moment is shot through a preset satellite, and the spatial resolution of a far infrared channel of the preset satellite is higher than that of a middle infrared channel; determining a target area where the fire point image element in the satellite image is located according to the intermediate infrared channel; and determining an open fire region in the target region according to the far infrared channel. The satellite remote sensing system can improve the position precision of the open fire area with strong fire in the large-range forest grassland fire by satellite remote sensing at the morning and evening time, can provide more accurate positioning for fire monitoring, and provides powerful technical support for fire suppression and relief.
Drawings
Fig. 1 is a schematic flow chart of the method for determining an open fire region in a satellite image according to the present invention;
FIG. 2 is a schematic diagram of an image element of the present invention initially identified as a fire point through a mid-infrared channel;
FIG. 3 is a schematic diagram of the present invention, in which the middle infrared channel is initially identified as the fire point image element, the accurate position of the fire point is further identified by using the far infrared channel;
fig. 4 is a block diagram of the determination device of the open fire region in the satellite image according to the present invention.
Detailed Description
Exemplary embodiments of the present disclosure will be described in more detail below with reference to the accompanying drawings. While exemplary embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure may be embodied in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the disclosure to those skilled in the art.
As shown in fig. 1, an embodiment of the present invention provides a method for determining an open fire region in a satellite image, including:
and 13, determining an open fire area in the target area according to the far infrared channel.
In this embodiment, the preset satellite may be an FY-3E satellite (wind cloud three E star), and the FY-3E satellite provides observation at the morning and evening time, so as to increase the monitoring accuracy, frequency and timeliness of forest and grassland fires at the morning and evening time;
the FY-3E satellite has two mid-infrared channels with one-kilometer resolution, the wavelengths are respectively 3.8 micrometers and 4.05 micrometers, fire point information can be monitored, and the fire point positioning accuracy can be one kilometer;
the FY-3E satellite has two far infrared channels with a resolution of 250 meters, with wavelengths of 10.8 μm and 12 μm, respectively, the resolution of the far infrared channels being 4 times higher than the resolution of the mid infrared channels.
In the embodiment of the invention, firstly, the target area where the fire point image element is located in the satellite image shot by the preset satellite at the morning and evening time is determined through the intermediate infrared channel of the preset satellite, and then the open fire area is further determined in the target area according to the far infrared channel of the preset satellite, so that the position accuracy of the open fire area with strong fire behavior in the satellite remote sensing of the large-range forest and grassland fire at the morning and evening time can be improved, more accurate positioning can be provided for fire monitoring, and powerful technical support can be provided for fire fighting and disaster relief.
In an optional embodiment of the present invention, step 12 may include:
step 121, determining the mixed pixel brightness temperature of the intermediate infrared channel, at least one pixel meeting a first preset condition as a fire point pixel in the satellite image; the first preset condition is determined according to the average value of the background brightness temperature of the intermediate infrared channel, the standard deviation of the background brightness temperature of the intermediate infrared channel, the brightness temperature difference of the intermediate infrared channel and the far infrared channel, the standard deviation of the background brightness temperature difference of the intermediate infrared channel and the far infrared channel and a first preset coefficient;
specifically, the step 121 may include:
step 1211, to satisfy the first preset condition: t is a unit of 1 ≥(T 1BG +a×δT 1BG ) And T is 2 ≥(T 2BG +a×δT 2BG ) Is determined as a fire pixel in the satellite image;
wherein, T 1 Is the mixed pixel brightness temperature of the middle infrared channel;
T 1BG the average value of the background brightness temperature of the intermediate infrared channel is obtained;
δT 1BG standard deviation of background brightness temperature of the intermediate infrared channel;
T 2 the bright temperature difference of the middle infrared channel and the far infrared channel;
T 2BG the average value of the background bright temperature difference of the middle infrared channel and the far infrared channel is obtained;
δT 2BG the standard deviation of the background bright temperature difference of the middle infrared channel and the far infrared channel is obtained;
a is a first preset coefficient, and a is set as a first value for a region with solar radiation in the satellite image; and a is a second value for the area without the sun irradiation in the satellite image, and the first value is larger than the second value.
Wherein if δ T 1BG When less than 2K, will delta T 1BG Is set to be 2K;
if delta T 1BG When greater than 3K, will δ T 1BG Set to 3K;
if delta T 2BG When less than 2K, will delta T 2BG Is set to be 2K;
if delta T 2BG When greater than 3K, will δ T 2BG Set to 3K;
the first value is 4 and the second value is 3.
And step 122, determining a region formed by at least one fire point pixel in the satellite image as the target region.
It should be noted that, since the FY-3E satellite is a morning and evening satellite, one part of the image is irradiated by the sun, the other part of the image is not irradiated by the sun, and the fire point identification of two different areas with or without solar radiation is different in calculation method, it is determined whether the pixel can receive solar radiation, mainly according to the solar zenith angle, and when the solar zenith angle is greater than a preset angle (90 degrees), the pixel is an area without solar radiation.
In this embodiment, the mixed pixel brightness temperature of the mid-infrared channel satisfies a first preset condition: t is 1 ≥(T 1BG +a×δT 1BG ) And T 2 ≥(T 2BG +a×δT 2BG ) When the satellite image is in a hot state, the pixel corresponding to the mixed pixel brightness temperature meeting the first preset condition is used as a fire pixel in the satellite image, and the area formed by at least one fire pixel in the satellite image is determined as the target area;
wherein, the T is 1 The mixed pixel brightness temperature of the middle infrared channel is set;
T 1BG the average value of the background brightness temperature of the intermediate infrared channel is obtained;
δT 1BG standard deviation of background brightness temperature of the intermediate infrared channel;
T 2 the bright temperature difference of the middle infrared channel and the far infrared channel;
T 2BG the background bright temperature difference of the middle infrared channel and the far infrared channel;
δT 2BG the standard difference of the background bright temperature difference of the middle infrared channel and the far infrared channel is shown;
a is a first preset coefficient, and for the satellite image, the coefficient isA region irradiated by the sun, a being set to a first value; for the area without the sun irradiation in the satellite image, a is a second value, and the first value is larger than the second value; if delta T 1BG When less than 2K, will delta T 1BG Set to 2K; if delta T 1BG When greater than 3K, will delta T 1BG Set to 3K; if delta T 2BG When less than 2K, will delta T 2BG Set to 2K; if delta T 2BG When greater than 3K, will delta T 2BG Set to 3K; the first value is preferably 4 and the second value is preferably 3.
Here, the average value of the background brightness temperature is an average value of the brightness temperatures of the pixels without open fire around the detection pixel;
calculating that the pixel selected by the background brightness temperature can be the same as the underlay surface type and vegetation coverage of the detection pixel as much as possible;
in the daytime, the difference between the bright temperature difference of the mid infrared and the far infrared caused by the specular reflection of solar radiation on the cloud surface and the water body can reach dozens of K;
although the reflection of solar radiation on a middle infrared channel on an exposed earth surface is relatively weak, due to the fact that the requirement on fire point judgment sensitivity is high, even if the difference is only a few degrees, pixels which are not fire points originally can reach a heating threshold value of fire point judgment, and therefore in the area with solar radiation, factors such as a cloud area, solar flares, suspected fire point pixels, a water body, a desert area, areas with different vegetation coverage degrees and the like are removed in background brightness temperature calculation.
Cloud detection products can be used for cloud area pixel judgment;
judging usable land utilization data of water bodies and desert regions; the solar flare area can use a flare angle value, and when the flare angle value is smaller than a first preset value, fire point judgment is not made, wherein the first preset value can be 15 degrees;
if the clear sky vegetation pixels in the neighborhood around the detection pixel meet the following conditions, the clear sky vegetation pixels are taken as suspected fire point pixels:
T MIR ＞(T MIR_AVG +ΔT MIR ) And T M-F ＞(T M-F_AVG + 8K), or T MIR ＞T MIR_WM ；
Wherein, theT is described MIR_AVG The average value of the brightness and the temperature of vegetation pixels in clear sky in r x r pixels around an infrared channel in a detection pixel is T M-F_AVG The mean value of difference between the bright temperature of the infrared channel and the bright temperature of the far infrared channel in the clear sky vegetation pixel in the r x r pixel around the detection pixel is delta T MIR The reference value was 10K, T MIR_WM The reference value is 330K, the initial value of r is set to 7, and the non-suspected high-temperature pixels in the neighborhood of the r × r pixels around the detection pixel are used as the pixels in the background area. The number of pixels in the background area is less than 20% of the number of pixels in the neighborhood, the neighborhood is expanded to 9 multiplied by 9 pixels, 11 multiplied by 11 pixels, \ 8230, and the number of pixels is up to 49 multiplied by 49 pixels (reference value); if T still does not reach 20% MIR If T is less than 345k, the determination process for the detected pixel can be abandoned, if T is MIR If the value is greater than or equal to 345k, the value of r can be increased continuously until 81 x 81 pixels (reference value) still reaches less than 20%, and the judgment processing of the detection pixel can be abandoned.
Because the FY-3E satellite is a morning and evening satellite, in the area without sun irradiation, the influence of factors such as cloud areas, solar flares, water bodies, desert areas, different vegetation coverage mixed areas and the like is not considered in the calculation of the background brightness temperature;
as shown in fig. 2, the thick solid line boxes a and B represent that the mid-infrared channel passes through a target area where the fire point image element is initially determined to be located, and the spatial resolution of the image element is 1km; the thin solid line boxes represent the picture elements of the far infrared channel, the resolution of which is 250m.
Therefore, according to the brightness temperature of the mixed pixels of the intermediate infrared channel, the area formed by at least one fire point pixel in the satellite image can be determined. This enables the approximate area in which the fire pixels are located to be locked, reducing the reference range for determining the location of a particular fire point.
In the above-described embodiments of the invention, in determining the fire pixels in the satellite images,
the open fire is obviously higher than the temperature and radiance of other ground features on the ground surface, the pixel radiance containing the open fire can be regarded as the linear combination of an open fire area and a non-open fire area on the basis of a linear mixed pixel decomposition technology:
N imix ＝P*N ihi +(1-P)*N ibg ；
wherein, the N is imix For fire point image element channel i radiance, the N ihi Is the radiance of the open fire region, N ibg The radiance of the non-open-fire area, and P is the proportion of the area of the sub-pixel open-fire area in the pixel to the area of the pixel.
The temperature difference Δ T between the pixels containing fire and non-fire may be calculated by the following formula:
wherein, delta T i The difference between the lighting pixel and the background brightness temperature is T imix Channel i light temperature, said T iB For background brightness, N imix For fire pixel channel i radiance, said N iB For background radiance, said V i Is the center wave number, C, of channel i 1 ＝1.1910659×10 -5 mW(m 2 ，sr，cm -4 )，C 2 ＝1.438833K/cm -1 。
Based on the central wave number of the middle infrared channel and the central wave number of the far infrared channel, supposing that the background temperature is 290K, the temperature of the open fire region is 750K, the proportion of the area of the sub-pixel fire point to the area of the pixel fire point is changed from 0.0001 to 0.005, and the parameters are substituted into
Under the condition of the same spatial resolution, the brightness temperature increment of the open fire area caused by the middle infrared channel and the far infrared channel is obviously different, and the brightness temperature increment caused by the middle infrared channel is obviously higher than that caused by the far infrared channel;
because the resolution ratio of a far infrared channel of the FY-3E satellite is 1km, the resolution ratio of a middle infrared channel is 250m, the resolution ratio of the far infrared channel is 4 times of that of the middle infrared channel, and the area ratio of the fire point area of the same sub-pixel in the far infrared channel pixel is 16 times of that of the middle infrared channel, the difference between the brightness temperature increment of the far infrared channel mixed pixel and the middle infrared channel is obviously reduced;
when the area proportion of the open fire region reaches 0.0004, the brightness temperature increment of the middle infrared channel is about 19K, the brightness temperature increment of the far infrared channel is about 7K, and when the area proportion of the open fire region is 0.005, the brightness temperature increment of the middle infrared channel is about 78K, the brightness temperature increment of the far infrared channel is about 71K, and the increment of the far infrared channel is very close to that of the middle infrared channel, so that a basis is provided for judging and identifying fire point pixels of the far infrared channel with higher spatial resolution than that of the middle infrared channel, and the positioning accuracy is improved.
In another optional embodiment of the present invention, step 13 may include:
step 131, in the target area, determining the mixed pixel brightness temperature of the far infrared channel and at least one pixel meeting a second preset condition as an open fire pixel in the satellite image; the second preset condition is determined according to the average value of the background brightness temperature of the far infrared channel, the standard deviation of the background brightness temperature of the far infrared channel and a second preset coefficient;
specifically, the second preset condition is to be satisfied: t is 3 ≥T 3BG +q×δT 3BG The at least one pixel is determined as an open fire pixel in the satellite image;
wherein, T 3 The mixed pixel which is a far infrared channel is bright;
T 3BG the average value of the background brightness temperature of the far infrared channel;
δT 3BG standard deviation of background brightness temperature of a far infrared channel;
q is a second preset coefficient, and for the area with the sun irradiation in the satellite image, q is set as a third value; and q is a fourth value for the area without the sun irradiation in the satellite image, and the third value is larger than the fourth value.
Wherein when delta T 3BG When less than 1K, will δ T 3BG Is set to 1K; when delta T 3BG When greater than 2K, will δ T 3BG Set to 2K;
the third value is 3 and the fourth value is 2.
And 132, determining the area formed by at least one open fire point pixel.
It should be noted that the far infrared channel has a resolution obviously higher than 1km for wild fire, and fire identification information of the mid-infrared channel needs to be combined when the far infrared channel is used for identifying the fire.
In the background brightness temperature estimation of the fire point of the far infrared channel, the far infrared channel pixel in the fire point coverage range of the medium infrared channel can be used, the medium infrared channel with the resolution of 1km identifies and marks the pixel of the fire point, and the positioning precision of the fire point is kilometer level. Fire point information of the mid-infrared channel identification provides necessary conditions for identifying far-infrared channel fire point pixels in the coverage range of the mid-infrared channel, namely one or more pixels are bound to contain open fire in 16 far-infrared fire point pixels in the coverage range of the mid-infrared fire point pixels;
and (3) taking the coverage range of the fire point pixel of the intermediate infrared channel identification as the area of the far infrared channel fire point identification with the resolution of 250m and the statistical calculation of the background brightness temperature, and judging whether the far infrared channel pixel in the area is the fire point pixel by pixel.
In this embodiment, in the target area, the mixed pixel brightness temperature of the far infrared channel satisfies a second preset condition: t is 3 ≥T 3BG +q×δT 3BG And then, taking the pixel corresponding to the mixed pixel brightness temperature meeting the second preset condition as an open fire pixel in the satellite image, and determining the open fire region by using the region formed by at least one open fire pixel.
Wherein, T 3 The mixed pixel which is a far infrared channel is bright; t is a unit of 3BG The average value of the background brightness temperature of the far infrared channel is obtained;
the area of the background brightness temperature statistics is the fire point pixel coverage area judged and identified by the middle infrared channel, and if no background pixel meeting the condition is found, the far infrared neighborhood range can be expanded to calculate the pixel brightness temperature average value as the background brightness temperature;
δT 3BG standard deviation of background brightness temperature of a far infrared channel; q is a second preset coefficient, and for the area with the sun irradiation in the satellite image, q is set as a third value; for the area without sun irradiation in the satellite image, q is a fourth value, and the third value is greater than the fourth value(ii) a When delta T 3BG When less than 1K, will delta T 3BG Is set to 1K; when delta T 3BG When greater than 2K, will delta T 3BG Set to 2K;
the third value is preferably 3 and the fourth value is preferably 2.
As shown in fig. 3, the thick solid line boxes a and B are two pixels of which the intermediate infrared channel is primarily determined as a fire point, 5 areas a, B, c, d and e are respectively arranged in the target areas a and B of which the intermediate infrared channel is primarily determined as a fire point pixel, the pixels are pixels of which the far infrared channel is determined as an open fire point area, the resolution of the pixels is 250 meters, and the position accuracy of the pixels in the 5 areas is further improved compared with the pixels a and B. The fire point pixel position of the far infrared channel of 250 meters is determined, and fire point position information with higher positioning precision than that of the middle infrared channel with 1km spatial resolution is finally realized.
In the embodiment of the invention, by the method for determining the open fire region in the satellite image, the position accuracy of the strong fire region in a satellite remote sensing large-range fire scene at the morning and evening can be improved, more accurate positioning can be provided for fire monitoring, and powerful technical support can be provided for fire fighting and relief.
As shown in fig. 4, an embodiment of the present invention further provides an apparatus 40 for determining an open fire region in a satellite image, where the apparatus 40 includes:
an obtaining module 41, configured to obtain a satellite image of a preset satellite at a morning and evening time, where a spatial resolution of a far infrared channel of the preset satellite is higher than a spatial resolution of a mid infrared channel;
the processing module 42 is configured to determine a target area where the fire pixel in the satellite image is located according to the mid-infrared channel; and determining an open fire area in the target area according to the far infrared channel.
Optionally, determining a target area where the fire pixel in the satellite image is located according to the mid-infrared channel includes:
determining the mixed pixel brightness temperature of the intermediate infrared channel and at least one pixel meeting a first preset condition as a fire point pixel in the satellite image; the first preset condition is determined according to the average value of the background brightness temperature of the middle infrared channel, the standard deviation of the background brightness temperature of the middle infrared channel, the brightness temperature difference of the middle infrared channel and the far infrared channel, the standard deviation of the background brightness temperature difference of the middle infrared channel and the far infrared channel and a first preset coefficient;
and determining a region formed by at least one fire point pixel in the satellite image as the target region.
Optionally, determining the mixed pixel brightness temperature of the mid-infrared channel and at least one pixel meeting a first preset condition as a fire pixel in the satellite image, including:
the first preset condition will be satisfied: t is 1 ≥(T 1BG +a×δT 1BG ) And T 2 ≥(T 2BG +a×δT 2BG ) Is determined as a fire pixel in the satellite image;
wherein, T 1 The mixed pixel brightness temperature of the middle infrared channel is set;
T 1BG the average value of the background brightness temperature of the middle infrared channel is obtained;
δT 1BG standard deviation of background brightness temperature of the intermediate infrared channel;
T 2 the bright temperature difference of the middle infrared channel and the far infrared channel;
T 2BG the average value of the background bright temperature difference of the middle infrared channel and the far infrared channel is obtained;
δT 2BG the standard deviation of the background bright temperature difference of the middle infrared channel and the far infrared channel is obtained;
a is a first preset coefficient, and a is set as a first value for a region with solar radiation in the satellite image; for the area without sun irradiation in the satellite image, a is a second value, and the first value is larger than the second value.
Alternatively, if δ T 1BG When less than 2K, will delta T 1BG Set to 2K;
if delta T 1BG When greater than 3K, will delta T 1BG Set to 3K;
if delta T 2BG When less than 2K, will delta T 2BG Set to 2K;
if delta T 2BG When greater than 3K, will delta T 2BG Set to 3K;
the first value is 4 and the second value is 3.
Optionally, determining an open fire region in the target region according to the far infrared channel, including:
in the target area, the brightness temperature of the mixed pixel of the far infrared channel and at least one pixel meeting a second preset condition are determined as an open fire pixel in the satellite image; the second preset condition is determined according to the average value of the background brightness temperature of the far infrared channel, the standard deviation of the background brightness temperature of the far infrared channel and a second preset coefficient;
and determining the open fire area according to the area formed by at least one open fire point pixel.
Optionally, the second preset condition is to be satisfied: t is a unit of 3 ≥T 3BG +q×δT 3BG Is determined as the pixel of the open fire point in the satellite image;
wherein, T 3 The mixed pixel which is a far infrared channel is bright;
T 3BG the average value of the background brightness temperature of the far infrared channel;
δT 3BG standard deviation of background brightness temperature of a far infrared channel;
q is a second preset coefficient, and for the area with the sun irradiation in the satellite image, q is set as a third value; and q is a fourth value for the area without the sun irradiation in the satellite image, and the third value is larger than the fourth value.
Alternatively, when δ T 3BG When less than 1K, will delta T 3BG Is set as 1K; when delta T 3BG When greater than 2K, will delta T 3BG Is set to be 2K;
the third value is 3 and the fourth value is 2.
It should be noted that the apparatus is an apparatus corresponding to the above method, and all the implementations in the above method embodiment are applicable to the embodiment of the apparatus, and the same technical effects can be achieved.
Embodiments of the present invention also provide a computing device, comprising: a processor, a memory storing a computer program which, when executed by the processor, performs the method as claimed above. All the implementation manners in the above method embodiment are applicable to the embodiment, and the same technical effect can be achieved
Embodiments of the present invention also provide a computer-readable storage medium, including: storing instructions that, when executed on a computer, cause the computer to perform the method as described above. All the implementation manners in the above method embodiment are applicable to this embodiment, and the same technical effect can be achieved.
Those of ordinary skill in the art will appreciate that the various illustrative elements and algorithm steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware or combinations of computer software and electronic hardware. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the implementation. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present invention.
It can be clearly understood by those skilled in the art that, for convenience and simplicity of description, the specific working processes of the above-described systems, apparatuses and units may refer to the corresponding processes in the foregoing method embodiments, and are not described herein again.
In the embodiments provided in the present invention, it should be understood that the disclosed apparatus and method may be implemented in other manners. For example, the above-described apparatus embodiments are merely illustrative, and for example, the division of the units is only one logical division, and other divisions may be realized in practice, for example, a plurality of units or components may be combined or integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection through some interfaces, devices or units, and may be in an electrical, mechanical or other form.
The units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of the embodiment.
In addition, functional units in the embodiments of the present invention may be integrated into one processing unit, or each unit may exist alone physically, or two or more units are integrated into one unit.
The functions may be stored in a computer-readable storage medium if they are implemented in the form of software functional units and sold or used as separate products. Based on such understanding, the technical solution of the present invention or a part thereof which substantially contributes to the prior art may be embodied in the form of a software product, which is stored in a storage medium and includes several instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U disk, a removable hard disk, a ROM, a RAM, a magnetic disk or an optical disk, and various media capable of storing program codes.
Furthermore, it is to be noted that in the device and method of the invention, it is obvious that the individual components or steps can be decomposed and/or recombined. These decompositions and/or recombinations are to be regarded as equivalents of the present invention. Also, the steps of performing the series of processes described above may naturally be performed chronologically in the order described, but need not necessarily be performed chronologically, and some steps may be performed in parallel or independently of each other. It will be understood by those skilled in the art that all or any of the steps or elements of the method and apparatus of the present invention may be implemented in any computing device (including processors, storage media, etc.) or network of computing devices, in hardware, firmware, software, or any combination thereof, which can be implemented by those skilled in the art using their basic programming skills after reading the description of the present invention.
The object of the invention is thus also achieved by a program or a set of programs running on any computing device. The computing device may be a general purpose device as is well known. The object of the invention is thus also achieved solely by providing a program product containing program code for implementing the method or device. That is, such a program product also constitutes the present invention, and a storage medium storing such a program product also constitutes the present invention. It is to be understood that the storage medium may be any known storage medium or any storage medium developed in the future. It is further noted that in the apparatus and method of the present invention, it is apparent that each component or step can be decomposed and/or recombined. These decompositions and/or recombinations are to be regarded as equivalents of the present invention. Also, the steps of executing the series of processes described above may naturally be executed chronologically in the order described, but need not necessarily be executed chronologically. Some steps may be performed in parallel or independently of each other.
While the foregoing is directed to the preferred embodiment of the present invention, it will be understood by those skilled in the art that various changes and modifications may be made without departing from the spirit and scope of the invention as defined in the appended claims.
Claims (10)
1. A method for determining an open fire region in a satellite image, comprising:
acquiring a satellite image shot by a preset satellite at morning and evening, wherein the spatial resolution of a far infrared channel of the preset satellite is higher than the spatial resolution of a middle infrared channel;
determining a target area where the fire point image element in the satellite image is located according to the intermediate infrared channel;
and determining an open fire area in the target area according to the far infrared channel.
2. The method for determining the open fire region in the satellite image according to claim 1, wherein determining the target region in which the fire pixel is located in the satellite image according to the mid-infrared channel comprises:
determining the mixed pixel brightness temperature of the intermediate infrared channel and at least one pixel meeting a first preset condition as a fire point pixel in the satellite image; the first preset condition is determined according to the average value of the background brightness temperature of the intermediate infrared channel, the standard deviation of the background brightness temperature of the intermediate infrared channel, the brightness temperature difference of the intermediate infrared channel and the far infrared channel, the standard deviation of the background brightness temperature difference of the intermediate infrared channel and the far infrared channel and a first preset coefficient;
and determining a region formed by at least one fire point pixel in the satellite image as the target region.
3. The method for determining the open fire region in the satellite image according to claim 2, wherein the step of determining the mixed pixel brightness temperature of the intermediate infrared channel and at least one pixel meeting a first preset condition as the fire pixel in the satellite image comprises the following steps:
the first preset condition will be satisfied: t is 1 ≥(T 1BG +a×δT 1BG ) And T 2 ≥(T 2BG +a×δT 2BG ) Is determined as a fire pixel in the satellite image;
wherein, T 1 Is the mixed pixel brightness temperature of the middle infrared channel;
T 1BG the average value of the background brightness temperature of the middle infrared channel is obtained;
δT 1BG standard deviation of background brightness temperature of the intermediate infrared channel;
T 2 the bright temperature difference of the middle infrared channel and the far infrared channel;
T 2BG the average value of the background bright temperature difference of the middle infrared channel and the far infrared channel is obtained;
δT 2BG standard for background bright temperature difference of middle infrared channel and far infrared channelDifference;
a is a first preset coefficient, and a is set as a first value for a region with solar radiation in the satellite image; and a is a second value for the area without the sun irradiation in the satellite image, and the first value is larger than the second value.
4. The method for determining an open flame region in a satellite image according to claim 3,
if delta T 1BG When less than 2K, will delta T 1BG Is set to be 2K;
if delta T 1BG When greater than 3K, will delta T 1BG Set to 3K;
if delta T 2BG When less than 2K, will delta T 2BG Set to 2K;
if delta T 2BG When greater than 3K, will δ T 2BG Set to 3K;
the first value is 4 and the second value is 3.
5. The method for determining an open fire region in a satellite image according to claim 1, wherein determining an open fire region in the target region based on the far-infrared channel comprises:
in the target area, the brightness temperature of the mixed pixel of the far infrared channel and at least one pixel meeting a second preset condition are determined as an open fire pixel in the satellite image; the second preset condition is determined according to the average value of the background brightness temperature of the far infrared channel, the standard deviation of the background brightness temperature of the far infrared channel and a second preset coefficient;
and determining the open fire area by using the area formed by at least one open fire pixel.
6. The method for determining an open flame region in a satellite image according to claim 5,
the second preset condition will be satisfied: t is 3 ≥T 3BG +q×δT 3BG Is determined as the pixel of the open fire point in the satellite image;
wherein, T 3 The mixed pixel which is a far infrared channel is bright;
T 3BG is the average value of background brightness temperature of far infrared channel
δT 3BG Standard deviation of background brightness temperature of a far infrared channel;
q is a second preset coefficient, and for the area with the sun irradiation in the satellite image, q is set as a third value; and q is a fourth value for the area without the sun irradiation in the satellite image, and the third value is larger than the fourth value.
7. The method for determining an open fire region in a satellite image according to claim 6,
when delta T 3BG When less than 1K, will delta T 3BG Is set to 1K;
when delta T 3BG When greater than 2K, will delta T 3BG Is set to be 2K;
the third value is 3 and the fourth value is 2.
8. An apparatus for determining an area of an open flame in a satellite image, comprising:
the system comprises an acquisition module, a processing module and a processing module, wherein the acquisition module is used for acquiring a satellite image at morning and evening time shot by a preset satellite, and the spatial resolution of a far infrared channel of the preset satellite is higher than the spatial resolution of a middle infrared channel;
the processing module is used for determining a target area where the fire point image element in the satellite image is located according to the intermediate infrared channel; and determining an open fire region in the target region according to the far infrared channel.
9. A computing device, comprising: a processor, a memory storing a computer program which, when executed by the processor, performs the method of any of claims 1 to 7.
10. A computer-readable storage medium storing instructions that, when executed on a computer, cause the computer to perform the method of any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210986567.5A CN115471743A (en) | 2022-08-17 | 2022-08-17 | Method, device and equipment for determining open fire area in satellite image |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210986567.5A CN115471743A (en) | 2022-08-17 | 2022-08-17 | Method, device and equipment for determining open fire area in satellite image |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115471743A true CN115471743A (en) | 2022-12-13 |
Family
ID=84368170
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210986567.5A Pending CN115471743A (en) | 2022-08-17 | 2022-08-17 | Method, device and equipment for determining open fire area in satellite image |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115471743A (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080027649A1 (en) * | 2004-11-30 | 2008-01-31 | Deutsches Zentrum Fur Luft-Und Raumfahrt E.V. | Method for Recognising Hot Targets on the Earth |
| US20090262978A1 (en) * | 2006-07-28 | 2009-10-22 | Telespazio S.P.A. | Automatic Detection Of Fires On Earth's Surface And Of Atmospheric Phenomena Such As Clouds, Veils, Fog Or The Like, Using A Satellite System |
| CN113340432A (en) * | 2021-06-09 | 2021-09-03 | 广东电网有限责任公司 | Fire monitoring method and system based on stationary meteorological satellite |
| CN114550006A (en) * | 2022-02-23 | 2022-05-27 | 黑龙江省生态气象中心(东北卫星气象数据中心) | Polar-orbit meteorological satellite fire point identification system, storage medium and equipment |
| CN114663776A (en) * | 2022-02-25 | 2022-06-24 | 星图智绘(西安)数字科技有限公司 | GK-2A-based high-frequency space self-adaptive threshold fire remote sensing monitoring method |
-
2022
- 2022-08-17 CN CN202210986567.5A patent/CN115471743A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080027649A1 (en) * | 2004-11-30 | 2008-01-31 | Deutsches Zentrum Fur Luft-Und Raumfahrt E.V. | Method for Recognising Hot Targets on the Earth |
| US20090262978A1 (en) * | 2006-07-28 | 2009-10-22 | Telespazio S.P.A. | Automatic Detection Of Fires On Earth's Surface And Of Atmospheric Phenomena Such As Clouds, Veils, Fog Or The Like, Using A Satellite System |
| CN113340432A (en) * | 2021-06-09 | 2021-09-03 | 广东电网有限责任公司 | Fire monitoring method and system based on stationary meteorological satellite |
| CN114550006A (en) * | 2022-02-23 | 2022-05-27 | 黑龙江省生态气象中心(东北卫星气象数据中心) | Polar-orbit meteorological satellite fire point identification system, storage medium and equipment |
| CN114663776A (en) * | 2022-02-25 | 2022-06-24 | 星图智绘(西安)数字科技有限公司 | GK-2A-based high-frequency space self-adaptive threshold fire remote sensing monitoring method |
Non-Patent Citations (3)
| Title |
|---|
| 胡秀清等: "FY-1C/1D全球海上气溶胶业务反演算法研究", 《海洋学报(中文版)》, pages 56 - 65 * |
| 郑伟等: "基于FY-3D/MERSI-Ⅱ远红外数据的火情监测研究", 《红外与毫米波学报》, pages 120 - 127 * |
| 陈洁等: "Himawari-8静止气象卫星草原火监测分析", 《自然灾害学报》, pages 197 - 204 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Hua et al. | The progress of operational forest fire monitoring with infrared remote sensing | |
| CN112435207B (en) | Forest fire monitoring and early warning method based on sky-ground integration | |
| Lee et al. | Improved detection of hotspots using the AVHRR 3.7-um channel | |
| CN101527073B (en) | Fire detection system and method thereof | |
| CN108648400A (en) | One kind is based on multispectral transmission line forest fire exploration prior-warning device and method for early warning | |
| CN108682105B (en) | One kind is based on multispectral transmission line forest fire exploration prior-warning device and method for early warning | |
| CN111982291A (en) | Fire point positioning method, device and system based on unmanned aerial vehicle | |
| Grammalidis et al. | A multi-sensor network for the protection of cultural heritage | |
| Ganci et al. | HOTSAT: a multiplatform system for the thermal monitoring of volcanic activity using satellite data | |
| CN113340432B (en) | Fire monitoring method and system based on stationary meteorological satellite | |
| Werkmeister et al. | Comparing satellite-to ground-based automated and manual cloud coverage observations–a case study | |
| CN112862150A (en) | Forest fire early warning method based on image and video multi-model | |
| CN115507959A (en) | Infrared radiation characteristic analysis method for target detection | |
| CN109945969B (en) | Method and device for determining earth radiation balance based on meteorological satellite observation | |
| Poduri et al. | Visibility monitoring using mobile phones | |
| CN114112065A (en) | Method for judging and recognizing fire danger by satellite remote sensing | |
| CN115471743A (en) | Method, device and equipment for determining open fire area in satellite image | |
| CN116597592A (en) | Method and system for monitoring fire points of polar orbit meteorological satellite in morning and evening time | |
| CN201435130Y (en) | Fire detection system | |
| Simpson et al. | Improved cloud top height retrieval under arbitrary viewing and illumination conditions using AVHRR data | |
| Bosch et al. | A ground system for early forest fire detection based on infrared signal processing | |
| CN114550006A (en) | Polar-orbit meteorological satellite fire point identification system, storage medium and equipment | |
| CN110133668B (en) | Infrared detection method and device for offshore targets | |
| Ramaswamy et al. | An adaptive tracker for ShipIR/NTCS | |
| Valero et al. | Flame filtering and perimeter localization of wildfires using aerial thermal imagery |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |