Конструкция и способ для обнаружения и/или определения местонахождения магнитного материала в области воздействия

Конструкция и способ для обнаружения и/или определения местонахождения магнитного материала в области воздействия

Иллюстрации

Показать все

Реферат

ОБЛАСТЬ ТЕХНИКИ

Настоящее изобретение относится к конструкции для обнаружения и/или определения местонахождения магнитного материала в области воздействия. Настоящее изобретение относится дополнительно к соответствующему способу и к компьютерной программе для управления конструкцией.

УРОВЕНЬ ТЕХНИКИ ИЗОБРЕТЕНИЯ

Конструкция такого рода известна из немецкой патентной заявки DE 10151778 A1. Важно, что в описанной в этой публикации конструкции пространственное распределение напряженности магнитного поля создается таким, что в зоне обследования формируется первая подзона, имеющая относительно низкую напряженность магнитного поля, и вторая подзона, имеющая относительно высокую напряженность магнитного поля. Положение в пространстве подзон в зоне обследования затем сдвигается так, что намагниченность частиц в зоне обследования локально изменяется. Регистрируются сигналы, которые зависят от намагниченности в зоне обследования, причем намагниченность оказывается под влиянием сдвига пространственного положения подзон, и информация относительно пространственного распределения магнитных частиц в зоне обследования извлекается из этих сигналов так, чтобы могло быть сформировано изображение зоны обследования. Такая конструкция имеет преимущество в том, что она может быть использована для исследования произвольных объектов обследования - например, человеческих тел - неразрушающим образом, и не приводя к каким-либо повреждениям, и с высоким пространственным разрешением, и близко к поверхности, и удаленно от поверхности обследуемого объекта.

Подобная конструкция и способ известны из В. Gleich и J. Weizenecker (2005), "Томографическая визуализация с использованием нелинейного отклика магнитных частиц" в Nature, vol.435, 30 June 2005, pp.1214-1217. Описанные в этой публикации конструкция и способ для визуализации магнитных частиц (MPI) используют преимущество кривой нелинейной намагниченности малых магнитных частиц.

MPI представляет собой способ для визуализации распределения магнитных наночастиц, который объединяет высокую чувствительность с возможностью быстрой динамической визуализации, что делает его перспективным для применений в медицинской визуализации. MPI применяет новый способ кодирования сигнала, основанный на динамическом смещении процесса локализованного возбуждения, и позволяет быструю объемную визуализацию. Однако в отличие от установившихся методов визуализации, таких как MRI и СТ, пока еще не подобрано простое математическое преобразование для восстановления изображения по полученным данным. Поэтому восстановление изображения MPI требует сведений о "системной функции", отесывающей отклик системы на данное пространственное распределение частиц, то есть отображения положения частицы в частотном отклике. Для решения проблемы восстановления системная функция должна быть инвертирована, что обычно требует некоторой конструкции регуляризации.

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

СУЩНОСТЬ ИЗОБРЕТЕНИЯ

Цель настоящего изобретения заключается в предоставлении конструкции и способа для обнаружения и/или определения местонахождения магнитного материала в области воздействия, при которых требуется меньше времени для сбора данных о системной функции и/или для восстановления изображения и/или при которых требуется меньший объем памяти для сохранения данных о системной функции.

В первом объекте настоящего изобретения представлена конструкция, как определено в п.1 формулы, содержащая:

- средство выборки для создания выборочного магнитного поля, имеющего такую пространственную конфигурацию напряженности магнитного поля, что в области воздействия формируются первая подзона, имеющая низкую напряженность магнитного поля, и вторая подзона, имеющая более высокую напряженность магнитного поля,

- средство управления для изменения пространственного положения двух подзон в области воздействия посредством управляющего магнитного поля так, чтобы намагниченность магнитного материала изменилась локально,

- средство приема для получения сигналов обнаружения, причем сигналы обнаружения зависят от намагниченности в области воздействия, причем на намагниченность влияет изменение пространственного положения первой и второй подзон,

- средство сохранения для сохранения подмножества данных системной функции конструкции, причем упомянутая системная функция содержит множество данных системной функции, описывающих соотношение между пространственным положением магнитного материала и характеристикой системы для упомянутой конструкции и траектории, вдоль которой упомянутая первая подзона перемещается для сбора упомянутых данных системной функции, и

- средство восстановления для восстановления пространственного распределения магнитного материала в области воздействия по сигналам обнаружения и сохраненному подмножеству данных системной функции, используя дополнительные сведения о структуре системной функции, причем упомянутые средства восстановления приспособлены для создания сокращенного представления системной функции посредством применения преобразования Чебышева или косинусного типа вдоль пространственного измерения и решения задачи намного более разреженного восстановления, возникающей после преобразования.

В дополнительном объекте настоящего изобретения представлен соответствующий способ, как определено в п.11 формулы.

Еще в одном объекте настоящего изобретения представлена компьютерная программа, содержащая средство программного кода для того, чтобы компьютер управлял конструкцией, как заявлено в п.1 формулы, для выполнения этапов способа, как заявлено в п.11 формулы, когда упомянутая компьютерная программа выполняется на компьютере.

Предпочтительные варианты реализации изобретения определены в зависимых формулах. Следует понимать, что заявленный способ имеет подобные и/или идентичные предпочтительные варианты реализации, как и заявленное устройство, и как определено в зависимых формулах.

Настоящее изобретение основано на идее о том, что нет необходимости в получении и сохранении полной системной функции. Из теоретического понимания процесса кодирования сигнала возникло понимание структуры системной функции, что было использовано, чтобы ускорить получение системной функции и/или даже смоделировать ее части или всю системную функцию. Информация о матричной структуре также была использована, чтобы найти более компактные представления системной функции, помогая снизить требования к объему и конфигурации памяти и ускорить восстановление. Наконец, идентификация математического преобразования, ведущего к изображению из данных, может быть использована для ускорения процесса восстановления.

В соответствии с первым объектом настоящего изобретения эта общая идея используется, чтобы сократить область памяти для сохранения данных системной функции. В частности, сохраняется только подмножество данных системной функции, и только это подмножество сохраненных данных системной функции используется при восстановлении пространственного распределения магнитного материала в области воздействия, например, при восстановлении изображения по сигналам обнаружения, получаемым средством приема. Кроме того, сведения о структуре системной функции используются при восстановлении. Эти сведения могут содержать, как определено в предпочтительных вариантах реализации, простые пространственные симметрии, возникающие в частотных компонентах системной функции, представлениях пространственных зависимостей частотных компонент аналитическими функциями или информации о спектральной избыточности в пространственных компонентах. Эти сведения могут быть получены из аналитического анализа конструкции кодирования MPI или моделирования процесса кодирования, и для данной комбинации конструкции, то есть установки сканера, траектории движения первой подзоны, имеющей низкую напряженность магнитного поля (также обозначаемую как "точка, свободная от поля", FFP), и для данного магнитного материала.

Использование сведений о пространственных свойствах системной функции позволяют получить только пространственно сокращенную системную функцию, причем эти сведения используются в соответствии со вторым объектом настоящего изобретения. Тем самым может быть также ускорена калибровочная процедура.

Обычно использование сведений о спектральных свойствах системной функции позволяет сократить размер матрицы системной функции, тем самым ускоряя сбор данных, процесс обработки и восстановления.

В соответствии с предпочтительным вариантом реализации средство приема приспособлено для сбора только упомянутого подмножества данных системной функции конструкции посредством обнаружения сигналов, в то время как зонд упомянутого магнитного материала последовательно помещается во множество различных положений в упомянутой области воздействия. Этот вид сбора данных не очень отличается, в частности, относительно движения FFP через область воздействия, то есть относительно используемой траектории, от сбора сигналов от (реального) обследуемого объекта, например пациента, кроме того факта, что в соответствии с настоящим изобретением зонд не перемещается по всем положениям сетки, покрывающей область воздействия, но только по сокращенному множеству ее точек.

Много различных опций существует для положений, в которых зонд помещается для получения сигналов для сбора упомянутого подмножества данных системной функции. В соответствии с преимущественными вариантами реализации, множество положений располагается в квадранте или октанте упомянутой области воздействия или распределяется чередующимся образом по области воздействия.

Предпочтительно, средства восстановления приспособлены для первого восстановления полной системной функции по сохраненному подмножеству из множества данных системной функции и для последующего восстановления пространственного распределения магнитного материала по сигналам обнаружения, используя восстановленное множество данных полной системной функции. Полная системная функция может, после ее восстановления, также быть сохранена в средстве сохранения для использования в последующих восстановлениях.

Альтернативно, средства восстановления приспособлены для восстановления данных системной функции, которые требуются для восстановления пространственного распределения магнитного материала, но не содержатся в сохраняемом подмножестве множества данных системной функции, непрерывно, в течение восстановления пространственного распределения магнитного материала.

Хотя первый вариант реализации этих двух альтернатив имеет преимущество в том, что фактическое восстановление происходит быстрее, что особенно справедливо для последующих восстановлений, последний вариант реализации имеет преимущество потому, что требуется меньшая занимаемая область памяти.

Как рассмотрено выше, дополнительные сведения используются для восстановления полных данных системной функции (или, по меньшей мере, данных системной функции, требуемых для восстановления пространственного распределения магнитного материала). В предпочтительном варианте реализации используются пространственные симметрии, в частности пространственные зеркальные симметрии, имеющиеся в частотных компонентах системной функции. Такие зеркальные симметрии могут быть получены не только практическими измерениями, но также и из теоретических соображений, используя сведения о траектории (и ее симметрии), кривой намагниченности частиц и монтаже конструкции (и ее симметрии), в частности о ее средстве выборки, средстве управления и приемном средстве. Например, для траектории Лиссажу с частицами Ланжевена наблюдается пространственная зеркальная симметрия, которая имеет вполне определенную частотно-зависимую четность. Для случая 2D визуализации это означает, что только одна четверть полной системной функции должна быть измерена и сохранена.

В соответствии с другим вариантом реализации средство восстановления приспособлено для восстановления пространственного распределения магнитного материала в области воздействия, используя системную функцию с пониженным числом спектральных компонент, полученных из суммирования частотных компонент сохраненного подмножества данных системной функции с подобной пространственной информацией. Поскольку множественные частоты управляющего поля включены в движение свободной от поля точки, сигнал, создаваемый нелинейным откликом частиц, содержит суммы и разности этих частот. Соответственные частотные компоненты имеют одну и ту же пространственную конфигурацию и поэтому могут быть объединены без потери информации. Это приводит к уменьшенному числу строк в матрице системной функции, что означает, что только подмножество данных полной системной функции требуется и изначально сохраняется.

В соответствии с настоящим изобретением средства восстановления приспособлены для создания сокращенного представления системной функции, применяя математическое преобразование по пространственному измерению и решая задачу более разреженного восстановления, возникающую после преобразования. Тип операции преобразования основан на сведениях о функциональной форме пространственной конфигурации спектральных компонент в системной функции. Для траектории Лиссажу FFP спектральные компоненты тесно связаны с функциями Чебышева 2-го рода. Поэтому преобразование Чебышева адекватно для достижения разреженного представления матрицы системной функции и проблемы восстановления, которая сохраняет память и ускоряет восстановление. Иначе говоря, вследствие этого близкого соотношения каждая частотная компонента может быть аппроксимирована только несколькими коэффициентами соответствующего разложения. Альтернативно, может быть использовано преобразование косинусного типа.

Предпочтительно, упомянутое сокращенное представление системной функции создается преобразованием полной системной функции. В соответствии с практическим воплощением этого варианта реализации первыми получаются данные полной системной функции, и из них определяются коэффициенты сокращенного представления, например, функциональное описание с использованием полиномов. Это экономит область памяти и усилия при восстановлении. В соответствии с другим практическим воплощением этого варианта реализации первым получается только подмножество функциональных данных, из которого полные данные системной функции воспроизводятся для определения функционального описания. Альтернативно, функциональное описание, например полиномиальное, может быть приспособлено к полученному подмножеству, как рассматривается ниже.

Еще в одном варианте реализации средство приема приспособлено для непосредственного измерения упомянутого сокращенного представления системной функции в выбранных точках области воздействия, и средства восстановления приспособлены для определения коэффициентов упомянутых аналитических функций, которые используются для описания пространственной зависимости из упомянутого измеренного сокращенного представления системной функции. Точки подмножества выбираются так, что коэффициенты функционального описания могут быть определены по возможности просто.

Следует отметить, что изобретение реализуется не только в заявленных конструкциях, способах и компьютерных программах, но также и в средстве обработки, которое может быть реализовано в аппаратных средствах, программном обеспечении или их комбинации, которые действуют как средство восстановления в конструкции в соответствии с настоящим изобретением. Таким образом, изобретение также относится к средству обработки для восстановления пространственного распределения магнитного материала в области воздействия по сигналам обнаружения и сохраняемому подмножеству данных системной функции, используя дополнительные сведения о структуре системной функции. Дополнительно, настоящее изобретение относится также к соответствующему способу обработки и соответственной компьютерной программе для выполнения упомянутого способа.

КРАТКОЕ ОПИСАНИЕ ЧЕРТЕЖЕЙ

Эти и другие объекты изобретения будут очевидны из рассмотрения описываемого ниже варианта(ов) реализации. На нижеследующих чертежах:

Фиг.1 изображает конструкцию в соответствии с настоящим изобретением,

Фиг.2 - пример конфигурации линий поля, произведенной конструкцией в соответствии с настоящим изобретением,

Фиг.3 - увеличенный вид магнитной частицы, присутствующей в области воздействия,

Фиг.4а и 4b - характеристики намагниченности таких частиц,

Фиг.5 - блок-конструкция аппарата в соответствии с настоящим изобретением,

Фиг.6 - детектируемый временной сигнал s(t) и его спектр S_n,

Фиг.7 - характеристика M(t) намагниченности частиц, получаемый временной сигнал s(t), и значение спектральных компонент S_n для различных управляющих полей и кривых намагниченности частиц,

фиг.8 - диаграммы, иллюстрирующие соотношение между идеальным откликом частиц и смещением выборочного поля,

фиг.9 - диаграммы, иллюстрирующие пространственную зависимость спектральных компонент сигнала для гармонического управляющего поля в комбинации с постоянным градиентом выборочного поля,

фиг.10 - диаграммы, иллюстрирующие производную намагниченности идеальных частиц относительно компоненты поля H_z,

фиг.11 - диаграммы, иллюстрирующие системную функцию идеальных частиц при различных гармониках для 1D движения FFP по центральной линии в z направлении,

фиг.12 - диаграммы, иллюстрирующие последовательные частотные компоненты системной функции идеальных частиц для 2D Лиссажу FFP движения с отношением частот x/z, составляющим 24/25,

фиг.13 - графическая таблица 2D функций Чебышева,

фиг.14 - графики ортогональности для первых 256 компонент базисного набора,

фиг.15 - образец изображения 64×64 и восстановление из разложения на компоненты Чебышева и системной функции, и

фиг.16-19 - использование матрицы системной функции или ее частей для восстановления в соответствии с известным способом и в соответствии с вариантами реализации изобретения.

ПОДРОБНОЕ ОПИСАНИЕ ИЗОБРЕТЕНИЯ

На фиг.1 показан произвольный объект, исследуемый посредством конструкции 10 в соответствии с настоящим изобретением. Номер 350 на фиг.1 обозначает объект, в данном случае человек или животное, размещаемые на столе для пациента, для которого показана только верхняя часть. До применения способа в соответствии с настоящим изобретением магнитные частицы 100 (не показаны на фиг.1) размещались в области воздействия 300 конструкции 10 по изобретению. В частности, перед терапевтической и/или диагностической обработкой, например, опухоли, магнитные частицы 100 помещаются в область воздействия 300, например, посредством жидкости (не показана), содержащей магнитные частицы 100, которая вводится в тело пациента 350.

Как пример варианта реализации настоящего изобретения, на фиг.2 показана конструкция 10, содержащая множество катушек, формирующих средство 210 выборки, диапазон которого определяет область воздействия 300, которая также называется областью обработки 300. Например, средство 210 выборки, размещено выше и ниже пациента 350 или выше и ниже вершины стола. Например, средство 210 выборки содержит первую пару катушек 210′, 210′′, каждая из которых содержит две идентичных обмотки 210′ и 210′′, которые расположены коаксиально выше и ниже пациента 350 и по которым протекают равные токи, но в противоположных направлениях. Первая пара катушек 210′, 210′′, обе вместе, далее называются средством 210 выборки. Предпочтительно, в этом случае используются постоянные токи. Средство 210 выборки создает выборочное магнитное поле 211, которое в целом является градиентным магнитным полем, представленным на фиг.2 полевыми линиями. Оно имеет по существу постоянный градиент в направлении (например, вертикальном) оси пары катушки средства 210 выборки и достигает нулевого значения в точке на этой оси. Начиная с этой, свободной от поля точки (не показанной отдельно на фиг.2) напряженность выборочного магнитного поля 211 увеличивается во всех трех пространственных направлениях при увеличении расстояния от точки без поля. В первой подзоне 301, или области 301, которая обозначена пунктирной линией вокруг точки без поля, напряженность настолько мала, что намагниченность частиц 100, присутствующих в этой первой подзоне 301, не насыщается, тогда как намагниченность частиц 100, присутствующих во второй подзоне 302 (вне области 301), отвечает состоянию насыщения. Точка без поля, или первая подзона 301 области воздействия 300, предпочтительно является пространственно когерентной зоной; она может также быть точечной зоной, или линией, или плоской зоной. Во второй подзоне 302 (то есть в остающейся части области воздействия 300 с внешней стороны первой подзоны 301) напряженность магнитного поля достаточно велика, чтобы сохранить частицы 100 в состоянии насыщения. Изменяя положения двух подзон 301, 302 в пределах области воздействия 300, намагниченность (полная) в области воздействия 300 изменяется. Измеряя намагниченность в области воздействия 300, или физические параметры, связанные с намагниченностью, может быть получена информация о пространственном распределении магнитных частиц в области воздействия. Для изменения относительного пространственного положения двух подзон 301, 302 в области воздействия 300 дополнительное магнитное поле, так называемое магнитное поле 221 управления, накладывается на выборочное поле 211 в области воздействия 300 или, по меньшей мере, в части области воздействия 300.

На фиг.3 показан пример магнитной частицы 100 такого рода, как используется вместе с конструкцией 10 настоящего изобретения. Она содержит, например, сферическую подложку 101, например, из стекла, на которой имеется мягкомагнитный слой 102 толщиной, например, 5 нм, и который состоит, например, из железоникелевого сплава (например, Пермаллоя). Этот слой может быть покрыт, например, покровным слоем 103, который защищает частицу 100 от химически и/или физически агрессивных сред, например кислот. Напряженность выборочного магнитного поля 211, требуемая для насыщения намагниченности таких частиц 100, зависит из различных параметров, например диаметра частиц 100, используемого магнитного материала для магнитного слоя 102, и от других параметров.

В случае, например, диаметра 10 мкм требуется магнитное поле приблизительно 800 А/м (соответствующее приблизительно плотности потока 1 мТ), тогда как в случае диаметра 100 мкм достаточно магнитного поля напряженностью 80 А/м. Еще меньшие значения получаются, когда выбирается покрытие 102 из материала, имеющего более низкую насыщенную намагниченность, или когда толщина слоя 102 уменьшена.

Для дополнительных особенностей предпочтительных магнитных частиц 100 соответствующие части патента DE 10151778, тем самым, включены здесь посредством ссылки, особенно параграфы 16-20 и параграфы 57-61 из патента ЕР 1304542 А2, заявляя приоритет относительно патента DE 10151778.

Размер первой подзоны 301 зависит, с одной стороны, от градиента напряженности выборочного магнитного поля 211 и, с другой стороны, от напряженности магнитного поля, требуемого для насыщения. Для достаточной насыщенности магнитных частиц 100 при напряженности магнитного поля 80 А/м и градиенте (в данном направлении пространства) напряженности выборочного магнитного поля 211, составляющем 160·10³ А/м², первая подзона 301, в которой намагниченность частиц 100 не насыщается, имеет размеры приблизительно 1 мм (в данном направлении пространства).

Когда дополнительное магнитное поле, далее называемое управляющим магнитным полем 221, накладывается на выборочное магнитное поле 210 (или градиентное магнитное поле 210) в области воздействия 300, первая подзона 301 сдвигается относительно второй подзоны 302 в направлении этого управляющего магнитного поля 221; степень этого сдвига увеличивается при увеличении напряженности управляющего магнитного поля 221. Когда наложенное управляющее магнитное поле 221 является переменным полем, положение первой подзоны 301 изменяется соответственно во времени и в пространстве. Выгодно принять или продетектировать сигналы от магнитных частиц 100, расположенных в первой подзоне 301, в другом диапазоне частот (сдвинутом к более высоким частотам), нежели диапазон частот изменения управляющего магнитного поля 221. Это возможно, поскольку частотные компоненты высших гармоник частоты управляющего магнитного поля 221 возникают вследствие изменения намагниченности магнитных частиц 100 в области воздействия 300 в результате нелинейности характеристик намагниченности.

Для создания этих магнитных управляющих полей 221 для любого данного направления в пространстве предоставляются три дополнительных пары катушек, а именно вторая пара 220′ катушек, третья пара 220′′ катушек и четвертая пара 220′′′ катушек, которые вместе далее называются средством 220 управления. Например, вторая пара 220′ катушек создает компоненту управляющего магнитного поля 221, которая ориентирована в направлении оси катушки первой пары 210′, 210′′ катушек, или средства 210 выборки, то есть, например, вертикально. С этой целью по обмоткам второй пары 220′ катушек протекают одинаковые токи в одном и том же направлении. Эффект, который может быть достигнут посредством второй пары катушки 220′, может быть в принципе также достигнут суперпозицией токов в том же самом направлении с противоположными, равными токами в первой паре 210′, 210′′ катушек так, чтобы ток уменьшался в одной катушке и увеличивался в другой катушке. Однако может быть выгодно, и особенно с целью интерпретации сигнала с большим отношением сигнала к шуму, когда постоянное во времени (или квазипостоянное) выборочное поле 211 (также называемое градиентным магнитным полем) и переменное во времени вертикальное управляющее магнитное поле создаются отдельными парами катушек средства 210 выборки и средства 220 управления.

Две дополнительные пары 220′′, 220′′′ катушек предоставляются для создания компонент управляющего магнитного поля 221, которые ориентированы в другом направлении в пространстве, например, горизонтально в продольном направлении области воздействия 300 (или пациента 350) и в перпендикулярном к нему направлении. Если бы третья и четвертая пары 220′′, 220′′′ катушек, типа катушек Гельмгольца (подобно парам катушек для средства 210 выборки и средства 220 управления), использовались с этой целью, то эти пары катушек должны были быть размещены слева и справа от области обработки, или спереди и позади этой области, соответственно. Это могло бы влиять на доступность области воздействия 300, или области обработки 300. Поэтому третья и/или четвертая пары магнитных катушек, или катушки 220′′, 220′′′, также размещаются выше и ниже области воздействия 300 и поэтому конфигурация их обмоток должна отличаться от таковой для второй пары к 220′ катушек. Катушки этого вида хорошо известны в магниторезонансных аппаратах с открытыми магнитами (открытое MRI), в которых пара радиочастотных (RF) катушек располагается выше и ниже области обработки, причем упомянутая пара RF катушек может создавать горизонтальное переменное магнитное поле. Поэтому конструкция таких катушек не нуждается в дополнительных пояснениях в данном случае.

Конструкция 10 в соответствии с настоящим изобретением дополнительно содержит средства 230 приема, которые только конструкциятично показаны на фиг.1. Средство 230 приема обычно содержит катушки, которые способны детектировать сигналы, индуцированные конфигурацией намагниченности магнитных частиц 100 в области воздействия 300. Катушки этого вида хорошо известны в магниторезонансных аппаратах, в которых, например, пара радиочастотных (RF) катушек располагается вокруг области воздействия 300, чтобы иметь настолько большое отношение сигнала к шуму, насколько это возможно. Поэтому конструкция таких катушек не нуждается в дополнительных пояснениях в данном случае.

В альтернативном варианте реализации для показанного на фиг.1 средства 210 выборки постоянные магниты (не показаны) могут быть использованы для создания градиентного выборочного магнитного поля 211. В пространстве между двумя полюсами таких (противоположных) постоянных магнитов (не показаны) формируется магнитное поле, которое подобно таковому на фиг.2, то есть когда противоположные полюса имеют одинаковую полярность. В другом альтернативном варианте реализации конструкции в соответствии с настоящим изобретением средство 210 выборки содержит и, по меньшей мере, один постоянный магнит, и, по меньшей мере, одну катушку 210′, 210′′, как изображено на фиг.2.

Частотные диапазоны, обычно используемые для или в различных компонентах средства 210 выборки, средства 220 управления и средства 230 приема примерно следующие: магнитное поле, создаваемое средством 210 выборки, или вовсе не изменяется во времени, или изменение сравнительно медленное, предпочтительно приблизительно от 1 Гц и приблизительно до 100 Гц. Магнитное поле, создаваемое средством 220 управления, изменяется предпочтительно приблизительно от 25 кГц и приблизительно до 100 кГц. Изменения магнитного поля, к которому должно быть чувствительно средство приема, находится предпочтительно в частотном диапазоне приблизительно от 50 кГц и приблизительно до 10 МГц.

На фиг.4а и 4b показаны характеристики намагниченности, то есть изменение намагниченности М частицы 100 (не показана на фиг.4а и 4b) как функция напряженности Н поля в местоположении этой частицы 100, в дисперсной системе таких частиц. Видно, что намагниченность М не изменяется при напряженности поля более +Н_с и при напряженности поля менее -Н_с, что означает достижение насыщенной намагниченности. Намагниченность М не насыщается между значениями +Н_с и -Н_с.

На фиг.4а показано действие синусоидального магнитного поля H(t) в местоположении частицы 100, где абсолютные значения получающегося синусоидального магнитного поля H(t) (то есть "фиксируемые частицей 100") оказываются ниже напряженности магнитного поля, требуемого для магнитного насыщения частицы 100, то есть в случае, когда дополнительное магнитное поле не активно. Намагниченность частицы 100, или частиц 100, в этих условиях изменяется в пределах значений своего насыщения с частотой изменения магнитного поля H(t). Результирующее изменение намагниченности во времени отмечено обозначением M(t) справа на фиг.4а. Видно, что намагниченность также изменяется периодически и что намагниченность такой частицы периодически полностью обращается.

Пунктирная часть линии в центре кривой обозначает приблизительное среднее изменение намагниченности M(t) как функции напряженности синусоидального магнитного поля H(t). Как отклонение от этой центральной линии намагниченность простирается несколько вправо, когда магнитное поле Н увеличивается от -Н_с до +Н_с, и несколько влево, когда магнитное поле Н уменьшается от +Н_с до -Н_с. Этот известный эффект называется эффектом магнитного гистерезиса, который лежит в основе механизма для производства тепла. Площадь петли гистерезиса, которая сформирована между линиями кривой, форма и размер которой зависят от материала, является мерой для производства тепла при изменении намагниченности.

На фиг.4b показано действие синусоидального магнитного поля H(t), на которое наложено статическое магнитное поле H₁. Поскольку намагниченность находится в насыщенном состоянии, она фактически не подвержена влиянию синусоидального магнитного поля H(t). Намагниченность M(t) остается постоянной во времени при этой площади. Следовательно, магнитное поле H(t) не вызывает изменения состояния намагниченности.

На фиг.5 показана блок-конструкция аппарата 10, показанного на фиг.1. Пара 210′, 210′′ катушек показана конструкциятично на фиг.5 и имеет для ясности обозначение 210. Пара (первое магнитное средство) 210 катушек питается постоянным током от управляемого источника 32 тока, причем упомянутый источник тока управляется блоком 76 управления. Блок 76 управления соединен с компьютером 12, который связан с монитором 13 для отображения распределения магнитных частиц в обследуемой области и с блоком 14 ввода, например с клавиатурой 14.

Пары (второе магнитное средство) 220′, 220′′, 220′′′ катушек соединены с усилителями 41, 51, 61 тока, от которого они получают питающий ток. Усилители 41, 51, 61 тока, в свою очередь, в каждом случае соединяются с источниками 42, 52, 62, переменного тока, который задает временное изменение усиливаемых токов Ix, Iy, Iz. Источники 42, 52, 62 переменного тока управляются блоком 76 управления.

Приемная катушка (средство приема) 230 также показана конструкциятично на фиг.5. Сигналы, индуцированные в приемной катушке 230, подаются на блок 71 фильтрации, посредством которого сигналы фильтруются. Цель этой фильтрации состоит в том, чтобы отделить измеренные значения, которые вызваны намагниченностью в обследуемой области, которая подвержена влиянию изменения положения этих двух частичных областей (301, 302), от других, мешающих сигналов. С этой целью блок 71 фильтрации может быть сконструирован, например, так, чтобы сигналы, которые имеют временные частоты, меньшие, чем временные частоты, на которых работают пары 220′, 220′′, 220′′′ катушек, или меньшие, чем вдвое этих временных частот, не проходили бы блок 71 фильтрации. Сигналы затем передаются через блок 72 усилителя на аналого-цифровой преобразователь 73 (ADC). Оцифрованные сигналы, произведенные аналого-цифровым преобразователем 73, подаются на блок 74 обработки изображения (также называемый средством восстановления), который восстанавливает пространственное распределение магнитных частиц из этих сигналов и соответственного положения, которое предполагалось для первой частичной области 301 первого магнитного поля в обследуемой области в течение приема соответственного сигнала и которое блок 74 обработки изображения получает из блока 76 управления. Восстановленное пространственное распределение магнитных частиц, наконец, передается через блок 76 управления на компьютер 12, который отображает его на мониторе 13.

Аппарат содержит также средство 75 сохранения, например жесткий диск или полупроводниковую память, связанную с блоком 74 обработки изображения, для сохранения получаемых сигналов обнаружения и данных системной функции аппарата 10. Средства 75 сохранения предоставляются для сохранения подмножества данных системной функции конструкции. Средства 74 восстановления (блок 74 обработки изображения) предоставляются для восстановления пространственного распределения магнитного материала 100 в области воздействия 300 по сигналам обнаружения и сохраненному подмножеству данных системной функции, используя дополнительные сведения о структуре системной функции.

Как рассмотрено выше, MPI использует новый способ кодирования сигнала, основанный на динамическом смещении процесса локализованного возбуждения, и позволяет быструю объемную визуализацию. Однако в отличие от установленных методов визуализации, подобных MRI и СТ, до сих пор не было простого математического преобразования для восстановления изображения из получаемых данных. Поэтому восстановление изображения MPI требует сведений о "системной функции", описывающей отклик системы для данного пространственного распределения частиц, то есть отображение положения частицы на частотной характеристике. Чтобы решить задачу восстановления, системная функция должна быть инвертирована, обычно при некоторой конструкции регуляризации.

До настоящего времени, системная функция определ