CN104287777B - 超音波扫描方法及超音波扫描*** - Google Patents
超音波扫描方法及超音波扫描*** Download PDFInfo
- Publication number
- CN104287777B CN104287777B CN201410553268.8A CN201410553268A CN104287777B CN 104287777 B CN104287777 B CN 104287777B CN 201410553268 A CN201410553268 A CN 201410553268A CN 104287777 B CN104287777 B CN 104287777B
- Authority
- CN
- China
- Prior art keywords
- image
- binaryzation
- ultrasonic
- images
- noise
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8995—Combining images from different aspect angles, e.g. spatial compounding
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—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/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
- A61B8/5238—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52046—Techniques for image enhancement involving transmitter or receiver
- G01S7/52047—Techniques for image enhancement involving transmitter or receiver for elimination of side lobes or of grating lobes; for increasing resolving power
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Animal Behavior & Ethology (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biophysics (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Acoustics & Sound (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
本发明提供一种超音波扫描方法,包含下列步骤:对标的物朝N个不同角度发出N组超音波讯号且接收自该标的物反射及/或散射的该N组超音波讯号,其中每一组该超音波讯号分别对应该N个不同角度的其中之一,且N为大于1的正整数;将该N组超音波讯号转换为N张超音波输入影像;对该N张超音波输入影像执行二值化演算法,以得到N张二值化影像;对该N张二值化影像进行比较程序,以判断该N张二值化影像的至少其中之一是否存在杂讯;当该N张二值化影像中的第i张二值化影像存在该杂讯时,将该杂讯自对应该第i张二值化影像的第i张超音波输入影像中滤除,其中i为小于或等于N的正整数;以及将该N张超音波输入影像合成为超音波输出影像。
Description
技术领域
本发明关于一种超音波扫描方法及超音波扫描***,尤指一种可有效滤除杂讯的超音波扫描方法及超音波扫描***。
背景技术
由于超音波扫描具有不破坏材料结构以及人体细胞的特性,因而普遍地被应用于材料领域以及临床医学检测。由于复合影像可提供较佳的解析度,复合影像已广泛地运用于超音波扫描***上。复合影像的技术原理必须先自多个不同角度发射与接收超音波讯号,再将对应此多个不同角度的超音波影像合成为一张超音波复合影像。然而,从斜角发射与接收的超音波讯号所产生的超音波影像常会产生旁瓣(side lobe)、假影(artifact)或其它杂讯,进而影响超音波复合影像的准确性。因此,如何有效滤除杂讯便成为超音波扫描技术中的一个重要的研究课题。
发明内容
本发明的目的之一在于提供一种可有效滤除杂讯的超音波扫描方法及超音波扫描***,以解决上述问题。
一方面,本发明提供一种超音波扫描方法,包含下列步骤:
对标的物朝N个不同角度发出N组超音波讯号且接收自该标的物反射及/或散射的该N组超音波讯号,其中每一组该超音波讯号分别对应该N个不同角度的其中之一,且N为大于1的正整数;
将该N组超音波讯号转换为N张超音波输入影像;
对该N张超音波输入影像执行二值化演算法,以得到N张二值化影像;
对该N张二值化影像进行比较程序,以判断该N张二值化影像的至少其中之一是否存在杂讯;
当该N张二值化影像中的第i张二值化影像存在该杂讯时,将该杂讯自对应该第i张二值化影像的第i张超音波输入影像中滤除,其中i为小于或等于N的正整数;以及
将该N张超音波输入影像合成为超音波输出影像。
较佳的,对该N张超音波输入影像执行二值化演算法,以得到N张二值化影像的步骤包含:
选取该N张超音波输入影像的其中之一作为参考影像;
将每一张该超音波输入影像分割成对应的M个区块,其中M为正整数;
将该参考影像中各该区块的所有像素值平均,以得到M个二值化阀值;以及
以对应各该区块的该二值化阀值对该N张超音波输入影像中对应的各该区块执行该二值化演算法,以得到该N张二值化影像;
其中,该比较程序包含下列步骤:
将除了该参考影像外的每一张该超音波输入影像的该二值化影像中的M个区块分别与该参考影像的该二值化影像中对应的M个区块进行比较;以及
当两对应区块的像素差异量大于杂讯滤除阀值时,判断非该参考影像的该超音波输入影像中的该区块存在该杂讯。
较佳的,对应该参考影像的该超音波讯号是朝该N个不同角度其中的最小角度发出。
较佳的,该最小角度为0度。
较佳的,对该N张超音波输入影像执行二值化演算法,以得到N张二值化影像的步骤包含:
将每一张该超音波输入影像分割成对应的M个区块,其中M为正整数;
将该N张超音波输入影像中对应的各该区块的所有像素值平均,以得到M个二值化阀值;以及
以对应各该区块的该二值化阀值对该N张超音波输入影像中对应的各该区块执行该二值化演算法,以得到该N张二值化影像;
其中,该比较程序包含下列步骤:
将该N张二值化影像中对应的该M个区块相互进行比较;以及
当该M个区块中的特定区块与其它对应的区块的像素差异量皆大于杂讯滤除阀值时,判断该特定区块存在该杂讯。
另一方面,本发明提供一种超音波扫描***,包含超音波探头以及超音波扫描仪。超音波探头用以对标的物朝N个不同角度发出N组超音波讯号且接收自该标的物反射及/或散射的该N组超音波讯号,其中每一组该超音波讯号分别对应该N个不同角度的其中之一,且N为大于1的正整数。超音波扫描仪用以与该超音波探头形成通讯。该超音波扫描仪包含:转换单元、二值化单元、比较单元、杂讯滤除单元以及影像合成单元。该转换单元用以将该N组超音波讯号转换为N张超音波输入影像;该二值化单元用以对该N张超音波输入影像执行二值化演算法,以得到N张二值化影像;该比较单元用以对该N张二值化影像进行比较程序,以判断该N张二值化影像的至少其中之一是否存在杂讯;该杂讯滤除单元,用以当该N张二值化影像中的第i张二值化影像存在该杂讯时,将该杂讯自对应该第i张二值化影像的第i张超音波输入影像中滤除,其中i为小于或等于N的正整数;该影像合成单元用以将该N张超音波输入影像合成为超音波输出影像。
较佳的,该二值化单元选取该N张超音波输入影像的其中之一作为参考影像,将每一张该超音波输入影像分割成对应的M个区块,将该参考影像中各该区块的所有像素值平均,以得到M个二值化阀值,且以对应各该区块的该二值化阀值对该N张超音波输入影像中对应的各该区块执行该二值化演算法,以得到该N张二值化影像,其中M为正整数;
该比较单元将除了该参考影像外的每一张该超音波输入影像的该二值化影像中的M个区块分别与该参考影像的该二值化影像中对应的M个区块进行比较;以及
当两对应区块的像素差异量大于杂讯滤除阀值时,该比较单元判断非该参考影像的该超音波输入影像中的该区块存在该杂讯,以判断该N张二值化影像的至少其中之一是否存在该杂讯。
较佳的,对应该参考影像的该超音波讯号是朝该N个不同角度其中的最小角度发出。
较佳的,该最小角度为0度。
较佳的,该二值化单元将每一张该超音波输入影像分割成对应的M个区块,将该N张超音波输入影像中对应的各该区块的所有像素值平均,以得到M个二值化阀值,且以对应各该区块的该二值化阀值对该N张超音波输入影像中对应的各该区块执行该二值化演算法,以得到该N张二值化影像,其中M为正整数;
该比较单元将该N张二值化影像中对应的该M个区块相互进行比较;以及
当该M个区块中的特定区块与其它对应的区块的像素差异量皆大于杂讯滤除阀值时,该比较单元判断该特定区块存在该杂讯,以判断该N张二值化影像的至少其中之一是否存在该杂讯。
与现有技术相比,本发明提供的超音波扫描方法及超音波扫描***,在得到对应多个不同角度的多张超音波输入影像后,先将此多张超音波输入影像转换为多张二值化影像,再对此多张二值化影像进行比较程序,以判断二值化影像是否存在杂讯。当二值化影像存在杂讯时,本发明即将杂讯自对应此二值化影像的超音波输入影像中滤除。在将杂讯自超音波输入影像中滤除后,再将多张超音波输入影像合成为单一超音波输出影像。由于本发明在合成超音波输出影像前先将杂讯自超音波输入影像中滤除,因此,合成后的超音波输出影像便不会受到杂讯影响,可有效提高超音波输出影像的准确性。
附图说明
图1为根据本发明一实施例的超音波扫描***的功能方块图;
图2为以图1中的超音波探头对标的物进行超音波扫描的示意图;
图3为图1中的超音波扫描仪根据图2中的扫描结果所产生的超音波输入影像的示意图;
图4为图1中的超音波扫描仪将图3中的超音波输入影像转换为二值化影像的示意图;
图5为图3中的超音波输入影像滤除杂讯后合成为超音波输出影像的示意图;
图6为根据本发明一实施例的超音波扫描方法的流程图;
图7为图3中的每一个超音波输入影像被分割成对应的九个区块的示意图;
图8为图6中的步骤S14、S16的一实施例的详细流程图;
图9为图6中的步骤S14、S16的另一实施例的详细流程图。
具体实施方式
为使对本发明的目的、构造、特征、及其功能有进一步的了解,兹配合实施例详细说明如下。
请参阅图1至图6,图1为根据本发明一实施例的超音波扫描***1的功能方块图,图2为以图1中的超音波探头10对标的物3进行超音波扫描的示意图,图3为图1中的超音波扫描仪12根据图2中的扫描结果所产生的超音波输入影像UI1、UI2、UI3的示意图,图4为图1中的超音波扫描仪12将图3中的超音波输入影像UI1、UI2、UI3转换为二值化影像BI1、BI2、BI3的示意图,图5为图3中的超音波输入影像UI1、UI2、UI3滤除杂讯后合成为超音波输出影像UI4的示意图,图6为根据本发明一实施例的超音波扫描方法的流程图。图6中的超音波扫描方法适用于图1中的超音波扫描***1。
如图1所示,超音波扫描***1包含超音波探头10以及超音波扫描仪12,其中超音波扫描仪12与超音波探头10形成通讯。超音波扫描仪12可为电脑或其它具有资料运算、处理与显示功能的电子装置。超音波扫描仪12与超音波探头10可藉由缆线(cable)形成通讯,使得超音波探头10在对标的物进行超音波扫描时,可于超音波扫描仪12的显示器上显示超音波扫描影像。
于此实施例中,超音波扫描仪12包含转换单元120、二值化单元122、比较单元124、杂讯滤除单元126以及影像合成单元128,其中二值化单元122电性连接于转换单元120,比较单元124电性连接于二值化单元122,杂讯滤除单元126电性连接于比较单元124,且影像合成单元128电性连接于转换单元120与杂讯滤除单元126。
如图2所示,在以超音波探头10对标的物3进行超音波扫描时,操作人员可操作超音波探头10对标的物3朝N个不同角度发出N组超音波讯号且接收自标的物3反射及/或散射的N组超音波讯号(图6中的步骤S10),其中每一组超音波讯号分别对应N个不同角度的其中之一,且N为大于1的正整数。每一组超音波讯号包含一束至多束的超音波讯号,而每一组超音波讯号所包含的超音波讯号个数视超音波探头10的规格而定。以图2所绘示的实施例为例,超音波探头10对标的物3朝三个不同角度θ1、θ2、θ3发出三组超音波讯号S1、S2、S3且接收自标的物3反射及/或散射的三组超音波讯号S1、S2、S3(亦即,N=3),其中角度θ1可为正30度,角度θ2可为0度,且角度θ3可为负30度,但不以此为限。需说明的是,超音波讯号的发射角度与数量可根据实际应用来决定,不以上述实施例为限。
接着,超音波探头10将三组超音波讯号S1、S2、S3传送至超音波扫描仪12。超音波扫描仪12在接收三组超音波讯号S1、S2、S3后,超音波扫描仪12的转换单元120会将三组超音波讯号S1、S2、S3转换为三张超音波输入影像UI1、UI2、UI3(图6中的步骤S12),如图3所示。于此实施例中,超音波输入影像UI1是由对应正30度的角度的超音波讯号S1转换得到,超音波输入影像UI2是由对应0度的角度的超音波讯号S2转换得到,且超音波输入影像UI3是由对应负30度的角度的超音波讯号S3转换得到。
接着,超音波扫描仪12的二值化单元122对三张超音波输入影像UI1、UI2、UI3执行二值化演算法,以得到N张二值化影像BI1、BI2、BI3(图6中的步骤S14),如图4所示。
接着,超音波扫描仪12的比较单元124对三张二值化影像BI1、BI2、BI3进行比较程序,以判断三张二值化影像BI1、BI2、BI3的至少其中之一是否存在杂讯(图6中的步骤S16)。当三张二值化影像BI1、BI2、BI3中的第i张二值化影像存在杂讯时,超音波扫描仪12的杂讯滤除单元126将杂讯自对应第i张二值化影像的第i张超音波输入影像中滤除(图6中的步骤S18),其中i为小于或等于N的正整数。于此实施例中,超音波扫描仪12的比较单元124判断三张二值化影像BI1、BI2、BI3中的第1张二值化影像BI1存在杂讯N1,且第3张二值化影像BI3存在杂讯N3,如图4所示。因此,超音波扫描仪12的杂讯滤除单元126会将杂讯N1自对应第1张二值化影像的第1张超音波输入影像UI1中滤除,且将杂讯N3自对应第3张二值化影像的第3张超音波输入影像UI3中滤除。
在将杂讯N1、N3分别自超音波输入影像UI1、UI3中滤除后,超音波扫描仪12的影像合成单元128即会将三张超音波输入影像UI1、UI2、UI3合成为一超音波输出影像UI4(图6中的步骤S20),如图5所示。最后,超音波扫描仪12的显示器上即会显示超音波输出影像UI4。由于本发明在合成超音波输出影像UI4前即先将杂讯N1、N3分别自超音波输入影像UI1、UI3中滤除,因此,合成后的超音波输出影像UI4便不会受到杂讯N1、N3影响,可有效提高超音波输出影像UI4的准确性。需说明的是,上述的杂讯N1、N3可为超音波影像中常见的旁瓣(sidelobe)、假影(artifact)、白色杂斑或其它杂讯。
请参阅图7以及图8,图7为图3中的每一个超音波输入影像UI1、UI2、UI3被分割成对应的九个区块的示意图,图8为图6中的步骤S14、S16的一实施例的详细流程图。于此实施例中,图6中的步骤S14可包含图8中的步骤S30至S36,且图6中的步骤S16中的比较程序可包含图8中的步骤S38至S40。
在超音波扫描仪12的转换单元120将三组超音波讯号S1、S2、S3转换为三张超音波输入影像UI1、UI2、UI3后(图6中的步骤S12),超音波扫描仪12的二值化单元122可先选取三张超音波输入影像UI1、UI2、UI3的其中之一作为参考影像(图8中的步骤S30)。以此实施例中,本发明可以朝三个不同角度θ1、θ2、θ3其中的最小角度发出的超音波讯号转换得到的超音波输入影像作为上述的参考影像。换言之,对应参考影像的超音波讯号朝三个不同角度θ1、θ2、θ3其中的最小角度发出。以图2所绘示的实施例为例,三个不同角度θ1、θ2、θ3其中的最小角度为等于0度的θ2,因此,超音波扫描仪12的二值化单元122会选取超音波输入影像UI2作为参考影像。
接着,超音波扫描仪12的二值化单元122可将每一张超音波输入影像UI1、UI2、UI3分割成对应的M个区块(图8中的步骤S32),其中M为正整数。如图7所示,超音波扫描仪12的二值化单元122可将每一张超音波输入影像UI1、UI2、UI3分割成对应的九个区块(亦即,M=9),但不以此为限。超音波扫描仪12的二值化单元122可根据实际应用的需求将每一张超音波输入影像UI1、UI2、UI3分割成对应的一个或多个区块。如图7所示,超音波输入影像UI1、UI2、UI3的区块C11-C19、C21-C29、C31-C39分别相互对应。
接着,超音波扫描仪12的二值化单元122可将参考影像UI2中各区块C21-C29的所有像素值平均,以得到九个二值化阀值(图8中的步骤S34)。进一步来说,二值化单元122可将区块C21的所有像素值平均而得到一个二值化阀值,二值化单元122可将区块C22的所有像素值平均而得到另一个二值化阀值,以此类推,进而得到九个二值化阀值。
接着,超音波扫描仪12的二值化单元122可以对应各区块C21-C29的二值化阀值对三张超音波输入影像UI1、UI2、UI3中对应的各区块C11-C19、C21-C29、C31-C39执行二值化演算法,以得到三张二值化影像BI1、BI2、BI3(图8中的步骤S36)。进一步来说,二值化单元122可以对应区块C21的二值化阀值对三张超音波输入影像UI1、UI2、UI3中对应的各区块C11、C21、C31执行二值化演算法,以对应区块C22的二值化阀值对三张超音波输入影像UI1、UI2、UI3中对应的各区块C12、C22、C32执行二值化演算法,以此类推,进而得到三张二值化影像BI1、BI2、BI3。
接着,超音波扫描仪12的比较单元124可将除了参考影像UI2外的每一张超音波输入影像UI1、UI3的二值化影像BI1、BI3中的区块C11-C19、C31-C39分别与参考影像UI2的二值化影像BI2中对应的区块C21-C29进行比较(图8中的步骤S38)。当两对应区块的像素差异量大于杂讯滤除阀值时,超音波扫描仪12的比较单元124即会判断非参考影像UI2的超音波输入影像中的区块存在杂讯(图8中的步骤S40)。
以图7所绘示的实施例为例,假设杂讯滤除阀值设定为10,且二值化影像BI1与二值化影像BI2的两对应区块C11、C21的像素差异量为30,超音波扫描仪12的比较单元124即会判断超音波输入影像UI1中的区块C11存在杂讯N1;同理,假设杂讯滤除阀值设定为10,且二值化影像BI3与二值化影像BI2的两对应区块C33、C23的像素差异量为20,超音波扫描仪12的比较单元124即会判断超音波输入影像UI3中的区块C33存在杂讯N3。此外,由于二值化影像BI1与二值化影像BI2的其它对应区块C12-C19、C22-C29的像素差异量皆为0或小于杂讯滤除阀值,则超音波扫描仪12的比较单元124即会判断超音波输入影像UI1中的其它区块C12-C19皆不存在杂讯;同理,由于二值化影像BI3与二值化影像BI2的其它对应区块C31-C32、C34-C39、C21-C22、C24-C29的像素差异量皆为0或小于杂讯滤除阀值,则超音波扫描仪12的比较单元124即会判断超音波输入影像UI3中的其它区块C31-C32、C34-C39皆不存在杂讯。
接着,再执行图6中的步骤S18、S20,超音波扫描仪12的影像合成单元128即会将滤除杂讯N1、N3后的三张超音波输入影像UI1、UI2、UI3合成为超音波输出影像UI4。
图9,图9为图6中的步骤S14、S16的另一实施例的详细流程图。于此实施例中,图6中的步骤S14可包含图9中的步骤S50至S54,且图6中的步骤S16中的比较程序可包含图9中的步骤S56至S58。以下同样利用图7搭配图9来进行说明。
在超音波扫描仪12的转换单元120将三组超音波讯号S1、S2、S3转换为三张超音波输入影像UI1、UI2、UI3后(图6中的步骤S12),超音波扫描仪12的二值化单元可将每一张超音波输入影像UI1、UI2、UI3分割成对应的M个区块(图9中的步骤S50),其中M为正整数。如图7所示,超音波扫描仪12的二值化单元122可将每一张超音波输入影像UI1、UI2、UI3分割成对应的九个区块(亦即,M=9),但不以此为限。超音波扫描仪12的二值化单元122可根据实际应用的需求将每一张超音波输入影像UI1、UI2、UI3分割成对应的一个或多个区块。如图7所示,超音波输入影像UI1、UI2、UI3的区块C11-C19、C21-C29、C31-C39分别相互对应。
接着,超音波扫描仪12的二值化单元122可将三张超音波输入影像UI1、UI2、UI3中对应的各区块C11-C19、C21-C29、C31-C39的所有像素值平均,以得到九个二值化阀值(图8中的步骤S52)。进一步来说,二值化单元122可将三张超音波输入影像UI1、UI2、UI3中对应的区块C11、C21、C31的所有像素值平均而得到一个二值化阀值,二值化单元122可将三张超音波输入影像UI1、UI2、UI3中对应的区块C12、C22、C32的所有像素值平均而得到另一个二值化阀值,以此类推,进而得到九个二值化阀值。
接着,超音波扫描仪12的二值化单元122可以对应各区块C11-C19、C21-C29、C31-C39的二值化阀值对三张超音波输入影像UI1、UI2、UI3中对应的各区块C11-C19、C21-C29、C31-C39执行二值化演算法,以得到三张二值化影像BI1、BI2、BI3(图9中的步骤S54)。进一步来说,二值化单元122可以对应区块C11、C21、C31的二值化阀值对三张超音波输入影像UI1、UI2、UI3中对应的各区块C11、C21、C31执行二值化演算法,以对应区块C12、C22、C32的二值化阀值对三张超音波输入影像UI1、UI2、UI3中对应的各区块C12、C22、C32执行二值化演算法,以此类推,进而得到三张二值化影像BI1、BI2、BI3。
接着,超音波扫描仪12的比较单元124可将二值化影像BI1、BI2、BI3中对应的区块C11-C19、C21-C29、C31-C39相互进行比较(图9中的步骤S56)。当区块C11-C19、C21-C29、C31-C39中的特定区块与其它对应的区块的像素差异量皆大于杂讯滤除阀值时,超音波扫描仪12的比较单元124即会判断此特定区块存在杂讯(图9中的步骤S58)。
以图7所绘示的实施例为例,假设杂讯滤除阀值设定为10,二值化影像BI1与二值化影像BI2的两对应区块C11、C21的像素差异量为30,二值化影像BI1与二值化影像BI3的两对应区块C11、C31的像素差异量为30,且二值化影像BI2与二值化影像BI3的两对应区块C21、C31的像素差异量为0或小于杂讯滤除阀值,由于区块C11与其它对应的区块C21、C31的像素差异量皆大于杂讯滤除阀值,超音波扫描仪12的比较单元124即会判断超音波输入影像UI1中的特定区块C11存在杂讯N1;同理,假设杂讯滤除阀值设定为10,二值化影像BI1与二值化影像BI2的两对应区块C13、C23的像素差异量为0或小于杂讯滤除阀值,二值化影像BI1与二值化影像BI3的两对应区块C13、C33的像素差异量为20,且二值化影像BI2与二值化影像BI3的两对应区块C23、C33的像素差异量为20,由于区块C33与其它对应的区块C13、C23的像素差异量皆大于杂讯滤除阀值,超音波扫描仪12的比较单元124即会判断超音波输入影像UI3中的特定区块C33存在杂讯N3。此外,由于二值化影像BI1、BI2、BI3的其它对应区块C12、C14-C19、C22、C24-C29、C32、C34-C39的像素差异量皆为0或小于杂讯滤除阀值,则超音波扫描仪12的比较单元124即会判断超音波输入影像UI1、UI2、UI3中的其它区块C12、C14-C19、C22、C24-C29、C32、C34-C39皆不存在杂讯。
接着,再执行图6中的步骤S18、S20,超音波扫描仪12的影像合成单元128即会将滤除杂讯N1、N3后的三张超音波输入影像UI1、UI2、UI3合成为超音波输出影像UI4。
需说明的是,图6所示的步骤S10至S20、图8所示的步骤S30至S40以及图9所示的步骤S50至S58的控制逻辑可以软件设计来实现。当然,控制逻辑中的各个部分或功能皆可通过软件、硬件或软硬件的组合来实现。此外,图1所示的转换单元120、二值化单元122、比较单元124、杂讯滤除单元126与影像合成单元128可以电路设计来实现。
综上所述,本发明在得到对应多个不同角度的多张超音波输入影像后,先将此多张超音波输入影像转换为多张二值化影像,再对此多张二值化影像进行比较程序,以判断二值化影像是否存在杂讯。当二值化影像存在杂讯时,本发明即将杂讯自对应此二值化影像的超音波输入影像中滤除。在将杂讯自超音波输入影像中滤除后,再将多张超音波输入影像合成为单一超音波输出影像。由于本发明在合成超音波输出影像前即先将杂讯自超音波输入影像中滤除,因此,合成后的超音波输出影像便不会受到杂讯影响,可有效提高超音波输出影像的准确性。于上述的比较程序中,本发明可根据实际应用的需求,自多张超音波输入影像中选取一张超音波输入影像作为比较基准的参考影像,或是以多张超音波输入影像两两相互比较,以判断二值化影像是否存在杂讯。
本发明已由上述相关实施例加以描述,然而上述实施例仅为实施本发明的范例。必需指出的是,已揭露的实施例并未限制本发明的范围。相反地,在不脱离本发明的精神和范围内所作的更动与润饰,均属本发明的专利保护范围。
Claims (6)
1.一种超音波扫描方法,其特征在于,包含下列步骤:
对标的物朝N个不同角度发出N组超音波讯号且接收自该标的物反射及/或散射的该N组超音波讯号,其中每一组该超音波讯号分别对应该N个不同角度的其中之一,且N为大于1的正整数;
将该N组超音波讯号转换为N张超音波输入影像;
对该N张超音波输入影像执行二值化演算法,以得到N张二值化影像;
对该N张二值化影像进行比较程序,以判断该N张二值化影像的至少其中之一是否存在杂讯;
当该N张二值化影像中的第i张二值化影像存在该杂讯时,将该杂讯自对应该第i张二值化影像的第i张超音波输入影像中滤除,其中i为小于或等于N的正整数;以及
将该N张超音波输入影像合成为超音波输出影像;
对该N张超音波输入影像执行二值化演算法,以得到N张二值化影像的步骤包含:
选取该N张超音波输入影像的其中之一作为参考影像;
将每一张该超音波输入影像分割成对应的M个区块,其中M为正整数;
将该参考影像中各该区块的所有像素值平均,以得到M个二值化阀值;以及
以对应各该区块的该二值化阀值对该N张超音波输入影像中对应的各该区块执行该二值化演算法,以得到该N张二值化影像;
其中,该比较程序包含下列步骤:
将除了该参考影像外的每一张该超音波输入影像的该二值化影像中的M个区块分别与该参考影像的该二值化影像中对应的M个区块进行比较;以及
当两对应区块的像素差异量大于杂讯滤除阀值时,判断非该参考影像的该超音波输入影像中的该区块存在该杂讯。
2.如权利要求1所述的超音波扫描方法,其特征在于,对应该参考影像的该超音波讯号是朝该N个不同角度其中的最小角度发出。
3.如权利要求2所述的超音波扫描方法,其特征在于,该最小角度为0度。
4.一种超音波扫描***,其特征在于,包含:
超音波探头,用以对标的物朝N个不同角度发出N组超音波讯号且接收自该标的物反射及/或散射的该N组超音波讯号,其中每一组该超音波讯号分别对应该N个不同角度的其中之一,且N为大于1的正整数;以及
超音波扫描仪,用以与该超音波探头形成通讯,该超音波扫描仪包含:
转换单元,用以将该N组超音波讯号转换为N张超音波输入影像;
二值化单元,用以对该N张超音波输入影像执行二值化演算法,以得到N张二值化影像;
比较单元,用以对该N张二值化影像进行比较程序,以判断该N张二值化影像的至少其中之一是否存在杂讯;
杂讯滤除单元,用以当该N张二值化影像中的第i张二值化影像存在该杂讯时,将该杂讯自对应该第i张二值化影像的第i张超音波输入影像中滤除,其中i为小于或等于N的正整数;以及
影像合成单元,用以将该N张超音波输入影像合成为超音波输出影像;
该二值化单元选取该N张超音波输入影像的其中之一作为参考影像,将每一张该超音波输入影像分割成对应的M个区块,将该参考影像中各该区块的所有像素值平均,以得到M个二值化阀值,且以对应各该区块的该二值化阀值对该N张超音波输入影像中对应的各该区块执行该二值化演算法,以得到该N张二值化影像,其中M为正整数;
该比较单元将除了该参考影像外的每一张该超音波输入影像的该二值化影像中的M个区块分别与该参考影像的该二值化影像中对应的M个区块进行比较;以及
当两对应区块的像素差异量大于杂讯滤除阀值时,该比较单元判断非该参考影像的该超音波输入影像中的该区块存在该杂讯,以判断该N张二值化影像的至少其中之一是否存在该杂讯。
5.如权利要求4所述的超音波扫描***,其特征在于,对应该参考影像的该超音波讯号是朝该N个不同角度其中的最小角度发出。
6.如权利要求5所述的超音波扫描***,其特征在于,该最小角度为0度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410553268.8A CN104287777B (zh) | 2014-10-17 | 2014-10-17 | 超音波扫描方法及超音波扫描*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410553268.8A CN104287777B (zh) | 2014-10-17 | 2014-10-17 | 超音波扫描方法及超音波扫描*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104287777A CN104287777A (zh) | 2015-01-21 |
CN104287777B true CN104287777B (zh) | 2017-09-26 |
Family
ID=52307918
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410553268.8A Expired - Fee Related CN104287777B (zh) | 2014-10-17 | 2014-10-17 | 超音波扫描方法及超音波扫描*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104287777B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105806473B (zh) * | 2016-03-14 | 2018-12-11 | 苏州佳世达电通有限公司 | 超音波探头的参数设定方法 |
TWI644526B (zh) * | 2017-02-21 | 2018-12-11 | 佳世達科技股份有限公司 | 超音波系統及雜訊消除方法 |
CN113311440B (zh) * | 2020-02-26 | 2024-06-14 | 苏州佳世达电通有限公司 | 用于改善超声波影像质量的方法及其探测*** |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH05192336A (ja) * | 1992-01-23 | 1993-08-03 | Toshiba Corp | 超音波診断装置 |
JPH09220228A (ja) * | 1996-02-19 | 1997-08-26 | Ge Yokogawa Medical Syst Ltd | 血流情報画像表示方法および超音波診断装置 |
JP4473981B2 (ja) * | 1999-07-21 | 2010-06-02 | 株式会社日立メディコ | 超音波診断装置 |
CN101854537A (zh) * | 2009-03-31 | 2010-10-06 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声图像数据优化及其造影定量分析方法和*** |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SG120301A1 (en) * | 2004-08-27 | 2006-03-28 | Univ Nanyang | Home ultrasound system |
US8617075B2 (en) * | 2005-02-09 | 2013-12-31 | Hitachi Medical Corporation | Ultrasonic diagnostic apparatus and ultrasonic imaging method |
US20110196236A1 (en) * | 2010-02-08 | 2011-08-11 | Gokul Swamy | System and method of automated gestational age assessment of fetus |
US8712505B2 (en) * | 2010-11-11 | 2014-04-29 | University Of Pittsburgh-Of The Commonwealth System Of Higher Education | Automated macular pathology diagnosis in three-dimensional (3D) spectral domain optical coherence tomography (SD-OCT) images |
EP2797514A1 (en) * | 2011-12-30 | 2014-11-05 | Koninklijke Philips N.V. | Method and apparatus for needle visualization enhancement in ultrasound imaging |
JP5272084B2 (ja) * | 2012-01-11 | 2013-08-28 | 日立アロカメディカル株式会社 | 超音波診断装置 |
US9642592B2 (en) * | 2013-01-03 | 2017-05-09 | Siemens Medical Solutions Usa, Inc. | Needle enhancement in diagnostic ultrasound imaging |
-
2014
- 2014-10-17 CN CN201410553268.8A patent/CN104287777B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH05192336A (ja) * | 1992-01-23 | 1993-08-03 | Toshiba Corp | 超音波診断装置 |
JPH09220228A (ja) * | 1996-02-19 | 1997-08-26 | Ge Yokogawa Medical Syst Ltd | 血流情報画像表示方法および超音波診断装置 |
JP4473981B2 (ja) * | 1999-07-21 | 2010-06-02 | 株式会社日立メディコ | 超音波診断装置 |
CN101854537A (zh) * | 2009-03-31 | 2010-10-06 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声图像数据优化及其造影定量分析方法和*** |
Also Published As
Publication number | Publication date |
---|---|
CN104287777A (zh) | 2015-01-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6263094B1 (en) | Acoustic data acquisition/playback system and method | |
CN104287777B (zh) | 超音波扫描方法及超音波扫描*** | |
Alzahrani et al. | A real-time edge detector: algorithm and VLSI architecture | |
CN106680825B (zh) | 一种声学阵列成像***与方法 | |
JP5121389B2 (ja) | 対象体の大きさを測定するための超音波診断装置及び方法 | |
KR101911734B1 (ko) | 소프트웨어 기반의 초음파 이미징 시스템 | |
TWI653033B (zh) | 超音波成像系統及超音波成像方法 | |
CN102727255A (zh) | 一种超声图像空间复合成像的方法及装置 | |
CN108523922A (zh) | 多频成像方法、装置及设备 | |
US20070167795A1 (en) | Image processing system and method for controlling gains for color flow images | |
CN108156452A (zh) | 一种检测传感器的方法、装置、设备及存储介质 | |
JP2006508729A (ja) | 高フレームレート3次元超音波画像化器 | |
Petrucci et al. | A development of a new image analysis technique for detecting the flame front evolution in spark ignition engine under lean condition | |
TWI521225B (zh) | 超音波掃描方法及超音波掃描系統 | |
CN108389174B (zh) | 超音波成像***及超音波成像方法 | |
TW200937342A (en) | Line drawing processor, program and line drawing processing method | |
KR101627319B1 (ko) | 의료진단을 위한 의료영상 처리장치 및 그 방법 | |
CN115439509A (zh) | 一种多目标跟踪方法、装置、计算机设备及存储介质 | |
US6370264B1 (en) | Method and apparatus for ultrasonic color flow imaging | |
JP7305509B2 (ja) | 情報処理装置、その制御方法、プログラム、及び記憶媒体 | |
CN109598718B (zh) | 设备更换需求分析机构 | |
CN112401934A (zh) | 超声成像方法和*** | |
JP2008036255A (ja) | 超音波診断装置 | |
CN116910758B (zh) | 恶意软件检测方法、装置、电子设备、芯片及存储介质 | |
TWI453404B (zh) | 超音波成像系統及其影像處理方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170926 Termination date: 20201017 |
|
CF01 | Termination of patent right due to non-payment of annual fee |