CN106175839A - 用于超声相控阵的扇形扫描成像的扫描检测方法和装置 - Google Patents
用于超声相控阵的扇形扫描成像的扫描检测方法和装置 Download PDFInfo
- Publication number
- CN106175839A CN106175839A CN201510219036.3A CN201510219036A CN106175839A CN 106175839 A CN106175839 A CN 106175839A CN 201510219036 A CN201510219036 A CN 201510219036A CN 106175839 A CN106175839 A CN 106175839A
- Authority
- CN
- China
- Prior art keywords
- point
- interpolation
- coordinate
- scanning
- scan line
- 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
Landscapes
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
一种用于超声相控阵的扇形扫描成像的扫描检测方法和装置,该方法包括:坐标变换步骤,对扫描区域内多个采样点的极坐标数据进行坐标变换,得到成像坐标系上对应像素点的每个该采样点的直角坐标数据,并存储该直角坐标数据;对坐标转换后的扫描线进行插值步骤,通过判断待插值点与其周围的采样点的直线位置关系,找到与该待插值点距离最近的四个回波信号采样点,根据该待插值点与该四个回波信号采样点的距离关系,按权重进行插值;以及存储该插值结果。本发明还提供了用于超声相控阵的扇形扫描成像的扫描检测装置,以及该装置的使用方法。
Description
技术领域
本发明涉及一种超声相控阵检测技术,特别是一种适用于等声程、等深度、等跨度三种聚焦形式,也同样适用于线性扫描及任意聚焦点组合的基于斜探头的扇形扫描成像的扫描检测方法和装置。
背景技术
相控阵超声检测技术是通过控制阵列探头各个阵元的接收与发射的延迟时间,形成合成声束的聚焦、扫描等,从而实现超声波束的偏转、聚焦等各种扫描效果,最终在扫描范围内实现高分辨率的超声缺陷成像。在相控阵超声发射状态下,阵列换能器中各阵元按一定延迟规律顺序激发,产生的超声发射子波束在空间合成,形成聚焦点和指向性。改变各阵元激发的延迟规律,可以改变焦点位置的波束指向,形成在一定空间范围内的扫描聚焦。和常规的扫描成像相比,相控阵超声成像技术采用阵列式的探头,只需改变探头中各个阵元激励的延迟时序,无需移动探头就可以方便地实现全方位和多角度扫查。
在相控阵超声接收状态下,列阵换能器的各个阵元接收回波信号,按不同延时值进行延时,然后加权求和作为输出。通过设定一定的延时规律,可以实现对声场中的指定物点进行聚焦接收,采用不同的延时规律,即可实现对不同点和不同方向上的接收聚焦和扫描。例如申请号为“200580020418.8”,名称为“阵列式超声换能器”的中国发明专利申请,申请号为“200580028604.6”,名称为“二维超声换能器阵列”的中国发明专利申请,申请号为“200680029122.7”,名称为“超声换能器阵列”的中国发明专利申请所公开的各种换能器,均为应用超声成像,进行非侵入地探查固、液、气体以进行相应检测。
现有技术的超声检测过程需要用楔块调整检测的范围,一般情况下超声相控阵可以通过改变延迟时间的方法来改变检测角度,但是实际测试中由于超声旁瓣等的干扰,当相控阵偏转角度过大时就会影响检测结果,因此常常在超声检测中使用楔块,通过加楔块增加超声相控阵的扫查范围。
参见图1,图1为超声扫描原理图。如图所示,设要扫描的范围是β1~βn,间隔角度为Δβ°,则共有N=(βn-β1)/Δβ条扫描线。根据Snell定律(公式1),反射角不相等,对应的入射角也不相等,而这N条扫描线都是从阵列中心点发射出来的,因此,这N条扫描线在楔块1与待测工件2的交界处会产生N个不同的折射点。这种现象即为入射点漂移。
在相控阵超声发射状态下,阵列换能器中各阵元按一定延迟规律顺序激发,产生的超声发射子波束在空间合成,形成聚焦点和指向性。由于入射点漂移,N条扫描线不是以同一点进入工件中,所以成像后的扇扫图形也不是以一点为圆心的扇形,对检测精度造成一定影响。
发明内容
本发明所要解决的技术问题是针对现有技术的上述问题,提供一种适用于等声程、等深度、等跨度三种聚焦形式,也适用于线性扫描及任意聚焦点组合的用于超声相控阵的扇形扫描成像的扫描检测方法和装置。
为了实现上述目的,本发明提供了一种用于超声相控阵的扇形扫描成像的扫描检测方法,其中,包括如下步骤:
S100、坐标变换,对扫描区域内多个采样点的极坐标数据进行坐标变换,得到成像坐标系上对应像素点的每个该采样点的直角坐标数据,并存储该直角坐标数据;
S200、对坐标转换后的扫描线进行插值,通过判断待插值点与其周围的采样点的直线位置关系,找到与该待插值点距离最近的四个回波信号采样点,根据该待插值点与该四个回波信号采样点的距离关系,按权重进行插值;以及
S300、存储该插值结果。
上述的超声相控阵的扇形扫描成像的扫描检测方法,其中,所述步骤S100进一步包括:
S101、生成第一数据存储表,根据扫描范围内N条扫描线角度、声程方向采样点个数M和直角坐标系上的图像分辨率,计算出M*N个采样点的极坐标对应该直角坐标系的直角坐标,并存储在第一数据存储表中;
S102、生成第二数据存储表,根据该声程方向的每个该采样点的直角坐标,计算出每条扫描线相邻两点的直线方程,并将该直线方程对应的直线的斜率和该直线在Y轴上的截距存储在第二数据存储表中;
S103、生成第三数据存储表,根据在同一声程上的N个采样点坐标,计算相邻两个扫描点的直线方程,将该直线方程对应的直线的斜率和该直线在Y轴上的截距存储在第三数据存储表中。
上述的超声相控阵的扇形扫描成像的扫描检测方法,其中,所述步骤S101进一步包括:
设声波在楔块中传播速度为c1,在工件中传播速度为c2,阵列中心点的高度为h,声波入射角依次为α1~αn,折射角依次为β1~βn,则极坐标系中各入射点位置为:
将以不同入射点为起始点的采样点转换到直角坐标上,转换公式如下:
上述的超声相控阵的扇形扫描成像的扫描检测方法,其中,所述步骤S200进一步包括:
S201、邻近点查找,通过判断待插值点与临近采样点确定的直线方程的位置关系,确定待插值点的四个邻近采样点坐标值;
S202,插值,利用该四个邻近采样点对该待插值点进行线性插值。
上述的超声相控阵的扇形扫描成像的扫描检测方法,其中,所述步骤S201进一步包括:
S2011、确定该四个邻近采样点的纵坐标,遍历直角坐标系上各个待插值点P,从第一条声程方向的扫描线开始,判断点P与扫描线上采样点确定的直线的位置关系,如果点P在扫描线上方,继续判断下一条扫描线,直到点P位于当前扫描线的下方为止,则确定点P位于当前扫描线与前一扫描线之间,由此获得点P的最邻近四个采样点的纵坐标;
S2012、确定该四个邻近采样点的横坐标,从已找到的第n条线和第n+1条扫描线中,判断点P与声程方向的采样点确定的直线位置关系,如果点P在第一条直线下方,继续判断下一条直线,直到点P位于当前直线的上方位置,则确定点P位于当前直线与前一条直线之间,由此获得点P的最邻近四个采样点的横坐标。
上述的超声相控阵的扇形扫描成像的扫描检测方法,其中,在步骤S2011中,如果点P位于第一条扫描线下方或最后一条扫描线上方,则不予以插值。
上述的超声相控阵的扇形扫描成像的扫描检测方法,其中,在步骤S2012中,如果点P位于第一条直线的上方或最后一条直线的下方,则不予以插值。
上述的超声相控阵的扇形扫描成像的扫描检测方法,其中,所述步骤S202进一步包括:
S2021、找到该待插值点两侧线段上满足于一设定比例关系且与该待插值点在同一直线上的两点,并对该两点进行线性插值;
S2022、利用该两点的灰度值对该待插值点进行线性插值,或者在角度方向按探头指向性函数对该待插值点进行插值。
为了更好地实现上述目的,本发明还提供了一种用于超声相控阵的扇形扫描成像的扫描检测装置,其中,包括:
计算模块,用于对扫描区域内多个采样点的极坐标数据进行坐标变换,得到成像坐标系上对应像素点的每个该采样点的直角坐标数据,并存储该直角坐标数据;并且,对坐标转换后的扫描线进行插值,通过判断待插值点与其周围的采样点的直线位置关系,找到与该待插值点距离最近的四个回波信号采样点,根据该待插值点与该四个回波信号采样点的距离关系,按权重进行插值;
存储模块,与所述计算模块连接,用于存储所述计算模块的坐标变换结果及插值结果;
输出显示模块,与所述存储模块连接,用于输出并显示计算或检测结果;以及
数据输入模块,与所述计算模块连接,用于扫描检测待测物体并将扫描数据传输给所述计算模块。
为了更好地实现上述目的,本发明还提供了一种上述超声相控阵的扇形扫描成像的扫描检测装置的使用方法。
本发明的技术效果在于:
本发明用于斜探头超声相控阵成像的扫描转换,能把以极坐标方式存储的图像数据即回波信号采样点,以直角坐标的方式准确、直观地显示在显示设备上。
以下结合附图和具体实施例对本发明进行详细描述,但不作为对本发明的限定。
附图说明
图1为超声扫描原理图;
图2为本发明的检测原理示意图;
图3为本发明一实施例的方法流程图;
图4为本发明另一实施例的方法流程图;
图5为本发明一实施例的装置结构框图。
其中,附图标记
1 楔块
2 待测工件
3 阵列换能器
4 计算模块
5 存储模块
6 输出显示模块
7 数据输入模块
A 折射点
S100-S300 步骤
具体实施方式
下面结合附图对本发明的结构原理和工作原理作具体的描述:
扇形扫描显示是一种灰度调制显示,它将A扫描(Amplitude scan)显示的回波幅度进行灰度调制。通过改变探头的角度,从而使超声指向快速变化,使每隔一定小角度,被探测方向不同深度所有界面的反射回波,都以灰度点的形式显示在对应的扫描线上,形成一条A扫,所有A扫可形成一幅由探头探测角度和声程决定的二维超声断层图像,成为扇形扫描图像。
在相控阵超声发射状态下,阵列换能器3中各阵元按一定延迟规律顺序激发,产生的超声发射子波束在空间合成,形成聚焦点和指向性。改变各阵元激发的延迟规律,可以改变焦点位置的波束指向,形成在一定空间范围内的扫描聚焦。在相控阵超声接收状态下,列阵换能器的各个阵元接收回波信号,按不同延时值进行延时,然后加权求和作为输出。通过设定一定的延时规律,可以实现对声场中的指定物点进行聚焦接收,采用不同的延时规律,即可实现对不同点和不同方向上的接收聚焦和扫描。
超声检测过程需要用楔块1调整检测的范围,一般情况下超声相控阵可以通过改变延迟时间的方法来改变检测角度,但是实际测试中由于超声旁瓣等的干扰,当相控阵偏转角度过大时就会影响检测结果,因此常常在超声检测中使用楔块1。通过楔块1增加超声相控阵的扫查范围。相控阵超声通过斜探头扫描形成的区域为扇形,相应的成像区域也为扇形。由于检测时探头是以极坐标方式扫查的,回波信号的存储也是以极坐标的方式存储,然而在显示器端,图像是以直角坐标方式显示的,这样就需要将以极坐标方式存储的数据通过坐标转换以直角坐标的方式存储。
参见图2,图2为本发明的检测原理示意图。该图显示为全部在待测工件2中,A为折射点,上边缘即为待测工件2上表面,即图1待测工件2中的部分,也就是所谓的“二维超声断层”。图中的每条斜线即为声程方向,该声程上的信号通过AD转换,导入计算机就是A扫。其中,θ角具体取值为折射角,α角具体取值为折射角,图中每个竖线与横线交叉点为灰度点。如图2所示,在坐标转换过程中,由于转换后的扫描线落在直角坐标时,在临近的两条扫描线之间的距离随着声程方向的增大而逐渐增大,不可能覆盖到直角坐标系的所有像素点上;并且,经过坐标转换后的采样点在直角坐标系的横坐标与纵坐标也不一定是整数,相邻的两个采样点的距离也会因为图像比例改变的原因,在它们之间会有没有覆盖到的像素点,这样就需要对坐标转换后的扫描线进行插值。
参见图3,图3为本发明一实施例的方法流程图。本发明的用于超声相控阵的扇形扫描成像的扫描检测方法,包括如下步骤:
步骤S100、坐标变换,对扫描区域内多个采样点的极坐标数据进行坐标变换,得到成像坐标系上对应像素点的每个该采样点的直角坐标数据,并存储该直角坐标数据;
步骤S200、对坐标转换后的扫描线进行插值,通过判断待插值点与其周围的采样点的直线位置关系,找到与该待插值点距离最近的四个回波信号采样点,根据该待插值点与该四个回波信号采样点的距离关系,按权重进行插值;
步骤S300、存储该插值结果。
参见图4,图4为本发明另一实施例的方法流程图。其中,所述步骤S100还进一步求同声程角度方向上的采样点确定的直线方程,或求声程方向上采样点确定的直线方程,具体可包括:
步骤S101、生成第一数据存储表,根据扫描范围内N条扫描线角度、声程方向采样点个数M和直角坐标系上的图像分辨率,计算出M*N个采样点的极坐标对应该直角坐标系的直角坐标,并存储在第一数据存储表中;
步骤S102、生成第二数据存储表,根据该声程方向的每个该采样点的直角坐标,计算出每条扫描线相邻两点的直线方程,并将该直线方程对应的直线的斜率和该直线在Y轴上的截距存储在第二数据存储表中;
步骤S103、生成第三数据存储表,根据在同一声程上的N个采样点坐标,计算相邻两个扫描点的直线方程,将该直线方程对应的直线的斜率和该直线在Y轴上的截距存储在第三数据存储表中。
其中,所述步骤S101进一步包括:
参见图2,设声波在楔块1中传播速度为c1,在工件2中传播速度为c2,阵列中心点的高度为h,声波入射角依次为α1~αn,折射角依次为β1~βn,则极坐标系中各入射点位置为:
其中,x为每条扫描线在工件2上入射点的x坐标。
将以不同入射点为起始点的采样点转换到直角坐标上,转换公式如下:
其中,x,y为每条扫描线上的采样点的直角坐标,r,θ为相应采样点的极坐标。
其中,所述步骤S200进一步包括对于直角坐标系上任意一个待插值点,判断该点是否需要插值,若判断结果为否,则直接对该点的灰度值赋值为0;若判断结果为是,则判断该点位于哪两条扫描线之间,然后判断该点位于声程方向的哪两条直线之间,并根据确定的四个点的坐标在两侧的扫描线上找到符合条件的两个点,分别给符合条件的两个点插值,最后根据待插值点与找到的两个点的距离进行插值。具体可包括:
步骤S201、邻近点查找,通过判断待插值点与临近采样点确定的直线方程的位置关系,确定待插值点的四个邻近采样点坐标值;
步骤S202,插值,利用该四个邻近采样点对该待插值点进行线性插值。
其中,所述步骤S201进一步包括:
步骤S2011、确定该四个邻近采样点的纵坐标,遍历直角坐标系上各个待插值点P,从第一条声程方向的扫描线开始,判断点P与扫描线上采样点确定的直线的位置关系,如果点P在扫描线上方,继续判断下一条扫描线,直到点P位于当前扫描线的下方为止,则确定点P位于当前扫描线与前一扫描线之间,由此获得点P的最邻近四个采样点的纵坐标;
步骤S2012、确定该四个邻近采样点的横坐标,从已找到的第n条线和第n+1条扫描线中,判断点P与声程方向的采样点确定的直线位置关系,如果点P在第一条直线下方,继续判断下一条直线,直到点P位于当前直线的上方位置,则确定点P位于当前直线与前一条直线之间,由此获得点P的最邻近四个采样点的横坐标。
其中,在步骤S2011中,如果点P位于第一条扫描线下方或最后一条扫描线上方,则不予以插值。在步骤S2012中,如果点P位于第一条直线的上方或最后一条直线的下方,则不予以插值。
其中,所述步骤S202进一步包括:
S2021、找到该待插值点两侧线段上满足于一设定比例关系且与该待插值点在同一直线上的两点,并对该两点进行线性插值;
S2022、利用该两点的灰度值对该待插值点进行线性插值,或者在角度方向按探头指向性函数对该待插值点进行插值。
参见图5,图5为本发明一实施例的装置结构框图。本发明的用于超声相控阵的扇形扫描成像的扫描检测装置,包括:
计算模块4,用于对扫描区域内多个采样点的极坐标数据进行坐标变换,得到成像坐标系上对应像素点的每个该采样点的直角坐标数据,并存储该直角坐标数据;并且,对坐标转换后的扫描线进行插值,通过判断待插值点与其周围的采样点的直线位置关系,找到与该待插值点距离最近的四个回波信号采样点,根据该待插值点与该四个回波信号采样点的距离关系,按权重进行插值;
存储模块5,与所述计算模块4连接,用于存储所述计算模块4的坐标变换结果及插值结果;
输出显示模块6,与所述存储模块5连接,用于输出并显示计算或检测结果;以及
数据输入模块7,与所述计算模块4连接,用于扫描检测待测物体并将扫描数据传输给所述计算模块4。
本发明还提供了上述超声相控阵的扇形扫描成像的扫描检测方法或装置的使用方法。具体包括,第一步,针对具体的待测工件2,检测人员设定好扫描线的角度和检测深度,计算模块4计算所有采样点的极坐标和直角坐标,并存入存储模块5;第二步,计算模块4计算显示器像素点与采样点之间的关系,算出像素点与采样点信号之间的插值权重关系,并存入存储模块5。第三步,数据输入模块7得到各采样点的采样数据存入存储模块5。第四步,计算模块4根据各采样点的采样数据插值得到像素点显示数据,存入存储模块5,同时输出到输出显示模块6,输出显示模块6显示检测结果,至此完成检测过程。
下面以具体数据对本发明的超声相控阵的扇形扫描成像的扫描检测方法予以进一步说明:
首先,根据扫描范围内N条扫描线角度、声程方向采样点个数M、建立在直角坐标系上的图像分辨率,计算出M*N个点在直角坐标系上对应的坐标,存储在第一数据存储表中。计算公式如下:
1)如图2,声波在楔块1中传播速度为c1,在工件2中传播速度为c2,阵列中心点的高度为h,声波入射角依次为α1~βn,折射角依次为β1~βn,则极坐标系中各入射点位置是:
其中,x为每条扫描线在待测工件2上入射点的x坐标。
2)通过公式(3)将以不同入射点为起始点的采样点转换到直角坐标系:
其中,x,y为每条扫描线上的采样点的直角坐标,r,θ为相应采样点的极坐标。
其次,根据声程方向的采样点坐标,计算出每条扫描线的相邻两点确定的直线的方程,将直线的斜率,和直线在Y轴上的截距存储在第二数据存储表中。由于一条扫描线上的所有采样点在一条直线上,所以每条直线只需用前两个点的坐标就可以求出这条扫描线上所有相邻点确定的直线方程,共N个方程;
然后,根据在同一声程上的N个采样点坐标,计算相邻两个扫描点确定的直线的方程,将直线的斜率,和直线在Y轴上的截距存储在第三数据存储表中。同一声程上的采样点,共(N-1)*M个方程;
本发明的插值方法是找到距待插值点最近的四个点,采用线性插值的思想,根据待插值点与周围四个点的位置大小关系分配权重,将周围四个点的灰度值按权重分配到待插值点中。
具体插值步骤如下:
遍历直角坐标系上各个待插值点,并进行如下步骤,参见图2,对图2中的点P(j,i)(假设该点P(10,13))进行插值。
1、从第一条声程方向的扫描线开始,判断点P与扫描线上采样点确定的直线的位置关系,如果点P在扫描线上方,继续判断下一条扫描线,直到点P位于当前扫描线的下方为止,则确定点P位于当前扫描线与前一扫描线之间。由此获得点P的最邻近四个点的纵坐标。如果点P位于第一条扫描线下方或最后一条扫描线上方,则不予以插值。设直线方程为y=k×x+b判断公式如下:
公式(4)
图2中,n条直线方程为:y=x+1;
n+1条直线方程为:y=1.3×x+1.7,带入i=13,j=10;
i(13)≥1×j(10)+1且i(13)<1.3×j(10)+1.7,
所以判断P在直线n与n+1之间。
2、从已找到的第n条线和第n+1条扫描线中,判断点P与声程方向的采样点确定的直线的位置关系,如果点P在第一条直线下方,继续判断下一条直线,直到点P位于当前直线的上方位置,则确定点P位于当前直线与前一条直线之间。由此获得点P的最邻近四个点的横坐标。如果点P位于第一条直线的上方或最后一条直线的下方,则不予以插值。判断公式参见公式(4)。
如此找到a点(9.8,10.8),b点(11.7,12.7),c点(8.5,12.75),d点(10.3,15.09),其灰度值分别设为1,2,3,4;
3、找到待插值点P两侧扫描线方向上且满足如下两个条件的两个点,如图2中的点e(x1,y1)和点f(x2,y2):
1)点a和点e的距离与点b和点e的距离之比等于点c和点f的距离与点d和点f的距离之比;
(x1-9.8)/(10.7-x1)=(x2-8.5)/(10.3-x2);
2)点e、点f、点P在同一直线上。
则:(10-x1)/(13-y1)=(x2-10)/(y2-13)
求得e点(10.76,11.76),f点(9.42,13.94)
4、点P的灰度值可基于以下原则获得:
1)点a和点e的距离除以点a和点b的距离之比(0.51)再乘以点b的灰度值(2),点b和点e的距离除以点a和点b的距离之比(0.49)乘以点a的灰度值(1),将两项值相加,求得点e的灰度值0.51×2+0.49×1=1.51;
2)点c和点f的距离除以点c和点d的距离之比(0.51)再乘以点d的灰度值(4),点d和点f的距离除以点c和点d的距离之比(0.49)乘以点c的灰度值(3),将两项值相加,求得点f的灰度值0.51×4+0.49×3=3.51;
3)点e和点P的距离除以点e和点f的距离之比(0.57)再乘以点f的灰度值(3.51),点f和点P的距离除以点e和点f的距离之比(0.43)乘以点e的灰度值(1.51),将两项值相加,求得点P的灰度值0.57×3.51+0.43×1.51=2.65。
本发明用于斜探头超声相控阵成像的扫描转换,具体包括坐标变换、采样点坐标及相关数据存储、插值等步骤。先将扫描坐标系(极坐标系)的采样点进行坐标变换,得到成像坐标系(直角坐标系)上的坐标值,能把以极坐标方式存储的图像数据即回波信号采样点,以直角坐标的方式准确、直观地显示在显示设备上。
当然,本发明还可有其它多种实施例,在不背离本发明精神及其实质的情况下,熟悉本领域的技术人员当可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。
Claims (10)
1.一种用于超声相控阵的扇形扫描成像的扫描检测方法,其特征在于,包括如下步骤:
S100、坐标变换,对扫描区域内多个采样点的极坐标数据进行坐标变换,得到成像坐标系上对应像素点的每个该采样点的直角坐标数据,并存储该直角坐标数据;
S200、对坐标转换后的扫描线进行插值,通过判断待插值点与其周围的采样点的直线位置关系,找到与该待插值点距离最近的四个回波信号采样点,根据该待插值点与该四个回波信号采样点的距离关系,按权重进行插值;以及
S300、存储该插值结果。
2.如权利要求1所述的超声相控阵的扇形扫描成像的扫描检测方法,其特征在于,所述步骤S100进一步包括:
S101、生成第一数据存储表,根据扫描范围内N条扫描线角度、声程方向采样点个数M和直角坐标系上的图像分辨率,计算出M*N个采样点的极坐标对应该直角坐标系的直角坐标,并存储在第一数据存储表中;
S102、生成第二数据存储表,根据该声程方向的每个该采样点的直角坐标,计算出每条扫描线相邻两点的直线方程,并将该直线方程对应的直线的斜率和该直线在Y轴上的截距存储在第二数据存储表中;
S103、生成第三数据存储表,根据在同一声程上的N个采样点坐标,计算相邻两个扫描点的直线方程,将该直线方程对应的直线的斜率和该直线在Y轴上的截距存储在第三数据存储表中。
3.如权利要求2所述的超声相控阵的扇形扫描成像的扫描检测方法,其特征在于,所述步骤S101进一步包括:
设声波在楔块中传播速度为c1,在工件中传播速度为c2,阵列中心点的高度为h,声波入射角依次为α1~αn,折射角依次为β1~βn,则极坐标系中各入射点位置为:
其中,x为每条扫描线在待测工件上入射点的x坐标;
将以不同入射点为起始点的采样点转换到直角坐标上,转换公式如下:
其中,x,y为每条扫描线上的采样点的直角坐标,r,θ为相应采样点的极坐标。
4.如权利要求1、2或3所述的超声相控阵的扇形扫描成像的扫描检测方法,其特征在于,所述步骤S200进一步包括:
S201、邻近点查找,通过判断待插值点与临近采样点确定的直线方程的位置关系,确定待插值点的四个邻近采样点坐标值;
S202,插值,利用该四个邻近采样点对该待插值点进行线性插值。
5.如权利要求4所述的超声相控阵的扇形扫描成像的扫描检测方法,其特征在于,所述步骤S201进一步包括:
S2011、确定该四个邻近采样点的纵坐标,遍历直角坐标系上各个待插值点P,从第一条声程方向的扫描线开始,判断点P与扫描线上采样点确定的直线的位置关系,如果点P在扫描线上方,继续判断下一条扫描线,直到点P位于当前扫描线的下方为止,则确定点P位于当前扫描线与前一扫描线之间,由此获得点P的最邻近四个采样点的纵坐标;
S2012、确定该四个邻近采样点的横坐标,从已找到的第n条线和第n+1条扫描线中,判断点P与声程方向的采样点确定的直线位置关系,如果点P在第一条直线下方,继续判断下一条直线,直到点P位于当前直线的上方位置,则确定点P位于当前直线与前一条直线之间,由此获得点P的最邻近四个采样点的横坐标。
6.如权利要求5所述的超声相控阵的扇形扫描成像的扫描检测方法,其特征在于,在步骤S2011中,如果点P位于第一条扫描线下方或最后一条扫描线上方,则不予以插值。
7.如权利要求5所述的超声相控阵的扇形扫描成像的扫描检测方法,其特征在于,在步骤S2012中,如果点P位于第一条直线的上方或最后一条直线的下方,则不予以插值。
8.如权利要求4所述的超声相控阵的扇形扫描成像的扫描检测方法,其特征在于,所述步骤S202进一步包括:
S2021、找到该待插值点两侧线段上满足于一设定比例关系且与该待插值点在同一直线上的两点,并对该两点进行线性插值;
S2022、利用该两点的灰度值对该待插值点进行线性插值,或者在角度方向按探头指向性函数对该待插值点进行插值。
9.一种用于超声相控阵的扇形扫描成像的扫描检测装置,其特征在于,包括:
计算模块,用于对扫描区域内多个采样点的极坐标数据进行坐标变换,得到成像坐标系上对应像素点的每个该采样点的直角坐标数据,并存储该直角坐标数据;并且,对坐标转换后的扫描线进行插值,通过判断待插值点与其周围的采样点的直线位置关系,找到与该待插值点距离最近的四个回波信号采样点,根据该待插值点与该四个回波信号采样点的距离关系,按权重进行插值;
存储模块,与所述计算模块连接,用于存储所述计算模块的坐标变换结果及插值结果;
输出显示模块,与所述存储模块连接,用于输出并显示计算或检测结果;以及
数据输入模块,与所述计算模块连接,用于扫描检测待测物体并将扫描数据传输给所述计算模块。
10.一种上述权利要求9所述的超声相控阵的扇形扫描成像的扫描检测装置的使用方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510219036.3A CN106175839A (zh) | 2015-04-30 | 2015-04-30 | 用于超声相控阵的扇形扫描成像的扫描检测方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510219036.3A CN106175839A (zh) | 2015-04-30 | 2015-04-30 | 用于超声相控阵的扇形扫描成像的扫描检测方法和装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106175839A true CN106175839A (zh) | 2016-12-07 |
Family
ID=57458567
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510219036.3A Pending CN106175839A (zh) | 2015-04-30 | 2015-04-30 | 用于超声相控阵的扇形扫描成像的扫描检测方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106175839A (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107255674A (zh) * | 2017-06-21 | 2017-10-17 | 景德镇陶瓷大学 | 一种圆柱壳内壁涂层界面缺陷的超声检测扫描成像装置及其应用方法 |
CN107356670A (zh) * | 2017-07-14 | 2017-11-17 | 武汉工程大学 | 一种基于斜入射的超声相控阵焊缝缺陷检测方法 |
CN108428251A (zh) * | 2018-03-09 | 2018-08-21 | 深圳市中捷视科科技有限公司 | 一种基于机器视觉技术激光结构光自动标定方法 |
CN110415316A (zh) * | 2019-07-02 | 2019-11-05 | 中国航空工业集团公司雷华电子技术研究所 | 一种提高气象视频显示准确性和连续性的插值方法 |
CN111175381A (zh) * | 2019-12-09 | 2020-05-19 | 中北大学 | 基于全矩阵数据的复合构件界面快速成像定量检测方法 |
CN112665528A (zh) * | 2020-12-31 | 2021-04-16 | 河南卫华重型机械股份有限公司 | 一种激光扫描三维成像的矫正方法 |
CN113240587A (zh) * | 2021-07-12 | 2021-08-10 | 深圳华声医疗技术股份有限公司 | 超分辨率扫描变换方法、装置、超声设备及存储介质 |
CN113848375A (zh) * | 2021-10-20 | 2021-12-28 | 国网湖南省电力有限公司 | 一种油绝缘变压器内部器件带电检测装置及其应用方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1073791A (zh) * | 1991-12-20 | 1993-06-30 | 菲利浦光灯制造公司 | 图象处理设备 |
CN101390756A (zh) * | 2007-09-20 | 2009-03-25 | 深圳迈瑞生物医疗电子股份有限公司 | 超声成像系统中的扫描线插线处理方法 |
CN101441401A (zh) * | 2007-11-20 | 2009-05-27 | 深圳迈瑞生物医疗电子股份有限公司 | 快速确定超声成像系统中成像区域的方法和装置 |
CN101449983A (zh) * | 2007-11-29 | 2009-06-10 | 深圳迈瑞生物医疗电子股份有限公司 | 超声线阵偏转成像的扫描变换插值方法和装置 |
CN101729893A (zh) * | 2008-08-15 | 2010-06-09 | 北京北大众志微系统科技有限责任公司 | 基于软硬件协同处理的mpeg多格式兼容解码方法及其装置 |
CN102025960A (zh) * | 2010-12-07 | 2011-04-20 | 浙江大学 | 一种基于自适应插值的运动补偿去隔行方法 |
CN102135615A (zh) * | 2010-03-22 | 2011-07-27 | 上海通号轨道交通工程技术研究中心有限公司 | 一种测距定标方法 |
CN102222366A (zh) * | 2011-07-29 | 2011-10-19 | 电子科技大学 | 复杂空间曲面的拟合方法 |
CN103837608A (zh) * | 2014-03-12 | 2014-06-04 | 深圳市神视检验有限公司 | 一种相控阵接收动态聚焦补偿方法及系统 |
CN104359982A (zh) * | 2014-11-19 | 2015-02-18 | 哈尔滨工业大学 | 一种薄板焊缝超声相控阵检测楔块的方法 |
-
2015
- 2015-04-30 CN CN201510219036.3A patent/CN106175839A/zh active Pending
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1073791A (zh) * | 1991-12-20 | 1993-06-30 | 菲利浦光灯制造公司 | 图象处理设备 |
CN101390756A (zh) * | 2007-09-20 | 2009-03-25 | 深圳迈瑞生物医疗电子股份有限公司 | 超声成像系统中的扫描线插线处理方法 |
CN101441401A (zh) * | 2007-11-20 | 2009-05-27 | 深圳迈瑞生物医疗电子股份有限公司 | 快速确定超声成像系统中成像区域的方法和装置 |
CN101449983A (zh) * | 2007-11-29 | 2009-06-10 | 深圳迈瑞生物医疗电子股份有限公司 | 超声线阵偏转成像的扫描变换插值方法和装置 |
CN101729893A (zh) * | 2008-08-15 | 2010-06-09 | 北京北大众志微系统科技有限责任公司 | 基于软硬件协同处理的mpeg多格式兼容解码方法及其装置 |
CN102135615A (zh) * | 2010-03-22 | 2011-07-27 | 上海通号轨道交通工程技术研究中心有限公司 | 一种测距定标方法 |
CN102025960A (zh) * | 2010-12-07 | 2011-04-20 | 浙江大学 | 一种基于自适应插值的运动补偿去隔行方法 |
CN102222366A (zh) * | 2011-07-29 | 2011-10-19 | 电子科技大学 | 复杂空间曲面的拟合方法 |
CN103837608A (zh) * | 2014-03-12 | 2014-06-04 | 深圳市神视检验有限公司 | 一种相控阵接收动态聚焦补偿方法及系统 |
CN104359982A (zh) * | 2014-11-19 | 2015-02-18 | 哈尔滨工业大学 | 一种薄板焊缝超声相控阵检测楔块的方法 |
Non-Patent Citations (2)
Title |
---|
张红旭: "医学超声图像数字扫描变换方法研究及硬件实现", 《天津大学硕士学位论文》 * |
许药林: "超声相控阵成像算法研究及软件系统设计", 《南京航空航天大学硕士学位论文》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107255674A (zh) * | 2017-06-21 | 2017-10-17 | 景德镇陶瓷大学 | 一种圆柱壳内壁涂层界面缺陷的超声检测扫描成像装置及其应用方法 |
CN107356670A (zh) * | 2017-07-14 | 2017-11-17 | 武汉工程大学 | 一种基于斜入射的超声相控阵焊缝缺陷检测方法 |
CN108428251A (zh) * | 2018-03-09 | 2018-08-21 | 深圳市中捷视科科技有限公司 | 一种基于机器视觉技术激光结构光自动标定方法 |
CN110415316A (zh) * | 2019-07-02 | 2019-11-05 | 中国航空工业集团公司雷华电子技术研究所 | 一种提高气象视频显示准确性和连续性的插值方法 |
CN110415316B (zh) * | 2019-07-02 | 2023-05-05 | 中国航空工业集团公司雷华电子技术研究所 | 一种提高气象视频显示准确性和连续性的插值方法 |
CN111175381A (zh) * | 2019-12-09 | 2020-05-19 | 中北大学 | 基于全矩阵数据的复合构件界面快速成像定量检测方法 |
CN111175381B (zh) * | 2019-12-09 | 2022-12-06 | 中北大学 | 基于全矩阵数据的复合构件界面快速成像定量检测方法 |
CN112665528A (zh) * | 2020-12-31 | 2021-04-16 | 河南卫华重型机械股份有限公司 | 一种激光扫描三维成像的矫正方法 |
CN112665528B (zh) * | 2020-12-31 | 2023-10-20 | 河南卫华重型机械股份有限公司 | 一种激光扫描三维成像的矫正方法 |
CN113240587A (zh) * | 2021-07-12 | 2021-08-10 | 深圳华声医疗技术股份有限公司 | 超分辨率扫描变换方法、装置、超声设备及存储介质 |
CN113848375A (zh) * | 2021-10-20 | 2021-12-28 | 国网湖南省电力有限公司 | 一种油绝缘变压器内部器件带电检测装置及其应用方法 |
CN113848375B (zh) * | 2021-10-20 | 2023-11-28 | 国网湖南省电力有限公司 | 一种油绝缘变压器内部器件带电检测装置及其应用方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106175839A (zh) | 用于超声相控阵的扇形扫描成像的扫描检测方法和装置 | |
EP2639581B1 (en) | Ultrasonic inspection method | |
CN101449983B (zh) | 超声线阵偏转成像的扫描变换插值方法和装置 | |
EP2356941B1 (en) | Ultrasonographic system | |
US5172343A (en) | Aberration correction using beam data from a phased array ultrasonic scanner | |
US6923066B2 (en) | Ultrasonic transmitting and receiving apparatus | |
US5415173A (en) | Ultrasound diagnosis system | |
CN109490419A (zh) | 一种全聚焦成像的声束校准方法 | |
GB1579939A (en) | Apparatus for examining at least one object by acoustic wave energy | |
Ranganathan et al. | A novel beamformer design method for medical ultrasound. Part I: Theory | |
JPH05228147A (ja) | 流動する反射体の大きさ及び方向の表示を求める方法及びコヒーレント作像システム | |
JPH1156837A (ja) | 超音波イメージング・システム及び超音波イメージング・システムの開口をアポダイズする方法 | |
US10197670B2 (en) | Ultrasonic measurement apparatus, ultrasonic diagnostic apparatus, and ultrasonic measurement method | |
KR0180057B1 (ko) | 초음파시스템의 3차원 영상 획득장치 | |
CN108226304A (zh) | 一种基于测量模型的超声相控阵线扫描灵敏度计算方法 | |
JP5593208B2 (ja) | 超音波探傷装置および超音波探傷方法 | |
US20020193688A1 (en) | Three-dimensional ultrasound imaging system for performing receive-focusing at voxels corresponding to display pixels | |
US20170032557A1 (en) | Ultrasound Focal Zone System and Method | |
US10042052B2 (en) | Method and system for determining a location of a reflecting scatterer in a medium | |
Kim et al. | High volume rate 3-D ultrasound imaging using cross array based on synthetic transmit focusing | |
JP2013102959A (ja) | 超音波診断装置及び方法 | |
JP6047041B2 (ja) | 超音波診断装置、超音波診断装置の信号処理方法およびプログラム | |
CN107340331A (zh) | 用于板状结构检测的真时延无频散sh0波相控阵系统 | |
Cosarinsky et al. | Fast and automatic array tilt compensation for 1.5 D array transducers | |
KR100543736B1 (ko) | 초음파 영상 표시 방법 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20161207 |