实用新型内容
本实用新型要解决的技术问题是提供一种数字扫描变换装置,克服现有技术的DSC插值装置不能根据图像特点调整选择不同插值方法的缺陷。
本实用新型为解决上述技术问题所采用的技术方案为:
一种数字扫描变换装置,包括线数据缓存器、帧图像存储器、扫描变换输入控制器、扫描变换输出控制器、坐标变换器,所述线数据缓存器与所述帧图像存储器相连,所述坐标变换器分别与所述扫描变换输出控制器和所述帧图像存储器相连,所述扫描变换输入控制器与所述帧图像存储器相连,还包括选择器和至少两个插值子模块,所述选择器的输入端分别与所述坐标变换器和所述帧图像存储器相连,所述选择器的输出端分别与每个插值子模块相连,由选择器控制信号控制所述选择器将所述坐标变换器的输出数据和所述帧图像存储器的输出数据送往其中一个插值子模块。
所述的数字扫描变换装置,其中所述插值子模块包括双线性插值子模块。
所述的数字扫描变换装置,其中所述插值子模块包括最邻近插值子模块。
所述的数字扫描变换装置,其中所述插值子模块包括拓展线性插值子模块。
所述的数字扫描变换装置,其中所述选择器设为三选一选择器。
所述的数字扫描变换装置,其中所述拓展线性插值子模块内存储插值平滑系数,由输入到所述拓展线性插值子模块的选择信号选择不同的插值平滑系数。
所述的数字扫描变换装置,其中所述选择器和所述插值子模块由现场可编程逻辑门阵列FPGA实现。
本实用新型的有益效果:本实用新型数字扫描变换装置对插值模块进行了改进,可以根据图像的特点来选择不同的数字扫描变换插值方法,优化了数字扫描变换的效果。
具体实施方式
下面根据附图和实施例对本实用新型作进一步详细说明:
如图3所示,本实用新型数字扫描变换装置,包括线数据缓存器、帧图像存储器、扫描变换输入控制器、扫描变换输出控制器、坐标变换器,所述线数据缓存器与所述帧图像存储器相连,所述坐标变换器分别与所述扫描变换输出控制器和所述帧图像存储器相连,所述扫描变换输入控制器与所述帧图像存储器相连,还包括选择器和至少两个插值子模块,所述选择器的输入端分别与所述坐标变换器和所述帧图像存储器相连,所述选择器的输出端分别与每个插值子模块相连,由选择器控制信号控制所述选择器将所述坐标变换器的输出数据和所述帧图像存储器的输出数据送往其中一个插值子模块。
如图4所示,在本实用新型的一个具体实施例中,插值子模块分别是双线性插值子模块、最邻近插值子模块和拓展线性插值子模块,选择器是三选一选择器,拓展线性插值子模块内存储插值平滑系数,由输入到所述拓展线性插值子模块的选择信号选择不同的插值平滑系数。
如图5所示,DSC中采样数据的存储模式共有Sr行,Sc列。每一行上的采样数据均对应相同的扫描深度,从第一行到第Sr行,对应的扫描深度依次等间隔从零增加到设定的探头最大扫描深度。每一列上的采样数据均对应相同的扫查角度,从第一列到第Sc列,对应的扫查角度从-θ0增加到θ0.一般来说,Sr为512,而Sc在整节距扫描的时候大小为80(80阵元)到128(128阵元);在半节距扫描的时候为160(80阵元)到256(128阵元).一般而言,图像区的大小为512*512,所以DSC需要对原始采样数据进行插值,才能得到用来显示的数据.
图6和表1中显示了DSC计算中需要的各相关参数。
表1 DSC相关参数
参数名称 |
参数的定义 |
符号 |
备注 |
凸阵探头的扫查角度 |
凸阵探头扫查时,所成扇形的角度 |
2θ0 |
单位为弧度 |
凸阵探头的曲率半径 |
曲率半径为凸阵探头的物理扇形半径 |
r0 |
单位为物理单位毫米(mm) |
凸阵探头当前的扫描深度 |
与采样数据对应的凸阵探头的扫描深度 |
D0 |
单位为物理单位毫米(mm) |
凸阵探头能够设定的最小扫描深度 |
凸阵探头能够使用的扫描深度中的最小值 |
D0 |
单位为物理单位毫米(mm) |
采样数据的行数 |
|
Sr |
无量纲的参数 |
采样数据的列数 |
|
Sc |
无量纲的参数 |
显示在屏幕上的图像区宽度 |
进行DSC之后显示在屏幕上的图像区宽度 |
W0 |
以像素为单位 |
显示在屏幕上的图像区高度 |
进行DSC之后显示在屏幕上的图像区高度 |
H0 |
同W0 |
垂距 |
探头圆心到坐标系x-y坐标轴x的垂直距离 |
d |
单位为物理单位毫米(mm) |
垂距的计算公式:
d=r0×cosθ0 (1)
坐标平移公式:
而
由式(2)和式(3)可得
式(4)便是最终用于DSC的坐标平移公式。
参照图6,在各个参数都确定的情况下,设从x-y坐标系下的坐标(x,y)经过坐标平移之后对应的在u-v坐标系下的坐标为(u,v),而坐标(u,v)经过坐标变换后对应在R-θ坐标系下的坐标为(R,θ),可得坐标(u,v)到坐标R-θ的变换公式如式(5)所示:
极径量化因子为qr,极角量化因子为qθ,则量化后的极坐标为(Rq,θq):
如图7所示,为了让量化后的极坐标和采样数据的编号能够统一,需要给极径和极角去偏,量化去偏后的极坐标为(Rqf,θqf)。其中极径和极角都可以由整数部分及小数部分组成。设Rqf的整数部分为i,小数部分为α,设θqf的整数部分为j,小数部分为β。最邻近插值和双线性插值法都是基于被插值点周围的四个点Pi,j,Pi,j+1,Pi+1,j,Pi+1,j+1来进行插值的。
最邻近插值法中P的值为Pi,j,Pi,j+1,Pi+1,j,Pi+1,j+1四点中与P点距离最近的点的值。
双线性插值的计算公式为:
P=(1-α)×(1-β)×Pi,j+(1-α)×β×Pi,j+1+α×(1-β)×Pi+1,j+α×β×Pi+1,j+1 (9)
双线性插值的计算量适中,且效果也不错,只是会造成图像的模糊,且在某些情况下会出现锯齿.一般来说,Sr为512,而Sc最小为80;而图像区的大小为512*512.所以在图像的纵向,不需要再优化插值方法,而在图像的横向,因为声线数的不足,需要优化插值方法,在本实用新型中,使用了一种基于8点的拓展线性插值法,这是一种改良后的线性插值方法.根据超声声线数较少的特点,在横向采用了四点插值,而在纵向采用了两点插值,如图8所示.
在计算过程中,先根据纵向的两个点通过线性插值得到一个点,八个点可以得到四个点,分别是Pj-1,Pj,Pj+1,Pj+2,其计算公式如式(10)所示。再通过横向四个点的插值得到P点的值,其计算公式如式(11)所示。
Pj-1=(1-α)×Pi,j-1+α×Pi+1,j-1
Pj=(1-α)×Pi,j+α×Pi+1,j (10)
Pj+1=(1-α)×Pi,j+1+α×Pi+1,j+1
Pj+2=(1-α)×Pi,j+2+α×Pi+1,j+2
P=kj-1×Pj-1+kj×Pj+kj+1×Pj+1+kj+2×Pj+2 (11)
横向四点插值系数的计算公式为:
其中d为各点距P点的横向距离,而η是插值平滑系数,η介于0-1之间,η越大插值越锐利,η越小插值越平滑。Pj-1到P点的横向距离是1+β,Pj到P点的横向距离是β,Pj+1到P点的横向距离是1-β,Pj+2到P点的横向距离是2-β。
拓展线性插值子模块以如下方式来进行拓展线性插值:
1)根据插值平滑系数η,及公式(12),计算得出两个插值系数查找表;
2)因为Pj,Pj+1距P点的距离在0和1之间,而Pj-1,Pj+2距P点的距离在1和2之间,两两之间使用的计算公式相同;
3)根据显示坐标x,y求得i,j,α,β;
4)从数据存储器中读取所需要的8个点;
5)有8个点及α求得Pj-1,Pj,Pj+1,Pj+2;
6)根据β,由查找表得到插值系数kj-1,kj,kj+1,kj+2;
7)由Pj-1,Pj,Pj+1,Pj+2及kj-1,kj,kj+1,kj+2求得P点。
本实用新型数字扫描变换装置的实施例可由现场可编程逻辑门阵列FPGA实现:数据读地址及α,β由显示坐标求出。插值方法可由用户选择,控制信号DSCSEL负责插值方法的三选一,其中“00”代表双线性插值,“01”代表最邻近插值,“10”代表拓展线性插值。插值平滑系数η的默认值为0.2;用户可以选择0.1-0.5共五档,在硬件中由控制信号SMSEL选择η的值。每种插值方法都由不同的插值子模块实现,在选定一插值方法后,数据及α,β被输入到该插值子模块中,经过相应的插值计算后得到显示数据,最终输出到显示模块。
本领域技术人员不脱离本实用新型的实质和精神,可以有多种变形方案实现本实用新型,以上所述仅为本实用新型较佳可行的实施例而已,并非因此局限本实用新型的权利范围,凡运用本实用新型说明书及附图内容所作的等效结构变化,均包含于本实用新型的权利范围之内。