CN104434094A - 一种磁-热-声耦合成像的电导率图像重建方法 - Google Patents
一种磁-热-声耦合成像的电导率图像重建方法 Download PDFInfo
- Publication number
- CN104434094A CN104434094A CN201410767707.5A CN201410767707A CN104434094A CN 104434094 A CN104434094 A CN 104434094A CN 201410767707 A CN201410767707 A CN 201410767707A CN 104434094 A CN104434094 A CN 104434094A
- Authority
- CN
- China
- Prior art keywords
- mrow
- conductivity
- msup
- sigma
- conductive object
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 20
- 238000003325 tomography Methods 0.000 title abstract 2
- 230000005284 excitation Effects 0.000 claims abstract description 41
- 238000003384 imaging method Methods 0.000 claims description 24
- 230000005684 electric field Effects 0.000 claims description 9
- 239000007822 coupling agent Substances 0.000 claims description 7
- 230000008878 coupling Effects 0.000 claims description 6
- 238000010168 coupling process Methods 0.000 claims description 6
- 238000005859 coupling reaction Methods 0.000 claims description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 5
- 230000003321 amplification Effects 0.000 claims description 3
- 239000013256 coordination polymer Substances 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 230000006698 induction Effects 0.000 claims description 3
- 238000009413 insulation Methods 0.000 claims description 3
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 3
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 238000003860 storage Methods 0.000 claims 1
- 238000004364 calculation method Methods 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 description 6
- 238000001514 detection method Methods 0.000 description 3
- 230000005672 electromagnetic field Effects 0.000 description 3
- 238000002604 ultrasonography Methods 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 239000004020 conductor Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/053—Measuring electrical impedance or conductance of a portion of the body
- A61B5/0536—Impedance imaging, e.g. by tomography
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biomedical Technology (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Radiology & Medical Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
一种磁-热-声耦合成像的电导率图像重建方法,首先根据接收到的电磁超声信号利用时间反演法获取导电物体的热声源分布,然后根据热声源分布再重建得到导电物体的电导率分布。具体过程为,首先利用时间反演法获取某一段层面上的热声源,然后利用插值法获取导电物体的热声源,接着根据激励电流,利用毕奥-萨伐尔定律计算得到一次磁矢位的空间分量A1,对导电物体进行空间离散,给出电导率的初值[σ]0,然后根据电流连续性定理,已知[σ]0和A1求解标量电位空间分量[φ]1,再接着将标量电位的空间分量[φ]1代入热声源与电导率满足的关系式,得到更新后的电导率[σ]1;利用[σ]1替换[σ]0,重复以上描述的过程,直到电导率相对误差满足ε=||([σ]1-[σ]0)/[σ]0||2≤ε0,即可停止迭代过程。
Description
技术领域
本发明涉及一种电导率图像的重建方法,特别涉及一种磁-热-声耦合成像的电导率图像重建方法。
背景技术
目前传统电阻抗成像技术的灵敏度和空间分辨率不高,主要因为电阻抗成像通常采用频率较低的电磁波作为激励。由于波长远远大于成像体,导致电磁场探测对比度高,但分辨率低。毋庸置疑,单一场都有其物理局限性。因此多物理场成像技术受到越来越多的关注,即将一种物理场作用于生物组织,转换为另一种物理场进行检测,由一种物理场提供分辨率,另一种物理场提供对比度,实现对比度和分辨率的同时提高。电磁场和超声相结合的多物理场成像技术正是考虑到电磁场对人体组织电导率的高对比度和超声波探测的高分辨率特性,成为人们的研究热点,磁热声成像作为一种新型的多物理场成像技术在最近一年受到重视。
磁热声成像是由新加坡南洋理工大学在2013年首次提出的新型的电阻抗成像方法,通过对导电物体施加低于20MHz的交变磁场,在导电物体内部产生感应电场,进而产生焦耳热,激发热弹性的声信号,检测声信号进行成像。该方法是一种以交变磁场作为激励源,基于生物组织内部焦耳热吸收率的差异,以超声作为信息载体的无损生物医学影像技术。与微波热声成像技术相比,激励源的频率降低,可以深入到导电体的更深处,使磁热声图像扩展到人体组织的深层。由测量的超声信号到电导率的重建分为两个过程,首先由测量的超声信号重建热声源分布,然后利用热声源分布重建电导率分布,目前的相关文献和专利只重建了热声源(S=σE2(σ),这里E为电场强度的空间分量),而没有提及电导率σ的重建。显然,电场强度E与电导率σ的分布有关,从热声源S中重建出电导率σ是非常困难的。
发明内容
本发明的目的是克服现有的磁热声成像方法无法给出电导率分布的不足,提出一种基于线性泊松方程直接迭代的磁-热-声耦合成像的电导率图像重建方法。本发明可以精确的重建导电物体的电导率。
本发明基于磁热声成像的原理:利用激励线圈对导电物体施加MHz电流激励,导电物体在电磁激励下产生焦耳热,进而产生超声信号,利用超声换能器接收超声信号,并对接收到的超声信号利用电导率图像重建算法进行图像重建,即可获取导电物体的电导率图像分布。
本发明磁-热-声耦合成像的电导率图像重建方法包括五个步骤:第一步首先获取导电物体的磁热声信号;第二步根据声压波动方程获取导电物体的热声源分布;第三步求解标量电位空间分量;第四步利用求解的标量电位空间分量计算新的电导率;第五步设定的电导率与求解的电导率进行相对误差计算,如果满足要求则新的电导率即为求解的电导率,否则把新的电导率设定为电导率初值,重新进行步骤三和步骤四,直到设定电导率与求解电导率之间的相对误差满足要求,则停止迭代,得到最终的导电物体的电导率。
具体步骤如下:
第一步:获取导电物体的磁热声信号
为了获取高分辨率和高信噪比的磁热声信号,本发明采用准连续脉冲激励方法对导电物体施加电磁激励。所述的准连续脉冲激励的脉冲宽度在2ns-0.1us可调,连续脉冲个数可调。电磁激励源产生的准连续脉冲激励信号通过激励线圈作用到导电物体上,激励线圈与导电物体通过空气隔离,导电物体在感应电流作用下产生焦耳热,进而产生电磁热声信号,电磁热声信号通过耦合剂耦合到超声换能器,超声换能器接收到电磁超声信号后通过数据处理采集电路进行放大、滤波等处理,并存储。
第二步:获取热声源分布
已知磁热声成像的声压波动方程:
其中r为热声源位置坐标,p(r,t)是声压,cs为热声源在介质中的传播声速,CP为导电物体的比热容,β为导电物体的热膨胀系数,S(r)为热声源分布,δ(t)是狄拉克函数,t是时间项;
选取导电物体的某一断层面z=z0,根据声压波动方程可以求解该断层面上的热声源分布S(x,y,z0),并求解导电物体上的所有断层面上的热声源分布S(x,y,z)。
热声源S同时是电导率和电场强度的函数,因此也表示为:
S=σE2=σE·E (2)
其中σ为导电物体的电导率,E为激励电流在导电物体内产生的电场强度;
第三步:求解标量电位空间分量
根据电流连续性定理
其中,A1为一次磁矢位的空间分量,一次磁矢位与导电物体无关,是激励源在真空中产生的磁矢位,φ是标量电位的空间分量,为哈密顿算符。
对导电物体进行空间离散,给出电导率的初值[σ]0,同时根据激励电流大小,利用毕奥-萨伐尔定律计算得到一次磁矢位的空间分量A1,将[σ]0和A1代入公式(3),结合电绝缘边界条件,进行线性有限元求解,求解得到标量电位空间分量[φ]1。
第四步:求解电导率σ
由公式(2)可以得出:
将热声源S(x,y,z)和标量电位的空间分量[φ]1代入公式(4),得到更新后的电导率[σ]1;
第五步:迭代求解电导率
根据第四步的求解的电导率,计算是否满足以下数学表达式:
ε=||([σ]1-[σ]0)/[σ]0||2≤ε0
其中ε0为给定的最大相对误差,ε为相对误差,[σ]0给定的电导率初值。
如果满足则电导率[σ]1即为所求解的电导率;如果不满足要求,利用[σ]1替换[σ]0,重复进行第三步和第四步,直到满足要求为止。
附图说明
图1本发明重建方法所涉及的磁热声信号获取示意图;
图中:1激励线圈,2水槽,3导电物体,4超声换能器。
具体实施方式
以下结合附图和具体实施方式进一步说明本发明。
本发明基于磁热声成像原理,所述的磁-热-声耦合成像的电导率图像重建方法包括五个步骤:
第一步首先获取导电物体的磁热声信号;第二步根据声压波动方程获取导电物体的热声源分布;第三步求解标量电位空间分量;第四步利用求解的标量电位空间分量计算新的电导率;第五步设定的电导率与求解的电导率进行相对误差计算,如果满足要求则新的电导率即为求解的电导率,否则把新的电导率设定为电导率初值,重新进行步骤三和步骤四,直到设定电导率与求解电导率之间的相对误差满足要求,则停止迭代,得到最终的导电物体的电导率。
具体步骤如下:
第一步:获取导电物体磁热声信号
如图1所示,为了获取高分辨率和高信噪比的磁热声信号,本发明采用准连续脉冲激励方法对导电物体施加电磁激励。所述的准连续脉冲激励的脉冲宽度在2ns-0.1us可调,连续脉冲个数可调。电磁激励源产生的准连续脉冲激励信号通过激励线圈1作用到导电物体3上,激励线圈1与导电物体3通过空气隔离,激励线圈1位于水槽2的耦合剂之外,导电物体3和超声换能器4位于水槽的耦合剂之内,同时保证导电物体位于激励线圈1产生的激励磁场之中。导电物体3在感应电流作用下产生焦耳热,进而产生电磁热声信号,电磁热声信号通过耦合剂耦合到超声换能器4,超声换能器4接收到电磁超声信号后通过数据处理采集电路进行放大、滤波等处理,并存储。
第二步:获取热声源分布
已知磁热声成像的声压波动方程:
其中r为热声源位置坐标,p(r,t)是声压,cs为热声源在介质中的传播声速,CP为导电物体的比热容,β为导电物体的热膨胀系数,S(r)为热声源分布,δ(t)是狄拉克函数,t是时间项;
选取导电物体的某一断层面z=z0,根据声压波动方程可以求解该断层面上的热声源分布S(x,y,z0),并求解导电物体上的所有断层面上的热声源分布S(x,y,z)。
热声源同时是电导率和电场强度的函数,因此也表示为:
S=σE2=σE·E (2)
其中σ为导电物体的电导率,E为激励电流在导电物体内产生的电场强度。
第三步:求解标量电位空间分量
根据电流连续性定理:
其中A1为一次磁矢位的空间分量,一次磁矢位与导电物体无关,是激励源在真空中产生的磁矢位,φ是标量电位的空间分量,为哈密顿算符。对导电物体进行空间离散,给出电导率的初值[σ]0,同时根据激励电流大小,利用毕奥-萨伐尔定律计算得到一次磁矢位的空间分量A1,将[σ]0和A1代入公式(3),结合电绝缘边界条件,进行线性有限元求解,求解得到标量电位空间分量[φ]1;
第四步:求解电导率
有公式(2)得出:
将热声源S(x,y,z)和标量电位的空间分量[φ]1代入公式(4),得到更新后的电导率[σ]1;
第五步:迭代求解电导率
根据第四步的求解的电导率,计算是否满足:
ε=||([σ]1-[σ]0)/[σ]0||2≤ε0
其中,ε0为给定的最大相对误差,ε相对误差,[σ]0给定的电导率初值。
如果满足则电导率[σ]1即为所求解的电导率;如果不满足要求,利用[σ]1替换[σ]0,重复进行第三步和第四步,直到满足要求为止。
Claims (1)
1.一种磁-热-声耦合成像的电导率图像重建方法,其特征在于:所述的磁-热-声耦合成像的电导率图像重建方法包括以下步骤:
第一步:获取导电物体磁热声信号
采用准连续脉冲激励方法对导电物体施加电磁激励;所述的准连续脉冲激励的脉冲宽度在2ns-0.1us可调,连续脉冲个数可调,电磁激励源产生的准连续脉冲激励信号通过激励线圈(1)作用到导电物体(3)上,激励线圈(1)与导电物体(3)通过空气隔离,也即,激励线圈(1)位于水槽(2)的耦合剂之外,导电物体(3)和超声换能器(4)位于水槽的耦合剂之内,导电物体(3)位于激励线圈(1)产生的激励磁场之中;导电物体(3)在感应电流作用下产生焦耳热,进而产生电磁热声信号,电磁热声信号通过耦合剂耦合到超声换能器(4),超声换能器(4)接收到电磁超声信号后通过数据处理采集电路进行放大、滤波处理,并存储;
第二步:获取热声源分布
已知磁热声成像的声压波动方程:
其中r为热声源位置坐标,p(r,t)是声压,cs为热声源在介质中的传播声速,CP为导电物体(3)的比热容,β为导电物体(3)的热膨胀系数,S(r)为热声源分布,δ(t)是狄拉克函数,t是时间项;
选取导电物体的某一断层面z=z0,根据声压波动方程求解该断层面上的热声源分布S(x,y,z0),并求解导电物体上的所有断层面上的热声源分布S(x,y,z);
热声源同时是电导率和电场强度的函数,因此也表示为:
S=σE2=σE·E (2)
其中,σ为导电物体的电导率,E为激励电流在导电物体内产生的电场强度;
第三步:求解标量电位空间分量
根据电流连续性定理:
其中A1为一次磁矢位的空间分量,一次磁矢位与导电物体无关,是激励源在真空中产生的磁矢位,φ是标量电位的空间分量,为哈密顿算符;
对导电物体进行空间离散,给出电导率的初值[σ]0,同时根据激励电流大小,利用毕奥-萨伐尔定律计算得到一次磁矢位的空间分量A1,将[σ]0和A1代入公式(3),结合电绝缘边界条件,进行线性有限元求解,求解得到标量电位空间分量[φ]1;
第四步:求解电导率
有公式(2)可以得出:
将热声源S(x,y,z)和标量电位的空间分量[φ]1代入公式(4),得到更新后的电导率[σ]1;
第五步:迭代求解电导率
根据第四步的求解的电导率,计算是否满足ε=||([σ]1-[σ]0)/[σ]0||2≤ε0
其中ε0为给定的最大相对误差,ε为相对误差,[σ]0给定的电导率初值;如果满足则电导率[σ]1即为所求解的电导率;如果不满足要求,利用[σ]1替换[σ]0,重复进行第三步和第四步,直到满足要求为止。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410767707.5A CN104434094B (zh) | 2014-12-14 | 2014-12-14 | 一种磁‑热‑声耦合成像的电导率图像重建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410767707.5A CN104434094B (zh) | 2014-12-14 | 2014-12-14 | 一种磁‑热‑声耦合成像的电导率图像重建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104434094A true CN104434094A (zh) | 2015-03-25 |
CN104434094B CN104434094B (zh) | 2017-04-19 |
Family
ID=52881403
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410767707.5A Active CN104434094B (zh) | 2014-12-14 | 2014-12-14 | 一种磁‑热‑声耦合成像的电导率图像重建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104434094B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106481339A (zh) * | 2016-09-14 | 2017-03-08 | 中国石油大学(华东) | 基于电磁‑热‑声效应的天然气水合物随钻探测与模拟方法 |
CN106910171A (zh) * | 2017-02-20 | 2017-06-30 | 中国科学院电工研究所 | 基于直接迭代的注入电流式热声电导率图像重建方法 |
CN107411742A (zh) * | 2017-07-31 | 2017-12-01 | 电子科技大学 | 热声成像图像的构建方法及图像处理装置 |
CN109839582A (zh) * | 2019-02-28 | 2019-06-04 | 中国计量大学 | 一种集成电路三维电流的磁成像测试方法及装置 |
CN111419185A (zh) * | 2020-04-08 | 2020-07-17 | 国网山西省电力公司电力科学研究院 | 一种声速不均匀的磁声成像图像重建方法 |
CN115372416A (zh) * | 2022-02-23 | 2022-11-22 | 重庆文理学院 | 一种感应式脉冲压缩磁声检测方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080154128A1 (en) * | 2005-05-27 | 2008-06-26 | Cardiospectra, Inc. | Hemoglobin contrast in magneto-motive optical doppler tomography, optical coherence tomography, and ultrasound imaging methods and apparatus |
CN101343999A (zh) * | 2008-09-03 | 2009-01-14 | 中国科学院电工研究所 | 阵列磁声电导率成像测井方法及装置 |
US8036732B2 (en) * | 2006-10-18 | 2011-10-11 | Board Of Regents, The University Of Texas System | Hemoglobin contrast in magneto-motive optical doppler tomography, optical coherence tomography, and ultrasound imaging methods and apparatus |
CN102894974A (zh) * | 2012-10-16 | 2013-01-30 | 中国科学院电工研究所 | 一种磁声电成像系统及成像方法 |
-
2014
- 2014-12-14 CN CN201410767707.5A patent/CN104434094B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080154128A1 (en) * | 2005-05-27 | 2008-06-26 | Cardiospectra, Inc. | Hemoglobin contrast in magneto-motive optical doppler tomography, optical coherence tomography, and ultrasound imaging methods and apparatus |
US8036732B2 (en) * | 2006-10-18 | 2011-10-11 | Board Of Regents, The University Of Texas System | Hemoglobin contrast in magneto-motive optical doppler tomography, optical coherence tomography, and ultrasound imaging methods and apparatus |
CN101343999A (zh) * | 2008-09-03 | 2009-01-14 | 中国科学院电工研究所 | 阵列磁声电导率成像测井方法及装置 |
CN102894974A (zh) * | 2012-10-16 | 2013-01-30 | 中国科学院电工研究所 | 一种磁声电成像系统及成像方法 |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106481339A (zh) * | 2016-09-14 | 2017-03-08 | 中国石油大学(华东) | 基于电磁‑热‑声效应的天然气水合物随钻探测与模拟方法 |
CN106481339B (zh) * | 2016-09-14 | 2019-04-19 | 中国石油大学(华东) | 基于电磁-热-声效应的天然气水合物随钻探测与模拟方法 |
CN106910171A (zh) * | 2017-02-20 | 2017-06-30 | 中国科学院电工研究所 | 基于直接迭代的注入电流式热声电导率图像重建方法 |
CN107411742A (zh) * | 2017-07-31 | 2017-12-01 | 电子科技大学 | 热声成像图像的构建方法及图像处理装置 |
CN107411742B (zh) * | 2017-07-31 | 2019-10-01 | 电子科技大学 | 热声成像图像的构建方法及图像处理装置 |
CN109839582A (zh) * | 2019-02-28 | 2019-06-04 | 中国计量大学 | 一种集成电路三维电流的磁成像测试方法及装置 |
CN109839582B (zh) * | 2019-02-28 | 2022-06-28 | 中国计量大学 | 一种集成电路三维电流的磁成像测试方法及装置 |
CN111419185A (zh) * | 2020-04-08 | 2020-07-17 | 国网山西省电力公司电力科学研究院 | 一种声速不均匀的磁声成像图像重建方法 |
CN115372416A (zh) * | 2022-02-23 | 2022-11-22 | 重庆文理学院 | 一种感应式脉冲压缩磁声检测方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN104434094B (zh) | 2017-04-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104434094B (zh) | 一种磁‑热‑声耦合成像的电导率图像重建方法 | |
CN104434099B (zh) | 利用线性双旋度方程进行磁热声成像的电阻率重建方法 | |
Mariappan et al. | Magnetoacoustic tomography with magnetic induction: bioimepedance reconstruction through vector source imaging | |
CN107064302B (zh) | 一种注入电流式热声成像电导率重建方法 | |
CN104146708B (zh) | 一种基于声换能器接收特性的磁声图像重建方法及装置 | |
CN104473640B (zh) | 一种磁热声成像的电导率重建方法 | |
CN104473639A (zh) | 一种基于最优化迭代算法的磁热声成像电阻率重建方法 | |
CN111948291B (zh) | 一种纳秒脉冲电场诱导热声成像的电导率重建方法 | |
Gao et al. | A correlated microwave-acoustic imaging method for early-stage cancer detection | |
CN104013388A (zh) | 基于低频连续波的磁声耦合成像激励与检测方法及装置 | |
CN104458818B (zh) | 一种基于线性泊松方程的磁热声成像电导率重建方法 | |
Perez-Liva et al. | Full-wave attenuation reconstruction in the time domain for ultrasound computed tomography | |
Zhang et al. | A reconstruction algorithm for thermoacoustic tomography with compensation for acoustic speed heterogeneity | |
Labyed et al. | TR-MUSIC inversion of the density and compressibility contrasts of point scatterers | |
CN106885842A (zh) | 一种注入电流式热声成像电阻率重建方法 | |
CN104434100A (zh) | 一种磁热声成像的电阻率重建方法 | |
CN106580249B (zh) | 一种注入电流式热声成像方法 | |
CN107049315B (zh) | 基于最优化迭代方法的注入电流式热声电阻率图像重建方法 | |
CN103142216B (zh) | 一种基于光声成像技术的多层介质声速计算的方法 | |
CN115219594A (zh) | 一种高电导率非铁磁材料缺陷的磁声电检测方法与装置 | |
CN103605868A (zh) | 基于介质有限元叠加的磁声耦合成像声信号求解方法 | |
Zhang et al. | Lamb waves topological imaging combining with Green’s function retrieval theory to detect near filed defects in isotropic plates | |
CN106910171B (zh) | 基于直接迭代的注入电流式热声电导率图像重建方法 | |
Ali et al. | 2-D Slicewise Waveform Inversion of Sound Speed and Acoustic Attenuation for Ring Array Ultrasound Tomography Based on a Block LU Solver | |
Ramirez et al. | Sparsity constrained Born inversion for breast cancer detection |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |