CN108618799B - 一种基于空间相干性的超声ct成像方法 - Google Patents
一种基于空间相干性的超声ct成像方法 Download PDFInfo
- Publication number
- CN108618799B CN108618799B CN201810372759.0A CN201810372759A CN108618799B CN 108618799 B CN108618799 B CN 108618799B CN 201810372759 A CN201810372759 A CN 201810372759A CN 108618799 B CN108618799 B CN 108618799B
- Authority
- CN
- China
- Prior art keywords
- imaging
- ultrasonic
- spatial coherence
- signal
- array
- 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 45
- 238000013170 computed tomography imaging Methods 0.000 title claims abstract description 20
- 238000003384 imaging method Methods 0.000 claims abstract description 66
- 238000001914 filtration Methods 0.000 claims abstract description 31
- 238000004364 calculation method Methods 0.000 claims abstract description 5
- 238000007781 pre-processing Methods 0.000 claims abstract description 5
- 238000012545 processing Methods 0.000 claims abstract description 5
- 238000012805 post-processing Methods 0.000 claims abstract description 3
- XREKLQOUFWBSFH-UHFFFAOYSA-N dimethyl 2-acetylbutanedioate Chemical compound COC(=O)CC(C(C)=O)C(=O)OC XREKLQOUFWBSFH-UHFFFAOYSA-N 0.000 claims abstract 2
- 238000013507 mapping Methods 0.000 claims description 14
- 230000006835 compression Effects 0.000 claims description 13
- 238000007906 compression Methods 0.000 claims description 13
- 230000008569 process Effects 0.000 claims description 13
- 238000001514 detection method Methods 0.000 claims description 12
- 230000035945 sensitivity Effects 0.000 claims description 4
- 239000000284 extract Substances 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 abstract description 2
- 238000002604 ultrasonography Methods 0.000 description 8
- 230000003902 lesion Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 210000001519 tissue Anatomy 0.000 description 6
- 206010006187 Breast cancer Diseases 0.000 description 3
- 208000026310 Breast neoplasm Diseases 0.000 description 3
- 238000003325 tomography Methods 0.000 description 3
- 210000000481 breast Anatomy 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000003111 delayed effect Effects 0.000 description 2
- 238000013399 early diagnosis Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 210000003484 anatomy Anatomy 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000001125 extrusion Methods 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 210000005075 mammary gland Anatomy 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000012502 risk assessment Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 238000012285 ultrasound imaging Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/48—Diagnostic techniques
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5229—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
- A61B6/5247—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from an ionising-radiation diagnostic technique and a non-ionising radiation diagnostic technique, e.g. X-ray and ultrasound
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/13—Tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
- A61B8/5238—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image
- A61B8/5261—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image combining images from different diagnostic modalities, e.g. ultrasound and X-ray
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Physics & Mathematics (AREA)
- Heart & Thoracic Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Optics & Photonics (AREA)
- High Energy & Nuclear Physics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Pulmonology (AREA)
- Theoretical Computer Science (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明公开了一种基于空间相干性的超声CT成像方法，包括以下步骤：(1)采集数据，获得原始回波数据；(2)数据预处理；(3)成像区域网格化；(4)针对每个网格成像点的基于空间相干性的DMAS计算处理，得到网格成像点的初始成像信号；(5)数据后处理，最终得到超声CT图像。本发明基于滤波延时相乘叠加算法中虚拟接收信号形成的特点，通过计算组成滤波延时相乘叠加算法中虚拟接收信号的任意两个信号的空间相干性，对虚拟接收信号进行加权，与现有技术相比能够有效解决零相位滤波延时相乘叠加算法在超声CT反射成像领域应用时存在的问题，能够重建出高对比度，高对比度噪声比和低旁瓣水平的超声CT图像。
Description
技术领域
本发明属于超声断层成像技术领域，更具体地，涉及一种基于空间相干性的超声CT成像方法。
背景技术
超声CT成像，即超声断层成像，具有高分辨率、360°全方位扫描、无挤压、无辐射和低成本的优势，在乳腺癌筛查和早期诊断方面具有重要的临床价值和应用前景。
超声CT成像方式包括反射成像、声速成像和衰减成像。反射成像刻画了目标的解剖结构，声速和衰减成像显示了目标的功能变化，三者相互补充，不可或缺。
由美国Karmanos癌症中心研制的计算超声风险评估系统(Computed ultrasoundrisk evaluation，CURE)可以重建出亚毫米级高分辨率的超声CT反射图像，目前在乳腺癌早期检测领域处于领先地位。这种扫描器每次使用单阵元向环形阵列的中心发射超声波，所有阵元接收，直到每个阵元循环发射一次，因此可以采集到多个角度的超声信号，重建出斑点噪声少的高分辨率图像。CURE系统采用了延时叠加(Delay and sum)算法，然而这种方法重建的图像在对比度，噪声抑制等方面仍有待改进。零相位滤波(Zero-phase filtereddelay multiply and sum，ZPF-DMAS)延时相乘叠加算法近年来被引入到超声CT领域，尽管该方法能够一定程度上抑制旁瓣和噪声，还是发现该方法重建图像时仍存在以下问题：(1)旁瓣水平和对比度仍需改进；(2)背景不均匀，对比度噪声比(Contrast to noise ratio，CNR)较低。
发明内容
针对现有技术的以上缺陷或改进需求，本发明的目的在于提供一种基于空间相干性的超声CT成像方法，基于滤波延时相乘叠加算法中虚拟接收信号形成的特点，通过计算组成滤波延时相乘叠加算法中虚拟接收信号的两个信号的空间相干性，对虚拟接收信号进行加权，形成一种新的基于空间相干性的滤波延时相乘叠加波束形成算法，与现有技术相比能够有效解决零相位滤波延时相乘叠加算法在超声CT反射成像领域应用时存在的问题，能够重建出高对比度，高对比度噪声比和低旁瓣水平的超声CT图像。
为实现上述目的，按照本发明，提供了一种基于空间相干性的超声CT成像方法，其特征在于，包括以下步骤：
(1)采集数据：
对环形阵列中的各个阵元按顺时针方向或逆时针方向由1到N开始编号，所述环形阵列中的各个阵元均匀地分布在圆环上，所述N为该环形阵列中的阵元总数；接着，由编号为1的阵元开始发射超声信号直到编号为N的阵元；每个所述阵元发射超声信号时，该环形阵列中的各个阵元均接收并采集超声信号，获得原始回波数据；
(2)数据预处理：
对所述步骤(1)得到的所述原始回波数据进行滤波处理，滤除其中的一些干扰信号，得到一次滤波后的信号；
(3)成像区域网格化：将成像区域按照直角坐标或极坐标进行网格化，使成像区域分为G个网格成像点；
(4)针对每个所述网格成像点的基于空间相干性的DMAS计算处理，具体是针对每一个所述网格成像点：
基于ZPF-DMAS方法，记深度为n处，由第j号阵元和第(j+m)号阵元这两个阵元接收到的一次滤波后的信号分别为sj(n)和sj+m(n)，则由sj(n)和sj+m(n)组成虚拟接收信号为
其中，m为正整数，具体为这两个阵元中较大的阵元编号与较小的阵元编号之间的差值；
基于空间相干性计算权重因子
其中，N为接收阵元总数，深度n1和深度n2为预先设定的用于计算空间相干性的一段信号的起始和终止位置，参数a为预先设定的用于调整空间相干性的敏感度的参数；
基于空间相干性的权重因子加权后深度n处的信号为
该ySC[n]即为所述网格成像点的初始成像信号；
(5)数据后处理：对各个所述网格成像点的初始成像信号进行零相位滤波处理得到滤波成像信号，然后再依次进行包络检测、对数压缩和灰度映射，最终即可得到超声CT图像。
作为本发明的进一步优选，所述步骤(3)中，所述网格化是按照直角坐标进行的，具体是将所述成像区域分成p行q列，所述G满足G＝(p-1)×(q-1)，每个网格的中心取为成像点；
若被检测物为已知病变区域尺寸的病变组织时，任意一个成像网格的面积小于该病变组织内最小病变区域面积的1/4；
若被检测物为未知是否病变的组织时，任意一个成像网格的面积不超过(λ/2)×(λ/2)，其中λ为超声波波长。
作为本发明的进一步优选，所述步骤(3)中，所述网格化是按照极坐标进行的，具体是将所述成像区域以极坐标极点为中心均匀的划分出u条射线形成u个线段，再将每个线段均匀的划分成v个小线段，所述G满足G＝u×v+1；相邻两条所述射线之间的夹角不超过1°，任意一条所述小线段的长度不超过λ/2，其中λ为超声波波长。
作为本发明的进一步优选，所述步骤(4)中，对于所述预先设定的用于计算空间相干性的一段信号的起始和终止位置，这一段信号对应一个波长；所述参数a的取值范围为10～30。
作为本发明的进一步优选，所述步骤(2)中，所述滤波处理为带通滤波处理；优选的，该带通滤波处理是对所述原始回波数据进行6阶的频率在f1～f2范围内的巴特沃斯带通滤波处理；其中，f2不小于环形阵列换能器中心频率f0，
Q为滤波器的品质因数。
作为本发明的进一步优选，所述步骤(5)中，所述滤波处理是滤除直流信号和频率值超过环形阵列换能器中心频率值至少2倍的高频信号；优选的，滤除的所述高频信号其频率值超过环形阵列超声换能器中心频率值至少4倍。
作为本发明的进一步优选，所述步骤(5)中，所述包络检测，具体是对于各个成像点的所述滤波成像信号，检测其上包络，提取出其中的低频分量，即为被检测物信息数据；所述低频分量优选为频率为20KHz～2MHz的分量；
所述对数压缩，具体是对于获得的被检测物信息数据取对数进行压缩，控制压缩后的数据在40dB至60dB之间；
所述灰度映射，具体是对于获得的压缩后的数据，采用线性映射为0至255之间或0至511之间的灰度值，得到灰度值数据。
通过本发明所构思的以上技术方案，与现有技术相比(如ZPF-DMAS等)，通过计算组成滤波延时叠加算法中虚拟接收信号的两个信号的空间相干性，对虚拟接收信号进行加权，能够最终重建出高对比度，高对比度噪声比和低旁瓣水平的超声CT图像。
超声信号的空间协方差(Spatial covariance)表示了延时接收信号在空间中两点的相似性，是随接收信号的空间间隔变化的函数。它由Van Cittert-Zernike(VCZ)定理定义为发射波形和散射函数的内积的平方的傅里叶变换(可参见：Kakkad V,Dahl J,Ellestad S,et al.In vivo application of short-lag spatial coherence andharmonic spatial coherence imaging in fetal ultrasound[J].Ultrasonic Imaging,2015,37(2):101-116)。空间协方差函数除以归一化值(对应RF信号的方差)得到空间相干函数(Spatial coherence function)。由于来自散射点p处的回波信号在各向同性的均匀介质中传播时声学特性并未发生改变，因此经过延迟聚焦后信号一致性很好，空间相干性较高。其他没有散射点的成像区域对应的回波信号可看作是噪声，而噪声信号空间相干性较差，因此采用空间相干性对p点的波束加权可以有效抑制旁瓣和系统噪声。
本发明提出的波束形成方法从空间相干性角度对超声断层成像领域最新提出的ZPF-DMAS方法进行了优化。ZPF-DMAS方法中的虚拟接收信号的来源和阵元间的空间位置有一种潜在的关系，本发明率先发现并利用了这种关系。在本方法的实施过程中，由于高反射区域的回波信号加权空间的相干性后幅度会降低，导致重建图像中低回声区域和高回声区域的对比度被削弱，因此本发明优选提出了调整空间相干性敏感度的参数a，该参数可以放大高低回声区域的信号强度差异，提高重建图像中低回声区域和高回声区域的对比度。
本发明在采用基于空间相干性的零相位滤波延时相乘叠加(ZPF-SCDMAS)波束形成方法重建超声CT图像时，首先对环形阵列超声换能器采集到的数据进行预处理，其次将成像区域网格化，分成若干个成像点。在每个成像点处进行基于空间相干性的滤波延时相乘叠加波束形成，包括(a)计算每个成像点处组成滤波延时相乘叠加算法中虚拟接收信号的两个信号的空间相干性和(b)根据空间相干性加权虚拟接收信号得到每个成像点处的值。然后将每条线数据采用零相位滤波，保留2倍中心频率附近的频率分量。最后通过包络检测、对数压缩、灰度映射将图像显示出来。
本发明所提超声CT成像方法，根据两个阵元的空间相干性可以抑制噪声信号的干扰，提高轴上信号的权重，限制离轴信号的影响，从而能够极大地降低旁瓣和噪声水平，提高对比度(～7dB，即约7dB)和对比度噪声比(～5dB)，提高超声CT成像的质量，有助于临床医生对乳腺癌的早期诊断。
附图说明
图1是本发明基于空间相干性的超声CT成像方法(ZPF-SCDMAS)流程图。
图2中图(a)为传统DAS算法，图(b)为ZPF-DMAS算法，图(c)为ZPF-SCDMAS算法的重建的乳腺体模的超声CT图像(动态范围为60dB)。
图3中图(a)为传统DAS算法，图(b)为ZPF-DMAS算法，图(c)为ZPF-SCDMAS算法的重建的真实人体乳腺的超声CT图像(动态范围为60dB)。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
实施例1
本发明中基于空间相干性的超声CT成像方法，其步骤包括数据采集、数据预处理、成像区域网格化、根据基于空间相干性(SC)的延时相乘叠加(DMAS)方法计算每个成像点处的值、零相位带通滤波、图像显示步骤。
(1)数据采集:环形阵列的阵元均匀地分布在圆环上，当一个阵元发射超声信号时，所有阵元(包括发射该信号的发射阵元本身)，负责接收信号。每个阵元循环发射一次，方向可以按顺时针方向或逆时针方向。此外，可以按顺时针方向或逆时针方向对所有阵元进行编号，阵元编号依次为1、2、3、……N，N为阵元总数。
(2)数据预处理：对步骤1中采集到的信号进行滤波，滤除一些干扰信号。本系统可以使用带通滤波，具体来说，在进行乳腺体模成像前，由于环形阵列换能器的中心频率为2.5MHz，所以可以对数据进行了6阶的频率在0.5MHz-2.5MHz范围内的巴特沃斯带通滤波(当然也可以采用其他频率范围，如f1～f2，只要频率上限f2满足不小于超声换能器中心频率f0(该超声换能器中心频率即环形阵列换能器中心频率)，频率下限f1满足
其中Q为滤波器的品质因数)。排除一些高频的噪声信号干扰，提高信噪比。
(3)成像区域网格化：将成像区域网格化，例如可以分成p行q列，这样整个成像区域分为G个网格，G＝(p-1)*(q-1)。理论上剖分的网格尺寸可以无限接近0(即无限小)。虽然网格尺寸越小，求解出来的成像点处的值越精确，但是网格越密，计算量越大，所以最佳的G由成像区域和计算效率共同决定。例如，一般地，病变组织的尺寸较大时，剖分不要求那么密，网格数G可适当降低，只要满足网格尺寸小于最小病变组织尺寸的1/4即可；当未知成像区域中病变组织尺寸时，剖分的网格大小不超过(λ/2)*(λ/2)，λ为超声波波长。另外，成像区域也可以按照极坐标网格化(后面成像时需进行极坐标与直角坐标的转换)，例如可以均匀分成u条线，每条线上v个小线段，这样整个成像区域分为G个成像点，G＝u*v+1(加1即包括上极坐标原点)。例如，两条线之间的夹角一般不超过1°，每条线上两点间间隔不超过λ/2。
(4)每个成像点处基于空间相干性(SC)的DMAS计算：
在ZPF-DMAS方法中，深度为n处，由阵元j，j+m接收到的信号sj(n)，sj+m(n)组成虚拟接收信号为
其中，m为两阵元的阵元序号之间的间隔。用空间相干系数加权该虚拟接收信号，权重如下：
这里，N为接收阵元总数，深度n1和深度n2为所选的用于计算空间相干性的一段信号的起始和终止位置，起始和终止位置可以预先设定，例如，选取信号长度一般为一个波长。参数a为预先设定的用于调整空间相干性的敏感度的参数，取值范围可以为10～30的任意实数。因此，空间相干系数加权后深度n处的信号为
(5)零相位带通滤波：由于在步骤5中信号相乘会产生直流信号或高频信号，所以为了避免这些信号的干扰，将线数据经过零相位带通滤波，滤除直流信号和高频信号(高频信号例如可以是频率大于10MHz的信号)，例如可以保留2倍中心频率附近的频率分量。
(6)包络检测：包络检测检出零相位带通滤波后信号的上包络，提取出回波信号中携带的低频分量(低频分量例如可以是频率为20KHz～2MHz的分量)，即被测物体信息。本发明中使用的包络检测的方法是希尔伯特变换法，原始信号经过希尔伯特变换将得到原始信号的正交信号。以原始信号为实部，以希尔伯特变换得到的信号为虚部构造解析信号，这个解析信号的模就是所要求的实信号的包络。
(7)对数压缩：对数压缩将线数据信号归一化后，取以10为底的对数，并乘以20，单位是dB。取完对数之后可以调整回波的动态范围，以得到最好的成像效果，一般调整为40dB或60dB。具体地说，动态范围为40dB指的是将信号中的最大值映射到40dB，将比最大信号小40dB的信号以及更小的信号映射到0dB。
(8)灰度映射：灰度映射采用简单的线性映射(其他还有动态范围压缩法等)，即成比例地将最弱的信号映射到0，将最强的信号映射到255或511(分别对应图像显示中常用的256灰度级或512灰度级的图像)。
(9)图像显示：经过灰度映射之后的数据，直接调用MATLAB中的成像函数即可得到超声CT图像。
从图2和图3中我们可以看出本发明所提方法对超声CT图像的影响。本发明ZPF-SCDMAS方法重建出的超声CT图像具有高对比度，高对比度噪声比和低旁瓣水平。
当网格化是按照极坐标情况进行时，和上述按照直角坐标网格化情况不同的是：在步骤(6)包络检测之后，步骤(7)对数压缩前，通过插值将图像中每个像素点的值用极坐标系下的成像点的值来表示，从而实现了极坐标系向直角坐标系的转换，然后按照步骤(7)-(9)进行。
本发明中的数学符号(包括函数符号等)均满足数学领域的常规定义。除了上述实施例中所具体描述的包络检测、对数压缩和灰度映射处理过程外，也可采用现有技术中其他包络检测、对数压缩和灰度映射处理方式。
本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (8)
1.一种基于空间相干性的超声CT成像方法，其特征在于，包括以下步骤：
(1)采集数据：
对环形阵列中的各个阵元按顺时针方向或逆时针方向由1到N开始编号，所述环形阵列中的各个阵元均匀地分布在圆环上，所述N为该环形阵列中的阵元总数；接着，由编号为1的阵元开始发射超声信号直到编号为N的阵元；每个所述阵元发射超声信号时，该环形阵列中的各个阵元均接收并采集超声信号，获得原始回波数据；
(2)数据预处理：
对所述步骤(1)得到的所述原始回波数据进行滤波处理，滤除其中的一些干扰信号，得到一次滤波后的信号；
(3)成像区域网格化：将成像区域按照直角坐标或极坐标进行网格化，使成像区域分为G个网格成像点；
(4)针对每个所述网格成像点的基于空间相干性的DMAS计算处理，具体是针对每一个所述网格成像点：
基于ZPF-DMAS方法，记深度为n处，由第j号阵元和第(j+m)号阵元这两个阵元接收到的一次滤波后的信号分别为sj(n)和sj+m(n)，则由sj(n)和sj+m(n)组成虚拟接收信号为
其中，m为正整数，具体为这两个阵元中较大的阵元编号与较小的阵元编号之间的差值；
基于空间相干性计算权重因子
其中，N为接收阵元总数，深度n1和深度n2为预先设定的用于计算空间相干性的一段信号的起始和终止位置，参数a为预先设定的用于调整空间相干性的敏感度的参数；所述参数a的取值范围为10～30；
基于空间相干性的权重因子加权后深度n处的信号为
该ySC[n]即为所述网格成像点的初始成像信号；
(5)数据后处理：对各个所述网格成像点的初始成像信号进行零相位滤波处理得到滤波成像信号，然后再依次进行包络检测、对数压缩和灰度映射，最终即可得到超声CT图像。
2.如权利要求1所述基于空间相干性的超声CT成像方法，其特征在于，所述步骤(3)中，所述网格化是按照极坐标进行的，具体是将所述成像区域以极坐标极点为中心均匀的划分出u条射线形成u个线段，再将每个线段均匀的划分成v个小线段，所述G满足G＝u×v+1；相邻两条所述射线之间的夹角不超过1°，任意一条所述小线段的长度不超过λ/2，其中λ为超声波波长。
3.如权利要求1所述基于空间相干性的超声CT成像方法，其特征在于，所述步骤(4)中，对于所述预先设定的用于计算空间相干性的一段信号的起始和终止位置，这一段信号对应一个波长。
4.如权利要求1所述基于空间相干性的超声CT成像方法，其特征在于，所述步骤(2)中，所述滤波处理为带通滤波处理。
5.如权利要求4所述基于空间相干性的超声CT成像方法，其特征在于，所述带通滤波处理是对所述原始回波数据进行6阶的频率在f1～f2范围内的巴特沃斯带通滤波处理；其中，f2不小于环形阵列超声换能器中心频率f0，
Q为滤波器的品质因数。
6.如权利要求1所述基于空间相干性的超声CT成像方法，其特征在于，所述步骤(5)中，所述滤波处理是滤除直流信号和频率值超过环形阵列超声换能器中心频率值至少2倍的高频信号。
7.如权利要求6所述基于空间相干性的超声CT成像方法，其特征在于，滤除的所述高频信号其频率值超过环形阵列超声换能器中心频率值至少4倍。
8.如权利要求1所述基于空间相干性的超声CT成像方法，其特征在于，所述步骤(5)中，所述包络检测，具体是对于各个成像点的所述滤波成像信号，检测其上包络，提取出其中的低频分量，即为被检测物信息数据；所述低频分量为频率为20KHz～2MHz的分量；
所述对数压缩，具体是对于获得的被检测物信息数据取对数进行压缩，控制压缩后的数据在40dB至60dB之间；
所述灰度映射，具体是对于获得的压缩后的数据，采用线性映射为0至255之间或0至511之间的灰度值，得到灰度值数据。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810372759.0A CN108618799B (zh) | 2018-04-24 | 2018-04-24 | 一种基于空间相干性的超声ct成像方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810372759.0A CN108618799B (zh) | 2018-04-24 | 2018-04-24 | 一种基于空间相干性的超声ct成像方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108618799A CN108618799A (zh) | 2018-10-09 |
| CN108618799B true CN108618799B (zh) | 2020-06-02 |
Family
ID=63694254
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810372759.0A Active CN108618799B (zh) | 2018-04-24 | 2018-04-24 | 一种基于空间相干性的超声ct成像方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108618799B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110009032B (zh) * | 2019-03-29 | 2022-04-26 | 江西理工大学 | 一种基于高光谱成像的组装分类方法 |
| CN113647978B (zh) * | 2021-08-18 | 2023-11-21 | 重庆大学 | 一种带有截断因子的高鲁棒性符号相干系数超声成像方法 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102435992A (zh) * | 2011-09-26 | 2012-05-02 | 重庆博恩富克医疗设备有限公司 | 基于广义相干系数的合成聚焦成像方法 |
| CN102641136A (zh) * | 2011-02-21 | 2012-08-22 | 三星电子株式会社 | 超声波束成型的方法和用于该方法的设备 |
| CN102727256A (zh) * | 2012-07-23 | 2012-10-17 | 重庆博恩富克医疗设备有限公司 | 基于虚拟阵元的双聚焦波束合成方法及装置 |
| CN103536316A (zh) * | 2013-09-22 | 2014-01-29 | 华中科技大学 | 一种空时平滑相干因子类自适应超声成像方法 |
| CN105813572A (zh) * | 2013-12-09 | 2016-07-27 | 皇家飞利浦有限公司 | 基于图像信息的图像复合 |
| WO2016159395A1 (ko) * | 2015-03-27 | 2016-10-06 | 알피니언메디칼시스템 주식회사 | 공간 스무딩 연산이 간단한 빔포밍 장치, 초음파 이미징 장치 및 빔포밍 방법 |
| CN105997153A (zh) * | 2016-07-15 | 2016-10-12 | 华中科技大学 | 一种超声ct的成像方法 |
| CN108403148A (zh) * | 2018-04-16 | 2018-08-17 | 武汉维视医学影像有限公司 | 一种基于mv自适应波束形成的超声ct成像方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2011123529A1 (en) * | 2010-04-02 | 2011-10-06 | Duke University | Methods, systems and apparatuses for van-cittert zernike imaging |
| US10064602B2 (en) * | 2014-06-03 | 2018-09-04 | Siemens Medical Solutions Usa, Inc. | Coherence ultrasound imaging with broad transmit beams |
-
2018
- 2018-04-24 CN CN201810372759.0A patent/CN108618799B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102641136A (zh) * | 2011-02-21 | 2012-08-22 | 三星电子株式会社 | 超声波束成型的方法和用于该方法的设备 |
| CN102435992A (zh) * | 2011-09-26 | 2012-05-02 | 重庆博恩富克医疗设备有限公司 | 基于广义相干系数的合成聚焦成像方法 |
| CN102727256A (zh) * | 2012-07-23 | 2012-10-17 | 重庆博恩富克医疗设备有限公司 | 基于虚拟阵元的双聚焦波束合成方法及装置 |
| CN103536316A (zh) * | 2013-09-22 | 2014-01-29 | 华中科技大学 | 一种空时平滑相干因子类自适应超声成像方法 |
| CN105813572A (zh) * | 2013-12-09 | 2016-07-27 | 皇家飞利浦有限公司 | 基于图像信息的图像复合 |
| WO2016159395A1 (ko) * | 2015-03-27 | 2016-10-06 | 알피니언메디칼시스템 주식회사 | 공간 스무딩 연산이 간단한 빔포밍 장치, 초음파 이미징 장치 및 빔포밍 방법 |
| CN105997153A (zh) * | 2016-07-15 | 2016-10-12 | 华中科技大学 | 一种超声ct的成像方法 |
| CN108403148A (zh) * | 2018-04-16 | 2018-08-17 | 武汉维视医学影像有限公司 | 一种基于mv自适应波束形成的超声ct成像方法 |
Non-Patent Citations (1)
| Title |
|---|
| Spatial smoothing coherence factor for ultrasound computed tomography;Cuijuan lou等;《proceedings of SPIE》;20161231;第9790卷;979008-1至979008-8 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108618799A (zh) | 2018-10-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Hyun et al. | Beamforming and speckle reduction using neural networks | |
| US11030780B2 (en) | Ultrasound speckle reduction and image reconstruction using deep learning techniques | |
| KR101407425B1 (ko) | 최적 스펙클 추적을 위한 이중 경로 처리 | |
| US9081097B2 (en) | Component frame enhancement for spatial compounding in ultrasound imaging | |
| CN105411626B (zh) | 一种基于超声ct的合成孔径成像方法及系统 | |
| US20130343627A1 (en) | Suppression of reverberations and/or clutter in ultrasonic imaging systems | |
| US11408987B2 (en) | Ultrasonic imaging with multi-scale processing for grating lobe suppression | |
| CN108742705A (zh) | 一种实时检测肌肉形态参数的超声成像设备及方法 | |
| Morgan et al. | Multi-covariate imaging of sub-resolution targets | |
| Yoon et al. | Frequency equalized compounding for effective speckle reduction in medical ultrasound imaging | |
| Li et al. | Evaluational spatial compounding | |
| Gauss et al. | Wavefront estimation in the human breast | |
| CN108618799B (zh) | 一种基于空间相干性的超声ct成像方法 | |
| Wang et al. | Performance evaluation of coherence-based adaptive imaging using clinical breast data | |
| KR20150112057A (ko) | 공간 일관성 기초 초음파 신호 처리 모듈 및 그에 의한 초음파 신호 처리 방법 | |
| CN110974293B (zh) | 一种基于c型探头的合成孔径成像方法 | |
| AU2004264847B2 (en) | Method and apparatus for ultrasonic imaging | |
| TWI550268B (zh) | Method of Improving Sensitivity of Quantitative Tissue Characteristic of Ultrasonic | |
| CA2587137A1 (en) | Imaging and measuring bladder wall thickness | |
| Lou et al. | Zero-phase filtered delay multiply and sum in ultrasound computed tomography | |
| Wang et al. | Plane-wave ultrasound imaging based on compressive sensing with low memory occupation | |
| US20220022848A1 (en) | Ultrasound imaging system using coherence estimation of a beamformed signal | |
| Wang et al. | Synthetic aperture focusing technique for 3-D ultrasound computed tomography | |
| Lou et al. | A fast contrast improved zero-phase filtered delay multiply and sum in ultrasound computed tomography | |
| Wang et al. | Multiline acquisition beamforming for ultrasound computed tomography |
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 | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20220817 Address after: Floor 1-5, Building A, Building A, Phase I, Donghu High-tech Creative City, Wutong Lake New District, Ezhou City, Hubei Province, 431000 Patentee after: WUHAN WEISHI MEDICAL SCIENCE IMAGE Co.,Ltd. Address before: 430074 Hubei Province, Wuhan city Hongshan District Luoyu Road No. 1037 Patentee before: HUAZHONG University OF SCIENCE AND TECHNOLOGY |
|
| CP03 | Change of name, title or address | ||
| CP03 | Change of name, title or address |
Address after: Room 02, 18th Floor, Building 2, New Factory, No. 438, Gaoxin Avenue, Wuhan East Lake New Technology Development Zone, Wuhan, Hubei 430000, China Patentee after: Weishi Medical Imaging Co.,Ltd. Address before: Floor 1-5, Building A, Building A, Phase I, Donghu High-tech Creative City, Wutong Lake New District, Ezhou City, Hubei Province, 431000 Patentee before: WUHAN WEISHI MEDICAL SCIENCE IMAGE Co.,Ltd. |