CN114417579B - Different unit grid distance conversion method for troposphere electric wave environment numerical mode result - Google Patents

Different unit grid distance conversion method for troposphere electric wave environment numerical mode result Download PDF

Info

Publication number
CN114417579B
CN114417579B CN202111674069.9A CN202111674069A CN114417579B CN 114417579 B CN114417579 B CN 114417579B CN 202111674069 A CN202111674069 A CN 202111674069A CN 114417579 B CN114417579 B CN 114417579B
Authority
CN
China
Prior art keywords
grid
longitude
latitude
equal
unit
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202111674069.9A
Other languages
Chinese (zh)
Other versions
CN114417579A (en
Inventor
张玉生
郭相明
赵振维
李清亮
康士峰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Institute of Radio Wave Propagation CETC 22 Research Institute
Original Assignee
China Institute of Radio Wave Propagation CETC 22 Research Institute
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Institute of Radio Wave Propagation CETC 22 Research Institute filed Critical China Institute of Radio Wave Propagation CETC 22 Research Institute
Priority to CN202111674069.9A priority Critical patent/CN114417579B/en
Publication of CN114417579A publication Critical patent/CN114417579A/en
Application granted granted Critical
Publication of CN114417579B publication Critical patent/CN114417579B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

The invention discloses a method for converting different unit grid distances of troposphere electric wave environment numerical mode results, which is characterized in that on the basis of acquiring more fine troposphere electric wave environment parameter space-time grid data based on a mesoscale atmospheric numerical mode, directly output grid data set in a grid distance kilometer unit are taken as a source, and based on the steps of grid reconstruction, acquiring the area ratio of a single grid with equal longitude and latitude units formed by original kilometer unit grids and acquiring the component of the troposphere electric wave environment parameter of the original kilometer unit grids in the equal longitude and latitude units formed by the same longitude and latitude units, and the like, the equal grid distance data set in the grid distance longitude and latitude units are converted through a certain conversion relation.

Description

Different unit grid distance conversion method for troposphere electric wave environment numerical mode result
Technical Field
The invention belongs to the technical field of meteorological numerical simulation and radio, and particularly relates to a method for converting different unit grid distances of troposphere electric wave environment numerical mode results in the field.
Background
In modern electronic information technology, tropospheric radio wave environmental characteristics are important factors affecting effective performance of electronic information systems such as radar, communication, navigation and the like. The troposphere electric wave environmental characteristics greatly depend on the acquisition and statistics of troposphere electric wave environmental parameters, the acquisition and statistics of the parameters in the past can only depend on sparse and limited meteorological ground and sounding observation stations in an area for data acquisition and statistics, and for a vast ocean area, meteorological sounding stations have no or rare sites and cannot collect data. With the development of the mesoscale atmospheric numerical mode, finer troposphere electric wave environmental parameter space-time grid data can be obtained, and therefore the problems of acquisition of regional troposphere electric wave environmental data and statistical characteristic research can be solved.
Mesoscale atmospheric numerical models have evolved in the 80 s of the 20 th century, and in the 90 s some mesoscale models and simulation systems have been widely used worldwide, such as the sea-gas coupled mesoscale prediction system (COAMPS) in the united states, the british gas phase mesoscale business model, the canadian mesoscale compressible common model, the french mesoscale non-static model, the japanese regional spectral model, and so on. The wide mesoscale atmospheric model spread at home and abroad is the fifth generation mesoscale non-hydrostatic model MM5 and meteorological research and forecasting model (WRF) developed by the national atmospheric research center (NCAR) and the Pennsylvania State university in America in cooperation, and the WRF model is a numerical meteorological prediction and atmospheric simulation system for research and business applications, developed on the basis of MM5 and developed by the national atmospheric research center (NCAR), the national environmental prediction center (NCEP), the National Oceanic and Atmospheric Administration (NOAA), the national aerospace administration (NASA), the air force weather administration (AFAFA) and many organizations such as the Association and the WA. WRF has a flexible programming architecture and finds wide application in research, business, and teaching. The correlation of the mesoscale atmospheric numerical mode is described in detail in many publications and is not described in detail here.
In the process of acquiring finer tropospheric electric wave environment parameter space-time grid data based on the mesoscale atmospheric numerical mode, a user is required to perform research area presetting (the size of the research area is determined by the user as required) and horizontal space grid distance setting, the grid distance setting is set according to an equal kilometer unit (a fixed value is set by the user as required in advance, and d km is assumed here), so that the acquired tropospheric electric wave environment parameter horizontal space grid distance unit is km, in other words, the regional tropospheric electric wave environment data acquired based on the mesoscale atmospheric numerical mode are all formed by a plurality of grid data which are formed by taking d km (a fixed value set by the user in advance in the mode) as equal intervals, if the grid data are converted into a longitude and latitude unit, the grid data are not equidistant according to the longitude and latitude unit, and the projection mode of the numerical mode and the difference between north and south, and the latitude are in great relation.
The mode obtains the longitude and latitude unit grid data with unequal intervals or the kilometer unit grid data with equal intervals in the research area, and in practical application, the type data can be conveniently used in a small research area or at a single point. However, this type of data is inconvenient when studying the environmental characteristics of tropospheric electric waves in a large area, particularly in a global area, and it is more convenient to use mesh data of equal-spaced latitude and longitude units, for example, the mesh data used in the relevant standards in the radio wave propagation standard of the international telecommunication union (ITU-R) is reanalysis mesh data of equal-spaced latitude and longitude units (the data is the direct output result of a global atmospheric numerical mode), while the present example relates to an atmospheric mesoscale numerical mode, which considers the projection and coordinate modes, the representativeness of local climate and the depiction of the fine structure of the atmosphere, in order to obtain a more accurate and finer result than the previous data, but the output result of the mode is an equal-kilometer unit mesh, and the study of the environmental characteristics of global or large-area electric waves needs to be converted into mesh data of equal-spaced latitude and longitude units, but people rarely relate to a specific conversion method for different unit mesh data of the mode data.
Disclosure of Invention
The invention aims to provide a method for converting different unit grid distances of troposphere electric wave environment numerical mode results.
The invention adopts the following technical scheme:
in a method for transforming different unit cell pitches as a result of numerical patterns of tropospheric radiowave environments, the improvement comprising the steps of:
step 1, in a numerical mode, a user sets a research area (the research area is user-defined) and a horizontal grid distance d km in advance, wherein d is less than 60;
step 2, obtaining an equal kilometer unit grid data file which is output by a user in a mesoscale atmospheric numerical mode and takes d km as equal space, wherein the grid data file comprises longitude and latitude information of unequal space corresponding to each grid data in a research area and troposphere electric wave environmental parameter information to be researched;
step 3, determining four vertexes of a final equal-longitude and latitude unit grid to be converted in the research area according to a specific equal-kilometer unit grid distance d km of a troposphere electric wave environment numerical mode result and a grid distance c of a specific equal-longitude and latitude unit to be converted;
step 4, if the boundary of the final equal longitude and latitude unit grid to be converted exceeds the regional boundary of the original equal longitude and latitude unit grid in the research region, the grid point is not converted and assigned, the conversion work of the corresponding equal grid distance grid is terminated, the conversion work of the rest grid points in the research region is further carried out, the step 3 is repeatedly carried out, and all regions covered by the final equal longitude and latitude unit grid obtained later do not contain the grid which is terminated;
if the boundary of the final equal-longitude-latitude unit single grid to be converted does not exceed the regional boundary of the original equal-longitude-latitude unit grid, performing step 5;
step 5, acquiring the sub-area of the converted final equal longitude and latitude unit single grid formed by the original equal kilometer unit grids according to the four vertexes of the determined equal longitude and latitude unit grid to be converted;
sub-area S of original equal kilometer unit grid in final equal longitude and latitude unit single grid Score of origin The calculation is as follows:
S score of origin ＝d Original weft dividing ×d Original meridian point (1)
Wherein, d Original weft division The latitude distance, d, of the final equal longitude and latitude unit single grid of the original equal kilometer unit grid Original meridian branch The longitude distance of the area of the original equal-kilometer unit grid occupying the final equal-longitude and latitude unit single grid is used;
step 6, according to the obtained sub-area S Score of origin Obtaining the area ratio S of the single grid which is formed by the original grids of the same kilometer unit and forms the final unit of the same longitude and latitude Ratio of occupation of ；
The area ratio S of original equal kilometer unit single grids forming the final equal longitude and latitude unit grid is calculated based on the following formula (2) Ratio of occupation of ：
S Ratio of occupation of ＝(S Score of origin /S Finally, the product is processed )×100％ (2)
Wherein S is Finally, the product is processed The area of the single grid is the final unit with equal latitude and longitude;
obtaining S based on the following formula (3) Finally, the product is processed ：
S Finally, the product is processed ＝d Weft yarn ×d Warp beam (3)
Wherein d is Weft yarn The latitude distance of a single grid in final equal latitude and longitude units, d Warp beam The longitude distance of a single grid is the final equal longitude and latitude unit;
step 7, acquiring troposphere electric wave environment parameters W of the original kilometer unit grids Original source Single grid ratio S of final equal longitude and latitude unit Ratio of occupation of Component W of Account for ；
Reading the troposphere electric wave environment parameter value corresponding to each original kilometer unit grid according to the longitude and latitude from the grid data file as required, and acquiring the troposphere electric wave environment parameter W of the original kilometer unit grid according to the following formula (4) Original source Single grid ratio S of final equal longitude and latitude unit Ratio of the ingredients Component W of Account for ：
W Account for ＝W Original source ×S Ratio of occupation of (4)
Step 8, the component W obtained in the step 7 is processed Account for Adding the obtained final converted grid value W of one of the equal grid distance troposphere electric wave environmental parameters set in the unit of grid distance longitude and latitude Finally, the product is processed 。
Further, the number of grids in the research area set by the user is N multiplied by m, and each grid has N electric wave environment parameters; repeating the steps 7-8 to obtain other troposphere electric wave environment parameters of the final equal grid distance grid; and repeating the steps 3-9 to obtain the troposphere electric wave environmental parameters of other final equal grid distance grids in the research area.
The invention has the beneficial effects that:
although the tropospheric electric wave environment parameters acquired based on the mesoscale atmospheric numerical mode are finer, they are all composed of a lot of grid data formed at equal intervals by d km (fixed value set in advance by the user in the mode), if converted into latitude and longitude units, the grid data are not equidistant for the latitude and longitude units, and the data are inconvenient to use in large areas, especially in the global scope.
The invention discloses a method for converting different unit grid distances of troposphere electric wave environment numerical value mode results, which is based on the steps of grid reconstruction, obtaining the area ratio of a single grid with equal longitude and latitude units formed by original kilometer unit grids, obtaining the component of the troposphere electric wave environment parameters of the original kilometer unit grids in the single grid ratio with equal longitude and latitude units formed by the convection layer electric wave environment parameters and the like, and further processing the grid data output by an atmospheric mesoscale numerical value mode through a certain conversion relation, so that the kilometer unit grid data are converted into equidistant longitude and latitude unit grid data, and the method is convenient to practical application. The method can convert the longitude and latitude unit grid data which are not equidistant into the longitude and latitude unit grid data which are equidistant, and is convenient for the practical application of researching the troposphere electric wave environmental characteristics in a large area, particularly in a global range.
Drawings
FIG. 1 is a schematic diagram of grid spacing characteristics of grid data for a mesoscale numerical mode;
FIG. 2 is a diagram of a mesh data file format for tropospheric electric wave environment numerical pattern results;
FIG. 3 is a schematic diagram of grid conversion from equal kilometer units of grid distance to final equal longitude and latitude units of tropospheric electric wave environment numerical pattern results;
fig. 4 is a diagram showing an example of grid conversion from equal kilometer units of grid distance to final equal longitude and latitude units of a tropospheric electric wave environment numerical pattern result;
fig. 5 is a diagram showing an example of a mesh data file format of the tropospheric electric wave environment numerical pattern result.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention will be described in further detail below with reference to the accompanying drawings and examples. It should be understood that the specific embodiments described herein are merely illustrative of the invention and do not limit the invention.
The method comprises the steps of obtaining area ratio of single grids with equal longitude and latitude units formed by original kilometer unit grids, obtaining component of single grid ratio with equal longitude and latitude units formed by convection layer electric wave environment parameters of original kilometer unit grids and the like, and converting into equal grid distance data set in grid distance longitude and latitude units through a certain conversion relation on the basis of obtaining finer time-space grid data of convection layer electric wave environment parameters based on a mesoscale atmospheric numerical mode, taking directly output grid data set in grid distance kilometer units as a source, and based on the steps of grid reconstruction, obtaining the area ratio of single grids with equal longitude and latitude units formed by original kilometer unit grids and the like.
As is well known, tropospheric radio wave environment data in a research area acquired based on a mesoscale atmospheric numerical mode are all a lot of grid data formed by taking d km (a fixed value preset by a user in the mode) as an equal interval, and if the grid data are converted into latitude and longitude units, the grid data are not equidistant for the latitude and longitude units, because the projection mode of the numerical mode and the difference between north and south latitudes have great relations, for example, the latitude grid distance of a data grid point becomes larger along with the decreasing of latitude, and the longitude grid distance is basically unchanged along with the decreasing of longitude, as shown in fig. 1, each grid distance is a fixed number d km, but the difference between the latitudes of adjacent grid points is north latitude 4-north latitude 3< north latitude 3-north latitude 2< north latitude 2-north latitude 1, and obviously, the grid distances of latitude units are not equidistant. The use of the equidistant longitude and latitude unit grid data is more convenient when studying the troposphere electric wave environmental characteristics in a large area, particularly in a global range.
The method of the invention is that a user sets a horizontal grid distance d km in advance in a numerical mode, grid data with the d km (a fixed value set in advance by the user in the mode) being equal intervals output by a mesoscale atmospheric numerical mode are obtained, a grid data file contains unequal interval longitude and latitude information corresponding to each grid data in a research area and troposphere electric wave environment parameter information to be researched, and four vertexes of a certain grid with equal longitude and latitude to be converted are determined according to a specific equal kilometer unit grid distance (if the grid distance d km) and a specific equal longitude and latitude unit grid distance (if the grid distance c DEG) to be converted of a troposphere electric wave environment numerical mode result. After 4 vertexes of a grid with equal longitude and latitude units are determined, the grid is determined, then N electric wave environment parameters of the grid are determined most importantly, and the calculation method of each electric wave environment parameter is obtained by weighting and calculating the electric wave environment parameters of original several grids with equal kilometers according to the area size of the original several grids with equal kilometers forming the grid and the area formed by the electric wave environment parameters of the original several grids with equal kilometers forming the grid correspondingly. And analogizing the solving methods of other grid electric wave environment parameters in the research area in turn.
In particular, the method of the present invention is suitable for the case where the grid distance c > d/110, that is, the grid distance of the specific equal longitude and latitude unit to be converted is greater than the grid distance of the equal kilometer unit of the troposphere electric wave environment numerical value mode result, wherein 110 is a specific numerical value of how many kilometer units the longitude and latitude unit is converted into, c is the grid distance of the equal longitude and latitude unit, and d is the grid distance of the specific equal kilometer unit.
The method specifically comprises the following steps:
step 1, a research area (user-defined in the research area) and a horizontal grid distance d km are preset by a user in a numerical mode, then the research area and the horizontal grid distance are fixed values, the grid distance characteristics of grid data of a mesoscale numerical mode are reflected in fig. 1, and all grid distances in the research area are equal and are all equal to d km. Where d is less than 60, the value is too large because of non-uniform atmospheric levels, and the converted data is prone to large errors;
step 2, a user acquires the grid data with equal kilometers and equal intervals, wherein d km (a fixed value preset by the user in the mode) output by a mesoscale atmospheric numerical mode is equal kilometers, a grid data file contains unequal interval longitude and latitude information corresponding to each grid data in a research area and troposphere electric wave environment parameter information to be researched, the data are arranged in sequence as shown in fig. 2, latitude 1 and longitude 1 represent the longitude and latitude corresponding to the first initial grid (line 1 and column 1) of a numerical mode output result, latitude N and longitude m represent the longitude and latitude corresponding to the last grid (line N and column m) of the numerical mode output result, the number of grids is N multiplied by m, and each grid has N electric wave environment parameters; it should be noted that each grid in the grid data file may include a plurality of radio wave environment parameters, and the number of parameters and the data format are not limited, as long as the grid data file can store and acquire the longitude and latitude of the corresponding grid point and the radio wave environment parameters.
And 3, determining four vertexes of the final equal-longitude and latitude unit grid to be converted in the research area according to the specific equal-kilometer unit grid distance (grid distance d km) of the troposphere electric wave environment numerical value mode result and the grid distance (grid distance c DEG) of the specific equal-longitude and latitude unit to be converted.
The grid corresponding to latitude 1 and longitude 1 is the 1 st grid, identified as W in fig. 3. If the grid conversion from the grid of the equal kilometer unit to the grid of the equal latitude and longitude unit is performed in the grid in the figure, people generally use the grid for convenient processing, the starting longitude and latitude (also representing the longitude and latitude of the grid) at the lower left corner of the first grid point is generally an integer, the position of the starting longitude and latitude is within the coverage range of the original data area, the grid distance c ° of the grid of the equal latitude and longitude unit is also regular, generally, the starting longitude and latitude (also representing the longitude and latitude of the grid point) at the lower left corner of the first grid point is 0.25 ° or 0.5 ° or 0.75 ° or 1 ° and the like, therefore, the longitude and latitude of the other three vertexes of the grid point are also regular, and by analogy, the starting longitude and latitude (also representing the longitude and latitude of the grid point) at the lower left corner of the other grid point is also regular, for example, the grid of the longitude and latitude unit is to be converted into the grid of the longitude and latitude unit such as ABCD surrounded by 4 dotted lines shown in fig. 3.
Since the grid distance is c °, the coordinates of the longitude and latitude of the 4 vertices A, B, C, D of the grid can be given, where the coordinates of the longitude and latitude (lat 1, long 1) of the lower left corner B of the grid point are most easily determined, and the longitude and latitude (lat 1, long 1) is known as how many grids are spaced from the vertex of the lower left corner of the first grid point (also representing the longitude and latitude of the first grid point) in the east-west direction and how many grids are known in the north-south direction, i.e., the initial longitude and latitude of the lower left corner are respectively added by several grid distances in the spaced direction, and next, the coordinate position is A, C, D, and if displayed in the (latitude, longitude) format, the coordinate position is respectively (lat 1 ° + c °, long1 ° + c ° + lat1 ° + c °, long1 ° + c °).
And the calculation method of the 4 vertex coordinate positions of other grids is analogized in turn. As shown in fig. 3, if the equal longitude and latitude unit grid ABCD is used as the initial reference grid point, the longitude and latitude coordinates of the 4 vertexes A1, B1, C1, D1 enclosed by the first 4 dotted lines on the right side of the grid (ABCD grid in the figure) are also easily determined, wherein the positions of the vertex B1 on the lower left corner and the vertex A1 on the upper left corner are most easily determined, and the coordinates are the C point coordinate and the D point coordinate of the grid on the left side, respectively, (lat 1 ° + C °), and the coordinate positions of the vertex C1 on the lower right corner and the vertex D1 on the upper right corner are (lat 1 °, long1 ° +2C °), (lat 1 ° + C °, long1 ° +2 ° and + C °), and so on the left corner, the vertex on the lower right corner, the vertex on the left corner, the upper right corner, and the 4 vertexes on the right corner are (lat 1 ° +2 ° + C °), and (lat 1 ° +3 ° + C °).
As shown in fig. 3, it is easy to determine the longitude and latitude coordinates of the 4 vertices a ', B', C ', D' (fig. 3) of the grid surrounded by the first 4 dotted lines immediately above the A1B1C1D1 grid, wherein the positions of the lower left vertex B ', lower right vertex C' are most easily determined, and the coordinates are the A1 point coordinate and D1 point coordinate of the lower grid, respectively, (lat 1 ° + C °, long1 °), (lat 1 ° + C °, long1 ° + C °, and the coordinate positions of the upper left vertex and upper right vertex a ', D' are (lat 1 ° +2C °, long1 °), (lat 1 ° + C °, long1 ° + C °, and so on, the lower left vertex, upper left vertex, lower right vertex, upper right vertex, long 4 ° +2 ° coordinates of the second grid above vertically (lat 1 ° +2C °), (lat 1 ° +3 ° + long (long 1 ° + C °), and long1 ° + C °.
And determining the longitude and latitude coordinates of the four vertexes of other grids in the same way, and so on, and the description is omitted.
Since the resulting area of the mesoscale numerical mode output is not very large, no special cases are involved here, such as cross-equator, cross greenwich line, etc.
And 4, if the boundary of the final equal-longitude and latitude unit single grid to be converted exceeds the regional boundary of the original equal-longitude unit grid in the research region, the grid point is not converted and assigned, the conversion work of the corresponding equal-grid-distance grid is terminated, the conversion work of the rest grid points in the research region is further carried out, and the step 3 is repeatedly carried out. All areas covered by the final equal latitude and longitude unit grid acquired later will not contain the grid terminating the conversion.
And 5, if the boundary of the final equal-longitude-latitude unit single grid to be converted does not exceed the regional boundary of the original equal-longitude-latitude unit grid, performing the step.
And 5, acquiring the sub-areas of the converted equal-longitude and latitude unit single grid formed by the original equal-kilometer unit grids according to the four vertexes of the determined equal-longitude and latitude unit grid to be converted.
After determining four vertexes of a unit grid with equal longitude and latitude to be converted, the grid is determined, then the electric wave environment parameters of the converted grid are determined most importantly, and for the solution of each electric wave environment parameter, the method is obtained by weighting and calculating the electric wave environment parameters of the original several unit grids with equal kilometers according to the area sizes of the original several unit grids with equal kilometers forming the grid and the corresponding forming areas of the original several unit grids with equal kilometers.
After the four vertexes of the grid of the equal-longitude and latitude unit to be converted are determined, the new grid can be determined according to the longitude and latitude of the four vertexes, which are formed by the original equal-kilometer unit grids, because the longitude and latitude of the four vertexes of each original equal-kilometer unit grid can be read from the data file shown in fig. 2, according to the size and the range of the longitude and latitude, the new grid can be determined according to simple geometric relationship and inclusion relationship, the specific number is related to the grid distance of the grid and the converted grid distance, and the description is omitted.
For example, for the new grid of equal longitude and latitude units shown in (lat 1, long 1) in fig. 3, the vertex is A, B, C, D, the area of which is composed of partial areas of the original several equal kilometer unit grids, the specific number of the original several equal kilometer unit grids is related to the grid distance of the original grid and the converted grid distance, and the number of the original several equal kilometer unit grids composing the new grid can be determined according to the size and range of the longitude and latitude and according to the simple geometric relationship and the inclusion relationship.
For convenience of description, it is assumed that a partial area composition of 9 original grids of equal kilometers is determined as a schematic diagram (as shown in fig. 3), in which the original grids are (north latitude 2, longitude 1), (north latitude 3, longitude 2), (north latitude 4, longitude 2), (north latitude 2, longitude 3), (north latitude 3, longitude 3), (north latitude 4, longitude 3), (north latitude 2, longitude 4), (north latitude 3, longitude 4), and (north latitude 4, longitude 4), respectively, and the grids are identified according to a (latitude, longitude) format, wherein the latitude and the longitude are corresponding to the lower left vertex of the grids.
Sub-area S of original equal kilometer unit grid in final equal longitude and latitude unit single grid Score of origin The calculation is as follows:
S score of origin ＝d Original weft dividing ×d Original meridian point (1)
Wherein d is Original weft dividing Forming final equal longitude and latitude unit single grid area-divided latitude distance for original equal kilometer unit grid Original meridian point The longitude distance of the area of the original equal-kilometer unit grid occupying the final equal-longitude and latitude unit single grid is obtained. If the 9 original equal kilometer unit grids are divided into areas with latitude distances d Original weft dividing Distance d from longitude Original meridian point With (d) Original weft dividing ，d Original meridian branch ) Formally, latitude distance d of 1 to 9 parts area Original weft dividing Distance d from longitude Original meridian point Then it is:
(latitude 3-lat1, longitude 1-long 1), (latitude 4-latitude 3, longitude 3-long 1), (lat 1 ° + c ° -latitude 4, longitude 4-longitude 3), (latitude 4-latitude 3, longitude 4-longitude 3), (latitude 3-lat1, longitude 4-longitude 3), (lat 1 ° + c ° -latitude 4, long1+ c ° -longitude 4), (latitude 4-latitude 3, long1+ c ° -longitude 4), (latitude 3-lat1, long1+ c ° -longitude 4).
According to the formula (1), the following compounds can be obtained:
S score 1 = (latitude 3-lat 1) × (longitude 1-long 1)
S Score 2 = (latitude 4-latitude 3) × (longitude 3-long 1)
S Score 3 = (lat 1+ c-latitude 4) × (longitude 3-long 1)
S Score 4 = (lat 1+ c-latitude 4) × (longitude 4-longitude 3)
S Score 5 = (latitude 4-latitude 3) × (longitude 4-longitude 3)
S Score 6 = (latitude 3-lat 1) × (longitude 4-longitude 3)
S Score 7 = (lat 1+ c-latitude 4) × (long 1+ c-longitude 4)
S Score 8 = (latitude 4-latitude 3) × (long 1+ c-longitude 4)
S Score 9 = (latitude 3-lat 1) × (long 1+ c-longitude 4)
Step 6, according to the sub-area S of the single grid which is obtained and formed by the original grids of the same kilometer unit and has the same longitude and latitude unit Score of origin Obtaining the area ratio S of the original grids of the same kilometer unit to form the single grids of the same longitude and latitude unit Ratio of occupation of 。
Based on (2), the area ratio S of original single grids in the unit of equal kilometers and forming grids in the unit of equal longitude and latitude can be obtained Ratio of occupation of ：
S Ratio of occupation of ＝(S Score of origin /S Finally, the product is processed )×100％ (2)
Wherein S Score of origin Is the sub-area of the original equal kilometer unit grid occupying the final equal longitude and latitude unit single grid, S Finally, the product is processed The area of the single grid is the final unit of equal longitude and latitude.
Solving the area S of a single grid with equal longitude and latitude units based on the formula (3) Finally, the product is processed ：
S Finally, the product is processed ＝d Weft yarn ×d Warp beam (3)
Wherein d is Weft yarn The final latitude distance of a single grid with equal latitude and longitude units, d Warp beam The longitude distance of the single grid is the final equal longitude and latitude unit.
Then, the area ratio of the single grid which is formed by the original 9 grids of the equal kilometer units and forms the final equal longitude and latitude unit can be obtained according to the formula (2):
S ratio of 1 ＝(S Score 1 /S Finally, the product is processed )×100％
S Ratio of 2 ＝(S Score 2 /S Finally, the product is processed )×100％
S Ratio of 3 ＝(S Score 3 /S Finally, the product is processed )×100％
S Ratio of 4 ＝(S Score 4 /S Finally, the product is processed )×100％
S Ratio of 5 ＝(S Score 5 /S Finally, the product is processed )×100％
S Ratio of 6 ＝(S Score 6 /S Finally, the product is processed )×100％
S Ratio of 7 ＝(S Score 7 /S Finally, the product is processed )×100％
S Ratio of 8 ＝(S Score 8 /S Finally, the product is processed )×100％
S Ratio of 9 ＝(S Score 9 /S Finally, the product is processed )×100％
Step 7, forming the area ratio S of the single grid with the same longitude and latitude unit according to the obtained original grid with the same kilometer unit Ratio of occupation of Obtaining troposphere electric wave environment parameter W of original kilometer unit grids Original source Single grid ratio S of final equal longitude and latitude unit Ratio of occupation of Component W of Account for 。
And reading the troposphere electric wave environment parameter values corresponding to each original equal kilometer unit grid from the data file of the figure 2 according to the longitude and latitude. Acquiring troposphere electric wave environment parameters W of original kilometer unit grids according to formula (4) Original source Single grid ratio S of final equal longitude and latitude unit Ratio of occupation of Component W of Account for ：
W Account for ＝W Original source ×S Ratio of occupation of (4)
The troposphere electric wave environmental parameters 1 to N corresponding to the original kilometer unit grids can be read from the data file as required, so that convenience is brought to the readingHere, only one parameter is read, and each of the parameter values corresponding to 1 to 9 original kilometer unit grids is W Original 1 、W 2. Sup. St 、W Original 3 、W Original 4 、W Original 5 、W Original 6 、W Original 7 、W Original 8 、W Original 9 。
Then the components of the single grid ratio of the final equal longitude and latitude unit formed by the environmental parameters of certain tropospheric electric waves of 9 original equal kilometer unit grids are respectively:
W account for 1 ＝W Original 1 ×S Ratio of 1
W Account for 2 ＝W 2. Sup. St ×S Ratio of 2
W Account for 3 ＝W Original 3 ×S Ratio of 3
W Account for 4 ＝W Original 4 ×S Ratio of 4
W Account for 5 ＝W Original 5 ×S Ratio of 5
W Account for 6 ＝W Original 6 ×S Ratio of 6
W Account for 7 ＝W Original 7 ×S Ratio of 7
W Account for 8 ＝W Original 8 ×S Ratio of 8
W Account 9 ＝W 9 original ×S Ratio of 9
Step 8, according to the acquired environmental parameter value W of the troposphere electric wave of one of the original kilometer unit grids Original source Single grid ratio S of final equal longitude and latitude unit Ratio of occupation of Component W of Account for The sum of these components is the obtained final converted grid value W of one of the equal grid distance troposphere electric wave environment parameters set in the unit of grid distance longitude and latitude Finally, the product is processed 。
Obtaining a final grid value W of the electric wave environment parameter of the troposphere of a certain equal grid distance by the formula (5) Finally, the product is processed ：
W Finally, the product is processed ＝W Account for 1 +W Account for 2 +W Account 3 +W Account for 4 +W Account 5 +W Account for 6 +W Account for 7 +W Account for 8 +W Account 9 (5)
If other tropospheric electric wave environment parameters of the final equal-grid-distance grid are to be acquired, the steps 7-8 can be repeated, and the tropospheric electric wave environment parameters 1 to N are acquired according to the requirement. Repeating the steps 3-8, and obtaining the environmental parameters 1 to N of the troposphere electric wave of all final equal grid distance grids in the research area.
In embodiment 1, the 4 vertexes A, B, C, D of the new mesh in the final equal longitude and latitude unit enclosed by the 4 dotted lines in the figure (as shown in fig. 4) are determined, and two electric wave environment parameters of the air temperature and the atmospheric boundary layer height of the new mesh point are obtained.
Step 1, in a numerical mode, a user sets a horizontal grid distance d km in advance, wherein the horizontal grid distance d km is usually a fixed value as shown in fig. 1, and d is less than 60km and takes a value of 30km. The value is too large because of the non-uniform atmospheric level, and the converted data is easy to cause large errors;
step 2, the user obtains the 30km (fixed value preset by user in mode) output by the mesoscale atmospheric numerical mode as the equal kilometer unit grid data with equal space, the grid data file contains unequal space latitude and longitude information corresponding to each grid data and troposphere electric wave environment parameter information to be researched, the data sequence is arranged as figure 5, 1-9 columns of data in the data table are respectively latitude, longitude and altitude corresponding to the grid, each longitude number of the grid, grid sequence number in north and south direction under each longitude of the grid, air pressure value, air temperature value, water-vapor mixing ratio and air boundary layer height value, each grid in the data table does not calculate height information, 4 electric wave environment parameters are respectively air pressure value P and air temperature value t c Water-vapor mixing ratio q v And an atmospheric boundary layer height value PBLH. The first two columns 11.04795, 102.9431 in the first row represent the longitude and latitude corresponding to the first grid (row 1, column 1) of the start of the numerical mode output result, and so on, the first two columns 11.322, 102.9431 in the second row represent the second grid in the north-south direction, and so on. It should be noted that the grid data file may contain many electric wave environment parameters, and the number of parameters and data format are not limited, as long as the longitude and latitude and electric wave ring of the corresponding grid point can be obtainedThe environmental parameters are only needed.
And 3, determining four vertexes of a certain final equal-longitude and latitude unit grid to be converted according to the specific equal-kilometer unit grid distance (grid distance 30 km) of the troposphere electric wave environment numerical mode result and the specific equal-longitude and latitude unit grid distance (grid distance 0.5 degree) to be converted.
In the grid chart shown in FIG. 4, latitudes 1 to 7 are 11.04795 °, 11.322 °, 11.59579 °, 11.86932 °, 12.14256 °, 12.41553 °, 12.68821 °, and longitudes 1 to 7 are 102.9431 °,103.2225 °,103.5018 °,103.7812 °, 104.0606 °, respectively. The grid corresponding to latitude 1 and longitude 1 is the 1 st grid, as denoted by W in the figure. If the grid conversion from the grid of the same kilometer unit to the grid of the same longitude and latitude unit is performed in the grid in the figure, for example, the grid to be converted into the grid of the longitude and latitude (11.5 degrees and 103 degrees), and the grid distance is 0.5 degree, 4 vertexes A, B, C, D of the grid need to be obtained, wherein the position of the vertex B is most easily determined, the coordinates thereof are the north latitude 11.5 degrees, the east longitude 103 degrees, the A, C, D positions, and the coordinate positions thereof are respectively the north latitude 11.5 degrees +0.5 degrees, the east longitude 103 degrees +0.5 degrees, and the analogy is sequentially performed for the calculation methods of the 4 vertex coordinate positions of other grids. Since the resulting area of the mesoscale numerical mode output is not very large, special cases such as cross-equator, cross-greenwich mean line, etc. are not involved here.
And 4, if the boundary of the final equal-longitude and latitude unit single grid to be converted exceeds the regional boundary of the original equal-kilometer unit grid, the grid point is not converted and assigned, the conversion work of the corresponding equal-grid-distance grid is terminated, the conversion work of the rest grid points is further carried out, and the step 3 is repeatedly carried out. All areas covered by the final equal latitude and longitude unit grid acquired later will not contain the grid terminating the conversion.
And 5, if the boundary of the final equal-longitude-latitude unit single grid to be converted does not exceed the regional boundary of the original equal-longitude-latitude unit grid, performing the step.
And 5, acquiring the sub-areas of the converted final equal-longitude and latitude unit single grid formed by the original equal-kilometer unit grids according to the four vertexes of the determined final equal-longitude and latitude unit grid to be converted.
After determining four vertexes of a final equal-longitude-latitude unit grid to be converted, the grid is determined, then, most importantly, four electric wave environment parameters of the converted grid are determined, and for the solution of each electric wave environment parameter, the embodiment is obtained by weighting and calculating the electric wave environment parameters of the original several equal-longitude-latitude unit grids according to the area sizes of the original several equal-longitude-latitude unit grids forming the grid and the corresponding forming areas of the original several equal-longitude-latitude unit grids. For example, for the new mesh (11.5 °,103 °) shown in fig. 4, the vertices are A, B, C, D, and the coordinates of these vertices are obtained from the above steps, the area of the new mesh is composed of partial areas of the original several kilometer units of mesh, and it is known that it can be determined how many original kilometer units of mesh the new mesh is composed according to the size and range of the latitudes and longitudes, where a total of 6 partial areas of the original kilometer units of mesh are composed, and the original meshes are respectively (11.322 °,102.9431, (11.59579 °,102.9431 °), (11.86932 °,102.9431 °) (82 zxft 3434 °) (11.59579 °, 3638 °) (3724 zxft 374924 °) (374924 ° of mesh), where the vertices and longitudes of mesh are identified according to the left grid of the latitude and longitude is identified as 324924 ° of the mesh.
According to the formula (1), the sub-area S of the original equal kilometer unit grid occupying the final equal longitude and latitude unit single grid can be obtained Score of origin . In this embodiment, 6 original equal kilometer unit grids occupy the latitude distance d of the sub-area of the final equal longitude and latitude unit single grid Original weft dividing And a longitudinal distance d Original meridian branch Subtracting the nearest two adjacent latitude and longitude values, if the 6 latitude distances d of the division area of the original equal kilometer unit grid Original weft dividing Distance d from longitude Original meridian point With (d) Original weft division ，d Original meridian point ) Formally expressed, latitude distance d of 1 to 6 areas Original weft dividing Distance d from longitude Original meridian point Then it is: (11.59579 ° -11.5 °,103.2225 ° -103 °), (11.86932 ° -11.59579 °,103.2225 ° -103 °), (12.0 ° -11.86932 °,103.2225 ° -103 ° 3°)、(12.0°-11.86932°，103.5°-103.2225°)、(11.86932°-11.59579°，103.5°-103.2225°)、(11.59579°-11.5°，103.5°-103.2225°)。
According to the formula (1):
S score 1 ＝(11.59579-11.5)×(103.2225-103)＝0.0213
S Score 2 ＝(11.86932-11.59579)×(103.2225-103)＝0.0609
By analogy in the following way,
S score 3 ＝0.0291
S Score 4 ＝0.0363
S Original score 5 ＝0.0760
S Score 6 ＝0.0266
Step 6, according to the sub-area S of the single grid which is obtained and formed by the original grids of the same kilometer unit and has the same longitude and latitude unit Score of origin Obtaining the area ratio S of the single grid which is formed by the original grids of the same kilometer unit and forms the final unit of the same longitude and latitude Ratio of occupation of 。
Calculating the area S of the single grid with the final equal longitude and latitude unit based on the formula (3) Finally, the product is processed ：
S Finally, the product is processed ＝d Weft yarn ×d Warp beam ＝0.5×0.5＝0.25
Wherein d is Weft yarn The distance in latitude of the single grid in the final equal latitude and longitude unit is 0.5 DEG, d Warp beam The longitude distance of the single grid is the final equal latitude and longitude unit, here also 0.5 °.
Then, according to the formula (2), the area ratio of the single grid which is formed by the original 6 grids of the equal kilometers unit to the final equal longitude and latitude unit is obtained:
S ratio of 1 ＝(S Score 1 /S Finally, the product is processed )×100％＝8.52％
S Ratio of 2 ＝(S Score 2 /S Finally, the product is processed )×100％＝24.36％
S Ratio of 3 ＝(S Score 3 /S Finally, the product is processed )×100％＝11.64％
S Ratio of 4 ＝(S Score 4 /S Finally, the product is processed )×100％＝14.52％
S Ratio of 5 ＝(S Score 5 /S Finally, the product is processed )×100％＝30.4％
S Ratio of 6 ＝(S Score 6 /S Finally, the product is processed )×100％＝10.64％
Where S is Score 1 、S Score 2 、S Score 3 、S Score 4 、S Score 5 、S Score 6 、S Score 7 、S Score 8 、S Score 9 The area of the original kilometer unit grids from 1 to 6 in the same direction accounts for the final longitude and latitude unit single grid.
Step 7, forming the area ratio S of the single grid with the same longitude and latitude unit according to the obtained original grid with the same kilometer unit Ratio of occupation of Acquiring tropospheric electric wave environment parameters W of the original kilometer unit grid according to the formula (4) Original source Single grid ratio S of final equal longitude and latitude unit Ratio of occupation of Component W of Account for 。
The troposphere electric wave environmental parameters 1 to 4 corresponding to the original kilometer unit grids can be read from the data file of FIG. 5 as required, and the gas temperature value t is respectively read at the positions c The air temperature values corresponding to 1-6 original kilometer unit grids are respectively 28.8, 28.2, 26.5, 28.9, 29.1 and 29.3.
Acquiring troposphere electric wave environment parameters W of original kilometer unit grids according to formula (4) Original source Single grid ratio S of final equal longitude and latitude unit Ratio of the ingredients Component W of Account for Here, the air temperature value t c Component W of Account for ：
W Account for 1 ＝W Original 1 ×S Ratio of 1 ＝28.8×8.52％＝2.45°
W Account for 2 ＝W Original 2 ×S Ratio of 2 ＝28.2×24.36％＝6.87°
W Account 3 ＝W Original 3 ×S Ratio of 3 ＝26.5×11.64％＝3.08°
W Account for 4 ＝W Original 4 ×S Ratio of 4 ＝28.9×14.52％＝4.20°
W Account for 5 ＝W Original 5 ×S Ratio of 5 ＝29.1×30.4％＝8.85°
W Account for 6 ＝W Original 6 ×S Ratio of 6 ＝29.3×10.64％＝3.12°
Step 8, according to the acquired troposphere electric wave environmental parameter gas temperature value W of the original kilometer unit grids Original source Single grid ratio S of final equal longitude and latitude unit Ratio of occupation of Component W of Account for Adding the obtained components, i.e. the obtained final converted tropospheric electric wave environmental parameter temperature t of the final equal grid of latitude and longitude units set in equal grid distance latitude and longitude units c Grid value W of Finally, the product is processed 。
Obtaining the final grid value W of the troposphere electric wave environmental parameter of the final equal longitude and latitude unit grid distance by the formula (5) Finally, the product is processed ：
W Finally, the product is processed ＝W Account for 1 +W Account for 2 +W Account for 3 +W Account for 4 +W Account for 5 +W Account for 6 ＝28.6°
And if other tropospheric electric wave environment parameters of the final equal-grid-distance grid are to be acquired, repeating the steps 7-8, and acquiring different tropospheric electric wave environment parameters 1-N according to requirements.
Repeating the step 7, reading any parameter of troposphere electric wave environment parameters 1 to 4 corresponding to the original kilometer unit grids from the data file of fig. 5 as required, and reading the atmospheric boundary layer height value PBLH parameter at the position, wherein the atmospheric boundary layer height values corresponding to 1-6 original kilometer unit grids are respectively 98.9, 148.6, 456.6, 431.2, 165.2 and 99.8.
Obtaining troposphere electric wave environment parameter W of original equal kilometer unit grids according to formula (4) Original source Single grid ratio S of final equal longitude and latitude unit Ratio of occupation of Component W of Account for 。
Component W of atmospheric boundary layer height value PBLH Account for ：
W Account for 1 ＝W Original 1 ×S Ratio of 1 ＝98.9×8.52％＝8.43
W Account for 2 ＝W 2. Sup. St ×S Ratio of 2 ＝148.6×24.36％＝36.20
W Account for 3 ＝W Original 3 ×S Ratio of 3 ＝456.6×11.64％＝53.15
W Account for 4 ＝W Original 4 ×S Ratio of 4 ＝431.2×14.52％＝62.61
W Account 5 ＝W Original 5 ×S Ratio of 5 ＝165.2×30.4％＝50.22
W Account for 6 ＝W Original 6 ×S Ratio of 6 ＝99.8×10.64％＝10.62
Repeating the step 8, and acquiring the grid value W of the converted height value PBLH of the new grid atmospheric boundary layer according to the formula (5) Finally, the product is processed ：
W Finally, the product is processed ＝W Account for 1 +W Account for 2 +W Account for 3 +W Account for 4 +W Account 5 +W Account for 6 ＝221.2(m)
The tropospheric electric wave environment parameters 1 to N of the final equal-grid-distance grid can be obtained as required.
Repeating the steps 3-8, and obtaining the environmental parameters 1 to N of the troposphere electric wave of all final equal grid distance grids in the research area.
In the embodiment 1, the non-equidistant longitude and latitude unit grid data is converted into the final equidistant longitude and latitude unit grid data, and the new grid values corresponding to two electric wave environment parameters of the air temperature and the atmospheric boundary layer height of the new grid point are obtained, so that the method is applicable and practical, and the converted final equidistant longitude and latitude unit grid data is convenient for the practical application of researching the electric wave environment characteristics of the troposphere in a large area, particularly the global range.
Embodiment 2, the first 4 dotted lines next to the horizontal right side of the ABCD grid in the figure enclose 4 vertexes A1, B1, C1, D1 of a new grid with equal latitude and longitude units (as shown in fig. 4), and the air temperature electric wave environment parameter of the new grid point is obtained.
Step 1, in a numerical mode, a user sets a horizontal grid distance d km in advance, wherein the horizontal grid distance d km is usually a fixed value as shown in fig. 1, and d is less than 60km and takes a value of 30km. The value is too large because of the non-uniform atmospheric level, and the converted data is easy to cause large errors;
step 2, the user obtains the 30km (fixed value preset by the user in the mode) output by the mesoscale atmospheric numerical mode as the equi-kilometer unit grid data with equal intervals, the grid data file contains the longitude and latitude information with unequal intervals corresponding to each grid data in the research area and the troposphere electric wave environmental parameter information to be researched, the data sequence is arranged as shown in figure 5, 1-9 columns of data in the data table are respectively the latitude, longitude and altitude corresponding to the grid, each longitude number of the grid, the south-north grid column number under each longitude of the grid, the air pressure value, the air temperature value, the water-vapor mixing ratio and the air boundary layer height value, each grid in the data table does not calculate the height information, 4 electric wave environmental parameters are respectively the air pressure value P and the air temperature value t c Water-vapor mixing ratio q v And an atmospheric boundary layer height value PBLH. The first two columns 11.04795, 102.9431 in the first row represent the longitude and latitude corresponding to the first grid (row 1, column 1) of the start of the numerical mode output result, and so on, the first two columns 11.322, 102.9431 in the second row represent the second grid in the north-south direction, and so on. It should be noted that each grid in the grid data file may include many parameters of the radio wave environment, and the number of the parameters and the data format are not limited as long as the longitude and latitude and the radio wave environment parameters of the corresponding grid point can be obtained.
And 3, determining four vertexes of a certain final equal longitude and latitude unit grid to be converted according to the specific equal kilometer unit grid distance (the grid distance of the example is 30 km) of the troposphere electric wave environment numerical value mode result and the specific equal longitude and latitude unit grid distance to be converted (the grid distance of the example is 0.5 degrees).
In the grid chart shown in FIG. 4, latitudes 1 to 7 are 11.04795 °, 11.322 °, 11.59579 °, 11.86932 °, 12.14256 °, 12.41553 °, 12.68821 °, and longitudes 1 to 7 are 102.9431 °,103.2225 °,103.5018 °,103.7812 °, 104.0606 °, respectively. The grid corresponding to latitude 1 and longitude 1 is the 1 st grid, as denoted by W in the figure. If the grid in the figure is converted into a grid in equal kilometer units to equal latitude and longitude units, for example, if the grid to be converted into a grid in latitude and longitude (11.5 degrees and 103.5 degrees) has a grid distance of 0.5 degrees, 4 vertexes A1, B1, C1 and D1 of the grid need to be obtained, wherein the position of the vertex B1 is most easily determined, the coordinates of the vertex B1 are 11.5 degrees of north latitude and 103.5 degrees of east longitude representing the grid coordinates, the positions of the vertex A1, C1 and D1 are respectively 11.5 degrees of north latitude +0.5 degrees of east longitude 103.5 degrees, 11.5 degrees of north latitude, 103.5 degrees of east longitude +0.5 degrees of north longitude, 103.5 degrees of east longitude +0.5 degrees of east longitude, and the analogy is sequentially performed on the 4 vertex coordinate positions of other grids. Since the resulting area of the mesoscale numerical mode output is not very large, special cases such as cross-equator, cross-greenwich mean line, etc. are not involved here.
And 4, if the boundary of the final equal-longitude and latitude unit single grid to be converted exceeds the regional boundary of the original equal-kilometer unit grid, the grid point is not converted and assigned, the conversion work of the corresponding equal-grid-distance grid is terminated, the conversion work of the rest grid points is further carried out, and the step 3 is repeatedly carried out. All areas covered by the final equal latitude and longitude unit grid acquired later will not contain the grid terminating the conversion.
And 5, if the boundary of the final equal-longitude-latitude unit single grid to be converted does not exceed the regional boundary of the original equal-longitude-latitude unit grid, performing the step.
And 5, acquiring the sub-areas of the converted final equal-longitude and latitude unit single grid formed by the original equal-kilometer unit grids according to the four vertexes of the determined final equal-longitude and latitude unit grid to be converted.
After determining four vertexes of a final equal longitude and latitude unit grid to be converted, the grid is determined, then most importantly, the electric wave environment parameters of the converted grid are determined, and for the solution of each electric wave environment parameter, the embodiment calculates the electric wave environment parameters of the original equal kilometer unit grids according to the area size of the original equal kilometer unit grids forming the grid and the corresponding formed area weighting. For example, for the mesh (11.5 °,103.5 °) shown in fig. 4, the vertices are A1, B1, C1, and D1, the area of the vertices is composed of partial areas of the original several kilometer-equivalent unit meshes, the specific number of the original several kilometer-equivalent unit meshes is related to the mesh distance of the original mesh and the mesh distance after conversion, and according to the size and range of the longitude and latitude, it can be determined how many original kilometer-equivalent unit meshes the new mesh is composed according to a simple geometric relationship. A total of 9 partial areas of the original kilometer units of the grid are formed, the original grid being (11.322 °,103.2225 °), (11.59579 °,103.2225 °), (11.86932 °,103.2225 °), (11.86932 °,103.5018 °), (11.59579 °,103.5018 °), (11.322 °,103.5018 °), (11.86932 °,103.7812 °), (11.59579 °,103.7812 °), (11.322 °,103.7812 °), respectively, the grid being identified in a latitude (longitude) format, wherein the latitude and longitude are the left-lower-left vertices of the corresponding grid.
According to the formula (1), the sub-area S of the original equal kilometer unit grid in the final equal longitude and latitude unit single grid can be obtained Score of origin . In this example, 9 original kilometer unit grids occupy the latitude distance d of the sub-area of the final equal longitude and latitude unit single grid Original weft division And a longitudinal distance d Original meridian point Subtracting the nearest two adjacent latitude and longitude values, if the 9 latitude distances d of the original equal kilometer unit grid sub-area Original weft dividing Distance d from longitude Original meridian branch With (d) Original weft dividing ，d Original meridian point ) Formally, latitude distance d of 1 to 9 parts area Original weft dividing Distance d from longitude Original meridian point Then it is:
(11.59579°—11.5°，103.5018°—103.5°)、(11.86932°—11.59579°，103.5018°—103.5°)、(12.0°—11.86932°，103.5018°—103.5°)、(12.0°—11.86932°，103.7812°—103.5018°)、(11.86932°—11.59579°，103.7812°—103.5018°)、(11.59579°—11.5°，103.7812°—103.5018°)、(12.0°—11.86932°，104°—103.7812°)、(11.86932°—11.59579°，104°—103.7812°)、(11.59579°—11.5°，104°—103.7812°)。
according to formula (1):
S score 1 ＝(11.59579-11.5)×(103.5018-103.5)＝1.72422×10 -4
S Score 2 ＝(11.86932-11.59579)×(103.5018-103.5)＝4.9235×10 -4
By analogy in the following way,
S score 3 ＝2.3522×10 -4
S Score 4 ＝0.0365
S Score 5 ＝0.0764
S Score 6 ＝0.0268
S Score 7 ＝0.0286
S Score 8 ＝0.0598
S Score 9 ＝0.0210
Step 6, according to the sub-area S of the single grid which is obtained and formed by the original grids of the same kilometer unit and has the same longitude and latitude unit Score of origin Obtaining the area ratio S of the single grid which is formed by the original grids of the same kilometer unit and forms the final unit of the same longitude and latitude Ratio of occupation of 。
Calculating the area S of the single grid with the final equal longitude and latitude unit based on the formula (3) Finally, the product is processed ：
S Finally, the product is processed ＝d Weft yarn ×d Warp beam ＝0.5×0.5＝0.25
Wherein d is Weft yarn The distance in latitude of the single grid in the final equal latitude and longitude unit is 0.5 DEG, d Warp beam The longitude distance of the single grid is the final equal latitude and longitude unit, here also 0.5 °.
Then, the area ratio of the single grid which is formed by the original 9 grids of the equal kilometer units and forms the final equal longitude and latitude unit can be obtained according to the formula (2):
S ratio of 1 ＝(S Score 1 /S Finally, the product is processed )×100％＝0.0687％
S Ratio of 2 ＝(S Score 2 /S Finally, the product is processed )×100％＝0.1970％
S Ratio of 3 ＝(S Score 3 /S Finally, the product is processed )×100％＝0.0941％
S Ratio of 4 ＝(S Score 4 /S Finally, the product is processed )×100％＝14.6％
S Ratio of 5 ＝(S Score 5 /S Finally, the product is processed )×100％＝30.56％
S Ratio of 6 ＝(S Score 6 /S Finally, the product is processed )×100％＝10.72％
S Ratio of 7 ＝(S Score 7 /S Finally, the product is processed )×100％＝11.44％
S Ratio of 8 ＝(S Score 8 /S Finally, the product is processed )×100％＝23.92％
S Ratio of 9 ＝(S Score 9 /S Finally, the product is processed )×100％＝8.4％
Where S is Score 1 、S Score 2 、S Score 3 、S Score 4 、S Score 5 、S Score 6 、S Score 7 、S Score 8 、S Score 9 The area of the original kilometer unit grid from 1 to 9 in the same direction accounts for the final longitude and latitude unit single grid.
Step 7, forming the area ratio S of the single grid with the same longitude and latitude unit according to the obtained original grid with the same kilometer unit Ratio of the ingredients Acquiring tropospheric electric wave environment parameters W of the original kilometer unit grid according to the formula (4) Original source Single grid ratio S of final equal longitude and latitude unit Ratio of occupation of Component W of Account for 。
Reading tropospheric electric wave environment parameters 1 to 4 corresponding to the original kilometer unit grids from the data file as required, and reading a parameter gas temperature value t at the position c The corresponding air temperature values of 1-9 original kilometer unit grids are 29.3, 29.1, 28.9, 29.0, 29.2, 29.4, 28.7, 29.2 and 29.5 respectively, and the tropospheric electric wave environment parameter W of the original kilometer unit grids is obtained according to the formula (4) Original source Single grid ratio S of final equal longitude and latitude unit Ratio of occupation of Component W of Account for ：
W Account 1 ＝W Original 1 ×S Ratio of 1 ＝29.3×0.0687％＝0.020°
W Account for 2 ＝W 2. Sup. St ×S Ratio of 2 ＝29.1×0.1970％＝0.057°
W Account 3 ＝W Original 3 ×S Ratio of 3 ＝28.9×0.0941％＝0.027°
W Account for 4 ＝W Original 4 ×S Ratio of 4 ＝29.0×14.6％＝4.234°
W Account for 5 ＝W Original 5 ×S Ratio of 5 ＝29.2×30.56％＝8.924°
W Account for 6 ＝W Original 6 ×S Ratio of 6 ＝29.4×10.72％＝3.152°
W Account for 7 ＝W Original 7 ×S Ratio of 7 ＝28.7×11.44％＝3.283°
W Account for 8 ＝W Original 8 ×S Ratio of 8 ＝29.2×23.92％＝6.985°
W Account for 9 ＝W 9 original ×S Ratio of 9 ＝29.5×8.4％＝2.478°
Step 8, according to the acquired troposphere electric wave environmental parameter gas temperature value W of the original kilometer unit grid Original source Single grid ratio S of final equal longitude and latitude unit Ratio of the ingredients Component W of Account for Adding the obtained components to obtain a final converted tropospheric electric wave environmental parameter gas temperature grid value W of the final equal longitude and latitude unit grid set in equal grid distance longitude and latitude units Finally, the product is processed 。
Obtaining the grid value W of the troposphere electric wave environment parameter of the final equal longitude and latitude unit grid distance by the formula (5) Finally, the product is processed ：
W Finally, the product is processed ＝W Account for 1 +W Account for 2 +W Account for 3 +W Account for 4 +W Account for 5 +W Account for 6 +W Account for 7 +W Account for 8 +W Account for 9 ＝29.2°
And if other tropospheric electric wave environment parameters of the final equal-longitude-latitude unit grid are to be acquired, repeating the steps 7-8, and acquiring different tropospheric electric wave environment parameters from 1 to N according to requirements. Repeating the steps 3-8, and obtaining the environmental parameters 1 to N of the troposphere electric waves of the unit grids with equal longitude and latitude in all the research areas.
Embodiment 2 converts the grid data of the non-equidistant longitude and latitude units into new grid data of equidistant longitude and latitude units (the new grid of the converted equidistant longitude and latitude units is positioned on the horizontal right side of the new grid of the equidistant longitude and latitude units in embodiment 1), and obtains a new grid value corresponding to the electric wave environment parameter air temperature of the new grid point, which shows that the method of the invention is usable and practical, and the converted mesh data of the equidistant longitude and latitude units is convenient for the practical application of the large-area, especially the global range research on the troposphere electric wave environment characteristics.

Claims (2)

1. A method for converting different unit grid distances of tropospheric radiowave environment numerical patterns, comprising the steps of:
step 1, in a numerical mode, a user sets a research area and a horizontal grid distance d km in advance, wherein d is less than 60;
step 2, obtaining an equal kilometer unit grid data file which is output by a user in a mesoscale atmospheric numerical mode and takes d km as equal space, wherein the grid data file comprises longitude and latitude information of unequal space corresponding to each grid data in a research area and troposphere electric wave environmental parameter information to be researched;
step 3, determining four vertexes of a final equal-longitude and latitude unit grid to be converted in the research area according to a specific equal-kilometer unit grid distance d km of a troposphere electric wave environment numerical mode result and a grid distance c of a specific equal-longitude and latitude unit to be converted;
step 4, if the boundary of the final equal longitude and latitude unit grid to be converted exceeds the regional boundary of the original equal longitude and latitude unit grid in the research region, the grid point is not converted and assigned, the conversion work of the corresponding equal grid distance grid is terminated, the conversion work of the rest grid points in the research region is further carried out, the step 3 is repeatedly carried out, and all regions covered by the final equal longitude and latitude unit grid obtained later do not contain the grid which is terminated;
if the boundary of the final equal-longitude-latitude unit single grid to be converted does not exceed the regional boundary of the original equal-longitude-latitude unit grid, performing step 5;
step 5, acquiring the sub-area of the converted final equal longitude and latitude unit single grid formed by the original equal kilometer unit grids according to the four vertexes of the determined equal longitude and latitude unit grid to be converted;
single network with original equal kilometer unit grid occupying final equal longitude and latitude unitFractional area S of the grid Score of origin The calculation is as follows:
S score of origin ＝d Original weft division ×d Original meridian point (1)
Wherein, d Original weft dividing The latitude distance, d, of the final equal longitude and latitude unit single grid of the original equal kilometer unit grid Original meridian branch The longitude distance of the area of the original equal-kilometer unit grid occupying the final equal-longitude and latitude unit single grid is used;
step 6, according to the obtained sub-area S Score of origin Obtaining the area ratio S of the original equal kilometer unit grids to form the final equal longitude and latitude unit grids Ratio of occupation of ；
The area ratio S of the original equal-kilometer unit single grid forming the final equal-longitude and latitude unit grid is obtained based on the following formula (2) Ratio of the ingredients ：
S Ratio of occupation of ＝(S Score of origin /S Finally, the product is processed )×100％ (2)
Wherein S is Finally, the product is processed The area of the single grid is the final unit with equal longitude and latitude;
obtaining S based on the following formula (3) Finally, the product is processed ：
S Finally, the product is processed ＝d Weft yarn ×d Warp beam (3)
Wherein d is Weft yarn The final latitude distance of a single grid with equal latitude and longitude units, d Warp beam The longitude distance of a single grid is the final equal longitude and latitude unit;
step 7, acquiring troposphere electric wave environment parameters W of original equal kilometer unit grids Original source Single grid ratio S of final equal longitude and latitude unit Ratio of the ingredients Component W of Account for ；
Reading the troposphere electric wave environment parameter value corresponding to each original kilometer unit grid according to the longitude and latitude from the grid data file as required, and acquiring the troposphere electric wave environment parameter W of the original kilometer unit grid according to the following formula (4) Original source Single grid ratio S of final equal longitude and latitude unit Ratio of the ingredients Component W of Account for ：
W Account for ＝W Original source ×S Ratio of occupation of (4)
Step 8, the component W obtained in the step 7 is processed Account for Adding the obtained final converted grid value W of one of the equal grid distance troposphere electric wave environmental parameters set in the unit of grid distance longitude and latitude Finally, the product is processed 。
2. The method of converting different unit cell pitches of tropospheric electric wave environment numerical patterns as set forth in claim 1, wherein: setting the number of grids in a research area set by a user to be N multiplied by m, wherein each grid has N electric wave environment parameters; repeating the steps 7-8 to obtain other troposphere electric wave environment parameters of the final equal grid distance grid; and repeating the steps 3-9 to obtain the troposphere electric wave environmental parameters of other final equal grid distance grids in the research area.
CN202111674069.9A 2021-12-31 2021-12-31 Different unit grid distance conversion method for troposphere electric wave environment numerical mode result Active CN114417579B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111674069.9A CN114417579B (en) 2021-12-31 2021-12-31 Different unit grid distance conversion method for troposphere electric wave environment numerical mode result

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111674069.9A CN114417579B (en) 2021-12-31 2021-12-31 Different unit grid distance conversion method for troposphere electric wave environment numerical mode result

Publications (2)

Publication Number Publication Date
CN114417579A CN114417579A (en) 2022-04-29
CN114417579B true CN114417579B (en) 2023-01-24

Family

ID=81271480

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111674069.9A Active CN114417579B (en) 2021-12-31 2021-12-31 Different unit grid distance conversion method for troposphere electric wave environment numerical mode result

Country Status (1)

Country Link
CN (1) CN114417579B (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107945242A (en) * 2017-11-16 2018-04-20 中国科学院海洋研究所 It is a kind of towards IDL projection transform algorithms
CN109342658A (en) * 2018-10-30 2019-02-15 河南省环境保护科学研究院 The method that computer programing based on WRF model calculates atmospheric environment capacity

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107945242A (en) * 2017-11-16 2018-04-20 中国科学院海洋研究所 It is a kind of towards IDL projection transform algorithms
CN109342658A (en) * 2018-10-30 2019-02-15 河南省环境保护科学研究院 The method that computer programing based on WRF model calculates atmospheric environment capacity

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
GPS网与天文大地网坐标系统的转换方法;王刚等;《测绘学院学报》;20011230(第04期);全文 *

Also Published As

Publication number Publication date
CN114417579A (en) 2022-04-29

Similar Documents

Publication Publication Date Title
CN105095589B (en) A kind of mountain area power grid wind area is distributed drawing drawing method
CN111652975A (en) Method and system for evaluating available solar energy resources of urban building group
CN111292214A (en) Method for calculating high-precision Chinese offshore wave characteristic distribution
CN111652126B (en) Inversion radiation method based on satellite cloud image
CN112100922A (en) Wind resource prediction method based on WRF and CNN convolutional neural network
Yu et al. Global data assimilation experiments with scatterometer winds from SEASAT-A
McGregor et al. Nested simulations of perpetual January climate over the Australian region
CN114417579B (en) Different unit grid distance conversion method for troposphere electric wave environment numerical mode result
Yang et al. Regional spatiotemporal statistical database of evaporation ducts over the South China Sea for future long-range radio application
Šúri et al. Comparison of direct normal irradiation maps for Europe
CN109241212A (en) Based on mesoscale numerical value atmospheric model and high-resolution history rainfall inversion method
CN111582547B (en) Method for acquiring wind field distribution in different places by using wind field data set
CN112541620B (en) Typhoon storm water increasing prediction method and system with high prediction precision and efficiency
CN115841266A (en) Photovoltaic power generation potential evaluation method for photovoltaic power station site selection
Wang et al. Preliminary results of a new global ocean reanalysis
CN106446389A (en) Rapid data and image cutting method
CN112636893B (en) Method for improving eLoran system time service precision by using ASF grid and differential station
CN111582546A (en) Method and system for acquiring global rainfall information at different places by utilizing rainfall data set
CN117891962B (en) Graph database construction method and application of urban distributed photovoltaic system data
CN115238514B (en) Method and system for calculating satellite load observation simulation data
Zhou et al. The SIA method for spatial analysis of precipitation in the upper-middle reaches of the Yangtze River
Baghani Assessment of Rooftop Solar Power potential in Rural Areas using UAV Photogrammetry and GIS
Torres et al. Simulation of storm surge in northeast coast of the US; a closer look at the wind forcing
Dishon Determination of average ocean depths from bathymetric data
Sánchez et al. A study of the estimation of the photovoltaic potential at the urban level in tropical complex terrain

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