专利名称：反射层析激光雷达双折线反投影成像方法
技术领域：
本发明涉及反射层析激光雷达，特别是一种反射层析激光雷达双折线反投影成 像方法，通过把反投影波面近似为双折线，将所测得的球面反射投影信号反投影到 对应的双折线上，用双折线反投影信号完成相应距离和质量要求的图像重建。
背景技术：
在医学计算机辅助层析(简称为CT)中，根据一维透射投影信息实现两维图像 重建已经得到广泛的应用和发展。所得到的一维透射投影信息反应的是在多个不同 角度下，物体内部不同部位透射系数的差异。反投影成像的内容为重建图像某一点 的光强密度值可看作这一平面内所有经过该点的射线投影值之和。
如图l透射层析所示，/(x,力为代建图像，丄^为对应的直线^xcos^ + ysin^，
P0"一)为/0,力沿直线丄^的透射特征系数的积分，即角度^所对应的一维透射投 影
/"力A (1)
用反投影成像方法重建图像g(x,:v)，贝(L
g(X力=￡ P" cos《+ y sin《.A) A-
'=i (2)
其中，^为第z'次投影所对应的角度，A^为投影角度采样间隔，w为总的投影 个数。CT反投影成像把取自有限物体空间的射线均匀地反投影到射线所及的无限空 间的各点之上，使原先像素数值为零的点经反投影之后不再为零，产生星状伪迹。 为了减少重建图像星状伪迹，引入了滤波反投影算法，对获得的投影进行适度的调 整，调整之后的投影数据为-
《(r^)-i^(IA:l巧Cp(r^))) (3)
其中，符号尸i、 i^分别表示一维的傅里叶变换与傅里叶反变换，"为空间变量， A:为频率变量，则经过滤波反投影重建图像g^(义，力表达式如下<formula>formula see original document page 5</formula>
反射层析(基于距离)类似于上述CT滤波反投影重建图像的操作过程，其实质 是将到源点某一距离的有限表面区域反射所得的投影信息，反投影到隔源点相同距 离的所有区域，如图2反射层析激光雷达所示。所不同的是CT依据的透射特征系数， 只需要0-;r的透射投影信息，而反射层析依据的是反射特征系数，需要O-2;r的反 射投影信息。
随着激光的出现和发展，Charles L. Matson (参见文献R. M. Marino, R. N. Capes， W. E. Keicher et al… Tomographic image reconstruction from laser radar reflective projections [J] ， SPIE Laser Radar 111,1988， Vol. 999:248—263)， Parker. J. K.(参见文献Jeffrey K. Parker, E.B.Craig, D.I. Klick et al.. Reflective tomography:image from range resolved laser radar raearsurements[J], Applied Optics, 1988'27 (13) : 2642—2643) ， F. K. Knight (参 见文献Frederick L Knight, David Klick， Danette P. Ryan-Howaxd et al. . Laser radar reflective tomography utilizing a steak camera for precise range resolution [J]. Applied Optics, 1989, 28(12) :2196—2198)等很快的将反射层析 算法引入到激光雷达中。根据激光短脉冲探测和非相干接收器得到基于距离的一维 反射投影信息，激光雷达系统可以实现二维图像重建，具有很高的分辨率。当波阵 面发生弯曲时，如图3反射层析激光雷达所示，直接的球面反投影滤波方式计算量 大，需要的处理时间长。引入双折线的处理方法，在一定距离下，可以减少计算量, 实现图像的重建。
发明内容
本发明的目的在于为反射层析激光成像雷达提供一种反射层析激光雷达双折线 反投影成像方法，该方法的原理可靠，结构简单，易于实现。
本发明的具体技术解决方案如下
一种反射层析激光雷达双折线反投影成像方法，其特点是包括以下步骤 ①由探测物体的目标尺寸X皿、以探测器和光脉冲二者中较低的分辨率作为系 统的距离分辨率M、激光脉冲宽度r三个因素，根据奈奎斯特采样定律计算满足物 体完全重建所需要的最小投影角度采样间隔A^和单个投影最少采样点数iVp:A , 360 180A/
② 激光器发出光脉冲，用会聚透镜使光脉冲会聚后再次发散，通过调节透镜的 焦距来改变光斑的大小，使到达探测目标的光斑外轮廓涵盖整个探测目标，获得需
要光束最小发散角"^arcsinO,/2丄)，也就是所对应的双折线反投影的角度，其中 L是激光器与待测目标之间的距离根据激光器发出光脉冲的时刻、，最早接收到反射 信息的时刻^，最晚接收到反射信息的时刻,3，探测器到待测目标的距离丄为
铜),2一"，其中c为光速； 2
③ 将初始待测目标所对应的角度定义为^ = 0° ，激光器发出单脉冲照射到待测 目标上，探测器记录在角度^=0°下所对应的经过目标反射加宽后的反射投影信 息，即探测器所得到的光功率随时间的变化信息反映的是物体表面反射特征系数随
物体深度的变化信息
其中，/。(x,力为物体表面的反射特征系数，e代表到达物体表面的均匀光强， 户。"A)为角度A下获得的反射投影信息，^为物体表面积分面元；
④ 待测目标旋转角度A^后，探测角度变为《=^ +厶6激光器再次发出单个脉
冲，探测器记录相应的反射投影信息A(f,A);
⑤ 判断投影角度4是否大于36(r,若小于36(T，重复步骤④(i一l)次操作，待 测物体旋转込^， i为2以上的正整数，相应记录反射投影信息A"《)，直至《大
于360。；
⑥ 滤波处理利用公式r-"， C为光速，将所得到的随时间变化的投影信息 A《,《)，转换成随距离变化的投影信息，对每个角度获得投影信息进行滤波处理，
减少重建图像中的伪迹，具体滤波公式为^.(rW.hi^dwl^O^r,^))),其中"为
频域变量，符号《、巧—i分别表示信号一维傅里叶变换和傅里叶反变换；
⑦ 建立重建图像矩阵g^(x,力，其像元坐标间隔为所对应的系统距离分辨率
6A/ ，为获得完整的目标重建图像，图像矩阵g^(x，力行数、列数要大于，
并将矩阵所有像元初始值设为零；
⑧ 取图像矩阵g^(x，力第 一个像元所对应的点坐标(X，力)，按照角度
砍(/ = 0， 1,2…进行旋转，旋转之后的新坐标变为
x; = x！ cos(《.)+ 乂 sin(《),乂 = —& sin(《)+ k cos(《.)
判断新坐标点Ocp乂)到角度《、^所对应的双折线的距离是否满足小于或等于 像素间隔A/ 的V^/2倍，若满足以上条件，则图像矩阵元g刚(x"乂)累加上反射投 影值《,化.，A)A0 ，直至像元点坐标(^1,;;1)完成所有旋转角度^(/ = 0,1，2 7^)和距离 r(^0,l，2…Ag的计算。
具体重建公式如下
新坐标点(x;,力)到角度《、。所对应的双折线的距离" ^ I _x' sin a+1 _y' I (1 - cos <2) + - 。) sin a |
若"卓m ，则g刚",乃)=f;i>'(0,《)~
⑨ 将图像矩阵g^0c,"的所有其它像元，依次重复步骤⑧操作，最终将得到双
折线滤波反投影重建图像g^ 。
本发明的技术效果-
本文在充分考虑激光发散角，波阵面呈现球面的状态。为此我们对反射层析激 光雷达成像采用双折线滤波反投影的方法，将球面反射投影信息反投影到双折线上， 利用双折线信息实现图像的重建。改善了直接球面反投影计算量大、处理时间长的 问题。
图l透射层析示意图； 图2反射层析示意图； 图3反射层析波阵面弯曲示意图；图4反射层析反射投影探测示意图中l激光器；2会聚透镜；3探测目标；4光纤耦合头；5光纤；6探测器。
图5原正三角形图样，i 为物体所占区域半径；
图6双折线反投影原理示意图7双折线反投影重建图像二值化；
具体实施例方式
下面结合附图和实施例进一步详细说明本发明，但不应以此限制本发明的保护 范围。
一种反射层析激光雷达双折线反投影成像方法，包括以下步骤
(1) 由探测物体的目标尺寸、探测器距离分辨率、激光脉冲宽度r三个因素， 根据奈奎斯特采样定律计算满足物体完全重建所需要的最小投影角度采样间隔、 单个投影最少采样点数A^;
在实际选择中，距离分辨率要以探测器和光脉冲二者中较低的分辨率作为系统 距离分辨率，否则会造成信息的丢失。例如，所用探测器带宽为4GHz，采样率最大 为20GHz的示波器，脉冲激光器发出的脉冲宽度为8ps，可以计算示波器的距离分辨 率(7.5mm)低于光脉冲所对应的距离分辨率(1.2ram)，应以示波器分辨率作为系统
距离分辨率。假设系统距离分辨率为Ai ，转换到频域Fr化e"(工》sinc(Ai ")
其中，符号Fr表示傅里叶变换，符号"c/表示矩形函数，符号sinc表示sinc函 数。定义sinc函数第一个零点为图像频域截止频率^M: Mmax=l/Ai 。假设目标物
体最大尺寸为&a,应用标准的奈奎斯特采样定律，要想完全重建物体图像，傅里 叶频域采样间隔要小于l/j^^ ，可以得到反射层析旋转角度采样数
乂"訂皿"皿=2"^，投影角度采样间隔A「^ = l^。
每个投影上采样点数-
(2) 如图4所示，脉冲激光器发出光脉冲，用会聚透镜使光脉冲会聚后再次发 散，可以通过调节透镜焦距来改变光斑的大小，使到达目标物体的光斑外轮廓涵盖 整个探测目标；若物体距离脉冲激光器足够远，也可以不用透镜，光束的发散角导致远场光斑覆盖整个物体即可；假若探测目标最大尺寸为^a，双折线反投影所对 应的角度为c^arcsin(x^/2丄)，L为探测器到物体的距离，根据激光器发出光脉冲 的时刻^，最早接收到反射信息的时刻,2，最晚接收到反射信息的时刻G，探测器到
物体的距离￡为4"2+。/2_"，其中c为光速；
2
(3) 将初始测量目标所对应的角度定义为0。=0° ，激光器发出单脉冲照射到
目标物体上，探测器记录在角度0。=0°下，所对应的经过目标反射加宽后的反射投
影信息。在距离比较远的情况下，到达物体表面波面光强是均匀的，那么反射回来 的光功率信息的变化则是物体表面的反射特征系数引起的，即探测器所得到的光功
率随时间的变化信息反映的是物体表面反射特征系数随物体深度的变化信息
其中，/。(x,力为物体表面的反射特征系数，e代表到达物体表面的均匀光强，
A^，A)为角度A下获得的反射投影信息，^为物体表面积分面元；
(4) 当目标物体发生旋转角度/V)后，探测角度变为^ + Z^，激光器再次发出 单个脉冲，在角度^+Z^下探测器记录相应的反射投影信息A"《)，并更新当前角
(5) 判断投影角度^是否大于36(T，若小于360',多次重复步骤(4)操作，当 物体旋转込#' = 1，2,3……)时重新记录所对应的反射投影信息/7,.(/,论)，并更新角度 《=步。+ ，直至《大于360°;
(6) 利用公式^" = " (c为光速)，将所得到的随时间变化的投影信息，转换成 随距离变化的投影信息。对每个角度获得投影信息进行滤波处理，减少重建图像中
的伪迹，具体滤波公式为仏(r4.X(i 其中"为频域变量，符号
《、F「1分别表示信号一维傅里叶变换及傅里叶反变换；
(7) 建立图像g^(^,力，其像元坐标间隔为所对应的系统距离分辨率M。为
《
M
表示取整数上限，例如「2.2， = 3)，并将矩阵所有像元初始值设为零；
获得完整的目标重建图像，图像矩阵g^(x，力行数、列数要大于
(符号「]
(8)取图像矩阵g^O,力第一个像元所对应的点坐标(^乃)，按照角度
9A(Z = 0,1,2… )进行旋转，旋转之后的新坐标变为
x; = x！ cos(A) + x sin(《.),乂 = sin(^) +乃cos(《)
判断新坐标点(x;,力)到角度^ 、 ^所对应的双折线的距离是否满足小于或等于
像素间隔m的V^/2倍，若满足以上条件，则图像矩阵元g^(^，乃)累加上反射投 影值《化.，直至像元点坐标0cpk)完成所有旋转角度= 0， 1，2… )和距离 Kz、0,l,2…a^)的计算。 具体重建公式如下
新坐标点oc;，力)到角度《、g所对应的双折线的距离d:
I 一x' sin a+1 I (1 - cos a) + (i — ) sin |
<i =-, =-^-
>/2 — 2cos
若"幸m，则g则(^乂)-f;f;化,《)A^ (7)
(9)将图像矩阵g^O,y)的所有其它像元，依次重复步骤(8)操作，最终将 得到双折线滤波反投影重建图像g^ 。
计算机仿真采用边长为3m的正三角形(如图5所示)，距离分辨率采样m^. lm， 则物体完全重建所需要的投影角度间隔Av/《3.3度，A^2 34.64。模拟中Av/4度， 单投影采样点数iVp取满足采样条件的最小整数35，距离L-10米。依次按照具体实
施步骤进行计算机仿真，重建图像二值化结果如图7所示，可以证明本发明采用双 折线滤波反投影的方法对反射层析激光雷达进行图像重建是一种行之有效的办法。 计算量的减少可以直观的通过matlab运行时间来看，在同一外界条件下，仿真双折 线反投影需要2分16秒，而球面法需要2分49秒。故双折线滤波反投影可以在减 少计算量的前提下，实现相应距离、相应图像质量要求的图像重建。
权利要求
1、一种反射层析激光雷达双折线反投影成像方法，其特征在于包括以下步骤①由探测物体的目标尺寸xmax、以探测器和光脉冲二者中较低的分辨率作为系统的距离分辨率ΔR、激光脉冲宽度τ三个因素，根据奈奎斯特采样定律计算满足物体完全重建所需要的最小投影角度采样间隔Δφ和单个投影最少采样点数Np②激光器(1)发出光脉冲，用会聚透镜(2)使光脉冲会聚后再次发散，通过调节透镜的焦距来改变光斑的大小，使到达探测目标(3)的光斑外轮廓涵盖整个探测目标，获得需要光束最小发散角α＝arcsin(xmax/2L)，也就是所对应的双折线反投影的角度，其中L是激光器与待测目标(3)之间的距离根据激光器发出光脉冲的时刻t1，最早接收到反射信息的时刻t2，最晚接收到反射信息的时刻t3，探测器(6)到待测目标(3)的距离L为其中c为光速；③将初始待测目标所对应的角度定义为φ0＝0°，激光器(1)发出单脉冲照射到待测目标(3)上，探测器(6)记录在角度φ0＝0°下所对应的经过目标反射加宽后的反射投影信息，即探测器所得到的光功率随时间的变化信息反映的是物体表面反射特征系数随物体深度的变化信息p0(t，φ0)＝e∫f0(x，y)ds其中，f0(x，y)为物体表面的反射特征系数，e代表到达物体表面的均匀光强，p0(t，φ0)为角度φ0下获得的反射投影信息，ds为物体表面积分面元；④待测目标(3)发生旋转角度Δφ后，探测角度变为φ1＝φ0+Δφ，激光器(1)再次发出单个脉冲，探测器(6)记录相应的反射投影信息p1(t，φ1)；⑤判断投影角度φ1是否大于360°，若小于360°，重复步骤④i—1次操作，待测物体旋转iΔφ，i为2以上的正整数，相应记录反射投影信息pi(t，φi)，直至φi大于360°；⑥滤波处理利用公式r＝ct，c为光速，将所得到的随时间变化的投影信息pi(t，φi)，转换成随距离变化的投影信息，对每个角度获得投影信息进行滤波处理，减少重建图像中的伪迹，具体滤波公式为其中u为频域变量，符号F1、分别表示信号一维傅里叶变换和傅里叶反变换；⑦建立重建图像矩阵gFB1(x，y)，其像元坐标间隔为所对应的系统距离分辨率ΔR，为获得完整的目标重建图像，图像矩阵gFB1(x，y)行数、列数要大于并将矩阵所有像元初始值设为零；⑧取图像矩阵gFB1(x，y)第一个像元所对应的点坐标(x1，y1)，按照角度φi(i＝0，1，2…Nφ)进行旋转，旋转之后的新坐标变为判断新坐标点到角度φi、rj所对应的双折线的距离是否满足小于或等于像素间隔ΔR的倍，若满足以上条件，则图像矩阵元gFB1(x1，y1)累加上反射投影值qi(rj，φi)Δφ，直至像元点坐标(x1，y1)完成所有旋转角度φi(i＝0，1，2…Nφ)和距离r(i＝0，1，2…Np)的计算，具体重建公式如下新坐标点到角度φi、rj所对应的双折线的距离d若则⑨将图像矩阵gFB1(x，y)的所有其它像元，依次重复步骤⑧操作，最终将得到双折线滤波反投影重建图像gFB1。
全文摘要
一种反射层析激光雷达双折线反投影成像方法，实质是将反投影波阵面近似为双折线的方法，将激光雷达探测器所接收到的球面反射投影信号反投影到所对应的双折线上，通过这些双折线信号重建图像。本发明方法与球面直接反投影相比，双折线反投影算法计算量小，操作简单，易于实现。模拟仿真表明本发明方法可以替代球面反投影完成相应质量要求的图像重建。
文档编号G01S17/89GK101545976SQ20091005085
公开日2009年9月30日 申请日期2009年5月8日 优先权日2009年5月8日
发明者毅 严, 刘立人, 吴亚鹏, 煜 周, 孙建锋, 金晓峰 申请人:中国科学院上海光学精密机械研究所