EP3017316B1 - Mitigation of multipath distortions for tdoa-based geolocation - Google Patents

Mitigation of multipath distortions for tdoa-based geolocation Download PDF

Info

Publication number
EP3017316B1
EP3017316B1 EP14819980.5A EP14819980A EP3017316B1 EP 3017316 B1 EP3017316 B1 EP 3017316B1 EP 14819980 A EP14819980 A EP 14819980A EP 3017316 B1 EP3017316 B1 EP 3017316B1
Authority
EP
European Patent Office
Prior art keywords
receivers
measurements
signal
tdoa
transmitter
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
EP14819980.5A
Other languages
German (de)
English (en)
French (fr)
Other versions
EP3017316A1 (en
EP3017316A4 (en
Inventor
Eyal BOKER
Avinoam LEVY
Rami Goldreich
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Elit Systems Bmd And Land Ew - Elisra Ltd
Original Assignee
Elbit Systems Bmd and Land Ew Elisra Ltd
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 Elbit Systems Bmd and Land Ew Elisra Ltd filed Critical Elbit Systems Bmd and Land Ew Elisra Ltd
Priority to PL14819980T priority Critical patent/PL3017316T3/pl
Publication of EP3017316A1 publication Critical patent/EP3017316A1/en
Publication of EP3017316A4 publication Critical patent/EP3017316A4/en
Application granted granted Critical
Publication of EP3017316B1 publication Critical patent/EP3017316B1/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
    • G01S5/06—Position of source determined by co-ordinating a plurality of position lines defined by path-difference measurements
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
    • G01S5/0205—Details
    • G01S5/0218—Multipath in signal reception
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
    • G01S5/0205—Details
    • G01S5/0221—Receivers
    • G01S5/02213—Receivers arranged in a network for determining the position of a transmitter
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
    • G01S5/0273—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves using multipath or indirect path propagation signals in position determination

Definitions

  • the disclosed technique generally relates to radiolocation in a multipath environment.
  • Geolocation techniques attempt to provide an accurate determination of the current geographic location of a particular object.
  • a common form is radiolocation, which utilizes the transmission and reception of radio frequency (RF) signals and measurements of the RF signal characteristics.
  • RF radio frequency
  • the location of a transmitting source can be determined using a plurality of receivers that obtain measurements relating to the intensity (signal strength) or the phase (delay) of the transmitted RF signal.
  • One approach uses the received signal strength (RSS) to estimate the distance to the transmitting source, assuming that the transmitted signal strength is known.
  • the transmission distance can be estimated based on the time of arrival (TOA) at the receivers if the transmission time and speed of propagation are known.
  • TOA time of arrival
  • the angle of arrival (AOA) of the transmitted signal which can be obtained using a directional antenna, may be combined with the distance estimate to provide the transmitter position.
  • a further approach is known as time difference of arrival (TDOA), which is a form of multilateration that measures the differences between the arrival times at multiple spatially separated receivers situated at known locations. For each pair of receivers where their respective locations and the TDOA between them are provided, the set of possible locations for the transmitter that satisfy the measurements will form a hyperbolic curve. By utilizing additional TDOA measurements for other pairs of receivers, multiple hyperbolas can be constructed, and their intersection provides an indication for the likely location of the transmitter.
  • TDOA time difference of arrival
  • ⁇ 12 and ⁇ 12 represent the relative attenuation and the delay between the two sensors, respectively.
  • the nature of the multipath which is characterized by the number of reflections, the reflections power and the reflections delays is dependent on the nature of the environment. Dense urban environments are characterized by a large number of reflections which may be stronger than the direct path signals. Even if the multipath reflection signal is relatively weak (e.g., less than 40dB), it can still have a significant influence on the accuracy of the final location estimate.
  • U.S. Patent No. 5,999,129 to Rose entitled “Multiplatform ambiguous phase circle and TDOA protection emitter location", is directed to the geolocation of a stationary RF signal emitter from two or more moving observer aircraft.
  • the observers receive signals from the emitter, and the ambiguous phase difference between the signals is measured at corresponding update intervals.
  • the observers perform pulse TOA measurements of the received signals over a predetermined clock interval, from which the TDOA of corresponding same-pulse emitter signals are calculated.
  • a series of circular lines of position (LOP) are estimated for each observer based on the ambiguous phase differences measured and associated integer values.
  • Hyperbolic LOPs are also computed based on the TDOA calculations.
  • the emitter location is determined based on the intersection of the hyperbolic LOPs and the circular LOPs.
  • U.S. Patent No. 7,132,981 to Roberts entitled “Method of locating object using phase differences among multiple frequency beacons transmitted from spaced apart transmitter sites", is directed to a technique for geolocating a mobile object in an environment where other locating systems (e.g. GPS systems) may not be expected to operate successfully.
  • At least three transmitter sites whose geolocations are fixed and known transmit dual frequency beacons that are readily received by a mobile receiver within the environment of interest.
  • a receiver located with the object processes the three sets of received signals by measuring phase differences among respective pairs of beacons.
  • the object receiver processes the phase differences to obtain TDOA or TOA information, from which the location the object can be derived. Residual phase errors and frequency offsets may be calibrated out at the receiver.
  • U.S. Patent No. 7,139,583 to Yamasaki et al entitled “Positioning system and method based on time difference of arrival", is directed to a positioning system that measures the time difference of arrival between a terminal station and a plurality of access points, by accumulating a plurality of reception timings for measurement signals which have passed through different propagation environments, and performing a statistic process with respect to the reception timing samples.
  • each of the access points repeatedly measures a signal transmitted from the station, while varying parameters such as the time, position, or frequency. For example the measurements are performed while slightly changing the position of an antenna at an access point or the position of a station, or while the changing the frequency channel in use.
  • Delay profiles are created from the measurement signals, and a given timing is determined on each delay profile as a reception timing sample for the signal.
  • the reception timing samples are combined, and one reception timing is determined for each access point as a reference timing.
  • the TDOA between the station and access points is calculated based on the reception timings, from which the station position coordinates are calculated.
  • U.S. Patent No. 7,911,385 to Heuser entitled “RF transmitter geolocation system and related methods", is directed to a system for geolocating an RF transmitter in the presence of multipath interference.
  • the system includes a plurality of RF receivers arranged in spaced relation, and a controller configured to generate a plurality of measurements associated with the RF transmitter.
  • the measurements may be at least one of: frequency difference of arrival measurements, time of arrival measurements, time difference of arrival measurements, frequency of arrival measurements, and angle of arrival measurements.
  • the controller computes a plurality of ambiguity functions based upon the measurements and due to the multipath interference, and projects the ambiguity functions onto a common geo-referenced grid.
  • the controller detects a peak on the common geo-referenced grid indicative of a geolocation of the RF transmitter.
  • the common geo-referenced grid may be generated based upon the position, orientation, and/or the relative movement of the RF receivers.
  • U.S. Patent No. 8,049,668 to Gines et al entitled “Method and system for locating signal emitters using iterative residual weighting”, discloses determining the location of a signal emitting device using at least three sensors separated and spaced apart from each other. Estimated location data for the device is determined for each sensor or sensor pair, based on which an estimated location is determined. The estimated location data may be based on the angle of arrival, time of arrival, time difference of arrival or the relative signal powers of the received signals. Residual values for the estimated location data is determined, based on a difference between the estimated device location and the estimated location data for a corresponding sensor or sensor pair. The residual values are converted into corresponding weights for the estimated location data for each sensor or unique sensor pair, which are weighted accordingly. The estimated device location is then updated using the weighted estimated location data.
  • U.S. Patent No. 8,077,089 to Parker entitled “Precision geolocation of moving or fixed transmitters using multiple observers", discloses the geolocation of a moving transmitter using a plurality of fixed receiving devices and a moving receiver device.
  • the fixed receivers obtain wavelength-scaled phase-difference measurements, providing a shape of the transmitter trajectory.
  • the phase difference between the moving receiver and at least one fixed receiver is measured to obtain a phase error residual.
  • An estimated starting point of the transmitter is moved to obtain a best-fit residual.
  • the Time of Arrival (TOA), Time Difference of Arrival (TDOA), or Frequency Difference of Arrival (FDOA) for the signal at each receiver is calculated, the clock offset between receivers is determined based on the TOA/TDOA/FDOA of the signal within each receiver, and the receiver clocks are calibrated based on the clock offset.
  • a differential measurement error may be calculated by comparing the expected difference between signal parameter values and the actual difference in parameter measurements, and corrections may be applied to the location system equipment to minimize errors in subsequent measurements and/or applied to the calculated parameter measurements.
  • a method for determining a location of a radio frequency (RF) transmitter in the presence of multipath interference includes the procedure of providing a plurality of spatially separated RF receivers at known locations in a moderate multipath environment in the vicinity of the transmitter.
  • the method further includes the procedure of detecting an active RF signal associated with the transmitter received by at least one of the receivers and directing the receivers to acquire measurements of the detected RF signal.
  • the method further includes the procedure of repeatedly repositioning the receivers along a circular trajectory at a sufficiently rapid rate to obtain a sequence of measurements of the RF signal by each of the receivers at different positions along the circular trajectory at a diversity of relative phase differences between the direct-path and the multipath reflections of the RF signal.
  • the method further includes the procedure of calculating time difference of arrival (TDOA) measurements between pairs of the receivers based on the obtained measurements.
  • the method further includes the procedure of averaging the TDOA measurements for each of the pairs of receivers to provide a respective updated TDOA measurement value.
  • the method further includes the procedure of determining a location of the transmitter based on the updated TDOA measurement values.
  • the range of the circular trajectory may be at least half the wavelength of the detected RF signal.
  • the receivers may be repeatedly repositioned along a circular trajectory by simultaneously repositioning a plurality of sub-receivers associated with the receiver along respective sections of the circular trajectory.
  • the receivers may be repeatedly repositioned along a circular trajectory by simultaneously repositioning a plurality of receivers associated with a plurality of RF transmitters operating in different spectral ranges.
  • a system for determining a location of an RF transmitter in the presence of multipath interference includes a plurality of RF receivers, and a processor.
  • the receivers are spatially separated at known locations in a moderate multipath environment in the vicinity of the transmitter.
  • the processor detects an active RF signal associated with the transmitter received by at least one of the receivers, and directs the receivers to acquire measurements of the detected RF signal.
  • the processor further calculates TDOA measurements between pairs of the receivers based on a sequence of measurements of the RF signal obtained by each of the receivers at different positions along a circular trajectory by repeatedly repositioning the receivers along the circular trajectory at a sufficiently rapid rate to obtain measurements at a diversity of relative phase differences between the direct-path and the multipath reflections of the RF signal.
  • the processor further averages the TDOA measurements for each of pair of receivers to provide a respective updated TDOA measurement value.
  • the processor further determines a location of the transmitter based on the updated TDOA measurement values.
  • the range of the circular trajectory may be at least half the wavelength of the detected RF signal.
  • the system may further include repositioning means, for repeatedly repositioning at least one of the receivers to different positions along the circular trajectory.
  • the repositioning means may simultaneously reposition a plurality of sub-receivers associated with the receiver along respective sections of the circular trajectory.
  • the repositioning means may simultaneously reposition a plurality of receivers associated with a plurality of RF transmitters operating in different spectral ranges.
  • the disclosed technique overcomes the disadvantages of the prior art by providing a system and method for determining the location of a stationary transmitter with high accuracy via a TDOA-based radiolocation technique by mitigating the distortions resulting from multipath signal propagation.
  • the disclosed technique involves repositioning the receivers in a controlled pattern while receiving the transmitted signal, such as by moving the receivers along a circular trajectory, and collecting a series of different readings at different positions along the trajectory.
  • the subsequent TDOA measurements for the receiver pairs are then averaged, substantially cancelling out the multipath distortions in the signal readings, which tend to be distributed substantially uniformly for varying receiver positions.
  • the transmitter location is then calculated based on the averaged TDOA measurements.
  • location refers to any parameter relating to the global or relative position of a particular object with respect to any relevant reference coordinate system (e.g., latitude and longitude coordinates).
  • Figure 1 is a graph, generally referenced 110, illustrating the effect of multipath propagation on TDOA measurements, by means of a simulation.
  • the simulation modelled two FM modulated signals at a bandwidth of 25kHz for a pair of receivers located at identical ranges from the transmitter, hence the expected TDOA value is 0.
  • an additional shifted and attenuated signal replica (representing a single reflection) was added to the first receiver.
  • the reflection was delayed by 20 ⁇ s and attenuated by -10dB with respect to the direct-path signal.
  • the reflection was also multiplied by different phases ranging between 0 - 2 ⁇ .
  • Graph 110 depicts the TDOA measurement error (ns) as a function of the reflection phase (radians) for the receiver pair.
  • the text box in graph 110 presents (from left to right): the average TDOA error, the TDOA error standard deviation, and the TDOA root-mean-square error (RMSE).
  • TDOA error refers to the deviation of a certain measurement value from the expected value
  • TDOA bias refers to the deviation of the final averaged value (the average of all the measurements) from the expected value, also equivalent to the residual error after averaging.
  • the multipath effect is strongly dependent on the reflection phase, and introduces a maximum TDOA error of approximately 8 ⁇ s (peak value), which corresponds to a position error greater than 2km for the resultant transmitter position estimate.
  • the TDOA error as a function of the reflection phase is distributed non-uniformly about the expected value.
  • the average value is still biased relative to the expected value by 187ns (corresponding to the TDOA bias value), but is closer to the expected value than most of the individual measurement values taken independently.
  • FIG. 2 is a graph, generally referenced 120, illustrating the TDOA measurement bias dependence upon the delay and attenuation of a multipath reflection, by means of a simulation.
  • Graph 120 shows the TDOA root-mean-squared error (RMSE) introduced by one reflection in one receiver as a function of the reflection delay for various reflection coefficients.
  • the reflection coefficient represents the power loss (attenuation) of the reflection relative to the direct-path signal.
  • the reflection power relative to the direct-path power is a function of the relative path length and the reflection coefficient.
  • the TDOA RMSE increases with the reflection delay at lower delay ranges, but decreases with the reflection delay at higher delay ranges.
  • System 130 includes a radio frequency (RF) transmitter 132, a plurality of RF receivers 134, 136, 138, and a processor 140.
  • Processor 140 is coupled with receivers 134, 136, 138.
  • Receivers 134, 136, 138 are situated at separate locations but are each disposed within sufficient range to receive an RF signal transmission from transmitter 132.
  • Transmitter 132 may include standard components and/or circuitry associated with the transmission of a radio frequency signal, such as a signal generator, a modulator, an antenna, a power supply, and the like.
  • each of receivers 134, 136, 138 may include standard components and/or circuitry associated with the reception of a radio frequency signal, such as filters, amplifiers, an antenna, a power supply, and the like.
  • the disclosed technique is generally applicable to RF communication signals propagating through a moderate multipath environment in which there are relatively few multipath reflections and which entails line-of-sight propagation between the transmitter and the receivers (i.e., such that the direct-path signal is at least stronger than the multipath reflections).
  • the disclosed technique is applicable to RF signals having a bandwidth for which TDOA-based calculations would be feasible and would provide reasonably accurate results, such as bandwidths on the order of tens of kHz (higher bandwidths would generally require using TOA calculations).
  • the RF signal may be any frequency within the RF spectral range (e.g., VHF, UHF), and may be modulated.
  • transmitter 132 and receivers 134, 136, 138 are disposed at or near ground level, where multipath distortions may be particularly likely to occur. Nevertheless, the transmitter and receivers may alternatively be high above ground level, provided that they are in any environment subject to moderate multipath interference with line-of-sight propagation.
  • Transmitter 132 is situated at an unknown location and is substantially stationary, such that it does not shift its position (beyond a marginal amount) throughout the duration of the location determination process.
  • Receivers 134, 136, 138 are situated at separate known locations, and are movable to various positions for receiving a transmitted RF signal at a plurality of spatially separate locations, as will be elaborated upon further below.
  • Processor 140 provides synchronization between receivers 134, 136, 138, to enable the subsequent analysis of the time difference of arrival data. Processor 140 may also control the activation and operation of receivers 134, 136, 138, and performs signal processing. It is noted that any of the functionality of processor 140 may be implemented via a plurality of processors and/or may be partially or fully distributed among any of the other components of system 130. For example, at least one of receivers 134, 136, 138 may be integrated with a dedicated microprocessor, which may perform at least a portion of the signal processing associated with the respective component.
  • processor 140 instructs receivers 134, 136, 138 to collect data in a selected spectral range.
  • Processor 140 analyzes the received data to detect active RF signals that may be associated with transmitter 132. When a suitable active RF signal is detected as having been received by at least one receiver, processor 140 instructs each of receivers 134, 136, 138 to acquire measurements of the detected signal.
  • the detected signal is equivalent to the RF signal 142 transmitted by transmitter 132 convolved with the transmission function between transmitter 132 and the respective receiver 134, 136, 138 (the transmission function which may include direct-path and/or reflections of the RF signal 142).
  • Each of receivers 134, 136, 138 performs repeated measurements of RF signal 142 while shifting positions.
  • the receiver follows a predetermined path while acquiring a series of measurements of RF signal 142 at different points along the path.
  • receiver 134 at location R1 is moved along trajectory 144, and acquires measurements at each of trajectory points 144A, 144B, 144C, 144D and 144E.
  • receiver 136 at location R2 is moved along trajectory 146, and acquires measurements at each of trajectory points 146A, 146B, 146C, 146D and 146E.
  • receiver 138 at location R3 is moved along trajectory 148 and acquires measurements at each of trajectory points 148A, 148B, 148C, 148D and 148E.
  • the range of the movement pattern traversed by the receiver is determined as a function of the wavelength of the measurement signal 142, to ensure that measurements are obtained at various phase differences distributed between the entire range of 0-2 ⁇ .
  • the radius of trajectories 144, 146, 148 should be at least half the wavelength of RF signal 142 (or correspondingly, the minimum diameter is equal to the signal wavelength).
  • the movement pattern is along the plane at which the RF signals (direct-path and reflections) propagate from transmitter 132, to ensure that relative phase differences are generated between the direct path and reflections.
  • trajectories 144, 146 and 148 are depicted as substantially circular, which provides a uniform span of reflection delays regardless of the constellation of the transmitter and receiver. Other trajectory shapes or patterns may be less beneficial, but may be applied provided that a diversity of phase differences between the direct-path and reflection signals are obtained.
  • the receiver movement may be implemented using any suitable mechanism or technique for shifting the position at which the signal reception is performed by the respective receiver.
  • the receiver may be repositioned using an automatic mechanical device, such as a rotating device that moves the receiver along a circular trajectory.
  • a receiver may be coupled with a mobile platform, such as being mounted onto a vehicle or onto a robotic trolley, which can be maneuvered along a particular path to transport the coupled receiver to different positions along the path.
  • a receiver may be held or carried by a person, who can manually transport the receiver to different positions.
  • the rate at which the receiver is moved along its trajectory should be at least rapid enough to obtain sufficient coverage of relative phase differences between the direct-path and reflection signals.
  • the receiver movement may involve repeatedly following the same trajectory over multiple cycles, for example by shifting the receiver multiple times around a circular pattern (such as circular trajectories 144, 146 and 148). Repeatedly following the same trajectory path in a cyclic manner ensures that a sufficient number of measurements representing phase differences distributed between the entire range of 0-2 ⁇ are attained.
  • a rotating mechanism may be configured to repeatedly reposition the receiver along a circular trajectory at a rate of 60 rpm (i.e., corresponding to one full rotation each second).
  • the measurement data at includes at least the time difference of arrival (TDOA) of RF signal 142 between pairs of receivers 134, 136, 138 at each of the trajectory positions.
  • processor 140 calculates a series of TDOA measurements between receivers 134 and 136, between receivers 134 and 138, and between receivers 136 and 138, for each of the different positions of the receivers in the respective pairing. Subsequently, processor 140 calculates an average value of all the TDOA measurements for each receiver pair. Processor 140 then calculates the location of transmitter 132 based on the average TDOA values.
  • system 130 may include more than three receivers and localization can be performed based on TDOA data obtained for only three of the receiver pairings, for all possible receiver pairings, or for any number in between (the position estimation is more accurate as the number of receivers is increased).
  • Three receivers is the minimum number required to allow for localization, although it is possible to achieve the effect of a certain number of receivers using a smaller number of actual receivers (e.g., using two receivers, at least one of which is shifted to a substantially different location over the course of the localization, to effectively represent the position of a "third receiver").
  • receivers While in theory, if it is known in advance which receivers are influenced by multipath interference, then can reposition only the relevant receivers (rather than all of them) in order to obtain an improved location estimate. In practice, however, such information is not available, and so the general approach is to reposition all of the receivers.
  • a plurality of receivers are prearranged at different fixed positions along the trajectory, or alternatively, may be selectively deployed to fixed positions along the trajectory, e.g., by activating a switch that deploys all the receivers. In this manner, the need for physically moving or repositioning the receiver is obviated.
  • the sequence of measurements obtained at different positions along a trajectory by a receiver may be obtained in practice using multiple receivers (or receiving antennas), which are simultaneously repositioned along respective sections of the entire trajectory.
  • the measurements of receiver 134 along trajectory 144 may be obtained using a pair of "sub-receivers" 134A and 134B, where a first sub-receiver 134A is successively repositioned around one half of the circumference of circular trajectory 144 (e.g., passing through at least trajectory points 144A, 144B and 144C), while concurrently, the second sub-receiver 134B is successively repositioned around the other half of the circumference of circular trajectory 144 (e.g., passing through at least trajectory points 144D, 144E and 144A).
  • This pattern may be repeated multiple times (i.e., in a cyclic matter), as necessary.
  • four sub-receivers associated with receiver 134 may be concurrently repositioned simultaneously along respective quadrants of circular trajectory 144.
  • the repositioning mechanism may achieve a greater number of total measurements for a given duration.
  • a rotating device would obtain twice the number of measurements in the same period of time using two sub-receivers moving concurrently over half sectors of a circular trajectory, as compared to one receiver moving along the entire trajectory at the same rate.
  • a rotating device operating at 60 rpm with two sub-receivers would obtain twice the total number of measurements per second than when using a single receiver.
  • a single repositioning mechanism may be configured to reposition multiple receivers operative to detect different transmitters operating in different spectral domains.
  • a first set of RF receivers 134A, 136A, 138A are directed to determine the geolocation of a first RF transmitter 132A operative in a first RF spectral range domain (e.g., in the VHF range), while a second set of RF receivers 134B, 136B, 138B are directed to determine the geolocation of a second RF transmitter 132B operative in a second RF spectral range domain (e.g., in the UHF range), where both sets of transmitters/receivers are in the same vicinity.
  • a rotating device may be configured to rotate at least one receiver belonging to each set, together or separately, along the same or different circular trajectories.
  • the rotating device may rotate both receivers 134A and 134B along circular trajectory 144, such that receiver 134A obtains a sequence of RF measurements respective of transmitter 132A along different points of trajectory 144, while receiver 134B obtains a sequence of RF measurements respective of transmitter 132B along different points of trajectory 144.
  • FIG. 4 is a graph, generally referenced 160, illustrating the level of TDOA bias when the receivers are repositioned while acquiring measurements in accordance with the disclosed technique, based on the results of an experiment.
  • the experiment for graph 160 was performed using three receivers situated at the vertices of an approximately equilateral triangle with side lengths of -5km, and a transmitting RF signal with a 25kHz bandwidth at a frequency of 79MHz located inside the triangle.
  • Each of the receivers was repositioned along the perimeter of a 3m radius circular trajectory.
  • TDOA measurements were obtained for segments of 400ms every 200ms.
  • a segment represents a time instant in which the receivers are considered to be substantially stationary (it is assumed that the receiver does not move much over 400ms).
  • each segment represents a certain realization of the receiver position (which corresponds to a certain realization of relative phase).
  • Graph 160 depicts the TDOA measurement (in units of ns) as a function of the segment for the three receiver pairings.
  • the estimated TDOA readings are represented by the triangular, circular, and dotted symbols, respectively, and the corresponding expected TDOA values are represented by the dashed lines.
  • the text box in graph 160 provides statistical properties of the TDOA error for the TDOA measurements associated with each combination of receiver pairs, where the first (leftmost) column indicates the average TDOA error, the second column indicates the standard deviation of the TDOA error, and the third column indicates the TDOA RMSE.
  • Graph 160 demonstrates that the multipath effect is dominant for all of the receivers, as expressed by the high TDOA RMSE values of -600ns and maximum TDOA bias above 1 ⁇ s.
  • This level of TDOA bias may result in a localization error (i.e., the error in the resultant transmitter position estimate using the associated TDOA values) of above 500m.
  • the segments in which the TDOA error is significant correspond to receiver positions realizing phase differences that result in a high TDOA bias.
  • the TDOA measurements for each of the receiver pairs are distributed around the expected value. By averaging the TDOA measurements, the TDOA bias for the corresponding receiver pairs was reduced from approximately 1 ⁇ s (worst case) to approximately 250ns.
  • the resulting localization error is reduced from approximately 500m to approximately 38m.
  • FIG. 5 is a block diagram of a method for determining a location of an RF transmitter in the presence of multipath interference, operative in accordance with an embodiment of the disclosed technique.
  • procedure 182 a plurality of spatially separated RF receivers are provided at known locations in a moderate multipath environment in the vicinity of a stationary RF transmitter.
  • receivers 134, 136 and 138 are disposed at separate respective locations R1, R2 and R3, and are all within reception range of transmitter 132.
  • the environment is subject to multipath reflections and there is line-of-sight propagation between transmitter 132 and receivers 134, 136, 138
  • an active RF signal associated with the transmitter received by at least one receiver is detected, and the receivers are directed to acquire measurements of the detected RF signal.
  • processor 140 instructs receivers 134, 136, 138 to collect data in a selected spectral range, and analyzes the received data to detect active RF signals that may be associated with transmitter 132.
  • processor 140 instructs each of receivers 134, 136, 138 to acquire measurements of the detected signal, which may include the direct-path and/or reflections of an RF signal 142 being transmitted by transmitter 132.
  • a sequence of measurements of the RF signal is obtained at different positions along a trajectory be each of the receivers.
  • each of receivers 134, 136 and 138 are moved to various positions along a trajectory while obtaining readings of RF signal 142.
  • receiver 134 is sequentially shifted to positions 144A, 144B, 144C, 144D and 144E along circular trajectory 144, and obtains a respective reading at each position.
  • receiver 136 is sequentially shifted to various positions along trajectory 146
  • receiver 138 is sequentially shifted to various positions along trajectory 148.
  • the sequence of measurements may be obtained by multiple receivers (or receiving antennas) that are arranged at fixed positions along the trajectory, or that are selectively deployed to fixed positions along the trajectory.
  • TDOA measurements between pairs of receivers are calculated based on the obtained measurements.
  • processor 140 calculates TDOA measurements for different pairs of receivers 134, 136, 138, based on the obtained series of readings at different positions at each receiver.
  • TDOA measurements are obtained between receivers 134 and 136, between receivers 134 and 138, and between receivers 136 and 138, for each of the different positions of the receivers in the respective pairing.
  • processor 140 calculates an average value of all the TDOA measurements for each receiver pair. For example, processor 140 averages all the TDOA measurements between receivers 134 and 136 to obtain a final averaged TDOA value between receivers 134 and 136. Similarly, processor 140 obtains an averaged TDOA value between receivers 134 and 138, and between receivers 136 and 138.
  • a location of the transmitter is determined based on the updated TDOA measurement values.
  • processor 140 calculates the location of transmitter 132 based on the average TDOA values calculated for each receiver pair.
  • processor 140 uses the updated TDOA between each pair among receivers 134, 136, 138, to obtain a position estimate of transmitter 132.
  • the disclosed technique is applicable to geolocation for any purpose, and can be employed in a wide variety of applications, including both civilian and military applications.
  • the disclosed technique may be utilized in communications infrastructures, such as to locate potential pirate transmissions that are interfering with sanctioned broadcasts and other forms of spectrum management.
  • Another example is in conjunction with cellular networks, such as to provide accurate assessment of a mobile device, e.g., in addition to or instead of standard global positioning systems (GPS).
  • GPS global positioning systems
  • a plurality of cellular antennas may be configured to rotate along different planar trajectories in order to obtain receptions of a transmitted signal for which the direction of arrival is not known in advance.
  • a further potential application of the disclosed technique is for localization of missing or stolen vehicles.
EP14819980.5A 2013-07-01 2014-06-30 Mitigation of multipath distortions for tdoa-based geolocation Active EP3017316B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
PL14819980T PL3017316T3 (pl) 2013-07-01 2014-06-30 Łagodzenie zakłóceń wielodrogowych dla geolokacji opartej na TDOA

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
IL227285A IL227285A0 (en) 2013-07-01 2013-07-01 Multipath interference reduction for geo-location based on time-of-arrival
PCT/IL2014/050584 WO2015001552A1 (en) 2013-07-01 2014-06-30 Mitigation of multipath distortions for tdoa-based geolocation

Publications (3)

Publication Number Publication Date
EP3017316A1 EP3017316A1 (en) 2016-05-11
EP3017316A4 EP3017316A4 (en) 2016-12-21
EP3017316B1 true EP3017316B1 (en) 2019-01-02

Family

ID=49784201

Family Applications (1)

Application Number Title Priority Date Filing Date
EP14819980.5A Active EP3017316B1 (en) 2013-07-01 2014-06-30 Mitigation of multipath distortions for tdoa-based geolocation

Country Status (9)

Country Link
US (1) US9841489B2 (ko)
EP (1) EP3017316B1 (ko)
KR (1) KR101834093B1 (ko)
DK (1) DK3017316T3 (ko)
ES (1) ES2717947T3 (ko)
IL (1) IL227285A0 (ko)
PL (1) PL3017316T3 (ko)
SG (1) SG11201510708WA (ko)
WO (1) WO2015001552A1 (ko)

Families Citing this family (41)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11175375B2 (en) 2010-11-12 2021-11-16 Position Imaging, Inc. Position tracking system and method using radio signals and inertial sensing
US10416276B2 (en) 2010-11-12 2019-09-17 Position Imaging, Inc. Position tracking system and method using radio signals and inertial sensing
US9945940B2 (en) 2011-11-10 2018-04-17 Position Imaging, Inc. Systems and methods of wireless position tracking
US9933509B2 (en) * 2011-11-10 2018-04-03 Position Imaging, Inc. System for tracking an object using pulsed frequency hopping
US10269182B2 (en) 2012-06-14 2019-04-23 Position Imaging, Inc. RF tracking with active sensory feedback
US10180490B1 (en) 2012-08-24 2019-01-15 Position Imaging, Inc. Radio frequency communication system
US10234539B2 (en) 2012-12-15 2019-03-19 Position Imaging, Inc. Cycling reference multiplexing receiver system
US9482741B1 (en) 2013-01-18 2016-11-01 Position Imaging, Inc. System and method of locating a radio frequency (RF) tracking device using a calibration routine
US10856108B2 (en) 2013-01-18 2020-12-01 Position Imaging, Inc. System and method of locating a radio frequency (RF) tracking device using a calibration routine
US10634761B2 (en) 2013-12-13 2020-04-28 Position Imaging, Inc. Tracking system with mobile reader
US9497728B2 (en) 2014-01-17 2016-11-15 Position Imaging, Inc. Wireless relay station for radio frequency-based tracking system
US10200819B2 (en) 2014-02-06 2019-02-05 Position Imaging, Inc. Virtual reality and augmented reality functionality for mobile devices
US11132004B2 (en) 2015-02-13 2021-09-28 Position Imaging, Inc. Spatial diveristy for relative position tracking
US10642560B2 (en) 2015-02-13 2020-05-05 Position Imaging, Inc. Accurate geographic tracking of mobile devices
US10324474B2 (en) 2015-02-13 2019-06-18 Position Imaging, Inc. Spatial diversity for relative position tracking
US11501244B1 (en) 2015-04-06 2022-11-15 Position Imaging, Inc. Package tracking systems and methods
US11416805B1 (en) 2015-04-06 2022-08-16 Position Imaging, Inc. Light-based guidance for package tracking systems
US10853757B1 (en) 2015-04-06 2020-12-01 Position Imaging, Inc. Video for real-time confirmation in package tracking systems
US10148918B1 (en) 2015-04-06 2018-12-04 Position Imaging, Inc. Modular shelving systems for package tracking
US10200886B2 (en) * 2015-09-14 2019-02-05 Red Point Positioning Corporation Method to estimate and compensate for NLOS bias in time difference of arrival estimate
EP3637126A1 (en) 2016-01-05 2020-04-15 Elta Systems Ltd. Method of locating a transmitting source in multipath environment and system thereof
EP3193187A1 (en) * 2016-01-13 2017-07-19 ST SPORTSERVICE GmbH Method for calibrating a local positioning system based on time-difference-of-arrival measurements
US10444323B2 (en) 2016-03-08 2019-10-15 Position Imaging, Inc. Expandable, decentralized position tracking systems and methods
CN105929365B (zh) * 2016-04-27 2018-07-24 成都精位科技有限公司 一种uwb高精度定位系统及定位方法
KR101881868B1 (ko) * 2016-07-15 2018-07-25 국방과학연구소 펄스 신호 수집정보의 확률적 선택에 의한 위치탐지방법
US11436553B2 (en) 2016-09-08 2022-09-06 Position Imaging, Inc. System and method of object tracking using weight confirmation
JP6551615B2 (ja) * 2016-11-29 2019-07-31 株式会社村田製作所 位置推定装置及び位置推定方法
US10634503B2 (en) 2016-12-12 2020-04-28 Position Imaging, Inc. System and method of personalized navigation inside a business enterprise
US10634506B2 (en) 2016-12-12 2020-04-28 Position Imaging, Inc. System and method of personalized navigation inside a business enterprise
US10455364B2 (en) 2016-12-12 2019-10-22 Position Imaging, Inc. System and method of personalized navigation inside a business enterprise
US11120392B2 (en) 2017-01-06 2021-09-14 Position Imaging, Inc. System and method of calibrating a directional light source relative to a camera's field of view
IL256679B2 (en) * 2017-12-31 2023-02-01 Elta Systems Ltd Systems and methods for locating a signal source
US11802936B2 (en) * 2018-01-15 2023-10-31 Sony Network Communications Europe B.V. Determining a position of a device with respect to another device
US20210263128A1 (en) * 2018-09-14 2021-08-26 Telefonaktiebolaget Lm Ericsson (Publ) Network locationing rf planner
CA3111595A1 (en) 2018-09-21 2020-03-26 Position Imaging, Inc. Machine-learning-assisted self-improving object-identification system and method
CN109597026B (zh) * 2018-11-30 2022-03-18 四川九洲电器集团有限责任公司 一种多点定位综合仿真评估方法和系统
WO2020146861A1 (en) 2019-01-11 2020-07-16 Position Imaging, Inc. Computer-vision-based object tracking and guidance module
DE102020203238B4 (de) 2020-03-13 2023-02-16 Forschungszentrum Jülich GmbH Verfahren und System zur Bestimmung eines Mehrwegeeinflusses, Verfahren zur Positionsbestimmung eines Objekts und Vorrichtung zur Datenverarbeitung
CN111781556B (zh) * 2020-05-18 2023-08-18 中国民用航空总局第二研究所 基于目标信号tdoa滤波的分布式时差定位方法及装置
KR102384525B1 (ko) 2020-07-09 2022-04-08 엘지전자 주식회사 이동로봇의 제어방법
CN114615619B (zh) * 2022-02-25 2023-08-04 厦门大学 基于信号反射点的室内定位方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020196186A1 (en) * 2001-06-25 2002-12-26 Harris Corporation Method and system for calibrating wireless location systems

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5437055A (en) * 1993-06-03 1995-07-25 Qualcomm Incorporated Antenna system for multipath diversity in an indoor microcellular communication system
US6108555A (en) * 1996-05-17 2000-08-22 Ksi, Inc. Enchanced time difference localization system
US9288623B2 (en) * 2005-12-15 2016-03-15 Invisitrack, Inc. Multi-path mitigation in rangefinding and tracking objects using reduced attenuation RF technology
US7872583B1 (en) * 2005-12-15 2011-01-18 Invisitrack, Inc. Methods and system for multi-path mitigation in tracking objects using reduced attenuation RF technology
US8077089B2 (en) 2006-08-15 2011-12-13 Rincon Research Corporation Precision geolocation of moving or fixed transmitters using multiple observers
US7453400B2 (en) * 2007-02-02 2008-11-18 Bae Systems Information And Electronic Systems Integration Inc. Multiplatform TDOA correlation interferometer geolocation
US7911385B2 (en) 2009-02-27 2011-03-22 Harris Corporation RF transmitter geolocation system and related methods
CN101839971B (zh) 2009-03-17 2013-09-18 深圳市华颖锐兴科技有限公司 一种高精度gnss接收机
US8174444B2 (en) * 2009-09-26 2012-05-08 Rincon Research Corporation Method of correlating known image data of moving transmitters with measured radio signals
WO2011160698A1 (en) * 2010-06-25 2011-12-29 Innovationszentrum für Telekommunikationstechnik GmbH IZT System and method for determining a position
US8077091B1 (en) * 2010-07-06 2011-12-13 Intelligent Sciences, Ltd. System and method for determining a position of a mobile device within a surveillance volume in the presence of multipath interference

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020196186A1 (en) * 2001-06-25 2002-12-26 Harris Corporation Method and system for calibrating wireless location systems

Also Published As

Publication number Publication date
ES2717947T3 (es) 2019-06-26
EP3017316A1 (en) 2016-05-11
WO2015001552A1 (en) 2015-01-08
KR101834093B1 (ko) 2018-03-02
KR20160026991A (ko) 2016-03-09
PL3017316T3 (pl) 2019-06-28
US20160370453A1 (en) 2016-12-22
SG11201510708WA (en) 2016-01-28
EP3017316A4 (en) 2016-12-21
DK3017316T3 (en) 2019-04-15
US9841489B2 (en) 2017-12-12
IL227285A0 (en) 2013-12-31

Similar Documents

Publication Publication Date Title
EP3017316B1 (en) Mitigation of multipath distortions for tdoa-based geolocation
US9002378B2 (en) Method and system for estimation of mobile station velocity in a cellular system based on geographical data
US10557919B2 (en) Observed time difference of arrival angle of arrival discriminator
US7911385B2 (en) RF transmitter geolocation system and related methods
Porretta et al. Location, location, location
Gentner et al. Particle filter based positioning with 3GPP-LTE in indoor environments
US20030220117A1 (en) Location systems and methods
KR20050098901A (ko) 다수의 송신기들을 구비한 무선 시스템에서의 모바일 무선수신기 위치 찾기 시스템 및 방법
US11092665B2 (en) High frequency geo-location methods and systems
US20150241545A1 (en) Single Platform Doppler Geolocation
Dubois et al. GSM position tracking using a Kalman filter
Schmitz et al. TDOA fingerprinting for localization in non-line-of-sight and multipath environments
Wang et al. ToA and TDoA error models for NLoS propagation based on outdoor to indoor channel measurement
Garcia-Molina et al. Snapshot localisation of multiple jammers based on receivers of opportunity
Moschevikin et al. The impact of nlos components in time-of-flight networks for indoor positioning systems
Jain et al. HF source geolocation using an operational TDoA receiver network: Experimental results
Steffes et al. Direct position determination for TDOA-based single sensor localization
Faragher et al. Measurements of the effects of multipath interference on timing accuracy in a cellular radio positioning system
US20050176443A1 (en) Cellular phone geolocation system
Steffes et al. Direct single sensor TDOA localization using signal structure information
Weber et al. 3D RF emitter location estimation
Gabay et al. Radio frequency Echo mapping with cellular devices
Li et al. A Belief Propagation Algorithm for Multipath-based SLAM with Multiple Map Features: A mmWave MIMO Application
Abbas et al. Guessing the coordinates of targets in indoor environments using hyperbolic technique
Schmickmann et al. Bring Your Own Positioning System: An Infrastructure-free and Omnidirectional UWB-based Localization Approach

Legal Events

Date Code Title Description
PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

17P Request for examination filed

Effective date: 20160106

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

AX Request for extension of the european patent

Extension state: BA ME

DAX Request for extension of the european patent (deleted)
A4 Supplementary search report drawn up and despatched

Effective date: 20161118

RIC1 Information provided on ipc code assigned before grant

Ipc: G01S 5/06 20060101AFI20161114BHEP

Ipc: G01S 1/20 20060101ALI20161114BHEP

Ipc: G01S 3/46 20060101ALI20161114BHEP

Ipc: G01S 5/22 20060101ALI20161114BHEP

Ipc: G01S 5/02 20060101ALI20161114BHEP

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: EXAMINATION IS IN PROGRESS

17Q First examination report despatched

Effective date: 20170613

REG Reference to a national code

Ref country code: DE

Ref legal event code: R079

Ref document number: 602014039262

Country of ref document: DE

Free format text: PREVIOUS MAIN CLASS: G01S0005060000

Ipc: G01S0005020000

RIC1 Information provided on ipc code assigned before grant

Ipc: G01S 3/46 20060101ALI20180425BHEP

Ipc: G01S 5/06 20060101ALI20180425BHEP

Ipc: G01S 5/22 20060101ALI20180425BHEP

Ipc: G01S 5/02 20100101AFI20180425BHEP

Ipc: G01S 1/20 20060101ALI20180425BHEP

GRAP Despatch of communication of intention to grant a patent

Free format text: ORIGINAL CODE: EPIDOSNIGR1

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: GRANT OF PATENT IS INTENDED

INTG Intention to grant announced

Effective date: 20180723

RAP1 Party data changed (applicant data changed or rights of an application transferred)

Owner name: ELIT SYSTEMS BMD AND LAND EW - ELISRA LTD.

GRAS Grant fee paid

Free format text: ORIGINAL CODE: EPIDOSNIGR3

GRAA (expected) grant

Free format text: ORIGINAL CODE: 0009210

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE PATENT HAS BEEN GRANTED

AK Designated contracting states

Kind code of ref document: B1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

REG Reference to a national code

Ref country code: GB

Ref legal event code: FG4D

REG Reference to a national code

Ref country code: CH

Ref legal event code: EP

Ref country code: AT

Ref legal event code: REF

Ref document number: 1085093

Country of ref document: AT

Kind code of ref document: T

Effective date: 20190115

REG Reference to a national code

Ref country code: IE

Ref legal event code: FG4D

REG Reference to a national code

Ref country code: DE

Ref legal event code: R096

Ref document number: 602014039262

Country of ref document: DE

REG Reference to a national code

Ref country code: SE

Ref legal event code: TRGR

REG Reference to a national code

Ref country code: NL

Ref legal event code: FP

REG Reference to a national code

Ref country code: DK

Ref legal event code: T3

Effective date: 20190408

REG Reference to a national code

Ref country code: LT

Ref legal event code: MG4D

REG Reference to a national code

Ref country code: AT

Ref legal event code: MK05

Ref document number: 1085093

Country of ref document: AT

Kind code of ref document: T

Effective date: 20190102

REG Reference to a national code

Ref country code: ES

Ref legal event code: FG2A

Ref document number: 2717947

Country of ref document: ES

Kind code of ref document: T3

Effective date: 20190626

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: PT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190502

Ref country code: LT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190102

Ref country code: NO

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190402

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: BG

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190402

Ref country code: GR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190403

Ref country code: RS

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190102

Ref country code: IS

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190502

Ref country code: LV

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190102

Ref country code: HR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190102

REG Reference to a national code

Ref country code: DE

Ref legal event code: R097

Ref document number: 602014039262

Country of ref document: DE

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: IT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190102

Ref country code: RO

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190102

Ref country code: SK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190102

Ref country code: CZ

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190102

Ref country code: AL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190102

Ref country code: AT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190102

Ref country code: EE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190102

PLBE No opposition filed within time limit

Free format text: ORIGINAL CODE: 0009261

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: SM

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190102

26N No opposition filed

Effective date: 20191003

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: MC

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190102

REG Reference to a national code

Ref country code: CH

Ref legal event code: PL

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: SI

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190102

REG Reference to a national code

Ref country code: BE

Ref legal event code: MM

Effective date: 20190630

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: TR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190102

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: IE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20190630

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: LI

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20190630

Ref country code: BE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20190630

Ref country code: LU

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20190630

Ref country code: CH

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20190630

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: CY

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190102

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: HU

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT; INVALID AB INITIO

Effective date: 20140630

Ref country code: MT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190102

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: MK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190102

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: SE

Payment date: 20230314

Year of fee payment: 10

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: NL

Payment date: 20230515

Year of fee payment: 10

Ref country code: FR

Payment date: 20230510

Year of fee payment: 10

Ref country code: DK

Payment date: 20230613

Year of fee payment: 10

Ref country code: DE

Payment date: 20230502

Year of fee payment: 10

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: PL

Payment date: 20230403

Year of fee payment: 10

Ref country code: FI

Payment date: 20230615

Year of fee payment: 10

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: GB

Payment date: 20230511

Year of fee payment: 10

Ref country code: ES

Payment date: 20230707

Year of fee payment: 10