CN101559001A - 一种超声扫查方法和设备 - Google Patents
一种超声扫查方法和设备 Download PDFInfo
- Publication number
- CN101559001A CN101559001A CNA2009100816131A CN200910081613A CN101559001A CN 101559001 A CN101559001 A CN 101559001A CN A2009100816131 A CNA2009100816131 A CN A2009100816131A CN 200910081613 A CN200910081613 A CN 200910081613A CN 101559001 A CN101559001 A CN 101559001A
- Authority
- CN
- China
- Prior art keywords
- scanning
- image data
- scanning probe
- scan image
- data
- 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
Images
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
一种超声扫查方法和设备,由扫查探头进行扫描获得两组三维超声图像数据;对两组三维超声图像数据进行融合,重构出具有良好效果的三维超声图像数据。
Description
技术领域
本发明涉及一种超声扫查方法和设备,尤其是将同一被测对象的重合的扫查区域两组三维图像数据进行融合以获得一组高分辨率的三维超声图像数据的方法和设备。
背景技术
目前的医用超声波探头晶片基本上是一维矩阵,应用于超声诊断的探头一般是一维的,由于采用了动态聚焦技术,普通的一维探头在扫描线平面内(轴向和侧向)聚焦良好,然而在垂直于扫描平面(横向)上,声束聚集差,由于探头的晶阵本身存在宽度,导致在该平面内重构的图像分辨率较低,这在超声三维图像重构中影响整体图像质量。
目前国内外主要采取以下两种方法来解决这一问题:其一,将探头晶片表面横向制成弧形,或把探头表面的声透镜的横向上做出弧度,使得在横向上也有一定的机械聚焦。该方法除了工艺复杂外,只能作到横向上的固定焦点,对焦点外的其它深度,尤其是距该深度较远处得到的聚焦效果并不理想,因此,该方法不能从根本上提高三维图像特别是横向的分辨率;其二,采用二维晶片的探头(或所谓1.5维探头),在侧向和横向上都采用电子聚焦技术,可以在横向上实现良好聚焦。该方法从原理上可以解决横向分辨率低的问题,但是此种探头造价太高,同时,超声系统的发射/接收通道的造价也增加很多,不利于实际应用中的推广。
发明内容
为了克服背景技术存在的缺陷,本发明提供一种超声扫查方法和设备。
本发明的第一目的在于提供一种超声扫查方法,具体方法包括:
扫查探头沿第一方向移动进行扫查获得第一组扫查图像数据;
以Z轴为转轴,对扫查探头进行旋转90度,沿垂直于第一方向的第二方向移动进行扫查获得第二组扫查图像数据;
对两组扫查图像数据的影像重叠部分,根据扫查图像数据相关性,对两组扫查图像数据进行融合,获得第三组扫查图像数据。
其中所述的根据图像数据相关性,对两组扫查图像数据进行合成的基本原理是:
扫查探头Z轴垂直于被测对象的体表,沿扫查面移动,对被测对象进行扫查。
两组扫查图像数据所覆盖的为一被测对象的重合的扫查区域。
扫查探头沿第一方向进行扫查,其特征在于Z轴垂直于被测对象的体表,横向沿X轴,侧向沿Y轴放置,扫查探头沿X轴方向平行移动,对被测对象进行扫查。
扫查探头沿第二方向进行扫查,其特征在于Z轴垂直于被检测体的表面,横向沿Y轴,侧向沿X轴放置,扫查探头沿Y轴方向平行移动,对被测对象进行扫查。
扫查探头每间隔一定距离进行一次超声扫描,获得一帧超声扫描图像数据,S={X,Y,(x,y),s},其中包括获取该数据时的扫查探头X轴和Y轴摆放方向,该数据的位置坐标(x,y),以及该位置的一帧扫查图像数据s,该扫查图像数据为一矩阵,矩阵为p*q,p为扫查探头宽度方向的采样点个数,q为扫查的深度方向的采样点个数。
对被测对象的一扫查区域,扫查探头沿X轴方向平行移动,获得第一组扫查图像数据,CX={Sx1,Sx2,...Sxn},其中该组数据的y坐标相同,n为扫查探头沿X轴方向平行移动采样的个数。
对被测对象的该扫查区域,扫查探头沿Y轴方向平行移动,获得第二组扫查图像数据,CY={Sy1,Sy2,...Sym},其中该组数据的x坐标相同,m为扫查探头沿Y轴方向平行移动采样的个数。
对获得的第一组扫查图像数据,CX={Sx1,Sx2,...Sxn},进行三维合成转换,提取相同深度的数据,拼接成为平面,转换成为具有q个平面的扫查图像数据,TX={tx1,tx2,...txq},其中tx1,tx2,...txq为p*n的矩阵;
对获得的第二组扫查数据,CY={Sy1,Sy2,...Sym},进行三维合成转换,提取相同深度的数据,拼接成为平面,转换成为具有q个平面的扫查图像数据,TY={ty1,ty2,...tyq},其中ty1,ty2,...tyq为p*m的矩阵
对获得的三维图像数据TX、TY进行插值,插入值与最近邻的gij相等,插入个数与α相关,使tx1,tx2,...txq的任一矩阵
设置扫查探头横向步进的距离为αs,α为大于1的整数,s为扫查探头侧向的扫描分辨率。α的取值范围为,α≤(扫查探头横向最高分辨率÷扫查探头侧向分辨率)×2,α≥(扫查探头横向最高分辨率÷扫查探头侧向分辨率)÷2。
对获得的三维图像数据TX、TY进行插值,插入值与最近邻的gij相等,插入个数与α相关,使tx1,tx2,...txq的任一矩阵
成为
使txv′的任意相邻两点的X轴向和Y轴向物理间距一致,得到TX′;
使ty1,ty2,...tyq的任一矩阵,
成为
使tyv′的任意相邻两点的X轴向和Y轴向物理间距一致,得到TY′。
对所获得的三维图像数据TX′和TY′,对两组数据对应影像物理位置重叠的部分,对其每个相同深度的平面的数据{txv′,tyv′}分别进行,切分为多个正方矩阵Mu,根据图形相关性对每个正方矩阵Mu进行合成计算,其中Mu为:
其中,gij,i,j=1,2...α,为图像灰度矩阵中每个像素点的灰度值。
其中根据图形相关性对每个正方矩阵Mu进行合成计算,其特征在于包括:
最接近真实图像的图像灰度矩阵Mz,
txv′中对应Mz的扫查测量获得的平面真实灰度值矩阵Mtx为:
其中,gim i=1,2,3...α,表示图像第i行灰度的均值;;
tyv′中对应Mz的扫查测量获得的平面真实灰度值矩阵Mty为:
其中,gmj i=1,2,3...α,为图像第i列灰度的均值;
它们的均值关系为:
gi1+gi2+gi3+...+giα=α·gim
超声图像的相关性为:
gi1∶gi2∶gi3∶...∶giα=gm1∶gm2∶gm3∶...∶gmα
上式等同于:
将上式代入gi1+gi2+gi3+...+giα=α·gim得出:
k·gm1+k·gm2+k·gm3+...+k·gmα=α·gim
求得K值。
将K值代入:
即可求得gi1,gi2,gi3…giα的值。求出矩阵Mz的每个元素的灰度值,获得Mz矩阵的融合。
扫查探头对被测对象移动扫查包括至少一次的若干次平行扫描,相邻扫描区域有少量重叠,根据扫查数据的精确的物理位置,将多次扫描小区域所成图像拼合成一块无重叠的连续的大区域,即整个平面的图像,得到两幅无重叠的连续的三维图像。
本发明的第二目的是提供一种扫查设备,扫查设备具有一扫查平台,另有一中央处理器和一扫查数据存储器。
扫查平台内置有一扫查探头,扫查探头Z轴垂直于被测对象体表,沿第一方向和垂直于第一方向的第二方向移动进行扫查可以获得两组扫查图像数据,并将两组图像数据存储于存储器中。扫查探头沿第一方向,X轴方向进行扫查,横向沿X轴,侧向沿Y轴放置,平行移动对被测对象扫查;扫查探头沿第二方向,Y轴方向进行扫查,横向沿Y轴,侧向沿X轴放置,平行移动对被测对象扫查。
扫查平台内具有一旋转装置,旋转装置与扫查探头相连接,扫查探头固定于旋转装置前端。旋转装置以Z轴为转轴,驱动扫查探头进行旋转,其旋转范围包括正负90度,能控制扫查探头精确定位于正负90度的物理位置。
扫查平台具有一精确扫描位置定位模块,获取扫查探头扫查时的精确物理位置。位置定位模块驱动扫查探头沿X轴方向和Y轴方向平行移动,控制扫查探头位移至精确位置,在每间隔一定距离进行一次超声扫描,获得一帧扫查图像数据,S={X,Y,(x,y),s}。扫查探头横向每间隔一定的距离为αs,α为大于1的整数,s为扫查探头侧向的扫描分辨率;α的取值范围为,α≤(扫查探头横向最高分辨率÷扫查探头侧向分辨率)×2,α≥(扫查探头横向最高分辨率÷扫查探头侧向分辨率)÷2。
扫查设备中另有一中央处理器,由中央处理器控制,系统对获得的两组扫查图像数据,进行三维合成转换,提取相同的数据拼接成为一平面,转换成为具有多个平面的扫查图像数据。并根据两组扫查图像数据的相关性,对两组扫查图像数据进行合成计算,获得第三组扫查图像数据,并存储于存储器中。
本发明的第三目的是提供一种独立的扫查平台组件,扫查平台内置有一扫查探头,扫查探头具有扫查数据输出接口;扫查平台具有一旋转装置,旋转装置与扫查探头相连,旋转装置具有物理输入接口,接受控制命令;扫查平台具有一精确扫描位置定位模块,获取扫查探头扫查时的精确物理位置,并具有物理接口输出位置信息;扫查平台具有平行位移驱动模块,驱动扫查探头沿第一方向以及第二方向移动进行扫查,位移驱动模块具有物理输入接口,接受控制命令。这样,扫查平台可以通过上述接口和一计算机系统连接,其连接接口可以是常用的RS232,USB等接口,结合本发明的方法,加上设备扫查架等,组合成为一台扫查设备。
附图说明
图1是本发明设备的结构示意图
图2是图1沿C-C线的剖面视图
图3是扫查探头横向、侧向及轴向的示意图。
图4A是旋转扫查探头的驱动装置的示意图。
图4B是旋转扫查探头的驱动装置的流程图。
图5是扫查探头沿X轴方向和Y轴方向扫查的示意图。
图6A是扫查探头沿X轴方向移动获得第一组图像数据的示意图。
图6B是图6A获得的图像数据的分辨率的示意图。
图7A是扫查探头沿Y轴方向移动获得第二组图像数据的示意图。
图7B是图7A获得的图像数据的分辨率的示意图。
图8是将两组三维图像数据进行合成,重构出第三组图像数据的流程图。
图9是扫查探头沿X轴方向和Y轴方向进行一次扫查的示意图。
图10是扫查探头沿X轴方向进行至少一次的若干次扫描的示意图。
图11是扫查探头沿Y轴方向进行至少一次的若干次扫描的示意图。
图12是图10和图11所示,沿X轴和Y轴进行若干次扫描获得的图像数据,进行切分,重叠部分进行合成的流程图。
具体实施方式
以下结合附图对本发明的具体实施方式作进一步描述。请参阅图1和图2,其中图1是本发明设备的结构示意图,图2是图1沿C-C线的剖面视图。
本发明设备具有一扫查平台7,为一具有边框的矩形框架,采用坚固的塑料材质制成。该扫查平台7与被测对象紧密接触的一平面为一扫查面6,在进行扫查时,该扫查面6覆盖在待扫查的被测对象体表。扫查平台7框架的两侧边各装有一把手8供操作员扫查基本定位,从垂直方向压紧固定被测对象。扫查平台7的上部固定连接在可调节位置的摇臂架9的一端,该摇臂架9连接至超声扫查控制操作台10,在超声扫查控制操作台10上有一显示器13,该超声扫查控制操作台10内有一中央处理器(CPU)11和一存储器12。
扫查平台7内置有一扫查探头5,该扫查探头5可采用电子线阵探头,以满足不同的扫描范围要求。扫查平台7内还具有一旋转装置4,扫查探头5与旋转装置4相连接,固定于旋转装置4前端。扫查探头5Z轴垂直于被测对象体表,横向沿X轴,侧向沿Y轴放置,沿X轴方向移动,对被测对象进行扫查。以Z轴为转轴,旋转装置4驱动扫查探头5,进行旋转90度,扫查探头5Z轴垂直于被测对象体表,横向沿Y轴,侧向沿X轴放置,沿Y轴方向移动,对被测对象进行扫查。并将获得的两组扫查图像数据存储于存储器12中。
扫查平台7内还具有一精确扫描位置定位模块,用来获取扫查探头5扫查时的精确物理位置。一精确扫描位置定位模块为伺服电动机3,伺服电动机3与旋转装置4相连接。在中央处理器(CPU)11的控制下,伺服电动机3驱动扫查探头5沿X轴方向移动,在每间隔一定距离时进行一次超声扫描,获得一帧图像数据S={X,Y,(x,y),s},并获得第一组图像数据CX={Sx1,Sx2,...Sxn}。然后,旋转装置驱动扫查探头旋转90度,扫查探头沿Y轴方向移动,在每间隔一定距离时进行一次超声扫描,获得一帧图像数据S={X,Y,(x,y),s}。并获得同一被测对象的第二组图像数据CY={Sy1,Sy2,...Sym}。中央处理器(CPU)11根据两组图像数据相关性,对两组图像数据进行合成计算,获得第三组图像数据,并存储于存储器12中。
图3是扫查探头横向、侧向和轴向定义的示意图。
本发明采用的扫查探头5为电子线阵探头,包括横向14、侧向15和轴向16三个方向。扫查探头5对被测对象进行扫查,扫查探头5发送超声波在横向14上显示的平面为横向平面,在轴向16上显示的平面为成像平面。其中A为扫查探头侧向15的扫描超声束宽度,B为扫查探头横向14的扫描超声束宽度,C为扫查探头对被测对象进行扫查,形成的超声波在轴向16上聚集形状。
图4A是旋转扫查探头的驱动装置的示意图。
图4B旋转扫查探头的驱动装置的流程图。
以下结合图4A和图4B对旋转扫查探头的驱动装置的具体实施进行描述。扫查探头5与旋转装置4相连接,固定于旋转装置4前端,旋转装置4与伺服电动机3相连接,伺服电动机3与X轴运动导轨1和Y轴运动导轨2连接。
步骤401,中央处理器(CPU)通过导线与伺服电动机、旋转装置和扫查探头相连接。
步骤402,在中央处理器(CPU)控制下,伺服电动机3驱动扫查探头5,沿X轴运动导轨1,以步进间隔匀速沿X轴方向平行移动,对被测对象扫查。
步骤403,此时,以Z轴为转轴,旋转装置4驱动扫查探头5旋转正负90度,并控制扫查探头5精确定位于90度的物理位置。
步骤404,扫查探头5在伺服电动机3驱动下,沿Y轴运动导轨2,以步进间隔匀速沿Y轴方向平行移动,对被测对象扫查。
图5是扫查探头沿X轴方向和Y轴方向扫查的示意图。
在中央处理器(CPU)11控制下,将扫查平台7内的扫查探头5调整至扫查平台7的右上角,扫查探头5的起始位置,在伺服电动机3驱动下,沿X轴方向平行移动,对被测对象进行扫查,获取第一组图像数据,如图中(1)所示。以Z轴为转轴,旋转装置4驱动扫查探头5旋转90度,如图中(2)所示。旋转后的扫查探头5,以旋转后的第一位置为起始位置,在伺服电动机3的驱动下,沿Y轴方向平行移动,对同一被测对象的重合的扫查区域进行扫查,获取第二组图像数据,如图中(3)所示。
下面将结合示意图,以p=256,q=400,n=64,m=64,α=4,为一实施例说明本发明的方法和设备。
图6A是扫查探头沿X轴方向移动获得第一组图像数据的示意图。
将扫查探头5置于扫查平台7框架内的右上角,扫查探头5与扫查平台7框架内的右上角预留有一定的空间。在扫查平台7覆盖范围内,扫查探头5Z轴垂直于被测对象体表,横向沿X轴,侧向沿Y轴放置于被测对象的扫查区域内。
在中央处理器(CPU)11控制下,扫查探头5在伺服电动机3驱动下,沿X轴运动导轨1,以步进间隔均速沿X轴方向平行移动,并由伺服电动机控制扫查探头5位移至精确的物理位置,使扫查探头5在每间隔一定距离进行一次超声扫描,获得一帧图像数据,S={X,Y,(x,y),s},其中包括获取该数据时的扫查探头5X轴Y轴摆放方向,该数据的位置坐标(x,y),以及该位置的一帧图像数据s,该图像数据为一矩阵,矩阵p*q为256*400,p=256为扫查探头5宽度方向的采样点个数,q=400为扫查的深度方向的采样点个数。
扫查探头5对被测对象移动扫查,获得第一组系列二维图像数据,CX={Sx1,Sx2,...Sx64},其中该组数据的y坐标相同,n=64为扫查探头5沿X轴方向平行移动采样的个数。
系统在中央处理器(CPU)11控制下,将获得的第一组系列二维图像数据CX={Sx1,Sx2,...Sx64}进行三维合成转换,提取每个二维图像数据相同深度的数据,拼接成一平面,再提取每个二维图像数据相同深度的数据,再拼接成一平面,以此类推,将所有的第一组二维图像数据的任何深度的数据,都拼接成平面,转换成为具有q=400个平面的三维图像数据,如三维图像数据平面a,TX={tx1,tx2,...tx400},其中tx1,tx2,...tx400为p*n,256*64的矩阵。并将获得的第一组三维图像数据储存于存储器12中。
扫查探头5沿X轴方向移动扫查,获得第一组三维图像数据TX,TX为一矩阵。该矩阵行数为p=256,列数为n=64,由于扫查探头5侧向行间距与列间距不等,此矩阵不是一个正方矩阵,故要对此矩阵X轴方向,即行方向,任意相邻间隔两点间进行插值。设置扫查探头5横向步进距离为4s,α=4为列间距是行间距的倍数,s为扫查探头5侧向的扫描分辨率。如果扫查探头5步进距离大于该值4s,则合成图像的分辨率提高有限;如果扫查探头5步进距离小于该值4s,则运算量加大。因此,最优步进距离为接近横向最高分辨率/侧向分辨率,取为侧向相邻像素点距离的整数倍即可。
对三维图像数据TX进行插值,插入值与最近邻的gij相等,插入个数与α=4相关,使tx1,tx2,...tx400的任一矩阵,
成为
使txv′的任意相邻两点的X轴向和Y轴向物理间距一致,得到TX′。
上述所述的算法,并非仅限于所列举,凡是能够采用与本发明所述的类似或优化的算法都应该在本发明的保护范围之内。
图6B是图6A获得的图像数据的分辨率的示意图。
本发明所采用的扫查探头5晶片为一维线阵,扫查探头5沿X轴方向对被测对象进行扫查,在Y轴方向,由于是电子聚焦,所获得的图像数据分辨率高,而在X轴方向,由于是机械聚焦,分辨率低。因此图6A所获得的第一组三维图像数据在Y轴方向,分辨率高;而在X轴方向,分辨率低,只具有一维高分辨率图像数据。
图7A是扫查探头沿Y轴方向移动获得第二组图像数据的示意图。
在中央处理器(CPU)11控制下,扫查探头5在伺服电动机驱动下沿扫查面完成X轴方向平行移动扫查,获得第一组三维图像数据结束后,以Z轴为转轴,旋转装置4驱动扫查探头5进行旋转90度,旋转范围包括正负90度,并控制扫查探头5精确定位于90度的物理位置。
扫查探头5Z轴垂直于被测对象体表,横向沿Y轴,侧向沿X轴放置,在伺服电动机3驱动下,沿Y轴运动导轨2,以步进间隔匀速沿Y轴方向平行移动,对同一被测对象的重合的扫查区域进行扫查,并由伺服电动机3控制扫查探头5位移至精确的物理位置,使扫查探头5在每间隔一定距离进行一次超声扫描,获得一帧图像数据,S={X,Y,(x,y),s},其中包括获取该数据时的扫查探头5X轴Y轴摆放方向,该数据的位置坐标(x,y),以及该位置的一帧图像数据s,该图像数据为一矩阵,矩阵为p*q,256*400,p=256为扫查探头5宽度方向的采样点个数,q=400为扫查的深度方向的采样点个数。
扫查探头5对被测对象移动扫查,获得第二组系列二维图像数据,CY={Sy1,Sy2,...Sy64},其中该组数据的x坐标相同,m=64为扫查探头5沿Y轴方向平行移动采样的个数。
系统在中央处理器(CPU)11控制下,将获得的第二组系列二维图像数据CY={Sy1,Sy2,...Sy64},进行三维合成转换,提取每个二维图像数据相同深度的数据,拼接成一平面,再提取每个二维像数据相同深度的数据,再拼接成一平面,以此类推,将所有的第二组二维图像数据的任何深度的数据,都拼接成平面,转换成为具有q=400个平面的三维图像数据,如三维图像平面b,TY={ty1,ty2,...ty400},其中ty1,ty2,...ty400为p*m,256*64的矩阵。并将获得的第二组三维图像数据储存于存储器12中。
扫查探头5沿Y轴方向移动扫查,获得第二组三维图像数据TY,TY为一矩阵。该矩阵列数为p=256,行数为m=64,由于扫查探头侧向列间距与行间距不等,此矩阵不是一个正方矩阵,故要对此矩阵Y轴方向,即列方向,任意间隔的相邻两点间进行插值。设置扫查探头横向步进距离为4s,α=4为列间距是行间距的倍数,s为扫查探头侧向的扫描分辨率。如果扫查探头5步进距离大于该值4s,则合成图像的分辨率提高有限;如果扫查探头5步进距离小于该值4s,则运算量加大。因此,最优步进距离为接近横向最高分辨率/侧向分辨率,取为侧向相邻像素点距离的整数倍即可。
对三维图像数据TY进行插值,插入值与最近邻的gij相等,插入个数与α=4相关,使ty1,ty2,...ty400的任一矩阵,
成为
使tyv′的任意相邻两点的X轴向和Y轴向物理间距一致,得到TY′。
上述所述的算法,并非仅限于所列举,凡是能够采用与本发明所述的类似或优化的算法都应该在本发明的保护范围之内。
图7B是图7A获得的图像数据的分辨率的示意图。
本发明所采用的扫查探头5晶片为一维线阵,扫查探头5沿Y轴方向对被测对象进行扫查,在X轴方向,由于是电子聚焦,所获得的图像数据分辨率高,而在Y轴方向,由于是机械聚焦,分辨率低。因此图7A所获得的第二组三维图像数据在X轴方向,分辨率高;而在Y轴方向,分辨率低,只具有一维高分辨率图像数据。
图8是将两组三维图像数据进行合成,重构出第三组图像数据的流程图。
步骤801,通过对TX和TY进行插值获得TX′、TY′两组三维图像数据。
步骤802,根据两组图像数据的物理位置重叠部分,对两组每个相同深度的平面的数据{txv′,tyv′}分别进行,切分为多个正方矩阵Mu。任取它们某一相同深度的平面,设为ah、bh,根据图形不同方向的分辨率不同对每个正方矩阵Mu进行合成计算,其中Mu为:
其中,gij,i,j=1,2,3,4,为图像灰度矩阵中每个像素点的灰度值。
步骤803,平面ah在Y轴方向由于是电子聚焦,分辨率高,而在X轴方向是机械聚焦,分辨率低,假设散焦范围为4个像素,即可以近似用这4点的平均值代替真实值表示txv′中对应Mz的扫查测量获得的平面真实灰度值矩阵Mtx。Mz为最接近真实图像的图像灰度矩阵。
矩阵Mtx为:
其中,gim i=1,2,3,4,表示图像第i行灰度的均值;
步骤804,平面bh在X轴方向由于是电子聚集,分辨率高,而在Y轴方向是机械聚焦,分辨率低。假设散焦范围为4个像素,即可以近似用这4点的平均值代替真实值表示tyv′中对应Mz的扫查测量获得的平面真实灰度值矩阵Mty。Mz为最接近真实图像的图像灰度矩阵。
矩阵Mty为:
其中,gmi i=1,2,3,4,为图像第i列灰度的均值;
步骤805,可见矩阵Mtx与矩阵Mty各具有一维高分辨率分量,矩阵Mtx的高分辨率分量为(g164 g264 g364 g464)T,矩阵Mty的高分辨率分量为(g641 g642 g643 g644)。
它们的均值关系为:
gi1+gi2+gi3+gi4=4·gi64
超声图像的相关性为:
gi1∶gi2∶gi3∶gi4=g641∶g642∶g643∶g644
上式等同于:
将上式代入gi1+gi2+gi3+gi4=4·gi64得出:
k·g641+k·g642+k·g643+k·g644=4·gi64
求得K值。
将K值代入:
即可求得gi1,gi2,gi3,gi4的值。求出矩阵Mz的每个元素的灰度值,获得Mz矩阵的融合。
上述所述的算法,并非仅限于所列举,凡是能够采用与本发明所述的类似或优化的算法都应该在本发明的保护范围之内。
图9是扫查探头沿X轴和Y轴方向进行一次扫查的示意图。
在中央处理器(CPU)11控制下,扫查探头5横向沿X轴,侧向沿Y轴放置,Z轴垂直于被测对象体表,在伺服电动机3驱动下,沿X轴方向平行移动,对被测对象扫查。扫查范围覆盖整个被测对象,获取第一组三维图像数据a。
扫查探头5在完成一次沿X轴方向扫查后,以Z轴为转轴,在旋转装置4驱动下,旋转正90度,如(1)所示,并在其控制下,精确定位于正90度的物理位置上;或在旋转装置4驱动下,旋转负90度,如(2)所示,并在其控制下精确定位于负90度的物理位置上。此时,扫查探头横向沿Y轴,侧向沿X轴放置,Z轴垂直于被测对象体表。然后,扫查探头5在伺服电动机3驱动下,沿Y轴方向平行移动,对同一被测对象的重合的扫查区域进行扫查。扫查范围覆盖整个被测对象,获取第二组三维图像数据b。并将两组三维图像数据储存至存储器中。
图10是扫查探头沿X轴方向进行至少一次的若干次扫描的示意图。
扫查探头5沿X轴方向扫查被测对象,如果扫查探头5一次扫描行程不能覆盖整个被测对象体表,则可以分为若干次沿X轴方向平行扫描。扫查探头5在完成一次沿X轴方向扫描行程,获得一个小扫查区域的三维图像数据a1后,扫查探头5放置方向不变,在伺服电机3驱动下,做一次垂直于扫描行程间的切换运动。然后,扫查探头5在伺服电动机3驱下,再次沿X轴方向对被测对象另一个小扫查区域进行扫描,获得另一个三维图像数据a2。以此类推,完成下一个小扫查区域的扫描,获得三维图像数据a3,直到完成对整个被测对象的扫查。将获得的每个三维图像数据存储至储存器中。
由于相邻的扫描区域有部分重叠,根据扫查数据的精确的物理位置,任取a1、a2和a3图像数据某一相同深度的平面,在重叠的位置进行拼接,拼接成一被测对象完整的图像数据。同理,对a1、a2和a3图像数据每个深度的平面都在求出的重叠位置进行拼接。拼合成一块无重叠的连续的扫查区域,即整个平面a的三维图像数据,并将三维图像数据储存于存储器12中。
图11是扫查探头沿Y轴方向进行至少一次的若干次扫描的示意图。
扫查探头5沿X轴方向扫查被测对象,如果扫查探头5一次扫描行程不能覆盖整个被测对象体表,则可以分为若干次沿Y轴方向平行扫描。扫查探头5在完成一次沿Y轴方向扫描行程,获得一个小扫查区域的三维图像数据b1后,扫查探头放置方向不变,在伺服电机3驱动下,做一次垂直于扫描行程间的切换运动。然后,扫查探头5在伺服电动机驱下,再次沿Y轴方向对被测对象另一个小扫查区域进行扫描,获得另一个三维图像数据b2。以此类推,完成下一个小扫查区域的扫描,获得三维图像数据b3,直到完成对整个被测对象的扫查。将获得的每个三维图像存储至储存器中。
由于相邻的扫描区域有部分重叠,根据扫查的精确的物理位置,任取b1、b2和b3图像数据某一相同深度的平面,在重叠的位置进行拼接,拼接成一被测对象完整的图像数据。同理,对b1、b2和b3图像数据每个深度的平面都在重叠位置进行拼接。拼合成一块无重叠的连续的扫查区域,即整个平面b的三维图像数据,并将三维图像数据储存于存储器12中。
图12是图10和图11所示,沿X轴和Y轴进行若干次扫描获得的图像,进行切分,重叠部分进行合成的流程图。
步骤1201,扫查探头5沿X轴方向移动扫查,获得的第一组三维图像数据,即平面a,为一矩阵A和扫查探头5沿Y轴方向移动扫查,获得的第二组三维图像数据,即平面b,为另一矩阵B。
步骤1202,由于矩阵A和矩阵B行间距与列间距不等,不是一个正方形矩阵,因此对矩阵A的行方向和矩阵B列方向进行插值。使该两个矩阵的任意相邻两点的X轴方向和Y轴方向物理间距一致,得到矩阵A′和矩阵B′。
步骤1203,根据两组图像数据重叠的物理位置部分,对两组图像数据的每个相同深度的平面的数据分别进行切分为多个正方矩阵。任取它们任何一深度的平面,根据图形不同方向的分辨率不同,在重叠的位置进行拼接合成。
步骤1204,第一组图像数据,矩阵A,在其任何一深度的平面在Y轴方向由于是电子聚集,分辨率高,而在X轴方向是机械聚集,分辨率低。第二组图像数据,矩阵B,在其任何一深度的平面在X轴方向由于是电子聚集,分辨率高,而在Y轴方向是机械聚集,分辨率低。两组图像数据各具有一维高分辨率分量,将两组图像数据进行合成,即可重构出一组在X轴方向和Y轴方向上都具有高分辨率的三维图像数据。
虽然本发明是参考其优选实施例示出和描述的,但本领域的普通技术人员应该理解在不脱离附属的权利要求书所限定的本发明的精神和范围的情况下,可以进行形式和细节的各种改变。
Claims (23)
1、一种超声扫查方法,扫查探头对被测对象进行超声扫查,其特征在于包括:
扫查探头沿第一方向移动进行扫查获得第一组扫查图像数据;
旋转扫查探头,沿第二方向移动进行扫查获得第二组扫查图像数据;
对两组扫查图像数据的影像重叠部分的两组扫查图像数据进行合成,获得第三组扫查图像数据。
2、如权利要求1所述的方法,其特征在于扫查探头沿第一方向移动进行扫查获得第一组扫查图像数据;
以Z轴为转轴,对扫查探头进行旋转90度,沿第二方向移动进行扫查获得第二组扫查图像数据;
对两组扫查图像数据的影像重叠部分,根据扫查图像数据相关性,对两组扫查图像数据进行合成,获得第三组扫查图像数据。
3、如权利要求1所述的方法,其特征在于扫查探头Z轴垂直于被检测体的表面,沿扫查面移动,对被测对象进行扫查。
4、如权利要求1所述的方法,其特征在于两组扫查图像数据所覆盖的为一被测对象的重合的扫查区域。
5、如权利要求1所述的方法,其中扫查探头沿第一方向进行扫查,其特征在于Z轴垂直于被检测体的表面,横向沿X轴,侧向沿Y轴放置,扫查探头沿X轴方向平行移动,对被测对象进行扫查。
6、如权利要求1所述的方法,其中扫查探头沿第二方向进行扫查,其特征在于Z轴垂直于被检测体的表面,横向沿Y轴,侧向沿X轴放置,扫查探头沿Y轴方向平行移动,对被测对象进行扫查。
7、如权利要求1所述的方法,其特征在于扫查探头每间隔一定距离进行一次超声扫描,获得一帧超声扫描图像数据,S={X,Y,(x,y),s},其中包括获取该数据时的扫查探头X轴和Y轴摆放方向,该数据的位置坐标(x,y),以及该位置的一帧扫查图像数据s,该扫查图像数据为一矩阵,矩阵为p*q,p为扫查探头宽度方向的采样点个数,q为扫查的深度方向的采样点个数。
8、如权利要求1所述的方法,其特征在于,对被测对象的一扫查区域,扫查探头沿X轴方向平行移动,获得第一组扫查图像数据,CX={Sx1,Sx2,...Sxn},其中该组数据的y坐标相同,n为扫查探头沿X轴方向平行移动采样的个数;
对被测对象的该扫查区域,扫查探头沿Y轴方向平行移动,获得第二组扫查图像数据,CY={Sy1,Sy2,...Sym},其中该组数据的x坐标相同,m为扫查探头沿Y轴方向平行移动采样的个数。
9、如权利要求1或8所述的方法,其特征在于,对获得的第一组扫查图像数据,CX={Sx1,Sx2,...Sxn},进行三维合成转换,提取相同深度的数据,拼接成为平面,转换成为具有q个平面的扫查图像数据,TX={tx1,tx2,...txq},其中tx1,tx2,...txq为p*n的矩阵;
对获得的第二组扫查数据,CY={Sy1,Sy2,...Sym},进行三维合成转换,提取相同深度的数据,拼接成为平面,转换成为具有q个平面的扫查图像数据,TY={ty1,ty2,...tyq},其中ty1,ty2,...tyq为p*m的矩阵。
10、如权利要求1或9所述的方法,其特征在于设置扫查探头横向步进的距离为αs,α为大于1的整数,s为扫查探头侧向的扫描分辨率。
11、如权利要求10所述,其特征在于α的取值范围为,α≤(扫查探头横向最高分辨率÷扫查探头侧向分辨率)×2,α≥(扫查探头横向最高分辨率÷扫查探头侧向分辨率)÷2。
12、如权利要求1、9或10所述的方法,其特征在于对获得的三维图像数据TX、TY进行插值,插入值与最近邻的gij相等,插入个数与α相关,使tx1,tx2,...txq的任一矩阵
txv=
成为
txv’=
使txv′的任意相邻两点的X轴向和Y轴向物理间距一致,得到TX′;
使ty1,ty2,...tyq的任一矩阵,
tyv=
成为
tyv’=
使tyv′的任意相邻两点的X轴向和Y轴向物理间距一致,得到TY′。
13、如权利要求1或12所述的方法,其特征在于,对所获得的三维图像数据TX′和TY′,对两组数据对应影像物理位置重叠的部分,对其每个相同深度的平面的数据{txv′,tyv′}分别进行,切分为多个正方矩阵Mu,根据图形相关性对每个正方矩阵Mu进行合成计算,其中Mu为:
其中,gij,i,j=1,2...α,为图像灰度矩阵中每个像素点的灰度值。
14、如权利要求1或13所述的方法,其中根据图像相关性对每个正方矩阵Mu进行合成计算,其特征在于包括:
最接近真实图像的图像灰度矩阵Mz,
txv′中对应Mz的扫查测量获得的平面真实灰度值矩阵Mtx为:
其中,gim i=1,2,3...α,表示图像第i行灰度的均值;
tyv′中对应Mz的扫查测量获得的平面真实灰度值矩阵Mty为:
其中,gmi i=1,2,3...α,为图像第i列灰度的均值;
它们的均值关系为:
gi1+gi2+gi3+...+giα=α·gim
超声图像的相关性为:
gi1∶gi2∶gi3∶...∶giα=gm1∶gm2∶gm3∶...∶gmα
上式等同于:
将上式代入gi1+gi2+gi3+...+giα=α·gim得出:
k·gm1+k·gm2+k·gm3+...+k·gmα=α·gim
求得K值。
将K值代入:
求得gi1,gi2,gi3…giα的值。求出矩阵Mz的每个元素的灰度值,获得Mz矩阵的融合。
15、如权利要求1或8所述的方法,其特征在于,扫查探头对被检测体移动扫查包括至少一次的若干次平行扫描,相邻扫描区域有少量重叠,根据扫查数据的精确的物理位置,将多次扫描小区域所成图像拼合成一块无重叠的连续的大区域,即整个平面的图像,得到两幅无重叠的连续的三维图像。
16、一种超声扫查设备,具有一扫查平台,另有一中央处理器和一扫查数据存储器,其特征在于包括:
扫查平台内置有一超声扫查探头;
扫查平台具有一旋转装置,旋转装置与扫查探头相连;
扫查平台具有一精确扫描位置定位模块,获取扫查探头扫查时的精确物理位置;
扫查平台驱动扫查探头沿第一方向移动进行扫查获得第一组扫查图像数据,存储于存储器;
扫查平台驱动旋转扫查探头,沿垂直于第一方向的第二方向移动进行扫查获得第二组扫查图像数据,存储于存储器;
对两组扫查图像数据的影像重叠部分,中央处理器根据扫查图像数据相关性,对两组扫查图像数据进行合成计算,获得第三组扫查图像数据,存储于存储器。
17、如权利要求16所述的设备,其特征在于扫查探头Z轴垂直于被检测体的表面,沿扫查面移动,对被测对象进行扫查。
18、如权利要求17所述的设备,其特征在于旋转装置连接扫查探头,旋转装置以Z轴为转轴,驱动扫查探头进行旋转,其旋转范围包括正负90度,能精确定位于正负90度。
19、如权利要求16所述的设备,其中扫查探头沿第一方向进行扫查,其特征在于Z轴垂直于被检测体的表面,横向沿X轴,侧向沿Y轴放置,扫查探头沿X轴方向平行移动,对被测对象进行扫查;
其中扫查探头沿第二方向进行扫查,其特征在于Z轴垂直于被检测体的表面,横向沿Y轴,侧向沿X轴放置,扫查探头沿Y轴方向平行移动,对被测对象进行扫查。
20、如权利要求16所述的设备,其特征在于扫查平台驱动扫查探头每间隔一定距离进行一次超声扫描,获得一帧扫描图像数据,S={X,Y,(x,y),s},其中包括获取该数据时的扫查探头X轴和Y轴摆放方向,该数据的位置坐标(x,y),以及该位置的一帧扫查图像数据s,该扫查图像数据为一矩阵,矩阵为p*q,p为扫查探头宽度方向的采样点个数,q为扫查的深度方向的采样点个数。
21、如权利要求20所述,其特征在于扫查探头横向每间隔一定的距离为αs,α为大于1的整数,s为扫查探头侧向的扫描分辨率;α的取值范围为,α≤(扫查探头横向最高分辨率÷扫查探头侧向分辨率)×2,α≥(扫查探头横向最高分辨率÷扫查探头侧向分辨率)÷2。
22、如权利要求16所述的设备,其特征在于,对被测对象的一扫查区域,扫查探头沿X轴方向平行移动,获得第一组扫查图像数据,CX={Sx1,Sx2,...Sxn},其中该组数据的y坐标相同,n为扫查探头沿X轴方向平行移动采样的个数;
对被测对象的该扫查区域,扫查探头沿Y轴方向平行移动,获得第二组扫查图像数据,CY={Sy1,Sy2,...Sym},其中该组数据的x坐标相同,m为扫查探头沿Y轴方向平行移动采样的个数;
中央处理器对获得的第一组扫查图像数据,CX={Sx1,Sx2,...Sxn},进行三维合成转换,提取相同深度的数据,拼接成为平面,转换成为具有q个平面的扫查图像数据,TX={tx1,tx2,...txq},其中tx1,tx2,...txq为p*n的矩阵;
中央处理器对获得的第二组扫查图像数据,CY={Sy1,Sy2,...Sym},进行三维合成转换,提取相同深度的数据,拼接成为平面,转换成为具有q个平面的扫查图像数据,TY={ty1,ty2,...tyq},其中ty1,ty2,...tyq为p*m的矩阵;
中央处理器根据扫查图像数据相关性,对两组扫查图像数据进行合成计算,获得第三组扫查图像数据,存储于存储器。
23、一种扫查平台,其特征在于包括:
扫查平台内置有一扫查探头,扫查探头具有扫查数据输出接口;
扫查平台具有一旋转装置,旋转装置与扫查探头相连,旋转装置具有物理输入接口,接受控制命令;
扫查平台具有一精确扫描位置定位模块,获取扫查探头扫查时的精确物理位置,并具有物理接口输出位置信息;
扫查平台具有平行位移驱动模块,驱动扫查探头沿第一方向以及第二方向移动进行扫查,位移驱动模块具有物理输入接口,接受控制命令。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNA2009100816131A CN101559001A (zh) | 2009-04-07 | 2009-04-07 | 一种超声扫查方法和设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNA2009100816131A CN101559001A (zh) | 2009-04-07 | 2009-04-07 | 一种超声扫查方法和设备 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101559001A true CN101559001A (zh) | 2009-10-21 |
Family
ID=41218098
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNA2009100816131A Pending CN101559001A (zh) | 2009-04-07 | 2009-04-07 | 一种超声扫查方法和设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101559001A (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102144931A (zh) * | 2011-04-10 | 2011-08-10 | 陈莉 | 一种类显微三维超声成像法 |
CN102546031A (zh) * | 2010-12-16 | 2012-07-04 | 李百祺 | 无线功率传输系统、无线功率传送装置与无线功率接收装置 |
CN102743188A (zh) * | 2011-04-22 | 2012-10-24 | 李百祺 | 超声波自动扫描系统及其扫描方法 |
WO2013021574A1 (en) * | 2011-08-08 | 2013-02-14 | Canon Kabushiki Kaisha | Object information acquisition apparatus, object information acquisition system, display control method, display method, and program |
CN104688268A (zh) * | 2015-02-11 | 2015-06-10 | 戴守平 | 超声检查辅助装置 |
CN105147398A (zh) * | 2015-09-30 | 2015-12-16 | 深圳市人民医院 | 一种用于筛查和监测乳腺癌的多模影像学系统 |
CN105184760A (zh) * | 2014-05-30 | 2015-12-23 | 财团法人金属工业研究发展中心 | 牙体影像的接合方法 |
CN106691507A (zh) * | 2017-01-04 | 2017-05-24 | 沈阳东软医疗系统有限公司 | 一种4d超声探头扫描控制方法、装置和系统 |
CN109069131A (zh) * | 2016-04-18 | 2018-12-21 | 皇家飞利浦有限公司 | 用于乳腺组织成像的超声系统和方法 |
CN110327077A (zh) * | 2019-07-09 | 2019-10-15 | 深圳开立生物医疗科技股份有限公司 | 一种血流显示方法、装置及超声设备和存储介质 |
CN113812978A (zh) * | 2021-10-25 | 2021-12-21 | 深圳市德力凯医疗设备股份有限公司 | 数据采样方法、病灶部位检查方法及智能终端 |
-
2009
- 2009-04-07 CN CNA2009100816131A patent/CN101559001A/zh active Pending
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102546031A (zh) * | 2010-12-16 | 2012-07-04 | 李百祺 | 无线功率传输系统、无线功率传送装置与无线功率接收装置 |
CN102144931B (zh) * | 2011-04-10 | 2013-01-30 | 陈莉 | 一种类显微三维超声成像法 |
CN102144931A (zh) * | 2011-04-10 | 2011-08-10 | 陈莉 | 一种类显微三维超声成像法 |
CN102743188A (zh) * | 2011-04-22 | 2012-10-24 | 李百祺 | 超声波自动扫描系统及其扫描方法 |
CN102743188B (zh) * | 2011-04-22 | 2016-01-13 | 李百祺 | 超声波自动扫描系统及其扫描方法 |
WO2013021574A1 (en) * | 2011-08-08 | 2013-02-14 | Canon Kabushiki Kaisha | Object information acquisition apparatus, object information acquisition system, display control method, display method, and program |
CN105184760B (zh) * | 2014-05-30 | 2018-12-04 | 财团法人金属工业研究发展中心 | 牙体影像的接合方法 |
CN105184760A (zh) * | 2014-05-30 | 2015-12-23 | 财团法人金属工业研究发展中心 | 牙体影像的接合方法 |
CN104688268A (zh) * | 2015-02-11 | 2015-06-10 | 戴守平 | 超声检查辅助装置 |
CN105147398A (zh) * | 2015-09-30 | 2015-12-16 | 深圳市人民医院 | 一种用于筛查和监测乳腺癌的多模影像学系统 |
CN109069131A (zh) * | 2016-04-18 | 2018-12-21 | 皇家飞利浦有限公司 | 用于乳腺组织成像的超声系统和方法 |
CN106691507A (zh) * | 2017-01-04 | 2017-05-24 | 沈阳东软医疗系统有限公司 | 一种4d超声探头扫描控制方法、装置和系统 |
CN110327077A (zh) * | 2019-07-09 | 2019-10-15 | 深圳开立生物医疗科技股份有限公司 | 一种血流显示方法、装置及超声设备和存储介质 |
CN113812978A (zh) * | 2021-10-25 | 2021-12-21 | 深圳市德力凯医疗设备股份有限公司 | 数据采样方法、病灶部位检查方法及智能终端 |
CN113812978B (zh) * | 2021-10-25 | 2023-08-18 | 深圳市德力凯医疗设备股份有限公司 | 数据采样方法、病灶部位检查方法及智能终端 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101559001A (zh) | 一种超声扫查方法和设备 | |
CN101315341B (zh) | 分层摄影检测系统与方法 | |
CN101569541B (zh) | 三维超声波成像系统 | |
CN100582657C (zh) | 三维微观形貌斜扫描方法及装置 | |
CN1264062C (zh) | 一种多视角x射线立体成像的方法与系统 | |
CN102692214A (zh) | 一种狭窄空间双目视觉测量定位装置及方法 | |
CN1155799C (zh) | 激光线扫描三维测量双光刀虚拟网格映射标定方法及装置 | |
KR20080097218A (ko) | 현미경 매체기반 시료로부터 디지털 이미지데이터를 수집하기 위한 방법 및 장치와 컴퓨터 프로그램 제품 | |
JP2005533320A (ja) | 画像解析の方法、装置およびソフトウェア | |
CN108253939A (zh) | 可变视轴单目立体视觉测量方法 | |
CN104019829A (zh) | 一种基于pos系统的车载全景相机和线阵激光扫描仪的外参标定方法 | |
CN100998512A (zh) | 三维超声影像的重建方法 | |
CN108318581A (zh) | 一种不用装夹定位的圆弧面工件超声波c扫描自动检测方法 | |
CN101526337A (zh) | 三维图像扫描系统及三维图像扫描方法 | |
CN111189415A (zh) | 一种基于线结构光的多功能三维测量重建系统及方法 | |
CN104881864A (zh) | 人体头部三维扫描仪及三维建模方法 | |
EP3471063A1 (en) | Three-dimensional imaging method and system | |
CN110617763A (zh) | 一种旋转结构光式的h型钢尺寸测量方法 | |
CN116100204A (zh) | 一种基于坐标映射算法的定位引导方法 | |
CN102692421B (zh) | 一种具有计量转轴的高精度x射线显微镜扫描样品台 | |
CN216449449U (zh) | 一种表面检测装置 | |
CN103860197A (zh) | 自由式三维超声成像方法、装置和系统 | |
CN108322736B (zh) | 一种多个线阵相机绕视轴旋转角度校准的校准板和校准方法 | |
CN114359393B (zh) | 跨平台的视觉引导点胶引导方法 | |
CN102692422A (zh) | 一种计量型高精度x射线显微镜扫描样品台 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
DD01 | Delivery of document by public notice |
Addressee: Beijing United Imaging Systems Co., Ltd. Document name: Notification of an Office Action |
|
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Open date: 20091021 |