CN104237920B - 高反射环境中的交通工具定位 - Google Patents
高反射环境中的交通工具定位 Download PDFInfo
- Publication number
- CN104237920B CN104237920B CN201410279249.0A CN201410279249A CN104237920B CN 104237920 B CN104237920 B CN 104237920B CN 201410279249 A CN201410279249 A CN 201410279249A CN 104237920 B CN104237920 B CN 104237920B
- Authority
- CN
- China
- Prior art keywords
- main body
- signal
- speed
- transmitter
- estimation
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/21—Interference related issues ; Issues related to cross-correlation, spoofing or other methods of denial of service
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/14—Receivers specially adapted for specific applications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/22—Multipath-related issues
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/26—Acquisition or tracking or demodulation of signals transmitted by the system involving a sensor measurement for aiding acquisition or tracking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining 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/42—Determining position
- G01S19/421—Determining position by combining or switching between position solutions or signals derived from different satellite radio beacon positioning systems; by combining or switching between position solutions or signals derived from different modes of operation in a single system
- G01S19/426—Determining position by combining or switching between position solutions or signals derived from different satellite radio beacon positioning systems; by combining or switching between position solutions or signals derived from different modes of operation in a single system by combining or switching between position solutions or signals derived from different modes of operation in a single system
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining 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/42—Determining position
- G01S19/48—Determining position by combining or switching between position solutions derived from the satellite radio beacon positioning system and position solutions derived from a further system
- G01S19/49—Determining position by combining or switching between position solutions derived from the satellite radio beacon positioning system and position solutions derived from a further system whereby the further system is an inertial position system, e.g. loosely-coupled
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining 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/52—Determining velocity
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Navigation (AREA)
Abstract
公开了一种高反射环境中的交通工具定位。一种定位系统例如从诸如陀螺仪或加速度计之类的航向传感器获得要估计其位置或速度的主体的航向的估计。跟踪来自多个发射机的相应的信号,并且从所述被跟踪信号中的每个获得相应的多普勒测量结果。针对被跟踪信号中的每个,使用主体的航向的估计和相应的多普勒测量结果来估计主体的速度。然后确定估计的主体速度是否与已经沿着直接路径从发射机接收到的信号一致。然后当估计主体的位置或速度时忽视提供关于与已沿着直接路径从发射机接收到的信号不一致的主体速度的信息的信号。
Description
技术领域
本发明涉及一种定位系统,并且特别地涉及用于交通工具的定位系统。
背景技术
导航和定位系统被常常使用,例如使用可用全球导航卫星系统(GNSS)中的一个,诸如全球定位系统(GPS)。在此系统中,用户设备从多个人造卫星(satellite vehicle)接收信号,其中的每一个具有预定轨道。因此,用户设备能够对接收信号进行测量,并且能够使用这些来提供用户位置和速度的估计。
此系统的一个众所周知的问题是来自人造卫星的信号可能在被在人造卫星与用户设备之间的视线路径中或接近视线路径的障碍物反射出去或折射之后到达用户设备。这些障碍物还可能具有直接视线信号被显著衰减的效果,可能达到用户设备实际上不能将其接收的程度。
此问题的效果可能是用户设备不能进行用户位置和速度的良好估计。
US-7,702,459描述了供在其中系统可能由于反射而接收到多路径信号的环境中使用的定位系统。在所述系统中,将GPS测量结果与如当前基于较早GPS测量结果而估计的设备的位置和速度相比较。可以将与当前估计不一致的测量结果(例如由于其是基于反射信号)视为异常值(outlier),并且在生成设备的位置和速度的下一估计时将其忽视(disregard)。这假设系统可以不具有潜在障碍物的知识并且由此的接收信号上的反射的影响的知识。
发明内容
根据本发明的一个方面,提供了一种定位系统的操作方法。获得要估计其位置或速度的主体(body)的航向(heading)的估计。跟踪来自多个发射机相应的信号,并且从所述被跟踪信号中的每个获得相应的多普勒测量结果。针对被跟踪信号中的每个,使用主体的航向估计和相应的多普勒测量结果来估计主体的速度,并且确定主体的估计速度是否与已沿着直径路径从发射机接收到的信号一致。然后当估计主体的位置或速度时忽视提供关于与已沿着直接路径从发射机接收到的信号不一致的主体速度的信息的信号。
该方法可包括通过确定估计的主体速度是否暗示着经由反射而从发射机获得相应的信号来确定估计的主体速度是否与已沿着直接路径从发射机接收到的信号一致。
该方法可包括通过将在所述估计航向的方向上的估计的主体速度与阈值速度相比较来确定估计的主体速度是否与已沿着直接路径从发射机接收到的信号一致。
阈值速度可以是零,使得如果在所述估计航向的方向上估计的主体速度为负,则确定估计的主体速度并不与已沿着直接路径从发射机接收到的信号一致。
基于要估计其位置或速度的估计的主体航向与从发射机的直接路径之间的角度,阈值速度可以是可修改的。
主体航向的估计可基于从在主体上提供的航向传感器接收到的信号。
来自航向传感器的信号还可包括关于主体速度的信息,并且其中该方法包括通过将使用主体航向估计和相应的多普勒测量结果获得的所述主体的所述估计速度与关于从航向传感器获得的主体的速度的信息相比较来确定使用主体的航向估计和相应的多普勒测量结果获得的估计的主体速度是否与已沿着直接路径从发射机接收到的信号一致。
可在确定主体在被确定固定不动之后已开始移动之后立即获得主体航向的估计。
该方法还可以包括防止被忽视信号的进一步跟踪。
来自未被忽视的信号的测量结果可用来更新用于主体的位置/速度/时间解(timesolution)。
航向传感器可包括具有已知的与主体的位置关系的陀螺仪,用于指示主体的取向,并且然后还可包括温度传感器,用于对陀螺仪的输出进行校准。
航向传感器可包括用于指示主体的运动方向的加速度计或磁力计。
根据本发明的第二方面,提供了一种定位系统,具有用于生成要估计其位置或速度的主体的航向估计的航向传感器以及用于检测来自多个发射机的信号的接收机。该系统还包括用于执行根据本发明的第一方面的方法的处理器。
附图说明
现在将参考附图仅作为示例来描述优选实施例,在所述附图中:
图1是包括导航和定位系统的电子设备的示意性框图。
图2图示出图1的系统的可能部署。
图3图示出一个特殊情况下的由导航和定位系统接收到的信号上的反射效果。
图4图示出更一般情况下的由导航和定位系统接收到的信号上的反射效果。
图5图示出由导航和定位系统接收到的信号上的反射的更多的效果。
图6是图示出由导航和定位系统执行的方法的流程图。
具体实施方式
图1图示出在电子设备10中提供的导航和定位系统。作为示例,电子设备10可能是智能电话、平板电脑或笔记本电脑等。作为另一示例,电子设备10可能是提供有用于安装到交通工具的装置的专业化卫星导航设备。
导航和定位系统在图1中被示为被分配在多个单独处理设备之间。在一个示例中,电子设备10包括例如全球定位系统(GPS)芯片组12的全球导航卫星系统(GNSS)和也用来提供电子设备10的其它所需功能的主处理器14。然而,能够在单一集成电路中提供所有处理也是可能的。
电子设备10还包括至少一个附加传感器16,如下面更详细地描述的。
GNSS芯片组12包括用于接收由人造卫星发射的信号的射频(RF)前端电路22。接收信号被传递至以测量引擎24形式的软件,其可以被认为包括跟踪控制器26和测量处理器28。跟踪控制器26控制RF前端电路的操作,使得其从人造卫星接收预期信号,并且由测量处理器28提取用于确定设备的位置和速度所需的测量结果。
主处理器14运行以定位引擎32形式的软件,其包括控制器34和位置速度时间(PVT)处理器36。控制器34协调定位引擎32和测量引擎24,而PVT处理器36使用从GNSS芯片组12且从附加传感器16接收到的测量结果来计算设备的位置和速度的估计以及时间。虽然在本文中参考其在GNSS卫星导航和定位系统中的使用来描述本发明,但其能够同样地在其中设备从发射机接收信号的其它系统中使用,所述发射机的位置(和速度,如果适当的话)是已知的。例如,此发射机可能采取GNSS伪卫星(pseudolite)、室内通讯系统(IMES)发射机、蓝牙发射机、WiFi发射机、蜂窝式无线发射机或近场通信(NFC)发射机的形式。
如一般地理解的,控制RF前端电路22,使得其从多个发射机(诸如所示示例中的人造卫星)接收信号。由测量引擎24从这些信号中提取测量结果,并且供应给定位引擎32,其然后计算设备的位置和速度的估计以及时间。
通常,测量结果涉及由接收机接收信号所处的时间。如果接收机具有被已知完全与卫星时间同步的时钟,则接收信号所处的时间将指示接收机与人造卫星的距离,即范围。通过从三个或更多人造卫星获得范围测量结果,能够计算接收机的位置。实际上,大多数导航设备不具有此时钟,并且因此时间偏移是通过使用关于接收到信号所处时间的附加测量结果而寻求计算的参数中的一个。使用非同步接收机时钟计算的范围测量结果被称为伪距离。通常,从四个或更多人造卫星的伪距离测量结果连同接收机时钟时间偏移被用来计算接收机的三维位置。
测量结果还可以涉及从卫星接收到的每个信号的载波的相位。这提供能够在计算设备的位置和速度以及时间时使用的附加信息。
设备10的移动引起由接收机接收到的信号中的多普勒效应。也就是说,可以将接收信号的频率与RF前端电路22中的本地振荡器的频率相比较。接收信号的表观频率以取决于设备朝向或远离卫星的移动速度的分量的量而移位。其还取决于本地振荡器相对于卫星的时间帧的频率误差。本地振荡器的频率误差对多普勒测量结果的影响类似于接收机的时间偏移对范围测量结果的影响。然而,多普勒测量结果通常还被称为‘多普勒’而不是‘伪多普勒’。多普勒测量结果可以用来计算设备10的移动速度。通常,来自四个或更多人造卫星的多普勒测量结果连同接收机时钟频率偏移被用来计算接收机的三维速度。还可能使用重复多普勒测量结果来确定设备10的位置。例如,如果使用在一个时间下的多普勒测量结果来计算设备10的速度,则这可以被用来估计该时段中的设备移动,直至进行下一多普勒测量为止,并且可以使用此过程的连续重复来获得设备所遵循的路径的估计。
图2图示出其中可使用电子设备10的情况。
具体地,图2(a)是包含建筑物52、54的城市环境中的以速度S移动的交通工具50的侧视图,而图2(b)是从上面看的视图。在这里假设交通工具50包含图1中所示类型的电子设备10。出于某些目的,如下面更详细地描述的,假设电子设备10以其相对于交通工具的取向将一般地随着交通工具正在移动而保持相同的此方式固定于交通工具50是有用的。图2示出了用户50处在将从两个人造卫星60、62接收信号的位置中。将认识到的是通常必须从至少四个人造卫星接收信号以便能够获得设备的位置和速度的准确估计以及接收机时钟时间和频率偏移,但是为了明了起见在图2中仅示出了两个人造卫星,并且因为这足以举例说明本发明的操作。
图2示出了由交通工具50从第一人造卫星60接收视线信号70,并且还示出了由交通工具50在反射离开建筑物52之后从第一人造卫星60接收信号72。图2还示出了由交通工具50从第二人造卫星62接收视线信号74,并且还示出了由交通工具50在反射离开建筑物54之后从第二人造卫星62接收信号76。
如上所述,交通工具50正在以速度S在处于建筑物52与第一人造卫星62之间的直线上的方向上朝着第一人造卫星62移动。
交通工具50的移动在由GNSS芯片组12接收到的信号上引起多普勒效应,并且这些效应可以用来生成多普勒测量结果。
一般而言将多普勒测量结果理解成意指用户与人造卫星的伪矩离的变化速率。本讨论假设由人造卫星的移动引起的多普勒效应已被去除。然而,假设每个多普勒测量结果将仍受到由GNSS芯片组的RF前端电路中的振荡器所生成的信号的频率的偏移而引起的误差的影响。
图3图示出特殊情况下的多普勒测量结果上的多路径的影响。具体地,图3图示出在图2中所示的情况下能够由交通工具50从第一人造卫星60获得的多普勒测量结果中的误差。
如上所述,图2示出了其中交通工具50正在沿着在人造卫星60与反射器52之间延伸的直线以速度S移动的情况。(在示图中,S的正值表示朝向人造卫星60的移动,而S的负值表示远离人造卫星60朝向反射信号的建筑物52的移动。)多普勒测量结果涉及信号在连续测量之间行进的距离的变化。
因此,当交通工具50正在朝着人造卫星60、即在图2中向右移动时,交通工具50与第一人造卫星60之间的实际距离正在减小。然而,反射信号72所行进的距离正在增加,其正常地将暗示着交通工具50正在远离人造卫星60移动。更具体地,当交通工具50正在以速度S朝着人造卫星60移动时,多普勒测量结果指示交通工具50正在以速度S远离人造卫星60移动。因此在从多普勒测量结果中获得的速度符号方面存在由此的误差,并且在从多普勒测量结果中获得的速度量值方面存在2s的误差。
图3针对在从-1至+1范围内的速度(用适当单位表示)以图形方式图示出这种情况。针对任何给定速度,从多普勒测量结果中获得的速度误差具有相反符号,并具有两倍速度的量值。
图3适用于其中交通工具50正在沿着在人造卫星60与反射器52之间延伸的直线移动的特殊情况。在更一般情况下,反射器的位置将不是已知的。
图4图示出更一般情况下的可从多普勒测量结果获得的速度的误差。
图4假设电子设备10包括能够用来估计交通工具50的速度的至少一个附加传感器16。该附加传感器可包括加速度计,并且可另外包括陀螺仪和/或磁力计。因此,能够估计交通工具50的速度。
因此,图4示出了其中已知交通工具正在以速度S移动但未将交通工具的运动方向假设为已知的情况。具体地,未关于交通工具相对于人造卫星位置的移动方向或关于反射器相对于人造卫星与交通工具之间的视线的位置进行假设。(出于图示的目的,这对应于图2(b)中的情况,其中,交通工具正在从第二人造卫星62接收反射信号76。) 这还适用于其中反射器为三维和/或具有曲面的情况。
在图4中,水平轴线表示交通工具在朝着正在从其接收反射信号的人造卫星的方向上的速度的分量。因此,如果交通工具正在以速度S朝着人造卫星移动,则这表示+S的速度,而如果交通工具正在以速度S直接地远离人造卫星行进,则这表示-S的速度。在其它方向上的移动将表示中间速度。例如,如果交通工具正在垂直于人造卫星方向的方向上以速度S行进,则其在朝向人造卫星的方向上的速度将是零。
在图4中,垂直轴线表示由反射引起的多普勒测量结果中的误差。假设多普勒测量结果指示交通工具相对于人造卫星的位置的变化。然而,当其是正在被跟踪的反射信号时,多普勒测量结果实际上正在指示交通工具相对于反射器的位置的变化。因此,由这样的事实引起误差,即交通工具的移动引起交通工具相对于反射器的位置的变化,其不同于交通工具相对于人造卫星的位置的变化。
图4示出了用于交通工具的不同行进方向和用于反射器的不同位置的可能误差的范围。例如,线80对应于图3中所示的特殊情况,其中交通工具直接位于人造卫星与反射器之间,并且因此存在从多普勒测量结果获得的速度的符号方面的误差以及从多普勒测量结果获得的速度的量值方面的误差,并且该误差的量值等于在朝向人造卫星的方向上的速度的分量的两倍。
作为另一示例,线82对应于其中交通工具正在直接地朝着反射器移动的情况的范围,并且因此交通工具的速度具有在朝向人造卫星的方向上的分量(在从0至S的范围内),但是在朝向反射器的方向上具有S的速度,并且因此多普勒测量结果过高估计了交通工具的速度。
作为另一示例,线84对应于针对反射器的不同位置的其中交通工具正在直接地朝着人造卫星移动(且因此在该方向上具有速度S)的情况。因此,直接在交通工具“后面”(如从人造卫星看)的反射器引起多普勒测量结果中的-2S的误差,而直接邻近于人造卫星与交通工具之间的视线路径的反射器引起多普勒测量结果中的非常小的误差。
出于本文所述方法的目的,重要的是认识到当交通工具的速度为零时,反射器的存在不引起多普勒测量结果中的误差。(这忽略了由于由卫星运动产生的卫星-反射器-接收机的几何结构方面的变化而引起的小的效果。) 认识到无论交通工具的移动方向如何且无论反射器的形状和位置如何、多普勒测量结果中的误差落在图4中的阴影平行四边形内也是重要的。
在本发明的一个实施例中,假设电子设备10被以大体上固定的取向安装到交通工具50。还可以假设交通工具的运动方向一般地将基本上向前或基本上向后,并且侧向移动是不可能的。
因此,可以使用至少一个附加传感器16来估计交通工具50的运动方向,并且因此充当航向传感器。该附加传感器可包括加速度计,并且可另外包括陀螺仪和/或磁力计。该附加传感器因此可以用来提供交通工具50的速度的估计。
图5图示出使用从反射信号导出的多普勒测量结果来估计交通工具的速度的效果,其中,我们假设影响多普勒测量结果的所有其它误差是零,并且我们假设接收机时钟的频率偏移是零。在图5中,水平轴线表示交通工具相对于人造卫星的移动方向,并且0°表示直接地朝向人造卫星的移动且180°表示直接地远离人造卫星的移动。该效果关于视线是对称的,并且因此图5仅示出了用于0°和180°之间的方向的效果,因为180°和360°之间的方向产生相同效果。在图5中,垂直轴线表示能够根据反射器的位置基于多普勒测量结果来估计的用户速度。图5呈现了通过考虑处于非常低的高度的人造卫星的情况且通过假设我们知道交通工具的垂直速度是零的略微简化的分析。
在图5中,将用户的速度取为一个单位,如由线130所示。在这种情况下,阴影区域132指示根据反射器位置能够从用于交通工具的给定移动方向的单一多普勒测量结果导出的速度范围。
因此,如先前所述,如果交通工具正在直接地朝着人造卫星移动(即,处在水平轴线上的0°的位置处),并且直接在交通工具“后面”存在反射器(如从人造卫星看),则多普勒测量结果将导致-1单位的估计速度,而直接地邻近于人造卫星与交通工具之间的视线路径的反射器将仅导致多普勒测量结果中的非常小的误差。在其它位置上的反射器将导致在从-1至+1单位范围内的估计速度。
针对交通工具的其它移动方向(即,在水平轴线上的其它位置处),反射信号的使用能够引起估计速度中的甚至更大误差。特别地,当交通工具正在与到人造卫星的其视线路径成90°的方向上行进时,针对某些反射器位置,从反射信号获得的单一多普勒测量结果的使用能够产生趋向于无穷大的估计速度。
然而,将注意到的是在大多数情况下速度估计将是负的,或者其量值将被低估。因此,如果使用反射信号来估计用户的速度,并且又使用此估计来估计移动用户的位置,则这具有用户看起来“被推回”的效果,即看起来比实际情况已行进更短的距离,因为速度的估计量值低于实际值或者因为运动方向相反。
将特别地注意的是,对于许多反射器位置而言,使用从反射信号获得的多普勒测量结果估计的速度是负的(即其在与交通工具的实际移动方向相反的方向上)。
因此,如果例如根据在交通工具上提供的单独航向传感器从而交通工具的实际移动方向是已知的,则容易认识到使用此多普勒测量结果估计的速度与已知移动方向不一致。因此,能够根据图5确定多普勒测量结果可能是从反射信号获得的。
如果在交通工具上还提供了单独速度传感器,则可能识别其中能够认识到使用单一多普勒测量结果估计的速度与根据附加传感器估计的速度不一致的其它情况(对应于反射器位置的更宽范围)。在这种情况下,根据图5能够确定多普勒测量结果大概是从反射信号获得的。
图6是图示出根据本发明的实施例的方法的流程图。图6中所示的过程可以在测量引擎与定位引擎之间分裂,或者可以在单一设备中执行,或者可以以不同的方式在不同的设备之间分裂。例如,能够主要地在定位引擎中执行图6中所示的过程,并且测量引擎在定位引擎的控制下简单地向定位引擎提供所需测量结果。
该过程在其中确定交通工具是正在移动还是固定不动的步骤140中开始。此确定可以基于GNSS测量结果和/或从诸如陀螺仪、加速度计或磁力计之类的附加传感器接收到的信号。
如果确定交通工具正在移动,则GNSS测量结果和陀螺仪信号在顺序估计过程中被组合以提供交通工具航向的估计。如果传感器包括加速度计,则还可以在该估计过程中使用从那里接收到的信号。该测量结果还在步骤142中被用来更新各种功能或校准参数。例如,传感器还可包括温度传感器,在这种情况下,可以使用温度测量结果以便改善陀螺仪‘零级’偏移的温度相关模型的系数估计以及GNSS振荡器频率偏移的温度相关模型的系数。
更具体地,陀螺仪可能产生随温度而变的输出信号,并且能够对此进行建模的其中(在许多之中的)一个方式是:
其中
Rt是实际旋转速率,
G是陀螺仪输出,
T是温度,
S是缩放因数,
e是残余误差,以及
A0和A1是校准系数。
通过在不同温度下测量陀螺仪输出值,可能获得用于系数A0和A1的值,使得能够在后续测量的场合获得用于实际旋转速率的更准确值。
类似地,RF前端电路中的振荡器可以在随温度而变的频率下产生输出信号,并且能够对此进行建模的其中(在许多之中的)一个方式是:
其中
Ft是实际频率,
Fa是标称频率,
T是温度,
e是残余误差,以及
B0和B1是校准系数。
通过将输出频率值与不同温度下的接收信号的频率相比较,可能获得用于系数B0和B1的值,使得能够在后续测量场合获得用于实际频率的更准确值。
因此,在交通工具正在移动的时候,以这样的方式(即相对于交通工具对陀螺仪进行校准)来将GNSS和传感器测量结果组合。也就是说,估计参数,使得来自陀螺仪的数据单独地可以用来估计交通工具的航向。更具体地,系统获悉哪些陀螺仪信号与交通工具的移动相关联。如上文所讨论的,情况通常将是交通工具将正在相对于其自己的轴线向前或向后行进,具有移动的相对小的侧向分量。将不一定最初知道设备10被如何安装在交通工具中或其上面,并且因此陀螺仪的输出最初将不会与交通工具的任何特定方向上的移动相关联,但是这将在此校准之后变成可能。
如果在步骤140中确定交通工具是固定不动的,则过程转到步骤144,其中,利用此‘零速度’知识,执行陀螺仪零级偏移和GNSS振荡器频率偏移的‘快照(snapshot)’校准。
因此,通过在当已知实际旋转速率是零时的时间下测量陀螺仪输出值,可能获得如上文所定义的用于系数A0和A1的值,使得能够在后续测量场合获得用于实际旋转速率的更准确值。
类似地,通过将输出频率值与在当交通工具已知固定不动并且因此没有多普勒误差时的时间下接收信号的已知频率相比较,可能获得用于系数B0和B1的值,使得能够在后续测量场合获得用于实际频率的更准确值。
在此校准之后,过程转到步骤146,其中确定交通工具是否已开始移动。如果没有,则过程重复步骤144,直至检测到移动为止,并且过程转到步骤148。
在步骤148中,使用来自陀螺仪或诸如磁力计之类的其它航向传感器的数据连同来自交通工具航向和各种功能或校准参数的先前更新的估计来预测当前时间下的交通工具航向。还可以根据加速度计数据来估计交通工具速度。
另外,同时地或单独地,使用GNSS和(可选地)温度数据连同来自各种功能或校准参数的先前更新的估计来预测在当前时间下的GNSS振荡器频率偏移。来自先前更新的估计可以是来自在步骤144中执行的最近固定不动的快照校准的估计,但是其它实施方式是可能的。
因此,在当前更新时间下,进行用户航向、时钟漂移和可选地速度的估计,并且这些估计独立于或者仅微弱地取决于当前GNSS测量结果。
该过程然后转到步骤150,其中考虑最近GNSS测量结果。例如,在其中该处理被分裂在这两个模块之间的情况下,可以将这些测量结果从测量引擎发送到定位引擎。
如将显而易见的是,RF前端将从多个人造卫星接收信号。在反射是因素的情况下,RF前端将从至少一些人造卫星接收多于一个的信号。每个信号能够被用来生成多普勒测量结果。
如先前所讨论的,多普勒测量结果受到接收机的振荡器的假定频率中的误差的影响。通过检查当已知接收机为固定不动时接收到的信号的频率,能够对本地振荡器信号的频率进行校准。然后,可以假设本地振荡器的已知频率与当已知接收机正在移动时接收到的信号的频率之间的任何差异是由于该移动而引起的。因此,可以使用多普勒测量结果来估计水平面中的交通工具的速度在人造卫星方向上的分量,假设交通工具的任何垂直速度是小的。然后可以将其与用户相对于卫星方向的航向方向的知识组合,以便生成用户在其航向方向上的速度估计。将认识到的是由于多普勒测量结果中的误差的影响,此速度估计不可能是特别准确的。还将认识到的是此速度估计的准确度将高度取决于用户相对于卫星方向的航向的方向。然而,应强调的是此估计并不是接收机速度的最终估计。更确切地说,其是判定哪些测量结果可能已经从反射信号导出的手段的一部分,使得能够将此测量结果从接收机位置和速度的后续、最终估计中排除。
该过程然后转到步骤152,其中检查交通工具速度的这些单独估计中的每一个以确定其是否与已被反射的卫星信号一致。例如,这可以简单地涉及到确定速度估计是否是负的。也就是说,将速度估计与在步骤148中获得的运动方向相比较。如果当在步骤148中找到的运动方向具有远离人造卫星的分量(或者反之亦然)时发现从多普勒测量结果导出的速度估计具有朝向人造卫星的分量,则确定该测量结果可能已经从反射信号中导出。
更一般地,交通工具速度的每个单独估计的检查可能涉及到将该估计与阈值相比较。
因此,测试是否
S<x
其中:
S是基于用户相对于卫星方向的测量航向从单一多普勒测量结果获得的交通工具速度的估计,并且x是阈值值。
在上述的一个示例中,可将阈值值x设置成零,使得交通工具速度的单独估计被确定成与当速度估计为负时已经被反射的卫星信号一致。
在另一示例中,可将阈值值x设置成负常数,使得交通工具速度的单独估计被确定成与当速度估计比该常数更加为负时已经被反射的卫星信号一致。
在另一示例中,可将阈值值x设置成负常数与表示测量结果不确定性的值的和,使得基于测量结果的预期准确度来修改阈值。
作为另一示例,阈值值x可被修改以考虑这样的事实,即当交通工具正在以接近于垂直于其到卫星的方向的航向移动时,交通工具的速度估计可能是质量相对差的。这能够用来设置阈值以改变(增加或降低)交通工具速度的估计被确定为与已被反射的卫星信号一致的可能性。
如果发现从特定信号导出的多普勒测量结果与已沿着直接路径从卫星接收到的不一致(即其与已沿着反射路径接收到的更加一致),则该特定多普勒测量结果被标记(flag),以由所有后续处理拒绝。
具体地,当计算交通工具的位置和速度且可能还计算时间(即PVT解的分量)时,在后续步骤154中忽视该多普勒测量结果。该计算步骤一般是常规的,并且将从多个接收信号获得的测量结果组合,以便达到交通工具的位置和/或速度的估计,但是忽视从被确定为可能已经被反射的信号获得的测量结果。
另外,还可忽视从同一信号获得的伪距离测量结果。此外,定位引擎可指令测量引擎停止跟踪该信号,因此没有进一步的测量结果从反射信号中被导出。
因此,执行反射聚焦的异常值拒绝过程。已知当使用GNSS测量结果以及传感器输入等来计算PVT解时拒绝异常值测量结果。此已知过程通常考虑在一个时间接收到的一组测量结果,并且假设这些测量结果中的仅小的比例将是错误的。因此其假设可以通过与其它测量结果不一致来识别错误的测量结果。在这里所述的异常值拒绝过程中,与其它GNSS测量结果隔离地考虑每个GNSS测量结果,以确定其是否与从其它源、特别是从航向传感器获得的信息一致。如果测量结果通过此测试,则然后将其在PVT解的计算中使用,认识到此计算可包括另外的异常值拒绝步骤。
图6中所示的过程然后转到步骤156,其中考虑到所计算的PVT解,以与参考步骤144所述的相同方式来更新参数估计。
因此描述了特别是交通工具在其中GNSS信号高度易于受到来自建筑物等的反射的影响的环境中的情况下,允许以更大的精确度来确定用户的位置和/或速度的过程。
Claims (20)
1.一种定位系统的操作方法,该方法包括:
获得要估计其位置或速度的主体的航向的估计,
跟踪来自多个发射机的相应的信号;
从所述被跟踪信号中的每个中获得相应的多普勒测量结果;
针对所述被跟踪信号中的每个,使用主体的航向估计和相应的多普勒测量结果来估计所述主体的速度,确定估计的主体速度是否与已沿着直接路径从发射机接收到的信号一致;以及
当估计主体的位置或速度时,忽视提供关于与已沿着直接路径从发射机接收到的信号不一致的主体速度的信息的信号。
2.如权利要求1所述的方法,进一步包括:
通过确定估计的主体速度是否暗示着经由反射而从发射机获得相应的信号来确定估计的主体速度是否与已沿着直接路径从发射机接收到的信号一致。
3.如权利要求1所述的方法,进一步包括:
通过将在所述估计航向的方向上的估计的主体速度与阈值速度相比较来确定估计的主体速度是否与已沿着直接路径从发射机接收到的信号一致。
4.如权利要求3所述的方法,其中所述阈值速度是零,使得如果在所述估计航向的方向上估计的主体速度是负的,则确定估计的主体速度与已沿着直接路径从发射机接收到的信号不一致。
5.如权利要求3所述的方法,其中,基于要估计其位置或速度的估计的主体航向与从发射机的直接路径之间的角度,所述阈值速度是可修改的。
6.如权利要求1所述的方法,其中,主体航向的估计是基于从在主体上提供的航向传感器接收到的信号。
7.如权利要求6所述的方法,其中,来自航向传感器的信号还包括关于主体速度的信息,该方法进一步包括:通过将使用主体航向的估计和相应的多普勒测量结果获得的所述主体的估计的速度与从航向传感器获得的关于主体速度的信息相比较来确定使用主体航向的估计和相应多普勒测量结果获得的所述主体的所述估计的速度是否与已沿着直接路径从发射机接收到的信号一致。
8.如权利要求1所述的方法,进一步包括在确定主体已在被确定为固定不动之后开始移动之后立即获得主体的航向的估计。
9.如权利要求1所述的方法,还包括:
防止对被忽视信号的未来跟踪。
10.如权利要求1所述的方法,还包括:
使用来自未被忽视的信号的测量结果来更新用于主体的位置/速度/时间解。
11.如权利要求1所述的方法,其中,所述航向传感器包括陀螺仪,其具有已知的与主体的位置关系,用于指示主体的取向。
12.如权利要求11所述的方法,其中,所述航向传感器还包括用于对陀螺仪的输出进行校准的温度传感器。
13.如权利要求10所述的方法,还包括用于指示主体的运动方向的加速度计。
14.如权利要求1所述的方法,其中,航向传感器包括用于指示主体的运动方向的磁力计。
15.一种定位系统,包括:
航向传感器,用于生成要估计其位置或速度的主体的航向的估计;
接收机,用于跟踪来自多个发射机的相应的信号;以及
处理器,用于从所述被跟踪信号中的每个导出相应的多普勒测量结果,使用主体的航向估计和针对所述被跟踪信号中的每个的相应的多普勒测量结果来估计所述主体的速度,确定估计的主体速度是否与已沿着直接路径从发射机接收到的信号一致,以及当估计主体的位置或速度时,忽视提供关于与已沿着直接路径从发射机接收到的信号不一致的主体速度的信息的信号。
16.如权利要求15所述的定位系统,其中具有多普勒测量结果导出的与已沿着直接路径从所述多个发射机中的相应的发射机接收到的信号一致的速度的信号是与所述航向一致的。
17.如权利要求16所述的定位系统,其中具有多普勒测量结果导出的与已沿着直接路径从所述多个发射机中的相应的发射机接收到的信号不一致的速度的信号被忽视。
18.如权利要求15所述的定位系统,其中航向传感器包括陀螺仪,该陀螺仪具有已知的与所述主体的位置关系,用于指示主体的取向。
19.如权利要求18所述的定位系统,其中航向传感器进一步包括温度传感器,用于校准陀螺仪的输出。
20.如权利要求19所述的定位系统,其中航向传感器包括磁力计,用于指示主体的运动的方向。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP13173075.6A EP2816374B1 (en) | 2013-06-20 | 2013-06-20 | Vehicle positioning in high-reflection environments |
EP13173075.6 | 2013-06-20 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104237920A CN104237920A (zh) | 2014-12-24 |
CN104237920B true CN104237920B (zh) | 2017-04-12 |
Family
ID=48628570
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410279249.0A Expired - Fee Related CN104237920B (zh) | 2013-06-20 | 2014-06-20 | 高反射环境中的交通工具定位 |
Country Status (4)
Country | Link |
---|---|
US (1) | US20140375495A1 (zh) |
EP (1) | EP2816374B1 (zh) |
JP (1) | JP2015004678A (zh) |
CN (1) | CN104237920B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102015203616A1 (de) | 2015-02-27 | 2016-09-01 | Continental Teves Ag & Co. Ohg | Multipath Erkennung mittels Vergleich zweier unterschiedlicher GNSS Signalen |
DE102015218810A1 (de) | 2015-09-29 | 2017-03-30 | Continental Teves Ag & Co. Ohg | Verfahren zum Auswählen von Lokalisierungsalgorithmen in einem Fahrzeug |
CN105711483A (zh) * | 2016-04-20 | 2016-06-29 | 高碑店市博伟汽车配件有限公司 | 一种弯道辅助照明系统 |
US20180011200A1 (en) * | 2016-07-06 | 2018-01-11 | Qualcomm Incorporated | Satellite signal exclusion based on doppler information |
US10222472B2 (en) * | 2016-09-30 | 2019-03-05 | Veoneer Us, Inc. | System and method for detecting heading and velocity of a target object |
DE102016220729A1 (de) * | 2016-10-21 | 2018-04-26 | Robert Bosch Gmbh | Verfahren, System und Computerprogrammprodukt zur Verbesserung der Positionsbestimmung mittels Satellitennavigation |
JP6946660B2 (ja) * | 2017-02-23 | 2021-10-06 | 株式会社デンソー | 測位装置 |
DE112019006954T5 (de) * | 2019-03-01 | 2021-12-16 | Google Llc | Geschwindigkeitsbestimmung mit hilfe eines reflektierten positionierungssignals |
WO2022133294A1 (en) | 2020-12-17 | 2022-06-23 | Swift Navigation, Inc. | System and method for fusing dead reckoning and gnss data streams |
US11906640B2 (en) | 2022-03-01 | 2024-02-20 | Swift Navigation, Inc. | System and method for fusing sensor and satellite measurements for positioning determination |
WO2023167916A1 (en) | 2022-03-01 | 2023-09-07 | Swift Navigation, Inc. | System and method for detecting outliers in gnss observations |
US12019163B2 (en) | 2022-09-12 | 2024-06-25 | Swift Navigation, Inc. | System and method for GNSS correction transmission |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101088452A (zh) * | 2006-06-12 | 2007-12-19 | 奥林巴斯医疗株式会社 | 内窥镜插入形状检测装置 |
CN101606080A (zh) * | 2007-02-26 | 2009-12-16 | 丰田自动车株式会社 | 移动单元定位装置 |
CN102519450A (zh) * | 2011-12-12 | 2012-06-27 | 东南大学 | 一种用于水下滑翔器的组合导航装置及方法 |
Family Cites Families (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0666916A (ja) * | 1992-08-13 | 1994-03-11 | Fujitsu Ten Ltd | Gps受信機 |
JPH0843515A (ja) * | 1994-08-02 | 1996-02-16 | Nippondenso Co Ltd | Gps受信装置 |
JP3851376B2 (ja) * | 1996-05-15 | 2006-11-29 | 日本無線株式会社 | 測位システム衛星信号受信機 |
JP2001264409A (ja) * | 2000-03-21 | 2001-09-26 | Clarion Co Ltd | カーナビゲーション装置におけるマルチパス判定方法 |
US6693592B2 (en) * | 2000-12-22 | 2004-02-17 | The Charles Stark Draper Laboratory, Inc. | Geographical navigation using multipath wireless navigation signals |
JP2003057327A (ja) * | 2001-08-09 | 2003-02-26 | Matsushita Electric Ind Co Ltd | 航法衛星信号受信機 |
JP3705187B2 (ja) * | 2001-11-12 | 2005-10-12 | 株式会社デンソー | ナビゲーション装置 |
JP2004045126A (ja) * | 2002-07-10 | 2004-02-12 | Matsushita Electric Ind Co Ltd | 衛星信号受信機 |
US8239162B2 (en) * | 2006-04-13 | 2012-08-07 | Tanenhaus & Associates, Inc. | Miniaturized inertial measurement unit and associated methods |
JP4657050B2 (ja) * | 2005-07-29 | 2011-03-23 | 日本無線株式会社 | 速度ベクトル決定システム |
JP2008051572A (ja) * | 2006-08-23 | 2008-03-06 | Matsushita Electric Ind Co Ltd | ナビゲーション装置及び、その方法、並びにそのプログラム |
US7702459B2 (en) * | 2006-10-17 | 2010-04-20 | Alpine Electronics, Inc. | GPS accuracy adjustment to mitigate multipath problems for MEMS based integrated INS/GPS navigation systems |
JP2008139247A (ja) * | 2006-12-05 | 2008-06-19 | Alpine Electronics Inc | Gps測位装置 |
US8286463B2 (en) * | 2007-04-04 | 2012-10-16 | Qualcomm Incorporated | Method for determining the null point of a gyroscope |
US7821454B2 (en) * | 2007-11-20 | 2010-10-26 | Sirf Technology, Inc. | Systems and methods for detecting GPS measurement errors |
JP5164546B2 (ja) * | 2007-12-07 | 2013-03-21 | クラリオン株式会社 | 測位方法 |
US8044851B2 (en) * | 2009-04-17 | 2011-10-25 | Mstar Semiconductor, Inc. | Method for suppressing multipath errors in a satellite navigation receiver |
JP5352422B2 (ja) * | 2009-10-30 | 2013-11-27 | 株式会社豊田中央研究所 | 測位装置及びプログラム |
FR2962227B1 (fr) * | 2010-07-01 | 2012-07-27 | Inst Telecom Telecom Sudparis | Procede de reduction de l'eblouissement d'un recepteur recevant des signaux depuis des emetteurs |
FR2971857A1 (fr) * | 2011-02-17 | 2012-08-24 | Thales Sa | Procede et systeme de determination des parametres de navigation d'un aeronef |
US9903956B2 (en) * | 2011-09-12 | 2018-02-27 | Continental Teves Ag & Co. Ohg | Method for selecting a satellite |
-
2013
- 2013-06-20 EP EP13173075.6A patent/EP2816374B1/en not_active Not-in-force
-
2014
- 2014-06-19 JP JP2014126024A patent/JP2015004678A/ja active Pending
- 2014-06-20 US US14/309,950 patent/US20140375495A1/en not_active Abandoned
- 2014-06-20 CN CN201410279249.0A patent/CN104237920B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101088452A (zh) * | 2006-06-12 | 2007-12-19 | 奥林巴斯医疗株式会社 | 内窥镜插入形状检测装置 |
CN101606080A (zh) * | 2007-02-26 | 2009-12-16 | 丰田自动车株式会社 | 移动单元定位装置 |
CN102519450A (zh) * | 2011-12-12 | 2012-06-27 | 东南大学 | 一种用于水下滑翔器的组合导航装置及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104237920A (zh) | 2014-12-24 |
US20140375495A1 (en) | 2014-12-25 |
JP2015004678A (ja) | 2015-01-08 |
EP2816374B1 (en) | 2016-05-25 |
EP2816374A1 (en) | 2014-12-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104237920B (zh) | 高反射环境中的交通工具定位 | |
US8421670B2 (en) | Position estimation apparatus and computer readable medium storing position estimation program | |
US10267920B2 (en) | Positioning method | |
JPWO2014002211A1 (ja) | 測位装置 | |
JP2006126181A (ja) | 移動体姿勢検出装置 | |
JP2012203721A (ja) | 相対位置推定装置及びプログラム | |
JP2010256301A (ja) | マルチパス判定装置及びプログラム | |
JP6532126B2 (ja) | 位置特定装置、ナビゲーション装置、位置特定方法、および位置特定プログラム | |
JP2024016253A (ja) | 位置計測装置、位置計測方法、及びプログラム | |
JP2004069536A (ja) | データ検定装置および方法 | |
JP2012098185A (ja) | 方位角推定装置及びプログラム | |
WO2020149014A1 (ja) | 衛星選択装置、及びプログラム | |
US20220026583A1 (en) | Enhancing Sensitivity to Reflected GNSS Signals | |
WO2020110996A1 (ja) | 測位装置、速度測定装置、及びプログラム | |
JP2007225408A (ja) | 移動体の横滑り計測装置 | |
US9778369B2 (en) | Pedestrian positioning in high-reflection environments | |
JP4496891B2 (ja) | 位置検出方法及び位置検出装置 | |
JP2010060421A (ja) | 移動体用測位システム及びgnss受信装置 | |
JP7464113B2 (ja) | 位置計測装置、位置計測方法、及びプログラム | |
JP3019719B2 (ja) | 位置測定方法 | |
JPH07306253A (ja) | 移動情報出力装置 | |
US20220043166A1 (en) | Determining Velocity Using a Reflected Positioning Signal | |
JP2023019895A (ja) | 受信信号の到来方向と衛星方位が不一致となる衛星の検出方法および検出装置ならびにgnss受信機 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170412 Termination date: 20190620 |
|
CF01 | Termination of patent right due to non-payment of annual fee |