CN108627153A - A kind of rigid motion tracing system and its working method based on inertial sensor - Google Patents
A kind of rigid motion tracing system and its working method based on inertial sensor Download PDFInfo
- Publication number
- CN108627153A CN108627153A CN201810448583.2A CN201810448583A CN108627153A CN 108627153 A CN108627153 A CN 108627153A CN 201810448583 A CN201810448583 A CN 201810448583A CN 108627153 A CN108627153 A CN 108627153A
- Authority
- CN
- China
- Prior art keywords
- inertial sensor
- motion
- rigid body
- data
- sensor node
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/10—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
- G01C21/12—Navigation; 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/16—Navigation; 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/165—Navigation; 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
-
- G—PHYSICS
- G08—SIGNALLING
- G08C—TRANSMISSION SYSTEMS FOR MEASURED VALUES, CONTROL OR SIMILAR SIGNALS
- G08C17/00—Arrangements for transmitting signals characterised by the use of a wireless electrical link
- G08C17/02—Arrangements for transmitting signals characterised by the use of a wireless electrical link using a radio link
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Computer Networks & Wireless Communication (AREA)
- Navigation (AREA)
- User Interface Of Digital Computer (AREA)
Abstract
本发明公开了一种基于惯性传感器的刚体运动追踪系统及其工作方法,该系统包括可附着在刚体上的惯性传感器节点和数据采集处理终端。惯性传感器节点包含惯性传感器、计算模块和无线通信模块。数据采集处理终端包含无线通信模块和计算机系统。用户将惯性传感器节点布置在刚体上,节点实时采集传感器数据并无线传输至数据采集处理终端。当运动结束,计算机系统执行运动追踪和运动分析算法,对数据做预处理,获知惯性传感器节点的加速度,并计算其速度和运动轨迹,再基于惯性传感器节点的姿态以及刚体上其他点和传感器节点的相对位置关系对这些点的运动轨迹进行还原,最终可以呈现刚体选定点的运动轨迹并对运动过程中发生的平移和旋转进行分析。
The invention discloses a rigid body motion tracking system based on an inertial sensor and a working method thereof. The system includes an inertial sensor node and a data collection and processing terminal which can be attached to the rigid body. The inertial sensor node includes inertial sensors, computing modules and wireless communication modules. The data acquisition and processing terminal includes a wireless communication module and a computer system. The user arranges the inertial sensor node on the rigid body, and the node collects sensor data in real time and wirelessly transmits it to the data acquisition and processing terminal. When the motion is over, the computer system executes motion tracking and motion analysis algorithms, preprocesses the data, obtains the acceleration of the inertial sensor node, and calculates its velocity and motion trajectory, and then based on the attitude of the inertial sensor node and other points on the rigid body and sensor nodes The relative positional relationship of these points is restored to the trajectory of these points, and finally the trajectory of the selected point of the rigid body can be presented and the translation and rotation that occur during the movement can be analyzed.
Description
技术领域technical field
本发明涉及惯性传感器感知识别,运动追踪以及运动分析,具体是一种基于惯性传感器的刚体运动追踪系统及其工作方法。The invention relates to inertial sensor perception recognition, motion tracking and motion analysis, in particular to an inertial sensor-based rigid body motion tracking system and a working method thereof.
背景技术Background technique
运动追踪是指对目标对象的运动状态进行跟踪,主要包括测量、跟踪、记录目标对象在三维空间中的运动轨迹。当前运动追踪解决方案主要基于计算机图形学和图像处理技术,超声波、雷达、激光等定位技术以及惯性传感器。其中最主流的解决方案一般是基于计算机图形学和图像处理技术的。通过部署在三维空间中的多个摄像头,将运动物体(通常会布置标记节点)的运动状况以图像的方式记录下来,再由计算机系统对这些图像数据进行处理,最终得到运动物体在不同时刻的空间坐标(X,Y,Z)。基于计算机视觉的解决方案可以获得很高的精度,但其需要较大规模的部署,而且容易受视觉因素干扰,当运动物体受遮挡时精度会有很大偏差。Motion tracking refers to tracking the motion state of the target object, mainly including measuring, tracking, and recording the motion trajectory of the target object in three-dimensional space. Current motion tracking solutions are mainly based on computer graphics and image processing technologies, ultrasonic, radar, laser and other positioning technologies, and inertial sensors. The most mainstream solutions are generally based on computer graphics and image processing technology. Through multiple cameras deployed in the three-dimensional space, the motion status of the moving object (usually with marker nodes) is recorded in the form of images, and then the computer system processes these image data, and finally obtains the motion of the moving object at different times. Space coordinates (X,Y,Z). Solutions based on computer vision can achieve high accuracy, but they require large-scale deployment and are easily disturbed by visual factors. When moving objects are occluded, the accuracy will deviate greatly.
随着惯性传感器,可穿戴计算技术的发展,基于可穿戴惯性传感器以及智能终端的运动追踪技术也备受关注。借助可穿戴惯性传感器节点,用户可以实时获取惯性传感器数据,再由终端设备上的运动追踪和运动分析算法处理,可以实现对目标刚体的运动过程进行追踪和分析。目前一些基于惯性传感器的运动追踪工作主要关注于还原目标对象的运动姿态或者运动轨迹,缺少对运动过程中不同状态之间转换过程的分析,如何能够兼顾运动追踪与运动分析是一个值得探究的课题。With the development of inertial sensors and wearable computing technology, motion tracking technology based on wearable inertial sensors and smart terminals has also attracted much attention. With the help of wearable inertial sensor nodes, users can obtain inertial sensor data in real time, and then processed by the motion tracking and motion analysis algorithms on the terminal device, which can track and analyze the motion process of the target rigid body. At present, some motion tracking work based on inertial sensors mainly focuses on restoring the motion posture or motion trajectory of the target object, and lacks the analysis of the transition process between different states during the motion process. How to balance motion tracking and motion analysis is a topic worth exploring. .
发明内容Contents of the invention
为了解决上述技术问题,本发明提出一种基于惯性传感器的刚体运动追踪系统及工作方法,该系统和方法仅需要一个惯性传感器节点和数据采集处理终端,部署方便,操作简单,能够相对精准地捕捉附着惯性传感器节点的刚体的运动过程,获取各个时刻刚体的运动加速度和速度信息,还原刚体上多个选定点的运动轨迹,并为刚体运动过程建立数学模型,对运动过程中发生的平移和旋转进行可靠的分析。In order to solve the above technical problems, the present invention proposes a rigid body motion tracking system and working method based on inertial sensors. The system and method only need one inertial sensor node and a data acquisition and processing terminal, and are easy to deploy, easy to operate, and can relatively accurately capture The motion process of the rigid body attached to the inertial sensor node obtains the motion acceleration and velocity information of the rigid body at each moment, restores the motion trajectory of multiple selected points on the rigid body, and establishes a mathematical model for the motion process of the rigid body. Spin for reliable analysis.
为实现上述发明目的,本发明的技术方案如下:For realizing the above-mentioned purpose of the invention, the technical scheme of the present invention is as follows:
本发明的一种基于惯性传感器的刚体运动追踪系统,该系统包括:A kind of rigid body motion tracking system based on inertial sensor of the present invention, this system comprises:
惯性传感器节点:内置多种惯性传感器,具备一定的计算、存储和通信能力,可以方便地附着在刚体表面。Inertial sensor node: Built-in a variety of inertial sensors, with certain computing, storage and communication capabilities, can be easily attached to the surface of rigid bodies.
数据采集处理终端:具备计算、存储和通信能力的计算机系统,接收从惯性传感器节点传送过来的各种传感器数据,存储并处理,最终生成刚体运动的轨迹以及刚体运动过程中发生的平移和旋转的分析结果,并向用户展示。Data acquisition and processing terminal: a computer system with computing, storage and communication capabilities, receiving various sensor data transmitted from inertial sensor nodes, storing and processing, and finally generating the trajectory of rigid body motion and the translation and rotation during rigid body motion Analyze the results and present them to the user.
刚体是指在运动过程中以及受力作用后,形状和大小基本不变,且内部各点的相对位置不变或者变化程度小到可以忽略不计的物体。所述的刚体运动,是指刚体在三维空间中作旋转、平移的运动,可以理解为三维空间中一组三维坐标到另一组三维坐标的映射。此外,所述的运动追踪是指对目标对象的运动状态进行跟踪,主要包括测量、跟踪、记录目标对象在三维空间中的运动轨迹。A rigid body refers to an object whose shape and size are basically unchanged during the motion process and after the force is applied, and the relative position of each internal point is unchanged or the degree of change is negligible. The rigid body motion refers to the rotation and translation motion of a rigid body in three-dimensional space, which can be understood as the mapping from one set of three-dimensional coordinates to another set of three-dimensional coordinates in three-dimensional space. In addition, the motion tracking refers to tracking the motion state of the target object, which mainly includes measuring, tracking, and recording the motion trajectory of the target object in three-dimensional space.
所述惯性传感器节点固定安装于刚体表面,其内部包括:The inertial sensor node is fixedly installed on the surface of the rigid body, and its interior includes:
惯性传感器模块:包含加速度计,陀螺仪和磁力计,用于实时测量相应的惯性传感器数据。Inertial sensor module: contains accelerometer, gyroscope and magnetometer for real-time measurement of corresponding inertial sensor data.
计算模块:用于处理惯性传感器原始数据,通过内置的信号滤波算法实时获取惯性传感器节点的姿态信息,从而估计出重力加速度在惯性传感器自身参考坐标系三个轴上的分量,和原始加速度数据相减即可获取线性加速度。Calculation module: used to process the raw data of the inertial sensor, and obtain the attitude information of the inertial sensor node in real time through the built-in signal filtering algorithm, thereby estimating the components of the acceleration of gravity on the three axes of the inertial sensor's own reference coordinate system, which are compared with the original acceleration data Subtract to get the linear acceleration.
无线通信模块:用于与终端设备进行信息交互,主要用于接收数据采集处理终端的指令,并周期性地向终端设备发送传感器数据。Wireless communication module: used for information interaction with terminal equipment, mainly used to receive instructions from data acquisition and processing terminals, and periodically send sensor data to terminal equipment.
进一步的,所述的数据采集处理终端内部主要包括:Further, the inside of the data collection and processing terminal mainly includes:
无线通信模块:用于与惯性传感器节点进行信息交互。Wireless communication module: used for information interaction with inertial sensor nodes.
计算机系统:是指具备计算、存储和通信能力的计算机系统,通过无线通信模块接收惯性传感器数据作为输入,执行运动追踪和运动分析算法,获知刚体每一时刻的加速度和速度信息,计算刚体运动的轨迹,并分析刚体在运动过程中平移和旋转的情况。Computer system: refers to a computer system with computing, storage and communication capabilities. It receives inertial sensor data as input through a wireless communication module, executes motion tracking and motion analysis algorithms, obtains the acceleration and velocity information of the rigid body at each moment, and calculates the rigid body motion. trajectory, and analyze the translation and rotation of the rigid body during motion.
基于上述追踪系统,本发明还介绍工作方法,包括以下步骤:Based on the above-mentioned tracking system, the present invention also introduces a working method, comprising the following steps:
1)用户将单个惯性传感器节点布置在刚体上,准备触发刚体运动。1) The user arranges a single inertial sensor node on the rigid body, ready to trigger the rigid body motion.
2)用户为惯性传感器节点与数据采集处理终端建立无线连接。2) The user establishes a wireless connection between the inertial sensor node and the data acquisition and processing terminal.
3)用户在数据采集处理终端开启数据采集,惯性传感器节点开始实时采集各类传感器数据并周期性发送至数据采集处理终端,刚体运动开始。3) The user starts data collection at the data collection and processing terminal, and the inertial sensor node starts to collect various sensor data in real time and periodically sends them to the data collection and processing terminal, and the rigid body motion starts.
4)刚体运动结束,用户在数据采集处理终端终止数据采集。4) After the rigid body motion ends, the user terminates the data collection at the data collection and processing terminal.
5)数据采集处理终端通过计算机系统执行运动追踪和运动分析算法,获知刚体每一时刻的加速度和速度信息,生成刚体运动的轨迹以及刚体运动过程中发生的平移和旋转的分析结果,并向用户展示。5) The data acquisition and processing terminal executes the motion tracking and motion analysis algorithm through the computer system, obtains the acceleration and velocity information of the rigid body at each moment, generates the trajectory of the rigid body motion and the analysis results of the translation and rotation that occur during the motion of the rigid body, and reports to the user exhibit.
进一步的,其中,所述步骤5)中的运动追踪和运动分析算法,包括以下步骤:Further, wherein, the motion tracking and motion analysis algorithm in the step 5) includes the following steps:
步骤5.1数据预处理,先对原始传感器数据进行平滑操作,然后计算每一时刻的三轴线性加速度数据和陀螺仪数据的幅值,利用阈值判断的方法确定运动起始和终止的时刻。具体来讲,运动起始时刻,两类传感器数据的幅值会分别超过一个阈值,运动终止的时刻,两类传感器数据的幅值都会回到相应阈值之下。通过大量实验数据选定合理的阈值,即可准确并抽取出运动过程对应的数据。接着对抽取出的线性加速度数据进行参考坐标系的转换,这里的参考坐标系的转换是指利用重力加速度向量和磁力计数据向量构建一个地理坐标系,具体来讲,首先由于指向北方,所以水平面上指向地理东方的向量可以由下式求解:Step 5.1 Data preprocessing, first smoothing the original sensor data, then calculating the three-axis linear acceleration data and the amplitude of the gyroscope data at each moment, and using the threshold judgment method to determine the start and end moments of the movement. Specifically, at the beginning of the movement, the amplitudes of the two types of sensor data will respectively exceed a threshold, and at the moment of the end of the movement, the amplitudes of the two types of sensor data will return below the corresponding threshold. By selecting a reasonable threshold through a large amount of experimental data, the data corresponding to the motion process can be accurately extracted. Then, the extracted linear acceleration data is converted to the reference coordinate system. The conversion of the reference coordinate system here refers to the use of the gravitational acceleration vector and the magnetometer data vector To construct a geographic coordinate system, specifically, firstly due to points north, so a vector pointing geographically east on the horizon It can be solved by the following formula:
又由于不严格水平,所以按如下算式得到水平面上指向地地理北方的向量 And because of It is not strictly horizontal, so the vector pointing to the geographical north on the horizontal plane can be obtained according to the following formula
于是这三个正交的向量便可以构成地理坐标系,这三个向量各自规范化之后可以构建方向余弦矩阵DCM,这样便可以实现从惯性传感器节点自身坐标系到地理坐标系的转换:then These three orthogonal vectors can form a geographic coordinate system. After the three vectors are normalized, the direction cosine matrix DCM can be constructed, so that the transformation from the inertial sensor node's own coordinate system to the geographic coordinate system can be realized:
(Xg,Yg,Zg)’=DCM×(X,Y,Z)’(X g ,Y g ,Z g )'=DCM×(X,Y,Z)'
其中(Xg,Yg,Zg)为地理坐标系下的传感器数据,(X,Y,Z)为基于惯性传感器节点自身坐标系的原始数据。Where (X g , Y g , Z g ) is the sensor data in the geographic coordinate system, and (X, Y, Z) is the original data based on the inertial sensor node's own coordinate system.
根据具体应用场景的需要,可以利用同样的方法对地理坐标系再做一次变换,得到应用场景对应的参考坐标系,最后我们就可以得到运动过程中相对于固定参考坐标系的线性加速度数据。According to the needs of the specific application scenario, the geographic coordinate system can be transformed again using the same method to obtain the reference coordinate system corresponding to the application scenario. Finally, we can obtain the linear acceleration data relative to the fixed reference coordinate system during the motion.
步骤5.2执行运动追踪,基于惯性传感器节点上一时刻的运动状态和当前时刻的惯性传感器数据更新当前时刻惯性传感器节点的运动状态,从而获取惯性传感器节点每一时刻在三维空间中的相对位置和姿态信息。Step 5.2 Execute motion tracking, update the motion state of the inertial sensor node at the current moment based on the motion state of the inertial sensor node at the previous moment and the inertial sensor data at the current moment, so as to obtain the relative position and attitude of the inertial sensor node in three-dimensional space at each moment information.
步骤5.3执行轨迹还原,根据惯性传感器节点当前时刻在空间中的相对位置以及其姿态,结合刚体上选定点和惯性传感器节点的相对位置关系,对刚体上这些选定点的位置进行估计,从而还原出刚体上多个点的运动轨迹。Step 5.3 Execute trajectory restoration, according to the relative position of the inertial sensor node in space at the current moment and its attitude, combined with the relative position relationship between the selected point on the rigid body and the inertial sensor node, estimate the position of these selected points on the rigid body, so that Restore the trajectory of multiple points on the rigid body.
步骤5.4执行运动分析,选取刚体中的某个平面的运动轨迹,对选定时间片长度的运动过程进行分析,建立数学模型,将连续的运动等效分解为平移和旋转的组合,运动的旋转情况可以利用四元数的方法求解,最终输出运动过程中发生的平移和旋转的情况,平移由三维空间的一个向量来表示,旋转可以由一个四元数(包含旋转轴和旋转角)来表示。Step 5.4 performs motion analysis, selects the motion trajectory of a certain plane in the rigid body, analyzes the motion process of the selected time slice length, establishes a mathematical model, and decomposes the continuous motion equivalently into a combination of translation and rotation, and the rotation of motion The situation can be solved by using the quaternion method, and finally output the translation and rotation that occurs during the motion process. The translation is represented by a vector in three-dimensional space, and the rotation can be represented by a quaternion (including the rotation axis and rotation angle). .
进一步的,所述步骤5.2中惯性传感器节点的运动状态主要包括惯性传感器节点在空间中的速度、相对位置及其姿态,具体的运动状态更新步骤包括:Further, the motion state of the inertial sensor node in the step 5.2 mainly includes the speed, relative position and posture of the inertial sensor node in space, and the specific motion state update steps include:
步骤5.2.1相对位置信息更新:运动追踪算法以惯性传感器节点开始运动的位置作为应用场景对应的参考坐标系的原点,在时域上先对线性加速度进行一次积分,得到各个时刻的速度信息,再利用零速校正等误差修正手段,对速度进行补偿。接着对补偿后的速度在时域上继续做一次积分,求解相邻时刻发生的位移,实现在每一采样时刻对惯性传感器节点在空间中的相对位置(即三维坐标)进行更新。Step 5.2.1 Relative position information update: The motion tracking algorithm takes the position where the inertial sensor node starts to move as the origin of the reference coordinate system corresponding to the application scene, and first integrates the linear acceleration in the time domain to obtain the speed information at each moment. Then use error correction methods such as zero speed correction to compensate the speed. Then continue to integrate the compensated velocity in the time domain to solve the displacement occurring at adjacent moments, so as to update the relative position of the inertial sensor node in space (that is, the three-dimensional coordinates) at each sampling moment.
步骤5.2.2姿态信息更新:运动初始阶段由磁力计和重力加速度的瞬时数据构建包含惯性传感器节点姿态信息的矩阵作为惯性传感器节点的初始姿态,并记录下初始阶段节点局部坐标系至地理坐标系的旋转矩阵,后续阶段利用陀螺仪的数据在时域上的一次积分对初始姿态进行更新,并在一些合适的时间点利用磁力计和重力加速度测得的姿态进行校准,从而在每一采样时刻对惯性传感器节点的姿态信息进行更新。Step 5.2.2 Attitude information update: In the initial stage of the movement, a matrix containing the attitude information of the inertial sensor nodes is constructed from the instantaneous data of the magnetometer and the acceleration of gravity as the initial attitude of the inertial sensor nodes, and the local coordinate system of the initial stage is recorded to the geographic coordinate system In the subsequent stage, the initial attitude is updated by an integration of the gyroscope data in the time domain, and the attitude measured by the magnetometer and the acceleration of gravity is calibrated at some suitable time points, so that at each sampling moment Update the attitude information of the inertial sensor node.
本发明提供了一种基于惯性传感器的刚体运动追踪系统及其工作方法,相对于基于计算机视觉的运动追踪系统,更加轻量级,易部署,成本低廉且不受视觉因素限制。相较于其他同类型的基于惯性传感器的运动追踪系统,本系统能够还原刚体上多个点的运动轨迹,并且给出可靠的运动分析结果。The invention provides a rigid body motion tracking system based on an inertial sensor and a working method thereof. Compared with a motion tracking system based on computer vision, it is lighter in weight, easy to deploy, low in cost and not limited by visual factors. Compared with other motion tracking systems based on inertial sensors of the same type, this system can restore the motion trajectories of multiple points on the rigid body and give reliable motion analysis results.
本发明的有益效果具体表现为:The beneficial effects of the present invention are embodied as:
1、提出了一种新的对刚体在三维空间中运动过程进行分析的方法:在分析运动轨迹之外引入对运动过程中平移和旋转情况的分析。通过建立几何模型来将刚体在三维空间中的运动分解成平移和旋转的组合,使得运动分析的结果更为直观。1. A new method for analyzing the motion process of a rigid body in three-dimensional space is proposed: in addition to analyzing the motion trajectory, the analysis of the translation and rotation during the motion process is introduced. By establishing a geometric model to decompose the motion of a rigid body in three-dimensional space into a combination of translation and rotation, the results of motion analysis are more intuitive.
2、可以实现对刚体上多个选定点的运动轨迹进行还原:通过构建一个与应用场景对应的固定的参考坐标系,充分利用惯性传感器实时姿态信息以及刚体上选定点与惯性传感器节点的相对位置关系,实现对刚体上多个选定点的运动轨迹进行还原。2. It can realize the restoration of the motion trajectory of multiple selected points on the rigid body: by constructing a fixed reference coordinate system corresponding to the application scene, make full use of the real-time attitude information of the inertial sensor and the relationship between the selected point on the rigid body and the inertial sensor node The relative position relationship realizes the restoration of the motion trajectory of multiple selected points on the rigid body.
3、容易部署:不需要在目标周围提前部署大量摄像头以及线缆,只需要用户将惯性传感器节点布置于待追踪的刚体之上即可。3. Easy to deploy: It is not necessary to deploy a large number of cameras and cables around the target in advance, and only needs the user to arrange the inertial sensor nodes on the rigid body to be tracked.
4、成本相对低廉:不需要购置高成本的摄像头、交换机、线缆等设备,只需一台常见的数据采集处理终端和成本相对低廉的惯性传感器节点即可实现对刚体运动过程的追踪和分析。4. The cost is relatively low: there is no need to purchase high-cost cameras, switches, cables and other equipment, only a common data acquisition and processing terminal and a relatively low-cost inertial sensor node can realize the tracking and analysis of the rigid body motion process .
附图说明Description of drawings
图1为基于惯性传感器的刚体运动追踪系统架构图。Figure 1 is an architecture diagram of a rigid body motion tracking system based on inertial sensors.
图2为基于惯性传感器的刚体运动追踪系统的运动追踪和运动分析算法流程图。Fig. 2 is a flow chart of the motion tracking and motion analysis algorithm of the inertial sensor-based rigid body motion tracking system.
图3为基于惯性传感器的刚体运动追踪系统系统工作方法涉及的参考坐标系示意图。FIG. 3 is a schematic diagram of a reference coordinate system involved in the working method of the inertial sensor-based rigid body motion tracking system.
图4为刚体上选定点与惯性传感器节点的相对位置关系图。Figure 4 is a diagram of the relative position relationship between the selected point on the rigid body and the inertial sensor node.
图5为刚体运动示意图。Figure 5 is a schematic diagram of rigid body motion.
图6为刚体运动过程建模分析示意图。Figure 6 is a schematic diagram of the modeling and analysis of the rigid body motion process.
具体实施方式Detailed ways
为了便于本领域技术人员的理解,下面结合附图对本发明作进一步的说明。In order to facilitate the understanding of those skilled in the art, the present invention will be further described below in conjunction with the accompanying drawings.
图1是基于惯性传感器的刚体运动追踪系统架构图,主要由附着在刚体上的惯性传感器节点和数据采集处理终端两部分组成,如图中所示,两者主要通过蓝牙无线通信进行数据交互。Figure 1 is an architecture diagram of a rigid body motion tracking system based on inertial sensors. It is mainly composed of two parts: an inertial sensor node attached to a rigid body and a data acquisition and processing terminal. As shown in the figure, the two mainly exchange data through Bluetooth wireless communication.
所述惯性传感器节点固定安装于刚体表面,可以选择粘贴于刚体表面,确保惯性传感器模块与刚体的相对位置关系固定。惯性传感器节点中包含惯性传感器、计算模块和通信模块。惯性传感器包含加速度计,陀螺仪和磁力计,用于实时测量相应的惯性传感器数据;计算模块是用于处理原始传感器数据的,通过内置的信号滤波算法实时获取惯性传感器节点的姿态信息,从而估计出重力加速度在惯性传感器自身参考坐标系三个轴上的分量,和原始加速度数据相减即可获取线性加速度;无线通信模块用于与终端设备进行信息交互,主要是接收数据采集处理终端的指令,并周期性地向终端设备发送传感器数据。The inertial sensor node is fixedly installed on the surface of the rigid body, and can be optionally pasted on the surface of the rigid body to ensure that the relative positional relationship between the inertial sensor module and the rigid body is fixed. The inertial sensor node contains inertial sensors, computing modules and communication modules. The inertial sensor includes accelerometer, gyroscope and magnetometer, which are used to measure the corresponding inertial sensor data in real time; the calculation module is used to process the raw sensor data, and obtain the attitude information of the inertial sensor node in real time through the built-in signal filtering algorithm, so as to estimate The linear acceleration can be obtained by subtracting the components of the acceleration of gravity on the three axes of the inertial sensor's own reference coordinate system from the original acceleration data; the wireless communication module is used for information interaction with the terminal device, mainly to receive instructions from the data acquisition and processing terminal , and periodically send sensor data to the terminal device.
数据采集处理终端主要由无线通信模块和计算机系统组成。无线通信模块用于与惯性传感器节点进行信息交互;计算机系统是指具备计算、存储和通信能力的计算机系统,它通过无线通信模块接收惯性传感器数据作为输入,执行运动追踪和运动分析算法,计算刚体运动的轨迹,并分析刚体在运动过程中平移和旋转的情况。The data acquisition and processing terminal is mainly composed of a wireless communication module and a computer system. The wireless communication module is used for information interaction with the inertial sensor nodes; the computer system refers to a computer system with computing, storage and communication capabilities, which receives inertial sensor data as input through the wireless communication module, executes motion tracking and motion analysis algorithms, and calculates rigid body The trajectory of the movement, and analyze the translation and rotation of the rigid body during the movement.
图2是运动追踪和运动分析算法流程图,运动追踪和运动分析算法主要包括数据预处理,运动追踪,轨迹还原和运动分析。Figure 2 is a flow chart of the motion tracking and motion analysis algorithm. The motion tracking and motion analysis algorithm mainly includes data preprocessing, motion tracking, trajectory restoration, and motion analysis.
数据预处理部分主要是按时间顺序从前往后扫描线性加速度和陀螺仪传感器数据,准确寻找出运动开始和运动结束的时刻。首先对传感器数据进行平滑,可以用简单的均值平滑滤除一些高频噪声。然后分别计算三轴线性加速度和陀螺仪数据的幅值AccM和GyroM,即分别求解两个三维向量的2-范数。这里,通过前期大量实验数据的分析,我们分别对两种三轴传感器数据的幅值有一个阈值thresA和thresG,我们通过设置一个定长的滑动窗口,统计窗口内传感器数据大于和小于对应阈值的比率。当两种传感器数据的幅值小于设定阈值的比率均大于ρ(本系统设置为80%)时,则说明当前窗口运动状态不活跃,而当有一个传感器数据的幅值大于设定阈值的比率大于ρ时,则说明当前窗口运动状态活跃。我们找出第一个不活跃状态迁移到活跃状态的窗口,然后设定运动起始时间ts,再往后搜索由活跃状态向不活跃状态迁移的窗口,设定运动终止时间te,至此,我们只需截取ts和te之间的传感器数据,再对截取后的线性加速度数据进行参考坐标系的转换,作为算法往后部分的输入。这里的参考坐标系的转换是指利用重力加速度向量和磁力计数据向量构建一个地理坐标系,具体来讲,首先由于指向北方,所以水平面上指向地理东方的向量可以由下式求解:The data preprocessing part is mainly to scan the linear acceleration and gyroscope sensor data from front to back in time order, and accurately find out the moment when the movement starts and ends. First smooth the sensor data, some high frequency noise can be filtered out with simple mean smoothing. Then calculate the amplitudes AccM and GyroM of the three-axis linear acceleration and gyroscope data respectively, that is, solve the 2-norm of the two three-dimensional vectors respectively. Here, through the analysis of a large number of experimental data in the early stage, we have thresholds thresA and thresG for the amplitudes of the two three-axis sensor data respectively. We set a fixed-length sliding window to count the sensor data in the window greater than and less than the corresponding threshold. ratio. When the ratio of the amplitudes of the two sensor data is less than the set threshold is greater than ρ (this system is set to 80%), it means that the current window motion state is not active, and when the amplitude of one sensor data is greater than the set threshold When the ratio is greater than ρ, it means that the current window motion state is active. We find out the first window where the inactive state migrates to the active state, and then set the motion start time t s , and then search for the window that transitions from the active state to the inactive state, and set the motion end time t e , so far , we only need to intercept the sensor data between t s and t e , and then convert the intercepted linear acceleration data to the reference coordinate system as the input of the later part of the algorithm. The conversion of the reference coordinate system here refers to the use of the gravitational acceleration vector and the magnetometer data vector To construct a geographic coordinate system, specifically, firstly due to points north, so a vector pointing geographically east on the horizon It can be solved by the following formula:
又由于不严格水平,所以按如下算式得到水平面上指向地地理北方的向量 And because of It is not strictly horizontal, so the vector pointing to the geographical north on the horizontal plane can be obtained according to the following formula
于是这三个正交的向量便可以构成地理坐标系,这三个向量各自规范化之后可以构建方向余弦矩阵DCM,这样便可以实现从惯性传感器节点自身坐标系到地理坐标系的转换:then These three orthogonal vectors can form a geographic coordinate system. After the three vectors are normalized, the direction cosine matrix DCM can be constructed, so that the transformation from the inertial sensor node's own coordinate system to the geographic coordinate system can be realized:
(Xg,Yg,Zg)’=DCM×(X,Y,Z)’(X g ,Y g ,Z g )'=DCM×(X,Y,Z)'
其中(Xg,Yg,Zg)为地理坐标系下的传感器数据,(X,Y,Z)为基于惯性传感器节点自身坐标系的原始数据。Where (X g , Y g , Z g ) is the sensor data in the geographic coordinate system, and (X, Y, Z) is the original data based on the inertial sensor node's own coordinate system.
根据具体应用场景的需要,可以利用同样的方法对地理坐标系再做一次变换,得到应用场景对应的参考坐标系,最后我们就可以得到运动过程中相对于固定参考坐标系的线性加速度数据。According to the needs of the specific application scenario, the geographic coordinate system can be transformed again using the same method to obtain the reference coordinate system corresponding to the application scenario. Finally, we can obtain the linear acceleration data relative to the fixed reference coordinate system during the motion.
运动追踪部分主要是利用预处理阶段生成的基于应用场景对应的参考坐标系下的线性加速度数据和陀螺仪数据,对惯性传感器节点的运动轨迹进行更新,同时生成每一采样时刻的惯性传感器节点姿态信息。具体来说,轨迹更新部分,算法以惯性传感器节点开始运动的位置作为应用场景对应的参考坐标系的原点,在时域上先对线性加速度进行一次积分,得到各个时刻的速度信息,由于线性加速度数据存在噪声,积分会导致误差累积,我们利用零速校正等误差修正手段,对速度进行修正。接着对补偿后的速度在时域上继续做一次积分,求解相邻时刻发生的位移,实现对惯性传感器节点在空间中的相对位置(即三维坐标)进行更新。所述的零速校正主要是基于运动终止时刻,惯性传感器节点的速度为零的这样一个假设,这样在对线性加速度做第一次积分得到各时刻的速度的时候可以对速度做一个补偿,从而可以减少累积误差。姿态更新部分,运动初始阶段由磁力计和重力加速度的瞬时数据构建包含惯性传感器节点姿态信息的矩阵作为惯性传感器节点的初始姿态,并记录下初始阶段节点局部坐标系至地理坐标系的旋转矩阵,后续阶段利用陀螺仪的数据在时域上的一次积分对初始姿态进行更新,并在一些合适的时间点利用磁力计和重力加速度测得的姿态进行校准,以消除陀螺仪带来的累积误差,从而对惯性传感器节点的姿态信息进行更新。这里所说的合适的时间点主要由这些因素判定:The motion tracking part mainly uses the linear acceleration data and gyroscope data generated in the preprocessing stage based on the reference coordinate system corresponding to the application scene to update the motion trajectory of the inertial sensor node, and at the same time generate the attitude of the inertial sensor node at each sampling moment information. Specifically, in the trajectory update part, the algorithm uses the position where the inertial sensor node starts to move as the origin of the reference coordinate system corresponding to the application scene, and first integrates the linear acceleration in the time domain to obtain the speed information at each moment. There is noise in the data, and integration will lead to error accumulation. We use error correction methods such as zero-speed correction to correct the speed. Then continue to integrate the compensated velocity in the time domain to solve the displacement occurring at adjacent moments, so as to update the relative position of the inertial sensor node in space (that is, the three-dimensional coordinates). The zero-speed correction is mainly based on the assumption that the speed of the inertial sensor node is zero at the end of the motion, so that the speed can be compensated when the linear acceleration is integrated for the first time to obtain the speed at each moment, so that Cumulative error can be reduced. In the attitude update part, in the initial stage of the movement, a matrix containing the attitude information of the inertial sensor node is constructed from the instantaneous data of the magnetometer and the acceleration of gravity as the initial attitude of the inertial sensor node, and the rotation matrix from the local coordinate system of the node to the geographic coordinate system in the initial stage is recorded. In the subsequent stage, an integration of gyroscope data in the time domain is used to update the initial attitude, and at some suitable time points, the attitude measured by the magnetometer and the acceleration of gravity is used for calibration to eliminate the cumulative error caused by the gyroscope. In this way, the attitude information of the inertial sensor node is updated. The appropriate time point mentioned here is mainly determined by these factors:
1)实时的角速度信息,当角速度的值在一定阈值以下,我们认为此刻由磁力计和重力加速度测得的姿态信息是比较可靠的,可以用于校准。1) Real-time angular velocity information. When the value of the angular velocity is below a certain threshold, we believe that the attitude information measured by the magnetometer and the acceleration of gravity at this moment is relatively reliable and can be used for calibration.
2)陀螺仪计算得到的姿态变化与利用磁力计和重力加速度得到的姿态变化的相似性,当由两种手段测得的姿态变化比较相似时,我们认为此刻由磁力计和重力加速度测得的姿态信息也是比较可靠的,可以用于校准。2) The attitude change calculated by the gyroscope is similar to the attitude change obtained by using the magnetometer and the acceleration of gravity. When the attitude changes measured by the two methods are similar, we believe that the attitude change measured by the magnetometer and the acceleration of gravity Attitude information is also relatively reliable and can be used for calibration.
轨迹还原部分利用运动追踪部分生生成的惯性传感器节点轨迹信息和姿态信息,结合预先获知的节点自身参考坐标系下刚体上选定点与惯性传感器节点的相对位置关系,得到应用场景对应的固定参考坐标系下刚体上选定点与惯性传感器节点的相对位置关系,然后对刚体上选定点的运动轨迹进行估算还原。所述刚体上选定点与惯性传感器节点的相对位置关系具体可以参见下文对图4的说明。The trajectory restoration part uses the inertial sensor node trajectory information and attitude information generated by the motion tracking part, combined with the pre-known relative position relationship between the selected point on the rigid body and the inertial sensor node in the node's own reference coordinate system, to obtain the fixed reference corresponding to the application scene. The relative position relationship between the selected point on the rigid body and the inertial sensor node in the coordinate system, and then estimate and restore the trajectory of the selected point on the rigid body. For the relative position relationship between the selected point on the rigid body and the inertial sensor node, please refer to the description of FIG. 4 below for details.
运动分析部分以轨迹还原部分输出的刚体上选定点的轨迹作为输入,利用我们建立的刚体运动过程的数学模型进行分析。具体来讲,我们选定刚体上的一个平面上不在同一直线上的三个点,他们的运动轨迹是已知的,相当于这个平面的运动过程是已知的,这样利用我们建立的刚体运动过程数学模型可以将刚体在一个时刻至另一时刻的运动过程等效简化为一次平移和一次旋转的组合。最后我们给出运动的分析结果便是运动过程中所发生平移的向量旋转过程的旋转轴以及绕旋转轴逆时针转过的角度θequal。关于刚体运动过程的数学模型可以参见下文对图5,图6的说明。The motion analysis part takes the trajectory of the selected point on the rigid body output by the trajectory restoration part as input, and uses the mathematical model of the rigid body motion process we have established for analysis. Specifically, we select three points on a plane on the rigid body that are not on the same straight line, and their motion trajectories are known, which means that the motion process of this plane is known, so that using the rigid body motion we established The process mathematical model can equivalently simplify the motion process of a rigid body from one moment to another into a combination of a translation and a rotation. Finally, we give the analysis result of the motion is the vector of the translation during the motion Axis of rotation for the rotation process And the angle θ equal rotated counterclockwise around the axis of rotation. For the mathematical model of the rigid body motion process, please refer to the description of Fig. 5 and Fig. 6 below.
图3是基于惯性传感器的刚体运动追踪系统工作方法中涉及的几种参考坐标系的示意图。图中左半部分以点N为坐标系原点、XNYNZN为坐标轴的坐标系是惯性传感器节点自身的坐标系,惯性传感器节点产生的原始传感器数据都是基于这个坐标系的,而以点O为坐标系原点,XYZ为坐标轴的坐标系则是前面所述的应用场景对应的参考坐标系,这里只是以其中一种可能的情况作为示例,具体应用场景对应的参考坐标系可以有不同的定义。图中右半部所示的是地理坐标系,我们可以看到X轴指向地理东方,Y轴指向地理北方,Z轴指向重力的反方向。在实际的算法流程中,参考坐标系转换的过程便是将线性加速度数据由惯性传感器节点自身的坐标系转换至地理坐标系,再由地理坐标系转换至应用场景对应的参考坐标系。FIG. 3 is a schematic diagram of several reference coordinate systems involved in the working method of the inertial sensor-based rigid body motion tracking system. In the left half of the figure, the coordinate system with point N as the origin of the coordinate system and X N Y N Z N as the coordinate axis is the coordinate system of the inertial sensor node itself. The original sensor data generated by the inertial sensor node is based on this coordinate system. The coordinate system with point O as the origin and XYZ as the coordinate axis is the reference coordinate system corresponding to the above-mentioned application scenario. Here is just one of the possible situations as an example. The reference coordinate system corresponding to the specific application scenario There can be different definitions. The right half of the figure shows the geographic coordinate system. We can see that the X-axis points to the geographic east, the Y-axis points to the geographic north, and the Z-axis points to the opposite direction of gravity. In the actual algorithm process, the process of converting the reference coordinate system is to convert the linear acceleration data from the coordinate system of the inertial sensor node itself to the geographic coordinate system, and then convert from the geographic coordinate system to the reference coordinate system corresponding to the application scene.
图4是刚体上一些选定点与惯性传感器节点的相对位置关系的示意。图中点N代表惯性传感器节点,在刚体上选取不在同一平面的O、A、B、C四个点作为示例,图中给出的坐标系是惯性传感器节点自身的坐标系。简单起见,这里将点N作为坐标系原点。在这个坐标系下,由于刚体本身的形状大小以及惯性传感器的部署位置是已知的,所以点O、A、B、C、 N的坐标都是可以测算得到的,亦即图中标示出的四个向量都是可以测算得到的。由此我们便可以通过惯性传感器节点N在三维空间中的相对位置推算刚体上选定点 O、A、B、C在三维空间中的相对位置。具体来说,就是每一时刻都将选定点与惯性传感器节点的相对位置关系(即四个向量)由节点自身坐标系转换到地理坐标系,再转换至图3所示的应用场景对应的参考坐标系。而我们在运动追踪部分求解的惯性传感器节点 N的运动轨迹(即每一时刻惯性传感器节点在三维空间中的相对位置)是基于应用场景对应的参考坐标系的,这样就可以利用应用场景对应的参考坐标系下的四个向量求解点O、A、B、C在应用场景对应的参考坐标系下的坐标了,最终可以实现对刚体上选定点的轨迹还原。Figure 4 is a schematic diagram of the relative positional relationship between some selected points on the rigid body and the inertial sensor nodes. Point N in the figure represents the inertial sensor node. Four points O, A, B, and C that are not on the same plane are selected on the rigid body as an example. The coordinate system given in the figure is the coordinate system of the inertial sensor node itself. For simplicity, point N is taken as the origin of the coordinate system here. In this coordinate system, since the shape and size of the rigid body itself and the deployment position of the inertial sensor are known, the coordinates of points O, A, B, C, and N can be measured and obtained, that is, the points marked in the figure four vectors are all measurable. Therefore, we can estimate the relative positions of the selected points O, A, B, and C on the rigid body in the three-dimensional space through the relative position of the inertial sensor node N in the three-dimensional space. Specifically, the relative position relationship between the selected point and the inertial sensor node (that is, the four vectors ) is converted from the node's own coordinate system to the geographic coordinate system, and then converted to the reference coordinate system corresponding to the application scenario shown in FIG. 3 . The trajectory of the inertial sensor node N (that is, the relative position of the inertial sensor node in the three-dimensional space at each moment) that we solve in the motion tracking part is based on the reference coordinate system corresponding to the application scene, so that the application scene can be used. Four vectors in the reference coordinate system Solve the coordinates of points O, A, B, and C in the reference coordinate system corresponding to the application scene, and finally realize the restoration of the trajectory of the selected point on the rigid body.
图5是刚体运动过程的一个示例,不失一般性,这里用刚体上不在同一平面的四个点来代表一个刚体。我们可以看到图中刚体四个顶点的起始位置为O、A、B、C,经过运动之后,对应顶点的位置变换到了O′、A′、B′、C′,这就是刚体运动前后状态的一般化表示,如前文所述,相当于发生了一组三维坐标(O、A、B、C)到另一组三维坐标(O′、A′、B′、C′)的映射。Figure 5 is an example of the motion process of a rigid body. Without loss of generality, four points on the rigid body that are not on the same plane are used to represent a rigid body. We can see that the starting positions of the four vertices of the rigid body in the figure are O, A, B, and C. After the movement, the positions of the corresponding vertices are transformed to O', A', B', and C'. This is the movement of the rigid body before and after The generalized representation of state, as mentioned above, is equivalent to the mapping from one set of three-dimensional coordinates (O, A, B, C) to another set of three-dimensional coordinates (O', A', B', C').
图6是对刚体运动过程分析建模的一个示意图,刚体的运动过程是比较复杂的,但我们在分析过程中可以在保证运动前后状态不变的情况下对运动过程进行等效的变换以达到简化运动过程的目的。我们将相邻时间片之间刚体的运动过程分解为平移和旋转,这需要建立几何模型来描述。图6(a)是对运动过程进行分解的第一步,我们选取刚体上的一个点(O),其指向运动之后对应点(O′)的向量即是刚体作平移运动的向量。如图6(b),平移过后,刚体位置变换为O1、A1、B1、C1,其中O1与O′重合,在第二步中,我们需要将刚体的边O1A1绕过点O′的轴旋转至与O′A′重合,旋转轴的计算如下所示:Figure 6 is a schematic diagram of the analysis and modeling of the rigid body motion process. The motion process of the rigid body is relatively complicated, but we can perform equivalent transformations on the motion process to achieve The purpose of simplifying the exercise process. We decompose the motion process of a rigid body between adjacent time slices into translation and rotation, which requires a geometric model to describe. Figure 6(a) is the first step in decomposing the motion process. We select a point (O) on the rigid body, which points to the vector corresponding to the point (O′) after the motion That is, the vector for the translational motion of the rigid body. As shown in Figure 6(b), after the translation, the position of the rigid body is transformed into O 1 , A 1 , B 1 , and C 1 , where O 1 coincides with O′. In the second step, we need to convert the side O 1 A 1 of the rigid body The axis around the point O' is rotated to coincide with O'A', the axis of rotation The calculation of is as follows:
而旋转角度θ0则可以由向量与向量点乘获取。如图6(c),第一次旋转过后,刚体位置变换为O2、A2、B2、C2,边O2A2与边O′A′重合。此时,我们选择与O2、A2不在同一直线上的第三个点B2,在这一步中,刚体需要绕轴(转换为单位向量)旋转至平面O2A2B2与平面O′A′B′共面,旋转角θ1可以直接计算两个平面的夹角获取。具体地,这里是通过先求解图6(c) 中所示的两个平面的法向量和再求解两个法向量之间的夹角θ来获取两个平面之间的夹角θ1。最后,图(d)展示了经过第二次旋转之后的刚体位置O3、A3、B3、C3,亦即位置O′、A′、 B′、C′。图4展示的一次平移两次旋转的刚体运动分解过程是比较直观的,实际计算过程中,平移的变换过程利用平移过程向量的加减便能实现,旋转的变换过程则可以利用四元数的方法来求解。另外,我们可以利用四元数的方法将两次旋转进行复合,这样,刚体的运动过程就简化为一次平移和一次绕固定轴的旋转。我们以四元数和代表上述过程中的两次旋转,为旋转前点坐标的四元数表示,为第一次旋转后点坐标的四元数表示,第二次旋转后点坐标的四元数表示,具体的旋转复合过程如下所示:And the rotation angle θ 0 can be determined by the vector with vector Click to get. As shown in Figure 6(c), after the first rotation, the position of the rigid body is transformed into O 2 , A 2 , B 2 , and C 2 , and the side O 2 A 2 coincides with the side O'A'. At this point, we choose the third point B 2 which is not on the same straight line as O 2 and A 2 , in this step, the rigid body needs to go around the axis (converted to a unit vector ) is rotated until the plane O 2 A 2 B 2 is coplanar with the plane O'A'B', and the rotation angle θ 1 can be obtained by directly calculating the angle between the two planes. Specifically, here is by first solving the normal vectors of the two planes shown in Figure 6(c) and Then solve the angle θ between the two normal vectors to obtain the angle θ 1 between the two planes. Finally, figure (d) shows the rigid body positions O 3 , A 3 , B 3 , C 3 after the second rotation, ie positions O', A', B', C'. The rigid body motion decomposition process of one translation and two rotations shown in Figure 4 is relatively intuitive. In the actual calculation process, the transformation process of translation can be realized by adding and subtracting the translation process vector, and the transformation process of rotation can be realized by using the quaternion method to solve. In addition, we can use the quaternion method to compound the two rotations, so that the motion process of the rigid body is simplified to a translation and a rotation around a fixed axis. We take the quaternion and Representing the two rotations in the above process, is the quaternion representation of the point coordinates before rotation, is the quaternion representation of the point coordinates after the first rotation, The quaternion representation of the point coordinates after the second rotation, the specific rotation composite process is as follows:
就是复合之后的旋转四元数,由它可以获取旋转轴和旋转角θequal。最终的运动分析便是对选定的任意两个时刻之间的运动过程进行分析,给出简化运动过程中发生的平移的向量旋转的旋转轴以及逆时针旋转过的角度θequal。 It is the rotation quaternion after compounding, from which the rotation axis can be obtained and rotation angle θ equal . The final motion analysis is to analyze the motion process between any two selected moments, and give the vector that simplifies the translation that occurs during the motion process rotating axis of rotation And the angle θ equal rotated counterclockwise.
本发明有诸多应用途径,上述实施例仅是本发明的一种优选实施方式,因此本发明并不局限于上述实施方式。在不脱离本发明原理的情况下,本技术领域的技术人员可能在本发明的启示下设计出其他实施方式,这些也应视为本发明的保护范围。The present invention has many application ways, and the above-mentioned embodiment is only a preferred implementation of the present invention, so the present invention is not limited to the above-mentioned implementation. Without departing from the principles of the present invention, those skilled in the art may design other implementations under the inspiration of the present invention, and these should also be regarded as the protection scope of the present invention.
Claims (8)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810448583.2A CN108627153B (en) | 2018-05-11 | 2018-05-11 | Rigid body motion tracking system based on inertial sensor and working method thereof |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810448583.2A CN108627153B (en) | 2018-05-11 | 2018-05-11 | Rigid body motion tracking system based on inertial sensor and working method thereof |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108627153A true CN108627153A (en) | 2018-10-09 |
| CN108627153B CN108627153B (en) | 2022-06-07 |
Family
ID=63692625
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810448583.2A Active CN108627153B (en) | 2018-05-11 | 2018-05-11 | Rigid body motion tracking system based on inertial sensor and working method thereof |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108627153B (en) |
Cited By (23)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109343546A (en) * | 2018-10-10 | 2019-02-15 | 上海海洋大学 | A real-time measurement system for submersible six-degree-of-freedom motion |
| CN111307174A (en) * | 2020-02-05 | 2020-06-19 | 深圳前海达闼云端智能科技有限公司 | Calibration method of sensor, moving object and storage medium |
| CN111413995A (en) * | 2020-03-24 | 2020-07-14 | 北京科技大学 | Method and system for tracking relative position and synchronously controlling posture between double rigid body characteristic points |
| CN111536966A (en) * | 2020-04-07 | 2020-08-14 | 广州体育学院 | High-precision position detection sensor |
| CN111831959A (en) * | 2020-03-05 | 2020-10-27 | 北京嘀嘀无限科技发展有限公司 | Motion data processing method, motion data processing device, terminal and computer-readable storage medium |
| CN112248940A (en) * | 2020-08-27 | 2021-01-22 | 清华大学 | Sensor arrangement resolving method and system for vehicle body attitude control |
| CN112308110A (en) * | 2020-09-21 | 2021-02-02 | 哈尔滨工业大学(深圳) | Hand motion recognition method and system capable of realizing man-machine interaction |
| CN112489082A (en) * | 2020-12-03 | 2021-03-12 | 海宁奕斯伟集成电路设计有限公司 | Position detection method, position detection device, electronic equipment and readable storage medium |
| CN112557060A (en) * | 2020-12-23 | 2021-03-26 | 重庆交通大学 | Dynamic comfort human body three-dimensional sitting posture change detection system, method and terminal |
| CN112595338A (en) * | 2020-12-24 | 2021-04-02 | 中国联合网络通信集团有限公司 | Navigation method and navigation system |
| CN113358897A (en) * | 2020-03-04 | 2021-09-07 | 埃韦利克斯公司 | Sensor system for an actuator, actuator and method for moving an actuator part |
| CN113916228A (en) * | 2021-10-09 | 2022-01-11 | 台州学院 | A method for monitoring regional activity in young children |
| CN114139319A (en) * | 2021-12-13 | 2022-03-04 | 吉林大学 | Reconfigurable multifunctional numerical control machining module configuration analysis method |
| CN114562993A (en) * | 2022-02-28 | 2022-05-31 | 联想(北京)有限公司 | Track processing method and device and electronic equipment |
| CN114674310A (en) * | 2022-03-10 | 2022-06-28 | 北京理工大学 | An inertial positioning method based on near-zero velocity observation in running state |
| CN114796889A (en) * | 2021-01-21 | 2022-07-29 | 王全锋 | Precise radiotherapy system based on four-dimensional adaptive digital tracking conformal intensity modulated focusing medical electron accelerator |
| CN115131874A (en) * | 2022-06-29 | 2022-09-30 | 深圳市神州云海智能科技有限公司 | User behavior recognition prediction method and system and intelligent safety helmet |
| CN115372659A (en) * | 2022-08-05 | 2022-11-22 | 东莞市汇能电子有限公司 | A method for battery motion trajectory and speed calculation through a three-axis gyroscope |
| US20230045497A1 (en) * | 2021-08-03 | 2023-02-09 | Fujifilm Healthcare Corporation | Magnetic resonance imaging apparatus and control method thereof |
| CN116086442A (en) * | 2022-12-29 | 2023-05-09 | 北京元客视界科技有限公司 | Fusion positioning device of multiple positioning sensors comprising inertial sensor |
| CN116202515A (en) * | 2022-12-14 | 2023-06-02 | 北京诺亦腾科技有限公司 | Data detection method, device, equipment and medium based on compensatory action elimination |
| CN116228818A (en) * | 2023-03-27 | 2023-06-06 | 南京大学 | Linear rigid body three-dimensional space motion tracking system and method based on edge computing |
| CN116407822A (en) * | 2023-04-14 | 2023-07-11 | 河北师范大学 | Method, device and upper computer for processing information on skeet motion |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101367015A (en) * | 2007-06-18 | 2009-02-18 | 神奇游戏智能公司 | Self contained inertial navigation system for interactive control using movable controller |
| US20090087029A1 (en) * | 2007-08-22 | 2009-04-02 | American Gnc Corporation | 4D GIS based virtual reality for moving target prediction |
| CN101579238A (en) * | 2009-06-15 | 2009-11-18 | 吴健康 | Human motion capture three dimensional playback system and method thereof |
| CN101886927A (en) * | 2010-06-25 | 2010-11-17 | 武汉大学 | Three-dimensional motion tracking system and method based on inertial sensor and geomagnetic sensor |
| EP2419057B1 (en) * | 2009-04-13 | 2014-10-01 | DEKA Products Limited Partnership | System and apparatus for orientation control |
| CN105607760A (en) * | 2015-12-18 | 2016-05-25 | 上海开圣影视文化传媒股份有限公司 | Trace restoration method and system based on micro inertial sensor |
| CN106679649A (en) * | 2016-12-12 | 2017-05-17 | 浙江大学 | Hand movement tracking system and tracking method |
| CN106767770A (en) * | 2016-11-29 | 2017-05-31 | 西安交通大学 | A kind of detection of user's direction of travel and method for tracing based on portable intelligent equipment |
| CN107092882A (en) * | 2017-04-19 | 2017-08-25 | 南京大学 | A kind of Activity recognition system and its method of work perceived based on sub- action |
| CN107115114A (en) * | 2017-04-28 | 2017-09-01 | 王春宝 | Human Stamina evaluation method, apparatus and system |
| CN107421537A (en) * | 2017-09-14 | 2017-12-01 | 桂林电子科技大学 | Object motion attitude cognitive method and system based on inertial sensor rigid body grid |
| CN107560613A (en) * | 2017-08-15 | 2018-01-09 | 江苏科技大学 | Trajectory Tracking System and method in robot chamber based on nine axle inertial sensors |
-
2018
- 2018-05-11 CN CN201810448583.2A patent/CN108627153B/en active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101367015A (en) * | 2007-06-18 | 2009-02-18 | 神奇游戏智能公司 | Self contained inertial navigation system for interactive control using movable controller |
| US20090087029A1 (en) * | 2007-08-22 | 2009-04-02 | American Gnc Corporation | 4D GIS based virtual reality for moving target prediction |
| EP2419057B1 (en) * | 2009-04-13 | 2014-10-01 | DEKA Products Limited Partnership | System and apparatus for orientation control |
| CN101579238A (en) * | 2009-06-15 | 2009-11-18 | 吴健康 | Human motion capture three dimensional playback system and method thereof |
| CN101886927A (en) * | 2010-06-25 | 2010-11-17 | 武汉大学 | Three-dimensional motion tracking system and method based on inertial sensor and geomagnetic sensor |
| CN105607760A (en) * | 2015-12-18 | 2016-05-25 | 上海开圣影视文化传媒股份有限公司 | Trace restoration method and system based on micro inertial sensor |
| CN106767770A (en) * | 2016-11-29 | 2017-05-31 | 西安交通大学 | A kind of detection of user's direction of travel and method for tracing based on portable intelligent equipment |
| CN106679649A (en) * | 2016-12-12 | 2017-05-17 | 浙江大学 | Hand movement tracking system and tracking method |
| CN107092882A (en) * | 2017-04-19 | 2017-08-25 | 南京大学 | A kind of Activity recognition system and its method of work perceived based on sub- action |
| CN107115114A (en) * | 2017-04-28 | 2017-09-01 | 王春宝 | Human Stamina evaluation method, apparatus and system |
| CN107560613A (en) * | 2017-08-15 | 2018-01-09 | 江苏科技大学 | Trajectory Tracking System and method in robot chamber based on nine axle inertial sensors |
| CN107421537A (en) * | 2017-09-14 | 2017-12-01 | 桂林电子科技大学 | Object motion attitude cognitive method and system based on inertial sensor rigid body grid |
Non-Patent Citations (4)
| Title |
|---|
| TADEJ BERAVS等: "Development and validation of a wearable inertial measurement system for use with lower limb exoskeletons", 《2011 11TH IEEE-RAS INTERNATIONAL CONFERENCE ON HUMANOID ROBOTS》, 28 October 2011 (2011-10-28), pages 212 - 217, XP032070371, DOI: 10.1109/Humanoids.2011.6100914 * |
| 吕佐臣等: "《飞航导弹发射装置 下》", 29 February 1996, 宇航出版社, pages: 241 - 243 * |
| 陈靖等: "适用于户外增强现实系统的混合跟踪定位算法", 《计算机辅助设计与图形学学报》, no. 02, 15 February 2010 (2010-02-15) * |
| 马强等: "利用惯性传感器的人体运动空间轨迹追踪", 《单片机与嵌入式系统应用》, no. 01, 1 January 2018 (2018-01-01) * |
Cited By (32)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109343546A (en) * | 2018-10-10 | 2019-02-15 | 上海海洋大学 | A real-time measurement system for submersible six-degree-of-freedom motion |
| CN111307174A (en) * | 2020-02-05 | 2020-06-19 | 深圳前海达闼云端智能科技有限公司 | Calibration method of sensor, moving object and storage medium |
| US11896413B2 (en) | 2020-03-04 | 2024-02-13 | Ewellix AB | Sensor system for an actuator, actuator and method for moving an actuator part |
| CN113358897A (en) * | 2020-03-04 | 2021-09-07 | 埃韦利克斯公司 | Sensor system for an actuator, actuator and method for moving an actuator part |
| CN111831959A (en) * | 2020-03-05 | 2020-10-27 | 北京嘀嘀无限科技发展有限公司 | Motion data processing method, motion data processing device, terminal and computer-readable storage medium |
| CN111831959B (en) * | 2020-03-05 | 2024-06-25 | 北京嘀嘀无限科技发展有限公司 | Motion data processing method, device, terminal and computer readable storage medium |
| CN111413995A (en) * | 2020-03-24 | 2020-07-14 | 北京科技大学 | Method and system for tracking relative position and synchronously controlling posture between double rigid body characteristic points |
| CN111413995B (en) * | 2020-03-24 | 2021-03-05 | 北京科技大学 | Method and system for relative position tracking and attitude synchronization control between dual rigid body feature points |
| CN111536966A (en) * | 2020-04-07 | 2020-08-14 | 广州体育学院 | High-precision position detection sensor |
| CN112248940A (en) * | 2020-08-27 | 2021-01-22 | 清华大学 | Sensor arrangement resolving method and system for vehicle body attitude control |
| CN112308110B (en) * | 2020-09-21 | 2023-07-21 | 哈尔滨工业大学(深圳) | A hand movement recognition method and system capable of realizing human-computer interaction |
| CN112308110A (en) * | 2020-09-21 | 2021-02-02 | 哈尔滨工业大学(深圳) | Hand motion recognition method and system capable of realizing man-machine interaction |
| CN112489082A (en) * | 2020-12-03 | 2021-03-12 | 海宁奕斯伟集成电路设计有限公司 | Position detection method, position detection device, electronic equipment and readable storage medium |
| CN112557060A (en) * | 2020-12-23 | 2021-03-26 | 重庆交通大学 | Dynamic comfort human body three-dimensional sitting posture change detection system, method and terminal |
| CN112595338A (en) * | 2020-12-24 | 2021-04-02 | 中国联合网络通信集团有限公司 | Navigation method and navigation system |
| CN114796889A (en) * | 2021-01-21 | 2022-07-29 | 王全锋 | Precise radiotherapy system based on four-dimensional adaptive digital tracking conformal intensity modulated focusing medical electron accelerator |
| US12306281B2 (en) | 2021-08-03 | 2025-05-20 | Fujifilm Corporation | Magnetic resonance imaging apparatus and control method thereof |
| US20230045497A1 (en) * | 2021-08-03 | 2023-02-09 | Fujifilm Healthcare Corporation | Magnetic resonance imaging apparatus and control method thereof |
| US11927658B2 (en) * | 2021-08-03 | 2024-03-12 | Fujifilm Healthcare Corporation | Magnetic resonance imaging apparatus and control method thereof |
| CN113916228A (en) * | 2021-10-09 | 2022-01-11 | 台州学院 | A method for monitoring regional activity in young children |
| CN114139319A (en) * | 2021-12-13 | 2022-03-04 | 吉林大学 | Reconfigurable multifunctional numerical control machining module configuration analysis method |
| CN114139319B (en) * | 2021-12-13 | 2024-04-16 | 吉林大学 | Reconfigurable multifunctional numerical control processing module configuration analysis method |
| CN114562993A (en) * | 2022-02-28 | 2022-05-31 | 联想(北京)有限公司 | Track processing method and device and electronic equipment |
| CN114674310A (en) * | 2022-03-10 | 2022-06-28 | 北京理工大学 | An inertial positioning method based on near-zero velocity observation in running state |
| CN114674310B (en) * | 2022-03-10 | 2025-07-08 | 北京理工大学 | Inertial positioning method based on near zero speed observation under running state |
| CN115131874B (en) * | 2022-06-29 | 2023-10-17 | 深圳市神州云海智能科技有限公司 | User behavior recognition prediction method, system and intelligent safety helmet |
| CN115131874A (en) * | 2022-06-29 | 2022-09-30 | 深圳市神州云海智能科技有限公司 | User behavior recognition prediction method and system and intelligent safety helmet |
| CN115372659A (en) * | 2022-08-05 | 2022-11-22 | 东莞市汇能电子有限公司 | A method for battery motion trajectory and speed calculation through a three-axis gyroscope |
| CN116202515A (en) * | 2022-12-14 | 2023-06-02 | 北京诺亦腾科技有限公司 | Data detection method, device, equipment and medium based on compensatory action elimination |
| CN116086442A (en) * | 2022-12-29 | 2023-05-09 | 北京元客视界科技有限公司 | Fusion positioning device of multiple positioning sensors comprising inertial sensor |
| CN116228818A (en) * | 2023-03-27 | 2023-06-06 | 南京大学 | Linear rigid body three-dimensional space motion tracking system and method based on edge computing |
| CN116407822A (en) * | 2023-04-14 | 2023-07-11 | 河北师范大学 | Method, device and upper computer for processing information on skeet motion |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108627153B (en) | 2022-06-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108627153A (en) | A kind of rigid motion tracing system and its working method based on inertial sensor | |
| CN111462231B (en) | Positioning method based on RGBD sensor and IMU sensor | |
| CN109885080B (en) | Autonomous control system and autonomous control method | |
| CN107314778B (en) | Calibration method, device and system for relative attitude | |
| CN111091587B (en) | Low-cost motion capture method based on visual markers | |
| CN115479598A (en) | Positioning and mapping method based on multi-sensor fusion and tight coupling system | |
| CN105953796A (en) | Stable motion tracking method and stable motion tracking device based on integration of simple camera and IMU (inertial measurement unit) of smart cellphone | |
| CN107888828A (en) | Space-location method and device, electronic equipment and storage medium | |
| CN106708066A (en) | Autonomous landing method of unmanned aerial vehicle based on vision/inertial navigation | |
| CN110553648A (en) | method and system for indoor navigation | |
| JP6288858B2 (en) | Method and apparatus for estimating position of optical marker in optical motion capture | |
| CN114001736A (en) | Positioning method, positioning device, storage medium and electronic equipment | |
| CN114608554A (en) | Handheld SLAM equipment and robot instant positioning and mapping method | |
| CN115218906A (en) | Indoor SLAM-oriented visual inertial fusion positioning method and system | |
| CN119334346A (en) | A navigation and positioning method, medium and device for electric power inspection drone | |
| CN108917755B (en) | Imaging seeker line-of-sight angle zero error estimation method and device | |
| CN111811421B (en) | High-speed real-time deformation monitoring method and system | |
| CN112837374B (en) | A spatial positioning method and system | |
| CN111145267A (en) | IMU (inertial measurement unit) assistance-based 360-degree panoramic view multi-camera calibration method | |
| CN114964217B (en) | A method and device for estimating state information | |
| CN116678391A (en) | Calibration method and storage medium for mobile device and its multi-sensors | |
| JP6205387B2 (en) | Method and apparatus for acquiring position information of virtual marker, and operation measurement method | |
| CN116465396A (en) | A Cooperative Mapping Method for Multi-robots | |
| CN114429515A (en) | A point cloud map construction method, device and device | |
| CN117437290B (en) | Multi-sensor fusion type three-dimensional space positioning method for unmanned aerial vehicle in natural protection area |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |