CN115830472A - Urban complex underlying surface methane emission source identification method - Google Patents
Urban complex underlying surface methane emission source identification method Download PDFInfo
- Publication number
- CN115830472A CN115830472A CN202310040669.2A CN202310040669A CN115830472A CN 115830472 A CN115830472 A CN 115830472A CN 202310040669 A CN202310040669 A CN 202310040669A CN 115830472 A CN115830472 A CN 115830472A
- Authority
- CN
- China
- Prior art keywords
- methane
- data
- ratio
- waveband
- methane sensitive
- 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.)
- Granted
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a method for identifying a methane emission source of an urban complex underlying surface, which comprises the following steps: acquiring high-resolution satellite data of a target area; fitting the high-resolution satellite data to obtain reference methane sensitive waveband data; acquiring the ratio of the reference methane sensitive waveband data to the methane sensitive waveband data; and acquiring a methane emission source of the target area according to the ratio. The invention realizes the identification of the high space-time resolution of the methane emission source of the urban complex underlying surface through the characteristic wave band, and provides reference for the supervision of the methane emission source.
Description
Technical Field
The invention relates to the field of atmospheric environment air quality monitoring. In particular to a method for identifying a methane emission source of an urban complex underlying surface.
Background
Methane (CH) 4 ) As second only to CO 2 The second greenhouse gas of (2) has a contribution rate of about one fourth to global warming. It is estimated that the worldwide artificial methane emission is predicted to be 9390 million metric tons of CO in 2020 2 Equivalent, global methane action (GMI) locks about 54% of its sources of emissions, divided into five parts, including: agriculture (fecal management), coal mines, municipal solid waste, oil and gas systems, and sewage. Ergonomic CaliforniaThe institute jet propulsion laboratory scientist of Ruili Doen and colleagues research team has investigated more than 272000 infrastructure points, and through detecting, positioning and quantifying 564 outstanding methane point sources, the biggest source of methane point source emission is landfill, and then dairy and oil gas industries. It shows that although the monomer emission is small, the monomer emission has huge quantity base number and has considerable and significant contribution to methane emission, so that in addition to monitoring super emission sources with large emission, various types of CH with large quantity and scattered distribution 4 Emissions sources also need to be monitored, and identifying and supervising such emissions sources is of great significance for methane abatement.
Disclosure of Invention
The invention is provided based on the above requirements of the prior art, and mainly aims at solving the problem that the methane emission source cannot be efficiently, accurately and conveniently identified in the urban complex underlying surface area, and provides a method for identifying the methane emission source of the urban complex underlying surface.
In order to solve the problems, the invention is realized by adopting the following technical scheme:
a method for identifying a methane emission source of an urban complex underlying surface comprises the following steps:
acquiring high-resolution satellite data of a target area;
fitting the high-resolution satellite data to obtain reference methane sensitive band data;
acquiring the ratio of the reference methane sensitive waveband data to the methane sensitive waveband data;
and acquiring a methane emission source of the target area according to the ratio.
According to the design, the high-resolution satellite data is obtained, the methane sensitive waveband data and the non-methane sensitive waveband data are fitted to obtain the relation between the methane sensitive waveband data and the non-methane sensitive waveband data, data which do not meet the relation in the high-resolution satellite data are screened out, and the corresponding position, namely the area where the methane emission source is located, is found. The high space-time resolution of the methane emission source of the urban complex underlying surface is identified by setting the reference methane sensitive waveband as the characteristic waveband.
Optionally, before said fitting said high resolution satellite data to obtain reference methane sensitive band data, said method further comprises:
and based on a methane emission source database, removing the high-resolution satellite data of the area influenced by the methane emission source to obtain the removed high-resolution satellite data.
With the design, the fitting precision is improved through the background area after the methane emission source is removed, so that the identification result is more accurate.
Optionally, before said fitting said high resolution satellite data to obtain reference methane sensitive band data, said method further comprises:
and randomly screening the high-resolution satellite data to obtain the screened high-resolution satellite data.
With this design, the data volume can be reduced, and the computational efficiency is improved.
Optionally, said fitting said high resolution satellite data to obtain reference methane sensitive band data comprises:
performing regression fitting on the methane sensitive waveband data through the non-methane sensitive waveband data to determine a fitting function;
and obtaining reference methane sensitive waveband data according to the fitting function and the non-methane sensitive waveband data.
The regression fitting effect is good, so that the obtained fitting function parameter precision is higher, and the identification accuracy is improved.
Optionally, the obtaining a ratio of the reference methane sensitive band data to the methane sensitive band data includes:
and carrying out ratio operation on the reference methane sensitive waveband data based on the same geographical position information and the corresponding methane sensitive waveband data to obtain the ratio of the reference methane sensitive waveband data to the methane sensitive waveband data.
Optionally, the obtaining a ratio of the reference methane sensitive band data to the methane sensitive band data includes:
when the target area is positioned in the same scene image, determining the ratio according to the reference methane sensitive waveband representation reflectivity in the reference methane sensitive waveband data and the methane sensitive waveband representation reflectivity in the methane sensitive waveband data;
or
And determining the ratio according to the reference methane sensitive waveband radiance in the reference methane sensitive waveband data and the methane sensitive waveband radiance in the methane sensitive waveband data.
With the design, as the emission increment of methane in the methane emission area is increased, the radiance of a methane sensitive waveband is lower than that of the methane emission area without methane emission, so that the radiance variation and the emission increment of methane have a negative correlation relationship, and the methane emission area can be obtained through more accurate screening; when the target area is positioned in the same scene image, the radiance and the representation reflectivity of each position can be approximately in the same coefficient relation, and the data of the band representation reflectivity can be easily acquired, so that the calculation process can be reduced, and the calculation efficiency can be improved.
Optionally, the obtaining the methane emission source of the target area according to the ratio includes:
if the ratio is less than or equal to 1, no methane is discharged from the corresponding position in the target area;
if the ratio is larger than 1, methane is discharged from the corresponding position in the target area;
screening to obtain a methane emission source in the target area with the ratio being greater than 1;
or
Processing the ratio by using a median filtering or deep learning algorithm to obtain a high-value area;
and screening to obtain a methane emission source according to the high-value area.
Optionally, the formula for determining the ratio according to the reference methane sensitive waveband radiance in the reference methane sensitive waveband data and the methane sensitive waveband radiance in the methane sensitive waveband data is as follows:
r= F(B ref ) / F(B)
wherein r is the ratio of the radiance of a reference methane sensitive waveband to the radiance of the methane sensitive waveband; f (B) ref ) Is a reference methane sensitive band radiance; f (B) is the radiance of the methane sensitive waveband.
Optionally, before acquiring the methane emission source of the target area according to the ratio, rejecting the water body and the shadow area in the target area comprises:
and removing the water body with interference according to the analyzed water body index, and removing the shadow area by using a low-value area of a visible light wave band.
With this design, avoid the water because the characteristics of absorbing partial methane sensitive wave band, influence final result, improve the precision.
Optionally, the method further comprises:
displaying the high-value area on a high-definition map through a first corresponding relation to obtain a methane emission source of the target area; the first corresponding relation is the longitude and latitude relation between the longitude and latitude of the high-value area and the longitude and latitude of the corresponding high-definition map.
With this design, can more audio-visually see the high value region through the image show, compare the high definition map that corresponds the position and obtain latent methane emission source, let the user have better experience and feel.
Compared with the prior art, the method utilizes the high-resolution satellite data, obtains the relationship between the methane sensitive waveband data and the non-methane sensitive waveband data by fitting, screens out the data which do not meet the relationship in the high-resolution satellite data, and finds out the corresponding position, namely the area where the methane emission source is located. The high space-time resolution of the methane emission source of the urban complex underlying surface is identified by setting the reference methane sensitive waveband as the characteristic waveband, and a reference is provided for supervision of the methane emission source.
Drawings
In order to more clearly illustrate the embodiments of the present specification or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments described in the embodiments of the present specification, and other drawings can be obtained by those skilled in the art according to the drawings.
FIG. 1 is a flow chart of a method for quickly identifying a methane emission source of an urban complex underlying surface according to an embodiment of the present invention;
FIG. 2 is a diagram of image data of a Fengcai area 2021 year 1 month 17 day of a method for rapidly identifying a methane emission source of an urban complex underlying surface according to an embodiment of the present invention and random sampling points;
FIG. 3 is a diagram of a Fennel zone B12 of a method for identifying a source of methane emission from a complex underlying surface of a city according to an embodiment of the present invention ref A schematic distribution diagram of/B12;
FIG. 4 shows the identification result of a typical landfill I in a plateau area of a complex urban underlying surface methane emission source identification method according to an embodiment of the present invention;
fig. 5 shows the recognition result of a typical landfill II in a plateau area of the method for recognizing a methane emission source of a complex underlying surface of a city according to the embodiment of the present invention;
fig. 6 shows the identification result of a typical sewage treatment plant in a pool area of the identification method for the urban complex underlying surface methane emission source according to the embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
For the convenience of understanding of the embodiments of the present invention, the following description will be further explained with reference to specific embodiments, which are not to be construed as limiting the embodiments of the present invention.
Example 1
The traditional methane emission source mostly depends on statistical data of 'from bottom to top and from top to bottom', annual statistical data and manual investigation work are taken as main means, and huge manpower and material resources are consumed for dynamic tracking, defect checking and leakage repairing of the emission source, so that aiming at the problem which is not solved by the prior art, the embodiment provides the method for identifying the methane emission source of the complex underlying surface of the city, the flow of the method is shown in figure 1 and comprises the following steps:
s1: high resolution satellite data is acquired for a target area.
The target zone comprises a methane-emitting zone and a non-methane-emitting zone; the high resolution satellite data comprises methane sensitive band data and non-methane sensitive band data; that is, each region includes methane sensitive band data and non-methane sensitive band data. The methane sensitive waveband data comprises the radiance of a methane sensitive waveband or the apparent reflectivity of the methane sensitive waveband; the non-methane sensitive waveband data comprise non-methane sensitive waveband radiance or non-methane sensitive waveband expression reflectivity.
The embodiment of the invention obtains the data of the high-resolution satellite by using the high-resolution satellite with the characteristic wave band in the short wave infrared range, wherein the satellite can be a sentinel II satellite, a Landsat 8 satellite or a world view-3 satellite.
The formula for calculating the band radiance is as follows:
L λ =Gain*DN+Offset
wherein L is λ The equivalent radiance at the entrance pupil of the satellite loading channel is shown, gain is a calibration coefficient, DN is the pixel brightness value of the remote sensing image, and Offset is the Gain Offset.
The calculation formula of the band expressed reflectivity is as follows:
wherein ρ λ Reflectance, L, is expressed for the wavelength band λ Equivalent radiance at the entrance pupil of the satellite loading channel, d is the distance between the sun and the earth, ESUN λ Is the solar irradiance, theta is the solar altitude.
Preferably, before said fitting said high resolution satellite data to obtain reference methane sensitive band data, said method further comprises:
and based on a methane emission source database, removing the high-resolution satellite data of the area influenced by the methane emission source to obtain the removed high-resolution satellite data.
The methane emission source database is obtained through public data, the public data comprise POI data, data in the public database and government official document information, and high-resolution satellite data after elimination of oil and gas facilities, livestock and poultry farms, paddy fields, coal mines, sewage treatment plants, waste gas treatment facilities and the like which possibly relate to methane emission areas are obtained on the basis of the methane emission source database.
For the areas with no methane emission or almost negligible extremely low methane emission, such as buildings, in most cities, the areas can be directly used as background areas irrelevant to methane emission, and the high-resolution satellite data are fitted; preferably, in order to improve the accuracy of final identification, the fitting is performed after the area affected by the methane emission source is removed.
For the case that the methane emission source is widely distributed in cities, for example, a target city is a typical coal mine city or a service chemical city, the area affected by the methane emission source needs to be eliminated.
Preferably, before said fitting said high resolution satellite data to obtain reference methane sensitive band data, said method further comprises:
and randomly screening the high-resolution satellite data to obtain the screened high-resolution satellite data.
The specific process of the random screening in the step is as follows:
directly and randomly selecting n sample points after acquiring high-resolution satellite data of a target area to obtain high-resolution satellite data of the screened sample points;
or
And after the high-resolution satellite data after being removed is obtained, randomly selecting n sample points to obtain the high-resolution satellite data of the sample points after being screened.
In the target area, the point n is preferably greater than 5, and the greater the point n is, the more the effect of eliminating the point of the small methane discharge area which may be mixed, and the more the point n is preferably equal to or greater than 300.
S2: and fitting the high-resolution satellite data to obtain reference methane sensitive band data.
Preferably, said obtaining reference methane sensitive band data by fitting said high resolution satellite data comprises:
carrying out regression fitting on the methane sensitive waveband data through the non-methane sensitive waveband data to determine a fitting function;
and obtaining reference methane sensitive waveband data according to the fitting function and the non-methane sensitive waveband data.
And selecting the parameter of the function relation which best corresponds to the fitting result in the regression fitting process as a final fitting function. The non-methane sensitive band data may be the non-methane sensitive band data in the high resolution satellite data, the non-methane sensitive band data in the high resolution satellite data after being rejected, or the non-methane sensitive band data in the high resolution satellite data after being screened. The methane sensitive waveband data are methane sensitive waveband data corresponding to the non-methane sensitive waveband data based on the same geographical position information.
Preferably, the non-methane sensitive waveband data in the high-resolution satellite data is substituted into a function formula with the best fitting result, and the reference methane sensitive waveband data is obtained through calculation.
Through the regression fitting mode, the effect of establishing the relation between the two is good, and the accuracy of the result is improved.
The reference methane sensitive waveband data refers to methane sensitive waveband data of a region with low methane emission or no methane emission, which is obtained by using the fitting function and the non-methane sensitive waveband data.
According to the characteristic that methane has characteristic absorption in short wave infrared (about 1680nm and about 2300 nm) and the short wave infrared band range of a high-resolution satellite, the methane sensitive band corresponding to the satellite is determined, for example, the methane sensitive band is represented by B11 and B12 in a sentinel second satellite, B6 and B7 in Lansat 8 and B7 in WorldView-3, namely the methane sensitive band can be used as the methane emission sensitive band. And performing regression fitting on the methane sensitive waveband data by using the non-methane sensitive waveband data, and selecting a function relation formula which best corresponds to a fitting result.
The wavelength of the short-wave infrared mid-band 11 of the sentinel II is between-1560 and 1660 nanometers, the wavelength of the band 12 is between-2090 and 2290 nanometers, the pixel resolution of 20 meters is achieved, and according to the characteristic that the short-wave infrared mid-band has characteristic absorption around 1680nm and 2300nm, methane can absorb the SWIR band 11 and the band 12 of the methane sentinel II to a great extent, so that the two bands are selected as methane sensitive bands of the sentinel II.
When at least two sensitive wave bands for methane emission exist, establishing data relation between each methane sensitive wave band and non-methane sensitive wave band data respectively, fitting to obtain a plurality of functional relations, comparing fitting results, selecting the functional relation with the maximum correlation as the functional relation for calculation of the embodiment of the invention, wherein the methane sensitive wave band data corresponding to the functional relation are used for calculating the ratio of the reference methane sensitive wave band data to the methane sensitive wave band data in the next step, for example, the methane sensitive wave bands of sentinel second number data are B11 and B12, and other 10 wave bands (B1, B2, B3, B4, B5, B6, B7, B8 and B8A, B) are respectively fitted to B11 and B12, and the fitting result (correlation) of the B12 is superior to the B11 after the test, so the functional relation of the B12 is selected, and the B12 data is used for calculating the ratio of the reference methane sensitive wave band data to the methane sensitive wave band data.
S3: and acquiring the ratio of the reference methane sensitive waveband data to the methane sensitive waveband data.
Preferably, in this step, the method comprises:
and carrying out ratio operation on the reference methane sensitive waveband data based on the same geographical position information and the corresponding methane sensitive waveband data to obtain the ratio of the reference methane sensitive waveband data to the methane sensitive waveband data.
Preferably, the obtaining the ratio of the reference methane sensitive band data to the methane sensitive band data includes:
when the target area is positioned in the same scene image, determining the ratio according to the reference methane sensitive waveband representation reflectivity in the reference methane sensitive waveband data and the methane sensitive waveband representation reflectivity in the methane sensitive waveband data;
or
And determining the ratio according to the reference methane sensitive waveband radiance in the reference methane sensitive waveband data and the methane sensitive waveband radiance in the methane sensitive waveband data.
The principle of identifying methane by using the radiance ratio is that the radiance increment of a methane emission area is increased, so that the radiance of a methane sensitive waveband is lower than that of the methane emission area when no methane is emitted, so that the radiance variation and the radiance increment have a negative correlation, therefore, the ratio of the radiance of a reference methane sensitive waveband and the radiance of the methane sensitive waveband can represent the methane emission increment to a certain extent, and the formula for determining the ratio according to the radiance of the reference methane sensitive waveband in the reference methane sensitive waveband data and the radiance of the methane sensitive waveband in the methane sensitive waveband data is as follows:
r= F(B ref ) / F(B)
wherein r is the ratio of the radiance of a reference methane sensitive waveband to the radiance of the methane sensitive waveband; f (B) ref ) Reference methane sensitive band radiance; f (B) is the radiance of the methane sensitive waveband.
The same scene image refers to the range covered by one-time imaging completed by the satellite, and is determined by the satellite, and the coverage areas of different satellites are different.
Because the conversion relation exists between the apparent reflectivity and the radiance, when the target area is in the same scene image, the radiance and the expressed reflectivity of each position in the target area can be approximately regarded as rho to k.L for the same scene image, wherein L is the radiance, rho is the apparent reflectivity, and k is a coefficient, namely the radiance and the expressed reflectivity are approximately in a linear relation.
With the design, when the target area is positioned in the same scene image, because the data of the band representation reflectivity is easier to obtain, the calculation process can be reduced, and the calculation efficiency is improved.
S4: and acquiring a methane emission source of the target area according to the ratio.
Preferably, in this step, specifically including:
if the ratio is less than or equal to 1, no methane is discharged from the corresponding position in the target area;
if the ratio is larger than 1, methane gas exists at the corresponding position in the target area;
screening to obtain a methane emission source in the target area with the ratio being greater than 1;
or
Processing the ratio by using a median filtering or deep learning algorithm to obtain a high-value area;
and screening to obtain a methane emission source according to the high-value area.
The principle of determining whether methane is discharged or not by comparing the ratio with 1 is that for places without methane discharge, the non-methane sensitive waveband data is brought into a function obtained by fitting, the obtained result is the same as the corresponding methane sensitive waveband data, so that the ratio is equal to 1, or an error exists, the result is smaller than the corresponding methane sensitive waveband data, and the ratio is smaller than 1; for places with methane emission, because the increment of methane emission is in a negative correlation relation with the change of radiance, the data of the methane sensitive waveband is smaller than the data of the reference methane sensitive waveband obtained by calculation, the ratio is larger than 1, and in an area with the ratio larger than 1, a methane emission source can be obtained by selecting the position with the largest ratio.
In addition, the ratio of the reference methane sensitive waveband data to the methane sensitive waveband data can be displayed in an image mode, a high-value area of the ratio of the reference methane sensitive waveband data to the methane sensitive waveband data is identified through median filtering or deep learning, and a methane emission source is obtained through screening according to the high-value area.
Preferably, the method further comprises the following steps:
displaying the high-value area on a high-definition map through a first corresponding relation to obtain a methane emission source of the target area; the first corresponding relation is the longitude and latitude of the high-value area and the longitude and latitude of the corresponding high-definition map.
With this design, can more audio-visual high value region of seeing through the image show, obtain latent methane emission source to the high definition map of corresponding position, let the user have better experience sense.
Preferably, before acquiring the methane emission source of the target area according to the ratio, rejecting the water body and the shadow area in the target area comprises:
and removing the water body with interference according to the analyzed water body index, and removing the shadow area by using a low-value area of a visible light wave band.
As the water body also has absorption in the methane sensitive wave band, for example, B11/B12 of sentinel II, B6/B7 of lansat 8 and B7/B8 of WorldView-3 are partially absorbed by the water body according to the wave band range of a high-resolution satellite. Therefore, for the analysis result, certain interference exists in the water body and the shadow area, the interference water body and the shadow area are removed by analyzing the water body index (NDVI or MNDWI) and the low-value area (shadow) of the visible light wave band, and then the high-value distribution area is extracted from the result of removing the interference area, so that a more accurate methane emission source can be obtained.
By utilizing the provided method for identifying the methane emission source of the complex underlying surface of the city, the embodiment of the invention provides a specific embodiment, taking the data of a sentinel second number in the Fennel district in 1 month and 17 days of 2021 as an example, two solid waste landfill sites and one methane emission point of a sewage treatment plant in the Fennel district can be clearly identified through the wave band ratio, and the specific steps are as follows:
step one, acquiring high-resolution satellite cloud removing data of a platform area
And (3) acquiring apparent reflectivity data of a sentinel second satellite which is removed from the cloud in the Fengcai region 2021, 1, 17 days, and resampling all wave band data to 10 m resolution. And randomly selecting 1000 point positions from the resampled band data to extract 12 bands (B1, B2, B3, B4, B5, B6, B7, B8A, B, B11 and B12) of the sentinel II.
Step two, acquiring a reference methane sensitive wave band (B12) ref )
Since there is a conversion relationship between the apparent reflectance (TOA) and the radiance, the present example is simplified to direct fitting regression of TOA data.
The data of methane at sentinel number two are represented by B11 (center wavelength of 1610 nm) and B12 (center wavelength of 2190 nm). The other 11 wave bands are respectively fitted to B11 and B12, the fitting result of the B12 is tested to be superior to the B11 effect, and the fitting result of the B12 is selected in the example to obtain the apparent reflectivity (B12) of the reference methane sensitive wave band ref ) The equation fitted to the apparent reflectance of the non-methane sensitive band is as follows:
B12 ref =0.3478*(B2)-0.5747*(B3)+0.3478*(B4)+1.2105*(B5)-0.9974*(B6)-0.2449*(B7)+0.9947*(B8A)+0.0180
step three, calculating and obtaining the ratio of the reference methane sensitive waveband and the methane sensitive waveband of the plateau area
B12ref is calculated by utilizing the fitting formula according to the sentinel data of the whole wave band of the fairy ground, and the sentinel data is divided by the B12 wave band to obtain the ratio of the reference methane sensitive wave band and the methane sensitive wave band of the fairy ground.
Identifying a high-value area of the ratio of the pool area to obtain a methane emission source
And (3) removing the water body by using the water body index NDVI or the published water body distribution data, and removing shadows from the visible light wave bands B2, B3 and B4 which are lower than the former 5 percent.
The high-value area can be easily obtained by imaging the contrast value data, and the distribution of the methane emission point source is obtained.
As shown in fig. 2, the image data map of 17 days at 1 month and 1 year in the plateau region 2021 according to the embodiment of the present invention and the random sampling points;
FIG. 3 shows a plateau B12 in an embodiment of the present invention ref Schematic distribution of/B12;
as shown in figures 4-6 of the drawings,three B12 illustrating the Rich platform region ref The methane emission source corresponding to the high-value area of the/B12, wherein, the B12 of the typical refuse landfill I of the plateau area is shown in figure 4 ref a/B12 area, a corresponding sentinel second image and a high-definition map; FIG. 5 is B12 of a typical landfill II of a Fennel district ref a/B12 area, a corresponding sentinel second image and a high-definition map; FIG. 6 shows B12 of a typical sewage treatment plant in a plateau area ref a/B12 area, a corresponding sentinel second image and a high-definition map.
Compared with the prior art, the embodiment of the invention utilizes the high-resolution satellite data, obtains the relationship between the methane sensitive waveband data and the non-methane sensitive waveband data by fitting, screens out the data which do not meet the relationship in the high-resolution satellite data, and finds out the corresponding position, namely the area where the methane emission source is located. The high space-time resolution of the methane emission source of the urban complex underlying surface is identified by setting the reference methane sensitive waveband as the characteristic waveband, and a reference is provided for supervision of the methane emission source.
The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are merely exemplary embodiments of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.
Claims (10)
1. A method for identifying a methane emission source of an urban complex underlying surface is characterized by comprising the following steps:
acquiring high-resolution satellite data of a target area;
fitting the high-resolution satellite data to obtain reference methane sensitive waveband data;
acquiring the ratio of the reference methane sensitive waveband data to the methane sensitive waveband data;
and acquiring a methane emission source of the target area according to the ratio.
2. The identification method of claim 1, wherein prior to said fitting said high resolution satellite data to obtain reference methane sensitive band data, said method further comprises:
and based on a methane emission source database, removing the high-resolution satellite data of the area influenced by the methane emission source to obtain the removed high-resolution satellite data.
3. The identification method of claim 1 or 2, wherein prior to said fitting said high resolution satellite data to obtain reference methane sensitive band data, said method further comprises:
and randomly screening the high-resolution satellite data to obtain the screened high-resolution satellite data.
4. The identification method of claim 3, wherein said fitting said high resolution satellite data to obtain reference methane sensitive band data comprises:
performing regression fitting on the methane sensitive waveband data through the non-methane sensitive waveband data to determine a fitting function;
and obtaining reference methane sensitive waveband data according to the fitting function and the non-methane sensitive waveband data.
5. The identification method according to claim 4, wherein the obtaining the ratio of the reference methane sensitive band data to the methane sensitive band data comprises:
and carrying out ratio operation on the reference methane sensitive waveband data based on the same geographical position information and the corresponding methane sensitive waveband data to obtain the ratio of the reference methane sensitive waveband data to the methane sensitive waveband data.
6. The identification method according to claim 1, wherein the obtaining the ratio of the reference methane sensitive band data to the methane sensitive band data comprises:
when the target area is positioned in the same scene image, determining the ratio according to the reference methane sensitive waveband representation reflectivity in the reference methane sensitive waveband data and the methane sensitive waveband representation reflectivity in the methane sensitive waveband data;
or
And determining the ratio according to the reference methane sensitive waveband radiance in the reference methane sensitive waveband data and the methane sensitive waveband radiance in the methane sensitive waveband data.
7. The identification method according to claim 1 or 2, wherein the obtaining of the methane emission source of the target area according to the ratio comprises:
if the ratio is less than or equal to 1, no methane is discharged from the corresponding position in the target area;
if the ratio is larger than 1, methane is discharged from the corresponding position in the target area;
screening to obtain a methane emission source in the target area with the ratio being greater than 1;
or
Processing the ratio by using a median filtering or deep learning algorithm to obtain a high-value area;
and screening to obtain a methane emission source according to the high-value area.
8. The method of claim 6, wherein the ratio is determined based on the radiation in the reference methane sensitive band data and the radiation in the methane sensitive band data according to the following formula:
r= F(B ref ) / F(B)
wherein r is the ratio of the radiance of a reference methane sensitive waveband to the radiance of the methane sensitive waveband; f (B) ref ) Is a reference methane sensitive band radiance; f (B) is the radiance of the methane sensitive waveband.
9. The identification method according to claim 1, wherein before acquiring the methane emission source of the target area according to the ratio, eliminating water bodies and shadow areas in the target area comprises:
and removing the water body with interference according to the analyzed water body index, and removing the shadow area by using a low-value area of a visible light wave band.
10. The identification method according to claim 7, further comprising:
displaying the high-value area on a high-definition map through a first corresponding relation to obtain a methane emission source of the target area; the first corresponding relation is the longitude and latitude relation between the longitude and latitude of the high-value area and the longitude and latitude of the corresponding high-definition map.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310040669.2A CN115830472B (en) | 2023-01-12 | 2023-01-12 | Urban complex underlying surface methane emission source identification method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310040669.2A CN115830472B (en) | 2023-01-12 | 2023-01-12 | Urban complex underlying surface methane emission source identification method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115830472A true CN115830472A (en) | 2023-03-21 |
CN115830472B CN115830472B (en) | 2023-05-16 |
Family
ID=85520592
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310040669.2A Active CN115830472B (en) | 2023-01-12 | 2023-01-12 | Urban complex underlying surface methane emission source identification method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115830472B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117593658A (en) * | 2023-11-21 | 2024-02-23 | 岭南师范学院 | BP neural network-based earth surface high-resolution methane inversion method and application method |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20180039885A1 (en) * | 2016-08-04 | 2018-02-08 | International Business Machines Corporation | Satellite-based location identification of methane-emitting sites |
CN113361309A (en) * | 2020-03-18 | 2021-09-07 | 山东交通学院 | Urban multi-haze pollution parameter monitoring method based on medium-high resolution remote sensing image |
CN113791042A (en) * | 2021-08-15 | 2021-12-14 | 中国科学技术大学 | Satellite remote sensing atmospheric formaldehyde total amount of slope under low signal-to-noise ratio |
CN113869143A (en) * | 2021-09-09 | 2021-12-31 | 武汉大学 | Method for quickly identifying and extracting methane point source based on hyperspectral imager |
CN114839327A (en) * | 2022-07-01 | 2022-08-02 | 北京英视睿达科技股份有限公司 | Emission source CO of satellite monitoring target area 2 Hourly discharge method and system |
CN115205709A (en) * | 2022-09-15 | 2022-10-18 | 眉山环天智慧科技有限公司 | Forest fire point identification method based on satellite remote sensing |
CN115235997A (en) * | 2022-07-04 | 2022-10-25 | 南通智能感知研究院 | Soil texture inversion method based on satellite hyperspectral image |
-
2023
- 2023-01-12 CN CN202310040669.2A patent/CN115830472B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20180039885A1 (en) * | 2016-08-04 | 2018-02-08 | International Business Machines Corporation | Satellite-based location identification of methane-emitting sites |
CN113361309A (en) * | 2020-03-18 | 2021-09-07 | 山东交通学院 | Urban multi-haze pollution parameter monitoring method based on medium-high resolution remote sensing image |
CN113791042A (en) * | 2021-08-15 | 2021-12-14 | 中国科学技术大学 | Satellite remote sensing atmospheric formaldehyde total amount of slope under low signal-to-noise ratio |
CN113869143A (en) * | 2021-09-09 | 2021-12-31 | 武汉大学 | Method for quickly identifying and extracting methane point source based on hyperspectral imager |
CN114839327A (en) * | 2022-07-01 | 2022-08-02 | 北京英视睿达科技股份有限公司 | Emission source CO of satellite monitoring target area 2 Hourly discharge method and system |
CN115235997A (en) * | 2022-07-04 | 2022-10-25 | 南通智能感知研究院 | Soil texture inversion method based on satellite hyperspectral image |
CN115205709A (en) * | 2022-09-15 | 2022-10-18 | 眉山环天智慧科技有限公司 | Forest fire point identification method based on satellite remote sensing |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117593658A (en) * | 2023-11-21 | 2024-02-23 | 岭南师范学院 | BP neural network-based earth surface high-resolution methane inversion method and application method |
Also Published As
Publication number | Publication date |
---|---|
CN115830472B (en) | 2023-05-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Liu et al. | Improve ground-level PM2. 5 concentration mapping using a random forests-based geostatistical approach | |
Brunner et al. | Accounting for the vertical distribution of emissions in atmospheric CO< sub> 2</sub> simulations | |
Kumar et al. | Evaluation and utilization of MODIS and CALIPSO aerosol retrievals over a complex terrain in Himalaya | |
Naboureh et al. | An integrated object-based image analysis and CA-Markov model approach for modeling land use/land cover trends in the Sarab plain | |
Heenkenda et al. | Quantifying mangrove chlorophyll from high spatial resolution imagery | |
Chen et al. | Methane emissions from the Munich Oktoberfest | |
CN111415309B (en) | High-resolution remote sensing image atmosphere correction method based on minimum reflectivity method | |
Khan et al. | Monitoring land use land cover changes and its impacts on land surface temperature over Mardan and Charsadda Districts, Khyber Pakhtunkhwa (KP), Pakistan | |
Huang et al. | Estimating 2013–2019 NO2 exposure with high spatiotemporal resolution in China using an ensemble model | |
Hou et al. | A study on water quality parameters estimation for urban rivers based on ground hyperspectral remote sensing technology | |
CN116148188A (en) | Air-space-ground integrated lake water quality tracing method, system, equipment and storage medium | |
Bonannella et al. | Forest tree species distribution for Europe 2000–2020: mapping potential and realized distributions using spatiotemporal machine learning | |
Nickel et al. | Reorganisation of a long-term monitoring network using moss as biomonitor for atmospheric deposition in Germany | |
CN115830472B (en) | Urban complex underlying surface methane emission source identification method | |
Luther et al. | Observational constraints on methane emissions from Polish coal mines using a ground-based remote sensing network | |
CN114814167B (en) | Soil heavy metal content inversion method fusing multi-source environment variables and spectral information | |
CN114397277A (en) | Unmanned aerial vehicle water chlorophyll remote sensing detection system | |
Al-Shalan et al. | Methane emissions in Kuwait: Plume identification, isotopic characterisation and inventory verification | |
Ya'acob et al. | Temporal forest change detection and forest health assessment using Remote Sensing | |
Papachristopoulou et al. | Aerosol optical depth regime over megacities of the world | |
Guo et al. | Comparison of atmospheric CO2, CH4, and CO at two stations in the Tibetan Plateau of China | |
Ghurah et al. | Temporal change detection of land use/land cover using GIS and remote sensing techniques in South Ghor Regions, Al-Karak, Jordan | |
WO2023229705A1 (en) | System and method for automatically estimating gas emission parameters | |
Wang et al. | Ammonia nitrogen monitoring of urban rivers with UAV-borne hyperspectral remote sensing imagery | |
Mauya et al. | Modelling and Mapping Above Ground Biomass Using Sentinel 2 and Planet Scope Remotely Sensed Data in West Usambara Tropical Rainforests, Tanzania |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |