具体实施方式
为了在偏转凸阵数字扫描变换时实现一次坐标变换,而减少插值引入的噪声、减少运算复杂度,本发明首先查找至少两个与待求象素点最接近的邻近象素点;然后,获取邻近象素点的超声回波数据,利用该超声回波数据通过插值运算,获得待求象素点的超声回波数据,从而构建扫描变换后的超声图像。
上述插值运算过程中,可以以一距离因子和/或一角度因子作为权重,对邻近象素点的超声回波数据进行插值计算,该角度因子由所述待求像素点相对于所述邻近像素点在偏转极坐标系上的归一化夹角决定,该距离因子由所述待求像素点相对于所述邻近像素点在偏转极坐标系上的归一化距离决定。
为了实现上述查找功能,并减少插值运算所带入的误差和运算复杂度,本发明的查找过程中,通过计算待求象素点在探头圆弧上所对应的射线起始点,并以该起始点在极坐标系下的角度、以及该起始点与待求像素点的距离为参考值查找与待求象素点最接近的邻近象素点。其具体过程如下:
首先,对笛卡儿坐标系中的任一待求象素点,计算该待求象素点在探头圆弧上所对应的射线起始点,在极坐标系下的角度。在这一步骤中,为计算所述角度,可以采用如下方式:
A11、根据探头圆弧圆心、所述待求像素点及所述射线起始点在极坐标系下的位置关系,构建关于所述角度的三角函数关系式;
A12、结合笛卡儿坐标系中探头圆弧圆心和待求象素点的坐标,计算所述三角函数关系式中的中间量,其中的中间量可以参照公式(1)。
A13、将所述中间量的计算结果代入到三角函数关系式中,获得所述角度的值。
关于所述角度的三角函数计算关系式,参见以下公式(1)为:
上式中,(x,y)为笛卡尔坐标系中待求象素点的坐标;(x0,y0)为探头圆弧圆心的坐标,ROB为凸阵圆弧半径,θ为凸阵偏转角度,θB为所述射线起始点在极坐标下的角度。
然后,计算笛卡儿坐标系中待求象素点到所述射线起始点的距离。这一步过程中,利用三角函数关系获得的待求象素点到所述射线起始点的距离如下式(2):
上式中,(x,y)为笛卡尔坐标系中待求象素点的坐标;(x0,y0)为探头圆弧圆心的坐标,θ为凸阵偏转角度,θB为所述射线起始点在极坐标下的角度,如公式(1)。
再次,根据所述角度和距离的计算结果,在笛卡儿坐标系中查找至少两个与待求象素点距离最接近的邻近象素点,并且根据邻近象素点的个数多少需要满足上述条件。以下详细说明本发明的如何利用所述角度和距离的计算结果查找邻近象素点的,首先,根据扫描线密度,将所述角度和所述距离分别分解成整数部分和小数部分,所述角度的整数部分对应于扫描线线号,所述距离的整数部分对应于沿着扫描线的扫描深度上的采样点;然后,根据所述角度和所述距离的整数部分,查找与待求象素点最相近的邻近象素点。而对于所述角度和所述距离的小数部分可作为权重,参加插值运算。
下面结合附图,以一个具体实施例为例详细说明本发明上述计算过程的原理。
如图3所示,本实施例的具体过程如下:
1、假设笛卡儿坐标系的一个像素点A(x,y),定义该点为上述待求象素点;
2、计算极坐标系下像素点A(x,y)对应圆弧上射线起始点B的角度θB。如图4所示,给出了极坐标系下圆弧射线起始点计算示意图。射线OA、OC为极坐标系下凸阵正常扫描情况下的探头圆弧圆心O点所对应的射线,即其对应凸阵探头圆弧的法线方向。线段AB表示偏转凸阵的实际扫描线,所以,极坐标系下像素点A(x,y)对应圆弧上射线起始点B也就是偏转凸阵实际扫描线的起始点。射线OB为凸阵正常扫描情况下的一根扫描线,即其对应凸阵探头圆弧的法线方向。其中,图4给出的线段AB与射线OB之间的夹角θ即偏转凸阵实际扫描线偏转正常扫描线的角度,也就是通过控制发射脉冲的延时使发射扫描线偏转的偏转角度,属于已知的凸阵偏转扫描角度。
3、计算像素点A距起始点B的距离RAB。
4、根据扫描线密度将RAB和θB分解成整数部分(RIAB,θIB)和小数部分(RFAB,θFB)。
5、根据整数部分(RIAB,θIB)可以查找到和点A最相近的点,本发明以四点插值为例,根据小数部分(RFAB,θFB)对这四个点进行插值计算。
6、最后经过插值计算后的超声波信息值用于凸阵空间复合处理。下面给出算法的具体流程:
计算极坐标下对应圆弧上射线起始点B的角度θB可通过简单的三角函数正弦定理获得。如图4,给出了起始点B的角度θB及像素点A距离起始点B的距离RAB的计算示意图。如图4所示,A点为所要进行扫描变换像素点,O(x0,y0)点为凸阵探头圆弧的圆心,B点为所求的起始点位置。
利用三角函数正弦定理,有公式(3):
其中ROB为凸阵圆弧半径,ROA为像素点A(x,y)距离凸阵探头圆弧圆心的直线距离,即计算辅助线,可在笛卡儿坐标系中利用A点坐标(x,y)和O点的坐标(xo,yo)计算出来,如公式(4):
∠OBA是凸阵偏转角度θ的补角,即:∠OBA+θ=π,因此有公式(5):
sin(∠OBA)=sin(θ) (5)
这样,可通过下式(6)计算出∠OAB:
随之,可通过下式(7)计算∠AOB:
∠AOB=θ-∠OAB (7)
∠AOC又可通过点A和O的坐标而确定,如下式(8)
这样,通过∠AOC和∠AOB可以确定圆弧上射线起始点B的角度θB:
θB=∠BOC=∠AOC-∠AOB (9)
依次将上述中间量代入上述公式(9)计算θB可得上述公式(1),即,
像素点A距离起始点B的距离RAB也可以通过相应的正弦定理来确定,如下式(10):
即:
将上述各个公式的结果带入公式(11)可得上述公式(2),即,
如上,通过计算得到了起始点B的角度θB及像素点A距起始点B的距离RAB,将RAB和θB根据扫描线密度分解成整数部分(RIAB,θIB)和小数部分(RFAB,θFB)。整数部分θIB对应于扫描线线号,整数部分RIAB对应于沿着扫描线的扫描深度上的采样点。扫描变换器中包含一个存储器,其本质是一块内存器,用于存储超声回波数据。其内存的存储格式一般采用二维存储的方法,对应于实际超声成像时系统的采样。其中,一个方向是沿着扫描线线号变换的方向,另外一个方向是沿着扫描线上深度的变化。根据(RIAB,θIB)可以查到存储器中相邻的四个位置点所对应的超声回波数据。根据这四个相邻位置点的超声回波数据,经过插值计算处理,得到(RAB和θB)位置也就是点A对应的超声回波数据。插值处理时的系数和小数部分(RFAB,θFB)相关。小数部分(RFAB,θFB)通过一定的变换进行是归一化。两个相邻点的深度距离归一化为1,两条相邻扫描线之间的交角归一化为1。
如图5所示,不妨假设取得的四个邻近象素点分别是P1,P2,Q1,Q2。其中P1,P2的扫描深度相同,不妨假设其扫描深度是m。Q1,Q2的扫描深度相同,其扫描深度是m+1。P1和Q1,P2和Q2的扫描线号相同,不妨假设其对应的扫描线号分别是n,n+1,这样选邻近像素点的好处在于用得是最近的点差值,误差最小。这四个邻近象素点对应的超声回波信息(主要包含的是像素点的幅值信息)分别为XP1,XP2,XQ1,XQ2。目的是为了计算位置A的超声回波数据,其中位置点A对应的小数部分是(RFAB,θFB)。首先,根据小数部分RFAB进行两点线性插值,然后根据小数部分θFB进行两点线性插值。不妨假设位置点K1,K2的扫描深度相同,特别的与点A的扫描深度相同。位置点P1和K1,P2和K2的扫描线号相同。位置点K1距P1在R轴方向的归一化距离是RF,以P1和Q1的轴向距离为1,RF由RFAB决定,即RF=RFAB,那么K1距离Q1点在轴方向距离是1-RF。假设A点对应的偏转极坐标和扫描线n之间的归一化夹角是θF,以扫描线n和扫描线n+1之间的夹角为1,θF由θFB决定,即θF=θFB,那么A点对应的极坐标和扫描线n+1之间的归一化夹角是1-θF。邻近点内插值算法以归一化的夹角和归一化距离作为权重,进行插值运算。
首先,根据P1和Q1两点的超声波信息插出K1点的信息XK1,如下式(12):
XK1=XP1×(1-RF)+XQ1×RF (12)
然后,根据P2和Q2两点的超声波信息插出K2点的信息XK2,如下式(13)
XK2=XP2×(1-RF)+XQ2×RF (13)
那么,位置点A对应的超声回波数据XA如下式(14):
XA=XK1×(1-θF)+XK2×θF (14)
经过插值处理后的超声回波数据输出到空间复合等处理模块后,空间复合模块把经过偏转图象扫描变换的多角度图像进行复合处理,空间复合控制模块提供选定加权复合的图象加权系数及加权手段。
以上给出了最常用的四点插值运算的例子,而除了四点插值外,还可以采用四点的倍数个邻近像素点、两个邻近像素点、六个邻近像素点等偶数个邻近象素点方式来进行插值运算。
基于上述扫描变换的方法,本发明还提供了一种超声成像系统中的凸阵空间复合方法,其包括以下步骤:
S1、对接收到的超声回波数据进行数字扫描变换前的预处理,这里的预处理包括:放大、波束合成、滤波、包络检波。
S2、利用凸阵偏转角度,将偏转角度极坐标系下的超声回波数据转换为笛卡儿坐标系下的回波数据,该转换过程包括以下步骤:查找至少两个与待求象素点最接近的邻近象素点;获取所述邻近象素点的超声回波数据,利用该超声回波数据通过插值运算,获得待求象素点的超声回波数据;
S3、将笛卡儿坐标系下的超声回波数据进行复合处理,获得用以显示的图像数据。这里可以采用平均复合法对不同偏转角度的图像进行加权平均,获得复合后的图像数据。
图像的空间复合有平均、最大值、中值、最小均方根等方法,平均复合方法为多角度空间复合的传统方法,思路就是对不同角度的图像进行加权平均,在抑制斑点噪声的同时突出信号特征。图像的缺陷如斑点、声影等具有一定的随机性,因此通过多角度图像的像素平均可以有效地抑制这些缺陷同时保留有用的信号。假定复合成像使用扫描角度数量为N,经过偏转图像扫描变换后各帧图像为Ii,复合结果为CompIi,显然有:
实时的空间复合通过把前一次复合图像结果减去第一幅图像结果再叠加新的图像而实现,其处理方法如下:
CompIi+1=CompIi+Ii+1-Ii-N+1
由于不同区域参与复合的帧数不同,并且有边缘线权值的影响,会导致复合图像的不同区域(主要是边缘)的灰度不同,所以参加空间复合的图像加权系数应随图像复合帧数的变化而进行相应的调整。复合后的图像经灰度映射直接送至显示器。
本发明的方法主要应用于超声成像系统上,凸阵空间复合超声系统基本框架图如图2所示。凸阵空间复合超声系统包括:空间复合控制模块、时间增益波长放大器、波束合成模块、包络检波模块、扫描变换模块、空间复合模块、显示模块、发射/接收模块。空间复合模块控制着整个空间复合过程中的发射,接收,扫描变换及空间复合算法等。
空间复合成像时,发射/接收模块先被切换到发射状态,系统通过发射脉冲控制模块控制发射脉冲的形状、延时以及参与发射的阵元,使发射的超声波聚焦到预定扫描线上的预定焦点位置。特别的,空间复合控制模块通过控制发射脉冲延时使发射扫描线聚焦于指定的偏转角度。
接着,所述发射/接收模块切换到接收状态,超声回波经超声换能器各阵元接收转换成电信号。该电信号经时间增益波长放大器放大,以补偿不同深度下的超声波衰减;再送往波束合成模块,调整各阵元回波的延时并进行变迹,以提高当前接收扫描线回波信号的信噪比。特别地,空间复合控制模块会相应地控制波束合成的延时使接收扫描线偏转到与发射扫描线相同的角度。这样就构成了偏转角度图象的发射与接收全过程,形成了具有一定偏转角度的图像,为空间复合做准备。
经波束合成后的回波信号一般还经过动态滤波提高信噪比,之后经包络检波模块获得回波信号的包络。所述的包络检波模块可以采用对正交解调后的两路正交信号取模检测的方法。检波后的数据送至扫描变换模块,该扫描变换模块利用空间复合控制模块提供出相应的偏转角度,把各个偏转角度的扫描线数据转换为代表真实空间位置的图像数据,再送至空间复合模块,即实现偏转凸阵图像的坐标变换和凸阵空间复合。
从上述陈述的本发明的方法过程可以看出,为了使超声成像系统只完成一次坐标变换就可以为空间复合做好准备,本发明的扫描变换模块中包括:查找模块,和插值运算模块,如图6所示,查找模块用于查找至少两个与待求象素点最接近的邻近象素点;插值运算模块用于利用所述邻近象素点的超声回波数据进行插值运算,获得待求象素点的超声回波数据。其中,查找模块包括:角度计算单元、距离计算单元和查找单元;
角度计算单元用于计算待求象素点在探头圆弧上所对应的射线起始点在极坐标系下的角度;
距离计算单元用于计算所述待求象素点到所述射线起始点的距离;及
查找单元用于根据所述角度和距离的计算结果,在笛卡儿坐标系中查找至少两个与待求象素点距离最相近的邻近象素点。
上述各个模块均可以通过软件编程来实现,其原理可以参见上述相关说明。
综上所述,本发明针对超声凸阵偏转成像时的扫描变换系统,通过计算待求像素点在极坐标下对应圆弧上射线起始点的角度和待求像素点距起始点的距离,确定与像素点最相近的点,并进行插值计算,获取偏转情况下的笛卡尔坐标系的超声图像,用于凸阵的空间复合处理。可见,本发明直接利用扫描变换下像素点之间所具有的三角函数关系,将偏转的极坐标图像一次性转换到笛卡尔坐标系,然后再进行空间复合处理,这种方法仅进行了一次坐标变换,即一次插值处理,能够有效的减少因坐标变换引入的插值噪声。并且,本发明采用选取邻近象素点进行插值计算,从而缩小了插值误差。本发明的邻近象素点和插值系数是通过计算极坐标系下对应探头圆弧上射线起始点位置后,把角度和轴向距离分解成整数部分和小数部分获取的,这样做可以保证精确定位待求象素点地邻近象素点,减少插值误差,而且使得计算简单。凸阵探头的图像采用本发明的凸阵空间复合方法及扫描变换方法后,能够在保留图像信号的同时,有效的提高信噪比,降低回声失落,增强组织间的分辨能力,清晰地显示组织间的边界,提高图像质量。
上述各具体步骤的举例说明较为具体,并不能因此而认为是对本发明的专利保护范围的限制,本发明的专利保护范围应以所附权利要求为准。