CN111637972A - Remote sensing definition method for industrial heat island effect - Google Patents

Remote sensing definition method for industrial heat island effect Download PDF

Info

Publication number
CN111637972A
CN111637972A CN202010331105.0A CN202010331105A CN111637972A CN 111637972 A CN111637972 A CN 111637972A CN 202010331105 A CN202010331105 A CN 202010331105A CN 111637972 A CN111637972 A CN 111637972A
Authority
CN
China
Prior art keywords
surface temperature
industrial
area
response
layer
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.)
Granted
Application number
CN202010331105.0A
Other languages
Chinese (zh)
Other versions
CN111637972B (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.)
Aerospace Information Research Institute of CAS
Original Assignee
Aerospace Information Research Institute 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 Aerospace Information Research Institute of CAS filed Critical Aerospace Information Research Institute of CAS
Priority to CN202010331105.0A priority Critical patent/CN111637972B/en
Publication of CN111637972A publication Critical patent/CN111637972A/en
Application granted granted Critical
Publication of CN111637972B publication Critical patent/CN111637972B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • 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/48Thermography; Techniques using wholly visual means
    • G01J5/485Temperature profile
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • 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)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Algebra (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Optimization (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)
  • Radiation Pyrometers (AREA)

Abstract

The heat pollution generated by industrial activities has serious negative influence on urban ecological environment, but a quantitative measuring and calculating method and a flow for industrial area heat environment effect are lacked at present. The invention aims at the problems and discloses a remote sensing definition method of an industrial heat island effect, which comprises the following steps: step 1) framing the vector boundary of an industrial target area by utilizing Google Earth; step 2) carrying out data preprocessing on the Landsat8 image of the target industrial area and inverting the earth surface temperature; step 3) constructing 5 kilometer buffer areas on the periphery of an industrial area, and dividing 100 kilometer buffer areas equally to generate 100 layer-by-layer buffer areas of 50 meters; step 4) calculating the average earth surface temperature of each layer of buffer area, drawing a scatter diagram of the average earth surface temperature along with the change of the distance, and performing interpolation by using a B-spline curve; and 5) finding a first change critical point according to the surface temperature change curve, and further calculating a response range (km), a response temperature difference (DEG C) and response intensity (DEG C/km).

Description

Remote sensing definition method for industrial heat island effect
Technical Field
The invention relates to a remote sensing definition method of an industrial heat island effect, which provides a remote sensing measurement method of the heat island effect suitable for a heat emission ground object target aiming at the problem of the heat environment effect of an industrial park in a city.
Background
Worldwide, rapid urbanization has led to significant changes in urban climate and surface biochemical processes. Among them, the urban heat island effect (UHI effect) is one of the most typical urban ecological environmental problems, i.e., the phenomenon that the air temperature and the surface temperature (LST) of an urban area are higher than those of suburbs around the urban area. The rise of the urban temperature affects the habitability of urban residents, aggravates air pollution, increases the consumption of cooling energy and damages human health. In response to the above problems, urban researchers are eagerly evaluating strategies that can mitigate further warming of urban areas.
Current research has extensively summarized that human waste heat from industrial, transportation and building energy sources is one of the major causes of urban heat island effects (Kato and Yamaguchi 2005, Kotharkar and surarwar 2016, Papadopoulos and Moussiopoulos 2004, Han and Taylor 2015). In the course of resource development and utilization, the basic materials and energy consumption required for industrial production development raise a series of environmental problems (Rao et al 2018), causing redistribution of urban heat and fluctuation of surface energy (Chakraborty, Kant and Mitra 2015). According to previous studies, industrial sites play a dominant role in the formation of high surface temperatures (Li et al, 2011), and urban high-temperature regions are highly concentrated in industrial areas (Pearsall 2017, Tran et al 2017). Furthermore, if the influence of industrial production activities on the environment is neglected, the investment of manpower, material resources, and financial resources for alleviating environmental deterioration is even doubled (Wan, Shenand Choi 2017). In addition, artificial surfaces (i.e., man-made areas), such as commercial, residential and industrial areas, often study the changes in urban thermal environment as a whole (Li et al 2011, Chakraborty et al 2015). Therefore, the industrial park, as a separate research target, is dedicated to the study of the relationship between the thermal response and the thermal emissions of the industrial landscape target, and is yet to be deeply explored.
Previous studies have also focused on the impact of urban climate control targets on thermal environments, such as the cold island effect in urban greenhouses (large parks), urban waters, etc. (Yan, Wu and Dong 2018, Cao et al 2010, Bowler et al 2010, Bartesaghi-Koc, Osmond and Peters 2019). Early studies estimated the cold island effect range of influence of urban waters based on measured data of fixed point locations (Chang, Li and Chang 2007). In order to adopt a more convenient data acquisition mode, (Du et al.2016) utilizes Google Earth and landsat-8 satellite image data to quantitatively detect the cold island effect of the water body in Shanghai city. In summary, a large number of studies prove that dynamic thermal radiation information can be accurately extracted in real time by using a remote sensing technology, and a powerful support can be provided for regional thermal environment optimization (Wilson et al 2003, Stone and Rodgers2001, Meng et al 2018, Dos Santos et al 2017, Powers et al 2015, Feng et al 2019).
In summary, accurately quantifying the influence of the industrial park on the thermal environment is an urgent task of industrial thermal pollution control at present. Therefore, the invention focuses on the research in the large-scale industrial park in the city, and aims to establish a quantitative measuring and calculating method for the industrial heat island effect. Environmental departments and city planners can adjust the deployment of production tasks according to the intensity of industrial heat island effect in different seasons, and aim to reduce the adverse effect of industrial production activities on local geothermal environment.
Disclosure of Invention
The invention aims to provide a remote sensing definition method of an industrial heat island effect, which aims to solve the technical defect problem of the industrial heat island effect in the existing research and utilizes a thermal infrared remote sensing image to quantitatively measure and calculate the local geothermal environment effect of an urban industrial area.
The purpose of the invention is realized by the following technical steps:
step 1) framing the vector boundary of an industrial target area by utilizing Google Earth;
step 2) carrying out data preprocessing on the Landsat8 image of the target industrial area and inverting the earth surface temperature;
step 3) constructing 5 kilometer buffer areas on the periphery of an industrial area, and dividing 100 kilometer buffer areas equally to generate 100 layer-by-layer buffer areas of 50 meters;
step 4) calculating the average earth surface temperature of each layer of buffer area, drawing a scatter diagram of the average earth surface temperature along with the change of the distance, and performing interpolation by using a B-spline curve;
and 5) finding a first change critical point according to the surface temperature change curve, and further calculating a response range (km), a response temperature difference (DEG C) and response intensity (DEG C/km).
Drawings
FIG. 1 is an analytic graph based on the average surface temperature profile of a continuous buffer zone in an industrial zone;
FIG. 2 is a graph of inversion results of surface temperature of a sample steel plant;
FIG. 3 is a graph of the fitting results of the average surface temperature variation curves of the sample iron and steel plant;
Detailed Description
The invention 'a remote sensing definition method of industrial heat island effect' is further explained with the accompanying drawings.
(one) study region vector bounding box determination
Firstly, the basic position of the target area is determined by visual interpretation by using Google Earth, the vector boundary of the industrial target area is framed by using an 'adding polygon' tool and is stored in a kml or kmz format, and then vector data is converted into a shp format by using an 'rolling out from kml' tool of ArcGIS and is stored again.
(II) data preprocessing and surface temperature remote sensing inversion
Firstly, preprocessing of radiometric calibration and atmospheric correction is carried out on an original Landsat8 primary data product, and then the earth surface temperature is inverted by using Band10 data of Landsat8 based on a radiation transmission equation, wherein the radiation transmission equation is as follows:
Figure BDA0002464133200000031
in the above formula, TsRepresenting the surface temperature to be calculated; l issensorThe radiation brightness value after radiation calibration; b (T)s) Is at a temperature TsBlack body radiation value of time;
Figure BDA0002464133200000032
and
Figure BDA0002464133200000033
respectively refers to the uplink radiation and the downlink radiation of the atmosphere; τ is the atmospheric transmittance; the data preprocessing process can be realized by an ENVI radiation correction tool and an atmospheric correction tool, the calculation process of surface temperature inversion can be realized by an ENVI wave band operation tool, and the surface temperature inversion result of the sample area is shown in FIG. 2.
(III) layer-by-layer buffer construction
Constructing buffer areas of 5 kilometers at the periphery of an industrial area, and dividing 100 buffer areas equally to generate 100 buffer areas of 50 meters layer by layer; this step can be implemented in a batch process by the ArcGIS buffer tool.
(IV) mean surface temperature variation curve fitting
Calculating the average earth surface temperature of each layer of buffer area, drawing a scatter diagram of the average earth surface temperature changing along with the distance, and fitting by using a B-spline curve to obtain a smoothed earth surface temperature change curve (the abscissa is the distance from the boundary of the industrial area, and the ordinate is the earth surface temperature); the average surface temperature of the layer-by-layer buffer area can be obtained by using the attribute function of ArcGIS, and the generation of the surface temperature change curve is realized by using a nonlinear curve fitting tool after the ArcGIS is introduced into origin.
(V) quantitative index calculation
According to the surface temperature change curve (figure 1), a first turning critical point is found, and then a response range (km), a response temperature difference (DEG C) and response intensity (DEG C/km) are calculated. The turning critical point refers to the first point position with non-negative curve slope, and the abscissa of the point position is the response range (km); the ordinate of the turning critical point represents the surface temperature of the turning point, and the difference between the boundary surface temperature of the industrial area and the surface temperature of the turning point is calculated, so that the response temperature difference (DEG C) can be calculated; the response intensity (DEG C/km) is the ratio of the response temperature difference to the response range. The results of the sample zone mean surface temperature curve fit are shown in fig. 3.

Claims (6)

1. A remote sensing definition method of industrial heat island effect comprises the following steps:
step 1) framing the vector boundary of an industrial target area by utilizing Google Earth;
step 2) carrying out data preprocessing on the Landsat8 image of the target industrial area and inverting the earth surface temperature;
step 3) constructing 5 kilometer buffer areas on the periphery of an industrial area, and dividing 100 kilometer buffer areas equally to generate 100 layer-by-layer buffer areas of 50 meters;
step 4) calculating the average earth surface temperature of each layer of buffer area, drawing a scatter diagram of the average earth surface temperature along with the change of the distance, and performing interpolation by using a B-spline curve;
and 5) finding a first change critical point according to the surface temperature change curve, and further calculating a response range (km), a response temperature difference (DEG C) and response intensity (DEG C/km).
2. The method of claim 1, wherein step 1): the basic position of the target area is determined by visual interpretation by using Google Earth, the vector boundary of the industrial target area is framed by using an 'adding polygon' tool and is stored in a kml or kmz format, and then vector data is converted into a shp format by using an 'rolling out from kml' tool of ArcGIS.
3. The method of claim 1, wherein step 2): firstly, preprocessing radiation correction and atmospheric correction is carried out on an original Landsat8 primary data product, and then the earth surface temperature is inverted by using Band10 data of Landsat8 based on a radiation transmission equation, wherein the radiation transmission equation is as follows:
Figure FDA0002464133190000011
in the above formula, TsRepresenting the surface temperature to be calculated; l issensorThe radiation brightness value after radiation calibration; b (T)s) Is at a temperature TsBlack body radiation value of time;
Figure FDA0002464133190000012
and
Figure FDA0002464133190000013
respectively refers to the uplink radiation and the downlink radiation of the atmosphere; τ is the atmospheric transmittance; the data preprocessing process can be realized by an ENVI radiation calibration tool and an atmospheric correction tool, and the earth surface temperature inversion calculation process can be realized by an ENVI wave band operation tool.
4. The method of claim 1, wherein step 3): constructing buffer areas of 5 kilometers at the periphery of an industrial area, and dividing 100 buffer areas equally to generate 100 buffer areas of 50 meters layer by layer; this step can be implemented by the buffer tool of ArcGIS.
5. The method of claim 1, wherein step 4): calculating the average earth surface temperature of each layer of buffer area, drawing a scatter diagram of the average earth surface temperature changing along with the distance, and fitting by using a B-spline curve to obtain a smoothed earth surface temperature change curve (the abscissa is the distance from the boundary of the industrial area, and the ordinate is the earth surface temperature); the average surface temperature of the layer-by-layer buffer area can be obtained by using the 'attribute' function of ArcGIS, and the generation of the surface temperature change curve is realized by using a nonlinear curve fitting tool after the ArcGIS is introduced into origin.
6. The method of claim 1, wherein step 5): and finding a first turning critical point according to the surface temperature change curve, and further calculating a response range (km), a response temperature difference (DEG C) and response intensity (DEG C/km). The turning critical point refers to the first point position with non-negative curve slope, and the abscissa of the point position is the response range (km); the ordinate of the turning critical point represents the surface temperature of the turning point, and the difference between the boundary surface temperature of the industrial area and the surface temperature of the turning point is calculated, so that the response temperature difference (DEG C) can be calculated; the response intensity (DEG C/km) is the ratio of the response temperature difference to the response range.
CN202010331105.0A 2020-04-23 2020-04-23 Remote sensing definition method for industrial heat island effect Active CN111637972B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010331105.0A CN111637972B (en) 2020-04-23 2020-04-23 Remote sensing definition method for industrial heat island effect

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010331105.0A CN111637972B (en) 2020-04-23 2020-04-23 Remote sensing definition method for industrial heat island effect

Publications (2)

Publication Number Publication Date
CN111637972A true CN111637972A (en) 2020-09-08
CN111637972B CN111637972B (en) 2021-07-13

Family

ID=72328711

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010331105.0A Active CN111637972B (en) 2020-04-23 2020-04-23 Remote sensing definition method for industrial heat island effect

Country Status (1)

Country Link
CN (1) CN111637972B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115630535A (en) * 2022-12-19 2023-01-20 中国农业科学院农业资源与农业区划研究所 Dynamic quantification method and device for urban surface heat island strength and electronic equipment

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100305913A1 (en) * 2009-05-29 2010-12-02 Johnson Daniel P Method of modeling the socio-spatial dynamics of extreme urban heat events
CN102235911A (en) * 2010-04-20 2011-11-09 北京师范大学 Method for detecting thermal radiation directivity of urban goal
CN102944880A (en) * 2012-11-13 2013-02-27 中国科学院城市环境研究所 Urban industrial thermal pollution source investigation method based on remote sensing
CN103617565A (en) * 2013-12-06 2014-03-05 国网山东省电力公司 Method for establishing relationship between city micro weather and electric power air conditioner loads
CN105243266A (en) * 2015-09-18 2016-01-13 中国科学院遥感与数字地球研究所 Heat island intensity defining and spatial feature remote sensing analysis method based on variable threshold
CN108332859A (en) * 2018-01-18 2018-07-27 广州大学 A kind of extracting method and device of urban heat island range
CN108549858A (en) * 2018-04-08 2018-09-18 武汉理工大学 A kind of quantitative evaluation method of urban heat land effect
CN108776360A (en) * 2018-06-13 2018-11-09 华南农业大学 A kind of method of urban heat island strength Monitoring on Dynamic Change

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100305913A1 (en) * 2009-05-29 2010-12-02 Johnson Daniel P Method of modeling the socio-spatial dynamics of extreme urban heat events
CN102235911A (en) * 2010-04-20 2011-11-09 北京师范大学 Method for detecting thermal radiation directivity of urban goal
CN102944880A (en) * 2012-11-13 2013-02-27 中国科学院城市环境研究所 Urban industrial thermal pollution source investigation method based on remote sensing
CN103617565A (en) * 2013-12-06 2014-03-05 国网山东省电力公司 Method for establishing relationship between city micro weather and electric power air conditioner loads
CN105243266A (en) * 2015-09-18 2016-01-13 中国科学院遥感与数字地球研究所 Heat island intensity defining and spatial feature remote sensing analysis method based on variable threshold
CN108332859A (en) * 2018-01-18 2018-07-27 广州大学 A kind of extracting method and device of urban heat island range
CN108549858A (en) * 2018-04-08 2018-09-18 武汉理工大学 A kind of quantitative evaluation method of urban heat land effect
CN108776360A (en) * 2018-06-13 2018-11-09 华南农业大学 A kind of method of urban heat island strength Monitoring on Dynamic Change

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
SHI, YUAN: "Modelling the fine-scale spatiotemporal pattern of urban heat island effect using land use regression approach in a megacity", 《SCIENCE OF THE TOTAL ENVIRONMEN》 *
ZHANG, YU: "Research on the Contribution of Urban Land Surface Moisture to the Alleviation Effect of Urban Land Surface Heat Based on Landsat 8 Data", 《REMOTE SENSING》 *
刘勇洪: "京津冀城市群热岛定量评估", 《生态学报》 *
吴慧凤: "福州市城市热岛效应的转移矩阵与缓冲区", 《天津师范大学学报(自然科学版)》 *
胡远满: "RS与GIS在城市热岛效应研究中的应用", 《环境保护科学》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115630535A (en) * 2022-12-19 2023-01-20 中国农业科学院农业资源与农业区划研究所 Dynamic quantification method and device for urban surface heat island strength and electronic equipment

Also Published As

Publication number Publication date
CN111637972B (en) 2021-07-13

Similar Documents

Publication Publication Date Title
Georgakis et al. Studying the effect of “cool” coatings in street urban canyons and its potential as a heat island mitigation technique
Wu et al. Seasonal variation in the effects of urban environmental factors on land surface temperature in a winter city
Gao et al. Quantifying the nonlinear relationship between block morphology and the surrounding thermal environment using random forest method
Abbes et al. Wind resource estimation and wind park design in El-Kef region, Tunisia
Fu et al. Exploring the multitemporal surface urban heat island effect and its driving relation in the Beijing-Tianjin-Hebei urban agglomeration
Huang et al. Sensitivity of urban boundary layer simulation to urban canopy models and PBL schemes in Beijing
Huang et al. CFD analysis on traffic-induced air pollutant dispersion under non-isothermal condition in a complex urban area in winter
CN111637972B (en) Remote sensing definition method for industrial heat island effect
CN101936777A (en) Method for inversing air temperature of surface layer based on thermal infrared remote sensing
Lin et al. Dynamic influencing mechanism of traditional settlements experiencing urbanization: A case study of Chengzi Village
Veisi et al. Analysis of solar radiation towards optimization and location of the urban blocks in the neighborhood units
Liu et al. Impacts of high-albedo urban surfaces on outdoor thermal environment across morphological contexts: A case of Tianjin, China
Li et al. Modeling the Impact of Urban Three-Dimensional Expansion on Atmospheric Environmental Conditions in an Old Industrial District: A Case Study in Shenyang, China.
Yang et al. Transforming and validating urban microclimate data with multi-sourced microclimate datasets for building energy modelling at urban scale
Zubair et al. Assessment of photovoltaic capabilities in urban environments: A case study of Islamabad, Pakistan
Huo et al. Simulation of urban functional zone air temperature based on urban weather generator (UWG): a case study of Beijing, China
Xie et al. Traffic emission transportation in street canyons
Yuan et al. Urban sprawl and thermal response in the guangdong-hong kong-macao greater bay area, 2000–2020
CN107490363A (en) A kind of heat production enterprise production intensity monitoring method based on surface temperature data
Listyawati et al. Analysis of Urban Heat Island Phenomenon as A Global Warming Control Based on Remote Sensing in Jember Urban, Indonesia
Guan et al. A Quantitative Analysis of the Complex Response Relationship between Urban Green Infrastructure (Ugi) Structure/Spatial Pattern and Urban Thermal Environment in Shanghai
Park et al. Detecting and classifying rooftops with a CNN-based remote-sensing method for urban area cool roof application
Cao et al. Towards urban wind utilization: The spatial characteristics of wind energy in urban areas
Dave et al. Modelling the impact of urban heat island mitigation strategies on urban air quality
CN112464147A (en) Thermal landscape pattern measuring method based on thermal anomaly classification

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