CN113057676B - Ivus系统的图像降噪方法 - Google Patents
Ivus系统的图像降噪方法 Download PDFInfo
- Publication number
- CN113057676B CN113057676B CN202110309664.6A CN202110309664A CN113057676B CN 113057676 B CN113057676 B CN 113057676B CN 202110309664 A CN202110309664 A CN 202110309664A CN 113057676 B CN113057676 B CN 113057676B
- Authority
- CN
- China
- Prior art keywords
- pixel
- pixel point
- image
- gray value
- interference
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/12—Diagnosis using ultrasonic, sonic or infrasonic waves in body cavities or body tracts, e.g. by using catheters
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0891—Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of blood vessels
-
- 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/5269—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Vascular Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本公开描述一种IVUS系统的图像降噪方法,其包括:利用IVUS系统获取超声图像;选取第一像素区域,计算第一像素区域的平均灰度值;比较各个像素点的原始灰度值与第一像素区域的平均灰度值并获得比较结果,并根据比较结果标记出干扰像素点;选取第二像素区域,分别计算所选定的干扰像素点的原始灰度值、第二像素区域中各个像素点的原始灰度值的中位值与包括所选定的干扰像素点的第一像素区域的平均灰度值之间的相对偏差,若该中位值与该平均灰度值的相对偏差更小,则使用该中位值替换所选定的干扰像素点的原始灰度值。根据本公开,能够优化降噪处理过程并有效提高超声图像的图像质量。
Description
本申请是申请日为2019年10月18日、申请号为201910996383.5、发明名称为血管 内超声系统的图像降噪方法的专利申请的分案申请。
技术领域
本公开涉及一种IVUS系统的图像降噪方法。
背景技术
血管内超声(Intravascular ultrasound,IVUS)系统通过超声探头采集患者血管内待检区域的超声图像,从而辅助医生对待检区域是否存在病变以及存在何种病变进行诊断和治疗。具体而言,血管内超声系统具有可以发射超声声束的超声探头,通过使用超声探头例如在人体血管内发射超声声束以获取用于显示人体血管的组织结构和几何形态的超声图像。
血管内超声系统在工作过程中,容易受到手术室内其它设备例如呼吸机等的强电磁脉冲干扰,从而导致在超声图像上产生噪声,影响医生阅图诊断。为减少这类噪声对诊断的不利影响,需要对超声图像进行降噪处理以提高其图像质量。在现有技术中,通常是选取待处理图像和若干张与待处理图像相邻的图像,并基于该相邻的图像对待处理图像中的各个像素点进行中值滤波来进行降噪处理。
然而,上述现有技术对待处理图像的处理过程需要以多张图像作为处理对象,增加了系统的复杂性;并且采用了无差别的中值滤波方式,该滤波方式难以有效地对待处理图像中受到干扰的像素点进行降噪处理,从而不能有效地提高超声图像的图像质量。
发明内容
本公开是有鉴于上述现有技术的状况而提出的,其目的在于提供一种能够优化降噪处理过程并且有效地对超声图像中的干扰像素点进行降噪的血管内超声系统的图像降噪方法。
为此,本公开第一方面提供了一种血管内超声系统的图像降噪方法,其特征在于,利用血管内超声系统的超声探头在旋转的情况下沿着旋转方向和血管的长度方向运动来获取来自于血管内的超声图像;在所述超声图像中遍历各个像素点,并获取各个像素点的原始灰度值;选取包括该像素点的第一像素区域,计算所述第一像素区域的各个像素点的原始灰度值的平均值作为所述第一像素区域的平均灰度值;比较所述各个像素点的原始灰度值与包括该像素点的所述第一像素区域的平均灰度值获得比较结果,并根据所述比较结果标记出干扰像素点;并且在所述超声图像中遍历各个干扰像素点,并选取包括该干扰像素点的第二像素区域,计算所述第二像素区域中各个像素点的原始灰度值的中位值,并使用该中位值替换该干扰像素点的原始灰度值,其中,所述第二像素区域包括含有该干扰像素点且沿着旋转方向的规定像素行、以及含有该干扰像素点且沿着血管的径向方向的规定像素列。
在本公开的第一方面中,通过比较超声图像中各个像素点的原始灰度值与包括该像素点的第一像素区域的平均灰度值,能够准确地识别出超声图像中的干扰像素点,并且能够有针对性地对所识别的干扰像素点进行中值滤波处理以降低因干扰而造成的噪声,从而能够优化降噪处理并且有效地对超声图像中的干扰像素点进行降噪。
在本公开第一方面所涉及的图像降噪方法中,可选地,通过计算所述各个像素点的原始灰度值与包括该像素点的所述第一像素区域的平均灰度值的比值,来获得所述各个像素点的原始灰度值与包括该像素点的所述第一像素区域的平均灰度值的比较结果。由此,能够准确地识别超声图像中的干扰像素点。
在本公开第一方面所涉及的图像降噪方法中,可选地,当所述比值大于第一预定阈值或小于第二预定阈值时,将该像素点标记为所述干扰像素点,其中,所述第一预定阈值大于所述第二预定阈值。在这种情况下,通过将上述比值设定在预定区间,从而能够更方便地识别超声图像中的干扰像素点。
在本公开第一方面所涉及的图像降噪方法中,可选地,在选取包括该像素点的第一像素区域时,以该像素点作为所述第一像素区域的中心点。由此,能够使得该像素点与包括该像素点的第一像素区域的平均灰度值之间具有更强的相关性。
在本公开第一方面所涉及的图像降噪方法中,可选地,所述规定像素行的像素点的数量为3至9。在这种情况下,通过考虑超声图像在旋转方向上受到干扰的像素点来选取合适的像素行,从而能够更有效地降低超声图像在旋转方向上噪声。
在本公开第一方面所涉及的图像降噪方法中,可选地,所述第一像素区域与所述第二像素区域相同。在这种情况下,各个干扰像素点的平均灰度值和中值滤波基于相同的像素区域,从而能够更有效地对图像进行降噪处理。
在本公开第一方面所涉及的图像降噪方法中,可选地,当所述超声图像的像素点靠近边缘且无法选择所述第一像素区域时,以该边缘为对称轴进行镜像处理以选取所述第一像素区域。在这种情况下,通过考虑目标像素点和镜像像素点,从而能够在像素点靠近边缘时选取合适的第一像素区域。
在本公开第一方面所涉及的图像降噪方法中,可选地,当所述超声图像的干扰像素点靠近边缘且无法选择所述第二像素区域时,以该边缘为对称轴进行镜像处理以选取所述第二像素区域。在这种情况下,通过考虑目标像素点和镜像像素点,从而能够在干扰像素点靠近边缘时选取合适的第二像素区域。
在本公开第一方面所涉及的图像降噪方法中,可选地,在所述第二像素区域中,所述干扰像素点位于所述规定像素行的中央。由此,能够更准确地对各个干扰像素点进行中值滤波,从而能够更有效地对图像进行降噪处理。
本公开第二方面提供一种血管内超声系统,其特征在于,使用本公开第一方面提供的图像降噪方法对血管内超声采集的超声图像进行降噪处理。
根据本公开,能够准确地识别出超声图像中的干扰像素点,并且能够有针对性地对所识别的干扰像素点进行中值滤波处理以降低因干扰而造成的噪声,从而能够优化降噪处理过程并且有效地对超声图像中的干扰像素点进行降噪。
附图说明
现在将仅通过参考附图的例子进一步详细地解释本公开,其中:
图1是示出了本公开示例所涉及的血管内超声系统的应用场景示意图。
图2是示出了本公开示例所涉及的超声探头在血管内发射超声声束的示意图。
图3是示出了本公开示例所涉及的血管内超声系统1的框图示意图。
图4(a)是示出了本公开示例所涉及的超声图像的实时图示例,图4(b)是示出了图4(a)所示的实时图示例的展开图。
图5(a)是示出了图4(a)所示的实时图示例的像素点示意图,图5(b)是示出了图4(b)所示的展开图的像素点示意图。
图6是示出了本公开示例所涉及的图像降噪方法的流程图。
图7是示出了本公开示例所涉及的第一像素区域的例子的示意图。
图8是示出了本公开示例所涉及的超声图像进行镜像处理以选取第一像素区域的示意图。
图9是示出了本公开示例所涉及的第二像素区域的例子的示意图。
图10是示出了本公开示例所涉及的超声图像进行镜像处理以选取第二像素区域的示意图。
附图标记说明:
1…血管内超声系统,10…超声探头,20…合成模块,30…降噪模块,40…显示模块,2…血管,P1…第一像素区域,P2…第二像素区域。
具体实施方式
以下,参考附图,详细地说明本公开的优选实施方式。在下面的说明中,对于相同的部件赋予相同的符号,省略重复的说明。另外,附图只是示意性的图,部件相互之间的尺寸的比例或者部件的形状等可以与实际的不同。
需要说明的是,本公开中的术语“包括”和“具有”以及它们的任何变形,例如所包括或所具有的一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元,而是可以包括或具有没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
另外,在本公开的下面描述中涉及的小标题等并不是为了限制本公开的内容或范围,其仅仅是作为阅读的提示作用。这样的小标题既不能理解为用于分割文章的内容,也不应将小标题下的内容仅仅限制在小标题的范围内。
本公开涉及一种能够优化降噪处理并且有效地对超声图像中的干扰像素点进行降噪的血管内超声系统的图像降噪方法。血管内超声系统可以简称为IVUS系统,血管内超声系统的图像降噪方法可以简称为图像降噪方法。以下结合附图对图像降噪方法进行详细说明。
图1是示出了本公开示例所涉及的血管内超声系统1的应用场景示意图。图2是示出了本公开示例所涉及的超声探头10在血管2内发射超声声束的示意图。图3是示出了本公开示例所涉及的血管内超声系统1的框图示意图。
在本实施方式中,如图1、图2和图3所示,血管内超声系统1可以包括超声探头10、合成模块20、降噪模块30和显示模块40。超声探头10可以在血管2内同时或分时进行旋转和/或沿着血管2的长度方向进行平移,例如沿着如图2所示的旋转方向R进行旋转、沿着如图2所示的长度方向L进行平移。并且,超声探头10在进行旋转或平移时可以在血管2内发射超声声束并可以接收该超声声束在血管2内的反射波,例如超声声束在血管2内传播至管壁时的反射波。接着,合成模块20可以基于该反射波转换所形成的电信号生成用于显示血管2的组织结构和几何形态的超声图像。然后,降噪模块30可以使用本公开所涉及的图像降噪方法对该超声图像进行降噪处理,以降低例如因呼吸机所产生的脉冲干扰而造成的噪声从而有效地提高超声图像的图像质量。最后,显示模块40可以用于方便地显示该超声图像。
另外,在本实施方式中,超声探头10可以将电信号转换为声信号,也可以将声信号转换为电信号。具体而言,超声探头10可以将血管内超声系统1的电信号转换为声信号,从而在血管2内发射超声声束。超声探头10还可以接收该超声声束在血管2内的反射波,并将该反射波转换为电信号,血管内超声系统1可以进一步处理该电信号从而生成上述超声图像。
另外,在一些示例中,超声探头10发射的超声声束的中心频率可以为20MHz至60MHz。在一些示例中,优选地,超声探头10发射的超声声束的中心频率可以为40MHz至60MHz。由此,能够在血管2内高质量地获取超声图像。
图4(a)是示出了本公开示例所涉及的超声图像的实时图示例,图4(b)是示出了图4(a)所示的实时图示例的展开图。
在一些示例中,如图4(a)所示,超声探头10在血管2内所获取的超声图像的实时图例如显示模块40所显示的实时图可以是血管2的截面图,如图所示R方向即为超声探头10在血管2内进行旋转的旋转方向,如图所示D方向即为血管2的径向方向,该截面图包括血管2的管壁和管腔的截面(断面)信息。通过该截面图,医生等能够方便地获知血管2的断面形态和血流图形等信息。另外,在一些示例中,为了便于对图4(a)所示的截面图进行例如降噪处理,图4(a)所示的截面图可以展开以得到图4(b)所示的展开图。
图5(a)是示出了图4(a)所示的实时图示例的像素点示意图,图5(b)是示出了图4(b)所示的展开图的像素点示意图。
在一些示例中,如图5(a)所示,超声探头10在血管2内所获取的超声图像可以具有沿着超声探头10的旋转方向的像素行,例如沿着如图5(a)所示的旋转方向R而形成像素点B1、像素点B2和像素点B3。另外,上述超声图像还可以具有沿着血管2的径向方向的像素列,例如沿着如图5(a)所示的径向方向D而形成像素点B3、像素点B4和像素点B5,其中,血管2的径向方向是指沿着血管2管壁的一侧穿刺至相对的另一侧的方向。另外,在一些示例中,图5(a)所示的像素点示意图可以展开以得到如图5(b)所示的像素点示意图,从而便于对图4(a)所示的截面图进行例如降噪处理。
图6是示出了本公开示例所涉及的图像降噪方法的流程图。图7是示出了本公开示例所涉及的第一像素区域P1的例子的示意图。图8是示出了本公开示例所涉及的超声图像进行镜像处理以选取第一像素区域P1的示意图。图9是示出了本公开示例所涉及的第二像素区域P2的例子的示意图。图10是示出了本公开示例所涉及的超声图像进行镜像处理以选取第二像素区域P2的示意图。
如图6所示,本实施方式所涉及的图像降噪方法可以包括以下步骤:遍历超声图像中的各个像素点,并获取各个像素点的原始灰度值(步骤S100);选定任一像素点,选取包括该像素点的第一像素区域P1,计算第一像素区域P1的各个像素点的原始灰度值的平均值作为第一像素区域P1的平均灰度值(步骤S200);比较各个像素点的原始灰度值与包括该像素点的第一像素区域P1的平均灰度值以获得比较结果(步骤S300);根据比较结果标记出超声图像中的干扰像素点(步骤S400);在超声图像中遍历各个干扰像素点,并选取包括该干扰像素点的第二像素区域P2,计算第二像素区域P2中各个像素点的原始灰度值的中位值(步骤S500);使用该中位值替换该干扰像素点的原始灰度值(步骤S600),在步骤S500中,第二像素区域P2可以包括含有该干扰像素点且沿着旋转方向的规定像素行、以及含有该干扰像素点且沿着血管的径向方向的规定像素列。
在本实施方式中,通过比较超声图像中各个像素点的原始灰度值与包括该像素点的第一像素区域P1的平均灰度值,能够准确地识别出超声图像中的干扰像素点,并且能够有针对性地对所识别的干扰像素点进行中值滤波处理以降低因脉冲干扰而产生的噪声,从而能够优化降噪处理过程并且有效地对超声图像中的干扰像素点进行降噪。通过本实施方式的图像降噪方法对超声图像进行降噪处理,能够有效地提高超声图像的质量,从而有助于医生等对患者做出更准确的诊断。
在一些示例中,血管内超声系统1可以使用上述图像降噪方法,即通过上述步骤S100至步骤S600,以实现对超声图像的降噪处理。
(步骤S100)
在本实施方式中,如上所述,在步骤S100中,可以遍历超声图像中各个像素点,并获取各个像素点的原始灰度值。
在本实施方式中,超声图像通常是灰度图像。灰度图像是指图像中各个像素点的R、G、B三个通道的亮度级别相等,例如RGB(100,100,100),数值“100”即为该像素点的灰度值。另外,在本实施方式中,原始灰度值是指超声图像未经任何滤波处理时其像素点的灰度值。
在一些示例中,遍历超声图像中各个像素点时可以沿着旋转方向上的像素行逐行进行浏览,也可以沿着径向方向上的像素列逐列进行浏览。在这种情况下,通过规律地浏览,能够有效降低遗漏超声图像中像素点的可能性。
另外,在一些示例中,在获取各个像素点的原始灰度值之前可以预先获取超声图像中所包含的像素点的数量。在这种情况下,通过对比像素点的数量与所获取的原始灰度值的数量,能够简单方便地预判在获取原始灰度值时是否对超声图像中的像素点有所遗漏。
另外,在一些示例中,在获取各个像素点的原始灰度值之后可以将该原始灰度值记录在其对应的像素点上。由此,能够方便地读取各个像素点的原始灰度值。
接着,在记录超声图像中各个像素点的原始灰度值后,执行步骤S200。
(步骤S200)
在本实施方式中,如上所述,在步骤S200中,可以选定任一像素点,并选取包括该像素点的第一像素区域P1,计算第一像素区域P1的各个像素点的原始灰度值的平均值作为第一像素区域P1的平均灰度值。
在本实施方式中,如图7所示,在超声图像中选定任一像素点700,选取包括该像素点的第一像素区域P1。计算第一像素区域P1中各个像素点的原始灰度值的平均值,也即,计算像素点700、像素点701、像素点702、像素点703、像素点704、像素点705、像素点706、像素点707、像素点708的原始灰度值的平均值,使用该平均值作为第一像素区域P1的平均灰度值。
在一些示例中,第一像素区域P1可以是如图7所示的第一像素区域P1,即3×3像素点矩阵。但本实施方式的示例并不限于此,还可以使用其他方式选取第一像素区域P1,例如3×5像素点矩阵、5×5像素点矩阵、7×7像素点矩阵等。
另外,在一些示例中,优选地,如图7所示,在选取包括选定的像素点的第一像素区域P1时,该像素点可以是第一像素区域P1的中心点。在这种情况下,能够使得选定的像素点与包括该像素点的第一像素区域P1的平均灰度值之间具有更强的相关性,从而能够更准确地标记出超声图像中的干扰像素点。
具体而言,在本公开所涉及的图像降噪方法中,比较选定的像素点的原始灰度值和包括该像素点的第一像素区域P1的平均灰度值,若二者之间的相对偏差较大,则可以认定该像素点受到脉冲干扰的影响较大。基于这个理论,选定的像素点与包括该像素点的第一像素区域P1的平均灰度值之间的相关性越强,则对干扰像素点的标记就越准确。由此,在本实施方式中,优选地,在选取包括选定的像素点的第一像素区域P1时,以该像素点作为第一像素区域P1的中心点。
另外,在一些示例中,如图8所示,当选定的像素点800靠近超声图像的边缘且无法选择第一像素区域P1时,可以以该边缘为对称轴进行镜像处理以选取第一像素区域P1。在这种情况下,通过考虑选定的像素点和镜像像素点,从而能够在像素点靠近超声图像的边缘时选取合适的第一像素区域P1(参见图8)。
以下,结合图8对当像素点靠近超声图像的边缘时对第一像素区域P1的选取方式进行详细说明:
在一些示例中,如图8所示,当选定的像素点800靠近超声图像的边缘时,若以3×3像素点矩阵的方式选取第一像素区域P1,由于像素点800周围的像素区域不足以构成3×3像素点矩阵,则会导致无法选取第一像素区域P1。在这种情况下,通过镜像处理以合理扩大像素点800周围的像素区域,例如以如图8所示的A-A'为对称轴进行镜像处理后,有效且合理地扩大了像素点800周围的像素区域,像素点800的镜像像素点为800c。由此,能够选取合适的第一像素区域P1。但本实施方式的示例并不限于此,在另一些示例中,当选定的像素点靠近超声图像边缘时,还可以基于与该超声图像边缘相邻的超声图像以选取第一像素区域P1。
接着,在基于各个像素点选取合适的第一像素区域P1并计算第一像素区域P1的平均灰度值后,执行步骤S300。
(步骤S300)
在本实施方式中,如上所述,在步骤S300中,可以比较各个像素点的原始灰度值与包括该像素点的第一像素区域P1的平均灰度值以获得比较结果。
如上所述,在本实施方式所涉及的图像降噪方法中,比较选定的像素点的原始灰度值和包括该像素点的第一像素区域P1的平均灰度值,若选定的像素点的原始灰度值和包括该像素点的第一像素区域P1的平均灰度值之间的相对偏差较大,则可以认定该像素点受到脉冲干扰的影响较大。
在一些示例中,可以通过计算各个像素点的原始灰度值与包括该像素点的第一像素区域P1的平均灰度值的比值,来获得各个像素点的原始灰度值与包括该像素点的第一像素区域P1的平均灰度值的比较结果。由此,能够方便地计算各个像素点的原始灰度值与包括该像素点的第一像素区域P1的平均灰度值之间的相对偏差。
另外,在一些示例中,可以通过计算各个像素点的原始灰度值与包括该像素点的第一像素区域P1的平均灰度值的差值,来比较各个像素点的原始灰度值与包括该像素点的第一像素区域P1的平均灰度值,并获得比较结果。由此,能够方便地计算选定的像素点的原始灰度值和包括该像素点的第一像素区域P1的平均灰度值之间的绝对偏差。
另外,在本实施方式中,优选地,以各个像素点的原始灰度值与包括该像素点的第一像素区域P1的平均灰度值的比值作为比较结果。
接着,在获得选定的像素点的原始灰度值与包括该像素点的第一像素区域P1的平均灰度值的比较结果后,执行步骤S400。
(步骤S400)
在本实施方式中,如上所述,在步骤S400中,可以根据比较结果标记出超声图像中的干扰像素点。
在一些示例中,当用于表示比较结果的比值大于第一预定阈值或小于第二预定阈值时,将该像素点标记为干扰像素点,其中,第一预定阈值大于第二预定阈值。在这种情况下,通过将上述比值设定在预定区间,从而能够更方便地识别超声图像中的干扰像素点。
换言之,当选定的像素点的原始灰度值相较于包括该像素点的第一像素区域的平均灰度值的偏差较大时,即上述比值大于第一预定阈值或小于第二预定阈值时,可以认定该像素点受到脉冲干扰的影响较大,从而可以将该像素点标记为干扰像素点。
在一些示例中,干扰像素点通常是由于布置在手术室内的各种设备(例如呼吸机、显示器等)产生脉冲干扰而引起的噪声,这种噪声的特点常常是非持续性的,其频率与超声工作频率量级接近,在超声图像上经常表现为椒盐噪声(脉冲噪声)。而且这样的脉冲噪声容易在超声探头10的旋转方向上持续3-5个像素点,在血管2的径向方向持续若干个像素点,对灰度影响比较大。出于脉冲噪声对超声图像的灰度影响的考虑,在一些示例中,上述第一预定阈值可以为1.3,上述第二预定阈值可以为0.7。出于脉冲噪声选择性的考虑,例如考虑到对不同程度的脉冲噪声进行处理,优选地,上述第一预定阈值可以为1.1,上述第二预定阈值可以为0.9。在这种情况下,通过设定合理的预定区间,即具体的偏差范围,从而能够更准确地识别超声图像中的干扰像素点。
接着,在标记出超声图像中的干扰像素点后,执行步骤S500。
(步骤S500)
在本实施方式中,步骤S500可以包括在超声图像中遍历各个干扰像素点,并选取包括该干扰像素点的第二像素区域P2,计算第二像素区域P2中各个像素点的原始灰度值的中位值。
另外,在本实施方式中,第二像素区域P2可以包括含有该干扰像素点且沿着旋转方向的规定像素行、以及含有该干扰像素点且沿着血管的径向方向的规定像素列。
另外,在本实施方式中,如图9所示,在超声图像中选定任一干扰像素点900,选取包括该干扰像素点的第二像素区域P2。计算第二像素区域P2中各个像素点的原始灰度值的中位值。也即,将像素点900、像素点901、像素点902、像素点903、像素点904的原始灰度值按大小排序,选取位于中间的值,即中位值。
在一些示例中,如图9所示,第二像素区域P2可以是包括选定的干扰像素点的一个像素行。优选地,在第二像素区域P2中,选定的干扰像素点可以位于该像素行的中央。由此,能够更准确地对各个干扰像素点进行中值滤波,从而能够更有效地对超声图像进行降噪处理。
另外,在一些示例中,规定像素行的像素点的数量可以为3至9。在这种情况下,通过考虑超声图像在旋转方向上受噪声干扰影响的像素来选取合理数量的像素点,从而能够更有效地降低旋转方向上所形成的超声图像的噪声。
在临床实践中,通常情况下,超声图像因脉冲干扰而产生的噪声在旋转方向上可以持续3-5个像素点、在径向方向上可以持续若干个像素点。在这种情况下,通过本示例所描述的方式选取第二像素区域P2,能够有效地选取包括多个未受到干扰或受到干扰较小的像素点,从而能够使得计算所得的中位值更接近于未受到干扰或受到干扰较小的像素点的原始灰度值。
在一些示例中,第二像素区域P2可以与第一像素区域P1相同。在这种情况下,各个干扰像素点的平均灰度值和中值滤波基于相同的像素区域,从而能够更有效地对超声图像进行降噪处理。
另外,在一些示例中,对第二像素区域P2的选取可以进行多次操作,并以最后一次选取的第二像素区域P2为准。具体而言,第一次选取第二像素区域P2后,比较第一次选取的第二像素区域P2中干扰像素点的数量以及非干扰像素点的数量。若干扰像素点的数量大于或等于非干扰像素点的数量,则对第二像素区域P2进行第二次选取,例如,以沿着旋转方向将左右两端同步延长的方式适当扩大规定像素行中像素点的数量。比较第二次选取的第二像素区域P2中干扰像素点的数量和非干扰像素点的数量。以此类推,直至最后一次选取的第二像素区域P2中干扰像素点的数量小于非干扰像素点的数量。将上述最后一次选取的第二像素区域P2作为步骤S500中待处理的像素区域。
在一些示例中,如步骤S200中有关选取第一像素区域P1时描述的那样,如图10所示,当选定的干扰像素点1000靠近超声图像的边缘且无法选择第二像素区域P2时,可以以该边缘为对称轴进行镜像处理以选取第二像素区域P2。在这种情况下,通过考虑选定的干扰像素点和镜像像素点,从而能够在像素点靠近超声图像的边缘时选取合适的第二像素区域P2。
以下,结合图10对当像素点靠近超声图像的边缘时对第二像素区域P2的选取方式进行详细说明:
具体而言,如图10所示,当选定的干扰像素点1000靠近超声图像的边缘时,若以包括5个像素点的规定像素行的方式选取第二像素区域P2,由于像素点1000周围的像素区域不足以构成包括5个像素点的规定像素行,则会导致无法选取第二像素区域P2。在这种情况下,通过镜像处理以合理扩大像素点1000周围的像素区域,例如以如图10所示的B-B'为对称轴进行镜像处理后,有效且合理地扩大了像素点1000周围的像素区域,像素点1000的镜像像素点为1000c。由此,能够选取合适的第二像素区域P2。但本实施方式的示例并不限于此,在另一些示例中,当选定的像素点靠近超声图像边缘时,还可以基于与该超声图像边缘相邻的超声图像以选取第一像素区域P2。
基于各个干扰像素点选取合适的第二像素区域P2并计算第二像素区域P2中各个像素点的原始灰度值的中位值后,执行步骤S600。
(步骤S600)
在本实施方式中,步骤S600可以包括使用步骤S500计算所得的中位值替换该干扰像素点的原始灰度值。
在本实施方式中,通过选取包括选定的干扰像素点的合适的第二像素区域P2,计算第二像素区域P2中各个像素点的原始灰度值的中位值,并使用该中位值替代该干扰像素点的原始灰度值。由此,能够有效降低干扰像素点因脉冲干扰而产生的噪声。
具体而言,由于第二像素区域P2中的干扰像素点相较于第二像素区域P2中的其它像素点尤其是未受到干扰或受到干扰较小的像素点,该干扰像素点的原始灰度值存在较大偏差。由此,当将第二像素区域P2的各个像素点的原始灰度值按大小排序时,干扰像素点的原始灰度值不容易出现在中间位置。处于中间位置的原始灰度值属于某个未受到干扰或者受到干扰较小的像素点,相较于干扰像素点而言,其更适用于表示超声图像中的有效图像信息。
在本实施方式中,步骤S600结束后,本公开所涉及的图像降噪方法即完成对超声图像的降噪处理。
另外,在本实施方式中,对超声图像的降噪处理不限于上述步骤S100至步骤S600。在一些示例中,在步骤S400和步骤S500之间还可以包括步骤S410。步骤S410可以包括获取超声图像中干扰像素点的数量,并比较超声图像中像素点的数量与干扰像素点的数量,若超声图像中干扰像素点的数量在超声图像中像素点的数量中占有较高的比重,即干扰比重较高,则弃用该帧超声图像。
具体而言,通常情况下,超声图像中干扰像素点不能准确地表示血管2内的实际信息。在这种情况下,若超声图像中干扰像素点较多,则该超声图像所表示的信息与血管2的实际信息存在较大出入,医生等不能很好地基于该超声图像对患者做出诊断治疗。
另外,在一些示例中,若超声图像的干扰比重达到50%及以上,则弃用该帧超声图像。在另一些示例中,考虑到诊断的严谨性,若超声图像的干扰比重达到30%及以上,则弃用该帧超声图像。
另外,在一些示例中,在步骤S500和步骤S600之间还可以包括步骤S510。步骤S510包括选定任一干扰像素点,分别计算该干扰像素点的原始灰度值、以及包括该干扰像素点的第二像素区域P2的中位值与包括该干扰像素点的第一像素区域P1的平均灰度值之间的相对偏差。若二者中包括该干扰像素点的第二像素区域P2的中位值与包括该干扰像素点的第一像素区域P1的平均灰度值之间的相对偏差更大,则使用包括该干扰像素点的第一像素区域P1的平均灰度值替换该干扰像素点的原始灰度值。
具体而言,当包括该干扰像素点的第二像素区域P2的中位值与包括该干扰像素点的第一像素区域P1的平均灰度值之间的相对偏差更大时,则该中位值不适用于表示超声图像中的有效图像信息。
另外,在本实施方式中,在步骤S100至步骤S600中,使用本公开所涉及的图像降噪方法对超声图像进行降噪处理时,仅需以一帧超声图像作为处理对象。在这种情况下,血管内超声系统1中用于降噪处理的降噪模块30,即使不设置缓存装置也能够完成对超声图像的降噪处理。由此,能够有效优化血管内超声系统1的系统结构。
另外,在本实施方式中,使用本公开所涉及的图像降噪方法对超声图像进行降噪处理时,仅对超声图像中的干扰像素点进行降噪处理。由此,能够优化降噪处理。
根据本公开,能够准确地识别出超声图像中的干扰像素点,并且能够有针对性地对所识别的干扰像素点进行中值滤波处理以降低因干扰而造成的噪声,从而能够优化降噪处理并且有效地对超声图像中的干扰像素点进行降噪。
虽然以上结合附图和示例对本公开进行了具体说明,但是可以理解,上述说明不以任何形式限制本公开。本领域技术人员在不偏离本公开的实质精神和范围的情况下可以根据需要对本公开进行变形和变化,这些变形和变化均落入本公开的范围内。
以上在具体实施方式中描述了本公开的各种示例。尽管这些描述直接描述了上述示例,但是应该理解的是,本领域技术人员可以想到对这里示出和描述的特定示例的修改和/或变形。落入本说明书范围内的任何这样的修改和/或变形也包括在其中。除非特别指出,否则发明人的意图是说明书和权利要求书中的词语和短语被赋予普通技术人员的普通和习惯的含义。
Claims (10)
1.一种IVUS系统的图像降噪方法,其特征在于,
利用所述IVUS系统的超声探头在旋转的情况下沿着旋转方向和血管的长度方向运动来获取来自于血管内的超声图像,所述超声图像具有沿着所述超声探头的旋转方向的像素行以及沿着血管的径向方向的像素列,血管的径向方向为沿着血管管壁的一侧穿刺至另一侧的方向,所述超声图像为包括血管的管壁和管腔的截面信息的截面图像;在所述超声图像中遍历各个像素点,并获取各个像素点的原始灰度值;在所述超声图像中选定任一像素点并选取包括该像素点的第一像素区域,计算所述第一像素区域的各个像素点的原始灰度值的平均值作为所述第一像素区域的平均灰度值;比较各个像素点的原始灰度值与包括该像素点的第一像素区域的平均灰度值并获得比较结果,并根据所述比较结果标记出干扰像素点;在所述超声图像中选定任一干扰像素点并选取包括该干扰像素点的第二像素区域,分别计算所选定的干扰像素点的原始灰度值、所述第二像素区域中各个像素点的原始灰度值的中位值与包括所选定的干扰像素点的第一像素区域的平均灰度值之间的相对偏差,若该中位值与该平均灰度值的相对偏差更小,则使用该中位值替换所选定的干扰像素点的原始灰度值,若该中位值与该平均灰度值之间的相对偏差更大,则使用该平均灰度值替换所选定的干扰像素点的原始灰度值。
2.根据权利要求1所述的图像降噪方法,其特征在于:
所述第一像素区域为3×3像素点矩阵、3×5像素点矩阵、5×5像素点矩阵或7×7像素点矩阵。
3.根据权利要求1所述的图像降噪方法,其特征在于:
所述第二像素区域包括含有所选定的干扰像素点且沿着旋转方向的规定像素行、以及含有所选定的干扰像素点且沿着血管的径向方向的规定像素列,
若所述第二像素区域中标记为干扰像素点的数量不小于未标记为干扰像素点的数量,则对所述第二像素区域重新进行选取。
4.根据权利要求3所述的图像降噪方法,其特征在于:
通过沿着旋转方向同步延长所述规定像素行的左右两端,以扩大所述规定像素行中像素点的数量。
5.根据权利要求1所述的图像降噪方法,其特征在于:
通过计算各个像素点的原始灰度值与包括该像素点的第一像素区域的平均灰度值的比值,来获得各个像素点的原始灰度值与包括该像素点的第一像素区域的平均灰度值的所述比较结果。
6.根据权利要求5所述的图像降噪方法,其特征在于:
若所述比值大于第一预定阈值或小于第二预定阈值,则将该像素点标记为干扰像素点,所述第一预定阈值大于所述第二预定阈值。
7.根据权利要求1所述的图像降噪方法,其特征在于:
还包括在选取所述第一像素区域之前,在所述超声图像中遍历各个像素点,以获取所述超声图像中像素点的总数量以及各个像素点的原始灰度值,若所述超声图像中像素点的总数量与所获取的原始灰度值的总数量相等,则进行后续步骤。
8.根据权利要求7所述的图像降噪方法,其特征在于:
通过沿着旋转方向上的像素行逐行进行浏览或沿着径向方向的像素列逐列进行浏览以遍历所述超声图像中的各个像素点。
9.根据权利要求1所述的图像降噪方法,其特征在于:
当所述超声图像的像素点靠近边缘且无法选择所述第一像素区域时,以该边缘为对称轴进行镜像处理以选取所述第一像素区域。
10.根据权利要求1所述的图像降噪方法,其特征在于:
通过将所述截面图像展开以对所述超声图像进行降噪处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110309664.6A CN113057676B (zh) | 2019-10-18 | 2019-10-18 | Ivus系统的图像降噪方法 |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110309664.6A CN113057676B (zh) | 2019-10-18 | 2019-10-18 | Ivus系统的图像降噪方法 |
CN201910996383.5A CN110893109B (zh) | 2019-10-18 | 2019-10-18 | 血管内超声系统的图像降噪方法 |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910996383.5A Division CN110893109B (zh) | 2019-10-18 | 2019-10-18 | 血管内超声系统的图像降噪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113057676A CN113057676A (zh) | 2021-07-02 |
CN113057676B true CN113057676B (zh) | 2022-04-29 |
Family
ID=69785781
Family Applications (5)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910996383.5A Active CN110893109B (zh) | 2019-10-18 | 2019-10-18 | 血管内超声系统的图像降噪方法 |
CN202110309664.6A Active CN113057676B (zh) | 2019-10-18 | 2019-10-18 | Ivus系统的图像降噪方法 |
CN202110310915.2A Active CN113040826B (zh) | 2019-10-18 | 2019-10-18 | 用于对超声图像进行降噪处理的降噪模块 |
CN202110309682.4A Active CN113017699B (zh) | 2019-10-18 | 2019-10-18 | 用于对超声图像进行降噪的图像降噪方法 |
CN202110309683.9A Active CN113017700B (zh) | 2019-10-18 | 2019-10-18 | 血管内超声系统 |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910996383.5A Active CN110893109B (zh) | 2019-10-18 | 2019-10-18 | 血管内超声系统的图像降噪方法 |
Family Applications After (3)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110310915.2A Active CN113040826B (zh) | 2019-10-18 | 2019-10-18 | 用于对超声图像进行降噪处理的降噪模块 |
CN202110309682.4A Active CN113017699B (zh) | 2019-10-18 | 2019-10-18 | 用于对超声图像进行降噪的图像降噪方法 |
CN202110309683.9A Active CN113017700B (zh) | 2019-10-18 | 2019-10-18 | 血管内超声系统 |
Country Status (1)
Country | Link |
---|---|
CN (5) | CN110893109B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113570507B (zh) * | 2020-04-29 | 2024-03-26 | 浙江宇视科技有限公司 | 一种图像降噪方法、装置、设备及存储介质 |
US20230314379A1 (en) * | 2022-02-11 | 2023-10-05 | Halliburton Energy Services, Inc. | Efficient beam profile imaging for non-negligible wave properties and rotationally anisotropic geometries |
CN115644933B (zh) * | 2022-11-17 | 2023-08-22 | 深圳微创踪影医疗装备有限公司 | 导管冲刷控制方法、装置、计算机设备、存储介质 |
Family Cites Families (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020133077A1 (en) * | 2001-03-14 | 2002-09-19 | Edwardsen Stephen Dodge | Transesophageal ultrasound probe having a rotating endoscope shaft |
JP2006512938A (ja) * | 2002-09-16 | 2006-04-20 | イメージング セラピューティクス,インコーポレーテッド | 筋骨格疾患の画像化マーカ |
JP2006319793A (ja) * | 2005-05-13 | 2006-11-24 | Pioneer Electronic Corp | 画像信号処理回路、表示装置及び画像信号処理方法 |
CN100484479C (zh) * | 2005-08-26 | 2009-05-06 | 深圳迈瑞生物医疗电子股份有限公司 | 超声图像增强与斑点抑制方法 |
WO2009081410A2 (en) * | 2007-12-25 | 2009-07-02 | Medic Vision - Brain Technologies Ltd. | Noise reduction of images |
US8483488B2 (en) * | 2009-08-07 | 2013-07-09 | Medinol Ltd. | Method and system for stabilizing a series of intravascular ultrasound images and extracting vessel lumen from the images |
CN101908205B (zh) * | 2010-06-09 | 2011-11-30 | 河北师范大学 | 基于幻方编码的中值滤波方法 |
CN102136068B (zh) * | 2011-03-31 | 2012-11-21 | 中国科学院半导体研究所 | 基于平均灰度的距离选通图像有效信息区的提取方法 |
CN103845077B (zh) * | 2012-12-05 | 2016-01-20 | 深圳迈瑞生物医疗电子股份有限公司 | 超声图像增益优化方法及超声成像增益自动优化装置 |
US9501846B2 (en) * | 2013-01-21 | 2016-11-22 | Kowa Company, Ltd. | Image processing device, image processing method, image processing program, and recording medium storing said program for producing noise-eliminated image |
CN103177427A (zh) * | 2013-03-14 | 2013-06-26 | 哈尔滨工程大学 | 一种x波段雷达图像同频干扰去除方法 |
CN103268630B (zh) * | 2013-05-22 | 2015-11-18 | 北京工业大学 | 一种基于血管内超声影像的血管三维可视化方法 |
CN103871034B (zh) * | 2014-03-22 | 2017-03-22 | 四川大学 | 一种自适应的图像椒盐噪声滤除方法 |
CN104809701A (zh) * | 2015-04-16 | 2015-07-29 | 南京航空航天大学 | 一种基于迭代开关中均值的图像椒盐噪声去除方法 |
CN104778669B (zh) * | 2015-04-16 | 2017-12-26 | 北京邮电大学 | 快速图像去噪方法及装置 |
CN106558029A (zh) * | 2016-10-28 | 2017-04-05 | 成都西纬科技有限公司 | 一种图像滤波方法及装置 |
CN107093167A (zh) * | 2017-03-07 | 2017-08-25 | 北京环境特性研究所 | 一种用于紫外成像系统的自适应去噪算法 |
KR101859392B1 (ko) * | 2017-10-27 | 2018-05-18 | 알피니언메디칼시스템 주식회사 | 초음파 영상 기기 및 이를 이용한 클러터 필터링 방법 |
CN108093182A (zh) * | 2018-01-26 | 2018-05-29 | 广东欧珀移动通信有限公司 | 图像处理方法和装置、电子设备、计算机可读存储介质 |
EP3775986A4 (en) * | 2018-04-09 | 2022-01-05 | Butterfly Network, Inc. | ULTRASONIC SYSTEM CONFIGURATION METHODS AND APPARATUS USING MULTIPLE IMAGING PARAMETER VALUES |
CN108764325B (zh) * | 2018-05-23 | 2022-07-08 | 腾讯科技(深圳)有限公司 | 图像识别方法、装置、计算机设备和存储介质 |
CN109064418B (zh) * | 2018-07-11 | 2022-03-08 | 成都信息工程大学 | 一种基于非局部均值的非均匀噪声图像去噪方法 |
CN109191387B (zh) * | 2018-07-20 | 2021-09-24 | 河南师范大学 | 一种基于巴特沃斯滤波器的红外图像去噪方法 |
CN109549671A (zh) * | 2018-12-31 | 2019-04-02 | 深圳北芯生命科技有限公司 | 带无线通信模块的血管内超声系统 |
CN110163219B (zh) * | 2019-04-17 | 2023-05-16 | 安阳师范学院 | 基于图像边缘识别的目标检测方法 |
-
2019
- 2019-10-18 CN CN201910996383.5A patent/CN110893109B/zh active Active
- 2019-10-18 CN CN202110309664.6A patent/CN113057676B/zh active Active
- 2019-10-18 CN CN202110310915.2A patent/CN113040826B/zh active Active
- 2019-10-18 CN CN202110309682.4A patent/CN113017699B/zh active Active
- 2019-10-18 CN CN202110309683.9A patent/CN113017700B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN113057676A (zh) | 2021-07-02 |
CN113017699A (zh) | 2021-06-25 |
CN110893109A (zh) | 2020-03-20 |
CN113017700B (zh) | 2022-05-03 |
CN113017699B (zh) | 2022-05-03 |
CN110893109B (zh) | 2021-05-07 |
CN113040826A (zh) | 2021-06-29 |
CN113017700A (zh) | 2021-06-25 |
CN113040826B (zh) | 2022-05-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113057676B (zh) | Ivus系统的图像降噪方法 | |
JP5519891B2 (ja) | 有効な輪郭線作成法及びシステムを使用する血管内の超音波分析 | |
CN110678933B (zh) | 超声临床特征检测与相关联的设备、系统和方法 | |
EP1614387B1 (en) | Ultrasonic diagnostic apparatus, image processing apparatus and image processing method | |
US8094893B2 (en) | Segmentation tool for identifying flow regions in an image system | |
JP4921826B2 (ja) | 超音波診断装置及びその制御方法 | |
WO2015190180A1 (ja) | 医用診断装置、医用診断装置の作動方法および医用診断装置の作動プログラム | |
US20060079772A1 (en) | Ultrasonic image generating method | |
WO2017104627A1 (ja) | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム | |
WO2016098429A1 (ja) | 超音波観測装置 | |
US20190008483A1 (en) | Ultrasound observation apparatus, method of operating ultrasound observation apparatus, and computer readable recording medium | |
JP2004267584A (ja) | 超音波診断装置 | |
US20190142385A1 (en) | Ultrasonic observation apparatus, method of operating ultrasonic observation apparatus, and program for operating ultrasonic observation apparatus | |
JP4651379B2 (ja) | 超音波画像処理装置及び超音波画像処理方法、並びに、超音波画像処理プログラム | |
US12002201B2 (en) | Ultrasonic diagnostic device, operation method of ultrasonic diagnostic device, and computer readable recording medium | |
WO2023228564A1 (ja) | 画像切り出し支援装置、超音波診断装置および画像切り出し支援方法 | |
WO2022239530A1 (ja) | 画像処理装置、画像処理システム、画像処理方法、及び画像処理プログラム | |
EP4327750A1 (en) | Guided ultrasound imaging for point-of-care staging of medical conditions | |
JP7438038B2 (ja) | 超音波診断装置及び診断支援方法 | |
WO2021152745A1 (ja) | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム | |
WO2023189260A1 (ja) | コンピュータプログラム、情報処理装置及び情報処理方法 | |
JP5932189B1 (ja) | 超音波観測装置 | |
WO2021042242A1 (zh) | 一种超声成像设备及其超声回波信号的处理方法 | |
WO2024042044A1 (en) | Guided ultrasound imaging for point-of-care staging of medical conditions |
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 | ||
CB02 | Change of applicant information | ||
CB02 | Change of applicant information |
Address after: 518000 area E, 3rd floor, building 3, Tingwei Industrial Park, No.6 Liufang Road, Xin'an street, Shenzhen City, Guangdong Province Applicant after: Shenzhen Beixin Life Technology Co.,Ltd. Address before: 518000 area E, 3rd floor, building 3, Tingwei Industrial Park, No.6 Liufang Road, Xin'an street, Shenzhen City, Guangdong Province Applicant before: INSIGHT LIFETECH Co.,Ltd. (SHENZHEN ) |
|
GR01 | Patent grant | ||
GR01 | Patent grant |