CN113768547A - 冠状动脉成像方法及装置、存储介质及电子设备 - Google Patents
冠状动脉成像方法及装置、存储介质及电子设备 Download PDFInfo
- Publication number
- CN113768547A CN113768547A CN202111074706.9A CN202111074706A CN113768547A CN 113768547 A CN113768547 A CN 113768547A CN 202111074706 A CN202111074706 A CN 202111074706A CN 113768547 A CN113768547 A CN 113768547A
- Authority
- CN
- China
- Prior art keywords
- data
- coronary artery
- radio frequency
- groups
- data corresponding
- 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.)
- Granted
Links
- 210000004351 coronary vessel Anatomy 0.000 title claims abstract description 106
- 238000003384 imaging method Methods 0.000 title claims abstract description 55
- 238000001914 filtration Methods 0.000 claims abstract description 51
- 230000017531 blood circulation Effects 0.000 claims abstract description 27
- 238000002604 ultrasonography Methods 0.000 claims description 41
- 238000000034 method Methods 0.000 claims description 16
- 238000005096 rolling process Methods 0.000 claims description 16
- 238000004590 computer program Methods 0.000 claims description 8
- 125000004122 cyclic group Chemical group 0.000 claims description 7
- 230000000875 corresponding effect Effects 0.000 description 153
- 210000001519 tissue Anatomy 0.000 description 21
- 210000005242 cardiac chamber Anatomy 0.000 description 10
- 230000008569 process Effects 0.000 description 9
- 230000000747 cardiac effect Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000003759 clinical diagnosis Methods 0.000 description 3
- 238000002586 coronary angiography Methods 0.000 description 3
- 230000034994 death Effects 0.000 description 3
- 231100000517 death Toxicity 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 210000000038 chest Anatomy 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 208000031225 myocardial ischemia Diseases 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000012285 ultrasound imaging Methods 0.000 description 2
- RUDATBOHQWOJDD-UHFFFAOYSA-N (3beta,5beta,7alpha)-3,7-Dihydroxycholan-24-oic acid Natural products OC1CC2CC(O)CCC2(C)C2C1C1CCC(C(CCC(O)=O)C)C1(C)CC2 RUDATBOHQWOJDD-UHFFFAOYSA-N 0.000 description 1
- HPTJABJPZMULFH-UHFFFAOYSA-N 12-[(Cyclohexylcarbamoyl)amino]dodecanoic acid Chemical compound OC(=O)CCCCCCCCCCCNC(=O)NC1CCCCC1 HPTJABJPZMULFH-UHFFFAOYSA-N 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000002583 angiography Methods 0.000 description 1
- 230000003143 atherosclerotic effect Effects 0.000 description 1
- 238000010009 beating Methods 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 239000002872 contrast media Substances 0.000 description 1
- 208000029078 coronary artery disease Diseases 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 208000019622 heart disease Diseases 0.000 description 1
- 230000005865 ionizing radiation Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 230000006798 recombination Effects 0.000 description 1
- 238000005215 recombination Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 210000001562 sternum Anatomy 0.000 description 1
- RUDATBOHQWOJDD-UZVSRGJWSA-N ursodeoxycholic acid Chemical compound C([C@H]1C[C@@H]2O)[C@H](O)CC[C@]1(C)[C@@H]1[C@@H]2[C@@H]2CC[C@H]([C@@H](CCC(O)=O)C)[C@@]2(C)CC1 RUDATBOHQWOJDD-UZVSRGJWSA-N 0.000 description 1
- 230000002861 ventricular Effects 0.000 description 1
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/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/48—Diagnostic techniques
-
- 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
-
- 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/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
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)
- Vascular Medicine (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本申请提供了一种冠状动脉成像方法及装置、存储介质及电子设备，涉及超声成像技术领域。本申请实施例提供的冠状动脉成像方法，通过确定冠状动脉对应的第一IQ数据，提高了射频数据对应的帧率，其中，第一IQ数据对应M帧图像；通过基于预设组合信息，确定第一IQ数据对应的W组第二IQ数据的方式，提供了冠状动脉超声造影信号与心腔内血流超声造影信号、组织超声造影信号的分辨信息；通过分别对W组第二IQ数据进行滤波操作，得到W组第二IQ数据各自对应的第三IQ数据的方式，滤除了心腔内血流超声造影信号和组织超声造影信号；基于W组第二IQ数据各自对应的第三IQ数据，重建得到了冠状动脉对应的冠状动脉超声造影图像。
Description
技术领域
本申请涉及超声成像技术领域，具体涉及一种冠状动脉成像方法及装置、存储介质及电子设备。
背景技术
冠状动脉造影图像是诊断冠状动脉粥样硬化性心脏病(即冠心病)的一种常用而且有效的方式，现已广泛应用于临床诊断。
然而，在对冠状动脉超声造影成像的同时，由于心腔内血流超声造影信号和组织超声造影信号的干扰，冠状动脉超声造影信号会被遮盖，因此无法成像。并且，现有技术中的超快速超声冠状动脉成像方法均为开胸成像，会受到成像角度的限制，同时无法观察冠状动脉超声造影信号的细节。
发明内容
为了解决上述技术问题，提出了本申请。本申请实施例提供了一种冠状动脉成像方法及装置、存储介质及电子设备。
第一方面，本申请一实施例提供了一种冠状动脉成像方法，该冠状动脉成像方法包括：确定冠状动脉对应的第一IQ数据，其中，第一IQ数据对应M帧图像，M为大于1的正整数；基于预设组合信息，确定第一IQ数据对应的W组第二IQ数据，其中，W为小于M的正整数；分别对W组第二IQ数据进行滤波操作，得到W组第二IQ数据各自对应的第三IQ数据，其中，滤波操作用于滤除心腔内血流超声造影信号和组织超声造影信号；基于W组第二IQ数据各自对应的第三IQ数据，重建得到冠状动脉对应的冠状动脉超声造影图像。
结合第一方面，在第一方面的某些实现方式中，确定冠状动脉对应的第一IQ数据，包括：获取冠状动脉对应的射频数据，其中，射频数据基于N个不同角度的平面波信号得到，N为大于1的正整数；对循环复用N个不同角度的射频数据进行滚动重建操作，得到第一IQ数据，其中，滚动重建操作用于提高射频数据对应的帧率。
结合第一方面，在第一方面的某些实现方式中，对循环复用N个不同角度的射频数据进行滚动重建操作，得到第一IQ数据，包括：针对N个不同角度的射频数据中的第P个角度的射频数据，基于第P个角度的射频数据、第P个角度的射频数据之后的循环复用的N-1个不同角度各自对应的射频数据，确定第P个角度的射频数据对应的IQ数据，其中，P为小于或等于N的正整数；基于N个不同角度的射频数据各自对应的IQ数据，确定第一IQ数据。
结合第一方面，在第一方面的某些实现方式中，预设组合信息包括预设时间间隔信息和预设时间段信息，基于预设组合信息，确定第一IQ数据对应的W组第二IQ数据，包括：基于M帧图像对应的第一IQ数据中的第U帧图像对应的第一IQ数据、预设时间段信息，确定第1组第二IQ数据，其中，U为小于M的正整数；基于第1组第二IQ数据、预设时间间隔信息依次向后滑动，确定第2组至第W组第二IQ数据。
结合第一方面，在第一方面的某些实现方式中，分别对W组第二IQ数据进行滤波操作，得到W组第二IQ数据各自对应的第三IQ数据，包括：确定W组第二IQ数据各自对应的像素坐标集；基于W组第二IQ数据各自对应的像素坐标集，分别对W组第二IQ数据进行滤波操作，得到W组第二IQ数据各自对应的第三IQ数据。
结合第一方面，在第一方面的某些实现方式中，基于W组第二IQ数据各自对应的像素坐标集，分别对W组第二IQ数据进行滤波操作，得到W组第二IQ数据各自对应的第三IQ数据，包括：针对W组第二IQ数据中的每组第二IQ数据，基于第二IQ数据对应的像素坐标集包括的S个像素坐标，确定S个像素坐标各自对应的IQ数据集，其中，IQ数据集包括IQ数据集对应的像素坐标，S为正整数；基于预设频率阈值对S个像素坐标各自对应的IQ数据集进行滤波操作，得到每组第二IQ数据对应的第三IQ数据。
结合第一方面，在第一方面的某些实现方式中，冠状动脉超声造影图像的帧率大于或等于1千赫兹。
第二方面，本申请一实施例提供了一种冠状动脉成像装置，该冠状动脉成像装置包括：第一确定模块，用于确定冠状动脉对应的第一IQ数据，其中，第一IQ数据对应M帧图像，M为大于1的正整数；第二确定模块，用于基于预设组合信息，确定第一IQ数据对应的W组第二IQ数据，其中，W为小于M的正整数；滤波操作模块，用于分别对W组第二IQ数据进行滤波操作，得到W组第二IQ数据各自对应的第三IQ数据，其中，滤波操作用于滤除心腔内血流超声造影信号和组织超声造影信号；重建模块，用于基于W组第二IQ数据各自对应的第三IQ数据，重建得到冠状动脉对应的冠状动脉超声造影图像。
第三方面，本申请一实施例提供了一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序用于执行上述任一实施例提及的冠状动脉成像方法。
第四方面，本申请一实施例提供了一种电子设备，该电子设备包括：处理器；用于存储处理器可执行指令的存储器；该处理器用于执行上述任一实施例提及的冠状动脉成像方法。
本申请实施例提供的冠状动脉成像方法，通过确定冠状动脉对应的第一IQ数据的方式，提高了射频数据对应的帧率，其中，第一IQ数据对应M帧图像，M为大于1的正整数；通过基于预设组合信息，确定第一IQ数据对应的W组第二IQ数据的方式，提供了冠状动脉超声造影信号与心腔内血流超声造影信号、组织超声造影信号的分辨信息，其中，W为小于M的正整数；通过分别对W组第二IQ数据进行滤波操作，得到W组第二IQ数据各自对应的第三IQ数据的方式，滤除了心腔内血流超声造影信号和组织超声造影信号；基于W组第二IQ数据各自对应的第三IQ数据，重建得到了冠状动脉对应的冠状动脉超声造影图像。
附图说明
通过结合附图对本申请实施例进行更详细的描述，本申请的上述以及其他目的、特征和优势将变得更加明显。附图用来提供对本申请实施例的进一步理解，并且构成说明书的一部分，与本申请实施例一起用于解释本申请，并不构成对本申请的限制。在附图中，相同的参考标号通常代表相同部件或步骤。
图1所示为本申请一示例性实施例提供的冠状动脉成像方法的流程示意图。
图2所示为本申请一示例性实施例提供的确定冠状动脉对应的第一IQ数据的流程示意图。
图3所示为本申请一示例性实施例提供的对循环复用N个不同角度的射频数据进行滚动重建操作，得到第一IQ数据的流程示意图。
图4所示为本申请一示例性实施例提供的基于预设组合信息，确定第一IQ数据对应的W组第二IQ数据的流程示意图。
图5所示为本申请一示例性实施例提供的分别对W组第二IQ数据进行滤波操作，得到W组第二IQ数据各自对应的第三IQ数据的流程示意图。
图6所示为本申请一示例性实施例提供的基于W组第二IQ数据各自对应的像素坐标集，分别对W组第二IQ数据进行滤波操作，得到W组第二IQ数据各自对应的第三IQ数据的流程示意图。
图7所示为本申请一示例性实施例提供的冠状动脉成像装置的结构示意图。
图8所示为本申请一示例性实施例提供的电子设备的结构示意图。
具体实施方式
下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
缺血性心脏病是世界十大死因之首，因缺血性心脏病死亡的人数约占全球总死亡人数的16％。目前临床中，如DSA、CTA、MRA等冠状动脉造影方式，既不方便，还受到如分辨率、侵入性和电离辐射等限制。在科研领域中，基于CUDA或3D UDCA对冠状动脉进行超声多普勒成像，但均为开胸成像，会受到成像角度的限制，同时无法观察冠状动脉超声造影信号的细节。
另外，心脏的跳动导致造影剂持续时间短，超声帧率不足以捕捉心动周期的每一个阶段。心脏泵血过程中，心腔内高浓度的超声造影信号会遮盖冠状动脉超声造影信号，并且由于胸骨的存在，阻碍了传向心脏的超声声束。具体地，冠状动脉位于心室壁，在对冠状动脉成像的同时，也对心腔内的血流成像，但相对于冠状动脉内的血流而言，心腔内的血流流量大、流速高，因此心腔内血流超声造影信号会远远强于冠状动脉超声造影信号，冠状动脉超声造影信号因此被遮盖，无法成像。同时，冠状动脉成像需要滤除组织超声造影信号，保证冠状动脉超声造影信号不被组织超声造影信号遮盖。
图1所示为本申请一示例性实施例提供的冠状动脉成像方法的流程示意图。如图1所示，本申请实施例提供的冠状动脉成像方法包括如下步骤。
步骤10，确定冠状动脉对应的第一IQ数据。
具体地，第一IQ数据提高了射频数据对应的帧率，第一IQ数据对应M帧图像，M为大于1的正整数。
在本申请一实施例中，第一IQ数据可以是基于平面波相干角度复合技术(即滚动重建操作)，对射频数据进行预处理后得到。
步骤20，基于预设组合信息，确定第一IQ数据对应的W组第二IQ数据。
具体地，步骤20中提及的W组第二IQ数据中的每组第二IQ数据，提供了冠状动脉超声造影信号与心腔内血流超声造影信号、组织超声造影信号的分辨信息，基于分辨信息可以滤除心腔内血流超声造影信号和组织超声造影信号。另外，W为小于M的正整数。
步骤30，分别对W组第二IQ数据进行滤波操作，得到W组第二IQ数据各自对应的第三IQ数据。
具体地，步骤30中提及的滤波操作，滤除了心腔内血流超声造影信号和组织超声造影信号，保留了冠状动脉超声造影信号，进而得到冠状动脉超声造影图像。
步骤40，基于W组第二IQ数据各自对应的第三IQ数据，重建得到冠状动脉对应的冠状动脉超声造影图像。
在本申请一实施例中，冠状动脉超声造影图像的帧率大于或等于1千赫兹。本申请实施例提供的冠状动脉成像方法，帧率满足捕捉心动周期每一个阶段的条件，提高了临床诊断的准确性。
在实际应用过程中，首先通过确定冠状动脉对应的第一IQ数据的方式，提高了射频数据对应的帧率；继而通过基于预设组合信息，确定第一IQ数据对应的W组第二IQ数据的方式，提供了冠状动脉超声造影信号与心腔内血流超声造影信号、组织超声造影信号的分辨信息；然后通过分别对W组第二IQ数据进行滤波操作，得到W组第二IQ数据各自对应的第三IQ数据的方式，滤除了心腔内血流超声造影信号和组织超声造影信号；最后基于W组第二IQ数据各自对应的第三IQ数据，重建得到了冠状动脉对应的冠状动脉超声造影图像。本申请实施例得到的冠状动脉超声造影图像显示了心动周期的每一个阶段，提高了临床诊断的准确性。
图2所示为本申请一示例性实施例提供的确定冠状动脉对应的第一IQ数据的流程示意图。在本申请图1所示实施例的基础上延伸出本申请图2所示实施例，下面着重叙述图2所示实施例与图1所示实施例的不同之处，相同之处不再赘述。
如图2所示，在本申请实施例提供的冠状动脉成像方法中，确定冠状动脉对应的第一IQ数据步骤，包括如下步骤。
步骤11，获取冠状动脉对应的射频数据。
在本申请一实施例中，基于超快速超声成像设备获取冠状动脉对应的射频数据，具体地，射频数据基于N个不同角度的平面波信号得到，N为大于1的正整数。
步骤12，对循环复用N个不同角度的射频数据进行滚动重建操作，得到第一IQ数据。
在本申请一实施例中，滚动重建操作提高了射频数据对应的帧率。另外，由于不同角度下的平面波信号相关而噪声独立，使得滚动重建操作得到的第一IQ数据提高了冠状动脉超声造影图像的穿透力及横向分辨率。
本申请实施例提供的冠状动脉成像方法，通过首先获取冠状动脉对应的射频数据，继而对射频数据进行滚动重建操作，得到第一IQ数据的方式，提高了射频数据对应的帧率，并且第一IQ数据提高了冠状动脉超声造影图像的穿透力及横向分辨率。
图3所示为本申请一示例性实施例提供的对循环复用N个不同角度的射频数据进行滚动重建操作，得到第一IQ数据的流程示意图。在本申请图2所示实施例的基础上延伸出本申请图3所示实施例，下面着重叙述图3所示实施例与图2所示实施例的不同之处，相同之处不再赘述。
如图3所示，在本申请实施例提供的冠状动脉成像方法中，对循环复用N个不同角度的射频数据进行滚动重建操作，得到第一IQ数据步骤，包括如下步骤。
步骤121，针对N个不同角度的射频数据中的第P个角度的射频数据，基于第P个角度的射频数据、第P个角度的射频数据之后的循环复用的N-1个不同角度各自对应的射频数据，确定第P个角度的射频数据对应的IQ数据。
示例性地，N等于10，P等于3，针对10个不同角度的射频数据中的第3个角度的射频数据，基于第3个角度的射频数据至第3个角度的射频数据的下一个循环复用的共计10个不同角度各自对应的射频数据，确定第3个角度的射频数据对应的IQ数据。具体地，基于第3个角度的射频数据至第3个角度的射频数据的下一个循环复用的第2个角度的射频数据，确定第3个角度的射频数据对应的IQ数据。其中，P为小于或等于N的正整数。
步骤122，基于N个不同角度的射频数据各自对应的IQ数据，确定第一IQ数据。
示例性地，利用步骤121中提及的N等于10，P等于3，基于10个不同角度的射频数据各自对应的IQ数据，确定第一IQ数据。
需要说明的是，本申请实施例提供的确定第一IQ数据的方法，将N个不同角度的射频数据中的每个角度的射频数据复用N次，使得到的冠状动脉超声造影图像的帧率提升至原帧率的N倍，为之后找出冠状动脉超声造影信号与心腔内血流超声造影信号、组织超声造影信号的分辨信息提供了前提条件。
本申请实施例提供的冠状动脉成像方法，针对N个不同角度的射频数据中的第P个角度的射频数据，通过基于第P个角度的射频数据、第P个角度的射频数据之后的循环复用的N-1个不同角度各自对应的射频数据的方式，得到了第P个角度的射频数据对应的IQ数据，继而基于N个不同角度的射频数据各自对应的IQ数据，得到了第一IQ数据。另外，本申请实施例提供的确定第一IQ数据的方法，将N个不同角度的射频数据中的每个角度的射频数据复用N次，使得到的冠状动脉超声造影图像的帧率提升至原帧率的N倍，为之后找出冠状动脉超声造影信号与心腔内血流超声造影信号、组织超声造影信号的分辨信息提供了前提条件。
图4所示为本申请一示例性实施例提供的基于预设组合信息，确定第一IQ数据对应的W组第二IQ数据的流程示意图。在本申请图1所示实施例的基础上延伸出本申请图4所示实施例，下面着重叙述图4所示实施例与图1所示实施例的不同之处，相同之处不再赘述。
如图4所示，在本申请实施例提供的冠状动脉成像方法中，基于预设组合信息，确定第一IQ数据对应的W组第二IQ数据步骤，包括如下步骤。具体地，预设组合信息包括预设时间间隔信息和预设时间段信息。
步骤21，基于M帧图像对应的第一IQ数据中的第U帧图像对应的第一IQ数据、预设时间段信息，确定第1组第二IQ数据。
在本申请一实施例中，预设时间段信息可以是M帧图像对应的第一IQ数据中的连续Q帧图像对应的第一IQ数据，基于第U帧图像对应的第一IQ数据、第U帧图像对应的第一IQ数据之后的连续Q-1帧图像对应的第一IQ数据，即基于第U帧图像对应的第一IQ数据至第Q+U-1帧图像对应的第一IQ数据，确定第1组第二IQ数据。其中，Q为小于M的正整数。
示例性地，M等于100，U等于3，基于100帧图像对应的第一IQ数据中的第3帧图像对应的第一IQ数据、预设时间段信息，确定第1组第二IQ数据。具体地，预设时间段信息可以是100帧图像对应的第一IQ数据中的连续30帧图像对应的第一IQ数据，基于第3帧图像对应的第一IQ数据、第3帧图像对应的第一IQ数据之后的连续29帧图像对应的第一IQ数据，即基于第3帧图像对应的第一IQ数据至第32帧图像对应的第一IQ数据，确定第1组第二IQ数据。其中，U为小于M的正整数。
需要说明的是，预设时间段信息可以根据能够提供冠状动脉超声造影信号与心腔内血流超声造影信号、组织超声造影信号的分辨信息这一条件而设定。示例性地，分辨信息可以是，在一定时间内，冠状动脉超声造影信号与心腔内血流超声造影信号、组织超声造影信号对应的频率不同。
步骤22，基于第1组第二IQ数据、预设时间间隔信息依次向后滑动，确定第2组至第W组第二IQ数据。
在本申请一实施例中，预设时间间隔信息可以是M帧图像对应的第一IQ数据中的连续R帧图像对应的第一IQ数据，其中，R为小于M的正整数。示例性地，利用步骤21中提及的M等于100，U等于3，预设时间间隔信息可以是1帧图像对应的第一IQ数据，基于第1组第二IQ数据(即第3帧图像对应的第一IQ数据至第32帧图像对应的第一IQ数据)，向后滑动1帧图像对应的第一IQ数据(即预设时间间隔信息)，得到第2组第二IQ数据(即第4帧图像对应的第一IQ数据至第33帧图像对应的第一IQ数据)。基于相同的依次向后滑动原理，依次得到第3组至第W组第二IQ数据。
优选地，预设时间间隔信息小于或等于预设时间段信息，以将M帧图像对应的第一IQ数据经过滑动操作后，M帧图像对应的第一IQ数据中的任意一帧图像对应的第一IQ数据都被包含在与其对应的W组第二IQ数据中的任意一组第二IQ数据中。
本申请实施例提供的冠状动脉成像方法，基于M帧图像对应的第一IQ数据中的第U帧图像对应的第一IQ数据、预设时间段信息，得到了第1组第二IQ数据；基于第1组第二IQ数据、预设时间间隔信息依次向后滑动，得到了第2组至第W组第二IQ数据。其中，W组第二IQ数据中的每组第二IQ数据提供了冠状动脉超声造影信号与心腔内血流超声造影信号、组织超声造影信号的分辨信息，以便于之后滤除心腔内血流超声造影信号、组织超声造影信号。
图5所示为本申请一示例性实施例提供的分别对W组第二IQ数据进行滤波操作，得到W组第二IQ数据各自对应的第三IQ数据的流程示意图。在本申请图1所示实施例的基础上延伸出本申请图5所示实施例，下面着重叙述图5所示实施例与图1所示实施例的不同之处，相同之处不再赘述。
如图5所示，在本申请实施例提供的冠状动脉成像方法中，分别对W组第二IQ数据进行滤波操作，得到W组第二IQ数据各自对应的第三IQ数据步骤，包括如下步骤。
步骤31，确定W组第二IQ数据各自对应的像素坐标集。
示例性地，步骤31中提及的像素坐标集可以是W组第二IQ数据各自对应的像素坐标的集合。
步骤32，基于W组第二IQ数据各自对应的像素坐标集，分别对W组第二IQ数据进行滤波操作，得到W组第二IQ数据各自对应的第三IQ数据。
示例性地，分别对W组第二IQ数据各自对应的像素坐标集进行滤波操作，得到W组第二IQ数据各自对应的第三IQ数据。
本申请实施例提供的冠状动脉成像方法，首先确定W组第二IQ数据各自对应的像素坐标集，继而基于W组第二IQ数据各自对应的像素坐标集，分别对W组第二IQ数据进行滤波操作，得到了W组第二IQ数据各自对应的第三IQ数据，从而滤除了心腔内血流超声造影信号和组织超声造影信号，为得到冠状动脉超声造影图像提供了前提条件。
图6所示为本申请一示例性实施例提供的基于W组第二IQ数据各自对应的像素坐标集，分别对W组第二IQ数据进行滤波操作，得到W组第二IQ数据各自对应的第三IQ数据的流程示意图。在本申请图5所示实施例的基础上延伸出本申请图6所示实施例，下面着重叙述图6所示实施例与图5所示实施例的不同之处，相同之处不再赘述。
如图6所示，在本申请实施例提供的冠状动脉成像方法中，基于W组第二IQ数据各自对应的像素坐标集，分别对W组第二IQ数据进行滤波操作，得到W组第二IQ数据各自对应的第三IQ数据步骤，包括如下步骤。
步骤321，针对W组第二IQ数据中的每组第二IQ数据，基于第二IQ数据对应的像素坐标集包括的S个像素坐标，确定S个像素坐标各自对应的IQ数据集。
具体地，IQ数据集包括IQ数据集对应的像素坐标，S为正整数。
示例性地，针对S个像素坐标中的第1个像素坐标，确定第1个像素坐标对应的IQ数据集，进而针对S个像素坐标中的第2个像素坐标至第S个像素坐标，确定第2个像素坐标至第S个像素坐标各自对应的IQ数据集。
步骤322，基于预设频率阈值对S个像素坐标各自对应的IQ数据集进行滤波操作，得到每组第二IQ数据对应的第三IQ数据。
示例性地，预设频率阈值可以根据频率分辨信息而设定，该频率分辨信息可以是能够将冠状动脉超声造影信号与心腔内血流超声造影信号、组织超声造影信号对应的频率进行分辨的频率信息。例如，将IQ数据集中不满足预设频率阈值的元素滤除，以得到每组第二IQ数据对应的第三IQ数据。
示例性地，利用步骤321中提及的第1个像素坐标对应的IQ数据集，以及第2个像素坐标至第S个像素坐标各自对应的IQ数据集，基于预设频率阈值对第1个像素坐标对应的IQ数据集进行滤波操作，并对第2个像素坐标至第S个像素坐标各自对应的IQ数据集进行滤波操作，得到每组第二IQ数据对应的第三IQ数据。
本申请实施例提供的冠状动脉成像方法，针对W组第二IQ数据中的每组第二IQ数据，通过基于第二IQ数据对应的像素坐标集包括的S个像素坐标，得到了S个像素坐标各自对应的IQ数据集，继而基于预设频率阈值对S个像素坐标各自对应的IQ数据集进行滤波操作，得到了每组第二IQ数据对应的第三IQ数据，第三IQ数据只对应有冠状动脉超声造影信号，从而为得到冠状动脉超声造影图像提供了前提条件。
图7所示为本申请一示例性实施例提供的冠状动脉成像装置的结构示意图。如图7所示，本申请实施例提供的冠状动脉成像装置包括：
第一确定模块100，用于确定冠状动脉对应的第一IQ数据，其中，第一IQ数据对应M帧图像，M为大于1的正整数；
第二确定模块200，用于基于预设组合信息，确定第一IQ数据对应的W组第二IQ数据，其中，W为小于M的正整数；
滤波操作模块300，用于分别对W组第二IQ数据进行滤波操作，得到W组第二IQ数据各自对应的第三IQ数据，其中，滤波操作用于滤除心腔内血流超声造影信号和组织超声造影信号；
重建模块400，用于基于W组第二IQ数据各自对应的第三IQ数据，重建得到冠状动脉对应的冠状动脉超声造影图像。
在本申请一实施例中，第一确定模块100，还用于获取冠状动脉对应的射频数据，其中，射频数据基于N个不同角度的平面波信号得到，N为大于1的正整数；对循环复用N个不同角度的射频数据进行滚动重建操作，得到第一IQ数据，其中，滚动重建操作用于提高射频数据对应的帧率。
在本申请一实施例中，第一确定模块100，还用于针对N个不同角度的射频数据中的第P个角度的射频数据，基于第P个角度的射频数据、第P个角度的射频数据之后的循环复用的N-1个不同角度各自对应的射频数据，确定第P个角度的射频数据对应的IQ数据，其中，P为小于或等于N的正整数；基于N个不同角度的射频数据各自对应的IQ数据，确定第一IQ数据。
在本申请一实施例中，第二确定模块200，还用于基于M帧图像对应的第一IQ数据中的第U帧图像对应的第一IQ数据、预设时间段信息，确定第1组第二IQ数据，其中，U为小于M的正整数；基于第1组第二IQ数据、预设时间间隔信息依次向后滑动，确定第2组至第W组第二IQ数据。
在本申请一实施例中，滤波操作模块300，还用于确定W组第二IQ数据各自对应的像素坐标集；基于W组第二IQ数据各自对应的像素坐标集，分别对W组第二IQ数据进行滤波操作，得到W组第二IQ数据各自对应的第三IQ数据。
在本申请一实施例中，滤波操作模块300，还用于针对W组第二IQ数据中的每组第二IQ数据，基于第二IQ数据对应的像素坐标集包括的S个像素坐标，确定S个像素坐标各自对应的IQ数据集，其中，IQ数据集包括IQ数据集对应的像素坐标，S为正整数；基于预设频率阈值对S个像素坐标各自对应的IQ数据集进行滤波操作，得到每组第二IQ数据对应的第三IQ数据。
应当理解，图7提供的冠状动脉成像装置中的第一确定模块100、第二确定模块200、滤波操作模块300以及重建模块400的操作和功能可以参考上述图1至图6提供的冠状动脉成像方法，为了避免重复，在此不再赘述。
下面，参考图8来描述根据本申请实施例的电子设备。图8所示为本申请一示例性实施例提供的电子设备的结构示意图。
如图8所示，电子设备50包括一个或多个处理器501和存储器502。
处理器501可以是中央处理单元(CPU)或者具有数据处理能力和/或指令执行能力的其他形式的处理单元，并且可以控制电子设备50中的其他组件以执行期望的功能。
存储器502可以包括一个或多个计算机程序产品，所述计算机程序产品可以包括各种形式的计算机可读存储介质，例如易失性存储器和/或非易失性存储器。所述易失性存储器例如可以包括随机存取存储器(RAM)和/或高速缓冲存储器(cache)等。所述非易失性存储器例如可以包括只读存储器(ROM)、硬盘、闪存等。在所述计算机可读存储介质上可以存储一个或多个计算机程序指令，处理器501可以运行所述程序指令，以实现上文所述的本申请的各个实施例的冠状动脉成像方法以及/或者其他期望的功能。在所述计算机可读存储介质中还可以存储诸如包括射频数据、第一IQ数据、第二IQ数据、第三IQ数据、冠状动脉超声造影图像等各种内容。
在一个示例中，电子设备50还可以包括：输入装置503和输出装置504，这些组件通过总线系统和/或其他形式的连接机构(未示出)互连。
该输入装置503可以包括例如键盘、鼠标等等。
该输出装置504可以向外部输出各种信息，包括射频数据、第一IQ数据、第二IQ数据、第三IQ数据、冠状动脉超声造影图像等。该输出装置504可以包括例如显示器、扬声器、打印机、以及通信网络及其所连接的远程输出设备等等。
当然，为了简化，图8中仅示出了该电子设备50中与本申请有关的组件中的一些，省略了诸如总线、输入/输出接口等等的组件。除此之外，根据具体应用情况，电子设备50还可以包括任何其他适当的组件。
除了上述方法和设备以外，本申请的实施例还可以是计算机程序产品，其包括计算机程序指令，所述计算机程序指令在被处理器运行时使得所述处理器执行本说明书上述描述的根据本申请各种实施例的冠状动脉成像方法中的步骤。
所述计算机程序产品可以以一种或多种程序设计语言的任意组合来编写用于执行本申请实施例操作的程序代码，所述程序设计语言包括面向对象的程序设计语言，诸如Java、C++等，还包括常规的过程式程序设计语言，诸如“C”语言或类似的程序设计语言。程序代码可以完全地在用户计算设备上执行、部分地在用户设备上执行、作为一个独立的软件包执行、部分在用户计算设备上部分在远程计算设备上执行、或者完全在远程计算设备或服务器上执行。
此外，本申请的实施例还可以是计算机可读存储介质，其上存储有计算机程序指令，所述计算机程序指令在被处理器运行时使得所述处理器执行本说明书上述描述的根据本申请各种实施例的冠状动脉成像方法中的步骤。
所述计算机可读存储介质可以采用一个或多个可读介质的任意组合。可读介质可以是可读信号介质或者可读存储介质。可读存储介质例如可以包括但不限于电、磁、光、电磁、红外线、或半导体的系统、装置或器件，或者任意以上的组合。可读存储介质的更具体的例子(非穷举的列表)包括：具有一个或多个导线的电连接、便携式盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。
以上结合具体实施例描述了本申请的基本原理，但是，需要指出的是，在本申请中提及的优点、优势、效果等仅是示例而非限制，不能认为这些优点、优势、效果等是本申请的各个实施例必须具备的。另外，上述公开的具体细节仅是为了示例的作用和便于理解的作用，而非限制，上述细节并不限制本申请为必须采用上述具体的细节来实现。
本申请中涉及的器件、装置、设备、系统的方框图仅作为例示性的例子并且不意图要求或暗示必须按照方框图示出的方式进行连接、布置、配置。如本领域技术人员将认识到的，可以按任意方式连接、布置、配置这些器件、装置、设备、系统。诸如“包括”、“包含”、“具有”等等的词语是开放性词汇，指“包括但不限于”，且可与其互换使用。这里所使用的词汇“或”和“和”指词汇“和/或”，且可与其互换使用，除非上下文明确指示不是如此。这里所使用的词汇“诸如”指词组“诸如但不限于”，且可与其互换使用。
还需要指出的是，在本申请的装置、设备和方法中，各部件或各步骤是可以分解和/或重新组合的。这些分解和/或重新组合应视为本申请的等效方案。
提供所公开的方面的以上描述以使本领域的任何技术人员能够做出或者使用本申请。对这些方面的各种修改对于本领域技术人员而言是非常显而易见的，并且在此定义的一般原理可以应用于其他方面而不脱离本申请的范围。因此，本申请不意图被限制到在此示出的方面，而是按照与在此公开的原理和新颖的特征一致的最宽范围。
为了例示和描述的目的已经给出了以上描述。此外，此描述不意图将本申请的实施例限制到在此公开的形式。尽管以上已经讨论了多个示例方面和实施例，但是本领域技术人员将认识到其某些变型、修改、改变、添加和子组合。
Claims (10)
1.一种冠状动脉成像方法，其特征在于，包括：
确定冠状动脉对应的第一IQ数据，其中，所述第一IQ数据对应M帧图像，M为大于1的正整数；
基于预设组合信息，确定所述第一IQ数据对应的W组第二IQ数据，其中，W为小于M的正整数；
分别对所述W组第二IQ数据进行滤波操作，得到所述W组第二IQ数据各自对应的第三IQ数据，其中，所述滤波操作用于滤除心腔内血流超声造影信号和组织超声造影信号；
基于所述W组第二IQ数据各自对应的第三IQ数据，重建得到所述冠状动脉对应的冠状动脉超声造影图像。
2.根据权利要求1所述的冠状动脉成像方法，其特征在于，所述确定冠状动脉对应的第一IQ数据，包括：
获取所述冠状动脉对应的射频数据，其中，所述射频数据基于N个不同角度的平面波信号得到，N为大于1的正整数；
对循环复用N个不同角度的所述射频数据进行滚动重建操作，得到所述第一IQ数据，其中，所述滚动重建操作用于提高所述射频数据对应的帧率。
3.根据权利要求2所述的冠状动脉成像方法，其特征在于，所述对循环复用N个不同角度的所述射频数据进行滚动重建操作，得到所述第一IQ数据，包括：
针对所述N个不同角度的射频数据中的第P个角度的射频数据，基于所述第P个角度的射频数据、所述第P个角度的射频数据之后的循环复用的N-1个不同角度各自对应的射频数据，确定所述第P个角度的射频数据对应的IQ数据，其中，P为小于或等于N的正整数；
基于所述N个不同角度的射频数据各自对应的IQ数据，确定所述第一IQ数据。
4.根据权利要求1至3中的任一项所述的冠状动脉成像方法，其特征在于，所述预设组合信息包括预设时间间隔信息和预设时间段信息，所述基于预设组合信息，确定所述第一IQ数据对应的W组第二IQ数据，包括：
基于所述M帧图像对应的所述第一IQ数据中的第U帧图像对应的所述第一IQ数据、所述预设时间段信息，确定第1组第二IQ数据，其中，U为小于M的正整数；
基于所述第1组第二IQ数据、所述预设时间间隔信息依次向后滑动，确定第2组至第W组第二IQ数据。
5.根据权利要求1至3中的任一项所述的冠状动脉成像方法，其特征在于，所述分别对所述W组第二IQ数据进行滤波操作，得到所述W组第二IQ数据各自对应的第三IQ数据，包括：
确定所述W组第二IQ数据各自对应的像素坐标集；
基于所述W组第二IQ数据各自对应的像素坐标集，分别对所述W组第二IQ数据进行滤波操作，得到所述W组第二IQ数据各自对应的第三IQ数据。
6.根据权利要求5所述的冠状动脉成像方法，其特征在于，所述基于所述W组第二IQ数据各自对应的像素坐标集，分别对所述W组第二IQ数据进行滤波操作，得到所述W组第二IQ数据各自对应的第三IQ数据，包括：
针对所述W组第二IQ数据中的每组第二IQ数据，
基于所述第二IQ数据对应的像素坐标集包括的S个像素坐标，确定所述S个像素坐标各自对应的IQ数据集，其中，所述IQ数据集包括所述IQ数据集对应的像素坐标，S为正整数；
基于预设频率阈值对所述S个像素坐标各自对应的IQ数据集进行滤波操作，得到所述每组第二IQ数据对应的第三IQ数据。
7.根据权利要求1至3中的任一项所述的冠状动脉成像方法，其特征在于，所述冠状动脉超声造影图像的帧率大于或等于1千赫兹。
8.一种冠状动脉成像装置，其特征在于，包括：
第一确定模块，用于确定冠状动脉对应的第一IQ数据，其中，所述第一IQ数据对应M帧图像，M为大于1的正整数；
第二确定模块，用于基于预设组合信息，确定所述第一IQ数据对应的W组第二IQ数据，其中，W为小于M的正整数；
滤波操作模块，用于分别对所述W组第二IQ数据进行滤波操作，得到所述W组第二IQ数据各自对应的第三IQ数据，其中，所述滤波操作用于滤除心腔内血流超声造影信号和组织超声造影信号；
重建模块，用于基于所述W组第二IQ数据各自对应的第三IQ数据，重建得到所述冠状动脉对应的冠状动脉超声造影图像。
9.一种计算机可读存储介质，其特征在于，所述存储介质存储有计算机程序，所述计算机程序用于执行上述权利要求1至7中的任一项所述的冠状动脉成像方法。
10.一种电子设备，其特征在于，所述电子设备包括：
处理器；
用于存储所述处理器可执行指令的存储器；
所述处理器，用于执行上述权利要求1至7中的任一项所述的冠状动脉成像方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111074706.9A CN113768547B (zh) | 2021-09-14 | 2021-09-14 | 冠状动脉成像方法及装置、存储介质及电子设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111074706.9A CN113768547B (zh) | 2021-09-14 | 2021-09-14 | 冠状动脉成像方法及装置、存储介质及电子设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113768547A true CN113768547A (zh) | 2021-12-10 |
| CN113768547B CN113768547B (zh) | 2024-03-22 |
Family
ID=78843495
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111074706.9A Active CN113768547B (zh) | 2021-09-14 | 2021-09-14 | 冠状动脉成像方法及装置、存储介质及电子设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113768547B (zh) |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101625757A (zh) * | 2009-08-13 | 2010-01-13 | 华北电力大学（保定） | 一种抑制冠状动脉内超声图像序列中刚性运动伪影的方法 |
| CN107392994A (zh) * | 2017-06-30 | 2017-11-24 | 深圳大学 | 冠脉血管的三维重建方法、装置、设备及存储介质 |
| CN108135566A (zh) * | 2016-06-30 | 2018-06-08 | 北京深迈瑞医疗电子技术研究院有限公司 | 一种增强超声造影图像的方法、系统及超声造影成像设备 |
| WO2020077381A1 (en) * | 2018-10-17 | 2020-04-23 | The University Of Melbourne | Coronary evaluation system and process |
| CN111110277A (zh) * | 2019-12-27 | 2020-05-08 | 深圳开立生物医疗科技股份有限公司 | 超声成像方法、超声设备及存储介质 |
| CN111588410A (zh) * | 2020-05-18 | 2020-08-28 | 南京超维景生物科技有限公司 | 基于超声造影图像的图像处理方法及装置 |
| CN111728643A (zh) * | 2020-06-24 | 2020-10-02 | 南京超维景生物科技有限公司 | 单血管的超声血流灌注成像方法及其装置 |
| CN112184690A (zh) * | 2020-10-12 | 2021-01-05 | 推想医疗科技股份有限公司 | 冠脉血管走向的预测方法、预测模型的训练方法及装置 |
| CN112652052A (zh) * | 2020-12-15 | 2021-04-13 | 山东大学 | 一种基于血管分支配准的冠状动脉三维重建方法及系统 |
-
2021
- 2021-09-14 CN CN202111074706.9A patent/CN113768547B/zh active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101625757A (zh) * | 2009-08-13 | 2010-01-13 | 华北电力大学（保定） | 一种抑制冠状动脉内超声图像序列中刚性运动伪影的方法 |
| CN108135566A (zh) * | 2016-06-30 | 2018-06-08 | 北京深迈瑞医疗电子技术研究院有限公司 | 一种增强超声造影图像的方法、系统及超声造影成像设备 |
| CN107392994A (zh) * | 2017-06-30 | 2017-11-24 | 深圳大学 | 冠脉血管的三维重建方法、装置、设备及存储介质 |
| US20190117198A1 (en) * | 2017-06-30 | 2019-04-25 | Shenzhen University | Three-dimensional reconstruction method, apparatus and device and storage medium for coronary vessels |
| WO2020077381A1 (en) * | 2018-10-17 | 2020-04-23 | The University Of Melbourne | Coronary evaluation system and process |
| CN111110277A (zh) * | 2019-12-27 | 2020-05-08 | 深圳开立生物医疗科技股份有限公司 | 超声成像方法、超声设备及存储介质 |
| CN111588410A (zh) * | 2020-05-18 | 2020-08-28 | 南京超维景生物科技有限公司 | 基于超声造影图像的图像处理方法及装置 |
| CN111728643A (zh) * | 2020-06-24 | 2020-10-02 | 南京超维景生物科技有限公司 | 单血管的超声血流灌注成像方法及其装置 |
| CN112184690A (zh) * | 2020-10-12 | 2021-01-05 | 推想医疗科技股份有限公司 | 冠脉血管走向的预测方法、预测模型的训练方法及装置 |
| CN112652052A (zh) * | 2020-12-15 | 2021-04-13 | 山东大学 | 一种基于血管分支配准的冠状动脉三维重建方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113768547B (zh) | 2024-03-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Song et al. | Accelerated singular value-based ultrasound blood flow clutter filtering with randomized singular value decomposition and randomized spatial downsampling | |
| US20220142583A1 (en) | Methods and systems to quantify and remove asynchronous noise in biophysical signals | |
| US20180220997A1 (en) | System and method for accelerated clutter filtering in ultrasound blood flow imaging using randomized ultrasound data | |
| CN101743565B (zh) | 图像序列的滤波 | |
| WO2018134729A1 (en) | Sparsity-based ultrasound super-resolution imaging | |
| US10123761B2 (en) | Device and method for spatiotemporal reconstruction of a moving vascular pulse wave in the brain and other organs | |
| US20220296216A1 (en) | Image processing method and apparatus based on contrast-enhanced ultrasound images | |
| CN117858671A (zh) | 无电路心动周期确定 | |
| CN112581554B (zh) | 一种ct成像方法、装置、存储设备及医学成像系统 | |
| JP2022516163A (ja) | 脳病変情報提供装置及び方法 | |
| CN113768547B (zh) | 冠状动脉成像方法及装置、存储介质及电子设备 | |
| CN114027804A (zh) | 一种脉象诊断的方法、装置及可读存储介质 | |
| Niederhauser et al. | Graphics-processor-unit-based parallelization of optimized baseline wander filtering algorithms for long-term electrocardiography | |
| CN113616217B (zh) | 基线漂移曲线的生成方法和装置 | |
| Agrawal et al. | FPGA-based peak detection of ECG signal using histogram approach | |
| JP2015213714A (ja) | 心電波形検出装置、心電波形検出方法、心電波形検出プログラム、及び撮像装置 | |
| EP4059442B1 (en) | Image processing method and apparatus based on ultrasound contrast image | |
| Islam et al. | Denoising ECG signal using adaptive filter algorithms and cubic spline interpolation for regaining missing data points of ECG in telecardiology system | |
| CN111248894B (zh) | 使用ecg恢复脉搏波形及其他血液动力学的设备及方法 | |
| CN113842128B (zh) | 一种基于多重滤波与混合放大的非接触心率检测装置 | |
| JP2022528223A (ja) | 血液量パルス信号検出装置、血液量パルス信号検出方法、及びプログラム | |
| Farid et al. | Removal of power line interference from ECG signal using different adaptive filter algorithms and cubic spline interpolation for missing data points of EGC in telecardiology system | |
| EP4111953A1 (en) | Method and device for translating between stochastic signals | |
| CN117349600B (zh) | 一种基于双模态双输入的心音心电联合诊断方法及系统 | |
| CN112842264B (zh) | 多模态成像中数字滤波方法、装置和多模态成像技术系统 |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |