一种B型图像增益自动优化方法及装置制造方法及图纸

本发明专利技术提供一种B型图像增益自动优化方法及装置，包括如下步骤：S1：判断图像噪声区域及非噪声区域，S2：获取图像组织区域，S3：图像增益补偿曲线计算，S4：图像增益优化，提高了鲁棒性，本发明专利技术原理简单，计算复杂度低，能够针对不同个体不同部位自适应的计算图像的纵向和横向增益补偿。

An Automatic Optimizing Method and Device for B-Type Image Gain

【技术实现步骤摘要】
一种B型图像增益自动优化方法及装置
本专利技术是一种B型图像增益自动优化方法及装置，属于医疗领域。
技术介绍
现有技术中，随着医学技术和医疗诊断手段的不断发展，超声成像技术在临床诊断和医学科研等领域得到了广泛的应用。常见的超声成像B型系统处理流程包括发射控制模块按照指令激励探头，电信号转化为声信号。声波在人体组织中传播并发射，探头接收经人体组织反射的超声回波信号，将声信号转化为电信号。在模拟信号处理模块经过放大、滤波、时间增益补偿TGC(timegaincompensation,TGC)及模拟数字转换(ADC)等将模拟回波信号转换为数字回波信号。由于声波在发射、接收过程中能量会产生衰减，如果直接采用此回波信号进行处理，在不同深度图像的亮度表现并不一致，给用户诊断带来困难，因此在前端通常会进行时间增益补偿，由于是在模拟端处理也叫模拟时间增益补偿(analogtimegaincompensation,ATGC)，波束合成器将接收到的多路AD信号经过延时累加聚焦等处理得到一条RF信号。在获取整幅图像RF数据之后，数据经过信号处理模块包括IQ解调、低通滤波降采样、对数压缩，此时的对数压缩数据在深度方向亮度不均匀，不符合人眼观察习惯，需要用户手动调节电位器，调整不同深度下的增益，保证图像亮度一致。此时处理的信号是数字信号，因此也叫数字时间增益补偿(digitaltimegaincompensation,DTGC),同时需要手动调节旋钮即整体增益(globalgain)，调整合适的图像亮度。经过增益调整后的图像被送入图像处理模块，该模块包括动态范围调整、图像增强、数字扫描变换等最终送到显示器进行显示。在超声成像系统中增益调整对于临床诊断具有重要的意义，没有经过增益调整的图像亮度在超声设备上表现不均匀，目前超声厂家在出厂时已预设了增益，但用户在对患者进行疾病检查时，由于超声波的衰减特性，不同人体的声阻抗也有差异，厂家预设的增益并不能满足所有需求，这时需要用户手动调节DTGC、globalgain才能达到好的结果，给用户使用带来了负担，有时也无法调整出较好的结果。针对上述问题诞生了一种增益优化方法，该方法通过分析图像数据，计算出一条轴向增益补偿曲线，通过对轴向增益进行补偿，使图像亮度表现均匀。
技术实现思路
针对现有技术存在的不足，本专利技术目的是提供一种B型图像增益自动优化方法及装置，以解决上述
技术介绍
中提出的问题，本专利技术采用结构张量中特征值的判断方法，提高了鲁棒性。为了实现上述目的，本专利技术是通过如下的技术方案来实现：一种B型图像增益自动优化方法，包括如下步骤：S1：判断图像噪声区域及非噪声区域：获取超声对数压缩图像I(M*N，M为图像深度，N为图像宽度)根据超声噪声数据判断图像噪声区域及非噪声区域；S2：获取图像组织区域：二维图像的结构张量是对称且半正定的二维矩阵，因此I有两个特征值λ1、λ2，分别代表图像像素点最大和最小特征值：当λ1≈λ2≈0时，表示像素点灰度值在此邻域范围内变化较小，像素点位于平坦区域；当λ1＞＞λ2≈0时，表示在某一个方向上灰度值变化强，像素点位于边缘；当λ1≥λ2＞0时，表示在指定方向上特征向量存在很强的变化，像素点位于角点；利用上述特征值的表现形式，构建一个可靠函数P，可以用来判断超声图像中的组织区域；具体方法如下：a)为加快计算，分别对I进行横向、纵向平滑同时进行降采样处理，具体方法为对图像I进行二维高斯滤波，同时抽取奇数行、奇数列，获取图像I′；b)图像I′结构张量S(I)计算其中*代表卷积，δ为方差，取0.5，首先计算图像I′横向和纵向梯度、Ix、Iy，利用Ix、Iy计算结构张量，然后对结构张量做二维高斯平滑处理，为了加快计算，二维高斯平滑处理可以转化为一个一维纵向高斯平滑及一维横向高斯平滑；c)图像I′特征值计算d)图像组织区域TisReg判断根据上述特征值的表现形式，构建函数P：或P∈(0，1)，λ1、λ2属于非噪声区域在本专利技术中P取第二种方式，λ1、λ2属于非噪声区域：当λ1≈λ2≈0时，P接近于0，像素点位于平坦区域，认为是组织区域；当λ1＞＞λ2≈0或λ1≥λ2＞0时，认为是非组织区域；在此处提供一个判断组织区域的阈值TisThr，当P<TisThr时，认为是组织区域，此时TisReg为1，其余为非组织区域，TisReg为0；S3：图像增益补偿曲线计算判断是否参与轴向平均信号AxisMean、横向平均信号LateralMean计算，即：TisPreThr可以根据实际情况取值，在此方案中取0.9，当L取N时，计算轴向平均信号AxisMean，当L取M时，计算横向平均信号LateralMean，对0点采用相邻平均信号插值处理；最终轴向增益补偿曲线GainAxis为：GainAxis＝DstGain-xisMean，DstGain为期望增益，由于在纵向调整过程中，每一行进行了GainAxis的增益补偿，所以在此基础上要达到调整后亮度为DstGain的目标，横向增益补偿曲线GainLateral为：GainLateral＝LateralMean-∑LateralMean对GainAxis和GainLateral进行上采样处理，得到最终的增益补偿曲线；S4：图像增益优化根据公式Iout＝I+GainAxis+GainLateral，对图像进行横向及纵向增益优化，Iout为增益优化后图像。进一步地：应区分图像噪声区域及非噪声区域，具体通过关闭发射获取噪声数据，在横向取平均值获取一条纵向噪声曲线1-M，重复m次取其平均值，为减少系统的响应时间，噪声曲线预存在文件中，每次优化时直接使用噪声曲线数据处理，当对数压缩数据小于噪声曲线数据时，认为此处为非噪声区域，否则为噪声区域。一种B型图像增益装置，包括探头、发射控制模块、接收控制模块、模拟信号处理模块、波束合成器、信号处理模块、B型图像增益自动优化模块、图像处理模块、显示模块，具体来讲发射控制模块按照指令激励探头，电信号转为声信号，探头接收经组织反射的超声回波信号，将声信号转化为电信号，在模拟信号处理模块经过放大、滤波、时间增益补偿及模拟数字转换(ADC)等将模拟回波信号转换为数字回波信号；波束合成器将接收到的多路AD信号经过延时累加聚焦等处理得到一条RF信号；在获取RF数据之后，系统通过信号处理模块包括IQ解调、低通滤波降采样、对数压缩最后将数据送入B型图像增益自动优化模块，自动优化模块通过在键盘上设置一个快捷按钮进行启动，当启动该装置时通过该装置计算出合适的增益调节曲线，对B型图像进行增益调节，当关闭改按钮时则用系统默认的增益曲线进行调节；经过增益调整后的图像被送入图像处理模块，该模块包括动态范围调整、图像增强、数字扫描变换等最终送到显示器进行显示。B型图像增益自动优化模块，包括数据存储模块、图像增益曲线计算模块，图像增益曲线输出模块、图像增益优化处理模块：数据存储模块保存图像增益曲线计算模块需要用到的参数，包括噪声曲线数据及组织区域判断阈值，为了减少系统的响应时间，方便计算，噪声数据在图像增益优化前已经计算完成并保存在数据存储模块，在图像增益计算模块中直接使用，数据存储模块中的值可以根据实际情况进行修改；图像增益曲线计算模块本文档来自技高网...

【技术保护点】
1.一种B型图像增益自动优化方法，其特征在于包括如下步骤：S1：判断图像噪声区域及非噪声区域：获取超声对数压缩图像I(M*N，M为图像深度，N为图像宽度)根据超声噪声数据判断图像噪声区域及非噪声区域；S2：获取图像组织区域：二维图像的结构张量是对称且半正定的二维矩阵，因此I有两个特征值λ1、λ2，分别代表图像像素点最大和最小特征值：当λ1≈λ2≈0时，表示像素点灰度值在此邻域范围内变化较小，像素点位于平坦区域；当λ1＞＞λ2≈0时，表示在某一个方向上灰度值变化强，像素点位于边缘；当λ1≥λ2＞0时，表示在指定方向上特征向量存在很强的变化，像素点位于角点；利用上述特征值的表现形式，构建一个可靠函数P，可以用来判断超声图像中的组织区域；具体方法如下：a)为加快计算，分别对I进行横向、纵向平滑同时进行降采样处理，具体方法为对图像I进行二维高斯滤波，同时抽取奇数行、奇数列，获取图像I′；b)图像I′结构张量S(I)计算

【技术特征摘要】
1.一种B型图像增益自动优化方法，其特征在于包括如下步骤：S1：判断图像噪声区域及非噪声区域：获取超声对数压缩图像I(M*N，M为图像深度，N为图像宽度)根据超声噪声数据判断图像噪声区域及非噪声区域；S2：获取图像组织区域：二维图像的结构张量是对称且半正定的二维矩阵，因此I有两个特征值λ1、λ2，分别代表图像像素点最大和最小特征值：当λ1≈λ2≈0时，表示像素点灰度值在此邻域范围内变化较小，像素点位于平坦区域；当λ1＞＞λ2≈0时，表示在某一个方向上灰度值变化强，像素点位于边缘；当λ1≥λ2＞0时，表示在指定方向上特征向量存在很强的变化，像素点位于角点；利用上述特征值的表现形式，构建一个可靠函数P，可以用来判断超声图像中的组织区域；具体方法如下：a)为加快计算，分别对I进行横向、纵向平滑同时进行降采样处理，具体方法为对图像I进行二维高斯滤波，同时抽取奇数行、奇数列，获取图像I′；b)图像I′结构张量S(I)计算其中*代表卷积，δ为方差，取0.5，首先计算图像I′横向和纵向梯度、Ix、Iy，利用Ix、Iy计算结构张量，然后对结构张量做二维高斯平滑处理，为了加快计算，二维高斯平滑处理可以转化为一个一维纵向高斯平滑及一维横向高斯平滑；c)图像I′特征值计算d)图像组织区域TisReg判断根据上述特征值的表现形式，构建函数P：或λ1、λ2属于非噪声区域在本发明中P取第二种方式，λ1、λ2属于非噪声区域：当λ1≈λ2≈0时，P接近于0，像素点位于平坦区域，认为是组织区域；当λ1＞＞λ2≈0或λ1≥λ2＞0时，认为是非组织区域；在此处提供一个判断组织区域的阈值TisThr，当P<TisThr时，认为是组织区域，此时TisReg为1，其余为非组织区域，TisReg为0；S3：图像增益补偿曲线计算判断是否参与轴向平均信号AxisMean、横向平均信号LateralMean计算，即：TisPreThr可以根据实际情况取值，在此方案中取0.9，当L取N时，计算轴向平均信号AxisMean，当L取M时，计算横向平均信号LateralMean，对0点采用相邻平均信号插值处理；最终轴向增益补偿曲线GainAxis为：GainAxis＝DstGain-xisMean，DstGain为期望增益；由于在纵向调整过程中，每一行进行了GainAxis的增益补偿，所以在此基础上要达到调整后亮度为DstGain的目标，横向增益补偿曲线GainLateral为：GainLateral＝LateralMean-∑LateralMean对GainAx...

【专利技术属性】
技术研发人员：孙瑞超，黄帅，邢锐桐，陈晶，
申请(专利权)人：深圳蓝韵医学影像有限公司，
类型：发明
国别省市：广东,44