CN115630535B - Urban surface heat island strength dynamic quantification method and device and electronic equipment - Google Patents

Urban surface heat island strength dynamic quantification method and device and electronic equipment Download PDF

Info

Publication number
CN115630535B
CN115630535B CN202211632782.1A CN202211632782A CN115630535B CN 115630535 B CN115630535 B CN 115630535B CN 202211632782 A CN202211632782 A CN 202211632782A CN 115630535 B CN115630535 B CN 115630535B
Authority
CN
China
Prior art keywords
urban
dynamic
surface temperature
area
rural
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202211632782.1A
Other languages
Chinese (zh)
Other versions
CN115630535A (en
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 Agricultural Resources and Regional Planning of CAAS
Hebei GEO University
Original Assignee
Institute of Agricultural Resources and Regional Planning of CAAS
Hebei GEO University
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 Agricultural Resources and Regional Planning of CAAS, Hebei GEO University filed Critical Institute of Agricultural Resources and Regional Planning of CAAS
Priority to CN202211632782.1A priority Critical patent/CN115630535B/en
Publication of CN115630535A publication Critical patent/CN115630535A/en
Application granted granted Critical
Publication of CN115630535B publication Critical patent/CN115630535B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/24—Querying
    • G06F16/245—Query processing
    • G06F16/2458—Special types of queries, e.g. statistical queries, fuzzy queries or distributed queries
    • G06F16/2474—Sequence data queries, e.g. querying versioned data
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/29—Geographical information databases
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/10—Services
    • G06Q50/26—Government or public services
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08—Thermal analysis or thermal optimisation
    • 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

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • Databases & Information Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Tourism & Hospitality (AREA)
  • Educational Administration (AREA)
  • Health & Medical Sciences (AREA)
  • Software Systems (AREA)
  • Computational Linguistics (AREA)
  • Mathematical Physics (AREA)
  • Evolutionary Computation (AREA)
  • Remote Sensing (AREA)
  • Computer Hardware Design (AREA)
  • Development Economics (AREA)
  • Fuzzy Systems (AREA)
  • Geometry (AREA)
  • Probability & Statistics with Applications (AREA)
  • Economics (AREA)
  • General Health & Medical Sciences (AREA)
  • Human Resources & Organizations (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • General Business, Economics & Management (AREA)
  • Cooling Or The Like Of Semiconductors Or Solid State Devices (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Abstract

The invention belongs to the technical field of urban heat island analysis, and relates to a dynamic quantification method and device for urban surface heat island strength and electronic equipment, wherein the method comprises the following steps: extracting year-by-year dynamic city natural boundaries; extracting year-by-year dynamic city center points; extracting dynamic country natural boundaries year by year; respectively extracting the average instantaneous surface temperature of the urban representative area and the highest surface temperature of the urban and rural overall area by using the instantaneous surface temperature remote sensing data and fitting the Gaussian curved surface twice; determining an average instantaneous surface temperature of the rural area; constructing a dynamic instantaneous surface heat island intensity calculation model; and determining the dynamic surface heat island strength. The method can solve the problem of single urban boundary in the dynamic earth surface heat island intensity quantization model, realize the rapid and accurate extraction of the dynamic urban natural boundary and the dynamic rural natural boundary, and improve the accuracy of the average instantaneous earth surface temperature of the urban area and the uniformity of the dynamic earth surface heat island intensity quantization.

Description

Urban surface heat island strength dynamic quantification method and device and electronic equipment
Technical Field
The invention relates to the technical field of urban heat island analysis, in particular to a method and a device for dynamically quantifying the intensity of an earth surface heat island considering urban expansion and electronic equipment.
Background
With the acceleration of the global urban process, the urban heat island phenomenon becomes one of typical urban heat environment problems, and is a key factor affecting climate change, urban living environment and health.
The existing ground surface heat island intensity quantization model often ignores the expansion and merging process of the urban boundary, so that the selection of the background temperature reference value has great uncertainty, and the long-term change characteristics of the heat island intensity and range cannot be accurately reflected. In general, the surface heat island quantification method lacks unified standards, and the monitoring precision of the surface heat island phenomenon is low.
By means of relatively mature city clustering algorithm, researchers develop a method for defining fixed city and village range, and adopt a once temperature field model to quantify the intensity of the earth surface heat island. However, these models do not fully take into account the dynamic characteristics of the urban natural boundaries, thereby affecting the accuracy of extraction of the rural natural boundaries. In addition, the existing temperature field model is characterized by describing the distribution pattern characteristics of the surface temperature, constructing an ellipsoid to obtain temperature field parameters, and accordingly representing the surface heat island strength. However, the model has larger sensitivity to the selection of background temperature, the effect of dynamic urban and rural boundaries generated by urban expansion cannot be fully shown, the influence of urban edge transition zones cannot be eliminated by single ellipsoid parameters, and the quantification of the ground surface heat island strength should fully represent the maximum natural difference of urban and rural ground surfaces. Therefore, the traditional primary Gaussian field model based on the static urban and rural boundaries cannot effectively monitor the dynamic surface heat island strength.
Disclosure of Invention
In order to solve the technical problems in the current earth surface heat island intensity quantization method, the invention fully considers the dynamic change of urban and rural areas, and provides a dynamic urban and rural boundary extraction method, a device and electronic equipment which take urban expansion into consideration by quantizing the spatial distribution pattern of urban, suburban and rural earth surface temperatures. The method is used for improving the application mode of the Gaussian model in the earth surface temperature field, a dynamic earth surface heat island intensity quantization model is constructed through a two-time Gaussian temperature curved surface fitting process, the long-time accurate quantization of the instantaneous earth surface heat island intensity is realized, the method is suitable for the consistency quantization of the instantaneous earth surface heat island intensity of different areas at different moments, and further the monitoring precision and the application value of the dynamic earth surface heat island are improved.
In a first aspect, the invention provides a dynamic quantification method for urban surface heat island intensity, comprising the following steps:
constructing a long-time-sequence remote sensing land utilization classification database, and extracting dynamic city natural boundaries year by year;
based on the natural boundary of the city, extracting dynamic city center points year by using a discrete point group minimum maximum distance method;
acquiring expansion of a city natural boundary and displacement of a city center point, and extracting dynamic country natural boundaries year by year;
respectively extracting the average instantaneous surface temperature of the urban representative area and the highest surface temperature of the urban and rural overall area by using the instantaneous surface temperature remote sensing data and fitting the Gaussian curved surface twice;
determining the average instantaneous surface temperature of the rural area according to the highest surface temperature of the urban and rural overall area and the surface temperature of the urban central point;
based on the dynamic urban natural boundary and the dynamic rural natural boundary, constructing a dynamic instantaneous surface heat island intensity calculation model by utilizing the average instantaneous surface temperature of the urban representative area and the average instantaneous surface temperature of the rural area;
and determining the dynamic surface heat island intensity according to the dynamic instantaneous surface heat island intensity calculation model.
The invention provides a dynamic quantification device for urban surface heat island intensity, which comprises a first extraction unit, a second extraction unit, a third extraction unit, a fitting unit, a processing unit, a model construction unit and an output unit;
the first extraction unit is used for constructing a long-time-sequence remote sensing land utilization classification database and extracting dynamic city natural boundaries year by year;
the second extraction unit is used for extracting dynamic city center points year by utilizing a discrete point group minimum maximum distance method based on the city natural boundary;
the third extraction unit is used for obtaining the expansion of the natural boundary of the city and the displacement of the center point of the city, and extracting the dynamic natural boundary of the country year by year;
the fitting unit is used for respectively extracting the average instantaneous surface temperature of the urban representative area and the highest surface temperature of the urban and rural overall area by using the instantaneous surface temperature remote sensing data and fitting the Gaussian curved surface twice;
the processing unit is used for determining the average instantaneous surface temperature of the rural area according to the highest surface temperature of the urban and rural overall area and the surface temperature of the urban central point;
the model construction unit is used for constructing a dynamic instantaneous surface heat island intensity calculation model by utilizing the average instantaneous surface temperature of the urban representative area and the average instantaneous surface temperature of the rural area based on the dynamic urban natural boundary and the dynamic rural natural boundary;
and the output unit is used for determining the dynamic surface heat island intensity according to the dynamic instantaneous surface heat island intensity calculation model.
In a third aspect, the present invention provides an electronic device, comprising:
a processor and a memory;
the memory is used for storing computer operation instructions;
and the processor is used for executing the urban surface heat island strength dynamic quantification method by calling the computer operation instruction.
The beneficial effects of the invention are as follows:
(1) In the dynamic city range extraction, the dynamic change of the city range caused by the expansion of city pixels is fully considered, and the natural boundaries of the city are extracted year by year, so that the problem of single city boundary in the existing dynamic earth surface heat island intensity quantization model can be solved;
(2) The buffer area at the edge of the natural boundary of the city is used as a medium, the influence of slow transition of the physical properties of the natural earth surface of the peripheral buffer area of the city and the natural earth surface of the rural area is eliminated, the dynamic rural area range can be effectively extracted, and the natural boundary of the rural area with the corresponding physical meaning with the natural boundary of the city is obtained; the expansion of the natural boundaries of the cities and the urban center point displacement effect caused by the expansion are considered, the dynamic urban center point is extracted based on the urban range, and rural areas with the same area as the urban range are defined year by year, so that the rapid and accurate extraction of the natural boundaries of the dynamic cities and the dynamic rural natural boundaries is realized, and a practical and reliable monitoring range is provided for the ground heat island intensity quantization model;
(3) Aiming at the problem of larger uncertainty of a reference background value in a Gaussian model, starting from a quantization method for improving the average surface temperature of a city, acquiring the highest surface temperature of the city and the surface temperature of a annual city center point by combining a once Gaussian surface model, and improving the accuracy of the average instantaneous surface temperature of a representative city region;
(4) And carrying out Gaussian surface secondary fitting on a larger city overall range and different background temperatures, calculating the highest value of the surface temperature of the urban and rural overall area, acquiring the average surface temperature of the representative rural area by combining the surface temperature of the urban center year by year, and finally constructing an intuitive surface heat island intensity quantization model to realize the uniformity of the instantaneous dynamic surface heat island intensity quantization.
On the basis of the technical scheme, the invention can be improved as follows.
Further, the dynamic urban natural boundary is extracted, a city clustering method is adopted, and a breadth-first search algorithm is utilized to obtain the minimum convex hull, so that the dynamic urban natural boundary is obtained.
Further, the extracting of the city center point, based on the city natural boundary, uses a discrete point group minimum maximum distance method to extract a year-by-year dynamic city center point, including:
acquiring longitude and latitude coordinates of all pixel centers in the city natural boundary, acquiring distances between each pixel center and the rest of pixel center points by taking each pixel center as a starting point, determining a maximum distance set of each starting point and the rest of pixel center points, and screening out a minimum value from the maximum distance set, wherein the starting point where the minimum value is located is used as the city center point.
Further, the extracting of the natural boundaries of the country includes:
and extracting the minimum convex hull which takes the urban central point as the center and is outside the suburban area and the urban area year by adopting a buffer area calculation method as the rural natural area, and determining the rural natural boundary.
Further, by using the instantaneous surface temperature remote sensing data, the average instantaneous surface temperature of the urban representative area and the highest surface temperature of the urban and rural overall area are respectively extracted by fitting the Gaussian curved surface twice, and the method comprises the following steps:
performing first fitting on the urban natural boundary and the total Gaussian surface of the buffer area with the set range of the edge of the urban natural boundary to obtain the highest surface temperature of the city;
averaging the highest surface temperature with the surface temperature of the urban central point to obtain the average instantaneous surface temperature of the urban representative area;
and performing second fitting on the urban natural boundary, the urban natural boundary edge set range buffer area and the rural natural boundary area total Gaussian temperature field to obtain the highest surface temperature of the urban and rural total area.
Further, determining an average instantaneous surface temperature of the rural area based on the highest surface temperature of the urban and rural overall area and the urban central point surface temperature, comprising:
and calculating the difference between the highest surface temperature of the urban and rural overall area and the surface temperature of the urban central point to obtain the average instantaneous surface temperature of the rural area.
Further, based on the dynamic city natural boundary and the dynamic country natural boundary, constructing a dynamic instantaneous surface heat island intensity calculation model by using the average instantaneous surface temperature of the representative city region and the average instantaneous surface temperature of the country region, including: and calculating the difference between the average instantaneous surface temperature of the representative urban area of the year-by-year dynamic state and the average instantaneous surface temperature of the country area of the year-by-year dynamic state to obtain a dynamic instantaneous surface heat island intensity calculation model.
Drawings
FIG. 1 is a flow chart of a dynamic quantification method for urban surface heat island intensity provided in embodiment 1 of the present invention;
fig. 2 is a schematic diagram of a specific embodiment of a dynamic quantification method for urban surface heat island intensity provided in embodiment 1 of the present invention;
FIG. 3 is a schematic diagram of a dynamic quantification device for urban surface heat island intensity provided in embodiment 2 of the present invention;
fig. 4 is a schematic diagram of an electronic device according to embodiment 3 of the present invention.
Icon: 40-an electronic device; 410-a processor; 420-bus; 430-memory; 440-transceiver.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention, and it is apparent that the described embodiments are some embodiments of the present invention, but not all embodiments of the present invention. The components of the embodiments of the present invention generally described and illustrated in the figures herein may be arranged and designed in a wide variety of different configurations.
Example 1
As an embodiment, as shown in fig. 1, to solve the above technical problem, the present embodiment provides a method for dynamically quantifying the intensity of an urban ground heat island, including:
constructing a long-time-sequence remote sensing land utilization classification database, and extracting dynamic city natural boundaries year by year;
based on the natural boundary of the city, extracting dynamic city center points year by utilizing a discrete point group minimum maximum distance method;
acquiring expansion of a city natural boundary and displacement of a city center point, and extracting dynamic country natural boundaries year by year;
respectively extracting the average instantaneous surface temperature of the urban representative area and the highest surface temperature of the urban and rural overall area by using the instantaneous surface temperature remote sensing data and fitting the Gaussian curved surface twice;
determining the average instantaneous surface temperature of the rural area according to the highest surface temperature of the urban and rural overall area and the surface temperature of the urban central point;
based on the dynamic city natural boundary and the dynamic country natural boundary, constructing a dynamic instantaneous surface heat island intensity calculation model by utilizing the average instantaneous surface temperature of the representative region of the city and the average instantaneous surface temperature of the country region;
and determining the dynamic surface heat island intensity according to the dynamic instantaneous surface heat island intensity calculation model.
As shown in fig. 2, a schematic diagram of a specific embodiment of a dynamic quantification method for urban surface heat island intensity is shown.
Optionally, extracting the dynamic city natural boundary, adopting a city clustering method, and acquiring the minimum convex hull by using a breadth-first search algorithm to obtain the dynamic city natural boundary.
Optionally, the extracting of the city center point is based on the city natural boundary, and the extracting of the year-by-year dynamic city center point by using a discrete point group minimum maximum distance method comprises the following steps:
acquiring longitude and latitude coordinates of centers of all pixels in a natural boundary of the city, acquiring distances between each starting point and center points of all other pixels by taking each pixel center as a starting point, determining a maximum distance set of each starting point and center points of all other pixels, and screening out a minimum value from the maximum distance set, wherein the starting point where the minimum value is located is used as the center point of the city.
The specific calculation formula of the minimum maximum distance method is as follows:
Figure 14553DEST_PATH_IMAGE001
；
Figure 213585DEST_PATH_IMAGE002
；
wherein, the liquid crystal display device comprises a liquid crystal display device,
Figure 524480DEST_PATH_IMAGE003
is the first in the natural boundary of the city
Figure 654110DEST_PATH_IMAGE003
The center point of each pixel is defined by the center point,
Figure 999641DEST_PATH_IMAGE004
for all the number of picture elements.
Figure 618841DEST_PATH_IMAGE005
Is the first
Figure 620295DEST_PATH_IMAGE003
Center point (S) and (H)
Figure 366665DEST_PATH_IMAGE006
The straight line distance of the individual center points,
Figure 769965DEST_PATH_IMAGE007
taking all maximum distances for all linear distances
Figure 560066DEST_PATH_IMAGE008
Minimum value of (2)
Figure 48817DEST_PATH_IMAGE009
The pixel where it is
Figure 582566DEST_PATH_IMAGE003
Namely, the city center point.
Optionally, the extracting of the natural boundaries of the country includes:
and extracting the minimum convex hull which takes the urban central point as the center and is outside the suburban area and the urban area year by adopting a buffer area calculation method as a rural natural area to determine the rural natural boundary.
The urban boundary expansion and the urban central point displacement effect generated by the urban boundary expansion are considered, a buffer area calculation method is adopted on the basis of eliminating the edge of the urban natural boundary, which is the outer distance length (such as one kilometer) from year to year, as a suburban area, and the minimum convex hull which takes the dynamic urban central point as the center and is outside the suburban area and has the same area as the city is extracted year by year to be used as the dynamic rural natural area.
Optionally, the method for respectively extracting the average instantaneous surface temperature of the urban representative area and the highest surface temperature of the urban and rural overall area by using the instantaneous surface temperature remote sensing data through twice fitting the gaussian curved surface comprises the following steps:
performing first fitting on the urban natural boundary and the total Gaussian surface of the buffer area with the set range of the edge of the urban natural boundary to obtain the highest surface temperature of the city;
averaging the highest surface temperature with the surface temperature of the urban central point to obtain the average instantaneous surface temperature of the urban representative area;
and performing second fitting on the city natural boundary, the buffer area with the set range of the city natural boundary edge and the total Gaussian temperature field of the rural natural boundary area to obtain the highest surface temperature of the urban and rural total area.
In the practical application process, the first Gaussian surface fitting is to fit the natural boundary of the city and the total Gaussian temperature field of the buffer area of one kilometer at the edge of the natural boundary of the city, and the average value of all the surface temperatures of the buffer area is set as
Figure 840372DEST_PATH_IMAGE010
The center strength of the ellipsoid is
Figure 270216DEST_PATH_IMAGE011
The surface temperature of the urban central point is
Figure 308580DEST_PATH_IMAGE012
Fitting the shape and distance parameters of the Gaussian ellipsoid to obtain the center strength
Figure 318124DEST_PATH_IMAGE011
And the surface temperature of the central point of the city
Figure 263993DEST_PATH_IMAGE013
Combining, finally extracting the average instantaneous surface temperature of the representative area of city
Figure 395898DEST_PATH_IMAGE014
The formula is as follows:
Figure 859240DEST_PATH_IMAGE015
；
the second Gaussian surface fitting is to fit the city natural boundary, the buffer area of one kilometer at the edge of the city natural boundary and the total Gaussian temperature field of the country natural boundary area. Specifically, the average value of all the surface temperatures in the rural area is taken as the background temperature
Figure 734792DEST_PATH_IMAGE016
Fitting again to obtain new ellipsoidal central intensity as the highest surface temperature of urban and rural overall area
Figure 232769DEST_PATH_IMAGE017
。
Optionally, determining the average instantaneous surface temperature of the rural area according to the highest surface temperature of the urban and rural overall area and the surface temperature of the urban central point comprises:
and calculating the difference between the highest surface temperature of the urban and rural overall area and the surface temperature of the urban central point to obtain the average instantaneous surface temperature of the rural area.
Calculating the highest surface temperature of urban and rural overall areas
Figure 738837DEST_PATH_IMAGE017
And the surface temperature of the central point of the city
Figure 486213DEST_PATH_IMAGE018
Is the average instantaneous surface temperature of the rural area
Figure 103139DEST_PATH_IMAGE019
The formula is as follows:
Figure 3093DEST_PATH_IMAGE020
。
optionally, based on the dynamic city natural boundary and the dynamic country natural boundary, constructing a dynamic instantaneous surface heat island intensity calculation model by using the average instantaneous surface temperature of the representative region of the city and the average instantaneous surface temperature of the country region, including: calculating average instantaneous surface temperature for year-by-year dynamic urban representative regions
Figure 680062DEST_PATH_IMAGE021
Average instantaneous surface temperature with year-by-year dynamic rural areas
Figure 117997DEST_PATH_IMAGE022
And obtaining the dynamic instantaneous surface heat island intensity by the difference.
Let dynamic instantaneous surface heat island intensity be
Figure 335352DEST_PATH_IMAGE023
Then:
Figure 542342DEST_PATH_IMAGE024
。
by combining dynamic urban natural boundary and rural natural boundary extraction results, a dynamic earth surface heat island intensity quantization model is constructed, and urban and rural total area average highest earth surface temperature extraction and heat island intensity urban and rural temperature difference method are effectively combined, so that the method is an improvement of an application mode of a Gaussian curved surface in an earth surface heat island intensity quantization scene.
Is provided with
Figure 452529DEST_PATH_IMAGE025
Is the longitude and latitude of all pixels in the Gaussian surface temperature field,
Figure 112181DEST_PATH_IMAGE026
is a picture element
Figure 133226DEST_PATH_IMAGE025
Is set to be a surface temperature of the earth,
Figure 929144DEST_PATH_IMAGE027
is the surface temperature of the background pixel,
Figure 213495DEST_PATH_IMAGE028
is the center of the ellipsoid,
Figure 173492DEST_PATH_IMAGE029
is a half major axis of an ellipsoid,
Figure 935911DEST_PATH_IMAGE030
is the semi-minor axis of an ellipsoid,
Figure 914232DEST_PATH_IMAGE031
orientation of ellipsoid, coefficient
Figure 838325DEST_PATH_IMAGE032
Then the center strength of the ellipsoid is taken as the highest surface temperature of the city, and then:
Figure 534886DEST_PATH_IMAGE033
Figure 835417DEST_PATH_IMAGE034
Figure 933823DEST_PATH_IMAGE035
。
the above embodiment has the following advantages:
(1) In the dynamic city range extraction, the dynamic change of the city range caused by the expansion of city pixels is fully considered, and the natural boundaries of the city are extracted year by year, so that the problem of single city boundary in the existing dynamic earth surface heat island intensity quantization model can be solved;
(2) The buffer area at the edge of the natural boundary of the city is used as a medium, the influence of slow transition of the physical properties of the natural earth surface of the peripheral buffer area of the city and the natural earth surface of the rural area is eliminated, the dynamic rural area range can be effectively extracted, and the natural boundary of the rural area with the corresponding physical meaning with the natural boundary of the city is obtained; the expansion of the natural boundaries of the cities and the urban center point displacement effect caused by the expansion are considered, the dynamic urban center point is extracted based on the urban range, and rural areas with the same area as the urban range are defined year by year, so that the rapid and accurate extraction of the natural boundaries of the dynamic cities and the dynamic rural natural boundaries is realized, and a practical and reliable monitoring range is provided for the ground heat island intensity quantization model;
(3) Aiming at the problem of larger uncertainty of a reference background value in a Gaussian model, starting from a quantization method for improving the average surface temperature of a city, combining the highest surface temperature of the city and the surface temperature of a annual city center point obtained by a once Gaussian surface model, and improving the accuracy of the average instantaneous surface temperature of a city region;
(4) And carrying out Gaussian surface secondary fitting on a larger city overall range and different background temperatures, calculating the highest value of the earth surface temperature of the urban and rural overall area, acquiring the average earth surface temperature of the rural area by combining the annual urban central earth surface temperature, and finally constructing an intuitive earth surface heat island intensity quantization model to realize the uniformity of dynamic earth surface heat island intensity quantization.
Example 2
Based on the same principle as the method shown in the embodiment 1 of the invention, the embodiment of the invention also provides a dynamic quantification device for the urban surface heat island intensity, which comprises a first extraction unit, a second extraction unit, a third extraction unit, a fitting unit, a processing unit, a model construction unit and an output unit;
the first extraction unit is used for constructing a long-time-sequence remote sensing land utilization classification database and extracting dynamic city natural boundaries year by year;
the second extraction unit is used for extracting dynamic city center points year by utilizing a discrete point group minimum maximum distance method based on city natural boundaries;
the third extraction unit is used for obtaining the expansion of the urban natural boundary and the displacement of the urban central point and extracting the dynamic country natural region year by year;
the fitting unit is used for respectively extracting the average instantaneous surface temperature of the urban representative area and the highest surface temperature of the urban and rural overall area by utilizing the instantaneous surface temperature remote sensing data and fitting the Gaussian curved surface twice;
the processing unit is used for determining the average instantaneous surface temperature of the rural area according to the highest surface temperature of the urban and rural overall area and the surface temperature of the urban central point;
the model construction unit is used for constructing a dynamic instantaneous surface heat island intensity calculation model by utilizing the average instantaneous surface temperature of the urban representative area and the average instantaneous surface temperature of the rural area based on the dynamic urban natural boundary and the dynamic rural natural boundary;
and the output unit is used for determining the dynamic surface heat island intensity according to the dynamic instantaneous surface heat island intensity calculation model.
Optionally, the first extraction unit extracts dynamic natural boundaries of cities by adopting a city clustering method, and obtaining the minimum convex hull by using a breadth-first search algorithm to obtain the dynamic natural boundaries of cities.
Optionally, the second extracting unit extracts dynamic city center points year by year based on city natural boundaries and by using a discrete point group minimum maximum distance method, and the second extracting unit comprises:
acquiring longitude and latitude coordinates of centers of all pixels in a natural boundary of the city, acquiring distances between each starting point and center points of all other pixels by taking each pixel center as a starting point, determining a maximum distance set of each starting point and center points of all other pixels, and screening out a minimum value from the maximum distance set, wherein the starting point where the minimum value is located is used as the center point of the city.
Optionally, the third extraction unit extracts natural boundaries of the country, including:
and extracting the minimum convex hull which takes the urban central point as the center and is outside the suburban area and the urban area year by adopting a buffer area calculation method as a rural natural area to determine the rural natural boundary.
Optionally, the fitting unit respectively extracts the average instantaneous surface temperature of the urban representative area and the highest surface temperature of the urban and rural overall area by using the instantaneous surface temperature remote sensing data and fitting the gaussian curved surface twice, including:
performing first fitting on the urban natural boundary and the total Gaussian surface of the buffer area with the set range of the edge of the urban natural boundary to obtain the highest surface temperature of the city;
averaging the highest surface temperature with the surface temperature of the urban central point to obtain the average instantaneous surface temperature of the urban representative area;
and performing second fitting on the city natural boundary, the buffer area with the set range of the city natural boundary edge and the total Gaussian temperature field of the rural natural boundary area to obtain the highest surface temperature of the urban and rural total area.
Optionally, the processing unit determines an average instantaneous surface temperature of the rural area according to the highest surface temperature of the urban and rural overall area and the surface temperature of the urban central point, including:
and calculating the difference between the highest surface temperature of the urban and rural overall area and the surface temperature of the urban central point to obtain the average instantaneous surface temperature of the rural area.
Optionally, the model building unit builds the dynamic transient earth surface heat island intensity calculation model based on the dynamic city natural boundary and the dynamic country natural boundary by using the average transient earth surface temperature of the representative region of the city and the average transient earth surface temperature of the country region, and the model building unit comprises: and calculating the difference between the average instantaneous surface temperature of the year-by-year dynamic city representative region and the average instantaneous surface temperature of the year-by-year dynamic rural region to obtain a dynamic instantaneous surface heat island intensity calculation model.
Example 3
Based on the same principle as the method shown in the embodiment of the present invention, there is also provided an electronic device in the embodiment of the present invention, as shown in fig. 4, which may include, but is not limited to: a processor and a memory; a memory for storing a computer program; a processor for executing the method according to any of the embodiments of the invention by invoking a computer program.
In an alternative embodiment, an electronic device is provided, the electronic device 40 shown in fig. 4 comprising: a processor 410 and a memory 430. Processor 410 is coupled to memory 430, such as via bus 420.
Optionally, the electronic device 40 may further comprise a transceiver 440, and the transceiver 440 may be used for data interaction between the electronic device and other electronic devices, such as transmission of data and/or reception of data, etc. It should be noted that, in practical applications, the transceiver 440 is not limited to one, and the structure of the electronic device 40 is not limited to the embodiment of the present invention.
The processor 410 may be a CPU central processing unit, a general purpose processor, a DSP data signal processor, an ASIC specific integrated circuit, an FPGA field programmable gate array or other programmable logic device, a hardware component, or any combination thereof. Processor 410 may also be a combination that implements computing functionality, e.g., comprising one or more microprocessor combinations, a combination of a DSP and a microprocessor, etc.
Bus 420 may include a path to transfer information between the aforementioned components. Bus 420 may be a PCI peripheral component interconnect standard bus or an EISA extension industry standard architecture bus, among others. Bus 420 may be divided into a control bus, a data bus, an address bus, and the like. For ease of illustration, only one thick line is shown in fig. 4, but not only one bus or one type of bus.
Memory 430 may be, but is not limited to, ROM read-only memory or other type of static storage device that can store static information and instructions, RAM random-access memory or other type of dynamic storage device that can store information and instructions, EEPROM electrically erasable programmable read-only memory, CD-ROM read-only or other optical disk storage, optical disk storage (including optical disks, laser disks, compact disks, digital versatile disks, etc.), magnetic disk storage media, or any other medium that can be used to carry or store desired program code in the form of instructions or data structures and that can be accessed by a computer.
The memory 430 is used to store application program codes (computer programs) for executing the inventive arrangements and is controlled to be executed by the processor 410. The processor 410 is configured to execute application code stored in the memory 430 to implement what is shown in the foregoing method embodiments.
The above is only a preferred embodiment of the present invention, and is not intended to limit the present invention, but various modifications and variations can be made to the present invention by those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims (8)

1. The urban earth surface heat island strength dynamic quantification method is characterized by comprising the following steps of:
constructing a long-time-sequence remote sensing land utilization classification database, and extracting dynamic city natural boundaries year by year;
based on the natural boundary of the city, extracting dynamic city center points year by using a discrete point group minimum maximum distance method;
acquiring expansion of a city natural boundary and displacement of a city center point, and extracting dynamic country natural boundaries year by year;
respectively extracting the average instantaneous surface temperature of the urban representative area and the highest surface temperature of the urban and rural overall area by using the instantaneous surface temperature remote sensing data through twice fitting Gaussian curved surfaces, wherein the method comprises the following steps of: performing first fitting on the urban natural boundary and the total Gaussian surface of the buffer area with the set range of the edge of the urban natural boundary to obtain the highest surface temperature of the city; averaging the highest surface temperature with the surface temperature of the urban central point to obtain the average instantaneous surface temperature of the urban representative area; performing second fitting on the urban natural boundary, the urban natural boundary edge set range buffer area and the rural natural boundary area total Gaussian temperature field to obtain the highest surface temperature of the urban and rural total area;
determining the average instantaneous surface temperature of the rural area according to the highest surface temperature of the urban and rural overall area and the surface temperature of the urban central point;
based on the dynamic urban natural boundary and the dynamic rural natural boundary, constructing a dynamic instantaneous surface heat island intensity calculation model by utilizing the average instantaneous surface temperature of the urban representative area and the average instantaneous surface temperature of the rural area;
and determining the dynamic surface heat island intensity according to the dynamic instantaneous surface heat island intensity calculation model.
2. The urban surface heat island intensity dynamic quantification method according to claim 1, wherein the dynamic urban natural boundary is extracted by adopting a city clustering method, and a breadth-first search algorithm is utilized to obtain a minimum convex hull, so that the dynamic urban natural boundary is obtained.
3. The method for dynamically quantifying urban surface heat island intensity according to claim 1, wherein the extracting of the urban center point is based on the urban natural boundary, and the extracting of the annual dynamic urban center point by using a discrete point group minimum maximum distance method comprises:
acquiring longitude and latitude coordinates of all pixel centers in the city natural boundary, acquiring distances between each pixel center and the rest of pixel center points by taking each pixel center as a starting point, determining a maximum distance set of each starting point and the rest of pixel center points, and screening out a minimum value from the maximum distance set, wherein the starting point where the minimum value is located is used as the city center point.
4. The method for dynamically quantifying urban surface heat island intensity according to claim 1, wherein the extracting of the rural natural boundary comprises:
and extracting the minimum convex hull which takes the urban central point as the center and is outside the suburban area and the urban area year by adopting a buffer area calculation method as the rural natural area, and determining the rural natural boundary.
5. The method for dynamically quantifying urban surface heat island intensity according to claim 1, wherein determining a rural area average instantaneous surface temperature from the highest surface temperature of the urban and rural overall area and the urban center surface temperature comprises:
and calculating the difference between the highest surface temperature of the urban and rural overall area and the surface temperature of the urban central point to obtain the average instantaneous surface temperature of the rural area.
6. The method for dynamically quantifying urban surface heat island intensity according to claim 1, wherein constructing a dynamic instantaneous surface heat island intensity calculation model using an average instantaneous surface temperature of the urban representative region and an average instantaneous surface temperature of the rural region based on the dynamic urban natural boundary and the dynamic rural natural boundary, comprises: and calculating the difference between the average instantaneous surface temperature of the representative urban area of the year-by-year dynamic state and the average instantaneous surface temperature of the country area of the year-by-year dynamic state to obtain a dynamic instantaneous surface heat island intensity calculation model.
7. The urban surface heat island intensity dynamic quantification device is characterized by comprising a first extraction unit, a second extraction unit, a third extraction unit, a fitting unit, a processing unit, a model construction unit and an output unit;
the first extraction unit is used for constructing a long-time-sequence remote sensing land utilization classification database and extracting dynamic city natural boundaries year by year;
the second extraction unit is used for extracting dynamic city center points year by utilizing a discrete point group minimum maximum distance method based on the city natural boundary;
the third extraction unit is used for obtaining the expansion of the natural boundary of the city and the displacement of the center point of the city, and extracting the dynamic natural boundary of the country year by year;
the fitting unit is configured to respectively extract an average instantaneous surface temperature of a representative area of the city and a highest surface temperature of an overall area of the city and the country by using instantaneous surface temperature remote sensing data and fitting gaussian curved surfaces twice, and includes: performing first fitting on the urban natural boundary and the total Gaussian surface of the buffer area with the set range of the edge of the urban natural boundary to obtain the highest surface temperature of the city; averaging the highest surface temperature with the surface temperature of the urban central point to obtain the average instantaneous surface temperature of the urban representative area; performing second fitting on the urban natural boundary, the urban natural boundary edge set range buffer area and the rural natural boundary area total Gaussian temperature field to obtain the highest surface temperature of the urban and rural total area;
the processing unit is used for determining the average instantaneous surface temperature of the rural area according to the highest surface temperature of the urban and rural overall area and the surface temperature of the urban central point;
the model construction unit is used for constructing a dynamic instantaneous surface heat island intensity calculation model by utilizing the average instantaneous surface temperature of the urban representative area and the average instantaneous surface temperature of the rural area based on the dynamic urban natural boundary and the dynamic rural natural boundary;
and the output unit is used for determining the dynamic surface heat island intensity according to the dynamic instantaneous surface heat island intensity calculation model.
8. An electronic device, comprising:
a processor and a memory;
the memory is used for storing computer operation instructions;
the processor is configured to perform the method of any one of claims 1 to 6 by invoking the computer operating instructions.
CN202211632782.1A 2022-12-19 2022-12-19 Urban surface heat island strength dynamic quantification method and device and electronic equipment Active CN115630535B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211632782.1A CN115630535B (en) 2022-12-19 2022-12-19 Urban surface heat island strength dynamic quantification method and device and electronic equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211632782.1A CN115630535B (en) 2022-12-19 2022-12-19 Urban surface heat island strength dynamic quantification method and device and electronic equipment

Publications (2)

Publication Number Publication Date
CN115630535A CN115630535A (en) 2023-01-20
CN115630535B true CN115630535B (en) 2023-05-12

Family

ID=84909705

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211632782.1A Active CN115630535B (en) 2022-12-19 2022-12-19 Urban surface heat island strength dynamic quantification method and device and electronic equipment

Country Status (1)

Country Link
CN (1) CN115630535B (en)

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103473613B (en) * 2013-09-09 2017-01-11 武汉理工大学 Landscape structure-surface temperature-electricity consumption coupling model and application thereof
CN107678075A (en) * 2017-11-13 2018-02-09 深圳先进技术研究院 A kind of urban heat land effect monitoring method and system based on domestic satellite
CN109959970B (en) * 2017-12-26 2020-10-09 中国科学院地理科学与资源研究所 Sky hemisphere thermal infrared atmospheric downlink radiation ground measurement method
CN108776360A (en) * 2018-06-13 2018-11-09 华南农业大学 A kind of method of urban heat island strength Monitoring on Dynamic Change
CN111637972B (en) * 2020-04-23 2021-07-13 中国科学院空天信息创新研究院 Remote sensing definition method for industrial heat island effect

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
史姝姝 ; 匡文慧 ; 董斯齐 ; .21世纪以来西安城乡梯度土地覆盖变化及对城市热岛影响时空特征.遥感技术与应用.2020,(第03期),全文. *

Also Published As

Publication number Publication date
CN115630535A (en) 2023-01-20

Similar Documents

Publication Publication Date Title
Zhang et al. Vectorized rooftop area data for 90 cities in China
Willmott et al. On the use of dimensioned measures of error to evaluate the performance of spatial interpolators
CN104619009B (en) A kind of method of adjustment, device and the mobile terminal in location data sampling period
Na et al. Bidirectional DEM relief shading method for extraction of gully shoulder line in loess tableland area
Wei et al. A sub-pixel edge detection algorithm based on Zernike moments
Li et al. Real-time flood simulations using CA model driven by dynamic observation data
CN112000893B (en) Resident region prediction method, device, equipment and storage medium
CN111080682A (en) Point cloud data registration method and device
Guo et al. Using VIIRS-DNB and landsat data for impervious surface area mapping in an arid/semiarid region
KR20140056828A (en) Apparatus, method and computer readable recording medium for analyzing a floating population using a user terminal
CN115630535B (en) Urban surface heat island strength dynamic quantification method and device and electronic equipment
Xue et al. Bias estimation and correction for triangle-based surface area calculations
Ruiz-Lendínez et al. A point-based methodology for the automatic positional accuracy assessment of geospatial databases
Wang et al. Edge-directed interpolation-based sub-pixel mapping
Ying et al. Terrestrial surface-area increment: the effects of topography, DEM resolution, and algorithm
CN113849531B (en) Query method and device
CN112669426B (en) Three-dimensional geographic information model rendering method and system based on generation countermeasure network
Li et al. Exploring the performance of spatio-temporal assimilation in an urban cellular automata model
Xue et al. Surface area calculation for DEM-based terrain model
Hyväluoma Reducing the grid orientation dependence of flow routing on square-grid digital elevation models
CN107133295B (en) Accelerated inverse address resolution method and device based on recorder management platform
Ge et al. Rough set-derived measures in image classification accuracy assessment
CN113655483B (en) Method, system, equipment and medium for constructing weather radar reflectivity jigsaw data set
Ma et al. Establishment of regional tropospheric delay model in Australia
CN117195016B (en) Sewage treatment mode determining method and device, computer equipment and storage 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
GR01 Patent grant
GR01 Patent grant