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 PDF

Info

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
area
surface heat
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
Application number
CN202211259580.7A
Other languages
Chinese (zh)
Inventor
吝涛
张浚茂
林美霞
姚霞
陈媛
张喻魁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Urban Environment of CAS
Original Assignee
Institute of Urban Environment of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Institute of Urban Environment of CAS filed Critical Institute of Urban Environment of CAS
Priority to CN202211259580.7A priority Critical patent/CN115628814A/en
Publication of CN115628814A publication Critical patent/CN115628814A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00Radiation pyrometry, e.g. infrared or optical thermometry
    • G01J5/007Radiation pyrometry, e.g. infrared or optical thermometry for earth observation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/10Terrestrial scenes
    • G06V20/13Satellite images
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A30/00Adapting or protecting infrastructure or their operation
    • Y02A30/60Planning 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

本发明涉及一种城市地表热效应时空测度方法、终端设备及存储介质,该方法中包括:S1:基于热红外遥感数据,反演得到研究区范围空间连续的地表温度栅格数据;S2:计算研究区范围地表温度的平均值和标准差,并根据平均值和标准差所设定的地表温度区间范围,识别研究区的高温区和次高温区;S3:基于识别到的高温区和次高温区的面积,计算研究区的地表热幅度指数和地表热强度指数;S4:对研究区的地表热幅度指数与地表热强度指数进行加权求和,得到研究区的地表热效应指数。本发明以标准化指数对区域地表热影响的范围(热幅度)和程度(热强度)进行了定量化评估,能够生成空间连续的热力指数图层,从而更全面地反映城市地表热环境的空间异质性特征。

Figure 202211259580

The present invention relates to a time-space measurement method of urban surface thermal effect, terminal equipment and storage medium. The method includes: S1: Based on thermal infrared remote sensing data, inversion obtains spatially continuous surface temperature grid data in the research area; S2: Computational research The average and standard deviation of the surface temperature in the area, and according to the range of the surface temperature set by the average and standard deviation, identify the high temperature area and sub-high temperature area in the study area; S3: Based on the identified high temperature area and sub-high temperature area Calculate the surface heat magnitude index and surface heat intensity index of the study area; S4: The weighted sum of the surface heat magnitude index and the surface heat intensity index of the study area is obtained to obtain the surface heat effect index of the study area. The present invention quantitatively evaluates the scope (heat range) and degree (heat intensity) of regional surface heat influence by standardized index, and can generate a spatially continuous thermal index layer, thereby more comprehensively reflecting the spatial variation of the urban surface heat environment. qualitative features.

Figure 202211259580

Description

一种城市地表热效应时空测度方法、终端设备及存储介质A spatio-temporal measurement method, terminal equipment, and storage medium for urban surface thermal effects

技术领域technical field

本发明涉及城市生态学中遥感地理信息技术应用领域,尤其涉及一种城市地表热效应时空测度方法、终端设备及存储介质。The invention relates to the application field of remote sensing geographic information technology in urban ecology, in particular to a method for measuring urban surface thermal effects in time and space, a terminal device and a storage medium.

背景技术Background technique

城市地表热环境是能够对人体的冷暖感受、健康水平和居民生存发展等产生重要影响的与热有关的物理环境,通常以地表温度表征。随着城市化进程的不断推进以及城市人口的不断增长,城市内部大量的自然地表被人工地表所取代,下垫面的更新,特别是不透水面的剧增,通过改变地表热辐射传输的平衡提升了城市局部或整体的地表温度,加剧城市地表热效应。城市地表热效应的加剧是城市地表热环境恶化的主要原因和重要表现,其中最为典型的热岛效应与城市居民健康、社会经济发展、自然生态环境保护等息息相关。The urban surface thermal environment is a heat-related physical environment that can have an important impact on the human body's cold and warm feelings, health levels, and the survival and development of residents. It is usually characterized by surface temperature. With the continuous advancement of the urbanization process and the continuous growth of the urban population, a large number of natural surfaces in the city are replaced by artificial surfaces. It increases the local or overall surface temperature of the city and intensifies the thermal effect of the urban surface. The intensification of urban surface thermal effect is the main reason and important manifestation of the deterioration of urban surface thermal environment. The most typical heat island effect is closely related to the health of urban residents, social and economic development, and natural ecological environment protection.

城市地表热效应的量化评估是城市地表热环境研究的重要基础,已有学者对城市地表热效应的定量化测度进行了一定的尝试,如基于城市下垫面热力性质差异建立的三维动态热效应数值模型,或者利用平均地表温度差值计算的不同地表类型的热效应,但目前的地表热效应相关量化指数仍存在以下不足之处:时间维度上,由于遥感数据获取热辐射信息的“瞬时性”特征以及“高时、空间分辨率难以统一”的问题,直接基于地表温度计算得到的热效应指数存在极大的不确定性,从而缺乏时序代表性,这使得城市地表热效应的长时间动态演变无法得到有效分析;尽管有部分指数在计算过程中利用平均地表温度消除了时间序列上的不确定性,但大多是基于土地利用或覆被的分类统计,因此只能得到相应地类整体的热效应值,从而在空间维度上,既无法充分识别城市地表热效应的空间异质性,导致同一城市内不同区域之间地表热效应对比未能得到有效的度量和分析,又因不同城市的地表结构差异而不利于进行城市间的对比分析,这为充分理解城市热环境特征及其影响因素的内外差异带来了一定的困难。Quantitative evaluation of urban surface thermal effects is an important basis for the study of urban surface thermal environment. Scholars have made some attempts to quantify urban surface thermal effects, such as the 3D dynamic thermal effect numerical model established based on the thermal properties of urban underlying surfaces. Or use the average surface temperature difference to calculate the thermal effects of different surface types, but the current quantitative index of surface thermal effects still has the following shortcomings: In the time dimension, due to the "instantaneous" characteristics of thermal radiation information obtained from remote sensing data and the "high The thermal effect index calculated directly based on the surface temperature has great uncertainty, so it lacks time-series representativeness, which makes it impossible to effectively analyze the long-term dynamic evolution of the urban surface thermal effect; although Some indexes use the average surface temperature in the calculation process to eliminate the uncertainty in the time series, but most of them are based on the classification statistics of land use or cover, so only the thermal effect value of the corresponding land type can be obtained, so that in the spatial dimension In fact, the spatial heterogeneity of urban surface thermal effects cannot be fully identified, resulting in the ineffective measurement and analysis of the comparison of surface thermal effects between different regions in the same city, and the differences in the surface structure of different cities are not conducive to inter-city comparisons. Comparative analysis has brought certain difficulties to fully understand the characteristics of urban thermal environment and the internal and external differences of its influencing factors.

总体而言,尚未见有充分反映地表热效应含义的通用指数,且相关指数的代表性和使用有效性有待提高。因此,亟需发展一种不受时空间量纲差异影响且高效通用的城市地表热效应多尺度时空测度方法。In general, there is no general index that fully reflects the meaning of surface thermal effect, and the representativeness and effectiveness of related indices need to be improved. Therefore, there is an urgent need to develop a multi-scale spatio-temporal measurement method for urban surface thermal effects that is not affected by the difference in spatio-temporal dimension and is efficient and general.

发明内容Contents of the invention

为了解决上述问题,本发明提出了一种城市地表热效应时空测度方法、终端设备及存储介质。In order to solve the above problems, the present invention proposes a spatio-temporal measurement method, terminal equipment and storage medium for urban surface thermal effects.

具体方案如下:The specific plan is as follows:

一种城市地表热效应时空测度方法,包括以下步骤:A spatio-temporal measurement method for urban surface thermal effects, comprising the following steps:

S1:基于热红外遥感数据,反演得到研究区范围空间连续的地表温度栅格数据;S1: Based on the thermal infrared remote sensing data, the spatially continuous surface temperature grid data of the study area is retrieved;

S2:计算研究区范围地表温度的平均值和标准差,并根据平均值和标准差所设定的地表温度区间范围,识别研究区的高温区和次高温区;S2: Calculate the average value and standard deviation of the surface temperature in the study area, and identify the high temperature area and sub-high temperature area in the study area according to the range of surface temperature set by the average value and standard deviation;

S3:基于识别到的高温区和次高温区的面积,计算研究区的地表热幅度指数和地表热强度指数;S3: Calculate the surface heat amplitude index and surface heat intensity index of the study area based on the area of the identified high temperature area and sub-high temperature area;

S4:对研究区的地表热幅度指数与地表热强度指数进行加权求和,得到研究区的地表热效应指数。S4: Weighted and summed the surface heat amplitude index and surface heat intensity index in the study area to obtain the surface heat effect index in the study area.

进一步的,步骤S1中所采用的遥感影像数据,其研究区范围应无云遮挡。Further, the remote sensing image data used in step S1 should have no cloud cover in the research area.

进一步的,步骤S2中,高温区对应的地表温度区间范围为:Ts>μ+std;次高温区对应的地表温度区间范围为:μ+0.5std<Ts≤μ+std;其中,Ts表示地表温度,μ表示地表温度的平均值,std表示地表温度的标准差。Further, in step S2, the range of the surface temperature range corresponding to the high temperature zone is: T s >μ+std; the range of the surface temperature range corresponding to the sub-high temperature zone is: μ+0.5std<T s ≤μ+std; where, T s represents the surface temperature, μ represents the average value of the surface temperature, and std represents the standard deviation of the surface temperature.

进一步的,步骤S3中的地表热幅度指数包括对应研究区整体的综合地表热幅度指数和对应研究区中每一个栅格像元的空间热幅度指数,地表热强度指数包括对应研究区整体的综合地表热强度指数和对应研究区中每一个栅格像元的空间热强度指数;在空间热幅度指数和空间热强度指数的计算中,通过移动窗口法逐像元扫描计算,将正方形窗口的整体计算值赋值给窗口中心像元;Further, the surface heat magnitude index in step S3 includes the comprehensive surface heat magnitude index corresponding to the whole research area and the spatial heat magnitude index corresponding to each grid cell in the research area, and the surface heat intensity index includes the comprehensive comprehensive index corresponding to the whole research area The surface thermal intensity index and the spatial thermal intensity index corresponding to each grid pixel in the study area; in the calculation of the spatial thermal amplitude index and the spatial thermal intensity index, the moving window method is used to scan and calculate pixel by pixel, and the whole square window The calculated value is assigned to the center pixel of the window;

综合地表热幅度指数的计算公式为:

Figure BDA0003890942060000031
The formula for calculating the comprehensive surface heat range index is:
Figure BDA0003890942060000031

综合地表热强度指数的计算公式为:

Figure BDA0003890942060000032
The formula for calculating the comprehensive surface heat intensity index is:
Figure BDA0003890942060000032

空间地表热幅度指数的计算公式为:

Figure BDA0003890942060000033
The calculation formula of the spatial surface heat amplitude index is:
Figure BDA0003890942060000033

空间地表热强度指数的计算公式为:

Figure BDA0003890942060000034
The calculation formula of space surface heat intensity index is:
Figure BDA0003890942060000034

其中,HAI表示整个研究区的综合地表热幅度指数,HII表示整个研究区的综合地表热强度指数,Shtz和Sshtz分别表示研究区的高温区和次高温区的面积,S表示研究区的总面积,空间地表热幅度指数HAIi表示第i个像元的地表热幅度指数,空间地表热强度指数HIIi表示第i个像元的地表热强度指数,

Figure BDA0003890942060000035
表示以像元i为中心的窗口包含的高温区的面积,
Figure BDA0003890942060000036
表示以像元i为中心的窗口包含的次高温区的面积,Si表示以像元i为中心的窗口的总面积。Among them, HAI represents the comprehensive surface heat amplitude index of the entire study area, HII represents the comprehensive surface heat intensity index of the entire study area, Shtz and S shtz represent the areas of high temperature and sub-high temperature regions in the study area, respectively, and S represents the area of the study area The total area, the spatial surface heat amplitude index HAI i represents the surface heat amplitude index of the i-th pixel, and the spatial surface heat intensity index HII i represents the surface heat intensity index of the i-th pixel,
Figure BDA0003890942060000035
Indicates the area of the high-temperature region contained in the window centered on pixel i,
Figure BDA0003890942060000036
Indicates the area of the sub-high temperature zone contained in the window centered on pixel i, and S i indicates the total area of the window centered on pixel i.

进一步的,步骤S4中的地表热效应指数包括对应研究区整体的综合地表热效应指数和对应研究区中每一个栅格像元的空间地表热效应指数;Further, the surface thermal effect index in step S4 includes the comprehensive surface thermal effect index corresponding to the whole research area and the spatial surface thermal effect index corresponding to each grid cell in the research area;

综合地表热效应指数的计算公式为:HEI=αHAI+βHII;The calculation formula of comprehensive surface thermal effect index is: HEI=αHAI+βHII;

空间地表热效应指数的计算公式为:HEIi=αHAIi+βHIIiThe calculation formula of space surface thermal effect index is: HEI i = αHAI i + βHII i ;

其中,HEI表示整个研究区的综合地表热效应指数,空间地表热效应指数HEIi表示第i个像元的地表热效应指数,α和β均表示权重系数且α+β=1。Among them, HEI represents the comprehensive surface heating effect index of the entire study area, the spatial surface heating effect index HEI i represents the surface heating effect index of the i-th pixel, α and β both represent weight coefficients and α+β=1.

一种城市地表热效应时空测度终端设备,包括处理器、存储器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现本发明实施例上述的方法的步骤。A terminal device for space-time measurement of urban surface thermal effects, including a processor, a memory, and a computer program stored in the memory and operable on the processor, the processor implements the embodiment of the present invention when executing the computer program steps of the method described above.

一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时实现本发明实施例上述的方法的步骤。A computer-readable storage medium. The computer-readable storage medium stores a computer program. When the computer program is executed by a processor, the steps of the above-mentioned method in the embodiment of the present invention are implemented.

本发明采用如上技术方案,弥补了现有相关指数对区域地表热效应空间分布特征刻画的不足,同时解决区域地表热效应时空间对比分析困难的问题,实现对城市地表热效应的科学定量有效测度,直观地揭示城市地表热环境的时间演变及空间分布特征。The present invention adopts the above technical scheme, which makes up for the deficiency of the existing correlation index in describing the spatial distribution characteristics of the regional surface thermal effect, and at the same time solves the problem of difficult time-space comparative analysis of the regional surface thermal effect, and realizes the scientific, quantitative and effective measurement of the urban surface thermal effect, intuitively Reveal the temporal evolution and spatial distribution characteristics of urban surface thermal environment.

附图说明Description of drawings

图1所示为本发明实施例一的流程图。FIG. 1 is a flowchart of Embodiment 1 of the present invention.

图2所示为本发明实施例一的地表温度反演结果图。Fig. 2 is a graph showing the inversion results of the surface temperature in Embodiment 1 of the present invention.

图3所示为本发明实施例一的地表温度分区结果图。Fig. 3 is a diagram showing the results of the land surface temperature divisions in Embodiment 1 of the present invention.

图4所示为本发明实施例一的地表热幅度空间分布结果图。FIG. 4 is a diagram showing the result of spatial distribution of surface heat amplitude in Embodiment 1 of the present invention.

图5所示为本发明实施例一的地表热强度空间分布结果图。FIG. 5 is a diagram showing the spatial distribution results of surface thermal intensity in Embodiment 1 of the present invention.

图6所示为本发明实施例一的地表热效应空间分布结果图。FIG. 6 is a diagram showing the spatial distribution results of surface thermal effects in Embodiment 1 of the present invention.

具体实施方式Detailed ways

为进一步说明各实施例,本发明提供有附图。这些附图为本发明揭露内容的一部分,其主要用以说明实施例,并可配合说明书的相关描述来解释实施例的运作原理。配合参考这些内容,本领域普通技术人员应能理解其他可能的实施方式以及本发明的优点。To further illustrate the various embodiments, the present invention is provided with accompanying drawings. These drawings are a part of the disclosure of the present invention, which are mainly used to illustrate the embodiments, and can be combined with related descriptions in the specification to explain the operating principles of the embodiments. With reference to these contents, those skilled in the art should understand other possible implementations and advantages of the present invention.

现结合附图和具体实施方式对本发明作进一步说明。The present invention will be further described in conjunction with the accompanying drawings and specific embodiments.

实施例一:Embodiment one:

本发明实施例提供了一种城市地表热效应时空测度方法,如图1所示,所述方法包括以下步骤:An embodiment of the present invention provides a method for measuring urban surface thermal effects in space and time, as shown in Figure 1, the method includes the following steps:

S1:基于热红外遥感数据,反演得到研究区范围空间连续的地表温度栅格数据。S1: Based on the thermal infrared remote sensing data, the spatially continuous surface temperature grid data of the study area is obtained by inversion.

该实施例中以面积规模、地理位置、气候条件均相同,但城市化水平差异显著的厦门岛和金门岛为研究区,获取了研究区范围2000-2020年间7期晴朗无云的Landsat系列遥感影像数据,时间分辨率为3或4年,空间分辨率为30m。基于Landsat TM6/TIRS10数据,使用单窗算法反演得到两个有居民海岛的地表温度,如图2所示。In this example, Xiamen Island and Kinmen Island, which have the same area size, geographical location, and climate conditions, but with significant differences in urbanization levels, were used as the research area, and Landsat series remote sensing data of seven periods of clear and cloudless periods in the study area from 2000 to 2020 were obtained. Image data with a temporal resolution of 3 or 4 years and a spatial resolution of 30m. Based on the Landsat TM6/TIRS10 data, the surface temperature of two inhabited islands was obtained by using the single-window algorithm, as shown in Figure 2.

S2:计算研究区范围地表温度的平均值和标准差,并根据平均值和标准差所设定的地表温度区间范围,识别研究区的高温区和次高温区。S2: Calculate the average and standard deviation of the surface temperature in the study area, and identify the high temperature area and sub-high temperature area in the study area according to the range of surface temperature set by the average and standard deviation.

该实施例中,在计算两个海岛各自地表温度的平均值与标准差的基础上,通过表1所示的区间划分原则对厦门岛和金门岛的地表温度进行分区,识别高温区和次高温区,如图3所示。In this example, on the basis of calculating the average and standard deviation of the respective surface temperatures of the two islands, the surface temperatures of Xiamen Island and Jinmen Island are divided according to the principle of interval division shown in Table 1, and high temperature areas and sub-high temperature areas are identified area, as shown in Figure 3.

表1Table 1

Figure BDA0003890942060000051
Figure BDA0003890942060000051

Figure BDA0003890942060000061
Figure BDA0003890942060000061

其中,Ts表示地表温度,μ表示地表温度的平均值,std表示地表温度的标准差。Among them, T s represents the surface temperature, μ represents the average value of the surface temperature, and std represents the standard deviation of the surface temperature.

S3:基于识别到的高温区和次高温区的面积,计算研究区的地表热幅度指数和地表热强度指数。S3: Calculate the surface heat amplitude index and surface heat intensity index of the study area based on the areas of the identified high temperature areas and sub-high temperature areas.

该实施例中,地表热幅度指数包括对应研究区整体的综合地表热幅度指数和对应研究区中每一个栅格像元的空间热幅度指数,地表热强度指数包括对应研究区整体的综合地表热强度指数和对应研究区中每一个栅格像元的空间热强度指数。In this embodiment, the surface heat magnitude index includes the comprehensive surface heat magnitude index corresponding to the whole research area and the spatial heat magnitude index corresponding to each grid cell in the research area, and the surface heat intensity index includes the comprehensive surface heat magnitude index corresponding to the whole research area. The intensity index and the spatial heat intensity index of each grid cell in the corresponding study area.

综合地表热幅度指数的计算公式为:

Figure BDA0003890942060000062
The formula for calculating the comprehensive surface heat range index is:
Figure BDA0003890942060000062

综合地表热强度指数的计算公式为:

Figure BDA0003890942060000063
The formula for calculating the comprehensive surface heat intensity index is:
Figure BDA0003890942060000063

其中,HAI(0<HAI<1)为整个研究区的综合地表热幅度指数,HII(0<HII<1)为整个研究区的综合地表热强度指数,Shtz和Sshtz分别表示研究区的高温区和次高温区的面积,S表示研究区的总面积,空间地表热幅度指数HAIi(0≤HAIi≤1)表示第i个像元的地表热幅度指数,空间地表热强度指数HIIi(0≤HIIi≤1)表示第i个像元的地表热强度指数,

Figure BDA0003890942060000064
表示以像元i为中心的窗口包含的高温区的面积,
Figure BDA0003890942060000065
表示以像元i为中心的窗口包含的次高温区的面积,Si表示以像元i为中心的窗口的总面积。可见,地表热强度指数为区域高温区面积占比,强调研究区的热力属性性质差异;地表热幅度指数为区域热区(高温区和次高温区)面积占比,强调研究区的热力影响空间范围差异。Among them, HAI (0<HAI<1) is the comprehensive surface heat amplitude index of the whole study area, HII (0<HII<1) is the comprehensive surface heat intensity index of the whole study area, Shtz and S shtz respectively represent the The area of the high temperature area and the sub-high temperature area, S represents the total area of the study area, the spatial surface heat amplitude index HAI i (0≤HAI i ≤1) represents the surface heat amplitude index of the i-th pixel, and the spatial surface heat intensity index HII i (0≤HII i ≤1) represents the surface heat intensity index of the i-th pixel,
Figure BDA0003890942060000064
Indicates the area of the high-temperature region contained in the window centered on pixel i,
Figure BDA0003890942060000065
Indicates the area of the sub-high temperature zone contained in the window centered on pixel i, and S i indicates the total area of the window centered on pixel i. It can be seen that the surface heat intensity index is the area ratio of the regional high temperature area, emphasizing the differences in the thermal properties of the study area; the surface heat amplitude index is the area proportion of the regional hot area (high temperature area and sub-high temperature area), emphasizing the thermal influence space of the study area range difference.

对于本实施例中所研究的两个有居民海岛,其2000-2020年的地表热幅度指数和地表热强度指数计算结果如下表2所示。结果表明,时序演变上,厦门岛的地表热幅度波动减小,但地表热强度波动增大;与此相反,金门岛的地表热幅度在波动中增大,但地表热强度在波动中减小。空间对比上,厦门岛的多年平均地表热幅度略大于金门岛,但多年平均地表热强度小于金门岛。For the two inhabited islands studied in this example, the calculation results of the surface heat amplitude index and surface heat intensity index from 2000 to 2020 are shown in Table 2 below. The results show that, in terms of time series evolution, the fluctuation of surface heat amplitude in Xiamen Island decreases, but the fluctuation of surface heat intensity increases; on the contrary, the amplitude of surface heat in Jinmen Island increases in fluctuation, but the surface heat intensity decreases in fluctuation . In terms of spatial comparison, the annual average surface heat amplitude of Xiamen Island is slightly larger than that of Kinmen Island, but the annual average surface heat intensity is smaller than that of Kinmen Island.

表2Table 2

Figure BDA0003890942060000071
Figure BDA0003890942060000071

进一步的,使用移动窗口法将窗口整体计算值赋值给窗口中心像元,并通过逐像元扫描计算,实现研究区范围内地表热强度和地表热幅度的空间化,分别定义为空间地表热幅度指数和空间地表热强度指数。空间地表热幅度指数的计算公式为:

Figure BDA0003890942060000072
空间地表热强度指数的计算公式为:
Figure BDA0003890942060000081
其中,空间地表热幅度指数HAIi(0≤HAIi≤1)表示第i个像元的地表热幅度指数,空间地表热强度指数HIIi(0≤HIIi≤1)表示第i个像元的地表热强度指数,
Figure BDA0003890942060000082
表示以像元i为中心的窗口包含的高温区的面积,
Figure BDA0003890942060000083
表示以像元i为中心的窗口包含的次高温区的面积,Si表示以像元i为中心的窗口的总面积。Furthermore, the moving window method is used to assign the overall calculation value of the window to the center pixel of the window, and through pixel-by-pixel scanning calculation, the spatialization of the surface heat intensity and surface heat amplitude within the study area is realized, which are respectively defined as the spatial surface heat amplitude Index and Spatial Surface Heat Intensity Index. The calculation formula of the spatial surface heat amplitude index is:
Figure BDA0003890942060000072
The calculation formula of space surface heat intensity index is:
Figure BDA0003890942060000081
Among them, the spatial surface heat amplitude index HAI i (0≤HAI i ≤1) represents the surface heat amplitude index of the i-th pixel, and the spatial surface heat intensity index HII i (0≤HII i ≤1) represents the i-th pixel The surface heat intensity index of
Figure BDA0003890942060000082
Indicates the area of the high-temperature region contained in the window centered on pixel i,
Figure BDA0003890942060000083
Indicates the area of the sub-high temperature zone contained in the window centered on pixel i, and S i indicates the total area of the window centered on pixel i.

移动窗口一般为正方形格网,其大小会直接影响其中心像元的具体计算值,且过小的移动窗口计算得到的地表热效应空间分布结果,其数值分布过于离散,缺乏属性连续性;而过大的移动窗口计算得到的热效应空间分布结果,其数值的空间粒度较大,容易丢失空间细节。在实际应用过程中,为了保证所得数值属性的连续性,又能充分反映热指数的空间分布细节,需要同时考虑地表温度栅格数据的空间分辨率和实际分析需求来确定移动窗口的大小。对于本实施例中所研究的两个有居民海岛,以210×210m的网格大小计算其空间热幅度指数和空间热强度指数,得到厦门岛和金门岛的地表热幅度和地表热强度空间分布结果如图4和图5所示。The moving window is generally a square grid, and its size will directly affect the specific calculation value of its central pixel, and the spatial distribution of the surface thermal effect calculated by the moving window that is too small, its numerical distribution is too discrete and lacks attribute continuity; The spatial distribution of the thermal effect calculated with a large moving window has a large spatial granularity, and the spatial details are easily lost. In the actual application process, in order to ensure the continuity of the obtained numerical attributes and fully reflect the details of the spatial distribution of the heat index, it is necessary to consider the spatial resolution of the surface temperature raster data and the actual analysis requirements to determine the size of the moving window. For the two inhabited islands studied in this example, the spatial thermal amplitude index and spatial thermal intensity index are calculated with a grid size of 210×210m to obtain the spatial distribution of the surface thermal amplitude and surface thermal intensity of Xiamen Island and Kinmen Island The results are shown in Figure 4 and Figure 5.

S4:对研究区的地表热幅度指数与地表热强度指数进行加权求和,得到研究区的地表热效应指数。S4: Weighted and summed the surface heat amplitude index and surface heat intensity index in the study area to obtain the surface heat effect index in the study area.

该实施例中地表热效应指数包括对应研究区整体的综合地表热效应指数和对应研究区中每一个栅格像元的空间地表热效应指数。综合地表热效应指数HEI(0<HEI<1)的计算公式为:HEI=αHAI+βHII。其中,α和β均为权重系数且α+β=1。可见,综合地表热效应指数综合了区域地表热幅度和地表热强度两方面的含义,既考虑了热力属性性质差异,又考虑了热力影响空间范围差异,并且可按照实际研究目标和需求,根据两者的贡献,调整权重系数以生成合适的针对性热效应指数。特别地,在认为地表热幅度和地表热强度对地表热效应的贡献相同的情况下,取α=β=0.5。In this embodiment, the surface thermal effect index includes a comprehensive surface thermal effect index corresponding to the whole research area and a spatial surface thermal effect index corresponding to each grid cell in the research area. The calculation formula of comprehensive surface heat effect index HEI (0<HEI<1) is: HEI=αHAI+βHII. Wherein, α and β are weight coefficients and α+β=1. It can be seen that the comprehensive surface thermal effect index combines the meanings of the regional surface thermal range and surface thermal intensity, not only considers the difference in thermal properties, but also considers the difference in the spatial range of thermal influence, and can be based on the actual research goals and needs. , the weight coefficients are adjusted to generate a suitable targeted thermal effect index. In particular, α = β = 0.5 under the assumption that the magnitude of the surface heat and the intensity of the surface heat have the same contribution to the surface heat effect.

对于本实施例中所研究的两个有居民海岛,认为地表热幅度和地表热强度对地表热效应的贡献相同,取α=β=0.5,计算得到其综合地表热效应指数如下表3所示。结果表明,2000-2020年,两个海岛的综合地表热效应均在波动中有所加剧,且金门岛的多年平均综合地表热效应略强于厦门岛。For the two inhabited islands studied in this example, it is considered that the surface heat range and surface heat intensity have the same contribution to the surface heat effect, and α = β = 0.5, and the calculated comprehensive surface heat effect index is shown in Table 3 below. The results show that, from 2000 to 2020, the comprehensive surface thermal effect of the two islands fluctuated and intensified, and the multi-year average comprehensive surface thermal effect of Kinmen Island was slightly stronger than that of Xiamen Island.

表3table 3

年份years 20002000 20032003 20062006 20102010 20142014 20172017 20202020 多年平均Multi-year average 厦门岛Xiamen Island 22.4822.48 22.0222.02 23.7523.75 22.1222.12 21.3321.33 22.7422.74 22.7922.79 22.4622.46 金门岛Golden Gate Island 21.7521.75 24.4124.41 22.3622.36 23.0023.00 23.1623.16 22.5922.59 22.9522.95 22.8922.89

使用空间地表热幅度指数和空间地表热强度指数计算空间地表热效应指数,实现地表热效应的空间化,计算公式为:HEIi=αHAIi+βHIIi。其中,HEIi(0≤HEIi≤1)为空间热效应指数,即像元i的地表热效应指数。基于两个海岛的地表热幅度和地表热强度空间化结果,计算得到两者的地表热效应空间分布结果如图6所示。可以看到,受地表热强度和地表热效应空间分布的共同影响,厦门岛强地表热效应区域主要集中分布在海岛北部和东南部沿岸,金门岛强地表热效应区域则主要集中分布在北部沿岸及南部海湾。The space surface heat effect index is calculated by using the space surface heat amplitude index and the space surface heat intensity index to realize the spatialization of the surface heat effect. The calculation formula is: HEI i =αHAI i +βHII i . Among them, HEI i (0≤HEI i ≤1) is the spatial thermal effect index, that is, the surface thermal effect index of pixel i. Based on the spatialization results of the surface heat amplitude and surface heat intensity of the two islands, the spatial distribution results of the surface heat effects of the two islands are shown in Figure 6. It can be seen that due to the joint influence of the surface heat intensity and the spatial distribution of the surface heat effect, the areas with strong surface heat effects in Xiamen Island are mainly concentrated in the northern and southeastern coasts of the island, while the areas with strong surface heat effects in Kinmen Island are mainly distributed in the northern coast and the southern bay. .

综上,本实施例提供的城市地表热效应时空测度方法,通过结合具体的地表温度分区格局和研究需求,以标准化的指数对区域地表热影响的范围(地表热幅度)和程度(地表热强度)进行了定量化评估,进而生成可同时用于纵向时间动态演变探究和横向区域空间对比分析的地表热效应指数,充分表征区域地表热环境状况。特别地,进一步利用空间分析技术,将对应的热力指数应用于地温分区栅格图层的逐像元计算,生成空间连续的热力指数图层,从而更全面地反映城市地表热环境的空间异质性特征。In summary, the spatio-temporal measurement method of urban surface heat effect provided by this example combines the specific surface temperature zoning pattern and research needs, and uses a standardized index to affect the range (surface heat range) and degree (surface heat intensity) of the regional surface heat. Quantitative evaluation was carried out to generate a surface thermal effect index that can be used for both longitudinal temporal dynamic evolution exploration and horizontal regional spatial comparative analysis, fully characterizing the regional surface thermal environment. In particular, the spatial analysis technology is further used to apply the corresponding thermal index to the pixel-by-pixel calculation of the geothermal zonal grid layer to generate a spatially continuous thermal index layer, thereby more comprehensively reflecting the spatial heterogeneity of the urban surface thermal environment sexual characteristics.

本发明提供的具体实施例一同时表明,当研究区涉及多个城市(或其他基本分析单元)以及长时间序列分析时,本实施例提供的城市地表热效应时空测度方法不仅能够反映同一城市单元内部各景观类型在不同时期所呈现出的热环境特征差异,还能比较同一时期的不同城市单元之间,城市化发展进程对区域地表热效应造成的影响。The specific embodiment provided by the present invention shows that at the same time, when the research area involves multiple cities (or other basic analysis units) and long-term series analysis, the urban surface thermal effect time-space measurement method provided by this embodiment can not only reflect the same city unit. The differences in thermal environment characteristics of each landscape type in different periods can also compare the impact of urbanization development process on regional surface thermal effects among different urban units in the same period.

实施例二:Embodiment two:

本发明还提供一种城市地表热效应时空测度终端设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现本发明实施例一的上述方法实施例中的步骤。The present invention also provides a terminal device for space-time measurement of urban surface thermal effects, including a memory, a processor, and a computer program stored in the memory and operable on the processor. When the processor executes the computer program, the Steps in the above method embodiment of Embodiment 1 of the present invention.

进一步地,作为一个可执行方案,所述城市地表热效应时空测度终端设备可以是桌上型计算机、笔记本、掌上电脑及云端服务器等计算设备。所述城市地表热效应时空测度终端设备可包括,但不仅限于,处理器、存储器。本领域技术人员可以理解,上述城市地表热效应时空测度终端设备的组成结构仅仅是城市地表热效应时空测度终端设备的示例,并不构成对城市地表热效应时空测度终端设备的限定,可以包括比上述更多或更少的部件,或者组合某些部件,或者不同的部件,例如所述城市地表热效应时空测度终端设备还可以包括输入输出设备、网络接入设备、总线等,本发明实施例对此不做限定。Furthermore, as an executable solution, the terminal device for measuring urban surface thermal effects in space and time can be computing devices such as desktop computers, notebooks, palmtop computers, and cloud servers. The spatio-temporal measurement terminal device of urban surface thermal effect may include, but not limited to, a processor and a memory. Those skilled in the art can understand that the composition and structure of the above-mentioned space-time measurement terminal equipment of urban surface thermal effect is only an example of the time-space measurement terminal equipment of urban surface thermal effect, and does not constitute a limitation on the time-space measurement terminal equipment of urban surface thermal effect, and may include more than the above or fewer components, or a combination of certain components, or different components, for example, the terminal device for measuring the thermal effect of urban surface space-time may also include input and output devices, network access devices, buses, etc., which are not covered by the embodiment of the present invention limited.

进一步地,作为一个可执行方案,所称处理器可以是中央处理单元(CentralProcessing Unit,CPU),还可以是其他通用处理器、数字信号处理器(Digital SignalProcessor,DSP)、专用集成电路(Application Specific Integrated Circuit,ASIC)、现场可编程门阵列(Field-Programmable Gate Array,FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等,所述处理器是所述城市地表热效应时空测度终端设备的控制中心,利用各种接口和线路连接整个城市地表热效应时空测度终端设备的各个部分。Further, as an executable solution, the so-called processor can be a central processing unit (Central Processing Unit, CPU), and can also be other general-purpose processors, digital signal processors (Digital Signal Processor, DSP), application specific integrated circuits (Application Specific Integrated Circuit, ASIC), Field-Programmable Gate Array (Field-Programmable Gate Array, FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. The general-purpose processor can be a microprocessor or the processor can also be any conventional processor, etc., and the processor is the control center of the terminal equipment for measuring the thermal effect of the urban surface, and connects the entire urban surface with various interfaces and lines. Spatiotemporal measurement of thermal effects on various parts of the terminal equipment.

所述存储器可用于存储所述计算机程序和/或模块,所述处理器通过运行或执行存储在所述存储器内的计算机程序和/或模块,以及调用存储在存储器内的数据,实现所述城市地表热效应时空测度终端设备的各种功能。所述存储器可主要包括存储程序区和存储数据区,其中,存储程序区可存储操作系统、至少一个功能所需的应用程序;存储数据区可存储根据手机的使用所创建的数据等。此外,存储器可以包括高速随机存取存储器,还可以包括非易失性存储器,例如硬盘、内存、插接式硬盘,智能存储卡(Smart Media Card,SMC),安全数字(Secure Digital,SD)卡,闪存卡(Flash Card)、至少一个磁盘存储器件、闪存器件、或其他易失性固态存储器件。The memory can be used to store the computer programs and/or modules, and the processor realizes the city by running or executing the computer programs and/or modules stored in the memory and calling the data stored in the memory. Various functions of terminal equipment for space-time measurement of surface thermal effect. The memory may mainly include a program storage area and a data storage area, wherein the program storage area may store an operating system and at least one application required by a function; the data storage 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 hard disk, internal memory, plug-in hard disk, smart memory card (Smart Media Card, SMC), secure digital (Secure Digital, SD) card , a flash memory card (Flash Card), at least one magnetic disk storage device, flash memory device, or other volatile solid-state storage devices.

本发明还提供一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时实现本发明实施例上述方法的步骤。The present invention also provides a computer-readable storage medium, where the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the steps of the above method in the embodiment of the present invention are implemented.

所述城市地表热效应时空测度终端设备集成的模块/单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明实现上述实施例方法中的全部或部分流程,也可以通过计算机程序来指令相关的硬件来完成,所述的计算机程序可存储于一计算机可读存储介质中,该计算机程序在被处理器执行时,可实现上述各个方法实施例的步骤。其中,所述计算机程序包括计算机程序代码,所述计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。所述计算机可读介质可以包括:能够携带所述计算机程序代码的任何实体或装置、记录介质、U盘、移动硬盘、磁碟、光盘、计算机存储器、只读存储器(ROM,Read-OnlyMemory)、随机存取存储器(RAM,Random Access Memory)以及软件分发介质等。If the integrated modules/units of the spatio-temporal measurement terminal equipment for urban surface thermal effects are realized in the form of software function units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the present invention realizes all or part of the processes in the methods of the above embodiments, and can also be completed by instructing related hardware through a computer program. The computer program can be stored in a computer-readable storage medium, and the computer When the program is executed by the processor, the steps in the above-mentioned various method embodiments can be realized. Wherein, the computer program includes computer program code, and the computer program code may be in the form of source code, object code, executable file or some intermediate form. The computer-readable medium may include: any entity or device capable of carrying the computer program code, a recording medium, a U disk, a removable hard disk, a magnetic disk, an optical disk, a computer memory, a read-only memory (ROM, Read-OnlyMemory), Random Access Memory (RAM, Random Access Memory) and software distribution media, etc.

尽管结合优选实施方案具体展示和介绍了本发明,但所属领域的技术人员应该明白,在不脱离所附权利要求书所限定的本发明的精神和范围内,在形式上和细节上可以对本发明做出各种变化,均为本发明的保护范围。Although the present invention has been particularly shown and described in connection with preferred embodiments, it will be understood by those skilled in the art that changes in form and details may be made to the present invention without departing from the spirit and scope of the invention as defined by the appended claims. Making various changes is within the protection scope of the present invention.

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:
Figure FDA0003890942050000021
the calculation formula of the comprehensive surface heat intensity index is as follows:
Figure FDA0003890942050000022
the calculation formula of the spatial surface heat amplitude index is as follows:
Figure FDA0003890942050000023
the calculation formula of the spatial surface heat intensity index is as follows:
Figure FDA0003890942050000024
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,
Figure FDA0003890942050000025
indicates the area of the high-temperature region contained in the window with the pixel i as the center,
Figure FDA0003890942050000026
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.
CN202211259580.7A 2022-10-14 2022-10-14 Urban surface thermal effect space-time measurement method, terminal device and storage medium Pending CN115628814A (en)

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)

* Cited by examiner, † Cited by third party
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

Cited By (5)

* Cited by examiner, † Cited by third party
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
CN115628814A (en) Urban surface thermal effect space-time measurement method, terminal device and storage medium
Zhou et al. Deceleration of China’s human water use and its key drivers
Li et al. The surface urban heat island response to urban expansion: A panel analysis for the conterminous United States
Tan et al. Urban land expansion and arable land loss in China—a case study of Beijing–Tianjin–Hebei region
Liu et al. Assessing the coordination between economic growth and urban climate change in China from 2000 to 2015
Chen et al. Precipitation dominants synergies and trade-offs among ecosystem services across the Qinghai-Tibet Plateau
Sun et al. Mapping urban socio-economic vulnerability related to heat risk: A grid-based assessment framework by combing the geospatial big data
Geng et al. Spatiotemporal patterns and driving forces of remotely sensed urban agglomeration heat islands in South China
Zeng et al. Integrating anthropogenic heat emissions and cooling accessibility to explore environmental justice in heat-related health risks in Shanghai, China
Sigauke et al. Modelling non-stationary time series using a peaks over threshold distribution with time varying covariates and threshold: An application to peak electricity demand
CN113344291A (en) Urban inland inundation range forecasting method, device, medium and equipment
Liu et al. China Building Rooftop Area: the first multi-annual (2016–2021) and high-resolution (2.5 m) building rooftop area dataset in China derived with super-resolution segmentation from Sentinel-2 imagery
Li et al. Identifying and characterizing frequency and maximum durations of surface urban heat and cool island across global cities
CN110119556A (en) A spatiotemporal evolution analysis method of regional water conservation function
Yin et al. Changes in global heat waves and its socioeconomic exposure in a warmer future
Wang et al. Regulation of urban morphology on thermal environment across global cities
Yuan et al. Using street view images and a geographical detector to understand how street-level built environment is associated with urban poverty: A case study in Guangzhou
Liang et al. Seasonal impacts of built environment and its interactions on urban park cooling effects in Nanjing, China
CN116934073A (en) Refined measurement method of urban disaster resilience based on spatiotemporal activity analysis
Li et al. Spatial evolution of cultivated land in the Heilongjiang Province in China from 1980 to 2015
Zhang et al. Estimating 1 km gridded daily air temperature using a spatially varying coefficient model with sign preservation
Sharma et al. Environmental change and groundwater variability in South Bihar, India
CN118248035A (en) Method and device for generating rain and flood disaster risk map by coupling intelligent algorithm and multi-source data
Liu et al. CBRA: The first multi-annual (2016–2021) and high-resolution (2.5 m) building rooftop area dataset in China derived with Super-resolution Segmentation from Sentinel-2 imagery
Zhang et al. Trade-offs and synergies of ecosystem services and their threshold effects in the largest tableland of the Loess Plateau

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