计算机仿真论文(计算机仿真 论文)水下光电成像系统计算机仿真

1、计算机仿真论文（计算机仿真 论文）：水下光电成像系统计算机仿真摘要:为进行水下目标探测,该文针对Jaffe模型的不足,利用小角度近似原理,建立了水下光电成像系统的理论模型。考虑散射光的影响,给出了给定距离处的水下成像面各像素的信号辐射量的计算方法和对比度模型。对常规激光助视水下成像探测系统进行了仿真,研究了光源系统和接受系统的距离与成像对比度的关系。与Jaffe模型相比,该文所建模型中变量的物理意义更加明确,而仿真结果与其类似,证明该文所述方法的正确性。关键词:水下成像;小角度近似;对比度模型;计算机仿真Computer Simulation forUnderwater Photoelectr

2、ic Imaging SystemAbstract:In view ofthe shortcomingsofthe Jaffemode,l a theoreticalmodelofunderwaterphotoelectricmi aging system is established forunderwater targetdetection using the small angle approxmi ation princi-ple. Considering the influence ofscatted ligh,t the calculationmethod ofsignal rad

3、iation quantity and thecontrastmodel ofpixels in underwater mi agingplane ofa given distance isproposed. An underwater mi a-ging detection system with extra light source is smi ulated. The relationship of the distances between thelightsource system and the acceptance system and the mi age contrast i

4、s researched. Comparedwith theJaffemode,l the physicalmeaning of the variables in the proposedmodel is clearer, and the smi ulationresults are smi ilar, which proves the correctness of the proposedmethod.Keywords:underwater imaging; small angle approximation; contrastmodels; computer simulation水下光电成

5、像探测系统在海洋探测、水下兵器试验等方面有着越来越重要的应用1, 2。水下光电成像探测系统的主要类型包括:摄像机-光源紧凑型、摄像机-光源分开型、距离门选通型和激光同步扫描型。这些不同类型的水下成像探测系统在不同的场合均有不同的应用。符合实际的水下成像系统模型有广泛的应用,它可以为系统设计、后续图像处理算法研究、水下图像识别等提供重要的指导3;也可以为海水光学特性研究和测量提供更加便利和快捷的方法,例如人们研究发现海水的相位函数和传递函数之间存在傅里叶变换的关系,从而为相位函数的计算提供了便利4。研究人员对水下光电成像探测系统的系统建模作了大量工作:Tan等人5研究了短距离工作时的距离门选通水

6、下光电成像探测系统的计算模型,黄有为等人6研究了基于光束空间展宽的水下前向散射成像模型, Jaffe7分析了同步扫描水下成像的性能极限,Rampolla基于后向散射效应给出了激光助视光电成像系统的理论模型8。在这些模型中,除Jaffe外,均没有考虑整个系统的成像过程,虽然Jaffe考虑小角度散射原理所建立的模型最为贴近实际,但在其模型中,并没有给出各种参量的具体计算方法。本文针对Jaffe模型的不足利用小角度近似原理,建立了水下光电成像系统详细的数学计算模型,给出了仿真步骤,得到仿真结果。与Jaffe模型相比,本文所建模型中各项的物理意义更加明确,而仿真结果与其模型类似,表明本文所述方法的正确

7、性。1信号辐射计算模型1. 1接受系统的信号分类水下光学成像系统的建模反映了从发射系统发出的光线到达接受系统的过程。如图1所示,从光源发出的光可以有三种到达接受系统的途径:直接传输、前向散射和后向散射。光源发出的光直接经目标反射到达接收系统称为直接传输;经水的散射再经目标发射到达接收系统称为前向散射,这部分光称为前向散射光;没有经过目标发射而是经过水的散射进入接收系统称为后向散射,这部分光称为后向散射光。假设成像系统接收的总辐射光通量ET是这3个元件的线性组合ET=Ed+Efs+Ebs式中:Ebs为后向散射;Efs为前向散射;Ed为直接传输。1.2系统结构模型考虑发射和接收系统间的距离,成像系

8、统框图如图2所示,假定接收系统与发射系统之间的间距为d。距接收系统Rm远处与接收系统全视场对应的物面宽度为Wm、高度为Hm。以物面中心(Wm/2,Hm/2)处为视场坐标原点,点(m,n)处极限可分辨单元的角坐标为(R,R),对应的直角坐标为(H,W),已知物面处的标准像素间距为,则它们之间有如下关系m=Hm2-H /n=Wm2-W /R=arctanWRmR=arctanH(W2+R2m)1/2(1)在该物面与接收系统之间的辐射传播空间内, (m,n)成像主光线路径上任选的至接收系统OR距离为Ra的空间点(R,R,Ra)相对于发射系统Ol的极坐标为(1,1,R1),其值可由式(2)决定1=ar

9、ctan tanR+dRacosRcosR1=arctancos1tanRcosRR1=RasinR/sin1(2)1.3信号和前向散射辐射量计算为计算景物面的直接反射量,必须首先计算到物面的照射光强。假设物面(m,n)相对于发射系统O1的极坐标为(s,s,Rs),激光束沿着主轴Rm方向照射,其强度为I,在考虑的物面单元(m,n)方向的辐射强度为I(s,s),经光强被衰减而且球面展开,有EI(m,n,s,s)=I(s,s)cosexp(-Rs) /R2s(3)397南京理工大学学报(自然科学版)第34卷第3期为O1与物点(m, n)的连线与过点(m, n)且垂直于物面的直线的夹角。显见,根据图

10、2,Rs和的大小为Rs=H2+R2m+(d+W)2(4)=arccosRmH2+R2m+(d+W)2(5)对反射光线而言,点(m,n)可视为开始点,入射辐射光通量用EI(m,n,0)表示,且有EI(m,n,0)=EI(m, n,s,s)。为精确考虑入射光通量,必须考虑由于小角度前向散射部分形成的照射,这可以通过引入点扩散函数g得到EI(m,n,0)=EI(m,n,0)g(m,n|Rs,G,c,B)+EI(m,n,0) (6)式中:g(m,n|Rs,G, c,B)=exp(-GRs)-exp(-Rs)F-1exp(-BRsf) (7)G是实验常数|G|,B是实验衰减因子,f是空间频率。现计算反射

11、辐射,它等于入射辐射乘以点(m,n)处的反射率(m,n)1,典型海洋地图中0.02(m,n)0.1。同时考虑郎伯定律,如果出射光线与点(m,n)的连线与穿过该点垂直于物面的直线夹角为,必须引入1个几何因子cos。考虑光学系统的孔径、反射波的球形展开,得到Ed(m,n)=EI(m,n,0)exp(-RR)(m,n)4f#2Tlcos4(8)式中:f#是接收系统的F数,Tl是透镜的透过率。显见,根据图2,cos的计算公式如下cos=RmR2m+W2+H2(9)将直接辐射量与点扩散函数g相卷积,即可得到前向散射辐射量Efs(m,n)=Ed(m,n,0)g(m,n|RR,G,B)(10)1.4后向散射

12、辐射量计算由于后向散射从大角度范围内进入接收系统,所以后向散射是水下成像系统模型计算中最费时的一部分。目前通用的方法是将三维空间切成平行于接收系统像面厚度为Zi的一系列平面,分别计算入射到这些平面的辐射量,方法与前面一样,即分别计算直接辐射和小角度散射引起的附加辐射。然后计算由于分割的切片对入射光向着接收系统方向散射(即后向散射)而到达接受系统的辐射量。从发射系统发出的光传输到三维空间的入射光通量用Es(x,y, z)表示,根据以上讨论,入射到体积元的入射量是直接投射和前向散射之和Es(m,n,R1)=Es, d(m,n,R1)+Es, fs(m,n,R1)(11)Es, d(m,n,R1)=

13、I(1,1)exp(-R1)R21(12)Efs(m,n,R1)=Es, d(m,n,R1)g(m,n|R1,G,B) (13)由于散射作用,切片将部分入射光散射转换到接受系统方向的信号通路中,具体转换大小由介质的内在特性体积散射函数决定。在体积V的水域对入射光向着接受系统单元(x,y)散射的辐射强度为Hbs(,x,y)=(bs)Es(m,n,R1)V(14)式中:()为体积散射函数,V是三维空间中一个递增的体积,bs是计算的体积元和光源的连线与体积元和接受系统的连线的夹角。(bs)=(R-1,R-1)表示从(1,1)方向的入射光经切片单元散射转向(R,R)方向从而进入信号通道的能力。最后对所

14、有的切片形成的辐射量进行相加,并考虑光学系统的影响,后向散射的直接部分为Ebs, d(m,n) =Ni=1exp(-cRai)(bs)Es(m,n,Rli)Rai4f#2cos3Tl(15)式中:Rai是后向散射体积V的厚度,Rai是接受系统中某个像素到水体切片中心的距离,N是切片的总个数,同样考虑后向散射光向接受系统传输时前向散射的影响,由于后向散射引起的总辐射量为Ebs(m,n)=Ebs, d(m,n)+Ebs, d(m,n)g(m,n|RR,G, c,B) (16)2计算机仿真分析2.1对比度计算模型背向与前向散射导致成像对比度恶化,相邻两单元,原始对比度为C0=A-BA=EdA-EdB

15、EdA(17)398总第172期孔捷张保民水下光电成像系统计算机仿真由后向散射与前向散射对成像辐射的叠加,使成像对比度衰降C=EdA-EdBEdA+Ebs+Efs=C0(18)式中:=(1+)-1。背向与前向散射造成成像对比衰降的恶化系数=(Ebs+Efs)EdA(19)为评价散射对对比度的影响,采用评价参数对比度透过率C(m, n),它是直接光分量与总的光通量的比值C(m,n)=Ed(m,n)ET(m,n)(20)2.2信号光计算步骤(1)假设探测目标到接受系统的垂直距离为Rm。接受系统的探测器单元像素尺寸为dxdy,光学系统的焦距为f,从而可计算出物面处的标准像素间距为和物面第(m, n)

16、个探测单元的位置为x=Rmfdxy=RmfdyH=mxW=my(21)(2)利用式(4)和(5)计算Rs和,再利用式(4)计算入射到物面单元(m,n)的光通量;(3)在计算出所有物面单元的入射光通量分布后,考虑光波传输中的前向散射并利用小角度散射原理,利用式(6)和(7)计算最终的入射光通量分布;(4)利用式(9)、(10)和(11)计算最终接受系统得到的信号光通量分布。2.3后向散射光计算步骤参见图2,本文对某个接受单元的后向散射光的计算主要针对其对应接受光主轴进行,即沿着连接物面单元点(m,n)和接受系统OR的连线所在区域进行,在实际系统中,还应考虑其它方向的散射光。具体算法如下:(1)将

17、物面单元点(m, n)和接受系统OR的连线所在区域分成M等份,编号从接近物面单元点(m,n)的区域开始,并利用公式(1)(3)计算每个区域相对发射系统Ol的角坐标。(2)根据发射系统发射光的最大偏角光线和物面单元点(m, n)与接受系统OR的连线的交点,决定散射的最大区域编号,设为N。(3)选定1个体积相位函数模型,利用公式(12)(17)计算整个区域散射到像面的光通量。2.4介质的相位函数为了计算后向散射光,必须首先确定相位函数。相位函数是描述1个散射介质对入射光向各个方向进行散射的能力。目前已有多种相位函数模型,其中,比较通用的是Hengey-Greenstein(OTHG)相位函数,它是

18、1个单参数函数,而且数学表示也比较简单OTHG()=141-2(1+2-2cos)3/2(22)式中:0,1,称为前向散射参数。对上述入射光方向为(l,l),出射光方向为(R,R),有cos=+(1-)(1-)cos(l-R) (23)式中:=cosl,=cosR。2.5仿真分析仿真工作主要针对一般激光助视水下成像系统,主要研究激光发射系统和接受系统的距离与对比度的关系。假设光源系统与接受系统平行放置,为简化计算,且假设光源向各个方向发出的光强相同。仿真中,前向散射系数为0.5。仿真结果如图3所示,可知:两种方法得到的结果类似;适当拉开光源系统和接受系统的距离,有助于降低散射光特别是后向散射光

19、对系统对比度的影响,从而提高系统的成像性能。从物理结构上这也容易解释,拉开距离意味着增大了入射和散射角,有利于降低散射作用,同时由于激光主轴仍然对着探测物面,因此信号光并没有受到太大的影响。图3对比度与发射和接收系统间距离的关系仿真3结束语利用小角度散射近似,考虑水下激光成像系统的成像过程,建立了水下成像探测系统理论模型。考虑散射光的影响,给出了给定距离处的水399南京理工大学学报(自然科学版)第34卷第3期下成像面各像素的信号辐射量的计算方法和对比度计算模型。最后对在实际应用中比较普遍的常规激光助视水下成像探测系统进行了仿真,研究了光源系统和接受系统的距离与成像对比度的关系。结果表明这些模型的建立可以为水下成像探测系统的设计和优化提供指导。下一步的工作将采用该模型,计算距离门选通激光成像系统的成像情景。参考文献:1徐啟阳,杨坤涛,王新兵,等.蓝绿激光雷达海洋探测M.北京:国防工业出版社, 2002.2张毅,柏连发,陈钱,等.布里渊散射在水下激光成像探测中的应用 J.红外与激光工程, 2009,38(4): 665-668.3 Jaffe JS. Computermodelin