一种超声图像清晰度提升方法技术

本发明专利技术公开了一种超声图像清晰度提升方法，设计了一种基于超声成像中散射与反射模型，结合自然图像降质模型中提高对比度方法求解，恢复强边缘或高密度区的反射信号，对于不同亮度区域，选择不同等级曲线控制散射影响，从而获取更高对比度的清晰图像。本发明专利技术提供一种具有更高清晰度的医学超声图像显示方法，尤其对于血管，胎儿羊水的等液性物质与组织的分界处，加强边界轮廓及增强组织内部小细节的识别，通过本发明专利技术方法处理提供高分辨率的图像，从而更好地辅助医生诊断。

【技术实现步骤摘要】

本专利技术属于医学超声成像显示
，涉及，具 体地说，涉及一种医学超声图像质量的对比度与细节分辨率提升的方法。
技术介绍
医学超声影像因其对人体无损伤、非侵入、成像方式多种多样、检查费用低、实时 性好等优点，已广泛应用于医学诊断的各个领域；特别是二维B模式影像，可提供人体组织 结构在不同病理状态下的切面断层信息，能较为准确的发现人体组织中较深处微小病变， 已成为医学检查的首选项目。超声成像是基于超声波的成像技术，以声学为物理基础。受 声学物理特性的影响，如散射及干涉、组织声衰减及不同传播速度，加之不可避免的系统电 子噪声等多种因素共同作用，超声影像相比与自然图像表现出较差的对比分辨率、细节分 辨率及亮度均匀一致性，给超声医师的临床诊断工作带来困难。为了更好地提供清晰图像，对比度提升，图像增强等后处理操作是必须的，按处理 方法可以分为线性滤波和非线性滤波，按处理域可以分为空域方法和频域方法。空域方法 是直接在图像上对象元进行处理；频域方法的基础是卷积定理。对比度小的图像，表现为灰 度动态范围小，对某一段范围的灰度进行拉伸变换，增加灰度动态范围。图像的灰度动态范 围很大，而显示设备的显示能力不够，常用的方法是对数压缩，调整图像的动态范围，或为 突出某一灰度级采用灰度级切片。直方图均衡化方法单纯从像素上提升，寻找一个灰度变 换，使得变换后图像在每一灰度级的像素个数都相同不能保证组织的真实性。自然图像对 比度增强中，采用一系列技术去改善图像的视觉效果，或将图像转换成一种更适合于人或 机器进行分析处理的形式。结合图像对比度提升，与超声医学影像成像特点，我们知道，医学超声成像中探头 接收的回波信号主要由组织对探头所发射的超声波反射与散射信号等所形成，利用接收到 的信号形成图像，从而判别组织结构及其中可能的病变。复杂的医学超声图像被各种模型 与算法所描述，但由于组织本来的复杂性，噪声的影响，没有一种方法能完全准确的描述身 体内组织情况，我们只能无限去接近它，噪声对图像的影响就好比给图像蒙上一层纱布，从 而使图像对比灵敏度下降。医学超声图像的显示应用中，现有的医学超声诊断仪设备的后 处理中，常常加入时间复合，空间复合，频率复合，平滑，声速衰减补偿等相关后处理，在提 高显示图像信噪比的同时，会使得超声图像蒙上一层更重的纱布似的“雾”，从而影响对比 度与细节分辨率。噪声和边缘都表现为图像上灰度变化剧烈的区域。平滑噪声的时候也模 糊了边缘，增强边缘的同时也强化了噪声，这是一个基本的矛盾。为提高显示对比度，非线 性对比拉升技术、对数函数,指数函数,伽马函数，Sigmoid函数(S形状)等自然图像对比 度提高的处理方法常用于优化医学超声图像，或是一些从人眼视觉感知上提升图像对比度 的方法，总体而言，这些方法仅是从像素上去提高图像质量，几乎没涉及到医学超声成像原 理。当目标物远小于波长时，微小结构中散射起主要作用，生物组织内的散射主要由组织结构的大小和声学特性决定。在医学超声图像的显示中，通常自适应调节最强的边界或高强度值接近饱和值(例如，在8位数字表达中，使用255值)，这样能获得更丰富的对比度且便于医生区分组织中病变。
技术实现思路
本专利技术的目的在于克服上述技术缺陷，提供，该方法基于超声成像中散射与反射模型，借助自然图像中降质模型中提高对比度方法进行求解，恢复强边缘或高密度区的反射信号，对于不同亮度区域，选择不同等级曲线降低散射影响，从而获取更高对比度的清晰图像，提高诊断准确性。本专利技术尤其对于血管，胎儿羊水的等液性物质与组织的分界处，加强边界轮廓及增强组织内部小细节的识别，通过本专利技术方法处理提供高分辨率的图像，从而更好地辅助医生诊断。其技术方案为，包括如下步骤I)根据当前图像性能分析，进行对比度判定，是否需要清晰度提升处理，其中包括计算主要图像区域的均值与边缘数计算，对比度判定包括以下步骤a)计算出主要组织部位区域的像素均值；b)统计表明一幅图像如果有更强的对比度，将会产生更多的边缘数，计算出组织部位区域的边缘数；c)综合均值与边缘数的信息得出对比度判定2)根据超声图像形成特性，设定窗口的大小；根据采集当前图像系统设定的深度，超声频率，选择合理窗口大小Scx = O. 87 λ z/DScz = O. 9IC0/Af式中，S。,，Sra，分别为横向单元大小与纵向单元的大小，λ是波长，(λ = C0/Fc), z 为采集图像深度，D是超声波束形成时所激活的探头孔径，C0是声速，(例如，1. 54mm/ μ s)， Fc是中心频率，Δ f,通过FeXbandwidth,“bandwidth”是带宽，以探头发射频率为基础的一个分数(例如80% );3)按照设定窗口大小对整个图像进行遍历，在窗口内取最小值，获取块区域内暗元素；4)形成块区域内暗元素图，便于查看特性，并进行优化；5)通过块区域内暗元素图，计算出散射与强反射分界的A值；6)构建与区域像素亮度对应的减函数等级，由此得到ω ；7)通过前面步骤2)，3)，5)，6)操作，利用公式本文档来自技高网...

【技术保护点】
一种超声图像清晰度提升方法，其特征在于，包括如下步骤：1)根据当前图像性能分析，进行对比度判定，是否需要清晰度提升处理，其中包括计算主要图像区域的均值与边缘数计算，对比度判定包括以下步骤：a)计算出主要组织部位区域的像素均值；b)统计表明一幅图像如果有更强的对比度，将会产生更多的边缘数，计算出组织部位区域的边缘数；c)综合均值与边缘数的信息得出对比度判定：2)根据超声图像形成特性，设定窗口的大小；根据采集当前图像系统设定的深度，超声频率，选择合理窗口大小：Scx＝0.87λz/DScz＝0.91C0/Δf式中，Scx，Scz，分别为横向单元大小与纵向单元的大小，λ是波长，(λ＝C0/Fc)，z为采集图像深度，D是超声波束形成时所激活的探头孔径，C0是声速，Fc是中心频率，Δf，通过Fc×bandwidth，“bandwidth”是带宽，以探头发射频率为基础的一个分数；3)按照设定窗口大小对整个图像进行遍历，在窗口内取最小值，获取块区域内暗元素；4)形成块区域内暗元素图，便于查看特性，并进行优化；5)通过块区域内暗元素图，计算出散射与强反射分界的A值；6)构建与区域像素亮度对应的减函数等级，由此得到ω；7)通过前面步骤2)，3)，5)，6)操作，利用公式t(z)=1-ωminy∈Ω(z)(R(y)A)得出传输率；式中，t(z)，是强反射或高密度区信号到达探头处的信号强度，R(y)为处理前的图，ω为根据区域亮度调整的减函数，A为散射与强反射或高密度区的分界值，Ω为深度z处的步骤2)所得的一个块大小区域；8)获取传输率图，并进行适当优化；9)利用公式S(z)=R(z)-Amax(t0,t(z))+A获得按不同等级减少散射影响的强反射与高密度区图像，从而提升医学影像图的对比度与细节分辨率，其中，max为求取较大值，R(z)为探头接收信号强度图，S(z)为提升对比度后恢复出来的图，t0为一个下限值。...

【技术特征摘要】
1.一种超声图像清晰度提升方法，其特征在于，包括如下步骤 1)根据当前图像性能分析，进行对比度判定，是否需要清晰度提升处理，其中包括计算主要图像区域的均值与边缘数计算，对比度判定包括以下步骤 a)计算出主要组织部位区域的像素均值； b)统计表明一幅图像如果有更强的对比度，将会产生更多的边缘数，计算出组织部位区域的边缘数； c)综合均值与边缘数的信息得出对比度判定 2)根据超声图像形成特性，设定窗口的大小； 根据采集当前图像系统设定的深度，超声频率，选择合理窗口大小Scx = O. 87 λ z/DScz = 0.91C0/Af 式中，Scx, Sra，分别为横向单元大小与纵向单元的大小，λ是波长，(λ = C0/Fc), z为采集图像深度，D是超声波束形...

【专利技术属性】
技术研发人员：郭德全，佘跃，王立，刘东权，
申请(专利权)人：声泰特成都科技有限公司，
类型：发明
国别省市：