CN113887780B - Satellite remote sensing earth surface temperature estimation method, device and equipment - Google Patents
Satellite remote sensing earth surface temperature estimation method, device and equipment Download PDFInfo
- Publication number
- CN113887780B CN113887780B CN202110991141.4A CN202110991141A CN113887780B CN 113887780 B CN113887780 B CN 113887780B CN 202110991141 A CN202110991141 A CN 202110991141A CN 113887780 B CN113887780 B CN 113887780B
- Authority
- CN
- China
- Prior art keywords
- land
- estimated
- temperature
- land cover
- temperature difference
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000005259 measurement Methods 0.000 claims abstract description 34
- 101100533803 Homo sapiens SLCO1B3 gene Proteins 0.000 claims description 87
- 102100027239 Solute carrier organic anion transporter family member 1B3 Human genes 0.000 claims description 87
- 101150026176 zfyve28 gene Proteins 0.000 claims description 87
- 101001065658 Homo sapiens Leukocyte-specific transcript 1 protein Proteins 0.000 claims description 27
- 102100027233 Solute carrier organic anion transporter family member 1B1 Human genes 0.000 claims description 27
- 238000004891 communication Methods 0.000 claims description 17
- 101000836287 Homo sapiens Putative solute carrier organic anion transporter family member 1B7 Proteins 0.000 claims description 12
- 102100027230 Putative solute carrier organic anion transporter family member 1B7 Human genes 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 6
- 238000012544 monitoring process Methods 0.000 abstract description 12
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 13
- 230000008859 change Effects 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000006978 adaptation Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000012417 linear regression Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000005215 recombination Methods 0.000 description 2
- 230000006798 recombination Effects 0.000 description 2
- 244000025254 Cannabis sativa Species 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000010792 warming Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J5/007—Radiation pyrometry, e.g. infrared or optical thermometry for earth observation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J2005/0077—Imaging
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A30/00—Adapting or protecting infrastructure or their operation
- Y02A30/60—Planning or developing urban green infrastructure
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Strategic Management (AREA)
- Environmental & Geological Engineering (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Game Theory and Decision Science (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Geology (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Radar Systems Or Details Thereof (AREA)
- Radiation Pyrometers (AREA)
Abstract
The embodiment of the invention provides a satellite remote sensing earth surface temperature estimation method, device and equipment, wherein the method comprises the following steps: obtaining measured surface temperatures of various land coverage types; obtaining a temperature difference relation between a first land cover type and the rest land cover types in at least one land cover type according to the measured surface temperatures of the plurality of land cover types; obtaining a first estimated temperature difference between the mixed pixel temperature and the estimated surface temperature of the first earth coverage type by satellite measurement; obtaining a second estimated temperature difference between the estimated surface temperature of the other land coverage types and the estimated surface temperature of the first land coverage type according to the temperature difference relation and the first estimated temperature difference; obtaining target surface temperatures of the plurality of land coverage types according to the second estimated temperature difference; the earth surface temperature information with high observation frequency and high spatial resolution can be obtained, and the timeliness and the precision of satellite remote sensing urban heat island monitoring are improved.
Description
Technical Field
The invention relates to the technical field of earth surface temperature monitoring, in particular to a satellite remote sensing earth surface temperature estimation method, device and equipment.
Background
With the acceleration of urban development and climate warming, the aggravation of urban heat island effect becomes an important factor of concern in urban environment monitoring, urban planning and the like. The meteorological satellite has the characteristics of wide observation range, high imaging frequency, sensitivity to ground surface temperature information and the like, and can play an important role in urban heat island monitoring. The land resource, high-grade satellite series has high spatial resolution and can finely invert the capability of the urban underlying land coverage information.
Disclosure of Invention
The invention provides a satellite remote sensing earth surface temperature estimation method, device and equipment. The method can obtain the refined underlying surface temperature information with high observation frequency and high spatial resolution, provides support for urban heat island monitoring and evaluation, and improves the timeliness and precision of satellite remote sensing urban heat island monitoring.
In order to solve the above technical problems, the embodiments of the present invention provide the following solutions:
a satellite remote sensing earth surface temperature estimation method comprises the following steps:
obtaining measured surface temperatures of various land coverage types;
obtaining a temperature difference relation between a first land cover type and the rest land cover types in at least one land cover type according to the measured surface temperatures of the plurality of land cover types;
Obtaining a first estimated temperature difference between the mixed pixel temperature and the estimated surface temperature of the first earth coverage type by satellite measurement;
obtaining a second estimated temperature difference between the estimated surface temperature of the other land coverage types and the estimated surface temperature of the first land coverage type according to the temperature difference relation and the first estimated temperature difference;
and obtaining the target surface temperatures of the plurality of land coverage types according to the second estimated temperature difference.
Optionally, obtaining a temperature difference relationship between a first land cover type and the remaining land cover types in the at least one land cover type according to the surface temperatures of the plurality of land cover types includes:
obtaining a first measured temperature difference between the measured temperature of the first land cover type and the measured temperature of the second land cover type according to the measured temperature of the first land cover type and the measured temperature of the second land cover type;
obtaining a second measured temperature difference between the measured temperature of the first land cover type and the measured temperature of the remaining land cover type according to the measured temperature of the first land cover type and the measured temperature of the remaining land cover type;
And obtaining a temperature difference relation between a first land cover type and the rest land cover types in at least one land cover type according to the first measured temperature difference and the second measured temperature difference.
Optionally, the temperature difference relationship is: deltaLST 1 The rest are ＝a×△LST1 Second one +b；
Wherein DeltaLST 1 The rest are For the second measured temperature difference, deltaLST 1 Second one For the first measured temperature difference, a and b are relationship constants.
Optionally, obtaining a first estimated temperature difference between the temperature of the mixed pixels and the estimated surface temperature of the first earth coverage type by satellite measurement comprises:
obtaining an estimated surface temperature and a mixed pixel temperature of a first land coverage type through satellite measurement;
according to the estimated surface temperature and the mixed pixel temperature of the first land coverage type, the following formula is adopted:
△LST2 first one ＝LST2 Mixed pixel -LST2 First land cover type Obtaining a first estimated temperature difference;
wherein DeltaLST 2 First one For a first estimated temperature difference, LST2 Mixed pixel To mix pixel temperature, LST2 First land cover type An estimated surface temperature for a first earth coverage type.
Optionally, the mixed pixel temperature is the pixel with the highest proportion of the area of the second land coverage type in the area of one pixel.
Optionally, obtaining a second estimated temperature difference between the estimated surface temperature of the remaining land cover types and the estimated surface temperature of the first land cover type according to the temperature difference relationship and the first estimated temperature difference, including:
based on the temperature difference relation and the first estimated temperature difference, the method passes through the formula delta LST2 Second one ＝a×△LST2 First one +b obtaining a second estimated temperature difference;
wherein DeltaLST 2 Second one For the second estimated temperature difference, ΔLST2 First one For the first estimated temperature difference, a and b are relationship constants.
Optionally, obtaining the estimated surface temperature of the remaining land coverage type according to the second estimated temperature difference includes:
according to the second estimated temperature difference, the following formula is adopted: LST2 Type of remaining land cover ＝LST2 First land cover type +△LST2 Second one Obtaining an estimated surface temperature for the remaining land coverage types;
wherein LST2 Type of remaining land cover LST2 for estimating surface temperature for remaining land cover types First land cover type For an estimated surface temperature of a first land cover type ΔLST2 Second one For a second estimated temperature difference;
according to the formula:
obtaining a target surface temperature of the second land cover type;
wherein LST3 Second land cover type LST2 for target surface temperature of second land cover type Mixed pixel To mix pixel temperature, LST2 Type of remaining land cover LST2 for estimating surface temperature for remaining land cover types First land cover type For the estimated surface temperature of the first earth coverage type, R The rest are R is the ratio of the area of the residual land coverage type to the area of the mixed pixel Second one R is the ratio of the area of the second land cover type to the area of the mixed pixel First one A ratio of an area of the first land cover type to an area of the hybrid pixel;
and obtaining the target surface temperature of the other land coverage types according to the target surface temperature of the second land coverage type and the estimated surface temperature of the first land coverage type.
The invention also provides a satellite remote sensing earth surface temperature estimation device, which comprises:
the acquisition module is used for acquiring actual earth surface temperatures of various land coverage types;
the processing module is used for obtaining the temperature difference relation between the first land coverage type and the rest land coverage types in at least one land coverage type according to the actually measured surface temperatures of the plurality of land coverage types; obtaining a first estimated temperature difference between the mixed pixel temperature and the estimated surface temperature of the first earth coverage type by satellite measurement; obtaining a second estimated temperature difference between the estimated surface temperature of the other land coverage types and the estimated surface temperature of the first land coverage type according to the temperature difference relation and the first estimated temperature difference; and obtaining the target surface temperatures of the plurality of land coverage types according to the second estimated temperature difference.
The present invention provides a computing device comprising: the device comprises a processor, a memory, a communication interface and a communication bus, wherein the processor, the memory and the communication interface complete communication with each other through the communication bus;
the memory is used for storing at least one executable instruction, and the executable instruction enables the processor to execute the operation corresponding to the satellite remote sensing earth surface temperature estimation method.
The present invention provides a computer readable storage medium storing instructions that, when executed on a computer, cause the computer to perform a method of estimating satellite remote sensing earth surface temperature as described above.
The scheme of the invention at least comprises the following beneficial effects:
according to the scheme, the actual measured surface temperature of various land coverage types is obtained; obtaining a temperature difference relation between a first land cover type and the rest land cover types in at least one land cover type according to the measured surface temperatures of the plurality of land cover types; obtaining a first estimated temperature difference between the mixed pixel temperature and the estimated surface temperature of the first earth coverage type by satellite measurement; obtaining a second estimated temperature difference between the estimated surface temperature of the other land coverage types and the estimated surface temperature of the first land coverage type according to the temperature difference relation and the first estimated temperature difference; and obtaining the target surface temperatures of the plurality of land coverage types according to the second estimated temperature difference. The satellite remote sensing refined underlying surface temperature with high observation frequency and high spatial resolution can be obtained, and the timeliness and the precision of satellite remote sensing urban heat island monitoring are improved.
Drawings
FIG. 1 is a flow chart of a method for estimating satellite remote sensing earth surface temperature according to an embodiment of the invention;
fig. 2 is a schematic structural diagram of an estimation device for satellite remote sensing earth surface temperature according to an embodiment of the present invention.
Detailed Description
Exemplary embodiments of the present invention will be described in more detail below with reference to the accompanying drawings. While exemplary embodiments of the present invention are shown in the drawings, it should be understood that the present invention may be embodied in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the invention to those skilled in the art.
As shown in fig. 1, the invention provides a satellite remote sensing earth surface temperature estimation method, which comprises the following steps:
step 11, obtaining actual measured surface temperatures of various land coverage types;
step 12, obtaining a temperature difference relation between a first land cover type and the rest land cover types in at least one land cover type according to the measured surface temperatures of the land cover types;
step 13, obtaining a first estimated temperature difference between the temperature of the mixed pixels and the estimated surface temperature of the first land coverage type through satellite measurement;
Step 14, obtaining a second estimated temperature difference between the estimated surface temperature of the other land coverage types and the estimated surface temperature of the first land coverage type according to the temperature difference relation and the first estimated temperature difference;
and step 15, obtaining the target surface temperature of the other land coverage types according to the second estimated temperature difference.
In the embodiment, ground surface temperatures of various land coverage types are obtained based on ground actual measurement, a temperature difference relation between the first land coverage type and the rest land coverage type is obtained according to the ground surface temperatures of the first land coverage type, the second land coverage type and the rest land coverage type, when a satellite monitors the ground surface temperature of a region in a pixel, a first estimated temperature difference is obtained through a satellite measurement mode, a second estimated temperature difference can be obtained according to the temperature difference relation and the first estimated temperature difference, and a target ground surface temperature of the rest land coverage type is obtained according to the second estimated temperature difference; the satellite remote sensing refined underlying surface temperature information with high observation frequency and high spatial resolution can be obtained, and the timeliness and the precision of satellite remote sensing urban heat island monitoring are improved.
In an alternative embodiment of the present invention, step 12 includes:
Step 121, obtaining a first measured temperature difference between the measured temperature of the first land cover type and the measured temperature of the second land cover type according to the measured temperature of the first land cover type and the measured temperature of the second land cover type;
step 122, obtaining a second measured temperature difference between the measured temperature of the first land cover type and the measured temperature of the other land cover type according to the measured temperature of the first land cover type and the measured temperature of the other land cover type;
and step 123, obtaining a temperature difference relation between the first land cover type and the rest land cover types in at least one land cover type according to the first measured temperature difference and the second measured temperature difference.
In this embodiment, at least one of a plurality of land cover types is included in one pixel, the measured surface temperature of the plurality of land cover types can be measured on the ground preferably periodically by an infrared thermometer, a first measured temperature difference and a second measured temperature difference are established between the first land cover type and the second land cover type and between the first land cover type and the rest of land cover types respectively, and a regression formula of the second measured temperature difference along with the change of the first measured temperature difference is established according to the first measured temperature difference and the second measured temperature difference, wherein the regression formula is the temperature difference relation between the first land cover type and the rest of land cover types;
It should be noted that the plurality of land cover types preferably include a first land cover type, a second land cover type, and other land cover types, wherein the second land cover type preferably includes an urban area, the first land cover type may select a land cover type that generally appears in a large area around the urban area, the first land cover type is preferably a farmland, and the remaining land cover types preferably include at least one of a forest land, a grass land, a water body, and a bare land other than the first land cover type and the second land cover type, in order to facilitate an algorithm process for satellite remote sensing of a surface temperature.
In an alternative embodiment of the present application, the temperature difference relationship is:
△LST1 the rest are ＝a×△LST1 Second one +b；
Wherein DeltaLST 1 The rest are For the second measured temperature difference, deltaLST 1 Second one For the first measured temperature difference, a and b are relationship constants.
In the embodiment, the regression formula corresponding to the temperature difference relation of the second measured temperature difference along with the change of the first measured temperature difference is a first-order linear regression formula, but the application is not limited thereto, and the regression formula can be obtained as the temperature difference relation according to the actual situation of the accuracy requirement of the second measured temperature difference along with the change of the first measured temperature difference.
In an alternative embodiment of the present invention, step 13 includes:
step 131, obtaining an estimated surface temperature and a mixed pixel temperature of the first land coverage type through satellite measurement;
step 132, according to the estimated surface temperature and the mixed pixel temperature of the first land coverage type, the following formula is adopted: deltaLST 2 First one ＝LST2 Mixed pixel -LST2 First land cover type Obtaining a first estimated temperature difference;
wherein DeltaLST 2 First one For a first estimated temperature difference, LST2 Mixed pixel To mix pixel temperature, LST2 First land cover type An estimated surface temperature for a first earth coverage type.
In addition, the mixed pixel is a pixel with the highest proportion of the area of the second land coverage type in the area of one pixel.
In this embodiment, an estimated surface temperature of a first land coverage type and a mixed pixel temperature are obtained through satellite measurement, where the estimated surface temperature of the first land coverage type is measured, and all the land coverage types in the measured pixels are surface temperatures corresponding to the pixels of the first land coverage type, and the mixed pixel is a pixel with the highest area ratio of the second land coverage type in the pixels of the second land coverage type, and it is noted that the mixed pixel includes the second land coverage type and may also include the remaining land coverage types;
In a specific embodiment 1, when all pixels in a unit pixel set are SP1, SP2, SP3, and SP4, respectively, wherein SP1 includes 100% of the first land cover type, SP2 includes 90% of the first land cover type and 10% of the second land cover type, SP3 includes 80% of the second land cover type, 10% of the first land cover type, and 10% of the remaining land cover type, and SP4 includes 80% of the remaining land cover type, 10% of the first land cover type, and 10% of the second land cover type; therefore, when the satellite remotely measures the estimated surface temperature of the first land coverage type, the sub-pixel SP1 is selected to measure the estimated surface temperature of the first land coverage type, and when the satellite measures the mixed pixel temperature, the SP3 is selected as the mixed pixel to measure the mixed pixel temperature;
optionally, satellite measuring the estimated surface temperature of the at least one earth coverage type may include:
acquiring first radiance of at least one land cover type through a sensor of a satellite;
according to the radiance of at least one land coverage type, respectively calculating the surface temperature corresponding to the first radiance of the at least one land coverage type by utilizing a Planck formula;
Processing the surface temperature corresponding to the first radiance to obtain estimated surface temperature of at least one land coverage type;
in this embodiment, the processing the surface temperature corresponding to the first radiance specifically includes:
by the formulaObtaining inversion surface temperature by using a far infrared channel and an infrared splitting window channel on a satellite; the inversion surface temperature is the estimated surface temperature of at least one land coverage type; wherein T is s To invert the surface temperature, T 4 Is the radiance of a far infrared channel, T 5 Is the radiance of an infrared splitting window channel, A 0 P is the average emissivity of the far infrared channel, and M is the emissivity difference of the infrared splitting window channel;
specifically, P and M are:
wherein,△ε＝ε 4 -ε 5 ，ε 4 emissivity of far infrared channel epsilon 5 For the emissivity of the infrared splitting window channel, ε is the average emissivity of the far infrared channel and the infrared splitting window channel, Δε is the difference in emissivity between the far infrared channel and the infrared splitting window channel, A 0 The values of α, β, γ ', α ' and β ' are related to the relevant parameters of the satellite used for the measurement, and may depend on the satellite used in the actual measurement.
In an alternative embodiment of the present invention, step 14 includes:
Step 141, according to the temperature difference relation and the first estimated temperature difference, passing through the formula DeltaLST 2 Second one ＝a×△LST2 First one +b obtaining a second estimated temperature difference;
wherein DeltaLST 2 Second one For the second estimated temperature difference, ΔLST2 First one For the first estimated temperature difference, a and b are relationship constants.
In this embodiment, ΔLST2 is calculated by the formula Second one ＝a×△LST2 First one And +b, obtaining a second estimated temperature difference, wherein the second estimated temperature difference is calculated based on the first estimated temperature difference and a regression formula corresponding to the temperature difference relation obtained according to actual measurement, and if the estimated surface temperature of the corresponding land coverage type is obtained by directly using the second estimated temperature difference through calculation, a great amount of deviation between the surface temperature estimated result of satellite remote sensing monitoring and the actual measurement result of the ground can occur, so that the second estimated temperature difference needs to be processed in step 15.
In an alternative embodiment of the present invention, step 15 includes:
step 151, according to the second estimated temperature difference, passing through the formula: LST2 Type of remaining land cover ＝LST2 First land cover type +△LST2 Second one Obtaining an estimated surface temperature for the remaining land coverage types;
wherein LST2 Type of remaining land cover LST2 for estimating surface temperature for remaining land cover types First land cover type For an estimated surface temperature of a first land cover type ΔLST2 Second one For a second estimated temperature difference;
step 152, according to the formula:
obtaining a target surface temperature of the second land cover type;
wherein LST3 Second land cover type LST2 for target surface temperature of second land cover type Mixed pixel To mix pixel temperature, LST2 Type of remaining land cover LST2 for estimating surface temperature for remaining land cover types First land cover type For the estimated surface temperature of the first earth coverage type, R The rest are R is the ratio of the area of the residual land coverage type to the area of the mixed pixel Second one R is the ratio of the area of the second land cover type to the area of the mixed pixel First one A ratio of an area of the first land cover type to an area of the hybrid pixel;
and step 153, obtaining the target surface temperature of the other land coverage types according to the target surface temperature of the second land coverage type and the estimated surface temperature of the first land coverage type.
In this embodiment, according to the second estimated temperature difference, the estimated surface temperature of the remaining land coverage type is calculated, and then according to the formula in step 152, the target surface temperature of the second land coverage type is obtained, and further the target surface temperature of the remaining land coverage type is obtained;
Wherein step 152 may comprise: acquiring the ratio R of the area of the first land coverage type in the mixed pixel to the area of the mixed pixel First one Ratio R of area of second land cover type in mixed pixel area Second one Ratio R of area of other land cover type in mixed pixel area The rest are The method comprises the steps of carrying out a first treatment on the surface of the Wherein R is The rest are +R First one +R Second one =1; the radiance of the hybrid pixel is described as follows:
L mixed pixel ＝R First one ×L First one +R Second one ×L Second one +R The rest are ×L The rest are ；
Wherein L is Mixed pixel For the radiance of the mixed pixels, L First one For the first land cover type of radiance, L Second one Is the radiance of the second land cover type, L The rest are For the radiance of the remaining land cover type, R The rest are R is the ratio of the area of the residual land coverage type to the area of the mixed pixel Second one R is the ratio of the area of the second land cover type to the area of the mixed pixel First one Is the firstThe ratio of the area of the land cover type in the area of the mixed pixel;
LST is calculated by boltzmann formula l=σ 4 Wherein σ=5.6704×10 8 The radiance of the hybrid pixel is described as follows:
(LST2 mixed pixel ) 4 ＝R First one ×(LST3 First type of land ) 4 +R Second one ×(LST2 Second type of land ) 4 +R The rest are ×(LST2 Other types of land ) 4 And then the formula is obtained:
the target surface temperature of the second land coverage type can be calculated by the formula;
wherein LST3 Second land cover type LST2 for target surface temperature of second land cover type Mixed pixel To mix pixel temperature, LST2 Type of remaining land cover LST2 for estimating surface temperature for remaining land cover types First land cover type For the estimated surface temperature of the first earth coverage type, R The rest are R is the ratio of the area of the residual land coverage type to the area of the mixed pixel Second one R is the ratio of the area of the second land cover type to the area of the mixed pixel First one A ratio of an area of the first land cover type to an area of the hybrid pixel;
when the first land cover type is used for satellite measurement to select the pixels, the pixels with the first land cover type are selected, so that the estimated surface temperature of the first land cover type can be regarded as the target surface temperature of the first land cover type, and the target surface temperatures of the other land cover types are obtained according to the estimated surface temperatures of the second land cover type and the first land cover type.
In a specific embodiment 2, the high-resolution satellite remotely measures the ground temperature of urban areas with a range of land coverage types, including various land coverage types, and the high-resolution satellite is used to more accurately measure the specific land coverage types in the pixels, including woodland, grassland, farmland, water, bare land (sparse vegetation) and urban areas (roads, squares, buildings, etc.);
step 1: obtaining the actual measurement temperature of each land coverage type in the area range, wherein the actual measurement temperature of farmland is LST Farm land The actual measurement temperature of urban areas is LST Urban area The actual measurement temperature of the forest land is LST Woodlands The measured temperature of the grassland is LST Grassland The measured temperature of the water body is LST Water body The measured temperature of bare land is LST Bare land ；
Step 2: according to the measured temperature of each land cover type, the temperature difference between the measured temperature of the land cover type except the farmland and the measured temperature of the farmland is obtained respectively, specifically:
actual measurement temperature difference delta LST1 of urban area Urban area ＝LST1 Urban area -LST1 Farm land ；
Actual measurement of temperature difference delta LST1 in forest land Woodlands ＝LST1 Woodlands -LST1 Farm land ；
Grassland actual measurement temperature difference delta LST1 Grassland ＝LST1 Grassland -LST1 Farm land ；
Actual measurement of temperature difference delta LST1 of water body Water body ＝LST1 Water body -LST1 Farm land ；
Actually measured temperature difference delta LST1 of bare land Bare land ＝LST1 Bare land -LST1 Farm land ；
Step 3: obtaining DeltaLST 1 according to the temperature difference between the measured temperature of the land cover type except the farmland and the measured temperature of the farmland Urban area Respectively with delta LST1 Woodlands 、△LST1 Grassland 、△LST1 Water body DeltaLST 1 Bare land In this embodiment, the temperature difference relationship is a first-order linear regression formula, and specifically:
△LST1 grassland ＝a1*△LST1 Urban area +b1；
△LST1 Woodlands ＝a2*△LST1 Urban area +b2；
△LST1 Bare land ＝a3*△LST1 Urban area +b3；
△LST1 Water body ＝a4*△LST1 Urban area +b4；
Wherein a1, a2, a3, a4, b1, b2, b3 and b4 are relationship constants;
step 4, measuring the earth surface temperature of each pixel in the area range through satellite remote sensing, firstly measuring the earth surface temperature of the pixel with the full land coverage type as a farmland, selecting a mixed pixel with the largest urban area occupation area in the pixel, taking the earth surface temperature of the mixed pixel as the estimated earth surface temperature of the satellite remote sensing, and obtaining the earth surface temperature of the pixel with the full land coverage type as a farmland through the formula delta LST2 Urban area ＝LST2 Urban area -LST2 Farm land Calculating to obtain the temperature difference between the estimated surface temperature of the urban area and the surface temperature of the farmland;
step 5, since the mixed pixel includes a plurality of land cover types, the boltzmann formula l=σ×lst is combined 4 Wherein σ=5.6704×10 8 Can be calculated by the formula
LST3 Urban area 4 ×R Urban area ＝(LST2 Urban area ) 4 -R First one ×(LST2 Farm land ) 4 -R Woodlands ×(LST2 Woodlands ) 4
-R Grassland ×(LST2 Grassland ) 4 -R Water body ×(LST2 Water body ) 4
-R Bare land ×(LST2 Bare land ) 4
The earth surface temperature of the urban area measured by the corrected satellite remote sensing is obtained through calculation, and according to the earth surface temperature of the urban area and the temperature difference relation in the step 3, the temperature difference between the estimated earth surface temperature of the urban area and the earth surface temperature of the farmland is combined to obtain the temperature difference between the farmland and the forest land, the bare land, the water body and the grassland respectively:
△LST2 grassland ＝a1*△LST3 Urban area +b1；
△LST2 Woodlands ＝a2*△LST3 Urban area +b2；
△LST2 Bare land ＝a3*△LST3 Urban area +b3；
△LST2 Water body ＝a4*△LST3 Urban area +b4；
Through the steps 1 to 5, the satellite remote sensing technology with high resolution is utilized, so that the surface temperature of the refined underlying surface can be observed, and particularly the method can be used for improving the urban heat island monitoring efficiency, so that the urban heat island monitoring result is more accurate.
In the embodiment of the invention, the actual measured surface temperature of various land coverage types is obtained; obtaining a temperature difference relation between a first land cover type and the rest land cover types in at least one land cover type according to the measured surface temperatures of the plurality of land cover types; obtaining a first estimated temperature difference between the mixed pixel temperature and the estimated surface temperature of the first earth coverage type by satellite measurement; obtaining a second estimated temperature difference between the estimated surface temperature of the other land coverage types and the estimated surface temperature of the first land coverage type according to the temperature difference relation and the first estimated temperature difference; and obtaining the target surface temperatures of the plurality of land coverage types according to the second estimated temperature difference. The satellite remote sensing refined underlying surface temperature information with high observation frequency and high spatial resolution can be obtained, and the timeliness and the precision of satellite remote sensing earth ground heat island monitoring are improved.
As shown in fig. 2, the present invention further provides an apparatus 20 for estimating satellite remote sensing earth surface temperature, the apparatus comprising:
an acquisition module 21 for acquiring measured surface temperatures of various land cover types;
a processing module 22, configured to obtain a temperature difference relationship between a first land cover type and the remaining land cover types in the at least one land cover type according to the measured surface temperatures of the plurality of land cover types; obtaining a first estimated temperature difference between the mixed pixel temperature and the estimated surface temperature of the first earth coverage type by satellite measurement; obtaining a second estimated temperature difference between the estimated surface temperature of the other land coverage types and the estimated surface temperature of the first land coverage type according to the temperature difference relation and the first estimated temperature difference; and obtaining the target surface temperatures of the plurality of land coverage types according to the second estimated temperature difference.
Optionally, obtaining a temperature difference relationship between a first land cover type and the remaining land cover types in the at least one land cover type according to the surface temperatures of the plurality of land cover types includes:
obtaining a first measured temperature difference between the measured temperature of the first land cover type and the measured temperature of the second land cover type according to the measured temperature of the first land cover type and the measured temperature of the second land cover type;
Obtaining a second measured temperature difference between the measured temperature of the first land cover type and the measured temperature of the remaining land cover type according to the measured temperature of the first land cover type and the measured temperature of the remaining land cover type;
and obtaining a temperature difference relation between a first land cover type and the rest land cover types in at least one land cover type according to the first measured temperature difference and the second measured temperature difference.
Optionally, the temperature difference relationship is: deltaLST 1 The rest are ＝a×△LST1 Second one +b；
Wherein DeltaLST 1 The rest are For the second measured temperature difference, deltaLST 1 Second one For the first measured temperature difference, a and b are relationship constants.
Optionally, obtaining a first estimated temperature difference between the temperature of the mixed pixels and the estimated surface temperature of the first earth coverage type by satellite measurement comprises:
obtaining an estimated surface temperature and a mixed pixel temperature of a first land coverage type through satellite measurement;
according to the estimated surface temperature and the mixed pixel temperature of the first land coverage type, the following formula is adopted:
△LST2 first one ＝LST2 Mixed pixel -LST2 First land cover type Obtaining a first estimated temperature difference;
wherein DeltaLST 2 First one For a first estimated temperature difference, LST2 Mixed pixel To mix pixel temperature, LST2 First land cover type An estimated surface temperature for a first earth coverage type.
Optionally, the mixed pixel temperature is the pixel with the highest proportion of the area of the second land coverage type in the area of one pixel.
Optionally, obtaining a second estimated temperature difference between the estimated surface temperature of the remaining land cover types and the estimated surface temperature of the first land cover type according to the temperature difference relationship and the first estimated temperature difference, including:
based on the temperature difference relation and the first estimated temperature difference, the method passes through the formula delta LST2 Second one ＝a×△LST2 First one +b obtaining a second estimated temperature difference;
wherein DeltaLST 2 Second one For the second estimated temperature difference, ΔLST2 First one For the first estimated temperature difference, a and b are relationship constants.
Optionally, obtaining the estimated surface temperature of the remaining land coverage type according to the second estimated temperature difference includes: according to the second estimated temperature difference, the following formula is adopted: LST2 Type of remaining land cover ＝LST2 First land cover type +△LST2 Second one Obtaining an estimated surface temperature for the remaining land coverage types;
wherein LST2 Type of remaining land cover LST2 for estimating surface temperature for remaining land cover types First land cover type For an estimated surface temperature of a first land cover type ΔLST2 Second one For a second estimated temperature difference;
according to the formula:
obtaining a target surface temperature of the second land cover type;
wherein LST3 Second land cover type LST2 for target surface temperature of second land cover type Mixed pixel To mix pixel temperature, LST2 The rest of the land is covered withCover type LST2 for estimating surface temperature for remaining land cover types First land cover type For the estimated surface temperature of the first earth coverage type, R The rest are R is the ratio of the area of the residual land coverage type to the area of the mixed pixel Second one R is the ratio of the area of the second land cover type to the area of the mixed pixel First one A ratio of an area of the first land cover type to an area of the hybrid pixel;
and obtaining the target surface temperature of the other land coverage types according to the target surface temperature of the second land coverage type and the estimated surface temperature of the first land coverage type.
It should be noted that, the apparatus 20 is an apparatus corresponding to the above method, and all implementation manners in the embodiment of the satellite remote sensing earth surface temperature estimation method are applicable to the embodiment of the apparatus, so that the same technical effects can be achieved.
Embodiments of the present invention provide a computing device comprising: the device comprises a processor, a memory, a communication interface and a communication bus, wherein the processor, the memory and the communication interface complete communication with each other through the communication bus;
the memory is used for storing at least one executable instruction, and the executable instruction enables the processor to execute the operation corresponding to the satellite remote sensing earth surface temperature estimation method.
It should be noted that, the computing device is a device corresponding to the above method, and all implementation manners in the embodiment of the above method for estimating the surface temperature are applicable to the embodiment of the computing device, so that the same technical effects can be achieved.
Embodiments of the present invention also provide a computer-readable storage medium storing instructions that, when executed on a computer, cause the computer to perform the method of estimating satellite remote sensing surface temperature as described above.
Those of ordinary skill in the art will appreciate that the various illustrative elements and algorithm steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware, or combinations of computer software and electronic hardware. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the solution. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present invention.
It will be clear to those skilled in the art that, for convenience and brevity of description, specific working procedures of the above-described systems, apparatuses and units may refer to corresponding procedures in the foregoing method embodiments, and are not repeated herein.
In the embodiments provided in the present invention, it should be understood that the disclosed apparatus and method may be implemented in other manners. For example, the apparatus embodiments described above are merely illustrative, e.g., the division of the units is merely a logical function division, and there may be additional divisions when actually implemented, e.g., multiple units or components may be combined or integrated into another system, or some features may be omitted or not performed. Alternatively, the coupling or direct coupling or communication connection shown or discussed with each other may be an indirect coupling or communication connection via some interfaces, devices or units, which may be in electrical, mechanical or other form.
The units described as separate units may or may not be physically separate, and units shown as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.
In addition, each functional unit in the embodiments of the present invention may be integrated in one processing unit, or each unit may exist alone physically, or two or more units may be integrated in one unit.
The functions, if implemented in the form of software functional units and sold or used as a stand-alone product, may be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention may be embodied essentially or in a part contributing to the prior art or in a part of the technical solution, in the form of a software product stored in a storage medium, comprising several instructions for causing a computer device (which may be a personal computer, a server, a network device, etc.) to perform all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a usb disk, a removable hard disk, a ROM, a RAM, a magnetic disk, or an optical disk, etc.
Furthermore, it should be noted that in the apparatus and method of the present invention, it is apparent that the components or steps may be disassembled and/or assembled. Such decomposition and/or recombination should be considered as equivalent aspects of the present invention. Also, the steps of performing the series of processes described above may naturally be performed in chronological order in the order of description, but are not necessarily performed in chronological order, and some steps may be performed in parallel or independently of each other. It will be appreciated by those of ordinary skill in the art that all or any of the steps or components of the methods and apparatus of the present invention may be implemented in hardware, firmware, software, or a combination thereof in any computing device (including processors, storage media, etc.) or network of computing devices, as would be apparent to one of ordinary skill in the art after reading this description of the invention.
The object of the invention can thus also be achieved by running a program or a set of programs on any computing device. The computing device may be a well-known general purpose device. The object of the invention can thus also be achieved by merely providing a program product containing program code for implementing said method or apparatus. That is, such a program product also constitutes the present invention, and a storage medium storing such a program product also constitutes the present invention. It is apparent that the storage medium may be any known storage medium or any storage medium developed in the future. It should also be noted that in the apparatus and method of the present invention, it is apparent that the components or steps may be disassembled and/or assembled. Such decomposition and/or recombination should be considered as equivalent aspects of the present invention. The steps of executing the series of processes may naturally be executed in chronological order in the order described, but are not necessarily executed in chronological order. Some steps may be performed in parallel or independently of each other.
While the foregoing is directed to the preferred embodiments of the present invention, it will be appreciated by those skilled in the art that various modifications and adaptations can be made without departing from the principles of the present invention, and such modifications and adaptations are intended to be comprehended within the scope of the present invention.
Claims (5)
1. The method for estimating the satellite remote sensing earth surface temperature is characterized by comprising the following steps of:
obtaining measured surface temperatures of various land coverage types;
obtaining a temperature difference relation between a first land cover type and the rest land cover types in at least one land cover type according to the measured surface temperatures of the plurality of land cover types;
obtaining a first estimated temperature difference between the mixed pixel temperature and the estimated surface temperature of the first earth coverage type by satellite measurement;
obtaining a second estimated temperature difference between the estimated surface temperature of the other land coverage types and the estimated surface temperature of the first land coverage type according to the temperature difference relation and the first estimated temperature difference;
obtaining target surface temperatures of the plurality of land coverage types according to the second estimated temperature difference;
according to the surface temperatures of the plurality of land cover types, obtaining a temperature difference relation between a first land cover type and the rest of land cover types in at least one land cover type comprises the following steps:
obtaining a first measured temperature difference between the measured temperature of the first land cover type and the measured temperature of the second land cover type according to the measured temperature of the first land cover type and the measured temperature of the second land cover type;
Obtaining a second measured temperature difference between the measured temperature of the first land cover type and the measured temperature of the remaining land cover type according to the measured temperature of the first land cover type and the measured temperature of the remaining land cover type;
obtaining a temperature difference relation between a first land cover type and the rest land cover types in at least one land cover type according to the first measured temperature difference and the second measured temperature difference;
wherein, the temperature difference relation is: deltaLST 1 The rest are ＝a×△LST1 Second one +b；
Wherein DeltaLST 1 The rest are For the second measured temperature difference, deltaLST 1 Second one A and b are relationship constants for a first measured temperature difference;
wherein obtaining a first estimated temperature difference between the temperature of the hybrid pel and the estimated surface temperature of the first earth coverage type by satellite measurements comprises:
obtaining an estimated surface temperature and a mixed pixel temperature of a first land coverage type through satellite measurement;
according to the estimated surface temperature and the mixed pixel temperature of the first land coverage type, the following formula is adopted:
△LST2 first one ＝LST2 Mixed pixel -LST2 First land cover type Obtaining a first estimated temperature difference;
wherein DeltaLST 2 First one For a first estimated temperature difference, LST2 Mixed pixel To mix pixel temperature, LST2 First land cover type Estimating a surface temperature for a first land cover type;
wherein obtaining a second estimated temperature difference between the estimated surface temperature of the remaining land cover types and the estimated surface temperature of the first land cover type according to the temperature difference relationship and the first estimated temperature difference comprises:
based on the temperature difference relation and the first estimated temperature difference, the method passes through the formula delta LST2 Second one ＝a×△LST2 First one +b obtaining a second estimated temperature difference;
wherein DeltaLST 2 Second one For the second estimated temperature difference, ΔLST2 First one For the first estimated temperature difference, a and b are relationship constants;
wherein obtaining the estimated surface temperature of the remaining land coverage type according to the second estimated temperature difference comprises:
according to the second estimated temperature difference, the following formula is adopted: LST2 Type of remaining land cover ＝LST2 First land cover type +△LST2 Second one Obtaining an estimated surface temperature for the remaining land coverage types;
wherein LST2 Type of remaining land cover LST2 for estimating surface temperature for remaining land cover types First land cover type For an estimated surface temperature of a first land cover type ΔLST2 Second one For a second estimated temperature difference;
According to the formula:
obtaining a target surface temperature of the second land cover type;
wherein LST3 Second land cover type LST2 for target surface temperature of second land cover type Mixed pixel To mix pixel temperature, LST2 Type of remaining land cover LST2 for estimating surface temperature for remaining land cover types First land cover type For the estimated surface temperature of the first earth coverage type, R The rest are R is the ratio of the area of the residual land coverage type to the area of the mixed pixel Second one R is the ratio of the area of the second land cover type to the area of the mixed pixel First one A ratio of an area of the first land cover type to an area of the hybrid pixel;
and obtaining the target surface temperature of the other land coverage types according to the target surface temperature of the second land coverage type and the estimated surface temperature of the first land coverage type.
2. The method of estimating a satellite remote sensing earth's surface temperature according to claim 1, wherein the mixed pixel temperature is a pixel having a highest proportion of an area of the second earth coverage type in an area of one pixel.
3. An apparatus for estimating satellite remote sensing earth surface temperature, the apparatus comprising:
The acquisition module is used for acquiring actual earth surface temperatures of various land coverage types;
the processing module is used for obtaining the temperature difference relation between the first land coverage type and the rest land coverage types in at least one land coverage type according to the actually measured surface temperatures of the plurality of land coverage types; obtaining a first estimated temperature difference between the mixed pixel temperature and the estimated surface temperature of the first earth coverage type by satellite measurement; obtaining a second estimated temperature difference between the estimated surface temperature of the other land coverage types and the estimated surface temperature of the first land coverage type according to the temperature difference relation and the first estimated temperature difference; obtaining target surface temperatures of the plurality of land coverage types according to the second estimated temperature difference; according to the surface temperatures of the plurality of land cover types, obtaining a temperature difference relation between a first land cover type and the rest of land cover types in at least one land cover type comprises the following steps:
obtaining a first measured temperature difference between the measured temperature of the first land cover type and the measured temperature of the second land cover type according to the measured temperature of the first land cover type and the measured temperature of the second land cover type;
Obtaining a second measured temperature difference between the measured temperature of the first land cover type and the measured temperature of the remaining land cover type according to the measured temperature of the first land cover type and the measured temperature of the remaining land cover type;
obtaining a temperature difference relation between a first land cover type and the rest land cover types in at least one land cover type according to the first measured temperature difference and the second measured temperature difference;
wherein the temperature differenceThe relation is: deltaLST 1 The rest are ＝a×△LST1 Second one +b；
Wherein DeltaLST 1 The rest are For the second measured temperature difference, deltaLST 1 Second one A and b are relationship constants for a first measured temperature difference;
wherein obtaining a first estimated temperature difference between the temperature of the hybrid pel and the estimated surface temperature of the first earth coverage type by satellite measurements comprises:
obtaining an estimated surface temperature and a mixed pixel temperature of a first land coverage type through satellite measurement;
according to the estimated surface temperature and the mixed pixel temperature of the first land coverage type, the following formula is adopted:
△LST2 first one ＝LST2 Mixed pixel -LST2 First land cover type Obtaining a first estimated temperature difference;
wherein DeltaLST 2 First one For a first estimated temperature difference, LST2 Mixed pixel To mix pixel temperature, LST2 First land cover type Estimating a surface temperature for a first land cover type;
wherein obtaining a second estimated temperature difference between the estimated surface temperature of the remaining land cover types and the estimated surface temperature of the first land cover type according to the temperature difference relationship and the first estimated temperature difference comprises:
based on the temperature difference relation and the first estimated temperature difference, the method passes through the formula delta LST2 Second one ＝a×△LST2 First one +b obtaining a second estimated temperature difference;
wherein DeltaLST 2 Second one For the second estimated temperature difference, ΔLST2 First one For the first estimated temperature difference, a and b are relationship constants;
wherein obtaining the estimated surface temperature of the remaining land coverage type according to the second estimated temperature difference comprises:
according to the second estimated temperature difference, the following formula is adopted: LST2 Type of remaining land cover ＝LST2 First land cover type +△LST2 Second one Obtaining an estimated surface temperature for the remaining land coverage types;
wherein LST2 Type of remaining land cover LST2 for estimating surface temperature for remaining land cover types First land cover type For an estimated surface temperature of a first land cover type ΔLST2 Second one For a second estimated temperature difference;
According to the formula:
obtaining a target surface temperature of the second land cover type;
wherein LST3 Second land cover type LST2 for target surface temperature of second land cover type Mixed pixel To mix pixel temperature, LST2 Type of remaining land cover LST2 for estimating surface temperature for remaining land cover types First land cover type For the estimated surface temperature of the first earth coverage type, R The rest are R is the ratio of the area of the residual land coverage type to the area of the mixed pixel Second one R is the ratio of the area of the second land cover type to the area of the mixed pixel First one A ratio of an area of the first land cover type to an area of the hybrid pixel;
and obtaining the target surface temperature of the other land coverage types according to the target surface temperature of the second land coverage type and the estimated surface temperature of the first land coverage type.
4. A computing device, comprising: the device comprises a processor, a memory, a communication interface and a communication bus, wherein the processor, the memory and the communication interface complete communication with each other through the communication bus;
the memory is configured to store at least one executable instruction, where the executable instruction causes the processor to perform the operations corresponding to the method for estimating satellite remote sensing earth surface temperature according to any one of claims 1 to 2.
5. A computer readable storage medium storing instructions which, when run on a computer, cause the computer to perform the method of estimating satellite remote sensing surface temperature according to any of claims 1 to 2.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110991141.4A CN113887780B (en) | 2021-08-26 | 2021-08-26 | Satellite remote sensing earth surface temperature estimation method, device and equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110991141.4A CN113887780B (en) | 2021-08-26 | 2021-08-26 | Satellite remote sensing earth surface temperature estimation method, device and equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113887780A CN113887780A (en) | 2022-01-04 |
| CN113887780B true CN113887780B (en) | 2023-11-24 |
Family
ID=79011187
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110991141.4A Active CN113887780B (en) | 2021-08-26 | 2021-08-26 | Satellite remote sensing earth surface temperature estimation method, device and equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113887780B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101187630A (en) * | 2007-12-05 | 2008-05-28 | 北京大学 | Agricultural drought monitoring method |
| CN101629850A (en) * | 2009-08-24 | 2010-01-20 | 中国农业科学院农业资源与农业区划研究所 | Method for inversing land surface temperature from MODIS data |
| CN103955606A (en) * | 2014-04-23 | 2014-07-30 | 北京大学 | Remote sensing technology-based grassland locust plague progressive prediction method |
| CN110208878A (en) * | 2019-06-14 | 2019-09-06 | 广西海佩智能科技有限公司 | Green Roof weather monitoring and tropical island effect impact evaluation method |
| CN111191673A (en) * | 2019-11-29 | 2020-05-22 | 广州地理研究所 | Ground surface temperature downscaling method and system |
| CN112487346A (en) * | 2020-10-26 | 2021-03-12 | 中国农业科学院农业资源与农业区划研究所 | Mountain land surface temperature remote sensing retrieval method |
| CN112857583A (en) * | 2021-01-14 | 2021-05-28 | 中国科学院地理科学与资源研究所 | Vegetation temperature and bare ground temperature estimation method based on satellite remote sensing data |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20210209803A1 (en) * | 2020-01-06 | 2021-07-08 | Quantela Inc | Computer-based method and system for geo-spatial analysis |
-
2021
- 2021-08-26 CN CN202110991141.4A patent/CN113887780B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101187630A (en) * | 2007-12-05 | 2008-05-28 | 北京大学 | Agricultural drought monitoring method |
| CN101629850A (en) * | 2009-08-24 | 2010-01-20 | 中国农业科学院农业资源与农业区划研究所 | Method for inversing land surface temperature from MODIS data |
| CN103955606A (en) * | 2014-04-23 | 2014-07-30 | 北京大学 | Remote sensing technology-based grassland locust plague progressive prediction method |
| CN110208878A (en) * | 2019-06-14 | 2019-09-06 | 广西海佩智能科技有限公司 | Green Roof weather monitoring and tropical island effect impact evaluation method |
| CN111191673A (en) * | 2019-11-29 | 2020-05-22 | 广州地理研究所 | Ground surface temperature downscaling method and system |
| CN112487346A (en) * | 2020-10-26 | 2021-03-12 | 中国农业科学院农业资源与农业区划研究所 | Mountain land surface temperature remote sensing retrieval method |
| CN112857583A (en) * | 2021-01-14 | 2021-05-28 | 中国科学院地理科学与资源研究所 | Vegetation temperature and bare ground temperature estimation method based on satellite remote sensing data |
Non-Patent Citations (3)
| Title |
|---|
| 土地利用/覆盖变化对城市地表参数的影响;钱乐祥;崔海山;千怀遂;夏丽华;;广州大学学报(自然科学版)(第05期);全文 * |
| 基于遥感数据的城市地表温度与土地覆盖定量研究;夏俊士;杜培军;张海荣;刘培;;遥感技术与应用(第01期);全文 * |
| 西安市城市热岛效应卫星遥感分析;王娟敏;孙娴;毛明策;杨联安;;陕西气象(第03期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113887780A (en) | 2022-01-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Konings et al. | Vegetation optical depth and scattering albedo retrieval using time series of dual-polarized L-band radiometer observations | |
| Liang et al. | Multi-factor modeling of above-ground biomass in alpine grassland: A case study in the Three-River Headwaters Region, China | |
| Su et al. | Inter-comparison of microwave satellite soil moisture retrievals over the Murrumbidgee Basin, southeast Australia | |
| Wu et al. | Evaluation of the SMAP brightness temperature downscaling algorithm using active–passive microwave observations | |
| Seguin et al. | IRSUTE: a minisatellite project for land surface heat flux estimation from field to regional scale | |
| Hu et al. | A Landsat-5 atmospheric correction based on MODIS atmosphere products and 6S model | |
| Mundava et al. | Evaluation of vegetation indices for rangeland biomass estimation in the Kimberley area of Western Australia | |
| Mokarram et al. | Modeling the behavior of Vegetation Indices in the salt dome of Korsia in North-East of Darab, Fars, Iran | |
| Zhang et al. | Evaluation of the visible and shortwave infrared drought index in China | |
| Vuolo et al. | Retrieval of Leaf Area Index from CHRIS/PROBA data: an analysis of the directional and spectral information content | |
| CN107145728A (en) | A kind of desert steppe non-acoustic detection withered grass weight remote sensing estimation method | |
| Junfeng et al. | Changes in glacier volume in the north bank of the Bangong Co Basin from 1968 to 2007 based on historical topographic maps, SRTM, and ASTER stereo images | |
| Yang et al. | Estimating sub-pixel temperatures using the triangle algorithm | |
| Harfenmeister et al. | Detecting phenological development of winter wheat and winter barley using time series of sentinel-1 and sentinel-2 | |
| Zhang et al. | Improving hydrological and vegetation modelling using regional model calibration schemes together with remote sensing data | |
| Merchant et al. | Bias correction and covariance parameters for optimal estimation by exploiting matched in-situ references | |
| JP4810604B2 (en) | Paddy rice planting situation grasping system, paddy rice cultivation situation grasping method, and paddy rice cultivation situation grasping program | |
| CN105841819A (en) | Estimation method and device of land surface temperature under cloudy condition | |
| CN113887780B (en) | Satellite remote sensing earth surface temperature estimation method, device and equipment | |
| Graf et al. | Propagating sentinel-2 top-of-atmosphere radiometric uncertainty into land surface phenology metrics using a monte carlo framework | |
| Liu et al. | Comparison of RO tropopause height based on different tropopause determination methods | |
| Yu et al. | Hyperspectral database prediction of ecological characteristics for grass species of alpine grasslands | |
| CN113641958B (en) | Method, device and equipment for processing fire information of stationary meteorological satellites | |
| Peng et al. | Uncertainties in estimating normalized difference temperature index from TOA radiances | |
| He et al. | Discovering optimal triplets for assessing the uncertainties of satellite-derived evapotranspiration products |
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 |