WO2021163846A1 - Target tracking method and target tracking apparatus - Google Patents

Target tracking method and target tracking apparatus Download PDF

Info

Publication number
WO2021163846A1
WO2021163846A1 PCT/CN2020/075556 CN2020075556W WO2021163846A1 WO 2021163846 A1 WO2021163846 A1 WO 2021163846A1 CN 2020075556 W CN2020075556 W CN 2020075556W WO 2021163846 A1 WO2021163846 A1 WO 2021163846A1
Authority
WO
WIPO (PCT)
Prior art keywords
time
state quantity
observation
updated
update
Prior art date
Application number
PCT/CN2020/075556
Other languages
French (fr)
Chinese (zh)
Inventor
周鹏
冯源
张欢
李选富
吴祖光
李维
Original Assignee
华为技术有限公司
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 华为技术有限公司 filed Critical 华为技术有限公司
Priority to PCT/CN2020/075556 priority Critical patent/WO2021163846A1/en
Priority to CN202080095340.0A priority patent/CN115039095A/en
Publication of WO2021163846A1 publication Critical patent/WO2021163846A1/en

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/17Function evaluation by approximation methods, e.g. inter- or extrapolation, smoothing, least mean square method

Definitions

  • the embodiments of the present application relate to the field of artificial intelligence, and more specifically, to a target tracking method and a target tracking device.
  • the target tracking method can be applied to many scenarios, such as target detection, lane line detection, positioning and so on.
  • Common target tracking methods need to use sensor measurement data combined with the state transition model of the tracked object to perform Bayesian estimation.
  • Bayes filters can be used to achieve target tracking.
  • the measurement data of the sensor is input to the Bayesian filter to obtain an estimate of the state of the tracked object.
  • Fig. 1 is a schematic diagram of the measurement data input to the Bayesian filter appearing out of order in time. The abscissa in Figure 1 shows the real time.
  • the data transmitted by the sensor includes a timestamp (header: Xms) and measurement data (data).
  • the timestamp is used to indicate that the measurement or collection time of the measurement data is Xms. time.
  • the measurement data measured by sensor 0 at 1ms is input into the Bayesian filter at 3ms, and the measurement data measured by sensor 1 at 0ms is input into the Bayesian filter at 4ms, resulting in the time to input the measurement data of the Bayesian filter. Out of order.
  • the Bayesian filter first outputs the state quantity calculated based on the measurement data measured at 1ms with a time stamp of 1ms, and then outputs the state calculated based on the measurement data measured at 0ms with a time stamp of 0ms As a result, the time stamp of the output state quantity jumps with the change of the time stamp of the measurement data.
  • One method is to determine whether the time stamp of the measurement data is later than the time stamp of the latest updated state quantity in the Bayesian filter, and only the measurement data whose collection time is later than the update time of the latest state quantity is input to the Bayesian filter. Yes filter.
  • this method will directly cause the lack of observations, leading to a decrease in the confidence of the output results.
  • Another method is to extrapolate the late observation data to the current moment. However, extrapolation in the observation space relies on strong assumptions and tends to enlarge errors.
  • the present application provides a target tracking method and a target tracking device, which can realize target tracking when the acquisition time of the observation is out of order, and avoid the out-of-order jump of the time stamp of the state quantity with the out-of-order change of the time stamp of the observation. .
  • a target tracking method which includes: acquiring the actual observation of the moving object collected by the sensor; in the case that the acquisition time of the actual observation is earlier than the update time of the current state of the moving object, according to the moving object
  • the corresponding state transition model obtains the estimated observation; the current state quantity is updated according to the actual observation and the estimated observation, and the updated state quantity is used as the state quantity corresponding to the update moment.
  • Observation can include the speed of the moving object, the acceleration of the moving object, or the position of the moving object.
  • the observation may include the position and/or speed of the moving object measured by the millimeter wave radar.
  • the state quantity may include the speed of the moving object, the acceleration of the moving object, or the position of the moving object.
  • the update time of the current state quantity of the moving object can be understood as the time when the state quantity of the moving object was updated last time.
  • the state transition model of the moving object includes a kinematics model of the moving object.
  • the kinematics model may include a uniform linear motion model and the like.
  • the state transition is performed according to the state transition model corresponding to the moving object, and then the current state quantity is updated, so that a more accurate state quantity corresponding to the latest update time can be obtained.
  • the updated state quantity is the state quantity corresponding to the latest update time, rather than the state quantity corresponding to the collection time, so as to avoid the out-of-order jump of the output state quantity in time.
  • obtaining estimated observations according to the state transition model corresponding to the moving object includes: obtaining the likelihood of the observation corresponding to the acquisition time according to the state transition model corresponding to the moving object Function, the estimated observation is determined according to the likelihood function of the observation corresponding to the acquisition time, the state transition model is determined according to the first probability density function, the first probability density function is to transfer the current state quantity to the corresponding acquisition time The probability density function of the state quantity.
  • the current state quantity is updated based on actual observations and estimated observations, and the updated state quantity is used as the state quantity corresponding to the update time, including: The likelihood function of the observation corresponding to the time updates the probability density function of the current state quantity; the state quantity corresponding to the update time is determined according to the probability density function of the updated state quantity.
  • the likelihood function of the observation corresponding to the acquisition time satisfies:
  • X k ') represents the likelihood function of the observation z k-1 corresponding to the time t k-1
  • X represents the state quantity updated at the time t k-1
  • X) Z represents Z k-1 k-1 is used on the time t where k-1 is updated state quantity
  • X k ') is the first probability density function, the first probability density function is to transfer the state quantity X k ' updated at t k to update at t k-1 The probability density function of the state quantity.
  • the probability density function of the updated state quantity satisfies:
  • z k ) represents the probability density function of the state quantity X k updated at time t k
  • z k ,z k-2 ) represents the absence of z
  • z k represents the observations collected at time t k
  • z k represents the set of observations collected at k time ⁇ z 1 , z 2 ,...,z k-2 ,z k-1 ,z k ⁇
  • z k-2 represents the set of observations collected at k-2 moments ⁇ z 1 ,z 2 ,...,z k-2 ⁇ .
  • the current state quantity is updated based on the actual observations and estimated observations, and the updated state quantity is used as the state quantity corresponding to the update time, including:
  • the expectation of the state quantity corresponding to the update time is regarded as the state quantity corresponding to the update time
  • the expectation of the state quantity corresponding to the update time is related to the Kalman gain value.
  • the Kalman gain value is related to the first covariance, the observation matrix at the acquisition time, the covariance of the observation matrix, and the variance of the observation at the acquisition time.
  • a covariance refers to the covariance of the state quantity transferred from the update time to the collection time.
  • the Kalman gain value satisfies:
  • k represents t K transferred from the update time to the acquisition time t covariance k-1 state quantity
  • k represents the lack Concept t k-1 time captured measurement Z k-1
  • the covariance of the state quantity updated at time t k F k-1
  • Q k represents the covariance of the prediction matrix
  • H k-1 Represents the observation matrix at time t k-1
  • Var(z k-1 ) represents the variance of the observation z k-1 collected at time t k-1
  • R k-1 represents the covariance of the observation matrix.
  • the expectation of the state quantity corresponding to the update time meets:
  • k represents the expectation of the state quantity updated at time t k
  • the estimated value of the observation corresponding to the collection time satisfies:
  • H k-1 represents a t the observation matrix k-1 time
  • k indicates a transition from t K time to t k-1 time of the state transition matrix
  • k represents the concept of a lack of time t k-1 measurement state quantity acquired in the case where time t k z k-1 update is desired.
  • the estimated observation in the case where the acquisition time of the actual observation is earlier than the update time of the current state of the moving object, the estimated observation is obtained according to the state transition model corresponding to the moving object , Including: when the acquisition time of the observation is earlier than the update time of the current state quantity of the moving object, and the time difference between the acquisition time and the update time is less than or equal to the threshold value, the estimated observation is obtained according to the state transition model corresponding to the moving object quantity.
  • the state transition is performed only when the time difference between the collection time and the update time is less than or equal to the threshold, and then the state quantity is updated, which can avoid the use of This observation updates the state quantity, which reduces the confidence of the updated state quantity.
  • a target tracking device which includes an acquisition module and a processing module, wherein the acquisition module is used to: acquire the actual observation of a moving object collected by the sensor; the processing module is used to: early in the acquisition time of the actual observation In the case of the update time of the current state quantity of the moving object, the estimated observation quantity is obtained according to the state transition model corresponding to the moving object; the current state quantity is updated according to the actual observation quantity and the estimated observation quantity, and the updated state quantity is used as the update The state quantity corresponding to the moment.
  • the state transition is performed according to the state transition model corresponding to the moving object, and then the current state quantity is updated, so that a more accurate state quantity corresponding to the latest update time can be obtained.
  • the updated state quantity is the state quantity corresponding to the latest update time, rather than the state quantity corresponding to the collection time, so as to avoid the out-of-order jump of the output state quantity in time.
  • the processing module is used to obtain the likelihood function of the observation corresponding to the acquisition time according to the state transition model corresponding to the moving object, and the estimated observation is based on the acquisition time
  • the likelihood function of the corresponding observation is determined, and the state transition model is determined according to the first probability density function.
  • the first probability density function is the probability density function that transfers the current state quantity to the state quantity corresponding to the collection time.
  • the processing module is used to: update the probability density function of the current state quantity according to the likelihood function of the observation corresponding to the acquisition time; according to the updated state quantity The probability density function of determines the state quantity corresponding to the update moment.
  • the likelihood function of the observation corresponding to the acquisition time satisfies:
  • X k ') represents the likelihood function of the observation z k-1 corresponding to the time t k-1
  • X represents the state quantity updated at the time t k-1
  • X) Z represents Z k-1 k-1 is used on the time t where k-1 is updated state quantity
  • X k ') is the first probability density function, the first probability density function is to transfer the state quantity X k ' updated at t k to update at t k-1 The probability density function of the state quantity.
  • the probability density function of the updated state quantity satisfies:
  • z k ) represents the probability density function of the state quantity X k updated at time t k
  • z k ,z k-2 ) represents the absence of z
  • z k represents the observations collected at time t k
  • z k represents the set of observations collected at k time ⁇ z 1 , z 2 ,...,z k-2 ,z k-1 ,z k ⁇
  • z k-2 represents the set of observations collected at k-2 moments ⁇ z 1 ,z 2 ,...,z k-2 ⁇ .
  • the processing module is used to: determine the expectation of the state quantity corresponding to the update time according to the actual observation and the estimated observation; and use the expectation of the state quantity corresponding to the update time as the update
  • the state quantity corresponding to the time; among them, the expectation of the state quantity corresponding to the update time is related to the Kalman gain value, the Kalman gain value and the first covariance, the observation matrix at the acquisition time, the covariance of the observation matrix and the observation at the acquisition time
  • the variance of the quantity is related, and the first covariance refers to the covariance of the state quantity transferred from the update time to the collection time.
  • the Kalman gain value satisfies:
  • k represents t K transferred from the update time to the acquisition time t covariance k-1 state quantity
  • k represents the lack Concept t k-1 time captured measurement Z k-1
  • the covariance of the state quantity updated at time t k F k-1
  • Q k represents the covariance of the prediction matrix
  • H k-1 Represents the observation matrix at time t k-1
  • Var(z k-1 ) represents the variance of the observation z k-1 collected at time t k-1
  • R k-1 represents the covariance of the observation matrix.
  • the expectation of the state quantity corresponding to the update time meets:
  • k represents the expectation of the state quantity updated at time t k
  • the estimated value of the observation corresponding to the collection time satisfies:
  • H k-1 represents a t the observation matrix k-1 time
  • k indicates a transition from t K time to t k-1 time of the state transition matrix
  • k represents the concept of a lack of time t k-1 measurement state quantity acquired in the case where time t k z k-1 update is desired.
  • the processing module is used to: the actual observation measurement acquisition time is earlier than the update time of the current state quantity of the moving object, and the time difference between the acquisition time and the update time When the value is less than or equal to the threshold, the estimated observation is obtained according to the state transition model corresponding to the moving object.
  • a target tracking device in a third aspect, includes: a memory for storing a program; a processor for executing the program stored in the memory, and when the program stored in the memory is executed, the processor is configured to execute the first aspect In the method.
  • a computer program product includes: computer program code, which when the computer program product runs on a computer, causes the computer to execute the method in the first aspect.
  • a computer-readable storage medium stores a computer program, and when the computer program runs on a computer, the computer executes the method in the first aspect.
  • the method of the first aspect may specifically refer to the first aspect and a method in any one of the various implementation manners of the first aspect.
  • Figure 1 is a schematic diagram of measurement data input into Bayesian filtering
  • Fig. 2 is a schematic structural diagram of a vehicle provided by an embodiment of the present application.
  • Fig. 3 is a schematic structural diagram of a computer system provided by an embodiment of the present application.
  • FIG. 4 is a schematic diagram of the application of a cloud-side command automatic driving vehicle provided by an embodiment of the present application
  • FIG. 5 is a schematic structural diagram of a target tracking device provided by an embodiment of the present application.
  • FIG. 6 is a schematic flowchart of a target tracking method provided by an embodiment of the present application.
  • FIG. 7 is a schematic structural diagram of another target tracking device provided by an embodiment of the present application.
  • FIG. 8 is a schematic structural diagram of another target tracking device provided by an embodiment of the present application.
  • Fig. 2 is a functional block diagram of a vehicle 100 provided by an embodiment of the present application.
  • the vehicle 100 may be a manually driven vehicle, or the vehicle 100 may be configured in a fully or partially automatic driving mode.
  • the vehicle 100 can control its own vehicle while in the automatic driving mode, and can determine the current state of the vehicle and its surrounding environment through human operations, determine the possible behavior of at least one other vehicle in the surrounding environment, and The confidence level corresponding to the possibility of other vehicles performing possible behaviors is determined, and the vehicle 100 is controlled based on the determined information.
  • the vehicle 100 can be placed to operate without human interaction.
  • the vehicle 100 may include various subsystems, such as a traveling system 110, a sensing system 120, a control system 130, one or more peripheral devices 140 and a power supply 160, a computer system 150, and a user interface 170.
  • a traveling system 110 a sensing system 120
  • a control system 130 a control system 130
  • peripheral devices 140 and a power supply 160 a computer system 150
  • a user interface 170 a user interface 170.
  • the vehicle 100 may include more or fewer subsystems, and each subsystem may include multiple elements.
  • each of the subsystems and elements of the vehicle 100 may be wired or wirelessly interconnected.
  • the travel system 110 may include components for providing power movement to the vehicle 100.
  • the travel system 110 may include an engine 111, a transmission 112, an energy source 113, and wheels 114/tires.
  • the engine 111 may be an internal combustion engine, an electric motor, an air compression engine, or other types of engine combinations; for example, a hybrid engine composed of a gasoline engine and an electric motor, or a hybrid engine composed of an internal combustion engine and an air compression engine.
  • the engine 111 can convert the energy source 113 into mechanical energy.
  • the energy source 113 may include gasoline, diesel, other petroleum-based fuels, propane, other compressed gas-based fuels, ethanol, solar panels, batteries, and other power sources.
  • the energy source 113 may also provide energy for other systems of the vehicle 100.
  • the transmission device 112 may include a gearbox, a differential, and a drive shaft; wherein, the transmission device 112 may transmit mechanical power from the engine 111 to the wheels 114.
  • the transmission device 112 may also include other devices, such as a clutch.
  • the drive shaft may include one or more shafts that can be coupled to one or more wheels 114.
  • the sensing system 120 may include several sensors that sense information about the environment around the vehicle 100.
  • the sensing system 120 may include a positioning system 121 (for example, a GPS system, a Beidou system or other positioning systems), an inertial measurement unit 122 (IMU), a radar 123, a laser rangefinder 124, and a camera 125.
  • the sensing system 120 may also include sensors of the internal system of the monitored vehicle 100 (for example, an in-vehicle air quality monitor, a fuel gauge, an oil temperature gauge, etc.). Sensor data from one or more of these sensors can be used to detect objects and their corresponding characteristics (position, shape, direction, speed, etc.). Such detection and identification are key functions for the safe operation of the autonomous vehicle 100.
  • the positioning system 121 can be used to estimate the geographic location of the vehicle 100.
  • the IMU 122 may be used to sense changes in the position and orientation of the vehicle 100 based on inertial acceleration.
  • the IMU 122 may be a combination of an accelerometer and a gyroscope.
  • the radar 123 may use radio signals to sense objects in the surrounding environment of the vehicle 100. In some embodiments, in addition to sensing the object, the radar 123 may also be used to sense the speed and/or direction of the object.
  • the laser rangefinder 124 may use laser light to sense objects in the environment where the vehicle 100 is located.
  • the laser rangefinder 124 may include one or more laser sources, laser scanners, and one or more detectors, as well as other system components.
  • the camera 125 may be used to capture multiple images of the surrounding environment of the vehicle 100.
  • the camera 125 may be a still camera or a video camera.
  • control system 130 controls the operation of the vehicle 100 and its components.
  • the control system 130 may include various elements, such as a steering system 131, a throttle 132, a braking unit 133, a computer vision system 134, a route control system 135, and an obstacle avoidance system 136.
  • the steering system 131 may be operated to adjust the forward direction of the vehicle 100.
  • it may be a steering wheel system in one embodiment.
  • the throttle 132 may be used to control the operating speed of the engine 111 and thereby control the speed of the vehicle 100.
  • the braking unit 133 may be used to control the deceleration of the vehicle 100; the braking unit 133 may use friction to slow down the wheels 114. In other embodiments, the braking unit 133 may convert the kinetic energy of the wheels 114 into electric current. The braking unit 133 may also take other forms to slow down the rotation speed of the wheels 114 to control the speed of the vehicle 100.
  • the computer vision system 134 may be operable to process and analyze the images captured by the camera 125 in order to identify objects and/or features in the surrounding environment of the vehicle 100.
  • the aforementioned objects and/or features may include traffic signals, road boundaries and obstacles.
  • the computer vision system 134 may use object recognition algorithms, structure from motion (SFM) algorithms, video tracking, and other computer vision technologies.
  • the computer vision system 134 may be used to map the environment, track objects, estimate the speed of objects, and so on.
  • the route control system 135 may be used to determine the travel route of the vehicle 100.
  • the route control system 135 may combine data from sensors, GPS, and one or more predetermined maps to determine a travel route for the vehicle 100.
  • the obstacle avoidance system 136 may be used to identify, evaluate, and avoid or otherwise cross potential obstacles in the environment of the vehicle 100.
  • control system 130 may additionally or alternatively include components other than those shown and described. Alternatively, a part of the components shown above may be reduced.
  • the vehicle 100 can interact with external sensors, other vehicles, other computer systems, or users through a peripheral device 140; wherein, the peripheral device 140 can include a wireless communication system 141, an onboard computer 142, a microphone 143 and/ Or speaker 144.
  • the peripheral device 140 can include a wireless communication system 141, an onboard computer 142, a microphone 143 and/ Or speaker 144.
  • the peripheral device 140 may provide a means for the vehicle 100 to interact with the user interface 170.
  • the onboard computer 142 may provide information to the user of the vehicle 100.
  • the user interface 116 can also operate the onboard computer 142 to receive user input; the onboard computer 142 can be operated through a touch screen.
  • the peripheral device 140 may provide a means for the vehicle 100 to communicate with other devices located in the vehicle.
  • the microphone 143 may receive audio (eg, voice commands or other audio input) from the user of the vehicle 100.
  • the speaker 144 may output audio to the user of the vehicle 100.
  • the wireless communication system 141 may wirelessly communicate with one or more devices directly or via a communication network.
  • the wireless communication system 141 can use 3G cellular communication; for example, code division multiple access (CDMA), EVD0, global system for mobile communications (GSM)/general packet radio service (general packet radio service) packet radio service, GPRS), or 4G cellular communication, such as long term evolution (LTE); or, 5G cellular communication.
  • CDMA code division multiple access
  • EVD0 global system for mobile communications
  • GSM global system for mobile communications
  • general packet radio service general packet radio service
  • GPRS general packet radio service
  • 4G cellular communication such as long term evolution (LTE)
  • LTE long term evolution
  • 5G cellular communication 5G cellular communication.
  • the wireless communication system 141 can communicate with a wireless local area network (WLAN) by using wireless Internet access (WiFi).
  • WiFi wireless Internet access
  • the wireless communication system 141 may directly communicate with the device using an infrared link, Bluetooth, or ZigBee; other wireless protocols, such as various vehicle communication systems, for example, the wireless communication system 141 may include one or Multiple dedicated short range communications (DSRC) devices, these devices may include public and/or private data communications between vehicles and/or roadside stations.
  • DSRC dedicated short range communications
  • the power supply 160 may provide power to various components of the vehicle 100.
  • the power source 160 may be a rechargeable lithium ion or lead-acid battery.
  • One or more battery packs of such batteries may be configured as a power source to provide power to various components of the vehicle 100.
  • the power source 160 and the energy source 113 may be implemented together, such as in some all-electric vehicles.
  • part or all of the functions of the vehicle 100 may be controlled by the computer system 150, where the computer system 150 may include at least one processor 151, and the processor 151 is executed in a non-transitory computer readable medium stored in the memory 152, for example.
  • the computer system 150 may also be multiple computing devices that control individual components or subsystems of the vehicle 100 in a distributed manner.
  • the processor 151 may be any conventional processor, such as a commercially available CPU.
  • the processor may be a dedicated device such as an ASIC or other hardware-based processor.
  • FIG. 2 functionally illustrates the processor, the memory, and other elements of the computer in the same block, those of ordinary skill in the art should understand that the processor, computer, or memory may or may not actually include Multiple processors, computers or memories in the same physical enclosure.
  • the memory may be a hard disk drive or other storage medium located in a housing other than the computer. Therefore, a reference to a processor or computer will be understood to include a reference to a collection of processors or computers or memories that may or may not operate in parallel. Rather than using a single processor to perform the steps described here, some components such as steering components and deceleration components may each have its own processor that only performs calculations related to component-specific functions .
  • the processor may be located away from the vehicle and wirelessly communicate with the vehicle.
  • some of the processes described herein are executed on a processor disposed in the vehicle and others are executed by a remote processor, including taking the necessary steps to perform a single manipulation.
  • the memory 152 may contain instructions 153 (eg, program logic), which may be executed by the processor 151 to perform various functions of the vehicle 100, including those functions described above.
  • the memory 152 may also contain additional instructions, for example, including sending data to, receiving data from, interacting with, and/or performing data to one or more of the traveling system 110, the sensing system 120, the control system 130, and the peripheral device 140. Control instructions.
  • the memory 152 may also store data, such as road maps, route information, the position, direction, and speed of the vehicle, and other such vehicle data, as well as other information. Such information may be used by the vehicle 100 and the computer system 150 during the operation of the vehicle 100 in autonomous, semi-autonomous, and/or manual modes.
  • the user interface 170 may be used to provide information to or receive information from a user of the vehicle 100.
  • the user interface 170 may include one or more input/output devices in the set of peripheral devices 140, for example, a wireless communication system 141, a car computer 142, a microphone 143, and a speaker 144.
  • the computer system 150 may control the functions of the vehicle 100 based on inputs received from various subsystems (for example, the traveling system 110, the sensing system 120, and the control system 130) and from the user interface 170.
  • the computer system 150 may use input from the control system 130 in order to control the braking unit 133 to avoid obstacles detected by the sensing system 120 and the obstacle avoidance system 136.
  • the computer system 150 is operable to provide control of many aspects of the vehicle 100 and its subsystems.
  • one or more of these components described above may be installed or associated with the vehicle 100 separately.
  • the storage 152 may exist partially or completely separately from the vehicle 100.
  • the above-mentioned components may be communicatively coupled together in a wired and/or wireless manner.
  • FIG. 2 should not be construed as a limitation to the embodiment of the present application.
  • the vehicle 100 may be an autonomous vehicle traveling on a road, and may recognize objects in its surrounding environment to determine the adjustment to the current speed.
  • the object may be other vehicles, traffic control equipment, or other types of objects.
  • each recognized object can be considered independently, and based on the respective characteristics of the object, such as its current speed, acceleration, distance from the vehicle, etc., can be used to determine the speed to be adjusted by the self-driving car.
  • the vehicle 100 or a computing device associated with the vehicle 100 may be based on the characteristics of the identified object and the state of the surrounding environment (for example, traffic, Rain, ice on the road, etc.) to predict the behavior of the identified object.
  • each recognized object depends on each other's behavior. Therefore, all recognized objects can also be considered together to predict the behavior of a single recognized object.
  • the vehicle 100 can adjust its speed based on the predicted behavior of the identified object.
  • the self-driving car can determine based on the predicted behavior of the object that the vehicle will need to be adjusted (e.g., accelerate, decelerate, or stop) to a stable state.
  • other factors may also be considered to determine the speed of the vehicle 100, such as the lateral position of the vehicle 100 on the road on which it is traveling, the curvature of the road, the proximity of static and dynamic objects, and so on.
  • the computing device can also provide instructions to modify the steering angle of the vehicle 100 so that the self-driving car follows a given trajectory and/or maintains an object near the self-driving car (for example, , The safe horizontal and vertical distances of cars in adjacent lanes on the road.
  • the above-mentioned vehicle 100 may be a car, truck, motorcycle, bus, boat, airplane, helicopter, lawn mower, recreational vehicle, playground vehicle, construction equipment, tram, golf cart, train, and trolley, etc.
  • the application examples are not particularly limited.
  • the vehicle 100 shown in FIG. 2 may be an automatic driving vehicle, and the automatic driving system will be described in detail below.
  • Fig. 3 is a schematic diagram of an automatic driving system provided by an embodiment of the present application.
  • the automatic driving system shown in FIG. 3 includes a computer system 201, where the computer system 201 includes a processor 203, and the processor 203 is coupled to a system bus 205.
  • the processor 203 may be one or more processors, where each processor may include one or more processor cores.
  • the display adapter 207 (video adapter) can drive the display 209, and the display 209 is coupled to the system bus 205.
  • the system bus 205 may be coupled to an input/output (I/O) bus 213 through a bus bridge 211, and an I/O interface 215 is coupled to an I/O bus.
  • I/O input/output
  • the I/O interface 215 communicates with a variety of I/O devices, such as input devices 217 (such as keyboard, mouse, touch screen, etc.), media tray 221, (such as CD-ROM, multimedia interface, etc.) .
  • the transceiver 223 can send and/or receive radio communication signals, and the camera 255 can capture landscape and dynamic digital video images.
  • the interface connected to the I/O interface 215 may be the USB port 225.
  • the processor 203 may be any traditional processor, such as a reduced instruction set computer (RISC) processor, a complex instruction set computer (CISC) processor, or a combination of the foregoing.
  • RISC reduced instruction set computer
  • CISC complex instruction set computer
  • the processor 203 may be a dedicated device such as an application specific integrated circuit (ASIC); the processor 203 may be a neural network processor or a combination of a neural network processor and the foregoing traditional processors.
  • ASIC application specific integrated circuit
  • the computer system 201 may be located far away from the autonomous driving vehicle, and may wirelessly communicate with the autonomous driving vehicle.
  • some of the processes described herein are executed on a processor provided in an autonomous vehicle, and others are executed by a remote processor, including taking actions required to perform a single manipulation.
  • the computer system 201 can communicate with the software deployment server 249 through the network interface 229.
  • the network interface 229 may be a hardware network interface, such as a network card.
  • the network 227 may be an external network, such as the Internet, or an internal network, such as an Ethernet or a virtual private network (VPN).
  • the network 227 may also be a wireless network, such as a wifi network, a cellular network, and so on.
  • the hard disk drive interface is coupled with the system bus 205
  • the hardware drive interface 231 can be connected with the hard drive 233
  • the system memory 235 is coupled with the system bus 205.
  • the data running in the system memory 235 may include an operating system 237 and application programs 243.
  • the operating system 237 may include a parser 239 (shell) and a kernel 241 (kernel).
  • the shell 239 is an interface between the user and the kernel of the operating system.
  • the shell can be the outermost layer of the operating system; the shell can manage the interaction between the user and the operating system, for example, waiting for the user's input, interpreting the user's input to the operating system, and processing various operating systems The output result.
  • the kernel 241 may be composed of those parts of the operating system that are used to manage memory, files, peripherals, and system resources. Directly interact with the hardware.
  • the operating system kernel usually runs processes and provides inter-process communication, providing CPU time slice management, interrupts, memory management, IO management, and so on.
  • Application programs 243 include programs that control auto-driving cars, such as programs that manage the interaction between autonomous vehicles and obstacles on the road, programs that control the route or speed of autonomous vehicles, and programs that control interaction between autonomous vehicles and other autonomous vehicles on the road. .
  • the application program 243 also exists on the system of the software deployment server 249. In one embodiment, the computer system 201 may download the application program from the software deployment server 249 when the automatic driving-related program 247 needs to be executed.
  • the senor 253 may be associated with the computer system 201, and the sensor 253 may be used to detect the environment around the computer 201.
  • the sensor 253 can detect animals, cars, obstacles, and crosswalks. Further, the sensor can also detect the surrounding environment of the above-mentioned animals, cars, obstacles, and crosswalks, such as: the environment around the animals, for example, when the animals appear around them. Other animals, weather conditions, the brightness of the surrounding environment, etc.
  • the senor may be a camera, an infrared sensor, a chemical detector, a microphone, etc.
  • Multiple sensors can be used to detect the location of obstacles around the vehicle, and the location of the obstacles can be obtained based on the data obtained by the multiple sensors.
  • obtaining the location of the obstacle based on the data acquired by multiple sensors can be implemented by the target tracking method in the embodiment of the present application.
  • the computer system 150 shown in FIG. 2 may also receive information from other computer systems or transfer information to other computer systems.
  • the sensor data collected from the sensor system 120 of the vehicle 100 may be transferred to another computer to process the data.
  • data from the computer system 312 may be transmitted to the server 320 on the cloud side via the network for further processing.
  • the network and intermediate nodes can include various configurations and protocols, including the Internet, World Wide Web, Intranet, virtual private network, wide area network, local area network, private network using one or more company’s proprietary communication protocols, Ethernet, WiFi and HTTP, And various combinations of the foregoing; this communication can be by any device capable of transferring data to and from other computers, such as modems and wireless interfaces.
  • the server 320 may include a server with multiple computers, such as a load balancing server group, which exchanges information with different nodes of the network for the purpose of receiving, processing, and transmitting data from the computer system 312.
  • the server may be configured similarly to the computer system 312, with a processor 330, a memory 340, instructions 350, and data 360.
  • the data 360 of the server 320 may include information related to road conditions around the vehicle.
  • the server 320 may receive, detect, store, update, and transmit information related to the road conditions of the vehicle.
  • State space refers to a collection of state variables that describe all possible states of a system.
  • the car can be regarded as a system, and the operation of the car by the user can be regarded as input variables. When there is an input of an operation signal, it will have a clear impact on the speed, acceleration, angular velocity and other variables of the car. These affected variables can all be regarded as the system. The component of the state variable.
  • Observation refers to the process of obtaining state variable estimates directly or indirectly through a certain measurement method.
  • Bayesian estimation means that starting from any time k-1, the probability distribution of the state a priori estimate at the next time k is calculated, which is called prediction; then after obtaining the observation value at time k, the prediction link is obtained The a priori estimate of is revised, and the posterior estimate of the state at time k is obtained, which is called update.
  • Bayesian recursive estimation is called Bayesian filtering when used in actual engineering.
  • the Bayesian filter is divided into two parts: prediction and update.
  • X is the state quantity, that is, the output quantity of the Bayesian filter, and the space where the state quantity is located is the state space.
  • z is the observation, that is, the input of the Bayesian filter, and the space where the observation is located is the observation space.
  • f(*) is the probability density function.
  • Fig. 5 is a schematic diagram of a target tracking device according to an embodiment of the present application.
  • the target tracking device 402 can be applied to the computer system 401.
  • the target tracking device includes a prediction unit 410 and an update unit 420.
  • the prediction unit 410 includes a judgment module 411, and the update unit 420 includes an observation transfer module 422.
  • the prediction unit 410 may also include a state transition module 412.
  • the update unit 420 may also include an update module 421.
  • the judging module 411 is used to judge whether the collection time of the observation is earlier than the update time of the current state quantity. That is to judge whether the observation is late.
  • the current state quantity refers to the state quantity of the latest update. Specifically, when the collection time of the observation is earlier than the latest update time of the state quantity, the observation is a late observation.
  • the observation collected by the sensor is the actual observation.
  • the state transition module 412 is configured to obtain the predicted value of the state amount of the current update according to the state amount of the previous update when the collection time of the observation is not earlier than the update time of the current state amount.
  • the specific process satisfies formula (1).
  • the update module 421 is configured to update the state quantity according to the predicted value of the current state quantity obtained by the state transition module 412.
  • the specific process satisfies formula (2) and formula (3).
  • the observation transfer unit 422 is configured to perform state transition on the observation when the acquisition time of the observation is earlier than the update time of the current state quantity to obtain the estimated observation.
  • the observation transfer unit 422 is further configured to update the current state quantity according to the actual observation and the estimated observation, and use the updated state quantity as the state quantity corresponding to the update time.
  • FIG. 6 is a schematic flowchart of a target tracking method 500 according to an embodiment of the present application.
  • the method shown in FIG. 6 may be executed by the target tracking device in the embodiment of the present application.
  • the method 500 includes steps S510 to S560. Steps S510 to S560 will be described in detail below.
  • S510 Obtain actual observations of the moving object collected by the sensor.
  • the observations collected by the sensors are the actual observations, and the actual observations may also be referred to as "observations" in the embodiments of the present application.
  • Observation can include the speed of the moving object, the acceleration of the moving object, or the position of the moving object.
  • the observation may include the position and/or speed of the moving object measured by the millimeter wave radar.
  • S520 Determine whether the collection time of the observation is earlier than the update time of the current state quantity of the moving object.
  • step S530 If the collection time of the observation is not earlier than the update time of the current state quantity of the moving object, step S530 is executed. If the collection time of the observation is earlier than the update time of the current state quantity of the moving object, step S540 is executed.
  • the acquisition time of the observation is earlier than the update time of the current state quantity of the moving object. It can be understood that the observation can be used to update the current state quantity of the moving object, but it cannot actually be used for the current state quantity of the moving object. Update.
  • the collection time of the observation can be indicated by the time stamp of the observation.
  • the state quantity may include the speed of the moving object, the acceleration of the moving object, or the position of the moving object.
  • the state quantity may be the result of the tracking, and the result of the tracking may be the position of the moving object.
  • the update time of the current state quantity of the moving object can be understood as the time when the state quantity of the moving object was updated last time.
  • the state quantity of the moving object and the update time of the state quantity can be stored in the tracking list in the Bayesian filter.
  • judging whether the collection time of the observation is earlier than the update time of the current state quantity of the moving object may be judging whether the time stamp of the observation is earlier than the latest update time in the tracking list.
  • the current state quantity of the moving object can be updated through the Bayesian filter.
  • the current state quantity of the moving object can be updated according to the above formula (1), formula (2) and formula (3).
  • S540 Determine whether the time difference between the collection time of the observation and the update time of the current state quantity of the moving object is greater than a first threshold.
  • the observation is discarded. That is, the input observation is not used to update the state quantity. In this way, when the time difference is too large, using the observation to update the state quantity may reduce the confidence of the updated state quantity.
  • step S550 is executed.
  • step S540 may also be to determine whether the time difference between the collection time of the observation and the update time of the current state quantity of the moving object is greater than or equal to the first threshold.
  • step S550 is executed.
  • step S520, step S530, and step S540 are optional steps, and the method 500 of the embodiment of the present application may execute step S550 after step S510.
  • the state transition model may be a kinematics model of a moving object.
  • the kinematics model may include a uniform linear motion model and the like.
  • the likelihood function of the observation corresponding to the acquisition time is obtained according to the state transition model corresponding to the moving object.
  • the estimated observation is determined based on the likelihood function of the observation corresponding to the acquisition time.
  • the state transition model is based on the first The probability density function is determined, and the first probability density function is a probability density function that transfers the current state quantity to the state quantity corresponding to the collection time.
  • step S550 will be described below in conjunction with a formula.
  • the observation z k collected at the time t k is input to the Bayesian filter, and the Bayesian filter updates the state quantity according to the observation z k , and uses the obtained state quantity X k 'as the current state quantity X k '.
  • the latest update of the state quantity is time t k , that is, the update time of the current state quantity X k ′ is time t k .
  • t k-1 time View captured measurement z k-1 input Bayesian filter.
  • the acquisition time of the above observation is earlier than the update time of the current state quantity. It can be understood that the Bayesian filter receives the observation z k and updates the state quantity and then receives the observation z k-1 .
  • the state transition model from time t k to time t k-1 is f k-1
  • this state transition is Markov transition.
  • X k ') is the above-mentioned first probability density function, and the first probability density function is to transfer the state quantity X k ' updated at t k to the state quantity updated at t k-1 The probability density function.
  • X k ') represents the likelihood function of the observation z k-1 corresponding to the time t k-1
  • X represents the state quantity updated at the time t k-1
  • X) Z represents Z k-1 k-1 is used on the time t where k-1 is updated state quantity Likelihood function.
  • S560 Update the current state quantity according to the actual observation and the estimated observation, and use the updated state quantity as the state quantity corresponding to the update time.
  • the probability density function of the current state quantity may be updated according to the likelihood function of the observation corresponding to the acquisition time.
  • the state quantity corresponding to the update time is determined according to the probability density function of the updated state quantity.
  • step S560 will be described below in conjunction with a formula.
  • z k ) represents the probability density function of the state quantity X k updated at time t k
  • z k ,z k-2 ) represents the absence of z
  • k-1 the probability density function of the state quantity X k 'updated at time t k. That is, the posterior probability corresponding to the current state quantity X k'.
  • z k represents the observations collected at time t k
  • z k represents the set of observations collected at k time ⁇ z 1 ,z 2 ,...,z k-2 ,z k-1 ,z k ⁇
  • z k -2 represents the set of observations collected at k-2 moments ⁇ z 1 , z 2 ,..., z k-2 ⁇ .
  • the Kalman filter is an implementation form of the Bayes filter. Taking the Kalman filter as an example, step S560 will be described.
  • the above step S560 may include determining the expectation of the state quantity corresponding to the update time according to the actual observation and the estimated observation; and taking the expectation of the state quantity corresponding to the update time as the state quantity corresponding to the update time.
  • the expectation of the state quantity corresponding to the update time is related to the Kalman gain value.
  • the Kalman gain value is related to the first covariance, the observation matrix at the time of collection, the covariance of the observation matrix, and the variance of the observation at the time of collection.
  • the first covariance refers to the covariance of the state quantity transferred from the update time to the collection time.
  • the Kalman gain value satisfies:
  • k represents t K transferred from the update time to the acquisition time t covariance k-1 state quantity
  • k represents the lack Concept t k-1 time captured measurement Z k-1
  • the covariance of the state quantity updated at time t k F k-1
  • Q k represents the covariance of the prediction matrix
  • H k-1 Represents the observation matrix at time t k-1
  • Var(z k-1 ) represents the variance of the observation z k-1 collected at time t k-1
  • R k-1 represents the covariance of the observation matrix.
  • k represents the expectation of the state quantity updated at time t k
  • k corresponding to the update time is taken as the state quantity corresponding to the update time.
  • the estimated value of the observation corresponding to the collection time satisfies:
  • H k-1 represents a t the observation matrix k-1 time
  • k indicates a transition from t K time to t k-1 time of the state transition matrix
  • k t k represents a time measurement in the case where z k-1 is updated in a desired state quantity of the absence of the concept of time t k-1 collected.
  • the estimated observation at time t k-1 can be understood as the state transition matrix F k-1
  • the covariance of the updated state quantity at time t k is calculated.
  • the covariance of the updated state at t k satisfies:
  • k represents the covariance of the state quantity updated at time t k.
  • the state transition is performed according to the state transition model corresponding to the moving object, and then the current state quantity is updated, and the more accurate state quantity corresponding to the latest update time can be obtained, instead of the state quantity corresponding to the collection time. , To avoid out-of-order jumps of the output state quantity in time.
  • step S550 and step S560 will be described.
  • R represents the distance between the sensor and the target.
  • k t k represents a time measurement in the case where z k-1 is updated in a desired state quantity of the absence of the concept of time t k-1 collected.
  • k represents the missing covariance Concept time t k-1 measured state quantity acquisition case where z k-1 t k is the time updates.
  • k can be 3
  • k may be 1
  • t k-1 time View captured measurement z k-1 may be four.
  • k satisfies:
  • the Kalman gain value satisfies:
  • k can be used as the output result, that is, the state quantity corresponding to the time t k.
  • the target tracking method of the embodiment of the present application is described in detail above with reference to FIG. 6, and the device embodiment of the present application will be described in detail below with reference to FIG. 7 to FIG. 8. It should be understood that the target tracking device in the embodiment of the present application can execute the target tracking method of the foregoing embodiment of the present application, that is, the specific working process of the following various products can refer to the corresponding process in the foregoing method embodiment.
  • Fig. 7 is a schematic block diagram of a target tracking device provided by an embodiment of the present application. It should be understood that the target tracking device 1000 can execute the target tracking method shown in FIG. 6.
  • the target tracking device 1000 includes: an acquiring unit 1010 and a processing unit 1020.
  • the acquiring unit 1010 is used to acquire the actual observation of the moving object collected by the sensor.
  • the processing unit 1020 is used to obtain the estimated observation according to the state transition model corresponding to the moving object when the acquisition time of the actual observation is earlier than the update time of the current state quantity of the moving object;
  • the state quantity is updated, and the updated state quantity is used as the state quantity corresponding to the update time.
  • the processing unit 1020 is configured to: obtain the likelihood function of the observation corresponding to the acquisition time according to the state transition model corresponding to the moving object, and the estimated observation is determined according to the likelihood function of the observation corresponding to the acquisition time,
  • the state transition model is determined according to the first probability density function, and the first probability density function is the probability density function that transfers the current state quantity to the state quantity corresponding to the collection time.
  • the processing unit 1020 is configured to: update the probability density function of the current state quantity according to the likelihood function of the observation corresponding to the acquisition time; determine the state quantity corresponding to the update time according to the probability density function of the updated state quantity .
  • the likelihood function of the observation corresponding to the acquisition moment satisfies:
  • X k ') represents the likelihood function of the observation z k-1 corresponding to the time t k-1
  • X represents the state quantity updated at the time t k-1
  • X) Z represents Z k-1 k-1 is used on the time t where k-1 is updated state quantity
  • X k ') is the first probability density function, the first probability density function is to transfer the state quantity X k ' updated at t k to update at t k-1 The probability density function of the state quantity.
  • the probability density function of the updated state quantity satisfies:
  • z k ) represents the probability density function of the state quantity X k updated at time t k
  • z k ,z k-2 ) represents the absence of z
  • z k represents the observations collected at time t k
  • z k represents the set of observations collected at k time ⁇ z 1 , z 2 ,...,z k-2 ,z k-1 ,z k ⁇
  • z k-2 represents the set of observations collected at k-2 moments ⁇ z 1 ,z 2 ,...,z k-2 ⁇ .
  • the processing unit 1020 is configured to: determine the expectation of the state quantity corresponding to the update time according to actual observations and estimated observations; use the expectation of the state quantity corresponding to the update time as the state quantity corresponding to the update time; wherein, the update time corresponds to The expectation of the state quantity is related to the Kalman gain value.
  • the Kalman gain value is related to the first covariance, the observation matrix at the acquisition time, the covariance of the observation matrix, and the variance of the observation at the acquisition time.
  • the first covariance refers to It is the covariance of the state amount from the update time to the collection time.
  • the Kalman gain value satisfies:
  • k represents t K transferred from the update time to the acquisition time t covariance k-1 state quantity
  • k represents the lack Concept t k-1 time captured measurement Z k-1
  • the covariance of the state quantity updated at time t k F k-1
  • Q k represents the covariance of the prediction matrix
  • H k-1 Represents the observation matrix at time t k-1
  • Var(z k-1 ) represents the variance of the observation z k-1 collected at time t k-1
  • R k-1 represents the covariance of the observation matrix.
  • the expectation of the state quantity corresponding to the update time meets:
  • k represents the expectation of the state quantity updated at time t k
  • the estimated value of the observation corresponding to the collection time satisfies:
  • H k-1 represents a t the observation matrix k-1 time
  • k indicates a transition from t K time to t k-1 time of the state transition matrix
  • k represents the concept of a lack of time t k-1 measurement state quantity acquired in the case where time t k z k-1 update is desired.
  • the processing unit 1020 is configured to: in the case where the acquisition time of the observation is earlier than the update time of the current state quantity of the moving object, and the time difference between the acquisition time and the update time is less than or equal to a threshold, corresponding to the moving object
  • the state transition model is estimated to be observed.
  • target tracking device 1000 is embodied in the form of a functional unit.
  • unit herein can be implemented in the form of software and/or hardware, which is not specifically limited.
  • a "unit” may be a software program, a hardware circuit, or a combination of the two that realizes the above-mentioned functions.
  • the hardware circuit may include an application specific integrated circuit (ASIC), an electronic circuit, and a processor for executing one or more software or firmware programs (such as a shared processor, a dedicated processor, or a group processor). Etc.) and memory, merged logic circuits and/or other suitable components that support the described functions.
  • the units of the examples described in the embodiments of the present application can be implemented by electronic hardware or a combination of computer software and electronic hardware. Whether these functions are executed by hardware or software depends on the specific application and design constraint conditions of the technical solution. Professionals and technicians can use different methods for each specific application to implement the described functions, but such implementation should not be considered beyond the scope of this application.
  • FIG. 8 is a schematic diagram of the hardware structure of the target tracking device provided by an embodiment of the present application.
  • the target tracking apparatus 1200 (the target tracking 1200 may specifically be a computer device) includes a memory 1201, a processor 1202, a communication interface 1203, and a bus 1204. Among them, the memory 1201, the processor 1202, and the communication interface 1203 implement communication connections between each other through the bus 1204.
  • the memory 1201 may be a read only memory (ROM), a static storage device, a dynamic storage device, or a random access memory (RAM).
  • the memory 1201 may store a program.
  • the processor 1202 is configured to execute each step of the target tracking method of the embodiment of the present application, for example, execute each step shown in FIG. 6.
  • the target tracking device shown in the embodiment of the present application may be a server, for example, it may be a server in the cloud, or may also be a chip configured in a server in the cloud.
  • the processor 1202 may adopt a general central processing unit (CPU), a microprocessor, an application specific integrated circuit (ASIC), or one or more integrated circuits for executing related programs to realize this The target tracking method of the application method embodiment.
  • CPU central processing unit
  • ASIC application specific integrated circuit
  • the processor 1202 may also be an integrated circuit chip with signal processing capability.
  • each step of the target tracking method of the present application can be completed by an integrated logic circuit of hardware in the processor 1202 or instructions in the form of software.
  • the above-mentioned processor 1202 may also be a general-purpose processor, a digital signal processing (digital signal processing, DSP), an application-specific integrated circuit (ASIC), an off-the-shelf programmable gate array (field programmable gate array, FPGA) or other programmable logic devices, Discrete gates or transistor logic devices, discrete hardware components.
  • DSP digital signal processing
  • ASIC application-specific integrated circuit
  • FPGA field programmable gate array
  • the general-purpose processor may be a microprocessor or the processor may also be any conventional processor or the like.
  • the steps of the method disclosed in the embodiments of the present application may be directly embodied as being executed and completed by a hardware decoding processor, or executed and completed by a combination of hardware and software modules in the decoding processor.
  • the software module can be located in a mature storage medium in the field, such as random access memory, flash memory, read-only memory, programmable read-only memory, or electrically erasable programmable memory, registers.
  • the storage medium is located in the memory 1201, and the processor 1202 reads the information in the memory 1201, and combines its hardware to complete the functions required by the units included in the target tracking device shown in FIG. 7 in the implementation of this application, or execute the method of this application The target tracking method shown in FIG. 6 of the embodiment.
  • the communication interface 1203 uses a transceiving device such as but not limited to a transceiver to implement communication between the target tracking device 1200 and other devices or a communication network.
  • a transceiving device such as but not limited to a transceiver to implement communication between the target tracking device 1200 and other devices or a communication network.
  • the bus 1204 may include a path for transferring information between various components of the target tracking device 1200 (for example, the memory 1201, the processor 1202, and the communication interface 1203).
  • target tracking device 1200 only shows a memory, a processor, and a communication interface, in the specific implementation process, those skilled in the art should understand that the target tracking device 1200 may also include other necessary for normal operation. Device. At the same time, according to specific needs, those skilled in the art should understand that the above-mentioned target tracking device 1200 may also include hardware devices that implement other additional functions.
  • target tracking device 1200 may also only include the necessary components to implement the embodiments of the present application, and not necessarily include all the components shown in FIG. 8.
  • the disclosed system, device, and method can be implemented in other ways.
  • the device embodiments described above are only illustrative.
  • the division of the units is only a logical function division, and there may be other divisions in actual implementation, for example, multiple units or components may be combined or It can be integrated into another system, or some features can be ignored or not implemented.
  • the displayed or discussed mutual coupling or direct coupling or communication connection may be indirect coupling or communication connection through some interfaces, devices or units, and may be in electrical, mechanical or other forms.
  • the units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in one place, or they may be distributed on multiple network units. Some or all of the units may be selected according to actual needs to achieve the objectives of the solutions of the embodiments.
  • the functional units in the various embodiments of the present application may be integrated into one processing unit, or each unit may exist alone physically, or two or more units may be integrated into one unit.
  • the function is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer readable storage medium.
  • the technical solution of the present application essentially or the part that contributes to the existing technology or the part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including Several instructions are used to make a computer device (which may be a personal computer, a server, or an access network device, etc.) execute all or part of the steps of the methods described in the various embodiments of the present application.
  • the aforementioned storage media include: U disk, mobile hard disk, read-only memory (read-only Memory, ROM), random access memory (random access memory, RAM), magnetic disk or optical disk and other media that can store program code .

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Traffic Control Systems (AREA)

Abstract

A target tracking method and a target tracking apparatus in the field of artificial intelligence. The method comprises: an actual observation quantity of a moving body captured by a sensor is obtained; if a capture time of the actual observation quantity is earlier than an update time of a current state quantity of the moving body, an estimated observation quantity is obtained according to a state transition model for the moving body; updating is performed on the current state quantity according to the actual observation quantity and the estimated observation quantity, and the updated state quantity serves as a state quantity corresponding to the update time. The present invention is able to implement target tracking when capture times of observation quantities are in the wrong order, preventing state quantity time stamps from jumping around in an incorrect order due to observation quantity time stamps becoming out of order.

Description

目标跟踪方法以及目标跟踪装置Target tracking method and target tracking device 技术领域Technical field
本申请实施例涉及人工智能领域,并且更具体地,涉及一种目标跟踪方法以及目标跟踪装置。The embodiments of the present application relate to the field of artificial intelligence, and more specifically, to a target tracking method and a target tracking device.
背景技术Background technique
目标跟踪方法能够应用于很多场景,如目标检测、车道线检测、定位等。常见的目标跟踪方法,需要利用传感器(sensor)的测量数据,结合所跟踪的对象自身的状态转移模型进行贝叶斯估计(Bayesian estimation)。例如可以采用贝叶斯滤波器(Bayes filters)实现目标跟踪。具体地,将传感器的测量数据输入贝叶斯滤波器,得到对所跟踪的对象的状态的估计。然而,在存在多个传感器的情况下,由于传输路径等原因,可能存在一个传感器的测量数据到达贝叶斯滤波器的时刻晚于其本应到达贝叶斯滤波器的时刻。图1是一种输入贝叶斯滤波器的测量数据在时间上出现乱序的示意图。图1中横坐标示出了真实时间(real time),传感器传递的数据包括时间戳(header:Xms)和测量数据(data),时间戳用于表示该测量数据的测量或采集的时刻为Xms时刻。传感器0于1ms时刻测量的测量数据在3ms时刻输入贝叶斯滤波器,传感器1于0ms时刻测量的测量数据在4ms时刻输入贝叶斯滤波器,造成输入贝叶斯滤波器的测量数据的时间上的乱序。在该情况下,贝叶斯滤波器先输出基于1ms时刻测量的测量数据计算得到的时间戳为1ms时刻的状态量,然后输出基于0ms时刻测量的测量数据计算得到的时间戳为0ms时刻的状态量,导致输出状态量的时间戳随着测量数据的时间戳的变化产生跳动。The target tracking method can be applied to many scenarios, such as target detection, lane line detection, positioning and so on. Common target tracking methods need to use sensor measurement data combined with the state transition model of the tracked object to perform Bayesian estimation. For example, Bayes filters can be used to achieve target tracking. Specifically, the measurement data of the sensor is input to the Bayesian filter to obtain an estimate of the state of the tracked object. However, in the case where there are multiple sensors, due to transmission paths and other reasons, there may be a time when the measurement data of one sensor reaches the Bayes filter later than the time when it should reach the Bayes filter. Fig. 1 is a schematic diagram of the measurement data input to the Bayesian filter appearing out of order in time. The abscissa in Figure 1 shows the real time. The data transmitted by the sensor includes a timestamp (header: Xms) and measurement data (data). The timestamp is used to indicate that the measurement or collection time of the measurement data is Xms. time. The measurement data measured by sensor 0 at 1ms is input into the Bayesian filter at 3ms, and the measurement data measured by sensor 1 at 0ms is input into the Bayesian filter at 4ms, resulting in the time to input the measurement data of the Bayesian filter. Out of order. In this case, the Bayesian filter first outputs the state quantity calculated based on the measurement data measured at 1ms with a time stamp of 1ms, and then outputs the state calculated based on the measurement data measured at 0ms with a time stamp of 0ms As a result, the time stamp of the output state quantity jumps with the change of the time stamp of the measurement data.
通常可以采用两种方法解决上述问题。一种方法为判断测量数据的时间戳是否晚于贝叶斯滤波器中最新一次更新的状态量的时间戳,仅将采集时刻晚于最近一次更新的状态量的更新时刻的测量数据输入至贝叶斯滤波器。但该方法会直接造成观测量的缺失,导致输出结果的置信度下降。另一种方法为将迟到的观测数据外推至当前时刻。然而,在观测空间中做外推依赖于强假设,容易扩大误差。There are usually two ways to solve the above problems. One method is to determine whether the time stamp of the measurement data is later than the time stamp of the latest updated state quantity in the Bayesian filter, and only the measurement data whose collection time is later than the update time of the latest state quantity is input to the Bayesian filter. Yes filter. However, this method will directly cause the lack of observations, leading to a decrease in the confidence of the output results. Another method is to extrapolate the late observation data to the current moment. However, extrapolation in the observation space relies on strong assumptions and tends to enlarge errors.
发明内容Summary of the invention
本申请提供一种目标跟踪方法以及目标跟踪装置,能够在观测量的采集时刻出现乱序的情况下实现目标跟踪,避免状态量的时间戳随观测量的时间戳的乱序变化发生乱序跳动。The present application provides a target tracking method and a target tracking device, which can realize target tracking when the acquisition time of the observation is out of order, and avoid the out-of-order jump of the time stamp of the state quantity with the out-of-order change of the time stamp of the observation. .
第一方面,提供了一种目标跟踪方法,包括:获取传感器采集的运动物体的实际观测量;在实际观测量的采集时刻早于运动物体的当前状态量的更新时刻的情况下,根据运动物体对应的状态转移模型得到估计观测量;根据实际观测量和估计观测量对当前状态量进行更新,将更新后的状态量作为该更新时刻对应的状态量。In the first aspect, a target tracking method is provided, which includes: acquiring the actual observation of the moving object collected by the sensor; in the case that the acquisition time of the actual observation is earlier than the update time of the current state of the moving object, according to the moving object The corresponding state transition model obtains the estimated observation; the current state quantity is updated according to the actual observation and the estimated observation, and the updated state quantity is used as the state quantity corresponding to the update moment.
观测量可以包括运动物体的速度、运动物体的加速度或运动物体的位置等。例如,该观测量可以包括毫米波雷达测量的运动物体的位置和/或速度等。Observation can include the speed of the moving object, the acceleration of the moving object, or the position of the moving object. For example, the observation may include the position and/or speed of the moving object measured by the millimeter wave radar.
状态量可以包括运动物体的速度、运动物体的加速度或运动物体的位置等。The state quantity may include the speed of the moving object, the acceleration of the moving object, or the position of the moving object.
运动物体的当前状态量的更新时刻,可以理解为,运动物体的状态量最新一次更新的时刻。The update time of the current state quantity of the moving object can be understood as the time when the state quantity of the moving object was updated last time.
可选地,运动物体的状态转移模型包括运动物体的运动学模型。例如,该运动学模型可以包括匀速直线运动模型等。Optionally, the state transition model of the moving object includes a kinematics model of the moving object. For example, the kinematics model may include a uniform linear motion model and the like.
根据本申请实施例的方案,根据运动物体对应的状态转移模型进行状态转移,进而对当前状态量进行更新,能够得到较准确的最新的更新时刻对应的状态量。同时,更新后的状态量为最新的更新时刻对应的状态量,而不是采集时刻对应的状态量,避免输出的状态量在时间上出现乱序跳动。According to the solution of the embodiment of the present application, the state transition is performed according to the state transition model corresponding to the moving object, and then the current state quantity is updated, so that a more accurate state quantity corresponding to the latest update time can be obtained. At the same time, the updated state quantity is the state quantity corresponding to the latest update time, rather than the state quantity corresponding to the collection time, so as to avoid the out-of-order jump of the output state quantity in time.
结合第一方面,在第一方面的某些实现方式中,根据运动物体对应的状态转移模型得到估计观测量,包括:根据运动物体对应的状态转移模型得到关于采集时刻对应的观测量的似然函数,估计观测量是根据关于采集时刻对应的观测量的似然函数确定的,状态转移模型是根据第一概率密度函数确定的,第一概率密度函数是将当前状态量转移至采集时刻对应的状态量的概率密度函数。In combination with the first aspect, in some implementations of the first aspect, obtaining estimated observations according to the state transition model corresponding to the moving object includes: obtaining the likelihood of the observation corresponding to the acquisition time according to the state transition model corresponding to the moving object Function, the estimated observation is determined according to the likelihood function of the observation corresponding to the acquisition time, the state transition model is determined according to the first probability density function, the first probability density function is to transfer the current state quantity to the corresponding acquisition time The probability density function of the state quantity.
结合第一方面,在第一方面的某些实现方式中,根据实际观测量和估计观测量对当前状态量进行更新,将更新后的状态量作为更新时刻对应的状态量,包括:根据关于采集时刻对应的观测量的似然函数对当前状态量的概率密度函数进行更新;根据更新后的状态量的概率密度函数确定更新时刻对应的状态量。In combination with the first aspect, in some implementations of the first aspect, the current state quantity is updated based on actual observations and estimated observations, and the updated state quantity is used as the state quantity corresponding to the update time, including: The likelihood function of the observation corresponding to the time updates the probability density function of the current state quantity; the state quantity corresponding to the update time is determined according to the probability density function of the updated state quantity.
结合第一方面,在第一方面的某些实现方式中,关于采集时刻对应的观测量的似然函数满足:In combination with the first aspect, in some implementations of the first aspect, the likelihood function of the observation corresponding to the acquisition time satisfies:
g(z k-1|X k')=∫g(z k-1|X)f k-1|k(X|X k')dX g(z k-1 |X k ')=∫g(z k-1 |X)f k-1|k (X|X k ')dX
其中,g(z k-1|X k')表示关于t k-1时刻对应的观测量z k-1的似然函数,X表示t k-1时刻更新的状态量,X k'表示在缺少z k-1的情况下t k时刻更新的状态量,g(z k-1|X)表示在z k-1用于t k-1时刻更新状态量的情况下关于z k-1的似然函数,f k-1|k(X|X k')为第一概率密度函数,第一概率密度函数是将t k时刻更新的状态量X k'转移至到t k-1时刻更新的状态量的概率密度函数。 Among them, g(z k-1 |X k ') represents the likelihood function of the observation z k-1 corresponding to the time t k-1 , X represents the state quantity updated at the time t k-1 , and X k 'represents the missing state quantity update time t K in the case of Z k-1, g (z k-1 | X) Z represents Z k-1 k-1 is used on the time t where k-1 is updated state quantity Likelihood function, f k-1|k (X|X k ') is the first probability density function, the first probability density function is to transfer the state quantity X k ' updated at t k to update at t k-1 The probability density function of the state quantity.
结合第一方面,在第一方面的某些实现方式中,更新后的状态量的概率密度函数满足:In combination with the first aspect, in some implementations of the first aspect, the probability density function of the updated state quantity satisfies:
Figure PCTCN2020075556-appb-000001
Figure PCTCN2020075556-appb-000001
其中,f k|k(X k|z k)表示t k时刻更新的状态量X k的概率密度函数,f' k|k(X k'|z k,z k-2)表示在缺少z k-1的情况下t k时刻更新的状态量X k'的概率密度函数,z k表示t k时刻采集的观测量,z k表示k个时刻采集的观测量的集合{z 1,z 2,...,z k-2,z k-1,z k},z k-2表示k-2个时刻采集的观测量的集合{z 1,z 2,...,z k-2}。 Among them, f k|k (X k |z k ) represents the probability density function of the state quantity X k updated at time t k , and f'k|k (X k '|z k ,z k-2 ) represents the absence of z In the case of k-1 , the probability density function of the state quantity X k 'updated at time t k , z k represents the observations collected at time t k , and z k represents the set of observations collected at k time {z 1 , z 2 ,...,z k-2 ,z k-1 ,z k }, z k-2 represents the set of observations collected at k-2 moments {z 1 ,z 2 ,...,z k-2 }.
结合第一方面,在第一方面的某些实现方式中,根据实际观测量和估计观测量对当前状态量进行更新,将更新后的状态量作为更新时刻对应的状态量,包括:In combination with the first aspect, in some implementations of the first aspect, the current state quantity is updated based on the actual observations and estimated observations, and the updated state quantity is used as the state quantity corresponding to the update time, including:
根据实际观测量和估计观测量确定更新时刻对应的状态量的期望;Determine the expectation of the state quantity corresponding to the update time according to the actual observations and estimated observations;
将更新时刻对应的状态量的期望作为更新时刻对应的状态量;The expectation of the state quantity corresponding to the update time is regarded as the state quantity corresponding to the update time;
其中,更新时刻对应的状态量的期望与由卡尔曼增益值有关,卡尔曼增益值与第一协 方差、采集时刻的观测矩阵、观测矩阵的协方差和采集时刻的观测量的方差有关,第一协方差指的是由更新时刻转移至采集时刻的状态量的协方差。Among them, the expectation of the state quantity corresponding to the update time is related to the Kalman gain value. The Kalman gain value is related to the first covariance, the observation matrix at the acquisition time, the covariance of the observation matrix, and the variance of the observation at the acquisition time. A covariance refers to the covariance of the state quantity transferred from the update time to the collection time.
结合第一方面,在第一方面的某些实现方式中,卡尔曼增益值满足:In combination with the first aspect, in some implementations of the first aspect, the Kalman gain value satisfies:
Figure PCTCN2020075556-appb-000002
Figure PCTCN2020075556-appb-000002
Var(z k-1)满足: Var(z k-1 ) satisfies:
Figure PCTCN2020075556-appb-000003
Figure PCTCN2020075556-appb-000003
P k-1|k满足: P k-1|k satisfies:
Figure PCTCN2020075556-appb-000004
Figure PCTCN2020075556-appb-000004
其中,P k-1|k表示由更新时刻t k转移到采集时刻t k-1的状态量的协方差,P' k|k表示在缺少t k-1时刻采集的观测量z k-1的情况下t k时刻更新的状态量的协方差,F k-1|k表示从t k时刻转移到t k-1时刻的状态转移矩阵,Q k表示预测矩阵的协方差,H k-1表示t k-1时刻的观测矩阵,Var(z k-1)表示t k-1时刻采集的观测量z k-1的方差,R k-1表示观测矩阵的协方差。 Wherein, P k-1 | k represents t K transferred from the update time to the acquisition time t covariance k-1 state quantity, P 'k | k represents the lack Concept t k-1 time captured measurement Z k-1 In the case of t k , the covariance of the state quantity updated at time t k, F k-1|k represents the state transition matrix from t k to t k-1 , Q k represents the covariance of the prediction matrix, H k-1 Represents the observation matrix at time t k-1 , Var(z k-1 ) represents the variance of the observation z k-1 collected at time t k-1 , and R k-1 represents the covariance of the observation matrix.
结合第一方面,在第一方面的某些实现方式中,更新时刻对应的状态量的期望满足:In combination with the first aspect, in some implementations of the first aspect, the expectation of the state quantity corresponding to the update time meets:
Figure PCTCN2020075556-appb-000005
Figure PCTCN2020075556-appb-000005
其中,x k|k表示t k时刻更新的状态量的期望,
Figure PCTCN2020075556-appb-000006
表示t k-1时刻的估计观测量。
Among them, x k|k represents the expectation of the state quantity updated at time t k,
Figure PCTCN2020075556-appb-000006
Represents the estimated observation at time t k-1.
结合第一方面,在第一方面的某些实现方式中,采集时刻对应的观测量的估计值满足:In combination with the first aspect, in some implementations of the first aspect, the estimated value of the observation corresponding to the collection time satisfies:
Figure PCTCN2020075556-appb-000007
Figure PCTCN2020075556-appb-000007
其中,
Figure PCTCN2020075556-appb-000008
表示t k-1时刻的估计观测量,H k-1表示t k-1时刻的观测矩阵,F k-1|k表示从t k时刻转移到t k-1时刻的状态转移矩阵,x' k|k表示在缺少t k-1时刻采集的观测量z k-1的情况下的t k时刻更新的状态量的期望。
in,
Figure PCTCN2020075556-appb-000008
Represents t estimated observations k-1 time point, H k-1 represents a t the observation matrix k-1 time, F k-1 | k indicates a transition from t K time to t k-1 time of the state transition matrix, x ' k | k represents the concept of a lack of time t k-1 measurement state quantity acquired in the case where time t k z k-1 update is desired.
结合第一方面,在第一方面的某些实现方式中,在实际观测量的采集时刻早于运动物体的当前状态量的更新时刻的情况下,根据运动物体对应的状态转移模型得到估计观测量,包括:在观测量的采集时刻早于运动物体的当前状态量的更新时刻,且采集时刻与更新时刻之间的时间差小于或等于阈值的情况下,根据运动物体对应的状态转移模型得到估计观测量。In combination with the first aspect, in some implementations of the first aspect, in the case where the acquisition time of the actual observation is earlier than the update time of the current state of the moving object, the estimated observation is obtained according to the state transition model corresponding to the moving object , Including: when the acquisition time of the observation is earlier than the update time of the current state quantity of the moving object, and the time difference between the acquisition time and the update time is less than or equal to the threshold value, the estimated observation is obtained according to the state transition model corresponding to the moving object quantity.
根据本申请实施例的方案,只有在采集时刻与更新时刻之间的时间差小于或等于阈值的情况下才进行状态转移,进而对状态量进行更新,这样可以避免在时间差过大的情况下,利用该观测量更新状态量而导致更新的状态量的置信度降低。According to the solution of the embodiment of the present application, the state transition is performed only when the time difference between the collection time and the update time is less than or equal to the threshold, and then the state quantity is updated, which can avoid the use of This observation updates the state quantity, which reduces the confidence of the updated state quantity.
第二方面,提供了一种目标跟踪装置,包括获取模块和处理模块,其中,获取模块用于:获取传感器采集的运动物体的实际观测量;处理模块用于:在实际观测量的采集时刻早于运动物体的当前状态量的更新时刻的情况下,根据运动物体对应的状态转移模型得到估计观测量;根据实际观测量和估计观测量对当前状态量进行更新,将更新后的状态量作为更新时刻对应的状态量。In a second aspect, a target tracking device is provided, which includes an acquisition module and a processing module, wherein the acquisition module is used to: acquire the actual observation of a moving object collected by the sensor; the processing module is used to: early in the acquisition time of the actual observation In the case of the update time of the current state quantity of the moving object, the estimated observation quantity is obtained according to the state transition model corresponding to the moving object; the current state quantity is updated according to the actual observation quantity and the estimated observation quantity, and the updated state quantity is used as the update The state quantity corresponding to the moment.
根据本申请实施例的方案,根据运动物体对应的状态转移模型进行状态转移,进而对当前状态量进行更新,能够得到较准确的最新的更新时刻对应的状态量。同时,更新后的状态量为最新的更新时刻对应的状态量,而不是采集时刻对应的状态量,避免输出的状态量在时间上出现乱序跳动。According to the solution of the embodiment of the present application, the state transition is performed according to the state transition model corresponding to the moving object, and then the current state quantity is updated, so that a more accurate state quantity corresponding to the latest update time can be obtained. At the same time, the updated state quantity is the state quantity corresponding to the latest update time, rather than the state quantity corresponding to the collection time, so as to avoid the out-of-order jump of the output state quantity in time.
结合第二方面,在第二方面的某些实现方式中,处理模块用于:根据运动物体对应的状态转移模型得到关于采集时刻对应的观测量的似然函数,估计观测量是根据关于采集时 刻对应的观测量的似然函数确定的,状态转移模型是根据第一概率密度函数确定的,第一概率密度函数是将当前状态量转移至采集时刻对应的状态量的概率密度函数。In conjunction with the second aspect, in some implementations of the second aspect, the processing module is used to obtain the likelihood function of the observation corresponding to the acquisition time according to the state transition model corresponding to the moving object, and the estimated observation is based on the acquisition time The likelihood function of the corresponding observation is determined, and the state transition model is determined according to the first probability density function. The first probability density function is the probability density function that transfers the current state quantity to the state quantity corresponding to the collection time.
结合第二方面,在第二方面的某些实现方式中,处理模块用于:根据关于采集时刻对应的观测量的似然函数对当前状态量的概率密度函数进行更新;根据更新后的状态量的概率密度函数确定更新时刻对应的状态量。With reference to the second aspect, in some implementations of the second aspect, the processing module is used to: update the probability density function of the current state quantity according to the likelihood function of the observation corresponding to the acquisition time; according to the updated state quantity The probability density function of determines the state quantity corresponding to the update moment.
结合第二方面,在第二方面的某些实现方式中,关于采集时刻对应的观测量的似然函数满足:In combination with the second aspect, in some implementations of the second aspect, the likelihood function of the observation corresponding to the acquisition time satisfies:
g(z k-1|X k')=∫g(z k-1|X)f k-1|k(X|X k')dX g(z k-1 |X k ')=∫g(z k-1 |X)f k-1|k (X|X k ')dX
其中,g(z k-1|X k')表示关于t k-1时刻对应的观测量z k-1的似然函数,X表示t k-1时刻更新的状态量,X k'表示在缺少z k-1的情况下t k时刻更新的状态量,g(z k-1|X)表示在z k-1用于t k-1时刻更新状态量的情况下关于z k-1的似然函数,f k-1|k(X|X k')为第一概率密度函数,第一概率密度函数是将t k时刻更新的状态量X k'转移至到t k-1时刻更新的状态量的概率密度函数。 Among them, g(z k-1 |X k ') represents the likelihood function of the observation z k-1 corresponding to the time t k-1 , X represents the state quantity updated at the time t k-1 , and X k 'represents the missing state quantity update time t K in the case of Z k-1, g (z k-1 | X) Z represents Z k-1 k-1 is used on the time t where k-1 is updated state quantity Likelihood function, f k-1|k (X|X k ') is the first probability density function, the first probability density function is to transfer the state quantity X k ' updated at t k to update at t k-1 The probability density function of the state quantity.
结合第二方面,在第二方面的某些实现方式中,更新后的状态量的概率密度函数满足:In combination with the second aspect, in some implementations of the second aspect, the probability density function of the updated state quantity satisfies:
Figure PCTCN2020075556-appb-000009
Figure PCTCN2020075556-appb-000009
其中,f k|k(X k|z k)表示t k时刻更新的状态量X k的概率密度函数,f' k|k(X k'|z k,z k-2)表示在缺少z k-1的情况下t k时刻更新的状态量X k'的概率密度函数,z k表示t k时刻采集的观测量,z k表示k个时刻采集的观测量的集合{z 1,z 2,...,z k-2,z k-1,z k},z k-2表示k-2个时刻采集的观测量的集合{z 1,z 2,...,z k-2}。 Among them, f k|k (X k |z k ) represents the probability density function of the state quantity X k updated at time t k , and f'k|k (X k '|z k ,z k-2 ) represents the absence of z In the case of k-1 , the probability density function of the state quantity X k 'updated at time t k , z k represents the observations collected at time t k , and z k represents the set of observations collected at k time {z 1 , z 2 ,...,z k-2 ,z k-1 ,z k }, z k-2 represents the set of observations collected at k-2 moments {z 1 ,z 2 ,...,z k-2 }.
结合第二方面,在第二方面的某些实现方式中,处理模块用于:根据实际观测量和估计观测量确定更新时刻对应的状态量的期望;将更新时刻对应的状态量的期望作为更新时刻对应的状态量;其中,更新时刻对应的状态量的期望与由卡尔曼增益值有关,卡尔曼增益值与第一协方差、采集时刻的观测矩阵、观测矩阵的协方差和采集时刻的观测量的方差有关,第一协方差指的是由更新时刻转移至采集时刻的状态量的协方差。In combination with the second aspect, in some implementations of the second aspect, the processing module is used to: determine the expectation of the state quantity corresponding to the update time according to the actual observation and the estimated observation; and use the expectation of the state quantity corresponding to the update time as the update The state quantity corresponding to the time; among them, the expectation of the state quantity corresponding to the update time is related to the Kalman gain value, the Kalman gain value and the first covariance, the observation matrix at the acquisition time, the covariance of the observation matrix and the observation at the acquisition time The variance of the quantity is related, and the first covariance refers to the covariance of the state quantity transferred from the update time to the collection time.
结合第二方面,在第二方面的某些实现方式中,卡尔曼增益值满足:In combination with the second aspect, in some implementations of the second aspect, the Kalman gain value satisfies:
Figure PCTCN2020075556-appb-000010
Figure PCTCN2020075556-appb-000010
Var(z k-1)满足: Var(z k-1 ) satisfies:
Figure PCTCN2020075556-appb-000011
Figure PCTCN2020075556-appb-000011
P k-1|k满足: P k-1|k satisfies:
Figure PCTCN2020075556-appb-000012
Figure PCTCN2020075556-appb-000012
其中,P k-1|k表示由更新时刻t k转移到采集时刻t k-1的状态量的协方差,P' k|k表示在缺少t k-1时刻采集的观测量z k-1的情况下t k时刻更新的状态量的协方差,F k-1|k表示从t k时刻转移到t k-1时刻的状态转移矩阵,Q k表示预测矩阵的协方差,H k-1表示t k-1时刻的观测矩阵,Var(z k-1)表示t k-1时刻采集的观测量z k-1的方差,R k-1表示观测矩阵的协方差。 Wherein, P k-1 | k represents t K transferred from the update time to the acquisition time t covariance k-1 state quantity, P 'k | k represents the lack Concept t k-1 time captured measurement Z k-1 In the case of t k , the covariance of the state quantity updated at time t k, F k-1|k represents the state transition matrix from t k to t k-1 , Q k represents the covariance of the prediction matrix, H k-1 Represents the observation matrix at time t k-1 , Var(z k-1 ) represents the variance of the observation z k-1 collected at time t k-1 , and R k-1 represents the covariance of the observation matrix.
结合第二方面,在第二方面的某些实现方式中,更新时刻对应的状态量的期望满足:In combination with the second aspect, in some implementations of the second aspect, the expectation of the state quantity corresponding to the update time meets:
Figure PCTCN2020075556-appb-000013
Figure PCTCN2020075556-appb-000013
其中,x k|k表示t k时刻更新的状态量的期望,
Figure PCTCN2020075556-appb-000014
表示t k-1时刻的估计观测量。
Among them, x k|k represents the expectation of the state quantity updated at time t k,
Figure PCTCN2020075556-appb-000014
Represents the estimated observation at time t k-1.
结合第二方面,在第二方面的某些实现方式中,采集时刻对应的观测量的估计值满足:In combination with the second aspect, in some implementations of the second aspect, the estimated value of the observation corresponding to the collection time satisfies:
Figure PCTCN2020075556-appb-000015
Figure PCTCN2020075556-appb-000015
其中,
Figure PCTCN2020075556-appb-000016
表示t k-1时刻的估计观测量,H k-1表示t k-1时刻的观测矩阵,F k-1|k表示从t k时刻转移到t k-1时刻的状态转移矩阵,x' k|k表示在缺少t k-1时刻采集的观测量z k-1的情况下的t k时刻更新的状态量的期望。
in,
Figure PCTCN2020075556-appb-000016
Represents t estimated observations k-1 time point, H k-1 represents a t the observation matrix k-1 time, F k-1 | k indicates a transition from t K time to t k-1 time of the state transition matrix, x ' k | k represents the concept of a lack of time t k-1 measurement state quantity acquired in the case where time t k z k-1 update is desired.
结合第二方面,在第二方面的某些实现方式中,处理模块用于:在实际观测量的采集时刻早于运动物体的当前状态量的更新时刻,且采集时刻与更新时刻之间的时间差小于或等于阈值的情况下,根据运动物体对应的状态转移模型得到估计观测量。With reference to the second aspect, in some implementations of the second aspect, the processing module is used to: the actual observation measurement acquisition time is earlier than the update time of the current state quantity of the moving object, and the time difference between the acquisition time and the update time When the value is less than or equal to the threshold, the estimated observation is obtained according to the state transition model corresponding to the moving object.
应理解,在上述第一方面中对相关内容的扩展、限定、解释和说明也适用于第二方面中相同的内容。It should be understood that the expansion, limitation, explanation and description of the related content in the above-mentioned first aspect are also applicable to the same content in the second aspect.
第三方面,提供了一种目标跟踪装置,该装置包括:存储器,用于存储程序;处理器,用于执行存储器存储的程序,当存储器存储的程序被执行时,处理器用于执行第一方面中的方法。In a third aspect, a target tracking device is provided. The device includes: a memory for storing a program; a processor for executing the program stored in the memory, and when the program stored in the memory is executed, the processor is configured to execute the first aspect In the method.
第四方面,提供一种计算机程序产品,所述计算机程序产品包括:计算机程序代码,当该计算机程序产品在计算机上运行时,使得计算机执行上述第一方面中的方法。In a fourth aspect, a computer program product is provided. The computer program product includes: computer program code, which when the computer program product runs on a computer, causes the computer to execute the method in the first aspect.
第五方面,提供了一种计算机可读存储介质,该计算机可读存储介质存储有计算机程序,当该计算机程序在计算机上运行时,使得计算机执行上述第一方面中的方法。In a fifth aspect, a computer-readable storage medium is provided, the computer-readable storage medium stores a computer program, and when the computer program runs on a computer, the computer executes the method in the first aspect.
应理解,本申请中,第一方面的方法具体可以是指第一方面以及第一方面中各种实现方式中的任意一种实现方式中的方法。It should be understood that, in this application, the method of the first aspect may specifically refer to the first aspect and a method in any one of the various implementation manners of the first aspect.
附图说明Description of the drawings
图1是测量数据输入贝叶斯滤波的示意图;Figure 1 is a schematic diagram of measurement data input into Bayesian filtering;
图2是本申请实施例提供的一种车辆的结构示意图;Fig. 2 is a schematic structural diagram of a vehicle provided by an embodiment of the present application;
图3是本申请实施例提供的一种计算机系统的结构示意图;Fig. 3 is a schematic structural diagram of a computer system provided by an embodiment of the present application;
图4是本申请实施例提供的一种云侧指令自动驾驶车辆的应用示意图;FIG. 4 is a schematic diagram of the application of a cloud-side command automatic driving vehicle provided by an embodiment of the present application;
图5是本申请实施例提供的一种目标跟踪装置的结构示意图;FIG. 5 is a schematic structural diagram of a target tracking device provided by an embodiment of the present application;
图6是本申请实施例提供的一种目标跟踪方法的示意性流程图;FIG. 6 is a schematic flowchart of a target tracking method provided by an embodiment of the present application;
图7是本申请实施例提供的另一种目标跟踪装置的结构示意图;FIG. 7 is a schematic structural diagram of another target tracking device provided by an embodiment of the present application;
图8是本申请实施例提供的又一种目标跟踪装置的结构示意图。FIG. 8 is a schematic structural diagram of another target tracking device provided by an embodiment of the present application.
具体实施方式Detailed ways
下面将结合附图,对本申请中的技术方案进行描述。The technical solution in this application will be described below in conjunction with the accompanying drawings.
图2是本申请实施例提供的车辆100的功能框图。Fig. 2 is a functional block diagram of a vehicle 100 provided by an embodiment of the present application.
其中,车辆100可以是人工驾驶车辆,或者可以将车辆100配置可以为完全或部分地自动驾驶模式。Among them, the vehicle 100 may be a manually driven vehicle, or the vehicle 100 may be configured in a fully or partially automatic driving mode.
在一个示例中,车辆100可以在处于自动驾驶模式中的同时控制自车,并且可通过人为操作来确定车辆及其周边环境的当前状态,确定周边环境中的至少一个其他车辆的可能行为,并确定其他车辆执行可能行为的可能性相对应的置信水平,基于所确定的信息来控制车辆100。在车辆100处于自动驾驶模式中时,可以将车辆100置为在没有和人交互的 情况下操作。In one example, the vehicle 100 can control its own vehicle while in the automatic driving mode, and can determine the current state of the vehicle and its surrounding environment through human operations, determine the possible behavior of at least one other vehicle in the surrounding environment, and The confidence level corresponding to the possibility of other vehicles performing possible behaviors is determined, and the vehicle 100 is controlled based on the determined information. When the vehicle 100 is in the automatic driving mode, the vehicle 100 can be placed to operate without human interaction.
车辆100中可以包括各种子系统,例如,行进系统110、传感系统120、控制系统130、一个或多个外围设备140以及电源160、计算机系统150和用户接口170。The vehicle 100 may include various subsystems, such as a traveling system 110, a sensing system 120, a control system 130, one or more peripheral devices 140 and a power supply 160, a computer system 150, and a user interface 170.
可选地,车辆100可以包括更多或更少的子系统,并且每个子系统可包括多个元件。另外,车辆100的每个子系统和元件可以通过有线或者无线互连。Alternatively, the vehicle 100 may include more or fewer subsystems, and each subsystem may include multiple elements. In addition, each of the subsystems and elements of the vehicle 100 may be wired or wirelessly interconnected.
示例性地,行进系统110可以包括用于向车辆100提供动力运动的组件。在一个实施例中,行进系统110可以包括引擎111、传动装置112、能量源113和车轮114/轮胎。其中,引擎111可以是内燃引擎、电动机、空气压缩引擎或其他类型的引擎组合;例如,汽油发动机和电动机组成的混动引擎,内燃引擎和空气压缩引擎组成的混动引擎。引擎111可以将能量源113转换成机械能量。Illustratively, the travel system 110 may include components for providing power movement to the vehicle 100. In one embodiment, the travel system 110 may include an engine 111, a transmission 112, an energy source 113, and wheels 114/tires. The engine 111 may be an internal combustion engine, an electric motor, an air compression engine, or other types of engine combinations; for example, a hybrid engine composed of a gasoline engine and an electric motor, or a hybrid engine composed of an internal combustion engine and an air compression engine. The engine 111 can convert the energy source 113 into mechanical energy.
示例性地,能量源113可以包括汽油、柴油、其他基于石油的燃料、丙烷、其他基于压缩气体的燃料、乙醇、太阳能电池板、电池和其他电力来源。能量源113也可以为车辆100的其他系统提供能量。Exemplarily, the energy source 113 may include gasoline, diesel, other petroleum-based fuels, propane, other compressed gas-based fuels, ethanol, solar panels, batteries, and other power sources. The energy source 113 may also provide energy for other systems of the vehicle 100.
示例性地,传动装置112可以包括变速箱、差速器和驱动轴;其中,传动装置112可以将来自引擎111的机械动力传送到车轮114。Exemplarily, the transmission device 112 may include a gearbox, a differential, and a drive shaft; wherein, the transmission device 112 may transmit mechanical power from the engine 111 to the wheels 114.
在一个实施例中,传动装置112还可以包括其他器件,比如离合器。其中,驱动轴可以包括可耦合到一个或多个车轮114的一个或多个轴。In an embodiment, the transmission device 112 may also include other devices, such as a clutch. Among them, the drive shaft may include one or more shafts that can be coupled to one or more wheels 114.
示例性地,传感系统120可以包括感测关于车辆100周边的环境的信息的若干个传感器。Exemplarily, the sensing system 120 may include several sensors that sense information about the environment around the vehicle 100.
例如,传感系统120可以包括定位系统121(例如,GPS系统、北斗系统或者其他定位系统)、惯性测量单元122(inertial measurement unit,IMU)、雷达123、激光测距仪124以及相机125。传感系统120还可以包括被监视车辆100的内部系统的传感器(例如,车内空气质量监测器、燃油量表、机油温度表等)。来自这些传感器中的一个或多个的传感器数据可用于检测对象及其相应特性(位置、形状、方向、速度等)。这种检测和识别是自主车辆100的安全操作的关键功能。For example, the sensing system 120 may include a positioning system 121 (for example, a GPS system, a Beidou system or other positioning systems), an inertial measurement unit 122 (IMU), a radar 123, a laser rangefinder 124, and a camera 125. The sensing system 120 may also include sensors of the internal system of the monitored vehicle 100 (for example, an in-vehicle air quality monitor, a fuel gauge, an oil temperature gauge, etc.). Sensor data from one or more of these sensors can be used to detect objects and their corresponding characteristics (position, shape, direction, speed, etc.). Such detection and identification are key functions for the safe operation of the autonomous vehicle 100.
其中,定位系统121可以用于估计车辆100的地理位置。IMU122可以用于基于惯性加速度来感测车辆100的位置和朝向变化。在一个实施例中,IMU 122可以是加速度计和陀螺仪的组合。Among them, the positioning system 121 can be used to estimate the geographic location of the vehicle 100. The IMU 122 may be used to sense changes in the position and orientation of the vehicle 100 based on inertial acceleration. In an embodiment, the IMU 122 may be a combination of an accelerometer and a gyroscope.
示例性地,雷达123可以利用无线电信号来感测车辆100的周边环境内的物体。在一些实施例中,除了感测物体以外,雷达123还可用于感测物体的速度和/或前进方向。Exemplarily, the radar 123 may use radio signals to sense objects in the surrounding environment of the vehicle 100. In some embodiments, in addition to sensing the object, the radar 123 may also be used to sense the speed and/or direction of the object.
示例性地,激光测距仪124可以利用激光来感测车辆100所位于的环境中的物体。在一些实施例中,激光测距仪124可以包括一个或多个激光源、激光扫描器以及一个或多个检测器,以及其他系统组件。Exemplarily, the laser rangefinder 124 may use laser light to sense objects in the environment where the vehicle 100 is located. In some embodiments, the laser rangefinder 124 may include one or more laser sources, laser scanners, and one or more detectors, as well as other system components.
示例性地,相机125可以用于捕捉车辆100的周边环境的多个图像。例如,相机125可以是静态相机或视频相机。Exemplarily, the camera 125 may be used to capture multiple images of the surrounding environment of the vehicle 100. For example, the camera 125 may be a still camera or a video camera.
如图2所示,控制系统130为控制车辆100及其组件的操作。控制系统130可以包括各种元件,比如可以包括转向系统131、油门132、制动单元133、计算机视觉系统134、路线控制系统135以及障碍规避系统136。As shown in FIG. 2, the control system 130 controls the operation of the vehicle 100 and its components. The control system 130 may include various elements, such as a steering system 131, a throttle 132, a braking unit 133, a computer vision system 134, a route control system 135, and an obstacle avoidance system 136.
示例性地,转向系统131可以操作来调整车辆100的前进方向。例如,在一个实施例 中可以为方向盘系统。油门132可以用于控制引擎111的操作速度并进而控制车辆100的速度。Illustratively, the steering system 131 may be operated to adjust the forward direction of the vehicle 100. For example, it may be a steering wheel system in one embodiment. The throttle 132 may be used to control the operating speed of the engine 111 and thereby control the speed of the vehicle 100.
示例性地,制动单元133可以用于控制车辆100减速;制动单元133可以使用摩擦力来减慢车轮114。在其他实施例中,制动单元133可以将车轮114的动能转换为电流。制动单元133也可以采取其他形式来减慢车轮114转速从而控制车辆100的速度。Illustratively, the braking unit 133 may be used to control the deceleration of the vehicle 100; the braking unit 133 may use friction to slow down the wheels 114. In other embodiments, the braking unit 133 may convert the kinetic energy of the wheels 114 into electric current. The braking unit 133 may also take other forms to slow down the rotation speed of the wheels 114 to control the speed of the vehicle 100.
如图2所示,计算机视觉系统134可以操作来处理和分析由相机125捕捉的图像以便识别车辆100周边环境中的物体和/或特征。上述物体和/或特征可以包括交通信号、道路边界和障碍物。计算机视觉系统134可以使用物体识别算法、运动中恢复结构(structure from motion,SFM)算法、视频跟踪和其他计算机视觉技术。在一些实施例中,计算机视觉系统134可以用于为环境绘制地图、跟踪物体、估计物体的速度等等。As shown in FIG. 2, the computer vision system 134 may be operable to process and analyze the images captured by the camera 125 in order to identify objects and/or features in the surrounding environment of the vehicle 100. The aforementioned objects and/or features may include traffic signals, road boundaries and obstacles. The computer vision system 134 may use object recognition algorithms, structure from motion (SFM) algorithms, video tracking, and other computer vision technologies. In some embodiments, the computer vision system 134 may be used to map the environment, track objects, estimate the speed of objects, and so on.
示例性地,路线控制系统135可以用于确定车辆100的行驶路线。在一些实施例中,路线控制系统135可结合来自传感器、GPS和一个或多个预定地图的数据以为车辆100确定行驶路线。Illustratively, the route control system 135 may be used to determine the travel route of the vehicle 100. In some embodiments, the route control system 135 may combine data from sensors, GPS, and one or more predetermined maps to determine a travel route for the vehicle 100.
如图2所示,障碍规避系统136可以用于识别、评估和避免或者以其他方式越过车辆100的环境中的潜在障碍物。As shown in FIG. 2, the obstacle avoidance system 136 may be used to identify, evaluate, and avoid or otherwise cross potential obstacles in the environment of the vehicle 100.
在一个实例中,控制系统130可以增加或替换地包括除了所示出和描述的那些以外的组件。或者也可以减少一部分上述示出的组件。In one example, the control system 130 may additionally or alternatively include components other than those shown and described. Alternatively, a part of the components shown above may be reduced.
如图2所示,车辆100可以通过外围设备140与外部传感器、其他车辆、其他计算机系统或用户之间进行交互;其中,外围设备140可包括无线通信系统141、车载电脑142、麦克风143和/或扬声器144。As shown in FIG. 2, the vehicle 100 can interact with external sensors, other vehicles, other computer systems, or users through a peripheral device 140; wherein, the peripheral device 140 can include a wireless communication system 141, an onboard computer 142, a microphone 143 and/ Or speaker 144.
在一些实施例中,外围设备140可以提供车辆100与用户接口170交互的手段。例如,车载电脑142可以向车辆100的用户提供信息。用户接口116还可操作车载电脑142来接收用户的输入;车载电脑142可以通过触摸屏进行操作。在其他情况中,外围设备140可以提供用于车辆100与位于车内的其它设备通信的手段。例如,麦克风143可以从车辆100的用户接收音频(例如,语音命令或其他音频输入)。类似地,扬声器144可以向车辆100的用户输出音频。In some embodiments, the peripheral device 140 may provide a means for the vehicle 100 to interact with the user interface 170. For example, the onboard computer 142 may provide information to the user of the vehicle 100. The user interface 116 can also operate the onboard computer 142 to receive user input; the onboard computer 142 can be operated through a touch screen. In other cases, the peripheral device 140 may provide a means for the vehicle 100 to communicate with other devices located in the vehicle. For example, the microphone 143 may receive audio (eg, voice commands or other audio input) from the user of the vehicle 100. Similarly, the speaker 144 may output audio to the user of the vehicle 100.
如图2所述,无线通信系统141可以直接地或者经由通信网络来与一个或多个设备无线通信。例如,无线通信系统141可以使用3G蜂窝通信;例如,码分多址(code division multiple access,CDMA))、EVD0、全球移动通信系统(global system for mobile communications,GSM)/通用分组无线服务(general packet radio service,GPRS),或者4G蜂窝通信,例如长期演进(long term evolution,LTE);或者,5G蜂窝通信。无线通信系统141可以利用无线上网(WiFi)与无线局域网(wireless local area network,WLAN)通信。As shown in FIG. 2, the wireless communication system 141 may wirelessly communicate with one or more devices directly or via a communication network. For example, the wireless communication system 141 can use 3G cellular communication; for example, code division multiple access (CDMA), EVD0, global system for mobile communications (GSM)/general packet radio service (general packet radio service) packet radio service, GPRS), or 4G cellular communication, such as long term evolution (LTE); or, 5G cellular communication. The wireless communication system 141 can communicate with a wireless local area network (WLAN) by using wireless Internet access (WiFi).
在一些实施例中,无线通信系统141可以利用红外链路、蓝牙或者紫蜂协议(ZigBee)与设备直接通信;其他无线协议,例如各种车辆通信系统,例如,无线通信系统141可以包括一个或多个专用短程通信(dedicated short range communications,DSRC)设备,这些设备可包括车辆和/或路边台站之间的公共和/或私有数据通信。In some embodiments, the wireless communication system 141 may directly communicate with the device using an infrared link, Bluetooth, or ZigBee; other wireless protocols, such as various vehicle communication systems, for example, the wireless communication system 141 may include one or Multiple dedicated short range communications (DSRC) devices, these devices may include public and/or private data communications between vehicles and/or roadside stations.
如图2所示,电源160可以向车辆100的各种组件提供电力。在一个实施例中,电源160可以为可再充电锂离子或铅酸电池。这种电池的一个或多个电池组可被配置为电源为 车辆100的各种组件提供电力。在一些实施例中,电源160和能量源113可一起实现,例如一些全电动车中那样。As shown in FIG. 2, the power supply 160 may provide power to various components of the vehicle 100. In one embodiment, the power source 160 may be a rechargeable lithium ion or lead-acid battery. One or more battery packs of such batteries may be configured as a power source to provide power to various components of the vehicle 100. In some embodiments, the power source 160 and the energy source 113 may be implemented together, such as in some all-electric vehicles.
示例性地,车辆100的部分或所有功能可以受计算机系统150控制,其中,计算机系统150可以包括至少一个处理器151,处理器151执行存储在例如存储器152中的非暂态计算机可读介质中的指令153。计算机系统150还可以是采用分布式方式控制车辆100的个体组件或子系统的多个计算设备。Exemplarily, part or all of the functions of the vehicle 100 may be controlled by the computer system 150, where the computer system 150 may include at least one processor 151, and the processor 151 is executed in a non-transitory computer readable medium stored in the memory 152, for example. The instruction 153. The computer system 150 may also be multiple computing devices that control individual components or subsystems of the vehicle 100 in a distributed manner.
例如,处理器151可以是任何常规的处理器,诸如商业可获得的CPU。For example, the processor 151 may be any conventional processor, such as a commercially available CPU.
可选地,该处理器可以是诸如ASIC或其它基于硬件的处理器的专用设备。尽管图2功能性地图示了处理器、存储器、和在相同块中的计算机的其它元件,但是本领域的普通技术人员应该理解该处理器、计算机、或存储器实际上可以包括可以或者可以不存储在相同的物理外壳内的多个处理器、计算机或存储器。例如,存储器可以是硬盘驱动器或位于不同于计算机的外壳内的其它存储介质。因此,对处理器或计算机的引用将被理解为包括对可以或者可以不并行操作的处理器或计算机或存储器的集合的引用。不同于使用单一的处理器来执行此处所描述的步骤,诸如转向组件和减速组件的一些组件每个都可以具有其自己的处理器,所述处理器只执行与特定于组件的功能相关的计算。Alternatively, the processor may be a dedicated device such as an ASIC or other hardware-based processor. Although FIG. 2 functionally illustrates the processor, the memory, and other elements of the computer in the same block, those of ordinary skill in the art should understand that the processor, computer, or memory may or may not actually include Multiple processors, computers or memories in the same physical enclosure. For example, the memory may be a hard disk drive or other storage medium located in a housing other than the computer. Therefore, a reference to a processor or computer will be understood to include a reference to a collection of processors or computers or memories that may or may not operate in parallel. Rather than using a single processor to perform the steps described here, some components such as steering components and deceleration components may each have its own processor that only performs calculations related to component-specific functions .
在此处所描述的各个方面中,处理器可以位于远离该车辆并且与该车辆进行无线通信。在其它方面中,此处所描述的过程中的一些在布置于车辆内的处理器上执行而其它则由远程处理器执行,包括采取执行单一操纵的必要步骤。In the various aspects described herein, the processor may be located away from the vehicle and wirelessly communicate with the vehicle. In other aspects, some of the processes described herein are executed on a processor disposed in the vehicle and others are executed by a remote processor, including taking the necessary steps to perform a single manipulation.
在一些实施例中,存储器152可包含指令153(例如,程序逻辑),指令153可以被处理器151执行来执行车辆100的各种功能,包括以上描述的那些功能。存储器152也可包含额外的指令,比如包括向行进系统110、传感系统120、控制系统130和外围设备140中的一个或多个发送数据、从其接收数据、与其交互和/或对其进行控制的指令。In some embodiments, the memory 152 may contain instructions 153 (eg, program logic), which may be executed by the processor 151 to perform various functions of the vehicle 100, including those functions described above. The memory 152 may also contain additional instructions, for example, including sending data to, receiving data from, interacting with, and/or performing data to one or more of the traveling system 110, the sensing system 120, the control system 130, and the peripheral device 140. Control instructions.
示例性地,除了指令153以外,存储器152还可存储数据,例如,道路地图、路线信息,车辆的位置、方向、速度以及其它这样的车辆数据,以及其他信息。这种信息可在车辆100在自主、半自主和/或手动模式中操作期间被车辆100和计算机系统150使用。Exemplarily, in addition to the instructions 153, the memory 152 may also store data, such as road maps, route information, the position, direction, and speed of the vehicle, and other such vehicle data, as well as other information. Such information may be used by the vehicle 100 and the computer system 150 during the operation of the vehicle 100 in autonomous, semi-autonomous, and/or manual modes.
如图2所示,用户接口170可以用于向车辆100的用户提供信息或从其接收信息。可选地,用户接口170可以包括在外围设备140的集合内的一个或多个输入/输出设备,例如,无线通信系统141、车载电脑142、麦克风143和扬声器144。As shown in FIG. 2, the user interface 170 may be used to provide information to or receive information from a user of the vehicle 100. Optionally, the user interface 170 may include one or more input/output devices in the set of peripheral devices 140, for example, a wireless communication system 141, a car computer 142, a microphone 143, and a speaker 144.
在本申请的实施例中,计算机系统150可以基于从各种子系统(例如,行进系统110、传感系统120和控制系统130)以及从用户接口170接收的输入来控制车辆100的功能。例如,计算机系统150可以利用来自控制系统130的输入以便控制制动单元133来避免由传感系统120和障碍规避系统136检测到的障碍物。在一些实施例中,计算机系统150可操作来对车辆100及其子系统的许多方面提供控制。In an embodiment of the present application, the computer system 150 may control the functions of the vehicle 100 based on inputs received from various subsystems (for example, the traveling system 110, the sensing system 120, and the control system 130) and from the user interface 170. For example, the computer system 150 may use input from the control system 130 in order to control the braking unit 133 to avoid obstacles detected by the sensing system 120 and the obstacle avoidance system 136. In some embodiments, the computer system 150 is operable to provide control of many aspects of the vehicle 100 and its subsystems.
可选地,上述这些组件中的一个或多个可与车辆100分开安装或关联。例如,存储器152可以部分或完全地与车辆100分开存在。上述组件可以按有线和/或无线方式来通信地耦合在一起。Optionally, one or more of these components described above may be installed or associated with the vehicle 100 separately. For example, the storage 152 may exist partially or completely separately from the vehicle 100. The above-mentioned components may be communicatively coupled together in a wired and/or wireless manner.
可选地,上述组件只是一个示例,实际应用中,上述各个模块中的组件有可能根据实际需要增添或者删除,图2不应理解为对本申请实施例的限制。Optionally, the above-mentioned components are only an example. In actual applications, the components in the above-mentioned modules may be added or deleted according to actual needs. FIG. 2 should not be construed as a limitation to the embodiment of the present application.
可选地,车辆100可以是在道路行进的自动驾驶汽车,可以识别其周围环境内的物体 以确定对当前速度的调整。物体可以是其它车辆、交通控制设备、或者其它类型的物体。在一些示例中,可以独立地考虑每个识别的物体,并且基于物体的各自的特性,诸如它的当前速度、加速度、与车辆的间距等,可以用来确定自动驾驶汽车所要调整的速度。Optionally, the vehicle 100 may be an autonomous vehicle traveling on a road, and may recognize objects in its surrounding environment to determine the adjustment to the current speed. The object may be other vehicles, traffic control equipment, or other types of objects. In some examples, each recognized object can be considered independently, and based on the respective characteristics of the object, such as its current speed, acceleration, distance from the vehicle, etc., can be used to determine the speed to be adjusted by the self-driving car.
可选地,车辆100或者与车辆100相关联的计算设备(如图2的计算机系统150、计算机视觉系统134、存储器152)可以基于所识别的物体的特性和周围环境的状态(例如,交通、雨、道路上的冰等等)来预测所述识别的物体的行为。Optionally, the vehicle 100 or a computing device associated with the vehicle 100 (such as the computer system 150, the computer vision system 134, and the memory 152 as shown in FIG. 2) may be based on the characteristics of the identified object and the state of the surrounding environment (for example, traffic, Rain, ice on the road, etc.) to predict the behavior of the identified object.
可选地,每一个所识别的物体都依赖于彼此的行为,因此,还可以将所识别的所有物体全部一起考虑来预测单个识别的物体的行为。车辆100能够基于预测的所述识别的物体的行为来调整它的速度。换句话说,自动驾驶汽车能够基于所预测的物体的行为来确定车辆将需要调整到(例如,加速、减速、或者停止)稳定状态。在这个过程中,也可以考虑其它因素来确定车辆100的速度,诸如,车辆100在行驶的道路中的横向位置、道路的曲率、静态和动态物体的接近度等等。Optionally, each recognized object depends on each other's behavior. Therefore, all recognized objects can also be considered together to predict the behavior of a single recognized object. The vehicle 100 can adjust its speed based on the predicted behavior of the identified object. In other words, the self-driving car can determine based on the predicted behavior of the object that the vehicle will need to be adjusted (e.g., accelerate, decelerate, or stop) to a stable state. In this process, other factors may also be considered to determine the speed of the vehicle 100, such as the lateral position of the vehicle 100 on the road on which it is traveling, the curvature of the road, the proximity of static and dynamic objects, and so on.
除了提供调整自动驾驶汽车的速度的指令之外,计算设备还可以提供修改车辆100的转向角的指令,以使得自动驾驶汽车遵循给定的轨迹和/或维持与自动驾驶汽车附近的物体(例如,道路上的相邻车道中的轿车)的安全横向和纵向距离。In addition to providing instructions to adjust the speed of the self-driving car, the computing device can also provide instructions to modify the steering angle of the vehicle 100 so that the self-driving car follows a given trajectory and/or maintains an object near the self-driving car (for example, , The safe horizontal and vertical distances of cars in adjacent lanes on the road.
上述车辆100可以为轿车、卡车、摩托车、公共汽车、船、飞机、直升飞机、割草机、娱乐车、游乐场车辆、施工设备、电车、高尔夫球车、火车、和手推车等,本申请实施例不做特别的限定。The above-mentioned vehicle 100 may be a car, truck, motorcycle, bus, boat, airplane, helicopter, lawn mower, recreational vehicle, playground vehicle, construction equipment, tram, golf cart, train, and trolley, etc. The application examples are not particularly limited.
在一种可能的实现方式中,上述图2所示的车辆100可以是自动驾驶车辆,下面对自动驾驶系统的进行详细描述。In a possible implementation manner, the vehicle 100 shown in FIG. 2 may be an automatic driving vehicle, and the automatic driving system will be described in detail below.
图3是本申请实施例提供的自动驾驶系统的示意图。Fig. 3 is a schematic diagram of an automatic driving system provided by an embodiment of the present application.
如图3所示的自动驾驶系统包括计算机系统201,其中,计算机系统201包括处理器203,处理器203和系统总线205耦合。处理器203可以是一个或者多个处理器,其中,每个处理器都可以包括一个或多个处理器核。显示适配器207(video adapter),显示适配器可以驱动显示器209,显示器209和系统总线205耦合。系统总线205可以通过总线桥211和输入输出(I/O)总线213耦合,I/O接口215和I/O总线耦合。I/O接口215和多种I/O设备进行通信,比如,输入设备217(如:键盘,鼠标,触摸屏等),媒体盘221(media tray),(例如,CD-ROM,多媒体接口等)。收发器223可以发送和/或接受无线电通信信号,摄像头255可以捕捉景田和动态数字视频图像。其中,和I/O接口215相连接的接口可以是USB端口225。The automatic driving system shown in FIG. 3 includes a computer system 201, where the computer system 201 includes a processor 203, and the processor 203 is coupled to a system bus 205. The processor 203 may be one or more processors, where each processor may include one or more processor cores. The display adapter 207 (video adapter) can drive the display 209, and the display 209 is coupled to the system bus 205. The system bus 205 may be coupled to an input/output (I/O) bus 213 through a bus bridge 211, and an I/O interface 215 is coupled to an I/O bus. The I/O interface 215 communicates with a variety of I/O devices, such as input devices 217 (such as keyboard, mouse, touch screen, etc.), media tray 221, (such as CD-ROM, multimedia interface, etc.) . The transceiver 223 can send and/or receive radio communication signals, and the camera 255 can capture landscape and dynamic digital video images. Wherein, the interface connected to the I/O interface 215 may be the USB port 225.
其中,处理器203可以是任何传统处理器,比如,精简指令集计算(reduced instruction set computer,RISC)处理器、复杂指令集计算(complex instruction set computer,CISC)处理器或上述的组合。The processor 203 may be any traditional processor, such as a reduced instruction set computer (RISC) processor, a complex instruction set computer (CISC) processor, or a combination of the foregoing.
可选地,处理器203可以是诸如专用集成电路(application specific integrated circuit,ASIC)的专用装置;处理器203可以是神经网络处理器或者是神经网络处理器和上述传统处理器的组合。Optionally, the processor 203 may be a dedicated device such as an application specific integrated circuit (ASIC); the processor 203 may be a neural network processor or a combination of a neural network processor and the foregoing traditional processors.
可选地,在本文所述的各种实施例中,计算机系统201可位于远离自动驾驶车辆的地方,并且可与自动驾驶车辆无线通信。在其它方面,本文所述的一些过程在设置在自动驾驶车辆内的处理器上执行,其它由远程处理器执行,包括采取执行单个操纵所需的动作。Optionally, in various embodiments described herein, the computer system 201 may be located far away from the autonomous driving vehicle, and may wirelessly communicate with the autonomous driving vehicle. In other respects, some of the processes described herein are executed on a processor provided in an autonomous vehicle, and others are executed by a remote processor, including taking actions required to perform a single manipulation.
计算机系统201可以通过网络接口229和软件部署服务器249通信。网络接口229可以是硬件网络接口,比如,网卡。网络227可以是外部网络,比如,因特网,也可以是内部网络,比如以太网或者虚拟私人网络(virtual private network,VPN)。可选地,网络227还可以是无线网络,比如wifi网络,蜂窝网络等。The computer system 201 can communicate with the software deployment server 249 through the network interface 229. The network interface 229 may be a hardware network interface, such as a network card. The network 227 may be an external network, such as the Internet, or an internal network, such as an Ethernet or a virtual private network (VPN). Optionally, the network 227 may also be a wireless network, such as a wifi network, a cellular network, and so on.
如图3所示,硬盘驱动接口和系统总线205耦合,硬件驱动器接口231可以与硬盘驱动器233相连接,系统内存235和系统总线205耦合。运行在系统内存235的数据可以包括操作系统237和应用程序243。其中,操作系统237可以包括解析器239(shell)和内核241(kernel)。shell 239是介于使用者和操作系统之内核(kernel)间的一个接口。Shell可以是操作系统最外面的一层;shell可以管理使用者与操作系统之间的交互,比如,等待使用者的输入,向操作系统解释使用者的输入,并且处理各种各样的操作系统的输出结果。内核241可以由操作系统中用于管理存储器、文件、外设和系统资源的那些部分组成。直接与硬件交互,操作系统内核通常运行进程,并提供进程间的通信,提供CPU时间片管理、中断、内存管理、IO管理等等。应用程序243包括控制汽车自动驾驶相关的程序,比如,管理自动驾驶的汽车和路上障碍物交互的程序,控制自动驾驶汽车路线或者速度的程序,控制自动驾驶汽车和路上其他自动驾驶汽车交互的程序。应用程序243也存在于软件部署服务器249的系统上。在一个实施例中,在需要执行自动驾驶相关程序247时,计算机系统201可以从软件部署服务器249下载应用程序。As shown in FIG. 3, the hard disk drive interface is coupled with the system bus 205, the hardware drive interface 231 can be connected with the hard drive 233, and the system memory 235 is coupled with the system bus 205. The data running in the system memory 235 may include an operating system 237 and application programs 243. The operating system 237 may include a parser 239 (shell) and a kernel 241 (kernel). The shell 239 is an interface between the user and the kernel of the operating system. The shell can be the outermost layer of the operating system; the shell can manage the interaction between the user and the operating system, for example, waiting for the user's input, interpreting the user's input to the operating system, and processing various operating systems The output result. The kernel 241 may be composed of those parts of the operating system that are used to manage memory, files, peripherals, and system resources. Directly interact with the hardware. The operating system kernel usually runs processes and provides inter-process communication, providing CPU time slice management, interrupts, memory management, IO management, and so on. Application programs 243 include programs that control auto-driving cars, such as programs that manage the interaction between autonomous vehicles and obstacles on the road, programs that control the route or speed of autonomous vehicles, and programs that control interaction between autonomous vehicles and other autonomous vehicles on the road. . The application program 243 also exists on the system of the software deployment server 249. In one embodiment, the computer system 201 may download the application program from the software deployment server 249 when the automatic driving-related program 247 needs to be executed.
示例性地,传感器253可以与计算机系统201关联,传感器253可以用于探测计算机201周围的环境。Exemplarily, the sensor 253 may be associated with the computer system 201, and the sensor 253 may be used to detect the environment around the computer 201.
举例来说,传感器253可以探测动物,汽车,障碍物和人行横道等,进一步传感器还可以探测上述动物,汽车,障碍物和人行横道等物体周围的环境,比如:动物周围的环境,例如,动物周围出现的其他动物,天气条件,周围环境的光亮度等。For example, the sensor 253 can detect animals, cars, obstacles, and crosswalks. Further, the sensor can also detect the surrounding environment of the above-mentioned animals, cars, obstacles, and crosswalks, such as: the environment around the animals, for example, when the animals appear around them. Other animals, weather conditions, the brightness of the surrounding environment, etc.
可选地,如果计算机系统201位于自动驾驶的汽车上,传感器可以是摄像头,红外线感应器,化学检测器,麦克风等。Optionally, if the computer system 201 is located in an auto-driving car, the sensor may be a camera, an infrared sensor, a chemical detector, a microphone, etc.
示例性地,传感器253可以为多个。多个传感器可以用于探测车辆周围障碍物的位置,基于多个传感器获取的数据得到障碍物的位置。具体地,基于多个传感器获取的数据得到障碍物的位置可以由本申请实施例的目标跟踪方法实现。Exemplarily, there may be multiple sensors 253. Multiple sensors can be used to detect the location of obstacles around the vehicle, and the location of the obstacles can be obtained based on the data obtained by the multiple sensors. Specifically, obtaining the location of the obstacle based on the data acquired by multiple sensors can be implemented by the target tracking method in the embodiment of the present application.
在一个示例中,图2所示的计算机系统150还可以从其它计算机系统接收信息或转移信息到其它计算机系统。或者,从车辆100的传感系统120收集的传感器数据可以被转移到另一个计算机对此数据进行处理。In an example, the computer system 150 shown in FIG. 2 may also receive information from other computer systems or transfer information to other computer systems. Alternatively, the sensor data collected from the sensor system 120 of the vehicle 100 may be transferred to another computer to process the data.
例如,如图4所示,来自计算机系统312的数据可以经由网络被传送到云侧的服务器320用于进一步的处理。网络以及中间节点可以包括各种配置和协议,包括因特网、万维网、内联网、虚拟专用网络、广域网、局域网、使用一个或多个公司的专有通信协议的专用网络、以太网、WiFi和HTTP、以及前述的各种组合;这种通信可以由能够传送数据到其它计算机和从其它计算机传送数据的任何设备,诸如调制解调器和无线接口。For example, as shown in FIG. 4, data from the computer system 312 may be transmitted to the server 320 on the cloud side via the network for further processing. The network and intermediate nodes can include various configurations and protocols, including the Internet, World Wide Web, Intranet, virtual private network, wide area network, local area network, private network using one or more company’s proprietary communication protocols, Ethernet, WiFi and HTTP, And various combinations of the foregoing; this communication can be by any device capable of transferring data to and from other computers, such as modems and wireless interfaces.
在一个示例中,服务器320可以包括具有多个计算机的服务器,例如负载均衡服务器群,为了从计算机系统312接收、处理并传送数据的目的,其与网络的不同节点交换信息。该服务器可以被类似于计算机系统312配置,具有处理器330、存储器340、指令350、和数据360。In one example, the server 320 may include a server with multiple computers, such as a load balancing server group, which exchanges information with different nodes of the network for the purpose of receiving, processing, and transmitting data from the computer system 312. The server may be configured similarly to the computer system 312, with a processor 330, a memory 340, instructions 350, and data 360.
示例性地,服务器320的数据360可以包括车辆周围道路情况的相关信息。例如,服务器320可以接收、检测、存储、更新、以及传送与车辆道路情况的相关信息。Exemplarily, the data 360 of the server 320 may include information related to road conditions around the vehicle. For example, the server 320 may receive, detect, store, update, and transmit information related to the road conditions of the vehicle.
下面对贝叶斯滤波器的相关内容做详细介绍。The following is a detailed introduction to the relevant content of the Bayesian filter.
状态空间是指描述一个系统所有可能出现的状态的状态变量集合。汽车即可以看作一个系统,用户对汽车的操作可以看作输入变量,有操作信号输入时,会对汽车的速度、加速度、角速度等变量产生明确的影响,这些受影响的变量都可以作为系统的状态变量分量。State space refers to a collection of state variables that describe all possible states of a system. The car can be regarded as a system, and the operation of the car by the user can be regarded as input variables. When there is an input of an operation signal, it will have a clear impact on the speed, acceleration, angular velocity and other variables of the car. These affected variables can all be regarded as the system. The component of the state variable.
观测指的是,通过某种量测手段直接或间接获取状态变量估计值的过程。Observation refers to the process of obtaining state variable estimates directly or indirectly through a certain measurement method.
贝叶斯估计指的是,从任意时刻k-1开始,先计算下一时刻k的状态先验估计的概率分布,称之为预测;然后在获取k时刻的观测值以后,对预测环节得到的先验估计进行修正,得到了k时刻状态的后验估计,称之为更新。Bayesian estimation means that starting from any time k-1, the probability distribution of the state a priori estimate at the next time k is calculated, which is called prediction; then after obtaining the observation value at time k, the prediction link is obtained The a priori estimate of is revised, and the posterior estimate of the state at time k is obtained, which is called update.
贝叶斯递推估计应用于实际工程中称为贝叶斯滤波。Bayesian recursive estimation is called Bayesian filtering when used in actual engineering.
贝叶斯滤波器分为预测和更新两个部分。The Bayesian filter is divided into two parts: prediction and update.
预测(predict)过程满足公式(1):The predict process satisfies formula (1):
f k+1|k(X|z k)=∫f k+1|k(X|X')f k|k(X'|z k)dX'  (1) f k+1|k (X|z k )=∫f k+1|k (X|X')f k|k (X'|z k )dX' (1)
更新(update)过程满足公式(2)和公式(3):The update process satisfies formula (2) and formula (3):
Figure PCTCN2020075556-appb-000017
Figure PCTCN2020075556-appb-000017
f k+1(z k+1|z k)=∫f k+1(z k+1|X)f k+1|k(X|z k)dX  (3) f k+1 (z k+1 |z k )=∫f k+1 (z k+1 |X)f k+1|k (X|z k )dX (3)
其中,X为状态量,即贝叶斯滤波器的输出量,状态量所在的空间即为状态空间。z为观测量,即贝叶斯滤波器的输入量,观测量所在的空间即为观测空间。z k+1为第k+1次的观测量,z k={z 0,z 1,...,z k}为前k次观测量的集合;f(*)为概率密度函数。 Among them, X is the state quantity, that is, the output quantity of the Bayesian filter, and the space where the state quantity is located is the state space. z is the observation, that is, the input of the Bayesian filter, and the space where the observation is located is the observation space. z k+1 is the k+1 observation, z k = {z 0 , z 1 ,..., z k } is the set of the previous k observations; f(*) is the probability density function.
图5是本申请实施例的目标跟踪装置的示意图。该目标跟踪装置402可以应用于计算机系统401中。该目标跟踪装置包括预测单元410和更新单元420。其中,预测单元410包括判断模块411,更新单元420包括观测量转移模块422。Fig. 5 is a schematic diagram of a target tracking device according to an embodiment of the present application. The target tracking device 402 can be applied to the computer system 401. The target tracking device includes a prediction unit 410 and an update unit 420. The prediction unit 410 includes a judgment module 411, and the update unit 420 includes an observation transfer module 422.
进一步地,预测单元410还可以包括状态转移模块412。进一步地,更新单元420还可以包括更新模块421。Further, the prediction unit 410 may also include a state transition module 412. Further, the update unit 420 may also include an update module 421.
为了更好地了解本申请实施例的目标跟踪方法的执行过程,下面先对图5中的各个模块的功能进行简单的描述。In order to better understand the execution process of the target tracking method in the embodiment of the present application, the function of each module in FIG. 5 is briefly described below.
判断模块411,用于判断观测量的采集时刻是否早于当前状态量的更新时刻。也就是判断该观测量是否为迟到的观测量。当前状态量指的是最新一次更新的状态量。具体地,当该观测量的采集时刻早于最新一次状态量的更新时刻,则该观测量为迟到的观测量。The judging module 411 is used to judge whether the collection time of the observation is earlier than the update time of the current state quantity. That is to judge whether the observation is late. The current state quantity refers to the state quantity of the latest update. Specifically, when the collection time of the observation is earlier than the latest update time of the state quantity, the observation is a late observation.
应理解,在本申请实施例中,由传感器采集得到的观测量即为实际观测量。It should be understood that, in the embodiments of the present application, the observation collected by the sensor is the actual observation.
状态转移模块412,用于在该观测量的采集时刻不早于当前状态量的更新时刻的情况下,根据前一次更新的状态量得到本次更新的状态量的预测值。具体过程满足公式(1)。The state transition module 412 is configured to obtain the predicted value of the state amount of the current update according to the state amount of the previous update when the collection time of the observation is not earlier than the update time of the current state amount. The specific process satisfies formula (1).
更新模块421,用于根据状态转移模块412得到的本次更新的状态量的预测值更新状态量。具体过程满足公式(2)和公式(3)。The update module 421 is configured to update the state quantity according to the predicted value of the current state quantity obtained by the state transition module 412. The specific process satisfies formula (2) and formula (3).
观测量转移单元422,用于在观测量的采集时刻早于当前状态量的更新时刻的情况下,对该观测量进行状态转移,得到估计观测量。观测量转移单元422还用于,根据实际观测量和估计观测量对当前状态量进行更新,将更新后的状态量作为更新时刻对应的状态量。The observation transfer unit 422 is configured to perform state transition on the observation when the acquisition time of the observation is earlier than the update time of the current state quantity to obtain the estimated observation. The observation transfer unit 422 is further configured to update the current state quantity according to the actual observation and the estimated observation, and use the updated state quantity as the state quantity corresponding to the update time.
下面结合图6对本申请实施例的目标跟踪方法500进行详细描述。图6是本申请实施例的目标跟踪方法500的示意性流程图。图6所示的方法可以由本申请实施例中的目标跟踪装置来执行。方法500包括步骤S510至步骤S560。下面对步骤S510至步骤S560进行详细说明。The target tracking method 500 of the embodiment of the present application will be described in detail below in conjunction with FIG. 6. FIG. 6 is a schematic flowchart of a target tracking method 500 according to an embodiment of the present application. The method shown in FIG. 6 may be executed by the target tracking device in the embodiment of the present application. The method 500 includes steps S510 to S560. Steps S510 to S560 will be described in detail below.
S510,获取传感器采集的运动物体的实际观测量。应理解,在本申请实施例中,由传感器采集得到的观测量即为实际观测量,该实际观测量在本申请实施例中也可以称为“观测量”。观测量可以包括运动物体的速度、运动物体的加速度或运动物体的位置等。例如,该观测量可以包括毫米波雷达测量的运动物体的位置和/或速度等。S510: Obtain actual observations of the moving object collected by the sensor. It should be understood that, in the embodiments of the present application, the observations collected by the sensors are the actual observations, and the actual observations may also be referred to as "observations" in the embodiments of the present application. Observation can include the speed of the moving object, the acceleration of the moving object, or the position of the moving object. For example, the observation may include the position and/or speed of the moving object measured by the millimeter wave radar.
S520,判断观测量的采集时刻是否早于运动物体的当前状态量的更新时刻。S520: Determine whether the collection time of the observation is earlier than the update time of the current state quantity of the moving object.
如果该观测量的采集时刻不早于运动物体的当前状态量的更新时刻,则执行步骤S530。如果该观测量的采集时刻早于运动物体的当前状态量的更新时刻,则执行步骤S540。If the collection time of the observation is not earlier than the update time of the current state quantity of the moving object, step S530 is executed. If the collection time of the observation is earlier than the update time of the current state quantity of the moving object, step S540 is executed.
观测量的采集时刻早于运动物体的当前状态量的更新时刻,可以理解为,该观测量原本可以用于运动物体的当前状态量的更新,但实际上未能用于运动物体的当前状态量的更新。The acquisition time of the observation is earlier than the update time of the current state quantity of the moving object. It can be understood that the observation can be used to update the current state quantity of the moving object, but it cannot actually be used for the current state quantity of the moving object. Update.
该观测量的采集时刻可以由该观测量的时间戳指示。The collection time of the observation can be indicated by the time stamp of the observation.
状态量可以包括运动物体的速度、运动物体的加速度或运动物体的位置等。例如,当方法500用于跟踪目标的位置,状态量可以为跟踪的结果,该跟踪的结果可以为运动物体的位置。The state quantity may include the speed of the moving object, the acceleration of the moving object, or the position of the moving object. For example, when the method 500 is used to track the position of a target, the state quantity may be the result of the tracking, and the result of the tracking may be the position of the moving object.
运动物体的当前状态量的更新时刻,可以理解为,运动物体的状态量最新一次更新的时刻。The update time of the current state quantity of the moving object can be understood as the time when the state quantity of the moving object was updated last time.
其中,运动物体的状态量以及状态量的更新时刻可以保存于贝叶斯滤波器中的跟踪列表中。Among them, the state quantity of the moving object and the update time of the state quantity can be stored in the tracking list in the Bayesian filter.
具体地,判断观测量的采集时刻是否早于运动物体的当前状态量的更新时刻,可以为,判断观测量的时间戳是否早于跟踪列表中最新一次更新的时刻。Specifically, judging whether the collection time of the observation is earlier than the update time of the current state quantity of the moving object may be judging whether the time stamp of the observation is earlier than the latest update time in the tracking list.
S530,更新运动物体的当前状态量,将更新后的状态量作为观测量的采集时刻对应的状态量。S530: Update the current state quantity of the moving object, and use the updated state quantity as the state quantity corresponding to the collection moment of the observation.
具体地,可以通过贝叶斯滤波器更新运动物体的当前状态量。例如,可以根据上述公式(1)、公式(2)和公式(3)更新运动物体的当前状态量。Specifically, the current state quantity of the moving object can be updated through the Bayesian filter. For example, the current state quantity of the moving object can be updated according to the above formula (1), formula (2) and formula (3).
S540,判断该观测量的采集时刻与运动物体的当前状态量的更新时刻之间的时间差是否大于第一阈值。S540: Determine whether the time difference between the collection time of the observation and the update time of the current state quantity of the moving object is greater than a first threshold.
如果该时间差大于第一阈值,则丢弃该观测量。也就是不利用该输入观测量进行状态量的更新。这样可以避免在时间差过大的情况下,利用该观测量更新状态量而导致更新的状态量的置信度降低。If the time difference is greater than the first threshold, the observation is discarded. That is, the input observation is not used to update the state quantity. In this way, when the time difference is too large, using the observation to update the state quantity may reduce the confidence of the updated state quantity.
如果该时间差小于或等于该第一阈值,则执行步骤S550。If the time difference is less than or equal to the first threshold, step S550 is executed.
可选地,步骤S540也可以为,判断该观测量的采集时刻与运动物体的当前状态量的更新时刻之间的时间差是否大于或等于第一阈值。Optionally, step S540 may also be to determine whether the time difference between the collection time of the observation and the update time of the current state quantity of the moving object is greater than or equal to the first threshold.
如果该时间差大于或等于第一阈值,则丢弃该观测量。也就是不利用该输入观测量进行状态量的更新。如果该时间差小于该第一阈值,则执行步骤S550。If the time difference is greater than or equal to the first threshold, then the observation is discarded. That is, the input observation is not used to update the state quantity. If the time difference is less than the first threshold, step S550 is executed.
需要说明的是,步骤S520、步骤S530和步骤S540为可选步骤,本申请实施例的方 法500可以在步骤S510之后执行步骤S550。It should be noted that step S520, step S530, and step S540 are optional steps, and the method 500 of the embodiment of the present application may execute step S550 after step S510.
S550,在观测量的采集时刻早于运动物体的当前状态量的更新时刻的情况下,根据运动物体对应的状态转移模型得到估计观测量。S550: In a case where the collection time of the observation is earlier than the update time of the current state quantity of the moving object, obtain the estimated observation according to the state transition model corresponding to the moving object.
当方法500用于跟踪目标的位置,状态转移模型可以为运动物体的运动学模型。例如,该运动学模型可以包括匀速直线运动模型等。When the method 500 is used to track the position of a target, the state transition model may be a kinematics model of a moving object. For example, the kinematics model may include a uniform linear motion model and the like.
具体地,根据运动物体对应的状态转移模型得到关于采集时刻对应的观测量的似然函数,估计观测量是根据关于采集时刻对应的观测量的似然函数确定的,状态转移模型是根据第一概率密度函数确定的,第一概率密度函数是将当前状态量转移至采集时刻对应的状态量的概率密度函数。Specifically, the likelihood function of the observation corresponding to the acquisition time is obtained according to the state transition model corresponding to the moving object. The estimated observation is determined based on the likelihood function of the observation corresponding to the acquisition time. The state transition model is based on the first The probability density function is determined, and the first probability density function is a probability density function that transfers the current state quantity to the state quantity corresponding to the collection time.
下面结合公式对步骤S550进行说明。The step S550 will be described below in conjunction with a formula.
t k时刻采集的观测量z k输入贝叶斯滤波器,贝叶斯滤波器根据该观测量z k更新状态量,将得到的状态量X k'作为当前状态量X k'。相应地,状态量最新一次的更新为t k时刻,也就是当前状态量X k'的更新时刻为t k时刻。在t k时刻之后,t k-1时刻采集的观测量z k-1输入贝叶斯滤波器。上述观测量的采集时刻早于当前状态量的更新时刻可以理解为,贝叶斯滤波器接收到观测量z k并更新了状态量之后又接收到观测量z k-1 The observation z k collected at the time t k is input to the Bayesian filter, and the Bayesian filter updates the state quantity according to the observation z k , and uses the obtained state quantity X k 'as the current state quantity X k '. Correspondingly, the latest update of the state quantity is time t k , that is, the update time of the current state quantity X k ′ is time t k . After the time t k, t k-1 time View captured measurement z k-1 input Bayesian filter. The acquisition time of the above observation is earlier than the update time of the current state quantity. It can be understood that the Bayesian filter receives the observation z k and updates the state quantity and then receives the observation z k-1 .
如果观测量z k-1输入贝叶斯滤波器的时刻早于观测量z k输入贝叶斯滤波器的时刻,则t k-1时刻的似然函数g(z k-1|X k-1)。 If the concept of time measurements z k-1 input to the Bayesian-filter concept of earlier measurement time Bayesian-filter input z k, t is the likelihood function g k-1 time (z k-1 | X k- 1 ).
如果观测量z k-1输入贝叶斯滤波器的时刻不早于观测量z k输入贝叶斯滤波器的时刻,则t k时刻到t k-1时刻的状态转移模型为f k-1|k(X|X k'),该状态转移即为马尔科夫转移。f k-1|k(X|X k')为上述第一概率密度函数,第一概率密度函数是将t k时刻更新的状态量X k'转移至到t k-1时刻更新的状态量的概率密度函数。 If the time when the observation z k-1 is input to the Bayesian filter is not earlier than the time when the observation z k is input to the Bayes filter, then the state transition model from time t k to time t k-1 is f k-1 |k (X|X k '), this state transition is Markov transition. f k-1|k (X|X k ') is the above-mentioned first probability density function, and the first probability density function is to transfer the state quantity X k ' updated at t k to the state quantity updated at t k-1 The probability density function.
关于采集时刻对应的观测量的似然函数满足:The likelihood function of the observation corresponding to the acquisition time satisfies:
g(z k-1|X k')=∫g(z k-1|X)f k-1|k(X|X k')dX g(z k-1 |X k ')=∫g(z k-1 |X)f k-1|k (X|X k ')dX
其中,g(z k-1|X k')表示关于t k-1时刻对应的观测量z k-1的似然函数,X表示t k-1时刻更新的状态量,X k'表示在缺少z k-1的情况下t k时刻更新的状态量,g(z k-1|X)表示在z k-1用于t k-1时刻更新状态量的情况下关于z k-1的似然函数。 Among them, g(z k-1 |X k ') represents the likelihood function of the observation z k-1 corresponding to the time t k-1 , X represents the state quantity updated at the time t k-1 , and X k 'represents the missing state quantity update time t K in the case of Z k-1, g (z k-1 | X) Z represents Z k-1 k-1 is used on the time t where k-1 is updated state quantity Likelihood function.
S560,根据实际观测量和估计观测量对当前状态量进行更新,将更新后的状态量作为更新时刻对应的状态量。S560: Update the current state quantity according to the actual observation and the estimated observation, and use the updated state quantity as the state quantity corresponding to the update time.
具体地,可以根据关于所述采集时刻对应的观测量的似然函数对当前状态量的概率密度函数进行更新。根据更新后的状态量的概率密度函数确定更新时刻对应的状态量。Specifically, the probability density function of the current state quantity may be updated according to the likelihood function of the observation corresponding to the acquisition time. The state quantity corresponding to the update time is determined according to the probability density function of the updated state quantity.
下面结合公式对步骤S560进行说明。The step S560 will be described below in conjunction with a formula.
更新后的状态量的概率密度函数满足:The probability density function of the updated state quantity satisfies:
Figure PCTCN2020075556-appb-000018
Figure PCTCN2020075556-appb-000018
其中,f k|k(X k|z k)表示t k时刻更新的状态量X k的概率密度函数,f' k|k(X k'|z k,z k-2)表示在缺少z k-1的情况下t k时刻更新的状态量X k'的概率密度函数。也就是当前状态量X k'对应的后验概率。z k表示t k时刻采集的观测量,z k表示k个时刻采集的观测量的集合{z 1,z 2,...,z k-2,z k-1,z k},z k-2表示k-2个时刻采集的观测量的集合{z 1,z 2,...,z k-2}。 Among them, f k|k (X k |z k ) represents the probability density function of the state quantity X k updated at time t k , and f'k|k (X k '|z k ,z k-2 ) represents the absence of z In the case of k-1 , the probability density function of the state quantity X k 'updated at time t k. That is, the posterior probability corresponding to the current state quantity X k'. z k represents the observations collected at time t k , z k represents the set of observations collected at k time {z 1 ,z 2 ,...,z k-2 ,z k-1 ,z k }, z k -2 represents the set of observations collected at k-2 moments {z 1 , z 2 ,..., z k-2 }.
进一步地,卡尔曼滤波为贝叶斯滤波器的一种实现形式,以卡尔曼滤波器为例,对步 骤S560进行说明。Further, the Kalman filter is an implementation form of the Bayes filter. Taking the Kalman filter as an example, step S560 will be described.
上述步骤S560,可以包括,根据实际观测量和估计观测量确定更新时刻对应的状态量的期望;将更新时刻对应的状态量的期望作为更新时刻对应的状态量。The above step S560 may include determining the expectation of the state quantity corresponding to the update time according to the actual observation and the estimated observation; and taking the expectation of the state quantity corresponding to the update time as the state quantity corresponding to the update time.
其中,更新时刻对应的状态量的期望与由卡尔曼增益值有关。卡尔曼增益值与第一协方差、采集时刻的观测矩阵、观测矩阵的协方差和采集时刻的观测量的方差有关。第一协方差指的是由更新时刻转移至采集时刻的状态量的协方差。Among them, the expectation of the state quantity corresponding to the update time is related to the Kalman gain value. The Kalman gain value is related to the first covariance, the observation matrix at the time of collection, the covariance of the observation matrix, and the variance of the observation at the time of collection. The first covariance refers to the covariance of the state quantity transferred from the update time to the collection time.
可选地,卡尔曼增益值满足:Optionally, the Kalman gain value satisfies:
Figure PCTCN2020075556-appb-000019
Figure PCTCN2020075556-appb-000019
Var(z k-1)满足: Var(z k-1 ) satisfies:
Figure PCTCN2020075556-appb-000020
Figure PCTCN2020075556-appb-000020
P k-1|k满足: P k-1|k satisfies:
Figure PCTCN2020075556-appb-000021
Figure PCTCN2020075556-appb-000021
其中,P k-1|k表示由更新时刻t k转移到采集时刻t k-1的状态量的协方差,P' k|k表示在缺少t k-1时刻采集的观测量z k-1的情况下t k时刻更新的状态量的协方差,F k-1|k表示从t k时刻转移到t k-1时刻的状态转移矩阵,Q k表示预测矩阵的协方差,H k-1表示t k-1时刻的观测矩阵,Var(z k-1)表示t k-1时刻采集的观测量z k-1的方差,R k-1表示观测矩阵的协方差。 Wherein, P k-1 | k represents t K transferred from the update time to the acquisition time t covariance k-1 state quantity, P 'k | k represents the lack Concept t k-1 time captured measurement Z k-1 In the case of t k , the covariance of the state quantity updated at time t k, F k-1|k represents the state transition matrix from t k to t k-1 , Q k represents the covariance of the prediction matrix, H k-1 Represents the observation matrix at time t k-1 , Var(z k-1 ) represents the variance of the observation z k-1 collected at time t k-1 , and R k-1 represents the covariance of the observation matrix.
更新时刻对应的状态量的期望满足:The expectation of the state quantity corresponding to the update time meets:
Figure PCTCN2020075556-appb-000022
Figure PCTCN2020075556-appb-000022
其中,x k|k表示t k时刻更新的状态量的期望,
Figure PCTCN2020075556-appb-000023
表示t k-1时刻的估计观测量。将更新时刻对应的状态量的期望x k|k作为更新时刻对应的状态量。
Among them, x k|k represents the expectation of the state quantity updated at time t k,
Figure PCTCN2020075556-appb-000023
Represents the estimated observation at time t k-1. The expected state quantity x k|k corresponding to the update time is taken as the state quantity corresponding to the update time.
可选地,采集时刻对应的观测量的估计值满足:Optionally, the estimated value of the observation corresponding to the collection time satisfies:
Figure PCTCN2020075556-appb-000024
Figure PCTCN2020075556-appb-000024
其中,
Figure PCTCN2020075556-appb-000025
表示t k-1时刻的估计观测量,H k-1表示t k-1时刻的观测矩阵,F k-1|k表示从t k时刻转移到t k-1时刻的状态转移矩阵,x' k|k表示在缺少所述t k-1时刻采集的观测量z k-1的情况下的t k时刻更新的状态量的期望。
in,
Figure PCTCN2020075556-appb-000025
Represents t estimated observations k-1 time point, H k-1 represents a t the observation matrix k-1 time, F k-1 | k indicates a transition from t K time to t k-1 time of the state transition matrix, x ' k | k t k represents a time measurement in the case where z k-1 is updated in a desired state quantity of the absence of the concept of time t k-1 collected.
Figure PCTCN2020075556-appb-000026
表示t k-1时刻的估计观测量可以理解为,通过状态转移矩阵F k-1|k将状态量的期望x' k|k转移至t k-1时刻,然后基于观测矩阵H k-1得到t k-1时刻的观测量的期望
Figure PCTCN2020075556-appb-000027
作为所述t k-1时刻的估计观测量。
Figure PCTCN2020075556-appb-000026
The estimated observation at time t k-1 can be understood as the state transition matrix F k-1|k to transfer the expected state quantity x'k|k to time t k-1 , and then based on the observation matrix H k-1 Obtain the expectation of the observation at t k-1
Figure PCTCN2020075556-appb-000027
As the estimated observation at the time t k-1.
进一步地,计算t k时刻更新的状态量的协方差。t k时刻更新的状态量的协方差满足: Further, the covariance of the updated state quantity at time t k is calculated. The covariance of the updated state at t k satisfies:
P k|k=(I-K' kH k-1)P k-1|k P k|k = (IK' k H k-1 )P k-1|k
其中,P k|k表示t k时刻更新的状态量的协方差。 Among them, P k|k represents the covariance of the state quantity updated at time t k.
根据本申请实施例的方案,根据运动物体对应的状态转移模型进行状态转移,进而对当前状态量进行更新,能够得到最新的更新时刻对应的较准确的状态量,而不是采集时刻对应的状态量,避免输出的状态量在时间上出现乱序跳动。According to the solution of the embodiment of the present application, the state transition is performed according to the state transition model corresponding to the moving object, and then the current state quantity is updated, and the more accurate state quantity corresponding to the latest update time can be obtained, instead of the state quantity corresponding to the collection time. , To avoid out-of-order jumps of the output state quantity in time.
下面以方法500应用于跟踪目标为例,对步骤S550和步骤S560进行说明。Hereinafter, taking the method 500 applied to tracking a target as an example, step S550 and step S560 will be described.
R表示传感器与目标之间的距离。x=(R)表示一个一维向量。R represents the distance between the sensor and the target. x=(R) represents a one-dimensional vector.
在缺少所述t k-1时刻采集的观测量z k-1的情况下的t k时刻更新的状态量x k满足
Figure PCTCN2020075556-appb-000028
x k~N(x;x′ k|k,P′ k|k)表示传感器与目标之间的距离为一个期望为x' k|k,方差为P′ k|k的随机数。P' k|k可以为1。
In the absence of the time t k-1 Concept captured measurement time t K in the case where Z k-1 is updated to meet the state quantity x k
Figure PCTCN2020075556-appb-000028
x k ~N(x; x′ k|k ,P′ k|k ) indicates that the distance between the sensor and the target is a random number with an expectation of x′ k|k and a variance of P′ k|k. P'k|k can be 1.
x' k|k表示在缺少所述t k-1时刻采集的观测量z k-1的情况下的t k时刻更新的状态量的期望。P' k|k表示在缺少t k-1时刻采集的观测量z k-1的情况下t k时刻更新的状态量的协方差。 x 'k | k t k represents a time measurement in the case where z k-1 is updated in a desired state quantity of the absence of the concept of time t k-1 collected. P 'k | k represents the missing covariance Concept time t k-1 measured state quantity acquisition case where z k-1 t k is the time updates.
例如,x' k|k可以为3,P' k|k可以为1,t k-1时刻采集的观测量z k-1可以为4。在该情况下,P k-1|k满足: For example, x 'k | k can be 3, P' k | k may be 1, t k-1 time View captured measurement z k-1 may be four. In this case, P k-1|k satisfies:
Figure PCTCN2020075556-appb-000029
Figure PCTCN2020075556-appb-000029
t k-1时刻采集的观测量z k-1的方差Var(z k-1)满足: The variance Var(z k-1 ) of the observation z k-1 collected at the time t k-1 satisfies:
Figure PCTCN2020075556-appb-000030
Figure PCTCN2020075556-appb-000030
卡尔曼增益值满足:The Kalman gain value satisfies:
Figure PCTCN2020075556-appb-000031
Figure PCTCN2020075556-appb-000031
t k-1时刻的估计观测量满足: The estimated observation at time t k-1 satisfies:
Figure PCTCN2020075556-appb-000032
Figure PCTCN2020075556-appb-000032
t k时刻更新的状态量的期望x k|k满足: k t state the amount of time to update the expectations x k | k satisfy:
Figure PCTCN2020075556-appb-000033
Figure PCTCN2020075556-appb-000033
其中,x k|k可以作为输出的结果,也就是t k时刻对应的状态量。 Among them, x k|k can be used as the output result, that is, the state quantity corresponding to the time t k.
t k时刻更新的状态量的协方差P k|k满足: t k renewed state quantity covariance P k | k satisfying:
P k|k=(I-K' kH k-1)P k-1|k=(1-1*1)*1=0 P k|k =(IK' k H k-1 )P k-1|k =(1-1*1)*1=0
应理解,上述举例是为了帮助本领域技术人员理解本申请实施例,而并非要将本申请实施例限于所例示的具体数值或具体场景。本领域技术人员根据所给出的例子,显然可以进行各种等价的修改或变化,这样的修改或变化也落入本申请实施例的范围内。It should be understood that the foregoing examples are intended to help those skilled in the art understand the embodiments of the present application, and are not intended to limit the embodiments of the present application to the specific numerical values or specific scenarios illustrated. Those skilled in the art can obviously make various equivalent modifications or changes based on the examples given, and such modifications or changes also fall within the scope of the embodiments of the present application.
上文结合图6,详细描述了本申请实施例的目标跟踪方法,下面将结合图7至图8,详细描述本申请的装置实施例。应理解,本申请实施例中的目标跟踪装置可以执行前述本申请实施例的目标跟踪方法,即以下各种产品的具体工作过程,可以参考前述方法实施例中的对应过程。The target tracking method of the embodiment of the present application is described in detail above with reference to FIG. 6, and the device embodiment of the present application will be described in detail below with reference to FIG. 7 to FIG. 8. It should be understood that the target tracking device in the embodiment of the present application can execute the target tracking method of the foregoing embodiment of the present application, that is, the specific working process of the following various products can refer to the corresponding process in the foregoing method embodiment.
图7是本申请实施例提供的目标跟踪装置的示意性框图。应理解,目标跟踪装置1000可以执行图6所示的目标跟踪方法。该目标跟踪装置1000包括:获取单元1010和处理单元1020。Fig. 7 is a schematic block diagram of a target tracking device provided by an embodiment of the present application. It should be understood that the target tracking device 1000 can execute the target tracking method shown in FIG. 6. The target tracking device 1000 includes: an acquiring unit 1010 and a processing unit 1020.
其中,获取单元1010用于获取传感器采集的运动物体的实际观测量。处理单元1020用于在实际观测量的采集时刻早于运动物体的当前状态量的更新时刻的情况下,根据运动物体对应的状态转移模型得到估计观测量;根据实际观测量和估计观测量对当前状态量进行更新,将更新后的状态量作为更新时刻对应的状态量。Wherein, the acquiring unit 1010 is used to acquire the actual observation of the moving object collected by the sensor. The processing unit 1020 is used to obtain the estimated observation according to the state transition model corresponding to the moving object when the acquisition time of the actual observation is earlier than the update time of the current state quantity of the moving object; The state quantity is updated, and the updated state quantity is used as the state quantity corresponding to the update time.
可选地,处理单元1020用于:根据运动物体对应的状态转移模型得到关于采集时刻对应的观测量的似然函数,估计观测量是根据关于采集时刻对应的观测量的似然函数确定的,状态转移模型是根据第一概率密度函数确定的,第一概率密度函数是将当前状态量转移至采集时刻对应的状态量的概率密度函数。Optionally, the processing unit 1020 is configured to: obtain the likelihood function of the observation corresponding to the acquisition time according to the state transition model corresponding to the moving object, and the estimated observation is determined according to the likelihood function of the observation corresponding to the acquisition time, The state transition model is determined according to the first probability density function, and the first probability density function is the probability density function that transfers the current state quantity to the state quantity corresponding to the collection time.
可选地,处理单元1020用于:根据关于采集时刻对应的观测量的似然函数对当前状态量的概率密度函数进行更新;根据更新后的状态量的概率密度函数确定更新时刻对应的状态量。Optionally, the processing unit 1020 is configured to: update the probability density function of the current state quantity according to the likelihood function of the observation corresponding to the acquisition time; determine the state quantity corresponding to the update time according to the probability density function of the updated state quantity .
可选地,关于采集时刻对应的观测量的似然函数满足:Optionally, the likelihood function of the observation corresponding to the acquisition moment satisfies:
g(z k-1|X k')=∫g(z k-1|X)f k-1|k(X|X k')dX g(z k-1 |X k ')=∫g(z k-1 |X)f k-1|k (X|X k ')dX
其中,g(z k-1|X k')表示关于t k-1时刻对应的观测量z k-1的似然函数,X表示t k-1时刻更新的状态量,X k'表示在缺少z k-1的情况下t k时刻更新的状态量,g(z k-1|X)表示在z k-1用于t k-1时刻更新状态量的情况下关于z k-1的似然函数,f k-1|k(X|X k')为第一概率密度函数,第一概率密度函数是将t k时刻更新的状态量X k'转移至到t k-1时刻更新的状态量的概率密度函数。 Among them, g(z k-1 |X k ') represents the likelihood function of the observation z k-1 corresponding to the time t k-1 , X represents the state quantity updated at the time t k-1 , and X k 'represents the missing state quantity update time t K in the case of Z k-1, g (z k-1 | X) Z represents Z k-1 k-1 is used on the time t where k-1 is updated state quantity Likelihood function, f k-1|k (X|X k ') is the first probability density function, the first probability density function is to transfer the state quantity X k ' updated at t k to update at t k-1 The probability density function of the state quantity.
可选地,更新后的状态量的概率密度函数满足:Optionally, the probability density function of the updated state quantity satisfies:
Figure PCTCN2020075556-appb-000034
Figure PCTCN2020075556-appb-000034
其中,f k|k(X k|z k)表示t k时刻更新的状态量X k的概率密度函数,f' k|k(X k'|z k,z k-2)表示在缺少z k-1的情况下t k时刻更新的状态量X k'的概率密度函数,z k表示t k时刻采集的观测量,z k表示k个时刻采集的观测量的集合{z 1,z 2,...,z k-2,z k-1,z k},z k-2表示k-2个时刻采集的观测量的集合{z 1,z 2,...,z k-2}。 Among them, f k|k (X k |z k ) represents the probability density function of the state quantity X k updated at time t k , and f'k|k (X k '|z k ,z k-2 ) represents the absence of z In the case of k-1 , the probability density function of the state quantity X k 'updated at time t k , z k represents the observations collected at time t k , and z k represents the set of observations collected at k time {z 1 , z 2 ,...,z k-2 ,z k-1 ,z k }, z k-2 represents the set of observations collected at k-2 moments {z 1 ,z 2 ,...,z k-2 }.
可选地,处理单元1020用于:根据实际观测量和估计观测量确定更新时刻对应的状态量的期望;将更新时刻对应的状态量的期望作为更新时刻对应的状态量;其中,更新时刻对应的状态量的期望与由卡尔曼增益值有关,卡尔曼增益值与第一协方差、采集时刻的观测矩阵、观测矩阵的协方差和采集时刻的观测量的方差有关,第一协方差指的是由更新时刻转移至采集时刻的状态量的协方差。Optionally, the processing unit 1020 is configured to: determine the expectation of the state quantity corresponding to the update time according to actual observations and estimated observations; use the expectation of the state quantity corresponding to the update time as the state quantity corresponding to the update time; wherein, the update time corresponds to The expectation of the state quantity is related to the Kalman gain value. The Kalman gain value is related to the first covariance, the observation matrix at the acquisition time, the covariance of the observation matrix, and the variance of the observation at the acquisition time. The first covariance refers to It is the covariance of the state amount from the update time to the collection time.
可选地,卡尔曼增益值满足:Optionally, the Kalman gain value satisfies:
Figure PCTCN2020075556-appb-000035
Figure PCTCN2020075556-appb-000035
Var(z k-1)满足: Var(z k-1 ) satisfies:
Figure PCTCN2020075556-appb-000036
Figure PCTCN2020075556-appb-000036
P k-1|k满足: P k-1|k satisfies:
Figure PCTCN2020075556-appb-000037
Figure PCTCN2020075556-appb-000037
其中,P k-1|k表示由更新时刻t k转移到采集时刻t k-1的状态量的协方差,P' k|k表示在缺少t k-1时刻采集的观测量z k-1的情况下t k时刻更新的状态量的协方差,F k-1|k表示从t k时刻转移到t k-1时刻的状态转移矩阵,Q k表示预测矩阵的协方差,H k-1表示t k-1时刻的观测矩阵,Var(z k-1)表示t k-1时刻采集的观测量z k-1的方差,R k-1表示观测矩阵的协方差。 Wherein, P k-1 | k represents t K transferred from the update time to the acquisition time t covariance k-1 state quantity, P 'k | k represents the lack Concept t k-1 time captured measurement Z k-1 In the case of t k , the covariance of the state quantity updated at time t k, F k-1|k represents the state transition matrix from t k to t k-1 , Q k represents the covariance of the prediction matrix, H k-1 Represents the observation matrix at time t k-1 , Var(z k-1 ) represents the variance of the observation z k-1 collected at time t k-1 , and R k-1 represents the covariance of the observation matrix.
可选地,更新时刻对应的状态量的期望满足:Optionally, the expectation of the state quantity corresponding to the update time meets:
Figure PCTCN2020075556-appb-000038
Figure PCTCN2020075556-appb-000038
其中,x k|k表示t k时刻更新的状态量的期望,
Figure PCTCN2020075556-appb-000039
表示t k-1时刻的估计观测量。
Among them, x k|k represents the expectation of the state quantity updated at time t k,
Figure PCTCN2020075556-appb-000039
Represents the estimated observation at time t k-1.
可选地,采集时刻对应的观测量的估计值满足:Optionally, the estimated value of the observation corresponding to the collection time satisfies:
Figure PCTCN2020075556-appb-000040
Figure PCTCN2020075556-appb-000040
其中,
Figure PCTCN2020075556-appb-000041
表示t k-1时刻的估计观测量,H k-1表示t k-1时刻的观测矩阵,F k-1|k表示从t k时刻转移到t k-1时刻的状态转移矩阵,x' k|k表示在缺少t k-1时刻采集的观测量z k-1的情况下的t k时刻更新的状态量的期望。
in,
Figure PCTCN2020075556-appb-000041
Represents t estimated observations k-1 time point, H k-1 represents a t the observation matrix k-1 time, F k-1 | k indicates a transition from t K time to t k-1 time of the state transition matrix, x ' k | k represents the concept of a lack of time t k-1 measurement state quantity acquired in the case where time t k z k-1 update is desired.
可选地,处理单元1020用于:在观测量的采集时刻早于运动物体的当前状态量的更新时刻,且采集时刻与更新时刻之间的时间差小于或等于阈值的情况下,根据运动物体对 应的状态转移模型得到估计观测量。Optionally, the processing unit 1020 is configured to: in the case where the acquisition time of the observation is earlier than the update time of the current state quantity of the moving object, and the time difference between the acquisition time and the update time is less than or equal to a threshold, corresponding to the moving object The state transition model is estimated to be observed.
需要说明的是,上述目标跟踪装置1000以功能单元的形式体现。这里的术语“单元”可以通过软件和/或硬件形式实现,对此不作具体限定。It should be noted that the above-mentioned target tracking device 1000 is embodied in the form of a functional unit. The term "unit" herein can be implemented in the form of software and/or hardware, which is not specifically limited.
例如,“单元”可以是实现上述功能的软件程序、硬件电路或二者结合。所述硬件电路可能包括应用特有集成电路(application specific integrated circuit,ASIC)、电子电路、用于执行一个或多个软件或固件程序的处理器(例如共享处理器、专有处理器或组处理器等)和存储器、合并逻辑电路和/或其它支持所描述的功能的合适组件。For example, a "unit" may be a software program, a hardware circuit, or a combination of the two that realizes the above-mentioned functions. The hardware circuit may include an application specific integrated circuit (ASIC), an electronic circuit, and a processor for executing one or more software or firmware programs (such as a shared processor, a dedicated processor, or a group processor). Etc.) and memory, merged logic circuits and/or other suitable components that support the described functions.
因此,在本申请的实施例中描述的各示例的单元,能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本申请的范围。Therefore, the units of the examples described in the embodiments of the present application can be implemented by electronic hardware or a combination of computer software and electronic hardware. Whether these functions are executed by hardware or software depends on the specific application and design constraint conditions of the technical solution. Professionals and technicians can use different methods for each specific application to implement the described functions, but such implementation should not be considered beyond the scope of this application.
图8是本申请实施例提供的目标跟踪装置的硬件结构示意图。FIG. 8 is a schematic diagram of the hardware structure of the target tracking device provided by an embodiment of the present application.
如图8所示,目标跟踪装置1200(目标跟踪1200具体可以是一种计算机设备)包括存储器1201、处理器1202、通信接口1203以及总线1204。其中,存储器1201、处理器1202、通信接口1203通过总线1204实现彼此之间的通信连接。As shown in FIG. 8, the target tracking apparatus 1200 (the target tracking 1200 may specifically be a computer device) includes a memory 1201, a processor 1202, a communication interface 1203, and a bus 1204. Among them, the memory 1201, the processor 1202, and the communication interface 1203 implement communication connections between each other through the bus 1204.
存储器1201可以是只读存储器(read only memory,ROM),静态存储设备,动态存储设备或者随机存取存储器(random access memory,RAM)。存储器1201可以存储程序,当存储器1201中存储的程序被处理器1202执行时,处理器1202用于执行本申请实施例的目标跟踪方法的各个步骤,例如,执行图6所示的各个步骤。The memory 1201 may be a read only memory (ROM), a static storage device, a dynamic storage device, or a random access memory (RAM). The memory 1201 may store a program. When the program stored in the memory 1201 is executed by the processor 1202, the processor 1202 is configured to execute each step of the target tracking method of the embodiment of the present application, for example, execute each step shown in FIG. 6.
应理解,本申请实施例所示的目标跟踪装置可以是服务器,例如,可以是云端的服务器,或者,也可以是配置于云端的服务器中的芯片。It should be understood that the target tracking device shown in the embodiment of the present application may be a server, for example, it may be a server in the cloud, or may also be a chip configured in a server in the cloud.
处理器1202可以采用通用的中央处理器(central processing unit,CPU),微处理器,应用专用集成电路(application specific integrated circuit,ASIC)或者一个或多个集成电路,用于执行相关程序以实现本申请方法实施例的目标跟踪方法。The processor 1202 may adopt a general central processing unit (CPU), a microprocessor, an application specific integrated circuit (ASIC), or one or more integrated circuits for executing related programs to realize this The target tracking method of the application method embodiment.
处理器1202还可以是一种集成电路芯片,具有信号的处理能力。在实现过程中,本申请的目标跟踪方法的各个步骤可以通过处理器1202中的硬件的集成逻辑电路或者软件形式的指令完成。The processor 1202 may also be an integrated circuit chip with signal processing capability. In the implementation process, each step of the target tracking method of the present application can be completed by an integrated logic circuit of hardware in the processor 1202 or instructions in the form of software.
上述处理器1202还可以是通用处理器、数字信号处理器(digital signal processing,DSP)、专用集成电路(ASIC)、现成可编程门阵列(field programmable gate array,FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本申请实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合本申请实施例所公开的方法的步骤可以直接体现为硬件译码处理器执行完成,或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器,闪存、只读存储器,可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器1201,处理器1202读取存储器1201中的信息,结合其硬件完成本申请实施中图7所示的目标跟踪装置中包括的单元所需执行的功能,或者,执行本申请方法实施例的图6所示的目标跟踪方法。The above-mentioned processor 1202 may also be a general-purpose processor, a digital signal processing (digital signal processing, DSP), an application-specific integrated circuit (ASIC), an off-the-shelf programmable gate array (field programmable gate array, FPGA) or other programmable logic devices, Discrete gates or transistor logic devices, discrete hardware components. The methods, steps, and logical block diagrams disclosed in the embodiments of the present application can be implemented or executed. The general-purpose processor may be a microprocessor or the processor may also be any conventional processor or the like. The steps of the method disclosed in the embodiments of the present application may be directly embodied as being executed and completed by a hardware decoding processor, or executed and completed by a combination of hardware and software modules in the decoding processor. The software module can be located in a mature storage medium in the field, such as random access memory, flash memory, read-only memory, programmable read-only memory, or electrically erasable programmable memory, registers. The storage medium is located in the memory 1201, and the processor 1202 reads the information in the memory 1201, and combines its hardware to complete the functions required by the units included in the target tracking device shown in FIG. 7 in the implementation of this application, or execute the method of this application The target tracking method shown in FIG. 6 of the embodiment.
通信接口1203使用例如但不限于收发器一类的收发装置,来实现目标跟踪装置1200与其他设备或通信网络之间的通信。The communication interface 1203 uses a transceiving device such as but not limited to a transceiver to implement communication between the target tracking device 1200 and other devices or a communication network.
总线1204可包括在目标跟踪装置1200各个部件(例如,存储器1201、处理器1202、通信接口1203)之间传送信息的通路。The bus 1204 may include a path for transferring information between various components of the target tracking device 1200 (for example, the memory 1201, the processor 1202, and the communication interface 1203).
应注意,尽管上述目标跟踪装置1200仅仅示出了存储器、处理器、通信接口,但是在具体实现过程中,本领域的技术人员应当理解,目标跟踪装置1200还可以包括实现正常运行所必须的其他器件。同时,根据具体需要本领域的技术人员应当理解,上述目标跟踪装置1200还可包括实现其他附加功能的硬件器件。It should be noted that although the above-mentioned target tracking device 1200 only shows a memory, a processor, and a communication interface, in the specific implementation process, those skilled in the art should understand that the target tracking device 1200 may also include other necessary for normal operation. Device. At the same time, according to specific needs, those skilled in the art should understand that the above-mentioned target tracking device 1200 may also include hardware devices that implement other additional functions.
此外,本领域的技术人员应当理解,上述目标跟踪装置1200也可仅仅包括实现本申请实施例所必须的器件,而不必包括图8中所示的全部器件。In addition, those skilled in the art should understand that the above-mentioned target tracking device 1200 may also only include the necessary components to implement the embodiments of the present application, and not necessarily include all the components shown in FIG. 8.
本领域普通技术人员可以意识到,结合本文中所公开的实施例描述的各示例的单元及算法步骤,能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本申请的范围。A person of ordinary skill in the art may realize that the units and algorithm steps of the examples described in combination with the embodiments disclosed herein can be implemented by electronic hardware or a combination of computer software and electronic hardware. Whether these functions are executed by hardware or software depends on the specific application and design constraint conditions of the technical solution. Professionals and technicians can use different methods for each specific application to implement the described functions, but such implementation should not be considered beyond the scope of this application.
所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,上述描述的系统、装置和单元的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。Those skilled in the art can clearly understand that, for the convenience and conciseness of description, the specific working process of the system, device and unit described above can refer to the corresponding process in the foregoing method embodiment, which will not be repeated here.
在本申请所提供的几个实施例中,应该理解到,所揭露的系统、装置和方法,可以通过其它的方式实现。例如,以上所描述的装置实施例仅仅是示意性的,例如,所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或单元的间接耦合或通信连接,可以是电性,机械或其它的形式。In the several embodiments provided in this application, it should be understood that the disclosed system, device, and method can be implemented in other ways. For example, the device embodiments described above are only illustrative. For example, the division of the units is only a logical function division, and there may be other divisions in actual implementation, for example, multiple units or components may be combined or It can be integrated into another system, or some features can be ignored or not implemented. In addition, the displayed or discussed mutual coupling or direct coupling or communication connection may be indirect coupling or communication connection through some interfaces, devices or units, and may be in electrical, mechanical or other forms.
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。The units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in one place, or they may be distributed on multiple network units. Some or all of the units may be selected according to actual needs to achieve the objectives of the solutions of the embodiments.
另外,在本申请各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。In addition, the functional units in the various embodiments of the present application may be integrated into one processing unit, or each unit may exist alone physically, or two or more units may be integrated into one unit.
所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者接入网设备等)执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(read-only Memory,ROM)、随机存取存储器(random access memory,RAM)、磁碟或者光盘等各种可以存储程序代码的介质。If the function is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer readable storage medium. Based on this understanding, the technical solution of the present application essentially or the part that contributes to the existing technology or the part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including Several instructions are used to make a computer device (which may be a personal computer, a server, or an access network device, etc.) execute all or part of the steps of the methods described in the various embodiments of the present application. The aforementioned storage media include: U disk, mobile hard disk, read-only memory (read-only Memory, ROM), random access memory (random access memory, RAM), magnetic disk or optical disk and other media that can store program code .
以上所述,仅为本申请的具体实施方式,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本申请的保护范围之内。因此,本申请的保护范围应以所述权利要求的保护范围为准。The above are only specific implementations of this application, but the protection scope of this application is not limited to this. Any person skilled in the art can easily think of changes or substitutions within the technical scope disclosed in this application. Should be covered within the scope of protection of this application. Therefore, the protection scope of this application should be subject to the protection scope of the claims.

Claims (22)

  1. 一种目标跟踪方法,其特征在于,包括:A target tracking method, characterized in that it comprises:
    获取传感器采集的运动物体的实际观测量;Obtain the actual measurement of the moving object collected by the sensor;
    在所述实际观测量的采集时刻早于所述运动物体的当前状态量的更新时刻的情况下,根据所述运动物体对应的状态转移模型得到估计观测量;In the case that the acquisition time of the actual observation is earlier than the update time of the current state quantity of the moving object, obtaining the estimated observation according to the state transition model corresponding to the moving object;
    根据所述实际观测量和所述估计观测量对所述当前状态量进行更新,将更新后的状态量作为所述更新时刻对应的状态量。The current state quantity is updated according to the actual observation and the estimated observation, and the updated state quantity is used as the state quantity corresponding to the update moment.
  2. 如权利要求1所述的方法,其特征在于,所述根据所述运动物体对应的状态转移模型得到估计观测量,包括:The method according to claim 1, wherein the obtaining the estimated observations according to the state transition model corresponding to the moving object comprises:
    根据所述运动物体对应的状态转移模型得到关于所述采集时刻对应的观测量的似然函数,所述估计观测量是根据所述关于所述采集时刻对应的观测量的似然函数确定的,所述状态转移模型是根据第一概率密度函数确定的,所述第一概率密度函数是将所述当前状态量转移至所述采集时刻对应的状态量的概率密度函数。Obtaining the likelihood function of the observation corresponding to the acquisition time according to the state transition model corresponding to the moving object, and the estimated observation is determined according to the likelihood function of the observation corresponding to the acquisition time, The state transition model is determined according to a first probability density function, and the first probability density function is a probability density function that transfers the current state quantity to the state quantity corresponding to the collection time.
  3. 如权利要求2所述的方法,其特征在于,所述根据所述实际观测量和所述估计观测量对所述当前状态量进行更新,将更新后的状态量作为所述更新时刻对应的状态量,包括:The method according to claim 2, wherein the current state quantity is updated according to the actual observation and the estimated observation, and the updated state quantity is used as the state corresponding to the update time Quantity, including:
    根据所述关于所述采集时刻对应的观测量的似然函数对所述当前状态量的概率密度函数进行更新;Updating the probability density function of the current state quantity according to the likelihood function of the observation corresponding to the acquisition time;
    根据更新后的状态量的概率密度函数确定更新时刻对应的状态量。The state quantity corresponding to the update time is determined according to the probability density function of the updated state quantity.
  4. 如权利要求2或3所述的方法,其特征在于,所述关于所述采集时刻对应的观测量的似然函数满足:The method according to claim 2 or 3, wherein the likelihood function of the observation corresponding to the acquisition time satisfies:
    g(z k-1|X k')=∫g(z k-1|X)f k-1|k(X|X k')dX; g(z k-1 |X k ')=∫g(z k-1 |X)f k-1|k (X|X k ')dX;
    其中,g(z k-1|X k')表示关于t k-1时刻对应的观测量z k-1的似然函数,X表示t k-1时刻更新的状态量,X k'表示在缺少z k-1的情况下t k时刻更新的状态量,g(z k-1|X)表示在z k-1用于t k-1时刻更新状态量的情况下关于z k-1的似然函数,f k-1|k(X|X k')为所述第一概率密度函数,所述第一概率密度函数是将t k时刻更新的状态量X k'转移至到t k-1时刻更新的状态量的概率密度函数。 Among them, g(z k-1 |X k ') represents the likelihood function of the observation z k-1 corresponding to the time t k-1 , X represents the state quantity updated at the time t k-1 , and X k 'represents the missing state quantity update time t K in the case of Z k-1, g (z k-1 | X) Z represents Z k-1 k-1 is used on the time t where k-1 is updated state quantity likelihood function, f k-1 | k ( X | X k ') of the first probability density function, the first probability density function is a t k renewed state quantity X k' is transferred to to t k -1 The probability density function of the updated state quantity.
  5. 如权利要求4所述的方法,其特征在于,所述更新后的状态量的概率密度函数满足:The method according to claim 4, wherein the probability density function of the updated state quantity satisfies:
    Figure PCTCN2020075556-appb-100001
    Figure PCTCN2020075556-appb-100001
    其中,f k|k(X k|z k)表示t k时刻更新的状态量X k的概率密度函数,f' k|k(X k'|z k,z k-2)表示在缺少z k-1的情况下t k时刻更新的状态量X k'的概率密度函数,z k表示t k时刻采集的观测量,z k表示k个时刻采集的观测量的集合{z 1,z 2,...,z k-2,z k-1,z k},z k-2表示k-2个时刻采集的观测量的集合{z 1,z 2,...,z k-2}。 Among them, f k|k (X k |z k ) represents the probability density function of the state quantity X k updated at time t k , and f'k|k (X k '|z k ,z k-2 ) represents the absence of z In the case of k-1 , the probability density function of the state quantity X k 'updated at time t k , z k represents the observations collected at time t k , and z k represents the set of observations collected at k time {z 1 , z 2 ,...,z k-2 ,z k-1 ,z k }, z k-2 represents the set of observations collected at k-2 moments {z 1 ,z 2 ,...,z k-2 }.
  6. 如权利要求1所述的方法,其特征在于,所述根据所述实际观测量和所述估计观测量对所述当前状态量进行更新,将更新后的状态量作为所述更新时刻对应的状态量,包 括:The method of claim 1, wherein the current state quantity is updated according to the actual observation and the estimated observation, and the updated state quantity is used as the state corresponding to the update time Quantity, including:
    根据所述实际观测量和所述估计观测量确定更新时刻对应的状态量的期望;Determine the expectation of the state quantity corresponding to the update time according to the actual observation and the estimated observation;
    将所述更新时刻对应的状态量的期望作为所述更新时刻对应的状态量;Taking the expectation of the state quantity corresponding to the update time as the state quantity corresponding to the update time;
    其中,所述更新时刻对应的状态量的期望与由卡尔曼增益值有关,所述卡尔曼增益值与第一协方差、采集时刻的观测矩阵、所述观测矩阵的协方差和采集时刻的观测量的方差有关,所述第一协方差指的是由所述更新时刻转移至所述采集时刻的状态量的协方差。Wherein, the expectation of the state quantity corresponding to the update time is related to the Kalman gain value, the Kalman gain value and the first covariance, the observation matrix at the acquisition time, the covariance of the observation matrix and the observation at the acquisition time The variance of the quantity is related, and the first covariance refers to the covariance of the state quantity transferred from the update time to the collection time.
  7. 如权利要求6所述的方法,其特征在于,所述卡尔曼增益值满足:The method according to claim 6, wherein the Kalman gain value satisfies:
    Figure PCTCN2020075556-appb-100002
    Figure PCTCN2020075556-appb-100002
    Var(z k-1)满足: Var(z k-1 ) satisfies:
    Figure PCTCN2020075556-appb-100003
    Figure PCTCN2020075556-appb-100003
    P k-1|k满足: P k-1|k satisfies:
    Figure PCTCN2020075556-appb-100004
    Figure PCTCN2020075556-appb-100004
    其中,P k-1|k表示由更新时刻t k转移到采集时刻t k-1的状态量的协方差,P' k|k表示在缺少t k-1时刻采集的观测量z k-1的情况下t k时刻更新的状态量的协方差,F k-1|k表示从t k时刻转移到t k-1时刻的状态转移矩阵,Q k表示预测矩阵的协方差,H k-1表示t k-1时刻的观测矩阵,Var(z k-1)表示t k-1时刻采集的观测量z k-1的方差,R k-1表示所述观测矩阵的协方差。 Wherein, P k-1 | k represents t K transferred from the update time to the acquisition time t covariance k-1 state quantity, P 'k | k represents the lack Concept t k-1 time captured measurement Z k-1 In the case of t k , the covariance of the state quantity updated at time t k, F k-1|k represents the state transition matrix from t k to t k-1 , Q k represents the covariance of the prediction matrix, H k-1 Represents the observation matrix at time t k-1 , Var(z k-1 ) represents the variance of the observation z k-1 collected at time t k-1 , and R k-1 represents the covariance of the observation matrix.
  8. 如权利要求7所述的方法,其特征在于,所述更新时刻对应的状态量的期望满足:The method according to claim 7, wherein the expectation of the state quantity corresponding to the update time meets:
    Figure PCTCN2020075556-appb-100005
    Figure PCTCN2020075556-appb-100005
    其中,x k|k表示t k时刻更新的状态量的期望,
    Figure PCTCN2020075556-appb-100006
    表示t k-1时刻的估计观测量。
    Among them, x k|k represents the expectation of the state quantity updated at time t k,
    Figure PCTCN2020075556-appb-100006
    Represents the estimated observation at time t k-1.
  9. 如权利要求6至8中任一项所述的方法,其特征在于,所述估计观测量满足:The method according to any one of claims 6 to 8, wherein the estimated observation quantity satisfies:
    Figure PCTCN2020075556-appb-100007
    Figure PCTCN2020075556-appb-100007
    其中,
    Figure PCTCN2020075556-appb-100008
    表示t k-1时刻的估计观测量,H k-1表示t k-1时刻的观测矩阵,F k-1|k表示从t k时刻转移到t k-1时刻的状态转移矩阵,x' k|k表示在缺少所述t k-1时刻采集的观测量z k-1的情况下的t k时刻更新的状态量的期望。
    in,
    Figure PCTCN2020075556-appb-100008
    Represents t estimated observations k-1 time point, H k-1 represents a t the observation matrix k-1 time, F k-1 | k indicates a transition from t K time to t k-1 time of the state transition matrix, x ' k | k t k represents a time measurement in the case where z k-1 is updated in a desired state quantity of the absence of the concept of time t k-1 collected.
  10. 如权利要求1至9中任一项所述的方法,其特征在于,所述在所述观测量的采集时刻早于所述运动物体的当前状态量的更新时刻的情况下,根据所述运动物体对应的状态转移模型得到估计观测量,包括:The method according to any one of claims 1 to 9, characterized in that, in the case that the acquisition time of the observation is earlier than the update time of the current state quantity of the moving object, the The state transition model corresponding to the object gets estimated observations, including:
    在所述实际观测量的采集时刻早于所述运动物体的当前状态量的更新时刻,且所述采集时刻与所述更新时刻之间的时间差小于或等于阈值的情况下,根据所述运动物体对应的状态转移模型得到估计观测量。In the case that the acquisition time of the actual observation is earlier than the update time of the current state quantity of the moving object, and the time difference between the acquisition time and the update time is less than or equal to a threshold value, according to the moving object The corresponding state transition model gets estimated observations.
  11. 一种目标跟踪装置,其特征在于,包括获取模块和处理模块,其中,A target tracking device, characterized in that it comprises an acquisition module and a processing module, wherein,
    所述获取模块用于:获取传感器采集的运动物体的实际观测量;The acquisition module is used to: acquire the actual observation of the moving object collected by the sensor;
    所述处理模块用于:The processing module is used for:
    在所述实际观测量的采集时刻早于所述运动物体的当前状态量的更新时刻的情况下,根据所述运动物体对应的状态转移模型得到估计观测量;In the case that the acquisition time of the actual observation is earlier than the update time of the current state quantity of the moving object, obtaining the estimated observation according to the state transition model corresponding to the moving object;
    根据所述实际观测量和所述估计观测量对所述当前状态量进行更新,将更新后的状态量作为所述更新时刻对应的状态量。The current state quantity is updated according to the actual observation and the estimated observation, and the updated state quantity is used as the state quantity corresponding to the update moment.
  12. 如权利要求11所述的装置,其特征在于,所述处理模块用于:The device of claim 11, wherein the processing module is configured to:
    根据所述运动物体对应的状态转移模型得到关于所述采集时刻对应的观测量的似然 函数,所述估计观测量是根据所述关于所述采集时刻对应的观测量的似然函数确定的,所述状态转移模型是根据第一概率密度函数确定的,所述第一概率密度函数是将所述当前状态量转移至所述采集时刻对应的状态量的概率密度函数。Obtaining the likelihood function of the observation corresponding to the acquisition time according to the state transition model corresponding to the moving object, and the estimated observation is determined according to the likelihood function of the observation corresponding to the acquisition time, The state transition model is determined according to a first probability density function, and the first probability density function is a probability density function that transfers the current state quantity to the state quantity corresponding to the collection time.
  13. 如权利要求12所述的装置,其特征在于,所述处理模块用于:The device of claim 12, wherein the processing module is configured to:
    根据所述关于所述采集时刻对应的观测量的似然函数对所述当前状态量的概率密度函数进行更新;Updating the probability density function of the current state quantity according to the likelihood function of the observation corresponding to the acquisition time;
    根据更新后的状态量的概率密度函数确定更新时刻对应的状态量。The state quantity corresponding to the update time is determined according to the probability density function of the updated state quantity.
  14. 如权利要求12或13所述的装置,其特征在于,所述关于所述采集时刻对应的观测量的似然函数满足:The device according to claim 12 or 13, wherein the likelihood function of the observation corresponding to the acquisition time satisfies:
    g(z k-1|X k')=∫g(z k-1|X)f k-1|k(X|X k')dX; g(z k-1 |X k ')=∫g(z k-1 |X)f k-1|k (X|X k ')dX;
    其中,g(z k-1|X k')表示关于t k-1时刻对应的观测量z k-1的似然函数,X表示t k-1时刻更新的状态量,X k'表示在缺少z k-1的情况下t k时刻更新的状态量,g(z k-1|X)表示在z k-1用于t k-1时刻更新状态量的情况下关于z k-1的似然函数,f k-1|k(X|X k')为所述第一概率密度函数,所述第一概率密度函数是将t k时刻更新的状态量X k'转移至到t k-1时刻更新的状态量的概率密度函数。 Wherein, g (z k-1 | X k ') denotes Views on T k-1 corresponding to the time measurement z likelihood function k-1, X is t k 1-time updating of the state quantity, X k' represents missing state quantity update time t K in the case of Z k-1, g (z k-1 | X) Z represents Z k-1 k-1 is used on the time t where k-1 is updated state quantity likelihood function, f k-1 | k ( X | X k ') of the first probability density function, the first probability density function is a t k renewed state quantity X k' is transferred to to t k -1 The probability density function of the updated state quantity.
  15. 如权利要求14所述的装置,其特征在于,所述更新后的状态量的概率密度函数满足:The device according to claim 14, wherein the probability density function of the updated state quantity satisfies:
    Figure PCTCN2020075556-appb-100009
    Figure PCTCN2020075556-appb-100009
    其中,f k|k(X k|z k)表示t k时刻更新的状态量X k的概率密度函数,f' k|k(X k'|z k,z k-2)表示在缺少z k-1的情况下t k时刻更新的状态量X k'的概率密度函数,z k表示t k时刻采集的观测量,z k表示k个时刻采集的观测量的集合{z 1,z 2,...,z k-2,z k-1,z k},z k-2表示k-2个时刻采集的观测量的集合{z 1,z 2,...,z k-2}。 Among them, f k|k (X k |z k ) represents the probability density function of the state quantity X k updated at time t k , and f'k|k (X k '|z k ,z k-2 ) represents the absence of z In the case of k-1 , the probability density function of the state quantity X k 'updated at time t k , z k represents the observations collected at time t k , and z k represents the set of observations collected at k time {z 1 , z 2 ,...,z k-2 ,z k-1 ,z k }, z k-2 represents the set of observations collected at k-2 moments {z 1 ,z 2 ,...,z k-2 }.
  16. 如权利要求11所述的装置,其特征在于,所述处理模块用于:The device of claim 11, wherein the processing module is configured to:
    根据所述实际观测量和所述估计观测量确定更新时刻对应的状态量的期望;Determine the expectation of the state quantity corresponding to the update time according to the actual observation and the estimated observation;
    将所述更新时刻对应的状态量的期望作为所述更新时刻对应的状态量;Taking the expectation of the state quantity corresponding to the update time as the state quantity corresponding to the update time;
    其中,所述更新时刻对应的状态量的期望与由卡尔曼增益值有关,所述卡尔曼增益值与第一协方差、采集时刻的观测矩阵、所述观测矩阵的协方差和采集时刻的观测量的方差有关,所述第一协方差指的是由所述更新时刻转移至所述采集时刻的状态量的协方差。Wherein, the expectation of the state quantity corresponding to the update time is related to the Kalman gain value, the Kalman gain value and the first covariance, the observation matrix at the acquisition time, the covariance of the observation matrix and the observation at the acquisition time The variance of the quantity is related, and the first covariance refers to the covariance of the state quantity transferred from the update time to the collection time.
  17. 如权利要求16所述的装置,其特征在于,所述卡尔曼增益值满足:The apparatus according to claim 16, wherein the Kalman gain value satisfies:
    Figure PCTCN2020075556-appb-100010
    Figure PCTCN2020075556-appb-100010
    Var(z k-1)满足: Var(z k-1 ) satisfies:
    Figure PCTCN2020075556-appb-100011
    Figure PCTCN2020075556-appb-100011
    P k-1|k满足: P k-1|k satisfies:
    Figure PCTCN2020075556-appb-100012
    Figure PCTCN2020075556-appb-100012
    其中,P k-1|k表示由更新时刻t k转移到采集时刻t k-1的状态量的协方差,P' k|k表示在缺少t k-1时刻采集的观测量z k-1的情况下t k时刻更新的状态量的协方差,F k-1|k表示从t k时刻转移到t k-1时刻的状态转移矩阵,Q k表示预测矩阵的协方差,H k-1表示t k-1时刻的观测矩阵, Var(z k-1)表示t k-1时刻采集的观测量z k-1的方差,R k-1表示所述观测矩阵的协方差。 Wherein, P k-1 | k represents t K transferred from the update time to the acquisition time t covariance k-1 state quantity, P 'k | k represents the lack Concept t k-1 time captured measurement Z k-1 In the case of t k , the covariance of the state quantity updated at time t k, F k-1|k represents the state transition matrix from t k to t k-1 , Q k represents the covariance of the prediction matrix, H k-1 Represents the observation matrix at time t k-1 , Var(z k-1 ) represents the variance of the observation z k-1 collected at time t k-1 , and R k-1 represents the covariance of the observation matrix.
  18. 如权利要求17所述的装置,其特征在于,所述更新时刻对应的状态量的期望满足:The device according to claim 17, wherein the expectation of the state quantity corresponding to the update time meets:
    Figure PCTCN2020075556-appb-100013
    Figure PCTCN2020075556-appb-100013
    其中,x k|k表示t k时刻更新的状态量的期望,
    Figure PCTCN2020075556-appb-100014
    表示t k-1时刻的估计观测量。
    Among them, x k|k represents the expectation of the state quantity updated at time t k,
    Figure PCTCN2020075556-appb-100014
    Represents the estimated observation at time t k-1.
  19. 如权利要求16至18中任一项所述的装置,其特征在于,所述估计观测量满足:The device according to any one of claims 16 to 18, wherein the estimated observation quantity satisfies:
    Figure PCTCN2020075556-appb-100015
    Figure PCTCN2020075556-appb-100015
    其中,
    Figure PCTCN2020075556-appb-100016
    表示t k-1时刻的估计观测量,H k-1表示t k-1时刻的观测矩阵,F k-1|k表示从t k时刻转移到t k-1时刻的状态转移矩阵,x' k|k表示在缺少所述t k-1时刻采集的观测量z k-1的情况下的t k时刻更新的状态量的期望。
    in,
    Figure PCTCN2020075556-appb-100016
    Represents t estimated observations k-1 time point, H k-1 represents a t the observation matrix k-1 time, F k-1 | k indicates a transition from t K time to t k-1 time of the state transition matrix, x ' k | k t k represents a time measurement in the case where z k-1 is updated in a desired state quantity of the absence of the concept of time t k-1 collected.
  20. 如权利要求11至19中任一项所述的装置,其特征在于,所述处理模块用于:The device according to any one of claims 11 to 19, wherein the processing module is configured to:
    在所述实际观测量的采集时刻早于所述运动物体的当前状态量的更新时刻,且所述采集时刻与所述更新时刻之间的时间差小于或等于阈值的情况下,根据所述运动物体对应的状态转移模型得到估计观测量。In the case that the acquisition time of the actual observation is earlier than the update time of the current state quantity of the moving object, and the time difference between the acquisition time and the update time is less than or equal to a threshold value, according to the moving object The corresponding state transition model gets estimated observations.
  21. 一种目标跟踪装置,其特征在于,包括至少一个处理器和存储器,所述至少一个处理器与所述存储器耦合,用于读取并执行所述存储器中的指令,以执行如权利要求1至10中任一项所述的方法。A target tracking device, characterized by comprising at least one processor and a memory, the at least one processor is coupled to the memory, and is configured to read and execute instructions in the memory to execute as claimed in claims 1 to 10. The method of any one of 10.
  22. 一种计算机可读存储介质,其特征在于,所述计算机可读存储介质上存储有计算机程序,当所述计算机程序在计算机上运行时,使得所述计算机执行如权利要求1至10中任一项所述的方法。A computer-readable storage medium, characterized in that a computer program is stored on the computer-readable storage medium, and when the computer program runs on a computer, the computer executes any one of claims 1 to 10 The method described in the item.
PCT/CN2020/075556 2020-02-17 2020-02-17 Target tracking method and target tracking apparatus WO2021163846A1 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
PCT/CN2020/075556 WO2021163846A1 (en) 2020-02-17 2020-02-17 Target tracking method and target tracking apparatus
CN202080095340.0A CN115039095A (en) 2020-02-17 2020-02-17 Target tracking method and target tracking device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/CN2020/075556 WO2021163846A1 (en) 2020-02-17 2020-02-17 Target tracking method and target tracking apparatus

Publications (1)

Publication Number Publication Date
WO2021163846A1 true WO2021163846A1 (en) 2021-08-26

Family

ID=77390277

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2020/075556 WO2021163846A1 (en) 2020-02-17 2020-02-17 Target tracking method and target tracking apparatus

Country Status (2)

Country Link
CN (1) CN115039095A (en)
WO (1) WO2021163846A1 (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113986787A (en) * 2021-10-09 2022-01-28 许继集团有限公司 Multi-CPU communication data detection method and system
CN115442762A (en) * 2022-08-22 2022-12-06 浙江工业大学 Target tracking method based on distributed consistency filtering of wireless sensor network
CN115131426B (en) * 2022-07-28 2024-03-22 苏州轻棹科技有限公司 Processing method for estimating center point of rear axle of vehicle

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104331902A (en) * 2014-10-11 2015-02-04 深圳超多维光电子有限公司 Target tracking method, target tracking device, 3D display method and 3D display device
CN104994470A (en) * 2015-06-26 2015-10-21 北京航空航天大学 Wireless sensor network collaborative tracking method based on TOA and RSS fusion
CN106021194A (en) * 2016-05-19 2016-10-12 哈尔滨工业大学 Multi-sensor multi-target tracking error estimation method
US20180095103A1 (en) * 2016-09-30 2018-04-05 Panasonic Intellectual Property Management Co., Ltd. State calculation apparatus, state calculation method, and recording medium storing program for moving object

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104331902A (en) * 2014-10-11 2015-02-04 深圳超多维光电子有限公司 Target tracking method, target tracking device, 3D display method and 3D display device
CN104994470A (en) * 2015-06-26 2015-10-21 北京航空航天大学 Wireless sensor network collaborative tracking method based on TOA and RSS fusion
CN106021194A (en) * 2016-05-19 2016-10-12 哈尔滨工业大学 Multi-sensor multi-target tracking error estimation method
US20180095103A1 (en) * 2016-09-30 2018-04-05 Panasonic Intellectual Property Management Co., Ltd. State calculation apparatus, state calculation method, and recording medium storing program for moving object

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113986787A (en) * 2021-10-09 2022-01-28 许继集团有限公司 Multi-CPU communication data detection method and system
CN113986787B (en) * 2021-10-09 2023-06-30 河南源网荷储电气研究院有限公司 Multi-CPU communication data detection method and system
CN115131426B (en) * 2022-07-28 2024-03-22 苏州轻棹科技有限公司 Processing method for estimating center point of rear axle of vehicle
CN115442762A (en) * 2022-08-22 2022-12-06 浙江工业大学 Target tracking method based on distributed consistency filtering of wireless sensor network
CN115442762B (en) * 2022-08-22 2024-05-03 浙江工业大学 Target tracking method based on distributed consistency filtering of wireless sensor network

Also Published As

Publication number Publication date
CN115039095A (en) 2022-09-09

Similar Documents

Publication Publication Date Title
CN109901574B (en) Automatic driving method and device
EP4071661A1 (en) Automatic driving method, related device and computer-readable storage medium
WO2021102955A1 (en) Path planning method for vehicle and path planning apparatus for vehicle
WO2022001773A1 (en) Trajectory prediction method and apparatus
WO2021196879A1 (en) Method and device for recognizing driving behavior of vehicle
EP3835908A1 (en) Automatic driving method, training method and related apparatuses
US20220080972A1 (en) Autonomous lane change method and apparatus, and storage medium
WO2021163846A1 (en) Target tracking method and target tracking apparatus
CN112534483B (en) Method and device for predicting vehicle exit
CN113498529B (en) Target tracking method and device
CN113835421B (en) Method and device for training driving behavior decision model
WO2022062825A1 (en) Vehicle control method, device, and vehicle
US20230048680A1 (en) Method and apparatus for passing through barrier gate crossbar by vehicle
WO2022156309A1 (en) Trajectory prediction method and apparatus, and map
CN111950726A (en) Decision method based on multi-task learning, decision model training method and device
EP4307251A1 (en) Mapping method, vehicle, computer readable storage medium, and chip
WO2022051951A1 (en) Lane line detection method, related device, and computer readable storage medium
CN112810603B (en) Positioning method and related product
WO2022089577A1 (en) Pose determination method and related device thereof
CN113859265B (en) Reminding method and device in driving process
WO2022017307A1 (en) Autonomous driving scenario generation method, apparatus and system
WO2022022284A1 (en) Target object sensing method and apparatus
WO2021159397A1 (en) Vehicle travelable region detection method and detection device
WO2022061725A1 (en) Traffic element observation method and apparatus
CN114556251B (en) Method and device for determining a passable space for a vehicle

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 20919432

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 20919432

Country of ref document: EP

Kind code of ref document: A1