CN114076946A - 一种运动估计方法及装置 - Google Patents

一种运动估计方法及装置 Download PDF

Info

Publication number
CN114076946A
CN114076946A CN202010831971.6A CN202010831971A CN114076946A CN 114076946 A CN114076946 A CN 114076946A CN 202010831971 A CN202010831971 A CN 202010831971A CN 114076946 A CN114076946 A CN 114076946A
Authority
CN
China
Prior art keywords
velocity vector
sensor
carrier
estimate
ith
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.)
Pending
Application number
CN202010831971.6A
Other languages
English (en)
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.)
Huawei Technologies Co Ltd
Original Assignee
Huawei Technologies 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 Huawei Technologies Co Ltd filed Critical Huawei Technologies Co Ltd
Priority to CN202010831971.6A priority Critical patent/CN114076946A/zh
Priority to PCT/CN2021/108627 priority patent/WO2022037370A1/zh
Publication of CN114076946A publication Critical patent/CN114076946A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/93Radar or analogous systems specially adapted for specific applications for anti-collision purposes
    • G01S13/931Radar or analogous systems specially adapted for specific applications for anti-collision purposes of land vehicles
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/10Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
    • G01C21/12Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
    • G01C21/16Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/10Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
    • G01C21/12Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
    • G01C21/16Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation
    • G01C21/165Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation combined with non-inertial navigation instruments
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/93Sonar systems specially adapted for specific applications for anti-collision purposes
    • G01S15/931Sonar systems specially adapted for specific applications for anti-collision purposes of land vehicles
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42Determining position
    • G01S19/45Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement
    • G01S19/46Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement the supplementary measurement being of a radio-wave signal type
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42Determining position
    • G01S19/45Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement
    • G01S19/47Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement the supplementary measurement being an inertial measurement, e.g. tightly coupled inertial

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Electromagnetism (AREA)
  • Automation & Control Theory (AREA)
  • Acoustics & Sound (AREA)
  • Navigation (AREA)
  • Traffic Control Systems (AREA)

Abstract

本申请公开了一种运动估计方法及装置,属于传感器技术领域,可用于辅助驾驶和自动驾驶。其中方法包括:获取M个第一传感器的转动角速度矢量估计值和N个第二传感器的瞬时速度矢量估计值,M≥1,N≥1,根据N个第二传感器的瞬时速度矢量估计值和N个第二传感器所在的载体的第一转动角速度矢量估计值,确定载体的第一平动速度矢量估计值,其中,第一转动角速度矢量估计值是根据M个第一传感器的转动角速度矢量估计值确定的。该方法可用于传感器感知的数据处理过程中,该方法用于准确估计传感器所在载体的运动,可以应用于车联网,如车辆外联V2X、车间通信长期演进技术LTE‑V、车辆‑车辆V2V等。

Description

一种运动估计方法及装置
技术领域
本申请涉及传感器技术领域,尤其涉及一种运动估计方法及装置。
背景技术
先进辅助驾驶系统(advanced driver assistant system,ADAS)或者自动驾驶(autonomous driving,AD)系统会配置多种传感器,例如毫米波雷达、激光雷达、超声波传感器如声呐、视觉传感器如相机或者摄像头等传感器,用于感知周边环境信息,周边环境信息包括运动目标和静止目标。对于运动目标和静止目标,通常采用不同的方法分析处理,示例性的,对运动目标(比如车辆、行人)进行分类、识别和跟踪,对静止目标(比如障碍物、护栏、路沿)分类和识别。通过上述方式,可以为自动驾驶提供额外信息,如规避障碍物、提供可行驶区域等。
传感器通常可以安装于载体,传感器跟随传感器所在载体运动。一方面,传感器所在载体的运动导致运动目标和静止目标无法独立分析,因此,需要估计传感器所在载体的运动,从而实现运动目标和静止目标的分离。另一方面,运动目标的跟踪通常基于运动模型,如常速度(constant velocity,CV)/常加速度(constant acceleration,CA)/匀速圆周运动(coordinated turn,CT)等模型,且模型通常假定相对地面或者大地坐标系,传感器所在载体的运动将导致上述模型失效或者跟踪性能下降,因此,需要对传感器所在载体的运动进行补偿。
综上,准确估计传感器所在载体的运动,是本领域人员正在解决的技术问题。
发明内容
本申请提供一种运动估计方法及装置,用于准确估计传感器所在载体的运动。
第一方面,本申请提供一种运动估计方法,该方法包括:
获取M个第一传感器的转动角速度矢量估计值和N个第二传感器的瞬时速度矢量估计值;其中,M≥1,N≥1;
根据所述N个第二传感器的瞬时速度矢量估计值和所述N个第二传感器所在的载体的第一转动角速度矢量估计值,确定所述载体的第一平动速度矢量估计值,其中,所述第一转动角速度矢量估计值是根据所述M个第一传感器的转动角速度矢量估计值确定的。
上述技术方案中,由于第一传感器的转动角速度矢量估计值和第二传感器的瞬时速度矢量估计值较为准确,根据至少一个第一传感器的转动角速度矢量估计值,确定载体的第一转动角速度矢量估计值,再结合第二传感器的瞬时速度矢量估计值,确定载体的第一平动速度矢量估计值,有助于得到较为准确的传感器所在载体的运动。
在一种可能的实现方式中,所述第一平动速度矢量估计值基于如下关系式确定:
Figure BDA0002638325200000011
其中,
Figure BDA0002638325200000012
为第一平动速度矢量估计值,ω为第一转动角速度矢量估计值,v2,j为第j个第二传感器的瞬时速度矢量估计值,r2,j为第j个第二传感器的坐标系相对于载体坐标系的位置平移矢量。
上述关系式基于刚体的平动速度矢量、瞬时速度矢量、转动角速度以及位置平移矢量的关系得到,且上述关系式可以存在多种变形。根据上述关系式可以较准确确定出第一平动速度矢量估计值。
在一种可能的实现方式中,所述第一平动速度矢量估计值满足如下关系式:
Figure BDA0002638325200000021
其中,
Figure BDA0002638325200000022
为第一平动速度矢量估计值,ω为第一转动角速度矢量估计值,w2,j为第j个第二传感器的加权系数,v2,j为第j个第二传感器的瞬时速度矢量估计值,r2,j为第j个第二传感器的坐标系相对于载体坐标系的位置平移矢量。
在一种可能的实现方式中,所述第一转动角速度矢量估计值满足如下关系式:
Figure BDA0002638325200000023
其中,ω为第一转动角速度矢量估计值,w1,i为第i个第一传感器的加权系数,ω1,i为第i个第一传感器的转动角速度矢量估计值。
在一种可能的实现方式中,还包括:
获取所述M个第一传感器中M′个第一传感器的归一化的平动速度矢量估计值,其中,1≤M′≤M;
根据所述第一平动速度矢量估计值、所述M′个第一传感器的归一化的平动速度矢量估计值,确定所述载体的第二平动速度矢量估计值。
上述技术方案中,第一传感器评估自身运动得到的归一化的平动速度矢量估计值较准确,该将第一传感器的归一化的平动速度矢量估计值和第一平动速度矢量估计值进行融合,得到载体的第二平动速度矢量估计值,可进一步提高载体的平动速度矢量估计值的精确度。
在一种可能的实现方式中,所述第二平动速度矢量估计值基于如下关系式确定:
Figure BDA0002638325200000024
其中,
Figure BDA0002638325200000025
为第二平动速度矢量估计值,ω为第一转动角速度矢量估计值,
Figure BDA0002638325200000026
为第i个第一传感器的归一化的平动速度矢量估计值,r1,i为第i个第一传感器的坐标系相对于载体坐标系的位置平移矢量,si为第i个第一传感器的平动速度矢量的归一化参数或者尺度伸缩因子,si由第一平动速度矢量估计值确定。
上述关系式基于刚体的平动速度矢量、瞬时速度矢量、转动角速度以及位置平移矢量的关系得到,且上述关系式可以存在多种变形。根据上述关系式可以较准确确定出第二平动速度矢量估计值。
在一种可能的实现方式中,所述第二平动速度矢量估计值满足如下关系式:
Figure BDA0002638325200000027
其中,tK为第二平动速度矢量估计值
Figure BDA0002638325200000028
tk为第k次迭代中载体的平动速度矢量估计值,ω为第一转动角速度矢量估计值,
Figure BDA0002638325200000029
为第i个第一传感器的归一化的平动速度矢量估计值,r1,i为第i个第一传感器的坐标系相对于载体坐标系的位置平移矢量,si,k为第k次迭代中第i个第一传感器的平动速度矢量的归一化参数或者尺度伸缩因子,si,k由第一平动速度矢量估计值或第k-1次迭代中载体的平动速度矢量估计值确定。
在一种可能的实现方式中,所述第二平动速度矢量估计值满足如下关系式:
Figure BDA0002638325200000031
其中,tk为第二平动速度矢量估计值
Figure BDA0002638325200000032
tk为第k次迭代中载体的平动速度矢量估计值,w′1,i,k为第k次迭代中第i个第一传感器的加权系数,ω为第一转动角速度矢量估计值,
Figure BDA0002638325200000033
为第i个第一传感器的归一化的平动速度矢量估计值,r1,i为第i个第一传感器的坐标系相对于载体坐标系的位置平移矢量,si,k为第k次迭代中第i个第一传感器的平动速度矢量的归一化参数或者尺度伸缩因子,si,k由第一平动速度矢量估计值或第k-1次迭代中载体的平动速度矢量估计值确定。
在一种可能的实现方式中,第k次迭代中第i个第一传感器的平动速度矢量的归一化参数或者尺度伸缩因子满足如下关系式:
si,k=‖tk-1+ω×r1,i
其中,tk-1为第k-1次迭代中载体的平动速度矢量估计值,t0为第一平动速度矢量估计值。
上述技术方案中,每次迭代中,根据多个第一传感器的归一化的平动速度矢量估计值,确定载体的平动速度矢量估计值。将多个第一传感器的参数作为一次迭代的输入,有助于提高第二平动速度矢量的估计精度。
在一种可能的实现方式中,所述第二平动速度矢量估计值满足如下关系式:
Figure BDA0002638325200000034
其中,
Figure BDA0002638325200000035
为第二平动速度矢量估计值
Figure BDA0002638325200000036
Figure BDA0002638325200000037
为第l轮迭代中第i个第一传感器对应的载体的平动速度矢量估计值,ω为第一转动角速度矢量估计值,
Figure BDA0002638325200000038
为第i个第一传感器的归一化的平动速度矢量估计值,r1,i第i个第一传感器的坐标系相对于载体坐标系的位置平移矢量,si,l为第l轮迭代中第i个第一传感器的平动速度矢量的归一化参数或者尺度伸缩因子,si,l由第一平动速度矢量估计值或
Figure BDA0002638325200000039
Figure BDA00026383252000000310
确定。
在一种可能的实现方式中,所述第l轮迭代中第i个第一传感器的平动速度矢量的归一化参数或者尺度伸缩因子满足如下关系式:
Figure BDA00026383252000000311
其中,
Figure BDA00026383252000000312
为第一平动速度矢量估计值;
Figure BDA00026383252000000313
上述技术方案中,针对每个第一传感器的归一化的平动速度矢量估计值,确定每个第一传感器对应的载体的平动速度矢量估计值,相当于,针对每个第一传感器进行一次迭代,将一个第一传感器的参数作为一次迭代的输入,可以快速得到较为准确的第二平动速度矢量估计值。
在一种可能的实现方式中,还包括:
基于如下关系式确定所述载体的第二转动角速度矢量估计值:
Figure BDA00026383252000000314
其中,ω′为第二转动角速度矢量估计值,v2,j为第j个第二传感器的瞬时速度矢量估计值,
Figure BDA0002638325200000041
为第二平动速度矢量估计值,
Figure BDA0002638325200000042
为[r2,j]×的逆矩阵,[r2,j]×为与r2,j对应的反对称矩阵,r2,j为第j个第二传感器的坐标系相对于载体坐标系的位置平移矢量。
上述技术方案中,根据较准确的第二平动速度矢量估计值进一步确定较准确的第二转动角速度矢量估计值。
第二方面,本申请提供一种运动估计装置,该装置包括:
获取单元和处理单元;
所述获取单元用于获取M个第一传感器的转动角速度矢量估计值和N个第二传感器的瞬时速度矢量估计值;其中,M≥1,N≥1;
所述处理单元用于根据所述N个第二传感器的瞬时速度矢量估计值和所述N个第二传感器所在的载体的第一转动角速度矢量估计值,确定所述载体的第一平动速度矢量估计值,其中,所述第一转动角速度矢量估计值是根据所述M个第一传感器的转动角速度矢量估计值确定的。
在一种可能的实现方式中,所述处理单元具体用于基于如下关系式确定所述第一平动速度矢量估计值:
Figure BDA0002638325200000043
其中,
Figure BDA0002638325200000044
为第一平动速度矢量估计值,ω为第一转动角速度矢量估计值,v2,j为第j个第二传感器的瞬时速度矢量估计值,r2,j为第j个第二传感器的坐标系相对于载体坐标系的位置平移矢量。
在一种可能的实现方式中,所述第一平动速度矢量估计值满足如下关系式:
Figure BDA0002638325200000045
其中,
Figure BDA0002638325200000046
为第一平动速度矢量估计值,ω为第一转动角速度矢量估计值,w2,j为第j个第二传感器的加权系数,v2,j为第j个第二传感器的瞬时速度矢量估计值,r2,j为第j个第二传感器的坐标系相对于载体坐标系的位置平移矢量。
在一种可能的实现方式中,所述第一转动角速度矢量估计值满足如下关系式:
Figure BDA0002638325200000047
其中,ω为第一转动角速度矢量估计值,w1,i为第i个第一传感器的加权系数,ω1,i为第i个第一传感器的转动角速度矢量估计值。
在一种可能的实现方式中,所述获取单元还用于获取所述M个第一传感器中M′个第一传感器的归一化的平动速度矢量估计值,其中,1≤M′≤M;所述处理单元还用于根据所述第一平动速度矢量估计值、所述M′个第一传感器的归一化的平动速度矢量估计值,确定所述载体的第二平动速度矢量估计值。
在一种可能的实现方式中,所述处理单元具体用于基于如下关系式确定所述第二平动速度矢量估计值:
Figure BDA0002638325200000048
其中,
Figure BDA0002638325200000051
为第二平动速度矢量估计值,ω为第一转动角速度矢量估计值,
Figure BDA0002638325200000052
为第i个第一传感器的归一化的平动速度矢量估计值,r1,i为第i个第一传感器的坐标系相对于载体坐标系的位置平移矢量,si为第i个第一传感器的平动速度矢量的归一化参数或者尺度伸缩因子,si由第一平动速度矢量估计值确定。
在一种可能的实现方式中,所述第二平动速度矢量估计值满足如下关系式:
Figure BDA0002638325200000053
其中,tK为第二平动速度矢量估计值
Figure BDA0002638325200000054
tk为第k次迭代中载体的平动速度矢量估计值,ω为第一转动角速度矢量估计值,
Figure BDA0002638325200000055
为第i个第一传感器的归一化的平动速度矢量估计值,r1,i为第i个第一传感器的坐标系相对于载体坐标系的位置平移矢量,si,k为第k次迭代中第i个第一传感器的平动速度矢量的归一化参数或者尺度伸缩因子,si,k由第一平动速度矢量估计值或第k-1次迭代中载体的平动速度矢量估计值确定。
在一种可能的实现方式中,所述第二平动速度矢量估计值满足如下关系式:
Figure BDA0002638325200000056
其中,tK为第二平动速度矢量估计值
Figure BDA0002638325200000057
tk为第k次迭代中载体的平动速度矢量估计值,w′1,i,k为第k次迭代中第i个第一传感器的加权系数,ω为第一转动角速度矢量估计值,
Figure BDA0002638325200000058
为第i个第一传感器的归一化的平动速度矢量估计值,r1,i为第i个第一传感器的坐标系相对于载体坐标系的位置平移矢量,si,k为第k次迭代中第i个第一传感器的平动速度矢量的归一化参数或者尺度伸缩因子,si,k由第一平动速度矢量估计值或第k-1次迭代中载体的平动速度矢量估计值确定。
在一种可能的实现方式中,第k次迭代中第i个第一传感器的平动速度矢量的归一化参数或者尺度伸缩因子满足如下关系式:
si,k=‖tk-1+ω×r1,i
其中,tk-1为第k-1次迭代中载体的平动速度矢量估计值,t0为第一平动速度矢量估计值。
在一种可能的实现方式中,所述第二平动速度矢量估计值满足如下关系式:
Figure BDA0002638325200000059
其中,
Figure BDA00026383252000000510
为第二平动速度矢量估计值
Figure BDA00026383252000000511
Figure BDA00026383252000000512
为第l轮迭代中第i个第一传感器对应的载体的平动速度矢量估计值,ω为第一转动角速度矢量估计值,
Figure BDA00026383252000000513
为第i个第一传感器的归一化的平动速度矢量估计值,r1,i第i个第一传感器的坐标系相对于载体坐标系的位置平移矢量,si,l为第l轮迭代中第i个第一传感器的平动速度矢量的归一化参数或者尺度伸缩因子,si,l由第一平动速度矢量估计值或
Figure BDA00026383252000000514
Figure BDA00026383252000000515
确定。
在一种可能的实现方式中,所述第l轮迭代中第i个第一传感器的平动速度矢量的归一化参数或者尺度伸缩因子满足如下关系式:
Figure BDA00026383252000000516
其中,
Figure BDA00026383252000000517
为第一平动速度矢量估计值;
Figure BDA00026383252000000518
在一种可能的实现方式中,所述处理单元还用于基于如下关系式确定所述载体的第二转动角速度矢量估计值:
Figure BDA00026383252000000519
其中,ω′为第二转动角速度矢量估计值,v2,j为第j个第二传感器的瞬时速度矢量估计值,
Figure BDA0002638325200000061
为第二平动速度矢量估计值,
Figure BDA0002638325200000062
为[r2,j]×的逆矩阵,[r2,j]×为与r2,j对应的反对称矩阵,r2,j为第j个第二传感器的坐标系相对于载体坐标系的位置平移矢量。
第三方面,本申请提供一种通信装置,包括至少一个处理器和通信接口,所述通信接口用于接收来自所述通信装置之外的其它通信装置的信号并传输至所述至少一个处理器或将来自所述至少一个处理器的信号发送给所述通信装置之外的其它通信装置,所述至少一个处理器通过逻辑电路或执行代码指令用于实现上述第一方面或第一方面的任意可能的实现方式中的方法。
第四方面,本申请提供一种计算机可读存储介质,所述计算机可读存储介质中存储有计算机程序或指令,当所述计算机程序或指令被通信装置执行时,实现上述第一方面或第一方面的任意可能的实现方式中的方法。
第五方面,本申请提供一种计算机程序产品,所述计算机程序产品包括计算机程序或指令,当所述计算机程序或指令被通信装置执行时,实现上述第一方面或第一方面的任意可能的实现方式中的方法。
第六方面,本申请提供一种芯片,包括至少一个处理器和接口;
所述接口,用于为所述至少一个处理器提供程序指令或者数据;
所述至少一个处理器用于执行所述程序行指令,以实现上述第一方面或第一方面的任意可能的实现方式中的方法。
第七方面,本申请提供一种终端,所述终端包含上述第二方面提供的任意运动估计装置、上述第三方面提供的任意通信装置、或者上述第四方面提供的任意计算机可读存储介质。进一步可选的,所述终端可以为车辆、无人机、机器人、智能家居设备或者卫星等。
上述第二方面至第七方面中任一方面可以达到的技术效果可以参照上述第一方面中有益效果的描述,此处不再重复赘述。
附图说明
图1为本申请提供的一种自运动估计系统的结构示意图;
图2为本申请提供的一种运动估计方法的流程示意图;
图3为本申请提供的一种车载系统中多传感器配置示意图;
图4为本申请提供的一种转动角速度矢量的示意图;
图5为本申请提供的一种位置平移矢量的示意图;
图6为本申请提供的一种载体坐标系与传感器坐标系之间的变换关系的示意图;
图7为本申请提供的一种运动估计装置的结构示意图;
图8为本申请提供的一种芯片的结构示意图。
具体实施方式
下面将结合附图,对本申请实施例进行详细描述。
请参见图1,图1是本申请实施例提供的一种自运动估计系统的结构示意图,该系统包括第一传感器1010、第一运动感知模块1011、第二传感器1020、第二运动感知模块1021和数据处理模块1030。
第一传感器1010可以为视觉传感器,例如,相机或者摄像头、红外热成像传感器等。第一传感器1010可以提供视觉测量数据,例如图像或者视频。第一运动感知模块1011用于根据第一传感器1010提供的测量数据确定运动测量数据,如传感器运动的转动角速度矢量和/或归一化或尺度伸缩(scaled)的平动速度矢量等或者带有完整尺度信息的平动速度矢量估计值。第二传感器1020可以为雷达传感器、超声波传感器、惯性测量传感器或者定位传感器等,例如,毫米波雷达、声呐、激光雷达、惯性测量单元(inertial measurementunit,IMU)或者全球导航卫星系统(global navigation satellite system,GNSS)等。第二传感器1020用于提供位置测量数据和/或速度测量数据,例如位置和/或径向速度或者速度投影分量的测量数据。第二运动感知模块1021用于根据第二传感器1020提供的位置测量数据和/或速度测量数据,确定运动测量数据,如传感器运动的瞬时平动速度矢量等。数据处理模块1030用于处理第一运动感知模块1011和第二运动感知模块1021提供的运动测量数据。本申请中,运动测量数据也可以称为是运动感知数据。
第一传感器1010、第一运动感知模块1011、第二传感器1020、第二运动感知模块1021和数据处理模块1030可以通过有线或者无线方式等连接在一起;第一传感器1010和第二传感器1020可以分布于载体的相同或者不同的位置;第一运动感知模块1011和第二运动感知模块1021可以分别与第一传感器1010和第二传感器1020集成在一起;也可以与数据处理模块1030集成在一起;也可以分别独立于其它模块存在,本申请不做限定。
在一种示例中,第一传感器1010、第一运动感知模块1011、第二传感器1020、第二运动感知模块1021和数据处理模块1030部署于一个处理器系统上。在又一种示例中,第一传感器1010、第二传感器1020分别部署于一个处理器系统上;第一运动感知模块1011、第二运动感知模块1021和数据处理模块1030部署于一个处理器系统上。
本申请实施例可适用各种载体中的多传感器系统,其中载体如车载(如汽车、摩托车或自行车等)、机载(如无人机、直升机或喷气式飞机或者气球)、船载(轮船、汽艇或者舰艇等)、星载(如卫星)或者智能体(如机器人等)等。
示例性的,载体为车载,车载上可以承载有至少一个第一传感器和至少一个第二传感器。比如,车载安装有1个第一传感器和1个第二传感器,或者,车载安装有1个第一传感器和5个第二传感器,或者,车载安装有6个第一传感器和5个第二传感器等。
目前,获得传感器运动速度或者载体运动速度估计的方法有多种。以车载平台为例,车载的运动估计中,可以是采用如下三种传感器。
一、IMU
IMU是测量物体三轴姿态角(或角速度)以及加速度的装置。一般地,一个IMU内会装有三轴的陀螺仪和三个方向的加速度计,来测量物体在三维空间中的角速度和加速度,并以此可以解算出物体运动速度和姿态。
二、雷达传感器
雷达传感器,通常可以提供距离、方位角和径向速度测量数据。基于静止目标物体的方位角和径向速度分量测量数据,利用最小二乘方法或者其他的方式,可以获得传感器相对于大地的瞬时速度,特别的,雷达传感器可以获取到较为精确的纵向速度估计值。此外,基于上述速度估计值,可以进一步得到传感器的横摆角速度(yaw rate)估计值。
三、视觉传感器
视觉传感器通常可以提供连续的两帧或者多帧图像。基于上述两帧或者多帧图像,利用光流法或者特征点对应的方法或者直接优化光亮度(intensity)的目标物体函数的方法,可以获得传感器尺度伸缩的平动速度估计值和转动速度估计值。
但上述三种传感器均存在各自缺陷:
(1)IMU的运动速度估计是基于加速度计积累得到,测量误差会随时间累积,因而存在误差累积问题,需要利用其它传感器进行额外的校准,而且,一般用于车辆中的IMU的精度太低,若选用高精度的IMU,成本高昂。
(2)雷达传感器获取到的传感器的横向速度估计值的精确度较低,且不能获取传感器的俯仰角速度估计值(pitch rate)和滚转角速度(roll rate)估计值。
(3)视觉传感器存在尺度伸缩问题,深度信息与平动速度的各个分量耦合在一起,通常无法获得精确的深度估计,无法精确获得平动速度估计值,或者仅能获得平动速度的一个尺度伸缩的估计值。
为解决上述问题,本申请实施例提供一种运动估计的方法,用于较精确的确定载体的平动速度矢量估计值和转动角速度矢量估计值,从而对载体实现精确的运动估计。
需要指出的是,本申请中,平动速度矢量估计值也可以是平动位移矢量估计值,平动位移矢量估计值可以是两帧间的位置偏移矢量的估计值,也可以是两帧间的时间差和平动速度矢量估计值的乘积。
请参见图2,图2是本申请实施例提供的一种运动估计方法的流程示意图,该方法的执行主体可以是传感器系统或者融合感知系统或者集成上述系统的规划/控制系统如辅助驾驶或者自动驾驶系统等。或者,该方法的执行主体也可以是软件或者硬件(如与相应传感器通过无线或者有线连接或者集成在一起的数据处理装置)。
以下不同的执行步骤可以集中式实现,或者,以下不同的执行步骤也可以分布式实现。方法包括但不限于如下步骤:
步骤201,获取M个第一传感器的转动角速度矢量估计值和N个第二传感器的瞬时速度矢量估计值。
其中,M个第一传感器和N个第二传感器承载于一个载体上,第一传感器的数量M≥1,第二传感器的数量N≥1。
具体的,第一传感器可以为视觉传感器如相机、摄像头、红外传感器或其它成像传感器等,或惯性测量传感器如IMU,第二传感器可以为雷达传感器如毫米波雷达、激光雷达等,或超声波传感器如声呐等。
一种实现中,M(M>1)个第一传感器可以是相同类型或不同类型,示例性的,M个第一传感器中包括M1个视觉传感器如相机或者摄像头和M2个惯性测量传感器如IMU,其中,M1+M2=M,且M1≥0,M2≥0。和/或,N(N>1)个第二传感器可以是相同类型或不同类型,示例性的,N个第二传感器中包括N1个雷达传感器如毫米波雷达、激光雷达和N2个超声波传感器如声呐,其中,N1+N2=N,且N1≥0,N 2≥0。
M个第一传感器和N个第二传感器可以安装于载体的相同位置或者不同位置。下面以车载为例说明:
例子1,车载前端安装1个相机和1个毫米波雷达;
例子2,车载前端安装1个相机和1个毫米波雷达,车载4个角位置另外安装4个毫米波雷达;
例子3,车载前端安装1个毫米波雷达,4个角位置另外安装4个毫米波雷达,另外在车载上均匀安装6个相机。
此外,车载上还可以进一步安装IMU或者GNSS等。
示例性的,一种车载系统中多传感器配置示意图如图3所示,车载系统的多传感器可以包括1个相机、5个毫米波雷达和1个IMU,其中IMU的安装位置可以接近车载坐标系(也可以称为是车体坐标系、载体坐标系)的原点位置,车载坐标系的原点可以位于车体后轴中心。
上述例子1至例子3中,毫米波雷达也可以替换成激光雷达或超声波传感器等,相机也可以替换成摄像头或红外传感器等,也可以在原基础上增加至少一个激光雷达。
比如上述例子3中,可以将5个毫米波雷达的其中2个毫米波雷达替换为激光雷达,又或者将5个毫米波雷达全部替换为激光雷达,此外,还可以在原5个毫米波雷达的基础上增加1至3个激光雷达。
应理解,此处只是本申请示例性提供的车载上安装第一传感器和第二传感器的实现方式,其并不对本申请做任何限定。
具体的,所述获取M个第一传感器的转动角速度矢量估计值,可以是从传感器通过有线或者无线接口直接得到,其中转动角速度矢量估计值可以基于传感器的测量数据通过运动或者估计算法得到或者是由传感器直接测量得到;
或者,所述获取M个第一传感器的转动角速度矢量估计值,可以是从传感器通过有线或者无线接口直接得到传感器的测量数据,转动角速度矢量估计值根据传感器的测量数据通过运动或者估计算法得到或者是由传感器测量数据直接得到。
作为一种实现方式,转动角速度矢量估计值可以根据第一传感器的测量数据通过估计得到。此时,第一传感器的测量数据并不直接包含第一传感器的运动测量的测量值。
示例性的,第一传感器为视觉传感器如相机,相机的原始测量数据为视觉测量数据,比如图像或者视频,可以根据图像或者视频中的特征点、线、平面或者区域,基于数据的光学特性或者几何特性确定相机的转动角速度矢量估计值。例如,可以基于8点法或者5点法或者单应性(Homography)或者光流法等,得到相机的转动角速度矢量估计值。基于图像或者视频,得到传感器的转动角速度为现有技术,此处不赘述。
作为另一种实现方式,转动角速度矢量估计值可以从第一传感器的测量数据中直接得到。此时,第一传感器可以直接测量包含转动角速度矢量。
示例性的,第一传感器为惯性测量传感器如IMU,IMU可以直接测量转动角速度矢量。
需要指出的是,上述M个第一传感器可以包含相同类型或不同类型的第一传感器,示例性的,M个第一传感器可以包括M1个视觉传感器例如相机或者摄像头和M2个惯性测量传感器如IMU,其中,M1+M2=M,且M1≥0,M2≥0,相应的,对M1个视觉传感器的测量数据处理得到M1个视觉传感器的转动角速度矢量估计值,以及从M2个惯性测量传感器的测量数据中直接得到读取M2个惯性测量传感器的转动角速度矢量估计值。
可选地,作为一种实现方式,转动角速度矢量可以为三维矢量,ω=[ωx ωy ωz]T,如图4所示。
可选地,作为另一种实现方式,第一传感器或者载体在平面内运动,平面如地面或者平面轨道。转动角速度矢量可以表示为ω=[0 0 ωz]T,此时,转动角速度矢量可以简化为用ωz表示。
需要指出的是,上述从传感器的测量数据直接得到或者基于传感器的测量数据通过运动估计得到转动角速度矢量估计值,可以根据传感器坐标系相对于载体坐标的变换关系,得到传感器在载体坐标系的转动角速度矢量估计值。
具体的,所述获取N个第二传感器的瞬时速度矢量估计值,可以是从传感器通过有线或者无线接口直接得到,其中瞬时速度矢量估计值可以基于传感器的测量数据通过运动或者估计算法得到;
或者,所述获取N个第二传感器的瞬时速度矢量估计值,可以是从传感器通过有线或者无线接口直接得到传感器的测量数据,瞬时速度矢量估计值根据传感器的测量数据通过运动或者估计算法得到。
作为一种实现方式,瞬时速度矢量估计值可以根据第二传感器的测量数据通过估计得到。此时,第二传感器的测量数据并不直接包含第二传感器的运动测量的测量值。
示例性的,第二传感器为毫米波雷达或者激光雷达或者超声波传感器如声呐,第二传感器的测量数据可以包括位置和径向速度,或者包括角度和径向速度。可以根据其中静止目标的测量数据,基于最小二乘法或者正交距离回归法或者最小均方误差准则等估计方法,确定瞬时速度矢量估计值。此外,也可以根据第二传感器的多次位置测量数据,根据其中静止目标的测量数据,确定瞬时速度矢量估计值。本申请实施例不做限定。
需要指出的是,载体获取第一传感器的转动角速度矢量估计值的数量可以小于或等于实际上载体中第一传感器的数量,载体获取第二传感器的瞬时速度矢量估计值的数量可以小于或等于实际上载体中第二传感器的数量。
第一个示例性中,载体上承载有6个第一传感器和3个第二传感器,载体可以获取6个第一传感器的转动角速度矢量估计值和3个第二传感器的瞬时速度矢量估计值。
第二个示例性中,载体上承载有6个第一传感器和3个第二传感器,载体可以获取4个第一传感器的转动角速度矢量估计值和2个第二传感器的瞬时速度矢量估计值。
步骤202,根据N个第二传感器的瞬时速度矢量估计值、N个第二传感器的外部参数和N个第二传感器所在的载体的第一转动角速度矢量估计值,确定载体的第一平动速度矢量估计值。
其中,N个第二传感器的外部参数可以包括N个第二传感器相对于载体坐标系的位置平移矢量,或者,包括N个第二传感器的坐标系原点相对于载体坐标系原点的位置平移矢量。示例性的,对于任一个第二传感器来说,第二传感器相对于载体坐标系的位置平移矢量用于将该第二传感器的坐标系原点平移至与载体坐标系原点一致。
以车载上安装5个第二传感器为例,该5个第二传感器可以包括毫米波雷达或者激光雷达或者超声波传感器,参照图5所示,5个第二传感器分别位于车载的不同位置,该5个第二传感器相对于车载坐标系原点的位置平移矢量分别为r21,r22,…,r25,该5个第二传感器的瞬时速度矢量估计值分别为v21,v22,…,v25。通常r21,r22,…,r25互不相同时,v21,v22,…,v25也互不相同。
具体的,根据N个第二传感器的瞬时速度矢量估计值和N个第二传感器的外部参数以及载体的第一转动角速度矢量估计值,确定载体的第一平动速度矢量估计值,可以是,载体的第一平动速度矢量估计值基于刚体的平动速度矢量、瞬时速度矢量、转动角速度以及位置平移矢量的关系得到,其中瞬时速度矢量以及位置平移矢量从N个第二传感器的瞬时速度矢量估计值和N个第二传感器的外部参数确定,转动角速度从载体的第一转动角速度矢量估计值确定。
具体的,刚体的平动速度矢量、瞬时速度矢量、转动角速度矢量以及位置平移矢量的关系式为t=v-ω×r,或者t=v+r×ω,其中,t为刚体的平动速度矢量,v为刚体的瞬时速度矢量,ω为刚体的转动角速度矢量,r为位置平移矢量,×表示矢量的叉积。
基于刚体的平动速度矢量、瞬时速度矢量、转动角速度矢量以及位置平移矢量的关系式,可以存在多种变形关系式,该多种变形关系式均可以得到载体的第一平动速度矢量估计值。
示例性的,确定第一平动速度矢量估计值可以基于如下关系式得到:
Figure BDA0002638325200000111
或者
Figure BDA0002638325200000112
其中,
Figure BDA0002638325200000113
为载体的第一平动速度矢量估计值,ω为载体的第一转动角速度矢量估计值,v2为第二传感器的瞬时速度矢量估计值,r2为第二传感器相对于载体坐标系原点的位置平移矢量,×表示矢量的叉积。
具体的,上述关系式可以为
Figure BDA0002638325200000114
其中,ω=[ωx ωy ωz]T
具体的,上述关系式也可以为
Figure BDA0002638325200000115
其中,r2=[rx,2 ry,2 rz,2]T
当然也可以根据其他的方式,本申请实施例不做限定。
具体的,r2可以为第二传感器的外部参数,可以为第二传感器的坐标系原点相对于载体坐标系原点的位置平移矢量。
作为一种实现方式,第一平动速度矢量估计值可以是基于如下关系式确定:
Figure BDA0002638325200000116
其中,
Figure BDA0002638325200000117
为载体的第一平动速度矢量估计值,t2,j为根据第j个第二传感器的瞬时速度矢量估计值及其外部参数和载体的第一转动角速度矢量估计值确定的载体的平动速度矢量估计值。
作为一种实现方式,基于最小均方误差(minimum mean sqaured error,MMSE)或者最小二乘(least square,LS)法等确定第一平动速度矢量估计值。
第一种实现方式中,第一平动速度矢量估计值是载体的N个平动速度矢量估计值的加权和,其中,载体的N个平动速度矢量估计值分别根据N个第二传感器的瞬时速度矢量估计值及其外部参数和载体的第一转动角速度矢量估计值确定。
一种具体实现中,载体的第一平动速度矢量估计值满足如下关系式:
Figure BDA0002638325200000118
其中,
Figure BDA0002638325200000119
为载体的第一平动速度矢量估计值,t2,j为根据第j个第二传感器的瞬时速度矢量估计值及其外部参数和载体的第一转动角速度矢量估计值确定的载体的平动速度矢量估计值,w2,j为与第j个第二传感器对应的加权系数或者加权系数矩阵。
具体的,加权系数w2,j可以根据t2,j的估计误差或测量误差的概率密度函数或统计特性或协方差矩阵确定。例如,w2,j可以根据t2,j的估计误差或测量误差的协方差矩阵确定,
Figure BDA0002638325200000121
其中,
Figure BDA0002638325200000122
P2,j为t2,j的估计误差或测量误差的协方差。
示例性的,根据第j个第二传感器的瞬时速度矢量估计值及其外部参数和载体的第一转动角速度矢量估计值确定载体的平动速度矢量估计值,可以是基于如下关系式确定:
t2,j=v2,j-ω×r2,j
相应的,作为第一种实现方式的一种具体示例,可以根据如下关系式确定第一平动速度矢量估计值:
Figure BDA0002638325200000123
其中,
Figure BDA0002638325200000124
为载体的第一平动速度矢量估计值,ω为载体的第一转动角速度矢量估计值,v2,j为第j个第二传感器的瞬时速度矢量估计值,r2,j为第j个第二传感器的坐标系相对于载体坐标系的位置平移矢量,w2,j为与第j个第二传感器对应的加权系数或者加权系数矩阵。
具体的,加权系数w2,j可以根据v2,j和/或r2,j测量误差或者估计误差确定。
例如,w2,j根据v2,j和r2,j测量误差或者估计误差确定,例如,w2,j满足如下关系式:
Figure BDA0002638325200000125
其中,
Figure BDA0002638325200000126
其中,
Figure BDA0002638325200000127
为v2,j的测量误差或者估计误差的协方差矩阵,
Figure BDA0002638325200000128
为r2,j的测量误差或者估计误差的协方差矩阵。Ω根据载体的转动角速度矢量得到,具体为,
Figure BDA0002638325200000129
其中,ω=[ωx ωy ωz]T
第二种实现方式中,第一平动速度矢量估计值是载体的N个平动速度矢量估计值的均值,其中,载体的N个平动速度矢量估计值分别根据N个第二传感器的瞬时速度矢量估计值及其外部参数和载体的第一转动角速度矢量估计值确定。
一种具体实现中,载体的第一平动速度矢量估计值满足如下关系式:
Figure BDA00026383252000001210
其中,
Figure BDA00026383252000001211
为载体的第一平动速度矢量估计值,t2,j为根据第j个第二传感器的瞬时速度矢量估计值及其外部参数和载体的第一转动角速度矢量估计值确定的载体的平动速度矢量估计值。
示例性的,根据第j个第二传感器的瞬时速度矢量估计值及其外部参数和载体的第一转动角速度矢量估计值确定载体的平动速度矢量估计值,可以是基于t2,j=v2,j-ω×r2,j。相应的,作为第二种实现方式的一种具体示例,第一平动速度矢量估计值满足如下关系式:
Figure BDA00026383252000001212
其中,
Figure BDA00026383252000001213
为载体的第一平动速度矢量估计值,ω为载体的第一转动角速度矢量估计值,v2,j为第j个第二传感器的瞬时速度矢量估计值,r2,j为第j个第二传感器的坐标系相对于载体坐标系的位置平移矢量。
此处,需要说明的是,第一种实现方式中,第一平动速度矢量估计值是载体的N个平动速度矢量估计值的加权和,而第二种实现方式中,第一平动速度矢量估计值是载体的N个平动速度矢量估计值的均值,也可以理解,第二种实现方式是第一种实现方式的一种特殊形式,若第一种实现方式中各第二传感器的加权系数相同,则可以为第二种实现方式。进一步的,该说明同样适用于其它实现方式中多个传感器估计值的加权和与多个传感器估计值的均值之间的关系。
作为另一种实现方式,N等于1,根据第一转动角速度矢量估计值、第二传感器的瞬时速度矢量估计值、第二传感器的坐标系相对于载体坐标系的位置平移矢量,确定第一平动速度矢量估计值,具体可参照上述关系式得到。
此外,还可以根据载体的第一转动角速度矢量估计值、第j个第二传感器的瞬时速度矢量估计值和第j个第二传感器的坐标系相对于载体坐标系的位置平移矢量,通过正交距离回归(orthogonal distance regression,ODR)方法,确定载体的平动速度矢量估计值t2,j
本申请实施例中,第一转动角速度矢量估计值是根据M个第一传感器的转动角速度矢量估计值确定的,其中M≥1。
作为一种实现方式,第一转动角速度矢量估计值,可以基于如下关系式确定:
ω=ω1,i,i=1,…,M,M≥1
其中,ω为第一转动角速度矢量估计值,ω1,i为第i个第一传感器的转动角速度矢量估计值。
作为一种实现方式,基于最小均方误差或者最小二乘法等确定第一转动角速度矢量估计值。
第1种实现方式中,第一转动角速度矢量估计值是载体的M个转动角速度矢量估计值的加权和,其中,载体的M个转动角速度矢量估计值分别根据M个第一传感器的转动角速度矢量估计值确定。
一种具体实现中,载体的第一转动角速度矢量估计值满足如下关系式:
Figure BDA0002638325200000131
其中,ω为载体的第一转动角速度矢量估计值,ω1,i为根据第i个第一传感器的转动角速度矢量估计值得到载体的转动角速度矢量估计值,w1,i为与第i个第一传感器对应的加权系数或者加权系数矩阵。
具体的,加权系数w1,i可以根据ω1,i的估计误差或测量误差的概率密度函数或统计特性或协方差矩阵确定。例如,w1,i可以根据ω1,i的估计误差或测量误差的协方差矩阵确定,具体为
Figure BDA0002638325200000132
其中,
Figure BDA0002638325200000133
P1,i为ω1,i的估计误差或测量误差的协方差。
第2种实现方式中,第一转动角速度矢量估计值是载体的M个转动角速度矢量估计值的均值,其中,载体的M个转动角速度矢量估计值分别根据M个第一传感器的转动角速度矢量估计值确定。
一种具体实现中,载体的第一转动角速度矢量估计值满足如下关系式:
Figure BDA0002638325200000134
其中,ω为载体的第一转动角速度矢量估计值,ω1,i为根据第i个第一传感器的转动角速度矢量估计值得到的载体的转动角速度矢量估计值。
作为另一种实现方式,M等于1,根据第一传感器的转动角速度矢量估计值,确定第一转动角速度矢量估计值,具体可参照上述关系式得到。
需要指出的是,以上,第一传感器的转动角速度矢量估计值、第二传感器的瞬时速度矢量估计值、载体的第一转动角速度矢量估计值、载体的第一平动速度矢量估计值等均为相对于载体坐标系定义的。而实际应用中,传感器的测量数据往往相对于传感器坐标系定义,因此,从第一传感器或者第二传感器得到运动速度矢量,包括转动角速度矢量和平动速度矢量或者瞬时速度矢量,往往相对于所在传感器定义比较方便。此时,需要基于传感器相对于载体坐标系的外部参数得到相对于载体坐标系的转动角速度矢量和平动速度矢量或者瞬时速度矢量。
不失一般性的,如图6所示,载体坐标系与传感器坐标系之间的变换关系通常可以用传感器的外部参数确定,外部参数可以包括旋转参数和传感器坐标系相对于载体坐标系的位置平移矢量。基于旋转参数,可将传感器坐标系的方向旋转至与载体坐标系的方向一致。基于位置平移矢量,可将传感器坐标系的原点平移至与载体坐标系的原点一致。
传感器坐标系相对于载体坐标系的位置平移矢量可以是如图6中矢量r,或者可以是如图5中的r21,r22,…,r25;再如前述r2,j,j=1,…,N,N≥1。
旋转参数用于表示载体坐标系和传感器坐标系之间的旋转;具体的,可以是用四元数、旋转矩阵、欧拉角等表示。其中,四元数、旋转矩阵、欧拉角等可以互相转化。示例性的,旋转矩阵可以根据四元数得到,或旋转矩阵可以根据欧拉角得到。示例性的,传感器坐标系与载体坐标系的方向可以一致,此时,旋转参数为单位矩阵。
本申请实施例,进一步基于传感器的外部参数和相对于传感器坐标系的运动速度矢量,得到相对于载体坐标系的运动速度矢量,其中,运动速度矢量可以包括转动角速度、平动速度矢量、瞬时运动速度矢量中的一个或多个。
具体的,可以基于第二传感器的外部参数和相对于第二传感器坐标系的瞬时速度矢量,得到相对于载体坐标系的瞬时速度矢量,其中外部参数包括旋转参数。
例如,旋转参数可以是旋转矩阵,相对于载体坐标系的瞬时平动速度矢量可以由以下关系式确定:
v2,j=R2,jv′2,j
其中,v2,j为第j个第二传感器相对于载体坐标系的瞬时速度矢量估计值,v′2,j第j个第二传感器相对于传感器坐标系的瞬时速度矢量估计值,R2,j为第二传感器的传感器坐标系到载体坐标系的旋转参数。
具体的,可以基于第一传感器的外部参数和相对于第一传感器坐标系的转动角速度矢量,得到相对于载体坐标系的转动角速度矢量,其中外部参数包括旋转参数。
例如,旋转参数可以是旋转矩阵,相对于载体坐标系的转动角速度矢量可以由以下关系式确定:
ω1,i=R1,iω′1,i
其中,ω1,i为第i个第一传感器相对于载体坐标系的转动角速度矢量估计值,ω′1,i为第i个第一传感器相对于该第一传感器的坐标系的转动角速度矢量估计值,R1,i为第一传感器的传感器坐标系到载体坐标系的旋转矩阵。
需要指出的是,旋转参数可以是固定值,也可以是运动估计中根据在线算法估计得到;位置平移矢量可以是固定值,也可以是运动估计中根据在线算法估计得到。本申请实施例对此不限定。
基于上述变换关系,可以根据传感器的外部参数、M个第一传感器的转动角速度矢量估计值,确定载体的第一转动角速度矢量估计值。其中传感器的外部参数可以包括传感器的旋转参数。
具体的,根据M个第一传感器的转动角速度矢量估计值,确定载体的第一转动角速度矢量估计值。
作为一种实现方式,可以基于如下关系式,确定第一转动角速度矢量估计值。
ω=R1,iω′1,i
其中,ω为第一转动角速度矢量估计值,ω′1,i为第i个第一传感器相对于传感器坐标系的转动角速度矢量估计值,R1,i为第i个第一传感器的坐标系相对于载体坐标系的旋转参数。
作为一种实现方式,根据M个第一传感器中各第一传感器相对于传感器坐标系的转动角速度矢量估计值、各第一传感器的旋转参数,确定载体的第一转动角速度矢量估计值。
在一个示例中,基于如下关系式确定第一转动角速度矢量估计值:
Figure BDA0002638325200000151
其中,ω为载体的第一转动角速度矢量估计值,ω′1,i为第i个第一传感器相对于传感器坐标系的转动角速度矢量估计值,w1,i为与第i个第一传感器对应的加权系数或者加权系数矩阵,R1,i为第i个第一传感器的坐标系相对于载体坐标系的旋转参数。
在又一个示例中,基于如下关系式确定第一转动角速度矢量估计值:
Figure BDA0002638325200000152
其中,ω为载体的第一转动角速度矢量估计值,ω′1,i为第i个第一传感器相对于传感器坐标系的转动角速度矢量估计值,R1,i为第i个第一传感器的坐标系相对于载体坐标系的旋转参数。
本申请实施例中,根据至少一个第一传感器的转动角速度矢量估计值,可以准确确定载体的第一转动角速度矢量估计值。结合第二传感器的瞬时速度矢量估计值,可以准确得到载体的第一平动速度矢量估计值。因此,采用本申请实施例的方法,可以实现载体的运动进行补偿,有助于实现运动目标和静止目标的分离,同时有助于实现对载体运动的定位和跟踪。
可选地,为了进一步提升载体的平动速度矢量估计值的准确性,本申请实施例可以进一步包括如下步骤203和步骤204。
步骤203,获取M个第一传感器中M′个第一传感器的归一化的平动速度矢量估计值。
具体的,M个第一传感器中包括M′个第一传感器,其中M′≤M,该M′个第一传感器可以是视觉传感器例如相机或者摄像头。
一个具体示例中,M个第一传感器包括M′个视觉传感器例如相机或者摄像头和M-M′个惯性测量传感器例如IMU。
作为一种实现方式,M′个第一传感器可以是视觉传感器例如相机或者摄像头,获取M′个第一传感器的归一化的平动速度矢量估计值,可以是根据相机或者摄像头得到的图像或者视频,根据其中的特征点、线、或者平面或者区域,基于数据的光学特性或者几何特性确定,例如基于8点法或者5点法或者单应性(Homography)或者光流法等方法得到,本申请实施例不做限定。
第一传感器的归一化的平动速度矢量估计值,理解为,根据第一传感器确定的尺度伸缩的平动速度矢量估计值,也即,第一传感器的归一化的平动速度矢量估计值与第一传感器的平动速度矢量估计值成比例。应理解,第一传感器的平动速度矢量估计值为第一传感器实际运动的平动速度矢量估计值。
示例性的,在第一传感器的坐标系中,第一传感器的归一化的平动速度矢量估计值可以表示为
Figure BDA0002638325200000161
第一传感器的平动速度矢量估计值可以表示为v1′,二者之间符合关系式:
Figure BDA0002638325200000162
其中,v1′=[v′1x,v′1y,v′1z]T
Figure BDA0002638325200000163
s为第一传感器的平动速度矢量估计值的归一化参数或者尺度伸缩因子。
具体的,归一化参数或者尺度伸缩因子可以是平动速度矢量估计值的幅度或者范数或者模,或平动速度矢量估计值的某一分量,如z轴分量。
以光流法为例,可以基于以下关系式确定归一化的平动速度矢量估计值
Figure BDA0002638325200000164
Figure BDA0002638325200000165
Figure BDA0002638325200000166
其中u,v为图像平面上的光流分量,s1=[-f 0 x]T,s2=[0 -f y]T
Figure BDA0002638325200000167
×表示矢量的叉积,f为相机焦距,x,y为图像平面的像素位置,x∈[px-wx,px+wx],y∈[py-wy,py+wy];其中(px,py)为中心位置,wx和wy为非负整数,wx=0,1,2,3,4…;wy=0,1,2,3,4…。Z′为像素对应的目标点的相对深度,
Figure BDA0002638325200000168
和ω为相对于传感器坐标系的归一化平动速度矢量和转动角速度矢量。
其中归一化的平动速度矢量估计值
Figure BDA0002638325200000169
和相对深度Z′,满足如下关系
Figure BDA00026383252000001610
或者,等价地
Figure BDA00026383252000001611
Z=sZ′
其中t′为相对于传感器坐标系的绝对平动速度矢量,Z为像素对应的目标点的绝对深度,t′z为t′的z轴分量,s=t′z为尺度伸缩因子。需要指出的是,尺度伸缩因子s不限于t′z,可以根据需要选择其它值,例如尺度伸缩因子s为t′的范数或者幅度等,如s=‖t′‖。
需要指出的是,获取M′个第一传感器的归一化的平动速度矢量估计值,可以进一步包括,根据第一传感器的外部参数得到相对于载体坐标系的归一化的平动速度矢量估计值,其中第一传感器的外部参数包括第一传感器坐标系相对于载体坐标系的旋转参数。
例如,相对于载体坐标系的归一化的平动速度矢量估计值
Figure BDA00026383252000001612
Figure BDA00026383252000001613
其中,R1为第一传感器的坐标系相对于载体坐标系的旋转矩阵。
需要指出的是,上述旋转变换为正交变换矩阵,并不改变归一化参数。下面,为方便描述,可以以归一化参数为例说明。
步骤204,根据第一平动速度矢量估计值、M′个第一传感器的归一化的平动速度矢量估计值,确定载体的第二平动速度矢量估计值。
本申请中,在确定出第一平动速度矢量估计值的基础上,可以基于M′个第一传感器的归一化的平动速度矢量估计值,进一步提高载体的平动速度矢量的估计精度,此处将进一步更新得到的载体的平动速度矢量估计值称为第二平动速度矢量估计值。
具体的,可以根据如下归一化的第一传感器的平动速度矢量、第一传感器的位置平移矢量、载体的平动速度矢量和载体的转动角速度矢量的关系,确定载体的第二平动速度矢量估计值:
Figure BDA0002638325200000171
其中,t为载体的平动速度矢量,ω为载体的转动角速度矢量估计值,
Figure BDA0002638325200000172
为第一传感器的归一化的平动速度矢量估计值,r1为第一传感器相对于载体坐标系的位置平移矢量,s为第一传感器的平动速度矢量的归一化参数。
作为一种实现方式,可以根据载体的第一平动速度矢量、载体的第一转动角速度矢量、第一传感器的归一化的平动速度矢量以及第一传感器相对于载体坐标系的位置平移矢量的关系,确定载体的第二平动速度矢量估计值:
具体的,载体的第二平动速度矢量估计值根据以下关系式确定:
Figure BDA0002638325200000173
其中,
Figure BDA0002638325200000174
为载体的第二平动速度矢量估计值,ω为载体的第一转动角速度矢量估计值,
Figure BDA0002638325200000175
为第i个第一传感器的归一化的平动速度矢量估计值,r1,i为第i个第一传感器相对于载体坐标系的位置平移矢量,si为第i个第一传感器的平动速度矢量的归一化参数。
具体的,载体的第二平动速度矢量估计值根据以下关系式确定:
Figure BDA0002638325200000176
其中,
Figure BDA0002638325200000177
为载体的第二平动速度矢量估计值,w′1,i为与第i个第一传感器对应的加权系数矩阵,ω为载体的第一转动角速度矢量估计值,
Figure BDA0002638325200000178
为第i个第一传感器的归一化的平动速度矢量估计值,r1,i为第i个第一传感器相对于载体坐标系的位置平移矢量,si为第i个第一传感器的平动速度矢量的归一化参数。加权系数矩阵w′1,i可以根据
Figure BDA0002638325200000179
的估计误差的协方差矩阵确定,与之前方法类似,此处不详述。
具体的,载体的第二平动速度矢量估计值也根据以下关系式确定:
Figure BDA00026383252000001710
其中ω为载体的第一转动角速度矢量估计值,
Figure BDA00026383252000001711
为第i个第一传感器的归一化的平动速度矢量估计值,r1,i为第i个第一传感器相对于载体坐标系的位置平移矢量,si为第i个第一传感器的平动速度矢量的归一化参数。
可选的,si具体可以根据下述关系式确定:
Figure BDA00026383252000001712
其中
Figure BDA00026383252000001713
为载体的第一平动速度矢量估计值。
作为另一种实现方式,可以根据上述载体的平动速度矢量、载体的第一转动角速度矢量和第一传感器的平动速度矢量以及第一传感器相对于载体坐标系的位置平移矢量的关系,确定载体的第二平动速度矢量估计值,符合以下关系式:
Figure BDA0002638325200000181
其中,载体的第二平动速度矢量估计值由上述
Figure BDA0002638325200000182
确定。ω为载体的第一转动角速度矢量估计值,
Figure BDA0002638325200000183
为第i个第一传感器的归一化的平动速度矢量估计值,r1,i为第i个第一传感器相对于载体坐标系的位置平移矢量,si为第i个第一传感器的平动速度矢量的归一化参数,si具体可以根据下述关系式确定:
Figure BDA0002638325200000184
其中
Figure BDA0002638325200000185
Figure BDA0002638325200000186
为载体的第一平动速度矢量估计值。
进一步的,基于上述两种实现方式中任一种实现方式,可以通过迭代方式确定第二平动速度矢量估计值。具体的,每次迭代可以得到一次载体的平动速度矢量估计值,最后一次迭代得到的载体的平动速度矢量估计值作为载体的第二平动速度矢量估计值。迭代实现方式可以进一步利用传感器的平动速度矢量估计值和位置平移矢量,从而提高载体的平动速度矢量的估计精度。
具体的,在第一种迭代方式中,可以是根据如下关系式得到载体的第二平动速度矢量估计值:
Figure BDA0002638325200000187
其中,tk为第k次迭代中得到的载体的平动速度矢量估计值,也可以理解,第K次迭代(最后一次迭代)得到t作为载体的第二平动速度矢量估计值
Figure BDA0002638325200000188
ω为载体的第一转动角速度矢量估计值,
Figure BDA0002638325200000189
为第一传感器的归一化的平动速度矢量估计值,r1为第一传感器相对于载体坐标系的位置平移矢量,sk为第k次迭代中第一传感器的平动速度矢量的归一化参数。具体的,sk可以是根据如下关系式得到:
sk=‖tk-1+ω×r1
其中,tk-1为第k-1次迭代中得到的载体的平动速度矢量估计值,ω为第一转动角速度矢量估计值,r1为第一传感器相对于载体坐标系的位置平移矢量。
作为一种具体实现,每次迭代可以是,根据M′个第一传感器的归一化的平动速度矢量估计值、第k-1次迭代中M′个第一传感器的平动速度矢量的归一化参数、M′个第一传感器相对于载体坐标系的位置平移矢量和第一转动角速度矢量估计值,确定第k次迭代载体的平动速度矢量估计值。其中,第k次迭代中M′个第一传感器的平动速度矢量的归一化参数可以根据第k-1次迭代中得到的载体的平动速度矢量估计值确定。通过迭代过程,将第K次迭代(最后一轮迭代)中得到符合预设条件的载体的平动速度矢量估计值作为载体的第二平动速度矢量估计值。
作为一种实现方式,可以根据上述第一传感器的平动速度矢量、第一传感器的位置平移矢量与载体的平动速度矢量和载体的转动角速度矢量的关系,确定载体的第二平动速度矢量估计值,具体包括:
载体的第二平动速度矢量估计值根据以下关系式确定:
Figure BDA00026383252000001810
其中,tk为第k次迭代中载体的平动速度矢量估计值,tK为第二平动速度矢量估计值
Figure BDA00026383252000001811
ω为第一转动角速度矢量估计值,
Figure BDA00026383252000001812
为第i个第一传感器的归一化的平动速度矢量估计值,r1,i为第i个第一传感器相对于载体坐标系的位置平移矢量,si,k为第k次迭代中第i个第一传感器的平动速度矢量的归一化参数,si,k由第一平动速度矢量估计值或第k-1次迭代中载体的平动速度矢量估计值确定,w′1,i,k为第k次迭代中第i个第一传感器的加权系数。
w′1,i,k在每次迭代中可以取固定值或者根据预设算法确定。示例性的,w′1,i,k可以根据si·v1,i-ω×r1,i的估计误差的协方差矩阵确定,与之前方法类似,此处不详述。
或者,载体的第二平动速度矢量估计值
Figure BDA0002638325200000191
根据以下关系式确定:
Figure BDA0002638325200000192
其中,tk为第k次迭代中载体的平动速度矢量估计值,tK为第二平动速度矢量估计值
Figure BDA0002638325200000193
ω为第一转动角速度矢量估计值,
Figure BDA0002638325200000194
为第i个第一传感器的归一化的平动速度矢量估计值,r1,i为第i个第一传感器相对于载体坐标系的位置平移矢量,si,k为第k次迭代中第i个第一传感器的平动速度矢量的归一化参数,si,k由第一平动速度矢量估计值或第k-1次迭代中载体的平动速度矢量估计值确定。
si,k是根据第一平动速度矢量估计值或第k-1次迭代中载体的平动速度矢量估计值确定,可以有如下两种情况:
k等于1,也即第一次迭代中,si,1可以由第一平动速度矢量估计值根据如下关系式确定:
Figure BDA0002638325200000195
其中,
Figure BDA0002638325200000196
为第一平动速度矢量估计值,ω为第一转动角速度矢量估计值,r1,i为第i个第一传感器的坐标系相对于载体坐标系的位置平移矢量。
k大于1,也即之后的几次迭代中,si,k可以由第k-1次迭代中载体的平动速度矢量估计值根据如下关系式确定:
si,k=‖tk-1+ω×r1,i
其中,tk-1为第k-1次迭代中载体的平动速度矢量估计值,ω为第一转动角速度矢量估计值,r1,i为第i个第一传感器的坐标系相对于载体坐标系的位置平移矢量。
可选的,作为一种具体的实现方式,M′等于1,可以根据第一传感器的归一化的平动速度矢量估计值、第k-1次迭代中第一传感器的平动速度矢量的归一化参数、第一传感器的坐标系相对于载体坐标系的位置平移矢量和第一转动角速度矢量估计值,确定第k次迭代中载体的平动速度矢量估计值,具体可参照上述关系式得到。其中,第k次迭代中第一传感器的平动速度矢量的归一化参数是根据第一平动速度矢量估计值或第k-1次迭代中载体的平动速度矢量估计值确定。
需要说明的是,可以将第一种迭代方式的迭代终止条件设置为,第二平动速度矢量估计值和第一平动速度矢量估计值之间的矢量距离不大于第一预设的阈值或者门限。相当于,确定第k次迭代得到的载体的平动速度矢量估计值与第一平动速度矢量估计值之间的矢量距离,若该矢量距离大于第一预设的阈值或者门限,则进一步执行第k+1次迭代,若该矢量距离不大于第一预设的阈值或者门限,则确定迭代终止,此时,第k次迭代(也即第K次迭代,或最后一次迭代)得到的载体的平动速度矢量估计值可以称为是第二平动速度矢量估计值。
示例性的,第二平动速度矢量估计值和第一平动速度矢量估计值满足如下关系式:
Figure BDA0002638325200000197
其中,Threshold1为第一预设阈值或者门限。
另外,还可以将第一种迭代方式的迭代终止条件设置为,达到最大迭代次数。相当于,设置最大迭代次数为K,也即,共计执行K次迭代,将第K次迭代(也即最后一次迭代)得到的载体的平动速度矢量估计值称为是第二平动速度矢量估计值。示例性的,可以设置最大迭代次数K的取值为20。
上述第一种迭代方式中,每次迭代中,根据多个第一传感器的归一化的平动速度矢量估计值,确定载体的平动速度矢量估计值。将多个第一传感器的参数作为一次迭代的输入,有助于提高第二平动速度矢量的估计精度。
具体的,第二种迭代方式中,根据第i个第一传感器对应的载体的平动速度矢量估计值、第i个第一传感器的平动速度矢量的归一化参数、第i个第一传感器的坐标系相对于载体坐标系的位置平移矢量,确定第i个第一传感器对应的载体的平动速度矢量估计值。其中,第i个第一传感器的平动速度矢量的归一化参数是根据第i-1个第一传感器对应的载体的平动速度矢量估计值确定的。
等价地,第二种迭代方式为根据前一个第一传感器对应的载体的平动速度矢量估计值,确定后一个第一传感器的平动速度矢量的归一化参数,进而确定后一个第一传感器对应的载体的平动速度矢量估计值,其主要适用于M′大于1的情况。
可以是根据如下关系式得到的载体的平动速度矢量估计值:
Figure BDA0002638325200000201
其中,
Figure BDA0002638325200000202
为第i个第一传感器对应的载体的平动速度矢量估计值,
Figure BDA0002638325200000203
为载体的第二平动速度矢量估计值
Figure BDA0002638325200000204
ω为第一转动角速度矢量估计值,
Figure BDA0002638325200000205
为第i个第一传感器的归一化的平动速度矢量估计值,r1,i第i个第一传感器的坐标系相对于载体坐标系的位置平移矢量,si为第i个第一传感器的平动速度矢量的归一化参数,si由第一平动速度矢量估计值或
Figure BDA0002638325200000206
确定。
具体的,si由可以是根据如下关系式得到。
Figure BDA0002638325200000207
其中,
Figure BDA0002638325200000208
为第i-1个第一传感器对应的载体的平动速度矢量估计值,
Figure BDA0002638325200000209
为第一平动速度矢量估计值;ω为第一转动角速度矢量估计值,r1,i为第i个第一传感器的坐标系相对于载体坐标系的位置平移矢量。
此外,本申请实施例可以执行多轮迭代,每轮迭代包括M′个第一传感器之间的迭代。可以是根据如下关系式得到的载体的平动速度矢量估计值:
Figure BDA00026383252000002010
其中,
Figure BDA00026383252000002011
为第l轮迭代中第i个第一传感器对应的载体的平动速度矢量估计值,也可以理解,第L轮迭代(最后一轮迭代)中,ti,L′为第二平动速度矢量估计值
Figure BDA00026383252000002012
ω为第一转动角速度矢量估计值,
Figure BDA00026383252000002013
为第i个第一传感器的归一化的平动速度矢量估计值,r1,i第i个第一传感器的坐标系相对于载体坐标系的位置平移矢量,si,l为第l轮迭代中第i个第一传感器的平动速度矢量的归一化参数,si,l由第一平动速度矢量估计值或
Figure BDA00026383252000002014
Figure BDA00026383252000002015
确定。
具体的,si,l由可以是根据如下关系式得到。
Figure BDA00026383252000002016
其中,
Figure BDA0002638325200000211
为第l轮迭代中第i-1个第一传感器对应的载体的平动速度矢量估计值,
Figure BDA0002638325200000212
为第一平动速度矢量估计值;
Figure BDA0002638325200000213
ω为第一转动角速度矢量估计值,r1,i为第i个第一传感器的坐标系相对于载体坐标系的位置平移矢量。
需要说明的是,前一轮迭代中最后一次迭代是后一轮迭代中第一次迭代的前一次迭代,或者说,第l-1轮迭代中第M′次迭代是第l轮迭代中第1次迭代的前一次迭代。当然,也可以理解,第k次迭代可以是第l-1轮迭代中第M′次迭代,第k+1次迭代可以是第l轮迭代中第1次迭代。进一步的,后一轮迭代的参数可以由前一轮迭代的参数确定,具体的,第l轮迭代中第1个第一传感器的平动速度矢量的归一化参数由第l-1轮迭代中第M′个第一传感器对应的载体的平动速度矢量估计值确定。
si,l是根据第一平动速度矢量估计值,或第l轮迭代中第i-1个第一传感器对应的载体的平动速度矢量估计值,或第l-1轮迭代中第M′个第一传感器对应的载体的平动速度矢量估计值确定,可以有如下三种情况:
i等于1,l等于1,si,l是根据第一平动速度矢量估计值确定,可以参照如下关系式。
Figure BDA0002638325200000214
其中,
Figure BDA0002638325200000215
为第一平动速度矢量估计值,ω为第一转动角速度矢量估计值,r1,1为第1个第一传感器的坐标系相对于载体坐标系的位置平移矢量。
i大于1,si,l是根据第l轮迭代中第i-1个第一传感器对应的载体的平动速度矢量估计值确定,可以参照如下关系式。
Figure BDA0002638325200000216
其中,
Figure BDA0002638325200000217
为第l轮迭代中第i-1个第一传感器对应的载体的平动速度矢量估计值,ω为第一转动角速度矢量估计值,r1,i为第i个第一传感器的坐标系相对于载体坐标系的位置平移矢量。
i等于1,l大于1,s1,l是根据第l-1轮迭代中第M′个第一传感器对应的载体的平动速度矢量估计值确定,可以参照如下关系式。
Figure BDA0002638325200000218
其中,
Figure BDA0002638325200000219
为第l-1轮迭代中第M′个第一传感器对应的载体的平动速度矢量估计值,ω为第一转动角速度矢量估计值,r1,i为第i个第一传感器的坐标系相对于载体坐标系的位置平移矢量。
需要说明的是,可以将第二种迭代方式的迭代终止条件设置为,第二平动速度矢量估计值和第一平动速度矢量估计值之间的矢量距离不大于第二预设阈值或门限。相当于,确定第l轮迭代中第i-1个第一传感器对应的载体的平动速度矢量估计值与第一平动速度矢量估计值之间的矢量距离,若该矢量距离大于第二预设阈值或门限,则进一步确定第l轮迭代中第i个第一传感器对应的载体的平动速度矢量估计值,若该矢量距离不大于第二预设阈值或门限,则确定迭代终止,此时,可以将第l轮迭代中第i-1个第一传感器对应的载体的平动速度矢量估计值称为是第二平动速度矢量估计值。
示例性的,第二平动速度矢量估计值和第一平动速度矢量估计值满足如下关系式:
Figure BDA00026383252000002110
其中,Threshold2为第二预设阈值或门限。
另外,还可以将第二种迭代方式的迭代终止条件设置为,达到最大迭代轮数。相当于,设置最大迭代轮数为L,也即,共计执行L轮迭代,将第L轮迭代(也即最后一轮迭代)中某个第一传感器对应的载体的平动速度矢量估计值称为是第二平动速度矢量估计值。示例性的,设置最大迭代轮数L的取值为5,设置第5轮中第M′个第一传感器对应的载体的平动速度矢量估计值为第二平动速度矢量估计值。
上述第二种迭代方式中,针对每个第一传感器的归一化的平动速度矢量估计值,确定每个第一传感器对应的载体的平动速度矢量估计值,相当于,针对每个第一传感器进行一次迭代,将一个第一传感器的参数作为一次迭代的输入,可以快速得到较为准确的第二平动速度矢量估计值。
本申请实施例中,第一传感器评估自身运动得到的归一化的平动速度矢量估计值较准确,该将第一传感器的归一化的平动速度矢量估计值和第一平动速度矢量估计值进行融合,得到载体的第二平动速度矢量估计值,相比于第一平动速度矢量估计值,可进一步提高载体的平动速度矢量估计值的精确度。采用载体的第一转动角速度矢量估计值和载体的第二平动速度矢量估计值对载体的运动进行补偿,有助于实现运动目标和静止目标的分离,同时有助于实现对载体运动的定位和跟踪。
本申请实施例中,M'个第一传感器中的第一传感器可以理解是可以获得归一化的平动速度矢量估计值的第一传感器,M'个第一传感器可以是视觉传感器,例如相机、摄像头、红外传感器等。一个示例中,M'个第一传感器可以是M1'个相机和M2'个红外传感器,其中,M1'+M2'=M',且M1'≥0,M2'≥0。
此外,本申请实施例还可以进一步结合M-M'个可以获取平动速度矢量估计值的第一传感器,其中M-M′个第一传感器获取的平动速度矢量估计值包含速度矢量的尺度信息,例如各个速度分量的完整信息,而不仅仅是归一化值或者方向信息,从而提高载体的第二平动速度矢量估计值的准确性。
也即,本申请提供再一种确定第二平动速度矢量估计值的方式如下。
在步骤203中可以是获取M′个第一传感器如视觉传感器的归一化的平动速度矢量估计值和M-M′个第一传感器如惯性传感器(如IMU)的平动速度矢量估计值。
在步骤204中,先根据M-M′个第一传感器如惯性传感器的平动速度矢量估计值,确定M-M′个第一传感器如惯性传感器对应的载体的平动速度矢量估计值,
示例性的,可基于关系式
Figure BDA0002638325200000221
Figure BDA0002638325200000222
确定,关系式的具体描述可参照上述实现方式。
作为一种直接实现方式,可以基于关系式
Figure BDA0002638325200000223
Figure BDA0002638325200000224
Figure BDA0002638325200000225
确定视觉传感器对应的载体的平动速度矢量估计值,将视觉传感器对应的载体的平动速度矢量估计值、惯性传感器对应的载体的平动速度矢量估计值进行加权,得到第二平动速度矢量估计值。
作为一种迭代的实现方式,可以将惯性传感器对应的载体的平动速度矢量估计值作为第一次迭代中的载体的平动速度矢量估计值,也即,k等于1中对应的
Figure BDA0002638325200000226
中,
Figure BDA0002638325200000227
替换为惯性传感器对应的载体的平动速度矢量估计值。基于步骤204中的迭代方式进行迭代,直至符合迭代条件,将符合迭代条件的第k次迭代中的载体的平动速度矢量估计值作为第二平动速度矢量估计值。
进一步的,在本申请中,还可能存在M′=0,也即,在步骤203中获取M个惯性传感器的平动速度矢量估计值。
相应的,在步骤204中,可以根据M个惯性传感器的平动速度矢量估计值,确定惯性传感器对应的载体的平动速度矢量估计值,示例性的,可基于关系式
Figure BDA0002638325200000231
Figure BDA0002638325200000232
Figure BDA0002638325200000233
确定惯性传感器对应的载体的平动速度矢量估计值,关系式的具体描述可参照上述实现方式。将惯性传感器对应的载体的平动速度矢量估计值和第一平动速度矢量估计值进行加权得到第二平动速度矢量估计值。
可选地,根据上述第二平动速度矢量估计值,可以进一步更新载体的转动角速度矢量估计值。此处将更新的该转动角速度矢量估计值称为是载体的第二转动角速度矢量估计值。
可选的,步骤205,根据载体的第二平动速度矢量估计值,确定载体的第二转动角速度矢量估计值。
可以基于如下关系式,确定第二转动角速度矢量估计值:
Figure BDA0002638325200000234
或者,等价地
Figure BDA0002638325200000235
或者,等价地
Figure BDA0002638325200000236
其中,ω′为第二转动角速度矢量估计值,v2为第二传感器的瞬时速度矢量估计值,
Figure BDA0002638325200000237
为第二平动速度矢量估计值,r2为第二传感器的坐标系相对于载体坐标系的位置平移矢量,×表示矢量的叉积。其中[r2]×为与r2对应的反对称矩阵。
基于上述关系式,第二转动角速度矢量估计值的表达式可以为
Figure BDA0002638325200000238
作为一种实现方式,可以基于如下关系式,确定第二转动角速度矢量估计值:
Figure BDA0002638325200000239
其中,ω′为第二转动角速度矢量估计值,v2,j为第j个第二传感器的瞬时速度矢量估计值,
Figure BDA00026383252000002310
为第二平动速度矢量估计值,
Figure BDA00026383252000002311
为[r2,j]×的逆矩阵,[r2,j]×为与r2,j对应的反对称矩阵,r2,j为第j个第二传感器的坐标系相对于载体坐标系的位置平移矢量。
其中,r2,j=[rx,2,j ry,2,j rz,2,j]T
Figure BDA00026383252000002312
作为一种实现方式,根据第二转动角速度矢量估计值、N个第二传感器的瞬时速度矢量估计值、N个第二传感器中各第二传感器的坐标系相对于载体坐标系的位置平移矢量,确定第二转动角速度矢量估计值,其中N≥1。
在一个示例中,基于如下关系式确定第二转动角速度矢量估计值:
Figure BDA00026383252000002313
其中,ω′为第二转动角速度矢量估计值,w2,j为第j个第二传感器的加权系数,v2,j为第j个第二传感器的瞬时速度矢量估计值,
Figure BDA00026383252000002314
为第二平动速度矢量估计值,
Figure BDA00026383252000002315
为[r2,j]×的逆矩阵,[r2,j]×为与r2,j对应的反对称矩阵,r2,j为第j个第二传感器的坐标系相对于载体坐标系的位置平移矢量。
在又一个示例中,基于如下关系式确定第一转动角速度矢量估计值:
Figure BDA00026383252000002316
其中,ω′为第二转动角速度矢量估计值,v2,j为第j个第二传感器的瞬时速度矢量估计值,
Figure BDA0002638325200000241
为第二平动速度矢量估计值,
Figure BDA0002638325200000242
为[r2,j]×的逆矩阵,[r2,j]×为与r2,j对应的反对称矩阵,r2,j为第j个第二传感器的坐标系相对于载体坐标系的位置平移矢量。
作为另一种实现方式中,N等于1,根据第二转动角速度矢量估计值、第二传感器的瞬时速度矢量估计值、第二传感器的坐标系相对于载体坐标系的位置平移矢量,确定第二转动角速度矢量估计值,具体可参照上述关系式得到。
本申请中,可以将载体的第二平动速度矢量估计值和载体的第二转动角速度矢量估计值,分别作为最终确定出的载体的平动速度矢量估计值和载体的转动角速度矢量估计值。
此外,还可以将第一转动角速度矢量估计值和第二转动角速度矢量估计值进一步融合,将融合后的转动角速度矢量估计值作为最终确定的载体的转动角速度矢量估计值。具体融合中,可以是确定载体的第一转动角速度矢量估计值和载体的第二转动角速度矢量估计值的平均值。或者,基于最小均方误差进行加权组合,将载体的第一转动角速度矢量估计值和载体的第二转动角速度矢量估计值进行融合。
本申请实施例中,根据载体的第二平动速度矢量估计值,确定载体的第二转动角速度矢量估计值,相比于第一转动角速度矢量估计值,可进一步提高载体的转动角速度矢量估计值的精确度。采用载体的第二转动角速度矢量估计值和载体的第二平动速度矢量估计值对载体的运动进行补偿,有助于实现运动目标和静止目标的分离,同时有助于实现对载体运动的定位和跟踪。
本文中描述的各个实施例可以为独立的方案,也可以根据内在逻辑进行组合,这些方案都落入本申请的保护范围中。
本申请实施例中对模块的划分是示意性的,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式。另外,在本申请各个实施例中的各功能模块可以集成在一个处理器中,也可以是单独物理存在,也可以两个或两个以上模块集成在一个模块中。上述集成的模块既可以采用硬件的形式实现,也可以采用软件功能模块的形式实现。
与上述构思相同,本申请实施例还提供一种运动估计装置用于实现上述方法。
示例性的,该运动估计装置可以是传感器系统或者融合感知系统或者集成上述系统的规划/控制系统如辅助驾驶或者自动驾驶系统等。或者,该运动估计装置还可以是软件或者硬件(如与相应传感器通过无线或者有线连接或者集成在一起的数据处理装置)。
该运动估计装置可为具有运动估计功能的车辆,或者为具有运动估计功能的其他部件。该运动估计装置包括但不限于:车载终端、车载控制器、车载模块、车载模组、车载部件、车载芯片、车载单元、车载雷达或车载摄像头等其他传感器,车辆可通过该车载终端、车载控制器、车载模块、车载模组、车载部件、车载芯片、车载单元、车载雷达或摄像头,实施本申请提供的方法。
该运动估计装置还可以为除了车辆之外的其他具有运动估计功能的智能终端,或设置在除了车辆之外的其他具有运动估计功能的智能终端中,或设置于该智能终端的部件中。该智能终端可以为智能运输设备、智能家居设备、机器人等其他终端设备。该运动估计装置包括但不限于智能终端或智能终端内的控制器、芯片、雷达或摄像头等其他传感器、以及其他部件等。
该运动估计装置可以是一个通用设备或者是一个专用设备。在具体实现中,该运动估计装置还可以台式机、便携式电脑、网络服务器、掌上电脑(personal digitalassistant,PDA)、移动手机、平板电脑、无线终端设备、嵌入式设备或其他具有处理功能的设备。本申请实施例不限定该运动估计装置的类型。
该运动估计装置还可以是具有处理功能的芯片或处理器,该运动估计装置可以包括多个处理器。处理器可以是一个单核(single-CPU)处理器,也可以是一个多核(multi-CPU)处理器。该具有处理功能的芯片或处理器可以设置在传感器中,也可以不设置在传感器中,而设置在传感器输出信号的接收端。
如图7为本申请示例性提供的一种运动估计装置700。该运动估计装置700可以包括:获取单元701和处理单元702。应理解,装置实施例的描述与方法实施例的描述相互对应,因此,未详细描述的内容可以参见上文方法实施例,为了简洁,这里不再赘述。
示例性的,所述获取单元701用于获取M个第一传感器的转动角速度矢量估计值和N个第二传感器的瞬时速度矢量估计值;其中,M≥1,N≥1;
所述处理单元702用于根据所述N个第二传感器的瞬时速度矢量估计值和所述N个第二传感器所在的载体的第一转动角速度矢量估计值,确定所述载体的第一平动速度矢量估计值,其中,所述第一转动角速度矢量估计值是根据所述M个第一传感器的转动角速度矢量估计值确定的。
在一种可能的实现方式中,所述处理单元702具体用于基于如下关系式确定所述第一平动速度矢量估计值:
Figure BDA0002638325200000251
其中,
Figure BDA0002638325200000252
为第一平动速度矢量估计值,ω为第一转动角速度矢量估计值,v2,j为第j个第二传感器的瞬时速度矢量估计值,r2,j为第j个第二传感器的坐标系相对于载体坐标系的位置平移矢量。
在一种可能的实现方式中,所述第一平动速度矢量估计值满足如下关系式:
Figure BDA0002638325200000253
其中,
Figure BDA0002638325200000254
为第一平动速度矢量估计值,ω为第一转动角速度矢量估计值,w2,j为第j个第二传感器的加权系数,w2,j为第j个第二传感器的瞬时速度矢量估计值,r2,j为第j个第二传感器的坐标系相对于载体坐标系的位置平移矢量。
在一种可能的实现方式中,所述第一转动角速度矢量估计值满足如下关系式:
Figure BDA0002638325200000255
其中,ω为第一转动角速度矢量估计值,w1,i为第i个第一传感器的加权系数,ω1,i为第i个第一传感器的转动角速度矢量估计值。
在一种可能的实现方式中,所述获取单元701还用于获取所述M个第一传感器中M′个第一传感器的归一化的平动速度矢量估计值,其中,1≤M′≤M;所述处理单元702还用于根据所述第一平动速度矢量估计值、所述M′个第一传感器的归一化的平动速度矢量估计值,确定所述载体的第二平动速度矢量估计值。
在一种可能的实现方式中,所述处理单元702具体用于基于如下关系式确定所述第二平动速度矢量估计值:
Figure BDA0002638325200000261
其中,
Figure BDA0002638325200000262
为第二平动速度矢量估计值,ω为第一转动角速度矢量估计值,
Figure BDA0002638325200000263
为第i个第一传感器的归一化的平动速度矢量估计值,r1,i为第i个第一传感器的坐标系相对于载体坐标系的位置平移矢量,si为第i个第一传感器的平动速度矢量的归一化参数或者尺度伸缩因子,si由第一平动速度矢量估计值确定。
在一种可能的实现方式中,所述第二平动速度矢量估计值满足如下关系式:
Figure BDA0002638325200000264
其中,tK为第二平动速度矢量估计值
Figure BDA0002638325200000265
tk为第k次迭代中载体的平动速度矢量估计值,ω为第一转动角速度矢量估计值,
Figure BDA0002638325200000266
为第i个第一传感器的归一化的平动速度矢量估计值,r1,i为第i个第一传感器的坐标系相对于载体坐标系的位置平移矢量,si,k为第k次迭代中第i个第一传感器的平动速度矢量的归一化参数或者尺度伸缩因子,si,k由第一平动速度矢量估计值或第k-1次迭代中载体的平动速度矢量估计值确定。
在一种可能的实现方式中,所述第二平动速度矢量估计值满足如下关系式:
Figure BDA0002638325200000267
其中,tK为第二平动速度矢量估计值
Figure BDA0002638325200000268
tk为第k次迭代中载体的平动速度矢量估计值,w′1,i,k为第k次迭代中第i个第一传感器的加权系数,ω为第一转动角速度矢量估计值,
Figure BDA0002638325200000269
为第i个第一传感器的归一化的平动速度矢量估计值,r1,i为第i个第一传感器的坐标系相对于载体坐标系的位置平移矢量,si,k为第k次迭代中第i个第一传感器的平动速度矢量的归一化参数或者尺度伸缩因子,si,k由第一平动速度矢量估计值或第k-1次迭代中载体的平动速度矢量估计值确定。
在一种可能的实现方式中,第k次迭代中第i个第一传感器的平动速度矢量的归一化参数或者尺度伸缩因子满足如下关系式:
si,k=‖tk-1+ω×r1,i
其中,tk-1为第k-1次迭代中载体的平动速度矢量估计值,t0为第一平动速度矢量估计值。
在一种可能的实现方式中,所述第二平动速度矢量估计值满足如下关系式:
Figure BDA00026383252000002610
其中,
Figure BDA00026383252000002611
为第二平动速度矢量估计值
Figure BDA00026383252000002612
Figure BDA00026383252000002613
为第l轮迭代中第i个第一传感器对应的载体的平动速度矢量估计值,ω为第一转动角速度矢量估计值,
Figure BDA00026383252000002614
为第i个第一传感器的归一化的平动速度矢量估计值,r1,i第i个第一传感器的坐标系相对于载体坐标系的位置平移矢量,si,l为第l轮迭代中第i个第一传感器的平动速度矢量的归一化参数或者尺度伸缩因子,si,l由第一平动速度矢量估计值或
Figure BDA00026383252000002615
Figure BDA00026383252000002616
确定。
在一种可能的实现方式中,所述第l轮迭代中第i个第一传感器的平动速度矢量的归一化参数或者尺度伸缩因子满足如下关系式:
Figure BDA00026383252000002617
其中,
Figure BDA00026383252000002618
为第一平动速度矢量估计值;
Figure BDA00026383252000002619
在一种可能的实现方式中,所述处理单元702还用于基于如下关系式确定所述载体的第二转动角速度矢量估计值:
Figure BDA0002638325200000271
其中,ω′为第二转动角速度矢量估计值,v2,j为第j个第二传感器的瞬时速度矢量估计值,
Figure BDA0002638325200000272
为第二平动速度矢量估计值,
Figure BDA0002638325200000273
为[r2,j]×的逆矩阵,[r2,j]×为与r2,j对应的反对称矩阵,r2,j为第j个第二传感器的坐标系相对于载体坐标系的位置平移矢量。
与上述构思相同,如图8所示,为本申请实施例提供的一种芯片的结构示意图。
芯片800包括一个或多个处理器801以及接口电路802。可选的,所述芯片800还可以包含总线803。
其中,处理器801可能是一种集成电路芯片,具有信号的处理能力。在实现过程中,上述方法的各步骤可以通过处理器801中的硬件的集成逻辑电路或者软件形式的指令完成。上述的处理器801可以是通用处理器、数字通信器(DSP)、专用集成电路(ASIC)、现场可编程门阵列(FPGA)或者其它可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本申请实施例中的公开的各方法、步骤。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。
接口电路802可以用于数据、指令或者信息的发送或者接收,处理器801可以利用接口电路802接收的数据、指令或者其它信息,进行加工,可以将加工完成信息通过接口电路802发送出去。
可选的,芯片还包括存储器,存储器可以包括只读存储器和随机存取存储器,并向处理器提供操作指令和数据。存储器的一部分还可以包括非易失性随机存取存储器(NVRAM)。
可选的,存储器存储了可执行软件模块或者数据结构,处理器可以通过调用存储器存储的操作指令(该操作指令可存储在操作系统中),执行相应的操作。
可选的,芯片可以使用在本申请实施例涉及的通信装置(包括主节点和从节点)中。可选的,接口电路802可用于输出处理器801的执行结果。关于本申请的一个或多个实施例提供的数据传输方法可参考前述各个实施例,这里不再赘述。
需要说明的,处理器801、接口电路802各自对应的功能既可以通过硬件设计实现,也可以通过软件设计来实现,还可以通过软硬件结合的方式来实现,这里不作限制。
本申请实施例还提供一种雷达系统,用于为车辆提供运动估计功能。其包含至少一个本申请上述实施例提到的运动估计装置,该系统内的至少一个运动估计装置可以集成为一个整机或设备,或者该系统内的至少一个运动估计装置也可以独立设置为元件或装置。
本申请实施例还提供一种传感器系统,用于为车辆提供运动估计功能。其包含至少一个本申请上述实施例提到的运动估计装置,以及,摄像头或雷达等传感器中的至少一个,该系统内的至少一个传感器装置可以集成为一个整机或设备,或者该系统内的至少一个传感器装置也可以独立设置为元件或装置。
本申请实施例还提供一种系统,应用于无人驾驶或智能驾驶中,其包含至少一个本申请上述实施例提到的运动估计装置、摄像头、雷达等传感器中的至少一个,该系统内的至少一个装置可以集成为一个整机或设备,或者该系统内的至少一个装置也可以独立设置为元件或装置。
进一步,上述任一系统可以与车辆的中央控制器进行交互,为所述车辆驾驶的决策或控制提供探测和/或融合信息。
本申请实施例还提供一种车辆,所述车辆包括至少一个本申请上述实施例提到的运动估计装置或上述任一系统。
本申请实施例还提供一种通信装置,包括处理器和通信接口,所述通信接口用于接收来自所述通信装置之外的其它通信装置的信号并传输至所述处理器或将来自所述处理器的信号发送给所述通信装置之外的其它通信装置,所述处理器通过逻辑电路或执行代码指令用于实现上述如图2中的方法。
本申请实施例还一种计算机可读存储介质,所述计算机可读存储介质中存储有计算机程序或指令,当所述计算机程序或指令被通信装置执行时,实现上述如图2中的方法。
本申请实施例还提供一种计算机程序产品,所述计算机程序产品包括计算机程序或指令,当所述计算机程序或指令被通信装置执行时,实现上述如图2中的方法。
本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
显然,本领域的技术人员可以对本申请进行各种改动和变型而不脱离本申请的范围。这样,倘若本申请的这些修改和变型属于本申请权利要求及其等同技术的范围之内,则本申请也意图包含这些改动和变型在内。

Claims (26)

1.一种运动估计的方法,其特征在于,包括:
获取M个第一传感器的转动角速度矢量估计值和N个第二传感器的瞬时速度矢量估计值;其中,M≥1,N≥1;
根据所述N个第二传感器的瞬时速度矢量估计值和所述N个第二传感器所在的载体的第一转动角速度矢量估计值,确定所述载体的第一平动速度矢量估计值,其中,所述第一转动角速度矢量估计值是根据所述M个第一传感器的转动角速度矢量估计值确定的。
2.如权利要求1所述的方法,其特征在于,所述第一平动速度矢量估计值基于如下关系式确定:
Figure FDA0002638325190000011
其中,
Figure FDA0002638325190000012
为第一平动速度矢量估计值,ω为第一转动角速度矢量估计值,v2,j为第j个第二传感器的瞬时速度矢量估计值,r2,j为第j个第二传感器的坐标系相对于载体坐标系的位置平移矢量。
3.如权利要求1所述的方法,其特征在于,所述第一平动速度矢量估计值满足如下关系式:
Figure FDA0002638325190000013
其中,
Figure FDA0002638325190000014
为第一平动速度矢量估计值,ω为第一转动角速度矢量估计值,w2,j为第j个第二传感器的加权系数,v2,j为第j个第二传感器的瞬时速度矢量估计值,r2,j为第j个第二传感器的坐标系相对于载体坐标系的位置平移矢量。
4.如权利要求1至3任一项所述的方法,其特征在于,所述第一转动角速度矢量估计值满足如下关系式:
Figure FDA0002638325190000015
其中,ω为第一转动角速度矢量估计值,w1,i为第i个第一传感器的加权系数,ω1,i为第i个第一传感器的转动角速度矢量估计值。
5.如权利要求1至4任一项所述的方法,其特征在于,还包括:
获取所述M个第一传感器中M′个第一传感器的归一化的平动速度矢量估计值,其中,1≤M′≤M;
根据所述第一平动速度矢量估计值、所述M′个第一传感器的归一化的平动速度矢量估计值,确定所述载体的第二平动速度矢量估计值。
6.如权利要求5所述的方法,其特征在于,所述第二平动速度矢量估计值基于如下关系式确定:
Figure FDA0002638325190000016
其中,
Figure FDA0002638325190000017
为第二平动速度矢量估计值,ω为第一转动角速度矢量估计值,
Figure FDA0002638325190000018
为第i个第一传感器的归一化的平动速度矢量估计值,r1,i为第i个第一传感器的坐标系相对于载体坐标系的位置平移矢量,si为第i个第一传感器的平动速度矢量的归一化参数或者尺度伸缩因子,si由第一平动速度矢量估计值确定。
7.如权利要求5所述的方法,其特征在于,所述第二平动速度矢量估计值满足如下关系式:
Figure FDA0002638325190000021
其中,tK为第二平动速度矢量估计值
Figure FDA0002638325190000022
tk为第k次迭代中载体的平动速度矢量估计值,ω为第一转动角速度矢量估计值,
Figure FDA0002638325190000023
为第i个第一传感器的归一化的平动速度矢量估计值,r1,i为第i个第一传感器的坐标系相对于载体坐标系的位置平移矢量,si,k为第k次迭代中第i个第一传感器的平动速度矢量的归一化参数或者尺度伸缩因子,si,k由第一平动速度矢量估计值或第k-1次迭代中载体的平动速度矢量估计值确定。
8.如权利要求5所述的方法,其特征在于,所述第二平动速度矢量估计值满足如下关系式:
Figure FDA0002638325190000024
其中,tK为第二平动速度矢量估计值
Figure FDA0002638325190000025
tk为第k次迭代中载体的平动速度矢量估计值,w′1,i,k为第k次迭代中第i个第一传感器的加权系数,ω为第一转动角速度矢量估计值,
Figure FDA0002638325190000026
为第i个第一传感器的归一化的平动速度矢量估计值,r1,i为第i个第一传感器的坐标系相对于载体坐标系的位置平移矢量,si,k为第k次迭代中第i个第一传感器的平动速度矢量的归一化参数或者尺度伸缩因子,si,k由第一平动速度矢量估计值或第k-1次迭代中载体的平动速度矢量估计值确定。
9.如权利要求7或8所述的方法,其特征在于,第k次迭代中第i个第一传感器的平动速度矢量的归一化参数或者尺度伸缩因子满足如下关系式:
si,k=||tk-1+ω×r1,i||
其中,tk-1为第k-1次迭代中载体的平动速度矢量估计值,t0为第一平动速度矢量估计值。
10.如权利要求5所述的方法,其特征在于,所述第二平动速度矢量估计值满足如下关系式:
Figure FDA0002638325190000027
其中,
Figure FDA0002638325190000028
为第二平动速度矢量估计值
Figure FDA0002638325190000029
Figure FDA00026383251900000210
为第l轮迭代中第i个第一传感器对应的载体的平动速度矢量估计值,ω为第一转动角速度矢量估计值,
Figure FDA00026383251900000211
为第i个第一传感器的归一化的平动速度矢量估计值,r1,i第i个第一传感器的坐标系相对于载体坐标系的位置平移矢量,si,l为第l轮迭代中第i个第一传感器的平动速度矢量的归一化参数或者尺度伸缩因子,si,l由第一平动速度矢量估计值或
Figure FDA00026383251900000212
Figure FDA00026383251900000213
确定。
11.如权利要求10所述的方法,其特征在于,所述第l轮迭代中第i个第一传感器的平动速度矢量的归一化参数或者尺度伸缩因子满足如下关系式:
Figure FDA00026383251900000214
其中,
Figure FDA00026383251900000215
为第一平动速度矢量估计值;
Figure FDA00026383251900000216
12.如权利要求5至11任一项所述的方法,其特征在于,还包括:
基于如下关系式确定所述载体的第二转动角速度矢量估计值:
Figure FDA0002638325190000031
其中,ω′为第二转动角速度矢量估计值,v2,j为第j个第二传感器的瞬时速度矢量估计值,
Figure FDA0002638325190000032
为第二平动速度矢量估计值,
Figure FDA0002638325190000033
为[r2,j]×的逆矩阵,[r2,j]×为与r2,j对应的反对称矩阵,r2,j为第j个第二传感器的坐标系相对于载体坐标系的位置平移矢量。
13.一种运动估计的装置,其特征在于,包括:
获取单元和处理单元;
所述获取单元用于获取M个第一传感器的转动角速度矢量估计值和N个第二传感器的瞬时速度矢量估计值;其中,M≥1,N≥1;
所述处理单元用于根据所述N个第二传感器的瞬时速度矢量估计值和所述N个第二传感器所在的载体的第一转动角速度矢量估计值,确定所述载体的第一平动速度矢量估计值,其中,所述第一转动角速度矢量估计值是根据所述M个第一传感器的转动角速度矢量估计值确定的。
14.如权利要求13所述的装置,其特征在于,所述处理单元具体用于基于如下关系式确定所述第一平动速度矢量估计值:
Figure FDA0002638325190000034
其中,
Figure FDA0002638325190000035
为第一平动速度矢量估计值,ω为第一转动角速度矢量估计值,v2,j为第j个第二传感器的瞬时速度矢量估计值,r2,j为第j个第二传感器的坐标系相对于载体坐标系的位置平移矢量。
15.如权利要求13所述的装置,其特征在于,所述第一平动速度矢量估计值满足如下关系式:
Figure FDA0002638325190000036
其中,
Figure FDA0002638325190000037
为第一平动速度矢量估计值,ω为第一转动角速度矢量估计值,w2,j为第j个第二传感器的加权系数,v2,j为第j个第二传感器的瞬时速度矢量估计值,r2,j为第j个第二传感器的坐标系相对于载体坐标系的位置平移矢量。
16.如权利要求13至15任一项所述的装置,其特征在于,所述第一转动角速度矢量估计值满足如下关系式:
Figure FDA0002638325190000038
其中,ω为第一转动角速度矢量估计值,w1,i为第i个第一传感器的加权系数,ω1,i为第i个第一传感器的转动角速度矢量估计值。
17.如权利要求13至16任一项所述的装置,其特征在于,所述获取单元还用于获取所述M个第一传感器中M′个第一传感器的归一化的平动速度矢量估计值,其中,1≤M′≤M;所述处理单元还用于根据所述第一平动速度矢量估计值、所述M′个第一传感器的归一化的平动速度矢量估计值,确定所述载体的第二平动速度矢量估计值。
18.如权利要求17所述的装置,其特征在于,所述处理单元具体用于基于如下关系式确定所述第二平动速度矢量估计值:
Figure FDA0002638325190000041
其中,
Figure FDA0002638325190000042
为第二平动速度矢量估计值,ω为第一转动角速度矢量估计值,
Figure FDA0002638325190000043
为第i个第一传感器的归一化的平动速度矢量估计值,r1,i为第i个第一传感器的坐标系相对于载体坐标系的位置平移矢量,si为第i个第一传感器的平动速度矢量的归一化参数或者尺度伸缩因子,si由第一平动速度矢量估计值确定。
19.如权利要求17所述的装置,其特征在于,所述第二平动速度矢量估计值满足如下关系式:
Figure FDA0002638325190000044
其中,tk为第二平动速度矢量估计值
Figure FDA0002638325190000045
tk为第k次迭代中载体的平动速度矢量估计值,ω为第一转动角速度矢量估计值,
Figure FDA0002638325190000046
为第i个第一传感器的归一化的平动速度矢量估计值,r1,i为第i个第一传感器的坐标系相对于载体坐标系的位置平移矢量,si,k为第k次迭代中第i个第一传感器的平动速度矢量的归一化参数或者尺度伸缩因子,si,k由第一平动速度矢量估计值或第k-1次迭代中载体的平动速度矢量估计值确定。
20.如权利要求17所述的装置,其特征在于,所述第二平动速度矢量估计值满足如下关系式:
Figure FDA0002638325190000047
其中,tK为第二平动速度矢量估计值
Figure FDA0002638325190000048
tk为第k次迭代中载体的平动速度矢量估计值,w′1,i,k为第k次迭代中第i个第一传感器的加权系数,ω为第一转动角速度矢量估计值,
Figure FDA0002638325190000049
为第i个第一传感器的归一化的平动速度矢量估计值,r1,i为第i个第一传感器的坐标系相对于载体坐标系的位置平移矢量,si,k为第k次迭代中第i个第一传感器的平动速度矢量的归一化参数或者尺度伸缩因子,si,k由第一平动速度矢量估计值或第k-1次迭代中载体的平动速度矢量估计值确定。
21.如权利要求19或20所述的装置,其特征在于,第k次迭代中第i个第一传感器的平动速度矢量的归一化参数或者尺度伸缩因子满足如下关系式:
si,k=||tk-1+ω×r1,i||
其中,tk-1为第k-1次迭代中载体的平动速度矢量估计值,t0为第一平动速度矢量估计值。
22.如权利要求17所述的装置,其特征在于,所述第二平动速度矢量估计值满足如下关系式:
Figure FDA00026383251900000410
其中,
Figure FDA00026383251900000411
为第二平动速度矢量估计值
Figure FDA00026383251900000412
Figure FDA00026383251900000413
为第l轮迭代中第i个第一传感器对应的载体的平动速度矢量估计值,ω为第一转动角速度矢量估计值,
Figure FDA00026383251900000414
为第i个第一传感器的归一化的平动速度矢量估计值,r1,i第i个第一传感器的坐标系相对于载体坐标系的位置平移矢量,si,l为第l轮迭代中第i个第一传感器的平动速度矢量的归一化参数或者尺度伸缩因子,si,l由第一平动速度矢量估计值或
Figure FDA00026383251900000415
Figure FDA00026383251900000416
确定。
23.如权利要求22所述的装置,其特征在于,所述第l轮迭代中第i个第一传感器的平动速度矢量的归一化参数或者尺度伸缩因子满足如下关系式:
Figure FDA00026383251900000417
其中,
Figure FDA0002638325190000051
为第一平动速度矢量估计值;
Figure FDA0002638325190000052
24.如权利要求17至23任一项所述的装置,其特征在于,所述处理单元还用于基于如下关系式确定所述载体的第二转动角速度矢量估计值:
Figure FDA0002638325190000053
其中,ω′为第二转动角速度矢量估计值,v2,j为第j个第二传感器的瞬时速度矢量估计值,
Figure FDA0002638325190000054
为第二平动速度矢量估计值,
Figure FDA0002638325190000055
为[r2,j]×的逆矩阵,[r2,j]×为与r2,j对应的反对称矩阵,r2,j为第j个第二传感器的坐标系相对于载体坐标系的位置平移矢量。
25.一种芯片,其特征在于,包括至少一个处理器和接口;
所述接口,用于为所述至少一个处理器提供程序指令或者数据;
所述至少一个处理器用于执行所述程序行指令,以实现如权利要求1至12中任一项所述的方法。
26.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质中存储有计算机程序或指令,当所述计算机程序或指令被通信装置执行时,实现如权利要求1至12中任一项所述的方法。
CN202010831971.6A 2020-08-18 2020-08-18 一种运动估计方法及装置 Pending CN114076946A (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202010831971.6A CN114076946A (zh) 2020-08-18 2020-08-18 一种运动估计方法及装置
PCT/CN2021/108627 WO2022037370A1 (zh) 2020-08-18 2021-07-27 一种运动估计方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010831971.6A CN114076946A (zh) 2020-08-18 2020-08-18 一种运动估计方法及装置

Publications (1)

Publication Number Publication Date
CN114076946A true CN114076946A (zh) 2022-02-22

Family

ID=80281546

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010831971.6A Pending CN114076946A (zh) 2020-08-18 2020-08-18 一种运动估计方法及装置

Country Status (2)

Country Link
CN (1) CN114076946A (zh)
WO (1) WO2022037370A1 (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115575923B (zh) * 2022-12-08 2023-03-28 千巡科技(深圳)有限公司 一种地面机器人静止判断方法、系统、装置以及存储介质

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3862015B2 (ja) * 2002-10-25 2006-12-27 オムロン株式会社 車載用レーダ装置
US8855848B2 (en) * 2007-06-05 2014-10-07 GM Global Technology Operations LLC Radar, lidar and camera enhanced methods for vehicle dynamics estimation
JP2010085196A (ja) * 2008-09-30 2010-04-15 Toyota Motor Corp 物体挙動推定方法および物体挙動推定装置
US10094919B2 (en) * 2015-10-06 2018-10-09 GM Global Technology Operations LLC Radar-vision fusion for target velocity estimation
EP3252501B1 (en) * 2016-06-03 2022-06-08 Veoneer Sweden AB Enhanced object detection and motion state estimation for a vehicle environment detection system
US10599931B2 (en) * 2017-08-21 2020-03-24 2236008 Ontario Inc. Automated driving system that merges heterogenous sensor data
US10908273B2 (en) * 2018-08-10 2021-02-02 GM Global Technology Operations LLC Object velocity vector estimation by multiple radars with different observation angles

Also Published As

Publication number Publication date
WO2022037370A1 (zh) 2022-02-24

Similar Documents

Publication Publication Date Title
CN109887057B (zh) 生成高精度地图的方法和装置
US10788830B2 (en) Systems and methods for determining a vehicle position
CN110148185B (zh) 确定成像设备坐标系转换参数的方法、装置和电子设备
AU2018282302B2 (en) Integrated sensor calibration in natural scenes
WO2020093378A1 (en) Vehicle positioning system using lidar
WO2021143286A1 (zh) 车辆定位的方法、装置、控制器、智能车和系统
US11747144B2 (en) Real time robust localization via visual inertial odometry
JP5992184B2 (ja) 画像データ処理装置、画像データ処理方法および画像データ処理用のプログラム
US10094919B2 (en) Radar-vision fusion for target velocity estimation
CN110411457B (zh) 基于行程感知与视觉融合的定位方法、系统、终端和存储介质
JP2019528501A (ja) マルチカメラシステムにおけるカメラ位置合わせ
CN112005079A (zh) 用于更新高清地图的系统和方法
CN110458885B (zh) 基于行程感知与视觉融合的定位系统和移动终端
US11866056B2 (en) Ballistic estimation of vehicle data
CN114730472A (zh) 车载相机的外部参数的标定方法及相关装置
CN115236714A (zh) 多源数据融合定位方法、装置、设备及计算机存储介质
Xian et al. Fusing stereo camera and low-cost inertial measurement unit for autonomous navigation in a tightly-coupled approach
WO2022037370A1 (zh) 一种运动估计方法及装置
WO2021056283A1 (en) Systems and methods for adjusting a vehicle pose
CN111207688B (zh) 在载运工具中测量目标对象距离的方法、装置和载运工具
CN113034538B (zh) 一种视觉惯导设备的位姿跟踪方法、装置及视觉惯导设备
CN112050830B (zh) 一种运动状态估计方法及装置
CN115508836A (zh) 一种运动状态确定方法及装置
CN113470342B (zh) 一种自运动估计的方法及装置
WO2022033139A1 (zh) 一种自运动估计方法及相关装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination