CN114329854B - 一种基于多目标约束的二维空间视觉传感器布局优化方法 - Google Patents

一种基于多目标约束的二维空间视觉传感器布局优化方法 Download PDF

Info

Publication number
CN114329854B
CN114329854B CN202011228559.1A CN202011228559A CN114329854B CN 114329854 B CN114329854 B CN 114329854B CN 202011228559 A CN202011228559 A CN 202011228559A CN 114329854 B CN114329854 B CN 114329854B
Authority
CN
China
Prior art keywords
model
coverage
point
space
sensor
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.)
Active
Application number
CN202011228559.1A
Other languages
English (en)
Other versions
CN114329854A (zh
Inventor
全权
邓恒
吴楠
蔡开元
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihang University
Original Assignee
Beihang University
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 Beihang University filed Critical Beihang University
Priority to CN202011228559.1A priority Critical patent/CN114329854B/zh
Publication of CN114329854A publication Critical patent/CN114329854A/zh
Application granted granted Critical
Publication of CN114329854B publication Critical patent/CN114329854B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/10Internal combustion engine [ICE] based vehicles
    • Y02T10/40Engine management systems

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Testing Or Calibration Of Command Recording Devices (AREA)

Abstract

本发明公开一种基于多目标约束的二维空间视觉传感器布局优化方法,具体如下:步骤一、建立视觉传感网络模型,所述的模型包括感知模型、空间模型、覆盖模型和重建误差模型;步骤二、定义覆盖率和重建品质;步骤三、总体设计多目标优化。本发明提出的一种基于多目标约束的二维空间视觉传感器布局优化方法,该方法能够在给定的二维空间中,根据覆盖范围大和重建误差小的目标要求,设计多目标优化算法,提供一种合理的摄像机布局方案,指导摄像机实际安装调试,能够大大节省人力物力、加快网络部署,降维搜索的方法可以解决任意形状空间的传感器布局问题,也能减少优化变量个数,提高优化速度。

Description

一种基于多目标约束的二维空间视觉传感器布局优化方法
技术领域
本发明涉及一种基于多目标约束的二维空间视觉传感器布局优化方法,属于视觉传感器布局优化领域。
背景技术
基于视觉传感器的三维视觉跟踪技术在很多领域有着广泛应用,如自主导航、人机交互、空中加油、增强/虚拟现实等。以视觉导航为例:与全球定位系统GPS、惯性导航系统等传统方法相比,基于视觉的定位方式精度高、信息丰富、抗干扰能力强、适用于室内环境,近年来受到国内外广泛关注。如何部署视觉传感器是视觉传感器网络应用中的重要问题之一,它直接影响着网络的覆盖范围、连通性和能源利用效率。目前的实践中,大多数摄像机布局问题都是通过人工反复试验得到解决,即不断地调整布局然后试验判断结果。相比之下,摄像机布局的自动决策问题会更有意义,一方面可以节约摄像机反复布局安装的时间和人力物力成本,一方面通过数值计算可以获得更加鲁棒的解决方法。同时,针对视觉定位任务,摄像机布局不仅需要满足覆盖范围要求,而且需要满足目标至少被两个摄像机同时覆盖的约束条件,从而实现三维重建。因此,本发明提出了一种基于多目标约束的二维空间视觉传感器布局优化方法,该方法能够在给定的二维空间中,根据覆盖范围大和重建误差小的目标要求,设计多目标优化算法,提供一种合理的摄像机布局方案(包括摄像机位置和方位角)。该方法直接在视觉传感网络布置之前给出摄像机部署方案,指导摄像机实际安装调试,能够大大节省人力物力、加快网络部署。
发明内容
本发明提供了一种基于多目标约束的二维空间视觉传感器布局优化方法,它同时考虑了视觉传感网络覆盖范围大和定位误差小的要求,解决了广域空间中视觉传感器部署困难的问题。该方法以简单的二维空间出发,适用于大多数方向传感器的布局优化问题,它能够解决大范围的布局问题,不需要额外增加复杂的测量工具,方便快捷,对实际工程部署有很强的指导意义。
本发明中的问题可以简单描述为:给定规则的二维矩形空间、摄像机模型和摄像机总数,优化摄像机布局方案(包括摄像机位置和方位角)使得目标空间覆盖率大的同时满足重建误差小。将上述问题抽象为一个优化问题。令
Figure BDA0002764406840000021
为任意连通的局部二维区域。
Figure BDA0002764406840000022
为给定任务,
Figure BDA0002764406840000023
为包含
Figure BDA0002764406840000024
所需的所有约束的集合。这些约束可能包括二维空间(即覆盖范围)
Figure BDA0002764406840000025
的位置约束,传感器的约束(传感器类型数目等)。传感器布局优化问题可以定义:给定区域
Figure BDA0002764406840000026
找到一组满足约束条件
Figure BDA0002764406840000027
的传感器位置,并最大化给定的目标函数G(或者最小化损失函数)。用数学的形式表示为:
Figure BDA0002764406840000028
其中,Π={π12,…,πN}为所有摄像机的布局方式,N为优化后的总摄像机数目。
注意到上述定义(1)只是一个抽象的表达方式,在具体问题面前,需要具体化每一部分,即具体确定优化的空间区域
Figure BDA0002764406840000029
不同的任务
Figure BDA00027644068400000210
不同的约束条件
Figure BDA00027644068400000211
以及最终的目标函数G,在本方法中,给出了一种具体的模型建立以及多目标优化问题解决方法,将最大化覆盖范围和最小化重建误差作为目标函数,提出一种加权的目标函数综合考虑它们。
一种基于多目标约束的二维空间视觉传感器布局优化方法,为离线方法,在视觉传感器实际部署之前,对整个视觉网络的传感器进行布局设计,具体如下:
步骤一、建立视觉传感网络模型,所述的模型包括感知模型、空间模型、覆盖模型和重建误差模型;
其中,感知模型采用等腰三角形模型。该等腰三角形的参数由摄像机的参数决定:其顶点为摄像机的二维坐标(cx,cy),底边上的高为摄像机的有效感知距离dof,顶角为摄像机的视场角fov,底边上的高与x轴的夹角为摄像机的方位角
Figure BDA00027644068400000212
则底边a=2doftan(fov/2)。这样,摄像机的位姿参数记为坐标Cp=(cx,cy)和方位角
Figure BDA00027644068400000213
则一个已知类型的摄像机的状态可用参数
Figure BDA00027644068400000214
表示。将待测目标点的位置坐标设为(x,y),经过数学推导,可将该特定摄像机的感知区域用以下三个不等式表示:
Figure BDA0002764406840000031
空间模型:将目标空间
Figure BDA0002764406840000032
分为两部分:边界
Figure BDA0002764406840000033
和测量区域
Figure BDA0002764406840000034
分别对应房间的墙壁和内部空间。将目标空间内的一个点称为工作点,用Wp表示。由于只关注测量区域内的覆盖和重建问题,有
Figure BDA0002764406840000035
同理,由于传感器只能放置在边界上,有
Figure BDA0002764406840000036
将目标空间抽象为无孔洞的矩形,则边界
Figure BDA0002764406840000037
为矩形的边,测量区域
Figure BDA0002764406840000038
为矩形的内部区域。采用均匀网格采样模型,边界
Figure BDA0002764406840000039
和测量区域
Figure BDA00027644068400000310
使用不同采样间隔,分别记为Δb和Δm,同时有Δb≤Δm。同样采用均匀采样的方法将圆周分为
Figure BDA00027644068400000311
份,即采样间隔为
Figure BDA00027644068400000312
确定采样间隔后,计算目标空间边界区域在x轴和y轴上采样点个数分别为sxb,syb,测量区域在x轴和y轴上采样点个数分别为sxm,sym,系统中方位角只能取
Figure BDA00027644068400000313
个值。
覆盖模型采用普遍使用的0/1覆盖模型,处在传感器感知范围内的点状态为被覆盖,反之则为不被覆盖。另外,考虑到图像重建的理论要求,本方法考虑k-覆盖问题,即空间中任一点至少被k个传感器节点覆盖。
重建误差模型使用测量方程的雅可比行列式定义的几何精度因子GDOP(Geometric Dilution of Precision)。GDOP针对工作空间中的每个点,定义观察它的两个传感器导致的不确定性的度量。
步骤二、定义覆盖率和重建品质
1)定义覆盖率
定义以下二进制变量:
Figure BDA00027644068400000314
其中,i1∈[0,sxb-1],j1∈[0,syb-1],i2∈[0,sxm-1],j2∈[0,sym-1],
Figure BDA00027644068400000315
将以上离散值转换为连续目标空间中的坐标和角度后,
Figure BDA00027644068400000316
的值根据公式(2)计算得出,并将所有的值提前储存在矩阵中,以便随时调用。
确定参数组合
Figure BDA0002764406840000041
和一组传感器的布局方案后,可以定义二进制可见性矩阵V=(vij),其中元素为:
Figure BDA0002764406840000042
其中,任意vij的值可以通过查询式(3)得出。那么,用二进制变量bj表示点j是否满足k-覆盖:
Figure BDA0002764406840000043
则目标空间的覆盖率定义为:满足k-覆盖的采样点个数与总采样点个数之比:
Figure BDA0002764406840000044
2)定义重建品质
在二维空间中视觉传感网络中,两个视觉传感器坐标为(xi,yi),i=1,2,分别测量目标点相对于自身的方位角θi,i=1,2,而目标点位置必须是同时满足两个角度约束的点,即两条射线的交点。这里将θi定义为从x轴正半轴开始逆时针旋转过的角度,取值范围为[0,2π]。ri,i=1,2为两个视觉传感器相对空间点的距离向量,且其模长为|ri|=ri,i=1,2。经过数学推导,可以得到GDOP等于目标与两个传感器之间的距离和两个传感器之间的相对方位角θ的余割的乘积,具体表达式为:
Figure BDA0002764406840000045
几何精度因子(7)表示重建过程将传感器测量误差放大为重建误差的程度,则其倒数可以表示重建结果的品质,再选择合适的比例系数dscale,将品质函数限制在[0,1],可推导出目标点j的重建品质函数:
Figure BDA0002764406840000046
上式(8)给出传感器i和传感器i′对于目标点j的重建品质。其中,vij,vi′j的表达式见式(4),乘数vijvi′j保证了当采样点
Figure BDA0002764406840000051
没有被2-覆盖时,重建品质qii′j=0。
遍历目标空间中的每个采样点,求取它们的重建品质和并求平均,即可得到整个空间的重建品质:
Figure BDA0002764406840000052
步骤三、总体设计多目标优化
本方法用空间覆盖率和空间重建品质的加权和作为综合评价传感器系统性能的量化标准,即:
e=σr+(1-σ)q         (10)
其中,r,q的表达式见公式(6)和公式(9),加权系数σ∈[0,1]表示评价标准倾向于覆盖率的程度。当σ=1时,评价标准为单一的最大化覆盖率;σ=0时,则为最大化重建品质,即最小化重建误差。
综合上述分析,二维优化问题可整理为以下模型:
给定目标区域长和宽,采样参数组合
Figure BDA0002764406840000053
传感器数量N,优化以下问题:
Figure BDA0002764406840000054
其中,
Figure BDA0002764406840000055
G=r,q,or e。
针对离散优化问题(11),本发明采用遗传算法(GA,Genetic Algorithm)和粒子群算法(PSO,Particle Swarm Optimization)的组合算法来求解,在提高运行速度的同时,尽可能避免算法陷入局部最优解。组合算法流程如下:
a)求取初始种群:将最大迭代次数设为50,运行粒子群算法,得到初步优化后的种群P;
b)筛选初始种群:保留种群P中最佳的9个个体,以及整个优化过程中的最佳个体,共10个个体作为部分初始种群,不足的部分则随机生成,两部分共同构成初始种群P(0);
c)迭代:以P(0)为初始种群,运行遗传算法;
d)判断终止条件:达到终止条件后,整个运算过程中适应度最大的个体为最优解,终止算法。
为简化优化参数和约束条件,本方法提出了一种用一维参数表示传感器位置信息的降维搜索方法。若将直角坐标系原点设为起点,逆时针沿目标空间的边界运动,则边界上的每个离散点唯一对应一个运动的前向距离d。以此前向距离代表传感器的位置信息,即可将二维参数(cx,cy)降为一维参数d。这样,针对有N个传感器的二维网络,其中第i个传感器的位姿信息由参数组合
Figure BDA0002764406840000061
简化为
Figure BDA0002764406840000062
整个网络的决策变量即为
Figure BDA0002764406840000063
本发明提出的一种基于多目标约束的二维空间视觉传感器布局优化方法,该方法能够在给定的二维空间中,根据覆盖范围大和重建误差小的目标要求,设计多目标优化算法,提供一种合理的摄像机布局方案(包括摄像机位置和方位角),指导摄像机实际安装调试,能够大大节省人力物力、加快网络部署,降维搜索的方法可以解决任意形状空间的传感器布局问题,也能减少优化变量个数,提高优化速度。
附图说明
图1是等腰三角形感知模型示意图;
图2是二维空间均匀采样示意图;
图3是二维空间重建原理示意图;
图4是本方法的流程图;
图5是6台相机优化布局实验结果示意图;
图中符号说明如下:
图1中的符号说明:oxy表示二维空间坐标轴,顶点表示摄像机的坐标(cx,cy),底边上的高表示摄像机的感知距离dof,顶角为摄像机的视场角fov,底边上的高与x轴的夹角为摄像机的方位角
Figure BDA0002764406840000064
图3中的符号说明:oxy表示二维空间坐标轴,两个视觉传感器坐标为(xi,yi),分别测量目标点r=(x,y)相对于自身的方位角θi,θ为两个视觉传感器之间的相对方位角,ri为两个视觉传感器相对空间点的距离。
具体实施方式
以下结合附图和实施例,对本发明的技术方案做进一步的说明。
本发明提出了一种基于多目标约束的二维空间视觉传感器布局优化方法,是一种离线方法,在视觉传感器实际部署之前,就需要对整个视觉网络的传感器进行合理的布局设计。整个离线优化过程的实现步骤具体如下:
(1)视觉传感网络模型建立
本部分的主要目的是需要将真实的摄像机和物理空间抽象为数学模型,用数学语言描述了传感器布局优化问题,为建立优化目标函数打下基础,以便量化地分析与解决问题。所涉及到的模型包括感知模型、空间模型、覆盖模型和重建误差模型。
1)感知模型建立
在二维空间问题中,摄像机感知模型采用等腰三角形模型(如图1)。该等腰三角形的参数由摄像机的参数决定:其顶点为摄像机的二维坐标(cx,cy),底边上的高为摄像机的有效感知距离dof,顶角为摄像机的视场角fov,底边上的高与x轴的夹角为摄像机的方位角
Figure BDA0002764406840000071
则底边a=2doftan(fov/2)。这样,摄像机的位姿参数记为坐标Cp=(cx,cy)和方位角
Figure BDA0002764406840000072
则一个已知类型的摄像机的状态可用参数
Figure BDA0002764406840000073
表示。将待测目标点的位置坐标设为(x,y),经过数学推导,可将该特定摄像机的感知区域用以下三个不等式表示:
Figure BDA0002764406840000074
2)空间模型建立
空间模型描述了传感器网络需要覆盖的区域,或称之为目标空间。在真实世界里,一个房间的四周是墙壁,摄像机可以放置在墙壁上,用于测量处于房间内部的目标物体的位姿信息,以达到视觉跟踪的目的。故将目标空间
Figure BDA0002764406840000075
分为两部分:边界
Figure BDA0002764406840000076
和测量区域
Figure BDA0002764406840000077
分别对应房间的墙壁和内部空间。将目标空间内的一个点称为工作点,用Wp表示。由于只关注测量区域内的覆盖和重建问题,有
Figure BDA0002764406840000078
同理,由于传感器只能放置在边界上,有
Figure BDA0002764406840000079
将目标空间抽象为无孔洞的矩形,则边界
Figure BDA0002764406840000081
为矩形的边,测量区域
Figure BDA0002764406840000082
为矩形的内部区域。如图2所示,采用均匀网格采样模型,边界
Figure BDA0002764406840000083
和测量区域
Figure BDA0002764406840000084
使用不同采样间隔,分别记为Δb和Δm。为保证最终布局方案的精度,一般情况下有:Δb≤Δm。同样采用均匀采样的方法,将圆周分为
Figure BDA0002764406840000085
份,即采样间隔为
Figure BDA0002764406840000086
确定采样间隔后,可计算出目标空间边界区域在x轴和y轴上采样点个数分别为sxb,syb,测量区域在x轴和y轴上采样点个数分别为sxm,sym,系统中方位角只能取
Figure BDA0002764406840000087
个值。
这样,C和Wp的取值范围均为有限离散点,连续问题就转化为离散问题。当采样参数组合
Figure BDA0002764406840000088
选择合理时,可认为离散问题的解近似于连续问题的解。
3)覆盖模型建立
针对覆盖问题,从传感器特性的角度,需要描述一个特定的传感器可以覆盖的范围。本方法采用普遍使用的0/1覆盖模型,即一个空间点只可能有两种状态:被覆盖和不被覆盖,处在传感器感知范围内的点状态为被覆盖,反之则为不被覆盖。另外,从目标空间的角度,还可以给出覆盖度的概念:若一个区域中任一点被k个传感器节点覆盖,则称这片区域的覆盖度为k。如何保证一个区域的覆盖度等于k,是提高传感器网络覆盖质量的关键,这个问题也被称为“k-覆盖问题”。考虑到图像重建的理论要求,本方法考虑k-覆盖问题。
4)重建误差模型
为了计算仅由传感器放置的几何特性引起的不确定性,本方法使用测量方程的雅可比行列式定义一种宽松的几何精度因子(GDOP,Geometric Dilution of Precision)。GDOP针对工作空间中的每个点,定义观察它的两个传感器导致的不确定性的度量。
(2)覆盖率和重建品质定义
1)覆盖率定义
根据连续的覆盖模型,位于(cx,cy),方位角为
Figure BDA0002764406840000089
的传感器当且仅当公式(2)成立时能覆盖点(x,y)。分析离散的目标空间,使用离散化后的传感器位姿参数和测量点坐标,我们可以定义以下二进制变量:
Figure BDA0002764406840000091
任取i1∈[0,sxb-1],j1∈[0,syb-1],i2∈[0,sxm-1],j2∈[0,sym-1],
Figure BDA0002764406840000092
将以上离散值转换为连续目标空间中的坐标和角度后,c的值仍可以通过式公式(2)计算得出。在优化过程中,将遍历计算所有c的值并提前储存在矩阵中,以便随时调用。
确定参数组合
Figure BDA0002764406840000093
和一组传感器的布局方案后,可以定义二进制可见性矩阵V=(vij),其中元素为:
Figure BDA0002764406840000094
其中,任意vij的值可以通过查询式(3)得出。那么,能够覆盖空间点j的传感器数量为:
nj=∑ivij                                 (5)
用二进制变量表示点j是否满足k-覆盖:
Figure BDA0002764406840000095
测量区域内被k个以上的传感器覆盖的所有采样点数量为:
Nc=∑jbj                                 (7)
则目标空间的覆盖率定义为:满足k-覆盖的采样点个数与总采样点个数之比:
Figure BDA0002764406840000096
这样覆盖率r便能够衡量目标空间被覆盖情况的量化标准。为保证目标空间被尽可能覆盖,需要最大化该函数。
2)重建品质函数定义
图3解释了二维空间中视觉传感网络的重建原理,两个视觉传感器坐标为(xi,yi),i=1,2,分别测量目标点相对于自身的方位角θi,i=1,2,而目标点位置必须是同时满足两个角度约束的点,即两条射线的交点。这里将θi定义为从x轴正半轴开始逆时针旋转过的角度,取值范围为[0,2π]。ri,i=1,2为两个视觉传感器相对空间点的距离向量,且其模长为ri=ri,i=1,2。经过数学推导,可以得到GDOP等于目标与两个传感器之间的距离和两个传感器之间的相对方位角θ的余割的乘积,具体表达式为:
Figure BDA0002764406840000101
对于这种情况,当目标位于两个传感器的连线上时(θ=0),GDOP趋于无穷大;当两个传感器之间的角度变大时,GDOP降低。当传感器之间的角度为直角且目标与传感器接近时,GDOP最佳。
几何精度因子(9)表示重建过程将传感器测量误差放大为重建误差的程度,则其倒数可以表示重建结果的品质,再选择合适的比例系数dscale,将品质函数限制在[0,1],可推导出品质函数:
Figure BDA0002764406840000102
上式(10)给出传感器i和传感器i′对于目标点j的重建品质。其中,vij,vi′j的表达式见式(4),乘数vijvi′j保证了当采样点
Figure BDA0002764406840000103
没有被2-覆盖时,重建品质qii′j=0。
遍历目标空间中的每个采样点,求取它们的重建品质和并求平均,即可得到整个空间的平均重建品质:
Figure BDA0002764406840000104
平均重建品质与平均重建误差成反比例关系,平均重建品质越高,则平均重建误差越小,两者均可以衡量系统的测量精度。为与覆盖范围函数统一优化方法,本方法中使用最大化平均重建品质作为优化目标。
(3)多目标优化总体设计
本部分的主要目的是利用上部分建立的数学模型,以及二维空间覆盖率r和重建品质q,再根据这两个目标函数的特点,提出一种多目标优化函数的建立方法。同时设计合适的多目标优化问题解决方法,以解决传感器网络布局优化问题。可以看出,上述两个评价指标值均在[0,1]之间,数量级相同。这样,本方法用它们的加权和作为综合评价传感器系统性能的量化标准,即:
e=σr+(1-σ)q         (12)
其中,r,q的表达式见公式(8)和公式(11),加权系数σ∈[0,1]表示评价标准倾向于覆盖率的程度。当σ=1时,评价标准为单一的最大化覆盖率;σ=0时,则为最大化重建品质,即最小化重建误差。
综合上述分析,二维优化问题可整理为以下模型:
给定目标区域长和宽,采样参数组合
Figure BDA0002764406840000111
传感器数量N,优化以下问题:
Figure BDA0002764406840000112
其中,
Figure BDA0002764406840000113
G=r或q或e。
针对离散优化问题(13),本方法采用遗传算法(GA,Genetic Algorithm)和粒子群算法(PSO,Particle Swarm Optimization)的组合算法来求解。遗传算法较为复杂,运行速度慢,但不容易陷入局部最优解;而粒子群算法较为简单,运行速度快,但容易出现过早收敛到局部最优解的情况。故本方法中组合算法的想法是:用粒子群算法的结果作为遗传算法的初值,以结合两者优点,在提高运行速度的同时,尽可能避免算法陷入局部最优解。组合算法流程如下:
a)求取初始种群:将最大迭代次数设为50,运行粒子群算法,得到初步优化后的种群P;
b)筛选初始种群:保留种群P中最佳的9个个体,以及整个优化过程中的最佳个体,共10个个体作为部分初始种群,不足的部分则随机生成,两部分共同构成初始种群P(0);
c)迭代:以P(0)为初始种群,运行遗传算法;
d)判断终止条件:达到终止条件后,整个运算过程中适应度最大的个体为最优解,终止算法。
为简化优化参数和约束条件,本方法提出了一种用一维参数表示传感器位置信息的降维搜索方法。若将直角坐标系原点设为起点,逆时针沿目标空间的边界运动,则边界上的每个离散点唯一对应一个运动的前向距离d。以此前向距离代表传感器的位置信息,即可将二维参数(cx,cy)降为一维参数d。这样,针对有N个传感器的二维网络,其中第i个传感器的位姿信息由参数组合
Figure BDA0002764406840000121
简化为
Figure BDA0002764406840000122
整个网络的决策变量即为
Figure BDA0002764406840000123
降维搜索方法不仅减少了优化参数,还直接将搜索范围限定在矩形的边界,省去了额外定义的约束条件,可避免搜索无用解空间、简化搜索过程、提高运行速度。同时,该方法理论上还可以解决任意形状的空间内的传感器布局问题。由于只通过前向距离d一个参数定义传感器位置,针对不规则形状的空间,只要已知空间边界的数学表达式,即可使用该降维搜索方法。
至此,一种基于多目标约束的二维空间视觉传感器布局优化方法已经完成,简要流程图见图4。
本发明提出了基于多目标优化的二维空间视觉传感器布局优化方法,具体仿真实施例如下:
仿真实验所用的电脑配置如下:Intel(R)Core(TM)i7-6700HQ CPU,主频2.6GHz,内存8GB;代码在MATLAB R2017b上编写并运行。根据实际需求,选择合理的环境参数、传感器参数和优化目标函数,运行仿真程序并得到布局方案。通过经验判断和理论计算,分析优化后的布局方案是否满足需求。
首先,考虑需要测量空间中的目标物体位姿参数的情况,即目标区域内采样点需要尽可能被两个以上传感器覆盖,且同时优化网络覆盖率和目标物体重建精度的情况。假设利用6台相机覆盖空间(60,60)m,相机视场角37°,视距60m,Δb=5,Δm=5,
Figure BDA0002764406840000125
同时最大化覆盖率和重建品质的优化结果如图5所示,可以看出整个目标区域覆盖情况良好(网络覆盖率94%),且传感器分布均匀合理,结果符合预期。
接着,讨论不同优化算法的对比试验,包括粒子群算法、遗传算法和本发明所述的组合算法。设计3组不同参数(见表1),分别使用3种算法优化,记录实验结果,以对比算法的性能。实验结果如表2所示,可以看出,在同一组参数设置下,粒子群算法的运行时间明显比遗传算法短,而优化结果略差于遗传算法,组合优化算法运行时间居于遗传算法和粒子群算法中间,且优化结果最好。
Figure BDA0002764406840000124
Figure BDA0002764406840000131
表1不同优化算法对比实验参数
Figure BDA0002764406840000132
表2不同优化算法对比实验结果
最后,讨论不同搜索方法的对比实验,包括常规搜索方法和本发明所述的降维搜索方法。选取遗传算法为最终的优化算法,实验选取3组不同的环境参数和传感器参数(见表3),分别用常规方法和降维方法开展实验,结果如表4。可以看到,降维搜索方法在显著降低运行时间的同时,还提高了算法的优化结果,说明其有同时提升算法速度和准确性的优点。
Figure BDA0002764406840000133
表3不同搜索方法的对比实验参数
Figure BDA0002764406840000134
Figure BDA0002764406840000141
表4不同搜索方法的对比实验结果。

Claims (3)

1.一种基于多目标约束的二维空间视觉传感器布局优化方法,其特征在于:为离线方法,在视觉传感器实际部署之前,对整个视觉网络的传感器进行布局设计,具体如下:
步骤一、建立视觉传感网络模型,所述的模型包括感知模型、空间模型、覆盖模型和重建误差模型;
其中,感知模型采用等腰三角形模型;该等腰三角形的参数由摄像机的参数决定:其顶点为摄像机的二维坐标(cx,cy),底边上的高为摄像机的有效感知距离dof,顶角为摄像机的视场角fov,底边上的高与x轴的夹角为摄像机的方位角
Figure FDA00027644068300000114
则底边a=2dof tan(fov/2);这样,摄像机的位姿参数记为坐标Cp=(cx,cy)和方位角
Figure FDA00027644068300000115
则一个已知类型的摄像机的状态可用参数
Figure FDA00027644068300000116
表示;将待测目标点的位置坐标设为(x,y),经过数学推导,可将该特定摄像机的感知区域用以下三个不等式表示:
Figure FDA0002764406830000011
空间模型:将目标空间
Figure FDA0002764406830000012
分为两部分:边界
Figure FDA0002764406830000013
和测量区域
Figure FDA0002764406830000014
分别对应房间的墙壁和内部空间;将目标空间内的一个点称为工作点,用Wp表示;由于只关注测量区域内的覆盖和重建问题,有
Figure FDA0002764406830000015
同理,由于传感器只能放置在边界上,有
Figure FDA0002764406830000016
将目标空间抽象为无孔洞的矩形,则边界
Figure FDA0002764406830000017
为矩形的边,测量区域
Figure FDA0002764406830000018
为矩形的内部区域;采用均匀网格采样模型,边界
Figure FDA0002764406830000019
和测量区域
Figure FDA00027644068300000110
使用不同采样间隔,分别记为Δb和Δm,同时有Δb≤Δm;同样采用均匀采样的方法将圆周分为
Figure FDA00027644068300000111
份,即采样间隔为
Figure FDA00027644068300000112
确定采样间隔后,计算目标空间边界区域在x轴和y轴上采样点个数分别为sxb,syb,测量区域在x轴和y轴上采样点个数分别为sxm,sym,系统中方位角只能取
Figure FDA00027644068300000113
个值;
覆盖模型采用普遍使用的0/1覆盖模型,处在传感器感知范围内的点状态为被覆盖,反之则为不被覆盖;
重建误差模型使用测量方程的雅可比行列式定义的几何精度因子GDOP;GDOP针对工作空间中的每个点,定义观察它的两个传感器导致的不确定性的度量;
步骤二、定义覆盖率和重建品质
1)定义覆盖率
定义以下二进制变量:
Figure FDA0002764406830000021
其中,i1∈[0,sxb-1],j1∈[0,syb-1],i2∈[0,sxm-1],j2∈[0,sym-1],
Figure FDA0002764406830000022
将以上离散值转换为连续目标空间中的坐标和角度后,
Figure FDA0002764406830000023
的值根据公式(2)计算得出,并将所有的值提前储存在矩阵中,以便随时调用;
确定参数组合
Figure FDA0002764406830000024
和一组传感器的布局方案后,可以定义二进制可见性矩阵V=(vij),其中元素为:
Figure FDA0002764406830000025
其中,任意vij的值可以通过查询式(3)得出;那么,用二进制变量bj表示点j是否满足k-覆盖:
Figure FDA0002764406830000026
则目标空间的覆盖率定义为:满足k-覆盖的采样点个数与总采样点个数之比:
Figure FDA0002764406830000027
2)定义重建品质
在二维空间中视觉传感网络中,两个视觉传感器坐标为(xi,yi),i=1,2,分别测量目标点相对于自身的方位角θi,i=1,2,而目标点位置必须是同时满足两个角度约束的点,即两条射线的交点;这里将θi定义为从x轴正半轴开始逆时针旋转过的角度,取值范围为[0,2π];ri,i=1,2为两个视觉传感器相对空间点的距离向量,且其模长为|ri|=ri,i=1,2;经过数学推导,可以得到GDOP等于目标与两个传感器之间的距离和两个传感器之间的相对方位角θ的余割的乘积,具体表达式为:
Figure FDA0002764406830000031
几何精度因子(7)表示重建过程将传感器测量误差放大为重建误差的程度,则其倒数可以表示重建结果的品质,再选择合适的比例系数dscale,将品质函数限制在[0,1],可推导出目标点j的重建品质函数:
Figure FDA0002764406830000032
上式(8)给出传感器i和传感器i′对于目标点j的重建品质;其中,vij,vi′j的表达式见式(4),乘数vijvi′j保证了当采样点
Figure FDA0002764406830000036
没有被2-覆盖时,重建品质qii′j=0;
遍历目标空间中的每个采样点,求取它们的重建品质和并求平均,即可得到整个空间的重建品质:
Figure FDA0002764406830000033
步骤三、总体设计多目标优化
本方法用空间覆盖率和空间重建品质的加权和作为综合评价传感器系统性能的量化标准,即:
e=σr+(1-σ)q         (10)
其中,r,q的表达式见公式(6)和公式(9),加权系数σ∈[0,1]表示评价标准倾向于覆盖率的程度;当σ=1时,评价标准为单一的最大化覆盖率;σ=0时,则为最大化重建品质,即最小化重建误差;
综合上述分析,二维优化问题可整理为以下模型:
给定目标区域长和宽,采样参数组合
Figure FDA0002764406830000034
传感器数量N,优化以下问题:
Figure FDA0002764406830000035
其中,
Figure FDA0002764406830000041
G=r,q,or e;
针对离散优化问题(11),采用遗传算法和粒子群算法的组合算法来求解。
2.根据权利要求1所述的一种基于多目标约束的二维空间视觉传感器布局优化方法,其特征在于:所述的遗传算法和粒子群算法的组合算法,流程如下:
a)求取初始种群:将最大迭代次数设为50,运行粒子群算法,得到初步优化后的种群P;
b)筛选初始种群:保留种群P中最佳的9个个体,以及整个优化过程中的最佳个体,共10个个体作为部分初始种群,不足的部分则随机生成,两部分共同构成初始种群P(0);
c)迭代:以P(0)为初始种群,运行遗传算法;
d)判断终止条件:达到终止条件后,整个运算过程中适应度最大的个体为最优解,终止算法。
3.根据权利要求1所述的一种基于多目标约束的二维空间视觉传感器布局优化方法,其特征在于:所述步骤三具体的优化算法实施过程中,为简化优化参数和约束条件,进一步提出了一种用一维参数表示传感器位置信息的降维搜索方法;若将直角坐标系原点设为起点,逆时针沿目标空间的边界运动,则边界上的每个离散点唯一对应一个运动的前向距离d;以此前向距离代表传感器的位置信息,即可将二维参数(cx,cy)降为一维参数d;这样,针对有N个传感器的二维网络,其中第i个传感器的位姿信息由参数组合
Figure FDA0002764406830000042
简化为
Figure FDA0002764406830000043
整个网络的决策变量即为
Figure FDA0002764406830000044
CN202011228559.1A 2020-11-06 2020-11-06 一种基于多目标约束的二维空间视觉传感器布局优化方法 Active CN114329854B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011228559.1A CN114329854B (zh) 2020-11-06 2020-11-06 一种基于多目标约束的二维空间视觉传感器布局优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011228559.1A CN114329854B (zh) 2020-11-06 2020-11-06 一种基于多目标约束的二维空间视觉传感器布局优化方法

Publications (2)

Publication Number Publication Date
CN114329854A CN114329854A (zh) 2022-04-12
CN114329854B true CN114329854B (zh) 2023-05-12

Family

ID=81031899

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011228559.1A Active CN114329854B (zh) 2020-11-06 2020-11-06 一种基于多目标约束的二维空间视觉传感器布局优化方法

Country Status (1)

Country Link
CN (1) CN114329854B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116894400B (zh) * 2023-09-11 2023-12-12 南京邮电大学 一种基于视觉传感网络的空间组合选址方法
CN117010132B (zh) * 2023-09-27 2023-12-26 中国船舶集团有限公司第七一九研究所 一种水下多基地声纳系统的空间阵位优化方法及系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102722103A (zh) * 2012-05-18 2012-10-10 镇江雅迅软件有限责任公司 基于模糊分析和可视化技术的多目标参数寻优方法
CN103630072A (zh) * 2013-10-25 2014-03-12 大连理工大学 双目视觉测量系统中摄像机的布局优化方法
CN105069240A (zh) * 2015-08-17 2015-11-18 湖北工业大学 一种空间测量定位系统测站布局智能优化方法
CN107066747A (zh) * 2017-04-24 2017-08-18 哈尔滨理工大学 一种视觉测量网络组网规划方法
WO2020000248A1 (zh) * 2018-06-27 2020-01-02 大连理工大学 一种基于空间重构的航空发动机过渡态加速过程关键性能参数预测方法
CN111460597A (zh) * 2020-04-02 2020-07-28 电子科技大学 一种基于改进的多目标粒子群优化算法的雷达布站方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102722103A (zh) * 2012-05-18 2012-10-10 镇江雅迅软件有限责任公司 基于模糊分析和可视化技术的多目标参数寻优方法
CN103630072A (zh) * 2013-10-25 2014-03-12 大连理工大学 双目视觉测量系统中摄像机的布局优化方法
CN105069240A (zh) * 2015-08-17 2015-11-18 湖北工业大学 一种空间测量定位系统测站布局智能优化方法
CN107066747A (zh) * 2017-04-24 2017-08-18 哈尔滨理工大学 一种视觉测量网络组网规划方法
WO2020000248A1 (zh) * 2018-06-27 2020-01-02 大连理工大学 一种基于空间重构的航空发动机过渡态加速过程关键性能参数预测方法
CN111460597A (zh) * 2020-04-02 2020-07-28 电子科技大学 一种基于改进的多目标粒子群优化算法的雷达布站方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
乔玉晶 ; 谭世征 ; 姜金刚 ; .一种多视觉测量组网规划策略.光学学报.2018,(第05期),全文. *
宋婷婷 ; 张达敏 ; 王依柔 ; 徐航 ; 樊英 ; 王栎桥 ; .基于改进鲸鱼优化算法的WSN覆盖优化.传感技术学报.2020,(第03期),全文. *

Also Published As

Publication number Publication date
CN114329854A (zh) 2022-04-12

Similar Documents

Publication Publication Date Title
CN107179768B (zh) 一种障碍物识别方法及装置
CN114329854B (zh) 一种基于多目标约束的二维空间视觉传感器布局优化方法
CN111429574A (zh) 基于三维点云和视觉融合的移动机器人定位方法和系统
CN103984981B (zh) 基于高斯过程模型的建筑物环境传感器测点优化方法
US9111172B2 (en) Information processing device, information processing method, and program
Fatlawi et al. Optimal camera placement using sine-cosine algorithm
CN115375868B (zh) 地图显示和遥感地图显示方法、计算设备以及存储介质
CN115563732B (zh) 基于虚拟现实的喷涂轨迹仿真优化方法及装置
CN116824396A (zh) 一种多卫星数据融合自动化解译的方法
CN114329855B (zh) 一种无线视觉传感网络的传感器布局优化与快速部署方法
US20070192068A1 (en) Method For Planning A Security Array Of Sensor Units
CN115983478A (zh) 分布式光伏发电功率预测分析方法、系统、终端及介质
CN116341231A (zh) 人工增雨催化剂扩散轨迹的模拟方法、装置、服务器和存储介质
CN111413679B (zh) 一种计算无人机蜂群rcs的方法
CN115018983A (zh) 移相变压器选址方法、装置、电子设备及存储介质
CN111489440B (zh) 对非标零件的三维扫描建模方法
CN112929818A (zh) 一种基于卡尔曼滤波和泛克里金插值的室内定位方法
CN108205683B (zh) 基于可见光卫星图片的红外目标识别系统阈值优化方法
CN117689846B (zh) 线状目标的无人机摄影重建多交向视点生成方法及装置
Kim et al. Optimal Pre-processing of Laser Scanning Data for Indoor Scene Analysis and 3D Reconstruction of Building Models
CN117745977B (zh) 基于3dgis的山地自然风场分析方法、系统及终端
CN116834003B (zh) 光伏组件的智能安装方法及系统
Ivanov et al. Review of Drone Swarms Usage for 3D Reconstruction
CN117058358B (zh) 一种场景边界检测方法和移动平台
CN117521172A (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