Reference definition method, device and equipment suitable for regional observation station and storage medium
Patent Information
- Application Number
- CN202510722283.9
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2025-05-30
- Publication Date
- 2025-10-17
Smart Images
Figure CN120804472A_ABST
Abstract
Description
TECHNICAL FIELD
[0001] The present application relates to the field of geodetic surveying, and in particular to a reference definition method and device suitable for regional survey stations, equipment and storage medium. BACKGROUND
[0002] The reference definition problem is a key problem in space geodetic surveying, which refers to the system rank deficiency problem when estimating the station coordinates and the earth reference frame parameters simultaneously. Due to the seven degrees of freedom of the earth reference frame, including three translations, three rotations and one scale, the normal equation set is not full rank and cannot be directly solved. The significance of reference definition is to solve this mathematical ill-posed problem, establish a stable and consistent earth reference frame, and ensure the continuity and consistency of coordinate results obtained at different times and by different technologies. In addition, correct reference definition is crucial for separating different physical processes in the earth system (such as crustal movement and earth movement), and provides a foundation for global navigation, crustal deformation monitoring and geodynamics research. Current reference definition methods mainly focus on global stations, and common methods include loose constraint method, tight constraint method and minimum constraint method.
[0003] The minimum constraint method is currently the preferred method for maintaining the International Terrestrial Reference Frame (ITRF), which only applies necessary degree of freedom constraints, introduces an equal number of constraint conditions as the rank deficiency of the normal equation, and does not interfere with the reference information of the ground station network itself.
[0004] Existing reference definition methods are mainly suitable for global station networks. The reference definition model is highly sensitive to station distribution, and has outstanding problems when applied to regional networks. Although the loose constraint method is simple to implement, the subjective weight selection leads to insufficient stability of the framework, making it difficult to meet the needs of modern high-precision earth science research. The tight constraint method excessively restricts the natural changes of stations, and when there are unidentified abnormalities or local deformations in the reference stations, the errors will be propagated to the entire network, distorting the true crustal deformation signal. The widely recognized minimum constraint method, although theoretically the most reasonable, its effective implementation is heavily dependent on the spatial distribution, number and quality of reference stations, which is often difficult to meet in regional networks.
[0005] There are essential differences between regional station network and global station network in many aspects, which make it face serious challenges when directly applying the reference definition method of global framework to regional network. First of all, the geometric configuration difference: global network stations are distributed on the entire earth's surface, providing all-around geometric coverage, while regional network (for example, A country region) only covers a limited geographical range, the geometric configuration is poor, resulting in weaker direction constraints in reference frame definition. Therefore, when applying the reference definition method of global station network to regional station network, there is a mismatch problem, which leads to poor accuracy of regional station network. SUMMARY
[0006] The application provides a reference definition method, device and equipment suitable for regional stations and a storage medium to solve the problem of poor accuracy of regional station network in the prior art and improve the accuracy of regional station network.
[0007] The application provides a reference definition method suitable for regional stations, which comprises the following steps: based on the first historical calculation result of the global station network and the second historical calculation result of the regional station network, extracting the transformation characteristics of the preset Helmert parameter, and the transformation characteristics representing the conversion deviation between the first historical calculation result and the second historical calculation result; based on the preset Helmert parameter, the transformation characteristics and the minimum constraint condition of the regional station network, constructing a target virtual observation equation set of the regional station network; based on the actual observation data of the regional station network and the target virtual observation equation set, obtaining a target reference definition suitable for the regional station, and based on the target reference definition, obtaining the earth reference frame parameter of the corresponding region of the regional station network, and the target virtual observation equation set is used to make the minimum constraint condition match the actual statistics of the parameters of the regional station network.
[0008] According to the reference definition method suitable for regional stations provided by the application, the first historical calculation result comprises first historical calculation results at different times, and the second historical calculation result comprises second historical calculation results at different times. Based on the first historical calculation result of the global station network and the second historical calculation result of the regional station network, the transformation characteristics of the preset Helmert parameter are extracted, which comprises the following steps: based on the preset Helmert parameter, constructing a Helmert transformation equation between the first historical calculation result and the second historical calculation result; substituting the first historical calculation result at different times and the second historical calculation result at different times into the Helmert transformation equation, and solving the preset Helmert parameter at different times; calculating the average value, variance, covariance and time-varying characteristics of the preset Helmert parameter at different times to obtain the transformation characteristics.
[0009] The method for defining the reference of the regional survey station provided by the application comprises the following steps: obtaining a virtual observation parameter corresponding to a minimum constraint condition; constructing a first target virtual observation equation based on the minimum constraint condition and the virtual observation parameter; constructing a second target virtual observation equation based on a constraint condition of a preset Helmert parameter, the preset Helmert parameter and an average value of the preset Helmert parameter; and obtaining a target virtual observation equation set based on the first target virtual observation equation and the second target virtual observation equation.
[0010] The method for defining the reference of the regional survey station provided by the application further comprises the following steps after the target virtual observation equation set of the regional survey station network is constructed based on the preset Helmert parameter, the transformation feature and the minimum constraint condition of the regional survey station network: solving the target virtual observation equation set by using a constrained least square method to obtain a station coordinate of the virtual observation parameter and a value of the preset Helmert parameter; and evaluating the accuracy and stability of the station coordinate of the virtual observation parameter and the value of the preset Helmert parameter.
[0011] The preset Helmert parameter comprises an X-direction translation parameter, a Y-direction translation parameter, a Z-direction translation parameter, an X-direction rotation parameter, a Y-direction rotation parameter and a Z-direction rotation parameter.
[0012] The method for defining the reference of the regional survey station provided by the application further comprises the following steps after the target virtual observation equation set of the regional survey station network is constructed based on the preset Helmert parameter, the transformation feature and the minimum constraint condition of the regional survey station network: solving the target virtual observation equation set by using a constrained least square method to obtain a station coordinate of the virtual observation parameter and a value of the preset Helmert parameter; and evaluating the accuracy and stability of the station coordinate of the virtual observation parameter and the value of the preset Helmert parameter.
[0013] The minimum constraint condition comprises a no overall rotation condition constraint and a no overall translation condition constraint.
[0014] The application further provides a reference definition device suitable for a regional station, comprising: a feature extraction module, configured to extract a transformation feature of a preset Helmert parameter based on a first historical calculation result of a global station network and a second historical calculation result of a regional station network, the transformation feature representing a conversion deviation between the first historical calculation result and the second historical calculation result; a construction module, configured to construct a target virtual observation equation set of the regional station network based on the preset Helmert parameter, the transformation feature and a minimum constraint condition of the regional station network; and an observation module, configured to obtain a target reference definition suitable for the regional station based on actual observation data of the regional station network and the target virtual observation equation set, and obtain an earth reference frame parameter corresponding to a region of the regional station network based on the target reference definition, the target virtual observation equation set being used to make the minimum constraint condition match actual statistics of parameters of the regional station network.
[0015] The application further provides an electronic device, comprising a memory, a processor and a computer program stored in the memory and executable on the processor, and the processor implements the reference definition method suitable for the regional station according to any one of the above when executing the computer program.
[0016] The application further provides a non-transitory computer readable storage medium, which stores a computer program, and the computer program is executed by a processor to implement the reference definition method suitable for the regional station according to any one of the above.
[0017] The application further provides a computer program product, comprising a computer program, and the computer program is executed by a processor to implement the reference definition method suitable for the regional station according to any one of the above.
[0018] The reference definition method, device, equipment and storage medium suitable for the regional station provided by the application introduce historical information of a global station network into a regional station network through a transformation feature of a preset Helmert parameter, make up for the defect of poor geometry configuration of the regional station network, improve the matching of a reference definition model applied to the regional station network, and further improve the observation accuracy of the regional station network. BRIEF DESCRIPTION OF DRAWINGS
[0019] In order to more clearly illustrate the technical solutions of the present application or the prior art, the following will briefly introduce the drawings needed to be used in the embodiments or prior art description. Obviously, the drawings described below are only some embodiments of the present application, and the ordinary skilled in the art can obtain other drawings according to these drawings without any creative effort.
[0020] Figure 1 is a flowchart of the reference definition method suitable for regional stations provided by the present application.
[0021] Figure 2 is a flowchart of constructing a target virtual observation equation set provided by the present application.
[0022] Figure 3 is a structural schematic diagram of the reference definition device suitable for regional stations provided by the present application.
[0023] Figure 4 is a structural schematic diagram of the electronic device provided by the present application. DETAILED DESCRIPTION
[0024] In order to make the objects, technical solutions and advantages of the present application clearer, the technical solutions of the present application will be described clearly and completely below with reference to the drawings in the present application. Obviously, the described embodiments are only some embodiments of the present application, not all embodiments. Based on the embodiments in the present application, all other embodiments obtained by the ordinary skilled in the art without any creative effort belong to the protection scope of the present application.
[0025] The minimum constraint method is the preferred method in the maintenance of the International Terrestrial Reference Frame (ITRF) at present, which only applies necessary degree of freedom constraints, introduces constraint conditions equivalent to the rank deficiency of normal equation, and does not interfere with the reference information of the ground station network itself. The minimum constraint conditions mainly include: No Net Rotation (NNR), No Net Translation (NNT) and No Net Scale (NNS), which are used to define the rotation, translation and scale reference of the ground station network, respectively. In the maintenance of the Earth Reference Frame, NNT and NNR constraints are often added when estimating the Earth rotation parameters, and NNT constraints are often added when estimating the Earth motion. The advantage of the minimum constraint method is to preserve the internal geometry of the station network, maximize the reflection of the true crustal deformation information, and avoid the distortion of the framework caused by excessive constraints. This method is realized by constructing virtual observation equations, and is the standard method for high-precision Earth Reference Frame maintenance.
[0026] When applying minimum constraints, it is necessary to select core stations in the station network. The core station network method is an important strategy for implementing benchmark definition. This method is based on the construction of a benchmark network by selecting core stations with uniform distribution, high observation quality, and stable time series. For different space geodetic technologies, the corresponding international service agencies will recommend a list of high-quality stations suitable for benchmark definition. These core stations should avoid stations affected by earthquakes, instrument anomalies, or local deformations. In practical applications, virtual observation equations are constructed for these core stations, and NNR or NNT is applied. The benchmark definition achieved through these core stations not only ensures the stability of the reference frame, but also effectively reflects the dynamic changes of the earth. The selection of the core station network is crucial and directly affects the quality and reliability of the reference frame. Therefore, the performance of the core stations needs to be evaluated regularly, and the core station list needs to be updated as needed.
[0027] As the primary means of implementing benchmark definitions, the effectiveness of the core station network approach is highly dependent on the quality and distribution of core stations. The core stations recommended by international service agencies are relatively evenly distributed at the global scale, but significant deviations may occur at the regional scale. In particular, when high-quality stations are scarce or unevenly distributed within a region, relying solely on a limited number of core stations within the region to implement benchmark definitions can easily lead to systematic biases and instabilities in the solution. While combining regional and global networks can leverage the stability of the global framework, it also introduces errors and uncertainties from the global framework and makes it difficult to highlight regionally specific deformation characteristics.
[0028] The following combination Figures 1-4 The present invention describes a datum definition method, device and electronic equipment applicable to regional survey stations.
[0029] Figure 1 This is a flow chart of the benchmark definition method applicable to regional stations provided by the present invention. Figure 1 As shown, the benchmark definition method applicable to regional stations includes steps S100 to S300, and each step is described in detail as follows.
[0030] S100: extracting transformation features of preset Helmert parameters based on a first historical solution result of a global station network and a second historical solution result of a regional station network.
[0031] The transformation feature represents a conversion deviation between the first historical solution result and the second historical solution result.
[0032] The region of the present invention includes a country or a province. The global station network is used to observe the movement of the entire earth. The regional station network is used to observe and solve the parameters of the earth reference frame, including the earth's rotation parameters, the motion of the center of the earth, etc.
[0033] The first historical calculation result of a global station network and the second historical calculation result of a regional station network are obtained. For example, the Solution Independent Exchange Format (SINEX) files of the global network weekly or daily calculation results officially released by the International Global Navigation Satellite System Service (IGS) are collected to obtain the first historical calculation result. The time span of the first historical calculation result is not less than 5 years to fully cover the seasonal change cycle. At the same time, the calculation results of the selected regional (for example, domestic) and surrounding regional network in the corresponding format are collected. In the data preprocessing link, all SINEX files are unified in the coordinate system (for example, converted to International Terrestrial Reference Frame (ITRF2014 or ITRF2020)), epoch reduction, antenna phase center consistency check, special attention is paid to the coordinate jump caused by major earthquakes and the systematic deviation caused by equipment replacement, and obvious abnormal values are removed to obtain the second historical calculation result. Subsequently, a high-quality core station set is selected based on strict standards: the observation history length is more than 5 years, the annual effective observation week number is not less than 40 weeks, the weekly solution coordinate repeatability is better than 3mm (plane) and 5mm (elevation), the linear trend is clear and there is no significant nonlinear change. In the global station network, the spatial distribution of the core station should cover the main tectonic plates, especially to ensure that the corresponding plate of the regional station network has enough representative stations. In the regional station network, the core station should be appropriately densified to ensure that the main tectonic units (for example, plateau, plain, coastal area) of the corresponding plate of the regional station network have representative stations.
[0034] For the global core station network and the regional station network after screening, the preset Helmert transformation of each SINEX calculation result relative to the ITRF reference frame is calculated, and the transformation features of the Helmert parameters are extracted. The transformation features represent the conversion deviation between the first historical calculation result and the second historical calculation result.
[0035] S200: Based on the preset Helmert parameters, the transformation features, and the minimum constraint conditions of the regional station network, a set of target virtual observation equations of the regional station network is constructed.
[0036] The minimum constraint conditions include: no overall rotation condition constraint and no overall translation condition constraint.
[0037] The minimum constraint conditions include a no net rotation (NNR) constraint and a no net translation (NNT) constraint. A first target virtual observation equation is constructed according to the minimum constraint conditions, and a second target virtual observation equation is constructed according to preset Helmert parameters and transformation characteristics. A target virtual observation equation set is obtained according to the first target virtual observation equation and the second target virtual observation equation.
[0038] According to the no net rotation constraint and the no net translation constraint, the minimum constraint conditions are determined, the overall rotation deviation and the overall translation deviation of the reference definition of the global station network applied to the regional station network are mainly considered, and the stability of the regional station network is ensured.
[0039] S300: Based on the actual observation data of the regional station network and the target virtual observation equation set, a target reference definition suitable for the regional station is obtained, and the earth reference frame parameters of the corresponding region of the regional station network are obtained based on the target reference definition.
[0040] The target virtual observation equation set is used to match the minimum constraint conditions with the actual statistics of the parameters of the regional station network.
[0041] The target virtual observation equation set is solved to obtain constraint information of the virtual observation parameters. Station coordinates of the actual observation data are obtained. Based on the station coordinates of the actual observation data and the constraint information (target reference definition) of the virtual observation parameters, the earth reference frame parameters (for example, the earth rotation parameters) of the corresponding region of the regional station network are solved.
[0042] The reference definition method suitable for the regional station provided by the embodiment introduces the historical information of the global station network into the regional station network through the transformation characteristics of the preset Helmert parameters, makes up for the defects of the poor geometry of the regional station network, improves the matching of the reference definition model applied to the regional station network, and further improves the observation accuracy of the regional station network. The preset Helmert parameter constraint is further considered on the basis of the conventional minimum constraint condition, the target virtual observation equation set is constructed, the systematic error caused by the global crustal change that cannot be reflected by the regional station network is eliminated, and more reliable reference constraint information is provided for solving the earth rotation parameters of the regional station.
[0043] Based on the above embodiment, the first historical calculation result includes first historical calculation results at different times, and the second historical calculation result includes second historical calculation results at different times. Based on the first historical calculation result of the global station network and the second historical calculation result of the regional station network, the transformation characteristics of the preset Helmert parameters are extracted, including the following steps: Constructing a Helmert transformation equation between the first historical solution result and the second historical solution result based on preset Helmert parameters; Substitute the first historical solution results at different times and the second historical solution results at different times into the Helmert transformation equation to solve the preset Helmert parameters at different times; The mean, variance, covariance and time-varying characteristics of the preset Helmert parameters at different times are calculated to obtain the transformation features.
[0044] In one embodiment, the preset Helmert parameters include an X-direction translation parameter, a Y-direction translation parameter, a Z-direction translation parameter, an X-direction rotation parameter, a Y-direction rotation parameter, and a Z-direction rotation parameter.
[0045] The Helmert transformation equation is shown in formula (1).
[0046] Formula (1); in, 、 and The first historical solution for the global station network, The first in the global network of stations The horizontal coordinates of the measuring stations, The first in the global network of stations The vertical coordinate of each measuring station, The first in the global station network of a measuring station Axis coordinates, 、 and is the second historical solution result of the regional station network, The regional station network The horizontal coordinates of the measuring stations, The regional station network The vertical coordinate of each measuring station, The regional station network of a measuring station Axis coordinates, is the X-direction translation parameter, is the Y-direction translation parameter, is the Z-direction translation parameter, is the X-axis rotation parameter, is the Y-axis rotation parameter, The Z-axis rotation parameter.
[0047] According to the Helmert transformation equation, the complete six preset parameters Helmert transformation of the SINEX solution results (including the first historical solution and the second historical solution) of different epochs relative to the ITRF reference frame are calculated, the weighted least square estimation method is used, the SINEX coordinates of the core station are used as the observation values, and the preset Helmert parameters at different times are solved. Wherein, the weight matrix (the weight of the observation value) is constructed based on the variance-covariance information in the SINEX file. Then, statistical analysis is performed on the preset Helmert parameters at different times, the expected mean and standard deviation of the preset Helmert parameters at different times are calculated, the correlation coefficient matrix between the preset Helmert parameters is analyzed, the stability and seasonal variation characteristics of the preset Helmert parameter sequence at different times are evaluated, and finally the transformation characteristics of the preset Helmert parameters are obtained, including the average value of the preset Helmert parameters, the complete variance-covariance matrix (including the variance and covariance of the preset Helmert parameters) and the time-varying characteristics.
[0048] According to the first historical solution and the second historical solution, the transformation characteristics of the preset Helmert parameters are extracted, and then the target virtual observation equation set is constructed, which is beneficial to overcome the time series instability defects and the poor geometric configuration defects of the regional station network.
[0049] Based on the above embodiment, as shown in Figure 2 , based on the preset Helmert parameters, the transformation characteristics and the minimum constraint condition of the regional station network, the target virtual observation equation set of the regional station network is constructed, including the following steps: Obtaining the virtual observation parameters corresponding to the minimum constraint condition; Based on the minimum constraint condition and the virtual observation parameters, the first target virtual observation equation is constructed; Based on the constraint condition of the preset Helmert parameters, the preset Helmert parameters and the average value of the preset Helmert parameters, the second target virtual observation equation is constructed; Based on the first target virtual observation equation and the second target virtual observation equation, the target virtual observation equation set is obtained.
[0050] The target virtual observation equation set is as shown in formula (2).
[0051] Formula (2); Wherein, , , ,..., , and are virtual observation parameters, , , , , and is a preset Helmert parameter, , , , , and is an average value of the corresponding preset Helmert parameter, is a first target virtual observation equation, is a second target virtual observation equation.
[0052] Further, the weight of each virtual observation parameter is determined by the accuracy of the historical solution results (including the first historical solution result and the second historical solution result). The weight of the average value of the preset Helmert parameter can be dynamically adjusted according to the change of time.
[0053] According to the minimum constraint condition and the virtual observation parameter, the first target virtual observation equation is constructed; according to the constraint condition of the preset Helmert parameter, the preset Helmert parameter and the average value of the preset Helmert parameter, the second target virtual observation equation is constructed, on the basis of the traditional minimum constraint condition, the constraint of the preset Helmert parameter is further considered, the system error caused by the regional station network unable to reflect the global crustal change is eliminated, and a more reliable spatial reference is provided for the regional crustal deformation research.
[0054] Based on the above embodiment, the earth reference frame parameters of the corresponding region of the regional station network are obtained based on the target reference definition, including the following steps: The actual observation data are strictly preprocessed, including cycle slip repair, multi-path effect weakening, troposphere delay estimation and the like. The target virtual observation equation set is solved to obtain the constraint information of the virtual observation parameter. Based on the preprocessed actual observation data, the constraint information of the virtual observation parameter and the preset weight of the station coordinates of the virtual observation parameter, the earth reference frame parameters are solved.
[0055] The actual observation data of the regional station network are strictly preprocessed, including cycle slip repair, multi-path effect weakening, troposphere delay estimation and the like.
[0056] The preprocessed actual observation data are solved to obtain the station coordinates of the actual observation data. The target virtual observation equation set is solved by using the constrained least square method to realize the joint optimization estimation of the virtual observation parameter and the preset Helmert parameter, and the constraint information of the virtual observation parameter is obtained.
[0057] According to the pre-processed actual observation data, the constraint information of the virtual observation parameters, and the preset weight of the station coordinates of the virtual observation parameters, the regional station network is used to solve the earth reference frame parameters (for example, the earth rotation parameters).
[0058] According to the pre-processed actual observation data, the constraint information of the virtual observation parameters, and the preset weight of the station coordinates of the virtual observation parameters, the earth reference frame parameters are solved, the problems of poor geometry of the regional station network, uneven distribution of the stations, and complex construction environment are solved, and the stability and reliability of the regional station network are improved.
[0059] Based on the above embodiment, based on the preset Helmert parameters, the transformation characteristics, and the minimum constraint conditions of the regional station network, after constructing the target virtual observation equation set of the regional station network, the following steps are further included: The least square method with constraints is used to solve the target virtual observation equation set, and the station coordinates of the virtual observation parameters and the values of the preset Helmert parameters are obtained. The accuracy and stability of the station coordinates of the virtual observation parameters and the values of the preset Helmert parameters are evaluated.
[0060] The least square method with constraints is used to solve the target virtual observation equation set, and the station coordinates of the virtual observation parameters and the values of the preset Helmert parameters are obtained.
[0061] The quality of the calculation results (including the station coordinates of the virtual observation parameters and the values of the preset Helmert parameters) is comprehensively evaluated, for example, the consistency of the estimated preset Helmert parameters with the prior constraints is tested, the stability and noise characteristics of the station coordinate time series in repeated observations are analyzed, the calculation results of the present application are compared with the results of the traditional minimum constraint (NNR / NNT) method, the framework stability and parameter accuracy improvement degree of the regional station network are evaluated, and the fine expression ability of the framework of the regional station network of the present application to the earth reference frame parameters such as the earth rotation parameters is evaluated by comparing with known geodynamic phenomena such as earthquakes and tectonic deformation, and finally a high-quality regional reference frame which maintains the high consistency of the ITRF benchmark and accurately reflects the regional crustal deformation characteristics is formed.
[0062] The reference definition method suitable for regional stations of the present application is simple, adaptable, highly operable, and easy to integrate and apply in existing processing software and processes. Compared with new methods that need complex algorithm modification, the reference definition method suitable for regional stations of the present application only needs to adjust parameters on the basis of existing space geodetic measurement software, and can be realized. The increased calculation burden is very small (only 6-7 global parameters are added), and the processing efficiency is basically not affected. The core prior constraint information of the present application is derived from long-term accumulated historical data statistics, and is objective and reliable, avoiding the uncertainty caused by subjective parameter setting. At the same time, these prior constraints can be flexibly updated according to the new framework release or processing strategy changes, so that the reference definition method suitable for regional stations of the present application has good time adaptability. In addition, the reference definition method suitable for regional stations of the present application has good compatibility with other technical means (such as multi-technology combination, geophysical model constraint, etc.), can be used in cooperation with other advanced methods while maintaining its own advantages, and has good expansion potential. In the processing of abnormal events (such as earthquakes, instrument failures, etc.), compared with the core station fixed method, the present application is less affected by single-point abnormalities, runs more stably, has low maintenance cost, and is particularly suitable for frequent tectonic activity environments.
[0063] The reference definition device suitable for regional stations provided by the present application is described below, and the reference definition device suitable for regional stations described below can be correspondingly referred to the reference definition method suitable for regional stations described above.
[0064] As shown in Figure 3 A reference definition device suitable for regional stations includes a feature extraction module 301, a construction module 302, and an observation module 303.
[0065] The feature extraction module 301 is configured to extract transformation features of preset Helmert parameters based on first historical calculation results of a global station network and second historical calculation results of a regional station network, the transformation features characterizing conversion biases between the first historical calculation results and the second historical calculation results.
[0066] The construction module 302 is configured to construct a target virtual observation equation set of the regional station network based on preset Helmert parameters, transformation features, and minimum constraint conditions of the regional station network.
[0067] The observation module 303 is configured to obtain a target reference definition suitable for regional stations based on actual observation data of the regional station network and the target virtual observation equation set, and obtain earth reference frame parameters of a corresponding region of the regional station network based on the target reference definition, the target virtual observation equation set being configured to match the minimum constraint conditions with actual statistics of parameters of the regional station network.
[0068] The base definition device suitable for the regional station provided by the embodiment of the application introduces the historical information of the global station network into the regional station network through the transformation characteristics of the preset Helmert parameter, makes up for the defect of the poor geometry configuration of the regional station network, improves the matching of the base definition model applied to the regional station network, and further improves the observation accuracy of the regional station network. On the basis of the traditional minimum constraint condition, the application further considers the constraint of the preset Helmert parameter, constructs the target virtual observation equation set, eliminates the systematic error caused by the global crustal change which cannot be reflected by the regional station network, and provides more reliable base constraint information for solving the earth rotation parameter of the regional station.
[0069] In one embodiment, the first historical calculation result includes first historical calculation results at different time instants, the second historical calculation result includes second historical calculation results at different time instants, and the feature extraction module 301 is configured to: construct a Helmert transformation equation between the first historical calculation result and the second historical calculation result based on the preset Helmert parameter; substitute the first historical calculation result at different time instants and the second historical calculation result at different time instants into the Helmert transformation equation, and solve the preset Helmert parameter at different time instants; and calculate the average value, the variance, the covariance and the time-varying characteristics of the preset Helmert parameter at different time instants to obtain the transformation characteristics.
[0070] In one embodiment, the construction module 302 is configured to: obtain a virtual observation parameter corresponding to the minimum constraint condition; construct a first target virtual observation equation based on the minimum constraint condition and the virtual observation parameter; construct a second target virtual observation equation based on the constraint condition of the preset Helmert parameter, the preset Helmert parameter and the average value of the preset Helmert parameter; and obtain the target virtual observation equation set based on the first target virtual observation equation and the second target virtual observation equation.
[0071] In one embodiment, the construction module 302 is configured to: solve the target virtual observation equation set by using the least square method with constraint to obtain the station coordinates of the virtual observation parameter and the value of the preset Helmert parameter; and evaluate the accuracy and stability of the station coordinates of the virtual observation parameter and the value of the preset Helmert parameter.
[0072] In one embodiment, the preset Helmert parameter includes an X-direction translation parameter, a Y-direction translation parameter, a Z-direction translation parameter, an X-direction rotation parameter, a Y-direction rotation parameter and a Z-direction rotation parameter.
[0073] In one embodiment, the observation module 303 is configured to: perform cycle slip repair, multipath mitigation and troposphere delay estimation on the actual observation data to obtain preprocessed actual observation data; solve the preprocessed actual observation data; solve a set of target virtual observation equations to obtain constraint information of the virtual observation parameters; and solve the Earth reference frame parameters based on the preprocessed actual observation data, the constraint information of the virtual observation parameters and a preset weight of station coordinates of the virtual observation parameters.
[0074] In one embodiment, the minimum constraint condition comprises: a no overall rotation condition constraint and a no overall translation condition constraint.
[0075] Figure 4 An example of a schematic diagram of a physical structure of an electronic device is shown in Figure 4 As shown, the electronic device can include a processor 410, a communications interface 420, a memory 430 and a communications bus 440, wherein the processor 410, the communications interface 420 and the memory 430 can communicate with each other through the communications bus 440. The processor 410 can invoke a logical instruction in the memory 430 to execute a reference definition method suitable for regional stations, the method comprising: based on a first historical solution result of a global station network and a second historical solution result of a regional station network, extracting a transformation feature of a preset Helmert parameter, the transformation feature representing a conversion bias between the first historical solution result and the second historical solution result; based on the preset Helmert parameter, the transformation feature and a minimum constraint condition of the regional station network, constructing a set of target virtual observation equations of the regional station network; based on actual observation data of the regional station network and the set of target virtual observation equations, obtaining a target reference definition suitable for the regional stations, and based on the target reference definition, obtaining an Earth reference frame parameter corresponding to a region of the regional station network, the set of target virtual observation equations being used to make the minimum constraint condition match actual statistics of parameters of the regional station network.
[0076] In addition, the logic instructions in the memory 430 described above can be implemented in the form of a software function unit and sold or used as an independent product, and can be stored in a computer readable storage medium. Based on such understanding, the technical solutions of the present application essentially or the part that contributes to the prior art or part of the technical solutions can be embodied in the form of a software product. The computer software product is stored in a storage medium, and includes a plurality of instructions for causing a computer device (which can be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of the present application. The aforementioned storage medium includes: a U disk, a mobile hard disk, a read-only memory (ROM, Read-Only Memory), a random access memory (RAM, Random Access Memory), a magnetic disk or an optical disk, and various media that can store program codes.
[0077] In another aspect, the present application also provides a non-transitory computer readable storage medium, which stores a computer program, and the computer program is executed by a processor to implement the reference definition method for regional stations provided by the above-mentioned methods. The method comprises: based on a first historical solution result of a global station network and a second historical solution result of a regional station network, extracting a transformation feature of a preset Helmert parameter, the transformation feature characterizing a conversion bias between the first historical solution result and the second historical solution result; based on the preset Helmert parameter, the transformation feature and a minimum constraint condition of the regional station network, constructing a target virtual observation equation set of the regional station network; based on actual observation data of the regional station network and the target virtual observation equation set, obtaining a target reference definition suitable for the regional stations, and based on the target reference definition, obtaining an earth reference frame parameter corresponding to a region of the regional station network, the target virtual observation equation set is used to make the minimum constraint condition match the actual statistics of the parameters of the regional station network.
[0078] The device embodiments described above are only schematic, wherein the units illustrated as separate components can or can not be physically separated, and the components illustrated as units can or can not be physical units, i.e., they can be located in one place, or distributed on multiple network units. Part or all of the modules can be selected according to actual needs to achieve the purpose of the present embodiment. Those skilled in the art can understand and implement without creative labor.
[0079] Those skilled in the art can clearly understand the technical solutions of the various embodiments from the above description of the embodiments, and the various embodiments can be implemented by means of software with the necessary general hardware platforms, and of course, can also be implemented by hardware. Based on such understanding, the above technical solutions, essentially or in other words, the part of the prior art that makes a contribution, can be embodied in the form of a software product, which can be stored in a computer readable storage medium, such as a ROM / RAM, a magnetic disk, an optical disk, and the like, and includes a number of instructions for causing a computer device (which can be a personal computer, a server, or a network device, etc.) to execute the methods described in the various embodiments or some parts of the embodiments.
[0080] Finally, it should be noted that: the above embodiments are only used to illustrate the technical solutions of the present application, rather than limit them; although the present application has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that: it can still modify the technical solutions recorded in the foregoing embodiments, or make equivalent replacement for some technical features therein; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present application.
Claims
1. A datum definition method applicable to regional stations, characterized in that: include: Extracting a transformation feature of a preset Helmert parameter based on a first historical solution result of a global station network and a second historical solution result of a regional station network, wherein the transformation feature represents a conversion deviation between the first historical solution result and the second historical solution result; Constructing a target virtual observation equation set of the regional station network based on preset Helmert parameters, the transformation characteristics, and minimum constraints of the regional station network; Based on actual observation data of the regional station network and the target virtual observation equation set, a target datum definition applicable to the regional station is obtained, and based on the target datum definition, the earth reference frame parameters of the area corresponding to the regional station network are obtained. The target virtual observation equation set is used to match the minimum constraint conditions with the actual statistics of the parameters of the regional station network.
2. The method for defining a benchmark for a regional station according to claim 1, wherein: The first historical solution results include first historical solution results at different times, the second historical solution results include second historical solution results at different times, and the first historical solution results based on the global station network and the second historical solution results based on the regional station network are used to extract transformation features of preset Helmert parameters, including: constructing a Helmert transformation equation between the first historical solution result and the second historical solution result based on the preset Helmert parameters; Substituting the first historical solution results at different times and the second historical solution results at different times into the Helmert transformation equation to solve the preset Helmert parameters at different times; The mean, variance, covariance and time-varying characteristics of the preset Helmert parameters at different moments are calculated to obtain the transformation characteristics.
3. The method for defining a benchmark for a regional station according to claim 2, wherein: The target virtual observation equation set of the regional station network is constructed based on the preset Helmert parameters, the transformation characteristics and the minimum constraint conditions of the regional station network, including: Obtaining virtual observation parameters corresponding to the minimum constraint condition; Constructing a first target virtual observation equation based on the minimum constraint condition and the virtual observation parameter; Constructing a second target virtual observation equation based on a constraint condition of a preset Helmert parameter, the preset Helmert parameter, and an average value of the preset Helmert parameter; The target virtual observation equation set is obtained based on the first target virtual observation equation and the second target virtual observation equation.
4. The method for defining a benchmark for a regional station according to claim 3, wherein: After constructing the target virtual observation equation set of the regional station network based on the preset Helmert parameters, the transformation characteristics and the minimum constraint conditions of the regional station network, the method further includes: Solving the target virtual observation equation set using a constrained least squares method to obtain the station coordinates of the virtual observation parameters and the values of the preset Helmert parameters; Evaluate the accuracy and stability of the station coordinates of the virtual observation parameters and the values of the preset Helmert parameters.
5. The method for defining a benchmark for a regional station according to claim 1, wherein: The preset Helmert parameters include an X-direction translation parameter, a Y-direction translation parameter, a Z-direction translation parameter, an X-direction rotation parameter, a Y-direction rotation parameter, and a Z-direction rotation parameter.
6. The method for defining a benchmark for a regional station according to claim 1, wherein: The obtaining of earth reference frame parameters of the area corresponding to the regional station network based on the target datum definition includes: Performing cycle slip repair, multipath effect reduction, and tropospheric delay estimation on the actual observation data to obtain preprocessed actual observation data; Solving the target virtual observation equation set to obtain constraint information of virtual observation parameters; The earth reference frame parameters are solved based on the pre-processed actual observation data, the constraint information of the virtual observation parameters and the weights of the preset station coordinates of the virtual observation parameters.
7. The method for defining a benchmark for a regional station according to claim 1, wherein: The minimum constraint conditions include: no overall rotation constraint and no overall translation constraint.
8. A datum definition device suitable for regional survey stations, characterized in that: include: a feature extraction module, configured to extract a transformation feature of a preset Helmert parameter based on a first historical solution result of a global station network and a second historical solution result of a regional station network, wherein the transformation feature represents a conversion deviation between the first historical solution result and the second historical solution result; A construction module, configured to construct a target virtual observation equation set of the regional station network based on preset Helmert parameters, the transformation characteristics, and minimum constraints of the regional station network; An observation module is configured to obtain a target datum definition applicable to a regional station based on actual observation data of the regional station network and the target virtual observation equation set, and to obtain earth reference frame parameters of a region corresponding to the regional station network based on the target datum definition, wherein the target virtual observation equation set is configured to match the minimum constraint condition with actual statistics of the parameters of the regional station network.
9. An electronic device comprising a memory, a processor, and a computer program stored in the memory and running on the processor, characterized in that: When the processor executes the computer program, the benchmark definition method applicable to regional stations according to any one of claims 1 to 7 is implemented.
10. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the benchmark definition method applicable to regional stations according to any one of claims 1 to 7 is implemented.