CN114323308B - Correction method for remote sensing city ground surface temperature angle effect - Google Patents
Correction method for remote sensing city ground surface temperature angle effect Download PDFInfo
- Publication number
- CN114323308B CN114323308B CN202210007809.1A CN202210007809A CN114323308B CN 114323308 B CN114323308 B CN 114323308B CN 202210007809 A CN202210007809 A CN 202210007809A CN 114323308 B CN114323308 B CN 114323308B
- Authority
- CN
- China
- Prior art keywords
- angle
- model
- remote sensing
- ground surface
- guta
- 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
Images
Classifications
-
- 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
- Radiation Pyrometers (AREA)
- Air Conditioning Control Device (AREA)
Abstract
A correction method for the angle effect of the ground surface temperature of a remote sensing city comprises the following steps: s1, calculating a in GUTA-T model based on computer model simulation data sw‑ig 、a w 、a sg‑ig The method comprises the steps of carrying out a first treatment on the surface of the S2, calculating a in the GUTA-T model in a satellite data inversion mode sw‑ig 、a w 、a sg‑ig The method comprises the steps of carrying out a first treatment on the surface of the S3, determining a specific form of a satellite angle effect parameterized model; s4, correcting the angle effect of the ground surface temperature of the remote sensing city. The method is suitable for remote sensing pixel scale heat radiation directivity simulation, and can effectively correct the angle effect of the urban ground surface temperature of remote sensing observation.
Description
Technical Field
The invention relates to the field of remote sensing measurement, in particular to a correction method for a remote sensing city ground surface temperature angle effect.
Background
The urban surface temperature is an important parameter reflecting urban underlying surface energy balance and local climate change, and the acquisition of the urban surface temperature information in a large range by using a thermal infrared remote sensing technology is a very natural choice. However, the heterogeneity of the three-dimensional structure and the surface temperature distribution of the urban earth surface is remarkable, so that the thermal radiation information obtained by the remote sensing sensor from different angles is greatly different, and the problem of angle effect is caused. Different types of satellite sensors may have different observation angles (for example, landsat satellites acquire ground surface temperature observed at vertical angles, MODIS sensors acquire ground surface temperature observed at vertical or inclined angles, the zenith angle of observation can reach 65 degrees), the observation angles of the same sensor passing through the same area can also change along with the passing time (for example, the sensors MODIS, AVHRR and the like carried by polar orbit satellites, and the zenith angle of observation can change from-65 degrees to 65 degrees). According to the existing research results, the angle effect of MODIS satellites reaches 8.8K and 6.5K in the building dense areas of the centers of Chicago and New York City in the United states, respectively. In the old urban area of Mantzz, france, the angle effect reaches 6.6K. This angular effect hampers inversion and application of satellite surface temperature products, mainly in: (1) reduced comparability between different remote sensing products; (2) Space-time characterization of urban surface temperature using polar orbit satellite data was hampered (Hu et al 2016; lagouarde et al 2004). Therefore, it is necessary to correct for the angular effects of the remote sensing city surface temperature.
The parameterized model is an effective means for realizing remote sensing observation angle effect correction. Compared with a complex forward simulation model, such as a computer model, a radiation transmission model and the like, the parameterized model simplifies the physical process to a certain extent, reduces input parameters and has simpler form. The parameterized model achieves good effect on natural earth surface heat radiation directivity reconstruction (Cao et al 2021; duffour et al 2016; erida et al 2018), but has shortcomings in urban earth surface application (Lagouarde and Irvine 2008; sun et al 2015). The main reason is that: urban building surfaces and vegetation surfaces have distinct geometric structures and physical properties, and the directional characteristics of heat radiation of the urban building surfaces and the vegetation surfaces are different. Considering building geometry is critical in urban heat radiation directivity simulation, as Voogt (2008) reveals that urban heat radiation directivity is closely related to building geometry and component temperature. Krayenhoff and Voogt (2016) quantitatively analyzes the relationship of the directional intensity of heat radiation to the geometry of the building based on the simulated data. The city building is simplified into cuboid which is identical in size and faces to random distribution by referring to modeling thought of building surface two-way reflection (Roujean et al 1992), and according to the geometrical relationship of a solar building-sensor, the geometrical model of building scene heat radiation directivity is developed based on geometrical optics principle by considering different component temperature differences, wherein the simulation sparsity (GUTA-spark; wang et al 2018 b), shielding (GUTA-osg; wang et al 2018 a) and dense (GUTA-dense; wang and Chen 2019) and the heat radiation directivity is expressed as a function of a solar angle and a building geometrical structure. Based on the GUTA model, the component temperature difference intra-day circulation mode is coupled by considering the time phase change of the urban heat radiation directivity, and the urban heat radiation directivity time expansion geometric model GUTA-T (Wang et al 2020) is constructed, so that the simulation of the urban heat radiation directivity intra-day and seasonal change is realized. Although the GUTA-T model gets rid of the dependence of the GUTA model on instantaneous/on-time multi-angle temperature observation, the solar zenith angle is taken as a time factor to replace a component temperature difference parameter which changes rapidly with time in the GUTA model, so that the GUTA-T model has the potential of being applied to satellite angle effect correction, and the applicability of the GUTA-T model on the scale of remote sensing pixels still needs to be discussed. Therefore, an effective correction method for remotely sensing and observing the urban earth surface temperature angle effect is still lacking at present.
Disclosure of Invention
The invention provides a correction method of a remote sensing urban ground surface temperature angle effect, which is suitable for remote sensing pixel scale heat radiation directivity simulation and can effectively correct the remote sensing urban ground surface temperature angle effect.
The technical scheme of the invention is as follows:
a correction method for the angle effect of the ground surface temperature of a remote sensing city comprises the following steps:
s1, calculating a in GUTA-T model based on computer model simulation data sw-ig 、a w 、a sg-ig Wherein the form of the GUTA-T model is:
wherein DeltaT represents the difference between the temperature in the oblique direction and the temperature in the vertical direction, k iso 、k ori And k shw Respectively represent the contribution of wall surface, orientation and ground shadow direction, a sw-ig 、a w 、a sg-ig Respectively representing the response degree of the temperature difference between the shadow wall surface and the illumination ground surface, the temperature difference between the different facing wall surfaces and the temperature difference between the shadow ground surface and the illumination ground surface to the zenith angle of the sun under the clear sky condition, and theta v Representing the zenith angle, θ observed by the sensor s Represents the zenith angle of the sun,
representing the relative azimuth angles of the sensor and the sun;
s2, calculating a in the GUTA-T model in a satellite data inversion mode sw-ig 、a w 、a sg-ig ；
S3, determining a specific form of a satellite angle effect parameterized model; and
s4, correcting the angle effect of the ground surface temperature of the remote sensing city.
Preferably, in the method for correcting the remote sensing urban surface temperature angle effect, in step S1, building geometry, physical properties, meteorological data are input into an urban surface energy balance model, three-dimensional surface temperature of a scene is simulated, average temperatures of a shadow wall surface, an illumination ground surface and the shadow ground surface are counted, and a temperature difference Δt between the shadow wall surface and the illumination ground surface is calculated sw-ig Temperature difference delta T of wall bodies with different directions w And the temperature difference delta T between the shadow ground and the illumination ground sg-ig Then to DeltaT sw-ig 、ΔT w 、ΔT sg-ig Along with the zenith angle theta of the sun s Is fitted to the variation of (a) and a is determined from the intercept of the fitting line on the y-axis sw-ig And a sg-ig Determining a based on the maximum value of the fit line in the y-direction w 。
Preferably, in the method for correcting the earth surface temperature angle effect of the remote sensing city, in step S2, temperature observation data of a plurality of moments of a research city are input into a GUTA-T model, the change of the satellite angle effect along with the observation angle is calculated, and the GUTA-T model is calibrated by utilizing remote sensing multi-angle temperature observation; simultaneously inputting solar and observation angle data and building geometric structure parameters into a GUTA-T model; finally, a is obtained by fitting a GUTA-T model sw-ig 、a w 、a sg-ig Is a value of (2).
Preferably, in the correction method of the remote sensing city ground surface temperature angle effect, in step S2, the sun and observation angle data includes an average sun zenith angle, an observation zenith angle, and a relative azimuth angle within each observation zenith angle interval; and building geometry parameters include the ratio of roof to building average footprint and the ratio of building height to street width.
Preferably, in the method for correcting the angle effect of the ground surface temperature of the remote sensing city, in step S3, the coupling step S1 calculates based on the computer model simulation dataThe a obtained sw-ig 、a w 、a sg-ig The value and a calculated by means of satellite data inversion in step S2 sw-ig 、a w 、a sg-ig The value is optimized to obtain the final a sw-ig 、a w 、a sg-ig And (5) determining the value, further determining the specific form of the satellite angle effect parameterized model, and obtaining an optimized GUTA-T model.
Preferably, in the method for correcting the angle effect of the ground surface temperature of the remote sensing city, in step S4, multi-angle remote sensing data at a plurality of times is input into the satellite angle effect parameterized model obtained in step S3, and the temperature difference between the remote sensing direction observation and the vertical direction is calculated, so that the remote sensing direction observation is corrected to the vertical direction, and the angle effect correction is realized.
According to the technical scheme of the invention, the beneficial effects are that:
the geometric model of the component of the method is suitable for simulating the heat radiation directivity of the remote sensing pixel scale, can be well suitable for simulating the earth surface temperature angle effect of a remote sensing city, and simulates the change characteristics of the angle effect along with the city latitude and the building geometric structure; the method is suitable for remote sensing inversion, and can realize the correction of the angle effect under the condition of known remote sensing multi-angle temperature observation.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below.
FIG. 1 is a flow chart of a method for correcting the angle effect of the surface temperature of a remote sensing city according to the present invention.
FIG. 2 is a flow chart of a specific embodiment of a method for calculating the satellite-observed urban surface temperature angle effect according to the present invention;
fig. 3 is a graph of the verification effect of the method of the present invention in different cities.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention more clear, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention.
The invention discloses an effective correction method for the angle effect of the ground surface temperature of a remote sensing city, which comprises the following steps:
s1, calculating a in GUTA-T model based on computer model simulation data sw-ig 、a w 、a sg-ig Wherein the form of the GUTA-T model is:
wherein DeltaT represents the difference between the temperature in the oblique direction and the temperature in the vertical direction, k iso 、k ori And k shw Is the "kernel" of the model, is a function of the sun-building-sensor geometry, representing the contributions of wall surface, orientation and ground shadow orientation, a sw-ig 、a w 、a sg-ig Representing the response degree of the temperature difference between the shadow wall surface and the illumination ground surface, the temperature difference between the different facing wall surfaces and the temperature difference between the shadow ground surface and the illumination ground surface to the zenith angle of the sun under the clear sky condition, and theta v Representing the zenith angle, θ observed by the sensor s Represents the zenith angle of the sun,
representing the relative azimuth angles of the sensor and the sun;
s2, calculating a in the GUTA-T model in a satellite data inversion mode sw-ig 、a w 、a sg-ig ；
S3, determining a specific form of a satellite angle effect parameterized model; and
s4, correcting the angle effect of the ground surface temperature of the remote sensing city.
The geometric model of the component of the method is suitable for simulating the heat radiation directivity of the remote sensing pixel scale, can be well suitable for simulating the earth surface temperature angle effect of a remote sensing city, and simulates the change characteristics of the angle effect along with the city latitude and the building geometric structure; the method is suitable for remote sensing inversion, and can realize the correction of the angle effect under the condition of known remote sensing multi-angle temperature observation.
The principle of the correction method of the remote sensing city surface temperature angle effect of the invention is as follows: on the basis of the GUTA-T model, the characteristics of model parameters on the remote sensing pixel scale are discussed, and a novel geometric model suitable for the remote sensing pixel scale heat radiation directivity simulation, namely an optimized GUTA-T model, is constructed. Wherein, the form of the GUTA-T model is as follows:
wherein DeltaT represents the difference between the temperature in the oblique direction and the temperature in the vertical direction, k iso 、k ori And k shw Is the "kernel" of the model, is a function of the solar-building-sensor geometry, and represents the contributions of wall surface, orientation and ground shadow orientation, θ, respectively v Representing the zenith angle, θ observed by the sensor s Represents the zenith angle of the sun,
representing the relative azimuth angle of the sensor and the sun, a sw-ig 、a w 、a sg-ig The response degree of component temperature differences (corresponding to the temperature differences of the shadow wall surface and the illumination ground, the temperature differences of the different facing wall surfaces and the temperature differences of the shadow ground and the illumination ground) to solar zenith angles under the clear sky condition is characterized, and the response degree is related to meteorological conditions, building geometric structures, vegetation coverage and material properties, and the values of three parameters are influenced by changing the shadow shielding relation and the temperature distribution between the surfaces of the building geometric structures and the vegetation coverage. The material properties affect the response of the material surface temperature to incident solar radiation. According to Krayenhoff and Voogt (2016), the building geometry has a greater effect on the directionality of heat radiation in urban areas than the material properties.
According to formula (1), if the influence of external weather conditions is ignored, a sw-ig 、a w 、a sg-ig In relation to the inherent properties of the earth's surface, it is assumed that the three parameter values remain unchanged for a short period of time. a, a sw-ig 、a w 、a sg-ig Characterization ofThe response degree of the component temperature difference to the zenith angle of the sun under the clear sky condition is calculated according to the following formula:
from the formula, a sw-ig 、a w 、a sg-ig Respectively equal to the solar zenith angle theta s At 0 DEG, the temperature difference delta T between the shadow wall surface and the illumination ground sw-ig Temperature difference delta T of wall bodies with different directions w And the temperature difference delta T between the shadow ground and the illumination ground sg-ig . Because the heat radiation directivity of the whole scene is the result of the temperature difference action of the average component in the scene, for a building scene with strong temperature distribution heterogeneity (such as non-unique illumination ground temperature), the average value of the temperature difference of the component in the scene is adopted to calculate a sw-ig 、a w 、a sg-ig Is a value of (2). If GUTA-T model is applied to remote sensing pixel scale, a sw-ig 、a w 、a sg-ig The average value of the component temperature differences corresponding to the pixel scale may be related to the average building geometry, vegetation coverage at the pixel scale. Therefore, a can be obtained by forward calculation based on the simulation result of the computer model sw-ig 、a w 、a sg-ig The method comprises the steps of carrying out a first treatment on the surface of the The GUTA-T model can be calibrated by utilizing remote sensing multi-angle temperature observation, and a is obtained by inversion sw-ig 、a w 、a sg-ig Is a value of (2). Thus, a is determined by two ways of computer simulation and multi-angle satellite data inversion sw-ig 、a w 、a sg-ig And determining a specific form of a satellite angle effect parameterized model applicable to remote sensing pixel scale angle effect simulation.
Based on the principle, the correction method of the remote sensing city surface temperature angle effect mainly comprises the following steps (see fig. 1 and 2):
s1, calculating a in GUTA-T model based on computer model simulation data sw-ig 、a w 、a sg-ig (forward process in FIG. 2), wherein the form of the GUTA-T model is as in equation (1) above. The computer model herein refers to an urban surface energy balance model for simulating urban surface temperature changes. Taking TUF3D model as an example, building geometry, physical properties and meteorological data are input into an energy balance model, three-dimensional surface temperature of a scene is simulated, average temperatures of a shadow wall surface, an illumination ground and the shadow ground are counted, and temperature difference delta T between the shadow wall surface and the illumination ground is calculated sw-ig Temperature difference delta T of wall bodies with different directions w And the temperature difference delta T between the shadow ground and the illumination ground sg-ig . Then to DeltaT sw-ig 、ΔT w 、ΔT sg-ig Along with the zenith angle theta of the sun s Is fitted to the variation of (a) and a is determined from the intercept of the fitting line on the y-axis sw-ig And a sg-ig Determining a based on the maximum value of the fit line in the y-direction w 。
S2, calculating a in the GUTA-T model in a satellite data inversion mode sw-ig 、a w 、a sg-ig (reverse process in fig. 2). Selecting a research city, inputting temperature observation data (at least two angles are needed to be known at each moment) at a plurality of moments into a GUTA-T model, calculating the change of satellite angle effect along with the observation angle by using an existing satellite angle effect calculation method (Hu et al 2016), and calibrating the GUTA-T model by using remote sensing multi-angle temperature observation. At the same time, known parameters are input into the GUTA-T model, wherein the known parameters comprise: (1) sun and observation angle data: average solar zenith angles, observation zenith angles and relative azimuth angles within each observation zenith angle interval; (2) Building geometry parameters (i.e., building structure parameters in fig. 2): roof to building average floor area ratio lambda p And the ratio of building height to street width h/w. Finally, a is obtained by fitting a GUTA-T model sw-ig 、a w 、a sg-ig Is a value of (2).
S3, determining a specific form of the satellite angle effect parameterized model. Coupling a calculated based on computer model simulation data and calculated by means of satellite data inversion sw-ig 、a w 、a sg-ig The value is optimized to obtain the final a sw-ig 、a w 、a sg-ig And (5) determining the value, further determining the specific form of the satellite angle effect parameterized model, and obtaining an optimized GUTA-T model.
S4, correcting the angle effect of the ground surface temperature of the remote sensing city. And (3) inputting the multi-angle remote sensing data at a plurality of moments into the optimized GUTA-T model (namely, the satellite angle effect parameterized model) obtained in the step (S3), calculating the temperature difference between the remote sensing direction observation and the vertical direction, and correcting the remote sensing direction observation to the vertical direction, so as to realize angle effect correction and obtain the hemispherical space temperature directivity at continuous moments.
Fig. 3 is a graph of the verification effect of the method of the present invention in different cities. Taking 6 cities as an example, from the south to the north, the latitudes are respectively (9.97 degrees, 53.62 degrees), the latitudes and longitudes are respectively (9.15 degrees, 45.7 degrees), the latitudes (114.33 degrees, 30.60 degrees), the latitudes (100.45 degrees, 13.56 degrees), the latitudes (40.30 degrees, the latitudes of-20.23 degrees), the latitudes (144.95 degrees and the latitudes of-37.86 degrees) in the northern hemisphere, and the latitudes are respectively from the top to the bottom in the northern hemisphere. The open low-rise building group (open low-rise) is most commonly distributed in a global city according to the local climate zone LCZ (local climate zone, stewart and ike 2012) dataset of the global city downloaded from WUDAPT. To obtain the pure pixels of LCZ, the MODIS data is masked by LCZ data, and MODIS pixels mainly comprising the wide and low building groups are screened out. The angle effect was extracted using MODIS LST data for 5-9 months of 10 years. According to (Stewart and like 2012), lambda is used p As average building structure of the open low building group, GUTA-T model was input =0.3 and h/w=0.5. Fitting the MODIS angle effect extraction result of Milan region with GUTA-T model to obtain a sw-ig 、a w And a sg-ig respectively-17.46K, 7.97K and-19.46K. The set of parameter values is input into a GUTA-T model in 6 cities, and the sensor and solar angle data of each city are combined to simulate MODIS LS of different citiesT angle effect. As shown in FIG. 3, the polarization diagram shows the relative geometrical relationship between the sensor observation angle and the sun angle during the satellite transit, and the relationship determines the angle effect of satellite observation and changes with latitude. In the summer of northern hemisphere, the closer the city is to the equator (such as Mangu), the closer the sun position at the satellite transit time is to the zenith direction, the closer the "hot spot" is to the zenith direction, the directivity is symmetrically distributed along with the observation angle, and the larger the influence of the angle effect on satellite observation is. 5-9 months is winter in southern hemisphere city, and the plane where the sensor is observed is vertical to the main plane of the sun. Compared with the northern hemisphere, the urban satellite observation in the southern hemisphere is less affected by the angle effect. In general, GUTA-T can accurately simulate the change rule of satellite observation angle effect along with latitude. The average root mean square error RMSD for 6 cities is 0.62K, while the average effective directivity is 3.31K.
As can be seen from the verification effect, the method provided by the invention can be well applied to simulation of the remote sensing city surface temperature angle effect, and the simulation angle effect is characterized by the change of city latitude and building geometry; the method is suitable for remote sensing inversion, and can realize the correction of the angle effect under the condition of known remote sensing multi-angle temperature observation.
The above examples are only specific embodiments of the present invention for illustrating the technical solution of the present invention, but not for limiting the scope of the present invention, and although the present invention has been described in detail with reference to the foregoing examples, it will be understood by those skilled in the art that the present invention is not limited thereto: any person skilled in the art may modify or easily conceive of the technical solution described in the foregoing embodiments, or perform equivalent substitution of some of the technical features, while remaining within the technical scope of the present disclosure; such modifications, changes or substitutions do not depart from the spirit and scope of the technical solutions of the embodiments of the present invention, and are intended to be included in the scope of the present invention. Therefore, the protection scope of the invention is subject to the protection scope of the claims.
Reference to the literature
Cao,B.,Roujean,J.-L.,Gastellu-Etchegorry,J.-P.,Liu,Q.,Du,Y.,Lagouarde,J.-P.,Huang,H.,Li,H.,Bian,Z.,Hu,T.,Qin,B.,Ran,X.,&Xiao,Q.(2021).A general framework of kernel-driven modeling in the thermal infrared domain.Remote Sensing of Environment,252,112157
Duffour,C.,Lagouarde,J.P.,&Roujean,J.L.(2016).A two parameter model to simulate thermal infrared directional effects for remote sensing applications.Remote Sensing of Environment,186,250-261
Ermida,S.L.,Trigo,I.F.,DaCamara,C.C.,&Roujean,J.-L.(2018).Assessing the potential of parametric models to correct directional effects on local to global remotely sensed LST.Remote Sensing of Environment,209,410-422
Hu,L.,Monaghan,A.,Voogt,J.A.,&Barlage,M.(2016).A first satellite-based observational assessment of urban thermal anisotropy.Remote Sensing of Environment,181,111-121
Krayenhoff,E.S.,&Voogt,J.A.(2016).Daytime thermal anisotropy of urban neighbourhoods:Morphological causation.Remote Sensing,8,108
Lagouarde,J.-P.,Moreau,P.,Irvine,M.,Bonnefond,J.-M.,Voogt,J.A.,&Solliec,F.(2004).Airborne experimental measurements of the angular variations in surface temperature over urban areas:case study of Marseille(France).Remote Sensing of Environment,93,443-462
Lagouarde,J.P.,&Irvine,M.(2008).Directional anisotropy in thermal infrared measurements over Toulouse city centre during the CAPITOUL measurement campaigns:first results.Meteorology and Atmospheric Physics,102,173-185
Roujean,J.L.,Leroy,M.,&Deschamps,P.Y.(1992).A bidirectional reflectance model of the Earth's surface for the correction of remote sensing data.Journal of Geophysical Research,97,20455-20468
Stewart,I.D.,&Oke,T.R.(2012).Local climate zones for urban temperature studies.Bulletin of the American Meteorological Society,93,1879-1900
Voogt,J.A.(2008).Assessment of an urban sensor wiew model for thermal anisotropy.Remote Sensing of Environment,112,482-495
Wang,D.,&Chen,Y.(2019).A Geometric Model to Simulate Urban Thermal Anisotropy in Simplified Dense Neighborhoods(GUTA-Dense).IEEE Transactions on Geoscience and Remote Sensing,57,6226-6239
Wang,D.,Chen,Y.,Cui,Y.,&Sun,H.(2018a).A Geometric Model to Simulate Urban Thermal Anisotropy for Simplified Neighborhoods.IEEE Transactions on Geoscience and Remote Sensing,56,4930-4944
Wang,D.,Chen,Y.,Voogt,J.A.,Krayenhoff,E.S.,Wang,J.,&Wang,L.(2020).An advanced geometric model to simulate thermal anisotropy time-series for simplified urban neighborhoods(GUTA-T).Remote Sensing of Environment,237,111547
Wang,D.,Chen,Y.,&Zhan,W.(2018b).A geometric model to simulate thermal anisotropy over a sparse urban surface(GUTA-sparse).Remote Sensing of Environment,209,263-274
Sun, chen Yunhao, wen Feng, wang Mengjie, ma Wei (2015) A model for driving directional variation of infrared equivalent emissivity of urban surface heat with high emissivity and uncertainty analysis thereof, and infrared and millimeter wave school report, 34,66-73
Claims (6)
1. A correction method for the angle effect of the ground surface temperature of a remote sensing city is characterized by comprising the following steps:
s1, calculating a in GUTA-T model based on computer model simulation data sw-ig 、a w 、a sg-ig Wherein the form of the GUTA-T model is:
wherein DeltaT represents the difference between the temperature in the oblique direction and the temperature in the vertical direction, k iso 、k ori And k shw Respectively represent the contribution of wall surface, orientation and ground shadow direction, a sw-ig 、a w 、a sg-ig Respectively representing the response degree of the temperature difference between the shadow wall surface and the illumination ground surface, the temperature difference between the different facing wall surfaces and the temperature difference between the shadow ground surface and the illumination ground surface to the zenith angle of the sun under the clear sky condition, and theta v Representing the zenith angle, θ observed by the sensor s Represents the zenith angle of the sun,
representing the relative azimuth angles of the sensor and the sun;
s2, calculating a in the GUTA-T model in a satellite data inversion mode sw-ig 、a w 、a sg-ig ；
S3, determining a specific form of a satellite angle effect parameterized model; and
s4, correcting the angle effect of the ground surface temperature of the remote sensing city.
2. The method according to claim 1, wherein in step S1, building geometry, physical properties, meteorological data are input into a city ground surface energy balance model, three-dimensional surface temperatures of a scene are simulated, average temperatures of a shadow wall surface, an illumination ground surface and the shadow ground surface are counted, and a temperature difference Δt between the shadow wall surface and the illumination ground surface is calculated sw-ig Temperature difference delta T of wall bodies with different directions w And the temperature difference delta T between the shadow ground and the illumination ground sg-ig Then to DeltaT sw-ig 、ΔT w 、ΔT sg-ig Along with the zenith angle theta of the sun s Is fitted to the variation of (a) and a is determined from the intercept of the fitting line on the y-axis sw-ig And a sg-ig Determining a based on the maximum value of the fit line in the y-direction w 。
3. The method according to claim 1, wherein in step S2, temperature observation data of a study city at several moments are input into the GUTA-T model, and calculatedThe satellite angle effect changes along with the observation angle, and the GUTA-T model is calibrated by utilizing remote sensing multi-angle temperature observation; simultaneously inputting solar and observation angle data and building geometric structure parameters into the GUTA-T model; finally, fitting the GUTA-T model to obtain a sw-ig 、a w 、a sg-ig Is a value of (2).
4. A method of correcting for the effects of ground surface temperature angles in a remote sensing city as claimed in claim 3, wherein the sun and observation angle data comprises an average sun zenith angle, an observation zenith angle and a relative azimuth angle within each observation zenith angle interval; and the building geometry parameters include a roof to building average footprint ratio and a building height to street width ratio.
5. The method for correcting the angle effect of the ground surface temperature of the remote sensing city according to claim 1, wherein in the step S3, a calculated based on the computer model simulation data in the step S1 is coupled sw-ig 、a w 、a sg-ig The value and a calculated in step S2 by means of inversion of said satellite data sw-ig 、a w 、a sg-ig The value is optimized to obtain the final a sw-ig 、a w 、a sg-ig And (3) determining the specific form of the satellite angle effect parameterized model to obtain an optimized GUTA-T model.
6. The method for correcting the angle effect of the ground surface temperature of the remote sensing city according to claim 1, wherein in the step S4, the multi-angle remote sensing data at a plurality of moments are input into the satellite angle effect parameterized model obtained in the step S3, the temperature difference between the remote sensing direction observation and the vertical direction is calculated, and the remote sensing direction observation is corrected to the vertical direction, so that the angle effect correction is realized.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210007809.1A CN114323308B (en) | 2022-01-06 | 2022-01-06 | Correction method for remote sensing city ground surface temperature angle effect |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210007809.1A CN114323308B (en) | 2022-01-06 | 2022-01-06 | Correction method for remote sensing city ground surface temperature angle effect |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114323308A CN114323308A (en) | 2022-04-12 |
| CN114323308B true CN114323308B (en) | 2023-05-02 |
Family
ID=81023978
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210007809.1A Active CN114323308B (en) | 2022-01-06 | 2022-01-06 | Correction method for remote sensing city ground surface temperature angle effect |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114323308B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2016106950A1 (en) * | 2014-12-30 | 2016-07-07 | 华中科技大学 | Zonal underground structure detection method based on sun illumination and shade compensation |
| CN108009392A (en) * | 2017-10-19 | 2018-05-08 | 桂林航天工业学院 | A kind of the Remote Sensing Reflectance model construction and Analysis method of dense vegetation earth's surface |
| CN108320285A (en) * | 2018-02-07 | 2018-07-24 | 中国地质大学（武汉） | Urban wetland tropical island effect analysis method based on multi-source Remote Sensing Images and system |
| WO2018145229A1 (en) * | 2017-02-10 | 2018-08-16 | 广西壮族自治区气象减灾研究所 | Accurate large-area inversion method for near-surface air temperature |
| CN109211960A (en) * | 2018-08-10 | 2019-01-15 | 北京师范大学 | A kind of dense city building scene thermal emission directionality strength calculation method |
| CN109446693A (en) * | 2018-11-08 | 2019-03-08 | 北京师范大学 | A kind of time-scale extension method of urban architecture scene thermal emission directionality intensity |
-
2022
- 2022-01-06 CN CN202210007809.1A patent/CN114323308B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2016106950A1 (en) * | 2014-12-30 | 2016-07-07 | 华中科技大学 | Zonal underground structure detection method based on sun illumination and shade compensation |
| WO2018145229A1 (en) * | 2017-02-10 | 2018-08-16 | 广西壮族自治区气象减灾研究所 | Accurate large-area inversion method for near-surface air temperature |
| CN108009392A (en) * | 2017-10-19 | 2018-05-08 | 桂林航天工业学院 | A kind of the Remote Sensing Reflectance model construction and Analysis method of dense vegetation earth's surface |
| CN108320285A (en) * | 2018-02-07 | 2018-07-24 | 中国地质大学（武汉） | Urban wetland tropical island effect analysis method based on multi-source Remote Sensing Images and system |
| CN109211960A (en) * | 2018-08-10 | 2019-01-15 | 北京师范大学 | A kind of dense city building scene thermal emission directionality strength calculation method |
| CN109446693A (en) * | 2018-11-08 | 2019-03-08 | 北京师范大学 | A kind of time-scale extension method of urban architecture scene thermal emission directionality intensity |
Non-Patent Citations (2)
| Title |
|---|
| 城市典型下垫面完全表面温度遥感估算研究;姜璐;《中国优秀硕士学位论文全文数据库》;全文 * |
| 热红外核驱动模型在垄行植被热辐射方向性特征拟合中的精度分析;曹彪;《遥感学报》;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114323308A (en) | 2022-04-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Lagouarde et al. | Airborne experimental measurements of the angular variations in surface temperature over urban areas: case study of Marseille (France) | |
| Wang et al. | Comparing the diurnal and seasonal variabilities of atmospheric and surface urban heat islands based on the Beijing urban meteorological network | |
| Lagouarde et al. | Modelling daytime thermal infrared directional anisotropy over Toulouse city centre | |
| Francois et al. | Analytical parameterization of canopy directional emissivity and directional radiance in the thermal infrared. Application on the retrieval of soil and foliage temperatures using two directional measurements | |
| Guillevic et al. | Land Surface Temperature product validation using NOAA's surface climate observation networks—Scaling methodology for the Visible Infrared Imager Radiometer Suite (VIIRS) | |
| Hu et al. | How can we use MODIS land surface temperature to validate long‐term urban model simulations? | |
| Lagouarde et al. | Experimental characterization and modelling of the nighttime directional anisotropy of thermal infrared measurements over an urban area: Case study of Toulouse (France) | |
| Cheng et al. | Estimating the hemispherical broadband longwave emissivity of global vegetated surfaces using a radiative transfer model | |
| Kolassa et al. | Soil moisture retrieval from AMSR-E and ASCAT microwave observation synergy. Part 1: Satellite data analysis | |
| Duffour et al. | A two parameter model to simulate thermal infrared directional effects for remote sensing applications | |
| Wang et al. | Modeling the angular effect of MODIS LST in urban areas: A case study of Toulouse, France | |
| Zheng et al. | The retrieved urban LST in Beijing based on TM, HJ-1B and MODIS | |
| Ye et al. | Standing water effect on soil moisture retrieval from L-band passive microwave observations | |
| CN101936777A (en) | Method for inversing air temperature of surface layer based on thermal infrared remote sensing | |
| Yang et al. | A semi-empirical method for estimating complete surface temperature from radiometric surface temperature, a study in Hong Kong city | |
| Yang et al. | Assessing the impact of urban geometry on surface urban heat island using complete and nadir temperatures | |
| Wang et al. | An advanced geometric model to simulate thermal anisotropy time-series for simplified urban neighborhoods (GUTA-T) | |
| Firozjaei et al. | Impact of surface characteristics and their adjacency effects on urban land surface temperature in different seasonal conditions and latitudes | |
| Peischl et al. | Wheat canopy structure and surface roughness effects on multiangle observations at L-band | |
| Wang et al. | Urban thermal anisotropy: A comparison among observational and modeling approaches at different time scales | |
| Jensen et al. | Satellite remote sensing of urban heat islands: current practice and prospects | |
| Gao et al. | Spatial downscaling based on spectrum analysis for soil freeze/thaw status retrieved from passive microwave | |
| CN114323308B (en) | Correction method for remote sensing city ground surface temperature angle effect | |
| Huang et al. | Thermal emission hot-spot effect of crop canopies—Part I: Simulation | |
| Si et al. | An artificial neuron network with parameterization scheme for estimating net surface shortwave radiation from satellite data under clear sky—Application to simulated GF-5 data set |
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 |