WO2023135683A1 - 衛星信号受信装置、衛星信号選択方法、及びプログラム - Google Patents
衛星信号受信装置、衛星信号選択方法、及びプログラム Download PDFInfo
- Publication number
- WO2023135683A1 WO2023135683A1 PCT/JP2022/000721 JP2022000721W WO2023135683A1 WO 2023135683 A1 WO2023135683 A1 WO 2023135683A1 JP 2022000721 W JP2022000721 W JP 2022000721W WO 2023135683 A1 WO2023135683 A1 WO 2023135683A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- satellite
- signal
- satellite signal
- value
- satellite signals
- Prior art date
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/22—Multipath-related issues
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/28—Satellite selection
Definitions
- the present invention relates to technology for performing positioning and time synchronization with high accuracy by GNSS (Global Navigation Satellite System).
- GNSS Global Navigation Satellite System
- positioning and time synchronization processing is performed using GNSS satellite signals (hereinafter referred to as satellite signals) received by the GNSS antenna.
- satellite signals GNSS satellite signals
- the reception of satellite signals in the line-of-sight state may be blocked by structures, etc. that exist around the installation position of the GNSS antenna. In that case, the satellite signal is not received with the required signal strength at the GNSS antenna, or is received as an invisible satellite signal due to multipath reflection and diffraction by structures etc. existing around the installation position of the GNSS antenna. will be As a result, positioning performance and time synchronization performance by GNSS deteriorate.
- the present invention has been made in view of the above-mentioned points, and it is possible to appropriately select a satellite signal and perform positioning and time synchronization by GNSS with high accuracy even when the satellite signal reception environment is not good.
- the purpose is to provide technology.
- a signal selection unit that selects a predetermined number of satellite signals based on the reception quality of satellite signals received by a GNSS antenna; a measurement unit that performs positioning using a predetermined number of satellite signals selected by the signal selection unit and determines an initial estimated position;
- the signal selection unit is used for positioning or time synchronization based on a comparison result of comparing the reception time expected when the satellite signal is received as a direct wave at the initial estimated position and the actual reception time of the satellite signal.
- a satellite signal receiver is provided for selecting a plurality of satellite signals to receive.
- a technique that enables accurate positioning and time synchronization by GNSS even when the environment for receiving satellite signals is poor.
- FIG. 1 is a configuration diagram of a measuring device according to an embodiment of the present invention
- FIG. 4 is a flowchart of initial estimated position determination processing
- 4 is a flowchart of satellite selection processing
- 10 is a flowchart of a substitute satellite selection process
- FIG. 10 is a diagram illustrating an example of a processing procedure for selecting satellite signals
- FIG. 4 is a diagram showing an example of setting parameters
- It is a figure which shows the example of a setting of a GNSS bias value.
- FIG. 5 is a diagram showing an example of setting an elevation bias value; It is a figure which shows the example of a setting of a GNSS priority.
- FIG. 10 is a diagram showing a processing procedure for setting a bias value;
- FIG. 10 is a diagram showing a processing procedure for setting a bias value
- FIG. 4 is a diagram showing the maximum C/N 0 value of each group; It is a figure which shows the example of curve fitting.
- FIG. 5 is a diagram showing an example of bias value setting;
- FIG. 4 is a diagram showing how satellite signals are incident on and reflected from the walls of a building;
- FIG. 4 is a diagram showing how satellite signals are incident on and reflected from the walls of a building;
- FIG. 10 is a diagram showing an example of setting the dCN 0 value depending on the elevation angle;
- FIG. 10 is a diagram illustrating an example of a processing procedure for selecting satellite signals;
- FIG. 4 is a diagram showing an example of setting parameters;
- FIG. 10 is a diagram showing a processing procedure for setting a bias value;
- FIG. 5 is a diagram showing an example of bias value setting;
- FIG. 4 is a diagram showing an example of actual measurement of reception characteristics in the L1 band;
- FIG. 4 is a diagram showing an example of actual measurement of reception characteristics in the L2
- C A /N 0 mask scheme As a conventional method for eliminating invisible satellite signals, C A /N 0 mask scheme is known.
- the satellite signal C/ N0 value depends on the gain of the antenna, the reception sensitivity of the receiver, the cable loss between the antenna and the receiver, the type of satellite, etc., it is difficult to set the optimum threshold.
- Such interfering signals include intentionally generated GNSS jamming signals, equipment generated noise, and interfering signals from other communication systems.
- a satellite signal suitable for use in positioning and time synchronization is selected in a procedure described later assuming a multi-GNSS environment in which a large number of visible satellites can be secured.
- this procedure uses the strength (reception quality) of the received satellite signal as a basis for selection, it enables the selection of a satellite signal that eliminates the influence of the individual characteristics of the antenna and the receiver and the interference signal.
- visible satellite signals are preferentially selected, and if the number of visible satellites is small, in addition to visible satellites, invisible satellites with a small impact on accuracy deterioration and small propagation delay are selected.
- the C/N 0 value is used as the index of reception quality, but a reception quality index other than the C/N 0 value may be used.
- the "satellite signal" when "selecting the satellite signal” is assumed to be associated with the GNSS satellite that is the transmission source of the satellite signal. For example, if GNSS satellite A, GNSS satellite B, and GNSS satellite C are three different arbitrary GNSS satellites, selecting the three satellite signals means the satellite signal from GNSS satellite A, the satellite signal from GNSS satellite B satellite signals, and satellite signals from GNSS satellite C.
- FIG. 1 shows a configuration example of a measuring device 100 according to an embodiment of the present invention.
- the measurement device 100 includes a GNSS antenna 110, a signal reception unit 120, a signal selection unit 130, a measurement unit 140, an output unit 150, a signal data storage unit 160, a bias value setting unit 170, and a bias value storage unit 180.
- the measuring device 100 is a device that receives and processes satellite signals, and may be called a "satellite signal receiving device.”
- the GNSS antenna 110 receives radio waves transmitted from GNSS satellites in orbit and converts the radio waves into electrical signals.
- This electrical signal may be referred to as a "satellite signal”.
- the GNSS antenna 110 and the signal receiver 120 are connected by a cable, and satellite signals are sent to the signal receiver 120 by the cable.
- An amplifier may be provided between the GNSS antenna 110 and the signal receiver 120 when the distance between the GNSS antenna 110 and the signal receiver 120 is long.
- the signal receiving unit 120 receives a satellite signal, measures the C/N 0 value, and identifies the type of the GNSS satellite that is the transmission source of the received satellite signal. It also uses satellite orbital information (eg almanac, ephemeris) to measure the elevation angle. Satellite orbit information may be obtained from navigation messages of satellite signals, or may be obtained from other means (eg, a server on a network).
- the signal receiving section 120 sends the identification information (code such as the PRN number) of the received satellite signal, the elevation angle of the satellite signal, the C/N 0 value, and the satellite type to the signal selecting section 130 .
- the signal receiving unit 120 stores the identification information, elevation angle, C/N 0 value, and satellite type for each received satellite signal in the signal data storage unit 160 .
- the elevation angle is the angle between the line of sight and the horizontal plane when viewing the GNSS satellite that is the source of the satellite signal from the satellite signal reception point (that is, the GNSS antenna). For example, if a GNSS satellite is at the zenith, its elevation angle is 90°.
- GNSS satellites targeted in this embodiment are GPS, GLONASS, Galileo, BeiDou, QZSS, and NAVIC. However, these are examples and may be more or less than these types.
- the signal selection unit 130 selects satellite signals to be used for positioning and time synchronization from among the plurality of received satellite signals. The selection procedure will be described later.
- the measurement unit 140 performs time synchronization using a satellite signal transmitted from a GNSS satellite equipped with an atomic clock whose time is precisely managed with respect to absolute time, thereby performing time synchronization with high accuracy with respect to absolute time. Calculate the time information.
- the absolute time here is, for example, Coordinated Universal Time (UTC). Note that the measurement unit 140 may perform only one of positioning and time synchronization.
- the measurement unit 140 uses satellite signals from, for example, four or more GNSS satellites to code position the four parameters of the three-dimensional coordinate information (x, y, z) of the reception position and the time offset ( ⁇ t). Positioning and time synchronization are performed at the same time.
- the measurement unit 140 may perform carrier phase positioning in addition to code positioning.
- the measurement unit 140 outputs the time information based on this absolute time and the position information as the positioning result through the output unit 150 .
- the base station uses the time information synchronized with the absolute time, for example, the adjacent base station (synchronized with the absolute time) and TDD ( Time Division Duplex)
- TDD Time Division Duplex
- the bias value setting unit 170 sets (calculates) a bias value using the satellite signal data stored in the signal data storage unit 160 and stores the set bias value in the bias value storage unit 180 .
- the bias values stored in the bias value storage section 180 are used for satellite signal selection processing in the signal selection section 130 . Details of the bias value setting operation by the bias value setting unit 170 will be described later.
- the measuring device 100 may be a single device that is physically integrated, or may have some functional units that are physically separated, and a plurality of separated functional units are connected by a network. It may be a device that has been installed.
- the measurement device 100 may include all the functions shown in FIG. may be installed in the measuring device 100 and used.
- the satellite signal Selection and positioning operations may be performed on the cloud.
- the measurement unit 140 on the cloud returns the positioning calculation result to the output unit 150 .
- a device consisting of "the signal data storage unit 160 and the bias value setting unit 170" in the measurement device 100 may be provided on a network (for example, on a cloud), and the remaining functions may be installed in the measurement device 100 and used. .
- observation data is output from the signal receiving unit 120 provided in the measuring device 100, the observation data is stored in the signal data storage unit 160 provided on the cloud, and the bias value setting unit 170 provided on the cloud A bias value is set using the stored data.
- the bias value is returned from the bias value setting unit 170 on the cloud to the bias value storage unit 180 .
- the measuring device 100 performs satellite signal selection processing in the order of (1) initial estimated position determination processing, (2) satellite selection processing, and (3) substitute satellite selection processing.
- the outline of each process is as follows.
- Initial estimated position determination processing is processing for determining an initial estimated position that serves as a reference when selecting satellites.
- Satellite selection processing selects NLOS (Non Line Of Sight) satellite signals with short delays in addition to LOS (Line Of Sight) satellite signals by comparing the estimated reception time at the initial estimated position with the actual reception time. processing.
- Substitute satellite selection processing is processing to select additional satellites in order to suppress deterioration of positioning and time synchronization accuracy due to deterioration of DOP (Dilution Of Precision) value when the number of selected satellite signals is small. be.
- DOP Different Of Precision
- dCN 0 maximum value of the difference between the maximum C/N 0 value of the received satellite signal and the C/N 0 value of the satellite signal selected in the initial estimated position determination process N 0 : the lowest value selected in the initial estimated position determination process Number of satellite signals dT: Threshold value of delay time from reception time expected in initial estimated position in satellite selection process Ns : Maximum number of times initial estimated position is updated in satellite selection process dp: Completion judgment in satellite selection process N 1 : Minimum number of satellite signals selected in the process of selecting subsatellites ⁇ (1) Initial estimated position determination process> The initial estimated position determination process will be described with reference to the flowchart of FIG.
- the signal selector 130 extracts the satellite signal with the largest C/ N0 value from the plurality of satellite signals received by the signal receiver 120.
- the signal selection unit 130 selects satellite signals whose C/N 0 value difference from the maximum C/N 0 value extracted in S1-1 is equal to or less than a preset value (dCN 0 ). , to select from the plurality of satellite signals.
- the signal selection unit 130 determines the C/N 0 value. Satellite signals are interpolated from the larger (next-highest) satellite signal until the minimum number of satellite signals (N 0 ) is reached.
- the measuring unit 140 determines the initial estimated position of the measuring device 100 by code positioning using the satellite signals selected in S1-1 to S1-3.
- a positioning method other than code positioning may be used as the positioning method. The same applies to the step of performing code positioning in subsequent processing.
- Satellite selection processing will be described with reference to the flowchart of FIG.
- the signal selection unit 130 converts each satellite signal to the initial estimated position on the premise that the clock bias of the receiver (measuring device 100) and the time bias between satellite systems have been corrected in the code positioning process of S1-4. , the expected reception time and the actual reception time when received as a direct wave by the signal receiving unit 120 are compared.
- the expected reception time when a satellite signal is received as a direct wave can be calculated from the straight-line distance between the source satellite position of the satellite signal calculated from the orbital information and the initial estimated position.
- the actual reception time may be measured by any of the signal receiving section 120, the signal selecting section 130, and the measuring section 140.
- the signal selection unit 130 rejects satellite signals received with a delay of dT or more from the reception time at the initial estimated position where the actual reception time is expected, based on the comparison in S2-1.
- dT be X [ns]. Assuming that the expected reception time of satellite signal A is T and the actual reception time of satellite signal A is T+X+1 [ns], satellite signal A is rejected because the delay time is X+1 (dT or longer).
- the measuring unit 140 performs code positioning using a plurality of satellite signals received by the signal receiving unit 120, excluding the satellite signals rejected in S2-2. Update estimated location.
- the dT value may be decreased for each repetition.
- dT X [ns]
- n an integer equal to or greater than 2
- the signal selector 130 compensates for reaching this number.
- delay times are compared among a plurality of non-selected satellite signals that are candidates for selection.
- the delay time is "the difference between the expected reception time and the actual reception time when the satellite signal is received as a direct wave at the initial estimated position" to be compared with dT in S2-2.
- the value measured and calculated in S2-2 at the end of the iteration may be used, or the initial estimated position at the end of the calculation in S2-4 may be used in this S3- A value newly measured and calculated at time 1 may be used.
- the signal selection unit 130 selects the satellite signal with the smallest delay time in S3-1 and the satellite signal with the next smallest delay time in addition to the satellite signal selected at the end of the calculation in S2-4.
- the DOP values are compared in each case, and the satellite signal with the smaller cost value calculated from the delay time and the DOP value is selected as the satellite signal to be added.
- satellite signal A, satellite signal B, and satellite signal C are selected at the end of the calculation of S2-4.
- satellite signal D is selected as the satellite signal with the shortest delay time
- satellite signal E is selected as the satellite signal with the next shortest delay time in S3-2.
- the delay time of satellite signal D is 2
- the delay time of satellite signal E is 3
- the DOP value of "satellite signal A, satellite signal B, satellite signal C, satellite signal D” is 7, and "satellite signal A, satellite
- the DOP value of signal B, satellite signal C, and satellite signal E is 4.
- a satellite type premium may be set for the cost value.
- the higher the priority the smaller the numerical value of the priority
- the premium of priority 1 is set to 1
- the premium of priority 2 is set to 2
- the premium of priority 3 is set to 3
- the premium of priority 4 is set to 4
- the premium of priority 5 is set to 5.
- satellite signal D is added. selected as the satellite signal to be used.
- the signal selection unit 130 determines whether or not the number of satellite signals reaches a preset number (N 1 ). repeat. During this repetition, the satellite signals other than the satellite signals that have already been determined to be added are processed.
- elevation angle dependency and GNSS type/frequency band dependency are taken into account in normalizing the C/ N0 value.
- the GNSS type means the type of navigation satellite system such as GPS and GLONASS.
- the reason for considering the elevation angle dependence is that the smaller the elevation angle of the satellite, the longer the propagation path in the troposphere near the surface of the earth, and the more the satellite signal tends to be attenuated.
- the reason why the GNSS type dependency is considered is that the signal frequency and transmission power differ depending on the GNSS type, causing a difference in the C/N 0 value.
- the first embodiment and the second embodiment will be described as specific examples of the satellite selection procedure in the initial estimated position determination process.
- the second embodiment mainly the parts different from the first embodiment will be described.
- CN 0max is the maximum C/N 0 value of all received satellite signals in the L1 band.
- dCN 0 is a parameter that determines the selection range of satellite signals. Maximum difference.
- N0 is the number of selected satellite signals, specifically the minimum number of satellite signals selected in the initial position estimation process. In the first embodiment, reception is performed in the L1 band, but reception in the L1 band is an example.
- the signal selection unit 130 normalizes the C/N 0 values of all satellite signals received in the L1 band, considering the GNSS type and elevation angle dependency. Specifically, normalization is performed by adding a GNSS bias value and an elevation angle bias value preset by the bias value setting unit 170 to the C/N 0 value obtained by observation.
- FIG. 7 shows a setting example of the GNSS bias value
- FIG. 8 shows a setting example of the elevation bias value. These bias values are stored in the bias value storage unit 180 .
- the C/N 0 value means the C/N 0 value after normalization.
- the signal selector 130 selects the satellite signal with the largest C/N 0 value from all the satellite signals received in the L1 band, and records the C/N 0 value as CN 0max .
- the signal selector 130 selects the satellite signal with the largest C/N 0 value from all the satellite signals received in the L1 band, and records the C/N 0 value as CN 0max .
- the signal selection unit 130 sets a value smaller than CN 0max by dCN 0 (eg, 10 dB) as the lower limit of C/N 0 with respect to the C/N 0 value (CN 0 ) of the satellite signal selected in S102. Select the satellite signal that satisfies the conditions from the satellite signals that That is, the signal selector 130 selects all satellite signals that satisfy the C/N 0 value CN 0max ⁇ dCN 0 ⁇ C/N 0 ⁇ CN 0max from all the satellite signals received in the L1 band.
- dCN 0 eg, 10 dB
- the signal selection unit 130 determines whether or not the number of satellite signals selected in S102 and S103 is equal to or greater than the preset minimum number of selected satellite signals (N 0 ). If the determination result of S104 is Yes, the signal selection process by the signal selection unit 130 is terminated. The signal selection unit 130 notifies the identification information (code such as the PRN number) of the selected satellite signal to the measurement unit 140, so that the measurement unit 140 performs positioning and time synchronization using the selected satellite signal. It can be carried out.
- the identification information code such as the PRN number
- the signal selection unit 130 selects satellite signals in descending order of the C/N 0 value "CN 0max ⁇ dCN 0 " or less and the next highest C/N 0 value based on the preset GNSS type priority. Signals are selected and supplemented so that the total number of selected satellite signals is N0 .
- FIG. 9 shows an example of setting the priority of GNSS types.
- the setting values of the priority levels of the GNSS types are also stored in the bias value storage unit 180 , and the signal selection unit 130 refers to the setting values stored in the bias value storage unit 180 .
- FIG. 9 shows that GPS has the highest priority and GLO (GLONASS) has the lowest priority.
- GPS and QZSS as navigation satellite systems, are completely synchronized in time and have a small clock bias. Based on such categorization, priorities are set as shown in FIG.
- a premium additional value, unit is dB
- the required number is selected in order from the satellite signal with the highest C/ N0 value.
- the premium of priority 1 is set to 5
- the premium of priority 2 is set to 4
- the premium of priority 3 is set to 3
- the premium of priority 4 is set to 2
- the premium of priority 5 is set to 1.
- N0 is 5 and three satellite signals are selected in S102 and S103.
- Select satellite signal 3 and satellite signal 2 whose values are 29 dB-Hz and 28 dB-Hz.
- the signal receiving unit 120 continuously collects satellite signal data. Concerning the length of time for collection, continuous 24-hour collection is sufficient in an open-sky environment. Other reception environments require longer-term continuous acquisition. Data may be collected at any time to update the bias value.
- the collected satellite signal data are stored in the signal data storage unit 160 as a set of (GNSS type, elevation angle, C/N 0 value).
- the bias value setting unit 170 groups the data of the same GNSS type for each elevation angle range, and calculates the C/ Extract the maximum of the N0 values.
- FIG. 11 shows an example of the processing of S203 in a certain GNSS type.
- the elevation angles are grouped into 0° to 15°, 15° to 30°, 30° to 45°, 45° to 60°, 60° to 75°, and 75° to 90°.
- the maximum C/N 0 value of is extracted.
- the bias value setting unit 170 applies curve fitting to the extracted maximum value data using, for example, the nonlinear least squares method.
- the bias value setting unit 170 repeats curve fitting several times, excluding the largest outlier.
- FIG. 12 shows an example of S204 and S205 for the GNSS types shown in FIG.
- the bias value setting unit 170 generates a fitting function for each GNSS type, and in S207, sets the GNSS type/elevation angle bias value using the fitting function for each GNSS type.
- An example of S206 and S207 is shown in FIG. As shown in FIG. 13, in any GNSS type, a larger bias value is set as the elevation angle decreases. In addition, in the example of FIG. 13, between the GNSS types, bias values are set in the order of GNSS-C>GNSS-B>GNSS-A.
- the dCN 0 value is a parameter that determines the range of C/N 0 values for selecting satellite signals, as described in S103 of FIG.
- the dCN 0 value may be a fixed value that does not depend on the elevation angle of the satellite signal, an example in which the dCN 0 value is determined depending on the elevation angle of the satellite signal will be described below.
- the example described here assumes a case where the reflecting surface of the satellite signal is the vertical wall surface (concrete or glass) of a building, such as in an urban area.
- FIG. 14 shows how a satellite signal with a high elevation angle is incident on a vertical wall surface of a building and is reflected
- FIG. 15 shows how a satellite signal with a low elevation angle is incident on a vertical wall surface of a building and is reflected. showing. As shown in FIGS. 14 and 15, the angle of incidence of the low elevation satellite signals on the vertical walls of the building is greater than the angle of incidence of the high elevation satellite signals on the vertical walls of the building.
- the satellite signal with a low elevation angle has a relatively high reflectance (the signal strength of the reflected wave is large) compared to the satellite signal with a high elevation angle. )It is expected.
- FIG. 16 shows a setting example of the dCN 0 value given elevation angle dependence. As shown in FIG. 16, the dCN 0 value is set to increase as the elevation angle of the satellite signal increases.
- Such elevation angle dependent setting values may be stored in the bias value storage unit 180 in the form of functions corresponding to the curves in FIG . may be stored in the bias value storage unit 180 in the form of a table holding .
- the signal selection unit 130 stores the bias value. Referring to section 180, obtain the dCN 0 value corresponding to the elevation angle of the satellite signal, and use the dCN 0 value to determine whether "CN 0max - dCN 0 ⁇ C/N 0 ⁇ CN 0max " is satisfied. to decide.
- the signal selection unit 130 determines whether the C/N 0 value of a certain satellite signal is equal to or less than "CN 0max -dCN 0 ". Referring to 180, obtain the dCN 0 value corresponding to the elevation angle of the satellite signal and use the dCN 0 value to determine if it is less than or equal to "CN 0max - dCN 0 ".
- the dCN 0 value of a satellite signal with a low elevation angle is smaller than that of a satellite signal with a high elevation angle.
- a satellite signal with a low elevation angle has a narrower range of "CN 0max -dCN 0 ⁇ C/N 0 ⁇ CN 0max " than a satellite signal with a high elevation angle. That is, lower elevation satellite signals are filtered more severely than higher elevation satellite signals. The reason for making the dCN 0 value elevation dependent in this way so that the lower elevation satellite signals are filtered more severely than the higher elevation satellite signals is explained below.
- the satellite signal at a low elevation angle is in a state close to total reflection.
- the difference from the signal intensity (the reference signal intensity normalized by the bias value in FIGS. 7 and 8) when received as a direct wave without an obstacle is small.
- satellite signals at low elevation angles have a longer optical path length in media that attenuate the signal strength, such as the ionosphere and troposphere, and the signal strength when received as a direct wave decreases. Since the drop in signal strength is small, it is necessary to reduce the dCN 0 value and perform strict filtering in order to remove multipath signals (reflected waves) of invisible satellite signals. High elevation satellites do the opposite, widening the range CN 0max -dCN 0 ⁇ C/N 0 ⁇ CN 0max , making them more likely to be selected.
- the dCN 0 value increases as the elevation angle of the satellite signal increases.
- the dCN 0 value may be given a different elevation dependence.
- the technique according to the present invention assumes that at least one visible satellite exists.
- the reference C/N 0 value CN 0max
- the probability that at least one visible satellite exists is improved.
- the second embodiment differs from the first embodiment in that the measuring device 100 selects satellite signals for each frequency band. That is, in the first embodiment, satellite signals were selected only for the L1 band as an example, but in the second embodiment, for each of a plurality of frequency bands output from each satellite, Select a satellite signal.
- the second embodiment is a variation of the embodiment of the invention.
- the reason for selecting satellite signals for each frequency band in the second embodiment is as follows.
- each satellite outputs signals in multiple frequency bands, from the viewpoint of selecting satellites suitable for positioning according to the position of the satellite (visible/invisible), satellite signals in any one frequency band can / If invisibleness can be determined accurately, there is no need to select satellites using signals in each of a plurality of frequency bands.
- a different satellite signal can be used for each frequency band in the positioning calculation. Since the frequency bands supported by satellites differ (for example, the GPS L5 frequency band is only supported by some satellites), by selecting satellite signals individually for each frequency band, variations in positioning calculation policy settings (for each frequency band) can be achieved. , etc. ) can be broadened.
- the device configuration of the measuring device 100 in the second embodiment is the same as the device configuration in the first embodiment, as shown in FIG.
- the operation of each unit is basically the same as in the first embodiment, but differs from the first embodiment in that operations for selecting satellite signals for each frequency band are performed.
- the signal receiving unit 120 sends, for each frequency band, the identification information (code such as the PRN number) of the received satellite signal, the elevation angle of the satellite signal, the C/N 0 value, and the satellite type to the signal selection unit 130. . Further, the signal receiving unit 120 stores identification information, elevation angle, C/N 0 value, and satellite type for each received satellite signal in the signal data storage unit 160 for each frequency band.
- the L1 band and the L2 band are targeted as a plurality of frequency bands.
- the L5 band may be used, or a frequency band other than the L1 band, the L2 band, and the L5 band may be used.
- FIG. 17 is a flow chart showing the operation of the signal selector 130. As shown in FIG. The flow is basically the same as the flow in the first embodiment shown in FIG. 5, but in the second embodiment, the flow in FIG. 17 is repeated for each frequency band, response), it is different from the first embodiment in that it is determined whether or not the minimum C/N 0 value in the frequency band being processed is satisfied. Note that FIG. 17 shows, as an example, processing for the L1 band in repetition for each frequency band.
- CN 0L1 is the lowest C/N 0 value of the selected satellite in the L1 band.
- dCN 0L1 is a parameter that determines the selection range of satellite signals in the L1 band (corresponding to dCN 0 described above).
- N0L1 is the number of selected satellite signals in the L1 band. Similar parameters are set for the L2 band. When using other frequency bands, parameters may be set for each frequency band. For example, if it is the L5 band, CN 0L5 or the like is set.
- the processing of the flow of FIG. 17 is executed for the L1 band.
- the processing of S101 and S102 is the same as in the first embodiment.
- the bias values shown in FIGS. 7 and 8 are set for each frequency band. normalize using .
- the satellite signal with the largest C/N 0 value is selected from all the satellite signals received in the corresponding frequency band (first L1 band), and the selected C/N 0 value is recorded as CN 0max .
- the signal selection unit 130 sets a value smaller than CN 0max by dCN 0L1 (eg, 10 dB) as the lower limit of C/N 0 with respect to the C/N 0 value (CN 0max ) of the satellite signal selected in S102. Select the satellite signal that satisfies the conditions from the satellite signals that Here, the signal selection unit 130 determines that the C/N 0 value satisfies CN 0max ⁇ dCN 0L1 ⁇ C/N 0 ⁇ CN 0max and CN 0L1 ⁇ C/N 0 from all satellite signals received in the L1 band. Select all satellite signals that satisfy
- the signal selection unit 130 determines whether or not the number of satellite signals selected in S102 and S103 is equal to or greater than the preset minimum number of selected satellite signals (N 0L1 ). If the determination result of S104 is Yes, the signal selection processing by the signal selection unit 130 for the L1 band is ended. The signal selection unit 130 notifies the identification information (code such as the PRN number) of the selected satellite signal to the measurement unit 140, so that the measurement unit 140 performs positioning and time synchronization using the selected satellite signal. It can be carried out.
- the identification information code such as the PRN number
- the signal selection unit 130 executes the processing of the flow in FIG. 17 for the next frequency band (for example, the L2 band).
- the processing of S105 may be the same as the processing described in the first embodiment. That is, in S105, the signal selection unit 130, based on the preset GNSS type priority (example: FIG. 9), the C / N 0 value is "CN 0max - dCN 0L1 " or less C / N 0 value is large Satellite signals are selected in order from the next satellite signal, and compensation is performed so that the total number of selected satellite signals becomes N0L1 .
- the setting of priorities as shown in FIG. 9 may be determined for each frequency band.
- the signal selector 130 selects a satellite signal using the priority corresponding to the frequency band being processed.
- a premium additional value
- the premium of priority 1 is set to 5
- the premium of priority 2 is set to 4
- the premium of priority 3 is set to 3
- the premium of priority 4 is set to 2
- the premium of priority 5 is set to 1.
- N 0L1 is 5 and three satellite signals are selected in S102 and S103.
- Select satellite signal 3 and satellite signal 2 whose values are 29 dB-Hz and 28 dB-Hz.
- the signal selection unit 130 When S105 ends, the signal selection unit 130 notifies the identification information (code such as the PRN number) of the selected satellite signal to the measurement unit 140, so that the measurement unit 140 can use the selected satellite signal. Positioning and time synchronization can be performed.
- the signal selection unit 130 executes the processing of the flow in FIG. 17 for the next frequency band (for example, the L2 band).
- positioning and time synchronization are performed using satellite signals selected for each frequency band. may be used for positioning and time synchronization.
- satellite signals 1, 2, 3, and 4 are selected in the L1 band, and satellite signals 5, 6, 7, and 8 are selected in the L2 band.
- the DOP (Dilution of Precision) values of satellite signals 1, 2, 3, and 4 are compared with the DOP values of satellite signals 5, 6, 7, and 8, and the frequency band with the smaller DOP value is selected. may be selected to perform positioning and time synchronization using L1 band and L2 band signals.
- the signal selector 130 may select a backup satellite signal by considering the DOP value. For example, the signal selection unit 130 selects the satellite signals A, B, and C in descending order of the C/N 0 value of "CN 0max ⁇ dCN 0 " or less and the next highest C/N 0 value as the substitute satellite signal candidate. , calculate the DOP value when the candidate satellites A, B, and C are added to the already selected satellite signal, and select the satellite signal with the smallest DOP value.
- satellite signals already selected at S104 are satellite signals 1, 2, and 3, "satellite signals 1, 2, 3, A”, “satellite signals 1, 2, 3, B", Calculate the DOP value for each of "satellite signals 1, 2, 3, C”. If "satellite signal 1, 2, 3, A” has the lowest DOP value, select "satellite signal 1, 2, 3, A”. If the number of satellites to be selected is greater than 4, the above process should be repeated until that number of satellites is reached.
- selection method 1 the alternate satellite signal selection method described in the first embodiment
- selection method 2 the above method using the DOP value
- the signal selector 130 may select the sub-satellite signal by a combination of selection methods 1 and 2.
- selection method 1 is used to select an alternate satellite signal
- selection method 2 is performed on each of the selected satellite signals, and a satellite signal with a small DOP value is selected.
- a cost value is set as an expected value of the degree of improvement in positioning accuracy by satellite signals selected based on selection methods 1 and 2, and the total cost value (evaluation value) of selection methods 1 and 2 is set.
- satellite signal A and satellite signal B are selected as substitute satellite signal candidates by a combination of selection methods 1 and 2 .
- the C/N 0 value of satellite signal A is 30 dB-Hz
- the C/N 0 value of satellite signal B is 28 dB-Hz
- the DOP value when satellite signal A is selected is 5, and satellite signal B is selected.
- the DOP value is 4 when the total cost value (evaluation value) of selection methods 1 and 2 is set as "DOP value/(C/N 0 value)", the cost values of satellite signals A and B are 1/6 and 1/7, respectively. , satellite signal B is smaller than satellite signal A, so satellite signal B is selected as the alternate satellite signal.
- the bias value setting operation executed by the bias value setting unit 170 is basically the same as the bias value setting operation in the first embodiment, but the second embodiment This differs from the first embodiment in that a bias value is set for each frequency band.
- FIG. 19 shows a flowchart of the bias value setting operation in the second embodiment.
- the signal receiver 120 continuously collects satellite signal data for each frequency band.
- the collected satellite signal data is stored in the signal data storage unit 160 as a set of (GNSS type, frequency band, elevation angle, C/N 0 value).
- the bias value setting unit 170 divides the data of the same GNSS type and frequency band into groups for each elevation angle range. Extract the maximum value of the C/N 0 value of An example of the processing is as described with reference to FIG.
- the bias value setting unit 170 applies curve fitting to the extracted maximum value data using a nonlinear least squares method or the like.
- the bias value setting unit 170 repeats curve fitting several times, excluding the largest outlier. Examples of S204 and S205 for the GNSS types shown in FIG. 11 are as shown in FIG.
- the bias value setting unit 170 In S206, the bias value setting unit 170 generates a fitting function for each GNSS type, and in S207, sets the GNSS type/elevation angle bias value using the fitting function for each GNSS type.
- FIG. 20 shows an example of S206 and S207 for the L1 band.
- a bias value is set for each frequency band of each satellite signal (for example, L1 band, L2 band, and L5 band in the case of GPS) for each GNSS satellite type. This is because the satellite signal reception characteristics depend not only on the GNSS satellite type and elevation angle, but also on the frequency band of the satellite signal. A bias value based only on the GNSS satellite type and frequency band may be set without using the elevation angle.
- FIGS. 21 and 22 show actual measurement examples of differences in reception characteristics depending on frequency bands for the same combination of GNSS antenna and GNSS receiver.
- FIG. 21 shows the GPS L1 signal
- FIG. 22 shows the GPS L2 signal.
- the horizontal axis is the elevation angle (°)
- the vertical axis is the C/N 0 value (dB-Hz).
- the transmission signal output may differ depending on individual satellites.
- the transmission signal output may differ depending on the orbit (GEO/IGSO/MEO) of the satellite.
- bias values may be set for individual satellites.
- the satellite signal is corrected by adding an individual bias value to the reception quality during normalization.
- the separate bias values apply in addition to the GNSS and elevation bias values shown in FIGS.
- only the separate bias values may be applied without applying the GNSS bias value and the elevation bias value.
- elevation bias values and individual bias values may be applied without applying GNSS bias values.
- the received signal strength is measured in advance for each GNSS type and for each satellite, and the measured value is stored in the signal data storage unit 160.
- the bias value setting unit 170 selects a satellite in which an event similar to that of satellite A described above occurs, and sets an individual bias value for the selected satellite.
- an individual bias value for a specific satellite may be applied to cases other than the case related to the transmission signal output such as satellite A described above.
- dCN 0 (Regarding the setting value of dCN 0 ) Regarding the dCN 0 value for each frequency band (that is, dCN 0L1 value, dCN 0L2 value, etc.) in the second embodiment, as described with reference to FIGS.
- a value of 0 may be set.
- a dCN 0 value given elevation angle dependence as shown in FIG. 16 is set for each frequency band.
- FIG. 23 is a diagram showing a hardware configuration example of a computer that can be used as the measuring device 100 in this embodiment.
- the computer may be a computer as a physical device, or may be a virtual machine on the cloud.
- the computer of FIG. 23 has a drive device 1000, an auxiliary storage device 1002, a memory device 1003, a CPU 1004, an interface device 1005, a display device 1006, an input device 1007, an output device 1008, etc., which are connected to each other via a bus B. .
- a bus B a bus connecting a drive device 1000, an auxiliary storage device 1002, a memory device 1003, a CPU 1004, an interface device 1005, a display device 1006, an input device 1007, an output device 1008, etc., which are connected to each other via a bus B.
- the GNSS antenna 110 is not shown in FIG. GNSS antenna 110 is, for example, connected to interface device 1005 .
- a program that implements the processing in the computer is provided by a recording medium 1001 such as a CD-ROM or memory card, for example.
- a recording medium 1001 such as a CD-ROM or memory card
- the program is installed from the recording medium 1001 to the auxiliary storage device 1002 via the drive device 1000 .
- the program does not necessarily need to be installed from the recording medium 1001, and may be downloaded from another computer via the network.
- the auxiliary storage device 1002 stores installed programs, as well as necessary files and data.
- the memory device 1003 reads and stores the program from the auxiliary storage device 1002 when a program activation instruction is received.
- the CPU 1004 implements functions related to the measuring device 100 according to programs stored in the memory device 1003 .
- Interface device 1005 is used as an interface to connect to GNSS antenna 110 .
- a display device 1006 displays a GUI (Graphical User Interface) or the like by a program.
- An input device 1007 is composed of a keyboard, a mouse, buttons, a touch panel, or the like, and is used to input various operational instructions.
- the output device 1008 outputs the calculation result.
- visible satellite signals can be selected from received satellite signals with high accuracy without depending on the characteristics of the receiver and antenna. This makes it possible to improve positioning accuracy in poor reception environments. Also, the influence of interfering signals on satellite selection can be reduced.
- At least a measuring device, a measuring method, and a program described in the following items are provided.
- (Section 1) a signal selection unit that selects a predetermined number of satellite signals based on the reception quality of the satellite signals received by the GNSS antenna; a measurement unit that performs positioning using a predetermined number of satellite signals selected by the signal selection unit and determines an initial estimated position; The signal selection unit is used for positioning or time synchronization based on a comparison result of comparing the reception time expected when the satellite signal is received as a direct wave at the initial estimated position and the actual reception time of the satellite signal.
- a satellite signal receiver that selects a plurality of satellite signals to receive.
- the signal selection unit rejects the satellite signal actually received at a delay time equal to or greater than a threshold from the expected reception time when received as a direct wave at the initial estimated position based on the comparison result. Satellite signal receiver as described. (Section 3) 3. The satellite signal receiving device according to claim 2, wherein the measurement unit performs positioning using a plurality of satellite signals excluding the rejected satellite signals, and updates the initial estimated position with the obtained coordinate values. (Section 4) The signal selection unit updates the initial estimated position, repeats the process based on the comparison result using the updated initial estimated position, and terminates the repeated process when a predetermined termination condition is satisfied. Satellite signal receiving device according to .
- the signal selector compares the delay times among the plurality of non-selected satellite signals, and delays the signals. Satellite signal reception according to claim 4, wherein the non-selected satellite signals to be added to the selected satellite signal are determined based on the DOP value when the non-selected satellite signals with shorter times are added to the selected satellite signal.
- a satellite signal selection method performed by a satellite signal receiving device, comprising: selecting a predetermined number of satellite signals based on the reception quality of the satellite signals received by the GNSS antenna; performing positioning using the predetermined number of satellite signals to determine an initial position estimate; Select multiple satellite signals to be used for positioning or time synchronization based on the result of comparison between the expected reception time when the satellite signal is received as a direct wave at the initial estimated position and the actual reception time of the satellite signal.
- a satellite signal selection method comprising the step of selecting.
- (Section 7) A program for causing a computer to function as each unit in the satellite signal receiving apparatus according to any one of items 1 to 5.
- bias value setting unit 180 bias value storage unit 1000 drive device 1001 recording medium 1002 auxiliary storage device 1003 memory device 1004 CPU 1005 interface device 1006 display device 1007 input device 1008 output device
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
GNSSアンテナにより受信した衛星信号の受信品質に基づいて、所定数の衛星信号を選択する信号選択部と、前記信号選択部により選択された所定数の衛星信号を用いて測位を実行し、初期推定位置を決定する計測部と、を備え、前記信号選択部は、衛星信号を初期推定位置において直接波として受信した場合に期待される受信時刻と当該衛星信号の実際の受信時刻とを比較した比較結果に基づいて、測位又は時刻同期に使用する複数の衛星信号を選択する衛星信号受信装置。
Description
本発明は、GNSS(Global Navigation Satellite System:全地球航法衛星システム)による測位及び時刻同期を高精度に行う技術に関連するものである。
近年、GNSSによる測位及び時刻同期が幅広いアプリケーションにおいて活用されている。
GNSSによる測位及び時刻同期においては、GNSSアンテナにより受信したGNSS衛星信号(以降、衛星信号)を用いて測位及び時刻同期の処理が実行される。
GNSSアンテナの設置位置の周囲に存在する構造物等により衛星信号の見通し状態での受信が遮られる場合がある。その場合、当該衛星信号は、GNSSアンテナにおいて必要な信号強度で受信されないか、又は、GNSSアンテナの設置位置の周囲に存在する構造物等により反射・回折するマルチパスにより、不可視衛星信号として受信されることになる。その結果、GNSSによる測位性能及び時刻同期性能が劣化する。
塚越他"オープンソースGNSSライブラリを用いた遊歩道環境での自律移動ロボットナビゲーションのための位置推定"計測自動制御学会論文集Vol.52, No.5, 276/283(2016)
GNSSによる測位及び時刻同期の精度を向上させる上では見通し状態で受信可能な可視衛星信号を多く受信し、精度の劣化に大きな影響を及ぼす、見通し状態で受信できない不可視衛星信号を測位及び時刻同期で使用する衛星信号から有効に排除することが重要である。
本発明は上記の点に鑑みてなされたものであり、衛星信号の受信環境が良くない場合でも、衛星信号を適切に選択し、精度良く、GNSSによる測位及び時刻同期を行うことを可能とする技術を提供することを目的とする。
開示の技術によれば、GNSSアンテナにより受信した衛星信号の受信品質に基づいて、所定数の衛星信号を選択する信号選択部と、
前記信号選択部により選択された所定数の衛星信号を用いて測位を実行し、初期推定位置を決定する計測部と、を備え、
前記信号選択部は、衛星信号を初期推定位置において直接波として受信した場合に期待される受信時刻と当該衛星信号の実際の受信時刻とを比較した比較結果に基づいて、測位又は時刻同期に使用する複数の衛星信号を選択する
衛星信号受信装置が提供される。
前記信号選択部により選択された所定数の衛星信号を用いて測位を実行し、初期推定位置を決定する計測部と、を備え、
前記信号選択部は、衛星信号を初期推定位置において直接波として受信した場合に期待される受信時刻と当該衛星信号の実際の受信時刻とを比較した比較結果に基づいて、測位又は時刻同期に使用する複数の衛星信号を選択する
衛星信号受信装置が提供される。
開示の技術によれば、衛星信号の受信環境が良くない場合でも、精度良く、GNSSによる測位及び時刻同期を行うことを可能とする技術が提供される。
以下、図面を参照して本発明の実施の形態(本実施の形態)を説明する。以下で説明する実施の形態は一例に過ぎず、本発明が適用される実施の形態は、以下の実施の形態に限定されるわけではない。
(課題の詳細、実施の形態の概要)
近年、GPS以外の航法衛星システムとしてGLONASS、Galileo、BeiDou、QZSS、NAVICなどが利用できるようになり、衛星数が増加している。
近年、GPS以外の航法衛星システムとしてGLONASS、Galileo、BeiDou、QZSS、NAVICなどが利用できるようになり、衛星数が増加している。
前述したように、GNSSによる測位及び時刻同期の精度を向上させる上では見通し状態で受信可能な可視衛星信号を多く受信し、精度の劣化に大きな影響を及ぼす、見通し状態で受信できない不可視衛星信号を測位及び時刻同期で使用する衛星信号から有効に排除することが重要である。
不可視衛星信号を排除する従来の方法として、受信した衛星信号から予め設定した閾値以下のC/N0(Carrier-To-Noise Power Density:搬送波対雑音密度比)値の衛星信号を排除する、C/N0マスク方式が知られている。
しかし、衛星信号のC/N0値はアンテナのゲイン、受信機の受信感度、アンテナ‐受信機間のケーブル損失、衛星種別等に依存するため最適な閾値を設定することは困難である。
また、C/N0マスク方式では衛星信号の信号帯域幅に干渉信号が混入した場合、衛星信号のC/N0値が全体的に低下し、C/N0マスクにより衛星信号をロストする結果、測位及び時刻同期ができなくなるリスクがある。このような干渉信号として、意図的に発生するGNSSジャミング(妨害)信号の他、機器が発生するノイズ、他の通信システムからの干渉信号がある。
本実施の形態では、見通し状態で受信可能な可視衛星信号を多く受信し、精度の劣化に大きな影響を及ぼす、見通し状態で受信できない不可視衛星信号を測位及び時刻同期で使用する衛星信号から有効に排除するために、可視衛星数が多く確保可能なマルチGNSS環境を想定した後述する手順で、測位及び時刻同期の使用に適した衛星信号を選択することとしている。当該手順は、受信した衛星信号の強度(受信品質)を選択の根拠に使用するもののアンテナ及び受信機の個体特性、干渉信号の影響を排除した衛星信号の選択を可能としている。
また、当該手順では可視衛星信号を優先的に選択し、可視衛星数が少ない場合は、可視衛星に加えて精度劣化への影響の小さい、伝搬遅延の小さい不可視衛星を選択する。
以下、本発明の実施の形態における構成及び動作の例を詳細に説明する。なお、以下で説明する処理では、受信品質の指標としてC/N0値を使用するが、C/N0値以外の受信品質の指標を使用してもよい。また、本実施の形態において、"衛星信号を選択する"際の"衛星信号"は、その衛星信号の送信元のGNSS衛星と紐づいているものとする。例えば、GNSS衛星A、GNSS衛星B、GNSS衛星Cを3機の異なる任意のGNSS衛星であるとすると、3個の衛星信号を選択するとは、GNSS衛星Aからの衛星信号、GNSS衛星Bからの衛星信号、及びGNSS衛星Cからの衛星信号を選択することを意味する。
(装置構成)
図1に、本発明の実施の形態における計測装置100の構成例を示す。本実施の形態における計測装置100は、GNSSアンテナ110、信号受信部120、信号選択部130、計測部140、出力部150、信号データ格納部160、バイアス値設定部170、バイアス値格納部180を有する。なお、計測装置100は、衛星信号を受信して処理する装置であり、これを「衛星信号受信装置」と呼んでもよい。
図1に、本発明の実施の形態における計測装置100の構成例を示す。本実施の形態における計測装置100は、GNSSアンテナ110、信号受信部120、信号選択部130、計測部140、出力部150、信号データ格納部160、バイアス値設定部170、バイアス値格納部180を有する。なお、計測装置100は、衛星信号を受信して処理する装置であり、これを「衛星信号受信装置」と呼んでもよい。
GNSSアンテナ110は、軌道上のGNSS衛星から送信される電波を受信し、電波を電気信号に変換する。この電気信号を「衛星信号」と呼んでもよい。
GNSSアンテナ110と信号受信部120とはケーブルで接続され、衛星信号はケーブルにより信号受信部120に送られる。GNSSアンテナ110と信号受信部120との間の距離が長い場合には、GNSSアンテナ110と信号受信部120との間に増幅器が備えられてもよい。
信号受信部120は、衛星信号を受信し、C/N0値を計測するとともに、受信した衛星信号の送信元のGNSS衛星の種別を識別する。また、衛星の軌道情報(例:アルマナック、エフェメリス)を使用し、仰角を計測する。衛星の軌道情報は衛星信号の航法メッセージから取得してもよいし、他の手段(例:ネットワーク上のサーバ)から取得してもよい。信号受信部120は、受信した衛星信号の識別情報(PRN番号等のコード)、当該衛星信号の仰角、C/N0値、及び衛星種別を信号選択部130に送る。また、信号受信部120は、受信した衛星信号毎の識別情報、仰角、C/N0値、衛星種別を信号データ格納部160に格納する。なお、仰角とは、衛星信号の受信地点(すなわちGNSSアンテナ)から衛星信号の送信元のGNSS衛星を見る場合の視線と水平面とのなす角度である。例えば、GNSS衛星が天頂にある場合、その仰角は90°である。
本実施の形態において対象とするGNSS衛星の種別は、GPS、GLONASS、Galileo、BeiDou、QZSS、及びNAVICである。ただし、これらは例であり、これらの種別よりも多くてもよいし、少なくてもよい。
信号選択部130は、受信した複数の衛星信号の中から、測位と時刻同期に使用する衛星信号を選択する。選択の手順については後述する。
計測部140は、絶対時刻に対して精密に時刻が管理される原子時計を搭載したGNSS衛星から送信される衛星信号を用いて時刻同期を行うことにより、絶対時刻に対して高精度に時刻同期した時刻情報を算出する。ここでの絶対時刻とは、例えば、協定世界時(UTC:Coordinated Universal Time)である。なお、計測部140は、測位と時刻同期のうちの一方のみを行うこととしてもよい。
受信した衛星信号から、当該衛星信号がGNSS衛星から送信された絶対時刻を知ることができるが、GNSS衛星からGNSSアンテナ110の位置に衛星信号が到達するまでの伝搬時間を計測し、計測部140の時刻と衛星の時刻との間の時刻オフセット値Δtを補正しなければ、受信位置において正確な絶対時刻は得られない。
そこで、計測部140は、例えば4機以上のGNSS衛星からの衛星信号を用いて、受信位置の3次元座標情報(x,y,z)、及び時刻オフセット(Δt)の4つのパラメータをコード測位により算出することにより、測位と時刻同期を同時に行う。計測部140はコード測位の他に搬送波位相測位を行ってもよい。
計測部140は、出力部150を介してこの絶対時刻に基づく時刻情報と測位結果である位置情報を出力する。例えば、計測装置100がモバイルネットワークにおける基地局であるとすると、当該基地局が絶対時刻に同期した時刻情報を利用して、例えば、(絶対時刻に同期している)隣接する基地局とTDD(Time Division Duplex:時分割複信)信号フレームの上り、下り信号のタイムスロット構成(並び)を一致させた上で、信号フレームの送信タイミングを同期させることで、隣接する基地局と互いに干渉しないようにTDD信号を送信することができる。
バイアス値設定部170は、信号データ格納部160に格納されている衛星信号データを用いて、バイアス値を設定(算出)し、設定したバイアス値をバイアス値格納部180に格納する。バイアス値格納部180に格納されたバイアス値は、信号選択部130における衛星信号の選択処理に用いられる。バイアス値設定部170によるバイアス値設定動作の詳細については後述する。
本実施の形態における計測装置100は、物理的にまとまった1つの装置であってもよいし、いくつかの機能部が物理的に分離していて、分離された複数の機能部がネットワークにより接続された装置であってもよい。
また、計測装置100は、図1に示す機能を全て含むこととしてもよいし、一部の機能(例えば、信号選択部130と計測部140)がネットワーク上(例えばクラウド上)に備えられ、残りの機能が計測装置100に搭載されて使用されてもよい。
例えば、計測装置100に備えられた信号受信部120から観測データを出力し、当該観測データをクラウド上に設けた「信号選択部130と計測部140」からなる装置に送信することで、衛星信号選択及び測位演算をクラウド上で実施してもよい。この場合、クラウド上の計測部140から、出力部150へ測位演算結果が返される。
また、計測装置100における「信号データ格納部160とバイアス値設定部170」からなる装置がネットワーク上(例えばクラウド上)に備えられ、残りの機能が計測装置100に搭載されて使用されてもよい。
例えば、計測装置100に備えられた信号受信部120から観測データを出力し、当該観測データをクラウド上に設けた信号データ格納部160に格納し、クラウド上に設けたバイアス値設定部170が、格納したデータを用いてバイアス値の設定を行う。この場合、クラウド上のバイアス値設定部170から、バイアス値格納部180へバイアス値が返される。
(衛星信号の選択のための処理内容)
次に、衛星信号の選択のための計測装置100の動作(特に信号選択部130及び計測部140の動作)について説明する。計測装置100は、(1)初期推定位置決定処理、(2)衛星選択処理、(3)補欠衛星選択処理の手順で衛星信号の選択処理を行う。各処理の概要は下記のとおりである。
次に、衛星信号の選択のための計測装置100の動作(特に信号選択部130及び計測部140の動作)について説明する。計測装置100は、(1)初期推定位置決定処理、(2)衛星選択処理、(3)補欠衛星選択処理の手順で衛星信号の選択処理を行う。各処理の概要は下記のとおりである。
(1)初期推定位置決定処理は、衛星選択を行うにあたり基準となる初期推定位置を決定する処理である。
(2)衛星選択処理は、初期推定位置における推定受信時刻と実際の受信時刻との比較によりLOS(Line Of Sight)衛星信号に加え、遅延の小さいNLOS(Non Line Of Sight)衛星信号を選択する処理である。
(3)補欠衛星選択処理は、選択された衛星信号数が少ない場合、DOP(Dilution Of Precision)値の劣化に伴う測位及び時刻同期精度の劣化を抑制するために追加で衛星を選択する処理である。
以下、(1)初期推定位置決定処理、(2)衛星選択処理、(3)補欠衛星選択処理のそれぞれを、フローチャートを参照して説明する。
ここで使用するパラメータの記号の意味は、下記のとおりである。なお、初期推定位置決定処理における衛星信号の選択処理の具体例を後述するが、その具体例で使用するパラメータについては、その具体例において改めて説明する。
dCN0:受信した衛星信号のC/N0値の最大値と初期推定位置処理で選択する衛星信号のC/N0値の差分の最大値
N0:初期推定位置決定処理において選択される最低衛星信号数
dT:衛星選択処理における初期推定位置において期待される受信時刻からの遅延時間の閾値
Ns:衛星選択処理において初期推定位置が更新される最大回数
dp:衛星選択処理において完了判定をするための更新前後の初期推定位置の差分の閾値
N1:補欠衛星選択プロセスにおいて選択される最低衛星信号数
<(1)初期推定位置決定処理>
図2のフローチャートを参照して、初期推定位置決定処理を説明する。
N0:初期推定位置決定処理において選択される最低衛星信号数
dT:衛星選択処理における初期推定位置において期待される受信時刻からの遅延時間の閾値
Ns:衛星選択処理において初期推定位置が更新される最大回数
dp:衛星選択処理において完了判定をするための更新前後の初期推定位置の差分の閾値
N1:補欠衛星選択プロセスにおいて選択される最低衛星信号数
<(1)初期推定位置決定処理>
図2のフローチャートを参照して、初期推定位置決定処理を説明する。
S1-1において、信号選択部130は、信号受信部120により受信した複数の衛星信号から最もC/N0値の大きい衛星信号を抽出する。
S1-2において、信号選択部130は、S1-1で抽出された最大のC/N0値とのC/N0値の差分が予め設定された値(dCN0)以下である衛星信号を、上記複数の衛星信号から選択する。
S1-3において、信号選択部130は、S1-1とS1-2で選択された衛星信号の数が予め設定された最低衛星信号数(N0)に満たない場合はC/N0値の大きい(次点の)衛星信号から、最低衛星信号数(N0)になるまで衛星信号を補完する。
S1-4において、計測部140は、S1-1~S1-3で選択された衛星信号を使用したコード測位により、計測装置100の初期推定位置を決定する。なお、測位方法として、コード測位以外の測位方法を用いてもよい。以降の処理のコード測位を行うステップにおいても同様である。
なお、S1-1~S1-3については、より詳細な例を第1の実施の形態及び第2の実施の形態として後述する。
<(2)衛星選択処理>
次に、図3のフローチャートを参照して、衛星選択処理を説明する。
次に、図3のフローチャートを参照して、衛星選択処理を説明する。
S2-1において、信号選択部130は、S1-4のコード測位の過程で受信機(計測装置100)のクロックバイアス、衛星システム間の時刻バイアスを補正した前提で、各衛星信号を初期推定位置において信号受信部120により直接波として受信した場合に期待される受信時刻と実際の受信時刻を比較する。
衛星信号を直接波として受信した場合に期待される受信時刻については、軌道情報から算出される当該衛星信号の送信元の衛星位置と初期推定位置との間の直線距離から算出することができる。実際の受信時刻は、信号受信部120、信号選択部130、計測部140のいずれで測定してもよい。
S2-2において、信号選択部130は、S2-1の比較により、実際の受信時刻が期待される初期推定位置での受信時刻よりもdT以上の時間遅れて受信された衛星信号を棄却する。
例えば、dTをX[ns]とする。衛星信号Aの期待される受信時刻をTとし、衛星信号Aの実際の受信時刻がT+X+1[ns]であったとすると、遅延時間がX+1(dT以上)なので、衛星信号Aは棄却される。
S2-3において、計測部140は、S2-2で棄却された衛星信号を除く、信号受信部120により受信した複数の衛星信号でコード測位を行い、得られた座標値で計測装置100の初期推定位置を更新する。
S2-4において、(1)S2-3で選択された衛星信号が4個未満の場合(言い換えると、選択された衛星が4機以下の場合)、(2)更新前後の初期推定位置のベクトル差が予め設定した値(dp)以下の場合、(3)予め設定した計算回数(Ns)に達した場合、のいずれかに該当する場合には計算を終了する。(1)~(3)のいずれでもない場合には、S2-1~S2-3を繰り返す。なお、計算を終了する条件は(1)~(3)に限定されない。
上記の(2)について、例えば、n回目の初期推定位置の更新後の位置がxn(位置を表すベクトル)であり、n+1回目の初期推定位置の更新後の位置がxn+1である場合において、|xn-xn+1|がdp[m]以下であれば、計算を終了する。
S2-1~S2-3の繰り返し処理において、繰り返し毎に、dT値を小さくしてもよい。例えば、1回目のS2-2の計算において、dT=X[ns]とし、n(nは2以上の整数))回目のS2-2の計算において、dT=X-(n-1)×Δ(Δは正の実数、但しdT>0)とすることとしてもよい。このような処理により、繰り返し毎にdTが小さくなるので、棄却される衛星信号が増加し、より期待される受信時刻からの遅延の小さい衛星信号が選択され易くなる。
<(3)補欠衛星選択処理>
次に、図4のフローチャートを参照して、補欠衛星選択処理を説明する。
次に、図4のフローチャートを参照して、補欠衛星選択処理を説明する。
S3-1において、信号選択部130は、S2-4の計算終了の時点で選択された衛星信号の数が予め設定された数(N1)よりも小さい場合にはこの数に達するために補う補欠の衛星信号を選択するため、選択候補となる複数の非選択衛星信号の間で遅延時間を比較する。遅延時間とは、S2-2でdTと比較の対象となる「衛星信号を初期推定位置において直接波として受信した場合に期待される受信時刻と実際の受信時刻との差分」である。各非選択衛星信号の遅延時間については、繰り返しの最後のS2-2で測定・算出した値を使用してもよいし、S2-4の計算終了の時点における初期推定位置を用いてこのS3-1の時点で新たに測定・算出した値を使用してもよい。
S3-2において、信号選択部130は、S3-1で遅延時間の最も小さい衛星信号とその次に遅延時間の小さい衛星信号に対し、S2-4の計算終了時点で選択された衛星信号に加えた場合のDOP値を比較し、遅延時間とDOP値で計算されるコスト値が小さくなる方の衛星信号を、追加する衛星信号として選択する。
例えば、S2-4の計算終了時点で衛星信号A、衛星信号B、衛星信号Cが選択されたとする。S3-2において、遅延時間の最も小さい衛星信号として衛星信号Dが選択され、その次に遅延時間の小さい衛星信号として衛星信号Eが選択されたとする。
ここで、衛星信号Dの遅延時間が2、衛星信号Eの遅延時間が3、「衛星信号A、衛星信号B、衛星信号C、衛星信号D」のDOP値が7、「衛星信号A、衛星信号B、衛星信号C、衛星信号E」のDOP値が4になったとする。コスト値を「遅延時間×DOP値」であるとすると、「衛星信号A、衛星信号B、衛星信号C、衛星信号D」のコスト値は、2×7=14になり、「衛星信号A、衛星信号B、衛星信号C、衛星信号E」のコスト値は、3×4=12になるので、この場合、衛星信号Eが追加する衛星信号として選択される。
また、コスト値に衛星種別のプレミアムを設定してもよい。例えば、後述する図9に示すように衛星の優先順位が設定されている場合において、優先順位が高いほど(優先順位の数値が小さいほど)、コスト値のプレミアムとして小さな値を付与する。
例えば、図9の例において、優先順位1のプレミアムを1とし、優先順位2のプレミアムを2とし、優先順位3のプレミアムを3とし、優先順位4のプレミアムを4とし、優先順位5のプレミアムを5とする。
上記の衛星信号D、Eの例において、衛星信号Dのプレミアムが1、衛星信号Eのプレミアムが4であるとすると、「衛星信号A、衛星信号B、衛星信号C、衛星信号D」のコスト値は、2×7+1=15になり、「衛星信号A、衛星信号B、衛星信号C、衛星信号E」のコスト値は、3×4+4=16になるので、この場合、衛星信号Dが追加する衛星信号として選択される。
S3-3において、信号選択部130は、衛星信号の数が予め設定された数(N1)に達するか否かを判定し、達していない場合にはS3-1、S3-2のプロセスを繰り返す。この繰り返しの際には、既に追加を決定した衛星信号を除く衛星信号に対して処理を行う。
S3-3において、衛星信号の数が予め設定された数(N1)に達した場合、計算を終了し、例えば、選択された衛星信号を用いた測位及び時刻同期を行う。
(初期推定位置決定処理における処理の具体例)
以下、計測装置100の具体的な動作例として、前述した初期推定位置決定処理における衛星信号の選択の手順の具体例を説明する。
以下、計測装置100の具体的な動作例として、前述した初期推定位置決定処理における衛星信号の選択の手順の具体例を説明する。
以下の具体例の説明において、C/N0値の正規化にあたり、仰角依存性とGNSS種別・周波数帯依存性を考慮している。ここでGNSS種別とはGPSやGLONASS等の航法衛星システムの種別を意味する。仰角依存性を考慮する理由は衛星の仰角が小さいほど地表に近い対流圏における伝搬経路が長くなり、衛星信号がより減衰する傾向にあるためである。GNSS種別依存性を考慮する理由はGNSS種別により信号周波数や送信電力が異なり、C/N0値に差分を生じるためである。さらに同一のGNSS種別でも信号の周波数帯(GPSの場合はL1帯、L2帯、L5帯、等)によりC/N0値に差分を生じる。なお、仰角依存性とGNSS種別・周波数帯依存性のうちのいずれか一方を考慮することとしてもよい。
以下では、初期推定位置決定処理における衛星選択の手順の具体例として、第1の実施の形態と、第2の実施の形態を説明する。第2の実施の形態については、第1の実施の形態と異なる部分について主に説明する。
―――――――――――第1の実施の形態―――――――――――
(第1の実施の形態の信号選択部130の動作例)
第1の実施の形態における信号選択部130の動作例を、図5に示すフローチャートの手順に沿って詳細に説明する。手順の説明の中で、図6~図9も参照する。
第1の実施の形態における信号選択部130の動作例を、図5に示すフローチャートの手順に沿って詳細に説明する。手順の説明の中で、図6~図9も参照する。
まず、図6を参照して、手順において使用される設定パラメータについて説明する。図6に示すとおり、CN0maxは、L1帯の全受信衛星信号のC/N0値の最大値である。dCN0は、衛星信号の選択範囲を決めるパラメータであり、具体的には、受信した衛星信号のC/N0値の最大値と初期推定位置処理で選択する衛星信号のC/N0値の差分の最大値である。N0は、選択衛星信号数であり、具体的には、初期推定位置決定処理において選択される最低衛星信号数である。なお、第1の実施の形態では、L1帯で受信を行っているが、L1帯で受信を行うことは一例である。
図5のS101において、信号選択部130は、L1帯で受信した全衛星信号のC/N0値をGNSS種別及び仰角依存性を考慮して正規化する。具体的には、観測で得られたC/N0値に、バイアス値設定部170により予め設定されたGNSSバイアス値と仰角バイアス値を加えることにより正規化を行う。
図7にGNSSバイアス値の設定例を示し、図8に仰角バイアス値の設定例を示す。これらのバイアス値はバイアス値格納部180に格納されている。
例えば、ある衛星信号の観測で得られたC/N0値が30dB-Hzであり、仰角が30°であり、衛星種別がGLO(GLONASS)であるとすると、信号選択部130は、当該衛星信号の補正後(正規化後)のC/N0値を30+4+2=36dB-Hzとする。以降、C/N0値は正規化後のC/N0値であることを意味する。
図5のS102において、信号選択部130は、L1帯で受信した全衛星信号の中からC/N0値の最も大きい衛星信号を選択し、そのC/N0値をCN0maxとして記録する。なお、ここでは、前提条件として可視衛星信号が少なくとも1つ存在することを想定している。
S103において、信号選択部130は、S102で選択した衛星信号のC/N0値(CN0)に対し、CN0maxよりdCN0(例:10dB)小さい値をC/N0の下限とし、受信した衛星信号から条件を満たす衛星信号を選択する。すなわち、信号選択部130は、L1帯で受信した全衛星信号から、C/N0値がCN0max-dCN0<C/N0<CN0maxを満たす全ての衛星信号を選択する。
S104において、信号選択部130は、S102とS103で選択された衛星信号の数が、予め設定した最低選択衛星信号数(N0)以上であるか否かを判定する。S104の判定結果がYesであれば、信号選択部130による信号選択処理を終了する。信号選択部130は、選択された衛星信号の識別情報(PRN番号等のコード)を計測部140に通知することで、計測部140は、選択された衛星信号を使用して測位及び時刻同期を行うことができる。
S104の判定結果がNoである場合、すなわち、S102とS103で選択された衛星信号の数が、予め設定した最低選択衛星信号数(N0)未満である場合、S105に進む。
S105において、信号選択部130は、予め設定したGNSS種別の優先順位に基づき、C/N0値が「CN0max-dCN0」以下でC/N0値の大きい次点の衛星信号から順に衛星信号を選択し、トータルの選択衛星信号数がN0となるよう補填する。
図9に、GNSS種別の優先順位の設定例を示す。GNSS種別の優先順位の設定値についてもバイアス値格納部180に格納されており、信号選択部130は、バイアス値格納部180に格納された設定値を参照する。図9において、GPSの優先順位が最も高く、GLO(GLONASS)の優先順位が最も低いことが示されている。
ここで、GNSS種別の優先順位に基づく選択について説明する。GNSS衛星の種別毎に、衛星を運用している時計の絶対時刻を基準とした時刻精度に差分(クロックバイアス)が存在する。S105で補欠の衛星信号を選択する際にはクロックバイアスを含む、GNSSの信頼度を考慮して選択することとしている。
例えば、GPSとQZSSは航法衛星システムとして時刻が互いに完全に同期しており、クロックバイアスが小さいのでカテゴリー1、カテゴリー2がGalileo、カテゴリー3がGLONASS及びBeiDouといった形で分類することができる。このようなカテゴリー分けに基づいて、図9に示すような優先順位の設定がなされる。
上記のような信頼度に基づく優先順位を考慮した補欠衛星信号の選択方法としては、一例として、優先順位(又は信頼度のカテゴリー)に応じてC/N0値にプレミアム(付加する値、単位はdB)を設定し、C/N0値が最高値の衛星信号から必要数を順に選択する方法がある。
例えば、図9の例において、優先順位1のプレミアムを5とし、優先順位2のプレミアムを4とし、優先順位3のプレミアムを3とし、優先順位4のプレミアムを2とし、優先順位5のプレミアムを1とする。
一例として、N0が5であるとし、S102、S103で3個の衛星信号が選択されたとする。また、「CN0max-dCN0」以下のC/N0値を持つ衛星信号として、衛星信号1(C/N0値=26dB-Hz、プレミアム=1)、衛星信号2(C/N0値=25dB-Hz、プレミアム=3)、衛星信号3(C/N0値=24dB-Hz、プレミアム=5)があるとすると、S105において、信号選択部130は、プレミアムを加えたC/N0値が、29dB-Hz、28dB-Hzとなる衛星信号3と衛星信号2を選択する。
(バイアス値設定に関する動作例)
次に、バイアス値を設定するための動作例を、図10に示すフローチャートの手順に沿って詳細に説明する。手順の説明の中で、図11~図13も参照する。
次に、バイアス値を設定するための動作例を、図10に示すフローチャートの手順に沿って詳細に説明する。手順の説明の中で、図11~図13も参照する。
S201において、信号受信部120により、衛星信号データを継続的に収集する。収集する時間長に関して、オープンスカイ環境であれば、連続で24時間収集すれば十分である。それ以外の受信環境では更に長期連続収集が必要である。随時データを収集し、バイアス値を更新してもよい。
S202において、収集した衛星信号データを(GNSS種別、仰角、C/N0値)のセットとして信号データ格納部160に格納する。
S203において、バイアス値設定部170は、信号データ格納部160に格納された衛星信号データに基づいて、同一GNSS種別のデータに対し、仰角の範囲毎にデータをグループ分けし、各グループのC/N0値の最大値を抽出する。
図11に、あるGNSS種別におけるS203の処理の例を示す。図11の例では、仰角が0°~15°、15°~30°、30°~45°、45°~60°、60°~75°、75°~90°にグループ分けされ、各グループのC/N0値の最大値が抽出されている。
S204において、バイアス値設定部170は、抽出した最大値データに対し例えば非線形最小二乗法によりカーブフィッティングを適用する。S205において、バイアス値設定部170は、最も大きい外れ値を除いたカーブフィッティングを数回繰り返す。図11に示したGNSS種別に対するS204、S205の例を図12に示す。
S206において、バイアス値設定部170は、GNSS種別毎にフィッティング関数を生成し、S207において、各GNSS種別のフィッティング関数によりGNSS種別・仰角のバイアス値を設定する。S206、S207の例を図13に示す。図13に示すように、いずれのGNSS種別においても、仰角が小さいほど大きなバイアス値が設定される。また、図13の例において、GNSS種別間では、GNSS-C>GNSS-B>GNSS-Aの順の大きさのバイアス値が設定される。
(dCN0の設定値について)
次に、dCN0の設定値(dCN0値と呼ぶ)について説明する。dCN0値は、図5のS103で説明したとおり、衛星信号を選択するためのC/N0値の範囲を決めるパラメータである。dCN0値は、衛星信号の仰角に依らない固定値であってもよいが、以下では、衛星信号の仰角に依存してdCN0値を決める例について説明する。ここで説明する例は、都市部のように、衛星信号の反射面が建造物の鉛直方向の壁面(コンクリートorガラス)であるような場合を想定した例である。
次に、dCN0の設定値(dCN0値と呼ぶ)について説明する。dCN0値は、図5のS103で説明したとおり、衛星信号を選択するためのC/N0値の範囲を決めるパラメータである。dCN0値は、衛星信号の仰角に依らない固定値であってもよいが、以下では、衛星信号の仰角に依存してdCN0値を決める例について説明する。ここで説明する例は、都市部のように、衛星信号の反射面が建造物の鉛直方向の壁面(コンクリートorガラス)であるような場合を想定した例である。
図14は、仰角の高い衛星信号が建造物の垂直な壁面に入射して反射する様子を示し、図15は、仰角の低い衛星信号が建造物の垂直な壁面に入射して反射する様子を示している。図14、図15に示されるとおり、仰角の低い衛星信号が建造物の垂直な壁面に入射する入射角は、仰角の高い衛星信号が建造物の垂直な壁面に入射する入射角よりも大きい。
衛星信号の建造物の垂直な壁面による反射率は入射角に依存するため、仰角の低い衛星信号は仰角の高い衛星信号と比較して相対的に反射率が大きい(反射波の信号強度が大きい)ことが期待される。
そこで、dCN0値に仰角依存性を与えることが可視/不可視の衛星選択においては有効となる。図16に、仰角依存性を与えたdCN0値の設定例を示す。図16に示すように、衛星信号の仰角が大きくなるとdCN0値も大きくなるように設定する。このような仰角依存性を持つ設定値は、例えば、図16の曲線に相当する関数の形でバイアス値格納部180に保存されてもよいし、各仰角(例えば5°刻み)に対するdCN0値を保持するテーブルの形でバイアス値格納部180に保存されてもよい。
信号選択部130は、前述したS103において、ある衛星信号のC/N0値が「CN0max-dCN0<C/N0<CN0max」を満たすか否かを判断する際に、バイアス値格納部180を参照して、その衛星信号の仰角に対応するdCN0値を取得し、そのdCN0値を用いて「CN0max-dCN0<C/N0<CN0max」を満たすか否かを判断する。
また、前述したS105の補欠の衛星信号選択において、信号選択部130は、ある衛星信号のC/N0値が「CN0max-dCN0」以下か否かを判断する際に、バイアス値格納部180を参照して、その衛星信号の仰角に対応するdCN0値を取得し、そのdCN0値を用いて「CN0max-dCN0」以下か否かを判断する。
「CN0max-dCN0<C/N0<CN0max」により衛星信号を選択するか否かの判断において、仰角の低い衛星信号のほうが仰角の高い衛星信号よりもdCN0値が小さくなるので、仰角の低い衛星信号のほうが仰角の高い衛星信号よりも「CN0max-dCN0<C/N0<CN0max」の範囲が狭くなる。すなわち、仰角の低い衛星信号のほうが仰角の高い衛星信号よりも厳しめのフィルタリングがなされる。このように、仰角の低い衛星信号のほうが仰角の高い衛星信号よりも厳しめのフィルタリングがなされるようにdCN0値に仰角依存性を持たせる理由を以下に説明する。
都市部において衛星信号の反射面が建造物の鉛直方向の壁面(コンクリートorガラス)であると想定すると低仰角の衛星信号は全反射に近い状態になり、反射した衛星信号の信号強度と、仮に障害物が存在せずに直接波として受信した場合の信号強度(図7、図8のバイアス値により正規化される基準信号強度)との差分は小さくなる。
つまり、低仰角の衛星信号は電離層や対流圏といった信号強度を減衰させる媒質の光路長が長くなる分、直接波として受信した場合の信号強度は小さくなるが、一方で、建造物で反射した際の信号強度の低下は小さいため、不可視衛星信号のマルチパス信号(反射波)を除去する上ではdCN0値を小さくし、厳しめにフィルタリングする必要がある。高仰角衛星はその逆になり、「CN0max-dCN0<C/N0<CN0max」の範囲を広げ、選択されやすくする。
なお、衛星信号の仰角が大きくなるとdCN0値も大きくなるように仰角依存性を与えることは例である。環境によっては、上記とは異なる仰角依存性をdCN0値に与えることとしてもよい。
なお、各衛星種別単位での衛星選択を行うことも考えられるが、それは行わない。その理由は下記のとおりである。
本発明に係る技術では、少なくとも1つの可視衛星が存在することが前提となる。各衛星種別単位で衛星選択を行った場合、ある衛星種別で可視衛星が存在しなかった場合には基準C/N0値(CN0max)が不適切な値となり、衛星選択の精度が劣化する可能性がある。第1及び第2の実施の形態のように、全衛星種別を対象とすれば、少なくとも1つの可視衛星が存在する確率が向上する。
―――――――――――第2の実施の形態―――――――――――
次に、第2の実施の形態について説明する。第2の実施の形態は、計測装置100が、周波数帯毎に衛星信号の選択を行う点が第1の実施の形態と異なる。すなわち、第1の実施の形態では、例としてL1帯のみを対象として衛星信号の選択を行っていたが、第2の実施の形態では、各衛星の出力する複数の周波数帯のそれぞれに対し、衛星信号の選択を行う。
なお、第1の実施の形態で説明した技術により、発明の効果を発揮することができる。第2の実施の形態は、発明の実施形態のバリエーションである。第2の実施の形態において、周波数帯毎に衛星信号の選択を行う理由は下記のとおりである。
各衛星は複数の周波数帯の信号を出力しているため、衛星の位置(可視/不可視)により、測位に適する衛星を選択するという観点では、何れか1つの周波数帯の衛星信号により衛星の可視/不可視が正確に判定できるのであれば、複数の周波数帯それぞれの信号で衛星選択を行う必要性はない。
しかし、実際には可視/不可視を100%正確に判定できることは保証されない。また、周波数帯毎にアンテナ・レシーバの受信特性や干渉信号の混入状態が異なる可能性があり、複数の周波数帯でそれぞれ衛星信号の選択を行うことで測位演算により適した衛星信号の組み合わせを選択する可能性がある。
測位演算においては周波数帯毎に異なる衛星信号を使用することができる。衛星によりサポートする周波数帯が異なる(例えばGPSのL5周波数帯は一部の衛星のみ対応)ため、周波数帯毎に衛星信号を個別に選択することによって、測位演算のポリシー設定のバリエーション(周波数帯毎にN0値を変える等)を広げることができる。
第2の実施の形態における計測装置100の装置構成は第1の実施の形態における装置構成と同じであり、図1に示したとおりである。各部の動作も基本的に第1の実施の形態と同じであるが、周波数帯毎に衛星信号の選択を行うための動作を行う点が第1の実施の形態と異なる。
すなわち、信号受信部120は、周波数帯毎に、受信した衛星信号の識別情報(PRN番号等のコード)、当該衛星信号の仰角、C/N0値、及び衛星種別を信号選択部130に送る。また、信号受信部120は、周波数帯毎に、受信した衛星信号毎の識別情報、仰角、C/N0値、衛星種別を信号データ格納部160に格納する。本実施の形態では、複数の周波数帯としてL1帯とL2帯を対象とする。ただし、L1帯とL2帯を使用することは例であり、これらに加えてL5帯を使用してもよいし、L1帯、L2帯、L5帯以外の周波数帯を使用してもよい。
(第2の実施の形態における信号選択部130の動作例)
次に、第2の実施の形態における信号選択部130の動作例を説明する。図17は、信号選択部130の動作を示すフローチャートである。基本的に図5に示した第1の実施の形態でのフローと同じであるが、第2の実施の形態では、図17のフローを周波数帯毎に繰り返すとともに、S113(図5のS103に対応)において、処理中の周波数帯における最低C/N0値を満たすかどうかの判定を行う点が第1の実施の形態と異なる。なお、図17は、例として、周波数帯毎の繰り返しにおけるL1帯についての処理を示している。
次に、第2の実施の形態における信号選択部130の動作例を説明する。図17は、信号選択部130の動作を示すフローチャートである。基本的に図5に示した第1の実施の形態でのフローと同じであるが、第2の実施の形態では、図17のフローを周波数帯毎に繰り返すとともに、S113(図5のS103に対応)において、処理中の周波数帯における最低C/N0値を満たすかどうかの判定を行う点が第1の実施の形態と異なる。なお、図17は、例として、周波数帯毎の繰り返しにおけるL1帯についての処理を示している。
まず、図18を参照して、第2の実施の形態の手順において使用される設定パラメータについて説明する。図18に示すとおり、CN0L1は、L1帯の選択衛星の最低C/N0値である。dCN0L1は、L1帯の衛星信号の選択範囲を決めるパラメータである(前述したdCN0に相当)。N0L1は、L1帯の選択衛星信号数である。L2帯の関しても同様のパラメータが設定されている。なお、他の周波数帯を使う場合、周波数帯毎にパラメータを設定すればよい。例えば、L5帯であればCN0L5等を設定する。
まず、L1帯について、図17のフローの処理を実行する。S101、S102の処理は第1の実施の形態と同じである。ただし、第2の実施の形態では、周波数帯毎に、図7、図8に示したようなバイアス値が設定されており、S101の正規化処理では、処理中の周波数帯に対応するバイアス値を用いて正規化を行う。S102では、該当周波数帯(最初はL1帯)で受信した全衛星信号の中からC/N0値の最も大きい衛星信号を選択し、そのC/N0値をCN0maxとして記録する。
S113において、信号選択部130は、S102で選択した衛星信号のC/N0値(CN0max)に対し、CN0maxよりdCN0L1(例:10dB)小さい値をC/N0の下限とし、受信した衛星信号から条件を満たす衛星信号を選択する。ここでは、信号選択部130は、L1帯で受信した全衛星信号から、C/N0値がCN0max-dCN0L1<C/N0<CN0maxを満たし、かつ、CN0L1<C/N0を満たす全ての衛星信号を選択する。
S104において、信号選択部130は、S102とS103で選択された衛星信号の数が、予め設定した最低選択衛星信号数(N0L1)以上であるか否かを判定する。S104の判定結果がYesであれば、L1帯についての信号選択部130による信号選択処理を終了する。信号選択部130は、選択された衛星信号の識別情報(PRN番号等のコード)を計測部140に通知することで、計測部140は、選択された衛星信号を使用して測位及び時刻同期を行うことができる。
信号選択部130は、次の周波数帯(例えばL2帯)について、図17のフローの処理を実行する。
S104の判定結果がNoである場合、すなわち、S102とS103で選択された衛星信号の数が、予め設定した最低選択衛星信号数(N0L1)未満である場合、S105に進む。
S105の処理は、第1の実施の形態で説明した処理と同じでよい。すなわち、S105において、信号選択部130は、予め設定したGNSS種別の優先順位(例:図9)に基づき、C/N0値が「CN0max-dCN0L1」以下でC/N0値の大きい次点の衛星信号から順に衛星信号を選択し、トータルの選択衛星信号数がN0L1となるよう補填する。
なお、図9に示したような優先順位の設定は、周波数帯毎に定められてもよい。その場合、信号選択部130は、処理中の周波数帯に対応する優先順位を用いて、衛星信号の選択を行う。
第1の実施の形態と同様に、信頼度に基づく優先順位を考慮した補欠衛星信号の選択方法として、優先順位(又は信頼度のカテゴリー)に応じてC/N0値にプレミアム(付加する値)を設定し、C/N0値が最高値の衛星信号から必要数を順に選択する方法を用いることができる。
例えば、図9の例において、優先順位1のプレミアムを5とし、優先順位2のプレミアムを4とし、優先順位3のプレミアムを3とし、優先順位4のプレミアムを2とし、優先順位5のプレミアムを1とする。
一例として、N0L1が5であるとし、S102、S103で3個の衛星信号が選択されたとする。また、「CN0max-dCN0L1」以下のC/N0値を持つ衛星信号として、衛星信号1(C/N0値=26dB-Hz、プレミアム=1)、衛星信号2(C/N0値=25dB-Hz、プレミアム=3)、衛星信号3(C/N0値=24dB-Hz、プレミアム=5)があるとすると、S105において、信号選択部130は、プレミアムを加えたC/N0値が、29dB-Hz、28dB-Hzとなる衛星信号3と衛星信号2を選択する。
S105が終了すると、信号選択部130は、選択された衛星信号の識別情報(PRN番号等のコード)を計測部140に通知することで、計測部140は、選択された衛星信号を使用して測位及び時刻同期を行うことができる。
信号選択部130は、次の周波数帯(例えばL2帯)について、図17のフローの処理を実行する。
なお、上記の例では、周波数帯毎に、選択された衛星信号を用いて測位及び時刻同期を行うこととしているが、特定の周波数帯で選択された衛星信号に基づき複数の周波数帯の衛星信号を使用する測位及び時刻同期を行うこととしてもよい。
例えば、L1帯とL2帯のそれぞれで図14のフローを実行することで、L1帯において衛星信号1、2、3、4が選択され、L2帯において衛星信号5、6、7、8が選択された場合において、衛星信号1、2、3、4のDOP(Dilution of Precision)値と、衛星信号5、6、7、8のDOP値とを比較して、DOP値が小さいほうの周波数帯の衛星信号を選択して、L1帯とL2帯の信号を使用した測位及び時刻同期を行うこととしてもよい。
(補欠衛星信号の選択方法のバリエーション)
第1の実施の形態及び第2の実施の形態において、S105における補欠衛星信号の選択についての他の例を説明する。信号選択部130は、DOP値を考慮して、補欠の衛星信号を選択してもよい。例えば、信号選択部130は、C/N0値が「CN0max-dCN0」以下でC/N0値の大きい次点の衛星信号から順に衛星信号A、B、Cを、補欠衛星信号候補として選択し、既に選択されている衛星信号に補欠衛星候補A、B、Cをそれぞれ加えた場合のDOP値を計算し、DOP値が最も小さくなる衛星信号を選択する。
第1の実施の形態及び第2の実施の形態において、S105における補欠衛星信号の選択についての他の例を説明する。信号選択部130は、DOP値を考慮して、補欠の衛星信号を選択してもよい。例えば、信号選択部130は、C/N0値が「CN0max-dCN0」以下でC/N0値の大きい次点の衛星信号から順に衛星信号A、B、Cを、補欠衛星信号候補として選択し、既に選択されている衛星信号に補欠衛星候補A、B、Cをそれぞれ加えた場合のDOP値を計算し、DOP値が最も小さくなる衛星信号を選択する。
例えば、S104の時点で既に選択されている衛星信号が衛星信号1、2、3であるとすると、「衛星信号1、2、3、A」、「衛星信号1、2、3、B」、「衛星信号1、2、3、C」のそれぞれのDOP値を計算する。もしも「衛星信号1、2、3、A」のDOP値が最小であれば、「衛星信号1、2、3、A」を選択する。選択すべき衛星数が、4よりも大きければ、上記の処理をその衛星数になるまで繰り返せばよい。
ここで、第1の実施の形態で説明した補欠衛星信号選択方法を選択方法1とし、DOP値を用いる上記の方法を選択方法2とする。信号選択部130は、選択方法1、2の組み合わせで補欠衛星信号を選択してもよい。
組み合わせの例としては、選択方法1により、補欠の衛星信号を選択し、その選択された衛星信号のそれぞれで選択方法2を実施して、DOP値が小さい衛星信号を選択する。
また、例えば、選択方法1、2に基づき選択された衛星信号による測位精度の改善度合いを期待値としてコスト値(評価値)を設定し、選択方法1と2のトータルのコスト値(評価値)が最小となる衛星信号を補欠衛星信号として選択してもよい。例えば、選択方法1、2の組み合わせにより補欠衛星信号の候補として衛星信号Aと衛星信号Bが選択されたとする。例えば、衛星信号AのC/N0値が30dB-Hz、衛星信号BのC/N0値が28dB-Hzであり、衛星信号Aを選択した場合のDOP値が5、衛星信号Bを選択した場合のDOP値が4であったとする。選択方法1と2のトータルのコスト値(評価値)を「DOP値÷(C/N0値)」として設定した場合、衛星信号A、Bのコスト値はそれぞれ1/6、1/7となり、衛星信号Aよりも衛星信号Bのほうが小さいため、衛星信号Bを補欠衛星信号として選択する。
(バイアス値設定に関する動作例)
第2の実施の形態において、バイアス値設定部170により実行されるバイアス値設定動作は、基本的に第1の実施の形態でのバイアス値設定動作と同じであるが、第2の実施の形態では、周波数帯毎にバイアス値を設定する点が第1の実施の形態と異なる。
第2の実施の形態において、バイアス値設定部170により実行されるバイアス値設定動作は、基本的に第1の実施の形態でのバイアス値設定動作と同じであるが、第2の実施の形態では、周波数帯毎にバイアス値を設定する点が第1の実施の形態と異なる。
図19に、第2の実施の形態におけるバイアス値設定動作のフローチャートを示す。S201では、信号受信部120により、周波数帯毎に、衛星信号データを継続的に収集する。
S212において、収集した衛星信号データを(GNSS種別、周波数帯、仰角、C/N0値)のセットとして信号データ格納部160に格納する。
S213において、バイアス値設定部170は、信号データ格納部160に格納された衛星信号データに基づいて、同一GNSS種別・周波数帯のデータに対し、仰角の範囲毎にデータをグループ分けし、各グループのC/N0値の最大値を抽出する。その処理例は図8を参照して説明したとおりである。
S204において、バイアス値設定部170は、抽出した最大値データに対し非線形最小二乗法等によりカーブフィッティングを適用する。S205において、バイアス値設定部170は、最も大きい外れ値を除いたカーブフィッティングを数回繰り返す。図11に示したGNSS種別に対するS204、S205の例は、図12に示したとおりである。
S206において、バイアス値設定部170は、GNSS種別毎にフィッティング関数を生成し、S207において、各GNSS種別のフィッティング関数によりGNSS種別・仰角のバイアス値を設定する。L1帯に対する、S206、S207の例を図20に示す。
上述したとおり、第2の実施の形態では、GNSS衛星種別に対して各衛星信号の周波数帯(例えばGPSの場合はL1帯、L2帯、L5帯)毎にバイアス値を設定する。それは、衛星信号の受信特性がGNSS衛星種別、仰角以外に、衛星信号の周波数帯にも依存するためである。なお、仰角を用いずに、GNSS衛星種別・周波数帯のみに基づくバイアス値を設定してもよい。
同一のGNSSアンテナとGNSSレシーバの組み合わせに対する、周波数帯による受信特性の違いの実測例を図21、図22に示す。図21は、GPSのL1信号を示し、図22は、GPSのL2信号を示す。図21、図22ともに、横軸は仰角(°)であり、縦軸はC/N0値(dB-Hz)である。
(バイアス値設定に関するバリエーション)
同じGNSS種別でも個別の衛星により送信信号出力が異なるケースがある。例えば、衛星の軌道(GEO/IGSO/MEO)により送信信号出力が異なる場合がある。その場合、個別の衛星に対しバイアス値を設定してもよい。
同じGNSS種別でも個別の衛星により送信信号出力が異なるケースがある。例えば、衛星の軌道(GEO/IGSO/MEO)により送信信号出力が異なる場合がある。その場合、個別の衛星に対しバイアス値を設定してもよい。
例えば、衛星Aの送信信号強度が同じGNSS種別の他の衛星の信号強度より小さい場合、受信信号強度が小さくなり直接波として受信した場合でも選択されないケースがある。その場合に正規化の際に当該衛星信号については、受信品質に個別のバイアス値を加える補正を実施する。個別のバイアス値は、図7、図8に示したGNSSバイアス値と仰角バイアス値に追加して適用する。あるいは、個別のバイアス値を適用する衛星信号に対しては、GNSSバイアス値と仰角バイアス値を適用せずに、個別のバイアス値のみを適用してもよい。あるいは、個別のバイアス値を適用する衛星信号に対しては、GNSSバイアス値を適用せずに、仰角バイアス値と個別のバイアス値を適用してもよい。
どのGNSS種別のどの衛星に対して個別のバイアス値を設定するかについては、例えば、事前に、GNSS種別毎・衛星毎に、受信信号強度の測定を行って、測定値を信号データ格納部160に格納し、バイアス値設定部170が、上記の衛星Aと同様の事象が生じる衛星を選択し、選択した衛星に対して個別のバイアス値を設定する。
また、特定の衛星に対する個別のバイアス値を、上記の衛星Aのような送信信号出力に関わるケース以外のケースに適用してもよい。
(dCN0の設定値について)
第2の実施の形態における周波数帯毎のdCN0値(つまり、dCN0L1値、dCN0L2値など)に関しても、図14、図15を参照して説明したように、仰角依存性を与えたdCN0値を設定してもよい。この場合、例えば、図16に示したような仰角依存性を与えたdCN0値を周波数帯毎に設定する。
第2の実施の形態における周波数帯毎のdCN0値(つまり、dCN0L1値、dCN0L2値など)に関しても、図14、図15を参照して説明したように、仰角依存性を与えたdCN0値を設定してもよい。この場合、例えば、図16に示したような仰角依存性を与えたdCN0値を周波数帯毎に設定する。
(ハードウェア構成例)
図23は、本実施の形態における計測装置100として使用することができるコンピュータのハードウェア構成例を示す図である。当該コンピュータは、物理的な装置としてのコンピュータであってもよいし、クラウド上の仮想マシンであってもよい。
図23は、本実施の形態における計測装置100として使用することができるコンピュータのハードウェア構成例を示す図である。当該コンピュータは、物理的な装置としてのコンピュータであってもよいし、クラウド上の仮想マシンであってもよい。
図23のコンピュータは、それぞれバスBで相互に接続されているドライブ装置1000、補助記憶装置1002、メモリ装置1003、CPU1004、インタフェース装置1005、表示装置1006、入力装置1007、及び出力装置1008等を有する。なお、図20にはGNSSアンテナ110は示されていない。GNSSアンテナ110は、例えば、インタフェース装置1005に接続される。
当該コンピュータでの処理を実現するプログラムは、例えば、CD-ROM又はメモリカード等の記録媒体1001によって提供される。プログラムを記憶した記録媒体1001がドライブ装置1000にセットされると、プログラムが記録媒体1001からドライブ装置1000を介して補助記憶装置1002にインストールされる。但し、プログラムのインストールは必ずしも記録媒体1001より行う必要はなく、ネットワークを介して他のコンピュータよりダウンロードするようにしてもよい。補助記憶装置1002は、インストールされたプログラムを格納すると共に、必要なファイルやデータ等を格納する。
メモリ装置1003は、プログラムの起動指示があった場合に、補助記憶装置1002からプログラムを読み出して格納する。CPU1004は、メモリ装置1003に格納されたプログラムに従って、計測装置100に係る機能を実現する。インタフェース装置1005は、GNSSアンテナ110に接続するためのインタフェースとして用いられる。表示装置1006はプログラムによるGUI(Graphical User Interface)等を表示する。入力装置1007はキーボード及びマウス、ボタン、又はタッチパネル等で構成され、様々な操作指示を入力させるために用いられる。出力装置1008は演算結果を出力する。
(実施の形態の効果)
以上説明したように、本発明の実施の形態によれば、受信器やアンテナの特性に依存せず、受信した衛星信号から可視衛星信号を高い精度で選択することができる。これにより、良好ではない受信環境における測位精度を向上することができる。また、衛星の選択における干渉信号の影響を低減することができる。
以上説明したように、本発明の実施の形態によれば、受信器やアンテナの特性に依存せず、受信した衛星信号から可視衛星信号を高い精度で選択することができる。これにより、良好ではない受信環境における測位精度を向上することができる。また、衛星の選択における干渉信号の影響を低減することができる。
(付記)
本実施の形態において、少なくとも、下記の各項に記載された計測装置、計測方法、及びプログラムが提供される。
(第1項)
GNSSアンテナにより受信した衛星信号の受信品質に基づいて、所定数の衛星信号を選択する信号選択部と、
前記信号選択部により選択された所定数の衛星信号を用いて測位を実行し、初期推定位置を決定する計測部と、を備え、
前記信号選択部は、衛星信号を初期推定位置において直接波として受信した場合に期待される受信時刻と当該衛星信号の実際の受信時刻とを比較した比較結果に基づいて、測位又は時刻同期に使用する複数の衛星信号を選択する
衛星信号受信装置。
(第2項)
前記信号選択部は、前記比較結果に基づいて、初期推定位置において直接波として受信した場合に期待される受信時刻から閾値以上の遅延時間で実際に受信された衛星信号を棄却する
第1項に記載の衛星信号受信装置。
(第3項)
前記計測部は、棄却された前記衛星信号を除く複数の衛星信号を用いて測位を行い、得られた座標値で初期推定位置を更新する
第2項に記載の衛星信号受信装置。
(第4項)
前記信号選択部は、初期推定位置を更新し、更新した初期推定位置を用いた前記比較結果に基づく処理を繰り返し、所定の終了条件を満足した場合に、当該繰り返しの処理を終了する
第3項に記載の衛星信号受信装置。
(第5項)
前記繰り返しの処理を終了した時点での選択された衛星信号の数が予め設定した数よりも小さい場合において、前記信号選択部は、複数の非選択衛星信号の間で遅延時間を比較し、遅延時間の小さい非選択衛星信号を前記選択された衛星信号に加えた場合のDOP値に基づいて、前記選択された衛星信号に追加する非選択衛星信号を決定する
第4項に記載の衛星信号受信装置。
(第6項)
衛星信号受信装置が実行する衛星信号選択方法であって、
GNSSアンテナにより受信した衛星信号の受信品質に基づいて、所定数の衛星信号を選択するステップと、
前記所定数の衛星信号を用いて測位を実行し、初期推定位置を決定するステップと、
衛星信号を初期推定位置において直接波として受信した場合に期待される受信時刻と当該衛星信号の実際の受信時刻とを比較した比較結果に基づいて、測位又は時刻同期に使用する複数の衛星信号を選択するステップと
を備える衛星信号選択方法。
(第7項)
コンピュータを、第1項ないし第5項のうちいずれか1項に記載の衛星信号受信装置における各部として機能させるためのプログラム。
本実施の形態において、少なくとも、下記の各項に記載された計測装置、計測方法、及びプログラムが提供される。
(第1項)
GNSSアンテナにより受信した衛星信号の受信品質に基づいて、所定数の衛星信号を選択する信号選択部と、
前記信号選択部により選択された所定数の衛星信号を用いて測位を実行し、初期推定位置を決定する計測部と、を備え、
前記信号選択部は、衛星信号を初期推定位置において直接波として受信した場合に期待される受信時刻と当該衛星信号の実際の受信時刻とを比較した比較結果に基づいて、測位又は時刻同期に使用する複数の衛星信号を選択する
衛星信号受信装置。
(第2項)
前記信号選択部は、前記比較結果に基づいて、初期推定位置において直接波として受信した場合に期待される受信時刻から閾値以上の遅延時間で実際に受信された衛星信号を棄却する
第1項に記載の衛星信号受信装置。
(第3項)
前記計測部は、棄却された前記衛星信号を除く複数の衛星信号を用いて測位を行い、得られた座標値で初期推定位置を更新する
第2項に記載の衛星信号受信装置。
(第4項)
前記信号選択部は、初期推定位置を更新し、更新した初期推定位置を用いた前記比較結果に基づく処理を繰り返し、所定の終了条件を満足した場合に、当該繰り返しの処理を終了する
第3項に記載の衛星信号受信装置。
(第5項)
前記繰り返しの処理を終了した時点での選択された衛星信号の数が予め設定した数よりも小さい場合において、前記信号選択部は、複数の非選択衛星信号の間で遅延時間を比較し、遅延時間の小さい非選択衛星信号を前記選択された衛星信号に加えた場合のDOP値に基づいて、前記選択された衛星信号に追加する非選択衛星信号を決定する
第4項に記載の衛星信号受信装置。
(第6項)
衛星信号受信装置が実行する衛星信号選択方法であって、
GNSSアンテナにより受信した衛星信号の受信品質に基づいて、所定数の衛星信号を選択するステップと、
前記所定数の衛星信号を用いて測位を実行し、初期推定位置を決定するステップと、
衛星信号を初期推定位置において直接波として受信した場合に期待される受信時刻と当該衛星信号の実際の受信時刻とを比較した比較結果に基づいて、測位又は時刻同期に使用する複数の衛星信号を選択するステップと
を備える衛星信号選択方法。
(第7項)
コンピュータを、第1項ないし第5項のうちいずれか1項に記載の衛星信号受信装置における各部として機能させるためのプログラム。
以上、本実施の形態について説明したが、本発明はかかる特定の実施形態に限定されるものではなく、特許請求の範囲に記載された本発明の要旨の範囲内において、種々の変形・変更が可能である。
100 計測装置
110 GNSSアンテナ
120 信号受信部
130 信号選択部
140 計測部
150 出力部
160 信号データ格納部
170 バイアス値設定部
180 バイアス値格納部
1000 ドライブ装置
1001 記録媒体
1002 補助記憶装置
1003 メモリ装置
1004 CPU
1005 インタフェース装置
1006 表示装置
1007 入力装置
1008 出力装置
110 GNSSアンテナ
120 信号受信部
130 信号選択部
140 計測部
150 出力部
160 信号データ格納部
170 バイアス値設定部
180 バイアス値格納部
1000 ドライブ装置
1001 記録媒体
1002 補助記憶装置
1003 メモリ装置
1004 CPU
1005 インタフェース装置
1006 表示装置
1007 入力装置
1008 出力装置
Claims (7)
- GNSSアンテナにより受信した衛星信号の受信品質に基づいて、所定数の衛星信号を選択する信号選択部と、
前記信号選択部により選択された所定数の衛星信号を用いて測位を実行し、初期推定位置を決定する計測部と、を備え、
前記信号選択部は、衛星信号を初期推定位置において直接波として受信した場合に期待される受信時刻と当該衛星信号の実際の受信時刻とを比較した比較結果に基づいて、測位又は時刻同期に使用する複数の衛星信号を選択する
衛星信号受信装置。 - 前記信号選択部は、前記比較結果に基づいて、初期推定位置において直接波として受信した場合に期待される受信時刻から閾値以上の遅延時間で実際に受信された衛星信号を棄却する
請求項1に記載の衛星信号受信装置。 - 前記計測部は、棄却された前記衛星信号を除く複数の衛星信号を用いて測位を行い、得られた座標値で初期推定位置を更新する
請求項2に記載の衛星信号受信装置。 - 前記信号選択部は、初期推定位置を更新し、更新した初期推定位置を用いた前記比較結果に基づく処理を繰り返し、所定の終了条件を満足した場合に、当該繰り返しの処理を終了する
請求項3に記載の衛星信号受信装置。 - 前記繰り返しの処理を終了した時点での選択された衛星信号の数が予め設定した数よりも小さい場合において、前記信号選択部は、複数の非選択衛星信号の間で遅延時間を比較し、遅延時間の小さい非選択衛星信号を前記選択された衛星信号に加えた場合のDOP値に基づいて、前記選択された衛星信号に追加する非選択衛星信号を決定する
請求項4に記載の衛星信号受信装置。 - 衛星信号受信装置が実行する衛星信号選択方法であって、
GNSSアンテナにより受信した衛星信号の受信品質に基づいて、所定数の衛星信号を選択するステップと、
前記所定数の衛星信号を用いて測位を実行し、初期推定位置を決定するステップと、
衛星信号を初期推定位置において直接波として受信した場合に期待される受信時刻と当該衛星信号の実際の受信時刻とを比較した比較結果に基づいて、測位又は時刻同期に使用する複数の衛星信号を選択するステップと
を備える衛星信号選択方法。 - コンピュータを、請求項1ないし5のうちいずれか1項に記載の衛星信号受信装置における各部として機能させるためのプログラム。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2022/000721 WO2023135683A1 (ja) | 2022-01-12 | 2022-01-12 | 衛星信号受信装置、衛星信号選択方法、及びプログラム |
PCT/JP2022/034447 WO2023135863A1 (ja) | 2022-01-12 | 2022-09-14 | 衛星信号受信装置、衛星信号選択方法、及びプログラム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2022/000721 WO2023135683A1 (ja) | 2022-01-12 | 2022-01-12 | 衛星信号受信装置、衛星信号選択方法、及びプログラム |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2023135683A1 true WO2023135683A1 (ja) | 2023-07-20 |
Family
ID=87278619
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2022/000721 WO2023135683A1 (ja) | 2022-01-12 | 2022-01-12 | 衛星信号受信装置、衛星信号選択方法、及びプログラム |
PCT/JP2022/034447 WO2023135863A1 (ja) | 2022-01-12 | 2022-09-14 | 衛星信号受信装置、衛星信号選択方法、及びプログラム |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2022/034447 WO2023135863A1 (ja) | 2022-01-12 | 2022-09-14 | 衛星信号受信装置、衛星信号選択方法、及びプログラム |
Country Status (1)
Country | Link |
---|---|
WO (2) | WO2023135683A1 (ja) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006132949A (ja) * | 2004-11-02 | 2006-05-25 | Seiko Epson Corp | 測位装置、測位方法および測位プログラム |
JP2010151725A (ja) * | 2008-12-26 | 2010-07-08 | Toyota Motor Corp | Gnss受信装置及び測位方法 |
JP2010534849A (ja) * | 2007-07-31 | 2010-11-11 | テレ アトラス ベスローテン フエンノートシャップ | 位置を判定する方法及び装置 |
US20130321202A1 (en) * | 2012-05-30 | 2013-12-05 | Csr Technology Inc. | Methods and Apparatuses for Multipath Estimation and Correction in GNSS Navigation Systems |
JP2020112494A (ja) * | 2019-01-15 | 2020-07-27 | 株式会社豊田中央研究所 | 衛星選択装置、及びプログラム |
JP2020201282A (ja) * | 2018-07-13 | 2020-12-17 | 日本電信電話株式会社 | 航法衛星システム受信装置、その航法衛星信号処理方法及びプログラム |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4196781B2 (ja) * | 2003-08-01 | 2008-12-17 | セイコーエプソン株式会社 | 測位装置、測位装置の制御方法、及び、プログラム |
JP2008045896A (ja) * | 2006-08-11 | 2008-02-28 | Seiko Epson Corp | プログラム、測位回路及び電子機器 |
WO2008025151A1 (en) * | 2006-08-31 | 2008-03-06 | Sige Semiconductor Inc. | An apparatus and method for use in global position measurements |
EP2409289A1 (en) * | 2009-03-16 | 2012-01-25 | Tomtom Belgium N.V. | Method for updating digital maps |
JP2010256301A (ja) * | 2009-04-28 | 2010-11-11 | Toyota Central R&D Labs Inc | マルチパス判定装置及びプログラム |
US11899117B2 (en) * | 2018-03-09 | 2024-02-13 | Nec Solution Innovators, Ltd. | Moving body positioning system, method, and program |
-
2022
- 2022-01-12 WO PCT/JP2022/000721 patent/WO2023135683A1/ja unknown
- 2022-09-14 WO PCT/JP2022/034447 patent/WO2023135863A1/ja active Application Filing
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006132949A (ja) * | 2004-11-02 | 2006-05-25 | Seiko Epson Corp | 測位装置、測位方法および測位プログラム |
JP2010534849A (ja) * | 2007-07-31 | 2010-11-11 | テレ アトラス ベスローテン フエンノートシャップ | 位置を判定する方法及び装置 |
JP2010151725A (ja) * | 2008-12-26 | 2010-07-08 | Toyota Motor Corp | Gnss受信装置及び測位方法 |
US20130321202A1 (en) * | 2012-05-30 | 2013-12-05 | Csr Technology Inc. | Methods and Apparatuses for Multipath Estimation and Correction in GNSS Navigation Systems |
JP2020201282A (ja) * | 2018-07-13 | 2020-12-17 | 日本電信電話株式会社 | 航法衛星システム受信装置、その航法衛星信号処理方法及びプログラム |
JP2020112494A (ja) * | 2019-01-15 | 2020-07-27 | 株式会社豊田中央研究所 | 衛星選択装置、及びプログラム |
Non-Patent Citations (1)
Title |
---|
YOSHIDA, SEIJI ET AL.: "Improvement of the time synchronization precision in the multipath satellite signal reception environment by an algorithm for statistical satellite signals selection applied in the GPS receiver", IEICE TECHNICAL REPORT, vol. 116, no. 468, 16 February 2017 (2017-02-16), pages 7 - 12, XP055676183, ISSN: 0913-5685 * |
Also Published As
Publication number | Publication date |
---|---|
WO2023135863A1 (ja) | 2023-07-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20220404511A1 (en) | Method and system for calibrating a system parameter | |
RU2591953C2 (ru) | Навигационная система и способ разрешения целочисленных неоднозначностей с использованием ограничения неоднозначности двойной разности | |
US9304184B1 (en) | System and method for mitigating severe multipath interference for geolocation and navigation | |
CN102597701B (zh) | 用于补偿错误测量的***和方法 | |
EP1727300A1 (en) | Electric wave propagation characteristic estimation system, method thereof, and program | |
US20230251388A1 (en) | Navigation satellite system reception device, method for processing navigation satellite signal from same, and program | |
US9612337B2 (en) | Antenna beam forming for tracking a transmitter signal | |
EP1148344A1 (en) | Positioning of a wireless terminal with satellite positioning signals or base station signals | |
US11442173B2 (en) | GNSS-based timing synchronization | |
JP2008527364A (ja) | 位置決定方法及び装置 | |
WO2003023440A2 (en) | System and method to estimate the location of a receiver in a multi-path environment | |
KR20150130989A (ko) | 위치 추정에서 멀티패스 효과를 방지하기 위해 송신 다양화를 제공하는 시스템 및 방법 | |
CN111988730A (zh) | 利用无线信号进行定位的方法和设备 | |
KR100687243B1 (ko) | 다중 경로 오차 제거를 위한 코드추적루프 및 다중 경로오차 제거 방법 | |
CA2824680A1 (en) | Process and system to determine temporal changes in retransmission and propagation of signals used to measure distances, syncronize actuators and georeference applications | |
CN107003412B (zh) | 定位方法、定位***、校正信息生成方法、校正信息生成装置以及定位***中的中继站、终端 | |
US11340353B2 (en) | Multipath mitigation for GNSS receivers | |
JP2006258461A (ja) | 測位装置、測位方法および測位プログラム | |
KR101981225B1 (ko) | Gnss 측정치에 포함된 특이신호요소 정밀추정 방법 | |
WO2023135683A1 (ja) | 衛星信号受信装置、衛星信号選択方法、及びプログラム | |
KR100899545B1 (ko) | 항법 위성 전체 관측법을 이용한 시각 전송방법 | |
KR20050035900A (ko) | 간섭하는 다중경로 상태를 검출하는 처리절차 | |
JP7485081B2 (ja) | 衛星信号受信装置、衛星信号処理方法、及びプログラム | |
US8436773B2 (en) | Method for leveraging diversity for enhanced location determination | |
JP7315349B2 (ja) | 位置推定装置、及び位置推定方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 22920209 Country of ref document: EP Kind code of ref document: A1 |