CN113252793A - 基于掩膜压缩的单传感器三维超声成像装置及方法 - Google Patents
基于掩膜压缩的单传感器三维超声成像装置及方法 Download PDFInfo
- Publication number
- CN113252793A CN113252793A CN202110481195.6A CN202110481195A CN113252793A CN 113252793 A CN113252793 A CN 113252793A CN 202110481195 A CN202110481195 A CN 202110481195A CN 113252793 A CN113252793 A CN 113252793A
- Authority
- CN
- China
- Prior art keywords
- mask
- coding mask
- ultrasonic probe
- ultrasonic
- coding
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/06—Visualisation of the interior, e.g. acoustic microscopy
- G01N29/0654—Imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/22—Details, e.g. general constructional or apparatus details
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/01—Indexing codes associated with the measuring variable
- G01N2291/012—Phase angle
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/10—Number of transducers
- G01N2291/101—Number of transducers one transducer
Landscapes
- Physics & Mathematics (AREA)
- Analytical Chemistry (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Signal Processing (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
基于掩膜压缩的单传感器三维超声成像装置及方法,包括非聚焦水浸超声探头,非聚焦水浸超声探头固定于旋转平台的左端,旋转平台的右端固定有编码掩膜,旋转平台连接在旋转平台基座上,旋转平台基座装夹在卡盘上,卡盘固定于水箱底部;编码掩膜和样件相对,样件固定于样件基座顶部,样件基座固定于水箱底部;方法是非聚焦水浸超声探头激励超声波,编码掩膜进行差异化后传播至样件处;声学信号被样件反射,经由编码掩膜压缩后,单一的反射回波信号被非聚焦水浸超声探头接收;进行多次旋转编码掩膜,得到多组接收信号和压缩矩阵,进行目标图像的重构;本发明具有成本低、结构简单、成像分辨率高的优点。
Description
技术领域
本发明涉及超声成像技术领域,具体涉及基于掩膜压缩的单传感器三维超声成像装置及方法。
技术背景
超声作为一种成熟的检测和成像手段,在工业界应用十分广泛。目前主流的超声成像方式有扫描超声成像和相控阵超声成像两种。扫描超声成像检测方式灵活,硬件复杂度低,但是依赖于机械扫描,声束可达性差,耗费时间长;相控阵超声成像通过控制换能器阵列中各阵元发射脉冲的延迟时间,改变声波到达物体内某点时的相位关系,实现超声波的波束扫描、偏转和聚焦,然后由换能器阵列接收物体反射回的脉冲回波信号,通过延迟处理进行成像,超声相控阵检测效率高,具有多种显示模式,但其硬件复杂度和集成度很高,还要求比较高的数字控制和信号处理技术,成本过于高昂。
发明内容
为了克服上述现有技术的缺点,本发明的目的在于提供一种基于掩膜压缩的单传感器三维超声成像装置及方法,使用一个非聚焦水浸超声探头和一个简单的编码掩膜对目标样件进行成像,具有成本低、结构简单、成像分辨率高的优点。
为了达到上述目的,本发明采取的技术方案为:
一种基于掩膜压缩的单传感器三维超声成像装置,包括非聚焦水浸超声探头1,非聚焦水浸超声探头1固定于旋转平台2的左端,旋转平台2的右端固定有编码掩膜4,旋转平台2连接在旋转平台基座8上,旋转平台基座8装夹在卡盘7上,卡盘7固定于水箱9底部;
编码掩膜4和样件5相对,样件5固定于样件基座6顶部,样件基座6固定于水箱9底部。
所述的水箱9采用有机玻璃制成,内部放置有水。
所述的编码掩膜4固定于旋转平台2的中心孔中,编码掩膜4一个端面设有位置随机、直径1mm,深度随机分布在0-1mm范围的孔。
所述的旋转平台2包括旋转盘11,旋转盘11一周激光雕刻有分度为1°的刻度,能够精确旋转360°;旋转盘11的中心孔固定有编码掩膜4,中心孔外的旋转盘11表面连接有圆盘3;旋转盘11连接在基板13上,并通过锁定螺丝12锁定;旋转盘11和基板13之间设置有旋转千分尺14;在锁定螺丝12处于解锁状态时,旋转圆盘3带动旋转盘11精确旋转固定角度,实现编码掩膜4精确旋转;当锁定螺丝12处于锁定状态时,旋转千分尺14可以精细调节。
利用一种基于掩膜压缩的单传感器三维超声成像装置的方法,包括以下步骤:
1)非聚焦水浸超声探头1激励超声波s(t),编码掩膜4将非聚焦水浸超声探头1产生的相对均一的声束进行差异化,差异化后的声束继续传播至样件5处;
2)空间中相应像素点的声学信号up(t)被成像区域内的样件5反射,经由编码掩膜4压缩后,单一的反射回波信号y(t)被非聚焦水浸超声探头1接收;
3)进行多次旋转编码掩膜4,得到多组接收信号y和压缩矩阵A,进行目标图像的重构。
所述的步骤1)具体为:非聚焦水浸超声探头1在脉冲收发仪激励下发射超声波s(t),超声波的相位均匀性被编码掩膜4破坏,编码掩模4近似是一个点源的集合,这些点源位于编码掩模4各通道的表面,编码掩膜4表面的超声场认为是这些点源集合产生的声场的叠加,每个虚拟点源被激发的时间点是超声场从非聚焦水浸超声探头1表面位置通过编码掩膜4到达编码掩膜4末端的时间决定的,通道i被激发的时间点ti为:
其中di为通道i的高度,cacry为超声波在有机玻璃中的传播速度;
则通道i发出的超声信号fi(t)为:
空间中像素点p处的声学信号为up(t),每一个通道i对像素点p处的声场信号的贡献up,i为:
up,i=s(t-ti-tp) (3)
其中tp为从超声信号从编码掩膜通道i表面传播至像素点p处耗费的时间;
编码掩膜4共有I个通道,则像素点p处的声信号up(t)为:
所述的步骤2)具体为:空间中共有P个像素点存在散射体,则通道i接收得到的散射信号yi(t)是检测区域内所有散射体回波信号的叠加:
非聚焦水浸超声探头1接收到的信号是每个通道接收的信号的总和:
用yi表示通道i接收并输出到非聚焦水浸超声探头1的成像区域内所有散射体的散射信号;用ai,p表示散射强度为1时,像素点p散射回通道i的散射信号;用向量x表示成像区域内像素点散射强度组成的列向量,则
上式转化为矩阵形式
yi=[ai,1ai,2…ai,P]x=Aix (8)
非聚焦水浸超声探头1实际接收到的信号y近似是各通道信号的总和,
其中A被称为压缩矩阵。
所述的步骤3)具体为:旋转k次编码掩膜4得到了k组观测信号y(k)和k组压缩矩阵A(k),k组观测信号y(k)组合得到yassem,k组压缩矩阵组合得到Aassem,二者关系为:
本发明的有益效果为:
本发明使用非聚焦水浸超声探头1在脉冲收发仪激励下发射高频超声脉冲,抛弃阵列聚焦的方案会降低系统的硬件复杂度,同时又不存在超声扫描成像扫查时间过长的问题。本发明使用的编码掩膜4端面设有位置随机、深度在0-1mm范围内随机的孔,设计并加工完成后,掩膜结构4固定且已知。由于超声波在有机玻璃和水中的传播速度不同,掩膜结构4厚度的局部变化会导致固定的局部延迟,从而扰乱声场的相位。传统超声成像方法的相位均匀性被破坏后会导致图像出现伪影,使成像结果变得模糊,但在本发明中这种干扰模式是已知的,在超声波的发射和接收中破坏声场相位均匀性有益于降低像素点回波信号的相关性,改善图像分辨率。本发明具有结构简单、成本低、成像分辨高的优点。
附图说明
图1是本发明装置的结构示意图。
图2是本发明装置编码掩膜的结构示意图。
图3是本发明装置旋转平台的结构示意图。
图4是本发明装置圆盘的结构示意图。
图5是本发明装置样件的结构示意图。
图6是实施例样件的成像结果。
具体实施方式
下面结合附图和实施例对本发明作进一步的详细描述。
参照图1,一种基于掩膜压缩的单传感器三维超声成像装置,包括非聚焦水浸超声探头1,非聚焦水浸超声探头1固定于旋转平台2的左端,旋转平台2的右端固定有编码掩膜4,旋转平台2连接在旋转平台基座8上,旋转平台基座8装夹在卡盘7上,卡盘7固定于水箱9底部;
编码掩膜4和样件5相对,样件5固定于样件基座6顶部,样件基座6固定于水箱9底部。
所述的水箱9采用有机玻璃制成,内部放置有水。
参照图1及图2,所述的编码掩膜4固定于旋转平台2的中心孔中,编码掩膜4一个端面设有位置随机、直径1mm,深度随机分布在0-1mm范围的孔。
参照图1、图3及图4,所述的旋转平台2包括旋转盘11,旋转盘11一周激光雕刻有分度为1°的刻度,能够精确旋转360°,旋转盘11的中心孔固定有编码掩膜4,中心孔外的旋转盘11表面通过螺栓连接有圆盘3;旋转盘11连接在基板13上,并通过锁定螺丝12锁定;旋转盘11和基板13之间设置有旋转千分尺14;在锁定螺丝12处于解锁状态时,旋转圆盘3带动旋转盘11精确旋转固定角度,实现编码掩膜4的精确旋转;当锁定螺丝12处于锁定状态时,旋转千分尺14可以精细调节。
参照图1、图5,所述的样件5固定于样件基座6顶面的方槽中。
利用一种基于掩膜压缩的单传感器三维超声成像装置的方法,包括以下步骤:
1)非聚焦水浸超声探头1激励超声波s(t),编码掩膜4将非聚焦水浸超声探头1产生的相对均一的声束进行差异化,差异化后的声束继续传播至样件5处;具体为:
非聚焦水浸超声探头1在脉冲收发仪激励下发射超声波s(t),超声波的相位均匀性被编码掩膜4破坏,编码掩模4近似是一个点源的集合,这些点源位于编码掩模4各通道的表面,编码掩膜4表面的超声场可以认为是这些点源集合产生的声场的叠加,每个虚拟点源被激发的时间点是超声场从非聚焦水浸超声探头1表面位置通过编码掩膜4到达编码掩膜4末端的时间决定的,通道i被激发的时间点ti为:
其中di为通道i的高度,cacry为超声波在有机玻璃中的传播速度;
则通道i发出的超声信号fi(t)为:
空间中像素点p处的声学信号为up(t),每一个通道i对像素点p处的声场信号的贡献up,i为:
up,i=fi(t-tp)=s(t-ti-tp) (3)
其中tp为从超声信号从编码掩膜通道i表面传播至像素点p处耗费的时间;
编码掩膜4共有I个通道,则像素点p处的声学信号up(t)为:
2)空间中相应像素点的声学信号up(t)被成像区域内的样件5反射,经由编码掩膜4压缩后,单一的反射回波信号y(t)被非聚焦水浸超声探头1接收;具体为:
空间中共有P个像素点存在散射体,则通道i接收得到的散射信号yi(t)是检测区域内所有散射体回波信号的叠加:
非聚焦水浸超声探头1接收到的信号是每个通道接收的信号的总和:
用yi表示通道i接收并输出到非聚焦水浸超声探头1的成像区域内所有散射体的散射信号;用ai,p表示散射强度为1时,像素点p散射回通道i的散射信号;用向量x表示成像区域内像素点散射强度组成的列向量,则
上式转化为矩阵形式:
yi=[ai,1ai,2…ai,P]x=Aix (8)
由于编码掩膜4前的非聚焦水浸超声探头1可以有效地将整个超声场集成在其表面,因此非聚焦水浸超声探头1实际接收到的信号y近似是各通道信号的总和:
其中A被称为压缩矩阵;
3)进行多次旋转编码掩膜4,得到多组接收信号y和压缩矩阵A,进行目标图像的重构;具体为:
为了降低各像素回波之间的相关性,需要旋转非聚焦水浸超声探头1前的编码掩膜4,这样干涉声场随之改变,从而获得包含新信息的额外测量值;旋转k次编码掩膜4得到了k组观测信号y(k)和k组压缩矩阵A(k),k组观测信号y(k)组合得到yassem,k组压缩矩阵组合得到Aassem,二者关系为:
参照图6,(a)是样件5的原始图像,(b)是编码掩膜4不旋转时的成像结果,(c)是编码掩膜4每次旋转90°,旋转4次的成像结果,(d)是编码掩膜4每次旋转9°,旋转40次的成像结果。当编码掩膜4不旋转时,重建的图像完全无法辨识样件5,当编码掩膜4旋转多次时,可以较为清晰地观察到样件5形貌。相较于编码掩膜4旋转4次,编码掩膜4旋转40次可以更为准确的提取样件5信息,模糊区域减少,增加编码掩膜4的旋转次数可以有效提高图像重构的质量。
Claims (8)
1.一种基于掩膜压缩的单传感器三维超声成像装置,包括非聚焦水浸超声探头(1),其特征在于:非聚焦水浸超声探头(1)固定于旋转平台(2)的左端,旋转平台(2)的右端固定有编码掩膜(4),旋转平台(2)连接在旋转平台基座(8)上,旋转平台基座(8)装夹在卡盘(7)上,卡盘(7)固定于水箱(9)底部;
编码掩膜(4)和样件(5)相对,样件(5)固定于样件基座(6)顶部,样件基座(6)固定于水箱(9)底部。
2.根据权利要求1所述的一种基于掩膜压缩的单传感器三维超声成像装置,其特征在于:所述的水箱(9)采用有机玻璃制成,内部放置有水。
3.根据权利要求1所述的一种基于掩膜压缩的单传感器三维超声成像装置,其特征在于:所述的编码掩膜(4)固定于旋转平台(2)的中心孔中,编码掩膜(4)一个端面设有位置随机、直径1mm,深度随机分布在0-1mm范围的孔。
4.根据权利要求1所述的一种基于掩膜压缩的单传感器三维超声成像装置,其特征在于:所述的旋转平台(2)包括旋转盘(11),旋转盘(11)一周激光雕刻有分度为1°的刻度,能够精确旋转360°,旋转盘(11)的中心孔固定有编码掩膜(4),中心孔外的旋转盘(11)表面连接有圆盘(3);旋转盘(11)连接在基板(13)上,并通过锁定螺丝(12)锁定;旋转盘(11)和基板(13)之间设置有旋转千分尺(14);在锁定螺丝(12)处于解锁状态时,旋转圆盘(3)带动旋转盘(11)精确旋转固定角度,实现编码掩膜(4)精确旋转;当锁定螺丝(12)处于锁定状态时,旋转千分尺(14)实现精细调节。
5.利用一种基于掩膜压缩的单传感器三维超声成像装置的方法,其特征在于,包括以下步骤:
1)首先非聚焦水浸超声探头(1)激励超声波s(t),编码掩膜(4)将非聚焦水浸超声探头(1)产生的相对均一的声束进行差异化,差异化后的声束继续传播至样件(5)处;
2)然后空间中相应像素点的声学信号up(t)被成像区域内的样件(5)反射,经由编码掩膜(4)压缩后,单一的反射回波信号y(t)被非聚焦水浸超声探头(1)接收;
3)最后进行多次旋转编码掩膜(4),得到多组接收信号y和压缩矩阵A,进行目标图像的重构。
6.根据权利要求5所述的方法,其特征在于,所述的步骤1)具体为:非聚焦水浸超声探头(1)在脉冲收发仪激励下发射超声波s(t),超声波的相位均匀性被编码掩膜(4)破坏,编码掩模(4)近似是一个点源的集合,这些点源位于编码掩模(4)各通道的表面,编码掩膜(4)表面的超声场认为是这些点源集合产生的声场的叠加,每个虚拟点源被激发的时间点是超声场从非聚焦水浸超声探头(1)表面位置通过编码掩膜(4)到达编码掩膜(4)末端的时间决定的,通道i被激发的时间点ti为:
其中di为通道i的高度,cacry为超声波在有机玻璃中的传播速度;
则通道i发出的超声信号fi(t)为:
空间中像素点p处的声学信号为up(t),每一个通道i对像素点p处的声场信号的贡献up,i为:
up,i=s(t-ti-tp) (3)
其中tp为从超声信号从编码掩膜通道i表面传播至像素点p处耗费的时间;
编码掩膜(4)共有I个通道,则像素点p处的声信号up(t)为:
7.根据权利要求6所述的方法,其特征在于,所述的步骤2)具体为:空间中共有P个像素点存在散射体,则通道i接收得到的散射信号yi(t)是检测区域内所有散射体回波信号的叠加:
非聚焦水浸超声探头(1)接收到的信号是每个通道接收的信号的总和:
用yi表示通道i接收并输出到非聚焦水浸超声探头(1)的成像区域内所有散射体的散射信号;用ai,p表示散射强度为1时,像素点p散射回通道i的散射信号;用向量x表示成像区域内像素点散射强度组成的列向量,则
上式转化为矩阵形式
yi=[ai,1 ai,2 … ai,P]x=Aix (8)
非聚焦水浸超声探头(1)实际接收到的信号y近似是各通道信号的总和,
其中A被称为压缩矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110481195.6A CN113252793A (zh) | 2021-04-30 | 2021-04-30 | 基于掩膜压缩的单传感器三维超声成像装置及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110481195.6A CN113252793A (zh) | 2021-04-30 | 2021-04-30 | 基于掩膜压缩的单传感器三维超声成像装置及方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113252793A true CN113252793A (zh) | 2021-08-13 |
Family
ID=77223698
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110481195.6A Pending CN113252793A (zh) | 2021-04-30 | 2021-04-30 | 基于掩膜压缩的单传感器三维超声成像装置及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113252793A (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102068277A (zh) * | 2010-12-14 | 2011-05-25 | 哈尔滨工业大学 | 基于压缩感知的单阵元多角度观测光声成像装置及方法 |
CN106067162A (zh) * | 2016-06-03 | 2016-11-02 | 西安电子科技大学 | 集成成像超分辨率微单元图像阵列采集与重构方法 |
CN111521683A (zh) * | 2020-05-06 | 2020-08-11 | 北京工业大学 | 一种基于多阵元等幅同步激励的材料缺陷超声三维成像方法 |
US20210088640A1 (en) * | 2019-09-24 | 2021-03-25 | Technion Research & Development Foundation Limited | Coded aperture ultrasound detector |
-
2021
- 2021-04-30 CN CN202110481195.6A patent/CN113252793A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102068277A (zh) * | 2010-12-14 | 2011-05-25 | 哈尔滨工业大学 | 基于压缩感知的单阵元多角度观测光声成像装置及方法 |
CN106067162A (zh) * | 2016-06-03 | 2016-11-02 | 西安电子科技大学 | 集成成像超分辨率微单元图像阵列采集与重构方法 |
US20210088640A1 (en) * | 2019-09-24 | 2021-03-25 | Technion Research & Development Foundation Limited | Coded aperture ultrasound detector |
CN111521683A (zh) * | 2020-05-06 | 2020-08-11 | 北京工业大学 | 一种基于多阵元等幅同步激励的材料缺陷超声三维成像方法 |
Non-Patent Citations (2)
Title |
---|
PIETER KRUIZINGA等: "Compressive 3D ultrasound imaging using a single sensor", 《SCIENCE ADVANCES》 * |
PIM VAN DER MEULEN等: "Coding Mask Design for Single Sensor Ultrasound Imaging", 《IEEE TRANSACTIONS ON COMPUTATIONAL IMAGING》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Montaldo et al. | Coherent plane-wave compounding for very high frame rate ultrasonography and transient elastography | |
Lu | 2D and 3D high frame rate imaging with limited diffraction beams | |
Ylitalo et al. | Ultrasound synthetic aperture imaging: monostatic approach | |
Flax et al. | Phase-aberration correction using signals from point reflectors and diffuse scatterers: Basic principles | |
JP4717995B2 (ja) | 超音波ビーム経路の数値的最適化方式 | |
US5678554A (en) | Ultrasound transducer for multiple focusing and method for manufacture thereof | |
US8111810B2 (en) | Method for producing highly constrained ultrasound images | |
EP0066343B1 (en) | Method and apparatus for measuring ultrasonic attenuation characteristics | |
Harris et al. | Ultimate limits in ultrasonic imaging resolution | |
CN104898123B (zh) | 基于角域虚拟源的水浸超声合成孔径聚焦成像方法 | |
US9384530B2 (en) | Enhanced ultrasound image formation using qualified regions of overlapping transmit beams | |
Nikolov et al. | 3D synthetic aperture imaging using a virtual source element in the elevation plane | |
Li et al. | Evaluational spatial compounding | |
Hoeks et al. | Methods to evaluate the sample volume of pulsed Doppler systems | |
CN113252793A (zh) | 基于掩膜压缩的单传感器三维超声成像装置及方法 | |
Lu | Producing bowtie limited diffraction beams with synthetic array experiment | |
Hein | Triple-beam lens transducers for three-dimensional ultrasonic fluid flow estimation | |
Jensen et al. | Equipment and methods for synthetic aperture anatomic and flow imaging | |
Plotnick et al. | Circular synthetic aperture sonar imaging of simple objects illuminated by an evanescent wavefield | |
Foster et al. | Breast imaging with a conical transducer/annular array hybrid scanner | |
Ophir et al. | Correlation artifacts in speed of sound estimation in scattering media | |
Wilkinson | Principles of real-time two-dimensional B-scan ultrasonic imaging | |
CN116115260B (zh) | 高分辨高对比快速计算的发收时序同步超声被动空化成像方法及系统 | |
Halliwell | Ultrasonic imaging in medical diagnosis | |
Lu et al. | High frame rate imaging with a small number of array elements |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20210813 |
|
WD01 | Invention patent application deemed withdrawn after publication |