WO2014132618A1 - 状態検出方法、補正値処理装置、測位システム、および状態検出プログラム - Google Patents

状態検出方法、補正値処理装置、測位システム、および状態検出プログラム Download PDF

Info

Publication number
WO2014132618A1
WO2014132618A1 PCT/JP2014/000964 JP2014000964W WO2014132618A1 WO 2014132618 A1 WO2014132618 A1 WO 2014132618A1 JP 2014000964 W JP2014000964 W JP 2014000964W WO 2014132618 A1 WO2014132618 A1 WO 2014132618A1
Authority
WO
WIPO (PCT)
Prior art keywords
reference station
master
master reference
correction value
positioning position
Prior art date
Application number
PCT/JP2014/000964
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 日本電気株式会社
Priority to CN201480010469.1A priority Critical patent/CN105008956A/zh
Priority to US14/764,476 priority patent/US20150362596A1/en
Priority to JP2015502766A priority patent/JP6119838B2/ja
Publication of WO2014132618A1 publication Critical patent/WO2014132618A1/ja

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/03Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
    • G01S19/07Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing data for correcting measured positioning data, e.g. DGPS [differential GPS] or ionosphere corrections
    • 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/40Correcting position, velocity or attitude
    • G01S19/41Differential correction, e.g. DGPS [differential GPS]

Definitions

  • the present invention relates to a state detection method, a correction value processing device, a positioning system, and a state detection program, and more particularly to a method for detecting an abnormal state of a reference station in a differential GPS having a plurality of reference stations.
  • GPS Global Positioning System
  • differential GPS in which positioning accuracy by GPS is improved by installing a reference station fixed on the ground surface, has become widespread.
  • the positioning accuracy depends on the time and position accuracy of the reference station. For example, if an error occurs in the operation of the reference station, the GPS user receiver that affects the correction signal, receives the correction signal, corrects the signal received from the GPS satellite, and performs positioning, positions itself different from the correct position. It is estimated that
  • the requirement established as a GBAS (Ground Based Augmentation System) standard for differential GPS reference stations currently in operation is 10 ⁇ 5 integrity risk per 150 seconds.
  • the integrity risk is the probability of loss of integrity.
  • This requirement is established by ICAO (International Civil Aviation Organization) / RTCA (Radio Technical Commission for Aeronautics).
  • ICAO International Civil Aviation Organization
  • RTCA Radio Technical Commission for Aeronautics
  • category I In operation of the GBAS, category I (CAT-I) In the numerical integrity risk is 10-7 as a requirement may be ignored, category integrity risk is 10 -9 II and III (CAT-II / In III), the size cannot be ignored.
  • Patent Document 1 discloses a method for detecting a reference station in which an abnormality has occurred in a differential GPS (DGPS) using a plurality of reference stations.
  • Receiver / satellite specific differential correction adjusted to remove clock bias and reflect common measurement time when identifying N GPS satellites and K reference stations with indices n and k respectively. Let C n, k be the value.
  • the receiver / satellite specific operation correction value is generated by the differential correction processing means based on the GPS signal received from the GPS satellite. That is, first, the satellite-specific differential correction value is generated based on the received GPS signal and the known survey position of the GPS signal receiver of the reference station fixed on the ground.
  • GPS signal receivers have a clock time offset or bias with respect to the GPS satellite atomic clock time, and each receiver / satellite specific measurement has a different measurement time.
  • the GPS signal receiver generates a differential correction value without a clock bias from the receiver / satellite specific measurement value.
  • the GPS signal receiver then adjusts the receiver / satellite specific correction value according to the common synchronization time. Through these processes, the receiver / satellite specific differential correction value C n, k adjusted to reflect the common measurement time is obtained by removing the clock bias.
  • the receiver / satellite specific discrimination value Zn , k is used as an amount used for the verification.
  • the reference station is determined based on a predetermined detection threshold DT n, k . That is, the judgment formula:
  • > DT n, k To determine the reference station. Specifically, for a given reference station k, if all the GPS satellites (n 1,..., N) have a judgment formula:
  • Patent Documents 2 to 5 disclose related technologies.
  • Patent Document 1 can detect a reference station where an abnormality has occurred, but has the following problems.
  • the verification based on the above judgment formula is performed only in the range domain (measurement of the distance from the GPS user receiver (aircraft or the like) to the GPS satellite). Therefore, if the absolute value of the receiver / satellite specific discrimination value is less than the threshold value, it is not considered that an abnormality has occurred.
  • the sign of the receiver / satellite specific discrimination value is biased to the negative side or the positive side, a fatal position error may occur in the position domain (determination of the position of the GPS user receiver).
  • the receiver / satellite specific discrimination value used in the verification accurately reflects the determination of the presence or absence of abnormality. I can't say that. For example, for a reference station having an abnormality, the receiver / satellite specific discrimination value becomes small, and a detection failure may occur. In addition, the receiver / satellite specific discrimination value becomes large with respect to the reference station having no abnormality, and there is a possibility that erroneous detection occurs.
  • a state detection method of the present invention is a state detection method for detecting an abnormal state of a reference station in a positioning system having a satellite and a reference station, and is based on a signal received by the reference station from the satellite. Determining the first positioning position of the reference station, selecting a predetermined number of master reference stations in order from the smallest distance between the position information given in advance of the reference station and the first positioning position of the reference station, Based on the signal received by the master reference station from the satellite, a correction value for the master reference station is generated, and the reference station other than the master reference station determines the first positioning position of the reference station other than the master reference station as the correction of the master reference station.
  • the second positioning position of the reference station other than the master reference station is obtained by correcting with the value, and the reference other than the master reference station is obtained. Position by pre-given position information, a distance between the second measured position when larger than a first threshold value, the reference station other than the master reference station and determines that an abnormal state.
  • the state detection method of the present invention is a state detection method for detecting an abnormal state of the plurality of reference stations in a positioning system having a satellite and a plurality of reference stations, and is based on a signal received by the reference station from the satellites. Determining a first positioning position of the reference station, selecting a reference station in which a distance between the position information given in advance of the reference station and the first positioning position of the reference station is included in a predetermined range as a master reference station, Based on the signal received by the master reference station from the satellite, a correction value for the master reference station is generated, and the reference station other than the master reference station determines the first positioning position of the reference station other than the master reference station as the correction of the master reference station.
  • the second positioning position of the reference station other than the master reference station is obtained by correcting with the value, and the base position other than the master reference station is obtained.
  • a reference station other than the master reference station is determined to be in an abnormal state when a distance between a position based on position information given in advance of the station and the second positioning position is greater than a first threshold value. .
  • the correction value processing apparatus of the present invention includes a data receiving unit that receives a signal received from a satellite by a reference station, and a first positioning that calculates a first positioning position of the reference station based on the signal received by the data receiving unit.
  • Position calculating means including master reference station selecting means for selecting a predetermined number of master reference stations in order from the smallest distance between the position information given in advance of the reference station and the first positioning position of the reference station; and the master
  • a master reference station correction value calculating means for generating a correction value of the master reference station based on a signal received by the reference station from the satellite, and the reference station other than the master reference station is configured to determine the first positioning position of the reference station other than the master reference station.
  • the second positioning position of a reference station other than the master reference station is calculated by correcting with the correction value of the master reference station Second positioning position calculation means, and when the distance between the position based on the position information given in advance by a reference station other than the master reference station and the second positioning position is greater than a first threshold, An abnormality determining means for determining that the reference station is in an abnormal state.
  • the correction value processing apparatus of the present invention includes a data receiving unit that receives a signal received from a satellite by a reference station, and a first positioning that calculates a first positioning position of the reference station based on the signal received by the data receiving unit.
  • a master reference station correction value calculating means for generating a correction value of the master reference station based on a signal received by the reference station from the satellite, and the reference station other than the master reference station is configured to determine the first positioning position of the reference station other than the master reference station.
  • the second positioning position of a reference station other than the master reference station is determined.
  • the distance between the second positioning position calculating means and the position based on the position information given in advance of the reference station other than the master reference station and the second positioning position is larger than the first threshold value, other than the master reference station
  • an abnormality determining means for determining that the reference station is in an abnormal state.
  • the positioning system of the present invention is a positioning system including a satellite, a reference station, and a correction value processing device, and the correction value processing device is configured to perform a first operation of the reference station based on the signal received by the data receiving means.
  • a predetermined number of master reference stations are selected in ascending order of distance between first positioning position calculation means for calculating a positioning position, position information given in advance of the reference station, and the first positioning position of the reference station.
  • Master reference station selection means master reference station correction value calculation means for generating a correction value for the master reference station based on the signal received by the master reference station from the satellite, and the reference station other than the master reference station is the master reference station other than the master reference station
  • the distance between the second positioning position calculating means for calculating the second positioning position of the sub-station, the position based on the position information given in advance by the reference station other than the master reference station, and the second positioning position is less than the first threshold value.
  • an abnormality determining means for determining that a reference station other than the master reference station is in an abnormal state when it is larger.
  • the positioning system of the present invention is a positioning system including a satellite, a reference station, and a correction value processing device, and the correction value processing device is configured to perform a first operation of the reference station based on the signal received by the data receiving means.
  • First reference position calculating means for calculating a positioning position, a reference station in which a distance between the position information given in advance of the reference station and the first positioning position of the reference station is included in a predetermined range is selected as a master reference station.
  • Master reference station selection means master reference station correction value calculation means for generating a correction value for the master reference station based on the signal received by the master reference station from the satellite, and the reference station other than the master reference station is the master reference station other than the master reference station By correcting the first positioning position of the reference station with the correction value of the master reference station,
  • a second positioning position calculating means for calculating a second positioning position of the reference station, a position based on position information given in advance by a reference station other than the master reference station, and a distance between the second positioning position than the first threshold value
  • an abnormality determining means for determining that a reference station other than the master reference station is in an abnormal state when it is larger.
  • the storage medium of the present invention is a storage medium storing a program for detecting an abnormal state of the reference station in a positioning system having a satellite and a reference station, and based on the signal received by the data receiving means, the reference station A process of calculating the first positioning position of the reference station, a process of selecting a predetermined number of master reference stations in order from the smallest distance between the position information given in advance of the reference station and the first positioning position of the reference station; A process for generating a correction value for the master reference station based on a signal received from the satellite by the master reference station; and a reference station other than the master reference station determines the first positioning position of the reference station other than the master reference station as the master A second position for calculating a second positioning position of a reference station other than the master reference station is corrected by the correction value of the reference station.
  • a reference station other than the master reference station is abnormal when the distance between the position based on position position calculation processing and a position based on position information given in advance by a reference station other than the master reference station and the second positioning position is greater than a first threshold. And a process for determining that it is in a state.
  • the storage medium of the present invention is a storage medium storing a program for detecting an abnormal state of the reference station in a positioning system having a satellite and a reference station, and based on the signal received by the data receiving means, the reference station A process of calculating the first positioning position of the reference station, a process of selecting, as a master reference station, a reference station in which a distance between the position information given in advance of the reference station and a distance from the first positioning position of the reference station is included in a predetermined range; A process for generating a correction value for the master reference station based on a signal received from the satellite by the master reference station; and a reference station other than the master reference station determines the first positioning position of the reference station other than the master reference station as the master The second positioning position of a reference station other than the master reference station is calculated by correcting with the correction value of the reference station.
  • the reference station other than the master reference station When the distance between the two positioning position calculation process and the position based on the position information given in advance by the reference station other than the master reference station and the second positioning position is larger than the first threshold, the reference station other than the master reference station And a process for determining that the state is abnormal.
  • the abnormal state of the reference station in the satellite positioning system can be detected with high accuracy.
  • An example of the structure of the differential GPS which concerns on the 1st Embodiment of this invention is shown. It is a flowchart which shows an example of the procedure of the state detection method of the reference station of differential GPS which concerns on the 1st Embodiment of this invention. An example of the structure of the correction value processing apparatus which concerns on the 1st Embodiment of this invention is shown. An example of a structure of the data processing part which concerns on the 2nd Embodiment of this invention is shown. It is a flowchart which shows an example of the procedure of the state detection method of the reference station of differential GPS which concerns on the 2nd Embodiment of this invention. It is a flowchart which shows an example of the procedure of the state detection method of the reference station of differential GPS which concerns on the 3rd Embodiment of this invention.
  • FIG. 1 shows an example of the configuration of a differential GPS 101 according to the first embodiment of the present invention.
  • the receiving device 601 receives GPS signals from the GPS satellites 201 to 20n and estimates its own position.
  • the receiving device 601 is mounted on an aircraft, for example.
  • Each of the plurality of reference stations 301 to 304 receives GPS signals from the GPS satellites 201 to 20n, and transmits observation data to the correction value processing device 401.
  • the correction value processing device 401 calculates a correction value for positioning by receiving GPS signals from the GPS satellites 201 to 20n, and sends the correction value to the broadcast transmission device 501.
  • the broadcast transmission device is, for example, a VDB (VHF (Very High Frequency) Data Broadcast) transmission device.
  • VDB VHF (Very High Frequency) Data Broadcast
  • the receiving device 601 receives the correction value broadcast transmitted from the broadcast transmitting device 501, corrects the GPS signal received from the GPS satellites 201 to 20n, and estimates its position.
  • FIG. 3 shows an example of the configuration of the correction value processing apparatus 401 according to the present embodiment.
  • the correction value processing apparatus 401 includes a data reception unit 402 that receives observation data transmitted from the reference stations 301 to 304, and a correction value transmission unit 405 that transmits correction value information to the broadcast transmission apparatus 501 to provide the reception apparatus 601 with correction value information.
  • the correction value processing apparatus further includes a data processing unit 403 that processes the observation data received by the data receiving unit 402 to generate correction value information, and a data holding unit 404 that holds data regarding the reference stations 301 to 304.
  • the data holding unit 404 further temporarily holds the information processed by the data processing unit 403, and holds parameter values and the like required for the verification processing described below.
  • the number of reference stations is four, but the present invention is not limited to this.
  • the number of reference stations may be less than 4 or more than 4.
  • the correction value processing device 401 receives observation data from the reference stations 301 to 304 that have received GPS signals emitted from the GPS satellites 201 to 20n (step S201 in FIG. 2).
  • the correction value processing apparatus 401 measures (single positioning) the position of each of the reference stations 301 to 304 based on the observation data (step S202).
  • the single positioning is executed by using a pseudo distance from the four GPS satellites (distance from the GPS satellite to the receiver).
  • the clock built in the receiving device 601 is inferior in accuracy by several orders of magnitude compared to the atomic clock of a GPS satellite. Even if the transmission time of the radio wave carved by the atomic clock of the GPS satellite is accurate, the arrival time measured by the receiver may include a large error.
  • pseudo distances from four or more GPS satellites are used simultaneously. .
  • the positioning position of the reference station can be obtained by independent positioning. At this time, at least one of the following conditions may be imposed so that the positioning accuracy does not vary among the reference stations.
  • the correction value processing apparatus 401 selects a master reference station and a verification target reference station from the reference stations 301 to 304 (step S203).
  • Three-dimensional coordinate values are given in advance as high-precision numerical values for the positions of the reference stations fixed on the respective ground surfaces. Hereinafter, this is called a survey position. Coordinate values for the survey position are held in the data holding unit 404 of the correction value processing device 401.
  • positioning error is defined as a three-dimensional distance.
  • the positioning position obtained by independent positioning is compared with the surveying position, and the reference station with the smallest positioning error and the next smallest reference station are selected.
  • the reference station with the smallest positioning error is designated as master reference station A, and the next smallest reference station is designated as master reference station B.
  • two master reference stations are selected.
  • the number of master reference stations may be three or more.
  • the base station that is not selected as the master base station is the base station to be verified.
  • two reference stations are the test target reference station X and the test target reference station Y, respectively.
  • the correction value processing device 401 calculates the DGPS correction value at the master reference station based on the observation data obtained from the master reference station and the value of the surveying position (step S204).
  • the DGPS correction value at the master reference station A becomes the DGPS correction value A
  • the DGPS correction value at the master reference station B becomes the DGPS correction value B.
  • the DGPS correction value is obtained as follows.
  • the pseudo distance measured by the reference station includes an error due to the true distance between the GPS satellite and the reference station.
  • the pseudo distance includes an error due to internal noise such as a clock bias of the reference station (a difference between a time displayed by the reference station and a true time) or a clock bias of a GPS satellite.
  • the pseudorange includes errors due to tropospheric delay experienced by the reference station, ionospheric delay experienced by the reference station, multipath experienced by the reference station, and internal noise such as thermal noise experienced by the reference station.
  • the pseudo-range correction value is calculated from the position information obtained by surveying the position information of the GPS satellites broadcasted by the GPS satellites and the reference station, and from the geometric distance between the GPS satellites and the reference stations. It is defined as the value obtained by subtracting the pseudorange.
  • the geometric distance between the GPS satellite and the reference station is the true distance between the GPS satellite and the reference station, the error in the position information of the GPS satellite broadcast by the GPS satellite, and the distance obtained by surveying the reference station. The error is added.
  • the pseudo-range correction value derived from the above relational expression includes the satellite clock bias, the tropospheric delay amount experienced by the satellite, the ionospheric delay amount experienced by the satellite, and the error that can be recognized and separated as the error of the receiver. , Including other errors.
  • the obtained pseudo distance correction value by the reference station is applied to the correction of the pseudo distance of the receiving device.
  • the pseudo-range of the receiver is the true distance between the GPS satellite and the receiver, the error component that can be separated during the positioning calculation as the clock bias of the receiver, and the standard deviation that is unique or assumed for each observation. Contains the error component evaluated using.
  • the actual pseudorange correction value is averaged at multiple reference stations. Thereby, random errors such as multipath and noise of the receiving apparatus are suppressed. On the other hand, the ionospheric delay and tropospheric delay remain averaged.
  • step S205 the master base station selection verification is executed (step S205).
  • the two DGPS correction values are applied between the two master reference stations to perform DGPS positioning, and the positioning results are compared with the respective survey positions. Since the selection of the master reference station depends on the single positioning as described above, there is a possibility that the positioning error may be reduced accidentally even though there is an abnormality. The above verification of the selection of the master reference station eliminates this possibility.
  • the positioning position obtained from the result of independent positioning accurately reproduces the survey position of the reference station. That is, the reference station having a large positioning error is suspected of an operation failure.
  • the correction value generated from the reference station in which the failure has occurred may correct the GPS signal in a direction deviating from the correct position, and if included in the correction value information provided to the receiving device 601, it is a factor that increases the error. Become.
  • the correction value processing device 401 corrects the selection of the master reference station (step S207) and calculates the correction value of the newly selected master reference station. This is performed (step S204).
  • the correction value processing device 401 performs DGPS positioning by applying the DGPS correction value of the master reference station to the verification target reference station (step S208).
  • step S209 the verification of the verification target base station is executed.
  • the difference between the DGPS positioning result of the verification target reference station and the survey position of each reference station is obtained.
  • the difference value exceeds a predetermined threshold with respect to the DGPS correction values of the two master reference stations, it is determined that the verification target reference station is abnormal.
  • Table 1 shows the logic for determining pass / fail.
  • the threshold used for the above determination may be changed depending on the positioning accuracy.
  • the threshold value may be determined based on PDOP (Position Dilution of Precision) reflecting the arrangement status of GPS satellites and the number of satellites. This is because PDOP may change due to increase / decrease in the number of satellites, and positioning accuracy may change, and the above determination may be inaccurate with a fixed threshold value.
  • PDOP Position Dilution of Precision
  • the evaluation is performed based on the three-dimensional distance.
  • the Maxwell-Boltzmann distribution is applied to the test process based on this three-dimensional distance.
  • the DGPS correction values from each of a plurality of reference stations are averaged for the reference stations for the GPS signal from one GPS satellite. Then, the obtained average value is given to the receiving device.
  • the receiving device By averaging the DGPS correction values from a plurality of normally operating reference stations and providing them to the receiving device, variations for each reference station are suppressed, and the receiving device can determine its position with higher accuracy.
  • the data processing unit 403 of the correction value processing apparatus 401 generates correction value information by excluding observation data from the verification target reference station that is rejected as a result of the above-described verification. Then, the correction value transmission unit 405 transmits the correction value information to the broadcast transmission device 501.
  • the correction value processing apparatus 401 performs single positioning with respect to the reference stations 301 to 304. As described above, GPS signals from GPS satellites that may cause errors in positioning are excluded. For the single positioning, GPS signals from n GPS satellites 201 to 20n are used. The three-dimensional coordinate value of the positioning position of the reference station obtained by independent positioning is described as follows.
  • Reference station 301 ( xsa1 , ysa1 , zsa1 ); Reference station 302: (x sa2 , y sa2 , z sa2 ); Reference station 303: (x sa3 , y sa3 , z sa3 ); Reference station 304: ( xsa4 , ysa4 , zsa4 ).
  • the three-dimensional coordinate values of the surveying positions with respect to the reference stations 301 to 304 are described as follows.
  • Reference station 301 (x 1 , y 1 , z 1 ); Reference station 302: (x 2 , y 2 , z 2 ); Reference station 303: (x 3 , y 3 , z 3 ); Reference station 304: (x 4 , y 4 , z 4 ).
  • the difference between the single positioning position and the surveying position for each of the reference stations 301 to 304 is obtained as a three-dimensional distance.
  • Reference station 303: ⁇ sa3 ⁇ ((x sa3 ⁇ x 3 ) 2 + (y sa3 ⁇ y 3 ) 2 + (z sa3 ⁇ z 3 ) 2 );
  • Reference station 304: ⁇ sa4 ⁇ (( x sa4 -x 4) 2 + (y sa4 -y 4) 2 + (z sa4 -z 4) 2).
  • the reference station that gives the smallest difference and the reference station that gives the next smallest difference become master reference stations. That is, the reference station 302 becomes the master reference station A, and the reference station 304 becomes the master reference station B. A reference station that gives a larger difference than these is a test target reference station. In this embodiment, the reference station 303 becomes the verification target reference station X, and the reference station 301 becomes the verification target reference station Y.
  • the DGPS correction value obtained in the master reference station based on the GPS signals from the n GPS satellites 201 to 20n used for the single positioning is described as follows.
  • Positioning position by C A being tested reference station X (x dgpsXA, y dgpsXA , z dgpsXA); Positioning position by C B of test target reference station X: (x dgpsXB , y dgpsXB , z dgpsXB ); Positioning position by C A being tested reference station Y: (x dgpsYA, y dgpsYA , z dgpsYA); Positioning position by C B of test target reference station Y: (x dgpsYB , y dgpsYB , z dgpsYB ).
  • the obtained coordinate values are used for the verification of the reference station described below.
  • the respective DGPS positioning values are executed between the master reference station A and the master reference station B by applying the respective DGPS correction values.
  • the obtained positioning position is described as follows.
  • Positioning position by C B master reference station A (x dgpsAB, y dgpsAB , z dgpsAB); Positioning position by C A master reference station B: (x dgpsBA, y dgpsBA , z dgpsBA).
  • the difference between the positioning position and the surveying position for the master reference station obtained by the above processing is compared with a predetermined threshold value TH master . That is, Master reference station A: If ⁇ dgpsA ⁇ TH master , master reference station A passes; Master base station B: If ⁇ dgpsB ⁇ TH master , master base station B passes.
  • the threshold TH master may be appropriately determined so as to reflect the accuracy of positioning.
  • the threshold TH master may be determined based on the number of GPS satellites and PDOP (Position Dilution of Precision) reflecting the arrangement status of GPS satellites. Good.
  • r is a non-dimensional quantity obtained by dividing a three-dimensional error by a length typical for a target system.
  • the Maxwell-Boltzmann distribution which is a probability distribution function for the r, expressed by the following equation (1):
  • the threshold value TH master may be determined based on the maximum allowable false alarm probability, the maximum allowable detection failure probability, the standard deviation of the error at the time of fault free, and the like.
  • the master base station fails due to the above verification, the master base station that has failed and the base station to be tested are replaced in order to select an appropriate master base station.
  • This operation is processing for eliminating the possibility that the selection of the master reference station is based on the result of the single positioning, and the difference between the single positioning and the surveying position may be reduced by chance.
  • both master base stations fail, both master base stations and the two test target base stations are replaced.
  • This embodiment shows a DGPS configuration including four master reference stations. When the number of master reference stations is larger than 4, the test target reference station having the smallest difference between the result of independent positioning and the survey position and the next smallest test are shown. The target reference station is replaced with two master reference stations.
  • the DGPS correction value calculated based on the observation data also includes an error.
  • C A a DGPS correction value by the master reference station A. That is, the error included in the observation data from the master reference station A is reflected in the positioning position by C B of the master reference station A, and the DGPS correction value C A by the master reference station A is reflected in the positioning position by C A of the master reference station B. The error contained in is reflected.
  • the master base station A is replaced with one of the base stations to be tested, and the master base station is rejected again, the master base station A is returned.
  • the master reference station B may be replaced.
  • the difference between the DGPS positioning result and the survey position at the verification target reference station is calculated as a three-dimensional distance.
  • ⁇ dgpsXA ⁇ ((x dgpsXA ⁇ x 3 ) 2 + (y dgpsXA ⁇ y 3 ) 2 + (z dgpsXA ⁇ z 3 ) 2 );
  • ⁇ dgpsXB ⁇ ((x dgpsXB ⁇ x 3 ) 2 + (y dgpsXB ⁇ y 3 ) 2 + (z dgpsXB ⁇ z 3 ) 2 );
  • ⁇ dgpsYA ⁇ ((x dgpsYA ⁇ x 1 ) 2 + (y dgpsYA ⁇ y 1 ) 2 + (z dgpsYA ⁇ z 1 ) 2 );
  • test target base station X fails, otherwise it passes, If ( ⁇ dgpsYA ⁇ TH test ) and ( ⁇ dgpsYB ⁇ TH test ), the test target base station Y is rejected, and otherwise it is determined to be acceptable.
  • the threshold value TH test may be appropriately determined so as to reflect the accuracy of positioning as in the case of TH master. You may decide to.
  • r is a non-dimensional quantity obtained by dividing a three-dimensional error by a length typical for a target system.
  • the threshold value TH based on the maximum allowable false alarm probability, the maximum allowable detection failure probability, the error standard deviation at the time of fault free, etc. The test may be determined.
  • the correction value information transmitted from the correction value processing device 401 to the broadcast transmission device 501 includes the DGPS correction value of the reference station that has passed the test.
  • the occurrence of the abnormal state of the reference station that receives the GPS signal from the GPS satellite and corrects the GPS signal based on its own position information is highly accurate. It can be detected.
  • the reference base station that is referred to is a base station that operates normally.
  • an abnormal state can be detected based on observation data from a reference station that is determined to operate normally for a reference station that is likely to cause an abnormality. For this reason, even if an abnormality occurs in a plurality of reference stations, the abnormal state can be detected with high accuracy.
  • the correction value information in the differential GPS is used to detect the abnormal state of the reference station, but the present invention is not limited to this.
  • a positioning system that receives a signal from a signal source that periodically transmits its own time and position information and estimates the position of its own device, it is referred to in the process of generating correction information of the signal.
  • the present invention is preferably applied to detection of an abnormal state of a reference station.
  • the receiving device is mounted on an aircraft.
  • the present invention is not limited to this, and the receiving device may be mounted on a ship, or may be held by an automobile or a pedestrian.
  • the correction value information is provided to the receiving device as a VHF broadcast radio wave.
  • the correction value information is not limited to this, and may be radio waves in other frequency bands or may be transmitted by wire.
  • the present invention can be applied to an underwater positioning system.
  • the apparatus may be configured. good.
  • the computer-readable recording medium include a flexible disk, a CD-ROM (Compact Disk Read-Only Memory), and an MO (Magneto-Optical disk).
  • the configuration of the differential GPS 101 in the second embodiment of the present invention is the same as the configuration example of the differential GPS 101 in the first embodiment shown in FIG.
  • a configuration example of the correction value processing apparatus 401 in the second embodiment of the present invention is the same as the configuration example of the correction value processing apparatus 401 in the first embodiment shown in FIG.
  • FIG. 4 is a diagram illustrating a configuration example of the data processing unit 403 included in the correction value processing apparatus 401 according to the second embodiment of the present invention.
  • the data processing unit 403 includes a first positioning position calculation unit 406, a master reference station selection unit 407, a master reference station correction value calculation unit 408, a second positioning position calculation unit 409, and an abnormality determination unit 410.
  • the first positioning position calculation unit 406 obtains the first positioning position of the reference station based on the signal received from the satellite by the reference station.
  • the first positioning position calculation unit 406 performs independent positioning with respect to the reference stations 301 to 304.
  • the first positioning position calculation unit 406 uses GPS signals from the n GPS satellites 201 to 20n as the single positioning, and determines the positioning positions of the reference stations 301 to 304. The three-dimensional coordinate value of is obtained.
  • the master reference station selection unit 407 has a predetermined number of master reference stations in order from the smallest distance between the position information (surveying position) given in advance of the reference station and the first positioning position obtained by the first positioning position calculation unit 406. Select. As described in the first embodiment, the master reference station selection unit 407 obtains the difference between the positioning position and the surveying position for each of the reference stations 301 to 304 as a three-dimensional distance, and calculates the magnitude relationship for the difference. Next, based on the calculated magnitude relationship, the master reference station selection unit 407 selects, for example, the reference station that gives the smallest difference and the reference station that gives the next smallest difference as the master reference station.
  • the master reference station correction value calculation unit 408 generates a correction value for the master reference station based on the signal received by the master reference station from the satellite.
  • the correction value of the master reference station is, for example, the DGPS correction value in the first embodiment.
  • the second positioning position calculation unit 409 allows the reference station other than the master reference station to correct the first positioning position of the reference station other than the master reference station with the correction value of the master reference station, thereby determining the second positioning position of the reference station other than the master reference station. Ask.
  • the second positioning position calculation unit 409 applies each correction value of the master reference station to the first positioning position of the reference station other than the master reference station to obtain the second positioning position. That is, the second positioning position calculation unit 409 obtains the second positioning positions by the number of master reference stations.
  • the abnormality determination unit 410 is in an abnormal state when the distance between the position based on the position information given in advance by the reference station other than the master reference station and the second positioning position is larger than the first threshold. Is determined.
  • the abnormality determination unit 401 determines that a reference station other than the master reference station is abnormal, for example, based on the pass / fail determination logic of Table 1 in the first embodiment.
  • FIG. 5 is a flowchart showing an operation example of the correction value processing apparatus 401 in the second embodiment of the present invention.
  • the data receiving unit 402 of the correction value processing apparatus 401 receives observation data from the reference stations 301 to 304 that have received the GPS signals emitted from the GPS satellites 201 to 20n (step S501).
  • the data processing unit 403 of the correction value processing device 401 measures the position of each of the reference stations 301 to 304 (single positioning) based on the observation data, and obtains a first positioning distance (step S502).
  • the data processing unit 403 selects a predetermined number of master reference stations in order from the smallest distance between the position information (surveying position) given in advance of the reference station and the first positioning position (step S503).
  • the data processing unit 403 calculates a correction value for the master reference station based on the observation data obtained from the master reference station and the value of the surveying position (step S504).
  • the data processing unit 403 obtains the second positioning position of the reference station other than the master reference station by the reference station other than the master reference station correcting the first positioning position of the reference station other than the master reference station with the correction value of the master reference station ( Step S505).
  • the data processing unit 403 indicates that the reference station other than the master reference station is in an abnormal state when the distance between the position based on the predetermined position information of the reference station other than the master reference station and the second positioning position is larger than the first threshold. Is determined (step S506).
  • the method for detecting the state of the reference station in the differential GPS receiving the GPS signal from the GPS satellite and detecting the occurrence of the abnormal state of the reference station that corrects the GPS signal based on its own position information. Can do.
  • the state detection method can accurately determine whether the reference station is abnormal or normal by setting a threshold value in consideration of the number and arrangement of GPS satellites.
  • a reference station in which the distance between the position information (surveying position) given in advance of the reference station and the first positioning position of the reference station is included in a predetermined range is set as a master reference station. is there. That is, the third embodiment of the present invention selects the master reference station using the absolute value of the distance between the surveying position of the reference station and the first positioning position of the reference station.
  • the configuration of the differential GPS 101 in the third embodiment of the present invention is the same as the configuration example of the differential GPS 101 in the first embodiment shown in FIG.
  • the configuration example of the correction value processing apparatus 401 in the third embodiment of the present invention is the same as the configuration example of the correction value processing apparatus 401 in the first embodiment shown in FIG.
  • the configuration example of the correction value processing apparatus 401 in the third embodiment of the present invention is the same as the configuration example of the data processing unit 403 in the second embodiment shown in FIG.
  • the master reference station selection unit 407 selects a reference station in which the distance between the position information (surveying position) given in advance of the reference station and the first positioning position of the reference station is included in a predetermined range. Select as master reference station.
  • the three-dimensional coordinate value of the first positioning position of the reference station obtained by single positioning is described as follows.
  • Reference station 301 ( xsa1 , ysa1 , zsa1 ); Reference station 302: (x sa2 , y sa2 , z sa2 ); Reference station 303: (x sa3 , y sa3 , z sa3 ); Reference station 304: ( xsa4 , ysa4 , zsa4 ).
  • the three-dimensional coordinate values of the surveying positions with respect to the reference stations 301 to 304 are described as follows.
  • Reference station 301 (x 1 , y 1 , z 1 ); Reference station 302: (x 2 , y 2 , z 2 ); Reference station 303: (x 3 , y 3 , z 3 ); Reference station 304: (x 4 , y 4 , z 4 ).
  • the master reference station selection unit 407 obtains the distance between the first positioning position and the surveying position for each of the reference stations 301 to 304 as a three-dimensional distance from these coordinate values.
  • Reference station 303: ⁇ sa3 ⁇ ((x sa3 ⁇ x 3 ) 2 + (y sa3 ⁇ y 3 ) 2 + (z sa3 ⁇ z 3 ) 2 );
  • Reference station 304: ⁇ sa4 ⁇ (( x sa4 -x 4) 2 + (y sa4 -y 4) 2 + (z sa4 -z 4) 2).
  • the master reference station selection unit 407 determines whether or not each of the obtained distances ⁇ sa1 , ⁇ sa2 , ⁇ sa3, and ⁇ sa4 between the first positioning position and the surveying position is within a predetermined threshold (predetermined range). To do.
  • the master reference station selection unit 407 selects, as a master reference station , a reference station corresponding to a distance that is within a predetermined threshold value among ⁇ sa1 , ⁇ sa2 , ⁇ sa3, and ⁇ sa4 . For example, when ⁇ sa1 that is the distance of the reference station 301 and ⁇ sa2 that is the distance of the reference station 302 are within a predetermined range, the master reference station selection unit 407 selects the reference station 301 and the reference station 302 as the master reference station.
  • the master reference station selection unit 407 the result of the determination, delta sa1, delta sa2, is within a predetermined threshold value (predetermined range) of the delta sa3 and delta sa4 if the "distance" is not, increases the predetermined threshold value It is possible to make the determination again (expand the predetermined range).
  • the master reference station selection unit 407 may increase the predetermined threshold every time a determination is made (expand the predetermined range), and repeatedly execute the determination until the master reference station can be selected.
  • the master reference station selection unit 407 the result of the determination, delta sa1, delta sa2, if delta sa3 and delta predetermined threshold (predetermined range) of the sa4 within the "distance" is not, as in the first embodiment
  • the master reference station may be selected based on the magnitude relationship of the distance between the first positioning position and the surveying position.
  • the master reference station selection unit 407 the result of the determination, delta sa1, delta sa2, when the distance within the predetermined threshold (predetermined range) of the delta sa3 and delta sa4 absence may not select a master reference station. In this case, the correction value processing apparatus 401 may terminate the process without executing the subsequent process.
  • the master reference station selection unit 407 the result of the determination, delta sa1, delta sa2, if all delta sa3 and delta sa4 is within the predetermined threshold (predetermined range), by decreasing the predetermined threshold (predetermined range). The determination may be performed again. In this case, each time the determination is made, the master reference station selection unit 407 decreases the predetermined threshold value (narrows the predetermined range), and repeats the determination until there is a reference station other than the master reference station. Good.
  • the master reference station selection unit 407 the result of the determination, delta sa1, delta sa2, if all delta sa3 and delta sa4 is within the predetermined threshold (predetermined range), as in the first embodiment, the first The master reference station may be selected based on the magnitude relationship between the distance between the positioning position and the surveying position.
  • the master reference station selection unit 407 the result of the determination, delta sa1, delta sa2, if all delta sa3 and delta sa4 is within the predetermined threshold (predetermined range), to the free selection of all reference stations as a master reference station Good.
  • the correction value processing apparatus 401 does not execute processing for determining whether or not the reference station other than the master reference station is in an abnormal state.
  • the predetermined range that is a condition for selecting as a master reference station may be changed according to a request of a user (such as an administrator), for example.
  • a user such as an administrator
  • the predetermined range is set wide, even a reference station having a certain distance between the first positioning position and the surveying position is selected as a master reference station. In this case, the correction value of the master reference station increases, and as a result, the correction amount of the first positioning position of the reference station other than the master reference station also increases.
  • the second positioning position of the reference station other than the master reference station is a value obtained by correcting the first positioning position of the reference station other than the master reference station with the correction value of the master reference station. That is, if the correction value of the master reference station is large, the second positioning position of the reference station other than the master reference station is a value obtained by largely correcting the first positioning position.
  • the error correction amount in the DGPS positioning is increased, and therefore the first positioning position approaches the surveying position. Therefore, the possibility that the distance between the surveying position of a reference station other than the master reference station and the second positioning position is greater than the first threshold value is reduced.
  • the abnormality determination unit 410 is performed based on, for example, the pass / fail determination logic of Table 1 in the first embodiment. In this case, if the abnormality determination unit 410 passes the determination with the correction value in one master reference station, the abnormality determination unit 410 determines that the pass is a comprehensive determination. Therefore, when a reference station having a certain distance between the first positioning position and the surveying position is selected as a master reference station, the correction value increases, and the abnormality determination unit 410 increases the number of cases where a reference station other than the master reference station is determined to be acceptable. . As a result, the number of reference stations other than the master reference station that the abnormality determination unit 410 determines to be in an abnormal state is reduced. In other words, by widening the predetermined range, the abnormal state determination unit 410 determines that a reference station other than the master reference station is normal if there is some error.
  • the predetermined range is set to be narrow
  • the master reference station is not selected unless the distance between the first positioning position and the surveying position is small.
  • the number of reference stations other than the master reference station that the abnormality determination unit 410 determines to be in an abnormal state increases.
  • the abnormality determination unit 410 determines that a reference station other than the master reference station is abnormal even with a slight error.
  • the third embodiment of the present invention by changing a predetermined range that is a condition for selecting as a master reference station, it is possible to flexibly determine the accuracy of determining whether a reference station other than the master reference station is abnormal or normal. Can be adjusted.
  • FIG. 6 is a flowchart showing an operation example of the correction value processing apparatus 401 in the third embodiment of the present invention.
  • the data receiving unit 402 of the correction value processing device 401 receives observation data from the reference stations 301 to 304 that have received the GPS signals emitted from the GPS satellites 201 to 20n (step S601).
  • the data processing unit 403 of the correction value processing apparatus 401 measures the position of each of the reference stations 301 to 304 based on the observation data (single positioning), and obtains a first positioning distance (step S602).
  • the data processing unit 403 selects, from the reference stations 301 to 304, a reference station in which the distance between the position information (surveying position) given in advance of the reference station and the first positioning position of the reference station is included in the predetermined range as the master reference station. (Step S603) Next, the data processing unit 403 calculates a correction value for the master reference station based on the observation data obtained from the master reference station and the value of the surveying position (step S604).
  • the data processing unit 403 obtains the second positioning position of the reference station other than the master reference station by the reference station other than the master reference station correcting the first positioning position of the reference station other than the master reference station with the correction value of the master reference station ( Step S605).
  • the data processing unit 403 indicates that the reference station other than the master reference station is in an abnormal state when the distance between the position based on the predetermined position information of the reference station other than the master reference station and the second positioning position is larger than the first threshold. Is determined (step S606).
  • the reference station in which the distance between the position information (surveying position) given in advance of the reference station and the first positioning position of the reference station is included in the predetermined range.
  • the state detection method can adjust the conditions for selecting the master reference station by setting the predetermined range in consideration of the number and arrangement of reference stations. As a result, the state detection method can flexibly change the accuracy of determining whether a reference station other than the master reference station is abnormal or normal.
  • a state detection method for detecting an abnormal state of a reference station in a positioning system having a satellite and a reference station Based on the signal received by the reference station from the satellite, the first positioning position of the reference station is obtained, Select a predetermined number of master reference stations in order from the smallest distance between the position information given in advance of the reference station and the first positioning position of the reference station, Based on the signal received by the master reference station from the satellite, a correction value for the master reference station is generated, By correcting the first positioning position of the reference station other than the master reference station with the correction value of the master reference station, the reference station other than the master reference station obtains the second positioning position of the reference station other than the master reference station, A reference station other than the master reference station is determined to be in an abnormal state when a distance between a position based on position information given in advance by a reference station other than the master reference station and the second positioning position is greater than a first threshold.
  • a state detection method for detecting an abnormal state of the plurality of reference stations, Based on the signal received by the reference station from the satellite, the first positioning position of the reference station is obtained, Selecting a reference station in which a distance between a predetermined position information of the reference station and the first positioning position of the reference station is included in a predetermined range as a master reference station; Based on the signal received by the master reference station from the satellite, a correction value for the master reference station is generated, By correcting the first positioning position of the reference station other than the master reference station with the correction value of the master reference station, the reference station other than the master reference station obtains the second positioning position of the reference station other than the master reference station, A reference station other than the master reference station is determined to be in an abnormal state when a distance between a position based on position information given in advance by a reference station other than the master reference station and the second positioning position is greater than a first threshold.
  • the state detection method characterized
  • Appendix 3 The state detection method according to appendix 1 or 2, wherein the correction value of the master reference station includes an ionospheric delay amount and a tropospheric delay amount with respect to the master reference station.
  • the correction value of the master reference station is a value obtained by subtracting the first positioning position obtained by performing a smoothing process on the pseudo distance between the master reference station and the satellite from a predetermined position of the master reference station.
  • the positioning system includes a plurality of the master reference stations, Any one of Supplementary notes 1 to 4, wherein the correction value used for correcting the pseudo distance of a reference station other than the master reference station is an average value of the correction values generated by each of the plurality of master reference stations.
  • the positioning system has a first master reference station and a second master reference station, It was obtained by correcting the pseudo distance between the first master reference station and the satellite obtained by the first master reference station based on the signal received from the satellite by the correction value of the second master reference station. Determining a second positioning position of the first master reference station; When the distance between the position of the first master reference station based on the position information given in advance and the second positioning position of the first master reference station is larger than a second threshold, the first master reference station It is determined that there is an abnormal condition, Supplementary notes 1 to 5, wherein when the first master reference station is determined to be in an abnormal state, one of the plurality of master reference stations is selected as the first master reference station.
  • the state detection method in any one of.
  • Appendix 7 The state detection method according to any one of appendices 1 to 6, wherein a satellite that transmits a signal to be received is selected from the plurality of satellites.
  • the satellite that transmits the signal to be received is selected from among a plurality of the satellites, a satellite having an elevation angle larger than a predetermined angle, or a satellite that is visible from all reference stations included in the positioning system.
  • Appendix 9 The state detection method according to any one of appendices 1 to 8, wherein the first threshold value is determined with reference to a predetermined probability distribution function.
  • Appendix 10 The state detection method according to any one of appendices 1 to 9, wherein the first threshold value is determined based on a number of the satellites and an arrangement state of the satellites.
  • Data receiving means for receiving a signal received from a satellite by a reference station; First positioning position calculating means for calculating a first positioning position of the reference station based on the signal received by the data receiving means; Master reference station selection means for selecting a predetermined number of master reference stations in order from the smallest distance between the position information given in advance of the reference station and the first positioning position of the reference station; A master reference station correction value calculating means for generating a correction value of the master reference station based on a signal received from the satellite by the master reference station; The reference station other than the master reference station calculates a second positioning position of a reference station other than the master reference station by correcting the first positioning position of the reference station other than the master reference station with the correction value of the master reference station.
  • a reference station other than the master reference station is determined to be in an abnormal state when a distance between a position based on position information given in advance by a reference station other than the master reference station and the second positioning position is greater than a first threshold.
  • a correction value processing apparatus comprising:
  • [Appendix 12] Data receiving means for receiving a signal received from a satellite by a reference station; First positioning position calculating means for calculating a first positioning position of the reference station based on the signal received by the data receiving means; Master reference station selection means for selecting, as a master reference station, a reference station whose distance between the position information given in advance of the reference station and the first positioning position of the reference station is included in a predetermined range; A master reference station correction value calculating means for generating a correction value of the master reference station based on a signal received from the satellite by the master reference station; The reference station other than the master reference station calculates a second positioning position of a reference station other than the master reference station by correcting the first positioning position of the reference station other than the master reference station with the correction value of the master reference station.
  • a reference station other than the master reference station is determined to be in an abnormal state when a distance between a position based on position information given in advance by a reference station other than the master reference station and the second positioning position is greater than a first threshold.
  • a correction value processing apparatus comprising:
  • Appendix 13 13 The correction value processing apparatus according to appendix 11 or 12, wherein the correction value of the master reference station includes an ionospheric delay amount and a tropospheric delay amount with respect to the master reference station.
  • Additional correction means transmitting means for transmitting the correction value generated by the data processing means to a receiving device that estimates its own position based on a signal received from the satellite,
  • the correction value processing apparatus according to any one of the above.
  • a positioning system having a satellite, a reference station, and a correction value processing device includes: First positioning position calculating means for calculating a first positioning position of the reference station based on the signal received by the data receiving means; Master reference station selection means for selecting a predetermined number of master reference stations in order from the smallest distance between the position information given in advance of the reference station and the first positioning position of the reference station; A master reference station correction value calculating means for generating a correction value of the master reference station based on a signal received from the satellite by the master reference station; The reference station other than the master reference station calculates a second positioning position of a reference station other than the master reference station by correcting the first positioning position of the reference station other than the master reference station with the correction value of the master reference station.
  • a reference station other than the master reference station is determined to be in an abnormal state when a distance between a position based on position information given in advance by a reference station other than the master reference station and the second positioning position is greater than a first threshold.
  • a positioning system having a satellite, a reference station, and a correction value processing device includes: First positioning position calculating means for calculating a first positioning position of the reference station based on the signal received by the data receiving means; Master reference station selection means for selecting, as a master reference station, a reference station whose distance between the position information given in advance of the reference station and the first positioning position of the reference station is included in a predetermined range; A master reference station correction value calculating means for generating a correction value of the master reference station based on a signal received from the satellite by the master reference station; The reference station other than the master reference station calculates a second positioning position of a reference station other than the master reference station by correcting the first positioning position of the reference station other than the master reference station with the correction value of the master reference station.
  • a reference station other than the master reference station is determined to be in an abnormal state when a distance between a position based on position information given in advance by a reference station other than the master reference station and the second positioning position is greater than a first threshold.
  • Appendix 17 The positioning system according to appendix 15 or 16, wherein the correction value of the master reference station includes an ionospheric delay amount and a tropospheric delay amount with respect to the master reference station.
  • Appendix 18 Further comprising a receiving device for estimating its position based on a signal received from the satellite; The positioning system according to any one of appendices 15 to 17, wherein the correction value processing device further includes correction value transmitting means for transmitting the correction value generated by the data processing unit to the receiving device. .
  • a program for detecting an abnormal state of the reference station A process of calculating a first positioning position of the reference station based on the signal received by the data receiving means; A process of selecting a predetermined number of master reference stations in order from the smallest distance between the position information given in advance of the reference station and the first positioning position of the reference station; A process for generating a correction value for the master reference station based on a signal received by the master reference station from the satellite; The reference station other than the master reference station calculates a second positioning position of a reference station other than the master reference station by correcting the first positioning position of the reference station other than the master reference station with the correction value of the master reference station.
  • a reference station other than the master reference station is determined to be in an abnormal state when a distance between a position based on position information given in advance by a reference station other than the master reference station and the second positioning position is greater than a first threshold.
  • Appendix 20 In a positioning system having a satellite and a reference station, a program for detecting an abnormal state of the reference station, A process of calculating a first positioning position of the reference station based on the signal received by the data receiving means; A process of selecting, as a master reference station, a reference station in which a distance between a predetermined position information of the reference station and a distance between the first positioning position of the reference station is included in a predetermined range; A process for generating a correction value for the master reference station based on a signal received by the master reference station from the satellite; The reference station other than the master reference station calculates a second positioning position of a reference station other than the master reference station by correcting the first positioning position of the reference station other than the master reference station with the correction value of the master reference station.
  • a reference station other than the master reference station is determined to be in an abnormal state when a distance between a position based on position information given in advance by a reference station other than the master reference station and the second positioning position is greater than a first threshold.
  • Appendix 21 The storage medium storing the program according to appendix 19 or 20, wherein the correction value of the master reference station includes at least an ionospheric delay amount and a tropospheric delay amount with respect to the master reference station.
  • a state detection method for detecting an abnormal state of the reference station Obtaining a first positioning position of the reference station based on a pseudo distance between the reference station and the satellite obtained based on a signal received from the satellite by the reference station; Selecting a predetermined number of master reference stations in order from the smallest distance between the position of the reference station according to the position information given in advance and the first positioning position of the reference station; Based on the signal received by the master reference station from the satellite, a correction value for the master reference station is generated, By correcting the pseudo distance of a reference station other than the master reference station with the correction value of the master reference station, a second positioning position of a reference station other than the master reference station is obtained, A reference station other than the master reference station when a distance between a position based on position information given in advance by a reference station other than the master reference station and the second positioning position of a reference station other than the master reference station is greater than a first threshold value It
  • Appendix 23 The state detection method according to appendix 22, wherein the correction value of the master reference station includes an ionospheric delay amount and a tropospheric delay amount with respect to the master reference station.
  • the correction value of the master reference station is a value obtained by subtracting the first positioning position obtained by performing a smoothing process on the pseudo distance between the master reference station and the satellite from a predetermined position of the master reference station. 24.
  • the positioning system includes a plurality of the master reference stations, Any one of appendices 22 to 24, wherein the correction value used for correcting the pseudo distance of a reference station other than the master reference station is an average value of the correction values generated by each of the plurality of master reference stations.
  • the positioning system has a first master reference station and a second master reference station, It was obtained by correcting the pseudo distance between the first master reference station and the satellite obtained by the first master reference station based on the signal received from the satellite by the correction value of the second master reference station. Determining a second positioning position of the first master reference station; When the distance between the position of the first master reference station based on position information given in advance and the second positioning position of the first master reference station is larger than a second threshold, the first master reference station Is determined to be abnormal, Supplementary notes 22 to 25, wherein when determining that the first master reference station is in an abnormal state, one of the plurality of master reference stations is selected as the first master reference station.
  • the state detection method in any one of.
  • Appendix 27 The state detection method according to any one of appendices 22 to 26, wherein when obtaining the first positioning position of the reference station, a satellite that transmits a signal to be received is selected from the plurality of satellites.
  • the satellite that transmits the signal to be received is selected from among a plurality of the satellites, a satellite having an elevation angle larger than a predetermined angle, or a satellite that is visible from all reference stations included in the positioning system.
  • Appendix 30 The state detection method according to any one of appendices 22 to 29, wherein the first threshold value is determined based on the number of the satellites and an arrangement state of the satellites.
  • Data receiving means for receiving a signal received from a satellite by a reference station A first positioning position of the reference station is obtained based on a pseudo-range between the reference station and the satellite obtained based on the signal received by the data receiving means, and a position based on position information given in advance of the reference station; A predetermined number of master reference stations are selected in order from the smallest distance from the first positioning position of the reference station, and the correction value of the master reference station is determined based on the signal received by the master reference station from the satellite.
  • the second positioning position of the reference station other than the master reference station is obtained by correcting the pseudo distance of the reference station other than the master reference station by the correction value of the master reference station.
  • Appendix 32 32.
  • a positioning system having a satellite, a reference station, and a correction value processing device includes: The reference station obtains a first positioning position of the reference station based on a pseudo distance between the reference station and the satellite, which is obtained based on a signal received from the satellite, and a position based on position information given in advance of the reference station, and A predetermined number of master reference stations are selected in ascending order of distance from the first positioning position of the reference station, and a correction value for the master reference station is generated based on a signal received from the satellite by the master reference station.
  • the second positioning position of the reference station other than the master reference station is obtained by correcting the pseudo distance of the reference station other than the master reference station by the correction value of the master reference station, and given in advance by the reference station other than the master reference station.
  • Reference station other than the terpolymer reference station is determined to be abnormal state, and having a data processing unit, the positioning system.
  • Appendix 36 Further comprising a receiving device for estimating its position based on a signal received from the satellite; 36.
  • a program for detecting an abnormal state of the reference station Processing for obtaining a first positioning position of the reference station based on a pseudo distance between the reference station and the satellite obtained based on a signal received from the satellite by the reference station; A process of selecting a predetermined number of master reference stations in order from the smallest distance between the position of the reference station based on the position information given in advance and the first positioning position of the reference station; Based on the signal received from the satellite by the master reference station, processing to generate a correction value for the master reference station; A process of obtaining a second positioning position of a reference station other than the master reference station by correcting the pseudo distance of a reference station other than the master reference station with the correction value of the master reference station; A reference station other than the master reference station when a distance between a position based on position information given in advance by a reference station other than the master reference station and the second positioning position of a reference station other than the master reference station is greater than
  • Appendix 38 The storage medium storing the program according to appendix 37, wherein the correction value of the master reference station includes at least an ionospheric delay amount and a tropospheric delay amount with respect to the master reference station.
  • the present invention is not limited to the above embodiment, and can be suitably applied to a positioning system that uses a plurality of reference stations.

Abstract

[課題]衛星と基準局を有する測位システムにおいて、基準局の異常状態を高い精度で検出する状態検出方法、測位システム、及びプログラムを提供すること。 [解決手段] 衛星と基準局とを有する測位システムにおいて基準局の異常状態を検出する状態検出方法であって、前記基準局が前記衛星から受信した信号に基づき、前記基準局の第1測位位置を求め、前記基準局の予め与えられた位置情報と、前記基準局の前記第1測位位置との間の距離が小さいものから順に所定数のマスター基準局を選択し、前記マスター基準局が前記衛星から受信した信号に基づき、前記マスター基準局の補正値を生成し、前記マスター基準局以外の前記基準局が前記マスター基準局以外の前記基準局の前記第1測位位置を前記マスター基準局の前記補正値により補正することにより、前記マスター基準局以外の基準局の第2測位位置を求め、前記マスター基準局以外の基準局の予め与えられた位置情報による位置と、前記第2測位位置との距離が第1の閾値よりも大きいときに、前記マスター基準局以外の基準局が異常状態であると判定することを特徴とする。

Description

状態検出方法、補正値処理装置、測位システム、および状態検出プログラム
 本発明は、状態検出方法、補正値処理装置、測位システム、および状態検出プログラムに関し、特に、複数の基準局を有するディファレンシャルGPSにおいて、基準局の異常状態を検出する方法に関する。
 近年、GPS(Global Positioning System)を利用して測位を行う場面がますます増加しており、日常生活のGPSへの依存が高まっている。
 ここで、低仰角からのGPS衛星からの信号を受信する場合など、信号を搬送する電波が伝播する電離層及び大気圏の影響により、伝播に遅延が生じることがあり、GPSを利用した測位の精度には限界があった。
 この問題を解決するため、地表に固定された基準局を設置することにより、GPSによる測位の精度を向上させた、ディファレンシャルGPSが普及している。
 しかしながら、ディファレンシャルGPSは、基準局の時刻と位置の精度に、測位精度が依存する。例えば、基準局の動作に異常が生じると、補正信号に影響し、補正信号を受信してGPS衛星から受信した信号を補正して測位を行うGPSユーザ受信機は、正しい位置とは異なる位置を自身の位置と推定することになる。
 GPSによる測位の誤差は、日常生活への影響が甚大である。このため、ディファレンシャルGPSにおいては、基準局の動作の障害を早期に検出し、補正値への影響を回避することが求められている。
 現在運用されているディファレンシャルGPSの基準局に対して、GBAS(地上型補強システム:Ground Based Augmentation System)標準として定められている要件は、インテグリティリスクが150秒あたり10-5である。なお、インテグリティリスクは、完全性喪失確率のことである。この要件は、ICAO(国際民間航空機関:International Civil Aviation Organization)/RTCA(航空無線技術委員会:Radio Technical Comission for Aeronautics)が定めている。この要件を満たす基準局について、同時に2つの基準局が故障する確率は、10-10程度になる。GBASの運用において、要件としてのインテグリティリスクが10-7であるカテゴリI(CAT-I)ではこの数値は無視してよいが、インテグリティリスクが10-9であるカテゴリII及びIII(CAT-II/III)では無視できない大きさである。
 特許文献1は、複数の基準局を使用するディファレンシャルGPS(DGPS)について、異常が発生した基準局を検出する方法を開示する。N個のGPS衛星とK個の基準局にそれぞれインデックスnとkを付して特定するときに、クロックバイアスを除去し、共通測定時間を反映するように調整された受信機/衛星固有差動訂正値をCn,kとする。この受信機/衛星固有作動訂正値は、GPS衛星から受信したGPS信号をもとに、差動訂正処理手段により生成される。すなわち、まず、受信したGPS信号と、地上に固定された基準局のGPS信号受信機の既知の測量位置とを基に、衛星固有差動訂正値を生成する。GPS信号受信機は、GPS衛星の原子クロック時間に対してクロック時間オフセットまたはバイアスを有しており、受信機/衛星固有測定値はそれぞれ異なる測定時間を有する。GPS信号受信機は、受信機/衛星固有測定値より、クロックバイアスを伴わない差動訂正値を生成する。次に、GPS信号受信機は、共通同期時間に応じて受信機/衛星固有訂正値を調整する。これらの処理により、クロックバイアスを除去し、共通測定時間を反映するように調整された受信機/衛星固有差動訂正値Cn,kが得られる。
 次に、衛星固有差動訂正平均値が、上記のクロックバイアスを除去し共通測定時間を反映するように調整された受信機/衛星固有差動訂正値Cn,kを基準局の間で相加平均して得られた値<C>であるとする。すなわち、衛星固有差動訂正平均値は、次の関係式:
 <C>=(1/K)(Cn,1+Cn,2+・・・+Cn,K )
に従って求められる。ここで、検定に用いる量として、受信機/衛星固有弁別値Zn,kを用いる。この受信機/衛星固有弁別値Zn,kは、受信機/衛星固有差動訂正値Cn,kと、衛星固有差動訂正平均値<C>との差分であるとする。すなわち、次の関係式:
 Zn,k =Cn,k-<C
により求められる。このとき、予め定めた検出閾値DTn,kにより、基準局の判定を行う。すなわち、判定式:
 |Zn,k|>DTn,k
により基準局を判定する。具体的には、与えられた基準局kに対して、全てのGPS衛星(n=1,..,N)で判定式:|Zn,k|>DTn,kが成立する場合に、基準局kに異常があるものと判定する。
 また、特許文献2乃至5にも、関連する技術が開示されている。
米国特許5,600、329公報 特開2003-057327号公報 特開2007-010422号公報 特開2009-168804号公報 特開2011-043449号公報
 特許文献1は、異常の発生した基準局を検出することができるが、以下のような課題を有する。すなわち、上記の判定式による検定は、レンジドメイン(GPSユーザ受信機(航空機など)からGPS衛星までの距離測定)のみで行われる。従って、受信機/衛星固有弁別値の絶対値が閾値未満である場合、異常が発生しているとはみなされない。しかしながら、受信機/衛星固有弁別値の符号が負側或いは正側に偏っている場合は、ポジションドメイン(GPSユーザ受信機の位置の決定)では致命的な位置誤差を発生するおそれがある。また、基準局間平均値の算出には、異常が発生した基準局からの補正値も含まれるため、検定に用いた受信機/衛星固有弁別値が異常の有無の判定を正確に反映しているとはいえない。例えば、異常のある基準局に対しては、受信機/衛星固有弁別値が小さくなり、検出失敗が発生するおそれがある。また、異常のない基準局に対して受信機/衛星固有弁別値が大きくなり、誤検出が発生するおそれがある。
 上記の目的を達成するため、本発明の状態検出方法は、衛星と基準局とを有する測位システムにおいて基準局の異常状態を検出する状態検出方法であって、前記基準局が前記衛星から受信した信号に基づき、前記基準局の第1測位位置を求め、前記基準局の予め与えられた位置情報と、前記基準局の前記第1測位位置との間の距離が小さいものから順に所定数のマスター基準局を選択し、前記マスター基準局が前記衛星から受信した信号に基づき、前記マスター基準局の補正値を生成し、前記マスター基準局以外の前記基準局が前記マスター基準局以外の前記基準局の前記第1測位位置を前記マスター基準局の前記補正値により補正することにより、前記マスター基準局以外の基準局の第2測位位置を求め、前記マスター基準局以外の基準局の予め与えられた位置情報による位置と、前記第2測位位置との距離が第1の閾値よりも大きいときに、前記マスター基準局以外の基準局が異常状態であると判定することを特徴とする。
 また、本発明の状態検出方法は、衛星と複数の基準局とを有する測位システムにおいて、前記複数の基準局の異常状態を検出する状態検出方法であって、前記基準局が前記衛星から受信した信号に基づき、前記基準局の第1測位位置を求め、前記基準局の予め与えられた位置情報と、前記基準局の前記第1測位位置との距離が所定範囲に含まれている基準局をマスター基準局として選択し、前記マスター基準局が前記衛星から受信した信号に基づき、前記マスター基準局の補正値を生成し、前記マスター基準局以外の前記基準局が前記マスター基準局以外の前記基準局の前記第1測位位置を前記マスター基準局の前記補正値により補正することにより、前記マスター基準局以外の基準局の第2測位位置を求め、前記マスター基準局以外の基準局の予め与えられた位置情報による位置と、前記第2測位位置との距離が第1の閾値よりも大きいときに、前記マスター基準局以外の基準局が異常状態であると判定することを特徴とする。
 また、本発明の補正値処理装置は、基準局が衛星から受信した信号を受信するデータ受信手段と、前記データ受信手段が受信した前記信号基づき、前記基準局の第1測位位置を算出する第1測位位置算出手段と、前記基準局の予め与えられた位置情報と、前記基準局の前記第1測位位置との間の距離が小さいものから順に所定数のマスター基準局を選択するマスター基準局選択手段と、前記マスター基準局が前記衛星から受信した信号に基づき、前記マスター基準局の補正値を生成するマスター基準局補正値算出手段と、前記マスター基準局以外の前記基準局が前記マスター基準局以外の前記基準局の前記第1測位位置を前記マスター基準局の前記補正値により補正することにより、前記マスター基準局以外の基準局の第2測位位置を算出する第2測位位置算出手段と、前記マスター基準局以外の基準局の予め与えられた位置情報による位置と、前記第2測位位置との距離が第1の閾値よりも大きいときに、前記マスター基準局以外の基準局が異常状態であると判定する異常判定手段と、を有することを特徴とする。
 また、本発明の補正値処理装置は、基準局が衛星から受信した信号を受信するデータ受信手段と、前記データ受信手段が受信した前記信号基づき、前記基準局の第1測位位置を算出する第1測位位置算出手段と、前記基準局の予め与えられた位置情報と、前記基準局の前記第1測位位置との距離が所定範囲に含まれている基準局をマスター基準局として選択するマスター基準局選択手段と、前記マスター基準局が前記衛星から受信した信号に基づき、前記マスター基準局の補正値を生成するマスター基準局補正値算出手段と、前記マスター基準局以外の前記基準局が前記マスター基準局以外の前記基準局の前記第1測位位置を前記マスター基準局の前記補正値により補正することにより、前記マスター基準局以外の基準局の第2測位位置を算出する第2測位位置算出手段と、前記マスター基準局以外の基準局の予め与えられた位置情報による位置と、前記第2測位位置との距離が第1の閾値よりも大きいときに、前記マスター基準局以外の基準局が異常状態であると判定する異常判定手段と、を有することを特徴とする。
 また、本発明の測位システムは、衛星と、基準局と、補正値処理装置を有する測位システムであって、前記補正値処理装置は、前記データ受信手段が受信した前記信号基づき、前記基準局の第1測位位置を算出する第1測位位置算出手段と、前記基準局の予め与えられた位置情報と、前記基準局の前記第1測位位置との間の距離が小さいものから順に所定数のマスター基準局を選択するマスター基準局選択手段と、前記マスター基準局が前記衛星から受信した信号に基づき、前記マスター基準局の補正値を生成するマスター基準局補正値算出手段と、前記マスター基準局以外の前記基準局が前記マスター基準局以外の前記基準局の前記第1測位位置を前記マスター基準局の前記補正値により補正することにより、前記マスター基準局以外の基準局の第2測位位置を算出する第2測位位置算出手段と、前記マスター基準局以外の基準局の予め与えられた位置情報による位置と、前記第2測位位置との距離が第1の閾値よりも大きいときに、前記マスター基準局以外の基準局が異常状態であると判定する異常判定手段と、を有することを特徴とする。
 また、本発明の測位システムは、衛星と、基準局と、補正値処理装置を有する測位システムであって、前記補正値処理装置は、前記データ受信手段が受信した前記信号基づき、前記基準局の第1測位位置を算出する第1測位位置算出手段と、前記基準局の予め与えられた位置情報と、前記基準局の前記第1測位位置との距離が所定範囲に含まれている基準局をマスター基準局として選択するマスター基準局選択手段と、前記マスター基準局が前記衛星から受信した信号に基づき、前記マスター基準局の補正値を生成するマスター基準局補正値算出手段と、前記マスター基準局以外の前記基準局が前記マスター基準局以外の前記基準局の前記第1測位位置を前記マスター基準局の前記補正値により補正することにより、前記マスター基準局以外の基準局の第2測位位置を算出する第2測位位置算出手段と、前記マスター基準局以外の基準局の予め与えられた位置情報による位置と、前記第2測位位置との距離が第1の閾値よりも大きいときに、前記マスター基準局以外の基準局が異常状態であると判定する異常判定手段と、を有することを特徴とする。
 また、本発明の記憶媒体は、衛星と基準局とを有する測位システムにおいて、前記基準局の異常状態を検出するプログラムを記憶した記憶媒体であって、前記データ受信手段が受信した前記信号基づき、前記基準局の第1測位位置を算出する処理と、前記基準局の予め与えられた位置情報と、前記基準局の前記第1測位位置との間の距離が小さいものから順に所定数のマスター基準局を選択する処理と、前記マスター基準局が前記衛星から受信した信号に基づき、前記マスター基準局の補正値を生成する処理と、前記マスター基準局以外の前記基準局が前記マスター基準局以外の前記基準局の前記第1測位位置を前記マスター基準局の前記補正値により補正することにより、前記マスター基準局以外の基準局の第2測位位置を算出する第2測位位置算出処理と、前記マスター基準局以外の基準局の予め与えられた位置情報による位置と、前記第2測位位置との距離が第1の閾値よりも大きいときに、前記マスター基準局以外の基準局が異常状態であると判定する処理と、とを含むことを特徴とする。
 また、本発明の記憶媒体は、衛星と基準局とを有する測位システムにおいて、前記基準局の異常状態を検出するプログラムを記憶した記憶媒体であって、前記データ受信手段が受信した前記信号基づき、前記基準局の第1測位位置を算出する処理と、前記基準局の予め与えられた位置情報と、前記基準局の前記第1測位位置との距離が所定範囲に含まれている基準局をマスター基準局として選択する処理と、前記マスター基準局が前記衛星から受信した信号に基づき、前記マスター基準局の補正値を生成する処理と、前記マスター基準局以外の前記基準局が前記マスター基準局以外の前記基準局の前記第1測位位置を前記マスター基準局の前記補正値により補正することにより、前記マスター基準局以外の基準局の第2測位位置を算出する第2測位位置算出処理と、前記マスター基準局以外の基準局の予め与えられた位置情報による位置と、前記第2測位位置との距離が第1の閾値よりも大きいときに、前記マスター基準局以外の基準局が異常状態であると判定する処理と、とを含むことを特徴とする。
 本発明によれば、衛星測位システムにおける基準局の異常状態を高い精度で検出することができる。
本発明の第1の実施形態に係るディファレンシャルGPSの構成の一例を示す。 本発明の第1の実施形態に係るディファレンシャルGPSの基準局の状態検出方法の手順の一例を示すフローチャートである。 本発明の第1の実施形態に係る補正値処理装置の構成の一例を示す。 本発明の第2の実施形態に係るデータ処理部の構成の一例を示す。 本発明の第2の実施形態に係るディファレンシャルGPSの基準局の状態検出方法の手順の一例を示すフローチャートである。 本発明の第3の実施形態に係るディファレンシャルGPSの基準局の状態検出方法の手順の一例を示すフローチャートである。
 発明を実施するための最良の形態について図面を参照して詳細に説明する。ただし、本発明は以下に示す実施形態に限定されない。
 [第1の実施形態]
 図1は、本発明の第1の実施形態に係るディファレンシャルGPS101の構成の一例を示す。
 受信装置601は、GPS衛星201乃至20nからGPS信号を受信して、自身の位置を推定する。受信装置601は、例えば航空機に搭載される。
 複数の基準局301乃至304のそれぞれは、GPS衛星201乃至20nからのGPS信号を受信し、観測データを補正値処理装置401に送信する。補正値処理装置401は、GPS衛星201乃至20nからのGPS信号を受信して測位する際の補正値を算出し、放送送信装置501に送る。
 放送送信装置は、例えば、VDB(VHF(Very High Frequency) Data Broadcast)送信装置である。
 受信装置601は、放送送信装置501から送信された補正値放送を受信し、GPS衛星201乃至20nから受信したGPS信号を補正して、自身の位置を推定する。
 図3は、本実施形態に係る補正値処理装置401の構成の一例を示す。補正値処理装置401は、基準局301乃至304から送信された観測データを受信するデータ受信部402と、補正値情報を受信装置601に提供するため放送送信装置501に送信する補正値送信部405を有する。補正値処理装置は、さらに、データ受信部402で受信した観測データを処理して補正値情報を生成するデータ処理部403と、基準局301乃至304についてのデータを保持するデータ保持部404を有する。データ保持部404は、さらに、データ処理部403で処理された情報を一時的に保持し、以下に説明する検定処理で必要とされるパラメータの値などを保持する。
 次に、本発明の第1の実施形態に係るディファレンシャルGPSの基準局の状態検出方法の手順を、図1乃至3を参照して説明する。
 なお、本実施形態では、基準局の数を4にしたが、これに限定されない。基準局の数は4より少なくてもよいし、4より多くてもよい。
 補正値処理装置401は、GPS衛星201乃至20nから発せられたGPS信号を受信した基準局301乃至304から観測データを受信する(図2ステップS201)。
 次に、補正値処理装置401は、観測データを基に、基準局301乃至304のそれぞれについての位置を測位(単独測位)する(ステップS202)。
 単独測位は、4つのGPS衛星からの擬似距離(GPS衛星から受信機までの距離)を利用して実行される。受信装置601に内蔵された時計は、GPS衛星の原子時計と比較して精度が数桁劣る。GPS衛星の原子時計が刻む電波の発信時刻は、正確でも、受信機で測定された到達時刻には大きな誤差が含まれる可能性がある。GPS衛星からの電波を受信する3次元位置とともに受信装置601の時計の誤差も未知量に加えた4つの未知量を決定するために、同時に4個以上のGPS衛星からの擬似距離が利用される。
 単独測位により基準局の測位位置が得られる。このとき、基準局の間で測位の精度にばらつきが出ないように、次の条件のうちの少なくとも一つを課してもよい。
 (1)低仰角(例えば10度以下)の衛星を測位計算から排除し、電離層及び大気圏の影響或いはマルチパスによる測位誤差を低減する。
 (2)全基準局から可視である共通衛星のみを測位計算に使用する。
 (3)与えられた衛星に対して、同一のクロック・エフェメリス(時刻・軌道)パラメータを用いて測位計算を行う。すなわち、衛星から送信されたGPS信号に対して、使用するクロック・エフェメリスパラメータを統一する。
 次に、補正値処理装置401は、基準局301乃至304から、マスター基準局と検定対象基準局を選定する(ステップS203)。
 それぞれの地表に固定された基準局の位置について、3次元の座標値が高精度の数値として予め与えられている。以下これを測量位置とよぶ。この測量位置についての座標値が、補正値処理装置401のデータ保持部404に保持される。
 なお、測位位置と測量位置との差分を、以下では測位誤差とよぶ。測位誤差は3次元の距離として定義される。これらの演算は、補正値処理装置401で行われる。
 各基準局に対して、単独測位により得られた測位位置を測量位置と比較し、測位誤差が最小である基準局と次に小さい基準局を選択する。測位誤差が最小の基準局をマスター基準局Aとし、次に小さい基準局をマスター基準局Bとする。
 なお、本実施形態では、2つのマスター基準局を選定したが、マスター基準局の数が多ければ、マスター基準局を3以上にしてもよい。
 マスター基準局に選定されなかった基準局は、検定対象基準局となる。本実施形態では、2つの基準局が残るので、それぞれ検定対象基準局X及び検定対象基準局Yとなる。
 次に、補正値処理装置401は、マスター基準局から得た観測データと測量位置の値を基に、マスター基準局でのDGPS補正値を算出する(ステップS204)。マスター基準局AでのDGPS補正値がDGPS補正値Aとなり、マスター基準局BでのDGPS補正値がDGPS補正値Bとなる。
 DGPS補正値は次のようにして求められる。
 基準局が測定する擬似距離は、GPS衛星と基準局との真の距離による誤差を含む。また、当該疑似距離は、基準局のクロックバイアス(基準局の表示する時刻と真の時刻との差)、又は、GPS衛星のクロックバイアスなどの内部ノイズによる誤差を含む。さらに、当該疑似距離は、基準局の経験する対流圏遅延量、基準局の経験する電離層遅延量、基準局の経験するマルチパス、および、基準局の経験する熱雑音などの内部ノイズのそれぞれによる誤差も含む。
 擬似距離補正値は、GPS衛星が放送しているGPS衛星の位置情報と基準局を測量して求めた位置から算出した、GPS衛星と基準局との幾何学的距離から、スムージング処理を行った基準局の擬似距離を差し引いた値として定義される。このGPS衛星と基準局との幾何学的距離は、GPS衛星と基準局との真の距離に、GPS衛星が放送しているGPS衛星の位置情報の誤差と、基準局を測量して得られた距離の誤差を加えたものである。
 上記の関係式から導出された擬似距離補正値は、衛星のクロックバイアスと、衛星の経験する対流圏遅延量と、衛星の経験する電離層遅延量と、受信装置の誤差として認識され分離可能な誤差と、その他の誤差を含む。
 得られた基準局による擬似距離補正値は、受信装置の擬似距離の補正に適用される。受信装置の擬似距離は、GPS衛星と受信装置との真の距離と、受信装置のクロックバイアスとして測位計算の過程で分離可能な誤差成分と、それぞれの観測量に固有或いは仮定された標準偏差を使用して評価された誤差成分を含む。
 実際の擬似距離補正値は、複数の基準局で平均する。これにより、マルチパスや受信装置のノイズなどのランダムな誤差が抑圧される。一方、電離層遅延量及び対流圏遅延量は平均化して残る。
 次に、マスター基準局の選定の検定が実行される(ステップS205)。2つのマスター基準局の間で、お互いのDGPS補正値を適用して、DGPS測位を行い、測位の結果がそれぞれの測量位置と比較される。マスター基準局の選択は、上記のように単独測位に依存するため、異常があるにも関わらず、偶然に測位誤差が小さくなる可能性がある。上記のマスター基準局の選定の検定は、この可能性を排除する。
 基準局が正常に動作していれば、単独測位の結果により得られる測位位置は、基準局の測量位置を精度よく再現する。すなわち、測位誤差が大きい基準局は、動作の障害が疑われる。障害が発生している基準局から生成された補正値は、GPS信号を正しい位置から逸脱する方向に補正するおそれがあり、受信装置601に提供する補正値情報に含めると、誤差を増大させる要因となる。
 基準局に異常が発生し、該基準局によるDGPS補正値が正常な測位結果を与えないにもかかわらず、単独測位において偶然に測位誤差が小さかったために、該基準局をマスター基準局とする選択は適切な選択ではない。このように、マスター基準局の選択が適切でなければ(ステップS206NO)、補正値処理装置401は、マスター基準局の選択を修正し(ステップS207)、新たに選択されたマスター基準局の補正値の算出を行う(ステップS204)。
 マスター基準局の選択が適切であれば(ステップS206YES)、補正値処理装置401は、マスター基準局のDGPS補正値を、検定対象基準局へ適用して、DGPS測位を行う(ステップS208)。
 次に、検定対象基準局の検定が実行される(ステップS209)。
 検定対象基準局のDGPS測位結果とそれぞれの基準局の測量位置との差分が求められる。2つのマスター基準局のDGPS補正値に対して、この差分の値が所定の閾値を越えた場合に、当該検定対象基準局が異常であると判定する。この合否を判定するロジックを表1に示す。
 表1.合否判定ロジック
Figure JPOXMLDOC01-appb-I000001
 上記の判定に使用する閾値は、測位の精度により変更してもよい。例えば、当該閾値は、GPS衛星の配置状況を反映するPDOP(Position Dilution of Precision)と衛星数をもとに決定されてもよい。これは、衛星数の増減により、PDOPが変化し、測位精度が変化することがあり、固定した閾値では、上記の判定が不正確になるおそれがあるためである。
 DGPS測位結果と測量位置の比較においては、垂直方向成分と水平方向成分に分解すると、誤差の規模を正確に反映できないので,三次元での距離を基に評価が行われる。この三次元距離を基にした検定処理には、マクスウェル・ボルツマン分布を適用する。
 GBASの運用においては、1つのGPS衛星からのGPS信号について、複数の基準局のそれぞれからのDGPS補正値が、基準局について平均される。そして、得られた平均値が受信装置に与えられる。正常に動作する複数の基準局からのDGPS補正値を平均し受信装置に提供することにより、基準局ごとのばらつきが抑制され、受信装置は、より高い精度で自身の位置を測位することができる。
 補正値処理装置401のデータ処理部403は、上記の検定の結果、不合格とされた検定対象基準局からの観測データを除外して補正値情報を生成する。そして、補正値送信部405は、該補正値情報を放送送信装置501に送信する。
 次に、補正値処理装置401におけるデータの処理について、詳細に説明する。
 まず、補正値処理装置401は、基準局301乃至304に対して単独測位を行う。上述のように、測位に誤差を生じるおそれのあるGPS衛星からのGPS信号が除外される。単独測位には、GPS衛星201乃至20nのn個のGPS衛星からのGPS信号を使用する。単独測位によって得られた基準局の測位位置の三次元座標値は次のように記される。
 基準局301:(xsa1,ysa1,zsa1);
 基準局302:(xsa2,ysa2,zsa2);
 基準局303:(xsa3,ysa3,zsa3);
 基準局304:(xsa4,ysa4,zsa4)。
 また、基準局301乃至304に対する測量位置の三次元座標値は、次のように記される。
 基準局301:(x,y,z);
 基準局302:(x,y,z);
 基準局303:(x,y,z);
 基準局304:(x,y,z)。
 これらの座標値から、基準局301乃至304のそれぞれについての、単独測位位置と測量位置との差分は、三次元距離として求められる。
 基準局301:Δsa1=√((xsa1-x+(ysa1-y+(zsa1-z);
 基準局302:Δsa2=√((xsa2-x+(ysa2-y+(zsa2-z);
 基準局303:Δsa3=√((xsa3-x+(ysa3-y+(zsa3-z);
 基準局304:Δsa4=√((xsa4-x+(ysa4-y+(zsa4-z)。
 次に、上記の処理により得られた差分について、大小関係が調べられる。本実施形態では、例えば次のような大小関係が得られたとする。
 Δsa3>Δsa1>Δsa4>Δsa2
すなわち、基準局303についての単独測位位置と測量位置との差分が最も大きく、基準局302についての差分が最も小さくなったとする。
 この大小関係を基に、最も小さい差分を与える基準局と次に小さい差分を与える基準局がマスター基準局となる。すなわち、基準局302がマスター基準局Aとなり、基準局304をマスター基準局Bとなる。これらより大きな差分を与える基準局は、検定対象基準局とする。本実施例については、基準局303が検定対象基準局Xとなり、基準局301が検定対象基準局Yとなる。
 次に、マスター基準局でのDGPS補正値が求められる。
 単独測位に使用したn個のGPS衛星201乃至20nからのGPS信号をもとに、マスター基準局において求めたDGPS補正値が次のように記される。
 検定対象基準局XのCによる測位位置:(xdgpsXA,ydgpsXA,zdgpsXA);
 検定対象基準局XのCによる測位位置:(xdgpsXB,ydgpsXB,zdgpsXB);
 検定対象基準局YのCによる測位位置:(xdgpsYA,ydgpsYA,zdgpsYA);
 検定対象基準局YのCによる測位位置:(xdgpsYB,ydgpsYB,zdgpsYB)。
 得られた座標値は、下記で説明する基準局の検定に用いられる。
 同様にして、マスター基準局Aとマスター基準局Bとの間でそれぞれのDGPS補正値を適用して、お互いのDGPS測位が実行される。得られた測位位置は、次のように記される。
 マスター基準局AのCによる測位位置:(xdgpsAB,ydgpsAB,zdgpsAB);
 マスター基準局BのCによる測位位置:(xdgpsBA,ydgpsBA,zdgpsBA)。
 得られた測位位置の結果と測量位置との差分が三次元距離で求められる。すなわち、
 マスター基準局Aについて:
 ΔdgpsA=√((xdgpsAB-x+(ydgpsAB-y+(zdgpsAB-z);
 マスター基準局Bについて:
 ΔdgpsB=√((xdgpsBA-x+(ydgpsBA-y+(zdgpsBA-z
となる。
 次に、得られた結果を用いて、マスター基準局の検定が実行される。
 上記の処理により求められた、マスター基準局についての測位位置と測量位置との差分が、所定の閾値THmasterと比較される。すなわち、
 マスター基準局A:ΔdgpsA<THmasterであれば、マスター基準局Aは合格;
 マスター基準局B:ΔdgpsB<THmasterであれば、マスター基準局Bは合格。
 なお、閾値THmasterは、測位の精度を反映するように適宜定めてよいが、例えば、GPS衛星の配置状況を反映するPDOP(Position Dilution of Precision)とGPS衛星の個数を基に決定してもよい。
 ここで、rが、三次元誤差を対象とする系に典型的な長さで除して無次元化した量であるとする。この場合、下記の式(1)で示される、当該rについての確率分布関数であるマクスウェル・ボルツマン分布

を参照し、最大許容誤警報確率、最大許容検出失敗確率、フォルトフリー時の誤差の標準偏差などを基に、閾値の値THmasterが決定されてもよい。
 上記の検定により、マスター基準局が不合格になった場合、適切なマスター基準局を選定するために、不合格になったマスター基準局と検定対象基準局とが入れ替えられる。この操作は、マスター基準局の選定が、単独測位の結果に基づいており、単独測位と測量位置との差分が偶然に小さくなっている可能性があり、これを排除する処理である。
 いずれか一方のマスター基準局が不合格になった場合は、不合格になったマスター基準局と、検定対象基準局のうち、単独測位の結果と測量位置との差分が最も小さい検定対象基準局とが入れ替えられる。
 両方のマスター基準局が不合格になった場合は、両方のマスター基準局と2つの検定対象基準局とが入れ替えられる。本実施形態は、4つのマスター基準局を含むDGPSの構成を示すが、マスター基準局の数が4より大きい場合は、単独測位の結果と測量位置との差分が最も小さい検定対象基準局と次に小さい検定対象基準局が、2つのマスター基準局と入れ替えられる。
 マスター基準局を入れ替えたら、それぞれのマスター基準局でのDGPS補正値が算出され、上記の処理が繰り返される。
 なお、基準局に障害が発生し、観測データに誤りがある場合は、観測データを基に算出したDGPS補正値にも誤りが含まれる。このため、例えば、マスター基準局Aに障害が発生している場合は、マスター基準局AによるDGPS補正値であるCにも誤りが含まれる。すなわち、マスター基準局AのCによる測位位置には、マスター基準局Aからの観測データに含まれる誤りが反映され、マスター基準局BのCによる測位位置には、マスター基準局AによるDGPS補正値Cに含まれる誤りが反映される。
 従って、上記の検定によって、マスター基準局Aが不合格とされ、マスター基準局Aを検定対象基準局の一つと入れ替えて検定処理を行い、再びマスター基準局が不合格になった場合、マスター基準局Aが戻され、マスター基準局Bが入れ替えられてもよい。
 同様にして、両方のマスター基準局が不合格になった場合も、片方のマスター基準局のみに障害が発生している可能性があるので、マスター基準局の入れ替え後の検定処理の結果によっては、いずれか片方のマスター基準局が戻されてもよい。
 次に、検定対象基準局でのDGPS測位結果と測量位置との差分が三次元距離として算出される。
  検定対象基準局XのCによる測位位置と測量位置との差分:
 ΔdgpsXA=√((xdgpsXA-x+(ydgpsXA-y+(zdgpsXA-z);
 検定対象基準局XのCによる測位位置と測量位置との差分:
 ΔdgpsXB=√((xdgpsXB-x+(ydgpsXB-y+(zdgpsXB-z);
 検定対象基準局YのCによる測位位置と測量位置との差分:
 ΔdgpsYA=√((xdgpsYA-x+(ydgpsYA-y+(zdgpsYA-z);
 検定対象基準局XのCによる測位位置と測量位置との差分:
 ΔdgpsYB=√((xdgpsYB-x+(ydgpsYB-y+(zdgpsYB-z)。
 これらを用いて検定対象基準局の検定が実行される。
 すなわち、
 (ΔdgpsXA<THtest)かつ(ΔdgpsXB<THtest)であれば検定対象基準局Xは不合格、それ以外の場合は合格、
 (ΔdgpsYA<THtest)かつ(ΔdgpsYB<THtest)であれば検定対象基準局Yは不合格、それ以外の場合は合格
とする。
 閾値THtestは、THmasterと同様に、測位の精度を反映するように、適宜定めてよいが、例えば、GPS衛星の配置状況を反映するPDOP(Position Dilution of Precision)とGPS衛星の個数を基に決定してもよい。ここで、rが、三次元誤差を対象とする系に典型的な長さで除して無次元化した量とする。この場合に、rについての確率分布関数である、マクスウェル・ボルツマン分布を参照し、最大許容誤警報確率、最大許容検出失敗確率、フォルトフリー時の誤差の標準偏差などを基に、閾値の値THtestが決定されてもよい。
 上記の処理の結果、補正値処理装置401が放送送信装置501に送信する補正値情報は、検定に合格した基準局のDGPS補正値を含む。
 本実施形態に係るディファレンシャルGPSにおける基準局の状態検出方法によれば、GPS衛星からGPS信号を受信して、自身の位置情報をもとにGPS信号を補正する基準局の異常状態の発生が高い精度で検出できる。すなわち、マスター基準局自身の検定を行うことにより、参照される基準局が正常に動作する基準局であると信頼することができる。これにより、異常の発生のおそれのある基準局について、正常に動作すると判断される基準局からの観測データを基に、異常状態の検出ができる。このため、複数の基準局に異常が発生しても、異常状態が精度よく検出できる。また、正常に動作している基準局に異常が発生していると誤って判断するおそれもない。
 また、GPS衛星の数や配置を考慮して閾値を設定することにより、異常であるか正常であるかが精度よく判定できる。
 本発明の実施形態においては、ディファレンシャルGPSにおける補正値情報を使用して、基準局の異常状態を検出するが、これに限定されない。
 自身の時刻と位置の情報を定期的に発信する信号源からの信号を受信して、自装置の位置を推定する測位システムにおいて、該信号の補正情報の生成処理において参照される、固定された基準局の異常状態検出に、本発明は好適に適用される。
 すなわち、本実施形態では、受信装置が航空機に搭載される例を示したが、これに限定されず、受信装置は船舶に搭載されてもよいし、自動車や歩行者が保持してもよい。また、本実施形態では、補正値情報はVHF放送電波として受信装置に提供されるが、これに限定されず、この他の周波数帯の電波でもよいし、有線で伝送されてもよい。さらに、本発明は、水中での測位システムに適用することも可能である。
 なお、上述の処理動作を実行させるためのプログラムを、コンピュータ読取可能な記録媒体に格納して配布し、当該プログラムをコンピュータにインストールすることにより、上述の処理動作を実行する装置が構成されても良い。コンピュータに読取可能な記録媒体は、例えば、フレキシブルディスク、CD-ROM(Compact Disk Read-Only Memory)、MO(Magneto-Optical disk)などである。
 [第2の実施形態]
 本発明の第2の実施形態について、図面を参照して説明する。
 本発明の第2の実施形態におけるディファレンシャルGPS101の構成は、図1に示す第1の実施形態におけるディファレンシャルGPS101の構成の例と同様である。また、本発明の第2の実施形態における補正値処理装置401の構成例は、図3に示す第1の実施形態における補正値処理装置401の構成例と同様である。
 図4は、本発明の第2の実施形態における補正値処理装置401に含まれるデータ処理部403の構成例を示す図である。
 データ処理部403は、第1測位位置算出部406と、マスター基準局選定部407と、マスター基準局補正値算出部408と、第2測位位置算出部409と、異常判定部410とを備える。
 第1測位位置算出部406は、基準局が衛星から受信した信号に基づき、基準局の第1測位位置を求める。第1測位位置算出部406は、基準局301乃至304に対して単独測位を行う。第1測位位置算出部406は、単独測位として、第1の実施形態で説明したように、GPS衛星201乃至20nのn個のGPS衛星からのGPS信号を使用し、基準局301乃至304の測位位置の三次元座標値を求める。
 マスター基準局選定部407は、基準局の予め与えられた位置情報(測量位置)と、第1測位位置算出部406が求めた第1測位位置との間の距離が小さいものから順に所定数のマスター基準局を選択する。マスター基準局選定部407は、第1の実施形態で説明したように、基準局301乃至304それぞれについての測位位置と測量位置との差分を三次元距離として求め、当該差分について、大小関係を算出する。次に、マスター基準局選定部407は、算出した大小関係を基に、例えば最も小さい差分を与える基準局と次に小さい差分を与える基準局をマスター基準局として選択する。
 マスター基準局補正値算出部408は、マスター基準局が衛星から受信した信号に基づき、マスター基準局の補正値を生成する。マスター基準局の補正値は、例えば、第1の実施形態におけるDGPS補正値である。
 第2測位位置算出部409は、マスター基準局以外の基準局が、マスター基準局以外の基準局の第1測位位置を、マスター基準局の補正値により補正することにより、マスター基準局以外の基準局の第2測位位置を求める。第2測位位置算出部409は、マスター基準局以外の基準局の第1測位位置に、マスター基準局の各々の補正値を適用して、第2測位位置を求める。すなわち、第2測位位置算出部409は、マスター基準局の数だけ、第2測位位置を求める。
 異常判定部410は、マスター基準局以外の基準局の予め与えられた位置情報による位置と、第2測位位置との距離が第1の閾値よりも大きいときに、マスター基準局以外の基準局が異常状態であると判定する。異常判定部401が、マスター基準局以外の基準局が異常であるとの判定は、例えば、第1の実施形態における表1の合否判定ロジックに基づいて行われる。
 図5は、本発明の第2の実施形態における、補正値処理装置401の動作例を示すフローチャートである。
 補正値処理装置401のデータ受信部402は、GPS衛星201乃至20nから発せられたGPS信号を受信した基準局301乃至304から観測データを受信する(ステップS501)。
 次に、補正値処理装置401のデータ処理部403は、観測データを基に、基準局301乃至304のそれぞれについての位置を測位(単独測位)し、第1測位距離を求める(ステップS502)。
 次に、データ処理部403は、基準局の予め与えられた位置情報(測量位置)と、第1測位位置との間の距離が小さいものから順に所定数のマスター基準局を選択する(ステップS503)。
 データ処理部403は、マスター基準局から得た観測データと測量位置の値を基に、マスター基準局の補正値を算出する(ステップS504)。
 データ処理部403は、マスター基準局以外の基準局が、当該マスター基準局以外の基準局の第1測位位置をマスター基準局の補正値により補正することにより、当該マスター基準局以外の基準局の第2測位位置を求める(ステップS505)。
 データ処理部403は、マスター基準局以外の基準局の予め与えられた位置情報による位置と、第2測位位置との距離が第1の閾値よりも大きいときに、マスター基準局以外の基準局が異常状態であると判定する(ステップS506)。
 本実施形態に係るディファレンシャルGPSにおける基準局の状態検出方法によれば、GPS衛星からGPS信号を受信して、自身の位置情報をもとにGPS信号を補正する基準局の異常状態の発生を検出することができる。
 また、当該状態検出方法は、GPS衛星の数や配置を考慮して閾値を設定することにより、基準局が異常であるか正常であるかを精度よく判定できる。
 [第3の実施形態]
 本発明の第3の実施形態について、図面を参照して説明する。
 本発明の第3の実施形態は、基準局の予め与えられた位置情報(測量位置)と、基準局の第1測位位置との距離が所定範囲に含まれている基準局を、マスター基準局とするものである。すなわち、本発明の第3の実施形態は、基準局の測量位置と基準局の第1測位位置との距離の絶対値を使用して、マスター基準局を選定する。
 本発明の第3の実施形態におけるディファレンシャルGPS101の構成は、図1に示す第1の実施形態におけるディファレンシャルGPS101の構成の例と同様である。また、本発明の第3の実施形態における補正値処理装置401の構成例は、図3に示す第1の実施形態における補正値処理装置401の構成例と同様である。さらに、本発明の第3の実施形態における補正値処理装置401の構成例は、図4に示す第2の実施形態におけるデータ処理部403の構成例と同様である。
 本発明の第3の実施形態において、マスター基準局選定部407は、基準局の予め与えられた位置情報(測量位置)と、基準局の第1測位位置との距離が所定範囲に含まれている基準局をマスター基準局として選択する。
 単独測位によって得られた基準局の第1測位位置の三次元座標値は次のように記される。
 基準局301:(xsa1,ysa1,zsa1);
 基準局302:(xsa2,ysa2,zsa2);
 基準局303:(xsa3,ysa3,zsa3);
 基準局304:(xsa4,ysa4,zsa4)。
 また、基準局301乃至304に対する測量位置の三次元座標値は、次のように記される。
 基準局301:(x,y,z);
 基準局302:(x,y,z);
 基準局303:(x,y,z);
 基準局304:(x,y,z)。
 マスター基準局選定部407は、これらの座標値から、基準局301乃至304のそれぞれについての、第1測位位置と測量位置との距離を、三次元距離として求める。
 基準局301:Δsa1=√((xsa1-x+(ysa1-y+(zsa1-z);
 基準局302:Δsa2=√((xsa2-x+(ysa2-y+(zsa2-z);
 基準局303:Δsa3=√((xsa3-x+(ysa3-y+(zsa3-z);
 基準局304:Δsa4=√((xsa4-x+(ysa4-y+(zsa4-z)。
 その後、マスター基準局選定部407は、求めた第1測位位置と測量位置との距離Δsa1、Δsa2、Δsa3及びΔsa4のそれぞれが所定の閾値(所定範囲)以内であるか否かを判定する。
 マスター基準局選定部407は、判定の結果、Δsa1、Δsa2、Δsa3及びΔsa4のうち所定の閾値以内である距離に対応する基準局を、マスター基準局として選定する。マスター基準局選定部407は、例えば、基準局301の距離であるΔsa1と基準局302の当該距離であるΔsa2が所定の範囲以内である場合、基準局301と基準局302とを、マスター基準局として選定する。
 なお、マスター基準局選定部407は、判定の結果、Δsa1、Δsa2、Δsa3及びΔsa4のうち所定の閾値(所定範囲)以内である「距離」が無い場合には、所定の閾値を増加させて(所定範囲を拡げて)、再度判定を行ってもよい。この場合において、マスター基準局選定部407は、判定を行うごとに所定の閾値を増加させ(所定範囲を拡げ)、マスター基準局を選定できるまで、当該判定を繰り返し実行してもよい。
 また、マスター基準局選定部407は、判定の結果、Δsa1、Δsa2、Δsa3及びΔsa4のうち所定の閾値(所定範囲)以内の「距離」が無い場合、第1の実施形態と同様に、第1測位位置と測量位置との距離の大小関係に基づきマスター基準局を選定してもよい。
 また、マスター基準局選定部407は、判定の結果、Δsa1、Δsa2、Δsa3及びΔsa4のうち所定の閾値(所定範囲)以内の距離が無い場合、マスター基準局を選定しなくてもよい。この場合において、補正値処理装置401は、その後の処理を実行せず、処理を終了してもよい。
 一方、マスター基準局選定部407は、判定の結果、Δsa1、Δsa2、Δsa3及びΔsa4の全てが所定の閾値(所定範囲)以内である場合、所定の閾値を減少させて(所定範囲を狭めて)、再度判定を行ってもよい。この場合において、マスター基準局選定部407は、判定を行うごとに所定の閾値を減少させ(所定範囲を狭め)、マスター基準局以外の基準局が存在するようになるまで、当該判定を繰り返し実行してもよい。
 また、マスター基準局選定部407は、判定の結果、Δsa1、Δsa2、Δsa3及びΔsa4の全てが所定の閾値(所定範囲)以内である場合、第1の実施形態と同様に、第1測位位置と測量位置との距離の大小関係に基づき、マスター基準局を選定してもよい。
 また、マスター基準局選定部407は、判定の結果、Δsa1、Δsa2、Δsa3及びΔsa4の全てが所定の閾値(所定範囲)以内である場合、全ての基準局をマスター基準局として選定してもよい。この場合において、補正値処理装置401は、マスター基準局以外の基準局について、異常状態か否かを判定する処理を実行しない。
 ここで、マスター基準局選定部407において、マスター基準局として選択するための条件である所定範囲は、例えばユーザ(管理者など)の要求に応じて変更できてもよい。この場合において、ユーザ(管理者など)は、基準局の数や配置を考慮して所定範囲を設定することにより、マスター基準局の選定の条件を調整することができる。当該所定範囲を広く設定した場合、第1測位位置と測量位置との距離がある程度大きい基準局であっても、マスター基準局として選択される。この場合、当該マスター基準局の補正値は大きくなり、その結果、マスター基準局以外の基準局の第1測位位置の補正量も大きくなる。ここで、マスター基準局以外の基準局の第2測位位置は、当該マスター基準局以外の基準局の第1測位位置を、マスター基準局の補正値で補正した値である。すなわち、マスター基準局の補正値が大きければ、マスター基準局以外の基準局の第2測位位置は、第1測位位置を大きく補正したものとなる。
 第1測位位置を大きく補正していれば、DGPS測位における誤差の補正量が大きくなるため、第1測位位置は測量位置に近づく。したがって、マスター基準局以外の基準局の測量位置と、当該第2測位位置との距離が第1の閾値よりも大きくなる可能性は少なくなる。
 異常判定部410は、例えば、第1の実施形態における表1の合否判定ロジックに基づいて行われる。この場合において、異常判定部410は、1つのマスター基準局における補正値での判定で合格であれば、総合判定で合格と判定する。したがって、第1測位位置と測量位置との距離がある程度大きい基準局がマスター基準局として選定されると、補正値が大きくなり、異常判定部410が、マスター基準局以外の基準局を合格と判定する場合が増える。その結果、異常判定部410が異常状態と判定するマスター基準局以外の基準局の数が少なくなる。言い換えると、所定範囲を広くすることにより、異常状態判定部410は、ある程度の誤差であれば、マスター基準局以外の基準局を正常と判断するようになる。
 一方、当該所定範囲を狭く設定した場合、第1測位位置と測量位置との距離が小さい基準局でなければ、マスター基準局として選択されなくなる。その結果、異常判定部410が異常状態と判定するマスター基準局以外の基準局の数が多くなる。言い換えると、所定範囲を狭めることにより、異常判定部410は、少しの誤差であっても、マスター基準局以外の基準局が異常であると判断するようになる。
 すなわち、本発明の第3の実施形態において、マスター基準局として選択するための条件である所定範囲を変更することにより、マスター基準局以外の基準局が異常であるか正常であるかの判定の精度を柔軟に調整することができる。
 図6は、本発明の第3の実施形態における、補正値処理装置401の動作例を示すフローチャートである。
 補正値処理装置401のデータ受信部402は、GPS衛星201乃至20nから発せられたGPS信号を受信した基準局301乃至304から観測データを受信する(ステップS601)。
 次に、補正値処理装置401のデータ処理部403は、観測データを基に、基準局301乃至304のそれぞれについての位置を測位(単独測位)し、第1測位距離を求める(ステップS602)。
 次に、データ処理部403は、基準局301乃至304から、基準局の予め与えられた位置情報(測量位置)と、基準局の第1測位位置との距離が所定範囲に含まれている基準局をマスター基準局として選択する(ステップS603)
 次に、データ処理部403は、マスター基準局から得た観測データと測量位置の値を基に、マスター基準局の補正値を算出する(ステップS604)。
 データ処理部403は、マスター基準局以外の基準局が、当該マスター基準局以外の基準局の第1測位位置をマスター基準局の補正値により補正することにより、当該マスター基準局以外の基準局の第2測位位置を求める(ステップS605)。
 データ処理部403は、マスター基準局以外の基準局の予め与えられた位置情報による位置と、第2測位位置との距離が第1の閾値よりも大きいときに、マスター基準局以外の基準局が異常状態であると判定する(ステップS606)。
 上記のとおり、本発明の第3の実施形態の状態検出方法は、基準局の予め与えられた位置情報(測量位置)と、基準局の第1測位位置との距離が所定範囲に含まれている基準局を、マスター基準局とする。そして、当該状態検出方法は、基準局の数や配置を考慮して当該所定範囲を設定することにより、マスター基準局の選定の条件を調整することが可能となる。その結果、当該状態検出方法は、マスター基準局以外の基準局が異常であるか正常であるかの判定の精度を、柔軟に変更することが可能となる。
 上記の実施形態の一部又は全部は、以下の付記のようにも記載されうるが、以下には限られない。
 [付記1]
 衛星と基準局とを有する測位システムにおいて基準局の異常状態を検出する状態検出方法であって、
 前記基準局が前記衛星から受信した信号に基づき、前記基準局の第1測位位置を求め、
 前記基準局の予め与えられた位置情報と、前記基準局の前記第1測位位置との間の距離が小さいものから順に所定数のマスター基準局を選択し、
 前記マスター基準局が前記衛星から受信した信号に基づき、前記マスター基準局の補正値を生成し、
 前記マスター基準局以外の前記基準局が前記マスター基準局以外の前記基準局の前記第1測位位置を前記マスター基準局の前記補正値により補正することにより、前記マスター基準局以外の基準局の第2測位位置を求め、
 前記マスター基準局以外の基準局の予め与えられた位置情報による位置と、前記第2測位位置との距離が第1の閾値よりも大きいときに、前記マスター基準局以外の基準局が異常状態であると判定することを特徴とする状態検出方法。
 [付記2]
 衛星と複数の基準局とを有する測位システムにおいて、前記複数の基準局の異常状態を検出する状態検出方法であって、
 前記基準局が前記衛星から受信した信号に基づき、前記基準局の第1測位位置を求め、
 前記基準局の予め与えられた位置情報と、前記基準局の前記第1測位位置との距離が所定範囲に含まれている基準局をマスター基準局として選択し、
 前記マスター基準局が前記衛星から受信した信号に基づき、前記マスター基準局の補正値を生成し、
 前記マスター基準局以外の前記基準局が前記マスター基準局以外の前記基準局の前記第1測位位置を前記マスター基準局の前記補正値により補正することにより、前記マスター基準局以外の基準局の第2測位位置を求め、
 前記マスター基準局以外の基準局の予め与えられた位置情報による位置と、前記第2測位位置との距離が第1の閾値よりも大きいときに、前記マスター基準局以外の基準局が異常状態であると判定する
ことを特徴とする、状態検出方法。
 [付記3]
 前記マスター基準局の補正値は、前記マスター基準局に対する電離層遅延量と対流圏遅延量とを含むことを特徴とする、付記1又は2に記載の状態検出方法。
 [付記4]
 前記マスター基準局の補正値は、前記マスター基準局の予め与えられた位置から、前記マスター基準局と前記衛星との前記疑似距離にスムージング処理を実行して得た前記第1測位位置を差し引いた値であることを特徴とする付記1乃至3のいずれかに記載の状態検出方法。
 [付記5]
 前記測位システムは、複数の前記マスター基準局を含み、
 前記マスター基準局以外の基準局の前記疑似距離の補正に用いる前記補正値は、複数の前記マスター基準局のそれぞれで生成された前記補正値の平均値であることを特徴とする付記1乃至4のいずれかに記載の状態検出方法。
 [付記6]
 前記測位システムは、第1のマスター基準局と第2のマスター基準局を有し、
 前記第1のマスター基準局が前記衛星から受信した信号を基に求めた前記第1のマスター基準局と前記衛星との擬似距離を、前記第2のマスター基準局の補正値により補正することにより得られた、前記第1のマスター基準局の第2測位位置を求め、
 前記第1のマスター基準局の予め与えられた位置情報による位置と、前記第1のマスター基準局の前記第2測位位置との距離が第2の閾値よりも大きいときに、前記第1のマスター基準局が異常状態であると判定し、
 前記第1のマスター基準局が異常状態であると判定したときに、前記複数のマスター基準局以外の基準局のうちの1つを前記第1のマスター基準局に選択する
ことを特徴とする、付記1乃至5のいずれかに記載の状態検出方法。
 [付記7]
 複数の前記衛星から、受信する信号を送信する衛星を選択することを特徴とする、付記1乃至6のいずれかに記載の状態検出方法。
 [付記8]
 前記受信する信号を送信する衛星は、複数の前記衛星のうち、仰角が所定の角度よりも大きい衛星、又は、前記測位システムに含まれる全ての基準局から可視である衛星から選択されることを特徴とする、付記7に記載の状態検出方法。
 [付記9]
 前記第1の閾値を、所定の確率分布関数を参照して決定することを特徴とする、付記1乃至8のいずれかに記載の状態検出方法。
 [付記10]
 前記第1の閾値を、前記衛星の数と、前記衛星の配置状況とに基づいて決定することを特徴とする、付記1乃至9のいずれかに記載の状態検出方法。
 [付記11]
 基準局が衛星から受信した信号を受信するデータ受信手段と、
 前記データ受信手段が受信した前記信号基づき、前記基準局の第1測位位置を算出する第1測位位置算出手段と、
 前記基準局の予め与えられた位置情報と、前記基準局の前記第1測位位置との間の距離が小さいものから順に所定数のマスター基準局を選択するマスター基準局選択手段と、
 前記マスター基準局が前記衛星から受信した信号に基づき、前記マスター基準局の補正値を生成するマスター基準局補正値算出手段と、
 前記マスター基準局以外の前記基準局が前記マスター基準局以外の前記基準局の前記第1測位位置を前記マスター基準局の前記補正値により補正することにより、前記マスター基準局以外の基準局の第2測位位置を算出する第2測位位置算出手段と、
 前記マスター基準局以外の基準局の予め与えられた位置情報による位置と、前記第2測位位置との距離が第1の閾値よりも大きいときに、前記マスター基準局以外の基準局が異常状態であると判定する異常判定手段と、
を有することを特徴とする、補正値処理装置。
 [付記12]
 基準局が衛星から受信した信号を受信するデータ受信手段と、
 前記データ受信手段が受信した前記信号基づき、前記基準局の第1測位位置を算出する第1測位位置算出手段と、
 前記基準局の予め与えられた位置情報と、前記基準局の前記第1測位位置との距離が所定範囲に含まれている基準局をマスター基準局として選択するマスター基準局選択手段と、
 前記マスター基準局が前記衛星から受信した信号に基づき、前記マスター基準局の補正値を生成するマスター基準局補正値算出手段と、
 前記マスター基準局以外の前記基準局が前記マスター基準局以外の前記基準局の前記第1測位位置を前記マスター基準局の前記補正値により補正することにより、前記マスター基準局以外の基準局の第2測位位置を算出する第2測位位置算出手段と、
 前記マスター基準局以外の基準局の予め与えられた位置情報による位置と、前記第2測位位置との距離が第1の閾値よりも大きいときに、前記マスター基準局以外の基準局が異常状態であると判定する異常判定手段と、
を有することを特徴とする、補正値処理装置。
 [付記13]
 前記マスター基準局の補正値は、前記マスター基準局に対する電離層遅延量と対流圏遅延量とを含むことを特徴とする、付記11又は12に記載の補正値処理装置。
 [付記14]
 前記データ処理手段で生成した前記補正値を、前記衛星から受信した信号を基に自身の位置を推定する受信装置に送信する補正値送信手段をさらに備えることを特徴とする、付記11乃至13のいずれかに記載の補正値処理装置。
 [付記15]
 衛星と、基準局と、補正値処理装置を有する測位システムであって、
 前記補正値処理装置は、
 前記データ受信手段が受信した前記信号基づき、前記基準局の第1測位位置を算出する第1測位位置算出手段と、
 前記基準局の予め与えられた位置情報と、前記基準局の前記第1測位位置との間の距離が小さいものから順に所定数のマスター基準局を選択するマスター基準局選択手段と、
 前記マスター基準局が前記衛星から受信した信号に基づき、前記マスター基準局の補正値を生成するマスター基準局補正値算出手段と、
 前記マスター基準局以外の前記基準局が前記マスター基準局以外の前記基準局の前記第1測位位置を前記マスター基準局の前記補正値により補正することにより、前記マスター基準局以外の基準局の第2測位位置を算出する第2測位位置算出手段と、
 前記マスター基準局以外の基準局の予め与えられた位置情報による位置と、前記第2測位位置との距離が第1の閾値よりも大きいときに、前記マスター基準局以外の基準局が異常状態であると判定する異常判定手段と、
を有することを特徴とする、測位システム。
 [付記16]
 衛星と、基準局と、補正値処理装置を有する測位システムであって、
 前記補正値処理装置は、
 前記データ受信手段が受信した前記信号基づき、前記基準局の第1測位位置を算出する第1測位位置算出手段と、
 前記基準局の予め与えられた位置情報と、前記基準局の前記第1測位位置との距離が所定範囲に含まれている基準局をマスター基準局として選択するマスター基準局選択手段と、
 前記マスター基準局が前記衛星から受信した信号に基づき、前記マスター基準局の補正値を生成するマスター基準局補正値算出手段と、
 前記マスター基準局以外の前記基準局が前記マスター基準局以外の前記基準局の前記第1測位位置を前記マスター基準局の前記補正値により補正することにより、前記マスター基準局以外の基準局の第2測位位置を算出する第2測位位置算出手段と、
 前記マスター基準局以外の基準局の予め与えられた位置情報による位置と、前記第2測位位置との距離が第1の閾値よりも大きいときに、前記マスター基準局以外の基準局が異常状態であると判定する異常判定手段と、
を有することを特徴とする、測位システム。
 [付記17]
 前記マスター基準局の補正値は、前記マスター基準局に対する電離層遅延量と対流圏遅延量とを含むことを特徴とする、付記15又は16に記載の測位システム。
 [付記18]
 前記衛星から受信した信号を基に自身の位置を推定する受信装置をさらに有し、
 前記補正値処理装置は、前記データ処理部で生成した前記補正値を、前記受信装置に送信する補正値送信手段をさらに備えることを特徴とする、付記15乃至17のいずれかに記載の測位システム。
 [付記19]
 衛星と基準局とを有する測位システムにおいて、前記基準局の異常状態を検出するプログラムであって、
 前記データ受信手段が受信した前記信号基づき、前記基準局の第1測位位置を算出する処理と、
 前記基準局の予め与えられた位置情報と、前記基準局の前記第1測位位置との間の距離が小さいものから順に所定数のマスター基準局を選択する処理と、
 前記マスター基準局が前記衛星から受信した信号に基づき、前記マスター基準局の補正値を生成する処理と、
 前記マスター基準局以外の前記基準局が前記マスター基準局以外の前記基準局の前記第1測位位置を前記マスター基準局の前記補正値により補正することにより、前記マスター基準局以外の基準局の第2測位位置を算出する第2測位位置算出処理と、
 前記マスター基準局以外の基準局の予め与えられた位置情報による位置と、前記第2測位位置との距離が第1の閾値よりも大きいときに、前記マスター基準局以外の基準局が異常状態であると判定する処理と、
とを含むことを特徴とするプログラムを記憶した記憶媒体。
 [付記20]
 衛星と基準局とを有する測位システムにおいて、前記基準局の異常状態を検出するプログラムであって、
 前記データ受信手段が受信した前記信号基づき、前記基準局の第1測位位置を算出する処理と、
 前記基準局の予め与えられた位置情報と、前記基準局の前記第1測位位置との距離が所定範囲に含まれている基準局をマスター基準局として選択する処理と、
 前記マスター基準局が前記衛星から受信した信号に基づき、前記マスター基準局の補正値を生成する処理と、
 前記マスター基準局以外の前記基準局が前記マスター基準局以外の前記基準局の前記第1測位位置を前記マスター基準局の前記補正値により補正することにより、前記マスター基準局以外の基準局の第2測位位置を算出する第2測位位置算出処理と、
 前記マスター基準局以外の基準局の予め与えられた位置情報による位置と、前記第2測位位置との距離が第1の閾値よりも大きいときに、前記マスター基準局以外の基準局が異常状態であると判定する処理と、
とを含むことを特徴とするプログラムを記憶した記憶媒体。
 [付記21]
 前記マスター基準局の補正値は、少なくとも前記マスター基準局に対する電離層遅延量と対流圏遅延量とを含むことを特徴とする、付記19又は20に記載のプログラムを記憶した記憶媒体。
 [付記22]
 衛星と基準局とを有する測位システムにおいて、前記基準局の異常状態を検出する状態検出方法であって、
 前記基準局が前記衛星から受信した信号を基に求めた前記基準局と前記衛星との擬似距離を基に、前記基準局の第1の測位位置を求め、
 前記基準局の予め与えられた位置情報による位置と、前記基準局の前記第1の測位位置との間の距離が小さいものから順に所定数のマスター基準局を選択し、
 前記マスター基準局が前記衛星から受信した信号を基に、前記マスター基準局の補正値を生成し、
 前記マスター基準局以外の基準局の前記擬似距離を、前記マスター基準局の前記補正値により補正することにより、前記マスター基準局以外の基準局の第2の測位位置を求め、
 前記マスター基準局以外の基準局の予め与えられた位置情報による位置と、前記マスター基準局以外の基準局の前記第2の測位位置との距離が第1の閾値よりも大きいときに、前記マスター基準局以外の基準局が異常状態であると判定する
ことを特徴とする、状態検出方法。
 [付記23]
 前記マスター基準局の補正値は、前記マスター基準局に対する電離層遅延量と対流圏遅延量とを含むことを特徴とする、付記22に記載の状態検出方法。
 [付記24]
 前記マスター基準局の補正値は、前記マスター基準局の予め与えられた位置から、前記マスター基準局と前記衛星との前記疑似距離にスムージング処理を実行して得た前記第1の測位位置を差し引いた値であることを特徴とする付記22又は23に記載の状態検出方法。
 [付記25]
 前記測位システムは、複数の前記マスター基準局を含み、
 前記マスター基準局以外の基準局の前記疑似距離の補正に用いる前記補正値は、複数の前記マスター基準局のそれぞれで生成された前記補正値の平均値であることを特徴とする付記22乃至24のいずれかに記載の状態検出方法。
 [付記26]
 前記測位システムは、第1のマスター基準局と第2のマスター基準局を有し、
 前記第1のマスター基準局が前記衛星から受信した信号を基に求めた前記第1のマスター基準局と前記衛星との擬似距離を、前記第2のマスター基準局の補正値により補正することにより得られた、前記第1のマスター基準局の第2の測位位置を求め、
 前記第1のマスター基準局の予め与えられた位置情報による位置と、前記第1のマスター基準局の前記第2の測位位置との距離が第2の閾値よりも大きいときに、前記第1のマスター基準局が異常状態であると判定し、
 前記第1のマスター基準局が異常状態であると判定したときに、前記複数のマスター基準局以外の基準局のうちの1つを前記第1のマスター基準局に選択する
ことを特徴とする、付記22乃至25のいずれかに記載の状態検出方法。
 [付記27]
 前記基準局の第1の測位位置を求めるときに、複数の前記衛星から、受信する信号を送信する衛星を選択する
ことを特徴とする、付記22乃至26のいずれかに記載の状態検出方法。
 [付記28]
 前記受信する信号を送信する衛星は、複数の前記衛星のうち、仰角が所定の角度よりも大きい衛星、又は、前記測位システムに含まれる全ての基準局から可視である衛星から選択されることを特徴とする、付記27に記載の状態検出方法。
 [付記29]
 前記第1の閾値を、所定の確率分布関数を参照して決定することを特徴とする、付記22乃至28のいずれかに記載の状態検出方法。
 [付記30]
 前記第1の閾値を、前記衛星の数と、前記衛星の配置状況とに基づいて決定することを特徴とする、付記22乃至29のいずれかに記載の状態検出方法。
 [付記31]
 基準局が衛星から受信した信号を受信するデータ受信手段と、
 前記データ受信手段が受信した前記信号を基に求めた前記基準局と前記衛星との擬似距離を基に、前記基準局の第1の測位位置を求め、前記基準局の予め与えられた位置情報による位置と、前記基準局の前記第1の測位位置との間の距離が小さいものから順に所定数のマスター基準局を選択し、前記マスター基準局が前記衛星から受信した信号を基に、前記マスター基準局の補正値を生成し、前記マスター基準局以外の基準局の前記擬似距離を、前記マスター基準局の前記補正値により補正することにより、前記マスター基準局以外の基準局の第2の測位位置を求め、前記マスター基準局以外の基準局の予め与えられた位置情報による位置と、前記マスター基準局以外の基準局の前記第2の測位位置との距離が第1の閾値よりも大きいときに、前記マスター基準局以外の基準局が異常状態であると判定するデータ処理手段と
を有することを特徴とする、補正値処理装置。
 [付記32]
 前記マスター基準局の補正値は、前記マスター基準局に対する電離層遅延量と対流圏遅延量とを含むことを特徴とする、付記31に記載の補正値処理装置。
 [付記33]
 前記データ処理手段で生成した前記補正値を、前記衛星から受信した信号を基に自身の位置を推定する受信装置に送信する補正値送信手段をさらに備えることを特徴とする、付記31又は32に記載の補正値処理装置。
 [付記34]
 衛星と、基準局と、補正値処理装置を有する測位システムであって、
 前記補正値処理装置は、
 前記基準局が前記衛星から受信した信号を基に求めた前記基準局と前記衛星との擬似距離を基に前記基準局の第1の測位位置を求め、前記基準局の予め与えられた位置情報による位置と前記基準局の前記第1の測位位置との間の距離が小さいものから順に所定数のマスター基準局を選択し、前記マスター基準局が前記衛星から受信した信号を基に、前記マスター基準局の補正値を生成し、前記マスター基準局以外の基準局の前記擬似距離を、前記マスター基準局の前記補正値により補正することにより、前記マスター基準局以外の基準局の第2の測位位置を求め、前記マスター基準局以外の基準局の予め与えられた位置情報による位置と前記マスター基準局以外の基準局の前記第2の測位位置との距離が第1の閾値よりも大きいときに前記マスター基準局以外の基準局が異常状態であると判定する、データ処理手段
を有することを特徴とする、測位システム。
 [付記35]
 前記マスター基準局の補正値は、前記マスター基準局に対する電離層遅延量と対流圏遅延量とを含むことを特徴とする、付記34に記載の測位システム。
 [付記36]
 前記衛星から受信した信号を基に自身の位置を推定する受信装置をさらに有し、
 前記補正値処理装置は、前記データ処理部で生成した前記補正値を、前記受信装置に送信する補正値送信手段をさらに備えることを特徴とする、付記34又は35に記載の測位システム。
 [付記37]
 衛星と基準局とを有する測位システムにおいて、前記基準局の異常状態を検出するプログラムであって、
 前記基準局が前記衛星から受信した信号を基に求めた前記基準局と前記衛星との擬似距離を基に、前記基準局の第1の測位位置を求める処理と、
 前記基準局の予め与えられた位置情報による位置と、前記基準局の前記第1の測位位置との間の距離が小さいものから順に所定数のマスター基準局を選択する処理と、
 前記マスター基準局が前記衛星から受信した信号を基に、前記マスター基準局の補正値を生成する処理と、
 前記マスター基準局以外の基準局の前記擬似距離を、前記マスター基準局の前記補正値により補正することにより、前記マスター基準局以外の基準局の第2の測位位置を求める処理と、
 前記マスター基準局以外の基準局の予め与えられた位置情報による位置と、前記マスター基準局以外の基準局の前記第2の測位位置との距離が第1の閾値よりも大きいときに、前記マスター基準局以外の基準局が異常状態であると判定する処理
とを含むことを特徴とするプログラムを記憶した記憶媒体。
 [付記38]
 前記マスター基準局の補正値は、少なくとも前記マスター基準局に対する電離層遅延量と対流圏遅延量とを含むことを特徴とする、付記37に記載のプログラムを記憶した記憶媒体。
 本発明の状態検出方法等は、上記実施形態に基づき説明されているが、上記実施形態に限定されることなく、本発明の範囲内においてかつ本発明の基本的技術思想に基づき、上記実施形態に対し種々の変形、変更及び改良を含むことができることはいうまでもない。また、本発明の請求の範囲の枠内において、種々の開示要素の多様な組み合わせ・置換ないし選択が可能である。本発明のさらなる課題、目的及び展開形態は、請求の範囲を含む本発明の全開示事項からも明らかにされる。
 この出願は、2013年2月26日に出願された日本出願特願2013-035783を基礎とする優先権を主張し、その開示の全てをここに取り込む。
 本発明は、上記の実施形態に限定されるものではなく、複数の基準局を使用する測位システムに好適に適用できる。
 101  ディファレンシャルGPS
 201、20n  GPS衛星
 301、302、303、304  基準局
 401  補正値処理装置
 402  データ受信部
 403  データ処理部
 404  データ保持部
 405  補正値送信部
 406  第1測位位置算出部
 407  マスター基準局選定部
 408  マスター基準局補正値算出部
 409  第2測位位置算出部
 410  異常判定部
 501  放送送信装置
 601  受信装置

Claims (21)

  1.  衛星と基準局とを有する測位システムにおいて基準局の異常状態を検出する状態検出方法であって、
     前記基準局が前記衛星から受信した信号に基づき、前記基準局の第1測位位置を求め、
     前記基準局の予め与えられた位置情報と、前記基準局の前記第1測位位置との間の距離が小さいものから順に所定数のマスター基準局を選択し、
     前記マスター基準局が前記衛星から受信した信号に基づき、前記マスター基準局の補正値を生成し、
     前記マスター基準局以外の前記基準局が前記マスター基準局以外の前記基準局の前記第1測位位置を前記マスター基準局の前記補正値により補正することにより、前記マスター基準局以外の基準局の第2測位位置を求め、
     前記マスター基準局以外の基準局の予め与えられた位置情報による位置と、前記第2測位位置との距離が第1の閾値よりも大きいときに、前記マスター基準局以外の基準局が異常状態であると判定することを特徴とする状態検出方法。
  2.  衛星と複数の基準局とを有する測位システムにおいて、前記複数の基準局の異常状態を検出する状態検出方法であって、
     前記基準局が前記衛星から受信した信号に基づき、前記基準局の第1測位位置を求め、
     前記基準局の予め与えられた位置情報と、前記基準局の前記第1測位位置との距離が所定範囲に含まれている基準局をマスター基準局として選択し、
     前記マスター基準局が前記衛星から受信した信号に基づき、前記マスター基準局の補正値を生成し、
     前記マスター基準局以外の前記基準局が前記マスター基準局以外の前記基準局の前記第1測位位置を前記マスター基準局の前記補正値により補正することにより、前記マスター基準局以外の基準局の第2測位位置を求め、
     前記マスター基準局以外の基準局の予め与えられた位置情報による位置と、前記第2測位位置との距離が第1の閾値よりも大きいときに、前記マスター基準局以外の基準局が異常状態であると判定する
    を有することを特徴とする、状態検出方法。
  3.  前記マスター基準局の補正値は、前記マスター基準局に対する電離層遅延量と対流圏遅延量とを含むことを特徴とする、請求項1又は2に記載の状態検出方法。
  4.  前記マスター基準局の補正値は、前記マスター基準局の予め与えられた位置から、前記マスター基準局と前記衛星との前記疑似距離にスムージング処理を実行して得た前記第1測位位置を差し引いた値であることを特徴とする請求項1乃至3のいずれかに記載の状態検出方法。
  5.  前記測位システムは、複数の前記マスター基準局を含み、
     前記マスター基準局以外の基準局の前記疑似距離の補正に用いる前記補正値は、複数の前記マスター基準局のそれぞれで生成された前記補正値の平均値であることを特徴とする請求項1乃至4のいずれかに記載の状態検出方法。
  6.  前記測位システムは、第1のマスター基準局と第2のマスター基準局を有し、
     前記第1のマスター基準局が前記衛星から受信した信号を基に求めた前記第1のマスター基準局と前記衛星との擬似距離を、前記第2のマスター基準局の補正値により補正することにより得られた、前記第1のマスター基準局の第2測位位置を求め、
     前記第1のマスター基準局の予め与えられた位置情報による位置と、前記第1のマスター基準局の前記第2測位位置との距離が第2の閾値よりも大きいときに、前記第1のマスター基準局が異常状態であると判定し、
     前記第1のマスター基準局が異常状態であると判定したときに、前記複数のマスター基準局以外の基準局のうちの1つを前記第1のマスター基準局に選択する
    ことを特徴とする、請求項1乃至5のいずれかに記載の状態検出方法。
  7.  複数の前記衛星から、受信する信号を送信する衛星を選択することを特徴とする、請求項1乃至6のいずれかに記載の状態検出方法。
  8.  前記受信する信号を送信する衛星は、複数の前記衛星のうち、仰角が所定の角度よりも大きい衛星、又は、前記測位システムに含まれる全ての基準局から可視である衛星から選択されることを特徴とする、請求項7に記載の状態検出方法。
  9.  前記第1の閾値を、所定の確率分布関数を参照して決定することを特徴とする、請求項1乃至8のいずれかに記載の状態検出方法。
  10.  前記第1の閾値を、前記衛星の数と、前記衛星の配置状況とに基づいて決定することを特徴とする、請求項1乃至9のいずれかに記載の状態検出方法。
  11.  基準局が衛星から受信した信号を受信するデータ受信手段と、
     前記データ受信手段が受信した前記信号基づき、前記基準局の第1測位位置を算出する第1測位位置算出手段と、
     前記基準局の予め与えられた位置情報と、前記基準局の前記第1測位位置との間の距離が小さいものから順に所定数のマスター基準局を選択するマスター基準局選択手段と、
     前記マスター基準局が前記衛星から受信した信号に基づき、前記マスター基準局の補正値を生成するマスター基準局補正値算出手段と、
     前記マスター基準局以外の前記基準局が前記マスター基準局以外の前記基準局の前記第1測位位置を前記マスター基準局の前記補正値により補正することにより、前記マスター基準局以外の基準局の第2測位位置を算出する第2測位位置算出手段と、
     前記マスター基準局以外の基準局の予め与えられた位置情報による位置と、前記第2測位位置との距離が第1の閾値よりも大きいときに、前記マスター基準局以外の基準局が異常状態であると判定する異常判定手段と、
    を有することを特徴とする、補正値処理装置。
  12.  基準局が衛星から受信した信号を受信するデータ受信手段と、
     前記データ受信手段が受信した前記信号基づき、前記基準局の第1測位位置を算出する第1測位位置算出手段と、
     前記基準局の予め与えられた位置情報と、前記基準局の前記第1測位位置との距離が所定範囲に含まれている基準局をマスター基準局として選択するマスター基準局選択手段と、
     前記マスター基準局が前記衛星から受信した信号に基づき、前記マスター基準局の補正値を生成するマスター基準局補正値算出手段と、
     前記マスター基準局以外の前記基準局が前記マスター基準局以外の前記基準局の前記第1測位位置を前記マスター基準局の前記補正値により補正することにより、前記マスター基準局以外の基準局の第2測位位置を算出する第2測位位置算出手段と、
     前記マスター基準局以外の基準局の予め与えられた位置情報による位置と、前記第2測位位置との距離が第1の閾値よりも大きいときに、前記マスター基準局以外の基準局が異常状態であると判定する異常判定手段と、
    を有することを特徴とする、補正値処理装置。
  13.  前記マスター基準局の補正値は、前記マスター基準局に対する電離層遅延量と対流圏遅延量とを含むことを特徴とする、請求項11又は12に記載の補正値処理装置。
  14.  前記データ処理手段で生成した前記補正値を、前記衛星から受信した信号を基に自身の位置を推定する受信装置に送信する補正値送信手段をさらに備えることを特徴とする、請求項11乃至13のいずれかに記載の補正値処理装置。
  15.  衛星と、基準局と、補正値処理装置を有する測位システムであって、
     前記補正値処理装置は、
     前記データ受信手段が受信した前記信号基づき、前記基準局の第1測位位置を算出する第1測位位置算出手段と、
     前記基準局の予め与えられた位置情報と、前記基準局の前記第1測位位置との間の距離が小さいものから順に所定数のマスター基準局を選択するマスター基準局選択手段と、
     前記マスター基準局が前記衛星から受信した信号に基づき、前記マスター基準局の補正値を生成するマスター基準局補正値算出手段と、
     前記マスター基準局以外の前記基準局が前記マスター基準局以外の前記基準局の前記第1測位位置を前記マスター基準局の前記補正値により補正することにより、前記マスター基準局以外の基準局の第2測位位置を算出する第2測位位置算出手段と、
     前記マスター基準局以外の基準局の予め与えられた位置情報による位置と、前記第2測位位置との距離が第1の閾値よりも大きいときに、前記マスター基準局以外の基準局が異常状態であると判定する異常判定手段と、
    を有することを特徴とする、測位システム。
  16.  衛星と、基準局と、補正値処理装置を有する測位システムであって、
     前記補正値処理装置は、
     前記データ受信手段が受信した前記信号基づき、前記基準局の第1測位位置を算出する第1測位位置算出手段と、
     前記基準局の予め与えられた位置情報と、前記基準局の前記第1測位位置との距離が所定範囲に含まれている基準局をマスター基準局として選択するマスター基準局選択手段と、
     前記マスター基準局が前記衛星から受信した信号に基づき、前記マスター基準局の補正値を生成するマスター基準局補正値算出手段と、
     前記マスター基準局以外の前記基準局が前記マスター基準局以外の前記基準局の前記第1測位位置を前記マスター基準局の前記補正値により補正することにより、前記マスター基準局以外の基準局の第2測位位置を算出する第2測位位置算出手段と、
     前記マスター基準局以外の基準局の予め与えられた位置情報による位置と、前記第2測位位置との距離が第1の閾値よりも大きいときに、前記マスター基準局以外の基準局が異常状態であると判定する異常判定手段と、
    を有することを特徴とする、測位システム。
  17.  前記マスター基準局の補正値は、前記マスター基準局に対する電離層遅延量と対流圏遅延量とを含むことを特徴とする、請求項15又は16に記載の測位システム。
  18.  前記衛星から受信した信号を基に自身の位置を推定する受信装置をさらに有し、
     前記補正値処理装置は、前記データ処理部で生成した前記補正値を、前記受信装置に送信する補正値送信手段をさらに備えることを特徴とする、請求項15乃至17のいずれかに記載の測位システム。
  19.  衛星と基準局とを有する測位システムにおいて、前記基準局の異常状態を検出するプログラムであって、
     前記データ受信手段が受信した前記信号基づき、前記基準局の第1測位位置を算出する処理と、
     前記基準局の予め与えられた位置情報と、前記基準局の前記第1測位位置との間の距離が小さいものから順に所定数のマスター基準局を選択する処理と、
     前記マスター基準局が前記衛星から受信した信号に基づき、前記マスター基準局の補正値を生成する処理と、
     前記マスター基準局以外の前記基準局が前記マスター基準局以外の前記基準局の前記第1測位位置を前記マスター基準局の前記補正値により補正することにより、前記マスター基準局以外の基準局の第2測位位置を算出する第2測位位置算出処理と、
     前記マスター基準局以外の基準局の予め与えられた位置情報による位置と、前記第2測位位置との距離が第1の閾値よりも大きいときに、前記マスター基準局以外の基準局が異常状態であると判定する処理と、
    とを含むことを特徴とするプログラムを記憶した記憶媒体。
  20.  衛星と基準局とを有する測位システムにおいて、前記基準局の異常状態を検出するプログラムであって、
     前記データ受信手段が受信した前記信号基づき、前記基準局の第1測位位置を算出する処理と、
     前記基準局の予め与えられた位置情報と、前記基準局の前記第1測位位置との距離が所定範囲に含まれている基準局をマスター基準局として選択する処理と、
     前記マスター基準局が前記衛星から受信した信号に基づき、前記マスター基準局の補正値を生成する処理と、
     前記マスター基準局以外の前記基準局が前記マスター基準局以外の前記基準局の前記第1測位位置を前記マスター基準局の前記補正値により補正することにより、前記マスター基準局以外の基準局の第2測位位置を算出する第2測位位置算出処理と、
     前記マスター基準局以外の基準局の予め与えられた位置情報による位置と、前記第2測位位置との距離が第1の閾値よりも大きいときに、前記マスター基準局以外の基準局が異常状態であると判定する処理と、
    とを含むことを特徴とするプログラムを記憶した記憶媒体。
  21.  前記マスター基準局の補正値は、少なくとも前記マスター基準局に対する電離層遅延量と対流圏遅延量とを含むことを特徴とする、請求項19又は20に記載のプログラムを記憶した記憶媒体。
PCT/JP2014/000964 2013-02-26 2014-02-25 状態検出方法、補正値処理装置、測位システム、および状態検出プログラム WO2014132618A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN201480010469.1A CN105008956A (zh) 2013-02-26 2014-02-25 状态检测方法、校正值处理设备、定位系统和状态检测程序
US14/764,476 US20150362596A1 (en) 2013-02-26 2014-02-25 State detecting method, correction value processing device, positioning system, and state detection program
JP2015502766A JP6119838B2 (ja) 2013-02-26 2014-02-25 状態検出方法、補正値処理装置、測位システム、および状態検出プログラム

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2013035783 2013-02-26
JP2013-035783 2013-02-26

Publications (1)

Publication Number Publication Date
WO2014132618A1 true WO2014132618A1 (ja) 2014-09-04

Family

ID=51427895

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2014/000964 WO2014132618A1 (ja) 2013-02-26 2014-02-25 状態検出方法、補正値処理装置、測位システム、および状態検出プログラム

Country Status (4)

Country Link
US (1) US20150362596A1 (ja)
JP (1) JP6119838B2 (ja)
CN (1) CN105008956A (ja)
WO (1) WO2014132618A1 (ja)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106249264A (zh) * 2016-07-19 2016-12-21 深圳市安煋信息技术有限公司 一种定位方法、系统及导航装置
CN107290764A (zh) * 2016-04-13 2017-10-24 苏州宝时得电动工具有限公司 差分全球定位系统及其定位方法
JP2018084533A (ja) * 2016-11-25 2018-05-31 エヌ・ティ・ティ・データ・カスタマサービス株式会社 位置補正情報提供システム及び位置補正情報提供方法
WO2019049540A1 (ja) * 2017-09-07 2019-03-14 ヤンマー株式会社 測位システム
WO2019049541A1 (ja) * 2017-09-07 2019-03-14 ヤンマー株式会社 測位システム
JP2020012779A (ja) * 2018-07-20 2020-01-23 古野電気株式会社 測位装置、測位方法、および、測位プログラム
JP2020180857A (ja) * 2019-04-25 2020-11-05 富士通株式会社 測位システム、測位装置及び測位方法
CN113762558A (zh) * 2020-06-03 2021-12-07 千寻位置网络有限公司 基于地基增强系统的选站模型优化方法和选站方法
JP7457165B2 (ja) 2020-08-14 2024-03-27 中移(上海)信息通信科技有限公司 基地局決定方法、装置、機器及びコンピュータ記憶媒体

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9903952B2 (en) * 2015-03-18 2018-02-27 Amazon Technologies, Inc. GPS error correction via network of fixed point ground stations
DE102017206280A1 (de) * 2017-04-12 2018-10-18 Robert Bosch Gmbh Verfahren zum Überwachen einer Integrität von Referenzstationen eines Korrekturdienstsystems, Korrekturdienstsystem, Verfahren zum Betreiben eines satellitengestützten Navigationssystems und satellitengestütztes Navigationssystem
CN107247280B (zh) * 2017-04-21 2020-01-24 中国科学院光电研究院 一种定位认证和处理方法及装置
CN108802764B (zh) * 2017-04-28 2021-11-30 千寻位置网络有限公司 卫星地基增强系统的自检系统的构建方法和构建系统
EP3422780B1 (en) 2017-07-01 2020-05-13 Ruckus Wireless, Inc. Identifying a synchronization master for radio nodes
CN108020854B (zh) * 2017-10-26 2022-05-17 广州中南民航空管技术装备工程有限公司 一种场面目标态势显示方法及系统
CN108089214B (zh) * 2017-12-20 2021-06-15 北京卫星导航中心 一种卫星定位方法和卫星定位系统
KR102136725B1 (ko) * 2018-07-30 2020-07-23 전자부품연구원 다변측정 감지시스템의 중앙처리장치 및 이를 이용한 항공기 위치 산출방법
CN111279221A (zh) * 2018-12-20 2020-06-12 深圳市大疆创新科技有限公司 一种基准站的位置标定方法及装置
CN109633544B (zh) * 2018-12-26 2021-04-06 奇点新源国际技术开发(北京)有限公司 一种锚点坐标标定方法、锚点定位方法及装置
CN109597099B (zh) * 2018-12-26 2022-04-01 上海司南卫星导航技术股份有限公司 判断基准站接收机是否移动的方法、oem板卡及接收机
KR102488643B1 (ko) * 2019-09-11 2023-01-16 한국도로공사 측위 차이를 이용한 정밀 위치 보정 장치 및 그 방법

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09251068A (ja) * 1996-03-14 1997-09-22 Japan Radio Co Ltd ディファレンシャルgps基準局
JPH09274074A (ja) * 1996-04-08 1997-10-21 Nec Corp Dgps
JPH11513112A (ja) * 1995-06-30 1999-11-09 ハネウエル・インコーポレーテッド 整合性監視付き差動衛星位置把握システム地上局
JP2003219448A (ja) * 2002-01-25 2003-07-31 Hitachi Ltd 位置補正情報配信システム
US20060214844A1 (en) * 2005-02-15 2006-09-28 Fagan John E Navigation system using external monitoring
JP2006322753A (ja) * 2005-05-17 2006-11-30 Fuji Heavy Ind Ltd 車両の走行制御装置
JP2007500845A (ja) * 2003-07-30 2007-01-18 ナヴコム テクノロジー インコーポレイテッド 広域またはグローバルディファレンシャルgpsシステムのためのクロック補正値を生成するための方法
JP2010216814A (ja) * 2009-03-13 2010-09-30 Fujitsu Ltd 測位システムおよび測位方法

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5596328A (en) * 1994-08-23 1997-01-21 Honeywell Inc. Fail-safe/fail-operational differential GPS ground station system
US5638077A (en) * 1995-05-04 1997-06-10 Rockwell International Corporation Differential GPS for fleet base stations with vector processing mechanization
JP2917873B2 (ja) * 1995-10-12 1999-07-12 日本電気株式会社 差分自己位置標定方法及び基準局装置、ユーザ局装置
JP2003057327A (ja) * 2001-08-09 2003-02-26 Matsushita Electric Ind Co Ltd 航法衛星信号受信機
JP4723932B2 (ja) * 2005-06-29 2011-07-13 株式会社東芝 測位システム
JP2011043449A (ja) * 2009-08-24 2011-03-03 Seiko Epson Corp 衛星信号受信装置および衛星信号受信装置の制御方法
JP5424338B2 (ja) * 2010-03-18 2014-02-26 日本電気株式会社 衛星測位システムの異常値検出装置、異常値検出方法及び異常値検出プログラム
JP5609247B2 (ja) * 2010-05-07 2014-10-22 日本電気株式会社 監視局、制御方法、広域補強システム、および制御プログラム

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11513112A (ja) * 1995-06-30 1999-11-09 ハネウエル・インコーポレーテッド 整合性監視付き差動衛星位置把握システム地上局
JPH09251068A (ja) * 1996-03-14 1997-09-22 Japan Radio Co Ltd ディファレンシャルgps基準局
JPH09274074A (ja) * 1996-04-08 1997-10-21 Nec Corp Dgps
JP2003219448A (ja) * 2002-01-25 2003-07-31 Hitachi Ltd 位置補正情報配信システム
JP2007500845A (ja) * 2003-07-30 2007-01-18 ナヴコム テクノロジー インコーポレイテッド 広域またはグローバルディファレンシャルgpsシステムのためのクロック補正値を生成するための方法
US20060214844A1 (en) * 2005-02-15 2006-09-28 Fagan John E Navigation system using external monitoring
JP2006322753A (ja) * 2005-05-17 2006-11-30 Fuji Heavy Ind Ltd 車両の走行制御装置
JP2010216814A (ja) * 2009-03-13 2010-09-30 Fujitsu Ltd 測位システムおよび測位方法

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107290764B (zh) * 2016-04-13 2021-01-05 苏州宝时得电动工具有限公司 差分全球定位系统及其定位方法
CN107290764A (zh) * 2016-04-13 2017-10-24 苏州宝时得电动工具有限公司 差分全球定位系统及其定位方法
CN106249264A (zh) * 2016-07-19 2016-12-21 深圳市安煋信息技术有限公司 一种定位方法、系统及导航装置
JP2018084533A (ja) * 2016-11-25 2018-05-31 エヌ・ティ・ティ・データ・カスタマサービス株式会社 位置補正情報提供システム及び位置補正情報提供方法
WO2019049540A1 (ja) * 2017-09-07 2019-03-14 ヤンマー株式会社 測位システム
JP2019045457A (ja) * 2017-09-07 2019-03-22 ヤンマー株式会社 測位システム
JP2019045456A (ja) * 2017-09-07 2019-03-22 ヤンマー株式会社 測位システム
WO2019049541A1 (ja) * 2017-09-07 2019-03-14 ヤンマー株式会社 測位システム
JP2020012779A (ja) * 2018-07-20 2020-01-23 古野電気株式会社 測位装置、測位方法、および、測位プログラム
JP2020180857A (ja) * 2019-04-25 2020-11-05 富士通株式会社 測位システム、測位装置及び測位方法
JP7295400B2 (ja) 2019-04-25 2023-06-21 富士通株式会社 測位システム、測位装置及び測位方法
CN113762558A (zh) * 2020-06-03 2021-12-07 千寻位置网络有限公司 基于地基增强系统的选站模型优化方法和选站方法
CN113762558B (zh) * 2020-06-03 2023-07-11 千寻位置网络有限公司 基于地基增强系统的选站模型优化方法和选站方法
JP7457165B2 (ja) 2020-08-14 2024-03-27 中移(上海)信息通信科技有限公司 基地局決定方法、装置、機器及びコンピュータ記憶媒体

Also Published As

Publication number Publication date
CN105008956A (zh) 2015-10-28
JPWO2014132618A1 (ja) 2017-02-02
JP6119838B2 (ja) 2017-04-26
US20150362596A1 (en) 2015-12-17

Similar Documents

Publication Publication Date Title
JP6119838B2 (ja) 状態検出方法、補正値処理装置、測位システム、および状態検出プログラム
CN108169774B (zh) 支持rtppp和rtk的多模gnss单频周跳探测与修复方法
JP5305416B2 (ja) 衛星航法システムにおける電離圏異常を検出する方法及びその装置。
JP2013019893A (ja) バイアスがかかった測定値に基づいた衛星航法システムの誤り検出
EP2562561B1 (en) Method and apparatus for receiving positioning signals based on pseudorange corrections
KR101574819B1 (ko) 비정상적 의사-거리 측정에 대해 무선 내비게이션 수신기 사용자를 보호하는 방법
EP1952174B1 (en) Methods and apparatus to detect and correct integrity failures in satellite positioning system receivers
US8742984B2 (en) Integrity method for differential corrections
CN106066485B (zh) 一种周跳探测方法及装置
KR20120026998A (ko) 의사 거리 측정치를 선택함으로써 위치 추정치를 수정하기 위한 방법
JP2008298443A (ja) マルチパス検出装置、測位装置、姿勢方位標定装置、マルチパス検出方法およびマルチパス検出プログラム
CN110879407B (zh) 一种基于完好性风险模型的卫星导航观测量新息检测方法
JP6593879B2 (ja) 衛星測位システム、測位端末、測位方法、及びプログラム
WO2023020056A1 (zh) 定位授时方法、装置、终端设备以及计算机可读存储介质
KR20120131093A (ko) Gps 수신기의 위치 계산 방법, 그 방법을 구현하기 위한 프로그램이 기록된 기록매체 및 gps 수신기
CN114859387A (zh) 静态站周跳探测的方法、装置、设备及计算机存储介质
US20190004181A1 (en) Method and apparatus applicable to positioning in nlos environment
KR101503001B1 (ko) 지상설비 안테나와 기저선 길이 예측값을 이용한 위성항법시스템의 고장 판단 시스템 및 그 방법
KR20190050554A (ko) 시스템간 시각차이 보상을 통한 다중위성군 gnss 측위 시스템 및 그 방법
KR101040053B1 (ko) 의사거리 측정잡음 감시기반 위성전파항법시스템 위성시계 고장검출 및 고장위성 식별 방법
KR101040054B1 (ko) 수신기 시계오차 감시기반 위성전파항법시스템 위성시계 고장검출 및 고장위성 식별 방법
US20200116873A1 (en) Positioning method and positioning terminal
CN110941002B (zh) 一种自适应抗差的序贯最小二乘精密单点定位方法
RU2475776C2 (ru) Способ управления навигационной системой и навигационная система
JP2014044056A (ja) 測位装置、測位方法および測位プログラム

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 14756490

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2015502766

Country of ref document: JP

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 14764476

Country of ref document: US

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 14756490

Country of ref document: EP

Kind code of ref document: A1