反射投影成像投影图中心的对准方法.pdf

反射投影成像投影图中心的对准方法
    【技术领域】

    本发明涉及反投影成像系统，特别是一种反射投影成像投影图中心的对准方法。在反投影成像系统中，存在目标物体平动所产生的旋转中心偏移问题，最终将导致不同角度下反射投影中心失配，重建图像发生错位。本发明方法利用特征点的投影信息可以迅速快捷的实现不同角度反射投影中心的校准，操作简单，易于实现。

    背景技术

    在医学计算机辅助层析(简称为CT)中，根据一维透射投影信息实现两维图像重建已经得到广泛的应用和发展。所得到的一维透射投影信息反应的是在多个不同角度下，物体内部不同部位透射系数的差异。反投影成像的内容为重建图像某一点的光强密度值可看作这一平面内所有经过该点的射线投影值之和。

    如图1透射层析所示，f(x，y)为待建图像，L_r，φ为对应的直线r＝xcosφ+ysinφ，p(r，φ)为f(x，y)沿直线L_r，φ的透射特征系数的积分，即角度φ所对应的一维透射投影：

     p(r,φ)=∫Lr,φf(x,y)ds---(1)]]>

    用反投影成像方法重建图像g(x，y)，则：

     g(x,y)=Σi=1mp(xcosφi+ysinφi,φi)Δφ---(2)]]>

    其中，φ_i为第i次投影所对应的角度，Δφ为投影角度采样间隔，m为总的投影个数。

    反射层析(基于距离)类似于上述CT反投影重建图像的操作过程，其实质是将到源点某一距离的有限表面区域反射所得的投影信息，反投影到离源点同此距离的所有区域，如图2反射层析激光雷达所示。所不同的是CT依据的透射特征系数，只需要0-π的透射投影信息，而反射层析依据的是反射特征系数，需要0-2π的反射投影信息。

    随着激光的出现和发展，Charles L.Matson(参见文献R.M.Marino，R.N.Capes，W.E.Keicher et al..Tomographic image reconstruction from laser radarreflective projections[J]，SPIE Laser Radar III，1988，Vol.999：248-263)，Parker.J.K.(参见文献Jeffrey K.Parker，E.B.Cralg，D.I.Klick et al..Reflective tomography：image from range resolved laser radarmearsurements[J]，Applied Optics，1988，27(13)：2642-2643)，F.K.Knight(参见文献Frederick L.Knight，David Klick，Danette P.Ryan-Howard et al..Laserradar reflective tomography utilizing a steak camera for precise rangeresolution[J].Applied Optics，1989，28(12)：2196-2198)等很快的将反射层析算法引入到激光雷达中。在物体发生平动情况下，激光短脉冲探测和非相干接收器得到基于距离的多角度一维反射投影信息，所获得的旋转中心投影不再对应距离上同一位置。投影图对准就是将不同角度下的旋转中心投影排列在同一直线上，现有的投影图对准采用位相恢复迭代算法(参见文献Stephen D.Ford and CharlesL.Matson，Projection registration in reflective tomography[J]，SPIE，1999，Vol.3815：189-197)，存在计算量大，操作复杂的问题。

    本方法利用特征点投影信息可以迅速的实现含有特征点旋转中心投影的对准，可以避免复杂的数学迭代过程。

    【发明内容】

    本发明的目的在于为反射投影成像系统提供一种反射投影成像投影图中心的对准方法，该方法的原理可靠，易于实现。

    本发明的技术解决方案如下：

    一种反射投影成像投影图中心对准方法，其特点在于包括以下步骤：

    ①激光器发出光脉冲，用会聚透镜使光脉冲会聚后再次发散，通过调节透镜的焦距来改变光斑的大小，使到达探测目标的光斑外轮廓涵盖整个探测目标；

    ②将初始待测目标所对应的角度定义为，其中是在激光器，探测器和目标所确定平面内，激光器初始脉冲方向与水平方向夹角，如图3所示。激光器发出单脉冲照射到待测目标上，探测器记录在角度下所对应的经过目标反射加宽后的反射投影信息

    ③待测目标旋转角度后，探测角度变为激光器发出单个脉冲，探测器记录相应的反射投影信息待测目标旋转角度后，探测角度变为激光器发出单个脉冲，探测器记录相应的反射投影信息直至获得全部投影

    ④观察反射投影，提取反射投影中的特征点。特征点可以选取反射投影中反射强度明显高于周围区域的点(对应于反射投影中的峰值)，也可以是反射强度明显低于周围区域的点(对应于反射投影中的谷值)，如图4所示。值得注意的是要将特征点与随机噪声点分开，可以采用单个角度多个脉冲获得多个投影，然后对投影取平均去除随机噪声。将所提取的特征点定义为特征点A；

    ⑤观察特征点A的反射投影信息，将相邻角度投影Δt变化很小的角度范围定义为局域角度，其中Δt为相邻角度特征点A反射投影信息峰值上升沿(或谷值下降沿)时间延迟。如图5所示，Δt₁为角度所对应的特征点A反射投影信息上升沿时间延迟，Δt₂为角度所对应的特征点A反射投影信息上升沿时间延迟；

    ⑥局域角度确定特征点A到投影图中心的方法如图6所示：

    

    

     ΔL1=cΔt12---(3)]]>

     ΔL2=cΔt22---(4)]]>

    

    β₁＝90°-β₂                (6)

    

    α₃＝90°-α₁                (8)

    

    其中，ΔL₁，ΔL₂为时间延迟Δt₁，Δt₂所对应的距离延迟，为局域角度变量，R为特征点到旋转中心的距离。

    由(5)(6)(7)(8)(9)可得到：

    

    

    将(10)式代入(2)式，

    

    比较(1)(12)式，可以得到：

    

    将(13)式展开可以得到：

    

    

    

    将(14)式两边同时除以可以得到：

    

    将求出的α₂值代入(1)，从而可以求出特征点A到旋转中心距离R值和角度对应的投影角度θ_th：

    

    

    其中，c为光速，Δt₁为角度所对应的特征点A反射投影信息上升沿时间延迟，Δt₂为角度所对应的特征点A反射投影信息上升沿时间延迟；

    ⑦根据步骤④得到半径R和对应投影中心反投影角度θ_th，对所得到的含有特征点A反射信息的反射投影角度进行调整，具体调整公式为：

    

    其中i＝1，2……N，根据式(16)(18)可以得到含有特征点A反射投影信息的角度所对应的角度θ_i，在各个角度含有特征点A反射投影信息的相应投影p_i(t，θ_i)(i＝1，2……N)，特征点A反射投影信息到旋转中心投影距离ΔS_i，如图7所示：

    

    ⑧将单个角度θ_i特征点投影和旋转中心投影所对应的距离信息ΔS_i转换为时间信息Δτ_i：

    

    其中c是光速，由式(20)可以得到特征点投影和旋转中心投影的时间间隔。特征点投影和旋转中心投影的前后顺序，可以通过角度θ_i的正负来确定。θ_i为正，旋转中心投影在特征点投影前Δτ_i处；θ_i为负，旋转中心投影在特征点的后Δτ_i处，如图8所示。

    ⑨通过观测特征点投影和角度θ_i对应的时间延迟Δτ_i将多角度旋转中心投影排列在同一直线上，实现投影图的对准，如图9所示。

    本发明的技术效果：

    本发明通过观测特征点的投影信息变化，可以实现多角度下反投影中心的对准，从而解决了多角度反投影中心偏移所带来的重建图像错位问题。

    【附图说明】

    图1透射层析示意图；

    图2反射层析示意图；

    图3夹角φ₁与角度间隔Δφ示意图；

    图4特征点的提取(峰值或谷值)；

    图5不同角度下一维反射投影特征点反射投影信息时间延迟示意图；

    图6特征点确定反投影半径和角度示意图；

    图7特征点投影到旋转中心投影距离示意图；

    图8特征点投影和时间间隔Δτ_i确定单角度θ_i投影中心；

    图9多角度投影图对准示意图。

    【具体实施方式】

    下面结合附图和实施例进一步详细说明本发明，但不应以此限制本发明的保护范围。

    一种反射投影成像投影图中心对准方法，包括以下步骤：

    ①激光器发出光脉冲，用会聚透镜使光脉冲会聚后再次发散，通过调节透镜的焦距来改变光斑的大小，使到达探测目标的光斑外轮廓涵盖整个探测目标；

    ②将初始待测目标所对应的角度定义为其中是在激光器，探测器和目标所确定平面内，激光器初始脉冲方向与水平方向夹角，如图3所示。激光器发出单脉冲照射到待测目标上，探测器记录在角度下所对应的经过目标反射加宽后的反射投影信息

    ③待测目标旋转角度后，探测角度变为激光器发出单个脉冲，探测器记录相应的反射投影信息待测目标旋转角度后，探测角度变为激光器发出单个脉冲，探测器记录相应的反射投影信息直至获得全部投影

    ④观察反射投影，提取反射投影中的特征点。特征点可以选取反射投影中反射强度明显高于周围区域的点(对应于反射投影中的峰值)，也可以是反射强度明显低于周围区域的点(对应于反射投影中的谷值)，如图4所示。值得注意的是要将特征点与随机噪声点分开，可以采用的单个角度多个脉冲获得多个投影，然后对投影取平均去除随机噪声。将所提取的特征点定义为特征点A；

    ⑤观察特征点A的反射投影信息，将相邻角度投影Δt变化很小的角度范围定义为局域角度，其中Δt为相邻角度特征点A反射投影信息峰值上升沿(或谷值下降沿)时间延迟。如图5所示，Δt₁为角度所对应的特征点A反射投影信息上升沿时间延迟，Δt₂为角度所对应的特征点A反射投影信息上升沿时间延迟；

    ⑥局域角度确定特征点A到投影图中心的方法如图6所示：

    

    

     ΔL1=cΔt12---(3)]]>

     ΔL2=cΔt22---(4)]]>

    

    β₁＝90°-β₂          (6)

    

    α₃＝90°-α₁           (8)

    

    其中，ΔL₁，ΔL₂为时间延迟Δt₁，Δt₂所对应的距离延迟，为局域角度变量，R为特征点到旋转中心的距离。

    由(5)(6)(7)(8)(9)可得到：

    

    

    将(10)式代入(2)式，

    

    比较(1)(12)式，可以得到：

    

    将(13)式展开可以得到：

    

    

    

    将(14)式两边同时除以可以得到：

    

    将求出的α₂值代入(1)，从而可以求出特征点A到旋转中心距离R值和角度对应的投影角度θ_th：

    

    

    其中，c为光速，Δt₁为角度所对应的特征点A反射投影信息上升沿时间延迟，Δt₂为角度所对应的特征点A反射投影信息上升沿时间延迟；

    ⑦根据步骤④得到半径R和对应投影中心反投影角度θ_th，对所得到的含有特征点A反射信息的反射投影角度进行调整，具体调整公式为：

    

    其中i＝1，2……N，根据式(16)(18)可以得到含有特征点A反射投影信息的角度所对应的角度θ_i，在各个角度含有特征点A反射投影信息的相应投影p_i(t，θ_i)(i＝1，2……N)，特征点A反射投影信息到旋转中心投影距离ΔS_i，如图7所示：

    

    ⑧将单个角度θ_i特征点投影和旋转中心投影所对应的距离信息ΔS_i转换为时间信息Δτ_i：

    

    其中c是光速，由式(20)可以得到特征点投影和旋转中心投影的时间间隔。特征点投影和旋转中心投影的前后顺序，可以通过角度θ_i的正负来确定。θ_i为正，旋转中心投影在特征点投影前Δτ_i处；θ_i为负，旋转中心投影在特征点的后Δτ_i处，如图8所示。

    ⑨通过观测特征点投影和角度θ_i对应的时间延迟Δτ_i将多角度旋转中心投影排列在同一直线上，实现投影图的对准，如图9所示。