CN114417579A - 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
CN114417579A
CN114417579A CN202111674069.9A CN202111674069A CN114417579A CN 114417579 A CN114417579 A CN 114417579A CN 202111674069 A CN202111674069 A CN 202111674069A CN 114417579 A CN114417579 A CN 114417579A
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.)
Granted
Application number
CN202111674069.9A
Other languages
Chinese (zh)
Other versions
CN114417579B (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 depend on the acquisition and statistics of troposphere electric wave environmental parameters to a great extent, the acquisition and statistics of the parameters in the past can only depend on sparse and limited meteorological ground and exploration observation stations in an area for data acquisition and statistics, and for a vast ocean area, meteorological exploration stations are not available or rare 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 cooperation of the national atmospheric research center (NCAR) and the Pennsylvania State university, 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 cooperation of 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 (AFWA) and many organizations such as associations and universities. 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 environmental parameter space-time grid data based on the mesoscale atmospheric numerical mode, the user is required to perform study area pre-setting (study area size is user-defined as required) and horizontal space grid distance setting, the grid distance setting is set in units of equal kilometers (fixed value, assumed here as d km, pre-set by the user as required), the unit of the horizontal spatial grid distance of the acquired tropospheric electric wave environment parameter is km, in other words, the regional tropospheric electric wave environment data acquired based on the mesoscale atmospheric numerical mode are all composed of a plurality of grid data formed at equal intervals of d km (a fixed value set in advance by the user in the mode), if the grid data is converted into the latitude and longitude units, the grid data are not equidistant for the latitude and longitude units, and the grid data are greatly related to the projection mode of the numerical mode, the difference of north and south latitudes and the like.
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 characteristics of tropospheric radio wave environment 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 radio communication bureau (ITU-R) is reanalysis mesh data of equal-spaced latitude and longitude units (the data is the direct output result of the global atmospheric numerical mode), while the present example relates to the atmospheric mesoscale numerical mode, which considers the projection and coordinate modes, gives consideration to the representativeness of local climate and describes the fine structure of 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 characteristics of global or large-area radio wave environment needs to be converted into mesh data of equal-spaced latitude and longitude units, however, the specific conversion method of the different unit mesh data of the pattern data is rarely involved.
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 converting different unit cell distances as a result of numerical patterns of tropospheric radio 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;
original equal kilometer unit grid account final equal longitude and latitude singleSub-area S of bit-by-bit gridScore of originThe calculation is as follows:
Sscore of origin＝dOriginal weft dividing×dOriginal meridian point (1)
Wherein d isOriginal weft dividingThe latitude distance, d, of the final equal longitude and latitude unit single grid of the original equal kilometer unit gridOriginal meridian pointThe 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 SScore of originObtaining 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 latitudeRatio 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 occupation of：
SRatio of occupation of＝(SScore of origin/SFinally, the product is processed)×100％ (2)
Wherein S isFinally, the product is processedThe 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：
SFinally, the product is processed＝dWeft yarn×dWarp beam (3)
Wherein d isWeft yarnThe final latitude distance of a single grid with equal latitude and longitude units, dWarp beamThe 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 gridsOriginal sourceSingle grid ratio S of final equal longitude and latitude unitRatio of occupation ofComponent W ofAccount 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 sourceSingle grid ratio S of final equal longitude and latitude unitRatio of occupation ofComponent W ofAccount for：
WAccount for＝WOriginal source×SRatio of occupation of (4)
Step 8, the component W obtained in the step 7 is processedAccount forAdding 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 latitudeFinally, 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 longitude and latitude unit grid data which are not equidistant can be converted into the longitude and latitude unit grid data which are equidistant, and the practical application of researching the environmental characteristics of troposphere electric waves in a large area, particularly in a global range is facilitated.
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 are not intended to limit the invention.
The method comprises the steps of obtaining finer space-time grid data of the tropospheric electric wave environment parameters based on a mesoscale atmospheric numerical mode, taking directly output grid data set in a grid distance kilometer unit as a source, and converting the space-time grid data into equal grid distance data set in a grid distance latitude and longitude unit through a certain conversion relation based on the steps of reconstructing grids, obtaining the area occupation ratio of a single grid with equal latitude and longitude units formed by original kilometer unit grids and the like, obtaining the component of the single grid occupation ratio of the tropospheric electric wave environment parameters of the original kilometer unit grids with equal longitude and latitude units and the like.
As is well known, the tropospheric radio wave environment data in the research area acquired based on the mesoscale atmospheric numerical mode are all a lot of grid data formed at equal intervals by d km (a fixed value preset by a user in the mode), and if 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 a great relationship, for example, the latitude grid distance of a data grid point is increased 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 it is obviously seen that 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 comprises the steps that a horizontal grid distance d km is preset by a user in a numerical mode, grid data with the d km (a fixed value preset by the user in the mode) being equal intervals output by a medium-scale 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 environmental parameter information to be researched, and four vertexes of a certain unit grid with equal longitude and latitude to be converted are determined according to the specific equal kilometer unit grid distance (if the grid distance d km) and the specific equal longitude and latitude unit grid distance (if the grid distance c DEG) to be converted of the troposphere electric wave environmental numerical mode result. After 4 vertexes of a grid with the same longitude and latitude unit 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 a plurality of original electric wave environment parameters of the grid with the same kilometer unit according to the area size of the plurality of original electric wave environment parameters of the grid with the same kilometer unit. 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 converted data is prone to large errors due to the non-uniform atmospheric level, which is too large;
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 row m) corresponding to 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 4 vertices A, B, C, D of the grid can be given, wherein the coordinates of the longitude and latitude (also the longitude and latitude representing the grid point) of the bottom left corner B of the grid point (lat1, long1) are most easily determined, the longitude and latitude (lat1, long1) is known by how many grids the grid point is separated from the vertex of the bottom left corner of the first grid point (also the longitude and latitude representing the first grid point) in the east-west direction and how many grids the grid point is in the north-south direction, i.e. the initial longitude and latitude of the bottom left corner are respectively added by several grid distances apart, and secondly the position of A, C, D is determined, and the coordinate position of (lat1 ° + c °, long1 °), (lat1 °, long1 ° + c ° + lat1 ° + c °, long1 ° + c °) if displayed in (latitude, longitude) format.
And the calculation method of the 4 vertex coordinate positions of other grids is analogized in turn. As shown in FIG. 3, if the ABCD grid with equal latitude and longitude units is used as the initial reference grid point, the latitude and longitude coordinates of the 4 vertexes A1, B1, C1 and D1 enclosed by the first 4 dotted lines on the right and left sides of the grid (ABCD grid in the figure) are easily determined, wherein the coordinates of the vertex B1 on the lower left corner and the vertex A1 on the upper left corner are the C point coordinate and the D point coordinate of the grid on the left side, respectively, and the coordinates are (lat1, long1+ C, (lat1+ C, long1+ C), and the coordinates of the vertex C1 on the lower right corner and the vertex D1 on the lower right corner are (lat1, long1+ 2C), (lat1+ C, long1+ 2C), and so on the left vertex, upper left vertex, right vertex and right vertex of the second grid (lat1 long 2 + C) (lat 3629 long, (lat1 ° + c °, long1 ° +2c °), (lat1 °, long1 ° +3c °), and (lat1 ° + c °, long1 ° +3c °).
As shown in FIG. 3, it is easy to determine the longitude and latitude coordinates of the 4 vertexes A ', B', C ', D' (see FIG. 3) of the grid surrounded by the first 4 dotted lines immediately above the grid A1B1C1D1, wherein the positions of the lower left vertex B ', lower right vertex C' are most easily determined, and are the A1 point coordinates and D1 point coordinates of the lower grid, respectively, (lat1 ° + C °, long1 °), (lat1 ° + C °, long1 ° + C °), and the coordinate positions of the upper left vertex and upper right vertex A ', D' are (lat1 ° +2C °, long1 °), (lat1 ° + C °, long1 ° + C °), and so on, the coordinates of the lower left vertex, upper left vertex, right vertex, upper right vertex 4 vertexes of the vertically above grid are (lat1 ° + long1 ° + long1 ° + 8925 ° + long C ° + 36 ° + 898 ° + long C ° + 27 °, long1 ° + c °), (lat1 ° +3c °, 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, 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-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 with equal longitude and latitude units shown in (lat1, long1) of fig. 3, the vertex is A, B, C, D, the area of the new grid is composed of partial areas of the original several grids with equal kilometers, the specific number of the original several grids with equal kilometers is related to the grid distance of the original grid and the converted grid distance, and the number of the original several grids with equal kilometers forming the new grid can be determined according to the size and range of the longitude and latitude and the simple geometric relationship and 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 gridScore of originThe calculation is as follows:
Sscore of origin＝dOriginal weft dividing×dOriginal meridian point (1)
Wherein d isOriginal weft dividingForming final equal grade for original equal kilometer unit gridLatitude distance of the sub-area of a single grid of latitude and longitude units, dOriginal meridian pointThe 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 dOriginal weft dividingDistance d from longitudeOriginal meridian pointWith (d)Original weft dividing，dOriginal meridian point) Formally, latitude distance d of 1 to 9 parts areaOriginal weft dividingDistance d from longitudeOriginal meridian pointThen it is:
(latitude 3-lat1, longitude 1-long1), (latitude 4-latitude 3, longitude 3-long1), (lat1 ° + c ° -latitude 4, longitude 3-long1), (lat1 ° + c ° -latitude 4, longitude 4-long 3), (latitude 4-latitude 3, longitude 4-long 3), (latitude 3-lat1, longitude 4-long 3), (lat1 ° + c ° -latitude 4, long1+ c ° -long 4), (latitude 4-latitude 3, long1+ c ° -longitude 4), (latitude 3-lat1, long1+ c ° -longitude 4).
According to the formula (1):
Sscore 1(latitude 3-lat1) × (longitude 1-long1)
SScore 2(latitude 4-latitude 3) × (longitude 3-long1)
SScore 3(lat1+ c-latitude 4) × (longitude 3-long1)
SScore 4(lat1+ c-latitude 4) × (longitude 4-longitude 3)
SScore 5(latitude 4-latitude 3) × (longitude 4-longitude 3)
SScore 6(latitude 3-lat1) × (longitude 4-longitude 3)
SScore 7(lat1+ c-latitude 4) × (long1+ c-longitude 4)
SScore 8(latitude 4-latitude 3) × (long1+ c-longitude 4)
SScore 9(latitude 3-lat1) × (long1+ 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 unitScore of originObtaining the area ratio S of the single grids which are formed by the original grids of the same kilometer unit and have the same longitude and latitude unitRatio of occupation of。
Based on (2) the unit grids with equal longitude and latitude can be obtainedArea ratio S of original equal kilometer unit single gridRatio of occupation of：
SRatio of occupation of＝(SScore of origin/SFinally, the product is processed)×100％ (2)
Wherein SScore of originIs the sub-area of the original equal kilometer unit grid occupying the final equal longitude and latitude unit single grid, SFinally, the product is processedThe 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：
SFinally, the product is processed＝dWeft yarn×dWarp beam (3)
Wherein d isWeft yarnThe final latitude distance of a single grid with equal latitude and longitude units, dWarp beamThe 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):
Sratio of 1＝(SScore 1/SFinally, the product is processed)×100％
SRatio of 2＝(SScore 2/SFinally, the product is processed)×100％
SRatio of 3＝(SScore 3/SFinally, the product is processed)×100％
SRatio of 4＝(SScore 4/SFinally, the product is processed)×100％
SRatio of 5＝(SScore 5/SFinally, the product is processed)×100％
SRatio of 6＝(SScore 6/SFinally, the product is processed)×100％
SRatio of 7＝(SScore 7/SFinally, the product is processed)×100％
SRatio of 8＝(SScore 8/SFinally, the product is processed)×100％
SRatio of 9＝(SScore 9/SFinally, the product is processed)×100％
Step 7, forming an equal longitude and latitude list according to the obtained original equal kilometer unit gridsArea ratio S of single lattice of bitsRatio of occupation ofObtaining troposphere electric wave environment parameter W of original kilometer unit gridsOriginal sourceSingle grid ratio S of final equal longitude and latitude unitRatio of occupation ofComponent W ofAccount 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 sourceSingle grid ratio S of final equal longitude and latitude unitRatio of occupation ofComponent W ofAccount for：
WAccount for＝WOriginal source×SRatio of occupation of (4)
The environmental parameters 1 to N of troposphere electric wave corresponding to the original kilometer unit grids can be read from the data file as required, for the convenience of description, a certain parameter is read only, and certain parameter values corresponding to 1 to 9 original kilometer unit grids are respectively WOriginal 1、W2. sup. st、WOriginal 3、WOriginal 4、WOriginal 5、WOriginal 6、WOriginal 7、WOriginal 8、W9 original。
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:
Waccount for 1＝WOriginal 1×SRatio of 1
WAccount for 2＝W2. sup. st×SRatio of 2
WAccount for 3＝WOriginal 3×SRatio of 3
WAccount for 4＝WOriginal 4×SRatio of 4
WAccount for 5＝WOriginal 5×SRatio of 5
WAccount for 6＝WOriginal 6×SRatio of 6
WAccount for 7＝WOriginal 7×SRatio of 7
WAccount for 8＝WOriginal 8×SRatio of 8
WAccount for 9＝W9 original×SRatio of 9
Step 8, according to the acquired environmental parameter value W of the troposphere electric wave of one of the original kilometer unit gridsOriginal sourceSingle grid ratio S of final equal longitude and latitude unitRatio of occupation ofComponent W ofAccount forThe 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 latitudeFinally, 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：
WFinally, the product is processed＝WAccount for 1+WAccount for 2+WAccount for 3+WAccount for 4+WAccount for 5+WAccount for 6+WAccount for 7+WAccount for 8+WAccount for 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.
Embodiment 1, in this embodiment, 4 vertexes A, B, C, D of a new mesh in a final equal longitude and latitude unit (as shown in fig. 4) surrounded by 4 dotted lines in the figure are determined, and two electric wave environment parameters, namely 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 30 km. 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 acquires the 30km (fixed value preset by the user in the mode) output by the mesoscale atmospheric numerical mode as the equal kilometer unit grid data with equal intervals, and the grid data file contains the unequal interval longitude and latitude information corresponding to each grid data and the troposphere electric wave environment to be researchedParameter information, data sequence is arranged as figure 5, data in 1-9 columns in the data table are respectively latitude, longitude, altitude, longitude number of grid, south-north grid column number under each longitude of grid, air pressure value, air temperature value, water-steam mixing ratio and height value of atmospheric boundary layer, height information is not calculated for each grid in the data table, 4 electric wave environment parameters are provided, and are respectively air pressure value P and air temperature value tcWater-vapor mixing ratio qvAnd an atmospheric boundary layer height value PBLH. The first two columns 11.04795, 102.9431 in the first row represent the latitude and longitude corresponding to the first grid (row 1, column 1) of the numeric pattern 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 include many parameters of the electric wave environment, and the number of the parameters and the data format are not limited as long as the longitude and latitude of the corresponding grid point and the electric wave environment parameters 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 is 30km) of the troposphere electric wave environment numerical mode result and the specific equal longitude and latitude unit grid distance (the grid distance is 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 °, and 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 longitude and latitude units, for example, to convert into a grid in longitude and latitude (11.5 degrees and 103 degrees), and the grid distance is 0.5 degrees, 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 of the vertex B are 11.5 degrees of north latitude and 103 degrees of east longitude, and the coordinates of the vertex B are A, C, D degrees, and the coordinate positions of the vertex B are respectively 11.5 degrees +0.5 degrees of north latitude, 103 degrees of east longitude, 11.5 degrees of north latitude, 103 degrees of east longitude +0.5 degrees of north longitude, 11.5 degrees of north latitude +0.5 degrees of east longitude, 103 degrees of east longitude +0.5 degrees of east longitude, and the determination methods of the 4 vertex coordinate positions of other grids are analogized in sequence. 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, 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 well known that it can be determined how many of the 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 of the original meshes from 1 to 6, respectively (11.322 °, 102.9431 °), (11.59579 °, 102.9431 °), (11.86932 °, 102.9431 °) (11.86932 °, 103.2225 °) (11.59579 °, 103.2225 °) (11.322 °, 103.2225 °), where the meshes are identified according to the (latitude, longitude) format, where the latitude and longitude are the lower left vertex of the corresponding 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 obtainedScore 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 gridOriginal weft dividingAnd a longitudinal distance dOriginal meridian pointSubtracting the nearest two adjacent latitude and longitude values, if the 6 latitude distances d of the division area of the original equal kilometer unit gridOriginal weft dividingDistance d from longitudeOriginal meridian pointWith (d)Original weft dividing，dOriginal meridian point) Formally expressed, latitude distance d of 1 to 6 areasOriginal weft dividingDistance d from longitudeOriginal meridian pointThen it is: (11.59579 ° -11.5 °, 103.2225 ° -103 °), (11.86932 ° -11.59579 °, 103.2225 ° -103 °), (12.0 ° -11.86932 °, 103.2225 ° -103 °), (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):
Sscore 1＝(11.59579-11.5)×(103.2225-103)＝0.0213
SScore 2＝(11.86932-11.59579)×(103.2225-103)＝0.0609
By analogy in the following way,
Sscore 3＝0.0291
SScore 4＝0.0363
SScore 5＝0.0760
SScore 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 unitScore of originObtaining 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 latitudeRatio 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：
SFinally, the product is processed＝dWeft yarn×dWarp beam＝0.5×0.5＝0.25
Wherein d isWeft yarnThe distance in latitude of the single grid in the final equal latitude and longitude unit is 0.5 DEG, dWarp beamThe longitude distance of the single grid, which is the final equal longitude and latitude unit, is also 0.5 deg. here.
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:
Sratio of 1＝(SScore 1/SFinally, the product is processed)×100％＝8.52％
SRatio of 2＝(SScore 2/SFinally, the product is processed)×100％＝24.36％
SRatio of 3＝(SScore 3/SFinally, the product is processed)×100％＝11.64％
SRatio of 4＝(SScore 4/SFinally, the product is processed)×100％＝14.52％
SRatio of 5＝(SScore 5/SFinally, the product is processed)×100％＝30.4％
SRatio of 6＝(SScore 6/SFinally, the product is processed)×100％＝10.64％
Where S isScore 1、SScore 2、SScore 3、SScore 4、SScore 5、SScore 6、SScore 7、SScore 8、SScore 9The 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 unitRatio of occupation ofAcquiring tropospheric electric wave environment parameters W of the original kilometer unit grid according to the formula (4)Original sourceSingle grid ratio S of final equal longitude and latitude unitRatio of occupation ofComponent W ofAccount 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 positionscThe parameters, air temperature values corresponding to 1-6 original kilometer unit grids are respectively 28.8、28.2、26.5、28.9、29.1、29.3。
Acquiring troposphere electric wave environment parameters W of original kilometer unit grids according to formula (4)Original sourceSingle grid ratio S of final equal longitude and latitude unitRatio of occupation ofComponent W ofAccount forHere, the air temperature value tcComponent W ofAccount for：
WAccount for 1＝WOriginal 1×SRatio of 1＝28.8×8.52％＝2.45°
WAccount for 2＝W2. sup. st×SRatio of 2＝28.2×24.36％＝6.87°
WAccount for 3＝WOriginal 3×SRatio of 3＝26.5×11.64％＝3.08°
WAccount for 4＝WOriginal 4×SRatio of 4＝28.9×14.52％＝4.20°
WAccount for 5＝WOriginal 5×SRatio of 5＝29.1×30.4％＝8.85°
WAccount for 6＝WOriginal 6×SRatio 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 gridOriginal sourceSingle grid ratio S of final equal longitude and latitude unitRatio of occupation ofComponent W ofAccount forAdding 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 unitscGrid value W ofFinally, 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：
WFinally, the product is processed＝WAccount for 1+WAccount for 2+WAccount for 3+WAccount for 4+WAccount for 5+WAccount for 6＝28.6°
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 step 7, reading any parameter of the tropospheric 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 this point, wherein the atmospheric boundary layer height values corresponding to 1 to 6 original kilometer unit grids are 98.9, 148.6, 456.6, 431.2, 165.2 and 99.8 respectively.
Acquiring troposphere electric wave environment parameters W of original kilometer unit grids according to formula (4)Original sourceSingle grid ratio S of final equal longitude and latitude unitRatio of occupation ofComponent W ofAccount for。
Component W of atmospheric boundary layer height value PBLHAccount for：
WAccount for 1＝WOriginal 1×SRatio of 1＝98.9×8.52％＝8.43
WAccount for 2＝W2. sup. st×SRatio of 2＝148.6×24.36％＝36.20
WAccount for 3＝WOriginal 3×SRatio of 3＝456.6×11.64％＝53.15
WAccount for 4＝WOriginal 4×SRatio of 4＝431.2×14.52％＝62.61
WAccount for 5＝WOriginal 5×SRatio of 5＝165.2×30.4％＝50.22
WAccount for 6＝WOriginal 6×SRatio 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：
WFinally, the product is processed＝WAccount for 1+WAccount for 2+WAccount for 3+WAccount for 4+WAccount for 5+WAccount 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 embodiment is to determine the 4 vertexes a1, B1, C1 and D1 (as shown in fig. 4) of a new grid which is surrounded by the first 4 dashed lines next to the right side of the horizontal of the ABCD grid to form equal latitude and longitude units, and to obtain the air temperature electric wave environment parameters of the new grid point.
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 30 km. 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 tcWater-vapor mixing ratio qvAnd an atmospheric boundary layer height value PBLH. The first two columns 11.04795, 102.9431 in the first row represent the latitude and longitude corresponding to the first grid (row 1, column 1) of the numeric pattern 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 contain many electric wave environment parameters, without limitation to the number of parameters and data format, as long as the corresponding grid can be obtainedLongitude and latitude of the lattice points and electric wave environment parameters.
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 30km) 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 °, and 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 equal kilometer unit to the grid of equal longitude and latitude unit is performed in the grid in the figure, for example, if the grid to be converted into the grid of longitude and latitude (11.5 ° and 103.5 °) is to be a grid with a grid distance of 0.5 °, 4 vertices a1, B1, C1 and D1 of the grid need to be obtained, wherein the vertex B1 is most easily located, and its coordinates are 11.5 ° north latitude, 103.5 ° east longitude, a1, C1 and D1, which represent the coordinates of the grid, and their coordinate positions are 11.5 ° north latitude +0.5 °, 103.5 ° east longitude +0.5 ° north latitude, 103.5 ° east longitude, and 4 vertex coordinate positions of other grids are analogized in turn. 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 the most important thing is to determine the electric wave environment parameters of the converted grid, 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 mesh (11.5 °, 103.5 °) shown in fig. 4, the vertices are a1, B1, C1, and D1, the area of the mesh 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 converted mesh distance, and according to the size and range of the longitude and latitude, the number of kilometer-equivalent unit meshes of the new mesh can be determined according to a simple geometric relationship. The total area of 9 original kilometers unit grids is composed, the original grids are (11.322 degrees, 103.2225 degrees), (11.59579 degrees, 103.2225 degrees, (11.86932 degrees, 103.2225 degrees), (11.86932 degrees, 103.5018 degrees), (11.59579 degrees, 103.5018 degrees), (11.322 degrees, 103.5018 degrees), (11.86932 degrees, 103.7812 degrees, (11.59579 degrees, 103.7812 degrees), (11.322 degrees, 103.7812 degrees) from 1 to 9, the grids are identified according to (latitude, longitude) format, and the latitude and the longitude are corresponding to the lower left vertex of the grids.
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 obtainedScore 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 gridOriginal weft dividingAnd a longitudinal distance dOriginal meridian pointSubtracting the nearest two adjacent latitude and longitude values, if the 9 latitude distances d of the original equal kilometer unit grid sub-areaOriginal weft dividingDistance d from longitudeOriginal meridian pointWith (d)Original weft dividing，dOriginal meridian point) Form representation, from 1 toLatitude distance d of 9 divisionsOriginal weft dividingDistance d from longitudeOriginal meridian pointThen 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):
Sscore 1＝(11.59579-11.5)×(103.5018-103.5)＝1.72422×10-4
SScore 2＝(11.86932-11.59579)×(103.5018-103.5)＝4.9235×10-4
By analogy in the following way,
Sscore 3＝2.3522×10-4
SScore 4＝0.0365
SScore 5＝0.0764
SScore 6＝0.0268
SScore 7＝0.0286
SScore 8＝0.0598
SScore 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 unitScore of originObtaining 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 latitudeRatio 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：
SFinally, the product is processed＝dWeft yarn×dWarp beam＝0.5×0.5＝0.25
Wherein d isWeft yarnThe distance in latitude of the single grid in the final equal latitude and longitude unit is 0.5 DEG, dWarp beamFor final equal longitude and latitude unit single netThe longitudinal distance of the grid is 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):
Sratio of 1＝(SScore 1/SFinally, the product is processed)×100％＝0.0687％
SRatio of 2＝(SScore 2/SFinally, the product is processed)×100％＝0.1970％
SRatio of 3＝(SScore 3/SFinally, the product is processed)×100％＝0.0941％
SRatio of 4＝(SScore 4/SFinally, the product is processed)×100％＝14.6％
SRatio of 5＝(SScore 5/SFinally, the product is processed)×100％＝30.56％
SRatio of 6＝(SScore 6/SFinally, the product is processed)×100％＝10.72％
SRatio of 7＝(SScore 7/SFinally, the product is processed)×100％＝11.44％
SRatio of 8＝(SScore 8/SFinally, the product is processed)×100％＝23.92％
SRatio of 9＝(SScore 9/SFinally, the product is processed)×100％＝8.4％
Where S isScore 1、SScore 2、SScore 3、SScore 4、SScore 5、SScore 6、SScore 7、SScore 8、SScore 9The 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 unitRatio of occupation ofAcquiring tropospheric electric wave environment parameters W of the original kilometer unit grid according to the formula (4)Original sourceSingle grid ratio S of final equal longitude and latitude unitRatio of occupation ofComponent W ofAccount for。
The convection current corresponding to the original kilometer unit grids can be read from the data file as requiredLayer electric wave environment parameters 1 to 4, a parameter air temperature value t is read at the positioncThe 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 sourceSingle grid ratio S of final equal longitude and latitude unitRatio of occupation ofComponent W ofAccount for：
WAccount for 1＝WOriginal 1×SRatio of 1＝29.3×0.0687％＝0.020°
WAccount for 2＝W2. sup. st×SRatio of 2＝29.1×0.1970％＝0.057°
WAccount for 3＝WOriginal 3×SRatio of 3＝28.9×0.0941％＝0.027°
WAccount for 4＝WOriginal 4×SRatio of 4＝29.0×14.6％＝4.234°
WAccount for 5＝WOriginal 5×SRatio of 5＝29.2×30.56％＝8.924°
WAccount for 6＝WOriginal 6×SRatio of 6＝29.4×10.72％＝3.152°
WAccount for 7＝WOriginal 7×SRatio of 7＝28.7×11.44％＝3.283°
WAccount for 8＝WOriginal 8×SRatio of 8＝29.2×23.92％＝6.985°
WAccount for 9＝W9 original×SRatio 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 gridOriginal sourceSingle grid ratio S of final equal longitude and latitude unitRatio of occupation ofComponent W ofAccount forAdding 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 unitsFinally, 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：
WFinally, the product is processed＝WAccount for 1+WAccount for 2+WAccount for 3+WAccount for 4+WAccount for 5+WAccount for 6+WAccount for 7+WAccount for 8+WAccount 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 right side of the level 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 environmental parameter air temperature of the new grid point, which shows that the method of the invention is applicable and practical, and the converted data of the equidistant longitude and latitude units is convenient for the practical application of the electric wave environmental characteristics of the troposphere in large areas, especially the global range.

Claims (2)

1. A method for converting different unit grid distances of tropospheric radio 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;
sub-area S of original equal kilometer unit grid in final equal longitude and latitude unit single gridScore of originThe calculation is as follows:
Sscore of origin＝dOriginal weft dividing×dOriginal meridian point (1)
Wherein d isOriginal weft dividingThe latitude distance, d, of the final equal longitude and latitude unit single grid of the original equal kilometer unit gridOriginal meridian pointThe 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 SScore of originObtaining 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 latitudeRatio 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 occupation of：
SRatio of occupation of＝(SScore of origin/SFinally, the product is processed)×100％ (2)
Wherein S isFinally, the product is processedThe 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：
SFinally, the product is processed＝dWeft yarn×dWarp beam (3)
Wherein d isWeft yarnThe final latitude distance of a single grid with equal latitude and longitude units, dWarp beamThe 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 gridsOriginal sourceSingle grid ratio S of final equal longitude and latitude unitRatio of occupation ofComponent W ofAccount 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 sourceSingle grid ratio S of final equal longitude and latitude unitRatio of occupation ofComponent W ofAccount for：
WAccount for＝WOriginal source×SRatio of occupation of (4)
Step 8, the component W obtained in the step 7 is processedAccount forAdding 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 latitudeFinally, 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 true CN114417579A (en) 2022-04-29
CN114417579B 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网与天文大地网坐标系统的转换方法", 《测绘学院学报》 *

Also Published As

Publication number Publication date
CN114417579B (en) 2023-01-24

Similar Documents

Publication Publication Date Title
Pfenninger et al. Long-term patterns of European PV output using 30 years of validated hourly reanalysis and satellite data
CN111723482B (en) Satellite CO-based 2 Method for inverting surface carbon flux by column concentration observation
Nicks Stochastic generation of the occurrence, pattern, and location of maximum amount of daily rainfall
CN106054283B (en) A kind of method and device of inverting higher troposphere and lower stratosphere wind field
Bertotti et al. Wind and wave predictions in the Adriatic Sea
CN110909447B (en) High-precision short-term prediction method for ionization layer region
TW202105306A (en) A power generation amount prediction device, a power generation amount prediction method and program
CN105975791A (en) Sparse region rainfall estimation fusion method based on dual-smoothing method
CN112100922A (en) Wind resource prediction method based on WRF and CNN convolutional neural network
CN111292214A (en) Method for calculating high-precision Chinese offshore wave characteristic distribution
Costa et al. Development and validation of the Portuguese wind atlas
CN107067454B (en) Radar force diagram three-dimensional display method based on hypergraph component development platform
McGregor et al. Nested simulations of perpetual January climate over the Australian region
CN115203934A (en) Mountain area water-reducing downscaling method based on Logistic regression
CN109241212B (en) Retrieval method for historical rainfall based on mesoscale numerical atmospheric mode and high resolution
CN114417579B (en) Different unit grid distance conversion method for troposphere electric wave environment numerical mode result
Muis et al. Global projections of storm surges using high-resolution CMIP6 climate models: validation, projected changes, and methodological challenges
CN111582547B (en) Method for acquiring wind field distribution in different places by using wind field data set
He et al. Re-assessing vegetation carbon storage and emissions from land use change in China using surface area
CN116542029A (en) MDT manufacturing method and system based on ocean numerical mode and observation data assimilation
CN112541620B (en) Typhoon storm water increasing prediction method and system with high prediction precision and efficiency
Ko et al. A comparison of accuracy of the ocean thermal environments using the daily analysis data of the KMA NEMO/NEMOVAR and the US Navy HYCOM/NCODA
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
CN113743027A (en) Method and device for drawing wind resource map based on CFD technology

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