CN111246386A - Terminal positioning method and device - Google Patents
Terminal positioning method and device Download PDFInfo
- Publication number
- CN111246386A CN111246386A CN201811346262.8A CN201811346262A CN111246386A CN 111246386 A CN111246386 A CN 111246386A CN 201811346262 A CN201811346262 A CN 201811346262A CN 111246386 A CN111246386 A CN 111246386A
- Authority
- CN
- China
- Prior art keywords
- track
- points
- point
- data
- terminal
- 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
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W4/00—Services specially adapted for wireless communication networks; Facilities therefor
- H04W4/02—Services making use of location information
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W24/00—Supervisory, monitoring or testing arrangements
- H04W24/10—Scheduling measurement reports ; Arrangements for measurement reports
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W4/00—Services specially adapted for wireless communication networks; Facilities therefor
- H04W4/02—Services making use of location information
- H04W4/021—Services related to particular areas, e.g. point of interest [POI] services, venue services or geofences
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W4/00—Services specially adapted for wireless communication networks; Facilities therefor
- H04W4/02—Services making use of location information
- H04W4/029—Location-based management or tracking services
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W64/00—Locating users or terminals or network equipment for network management purposes, e.g. mobility management
Abstract
The embodiment of the invention provides a terminal positioning method and device. The method comprises the following steps: acquiring original measurement report data of a terminal to be measured, wherein the original measurement report data comprises track data of at least two first track points; according to the track data, second track points meeting the requirement of preset credibility probability are screened from the first track points, and the terminal moving trend direction is determined according to the track data of the second track points; according to the moving trend direction and a preset track fitting function, track correction is sequentially carried out on non-second track points in the original measurement report data, and a corrected reference point is obtained; and selecting a point which is within a preset distance range from the reference point and in the moving trend direction as a target track point of the non-second track point. The embodiment of the invention solves the problems of large positioning result error and disorder of single-user positioning tracks of the positioning algorithm based on MR data.
Description
Technical Field
The embodiment of the invention relates to the technical field of mobile communication, in particular to a terminal positioning method and device.
Background
With the development of mobile communication technology, Long Term Evolution (LTE) of universal mobile telecommunications technology has occupied a large amount of user market share with its superior characteristics. A time division Long Term Evolution (TD-LTE) system is a common system in LTE. Measurement is one of the important functions of the TD-LTE system, and the system needs to use the analysis result of the measurement data to realize the triggering of events such as cell selection, reselection and handover; the method can also be used for finding out problems in the network aiming at the statistical analysis of a large amount of measurement data. Compared with road test data, the measured data analysis result has the advantages of more comprehensively and completely reflecting network performance directly and the like.
Triggering of Measurement Reports (MRs) comprises two modes of event triggering and periodic triggering, and after Measurement is completed, periodic summary is carried out to generate Measurement report data files such as sample data (MRO) and statistical data (MRS). Typically, the period of the measurement report of the base station (eNodeB) or the terminal (UE) is set to 5120ms, i.e. periodic measurements are made every 5120 ms. The measurement report data mainly comes from the UE or eNodeB, and the measurement report generated in the radio resource management process.
In the network planning optimization work, considering that the UE periodically sends MR data, the actual distribution position of the UE has a corresponding relation with the position of a measuring point in the MR, so that MR positioning can be equivalent to UE positioning; by means of the UE positioning, the geographic rasterization presentation of MR distribution is realized, the current ubiquitous 4G weak coverage network problem can be positioned in a deep perspective mode by analyzing field information such as coverage performance contained in an MR, the outdoor and indoor deep coverage problems can be accurately positioned, and coverage evaluation and accurate planning are carried out according to the problems so as to support the later-stage 4G wireless network planning and fine optimization work.
The MR positioning accuracy is very important and is directly related to the precision of user behavior analysis under big data. Currently, in LTE, MR-based positioning methods mainly include the following:
1. a positioning method based on wireless signal propagation model positioning;
the position of a mobile station is estimated mainly according to the inverse proportion relation between the signal strength and the signal transmission distance in a Standard Propagation Model (SPM) in the process of wireless link budget. The specific implementation method comprises the steps of respectively detecting signals sent by transmitters through a plurality of receivers, collecting the signal intensity of the transmitters, calculating a signal transmission distance value by constructing a relational expression of the signal intensity and the transmission distance, and finally determining the position of the mobile station (terminal) by taking the distance value as a positioning parameter.
2. Timing Advance (TA) + Arrival direction (Angle of Arrival, AoA) positioning technology;
the TA + AOA is mainly for the receiving antenna to determine the direction of the mobile station according to the incident angle of the received signal, and then to estimate the position coordinate. Judging the arrival direction of the signals by using the incident angles of the detected multiple signals through the receiving antenna array; and estimating the position of the mobile station target according to the angle values of incoming waves in a plurality of directions.
3. Time Delay of Arrival (TDOA) location, i.e., Time Delay detection location:
the time delay detection positioning technology is a widely used technology in the current positioning system, and the principle of the time delay detection positioning technology is that a receiver establishes a positioning equation set by using a plurality of detected signal time delays, and then solution operation of a target position is carried out. The time delay detection is divided into two types, namely arrival time delay detection and arrival time delay difference detection, wherein the arrival time delay detection mode needs strict time synchronization of a network, and the complexity of system design is higher; the arrival time delay difference detection is an improvement on an arrival time delay detection mode, reduces the time synchronization requirement of a network, and is relatively easy for engineering realization.
4. An application service (Over The Top, OTT) positioning technology provided based on The internet;
a positioning technology for acquiring location information of a user based on application program (APP) software is called an OTT positioning technology, and the technical principle is to acquire longitude and latitude information when a terminal performs a HyperText Transfer Protocol (HTTP) service, and acquire the longitude and latitude information from a Uniform Resource Identifier (URI) of an HTTP Protocol header in information sent from the APP of the terminal. For example, the latitude and longitude information is obtained in a Payload (Payload) mode, and the Payload data sent to the terminal APP from the map server is decoded.
5. Hybrid positioning technology;
the hybrid location technology is to mix and use the above two or three location technologies, such as AOA, TDOA, location based on wireless signal propagation model, OTT location, etc., to detect and extract relevant location parameters for location integration operation.
6. Satellite positioning technology;
for example, GPS or AGPS requires a GPS/AGPS communication module for positioning, a GPS positioning function needs to be started in advance in a terminal, and the terminal measures and reports latitude and longitude information.
In the above positioning technologies, except the 6 th positioning technology that uses a satellite positioning technology to realize accurate positioning, other positioning methods or hybrid methods all realize MR positioning on the basis of terminal positioning, and evaluate the real situation of LTE whole-network wireless coverage by using MR data field information reported by UE.
However, the existing MR localization algorithm has the following disadvantages:
in the above mode 1, due to the complexity of the wireless signal propagation space, the uncertainty of the wireless channel is very strong, the distance estimation is prone to generate a deviation, and the positioning accuracy is difficult to meet the high-accuracy positioning requirement in the actual work.
Specifically, due to the complexity of the spatial propagation path of the wireless signal, such as various scattering, reflection, diffraction, and attenuation fading caused by blocking of various buildings, trees, and the like, the wireless signal of the base station rarely reaches the terminal through a straight propagation path, which results in an inevitable deviation of the distance calculated by using the wireless propagation model. Moreover, in terms of the propagation model formula, it is also difficult to find a set of parameter configurations suitable for all scenes, each scene has a certain difference, but even a slight difference will cause that it is difficult to achieve a high-precision effect in positioning. Due to the complexity of the actual environment, different propagation model parameters are suitable for different environments, so that the actual real environment needs to be combined when selecting the suitable propagation model and positioning method, including geographic terrain, building distribution, environmental characteristics, base station coverage density and other factors need to be comprehensively considered, and the propagation model and algorithm which are universal for each environment are difficult to select.
In the above 2 mode, the TA/AOA positioning method utilizes the TA value reported by the UE and the angle at which the wireless signal reaches the antenna to perform measurement and positioning. The AOA positioning technology requires the addition of an antenna array, which puts high requirements on the system design of the receiving end. Due to the complexity of the space environment, the blocking of buildings and the mobility of the terminal, the TA and AOA measurement are easy to deviate, and finally the MR positioning accuracy is greatly deviated. The positioning accuracy can theoretically reach 100-200 meters, however, due to the drift of base station signals, the single-point MR positioning has large random drift, which can cause hundreds of meters of positioning deviation, and can not meet the application requirement of high-precision positioning.
In the above 3 rd mode, since location information of APP needs to be acquired, on one hand, privacy of the user is unfavorable, and on the other hand, the APP is easily intercepted.
In the above 4 manner, there are three technical difficulties in OTT positioning at present, including: how to process massive operation data, resource investment and economic benefit need to be balanced; how to judge the used coordinate system needs to research a related matching judgment algorithm; and how to eliminate the interference information to ensure the accuracy of the longitude and latitude of the sampling point, and corresponding algorithms also need to be researched.
In the above-mentioned 5 th mode, the difficulty is increased by the difficulty of the above-mentioned 4 modes.
In the above-mentioned mode 6, although the GPS/AGPS positioning accuracy is the highest and can meet the working requirement, the terminal user often turns off the GPS positioning function of the terminal to avoid obtaining latitude and longitude information because of the requirement of protecting personal privacy.
In summary, the existing MR data-based localization algorithms have two significant drawbacks: firstly, the positioning result error caused by the complexity of a propagation space is large; secondly, the positioning azimuth deviation caused by signal drift has the characteristic of randomness, so that the single-user MR positioning track is out of sequence and disordered. Therefore, when the communication network construction planning and optimization is guided based on the result of the MR positioning algorithm, inaccurate description of coverage quality and coverage area is easily caused due to inaccurate positioning, which is likely to cause problems, and it is difficult to really and effectively guide the actual network planning and optimization work.
Disclosure of Invention
The embodiment of the invention provides a terminal positioning method and device, which are used for solving the problems of large positioning result error, and disorder of a single-user positioning track of a positioning algorithm based on MR data in the prior art.
In one aspect, an embodiment of the present invention provides a terminal positioning method, where the method includes:
acquiring original measurement report data of a terminal to be measured, wherein the original measurement report data comprises track data of at least two first track points;
according to the track data, second track points meeting the requirement of preset credibility probability are screened from the first track points, and the terminal moving trend direction is determined according to the track data of the second track points;
according to the moving trend direction and a preset track fitting function, track correction is sequentially carried out on non-second track points in the original measurement report data, and a corrected reference point is obtained;
and selecting a point with the distance between the reference point and the terminal in a preset distance range in the moving trend direction as a target track point of the non-second track point, wherein the target track point is a positioning point of the terminal.
In one aspect, an embodiment of the present invention provides a terminal positioning apparatus, where the apparatus includes:
the system comprises an acquisition module, a processing module and a processing module, wherein the acquisition module is used for acquiring original measurement report data of a terminal to be measured, and the original measurement report data comprises track data of at least two first track points;
the determining module is used for screening second track points meeting the requirement of preset credibility probability from the first track points according to the track data and determining the terminal moving trend direction according to the track data of the second track points;
the correction module is used for sequentially carrying out track correction on non-second track points in the original measurement report data according to the moving trend direction and a preset track fitting function to obtain a corrected reference point;
and the selection module is used for selecting the distance, the point of the reference point in a preset distance range in the moving trend direction is used as a target track point of the non-second track point, and the target track point is a positioning point of the terminal.
On the other hand, an embodiment of the present invention further provides an electronic device, which includes a memory, a processor, a bus, and a computer program stored in the memory and executable on the processor, where the processor implements the steps in the terminal positioning method when executing the program.
In still another aspect, an embodiment of the present invention further provides a non-transitory computer-readable storage medium, on which a computer program is stored, where the computer program, when executed by a processor, implements the steps in the above-mentioned terminal positioning method.
According to the terminal positioning method and device provided by the embodiment of the invention, original measurement report data of a terminal to be measured is obtained; according to the track data of the first track points of the original measurement report data, second track points meeting the requirement of preset credibility probability are screened from the first track points, and the terminal moving trend direction is determined according to the track data of the second track points; according to the moving trend direction and a preset track fitting function, track correction is sequentially carried out on non-second track points in the original measurement report data to obtain corrected reference points, and the deviation of the distance calculated based on a wireless propagation model is overcome; selecting a point which is within a preset distance range from the reference point and in the moving trend direction as a target track point of the non-second track point, wherein the target track point is a final positioning point of the terminal; comprehensively analyzing the terminal moving track from user dimension, time dimension, space dimension and data dimension through original measurement report data, and correcting track points by establishing a moving trend direction to overcome positioning result errors caused by propagation space complexity; through carrying out a lot of screening to the track point to establish credibility probability based on the orbit characteristic parameter between with the preceding track point, screen the track point through the parameter characteristic with adjacent track point, avoid because of the location position deviation that signal drift leads to, and because of the randomness of signal, thereby lead to the disorder and the confusion of single user MR location orbit.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and those skilled in the art can also obtain other drawings according to the drawings without creative efforts.
Fig. 1 is a schematic flowchart of a terminal positioning method according to an embodiment of the present invention;
FIG. 2 is one of schematic diagrams of a first example of an embodiment of the invention;
FIG. 3 is a second schematic diagram of a first example of the embodiment of the invention;
FIG. 4 is a schematic diagram of a second example of embodiment of the present invention;
fig. 5 is a schematic structural diagram of a terminal positioning device according to an embodiment of the present invention;
fig. 6 is a schematic structural diagram of an electronic device according to an embodiment of the present invention.
Detailed Description
In order to make the technical problems, technical solutions and advantages of the present invention more apparent, the following detailed description is given with reference to the accompanying drawings and specific embodiments. In the following description, specific details such as specific configurations and components are provided only to help the full understanding of the embodiments of the present invention. Thus, it will be apparent to those skilled in the art that various changes and modifications may be made to the embodiments described herein without departing from the scope and spirit of the invention. In addition, descriptions of well-known functions and constructions are omitted for clarity and conciseness.
It should be appreciated that reference throughout this specification to "an embodiment" or "an embodiment" means that a particular feature, structure or characteristic described in connection with the embodiment is included in at least one embodiment of the present invention. Thus, the appearances of the phrase "in an embodiment" or "in an embodiment" in various places throughout this specification are not necessarily all referring to the same embodiment. Furthermore, the particular features, structures, or characteristics may be combined in any suitable manner in one or more embodiments.
In various embodiments of the present invention, it should be understood that the sequence numbers of the following processes do not mean the execution sequence, and the execution sequence of each process should be determined by its function and inherent logic, and should not constitute any limitation to the implementation process of the embodiments of the present invention.
In the embodiments provided herein, it should be understood that "B corresponding to a" means that B is associated with a from which B can be determined. It should also be understood that determining B from a does not mean determining B from a alone, but may be determined from a and/or other information.
Fig. 1 shows a flowchart of a terminal positioning method according to an embodiment of the present invention.
As shown in fig. 1, the method for positioning a terminal provided in the embodiment of the present invention specifically includes the following steps:
The raw measurement report data mainly comes from the physical layer, Radio Link Control (RLC) layer, and measurement report generated by calculation in the Radio resource management process of the UE and eNodeB.
The raw measurement report data may include the following: basic core data such as an MME identification Code (Mme Code), an MME group identification Code (Mme group ID), a UE-in-MME identification Code (Mme UeS1apId), a TimeStamp (TimeStamp) and an LTE Cell identification (Cell Id); wherein, the MME is a network node (Mobility Management Entity). Sample data (MRO) may be included in the raw measurement report data, see table 1 for example, with table 1 showing an example of an MRO sample XML data file.
Table 1:
the original measurement report data comprises track data of at least two first track points, wherein the track data is MR data of the track points; the time interval of two adjacent first track points on the time stamp meets the preset requirement.
Optionally, as a first example, as shown in fig. 2, the raw measurement report data may be reported to a radio access Network element Management System (OMC-R) through statistics calculation (counter accumulation and statistics) by the eNodeB, where the OMC-R stores the raw measurement report data in a form of statistical data and reports the statistical data to a Network Management System (NMS) through the OMC-R, or as shown in fig. 3, the eNodeB directly reports the raw measurement report data to the OMC-R, and the OMC-R stores the raw measurement report data in a form of sample data.
And step 102, according to the track data, screening second track points meeting the requirement of preset credibility probability from the first track points, and determining the terminal moving trend direction according to the track data of the second track points.
Extracting track characteristic parameters of each track point from the track data, and determining the reliability probability of the track point according to the track characteristic parameters, namely determining the reliability probability according to the track characteristic parameters of the track point; the trajectory characteristic parameters may include displacement, time difference, moving speed, etc. from a previous trajectory point determined from the timestamp, in the time dimension.
The reliability probability is used for determining the reliability of the trace points; performing preset operation on the track characteristic parameters, and comprehensively considering multiple dimensions such as displacement, time difference, moving speed and the like to obtain reliability probability; and filtering the trace points which do not meet the requirement of the preset credibility probability in the first trace points to obtain second trace points.
Alternatively, the predetermined confidence probability requirement may be that the confidence probability is within a predetermined confidence interval or above a predetermined confidence threshold.
After the second track point is obtained through screening, determining the moving trend direction of the second track point according to a preset operation rule; alternatively, the movement tendency direction may be expressed in the form of a direction angle.
And 103, sequentially carrying out track correction on non-second track points in the original measurement report data according to the moving trend direction and a preset track fitting function to obtain a corrected reference point.
The preset track fitting function is used for fitting a track curve of the terminal; and according to the moving trend direction and a preset track fitting function, carrying out track correction on non-second track points in the original measurement report data, namely carrying out track correction on the original track points in the original measurement report data, so that the corrected reference points are positioned in the moving trend direction, and the deviation of the distance calculated based on the wireless propagation model is overcome.
And 104, selecting points which are within a preset distance range from the reference point and in the moving trend direction as target track points of the non-second track points, wherein the target track points are positioning points of the terminal.
After the reference point is determined, drawing a circle by taking the reference point as a circle center and a preset distance as a radius to obtain a circular distance range, and taking a point which is positioned in the preset distance range and in the moving trend direction as a target track point of the reference point, namely a target track point which is not the second track point, wherein the target track point is a final position point of the second track point; and after all the target track points are obtained, connecting all the target track points in series to obtain the moving track of the terminal.
In the above embodiment of the present invention, the original measurement report data of the terminal to be measured is obtained; according to the track data of the first track points of the original measurement report data, second track points meeting the requirement of preset credibility probability are screened from the first track points, and the terminal moving trend direction is determined according to the track data of the second track points; according to the moving trend direction and a preset track fitting function, track correction is sequentially carried out on non-second track points in the original measurement report data to obtain corrected reference points, and the deviation of the distance calculated based on a wireless propagation model is overcome; selecting a point which is within a preset distance range from the reference point and in the moving trend direction as a target track point of the non-second track point, wherein the target track point is a final positioning point of the terminal; comprehensively analyzing the terminal moving track from user dimension, time dimension, space dimension and data dimension through original measurement report data, and correcting track points by establishing a moving trend direction to overcome positioning result errors caused by propagation space complexity; by screening the track points for multiple times, establishing reliability probability based on the track characteristic parameters between the track points and the previous track point and screening the track points through the parameter characteristics of the adjacent track points, the positioning azimuth deviation caused by signal drift and the disorder and disorder of the single-user MR positioning track caused by the randomness of signals are avoided; the embodiment of the invention solves the problems of large positioning result error, and disorder of single-user positioning tracks of the positioning algorithm based on MR data in the prior art.
Optionally, in this embodiment of the present invention, the step of obtaining raw measurement report data of the terminal to be measured includes:
acquiring original measurement report data of a terminal to be measured, wherein the original measurement report data comprises track data of at least two original track points, and the track data at least comprises a timestamp;
and screening at least two first track points from the original track points according to the timestamp, and establishing a first track point set of the terminal.
Acquiring original measurement report data of a terminal to be measured, wherein the original measurement report data comprises track data of an original track point, and the track data comprises a timestamp; according to the 3GPP protocol specification, defining a terminal user unique identifier UserId, wherein the UserId is used for marking a resource identifier distributed by the current session of the terminal on a network side; alternatively, the UserId may be of the form:
UserId＝MmeGroupId+MmeCode+MmeUeS1apId；
wherein, the combination data of UserId is all digital type, and can form 8 bytes of an unsigned type of UINT64 by displacement, such as: UserId ═ MmeGroupId | mmeecode (upper 32-47bit) | mmeuees 1apId (lower 0-31 bit);
or combined into a concatenated string, such as:
UserId＝MmeGroupId|MmeCode|MmeUeS1apId。
after the UserId is determined, establishing a single-user LTE-MR data set (MRDataSet) with the UserId as an index, arranging the data sets according to a time stamp sequence, and carrying a sampling time stamp by track data of each original track point; each track data in the set is attached to a different cell CellId.
Traversing and positioning a data set with the terminal UserId as an index; and acquiring an LTE-MR data set of a user, and accessing each data sampling point according to a time sequence.
In the embodiment of the invention, the problem of multiplexing conflict in the range of the MME is avoided by setting the UserId, the attribution of the user MR sample of the multiplexing ID is judged according to the UserId and the MR acquisition sample period, the positioning error caused by the serial connection of different user MR data is avoided, and the positioning error caused by the overlapping of the user data is effectively avoided. And all MR acquisition samples are labeled with tissue data according to UserId while being ordered in time series for subsequent analysis.
And further, after all the original measurement report data of the terminal are collected, at least two first track points are screened from the original track points according to the time stamp, and a first track point set of the terminal is established.
Specifically, the step of screening at least two first track points from the original track points according to the timestamp and establishing a first track point set of the terminal includes:
if the original measurement report data comprises at least two adjacent region measurement report MR data of the original track point according to the timestamp, determining a corresponding first track point of the original track point according to the two adjacent region MR data and a preset positioning method;
and screening the first track points with the time intervals between adjacent points meeting the preset interval requirement according to the time stamps, and establishing a first track point set of the terminal.
If the original measurement report data comprises at least two adjacent area MR data aiming at a timestamp, determining longitude and latitude information of the original track point according to the two adjacent area MR data and a preset positioning method, namely, the track data of the first track point is more than the track data of the original track point by the longitude and latitude data;
specifically, the step of determining a first trace point corresponding to the original trace point according to the MR data of the two adjacent regions and a preset positioning method includes: and determining the longitude and latitude information of the first track point according to a TA/AOA positioning method or an arc positioning method.
Determining longitude and latitude information of a first track point according to an arc positioning method, wherein 3 distances between an original track point and a base station and between base stations of two adjacent areas are calculated according to a wireless signal propagation model (COST231-HATA model or SPM standard propagation model), and 3 circles are drawn by taking the points of the longitude and latitude of the 3 base stations as circle centers and the corresponding distances as radii respectively; and the intersection points of the 3 circles are the longitude and latitude information of the first track point. If the current original track point does not meet the positioning requirement, continuing to process the next LTE-MR data sample; and (4) processing each original track point in the set one by one for positioning, namely acquiring the longitude and latitude of each MR.
After determining the longitude and latitude information, performing time difference isolation on first track points, specifically, screening the first track points with the time interval between adjacent points meeting the preset interval requirement, including detecting the time stamp interval between the adjacent first track points;
as an example, obtaining the reporting MR period setting of the LTE-MR system:
Tmr{ms120，ms240，ms480，ms640，ms1024，ms2048，ms5120，ms10240，min1，min6，min12，min30，min60}；
the period is specified in 3GPP technical specification 36.133, and eNodeB or UE implements the period according to the technical specification requirement, and calculates the difference deltaT between Tdiff and MR period setting.
Ideally, the difference deltaT is 0; when the deltaT is greater than a preset timeout threshold (Tdiff _ MAX ═ 2 × Tmr can be set), judging that the LTE-MR data of the current user is finished, namely the connection of the current user is finished, indicating that the user reestablishes the connection, and entering the step of establishing the LTE-MR data set of the user when a new user LTE-MR starts;
that is, if the time interval between adjacent points is less than or equal to the preset timeout threshold, it indicates that the two first trace points are both trace points of the terminal. Through the step, the data conflict caused by multiplexing of the UserId in the MME equipment is solved, and the accuracy and reliability of LTE-MR data based on the user track are ensured; and traversing all original track point data of the user, and establishing a user track set UserTrackSet with the UserId as an identifier, namely a first track point set.
The circular arc positioning method carries out positioning analysis based on a positioning algorithm of wireless signal propagation model positioning, determines the reliability probability of the positioning result by combining a positioning result scene, and forms the user reference point positioning reference set by setting the reference point of which the threshold filtering part meets the threshold requirement.
Optionally, in the embodiment of the present invention, the step of screening, according to the trajectory data, second trajectory points that satisfy a preset reliability probability requirement from the first trajectory points includes:
determining a track characteristic parameter of the first track point according to the track data;
determining the reliability probability of each first track point according to the track characteristic parameters and a preset reliability calculation formula;
and screening the first track points with the credibility probability meeting the preset credibility probability requirement as second track points.
The track characteristic parameters are displacement, time difference, moving speed and the like between the track characteristic parameters and a previous track point, and the previous track point is determined according to the timestamp and is the previous track point on the time dimension.
And calculating the track characteristic parameters according to a preset reliability calculation formula, determining the reliability probability of each first track point, and screening the first track points with the reliability probabilities meeting the requirement of the preset reliability probability according to the reliability probability as second track points.
The preset reliability probability requirement may be that the reliability probability is within a preset reliability range or higher than a preset reliability threshold.
The preset confidence interval and the preset confidence threshold are determined according to the current moving state scene of the terminal.
Specifically, based on a first track point set, determining a CELL identifier of a current service CELL of each first track point by combining working parameter data of a base station to which a terminal belongs, and identifying whether the current service CELL is an indoor sub-base station or an outdoor macro-station; and judging a current movement state scene (CTXT) of the terminal by combining the movement speed and the base station type, wherein the CTXT comprises a room division scene, a fixed area, low-speed movement, urban vehicle-mounted movement, expressway movement and the like.
The scene displacement distance parameter Dn (unit: meter) is determined according to the moving state scene, such as the set room division D0, the fixed area D1, the low-speed moving D2, the urban vehicle-mounted moving D3 and the expressway moving D4.
In one aspect, a confidence interval [ Dn +/-Cn ] of the motion state scene CTXT at 95% confidence can be solved according to a preset confidence interval calculation formula, where Cn is a confidence interval threshold value.
On the other hand, a preset reliability threshold value PDmax can be set to be 0.5, the reliability probability of each first track point is calculated through a pd (x) function, high-reliability track points meeting the reliability requirement can be filtered according to the preset reliability threshold value, and a PDTrackPoints track point set, that is, a set of second track points, is obtained.
Specifically, the trajectory feature parameters at least include: displacement and time difference between the first track point and a previous first track point adjacent in time, and moving speed of the previous first track point;
the step of determining the reliability probability of each first track point according to the track characteristic parameters and a preset reliability calculation formula comprises the following steps:
determining the reliability probability of each first track point according to the following formula:
PD＝sigmoid{1/(|DS-UV0*TDS|/Cn)}；
wherein PD is the reliability probability, DS is the displacement, TDS is the time difference, UV0 is the moving speed, and Cn is a preset confidence interval threshold or a preset reliability threshold, which is a preset constant.
Calculating the displacement DS between two adjacent first track points according to the longitude and latitude data of the first track points, and calculating the time difference TDS according to the MR data; and calculating the moving speed UV of the current first track point according to the displacement DS and the time difference TDS.
And calculating the reliability probability according to the PD formula, wherein the sigmoid threshold function is defined as the following form:
if a preset reliability threshold value PDmax is set to be 0.5, calculating the reliability probability of each first track point through a PD (PD) (X) function, filtering high-reliability track points meeting the reliability requirement according to the set threshold value, and obtaining a PDTrackPoints track point set, namely a second track point set.
Optionally, in this embodiment of the present invention, the step of determining the terminal movement trend direction according to the trajectory data of the second trajectory point includes:
and presetting and sampling the second track points to obtain a second track point set, and determining the moving trend direction of the terminal on the second track points according to the timestamps and the displacements of the second track points in the second track point set.
Selecting a second track point as a user track point, sampling (sampling function S (DS, TDS, UV)) for continuous N adjacent sample points, ensuring uniform characteristics of sample time and space distribution, and establishing a second track point set; from the time stamps and the displacements, a movement trend direction DIR is determined at each sample point, which can be expressed in the form of a direction angle.
In the embodiment, the motion trail analysis model is established based on the MR data of the adjacent track points, the user motion and displacement model can be established due to the continuity of the terminal motion trail and the fixed characteristic of the MR sample report time interval, and the user motion direction and the motion displacement speed can be determined by combining the MR data context of the terminal trail in a specific sample analysis scene.
Further, in the embodiment of the present invention, the step of sequentially performing track correction on the non-second track points in the original measurement report data according to the moving trend direction and a preset track fitting function includes:
and according to a preset track fitting function, sequentially determining reference points of the non-second track points in the moving trend direction, wherein the reference points are the corrected reference points of the non-second track points.
According to a preset track fitting function, sequentially determining reference points of the non-second track points in the moving trend direction, wherein the reference points are corrected reference points of the non-second track points; optionally, for the non-second trace point MR data point positioning calculation, a least square interpolation calculation method can be adopted, and the following trace fitting function is established by combining the state data (DS, TDS, UV, CTXT, DIR) of the adjacent previous MR data point user:
TRACK＝K0+K1*UV+K2*CTXT+K3*log(DS)+K4*log(TDS)
wherein K0-K4 are preset parameters, and TRACK is the positioning result.
A reference point RefPoint is determined which combines displacement, velocity, direction in the direction of movement,
if K0-K4 is unknown, parameters K0, K1, K2, K3 and K4 can be solved by adopting a least square method aiming at the established TRACK fitting function TRACK, the sample data selected as a reference point is used as an input parameter to establish an equation for the TRACK function, and the value of K0, K1, K2, K3 and K4 can be solved by solving the equation set.
For sample data points which are not second TRACK points, solving a positioning result TRACK of adjacent non-reference points by using a TRACK fitting function, wherein DS0 (displacement), TDS0 (time difference), UV0 (speed), CTXT0 (scene) and DIR0 (direction) in a formula are all adjacent previous MR sample point state data; and solving the track result of each subsequent MR sample point in sequence to complete the subsequent non-reference point data calculation.
After the coordinates of the reference point are determined, a preset distance is determined as a radius, for example, a radius area of 20m, with the reference point as a center, as a positioning target area of a preset distance range, as a second example, as shown in fig. 4, a point P reference point, a circular area with the point P as a center, is shown as the preset distance range.
In the above embodiment of the present invention, the original measurement report data of the terminal to be measured is obtained; according to the track data of the first track points of the original measurement report data, second track points meeting the requirement of preset credibility probability are screened from the first track points, and the terminal moving trend direction is determined according to the track data of the second track points; according to the moving trend direction and a preset track fitting function, track correction is sequentially carried out on non-second track points in the original measurement report data to obtain corrected reference points, and the deviation of the distance calculated based on a wireless propagation model is overcome; selecting a point which is within a preset distance range from the reference point and in the moving trend direction as a target track point of the non-second track point, wherein the target track point is a final positioning point of the terminal; comprehensively analyzing the terminal moving track from user dimension, time dimension, space dimension and data dimension through original measurement report data, and correcting track points by establishing a moving trend direction to overcome positioning result errors caused by propagation space complexity; through carrying out a lot of screening to the track point to establish credibility probability based on the orbit characteristic parameter between with the preceding track point, screen the track point through the parameter characteristic with adjacent track point, avoid because of the location position deviation that signal drift leads to, and because of the randomness of signal, thereby lead to the disorder and the confusion of single user MR location orbit.
With the above description of the terminal positioning method according to the embodiment of the present invention, a terminal positioning apparatus according to the embodiment of the present invention will be described with reference to the accompanying drawings.
Referring to fig. 5, an embodiment of the present invention provides a terminal positioning apparatus, where the terminal positioning apparatus includes:
an obtaining module 501, configured to obtain original measurement report data of a terminal to be measured, where the original measurement report data includes track data of at least two first track points;
a determining module 502, configured to screen, according to the trajectory data, second trajectory points that meet a preset reliability probability requirement from the first trajectory points, and determine a terminal movement trend direction according to the trajectory data of the second trajectory points;
a correction module 503, configured to sequentially perform trajectory correction on non-second trajectory points in the original measurement report data according to the movement trend direction and a preset trajectory fitting function, so as to obtain a corrected reference point;
and the selecting module 504 is configured to select a point, which is within a preset distance range and in the moving trend direction, of the distance reference point as a target track point of the non-second track point, where the target track point is a positioning point of the terminal.
Optionally, in this embodiment of the present invention, the obtaining module 501 includes:
the system comprises an acquisition submodule and a processing submodule, wherein the acquisition submodule is used for acquiring original measurement report data of a terminal to be measured, the original measurement report data comprises track data of at least two original track points, and the track data at least comprises a timestamp;
and the establishing sub-module is used for screening at least two first track points from the original track points according to the timestamp and establishing a first track point set of the terminal.
Optionally, in this embodiment of the present invention, the establishing sub-module is configured to:
if it is determined that the original measurement report data includes at least two neighbor measurement report MR data of the original trace point according to the timestamp, determining a corresponding first trace point of the original trace point according to the two neighbor MR data and a preset positioning device;
and screening the first track points with the time intervals between adjacent points meeting the preset interval requirement according to the time stamps, and establishing a first track point set of the terminal.
Optionally, in this embodiment of the present invention, the determining module 502 includes:
the first determining submodule is used for determining a track characteristic parameter of the first track point according to the track data;
the second determining submodule is used for determining the reliability probability of each first track point according to the track characteristic parameters and a preset reliability calculation formula;
and the screening submodule is used for screening the first track point with the reliability probability meeting the requirement of the preset reliability probability as a second track point.
Optionally, in an embodiment of the present invention, the track characteristic parameters at least include: displacement and time difference between the first track point and a previous first track point adjacent in time, and moving speed of the previous first track point;
the second determination submodule is configured to:
determining the reliability probability of each first track point according to the following formula:
PD＝sigmoid{1/(|DS-UV0*TDS|/Cn)}；
wherein PD is the confidence probability, DS is the displacement, TDS is the time difference, and UV0 is the moving speed.
Optionally, in this embodiment of the present invention, the determining module 502 includes:
and the third determining submodule is used for presetting and sampling the second track points to obtain a second track point set, and determining the moving trend direction of the terminal on the second track points according to the time stamps and the displacements of the second track points in the second track point set.
Optionally, in this embodiment of the present invention, the modification module 503 is configured to:
and according to a preset track fitting function, sequentially determining reference points of the non-second track points in the moving trend direction, wherein the reference points are the corrected reference points of the non-second track points.
In the above embodiment of the present invention, the obtaining module 501 obtains the original measurement report data of the terminal to be measured; the determining module 502 is used for screening second track points meeting the requirement of preset credibility probability from the first track points according to the track data of the first track points of the original measurement report data, and determining the terminal moving trend direction according to the track data of the second track points; the correction module 503 sequentially performs track correction on non-second track points in the original measurement report data according to the moving trend direction and a preset track fitting function to obtain a corrected reference point, and overcomes the deviation of the distance calculated based on a wireless propagation model; the selection module 504 selects a point which is within a preset distance range from the reference point and in the moving trend direction as a target track point of the non-second track point, wherein the target track point is a final positioning point of the terminal; comprehensively analyzing the terminal moving track from user dimension, time dimension, space dimension and data dimension through original measurement report data, and correcting track points by establishing a moving trend direction to overcome positioning result errors caused by propagation space complexity; through carrying out a lot of screening to the track point to establish credibility probability based on the orbit characteristic parameter between with the preceding track point, screen the track point through the parameter characteristic with adjacent track point, avoid because of the location position deviation that signal drift leads to, and because of the randomness of signal, thereby lead to the disorder and the confusion of single user MR location orbit.
Fig. 6 is a schematic structural diagram of an electronic device according to yet another embodiment of the present invention.
As shown in fig. 6, the electronic device may include: a processor (processor)610, a communication Interface (Communications Interface)620, a memory (memory)630 and a communication bus 640, wherein the processor 610, the communication Interface 620 and the memory 630 communicate with each other via the communication bus 640. The processor 610 may call logic instructions in the memory 630 to perform the following method:
acquiring original measurement report data of a terminal to be measured, wherein the original measurement report data comprises track data of at least two first track points;
according to the track data, second track points meeting the requirement of preset credibility probability are screened from the first track points, and the terminal moving trend direction is determined according to the track data of the second track points;
according to the moving trend direction and a preset track fitting function, track correction is sequentially carried out on non-second track points in the original measurement report data, and a corrected reference point is obtained;
and selecting a point with the distance between the reference point and the terminal in a preset distance range in the moving trend direction as a target track point of the non-second track point, wherein the target track point is a positioning point of the terminal.
In addition, the logic instructions in the memory 630 may be implemented in software functional units and stored in a computer readable storage medium when the logic instructions are sold or used as independent products.
In another embodiment of the present invention, a non-transitory computer-readable storage medium is provided, where a computer program is stored on the non-transitory computer-readable storage medium, and when the computer program is executed by a processor, the steps in the method provided in the foregoing embodiment of the present invention are implemented, and details of the implementation are not repeated.
Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes.
The above-described embodiments of the apparatus are merely illustrative, and the units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment. One of ordinary skill in the art can understand and implement it without inventive effort.
Through the above description of the embodiments, those skilled in the art will clearly understand that each embodiment can be implemented by software plus a necessary general hardware platform, and certainly can also be implemented by hardware. With this understanding in mind, the above-described technical solutions may be embodied in the form of a software product, which can be stored in a computer-readable storage medium such as ROM/RAM, magnetic disk, optical disk, etc., and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute the methods described in the embodiments or some parts of the embodiments.
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.
Claims (10)
1. A terminal positioning method, characterized in that the method comprises:
acquiring original measurement report data of a terminal to be measured, wherein the original measurement report data comprises track data of at least two first track points;
according to the track data, second track points meeting the requirement of preset credibility probability are screened from the first track points, and the terminal moving trend direction is determined according to the track data of the second track points;
according to the moving trend direction and a preset track fitting function, track correction is sequentially carried out on non-second track points in the original measurement report data, and a corrected reference point is obtained;
and selecting a point with the distance between the reference point and the terminal in a preset distance range in the moving trend direction as a target track point of the non-second track point, wherein the target track point is a positioning point of the terminal.
2. The method of claim 1, wherein the step of obtaining raw measurement report data of the terminal under test comprises:
acquiring original measurement report data of a terminal to be measured, wherein the original measurement report data comprises track data of at least two original track points, and the track data at least comprises a timestamp;
and screening at least two first track points from the original track points according to the timestamp, and establishing a first track point set of the terminal.
3. The method according to claim 2, wherein the step of creating a first set of trace points of the terminal by selecting at least two first trace points from the original trace points according to the timestamp comprises:
if the original measurement report data comprises at least two adjacent region measurement report MR data of the original track point according to the timestamp, determining a corresponding first track point of the original track point according to the two adjacent region MR data and a preset positioning method;
and screening the first track points with the time intervals between adjacent points meeting the preset interval requirement according to the time stamps, and establishing a first track point set of the terminal.
4. The method according to claim 1, wherein the step of selecting, according to the trajectory data, second trajectory points that satisfy a preset confidence probability requirement from the first trajectory points comprises:
determining a track characteristic parameter of the first track point according to the track data;
determining the reliability probability of each first track point according to the track characteristic parameters and a preset reliability calculation formula;
and screening the first track points with the credibility probability meeting the preset credibility probability requirement as second track points.
5. The method according to claim 4, wherein the trajectory feature parameters include at least: displacement and time difference between the first track point and a previous first track point adjacent in time, and moving speed of the previous first track point;
the step of determining the reliability probability of each first track point according to the track characteristic parameters and a preset reliability calculation formula comprises the following steps:
determining the reliability probability of each first track point according to the following formula:
PD＝sigmoid{1/(|DS-UV0*TDS|/Cn)}；
wherein PD is the confidence probability, DS is the displacement, TDS is the time difference, UV0 is the movement speed, and Cn is a preset confidence interval threshold.
6. The method according to claim 1, wherein the step of determining the terminal movement trend direction according to the trajectory data of the second trajectory point comprises:
and presetting and sampling the second track points to obtain a second track point set, and determining the moving trend direction of the terminal on the second track points according to the timestamps and the displacements of the second track points in the second track point set.
7. The method according to claim 1, wherein the step of sequentially performing trajectory correction on the non-second trajectory points in the raw measurement report data according to the moving trend direction and a preset trajectory fitting function includes:
and according to a preset track fitting function, sequentially determining reference points of the non-second track points in the moving trend direction, wherein the reference points are the corrected reference points of the non-second track points.
8. A terminal positioning apparatus, characterized in that the apparatus comprises:
the system comprises an acquisition module, a processing module and a processing module, wherein the acquisition module is used for acquiring original measurement report data of a terminal to be measured, and the original measurement report data comprises track data of at least two first track points;
the determining module is used for screening second track points meeting the requirement of preset credibility probability from the first track points according to the track data and determining the terminal moving trend direction according to the track data of the second track points;
the correction module is used for sequentially carrying out track correction on non-second track points in the original measurement report data according to the moving trend direction and a preset track fitting function to obtain a corrected reference point;
and the selection module is used for selecting the distance, the point of the reference point in a preset distance range in the moving trend direction is used as a target track point of the non-second track point, and the target track point is a positioning point of the terminal.
9. An electronic device, comprising a memory, a processor, a bus and a computer program stored on the memory and executable on the processor, the processor implementing the steps in the terminal positioning method according to any of claims 1 to 7 when executing the program.
10. A non-transitory computer-readable storage medium having stored thereon a computer program, characterized in that: the program, when executed by a processor, implements the steps in a terminal positioning method as claimed in any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811346262.8A CN111246386B (en) | 2018-11-13 | 2018-11-13 | Terminal positioning method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811346262.8A CN111246386B (en) | 2018-11-13 | 2018-11-13 | Terminal positioning method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111246386A true CN111246386A (en) | 2020-06-05 |
| CN111246386B CN111246386B (en) | 2021-06-08 |
Family
ID=70863606
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811346262.8A Active CN111246386B (en) | 2018-11-13 | 2018-11-13 | Terminal positioning method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111246386B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111968398A (en) * | 2020-07-14 | 2020-11-20 | 深圳市综合交通运行指挥中心 | Method, device, terminal and medium for determining running state of public transport means |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101173980A (en) * | 2007-11-21 | 2008-05-07 | 湖南大学 | Indoor node locating algorithm based on ultra-broadband |
| US20100331664A1 (en) * | 2009-06-30 | 2010-12-30 | Joachim Graessner | Automatic positioning of a slice plane in mr angiography measurements |
| CN103776488A (en) * | 2012-04-02 | 2014-05-07 | 玛珂系统分析和开发有限公司 | Position determining device |
| CN108513251A (en) * | 2018-02-13 | 2018-09-07 | 北京天元创新科技有限公司 | A kind of localization method and system based on MR data |
-
2018
- 2018-11-13 CN CN201811346262.8A patent/CN111246386B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101173980A (en) * | 2007-11-21 | 2008-05-07 | 湖南大学 | Indoor node locating algorithm based on ultra-broadband |
| US20100331664A1 (en) * | 2009-06-30 | 2010-12-30 | Joachim Graessner | Automatic positioning of a slice plane in mr angiography measurements |
| CN103776488A (en) * | 2012-04-02 | 2014-05-07 | 玛珂系统分析和开发有限公司 | Position determining device |
| CN108513251A (en) * | 2018-02-13 | 2018-09-07 | 北京天元创新科技有限公司 | A kind of localization method and system based on MR data |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111968398A (en) * | 2020-07-14 | 2020-11-20 | 深圳市综合交通运行指挥中心 | Method, device, terminal and medium for determining running state of public transport means |
| CN111968398B (en) * | 2020-07-14 | 2022-07-29 | 深圳市综合交通运行指挥中心 | Method, device, terminal and medium for determining running state of public transport means |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111246386B (en) | 2021-06-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9829563B2 (en) | Positioning method and apparatus | |
| US9467803B2 (en) | Detecting multipath and determining positioning measurement uncertainty | |
| RU2390791C2 (en) | Positioning for wlan and other wireless networks | |
| CN108702585B (en) | Hybrid fingerprinting/OTDOA location techniques and systems | |
| US9137772B2 (en) | System and method for mass calibration of radio frequency fingerprint (RF-FP) location measurements | |
| EP2640116B1 (en) | Calibration method and device for coverage database | |
| CN102869038B (en) | Measurement compensation method and device in positioning of base station, server and system | |
| CN107925980B (en) | Method and device for positioning based on signal correlation function characteristic feedback | |
| JP2010511155A (en) | Method, software and system for determining the position of a user device | |
| CN103582115A (en) | Positioning method, control equipment and mobile communication system | |
| US9479902B1 (en) | Determining a propagation-time adjustment for a wireless coverage area, based on information provided by wireless terminals | |
| CN107197518B (en) | Method and device for determining terminal position | |
| US20180279251A1 (en) | Fingerprint Positioning For Mobile Terminals | |
| JP7382487B2 (en) | Method and device for on-demand positioning related application data | |
| CN103686998A (en) | Switching position positioning method and device based on AOA | |
| WO2016112689A1 (en) | Angle of arrival and time advance correction processing method and device | |
| US20100110181A1 (en) | Passive Positioning Information Of a Camera In large Studio Environment | |
| CN111246386B (en) | Terminal positioning method and device | |
| WO2018121439A1 (en) | Method and apparatus for determining direct path | |
| CN105338487B (en) | Techniques for multi-pass geolocation | |
| KR102416604B1 (en) | Mehtod and device for precise positioning for wireless communication system | |
| Wu et al. | DorFin: WiFi fingerprint-based localization revisited | |
| US10401468B2 (en) | Single positioning controller and positioning control system | |
| CN106796277B (en) | Location adjustment in a mobile communication network | |
| CN1753546A (en) | Location signaling strengthens between radio network controller and stand-alone service mobile positioning center |
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 |