技术领域technical field
本发明属于视觉测量领域,具体涉及利用两个高速摄像机测量开舱点位置的方法及装置。The invention belongs to the field of visual measurement, and in particular relates to a method and a device for measuring the position of a cabin opening point by using two high-speed cameras.
背景技术Background technique
子母弹用于对付集群目标,作为大纵深、大面积的压制兵器,其射击效力主要由子弹落点的分布和散布特性决定。而子弹落点的分布和散布取决于子母弹开舱点的位置及速度,因此有必要测量子母弹开舱点的空间位置。传统方法是用多台光电经纬仪跟踪母弹航迹,获得子母弹开舱点附近的图像,进而推算出子母弹的开舱点位置及速度。该方法因为开舱后目标的突变,致使经纬仪提取目标困难,很难获取精确的开舱点位置和速度,而且光电经纬仪的采样频率较低,导致开舱点的位置和速度精度较低。Submunitions are used to deal with cluster targets. As a large-depth, large-area suppression weapon, its shooting effectiveness is mainly determined by the distribution and dispersion characteristics of the bullets' impact points. The distribution and dispersion of bullet impact points depends on the position and speed of the opening point of the sub-cluster bomb, so it is necessary to measure the spatial position of the opening point of the sub-cluster bomb. The traditional method is to use multiple photoelectric theodolites to track the trajectory of the parent bomb, obtain images near the opening point of the sub-munition, and then calculate the position and speed of the opening point of the sub-munition. Due to the sudden change of the target after the opening of the cabin, it is difficult for the theodolite to extract the target, and it is difficult to obtain the precise position and velocity of the opening point. Moreover, the sampling frequency of the photoelectric theodolite is low, resulting in low accuracy of the position and velocity of the opening point.
发明内容Contents of the invention
为了快速、准确地获取较高精度的子母弹开舱点位置,弥补传统测量方法的不足,本发明提出了利用两个高速摄像机测量开舱点位置的方法及装置。In order to quickly and accurately obtain the position of the opening point of the sub-munitions with higher precision and make up for the shortcomings of the traditional measurement method, the present invention proposes a method and device for measuring the position of the opening point by using two high-speed cameras.
本发明的技术方案为:Technical scheme of the present invention is:
两个高速摄像机测量开舱点位置的方法,包括以下步骤,The method for two high-speed cameras measuring the position of the opening point comprises the following steps,
步骤一,第一高速摄像机的视场为ab、cd间的夹角,第二高速摄像机的视场为ac、bd间的夹角,摄像机靶标及理论开舱点位于四边形abdc中,两个高速摄像机的世界坐标Oj(xj,yj,zj)j=1,2,摄像机靶标的世界坐标O3(x3,y3,z3),三者相距一定距离,且摄像机靶标在第一高速摄像机和第二高速摄像机视场内成清晰的图像;Step 1, the field of view of the first high-speed camera is the angle between ab and cd, the field of view of the second high-speed camera is the angle between ac and bd, the camera target and the theoretical opening point are located in the quadrilateral abdc, the two high-speed The world coordinates of the camera Oj (xj ,yj ,zj )j=1,2, and the world coordinates of the camera target O3 (x3 ,y3 ,z3 ), the three are separated by a certain distance, and the camera target is at A clear image is formed in the fields of view of the first high-speed camera and the second high-speed camera;
步骤二,控制计算机利用网络接口采集第一高速摄像机拍摄的摄像机靶标的图像,并提取特征点,控制计算机利用网络接口采集第二高速摄像机拍摄的摄像机靶标的图像,并提取特征点,利用线性成像模型标定两个高速摄像机参数,利用高速摄像机与摄像机靶标的位置关系解算出摄像机与世界坐标系X轴的方位角α、与水平面的高低角λ;Step 2, the control computer utilizes the network interface to collect the image of the camera target taken by the first high-speed camera, and extracts feature points, and the control computer utilizes the network interface to collect the image of the camera target taken by the second high-speed camera, and extracts the feature points, and utilizes linear imaging The model calibrates the parameters of two high-speed cameras, and uses the positional relationship between the high-speed camera and the camera target to calculate the azimuth α between the camera and the X-axis of the world coordinate system, and the elevation angle λ between the horizontal plane;
其中,高速摄像机的图像坐标系和世界坐标系关系如式①,Among them, the relationship between the image coordinate system of the high-speed camera and the world coordinate system is as follows:
(au,av)为高速摄像机成像器件的像素尺寸;(au , av ) is the pixel size of the high-speed camera imaging device;
(u0,v0)为高速摄像机的光轴与图像平面的交点;(u0 , v0 ) is the intersection of the optical axis of the high-speed camera and the image plane;
R、T分别为摄像机坐标系与世界坐标系的旋转矩阵、平移矩阵;R and T are the rotation matrix and translation matrix of the camera coordinate system and the world coordinate system respectively;
利用高速摄像机与摄像机靶标的位置关系,根据式②解算出高速摄像机与世界坐标系X轴的方位角α、与水平面的高低角λ,Using the positional relationship between the high-speed camera and the camera target, the azimuth α between the high-speed camera and the X-axis of the world coordinate system and the elevation angle λ between the horizontal plane are calculated according to formula ②.
Oj(xj,yj,zj)j=1,2表示两个高速摄像机在世界坐标系中的坐标;Oj (xj ,yj ,zj )j=1,2 represents the coordinates of two high-speed cameras in the world coordinate system;
O3(x3,y3,z3)表示靶标在世界坐标系中的坐标;O3 (x3 , y3 , z3 ) represents the coordinates of the target in the world coordinate system;
步骤三,第一高速摄像机和第二高速摄像机在时统终端的同步脉冲作用下,通过网络接口获取目标的开舱点图像信息,经过图像处理得到开舱点的脱靶量信息;Step 3, the first high-speed camera and the second high-speed camera obtain the image information of the opening point of the target through the network interface under the synchronous pulse action of the time system terminal, and obtain the miss information of the opening point through image processing;
步骤四,控制计算机利用两个高速摄像机的参数及高速摄像机与世界坐标系X轴的方位角α、与水平面的高低角λ信息,再根据公式③计算出开舱点的位置信息(x,y,z),Step 4, the control computer uses the parameters of the two high-speed cameras, the azimuth α between the high-speed camera and the X-axis of the world coordinate system, and the height angle λ of the horizontal plane, and then calculates the position information of the opening point according to the formula ③ (x, y ,z),
Oj(xj,yj,zj)j=1,2表示两个高速摄像机在世界坐标系中的坐标;Oj (xj ,yj ,zj )j=1,2 represents the coordinates of two high-speed cameras in the world coordinate system;
(αj,λj)j=1,2表示两个高速摄像机与世界坐标系X轴的方位角、与水平面的高低角。(αj ,λj )j=1,2 represents the azimuth angle between the two high-speed cameras and the X-axis of the world coordinate system, and the elevation angle between them and the horizontal plane.
利用高速摄像机测量开舱点位置的方法所应用的测量装置,时统终端依靠串行接口连接第一高速摄像机,时统终端依靠串行接口连接第二高速摄像机;控制计算机依靠网络接口连接第一高速摄像机,控制计算机依靠网络接口连接第二高速摄像机;The measurement device used in the method of measuring the position of the opening point with a high-speed camera, the time system terminal is connected to the first high-speed camera by a serial interface, and the time system terminal is connected to the second high-speed camera by a serial interface; the control computer is connected to the first high-speed camera by a network interface. A high-speed camera, the control computer is connected to the second high-speed camera through the network interface;
摄像机的标定靶标,其特征点均能在两个高速摄像机的视场内成清晰的图像;The calibration target of the camera, its feature points can form a clear image in the field of view of the two high-speed cameras;
时统终端,其通过串行接口为第一高速摄像机提供同步脉冲信号和时间信息;其通过串行接口为第二高速摄像机提供同步脉冲信号和时间信息;A time system terminal, which provides a synchronous pulse signal and time information for the first high-speed camera through a serial interface; it provides a synchronous pulse signal and time information for the second high-speed camera through a serial interface;
控制计算机,其通过网络接口接收第一高速摄像机的靶标及子母弹开舱点的图像信息;其通过网络接口接收第二高速摄像机的靶标及子母弹开舱点的图像信息;控制计算机经过数据处理后交会出子母弹开舱点的位置信息。The control computer receives the image information of the target of the first high-speed camera and the opening point of the sub-cluster bomb through the network interface; it receives the image information of the target of the second high-speed camera and the opening point of the sub-cluster bomb through the network interface; the control computer rendezvous after data processing The position information of the opening point of the sub-cluster bomb is displayed.
本发明的有益效果是:利用两个高速摄像机采集靶标的特征点图像(两个高速摄像机及靶标位置已知),解算出摄像机参数、方位角和高低角;在时统终端的作用下,两个高速摄像机同时获取开舱点的图像信息,提取脱靶量信息,交会计算出开舱点的位置信息。此方法简单易行、测量精度较高,可有效降低测量成本。The beneficial effect of the present invention is: use two high-speed cameras to collect the feature point images of the target (the two high-speed cameras and the target position are known), and calculate the camera parameters, azimuth angle and elevation angle; under the action of the time system terminal, the two A high-speed camera acquires the image information of the opening point at the same time, extracts the miss information, and calculates the location information of the opening point by intersection. This method is simple and easy, has high measurement accuracy, and can effectively reduce the measurement cost.
附图说明Description of drawings
图1:本发明两个高速摄像机测量开舱点位置的方法所用装置的示意图。Fig. 1: the schematic diagram of the device used in the method for measuring the position of the opening point with two high-speed cameras of the present invention.
图2:本发明两个高速摄像机测量开舱点位置的方法的工作流程图。Fig. 2: The working flowchart of the method for measuring the position of the opening point by two high-speed cameras of the present invention.
图3:本发明中高速摄像机的线性成像模型示意图。Figure 3: Schematic diagram of the linear imaging model of the high-speed camera in the present invention.
具体实施方式Detailed ways
下面结合附图对本发明做进一步详细说明。The present invention will be described in further detail below in conjunction with the accompanying drawings.
如图1所示,两个高速摄像机测量开舱点位置的方法所应用的装置,包括时统终端3依靠串行接口5连接第一高速摄像机1,时统终端3依靠串行接口6连接第二高速摄像机2;控制计算机4依靠网络接口7连接第一高速摄像机1,控制计算机4依靠网络接口8连接第二高速摄像机2。As shown in Figure 1, the device that two high-speed cameras are used in the method for measuring the position of the opening point comprises that the time system terminal 3 is connected to the first high-speed camera 1 by means of the serial interface 5, and the time system terminal 3 is connected to the first high-speed camera by means of the serial interface 6. Two high-speed cameras 2; the control computer 4 is connected to the first high-speed camera 1 by means of the network interface 7, and the control computer 4 is connected to the second high-speed camera 2 by means of the network interface 8.
摄像机靶标9特征点均能在第一高速摄像机1和第二高速摄像机2的交叉视场内成清晰的图像。All the feature points of the camera target 9 can form a clear image in the intersection field of view of the first high-speed camera 1 and the second high-speed camera 2 .
时统终端3通过串行接口5为第一高速摄像机1提供同步脉冲信号和时间信息;时统终端3通过串行接口6为第二高速摄像机2提供同步脉冲信号和时间信息。Time system terminal 3 provides synchronization pulse signal and time information for first high-speed camera 1 through serial interface 5 ; time system terminal 3 provides synchronization pulse signal and time information for second high-speed camera 2 through serial interface 6 .
控制计算机4通过网络接口7接收第一高速摄像机1的靶标及子母弹开舱点的图像信息;控制计算机4通过网络接口8接收第二高速摄像机2的靶标及子母弹开舱点的图像信息;控制计算机4经过数据处理后交会出子母弹开舱点的位置信息。The control computer 4 receives the target of the first high-speed camera 1 and the image information of the opening point of the sub-cluster bomb through the network interface 7; the control computer 4 receives the target of the second high-speed camera 2 and the image information of the opening point of the sub-cluster bomb through the network interface 8; the control computer 4 After data processing, the location information of the opening point of the sub-munition is intersected.
如图2所示,两个高速摄像机测量开舱点位置的方法,实现步骤如下:As shown in Figure 2, the method for two high-speed cameras to measure the position of the opening point is implemented as follows:
步骤一,按图1所示,第一高速摄像机1的视场为ab、cd间的夹角,第二高速摄像机2的视场为ac、bd间的夹角,摄像机靶标9及理论开舱点位于四边形abdc中,在世界坐标系中,两个高速摄像机位于坐标Oj(xj,yj,zj)j=1,2,靶标位于坐标O3(x3,y3,z3),三者相距一定距离,且靶标在第一高速摄像机1和第二高速摄像机2视场内成清晰的图像。Step 1, as shown in Figure 1, the field of view of the first high-speed camera 1 is the angle between ab and cd, the field of view of the second high-speed camera 2 is the angle between ac and bd, the camera target 9 and the theoretical opening The point is located in the quadrilateral abdc. In the world coordinate system, the two high-speed cameras are located at the coordinates Oj (xj ,yj ,zj )j=1,2, and the target is located at the coordinates O3 (x3 ,y3 ,z3 ), the three are separated by a certain distance, and the target forms a clear image in the field of view of the first high-speed camera 1 and the second high-speed camera 2 .
步骤二,控制计算机4利用网络接口7接收第一高速摄像机1的摄像机靶标9的图像,并提取特征点,控制计算机4利用网络接口8接收第二高速摄像机2的摄像机靶标9的图像,并提取特征点。由图3可以得到摄像机的图像坐标系和世界坐标系关系如式①,Step 2, the control computer 4 utilizes the network interface 7 to receive the image of the camera target 9 of the first high-speed camera 1, and extracts feature points, and the control computer 4 utilizes the network interface 8 to receive the image of the camera target 9 of the second high-speed camera 2, and extracts Feature points. From Figure 3, the relationship between the image coordinate system of the camera and the world coordinate system can be obtained as formula ①,
(au,av)为高速摄像机成像器件的像素尺寸;(au , av ) is the pixel size of the high-speed camera imaging device;
(u0,v0)为高速摄像机的光轴与图像平面的交点;(u0 , v0 ) is the intersection of the optical axis of the high-speed camera and the image plane;
R、T分别为摄像机坐标系与世界坐标系的旋转矩阵、平移矩阵。R and T are the rotation matrix and translation matrix of the camera coordinate system and the world coordinate system, respectively.
利用高速摄像机与摄像机靶标的位置关系解算出摄像机与世界坐标系X轴的方位角α、与水平面的高低角λ,Use the positional relationship between the high-speed camera and the camera target to calculate the azimuth α between the camera and the X-axis of the world coordinate system, and the height angle λ between the horizontal plane,
Oj(xj,yj,zj)j=1,2表示两个高速摄像机在世界坐标系中的坐标;Oj (xj ,yj ,zj )j=1,2 represents the coordinates of two high-speed cameras in the world coordinate system;
O3(x3,y3,z3)表示摄像机靶标9在世界坐标系中的坐标。O3 (x3 , y3 , z3 ) represents the coordinates of the camera target 9 in the world coordinate system.
步骤三,第一高速摄像机1和第二高速摄像机2在时统终端3的同步脉冲作用下,通过网络接口获取目标的开舱点图像信息,经过图像处理得到开舱点的脱靶量信息。Step 3, the first high-speed camera 1 and the second high-speed camera 2 obtain the image information of the opening point of the target through the network interface under the action of the synchronous pulse of the time system terminal 3, and obtain the information on the amount of misses at the opening point through image processing.
步骤四,控制计算机4利用两个高速摄像机的参数及高速摄像机与世界坐标系X轴的方位角α、与水平面的高低角λ信息,根据公式③计算出开舱点的位置信息(x,y,z),Step 4, the control computer 4 utilizes the parameters of the two high-speed cameras and the azimuth α between the high-speed camera and the X-axis of the world coordinate system, and the height angle λ information of the horizontal plane, and calculates the position information (x, y) of the opening point according to the formula 3. ,z),
Oj(xj,yj,zj)j=1,2表示两个高速摄像机在世界坐标系中的坐标;Oj (xj ,yj ,zj )j=1,2 represents the coordinates of two high-speed cameras in the world coordinate system;
(αj,λj)j=1,2表示两个高速摄像机与世界坐标系X轴的方位角、与水平面的高低角。(αj ,λj )j=1,2 represents the azimuth angle between the two high-speed cameras and the X-axis of the world coordinate system, and the elevation angle between them and the horizontal plane.
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410122961.XACN103925909B (en) | 2014-03-28 | 2014-03-28 | Two high-speed camera measurements are opened the cabin the method for a position and device |
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410122961.XACN103925909B (en) | 2014-03-28 | 2014-03-28 | Two high-speed camera measurements are opened the cabin the method for a position and device |
| Publication Number | Publication Date |
|---|---|
| CN103925909Atrue CN103925909A (en) | 2014-07-16 |
| CN103925909B CN103925909B (en) | 2016-06-15 |
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410122961.XAExpired - Fee RelatedCN103925909B (en) | 2014-03-28 | 2014-03-28 | Two high-speed camera measurements are opened the cabin the method for a position and device |
| Country | Link |
|---|---|
| CN (1) | CN103925909B (en) |
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106643306A (en)* | 2016-12-30 | 2017-05-10 | 中国科学院长春光学精密机械与物理研究所 | High-speed imaging method and system for light screen target trajectory measuring system |
| CN107976146A (en)* | 2017-11-01 | 2018-05-01 | 中国船舶重工集团公司第七〇九研究所 | The self-calibrating method and measuring method of a kind of linear array CCD camera |
| CN117516481A (en)* | 2024-01-08 | 2024-02-06 | 北京奥博泰科技有限公司 | Dynamic image intersection measuring method and device |
| CN120232400A (en)* | 2025-05-30 | 2025-07-01 | 中国人民解放军95859部队 | A surface-to-surface intersection measurement method for terminal trajectory parameters based on linear motion constraints |
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000028332A (en)* | 1998-07-07 | 2000-01-28 | Mitsubishi Heavy Ind Ltd | Three-dimensional measuring device and method therefor |
| CN101970985A (en)* | 2008-02-29 | 2011-02-09 | 特林布尔公司 | Determining coordinates of a target in relation to a survey instrument having at least two cameras |
| CN103512559A (en)* | 2013-10-08 | 2014-01-15 | 北京理工大学 | Shot monocular video pose measurement method and target pattern |
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000028332A (en)* | 1998-07-07 | 2000-01-28 | Mitsubishi Heavy Ind Ltd | Three-dimensional measuring device and method therefor |
| CN101970985A (en)* | 2008-02-29 | 2011-02-09 | 特林布尔公司 | Determining coordinates of a target in relation to a survey instrument having at least two cameras |
| CN103512559A (en)* | 2013-10-08 | 2014-01-15 | 北京理工大学 | Shot monocular video pose measurement method and target pattern |
| Title |
|---|
| 吴能伟: "光电经纬仪目标异面交会测量的误差分析", 《中国仪器仪表》, no. 6, 30 June 2006 (2006-06-30)* |
| 赵立荣: "光电经纬仪单站空间余弦即多站面面交汇的飞机姿态测量", 《光学 精密工程》, vol. 17, no. 11, 30 November 2009 (2009-11-30)* |
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106643306A (en)* | 2016-12-30 | 2017-05-10 | 中国科学院长春光学精密机械与物理研究所 | High-speed imaging method and system for light screen target trajectory measuring system |
| CN106643306B (en)* | 2016-12-30 | 2018-07-06 | 中国科学院长春光学精密机械与物理研究所 | A kind of high speed imaging method and its system for light curtain target trajectory measurement system |
| CN107976146A (en)* | 2017-11-01 | 2018-05-01 | 中国船舶重工集团公司第七〇九研究所 | The self-calibrating method and measuring method of a kind of linear array CCD camera |
| CN107976146B (en)* | 2017-11-01 | 2019-12-10 | 中国船舶重工集团公司第七一九研究所 | Self-calibration method and measurement method of linear array CCD camera |
| CN117516481A (en)* | 2024-01-08 | 2024-02-06 | 北京奥博泰科技有限公司 | Dynamic image intersection measuring method and device |
| CN117516481B (en)* | 2024-01-08 | 2024-04-16 | 北京奥博泰科技有限公司 | Dynamic image intersection measuring method and device |
| CN120232400A (en)* | 2025-05-30 | 2025-07-01 | 中国人民解放军95859部队 | A surface-to-surface intersection measurement method for terminal trajectory parameters based on linear motion constraints |
| Publication number | Publication date |
|---|---|
| CN103925909B (en) | 2016-06-15 |
| Publication | Publication Date | Title |
|---|---|---|
| CN109272532B (en) | Model pose calculation method based on binocular vision | |
| CN103925909B (en) | Two high-speed camera measurements are opened the cabin the method for a position and device | |
| KR101898782B1 (en) | Apparatus for tracking object | |
| CN105785064B (en) | A kind of the Projectile velocity measurements device and its measurement method of any direction incidence | |
| CN103223236B (en) | Intelligent evaluation system for table tennis training machine | |
| CN108036776A (en) | It is a kind of to utilize radiographic measurement velocity of shot and the method for posture | |
| CN105373011B (en) | Detect the real-time emulation system and computer of electro-optical tracking device | |
| CN104330075B (en) | Rasterizing polar coordinate system object localization method | |
| CN103364579A (en) | Method for predicting ping-pong spin angle velocity of ping-pong robot | |
| CN103575167B (en) | Trajectory correction method for civil interceptor missiles | |
| CN105300345B (en) | Multi-target Tracking Method for Photoelectric Theodolite | |
| CN104482924A (en) | Revolution body object pose vision measurement method | |
| CN101839677A (en) | Acousto-optic automatic target reporting system | |
| CN117516481B (en) | Dynamic image intersection measuring method and device | |
| CN106643348A (en) | Semi-physical simulating device for guided missile | |
| CN110787433A (en) | Golf ball sensing system for golf simulator | |
| CN110866954B (en) | High-precision attitude measurement method of missile target under length constraint | |
| CN106156745A (en) | Pedestrian crossing traffic feature extraction method and device based on space-time track | |
| CN112710196A (en) | Target design method based on structural equivalence | |
| CN107990789B (en) | Method for determining position of projectile based on vibration signal | |
| CN102213567A (en) | Intelligent shooting scoring system | |
| CN108984808A (en) | Fuze warhead matching emulation intersection fragmentation point of impact calculation method of parameters | |
| Zhou et al. | Research on a measurement method for spatio-temporal distribution of fragments based on a high-speed camera network | |
| CN106225668A (en) | Wind-tunnel missile high speed pose measuring methods based on many refraction models | |
| CN106932605A (en) | A kind of gun bullet movement velocity measuring system |
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee | Granted publication date:20160615 Termination date:20180328 |