CN116549011A - 一种多普勒超声信号处理方法、系统和存储介质 - Google Patents
一种多普勒超声信号处理方法、系统和存储介质 Download PDFInfo
- Publication number
- CN116549011A CN116549011A CN202310708389.4A CN202310708389A CN116549011A CN 116549011 A CN116549011 A CN 116549011A CN 202310708389 A CN202310708389 A CN 202310708389A CN 116549011 A CN116549011 A CN 116549011A
- Authority
- CN
- China
- Prior art keywords
- line
- spectral lines
- spectral
- interpolation
- interpolated
- 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
- 238000003672 processing method Methods 0.000 title claims abstract description 9
- 230000003595 spectral effect Effects 0.000 claims abstract description 279
- 238000002604 ultrasonography Methods 0.000 claims abstract description 125
- 238000000034 method Methods 0.000 claims abstract description 34
- 230000009466 transformation Effects 0.000 claims description 38
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000007423 decrease Effects 0.000 claims description 3
- 239000000523 sample Substances 0.000 claims description 2
- 230000005236 sound signal Effects 0.000 description 28
- 238000001228 spectrum Methods 0.000 description 18
- 238000010586 diagram Methods 0.000 description 14
- 230000008569 process Effects 0.000 description 13
- 230000004048 modification Effects 0.000 description 9
- 238000012986 modification Methods 0.000 description 9
- 230000006870 function Effects 0.000 description 6
- 239000011800 void material Substances 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 4
- 230000006872 improvement Effects 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 230000000644 propagated effect Effects 0.000 description 3
- 239000008280 blood Substances 0.000 description 2
- 210000004369 blood Anatomy 0.000 description 2
- 230000000052 comparative effect Effects 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 238000012886 linear function Methods 0.000 description 2
- 238000012887 quadratic function Methods 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 230000001131 transforming effect Effects 0.000 description 2
- 230000007306 turnover Effects 0.000 description 2
- 241000579895 Chlorostilbon Species 0.000 description 1
- 206010028980 Neoplasm Diseases 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000000712 assembly Effects 0.000 description 1
- 238000000429 assembly Methods 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 238000013501 data transformation Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 239000010976 emerald Substances 0.000 description 1
- 229910052876 emerald Inorganic materials 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000010977 jade Substances 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 239000003607 modifier Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- ZLIBICFPKPWGIZ-UHFFFAOYSA-N pyrimethanil Chemical compound CC1=CC(C)=NC(NC=2C=CC=CC=2)=N1 ZLIBICFPKPWGIZ-UHFFFAOYSA-N 0.000 description 1
- 239000010979 ruby Substances 0.000 description 1
- 229910001750 ruby Inorganic materials 0.000 description 1
- 230000009469 supplementation Effects 0.000 description 1
- 230000001502 supplementing effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 238000012285 ultrasound imaging Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
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/488—Diagnostic techniques involving Doppler signals
-
- 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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Medical Informatics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biophysics (AREA)
- Physics & Mathematics (AREA)
- Heart & Thoracic Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本说明书公开了一种多普勒超声信号处理方法、系统和存储介质,所述方法包括:获取至少两段多普勒超声数据;对至少两段多普勒超声数据进行短时傅里叶变换,得到对应的至少两组频谱谱线;选择至少两组频谱谱线的至少部分谱线作为参考谱线,基于参考谱线和对应的插值参数,对至少两组频谱谱线进行插值处理。
Description
技术领域
本说明书涉及超声成像技术领域,特别涉及一种多普勒超声信号处理方法、系统和存储介质。
背景技术
多普勒超声是医用超声模式的一种,可以获取目标对象中血液的运动信息,并通过频谱图像和/或音频信号的方式呈现。在多模式工作时(例如,多普勒超声模式、B模式、color模式),多模式间通常采用间隔工作的方式,相应获取的多普勒超声数据存在缺失部分。因此,希望可以提供一种多普勒超声信号处理方法,通过插值处理,准确完整地呈现频谱图像和/或音频信号。
发明内容
本说明书的一个方面提供一种多普勒超声信号处理方法,所述方法包括:获取至少两段多普勒超声数据;对至少两段多普勒超声数据进行短时傅里叶变换,得到对应的至少两组频谱谱线;选择至少两组频谱谱线的至少部分谱线作为参考谱线,基于参考谱线和对应的插值参数,对至少两组频谱谱线进行插值处理。
本说明书的一个方面提供一种多普勒超声信号处理系统,所述系统包括:获取模块,用于获取至少两段多普勒超声数据;变换模块,用于对至少两段多普勒超声数据进行短时傅里叶变换,得到对应的至少两组频谱谱线;插值模块,用于选择至少两组频谱谱线的至少部分谱线作为参考谱线,基于参考谱线和对应的插值参数,对至少两组频谱谱线进行插值处理。
本说明书的一个方面提供一种多普勒超声信号处理系统,所述系统包括:超声设备,用于获取至少两段多普勒超声数据;至少一个处理器,用于对至少两段多普勒超声数据进行短时傅里叶变换,得到对应的至少两组频谱谱线;以及选择至少两组频谱谱线的至少部分谱线作为参考谱线,基于参考谱线和对应的插值参数,对至少两组频谱谱线进行插值处理。
本说明书的另一个方面提供一种计算机可读存储介质,所述存储介质存储计算机指令,当计算机读取存储介质中的计算机指令后,计算机执行多普勒超声信号处理方法。
本说明书的一些实施例对多普勒超声数据进行短时傅里叶变换得到频谱谱线,并基于频谱谱线获取参考谱线,根据参考谱线在频域对多普勒超声数据进行插值,可以减少频谱模糊的情况,从而可以准确地更具针对性地显示完整的频谱图像和播放不间断的音频信号。
本说明书的一些实施例针对不同的处理目标,确定不同的傅里叶变换参数,例如图像生成对应的傅里叶变换参数(例如,第一预设重叠长度)不同于音频生成对应的傅里叶变换参数(例如,第二预设重叠长度),从而可以针对频谱图像和音频信号的不同需求,获取不同数量的频谱谱线,一方面可以同时保证频谱图像和音频信号的质量,另一方面可以提高插值效率。
本说明书的一些实施例对小于预设变换长度的待变换多普勒超声数据进行数据填充,使得每段多普勒超声数据中的所有PW数据点都可以用于生成均匀的频谱谱线,从而减少频谱谱线的信息缺失。
附图说明
本说明书将以示例性实施例的方式进一步描述,这些示例性实施例将通过附图进行详细描述。这些实施例并非限制性的,在这些实施例中,相同的编号表示相同的结构,其中:
图1是根据本说明书的一些实施例所示的多普勒超声信号处理系统的应用场景示意图;
图2是根据本说明书的一些实施例所示的多普勒超声信号处理系统的示例性模块图;
图3是根据本说明书的一些实施例所示的多普勒超声信号处理方法的示例性流程图;
图4是根据本说明书的一些实施例所示的多普勒超声数据的示例性示意图;
图5是根据本说明书的一些实施例所示的短时傅里叶变换的示例性示意图;
图6是根据本说明书的一些实施例所示的对多普勒超声数据进行数据填充的示例性流程图;
图7是根据本说明书的一些实施例所示的插值处理的示例性示意图;
图8是根据本说明书的一些实施例所示的图像生成对应的插值处理过程的示例性示意图;
图9是根据本说明书的一些实施例所示的图像生成对应的频谱图像的示例性示意图;
图10是根据本说明书的一些实施例所示的音频生成对应的插值处理过程的示例性示意图;
图11是根据本说明书的一些实施例所示的经插值处理后生成的音频信号与未经插值处理而直接生成的音频信号的示例性对比示意图。
具体实施方式
为了更清楚地说明本说明书实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单的介绍。显而易见地,下面描述中的附图仅仅是本说明书的一些示例或实施例,对于本领域的普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图将本说明书应用于其它类似情景。除非从语言环境中显而易见或另做说明,图中相同标号代表相同结构或操作。
应当理解,本说明书中所使用的“系统”、“装置”、“单元”和/或“模组”是用于区分不同级别的不同组件、元件、部件、部分或装配的一种方法。然而,如果其他词语可实现相同的目的,则可通过其他表达来替换所述词语。
如本说明书和权利要求书中所示,除非上下文明确提示例外情形,“一”、“一个”、“一种”和/或“该”等词并非特指单数,也可包括复数。一般说来,术语“包括”与“包含”仅提示包括已明确标识的步骤和元素,而这些步骤和元素不构成一个排它性的罗列,方法或者设备也可能包含其它的步骤或元素。
本说明书中使用了流程图用来说明根据本说明书的实施例的系统所执行的操作。应当理解的是,前面或后面操作不一定按照顺序来精确地执行。相反,可以按照倒序或同时处理各个步骤。同时,也可以将其他操作添加到这些过程中,或从这些过程移除某一步或数步操作。
图1是根据本说明书的一些实施例所示的多普勒超声信号处理系统的应用场景示意图。在一些实施例中,如图1所示,多普勒超声信号处理系统100可以包括超声设备110、处理设备120、终端设备130、网络140和存储设备150。
在一些实施例中,多普勒超声信号处理系统100的组件可以以一种或多种各种方式连接。仅作为示例,如图1所示,超声设备110可以通过网络140连接到处理设备120。又例如,超声设备110可以直接连接到处理设备120(如连接超声设备110和处理设备120的虚线双向箭头所示)。作为进一步的示例,存储设备150可以直接或通过网络140连接到处理设备120。作为进一步的示例,终端设备130可以直接(如连接终端设备130和处理设备120的虚线双向箭头所示)与处理设备120连接和/或通过网络140与处理设备120连接。
超声设备110可以发射超声波至目标对象,以采集相应的超声数据。在一些实施例中,超声设备110可以通过超声探头向目标对象发射超声波,并从该目标对象获取反射超声波(例如,多普勒超声数据)。在一些实施例中,目标对象可以是人体、器官、机体、物体、损伤部位、肿瘤等或其任意组合。例如,目标对象在可以是心脏的一个或多个病变组织。在一些实施例中,超声设备110可以在多种模式(例如,B模式、color模式、脉冲多普勒模式)下工作,并采集相应模式下的超声数据。例如,以频谱多普勒模式为例,超声设备110可以发射超声波至目标对象,并采集目标对象反射的反射波。相应得到的多普勒超声数据可以体现目标对象中血液的运动信息(例如,血液流动速度和/或方向)。在一些实施例中,超声设备110可以在不同的时间段采集(例如,间隔采集)不同模式下的超声数据。为方便描述,后文将脉冲多普勒模式简称为PW模式,该模式下采集的超声数据称之为PW数据。可以理解,频谱多普勒模式还可以包括其他模式(例如,CW模式),相应还可以采集其他模式下的超声数据(例如,CW数据)。
处理设备120可以处理从超声设备110、终端设备130和/或存储设备150获取的数据和/或信息。例如,处理设备120可以对超声数据(例如,多普勒超声数据)执行短时傅里叶变换,以确定相应的频谱谱线。又例如,处理设备120可以对频谱谱线进行插值处理,以确定最终的频谱图像。再例如,处理设备120可以对频谱谱线进行短时傅里叶逆变换,以确定相应的音频信号。在一些实施例中,处理设备120可以包括中央处理单元(CPU)、数字信号处理器(DSP)等和/或其任意组合。在一些实施例中,处理设备120可以包括计算机、用户控制台、单个服务器或服务器组等。服务器组可以是集中式或分布式的。在一些实施例中,处理设备120可以是本地的或远程的。在一些实施例中,处理设备120可以在云平台上实现。仅作为示例,云平台可以包括私有云、公共云、混合云、社区云、分布式云、云间、多云等或其任意组合。在一些实施例中,处理设备120或处理设备120的一部分可以集成到超声设备110中。
终端设备130可以展示超声数据和/或超声数据的处理结果。在一些实施例中,终端设备130可以包含显示设备和/或音频设备。显示设备可以显示频谱图像,音频设备可以播放音频信号。在一些实施例中,终端设备130还可以实现用户与多普勒超声信号处理系统100的交互。在一些实施例中,终端设备130可以包括移动设备131、平板计算机132、笔记本计算机133等或其任意组合。在一些实施例中,终端设备130可以是处理设备120的一部分。
网络140可以包括促进多普勒超声信号处理系统100的信息和/或数据交换的任何合适的网络。在一些实施例中,一个或以上多普勒超声信号处理系统100的组件(例如,超声设备110、处理设备120、存储设备150、终端设备130)可以通过网络140与多普勒超声信号处理系统100的一个或以上其他组件通信信息和/或数据。例如,处理设备120可以经由网络从超声设备110接收多普勒超声数据。又例如,终端设备130可以经由网络140从处理设备120中获取频谱图像和/或音频信号。在一些实施例中,网络140可以是和/或包括公共网络、私有网络、广域网(WAN))、有线网络、无线网络、蜂窝网络、帧中继网络、虚拟专用网、卫星网络、电话网络、路由器、集线器、交换机等或其任意组合。
存储设备150可以存储数据、指令和/或任何其他信息。在一些实施例中,存储设备150可以存储从超声设备110、终端设备130和/或处理设备120获取的数据。在一些实施例中,存储设备150可以存储数据和/或指令,处理设备120可以执行或使用该数据和指令来执行本说明书中描述的示例性方法/系统。在一些实施例中,存储设备150可包括大容量存储器、可移动存储器、易失性读写内存、只读内存(ROM)等或其任意组合。在一些实施例中,存储设备150可在云平台上执行。在一些实施例中,存储设备150可以是处理设备120的一部分。
应该注意的是,上述描述仅出于说明性目的而提供,并不旨在限制本说明书的范围。对于本领域普通技术人员而言,在本说明书内容的指导下,可做出多种变化和修改。可以以各种方式组合本说明书描述的示例性实施例的特征、结构、方法和其他特征,以获取另外的和/或替代的示例性实施例。然而,这些变化与修改不会背离本说明书的范围。
图2是根据本说明书的一些实施例所示的多普勒超声信号处理系统的示例性模块图。
在一些实施例中,多普勒超声信号处理系统200可以包括获取模块210、变换模块220和/或插值模块230。
获取模块210可以用于获取至少两段多普勒超声数据。
变换模块220可以用于对至少两段多普勒超声数据进行短时傅里叶变换,得到对应的至少两组频谱谱线。在一些实施例中,变换模块220可以基于预设变换长度和预设重叠长度,执行短时傅里叶变换。在一些实施例中,若短时傅里叶变换中的待变换多普勒超声数据的长度小于预设变换长度,变换模块220可以对待变换多普勒超声数据进行数据填充。在一些实施例中,预设变换长度和预设重叠长度可以基于处理目标确定。在一些实施例中,处理目标可以包括图像生成和/或音频生成。在一些实施例中,预设重叠长度可以包括对应图像生成的第一预设重叠长度和/或所述音频生成的第二预设重叠长度。在一些实施例中,第一预设重叠长度可以大于第二预设重叠长度。
插值模块230可以用于选择至少两组频谱谱线的至少部分谱线作为参考谱线,基于参考谱线和对应的插值参数,对至少两组频谱谱线进行插值处理。在一些实施例中,插值参数可以与参考谱线和待插值谱线间的距离负相关。在一些实施例中,同一参考谱线,对应的插值参数可以随参考谱线和待插值谱线之间距离的增大而减小。在一些实施例中,一根待插值谱线可以基于至少一对参考谱线插值确定。在一些实施例中,对于同一待插值谱线,参考谱线对应的插值参数之和可以为预设常数,插值参数可以与对应的参考谱线和待插值谱线之间的距离负相关。在一些实施例中,处理目标可以为音频生成。在一些实施例中,待插值谱线可以包括第一根待插值谱线、最后一根待插值谱线和中间待插值谱线,其中,第一根待插值谱线和最后一根待插值谱线可以分别基于其对应的相邻的参考谱线确定;对于中间待插值谱线,插值参数可以与参考谱线和中间待插值谱线间的距离负相关。
图3是根据本说明书的一些实施例所示的多普勒超声信号处理方法的示例性流程图。在一些实施例中,流程300可以由多普勒超声信号处理系统100(例如,处理设备120)或多普勒超声信号处理系统200执行。例如,流程300可以以程序或指令的形式存储在存储设备(例如,存储设备150、系统的存储单元)中,当处理设备120或多普勒超声信号处理系统200执行指令时,可以实现流程300。下面呈现的流程300的操作示意图是说明性的。在一些实施例中,可以利用一个或以上未描述的附加操作和/或未讨论的一个或以上操作来完成该过程。另外,图3中示出的和下面描述的流程300的操作的顺序并非限制性的。
步骤310,获取至少两段多普勒超声数据。具体地,步骤310可以由获取模块210执行。
结合图1所述,超声设备110可以间隔地在多种模式(例如,B模式、C模式、PW模式)下工作,并采集相应模式下的超声数据。例如,图4是根据本说明书的一些实施例所示的多普勒超声数据的示例性示意图,如图4所示,超声设备110可以间隔地在PW模式和其他模式(例如,B模式或color模式)下工作,分别采集PW数据和其他超声数据(例如,B模式的超声数据或color模式的超声数据)。也就是说,间隔采集的PW数据间存在空缺部分,导致PW数据不连续。相应地,本说明书实施例所述的方法可以获取至少两段不连续的PW数据,对其进行插值处理,实现对空缺部分的插值补充。
如图4所示,获取模块210获取的至少两段多普勒超声数据可以简称为PW数据1、PW数据2、……、PW数据n,其中,每一段PW数据可以包括m个数据点。
步骤320,对至少两段多普勒超声数据进行短时傅里叶变换,得到对应的至少两组频谱谱线。具体地,步骤320可以由变换模块220执行。
在一些实施例中,变换模块220可以基于预设变换长度和预设重叠长度,执行短时傅里叶变换。
预设变换长度可以是进行一次短时傅里叶变换的数据长度。也就是说,预设变换长度可以理解为用于生成一根频谱谱线的PW数据点的数量。在一些实施例中,预设变换长度可以设置为2的n次方,例如,32、64、128等。在一些实施例中,预设变换长度可以基于每段多普勒超声数据包含的PW数据点的数量确定。每段多普勒超声数据包含的PW数据点的数量越多,预设变换长度可以越长。仅作为示例,假设每段多普勒超声数据包含的PW数据点的数量小于300,则预设变换长度可以为32;假设每段多普勒超声数据包含的PW数据点的数量大于等于300且小于600,则预设变换长度可以为64;假设每段多普勒超声数据包含的PW数据点的数量大于等于600,则预设变换长度可以为128。
预设重叠长度可以是相邻两次短时傅里叶变换对应的数据之间的重叠长度。也就是说,预设重叠长度可以理解为用于生成相邻两根频谱谱线的PW数据点的重叠数量。可以理解,预设重叠长度越长,则相邻两根频谱谱线的间隔越小,每段多普勒超声数据变换得到的频谱谱线数量越多。
短时傅里叶变换可以将时域的多普勒超声波数据转换为频域的频谱数据。仅作为示例,短时傅里叶变换可以对预设变换长度内的PW数据点进行傅里叶变换,然后基于预设重叠长度获取下一个预设变换长度内的PW数据点,重复上述操作,以获取每次短时傅里叶变换后的频谱谱线。
在一些实施例中,处理目标不同,所需执行的短时傅里叶变换的参数相应不同。处理目标可以指对多普勒超声数据进行数据处理的目的或目标。在一些实施例中,处理目标可以包括图像生成、音频生成等或其任意组合。图像生成可以理解为数据处理的目的是为了获取频谱图像(例如,在频率方向上连续的多普勒超声数据)。例如,可以将时域的多普勒超声数据变换为频域的频谱谱线,从而获取频谱图像。音频生成可以理解为数据处理的目的是为了获取音频信号(例如,在时间方向上连续的多普勒超声数据)。例如,可以将时域的多普勒超声数据变换为频域的频谱谱线,并对频谱谱线逆变换为时域的音频信号。关于频谱图像和音频信号的更多描述可以参见步骤330的相关描述,在此不再赘述。在一些实施例中,预设变化长度和预设重叠长度可以基于处理目标确定。
在一些实施例中,图像生成(即处理目标为图像生成)对应的预设变换长度和音频生成(即处理目标为音频生成)对应的预设变换长度可以不同。
在一些实施例中,预设重叠长度可以包括图像生成(即处理目标为图像生成)对应的第一预设重叠长度和/或音频生成(即处理目标为音频生成)对应的第二预设重叠长度。在一些实施例中,第一预设重叠长度与第二预设重叠长度不同。在一些实施例中,第一预设重叠长度可以大于第二预设重叠长度。也就是说,用于生成频谱图像的频谱谱线数量大于用于生成音频信号的频谱谱线数量。例如,第一预设重叠长度为32,第二预设重叠长度为16。
在本说明书的一些实施例中,图像生成对应的傅里叶变换参数(例如,第一预设重叠长度)不同于音频生成对应的傅里叶变换参数(例如,第二预设重叠长度),可以针对频谱图像和音频信号的不同需求,获取不同数量的频谱谱线,一方面可以同时保证频谱图像和音频信号的质量,另一方面可以提高插值效率。
在一些实施例中,对于每段PW数据,变换模块220可以取预设变换长度的PW数据点进行短时傅里叶变换,生成一根频谱谱线;然后取下一相同预设变换长度的PW数据点进行短时傅里叶变换,生成下一根频谱谱线,相邻两段预设变换长度的PW数据的重叠长度为预设重叠长度;以此类推,直到取完所有的PW数据点。
例如,图5是根据本说明书的一些实施例所示的短时傅里叶变换的示例性示意图,如图5所示,以一段PW数据为例,该段PW数据包含592个PW数据点,预设变换长度(如图中实线箭头所示)为64,预设重叠长度(如图中虚线箭头所示)为16,则变换模块220可以基于每段PW数据中的第1~64个PW数据点生成第1根频谱谱线,第49~112个PW数据点生成第2根频谱谱线,……,第529~592个PW数据点生成第12根频谱谱线。
在一些实施例中,若待变换多普勒超声数据的长度(例如,待变换多普勒超声数据中PW数据点的数量)小于预设变换长度,变换模块220可以对待变换多普勒超声数据进行数据填充。
仅作为示例,图6是根据本说明书的一些实施例所示的对多普勒超声数据进行数据填充的示例性流程图,如图6所示,以一段PW数据为例,该段PW数据包含592个PW数据点,预设变换长度(如图中实线箭头所示)为64,预设重叠长度(如图中虚线箭头所示)为32,变换模块220基于第1~64个PW数据点生成第1根频谱谱线后,待变换多普勒超声数据为剩余PW数据点及预设重叠长度对应的PW数据点,即待变换多普勒超声数据为第33~592个PW数据点。
在一些实施例中,若当前待变换多普勒超声数据的长度大于预设变换长度,则变换模块220继续从待变换多普勒超声数据中取预设变换长度的PW数据点执行短时傅里叶变换,生成下一根频谱谱线。例如,继续上述示例,待变换多普勒超声数据“第33~592个PW数据点”的长度为560,大于预设变换长度为64,则变换模块220可以从待变换多普勒超声数据“第33~592个PW数据点”中继续取64个PW数据点“第33~96个PW数据点”执行短时傅里叶变换,获取对应的第2根频谱谱线。
在一些实施例中,若待变换多普勒超声数据的长度小于预设变换长度,则变换模块220可以对待变换多普勒超声数据进行数据填充。例如,继续上述示例,变换模块220将第1~576个PW数据点变换为17根频谱谱线后,待变换多普勒超声数据为“第545~592个PW数据点”,其长度为48,小于预设变换长度64,则变换模块220可以对待变换多普勒超声数据进行数据填充。在一些实施例中,变换模块220可以在待变换多普勒超声数据的任意位置填充0数据,将待变换多普勒超声数据的长度填充为预设变换长度。仅作为示例,如图6所示,变换模块220可以用16个0数据在“第592个PW数据点”后填充第593~608个PW数据点,使得填充后的待变换多普勒超声数据“第545~608个PW数据点”的长度为64。
进一步地,变换模块220可以对填充后的待变换多普勒超声数据执行短时傅里叶变换,生成对应的频谱谱线。例如,继续上述示例,变换模块220可以对“第545~608个PW数据点”执行短时傅里叶变换,获取图像生成对应的第18根频谱谱线。
在本说明书的一些实施例中,对小于预设变换长度的待变换多普勒超声数据进行数据填充,使得每段多普勒超声数据中的所有PW数据点都可以用于生成均匀的频谱谱线,从而减少频谱谱线的信息缺失。
在一些情况下,多普勒超声数据可能包含目标对象的运动组织的运动信息(例如,患者移动引起的位移运动信息)。为了减少运动组织的运动信息对插值效果的影响,在一些实施例中,在对多勒普超声数据执行短时傅里叶变换前,变换模块220可以先对每段多勒普超声数据进行预处理。具体地,变换模块220可以对每段多勒普超声数据进行高通滤波,从而过滤掉多勒普超声数据中运动组织低频的运动信息。
仅作为示例,变换模块220可以对巴特沃斯滤波器进行投影初始化,获取m×m的二维方阵,其中,m为每段多普勒超声数据中包含的PW数据点。进一步地,变换模块220可以对m×m的二维方阵和m×1的每段PW数据进行矩阵乘法运算,获取过滤后的PW数据。
步骤330,选择至少两组频谱谱线的至少部分谱线作为参考谱线,基于参考谱线和对应的插值参数,对至少两组频谱谱线进行插值处理。具体地,步骤330可以由插值模块230执行。
结合前文所述,经过短时傅里叶变换后,可以得到至少两段多普勒超声数据分别对应的至少两组频谱谱线。在一些实施例中,每组频谱谱线可以包含多根频谱谱线。例如,图7是插值处理的示例性示意图,如图7所示,n段多普勒超声数据“PW数据1”、“PW数据2”、…“PW数据n”可以分别对应“第1组频谱谱线”、“第2组频谱谱线”、…“第n组频谱谱线”,每组频谱谱线可以包含i根频谱谱线。
结合步骤310所述,相邻两段多普勒超声数据之间存在空缺部分。相应地,经过短时傅里叶变换后,相邻两组频谱谱线之间也存在空缺部分。在一些实施例中,插值模块230可以选择至少部分频谱谱线作为参考谱线,对空缺部分进行插值补充。在一些实施例中,在对空缺部分进行插值补充时,每段空缺部分可以对应一个插值组,其中包括一根或多根经过插值确定的谱线。例如,如图7所示,“空缺1”、“空缺2”、…“空缺(n-1)”可以分别对应“第1个插值组”、“第2个插值组”、…“第(n-1)个插值组”。
在一些实施例中,插值模块230可以基于每段空缺部分对应的PW数据点数量和预设重叠长度,确定每个插值组包含的待插值谱线的数量。以一段空缺部分为例,该段空缺部分对应的PW数据点数量可以理解为在该段空缺部分对应的时间内超声设备110能够采集的PW数据点的数量。仅作为示例,结合图4,假设一段PW数据的数据长度(或数据点数量)为L,空缺部分的数据长度(或数据点数量)为D,将一段PW数据和随后的一段空缺部分作为一个数据长度为L+D的数据周期。假设预设变换长度为V,预设重叠长度为R,则一个数据周期内可以生成1+(L+D-V)/R根谱线,实际生成1+(L-V)/R根谱线,因此,每个插值组包含的待插值谱线的数量为D/R。具体例如,假设空缺部分的数据长度为384,若第一预设重叠长度为16,则对应的插值组包含的待插值谱线的数量为384/16=24根;若第二预设重叠长度为32,对应的插值组包含的待插值谱线的数量为384/32=12根。
在一些实施例中,对于一个插值组,插值模块230可以从该插值组前后相邻的两组频谱谱线中选择参考谱线。相应地,在前后相邻的两组频谱谱线中选择的参考谱线可以分别称之为前参考谱线组和后参考谱线组。例如,图8是根据本说明书的一些实施例所示的图像生成对应的插值处理过程的示例性示意图,如图8所示,以第1个插值组为例,插值模块230可以从第1个插值组之前相邻的“第1组频谱谱线”和之后相邻的“第2组频谱谱线”中分别选择参考谱线以分别形成第1组参考谱线和第2组参考谱线。在一些实施例中,插值模块230也可以进一步从插值组前后的其他组频谱谱线中选择参考谱线,在此不再赘述。
在一些实施例中,插值模块230可以将与插值组相邻的多条频谱谱线作为参考谱线。例如,将前参考谱线组中与插值组相邻的后几条频谱谱线和后参考谱线组中与插值组相邻的前几条频谱谱线作为参考谱线。在一些实施例中,插值模块230可以基于频谱谱线与插值组的距离,选择对应的参考谱线。在一些实施例中,插值模块230可以选择与插值组的距离在预设距离范围内的频谱谱线作为参考谱线。例如,如图8所示,插值模块230可以将第1组频谱谱线中最后4根频谱谱线(即参考谱线组b)和第2组频谱谱线中前4根频谱谱线(即参考谱线组a)作为第1个插值组的参考谱线。
在一些实施例中,插值模块230可以从频谱谱线中随机选择参考谱线。例如,仍以第1个插值组为例,插值模块230可以随机从第1组频谱谱线和第2组频谱谱线中选择参考谱线。
在一些实施例中,插值模块230可以按预设规则从频谱谱线中选择参考谱线。例如,仍以第1个插值组为例,插值模块230可以将第1组频谱谱线和第2组频谱谱线中的第1根、第4根、……、第7根频谱谱线作为第1个插值组的参考谱线。
在一些实施例中,插值模块230可以基于成像效率和/或成像质量,确定图像生成的插值组对应的参考谱线的数量。仅作为示例,插值模块230可以分别将插值组前后相邻的两组频谱谱线中2/3的频谱谱线作为对应的参考谱线,以获取较高的成像质量。作为又一示例,插值模块230可以分别将插值组前后相邻的两组频谱谱线中1/3的频谱谱线作为对应的参考谱线。例如,假设每组频谱谱线包含12根频谱谱线,插值模块230可以将第1组频谱谱线中的4根频谱谱线和第2组频谱谱线中的4根频谱谱线作为第1个插值组的参考谱线。与图像生成类似地,在一些实施例中,插值模块230可以基于音频生成效率和/或音频质量,确定音频生成的插值组对应的参考谱线的数量。
在一些实施例中,以一个插值组为例,插值模块230可以基于参考谱线和对应的插值参数,通过下述公式(1)确定待插值谱线:
其中,ji表示插值组中的第i根待插值谱线,c1表示该插值组对应的前参考谱线组中的参考谱线的数量,c2表示该插值组对应的后参考谱线组中的参考谱线数量,lb,k1表示前参考谱线组中的第k1根参考谱线,la,k2表示后参考谱线组中的第k2根参考谱线,表示参考谱线lb,k1对应的插值参数,/>表示参考谱线la,k2对应的插值参数,插值参数可以体现待插值谱线对应的参考谱线的权重。
在一些实施例中,插值参数可以与参考谱线和待插值谱线间的距离负相关。以一根待插值谱线和其对应的一个参考谱线组(例如,前参考谱线组)为例,插值模块230可以获取待插值谱线和每根参考谱线之间的距离,并基于待插值谱线和每根参考谱线之间的距离,通过预设函数确定每根参考谱线对应的插值参数。在一些实施例中,预设函数可以包括但不限于线性函数、二次函数、贝塞尔函数等。
在一些实施例中,对于同一参考谱线,其对应的插值参数随参考谱线和待插值谱线之间距离的增大而减小。结合前文,一个插值组中可以包括多根待插值谱线,不同的待插值谱线可以对应相同的参考谱线。也就是说,同一根参考谱线可以应用于不同的待插值谱线的确定。相应地,对于同一根参考谱线,其对应的插值参数为变量,随着待插值谱线的不同而不同。该参考谱线与待插值谱线间的距离越近,则对应的插值参数越大;反之,则对应的插值参数越小。
在本说明书的一些实施例中,插值参数与参考谱线和待插值谱线间的距离负相关,使得待插值谱线与其距离越近的参考谱线越相似,从而使得插值后的频谱谱线更加平滑。
在一些实施例中,一根待插值谱线可以基于至少一对参考谱线插值确定。由前述可知,对于一个插值组,插值模块230可以从该插值组前后相邻的两组频谱谱线中选择参考谱线。在前后相邻的两组频谱谱线中选择的参考谱线可以分别称之为前参考谱线组和后参考谱线组。相应地,前参考谱线组和后参考谱线组的参考谱线数量相同,且在数量序号上基于该插值组对称。在一些实施例中,插值模块230可以将插值组前后参考谱线组中基于该插值组对称的两条参考谱线作为一对参考谱线。例如,以前参考谱线组b和后参考谱线组a为例,前参考谱线组b中的第1根参考谱线可以和后参考谱线组a中的倒数第1根参考谱线作为一对参考谱线;前参考谱线组b中的第2根参考谱线可以和后参考谱线组a中的倒数第2根参考谱线作为一对参考谱线;……以此类推。
仅作为示例,插值模块230可以基于2对参考谱线和对应的4个插值参数,通过下述公式(2)确定待插值谱线:
其中,ji表示插值组中的第i根待插值谱线,lb,1和la,1分别表示第1对参考谱线中前参考谱线组中的第1根参考谱线和后参考谱线组中的第1根参考谱线,lb,2和la,2分别表示第2对参考谱线中前参考谱线组中的第2根参考谱线和后参考谱线组中的第2根参考谱线,分别表示参考谱线lb,1、lb,2、la,1、la,2对应的插值参数。
在一些实施例中,对于同一待插值谱线,参考谱线对应的插值参数之和可以为预设常数,且插值参数与对应的参考谱线和待插值谱线之间的距离负相关。例如,预设常数可以是1~2之间的任意数值。仅作为示例,预设常数为1,则公式(1)中的插值参数作为又一示例,预设常数为2,公式(2)中的插值参数
在一些实施例中,插值模块230可以基于待插值谱线和参考谱线之间的距离,通过预设函数确定一对参考谱线分别对应的插值参数。在一些实施例中,预设函数可以包括但不限于线性函数、二次函数、贝塞尔函数等。
在本说明书的一些实施例中,待插值谱线基于至少一对参考谱线插值确定,使得待插值谱线可以同时与前后两部分的频谱谱线平滑连接,减少频谱模糊、出现竖直白线等现象。
在一些实施例中,经过插值处理后,插值模块230可以基于频谱谱线和待插值谱线,获取在频率方向上连续的完整频谱谱线。进一步地,如图9所示,插值模块230可以将图像生成对应的完整频谱谱线作为图像生成对应的频谱图像。
在一些实施例中,频谱图像的格式可以包括Joint Photographic Experts Group(JPEG)图像格式、Tagged Image File Format(TIFF)图像格式、Graphics InterchangeFormat(GIF)图像格式、Kodak Flash PiX(FPX)图像格式、Digital Imaging andCommunications in Medicine(DICOM)图像格式等。
在一些实施例中,结合前文所述,不同的处理目标(例如,图像生成、音频生成)可以对应不同的处理参数(例如,傅里叶变换参数)。相应地,不同的处理目标也可以对应不同的插值处理条件。图10是根据本说明书的一些实施例所示的音频生成对应的插值处理过程的示例性示意图。
在一些实施例中,若处理目标为音频生成,插值组可以包括第一根待插值谱线、最后一根待插值谱线和中间待插值谱线。例如,如图10所示,第一根待插值谱线为j1,最后一根待插值谱线为j3,中间待插值谱线为j2。
在一些实施例中,第一根待插值谱线和最后一根待插值谱线可以分别基于其对应的相邻的参考谱线确定。
在一些实施例中,插值模块230可以基于第一根待插值谱线(或最后一根待插值谱线)和其相邻参考谱线之间的距离,确定对应的插值参数,然后将插值参数和该相邻参考谱线的乘积确定为第一根待插值谱线(或最后一根待插值谱线)。关于确定插值参数的详细描述可以参见前文,在此不再赘述。作为示例,如图10所示,lb,6是第一根待插值谱线j1的相邻参考谱线,插值模块230可以基于j1和lb,6之间的距离确定对应的插值参数,并基于相邻参考谱线与插值参数的乘积确定第一根待插值谱线j1。作为又一示例,la,1是最后一根待插值谱线j3的相邻参考谱线,插值模块230可以基于j3和la,1之间的距离确定对应的插值参数,并基于相邻参考谱线与插值参数的乘积确定最后一根待插值谱线j1。
在一些实施例中,插值模块230可以将第一根待插值谱线(或最后一根待插值谱线)的相邻参考谱线直接复制为第一根待插值谱线(或最后一根待插值谱线)。例如,如图10所示,lb,6是第一根待插值谱线j1的相邻参考谱线,插值模块230可以直接复制相邻参考谱线作为第一根待插值谱线,即j1=lb,6;la,1是最后一根待插值谱线j3的相邻参考谱线,插值模块230可以直接复制相邻参考谱线作为最后一根待插值谱线,即j3=la,1。
在一些实施例中,插值模块230可以基于至少一根参考谱线对应的至少一个插值参数,对至少一根参考谱线进行加权求和,从而获取对应的中间待插值谱线。在一些实施例中,对于中间待插值谱线,插值参数与参考谱线和中间待插值谱线间的距离负相关。例如,插值模块230可以基于公式(1)获取对应的中间待插值谱线。关于获取中间待插值谱线的详细描述可以参见前文,在此不再赘述。
在本说明书的一些实施例中,音频生成对应的第一根待插值谱线和最后一根待插值谱线分别基于其对应的相邻参考谱线确定,可以减少音频生成的谱线计算量。
在一些实施例中,完成针对音频生成的频谱谱线插值处理后,可以获取完整的频谱谱线。进一步地,插值模块230可以对完整频谱谱线进行傅里叶逆变换,以生成相应的音频信号。
在一些实施例中,音频信号可以包括左声道信号和右声道信号。在一些实施例中,插值模块230可以基于音频生成对应的完整频谱谱线获取正频信号和负频信号,并分别对正频部分和负频部分做翻转操作,获取对应的翻转负频和翻转正频。进一步地,将正频部分和翻转正频拼接作为正频信号,将负频部分和翻转负频拼接作为负频信号。更进一步地,插值模块230可以分别对频域的正频信号和负频信号做短时傅里叶逆变换,分别获取时域的左声道信号和右声道信号。
图11是根据本说明书的一些实施例所示的经插值处理后生成的音频信号与未经插值处理而直接生成的音频信号的示例性对比示意图。由此可见,经过插值处理后的音频信号更加连续饱满。
一般来说,人耳对音频信号的相位变化不敏感,使得音频信号听起来不够饱满真实,因此,在一些实施例中,插值模块230可以在获取音频信号之前,在音频生成对应的完整频谱谱线中添加低频段的随机相位噪声,使得声音听起来更加饱满真实。
在本说明书的一些实施例中,对多普勒超声数据进行短时傅里叶变换得到频谱谱线,并基于频谱谱线获取参考谱线,根据参考谱线在频域对多普勒超声数据进行插值,可以减少频谱模糊的情况,从而准确地更具针对性地显示完整的频谱图像和播放不间断的音频信号。
本说明书实施例可能带来的有益效果包括但不限于:(1)对多普勒超声数据进行短时傅里叶变换得到频谱谱线,并基于频谱谱线获取参考谱线,根据参考谱线在频域对多普勒超声数据进行插值,可以减少频谱模糊的情况,从而可以准确地更具针对性地显示完整的频谱图像和播放不间断的音频信号;(2)针对不同的处理目标,确定不同的傅里叶变换参数,例如图像生成对应的傅里叶变换参数(例如,第一预设重叠长度)不同于音频生成对应的傅里叶变换参数(例如,第二预设重叠长度),从而可以针对频谱图像和音频信号的不同需求,获取不同数量的频谱谱线,一方面可以同时保证频谱图像和音频信号的质量,另一方面可以提高插值效率;(3)对小于预设变换长度的待变换多普勒超声数据进行数据填充,使得每段多普勒超声数据中的所有PW数据点都可以用于生成均匀的频谱谱线,从而减少频谱谱线的信息缺失;(4)插值参数与参考谱线和待插值谱线间的距离负相关,使得待插值谱线与其距离越近的参考谱线越相似,从而使得插值后的频谱谱线更加平滑;(5)待插值谱线基于至少一对参考谱线插值确定,使得待插值谱线可以同时与前后两部分的频谱谱线平滑连接,减少频谱模糊、出现竖直白线等现象。需要说明的是,不同实施例可能产生的有益效果不同,在不同的实施例里,可能产生的有益效果可以是以上任意一种或几种的组合,也可以是其他任何可能获取的有益效果。
上文已对基本概念做了描述,显然,对于本领域技术人员来说,上述详细披露仅仅作为示例,而并不构成对本说明书的限定。虽然此处并没有明确说明,本领域技术人员可能会对本说明书进行各种修改、改进和修正。该类修改、改进和修正在本说明书中被建议,所以该类修改、改进、修正仍属于本说明书示范实施例的精神和范围。
同时,本说明书使用了特定词语来描述本说明书的实施例。如“一个实施例”、“一实施例”、和/或“一些实施例”意指与本说明书至少一个实施例相关的某一特征、结构或特点。因此,应强调并注意的是,本说明书中在不同位置两次或多次提及的“一实施例”或“一个实施例”或“一个替代性实施例”并不一定是指同一实施例。此外,本说明书的一个或多个实施例中的某些特征、结构或特点可以进行适当的组合。
此外,本领域技术人员可以理解,本说明书的各方面可以通过若干具有可专利性的种类或情况进行说明和描述,包括任何新的和有用的工序、机器、产品或物质的组合,或对他们的任何新的和有用的改进。相应地,本说明书的各个方面可以完全由硬件执行、可以完全由软件(包括固件、常驻软件、微码等)执行、也可以由硬件和软件组合执行。以上硬件或软件均可被称为“数据块”、“模块”、“引擎”、“单元”、“组件”或“系统”。此外,本说明书的各方面可能表现为位于一个或多个计算机可读介质中的计算机产品,该产品包括计算机可读程序编码。
计算机存储介质可能包含一个内含有计算机程序编码的传播数据信号,例如在基带上或作为载波的一部分。该传播信号可能有多种表现形式,包括电磁形式、光形式等,或合适的组合形式。计算机存储介质可以是除计算机可读存储介质之外的任何计算机可读介质,该介质可以通过连接至一个指令执行系统、装置或设备以实现通讯、传播或传输供使用的程序。位于计算机存储介质上的程序编码可以通过任何合适的介质进行传播,包括无线电、电缆、光纤电缆、RF、或类似介质,或任何上述介质的组合。
本说明书各部分操作所需的计算机程序编码可以用任意一种或多种程序语言编写,包括面向对象编程语言如Java、Scala、Smalltalk、Eiffel、JADE、Emerald、C++、C#、VB.NET、Python等,常规程序化编程语言如C语言、Visual Basic、Fortran2003、Perl、COBOL2002、PHP、ABAP,动态编程语言如Python、Ruby和Groovy,或其他编程语言等。该程序编码可以完全在用户计算机上运行、或作为独立的软件包在用户计算机上运行、或部分在用户计算机上运行部分在远程计算机运行、或完全在远程计算机或处理设备上运行。在后种情况下,远程计算机可以通过任何网络形式与用户计算机连接,比如局域网(LAN)或广域网(WAN),或连接至外部计算机(例如通过因特网),或在云计算环境中,或作为服务使用如软件即服务(SaaS)。
此外,除非权利要求中明确说明,本说明书所述处理元素和序列的顺序、数字字母的使用、或其他名称的使用,并非用于限定本说明书流程和方法的顺序。尽管上述披露中通过各种示例讨论了一些目前认为有用的发明实施例,但应当理解的是,该类细节仅起到说明的目的,附加的权利要求并不仅限于披露的实施例,相反,权利要求旨在覆盖所有符合本说明书实施例实质和范围的修正和等价组合。例如,虽然以上所描述的系统组件可以通过硬件设备实现,但是也可以只通过软件的解决方案得以实现,如在现有的处理设备或移动设备上安装所描述的系统。
同理,应当注意的是,为了简化本说明书披露的表述,从而帮助对一个或多个发明实施例的理解,前文对本说明书实施例的描述中,有时会将多种特征归并至一个实施例、附图或对其的描述中。但是,这种披露方法并不意味着本说明书对象所需要的特征比权利要求中提及的特征多。实际上,实施例的特征要少于上述披露的单个实施例的全部特征。
一些实施例中使用了描述成分、属性数量的数字,应当理解的是,此类用于实施例描述的数字,在一些示例中使用了修饰词“大约”、“近似”或“大体上”来修饰。除非另外说明,“大约”、“近似”或“大体上”表明所述数字允许有±20%的变化。相应地,在一些实施例中,说明书和权利要求中使用的数值参数均为近似值,该近似值根据个别实施例所需特点可以发生改变。在一些实施例中,数值参数应考虑规定的有效数位并采用一般位数保留的方法。尽管本说明书一些实施例中用于确认其范围广度的数值域和参数为近似值,在具体实施例中,此类数值的设定在可行范围内尽可能精确。
针对本说明书引用的每个专利、专利申请、专利申请公开物和其他材料,如文章、书籍、说明书、出版物、文档等,特此将其全部内容并入本说明书作为参考。与本说明书内容不一致或产生冲突的申请历史文件除外,对本说明书权利要求最广范围有限制的文件(当前或之后附加于本说明书中的)也除外。需要说明的是,如果本说明书附属材料中的描述、定义、和/或术语的使用与本说明书所述内容有不一致或冲突的地方,以本说明书的描述、定义和/或术语的使用为准。
最后,应当理解的是,本说明书中所述实施例仅用以说明本说明书实施例的原则。其他的变形也可能属于本说明书的范围。因此,作为示例而非限制,本说明书实施例的替代配置可视为与本说明书的教导一致。相应地,本说明书的实施例不仅限于本说明书明确介绍和描述的实施例。
Claims (11)
1.一种多普勒超声信号处理方法,其特征在于,所述方法包括:
获取至少两段多普勒超声数据;
对所述至少两段多普勒超声数据进行短时傅里叶变换,得到对应的至少两组频谱谱线;
选择所述至少两组频谱谱线的至少部分谱线作为参考谱线,基于所述参考谱线和对应的插值参数,对所述至少两组频谱谱线进行插值处理。
2.如权利要求1所述的方法,其特征在于,所述对所述至少两段多普勒超声数据进行短时傅里叶变换,包括:
基于预设变换长度和预设重叠长度,执行短时傅里叶变换;
若所述短时傅里叶变换中的待变换多普勒超声数据的长度小于所述预设变换长度,对所述待变换多普勒超声数据进行数据填充。
3.如权利要求2所述的方法,其特征在于,
所述预设变换长度和所述预设重叠长度基于处理目标确定,所述处理目标包括图像生成和/或音频生成,
所述预设重叠长度包括对应所述图像生成的第一预设重叠长度和/或对应所述音频生成的第二预设重叠长度,
所述第一预设重叠长度大于所述第二预设重叠长度。
4.如权利要求1所述的方法,其特征在于,所述插值参数与所述参考谱线和待插值谱线间的距离负相关。
5.如权利要求1所述的方法,其特征在于,对于同一参考谱线,对应的所述插值参数随所述参考谱线和待插值谱线之间距离的增大而减小。
6.如权利要求1所述的方法,其特征在于,一根待插值谱线基于至少一对参考谱线插值确定。
7.如权利要求1所述的方法,其特征在于,对于同一待插值谱线,所述参考谱线对应的插值参数之和为预设常数,所述插值参数与对应的所述参考谱线和所述待插值谱线之间的距离负相关。
8.如权利要求3所述的方法,其特征在于,所述处理目标为音频生成,待插值谱线包括第一根待插值谱线、最后一根待插值谱线和中间待插值谱线,其中,
所述第一根待插值谱线和所述最后一根待插值谱线分别基于其对应的相邻的所述参考谱线确定;
对于中间待插值谱线,所述插值参数与所述参考谱线和所述中间待插值谱线间的距离负相关。
9.一种多普勒超声信号处理系统,其特征在于,所述系统包括:
获取模块,用于获取至少两段多普勒超声数据;
变换模块,用于对所述至少两段多普勒超声数据进行短时傅里叶变换,得到对应的至少两组频谱谱线;
插值模块,用于选择所述至少两组频谱谱线的至少部分谱线作为参考谱线,基于所述参考谱线和对应的插值参数,对所述至少两组频谱谱线进行插值处理。
10.一种多普勒超声信号处理系统,其特征在于,所述系统包括:
超声探头,用于获取至少两段多普勒超声数据;
至少一个处理器,用于对所述至少两段多普勒超声数据进行短时傅里叶变换,得到对应的至少两组频谱谱线;以及选择所述至少两组频谱谱线的至少部分谱线作为参考谱线,基于所述参考谱线和对应的插值参数,对所述至少两组频谱谱线进行插值处理。
11.一种计算机可读存储介质,所述存储介质存储计算机指令,当计算机读取所述计算机指令时,所述计算机执行如权利要求1~8中任一项所述的多普勒超声信号处理方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310708389.4A CN116549011A (zh) | 2023-06-13 | 2023-06-13 | 一种多普勒超声信号处理方法、系统和存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310708389.4A CN116549011A (zh) | 2023-06-13 | 2023-06-13 | 一种多普勒超声信号处理方法、系统和存储介质 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116549011A true CN116549011A (zh) | 2023-08-08 |
Family
ID=87502013
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310708389.4A Pending CN116549011A (zh) | 2023-06-13 | 2023-06-13 | 一种多普勒超声信号处理方法、系统和存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116549011A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116989888A (zh) * | 2023-09-27 | 2023-11-03 | 之江实验室 | 声成像方法、装置、计算机设备和存储介质 |
-
2023
- 2023-06-13 CN CN202310708389.4A patent/CN116549011A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116989888A (zh) * | 2023-09-27 | 2023-11-03 | 之江实验室 | 声成像方法、装置、计算机设备和存储介质 |
CN116989888B (zh) * | 2023-09-27 | 2024-03-12 | 之江实验室 | 声成像方法、装置、计算机设备和存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11653900B2 (en) | Data augmentation for training deep learning models with ultrasound images | |
US8021301B2 (en) | Ultrasonic image processing apparatus, ultrasonic image processing method and ultrasonic image processing program | |
US20170109925A1 (en) | Utilizing depth from ultrasound volume rendering for 3d printing | |
CN102138808B (zh) | 超声波诊断装置及医用图像处理装置 | |
US20100310143A1 (en) | Filter Bank for Ultrasound Image Enhancement | |
JP2013505778A (ja) | 運動情報を用いた医用画像解析のためのコンピュータ可読媒体、システム、および方法 | |
CN116549011A (zh) | 一种多普勒超声信号处理方法、系统和存储介质 | |
JPH10155796A (ja) | 医療用超音波装置および該画像作成方法 | |
CN114376606B (zh) | 一种超声成像的滤波方法和系统 | |
US20200294230A1 (en) | Medical image processing apparatus, training data producing apparatus, and ultrasound diagnosis apparatus | |
KR20190044758A (ko) | 3차원 초음파 영상 복원 방법 및 그 초음파 영상 장치 | |
JP2011120881A (ja) | 副関心領域に基づいて3次元超音波映像を提供する超音波システムおよび方法 | |
US20110295116A1 (en) | Providing a motion image in an ultrasound system | |
CN115670520B (zh) | 血管内超声成像方法、装置、计算机设备、存储介质 | |
JPWO2014156269A1 (ja) | 超音波撮像装置及び超音波画像表示方法 | |
WO2020254159A1 (en) | Method and system for generating a synthetic elastrography image | |
CN111598989A (zh) | 一种图像渲染参数设置方法、装置、电子设备及存储介质 | |
JP7437192B2 (ja) | 医用画像処理装置 | |
Ketterling et al. | High-speed, high-frequency ultrasound, in utero vector-flow imaging of mouse embryos | |
JP2012050818A (ja) | カラードップラモード映像を提供する超音波システムおよび方法 | |
KR20160098010A (ko) | 초음파 진단 장치, 초음파 프로브 및 그 제어 방법 | |
JP6998477B2 (ja) | カラードップラー超音波イメージングを行うための方法及びシステム | |
JP7171625B2 (ja) | 超音波画像を処理するための方法及びシステム | |
US20220338839A1 (en) | Method, system, and storage medium for ultrasonic imaging | |
Song et al. | A novel design framework of synthetic radial aperture focusing for volumetric transrectal ultrasound imaging |
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 |