CN108896017A - 一种弹丸近炸破片群位置参数测量与计算方法 - Google Patents
一种弹丸近炸破片群位置参数测量与计算方法 Download PDFInfo
- Publication number
- CN108896017A CN108896017A CN201810438708.3A CN201810438708A CN108896017A CN 108896017 A CN108896017 A CN 108896017A CN 201810438708 A CN201810438708 A CN 201810438708A CN 108896017 A CN108896017 A CN 108896017A
- Authority
- CN
- China
- Prior art keywords
- light
- fragmentation
- field camera
- target
- sin
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000012634 fragment Substances 0.000 title claims abstract description 63
- 238000004364 calculation method Methods 0.000 title claims abstract description 29
- 238000005259 measurement Methods 0.000 title claims abstract description 28
- 238000013467 fragmentation Methods 0.000 claims abstract description 98
- 238000006062 fragmentation reaction Methods 0.000 claims abstract description 98
- 230000003287 optical effect Effects 0.000 claims abstract description 16
- 238000005516 engineering process Methods 0.000 claims abstract description 6
- 239000011159 matrix material Substances 0.000 claims description 30
- 238000012360 testing method Methods 0.000 claims description 14
- 238000001514 detection method Methods 0.000 claims description 12
- 238000003384 imaging method Methods 0.000 claims description 12
- 230000000007 visual effect Effects 0.000 claims description 12
- 238000000034 method Methods 0.000 claims description 8
- 230000000903 blocking effect Effects 0.000 claims description 4
- 230000000694 effects Effects 0.000 claims description 3
- 238000009499 grossing Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 239000000571 coke Substances 0.000 claims description 2
- 230000001360 synchronised effect Effects 0.000 claims description 2
- 239000004744 fabric Substances 0.000 claims 1
- 238000000926 separation method Methods 0.000 claims 1
- 238000004458 analytical method Methods 0.000 abstract description 5
- 238000001444 catalytic combustion detection Methods 0.000 description 12
- 238000010586 diagram Methods 0.000 description 7
- 238000000691 measurement method Methods 0.000 description 4
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000005622 photoelectricity Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C11/00—Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C11/00—Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
- G01C11/04—Interpretation of pictures
- G01C11/30—Interpretation of pictures by triangulation
- G01C11/34—Aerial triangulation
Landscapes
- Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Image Processing (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明提供一种弹丸近炸破片群位置参数测量与计算方法，其具体包括根据破片群的位置布置三个光场相机，建立近炸破片群散布解算模型；利用数字重聚焦技术，对遮挡、重合的多个破片目标进行分辨识别，获取破片目标图像坐标；结合破片目标经过数字重聚焦后的图像坐标，确定空间破片与各个相机光轴之间的夹角；基于双目视觉原理，对两两交汇的相机进行分析，根据空间几何关系计算出破片在预定平面的空间位置坐标参数等步骤。本发明给出比较形象的弹丸近炸破片群位置参数测量与计算方法，为新型目标毁伤计算与分析提供有效的理论依据。
Description
技术领域
本发明涉及光电测试技术领域，尤其涉及一种弹丸近炸破片群位置参数测量和计算方法。
背景技术
在实际作战环境中，破片是武器对目标的主要毁伤手段，破片能否有效毁伤目标与破片打击目标时的空间位置密切相关，因此，弹丸近炸破片的散布参数是目标毁伤计算的重要参数和分析毁伤评价体系的有效理论依据。
现有的研究也提出一些目标位置参数的装置与测量方法，具体包括光幕靶测试系统、声靶测试系统以及多CCD交汇测试系统等；光幕靶一般采用一体靶架式结构，采用阵列式光电发射接收原理对目标位置进行测量，当多目标经过探测区时，光幕靶很难进行匹配识别，另外，其装置受结构影响，靶面不能设置较大，因此难以对破片群散布参数进行分辨测试；多CCD交汇测试系统利用多个线阵CCD构成多个靶面，利用时空匹配的方法对多目标散布参数进行测试，但当测试多个破片目标时，多个线阵CCD无法分辨识别重合的破片目标，容易漏测，无法满足高精度测试散布参数的要求；声靶采用声阵列传感器对目标位置进行测量计算，当多个目标同时着靶时，多个目标的声波同时到达传感器，此时声靶无法从信号中分辨多目标的信息，因此同样无法满足破片群散布参数测试的要求；综上所述，现有的专利与文献中提供的目标位置测量装置与计算方法均无法分辨遮挡、重合的多个碎片，迫切需要一种满足近炸破片群散布参数测试要求的测量与计算方法。
发明内容
本发明提供一种弹丸近炸破片群位置参数测量和计算方法，用以解决现有技术中存在的技术问题。
本发明一方面提供一种弹丸近炸破片群位置参数测量与计算方法，其具体包括以下步骤：
步骤1：根据破片群的位置，布置三个光场相机，根据三个光场相机布置参数，建立基于三个面阵CCD相机的近炸空间破片散布解算模型；
步骤2：利用数字重聚焦技术，对遮挡、重合的多个破片目标进行分辨识别，获取破片目标的图像坐标；
步骤3：结合破片目标经过数字重聚焦后的图像坐标，确定空间破片与各个光场相机的光轴之间的夹角；
步骤4：当破片目标经过探测区域时，利用重聚焦解算目标图像坐标参数模型；基于建立的近炸空间破片散布解算模型，依据双目视觉原理，对两两交汇的光场相机进行分析，根据空间几何关系计算出破片在主坐标系下与从坐标系下的空间位置坐标参数；依据三个坐标系下的空间位置参数，根据同一破片目标的方位角与俯仰角在图像序列中不变的原理进行破片群匹配；
步骤5：取三个坐标系下的空间位置坐标平均值，获取破片目标的空间位置。
优选地，在步骤1中，包括以下步骤：
步骤S11，采用三个光场相机，且三个光场相机都具有满足测试要求的高帧频，将其布置到三角形的三个顶点，且使三个光场相机的视场交汇于探测区域；
步骤S12，对三个光场相机中两两一组进行校平；
步骤S13，测量每一组光场相机校平后的角度值αa1、αa3、αb1、αb2、αa2、αb3；
步骤S14，测量三个光场相机光学系统的俯仰角数值ε1、ε2、ε3；
步骤S15，测量每两个光场相机的光学镜头之间的距离d1、d2、d3；依据建立的破片空间位置解算模型，在炮口处设置外部触发源，当触发源探测火光时，三个光场相机同步触发采集图像,；
步骤S16，建立主坐标系，具体地，以第一光场相机和第三光场相机的中心连线S1S3为X轴，S1为原点，构成主坐标系XOY；
步骤S17，建立两个从坐标系，其中，以第一光场相机和第二光场相机的中心连线S1S2为X'轴，S1为原点O'，构成第一从坐标系；以第二光场相机和第三光场相机的中心连线S2S3为X”轴，S2为原点O”，构成第二从坐标系；其中，主坐标系与第一从坐标系的夹角为δ1，主坐标系与第二从坐标系的夹角为δ3，第一从坐标系与第二从坐标系的夹角为δ2，其中，δ1＝αb1-αa1，δ2＝αb2-αa2，δ3＝αb3-αa3。
优选地，在步骤2中，具体包括以下步骤：
步骤S21，对爆炸区域的光场信息利用光场相机内部的微透镜阵列进行处理，每一光场相机内的CCD成像探测面置于微透镜阵列的一倍焦距f处，设光线进入微透镜的位置坐标(ui,vi)，位置坐标矩阵U＝[u1,u2,…un]T，V＝[v1,v2,…vn]T，其中，i＝1,2,···,n；经过图像处理技术，获取空间中每一破片在CCD成像探测面上的图像坐标位置(xi′,yi′)，设图像坐标矩阵X′＝[x1′,x′2,…x′n]T，Y′＝[y1′,y′2,…y′n]T；当光场相机采集图像中出现多个破片目标重合或者遮挡的情况时，通过改变微透镜阵列参数ai，使重合遮挡破片目标区域的成像景深变浅，在还原被遮挡破片时，能够将前方的离焦遮挡破片严重虚化，实现透视效果；
步骤S22，将微透镜阵列平面视为u-v平面，探测器平面视为x-y平面，则重聚焦平面距微透镜平面的距离为aif，对遮挡、重合的破片目标分辨识别后，使得空间破片通过重新聚焦在CCD成像探测面上形成一个清晰的像，从而得到破片目标经过重聚焦后的图像坐标(xic,yic)，则Xc＝[x1c,x2c,…xnc]T，Yc＝[y1c,y2c,…ync]T，且另其中，i＝1,2,···,n，E为单位矩阵，且满足如下关系：
优选地，在步骤3中，具体包括以下步骤：
步骤S31，根据步骤2中的数字重聚焦方法获取每个破片目标的图像坐标(xic,yic)，计算每个破片与各个光场相机光轴之间的夹角βmi，m＝1,2,3，i＝1,2,···,n则有：
优选地，在步骤4中，包括以下步骤：
步骤S41，当破片群经过探测区域时，三个光场相机同步对破片群图像进行采集，依据采集的图像，调整微透镜阵列参数矩阵A进行分辨识别破片目标，解算图像坐标参数模型(Xc,Yc)；假设每台光场相机的焦距为fl，根据破片经过重聚焦后的图像坐标(xi′,yi′)，计算每个破片与各个光场相机光轴之间的夹角βmi，其中，i＝1,2,···,n，l＝1,2,3，m＝1,2,3。
步骤S42，分别对两两交汇的光场相机进行分析，根据几何关系计算第k幅光场相机采集的图像中第i个破片目标的空间位置(xki,yki,zki)，对应的空间位置矩阵(Xk,Yk,Zk)，其中，i＝1,2,…,n，k＝1,2,…，具体分为如下情况：
利用第一光场相机与第三光场相机视场交汇得到破片群的空间位置矩阵(X1,Y1,Z1)，具体如下：
X1＝d1·Q1 TW1E1-f1·R1 TM1
Y1＝d1·Q1 TW1E1
Z1＝d1·C1 TB1
其中，
Q1＝[sin(ε3-β31),sin(ε3-β32),…，sin(ε3-β3n)]T
W1＝[cos(ε1-β11),cos(ε1-β12),…，cos(ε1-β1n)]T
B1＝[cot(ε3+β31),cot(ε3+β32),…，cot(ε3+β3n)]T
R1＝[sin(β11),sin(β12),…，sin(β1n)]T
利用第一光场相机与第二光场相机视场交汇得到破片群的空间位置矩阵(X2,Y2,Z2)，具体如下：
X2＝(d2·Q2 TW2E2-f2·R1 TC2)·cosδ1-(d2·Q2 TE2C2·sinδ1)
Y2＝(d2·Q2 TW2E2-f2·R1 TC2)·sinδ1-(d2·Q2 TE2C2)·cosδ1
Z2＝d2·C1 TN1
其中，
N1＝[cot(ε2+β21),cot(ε2+β22),…，cot(ε2+β2n)]T
Q2＝[sin(ε2+β21),sin(ε2+β22),…，sin(ε2+β2n)]T
W2＝[cos(δ1i-ε1+β11),cos(δ1i-ε1+β12),…，cos(δ1i-ε1+β1n)]T
M2＝[sin(δ1-ε2+β11),sin(δ1-ε2+β12),…，sin(δ1-ε2+β1n)]T
利用第二光场相机与第三光场相机视场交汇得到破片群的空间位置矩阵(X3,Y3,Z3)，具体如下：
X3＝(d3·Q3 TW3E3-f3·R3 TZ3)·cosδ3-(d3·Q3 TZ3E3)·sinδ3+d2·cosδ1
Y3＝(d3·Q3 TW3E3-f3·R3 TZ3)·sinδ3-(d3·Q3 T·Z3E3)·cosδ3+d2·sinδ1
Z3＝d2C1 TN1
其中，
Q3＝[sin(δ3-ε3+β31),sin(δ3-ε3+β32),…，sin(δ3-ε3+β3n)]T
W3＝[cos(δ2-ε2-β21),cos(δ2-ε2-β22),…，cos(δ2-ε2-β2n)]T
R3＝[sin(β21),sin(β22),…，sin(β2n)]T
步骤S43，依据解算的空间位置矩阵，用相邻两幅图像待配准的破片目标群位置坐标进行作差，根据同一破片目标的方位角与俯仰角在图像序列中不变的原理进行匹配，若坐标差成比例，则待配准目标为同一目标。
优选地，在步骤5中，包括以下步骤：
步骤S51，取三个坐标系的空间位置平均值，从而得到破片目标的空间位置(xi,yi,zi)，其中，破片空间位置坐标矩阵X＝[x1,x2,…xn]T，Y＝[y1,y2,…yn]T，Z＝[z1,z2,…zn]T，且满足如下关系：
本发明利用了光场相机数字重聚焦原理，建立了破片群位置参数测量与计算模型。本发明可以分辨重合、遮挡的多个破片目标，给出了比较形象的基于光场相机的弹丸近炸破片群位置参数测量与计算方法，为新型目标毁伤计算提供科学依据。
附图说明
图1是本发明涉及的弹丸近炸破片群位置参数测量与计算方法的流程图；
图2是本发明涉及的弹丸近炸破片群位置参数测量与计算方法的空间几何结构示意图；
图3是本发明涉及的弹丸近炸破片群位置参数测量与计算方法中基于三相机布置参数建立解算模型的框图；
图4是本发明涉及的弹丸近炸破片群位置参数测量与计算方法中微透镜阵列重聚焦图像示意图；
图5是本发明涉及的弹丸近炸破片群位置参数测量与计算方法中破片目标飞行方向在相机成像几何关系示意图；
图6是本发明涉及的弹丸近炸破片群位置参数测量与计算方法中破片目标OZ轴方向坐标解算示意图。
图7是本发明涉及的弹丸近炸破片群位置参数测量与计算方法中主坐标系破片坐标解算示意图；
图8是本发明涉及的弹丸近炸破片群位置参数测量与计算方法中第一从坐标系破片坐标解算示意图；
图9是本发明涉及的弹丸近炸破片群位置参数测量与计算方法中第二从坐标系破片坐标解算示意图；
其中：
1、第一光场相机；2、第二光场相机；3、第三光场相机；4、破片群；
具体实施方式
本实施例涉及一种弹丸近炸破片群位置参数测量与计算方法，该弹丸近炸破片群位置参数测量与计算方法主要应用于多目标位置参数测试，可以识别分辨重合遮挡的破片目标，为目标毁伤计算和分析毁伤评价体系提供有效的理论依据。
本实施例所涉及的一种弹丸近炸破片群位置参数测量与计算方法，如图1-2所示，其包括以下步骤：
步骤1：针对破片群4的位置，布置三个光场相机，根据三个光场相机的布置参数，建立基于三个面阵CCD的近炸空间破片散布解算模型；
步骤2：利用数字重聚焦技术，对遮挡、重合的多个破片目标进行分辨识别，获取破片目标的图像坐标；
步骤3：结合破片目标经过数字重聚焦后的图像坐标，确定空间破片与各个光场相机的光轴之间的夹角；
步骤4：当破片目标经过探测区域时，利用重聚焦技术解算破片目标的图像坐标参数模型；基于建立的近炸空间破片散布解算模型，依据双目视觉原理，对两两交汇的相机进行分析，根据空间几何关系计算出破片在主坐标系下与从坐标系下的空间位置坐标参数；依据三个坐标系下的空间位置参数，根据同一破片目标的方位角与俯仰角在图像序列中不变的原理进行破片群匹配；
步骤5：为了提高破片群4的测量精度，取三个坐标系下的空间位置坐标平均值，获取破片目标的空间位置。
如图2和3所示，在上述步骤1中，具体地，包括以下步骤：
步骤S11，采用三个光场相机，分别是第一光场相机1、第二光场相机2和第三光场相机3，其中，三个光场相机都具有满足测试要求的高帧频，将三个光场相机布置到三角形的三个顶点，且使三个相机的视场交汇于探测区域；
步骤S12，对三个光场相机两两一组进行校平；
步骤S13，测量每一组光场相机校平后的角度值αa1、αa3、αb1、αb2、αa2、αb3；
步骤S14，测量三个光场相机光学系统的俯仰角数值ε1、ε2、ε3；
步骤S15，测量每两个光场相机的光学镜头之间的距离d1、d2、d3；依据建立的破片空间位置解算模型，在发出破片目标的炮口处设置外部触发源，当触发源探测火光时，三个光场相机同步触发采集图像；
步骤S16，建立主坐标系，具体地，以第一光场相机1和第三光场相机3的中心连线S1S3为X轴，S1为原点，构成主坐标系XOY；
步骤S17，建立两个从坐标系，其中，以第一光场相机1和第二光场相机2的中心连线S1S2为X'轴，S1为原点O'，构成第一从坐标系；以第二光场相机2和第三光场相机3的中心连线S2S3为X”轴，S2为原点O”，构成第二从坐标系；其中，主坐标系与第一从坐标系的夹角为δ1，主坐标系与第二从坐标系的夹角为δ3，第一从坐标系与第二从坐标系的夹角为δ2，其中，δ1＝αb1-αa1，δ2＝αb2-αa2，δ3＝αb3-αa3；第一从坐标系通过旋转角度δ1转换至主坐标系，第二从坐标系通过旋转角度δ3再平移转换至主坐标系；
如图4所示，在所述步骤2中，具体包括以下步骤：
步骤S21，对爆炸区域的光场信息利用每个光场相机内部的微透镜阵列进行处理，每一光场相机内的CCD成像探测面置于微透镜阵列的一倍焦距f处，设光线进入微透镜的位置坐标(ui,vi)，位置坐标矩阵为U＝[u1,u2,…un]T，V＝[v1,v2,…vn]T，其中，i＝1,2,···,n；经过图像处理技术，获取空间中每一破片在CCD成像探测面上的图像坐标位置(xi′,yi′)，其中，设定图像坐标矩阵为X′＝[x1′,x′2,…x′n]T，Y′＝[y1′,y′2,…y′n]T；当光场相机采集图像序列的某时刻出现多个破片目标重合或者遮挡的情况时，每个光场相机会采集对应时刻的一帧图像，根据三个光场相机的空间几何布置，在某一光场相机采集的该帧图像中重合遮挡的破片目标，在另一光场相机成像后会被分离；通过改变微透镜阵列参数ai，使得重合遮挡破片目标区域的成像景深变浅，在还原被遮挡破片时，能够将前方的离焦遮挡破片严重虚化，实现透视效果；
步骤S22，将微透镜阵列平面视为u-v平面，将探测器平面视为x-y平面，则重聚焦平面与微透镜平面的距离为aif，对遮挡、重合的破片目标分辨识别后，使得空间破片通过重新聚焦在CCD成像探测面上形成一个清晰的像，从而得到破片目标经过重聚焦后的图像坐标(xic,yic)，则Xc＝[x1c,x2c,…xnc]T，Yc＝[y1c,y2c,…ync]T，且另其中，i＝1,2,···,n，且满足如下关系：
如图5-6所示，在所述步骤3中，具体包括以下步骤：
步骤S31，根据步骤2中的数字重聚焦方法获取每个破片目标的图像坐标(xic,yic)，计算每个破片与各个光场相机光轴之间的夹角βmi，m＝1,2,3，i＝1,2,···,n则有：
如图7-9所示，在所述步骤4中，包括以下步骤：
步骤S41，当破片群4经过探测区域时，三个光场相机同步对破片群图像进行采集，依据采集的图像，调整微透镜阵列参数矩阵A进行分辨识别破片目标，解算图像坐标参数模型(Xc,Yc)；假设每台光场相机的焦距为fl，由步骤2可得到破片经过重聚焦后的图像坐标(xi′,yi′)，计算每个破片与各个相机光轴之间的夹角βmi，其中，i＝1,2,···,n，l＝1,2,3，m＝1,2,3；
步骤S42，分别对两两交汇的相机进行分析，根据几何关系计算第k幅光场相机采集的图像中第i个破片目标的空间位置(xki,yki,zki)，对应的空间位置矩阵(Xk,Yk,Zk)，其中，i＝1,2,…,n，k＝1,2,…，具体分为如下情况：
利用第一光场相机1与第三光场相机3视场交汇得到破片群的空间位置矩阵(X1,Y1,Z1)，具体如下：
X1＝d1·Q1 TW1E1-f1·R1 TM1
Y1＝d1·Q1 TW1E1
Z1＝d1·C1 TB1
其中，
Q1＝[sin(ε3-β31),sin(ε3-β32),…，sin(ε3-β3n)]T
W1＝[cos(ε1-β11),cos(ε1-β12),…，cos(ε1-β1n)]T
B1＝[cot(ε3+β31),cot(ε3+β32),…，cot(ε3+β3n)]T
R1＝[sin(β11),sin(β12),…，sin(β1n)]T
利用第一光场相机1与第二光场相机2视场交汇得到破片群的空间位置矩阵(X2,Y2,Z2)，具体如下：
X2＝(d2·Q2 TW2E2-f2·R1 TC2)·cosδ1-(d2·Q2 TE2C2·sinδ1)
Y2＝(d2·Q2 TW2E2-f2·R1 TC2)·sinδ1-(d2·Q2 TE2C2)·cosδ1
Z2＝d2·C1 TN1
其中，
N1＝[cot(ε2+β21),cot(ε2+β22),…，cot(ε2+β2n)]T
Q2＝[sin(ε2+β21),sin(ε2+β22),…，sin(ε2+β2n)]T
W2＝[cos(δ1i-ε1+β11),cos(δ1i-ε1+β12),…，cos(δ1i-ε1+β1n)]T
M2＝[sin(δ1-ε2+β11),sin(δ1-ε2+β12),…，sin(δ1-ε2+β1n)]T
利用第二光场相机2与第三光场相机3视场交汇得到破片群的空间位置矩阵(X3,Y3,Z3)，具体如下：
X3＝(d3·Q3 TW3E3-f3·R3 TZ3)·cosδ3-(d3·Q3 TZ3E3)·sinδ3+d2·cosδ1
Y3＝(d3·Q3 TW3E3-f3·R3 TZ3)·sinδ3-(d3·Q3 T·Z3E3)·cosδ3+d2·sinδ1
Z3＝d2C1 TN1
其中，
Q3＝[sin(δ3-ε3+β31),sin(δ3-ε3+β32),…，sin(δ3-ε3+β3n)]T
W3＝[cos(δ2-ε2-β21),cos(δ2-ε2-β22),…，cos(δ2-ε2-β2n)]T
R3＝[sin(β21),sin(β22),…，sin(β2n)]T
步骤S43，依据解算的空间位置矩阵，用相邻两幅图像待配准的破片目标群4位置坐标进行作差，根据同一破片目标的方位角与俯仰角在图像序列中不变的原理进行匹配，若坐标差成比例，则待配准目标为同一目标。
在所述步骤5中，包括以下步骤：
步骤S51，为了提高破片群4的测量精度，本实施例取三个坐标系的空间位置平均值，从而得到破片目标的空间位置(xi,yi,zi)，其中，破片空间位置坐标矩阵X＝[x1,x2,…xn]T，Y＝[y1,y2,…yn]T，Z＝[z1,z2,…zn]T，且满足如下关系：
显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。
Claims (6)
1.一种弹丸近炸破片群位置参数测量与计算方法，其具体包括以下步骤：
步骤1：根据破片群的位置，布置三个光场相机，根据三个光场相机布置参数，建立基于三个面阵CCD相机的近炸空间破片散布解算模型；
步骤2：利用数字重聚焦技术，对遮挡、重合的多个破片目标进行分辨识别，获取破片目标的图像坐标；
步骤3：结合破片目标经过数字重聚焦后的图像坐标，确定空间破片与各个光场相机的光轴之间的夹角；
步骤4：当破片目标经过探测区域时，利用重聚焦解算目标图像坐标参数模型；基于建立的近炸空间破片散布解算模型，依据双目视觉原理，对两两交汇的光场相机进行分析，根据空间几何关系计算出破片在主坐标系下与从坐标系下的空间位置坐标参数；依据三个坐标系下的空间位置参数，根据同一破片目标的方位角与俯仰角在图像序列中不变的原理进行破片群匹配；
步骤5：取三个坐标系下的空间位置坐标平均值，获取破片目标的空间位置。
2.根据权利要求1所述的弹丸近炸破片群位置参数测量与计算方法，其特征在于，在所述步骤1中，包括以下步骤：
步骤S11，采用三个光场相机，且三个光场相机都具有满足测试要求的高帧频，将其布置到三角形的三个顶点，且使三个光场相机的视场交汇于探测区域；
步骤S12，对三个光场相机中两两一组进行校平；
步骤S13，测量每一组光场相机校平后的角度值αa1、αa3、αb1、αb2、αa2、αb3；
步骤S14，测量三个光场相机光学系统的俯仰角数值ε1、ε2、ε3；
步骤S15，测量每两个光场相机的光学镜头之间的距离d1、d2、d3；依据建立的破片空间位置解算模型，在炮口处设置外部触发源，当触发源探测火光时，三个光场相机同步触发采集图像,；
步骤S16，建立主坐标系，具体地，以第一光场相机(1)和第三光场相机(3)的中心连线S1S3为X轴，S1为原点，构成主坐标系XOY；
步骤S17，建立两个从坐标系，其中，以第一光场相机(1)和第二光场相机(2)的中心连线S1S2为X'轴，S1为原点O'，构成第一从坐标系；以第二光场相机(2)和第三光场相机(3)的中心连线S2S3为X″轴，S2为原点O″，构成第二从坐标系；其中，主坐标系与第一从坐标系的夹角为δ1，主坐标系与第二从坐标系的夹角为δ3，第一从坐标系与第二从坐标系的夹角为δ2，其中，δ1＝αb1-αa1，δ2＝αb2-αa2，δ3＝αb3-αa3。
3.根据权利要求1所述的弹丸近炸破片群位置参数测量与计算方法，其特征在于，在所述步骤2中，具体包括以下步骤：
步骤S21，对爆炸区域的光场信息利用光场相机内部的微透镜阵列进行处理，每一光场相机内的CCD成像探测面置于微透镜阵列的一倍焦距f处，设光线进入微透镜的位置坐标(ui,vi)，位置坐标矩阵U＝[u1,u2,…un]T，V＝[v1,v2,…vn]T，其中，i＝1,2,…,n；经过图像处理技术，获取空间中每一破片在CCD成像探测面上的图像坐标位置(x′i,y′i)，设图像坐标矩阵X′＝[x′1,x′2,…x′n]T，Y′＝[y′1,y′2,…y′n]T；当光场相机采集图像中出现多个破片目标重合或者遮挡的情况时，通过改变微透镜阵列参数ai，使重合遮挡破片目标区域的成像景深变浅，在还原被遮挡破片时，能够将前方的离焦遮挡破片严重虚化，实现透视效果；
步骤S22，将微透镜阵列平面视为u-v平面，探测器平面视为x-y平面，则重聚焦平面距微透镜平面的距离为aif，对遮挡、重合的破片目标分辨识别后，使得空间破片通过重新聚焦在CCD成像探测面上形成一个清晰的像，从而得到破片目标经过重聚焦后的图像坐标(xic,yic)，则Xc＝[x1c,x2c,…xnc]T，Yc＝[y1c,y2c,…ync]T，且另其中，i＝1,2,…,n，E为单位矩阵，且满足如下关系：
4.根据权利要求1所述的弹丸近炸破片群位置参数测量与计算方法，其特征在于，在所述步骤3中，具体包括以下步骤：
步骤S31，根据步骤2中的数字重聚焦方法获取每个破片目标的图像坐标(xic,yic)，计算每个破片与各个光场相机光轴之间的夹角βmi，m＝1,2,3，i＝1,2,…,n则有：
5.根据权利要求1所述的弹丸近炸破片群位置参数测量与计算方法，其特征在于，在所述步骤4中，包括以下步骤：
步骤S41，当破片群经过探测区域时，三个光场相机同步对破片群图像进行采集，依据采集的图像，调整微透镜阵列参数矩阵A进行分辨识别破片目标，解算图像坐标参数模型(Xc,Yc)；假设每台光场相机的焦距为fl，根据破片经过重聚焦后的图像坐标(x′i,y′i)，计算每个破片与各个光场相机光轴之间的夹角βmi，其中，i＝1,2,…,n，l＝1,2,3，m＝1,2,3。
步骤S42，分别对两两交汇的光场相机进行分析，根据几何关系计算第k幅光场相机采集的图像中第i个破片目标的空间位置(xki,yki,zki)，对应的空间位置矩阵(Xk,Yk,Zk)，其中，i＝1,2,…,n，k＝1,2,…，具体分为如下情况：
利用第一光场相机(1)与第三光场相机(3)视场交汇得到破片群的空间位置矩阵(X1,Y1,Z1)，具体如下：
X1＝d1·Q1 TW1E1-f1·R1 TM1
Y1＝d1·Q1 TW1E1
Z1＝d1·C1 TB1
其中，
Q1＝[sin(ε3-β31),sin(ε3-β32),…，sin(ε3-β3n)]T
W1＝[cos(ε1-β11),cos(ε1-β12),…，cos(ε1-β1n)]T
B1＝[cot(ε3+β31),cot(ε3+β32),…，cot(ε3+β3n)]T
R1＝[sin(β11),sin(β12),…，sin(β1n)]T
利用第一光场相机(1)与第二光场相机(2)视场交汇得到破片群的空间位置矩阵(X2,Y2,Z2)，具体如下：
X2＝(d2·Q2 TW2E2-f2·R1 TC2)·cosδ1-(d2·Q2 TE2C2·sinδ1)
Y2＝(d2·Q2 TW2E2-f2·R1 TC2)·sinδ1-(d2·Q2 TE2C2)·cosδ1
Z2＝d2·C1 TN1
其中，
N1＝[cot(ε2+β21),cot(ε2+β22),…，cot(ε2+β2n)]T
Q2＝[sin(ε2+β21),sin(ε2+β22),…，sin(ε2+β2n)]T
W2＝[cos(δ1i-ε1+β11),cos(δ1i-ε1+β12),…，cos(δ1i-ε1+β1n)]T
M2＝[sin(δ1-ε2+β11),sin(δ1-ε2+β12),…，sin(δ1-ε2+β1n)]T
利用第二光场相机(2)与第三光场相机(3)视场交汇得到破片群的空间位置矩阵(X3,Y3,Z3)，具体如下：
X3＝(d3·Q3 TW3E3-f3·R3 TZ3)·cosδ3-(d3·Q3 TZ3E3)·sinδ3+d2·cosδ1
Y3＝(d3·Q3 TW3E3-f3·R3 TZ3)·sinδ3-(d3·Q3 T·Z3E3)·cosδ3+d2·sinδ1
Z3＝d2C1 TN1
其中，
Q3＝[sin(δ3-ε3+β31),sin(δ3-ε3+β32),…，sin(δ3-ε3+β3n)]T
W3＝[cos(δ2-ε2-β21),cos(δ2-ε2-β22),…，cos(δ2-ε2-β2n)]T
R3＝[sin(β21),sin(β22),…，sin(β2n)]T
步骤S43，依据解算的空间位置矩阵，用相邻两幅图像待配准的破片目标群位置坐标进行作差，根据同一破片目标的方位角与俯仰角在图像序列中不变的原理进行匹配，若坐标差成比例，则待配准目标为同一目标。
6.根据权利要求1所述的弹丸近炸破片群位置参数测量与计算方法，其特征在于，在所述步骤5中，包括以下步骤：
步骤S51，取三个坐标系的空间位置平均值，从而得到破片目标的空间位置(xi,yi,zi)，其中，破片空间位置坐标矩阵X＝[x1,x2,…xn]T，Y＝[y1,y2,…yn]T，Z＝[z1,z2,…zn]T，且满足如下关系：
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810438708.3A CN108896017B (zh) | 2018-05-09 | 2018-05-09 | 一种弹丸近炸破片群位置参数测量与计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810438708.3A CN108896017B (zh) | 2018-05-09 | 2018-05-09 | 一种弹丸近炸破片群位置参数测量与计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108896017A true CN108896017A (zh) | 2018-11-27 |
| CN108896017B CN108896017B (zh) | 2022-04-15 |
Family
ID=64342780
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810438708.3A Active CN108896017B (zh) | 2018-05-09 | 2018-05-09 | 一种弹丸近炸破片群位置参数测量与计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108896017B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109829945A (zh) * | 2018-11-28 | 2019-05-31 | 西安工业大学 | 一种近炸破片分布场弹目交汇的目标毁伤评估方法 |
| CN112435252A (zh) * | 2020-12-04 | 2021-03-02 | 西安工业大学 | 一种战斗部破片穿孔和凹坑检测方法 |
| CN113008089A (zh) * | 2021-02-18 | 2021-06-22 | 西安工业大学 | 一种双光场相机破片参数测试装置与测试方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2216996A (en) * | 1988-03-14 | 1989-10-18 | Sector Limited | System for detecting position of impact of a projectile. |
| CN102175149A (zh) * | 2011-01-25 | 2011-09-07 | 西安工业大学 | 一种飞行弹丸空间炸点三维坐标光电测量装置及测量方法 |
| CN102818561A (zh) * | 2012-07-09 | 2012-12-12 | 中国人民解放军国防科学技术大学 | 基于数字狭缝摄像技术的靶场弹丸运动参数测量方法 |
| EP2887642A2 (en) * | 2013-12-23 | 2015-06-24 | Nokia Corporation | Method, apparatus and computer program product for image refocusing for light-field images |
-
2018
- 2018-05-09 CN CN201810438708.3A patent/CN108896017B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2216996A (en) * | 1988-03-14 | 1989-10-18 | Sector Limited | System for detecting position of impact of a projectile. |
| CN102175149A (zh) * | 2011-01-25 | 2011-09-07 | 西安工业大学 | 一种飞行弹丸空间炸点三维坐标光电测量装置及测量方法 |
| CN102818561A (zh) * | 2012-07-09 | 2012-12-12 | 中国人民解放军国防科学技术大学 | 基于数字狭缝摄像技术的靶场弹丸运动参数测量方法 |
| EP2887642A2 (en) * | 2013-12-23 | 2015-06-24 | Nokia Corporation | Method, apparatus and computer program product for image refocusing for light-field images |
Non-Patent Citations (3)
| Title |
|---|
| 李翰山等: "基于摄像法测量弹丸的空间炸点位置", 《光学精密工程》 * |
| 李良威等: "高速摄影法测量弹丸的破片速度", 《高速摄影与光子学》 * |
| 桑晓月等: "弹丸近炸破片密度分布与运动规律研究", 《机械与电子》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109829945A (zh) * | 2018-11-28 | 2019-05-31 | 西安工业大学 | 一种近炸破片分布场弹目交汇的目标毁伤评估方法 |
| CN109829945B (zh) * | 2018-11-28 | 2022-11-18 | 西安工业大学 | 一种近炸破片分布场弹目交汇的目标毁伤评估方法 |
| CN112435252A (zh) * | 2020-12-04 | 2021-03-02 | 西安工业大学 | 一种战斗部破片穿孔和凹坑检测方法 |
| CN112435252B (zh) * | 2020-12-04 | 2023-05-09 | 西安工业大学 | 一种战斗部破片穿孔和凹坑检测方法 |
| CN113008089A (zh) * | 2021-02-18 | 2021-06-22 | 西安工业大学 | 一种双光场相机破片参数测试装置与测试方法 |
| CN113008089B (zh) * | 2021-02-18 | 2023-05-16 | 西安工业大学 | 一种双光场相机破片参数测试装置与测试方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108896017B (zh) | 2022-04-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107144240A (zh) | 一种检测玻璃面板表面缺陷的系统及方法 | |
| CN108896017A (zh) | 一种弹丸近炸破片群位置参数测量与计算方法 | |
| CN108426585B (zh) | 一种光场相机的几何标定方法 | |
| CN108986070B (zh) | 一种基于高速视频测量的岩石裂缝扩展实验监测方法 | |
| CN106553086B (zh) | 快速高精度的复杂曲面制孔点法矢量测量方法 | |
| CN109658457A (zh) | 一种激光与相机任意相对位姿关系的标定方法 | |
| CN103604371A (zh) | 移动终端及其物体测量方法 | |
| CN109341668B (zh) | 基于折射投影模型和光束追踪法的多相机测量方法 | |
| CN105866129A (zh) | 一种基于数字投影的产品表面质量在线检测方法 | |
| CN106052586A (zh) | 基于机器视觉的石材大板表面轮廓尺寸获取系统及方法 | |
| CN106875436B (zh) | 一种基于特征点密度由聚焦堆栈估计深度的方法和装置 | |
| CN107121061B (zh) | 一种基于椭圆锥约束的制孔点法矢量求解方法 | |
| CN103782232A (zh) | 投影仪及其控制方法 | |
| CN109089025A (zh) | 一种基于光场成像技术的影像仪数字聚焦方法 | |
| CN109377551A (zh) | 一种三维人脸重建方法、装置及其存储介质 | |
| CN105957096A (zh) | 一种用于三维数字图像相关的相机外参标定方法 | |
| CN107024829B (zh) | 多光谱相机像面装调方法 | |
| CN110487183A (zh) | 一种多目标光纤位置精确检测系统及使用方法 | |
| CN107167088B (zh) | 玻璃变形量的测量方法及装置 | |
| CN110146180A (zh) | 一种基于焦面哈特曼波前传感器的大视场图像清晰化装置及方法 | |
| CN103676487B (zh) | 一种工件高度测量装置及其校正方法 | |
| CN109900719A (zh) | 一种叶片表面刀纹的视觉检测方法 | |
| CN108872246A (zh) | 板面材料全视面缺陷检测系统 | |
| CN107084671A (zh) | 一种基于三线结构光的凹球直径测量系统及测量方法 | |
| CN108050955B (zh) | 基于结构光投影与数字图像相关的高温空气扰动滤除方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| EE01 | Entry into force of recordation of patent licensing contract | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20181127 Assignee: Shaanxi Yunrui Chuangzhi Technology Co.,Ltd. Assignor: XI'AN TECHNOLOGICAL University Contract record no.: X2023980044921 Denomination of invention: A Method for Measuring and Calculating the Position Parameters of the Near Burst Fragment Group of Projectiles Granted publication date: 20220415 License type: Common License Record date: 20231030 |