CN114909987B - 一种b扫描分割法抑制样品整体抖动的方法 - Google Patents

一种b扫描分割法抑制样品整体抖动的方法 Download PDF

Info

Publication number
CN114909987B
CN114909987B CN202210501036.2A CN202210501036A CN114909987B CN 114909987 B CN114909987 B CN 114909987B CN 202210501036 A CN202210501036 A CN 202210501036A CN 114909987 B CN114909987 B CN 114909987B
Authority
CN
China
Prior art keywords
scan
sample
scans
sub
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
Application number
CN202210501036.2A
Other languages
English (en)
Other versions
CN114909987A (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.)
Southeast University
Original Assignee
Southeast 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 Southeast University filed Critical Southeast University
Priority to CN202210501036.2A priority Critical patent/CN114909987B/zh
Publication of CN114909987A publication Critical patent/CN114909987A/zh
Application granted granted Critical
Publication of CN114909987B publication Critical patent/CN114909987B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00—Measuring instruments characterised by the use of optical techniques
    • G01B9/02—Interferometers
    • G01B9/0209—Low-coherence interferometers
    • G01B9/02091—Tomographic interferometers, e.g. based on optical coherence
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00—Measuring instruments characterised by the use of optical techniques
    • G01B9/02—Interferometers
    • G01B9/02083—Interferometers characterised by particular signal processing and presentation

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明提供了一种B扫描分割法抑制样品整体抖动的方法。将信息采集过程中的每个完整B扫描分为几个子B扫描，每个子B扫描在样本光束移动到下一个子B扫描之前重复多次，将所有子B扫描拼接成完整B扫描。分析连续B扫描之间的时间间隔对用最近提出的对数标度强度差标准差算法计算的微血管图像的影响，找到一个最佳的时间间隔，能有效提高图像信噪比.该技术可以应用于需要在同一位置进行多次扫描的光学相干断层扫描成像技术，如基于OCT的微血流图像和基于高分辨OCT的细胞动态成像技术。本发明可以有效抑制图像采集过程中的样品整体抖动，可以在总的数据采集时间几乎不变的前提下，提高成像质量，可广泛应用于光学相干断层扫描技术领域。

Description

一种B扫描分割法抑制样品整体抖动的方法
技术领域
本发明属于光学相干断层扫描(OCT)技术领域，尤其涉及一种B扫描分割法抑制样品整体抖动的方法。
背景技术
光学相干断层成像(OCT)是20世纪90年代提出的一种成像方式，具有高分辨率和非侵入性的优点。由于细胞中运动的细胞器以及组织内的红细胞(RBC)等都可以调制后向散射光信号，因此基于这种机制的多功能OCT被提出用于实现血流流速测量或细胞动态成像，如基于高分辨OCT的细胞动态成像和OCT血管成像(OCTA)。
基于高分辨OCT的细胞动态成像是最近提出的一种亚细胞动态成像技术，可以获得与组织切片相媲美的组织图像。OCTA是临床应用中最广泛使用的无创重建微血管图的技术之一，用于疾病诊断和治疗监测。现在已经提出了基于OCT振幅、相位或复杂信号的多种OCTA算法。
在动态成像技术中，为了提高动态信号检测的灵敏度，通常利用同一位置的多次B扫描(而不是A扫描)来计算动态信号，这会产生更大的时间间隔。由于在数据采集过程中，被测对象的心跳和呼吸以及动态细胞的组织弹性和环境振动等引起的样品整体抖动是不可避免的，并且可以用移动散射体的相同机制调制后向散射光信号，因此运动伪影总是出现在获得的动态图像上。克服这一问题的一种常用方法是：在后期数据处理过程中，使用图像配准算法来补偿体积偏移。然而，高精度的图像配准会大大增加数据处理的计算复杂度，这使得实时筛选特征具有挑战性。
发明内容
本发明目的在于提供一种B扫描分割法抑制样品整体抖动的方法,以解决高精度的图像配准会大大增加数据处理的计算复杂度，这使得实时筛选特征具有挑战性的技术问题。
为解决上述技术问题，本发明的具体技术方案如下：
一种B扫描分割法抑制样品整体抖动的方法，包括以下步骤：
步骤1、在使用光学相干断层扫描OCT系统对探查对象进行动态信号提取的过程中，用B扫描分割法进行扫描；
步骤2、基于B扫描分割法控制得到连续B扫描之间不同的时间间隔，测试不同时间间隔对于对数标度强度差标准差算法计算的微血管图像的对比噪声比影响；
步骤3、采用最优时间间隔扫描方案进行扫描，无需进行图像配准即可获得抑制了样品整体抖动的高对比度微血管图像和细胞动态图像。
进一步的，步骤1具体包括以下步骤：
在B扫描分段动态显微光学相干断层成像中，每个完整的B扫描被分为几个子B扫描，每个子B扫描在样本光束移动到下一个子B扫描之前重复多次；
当一个完整B扫描的所有子B扫描完成后，我们将所有子B扫描拼接到具有相同编号的完整B扫描中。
进一步的，当一个完整B扫描的所有子B扫描完成后，我们将所有子B扫描拼接到具有相同编号的完整B扫描中，具体包括以下步骤：
为了实现多个子B扫描图像的完美拼接，两个相邻子B扫描的振镜驱动电压要重叠，并在每个子B扫描的每个边缘采集8个额外的A扫描以补偿振镜机械惯性造成的模糊部分。
进一步的，步骤2具体包括以下步骤：
每个B扫描由若干个A扫描组成，因此通过跳过不同数量的A扫描即可得到不同的时间间隔；
计算对数标度强度差标准差信号与不同时间间隔的关系，得到对数标度强度差标准差的图像；
在对数标度强度差标准差图像中，选取两个区域标记用于对比噪声比的计算，得到对比噪声比与时间间隔的关系图。
进一步的，在对数标度强度差标准差图像中，选取两个区域标记用于对比噪声比计算，包括以下步骤：
要标记的两个区域，应分别在动态区域和背景区域中选取。
进一步的，所述的对数标度强度差标准差信号与不同时间间隔的关系根据下式计算得到：
其中DI＝log(Im)–log(Im+1)表示两个相邻A扫描的差分信号，N是深度方向上的窗口大小，是窗口内强度的平均值；
所述的微血管图像的对比噪声比根据下式计算得到：
其中是动态区域的平均值,/>是背景区域的平均值，Ibn是背景区域的标准差。
进一步的，步骤3具体包括以下步骤：
对比噪声比与时间间隔的关系图的峰值位置对应的时间间隔就是最佳时间间隔。
进一步的，步骤3具体包括以下步骤：
通过在快速扫描方向上的平均对数标度强度差标准差图像来定量评估抑制样品整体抖动的性能，选取不同时间间隔的扫描方案，计算得到在不同时间间隔下的平均对数标度强度差标准差值曲线图，平均对数标度强度差标准差值整体越低的扫描方案表明该方案的样品整体抖动抑制效果越好。
本发明的一种B扫描分割法抑制样品整体抖动的方法，具有以下优点：
1、本发明采用了B扫描分割法，两次完整B扫描之间的时间间隔可以减少数倍，同时总的数据采集时间几乎保持不变，从而有效抑制样品整体抖动和改善图像质量。
2、本发明通过控制每个完整B扫描划分成的子B扫描数量，调制时间间隔，采用最佳时间间隔扫描方案，进一步抑制了样品整体抖动，无需图像配准即可直接获得高对比度的微血管图像，有效提高了数据处理速度。
附图说明
图1为本发明实施例提供的一种B扫描分割法抑制样品整体抖动的方法步骤流程图；
图2为本发明的实施例提供的微血管图像对比噪声比(CNR)与时间间隔的关系图；
图3(a)为本发明实施例提供的在皮肤表面以下225-375μm深度范围内，以16.0ms的时间间隔拍摄的人体皮肤正面微血管图像；
图3(b)为本发明实施例提供的在皮肤表面以下225-375μm深度范围内，以8.2ms的时间间隔拍摄的人体皮肤正面微血管图像；
图3(c)为本发明实施例提供的在皮肤表面以下225-375μm深度范围内，以4.3ms的时间间隔扫描拍摄的人体皮肤正面微血管图像；
图3(d)为本发明实施例提供的在皮肤表面以下375-525μm深度范围内，以16.0ms的时间间隔扫描拍摄的人体皮肤正面微血管图像；
图3(e)为本发明实施例提供的在皮肤表面以下375-525μm深度范围内，以8.2ms的时间间隔扫描拍摄的人体皮肤正面微血管图像；
图3(f)为本发明实施例提供的在皮肤表面以下375-525μm深度范围内，以4.3ms的时间间隔扫描拍摄的人体皮肤正面微血管图像；
图4为本发明实施例提供的不同时间间隔下的平均对数标度强度差标准差曲线图；
图5(a)为本发明实施例提供的传统扫描方案获得的斑马鱼肠道的对比度图；
图5(b)为本发明实施例提供的B扫描分割法获得的斑马鱼肠道的对比度图；
图6为本发明实施例提供的两种方法的第一次B扫描和后续199次B扫描之间的互相关系数曲线图。
具体实施方式
为了更好地了解本发明的目的、结构及功能，下面结合附图，对本发明一种B扫描分割法抑制样品整体抖动的方法做进一步详细的描述。
在本发明的描述中，多个的含义是两个或两个以上，如果有描述到第一、第二只是用于区分技术特征为目的。而不能理解为指示或者暗示相对重要性或者隐含指明所指示的技术特征的数量或者隐含知名所指示的技术特征的先后关系。此外，除非另有定义，本文所使用的所有的技术和科学术语与本技术领域的技术人员通常理解的含义相同。
本发明的一种B扫描分割法抑制样品整体抖动的方法，包括以下步骤：
步骤1：在使用光学相干断层扫描OCT系统对探查对象进行动态信号提取的过程中，用B扫描分割法进行扫描；
在B扫描分段动态显微光学相干断层成像(BD-MOCT)中，每个完整的B扫描被分为几个子B扫描，每个子B扫描在样本光束移动到下一个子B扫描之前重复多次；
当一个完整B扫描的所有子B扫描完成后，我们将所有子B扫描拼接到具有相同编号的完整B扫描中：为了实现多个子B扫描图像的完美拼接，两个相邻子B扫描的振镜驱动电压要略微重叠，并在每个子B扫描的每个边缘采集8个额外的A扫描以补偿振镜机械惯性造成的模糊部分；
通过这种方式，两次连续拼接的B扫描之间的时间间隔可以减少几倍，从而有效抑制样品整体抖动。
步骤2：基于B扫描分割法控制得到连续B扫描之间不同的时间间隔，测试不同时间间隔对于对数标度强度差标准差(DSDLI)算法计算的微血管图像的对比噪声比(CNR)影响；
每个B扫描由若干个A扫描组成，因此通过跳过不同数量的A扫描即可得到不同的时间间隔；
计算对数标度强度差标准差(DSDLI)信号与不同时间间隔的关系：
其中DI＝log(Im)–log(Im+1)表示两个相邻A扫描的差分信号，N是深度方向上的窗口大小，是窗口内强度的平均值；
在DSDLI图像中，选取两个区域(分别在动态区域和背景区域中选取)标记用于CNR计算，得到CNR与时间间隔的关系图：所述的微血管图像的对比噪声比(CNR)根据下式计算得到：
其中是动态区域的平均值,/>和Ibn分别是背景区域的平均值和标准差。
步骤3：采用最优时间间隔扫描方案进行扫描，无需进行图像配准即可获得抑制了样品整体抖动的高对比度微血管图像。
CNR与时间间隔的关系图峰值位置对应的时间间隔就是最佳时间间隔；
还可以通过在快速扫描方向上的平均DSDLI图像来定量评估抑制样品整体抖动的性能，选取不同时间间隔的扫描方案，计算得到在不同时间间隔下的平均DSDLI值曲线图，平均DSDLI值整体越低的扫描方案表明该方案的样品整体抖动抑制效果越好。
通过计算不同时间间隔的B扫图像的相关系数来量化衡量样品抖动噪声对结果的影响，可以发现B扫分割方法可以更好地保持同一位置上连续B扫的相关性，意味着样品抖动噪声对结果有更小的影响。
下面是本发明的第一实施例：
参照图1，本发明第一实施例是其在微血流成像技术方面的应用，提供了一种B扫描分割法抑制样品整体抖动的方法，包括以下步骤：
步骤1：在使用光学相干断层扫描OCT系统对健康志愿者手指上1.5×1.5mm的区域进行动态信号提取的过程中，用B扫描分割法进行扫描；
每个完整的B扫描有600个有效A扫描组成，每个子B扫描在样本光束移动到下一个子B扫描之前重复多次；
当一个完整B扫描的所有子B扫描完成后，我们将所有子B扫描拼接到具有相同编号的完整B扫描中：为了实现多个子B扫描图像的完美拼接，两个相邻子B扫描的振镜驱动电压要略微重叠，并在每个子B扫描的每个边缘采集8个额外的A扫描以补偿振镜机械惯性造成的模糊部分。
步骤2：基于B扫描分割法控制得到(连续B扫描之间)不同的时间间隔，测试不同时间间隔对于对数标度强度差标准差(DSDLI)算法计算的微血管图像的对比噪声比(CNR)影响；
通过跳过不同数量的A扫描，以得到不同的时间间隔；
计算对数标度强度差标准差(DSDLI)信号与不同时间间隔的关系：
其中DI＝log(Im)–log(Im+1)表示两个相邻A扫描的差分信号，N是深度方向上的窗口大小，是窗口内强度的平均值；
在DSDLI图像中，选取两个区域(分别在动态区域和背景区域中选取)标记用于CNR计算，得到CNR与时间间隔的关系图：所述的微血管图像的对比噪声比(CNR)根据下式计算得到：
其中是动态区域的平均值,/>和Ibn分别是背景区域的平均值和标准差。
步骤3：采用最优时间间隔扫描方案进行扫描，无需进行图像配准即可获得抑制了样品整体抖动的高对比度微血管图像。
参照图2，本发明实施例中4.3ms接近峰值的间隔，即是最佳时间间隔；
图3(a)-(c)为在皮肤表面以下225-375μm深度范围内，分别以16.0ms、8.2ms和4.3ms的时间间隔拍摄的人类皮肤正面微血管图像；(d)-(f)为在皮肤表面以下375-525μm深度范围内，分别以16.0ms、8.2ms和4.3ms的时间间隔拍摄的人体皮肤正面微血管图像。
参照图4，本发明实施例还可以通过在快速扫描方向上的平均DSDLI图像来定量评估抑制样品整体抖动的性能，选取不同时间间隔的扫描方案，计算得到在不同时间间隔下的平均DSDLI值曲线图，平均DSDLI值整体越低的扫描方案表明该方案的样品整体抖动抑制效果越好。显然本发明实施例中，时间间隔为4.3ms的扫描方案对于样品整体抖动的抑制效果最佳。
下面是本发明的第二实施例：
本发明的第二实施例为其在基于高分辨OCT的细胞动态成像上的应用，提供了一种B扫描分割法抑制样品整体抖动的方法，包括以下步骤：
步骤一：提取一条斑马鱼的肠道，提取后立刻将肠浸入细胞培养基中；
步骤二：参照图5(b)，采用本发明的方法对斑马鱼肠道进行扫描成像，具体实施方式与上一实施例步骤1、步骤2、步骤3相同，在此就不再赘述；
步骤三：参照图5(a)，不采用本发明的方法，采用传统扫描方案对斑马鱼肠道进行扫描。对比图5(a)和(b)，可以看出采用本发明的扫描方案扫描，所成的像具有更丰富的动态信息；
步骤四：参照图6，分别计算采用和不采用本发明下，第一次B扫描与后续199次B扫描之间的相关系数以进行定量比较。可以看出，采用本发明方法的相关性下降比传统方法慢得多，这说明样品整体抖动被有效抑制了，并且利用B扫描分割法后获得的动态图像中不同组织细胞间具有更高的对比度。
可以理解，本发明是通过一些实施例进行描述的，本领域技术人员知悉的，在不脱离本发明的精神和范围的情况下，可以对这些特征和实施例进行各种改变或等效替换。另外，在本发明的教导下，可以对这些特征和实施例进行修改以适应具体的情况及材料而不会脱离本发明的精神和范围。因此，本发明不受此处所公开的具体实施例的限制，所有落入本申请的权利要求范围内的实施例都属于本发明所保护的范围内。

Claims (7)

1.一种B扫描分割法抑制样品整体抖动的方法，其特征在于，包括以下步骤：
步骤1、在使用光学相干断层扫描OCT系统对探查对象进行动态信号提取的过程中，用B扫描分割法进行扫描；
步骤2、基于B扫描分割法控制得到连续B扫描之间不同的时间间隔，测试不同时间间隔对于对数标度强度差标准差算法计算的微血管图像的对比噪声比影响；
步骤3、采用最优时间间隔扫描方案进行扫描，无需进行图像配准即可获得抑制了样品整体抖动的高对比度微血管图像和细胞动态图像；
所述步骤2具体包括以下步骤：
每个B扫描由若干个A扫描组成，因此通过跳过不同数量的A扫描即可得到不同的时间间隔；
计算对数标度强度差标准差信号与不同时间间隔的关系，得到对数标度强度差标准差的图像；
在对数标度强度差标准差图像中，选取两个区域标记用于对比噪声比的计算，得到对比噪声比与时间间隔的关系图。
2.根据权利要求1所述的B扫描分割法抑制样品整体抖动的方法，其特征在于，步骤1具体包括以下步骤：
在B扫描分段动态显微光学相干断层成像中，每个完整的B扫描被分为几个子B扫描，每个子B扫描在样本光束移动到下一个子B扫描之前重复多次；
当一个完整B扫描的所有子B扫描完成后，我们将所有子B扫描拼接到具有相同编号的完整B扫描中。
3.根据权利要求2所述的B扫描分割法抑制样品整体抖动的方法，其特征在于，当一个完整B扫描的所有子B扫描完成后，我们将所有子B扫描拼接到具有相同编号的完整B扫描中，具体包括以下步骤：
为了实现多个子B扫描图像的完美拼接，两个相邻子B扫描的振镜驱动电压要重叠，并在每个子B扫描的每个边缘采集8个额外的A扫描以补偿振镜机械惯性造成的模糊部分。
4.根据权利要求1所述的B扫描分割法抑制样品整体抖动的方法，其特征在于，在对数标度强度差标准差图像中，选取两个区域标记用于对比噪声比计算，包括以下步骤：
要标记的两个区域，应分别在动态区域和背景区域中选取。
5.根据权利要求1所述的B扫描分割法抑制样品整体抖动的方法，其特征在于，所述的对数标度强度差标准差信号与不同时间间隔的关系根据下式计算得到：
其中DI＝log(Im)–log(Im+1)表示两个相邻A扫描的差分信号，N是深度方向上的窗口大小，是窗口内强度的平均值；
所述的微血管图像的对比噪声比根据下式计算得到：
其中是动态区域的平均值,/>是背景区域的平均值，Ibn是背景区域的标准差。
6.根据权利要求1所述的B扫描分割法抑制样品整体抖动的方法，其特征在于，步骤3具体包括以下步骤：
对比噪声比与时间间隔的关系图的峰值位置对应的时间间隔就是最佳时间间隔。
7.根据权利要求1所述的B扫描分割法抑制样品整体抖动的方法，其特征在于，步骤3具体包括以下步骤：
通过在快速扫描方向上的平均对数标度强度差标准差图像来定量评估抑制样品整体抖动的性能，选取不同时间间隔的扫描方案，计算得到在不同时间间隔下的平均对数标度强度差标准差值曲线图，平均对数标度强度差标准差值整体越低的扫描方案表明该方案的样品整体抖动抑制效果越好。
CN202210501036.2A 2022-05-09 2022-05-09 一种b扫描分割法抑制样品整体抖动的方法 Active CN114909987B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210501036.2A CN114909987B (zh) 2022-05-09 2022-05-09 一种b扫描分割法抑制样品整体抖动的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210501036.2A CN114909987B (zh) 2022-05-09 2022-05-09 一种b扫描分割法抑制样品整体抖动的方法

Publications (2)

Publication Number Publication Date
CN114909987A CN114909987A (zh) 2022-08-16
CN114909987B true CN114909987B (zh) 2024-02-02

Family

ID=82766767

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210501036.2A Active CN114909987B (zh) 2022-05-09 2022-05-09 一种b扫描分割法抑制样品整体抖动的方法

Country Status (1)

Country Link
CN (1) CN114909987B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115496917B (zh) * 2022-11-01 2023-09-26 中南大学 一种GPR B-Scan图像中的多目标检测方法及装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104134199A (zh) * 2013-05-01 2014-11-05 佳能株式会社 处理偏振敏感光学相干断层成像的偏振数据的方法和设备
CN104523233A (zh) * 2014-12-29 2015-04-22 浙江大学 基于复数互相关的微血管光学造影及抖动补偿方法与系统
CN107452029A (zh) * 2017-07-31 2017-12-08 中国医学科学院生物医学工程研究所 一种光学微血管血流成像方法
CN110063714A (zh) * 2018-01-22 2019-07-30 株式会社多美 光学相干层析成像装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015120055A1 (en) * 2014-02-04 2015-08-13 University Of Southern California Optical coherence tomography (oct) system with phase-sensitive b-scan registration

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104134199A (zh) * 2013-05-01 2014-11-05 佳能株式会社 处理偏振敏感光学相干断层成像的偏振数据的方法和设备
CN104523233A (zh) * 2014-12-29 2015-04-22 浙江大学 基于复数互相关的微血管光学造影及抖动补偿方法与系统
CN107452029A (zh) * 2017-07-31 2017-12-08 中国医学科学院生物医学工程研究所 一种光学微血管血流成像方法
CN110063714A (zh) * 2018-01-22 2019-07-30 株式会社多美 光学相干层析成像装置

Also Published As

Publication number Publication date
CN114909987A (zh) 2022-08-16

Similar Documents

Publication Publication Date Title
CN109907731B (zh) 基于特征空间的光学相干层析的三维血流造影方法
US8480582B2 (en) Image processing apparatus and ultrasonic diagnosis apparatus
EP0880937B1 (en) Ultrasonic diagnostic imaging system with doppler assisted tracking of tissue motion
US8094893B2 (en) Segmentation tool for identifying flow regions in an image system
EP1884197B1 (en) Image diagnosing device
Suhling et al. Myocardial motion analysis from B-mode echocardiograms
US5533510A (en) Real time ultrasound endocardial displacement display
US7628755B2 (en) Apparatus and method for processing an ultrasound image
CN108670239B (zh) 一种基于特征空间的三维血流成像方法与系统
US20080095417A1 (en) Method for registering images of a sequence of images, particularly ultrasound diagnostic images
JP2014525328A (ja) 針を検出して追跡する方法
US8313437B1 (en) Vascular ultrasound intima-media thickness (IMT) measurement system
DE112004001861T5 (de) System und Verfahren zur Analyse örtlicher verformbarer Bewegungen
CN110111386B (zh) 一种肌肉组织中结构性角度的自动测量方法及系统
CN114909987B (zh) 一种b扫描分割法抑制样品整体抖动的方法
US12089989B2 (en) Analyzing apparatus and analyzing method
CN114748032A (zh) 一种基于oct血管成像技术的运动噪声补偿方法
CN107274405A (zh) 一种超声微血管血流成像方法
CN114565711A (zh) 基于深度学习的心脏图像重建方法及系统
Wolf et al. ROPES: A semiautomated segmentation method for accelerated analysis of three-dimensional echocardiographic data
EP1772825A1 (en) Method for registering images of a sequence of images, particularly ultrasound diagnostic images
CN112396622A (zh) 基于多维特征空间的微血流图像分割量化方法和系统
Sjoerdsma et al. A spatial near-field clutter reduction filter preserving tissue speckle in echocardiography
CN116138760A (zh) 一种自适应增强的激光散斑衬比血流成像方法
US7218767B2 (en) Method of improving the resolution of a medical nuclear image

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