CN107247690B - Estimate the method and service terminal of temperature - Google Patents
Estimate the method and service terminal of temperature Download PDFInfo
- Publication number
- CN107247690B CN107247690B CN201710434299.5A CN201710434299A CN107247690B CN 107247690 B CN107247690 B CN 107247690B CN 201710434299 A CN201710434299 A CN 201710434299A CN 107247690 B CN107247690 B CN 107247690B
- Authority
- CN
- China
- Prior art keywords
- air temperature
- value
- altitude
- temperature value
- temperature
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- 230000009467 reduction Effects 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 14
- 238000004088 simulation Methods 0.000 abstract 2
- 238000010586 diagram Methods 0.000 description 8
- 230000006870 function Effects 0.000 description 5
- 238000004422 calculation algorithm Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 238000004590 computer program Methods 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 229910001316 Ag alloy Inorganic materials 0.000 description 1
- BQCADISMDOOEFD-UHFFFAOYSA-N Silver Chemical compound [Ag] BQCADISMDOOEFD-UHFFFAOYSA-N 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000002844 melting Methods 0.000 description 1
- 230000008018 melting Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 238000010792 warming Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01W—METEOROLOGY
- G01W1/00—Meteorology
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01W—METEOROLOGY
- G01W1/00—Meteorology
- G01W2001/006—Main server receiving weather information from several sub-stations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Environmental & Geological Engineering (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Environmental Sciences (AREA)
- Algebra (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Ecology (AREA)
- Evolutionary Biology (AREA)
- Biodiversity & Conservation Biology (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Atmospheric Sciences (AREA)
- Radiation Pyrometers (AREA)
Abstract
The present invention relates to meteorological technical fields, more particularly to a kind of method and service terminal for estimating temperature, the method of the estimation temperature is applied to service terminal, and the method for the estimation temperature includes the temperature lapse rate of pixel where calculating according to satellite remote sensing inverting surface temperature and digital elevation model.According to the temperature lapse rate of meteorological site pixel where the temperature value and the meteorological site that practical height above sea level measures, the meteorological site is calculated in the first estimation temperature value of the predetermined height above sea level.Establish regression model, the meteorological site is obtained in the simulation temperature value of the predetermined height above sea level, according to the meteorological site in the first estimation temperature value and simulation temperature value of the predetermined height above sea level, the meteorological site is calculated in the second estimation temperature value of the predetermined height above sea level.The accurate estimation to temperature is realized by this programme, is avoided and is directly acquired temperature record bring error.
Description
Technical Field
The invention relates to the technical field of weather, in particular to a method for estimating air temperature and a service terminal.
Background
The air temperature regulation not only regulates a plurality of near-surface processes, but also influences various aspects such as snow melting in glaciers, frozen soil degradation, global climate warming and the like, and is an important near-surface meteorological parameter in models of various plants, hydrology, weather, environment and the like. At present, the spatial air temperature data is obtained indirectly mainly through spatial interpolation and remote sensing data of meteorological sites. The spatial interpolation is easily influenced by terrain, and the error is large; the remote sensing data adopts near infrared data to invert the air temperature with great difficulty, and the air temperature is estimated mainly by adopting a multi-parameter statistical relationship and a temperature vegetation index method. Therefore, it is necessary to provide an accurate air temperature estimation method.
Disclosure of Invention
The invention aims to provide a method for estimating air temperature, which is used for realizing accurate estimation of the air temperature and reducing the error of directly measuring air temperature data.
Another object of the present invention is to provide a service terminal to achieve accurate estimation of air temperature and reduce errors in directly measuring air temperature data.
In order to achieve the above purpose, the embodiment of the present invention adopts the following technical solutions:
in a first aspect, an embodiment of the present invention provides a method for estimating an air temperature, which is applied to a service terminal, and is configured to combine and analyze earth surface temperature data inverted on a satellite remote sensing image and an air temperature measured at a meteorological site to accurately estimate an air temperature, where the method includes:
calculating a first estimated air temperature value of the meteorological station at the preset altitude according to the air temperature value measured by the meteorological station at the actual altitude and the temperature reduction rate of the pixel of the meteorological station on the satellite remote sensing image;
establishing a regression model to obtain a simulated air temperature value of the meteorological station at the preset altitude;
and calculating a second estimated air temperature value of the meteorological station at the preset altitude according to the first estimated air temperature value and the simulated air temperature value of the meteorological station at the preset altitude.
In a second aspect, an embodiment of the present invention further provides a service terminal, where the service terminal includes:
a memory;
a processor; and a means for estimating air temperature, said means for estimating air temperature being installed in said memory and comprising one or more software functional modules executed by said processor, said means for estimating air temperature comprising:
the first calculation module is used for calculating a first estimated air temperature value of the meteorological station at the preset altitude according to the air temperature value measured by the meteorological station at the actual altitude and the temperature direct reduction rate of the pixel of the meteorological station on the satellite remote sensing image;
the model establishing module is used for establishing a regression model to obtain a simulated air temperature value of the meteorological station at the preset altitude;
and the second calculation module is used for calculating a second estimated air temperature value of the meteorological station at the preset altitude according to the first estimated air temperature value and the simulated air temperature value of the meteorological station at the preset altitude.
The method for estimating the air temperature is applied to the service terminal and comprises the step of calculating a first estimated air temperature value of the meteorological site at the preset altitude according to the air temperature value measured by the meteorological site at the actual altitude and the temperature reduction rate of the meteorological site. Establishing a regression model to obtain a simulated air temperature value of the meteorological site at the preset altitude, and calculating a second estimated air temperature value of the meteorological site at the preset altitude according to a first estimated air temperature value and the simulated air temperature value of the meteorological site at the preset altitude. According to the scheme, the temperature is accurately estimated, and errors caused by directly acquiring temperature data are avoided.
In order to make the aforementioned and other objects, features and advantages of the present invention comprehensible, preferred embodiments accompanied with figures are described in detail below.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings needed to be used in the embodiments will be briefly described below, it should be understood that the following drawings only illustrate some embodiments of the present invention and therefore should not be considered as limiting the scope, and for those skilled in the art, other related drawings can be obtained according to the drawings without inventive efforts.
Fig. 1 shows a schematic structural diagram of a service terminal according to an embodiment of the present invention.
Fig. 2 is a flow chart illustrating a method for estimating air temperature according to an embodiment of the present invention.
Fig. 3 is a flow chart illustrating sub-steps of a method for estimating an air temperature according to an embodiment of the present invention.
Fig. 4 is a flow chart illustrating another method for estimating air temperature according to an embodiment of the present invention.
Fig. 5 is a schematic diagram illustrating functional blocks of an apparatus for estimating an air temperature according to an embodiment of the present invention.
The figure is as follows: 300-a service terminal; 310-means for estimating air temperature; 320-a memory; 330-a processor; 311-a first calculation module; 312-model building module; 313-a second calculation module.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. The components of 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. Thus, the following detailed description of the embodiments of the present invention, presented in the figures, is not intended to limit the scope of the invention, as claimed, but is merely representative of selected embodiments of the invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments of the present invention without making any creative effort, shall fall within the protection scope of the present invention.
It should be noted that: like reference numbers and letters refer to like items in the following figures, and thus, once an item is defined in one figure, it need not be further defined and explained in subsequent figures. Meanwhile, in the description of the present invention, the terms "first", "second", and the like are used only for distinguishing the description, and are not to be construed as indicating or implying relative importance.
Referring to fig. 1, a schematic structural diagram of a service terminal 300 according to an embodiment of the present invention is shown, in the embodiment of the present invention, a method for estimating an air temperature is applied to the service terminal 300, and the service terminal 300 may be, but is not limited to, an intelligent electronic device such as a desktop computer. The service terminal 300 includes a device for estimating an air temperature 310, a memory 320, and a processor 330.
The elements of the memory 320 and the processor 330 are electrically connected to each other directly or indirectly to achieve data transmission or interaction. For example, the components may be electrically connected to each other via one or more communication buses or signal lines. The means 310 for estimating the air temperature includes at least one software function module which can be stored in the memory 320 in the form of software or firmware (firmware) or solidified in an Operating System (OS) of the service terminal 300. The processor 330 is used for executing executable modules stored in the memory 320, such as software functional modules and computer programs included in the device 310 for estimating air temperature.
The Memory 320 may be, but is not limited to, a Random Access Memory (RAM), a Read Only Memory (ROM), a Programmable Read-Only Memory (PROM), an Erasable Read-Only Memory (EPROM), an electrically Erasable Read-Only Memory (EEPROM), and the like. The memory 320 is used for storing programs, and the processor 330 executes the programs after receiving the execution instructions.
Referring to fig. 2, a method for estimating air temperature according to an embodiment of the present invention is applied to a service terminal 300, and the method for estimating air temperature includes: and step S110, calculating a first estimated air temperature value of the meteorological station at the preset altitude according to the air temperature value measured by the meteorological station at the actual altitude and the temperature direct reduction rate of the pixel of the meteorological station on the satellite remote sensing image.
By passingCalculating a first estimated air temperature value of the meteorological station at a preset altitude, wherein Ta，nThe measured air temperature value of the meteorological station at the actual altitude,and the temperature direct reduction rate of the pixel of the meteorological site on the satellite remote sensing image is represented.
And step S120, establishing a regression model to obtain a simulated air temperature value of the meteorological site at the preset altitude.
In the embodiment of the present invention, a regression model is established according to the longitude and latitude of the meteorological site and the normalized vegetation index, and the simulated temperature value of the meteorological site at the predetermined altitude is simulated according to the established regression model, but the present invention is not limited thereto, and the regression model may also be established according to other parameters. The regression model is Tm(i，j)＝a1·latitude+a2·lontitude+a3NDVI + b, wherein Tm(i，j)To simulate the gas temperature value, a1,a2,a3For regression coefficients, b is a constant, latitude is longitude, lottitude is latitude, and NDVI is the normalized vegetation index.
Step S130, calculating a second estimated air temperature value of the weather station at the predetermined altitude according to the first estimated air temperature value and the simulated air temperature value of the weather station at the predetermined altitude.
According to the first estimated air temperature value obtained in step S110 and the simulated air temperature value obtained in step S120, a second estimated air temperature value of the weather station at the predetermined altitude is calculated, where the second estimated air temperature value is a more accurate air temperature value of the weather station at the predetermined altitude. Specifically, please refer to fig. 3, which is a flowchart illustrating the sub-step of step S130 according to the embodiment of the present invention.
And S131, calculating an air temperature residual value according to the first estimated air temperature value and the simulated air temperature value of the meteorological station at the preset altitude, and calculating an air temperature residual interpolation value according to the air temperature residual value.
By obtaining TC，n＝TD，n-Tm，nResidual value of air temperature, wherein TC，nIs the residual value of air temperature, TD，nEstimating a first temperature value, T, for said meteorological site at said predetermined altitudem，nAnd simulating the air temperature value of the meteorological station at the preset altitude.
By means of kriging spatial interpolationObtaining an air temperature residual interpolation value, wherein TC(i，j)As interpolated values of air temperature residuals, TC，iIs the residual value of air temperature, deltaiAre observation weight coefficients.
Step S132, adding the air temperature residual error interpolation value and the simulated air temperature to obtain a second estimated air temperature value of the meteorological site at the preset altitude.
Through TD，a(i，j)＝TC(i，j)+Tm(i，j)Calculating by an algorithm to obtain a second estimated air temperature value, wherein TD，a(i，j)For the second estimated gas temperature value, TC(i，j)As interpolated values of air temperature residuals, Tm(i，j)To simulate the gas temperature value.
Referring to fig. 4, a flow chart of another method for estimating air temperature according to an embodiment of the present invention is shown, where the method for estimating air temperature is applied to a service terminal 300, and the method for estimating air temperature includes:
and S210, acquiring a plurality of earth surface temperature values from the satellite remote sensing image by inversion.
And (2) obtaining a plurality of surface temperature values of each pixel from the satellite remote sensing image by MODIS (moderate-resolution imaging spectrometer), wherein the surface temperature inversion values of each pixel comprise 4 inversion values in the day and at night.
Step S220, calculating an average value of a plurality of surface temperature values, wherein the average value of the plurality of surface temperature values is the surface temperature of each pixel.
Through TS＝(TTd+TAd+TTn+TDn) And/4 calculating the surface temperature of each pixel, wherein TTdAnd TAdRepresents the daytime surface temperature, T, of the pixelTnAnd TDnAnd (3) representing the nighttime surface temperature of the pixel, and obtaining the surface temperature of each pixel by inverting the average value of the 4 times surface temperatures of the pixel.
And step S230, calculating the temperature reduction rate of the pixels with different altitudes relative to the pixels with the preset altitude according to the surface temperature of the pixels at the actual altitude and inverted by the remote sensing image and the digital elevation model of the pixels corresponding to the area where the remote sensing image is located.
The method comprises the steps of obtaining the surface temperature of pixels at the actual altitude on a remote sensing image, and meanwhile calculating the temperature direct reduction rate of the pixels at different altitudes relative to the preset altitude by combining a digital elevation model of the pixels, wherein the preset altitude can be selected according to actual needs. In particular by
The algorithm calculates the temperature direct reduction rate of the pixels with different altitudes relative to the preset altitude, whereinDirect rate of temperature decrease, Hi,jAltitude, H, of other picture elementsDIs a predetermined altitude, TS(i，j)Surface temperature, T, of other picture elementsS，D(x，y)Is the surface temperature of the picture element at a predetermined altitude.
Step S240, calculating a first estimated air temperature value of the meteorological station at the preset altitude according to the air temperature value measured by the meteorological station at the actual altitude and the temperature reduction rate of the pixel where the meteorological station is located.
By passingCalculating a first estimated air temperature value of the meteorological station at a preset altitude, wherein Ta，nThe measured air temperature value of the meteorological station at the actual altitude,indicating the rate of temperature decrease at the meteorological site.
Step S250, establishing a regression model to obtain a simulated air temperature value of the meteorological site at the preset altitude.
In the embodiment of the invention, a regression model is established according to the longitude and latitude of the meteorological site and the normalized vegetation index, and the meteorological site is simulated to be in a preset state through the established regression modelThe simulated temperature value of the altitude is not limited to the above, and a regression model can be established through other parameters. The regression model is Tm(i，j)＝a1·latitude+a2·lontitude+a3NDVI + b, wherein Tm(i，j)To simulate the gas temperature value, ai,a2,a3For regression coefficients, b is a constant, latitude is longitude, lottitude is latitude, and NDVI is the normalized vegetation index.
Step S260, calculating a second estimated air temperature value of the weather station at the predetermined altitude according to the first estimated air temperature value and the simulated air temperature value of the weather station at the predetermined altitude.
According to the first estimated air temperature value obtained in step S240 and the simulated air temperature value obtained in step S250, a second estimated air temperature value of the weather station at the predetermined altitude is calculated, where the second estimated air temperature value is a more accurate air temperature value of the weather station at the predetermined altitude. In particular to a method for preparing a high-performance nano-silver alloy,
and subtracting the first estimated air temperature value of the meteorological station at the preset altitude from the simulated air temperature value to obtain an air temperature residual value, and calculating to obtain an air temperature residual interpolation value according to the air temperature residual value.
Through TC，n＝TD，n-Tm，nObtaining a residual air temperature value, wherein TC，nIs the residual value of air temperature, TD，nEstimating a first temperature value, T, for said meteorological site at said predetermined altitudem，nAnd simulating the air temperature value of the meteorological station at the preset altitude.
By means of kriging spatial interpolationObtaining an air temperature residual interpolation value, wherein TC(i，j)As interpolated values of air temperature residuals, TC，iIs the residual value of air temperature, deltaiAre observation weight coefficients.
And adding the air temperature residual error interpolation value and the simulated air temperature to obtain a second estimated air temperature value of the meteorological site at the preset altitude.
Through TD，a(i，j)＝TC(i，j)+Tm(i，j)Calculating by an algorithm to obtain a second estimated air temperature value, wherein TD，a(i，j)For the second estimated gas temperature value, TC(i，j)As interpolated values of air temperature residuals, Tm(i，j)To simulate the gas temperature value.
And step S270, calculating the air temperature values of the actual altitudes of the different pixels according to the temperature direct reduction rates of the pixels with different altitudes and the second estimated air temperature value of the meteorological station at the preset altitude.
The air temperature value at the actual altitude is calculated according to the temperature directly decreasing rate at different altitudes calculated in step S230 and the second estimated air temperature value at the predetermined altitude of the weather station calculated in step S260. In particular by
Calculating the air temperature value of the actual altitude by an algorithm, wherein Ta(i，j)Is the value of the temperature at actual altitude, TD，a(i，j)A second estimated air temperature value for the meteorological site at the predetermined altitude,the temperature decrease rate. Therefore, the air temperature values of other altitudes are obtained through the second estimated air temperature value of the meteorological station at the preset altitude, and the accurate estimation of the air temperature value of the actual altitude is realized.
Referring to fig. 5, which is a functional block diagram of an apparatus 310 for estimating an air temperature according to an embodiment of the present invention, the apparatus 310 for estimating an air temperature is applied to a service terminal 300, and the apparatus 310 for estimating an air temperature includes:
the first calculating module 311 is configured to calculate a first estimated air temperature value of the meteorological site at the predetermined altitude according to the air temperature value measured by the meteorological site at the actual altitude and the temperature directly decreasing rate of the pixel of the meteorological site on the satellite remote sensing image.
In the embodiment of the present invention, step S110, step S210 to step S240 may be executed by the first calculating module 311.
The model establishing module 312 is configured to establish a regression model to obtain a simulated air temperature value of the meteorological site at the predetermined altitude.
In an embodiment of the present invention, steps S120 and S250 may be performed by the model building module 312.
The second calculating module 313 is configured to calculate a second estimated air temperature value of the weather station at the predetermined altitude according to the first estimated air temperature value and the simulated air temperature value of the weather station at the predetermined altitude.
In the embodiment of the present invention, steps S130 to S132, step S260, and step S270 may be executed by the second calculation module 313.
Since the method of estimating the air temperature has been described in detail, it will not be described in detail.
In summary, the method for estimating air temperature according to the embodiments of the present invention is applied to a service terminal, and the method for estimating air temperature includes calculating a first estimated air temperature value of a weather station at a predetermined altitude according to an air temperature value measured at an actual altitude of the weather station and a temperature droop rate of the weather station. Establishing a regression model to obtain a simulated air temperature value of the meteorological site at the preset altitude, and calculating a second estimated air temperature value of the meteorological site at the preset altitude according to a first estimated air temperature value and the simulated air temperature value of the meteorological site at the preset altitude. According to the scheme, the temperature is accurately estimated, and errors caused by directly acquiring temperature data are avoided.
In the embodiments provided in the present application, it should be understood that the disclosed apparatus and method can be implemented in other ways. The apparatus embodiments described above are merely illustrative, and for example, the flowchart and block diagrams in the figures illustrate the architecture, functionality, and operation of possible implementations of apparatus, methods and computer program products according to various embodiments of the present invention. In this regard, each block in the flowchart or block diagrams may represent a module, segment, or portion of code, which comprises one or more executable instructions for implementing the specified logical function(s). It should also be noted that, in some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams and/or flowchart illustration, and combinations of blocks in the block diagrams and/or flowchart illustration, can be implemented by special purpose hardware-based systems which perform the specified functions or acts, or combinations of special purpose hardware and computer instructions.
In addition, the functional modules in the embodiments of the present invention may be integrated together to form an independent part, or each module may exist separately, or two or more modules may be integrated to form an independent part.
The functions, if implemented in the form of software functional modules and sold or used as a stand-alone product, may be stored in a computer readable storage medium. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes. It is noted that, herein, relational terms such as first and second, and the like may be used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Also, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other identical elements in a process, method, article, or apparatus that comprises the element.
The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention. It should be noted that: like reference numbers and letters refer to like items in the following figures, and thus, once an item is defined in one figure, it need not be further defined and explained in subsequent figures.
The above description is only for the specific embodiments of the present invention, but the scope of the present invention is not limited thereto, and any person skilled in the art can easily conceive of the changes or substitutions within the technical scope of the present invention, and all the changes or substitutions should be covered within the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.
Claims (8)
1. A method for estimating air temperature, which is applied to a service terminal and is used for combining and analyzing earth surface temperature data inverted on a satellite remote sensing image and air temperature measured by a meteorological site to accurately estimate the air temperature, and the method comprises the following steps:
calculating a first estimated air temperature value of the meteorological station at a preset altitude according to the air temperature value measured by the meteorological station at the actual altitude and the temperature reduction rate of the pixel of the meteorological station on the satellite remote sensing image;
establishing a regression model to obtain a simulated air temperature value of the meteorological station at the preset altitude;
calculating an air temperature residual value by using a first estimated air temperature value of the meteorological station at the preset altitude and the simulated air temperature value, and calculating an air temperature residual interpolation value according to the air temperature residual value;
adding the air temperature residual error interpolation value and the simulated air temperature value to obtain a second estimated air temperature value of the meteorological station at the preset altitude;
and calculating the gas temperature value of the actual altitude of the pixel according to the temperature direct reduction rate of the pixel on the satellite remote sensing image and the second estimated gas temperature value of the meteorological site at the preset altitude.
2. The method of estimating air temperature according to claim 1, characterized in that,
the calculation mode of calculating the air temperature residual error is as follows:
TC，n＝TD，n-Tm，n
wherein, TC，nIs the air temperature residual error, TD，nEstimating a first temperature value, T, for said meteorological site at said predetermined altitudem，nSimulating a temperature value for the meteorological site at the predetermined altitude;
the calculation mode of obtaining the air temperature residual space interpolation by calculating according to the air temperature residual is as follows:
wherein, TC(i，j)As interpolated values of air temperature residuals, TC，iIs the air temperature residual, δiAre observation weight coefficients.
3. The method of estimating air temperature according to claim 1, characterized in that,
the calculation mode of the actual altitude air temperature value of the pixel is as follows:
wherein, Ta(i，j)Is the value of the temperature at actual altitude, TD，a(i，j)A second estimated air temperature value for the meteorological site at the predetermined altitude,the temperature decrease rate.
4. The method of estimating air temperature according to claim 1, characterized in that it further comprises:
and calculating the temperature reduction rate of the pixels with different altitudes relative to the pixels at the preset altitude by combining the digital elevation model of the pixels in the corresponding area of the remote sensing image according to the surface temperature of the pixels at the actual altitude obtained by inversion on the satellite remote sensing image.
5. The method of estimating air temperature according to claim 4, characterized in that,
the calculation mode of the temperature direct reduction rate is as follows:
wherein,is the rate of direct decrease in temperature, Hi,jAltitude, H, of other picture elementsDIs a predetermined altitude, TS(i，j)Surface temperature, T, of other picture elementsS，D(x，y)Is the surface temperature of the picture element at a predetermined altitude.
6. The method of estimating air temperature according to claim 4, characterized in that it further comprises:
a plurality of earth surface temperature values of each pixel are obtained by inversion from the satellite remote sensing image;
and calculating the average value of a plurality of surface temperature values, wherein the average value of the plurality of surface temperature values is the surface temperature of each pixel.
7. The method of estimating air temperature according to claim 1, wherein said step of building a regression model to obtain a simulated air temperature value at said predetermined altitude for said meteorological site comprises:
establishing a regression model according to the longitude and the latitude of the meteorological site and the normalized vegetation index to obtain a simulated air temperature value of the meteorological site at the preset altitude;
the calculation mode for establishing the regression model is as follows:
Tm(i，j)＝a1·latitude+a2·lontitude+a3·NDVI+b
wherein, Tm(i，j)To simulate the gas temperature value, a1,a2,a3For regression coefficients, b is a constant, latitude is longitude, lottitude is latitude, and NDVI is the normalized vegetation index.
8. A service terminal, characterized in that the service terminal comprises:
a memory;
a processor; and a means for estimating air temperature, said means for estimating air temperature being installed in said memory and comprising one or more software functional modules executed by said processor, said means for estimating air temperature comprising:
the first calculation module is used for calculating a first estimated air temperature value of the meteorological station at a preset altitude according to the air temperature value measured by the meteorological station at the actual altitude and the temperature direct reduction rate of the pixel of the meteorological station on the satellite remote sensing image;
the model establishing module is used for establishing a regression model to obtain a simulated air temperature value of the meteorological station at the preset altitude;
the second calculation module is used for calculating an air temperature residual value according to the first estimated air temperature value and the simulated air temperature value of the meteorological station at the preset altitude and calculating an air temperature residual interpolation value according to the air temperature residual value; the second calculation module is further used for adding the air temperature residual error interpolation value and the simulated air temperature value to obtain a second estimated air temperature value of the meteorological site at the preset altitude; the second calculation module is further used for calculating the air temperature value of the actual altitude of the pixel according to the temperature reduction rate of the pixel on the satellite remote sensing image and the second estimated air temperature value of the meteorological site at the preset altitude.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710434299.5A CN107247690B (en) | 2017-06-09 | 2017-06-09 | Estimate the method and service terminal of temperature |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710434299.5A CN107247690B (en) | 2017-06-09 | 2017-06-09 | Estimate the method and service terminal of temperature |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107247690A CN107247690A (en) | 2017-10-13 |
| CN107247690B true CN107247690B (en) | 2019-10-01 |
Family
ID=60018215
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710434299.5A Active CN107247690B (en) | 2017-06-09 | 2017-06-09 | Estimate the method and service terminal of temperature |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107247690B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108961089A (en) * | 2018-05-28 | 2018-12-07 | 浙江大学 | Based on star the In Middle And Lower Reaches of Changjiang River rice Climatic regionalization remote-sensing monitoring method of multi-source data |
| CN109063684B (en) * | 2018-08-28 | 2022-03-29 | 中国科学院西北生态环境资源研究院 | Near-surface air temperature remote sensing monitoring method |
| CN109871637B (en) * | 2019-03-06 | 2022-11-04 | 成都信息工程大学 | Near-ground air temperature estimation method under cloud-sky condition |
| CN114117341A (en) * | 2020-09-01 | 2022-03-01 | 石河子大学 | Calculation method for annual average temperature in any region of Tianshan northern slope |
| CN114021371B (en) * | 2021-11-16 | 2023-03-03 | 中国科学院西北生态环境资源研究院 | Carbon reserve influence estimation method and device, electronic equipment and storage medium |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103530499A (en) * | 2013-08-29 | 2014-01-22 | 西南林业大学 | Method for building mountainous area surface temperature base line and application |
| CN106169014A (en) * | 2016-06-15 | 2016-11-30 | 中国水利水电科学研究院 | Region based on remotely-sensed data Surface sensible heat/latent heat flux inversion method and system |
-
2017
- 2017-06-09 CN CN201710434299.5A patent/CN107247690B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103530499A (en) * | 2013-08-29 | 2014-01-22 | 西南林业大学 | Method for building mountainous area surface temperature base line and application |
| CN106169014A (en) * | 2016-06-15 | 2016-11-30 | 中国水利水电科学研究院 | Region based on remotely-sensed data Surface sensible heat/latent heat flux inversion method and system |
Non-Patent Citations (3)
| Title |
|---|
| Determination of land surface temperature and its lapse rate in the Satluj River basin using NOAA data;Sanjay K. Jain 等;《International Journal of Remote Sensing》;20080630;第29卷(第11期);第3091-3103页 * |
| 基于GIS的日照市气温分布式模拟;王林林 等;《地理与地理信息科学》;20080915;第24卷(第5期);第47-50页 * |
| 基于MODIS遥感影像的云南山地气温垂直直减率场;王艳霞 等;《遥感技术与应用》;20140815;第29卷(第4期);第616-624页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107247690A (en) | 2017-10-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107247690B (en) | Estimate the method and service terminal of temperature | |
| Demuth et al. | Evaluation of Advanced Microwave Sounding Unit tropical-cyclone intensity and size estimation algorithms | |
| Xu et al. | A new satellite-based monthly precipitation downscaling algorithm with non-stationary relationship between precipitation and land surface characteristics | |
| Machado et al. | The convective system area expansion over Amazonia and its relationships with convective system life duration and high-level wind divergence | |
| CN108535731B (en) | Method and device for forecasting short rainfall | |
| JPWO2016031174A1 (en) | Simulation device, simulation method, and storage medium | |
| CN110390135B (en) | Method for improving forest fire spreading prediction precision | |
| CN110321784B (en) | Method, apparatus, electronic device and computer medium for soil moisture estimation | |
| Chabot-Couture et al. | An environmental data set for vector-borne disease modeling and epidemiology | |
| CN114868043B (en) | Method, apparatus, system and computer program product for estimating road surface temperature | |
| Pimentel et al. | Study of snow dynamics at subgrid scale in semiarid environments combining terrestrial photography and data assimilation techniques | |
| Jose et al. | Bias Correction and trend analysis of temperature data by a high-resolution CMIP6 Model over a Tropical River Basin | |
| Engelhardt et al. | Evaluation of gridded precipitation for norway using glacier mass‐balance measurements | |
| Balsa-Barreiro et al. | A new methodology to estimate the discrete-return point density on airborne LiDAR surveys | |
| West et al. | Exploring the sensitivity of coastal inundation modelling to DEM vertical error | |
| Holmlund et al. | Constraining 135 years of mass balance with historic structure-from-motion photogrammetry on Storglaciären, Sweden | |
| KR102580312B1 (en) | Observation data quality inspection apparatus and observation data quality inspection method using the same | |
| US10339691B2 (en) | Storage of time-varying data as parametrized functions | |
| Osejo et al. | Spatial distribution of precipitation and evapotranspiration estimates from Worldclim and Chelsa datasets: Improving long-term water balance at the watershed-scale in the Urabá region of Colombia | |
| CN110334381B (en) | Method, apparatus and computer storage medium for evapotranspiration estimation | |
| Xu et al. | Merging the MODIS and landsat terrestrial latent heat flux products using the multiresolution tree method | |
| US10317575B2 (en) | Method and apparatus for forecasting weather | |
| Lee et al. | Estimation of maximum daily fresh snow accumulation using an artificial neural network model | |
| Duffaut Espinosa et al. | Multifractal downscaling of rainfall using normalized difference vegetation index (NDVI) in the Andes Plateau | |
| Shi et al. | Integrating ASTER and GLASS broadband emissivity products using a multi-resolution Kalman filter |
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 | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20200703 Address after: 730000 No.318, Donggang West Road, Chengguan District, Lanzhou City, Gansu Province Patentee after: NORTHWEST INSTITUTE OF ECO-ENVIRONMENT AND RESOURCES, CAS Address before: Chengguan District of Gansu city of Lanzhou province Donggang West Road 730000 No. 320 Patentee before: Institute of environment and Engineering in cold and dry areas, Chinese Academy of Sciences |