[go: up one dir, main page]

RU2720161C2 - Электромагнитная томография с распознаванием картин интерференции - Google Patents

Электромагнитная томография с распознаванием картин интерференции Download PDF

Info

Publication number
RU2720161C2
RU2720161C2 RU2018117885A RU2018117885A RU2720161C2 RU 2720161 C2 RU2720161 C2 RU 2720161C2 RU 2018117885 A RU2018117885 A RU 2018117885A RU 2018117885 A RU2018117885 A RU 2018117885A RU 2720161 C2 RU2720161 C2 RU 2720161C2
Authority
RU
Russia
Prior art keywords
image
electromagnetic
electromagnetic interference
interference
perturbed
Prior art date
Application number
RU2018117885A
Other languages
English (en)
Other versions
RU2018117885A (ru
RU2018117885A3 (ru
Inventor
Сергей И. СЕМЕНОВ
Original Assignee
Эмтензор Гмбх
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 Эмтензор Гмбх filed Critical Эмтензор Гмбх
Publication of RU2018117885A publication Critical patent/RU2018117885A/ru
Publication of RU2018117885A3 publication Critical patent/RU2018117885A3/ru
Application granted granted Critical
Publication of RU2720161C2 publication Critical patent/RU2720161C2/ru

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0033Features or image-related aspects of imaging apparatus, e.g. for MRI, optical tomography or impedance tomography apparatus; Arrangements of imaging apparatus in a room
    • A61B5/004Features or image-related aspects of imaging apparatus, 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/0042Features or image-related aspects of imaging apparatus, 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R29/00Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
    • G01R29/08Measuring electromagnetic field characteristics
    • G01R29/0864Measuring electromagnetic field characteristics characterised by constructional or functional features
    • G01R29/0892Details 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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B34/00Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
    • A61B34/20Surgical navigation systems; Devices for tracking or guiding surgical instruments, e.g. for frameless stereotaxis
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
    • A61B5/053Measuring electrical impedance or conductance of a portion of the body
    • A61B5/0536Impedance imaging, e.g. by tomography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B34/00Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
    • A61B34/20Surgical navigation systems; Devices for tracking or guiding surgical instruments, e.g. for frameless stereotaxis
    • A61B2034/2046Tracking techniques
    • A61B2034/2051Electromagnetic tracking systems
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
    • A61B5/0507Detecting, 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)
  • Surgery (AREA)
  • Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Public Health (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Veterinary Medicine (AREA)
  • Animal Behavior & Ethology (AREA)
  • Pathology (AREA)
  • Biophysics (AREA)
  • Neurology (AREA)
  • Robotics (AREA)
  • Electromagnetism (AREA)
  • General Physics & Mathematics (AREA)
  • Optics & Photonics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

Изобретение относится к области медицины. Способ электромагнитной томографии с распознаванием картин интерференции (EMIPRT) для использования в системе реконструкции изображения содержит генерацию данных электромагнитного поля, соответствующих объекту в области визуализации, посредством системы электромагнитной томографии, причем данные электромагнитного поля измеряют на множестве приемников после его создания на множестве передатчиков и взаимодействия с объектом; и использование сгенерированных данных электромагнитного поля повторно, рекурсивным образом, с формированием невозмущенного изображения электромагнитной интерференции, формированием возмущенного изображения электромагнитной интерференции, основанного по меньшей мере частично на невозмущенном изображении электромагнитной интерференции, распознаванием картин электромагнитной интерференции в повторно сформированных возмущенных изображениях электромагнитной интерференции и формированием суперпозиционного изображения путем сведения к нулю или уменьшения распознанных картин электромагнитной интерференции от возмущенного изображения электромагнитной интерференции. Способ реконструкции изображения с использованием электромагнитной томографии с распознаванием картин интерференции содержит генерацию данных электромагнитного поля, соответствующих объекту в области визуализации, посредством системы электромагнитной томографии, причем данные электромагнитного поля измеряют на множестве приемников после его создания на множестве передатчиков и взаимодействия с объектом; и использование сгенерированных данных электромагнитного поля повторно, рекурсивным образом, с формированием невозмущенного изображения электромагнитной интерференции, формированием возмущенного изображения электромагнитной интерференции, основанного по меньшей мере частично на невозмущенном изображении электромагнитной интерференции, распознаванием картин электромагнитной интерференции в повторно сформированных возмущенных изображениях электромагнитной интерференции, и формированием суперпозиционного изображения путем сведения к нулю или уменьшения распознанных картин электромагнитной интерференции от возмущенного изображения электромагнитной интерференции. Технический результат заявленного изобретения состоит в увеличении точности представления изображений. 3 н. и 27 з.п. ф-лы, 8 ил.

Description

ПЕРЕКРЕСТНАЯ ССЫЛКА НА СВЯЗАННЫЕ ЗАЯВКИ
[0001] В отношении Соединенных Штатов настоящая заявка является обычной патентной заявкой США и испрашивает приоритет согласно 35 U.S.C. 119 (e) предварительной заявки на патент США с порядковым номером 62/242,915, поданной 16 октября 2015 года и озаглавленной «ELECTROMAGNETIC INTERFERENCE PATTERN RECOGNITION TOMOGRAPHY», причем заявка '915 в полном объеме явно включена в настоящий документ по ссылке. Кроме того, каждый из следующих патентов, патентных заявок и публикаций патентных заявок включен в настоящий документ по ссылки в полном объеме:
(a) патент США № 9414749 Семенова, поданный 16 августа 2016 года и ранее опубликованный 5 июня 2014 как публикация патентной заявки США № 2014/0155740 A1, которая предназначена по меньшей мере для предоставления информации об уровне техники и технической информации в отношении систем и условий изобретений настоящей патентной заявки; а также
(b) публикация патентной заявки США № 2012/0010493 A1, которая была опубликована 12 января 2012 года на основе патентной заявки № 13/173,078 Семенова, поданной 30 июня 2011 года и озаглавленной «SYSTEMS AND METHODS OF ELECTROMAGNETIC TOMOGRAPHY (EMT) DIFFERENTIAL (DYNAMIC) FUSED IMAGING», которая предназначена для предоставления информации об уровне техники и и технической информации в отношении 4D EMT визуализации.
ЗАЯВЛЕНИЕ ОБ АВТОРСКОМ ПРАВЕ
[0002] Все материалы в настоящей патентной заявке охраняются законодательством об авторских и смежных правах Соединенных Штатов Америки и других стран. Владелец авторских прав не возражает против факсимильного воспроизведения патентной заявки любым лицом или раскрытия патентной информации для отражения в государственных регистрационных документах, но в противном случае все прочие авторские права сохраняются за владельцем.
ПРЕДПОСЫЛКИ НАСТОЯЩЕГО ИЗОБРЕТЕНИЯ
Область настоящего изобретения
[0003] Настоящее изобретение в целом относится к электромагнитной томографии и, в частности, к использованию стратегий распознавания картин электромагнитной интерференции для устранения усиленных искажений, вызванных повторяющимися картинами интерференции.
Уровень техники
[0004] Электромагнитная томография (ЭМТ) - сравнительно новая модальность визуализации с большим потенциалом как для биомедицинских, так и для промышленных применений. Биомедицинские применения включают в себя, но не ограничиваются, неинвазивной оценкой функциональных и патологических состояний биологических тканей. Промышленные применения включают в себя, но не ограничиваются, разведкой нефти и газа, поиском и оценкой мин и оценкой потока в неметаллических трубах. Используя ЭМТ, различаются объекты, такие как биологические ткани и, следовательно, могут быть отображены на основе различий в диэлектрических свойствах таких объектов. Считается, что ЭМТ обладает высоким потенциалом для биомедицинских применений, что основано на недавней демонстрации зависимости диэлектрических свойств ткани от различных функциональных и патологических состояний ткани, таких как содержание крови и кислорода, ишемия и злокачественное развитие инфаркта, отеки и другие.
[0005] В течение последнего десятилетия или более были разработаны двумерные (2D), трехмерные (3D) и даже «четырехмерные» (4D) системы EMT и способы реконструкции изображений. Была продемонстрирована осуществимость технологии для различных биомедицинских и промышленных применений, например, для визуализации сердца и конечностей.
[0006] Как и при любой томографической визуализации классический сценарий визуализации ЭМТ состоит из циклов измерений сложных сигналов в зависимости от наличия исследуемого объекта, находящегося в так называемой области визуализации, создаваемых множеством передатчиков, расположенных в разных точках вокруг объекта и измеряемых на множестве приемников, расположенных в разных точках вокруг объекта. Это проиллюстрировано на ФИГ. 1. Местоположения передатчиков и приемников могут быть в области визуализации, на границе области визуализации или вне области визуализации. Как указано в другом месте настоящего документа, измеренная матрица ЭМ-сигналов затем может использоваться способами реконструкции изображения для реконструкции 3D распределения диэлектрических свойств объекта 19, т.е. для построения 3D изображения объекта. Компоненты оборудования, раскрытые и описанные в указанном выше патенте США № 9414749, могут использоваться для генерации необходимых EM-сигналов для реконструкции 3D изображения объекта 19.
[0007] Как правило, для реконструкции изображения очень важно точно описать распределение ЭМ поля по отношению к области 21 визуализации. Распределение ЭМ поля в камере визуализации - очень сложное явление, даже если внутри него нет исследуемого объекта.
[0008] Данное изобретение относится к использованию электромагнитных (ЭМ) полей для визуализации структуры объекта 19. Объект 19 может быть человеческим телом или частью человеческого тела, таким как голова, туловище, рука или тому подобное, но также может быть любым объектом без металлического экранирования. Использование ЭМ полей для визуализации внутри сильно экранированного объекта (но экранированного не посредством металла) является проблемой высокой сложности. Одним из примеров такого применения является визуализация человеческого мозга. Однако следует принимать во внимание, что другие такие приложения могут включать в себя визуализацию любой ткани человека, которая экранирована костной структурой. ЭМ визуализация мозга или других тканей, окруженных костью, представляет собой очень сложную проблему, учитывая высокий диэлектрический контраст. Задача состоит в том, чтобы реконструировать скрытые свойства глубоких тканей мозга, которые эффективно экранированы экраном с высоким диэлектрическим контрастом, в том числе черепом (с диэлектрическими свойствами в диапазоне 16+j5) и цереброспинальной жидкостью (с диэлектрическими свойствами в диапазоне 60+j30). Хотя это изобретение, как указано выше, применимо для визуализации любых объектов, можно полагать, что оно особенно применимо для визуализации внутри сильно экранированных объектов.
[0009] ЭМТ визуализация объектов с высоким диэлектрическим контрастом, в том числе биологических объектов, имеет очень сложную проблему так называемой «дифракционной томографии». Devaney A.J. «Current research topics in diffraction tomography», в Inverse Problmes in Scattering and Imaging, M. BNertero и E.R. Pike, Eds, New York: Adam Hilger, 1992, стр.47-58. Высокий диэлектрический контраст между тканями с высоким содержанием воды, такими как, мышечная ткань, но не ограничиваясь ею, и с низким содержанием воды, например, кость, но не ограничиваясь ею, является дополнительным осложнением при использовании ЭМ полей для визуализации. В попытке решить проблему дифракционной томографии были разработаны различные подходы в 2D и 3D-геометрии с использованием скалярных и векторных аппроксимаций. Смотри H. Harada, D. Wall, T. Takenaka, и T. Tanaka «Conjugate gradient method applied to inverse scattering problem», IEEE Trans. Antennas and Propagations, том. 43, 784-792, Aug. 1995; R.E. Kleinman, и P. M. van den Berg «A modified gradient method for two-dimensional problems in tomography», J.Comput. Appl. Math., том 42, стр. 17-35, Январь 1992; A. Abubakar, P.M. van den Berg, и J.J. Mallorqui «Imaging of Biomedical Data Using A Multiplicative Regularized Source Inversion Method» IEEE Trans. Microwave Theory and Techniques, т. 50, стр.1761-1771, Июль 2002; N. Joachimowicz, J.J. Mallorqui, J. Ch. Bolomey, и A. Brouguetas, « Convergence and stability assessment of Newton-Kantorovich reconstruction algorithms for microwave tomography», IEEE Trans. Medical Imaging, том 17, стр.562-570, Авг. 1998; P. Lobel, R. Kleinman, Ch. Pichot, L. Blanc-Fйraud, и M. Barlaud «Conjugate Gradient Method for Solving Inverse Scattering with Experimental Data», IEEE Antennas & Propagation Magazine, Том 38, стр. 48-51, Июнь. 1996; W.C. Chew, и Y.M. Wang «Reconstruction of two-dimensional permittivity distribution using the distorted Born iterative method», IEEE Trans. Medical Imaging, том 9, стр.218-225, Июнь 1990; P.M. Meaney, K.D. Paulsen, A. Hartov, и R.K. Crane «Microwave imaging for tissue assessment: Initial evaluation in multitarget tissue equivalent phantoms», IEEE Trans. Biomedical Engineering, том 43, стр. 878-890, Сент. 1996.
[0010] В предыдущих работах были разработаны математические алгоритмы и соответствующие им системы и программные реализации, которые оказались очень надежными и производят изображения объектов разного размера от нескольких сантиметров в вырезаемом собачьем сердце до полноразмерного тела в 2D, 3D и 3D-векторных случаях. Смотри A.E. Souvorov, A.E. Bulyshev, S.Y. Semenov, R.H. Svenson, A.G. Nazarov, Y.E. Sizov, и G.P. Tatsis «Microwave tomography: A two-dimensional Newton iterative scheme», IEEE Trans. Microwave Theory and Techniques, том. 46, стр. 1654-1659, Ноябрь. 1998; A.E. Bulyshev, A.E. Souvorov, S.Y. Semenov, R.H. Svenson, A.G. Nazarov, Y.E. Sizov, и G.P. Tatsis «Three-dimensional microwave tomography. Theory и computer experiments in scalar approximation», Inverse Problems, том. 16, стр. 863-875, Июнь 2000; Bulyshev A.E, Souvorov A.E., Semenov S.Y., Posukh V.G., Sizov Y.E. « Three-dimensional Vector Microwave 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 the Iteration Reconstruction Algorithm», IEEE Trans MTT, 1998, 46, 2, 133-141; Semenov S.Y., Bulyshev A.E., Abubakar A., Posukh V.G., Sizov Y.E., Souvorov A.E., Van den Berg P., Williams T. «Microwave tomographic imaging of the high dielectric contrast objects using different imaging approaches», IEEE Trans. MTT, т. 53, № 7, стр. 2284-2294, 2005; Semenov S.Y., Kellam J.F., Althausen P., Williams T.C., Abubakar A., Bulyshev A., Sizov Y. «Microwave tomographic imaging of the high dielectric contrast objects using different imaging approaches», 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 myocardial infarction. 1. Excised canine heartsʺ, Annals of Biomedical Engineering, 2003, 31, 262-270. Однако ни один из вышеперечисленных способов, приведенный в этих предыдущих работах, не оказался эффективным при визуализации внутри сильно экранированных объектов. Новый подход необходим для точного представления объектов при визуализации ЭМТ, имеющих высокий диэлектрический контрастный экран, таких как мозг человека, но не ограничиваясь им.
СУЩНОСТЬ НАСТОЯЩЕГО ИЗОБРЕТЕНИЯ
[0011] Некоторые примерные варианты осуществления настоящего изобретения могут преодолеть один или более из вышеупомянутых недостатков и других недостатков, не описанных выше, но настоящее изобретение не требуется для преодоления какого-либо особого недостатка, описанного выше, и некоторые примерные варианты осуществления настоящего изобретения могут не преодолевать какой-либо из недостатков, описанных выше.
[0012] В широком смысле настоящее изобретение в соответствии с одним из аспектов представляет собой способ использования электромагнитного поля так, что создается картина электромагнитной интерференции в области визуализации, выявляющая суперпозицию 3D диэлектрической структуры объекта вместе с картиной электромагнитной интерференции, в то время как картина электромагнитной интерференции была дополнительно распознана и применена к 3D электромагнитной суперпозиции, сводя на нет или уменьшая картину электромагнитной интерференции и выявляя 3D диэлектрическую структуру объекта.
[0013] В широком смысле настоящее изобретение в соответствии с другим аспектом представляет собой способ 4D динамической томографии с распознаванием объединенных электромагнитных картин.
[0014] В широком смысле настоящее изобретение в соответствии с другим аспектом представляет собой способ контроля жизнеспособности и функциональных состояний биологической ткани с использованием 4D динамической томографии с распознаванием объединенных электромагнитных картин.
[0015] В широком смысле настоящее изобретение в соответствии с другим аспектом представляет собой, как показано и описано, способ электромагнитной томографии с распознаванием картин интерференции (EMIPRT) для использования в системе реконструкции изображения.
[0016] В широком смысле настоящее изобретение в соответствии с другим аспектом представляет собой способ электромагнитной томографии с распознаванием картин интерференции (EMIPRT) для использования в системе реконструкции изображения, содержащий: генерацию данных электромагнитного поля, соответствующих объекту в области визуализации, посредством системы электромагнитной томографии, причем данные электромагнитного поля измеряют на множестве приемников после его создания на упомянутом множестве передатчиков и взаимодействия с объектом; и использование сгенерированных данных электромагнитного поля повторно, рекурсивным образом, с: формированием невозмущенного изображения электромагнитной интерференции, формированием возмущенного изображения электромагнитной интерференции, основанного по меньшей мере частично на невозмущенном изображении электромагнитной интерференции, распознаванием картин электромагнитной интерференции в повторно сформированных возмущенных изображениях электромагнитной интерференции и формированием суперпозиционного изображения путем сведения к нулю или уменьшения распознанных картин электромагнитной интерференции от возмущенного изображения электромагнитной интерференции.
[0017] В качестве признака этого аспекта этап формирования возмущенного изображения электромагнитной интерференции, основанного по меньшей мере частично на невозмущенном изображении электромагнитной интерференции, включает в себя формирование возмущенного изображения электромагнитной интерференции, основанного по меньшей мере частично на определении фактора объекта, который является функционалом разниц между экспериментальными электромагнитными полями и электромагнитными полями, вычисленными во время этапа формирования невозмущенного изображения электромагнитной интерференции. Более конкретно, фактор объекта определяется как
Figure 00000001
, где
Figure 00000002
- экспериментально смоделированное или измеренное значение соответственно z-компоненты электромагнитного поля, измеренного приемником j, когда передатчик i является источником электромагнитного поля, где
Figure 00000003
представлено в общем виде как
Figure 00000004
, где α, β и γ являются коэффициентами действительных ненулевых или нулевых значений, где Ω - оператор регуляризации, и где
Figure 00000005
является функцией ее аргумента; фактор объекта определяется как
Figure 00000006
, где
Figure 00000002
- экспериментально смоделированное или измеренное значение соответственно z-компоненты электромагнитного поля, измеренного приемником j, когда передатчик i является источником электромагнитного поля, и где
Figure 00000007
- максимальная норма экспериментально измеренной z-компоненты электромагнитного поля; фактор объекта определяется как
Figure 00000006
, где
Figure 00000002
- экспериментально смоделированное или измеренное значение соответственно z-компоненты электромагнитного поля, измеренного приемником j, когда передатчик i является источником электромагнитного поля, и где
Figure 00000008
является нормой экспериментально измеренной z-компоненты электромагнитного поля, измеренной приемником j, когда передатчик i является источником электромагнитного поля мощности θ; этап формирования возмущенного изображения электромагнитной интерференции включает в себя расчет
Figure 00000009
, где
Figure 00000010
и
Figure 00000011
представляют собой 3D распределения электромагнитных полей (x, y, z) от электромагнитных источников частоты f k , расположенных в местоположениях физических источников (от 1 до N) в системе электромагнитной томографии и в местоположении физических приемников (от 1 до M) соответственно, взятые как сопряженные значения, и где
Figure 00000012
- «ij»-ая составляющая фактора объекта от передатчика i к приемнику j; этап распознавания картин электромагнитной интерференции в повторно сформированных возмущенных изображений электромагнитной интерференции включает в себя вычисление сумм
Figure 00000009
, где
Figure 00000010
и
Figure 00000011
представляют собой 3D распределения электромагнитных полей (x, y, z) от электромагнитных источников частоты f k , расположенных в местоположениях физических источников (от 1 до N) в системе электромагнитной томографии и в местоположении физических приемников (от 1 до M) соответственно, взятых как сопряженные значения, и где
Figure 00000012
- «ij»-ая составляющая фактора объекта от передатчика i к приемнику j; и/или где этап распознавания картин электромагнитной интерференции в повторно сформированных возмущенных изображениях электромагнитной интерференции включает в себя вычисление, для итерации i> 1,
Figure 00000013
где для простоты частотные члены опущены, где
Figure 00000014
и
Figure 00000015
представляют собой 3D распределения электромагнитных полей (x,y,z) от электромагнитных источников, расположенных в местоположениях физических источников (от 1 до N) в системе электромагнитной томографии и в местоположениях физических приемников (от 1 до M) соответственно, взятых как сопряженные значения, и где
Figure 00000012
- «ij»-ая составляющая фактора объекта от передатчика i к приемнику j.
[0018] В другом признаке этого аспекта способ дополнительно включает в себя этап, выполняемый после каждого повторного этапа формирования суперпозиционного изображения, определения того, была ли достигнута целевая сходимость.
[0019] В другом признаке этого аспекта этот способ используется как часть способа генерации 4D дифференциальных (динамических) объединенных изображений. Более конкретно, генерация 4D дифференциальных (динамических) объединенных изображений включает в себя объединение по меньшей мере одного последовательно сформированного изображения, указывающего на относительное физиологическое изменение, с базовым анатомическим изображением для отображения в виде единого слитного изображения; и/или способ генерации 4D дифференциальных (динамических) объединенных изображений используется как часть способа контроля жизнеспособности и/или функциональных состояний биологической ткани, используя 4D динамическую электромагнитную томографию с распознаванием объединенных картин.
[0020] В другом признаке этого аспекта этапы формирования невозмущенного изображения электромагнитной интерференции, формирования возмущенного изображения электромагнитной интерференции, основанного по меньшей мере частично на невозмущенном изображении электромагнитной интерференции, распознавания картин электромагнитной интерференции в повторно сформированных возмущенных изображений электромагнитной интерференции и формирования суперпозиционного изображения путем сведения к нулю или уменьшения распознанных картин электромагнитной интерференции от возмущенного изображения электромагнитной интерференции выполняются последовательно.
[0021] В другом признаке этого аспекта, способ дополнительно включает в себя этап отображения суперпозиционного изображения посредством блока отображения.
[0022] В широком смысле настоящее изобретение в соответствии с другим аспектом представляет собой способ реконструкции изображения с использованием электромагнитной томографии с распознаванием картин интерференции, содержащий: генерацию данных электромагнитного поля, соответствующих объекту в области визуализации, посредством системы электромагнитной томографии, причем данные электромагнитного поля измеряют на множестве приемников после его создания на множестве передатчиков и взаимодействия с объектом; и использование сгенерированных данных электромагнитного поля повторно, рекурсивным образом, с: формированием невозмущенного изображения электромагнитной интерференции, формированием возмущенного изображения электромагнитной интерференции, основанного по меньшей мере частично на невозмущенном изображении электромагнитной интерференции, распознаванием картин электромагнитной интерференции в повторно сформированных возмущенных изображениях электромагнитной интерференции и формированием суперпозиционное изображения путем сведения к нулю или уменьшения распознанных картин электромагнитной интерференции от возмущенного изображения электромагнитной интерференции.
[0023] В качестве признака этого аспекта этап формирования возмущенного изображения электромагнитной интерференции, основанного по меньшей мере частично на невозмущенном изображении электромагнитной интерференции, включает в себя формирование возмущенного изображения электромагнитной интерференции, основанного по меньшей мере частично на определении фактора объекта, который является функционалом разниц между экспериментальными электромагнитными полями и электромагнитными полями, вычисленными во время этапа формирования невозмущенного изображения электромагнитной интерференции.
[0024] В другом признаке этого аспекта способ дополнительно включает в себя этап, выполняемый после каждого повторного этапа формирования суперпозиционного изображения, определения того, была ли достигнута целевая сходимость.
[0025] В другом признаке этого аспекта этот способ используется как часть способа генерации 4D дифференциальных (динамических) объединенных изображений. Более конкретно, генерация 4D дифференциальных (динамических) объединенных изображений включает в себя объединение по меньшей мере одного последовательно сформированного изображения, указывающего на относительное физиологическое изменение, с базовым анатомическим изображением для отображения в виде единого слитного изображения; и/или способ генерации 4D дифференциальных (динамических) объединенных изображений используется как часть способа контроля жизнеспособности и/или функциональных состояний биологической ткани, используя четырехмерную динамическую электромагнитную томографию с распознаванием объединенных картин.
[0026] В другом признаке этого аспекта этапы формирования невозмущенного изображения электромагнитной интерференции, формирования возмущенного изображения электромагнитной интерференции, основанного по меньшей мере частично на невозмущенном изображении электромагнитной интерференции, распознавания картин электромагнитной интерференции в повторно сформированных возмущенных изображений электромагнитной интерференции и формирования суперпозиционного изображения путем сведения к нулю или уменьшения распознанных картин электромагнитной интерференции от возмущенного изображения электромагнитной интерференции выполняются последовательно.
[0027] В другом признаке этого аспекта, способ дополнительно включает в себя этап отображения суперпозиционного изображения посредством блока отображения.
[0028] В широком смысле настоящее изобретение в соответствии с другим аспектом представляет собой систему реконструкции изображений с использованием электромагнитной томографии с распознаванием картин интерференции, содержащую: систему электромагнитной томографии, которая генерирует данные электромагнитного поля, соответствующие объекту в области визуализации, причем система электромагнитной томографии содержит множество передатчиков электромагнитных волн, множество приемников, которые измеряют электромагнитные данные после создания на упомянутом множестве передатчиков и взаимодействия с объектом, и пограничное устройство; и центр обработки, который используя сгенерированные данные электромагнитного поля повторно рекурсивным образом: выполняет этапы формирования невозмущенного изображения электромагнитной интерференции, формирования возмущенного изображения электромагнитной интерференции, основанного по меньшей мере частично на невозмущенном изображении электромагнитной интерференции, распознавания картин электромагнитной интерференции в повторно сформированных возмущенных изображениях электромагнитной интерференции и формирования суперпозиционного изображения путем сведения к нулю или уменьшения распознанных картин электромагнитной интерференции от возмущенного изображения электромагнитной интерференции.
[0029] В качестве признака этого аспекта этап формирования возмущенного изображения электромагнитной интерференции, основанного по меньшей мере частично на невозмущенном изображении электромагнитной интерференции, включает в себя формирование возмущенного изображения электромагнитной интерференции, основанного по меньшей мере частично на определении фактора объекта, который является функционалом разниц между экспериментальными электромагнитными полями и электромагнитными полями, вычисленными во время этапа формирования невозмущенного изображения электромагнитной интерференции.
[0030] В другом признаке этого аспекта центр обработки дополнительно выполняет этап, выполняемый после каждого повторного этапа формирования суперпозиционного изображения, определения того, была ли достигнута целевая сходимость.
[0031] В другом признаке этого аспекта этапы, выполняемые центром обработки, используются в качестве части способа генерации 4D дифференциальных (динамических) объединенных изображений. Более конкретно, генерация 4D дифференциальных (динамических) объединенных изображений включает в себя объединение по меньшей мере одного последовательно сформированного изображения, указывающего на относительное физиологическое изменение, с базовым анатомическим изображением для отображения в виде единого слитного изображения; и/или способ генерации 4D дифференциальных (динамических) объединенных изображений используется как часть способа контроля жизнеспособности и/или функциональных состояний биологической ткани, используя четырехмерную динамическую электромагнитную томографию с распознаванием объединенных картин.
[0032] В другом признаке этого аспекта этапы формирования невозмущенного изображения электромагнитной интерференции, формирования возмущенного изображения электромагнитной интерференции, основанного по меньшей мере частично на невозмущенном изображении электромагнитной интерференции, распознавания картин электромагнитной интерференции в повторно сформированных возмущенных изображений электромагнитной интерференции и формирования суперпозиционного изображения путем сведения к нулю или уменьшения распознанных картин электромагнитной интерференции от возмущенного изображения электромагнитной интерференции выполняются последовательно.
[0033] В другом признаке этого аспекта способ дополнительно включает в себя блок отображения, который отображает суперпозиционное изображение.
[0034] Прочие области возможного применения данного изобретения станут очевидными после ознакомления с подробным описанием, приведенным ниже. Подробное описание и конкретные примеры дают понимание наиболее предпочтительных областей использования изобретения, но приводятся исключительно в качестве иллюстрации и не призваны ограничивать объем изобретения.
КРАТКОЕ ОПИСАНИЕ ЧЕРТЕЖЕЙ
[0035] Дополнительные признаки, варианты осуществления и преимущества настоящего изобретения станут очевидными из следующего подробного описания со ссылкой на чертежи, на которых:
на ФИГ. 1 представлена упрощенная схематическая иллюстрация частей системы электромагнитной томографии (ЭМТ);
на ФИГ. 2 представлена блок-схема способа электромагнитной томографии с распознаванием картин интерференции (EMIPRT) в соответствии с одним или более предпочтительными вариантами осуществления настоящего изобретения;
на ФИГ. 3 представлено примерное невозмущенное изображение ЭМ интерференции для области визуализации, предназначенной для визуализации головы человека;
на ФИГ. 4А и 4В являются примерными возмущенными изображениями ЭМ интерференции головы человека для двух случаев;
на ФИГ. 5 представлена блок-схема одной стратегии распознавания картин и набор уравнений для блока распознавания картин, показанного на Фиг. 2;
на ФИГ. 7 представлена блок-схема одной возможной стратегии коррекции картины, которая может быть выполнена в блоке формирования суперпозиционного изображения для каждой из множества подобластей;
на ФИГ. 8A представлено примерное выходное изображение головы жертвы возможного инсульта из итерационного способа визуализации, аналогичного способу, показанному на ФИГ. 2, но где блок распознавания картин и его выходная информация не используются во время визуализации; и
на ФИГ. 8B представлено примерное выходное изображение той же головы жертвы возможного инсульта, что и на ФИГ. 8А, но где блок распознавания картин и его выходная информация используются блоком во время формирования суперпозиционного изображения.
ПОДРОБНОЕ ОПИСАНИЕ
[0036] Следует отметить, что любой специалист в данной области, имеющий общее представление о предмете вопроса, легко поймет то, что данное изобретение имеет большой спектр применения и потенциал. Кроме того, любой «предпочтительный» вариант рассматривается как часть наиболее оптимального способа использования, предусмотренного как инструмент для применения настоящего изобретения. Также в качестве дополнительных иллюстраций могут рассматриваться и другие варианты осуществления, что позволит предоставить полную и полезную с практической точки зрения характеристику изобретения. Следует понимать, что любой вариант осуществления может содержать только один или более вышеописанных аспектов изобретения и может дополнительно содержать только один или более вышеописанных признаков. Кроме того, многие варианты осуществления, такие как адаптации, вариации, модификации и аналогичные им изменения будут косвенно раскрыты в составе рассмотренных в настоящем документе вариантов осуществления, и они также входят в рамки настоящего изобретения.
[0037] Соответственно, хотя настоящее изобретение подробного описывается в данном документе в виде одного или более вариантов осуществления, понимается, что данная информация является иллюстративной, представляет собой примеры использования изобретения и приводится только для того, чтобы предоставить полное и полезное с практической точки зрения описание изобретения. Подробное раскрытие в настоящем документе информации в отношении одного или более вариантов осуществления не призвано и не может толковаться как ограничивающее объем патентной защиты данного изобретения, объем такой защиты определяется патентными формулами и с учетом их эквивалентов. Не предполагается, что объем патентной защиты будет определяться с ограничением какой-либо формулы положениями настоящего документа, которые не вошли непосредственно в саму патентную формулу.
[0038] Таким образом, например, любая очередность и/или временная последовательность этапов в составе различных процессов или способов, которые описаны в настоящем документе, является иллюстративной и не носит характер ограничения. Соответственно, следует исходить из того, что хотя этапы в рамках различных процессов или способов могут быть показаны и описаны в качестве некоторой последовательности или временного порядка, эти этапы в рамках такого процесса или способов не ограничиваются совершением их в какой-либо конкретной последовательности или порядке, если не указано иное. На самом деле, этапы в составе таких процессов или способов могут в целом осуществляться в различной очередности или последовательности, но все равно являться частью настоящего изобретения. Соответственно предполагается, что объем патентной защиты данного изобретения определяется по прилагаемым патентным формулам, а не по описанию, приведенному в настоящем документе.
[0039] Кроме того, важно отметить, что каждый термин, используемый в настоящем документе, должен быть понят специалистом в данной области, как значение такого термина в контексте его использования. Если значение какого-либо термина, используемого в настоящем документе (и понятого специалистом, с учетом его контекста), отличается каким-либо образом от конкретного словарного определения такого термина, то исходят из того, что приоритетный характер должно иметь значение, в котором этот термин понятен специалисту.
[0040] Что касается применимости 35 U.S.C. §112, 6, элемент пункта формулы изобретения не должен рассматриваться в соответствии с настоящим предусмотренным законом положением, если только явная фраза «означает для» или «этап для» фактически не используется в таком элементе пункта формулы изобретения, после чего это нормативное положение предназначено для применения в интерпретации такого элемента пункта формулы изобретения.
[0041] Кроме того, важно отметить, что при использовании в настоящем документе «какой-либо» и «какая-либо» в целом означают «по меньшей мере один/одна», и не исключают наличия множества таких объектов, если контекстом не продиктовано иное. Соответственно, ссылка на «корзину для пикника с яблоком внутри» является описанием «корзины для пикника с по меньшей мере одним яблоком внутри», а также «корзины для пикника с яблоками внутри». И наоборот, ссылка на «корзину для пикника с одним яблоком внутри» является описанием «корзины для пикника с только одним яблоком внутри».
[0042] При использовании в настоящем документе в качестве союза для объединения перечисляемых позиций «или» означает «по меньшей мере один из предметов», но не исключает допустимости наличия нескольких предметов из списка. Таким образом, ссылка на «корзину для пикника с сыром или крекерами внутри» является описанием «корзины для пикника с сыром, но без крекеров внутри», «корзины для пикника с крекерами, но без сыра внутри», и «корзины для пикника и с сыром, и с крекерами внутри». И наконец, при использовании для объединения перечисляемых позиций «и» означает «все предметы из списка». Таким образом, ссылка на «корзину для пикника с сыром и крекерами внутри» является описанием «корзины для пикника с сыром внутри, а также корзины для пикника, внутри которой есть еще и крекеры», и описанием «корзины для пикника с крекерами внутри, а также корзины для пикника, внутри которой есть еще и сыр».
[0043] Ниже рассматриваются схемы, на которых одинаковыми номерами обозначены одни и те же компоненты на различных изображениях, и описывается один или более вариантов осуществления данного изобретения. Приводимое далее описание одного или более предпочтительных вариантов осуществления является только примером по своему характеру и ни при каких обстоятельствах не должно ограничивать характер данного изобретения, его применение или формы использования.
[0044] На ФИГ. 1 представлена упрощенная схематическая иллюстрация частей системы 10 электромагнитной томографии (ЭМТ). В системе 10 пограничное устройство 12 окружает область 21 визуализации. Объект 19 (в этом примере, голова человека) помещается в область 21 визуализации. Множество ЭМ аппаратных устройств 30, обычно, но не обязательно всегда расположенных на пограничном устройстве 12, действуют как передатчики (источники) и/или приемники. (Как описано здесь, считается, что аппаратные устройства 30 обычно действуют и как передатчики, и как приемники, но будет понятно, что передатчики и приемники, которые отделены друг от друга, можно использовать дополнительно или альтернативно). Компьютерная система 28 служит в качестве центра обработки, где предварительно обрабатываются исходные данные, а изображения реконструируются и предварительно обрабатываются. По меньшей мере в некоторых вариантах осуществления ЭМТ система 10 дополнительно включает в себя фоновую или соответствующую текучую среду или другой носитель. Соответствующая среда представляет собой раствор или гель, который необходим или полезен внутри области визуализации, когда объект 19 визуализируется, для борьбы с электромагнитными проблемами согласования с телом и/или другими сложностями. По меньшей мере в некоторых вариантах осуществления компьютерная система 28 и ее функциональные средства обработки данных и программное обеспечение для визуализации непосредственно соединены с аппаратным устройствам 30 для передачи/приема ЭМ-полей, в то время как в других вариантах осуществления некоторые или все из компьютерной системы 28 дистанционно соединены посредством беспроводной технологии и/или высокоскоростных проводных соединений. Функционально большая часть работы ЭМТ системы 10 может быть аналогична работе, описанной в вышеупомянутом патенте США № 9414749, но здесь могут быть описаны различные конкретные варианты осуществления и признаки.
[0045] Как описано выше, ЭМТ визуализация объектов с высоким диэлектрическим контрастом, в том числе биологических объектов, включает в себя очень сложную проблему так называемой «дифракционной томографии». Высокий диэлектрический контраст между тканями с высоким содержанием воды, такими как, мышечная ткань, но не ограничиваясь ею, и с низким содержанием воды, например, кость, но не ограничиваясь ею, является дополнительным осложнением при использовании ЭМ полей для визуализации. На ФИГ. 2 представлена блок-схема способа 100 электромагнитной томографии с распознаванием картин интерференции (EMIPRT) в соответствии с одним или более предпочтительными вариантами осуществления настоящего изобретения. Способ 100 EMIPRT может быть выполнен с использованием данных ЭМ-полей, сгенерированных с использованием некоторых или всех элементов системы 10 ЭМТ, кратко описанных в отношении ФИГ.1.
[0046] Способ 100 EMIPRT является итерационным процессом, где проверка сходимости (показанная на этапе 180) происходит после каждой итерации по различным процессам визуализации до тех пор, пока не будут получены соответствующие результаты и предоставлены в качестве вывода 200 изображения. Первичными входными данными для способа 100 EMIPRT являются данные 110 области визуализации и экспериментально измеренные ЭМ данные 120. Данные 110 области визуализации включают в себя пространственное расположение источников ЭМ полевых аппаратных устройств 30 на множестве пространственных местоположений и с k ми частотами (k от 1 до K) (
Figure 00000016
и приемников ЭМ полевых аппаратных устройств 30 на множестве пространственных местоположений и с k ми частотами (k от 1 до K) (
Figure 00000016
по отношению к области 21 визуализации, где N - число источников ЭМ полевых устройств 30, а M - число приемников ЭМ полевых устройств 30. Данные 110 области визуализации также включают в себя диэлектрические свойства соответствующей среды (
Figure 00000017
) внутри области 21 визуализации, и также могут включать в себя другую физическую или относящуюся к окружающей среде информацию/данные. Экспериментально измеренные ЭМ данные 120 представляют собой матрицу, содержащую экспериментально измеренные ЭМ поля (например, как в форме амплитуды, фазы и поляризации, так и в форме реальной, мнимой и поляризационной форме) на точках приемников, представленных как
Figure 00000018
. Точки приемников ЭМ поля могут быть связаны с любым или всеми ЭМ полевыми устройствами 30.
[0047] В блоке 130 формируется «невозмущенное» изображение ЭМ интерференции. В начальном проходе формирование невозмущенного изображения ЭМ интерференции (показанного как вывод 135 блока 130 формирования невозмущенного изображения) получается только на основе данных 110 области визуализации способа 100 EMIPRT. В последующих итерациях формирование невозмущенного изображения интерференции предпочтительно также использует пространственное распределение диэлектрических свойств в области 21 визуализации. В зависимости от приложения пространственное распределение может быть 2D пространственным распределением или 3D пространственным распределением. По меньшей мере в некоторых вариантах осуществления 3D пространственное распределение диэлектрических свойств описывается как ɛ*(x,y,z) в декартовой системе координат и получается во время предыдущей итерации в блоке 170 формирования суперпозиционного изображения, описанном ниже. 3D пространственное распределение диэлектрических свойств ɛ*(x,y,z) можно охарактеризовать как изображение исследуемого объекта 19 по диэлектрической шкале.
[0048] Задача блока 130 формирования невозмущенного изображения ЭМ интерференции заключается в синтезе 2D изображения или 3D изображения (картины) интерференции 135 от упомянутого множества ЭМ источников/приемников 30, расположенных внутри или на границе области 21 визуализации. Во время процесса визуализации все ЭМ полевые устройства 30 (i=1 до N и j=1 до M) рассматриваются как источники ЭМ поля. Определение (3D) невозмущенного изображения 135 ЭМ интерференции генерируется из следующей функции:
(1)
Figure 00000019
,
где
Figure 00000020
и
Figure 00000021
представляют собой 3D распределение ЭМ полей (x,y,z) в области 21 визуализации от ЭМ полевых устройств 30 частоты f k зондирования, расположенных в местоположении физических источников (от 1 до N) и в местоположении физических приемников (от 1 до М) соответственно и где W k - это частотная весовая функция, которая учитывает различный вклад частотно-зависимых изображений ЭМ интерференции в функцию (1). На ФИГ. 3 представлено примерное невозмущенное изображение 135 ЭМ интерференции для области 21 визуализации, предназначенной для визуализации головы человека. Изображение 135 было получено на определенной частоте, которая в этом случае составляла 1 ГГц. Изображение представляет собой двумерное поперечное сечение в плоскости X-Y, где шкалы X и Y измеряются в сантиметрах.
[0049] Обратимся к ФИГ. 2, следующей важной операцией в способе 100 EMIPRT является определение фактора объекта, как показано на этапе 145, для применения к изображению ЭМ интерференции. Фактор объекта предпочтительно является функционалом различий между измеренным (из блока 120) и имитированным (из блока 130) ЭМ полями для каждой пары от передающих (i=от 1 до N) к принимающим (j=от 1 до M) устройств 30 при некоторой частоте зондирования. Таким образом, первичные вводы в блок 145 определения фактора объекта являются вторыми из двух основных вводов в способ 100 (экспериментально измеренные ЭМ поля в точках 120 приемников, описанных ранее) и входными данными 140 (рассчитанными или смоделированными ЭМ полями в точках 120 приемников), генерируемыми из блока 130 формирования невозмущенных изображений ЭМ интерференции. Эти два основных ввода в блок 145 определения фактора объекта являются матрицами данных. Первая матрица содержит экспериментально измеренные ЭМ поля, представленные как
Figure 00000018
, а вторая матрица содержит моделируемые ЭМ поля, представленные как
Figure 00000022
, где в обоих случаях
Figure 00000023
представляет собой комплексное векторное значение, представляющее электрические (E) и/или магнитные (H) компоненты ЭМ поля. Эти значения могут быть представлены как комплексное (действительная и мнимая части) значение каждой векторной компоненты E и/или H поля или скалярное значение E и/или H поля. Они также могут быть представлены как амплитуда и/или фаза каждой векторной компоненты E и/или H поля или скалярного значения E и/или H поля. По меньшей мере в одном варианте осуществления настоящего изобретения фактор объекта может быть представлен, вообще говоря, следующим образом:
(2)
Figure 00000001
,
где
Figure 00000002
представляет собой экспериментально смоделированное или измеренное значение соответственно z-компоненты ЭМ поля, измеренного приемником j, когда передатчик i является источником ЭМ поля, где
Figure 00000003
представляется в общем виде как:
(3)
Figure 00000024
,
где α, β и γ - коэффициенты вещественных ненулевых или нулевых значений, Ω - оператор регуляризации, а
Figure 00000025
- функция ее аргумента.
[0050] В других вариантах осуществления настоящего изобретения фактор объекта может иметь другой состав. Это будет влиять на изображение ЭМ интерференции и его картину. Один пример альтернативного фактора объекта представлен ниже для α=1, β=0 и γ=0:
(4)
Figure 00000006
,
где
Figure 00000002
представляет собой экспериментально смоделированное или измеренное значение соответственно z-компоненты ЭМ поля, измеренного приемником j, когда передатчик i является источником ЭМ поля, и где
Figure 00000007
- максимальная норма экспериментально измеренной z-компоненты ЭМ поля. Другой пример альтернативного фактора 145 объекта представлен ниже для α=1, β=0 и γ=0:
(5)
Figure 00000026
,
где
Figure 00000002
представляет собой экспериментально смоделированное или измеренное значение соответственно z-компоненты ЭМ поля, измеренного приемником j, когда передатчик i является источником ЭМ поля, и где
Figure 00000008
- норма экспериментально измеренной z-компоненты ЭМ поля, измеренного приемником j, когда передатчик i является источником электромагнитного поля мощности θ. Мощность θ может, например, иметь значение 1 или 2, или 3, или 1/2 или даже какое-то другое значение.
[0051] Вывод из блока 145 определения фактора объекта и вывод из блока 130 формирования невозмущенного изображения ЭМ интерференции (то есть 2D или 3D невозмущенного изображения 135 интерференции) используются для формирования возмущенного изображения ЭМ интерференции, как показано на блоке 150. Изображение, сгенерированное в блоке 150, представляет собой 2D или 3D изображение в зависимости от приложения. Формирование (3D) возмущенного изображения ЭМ интерференции в области 21 визуализации, как показано в блоке 150, включает в себя расчет сумм, как показано, например, в уравнении (6):
(6)
Figure 00000009
,
где
Figure 00000027
и
Figure 00000027
- 3D распределение ЭМ полей (x,y,z) от ЭМ источников 30 частоты f k , расположенных в местоположениях физических источников (от 1 до N) и в местоположении физических приемников (от 1 до M) соответственно, взятых как сопряженные значения, и
Figure 00000012
- это «ij»-ая составляющая выбранного фактора объекта от передатчика i к приемнику j, как определено в блоке 145.
[0052] Вывод блока 150 формирования возмущенного изображения ЭМ интерференции представляет собой возмущенное изображение 155 ЭМ интерференции (или данные, его представляющие). В этом отношении ФИГ. 4A и 4B представляют собой примерные возмущенные изображения 155 ЭМ интерференции головы человека для двух случаев. Каждое изображение 155 представляет собой двумерное поперечное сечение в плоскости X-Y, где шкалы X и Y измеряются в сантиметрах. (Понятно, что, хотя ФИГ. 4A и 4B представляют собой 2D иллюстрации, 3D изображения могут быть получены путем дальнейшего применения того же процесса.) Альтернативный фактор объекта, показанный на (4), использовался для этих двух примеров. При сравнении невозмущенного ЭМ изображения 135 с ФИГ. 3 и возмущенных ЭМ изображений 155 с ФИГ. 4А и 4В суперпозиция диэлектрической структуры объекта 410, 420 (поперечное сечение X-Y головы человека) вместе с картиной 415, 425 ЭМ интерференции хорошо видна на ФИГ. 4A и 4B, где наружное темное кольцо на каждой иллюстрации представляет собой человеческий череп, а топологически подобные картины внутри такого темного кольца являются остальными картинами ЭМ интерференции.
[0053] После формирования возмущенного изображения ЭМ интерференции в блоке 150 изображение (данные) 155 отправляется как в блок 160 распознавания картин, так и в блок 170 формирования суперпозиционного изображения. Распознавание картин интерференции и его применения выполняется как двухэтапный процесс, где этап 1 - само распознавание картин, а этап 2 представляет собой применение распознанной картины к 2D или 3D суперпозиционному ЭМ изображению, сводя на нет или уменьшая картину ЭМ интерференции и выявляя 2D или 3D диэлектрическую структуру объекта 19. Эти два этапа представлены на ФИГ. 2 блоком 160 распознавания картин и блоком 170 формирования суперпозиционного изображения.
[0054] Могут использоваться различные подходы к распознаванию картин. Например, на ФИГ. 5 представлена блок-схема одной стратегии 500 распознавания картин и набор уравнений для блока 160 распознавания картин, показанного на ФИГ. 2. На ФИГ. 5 используется 2D-пример, но будет понятно, что такое распознавание картин может быть применено к множеству поперечных сечений X-Y для создания 3D изображения. Что показано на ФИГ. 5.
[0055] Альтернативная стратегия распознавания картин, иногда называемая далее «2-ой T оптимизацией», преимущественно позволяет вырабатывать гораздо более резкие картины на начальной или первой итерации блока 170 формирования суперпозиционного изображения и, как следствие, упрощает распознавание картин на последующих итерациях. Ключевой особенностью 2-ой T оптимизации является использование обобщения (6) на итерации i> 1 для формирования суперпозиционного изображения (где для простоты частотные члены опущены):
(7)
Figure 00000028
В некоторых вариантах осуществления коэффициент α может быть независимым от номера итерации (и может быть найден методом проб), в то время как в других вариантах осуществления он может иметь зависимое от итерации значение, которое помогает ускорить процесс сходимости. В последнем случае одной из возможных стратегий является наличие значения коэффициента α на итерации i+1 как максимальной абсолютной разности между изображением на итерации i и изображением на итерации i - 1:
(8)
Figure 00000029
.
В других вариантах осуществления настоящего изобретения для поиска оптимального коэффициента α могут использоваться другие стратегии и/или алгоритм распознавания картин.
[0056] В некоторых вариантах осуществления одна и та же стратегия/процедура распознавания картин может выполняться на каждой итерации способа 100, тогда как в других вариантах осуществления стратегия/процедура распознавания картин может динамически обновляться после каждой последующей итерации.
[0057] Вывод из блока 160 распознавания картин и вывод из блока 150 формирования возмущенного изображения ЭМ интерференции (то есть возмущенного изображения 155 ЭМ интерференции) выдаются в блок 170 формирования суперпозиционного изображения. Формирование 170 суперпозиционного изображения может быть выполнено путем распознавания картин посредством блока 160 путем сравнения их с возмущенными изображениями 155 ЭМ интерференции, а затем выполнением корректировки картин. В блоке 170 данные анализируются для уменьшения или сведения к нулю картины ЭМ интерференции и выявления истинной 2D или 3D структуры объекта 19, представленного в качестве вывода 200 изображения.
[0058] Коррекция картины может быть выполнена следующим образом. Во-первых, по меньшей мере в некоторых вариантах осуществления способа 100 вся область 21 визуализации может подразделяться на ряд характерных подобластей визуализации. В этом отношении ФИГ. 6 - изображение 600, иллюстрирующее разбиение области 21 визуализации на подобласти. Более конкретно, используя пример головы человека в качестве объекта 19, область 21 визуализации может быть подразделена на наружную (внешнюю) область 605, область 610 черепа, область 615 цереброспинальной жидкости (CSF), область 620 серого вещества и область 625 белого вещества. На ФИГ. 7 представлена блок-схема одной возможной стратегии 700 коррекции картины, которая может быть выполнена в блоке 170 формирования суперпозиционного изображения для каждой из множества подобластей. На этапе 705 вычисляется нормализованное V {значение (6) или (7)}, где V - объем конкретной подобласти. На этапе 710, вычисляется нормализованное L - pattern {значение (6) или (7)}, где L - длина «ip»ой картины. Логика этапа 710 может быть легко понята при обращении к ФИГ. 3. Понятно, что имеются четкие круговые картины. Число этих круговых картин Np равно числу ячеек сетки по радиусу. Таким образом, для простых круговых картин, как на ФИГ. 3, интеграция на этапе 710 проводится по каждой картине (ʺipʺая картине) от 1 до Np по длине
Figure 00000030
. Аналогично, для более сложных картин объектов (например, нескольких кольцевых картин, представленных на ФИГ. 4A и 4B), этап 710 включает в себя вычисление нормализованного L - pattern {значения (6) или (7)}, где L - длина «ip»ой картины. Затем на этапе 715 изображение по ip ой картине корректируется путем установки нового значения (6) или (7), равного сумме старого значения и { V {(6) или (7)} - L -картина {(6) или (7)}. Этапы 710 и 715 повторяются для каждой ip ой картины.
[0059] Понятно, что стратегия коррекции картин, показанная на ФИГ. 7 - не единственная стратегия коррекции картин, которая может быть использована; другие стратегии коррекции картин могут дополнительно или альтернативно использоваться без отклонения от объема настоящего изобретения.
[0060] После каждой итерации посредством блока 170 результат оценивается «на предмет сходимости» на этапе 180. При достижении сходимости формируется итоговое выходное изображение 200.
[0061] Важное преимущество внедрения представленной здесь технологии распознавания картин можно понять следующим образом. По мере того, как продвигается итеративный процесс, представленный на ФИГ. 2, влияние картин ЭМ интерференции имеет тенденцию существенно искажать результирующие изображения целевой 2D или 3D диэлектрической структуры. Это связано с усилением такой картины с помощью аналогичных картин ЭМ интерференции, полученных в следующих итерациях, и в результирующего искажения в изображениях, произведенных ими. В этом отношении ФИГ. 8A является примерным выходным изображением головы 19 жертвы возможного инсульта из способа итеративной визуализации, аналогичного способу 100 на ФИГ. 2, но где блок 160 распознавания картин и его выходная информация не используются во время визуализации, тогда как на ФИГ. 8B - примерное выходное изображение 200 той же головы 19 жертвы возможного инсульта, что показано на ФИГ. 8А, но где блок 160 распознавания картин и его выходная информация используются в условиях блока 170 формирования суперпозиционного изображения. ФИГ. 8B ясно показывает область возможного инсульта 820 в правом нижнем углу изображения. Изображение на ФИГ. 8А, полученное без блока 160 распознавания картин интерференции и его применения в блоке 170, обнаруживается некоторое искажение изображения в области возможного инсульта 810, но без большей ясности можно только предположить, что искажение указывает область инсульта.
[0062] Также будет понятно, что, хотя описанная здесь методология (в том числе вариации и перестановки), была описана главным образом с точки зрения 2D и 3D изображений, считается, что она одинаково применима к технологии 4D ЭМТ (3D в пространстве плюс 1D во времени), где 4D динамические объединенные электромагнитные томографические изображения генерируются, как описано, например, в патентной заявке США 13/173 078, которая включена во всей полноте в настоящий документ по ссылке.
[0063] Исходя из изложенной выше информации, лицам, которые являются профессионалами в данной области, будет без труда понятно, что данное изобретение может широко применяться. Многочисленные варианты осуществления и адаптации данного изобретения, помимо тех, которые описаны в настоящей заявке, а также многочисленные вариации, модификации и эквивалентные действия, будут очевидны или разумно понятны на основе данного изобретения и приведенных выше его описаний, без отклонения от существа и объема данного изобретения.
[0064] Соответственно, несмотря на то, что данное изобретение было детально описано в настоящей заявке применительно к одному или более предпочтительным вариантам осуществления, следует понимать, что предоставленная информация носит исключительно иллюстративный характер и призвана дать представление о данном изобретении и приведена исключительно с целью предоставления полной и полезной с практической точки зрения информации об изобретении. Приведенная выше информация не должна толковаться как ограничивающая данное изобретение или иным образом исключающая любые другие варианты осуществления, адаптации, вариации, модификации или эквивалентные действия; причем настоящее изобретение ограничивается только прилагаемыми патентными формулами изобретения их эквивалентами.

Claims (57)

1. Способ электромагнитной томографии с распознаванием картин интерференции (EMIPRT) для использования в системе реконструкции изображения, содержащий:
генерацию данных электромагнитного поля, соответствующих объекту в области визуализации, посредством системы электромагнитной томографии, причем данные электромагнитного поля измеряют на множестве приемников после его создания на множестве передатчиков и взаимодействия с объектом; и
использование сгенерированных данных электромагнитного поля повторно, рекурсивным образом, с:
формированием невозмущенного изображения электромагнитной интерференции,
формированием возмущенного изображения электромагнитной интерференции, основанного по меньшей мере частично на невозмущенном изображении электромагнитной интерференции,
распознаванием картин электромагнитной интерференции в повторно сформированных возмущенных изображениях электромагнитной интерференции и
формированием суперпозиционного изображения путем сведения к нулю или уменьшения распознанных картин электромагнитной интерференции от возмущенного изображения электромагнитной интерференции.
2. Способ по п. 1, причем этап формирования возмущенного изображения электромагнитной интерференции, основанного по меньшей мере частично на невозмущенном изображении электромагнитной интерференции, включает в себя формирование возмущенного изображения электромагнитной интерференции, основанного по меньшей мере частично на определении фактора объекта, который является функционалом разниц между экспериментальными электромагнитными полями и электромагнитными полями, вычисленными во время этапа формирования невозмущенного изображения электромагнитной интерференции.
3. Способ по п. 2, причем фактор объекта определяют как
Figure 00000031
, где
Figure 00000032
- экспериментально смоделированное или измеренное значение соответственно z-компоненты электромагнитного поля, измеренного приемником j, когда передатчик i является источником электромагнитного поля, где
Figure 00000033
представлено в общем виде как
Figure 00000034
, где α, β и γ являются коэффициентами действительных ненулевых или нулевых значений, где Ω - оператор регуляризации, и где
Figure 00000035
является функцией ее аргумента.
4. Способ по п. 2, причем фактор объекта определяют как
Figure 00000036
, где
Figure 00000032
- экспериментально смоделированное или измеренное значение соответственно z-компоненты электромагнитного поля, измеренного приемником j, когда передатчик i является источником электромагнитного поля, и где
Figure 00000037
- максимальная норма экспериментально измеренной z-компоненты электромагнитного поля.
5. Способ по п. 2, причем фактор объекта определяют как
Figure 00000038
, где
Figure 00000032
- экспериментально смоделированное или измеренное значение соответственно z-компоненты электромагнитного поля, измеренного приемником j, когда передатчик i является источником электромагнитного поля, и где
Figure 00000039
- норма экспериментально измеренной z-компоненты электромагнитного поля, измеренного приемником j, когда передатчик i является источником электромагнитного поля мощности θ.
6. Способ по п. 2, причем этап формирования возмущенного изображения электромагнитной интерференции включает в себя расчет
Figure 00000040
где
Figure 00000041
и
Figure 00000042
представляют собой 3D распределения электромагнитных полей (x, y, z) от электромагнитных источников частоты ƒk, расположенных в местоположениях физических источников (от 1 до N) в системе электромагнитной томографии и в местоположении физических приемников (от 1 до M) соответственно, взятые как сопряженные значения, и причем
Figure 00000043
- «ij»-я составляющая фактора объекта от передатчика i к приемнику j.
7. Способ по п. 2, причем этап распознавания картин электромагнитной интерференции в повторно сформированных возмущенных изображениях электромагнитной интерференции включает в себя вычисление сумм
Figure 00000044
где
Figure 00000041
и
Figure 00000042
представляют собой 3D распределения электромагнитных полей (x, y, z) от электромагнитных источников частоты ƒk, расположенных в местоположениях физических источников (от 1 до N) в системе электромагнитной томографии и в местоположении физических приемников (от 1 до M) соответственно, взятых как сопряженные значения, и причем
Figure 00000043
- «ij»-я составляющая фактора объекта от передатчика i к приемнику j.
8. Способ по п. 2, причем этап распознавания картин электромагнитной интерференции в повторно сформированных возмущенных изображениях электромагнитной интерференции включает в себя вычисление, для итерации i>1,
Figure 00000045
где для простоты частотные члены опущены, где
Figure 00000046
и
Figure 00000047
представляют собой 3D распределения электромагнитных полей (x,y,z) от электромагнитных источников, расположенных в местоположениях физических источников (от 1 до N) в системе электромагнитной томографии и в местоположениях физических приемников (от 1 до M) соответственно, взятых как сопряженные значения, и причем
Figure 00000043
- «ij»-я составляющая фактора объекта от передатчика i к приемнику j.
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, дополнительно содержащая блок отображения, который отображает суперпозиционное изображение.
RU2018117885A 2015-10-16 2016-10-16 Электромагнитная томография с распознаванием картин интерференции RU2720161C2 (ru)

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 (3)

Publication Number Publication Date
RU2018117885A RU2018117885A (ru) 2019-11-21
RU2018117885A3 RU2018117885A3 (ru) 2020-02-17
RU2720161C2 true RU2720161C2 (ru) 2020-04-24

Family

ID=58518250

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2018117885A RU2720161C2 (ru) 2015-10-16 2016-10-16 Электромагнитная томография с распознаванием картин интерференции

Country Status (8)

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

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2783002C1 (ru) * 2021-07-02 2022-11-08 Государственное бюджетное учреждение здравоохранения города Москвы "Научно-практический клинический центр диагностики и телемедицинских технологий Департамента здравоохранения города Москвы" (ГБУЗ "НПКЦ ДиТ ДЗМ") Способ проведения магнитно-резонансной томографии предстательной железы у пациентов с металлоконструкциями тазобедренного сустава

Families Citing this family (6)

* 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
WO2017066731A1 (en) 2015-10-16 2017-04-20 Emtensor Gmbh Electromagnetic interference pattern recognition tomography
WO2018098387A1 (en) 2016-11-23 2018-05-31 Emtensor Gmbh Use of electromagnetic field for tomographic imaging of head
KR20220013537A (ko) * 2018-09-04 2022-02-04 이엠비전 메디칼 디바이시스 리미티드 의료 영상화 장치 및 방법
US12048518B2 (en) 2020-09-22 2024-07-30 Serguei Semenov Individually wearable electromagnetic sensing (iwEMS) system and method for non-invasive assessment of tissue blood and oxygen content
US20240078721A1 (en) * 2021-03-04 2024-03-07 EMvision Medical Devices Ltd Apparatus and Process for Electromagnetic Imaging

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020017905A1 (en) * 1998-04-13 2002-02-14 Ugo Conti Method and apparatus for measuring characteristics of geological formations
US20120083690A1 (en) * 2010-10-04 2012-04-05 Semenov Serguei Y System and method for electromagnetic imaging and therapeutics using specialized nanoparticles
US20140276012A1 (en) * 2013-03-15 2014-09-18 Emtensor Gmbh Wearable/man-portable electromagnetic tomographic imaging

Family Cites Families (89)

* 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
CA2290599A1 (en) 1997-05-23 1998-11-26 The Carolinas Heart Institute Electromagnetical imaging and therapeutic (emit) systems
EP1014853B1 (en) 1997-09-11 2011-12-14 WiSys Technology Foundation, Inc. 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
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
CA2451404C (en) 2001-07-06 2011-04-19 Wisconsin Alumni Research Foundation Space-time microwave imaging for cancer detection
WO2003031958A2 (en) 2001-10-12 2003-04-17 Wisconsin Alumni Research Foundation Microwave dielectric spectroscopy method and apparatus for assays of proteins and other biological molecules
AU2002351405A1 (en) 2001-12-18 2003-06-30 Mri Devices Corporation 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
EP1538972B1 (en) 2002-08-01 2013-11-13 California Institute Of Technology Remote-sensing method and device
WO2004052169A2 (en) 2002-09-27 2004-06-24 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
WO2005115235A1 (en) 2004-05-26 2005-12-08 Medical Device Innovations Limited Tissue detection and ablation apparatus and apparatus and method for actuating a tuner
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
DE602006020618D1 (de) 2005-12-22 2011-04-21 Visen Medical Inc Kombiniertes röntgen- und optisches tomographie-bildgebungssystem
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
CA2729166A1 (en) * 2008-06-27 2009-12-30 Wolfram R. Jarisch High efficiency computed tomography
EP2403401B1 (en) 2009-03-04 2017-05-10 Sensible Medical Innovations Ltd. System for monitoring intrabody tissues
WO2010143691A1 (ja) 2009-06-10 2010-12-16 国立大学法人静岡大学 診断装置
EP2457195B1 (en) 2009-07-23 2017-10-18 Medfields Diagnostics AB Classification of microwave scattering data
US8588284B2 (en) 2010-06-01 2013-11-19 Adeptence, Llc Systems and methods for networked wearable medical sensors
WO2011156810A2 (en) 2010-06-11 2011-12-15 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
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
WO2013005134A2 (en) 2011-07-01 2013-01-10 University Of Manitoba Imaging using probes
DE102011106405B4 (de) * 2011-07-02 2021-08-12 Drägerwerk AG & Co. KGaA Elektroimpedanztomographie-Gerät
IL319973A (en) 2011-12-20 2025-05-01 Sensible Medical Innovations Ltd Chest garment for placement of electromagnetic transducers and methods for using chest garment
US10743815B2 (en) 2012-01-19 2020-08-18 Cerebrotech Medical Systems, Inc. Detection and analysis of spatially varying fluid levels using magnetic signals
CA2863449C (en) 2012-01-19 2019-07-02 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
EP2922464B1 (en) 2012-11-21 2021-09-08 Emtensor GmbH Electromagnetic tomography solutions for scanning head
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
CA2960085A1 (en) 2014-09-03 2016-03-10 Cerebrotech Medical Systems, Inc. Detection and analysis of spatially varying fluid levels using magnetic signals
WO2017066731A1 (en) 2015-10-16 2017-04-20 Emtensor Gmbh Electromagnetic interference pattern recognition tomography
CN110177498B (zh) 2016-12-06 2024-03-15 麦德菲尔德诊断有限公司 用于检测主体中的不对称定位的内部对象的系统和方法
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
AU2018280336B2 (en) 2017-06-08 2024-03-21 EMvision Medical Devices Ltd A tomographic imaging system and process
EP3709868A1 (en) 2017-11-13 2020-09-23 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 (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020017905A1 (en) * 1998-04-13 2002-02-14 Ugo Conti Method and apparatus for measuring characteristics of geological formations
US20120083690A1 (en) * 2010-10-04 2012-04-05 Semenov Serguei Y System and method for electromagnetic imaging and therapeutics using specialized nanoparticles
US20140276012A1 (en) * 2013-03-15 2014-09-18 Emtensor Gmbh Wearable/man-portable electromagnetic tomographic imaging

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
L.MARIAPPAN et.al. Magnetoacoustic tomography with magnetic induction for high-resolution bioimepedance imaging through vector source reconstruction under the static field of MRI magnet., Med Phys. 2014 Feb; 41(2): 022902. *
L.MARIAPPAN et.al. Magnetoacoustic tomography with magnetic induction for high-resolution bioimepedance imaging through vector source reconstruction under the static field of MRI magnet., Med Phys. 2014 Feb; 41(2): 022902. S.Y.SEMENOV et al. Microwave-Tomographic Imaging of the High Dielectric-Contrast Objects Using Different Image-Reconstruction Approaches. IEEE transactions on microwave theory and techniques, vol. 53, no. 7, July 2005. *
S.Y.SEMENOV et al. Microwave-Tomographic Imaging of the High Dielectric-Contrast Objects Using Different Image-Reconstruction Approaches. IEEE transactions on microwave theory and techniques, vol. 53, no. 7, July 2005. *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2783002C1 (ru) * 2021-07-02 2022-11-08 Государственное бюджетное учреждение здравоохранения города Москвы "Научно-практический клинический центр диагностики и телемедицинских технологий Департамента здравоохранения города Москвы" (ГБУЗ "НПКЦ ДиТ ДЗМ") Способ проведения магнитно-резонансной томографии предстательной железы у пациентов с металлоконструкциями тазобедренного сустава

Also Published As

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

Similar Documents

Publication Publication Date Title
RU2720161C2 (ru) Электромагнитная томография с распознаванием картин интерференции
Ito et al. A direct sampling method to an inverse medium scattering problem
Song et al. Study on joint inversion algorithm of acoustic and electromagnetic data in biomedical imaging
Bhateja et al. An improved medical image fusion approach using PCA and complex wavelets
Rahmati et al. Level-set-based reconstruction algorithm for EIT lung images: first clinical results
Xiao et al. Multi-frequency difference method for intracranial hemorrhage detection by magnetic induction tomography
O'Shea et al. Temporal regularization of ultrasound‐based liver motion estimation for image‐guided radiation therapy
Lappas et al. Automatic contouring of normal tissues with deep learning for preclinical radiation studies
Irishina et al. Structural level set inversion for microwave breast screening
EP3472805B1 (en) A method and apparatus for mapping at least part of a structure in an image of at least part of a body of a subject
Xin et al. Time-efficient contrast-enhanced thermoacoustic imaging modality for 3-D breast cancer detection using compressive sensing
Ko et al. U-Net-based approach for automatic lung segmentation in electrical impedance tomography
Joshi et al. Multi-Modality Medical Image Fusion Using SWT & Speckle Noise Reduction with Bidirectional Exact Pattern Matching Algorithm
Qiu et al. Analysis of the magnetoacoustic tomography with magnetic induction
Zhang et al. Difference sensitivity matrix constructed for ultrasound modulated electrical resistance tomography
Abd-Elmoniem et al. Direct pixel to pixel principal strain mapping from tagging MRI using end to end deep convolutional neural network (DeepStrain)
Lindström et al. A hybrid finite element/finite difference method for reconstruction of dielectric properties of conductive objects
Movafagh et al. Huygens principle-based microwave brain imaging through finite difference time domain
Gou et al. Lung dynamic MRI deblurring using low‐rank decomposition and dictionary learning
JP2022059493A (ja) モデル生成方法、モデル生成装置、画像処理方法及び画像処理装置
Movafagh et al. Huygens Principle Imaging Method Powered by Deep Learning for Brain Stroke Classification
Otani et al. Assessing cardiac tissue function via action potential wave imaging using cardiac displacement data
Huy et al. Tomographic image reconstruction enhancement through median filtering and K-means clustering.
Rana et al. Medical image fusion using combined multi-resolution and multi-scaling transform
Hu et al. Score-Based diffusion model for three-dimensional electrical impedance tomography