CN115628814A - Urban surface thermal effect space-time measurement method, terminal device and storage medium - Google Patents
Urban surface thermal effect space-time measurement method, terminal device and storage medium Download PDFInfo
- Publication number
- CN115628814A CN115628814A CN202211259580.7A CN202211259580A CN115628814A CN 115628814 A CN115628814 A CN 115628814A CN 202211259580 A CN202211259580 A CN 202211259580A CN 115628814 A CN115628814 A CN 115628814A
- Authority
- CN
- China
- Prior art keywords
- index
- surface heat
- area
- temperature
- research area
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 230000000694 effects Effects 0.000 title claims abstract description 86
- 238000000691 measurement method Methods 0.000 title description 4
- 238000011160 research Methods 0.000 claims abstract description 65
- 238000000034 method Methods 0.000 claims abstract description 28
- 238000004364 calculation method Methods 0.000 claims description 30
- 238000004590 computer program Methods 0.000 claims description 17
- 238000005259 measurement Methods 0.000 claims description 15
- 238000004458 analytical method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 3
- 239000010931 gold Substances 0.000 description 3
- 229910052737 gold Inorganic materials 0.000 description 3
- 230000006870 function Effects 0.000 description 2
- 230000036541 health Effects 0.000 description 2
- 238000005192 partition Methods 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 238000012300 Sequence Analysis Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000012732 spatial analysis Methods 0.000 description 1
- 230000004083 survival effect Effects 0.000 description 1
- 238000013316 zoning Methods 0.000 description 1
Images
Classifications
-
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/13—Satellite images
-
- 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
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Astronomy & Astrophysics (AREA)
- Remote Sensing (AREA)
- Multimedia (AREA)
- Theoretical Computer Science (AREA)
- Radiation Pyrometers (AREA)
Abstract
The invention relates to a method for measuring time and space of urban earth surface thermal effect, a terminal device and a storage medium, wherein the method comprises the following steps: s1: based on the thermal infrared remote sensing data, performing inversion to obtain surface temperature grid data which are continuous in the range space of the research area; s2: calculating the average value and standard deviation of the surface temperature in the research area range, and identifying a high-temperature area and a secondary high-temperature area of the research area according to the surface temperature range set by the average value and standard deviation; s3: calculating the surface heat amplitude index and the surface heat intensity index of the research area based on the areas of the identified high-temperature area and the second high-temperature area; s4: and carrying out weighted summation on the earth surface heat amplitude index and the earth surface heat intensity index of the research area to obtain the earth surface heat effect index of the research area. The invention quantitatively evaluates the range (thermal amplitude) and the degree (thermal intensity) of the regional earth surface heat influence by the standardized index, and can generate a spatially continuous thermodynamic index map layer, thereby more comprehensively reflecting the spatial heterogeneity characteristics of the urban earth surface heat environment.
Description
Technical Field
The invention relates to the technical application field of remote sensing geographic information in urban ecology, in particular to a method for measuring the spatiotemporal measure of urban earth surface heat effect, terminal equipment and a storage medium.
Background
An urban surface thermal environment is a heat-related physical environment that can have a significant impact on the human body's cooling and heating feelings, health level, and resident survival development, and is generally characterized by a surface temperature. With the continuous promotion of the urbanization process and the continuous growth of urban population, a large number of natural ground surfaces inside cities are replaced by artificial ground surfaces, the updating of the underlying surfaces, particularly the sharp increase of impervious surfaces, the local or overall ground surface temperature of the cities is improved by changing the balance of ground surface heat radiation transmission, and the city ground surface heat effect is aggravated. The aggravation of urban surface heat effect is the main reason and important expression of urban surface heat environment deterioration, wherein the most typical heat island effect is closely related to urban resident health, social and economic development, natural ecological environment protection and the like.
The quantitative evaluation of the urban surface thermal effect is an important basis for the research of urban surface thermal environment, and scholars have made certain attempts to the quantitative measurement of the urban surface thermal effect, such as a three-dimensional dynamic thermal effect numerical model established based on the difference of thermal properties of urban underlying surfaces, or thermal effects of different types of tables calculated by using the difference of average surface temperature, but the current quantitative index related to the surface thermal effect still has the following disadvantages: in the time dimension, because of the problems of 'instantaneity' characteristic and 'high time and spatial resolution which are difficult to unify' of the thermal radiation information obtained by remote sensing data, a thermal effect index obtained directly on the basis of earth surface temperature calculation has great uncertainty, so that time sequence representativeness is lacked, and the long-time dynamic evolution of the urban earth surface thermal effect cannot be effectively analyzed; although some indexes utilize the average earth surface temperature to eliminate uncertainty in a time sequence in the calculation process, most of the indexes are based on classification statistics of land utilization or coverage, so that only the heat effect value of the whole corresponding ground class can be obtained, and therefore, in the spatial dimension, the spatial heterogeneity of the earth surface heat effect of the cities cannot be fully identified, so that the earth surface heat effect comparison between different areas in the same city cannot be effectively measured and analyzed, and the comparison and analysis among the cities are not facilitated due to the earth surface structure difference of different cities, which brings certain difficulty for fully understanding the internal and external differences of the urban heat environment characteristics and the influence factors thereof.
In general, no general index that sufficiently reflects the meaning of the terrestrial heat effect has been found, and the representativeness and the availability of the relevant index are to be improved. Therefore, it is urgently needed to develop an efficient and general multi-scale space-time measuring method for urban ground surface thermal effect, which is not affected by the difference of the time-space dimensions.
Disclosure of Invention
In order to solve the problems, the invention provides a method for measuring the urban ground surface thermal effect space-time, a terminal device and a storage medium.
The specific scheme is as follows:
a space-time measuring method for urban surface thermal effect comprises the following steps:
s1: based on the thermal infrared remote sensing data, performing inversion to obtain continuous ground surface temperature grid data in the range space of the research area;
s2: calculating the average value and standard deviation of the surface temperature in the research area range, and identifying a high-temperature area and a secondary high-temperature area of the research area according to the surface temperature range set by the average value and standard deviation;
s3: calculating the surface heat amplitude index and the surface heat intensity index of the research area based on the areas of the identified high-temperature area and the second high-temperature area;
s4: and carrying out weighted summation on the earth surface heat amplitude index and the earth surface heat intensity index of the research area to obtain the earth surface heat effect index of the research area.
Further, in the remote sensing image data used in the step S1, the research area range should be free of cloud occlusion.
Further, in step S2, the range of the surface temperature zone corresponding to the high temperature zone is: t is s Mu + std; surface temperature range corresponding to secondary high temperature areaComprises the following steps: mu +0.5std < T s Less than or equal to mu + std; wherein, T s Represents the surface temperature, μ represents the mean of the surface temperature, and std represents the standard deviation of the surface temperature.
Further, the surface heat amplitude index in step S3 includes a comprehensive surface heat amplitude index corresponding to the whole research area and a spatial heat amplitude index corresponding to each grid pixel in the research area, and the surface heat intensity index includes a comprehensive surface heat intensity index corresponding to the whole research area and a spatial heat intensity index corresponding to each grid pixel in the research area; in the calculation of the spatial heat amplitude index and the spatial heat intensity index, pixel-by-pixel scanning calculation is carried out by a moving window method, and the integral calculation value of the square window is assigned to a window center pixel;
the calculation formula of the comprehensive earth surface heat amplitude index is as follows:
the calculation formula of the comprehensive earth surface heat intensity index is as follows:
the calculation formula of the spatial surface heat amplitude index is as follows:
the formula for calculating the spatial surface heat intensity index is as follows:
wherein HAI represents the comprehensive surface heat amplitude index of the whole research area, HII represents the comprehensive surface heat intensity index of the whole research area, and S htz And S shtz Respectively showing the areas of the high-temperature area and the sub-high-temperature area of the research area, S showing the total area of the research area and the heat amplitude index HAI of the space earth surface i Representing the surface heat amplitude index and the spatial surface heat intensity index HII of the ith pixel element i Represents the surface heat intensity index of the ith pixel element,
indicates the area of the high-temperature region contained in the window with the pixel i as the center,
denotes the area of the sub-high temperature region contained in the window centered on the pixel i, S i Representing the total area of the window centred on the picture element i.
Further, the surface heat effect index in the step S4 includes a comprehensive surface heat effect index corresponding to the whole research area and a spatial surface heat effect index corresponding to each grid pixel in the research area;
the calculation formula of the comprehensive earth surface heat effect index is as follows: HEI = α HAI + β HII;
the calculation formula of the spatial earth surface heat effect index is as follows: HEI i ＝αHAI i +βHII i ；
Wherein HEI represents the comprehensive earth surface heat effect index of the whole research area, and the space earth surface heat effect index HEI i The surface thermal effect index of the ith pixel element is shown, alpha and beta both represent weight coefficients and alpha + beta =1.
The terminal equipment for measuring the urban ground surface thermal effect space-time comprises a processor, a memory and a computer program which is stored in the memory and can run on the processor, wherein the processor executes the computer program to realize the steps of the method of the embodiment of the invention.
A computer-readable storage medium, in which a computer program is stored which, when being executed by a processor, carries out the steps of the method as described above for embodiments of the invention.
By adopting the technical scheme, the method overcomes the defect that the existing correlation index is used for describing the spatial distribution characteristics of the regional earth surface heat effect, solves the problem that the spatial contrast analysis is difficult when the regional earth surface heat effect is used, realizes scientific and quantitative effective measurement of the urban earth surface heat effect, and visually reveals the time evolution and the spatial distribution characteristics of the urban earth surface heat environment.
Drawings
Fig. 1 is a flowchart illustrating a first embodiment of the present invention.
Fig. 2 is a diagram illustrating a surface temperature inversion result according to a first embodiment of the invention.
FIG. 3 is a chart showing the results of zoning of the surface temperature according to the first embodiment of the present invention.
Fig. 4 is a diagram showing the result of the spatial distribution of the surface heat amplitude according to the first embodiment of the present invention.
Fig. 5 is a diagram showing the results of the spatial distribution of the surface heat intensity according to the first embodiment of the present invention.
Fig. 6 is a diagram illustrating the result of the spatial distribution of the terrestrial thermal effect according to the first embodiment of the present invention.
Detailed Description
To further illustrate the various embodiments, the invention provides the accompanying drawings. The accompanying drawings, which are incorporated in and constitute a part of this disclosure, illustrate embodiments of the invention and, together with the description, serve to explain the principles of the embodiments. Those skilled in the art will appreciate still other possible embodiments and advantages of the present invention with reference to these figures.
The invention will now be further described with reference to the accompanying drawings and detailed description.
The first embodiment is as follows:
the embodiment of the invention provides a space-time measurement method for urban earth surface thermal effect, which comprises the following steps as shown in figure 1:
s1: and (4) performing inversion to obtain continuous earth surface temperature grid data in the range space of the research area based on the thermal infrared remote sensing data.
In the embodiment, a building island and a Jinmen island which have the same area scale, geographic position and climate condition but obvious urbanization level difference are taken as research areas, landsat series remote sensing image data of clear and cloudless 7-period within the range of the research areas from 2000 to 2020 is obtained, the time resolution is 3 or 4 years, and the spatial resolution is 30m. Based on Landsat TM6/TIRS10 data, the surface temperatures of two residential islands are obtained by single window algorithm inversion, as shown in FIG. 2.
S2: and calculating the average value and the standard deviation of the surface temperature in the research area range, and identifying the high-temperature area and the secondary high-temperature area of the research area according to the surface temperature interval range set by the average value and the standard deviation.
In this example, the average and standard deviation of the surface temperature of each of the two islands were calculated, and the surface temperatures of the xiamen island and jinmen island were divided into sections by the section division rule shown in table 1, and the high temperature zone and the sub-high temperature zone were identified, as shown in fig. 3.
TABLE 1
Wherein, T s Represents the surface temperature, μ represents the mean of the surface temperature, and std represents the standard deviation of the surface temperature.
S3: and calculating the surface heat amplitude index and the surface heat intensity index of the research area based on the areas of the identified high-temperature area and the second high-temperature area.
In this embodiment, the surface heat amplitude index includes a comprehensive surface heat amplitude index corresponding to the whole of the research area and a spatial heat amplitude index corresponding to each grid pixel in the research area, and the surface heat intensity index includes a comprehensive surface heat intensity index corresponding to the whole of the research area and a spatial heat intensity index corresponding to each grid pixel in the research area.
The calculation formula of the comprehensive earth surface heat amplitude index is as follows:
the calculation formula of the comprehensive earth surface heat intensity index is as follows:
wherein, HAI (0)<HAI<1) For comprehensive surface heat amplitude index, HII (0) of the entire study area<HII<1) For the whole grindingComprehensive surface heat intensity index, S, of the research area htz And S shtz Respectively showing the areas of the high-temperature area and the sub-high-temperature area of the research area, S showing the total area of the research area and the heat amplitude index HAI of the space earth surface i (0≤HAI i ≦ 1) represents the surface heat amplitude index of the ith pixel element, and the spatial surface heat intensity index HII i (0≤HII i ≦ 1) represents the surface heat intensity index of the ith pixel,
representing the area of the high temperature region contained by the window centered on the picture element i,
denotes the area of the sub-high temperature region contained in the window centered on the pixel i, S i Representing the total area of the window centred on the picture element i. Therefore, the surface heat intensity index is the area ratio of the high-temperature area of the area, and the difference of the thermal property and the property of the research area is emphasized; the surface heat amplitude index is the area ratio of the regional hot areas (high temperature area and secondary high temperature area), and emphasizes the difference of the thermal influence space range of the research area.
The results of the calculation of the surface heat amplitude index and the surface heat strength index in 2000-2020 for the two populated islands studied in this example are shown in table 2 below. The result shows that in the time sequence evolution, the fluctuation of the earth surface heat amplitude of the mansion island is reduced, but the fluctuation of the earth surface heat intensity is increased; in contrast, the surface heat amplitude of the gold door island increases in the wave, but the surface heat intensity decreases in the wave. In spatial comparison, the annual average ground surface heat amplitude of the mansion island is slightly larger than that of the Jinmen island, but the annual average ground surface heat intensity of the mansion island is smaller than that of the Jinmen island.
TABLE 2
Further, the window integral calculation value is assigned to a window center pixel by using a moving window method, and the surface heat intensity and the surface heat amplitude in the research area range are realized by scanning and calculating pixel by pixelThe degree spatialization is respectively defined as a space surface heat amplitude index and a space surface heat intensity index. The calculation formula of the spatial surface heat amplitude index is as follows:
the calculation formula of the spatial surface heat intensity index is as follows:
wherein the spatial surface heat amplitude index HAI i (0≤HAI i ≦ 1) represents the surface heat amplitude index of the ith pixel element, and the spatial surface heat intensity index HII i (0≤HII i ≦ 1) represents the surface heat intensity index of the ith pixel,
indicates the area of the high-temperature region contained in the window with the pixel i as the center,
denotes the area of the sub-high temperature region contained in the window centered on the pixel i, S i Representing the total area of the window centred on the picture element i.
The mobile window is generally a square grid, the size of the mobile window can directly influence the specific calculation value of the central pixel, and the numerical distribution of the surface thermal effect space distribution result obtained by the calculation of the undersized mobile window is too discrete and lacks of attribute continuity; the spatial granularity of the numerical value of the thermal effect spatial distribution result obtained by calculating the overlarge moving window is larger, and spatial details are easy to lose. In the practical application process, in order to ensure the continuity of the obtained numerical attributes and fully reflect the spatial distribution details of the thermal index, the size of the moving window needs to be determined by simultaneously considering the spatial resolution of the surface temperature grid data and the practical analysis requirements. For the two populated islands studied in this example, the spatial heat amplitude index and the spatial heat intensity index were calculated with a grid size of 210 × 210m, and the results of the spatial distribution of the surface heat amplitude and the surface heat intensity of the mansion islands and the golden islands are shown in fig. 4 and 5.
S4: and carrying out weighted summation on the earth surface heat amplitude index and the earth surface heat intensity index of the research area to obtain the earth surface heat effect index of the research area.
The surface thermal effect index in this embodiment includes a comprehensive surface thermal effect index corresponding to the entire study area and a spatial surface thermal effect index corresponding to each grid pixel in the study area. The calculation formula of the comprehensive earth surface heat effect index HEI (0-HEI-1) is as follows: HEI = α HAI + β HII. Wherein α and β are both weight coefficients and α + β =1. Therefore, the comprehensive earth surface heat effect index integrates the meanings of the area earth surface heat amplitude and the earth surface heat intensity, not only considers the property difference of the thermal property, but also considers the difference of the thermal influence space range, and the weight coefficient can be adjusted according to the actual research target and the demand and the contribution of the two to generate the appropriate targeted heat effect index. In particular, where the earth surface heat amplitude and the earth surface heat intensity are considered to contribute equally to the earth surface thermal effect, α = β =0.5 is taken.
For the two populated islands studied in this example, the earth surface heat magnitude and the earth surface heat intensity are considered to contribute the same to the earth surface heat effect, and a = β =0.5 is taken, and the comprehensive earth surface heat effect index is calculated as shown in table 3 below. The results show that in 2000-2020, the comprehensive earth surface heat effect of the two islands is greatly increased in fluctuation, and the annual average comprehensive earth surface heat effect of the Jinmen island is slightly stronger than that of the Xiamen island.
TABLE 3
| Year of year | 2000 | 2003 | 2006 | 2010 | 2014 | 2017 | 2020 | Mean over many years |
| Mansion gate island | 22.48 | 22.02 | 23.75 | 22.12 | 21.33 | 22.74 | 22.79 | 22.46 |
| Gold door island | 21.75 | 24.41 | 22.36 | 23.00 | 23.16 | 22.59 | 22.95 | 22.89 |
Calculating a spatial earth surface heat effect index by using the spatial earth surface heat amplitude index and the spatial earth surface heat intensity index, and realizing spatialization of the earth surface heat effect, wherein the calculation formula is as follows: HEI i ＝αHAI i +βHII i . Wherein, HEI i (0≤HEI i Less than or equal to 1) is a space thermal effect index, namely the pixel element iSurface thermal effect index of (a). The results of the spatial distribution of the surface thermal effects of the two islands are calculated based on the results of the spatialization of the surface thermal amplitude and the surface thermal intensity of the two islands are shown in fig. 6. It can be seen that under the common influence of the ground surface heat intensity and the ground surface heat effect spatial distribution, the mansion gate island strong ground surface heat effect area is mainly and intensively distributed on the north and south east coasts of the island, and the gold gate island strong ground surface heat effect area is mainly and intensively distributed on the north and south bays.
In summary, the urban surface thermal effect space-time measurement method provided in this embodiment quantitatively evaluates the range (surface thermal amplitude) and the degree (surface thermal intensity) of the regional surface thermal effect by using standardized indexes in combination with specific surface temperature partition patterns and research requirements, and further generates a surface thermal effect index that can be used for longitudinal time dynamic evolution research and transverse regional space contrast analysis, thereby sufficiently characterizing the regional surface thermal environment condition. Particularly, a spatial analysis technology is further utilized, the corresponding thermodynamic index is applied to pixel-by-pixel calculation of the ground temperature partition grid map layer, and a spatially continuous thermodynamic index map layer is generated, so that the spatial heterogeneity characteristics of the urban surface heat environment are reflected more comprehensively.
The specific embodiments provided by the invention show that when the research area relates to a plurality of cities (or other basic analysis units) and long-time sequence analysis, the method for measuring the urban terrestrial heat effect space-time measure provided by the embodiment can reflect the difference of thermal environment characteristics of each landscape type in the same urban unit in different periods, and can also compare the influence of an urban development process on the regional terrestrial heat effect among different urban units in the same period.
Example two:
the invention also provides terminal equipment for measuring the urban ground surface thermal effect space-time, which comprises a memory, a processor and a computer program stored in the memory and capable of running on the processor, wherein the processor executes the computer program to realize the steps of the method embodiment of the first embodiment of the invention.
Further, as an executable scheme, the terminal device for measuring the time-space degree of the urban earth surface thermal effect may be a desktop computer, a notebook computer, a palm computer, a cloud server and other computing devices. The terminal device for urban terrestrial heat effect space-time measurement can comprise, but is not limited to, a processor and a memory. It can be understood by those skilled in the art that the above-mentioned structure of the terminal device for urban terrestrial heat effect space-time measurement is only an example of the terminal device for urban terrestrial heat effect space-time measurement, and does not constitute a limitation on the terminal device for urban terrestrial heat effect space-time measurement, and may include more or less components than the above, or combine some components, or different components, for example, the terminal device for urban terrestrial heat effect space-time measurement may further include an input-output device, a network access device, a bus, and the like, which is not limited in this embodiment of the present invention.
Further, as an executable solution, the Processor may be a Central Processing Unit (CPU), other general purpose Processor, a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), a Field Programmable Gate Array (FPGA) or other Programmable logic device, a discrete Gate or transistor logic device, a discrete hardware component, and the like. The general processor can be a microprocessor or the processor can be any conventional processor and the like, the processor is a control center of the urban terrestrial heat effect space-time measurement terminal equipment, and various interfaces and lines are utilized to connect various parts of the whole urban terrestrial heat effect space-time measurement terminal equipment.
The memory can be used for storing the computer programs and/or modules, and the processor can realize various functions of the terminal equipment for measuring the urban surface thermal effect space-time by running or executing the computer programs and/or modules stored in the memory and calling the data stored in the memory. The memory can mainly comprise a program storage area and a data storage area, wherein the program storage area can store an operating system and an application program required by at least one function; the storage data area may store data created according to the use of the mobile phone, and the like. In addition, the memory may include high speed random access memory, and may also include non-volatile memory, such as a hard disk, a memory, a plug-in hard disk, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash memory Card (Flash Card), at least one magnetic disk storage device, a Flash memory device, or other volatile solid state storage device.
The present invention also provides a computer-readable storage medium, which stores a computer program, which, when executed by a processor, implements the steps of the above-mentioned method of an embodiment of the present invention.
The integrated module/unit of the terminal device for urban surface thermal effect space-time measurement can be stored in a computer readable storage medium if the module/unit is realized in the form of a software functional unit and is sold or used as an independent product. Based on such understanding, all or part of the flow of the method according to the embodiments of the present invention may also be implemented by a computer program, which may be stored in a computer-readable storage medium, and when the computer program is executed by a processor, the steps of the method embodiments may be implemented. Wherein the computer program comprises computer program code, which may be in the form of source code, object code, an executable file or some intermediate form, etc. The computer-readable medium may include: any entity or device capable of carrying the computer program code, recording medium, usb disk, removable hard disk, magnetic disk, optical disk, computer Memory, read-Only Memory (ROM), random Access Memory (RAM), software distribution medium, and the like.
While the invention has been particularly shown and described with reference to a preferred embodiment, it will be understood by those skilled in the art that various changes in form and detail may be made therein without departing from the spirit and scope of the invention as defined by the appended claims.
Claims (7)
1. A space-time measuring method for urban ground surface thermal effect is characterized by comprising the following steps:
s1: based on the thermal infrared remote sensing data, performing inversion to obtain surface temperature grid data which are continuous in the range space of the research area;
s2: calculating the average value and the standard deviation of the surface temperature in the research area range, and identifying a high-temperature area and a secondary high-temperature area of the research area according to the surface temperature interval range set by the average value and the standard deviation;
s3: calculating the surface heat amplitude index and the surface heat intensity index of the research area based on the areas of the identified high-temperature area and the second high-temperature area;
s4: and carrying out weighted summation on the earth surface heat amplitude index and the earth surface heat intensity index of the research area to obtain the earth surface heat effect index of the research area.
2. The method for spatiotemporal measurement of urban terrestrial thermoeffects according to claim 1, characterized in that: in the remote sensing image data adopted in the step S1, the research area range of the remote sensing image data is free from cloud shielding.
3. The method for spatiotemporal measurement of urban terrestrial thermoeffects according to claim 1, characterized in that: in step S2, the earth surface temperature interval range corresponding to the high temperature zone is: t is s Mu + std; the earth surface temperature interval range corresponding to the secondary high temperature area is as follows: mu +0.5std < T s Less than or equal to mu + std; wherein, T s Represents the surface temperature, μ represents the mean of the surface temperature, and std represents the standard deviation of the surface temperature.
4. The method for spatiotemporal measurement of urban terrestrial thermoeffects according to claim 1, characterized in that: the surface heat amplitude index in the step S3 comprises a comprehensive surface heat amplitude index corresponding to the whole research area and a space heat amplitude index corresponding to each grid pixel in the research area, and the surface heat intensity index comprises a comprehensive surface heat intensity index corresponding to the whole research area and a space heat intensity index corresponding to each grid pixel in the research area; in the calculation of the spatial heat amplitude index and the spatial heat intensity index, pixel-by-pixel scanning calculation is carried out by a moving window method, and the integral calculation value of the square window is assigned to a window center pixel;
the calculation formula of the comprehensive earth surface heat amplitude index is as follows:
the calculation formula of the comprehensive surface heat intensity index is as follows:
the calculation formula of the spatial surface heat amplitude index is as follows:
the calculation formula of the spatial surface heat intensity index is as follows:
wherein HAI represents the comprehensive surface heat amplitude index of the whole research area, HII represents the comprehensive surface heat intensity index of the whole research area, and S htz And S shtz Respectively representing the areas of a high-temperature area and a second high-temperature area of the research area, S representing the total area of the research area and the spatial earth surface heat amplitude index HAI i The surface heat amplitude index and the spatial surface heat intensity index HII of the ith pixel element are represented i Represents the surface heat intensity index of the ith pixel element,
indicates the area of the high-temperature region contained in the window with the pixel i as the center,
denotes the area of the sub-high temperature region contained in the window centered on the pixel i, S i Representing the total area of the window centred on the pixel element i.
5. The method for spatiotemporal measurement of urban terrestrial thermoeffects according to claim 4, characterized in that: the surface thermal effect index in the step S4 comprises a comprehensive surface thermal effect index corresponding to the whole research area and a space surface thermal effect index corresponding to each grid pixel in the research area;
the calculation formula of the comprehensive earth surface heat effect index is as follows: HEI = α HAI + β HII;
the calculation formula of the spatial earth surface heat effect index is as follows: HEI i ＝αHAI i +βHII i ；
Wherein HEI represents the comprehensive surface heat effect index of the whole research area, and the space surface heat effect index HEI i The surface thermal effect index of the ith pixel element is shown, alpha and beta both represent weight coefficients and alpha + beta =1.
6. A terminal device for urban earth surface thermal effect space-time measurement is characterized in that: comprising a processor, a memory and a computer program stored in said memory and running on said processor, said processor implementing the steps of the method according to any one of claims 1 to 5 when executing said computer program.
7. A computer-readable storage medium storing a computer program, characterized in that: the computer program realizing the steps of the method as claimed in any one of claims 1 to 5 when executed by a processor.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211259580.7A CN115628814A (en) | 2022-10-14 | 2022-10-14 | Urban surface thermal effect space-time measurement method, terminal device and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211259580.7A CN115628814A (en) | 2022-10-14 | 2022-10-14 | Urban surface thermal effect space-time measurement method, terminal device and storage medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115628814A true CN115628814A (en) | 2023-01-20 |
Family
ID=84904718
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211259580.7A Pending CN115628814A (en) | 2022-10-14 | 2022-10-14 | Urban surface thermal effect space-time measurement method, terminal device and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115628814A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116362567A (en) * | 2023-04-07 | 2023-06-30 | 生态环境部卫星环境应用中心 | Urban industrial heat source enterprise emergency response remote sensing evaluation method |
| CN116842343A (en) * | 2023-07-04 | 2023-10-03 | 南京林业大学 | Method for quantifying influence of urban forest on temperature based on satellite observation and space conversion |
| CN117670873A (en) * | 2024-01-30 | 2024-03-08 | 广东省科学院广州地理研究所 | Thermal comfort index detection method, thermal comfort index detection device, computer equipment and storage medium |
-
2022
- 2022-10-14 CN CN202211259580.7A patent/CN115628814A/en active Pending
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116362567A (en) * | 2023-04-07 | 2023-06-30 | 生态环境部卫星环境应用中心 | Urban industrial heat source enterprise emergency response remote sensing evaluation method |
| CN116362567B (en) * | 2023-04-07 | 2024-01-12 | 生态环境部卫星环境应用中心 | Urban industrial heat source enterprise emergency response remote sensing evaluation method |
| CN116842343A (en) * | 2023-07-04 | 2023-10-03 | 南京林业大学 | Method for quantifying influence of urban forest on temperature based on satellite observation and space conversion |
| CN117670873A (en) * | 2024-01-30 | 2024-03-08 | 广东省科学院广州地理研究所 | Thermal comfort index detection method, thermal comfort index detection device, computer equipment and storage medium |
| CN117670873B (en) * | 2024-01-30 | 2024-05-17 | 广东省科学院广州地理研究所 | Thermal comfort index detection method, thermal comfort index detection device, computer equipment and storage medium |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Sinha et al. | Assessing the spatial sensitivity of a random forest model: Application in gridded population modeling | |
| CN115628814A (en) | Urban surface thermal effect space-time measurement method, terminal device and storage medium | |
| Osipova et al. | Using step‐selection functions to model landscape connectivity for African elephants: Accounting for variability across individuals and seasons | |
| Robertson et al. | A fuzzy classification technique for predicting species’ distributions: applications using invasive alien plants and indigenous insects | |
| Fan et al. | Measuring the spatial arrangement of urban vegetation and its impacts on seasonal surface temperatures | |
| Parolo et al. | Toward improved species niche modelling: Arnica montana in the Alps as a case study | |
| Anderson | A framework for using niche models to estimate impacts of climate change on species distributions | |
| Mennis | Generating surface models of population using dasymetric mapping | |
| Rajasekar et al. | Spatio‐temporal modelling and analysis of urban heat islands by using Landsat TM and ETM+ imagery | |
| Wong et al. | Spatial variability of frontal area index and its relationship with urban heat island intensity | |
| Abdullah et al. | Investigating the impact of land use/land cover change on present and future land surface temperature (LST) of Chittagong, Bangladesh | |
| Surabuddin Mondal et al. | Modeling of spatio-temporal dynamics of land use and land cover in a part of Brahmaputra River basin using Geoinformatic techniques | |
| Thies et al. | Projecting land-use and land-cover changes in a tropical mountain forest of Southern Ecuador | |
| Kowe et al. | Effect of landscape pattern and spatial configuration of vegetation patches on urban warming and cooling in Harare metropolitan city, Zimbabwe | |
| Monahan et al. | Niche tracking and rapid establishment of distributional equilibrium in the house sparrow show potential responsiveness of species to climate change | |
| Biard et al. | Automated detection of weather fronts using a deep learning neural network | |
| Wisz et al. | Using potential distributions to explore determinants of Western Palaearctic migratory songbird species richness in sub‐Saharan Africa | |
| Pal et al. | Hurst exponent approach through rescaled range analysis to study the time series of summer monsoon rainfall over northeast India | |
| Su et al. | Geospatial assessment of agroecosystem health: development of an integrated index based on catastrophe theory | |
| Zhu et al. | Investigating the uncertainty and transferability of parameters in SWAT model under climate change | |
| Pan et al. | Understanding spatiotemporal evolution of the surface urban heat island in the Bangkok metropolitan region from 2000 to 2020 using enhanced land surface temperature | |
| Dost et al. | Quantification of water resource sustainability in response to drought risk assessment for Afghanistan river basins | |
| Desjardins et al. | A space–time parallel framework for fine-scale visualization of pollen levels across the Eastern United States | |
| Morrison et al. | Spatial scale affects novel and disappeared climate change projections in Alaska | |
| CN117078096A (en) | Urban activity comprehensive measurement method, system, equipment and medium |
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 |