CN108986084B - 一种新型oct图像显示方法 - Google Patents
一种新型oct图像显示方法 Download PDFInfo
- Publication number
- CN108986084B CN108986084B CN201810691602.4A CN201810691602A CN108986084B CN 108986084 B CN108986084 B CN 108986084B CN 201810691602 A CN201810691602 A CN 201810691602A CN 108986084 B CN108986084 B CN 108986084B
- Authority
- CN
- China
- Prior art keywords
- signal
- image
- time
- original data
- dimensional image
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10101—Optical tomography; Optical coherence tomography [OCT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20048—Transform domain processing
- G06T2207/20056—Discrete and fast Fourier transform, [DFT, FFT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30041—Eye; Retina; Ophthalmic
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明涉及一种新型OCT图像显示方法,包括有以下步骤:S1.将OCT设备所获得的原始数据进行傅里叶变换,转换成A‑扫描三维图像并实时显示输出到信号瀑布图中;S2.将OCT设备所获得的原始数据进行傅里叶变换后与手动转动的角度θ结合转换成B‑扫描二维图像并实时显示输出到扇形图中;S3.结合信号瀑布图和扇形图同时显示,让医护相关人员实现定点精准检测组织。
Description
技术领域
本发明涉及OCT图像处理领域,更具体地,涉及一种新型OCT图像显示方法。
背景技术
光学相干断层成像(Optical Coherence Tomography,简称OCT),已广泛应用在眼科诊断领域,这项技术是建立在光学、电子学以及计算机技术科学的基础上,是集光电、高速数据采集和图像处理等多项前沿学科为一体的新型成像技术,OCT凭借其具有高分辨率、高速成像等优点而备受人们的关注,并在生物医学与临床诊断领域开始得到重视和应用。
与现有的CT、超声、MRI等其他成像方式相比,OCT具有极高的分辨率,与传统的激光共聚焦显微镜相比,OCT的成像深度具有明显的优势。传统光学探头的核心技术大多采用光纤束进行光传导并进行成像,或者采用CCD技术进行成像,此类内窥探头仅能探测组织表面的病变,然而早期癌症的症状发生在表皮以下1-3毫米的深度,因此传统光学内窥探头就显得力不从心。目前也有通过超声原理进行医学成像的内窥探头,虽然可获得生物组织表层以下较深的组织信息,但分辨率仅为毫米量级,对早期的癌症易造成漏诊。因此,OCT技术的发展对医学领域有重大的意义。
目前广泛应用在人体体内的OCT技术中所利用的成像方法都是依靠电机驱动进入体内的OCT探头进行扇形遍历式扫描。采用此扫描运动方法的缺点有:1)对于探头有一定的要求;2)导致成像结果里存在着假象;3)由于微型电机的存在,硬件成本高;4)图像质量被严重压缩,有效信息被损坏。
为此,本发明提出一种低成本、速度快的定点精准检测的OCT图像显示方法。
发明内容
为了克服背景技术的不足,本发明提供了一种新型OCT图像显示方法,具有对定点精准测量、探头即插即用、实时三维成像、图像失真少和成本低等优点。本发明主要用于对人体内部组织的无创、实时、高分辨率三维成像,其分辨率为微米量级,可应用于人体组织狭窄空间和细小腔道的早期病变检测。
为实现以上发明目的,采用的技术方案是:
一种新型OCT图像显示方法,包括有以下步骤:
S1.将OCT设备所获得的原始数据进行傅里叶变换,转换成A-扫描三维图像并实时显示输出到信号瀑布图中;
S2.将OCT设备所获得的原始数据进行傅里叶变换后与手动转动的角度θ结合转换成B-扫描二维图像并实时显示输出到扇形图中;
S3.结合信号瀑布图和扇形图同时显示,让医护相关人员实现定点精准检测组织。
优选地,所述步骤S3将B-扫描二维图像进行归一化变换和坐标变换并以灰度值显示信号强弱,从而实时显示到扇形图中。
优选地,所述信号瀑布图为二维的图像,其Y轴代表实时深度,其灰度值代表实时干涉信号的强弱,其X轴代表时刻t,一个时刻t只有一条曲线。
优选地,所述信号瀑布图为三维的图像,其X轴代表实时深度,其Y轴代表实时干涉信号的强弱,其Z轴代表时刻t,一个时刻t只有一条曲线。
附图说明
图1为实施例1的新型OCT图像显示方法的流程图。
图2为实施例1的新型OCT图像显示方法的二维信号瀑布图。
图3为实施例1的新型OCT图像显示方法的三维信号瀑布图。
图4为实施例1的新型OCT图像显示方法的扇形图。
具体实施方式
附图仅用于示例性说明,不能理解为对本专利的限制;
以下结合附图和实施例对本发明做进一步的阐述。
实施例1
如图1所示,本发明提供了一种新型OCT图像显示方法,包括有以下步骤:
S1.将OCT设备所获得的原始数据进行傅里叶变换,转换成A-扫描三维图像并实时显示输出到信号瀑布图中;
S2.将OCT设备所获得的原始数据进行傅里叶变换后与手动转动的角度θ结合转换成B-扫描二维图像并实时显示输出到扇形图中;
S3.结合信号瀑布图和扇形图同时显示,让医护相关人员实现定点精准检测组织。
在一个优选地实施例中,信号瀑布图平面示意图如附图2所示,信号瀑布图是一个二维图像,信号瀑布图Y轴代表的数值为实时深度,信号瀑布图的灰度值数值为实时干涉信号的强弱,信号瀑布图X轴代表的数值代表着时刻t,一个时刻t只有一条曲线。信号瀑布图可以显示多个时间内OCT设备所获得的信号并进行对比比较。信号瀑布图可以实时显示并定点显示所测物体的深度信息。
扇形图与二维信号瀑布图存在坐标变换变换公式反映了第t时刻,此时对应展开的角度为θt,此时二维信号瀑布图所呈现横坐标xt,即为扇形图的角度θt。二维信号瀑布图的纵坐标yt即为扇形图中半径rt。二维信号瀑布图的灰度值对应的在扇形图中也以灰度值显示。结合信号瀑布图和扇形图同时显示,可以让使用的医护相关人员实现定点精准检测组织。
在一个优选的实施例中,信号瀑布图的平面示意图如附图3所示,信号瀑布图可以用三维图像表示,信号瀑布图X轴代表的数值为实时深度,信号瀑布图Y轴代表的数值为实时干涉信号的强弱,信号瀑布图Z轴代表的数值代表着时刻t,一个时刻t只有一条曲线。信号瀑布图可以显示多个时间内OCT设备所获得的信号并进行对比比较。信号瀑布图可以实时显示并定点显示所测物体的深度信息。
扇形图与三维信号瀑布图存在映射关系:映射关系反映了第t时刻,此时对应展开的角度为θt,此时三维图所呈现横坐标xt,即为扇形图当角度为θt时的半径rt。三维图的纵坐标yt在扇形图中则以灰度值大小呈现,纵坐标yt越大,则灰度值越大,反之亦然。结合信号瀑布图和扇形图同时显示,可以让使用的医护相关人员实现定点精准检测组织
扇形图的平面示意图如附图4所示,扇形图是一个二维图像,扇形图展开的角度θ即人工手动转动的角度θ,扇形图的每一个θ上的信息即为该θ数据下的A-扫描数据。B-扫描二维图像进行坐标变换、归一化、压缩处理后,将信号强弱以灰度值大小表示,从而得到扇形图。扇形图可以显示一定角度下的物体深度信息。扇形图可以针对异常部位手动反复测量;
压缩处理即是对原始数据进行坐标变换为扇形图时,对数据处理的一种方式。在变换时,以扫描点为圆心,工作距离为半径,随着工作距离越大,能显示的数据增多;反之,随着工作距离越短,能显示的数据减少。因此,距离圆心越远的即工作距离越长的数据点保留得越多,反之,越靠近圆心的即工作距离越短的数据点被舍弃得越多。
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明权利要求的保护范围之内。
Claims (3)
1.一种新型OCT图像显示方法,其特征在于:包括有以下步骤:
S1.将OCT设备所获得的原始数据进行傅里叶变换,转换成A-扫描三维图像并实时显示输出到信号瀑布图中;
S2.将OCT设备所获得的原始数据进行傅里叶变换后与手动转动的角度θ结合转换成B-扫描二维图像并实时显示输出到扇形图中;
S3.结合信号瀑布图和扇形图同时显示,让医护相关人员实现定点精准检测组织;
所述步骤S3将B-扫描二维图像进行归一化变换和坐标变换并以灰度值显示信号强弱,从而实时显示到扇形图中。
2.根据权利要求1所述的新型OCT图像显示方法,其特征在于:所述信号瀑布图为二维的图像,其Y轴代表实时深度,其灰度值代表实时干涉信号的强弱,其X轴代表时刻t,一个时刻t只有一条曲线。
3.根据权利要求1所述的新型OCT图像显示方法,其特征在于:所述信号瀑布图为三维的图像,其X轴代表实时深度,其Y轴代表实时干涉信号的强弱,其Z轴代表时刻t,一个时刻t只有一条曲线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810691602.4A CN108986084B (zh) | 2018-06-28 | 2018-06-28 | 一种新型oct图像显示方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810691602.4A CN108986084B (zh) | 2018-06-28 | 2018-06-28 | 一种新型oct图像显示方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108986084A CN108986084A (zh) | 2018-12-11 |
CN108986084B true CN108986084B (zh) | 2022-04-19 |
Family
ID=64538837
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810691602.4A Active CN108986084B (zh) | 2018-06-28 | 2018-06-28 | 一种新型oct图像显示方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108986084B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112348907B (zh) * | 2020-10-16 | 2024-03-26 | 佛山科学技术学院 | 一种用于光学相干层析成像三维数据的压缩方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104376549A (zh) * | 2014-11-20 | 2015-02-25 | 华北电力大学(保定) | 血管内超声与血管内oct图像的融合方法 |
CN106934761A (zh) * | 2017-02-15 | 2017-07-07 | 苏州大学 | 一种三维非刚性光学相干断层扫描图像的配准方法 |
CN106963338A (zh) * | 2017-03-31 | 2017-07-21 | 杭州电子科技大学 | 晶状体生物力学和光学特性无创在体成像系统及测量方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102835948B (zh) * | 2012-09-12 | 2016-03-30 | 无锡微奥科技有限公司 | 一种扫频光源oct实时图像显示方法及其系统 |
US9241626B2 (en) * | 2013-03-14 | 2016-01-26 | Carl Zeiss Meditec, Inc. | Systems and methods for improved acquisition of ophthalmic optical coherence tomography data |
CN104825120A (zh) * | 2015-05-08 | 2015-08-12 | 南京微创医学科技有限公司 | Oct内窥扫描成像系统中使用的光学时钟信号发生系统和方法 |
CN104794740A (zh) * | 2015-05-08 | 2015-07-22 | 南京微创医学科技有限公司 | 利用通用图像处理器处理oct信号的方法及系统 |
JP6702764B2 (ja) * | 2016-03-08 | 2020-06-03 | キヤノン株式会社 | 光干渉断層データの処理方法、該方法を実行するためのプログラム、及び処理装置 |
CN106073788B (zh) * | 2016-07-19 | 2023-03-17 | 杭州捷诺飞生物科技有限公司 | 基于oct的原位三维打印皮肤修复设备及其实现方法 |
CN107862661B (zh) * | 2017-11-06 | 2019-11-29 | 郑州轻工业学院 | 一种光学相干层析成像系统图像校正方法 |
-
2018
- 2018-06-28 CN CN201810691602.4A patent/CN108986084B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104376549A (zh) * | 2014-11-20 | 2015-02-25 | 华北电力大学(保定) | 血管内超声与血管内oct图像的融合方法 |
CN106934761A (zh) * | 2017-02-15 | 2017-07-07 | 苏州大学 | 一种三维非刚性光学相干断层扫描图像的配准方法 |
CN106963338A (zh) * | 2017-03-31 | 2017-07-21 | 杭州电子科技大学 | 晶状体生物力学和光学特性无创在体成像系统及测量方法 |
Non-Patent Citations (1)
Title |
---|
基于OCT图像的房角开放距离计算机辅助测量研究;武薇;《万方学位论文全文数据库》;20131008;第1-155页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108986084A (zh) | 2018-12-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Dima et al. | Non-invasive carotid imaging using optoacoustic tomography | |
JP6160487B2 (ja) | 超音波診断装置およびその制御方法 | |
Aspres et al. | Imaging the skin | |
CN103040429B (zh) | 一种用于口腔的光学影像检测装置和成像方法 | |
US10278663B2 (en) | Sensor coordinate calibration in an ultrasound system | |
Miranda et al. | Intrauterine photoacoustic and ultrasound imaging probe | |
CN113096056A (zh) | 一种基于区域互补的血管内图像融合方法 | |
CN108986084B (zh) | 一种新型oct图像显示方法 | |
CN107478414B (zh) | 一种oct成像回抽性能测试装置及方法 | |
CN203153692U (zh) | 一种用于口腔的光学影像检测装置 | |
CN102144931B (zh) | 一种类显微三维超声成像法 | |
CN106794000A (zh) | 超声图像处理方法及其超声成像装置 | |
CN111281426A (zh) | 表浅神经的超声成像及图像分析方法 | |
Nelson et al. | Multiscale label-free imaging of fibrillar collagen in the tumor microenvironment | |
CN210055992U (zh) | 一种血管内双模成像装置 | |
CN111493832A (zh) | 一种基于En face-OCT的内窥成像方法 | |
CN202426507U (zh) | 一体化oct硬质脑室镜系统 | |
Mahmoud et al. | High-resolution ultrasound imaging for jawbone surface | |
Jeng et al. | Ultrasound shear-wave computed tomography for elasticity imaging | |
Wells | Advances in ultrasound: from microscanning to telerobotics. | |
JP2009100996A (ja) | 超音波診断システム | |
Liu et al. | Seeing beyond skin deep: High‐resolution ultrasound in dermatology—A comprehensive review and future prospects | |
CN118134931B (zh) | 一种超声图像处理方法、设备及系统 | |
WO2022243714A1 (en) | Depth-surface imaging device for registering ultrasound images to each other and to surface images by using surface information | |
CN202437073U (zh) | 一体化oct硬质肛肠镜系统 |
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 |