CN105934786A - 异常行驶部位检测装置以及异常行驶部位检测方法 - Google Patents

异常行驶部位检测装置以及异常行驶部位检测方法 Download PDF

Info

Publication number
CN105934786A
CN105934786A CN201480073508.2A CN201480073508A CN105934786A CN 105934786 A CN105934786 A CN 105934786A CN 201480073508 A CN201480073508 A CN 201480073508A CN 105934786 A CN105934786 A CN 105934786A
Authority
CN
China
Prior art keywords
abnormal
turning
value
location detection
turning value
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.)
Granted
Application number
CN201480073508.2A
Other languages
English (en)
Other versions
CN105934786B (zh
Inventor
柏井忠大
根山亮
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.)
Toyota Motor Corp
Original Assignee
Toyota Motor Corp
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 Toyota Motor Corp filed Critical Toyota Motor Corp
Publication of CN105934786A publication Critical patent/CN105934786A/zh
Application granted granted Critical
Publication of CN105934786B publication Critical patent/CN105934786B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/01Detecting movement of traffic to be counted or controlled
    • G08G1/0104Measuring and analyzing of parameters relative to traffic conditions
    • G08G1/0125Traffic data processing
    • G08G1/0129Traffic data processing for creating historical data or processing based on historical data
    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/01Detecting movement of traffic to be counted or controlled
    • G08G1/0104Measuring and analyzing of parameters relative to traffic conditions
    • G08G1/0108Measuring and analyzing of parameters relative to traffic conditions based on the source of data
    • G08G1/0112Measuring and analyzing of parameters relative to traffic conditions based on the source of data from the vehicle, e.g. floating car data [FCD]
    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/01Detecting movement of traffic to be counted or controlled
    • G08G1/0104Measuring and analyzing of parameters relative to traffic conditions
    • G08G1/0137Measuring and analyzing of parameters relative to traffic conditions for specific applications
    • G08G1/0141Measuring and analyzing of parameters relative to traffic conditions for specific applications for traffic information dissemination

Landscapes

  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Traffic Control Systems (AREA)
  • Time Recorders, Dirve Recorders, Access Control (AREA)
  • Navigation (AREA)

Abstract

异常行驶部位检测装置检测异常行驶部位。异常行驶部位是发生向道路的宽度方向的异常行驶的部位。历史信息获取部从车辆获取包括转弯值的行驶历史信息。转弯值表示与所述车辆的位置及所述车辆的转弯有关的物理量。转弯值校正部关于多个所述行驶历史信息,校正这些行驶历史信息中包含的所述转弯值。所述转弯值被校正为与预先设定的每个运算对象地点关联的转弯值。基准转弯值运算部使用与相互相同的运算对象地点关联的多个所述转弯值来运算基准转弯值。基准转弯值是平均转弯值。异常行驶部位检测部关于作为异常行驶部位的检测对象的行驶历史信息,运算偏离所述基准转弯值的背离度。偏离所述基准转弯值的背离度大的部位被检测为所述异常行驶部位。

Description

异常行驶部位检测装置以及异常行驶部位检测方法
技术领域
本公开涉及检测为了躲避障碍物等而发生了向道路的宽度方向的异常的行驶的部位的异常行驶部位检测装置以及异常行驶部位检测方法。
背景技术
在近年来开展高速公路交通系统(ITS:Intelligent TransportSystems(智能交通系统))的过程中,特别是最近,例如还要求检测为了躲避障碍物等而发生了向道路的宽度方向的异常行驶的部位(异常行驶部位)等。通过这样检测异常行驶部位,例如还能够对用户提供与异常行驶部位有关的引导、以及反映出异常行驶部位的驾驶支援。
以往,作为为了检测异常行驶部位而使用的装置,有例如判断本车辆从行车道脱离的脱离判断装置(例如参照专利文献1)。脱离判断装置使用通过车载摄像装置拍摄道路而得到的图像数据、本车辆的转弯量,来判断本车辆从行车道脱离的可能性是否大。
专利文献1:日本特开2013-3913号公报
发明内容
但是,上述脱离判断装置所使用的车载摄像装置作为搭载于车辆的装置,由于比较昂贵等的理由而普及性低。因此,在从设置有脱离判断装置的不特定多个车辆收集行车道脱离信息来检测异常行驶部位的情况下,无法收集足够数量的数据,所以无法得到不依赖于个体差等一般性高的信息。要求开发能够使用普及性高的传感器来获取与异常行驶部位有关的信息的系统。
本公开的目的在于利用普及性高的传感器来获取与道路的异常行驶部位有关的信息。
与本公开的一个技术方案有关的异常行驶部位检测装置,是检测异常行驶部位的异常行驶部位检测装置,所述异常行驶部位是发生向道路的宽度方向的异常的行驶的部位,所述异常行驶部位检测装置具备:历史信息获取部,构成为从车辆获取包括转弯值的行驶历史信息,所述转弯值是表示与所述车辆的位置及所述车辆的转弯有关的物理量的值;转弯值校正部,构成为关于多个所述行驶历史信息,校正这些行驶历史信息中包含的所述转弯值,所述转弯值校正部构成为将所述转弯值校正为与预先设定的每个运算对象地点关联起来的转弯值;基准转弯值运算部,构成为使用与相互相同的运算对象地点关联起来的多个所述转弯值来运算基准转弯值,所述基准转弯值是平均转弯值;以及异常行驶部位检测部,构成为关于作为所述异常行驶部位的检测对象的行驶历史信息,运算偏离所述基准转弯值的背离度,所述异常行驶部位检测部构成为将偏离所述基准转弯值的背离度大的部位检测为所述异常行驶部位。
与本公开的另一技术方案有关的异常行驶部位检测方法,是检测异常行驶部位的异常行驶部位检测方法,所述异常行驶部位是发生向道路的宽度方向的异常的行驶的部位,所述异常行驶部位检测方法的特征在于,包括:从车辆获取包括转弯值的行驶历史信息,所述转弯值是表示与所述车辆的位置及所述车辆的转弯有关的物理量的值;关于多个行驶历史信息,校正这些行驶历史信息中包含的转弯值,所述校正将所述转弯值校正为与预先设定的每个运算对象地点关联起来的转弯值;使用与相互相同的运算对象地点关联起来的多个所述转弯值来运算基准转弯值,所述基准转弯值是平均转弯值;关于作为异常行驶部位的检测对象的行驶历史信息,运算偏离所述基准转弯值的背离度;以及将偏离所述基准转弯值的所述背离度大的部位检测为所述异常行驶部位。
根据上述结构或者方法,能够利用检测转弯值的传感器,根据基准转弯值与检测对象的行驶历史信息中包含的转弯值的背离度来检测异常行驶部位,其中基准转弯值是基于与每个运算对象地点关联起来的多个转弯值得到的平均转弯值。因此,能够利用例如偏航率传感器等普及性高的传感器来获得与道路的异常行驶部位有关的信息。
作为一个方案,所述异常行驶部位检测部构成为计算与相互相同的运算对象地点关联起来的多个转弯值的绝对偏差,所述异常行驶部位检测部进而构成为关于作为所述异常行驶部位的检测对象的所述行驶历史信息,将与所述运算对象地点关联起来的所述转弯值与所述基准转弯值的差分除以所述绝对偏差,从而计算所述背离度。
根据上述结构,在计算背离度时使用绝对偏差,所以能够提高背离度相对于偏移值的鲁棒性。
作为一个方案,所述异常行驶部位检测部构成为计算所述转弯值偏离所述基准转弯值的背离度的沿着所述运算对象地点的推移,所述异常行驶部位检测部进而构成为对计算出的所述背离度的推移实施平滑化处理。
根据上述结构,针对背离度的推移实施平滑化处理,所以能够抑制方向盘操作的摇晃等所致的转弯值的偏移值对异常行驶部位的检测造成的影响。
作为一个方案,所述异常行驶部位检测部构成为针对作为异常行驶部位的检测对象的每个行驶历史信息检测异常行驶部位,所述异常行驶部位检测部进而构成为在基于检测到所述异常行驶部位的多个行驶历史信息的相互相同的部位处的异常行驶部位的出现频度大的情况下,将所述部位最终判定为所述异常行驶部位。
在上述结构中,将每个行驶历史信息的异常行驶部位中的、出现频度大的异常行驶部位最终判定为异常行驶部位,所以能够得到不受个体限定的信息。
在一个方案中,所述行驶历史信息中包含的转弯值是车辆的偏航率。
根据上述结构,因为作为用于判断异常行驶部位的变量而使用偏航率,所以能够提高检测出的异常行驶部位的精度,所述偏航率通过与车速等信息并用而能够判断道路的宽度方向的轨迹。在行驶历史信息是从不特定多数的车辆收集到的信息的情况下,车辆发送包括偏航率的行驶历史信息即可,异常行驶部位检测装置对多个车辆的偏航率进行统计处理即可。在统计处理的对象的行驶历史信息是根据一台车辆的行驶历史收集到的信息的情况下,车辆获取包括偏航率的行驶历史信息即可。因此,相比于使用由车载摄像装置拍摄到的图像数据等来检测异常行驶部位的情况,能够降低成本,并且能够减轻搭载于车辆的控制装置、异常行驶部位检测装置所需要的运算负荷。
附图说明
图1是示出包括本公开的异常行驶部位检测装置的支援信息生成系统的概略结构的框图。
图2是示出从图1所示的车辆发送的探查信息的数据结构的图。
图3是示出图1所示的各行程的行驶轨迹的例子的示意图。
图4是示出图1所示的异常行驶部位检测装置的动作的流程图。
图5是说明在利用图1所示的异常行驶部位检测装置进行的异常行驶部位的检测中所使用的数据流的图。
图6是说明图2所示的行程中的偏航率的量化方法的图。
图7是说明计算图6所示的偏航率的平均值基准值和绝对偏差的方法的图。
图8是说明关于图8所示的1个行程信息而检测到的异常行驶部位的图。
图9是说明关于图8所示的行程信息,最终判定异常行驶部位的方法的图。
图10是关于异常行驶部位检测装置的变形例,说明在异常行驶部位的检测中所使用的数据流的图。
(符号说明)
10:异常行驶部位检测装置;11:支援系统;15:通信部;16:支援信息生成部;17:探查信息存储部;18:道路地图信息存储部;19:行程信息存储部;20:历史信息获取部;21:行程提取部;22:转弯值校正部;23:基准转弯值运算部;24:异常行驶部位检测部;50:历史发送系统;51:GPS接收部;52:车速传感器;53:偏航率传感器;54:车载通信部;55:车载控制部;100:车辆;101:探查信息;101:行程信息;N:网络。
具体实施方式
以下,说明对异常行驶部位检测装置以及异常行驶部位检测方法进行具体化的一个实施方式。在本实施方式中,将异常行驶部位检测装置具体化为构成探查车辆系统(probe car system)的一部分的装置而进行说明。探查车辆系统收集在不特定多数的车辆中生成的行驶历史信息来作为探查信息。
如图1所示,在本实施方式中,检测在道路的宽度方向上发生车辆的异常行驶的部位的异常行驶部位检测装置(以下称为检测装置10)构成支援信息生成系统11的一部分,该支援信息生成系统11生成用于提供道路信息或者驾驶支援等的支援信息。支援信息生成系统11经由网络N与搭载于车辆100的探查发送系统50连接。
说明搭载于车辆100的探查发送系统50。探查发送系统50具备接收从GPS(Global Positioning System:全球定位系统)卫星发送的电波的GPS接收部51、车速传感器52、偏航率传感器53、车载通信部54以及车载控制部55。车载控制部55根据从GPS接收部51输入的电波检测信号计算本车位置的纬度以及经度而作为绝对位置坐标。车载控制部55从车速传感器52输入车轮速度脉冲,从偏航率传感器53输入作为车辆100的转弯方向的角速度的偏航率。车载控制部55生成作为行驶历史信息的探查信息101,将生成的探查信息101经由车载通信部54而发送到支援信息生成系统11。
如图2所示,探查信息101包括车辆标识符102、时刻信息103、绝对位置坐标104、车速105、偏航率106。车辆标识符102是用于在支援信息生成系统11侧识别发送源的车辆100的标识符。时刻信息103表示生成探查信息101时的日期时间。绝对位置坐标104是根据GPS电波信号计算出的坐标。车速105是基于车载控制部55从车速传感器52输入的车轮速度脉冲而得到的车速,偏航率106是从偏航率传感器53输入的偏航率。
参照图1,说明支援信息生成系统11。支援信息生成系统11除了上述检测装置10以外,还具备用于接收探查信息101的通信部15以及支援信息生成部16。支援信息生成系统11具备探查信息存储部17、道路地图信息存储部18以及行程信息存储部19。
在本实施方式中,检测装置10具备CPU、RAM以及ROM等硬件以及用于检测异常行驶部位的程序,从而作为历史信息获取部20、行程提取部21、转弯值校正部22、基准转弯值运算部23、异常行驶部位检测部24而发挥功能。
历史信息获取部20经由通信部15获取从车辆100发送的探查信息101,将获取到的探查信息101储存到探查信息存储部17。
行程提取部21从所积蓄的探查信息101中,读入在通用的条件下生成的探查信息101。例如,行程提取部21读入在高速道路、干线道路的预定的行驶区域内生成的探查信息101、或者在预定的行驶区域内在特定的时间段生成的探查信息101等。在本实施方式中,在读入相互相同的行驶区域的探查信息101的情况下,即使是在相互不同的行车道上行驶的情况下,也作为在相互相同的行驶区域中行驶的信息处理。行程提取部21从所读入的探查信息101中的、从一台车辆100连续发送的探查信息101提取必要的数据,将把探查信息101以被发送的顺序按时系列排列而得到的信息作为行程信息110。行程是指一台车辆在某个目的的支配下,从起点(出发地)移动至终点(目的地)的单位。如果生成了行程信息110,则行程提取部21储存到行程信息存储部19中。
道路地图信息存储部18中储存有道路地图信息25。道路地图信息具有与链路以及设置在链路的两端的节点有关的信息,所述链路是将交叉路口或信号灯、岔路口等作为单位而划分出的区间。
支援信息生成部16根据从检测装置10输出的信息,生成道路信息、用于驾驶支援的信息。将所生成的信息,从通信部15经由网络N而发送到车辆100。
参照图3~图9,说明检测装置10的动作。
参照图3说明异常行驶部位。在存在与“nt”个行程Tj(1≤j≤nt)对应的车辆的行驶轨迹时,根据行驶环境,也有所有行驶轨迹都大致沿道路的区间。另一方面,在存在停车车辆、故障车辆等的躲避部位121时、在一个行车道上发生拥堵时等,如图3中区间122所示,车辆的行驶轨迹向道路的宽度方向的偏散变大,区间122是在道路120的长度方向上划分出的多个区间中的1个区间。在存在相对于沿道路的行驶轨迹的背离大的行程Tj、而且在相互相同的区间中存在多个背离大的行程Tj的情况下,将该区间检测为异常行驶部位。关于异常行驶部位,也可以并非如图3那样具有沿着道路的行进方向的长度的区间,而检测为地点。
依照图4说明用于检测异常行驶部位的主要的步骤。检测装置10获取在通用条件下生成的探查信息101(步骤S1),从探查信息101提取行程信息110(步骤2)。根据行程信息110,针对每个行程Tj进行车辆位置的辨识,从而提高车辆位置的精度(步骤S3)。
进而,由于行程信息110中包含的偏航率的测量位置并非恒定,所以检测装置10将这些偏航率与等间隔地设定的运算对象地点xi关联起来,从而进行量化(步骤S4)。进而,检测装置10使用量化后的偏航率,计算成为各运算对象地点xi的偏航率的平均的作为转弯基准值的偏航率基准值、以及偏航率的绝对偏差(步骤S5)。
接下来,检测装置10个别地读入为了计算平均以及绝对偏差而使用的行程信息110,处置为用于检测异常行驶部位的行程信息110。首先,针对每个行程信息110,计算偏离偏航率基准值的、量化后的偏航率的背离度(步骤S6)。针对每个运算对象地点xi计算背离度。检测装置10为了去掉方向盘操作的摇晃等所致的偏移值,对运算对象地点xi间的背离度的推移进行平滑化(步骤S7)。判断在表示平滑化后的背离度的推移的曲线中是否存在预先设定的阈值以上的区域,将阈值以上的区域作为行程的异常行驶部位(步骤S8)。在该阶段中,针对每个行程检测出的异常行驶部位成为异常行驶部位的候选。
检测装置10合并针对每个行程计算出的异常行驶部位,最终检测异常行驶部位(步骤S9)。
参照图5,说明在上述动作中使用的数据流。在计算偏航率基准值以及绝对偏差时,多个行程信息110被用作统计处理的对象。接下来,使用该偏航率基准值以及绝对偏差,针对每个行程信息110个别地计算偏离偏航率基准值的背离度。在本实施方式中,个别地读入用作统计处理的对象的行程信息110来计算背离度。在计算出每个行程信息110的背离度时,根据背离度,针对每个行程信息110检测异常行驶部位的候选。使用检测到异常行驶部位的候选的多个行程信息110,根据相互相同的部位的异常行驶部位的出现频度等,最终判定异常行驶部位。
接下来,详细阐述各步骤。在步骤S1中,如上所述,检测装置10从探查信息存储部17所积蓄的探查信息101中获取车辆100的行进方向相互相同、且在通用条件下生成的探查信息101,所述通用条件是指预定的行驶区域内或者预定的行驶区域内且预定的时间段等。
在步骤S2中,检测装置10生成从获取到的探查信息101中的、由一台车辆100连续发送的探查信息101提取必要的数据并按时系列排列而成的行程信息110。
在步骤S3中,检测装置10根据行程信息110中包含的绝对位置坐标104、时刻信息103以及车速105,辨识发送探查信息101的车辆100的位置。在本实施方式中,检测装置10验证绝对位置坐标104的误差,在误差大的情况下,通过对车速105进行时间积分而得到的车速积分值对绝对位置坐标间的车辆位置进行插值。通过该手法,即使在例如发生多径传播那样的GPS电波的接收状况恶化的情况下,也能够提高车辆位置的精度。
(偏航率的量化)
参照图6,说明对偏航率106进行量化的步骤S4。首先,检测装置10从道路地图信息存储部18中储存的道路地图信息25读出作为运算对象的行驶区域的链路信息,在一个链路L内,以例如每0.1m的固定间隔设定运算对象地点xi(1≤i≤np)。关于该运算对象地点xi,也可以使用与链路L关联起来预先设定的地点。
检测装置10读入1个行程信息110,获取行程信息110中的、在链路上测量到的多个偏航率Y1、Y2…以及表示检测到所述偏航率的位置的车辆位置。车辆位置是在步骤S3中辨识出的车辆的位置。检测装置10对检测到偏航率的车辆位置和运算对象地点xi进行比较,在它们一致时不校正偏航率,在不一致的情况下,通过线性插值、使用样条曲线的插值等公知的方法,对作为实测值的多个偏航率彼此间进行插值。根据进行插值而得到的值,计算与运算对象地点xi关联起来的偏航率yij。标号的“i”表示运算对象地点的识别编号,“j”表示行程的识别编号。运算对象地点xi的个数为“np”,行程Tj的数量为“nj”。
在针对一个行程Tj计算出与运算对象地点xi关联起来的偏航率yij时,针对其他行程Tj也计算偏航率yij。
(偏航率基准值以及绝对偏差的运算)
参照图7,说明运算偏航率基准值以及绝对偏差的步骤S5。检测装置10获取例如与运算对象地点x1关联起来的行程Tj即T1、T2…、Tnt的偏航率y1j也就是y11、y12…、y1np,计算中位值median(y1j)而作为偏航率基准值。关于运算对象地点x2…、xnp,也同样地,如median(y2j)、median(y3j)…、median(ynpj)那样,计算中位值median(yij)。
在本实施方式中,依照下述的式(1),计算使用中位值的中位绝对偏差(MAD:Median Absolute Deviation)。即,计算运算对象地点x1的偏航率y11、y12…、y1np与中位值median(x1)的差分的中位值,计算偏散的中位值。关于运算对象地点x2…、xnp,也同样地,计算MAD(x2)…、MAD(xnp)。通过这样使用中位值,能够提高针对偏移值的鲁棒性。
[式1]
MAD=1.4826·median(|yij-median(yij)|)…(1)
(背离度的运算)
接下来,说明步骤S6的背离度的运算。如上所述,在该步骤S6中,将一个行程信息110作为对象。检测装置10依照下述的式(2),从某个运算对象地点xi的偏航率yij减去所述运算对象地点xi的中位值median(yij),将由此得到的值除以中位绝对偏差MAD(yij),从而作为背离度y^ij。关于一个行程Tj,针对各运算对象地点xi的每一个计算该背离度y^ij。在关于一个行程Tj计算出背离度y^ij后,关于其他行程Tj,也针对各运算对象地点xi的每一个计算背离度y^ij。
[式2]
y ^ i j = y i j - m e d i a n ( y i j ) M A D ( y i j ) ... ( 2 )
(背离度的平滑化)
接下来,说明步骤S7中的背离度的平滑化。运算对象地点xi处的背离度y^ij中还包括方向盘操作的摇晃等所引起的偏移值。因此,检测装置10通过对背离度的推移曲线进行平滑化而去掉偏移值。在本实施方式中,根据行程信息110中包含的车速105,将预定的行驶区域内的行程Tj划分成从出发至停止的区间,在区间内计算移动平均。此时,既可以计算每固定距离的背离度y^ij的移动平均,也可以计算每固定时间的y^ij的移动平均。针对全部行程Tj进行背离度的平滑化。
(异常行驶部位的候选的检测)
接下来,参照图8,说明步骤S8中的异常行驶部位的候选的检测。检测装置10判断在平滑化了的背离度中是否存在超过预先设定的阈值yth的区域。关于阈值yth,经由试验,设定为异常行驶部位的检测中的适合率和再现率适当平衡的值。将背离度y^ij超过阈值yth的区域检测为所述行程Tj的异常行驶部位Zj。每个该行程Tj的异常行驶部位Zj成为用于最终判定异常行驶部位Z的候选。对全部行程Tj检测异常行驶部位Zj。
(异常行驶部位的判定)
接下来,参照图9说明步骤S9中的异常行驶部位的判定。在行程Tj(T1、T2…)彼此中,所检测出的异常行驶部位Zj未必一致或者重复,所以在相互相同的地点或者区间,以异常行驶部位Zj出现的频度比较大为基准,将所述部位确定为异常行驶部位Z。
即使行程Tj彼此间成为候选的异常行驶部位相互重复,异常行驶部位彼此的长度也未必一致,所以根据相互不同的行程Tj的异常行驶部位Zj的间隔、各异常行驶部位Zj的重复区域,最终决定异常行驶部位Z的长度。例如,也可以仅将异常行驶部位Zj重复的地点或者区间确定为异常行驶部位Z。在多个行程Tj的异常行驶部位Zj相互接近的情况下,也可以将包含相互接近的全部异常行驶部位Zj的区间确定为异常行驶部位Z。
在这样检测到异常行驶部位时,检测装置10将与异常行驶部位有关的信息输出到支援信息生成部16。支援信息生成部16根据与异常行驶部位有关的信息生成道路信息、驾驶支援信息,经由通信部15以及网络N发送到车辆100。
接下来,说明依照上述方法检测异常行驶部位时的作用。
在例如大多数的车辆在沿道路行驶的情况下,在该车群中存在变更行车道的车辆时,关于变更了行车道的车辆的行程,变更了行车道的区间被检测为异常行驶部位,但关于其他行程未被检测为异常行驶部位,所以在步骤S9中,最终该区间不被检测为异常行驶部位Z。
另一方面,在车辆例如停在道路的行驶行车道上的情况下,在行驶行车道上行驶的车辆中,既有在停车车辆的跟前50m处躲避的车辆,也有在跟前10m处躲避的车辆。因此,在所述行车道上行驶的车辆的行驶轨迹根据车辆而相互不同,停车位置附近的运算对象地点xi的偏航率yij的偏散变大。在该情况下,也存在所述运算对象地点xi的中位绝对偏差MAD(yij)变大的可能性,但偏离中位值median(yij)的背离度大的偏航率yij也变多。因此,如果适当地设定运算对象地点xi的间隔、进行平滑化的间隔,则在这些运算对象地点xi,关于停车位置附近的运算对象地点xi,偏航率yij为阈值以上的行程Tj变多,能够将停车车辆所存在的地点或者区间检测为异常行驶部位。关于被检测为异常行驶部位的地点或者区间,不仅是存在路上停车的部位,而且还有例如故障车辆、事故车辆、道路的遗撒物所存在的部位、特定的行车道的发生拥堵的部位、交叉路口跟前等。
如以上说明的那样,根据本实施方式,能够得到以下的效果。
(1)根据与每个运算对象地点xi关联起来的多个偏航率yij,运算偏航率的中位值。在检测异常行驶部位时,关于作为异常行驶部位的检测对象的行程信息110,将与运算对象地点xi关联起来的偏航率yij、和所述地点的偏离偏航率的中位值的背离度y^ij大的部位作为异常行驶部位。因此,能够利用车辆标准地搭载的普及性高的传感器,得到与道路的异常行驶部位有关的信息。因此,与异常行驶部位有关的信息成为不依赖于个体差、突发性的要因的信息,所以能够提高一般性。
(2)在计算与运算对象地点xi关联起来的偏航率yij、和所述地点的偏航率yij偏离中位值的背离度y^ij时,使用基于与所述地点xi关联起来的多个偏航率yij得到的中位绝对偏差MAD(yij)。因此,能够以运算对象地点xi的偏离中位值的平均的背离程度为基准,运算检测对象的偏航率yij偏离中位值的背离度y^ij,所以能够提高针对偏移值的背离度y^ij的鲁棒性。
(3)针对背离度y^ij沿着运算对象地点xi的推移计算移动平均值,所以能够抑制方向盘操作的摇晃等所致的转弯值的偏移值对异常行驶部位的检测造成的影响。
(4)在最终判定是否是异常行驶部位Z时,根据检测到异常行驶部位Zj的多个行程信息110,将出现频度大的部位最终判定为异常行驶部位,所以能够得到不受个体限定的信息。
(5)作为用于判断异常行驶部位Z的变量而使用偏航率yij,所以能够提高检测出的异常行驶部位Z、Zj的精度,所述偏航率yij通过与车速等信息并用而能够判断道路的宽度方向的轨迹。车辆100发送包括偏航率106的探查信息101即可,异常行驶部位检测装置10对多个车辆的偏航率进行统计处理即可。例如,相比于使用由车载摄像装置拍摄到的图像数据等来检测异常行驶部位的情况,本实施方式使用了普及性高的传感器,从而能够削减车辆侧所需要的成本,并且能够减轻车载控制部55、异常行驶部位检测装置10所需要的运算负荷等负荷。也无需在车载控制部55等中准备解析图像数据的应用。
(其他实施方式)
上述实施方式还能够以以下那样的方式实施。
·在上述实施方式中,异常行驶部位检测装置10构成为包括CPU、RAM、ROM等,但也可以是这以外的结构。例如,也可以由ASIC(特定用途集成电路)等构成。
·在上述实施方式中,与探查信息存储部17独立地设置了行程信息存储部19,但只要能够将行程信息110临时储存到任意的存储部,则也可以省略行程信息存储部。
·探查信息101中也可以不包括车辆标识符102,而包括例如对每个驾驶员赋予的用户标识符。探查信息101中也可以包括行进方向。在该情况下,能够省略利用异常行驶部位检测装置10进行的车辆100的行进方向的判断。探查信息101中也可以包括与车辆100的行驶行车道有关的信息。在该情况下,异常行驶部位检测装置10能够很容易地针对每个行驶行车道判断是否发生行车道的宽度方向的异常的行驶。
·在上述实施方式中,将异常行驶部位说明为在道路的长度方向上划分出的区间,但在车辆位置的精度高的情况下,也可以辨识道路的宽度方向的位置。
·在上述实施方式中,设为检测装置10提取行程,但检测装置10既可以获取由其他装置生成的行程信息110,也可以获取在车辆侧生成的行程信息110。
·在上述实施方式中,设为使用偏航率的中位值以及中位绝对偏差来计算偏航率的背离度,但也可以省略式(1)的系数“1.4826”。既可以使用中位值以外的平均,也可以使用中位绝对偏差以外的绝对偏差。例如,也可以使用与运算对象地点关联起来的多个偏航率的算术平均值以及平均绝对偏差来计算偏航率的背离度。在使用个数为“n”的偏航率yij的情况下,算术平均值为(yi1+yi2+…+yin)/n。平均绝对偏差是关于个数为“n”的全部偏航率,对从偏航率yij减去偏航率算术平均值Avr而得到的绝对值|yij-Avr|进行加法运算而得到总和,将总和|yi1-Avr|+|yi2-Avr|+…+|yin-Avr|除以“n”而得到的值、或者进一步乘以系数“1.253”而得到的值。
·在上述实施方式中,通过合并每个行程的异常行驶部位,最终判定了异常行驶部位,但也可以根据支援信息的目的,将每个行程的全部异常行驶部位判定为异常行驶部位。
·关于车辆位置的辨识,也可以通过用车速积分值对绝对位置坐标进行插值的方法以外的方法来进行。例如,也可以使用绝对位置坐标和道路地图信息25来进行车辆位置的辨识。或者,也可以使用车速积分值和道路地图信息25来进行车辆位置的辨识。也可以不由检测装置10进行车辆位置的辨识,而在车辆侧辨识车辆位置,异常行驶部位检测装置10使用已辨识的车辆位置。
·在上述实施方式中,等间隔地设定了运算对象地点xi,但也可以以不均等的间隔设定。例如,也可以在交叉路口附近增大运算对象地点的密度,在其以外的区域减小运算对象地点的密度。
·在上述实施方式中,通过移动平均进行了背离度的推移的平滑化,但也可以通过低通滤波器等其他公知的平滑化处理来去掉偏移值。
·如图10所示,根据检测异常行驶部位的目的,也可以使作为用于计算偏航率基准值的统计处理的对象的行程(T1、T2…、Tnt)和作为背离度的计算对象的行程(TA1、TA2…、TAn)为相互不同的行程。在该方案中,也可以例如使收集作为统计处理的对象的行程信息110的期间比收集背离度的计算对象的行程信息110的期间长。由此,根据作为统计处理的对象的行程信息110,能够得到更一般的信息。因此,通过作为统计处理的对象的行程信息110与作为背离度的计算对象的行程信息110的比较,能够在去掉操作摇晃等偏移值的同时检测异常的行驶,所以异常行驶部位的信息的一般性提高。
·在上述实施方式中,针对背离度的推移实施了平滑化处理,但也可以省略所述处理,进行背离度与阈值的比较。在该方案中,在例如检测到背离度的推移中的阈值以上的部位时,在阈值以上的部位的长度小于预定的长度的情况下,也可以不将所述部位检测为异常行驶部位。
·在上述实施方式中,根据上述式(2)计算了背离度。即,将与运算对象地点关联起来的偏航率和与该地点关联起来的多个偏航率的中位值的差分的绝对值,除以与该地点关联起来的多个偏航率的中位绝对偏差。其中,背离度的运算方法能够根据与异常行驶部位有关的信息的用途而变更,也可以通过这以外的手法计算背离度。例如,在求取还包括偏移值的背离度的情况下,也可以简单地将与运算对象地点关联起来的偏航率和算术平均值或者中位值等偏航率基准值的差分作为背离度。
·作为转弯值,探查信息101中除了偏航率106以外还包括横加速度、舵角或者基于陀螺仪的相对方向。也可以包括偏航率106、横加速度、舵角以及基于陀螺仪的相对方向中的多个。进而,也可以包括偏航率106、横加速度、舵角以及基于陀螺仪的相对方向中的至少一个、和表示转向指示器的ON或者OFF的信息。在该情况下,能够检测用户所意图的转弯方向。检测装置10根据探查信息101中包含的转弯值,检测车辆100在道路的宽度方向的举动。
·在上述实施方式中,说明为异常行驶部位检测装置10构成生成道路信息、驾驶支援信息的支援信息生成系统的一部分,但检测装置10也可以被用作通过简单地收集并统计探查信息101来检测异常行驶部位的装置。在该情况下,能够将与异常行驶部位有关的信息用于例如交通检查。
·在上述实施方式中,从搭载于车辆的偏航率传感器获取了转弯值,但也可以从车载传感器以外获取。例如,也可以从被带进车辆的智能手机、平板终端等便携信息终端中设置的陀螺仪、旋转矢量传感器等能够检测车辆的方向的传感器获取转弯值。
·在上述实施方式中,设为异常行驶部位检测装置10构成探查车辆系统的一部分,但异常行驶部位检测装置10也可以搭载于车辆100。在该情况下,异常行驶部位检测装置10收集并积蓄本车辆的行驶历史信息,将在相互不同的定时收集到的多个行驶历史信息作为统计处理的对象,在预定的行驶区域内计算偏航率基准值以及绝对偏差。将统计处理的对象的行驶历史信息之一或者最新的行驶历史信息作为背离度判定对象,计算偏离偏航率基准值的背离度,根据背离度的大小判定异常行驶部位。在该方案中,也可以在检测到异常行驶部位的情况下,在接下来在异常行驶部位行驶时,提醒驾驶员注意。也可以在检测到异常行驶部位的时间点,提醒驾驶员注意。

Claims (6)

1.一种异常行驶部位检测装置,检测异常行驶部位,所述异常行驶部位是发生向道路的宽度方向的异常的行驶的部位,所述异常行驶部位检测装置的特征在于,
所述异常行驶部位检测装置具备:
历史信息获取部,构成为从车辆获取包括转弯值的行驶历史信息,所述转弯值是表示与所述车辆的位置及所述车辆的转弯有关的物理量的值;
转弯值校正部,构成为关于多个所述行驶历史信息,校正这些行驶历史信息中包含的所述转弯值,所述转弯值校正部构成为将所述转弯值校正为与预先设定的每个运算对象地点关联起来的转弯值;
基准转弯值运算部,构成为使用与相互相同的运算对象地点关联起来的多个所述转弯值来运算基准转弯值,所述基准转弯值是平均转弯值;以及
异常行驶部位检测部,构成为关于作为所述异常行驶部位的检测对象的行驶历史信息,运算偏离所述基准转弯值的背离度,所述异常行驶部位检测部构成为将偏离所述基准转弯值的背离度大的部位检测为所述异常行驶部位。
2.根据权利要求1所述的异常行驶部位检测装置,其特征在于,所述异常行驶部位检测部构成为计算与相互相同的运算对象地点关联起来的多个转弯值的绝对偏差,
所述异常行驶部位检测部进而构成为关于作为所述异常行驶部位的检测对象的所述行驶历史信息,将与所述运算对象地点关联起来的所述转弯值与所述基准转弯值的差分除以所述绝对偏差,从而计算所述背离度。
3.根据权利要求1或者2所述的异常行驶部位检测装置,其特征在于,
所述异常行驶部位检测部构成为计算所述转弯值偏离所述基准转弯值的背离度的沿着所述运算对象地点的推移,
所述异常行驶部位检测部进而构成为针对计算的所述背离度的推移实施平滑化处理。
4.根据权利要求1~3中的任意一项所述的异常行驶部位检测装置,其特征在于,
所述异常行驶部位检测部构成为针对作为异常行驶部位的检测对象的每个行驶历史信息检测异常行驶部位,
所述异常行驶部位检测部进而构成为在基于检测到所述异常行驶部位的多个行驶历史信息的相互相同的部位处的异常行驶部位的出现频度大的情况下,将所述部位最终判定为所述异常行驶部位。
5.根据权利要求1~4中的任意一项所述的异常行驶部位检测装置,其特征在于,
所述行驶历史信息中包含的转弯值是车辆的偏航率。
6.一种异常行驶部位检测方法,检测异常行驶部位,所述异常行驶部位是发生向道路的宽度方向的异常的行驶的部位,所述异常行驶部位检测方法的特征在于,包括:
从车辆获取包括转弯值的行驶历史信息,所述转弯值是表示与所述车辆的位置及所述车辆的转弯有关的物理量的值;
关于多个行驶历史信息,校正这些行驶历史信息中包含的转弯值,所述校正将所述转弯值校正为与预先设定的每个运算对象地点关联起来的转弯值;
使用与相互相同的运算对象地点关联起来的多个所述转弯值来运算基准转弯值,所述基准转弯值是平均转弯值;
关于作为异常行驶部位的检测对象的行驶历史信息,运算偏离所述基准转弯值的背离度;以及
将偏离所述基准转弯值的所述背离度大的部位检测为所述异常行驶部位。
CN201480073508.2A 2014-01-21 2014-12-22 异常行驶部位检测装置以及异常行驶部位检测方法 Expired - Fee Related CN105934786B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2014008450A JP5929936B2 (ja) 2014-01-21 2014-01-21 特異走行箇所検出装置及び特異走行箇所検出方法
JP2014-008450 2014-01-21
PCT/JP2014/083888 WO2015111344A1 (ja) 2014-01-21 2014-12-22 特異走行箇所検出装置および特異走行箇所検出方法

Publications (2)

Publication Number Publication Date
CN105934786A true CN105934786A (zh) 2016-09-07
CN105934786B CN105934786B (zh) 2018-01-16

Family

ID=53681161

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201480073508.2A Expired - Fee Related CN105934786B (zh) 2014-01-21 2014-12-22 异常行驶部位检测装置以及异常行驶部位检测方法

Country Status (5)

Country Link
US (1) US10089867B2 (zh)
EP (1) EP3098797B1 (zh)
JP (1) JP5929936B2 (zh)
CN (1) CN105934786B (zh)
WO (1) WO2015111344A1 (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106952473A (zh) * 2017-04-01 2017-07-14 深圳市元征科技股份有限公司 道路故障检测方法及装置
CN108242145A (zh) * 2016-12-26 2018-07-03 高德软件有限公司 异常轨迹点检测方法和装置
CN108871357A (zh) * 2018-04-20 2018-11-23 长安大学 一种在电子地图上显示拥堵路段的事故车道的方法
CN111157018A (zh) * 2019-12-31 2020-05-15 华为技术有限公司 一种检测方法、检测装置及自动驾驶车辆
CN111325986A (zh) * 2019-04-04 2020-06-23 北京嘀嘀无限科技发展有限公司 异常停车监测方法、装置、电子设备以及存储介质
CN113744550A (zh) * 2020-05-15 2021-12-03 丰田自动车株式会社 信息处理装置以及信息处理系统

Families Citing this family (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
IL236234A0 (en) * 2014-12-14 2015-03-31 Google Llc Systems and methods for creating travel plans using location information
JP6410670B2 (ja) * 2015-06-17 2018-10-24 日本電信電話株式会社 通行難易度推定装置、方法およびプログラム
JP6657646B2 (ja) * 2015-08-06 2020-03-04 オムロン株式会社 障害物検知装置、障害物検知方法、および障害物検知プログラム
JP6699230B2 (ja) 2016-02-25 2020-05-27 住友電気工業株式会社 道路異常警告システム及び車載機
JP6896374B2 (ja) * 2016-04-22 2021-06-30 株式会社東芝 プローブ情報処理サーバ、およびプローブ情報処理方法
JP6776694B2 (ja) * 2016-07-25 2020-10-28 日産自動車株式会社 情報提供方法および情報提供装置
JP6718811B2 (ja) * 2016-12-28 2020-07-08 本田技研工業株式会社 情報処理システム、および情報処理方法
US20190362433A1 (en) * 2017-01-26 2019-11-28 Sony Corporation Information processing apparatus and method, program, and information processing system
WO2018225347A1 (ja) * 2017-06-09 2018-12-13 ボッシュエンジニアリング株式会社 走行障害検出装置及び車両ナビゲーションシステム
US10482402B2 (en) * 2017-08-22 2019-11-19 Caterpillar Inc. Method and system for providing notifications related to service information of machine
JP6446116B1 (ja) * 2017-11-27 2018-12-26 トーヨーカネツソリューションズ株式会社 搬送システム検査装置
US10448356B1 (en) * 2018-03-30 2019-10-15 AVAST Software s.r.o. Mobile device location anomaly detection based on non-location information
JP2019185109A (ja) * 2018-04-02 2019-10-24 パイオニア株式会社 検知装置、検知方法及びプログラム
JP7129282B2 (ja) * 2018-08-29 2022-09-01 日産自動車株式会社 障害物区間推定方法及び障害物区間推定装置
JP2020154623A (ja) * 2019-03-19 2020-09-24 株式会社Subaru 交通制御システム
CN112820111B (zh) * 2021-02-03 2022-11-01 腾讯科技(深圳)有限公司 一种异常车辆识别方法、装置和电子设备
US20230104188A1 (en) * 2021-09-28 2023-04-06 Here Global B.V. Method, apparatus, and system for calibrating vehicle motion data based on mobile device sensor data
CN118413823B (zh) * 2024-07-03 2024-09-03 大连海大赢海科技有限公司 一种用于船端与岸基的多源异构数据传输方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1572597A (zh) * 2003-05-30 2005-02-02 丰田自动车株式会社 冲撞预测装置
JP2006313519A (ja) * 2005-04-04 2006-11-16 Sumitomo Electric Ind Ltd 障害物検出センター装置、障害物検出システム及び障害物検出方法
CN101287634A (zh) * 2005-10-13 2008-10-15 日产自动车株式会社 车辆驾驶辅助系统
CN102132335A (zh) * 2008-09-19 2011-07-20 株式会社日立制作所 行驶环境识别装置
CN102490721A (zh) * 2006-04-25 2012-06-13 丰田自动车株式会社 停车辅助装置

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5850193A (en) * 1995-03-30 1998-12-15 Sumitomo Electric Industries, Ltd. Apparatus for assisting driver in carefully driving
JP4367174B2 (ja) * 2004-02-25 2009-11-18 株式会社デンソー 車載送信装置および障害物検知システム
JP4945222B2 (ja) * 2006-11-28 2012-06-06 日立オートモティブシステムズ株式会社 突発事象解消判定システム
JP4845783B2 (ja) * 2007-03-16 2011-12-28 パイオニア株式会社 情報処理方法、車載装置および情報配信装置
JP4985150B2 (ja) 2007-06-29 2012-07-25 日産自動車株式会社 車両用運転支援装置および車両用運転支援装置を備える車両
JP4946903B2 (ja) * 2008-02-12 2012-06-06 トヨタ自動車株式会社 車両用背景画像収集装置
JP5774847B2 (ja) * 2010-12-20 2015-09-09 デジスパイス株式会社 車両走行再現評価装置
JP5389864B2 (ja) 2011-06-17 2014-01-15 クラリオン株式会社 車線逸脱警報装置
JP5905680B2 (ja) * 2011-08-12 2016-04-20 株式会社ナビタイムジャパン 探索結果生成システム、サーバ装置、端末装置、探索結果生成方法、および、プログラム
JP5741954B2 (ja) * 2012-02-17 2015-07-01 アイシン・エィ・ダブリュ株式会社 交通情報配信システム、ナビゲーションシステム、交通情報配信プログラム、及び交通情報配信方法
US9098956B2 (en) * 2013-09-26 2015-08-04 Lytx, Inc. Dynamic uploading protocol

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1572597A (zh) * 2003-05-30 2005-02-02 丰田自动车株式会社 冲撞预测装置
JP2006313519A (ja) * 2005-04-04 2006-11-16 Sumitomo Electric Ind Ltd 障害物検出センター装置、障害物検出システム及び障害物検出方法
CN101287634A (zh) * 2005-10-13 2008-10-15 日产自动车株式会社 车辆驾驶辅助系统
CN102490721A (zh) * 2006-04-25 2012-06-13 丰田自动车株式会社 停车辅助装置
CN102132335A (zh) * 2008-09-19 2011-07-20 株式会社日立制作所 行驶环境识别装置

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108242145A (zh) * 2016-12-26 2018-07-03 高德软件有限公司 异常轨迹点检测方法和装置
CN106952473A (zh) * 2017-04-01 2017-07-14 深圳市元征科技股份有限公司 道路故障检测方法及装置
CN108871357A (zh) * 2018-04-20 2018-11-23 长安大学 一种在电子地图上显示拥堵路段的事故车道的方法
CN108871357B (zh) * 2018-04-20 2021-09-07 长安大学 一种在电子地图上显示拥堵路段的事故车道的方法
CN111325986A (zh) * 2019-04-04 2020-06-23 北京嘀嘀无限科技发展有限公司 异常停车监测方法、装置、电子设备以及存储介质
CN111157018A (zh) * 2019-12-31 2020-05-15 华为技术有限公司 一种检测方法、检测装置及自动驾驶车辆
CN111157018B (zh) * 2019-12-31 2021-08-31 华为技术有限公司 一种检测方法、检测装置及自动驾驶车辆
CN113744550A (zh) * 2020-05-15 2021-12-03 丰田自动车株式会社 信息处理装置以及信息处理系统

Also Published As

Publication number Publication date
JP2015138316A (ja) 2015-07-30
EP3098797B1 (en) 2021-01-20
EP3098797A1 (en) 2016-11-30
US10089867B2 (en) 2018-10-02
CN105934786B (zh) 2018-01-16
JP5929936B2 (ja) 2016-06-08
WO2015111344A1 (ja) 2015-07-30
EP3098797A4 (en) 2017-03-01
US20160379484A1 (en) 2016-12-29

Similar Documents

Publication Publication Date Title
CN105934786A (zh) 异常行驶部位检测装置以及异常行驶部位检测方法
US10704920B2 (en) Traffic lane guidance system for vehicle and traffic lane guidance method for vehicle
CN104933856B (zh) 道路路况实时评估系统及方法
CN104567885B (zh) 导航装置以及高架上下道路判定方法
CN103164975A (zh) 警告系统、车载设备和服务器
CN104471626A (zh) 车道脱离判定装置、车道脱离警报装置以及使用它们的车辆控制系统
JP3970473B2 (ja) 監視手段付きgps装置
JP2009140008A (ja) 危険走行情報提供装置、危険走行判定プログラム及び危険走行判定方法
CN105612569A (zh) 驻车车辆检测装置、车辆管理系统以及控制方法
CN104192063B (zh) 车辆安全行驶警示系统及相应的警示方法
CN204055561U (zh) 车辆安全行驶警示系统
CN101958043A (zh) 路段旅行时间计算装置及其方法
CN109849785A (zh) 一种基于车联网的换道判定系统及方法
CN109387218B (zh) 车载设备及道路维护辅助管理系统
CN111183464B (zh) 基于车辆轨迹数据估计信号交叉路口的饱和流量的系统和方法
CN106846369A (zh) 基于双目视觉的车辆转向状态判别方法及装置
CN114599532A (zh) 载荷检测设备
US20230051377A1 (en) Mobility movemennt information acquiring method and mobility movement information acquiring apparatus
US11919535B2 (en) On-board detection and logging of driver smartphone use
CN103241243B (zh) 一种确定车辆状态的方法和装置
JP2016066144A (ja) 路面情報収集装置および路面情報解析システム
CN113793159B (zh) 逃票检测方法、装置及电子设备
US20230236020A1 (en) System and Method for Map Matching GNSS Positions of a Vehicle
JP5396988B2 (ja) 道路交通上の人件費算出システム、演算装置、及び、道路交通上のシステム導入効果評価方法
JP5162962B2 (ja) 物体検出装置

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

CF01 Termination of patent right due to non-payment of annual fee