CN105930633A - Method for forecasting urban heat island effect - Google Patents
Method for forecasting urban heat island effect Download PDFInfo
- Publication number
- CN105930633A CN105930633A CN201610207314.8A CN201610207314A CN105930633A CN 105930633 A CN105930633 A CN 105930633A CN 201610207314 A CN201610207314 A CN 201610207314A CN 105930633 A CN105930633 A CN 105930633A
- Authority
- CN
- China
- Prior art keywords
- matrix
- heat island
- calculated
- data
- moment
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 55
- 230000000694 effects Effects 0.000 title claims abstract description 35
- 239000011159 matrix material Substances 0.000 claims abstract description 104
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 7
- 238000012546 transfer Methods 0.000 claims description 26
- 230000005855 radiation Effects 0.000 claims description 8
- 230000007704 transition Effects 0.000 abstract description 6
- 238000011160 research Methods 0.000 description 11
- 230000008901 benefit Effects 0.000 description 5
- 239000002689 soil Substances 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000000116 mitigating effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
- 239000002918 waste heat Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2219/00—Indexing scheme relating to application aspects of data processing equipment or methods
- G06F2219/10—Environmental application, e.g. waste reduction, pollution control, compliance with environmental legislation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A30/00—Adapting or protecting infrastructure or their operation
- Y02A30/60—Planning or developing urban green infrastructure
Landscapes
- Radiation Pyrometers (AREA)
Abstract
The invention discloses a method for forecasting urban heat island effect. The method includes the steps of performing temperature inversion for remote sensing image data using a single window algorithm to obtain urban surface temperature data; grading the urban surface temperature data according to mean parameters and standard deviation parameters; conducting statistics of heat island area data of different grades; and calculating an initial state matrix S(0) based on the heat island region data of different grades at the first moment, calculating a transition matrix P based on the heat island region data of different grades at the second moment and the initial state matrix S(0), and calculating the evolution state data of heat island regions in the subsequent required time based on the initial state matrix S(0) and the transition matrix P. The invention constructs the state transition matrix and the initial state matrix by utilizing the temperature inversion grading results, and forecasts the next evolution trend of the heat island based on Markov chain transition probabilities. The invention is high in forecasting accuracy and has a great practical significance.
Description
Technical field
The present invention relates to a kind of Remote Sensing Data Processing method, particularly relate to the prediction of a kind of urban heat land effect
Method.
Background technology
Tropical island effect refers to when urban development to certain scale, due to the change of city underlying surface character, air
Pollute and artificial waste heat discharge etc. makes city temperature apparently higher than suburb, forms showing of similar high temperature isolated island
As.Underlying surface refers to the earth surface directly contacted with air lower floor, it include landform, geology, soil,
River and vegetation etc., be one of key factor affecting weather.Have from literature record from tropical island effect,
Just there are many scientists that many big and small cities, the whole world are observed, find according to statistics, to 2012
Till Nian, having in the world in the city of more than 1000 different scales and there is tropical island effect, scope is throughout south
Different latitude position, the Northern Hemisphere, different altitude height areas.Nowadays, tropical island effect has become as one
City public hazards, healthy to social economy, urban ecology and resident works the mischief.
In order to deeply probe into the Forming Mechanism of tropical island effect, alleviate tropical island effect adverse effect,
Lot of domestic and foreign scholar has carried out the research for urban heat island, and urban heat land effect has become current city
City's weather and one of highly important research contents in environmental studies.According to basic data source difference, city
The research method of city's thermal environment effect can be divided into four classes: meteorological data method, observation method of layouting, numerical value
Simulation method and remote sensing method.Owing to remotely-sensed data is for the abundant expression forms of urban thermal environment effect, can
Respond well depending on changing, beneficially the dynamic monitoring of city thermal field and Inner construction analysis, became city in recent years
The main stream approach of city's thermal environment effect study.At present, remote sensing technology research heat island, three kinds of sides can be divided into again
Method:
One, monitoring method based on vegetation index.
Vegetation index can react the density state of the ground vegetation extracted, and is widely used in qualitative and fixed
In the research evaluating vegetation coverage and vigor of amount, city's heat island is had by substantial amounts of research display surface vegetation
There is mitigation.Additionally one of shortcoming of vegetation index be requirement height above sea level within the specific limits, if high
Path difference reaches more than 500m such as hilly country, and result will be far from each other, its two be to ground cover
Type has certain requirement, vegetation is not enough in the winter time, be difficult to extract in the case of, result of study will lack
Weary accuracy.
Its two, the radiation temperature value method of inversion.
The bright temperature of satellite that thermal infrared remote sensing detects, by remotely-sensed data inverting surface temperature and temperature three
There is substantial connection, utilize bright temperature or surface temperature that urban Heat Environment is analyzed research and there is its conjunction
Rationality.Owing to the method difference for the treatment of temperature can be divided into again: monitoring method based on bright temperature and based on ground temperature
Monitoring method.The former advantage is take into account the factor impact of radiating surface and air etc., and shortcoming is
Calculating process is complicated, and X factor impact is too many.In view of the hardware condition of present stage satellite itself, base
Inverting in surface temperature there be difficulties involved when, currently mainly by trying to achieve atmospheric parameter, emissivity
Inverting surface temperature is carried out etc. parameter.
Its three, detection method based on " thermal landscape ".
Utilize the viewpoint of landscape esthetics, under the support of GIS and remote sensing technology, for thermal environment spatial framework with
Process sets up a set of research method and evaluation index.The method by multiformity, separating degree, dominance, point
The many factors such as dimension, shape index composition evaluation index, the advantage of the method is can be quickly and easily
Obtain the key elements such as the land use pattern in spatial dimension, ground temperature, NDVI;Shortcoming is exactly Time Continuous
Deficiency in property.
In the above four kinds of method.The advantage of conventional method is the description to urban heat land effect and announcement
Simple and brief, but it is substantially the research of on point scale or pure concept, Points replacing surfaces, is difficult to true
With effectively extend to and on, the feature such as layout, internal structure thus for the flat of research urban heat island
There is the biggest difficulty." urban heat island " that although RS obtains and temperature " urban heat island " traditionally
There is marked difference, " urban heat island " that obtain such as RS is the strongest, spatial variations rate is maximum,
Antithesis, but the two distribution on space-time is similar to this " urban heat island " characterized with temperature, and
And RS can carry out surface temperature mensuration timely, objective, in large area, describe urban heat island time
Space division cloth has the advantage that conventional method is incomparable.Especially with the quickening of remote sensing quantification paces,
Increasing Remote Sensing Parameters can be obtained by remotely-sensed data inverting, this more remote sensing technology wide
A good basis has been established in general application.
Urban heat land effect is monitored, analysis of law, extraction factor of influence simultaneously, use a certain
The research that tropical island effect is predicted by kind forecast model is also little, calculates heat further based on forecast model
Island effect tend to be steady state time distribution characteristics the most less.And study urban heat land effect, not only need
Deeply to probe into current tropical island effect distribution characteristics and inherent law, predict tropical island effect based on current state
Next step regularity of distribution is to improving urban environment and city decision-making is the most significant.
Summary of the invention
For drawbacks described above, the invention provides the Forecasting Methodology of a kind of urban heat land effect.
The Forecasting Methodology of the urban heat land effect that the present invention provides, including:
Step 1, uses mono window algorithm to carry out temperature retrieval remote sensing image data, it is thus achieved that urban surface temperature
Degrees of data;
Step 2, carries out classification to urban ground temperature data according to Mean Parameters and standard deviation criteria;
Step 3, adds up the heat island area data of each rank, described heat island area data include position,
Area, account for total area ratio;
Step 4, is calculated original state square according to the heat island area data of each rank in the first moment
Battle array S (0), according to heat island area data and the described original state matrix S of each rank in the second moment
(0) it is calculated transfer matrix P, counts according to described original state matrix S (0) and transfer matrix P
Calculate the described heat island region evolving state data in the follow-up required moment.
Above-mentioned Forecasting Methodology can also have the following characteristics that
In step 1, temperature retrieval specifically includes:
It is calculated radiation brightness according to remotely-sensed data, is calculated ground brightness according to described radiation brightness
Temperature value;
Calculate normalized differential vegetation index, calculate vegetation coverage according to described normalized differential vegetation index, according to
Described vegetation coverage calculates earth's surface emissivity;
It is calculated described surface temperature number according to described ground brightness temperature value and described earth's surface emissivity
According to.
Above-mentioned Forecasting Methodology can also have the following characteristics that
Described step 2 carries out classification according to below table:
Above-mentioned Forecasting Methodology can also have the following characteristics that
Described transfer matrix P meets following condition: the value of each element is positive number, and each row vector
All elements sum be 1.
Above-mentioned Forecasting Methodology can also have the following characteristics that
The temperature data in the second moment according to described heat island region and described original state in described step 4
Matrix S (0) is calculated transfer matrix P and includes:
Described second moment is k1 with the discrete time difference in described first moment, according to described heat island region
The temperature data in the second moment is calculated state matrix S (k1) corresponding to the second moment, according to S (k1)=S
(0)Pk1It is calculated P;
Original state matrix S described in described step 4 (0) and transfer matrix P calculate described heat island region
Evolving state in the follow-up required moment includes:
The described required moment is k2 with the discrete time difference in described first moment, is calculated institute according to following formula
Take and carve corresponding state matrix S (k2): S (k2)=S (0) Pk2。
Above-mentioned Forecasting Methodology can also have the following characteristics that
Described method also includes: use iterative method pre-according to described original state matrix and described transfer matrix
Survey the heat island region of each rank evolving state under plateau.
Above-mentioned Forecasting Methodology can also have the following characteristics that
Described use iterative method predicts heat island region according to described original state matrix and described transfer matrix
Evolving state under plateau includes:
The difference calculating the transposed matrix P ' of described transfer matrix P and the E of unit matrix obtains matrix A, root
It is calculated matrix X according to AX=0 and is the matrix of heat island region evolving state under plateau.
Above-mentioned Forecasting Methodology can also have the following characteristics that
Described it be calculated matrix X according to AX=0 and include: iteration precision is set, arranges at the beginning of matrix X
Initial value, calculates the value of the X after iteration according to following formula, and after carrying out successive ignition calculating, obtain changes
The error of the respective element of the X before the value of each element of the X after Dai and iteration is respectively less than described iteration
During precision, using the X that finally gives as calculated matrix X;
Wherein, the line number of matrix A and columns are the value that n, i and j are between 0 to n, aijFor matrix
The element of the i-th row jth row in A,For the element value before iteration,For the element value after iteration.
Above-mentioned Forecasting Methodology can also have the following characteristics that
When the initial value of matrix X is set, it is set to random value, or below the meeting of matrix X is set
Condition: the value of each element is all higher than 0 and less than 1, and the element sum often gone is 1.
The present invention utilizes temperature retrieval classification results to construct state-transition matrix and original state matrix, root
Walk transition probability according to Markov Chain K, next step evolving trend of heat island is predicted, meanwhile, according to
The thought of the steady chain of Markov, the problem using the iterative algorithm infinite solution of solving matrix equation, and then solve
The certainly calculating of the steady chain of Markov.The prediction accuracy of the present invention is high, has bigger Practical significance.
Accompanying drawing explanation
Fig. 1 is the flow chart of the Forecasting Methodology of urban heat land effect in the present invention.
Specific embodiment
For making the purpose of the embodiment of the present invention, technical scheme and advantage clearer, below in conjunction with this
Accompanying drawing in bright embodiment, is clearly and completely described the technical scheme in the embodiment of the present invention,
Obviously, described embodiment is a part of embodiment of the present invention rather than whole embodiments.Based on
Embodiment in the present invention, those of ordinary skill in the art are obtained under not making creative work premise
The every other embodiment obtained, broadly falls into the scope of protection of the invention.It should be noted that do not conflicting
In the case of, the embodiment in the application and the feature in embodiment can mutual combination in any.
Fig. 1 is the flow chart of the Forecasting Methodology of urban heat land effect in the present invention.The method includes:
Step 1, uses mono window algorithm to carry out temperature retrieval remote sensing image data, it is thus achieved that urban surface temperature
Degrees of data;
Step 2, carries out classification to urban ground temperature data according to Mean Parameters and standard deviation criteria;
Step 3, adds up the heat island area data of each rank, heat island area data include position, area,
Account for total area ratio;
Step 4, is calculated original state square according to the heat island area data of each rank in the first moment
Battle array S (0), according to heat island area data and the original state matrix S of each rank in the second moment
(0) it is calculated transfer matrix P, calculates heat according to original state matrix S (0) and transfer matrix P
Region, island is in the evolving state data in follow-up required moment.
The following detailed description of this method:
In step 1, temperature retrieval specifically includes:
Step 101: be calculated radiation brightness according to remotely-sensed data, is calculated ground according to radiation brightness
Face brightness temperature value.
Concrete, according to Lλ=(Gain*DN)+offset following formula is calculated amount of radiation temperature Lλ, wherein, Gain
Representing gain coefficient, DN represents the size of each pixel, and offset represents deviation ratio.
Bright temperature T in ground is calculated according to following formula.
T is bright temperature (unit: DEG C), K1、K2For revising coefficient.
Step 102: calculate normalized differential vegetation index, calculates vegetation coverage according to normalized differential vegetation index,
Earth's surface emissivity is calculated according to vegetation coverage;
Normalized differential vegetation index NDVI is calculated according to following formula
V in formulaNIR、VRIt is respectively near infrared band, the pixel value of red spectral band.
According to following formula calculating vegetation coverage:
Wherein VNDVIRepresent NDVI value, i.e. normalized differential vegetation index, VNDVI_minAnd VNDVI_maxRepresent respectively
The maximum of NDVI and minima, can respectively be taken as 0.5 and 0.2.
Step 103: be calculated surface temperature data according to ground brightness temperature value and earth's surface emissivity.
Surface temperature Ts is calculated according to following formula
Wherein, λ represents the centre wavelength of Thermal infrared bands, ρ=h c/ σ 1.438 × 10‐2M K, wherein
Planck's constant h=6.626 × 10-34J s, Boltzmann constant σ=1.38 × 10-23J/K, the light velocity
C=2.998 × 108m/s。
Wherein, ε is earth's surface emissivity.
ε=εsoil, work as VNDVI< 0.2
ε=εveg, work as VNDVI> 0.5
ε=εsoilPV+ ε=εVEG(1-PV), as 0.2≤VNDVI≤0.5
Wherein, ε in formulasoilFor soil emissivity, εvegFor vegetation emissivity, PVFor vegetation coverage.
The empirical value that can also directly use Nichol to propose in 1994 in this method, i.e. has vegetation
Ground surface ε=0.95, unvegetated ground surface ε=0.92.
Described step 2 carries out classification according to below table:
Wherein, μ is Mean Parameters, and std is standard deviation criteria, and Ts is temperature parameter.
Transfer matrix P meets following condition: the value of each element is positive number, and the institute of each row vector
Having element sum is 1.
The temperature data in the second moment according to heat island region and original state matrix S (0) meter in step 4
Calculation obtains transfer matrix P and includes: the second moment was k1 with the discrete time difference in the first moment, according to heat island
The temperature data in second moment in region is calculated state matrix S (k1) corresponding to the second moment, according to
S (k1)=S (0) Pk1It is calculated P;
In step 4, original state matrix S (0) and transfer matrix P calculating heat island region is taken follow-up
The evolving state carved includes: the required moment is k2 with the discrete time difference in the first moment, calculates according to following formula
Obtain state matrix S (k2): S (k2)=S (0) P corresponding to required momentk2。
Method also includes: use iterative method to predict each rank according to original state matrix and transfer matrix
Heat island region evolving state under plateau.
Iterative method is used to predict that heat island region is under plateau according to original state matrix and transfer matrix
Evolving state include: the difference of the E calculating transposed matrix P ' and the unit matrix of transfer matrix P obtains square
Battle array A, is calculated matrix X according to AX=0 and is heat island region evolving state under plateau
Matrix.
It is calculated matrix X according to AX=0 to include: arrange iteration precision, the initial value of matrix X is set,
The value of the X after iteration is calculated, after carrying out successive ignition calculating, after the iteration obtained according to following formula
When the error of the respective element of the X before the value of each element of X and iteration is respectively less than iteration precision, will be
The X obtained eventually is as calculated matrix X;
Wherein, the line number of matrix A and columns are the value that n, i and j are between 0 to n, aijFor matrix
The element of the i-th row jth row in A,For the element value before iteration,For the element value after iteration.
When the initial value of matrix X is set, it is set to random value, or below the meeting of matrix X is set
Condition: the value of each element is all higher than 0 and less than 1, and the element sum often gone is 1.
Specific embodiment
Obtain the history remotely-sensed data of city A, use mono window algorithm to carry out temperature retrieval, it is thus achieved that city
Surface temperature data, carry out classification to urban ground temperature data according to Mean Parameters and standard deviation criteria,
It is divided into low-temperature space, secondary middle warm area, middle warm area, secondary high-temperature region, high-temperature region, extra-high warm area;Add up each
The heat island area data of rank, heat island area data includes position, area, accounts for total area ratio.According to
The heat island area data of each rank in October, 2009 is calculated original state matrix S (0), root
It is calculated current matrix S (K) according to the heat island area data of each rank in October, 2014.Root
Relation S (K)=S (0) P according to S (0) and S (K)K, obtain transfer matrix P.According to S (2K)
=S (0) P2KRelation be calculated the data in October, 2019.Calculate the transposition square of transfer matrix P
The difference of battle array P ' and the E of unit matrix obtains matrix A, is calculated matrix X according to AX=0 and is city
The matrix data of A evolving state under plateau.
The present invention utilizes temperature retrieval classification results to construct state-transition matrix and original state matrix, root
Walk transition probability according to Markov Chain K, next step evolving trend of heat island is predicted, meanwhile, according to
The thought of the steady chain of Markov, the problem using the iterative algorithm infinite solution of solving matrix equation, and then solve
The certainly calculating of the steady chain of Markov.The prediction accuracy of the present invention is high, has bigger Practical significance.
Descriptions above can combine enforcement individually or in every way, and these modification
Mode is all within protection scope of the present invention.
It should be noted that in this article, term " include ", " comprising " or its any other variant meaning
Containing comprising of nonexcludability, so that include that the article of a series of key element or equipment not only include
Those key elements, but also include other key elements being not expressly set out, or also include for this article
Or the key element that equipment is intrinsic.In the case of there is no more restriction, statement " including ... " limit
Key element, it is not excluded that in the article including key element or equipment, there is also other identical element.
Above example is only in order to illustrate technical scheme and unrestricted, reference only to preferably implementing
The present invention has been described in detail by example.It will be understood by those within the art that, can be to this
Bright technical scheme is modified or equivalent, without deviating from spirit and the model of technical solution of the present invention
Enclose, all should contain in the middle of scope of the presently claimed invention.
Claims (9)
1. the Forecasting Methodology of a urban heat land effect, it is characterised in that including:
Step 1, uses mono window algorithm to carry out temperature retrieval remote sensing image data, it is thus achieved that urban surface temperature
Degrees of data;
Step 2, carries out classification to urban ground temperature data according to Mean Parameters and standard deviation criteria;
Step 3, adds up the heat island area data of each rank, described heat island area data include position,
Area, account for total area ratio;
Step 4, is calculated original state square according to the heat island area data of each rank in the first moment
Battle array S (0), according to heat island area data and the described original state matrix S of each rank in the second moment
(0) it is calculated transfer matrix P, counts according to described original state matrix S (0) and transfer matrix P
Calculate the described heat island region evolving state data in the follow-up required moment.
2. the Forecasting Methodology of urban heat land effect as claimed in claim 1, it is characterised in that
In step 1, temperature retrieval specifically includes:
It is calculated radiation brightness according to remotely-sensed data, is calculated ground brightness according to described radiation brightness
Temperature value;
Calculate normalized differential vegetation index, calculate vegetation coverage according to described normalized differential vegetation index, according to
Described vegetation coverage calculates earth's surface emissivity;
It is calculated described surface temperature number according to described ground brightness temperature value and described earth's surface emissivity
According to.
3. the Forecasting Methodology of urban heat land effect as claimed in claim 1, it is characterised in that
Described step 2 carries out classification according to below table:
Wherein, μ is Mean Parameters, and std is standard deviation criteria, and Ts is temperature parameter.
4. the Forecasting Methodology of urban heat land effect as claimed in claim 1, it is characterised in that
Described transfer matrix P meets following condition: the value of each element is positive number, and each row vector
All elements sum be 1.
5. the Forecasting Methodology of urban heat land effect as claimed in claim 1, it is characterised in that
The temperature data in the second moment according to described heat island region and described original state in described step 4
Matrix S (0) is calculated transfer matrix P and includes:
Described second moment is k1 with the discrete time difference in described first moment, according to described heat island region
The temperature data in the second moment is calculated state matrix S (k1) corresponding to the second moment, according to S (k1)=S
(0)Pk1It is calculated P;
Original state matrix S described in described step 4 (0) and transfer matrix P calculate described heat island region
Evolving state in the follow-up required moment includes:
The described required moment is k2 with the discrete time difference in described first moment, is calculated institute according to following formula
Take and carve corresponding state matrix S (k2): S (k2)=S (0) Pk2。
6. the Forecasting Methodology of urban heat land effect as claimed in claim 1, it is characterised in that
Described method also includes: use iterative method pre-according to described original state matrix and described transfer matrix
Survey the heat island region of each rank evolving state under plateau.
7. the Forecasting Methodology of urban heat land effect as claimed in claim 6, it is characterised in that
Described use iterative method predicts heat island region according to described original state matrix and described transfer matrix
Evolving state under plateau includes:
The difference calculating the transposed matrix P ' of described transfer matrix P and the E of unit matrix obtains matrix A, root
It is calculated matrix X according to AX=0 and is the matrix of heat island region evolving state under plateau.
8. the Forecasting Methodology of urban heat land effect as claimed in claim 7, it is characterised in that
Described it be calculated matrix X according to AX=0 and include: iteration precision is set, arranges at the beginning of matrix X
Initial value, calculates the value of the X after iteration according to following formula, and after carrying out successive ignition calculating, obtain changes
The error of the respective element of the X before the value of each element of the X after Dai and iteration is respectively less than described iteration
During precision, using the X that finally gives as calculated matrix X;
Wherein, the line number of matrix A and columns are the value that n, i and j are between 0 to n, aijFor matrix
The element of the i-th row jth row in A,For the element value before iteration,For the element value after iteration.
9. the Forecasting Methodology of urban heat land effect as claimed in claim 8, it is characterised in that
When the initial value of matrix X is set, it is set to random value, or below the meeting of matrix X is set
Condition: the value of each element is all higher than 0 and less than 1, and the element sum often gone is 1.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610207314.8A CN105930633B (en) | 2016-04-05 | 2016-04-05 | A kind of prediction technique of urban heat land effect |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610207314.8A CN105930633B (en) | 2016-04-05 | 2016-04-05 | A kind of prediction technique of urban heat land effect |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105930633A true CN105930633A (en) | 2016-09-07 |
CN105930633B CN105930633B (en) | 2018-08-17 |
Family
ID=56840538
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610207314.8A Expired - Fee Related CN105930633B (en) | 2016-04-05 | 2016-04-05 | A kind of prediction technique of urban heat land effect |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105930633B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107391910A (en) * | 2017-07-03 | 2017-11-24 | 国网湖南省电力公司 | Computational methods and its system of a kind of thermal power plant's hot driving to mixed layer atmospheric heating |
CN109211791A (en) * | 2018-08-10 | 2019-01-15 | 北京观微科技有限公司 | Crop condition monitoring method and system |
CN109612587A (en) * | 2018-12-18 | 2019-04-12 | 广州大学 | A kind of urban Heat Environment cause diagnosis method and system |
CN110208878A (en) * | 2019-06-14 | 2019-09-06 | 广西海佩智能科技有限公司 | Green Roof weather monitoring and tropical island effect impact evaluation method |
CN112200349A (en) * | 2020-09-16 | 2021-01-08 | 平衡机器科技(深圳)有限公司 | Remote sensing image heat island effect prediction method based on single window algorithm and PredRNN |
-
2016
- 2016-04-05 CN CN201610207314.8A patent/CN105930633B/en not_active Expired - Fee Related
Non-Patent Citations (9)
Title |
---|
DAVID R. STREUTKER: ""A remote sensing study of the urban heat island of Houston, Texas"", 《INTERNATIONAL JOURNAL OF REMOTE SENSING》 * |
ZHU DONGJIN: ""The stationary distribution of markov chains in random environments"", 《CHINESE JOURNAL OF APPLIED PROBABILITY AND STATISTICS》 * |
刘文渊等: ""不同地表参数变化的上海市热岛效应时空分析"", 《遥感技术与应用》 * |
夏安邦: "《系统建模理论与方法》", 30 September 2008, 北京:机械工业出版社 * |
张殿江等: ""基于马尔科夫链模型的城市热岛扩散趋势预测——以天津滨海新区为例"", 《中国人口·资源与环境》 * |
王文杰等: ""基于遥感的北京市城市化发展与城市热岛效应变化关系研究"", 《环境科学研究》 * |
范新岗等: ""气候系统可预报性的全局研究"", 《气象学报》 * |
许辉熙: ""成都平原中等城市的热岛效应动态特征对比研究"", 《测绘与空间地理信息》 * |
许辉熙: ""资阳市城市热岛效应动态评价"", 《四川建筑》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107391910A (en) * | 2017-07-03 | 2017-11-24 | 国网湖南省电力公司 | Computational methods and its system of a kind of thermal power plant's hot driving to mixed layer atmospheric heating |
CN109211791A (en) * | 2018-08-10 | 2019-01-15 | 北京观微科技有限公司 | Crop condition monitoring method and system |
CN109612587A (en) * | 2018-12-18 | 2019-04-12 | 广州大学 | A kind of urban Heat Environment cause diagnosis method and system |
CN110208878A (en) * | 2019-06-14 | 2019-09-06 | 广西海佩智能科技有限公司 | Green Roof weather monitoring and tropical island effect impact evaluation method |
CN112200349A (en) * | 2020-09-16 | 2021-01-08 | 平衡机器科技(深圳)有限公司 | Remote sensing image heat island effect prediction method based on single window algorithm and PredRNN |
CN112200349B (en) * | 2020-09-16 | 2022-06-03 | 平衡机器科技(深圳)有限公司 | Remote sensing image heat island effect prediction method based on single window algorithm and PredRNN |
Also Published As
Publication number | Publication date |
---|---|
CN105930633B (en) | 2018-08-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li et al. | Quantifying 3D building form effects on urban land surface temperature and modeling seasonal correlation patterns | |
Zhang et al. | An urban heat island study in Nanchang City, China based on land surface temperature and social-ecological variables | |
CN105930633A (en) | Method for forecasting urban heat island effect | |
Zeng et al. | Daily global solar radiation in China estimated from high‐density meteorological observations: a random forest model framework | |
CN105243435B (en) | A kind of soil moisture content prediction technique based on deep learning cellular Automation Model | |
CN109635309A (en) | A kind of surface temperature space NO emissions reduction method | |
CN107678075A (en) | A kind of urban heat land effect monitoring method and system based on domestic satellite | |
CN102507586B (en) | Remote sensing monitoring method for carbon emission | |
Mishra et al. | Spatial representation of organic carbon and active-layer thickness of high latitude soils in CMIP5 earth system models | |
Wang et al. | Spatial distribution and influencing factors on urban land surface temperature of twelve megacities in China from 2000 to 2017 | |
Dai et al. | Thermal impacts of greenery, water, and impervious structures in Beijing’s Olympic area: A spatial regression approach | |
Yang et al. | The impact analysis of water body landscape pattern on urban heat island: a case study of Wuhan City | |
Shi et al. | Influence of landscape features on urban land surface temperature: Scale and neighborhood effects | |
Han et al. | Seasonal effects of urban morphology on land surface temperature in a three-dimensional perspective: A case study in Hangzhou, China | |
Xie et al. | A backpropagation neural network improved by a genetic algorithm for predicting the mean radiant temperature around buildings within the long-term period of the near future | |
Wu et al. | Nonlinear forces in urban thermal environment using Bayesian optimization-based ensemble learning | |
Jiang et al. | A Municipal PM2. 5 Forecasting Method Based on Random Forest and WRF Model. | |
CN114005048A (en) | Multi-temporal data-based land cover change and thermal environment influence research method | |
Chen et al. | Integrating weather observations and local-climate-zone-based landscape patterns for regional hourly air temperature mapping using machine learning | |
Chen et al. | Quantitative analysis of the building-level relationship between building form and land surface temperature using airborne LiDAR and thermal infrared data | |
Lin et al. | Exploring the relationship between thermal environmental factors and land surface temperature of a “furnace city” based on local climate zones | |
Wu et al. | Simulating spatiotemporal land use change in middle and high latitude regions using multiscale fusion and cellular automata: The case of Northeast China | |
Pouyan et al. | Spatial and seasonal modeling of the land surface temperature using random forest | |
Chang et al. | Impact of refined land surface properties on the simulation of a heavy convective rainfall process in the Pearl River Delta region, China | |
Yu et al. | Downscaling mapping method for local climate zones from the perspective of deep learning |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180817 |
|
CF01 | Termination of patent right due to non-payment of annual fee |