CN104287777A - 超音波扫描方法及超音波扫描系统 - Google Patents

超音波扫描方法及超音波扫描系统 Download PDF

Info

Publication number
CN104287777A
CN104287777A CN201410553268.8A CN201410553268A CN104287777A CN 104287777 A CN104287777 A CN 104287777A CN 201410553268 A CN201410553268 A CN 201410553268A CN 104287777 A CN104287777 A CN 104287777A
Authority
CN
China
Prior art keywords
image
ultrasound
binaryzation
block
input image
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
Application number
CN201410553268.8A
Other languages
English (en)
Other versions
CN104287777B (zh
Inventor
陈翔洲
郭富彦
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Qisda Suzhou Co Ltd
Qisda Corp
Original Assignee
Qisda Suzhou Co Ltd
Qisda Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Qisda Suzhou Co Ltd, Qisda Corp filed Critical Qisda Suzhou Co Ltd
Priority to CN201410553268.8A priority Critical patent/CN104287777B/zh
Publication of CN104287777A publication Critical patent/CN104287777A/zh
Application granted granted Critical
Publication of CN104287777B publication Critical patent/CN104287777B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8995Combining images from different aspect angles, e.g. spatial compounding
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5238Devices 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52046Techniques for image enhancement involving transmitter or receiver
    • G01S7/52047Techniques for image enhancement involving transmitter or receiver for elimination of side lobes or of grating lobes; for increasing resolving power

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可为超音波影像中常见的旁瓣(side lobe)、假影(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 (10)

1.一种超音波扫描方法,其特征在于,包含下列步骤:
对标的物朝N个不同角度发出N组超音波讯号且接收自该标的物反射及/或散射的该N组超音波讯号,其中每一组该超音波讯号分别对应该N个不同角度的其中之一,且N为大于1的正整数;
将该N组超音波讯号转换为N张超音波输入影像;
对该N张超音波输入影像执行二值化演算法,以得到N张二值化影像;
对该N张二值化影像进行比较程序,以判断该N张二值化影像的至少其中之一是否存在杂讯;
当该N张二值化影像中的第i张二值化影像存在该杂讯时,将该杂讯自对应该第i张二值化影像的第i张超音波输入影像中滤除,其中i为小于或等于N的正整数;以及
将该N张超音波输入影像合成为超音波输出影像。
2.如权利要求1所述的超音波扫描方法,其特征在于,对该N张超音波输入影像执行二值化演算法,以得到N张二值化影像的步骤包含:
选取该N张超音波输入影像的其中之一作为参考影像;
将每一张该超音波输入影像分割成对应的M个区块,其中M为正整数;
将该参考影像中各该区块的所有像素值平均,以得到M个二值化阀值;以及
以对应各该区块的该二值化阀值对该N张超音波输入影像中对应的各该区块执行该二值化演算法,以得到该N张二值化影像;
其中,该比较程序包含下列步骤:
将除了该参考影像外的每一张该超音波输入影像的该二值化影像中的M个区块分别与该参考影像的该二值化影像中对应的M个区块进行比较;以及
当两对应区块的像素差异量大于杂讯滤除阀值时,判断非该参考影像的该超音波输入影像中的该区块存在该杂讯。
3.如权利要求2所述的超音波扫描方法,其特征在于,对应该参考影像的该超音波讯号是朝该N个不同角度其中的最小角度发出。
4.如权利要求3所述的超音波扫描方法,其特征在于,该最小角度为0度。
5.如权利要求1所述的超音波扫描方法,其特征在于,对该N张超音波输入影像执行二值化演算法,以得到N张二值化影像的步骤包含:
将每一张该超音波输入影像分割成对应的M个区块,其中M为正整数;
将该N张超音波输入影像中对应的各该区块的所有像素值平均,以得到M个二值化阀值;以及
以对应各该区块的该二值化阀值对该N张超音波输入影像中对应的各该区块执行该二值化演算法,以得到该N张二值化影像;
其中,该比较程序包含下列步骤:
将该N张二值化影像中对应的该M个区块相互进行比较;以及
当该M个区块中的特定区块与其它对应的区块的像素差异量皆大于杂讯滤除阀值时,判断该特定区块存在该杂讯。
6.一种超音波扫描系统,其特征在于,包含:
超音波探头,用以对标的物朝N个不同角度发出N组超音波讯号且接收自该标的物反射及/或散射的该N组超音波讯号,其中每一组该超音波讯号分别对应该N个不同角度的其中之一,且N为大于1的正整数;以及
超音波扫描仪,用以与该超音波探头形成通讯,该超音波扫描仪包含:
转换单元,用以将该N组超音波讯号转换为N张超音波输入影像;
二值化单元,用以对该N张超音波输入影像执行二值化演算法,以得到N张二值化影像;
比较单元,用以对该N张二值化影像进行比较程序,以判断该N张二值化影像的至少其中之一是否存在杂讯;
杂讯滤除单元,用以当该N张二值化影像中的第i张二值化影像存在该杂讯时,将该杂讯自对应该第i张二值化影像的第i张超音波输入影像中滤除,其中i为小于或等于N的正整数;以及
影像合成单元,用以将该N张超音波输入影像合成为超音波输出影像。
7.如权利要求6所述的超音波扫描系统,其特征在于,该二值化单元选取该N张超音波输入影像的其中之一作为参考影像,将每一张该超音波输入影像分割成对应的M个区块,将该参考影像中各该区块的所有像素值平均,以得到M个二值化阀值,且以对应各该区块的该二值化阀值对该N张超音波输入影像中对应的各该区块执行该二值化演算法,以得到该N张二值化影像,其中M为正整数;
该比较单元将除了该参考影像外的每一张该超音波输入影像的该二值化影像中的M个区块分别与该参考影像的该二值化影像中对应的M个区块进行比较;以及
当两对应区块的像素差异量大于杂讯滤除阀值时,该比较单元判断非该参考影像的该超音波输入影像中的该区块存在该杂讯,以判断该N张二值化影像的至少其中之一是否存在该杂讯。
8.如权利要求7所述的超音波扫描系统,其特征在于,对应该参考影像的该超音波讯号是朝该N个不同角度其中的最小角度发出。
9.如权利要求8所述的超音波扫描系统,其特征在于,该最小角度为0度。
10.如权利要求6所述的超音波扫描系统,其特征在于,该二值化单元将每一张该超音波输入影像分割成对应的M个区块,将该N张超音波输入影像中对应的各该区块的所有像素值平均,以得到M个二值化阀值,且以对应各该区块的该二值化阀值对该N张超音波输入影像中对应的各该区块执行该二值化演算法,以得到该N张二值化影像,其中M为正整数;
该比较单元将该N张二值化影像中对应的该M个区块相互进行比较;以及
当该M个区块中的特定区块与其它对应的区块的像素差异量皆大于杂讯滤除阀值时,该比较单元判断该特定区块存在该杂讯,以判断该N张二值化影像的至少其中之一是否存在该杂讯。
CN201410553268.8A 2014-10-17 2014-10-17 超音波扫描方法及超音波扫描系统 Expired - Fee Related CN104287777B (zh)

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 true CN104287777A (zh) 2015-01-21
CN104287777B 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)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105806473A (zh) * 2016-03-14 2016-07-27 苏州佳世达电通有限公司 超音波探头的参数设定方法
TWI644526B (zh) * 2017-02-21 2018-12-11 佳世達科技股份有限公司 超音波系統及雜訊消除方法
CN113311440A (zh) * 2020-02-26 2021-08-27 苏州佳世达电通有限公司 用于改善超声波影像质量的方法及其探测系统

Citations (11)

* Cited by examiner, † Cited by third party
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 血流情報画像表示方法および超音波診断装置
CN101102721A (zh) * 2005-02-09 2008-01-09 株式会社日立医药 超声波诊断装置以及超声波摄像方法
JP4473981B2 (ja) * 1999-07-21 2010-06-02 株式会社日立メディコ 超音波診断装置
CN101854537A (zh) * 2009-03-31 2010-10-06 深圳迈瑞生物医疗电子股份有限公司 一种超声图像数据优化及其造影定量分析方法和系统
CN102247172A (zh) * 2010-02-08 2011-11-23 通用电气公司 自动化胎儿胎龄评估的系统和方法
US20110301464A1 (en) * 2004-08-27 2011-12-08 Nanyang Technological University Home ultrasound system
US20120184845A1 (en) * 2010-11-11 2012-07-19 University Of Pittsburgh - Of The Commonwealth System Of Higher Education Automated macular pathology diagnosis in threedimensional (3d) spectral domain optical coherence tomography (sd-oct) images
WO2013098696A1 (en) * 2011-12-30 2013-07-04 Koninklijke Philips Electronics N.V. Method and apparatus for needle visualization enhancement in ultrasound imaging
CN103908300A (zh) * 2013-01-03 2014-07-09 西门子公司 诊断超声成像中的针增强
CN104039235A (zh) * 2012-01-11 2014-09-10 日立阿洛卡医疗株式会社 超声波诊断装置

Patent Citations (11)

* Cited by examiner, † Cited by third party
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 株式会社日立メディコ 超音波診断装置
US20110301464A1 (en) * 2004-08-27 2011-12-08 Nanyang Technological University Home ultrasound system
CN101102721A (zh) * 2005-02-09 2008-01-09 株式会社日立医药 超声波诊断装置以及超声波摄像方法
CN101854537A (zh) * 2009-03-31 2010-10-06 深圳迈瑞生物医疗电子股份有限公司 一种超声图像数据优化及其造影定量分析方法和系统
CN102247172A (zh) * 2010-02-08 2011-11-23 通用电气公司 自动化胎儿胎龄评估的系统和方法
US20120184845A1 (en) * 2010-11-11 2012-07-19 University Of Pittsburgh - Of The Commonwealth System Of Higher Education Automated macular pathology diagnosis in threedimensional (3d) spectral domain optical coherence tomography (sd-oct) images
WO2013098696A1 (en) * 2011-12-30 2013-07-04 Koninklijke Philips Electronics N.V. Method and apparatus for needle visualization enhancement in ultrasound imaging
CN104039235A (zh) * 2012-01-11 2014-09-10 日立阿洛卡医疗株式会社 超声波诊断装置
CN103908300A (zh) * 2013-01-03 2014-07-09 西门子公司 诊断超声成像中的针增强

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105806473A (zh) * 2016-03-14 2016-07-27 苏州佳世达电通有限公司 超音波探头的参数设定方法
CN105806473B (zh) * 2016-03-14 2018-12-11 苏州佳世达电通有限公司 超音波探头的参数设定方法
TWI644526B (zh) * 2017-02-21 2018-12-11 佳世達科技股份有限公司 超音波系統及雜訊消除方法
US10813627B2 (en) 2017-02-21 2020-10-27 Qisda Corporation Ultrasound system and noise eliminating method
CN113311440A (zh) * 2020-02-26 2021-08-27 苏州佳世达电通有限公司 用于改善超声波影像质量的方法及其探测系统

Also Published As

Publication number Publication date
CN104287777B (zh) 2017-09-26

Similar Documents

Publication Publication Date Title
US8600122B2 (en) Apparatus and method for culling substantially redundant data in fingerprint sensing circuits
CN111247525A (zh) 一种车道检测方法、装置及车道检测设备、移动平台
US8717228B2 (en) Method and device for detecting target object, and radar apparatus
CN104287777A (zh) 超音波扫描方法及超音波扫描系统
US11372087B2 (en) Radar system and method for receiving and analyzing radar signals
Rigau et al. Medical image segmentation based on mutual information maximization
CN102540160B (zh) 信号处理装置、雷达装置、信号处理方法
US11657485B2 (en) Method for expanding image depth and electronic device
WO2023071992A1 (zh) 多传感器信号融合的方法、装置、电子设备及存储介质
US20190282205A1 (en) Ultrasound imaging system and ultrasound imaging method
CN102034109A (zh) 一种基于统计特性的合成孔径声纳目标检测方法
JP2019074478A (ja) 識別装置、識別方法およびプログラム
CN102540147A (zh) 信号处理装置、雷达装置、信号处理方法及信号处理程序
US9349069B2 (en) Dynamic line-detection system for processors having limited internal memory
US20100128975A1 (en) Method and system for object recognition based on a trainable dynamic system
TWI521225B (zh) 超音波掃描方法及超音波掃描系統
El Hajjouji et al. FPGA based real-time lane detection and tracking implementation
US9916663B2 (en) Image processing method and process simulation apparatus
CN115439509B (zh) 一种多目标跟踪方法、装置、计算机设备及存储介质
CN103051829A (zh) 基于fpga平台的原始图像数据降噪系统及降噪方法
CN114581944B (zh) 一种毫米波图像处理方法、装置及电子设备
Geninatti et al. Real-Time Reconfigurable Processor to Detect Similarities in Compressed Video Using Generalized Hough Transformation
CN108389174B (zh) 超音波成像系统及超音波成像方法
CN115407281A (zh) 点云去噪方法、装置、探测雷达和存储介质
JP4878239B2 (ja) 超音波診断装置

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