CN103308452A - 一种基于景深融合的光学投影断层成像图像获取方法 - Google Patents

一种基于景深融合的光学投影断层成像图像获取方法 Download PDF

Info

Publication number
CN103308452A
CN103308452A CN2013101999831A CN201310199983A CN103308452A CN 103308452 A CN103308452 A CN 103308452A CN 2013101999831 A CN2013101999831 A CN 2013101999831A CN 201310199983 A CN201310199983 A CN 201310199983A CN 103308452 A CN103308452 A CN 103308452A
Authority
CN
China
Prior art keywords
image
matrix
dimensional image
pixel value
image sequence
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.)
Granted
Application number
CN2013101999831A
Other languages
English (en)
Other versions
CN103308452B (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.)
Institute of Automation of Chinese Academy of Science
Original Assignee
Institute of Automation of Chinese Academy of Science
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 Institute of Automation of Chinese Academy of Science filed Critical Institute of Automation of Chinese Academy of Science
Priority to CN201310199983.1A priority Critical patent/CN103308452B/zh
Publication of CN103308452A publication Critical patent/CN103308452A/zh
Application granted granted Critical
Publication of CN103308452B publication Critical patent/CN103308452B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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.按权利要求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所述的方法,其特征在于,比较完最后一幅图像后,得到的目标合成图像像素矩阵对应的图像即为最终得到的全景深图。
CN201310199983.1A 2013-05-27 2013-05-27 一种基于景深融合的光学投影断层成像图像获取方法 Active CN103308452B (zh)

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 true CN103308452A (zh) 2013-09-18
CN103308452B 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)

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104065859A (zh) * 2014-06-12 2014-09-24 青岛海信电器股份有限公司 一种全景深图像的获取方法及摄像装置
CN105118088A (zh) * 2015-08-06 2015-12-02 曲阜裕隆生物科技有限公司 一种基于病理切片扫描装置的3d成像及融合的方法
CN105407280A (zh) * 2015-11-11 2016-03-16 广东欧珀移动通信有限公司 全景图像合成方法和***
CN105451010A (zh) * 2014-08-18 2016-03-30 惠州友华微电子科技有限公司 一种景深获取装置及其获取的方法
CN106339476A (zh) * 2016-08-30 2017-01-18 北京寺库商贸有限公司 一种图像处理方法和***
CN107038724A (zh) * 2015-10-28 2017-08-11 舆图行动股份有限公司 全景鱼眼相机影像校正、合成与景深重建方法与***
CN108495032A (zh) * 2018-03-26 2018-09-04 广东欧珀移动通信有限公司 图像处理方法、装置、存储介质及电子设备
CN108735283A (zh) * 2017-04-13 2018-11-02 佳能株式会社 信息处理装置、***、方法和存储介质
CN110160956A (zh) * 2019-05-27 2019-08-23 广州英特美迪科技有限公司 病理切片扫描影像分析***及其扫描方法
CN111242880A (zh) * 2019-12-30 2020-06-05 广州市明美光电技术有限公司 一种用于显微镜的多景深图像叠加方法、设备及介质
WO2020135040A1 (zh) * 2018-12-29 2020-07-02 中国科学院深圳先进技术研究院 图像三维信息提取方法、对象成像方法、装置及***
CN111861915A (zh) * 2020-07-08 2020-10-30 北京科技大学 一种显微成像场景下失焦扩散效应消除方法及装置
CN114881907A (zh) * 2022-06-30 2022-08-09 江苏集萃苏科思科技有限公司 光学显微图像多景深焦点合成方法、***及图像处理方法
WO2023102724A1 (zh) * 2021-12-07 2023-06-15 宁德时代新能源科技股份有限公司 图像的处理方法和***

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101615289A (zh) * 2009-08-05 2009-12-30 北京优纳科技有限公司 切片组织的三维采集及多层图像融合方法
US20110193953A1 (en) * 2010-02-05 2011-08-11 Applied Vision Company, Llc System and method for estimating the height of an object using tomosynthesis-like techniques
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射线发光断层成像***及其方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101615289A (zh) * 2009-08-05 2009-12-30 北京优纳科技有限公司 切片组织的三维采集及多层图像融合方法
US20110193953A1 (en) * 2010-02-05 2011-08-11 Applied Vision Company, Llc System and method for estimating the height of an object using tomosynthesis-like techniques
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)

* Cited by examiner, † Cited by third party
Title
王会峰: "旋转扫描序列图像的全景图拼接方法", 《光学学报》 *

Cited By (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104065859A (zh) * 2014-06-12 2014-09-24 青岛海信电器股份有限公司 一种全景深图像的获取方法及摄像装置
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成像及融合的方法
CN107038724A (zh) * 2015-10-28 2017-08-11 舆图行动股份有限公司 全景鱼眼相机影像校正、合成与景深重建方法与***
CN105407280A (zh) * 2015-11-11 2016-03-16 广东欧珀移动通信有限公司 全景图像合成方法和***
CN106339476B (zh) * 2016-08-30 2019-10-29 北京寺库商贸有限公司 一种图像处理方法和***
CN106339476A (zh) * 2016-08-30 2017-01-18 北京寺库商贸有限公司 一种图像处理方法和***
CN108735283A (zh) * 2017-04-13 2018-11-02 佳能株式会社 信息处理装置、***、方法和存储介质
CN108495032A (zh) * 2018-03-26 2018-09-04 广东欧珀移动通信有限公司 图像处理方法、装置、存储介质及电子设备
CN108495032B (zh) * 2018-03-26 2020-08-04 Oppo广东移动通信有限公司 图像处理方法、装置、存储介质及电子设备
CN111381357A (zh) * 2018-12-29 2020-07-07 中国科学院深圳先进技术研究院 图像三维信息提取方法、对象成像方法、装置及***
WO2020135040A1 (zh) * 2018-12-29 2020-07-02 中国科学院深圳先进技术研究院 图像三维信息提取方法、对象成像方法、装置及***
CN111381357B (zh) * 2018-12-29 2021-07-20 中国科学院深圳先进技术研究院 图像三维信息提取方法、对象成像方法、装置及***
CN110160956A (zh) * 2019-05-27 2019-08-23 广州英特美迪科技有限公司 病理切片扫描影像分析***及其扫描方法
CN111242880A (zh) * 2019-12-30 2020-06-05 广州市明美光电技术有限公司 一种用于显微镜的多景深图像叠加方法、设备及介质
CN111861915A (zh) * 2020-07-08 2020-10-30 北京科技大学 一种显微成像场景下失焦扩散效应消除方法及装置
CN111861915B (zh) * 2020-07-08 2021-08-13 北京科技大学 一种显微成像场景下失焦扩散效应消除方法及装置
WO2023102724A1 (zh) * 2021-12-07 2023-06-15 宁德时代新能源科技股份有限公司 图像的处理方法和***
US11948287B2 (en) 2021-12-07 2024-04-02 Contemporary Amperex Technology Co., Limited Image processing method and system
CN114881907A (zh) * 2022-06-30 2022-08-09 江苏集萃苏科思科技有限公司 光学显微图像多景深焦点合成方法、***及图像处理方法
CN114881907B (zh) * 2022-06-30 2022-09-23 江苏集萃苏科思科技有限公司 光学显微图像多景深焦点合成方法、***及图像处理方法

Also Published As

Publication number Publication date
CN103308452B (zh) 2015-05-06

Similar Documents

Publication Publication Date Title
CN103308452B (zh) 一种基于景深融合的光学投影断层成像图像获取方法
CN101896856B (zh) 用于光学层析成像的景深扩展
CN104390981B (zh) 基于图像采集装置的快速准确对焦扫描病理切片组织的方法
CN109413407B (zh) 高空间分辨率光场采集装置与图像生成方法
US9088729B2 (en) Imaging apparatus and method of controlling same
CN102980664B (zh) 超像素微扫描方法及相应的红外超分辨率实时成像装置
CN101441320B (zh) 一种基于显微成像检测的高动态图像获取装置及其方法
CN105158893B (zh) 基于lcd液晶面板的可编程孔径显微镜***的光场成像方法
CN107065159A (zh) 一种基于大照明数值孔径的大视场高分辨率显微成像装置及迭代重构方法
CN105704371B (zh) 一种光场重聚焦方法
JP2017063414A (ja) ライト・フィールドを表すデータを生成する方法および装置
CN109118544A (zh) 基于透视变换的合成孔径成像方法
CN104168430B (zh) 一种tdi ccd焦平面编码超分辨率成像装置及方法
CN106908942A (zh) 基于微透镜阵列的高分辨并行显微成像仪
CN113267141B (zh) 一种显微三维信息采集设备
CN102822660A (zh) 断层扫描光场显微镜
KR20180053668A (ko) 화소 빔을 나타내는 데이터를 생성하기 위한 방법 및 장치
CN103744172A (zh) 一种具备空间光调制照明的共聚焦显微成像方法
CN108765285A (zh) 一种基于视频清晰度融合的大尺寸显微图像生成方法
CN108537862A (zh) 一种自适应降噪的傅里叶衍射扫描显微镜成像方法
RU2580870C2 (ru) Система получения изображения высокого разрешения
CN106530226A (zh) 获取高分辨率高清晰工业图像的实现方法
CN105651699B (zh) 一种基于面阵相机的动态跟焦方法
CN111080774B (zh) 一种应用深度采样进行光场重构的方法及***
CN104143179B (zh) 一种多线列时差扫描扩展采样的运动目标增强方法

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