Method, device, electronic device and storage medium for determining air temperature

Abstract

The present application provides an air temperature determination method, device, electronic device and storage medium, the method comprising: obtaining first air temperature data of multiple first sampling points and second air temperature data of multiple second sampling points; establishing a first air temperature semivariogram based on the first air temperature data; establishing a first surface temperature semivariogram based on the first surface temperature data; establishing multiple third sampling points and establishing a second surface temperature semivariogram based on the second surface temperature data; determining an unbiased second air temperature semivariogram based on the first air temperature semivariogram, the first surface temperature semivariogram and the second surface temperature semivariogram; determining an air temperature range of the third sampling point; determining the air temperature of the sampling point to be measured based on the sampling point to be measured, the second sampling point, the second air temperature data of the second sampling point, the third sampling point, the air temperature range of the third sampling point and the unbiased second air temperature semivariogram.

Description

Technical Field

[0001] The present application relates to the field of data processing, and in particular to a method, device, electronic device and storage medium for determining air temperature. Background Art

[0002] Spatial interpolation allows humans to estimate the population through samples and thus mathematically reflect the real world. Spatial interpolation methods are widely used in geography, environmental science, and engineering. However, traditional methods are difficult to meet complex real-world situations. For example, spatial stratification heterogeneity is widely present in nature, and the current mainstream method is Kriging, which relies heavily on the second-order stationary assumption of spatial objects. However, since this assumption ignores the existence of stratified heterogeneity and does not consider the quality of samples, it is difficult to achieve satisfactory results in high-precision spatial interpolation. In addition, other spatial interpolation methods such as inverse distance weighting (IDW) and nearest neighbor method (NN) are also directly affected by sample quality.

[0003] When the sample population is stratified and heterogeneous, and samples in some areas are extremely sparse or even missing, the observed samples are biased and cannot make accurate global or local statistical inferences. In this case, all of the above methods cannot obtain accurate spatial modeling, nor can they determine accurate air temperature data through interpolation. Summary of the invention

[0004] The present application provides a method, device, electronic device and storage medium for determining air temperature.

[0005] An embodiment of the present application provides, on one hand, a method for determining air temperature, the method comprising:

[0006] Obtaining first air temperature data of a plurality of first sampling points in a first region and second air temperature data of a plurality of second sampling points in a second region, wherein the first region and the second region belong to the same target region but belong to different climate zones, and the number of the second sampling points is less than a target ratio of the sum of the first sampling points and the second sampling points;

[0007] establishing a first air temperature semivariogram based on the first air temperature data;

[0008] Collecting the surface temperature data of the first sampling point by using a satellite to obtain the first surface temperature data of the first sampling point;

[0009] Establishing a first land surface temperature semivariogram based on the first land surface temperature data;

[0010] Establishing a plurality of third sampling points in the second area, and collecting surface temperature data of the third sampling points by using a satellite to obtain second surface temperature data of the third sampling points;

[0011] Establishing a second land surface temperature semivariogram based on the second land surface temperature data;

[0012] determining an unbiased second air temperature semivariogram based on the first air temperature semivariogram, the first land surface temperature semivariogram, and the second land surface temperature semivariogram;

[0013] Determining the air temperature range of the third sampling point using linear regression;

[0014] The air temperature of the sampling point to be measured is determined based on the sampling point to be measured, the second sampling point, the second air temperature data of the second sampling point, the third sampling point, the air temperature range of the third sampling point and the unbiased second air temperature semivariogram.

[0015] Wherein, establishing a first air temperature semivariogram based on the first air temperature data comprises:

[0016] Traverse multiple preset range values;

[0017] Determine at least one first sampling point pair whose distance is less than or equal to the current preset range value;

[0018] Determine a function value corresponding to the current preset range value according to at least one first sampling point pair of the current preset range value and first air temperature data of each first sampling point;

[0019] After the traversal is completed, the function value corresponding to each preset range value is obtained;

[0020] A first air temperature semivariogram is determined according to the plurality of preset range values ​​and a function value corresponding to each preset range value.

[0021] Wherein, determining the function value corresponding to the current preset range value according to at least one first sampling point pair of the current preset range value and the first air temperature data of each first sampling point includes:

[0022] Determine the square of the air temperature difference of each first sampling point pair according to the first air temperature data of each first sampling point;

[0023] The squares of the air temperature differences of all first sampling point pairs are summed, and the sum result is divided by twice the total number of first sampling point pairs to obtain the function value corresponding to the current preset range value.

[0024] The step of establishing a plurality of third sampling points in the second area includes:

[0025] constructing a window covering the second area;

[0026] A grid of windows is constructed based on a preset interval;

[0027] determining a plurality of candidate third sampling points based on the grid;

[0028] A candidate third sampling point in the second area is determined as the third sampling point.

[0029] Wherein, determining an unbiased second air temperature semivariogram based on the first air temperature semivariogram, the first land surface temperature semivariogram, and the second land surface temperature semivariogram comprises:

[0030] determining a semivariogram scale based on the first air temperature semivariogram and the first land surface temperature semivariogram;

[0031] An unbiased second air temperature semivariogram is determined based on the semivariogram scale and the second land surface temperature semivariogram.

[0032] Wherein, determining the air temperature range of the third sampling point by using linear regression includes:

[0033] collecting the land surface temperature data of the second sampling point by using a satellite to obtain third land surface temperature data of the second sampling point;

[0034] Establishing a linear regression function based on the first sampling point, the first air temperature data and the first ground surface temperature data of the first sampling point, the second sampling point and the second air temperature data and the third ground surface temperature data of the second sampling point, wherein the linear regression function represents the mapping relationship between the air temperature and the ground surface temperature in the target area;

[0035] An air temperature range of the third sampling point is determined based on the third sampling point, the second ground surface temperature data of the third sampling point, and the linear regression function.

[0036] Another aspect of the present application provides an air temperature determination device, the device comprising:

[0037] A collection module, configured to obtain first air temperature data of a plurality of first sampling points in a first area and second air temperature data of a plurality of second sampling points in a second area, wherein the first area and the second area belong to the same target area but belong to different climate zones, and the number of the second sampling points is less than a target ratio of the sum of the first sampling points and the second sampling points; collect surface temperature data of the first sampling points by using a satellite to obtain first surface temperature data of the first sampling points; and establish a plurality of third sampling points in the second area, and collect surface temperature data of the third sampling points by using a satellite to obtain second surface temperature data of the third sampling points;

[0038] A calculation module is used to establish a first air temperature semivariogram based on the first air temperature data; establish a first surface temperature semivariogram based on the first surface temperature data; establish a second surface temperature semivariogram based on the second surface temperature data; determine an unbiased second air temperature semivariogram based on the first air temperature semivariogram, the first surface temperature semivariogram and the second surface temperature semivariogram; determine the air temperature range of the third sampling point by linear regression; and determine the air temperature of the sampling point to be measured based on the sampling point to be measured, the second sampling point, the second air temperature data of the second sampling point, the third sampling point, the air temperature range of the third sampling point and the unbiased second air temperature semivariogram.

[0039] The acquisition module is further used to acquire the surface temperature data of the second sampling point by using a satellite to obtain the third surface temperature data of the second sampling point;

[0040] The calculation module is further used to establish a linear regression function based on the first sampling point, the first air temperature data and the first ground surface temperature data of the first sampling point, the second sampling point and the second air temperature data and the third ground surface temperature data of the second sampling point, wherein the linear regression function represents the mapping relationship between the air temperature and the ground surface temperature in the target area; and determine the air temperature range of the third sampling point based on the third sampling point, the second ground surface temperature data of the third sampling point and the linear regression function.

[0041] Another aspect of the present invention provides an electronic device, comprising: a processor and an input-output system;

[0042] The input-output system obtains first air temperature data of a plurality of first sampling points in the first area and second air temperature data of a plurality of second sampling points in the second area; collects surface temperature data of the first sampling points by using a satellite to obtain first surface temperature data of the first sampling points; and establishes a plurality of third sampling points in the second area, and collects surface temperature data of the third sampling points by using a satellite to obtain second surface temperature data of the third sampling points;

[0043] The processor establishes a first air temperature semivariogram based on the first air temperature data; establishes a first surface temperature semivariogram based on the first surface temperature data; establishes a second surface temperature semivariogram based on the second surface temperature data; determines an unbiased second air temperature semivariogram based on the first air temperature semivariogram, the first surface temperature semivariogram and the second surface temperature semivariogram; determines the air temperature range of the third sampling point using linear regression; and determines the air temperature of the sampling point to be measured based on the sampling point to be measured, the second sampling point, the second air temperature data of the second sampling point, the third sampling point, the air temperature range of the third sampling point and the unbiased second air temperature semivariogram.

[0044] On the other hand, the present invention provides a computer-readable storage medium, the storage medium storing a computer program, the computer program being used to obtain first air temperature data of multiple first sampling points in a first area and second air temperature data of multiple second sampling points in a second area; establish a first air temperature semivariogram based on the first air temperature data; use a satellite to collect the surface temperature data of the first sampling point to obtain the first surface temperature data of the first sampling point; establish a first surface temperature semivariogram based on the first surface temperature data; establish multiple third sampling points in the second area, and use a satellite to collect the surface temperature data of the third sampling point to obtain the second surface temperature data of the third sampling point; establish a second surface temperature semivariogram based on the second surface temperature data; determine an unbiased second air temperature semivariogram based on the first air temperature semivariogram, the first surface temperature semivariogram and the second surface temperature semivariogram; determine the air temperature range of the third sampling point using linear regression; and determine the air temperature of the sampling point to be measured based on the sampling point to be measured, the second sampling point, the second air temperature data of the second sampling point, the third sampling point, the air temperature range of the third sampling point and the unbiased second air temperature semivariogram.

[0045] It should be understood that the content described in this section is not intended to identify the key or important features of the embodiments of the present application, nor is it intended to limit the scope of the present application. Other features of the present application will become easily understood through the following description. BRIEF DESCRIPTION OF THE DRAWINGS

[0046] By reading the detailed description below with reference to the accompanying drawings, the above and other purposes, features and advantages of the exemplary embodiments of the present application will become readily understood. In the accompanying drawings, several embodiments of the present application are shown in an exemplary and non-limiting manner, wherein:

[0047] In the drawings, the same or corresponding reference numerals represent the same or corresponding parts.

[0048] Figure 1 A flow chart showing a method for determining air temperature according to an embodiment of the present application is shown;

[0049] Figure 2 A flow chart showing a method for determining air temperature according to another embodiment of the present application is shown;

[0050] Figure 3 A schematic diagram showing a first air temperature semivariogram according to an embodiment of the present application is shown;

[0051] Figure 4 A flow chart showing a method for determining air temperature according to another embodiment of the present application is shown;

[0052] Figure 5 A flow chart showing a method for determining air temperature according to another embodiment of the present application is shown;

[0053] Fig. 6A A schematic diagram of constructing a second area window according to an embodiment of the present application is shown;

[0054] Figure 6B A schematic diagram of constructing a second area window grid according to an embodiment of the present application is shown;

[0055] Figure 7 A flow chart showing a method for determining air temperature according to another embodiment of the present application is shown;

[0056] Figure 8 A flow chart showing a method for determining air temperature according to another embodiment of the present application is shown;

[0057] Fig. 9 A schematic structural diagram of an air temperature determination device according to an embodiment of the present application is shown;

[0058] Fig.10 A schematic diagram of the structure of an electronic device according to an embodiment of the present application is shown. DETAILED DESCRIPTION

[0059] In order to make the purpose, features, and advantages of the present application more obvious and easy to understand, the technical solutions in the embodiments of the present application will be clearly and completely described below in conjunction with the drawings in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative work are within the scope of protection of the present application.

[0060] In order to determine accurate air temperature data through spatial modeling and interpolation when there is stratified heterogeneity in the sample population and samples in some areas are extremely sparse or even missing, an embodiment of the present application provides an air temperature determination method, such as Figure 1 As shown, the method includes:

[0061] Step 101, obtaining first air temperature data of a plurality of first sampling points in a first area and second air temperature data of a plurality of second sampling points in a second area, wherein the first area and the second area belong to the same target area but to different climate zones, and the number of the second sampling points is less than the target ratio of the sum of the first sampling points and the second sampling points.

[0062] For example, the western part of the Qinghai-Tibet Plateau is dominated by a plateau subarctic climate, while the eastern and southern parts are dominated by a plateau temperate climate. There is significant stratified heterogeneity in these two climate regions. Most of the meteorological stations are concentrated in the east, which results in the collected sample data mainly reflecting the characteristics of the plateau temperate climate, and cannot reflect the characteristics of the plateau subarctic climate. A total of 51 sample data of the Qinghai-Tibet Plateau were obtained. Based on the classification of different climate zones, the east and south were divided into the first region, and the sample data of the first sampling point belonging to the first region were determined as the first air temperature data, totaling 48. The west was divided into the second region, and the sample data of the second sampling point belonging to the second region were determined as the second air temperature data, totaling 3.

[0063] In this embodiment, the target ratio is set to 10%. In other implementations, the target ratio can be set based on specific needs.

[0064] Step 102: establishing a first air temperature semivariogram based on the first air temperature data.

[0065] Semivariogram (also called semivariogram) is a calculation method used in mathematical modeling, and its expression is:

[0066]

[0067] Wherein, when establishing the first air temperature semivariogram, γ(h) is the function value of the semivariogram, h is the preset range value, N(h) is the number of first sampling point pairs whose distance is less than or equal to the range value h (i.e., point pairs composed of two different first sampling points, and the distance is the distance between the two first sampling points), Z(X α ) indicates that at the first sampling point X α The first air temperature data at Z(X α +h) indicates the current and first sampling point X α The first air temperature data at the first sampling point for calculation.

[0068] Each preset range value is traversed and substituted into the semivariogram, and the function value corresponding to each preset range value is obtained based on the distance between each two first sampling points and the first air temperature data of each first sampling point.

[0069] Fitting is performed based on multiple preset range values ​​and function values ​​corresponding to each preset range value to obtain a first air temperature semivariogram.

[0070] For example, for the 48 first sampling points in the first area of ​​the Qinghai-Tibet Plateau and the first air temperature data of each first sampling point, after calculation, the first air temperature semivariogram is:

[0071]

[0072] Among them, γ east (h) is the function value of the first air temperature semivariogram, and h is the distance between the first sampling points.

[0073] Step 103: Collect the land surface temperature data of the first sampling point by using a satellite to obtain the first land surface temperature data of the first sampling point.

[0074] The first land surface temperature data (LST data) of each first sampling point is obtained by using a remote sensing image of the land surface temperature at each first sampling point observed by a satellite.

[0075] Step 104: establishing a first land surface temperature semivariogram based on the first land surface temperature data.

[0076] A first land surface temperature semivariogram is established based on the observed first land surface temperature data at each first sampling point location.

[0077] Similarly, the expression based on the semivariogram is:

[0078]

[0079] When establishing the first surface temperature semivariogram, γ(h) is the function value of the semivariogram, h is the preset range value, N(h) is the number of first sampling point pairs whose distance is less than or equal to the range value h, and Z(X α ) indicates that at the first sampling point X α The first surface temperature data at Z(X α +h) indicates the current and first sampling point X α The first surface temperature data at the first sampling point for calculation.

[0080] Each preset range value is traversed and substituted into the semivariogram, and the function value corresponding to each preset range value is obtained based on the distance between each two first sampling points and the first surface temperature data of each first sampling point.

[0081] Fitting is performed based on multiple preset range values ​​and the function value corresponding to each preset range value to obtain a first surface temperature semivariogram.

[0082] For example, for the 48 first sampling points in the first area of ​​the Qinghai-Tibet Plateau and the first surface temperature data of each first sampling point, after calculation, the first surface temperature semivariogram is:

[0083]

[0084] Among them, γ east,LST (h) is the function value of the first surface temperature semivariogram, and h is the distance between the first sampling points.

[0085] Step 105: establish a plurality of third sampling points in the second area, and collect land surface temperature data of the third sampling points by using a satellite to obtain second land surface temperature data of the third sampling points.

[0086] Since the number of second sampling points in the second area is small, based on the rare second sampling points and the second air temperature data of the second sampling points, the surface temperature data obtained by satellite measurement of the second sampling points will also be small, and it is impossible to directly establish an accurate air temperature semivariogram or surface temperature semivariogram.

[0087] Therefore, multiple third sampling points are established in the second area, and multiple third sampling points in the second area can be established by randomly selecting in the second area or by constructing a grid. Then, the remote sensing image of the surface temperature at each third sampling point is observed by satellite to obtain the second surface temperature data of each third sampling point.

[0088] Step 106: Establish a second land surface temperature semivariogram based on the second land surface temperature data.

[0089] A second land surface temperature semivariogram is established based on the observed second land surface temperature data at each third sampling point location.

[0090] Similarly, the expression based on the semivariogram is:

[0091]

[0092] When establishing the second surface temperature semivariogram, γ(h) is the function value of the semivariogram, h is the preset range value, N(h) is the number of third sampling point pairs whose distance is less than or equal to the range value h, and Z(X α ) indicates that at the third sampling point X α The second surface temperature data at Z(X α +h) indicates the current and the third sampling point X αThe second surface temperature data at the third sampling point is calculated.

[0093] Each preset range value is traversed and substituted into the semivariogram, and the function value corresponding to each preset range value is obtained based on the distance between every two third sampling points and the second surface temperature data of each third sampling point.

[0094] Fitting is performed based on multiple preset range values ​​and the function value corresponding to each preset range value to obtain a second surface temperature semivariogram.

[0095] For example, for the three second sampling points in the second area of ​​the Qinghai-Tibet Plateau, the number is small, and 50 third sampling points in the second area are established by constructing a grid. Based on the 50 third sampling points in the second area and the second surface temperature data of each third sampling point, after calculation, the second surface temperature semivariogram is:

[0096]

[0097] Among them, γ west，LST (h) is the function value of the second surface temperature semivariogram, and h is the distance between the third sampling points.

[0098] Step 107: Determine an unbiased second air temperature semivariogram based on the first air temperature semivariogram, the first land surface temperature semivariogram, and the second land surface temperature semivariogram.

[0099] After obtaining the first air temperature semivariogram, the first surface temperature semivariogram, and the second surface temperature semivariogram, the unbiased second air temperature semivariogram can be determined by the following formula:

[0100]

[0101] For example, the first air temperature semivariogram is:

[0102]

[0103] The first surface temperature semivariogram is:

[0104]

[0105] The second surface temperature semivariogram is:

[0106]

[0107] Substituting the first air temperature semivariogram, the first surface temperature semivariogram, and the second surface temperature semivariogram into the above formula for determining the unbiased second air temperature semivariogram, the unbiased second air temperature semivariogram is obtained as:

[0108]

[0109] Step 108: Determine the air temperature range of the third sampling point by using linear regression.

[0110] A linear regression function of the target area is constructed based on the first air temperature data and the first surface temperature data of the first sampling point and the second air temperature data and the second surface temperature data of the second sampling point. The linear regression function can characterize the mapping relationship between the air temperature and the surface temperature of the target area.

[0111] Then, the air temperature range of the third sampling point is determined by using the linear regression function and the second ground surface temperature data of the third sampling point.

[0112] Step 109, determining the air temperature of the sampling point to be measured based on the sampling point to be measured, the second sampling point, the second air temperature data of the second sampling point, the third sampling point, the air temperature range of the third sampling point and the unbiased second air temperature semivariogram.

[0113] The sampling point to be measured, the second sampling point, the second air temperature data of the second sampling point, the third sampling point, the air temperature range of the third sampling point and the unbiased second air temperature semivariogram are used as input data, and the Bayesian maximum entropy calculation is performed using BMELIB2.0c (Bayesian Maximum Entropy Software Package) in matlab (a mathematical software used in data analysis, wireless communication, deep learning, image processing and computer vision, signal processing, quantitative finance and risk management, robots, control systems and other fields.) And through the calculation instruction: BMEintervalMode (targetcoors, tempcoors, softcoors, tempval, slow, sup, {'sphericalC', 'gaussianC'}, {paramwest1, paramwest2}, 5, 5, 400, -1). Through this calculation instruction, the air temperature of the sampling point to be measured is determined.

[0114] In the above scheme, since the air temperature can only be obtained through monitoring of the set meteorological stations, when there are sparse meteorological stations in different climate zones of the target area, it is impossible to establish an accurate model through the sample data of the location, and then determine the accurate air temperature of the corresponding location. The first air temperature semivariogram is established based on the first air temperature data of the first sampling point in the first area with more meteorological stations, and the first surface temperature data of the first sampling point is collected by satellite, and the first surface temperature semivariogram is established. Then, multiple third sampling points are established in the second area with sparse meteorological stations, and the second surface temperature data of the third sampling point is collected by satellite, and the second surface temperature semivariogram is established. Since there is a corresponding correlation between the air temperature and the surface temperature in the same area, the unbiased second air temperature semivariogram of the second area can be determined by the first air temperature semivariogram, the first surface temperature semivariogram and the second surface temperature semivariogram. Then, the air temperature range of the third sampling point is determined by linear regression. Finally, the accurate air temperature of the sampling point to be measured is determined based on the sampling point to be measured, the second sampling point, the second air temperature data of the second sampling point, the third sampling point, the air temperature range of the third sampling point and the unbiased second air temperature semivariogram.

[0115] In an example of the present application, a method for determining air temperature is also provided, such as Figure 2 As shown, the establishing of the first air temperature semivariogram based on the first air temperature data includes:

[0116] Step 201, traverse multiple preset range values.

[0117] In this embodiment, an initial range value and a final range value may be set to determine the range of the range value, and then a step size may be set to determine a plurality of range values.

[0118] For example, if the initial range value is set to 10, the final range value is set to 100, and the step length is 10, then 10, 20, 30, 40, 50, 60, 70, 80, 90 and 100 are obtained, a total of 10 preset range values.

[0119] In this embodiment, for the convenience of description, the preset range value currently traversed is referred to as the current preset range value.

[0120] Step 202: Determine at least one first sampling point pair whose distance is less than or equal to the current preset range value.

[0121] For example, there are first sampling points A, B and C in the first area, and the three first sampling points form three first sampling point pairs, namely, first sampling point pair AB, first sampling point pair AC and first sampling point pair BC. The distance corresponding to the first sampling point pair AB is 5 km, the distance between the first sampling point pair AC is 12 km, and the distance between the first sampling point pair BC is 8 km. If the current preset range value is 10 km, the first sampling point pair AB and the first sampling point pair BC are determined as the first sampling point pairs of the current preset range value.

[0122] Step 203: determining a function value corresponding to the current preset range value according to at least one first sampling point pair of the current preset range value and first air temperature data of each first sampling point.

[0123] The square of the difference between the first air temperature data of the two first sampling points in each first sampling point pair of the current preset range value is calculated and summed, and the sum result is divided by twice the number of first sampling point pairs corresponding to the current preset range value to obtain the function value corresponding to the current preset range value.

[0124] Step 204: After the traversal is completed, the function value corresponding to each preset range value is obtained.

[0125] Step 205: determining a first air temperature semivariogram according to the plurality of preset range values ​​and a function value corresponding to each preset range value.

[0126] like Figure 3 As shown, a coordinate system is established with the range value as the X-axis and the function value as the Y-axis, the preset range value is used as the horizontal axis, and the corresponding function value is used as the vertical axis to generate multiple coordinate points ( Figure 3 The black dots in the figure are placed in the coordinate system, and curve fitting is performed based on these multiple coordinate points ( Figure 3 The curve in is the fitting result), and the first air temperature semivariogram of the first area is obtained.

[0127] In an example of the present application, a method for determining air temperature is also provided, such as Figure 4 As shown, the determining of the function value corresponding to the current preset range value according to at least one first sampling point pair of the current preset range value and the first air temperature data of each first sampling point includes:

[0128] Step 301: Determine the square of the air temperature difference between each first sampling point pair according to the first air temperature data of each first sampling point.

[0129] For example, in the first sampling point pair AB, the first air temperature data of the first sampling point A is 5 degrees, and the first air temperature data of the first sampling point B is 0 degrees, then the square of the difference in the first air temperature data of the first sampling point pair AB is 25.

[0130] Step 302: sum the squares of the air temperature differences of all first sampling point pairs, and divide the sum by twice the total number of first sampling point pairs to obtain a function value corresponding to the current preset range value.

[0131] For example, there are a first sampling point pair AB, a first sampling point pair AC, and a first sampling point pair BC, and the squares of the first air temperature data differences are 27, 63, and 30, respectively. The squares of these first air temperature data differences are summed to obtain a sum result of 120, and the sum result is divided by twice the total number of first sampling point pairs, that is, 6, to obtain a function value of 20 corresponding to the current preset range value.

[0132] In an example of the present application, a method for determining air temperature is also provided, such as Figure 5 As shown, the establishing of a plurality of third sampling points in the second area comprises:

[0133] Step 401: construct a window covering the second area.

[0134] like Fig. 6A As shown, for the second area, a window capable of covering the second area is constructed.

[0135] Step 402: constructing a grid of the windows based on a preset interval.

[0136] In this embodiment, the preset interval is generally determined based on the size of the second area. In other implementations, it can also be set based on specific needs, such as based on the required number of third sampling points.

[0137] Step 403: determine a plurality of candidate third sampling points based on the grid.

[0138] like Figure 6B As shown, after constructing the grid of the window, multiple candidate third sampling points can be obtained based on the grid.

[0139] Step 404: determine the candidate third sampling point in the second area as the third sampling point.

[0140] In the above scheme, by constructing a window covering the second area, and then constructing a grid of the window based on a preset interval, multiple candidate third sampling points can be obtained, and the candidate third sampling points in the second area are determined as the third sampling points, and multiple third sampling points in the second area can be obtained. It is possible to subsequently obtain a sufficient amount of second surface temperature data based on the multiple third sampling points, and then establish a second surface temperature semivariogram, and finally determine the accurate air temperature of the sampling point to be measured in the target area.

[0141] In an example of the present application, a method for determining air temperature is also provided, such as Figure 7 As shown, the determining of an unbiased second air temperature semivariogram based on the first air temperature semivariogram, the first land surface temperature semivariogram and the second land surface temperature semivariogram comprises:

[0142] Step 501, determining a semivariogram ratio based on the first air temperature semivariogram and the first surface temperature semivariogram.

[0143] Specifically, the proportional function P of the semivariogram can be determined based on the following formula:

[0144]

[0145] Among them, γ east (h) is the first air temperature semivariogram, γ east，LST (h) is the first surface temperature semivariogram.

[0146] Step 501, determining an unbiased second air temperature semivariogram based on the semivariogram ratio and the second surface temperature semivariogram.

[0147] Specifically, the unbiased second air temperature semivariogram γ can be determined based on the following formula west，LST (h):

[0148] γ west，LST (h) = P·γ west，LST (h)

[0149] Among them, γ west，LST (h) is the second surface temperature semivariogram.

[0150] In an example of the present application, a method for determining air temperature is also provided, such as Figure 8 As shown, the method of determining the air temperature range of the third sampling point by linear regression includes:

[0151] Step 601: collect the land surface temperature data of the second sampling point by using a satellite to obtain third land surface temperature data of the second sampling point.

[0152] The remote sensing image of the surface temperature at each second sampling point is observed by satellite to obtain the third surface temperature data of each second sampling point.

[0153] Step 602: Establish a linear regression function based on the first sampling point, the first air temperature data and the first surface temperature data of the first sampling point, and the second sampling point and the second air temperature data and the third surface temperature data of the second sampling point, wherein the linear regression function represents the mapping relationship between the air temperature and the surface temperature in the target area.

[0154] The second surface temperature data of the third sampling point cannot directly represent the air temperature, and a simple linear transformation is required. First, a linear regression function is established using the first sampling point, the first air temperature data of the first sampling point and the first surface temperature data, the second sampling point, the second air temperature data of the second sampling point and the third surface temperature data. Since the surface temperature and air temperature in the same area are correlated, the linear regression function obtained by the above method can characterize the mapping relationship between the air temperature and the surface temperature in the target area.

[0155] Step 603: Determine the air temperature range of the third sampling point based on the third sampling point, the second ground surface temperature data of the third sampling point, and the linear regression function.

[0156] Substitute the second surface temperature data of the third sampling point into the linear regression function to obtain the air temperature data of the third sampling point. However, due to the uncertainty of remote sensing images collected by satellites, it is necessary to replace the air temperature data of the third sampling point with a 95% prediction interval to obtain the air temperature range of the third sampling point.

[0157] In an example of the present application, a method for determining air temperature is also provided, comprising:

[0158] Obtaining first air temperature data of a plurality of first sampling points in a first area, wherein the first area and the second area belong to the same target area but belong to different climate zones;

[0159] In the solution of the present application, the second air temperature data of the second area is not necessary. If there is no meteorological station in the second area (that is, there is no second sampling point) and no second air temperature data can be collected, the solution of the present application can still determine the accurate air temperature at any location in the target area.

[0160] establishing a first air temperature semivariogram based on the first air temperature data;

[0161] Collecting the surface temperature data of the first sampling point by using a satellite to obtain the first surface temperature data of the first sampling point;

[0162] Establishing a first land surface temperature semivariogram based on the first land surface temperature data;

[0163] Establishing a plurality of third sampling points in the second area, and collecting surface temperature data of the third sampling points by using a satellite to obtain second surface temperature data of the third sampling points;

[0164] Establishing a second land surface temperature semivariogram based on the second land surface temperature data;

[0165] determining an unbiased second air temperature semivariogram based on the first air temperature semivariogram, the first land surface temperature semivariogram, and the second land surface temperature semivariogram;

[0166] Determining the air temperature range of the third sampling point using linear regression;

[0167] When establishing a linear regression function, when there is no second sampling point and the second air temperature data of the second sampling point, a linear regression function can be established only based on the first sampling point, the first air temperature data of the first sampling point and the first surface temperature data. Similarly, a linear regression function that characterizes the mapping relationship between the air temperature and the surface temperature in the target area can be established.

[0168] The air temperature of the sampling point to be measured is determined based on the sampling point to be measured, the second sampling point, the second air temperature data of the second sampling point, the third sampling point, the air temperature range of the third sampling point and the unbiased second air temperature semivariogram.

[0169] In order to implement the above-mentioned air temperature determination method, as Fig. 9 As shown, an example of the present application provides an air temperature determination device, comprising:

[0170] The acquisition module 10 is used to obtain first air temperature data of a plurality of first sampling points in a first area and second air temperature data of a plurality of second sampling points in a second area, wherein the first area and the second area belong to the same target area but belong to different climate zones, and the number of the second sampling points is less than the target ratio of the sum of the first sampling points and the second sampling points; collect the surface temperature data of the first sampling points by using a satellite to obtain the first surface temperature data of the first sampling points; and establish a plurality of third sampling points in the second area, and collect the surface temperature data of the third sampling points by using a satellite to obtain the second surface temperature data of the third sampling points;

[0171] A calculation module 20 is used to establish a first air temperature semivariogram based on the first air temperature data; to establish a first surface temperature semivariogram based on the first surface temperature data; to establish a second surface temperature semivariogram based on the second surface temperature data; to determine an unbiased second air temperature semivariogram based on the first air temperature semivariogram, the first surface temperature semivariogram and the second surface temperature semivariogram; to determine the air temperature range of the third sampling point using linear regression; and to determine the air temperature of the sampling point to be measured based on the sampling point to be measured, the second sampling point, the second air temperature data of the second sampling point, the third sampling point, the air temperature range of the third sampling point and the unbiased second air temperature semivariogram.

[0172] The acquisition module 10 is further used to acquire the surface temperature data of the second sampling point by using a satellite to obtain the third surface temperature data of the second sampling point;

[0173] The calculation module 20 is further used to establish a linear regression function based on the first sampling point, the first air temperature data and the first ground surface temperature data of the first sampling point, the second sampling point and the second air temperature data and the third ground surface temperature data of the second sampling point, wherein the linear regression function represents the mapping relationship between the air temperature and the ground surface temperature in the target area; and determine the air temperature range of the third sampling point based on the third sampling point, the second ground surface temperature data of the third sampling point and the linear regression function.

[0174] Among them, the calculation module 20 is also used to traverse multiple preset range values; determine at least one first sampling point pair whose distance is less than or equal to the current preset range value; determine the function value corresponding to the current preset range value based on the at least one first sampling point pair of the current preset range value and the first air temperature data of each first sampling point; after the traversal is completed, obtain the function value corresponding to each preset range value; and determine the first air temperature semivariogram function based on the multiple preset range values ​​and the function value corresponding to each preset range value.

[0175] Among them, the calculation module 20 is also used to determine the square of the air temperature difference of each first sampling point pair according to the first air temperature data of each first sampling point; and sum the squares of the air temperature differences of all first sampling point pairs, and divide the sum by twice the total number of first sampling point pairs to obtain the function value corresponding to the current preset range value.

[0176] Among them, it also includes:

[0177] The processing module 30 is configured to construct a window covering the second area; construct a grid of the window based on a preset interval; determine a plurality of candidate third sampling points based on the grid; and determine the candidate third sampling point in the second area as the third sampling point.

[0178] Among them, the calculation module 20 is also used to determine the semivariogram ratio based on the first air temperature semivariogram and the first surface temperature semivariogram; and determine the unbiased second air temperature semivariogram based on the semivariogram ratio and the second surface temperature semivariogram.

[0179] According to an embodiment of the present disclosure, the present disclosure also provides an electronic device and a readable storage medium.

[0180] Fig.10 A schematic block diagram of an example electronic device 700 that can be used to implement an embodiment of the present disclosure is shown. The electronic device is intended to represent various forms of digital computers, such as laptop computers, desktop computers, workstations, personal digital assistants, servers, blade servers, mainframe computers, and other suitable computers. The electronic device can also represent various forms of mobile devices, such as personal digital processing, cellular phones, smart phones, wearable devices, and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely examples and are not intended to limit the implementation of the present disclosure described and / or required herein.

[0181] like Fig.10 As shown, the device 700 includes a computing unit 701, which can perform various appropriate actions and processes according to a computer program stored in a read-only memory (ROM) 702 or a computer program loaded from a storage unit 708 into a random access memory (RAM) 703. In the RAM 703, various programs and data required for the operation of the device 700 can also be stored. The computing unit 701, the ROM 702, and the RAM 703 are connected to each other via a bus 704. An input / output (I / O) interface 705 is also connected to the bus 704.

[0182] A number of components in the device 700 are connected to the I / O interface 705, including: an input unit 706, such as a keyboard, a mouse, etc.; an output unit 707, such as various types of displays, speakers, etc.; a storage unit 708, such as a disk, an optical disk, etc.; and a communication unit 709, such as a network card, a modem, a wireless communication transceiver, etc. The communication unit 709 allows the device 700 to exchange information / data with other devices through a computer network such as the Internet and / or various telecommunication networks.

[0183] The computing unit 701 may be a variety of general and / or special processing components with processing and computing capabilities. Some examples of the computing unit 701 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various dedicated artificial intelligence (AI) computing chips, various computing units running machine learning model algorithms, digital signal processors (DSPs), and any appropriate processors, controllers, microcontrollers, etc. The computing unit 701 performs the various methods and processes described above, such as the air temperature determination method. For example, in some embodiments, the air temperature determination method may be implemented as a computer software program, which is tangibly contained in a machine-readable medium, such as a storage unit 708. In some embodiments, part or all of the computer program may be loaded and / or installed on the device 700 via the ROM 702 and / or the communication unit 709. When the computer program is loaded into the RAM 703 and executed by the computing unit 701, one or more steps of the air temperature determination method described above may be performed. Alternatively, in other embodiments, the computing unit 701 may be configured to perform the air temperature determination method in any other appropriate manner (e.g., by means of firmware).

[0184] Various implementations of the systems and techniques described above herein can be implemented in digital electronic circuit systems, integrated circuit systems, field programmable gate arrays (FPGAs), application specific integrated circuits (ASICs), application specific standard products (ASSPs), integrated systems on a chip (SOCs), complex programmable logic devices (CPLDs), computer hardware, firmware, software, and / or combinations thereof. These various implementations can include: being implemented in one or more computer programs that can be executed and / or interpreted on a programmable system including at least one programmable processor, which can be a special purpose or general purpose programmable processor that can receive data and instructions from a storage system, at least one input device, and at least one output device, and transmit data and instructions to the storage system, the at least one input device, and the at least one output device.

[0185] The program code for implementing the method of the present disclosure may be written in any combination of one or more programming languages. These program codes may be provided to a processor or controller of a general-purpose computer, a special-purpose computer, or other programmable data processing device, so that the program code, when executed by the processor or controller, enables the functions / operations specified in the flow chart and / or block diagram to be implemented. The program code may be executed entirely on the machine, partially on the machine, partially on the machine and partially on a remote machine as a stand-alone software package, or entirely on a remote machine or server.

[0186] In the context of the present disclosure, a machine-readable medium may be a tangible medium that may contain or store a program for use by or in conjunction with an instruction execution system, device, or equipment. A machine-readable medium may be a machine-readable signal medium or a machine-readable storage medium. A machine-readable medium may include, but is not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, device, or equipment, or any suitable combination of the foregoing. A more specific example of a machine-readable storage medium may include an electrical connection based on one or more lines, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing.

[0187] To provide interaction with a user, the systems and techniques described herein can be implemented on a computer having: a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to the user; and a keyboard and pointing device (e.g., a mouse or trackball) through which the user can provide input to the computer. Other types of devices can also be used to provide interaction with the user; for example, the feedback provided to the user can be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user can be received in any form (including acoustic input, voice input, or tactile input).

[0188] The systems and techniques described herein may be implemented in a computing system that includes back-end components (e.g., as a data server), or a computing system that includes middleware components (e.g., an application server), or a computing system that includes front-end components (e.g., a user computer with a graphical user interface or a web browser through which a user can interact with implementations of the systems and techniques described herein), or a computing system that includes any combination of such back-end components, middleware components, or front-end components. The components of the system may be interconnected by any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include: a local area network (LAN), a wide area network (WAN), and the Internet.

[0189] A computer system may include a client and a server. The client and the server are generally remote from each other and usually interact through a communication network. The relationship of client and server is generated by computer programs running on respective computers and having a client-server relationship with each other. The server may be a cloud server, a server of a distributed system, or a server combined with a blockchain.

[0190] It should be understood that the various forms of processes shown above can be used to reorder, add or delete steps. For example, the steps recorded in this disclosure can be executed in parallel, sequentially or in different orders, as long as the desired results of the technical solutions disclosed in this disclosure can be achieved, and this document does not limit this.

[0191] In addition, the terms "first" and "second" are used for descriptive purposes only and should not be understood as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include at least one of the features. In the description of the present disclosure, the meaning of "plurality" is two or more, unless otherwise clearly and specifically defined.

[0192] The above is only a specific embodiment of the present disclosure, but the protection scope of the present disclosure is not limited thereto. Any person skilled in the art who is familiar with the technical field can easily think of changes or substitutions within the technical scope disclosed in the present disclosure, which should be included in the protection scope of the present disclosure. Therefore, the protection scope of the present disclosure should be based on the protection scope of the claims.

Claims

1. A method for determining air temperature, the method comprising: Obtaining first air temperature data of a plurality of first sampling points in a first region and second air temperature data of a plurality of second sampling points in a second region, wherein the first region and the second region belong to the same target region but belong to different climate zones, and the number of the second sampling points is less than a target ratio of the sum of the first sampling points and the second sampling points; establishing a first air temperature semivariogram based on the first air temperature data; Collecting the surface temperature data of the first sampling point by using a satellite to obtain the first surface temperature data of the first sampling point; Establishing a first land surface temperature semivariogram based on the first land surface temperature data; Establishing a plurality of third sampling points in the second area, and collecting surface temperature data of the third sampling points by using a satellite to obtain second surface temperature data of the third sampling points; Establishing a second land surface temperature semivariogram based on the second land surface temperature data; determining an unbiased second air temperature semivariogram based on the first air temperature semivariogram, the first land surface temperature semivariogram, and the second land surface temperature semivariogram; Determining the air temperature range of the third sampling point using linear regression; Determine the air temperature of the sampling point to be measured based on the sampling point to be measured, the second sampling point, the second air temperature data of the second sampling point, the third sampling point, the air temperature range of the third sampling point and the unbiased second air temperature semivariogram; Wherein, determining an unbiased second air temperature semivariogram based on the first air temperature semivariogram, the first land surface temperature semivariogram, and the second land surface temperature semivariogram comprises: determining a semivariogram scale based on the first air temperature semivariogram and the first land surface temperature semivariogram; determining an unbiased second air temperature semivariogram based on the semivariogram ratio and the second land surface temperature semivariogram; Wherein, determining the air temperature range of the third sampling point by using linear regression includes: collecting the land surface temperature data of the second sampling point by using a satellite to obtain third land surface temperature data of the second sampling point; Establishing a linear regression function based on the first sampling point, the first air temperature data and the first ground surface temperature data of the first sampling point, the second sampling point and the second air temperature data and the third ground surface temperature data of the second sampling point, wherein the linear regression function represents the mapping relationship between the air temperature and the ground surface temperature in the target area; determining an air temperature range of the third sampling point based on the third sampling point, the second ground surface temperature data of the third sampling point, and the linear regression function; Wherein, determining the air temperature of the sampling point to be measured based on the sampling point to be measured, the second sampling point, the second air temperature data of the second sampling point, the third sampling point, the air temperature range of the third sampling point and the unbiased second air temperature semivariogram comprises: The sampling point to be measured, the second sampling point, the second air temperature data of the second sampling point, the third sampling point, the air temperature range of the third sampling point and the unbiased second air temperature semivariogram are used as input data, input into matlab using BMELIB2.0c to perform Bayesian maximum entropy calculation, and the air temperature of the sampling point to be measured is determined through calculation instructions.

2. The method according to claim 1, wherein establishing a first air temperature semivariogram based on the first air temperature data comprises: Traverse multiple preset range values; Determine at least one first sampling point pair whose distance is less than or equal to the current preset range value; Determine a function value corresponding to the current preset range value according to at least one first sampling point pair of the current preset range value and first air temperature data of each first sampling point; After the traversal is completed, the function value corresponding to each preset range value is obtained; A first air temperature semivariogram is determined according to the plurality of preset range values ​​and a function value corresponding to each preset range value.

3. The method according to claim 2, wherein determining the function value corresponding to the current preset range value according to at least one first sampling point pair of the current preset range value and the first air temperature data of each first sampling point comprises: Determine the square of the air temperature difference of each first sampling point pair according to the first air temperature data of each first sampling point; The squares of the air temperature differences of all first sampling point pairs are summed, and the sum result is divided by twice the total number of first sampling point pairs to obtain the function value corresponding to the current preset range value.

4. The method according to claim 1, wherein establishing a plurality of third sampling points in the second area comprises: constructing a window covering the second area; A grid of windows is constructed based on a preset interval; determining a plurality of candidate third sampling points based on the grid; A candidate third sampling point in the second area is determined as the third sampling point.

5. An air temperature determination device, the device comprising: A collection module is used to obtain first air temperature data of a plurality of first sampling points in a first area and second air temperature data of a plurality of second sampling points in a second area, wherein the first area and the second area belong to the same target area but belong to different climate zones, and the number of the second sampling points is less than a target ratio of the sum of the first sampling points and the second sampling points; and to collect the surface temperature data of the first sampling points by using a satellite to obtain the first surface temperature data of the first sampling points; and establishing a plurality of third sampling points in the second area, and collecting surface temperature data of the third sampling points by using satellites to obtain second surface temperature data of the third sampling points; A calculation module, used for establishing a first air temperature semivariogram based on the first air temperature data; Establishing a first land surface temperature semivariogram based on the first land surface temperature data; Establish a second surface temperature semivariogram based on the second surface temperature data; determine an unbiased second air temperature semivariogram based on the first air temperature semivariogram, the first surface temperature semivariogram and the second surface temperature semivariogram; determine the air temperature range of the third sampling point using linear regression; and determine the air temperature of the sampling point to be measured based on the sampling point to be measured, the second sampling point, the second air temperature data of the second sampling point, the third sampling point, the air temperature range of the third sampling point and the unbiased second air temperature semivariogram; The calculation module is further used to determine a semivariogram ratio based on the first air temperature semivariogram and the first land surface temperature semivariogram; and determine an unbiased second air temperature semivariogram based on the semivariogram ratio and the second land surface temperature semivariogram; The acquisition module is further used to acquire the surface temperature data of the second sampling point by using a satellite to obtain the third surface temperature data of the second sampling point; The calculation module is further used to establish a linear regression function based on the first sampling point, the first air temperature data and the first ground surface temperature data of the first sampling point, and the second sampling point and the second air temperature data and the third ground surface temperature data of the second sampling point, wherein the linear regression function represents the mapping relationship between the air temperature and the ground surface temperature in the target area; and determining an air temperature range of the third sampling point based on the third sampling point, the second ground surface temperature data of the third sampling point and the linear regression function; The calculation module is also used to take the sampling point to be measured, the second sampling point, the second air temperature data of the second sampling point, the third sampling point, the air temperature range of the third sampling point and the unbiased second air temperature semivariogram as input data, input them into matlab using BMELIB2.0c to perform Bayesian maximum entropy calculation, and determine the air temperature of the sampling point to be measured through calculation instructions.

6. An electronic device comprising: processor and input-output systems; The input-output system obtains first air temperature data of a plurality of first sampling points in a first region and second air temperature data of a plurality of second sampling points in a second region, the first region and the second region belong to the same target region but belong to different climate zones, and the number of the second sampling points is less than a target ratio of the sum of the first sampling points and the second sampling points; The processor establishes a first air temperature semivariogram based on the first air temperature data; The input-output system collects the surface temperature data of the first sampling point by using a satellite to obtain the first surface temperature data of the first sampling point; The processor establishes a first land surface temperature semivariogram based on the first land surface temperature data; The input-output system establishes a plurality of third sampling points in the second area, and uses satellites to collect surface temperature data of the third sampling points to obtain second surface temperature data of the third sampling points; The processor establishes a second land surface temperature semivariogram based on the second land surface temperature data; The processor determines an unbiased second air temperature semivariogram based on the first air temperature semivariogram, the first land surface temperature semivariogram, and the second land surface temperature semivariogram; The processor determines the air temperature range of the third sampling point using linear regression; The processor determines the air temperature of the sampling point to be measured based on the sampling point to be measured, the second sampling point, the second air temperature data of the second sampling point, the third sampling point, the air temperature range of the third sampling point and the unbiased second air temperature semivariogram; The processor determines a semivariogram scale based on the first air temperature semivariogram and the first land surface temperature semivariogram; The processor determines an unbiased second air temperature semivariogram based on the semivariogram scale and the second land surface temperature semivariogram; The input-output system collects the surface temperature data of the second sampling point by using a satellite to obtain the third surface temperature data of the second sampling point; The processor establishes a linear regression function based on the first sampling point, the first air temperature data and the first ground surface temperature data of the first sampling point, the second sampling point and the second air temperature data and the third ground surface temperature data of the second sampling point, wherein the linear regression function represents the mapping relationship between the air temperature and the ground surface temperature in the target area; The processor determines an air temperature range of the third sampling point based on the third sampling point, the second ground surface temperature data of the third sampling point, and the linear regression function; The processor takes the sampling point to be measured, the second sampling point, the second air temperature data of the second sampling point, the third sampling point, the air temperature range of the third sampling point and the unbiased second air temperature semivariogram as input data, inputs them into matlab using BMELIB2.0c to perform Bayesian maximum entropy calculation, and determines the air temperature of the sampling point to be measured through calculation instructions.

7. A computer-readable storage medium storing a computer program, the computer program being used to obtain first air temperature data of a plurality of first sampling points in a first region and second air temperature data of a plurality of second sampling points in a second region, the first region and the second region being both target regions but belonging to different climate zones, and the number of the second sampling points being less than a target ratio of the sum of the first sampling points and the second sampling points; Establishing a first air temperature semivariogram based on the first air temperature data; using a satellite to collect the surface temperature data of the first sampling point to obtain the first surface temperature data of the first sampling point; Establishing a first land surface temperature semivariogram based on the first land surface temperature data; Establishing a plurality of third sampling points in the second area, and collecting surface temperature data of the third sampling points by using a satellite to obtain second surface temperature data of the third sampling points; Establish a second surface temperature semivariogram based on the second surface temperature data; determine an unbiased second air temperature semivariogram based on the first air temperature semivariogram, the first surface temperature semivariogram and the second surface temperature semivariogram; determine the air temperature range of the third sampling point using linear regression; determine the air temperature of the sampling point to be measured based on the sampling point to be measured, the second sampling point, the second air temperature data of the second sampling point, the third sampling point, the air temperature range of the third sampling point and the unbiased second air temperature semivariogram; The computer program is further configured to determine a semivariogram ratio based on the first air temperature semivariogram and the first land surface temperature semivariogram; determine an unbiased second air temperature semivariogram based on the semivariogram ratio and the second land surface temperature semivariogram; The computer program is also used to collect the surface temperature data of the second sampling point by using a satellite to obtain the third surface temperature data of the second sampling point; based on the first sampling point, the first air temperature data and the first surface temperature data of the first sampling point, the second sampling point and the second air temperature data and the third surface temperature data of the second sampling point, a linear regression function is established, wherein the linear regression function represents the mapping relationship between the air temperature and the surface temperature in the target area; determining an air temperature range of the third sampling point based on the third sampling point, the second ground surface temperature data of the third sampling point, and the linear regression function; The computer program is also used to use the sampling point to be measured, the second sampling point, the second air temperature data of the second sampling point, the third sampling point, the air temperature range of the third sampling point and the unbiased second air temperature semivariogram as input data, input them into matlab using BMELIB2.0c to perform Bayesian maximum entropy calculation, and determine the air temperature of the sampling point to be measured through calculation instructions.

Images