发明名称 一种合成聚焦的便携式B型超声成像方法
摘要 一种合成聚焦的便携式B型超声成像方法,属超声医学影像学技术领域,便携式B超设备的射频回波信号在模/数转换处理中被随机采样,然后进行波束形成、解调检波处理而形成欠采样的扫描线,最后采用压缩感知的信号重构方法对欠采样的扫描线进行信号恢复后显像。本发明方法的步骤为:预处理、产生伪随机脉冲序列和测量矩阵、换能器的一个阵元发射超声波、换能器的所有阵元都接收超声射频回波信号、随机采样的模/数转换、波束形成、解调检波、对数压缩、压缩感知的信号重构、帧缓存、显像。本发明解决了传统的合成聚焦波束形成方法的系统复杂度高、数据规模庞大、难以在便携式B超设备中得到实际应用的问题,实现了合成聚焦的便携式B超的成像方法。
申请公布号 CN102920478B 申请公布日期 2014.04.16
申请号 CN201210480709.7 申请日期 2012.11.22
申请人 山东大学 发明人 孙丰荣;杜衍震;李凯一;金鑫
分类号 A61B8/00(2006.01)I 主分类号 A61B8/00(2006.01)I
代理机构 济南金迪知识产权代理有限公司 37219 代理人 许德山
主权项 1.一种合成聚焦的便携式B型超声成像方法,便携式B型超声设备的射频回波信号在模/数转换处理中被随机采样,然后进行波束形成、解调检波处理而形成欠采样的扫描线,最后采用压缩感知的信号重构方法对欠采样的扫描线进行信号恢复后显像;假定便携式B型超声设备有定点聚焦和动态聚焦两种聚焦模式,探查深度为<img file="DEST_PATH_IMAGE001.GIF" wi="18" he="18" />,每条扫描线包含<img file="479033DEST_PATH_IMAGE002.GIF" wi="22" he="18" />个像素,换能器的阵元个数为<img file="DEST_PATH_IMAGE003.GIF" wi="20" he="20" />,编号依次为<img file="780701DEST_PATH_IMAGE004.GIF" wi="68" he="22" />;记选通阵元索引指示器为<img file="DEST_PATH_IMAGE005.GIF" wi="14" he="18" />,其值为换能器的阵元编号;本方法步骤如下:S1)预处理确定便携式B型超声设备的聚焦模式,并据此确定模/数转换器对原始超声射频回波信号开始采样的延时;确定欠采样率<img file="954325DEST_PATH_IMAGE006.GIF" wi="14" he="18" />;确定压缩感知的信号重构方法中的稀疏化变换<img file="DEST_PATH_IMAGE007.GIF" wi="20" he="18" />,该变换为离散的小波变换,或者为离散的余弦变换;初始化选通阵元索引指示器<img file="279127DEST_PATH_IMAGE005.GIF" wi="14" he="18" />的值为1;S2)产生伪随机脉冲序列和测量矩阵便携式B型超声设备的主控芯片产生伪随机脉冲序列,记作<img file="648928DEST_PATH_IMAGE008.GIF" wi="26" he="28" />,其第<img file="726475DEST_PATH_IMAGE010.GIF" wi="9" he="18" />个脉冲的到达时刻记作<img file="DEST_PATH_IMAGE011.GIF" wi="16" he="25" />;测量矩阵为压缩感知的信号重构方法中所使用的测量矩阵,记作<img file="523529DEST_PATH_IMAGE012.GIF" wi="18" he="17" />,其结构有如下特征:(1)测量矩阵<img file="70048DEST_PATH_IMAGE012.GIF" wi="18" he="17" />的列数为<img file="243541DEST_PATH_IMAGE002.GIF" wi="22" he="18" />,行数记为<img file="DEST_PATH_IMAGE013.GIF" wi="18" he="18" />,<img file="942638DEST_PATH_IMAGE014.GIF" wi="68" he="22" />,(2)<img file="910594DEST_PATH_IMAGE012.GIF" wi="18" he="17" />的每行只有一个元素的值为1,其余元素的值都为0,且<img file="6726DEST_PATH_IMAGE012.GIF" wi="18" he="17" />的第一行、第一列元素的值为1,(3)<img file="656013DEST_PATH_IMAGE012.GIF" wi="18" he="17" />的第<img file="DEST_PATH_IMAGE015.GIF" wi="14" he="21" />行值为1的元素所在的列为<img file="521201DEST_PATH_IMAGE016.GIF" wi="113" he="49" />,其中<img file="DEST_PATH_IMAGE017.GIF" wi="18" he="26" />为伪随机脉冲序列第<img file="847009DEST_PATH_IMAGE015.GIF" wi="14" he="21" />个脉冲的到达时刻,<img file="430437DEST_PATH_IMAGE018.GIF" wi="16" he="25" />为伪随机脉冲序列第1个脉冲的到达时刻,<img file="DEST_PATH_IMAGE019.GIF" wi="14" he="16" />为超声波在探测介质中的传播速度,<img file="617836DEST_PATH_IMAGE001.GIF" wi="18" he="18" />为便携式B型超声设备的探查深度,<img file="603109DEST_PATH_IMAGE002.GIF" wi="22" he="18" />为每条扫描线包含的像素个数;S3)换能器的一个阵元发射超声波便携式B型超声设备换能器的一个阵元在高压脉冲的激励下发射超声波,该阵元的编号为<img file="912868DEST_PATH_IMAGE005.GIF" wi="14" he="18" />的值;上述高压脉冲是便携式B型超声设备的主控芯片所产生的高频激励脉冲经高压脉冲产生电路得到;S4)换能器的所有阵元都接收超声射频回波信号便携式B型超声设备换能器的所有阵元都接收超声射频回波信号;对每个阵元所接收到的超声射频回波信号都进行前置放大、时间增益补偿前端处理,获得原始超声射频回波信号集<img file="660289DEST_PATH_IMAGE020.GIF" wi="112" he="28" />,其中<img file="DEST_PATH_IMAGE021.GIF" wi="121" he="25" />为第<img file="713695DEST_PATH_IMAGE022.GIF" wi="14" he="18" />个阵元接收到的、经上述前端处理后的原始超声射频回波信号;S5)随机采样的模/数转换对原始超声射频回波信号<img file="DEST_PATH_IMAGE023.GIF" wi="97" he="25" />进行模/数转换处理;该模/数转换处理是基于随机采样的;该随机采样是由步骤S2)产生的伪随机脉冲序列<img file="491158DEST_PATH_IMAGE008.GIF" wi="26" he="28" />控制模/数转换器对每个原始超声射频回波信号的采样时刻实现的,具体为:模/数转换器对第<img file="971818DEST_PATH_IMAGE010.GIF" wi="9" he="18" />个原始超声射频回波信号<img file="451210DEST_PATH_IMAGE024.GIF" wi="21" he="25" />的第<img file="DEST_PATH_IMAGE025.GIF" wi="14" he="20" />个采样点的采样时刻为<img file="308308DEST_PATH_IMAGE026.GIF" wi="42" he="21" />,其中<img file="DEST_PATH_IMAGE027.GIF" wi="18" he="25" />为伪随机脉冲序列<img file="940277DEST_PATH_IMAGE008.GIF" wi="26" he="28" />第<img file="591838DEST_PATH_IMAGE025.GIF" wi="14" he="20" />个脉冲的到达时刻,<img file="325570DEST_PATH_IMAGE028.GIF" wi="17" he="21" />为模/数转换器对第<img file="720780DEST_PATH_IMAGE010.GIF" wi="9" he="18" />个原始超声射频回波信号<img file="269573DEST_PATH_IMAGE024.GIF" wi="21" he="25" />开始采样的延时;原始超声射频回波信号经随机采样的模/数转换处理后,获得欠采样的超声射频回波信号集<img file="DEST_PATH_IMAGE029.GIF" wi="96" he="28" />,其中信号<img file="295298DEST_PATH_IMAGE030.GIF" wi="308" he="26" />为第<img file="562331DEST_PATH_IMAGE010.GIF" wi="9" he="18" />个原始超声射频回波信号<img file="948182DEST_PATH_IMAGE024.GIF" wi="21" he="25" />经上述随机采样的模/数转换处理后获得的欠采样的超声射频回波信号,<img file="DEST_PATH_IMAGE031.GIF" wi="148" he="26" />为该欠采样的超声射频回波信号的第<img file="351481DEST_PATH_IMAGE032.GIF" wi="14" he="16" />个离散信号点,<img file="610424DEST_PATH_IMAGE013.GIF" wi="18" he="18" />为步骤S2)所述测量矩阵<img file="302437DEST_PATH_IMAGE012.GIF" wi="18" he="17" />的行数;S6)波束形成对欠采样的超声射频回波信号<img file="DEST_PATH_IMAGE033.GIF" wi="81" he="25" />进行叠加求和,得到欠采样的超声射频信号线<img file="39448DEST_PATH_IMAGE034.GIF" wi="92" he="28" />,其中<img file="DEST_PATH_IMAGE035.GIF" wi="212" he="46" />为欠采样的超声射频信号线<img file="251249DEST_PATH_IMAGE036.GIF" wi="16" he="18" />的第<img file="681094DEST_PATH_IMAGE015.GIF" wi="14" he="21" />个离散信号点,<img file="594823DEST_PATH_IMAGE013.GIF" wi="18" he="18" />为上述测量矩阵<img file="135526DEST_PATH_IMAGE012.GIF" wi="18" he="17" />的行数;S7)解调检波、对数压缩对欠采样的射频信号线<img file="778997DEST_PATH_IMAGE036.GIF" wi="16" he="18" />进行解调检波处理得到到欠采样的扫描线,也称作欠采样的声线,或欠采样的视频信号线;然后对欠采样的扫描线进行对数压缩以减小输出动态范围,经此对数压缩处理的欠采样的扫描线记作<img file="301114DEST_PATH_IMAGE038.GIF" wi="9" he="20" />;S8)压缩感知的信号重构采用压缩感知的信号重构方法,对欠采样的扫描线<img file="30035DEST_PATH_IMAGE038.GIF" wi="9" he="20" />进行信号恢复;压缩感知的信号重构方法中的测量矩阵是上述矩阵<img file="108850DEST_PATH_IMAGE012.GIF" wi="18" he="17" />,稀疏化变换为<img file="544510DEST_PATH_IMAGE007.GIF" wi="20" he="18" />;待恢复扫描线记作<img file="DEST_PATH_IMAGE039.GIF" wi="14" he="16" />,在稀疏化变换下的稀疏化表示记作<img file="50578DEST_PATH_IMAGE040.GIF" wi="13" he="16" />,其表达式为<img file="DEST_PATH_IMAGE041.GIF" wi="50" he="20" />;首先通过对欠采样的扫描线<img file="955211DEST_PATH_IMAGE038.GIF" wi="9" he="20" />的表达公式<img file="837717DEST_PATH_IMAGE042.GIF" wi="60" he="20" />来求解得到待恢复扫描线<img file="127884DEST_PATH_IMAGE039.GIF" wi="14" he="16" />的稀疏化表示<img file="70432DEST_PATH_IMAGE040.GIF" wi="13" he="16" />,然后再利用稀疏化表达式<img file="508366DEST_PATH_IMAGE041.GIF" wi="50" he="20" />得到待恢复扫描线<img file="115934DEST_PATH_IMAGE039.GIF" wi="14" he="16" />;待恢复扫描线<img file="588504DEST_PATH_IMAGE039.GIF" wi="14" he="16" />是<img file="436374DEST_PATH_IMAGE002.GIF" wi="22" he="18" />维的;S9)帧缓存将恢复后的扫描线<img file="361605DEST_PATH_IMAGE039.GIF" wi="14" he="16" />送入帧缓存进行存储;令选通阵元索引指示器<img file="523596DEST_PATH_IMAGE005.GIF" wi="14" he="18" />的值加1;判断选通阵元索引指示器<img file="850672DEST_PATH_IMAGE005.GIF" wi="14" he="18" />的值是否为<img file="DEST_PATH_IMAGE043.GIF" wi="38" he="20" />,若选通阵元索引指示器<img file="823438DEST_PATH_IMAGE005.GIF" wi="14" he="18" />的值为<img file="235965DEST_PATH_IMAGE043.GIF" wi="38" he="20" />,则转到步骤S10),否则返回步骤S2);S10)显像便携式B型超声设备显示一帧图像。
地址 250100 山东省济南市历城区山大南路27号