CN103308452B - 一种基于景深融合的光学投影断层成像图像获取方法 - Google Patents
一种基于景深融合的光学投影断层成像图像获取方法 Download PDFInfo
- Publication number
- CN103308452B CN103308452B CN201310199983.1A CN201310199983A CN103308452B CN 103308452 B CN103308452 B CN 103308452B CN 201310199983 A CN201310199983 A CN 201310199983A CN 103308452 B CN103308452 B CN 103308452B
- Authority
- CN
- China
- Prior art keywords
- image
- matrix
- dimensional image
- pixel value
- pixel
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 230000003287 optical effect Effects 0.000 title claims abstract description 28
- 230000004927 fusion Effects 0.000 title claims abstract description 12
- 238000003325 tomography Methods 0.000 title abstract description 7
- 239000011159 matrix material Substances 0.000 claims description 39
- 238000003384 imaging method Methods 0.000 claims description 28
- 230000008569 process Effects 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 2
- 238000006467 substitution reaction Methods 0.000 claims description 2
- 230000002194 synthesizing effect Effects 0.000 abstract description 2
- 238000002591 computed tomography Methods 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 241000382353 Pupa Species 0.000 description 2
- 230000003321 amplification Effects 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Landscapes
- Image Processing (AREA)
- Microscoopes, Condenser (AREA)
Abstract
本发明公开了一种基于景深融合的光学投影断层成像图像获取方法,其包括:步骤1、通过在同一旋转角度下,移动样本与显微装置的相对位置,采集不同聚焦位置上的二维图像序列;步骤2、在所述二维图像序列中,选取每个像素位置在该二维图像序列中最清晰的点进行合成,生成全景深图。本发明可实现光学投影断层成像***采集的不同聚焦位置图像融合,得到的全景深图像用于三维重建,提高了三维图像的分辨率。
Description
技术领域
本发明涉及一种光学投影断层成像技术领域,尤其涉及一种基于景深融合的光学投影断层成像图像获取方法。
背景技术
光学投影断层成像(Optical Projection Tomography,OPT)技术可以实现1-10毫米尺度生物样本的结构和分子特异性成像。OPT成像过程中样本比较薄且相对透明,可见光穿过样本时散射效应可以忽略不计,样本对光子主要表现为吸收作用,因此可认为光线沿直线传播穿过样本。OPT扫描时将样本固定在一个高精度旋转台上,采用高稳定的LED或激光器作为照射源,照射在样本上,穿过样本的光信号通过显微装置放大后传输到高灵敏度的电子倍增式CCD探测器,由CCD探测器采集并成像,***每采集一幅二维投影图,旋转台旋转一个固定角度,***采集下一幅图像,依次共采集360°范围内的投影数据。因此,光学投影断层成像***扫描采集到的数据是一系列不同角度下光线穿过样本的二维投影图像。将所有投影图像运用滤波反投影方法进行重建,即可得到三维图像。二维投影图的清晰与否,直接关系到三维图像分辨率的高低。
光学投影断层成像技术属于显微成像。在显微光学成像中,低倍显微镜工作距离长,景深比较大,由于显微镜物镜焦深范围小,随着放大倍数的增大,景深会相应减小。只有那些在聚焦平面或其附近的结构才是可见的,这使得即便是结构最简单的、三维深度相对平坦的物体也不可能在一幅图像中完全聚焦清晰。而在光学投影断层成像中,样本厚度一般为2~3mm,当高倍数放大后景深问题变得尤为严重。而观察样本时,要求显微镜成像既要有更高的分辨率又要有足够的景深,这是传统光学硬件的矛盾。然而对于光学投影断层成像而言,这又是不得不解决的问题,所以景深问题是制约光学投影断层成像获得高分辨率图像的一个瓶颈。
发明内容
为解决上述问题,本发明提供了一种基于景深融合的光学投影断层成像图像获取方法。该方法包括:
步骤1、通过在同一旋转角度下,移动样本与显微装置的相对位置,采集不同聚焦位置上的二维图像序列;
步骤2、在所述二维图像序列中,选取每个像素位置在该二维图像序列中最清晰的点进行合成,生成全景深图。
本发明可实现光学投影断层成像***采集的不同聚焦位置图像融合,得到的全景深图像用于三维重建,提高了三维图像的分辨率。
附图说明
图1为本发明实施例中基于景深融合的光学投影断层成像图像获取方法中成像***数据采集的原理示意图;
图2为本发明实施例中基于景深融合的光学投影断层成像图像获取方法中图像合成的步骤流程图;
图3为本发明实施例中基于景深融合的光学投影断层成像图像获取方法的具体实验结果图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明进一步详细说明。虽然本文可提供包含特定值的参数的示范,但应了解,参数无需确切等于相应的值,而是在可接受的误差容限或设计约束内近似于所述值。
本发明公开了一种基于景深融合的光学投影断层成像图像获取方法。本发明专门针对光学投影断层成像在某一固定旋转角度下、不同聚焦位置图像进行景深融合,选取每个像素位置在图像序列上最清晰点的像素值作为合成图像上该像素的灰度值进行图像融合,生成全景深图,从而提高三维图像的分辨率。
本发明公开的上述方法分为两个主要步骤:数据获取、图像融合。其中,数据获取步骤用于利用光学投影断层成像采集某一固定旋转角度下、不同聚焦位置的二维序列图像;图像融合步骤用于根据所获取到的二维图像序列中每个像素位置在该二维图像序列上最清晰的点进行合成,生成全景深图,用于三维重建,从而提高三维重建体分辨率。
本发明公开的上述方法的详细步骤如下:
步骤S1:通过在同一旋转角度下,移动样本与显微装置的相对位置,采集到不同聚焦位置上的二维图像序列。
光学投影断层成像技术属于显微成像,当样本尺寸非常小时,需要采用高放大倍数的显微装置,进行光信号放大。而高倍数放大时,景深随之变得非常小,这将导致在某一固定角度下不能采集到每个位置都很清晰的图像。但是在实际使用光学投影断层成像***时,要求成像既要有高的分辨率又要有足够的景深,这是传统光学硬件的矛盾。解决此问题的一个有效方法是利用数字图像处理技术,对图像序列进行景深融合。
图1示出了本发明中成像***进行数据采集的原理示意图。如图1所示,显微装置和样本旋转中心在同一轴线上。通过径向移动样本与显微装置的相对位置,并采集数据,使得样本的所有清晰部分都包含在所采集到的图像序列中。水平平移台步距选择通过单方向移动平移台,使得样本在探测器上的图像由模糊到清晰,再到模糊。整个过程中,在第一次模糊后,出现第一个聚焦清晰的位置为该次图像序列采集的起点,逐步等距移动样本,使得图像再次模糊后,前一次聚焦清晰位置为终点,记录下所移动的步数和步距,重复下一个角度的数据采集,直到采集完360°范围内不同旋转角度下的所有数据。其中,在采集图像前需要调整样本位置处于探测器的成像区域中,被采集到的二维图像序列之间Y轴和Z轴位置不发生变化,而它们之间的步距(即X轴方向上的位置)影响最终生成的全景深图像的清晰度;水平平移台的运动通过计算机软件控制,从而可获得步距精确且均匀的二维图像序列。
步骤S2:在所述二维图像序列中,选取每个像素位置在图像序列上最清晰点的像素值作为合成图像上该像素的灰度值进行图像融合,生成全景深图。即针对步骤S1得到的某一角度下的多张不同聚焦位置的二维序列图像,进行景深融合。
图2示出了本发明中对不同聚焦位置的序列图像进行景深融合的步骤流程图。如图2所示,该步骤具体包括:
第一步,根据所采集序列图像的大小,设定目标合成图像像素矩阵M0和用于中间数据处理的矩阵T的大小,其大小等于所采集图像的像素大小,并将M0和T初始化为零矩阵,同时,初始化比较图层序号i=1;
第二步,输入待比较的图层数N,即一共有N幅不同聚焦位置的图像;
第三步,判断待比较图像是否已经完成比较,即i是否等于N,若未完成比较,则i=i+1,并转向第四步,否则转向第六步;
第四步,将待比较图像像素矩阵Mi与目标合成图像像素矩阵M0进行比较,生成中间矩阵T,T中全部元素均为0或1,1代表在相应位置上Mi的元素大于M0中的元素;0代表在相应位置上Mi的元素不大于M0中的元素;其中,待比较图像像素矩阵Mi中的元素为所采集的序列图像中待比较图像的像素值。
第五步,根据生成的中间矩阵T更新目标合成图像像素矩阵M0,即在中间矩阵中元素为1的位置,使用Mi中的元素替代M0相应位置的元素,使用图像Mi中的像素值替代图像M0相应位置的像素值,至此,完成一层图像的融合;转向第三步继续判断,若未完成比较,则重复第四、第五步,否则转向第六步;
第六步,完成图像融合,输出目标合成图像像素矩阵M0;
第七步,某一角度的景深融合结束,M0对应的图像即为该角度下的景深融合图像。
图3示出了本发明中利用上述方法对果蝇蛹活体进行成像实验的结果图。实验中采用活体果蝇蛹进行景深融合处理的光学投影断层成像,共采集360°下的二维图像,每个角度采集6张不同聚焦位置的图像,即水平平移台移动5次。
本实验选择平移台步进0.25毫米,移动5次,共采集6张图像。其结果如图3中1、2、3、4、5、6所示,分别表示不同步距下采集到的二维图像。可见,图像序列中图像由模糊到清晰再到模糊,且每张都有清晰的部分。采集到的图像按顺序依次保存,便于图像融合。其中7所示为该角度下运用景深融合得到的全景深图像,由图中可见,该图像集合了图3中1~6幅图中的清晰位置。这将明显改善单一图像的不足,提高结果图像的清晰度及信息包含量,为进一步的科研分析提供更多的样本信息。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (12)
1.一种基于景深融合的光学投影断层成像图像获取方法,包括:
步骤1、通过在同一旋转角度下,移动样本与显微装置的相对位置,采集不同聚焦位置上的二维图像序列;
步骤2、在所述二维图像序列中,选取每个像素位置在该二维图像序列中最清晰的点进行合成,生成全景深图;
其中,步骤2具体包括:
第一步,根据所采集二维图像序列中图像的大小,设定目标合成图像像素矩阵M0和用于中间数据处理的矩阵T的大小,其大小等于所采集图像的像素大小,并将M0和T初始化为零矩阵,同时,初始化比较图像序号i=1;
第二步,输入待比较的图像数N;
第三步,判断待比较图像是否已经完成比较,即i是否等于N,若未完成比较,则i=i+1,并转向第四步,否则转向第六步;
第四步,将待比较图像像素矩阵Mi与目标合成图像像素矩阵M0进行比较,生成中间矩阵T,T中全部元素均为0或1,1代表在相应位置上Mi的元素大于M0中的元素;0代表在相应位置上Mi的元素不大于M0中的元素;其中,待比较图像像素矩阵Mi中的元素为所述二维图像序列中待比较图像的像素值;
第五步,根据生成的中间矩阵T更新目标合成图像像素矩阵M0,即在中间矩阵中元素为1的位置,使用Mi中的元素替代M0相应位置的元素,使用图像Mi中的像素值替代图像M0相应位置的像素值,至此,完成一层图像的融合;转向第三步继续判断,若未完成比较,则重复第四、第五步,否则转向第六步;
第六步,完成图像融合,输出目标合成图像像素矩阵M0;
第七步,某一角度的景深融合结束,M0对应的图像即为该角度下的景深融合图像。
2.按权利要求1所述的方法,其特征在于,该方法采集360°范围内不同旋转角度下的二维图像序列,其先采集某一角度下的不同聚焦位置的二维图像序列后,再旋转至下一个角度后采集该下一个角度的不同聚焦位置的二维图像序列。
3.按权利要求1所述的方法,其特征在于,所述移动样本与显微装置的相对位置是通过水平平移台在X轴方向上移动样本的位置实现。
4.按权利要求3所述的方法,其特征在于,采集图像前调整样本位置处于探测器的成像区域中。
5.按权利要求1所述的方法,其特征在于,被采集到的二维图像序列之间Y轴和Z轴位置不发生变化,而它们之间的步距影响最终生成的全景深图像的清晰度。
6.按权利要求1所述的方法,其特征在于,水平平移台的运动通过计算机控制,从而可获得步距精确且均匀的二维图像序列。
7.按权利要求1所述的方法,其特征在于,步骤2中不同聚焦位置上的二维图像根据每幅图像的像素值,选取在该位置上最清晰点的像素值作为合成图像的像素值。
8.按权利要求7所述的方法,其特征在于,步骤2中图像融合通过逐幅比较二维图像序列中每幅图像的相应像素值实现。
9.按权利要求8所述的方法,其特征在于,逐幅比较时,先初始化目标合成图像像素矩阵、中间矩阵和待比较二维图像序列序号,利用目标合成图像像素矩阵与待比较图像像素矩阵生成中间矩阵,并通过中间矩阵更新目标合成图像像素矩阵中的元素像素值实现。
10.按权利要求9所述的方法,其特征在于,中间矩阵为0、1矩阵,若待比较图像像素矩阵中的元素像素值大于目标合成图像像素矩阵的相应元素像素值,则中间矩阵的相应元素置为1,否则为0。
11.按权利要求9所述的方法,其特征在于,中间矩阵用于目标合成图像像素矩阵中元素像素值的更新,其中将中间矩阵中元素值为1的元素位置对应的目标合成图像像素的元素像素值替换为待比较图像像素矩阵中相应的元素像素值。
12.按权利要求9所述的方法,其特征在于,比较完最后一幅图像后,得到的目标合成图像像素矩阵对应的图像即为最终得到的全景深图。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310199983.1A CN103308452B (zh) | 2013-05-27 | 2013-05-27 | 一种基于景深融合的光学投影断层成像图像获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310199983.1A CN103308452B (zh) | 2013-05-27 | 2013-05-27 | 一种基于景深融合的光学投影断层成像图像获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103308452A CN103308452A (zh) | 2013-09-18 |
CN103308452B true CN103308452B (zh) | 2015-05-06 |
Family
ID=49133911
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310199983.1A Active CN103308452B (zh) | 2013-05-27 | 2013-05-27 | 一种基于景深融合的光学投影断层成像图像获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103308452B (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104065859B (zh) * | 2014-06-12 | 2017-06-06 | 青岛海信电器股份有限公司 | 一种全景深图像的获取方法及摄像装置 |
CN105451010A (zh) * | 2014-08-18 | 2016-03-30 | 惠州友华微电子科技有限公司 | 一种景深获取装置及其获取的方法 |
CN105118088A (zh) * | 2015-08-06 | 2015-12-02 | 曲阜裕隆生物科技有限公司 | 一种基于病理切片扫描装置的3d成像及融合的方法 |
TWI555378B (zh) * | 2015-10-28 | 2016-10-21 | 輿圖行動股份有限公司 | 一種全景魚眼相機影像校正、合成與景深重建方法與其系統 |
CN105407280B (zh) * | 2015-11-11 | 2019-02-12 | Oppo广东移动通信有限公司 | 全景图像合成方法和*** |
CN106339476B (zh) * | 2016-08-30 | 2019-10-29 | 北京寺库商贸有限公司 | 一种图像处理方法和*** |
JP7039179B2 (ja) * | 2017-04-13 | 2022-03-22 | キヤノン株式会社 | 情報処理装置、情報処理システム、情報処理方法及びプログラム |
CN108495032B (zh) * | 2018-03-26 | 2020-08-04 | Oppo广东移动通信有限公司 | 图像处理方法、装置、存储介质及电子设备 |
CN111381357B (zh) * | 2018-12-29 | 2021-07-20 | 中国科学院深圳先进技术研究院 | 图像三维信息提取方法、对象成像方法、装置及*** |
CN110160956B (zh) * | 2019-05-27 | 2022-04-01 | 广州英特美迪科技有限公司 | 病理切片扫描影像分析***及其扫描方法 |
CN111242880B (zh) * | 2019-12-30 | 2023-05-02 | 广州市明美光电技术有限公司 | 一种用于显微镜的多景深图像叠加方法、设备及介质 |
CN111861915B (zh) * | 2020-07-08 | 2021-08-13 | 北京科技大学 | 一种显微成像场景下失焦扩散效应消除方法及装置 |
WO2023102724A1 (zh) * | 2021-12-07 | 2023-06-15 | 宁德时代新能源科技股份有限公司 | 图像的处理方法和*** |
CN114881907B (zh) * | 2022-06-30 | 2022-09-23 | 江苏集萃苏科思科技有限公司 | 光学显微图像多景深焦点合成方法、***及图像处理方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101615289A (zh) * | 2009-08-05 | 2009-12-30 | 北京优纳科技有限公司 | 切片组织的三维采集及多层图像融合方法 |
CN102727188A (zh) * | 2012-07-26 | 2012-10-17 | 中国科学院自动化研究所 | 一种基于拼合螺旋扫描方式的光学投影断层成像方法 |
CN102973246A (zh) * | 2012-12-03 | 2013-03-20 | 中国科学院自动化研究所 | 自动处理景深的光学投影断层成像方法 |
CN103065359A (zh) * | 2013-01-14 | 2013-04-24 | 厦门大学 | 一种光学成像的三维轮廓重建***及其重建方法 |
CN103110425A (zh) * | 2013-03-11 | 2013-05-22 | 西安电子科技大学 | 一种联合成像的锥束x射线发光断层成像***及其方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8508591B2 (en) * | 2010-02-05 | 2013-08-13 | Applied Vision Corporation | System and method for estimating the height of an object using tomosynthesis-like techniques |
-
2013
- 2013-05-27 CN CN201310199983.1A patent/CN103308452B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101615289A (zh) * | 2009-08-05 | 2009-12-30 | 北京优纳科技有限公司 | 切片组织的三维采集及多层图像融合方法 |
CN102727188A (zh) * | 2012-07-26 | 2012-10-17 | 中国科学院自动化研究所 | 一种基于拼合螺旋扫描方式的光学投影断层成像方法 |
CN102973246A (zh) * | 2012-12-03 | 2013-03-20 | 中国科学院自动化研究所 | 自动处理景深的光学投影断层成像方法 |
CN103065359A (zh) * | 2013-01-14 | 2013-04-24 | 厦门大学 | 一种光学成像的三维轮廓重建***及其重建方法 |
CN103110425A (zh) * | 2013-03-11 | 2013-05-22 | 西安电子科技大学 | 一种联合成像的锥束x射线发光断层成像***及其方法 |
Non-Patent Citations (1)
Title |
---|
旋转扫描序列图像的全景图拼接方法;王会峰;《光学学报》;20090531;第29卷(第5期);第1211-1215页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103308452A (zh) | 2013-09-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103308452B (zh) | 一种基于景深融合的光学投影断层成像图像获取方法 | |
CN102980664B (zh) | 超像素微扫描方法及相应的红外超分辨率实时成像装置 | |
CA2904111C (en) | High-resolution computed tomography | |
CN105158893B (zh) | 基于lcd液晶面板的可编程孔径显微镜***的光场成像方法 | |
CN104390981B (zh) | 基于图像采集装置的快速准确对焦扫描病理切片组织的方法 | |
CN107709968A (zh) | 荧光活检样本成像仪及方法 | |
CN103561655B (zh) | 医用图像诊断装置、医用图像处理装置以及方法 | |
CN109118544A (zh) | 基于透视变换的合成孔径成像方法 | |
CN113267141B (zh) | 一种显微三维信息采集设备 | |
CN104459971A (zh) | 一种基于集成光子芯片的结构光照明显微成像*** | |
CN102822660A (zh) | 断层扫描光场显微镜 | |
CN105704371A (zh) | 一种光场重聚焦方法 | |
CN106908942A (zh) | 基于微透镜阵列的高分辨并行显微成像仪 | |
CN107145052B (zh) | 基于数字插值和相位迭代的全息显微成像方法 | |
CN102512140B (zh) | 一种定位光学投影断层成像旋转中心的方法 | |
CN104713864B (zh) | 基于仿体的光学投影断层成像***的几何校正方法 | |
CN109856169A (zh) | 一种高分辨显微能谱ct成像方法及*** | |
CN104138267B (zh) | 一种x射线成像方法和装置 | |
CN105651699B (zh) | 一种基于面阵相机的动态跟焦方法 | |
CN106934110A (zh) | 一种由聚焦堆栈重建光场的滤波反投影方法和装置 | |
CN106447717B (zh) | 一种基于多角度的选择性光片照明显微成像的重建方法 | |
CN102727188B (zh) | 一种基于拼合螺旋扫描方式的光学投影断层成像方法 | |
US7813529B1 (en) | Optical 3-d surface tomography using depth from focus of partially overlapping 2-d images | |
Zhang et al. | Computational Super-Resolution Imaging With a Sparse Rotational Camera Array | |
Tan et al. | Sharpness preserved sinogram synthesis using convolutional neural network for sparse-view CT imaging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |