WO2019008845A1 - 測位装置 - Google Patents

測位装置 Download PDF

Info

Publication number
WO2019008845A1
WO2019008845A1 PCT/JP2018/012620 JP2018012620W WO2019008845A1 WO 2019008845 A1 WO2019008845 A1 WO 2019008845A1 JP 2018012620 W JP2018012620 W JP 2018012620W WO 2019008845 A1 WO2019008845 A1 WO 2019008845A1
Authority
WO
WIPO (PCT)
Prior art keywords
positioning
unit
satellite
environment
category
Prior art date
Application number
PCT/JP2018/012620
Other languages
English (en)
French (fr)
Inventor
正剛 隈部
ワタンユー チットマン
Original Assignee
株式会社デンソー
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 株式会社デンソー filed Critical 株式会社デンソー
Publication of WO2019008845A1 publication Critical patent/WO2019008845A1/ja
Priority to US16/708,707 priority Critical patent/US11221415B2/en

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/24Acquisition or tracking or demodulation of signals transmitted by the system
    • G01S19/28Satellite selection
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/03Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
    • G01S19/05Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing aiding data
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/22Multipath-related issues
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/24Acquisition or tracking or demodulation of signals transmitted by the system
    • G01S19/25Acquisition or tracking or demodulation of signals transmitted by the system involving aiding data received from a cooperating element, e.g. assisted GPS
    • G01S19/254Acquisition or tracking or demodulation of signals transmitted by the system involving aiding data received from a cooperating element, e.g. assisted GPS relating to Doppler shift of satellite signals
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/03Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
    • G01S19/10Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing dedicated supplementary positioning signals
    • G01S19/11Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing dedicated supplementary positioning signals wherein the cooperating elements are pseudolites or satellite radio beacon positioning system signal repeaters
    • G01S19/115Airborne or satellite based pseudolites or repeaters
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42Determining position
    • G01S19/51Relative positioning

Definitions

  • the present disclosure relates to a positioning device that performs positioning using a positioning signal transmitted from a positioning satellite.
  • the calculation result is obtained because the positioning signal transmitted from the positioning satellite is reflected by a building or the like and received by the positioning device (so-called multi-path) It is known that the position information as can contain an error.
  • GNSS Global Positioning Satellite System
  • a residual which is hereinafter referred to as a pseudo-range residual
  • the pseudorange residual is calculated from the distance (coordinate origin value) determined based on the coordinates of the current position of the positioning device to be estimated and the position coordinates of the positioning satellite, and the reception result of the positioning signal. Is the difference between the pseudoranges as observed values.
  • pseudo-range residuals are sequentially calculated for each captured satellite, and positioning operation processing is performed using only observation data of captured satellites for which the pseudo-range residual is less than a predetermined level at the previous time.
  • the threshold for the pseudorange residual acts as a reference when selecting the capture satellites used for the positioning operation processing.
  • the threshold value for the pseudorange residual which is set in the assumed configuration, be larger in a multipath environment such as an urban area.
  • the threshold for the pseudorange residual for selecting the acquisition satellites used for the positioning operation processing is too small, the number of acquisition satellites usable for the positioning operation processing decreases, and the positioning accuracy is improved. It is because it becomes difficult to obtain an effect.
  • the quality of positioning signals from individual captured satellites is better than in a multipath environment, so even if the number of captured satellites used in positioning processing is relatively high, it is sufficiently high. Accuracy positioning is possible.
  • the threshold for the pseudorange residual is set to a small value, it is possible to suppress the possibility that observation data having a large error such as the pseudorange will be mixed in with the observation data used for the positioning operation processing, further improving the positioning accuracy. I can expect it. That is, under the open sky environment, it is preferable that the threshold for the pseudorange residual be smaller.
  • conditions for selecting and acquiring captured satellites (in other words, observation data) used for positioning calculation processing which should be set to obtain positioning results with higher accuracy, differ depending on the surrounding environment. .
  • the present disclosure has been made based on the above circumstances, and an object of the present disclosure is to achieve both improvement in positioning accuracy in an open sky environment and improvement in positioning accuracy in a multipath environment.
  • a positioning device is a transmission source of a received positioning signal based on a signal receiving unit that receives positioning signals transmitted from each of a plurality of positioning satellites, and a positioning signal received by the signal receiving unit.
  • the environment index value calculation unit sequentially calculates at least one of the amount of Doppler shift for the positioning satellite and the pseudorange as an index value for determining whether or not under multipath environment, and positioning by the signal reception unit
  • a scatter degree calculation unit that calculates, for each captured satellite that is a positioning satellite that can receive a signal, an index value scatter degree that indicates the degree of variation of the index value calculated by the environment index value calculation unit for the captured satellite;
  • a primary selection unit that determines, among the plurality of captured satellites, a non-use satellite that is a captured satellite that is not used for positioning calculation processing, based on the index value dispersion degree for each captured satellite calculated by the degree calculation unit;
  • Environment determines whether the surrounding environment is an environment where multipath is likely to occur, based on the deletion rate, which is the ratio of
  • the condition change unit dynamically changes the condition (hereinafter referred to as the adoption condition) of the captured satellite used for the positioning operation process according to the determination result of the environment determination unit. Therefore, for example, when it is determined by the environment determination unit that the surrounding environment is an environment in which multipath does not easily occur, the condition change unit changes the adoption condition to a relatively severe condition, thereby performing positioning operation processing.
  • the acquisition satellites used can be narrowed down to acquisition satellites with good signal quality to improve the positioning accuracy.
  • the adoption condition is relaxed to increase the number of captured satellites used for positioning calculation processing, thereby improving the positioning accuracy. It can also be illustrated.
  • condition change unit dynamically changes the adoption condition according to the surrounding environment, so that the improvement of the positioning accuracy under the open sky environment and the improvement of the positioning accuracy under the multipath environment Can be compatible.
  • FIG. 1 is a diagram showing an example of usage of a positioning device
  • FIG. 2 is a block diagram showing an example of a schematic configuration of a positioning device
  • FIG. 3 is a diagram for explaining the operation of the environment determination unit
  • FIG. 4 is a flowchart for explaining the secondary filtering process performed by the secondary filter unit.
  • FIG. 5 is a flowchart for explaining exception processing performed by the secondary filter unit.
  • FIG. 6 is a block diagram showing an example of a schematic configuration of a second-order filter unit of the first modification
  • FIG. 7 is a diagram showing an operation example of the threshold value changing unit of the first modification;
  • FIG. 1 is a diagram showing an example of usage of a positioning device
  • FIG. 2 is a block diagram showing an example of a schematic configuration of a positioning device
  • FIG. 3 is a diagram for explaining the operation of the environment determination unit
  • FIG. 4 is a flowchart for explaining the secondary filtering process performed by the secondary filter unit.
  • FIG. 5 is a flowchar
  • FIG. 8 is a block diagram showing an example of a schematic configuration of a positioning device in a second modification
  • FIG. 9 is a diagram for explaining the operation of the environment determination unit in the second modification
  • FIG. 10 is a diagram for explaining the operation of the environment determination unit in the second modification.
  • FIG. 1 is a diagram showing a schematic configuration of a satellite positioning system 100 including a positioning device (POSI APPA) 1 of the present embodiment.
  • the positioning device 1 is a device that sequentially receives current positioning signals by receiving positioning signals transmitted from the plurality of positioning satellites 2 and the quasi-zenith satellite 3, respectively.
  • the positioning device 1 is mounted on a vehicle Hv, and the positioning result (that is, current position information) of the positioning device 1 is used by any application software such as automatic driving and route guidance.
  • Vehicle Hv may be a four-wheeled vehicle, a two-wheeled vehicle, a three-wheeled vehicle or the like.
  • motorcycles also include motor bikes.
  • the vehicle Hv may be, for example, a construction vehicle, an agricultural work vehicle, or the like whose main purpose is other than the movement of the occupant.
  • the positioning device 1 may be mounted other than a vehicle.
  • the positioning device 1 may be mounted on a communication terminal (hereinafter referred to as a portable terminal) carried by a user, such as a smartphone, a tablet terminal, or a wearable terminal.
  • a communication terminal hereinafter referred to as a portable terminal
  • the positioning device 1 is used in any area on the ground. Here, as an example, it is assumed that the positioning device 1 is used in Japan. The details of the configuration, functions and the like of the positioning device 1 will be described later separately.
  • the positioning satellite 2 is an artificial satellite that constitutes a predetermined global positioning satellite system (hereinafter, GNSS: Global Navigation Satellite System).
  • GNSS Global Navigation Satellite System
  • GPS Global Navigation Satellite System
  • Galileo GPS
  • GLONASS Global Navigation Satellite System
  • BeiDou BeiDou
  • the positioning satellite 2 will be described as the artificial satellite forming the GPS (Global Positioning System) which is one of GNSS, but the configuration is not limited to this.
  • the positioning satellite 2 may be a satellite configuring GLONASS, or may be a satellite configuring BeiDou.
  • the satellite positioning system 100 may be a system combining multiple types of GNSS.
  • the positioning device 1 corresponds to each of a plurality of GNSSs such as GPS, Galileo, GLONASS, BeiDou, etc., and is configured to be able to execute positioning operation processing using a combination of positioning signals of a plurality of GNSSs. It may be done. For example, even if the positioning device 1 performs positioning operation processing using both the positioning signal transmitted from the positioning satellite 2 configuring GPS and the positioning signal transmitted from the positioning satellite 2 configuring GLONASS good.
  • positioning satellites are shown in FIG. 1, there may actually be more. For example, there are 30 or more positioning satellites 2 constituting the GPS. In addition, there are 20 or more positioning satellites 2 constituting the GLONASS.
  • Each positioning satellite 2 sequentially (for example, every 50 milliseconds) transmits a signal (hereinafter, positioning signal) obtained by phase-modulating data indicating transmission time and the like using a C / A code unique to each positioning satellite.
  • the positioning signal indicates, for example, data indicating an error of the satellite clock, ephemeris data indicating the current position of the satellite itself, almanac data indicating approximate orbit information of all the positioning satellites, etc., in addition to the transmission time.
  • Various data are transmitted sequentially. Since the C / A code is unique to each positioning satellite 2, the C / A code functions as information indicating the transmission source.
  • the positioning satellite 2 for which the positioning device 1 can receive the positioning signal is also referred to as a captured satellite.
  • Quasi-Zenith Satellite 3 is an artificial satellite orbiting an orbit (so-called quasi-zenith orbit) staying for a long time over a specific area including the area where positioning device 1 is used (that is, Japan).
  • QZSS Quasi-Zenith Satellite System
  • the QZSS corresponds to an example of a regional navigation satellite system (RNSS).
  • RNSS regional navigation satellite system
  • the quasi-zenith satellite 3 is an optional element.
  • the quasi-zenith satellite 3 is configured to be able to operate integrally with GPS, and sequentially transmits a positioning signal similar to the positioning signal transmitted by the positioning satellite 2 as a GPS satellite.
  • the quasi-zenith satellite 3 also includes the quasi-zenith satellite 3 in the positioning satellites 2 described below in order to provide the same function as the positioning satellite 2.
  • the positioning satellites of the present disclosure include quasi-zenith satellites 3 in addition to positioning satellites 2 as GPS satellites.
  • the quasi-zenith satellite 3 transmits a correction signal indicating correction information on a predetermined positioning satellite 2 using a satellite channel, in addition to the positioning signal.
  • the correction information here is, for example, information indicating precise satellite coordinates of the target positioning satellite 2, clock error, phase bias, ionospheric correction amount, tropospheric correction amount, and the like.
  • the ionospheric correction amount is a parameter that corrects the influence of the delay in the ionosphere, for example, an error of a pseudorange or carrier wave phase which occurs depending on the thickness of the ionosphere.
  • the troposphere correction amount is a parameter that corrects the influence of delay in the troposphere (for example, an error of pseudo range or carrier wave phase).
  • the correction signal transmitted by the quasi-zenith satellite 3 may be generated based on data generated and transmitted by a center (hereinafter, ground station) installed on the ground.
  • the correction signal may be the same signal as the SBAS augmentation message.
  • the positioning device 1 includes a carrier wave receiver 11, a reception processor 12, a dispersion degree calculator 13, an output unit 14, a primary filter unit 15, an environment determination unit 16, a secondary filter unit 17, and A positioning operation unit 18 is provided.
  • the reception processing unit 12 the dispersion degree calculation unit 13, the output unit 14, the primary filter unit 15, the environment determination unit 16, the secondary filter unit 17, and the positioning operation unit 18, a CPU (not shown) , Positioning program) can be realized.
  • the positioning device 1 includes a CPU that executes various arithmetic processing, a flash memory that is a non-volatile memory, a RAM that is a volatile memory, an I / O, and a bus line that connects these components. .
  • a program for causing a normal computer to function as the positioning device 1 may be stored in a non-transitory tangible storage medium.
  • the positioning program may be stored in a ROM, a flash memory or the like. Execution of the positioning program by the CPU corresponds to execution of a method corresponding to the positioning program.
  • the positioning device 1 receives the reception processing unit 12, the dispersion degree calculation unit 13, the output unit 14, the primary filter unit 15, and the environment determination unit as functional blocks that are expressed by execution of the positioning program by the CPU.
  • the second embodiment is provided with the second-order filter unit 17 and the positioning operation unit 18, the present invention is not limited to this.
  • Some or all of the functional blocks described above may be implemented as hardware.
  • the aspect realized as hardware also includes an aspect realized using one or more ICs and the like.
  • some or all of the above functional blocks may be realized by execution of software by the CPU and cooperation of hardware components.
  • the carrier wave receiver 11 is a communication module that receives a radio wave as a positioning signal transmitted from the positioning satellite 2, converts it into an electric signal, and outputs it.
  • the carrier wave receiver 11 includes an antenna (hereinafter referred to as a GNSS antenna) for receiving a positioning signal.
  • the carrier wave receiver 11 also includes a frequency conversion circuit for converting a received signal into a baseband signal, and an analog for converting an analog signal output from the GNSS antenna into a digital signal. It has a digital conversion unit and the like.
  • the carrier wave receiver 11 may include an amplifier circuit or the like.
  • the carrier wave receiver 11 corresponds to a signal receiver.
  • the carrier receiving unit 11 When the positioning device 1 is configured to be capable of performing positioning operation processing using positioning signals of a plurality of GNSSs, the carrier receiving unit 11 includes a plurality of GNSS antennas corresponding to the carrier frequencies of the respective GNSSs. It should just be.
  • the carrier wave receiver 11 may be configured to be able to receive the positioning signals of each positioning system to which the positioning device 1 should correspond.
  • the reception processing unit 12 is configured to generate observation data of the positioning satellite of the transmission source based on the positioning signal output from the carrier receiving unit 11.
  • the observation data is data on a captured satellite determined from the received positioning signal.
  • the reception processing unit 12 generates and outputs data indicating satellite number, observation time, Doppler shift amount, pseudo distance, satellite coordinates, carrier phase, and the like as observation data.
  • the satellite number is information for identifying each of the plurality of positioning satellites 2, and may be, for example, a PRN (Pseudo Random Noise) ID.
  • the observation time may be the time when the positioning signal for generating the observation data is received.
  • the amount of Doppler shift is a parameter that represents the difference between the carrier frequency and the reception frequency caused by the Doppler effect.
  • the pseudo distance is a distance determined from the time from when the positioning signal is transmitted from the positioning satellite 2 to when it is received by the reception processing unit 12 (hereinafter, signal propagation time) and the propagation speed of the radio wave. Note that the signal propagation time may be calculated from the difference between the transmission time and the reception time, or may be calculated based on the phase shift amount of the C / A code.
  • the satellite coordinates are information indicating the current position of the positioning satellite on the satellite orbit.
  • the reception processing unit 12 includes a pseudo distance calculation unit that calculates a pseudo distance, a Doppler shift amount calculation unit, and the like as sub functions for generating the observation data described above.
  • simulation distance calculation part and a Doppler shift amount calculation part is abbreviate
  • the reception processing unit 12 having the functions as the pseudo distance calculation unit and the Doppler shift amount calculation unit corresponds to the environment index value calculation unit.
  • the reception processing unit 12 generates the observation data described above each time the carrier receiving unit 11 receives a positioning signal. Thereby, observation data for each positioning satellite (that is, captured satellite) that can receive the positioning signal is sequentially generated. Then, the observation data for each captured satellite generated by the reception processing unit 12 is output to the output unit 14. Further, the pseudo distance generated by the reception processing unit 12 is associated with information (for example, a satellite number) indicating which of the captured satellites is a pseudo distance, and is output to the dispersion degree calculation unit 13.
  • the dispersion degree calculation unit 13 acquires pseudo distances for each captured satellite sequentially provided from the reception processing unit 12 and stores the pseudo distances separately for each captured satellite in a memory (not shown).
  • the memory may be realized using a rewritable storage medium such as a RAM.
  • a plurality of pseudoranges at different acquisition times for the same acquisition satellite may be sorted in chronological order and stored in memory, for example, so that the latest data is at the top. Also, data that has been stored for a certain period of time may be sequentially discarded.
  • the dispersion degree calculation unit 13 calculates, for each captured satellite, a pseudo distance dispersion degree that indicates the degree of variation of the pseudo distance calculated (in other words, observed) within the last predetermined time.
  • the pseudo distance scatter degree for a certain acquisition satellite may be a standard deviation with the pseudo distance calculated within the past constant time (for example, 5 seconds) from the current time as a population.
  • the dispersion degree may adopt another index indicating the degree of dispersion of data in statistics, such as dispersion or the difference between the maximum value and the minimum value.
  • the pseudo range scatter degree for a certain capture satellite may be calculated when a new pseudo range for the capture satellite is obtained.
  • the scattering degree calculation unit 13 may be configured to calculate the pseudo distance scattering degree for each captured satellite sequentially (for example, every 100 milliseconds) at a predetermined cycle.
  • the scattering degree calculation unit 13 is configured to calculate a shift amount scattering degree indicating the degree of variation of the Doppler shift amount calculated and observed within the latest predetermined time, instead of the pseudo distance scattering degree May be
  • the reception processing unit 12 may be configured to output the Doppler shift amount to the distribution degree calculation unit 13 instead of the pseudo distance.
  • the distribution degree calculation unit 13 may be configured to calculate both the pseudo distance distribution degree and the shift amount distribution degree for each captured satellite.
  • the reception processing unit 12 may be configured to sequentially output both the pseudo distance and the Doppler shift amount to the dispersion degree calculation unit 13.
  • the dispersion of the pseudo distance which is sequentially calculated becomes large, so the pseudo distance is used as an index value of whether or not the surrounding environment is under the multipath environment.
  • the arrival direction of the positioning signal may change successively. If the direction of arrival changes, the amount of Doppler shift also changes.
  • the variation in the Doppler shift amount that is sequentially calculated becomes large, so the Doppler shift amount that is sequentially calculated is also the multipath environment in the surrounding environment. It functions as an indicator value of whether or not there is any. That is, the pseudo distance and the Doppler shift amount correspond to the index value. Further, the pseudo distance distribution degree and the shift amount distribution degree correspond to the index value distribution degree.
  • the pseudo distance scatter degree calculated by the scatter degree calculation unit 13 is output to the output unit 14 in association with information (for example, a satellite number) indicating which of the captured satellites the pseudo distance scatter degree is.
  • the output unit 14 outputs, to the primary filter unit 15, the observation data for each captured satellite input from the reception processing unit 12 and the pseudo distance distribution degree for each captured satellite input from the distribution calculation unit 13.
  • the primary filter unit 15 deletes observation data of captured satellites having a pseudo distance distribution degree equal to or greater than a predetermined threshold (hereinafter, a dispersion degree threshold), among observation data of each captured satellite. Then, the primary filter unit 15 outputs, to the secondary filter unit 17, observation data remaining after deleting observation data of captured satellites for which the pseudo distance distribution degree is equal to or greater than the distribution degree threshold value. That is, the primary filter unit 15 uses the observation data for each captured satellite so that only observation data for captured satellites having a relatively small degree of pseudo-dispersion are used in the secondary filter unit 17 that is the configuration of the subsequent stage. It is the structure to filter.
  • a dispersion degree threshold a predetermined threshold
  • the primary filter unit 15 corresponds to a primary selection unit.
  • a captured satellite corresponding to observation data remaining as a result of the primary filtering process is also referred to as a primary residual satellite.
  • the primary filter unit 15 corresponds to a configuration for determining, from among the plurality of acquired satellites, non-use satellites that are acquired satellites that are not used for the positioning operation processing based on the pseudorange scatter degree for each acquired satellite. In other words, it corresponds to a configuration in which the acquisition satellites to be used for the positioning operation processing are selected from among the plurality of acquisition satellites based on the pseudo range scatter degree for each acquisition satellite.
  • the primary filter unit 15 calculates a deletion rate R that represents the ratio of the number of observation data (in other words, the total number of captured satellites) input from the output unit 14 and the number of observation data deleted.
  • the deletion rate R is represented by k / N.
  • the number M of observation data output to the secondary filter unit 17 is N ⁇ k. .
  • the primary filter unit 15 outputs data indicating the deletion rate R to the environment determination unit 16.
  • the environment determination unit 16 is based on the deletion rate R provided from the primary filter unit 15 in an environment where the surrounding environment of the positioning device 1 has many features that block and scatter signals from positioning satellites in the sky view (in other words, Then, it is determined whether or not the environment is likely to cause multipath. Whether multipath is likely to occur is determined by the good view of the sky above the positioning device 1.
  • the surrounding environment is classified into four categories in stages according to the goodness of the sky, and the environment determination unit 16 selects the current surrounding environment according to the deletion rate R. It shall be judged which category it corresponds to.
  • the four categories will be referred to as a first category C1, a second category C2, a third category C3 and a fourth category C4 in the order of good view of the sky.
  • the first category C1 has the best view over the sky
  • the fourth category C4 corresponds to the environment having the worst view over the sky.
  • the environment determination unit 16 determines that the surrounding environment corresponds to the category having the larger category number as the deletion rate R is larger.
  • the environment determining unit 16 determines that the surrounding environment corresponds to the first category C1.
  • a specific value of the first deletion rate Rth1 may be set appropriately by conducting a test. Here, 10% is taken as an example, but it may be 5% or 15%.
  • the deletion rate R is equal to or higher than the first deletion rate Rth1 and smaller than the predetermined second deletion rate Rth2, it is determined that the surrounding environment corresponds to the category C2.
  • a specific value of the second deletion rate Rth2 may be appropriately set by conducting a test in a range larger than the first deletion rate Rth1, and is 30% here as an example.
  • the environment determination unit 16 determines that the surrounding environment corresponds to the third category C3.
  • a specific value of the third deletion rate Rth3 may be appropriately set by performing a test larger than the second deletion rate Rth2, and is 60% as an example here. Then, when the deletion rate R is equal to or higher than the third deletion rate Rth3, it is determined that the surrounding environment corresponds to the fourth category C4.
  • the environment determination unit 16 determines the surrounding environment each time the deletion rate R is input from the primary filter unit 15. Naturally, the surrounding environment changes as the vehicle Hv travels, and the deletion rate R changes. That is, as shown in FIG. 3, as the deletion rate R changes, the determination result of the environment determination unit 16 also changes with time. The determination result of the environment determination unit 16 is sequentially output to the secondary filter unit 17.
  • the first category C1 is a category corresponding to a so-called open sky environment.
  • the sky ratio in the case of photographing the sky with a fisheye camera may be set to be 80% or more.
  • an environment in which the sky ratio is 75% or 85% or more may be defined as an open sky environment.
  • the first deletion rate Rth1 described above may be determined based on the test result of the deletion rate R in the open sky environment.
  • the open sky environment is an environment in which relatively accurate positioning results can be obtained.
  • the first category C1 corresponds to the open sky category.
  • the second category C2 is a category that indicates an area having relatively few buildings such as buildings in an urban area (in other words, an area in which multipath is less likely to occur).
  • the second category C2 may correspond to an environment in which the sky ratio is relatively low (for example, 65% or more and less than 80%).
  • the third category C3 is a category that indicates an area having a relatively large number of structures such as buildings in an urban area (in other words, an area where multipaths are likely to occur).
  • the criteria for separating the second category C2 and the third category C3 may be appropriately determined based on the test results.
  • the third category C3 may correspond to an environment in which the sky ratio is relatively low (for example, 50% or more and less than 65%).
  • the fourth category C4 is a category indicating an urban area where high-rise buildings and the like exist.
  • the fourth category C4 may be a category that assumes an environment in which the sky ratio is less than 50%.
  • the fourth category C4 corresponds to the urban area category.
  • the surrounding environment is classified into four stages of categories C1 to C4 according to the sky view, but the present invention is not limited to this.
  • the surrounding environment may be classified into two stages, three stages, or five stages. For example, it may be classified into two stages of whether or not it is an open sky environment.
  • the peripheral environment may be divided into three categories as a whole. Environment such as in the forest or inside the tunnel may be added as a category.
  • the secondary filter unit 17 removes observation data that is relatively strongly affected by multipath etc. from observation data for each primary residual satellite provided from the primary filter unit 15, and the remaining observation data Are output to the positioning operation unit 18. That is, the secondary filter unit 17 corresponds to a configuration for extracting observation data to be finally used for the positioning operation processing from the observation data for each primary residual satellite.
  • the secondary filter unit 17 includes a residual calculation unit 171 and a threshold change unit 172 as sub functions.
  • the secondary filter unit 17 including the residual calculating unit 171 and the threshold changing unit 172 operates roughly as follows.
  • the secondary filter unit 17 first calculates the pseudorange residual ⁇ D for each of the primary residual satellites extracted as a result of the primary filtering process by the primary filter unit 15. Then, only the observation data of the captured satellite whose pseudorange residual ⁇ D is less than the predetermined residual threshold Dth is output to the positioning operation unit 18.
  • the residual calculation unit 171 described above is configured to calculate the pseudorange residual ⁇ D for each primary residual satellite.
  • the pseudorange residual ⁇ D for a certain positioning satellite 2 is the linear distance from the current position of the positioning device 1 to the positioning satellite 2 and the time difference from the transmission of the positioning signal to the reception by the positioning device 1 and the radio wave
  • the absolute value of the difference with the distance (ie, pseudo distance) calculated based on the propagation velocity of The linear distance from the current positioning device 1 to the positioning satellite 2 can be calculated from the current position coordinates of the positioning device 1 and the position coordinates of the positioning satellite 2.
  • the pseudorange residual ⁇ D calculated as described above tends to take a large value when the positioning signal is strongly affected by multipath and the like and the pseudorange includes an error. That is, the pseudorange residual ⁇ D for a certain acquisition satellite functions as a parameter representing the possibility that the pseudorange calculated based on the positioning signal transmitted from the acquisition satellite contains an error. Therefore, the pseudorange residual ⁇ D corresponds to the error index value. Further, the residual calculation unit 171 corresponds to an error index value calculation unit.
  • the position coordinates as the current position of the positioning device 1 used when calculating the pseudo distance residual ⁇ D may be strictly the position coordinates calculated by the positioning operation unit 18 at the previous time. Since the positioning operation unit 18 basically performs the positioning operation process approximately every 50 to 100 milliseconds, the position coordinates calculated by the positioning operation unit 18 at the previous time are positioning 50 to 100 milliseconds before the present time It shows the position of the device 1. In calculating the pseudo distance residual ⁇ D, it is assumed that an error between the position coordinates calculated by the positioning operation unit 18 at the previous time and the actual current position can be ignored. Of course, by performing known dead reckoning based on the position coordinates calculated by the positioning operation unit 18 at the previous time, the current position may be estimated more strictly and the pseudo distance residual ⁇ D may be calculated.
  • secondary filtering processing a process of deleting observation data of a captured satellite whose pseudorange residual ⁇ D is equal to or more than a predetermined residual threshold among observation data for each primary residual satellite.
  • secondary residual satellite a captured satellite corresponding to observation data remaining as a result of the secondary filtering process.
  • the residual threshold Dth used in the secondary filtering process functions as a threshold for the pseudorange residual for selecting and selecting a captured satellite (in other words, observation data) used in the positioning operation process. Therefore, the residual threshold Dth corresponds to the filtering threshold.
  • the threshold changing unit 172 included in the secondary filter unit 17 is configured to change the value of the residual threshold Dth based on the category of the surrounding environment determined by the environment determining unit 16.
  • the threshold changing unit 172 sets the residual threshold Dth to a relatively small value when the environment determining unit 16 determines that the surrounding environment is an environment in which multipath is less likely to occur (that is, a category with a smaller number). Do. When it is determined by the environment determination unit 16 that the surrounding environment is an environment in which multipath is likely to occur (that is, a category with a large number), the residual threshold Dth is set to a relatively large value.
  • a plurality of values (in other words, candidate values) that can be set by the threshold value changing unit 172 as the residual threshold value Dth are set in advance. The multiple levels of candidate values are associated with each category of multiple levels. The details of the operation of the threshold value changing unit 172 will be described later with reference to FIG.
  • the residual threshold Dth is set to a relatively large value, the primary residual satellites are easily selected as positioning satellites to be used for the positioning operation processing, and if the residual threshold Dth is set to a relatively small value, The primary residual satellites are less likely to be selected as positioning satellites for use in the positioning operation process. That is, the residual threshold Dth functions as a parameter that defines the conditions of the captured satellites used for the positioning operation process.
  • the threshold value changing unit 172 is configured to loosen or tighten the conditions of the captured satellites used for the positioning calculation process based on the determination result of the environment determination unit 16. Therefore, the threshold change unit 172 corresponds to a condition change unit.
  • the secondary residual satellite described above corresponds to a captured satellite that satisfies the conditions of the captured satellite used for the positioning operation process set by the threshold value changing unit 172.
  • the positioning operation unit 18 uses the observation data provided by the secondary filter unit 17 to carry out a positioning operation process. That is, positioning operation processing is performed using observation data on secondary residual satellites remaining as a result of secondary filtering processing by the secondary filter unit 17.
  • PPP-AR more specifically, MADOCA-PPP
  • various known methods such as PPP-RTK can be adopted as a method of positioning operation processing (hereinafter, positioning method) itself. Since the specific procedure of each positioning system is well-known, it abbreviate
  • PPP is an abbreviation for Precise Point Positioning (that is, single carrier phase positioning)
  • AR is an abbreviation for Ambiguity Resolution
  • MADOCA is an abbreviation of Multi-gnss Advanced Demonstration tool for Orbit and Clock Analysis
  • RTK is an abbreviation of Real-Time Kinematic.
  • the positioning operation unit 18 outputs position information as a result of the positioning operation processing to a predetermined application.
  • An application that uses position information may be incorporated in the positioning device 1 or may be incorporated in an electronic device provided outside the positioning device 1.
  • the positioning operation unit 18 also outputs the result of the positioning operation process to the secondary filter unit 17.
  • the position information provided to the secondary filter unit 17 is used to calculate the next pseudo distance residual ⁇ D.
  • step S101 based on observation data (mainly pseudo distance and satellite coordinates) for each primary residual satellite and position information of the positioning device 1, the residual calculation unit 171 generates pseudo distance residuals for each primary residual satellite. After ⁇ D is calculated, the process proceeds to step S101.
  • step S102 the threshold changing unit 172 determines whether the determination result of the surrounding environment by the environment determining unit 16 is the first category C1.
  • the determination in step S102 is affirmative and the step S103 is executed.
  • a negative determination is made in step S102, and step S104 is executed.
  • the threshold changing unit 172 sets the residual threshold Dth to a predetermined first threshold Dth1 prepared in advance as a candidate value for the residual threshold Dth, and then proceeds to step S108.
  • the specific value of the first threshold Dth1 may be determined based on the test result of the pseudorange residual under the open sky environment. For example, the first threshold Dth1 may be set to a value of about 2 to 5 m.
  • the first threshold Dth1 corresponds to the candidate value associated with the first category C1 as the open sky category.
  • step S104 it is determined whether the determination result of the surrounding environment by the environment determination unit 16 is the second category C2.
  • the determination in step S104 is affirmative and the step S105 is executed.
  • a negative determination is made in step S104 and step S106 is performed.
  • the threshold changing unit 172 sets the residual threshold Dth to a predetermined second threshold Dth2, and proceeds to step S108.
  • a specific value of the second threshold Dth2 may be appropriately set in a range larger than the first threshold Dth1.
  • the second threshold Dth2 increases the number of secondary residual satellites by one or more compared to when the residual threshold Dth is set to the first threshold Dth1, and the number of positioning satellites used for the positioning operation processing
  • the value is set to a value at which an improvement in positioning accuracy can be obtained by the increase in the number of satellites used.
  • the specific value for obtaining the above effect may be determined based on the test result in the second category C2.
  • the second threshold Dth2 corresponds to the candidate value associated with the second category C2.
  • step S106 it is determined whether the determination result of the surrounding environment by the environment determination unit 16 is the third category C3.
  • the determination in step S106 is positive and the step S107 is executed.
  • the determination in step S106 is negative and the step S108 is executed.
  • the case where the determination result of the surrounding environment by the environment determining unit 16 is not the third category C3 corresponds to the case where the determination result of the surrounding environment by the environment determining unit 16 is the fourth category C4.
  • the threshold value changing unit 172 sets the residual threshold value Dth to a predetermined third threshold value Dth3, and then proceeds to step S108.
  • a specific value of the third threshold Dth3 may be appropriately set in a range larger than the second threshold Dth2.
  • the third threshold Dth3 increases the number of secondary residual satellites by one or more compared to when the residual threshold Dth is set to the second threshold Dth2, and the effect of improving the positioning accuracy by increasing the number of satellites used is It is preferable to set to the value of the level obtained.
  • the specific value for obtaining the above effect may be determined based on the test result in the third category C3.
  • the third threshold Dth3 corresponds to the candidate value associated with the third category C3.
  • the threshold value changing unit 172 sets the residual threshold value Dth to a predetermined fourth threshold value Dth4, and then proceeds to step S108.
  • a specific value of the fourth threshold Dth4 may be appropriately set in a range larger than the third threshold Dth3.
  • the fourth threshold Dth4 increases the number of secondary remaining satellites by 1 to 2 or more as compared to the case where the residual threshold Dth is set to the third threshold Dth3, and the effect of improving the positioning accuracy by increasing the number of satellites used is It is preferable to set to the value of the level obtained.
  • the specific value for obtaining the above effect may be determined based on the test result in the fourth category C4.
  • the fourth threshold Dth4 which corresponds to the maximum value of candidate values that can be set as the residual threshold, is set to a sufficiently large value equivalent to infinity so that all primary residual satellites can be used for positioning processing. It may be set.
  • the fourth threshold value Dth4 corresponds to a candidate value associated with the fourth category C4 as the urban area tegor.
  • step S109 among the observation data for each captured satellite input from the primary filter unit 15, the secondary filter unit 17 has the pseudo distance residual ⁇ D smaller than the residual threshold Dth set in the above processing.
  • the observation data of the capture satellite is extracted, and the process proceeds to step S110.
  • step S110 the secondary filter unit 17 determines whether the positioning operation unit 18 can execute the positioning operation processing based on the observation data extracted in step S109. For example, when the number of observation data extracted in step S109 (in other words, the number of secondary remaining satellites) is four or more, the positioning operation unit 18 determines that the positioning operation is possible. On the other hand, when the number of observation data extracted in step S109 is less than 4, it is determined that the positioning operation unit 18 can not perform the positioning operation.
  • step S110 When the number of observation data extracted in step S109 is sufficiently large that the positioning operation unit 18 can perform positioning operation, the determination in step S110 is positive and step S111 is executed. On the other hand, when the number of observation data extracted in step S109 is so small that the positioning operation unit 18 can not perform positioning operation, the determination in step S110 is negative and step S112 is executed.
  • the number of observation data to be used as a determination criterion of whether or not the positioning operation unit 18 can perform positioning operation is 3 It may be
  • step S111 the secondary filter unit 17 outputs the observation data extracted in step S109 to the positioning operation processing, and the present flow ends.
  • the positioning operation unit 18 performs positioning operation processing using the observation data for each captured satellite provided from the secondary filter unit 17, and identifies the current position. Then, the position information as the positioning result is output to the predetermined application and the secondary filter unit 17.
  • step S112 the secondary filter unit 17 performs predetermined exception processing.
  • the exception processing is processing in the case where the positioning operation unit 18 can not carry out the positioning operation processing on the observation data extracted by the series of processing in S101 to S109. This exception processing will be described using the flowchart shown in FIG. The flowchart illustrated in FIG. 5 is executed when the process proceeds to step S112 in FIG.
  • the threshold value changing unit 172 sets the value of the residual threshold value Dth to a value one step larger than the current value, and proceeds to step S202. Specifically, when the current residual threshold Dth is set to the first threshold Dth1, the residual threshold Dth is set to the second threshold Dth2. When the current residual threshold Dth is set to the second threshold Dth2, the residual threshold Dth is set to the third threshold Dth3. Similarly, when the current residual threshold Dth is set to the third threshold Dth3, the residual threshold Dth is set to the fourth threshold Dth4.
  • step S202 of the observation data for each captured satellite input from the primary filter unit 15, the observation about the captured satellite whose pseudorange residual ⁇ D is less than the residual threshold Dth newly set in step S201. Data is extracted and it moves to step S203.
  • step S203 it is determined based on the observation data extracted in step S202 whether the positioning operation unit 18 can execute the positioning operation process.
  • the substantive process contents of step S203 are the same as those of step S110 described above, and therefore the description thereof is omitted.
  • the determination in step S203 is positive and step S204 is executed.
  • the determination in step S203 is negative and the process returns to step S201.
  • step S204 after prohibiting the positioning operation unit 18 from outputting the result of the positioning operation process to the external application, the observation data extracted in step S202 is output, and the process proceeds to step S205.
  • the positioning operation unit 18 performs a positioning operation process using the observation data provided from the secondary filter unit 17, and returns the result to the secondary filter unit 17. Further, based on the instruction from the secondary filter unit 17, the positioning result is not output to an external application. That is, when the positioning calculation processing can not be performed using only observation data that satisfies the residual threshold value Dth according to the determination result of the environment determination unit 16, the positioning result is discarded externally.
  • step S205 the result of the positioning operation process in the positioning operation unit 18 is acquired, and the present flow ends.
  • the position information acquired in step S205 is used to calculate the next pseudo distance residual.
  • a positioning device that does not have a configuration corresponding to the primary filter unit 15 and selects and uses captured satellites used for positioning operation processing using a constant residual threshold regardless of the surrounding environment To introduce the effects of the present embodiment.
  • the quality of positioning signals from individual captured satellites tends to be better than in multi-path environments such as urban areas, so the number of observation data (in other words, captured satellites) used You can get enough information even if you That is, in the open sky environment, sufficiently high precision positioning can be performed even if the number of captured satellites used for the positioning operation processing is relatively small.
  • the positioning accuracy does not necessarily improve. If the residual threshold value is set to a relatively large value in the first comparison configuration, there is a high possibility that observation data having a large error such as pseudo distance may be mixed with observation data used in the positioning operation process. There is a concern that the accuracy may deteriorate.
  • the residual threshold be large to some extent. If the residual threshold value is set to a relatively small value in the first comparison configuration, the number of captured satellites used for the positioning operation processing is reduced, and the effect of improving the accuracy by increasing the number of satellites used is obtained. Because it is difficult to
  • the threshold changing unit 172 determines that the environment is the environment where multipath is less likely to occur (that is, a category with a smaller number)
  • the residual is left
  • the threshold Dth is set to a relatively small value.
  • the threshold changing unit 172 has a relatively large value for the residual threshold Dth.
  • the adoption condition is eased, and the number of captured satellites used for the positioning operation processing can be increased.
  • the number of captured satellites used for positioning operation processing is high, the positioning result tends to be more accurate. Therefore, the number of captured satellites used for positioning operation processing is increased by the above processing. Higher precision can be achieved.
  • the positioning device 1 of the present embodiment can be expected to further improve the positioning accuracy by providing the primary filter unit 15. Specifically, it is as follows. In a multipath environment, the fluctuation of the pseudorange is relatively large. Naturally, the pseudorange residuals are also highly volatile because the pseudoranges of captured satellites that are strongly affected by multipath are highly volatile. That is, even if the pseudorange residual of a certain acquisition satellite is small at one time, it is not necessarily small at the next time. Therefore, in the first comparison configuration, there remains a possibility that observation data having a large error such as pseudo distance may be mixed in with the observation data used for the positioning operation process.
  • the primary filter unit 15 has only the observation data of captured satellites having a relatively small degree of pseudo-distance distribution among the observation data of each captured satellite as the configuration of the subsequent stage. It is output to the secondary filter unit 17. Therefore, since observation data on captured satellites with large fluctuations in pseudorange are removed by primary filtering processing, it is possible to suppress the possibility that observation data with large errors such as pseudorange will be mixed in with observation data used for positioning operation processing. be able to. As described above, according to the configuration of the present embodiment, it is possible to cope with the case where the pseudorange residual instantaneously falls below the residual threshold Dth, and further improvement of the positioning accuracy can be expected.
  • the category of the surrounding environment is determined based on the deletion rate R determined from the reception status of the positioning signal.
  • second comparison configuration a configuration using map data may be considered.
  • the positioning device needs to have a function of acquiring map data by communicating with a storage device in which map data is stored or an external server, and the cost for that amount. Will increase.
  • the environment determination unit 16 determines the category of the surrounding environment based on the deletion rate R determined from the reception status of the positioning signal. It does not have to have the ability to communicate with it. That is, according to the present embodiment, the manufacturing cost can be suppressed more than in the second comparison configuration.
  • the present invention can also be applied to devices that do not have a wide area communication module.
  • the first case when five observation data are deleted when the number of acquisition satellites is 10 (hereinafter, the first case) and when five observation data are deleted when the number of acquisition satellites is 30, (hereinafter, the second case)
  • the first case suggests that it is an environment in which multipath is more likely to occur (in other words, an environment with many structures around it).
  • the surrounding environment there is a relatively high possibility that the surrounding environment may be erroneously determined.
  • the environment determination unit 16 determines the category of the surrounding environment not by the number of deleted observation data but by the deleted ratio (that is, the deletion rate R) as in the configuration of the present embodiment, the surrounding environment is erroneous. It is possible to suppress the risk of judging. In other words, by using the deletion rate R instead of the number of deletions of observation data as an index for determining the surrounding environment, it is possible to more accurately determine the surrounding environment.
  • the threshold changing unit 172 may determine the residual threshold Dth depending on whether or not the quasi-zenith satellite 3 has been captured, in addition to the determination result of the environment determining unit 16. For example, when the quasi-zenith satellite 3 can be captured, the residual threshold Dth may be set to a smaller value than when the quasi-zenith satellite 3 can not be captured.
  • a configuration based on the above-described idea will be described as a first modification.
  • the secondary filter unit 17 in the modification 1 includes a QZS capture determination unit 173 in addition to the residual calculation unit 171 and the threshold change unit 172.
  • the QZS acquisition determination unit 173 is configured to determine whether the positioning signal from the quasi-zenith satellite 3 can be received (in other words, whether or not acquisition can be performed). Whether or not the quasi-zenith satellite 3 can be captured can be determined by the elevation angle of the captured satellite or the transmission source information indicated by the positioning signal.
  • the QZS capture determination unit 173 may determine that the quasi-zenith satellite 3 can be captured.
  • the quasi-zenith angle is an angle for separating the positioning satellite 2 (for example, a GPS satellite) taking an orbit that orbits the earth once, and the quasi-zenith satellite 3.
  • the quasi-zenith angle may be set to the lower limit (for example, 70 °) of the angle that can be regarded as the zenith direction.
  • the elevation angle is an angle formed by a straight line connecting the captured satellite and the positioning device 1 with respect to the ground plane relative to the current position of the positioning device 1. The elevation angle can be calculated based on the position coordinates of each of the positioning device 1 and the captured satellite.
  • the QZS capture determination unit 173 may determine based on transmission source information (for example, PRN-ID) indicated by the positioning signal. Furthermore, when the correction signal transmitted by the quasi-zenith satellite 3 can be received, it may be determined that the quasi-zenith satellite 3 can be captured.
  • the QZS capture determination unit 173 corresponds to a quasi-zenith capture determination unit.
  • the threshold change unit 172 of the modification 1 sets the residual threshold Dth to the category determined by the environment determination unit 16. Set to the associated value.
  • the residual threshold Dth is higher than the value associated with the category determined by the environment determination unit 16. Set to a small value. For example, as shown in FIG. 7, when the QZS capture determination unit 173 determines that the quasi-zenith satellite 3 can be captured, it is determined that the residual threshold Dth can not be captured. In this case, the value to be applied may be set to 0.8 times the value to be applied.
  • the magnification here is not limited to 0.8, but may be 0.9, 0.7, or the like.
  • the residual threshold Dth may be reduced by subtracting a predetermined value (e.g., 0.5 m) instead of reducing the residual threshold Dth by multiplying the coefficient by a predetermined magnification.
  • FIG. 7 discloses an aspect in which the residual threshold Dth when the quasi-zenith satellite 3 can be captured is set to a value smaller than the value when the quasi-zenith satellite 3 is not captured in all categories. But it is not limited to this.
  • the residual threshold Dth when the quasi-zenith satellite 3 can be captured only when the category determined by the environment determination unit 16 corresponds to C1 to C2, when the quasi-zenith satellite 3 can not be captured It may be set to a value smaller than the value.
  • the quasi-zenith satellite 3 When the quasi-zenith satellite 3 can be captured, the pseudoranges for each positioning signal are made more accurate by the ionospheric correction amount and the troposphere correction amount provided from the quasi-zenith satellite 3. Hence, the pseudorange residual itself of each acquisition satellite can be generally small. As described above, since the average level of the pseudorange residual differs depending on whether or not the quasi-zenith satellite 3 has been acquired, the same value is temporarily set in the case where the quasi-zenith satellite 3 is acquired and not acquired. If applied, if the quasi-zenith satellite 3 can be captured, the possibility that the observation data having a large error derived from the multipaths may be mixed with the observation data used for the positioning operation processing is increased.
  • the residual threshold Dth is set to a smaller value than in the case where the quasi-zenith satellite 3 is not captured. It is possible to further suppress the possibility that the observation data having a large error may be mixed in with the observation data used for the positioning operation process. As a result, further improvement of the positioning accuracy can be expected.
  • the QZS acquisition determination unit 173 determines whether or not the positioning signal from the quasi-zenith satellite 3 can be received, but the present invention is not limited to this.
  • the QZS acquisition determination unit 173 may determine whether or not the reinforcement signal from the quasi-zenith satellite 3 can be received.
  • the threshold changing unit 172 performs local correction information.
  • the same processing as that of the above-described first modification may be performed depending on whether or not acquisition is possible. That is, when the local correction information can be acquired, the residual threshold Dth is set to a smaller value than when the local correction information can not be acquired. According to such a configuration, the same effect as that of the first modification can be obtained.
  • short-range communication direct wireless communication
  • the determination result of the surrounding environment in another device is acquired via the short-range communication unit 4, and the current surrounding environment category It may be determined finally.
  • the positioning device 1 can correct the error when only the device itself misjudges the category of the surrounding environment for some reason. That is, the accuracy and reliability of the category finally used can be increased.
  • the above-mentioned narrow area communication is communication which does not go through a wide area communication network, in other words.
  • the short-range communication unit 4 performs communication directly (that is, inter-vehicle communication) with other vehicles existing around the vehicle Hv.
  • the frequency band used for inter-vehicle communication is, for example, the 760 MHz band.
  • 2.4 GHz and 5.9 GHz bands can also be used.
  • a communication standard for realizing inter-vehicle communication can adopt any one.
  • the standard of WAVE (Wireless Access in Vehicular Environment) disclosed in IEEE 1609 or the like can be adopted.
  • the short-range communication unit 4 as a communication module for inter-vehicle communication is configured to be communicable with other vehicles within a radius of several hundred meters centering on the host vehicle Hv.
  • the short-range communication part 4 is set as a communication module which provides the wireless communication function based on the inter-vehicle communication specification which becomes about several hundreds m communication range as an example here, it is not restricted to this.
  • the short range communication unit 4 is a communication module for performing communication (hereinafter referred to as short distance communication) conforming to a predetermined short distance wireless communication standard in which the communication range is, for example, several tens of meters at the maximum. It is good.
  • Bluetooth Low Energy Bluetooth is a registered trademark
  • Wi-Fi registered trademark
  • ZigBee registered trademark
  • the short range communication unit 4 determines a determination result packet indicating category information (for example, category number) indicating a category determined by the environment determination unit 16 of the own device based on a request from the positioning device 1 (in particular, the environment determination unit 16). Broadcast. Also, when the determination result packet transmitted from the other vehicle is received, the data indicated in the determination result packet is provided to the positioning device 1. As a premise, the environment determination unit 16 in the present modification 2 generates a determination result packet each time the determination of the category of the surrounding environment based on the deletion rate input from the primary filter unit 15 is performed, It is output to the communication unit 4 and broadcasted.
  • category information for example, category number
  • the determination result packet is a communication packet including information (for example, a category number) indicating the determination result of the surrounding environment and position information as the calculation result of the positioning operation unit 18.
  • information for example, a category number
  • the determination result packet includes information indicating the transmission time of the communication packet, the traveling speed, and the traveling direction.
  • Environment determination unit 16 in the present modification determines the determination results of other vehicles (hereinafter, peripheral vehicles) present within predetermined distance D1 from own vehicle Hv. Obtain one by one.
  • a majority determination unit 161 is provided which finally determines to which category the surrounding environment belongs by majority determination with the determination results of the surrounding vehicles and the host vehicle Hv as the population.
  • the majority decision unit 161 determines that the surrounding environment is the first category C1 as a result of the determination of the surrounding vehicle in the situation where it is determined to be the second category C2 itself. If it is a group, it is finally determined that the surrounding environment is the first category C1. That is, it discards its own determination result that the surrounding environment is the second category C2, and determines that the category of the surrounding environment is the first category.
  • the distance D1 is preferably set to a sufficiently small value such as, for example, about several tens of meters (for example, 15 m) as compared with the distance in which inter-vehicle communication can be performed. This is because the surrounding environment for the other vehicle traveling in a place away from the host vehicle Hv is likely to be different from the surrounding environment of the host vehicle Hv.
  • the other vehicle corresponding to the surrounding vehicle may be extracted using the current position of the host vehicle Hv and the current position of the transmission source vehicle.
  • a transmission source vehicle whose current position is within a distance D1 from the host vehicle Hv may be adopted as a surrounding vehicle.
  • the positioning device 1 mounted on the host vehicle Hv corresponds to the own device, and the positioning device 1 mounted on the nearby vehicle corresponds to the peripheral device.
  • the positioning device 1 corrects the error when only the own device misjudges the category of the surrounding environment for some reason. be able to. As a result, the accuracy and reliability of the finally used category can be increased.
  • the mode category which is the determination result with the largest number of appearances (in other words, the majority) among the determination results of the surrounding vehicles and the determination result of the host vehicle Hv, exceeds the majority of the whole.
  • the peripheral environment may adopt the majority judgment result only when it is present.
  • the fact that the most frequent category does not exceed a majority means that the judgment result for each vehicle (more specifically, the positioning device 1) is broken, and the majority judgment result is not necessarily the correct judgment result. It is because there is not. That is, by adopting the majority judgment result only when the majority judgment result exceeds the majority, it is possible to improve the judgment accuracy of the surrounding environment. In addition, despite the fact that oneself has correctly determined the surrounding environment, it is also possible to reduce the possibility that an erroneous determination result will be adopted as a final category by being pulled by the determination result of the surrounding vehicle.
  • the present invention is not limited thereto.
  • the determination result of the other vehicle which is ahead of the own vehicle Hv and is within the distance D2 from the own vehicle Hv is used (in other words, (As a population), you may decide the final surrounding environment.
  • the deletion rate R for determining the category changes according to the pseudo distance distribution, which is a statistical index determined based on observation data within a predetermined time in the past. Even if the actual surrounding environment changes, it does not appear as a change in the deletion rate R until the pseudo distance distribution degree exceeds the predetermined distribution degree threshold. In other words, it takes time according to the number of data for calculating the pseudo distance distribution degree until the change of the actual surrounding environment is reflected in the determination result of the environment determination unit 16.
  • the determination result of the other vehicle traveling ahead of the own vehicle Hv is the environment surrounding the own vehicle Hv than the determination result of the other vehicle traveling behind the own vehicle Hv.
  • the category is determined by configuring the final peripheral environment category using the determination result of the other vehicle existing ahead of the host vehicle Hv and within the distance D2 from the host vehicle Hv. Accuracy can be improved.
  • the distance D2 may be approximately the same as the distance D1.
  • the determination result of the other vehicle whose traveling direction is the same as the host vehicle Hv and the determination result of the host vehicle Hv The category may be determined.
  • the surrounding environment is different between on an expressway and on a general road (hereinafter, a high-speed roadway) disposed below the expressway.
  • a general road hereinafter, a high-speed roadway
  • the view on the highway is better on the highway, and the view on the highway is worse on the high speed road.
  • final determination using the determination result of the host vehicle Hv of the other vehicle traveling on the same type of road as the host vehicle Hv It is preferable to determine the category of the surrounding environment. Whether the vehicle is traveling on a road of the same type as the host vehicle Hv may be specified by sharing the type of the traveling path by inter-vehicle communication, or determined based on height information included in the position information. May be In addition, since a significant difference occurs in the traveling speed between the expressway and the general road, it can be determined depending on whether the traveling speed is the same level (for example, within ⁇ 10 km).
  • the error indicator value may be a noise content rate that represents the ratio of noise to the positioning signal.
  • the noise content rate may be a value obtained by dividing the variance of the noise by the variance of the signal, or a value obtained by taking the logarithm thereof.
  • the noise content rate is a parameter that means that the smaller the value, the smaller the possibility that the pseudo distance calculated based on the positioning signal contains an error.
  • the noise content rate may be calculated as the reciprocal of the S / N ratio of the positioning signal.
  • the error index value may be an elevation residual determined by subtracting the elevation from 90 °. This is because the pseudorange residual ⁇ D tends to be smaller as the elevation angle is closer to 90 °, that is, as the elevation residual is smaller.
  • the inverse of the elevation angle may be adopted as the elevation residual.
  • the elevation residual may be a parameter having a smaller value as the elevation is larger.
  • the pseudorange residual ⁇ D directly represents the magnitude of the possibility that the pseudorange contains an error more than the noise content rate or the elevation residual. Therefore, it is preferable to use the pseudorange residual ⁇ D as the error index value.
  • the secondary filtering process may be performed using any two or all of the three indices such as the pseudorange residual ⁇ D, the noise content rate, and the elevation residual, in combination.
  • each unit is expressed as S101, for example.
  • each part can be divided into a plurality of sub-parts, while a plurality of parts can be combined into one part.
  • each part configured in this way can be referred to as a circuit, a device, a module, or a means.

Abstract

測位装置(1)は、信号受信部で測位信号を受信できている測位衛星である捕捉衛星毎に、その捕捉衛星に対して指標値散布度を算出する散布度算出部(13)と、捕捉衛星毎の指標値散布度に基づいて、複数の捕捉衛星のうち、測位演算処理に使用しない捕捉衛星である不使用衛星を決定する1次選択部(15)と、不使用衛星に設定された捕捉衛星の数の比である削除率に基づいて、周辺環境はマルチパスが発生しやすい環境であるか否かを判定する環境判定部(16)と、周辺環境の判定結果に基づいて、測位演算処理に用いる捕捉衛星の条件を変更する条件変更部(172)と、設定されている条件を充足している測位衛星からの測位信号を用いて測位演算処理を実施する測位演算部(18)と、を備える。

Description

測位装置 関連出願の相互参照
 本出願は、2017年7月7日に出願された日本特許出願番号2017-133972号に基づくもので、ここにその記載内容を援用する。
 本開示は、測位衛星から送信される測位信号を用いて測位を行う測位装置に関する。
 従来、全球測位衛星システム(GNSS:Global Navigation Satellite System)の技術分野では、測位衛星から送信された測位信号が建物などで反射されて測位装置で受信されること(いわゆるマルチパス)によって、演算結果としての位置情報に誤差が含まれうることが知られている。
 特許文献1では、そのようなマルチパスに由来する測位誤差を抑制するための構成として、疑似距離の誤差の大きさを示すパラメータである残差(以降、擬似距離残差とする)を捕捉衛星毎に算出し、擬似距離残差が予め定められている閾値を超過している捕捉衛星からの測位信号を測位演算処理に用いない構成が提案されている。なお、ここでの擬似距離残差とは、推定される測位装置の現在位置の座標と測位衛星の位置座標からとに基づいて定まる距離(座標由来値)と、測位信号の受信結果から算出される観測値としての擬似距離の差である。
特開2013-108961号公報
 特許文献1のように捕捉衛星毎に擬似距離残差を逐次算出し、前時刻において擬似距離残差が所定のレベル以下となっている捕捉衛星の観測データのみを用いて測位演算処理を実施する構成(以降、想定構成)では、擬似距離残差に対する閾値が、測位演算処理に用いる捕捉衛星を取捨選択する際の基準として作用する。
 一般的に、都市部等のマルチパス環境下では、オープンスカイ環境下に比べて個々の捕捉衛星からの測位信号の品質が劣化するため、測位演算処理に用いる測位衛星の数は多いほど、測位精度の向上が期待できる。そのような事情を鑑みると、都市部等のマルチパス環境下では、想定構成で設定される、擬似距離残差に対する閾値は大きい方が好ましい。仮に上記想定構成において、測位演算処理に用いる捕捉衛星を取捨選択するための擬似距離残差に対する閾値が小さ過ぎると、測位演算処理に使用可能な捕捉衛星が少なくなってしまい、測位精度を向上させる効果が得られにくくなってしまうためである。
 一方、オープンスカイ環境下では、マルチパス環境下に比べて個々の捕捉衛星からの測位信号の品質は良好であるため、測位演算処理に用いる捕捉衛星が相対的に少数であっても十分に高精度測位が可能である。また、擬似距離残差に対する閾値を小さい値に設定すれば、擬似距離等の誤差が大きい観測データが測位演算処理に用いる観測データに混ざり込む恐れを抑制することができ、測位精度のさらなる向上が期待できる。つまり、オープンスカイ環境下では、擬似距離残差に対する閾値は小さいほうが好ましい。
 このように、測位演算処理に用いる捕捉衛星(換言すれば観測データ)を取捨選択するための条件であって、より精度がよい測位結果を得るために設定されるべき条件は、周辺環境によって異なる。
 本開示は、この事情に基づいて成されたものであり、その目的とするところは、オープンスカイ環境下での測位精度の向上と、マルチパス環境下での測位精度の向上とを両立させることができる測位装置を提供することにある。
 本開示の一態様による測位装置は、複数の測位衛星のそれぞれから送信される測位信号を受信する信号受信部と、信号受信部が受信した測位信号に基づいて、受信した測位信号の送信元としての測位衛星についてのドップラーシフト量及び擬似距離の少なくとも何れか一方を、マルチパス環境下であるか否かを判定するための指標値として逐次算出する環境指標値算出部と、信号受信部で測位信号を受信できている測位衛星である捕捉衛星毎に、その捕捉衛星に対して環境指標値算出部が算出した指標値のばらつき度合いを示す指標値散布度を算出する散布度算出部と、散布度算出部が算出した捕捉衛星毎の指標値散布度に基づいて、複数の捕捉衛星のうち、測位演算処理に使用しない捕捉衛星である不使用衛星を決定する1次選択部と、捕捉衛星の総数に対する1次選択部によって不使用衛星に設定された捕捉衛星の数の比である削除率に基づいて、周辺環境はマルチパスが発生しやすい環境であるか否かを判定する環境判定部と、環境判定部の判定結果に基づいて、測位演算処理に用いる捕捉衛星の条件を変更する条件変更部と、条件変更部によって設定されている条件を充足している測位衛星からの測位信号を用いて測位演算処理を実施する測位演算部と、を備える。
 上記の構成では、条件変更部が、環境判定部の判定結果に応じて、測位演算処理に用いる捕捉衛星の条件(以降、採用条件)を動的に変更する。故に、例えば、環境判定部によって周辺環境はマルチパスが生じにくい環境であると判定されている場合には、条件変更部が採用条件を相対的に厳しい条件に変更することで、測位演算処理に使用する捕捉衛星を信号品質が良い捕捉衛星に絞り、測位精度の向上を図ることができる。また、環境判定部によって周辺環境はマルチパスを生じやすい環境であると判定されている場合には採用条件を緩和することで測位演算処理に使用する捕捉衛星の数を増やし、測位精度の向上を図ることもできる。
 つまり上記の構成によれば、条件変更部が採用条件を周辺環境に応じて動的に変更するため、オープンスカイ環境下での測位精度の向上と、マルチパス環境下での測位精度の向上とを両立させることができる。
 本開示についての上記目的およびその他の目的、特徴や利点は、添付の図面を参照しながら下記の詳細な記述により、より明確になる。その図面は、
図1は、測位装置の利用態様の一例を示した図であり、 図2は、測位装置の概略的な構成の一例を示すブロック図であり、 図3は、環境判定部の作動を説明するための図であり、 図4は、2次フィルタ部が実施する2次フィルタリング処理を説明するためのフローチャートであり、 図5は、2次フィルタ部が実施する例外処理を説明するためのフローチャートであり、 図6は、変形例1の2次フィルタ部の概略的な構成の一例を示すブロック図であり、 図7は、変形例1の閾値変更部の作動例を示す図であり、 図8は、変形例2における測位装置の概略的な構成の一例を示すブロック図であり、 図9は、変形例2における環境判定部の作動を説明するための図であり、 図10は、変形例2における環境判定部の作動を説明するための図である。
 本開示の実施形態としての測位装置1について、図を用いて説明する。図1は、本実施形態の測位装置(POSI APPA)1を含む衛星測位システム100の概略的な構成を示す図である。測位装置1は、概略的に、複数の測位衛星2、及び準天頂衛星3のそれぞれから送信される測位信号を受信することによって現在位置を逐次する装置である。測位装置1は、車両Hvに搭載されており、測位装置1の測位結果(つまり現在位置情報)は、例えば自動運転や経路案内等といった任意のアプリケーションソフトウェアで利用される。
 なお、車両Hvは、四輪自動車であってもよいし、二輪自動車や、三輪自動車等であってもよい。二輪自動車には原動機付き自転車も含まれる。車両Hvは例えば建設車両や農業用作業車など、乗員の移動以外を主目的とする車両であってもよい。また、測位装置1は車両以外に搭載されていても良い。例えば測位装置1は、スマートフォンやタブレット端末、ウェアラブル端末などといった、ユーザによって携帯される通信端末(以降、携帯端末)に搭載されていても良い。
 この測位装置1は、地上の任意の地域で使用される。ここでは一例として測位装置1は日本国内で使用される場合を想定して説明する。測位装置1の構成及び機能等についての詳細は別途後述する。
 測位衛星2は、所定の全球測位衛星システム(以降、GNSS:Global Navigation Satellite System)を構成する人工衛星である。GNSSとしては、GPSや、Galileo、GLONASS、BeiDou等がある。本実施形態では一例として、測位衛星2はGNSSの1つであるGPS(Global Positioning System)を構成する人工衛星であるものとして以下の構成を説明するが、これに限らない。測位衛星2は、GLONASSを構成する人工衛星であってもよいし、BeiDouを構成する人工衛星であってもよい。
 また、衛星測位システム100は、複数種類のGNSSを組み合わせたシステムであっても良い。換言すれば、測位装置1は、GPSや、Galileo、GLONASS、BeiDouなどの、複数のGNSSのそれぞれに対応しており、複数のGNSSの測位信号を組み合わせて用いて測位演算処理を実行可能に構成されていてもよい。たとえば測位装置1は、GPSを構成する測位衛星2から送信される測位信号と、GLONASSを構成する測位衛星2から送信される測位信号の両方を用いて測位演算処理を実施するものであっても良い。
 なお、図1では測位衛星を4機しか示していないが、実際にはもっと多数存在しうる。例えば、GPSを構成する測位衛星2は30機以上存在する。また、GLONASSを構成する測位衛星2は、20機以上存在する。
 各測位衛星2は、送信時刻等を示すデータを、測位衛星毎に固有のC/Aコードを用いて位相変調した信号(以降、測位信号)を逐次(例えば50ミリ秒毎に)送信する。測位信号は、送信時刻の他に、例えば、衛星時計の誤差を示すデータや、衛星自身の現在位置を示すエフェメリスデータ、全測位衛星の概略的な軌道情報を示すアルマナックデータなどを示す。種々のデータは順次送信される。C/Aコードは測位衛星2毎に固有であるため、C/Aコードは、送信元を示す情報として機能する。便宜上、測位装置1が測位信号を受信できている測位衛星2のことを、捕捉衛星とも記載する。
 準天頂衛星3は、測位装置1が使用される地域(つまり日本)を含む、特定の地域の上空に長時間とどまる軌道(いわゆる準天頂軌道)を周回している人工衛星であって、準天頂衛星システム(以降、QZSS:Quasi-Zenith Satellite System)を構成する。QZSSは、地域航法衛星システム(いわゆるRNSS:Regional Navigation Satellite System)の一例に相当する。なお、準天頂衛星3は任意の要素である。
 準天頂衛星3は、GPSと一体運用可能に構成されており、GPS衛星としての測位衛星2が送信する測位信号と同様の測位信号を逐次送信する。このように準天頂衛星3は、測位衛星2と同様の機能を提供するため、以降に記載する測位衛星2には準天頂衛星3も含まれるものとする。また、本開示の測位衛星は、GPS衛星としての測位衛星2以外、準天頂衛星3も含まれる。
 準天頂衛星3は、測位信号の他に、衛星回線を用いて所定の測位衛星2についての補正情報を示す補正信号を送信する。ここでの補正情報とは、例えば、対象とする測位衛星2の精密衛星座標、時計誤差、位相バイアス、電離層補正量、対流圏補正量などを示す情報である。電離層補正量は、例えば電離層の厚みに応じて生じる擬似距離や搬送波位相の誤差といった、電離層での遅延の影響を補正するパラメータである。対流圏補正量は、対流圏での遅延の影響(例えば擬似距離や搬送波位相の誤差)を補正するパラメータである。準天頂衛星3が送信する補正信号は、地上に設置されたセンタ(以降、地上局)で生成されて送信されてきたデータに基づいて生成されれば良い。補正信号は、SBAS補強メッセージと同様の信号とすることができる。
 <測位装置1の構成及び機能について>
 次に、図2を用いて測位装置1の構成及び機能について説明する。測位装置1は、図2に示すように、搬送波受信部11、受信処理部12、散布度算出部13、出力部14、1次フィルタ部15、環境判定部16、2次フィルタ部17、及び測位演算部18を備える。なお、受信処理部12、散布度算出部13、出力部14、1次フィルタ部15、環境判定部16、2次フィルタ部17、及び測位演算部18は、図示しないCPUが所定のプログラム(以降、測位プログラム)を実行することによって実現することができる。その場合、測位装置1は、種々の演算処理を実行するCPU、不揮発性のメモリであるフラッシュメモリ、揮発性のメモリであるRAM、I/O、及びこれらの構成を接続するバスラインなどを備える。
 通常のコンピュータを測位装置1として機能させるためのプログラム(つまり測位プログラム)は、非遷移的実体的記録媒体(non- transitory tangible storage medium)に格納されていればよい。例えば測位プログラムはROMやフラッシュメモリ等に保存されていれば良い。CPUが測位プログラムを実行することは、測位プログラムに対応する方法が実行されることに相当する。
 本実施形態では一例として測位装置1は、CPUが測位プログラムを実行することによって発現される機能ブロックとして受信処理部12、散布度算出部13、出力部14、1次フィルタ部15、環境判定部16、2次フィルタ部17、及び測位演算部18を備えているものとするが、これに限らない。上述した機能ブロックの一部又は全部は、ハードウェアとして実現されてもよい。ハードウェアとして実現する態様には、一つ或いは複数のIC等を用いて実現する態様も含まれる。さらに、上記の機能ブロックの一部又は全部は、CPUによるソフトウェアの実行とハードウェア部材の協働によって実現されてもよい。
 搬送波受信部11は、測位衛星2から送信される測位信号としての電波を受信し、電気信号に変換して出力する通信モジュールである。搬送波受信部11は、測位信号を受信するためのアンテナ(以降、GNSSアンテナ)を備える。また、搬送波受信部11は、少なくとも1つのGNSSアンテナの他に、受信信号をベースバンド帯の信号に変換するための周波数変換回路や、GNSSアンテナから出力されるアナログ信号をデジタル信号に変換するアナログデジタル変換部などを備える。搬送波受信部11は、増幅回路等を備えていても良い。搬送波受信部11が信号受信部に相当する。
 なお、測位装置1が複数のGNSSの測位信号を用いて測位演算処理を実施可能に構成されている場合には、搬送波受信部11は、各GNSSの搬送周波数に対応する複数のGNSSアンテナを備えていればよい。搬送波受信部11は測位装置1が対応すべき各測位システムの測位信号をそれぞれ受信可能に構成されていればよい。
 受信処理部12は、搬送波受信部11から出力される測位信号に基づいて、送信元の測位衛星についての観測データを生成する構成である。観測データは、受信した測位信号から定まる捕捉衛星についてのデータである。ここでは一例として、受信処理部12は、観測データとして、衛星番号、観測時刻、ドップラーシフト量、擬似距離、衛星座標、搬送波位相等を示すデータを生成して出力するものとする。衛星番号は、複数の測位衛星2のそれぞれを区別するための情報であって、例えば、PRN(Pseudo Random Noise)IDとすればよい。観測時刻は当該観測データを生成するための測位信号を受信した時刻とすればよい。
 ドップラーシフト量は、ドップラー効果によって生じる搬送波周波数と受信周波数の差を表すパラメータである。擬似距離は、測位信号が測位衛星2から送信されてから受信処理部12で受信されるまでの時間(以降、信号伝搬時間)と、電波の伝搬速度とから定まる距離である。なお、信号伝搬時間は、送信時刻と受信時刻との差から算出されても良いし、C/Aコードの位相のずれ量に基づいて算出されてもよい。衛星座標は、測位衛星の衛星軌道上における現在位置を示す情報である。
 なお、観測データとしては上述した全ての情報を含んでいる必要はなく、観測データが含むべき具体的な項目は適宜設計されれば良い。また、観測データには、アルマナックデータや、衛星座標の算出に用いたエフェメリスデータ等が含まれていても良い。観測データは、測位演算処理に必要な情報を含んでいればよい。受信処理部12は、上述した観測データを生成するためのサブ機能として、擬似距離を算出する擬似距離算出部や、ドップラーシフト量算出部などを備える。なお、擬似距離算出部や、ドップラーシフト量算出部の図示は省略している。このように擬似距離算出部やドップラーシフト量算出部としての機能を備える受信処理部12が環境指標値算出部に相当する。
 受信処理部12は、搬送波受信部11が測位信号を受信する度に上述した観測データを生成する。これにより、測位信号を受信できている測位衛星(つまり捕捉衛星)毎の観測データが逐次生成される。そして、受信処理部12が生成した捕捉衛星ごとの観測データは出力部14に出力される。また、受信処理部12が生成した擬似距離は、何れの捕捉衛星についての擬似距離であるかを示す情報(たとえば衛星番号)と対応づけられて散布度算出部13に出力される。
 散布度算出部13は、受信処理部12から逐次提供される捕捉衛星毎の擬似距離を取得し、捕捉衛星毎に区別して図示しないメモリに保存していく。メモリはRAM等の書き換え可能な記憶媒体を用いて実現されれば良い。同一捕捉衛星についての取得時点が異なる複数の疑似距離は、例えば、最新のデータが先頭となるように時系列順にソートされてメモリに保存されれば良い。また、保存されてから一定時間経過したデータは順次破棄されていけば良い。
 散布度算出部13は、捕捉衛星毎に、直近一定時間以内において算出(換言すれば観測)された擬似距離のばらつき度合いを示す擬似距離散布度を算出する。或る捕捉衛星についての擬似距離散布度は、現時点から過去一定時間(例えば5秒)以内において算出された擬似距離を母集団とする標準偏差とすれば良い。なお、他の態様として散布度は、分散や、最大値と最小値の差など、統計学においてデータのばらつき度合いを示す他の指標を採用することもできる。
 或る捕捉衛星についての擬似距離散布度は、当該捕捉衛星についての新たな擬似距離を取得した場合に算出されれば良い。また、散布度算出部13は、所定の周期で各捕捉衛星についての擬似距離散布度を逐次(例えば100ミリ秒毎に)算出するように構成されていても良い。
 なお、他の態様として、散布度算出部13は、擬似距離散布度の代わりに、直近一定時間以内において算出観測されたドップラーシフト量のばらつき度合いを示すシフト量散布度を算出するように構成されていても良い。その場合、受信処理部12は、擬似距離の代わりにドップラーシフト量を散布度算出部13に出力するように構成されていればよい。
 また、他の態様として、散布度算出部13は、擬似距離散布度とシフト量散布度の両方を捕捉衛星毎に算出するように構成されていても良い。その場合、受信処理部12は、擬似距離とドップラーシフト量の両方を散布度算出部13に逐次出力するように構成されていればよい。
 なお、車両Hvの周辺環境がマルチパス環境下である場合には、逐次算出される擬似距離のばらつきは大きくなるため、擬似距離は周辺環境がマルチパス環境下であるか否かの指標値として機能する。また、車両Hvの周辺環境がマルチパス環境下である場合には、測位信号の到来方向が逐次変化し得る。到来方向が変わればドップラーシフト量も変化する。つまり、車両Hvの周辺環境がマルチパス環境下である場合には、逐次算出されるドップラーシフト量のばらつきは大きくなるため、逐次算出されるドップラーシフト量もまた、周辺環境がマルチパス環境下であるか否かの指標値として機能する。すなわち、擬似距離及びドップラーシフト量が指標値に相当する。また、擬似距離散布度及びシフト量散布度が指標値散布度に相当する。
 散布度算出部13が算出した擬似距離散布度は、何れの捕捉衛星についての擬似距離散布度であるかを示す情報(たとえば衛星番号)と対応づけられて出力部14に出力される。出力部14は、受信処理部12から入力された捕捉衛星毎の観測データ、及び、散布度算出部13から入力された捕捉衛星毎の擬似距離散布度を1次フィルタ部15に出力する。
 1次フィルタ部15は、捕捉衛星毎の観測データのうち、擬似距離散布度が所定の閾値(以降、散布度閾値)以上となっている捕捉衛星の観測データを削除する。そして、1次フィルタ部15は、擬似距離散布度が散布度閾値以上となっている捕捉衛星の観測データを削除した後に残っている観測データを2次フィルタ部17に出力する。つまり、1次フィルタ部15は、擬似距離散布度が相対的に小さい捕捉衛星についての観測データのみが後段の構成である2次フィルタ部17で利用されるように、捕捉衛星毎の観測データをフィルタリングする構成である。
 便宜上、捕捉衛星毎の観測データのうち、擬似距離散布度が散布度閾値以上となっている捕捉衛星の観測データを削除する処理のことを一次フィルタリング処理と称する。1次フィルタ部15が1次選択部に相当する。また、一次フィルタリング処理の結果として残った観測データに対応する捕捉衛星のことを一次残存衛星とも称する。
 なお、或る捕捉衛星についての観測データを後段処理に出力しない(つまり削除)するということは、当該捕捉衛星を測位演算処理に使用しないことに相当する。つまり、1次フィルタ部15は、捕捉衛星毎の擬似距離散布度に基づいて、複数の捕捉衛星のなかから測位演算処理に使用しない捕捉衛星である不使用衛星を決定する構成に相当する。換言すれば、捕捉衛星毎の擬似距離散布度に基づいて、複数の捕捉衛星のなかから測位演算処理に使用する捕捉衛星を取捨選択する構成に相当する。
 また、1次フィルタ部15は、出力部14から入力された観測データの数(換言すれば捕捉衛星の総数)と、削除した観測データの数の比を表す削除率Rを算出する。出力部14から入力された観測データの数をN、削除した観測データの数をkとした場合、削除率Rは、k/Nによって表される。なお、出力部14から入力された観測データの数をN、削除した観測データの数をkとした場合に、2次フィルタ部17に出力される観測データの数Mは、N-kである。1次フィルタ部15は、削除率Rを示すデータを環境判定部16に出力する。
 環境判定部16は、1次フィルタ部15から提供される削除率Rに基づいて、測位装置1の周辺環境が、上空視界において測位衛星からの信号を遮断及び散乱させる地物が多い環境(換言すればマルチパスが発生しやすい環境)であるか否かを判定する構成である。マルチパスが発生しやすいかどうかは、測位装置1上空の見晴らしの良さによって定まる。
 本実施形態では一例として、周辺環境を上空の見晴らしの良さに応じて段階的に4つのカテゴリに分類しており、環境判定部16は削除率Rに応じて現在の周辺環境がカテゴリのうちの何れのカテゴリに該当するのかを判定するものとする。以降では4つのカテゴリを上空の見晴らしが良い順に、第1カテゴリC1、第2カテゴリC2、第3カテゴリC3、第4カテゴリC4と称する。第1カテゴリC1が最も上空の見晴らしが良く、第4カテゴリC4が最も上空の見晴らしが悪い環境に相当する。
 上空の見晴らしが悪いほど、ビル等の遮蔽物が多く、マルチパスが発生しやすい環境であることを意味する。また、マルチパスが発生しやすいほど、擬似距離散布度は相対的に大きい値となり、その結果として、削除率Rも相対的に大きい値へとなる。故に、環境判定部16は、削除率Rが大きいほど、周辺環境はカテゴリ番号が大きいカテゴリに該当すると判定する。
 具体的には、環境判定部16は、削除率Rが所定の第1削除率Rth1未満である場合に、周辺環境は第1カテゴリC1に該当すると判定する。第1削除率Rth1の具体的な値は試験を実施して適宜設定されれば良い。ここでは一例として10%とするが、その他、5%や15%などであってもよい。また、削除率Rが、第1削除率Rth1以上、且つ、所定の第2削除率Rth2未満である場合には周辺環境はカテゴリC2に該当すると判定する。第2削除率Rth2の具体的な値は、第1削除率Rth1よりも大きい範囲において、試験を実施して適宜設定されれば良く、ここでは一例として30%とする。
 さらに、環境判定部16は、削除率Rが、第2削除率Rth2以上、且つ、所定の第3削除率Rth3未満である場合には周辺環境は第3カテゴリC3に該当すると判定する。第3削除率Rth3の具体的な値は、第2削除率Rth2よりも大きい試験を実施して適宜設定されれば良く、ここでは一例として60%とする。そして、削除率Rが、第3削除率Rth3以上である場合には周辺環境は第4カテゴリC4に該当すると判定する。
 環境判定部16は、1次フィルタ部15から削除率Rが入力される度に、周辺環境の判定を実施する。当然、車両Hvが走行するにつれて周辺環境は遷移し、削除率Rを変化する。すなわち、図3に示すように、削除率Rの推移に伴って環境判定部16の判定結果もまた経時的に変化していく。環境判定部16の判定結果は2次フィルタ部17に逐次出力される。
 なお、各カテゴリが指す環境は、具体的には次の通りである。第1カテゴリC1は、いわゆるオープンスカイ環境に対応するカテゴリである。オープンスカイ環境は、定量的には、例えば魚眼カメラにて上空を撮影した場合の天空比(換言すれば天空率)が80%以上となる環境とすればよい。もちろん、天空比が75%や85%以上となる環境をオープンスカイ環境と定義しても良い。測位装置1がオープンスカイ環境下に存在する場合、マルチパスが発生しにくいため、擬似距離散布度は相対的に小さい値となることが期待され、その結果として、削除率Rも相対的に小さい値へとなることが期待できる。前述の第1削除率Rth1はオープンスカイ環境下での削除率Rの試験結果に基づいて決定されれば良い。なお、オープンスカイ環境下では、各測位衛星からの測位信号が、ビル等で反射されること無く、直接的に受信できる可能性が高い。そのため、オープンスカイ環境は、相対的に高精度な測位結果が得られる環境である。第1カテゴリC1がオープンスカイカテゴリに相当する。
 第2カテゴリC2は、市街地の中でも相対的にビル等の構造物が少ない(換言すればマルチパスが発生しにくい)エリアを指すカテゴリである。例えば第2カテゴリC2は、天空比が相対的に低い値(例えば65%以上、80%未満)となる環境に相当するものとすればよい。
 第3カテゴリC3は、市街地の中でも相対的にビル等の構造物が多い(換言すればマルチパスが発生しやすい)エリアを指すカテゴリである。第2カテゴリC2と第3カテゴリC3とを切り分ける基準は試験結果に基づいて適宜決定されれば良い。例えば第3カテゴリC3は、天空比が相対的に低い値(例えば50%以上、65%未満)となる環境に相当するものとすればよい。第4カテゴリC4は、高層ビル等が存在する都市部を指すカテゴリである。例えば第4カテゴリC4は、天空比が50%未満となるような環境を想定したカテゴリとすればよい。第4カテゴリC4が都市部カテゴリに相当する。
 なお、本実施形態では、周辺環境を上空視界に応じてカテゴリC1~C4の4段階に分類しているがこれに限らない。周辺環境は2段階や3段階、5段階に分類しても良い。例えばオープンスカイ環境であるか否かの2段階に分類分けされていても良い。また、本実施形態の第2カテゴリC2と第3カテゴリC3とを統合して1つのカテゴリとすることで、全体として周辺環境を3つのカテゴリに分けても良い。森の中などの環境やトンネル内部などをカテゴリとして追加してもよい。
 2次フィルタ部17は、1次フィルタ部15から提供される一次残存衛星毎の観測データから、相対的にマルチパス等の影響を強く受けている観測データを削除して、残っている観測データを測位演算部18に出力する構成である。つまり、2次フィルタ部17は、一次残存衛星毎の観測データの中から、最終的に測位演算処理に用いられる観測データを抽出する構成に相当する。2次フィルタ部17は、サブ機能として、残差算出部171及び閾値変更部172を備える。
 この残差算出部171、閾値変更部172を含む2次フィルタ部17は、概略的に、次のように作動する。2次フィルタ部17はまず、1次フィルタ部15による一次フィルタリング処理の結果として抽出される一次残存衛星のそれぞれについての擬似距離残差ΔDを算出する。そして、擬似距離残差ΔDが所定の残差閾値Dth未満となっている捕捉衛星の観測データのみを測位演算部18に出力する。前述の残差算出部171は、各一次残存衛星についての擬似距離残差ΔDを算出する構成である。
 或る測位衛星2についての擬似距離残差ΔDとは、測位装置1の現在位置から測位衛星2までの直線距離と、測位信号が送信されてから測位装置1で受信されるまでの時間差と電波の伝搬速度に基づいて算出される距離(すなわち擬似距離)との差の絶対値である。現在の測位装置1から測位衛星2までの直線距離は、現在の測位装置1の位置座標と測位衛星2の位置座標とから算出することができる。
 以上のようにして算出される擬似距離残差ΔDは、測位信号がマルチパス等の影響を強く受け、擬似距離が誤差を含んでいる場合には、大きい値をとる傾向がある。つまり、或る捕捉衛星についての擬似距離残差ΔDは当該捕捉衛星から送信された測位信号に基づいて算出された擬似距離が誤差を含んでいる可能性の大きさを表すパラメータとして機能する。故に、擬似距離残差ΔDは誤差指標値に相当する。また、残差算出部171が誤差指標値算出部に相当する。
 なお、擬似距離残差ΔDを算出する際に用いる測位装置1の現在位置としての位置座標は、厳密には、前時刻において測位演算部18が算出した位置座標が使用されれば良い。測位演算部18は基本的に約50~100ミリ秒毎に測位演算処理を実施するため、前時刻において測位演算部18が算出した位置座標とは、現在よりも50~100ミリ秒前の測位装置1の位置を示すものである。擬似距離残差ΔDを算出する上では、前時刻において測位演算部18が算出した位置座標と実際の現在位置との誤差は無視できるものとする。もちろん、前時刻において測位演算部18が算出した位置座標を元に周知のデッドレコニングを実施することによって、より厳密な現在位置を推定し、擬似距離残差ΔDを算出しても良い。
 便宜上、一次残存衛星毎の観測データのうち、擬似距離残差ΔDが所定の残差閾値以上となっている捕捉衛星の観測データを削除する処理のことを2次フィルタリング処理と称する。また、2次フィルタリング処理の結果として残った観測データに対応する捕捉衛星のことを2次残存衛星とも称する。
 2次フィルタリング処理で使用される残差閾値Dthは、測位演算処理に用いる捕捉衛星(換言すれば観測データ)を取捨選択するための擬似距離残差に対する閾値として機能する。故に残差閾値Dthがフィルタリング閾値に相当する。2次フィルタ部17が備える閾値変更部172は、環境判定部16が判定した周辺環境のカテゴリに基づいて残差閾値Dthの値を変更する構成である。
 閾値変更部172は、環境判定部16によって周辺環境はマルチパスが生じにくい環境(つまり番号が小さいカテゴリ)であると判定されている場合には、残差閾値Dthを相対的に小さい値に設定する。また、環境判定部16によって周辺環境はマルチパスを生じやすい環境(つまり番号が大きいカテゴリ)であると判定されている場合には、残差閾値Dthを相対的に大きい値に設定する。閾値変更部172が残差閾値Dthとして設定可能な複数段階の値(換言すれば候補値)は、予め設定されている。複数段階の候補値は、複数段階の各カテゴリに対応づけられている。この閾値変更部172の作動の詳細については別途図4等を用いて後述する。
 なお、残差閾値Dthを相対的に大きい値に設定すれば、一次残存衛星は測位演算処理に使用する測位衛星として選ばれやすくなり、残差閾値Dthを相対的に小さい値に設定すれば、一次残存衛星は測位演算処理に使用する測位衛星として選ばれにくくなる。つまり、残差閾値Dthは、測位演算処理に使用する捕捉衛星の条件を規定するパラメータとして機能する。
 すなわち、閾値変更部172は、環境判定部16の判定結果に基づいて測位演算処理に使用される捕捉衛星の条件を緩めたり厳しくしたりする構成である。故に、この閾値変更部172が、条件変更部に相当する。前述の2次残存衛星とは、閾値変更部172によって設定されている測位演算処理に使用する捕捉衛星の条件を充足している捕捉衛星に相当する。
 測位演算部18は、2次フィルタ部17から提供された観測データを用いて測位演算処理を実施する。すなわち、2次フィルタ部17による2次フィルタリング処理の結果として残った2次残存衛星についての観測データを用いて測位演算処理を実施する。測位演算処理の方式(以降、測位方式)自体は、例えばPPP-AR(より具体的にはMADOCA-PPP)や、PPP-RTK等、周知の種々の方式を採用することができる。各測位方式の具体的な手順は公知であるためここではその詳細については省略する。
 なお、PPPは、Precise Point Positioning(つまり単独搬送波位相測位)の略であり、ARは、Ambiguity Resolutionの略である。MADOCAは、Multi-gnss Advanced Demonstration tool for Orbit and Clock Analysisの略である。RTKは、Real-Time Kinematicの略である。
 測位演算部18は、測位演算処理の結果としての位置情報を所定のアプリケーションに出力する。位置情報を利用するアプリケーションは、測位装置1に組み込まれていても良いし、測位装置1の外部に設けられた電子機器に組み込まれていても良い。また測位演算部18は、測位演算処理の結果を2次フィルタ部17にも出力する。2次フィルタ部17に提供された位置情報は、次回の擬似距離残差ΔDの算出に供される。
 <2次フィルタ部17の作動について>
 次に図4に示すフローチャートを用いて、2次フィルタ部17が実施する2次フィルタリング処理について説明する。図4に示すフローチャートは、1次フィルタ部15から観測データが入力される度に開始されれば良い。
 まずステップS101では残差算出部171が、1次残存衛星毎の観測データ(主として擬似距離及び衛星座標)と、測位装置1の位置情報とに基づいて、1次残存衛星毎の擬似距離残差ΔDを算出してステップS101に移る。
 ステップS102では閾値変更部172が、環境判定部16による周辺環境の判定結果が第1カテゴリC1であるか否かを判定する。環境判定部16による周辺環境の判定結果が第1カテゴリC1である場合には、ステップS102を肯定判定してステップS103を実行する。一方、環境判定部16による周辺環境の判定結果が第1カテゴリC1ではない場合には、ステップS102を否定判定してステップS104を実行する。
 ステップS103では閾値変更部172が、残差閾値Dthを、残差閾値Dthの候補値として予め用意されている所定の第1閾値Dth1に設定してステップS108に移る。第1閾値Dth1の具体的な値は、オープンスカイ環境下における擬似距離残差の試験結果に基づいて決定されれば良い。例えば第1閾値Dth1は2~5m程度の値に設定されれば良い。第1閾値Dth1がオープンスカイカテゴリとしての第1カテゴリC1に対応付けられている候補値に相当する。
 ステップS104では、環境判定部16による周辺環境の判定結果が第2カテゴリC2であるか否かを判定する。環境判定部16による周辺環境の判定結果が第2カテゴリC2である場合には、ステップS104を肯定判定してステップS105を実行する。一方、環境判定部16による周辺環境の判定結果が第2カテゴリC2ではない場合には、ステップS104を否定判定してステップS106を実行する。
 ステップS105では閾値変更部172が、残差閾値Dthを所定の第2閾値Dth2に設定してステップS108に移る。第2閾値Dth2の具体的な値は、第1閾値Dth1よりも大きい範囲において適宜設定されれば良い。第2閾値Dth2は、残差閾値Dthが第1閾値Dth1に設定されている場合に比べて2次残存衛星の数が1~2以上増加し、測位演算処理に使用する測位衛星の数(以降、使用衛星数)の増加による測位精度の向上効果が得られるレベルの値に設定されていることが好ましい。上記効果が得られる具体的な値は、第2カテゴリC2における試験結果に基づいて決定されれば良い。第2閾値Dth2が第2カテゴリC2に対応付けられている候補値に相当する。
 ステップS106では、環境判定部16による周辺環境の判定結果が第3カテゴリC3であるか否かを判定する。環境判定部16による周辺環境の判定結果が第3カテゴリC3である場合には、ステップS106を肯定判定してステップS107を実行する。一方、環境判定部16による周辺環境の判定結果が第3カテゴリC3ではない場合には、ステップS106を否定判定してステップS108を実行する。なお、環境判定部16による周辺環境の判定結果が第3カテゴリC3ではない場合とは、環境判定部16による周辺環境の判定結果が第4カテゴリC4である場合に相当する。
 ステップS107では閾値変更部172が、残差閾値Dthを所定の第3閾値Dth3に設定してステップS108に移る。第3閾値Dth3の具体的な値は、第2閾値Dth2よりも大きい範囲において適宜設定されれば良い。第3閾値Dth3は、残差閾値Dthが第2閾値Dth2に設定されている場合に比べて2次残存衛星の数が1~2以上増加し、使用衛星数の増加による測位精度の向上効果が得られるレベルの値に設定されていることが好ましい。上記効果が得られる具体的な値は、第3カテゴリC3における試験結果に基づいて決定されれば良い。第3閾値Dth3が第3カテゴリC3に対応付けられている候補値に相当する。
 ステップS108では閾値変更部172が、残差閾値Dthを所定の第4閾値Dth4に設定してステップS108に移る。第4閾値Dth4の具体的な値は、第3閾値Dth3よりも大きい範囲において適宜設定されれば良い。第4閾値Dth4は、残差閾値Dthが第3閾値Dth3に設定されている場合に比べて2次残存衛星の数が1~2以上増加し、使用衛星数の増加による測位精度の向上効果が得られるレベルの値に設定されていることが好ましい。上記効果が得られる具体的な値は、第4カテゴリC4における試験結果に基づいて決定されれば良い。なお、残差閾値として設定可能な候補値の最大値に相当する第4閾値Dth4に限っては、1次残存衛星を全て測位演算処理に使用できるように、無限大に相当する十分大きい値に設定されていても良い。第4閾値Dth4が都市部テゴリとしての第4カテゴリC4に対応付けられている候補値に相当する。
 ステップS109では2次フィルタ部17が、1次フィルタ部15から入力された捕捉衛星毎の観測データのうち、擬似距離残差ΔDが上記処理で設定されている残差閾値Dth未満となっている捕捉衛星についての観測データを抽出してステップS110に移る。
 ステップS110では2次フィルタ部17が、ステップS109で抽出された観測データに基づいて測位演算部18が測位演算処理を実行可能であるか否かを判定する。例えば、ステップS109で抽出された観測データの数(換言すれば2次残存衛星の数)が4以上である場合に、測位演算部18が測位演算可能であると判定する。一方、ステップS109で抽出された観測データの数が4未満である場合には、測位演算部18が測位演算不可であると判定する。
 ステップS109で抽出された観測データの数が、測位演算部18が測位演算可能なほど十分に多い場合には、S110を肯定判定してステップS111を実行する。一方、ステップS109で抽出された観測データの数が、測位演算部18が測位演算不可能なほど少ない場合には、ステップS110を否定判定してステップS112を実行する。なお、測位装置1の高さは前時刻における高さと同じであるという推定を利用する場合には、測位演算部18が測位演算可能であるか否かの判定基準とする観測データの数は3であっても良い。
 ステップS111では2次フィルタ部17が、ステップS109で抽出した観測データを測位演算処理に出力し、本フローを終了する。なお、ステップS111以降の処理としては、測位演算部18は2次フィルタ部17から提供された捕捉衛星毎の観測データを用いて測位演算処理を実施し、現在位置を特定する。そして、測位結果としての位置情報は所定のアプリケーション及び2次フィルタ部17へと出力される。
 ステップS112では2次フィルタ部17が、所定の例外処理を実施する。例外処理はS101~S109での一連の処理によって抽出された観測データでは、測位演算部18が測位演算処理を実施できない場合の処理である。この例外処理については図5に示すフローチャートを用いて説明する。図5に示すフローチャートは、図4のステップS112に移った際に実行される。
 まずステップS201では閾値変更部172が、残差閾値Dthの値を、現在の値よりも1段階大きい値に設定してステップS202に移る。具体的には、現在の残差閾値Dthが第1閾値Dth1に設定されている場合には残差閾値Dthを第2閾値Dth2に設定する。また、現在の残差閾値Dthが第2閾値Dth2に設定されている場合には残差閾値Dthを第3閾値Dth3に設定する。同様に、現在の残差閾値Dthが第3閾値Dth3に設定されている場合には残差閾値Dthを第4閾値Dth4に設定する。
 なお、現在の残差閾値Dthが第4閾値Dth4に設定されている場合には、測位演算不可として本フローを終了する。その場合、最後に測位演算処理を実行できた時の測位結果を起点としたデッドレコニングを実施し、仮の現在位置を推定するものとする。
 ステップS202では、1次フィルタ部15から入力された捕捉衛星毎の観測データのうち、擬似距離残差ΔDがステップS201で新たに設定された残差閾値Dth未満となっている捕捉衛星についての観測データを抽出してステップS203に移る。
 ステップS203では、ステップS202で抽出された観測データに基づいて測位演算部18が測位演算処理を実行可能であるか否かを判定する。このステップS203の実体的な処理内容は、前述のステップS110と同様であるため説明は省略する。ステップS202で抽出された観測データの数が、測位演算部18が測位演算可能なほど多い場合には、ステップS203を肯定判定してステップS204を実行する。一方、ステップS202で抽出された観測データの数が、測位演算部18が測位演算不可能なほど少ない場合には、ステップS203を否定判定してステップS201に戻る。
 ステップS204では、測位演算部18に対して測位演算処理の結果を外部アプリケーションに出力することを禁止した上で、ステップS202で抽出した観測データを出力し、ステップS205に移る。なお、このステップS204を受けて測位演算部18は、2次フィルタ部17から提供された観測データを用いて測位演算処理を実施し、その結果を2次フィルタ部17に返送する。また、2次フィルタ部17からの指示に基づき、外部のアプリケーションには測位結果を出力しない。すなわち、環境判定部16の判定結果に応じた残差閾値Dthを充足する観測データのみを用いて測位演算処理を実行できなかった場合には、対外的にはその測位結果を破棄する。
 ステップS205では、測位演算部18での測位演算処理の結果を取得して本フローを終了する。このステップS205で取得した位置情報は、次回の擬似距離残差を算出する際に使用される。
 <実施形態の効果について>
 ここでは第1比較構成として、1次フィルタ部15に相当する構成を備えず、かつ、周辺環境によらずに一定の残差閾値を用いて測位演算処理に用いる捕捉衛星を取捨選択する測位装置を導入し、本実施形態の効果について説明する。
 一般的に、オープンスカイ環境では、都市部等のマルチパス環境下に比べて個々の捕捉衛星からの測位信号の品質が良い傾向があるため、使用する観測データ(換言すれば捕捉衛星)の数を絞っても十分な情報が得られる。つまり、オープンスカイ環境下では、測位演算処理に用いる捕捉衛星が相対的に少数であっても十分に高精度測位が可能である。また、オープンスカイ環境下であっても、仰角が低いなどの事情によって相対的に誤差が大きい擬似距離を提供する捕捉衛星も存在しうる。
 そのため、オープンスカイ環境下では、観測データの数を増やせば増やすほど、必ずしも測位精度が向上するわけではない。仮に第1比較構成において残差閾値が相対的に大きい値に設定されている場合には、擬似距離等の誤差が大きい観測データが測位演算処理に用いる観測データに混ざり込む恐れが高まるため、かえって精度が劣化してしまう懸念がある。
 一方、マルチパス環境下では、オープンスカイ環境下に比べて個々の捕捉衛星からの測位信号の品質が劣化するため、測位演算処理に用いる測位衛星の数は多いほど、測位精度の向上が期待できる。すなわち、残差閾値はある程度大きいことが好ましい。仮に第1比較構成において残差閾値が相対的に小さい値に設定されている場合には、測位演算処理に用いる捕捉衛星の数を減ってしまい、使用衛星数の増加による精度向上の効果が得られにくくなるからである。
 このようにより精度がよい測位結果が得るために設定されるべき、測位演算処理に用いる捕捉衛星(換言すれば観測データ)を取捨選択する条件(以降、採用条件)は、周辺環境に応じて異なる。そのため、第1比較構成では、オープンスカイ環境下での測位精度の向上と、マルチパス環境下での測位精度の向上とを両立させることが困難である。
 対して、本実施形態の構成では、閾値変更部172が、環境判定部16によって周辺環境はマルチパスが生じにくい環境(つまり番号が小さいカテゴリ)であると判定されている場合には、残差閾値Dthを相対的に小さい値に設定する。これにより、信号品質がよい測位衛星からの観測データのみが測位演算処理に使用されることとなる。オープンスカイ環境下では測位演算処理に使用する捕捉衛星の数が少なくても十分に高精度測位が可能である為、上記処理によって測位演算処理に使用する捕捉衛星を信号品質が良い捕捉衛星に絞ることで、より一層高精度化を図ることができる。
 また、閾値変更部172は、環境判定部16によって周辺環境はマルチパスを生じやすい環境(つまり番号が大きいカテゴリ)であると判定されている場合には、残差閾値Dthを相対的に大きい値に設定する。これにより、採用条件が緩和され、測位演算処理に使用する捕捉衛星の数を増やすことができる。マルチパス環境下では測位演算処理に使用する捕捉衛星の数が多いほうが、測位結果が高精度化する傾向がある為、上記処理によって測位演算処理に使用する捕捉衛星の数を増やすことで、より一層高精度化を図ることができる。
 つまり本実施形態の構成によれば、採用条件を周辺環境に応じて動的に変更するため、オープンスカイ環境下での測位精度の向上と、マルチパス環境下での測位精度の向上とを両立させることができる。
 また、本実施形態の測位装置1は、1次フィルタ部15を備えることによって、より一層測位精度を向上させることが期待できる。具体的には次の通りである。マルチパス環境下では、相対的に擬似距離の変動が大きい。当然、マルチパスの影響を強く受けている捕捉衛星の擬似距離は変動が大きいため、擬似距離残差もまた変動が激しい。つまり、或る捕捉衛星の擬似距離残差が、或る時点においては小さかったとしても、次の時点でも小さいとは限らない。故に、第1比較構成では、擬似距離等の誤差が大きい観測データが、測位演算処理に用いる観測データに混ざり込む恐れが残っている。
 対して、本実施形態の構成によれば、1次フィルタ部15が、捕捉衛星毎の観測データのうち、擬似距離散布度が相対的に小さい捕捉衛星についての観測データのみを後段の構成である2次フィルタ部17に出力する。故に、擬似距離の変動が大きい捕捉衛星についての観測データは1次フィルタリング処理で除去されるため、擬似距離等の誤差が大きい観測データが、測位演算処理に用いる観測データに混ざり込む恐れを抑制することができる。このように本実施形態の構成によれば、瞬時的に擬似距離残差が残差閾値Dthを下回るようなケースにも対応できており、測位精度のさらなる向上が期待できる。
 さらに、上記構成では測位信号の受信状況から定まる削除率Rに基づいて周辺環境のカテゴリを判定する。周辺環境のカテゴリを判定する他の構成(以降、第2比較構成)としては、地図データを用いる構成も考えられる。しかしながら、そのような第2比較構成では、測位装置が、地図データが格納された記憶装置、又は、外部サーバと通信して地図データを取得する機能を備えている必要があり、その分だけコストが増大してしまう。対して本実施形態の環境判定部16は、測位信号の受信状況から定まる削除率Rに基づいて周辺環境のカテゴリを判定するため、測位装置1が、地図データを記憶している装置や外部サーバと通信する機能を備えている必要はない。つまり、本実施形態によれば第2比較構成よりも製造コストを抑制することができる。また、広域通信用モジュールを備えないデバイスにも適用できる。
 また、捕捉衛星数が10のときに5つの観測データを削除した場合(以降、第1ケース)と、捕捉衛星数が30のときに5つの観測データを削除した場合(以降、第2ケース)とでは、観測データの削除数は同じであっても、マルチパスの発生度合いは異なる。当然、第1ケースのほうが、マルチパスの発生し易い環境(換言すれば周囲に構造物が多い環境)であることを示唆している。そのような事情を鑑みれば、仮に観測データの削除数で周辺環境を判別しようとすると、周辺環境を誤判定してしまう恐れが相対的に高い。対して、本実施形態の構成のように、環境判定部16は、観測データを削除した数ではなく、削除した比率(つまり削除率R)によって周辺環境のカテゴリを判別すれば、周辺環境を誤判定してしまう恐れを抑制することができる。換言すれば、周辺環境を判別するための指標として観測データの削除数ではなく削除率Rを用いることにより、より精度良く周辺環境を判別することができる。
 以上、本開示の実施形態の一例について説明したが、本開示は上述の実施形態に限定されるものではなく、以降で述べる種々の変形例も本開示の技術的範囲に含まれ、さらに、下記以外にも要旨を逸脱しない範囲内で種々変更して実施することができる。もちろん、種々の変形例は適宜組み合わせて実施することができる。
 なお、前述の実施形態で述べた部材と同一の機能を有する部材については、同一の符号を付し、その説明を省略する。また、構成の一部のみに言及している場合、他の部分については先に説明した実施形態の構成を適用することができる。
 (変形例1)
 閾値変更部172は、環境判定部16の判定結果に加えて、準天頂衛星3を捕捉できているか否かによって残差閾値Dthを決定しても良い。例えば、準天頂衛星3を捕捉できている場合には、準天頂衛星3を捕捉できていない場合よりも残差閾値Dthを小さい値に設定しても良い。ここでは上記の思想に基づく構成の一例について、変形例1として説明する。
 変形例1における2次フィルタ部17は、図6に示すように、残差算出部171、閾値変更部172に加えて、QZS捕捉判定部173を備える。QZS捕捉判定部173は、準天頂衛星3からの測位信号を受信できているか否か(換言すれば捕捉できているか否か)を判定する構成である。準天頂衛星3を捕捉できているか否かは、捕捉衛星の仰角や測位信号が示す送信元情報によって判定することができる。
 例えばQZS捕捉判定部173は、仰角が所定の準天頂角度以上となっている捕捉衛星が存在する場合には準天頂衛星3を捕捉できていると判定すればよい。ここでの準天頂角度とは、地球を一周回る軌道をとっている測位衛星2(例えばGPS衛星)と、準天頂衛星3とを切り分けるための角度である。準天頂角度は天頂方向と見なすことができる角度の下限値(例えば70°)に設定されれば良い。なお、仰角は、捕捉衛星と測位装置1とを結ぶ直線が、測位装置1の現在位置を基準とする地平面に対してなす角度である。仰角は、測位装置1と捕捉衛星のそれぞれの位置座標に基づき算出することができる。
 また、QZS捕捉判定部173は、測位信号が示す送信元情報(例えばPRN-ID)に基づいて判定しても良い。さらに、準天頂衛星3が送信する補正信号を受信できている場合に、準天頂衛星3を捕捉できていると判定しても良い。QZS捕捉判定部173が準天頂捕捉判定部に相当する。
 変形例1の閾値変更部172は、QZS捕捉判定部173によって準天頂衛星3を捕捉できていないと判定されている場合には、残差閾値Dthを環境判定部16によって判定されているカテゴリに対応づけられている値に設定する。一方、QZS捕捉判定部173によって準天頂衛星3を捕捉できていると判定されている場合には、残差閾値Dthを環境判定部16によって判定されているカテゴリに対応づけられている値よりも小さい値に設定する。例えば図7に示すように、QZS捕捉判定部173によって準天頂衛星3を捕捉できていると判定されている場合には、残差閾値Dthを、準天頂衛星3を捕捉できていないと判定されている場合に適用される値を0.8倍した値に設定すればよい。
 なお、ここでの倍率は、0.8に限らず、0.9や、0.7などであっても良い。また、所定の倍率の係数を乗じることによって残差閾値Dthを低減するのではなく、所定の値(例えば0.5m)を減算することによって残差閾値Dthを低減しても良い。さらに、図7では全てのカテゴリにおいて、準天頂衛星3を捕捉できている場合の残差閾値Dthを、準天頂衛星3を捕捉できていない場合の値よりも小さい値に設定する態様を開示したが、これに限らない。環境判定部16によって判定されているカテゴリが、C1~C2に該当する場合にのみ、準天頂衛星3を捕捉できている場合の残差閾値Dthを、準天頂衛星3を捕捉できていない場合の値よりも小さい値に設定してもよい。
 <変形例1の効果>
 準天頂衛星3を捕捉できている場合には、準天頂衛星3から提供される電離層補正量、対流圏補正量によって、各測位信号についての擬似距離が高精度化される。故に、各捕捉衛星の擬似距離残差自体が全般的に小さくなりうる。このように、準天頂衛星3を捕捉できているか否かによって擬似距離残差の平均レベルが異なるため、仮に準天頂衛星3を捕捉できている場合と捕捉できていない場合とで同一の値を適用してしまうと、準天頂衛星3を捕捉できている場合には、マルチパスに由来する誤差が大きい観測データが測位演算処理に用いる観測データに混ざり込む恐れが高まってしまう。
 対して、本変形例1の構成によれば、準天頂衛星3を捕捉できている場合には捕捉できていない場合よりも残差閾値Dthを小さい値に設定するため、相対的に擬似距離の誤差が大きい観測データが測位演算処理に用いる観測データに混ざり込む恐れをより一層抑制することができる。その結果、測位精度のさらなる向上を期待できる。
 なお、以上では、QZS捕捉判定部173は、準天頂衛星3からの測位信号を受信できているか否かを判定するものとしたが、これに限らない。QZS捕捉判定部173は、準天頂衛星3からの補強信号を受信できているか否かを判定するものとしても良い。
 また、測位装置1が広域通信網等を介して地上局から電離層補正量、対流圏補正量等のローカル補正情報を取得可能に構成されている場合には、閾値変更部172は、ローカル補正情報を取得できているか否かによって上述した変形例1と同様の処理を実施しても良い。つまり、ローカル補正情報を取得できている場合には、ローカル補正情報を取得できていない場合よりも残差閾値Dthを小さい値に設定する。そのような構成によっても変形例1と同様の効果を奏する。
 (変形例2)
 図8に示すように測位装置1が、所定の周波数帯の電波を用いて測位装置1の周辺に存在する他装置と直接的な無線通信(以降、狭域通信)を実施するための通信モジュールである狭域通信部4と相互通信可能に接続されている場合には、狭域通信部4を介して他装置での周辺環境の判定結果を取得し、多数決によって現在の周辺環境のカテゴリを最終的に決定してもよい。そのような構成によれば、測位装置1は、何かしらの理由で自装置のみが周辺環境のカテゴリを誤判定している場合に、当該誤りを修正することができる。つまり、最終的に使用されるカテゴリの精度や信頼度を高めることができる。
 以下、上記の思想に基づく構成を変形例2として説明する。なお、上述の狭域通信とは、換言すれば、広域通信網を介さない通信である。ここでは一例として、測位装置1が車両に搭載されていることを踏まえ、狭域通信部4は車両Hvの周辺に存在する他車両と直接的な(つまり車車間通信)を実施するための通信モジュールとする。
 車車間通信に用いられる周波数帯は、たとえば、760MHz帯である。その他、2.4GHz、5.9GHz帯などを用いることもできる。車車間通信を実現するための通信規格は任意のものを採用することができる。たとえば、IEEE1609等にて開示されているWAVE(Wireless Access in Vehicular Environment)の規格を採用することができる。
 車車間通信用の通信モジュールとしての狭域通信部4は、自車両Hvを中心とする半径数百m以内に存在する他車両と通信可能に構成されている。なお、ここでは一例として狭域通信部4は、通信範囲が数百m程度となる車車間通信規格に準拠した無線通信機能を提供する通信モジュールとするが、これに限らない。他の態様として狭域通信部4は、通信範囲が例えば最大でも数十メートル程度となる所定の近距離無線通信規格に準拠した通信(以降、近距離通信とする)を実施する通信モジュールであっても良い。例えばBluetooth Low Energy(Bluetoothは登録商標)や、Wi-Fi(登録商標)、ZigBee(登録商標)等が上述した近距離無線通信規格に該当する。
 狭域通信部4は、測位装置1(特に環境判定部16)からの要求に基づいて、自装置の環境判定部16が判定したカテゴリを示すカテゴリ情報(例えばカテゴリ番号)を示す判定結果パケットを同報送信する。また、他車両から送信される判定結果パケットを受信すると、当該判定結果パケットに示されるデータを測位装置1に提供する。なお、前提として本変形例2における環境判定部16は、1次フィルタ部15から入力された削除率に基づいた周辺環境のカテゴリの判定を実施する度に、判定結果パケットを生成し、狭域通信部4に出力して同報送信させる。
 判定結果パケットは、周辺環境の判定結果を示す情報(例えばカテゴリ番号)と、測位演算部18の演算結果として位置情報とを含む通信パケットである。また、ここではより好ましい態様として、判定結果パケットには、カテゴリ番号や現在位置情報に加えて、当該通信パケットの送信時刻や、走行速度、進行方向を示す情報が含まれているものとする。
 本変形例における環境判定部16は、狭域通信部4と協働して図9に示すように自車両Hvから所定の距離D1以内に存在する他車両(以降、周辺車両)の判定結果を逐次取得する。また、周辺車両と自車両Hvのそれぞれの判定結果を母集団とした多数決によって周辺環境が何れのカテゴリに該当するかを最終的に決定する多数決部161を備える。
 この多数決部161は、例えば図9に示すように、自分自身は第2カテゴリC2であると判定している状況において周辺車両の判定結果として周辺環境は第1カテゴリC1であるという判定結果が多数派である場合には、最終的に周辺環境は第1カテゴリC1であると判定する。つまり、周辺環境は第2カテゴリC2であるという自分自身の判定結果を破棄して、周辺環境のカテゴリを第1カテゴリであると判定する。
 距離D1は、例えば数十メートル程度(例えば15m)など、車車間通信可能な距離に比べて十分に小さい値に設定されていることが好ましい。自車両Hvから離れている場所を走行している他車両にとっての周辺環境は、自車両Hvの周辺環境とは異なる可能性が高いためである。車車間通信をしている他車両のうち、周辺車両に該当する他車両は、自車両Hvの現在位置と送信元車両の現在位置とを用いて抽出すれば良い。現在位置が、自車両Hvから距離D1以内となっている送信元車両を周辺車両として採用すれば良い。自車両Hvに搭載されている測位装置1が自装置に相当し、周辺車両に搭載されている測位装置1が周辺装置に相当する。
 このような変形例2の構成によれば、冒頭にも記載の通り、測位装置1は、何かしらの理由で自装置のみが周辺環境のカテゴリを誤判定している場合に、当該誤りを修正することができる。その結果、最終的に使用されるカテゴリの精度や信頼度を高めることができる。
 なお、多数決部161は、周辺車両の判定結果と自車両Hvの判定結果の中で最も出現数が多い(換言すれば多数派の)判定結果である最頻カテゴリが、全体の過半数を超えている場合にのみ、周辺環境は多数派の判定結果を採用するようにしてもよい。最頻カテゴリが過半数を超えていないということは、車両(より具体的には測位装置1)毎の判定結果が割れていることを意味し、多数派の判定結果が必ずしも正しい判定結果とは限らないためである。つまり、多数派の判定結果が過半数を超えている場合にのみ多数派の判定結果を採用することで、周辺環境の判定精度を高めることができる。また、自分自身が周辺環境を正しく判定できていたにも関わらず、周辺車両の判定結果に引っ張られて誤った判定結果を最終的なカテゴリとして採用してしまう恐れも低減できる。
 なお、上述した構成では、自車両Hvから距離D1以内に位置する他車両の判定結果を用いて、最終的な周辺環境を決定する態様を開示したが、これに限らない。図10に示すように、車車間通信を実施している車両のうち、自車両Hvの前方であって、且つ、自車両Hvから距離D2以内に存在する他車両の判定結果を用いて(換言すれば母集団として)、最終的な周辺環境を決定してもよい。
 これは次の理由による。カテゴリを決定するための削除率Rは、過去一定時間以内の観測データを元に定まる統計的な指標である擬似距離散布度によって変化する。実際の周辺環境が変化したとしても、擬似距離散布度が所定の散布度閾値を超えるまでは、削除率Rの変化としては表れない。つまり、実際の周辺環境の変化が環境判定部16の判定結果に反映されるまでには擬似距離散布度を算出するためのデータ数に応じた時間がかかる。
 そのような性質を鑑みると、自車両Hvの後方を走行している他車両の判定結果よりも、自車両Hvの前方を走行している他車両の判定結果のほうが、自車両Hvの周辺環境を正しく判定している可能性が高い。故に、自車両Hv前方であって、且つ、自車両Hvから距離D2以内に存在する他車両の判定結果を用いて最終的な周辺環境のカテゴリを決定するように構成することで、カテゴリの判定精度を高めることができる。なお、距離D2は距離D1と同程度値とすればよい。
 また、例えば自車両Hvの前方を走行している他車両のうち、自車両Hvと進行方向を同一とする他車両の判定結果のみと自車両Hvの判定結果を用いて最終的な周辺環境のカテゴリを決定してもよい。
 さらに、高速道路上と当該高速道路の下に配置されている一般道路(以降、高速沿道)上とでは、周辺環境は異なるといえる。当然、高速道路上のほうが上空視界の見晴らしはよく、高速沿道では上空視界の見晴らしは悪い。そのような課題に着眼し、自車両Hvの周辺を走行している他車両のうち、自車両Hvと同じ種別の道路を走行している他車両の自車両Hvの判定結果を用いて最終的な周辺環境のカテゴリを決定することが好ましい。自車両Hvと同じ種別の道路を走行しているか否かは、車車間通信で走行路の種別を共有することによって特定しても良いし、位置情報に含まれる高さ情報に基づいて判別してもよい。また、高速道路と一般道路では走行速度に有意の差が生じるため、走行速度が同じレベル(例えば±10キロ以内)であるか否かによって判別することもできる。
 (変形例3)
 上述した実施形態では誤差指標値として、擬似距離残差ΔDを採用したが、これに限らない。誤差指標値は、測位信号に対するノイズの比率を表すノイズ含有率であってもよい。ノイズ含有率は、雑音の分散を信号の分散で割った値、又は、その対数を取った値とすればよい。ノイズ含有率は値が小さいほど、当該測位信号に基づいて算出された擬似距離が誤差を含んでいる可能性が小さいことを意味するパラメータである。ノイズ含有率は、測位信号のS/N比の逆数として算出されれば良い。
 また、誤差指標値は、90°から仰角を減算することによって定まる仰角残差であってもよい。仰角が90°に近いほど、つまり仰角残差が小さいほど、擬似距離残差ΔDが小さくなる傾向があるためである。なお、仰角の逆数を仰角残差として採用してもよい。仰角残差は仰角が大きいほど小さい値をとるパラメータとすればよい。
 ただし、擬似距離残差ΔDは、ノイズ含有率や仰角残差よりも擬似距離が誤差を含んでいる可能性の大きさを直接的に表す。故に、誤差指標値としては、擬似距離残差ΔDを用いることが好ましい。なお、擬似距離残差ΔD、ノイズ含有率、及び仰角残差といった3つの指標の任意の2つ又は全部を組み合わせて用いて2次フィルタリング処理を実施しても良い。
 本開示に記載されるフローチャート、あるいは、フローチャートの処理は、複数の部(あるいはステップと言及される)から構成され、各部は、たとえば、S101と表現される。さらに、各部は、複数のサブ部に分割されることができる、一方、複数の部が合わさって一つの部にすることも可能である。さらに、このように構成される各部は、サーキット、デバイス、モジュール、ミーンズとして言及されることができる。
 本開示は、実施例に準拠して記述されたが、本開示は当該実施例や構造に限定されるものではないと理解される。本開示は、様々な変形例や均等範囲内の変形をも包含する。加えて、様々な組み合わせや形態、さらには、それらに一要素のみ、それ以上、あるいはそれ以下、を含む他の組み合わせや形態をも、本開示の範畴や思想範囲に入るものである。

Claims (8)

  1.  複数の測位衛星のそれぞれから送信される測位信号を受信する信号受信部(11)と、
     前記信号受信部が受信した測位信号に基づいて、受信した測位信号の送信元としての測位衛星についてのドップラーシフト量及び擬似距離の少なくとも何れか一方を、マルチパス環境下であるか否かを判定するための指標値として逐次算出する環境指標値算出部(12)と、
     前記信号受信部で測位信号を受信できている測位衛星である捕捉衛星毎に、その捕捉衛星に対して前記環境指標値算出部が算出した前記指標値のばらつき度合いを示す指標値散布度を算出する散布度算出部(13)と、
     前記散布度算出部が算出した前記捕捉衛星毎の前記指標値散布度に基づいて、複数の前記捕捉衛星のうち、測位演算処理に使用しない前記捕捉衛星である不使用衛星を決定する1次選択部(15)と、
     前記捕捉衛星の総数に対する前記1次選択部によって不使用衛星に設定された前記捕捉衛星の数の比である削除率に基づいて、周辺環境はマルチパスが発生しやすい環境であるか否かを判定する環境判定部(16)と、
     前記環境判定部の判定結果に基づいて、測位演算処理に用いる前記捕捉衛星の条件を変更する条件変更部(172)と、
     前記条件変更部によって設定されている条件を充足している測位衛星からの測位信号を用いて測位演算処理を実施する測位演算部(18)と、を備える測位装置。
  2.  請求項1に記載の測位装置であって、
     前記環境判定部は、マルチパスの発生のしやすさに応じて段階的に予め設定されている複数のカテゴリのうちの何れに周辺環境は該当するかを、前記削除率に基づいて判定するものであって、
     複数の前記カテゴリには、少なくとも、マルチパスが最も発生しにくいオープンスカイ環境に相当するカテゴリであるオープンスカイカテゴリと、マルチパスが最も発生しやすい都市部に相当するカテゴリである都市部カテゴリとが含まれており、
     前記条件変更部は、前記環境判定部によって判定された周辺環境がよりマルチパスの生じにくいカテゴリであるほど、測位演算処理に用いる前記捕捉衛星の条件を厳しくするように構成されている測位装置。
  3.  請求項2に記載の測位装置であって、
     前記1次選択部によって前記不使用衛星に設定されなかった前記捕捉衛星である1次残存衛星毎に、当該捕捉衛星から送信された測位信号に基づいて算出された擬似距離が誤差を含んでいる可能性の大きさを表す誤差指標値を算出する誤差指標値算出部(171)を備え、
     前記測位演算処理に用いる前記捕捉衛星の条件は、前記誤差指標値が前記条件変更部によって設定されるフィルタリング閾値未満となっていることであって、
     前記条件変更部は、前記環境判定部によって判定された周辺環境がよりマルチパスの生じにくいカテゴリであるほど、前記フィルタリング閾値を小さい値に設定するように構成されている測位装置。
  4.  請求項3に記載の測位装置であって、
     前記誤差指標値は、前記測位装置の位置座標と前記捕捉衛星の位置座標とから定まる前記測位装置から前記捕捉衛星までの直線距離と、当該捕捉衛星の前記擬似距離との差の大きさを表す擬似距離残差である測位装置。
  5.  請求項3又は4に記載の測位装置であって、
     前記フィルタリング閾値の候補値として複数段階の前記候補値が設定されてあって、
     前記条件変更部は、前記環境判定部によって周辺環境は前記オープンスカイカテゴリであると判定されている場合には、前記フィルタリング閾値を、複数の前記候補値の中で最も小さい前記候補値に設定する一方、
     前記環境判定部によって周辺環境は前記都市部カテゴリであると判定されている場合には、前記フィルタリング閾値を、複数の前記候補値の中で最も大きい前記候補値に設定する測位装置。
  6.  請求項3から5の何れか1項に記載の測位装置であって、
     前記フィルタリング閾値の候補値として、複数の前記カテゴリのそれぞれと対応づけられた複数の前記候補値が設定されており、
     複数の前記測位衛星には、前記測位装置が使用される地域にとっての準天頂軌道を周回する準天頂衛星が少なくとも1つ含まれており、
     前記捕捉衛星毎の仰角又は測位信号に示されている送信元情報に基づいて、前記準天頂衛星を捕捉できているかを判定する準天頂捕捉判定部(173)を備え、
     前記条件変更部は、
     前記準天頂捕捉判定部によって前記準天頂衛星からの測位信号を捕捉できていないと判定されている場合には、前記フィルタリング閾値を、前記環境判定部によって判定されている前記カテゴリに対応づけられている前記候補値に設定する一方、
     前記準天頂捕捉判定部によって前記準天頂衛星からの測位信号を捕捉できていると判定されている場合には、前記フィルタリング閾値を、前記環境判定部によって判定されている前記カテゴリに対応づけられている前記候補値よりも小さい値に設定する測位装置。
  7.  請求項2から6の何れか1項に記載の測位装置であって、
     前記測位装置としての自装置から所定距離以内に存在する他の前記測位装置である周辺装置と直接的に通信するための狭域通信部(4)と相互通信可能に接続されており、
     前記環境判定部は、前記狭域通信部と協働して、前記自装置の前記環境判定部が前記削除率に基づいて判定したカテゴリを示すカテゴリ情報を前記周辺装置に同報送信するように構成されており、
     前記周辺装置のそれぞれが判定した周辺環境のカテゴリを前記狭域通信部と協働して取得するとともに、前記周辺装置と前記自装置の判定結果を用いた多数決によって、周辺環境が何れのカテゴリに該当するかを最終的に決定する多数決部(161)をさらに備える測位装置。
  8.  請求項7に記載の測位装置であって、
     前記多数決部は、前記周辺装置の判定結果と前記自装置の判定結果の中で最も出現数が多い判定結果である最頻カテゴリの出現数が過半数を超えている場合には、周辺環境は前記最頻カテゴリに該当すると判定する一方、前記最頻カテゴリの出現数が過半数を超えていない場合には前記自装置の判定結果を最終的な周辺環境として採用する測位装置。
PCT/JP2018/012620 2017-07-07 2018-03-28 測位装置 WO2019008845A1 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
US16/708,707 US11221415B2 (en) 2017-07-07 2019-12-10 Positioning device and positioning method

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2017133972A JP6743777B2 (ja) 2017-07-07 2017-07-07 測位装置
JP2017-133972 2017-07-07

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US16/708,707 Continuation US11221415B2 (en) 2017-07-07 2019-12-10 Positioning device and positioning method

Publications (1)

Publication Number Publication Date
WO2019008845A1 true WO2019008845A1 (ja) 2019-01-10

Family

ID=64950727

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2018/012620 WO2019008845A1 (ja) 2017-07-07 2018-03-28 測位装置

Country Status (3)

Country Link
US (1) US11221415B2 (ja)
JP (1) JP6743777B2 (ja)
WO (1) WO2019008845A1 (ja)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6733619B2 (ja) * 2017-07-07 2020-08-05 株式会社デンソー 測位装置
US20210247521A1 (en) * 2018-07-09 2021-08-12 Sony Semiconductor Solutions Corporation Receiver and control method of the same
CN113189629B (zh) * 2021-04-22 2023-07-07 和芯星通科技(北京)有限公司 一种卫星导航系统中观测值的管理方法
EP4084430A1 (en) * 2021-04-28 2022-11-02 Aptiv Technologies Limited Safety system for vehicle chassis sensors
WO2023234265A1 (ja) * 2022-05-31 2023-12-07 Necソリューションイノベータ株式会社 衛星測位装置、衛星測位方法、及びコンピュータ読み取り可能な記録媒体

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008292454A (ja) * 2007-04-23 2008-12-04 Seiko Epson Corp 測位方法、プログラム、記憶媒体、測位装置及び電子機器
JP2010032235A (ja) * 2008-07-25 2010-02-12 Seiko Epson Corp 測位装置、測位装置の制御方法及び測位装置の制御プログラム
JP2012233882A (ja) * 2011-04-27 2012-11-29 Mediatek Inc 全地球的航法衛星システム用受信機及びその制御方法
JP2016020812A (ja) * 2014-07-11 2016-02-04 トヨタ自動車株式会社 測位補正情報提供装置および衛星測位装置
WO2018009088A1 (en) * 2016-07-04 2018-01-11 Llc "Topcon Positioning Systems" Gnss positioning system and method using multiple processing threads
US20180091962A1 (en) * 2016-09-23 2018-03-29 Apple Inc. Context-dependent allocation of shared resources in a wireless communication interface

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3922585B2 (ja) 2004-05-13 2007-05-30 セイコーエプソン株式会社 測位装置、測位方法、測位プログラム、測位プログラムを記録したコンピュータ読み取り可能な記録媒体
US7898474B2 (en) * 2006-07-20 2011-03-01 Seiko Epson Corporation Positioning device, method of controlling positioning device, and recording medium having program for controlling positioning device recorded thereon
JP4910676B2 (ja) * 2006-12-15 2012-04-04 セイコーエプソン株式会社 演算回路、測位装置、指標算出方法及びプログラム
JP5186873B2 (ja) * 2007-10-10 2013-04-24 セイコーエプソン株式会社 測位方法、プログラム、測位装置及び電子機器
JP5186874B2 (ja) * 2007-10-10 2013-04-24 セイコーエプソン株式会社 測位方法、プログラム、測位装置及び電子機器
JP5228433B2 (ja) * 2007-10-15 2013-07-03 セイコーエプソン株式会社 測位方法、プログラム、測位装置及び電子機器
JP5347443B2 (ja) * 2008-11-11 2013-11-20 セイコーエプソン株式会社 位置算出方法及び位置算出装置
JP5994237B2 (ja) 2011-11-24 2016-09-21 株式会社豊田中央研究所 測位装置及びプログラム
JP6201762B2 (ja) * 2014-01-08 2017-09-27 株式会社デンソー 速度推定装置
JP6409346B2 (ja) * 2014-06-04 2018-10-24 株式会社デンソー 移動距離推定装置
US9798017B2 (en) * 2015-04-09 2017-10-24 Topcon Positioning Systems, Inc. Reducing time and increasing reliability of ambiguity resolution in GNSS
AU2015409237B2 (en) * 2015-09-17 2019-07-18 Mitsubishi Electric Corporation Positioning satellite selection device, positioning information transmitting device, and positioning system
JP6733619B2 (ja) * 2017-07-07 2020-08-05 株式会社デンソー 測位装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008292454A (ja) * 2007-04-23 2008-12-04 Seiko Epson Corp 測位方法、プログラム、記憶媒体、測位装置及び電子機器
JP2010032235A (ja) * 2008-07-25 2010-02-12 Seiko Epson Corp 測位装置、測位装置の制御方法及び測位装置の制御プログラム
JP2012233882A (ja) * 2011-04-27 2012-11-29 Mediatek Inc 全地球的航法衛星システム用受信機及びその制御方法
JP2016020812A (ja) * 2014-07-11 2016-02-04 トヨタ自動車株式会社 測位補正情報提供装置および衛星測位装置
WO2018009088A1 (en) * 2016-07-04 2018-01-11 Llc "Topcon Positioning Systems" Gnss positioning system and method using multiple processing threads
US20180091962A1 (en) * 2016-09-23 2018-03-29 Apple Inc. Context-dependent allocation of shared resources in a wireless communication interface

Also Published As

Publication number Publication date
US20200110182A1 (en) 2020-04-09
JP2019015636A (ja) 2019-01-31
US11221415B2 (en) 2022-01-11
JP6743777B2 (ja) 2020-08-19

Similar Documents

Publication Publication Date Title
US11221415B2 (en) Positioning device and positioning method
US9417327B2 (en) Selection method of satellites for RTK positioning calculation and a selection device of satellites for the same
US11275180B2 (en) Positioning device and positioning method
CN109313272B (zh) 使用速度积分的改进gnss接收器
CN105353391A (zh) 一种支持多类型定位终端的多网融合定位增强系统及方法
US20080180315A1 (en) Methods and systems for position estimation using satellite signals over multiple receive signal instances
CN111913200B (zh) 车辆群体差分定位方法、rsu设备、融合感知设备及系统
CN112327340B (zh) 终端定位精度评估方法、装置、设备以及介质
KR102188880B1 (ko) 단말, 기준국 및 위치 측위 방법
WO2015145718A1 (ja) 測位装置
JP4723932B2 (ja) 測位システム
CN113050142B (zh) 终端设备的定位方法、装置、电子设备及可读存储介质
US8638258B2 (en) Method and system for a virtual wide area GNSS reference network
CN115792980A (zh) 一种模型和数据双重驱动的gnss rtk定位选星方法及系统
Yuan et al. GPS multipath and NLOS mitigation for relative positioning in urban environments
US7184887B2 (en) Method and apparatus for calculating a figure of merit for GPS position using NMEA 0183 output
CN113093250B (zh) 伪距观测数据修复方法、定位信息确定方法及装置
Socharoentum et al. A machine learning approach to detect non-line-of-sight GNSS signals in Nav2Nav
JP2016020812A (ja) 測位補正情報提供装置および衛星測位装置
CN112904390B (zh) 定位方法、装置、计算机设备和存储介质
CN116626726A (zh) 确定关于车辆的基于gnss的定位的完整性信息的方法
CN113835109B (zh) 终端的定位方法、装置、电子设备、存储介质及程序产品
Tijero et al. Advanced GNSS Algorithms for Safe Autonomous Vehicles
Pasnikowski et al. Challenges for integrity in navigation of high precision
US20240085572A1 (en) Method for providing gnss-relevant route information

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: 18827337

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 18827337

Country of ref document: EP

Kind code of ref document: A1