CN102339143B - 轨迹平滑化的方法 - Google Patents

轨迹平滑化的方法 Download PDF

Info

Publication number
CN102339143B
CN102339143B CN201010228914.5A CN201010228914A CN102339143B CN 102339143 B CN102339143 B CN 102339143B CN 201010228914 A CN201010228914 A CN 201010228914A CN 102339143 B CN102339143 B CN 102339143B
Authority
CN
China
Prior art keywords
coordinate
mean number
speed
value
moving average
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
CN201010228914.5A
Other languages
English (en)
Other versions
CN102339143A (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.)
QIGAO ELECTRONIC CO Ltd
Original Assignee
QIGAO ELECTRONIC CO 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 QIGAO ELECTRONIC CO Ltd filed Critical QIGAO ELECTRONIC CO Ltd
Priority to CN201010228914.5A priority Critical patent/CN102339143B/zh
Publication of CN102339143A publication Critical patent/CN102339143A/zh
Application granted granted Critical
Publication of CN102339143B publication Critical patent/CN102339143B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Complex Calculations (AREA)

Abstract

本发明提出一种轨迹平滑化的方法,包含:依据一物体的移动轨迹,产生多个轨迹点并取得该多个轨迹点的坐标;依据移动平均法,计算该多个轨迹点的移动平均坐标,其中该移动平均坐标以一平均个数作为分母进行平均计算;依据该物体移动的速度,适应性地调整该平均法的平均个数值,其中,当该物体移动的速度较快时,采用较低的平均个数值,当物体移动速度较慢时,采用较高的平均个数值。

Description

轨迹平滑化的方法
技术领域
本发明涉及一种轨迹平滑化的方法,特别是指一种可动态调整计算方式的轨迹平滑化方法。
背景技术
现有技术中,对于多个轨迹点平滑化的方法,是利用平均法来达成,以计算该多个轨迹点的移动平均坐标,其中该移动平均坐标以一平均个数(n)作为分母进行平均计算。请参考图1为平均个数的值与曲线平滑度的关系,其中C1、C2、C3分别代表利用不同平均个数计算移动平均坐标的曲线。曲线C1平均个数n为2;曲线C2的平均个数n为8;曲线C3的平均个数n为16。由图中可知,当平均个数n的值越高时,曲线的平滑效果越好。但平均个数n的值越高时,时间的延迟越久;请再参见图2,更进一步说明平均个数n与时间延迟的关系。曲线C4和曲线C5为利用不同平均个数n计算移动平均坐标的曲线,其中曲线C4的平均个数为4;曲线C5的平均个数为16。由图中可知,当平均个数n的值越高,平滑化的轨迹跟上实际的轨迹所需的时间越久。
综合以上所述,现有技术中,不管是以较多或较少的平均个数n作为移动中轨迹平滑化的基准,皆有其优缺点,无法让系统达到一个可以两边兼顾的状态。
发明内容
本发明目的之一在于克服现有技术的不足与缺陷,提出一种轨迹平滑化的方法。
为达上述目的,就其中一个观点而言,本发明提供了轨迹平滑化的方法,用以平滑化一物体的移动轨迹,包含:依据该物体的移动轨迹,产生多个轨迹点并取得该多个轨迹点的坐标;依据移动平均法,计算该多个轨迹点的移动平均坐标,其中该移动平均坐标以一平均个数(n)作为分母进行平均计算;以及依据该物体移动的速度,适应性地调整该平均个数值,其中,当该物体移动的速度较快时,采用较低的平均个数值,当物体移动速度较慢时,采用较高的平均个数值。
在其中一种较佳实施形式中,该移动平均法依以下公式产生移动平均坐标avgT
sumT=sum(T-1)-avg(T-1)+new;
avgT=sumT/n,
其中avgT和avg(T-1)分别代表目前时间T所产生出的移动平均坐标和前一时间T-1所产生出的移动平均坐标,sumT和sum(T-1)分别代表目前时间T的坐标值总和和前一时间T-1的坐标值总和,new代表目前时间T所得的轨迹点坐标,n代表平均个数。
在其中另一种较佳实施形式中,该依据该物体移动的速度,适应性地调整该平均个数值的步骤包括:提供至少一个速度临界值;判断该物体移动速度是否高于该速度临界值;以及根据判断结果,适应性地调整该平均个数值。
在其中另一种较佳实施形式中,该依据该物体移动的速度,适应性地调整该平均个数的数值的步骤包括:提供至少一组速度迟滞临界值,此速度迟滞临界值包含一高限值与一低限值;当该物体移动速度增加以致高于该高限值时,采用较低的平均个数值;以及当该物体移动速度降低以致低于该低限值时,采用较高的平均个数值。
在其中再一种较佳实施形式中,所取得的轨迹点坐标值为绝对坐标值。在其中再一种较佳实施形式中,所取得的轨迹点坐标值为相对坐标值,且方法更包含:将该相对坐标值转换为绝对坐标值。所述物体例如为遥控器或鼠标。
下面通过具体实施例详加说明,当更容易了解本发明的目的、技术内容、特点及其所达成的功效。
附图说明
图1标出多个轨迹点平滑化的方法中,不同平均个数的曲线的比较;
图2标不同平均个数n所计算出的平滑化轨迹曲线与时间延迟的关系图;
图3-5标出数个实施例,揭露依据物体移动速度决定平均个数n的方法;
图6标出依据物体移动速度对平均个数n作阶级调整的流程图;
图7标出阶级变换设有迟滞区间时的流程图;
图8标出阶级变换设有迟滞区间时,物体移动速度与平均个数n的关系。
图中符号说明
C1,C2,C3,C4,C5          曲线
n,n0,n1,n2               平均个数
stg0,stg1,stg2            阶级
S101~S107,S201~S207      步骤
t                           时间
VT1,VT2                    速度临界值
VT1_L,VT1_H,VT2_L,VT2_H  迟滞速度临界值
具体实施方式
如前所述,现有技术中,不管是以较多或较少的平均个数n作为轨迹平滑化的基准,皆有其缺点。有鉴于此,本发明即针对上述现有技术的不足,提出一种轨迹平滑化的方法,依据物体移动的速度,调整计算移动平均坐标的平均个数n。详言之,当物体移动速度较慢时,依据物体移动轨迹所产生的多个轨迹点的坐标容易受到噪声的干扰而导致轨迹有不平滑或锯齿状的情形发生。相对地,当物体移动较快时,由于位移量大,噪声所产生的效应较低,不平滑的情况较不显著。例如,当噪声的值为±1单位时,如果单位时间移动的距离为10单位,噪声的比例为10%;而当单位时间移动的距离为2单位时,噪声的比例则为50%。另一方面,当物体移动较快时,移动平均计算所造成的时间延迟效应较大;当物体移动速度较慢时,移动平均计算所造成的时间延迟效应则较不显著。因此,本发明利用实时侦测物体位移速度的快慢,调整计算移动平坐标的平均个数n,其中当物体移动的速度较快时,采用较少的平均个数n计算多个轨迹点的移动平均坐标avg,以提高系统运作的效率;当物体移动的速度较快时,采用较多的平均个数n计算多个轨迹点的移动平均坐标avg,以过滤移动轨迹的噪声,并提高物体移动轨迹的平滑度。
计算移动平均坐标avg的方法有很多,可依据使用者的需要代入不同的公式运算,产生移动平均坐标avg。举例而言,其中一种计算移动平均坐标avg的公式为:
sumT=sum(T-1)-avg(T-1)+new;
avgT=sum/n,
其中avgT和avg(T-1)分别代表目前时间T所产生出的移动平均坐标avg和前一时间T-1所产生出的移动平均坐标avg,sumT和sum(T-1)分别代表目前时间T的坐标值总和和前一时间T-1的坐标值总和,new代表目前时间T所得的轨迹点坐标,n代表平均个数。
根据本发明,依据物体移动速度调整平均个数n有多种实施方式,其中物体移动速度与平均个数n的关系例如可为:线性调整(如n=[-av+b],其中a,b为常数,n可取整数,请参见图3)、反比调整(如n=[a/v+b],其中a,b为常数,n可取整数,请参见图4)、或阶级调整(请参见图5)。以阶级调整的方法为例,请参见图6为阶级调整的流程图,并对照图5,其中n0>n1>n2,VT2>VT1:
在步骤S100中,系统开始初始化;接着,进入步骤S101。
在步骤S101中,以阶级stg0的预设平均个数值(n=n0),进行移动平均坐标avg运算;接着,进入步骤S102。
在步骤S102中,侦测物体的速度,如果物体的移动速度未超过第一速度临界值VT1,平均个数维持在起始阶级stg0的值n=n0,回到步骤S101;当物体移动速度增加到超过第一速度临界值VT1时,进入步骤S103。
在步骤S103中,以阶级stg1的预设平均个数值(n=n1),进行移动平均坐标avg运算;接着,进入步骤S104与步骤S107。
在步骤S104与步骤S107中,侦测物体的速度,如果物体的移动速度未低于第一速度临界值VT1、或未超过第二速度临界值VT2,平均个数维持在阶级stg1的值n=n1,回到步骤S103。当物体移动速度减低到低于第一速度临界值VT1时,回到步骤S101;当物体移动速度增加到超过第二速度临界值VT2时,进入步骤S105。其中,步骤S104与步骤S107可平行或以任意次序先后执行。
在步骤S105中,以阶级stg2的预设平均个数值(n=n2),进行移动平均坐标avg运算;接着,进入步骤S106。
在步骤S106中,侦测物体的速度,如果物体的移动速度持续超过第二速度临界值VT2,平均个数维持在阶级stg2的值n=n2,回到步骤S105;当物体移动速度减小到低于第二速度临界值VT2时,回到步骤S103。
上述阶级调整的优点为,平均个数n的值不需因物体移动速度的微幅改变而变动,有助于维持系统的稳定性并降低电路的复杂度。除上述外,阶级调整的方法尚有许多安排方式,本领域技术人员于了解本发明的内容后,应可思及各种变化,皆应包含在本发明的专利范围内。例如,在各阶变换处可安排一个迟滞区间,举例而言请参见图7并对照图8,其中n0>n1>n2,VT2_H>VT2_L>VT1_H>VT1_L:
图7中在步骤S200时,系统开始初始化;接着,进入步骤S201。
在步骤S201中,以阶级stg0的预设平均个数值(n=n0),进行移动平均坐标avg运算;接着,进入步骤S202。
在步骤S202中,侦测物体的速度,如果物体的移动速度未超过第一迟滞临界值的高限值VT1_H,平均个数维持在起始阶级stg0的值n0,回到步骤S201;当物体移动速度增加到超过第一迟滞临界值的高限值VT1_H时,进入步骤S203。
在步骤S203中,以阶级stg1的预设平均个数值(n=n1),进行移动平均坐标avg运算;接着,进入步骤S204与步骤S207。
在步骤S204与步骤S207中,侦测物体的速度;如果物体的移动速度未低于第一迟滞临界值的低限值VT1_L、或未超过第二迟滞临界值的高限值VT2_H,平均个数维持在阶级stg1的值n=n1,回到步骤S203。当物体移动速度减低到低于第一迟滞临界值的低限值VT1_L时,回到步骤S201;当物体移动速度增加到超过第二迟滞临界值的高限值VT2_H时,进入步骤S205。其中,步骤S204与步骤S207可平行或以任意次序先后执行。
在步骤S205中,以阶级stg2的预设平均个数值(n=n2),进行移动平均坐标avg运算;接着,进入步骤S206。
在步骤S206中,侦测物体的速度,如果物体的移动速度持续超过第二迟滞临界值的低限值VT2_L时,平均个数维持在阶级stg2的值n=n2,回到步骤S205;当物体移动速度减小到低于第二迟滞临界值的低限值VT2_L时,回到步骤S203。
上述图6或7的流程中,系统于起始后并不必须由阶级Stg0开始,而可由Stg0,Stg1或Stg2中任一阶级开始;平均个数n的起始默认值不必须为n0(最大值,对应于最低速度),而当然也可以从最小默认值n=n2(对应于最高速度)开始调整,或从任一中间步骤(如n=n1值)开始调整,亦即系统于起始后可衔接图6或7中任一步骤S101~S107或S201~S207。
上述实施例中,依据物体移动轨迹所产生的多个轨迹点的坐标可为绝对坐标,亦可为相对坐标;只是,当多个轨迹点的坐标为相对坐标时,计算移动平均坐标avg前,宜先转换相对坐标为绝对坐标,以利系统运算。所述物体可为在绝对坐标下移动的物体,亦可为在相对坐标下移动的物体。所述物体例如可为遥控器或鼠标等。
此外,还需说明:本发明以上所述的坐标计算,可适用于任意维度的坐标系统。虽然图1是以x-y二维坐标系统来说明,但本发明不限于应用在二维坐标系统,而可应用于一维至无限多维度。如维度等于或超过二维,则在运算时,每一维度的坐标可以分开平行运算、或依序分时运算、或两者综合,可视硬件成本与运算效率来做取舍。
以上已针对较佳实施例来说明本发明,只是以上所述,仅为使本领域技术人员易于了解本发明的内容,并非用来限定本发明的权利范围。凡此种种,在本发明的相同精神下,本领域技术人员可以思及各种等效变化,均应包含在本发明的范围之内。

Claims (6)

1.一种轨迹平滑化的方法,用以平滑化一物体的移动轨迹,其特征在于,包含:
依据该物体的移动轨迹,产生多个轨迹点并取得该多个轨迹点的坐标;
依据移动平均法,计算该多个轨迹点的移动平均坐标,其中该移动平均坐标是以一平均个数作为分母进行平均计算;以及
依据该物体移动的动态速度,适应性地调整该平均个数值,在此步骤中,当该物体移动的动态速度相对较快时,采用相对较低的平均个数值,当物体移动的动态速度相对较慢时,采用相对较高的平均个数值,且该步骤包括:
提供至少一个速度临界值;
判断该物体移动速度是否高于该速度临界值;以及
根据判断结果,适应性地调整该平均个数值。
2.如权利要求1所述的轨迹平滑化的方法,其中,该移动平均法依以下公式产生移动平均坐标avgT:
sumT=sum(T-1)-avg(T-1)+new
avgT=sumT/n
其中avgT和avg(T-1)分别代表目前时间T所产生出的移动平均坐标和前一时间T-1所产生出的移动平均坐标,sumT和sum(T-1)分别代表目前时间T的坐标值总和和前一时间T-1的坐标值总和,new代表目前时间T所得的轨迹点坐标,n代表平均个数。
3.如权利要求1所述的轨迹平滑化的方法,其中,该依据该物体移动的速度,适应性地调整该平均个数的值的步骤包括:
提供至少一组速度迟滞临界值,此速度迟滞临界值包含一高限值与一低限值;
当该物体移动速度增加以致高于该高限值时,采用较低的平均个数值;以及
当该物体移动速度降低以致低于该低限值时,采用较高的平均个数值。
4.如权利要求1所述的轨迹平滑化的方法,其中,所取得的轨迹点坐标值为绝对坐标值。
5.如权利要求1所述的轨迹平滑化的方法,其中,所取得的轨迹点坐标值为相对坐标值,且方法还包含:将该相对坐标值转换为绝对坐标值。
6.如权利要求1所述的轨迹平滑化的方法,其中,该物体为一遥控器或鼠标。
CN201010228914.5A 2010-07-14 2010-07-14 轨迹平滑化的方法 Active CN102339143B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201010228914.5A CN102339143B (zh) 2010-07-14 2010-07-14 轨迹平滑化的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201010228914.5A CN102339143B (zh) 2010-07-14 2010-07-14 轨迹平滑化的方法

Publications (2)

Publication Number Publication Date
CN102339143A CN102339143A (zh) 2012-02-01
CN102339143B true CN102339143B (zh) 2015-01-14

Family

ID=45514908

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201010228914.5A Active CN102339143B (zh) 2010-07-14 2010-07-14 轨迹平滑化的方法

Country Status (1)

Country Link
CN (1) CN102339143B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106289320B (zh) * 2016-08-01 2019-03-05 成都天合世纪科技有限责任公司 一种误差校准方法及装置
TW201926149A (zh) * 2017-11-22 2019-07-01 財團法人資訊工業策進會 工段歷程追蹤系統及工段歷程追蹤方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5661502A (en) * 1996-02-16 1997-08-26 Ast Research, Inc. Self-adjusting digital filter for smoothing computer mouse movement
CN1652143A (zh) * 2004-02-06 2005-08-10 原相科技股份有限公司 指针装置及轨迹平滑方法
CN101739158A (zh) * 2008-11-24 2010-06-16 原相科技股份有限公司 轨迹平滑化方法及使用所述方法的导航装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5661502A (en) * 1996-02-16 1997-08-26 Ast Research, Inc. Self-adjusting digital filter for smoothing computer mouse movement
CN1652143A (zh) * 2004-02-06 2005-08-10 原相科技股份有限公司 指针装置及轨迹平滑方法
CN101739158A (zh) * 2008-11-24 2010-06-16 原相科技股份有限公司 轨迹平滑化方法及使用所述方法的导航装置

Also Published As

Publication number Publication date
CN102339143A (zh) 2012-02-01

Similar Documents

Publication Publication Date Title
CN108845973B (zh) 一种基于改进Quinn算法的多普勒频率估计方法
CN102436438A (zh) 基于gpu的稀疏矩阵数据存储方法
EP2348685A3 (en) Method and apparatus for reducing peak to average power ratio using peak windowing
CN102339143B (zh) 轨迹平滑化的方法
CN107274422A (zh) 一种基于法线信息和k邻域搜索结合的点云边缘提取方法
CN103675447B (zh) 一种电气化铁路的高精度实时谐波分析方法
EP2811417A3 (en) Systems and methods for synchronizing geographic information system (GIS) network models
CN110703684B (zh) 一种端点速度不限的轨迹规划方法及装置
CN103337084B (zh) 一种基于人工地物特征的地物图斑自动生成方法
WO2014121573A1 (zh) 基于非局部均值的susan边缘检测方法及系统
CN106202002A (zh) 一种用于检测水文系列参数是否变异的方法
TWI448947B (zh) Capacitive touchpad detection method
CN103438827B (zh) 螺丝的螺纹区域峰值线的检测方法
CN103473607B (zh) 按轨迹特征优化和匹配外推模型的风电超短期预测方法
CN103488844B (zh) 一种基于厂站地理信息的电网网架图自动布线方法
CN104376578A (zh) 一种应用于直录播系统中的运动目标检测处理方法及设备
CN107885696A (zh) 一种利用观测序列相似性实现缺失数据修复的方法
CN112600210A (zh) 一种应用于低电压治理装置的电压自适应调整方法
CN1399443A (zh) 网络节点拓扑结构图示方法
EP2533166A3 (en) Automatic extremum detection on a surface mesh of a component
CN104238457B (zh) 一种计算复杂度自适应的nurbs曲线插补方法
CN102801983A (zh) 一种基于dct的去噪方法及装置
CN103778595B (zh) 一种图像缩小处理的方法及装置
CN106648395B (zh) 一种书写笔迹分段平滑方法和系统
CN106655799B (zh) Buck-Boost矩阵变换器稳定性判定方法及装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant