CN111358492A - 一种四维造影图像生成方法、装置、设备及存储介质 - Google Patents
一种四维造影图像生成方法、装置、设备及存储介质 Download PDFInfo
- Publication number
- CN111358492A CN111358492A CN202010129489.8A CN202010129489A CN111358492A CN 111358492 A CN111358492 A CN 111358492A CN 202010129489 A CN202010129489 A CN 202010129489A CN 111358492 A CN111358492 A CN 111358492A
- Authority
- CN
- China
- Prior art keywords
- ultrasonic probe
- contrast data
- contrast
- dimensional
- swing
- 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.)
- Pending
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/48—Diagnostic techniques
- A61B8/481—Diagnostic techniques involving the use of contrast agent, e.g. microbubbles introduced into the bloodstream
-
- 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/5223—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
-
- 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)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Veterinary Medicine (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Hematology (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本申请公开了一种四维造影图像生成方法,该方法可以包括以下步骤:获得目标部位当前帧的造影数据;根据超声探头的定位信息,确定超声探头的摆动是否达到预设修正条件;如果达到,则对造影数据进行修正,并基于修正后的造影数据,生成四维时间渲染图像;如果未达到,则基于造影数据,生成四维时间渲染图像。应用本申请实施例所提供的技术方案,可以有效避免超声探头的摆动造成图像偏移,从而导致的时间成像的渲染结果颜色错乱的问题,同时,无需进行耗时低效的三维特征匹配操作来实现造影图像偏移的调整,提高了四维造影图像生成效率。本申请还公开了一种四维造影图像生成装置、设备及存储介质,具有相应技术效果。
Description
技术领域
本申请涉及图像处理技术领域,特别是涉及一种四维造影图像生成方法、装置、设备及存储介质。
背景技术
随着医学技术的进步,造影技术逐渐发展起来。造影,是在放射诊断学中,通过摄入含原子序数高的元素的物质,然后在待诊断体内部位摄取放射照片以供医学诊断。对缺乏自然对比的结构或器官,可将密度高于或低于该结构或器官的物质即造影剂引入器官内或其周围间隙,使之产生对比显影。
为了方便直观观察造影剂流动的方法与过程,可以在成像过程中,采用多段颜色伪彩,将造影剂随时间的流动与颜色一一对应。但是,如果在图像采集过程中,超声探头的摆动过大,则将容易出现难以对三维图像进行有效匹配,渲染颜色错乱的问题。
发明内容
本申请的目的是提供一种四维造影图像生成方法、装置、设备及存储介质,以避免生成的四维造影图像出现颜色错乱问题。
为解决上述技术问题,本申请提供如下技术方案:
一种四维造影图像生成方法,包括:
获得目标部位当前帧的造影数据;
根据超声探头的定位信息,确定所述超声探头的摆动是否达到预设修正条件;
如果确定所述超声探头的摆动达到所述预设修正条件,则对所述造影数据进行修正,并基于修正后的所述造影数据,生成四维时间渲染图像;
如果确定所述超声探头的摆动未达到所述预设修正条件,则基于所述造影数据,生成所述四维时间渲染图像。
在本申请的一种具体实施方式中,所述超声探头上设置有定位装置,所述根据超声探头的定位信息,确定所述超声探头的摆动是否达到预设修正条件,包括:
根据超声探头的定位信息,确定当前帧所述超声探头的空间位置;
将所述当前帧超声探头的空间位置与预先获得的首帧所述超声探头的空间位置进行比对,确定所述超声探头的空间旋转角度和平移距离;
根据所述空间旋转角度和平移距离,确定所述超声探头的摆动是否达到预设修正条件。
在本申请的一种具体实施方式中,所述根据所述空间旋转角度和平移距离,确定所述超声探头的摆动是否达到预设修正条件,包括:
如果所述空间旋转角度大于预设角度阈值和/或所述平移距离大于预设距离阈值,则确定所述超声探头的摆动达到预设修正条件。
在本申请的一种具体实施方式中,所述对所述造影数据进行修正,包括:
将电磁定位坐标系确定为世界坐标系;
将超声图像坐标系转换到所述电磁定位坐标系;
确定所述超声探头摆动的仿射矩阵;
结合所述仿射矩阵,在所述电磁定位坐标系上对所述造影数据进行坐标变换,并转换回所述超声图像坐标系,获得修正后的所述造影数据。
在本申请的一种具体实施方式中,在所述获得目标部位当前帧的造影数据之后、所述根据超声探头的定位信息,确定所述超声探头的摆动是否达到预设修正条件之前,还包括:
确定所述造影数据中每个体素是否为噪声点;
在所述造影数据中剔除所述噪声点。
在本申请的一种具体实施方式中,所述确定所述造影数据中每个体素是否为噪声点,包括:
确定所述造影数据中的每个体素对应的灰度均值;
将对应的灰度均值小于预设的灰度阈值的体素确定为噪声点。
在本申请的一种具体实施方式中,所述确定所述造影数据中的每个体素对应的灰度均值,包括:
确定所述造影数据中的每个体素对应的设定大小的区域中体素的灰度和;
将每个体素对应的灰度和与相应区域中体素个数的比值,确定为每个体素对应的灰度均值。
在本申请的一种具体实施方式中,在对已获得的所述目标部位的完整造影数据进行后处理的情况下,在所述生成四维时间渲染图像之前,还包括:
基于所述目标部位的组织特征,根据三维区域生长法在所述目标部位的完整造影数据中分割出所述目标部位。
一种四维造影图像生成装置,包括:
造影数据获得模块,用于获得目标部位当前帧的造影数据;
判断模块,用于根据超声探头的定位信息,确定所述超声探头的摆动是否达到预设修正条件;
造影数据修正模块,用于在确定所述超声探头的摆动达到所述预设修正条件时,对所述造影数据进行修正;
造影图像生成模块,用于基于修正后的所述造影数据,生成四维时间渲染图像;
所述造影图像生成模块,还用于在确定所述超声探头的摆动未达到所述预设修正条件时,基于所述造影数据,生成所述四维时间渲染图像。
一种四维造影图像生成设备,包括:
存储器,用于存储计算机程序;
处理器,用于执行所述计算机程序时实现上述任一项所述四维造影图像生成方法的步骤。
一种计算机可读存储介质,所述计算机可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时实现上述任一项所述四维造影图像生成方法的步骤。
应用本申请实施例所提供的技术方案,获得目标部位当前帧的造影数据后,根据当前帧超声探头的空间位置,确定超声探头的摆动是否达到预设修正条件,如果达到,则可以对造影数据进行修正,并基于修正后的造影数据,生成四维造影图像,如果未达到,则基于造影数据,生成四维造影图像。在超声探头的摆动达到修正条件时,对当前帧的造影数据进行修正后再生成四维造影图像,这样可以有效避免超声探头的摆动造成图像偏移,从而导致的时间成像的渲染结果颜色错乱的问题,同时,无需进行耗时低效的三维特征匹配操作来实现造影图像偏移的调整,提高了四维造影图像生成效率。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本申请实施例中一种四维造影图像生成方法的实施流程图;
图2为本申请实施例中三维仿射变换前后图像对比图;
图3为本申请实施例中对造影数据进行预处理前后对比图;
图4为本申请实施例中目标部位分割前后对比图;
图5为本申请实施例中一种四维造影图像生成装置的结构示意图;
图6为本申请实施例中一种四维造影图像生成设备的结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请方案,下面结合附图和具体实施方式对本申请作进一步的详细说明。显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。
参见图1所示,为本申请实施例所提供的一种四维造影图像生成方法的实施流程图,该方法可以包括以下步骤:
S110:获得目标部位当前帧的造影数据。
目标部位为待诊断体内部位,如子宫输卵管等部位。可以向目标部位进行造影剂推注进而获得造影数据。在向目标部位进行造影剂的推注过程中,可以对目标部位的造影进行扫查,获取当前帧的造影数据。该造影数据可以是三维造影图像。
S120:根据超声探头的定位信息,确定超声探头的摆动是否达到预设修正条件。
造影数据可以通过超声探头采集得到。在获得目标部位当前帧的造影数据的同时,可以获得超声探头的定位信息。具体的,可以采用探头定位技术得到超声探头的定位信息。通过定位信息,可以确定当前帧超声探头的空间位置。
目前常用的空间定位技术主要有光学定位和电磁定位,电磁定位技术具有精度高、操作方便灵活、无视线遮挡问题等优点。电磁定位***,主要包括电磁发射器和探头接收器,电磁发射器通过探头接收器的反馈信号,采集得到超声探头的定位信息。
根据超声探头的定位信息,可以确定超声探头的摆动是否过大,是否达到预设修正条件。该预设修正条件可以根据实际情况进行设定和调整。
因为超声探头的摆动会造成图像偏移,从而导致时间成像的渲染结果颜色错乱。所以,本申请实施例中,当超声探头的摆动达到预设修正条件时,可以继续执行步骤S130的操作,对当前帧的造影数据进行修正处理,如果超声探头的摆动未达到预设修正条件,则对当前帧的造影数据可以不做修正处理,执行步骤S140的操作,直接利用当前帧的造影数据进行四维时间渲染图像的生成。
S130:如果确定超声探头的摆动达到预设修正条件,则对造影数据进行修正,并基于修正后的造影数据,生成四维时间渲染图像。
根据超声探头的定位信息,可以确定超声探头的摆动是否过大,是否达到预设修正条件。如果确定超声探头的摆动达到预设修正条件,则可以对当前帧的造影数据进行修正。如可以通过空间坐标变换等方式对当前帧的造影数据进行修正。
在对造影数据进行修正后,可以基于修正后的造影数据,生成四维时间造影图像。具体的,可以将当前帧修正后的造影数据与前一帧的造影数据进行比对,得到新的造影区域,在新出现的造影区域内填充相应颜色,生成四维时间渲染图像。这里前一帧的造影数据可以是前一帧的不需要修正处理的造影数据还可以是前一帧的经过修正处理的造影数据。
S140:如果确定超声探头的摆动未达到修正条件,则基于造影数据,生成四维时间渲染图像。
在确定超声探头的摆动未达到修正条件时,不需要对当前帧的造影数据做修正处理,可以直接基于当前帧的造影数据,生成四维时间渲染图像。同样,可以将当前帧的造影数据与前一帧的造影数据进行比对,得到新的造影区域,在新出现的造影区域内填充相应颜色,生成四维时间渲染图像。不同时间对应的造影区域可以对应不同颜色。这里前一帧的造影数据可以是前一帧的不需要修正处理的造影数据还可以是前一帧的经过修正处理的造影数据。
在生成四维时间渲染图像过程中,可以采用多段颜色伪彩,将造影剂随时间的流动与颜色一一对应,通过四维时间渲染图像可以方便直观观察造影剂流动的方向与过程。
对于目标部位的每一帧造影数据均按照上述步骤进行,可以得到最终的目标部位的四维时间渲染图像。
当然,在实际应用中,可以先针对于每一帧造影数据,确定是否对其进行修正处理,并在确定需要对其进行修正处理时,对其进行修正处理。最后再基于最终获得的针对目标部位的造影数据,生成四维时间渲染图像。
应用本申请实施例所提供的方法,获得目标部位当前帧的造影数据后,根据超声探头的定位信息,确定超声探头的摆动是否达到预设修正条件,如果达到,则可以对造影数据进行修正,并基于修正后的造影数据,生成四维时间渲染图像,如果未达到,则基于造影数据,生成四维时间渲染图像。因为超声探头的摆动会造成图像偏移,如果超声探头的摆动过大,则将使得图像偏移幅度较大,从而会导致时间成像的渲染结果出现颜色错乱的情况。在超声探头的摆动过大,达到预设修正条件时,对当前帧的造影数据进行修正后,再基于修正后的造影数据生成四维时间渲染图像,可以有效避免超声探头的摆动造成图像偏移,从而导致的时间成像的渲染结果颜色错乱的问题,同时,不需要进行耗时低效的三维特征匹配操作来实现造影图像偏移的调整,提高了四维时间渲染图像的生成效率。
在实际应用中,可以在向目标部位推注造影剂的过程中,实时获得目标部位的每一帧的造影数据,实时对于每一帧的造影数据进行相应处理,生成四维时间渲染图像。还可以在向目标部位推注造影剂的过程中,记录相关数据,如造影数据、超声探头的定位信息等,在后期需要进行四维时间渲染图像的生成时,基于已记录的相关数据,获得目标部位每一帧的造影数据,并针对每一帧的造影数据执行本申请步骤,生成四维时间渲染图像。还可以在向目标部位推注造影剂的过程中,实时获得目标部位的每一帧的造影数据,实时对于每一帧的造影数据进行相应处理,并记录相关数据,如不需要处理的造影数据、进行过处理的造影数据等,在后期需要进行四维时间渲染图像的生成时,基于已记录的相关数据,生成四维时间渲染图像。
在本申请的一个实施例中,超声探头上设置有定位装置,步骤S120可以包括以下步骤:
步骤一:根据超声探头的定位信息,确定当前帧超声探头的空间位置;
步骤二:将当前帧超声探头的空间位置与预先获得的首帧超声探头的空间位置进行比对,确定超声探头的空间旋转角度和平移距离;
步骤三:根据空间旋转角度和平移距离,确定超声探头的摆动是否达到预设修正条件。
为便于描述,将上述三个步骤结合起来进行说明。
在本申请实施例中,超声探头上可以设置有定位装置,该定位装置可以是电磁定位***中的探头接收器,电磁定位***的电磁发射器通过探头接收器的反馈信号,可以采集得到超声探头的定位信息。
在实际应用中,可以通过电磁定位***预先获得首帧超声探头的空间位置。在获得目标部位当前帧的造影数据时,可以根据超声探头的定位信息,确定当前帧超声探头的空间位置。
将当前帧超声探头的空间位置与首帧超声探头的空间位置进行比对,根据比对结果,可以确定超声探头的空间旋转角度和平移距离。即以首帧超声探头的空间位置为基准,确定当前帧超声探头的空间旋转角度和平移距离。
根据确定的空间旋转角度和平移距离,可以确定超声探头的摆动是否达到预设修正条件。具体的,如果空间旋转角度大于预设角度阈值和/或平移距离大于预设距离阈值,则可以确定超声探头的摆动达到预设修正条件。角度阈值和距离阈值可以根据实际情况进行设定和调整。
即在超声探头的空间旋转角度大于预设角度阈值和/或平移距离大于预设距离阈值时,可以认为超声探头的摆动过大,会造成图像偏移,导致时间成像的渲染结果颜色错乱,在这种情况下,可以确定超声探头的摆动达到了预设修正条件。
通过将当前帧超声探头的空间位置与预先获得的首帧超声探头的空间位置进行比对,确定超声探头的空间旋转角度和平移距离,从而根据空间旋转角度和平移距离,可以准确确定超声探头的摆动是否达到预设修正条件。
在本申请的一个实施例中,步骤S130对造影数据进行修正,可以包括以下步骤:
第一个步骤:将电磁定位坐标系确定为世界坐标系;
第二个步骤:将超声图像坐标系转换到电磁定位坐标系;
第三个步骤:确定超声探头摆动的仿射矩阵;
第四个步骤:结合仿射矩阵,在电磁定位坐标系上对造影数据进行坐标变换,并转换回超声图像坐标系,获得修正后的造影数据。
为便于描述,将上述四个步骤结合起来进行说明。
在本申请实施例中,可以使用电磁定位***进行超声探头的定位信息。电磁定位***中的电磁发射器位置是不动的,可以将电磁定位坐标系确定为世界坐标系。在确定超声探头的摆动达到预设修正条件时,可以对造影数据进行修正。超声探头的摆动在电磁定位坐标系,可以先将超声图像坐标系转换到电磁定位坐标系,然后确定超声探头摆动的仿射矩阵。
结合仿射矩阵,在电磁定位坐标系上对造影数据进行坐标变换,并转换回超声图像坐标系,即可获得修正后的造影数据。
举例而言,通过计算得到的超声探头的三维仿射矩阵如下:
像素P1(x1,y1,z1)和P2(x2,y2,z2)变换关系为:P2=M·P1。变换前后的图像对比如图2所示。
在本申请的一个实施例中,在步骤S110获得目标部位当前帧的造影数据之后、步骤S120根据超声探头的定位信息,确定超声探头的摆动是否达到预设修正条件之前,该方法还可以包括以下步骤:
步骤一:确定造影数据中每个体素是否为噪声点;
步骤二:在造影数据中剔除噪声点。
为便于描述,将上述两个步骤结合起来进行说明。
在本申请实施例中,获得目标部位当前帧的造影数据之后,可以先对造影数据进行降噪等预处理。
在对造影数据进行预处理时,可以先确定造影数据中每个体素是否为噪声点。具体的,可以先确定造影数据中的每个体素对应的灰度均值,将对应的灰度均值小于预设的灰度阈值的体素确定为噪声点。
即针对造影数据中的每个体素,可以先确定该体素对应的灰度均值,如果该体素对应的灰度均值大于或等于预设的灰度阈值,则可以确定该体素为非噪声点,如果该体素对应的灰度均值小于该预设的灰度阈值,则可以确定该体素为噪声点。针对于造影数据中的每个体素都可以按照这样的方式进行噪声点的确定。灰度阈值可以根据实际情况进行设定和调整,本申请实施例对此不做限制。
在确定造影数据中的每个体素对应的灰度均值时,可以先确定造影数据中的每个体素对应的设定大小的区域中体素的灰度和,然后将每个体素对应的灰度和与相应区域中体素个数的比值,确定为每个体素对应的灰度均值。
即针对于造影数据中的每个体素,确定该体素对应的设定大小的区域中体素的灰度和,将该灰度和与该区域中体素个数的比值,确定为该体素对应的灰度均值。设定大小可以根据实际情况进行设定和调整,在进行噪声点的确定时均使用统一标准即可。
每个体素具有各自对应的灰度均值,再基于灰度均值确定每个体素是否为噪声点。
在对造影数据中的每个体素进行噪声点的确定之后,可以在造影数据中剔除噪声点,得到降噪后的造影数据,并对降噪后的造影数据进行进一步操作,如修正等操作,最后生成四维时间渲染图像。这种降噪预处理方式不同于传统滤波运算,不会影响感兴趣图像区域的分辨率。如图3所示即为预处理前后对比图。
为便于理解,举例说明。
假设当前帧的造影数据I(x,y,z),设定大小的区域为(a,b,c),对于造影数据中的每个体素I(x0,y0,z0),可以计算得到该体素对应的设定大小的区域中体素的灰度和:
再计算该体素对应的灰度和与该区域中体素个数的比值,将该比值确定为该体素对应的灰度均值。
如果该灰度均值大于或等于预设的灰度阈值,则可以确定该体素为非噪声点,I(x0,y0,z0)可以保持不变,如果该灰度均值小于该预设的灰度阈值,则可以确定该体素为噪声点,可以令I(x0,y0,z0)=0。
在对造影数据进行修正、进行四维时间渲染图像的生成之前,先对造影数据进行降噪等预处理,可以减少造影数据中的噪声,避免出现最后生成的四维时间渲染图像的渲染结果颜色杂乱的问题。
在本申请的一个实施例中,在对已获得目标部位的完整造影数据进行后处理的情况下,在生成四维时间渲染图像之前,该方法还可以包括以下步骤:
基于目标部位的组织特征,根据三维区域生长法在目标部位的完整造影数据中分割出目标部位。
在本申请实施例中,可以先获得目标部位的完整造影数据,然后在对完整造影数据进行回调,生成四维时间渲染图像。在对已获得的目标部位的完整造影数据进行后处理的情况下,在生成四维时间渲染图像之前,可以基于目标部位的组织特征,在造影数据中分割出目标部位。具体的,可以根据三维区域生长法在目标部位的完整造影数据中分割出目标部位。这样可以有效去除噪声,加快后续操作的运算速度。
以目标部位为子宫输卵管为例,由于子宫输卵管形状特征和灰度特征等组织特征比较明显,可使用三维区域生长法分割出子宫输卵管。使用三维区域生长法对子宫输卵管进行分割主要采用的是三维连通域标记法,这种方法计算快速、便捷。具体的,可以根据连通区域大小进行排序,得到区域体积最大的前n个区域。由组织特征可知,子宫输卵管一般是一个或几个大的连通区域,而噪声是分布零散的小区域,因此子宫输卵管与噪声可以通过连通区域的体积大小进行区分。其中,连通区域体积大小可以通过连通区域体素个数计算得到。在实际应用中,n值大小可以通过后续临床图像得到一个较为合适的值,也可以根据实际情况进行调节。
在造影数据中分割出目标部位,也就是在造影数据中去除了非目标部位的区域,进一步减少了造影数据中的噪声,如图4所示为目标部位分割前后对比图,后续只需要对分割出的目标部位的造影数据进行四维时间渲染图像的生成操作,可以加快运算速度。
相应于上面的方法实施例,本申请实施例还提供了一种四维造影图像生成装置,下文描述的一种四维造影图像生成装置与上文描述的一种四维造影图像生成方法可相互对应参照。
参见图5所示,该装置可以包括以下模块:
造影数据获得模块510,用于获得目标部位当前帧的造影数据;
判断模块520,用于根据超声探头的定位信息,确定超声探头的摆动是否达到预设修正条件;
造影数据修正模块530,用于在确定超声探头的摆动达到预设修正条件时,对造影数据进行修正;
造影图像生成模块540,用于基于修正后的造影数据,生成四维时间渲染图像;
造影图像生成模块540,还用于在确定超声探头的摆动未达到预设修正条件时,基于造影数据,生成四维时间渲染图像。
应用本申请实施例所提供的装置,获得目标部位当前帧的造影数据后,根据超声探头的定位信息,确定超声探头的摆动是否达到预设修正条件,如果达到,则可以对造影数据进行修正,并基于修正后的造影数据,生成四维时间渲染图像,如果未达到,则基于造影数据,生成四维时间渲染图像。因为超声探头的摆动会造成图像偏移,如果超声探头的摆动过大,则将使得图像偏移幅度较大,从而会导致时间成像的渲染结果出现颜色错乱的情况。在超声探头的摆动过大,达到预设修正条件时,对当前帧的造影数据进行修正后,再基于修正后的造影数据生成四维时间渲染图像,可以有效避免超声探头的摆动造成图像偏移,从而导致的时间成像的渲染结果颜色错乱的问题,同时,不需要进行耗时低效的三维特征匹配操作来实现造影图像偏移的调整,提高了四维时间渲染图像的生成效率。
在本申请的一种具体实施方式中,超声探头上设置有定位装置,判断模块520,用于:
根据超声探头的定位信息,确定当前帧超声探头的空间位置;
将当前帧超声探头的空间位置与预先获得的首帧超声探头的空间位置进行比对,确定超声探头的空间旋转角度和平移距离;
根据空间旋转角度和平移距离,确定超声探头的摆动是否达到预设修正条件。
在本申请的一种具体实施方式中,判断模块520,用于:
如果空间旋转角度大于预设角度阈值和/或平移距离大于预设距离阈值,则确定超声探头的摆动达到预设修正条件。
在本申请的一种具体实施方式中,造影数据修正模块530,用于:
将电磁定位坐标系确定为世界坐标系;
将超声图像坐标系转换到电磁定位坐标系;
确定超声探头摆动的仿射矩阵;
结合仿射矩阵,在电磁定位坐标系上对造影数据进行坐标变换,并转换回超声图像坐标系,获得修正后的造影数据。
在本申请的一种具体实施方式中,还包括预处理模块,用于:
在获得目标部位当前帧的造影数据之后、根据超声探头的定位信息,确定超声探头的摆动是否达到预设修正条件之前,确定造影数据中每个体素是否为噪声点;
在造影数据中剔除噪声点。
在本申请的一种具体实施方式中,预处理模块,用于:
确定造影数据中的每个体素对应的灰度均值;
将对应的灰度均值小于预设的灰度阈值的体素确定为噪声点。
在本申请的一种具体实施方式中,预处理模块,用于:
确定造影数据中的每个体素对应的设定大小的区域中体素的灰度和;
将每个体素对应的灰度和与相应区域中体素个数的比值,确定为每个体素对应的灰度均值。
在本申请的一种具体实施方式中,还包括分割模块,用于:
在对已获得的目标部位的完整造影数据进行后处理的情况下,在生成四维时间渲染图像之前,基于目标部位的组织特征,根据三维区域生长法在目标部位的完整造影数据中分割出目标部位。
相应于上面的方法实施例,本申请实施例还提供了一种四维造影图像生成设备,包括:
存储器,用于存储计算机程序;
处理器,用于执行计算机程序时实现上述四维造影图像生成方法的步骤。
如图6所示,为四维造影图像生成设备的组成结构示意图,四维造影图像生成设备可以包括:处理器10、存储器11、通信接口12和通信总线13。处理器10、存储器11、通信接口12均通过通信总线13完成相互间的通信。
在本申请实施例中,处理器10可以为中央处理器(Central Processing Unit,CPU)、特定应用集成电路、数字信号处理器、现场可编程门阵列或者其他可编程逻辑器件等。
处理器10可以调用存储器11中存储的程序,具体的,处理器10可以执行四维造影图像生成方法的实施例中的操作。
存储器11中用于存放一个或者一个以上程序,程序可以包括程序代码,程序代码包括计算机操作指令,在本申请实施例中,存储器11中至少存储有用于实现以下功能的程序:
获得目标部位当前帧的造影数据;
根据超声探头的定位信息,确定超声探头的摆动是否达到预设修正条件;
如果确定超声探头的摆动达到预设修正条件,则对造影数据进行修正,并基于修正后的造影数据,生成四维时间渲染图像;
如果确定超声探头的摆动未达到预设修正条件,则基于造影数据,生成四维时间渲染图像。
在一种可能的实现方式中,存储器11可包括存储程序区和存储数据区,其中,存储程序区可存储操作***,以及至少一个功能(比如数据处理功能、图像生成功能)所需的应用程序等;存储数据区可存储使用过程中所创建的数据,如造影数据、修正数据等。
此外,存储器11可以包括高速随机存取存储器,还可以包括非易失性存储器,例如至少一个磁盘存储器件或其他易失性固态存储器件。
通信接口13可以为通信模块的接口,用于与其他设备或者***连接。
当然,需要说明的是,图6所示的结构并不构成对本申请实施例中四维造影图像生成设备的限定,在实际应用中四维造影图像生成设备可以包括比图6所示的更多或更少的部件,或者组合某些部件。
相应于上面的方法实施例,本申请实施例还提供了一种计算机可读存储介质,计算机可读存储介质上存储有计算机程序,计算机程序被处理器执行时实现上述四维造影图像生成方法的步骤。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其它实施例的不同之处,各个实施例之间相同或相似部分互相参见即可。
专业人员还可以进一步意识到,结合本文中所公开的实施例描述的各示例的单元及算法步骤,能够以电子硬件、计算机软件或者二者的结合来实现,为了清楚地说明硬件和软件的可互换性,在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本申请的范围。
结合本文中所公开的实施例描述的方法或算法的步骤可以直接用硬件、处理器执行的软件模块,或者二者的结合来实施。软件模块可以置于随机存储器(RAM)、内存、只读存储器(ROM)、电可编程ROM、电可擦除可编程ROM、寄存器、硬盘、可移动磁盘、CD-ROM、或技术领域内所公知的任意其它形式的存储介质中。
本文中应用了具体个例对本申请的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本申请的技术方案及其核心思想。应当指出,对于本技术领域的普通技术人员来说,在不脱离本申请原理的前提下,还可以对本申请进行若干改进和修饰,这些改进和修饰也落入本申请权利要求的保护范围内。
Claims (11)
1.一种四维造影图像生成方法,其特征在于,包括:
获得目标部位当前帧的造影数据;
根据超声探头的定位信息,确定所述超声探头的摆动是否达到预设修正条件;
如果确定所述超声探头的摆动达到所述预设修正条件,则对所述造影数据进行修正,并基于修正后的所述造影数据,生成四维时间渲染图像;
如果确定所述超声探头的摆动未达到所述预设修正条件,则基于所述造影数据,生成所述四维时间渲染图像。
2.根据权利要求1所述的方法,其特征在于,所述超声探头上设置有定位装置,所述根据超声探头的定位信息,确定所述超声探头的摆动是否达到预设修正条件,包括:
根据超声探头的定位信息,确定当前帧所述超声探头的空间位置;
将所述当前帧超声探头的空间位置与预先获得的首帧所述超声探头的空间位置进行比对,确定所述超声探头的空间旋转角度和平移距离;
根据所述空间旋转角度和平移距离,确定所述超声探头的摆动是否达到预设修正条件。
3.根据权利要求2所述的方法,其特征在于,所述根据所述空间旋转角度和平移距离,确定所述超声探头的摆动是否达到预设修正条件,包括:
如果所述空间旋转角度大于预设角度阈值和/或所述平移距离大于预设距离阈值,则确定所述超声探头的摆动达到预设修正条件。
4.根据权利要求1所述的方法,其特征在于,所述对所述造影数据进行修正,包括:
将电磁定位坐标系确定为世界坐标系;
将超声图像坐标系转换到所述电磁定位坐标系;
确定所述超声探头摆动的仿射矩阵;
结合所述仿射矩阵,在所述电磁定位坐标系上对所述造影数据进行坐标变换,并转换回所述超声图像坐标系,获得修正后的所述造影数据。
5.根据权利要求1所述的方法,其特征在于,在所述获得目标部位当前帧的造影数据之后、所述根据超声探头的定位信息,确定所述超声探头的摆动是否达到预设修正条件之前,还包括:
确定所述造影数据中每个体素是否为噪声点;
在所述造影数据中剔除所述噪声点。
6.根据权利要求5所述的方法,其特征在于,所述确定所述造影数据中每个体素是否为噪声点,包括:
确定所述造影数据中的每个体素对应的灰度均值;
将对应的灰度均值小于预设的灰度阈值的体素确定为噪声点。
7.根据权利要求6所述的方法,其特征在于,所述确定所述造影数据中的每个体素对应的灰度均值,包括:
确定所述造影数据中的每个体素对应的设定大小的区域中体素的灰度和;
将每个体素对应的灰度和与相应区域中体素个数的比值,确定为每个体素对应的灰度均值。
8.根据权利要求1至7之中任一项所述的方法,其特征在于,在对已获得的所述目标部位的完整造影数据进行后处理的情况下,在所述生成四维时间渲染图像之前,还包括:
基于所述目标部位的组织特征,根据三维区域生长法在所述目标部位的完整造影数据中分割出所述目标部位。
9.一种四维造影图像生成装置,其特征在于,包括:
造影数据获得模块,用于获得目标部位当前帧的造影数据;
判断模块,用于根据超声探头的定位信息,确定所述超声探头的摆动是否达到预设修正条件;
造影数据修正模块,用于在确定所述超声探头的摆动达到所述预设修正条件时,对所述造影数据进行修正;
造影图像生成模块,用于基于修正后的所述造影数据,生成四维时间渲染图像;
所述造影图像生成模块,还用于在确定所述超声探头的摆动未达到所述预设修正条件时,基于所述造影数据,生成所述四维时间渲染图像。
10.一种四维造影图像生成设备,其特征在于,包括:
存储器,用于存储计算机程序;
处理器,用于执行所述计算机程序时实现如权利要求1至8任一项所述四维造影图像生成方法的步骤。
11.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时实现如权利要求1至8任一项所述四维造影图像生成方法的步骤。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010129489.8A CN111358492A (zh) | 2020-02-28 | 2020-02-28 | 一种四维造影图像生成方法、装置、设备及存储介质 |
PCT/CN2020/096260 WO2021169086A1 (zh) | 2020-02-28 | 2020-06-16 | 一种四维造影图像生成方法、装置、设备及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010129489.8A CN111358492A (zh) | 2020-02-28 | 2020-02-28 | 一种四维造影图像生成方法、装置、设备及存储介质 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111358492A true CN111358492A (zh) | 2020-07-03 |
Family
ID=71198292
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010129489.8A Pending CN111358492A (zh) | 2020-02-28 | 2020-02-28 | 一种四维造影图像生成方法、装置、设备及存储介质 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN111358492A (zh) |
WO (1) | WO2021169086A1 (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112529804A (zh) * | 2020-12-14 | 2021-03-19 | 深圳开立生物医疗科技股份有限公司 | 一种图像去噪方法、装置、设备及计算机可读存储介质 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101102721A (zh) * | 2005-02-09 | 2008-01-09 | 株式会社日立医药 | 超声波诊断装置以及超声波摄像方法 |
CN101634551A (zh) * | 2009-08-18 | 2010-01-27 | 清华大学深圳研究生院 | 一种检测表面粗糙度的方法及其*** |
JP2010088585A (ja) * | 2008-10-06 | 2010-04-22 | Toshiba Corp | 超音波診断装置及び超音波プローブ |
CN102247171A (zh) * | 2010-03-30 | 2011-11-23 | 株式会社东芝 | 超声波诊断装置、超声波图像处理装置以及医用图像诊断装置 |
CN103473805A (zh) * | 2013-09-18 | 2013-12-25 | 内蒙古科技大学 | 基于改进区域生长算法测量三维重建肝脏模型体积的方法 |
CN104574329A (zh) * | 2013-10-09 | 2015-04-29 | 深圳迈瑞生物医疗电子股份有限公司 | 超声融合成像方法、超声融合成像导航*** |
CN106327487A (zh) * | 2016-08-18 | 2017-01-11 | 苏州润心医疗科技有限公司 | 基于x射线冠脉造影图像的冠状动脉血流储备分数计算方法 |
CN107440720A (zh) * | 2012-06-27 | 2017-12-08 | 东芝医疗***株式会社 | 超声波诊断装置以及图像数据的校正方法 |
CN109447974A (zh) * | 2018-10-31 | 2019-03-08 | 上海联影医疗科技有限公司 | 体数据处理方法、装置、图像处理工作站和可读存储介质 |
CN110313943A (zh) * | 2018-03-30 | 2019-10-11 | 佳能医疗***株式会社 | 医用诊断装置、医用图像处理装置以及图像处理方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4167162B2 (ja) * | 2003-10-14 | 2008-10-15 | アロカ株式会社 | 超音波診断装置 |
KR101563506B1 (ko) * | 2013-01-16 | 2015-10-27 | 삼성메디슨 주식회사 | 초음파 장치 및 초음파 장치의 정보 입력 방법 |
EP3142560B1 (en) * | 2014-05-12 | 2021-08-18 | Exact Imaging Inc. | Medical-imaging system and method thereof |
EP3517047B1 (en) * | 2016-09-26 | 2023-12-20 | FUJIFILM Corporation | Ultrasonic diagnostic apparatus and control method of ultrasonic diagnostic apparatus |
CN106725612B (zh) * | 2016-12-23 | 2020-10-23 | 深圳开立生物医疗科技股份有限公司 | 四维超声图像优化方法及*** |
CN109480900B (zh) * | 2018-11-01 | 2021-08-13 | 飞依诺科技(苏州)有限公司 | 4d超声探头扫查控制方法、装置及存储介质 |
-
2020
- 2020-02-28 CN CN202010129489.8A patent/CN111358492A/zh active Pending
- 2020-06-16 WO PCT/CN2020/096260 patent/WO2021169086A1/zh active Application Filing
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101102721A (zh) * | 2005-02-09 | 2008-01-09 | 株式会社日立医药 | 超声波诊断装置以及超声波摄像方法 |
JP2010088585A (ja) * | 2008-10-06 | 2010-04-22 | Toshiba Corp | 超音波診断装置及び超音波プローブ |
CN101634551A (zh) * | 2009-08-18 | 2010-01-27 | 清华大学深圳研究生院 | 一种检测表面粗糙度的方法及其*** |
CN102247171A (zh) * | 2010-03-30 | 2011-11-23 | 株式会社东芝 | 超声波诊断装置、超声波图像处理装置以及医用图像诊断装置 |
CN107440720A (zh) * | 2012-06-27 | 2017-12-08 | 东芝医疗***株式会社 | 超声波诊断装置以及图像数据的校正方法 |
CN103473805A (zh) * | 2013-09-18 | 2013-12-25 | 内蒙古科技大学 | 基于改进区域生长算法测量三维重建肝脏模型体积的方法 |
CN104574329A (zh) * | 2013-10-09 | 2015-04-29 | 深圳迈瑞生物医疗电子股份有限公司 | 超声融合成像方法、超声融合成像导航*** |
CN106327487A (zh) * | 2016-08-18 | 2017-01-11 | 苏州润心医疗科技有限公司 | 基于x射线冠脉造影图像的冠状动脉血流储备分数计算方法 |
CN110313943A (zh) * | 2018-03-30 | 2019-10-11 | 佳能医疗***株式会社 | 医用诊断装置、医用图像处理装置以及图像处理方法 |
CN109447974A (zh) * | 2018-10-31 | 2019-03-08 | 上海联影医疗科技有限公司 | 体数据处理方法、装置、图像处理工作站和可读存储介质 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112529804A (zh) * | 2020-12-14 | 2021-03-19 | 深圳开立生物医疗科技股份有限公司 | 一种图像去噪方法、装置、设备及计算机可读存储介质 |
Also Published As
Publication number | Publication date |
---|---|
WO2021169086A1 (zh) | 2021-09-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110570492A (zh) | 神经网络训练方法和设备、图像处理方法和设备以及介质 | |
CN110648337A (zh) | 髋关节分割方法、装置、电子设备和存储介质 | |
CN110998602A (zh) | 使用深度学习方法对3d牙颌面结构的分类和3d建模 | |
US20030174890A1 (en) | Image processing device and ultrasonic diagnostic device | |
KR20210005649A (ko) | 딥 러닝 기술을 이용하여 금속 영향을 받은 x-선 데이터의 복셀 표현을 자동 수정 | |
CN113888689A (zh) | 图像渲染模型训练、图像渲染方法及装置 | |
CN111598989B (zh) | 一种图像渲染参数设置方法、装置、电子设备及存储介质 | |
CN110176010B (zh) | 一种图像检测方法、装置、设备及存储介质 | |
JP2008511366A (ja) | 距離座標を使った特徴重み付け医療オブジェクト輪郭検出 | |
CN112598649B (zh) | 基于生成对抗网络的2d/3d脊椎ct非刚性配准方法 | |
CN111524200B (zh) | 在投影图像中分割金属对象的方法、设备和介质 | |
CN114493995A (zh) | 图像渲染模型训练、图像渲染方法及装置 | |
CN111242926A (zh) | 病灶检测方法、装置及电子设备 | |
CN113554742A (zh) | 一种三维图像的重建方法、装置、设备及存储介质 | |
JP2023515367A (ja) | モデルへの入力インスタンスの分布外検出 | |
CN113689337B (zh) | 一种基于生成对抗网络的超声图像超分辨率重建方法及*** | |
CN111358492A (zh) | 一种四维造影图像生成方法、装置、设备及存储介质 | |
US7068837B2 (en) | Image processing method | |
CN116363038B (zh) | 超声图像融合方法、装置、计算机设备及存储介质 | |
CN111681297B (zh) | 图像重建方法、计算机设备和存储介质 | |
CN113658284A (zh) | 用于训练结节检测***的来自ct图像的x射线图像合成 | |
Johari et al. | Metal artifact suppression in dental cone beam computed tomography images using image processing techniques | |
CN115439423B (zh) | 一种基于ct图像的识别方法、装置、设备及存储介质 | |
JP4679095B2 (ja) | 画像処理装置、画像処理方法、及びプログラム | |
CN113096238B (zh) | 一种x射线图模拟方法、装置、电子设备及存储介质 |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200703 |