CN107091810B - 基于线阵探测器的旋转式光学层析成像系统的成像方法 - Google Patents

基于线阵探测器的旋转式光学层析成像系统的成像方法 Download PDF

Info

Publication number
CN107091810B
CN107091810B CN201710167611.9A CN201710167611A CN107091810B CN 107091810 B CN107091810 B CN 107091810B CN 201710167611 A CN201710167611 A CN 201710167611A CN 107091810 B CN107091810 B CN 107091810B
Authority
CN
China
Prior art keywords
linear array
array detector
image
reflecting optics
optics group
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
CN201710167611.9A
Other languages
English (en)
Other versions
CN107091810A (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.)
Hefei Polytechnic University
Original Assignee
Hefei Polytechnic 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 Hefei Polytechnic University filed Critical Hefei Polytechnic University
Priority to CN201710167611.9A priority Critical patent/CN107091810B/zh
Publication of CN107091810A publication Critical patent/CN107091810A/zh
Application granted granted Critical
Publication of CN107091810B publication Critical patent/CN107091810B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/01—Arrangements or apparatus for facilitating the optical investigation
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/01—Arrangements or apparatus for facilitating the optical investigation
    • G01N2021/0106—General arrangement of respective parts
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N2021/178—Methods for obtaining spatial resolution of the property being measured
    • G01N2021/1785—Three dimensional
    • G01N2021/1787—Tomographic, i.e. computerised reconstruction from projective measurements
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2201/00—Features of devices classified in G01N21/00
    • G01N2201/06—Illumination; Optics
    • G01N2201/068—Optics, miscellaneous

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明属于光学成像领域，具体地讲涉及一种基于线阵探测器的旋转式光学层析成像系统，该系统包括沿着待测物待测面的正面方向依次顺序布置的反射镜片组、柱面镜和线阵探测器；所述线阵探测器连接计算机。本发明还相应地提供了一种利用上述系统进行成像的方法。本发明利用了基于线阵探测器的旋转式光学层析技术，同时采用旋转式光学系统，具有单一的旋转装置，结构比较简单，克服了圆锥扫描时的旋转偏心作用和光学系统离焦作用，减小了后期算法重构时的算法修正难度，提升了系统的稳定性。

Description

基于线阵探测器的旋转式光学层析成像系统的成像方法
技术领域
本发明属于光学成像领域，具体地讲涉及一种基于线阵探测器的旋转式光学层析成像系统的成像方法。
背景技术
随着光学技术的发展以及人们对图像质量的要求越来越高，各种新的成像系统也在不断出现。但是目前现有的成像方法主要是利用圆锥扫描式光学层析系统进行成像，圆锥扫描式光学层析系统主要是将圆锥扫描后的图像经过调制盘进行分割从而获得层析信号，该系统结构比较复杂，所需的光学元器件较多，加上受到圆锥扫描时的旋转偏心和光学系统离焦的作用，使后期图像重构时的算法和修正难度较大，系统的稳定性较差。
由于存在上述缺陷，现有的利用圆锥扫描式光学层析系统的成像系统显然已经不能满足人们对成像的需求。因此，提出一种结构简单、系统稳定的光学成像系统很有必要。
发明内容
根据现有技术中存在的问题，本发明的目的是提供了一种基于线阵探测器的旋转式光学层析成像系统的成像方法，本系统采用旋转式光学系统，具有单一的旋转装置，结构比较简单，克服了圆锥扫描时的旋转偏心作用和光学系统离焦作用，减小了后期算法重构时的算法修正难度，提升了系统的稳定性。
本发明采用以下技术方案：
一种基于线阵探测器的旋转式光学层析成像系统的成像方法，其特征在于，基于线阵探测器的旋转式光学层析成像系统包括沿着待测物待测面的正面方向依次顺序布置的反射镜片组、柱面镜和线阵探测器；所述线阵探测器连接计算机；所述反射镜片组的中心轴为射入和射出反射镜片组的待测物待测面反射光光场的中心线，反射镜片组可绕其中心轴旋转；所述待测物、柱面镜、线阵探测器的中心点与反射镜片组的旋转中心轴均处于一条直线上；所述反射镜片组由电机带动旋转，线阵探测器固定在平移台上，电机与平移台均与计算机电连接；所述电机采用步进的控制方式，电机每转动设定的周期角度带动反射镜片组转动一个角度，所述线阵探测器获取该角度下待测物的一维投影；所述电机和所述线阵探测器通过信号同步控制器实现同步；所述反射镜片组由镜片一、镜片二、镜片三组成；所述镜片三平行于所述中心轴，所述镜片一与镜片二关于镜片三的与中心轴垂直的中垂面对称布置，所述镜片一与镜片三的夹角为0°～45°；射入反射镜片组的待测物待测面反射光依次通过镜片一、镜片三、镜片二的反射射出反射镜片组；成像方法包括如下步骤：
S1，将标的图像置于待测物的位置，标的图像在光照下产生的反射光射入反射镜片组，经过反射镜片组内镜片的反射后，所射出的光线射在柱面镜的凸面上，经过柱面镜的汇聚作用，光线在线阵探测器的靶面上形成一条一维投影，线阵探测器采集一维投影信息，获得标的图像的一维投影函数；
S2，电机带动反射镜片组绕中心轴旋转，使标的图像在光照下产生的另一个角度的反射光射入反射镜片组，经过反射镜片组内镜片的反射后，所射出的光线射在柱面镜的凸面上，经过柱面镜的汇聚作用，光线在线阵探测器的靶面上形成另一条一维投影；线阵探测器采集一维投影信息，获得标的图像的一维投影函数；然后继续转动反射镜片组，直至获得标的图像在反射镜片组处于不同角度时的多个一维投影函数，即标的图像的多个采样数据；
S3，利用采样拟合圆曲线方法，对步骤S2中获得的标定图像的多个采样数据进行采样数据的拟合，得到拟合中心，所述拟合中心即反射镜片组的旋转中心；计算机控制平移台带动线阵探测器在其靶面所在的二维平面内移动，使线阵探测器的靶面中心与所述拟合中心重合，此时线阵探测器的靶面中心和反射镜片组的旋转中心初步校准；
S4，将步骤S1和S2中的标的图像换成待测物，对待测物再次进行步骤S1和S2，获得待测物待测面在反射镜片组处于不同角度时的多个一维投影，线阵探测器采集一维投影信息，获得多个一维投影函数，即待测物待测面的多个采样数据；
S5，对步骤S4中获得的待测物待测面的多个采样数据进行一维滤波处理，得到多个经过修正的线投影，再将此修正后的线投影做反投影运算，即进行反投影重建，得到多个投影图像，将此多个投影图像整合重建成一个初步重构图像；再运用投影正弦图法对所述初步重构图像确定投影中心；
S6，对步骤S5中获得的初步重构图像计算梯度值，并根据梯度优化判定法，选择步骤S5中得到的投影中心偏移若干像素距离的像素坐标点为多个新的投影中心，根据多个新的所述投影中心进行图像重构并计算各个新的重构图像梯度值，选择初步重构图像和各个新的重构图像中梯度值最大的对应的重构图像作为最终的重构图像。
优选的，所述步骤S6的具体过程如下：计算步骤S5中得到的初步重构图像的梯度值为A，记录存储下梯度值；对投影中心像素点向重构图像的左或右移动一个像素点；
若第一次向左移动一个像素点的距离作为新的投影中心，以新的投影中心重构图像后，计算其梯度值为B；若梯度值B比A值大，则继续重复刚才的步骤，即不断对中心像素点位置向同一方向移动，每移动一个像素点位置作为新的投影中心，以新的投影中心重构图像，计算其梯度值并与前一个重构图像的梯度值比较，直至重构图像的梯度值不再增大，此时最大梯度值对应的中心像素点位置确定为重构质量最好的投影中心，相应的重构图像为最佳重构图像；
若梯度值B比A值小，则对投影中心像素点位置向相反的右方向移动一个像素点作为新的投影中心，此时再以新的投影中心进行图像重构，并计算出重构后的图像梯度值C，并与梯度值A比较，若梯度值C比A值大，则也继续重复刚才的步骤，即对投影中心像素点位置继续向右方向移动，每移动一个像素点位置作为新的投影中心，重构图像后计算其梯度值并与前一个梯度值比较，直至重构图像的梯度值不再增大，此时最大梯度值对应的中心像素点位置确定为重构质量最好的投影中心，相应的重构图像为最佳重构图像；
若梯度值C比A值小，则可以确定此时梯度值A、B、C中A最大，则梯度值A对应的中心像素点位置就是重构图像质量最好的投影中心，相应的重构图像即为最佳重构图像；
优选的，所述标的图像为黑色背景，背景内有一个便于识别的白点。
优选的，在整个成像过程中，所述反射镜片组共旋转360度，反射镜片组每次旋转的角度均相同，且反射镜片组绕其中心轴旋转一周所需的次数为不小于75次的奇数次。
本发明的优点和有益效果在于：
1)本发明由待测物、反射镜片组、柱面镜、线阵探测器和计算机组成，本系统采用基于线阵探测器的旋转式光学层析成像技术得到待测物待测面的多个一维投影，并对多个一维投影进行图像重构，实现对待测物待测面的成像。由于本系统采用旋转式光学系统，具有单一的旋转装置，结构比较简单，克服了圆锥扫描时的旋转偏心作用和光学系统离焦作用，减小了后期算法重构时的算法修正难度，提升了系统的稳定性。
2)本发明旋转式光学层析系统获得层析信号是通过柱面镜的光学积分配合线阵探测器并行接收完成的，相比传统沿用的层析扫描系统通过调制盘分割并用单像素光电探测器接收的方式，具有高速、抗干扰的特性。
3)本发明引用采样拟合圆曲线方法对反射镜片组的旋转中心与线阵探测器的靶面中心进行标定和校准，利用滤波反投影算法对投影进行初步图像重构，运用投影正弦图法对初步重构图像确定投影中心，最后运用梯度优化判定法进行投影中心的最终校准和对重构图像最终校正。通过上述投影中心的再校准和图像的再重构，能够得到高质量的重构图像。
4)本发明采用标的图像对反射镜片组的旋转中心与线阵探测器的靶面中心进行校准，且标的图像为黑色背景、内有一个便于识别的白点，由于黑色不反射可见光，而白色可以反射可见光，因此利用仅有的一个白点的反射光可以简单、便捷的对投影中心进行标定，对反射镜片组的旋转中心与线阵探测器的靶面中心进行校准。
5)本发明反射镜片组绕其中心轴旋转一周所需的次数为奇数次，反射镜片组在旋转一周过程中避免了同一个像素点相对投影中心对称的位置上产生相同的扫描结果，也就避免了采样冗余，使成像效果达到最佳。
附图说明
图1为本发明结构示意图。
图2为本发明标的图像示意图。
图3a～3i为本发明将金山WPS图标作为待测物的多个重构图像。
附图标记：
10-待测物，20-反射镜片组，21-镜片一，22-镜片二，
23-镜片三，30-柱面镜，40-线阵探测器，50-计算机，60-平移台。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
如图1所示，本发明中的基于线阵探测器的旋转式光学层析成像系统，包括沿着待测物10待测面的正面方向依次顺序布置的反射镜片组20、柱面镜30和线阵探测器40；所述线阵探测器40连接计算机50；所述反射镜片组20的中心轴为射入和射出反射镜片组20反射光光场的中心线，反射镜片组20可绕其中心轴旋转；所述待测物10、柱面镜30、线阵探测器40的中心点与反射镜片组20的旋转中心轴均处于一条直线上；所述反射镜片组20由电机带动旋转，线阵探测器40固定在平移台60上，电机与平移台60均与计算机50电连接。
所述电机采用步进的控制方式，电机每转动设定的周期角度带动反射镜片组20转动一个角度，所述线阵探测器40便可获取该角度下待测物10的一维投影；所述电机和所述线阵探测器40通过信号同步控制器实现同步。
由于所述反射镜片组20由镜片一21、镜片二22、镜片三23组成，且所述镜片三23平行于所述中心轴，所述镜片一21与镜片二22关于镜片三23的与中心轴垂直的中垂面对称布置，所述镜片一21与镜片三23的夹角为0°～45°，射入反射镜片组20的待测物10待测面反射光依次通过镜片一21、镜片三23、镜片二22的反射射出反射镜片组20。因此，所述反射镜片组20就具备了如下特性：反射镜片组20旋转一个θ角度，所成的像就会旋转2θ角度。
反射镜片组20单次旋转角度不同，获得的一维投影个数不同，成像效果也就存在差异，因此旋转角度可根据实际情况进行选择。但是必须指出的是，这里设定反射镜片组20旋转一周所需的旋转次数应为奇数，由于反射镜片组20旋转一个θ角度，所成的像就会旋转2θ角度，若为偶数，反射镜片组20则在旋转一周过程中同一个像素点相对投影中心对称的位置上产生相同的投影结果，这样就导致了采样冗余，影响成像效果。因此，反射镜片组20旋转一周所需的旋转次数应为奇数。
本实施例中，设定反射镜片组20旋转一周所需的次数为225次，即每次旋转角度为1.6°，同时设步进电机与反射镜片组20传动比为1:30，即电机旋转30圈反射镜片组20转动一圈，反射镜片组20转动一周获得225个不同角度下的投影。旋转台每转动1.6°，等待3s供线阵探测器40进行采集投影信息，然后再转动，不断重复上述动作，转完一周为止。
所述柱面镜30采用平凸柱面透镜，柱面镜30的凸面正对反射镜片组20，柱面镜30的平面垂直于所述中心轴；让反射光从柱面镜30的凸面入射，从柱面镜30的平面射出，对通过柱面镜30的光线进行汇聚。
所述柱面镜30的所有焦点均位于线阵探测器40的靶面上，使光线通过柱面镜30的汇聚作用后全部落在线阵探测器40的靶面上形成一维投影直线，以便线阵探测器40进行采集数据。需要指出的是，由于初始设置时，所述柱面镜30的所有焦点均位于线阵探测器40的靶面上，因此后续通过调整平移台60使反射镜片组20的旋转中心与线阵探测器40的靶面中心重合时，所述平移台60只能在线阵探测器40靶面所在的二维平面内移动，这样也能同时保持线阵探测器40的靶面与柱面镜30的所有焦点始终处于重合的状态。
所述待测物10置于有自然光照射的环境下，以便产生自然反射光线进行后续成像工作。同样地，用于中心校准的标的图像也置于有自然光照射的环境下，该标的图像为黑色背景，背景内有一个便于识别的白点。
本发明实施例对各器件的型号、样式除做特殊说明的以外，其他器件的型号、样式不做限制，只要能完成上述功能的器件均可。
下面结合具体工作过程和附图，对本发明的基于线阵探测器的旋转式光学层析成像系统得成像方法进行详细说明：
S1，将标的图像置于待测物10的位置，标的图像在自然光的照射下产生的反射光射入反射镜片组20，经过反射镜片组20内镜片的反射后，所射出的光线射在柱面镜30的凸面上，经过柱面镜30的汇聚作用，光线在线阵探测器40的靶面上形成一条一维投影，线阵探测器40采集一维投影信息，获得标的图像的一维投影函数；
S2，电机带动反射镜片组20绕中心轴旋转，使标的图像在自然光的照射下产生的另一个角度的反射光射入反射镜片组20，经过反射镜片组20内镜片的反射后，所射出的光线射在柱面镜30的凸面上，经过柱面镜30的汇聚作用，光线在线阵探测器40的靶面上形成另一条一维投影；线阵探测器40采集一维投影信息，获得标的图像的一维投影函数；然后继续转动反射镜片组20，直至获得标的图像在反射镜片组20处于不同角度时的多个一维投影函数，即标的图像的多个采样数据；
S3，利用采样拟合圆曲线方法，对步骤S2中获得的标定图像的多个采样数据进行采样数据的拟合，得到拟合中心，所述拟合中心即反射镜片组20的旋转中心；计算机50控制平移台60带动线阵探测器40在其靶面所在的二维平面内移动，使线阵探测器40的靶面中心与所述拟合中心重合，此时线阵探测器40的靶面中心和反射镜片组20的旋转中心初步校准；
S4，将步骤S1和S2中的标的图像换成金山WPS图标，对金山WPS图标再次进行步骤S1和S2，获得金山WPS图标待测面在反射镜片组20处于不同角度时的多个一维投影，线阵探测器40采集一维投影信息，获得多个一维投影函数，即金山WPS图标待测面的多个采样数据；
S5，对步骤S4中获得的金山WPS图标待测面的多个采样数据进行一维滤波处理，得到多个经过修正的线投影，再将此修正后的线投影做反投影运算，即进行反投影重建，得到多个投影图像，将此多个投影图像整合重建成一个初步重构图像；再运用投影正弦图法对所述初步重构图像确定投影中心；
S6，对步骤S5中获得的初步重构图像计算梯度值，并根据梯度优化判定法，选择步骤S5中得到的投影中心偏移若干像素距离的像素坐标点为多个新的投影中心，根据所述多个新的投影中心进行图像重构并计算各个新的重构图像梯度值，选择初步重构图像和各个新的重构图像中梯度值最大的对应的重构图像作为最终的重构图像。具体过程如下：
计算步骤S5中得到的初步重构图像的梯度值为A，记录存储下梯度值；由于柱面镜30的汇聚作用，光线投影在纵方向上得到压缩，在线阵探测器40的靶面上形成一条一维投影，因此重构图像的质量主要取决于投影中心的横向左右偏移。因此，这里对投影中心像素点向左或右移动一个像素点；若第一次向左移动一个像素点的距离，以新的投影中心重构图像后，计算其梯度值为B；由于梯度峰值具有唯一性且呈现递减趋势，因此，若梯度值B比A值大，则继续重复刚才的步骤，即不断对中心像素点位置向同一方向移动，每移动一个像素点位置，以新的投影中心重构图像，计算其梯度值并与前一个重构图像的梯度值比较，直至重构图像的梯度值不再增大，此时最大梯度值对应的中心像素点位置确定为重构质量最好的投影中心，相应的重构图像为最佳重构图像；
若梯度值B比A值小，则对中心像素点位置向相反的方向移动一个像素点，此时再以新的投影中心进行图像重构，并计算出重构后的图像梯度值C，并与梯度值A比较，若梯度值C比A值大，则也继续重复刚才的步骤，即继续对中心像素点位置向同一方向移动，每移动一个像素点位置，重构图像后计算其梯度值并与前一个梯度值比较，直至重构图像的梯度值不再增大，此时最大梯度值对应的中心像素点位置确定为重构质量最好的投影中心，相应的重构图像为最佳重构图像；
若梯度值C比A值小，则可以确定此时梯度值A、B、C中A最大，则梯度值A对应的中心像素点位置就是重构图像质量最好的投影中心，相应的重构图像即为最佳重构图像。
如图3a～图3i所示，其中图3e为步骤S5中获得的金山WPS图标的初步重构图像，图3a～图3d和图3f～图3i为以周边中心像素点为投影中心重构出的重构图像，计算图3a～图3i对应重构图像的梯度值，如表1所示：
表1：
图像 图3a 图3b 图3c 图3d 图3e 图3f 图3g 图3h 图3i
梯度值 2619.9 2630.2 2721.7 2731.5 3219.9 2916.6 2910.0 2872.1 2862.2
通过表1可得知，步骤S5中获得的金山WPS图标的初步重构图像图3e对应的梯度值最大，则所述初步重构图像图3e对应的重构中心点为最终的投影中心，所述初步重构图像图3e为最终的重构图像。
综上所述，本发明针对物体成像难题，引入基于线阵探测器的旋转式光学层析的方法，通过结构简单的反射镜片组20的旋转获得不同角度下的投影，利用线阵探测器40采集投影信息，并利用计算机50进行采样拟合圆曲线法对反射镜片组20的旋转中心和线阵探测器40的投影中心进行初步校准，利用滤波反投影算法对投影进行初步图像重构，运用投影正弦图法对初步重构图像确定投影中心，最后运用梯度优化判定法进行投影中心的最终校准和对重构图像最终校正，实现了物体的快速稳定成像，满足了实际应用中的多种需要。

Claims (4)

1.一种基于线阵探测器的旋转式光学层析成像系统的成像方法，其特征在于，基于线阵探测器的旋转式光学层析成像系统包括沿着待测物(10)待测面的正面方向依次顺序布置的反射镜片组(20)、柱面镜(30)和线阵探测器(40)；所述线阵探测器(40)连接计算机(50)；所述反射镜片组(20)的中心轴为射入和射出反射镜片组(20)的待测物(10)待测面反射光光场的中心线，反射镜片组(20)可绕其中心轴旋转；所述待测物(10)、柱面镜(30)、线阵探测器(40)的中心点与反射镜片组(20)的旋转中心轴均处于一条直线上；所述反射镜片组(20)由电机带动旋转，线阵探测器(40)固定在平移台(60)上，电机与平移台(60)均与计算机(50)电连接；所述电机采用步进的控制方式，电机每转动设定的周期角度带动反射镜片组(20)转动一个角度，所述线阵探测器(40)获取该角度下待测物(10)的一维投影；所述电机和所述线阵探测器(40)通过信号同步控制器实现同步；所述反射镜片组(20)由镜片一(21)、镜片二(22)、镜片三(23)组成；所述镜片三(23)平行于所述中心轴，所述镜片一(21)与镜片二(22)关于镜片三(23)的与中心轴垂直的中垂面对称布置，所述镜片一(21)与镜片三(23)的夹角为0°～45°；射入反射镜片组(20)的待测物(10)待测面反射光依次通过镜片一(21)、镜片三(23)、镜片二(22)的反射射出反射镜片组(20)；包括以下步骤：
S1，将标的图像置于待测物(10)的位置，标的图像在光照下产生的反射光射入反射镜片组(20)，经过反射镜片组(20)内镜片的反射后，所射出的光线射在柱面镜(30)的凸面上，经过柱面镜(30)的汇聚作用，光线在线阵探测器(40)的靶面上形成一条一维投影，线阵探测器(40)采集一维投影信息，获得标的图像的一维投影函数；
S2，电机带动反射镜片组(20)绕中心轴旋转，使标的图像在光照下产生的另一个角度的反射光射入反射镜片组(20)，经过反射镜片组(20)内镜片的反射后，所射出的光线射在柱面镜(30)的凸面上，经过柱面镜(30)的汇聚作用，光线在线阵探测器(40)的靶面上形成另一条一维投影；线阵探测器(40)采集一维投影信息，获得标的图像的一维投影函数；然后继续转动反射镜片组(20)，直至获得标的图像在反射镜片组(20)处于不同角度时的多个一维投影函数，即标的图像的多个采样数据；
S3，利用采样拟合圆曲线方法，对步骤S2中获得的标定图像的多个采样数据进行采样数据的拟合，得到拟合中心，所述拟合中心即反射镜片组(20)的旋转中心；计算机(50)控制平移台(60)带动线阵探测器(40)在其靶面所在的二维平面内移动，使线阵探测器(40)的靶面中心与所述拟合中心重合，此时线阵探测器(40)的靶面中心和反射镜片组(20)的旋转中心初步校准；
S4，将步骤S1和S2中的标的图像换成待测物(10)，对待测物(10)再次进行步骤S1和S2，获得待测物(10)待测面在反射镜片组(20)处于不同角度时的多个一维投影，线阵探测器(40)采集一维投影信息，获得多个一维投影函数，即待测物(10)待测面的多个采样数据；
S5，对步骤S4中获得的待测物(10)待测面的多个采样数据进行一维滤波处理，得到多个经过修正的线投影，再将此修正后的线投影做反投影运算，即进行反投影重建，得到多个投影图像，将此多个投影图像整合重建成一个初步重构图像；再运用投影正弦图法对所述初步重构图像确定投影中心；
S6，对步骤S5中获得的初步重构图像计算梯度值，并根据梯度优化判定法，选择步骤S5中得到的投影中心偏移若干像素距离的像素坐标点为多个新的投影中心，根据多个新的所述投影中心进行图像重构并计算各个新的重构图像梯度值，选择初步重构图像和各个新的重构图像中梯度值最大的重构图像作为最终的重构图像。
2.根据权利要求1所述的一种基于线阵探测器的旋转式光学层析成像系统的成像方法，其特征在于，所述步骤S6的具体过程如下：
计算步骤S5中得到的初步重构图像的梯度值为A，记录存储下梯度值；对投影中心像素点向重构图像的左或右移动一个像素点；
若第一次向左移动一个像素点的距离作为新的投影中心，以新的投影中心重构图像后，计算其梯度值为B；若梯度值B比A值大，则继续重复刚才的步骤，即不断对中心像素点位置向同一方向移动，每移动一个像素点位置作为新的投影中心，以新的投影中心重构图像，计算其梯度值并与前一个重构图像的梯度值比较，直至重构图像的梯度值不再增大，此时最大梯度值对应的中心像素点位置确定为重构质量最好的投影中心，相应的重构图像为最佳重构图像；
若梯度值B比A值小，则对投影中心像素点位置向相反的右方向移动一个像素点作为新的投影中心，此时再以新的投影中心进行图像重构，并计算出重构后的图像梯度值C，并与梯度值A比较，若梯度值C比A值大，则也继续重复刚才的步骤，即对投影中心像素点位置继续向右方向移动，每移动一个像素点位置作为新的投影中心，重构图像后计算其梯度值并与前一个梯度值比较，直至重构图像的梯度值不再增大，此时最大梯度值对应的中心像素点位置确定为重构质量最好的投影中心，相应的重构图像为最佳重构图像；
若梯度值C比A值小，则可以确定此时梯度值A、B、C中A最大，则梯度值A对应的中心像素点位置就是重构图像质量最好的投影中心，相应的重构图像即为最佳重构图像。
3.根据权利要求1所述的一种基于线阵探测器的旋转式光学层析成像系统的成像方法，其特征在于：所述标的图像为黑色背景，背景内有一个便于识别的白点。
4.根据权利要求1所述的一种基于线阵探测器的旋转式光学层析成像系统的成像方法，其特征在于：在整个成像过程中，所述反射镜片组(20)共旋转360度，反射镜片组(20)每次旋转的角度均相同，且反射镜片组(20)绕其中心轴旋转一周所需的次数为不小于75次的奇数次。
CN201710167611.9A 2017-03-21 2017-03-21 基于线阵探测器的旋转式光学层析成像系统的成像方法 Active CN107091810B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710167611.9A CN107091810B (zh) 2017-03-21 2017-03-21 基于线阵探测器的旋转式光学层析成像系统的成像方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710167611.9A CN107091810B (zh) 2017-03-21 2017-03-21 基于线阵探测器的旋转式光学层析成像系统的成像方法

Publications (2)

Publication Number Publication Date
CN107091810A CN107091810A (zh) 2017-08-25
CN107091810B true CN107091810B (zh) 2019-10-11

Family

ID=59646227

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710167611.9A Active CN107091810B (zh) 2017-03-21 2017-03-21 基于线阵探测器的旋转式光学层析成像系统的成像方法

Country Status (1)

Country Link
CN (1) CN107091810B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110310338B (zh) * 2019-06-24 2022-09-06 西北工业大学 一种基于多中心投影模型的光场相机标定方法
CN112866529B (zh) * 2021-01-11 2022-04-08 合肥工业大学 单元探测器光学层析扫描时分调制成像系统
CN113126116B (zh) * 2021-04-08 2023-12-26 合肥工业大学 基于时分幅度调制的单像素层析扫描装置

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1766532A (zh) * 2005-12-09 2006-05-03 北京理工大学 空间调制干涉型计算层析成像光谱仪
CN101032390A (zh) * 2007-03-29 2007-09-12 浙江大学 用于在体光学活检的谱域oct内窥成像系统
CN201945702U (zh) * 2011-01-30 2011-08-24 福建中策光电科技有限公司 一种新型道威棱镜
CN102824185A (zh) * 2012-09-12 2012-12-19 北京大学 结合透声反光镜的光声层析成像系统及其成像方法
CN103728724A (zh) * 2013-12-19 2014-04-16 合肥工业大学 一种层析扫描系统及层析扫描方法
CN106112256A (zh) * 2016-07-28 2016-11-16 上海嘉强自动化技术有限公司 一种基于二分束棱镜与扫描振镜组合光学系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9737207B2 (en) * 2013-12-05 2017-08-22 Shenzhen Certainn Technology Co., Ltd. Method for quick switching to realize anterior and posterior eye segments imaging

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1766532A (zh) * 2005-12-09 2006-05-03 北京理工大学 空间调制干涉型计算层析成像光谱仪
CN101032390A (zh) * 2007-03-29 2007-09-12 浙江大学 用于在体光学活检的谱域oct内窥成像系统
CN201945702U (zh) * 2011-01-30 2011-08-24 福建中策光电科技有限公司 一种新型道威棱镜
CN102824185A (zh) * 2012-09-12 2012-12-19 北京大学 结合透声反光镜的光声层析成像系统及其成像方法
CN103728724A (zh) * 2013-12-19 2014-04-16 合肥工业大学 一种层析扫描系统及层析扫描方法
CN106112256A (zh) * 2016-07-28 2016-11-16 上海嘉强自动化技术有限公司 一种基于二分束棱镜与扫描振镜组合光学系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
光学层析成像系统研究;黄俊玲;《中国优秀硕士学位论文全文数据库信息科技辑》;20160515(第5期);第23-29、33页以及图4.2 *

Also Published As

Publication number Publication date
CN107091810A (zh) 2017-08-25

Similar Documents

Publication Publication Date Title
CN107091810B (zh) 基于线阵探测器的旋转式光学层析成像系统的成像方法
US9429665B2 (en) Radiation penetration system and calibration of the same
US7950849B2 (en) Method and device for geometry analysis and calibration of volumetric imaging systems
CN107314742B (zh) 一种旋转式光学层析成像系统及成像方法
CN102711613B (zh) 计算断层摄影成像方法及系统
US6140650A (en) Calibration of PET cameras
US7277522B2 (en) Method for focus adjustment in a CT apparatus
CN110533597B (zh) 伪影处理及旋转中心确定方法、装置与设备、存储介质
JPH0661326B2 (ja) らせん走査で断層像を作成する方法および装置
KR20040111005A (ko) 단층촬영장치
WO2004034908A1 (ja) X線断層撮影装置
CN101398397A (zh) 多次扫描模式的ct成像方法
US20110304856A1 (en) Lightwave interference measurement apparatus
CN110023998A (zh) 补偿运动模糊的从表面或切片扫描对对象的三维模型的重构
CN107631700B (zh) 三维扫描仪与全站仪相结合的三维数据测量方法
CN106821405B (zh) 一种x光机的参数标定方法、装置及系统
US8861829B2 (en) Method and system for reconstruction of tomographic images
CN105844594A (zh) 一种工业ct扫描系统旋转中心校正方法
US7596203B2 (en) Computer tomography method
WO2000062674A1 (en) Half field of view reduced-size ct detector
US9095259B2 (en) Method and system for high resolution nutated slice reconstruction using quarter detector offset
JPH03188833A (ja) 断層撮影像作成方法および装置
US11175242B2 (en) Geometric alignment, sample motion correction, and intensity normalization of computed tomography projections using pi-line optimization
WO2020054094A1 (ja) 三次元レーザー光走査装置の計測データ較正方法
CN1482433A (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