CN108366771B - 电磁干涉图案识别层析成像 - Google Patents

电磁干涉图案识别层析成像 Download PDF

Info

Publication number
CN108366771B
CN108366771B CN201680074323.2A CN201680074323A CN108366771B CN 108366771 B CN108366771 B CN 108366771B CN 201680074323 A CN201680074323 A CN 201680074323A CN 108366771 B CN108366771 B CN 108366771B
Authority
CN
China
Prior art keywords
image
electromagnetic interference
electromagnetic
forming
disturbed
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.)
Active
Application number
CN201680074323.2A
Other languages
English (en)
Other versions
CN108366771A (zh
Inventor
赛尔古依·Y·谢苗诺夫
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
EMTensor GmbH
Original Assignee
EMTensor GmbH
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by EMTensor GmbH filed Critical EMTensor GmbH
Publication of CN108366771A publication Critical patent/CN108366771A/zh
Application granted granted Critical
Publication of CN108366771B publication Critical patent/CN108366771B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
    • G01R29/08—Measuring electromagnetic field characteristics
    • G01R29/0864—Measuring electromagnetic field characteristics characterised by constructional or functional features
    • G01R29/0892—Details related to signal analysis or treatment; presenting results, e.g. displays; measuring specific signal features other than field strength, e.g. polarisation, field modes, phase, envelope, maximum value
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
    • A61B34/20—Surgical navigation systems; Devices for tracking or guiding surgical instruments, e.g. for frameless stereotaxis
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
    • A61B5/004—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part
    • A61B5/0042—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part for the brain
    • 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
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
    • A61B34/20—Surgical navigation systems; Devices for tracking or guiding surgical instruments, e.g. for frameless stereotaxis
    • A61B2034/2046—Tracking techniques
    • A61B2034/2051—Electromagnetic tracking systems
    • 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/0507—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  using microwaves or terahertz waves

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Surgery (AREA)
  • Medical Informatics (AREA)
  • Radiology & Medical Imaging (AREA)
  • General Health & Medical Sciences (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biomedical Technology (AREA)
  • Molecular Biology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Pathology (AREA)
  • Biophysics (AREA)
  • Neurology (AREA)
  • Electromagnetism (AREA)
  • Robotics (AREA)
  • General Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Optics & Photonics (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)

Abstract

用于图像重建***中的电磁干扰图案识别层析成像(EMIPRT)方法包括：经由电磁层析成像***生成与成像域中的对象相对应的电磁场数据，并且以递归的方式重复地使用所生成的电磁场数据，进行以下处理：形成未受干扰的电磁干涉图像，基于未受干扰的电磁干涉图像形成受干扰的电磁界面图像，识别重复形成的受干扰的电磁界面图像中的电磁干涉图案，以及通过使所识别的电磁干涉图案从受干扰的电磁界面图像中无效或衰减来形成叠加图像。形成受干扰的电磁界面图像还基于目标因子，该目标因子是在形成未受干扰的电磁干涉图像的步骤期间的实验电磁场与计算出的电磁场之间的差异的泛函。在形成叠加图像的每个重复步骤之后，该方法还包括确定是否已经达到收敛目标。

Description

电磁干涉图案识别层析成像
相关申请的交叉引用
就美国而言，本申请是于2015年10月16日提交的标题为“电磁干涉图案识别层析成像(ELECTROMAGNETIC INTERFERENCE PATTERN RECOGNITION TOMOGRAPHY)”的美国临时专利申请序列号62/242,915的美国非临时专利申请，并且根据35U.S.C.§119(e)要求其优先权，该‘915申请的全部内容通过引用明确合并于此。此外，以下专利、专利申请和专利申请公开的每一个的全部内容通过引用合并于此：
(a)授予Semenov的美国专利No.9,414,749，其于2016年8月16日发布，并且在之前于2014年6月5日作为美国专利申请公开No.2014/0155740 A1公开，其至少旨在提供关于当前专利申请的发明的***和环境的背景和技术信息；和
(b)于2012年1月12日公开的美国专利申请公开No.2012/0010493 A1，其基于授予Semenov的2011年6月30日提交的名称为“电磁层析成像(EMT)差分(动态)融合成像***和方法(SYSTEMS AND METHODS OF ELECTROMAGNETIC TOMOGRAPHY(EMT)DIFFERENTIAL(DYNAMIC)FUSED IMAGING)”美国专利申请序列号13/173,078，其旨在提供有关4D EMT成像的背景和技术信息。
版权声明
本专利文件中的所有材料受美国和其他国家版权法的保护。版权所有人不反对任何人对专利文件或专利公开进行传真复制，因为它出现在官方的政府记录中，但是，所有其他版权权利都被保留。
技术领域
本发明通常涉及电磁层析成像，并且具体地涉及使用电磁干涉图案识别策略来消除由重复干涉图案引起的放大失真。
背景技术
电磁层析成像(EMT)是一种相对较新的成像形式，具有生物医学和工业应用的巨大潜力。生物医学应用包括但不限于生物组织的功能和病理状况的非侵入性评估。工业应用包括但不限于石油和天然气勘探、矿山搜索和评估、以及非金属管道内的流量评估。使用EMT，可以区分生物组织等对象，并相应地根据这些对象的介电特性差异进行成像。基于最近证明的组织介电特性对组织的各种功能和病理状况(例如血液和氧气含量、缺血和梗塞恶性肿瘤、水肿等)的依赖性，EMT被认为具有生物医学应用的高潜力。
在过去的十年或更长的时间里，二维(2D)、三维(3D)乃至“四维”(4D)EMT***和图像重建方法已经得到发展。已经证明了用于各种生物医学和工业应用的技术的可行性，例如，用于心脏成像和四肢成像。
如在任何层析成像中那样，典型的EMT成像场景由复杂信号的测量周期组成，该复杂信号受到位于所谓的成像域内的研究对象的存在的影响、由位于所述对象周围的不同点处的多个发射机产生、并且在位于所述对象周围的不同点处的多个接收机上测量。这在图1中示出。发射机和接收机的位置可以位于成像域内、成像域的边界上或成像域外。如本文其他地方所述，然后可以在图像重建方法中使用所测量的EM信号矩阵，以便重建对象19的介电特性的3D分布，即构建对象的3D图像。在前述美国专利No.9,414,749中公开和描述的硬件组件可用于生成必要的EM信号以重建对象19的3D图像。
一般而言，精确地描述具有成像域21的EM场的分布对于图像重建是非常重要的。即使在内部没有感兴趣的对象的情况下，成像室内的EM场的分布也是非常复杂的现象。
本发明涉及使用电磁EM(场)对对象19的结构进行成像。对象19可以是人体或人体的一部分，例如头部、躯干、手臂等，但也可以是没有金属屏蔽的任何对象。使用EM场对被强屏蔽(但不是金属屏蔽)的对象的内部进行成像是一个高度复杂的问题。这种应用的一个示例是对人脑的成像。然而，应该理解的是，其他这样的应用可以包括对被骨结构屏蔽的任何人体组织的成像。对被骨骼所包围的大脑或其他组织进行EM成像呈现出非常复杂的高介电对比度问题。挑战在于重建被高介电对比度屏障有效屏蔽的深部脑组织的隐藏特性，该高介电对比度屏障包括颅骨(介电特性在16+j5范围内)和脑脊髓液(介电特性在60+j30范围内)。尽管如上所述，本发明适用于对任何对象成像，但据信其特别适用于对被强屏蔽对象的内部成像。
高介电对比度对象(包括生物对象)的EMT成像具有非常复杂的所谓“衍射层析成像”问题。(Devaney A.J.“Current research topics in diffraction tomography”,inInverse Problmes in Scattering and Imaging,M.BNertero and E.R.Pike,Eds,NewYork:Adam Hilger,1992,pp.47-58.)高水含量的组织(例如但不限于肌肉组织)和低水含量的组织(例如但不限于骨骼)之间的高介电对比度在使用EM场进行成像时呈现出额外的复杂性。为了尝试解决衍射层析成像问题，已经开发出了使用标量近似和矢量近似的2D和3D几何中的各种方法。(参见H.Harada,D.Wall,T.Takenaka,and T.Tanaka,“Conjugategradient method applied to inverse scattering problem”,IEEE Trans.Antennasand Propagations,vol.43,784-792,Aug.1995；R.E.Kleinman,and P.M.van den Berg,”Amodified gradient method for two-dimensional problems in tomography”,J.Comput.Appl.Math.,vol.42,pp.17-35,Jan.1992；A.Abubakar,P.M.van den Berg,andJ.J.Mallorqui,“Imaging of Biomedical Data Using A Multiplicative RegularizedSource Inversion Method,”IEEE Trans.Microwave Theory and Techniques,v.50,pp.1761-1771,Jul.2002；N.Joachimowicz,J.J.Mallorqui,J.Ch.Bolomey,andA.Brouguetas,“Convergence and stability assessment of Newton-Kantorovichreconstruction algorithms for microwave tomography”,IEEE Trans.MedicalImaging,vol.17,pp.562-570,Aug.1998；P.Lobel,R.Kleinman,Ch.Pichot,L.Blanc-Fйraud,and M.Barlaud“Conjugate Gradient Method for Solving Inverse Scatteringwith Experimental Data“,IEEE Antennas&Propagation Magazine,Vol.38,pp.48--51,Jun.1996；W.C.Chew,and Y.M.Wang,“Reconstruction of two-dimensionalpermittivity distribution using the distorted Born iterative method”,IEEETrans.Medical Imaging,vol.9,p.218-225,June 1990；P.M.Meaney,K.D.Paulsen,A.Hartov,and R.K.Crane,“Microwave imaging for tissue assessment:Initialevaluation in multitarget tissue equivalent phantoms”,IEEE Trans.BiomedicalEngineering,vol.43,pp.878-890,Sept.1996.)
在先前的作品中，已经开发出了数学算法及其相应的***和软件实现，其被证明是非常可靠的、并且已经实现了对不同尺寸(从离体的犬心脏的几厘米直到2D、3D和3D矢量案例中的全尺寸身体)的对象成像。(参见See A.E.Souvorov,A.E.Bulyshev,S.Y.Semenov,R.H.Svenson,A.G.Nazarov,Y.E.Sizov,and G.P.Tatsis,“Microwave tomography:A two-dimensional Newton iterative scheme”,IEEE Trans.Microwave Theory andTechniques,vol.46,pp.1654-1659,Nov.1998；A.E.Bulyshev,A.E.Souvorov,S.Y.Semenov,R.H.Svenson,A.G.Nazarov,Y.E.Sizov,and G.P.Tatsis,“Three-dimensional microwave tomography.Theory and computer experiments in scalarapproximation”,Inverse Problems,vol.16,pp.863-875,Jun.2000；Bulyshev A.E,Souvorov A.E.,Semenov S.Y.,Posukh V.G.,Sizov Y.E.“Three-dimensional VectorMicrowave Tomography.Theory and Computational experiments”,Inverse Problems,2004,20,4,1239-1259；Semenov S.Y.,Bulyshev A.E.,Souvorov A.E.,Svenson R.H.,Sizov Y.E.,Borisov V.Y.,Posukh V.G.,Kozlov I.M.,Nazarov A.G.,Tatsis G.P."Microwave Tomography:Theoretical and Experimental Investigation of theIteration Reconstruction Algorithm",IEEE Trans MTT,1998,46,2,133-141；SemenovS.Y.,Bulyshev A.E.,Abubakar A.,Posukh V.G.,Sizov Y.E.,Souvorov A.E.,Van denBerg P.,Williams T.“Microwave tomographic imaging of the high dielectriccontrast objects using different imaging approaches”,IEEE Trans.MTT,v.53,No7,pp 2284-2294,2005；Semenov S.Y.,Kellam J.F.,Althausen P.,Williams T.C.,Abubakar A.,Bulyshev A.,Sizov Y.“Microwave tomography for functional imagingof extremity soft tissues.Feasibility assessment”,Phys.Med.Biol.,2007,52,5705-5719；Semenov S.Y.,Bulyshev A.E.,Posukh V.G.,Sizov Y.E.,Williams T.C.,Souvorov A.E.“Microwave tomography for detection/imaging of myocardialinfarction.1.Excised canine hearts”,Annals of Biomedical Engineering,2003,31,262-270.)然而，这些先前的作品中所引用的上述方法都没有被证明在对被强屏蔽对象的内部成像时有效。需要新的方法对具有高介电对比度屏障的对象(例如但不限于人脑)的EMT成像进行精确表示。
发明内容
本发明的一些示例性实施例可以克服一个或多个上述缺点以及以上未描述的其他缺点，但是本发明不需要克服上述任何特定的缺点，并且本发明的一些示例性实施例可以不克服上述任何缺点。
广义地说，根据一个方面的本发明是以以下方式使用电磁场的方法：即在成像域内生成电磁干涉图片，显示对象的3D介电结构与电磁干涉图案的叠加，同时电磁干涉图案被进一步识别并被应用于3D电磁叠加，使电磁干涉图案无效或衰减并显示对象的3D介电结构。
广义地说，根据另一方面的本发明是一种4D动态融合电磁图案识别层析成像的方法。
广义地说，根据另一方面的本发明是利用4D动态融合电磁识别层析成像来监测生物组织的生存能力和功能状况的方法。
广义地说，根据另一方面的本发明是用于在所示和所述的图像重建***中使用的电磁干涉图案识别层析成像(EMIPRT)方法。
广义地说，根据另一方面的本发明是在图像重建***中使用的电磁干涉图案识别层析成像(EMIPRT)方法，其包括：经由电磁层析成像***，生成与成像域中的对象相对应的电磁场数据，其中所述电磁场数据在多个发射机处产生并在与所述对象相互作用之后在多个接收机处被测量；并且以递归的方式重复地使用所生成的电磁场数据：形成未受干扰的电磁干涉图像，至少部分地基于所述未受干扰的电磁干涉图像来形成受干扰的电磁干涉图像，识别重复形成的受干扰的电磁干涉图像中的电磁干涉图案，以及通过使所识别的电磁干涉图案从受干扰的电磁干涉图像中无效或衰减来形成叠加图像。
在该方面的特征中，所述至少部分地基于所述未受干扰的电磁干涉图像来形成受干扰的电磁干涉图像的步骤包括至少部分地基于对目标因子的确定来形成受干扰的电磁干涉图像，所述目标因子是在形成未受干扰的电磁干涉图像的步骤期间的实验电磁场与计算出的电磁场之间的差异的泛函。在进一步的特征中，所述目标因子被确定为
Figure GDA0001822578180000061
其中，当发射机i是所述电磁场的源时，
Figure GDA0001822578180000062
或
Figure GDA0001822578180000063
分别是由接收机j测得的所述电磁场的z分量的实验仿真值或测得值，其中
Figure GDA0001822578180000064
通常以
Figure GDA0001822578180000065
的形式表示，其中α、β和γ是系数，所述系数为非零实数或零值，其中Ω是正则化算子，并且其中
Figure GDA0001822578180000066
是其自变量的函数；所述目标因子被确定为
Figure GDA0001822578180000067
其中，当发射机i是所述电磁场的源时，
Figure GDA0001822578180000068
或
Figure GDA0001822578180000069
分别是由接收机j测得的所述电磁场的z分量的实验仿真值或测得值，并且其中
Figure GDA00018225781800000610
是所述电磁场的实验测得的z分量的最大范数；所述目标因子被确定为
Figure GDA00018225781800000611
其中，当发射机i是所述电磁场的源时，
Figure GDA00018225781800000612
或
Figure GDA00018225781800000613
分别是由接收机j测得的所述电磁场的z分量的实验仿真值或测得值，并且其中当发射机i是幂为θ的所述电磁场的源时，
Figure GDA00018225781800000614
是由接收机j测得的所述电磁场的实验测得的z分量的范数；所述形成受干扰的电磁干涉图像的步骤包括计算
Figure GDA00018225781800000615
其中
Figure GDA0001822578180000071
和
Figure GDA0001822578180000072
是来自于所述电磁层析成像***中的位于物理源(从1到N)的位置处和相应的物理接收机(从1到M)的位置处的频率为fk的电磁源的3D电磁场(x，y，z)分布，作为共轭值，并且其中，Object Factori，j是从发射机i到接收机j的所述目标因子的第“ij”个分量；所述识别重复形成的受干扰的电磁干涉图像中的电磁干涉图案的步骤包括计算总和
Figure GDA0001822578180000073
其中，
Figure GDA0001822578180000074
和
Figure GDA0001822578180000075
是来自于电磁层析成像***中的位于物理源(从1到N)的位置处和相应的物理接收机(从1到M)的位置处的频率为fk的电磁源的3D电磁场(x，y，z)分布，作为共轭值，并且其中，Object Factori，j是从发射机i到接收机j的所述目标因子的第“ij”个分量；和/或其中所述识别重复形成的受干扰的电磁干涉图像中的电磁干涉图案的步骤包括针对迭代i>1，计算
Figure GDA0001822578180000076
Figure GDA0001822578180000077
其中，为了简便省略了频率项，其中，
Figure GDA0001822578180000078
和
Figure GDA0001822578180000079
是来自于电磁层析成像***中的位于物理源(从1到N)的位置处和相应的物理接收机(从1到M)的位置处的电磁源的3D电磁场(x，y，z)分布，作为共轭值，并且其中，Object Factori，j是从发射机i到接收机j的所述目标因子的第“ij”个分量。
在该方面的另一个特征中，该方法还包括在形成叠加图像的每个重复步骤之后执行的确定是否已经达到收敛目标的步骤。
在该方面的另一个特征中，该方法用作生成4D差分(动态)融合图像的方法的一部分。在进一步的特征中，生成4D差分(动态)融合图像包括将指示相对生理变化的至少一个连续形成的图像与基线解剖图像相组合以显示为单个统一图像；和/或生成4D差分(动态)融合图像的方法被用作利用4D动态融合电磁图案识别层析成像来监测生物组织的生存能力和/或功能状况的方法的一部分。
在该方面的另一特征中，按顺序地执行以下步骤：形成未受干扰的电磁干涉图像，至少部分地基于未受干扰的电磁干涉图像形成受干扰的电磁干涉图像，识别重复形成的受干扰的电磁干涉图像中的电磁干涉图案，以及通过使所识别的电磁干涉图案从所述受干扰的电磁干涉图像中无效或衰减来形成叠加图像。
在该方面的另一特征中，该方法还包括经由显示单元显示叠加图像的步骤。
广义地说，根据另一方面的本发明是一种使用电磁干涉图案识别层析成像来重建图像的方法，包括：经由电磁层析成像***，生成与成像域中的对象相对应的电磁场数据，其中，所述电磁场数据在多个发射机处产生并与所述对象相互作用之后，在多个接收机处测量所述电磁场数据；并且以递归的方式重复地使用所生成的电磁场数据：形成未受干扰的电磁干涉图像，至少部分地基于未受干扰的电磁干涉图像来形成受干扰的电磁干涉图像，识别重复形成的受干扰的电磁干涉图像中的电磁干涉图案，以及通过使所识别的电磁干涉图案从受到干扰的电磁干涉图像中无效或衰减来形成叠加图像。
在该方面的特征中，所述至少部分地基于所述未受干扰的电磁干涉图像来形成受干扰的电磁干涉图像的步骤包括：至少部分地基于对目标因子的确定来形成受干扰的电磁干涉图像，所述目标因子是在形成未受干扰的电磁干涉图像的步骤期间的实验电磁场与计算出的电磁场之间的差异的泛函。
在该方面的另一个特征中，该方法还包括在形成叠加图像的每个重复步骤之后执行的确定是否已经达到收敛目标的步骤。
在该方面的另一个特征中，该方法用作生成4D差分(动态)融合图像的方法的一部分。在进一步的特征中，生成4D差分(动态)融合图像包括将指示相对生理变化的至少一个连续形成的图像与基线解剖图像相组合以显示为单个统一图像；和/或生成4D差分(动态)融合图像的方法被用作利用4D动态融合电磁图案识别层析成像监测生物组织的生存能力和/或功能状况的方法的一部分。
在该方面的另一特征中，按顺序地执行以下步骤：形成未受干扰的电磁干涉图像，至少部分地基于未受干扰的电磁干涉图像形成受干扰的电磁干涉图像，识别重复形成的受干扰的电磁干涉图像中的电磁干涉图案，以及通过使所识别的电磁干涉图案从所述受干扰的电磁干涉图像中无效或衰减来形成叠加图像。
在该方面的另一特征中，该方法还包括经由显示单元显示叠加图像的步骤。
广义地说，根据另一方面的本发明是一种使用电磁干涉图案识别层析成像的图像重建***，包括：电磁层析成像***，其生成与成像域中的对象相对应的电磁场数据，所述电磁层析成像***具有多个电磁发射机、测量在所述多个发射机处产生并与所述对象相互作用之后的电磁场数据的多个接收机、以及边界设备；以及处理中心，其以递归的方式重复地使用所生成的电磁场数据，执行以下步骤：形成未受干扰的电磁干涉图像，至少部分地基于未受干扰的电磁干涉图像来形成受干扰的电磁干涉图像，识别重复形成的受干扰的电磁干涉图像中的电磁干涉图案，以及通过使所识别的电磁干涉图案从所述受干扰的电磁干涉图像中无效或衰减来形成叠加图像。
在该方面的特征中，所述至少部分地基于未受干扰的电磁干涉图像来形成受干扰的电磁干涉图像的步骤包括：至少部分地基于对目标因子的确定来形成受干扰的电磁干涉图像，所述目标因子是在形成未受干扰的电磁干涉图像的步骤期间的实验电磁场与计算出的电磁场之间的差异的泛函。
在该方面的另一个特征中，所述处理中心还在形成叠加图像的每个重复步骤之后执行确定是否已经达到收敛目标的步骤。
在该方面的另一个特征中，由所述处理中心执行的步骤用作生成4D差分(动态)融合图像的方法的一部分。在进一步的特征中，生成4D差分(动态)融合图像包括将指示相对生理变化的至少一个连续形成的图像与基线解剖图像相组合以显示为单个统一图像；和/
或生成4D差分(动态)融合图像的方法被用作利用4D动态融合电磁图案识别层析成像监测生物组织的生存能力和/或功能状况的方法的一部分。
在该方面的另一特征中，按顺序地执行以下步骤：形成未受干扰的电磁干涉图像，至少部分地基于未受干扰的电磁干涉图像形成受干扰的电磁干涉图像，识别重复形成的受干扰的电磁干涉图像中的电磁干涉图案，以及通过使所识别的电磁干涉图案从所述受干扰的电磁干涉图像中无效或衰减来形成叠加图像。
在该方面的另一特征中，该方法还包括显示单元，其显示叠加图像。
根据下文提供的详细描述，本发明的其他应用领域将变得显而易见。应该理解的是，具体实施方式和具体示例虽然指示了本发明的优选实施例，但仅用于说明的目的，而不旨在限制本发明的范围。
附图说明
通过以下结合附图的详细描述，本发明的其他特征、实施例和优点将变得显而易见，其中：
图1是电磁层析成像(EMT)***的部分的简化示意图；
图2是根据本发明的一个或多个优选实施例在图像重建***中使用的EM干涉图案识别层析成像方法的流程图；
图3是用于针对人的头部成像设计的成像域的示例性未受干扰的EM干涉图像；
图4A和图4B是用于两个人的头部案例的示例性受干扰的EM干涉图像；
图5是用于图2中所示的图案识别框的一个图案识别策略和一组公式的流程图。
图7是可以在针对多个子域中的每一个的叠加图像形成框中执行的一种可能的图案校正策略的流程图；
图8A是来自于类似于图2的方法的迭代成像方法的虚拟的中风患者的头部的示例性输出图像，但是在图像形成期间未使用图案识别框及其输出信息；以及
图8B是与图8A相同的虚拟中风患者的头部的示例性输出图像，但是在叠加图像形成框期间使用了图案识别框及其输出信息。
具体实施方式
作为初步事项，相关领域中具有普通技术的人员(“普通技术人员”)容易理解，本发明具有广泛的用途和应用。此外，被讨论并被确定为“优选”的任何实施例被认为是用于执行本发明的预期的最佳模式的一部分。在提供完整且可实现的本发明的公开时，还可以讨论其他实施例以用于另外的说明目的。应该理解，任何实施例可以仅包含本发明的一个或多个上述公开的方面，并且可以进一步包含一个或多个上述公开的特征。此外，许多实施例(例如改编、变化、修改和等同布置)将由本文描述的实施例隐含地公开并且落入本发明的范围内。
因此，尽管本文中关于一个或多个实施例详细描述了本发明，但应理解的是，本公开是本发明的说明性和示例性，并且仅为了提供完整且可实现的本发明的公开目的。本文的一个或多个实施例的详细公开不旨在或不被解释为限制本发明所提供的专利保护范围，该范围由权利要求及其等同物限定。并不旨在通过将任何权利要求读成在权利要求本身没有明确出现的本文中找到的限制解释为限定本发明提供的专利保护范围。
因此，例如，本文描述的各种过程或方法的步骤的任何序列顺序和/或时间顺序是说明性的而非限制性的。因此，应该理解，尽管各种过程或方法的步骤可以被示出和描述为按照序列顺序或时间顺序，但是任何这样的过程或方法的步骤在没有指示的情况下不限于以任何特定序列或顺序来执行。实际上，这样的过程或方法中的步骤通常可以以各种不同的序列和顺序执行，而仍落入本发明的范围内。因此，本发明提供的专利保护范围旨在由所附权利要求限定，而不是由此处的说明书阐述。
另外，重要的是要注意，本文中使用的每个术语是指普通技术人员基于该术语在本文中的的上下文使用将会理解该术语的含义。当本文中所使用的术语的含义(正如由普通技术人员基于该术语的上下文使用所理解的)以任何方式与该术语的任何特定字典定义不同不同时，该术语的范围旨在以普通技术人员应理解为准。
关于35 U.S.C.§112,
Figure GDA0001822578180000121
的适用性，除非在权利要求要素中实际使用了“用于...的方法”或“用于...的步骤”这样的明确短语，否则不旨在根据该法定条款来阅读这种权利要求要素，因此，该法定条款旨在应用于对这种权利要求元素的解释。
此外，重要的是要注意，如本文所使用的，除非上下文使用另有规定，否则“一”和“一个”通常都表示“至少一个”，但不排除多个。因此，提到“具有苹果的野餐篮子”描述了“具有至少一个苹果的野餐篮子”以及“具有复数个苹果的野餐篮子”。相反，提到“具有单个苹果的野餐篮子”描述了“仅有一个苹果的野餐篮子。
当在本文中用于结合一系列项目时，“或”表示“至少一个项目”，但不排除所述系列中的多个项目。因此，提到“具有奶酪或饼干的野餐篮子”描述了“有奶酪而没有饼干的野餐篮子”、“有饼干而没有奶酪的野餐篮子”和“具有奶酪和饼干两者的野餐篮子”。最后，当在本文中用于结合一系列项目时，“和”表示“列表中的所有项目”。因此，提及“具有奶酪和饼干的野餐篮子”描述了“具有奶酪的野餐篮子，其中野餐篮子还具有饼干”以及描述了“具有饼干的野餐篮子，其中野餐篮子还有奶酪”。
现在参照附图，其中在几个示图中相同的附图标记表示相同的组件，接下来描述本发明的一个或多个优选实施例。以下对一个或多个优选实施例的描述本质上仅仅是示例性的，决不旨在限制本发明、其应用或用途。
图1是电磁层析成像(EMT)***10的部分的简化示意图。在***10中，边界设备12围绕成像域21。对象19(在该示例中为人的头部)放置在成像域21中。多个EM硬件装置30(通常但不一定总是设置在边界设备12上)用作发射机(源)和/或接收机。(如在本文中所述，硬件装置30通常被认为用作发射机和接收机两者，但是将理解，可以替代地或另外地利用彼此分离的发射机和接收机)。计算机***28用作处理中心，其中初始数据被预处理并且图像被重构和后处理。在至少一些实施例中，EMT***10还包括背景或匹配液或其他介质。当对象19成像时，匹配液是成像域内需要或有用的溶液或凝胶，以解决电磁体匹配问题和/或其他问题。在至少一些实施例中，计算机***28及其数据处理功能和成像软件直接连接到EM场发送/接收硬件装置30，而在其他实施例中，计算机***28中的一些或全部通过无线技术和/或高速有线连接而远程连接。在功能上，EMT***10的大部分操作可以类似于前述美国专利No.9,414,749中描述的操作，但是可以在本文中描述各种特定实施例和特征。
如上所述，高介电对比度对象(包括生物对象)的EMT成像涉及所谓的“衍射层析成像”的非常复杂的问题。在使用EM场进行成像时，高含水量的组织(例如但不限于肌肉组织)和低含水量的组织(例如但不限于骨骼)之间的高介电对比度呈现额外的复杂性。图2是根据本发明的一个或多个优选实施例的用于图像重建***的EM干涉图案识别层析成像(EMIPRT)方法100的流程图。可以利用使用图1中简要描述的EMT***10的一些或全部元件而生成的EM场数据来执行EMIPRT方法100。
EMIPRT方法100是迭代过程，其中在每次迭代之后通过各种图像形成处理进行收敛检查(在步骤180处示出)，直到获得合适的结果并将其提供为图像输出200。EMIPRT方法100的主要输入是成像域数据110和实验测得的EM数据120。成像域数据110包括与成像域21有关的位于多个空间位置处和第k个频率(k从1到K)(Ei＝1 to N(fk＝1 to K))处的源EM场硬件装置30的空间位置、和位于多个空间位置和第k个频率(k从1到K)(Ej＝1 to M(fk＝1 to K))处的接收机EM场硬件装置30的空间位置，其中N是源EM场装置30的数量，并且M是接收机EM场装置30的数量。成像域数据110还包括成像域21内的匹配介质(ε0)的介电特性，并且还可以包括其他的环境或物理信息/数据。实验测得的EM数据120是包含在接收机的点上通过实验测得的EM场(例如以振幅、相位和极化形式或以实部、虚部和极化形式)的矩阵，其表示为
Figure GDA0001822578180000141
EM场接收机的点可以与任何或全部的EM场装置30相关联。
在框130处，形成“未受干扰的”EM干涉图像。在初始处理中，仅基于EMIPRT方法100的成像域数据110来获得未受干扰的EM干涉图像(作为未受干扰的图像形成框130的输出135示出)的形成。在随后的迭代中，形成未受干扰的干涉图像优选地还使用成像域21内的介电特性的空间分布。根据应用，空间分布可以是2D空间分布或3D空间分布。在至少一些实施例中，介电特性的3D空间分布在笛卡尔坐标系中被描述为ε*(x，y，z)，并且在下面描述的重叠图像形成框170中在先前的迭代期间被获得。介电特性的3D空间分布ε*(x，y，z)在介电尺度中可以被表征为研究对象19的图像。
未受干扰的EM干涉图像形成框130的任务是合成来自于位于成像域21的内部或边界上的多个EM源/接收机30的2D或3D干涉图像(图片)135。图像形成过程期间，所有EM场装置30(i＝1到N和j＝1到M)被认为是EM场的源。由以下函数生成对(3D)未受干扰的EM干涉图像135的确定：
Figure GDA0001822578180000142
其中Ei＝1 to N(fk，x，y，z)和Ej＝1 to M(fk，x，y，z)是来自于位于物理源(从1到N)位置和相应的物理接收机(从1到M)位置的探测频率为fk的EM场装置的成像域21内的3DEM场(x，y，z)分布，并且Wk是频率权重函数，它考虑了将频率相关的EM干涉图像的不同输入输入到函数(1)中。图3是用于针对人的头部成像设计的成像域21的示例性未受干扰的EM干涉图像135。图像135是以特定的频率获得的，在这种情况下该频率是1GHz。图像是X-Y平面中的2D横截面视图，其中X和Y刻度以厘米为单位。
回到图2，EMIPRT方法100中的下一个主要操作是确定应用至EM干涉图像的目标因子，如框145所示。目标因子优选地是在特定探测频率处针对每对发送装置(i＝1到N)30到接收装置(j＝1到M)30而测得的EM场(来自框120)和仿真的EM场(来自框130)之间的差异的泛函。因此，目标因子确定框145的主输入是方法100的两个主要输入中的第二个(先前描述的在接收机的点上通过实验测得的EM场120)以及从未受干扰的EM干涉图像形成框130生成的输入数据140(在接收机的点上计算的或仿真的EM场)。输入到目标因子确定框145的这两个主要输入是数据的矩阵。第一个矩阵包含通过实验测得的EM场，表示为
Figure GDA0001822578180000151
第二个矩阵包含仿真的EM场，表示为
Figure GDA0001822578180000152
其中在这两种情况下，
Figure GDA0001822578180000153
是表示EM场的电(E)分量和/或磁(H)分量的复矢量值。这些值可以表示为E场和/或H场中的每个矢量分量的复数(实部和虚部)值、或表示为E场和/或H场的标量值。它们也可以表示为E场和/或H场中的每个矢量分量的幅度和/或相位、或者表示为E场和/或H场的标量值。在本发明的至少一个实施例中，目标因子可以以下面的一般术语来表示：
Figure GDA0001822578180000154
其中，当发射机i是EM场源时，
Figure GDA0001822578180000155
或
Figure GDA0001822578180000156
分别是由接收机j测得的EM场的z分量的实验仿真值或测得值，其中
Figure GDA0001822578180000157
以一般形式呈现为：
Figure GDA0001822578180000158
其中α、β和γ是系数，其为非零实数或零值，Ω是正则化算子，并且其中
Figure GDA0001822578180000159
是其自变量(argument)的函数。
在本发明的其他实施例中，目标因子可具有不同的组成。这将影响电磁干涉图像及其图案。对于α＝1、β＝0和γ＝0，下面给出了替代目标因子的一个示例：
Figure GDA00018225781800001510
其中，当发射机i是EM场的源时，
Figure GDA00018225781800001511
或
Figure GDA00018225781800001512
分别是由接收机j测得的EM场的z分量的实验仿真值或测得值，并且其中
Figure GDA0001822578180000161
是实验测得的EM场的z分量的最大范数。对于α＝1、β＝0和γ＝0，下面给出了替代目标因子145的另一个示例：
Figure GDA0001822578180000162
其中，当发射机i是EM场的源时，
Figure GDA0001822578180000163
或
Figure GDA0001822578180000164
分别是由接收机j测得的EM场的z分量的实验仿真值或测得值，并且其中当发射机i是幂为θ的EM场的源时，
Figure GDA0001822578180000165
是由接收机j测得的通过实验测得的EM场的z分量的范数。幂θ可以例如是以下值：1或2或3或1/2或甚至其他值。
目标因子确定框145的输出和未受干扰的EM干涉图像形成框130的输出(即，2D或3D的未受干扰的EM干涉图像135)被用于形成受干扰的EM干涉图像，如框150所示。根据应用，在框150中生成的图像是2D或3D图像。如框150所示，在成像域21中形成(3D)受干扰的EM干涉图像包括例如公式(6)中所示的和的计算：
Figure GDA0001822578180000166
其中
Figure GDA0001822578180000167
和
Figure GDA0001822578180000168
是来自于位于物理源(从1到N)的位置处和相应的物理接收机(从1到M)的位置处的频率为fk的EM源30的3D EM场(x，y，z)分布，作为共轭值，并且如在框145处确定的，Object Factori，j是从发射机i到接收机j的所选目标因子的第“ij”个分量。
受干扰的EM干涉图像形成框150的输出是受干扰的EM干涉图像155(或其代表的数据)。在这方面，图4A和4B是两个人的头部案例的示例性受干扰的EM干涉图像155。每个图像155是X-Y平面中的2D横截面视图，其中X和Y刻度以厘米为单位。(将会理解，尽管图4A和图4B是2D图示，但可通过对相同处理的进一步应用来实现3D图像。)在(4)处示出的替代性的目标因子被用于这两个示例。当比较图3的未受干扰的EM图像135和图4A和图4B的受干扰的EM图像155时，在图4A和4B中清楚地显示了对象410、420(人的头部的X-Y横截面)的介电结构与EM干涉图案415、425的叠加，其中每个图中的最外面的深色环代表人类颅骨，并且这样的深色环内的拓扑相似图案是剩余的EM干涉图案。
在框150中形成受干扰的EM干涉图像之后，将图像(数据)155发送到图案识别框160和叠加图像形成框170两者。干涉图案识别及其应用作为两阶段处理，其中阶段1是图案识别本身，阶段2是将识别的图案应用于2D或3D EM叠加图像，使EM干涉图案无效或衰减，并显示对象19的2D或3D介电结构。这两个阶段通过图案识别框160和叠加图像形成框170在图2中示出。
可以使用各种图案识别方法。例如，图5是用于图2中所示的图案识别框160的一个图案识别策略500和一组公式的流程图。在图5中使用2D示例，但可以理解的是，可以跨多个X-Y截面应用这种图案识别以创建3D图像，如图5所示。
在下文中有时被称为“第二T优化”的替代性图案识别策略有利地允许在叠加图像形成框170的初始迭代或第一次迭代中开发更清晰的图案，并且因此在后续迭代中使得图案识别更容易。第二T优化的关键特征是在迭代i>1时使用(6)的推广以形成叠加图像(为简单起见，频率项被省略)：
Figure GDA0001822578180000171
在一些实施例中，系数α可以独立于迭代次数(并且可以通过试算法找到)，而在其他实施例中，它可以具有迭代相关值，这有助于加速收敛过程。在后一种情况下，一种可能的策略是在迭代i+1时具有系数α的值作为迭代i时的图像与迭代i-1时的图像之间的最大绝对差异：
(8)max||ε(x，y，z)i-ε(x，y，z)i-1||.
在本发明的其他实施例中，可以使用其他图案识别策略和/或算法来找到最佳系数α。
在一些实施例中，可以在方法100的每次迭代期间执行相同的图案识别策略/过程，而在其他实施例中，图案识别策略/过程可以在每次连续迭代之后动态地更新。
将图案识别框160的输出和受干扰的EM干涉图像形成框150的输出(即，受干扰的EM干涉图像155)提供给叠加图像形成框170。可以通过经由框160识别图案、将它们与受干扰的EM干涉图像155进行比较、然后进行图案校正，来完成叠加图像形成170。在框170中，分析数据以使EM干涉图案衰减或无效并显示表示为图像输出200的对象19的真实2D或3D结构。
图案校正可以如下执行。首先，在方法100的至少一些实施例中，可以将整个成像域21细分为多个特征成像子域。在这方面，图6是示出将成像域21细分成子域的图像600。更具体地，使用人的头部的示例作为对象19，成像域21可以细分为外部(外面)域605、颅骨域610、脑脊髓液(CSF)域615、灰质域620、和白质域625。图7是可以针对每个子域在叠加图像形成框170中执行的一种可能的图案校正策略700的流程图。在步骤705，计算归一化的∫V{(6)或(7)的值}，其中V是特定子域的体积。在步骤710，计算归一化的∫L-pattern{(6)或(7)的值}，其中L是第“ip”个图案的长度。参考图3可以容易地理解步骤710的逻辑。可以理解的是，有清晰的圆形图案。那些圆形图案的数量Np等于半径上的网格的数量。因此，对于如图3中那些简单的圆形图案，步骤710的积分在从1到Np的每个图案(第“ip”个图案)上、在2πRNp的长度上进行。类似地，对于更复杂的对象图案(例如，图4A和4B中呈现的稍微环形的图案)，步骤710涉及计算归一化的∫L-pattern{(6)或(7)的值}，其中L是第“ip”个图案的长度。然后在步骤715处，通过将(6)或(7)的新值设置为等于旧值与{∫V{(6)或(7)}–∫L-pattern{(6)或(7)}的和来校正第ip个图案上的图像。然后针对每个第ip个图案重复步骤710和715。
应该理解的是，图7中所示的图案校正策略不是唯一可以利用的图案校正策略；在不脱离本发明的范围的情况下可以附加地或替代地使用其他图案校正策略。
在每次经过框170的迭代之后，在步骤180处针对“收敛目标”评估输出。当收敛目标时，生成最终输出图像200。
结合本文呈现的图案识别技术的重要益处可以理解如下。在进行如图2所示的迭代过程时，EM干涉图案的影响倾向于显著扭曲所得到的目标2D或3D介电结构的图像。这是因为通过在以下迭代中获得的相似EM干涉图案放大了这种图案，并由此产生了图像的失真。在这方面，图8A是来自于类似于图2的方法100的迭代成像方法的虚拟的中风患者的头部19的示例性输出图像，但是在图像形成期间图案识别框160及其输出信息未被使用，而图8B是与图8A相同的虚拟的中风患者的头部19的示例性输出图像200，但是在叠加图像形成框170期间使用了图案识别框160及其输出信息。图8B在图像右下方清楚地显示了虚拟中风区域820。在没有干涉图案识别框160及其在框170中的应用的情况下获得的图8A中的图像在虚拟中风区域810内显示了一些图像失真，但是在没有更清晰的情况下，只能推测失真代表了虚拟中风区域。
还将认识到，尽管本文描述的方法(包括变化和置换)已经主要根据2D和3D图像进行描述，但相信同样适用于4D EMT技术(空间3D加上时间上的1D)，其中例如在USPA 13/173,078中描述的那样生成4D动态融合电磁层析图像，其全部内容通过引用并入本文。
基于上述信息，本领域技术人员将容易理解，本发明易于广泛使用和应用。除了本文中具体描述的那些实施例以外，本发明的许多实施例和改编以及许多变化、修改和等同布置将从本发明及其前述描述显而易见或合理地提出，而不背离本发明的实质或范围。
因此，尽管已经在本文中关于一个或多个优选实施例详细描述了本发明，但应该理解的是，本公开仅仅是本发明的说明和示例，并且仅为了提供完整且可实现的本发明的公开的目的。前述公开内容不旨在解释为限制本发明或以其他方式排除任何此类其他实施例、改编、变化、修改或等同布置；本发明仅由所附权利要求及其等同物限制。

Claims (30)

1.一种在图像重建***中使用的电磁干涉图案识别层析成像方法，包括：
经由电磁层析成像***，生成与成像域中的对象相对应的电磁场数据，其中，所述电磁场数据在多个发射机处产生并且与所述对象相互作用之后，在多个接收机处测量所述电磁场数据；以及
以递归的方式重复地使用所生成的电磁场数据：
形成未受干扰的电磁干涉图像，
至少部分地基于所述未受干扰的电磁干涉图像来形成受干扰的电磁干涉图像，
识别重复形成的受干扰的电磁干涉图像中的电磁干涉图案，以及
通过使所识别的电磁干涉图案从所述受干扰的电磁干涉图像中无效或衰减来形成叠加图像。
2.根据权利要求1所述的方法，其中，所述至少部分地基于所述未受干扰的电磁干涉图像来形成受干扰的电磁干涉图像的步骤包括至少部分地基于对目标因子的确定来形成受干扰的电磁干涉图像，所述目标因子是在形成未受干扰的电磁干涉图像的步骤期间的实验电磁场与计算出的电磁场之间的差异的泛函。
3.根据权利要求2所述的方法，其中所述目标因子被确定为
Figure FDA0001822578170000011
其中，当发射机i是所述电磁场的源时，
Figure FDA0001822578170000012
或
Figure FDA0001822578170000013
分别是由接收机j测得的所述电磁场的z分量的实验仿真值或测得值，其中
Figure FDA0001822578170000014
一般以
Figure FDA0001822578170000015
的形式表示，其中α、β和γ是系数，所述系数为非零实数或零值，其中Ω是正则化算子，并且其中
Figure FDA0001822578170000016
是其自变量的函数。
4.根据权利要求2所述的方法，其中，所述目标因子被确定为
Figure FDA0001822578170000021
其中，当发射机i是所述电磁场的源时，
Figure FDA0001822578170000022
或
Figure FDA0001822578170000023
分别是由接收机j测得的所述电磁场的z分量的实验仿真值或测得值，并且其中
Figure FDA0001822578170000024
是所述电磁场的实验测得的z分量的最大范数。
5.根据权利要求2所述的方法，其中所述目标因子被确定为
Figure FDA0001822578170000025
其中，当发射机i是所述电磁场的源时，
Figure FDA0001822578170000026
或
Figure FDA0001822578170000027
分别是由接收机j测得的所述电磁场的z分量的实验仿真值或测得值，并且其中当发射机i是幂为θ的所述电磁场的源时，
Figure FDA0001822578170000028
是由接收机j测得的所述电磁场的实验测得的z分量的范数。
6.根据权利要求2所述的方法，其中，所述形成受干扰的电磁干涉图像的步骤包括计算
Figure FDA0001822578170000029
其中
Figure FDA00018225781700000210
和
Figure FDA00018225781700000211
是来自于所述电磁层析成像***中的位于物理源(从1到N)的位置处和相应的物理接收机(从1到M)的位置处的频率为fk的电磁源的3D电磁场(x，y，z)分布，作为共轭值，并且其中，Object Factori，j是从发射机i到接收机j的所述目标因子的第“ij”个分量。
7.根据权利要求2所述的方法，其中，所述识别重复形成的受干扰的电磁干涉图像中的电磁干涉图案的步骤包括计算总和
Figure FDA0001822578170000031
其中，
Figure FDA0001822578170000032
和
Figure FDA0001822578170000033
是来自于电磁层析成像***中的位于物理源(从1到N)的位置处和相应的物理接收机(从1到M)的位置处的频率为fk的电磁源的3D电磁场(x，y，z)分布，作为共轭值，并且其中，Object Factori，j是从发射机i到接收机j的所述目标因子的第“ij”个分量。
8.根据权利要求2所述的方法，其中，所述识别重复形成的受干扰的电磁干涉图像中的电磁干涉图案的步骤包括针对迭代i>1，计算
Figure FDA0001822578170000034
其中，为了简便省略了频率项，其中，
Figure FDA0001822578170000035
和
Figure FDA0001822578170000036
是来自于电磁层析成像***中的位于物理源(从1到N)的位置处和相应的物理接收机(从1到M)的位置处的电磁源的3D电磁场(x，y，z)分布，作为共轭值，并且其中，Object Factori，j是从发射机i到接收机j的所述目标因子的第“ij”个分量。
9.根据权利要求1所述的方法，还包括在形成叠加图像的每个重复步骤之后执行的确定是否已经达到收敛目标的步骤。
10.根据权利要求1所述的方法，其中，所述方法用作生成4D差分(动态)融合图像的方法的一部分。
11.根据权利要求10所述的方法，其中，生成4D差分(动态)融合图像包括将指示相对生理变化的至少一个连续形成的图像与基线解剖图像相组合以显示为单个统一图像。
12.根据权利要求11所述的方法，其中，生成4D差分(动态)融合图像的方法被用作利用4D动态融合电磁图案识别层析成像来监测生物组织的生存能力和/或功能状况的方法的一部分。
13.根据权利要求1所述的方法，其中按顺序地执行以下步骤：形成未受干扰的电磁干涉图像，至少部分地基于所述未受干扰的电磁干涉图像形成受干扰的电磁干涉图像，识别重复形成的受干扰的电磁干涉图像中的电磁干涉图案，以及通过使所识别的电磁干涉图案从所述受干扰的电磁干涉图像中无效或衰减来形成叠加图像。
14.根据权利要求1所述的方法，还包括经由显示单元显示所述叠加图像的步骤。
15.一种使用电磁干涉图案识别层析成像来重建图像的方法，包括：
经由电磁层析成像***，生成与成像域中的对象相对应的电磁场数据，其中，所述电磁场数据在多个发射机处产生并且与所述对象相互作用之后，在多个接收机处测量所述电磁场数据；以及
以递归的方式重复地使用所生成的电磁场数据：
形成未受干扰的电磁干涉图像，
至少部分地基于所述未受干扰的电磁干涉图像来形成受干扰的电磁干涉图像，
识别重复形成的受干扰的电磁干涉图像中的电磁干涉图案；以及
通过使所识别的电磁干涉图案从所述受干扰的电磁干涉图像中无效或衰减来形成叠加图像。
16.根据权利要求15所述的方法，其中，所述至少部分地基于所述未受干扰的电磁干涉图像来形成受干扰的电磁干涉图像的步骤包括至少部分地基于对目标因子的确定来形成受干扰的电磁干涉图像，所述目标因子是在形成未受干扰的电磁干涉图像的步骤期间的实验电磁场与计算出的电磁场之间的差异的泛函。
17.根据权利要求15所述的方法，还包括在形成叠加图像的每个重复步骤之后执行的确定是否已经达到收敛目标的步骤。
18.根据权利要求15所述的方法，其中，所述方法用作生成4D差分(动态)融合图像的方法的一部分。
19.根据权利要求18所述的方法，其中，生成4D差分(动态)融合图像包括将指示相对生理变化的至少一个连续形成的图像与基线解剖图像相组合以显示为单个统一图像。
20.根据权利要求19所述的方法，其中，生成4D差分(动态)融合图像的方法被用作利用4D动态融合电磁图案识别层析成像来监测生物组织的生存能力和/或功能状况的方法的一部分。
21.根据权利要求15所述的方法，其中按顺序地执行以下步骤：形成未受干扰的电磁干涉图像，至少部分地基于所述未受干扰的电磁干涉图像形成受干扰的电磁干涉图像，识别重复形成的受干扰的电磁干涉图像中的电磁干涉图案，以及通过使所识别的电磁干涉图案从所述受干扰的电磁干涉图像中无效或衰减来形成叠加图像。
22.根据权利要求15所述的方法，还包括经由显示单元显示所述叠加图像的步骤。
23.一种使用电磁干涉图案识别层析成像的图像重建***，包括：
电磁层析成像***，其生成与成像域中的对象相对应的电磁场数据，所述电磁层析成像***包括：
多个电磁发射机；
多个接收机，其测量在所述多个发射机处产生并与所述对象相互作用之后的电磁场数据；
以及边界设备；以及
处理中心，其以递归的方式重复地使用所生成的电磁场数据，执行以下步骤：
形成未受干扰的电磁干涉图像，
至少部分地基于所述未受干扰的电磁干涉图像来形成受干扰的电磁干涉图像，
识别重复形成的受干扰的电磁干涉图像中的电磁干涉图案，以及
通过使所识别的电磁干涉图案从所述受干扰的电磁干涉图像中无效或衰减来形成叠加图像。
24.根据权利要求23所述的图像重建***，其中，所述至少部分地基于未受干扰的电磁干涉图像来形成受干扰的电磁干涉图像的步骤包括至少部分地基于对目标因子的确定来形成受干扰的电磁干涉图像，所述目标因子是在形成未受干扰的电磁干涉图像的步骤期间的实验电磁场与计算出的电磁场之间的差异的泛函。
25.根据权利要求23所述的图像重建***，其中，所述处理中心还在形成叠加图像的每个重复步骤之后执行确定是否已经达到收敛目标的步骤。
26.根据权利要求23所述的图像重建***，其中，由所述处理中心执行的步骤用作生成4D差分(动态)融合图像的方法的一部分。
27.根据权利要求26所述的图像重建***，其中，生成4D差分(动态)融合图像包括将指示相对生理变化的至少一个连续形成的图像与基线解剖图像相组合以显示为单个统一图像。
28.根据权利要求27所述的图像重建***，其中，生成4D差分(动态)融合图像的方法被用作利用4D动态融合电磁图案识别层析成像来监测生物组织的生存能力和/或功能状况的方法的一部分。
29.根据权利要求23所述的图像重建***，其中，按顺序地执行以下步骤：形成未受干扰的电磁干涉图像，至少部分地基于未受干扰的电磁干涉图像形成受干扰的电磁干涉图像，识别重复形成的受干扰的电磁干涉图像中的电磁干涉图案，以及通过使所识别的电磁干涉图案从所述受干扰的电磁干涉图像中无效或衰减来形成叠加图像。
30.根据权利要求23所述的图像重建***，还包括显示所述叠加图像的显示单元。
CN201680074323.2A 2015-10-16 2016-10-16 电磁干涉图案识别层析成像 Active CN108366771B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201562242915P 2015-10-16 2015-10-16
US62/242,915 2015-10-16
PCT/US2016/057254 WO2017066731A1 (en) 2015-10-16 2016-10-16 Electromagnetic interference pattern recognition tomography

Publications (2)

Publication Number Publication Date
CN108366771A CN108366771A (zh) 2018-08-03
CN108366771B true CN108366771B (zh) 2021-11-02

Family

ID=58518250

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201680074323.2A Active CN108366771B (zh) 2015-10-16 2016-10-16 电磁干涉图案识别层析成像

Country Status (8)

Country Link
US (2) US10921361B2 (zh)
EP (1) EP3361955B1 (zh)
CN (1) CN108366771B (zh)
DK (1) DK3361955T3 (zh)
ES (1) ES2825898T3 (zh)
IL (1) IL258655B (zh)
RU (1) RU2720161C2 (zh)
WO (1) WO2017066731A1 (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140275944A1 (en) 2013-03-15 2014-09-18 Emtensor Gmbh Handheld electromagnetic field-based bio-sensing and bio-imaging system
DK3361955T3 (da) 2015-10-16 2020-10-26 Emtensor Gmbh Elektromagnetisk interferensmønstergenkendelsestomografi
CA3044844A1 (en) 2016-11-23 2018-05-31 Emtensor Gmbh Use of electromagnetic field for tomographic imaging of head
JP7397067B2 (ja) * 2018-09-04 2023-12-12 エムビジョン・メディカル・デバイシーズ・リミテッド 医療用イメージングのための装置及び処理
EP4301219A1 (en) * 2021-03-04 2024-01-10 Emvision Medical Devices Ltd Apparatus and process for electromagnetic imaging

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102132149A (zh) * 2008-06-27 2011-07-20 沃尔弗拉姆·R·雅里施 高效计算机断层摄影技术

Family Cites Families (91)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4157472A (en) 1976-09-16 1979-06-05 General Electric Company X-ray body scanner for computerized tomography
US4135131A (en) 1977-10-14 1979-01-16 The United States Of America As Represented By The Secretary Of The Army Microwave time delay spectroscopic methods and apparatus for remote interrogation of biological targets
US4247815A (en) 1979-05-22 1981-01-27 The United States Of America As Represented By The Secretary Of The Army Method and apparatus for physiologic facsimile imaging of biologic targets based on complex permittivity measurements using remote microwave interrogation
US4257278A (en) 1979-08-24 1981-03-24 General Electric Company Quantitative volume blood flow measurement by an ultrasound imaging system featuring a Doppler modality
US4638813A (en) 1980-04-02 1987-01-27 Bsd Medical Corporation Electric field probe
US4662222A (en) 1984-12-21 1987-05-05 Johnson Steven A Apparatus and method for acoustic imaging using inverse scattering techniques
DE3531893A1 (de) 1985-09-06 1987-03-19 Siemens Ag Verfahren zur bestimmung der verteilung der dielektrizitaetskonstanten in einem untersuchungskoerper sowie messanordnung zur durchfuehrung des verfahrens
DE3601983A1 (de) 1986-01-23 1987-07-30 Siemens Ag Verfahren und vorrichtung zur beruehrungslosen bestimmung der temperaturverteilung in einem untersuchungsobjekt
US4926868A (en) 1987-04-15 1990-05-22 Larsen Lawrence E Method and apparatus for cardiac hemodynamic monitor
US5069223A (en) 1990-02-14 1991-12-03 Georgetown University Method of evaluating tissue changes resulting from therapeutic hyperthermia
US5363050A (en) 1990-08-31 1994-11-08 Guo Wendy W Quantitative dielectric imaging system
US5233713A (en) 1991-03-27 1993-08-10 General Electric Company Head holder for nuclear imaging
US5222501A (en) 1992-01-31 1993-06-29 Duke University Methods for the diagnosis and ablation treatment of ventricular tachycardia
US5305748A (en) 1992-06-05 1994-04-26 Wilk Peter J Medical diagnostic system and related method
US5263050A (en) 1992-09-09 1993-11-16 Echelon Corporation Adaptive threshold in a spread spectrum communications system
US5405346A (en) 1993-05-14 1995-04-11 Fidus Medical Technology Corporation Tunable microwave ablation catheter
ZA948393B (en) 1993-11-01 1995-06-26 Polartechnics Ltd Method and apparatus for tissue type recognition
US5715819A (en) 1994-05-26 1998-02-10 The Carolinas Heart Institute Microwave tomographic spectroscopy system and method
AUPM851694A0 (en) 1994-09-30 1994-10-27 Barsamian, Sergei T New methods for diagnosis, detection of cell abnormalities and morphology of living systems
US6026173A (en) 1997-07-05 2000-02-15 Svenson; Robert H. Electromagnetic imaging and therapeutic (EMIT) systems
CN1261259A (zh) 1997-05-23 2000-07-26 卡罗莱纳心脏研究所 电子成像和治疗***
CA2303073C (en) 1997-09-11 2004-07-13 West Virginia University Research Corp. Electrical property enhanced tomography (epet) apparatus and method
US6697660B1 (en) 1998-01-23 2004-02-24 Ctf Systems, Inc. Method for functional brain imaging from magnetoencephalographic data by estimation of source signal-to-noise ratio
US6703838B2 (en) * 1998-04-13 2004-03-09 Schlumberger Technology Corporation Method and apparatus for measuring characteristics of geological formations
US6333087B1 (en) 1998-08-27 2001-12-25 Chevron Chemical Company Llc Oxygen scavenging packaging
US6233479B1 (en) 1998-09-15 2001-05-15 The Regents Of The University Of California Microwave hematoma detector
US6454711B1 (en) 1999-04-23 2002-09-24 The Regents Of The University Of California Microwave hemorrhagic stroke detector
US7373197B2 (en) 2000-03-03 2008-05-13 Intramedical Imaging, Llc Methods and devices to expand applications of intraoperative radiation probes
US6511427B1 (en) 2000-03-10 2003-01-28 Acuson Corporation System and method for assessing body-tissue properties using a medical ultrasound transducer probe with a body-tissue parameter measurement mechanism
US6481887B1 (en) 2000-04-12 2002-11-19 Ge Medical Systems Global Technology Company, Llc Emergency vehicle with medical image scanner and teleradiology system and method of operation
US6503203B1 (en) 2001-01-16 2003-01-07 Koninklijke Philips Electronics N.V. Automated ultrasound system for performing imaging studies utilizing ultrasound contrast agents
US6965661B2 (en) 2001-06-19 2005-11-15 Hitachi, Ltd. Radiological imaging apparatus and radiological imaging method
US7570063B2 (en) 2001-07-06 2009-08-04 Wisconsin Alumni Research Foundation Space-time microwave imaging for cancer detection
EP1436602A2 (en) 2001-10-12 2004-07-14 Wisconsin Alumni Research Foundation Microwave dielectric spectroscopy method and apparatus
US6865494B2 (en) 2001-12-18 2005-03-08 Mri Devices Corp. Method and apparatus for noise tomography
US7164105B2 (en) 2002-04-05 2007-01-16 Microwave Imaging Systems Technologies, Inc. Non-invasive microwave analysis systems
US8892189B2 (en) 2002-05-30 2014-11-18 Alcatel Lucent Apparatus and method for heart size measurement using microwave doppler radar
US7272431B2 (en) 2002-08-01 2007-09-18 California Institute Of Technology Remote-sensing method and device
AU2003302229A1 (en) 2002-09-27 2004-06-30 The Trustees Of Dartmouth College Imaging by magnetic resonance adsorption, elastography and tomography
US7239731B1 (en) 2002-11-26 2007-07-03 Emimaging Ltd System and method for non-destructive functional imaging and mapping of electrical excitation of biological tissues using electromagnetic field tomography and spectroscopy
WO2004062467A2 (en) 2002-12-31 2004-07-29 John Herbert Cafarella Multi-sensor breast tumor detection
WO2004073618A2 (en) 2003-02-14 2004-09-02 University Of Florida Breast cancer detection system
JP4263579B2 (ja) 2003-10-22 2009-05-13 アロカ株式会社 超音波診断装置
US20050135560A1 (en) 2003-12-17 2005-06-23 Ehud Dafni Portable computed tomography scanner and methods thereof
EP1748726B1 (en) 2004-05-26 2010-11-24 Medical Device Innovations Limited Tissue detection and ablation apparatus
US7782998B2 (en) 2004-12-21 2010-08-24 General Electric Company Method and apparatus for correcting motion in image reconstruction
US8253619B2 (en) 2005-02-15 2012-08-28 Techtronic Power Tools Technology Limited Electromagnetic scanning imager
US7962198B2 (en) 2005-04-27 2011-06-14 The Trustees Of Dartmouth College System and method for spectral-encoded high-rate hemodynamic tomography
DE102005026940A1 (de) 2005-06-06 2006-12-14 Schering Ag Röntgenanordnung zur Bilddarstellung eines Untersuchungsobjektes und Verwendung der Röntgenanordnung
US7387607B2 (en) 2005-06-06 2008-06-17 Intel Corporation Wireless medical sensor system
US20070025514A1 (en) 2005-06-06 2007-02-01 Ruediger Lawaczeck X-ray arrangement for graphic display of an object under examination and use of the x-ray arrangement
EP1968431B2 (en) 2005-12-22 2022-06-01 Visen Medical, Inc. Combined x-ray and optical tomographic imaging system
EP2020915B1 (en) 2006-05-22 2017-04-26 Medfields Diagnostics AB System relating to examination of an object
US9332922B2 (en) 2006-06-29 2016-05-10 Medfields Diagnostics Ab Solution for internal monitoring of body
US9709733B2 (en) 2006-09-28 2017-07-18 The Florida International University Board Of Trustees Hand-held optical probe based imaging system with 3D tracking facilities
US8423152B2 (en) 2007-05-14 2013-04-16 Bsd Medical Corporation Apparatus and method for selectively heating a deposit in fatty tissue in a body
US8089417B2 (en) 2007-06-01 2012-01-03 The Royal Institution For The Advancement Of Learning/Mcgill University Microwave scanning system and miniaturized microwave antenna
US7876114B2 (en) 2007-08-08 2011-01-25 Cascade Microtech, Inc. Differential waveguide probe
GB0721694D0 (en) 2007-11-05 2007-12-12 Univ Bristol Methods and apparatus for measuring the contents of a search volume
GB0721693D0 (en) 2007-11-05 2007-12-12 Univ Bristol Antenna for investigating structure of human or animal
US9672471B2 (en) 2007-12-18 2017-06-06 Gearbox Llc Systems, devices, and methods for detecting occlusions in a biological subject including spectral learning
US8708919B1 (en) 2008-01-25 2014-04-29 Raytheon Company System and method for remotely sensing vital signs
EP3181045B1 (en) 2009-03-04 2020-12-02 Sensible Medical Innovations Ltd. Methods and systems for monitoring intrabody tissues
WO2010143691A1 (ja) 2009-06-10 2010-12-16 国立大学法人静岡大学 診断装置
WO2011009945A2 (en) 2009-07-23 2011-01-27 Medfield Diagnostics Ab Classification of microwave scattering data
KR20130136907A (ko) 2010-06-01 2013-12-13 어뎁턴스, 엘엘씨 네트워크 연결된 착용가능한 의료용 센서를 위한 시스템 및 방법
US9635349B2 (en) 2010-06-11 2017-04-25 The Florida International University Board Of Trustees Second generation hand held optical imager
US9724010B2 (en) 2010-07-08 2017-08-08 Emtensor Gmbh Systems and methods of 4D electromagnetic tomographic (EMT) differential (dynamic) fused imaging
US20120083690A1 (en) * 2010-10-04 2012-04-05 Semenov Serguei Y System and method for electromagnetic imaging and therapeutics using specialized nanoparticles
US9079011B2 (en) 2011-01-03 2015-07-14 Wisconsin Alumni Research Foundation Microwave hyperthermia treatment system
US9649091B2 (en) 2011-01-07 2017-05-16 General Electric Company Wireless ultrasound imaging system and method for wireless communication in an ultrasound imaging system
US8376948B2 (en) 2011-02-17 2013-02-19 Vivant Medical, Inc. Energy-delivery device including ultrasound transducer array and phased antenna array
US9448187B2 (en) 2011-07-01 2016-09-20 University Of Manitoba Imaging using probes
DE102011106405B4 (de) * 2011-07-02 2021-08-12 Drägerwerk AG & Co. KGaA Elektroimpedanztomographie-Gerät
CN107174214B (zh) 2011-12-20 2021-03-30 明智医疗创新有限公司 定位电磁传感器的胸廓服装及使用所述胸廓服装的方法
US10743815B2 (en) 2012-01-19 2020-08-18 Cerebrotech Medical Systems, Inc. Detection and analysis of spatially varying fluid levels using magnetic signals
EP3571986A1 (en) 2012-01-19 2019-11-27 Cerebrotech Medical Systems, Inc. Diagnostic system for detection of fluid changes
DE102012205294B3 (de) 2012-03-30 2013-06-27 Siemens Aktiengesellschaft Verfahren und Steuereinrichtung zur Ansteuerung eines Magnetresonanzsystems
US9445729B2 (en) 2012-07-20 2016-09-20 Resmed Sensor Technologies Limited Range gated radio frequency physiology sensor
CA2936145C (en) 2012-11-21 2021-06-15 Emtensor Gmbh Electromagnetic tomography solutions for scanning head
US9072449B2 (en) 2013-03-15 2015-07-07 Emtensor Gmbh Wearable/man-portable electromagnetic tomographic imaging
US20140275944A1 (en) 2013-03-15 2014-09-18 Emtensor Gmbh Handheld electromagnetic field-based bio-sensing and bio-imaging system
US9357976B2 (en) * 2013-10-24 2016-06-07 General Electric Company System and method of noise deletion in computed tomography
US9420977B2 (en) 2014-03-19 2016-08-23 Tribogenics, Inc. Portable head CT scanner
CN112998685A (zh) 2014-09-03 2021-06-22 脑科技医疗***股份有限公司 利用磁信号检测和分析空间变化的流体水平
DK3361955T3 (da) 2015-10-16 2020-10-26 Emtensor Gmbh Elektromagnetisk interferensmønstergenkendelsestomografi
EP3551070A1 (en) 2016-12-06 2019-10-16 Medfields Diagnostics AB System and method for detecting an assymetrically positioned internal object in a body
WO2018127434A1 (en) 2017-01-09 2018-07-12 Medfield Diagnostics Ab Method and system for ensuring antenna contact and system function in applications of detecting internal dielectric properties in a body
US11551386B2 (en) 2017-06-08 2023-01-10 EMvision Medical Devices Ltd Tomographic imaging system and process
WO2019094877A1 (en) 2017-11-13 2019-05-16 Cerebrotech Medical Systems, Inc. Continuous autoregulation system
EP3797453A1 (en) 2018-05-23 2021-03-31 Medfield Diagnostics AB Solution for absorption of microwaves

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102132149A (zh) * 2008-06-27 2011-07-20 沃尔弗拉姆·R·雅里施 高效计算机断层摄影技术

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Microwave Tomography: Theoretical and Experimental Investigation of the Iteration Reconstruction Algorithm;Serguei Y. Semenov, et al;《IEEE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES》;19980228;第46卷(第2期);133-141 *
基于差分进化算法的多模态医学图像融合;许良凤, 林辉, 胡敏;《电子测量与仪器学报》;20130228;第27卷(第2期);110-114 *

Also Published As

Publication number Publication date
ES2825898T3 (es) 2021-05-17
EP3361955A4 (en) 2019-05-22
RU2720161C2 (ru) 2020-04-24
RU2018117885A3 (zh) 2020-02-17
DK3361955T3 (da) 2020-10-26
EP3361955B1 (en) 2020-09-02
IL258655A (en) 2018-06-28
IL258655B (en) 2021-05-31
US20180231594A1 (en) 2018-08-16
CN108366771A (zh) 2018-08-03
EP3361955A1 (en) 2018-08-22
WO2017066731A1 (en) 2017-04-20
RU2018117885A (ru) 2019-11-21
US20210181246A1 (en) 2021-06-17
US11892491B2 (en) 2024-02-06
WO2017066731A8 (en) 2018-05-31
US10921361B2 (en) 2021-02-16

Similar Documents

Publication Publication Date Title
CN108366771B (zh) 电磁干涉图案识别层析成像
US6373998B2 (en) Image processing electronic device for detecting dimensional variations
CN111771138A (zh) 具有用于根据欠采样超声数据产生图像的神经网络的超声***
US9542763B2 (en) Systems and methods for fast reconstruction for quantitative susceptibility mapping using magnetic resonance imaging
JP4527881B2 (ja) 磁気共鳴イメージング方法
Kunyansky et al. Rotational magneto-acousto-electric tomography (MAET): Theory and experimental validation
Galande et al. The art of medical image fusion: A survey
US20180180688A1 (en) System and Method for Localized Processing of Quantitative Susceptibility Maps in Magnetic Resonance Imaging
JP2013063272A (ja) Ｍｒｉの自動３次元分割方法
Prasad et al. DeepUCT: Complex cascaded deep learning network for improved ultrasound tomography
Liang et al. Coplanar electrical/ultrasonic dual-modality tomography for water continuous gas/oil/water three-phase distribution imaging
Alessandrini et al. Multiscale optical flow computation from the monogenic signal
Lu et al. Reconstruction of elasticity: a stochastic model-based approach in ultrasound elastography
Mamatjan Imaging of hemorrhagic stroke in magnetic induction tomography: An in vitro study
Lee et al. High-frequency conductivity at Larmor-frequency in human brain using moving local window multilayer perceptron neural network
Joshi et al. Multi-Modality Medical Image Fusion Using SWT & Speckle Noise Reduction with Bidirectional Exact Pattern Matching Algorithm
Hu et al. Data-driven system matrix manipulation enabling fast functional imaging and intra-image nonrigid motion correction in tomography
Hong et al. Accelerating compressed sensing reconstruction of subsampled radial k-space data using geometrically-derived density compensation
Gou et al. Lung dynamic MRI deblurring using low‐rank decomposition and dictionary learning
Otani et al. Assessing cardiac tissue function via action potential wave imaging using cardiac displacement data
Moerman et al. Evaluation of a validation method for MR imaging-based motion tracking using image simulation
De Tillieux et al. A volumetric integral equation formulation for magnetic induction tomography
Filist Two-Dimensional Walsh Spectral Transform in Problems of Automated Analysis of Ultrasound Images
Murillo et al. A review of motion estimation methods for non-invasive ultrasound motion and emerging strain imaging methods of carotid artery plaques
Grau et al. Motion-guided anisotropic filtering of ultrasound sequences

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
GR01 Patent grant
GR01 Patent grant