CN108304780A - A kind of crop straw burning fire point remote-sensing monitoring method based on No. three satellites of wind and cloud - Google Patents

A kind of crop straw burning fire point remote-sensing monitoring method based on No. three satellites of wind and cloud Download PDF

Info

Publication number
CN108304780A
CN108304780A CN201711471773.8A CN201711471773A CN108304780A CN 108304780 A CN108304780 A CN 108304780A CN 201711471773 A CN201711471773 A CN 201711471773A CN 108304780 A CN108304780 A CN 108304780A
Authority
CN
China
Prior art keywords
channel
image
fire point
cloud
satellites
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
Application number
CN201711471773.8A
Other languages
Chinese (zh)
Inventor
赵辉杰
吴锋
王玉玺
吴永欣
张帅
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
CETC 27 Research Institute
Original Assignee
CETC 27 Research Institute
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by CETC 27 Research Institute filed Critical CETC 27 Research Institute
Priority to CN201711471773.8A priority Critical patent/CN108304780A/en
Publication of CN108304780A publication Critical patent/CN108304780A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/10Terrestrial scenes
    • G06V20/13Satellite images
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/20Image preprocessing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/20Image preprocessing
    • G06V10/26Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion
    • G06V10/267Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion by performing operations on regions, e.g. growing, shrinking or watersheds
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/20Image preprocessing
    • G06V10/24Aligning, centring, orientation detection or correction of the image
    • G06V10/247Aligning, centring, orientation detection or correction of the image by affine transforms, e.g. correction due to perspective effects; Quadrilaterals, e.g. trapezoids

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Multimedia (AREA)
  • Theoretical Computer Science (AREA)
  • Astronomy & Astrophysics (AREA)
  • Remote Sensing (AREA)
  • Image Processing (AREA)

Abstract

The invention discloses a kind of crop straw burning fire point remote-sensing monitoring method based on No. three satellites of wind and cloud, the L1 grade image datas of the visible light infrared scanning radiometer in No. three satellites of wind and cloud are pre-processed first, it is then based on pretreated image data and carries out doubtful fire point extraction, obtain doubtful fire point distribution map, finally according to the land use pattern of fire point pixel and the pseudo- fiery point of the fixed heat source document weeding of inquiry, crop straw burning fire point is isolated on doubtful fire point distribution map;The present invention passes through the visible light infrared scanning radiometer fully studying the characteristic parameter of No. three satellites of wind and cloud and carry thereon(VIRR)Band setting, make it possible the remotely-sensed data using No. three satellites of domestic wind and cloud it is convenient, it is effective, be continuously completed crop straw burning fire point and monitor, and currently used external satellite crop straw burning Monitoring Result can be verified and be supplemented, the time frequency of China's crop straw burning monitoring is greatly improved, effectively enhances stalk and prohibits the ability to supervise burnt.

Description

A kind of crop straw burning fire point remote-sensing monitoring method based on No. three satellites of wind and cloud
Technical field
The present invention relates to crop straw burning fire point monitoring technical field more particularly to a kind of stalks based on No. three satellites of wind and cloud Burn fire point remote-sensing monitoring method.
Background technology
For a long time, crop straw burning prevention and cure of pollution are always one of important process of environmental protection;Open Burning Crop Residue mistake Cheng Zhonghui discharges a large amount of gaseous pollutants and particulate matter, causes centainly to pollute to atmospheric environment;On a large scale, the stalk of high intensity is burnt Burning event can cause the drastically decline of atmospheric visibility, increase haze weather odds, endanger health and the friendship of people Logical safety.
However, information-based monitoring can't be fully achieved mainly based on artificial means to the monitoring of crop straw burning at present Management level;In numerous rural areas of China, government, which mostly uses greatly, prohibits the measure for burning the tour of time value class to peasant's crop straw burning behavior It is supervised, however this examine on the spot will expend a large amount of manpower and materials, and have tunnel vision, be limited in scope, have larger Contingency.
With the high speed development of earth observation technology, the advantages that remote sensing is observed with its Large Area Synchronous, timeliness is strong, become Crop straw burning monitors efficient, economic, practical brand-new route.Due to the daily detection requirement high frequency time of crop straw burning, ensured sustained development, mesh Before be usually used in carrying out the light sensor data of this business mainly with the high grade satellite in orbit number of external hour of temporal resolution Based on Terra/MODIS, Aqua/MODIS, NOAA/AVHRR.However, depending on external remotely-sensed data unduly, there are certain wind Dangerous and limitation lacks independence to the spatio-temporal distribution of data, will greatly influence the progress of monitoring and warning work; In addition, the common satellite generally existing active time such as Terra, Aqua, NOAA is longer, the quality of data declines and Annual distribution is concentrated The problem of.
The transmitting of a new generation of China polar orbiting meteorological satellite wind and cloud No. three (FY-3) brings hope to solve this problem;Wind Cloud three first two experiment stars FY-3A and FY-3B went up to the air respectively at 2008 and 2010,2013, FY-3C transmittings Success, and formally put into service operation in June, 2014;FY-3A/C and FY-3B common networkings at present, form China's SSO (Sun Synchronous Orbit) gas As the business layout that the upper and lower noon star networking of satellite is observed, therefore, how using the remotely-sensed data of No. three satellites of wind and cloud come to stalk It burns fire point and carries out remote sensing monitoring, be a technical problem for being badly in need of solving.
Invention content
The object of the present invention is to provide a kind of crop straw burning fire point remote-sensing monitoring method based on No. three satellites of wind and cloud, can , continuously progress crop straw burning fire point monitoring convenient, effective using the remotely-sensed data of No. three satellites of domestic wind and cloud, effectively enhances straw Stalk prohibits the ability to supervise burnt.
The technical solution adopted by the present invention is:A kind of crop straw burning fire point remote sensing monitoring side based on No. three satellites of wind and cloud Method includes the following steps:
Step A:The data prediction of No. three satellites of wind and cloud;The radiation of visible light infrared scan is loaded on No. three satellites of wind and cloud Meter, the visible light infrared scanning radiometer includes channel one, channel two, channel three, channel four, channel five, channel six, channel Seven, the observation channel of channel eight, channel nine and channel ten totally ten one kilometer of resolution ratio, the channel one, channel two, channel Six, channel seven, channel eight, channel nine and channel ten are visible channel, and channel three, channel four and channel five are infrared channel; The image in channel one, channel two, channel three, channel four, channel five, channel six, channel seven, channel eight, channel nine and channel ten Data have collectively constituted the L1 grade image datas of visible light infrared scanning radiometer, by visible light infrared scanning radiometer L1 grades of image datas pre-processed to obtain the reflectivity image in channel one and channel two, channel three and channel four surface radiation The bright temperature image of the earth's surface in brightness image and channel three.
Step B:Doubtful fire point extraction is carried out based on the image data that step A is pre-processed;Specifically include step B1- B2;
Step B1:RGB pseudo color composings;For image on daytime, the surface radiation brightness image in channel three is assigned to red, The reflectivity image in channel two assigns green, and the reflectivity image in channel one assigns blue, on colored synthesis figure, burning Open fire area be denoted as cerise, cross flame range and be denoted as kermesinus;For night image, by the surface radiation brightness in channel three Image assigns red, and the surface radiation brightness image in channel four assigns green and blue, and open firing point is still cerise in the picture, Other topographical features are difficult identification, because not having visible channel participation in night composite diagram.
Step B2:Into row threshold division on the basis of step B1 carries out image processing;According to artificial visual result to logical The surface radiation brightness image in road three obtains the initial binary that target is the doubtful fire point extraction that 1 background is 0 into row threshold division Figure result;If indicating threshold value with T, which can be described with following formula:
Initial binary figure is compared with false colour composite image, threshold value T is repeatedly finely tuned, closer to most Good threshold value, to obtain doubtful fire point distribution map.
Step C:Fire point identification is carried out on doubtful fire point distribution map and is determined;According to the land use pattern of fire point pixel With the pseudo- fiery point of the fixed heat source document weeding of inquiry, crop straw burning fire point is isolated on doubtful fire point distribution map, stalk is obtained and burns Make a fire a distribution map.
Further the step A further comprises the data prediction of No. three satellites of wind and cloud:
Step A1:Sensor is calibrated;The purpose of the step is to obtain reflectivity image and the channel in channel one and channel two Three and channel four surface radiation brightness image.
First to the image data in channel one and channel two in the L1 grade image datas of visible light infrared scanning radiometer Visible light near-infrared calibration is carried out, the reflectivity image in channel one and channel two is obtained, calibration formula is as follows:
F=SDN+I
Wherein, F is channel albedo, and S is slope, and I is intercept, DN be visible channel earth observation count value, that is, it is bright Angle value;
Secondly to the image data in channel three and channel four in the L1 grade image datas of visible light infrared scanning radiometer Infrared channel calibration is carried out, obtains the surface radiation brightness image in channel three and channel four, infrared channel calibration includes that star is reached the standard grade Property calibration and spoke luminance non-linearity correct two steps.
The formula of linear scaled is on star:
NLIN=ScaleDN+Offset
N in formulaLINFor linear scaled spoke brightness value, unit is m W/m2·cm-1Sr, Scale are gain, and Offset is Intercept, DN are the earth observation count value of infrared channel.
Spoke luminance non-linearity is corrected, i.e. ground calibration, and formula is:
N=b0+(1+b1)NLIN+b2NLIN 2
N is revised calibration spoke brightness value, unit mW/m in formula2·cm-1Sr, i.e. surface radiation brightness value, b0、 b1、b2For correction coefficient.
Step A2:The bright temperature of earth's surface calculates;The purpose of the step is to calculate the bright temperature image of earth's surface in channel three;
Brightness temperature refers to when spectral radiance of the object under Same Wavelength is equal with blackbody spectrum radiation intensity Blackbody temperature;Channel three has obtained surface radiation brightness value by linear scaled and ground calibration on the star in step A1, utilizes Surface radiation brightness value can be converted into effective blackbody temperature by following formula:
In formulaFor effective blackbody temperature, C1=1.1910427 × 10-5mW/(m2·sr·cm-4), C2= 1.4387752cmK VCIt is the infrared channel center wave number that ground calibration obtains, N is the radiance after calibration.
Due in terms of sensor physics, acquired with above formulaMathematics must also be carried out according to sensor characteristics The actual detection temperature of wanted wave band can just be obtained by correcting, and formula is as follows:
T is the bright temperature of earth's surface of wanted wave band in formula, and A, B are constant, the respectively slope and intercept of temperature correction.
Further the step A further comprises step A3 to the data prediction of No. three satellites of wind and cloud;
Step A3:Geometric correction, using GLT geometric corrections method to the obtained channel one step A1 and step A2 and channel two Reflectivity image, channel three and the bright temperature image of the earth's surface in the surface radiation brightness image in channel four and channel three carry out geometry It corrects, the purpose of geometric correction is the geometric distortion for the L1 grade image datas for correcting visible light infrared scanning radiometer.
Pseudo- fire point is rejected according further to the land use pattern of fire point pixel is specially:By doubting in survey region The MCD12Q1 ground mulching categorical datas in L3 grades of products of distribution map and the MODIS in the same year, which are put, like fire carries out space overlapping, it will Fiery vertex type divides into forest fires, Grassland Fire point, urban heat island and crop straw burning fire point, to effectively reject outside crop land Heat anomaly point.
Be further advanced by the fixed heat source data of inquiry is specially to reject pseudo- fire point:Fixed heat source be divided into typical heat source and Two kinds of atypia heat source;Typical heat source is long-term high temperature dot, and main includes large-scale steel plant, steam power plant and Yao Chang, these works Factory possesses heat supply heat-removal equipment, is picked by searching map on the net or establishing data according to other existing fixed heat source data It removes;Atypia heat source is mainly some semiworks of construction near farmland, and distributed more widely and heat releasing source characteristic is unknown, by looking into The correlation experience data for asking long term monitoring accumulation rejects its part.
The present invention passes through the visible light infrared scan spoke fully studying the characteristic parameter of No. three satellites of wind and cloud and carry thereon The band setting of meter (VIRR) is penetrated, first to the L1 grade image datas of the visible light infrared scanning radiometer in No. three satellites of wind and cloud It is pre-processed, is then based on pretreated image data and carries out doubtful fire point extraction, obtain doubtful fire point distribution map, finally According to the land use pattern of fire point pixel and the pseudo- fiery point of the fixed heat source document weeding of inquiry, detached on doubtful fire point distribution map Go out crop straw burning fire point;That the invention enables the remotely-sensed datas using No. three satellites of domestic wind and cloud is convenient, effective, is continuously completed straw Stalk burns fire point monitoring and is possibly realized, and can currently used external satellite crop straw burning Monitoring Result be verified and be mended It fills, the time frequency of China's crop straw burning monitoring is greatly improved, effectively enhance stalk and prohibit the ability to supervise burnt.
Description of the drawings
Fig. 1 is the flow chart of the present invention.
Specific implementation mode
As shown in Figure 1, a kind of crop straw burning fire point remote-sensing monitoring method based on No. three satellites of wind and cloud, including walk as follows Suddenly:
Step A:The data prediction of No. three satellites (FY-3) of wind and cloud;It is loaded with outside visible red and sweeps on No. three satellites of wind and cloud Radiometer (VIRR) is retouched, the visible light infrared scanning radiometer includes channel one, channel two, channel three, channel four, channel Five, the observation channel of channel six, channel seven, channel eight, channel nine and channel ten totally ten one kilometer of resolution ratio, the channel One, channel two, channel six, channel seven, channel eight, channel nine and channel ten are visible channel, channel three, channel four and channel Five be infrared channel, and the centre wavelength in channel three, channel four and channel five is near 4 μm, 11 μm, 12 μm, to underlay surface radiation Variation is sensitive;Channel one, channel two, channel three, channel four, channel five, channel six, channel seven, channel eight, channel nine and channel Ten image data has collectively constituted the L1 grade image datas of visible light infrared scanning radiometer, by visible light infrared scan The L1 grade image datas of radiometer are pre-processed to obtain the reflectivity image in channel one and channel two, channel three and channel four The bright temperature image of the earth's surface in surface radiation brightness image and channel three.
Step A further comprises step A1, step A2 and step A3 to the data prediction of No. three satellites of wind and cloud.
Step A1:Sensor is calibrated;Sensor calibration is actually that the count value that satellite instrument detects is converted into physics Amount, the purpose of the step is the reflectivity image for obtaining channel one and channel two and the surface radiation brightness in channel three and channel four Image.
First to the image data in channel one and channel two in the L1 grade image datas of visible light infrared scanning radiometer Visible light near-infrared calibration is carried out, the reflectivity image in channel one and channel two is obtained, calibration formula is as follows:
F=SDN+I (1)
Wherein, F is channel albedo, and S is slope, and I is intercept, DN be visible channel earth observation count value, that is, it is bright Angle value;The numerical value of S and I is stored in the file attribute of L1 grade image datas of visible light infrared scanning radiometer.
Secondly to the image data in channel three and channel four in the L1 grade image datas of visible light infrared scanning radiometer Infrared channel calibration is carried out, obtains the surface radiation brightness image in channel three and channel four, infrared channel calibration includes that star is reached the standard grade Property calibration and spoke luminance non-linearity correct two steps.
The formula of linear scaled is on star:
NLIN=ScaleDN+Offset (2)
N in formulaLINFor linear scaled spoke brightness value, unit is m W/m2·cm-1Sr, Scale are gain, and Offset is Intercept, DN are the earth observation count value of infrared channel;Scale and Offset is stored in the L1 of visible light infrared scanning radiometer The science data of grade image data are concentrated.
Spoke luminance non-linearity is corrected, i.e. ground calibration, and formula is:
N=b0+(1+b1)NLIN+b2NLIN 2 (3)
N is revised calibration spoke brightness value, unit mW/m in formula2·cm-1Sr, i.e. surface radiation brightness value, b0、 b1、b2For correction coefficient;It is stored in the file attribute of L1 grade image datas of visible light infrared scanning radiometer.
Step A2:The bright temperature of earth's surface calculates;The purpose of the step is that the bright temperature image of earth's surface in channel three is calculated.
Brightness temperature refers to when spectral radiance of the object under Same Wavelength is equal with blackbody spectrum radiation intensity Blackbody temperature;Channel three has obtained surface radiation brightness value by linear scaled and ground calibration on the star in step A1, utilizes Surface radiation brightness value can be converted into effective blackbody temperature by following formula:
In formulaFor effective blackbody temperature, C1=1.1910427 × 10-5mW/(m2·sr·cm-4), C2= 1.4387752cmK VCIt is the infrared channel center wave number that ground calibration obtains, N is the radiance after calibration.
Due in terms of sensor physics, acquired with formula (4)It must also be according to sensor characteristics into line number Correct the actual detection temperature that can just obtain wanted wave band, and formula is as follows:
T is the bright temperature of earth's surface of wanted wave band in formula, and A, B are constant, the respectively slope and intercept of temperature correction.
Step A3:Geometric correction, using GLT geometric corrections method to the obtained channel one step A1 and step A2 and channel two Reflectivity image, channel three and the bright temperature image of the earth's surface in the surface radiation brightness image in channel four and channel three carry out geometry It corrects, the purpose of geometric correction is the geometric distortion for the L1 grade image datas for correcting visible light infrared scanning radiometer;Due to FY- 3VIRR spatial resolutions are low, have no longer been applicable in by the conventional geometric bearing calibration of ground control point, therefore have used GLT geometry school It executes and geometric correction is carried out to VIRR data.
Step B:Doubtful fire point extraction is carried out based on the image data that step A is pre-processed;Specifically include step B1- B2。
Step B1:RGB pseudo color composings;For image on daytime, the surface radiation brightness image in channel three is assigned to red, The reflectivity image in channel two assigns green, and the reflectivity image in channel one assigns blue, on colored synthesis figure, burning Open fire area be denoted as cerise, cross flame range and be denoted as kermesinus;For night orbital data, by the surface radiation in channel three Brightness image assigns red, and the surface radiation brightness image in channel four assigns green and blue, and open firing point is still fresh in the picture Red, other topographical features are difficult identification, because not having visible channel participation in night composite diagram.
Step B2:Into row threshold division on the basis of step B1 carries out pseudo color composing;According to artificial visual result pair The surface radiation brightness image in channel three obtains initial two that target is the doubtful fire point extraction that 1 background is 0 into row threshold division It is worth figure result;If indicating threshold value with T, which can be described with following formula:
Initial binary figure is compared with false colour composite image, threshold value T is repeatedly finely tuned, closer to most Good threshold value, to obtain doubtful fire point distribution map.
Step C:Fire point identification is carried out on doubtful fire point distribution map and is determined;According to the land use pattern of fire point pixel With the pseudo- fiery point of the fixed heat source document weeding of inquiry, crop straw burning fire point is isolated on doubtful fire point distribution map, stalk is obtained and burns Make a fire a distribution map.
Pseudo- fire point is rejected according to the land use pattern of fire point pixel is specially:By the doubtful fire point minute in survey region MCD12Q1 ground mulching categorical datas in L3 grades of products of Butut and the MODIS in the same year carry out space overlapping, by fiery vertex type Forest fires, Grassland Fire point, urban heat island and crop straw burning fire point are divided into, so as to effectively reject the heat outside crop land Abnormal point.
Pseudo- fire point is rejected by the fixed heat source data of inquiry is specially:Fixed heat source is divided into typical heat source and atypia heat Two kinds of source;Typical heat source is long-term high temperature dot, and main includes large-scale steel plant, steam power plant and Yao Chang, these factories possess pot The heat supplies heat-removal equipment such as stove, chimney, the heat discharged daily is relatively high, by searching map on the net or according to other existing Fixed heat source data is established data and is rejected;Atypia heat source is mainly some semiworks of construction near farmland, distribution Unknown compared with wide and heat releasing source characteristic, resulting pseudo- fire point is difficult to accomplish whole rejectings, can be by inquiring long term monitoring accumulation Correlation experience data its part is rejected;If in three months or more time, same position constantly or is frequently supervised Heat anomaly point is measured, which is probably fixed heat source, on-site inspection can be unfolded at this time, if being confirmed as heat source, you can Experiential basis of the correlation experience data as follow-up monitoring is added.
The invention enables simple using domestic No. three satellite remote sensing dates of wind and cloud, quickly and efficiently completion crop straw burning is fiery Point monitoring business is possibly realized, and can currently used external satellite crop straw burning Monitoring Result be verified and be supplemented, greatly The time frequency for improving the monitoring of China's crop straw burning, to avoid flame range farmland to a certain extent because rapid turn over is brought The influence failing to judge and judge by accident of mistake flame range, effectively enhance stalk and prohibit the ability to supervise burnt.
Finally it should be noted that:The above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations;Although Present invention has been described in detail with reference to the aforementioned embodiments, it will be understood by those of ordinary skill in the art that:It is still It can modify to the technical solution recorded in previous embodiment, either which part or all technical features are carried out etc. With replacement;And these modifications or replacements, technical solution of the embodiment of the present invention that it does not separate the essence of the corresponding technical solution Range.

Claims (5)

1. a kind of crop straw burning fire point remote-sensing monitoring method based on No. three satellites of wind and cloud, which is characterized in that include the following steps:
Step A:The data prediction of No. three satellites of wind and cloud;It is loaded with visible light infrared scanning radiometer on No. three satellites of wind and cloud, The visible light infrared scanning radiometer include channel one, channel two, channel three, channel four, channel five, channel six, channel seven, The observation channel of channel eight, channel nine and channel ten totally ten one kilometer of resolution ratio, it is the channel one, channel two, channel six, logical Road seven, channel eight, channel nine and channel ten are visible channel, and channel three, channel four and channel five are infrared channel;Channel One, the image data in channel two, channel three, channel four, channel five, channel six, channel seven, channel eight, channel nine and channel ten The L1 grade image datas for having collectively constituted visible light infrared scanning radiometer, pass through the L1 grades to visible light infrared scanning radiometer Image data is pre-processed to obtain the reflectivity image in channel one and channel two, channel three and the surface radiation brightness in channel four The bright temperature image of the earth's surface in image and channel three;
Step B:Doubtful fire point extraction is carried out based on the image data that step A is pre-processed;Specifically include step B1-B2;
Step B1:RGB pseudo color composings;For image on daytime, the surface radiation brightness image in channel three is assigned to red, channel Two reflectivity image assigns green, and the reflectivity image in channel one assigns blue, and on colored synthesis figure, what is burnt is bright Flame range is denoted as cerise, crosses flame range and is denoted as kermesinus;For night image, by the surface radiation brightness image in channel three Red is assigned, the surface radiation brightness image in channel four assigns green and blue, and open firing point is still cerise in the picture, other Topographical features are difficult identification, because not having visible channel participation in night composite diagram;
Step B2:Into row threshold division on the basis of step B1 carries out image processing;According to artificial visual result to channel three Surface radiation brightness image into row threshold division, obtain the initial binary figure knot that target is the doubtful fire point extraction that 1 background is 0 Fruit;If indicating threshold value with T, which can be described with following formula:
Initial binary figure is compared with false colour composite image, threshold value T is repeatedly finely tuned, closer to best threshold Value, to obtain doubtful fire point distribution map;
Step C:Fire point identification is carried out on doubtful fire point distribution map and is determined;According to the land use pattern of fire point pixel and look into The pseudo- fiery point of fixed heat source document weeding is ask, crop straw burning fire point is isolated on doubtful fire point distribution map, obtains crop straw burning fire Point distribution map.
2. the crop straw burning fire point remote-sensing monitoring method based on No. three satellites of wind and cloud as described in claim 1, it is characterised in that: The step A further comprises the data prediction of No. three satellites of wind and cloud:
Step A1:Sensor is calibrated;The purpose of the step is to obtain three He of reflectivity image and channel in channel one and channel two The surface radiation brightness image in channel four;
The image data in channel one and channel two in the L1 grade image datas of visible light infrared scanning radiometer is carried out first Visible light near-infrared is calibrated, and the reflectivity image in channel one and channel two is obtained, and calibration formula is as follows:
F=SDN+I
Wherein, F is channel albedo, and S is slope, and I is intercept, and DN is earth observation count value, that is, brightness of visible channel Value;
Secondly the image data in channel three and channel four in the L1 grade image datas of visible light infrared scanning radiometer is carried out Infrared channel is calibrated, and the surface radiation brightness image in channel three and channel four is obtained, and infrared channel calibration includes linearly fixed on star Mark and spoke luminance non-linearity correct two steps;
The formula of linear scaled is on star:
NLIN=ScaleDN+Offset
N in formulaLINFor linear scaled spoke brightness value, unit is m W/m2·cm-1Sr, Scale are gain, and Offset is intercept, DN is the earth observation count value of infrared channel;
Spoke luminance non-linearity is corrected, i.e. ground calibration, and formula is:
N=b0+(1+b1)NLIN+b2NLIN 2
N is revised calibration spoke brightness value, unit mW/m in formula2·cm-1Sr, i.e. surface radiation brightness value, b0、b1、b2 For correction coefficient;
Step A2:The bright temperature of earth's surface calculates;The purpose of the step is to calculate the bright temperature image of earth's surface in channel three;
Brightness temperature refers to black matrix when spectral radiance of the object under Same Wavelength is equal with blackbody spectrum radiation intensity Temperature;Channel three has obtained surface radiation brightness value, has utilized following formula by linear scaled and ground calibration on the star in step A1 Surface radiation brightness value can be converted into effective blackbody temperature:
In formulaFor effective blackbody temperature, C1=1.1910427 × 10-5mW/(m2·sr·cm-4), C2=1.4387752cm K, VCIt is the infrared channel center wave number that ground calibration obtains, N is the radiance after calibration;
Due in terms of sensor physics, acquired with above formulaMathematics must also be carried out according to sensor characteristics to correct The actual detection temperature of wanted wave band can be obtained, formula is as follows:
T is the bright temperature of earth's surface of wanted wave band in formula, and A, B are constant, the respectively slope and intercept of temperature correction.
3. the crop straw burning fire point remote-sensing monitoring method based on No. three satellites of wind and cloud as claimed in claim 2, it is characterised in that: The step A further comprises step A3 to the data prediction of No. three satellites of wind and cloud;
Step A3:Geometric correction, using GLT geometric corrections method to the anti-of the obtained channel one step A1 and step A2 and channel two The bright temperature image of earth's surface of the surface radiation brightness image and channel three of penetrating rate image, channel three and channel four carries out geometry and entangles Just, the purpose of geometric correction is the geometric distortion for the L1 grade image datas for correcting visible light infrared scanning radiometer.
4. the crop straw burning fire point remote-sensing monitoring method based on No. three satellites of wind and cloud, feature exist as claimed in claim 1 or 2 In:Pseudo- fire point is rejected according to the land use pattern of fire point pixel is specially:By the doubtful fire point distribution map in survey region Space overlapping is carried out with the MCD12Q1 ground mulching categorical datas in the L3 grades of products of MODIS in the same year, fiery vertex type is distinguished For forest fires, Grassland Fire point, urban heat island and crop straw burning fire point, to effectively reject the heat anomaly point outside crop land.
5. the crop straw burning fire point remote-sensing monitoring method based on No. three satellites of wind and cloud, feature exist as claimed in claim 1 or 2 In:Pseudo- fire point is rejected by the fixed heat source data of inquiry is specially:Fixed heat source is divided into typical heat source and atypia heat source two Kind;Typical heat source is long-term high temperature dot, and main includes large-scale steel plant, steam power plant and Yao Chang, these factories possess heat supply row Hot equipment is rejected by searching map on the net or establishing data according to other existing fixed heat source data;Atypia heat Source is mainly some semiworks of construction near farmland, and distributed more widely and heat releasing source characteristic is unknown, passes through and inquires long term monitoring The correlation experience data of accumulation rejects its part.
CN201711471773.8A 2017-12-29 2017-12-29 A kind of crop straw burning fire point remote-sensing monitoring method based on No. three satellites of wind and cloud Pending CN108304780A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711471773.8A CN108304780A (en) 2017-12-29 2017-12-29 A kind of crop straw burning fire point remote-sensing monitoring method based on No. three satellites of wind and cloud

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711471773.8A CN108304780A (en) 2017-12-29 2017-12-29 A kind of crop straw burning fire point remote-sensing monitoring method based on No. three satellites of wind and cloud

Publications (1)

Publication Number Publication Date
CN108304780A true CN108304780A (en) 2018-07-20

Family

ID=62867628

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711471773.8A Pending CN108304780A (en) 2017-12-29 2017-12-29 A kind of crop straw burning fire point remote-sensing monitoring method based on No. three satellites of wind and cloud

Country Status (1)

Country Link
CN (1) CN108304780A (en)

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109536231A (en) * 2018-12-21 2019-03-29 王武荣 Agricultural stalk and forestry waste comprehensively utilize monitoring system
CN110348421A (en) * 2019-07-18 2019-10-18 北华航天工业学院 Temperate Region in China forest fire detection method based on MODIS data
CN111008599A (en) * 2019-12-06 2020-04-14 国网湖南省电力有限公司 Mountain fire synchronous satellite monitoring and early warning method and system based on time scale feature evaluation
CN111380807A (en) * 2020-03-30 2020-07-07 中国科学院东北地理与农业生态研究所 Straw burning fire point information extraction method based on geostationary satellite remote sensing data
CN111696074A (en) * 2020-04-30 2020-09-22 中国资源卫星应用中心 Fire point monitoring system based on high score four and application
CN111753754A (en) * 2020-06-28 2020-10-09 三亚中科遥感研究所 Straw combustion fire point identification method based on heat source heavy industry area analysis
CN111858813A (en) * 2020-07-21 2020-10-30 云南电网有限责任公司带电作业分公司 Non-fire area eliminating method based on satellite technology
CN112834040A (en) * 2020-12-24 2021-05-25 合肥冬行信息科技有限公司 Ground thermal anomaly identification algorithm based on stationary satellite
CN112883857A (en) * 2021-02-05 2021-06-01 南京恩瑞特实业有限公司 Method for recognizing low-cloud large-fog-drop area by NRIET satellite
CN113177496A (en) * 2021-05-10 2021-07-27 浙江大华技术股份有限公司 Fire point detection method, device, equipment and storage medium
CN113657275A (en) * 2021-08-16 2021-11-16 中国科学院空天信息创新研究院 Automatic detection method for forest and grass fire points
CN113674166A (en) * 2021-07-30 2021-11-19 中国环境科学研究院 Accounting method for methane emission generated by waste incineration treatment
CN114689177A (en) * 2022-05-30 2022-07-01 安徽气象信息有限公司 System for effectively improving satellite detection fire point
CN116363047A (en) * 2022-08-23 2023-06-30 生态环境部卫星环境应用中心 Early warning method for pollution risk of straw incineration atmosphere
CN117789025A (en) * 2023-12-26 2024-03-29 北京工业大学 Open-air incineration emission estimation method for crop straws based on climate information

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104240429A (en) * 2014-08-27 2014-12-24 中国资源卫星应用中心 Method for monitoring forest fire by means of combination of infrared and multi-spectral wave bands
KR101484667B1 (en) * 2014-07-23 2015-01-27 대한민국 A method and system for fire detection using satellite observation

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101484667B1 (en) * 2014-07-23 2015-01-27 대한민국 A method and system for fire detection using satellite observation
CN104240429A (en) * 2014-08-27 2014-12-24 中国资源卫星应用中心 Method for monitoring forest fire by means of combination of infrared and multi-spectral wave bands

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
闫厚等: "气象卫星在森林火灾监测中的应用及其原理", 《森林防火》 *
陈鹏狮等: "辽宁森林火灾的遥感监测", 《辽宁气象》 *
马骏等: "基于多层次并行的风云三号气象卫星数据定标", 《微电子学与计算机》 *

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109536231A (en) * 2018-12-21 2019-03-29 王武荣 Agricultural stalk and forestry waste comprehensively utilize monitoring system
CN110348421A (en) * 2019-07-18 2019-10-18 北华航天工业学院 Temperate Region in China forest fire detection method based on MODIS data
CN111008599A (en) * 2019-12-06 2020-04-14 国网湖南省电力有限公司 Mountain fire synchronous satellite monitoring and early warning method and system based on time scale feature evaluation
CN111380807A (en) * 2020-03-30 2020-07-07 中国科学院东北地理与农业生态研究所 Straw burning fire point information extraction method based on geostationary satellite remote sensing data
CN111696074A (en) * 2020-04-30 2020-09-22 中国资源卫星应用中心 Fire point monitoring system based on high score four and application
CN111696074B (en) * 2020-04-30 2024-04-05 中国资源卫星应用中心 Fire monitoring system based on high-resolution fourth-order and application
CN111753754B (en) * 2020-06-28 2023-09-12 三亚中科遥感研究所 Straw burning fire point identification method based on heat source heavy industry area analysis
CN111753754A (en) * 2020-06-28 2020-10-09 三亚中科遥感研究所 Straw combustion fire point identification method based on heat source heavy industry area analysis
CN111858813A (en) * 2020-07-21 2020-10-30 云南电网有限责任公司带电作业分公司 Non-fire area eliminating method based on satellite technology
CN112834040A (en) * 2020-12-24 2021-05-25 合肥冬行信息科技有限公司 Ground thermal anomaly identification algorithm based on stationary satellite
CN112883857A (en) * 2021-02-05 2021-06-01 南京恩瑞特实业有限公司 Method for recognizing low-cloud large-fog-drop area by NRIET satellite
CN113177496A (en) * 2021-05-10 2021-07-27 浙江大华技术股份有限公司 Fire point detection method, device, equipment and storage medium
CN113674166A (en) * 2021-07-30 2021-11-19 中国环境科学研究院 Accounting method for methane emission generated by waste incineration treatment
CN113657275A (en) * 2021-08-16 2021-11-16 中国科学院空天信息创新研究院 Automatic detection method for forest and grass fire points
CN114689177A (en) * 2022-05-30 2022-07-01 安徽气象信息有限公司 System for effectively improving satellite detection fire point
CN116363047A (en) * 2022-08-23 2023-06-30 生态环境部卫星环境应用中心 Early warning method for pollution risk of straw incineration atmosphere
CN116363047B (en) * 2022-08-23 2024-01-30 生态环境部卫星环境应用中心 Early warning method for pollution risk of straw incineration atmosphere
CN117789025A (en) * 2023-12-26 2024-03-29 北京工业大学 Open-air incineration emission estimation method for crop straws based on climate information
CN117789025B (en) * 2023-12-26 2024-06-04 北京工业大学 Open-air incineration emission estimation method for crop straws based on climate information

Similar Documents

Publication Publication Date Title
CN108304780A (en) A kind of crop straw burning fire point remote-sensing monitoring method based on No. three satellites of wind and cloud
Wooster et al. Satellite remote sensing of active fires: History and current status, applications and future requirements
Xu et al. Major advances in geostationary fire radiative power (FRP) retrieval over Asia and Australia stemming from use of Himarawi-8 AHI
Rozenstein et al. Estimating cotton water consumption using a time series of Sentinel-2 imagery
Ichoku et al. Comparative analysis of daytime fire detection algorithms using AVHRR data for the 1995 fire season in Canada: Perspective for MODIS
Feng et al. Comparison of valid ocean observations between MODIS Terra and Aqua over the global oceans
CN102798470B (en) Method and device for monitoring range of heat abnormality in mining area
CN104615848A (en) Comprehensive application method for integrating forest fire danger forecasting and ground fire point monitoring
CN111505016B (en) Method for estimating straw burning discharge list
Dickinson et al. Measuring radiant emissions from entire prescribed fires with ground, airborne and satellite sensors–RxCADRE 2012
CN112102578B (en) Forest fire monitoring and early warning system, method, storage medium and computer equipment
CN113324656A (en) Unmanned aerial vehicle-mounted infrared remote sensing earth surface heat anomaly detection method and system
Pu et al. A dynamic algorithm for wildfire mapping with NOAA/AVHRR data
CN113221057A (en) Straw burning fire point monitoring method based on multi-temporal satellite images
CN112183319A (en) Forest and grassland fire behavior discrimination method based on remote sensing image and vegetation type map
Galindo et al. Real-time AVHRR forest fire detection in Mexico (1998-2000)
CN105096511A (en) Method for generating smoke area detection result automatically
Chen et al. Spatial distribution and pattern persistence of surface soil moisture and temperature over prairie from remote sensing
Halle et al. Firebird―Small Satellites for Wild Fire Assessment
Frantzova Remote sensing, GIS and Disaster management
Zhang et al. GF-4 Satellite Fire Detection with an improved Contextual algorithm
Ananth et al. Mapping of Burnt area and Burnt Severity using Landsat 8 Images: A Case Study of Bandipur forest Fire Region of Karnataka state India
Ismail et al. Forest fires hotspot detection in Indonesia using Himawari-8 satellite data
Bag et al. Chapter-4 Real-Time Monitoring of Crop Residue Burning through Satellite Remote Sensing
Kato et al. Detection of thermal anomaly using Sentinel-2A data

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
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20180720