RU2770544C2 - Коррекция увеличения жесткости рентгеновского излучения при томографической реконструкции с использованием модели ослабления альвареса-маковски - Google Patents
Коррекция увеличения жесткости рентгеновского излучения при томографической реконструкции с использованием модели ослабления альвареса-маковски Download PDFInfo
- Publication number
- RU2770544C2 RU2770544C2 RU2019111292A RU2019111292A RU2770544C2 RU 2770544 C2 RU2770544 C2 RU 2770544C2 RU 2019111292 A RU2019111292 A RU 2019111292A RU 2019111292 A RU2019111292 A RU 2019111292A RU 2770544 C2 RU2770544 C2 RU 2770544C2
- Authority
- RU
- Russia
- Prior art keywords
- model
- simplification
- sample
- atomic number
- density
- Prior art date
Links
- 238000000034 method Methods 0.000 claims abstract description 137
- 230000000694 effects Effects 0.000 claims abstract description 42
- 230000005855 radiation Effects 0.000 claims abstract description 40
- 238000001228 spectrum Methods 0.000 claims abstract description 22
- 230000009467 reduction Effects 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 10
- 230000002441 reversible effect Effects 0.000 abstract description 4
- 239000000126 substance Substances 0.000 abstract description 2
- 231100000289 photo-effect Toxicity 0.000 abstract 1
- 238000012800 visualization Methods 0.000 abstract 1
- 239000000523 sample Substances 0.000 description 63
- 239000000463 material Substances 0.000 description 37
- 238000002591 computed tomography Methods 0.000 description 26
- 238000012545 processing Methods 0.000 description 18
- 238000003384 imaging method Methods 0.000 description 17
- 230000008569 process Effects 0.000 description 15
- 238000004422 calculation algorithm Methods 0.000 description 14
- 238000012937 correction Methods 0.000 description 13
- 230000006870 function Effects 0.000 description 10
- 230000009977 dual effect Effects 0.000 description 9
- 238000004364 calculation method Methods 0.000 description 8
- 230000007246 mechanism Effects 0.000 description 8
- 230000003595 spectral effect Effects 0.000 description 8
- 238000003325 tomography Methods 0.000 description 7
- 239000000203 mixture Substances 0.000 description 6
- 238000010521 absorption reaction Methods 0.000 description 5
- 230000009471 action Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000005510 radiation hardening Methods 0.000 description 5
- 238000002474 experimental method Methods 0.000 description 4
- 238000001914 filtration Methods 0.000 description 4
- 230000033001 locomotion Effects 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 230000011218 segmentation Effects 0.000 description 4
- 238000013459 approach Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 230000003993 interaction Effects 0.000 description 3
- 238000002083 X-ray spectrum Methods 0.000 description 2
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 2
- 229910052782 aluminium Inorganic materials 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000013170 computed tomography imaging Methods 0.000 description 2
- 230000001066 destructive effect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000004927 fusion Effects 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 210000001519 tissue Anatomy 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 239000004809 Teflon Substances 0.000 description 1
- 229920006362 Teflon® Polymers 0.000 description 1
- 238000000862 absorption spectrum Methods 0.000 description 1
- NIXOWILDQLNWCW-UHFFFAOYSA-N acrylic acid group Chemical group C(C=C)(=O)O NIXOWILDQLNWCW-UHFFFAOYSA-N 0.000 description 1
- 229910021417 amorphous silicon Inorganic materials 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 239000012472 biological sample Substances 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 230000021615 conjugation Effects 0.000 description 1
- 238000002247 constant time method Methods 0.000 description 1
- 239000002872 contrast media Substances 0.000 description 1
- 210000003298 dental enamel Anatomy 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000000295 emission spectrum Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 239000003574 free electron Substances 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000004579 marble Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000010603 microCT Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000002040 relaxant effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 231100000027 toxicology Toxicity 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 238000004846 x-ray emission Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/06—Diaphragms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
- G01N23/046—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/03—Investigating materials by wave or particle radiation by transmission
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/10—Different kinds of radiation or particles
- G01N2223/101—Different kinds of radiation or particles electromagnetic radiation
- G01N2223/1016—X-ray
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/30—Accessories, mechanical or electrical features
- G01N2223/303—Accessories, mechanical or electrical features calibrating, standardising
- G01N2223/3035—Accessories, mechanical or electrical features calibrating, standardising phantom
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/40—Imaging
- G01N2223/401—Imaging image processing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/40—Imaging
- G01N2223/419—Imaging computed tomograph
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/40—Imaging
- G01N2223/423—Imaging multispectral imaging-multiple energy imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/60—Specific applications or type of materials
- G01N2223/616—Specific applications or type of materials earth materials
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/408—Dual energy
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/424—Iterative
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Pulmonology (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Immunology (AREA)
- Medical Informatics (AREA)
- Mathematical Analysis (AREA)
- Algebra (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Optics & Photonics (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
Группа изобретений относится к коррекции увеличения жесткости излучения для томографических реконструкций. Способы и устройства, раскрытые в данном документе, предусматривают коррекцию увеличения жесткости излучения в томографической реконструкции с использованием упрощенной модели ослабления Альвареса-Маковски. Типовой способ включает упрощение модели прямой проекции, при этом модель прямой проекции основана на модели ослабления Альвареса-Маковски (AM), причем упрощение модели прямой проекции приводит к упрощению модели ослабления AM за счет одного из допущений, включающих выбор фотоэффекта, постоянной плотности, постоянного атомного номера и плотности, пропорциональной атомному номеру, и выполнение итеративной реконструкции образца с использованием упрощенной модели прямой проекции, при этом итеративную реконструкцию взвешивают по первому спектру, причем полученные в результате измерений данные о визуализации образца, использованные в итеративной реконструкции, получают на первом уровне энергии, и при этом обратная операция итеративной реконструкции не является сопряженной с упрощенной моделью прямой проекции. Группа изобретений обеспечивает коррекцию увеличения жесткости излучения для томографических реконструкций. 2 н. и 18 з.п. ф-лы, 11 ил.
Description
ОБЛАСТЬ ТЕХНИКИ
[0001] Настоящая заявка в целом относится к компьютерной томографии, и конкретнее, к коррекции увеличения жесткости излучения для томографических реконструкций.
УРОВЕНЬ ТЕХНИКИ
[0002] Томография, или компьютерная томография (КТ), представляет собой процесс формирования трехмерной модели объекта путем объединения 2D проекций интересующего объекта, как правило, получаемых с использованием любого типа проникающих частиц или волн. Томография является быстро развивающейся технологией визуализации с широким спектром применения в различных областях, включающих среди прочего медицину, стоматологию, биологию, защиту окружающей среды, токсикологию, минералогию и электронику. В томографических процессах для получения информации различного типа, такой как, например, атомная структура и химический анализ образца, используются различный инструментарий, включая рентгеновские системы, просвечивающие электронные микроскопы (TEM) и/или сканирующие просвечивающие электронные микроскопы (STEM). Набор данных трехмерной (3D) томографии обычно получают путем обратного проецирования серии двумерных (2D) изображений внутренней области образца, полученных под разными углами.
[0003] В течение последних трех десятилетий КТ использовали в качестве неразрушающего метода визуализации для изучения геологических образцов, таких как ископаемые, коллекторные породы и почвы. Число возможных применений возросло с увеличением разрешающей способности и появлением методов визуализации компьютерной микротомографии (микро КТ), которые теперь способны к достижению субмикронного разрешения. Наряду с такими расширенными возможностями использования возникают новые задачи и перспективы, например в отношении коррекции того, что именуют эффектами увеличения жесткости излучения.
[0004] Эффекты увеличения жесткости излучения по меньшей мере частично обусловлены взаимодействием сканируемого материала со сканирующим пучком, которое может не учитываться должным образом в конкретно используемом способе КТ. Более конкретно, эффекты увеличения жесткости излучения могут быть обусловлены отсутствием учета спектрального состава сканирующего пучка и, дополнительно, трудностью обнаружения низкоэнергетических компонентов сканирующего пучка. Трудность в обнаружении низкоэнергетических компонентов сканирующего пучка может быть связана с большим ослаблением низкоэнергетических компонентов сканирующего пучка в сканируемом материале по сравнению с компонентами сканирующего пучка более высоких энергий.
[0005] Как правило, полихроматические данные ослабления, полученные из одного сканирования, содержат недостаточно информации для выполнения томографической реконструкции. Поэтому увеличение жесткости излучения представляет собой проблему, которая в отсутствие исходных данных об образце выглядит неразрешимой с помощью вычислений. Среди многих способов, ранее предложенных для конкретных случаев (таких как получение изображений костей и тканей) либо требующих значительных вычислительных затрат, ни один способ не продемонстрировал масштабного внедрения в технике визуализации для микро КТ. Поэтому, поскольку в последние годы произошло значительное сокращение других артефактов рентгеновской визуализации, увеличение жесткости излучения стало одним из основных источников артефактов визуализации во многих ситуациях и остается одной из ключевых проблем лабораторной рентгеновской томографии. В данном документе предлагаются способы коррекции увеличения жесткости излучения, которые работают в ситуациях, представляющих массовый интерес для визуализации с помощью микро КТ.
КРАТКОЕ ОПИСАНИЕ ГРАФИЧЕСКИХ МАТЕРИАЛОВ
[0006] Фиг. 1А и 1В представляют собой изображения системы КТ для получения сканирований КТ образцов и коррекции артефактов от увеличения жесткости излучения в соответствии с вариантом реализации изобретения по настоящему описанию.
[0007] Фиг. 2 представляет типовую блок-схему способа итеративной реконструкции.
[0008] Фиг. 3 представляет собой типовую блок-схему способа двухэнергетической итеративной реконструкции, использующего модель Альвареса-Маковски (AM) в соответствии с вариантом реализации изобретения по настоящему описанию.
[0009] Фиг. 4 представляет собой типовой способ томографической реконструкции с использованием упрощенной модели AM для получения скорректированных по увеличению жесткости излучения реконструкций в соответствии с вариантом реализации изобретения по настоящему описанию.
[0010] Фиг. 5 представляет собой типовой способ оценки значений параметров в соответствии с вариантом реализации изобретения по настоящему описанию.
[0011] Фиг. 6 представляет собой типовой способ выбора типа упрощения, используемого в способе по фиг.4 в соответствии с вариантом реализации изобретения по настоящему описанию.
[0012] Фиг. 7 представляет собой типовой двухэнергетический способ с использованием вариационных данных в соответствии с вариантом реализации изобретения по настоящему описанию.
[0013] Фиг. 8 представляет собой типовую функциональную блок-схему системы КT в соответствии с вариантом реализации изобретения по настоящему описанию.
[0014] Фиг. 9 и 10 иллюстрируют реконструкции с использованием методик, раскрытых в данном документе, вместе с типовым способом реконструкции без использования коррекции увеличения жесткости излучения.
[0015] Одинаковые номера позиций относятся к соответствующим частям для всех проекций графических элементов.
ПОДРОБНОЕ ОПИСАНИЕ СУЩНОСТИ ИЗОБРЕТЕНИЯ
[0016] Варианты реализации настоящего изобретения описаны ниже в тесной связи с устройством томографической визуализации для микроразмерной или наноразмерной компьютерной томографии малых объектов, в частности образцов цилиндрической формы, с использованием конусных пучков рентгеновского излучения и кольцевой либо спиральной траектории сканирования. Однако понятно, что способы, описанные в данном документе, общеприменимы к широкому кругу различных томографических способов и устройств, включая системы как с конусным, так и параллельным пучком излучения, и не ограничиваются каким-либо конкретным типом устройства, типом пучка излучения, типом объекта, линейным масштабом или траекторией сканирования
[0017] В контексте настоящей заявки и в пунктах формулы изобретения, формы единственного числа включают также и формы множественного числа за исключением случаев, когда в содержании явно подразумевается иное. Кроме того, термин «включает» означает «содержит». К тому же термин «сопряженный» не исключает наличия промежуточных звеньев между сопряженными элементами.
[0018] Никоим образом не следует считать, что данный документ ограничен системами, устройствами и способами, раскрытыми в нем. Наоборот, настоящее изобретение относится ко всем новым и неочевидным признакам и аспектам различных раскрытых вариантов реализации, как выбранным отдельно, так и в различных разноуровневых сочетаниях одного с другим. Раскрытые системы, способы и устройства не ограничены каким-либо конкретным аспектом или признаком, или сочетанием вышеуказанного, а также раскрытые системы, способы и устройства не требуют наличия какого-либо одного или нескольких конкретных преимуществ или решений проблемы. Любые описания принципов действия использованы для облегчения разъяснения, но раскрытые системы, способы и устройства не ограничены такими принципами действия.
[0019] Хотя для удобства изложения действия некоторых из раскрытых способов описаны в определенном последовательном порядке, следует понимать, что этот стиль описания включает в себя изменение порядка их выполнения за исключением случаев, когда нижеследующие конкретные формулировки требуют в описании использования определенной последовательности. Например, в некоторых случаях порядок выполнения действий, изложенных последовательно, может быть изменен либо такие действия могут выполняться одновременно. Помимо этого, для простоты изложения, прилагаемые фигуры могут не отображать различных вариантов использования раскрытых систем, способов и устройств в сочетании с иными системами, способами и устройствами. Кроме того, для описания раскрытых способов в описании иногда используются такие термины, как «изготавливать» и «получать». Эти термины представляют собой абстрактные формы высокого уровня для фактически выполняемых действий. Фактические действия, которые соответствуют этим терминам, изменяются в зависимости от конкретного варианта реализации и легко определяются специалистом в данной области техники.
[0020] В некоторых примерах значения, процессы или устройства обозначены как «самые низкие», «лучшие», «минимальные» или тому подобное. Понятно, что использование таких описаний указывает на возможность выбора из множества используемых функциональных альтернативных вариантов, и для выбранного варианта не требуется быть лучше, меньше или более предпочтительным по сравнению с другими.
[0021] Как было вкратце отмечено, 3D компьютерная томография является основным методом исследования внутренней структуры объектов, и ее привлекательность проистекает из ее удобства, универсальности и неразрушающего воздействия. Тем не менее, в процессе томографической реконструкции из набора 2D рентгенограмм на изображении появляются нежелательные признаки или артефакты. Эти артефакты являются результатом несовершенства инверсионных способов, используемых для построения объемных данных из рентгенограмм. Главным источником этих артефактов является ошибка точного моделирования физического процесса сканирования. Методы реконструкции, которые используют неправильную геометрию, некорректные шумовые модели и не предполагают компенсации относительного перемещения компонентов (или образца), как правило, приводят к негативным результатам. В лабораторной томографии в этой категории чаще всего встречаются артефакты, возникающие в результате неучтенной полихроматичности спектра, производимого источниками рентгеновского излучения. Наоборот, использование простой монохроматической модели приводит к артефактам увеличения жесткости излучения. При прохождении через образец увеличивается жесткость спектра рентгеновского излучения. Рентгеновское излучение более низкой энергии (или мягкое) испытывает более сильное ослабление по сравнению с рентгеновским излучением более высокой энергии (или жестким), и в результате средняя энергия рентгеновского излучения становится выше. Для более жесткого рентгеновского пучка образец создает меньшее ослабление, что приводит к получению несопоставимых результатов измерений. Если спектральные эффекты не учитывать при реконструкции, увеличение жесткости излучения приводит к появлению характерных артефактов в виде следов от банки (эффект чаши) или темных полос, см. Фиг. 9 и 10 в качестве примеров. Эти нежелательные особенности вызывают в реконструкциях изменения величины ослабления по шкале серого цвета. В результате визуально сложно отличить тип материала, из которого изготовлены оба объекта, а автоматическая пороговая сегментация становится невозможной.
[0022] Артефакты увеличения жесткости излучения нарушают линейность процесса проецирования, устанавливая зависимость параметров ослабления отдельного материала от характеристик его пространственного окружения. Это затрудняет простую и надежную идентификацию материалов, что является серьезным препятствием для использования рентгеновских изображений в таких приложениях, как: i) обнаружение неоднородности эмали для стоматологической КТ и ii) определение состава пористых пород в микро КТ для петрофизического моделирования. Понятно, что значительные усилия были направлены на устранение этих артефактов. Методики можно условно разделить на три категории: 1) линеаризация, 2) фильтрация спектра и 3) методы итеративной реконструкции.
[0023] Для монохроматического рентгеновского излучения при энергии E ослабление для индивидуального материала с коэффициентом ослабления μ(E) и толщиной t моделируется в соответствии с законом Ламберта-Бера (уравнение (1) приведено ниже). По закону Ламберта-Бера интенсивность прошедшего излучения I связана с интенсивностью исходного рентгеновского излучения I0 следующим соотношением:
[0024] Поэтому проецируемая величина ослабления -log(I/I0) линейно возрастает с увеличением толщины материала. Метод линеаризации основан на восстановлении этой линейной зависимости, например с использованием клина для калибровки по клину зависимости ослабления от толщины материала. Этот метод пригоден для индивидуального материала, но поскольку для каждого материала требуется уникальная кривая коррекции, он остается приближенным для образцов, состоящих из нескольких материалов. В результате коррекция для образцов, состоящих из нескольких материалов, может приводить к подавлению артефактов, но не способна полностью их устранить.
[0025] Физическая фильтрация спектра (для повышения его монохроматичности) достигается путем размещения некоторого известного материала постоянной толщины между источником рентгеновского излучения и образцом. С помощью этого можно удалить (то есть отфильтровать) большую часть рентгеновского излучения низкой энергии и обеспечить эффективное предварительное увеличение жесткости излучения до проведения визуализации. Эффекты увеличения жесткости излучения могут быть существенно уменьшены, но за счет потери величины потока рентгеновского излучения и отношения сигнал/шум.
[0026] И линеаризация, и фильтрация воздействуют только на внешние проявления проблемы, а не на причину ее возникновения, поскольку они не моделируют реальные физические процессы при построении проекционной модели и чрезмерно оптимистичны в проведении коррекции во всех возможных случаях без учета дополнительной физической информации о каждом конкретном сценарии.
[0027] Ограничения линеаризации и фильтрации мотивируют использование методов итеративной реконструкции, которые в процессе проецирования используют больше информации, с тем, чтобы более достоверно моделировать процесс получения изображений. Был разработан ряд способов с использованием рекуррентных схем, осуществляющих постепенное приближение фактического сканирования к уровню, признанному «достаточно хорошим». В этих способах используют различное количество априорной информации или собственные упрощения и дополнительные предположения.
[0028] Основным недостатком некоторых итеративных схем является их сложность. Способ линеаризации требует проведения оптимизации по большому количеству переменных на каждом итерационном этапе, что для каждой итерации требует выполнения целых восьми вычислительно затратных операций проецирования/обратного проецирования. Поэтому, хотя попытки использования правильной физической модели ослабления в процессе реконструкции необходимы, сложность создания верной с точки зрения физики и статистики физической модели означает, что этим способам еще предстоит найти более широкое применение. По тем же причинам способ приближения с помощью полиномов является неудовлетворительным, поскольку для обеспечения сходимости в нем требуется наличие точной сегментации. Поскольку в реальных материалах процессы сегментации затруднены именно из-за влияния увеличения жесткости излучения, то невозможность проведения идентификации материала во многих случаях является препятствием к его практическому использованию. В другом способе используют предположение, что за счет уменьшения энергии монохроматического рентгеновского излучения величину полихроматического ослабления можно линейно аппроксимировать относительно коэффициента ослабления для каждого полученного пикселя. Это монохроматическое рентгеновское излучение на выходе из материала будет иметь ту же интенсивность, что и полихроматическое рентгеновское излучение, проходящее тот же путь в материале. Это приближение может представлять собой приближенный метод расчета, но в целом, безусловно, неверный. В любом случае этот способ также привносит дополнительную сложность за счет добавления двух дополнительных параметров, не имеющих очевидного физического смысла и требующих внесения поправок в зависимости от образца.
[0029] Модель ослабления рентгеновского излучения Альвареса-Маковски (АМ) является функцией двух параметров материала: плотности материала ρ и атомного номера Z. В сочетании с двухэнергетическим способом визуализации и спектральными данными модель АМ можно использовать с полностью полихроматическими моделями ослабления. Поскольку в модели используют два неизвестных свойства, для получения реальных значений ослабления в качестве функции от Е необходимо иметь два независимых набора информации. В двухэнергетическом способе визуализации эти два набора данных представляют собой сканирования при разных уровнях энергии. Поскольку двухэнергетический способ визуализации учитывает полихроматическую модель ослабления, он по определению предусматривает меры по коррекции увеличения жесткости излучения в процессе реконструкции. Однако двухэнергетический способ требует удвоенного времени сканирования и по меньшей мере вдвое больше вычислений для получения реконструкции обычного сканирования, а также проведения согласования для сравнения и анализа двух сканирований.
[0030] Для решения этих проблем способы, раскрытые в данном документе, включают простые и понятные способы учета физической сущности полихроматического ослабления рентгеновского излучения в итеративных схемах. Для сведения двух неизвестных параметров материала, требуемых в модели AM ослабления рентгеновского излучения, к одному предлагается несколько допущений. При выборе одного из таких упрощений единственная физическая информация, необходимая для проведения коррекции, может содержаться в энергетическом спектре рентгеновского излучения. Используя эти допущения можно провести реконструкцию ослабления в экспериментальном образце при единственной заданной величине энергии рентгеновского излучения E0. Прямой способ моделирует полихроматический эксперимент с использованием одной или максимум двух операций проецирования. Используя допущение о принадлежности всех ошибок E0, можно пренебречь необходимостью использования сопряженной обратной операции, и при этом обратный процесс будет представлять собой только лишь обратное проецирование при одной величине энергии, усредненной по образцу.
[0031] В некоторых вариантах реализации изобретения модель АМ ослабления рентгеновского излучения может быть использована для итеративной реконструкции. Модель AM сдержит компонент, относящийся к ослаблению на основе фотоэлектрического эффекта, и компонент, относящийся к ослаблению на основе комптоновского рассеяния. Применительно к ослаблению рентгеновского излучения можно отметить, что существенное ослабление фотонов в области энергий 5-120 кэВ рентгеновского излучения в основном связано с влиянием фотоэлектрического эффекта и комптоновского рассеяния. Доля каждого компонента во всем образце изменяется в зависимости от материала и величины энергии. Фотоэлектрический эффект возникает при поглощении фотона орбитальным электроном и сопровождается вылетом электрона. Вероятность такого взаимодействия на единицу расстояния (τ) связана с числом атомов на единицу объема (N), атомным номером (Z) и энергией фотона (E), как указано в уравнении (2), приведенном ниже
[0032] Параметр O(Z) относится к функции второго порядка Z, которой можно пренебречь в упрощенной модели AM. Степени n и m являются постоянными и могут принимать значения от 3 до 5. В качестве принятого в модели AM значения в данном документе будет использоваться значение 3. Обратная кубическая зависимость от энергии фотонов означает, что благодаря фотоэлектрическому эффекту с увеличением энергии величина ослабления быстро уменьшается, поэтому даже небольшая разница в расстоянии ослабления приводит к значительному отличию в величинах ослабления при разных уровнях энергии, поэтому фотоэлектрический эффект ответственен за большую часть проявлений увеличения жесткости излучения.
[0033] Комптоновское рассеяние, напротив, остается гораздо более стабильным в рентгеновской области. Комптоновское рассеяние заключается в поглощении свободными электронами фотонов с высокой энергией и их повторной эмиссии, при этом путь фотона отклоняется на угол θ от его исходного направления. Вероятность взаимодействия определяется полным сечением Клейна-Нишины, которое имеет вид уравнения (3).
[0034] Поэтому любая достаточная модель полного ослабления фотонов должна учитывать как фотоэлектрический эффект, так и комптоновское рассеяние. Такая комбинация предложена в модели AM, которая принимает вид уравнения (4)
Здесь величина полного ослабления в материале на единицу длины μ при энергии E представлена в виде линейной комбинации двух явлений. K1 и K2 представляют собой два коэффициента, которые должны быть приведены в соответствие с фактической кривой ослабления, а n - это постоянная, обычно равная 4. Данное упрощение обоснованно пренебрегает учетом влияния краев полосы поглощения рентгеновского излучения, и, следовательно, не может служить точной моделью для элементов с высоким атомным номером, для которых к-края полосы поглощения находятся в обычной рентгеновской области. Однако оно очень хорошо моделирует общую форму спектра поглощения и является простым в том смысле, что приводит зависимость величины ослабления к двум отдельным параметрам материала: плотности ρ и атомного номера Z.
[0035] Фиг. 1А и 1В представляют собой изображения системы КТ 100 для получения сканирований КТ образцов и коррекции артефактов от увеличения жесткости излучения в соответствии с вариантом реализации изобретения по настоящему описанию. Система КТ 100 содержит по меньшей мере источник 102, держатель 104 и детектор 106. Источник рентгеновского излучения 102 выполняет эмиссию пучка рентгеновского излучения конусной формы 108, который передается через объект 110, установленный на держателе образца 104, на детектор 106. Конусный угол и веерный угол определяются соответственно как вертикальный и горизонтальный полууглы, стянутые и угол вентилятора определяются, соответственно, как вертикальные и горизонтальные половинные углы, стянутые детектором 106 в источнике 102. Детектор 106 содержит сцинтиллятор, который при облучении рентгеновскими лучами генерирует световое излучение в видимой области спектра, и светочувствительный датчик CCD, или плоский датчик на основе аморфного кремния, установленный за сцинтиллятором, который генерирует визуальные данные, представляющие собой двумерные отображения пространственного положения световых вспышек, производимых сцинтиллятором, при этом создавая двумерную карту интенсивности рентгеновского излучения, переданного через объект. Понятно, что каждое из этих изображений представляет информацию как о внешних, так и внутренних структурных особенностях объекта в виде проекции вдоль направления распространения рентгеновских лучей, передаваемых через объект на детектор 106. Данные изображения, отправленные детектором 106, получают в виде набора изображений, сохраненных в форме двоичных данных в компьютерном устройстве (не показано) системы для последующей обработки. Изображения получают последовательно, причем держатель образца 104 приводят в движение для поворота образца 110 на небольшой угол (а в случае сканирования по спиральной траектории - также для перемещения образца или объекта 110 на небольшое расстояние) между положениями для получения последовательных изображений, в результате получая разные геометрические проекции объекта по его внутреннему объему. Эти шаги повторяют до завершения полного оборота образца на угол, равный по меньшей мере 180°+веерный угол, и получения полного набора проекционных изображений. В случае спирального сканирования шаги повторяют до завершения образцом 110 линейного вертикального перемещения и поворота, достаточных для получения полной информации об области интереса образца. Хотя в описании для иллюстрации траектории образца используют движение образца 110 относительно источника 102 и детектора 106, движение считают относительным, и в некоторых вариантах реализации изобретения, наоборот, пару источник 102/детектор 106 перемещают относительно образца 110.
[0036] Рентгеновские изображения образца 110 получают последовательно, причем вращение и перемещение выполняют для незначительного вращения и перемещения образца между позициями получения последовательных изображений, в результате получая разные геометрические проекции по внутреннему объему образца 110. Как описано выше, вертикальное перемещение и вращение вокруг вертикальной оси можно сочетать для создания целого ряда возможных траекторий, так что соответствующий набор проекций, полученных по любой из этих траекторий, может предоставить полную информацию об образце 110. Такие траектории могут представлять собой, например, траектории седлового типа, варианты траекторий типа «круг-плюс-линия» и заполняющие пространство траектории. Заполняющие пространство траектории содержат некоторое количество проекционных углов, расположенных в упорядоченной совокупности точек вокруг по меньшей мере части образца 110.
[0037] Как было отмечено, обычные лабораторные системы КТ, такие как система 100, могут выполнять реконструкцию образцов, содержащих артефакты увеличения жесткости излучения, такие как следы от банок (эффект чаши) или темных полос. См. фиг. 9 в качестве примера. Эти артефакты, как правило, обусловлены отсутствием учета спектрального состава излучения от источника 102. Хотя учет спектрального состава источника рентгеновского излучения может приводить к получению более высокого отношение сигнал/шум, и следовательно более качественных изображений, ослабление образцом 110 рентгеновского излучения с более низкой энергией может препятствовать осуществлению требуемого учета с помощью традиционных моделей реконструкции. Например, типовая реконструкция может использовать итеративную схему, которая включает моделирование эксперимента с использованием расчетов по прямой и обратной проекции. Однако если моделирование прямого эксперимента не позволяет учитывать спектральный состав рентгеновского излучения, то качество реконструкции может быть утеряно.
[0038] Как было рассмотрено выше, способы томографической визуализации, раскрытые в данном документе, включают определение упрощенной модели прямой проекции для ее использования в итерационной схеме реконструкции. Определение упрощения может включать выполнение итеративной двухэнергетической реконструкции, при этом двеухэнергетические данные могут представлять собой разреженные низкоэнергетические данные, полученные из низкоэнергетического сканирования с высоким разрешением, и данные с низким разрешением. Данные с низким разрешением могут быть получены из отдельного высокоэнергетического сканирования или основаны на вариационных данных низкоэнергетического сканирования с высоким разрешением. Промежуточные значения, полученные из итеративной двухэнергетической реконструкции, могут быть использованы для определения усредненных параметров образца, таких как плотность ρ и атомный номер Z. Кроме того, при выборе лучшего упрощения можно использовать спектральную информацию источника 102 вместе с оценочными данными по ослаблению. После выбора упрощения выполняют одноэнергетическую итеративную схему вычислений с использованием упрощения для получения томографической реконструкции с коррекцией увеличения жесткости излучения (КЖИ).
[0039] Прямая проекция, использованная в схеме итеративной реконструкции, может представлять собой полную модель AM или упрощенные версии, как определено в данном документе. Как было отмечено, модель AM используют, поскольку эта модель учитывает основные механизмы ослабления рентгеновского излучения и к тому же является энергетически инвариантной по отношению к спектру источника 102. Применительно к ее использованию, система 100 может выполнять низкоэнергетическое сканирование образца 110 с высоким разрешением для выполнения одной или нескольких реконструкций образца. С целью устранения артефактов увеличения жесткости излучения для получения томографических реконструкций образца 110 с коррекцией увеличения жесткости излучения (КЖИ) система 100 может выполнять итеративную схему с использованием упрощения модели AM. За счет упрощения при расчете реконструкции требуются только данные сканирования при одной энергии, например низкоэнергетические данные сканирования с высоким разрешением. Однако, как будет подробнее рассмотрено ниже, при этом также могут быть определены различные другие параметры и проведен базовый выбор наилучшего упрощения для его последующего использования. Следует однако отметить, что различные другие параметры и базовый выбор наилучшего упрощения, могут являться необязательными, а вместо этого дополнительные алгоритмы может заменить первоначальная приближенная оценка.
[0040] Например, можно дать приближенную оценку упрощению, такому как выбор фотоэлектрического эффекта, в качестве единого механизма ослабления и использовать соответствующее упрощение для модели AM в итерационной схеме с одним значением энергией. Дополнительные начальные оценки для ослабления объема также могут быть сделаны на начальном этапе, но затем эти оценки улучшают по мере приближения алгоритма к данным измерения. Важной также является операция обратной проекции, использованная в итерационной схеме. Как было отмечено, в таких итерационных схемах принято использовать сопряженную обратную проекцию, но даже для упрощенной модели AM это может быть обременительно в вычислительном отношении. Поэтому ее можно заменить простым обратным проецированием при одном значении энергии, усредненной по образцу, нарушающим сопряженную связь. Тем не менее, было подтверждено, что такое изменение дает возможность получать требуемые реконструкции с КЖИ при низких вычислительных затратах. Различные алгоритмы и итерационные схемы будут подробнее рассмотрены ниже.
[0041] Фиг. 2 представляет типовую блок-схему способа 201. Способ 201 представляет собой иллюстрацию способа итеративной реконструкции, известного в данной области техники. Способ 201 может начинаться в блоке объединения 205, который суммирует либо величину оценки по объему, либо величину среднего ослабления, полученные из блока 203, с выходной величиной обратной проекции 213. Объединение данных может быть выполнено посредством сложения или умножения. Вначале выходной сигнал обратной проекции 213 может иметь нулевое значение, и значение 203 направляют для получения прямой проекции 207, которая проецирует значение 203 с целью моделирования эксперимента, то есть для имитации ослабления рентгеновского излучения образца. Выходной сигнал 207, т.е. имитация экспериментального значения, подают в блок объединения 209, который сравнивает результаты моделирования с данными измерений 211. Сравнение может быть выполнено путем вычитания или деления и может быть обратным процессу, выполняемому блоком объединения 205. Результат из блока объединения 209 передают в блок обратной проекции 213, который выполняет обратное проецирование. Обратная проекция может при необходимости являться сопряженной обратной проекцией, то есть выполнять операцию, обратную по отношению к выполняемой в блоке прямой проекции 207. Обратная проекция фактически обратно проецирует разностные данные для имитации источника, используемого для получения данных измерений. Выходные данные блока обратного проецирования 213 направляют в блок объединения 205 для объединения с величиной первоначальной оценки, то есть данными блока 203. Затем объединенные данные используют для выполнения следующей итерации по способу 201. Выполнение итераций по способу 201 можно повторять до получения сближения смоделированных данных с измеренными данными.
[0042] Способы итеративной реконструкции ориентированы на получение последовательно улучшенных оценок по объему, которые при проецировании максимально увеличивают вероятность соответствия данных измерения этой оценке. Значения оценки улучшают на каждой итерации путем вычисления разности (не обязательно путем вычитания) между проекцией оценки и измеренными данными и отправки этого значения обратно в поле объемной оценки для получения улучшенного значения для следующей итерации. Процедура осуществления любых методов итеративной реконструкции обобщена в алгоритме 1 (см. ниже). В качестве входных данных в способе используют данные измерений y и на выходе получают реконструкцию объема x. Структура четырех операций (строки 4-7), как правило, аналогична для всех итеративных методов, при этом различия заключаются в выборе и реализации этих четырех функций. Разностная функция образует парную функцию с функцией объединения, которая обычно представляет собой пару (сложение, вычитание) или (умножение, деление). И аналогичным образом функция обратной проекции, которая, как правило, сопряжена с функцией проекции, образует пару с функцией проекции.
1: procedure Iterative(y,n) ♦ Итеративная реконструкция сканирования y с использованием n итераций
2: x ← Initialize() ♦ Выбор начального значения оценки
3: for I=1, 2, 3 … n do
4: p ← P(x) ♦ Текущая оценка проекции
5: Δ ← d(p,y) ♦ Расчет разности между проецированной оценкой и сканированием
6: r ← B(Δ) ♦ Обратное проецирование разности
7: x ← a(r,x) ♦ Обединение с x
8: end for
9: end procedure
[0043] В некоторых вариантах реализации изобретения модель AM может быть использована в прямой проекции 207. Однако, поскольку модели AM содержит по меньшей мере две переменные величины, для полного расчета модели AM может потребоваться многократное сканирование и/или отдельные итерации по способу 201 при разных уровнях энергии. Кроме того, наличие сопряженной обратной проекции модели AM может дополнительно усложнять использование способа 201.
[0044] Фиг. 3 представляет собой типовую блок-схему способа двухэнергетической итеративной реконструкции 301, использующего модель AM в соответствии с вариантом реализации изобретения по настоящему описанию. Способ 301 может предусматривать реконструкцию образца с использованием данных измерения на двух разных уровнях энергии источника, дополнительно взвешенных по спектру, соответствующему каждому из этих энергетических уровней источника. Способ 301 может выполняться по алгоритму 1, рассмотренному выше, с изменениями, внесенными в прямое проецирование, добавлением второго уровня энергии источника и несопряженной обратной проекцией. Несмотря на это, конечным результат может представлять собой реконструкцию образца с коррекцией увеличения жесткости излучения (КЖИ).
[0045] Хотя для реализации полной модели AM можно использовать двухэнергетическую итеративную реконструкцию, для такой модели требуются данные двух низкоэнергетических сканирований с высоким разрешением, что может потребовать дополнительного времени сканирования, дополнительных вычислительных затрат и более сложного обратного проецирования. Однако в двухэнергетической итеративной реконструкции можно использовать данные низкоэнергетического сканирования с высоким разрешением в сочетании с данными низкого разрешения, для получения которых может потребоваться меньше времени. Данные низкого разрешения могут быть получены из высокоэнергетического сканирования, например сканирования с использованием рентгеновских излучения высокой энергии, имеющих по определению изначально более высокую жесткость, или вариационных данных низкого разрешения. В любом случае для пересчета прямых проекций с учетом спектрального состава рентгеновского излучения используют рентгеновские спектры обоих сканирований. Хотя способ двухэнергетической итерационной реконструкции можно использовать для получения реконструкций, при использовании данных низкого разрешения или вариационных данных такие реконструкции могут не являться реконструкциями с полной коррекцией увеличения жесткости излучения. Однако полученные из двухэнергетической реконструкции промежуточные параметры, такие как ρZ3 и ρ, можно использовать в других способах, раскрытых в данном документе, для определения средних значений параметров и упрощения модели AM для получения реконструкции образца с КЖИ.
[0046] Способ 301 может начинаться с использования начальных оценок 303 и 305, которые представляют собой величины объемных оценок или оценок ослабления образца при соответствующих уровнях энергии источника E0 и E1. Начальные значения из блоков 303, 305 отправляют в блоки объединения данных 307, 309 соответственно, в которых происходит либо суммирование, либо умножение данных обратного проецирования 343 и 345 соответственно. Вначале блоки объединения данных 307 и 309 могут не принимать данные обратных проекций 343, 345 для объединения с величинами оценок из блоков 303 и 305, но последующие шаги итерации по способу 301 включают объединение данных с данными из блоков обратного проецирования 343, 345. Выходные данные блоков объединения 307 и 309 подают на блок обработки 311, который выполняет расчет обратной модели AM при уровнях энергии E0 и E1. Например, обратная модель AM может соответствовать обратному уравнению 4. На основе входных данных блок обработки 311 определяет значение плотности (ρ) и атомного номера (Z) для образца. Значения ρ и Z, которые могут называться промежуточными параметрами образца, направляют в блок 313, а значение ρ направляют в блок 317. В блоке 313 определяют значение для ρ * Z3, которое используют в компоненте фотоэлектрического эффекта модели AM, при этом блок 317 используют для сохранения значений ρ, которое используют в компоненте комптоновского рассеяния модели AM.
[0047] Значения ρZ3 и ρ отдельно проецируют с использованием модели AM (уравнение 4) в блоках прямого проецирования 315 и 319 соответственно, при этом отдельные проекции объединяют для моделирования проекции при каждом уровне энергии E с использованием модели AM в блоке обработки 321. Результат выполнения обработки данных в блоке 321 направляют в оба блока объединения данных 325 и 329, где производят их взвешивание по спектрам 0 и 1 соответственно, полученным из 323 и 327. Спектры 0 и 1 относятся к уровням энергии источника E0 и E1 соответственно и могут быть измерены с использованием системы 100 без образца в камере. На выходе блоков объединения 325, 329 получают значения оценки проекции по спектрам 0 и 1 соответственно с использованием модели AM, как было отмечено, в блоках 331 и 333. Эти оценки затем направляют в блоки объединения данных 335 и 337 соответственно, где проводят сравнение значений оценки с данными измерений по тем же спектрам 0 и 1. Блоки объединения данных 335 и 337 могут выполнять обратные операции, такие как выполняют блоки объединения 307 и 309, соответственно, как было отмечено выше. Разницу или остаток после сравнения затем направляют в блоки обратных проекций 343 и 345 соответственно, которые выполняют несопряженное обратное проецирование. Несопряженные обратные проекции могут представлять собой простую обратную проекцию при одном уровне энергии, а не по соответствующим спектрам 0 и 1. Выходные данные блоков обратных проекций 343 и 345 затем направляют в блоки объединения данных 307 и 309 соответственно для выполнения следующей итерации по способу 301. Выполнение итераций по способу 301 можно продолжать до получения сближения данных.
[0048] Типовая итеративная схема, обобщенная в алгоритме 1, принимает вид, представленный, например, на фиг. 2 и 3. Здесь оператор проецирования представлен в виде интегрирования коэффициента затухания объекта μ по линии L, в соответствии со следующим уравнением:
P(μ)(L)=∫L μ(s)ds. (5)
По закону Ламберта-Бера зависимость для интенсивности принимает вид
I(L)=I0e-P(μ)(L). (6)
[0049] Для проведения коррекции увеличения жесткости излучения этот расчет монохроматической интенсивности для полихроматического излучения изменяют следующим образом:
I(L)=I0 ∑ S(E)e-P(μ(E))(L), (7)
где суммирование выполняют по E.
Здесь μ(E) представляет собой ослабление при уровне энергии E согласно модели AM. Получают проекцию величины ослабления при каждом уровне энергии, преобразуют ее в интенсивность и взвешивают по интенсивности спектра для аналогичного уровня энергии S(E) - см., например, блоки 315, 319-325, 329 по способу 301 в качестве примера. Суммирование этих величин по всем уровням энергии E дает в результате окончательное значение интенсивности по линии L. Из уравнения (7) представляется, что оператор проецирования необходимо использовать при каждом уровне энергии E. Однако на практике это может не потребоваться, если проецирование соответствующих параметров материала выполняют отдельно (см. блоки 313-315 и 317-319). Затем из этих данных можно получить значения ослабления при каждом уровне энергии с помощью простых операций сложения и масштабирования. С точки зрения упрощения модели AM для одного из упрощений требуется наличие только одной проекции и масштабирование с коэффициентом 1/E3 (учет только фотоэлектрического эффекта), тогда как для другого (Z пропорционально ρ) требуются две проекции с коэффициентами масштабирования 1/E3 и fKN (E).
[0050] Следует выделить, что использование полной модели AM можно выполнить только при наличии достаточного объема информации, полученной при сканировании при двух уровнях энергии. Именно поэтому и требуются упрощения. Существует также вопрос выбора основного уровня энергии E0, который можно использовать для масштабирования величины ослабления для получения значений ослабления для всех уровней энергии.
[0051] Модель AM предоставляет необходимую информацию об ослаблении рентгеновского излучения в зависимости от энергии рентгеновского излучения и о параметрах материала, таких как атомный номер (Z), атомная масса (A) и электронная плотность (ρ). Модель включает в себя комптоновское рассеяние и фотоэлектрический эффект. Допущение о равенстве A=2Z, сводит модель к двум переменным. При выполнении сканирования образцов, состоящих из одного материала, для которого такие параметры могут быть получены, модель AM можно использовать непосредственно в алгоритме реконструкции. Однако это несколько ограничивает область применения модели, особенно в тех случаях, когда основной задачей томографии является идентификация сканирования, содержащего несколько неизвестных материалов, для которых плотность и эффективный атомный номер не только неизвестны, но могут попиксельно изменяться. Необходимы способы, позволяющие неким образом различать эти параметры материала во время выполнения реконструкции. Это изначально некорректно поставленная задача. Томограмма содержит только один тип информации, которая представляет собой коэффициент ослабления для каждого пикселя, поэтому имеющихся данных недостаточно для использования двух компонентов ослабления. Это связано с тем, что плотность и атомный номер материала являются независимыми величинами. Следовательно, получение необходимых значений для использования полной модели AM требует определенных предположений.
[0052] Один из раскрытых в данном документе способов, то есть упрощение модели AM, включает наличие простой взаимосвязи между двумя величинами: линейной зависимости между Z и ρ, например Z=cρ. Хотя это упрощение дает возможность осуществлять моделирование процесса физического ослабления полихроматического излучения в итерационной схеме, уровень его технологичности остается под вопросом. Поскольку вклады фотоэлектрического поглощения и комптоновского рассеяния в величину ослабления необходимо рассчитывать, а затем выполнять их проецирование и отдельное обратное, сложность алгоритма неизбежно возрастает. Однако, ослабив требования к сопряженной операции обратного проецирования и допуская наличие простых соотношений свойств материала, алгоритм может быть значительно упрощен.
[0053] Фиг. 4 представляет собой типовой способ 401 томографической реконструкции с использованием упрощенной модели AM для получения скорректированных по увеличению жесткости излучения реконструкций в соответствии с вариантом реализации изобретения по настоящему описанию. Способ 401 использует одно из нескольких упрощений модели AM, с тем, чтобы скорректированные по увеличению жесткости излучения реконструкции могли быть получены из сканирования при одном уровне энергии вместо сканирования при двух уровнях энергии. Способ 401 представляет собой упрощение способа 301 и включает упрощение модели AM, которая в свою очередь устанавливает переменные величины, используемые при выполнении способа 401. Например, для упрощения, которое заключается в использовании в модели АМ только фотоэлектрического эффекта, дополнительные блоки 407, 409 и 413 можно исключить. Следует отметить, что дополнительные технологические блоки могут понадобиться только тогда, когда упрощение заключается в пропорциональности Z к ρ, что требует расчета значений c3ρ4 и ρ. В некоторых вариантах реализации изобретения способ 401 можно выполнять с использованием предположений относительно начальных значений блока 403, и при необходимости, блоков 409 и 413, при этом случайным образом осуществляя выбор упрощения модели AM, например, упрощения пропорциональности или только упрощение выбора фотоэлектрического эффекта. Однако в других вариантах реализации изобретения, на этапах, предшествующих выполнению способа 401, можно выполнять способы 301, 501 и 601 с целью определения начальных значений c, ρ и Z и дополнительного определения наиболее пригодного упрощения, используемого в способе 401. Следует отметить, что в данном документе предусмотрены обе последовательности операций, при этом изначальное определение c, ρ, Z и выбор упрощения не являются необходимыми для способа, раскрытого в данном документе, который заключается в упрощении модели AM для прямого проецирования и использования несопряженной обратной проекции в итерационной схеме.
[0054] Как и в случае способов 201 и 301, способ 401 может начинаться в блоке объединения данных 405, который объединяет начальные значения из 403 с выходными данными обратного проецирования 429. Тем не менее, вначале выходные данные обратного проецирования 429 могут быть равными нулю. Действие, выполняемое в блоке объединения данных 405, может представлять собой сложение или умножение. В зависимости от упрощения, как было отмечено выше, выходные данные из блока объединения 405 могут направляться непосредственно в блок обработки 417 или в блок обработки 407. Например, если упрощение основано на постоянной плотности, постоянном атомном номере, только фотоэлектрическом эффекте или только комптоновском рассеянии, то выходные данные блока обработки 405 направляют в блок обработки 417. С другой стороны, в случае выбора упрощения пропорциональности атомного номера величине выходные данные из блока объединения 405 подают в блок обработки 407.
[0055] Поскольку способ 401 является упрощенным вариантом способа 301, детали последовательности операций не рассматриваются для краткости. Тем не менее, влияние упрощений и способы их реализации будет рассмотрено. Следует также повторить, что обратное проецирование, выполняемое в блоке обратного проецирования 429, является несопряженным проецированием, например, эта операция не является обратной для операции проецирования, выполняемой в блоках 411 и/или 417. В целом предлагается семейство алгоритмов для определения томграфической карты ослабления при одном уровне энергии рентгеновского излучения. Алгоритмы образованы с использованием различных упрощений модели AM (для сведения двух параметров материала к одной переменной) в сочетании с несколькими способами итеративной реконструкции по максимальной вероятности. Во всех упрощениях требуется создание простой модели, позволяющей точнее моделировать результаты измерения измерение физических данных, при этом операция обратного проецирования не отличается от типовой операции обратного проецирования в допущении «монохроматического» излучения.
[0056] Одно из упрощений заключается в допущении постоянной величины плотности ρ, что может быть полезным для разделения в образце различных материалов с одинаковой плотностью, таких как биологические образцы с контрастными агентами, поглощенными различными тканями образца. При использовании этого упрощения важная задача относится к успешности выбора постоянного значения плотности.
[0057] Другое упрощение заключается в допущении постоянного атомного номера Z, что может быть применимо, когда образец состоит из одного материала, но с различной степенью пористости (например, плотности). В этом допущении масштабирование по уровню энергии учитывает как фотоэлектрический эффект, так и комптоновское рассеяние и благодаря упрощению может требовать выполнения только одной операции проецирования на шаг итерации. Выбор значения Z может быть сделан на основе априорных знаний.
[0058] Другое упрощение состоит в допущении выбора фотоэлектрического эффекта в качестве единственного механизма ослабления. Хотя этот подход может показаться огромным упрощением, такое упрощение приводит к действенной качественной коррекции артефактов увеличения жесткости излучения, поскольку большая часть несоответствий возникает при низких энергиях, когда фотоэлектрическое поглощение дает максимальный вклад в величину ослабления. Упрощение включает масштабирование фотоэлектрического эффекта 1/E3, что указывает на более выраженную разницу в ослаблении между двумя различными уровнями энергии для фотоэлектрического компонента модели AM по сравнению с масштабированием fKN(E) комптоновского компонента модели. Масштабирование может подчеркивать чувствительность фотоэлектрического компонента модели AM к спектру зондирующего рентгеновского пучка и эффекту увеличения жесткости излучения. Дополнительное преимущество может заключаться в необходимости выполнения только одной операции проецирования на шаг итерации по способу 401.
[0059] Четвертое упрощение заключается в допущении выбора комптоновского рассеяния в качестве единственного механизма ослабления, которое может являться обратным упрощением по отношению к фотоэлектрическому эффекту и может пренебречь только вкладом фотоэлектрического эффекта. Это упрощение может иметь меньше применений по сравнению с другими, перечисленными в данном документе, но включено в него с тем, чтобы получить полное представление о возможных упрощениях. Это упрощение может быть использовано в томографии, использующей рентгеновские лучи очень жесткого спектра, например, только с высокой энергией, и/или образцы, состоящие полностью из очень легких элементов. Использование этого упрощения может также потребовать выполнения только одной операции проецирования на шаг итерации по способу 401.
[0060] Последнее упрощение заключается в допущении пропорциональности атомного номера Z величине плотности ρ. Например, Z может быть пропорциональным величине cρ, при этом c является постоянной пропорциональности. Примерное значение c может составлять около 4,88, что дает возможность построения приемлемой модели для материалов из середины периодической таблицы. Вместо установления соотношения между компонентами фотоэлектрического эффекта и комптоновского рассеяния в модели AM, предлагается соотношение между коэффициентами фотоэлектрического эффекта и комптоновского рассеяния ρZ3=c3ρ4 и ρ соответственно. Отметим, что для этого допущения требуются две операции проецирования в прямой модели: одна для ρ4 и одна для ρ, поскольку ∫Lρ(x)4ds ≠ (∫Lρ(x)ds)4.
[0061] Для дополнительной иллюстрации вышеупомянутые упрощения будут подробно рассмотрены по отношению к модели AM и их использования по меньшей мере в способе 401. Эти упрощения используют для изменения прямой модели в каждой из итерационных схем. Конечно, изменение только прямой модели приводит к нарушению сопряженности с обратной операцией (в данном случае - обратным проецированием), которая, как правило, предусматривает быструю сходимость. Наблюдалось и, что, несмотря на эту теоретическую несовместимость, все эти способы показали устойчиво работу с имитированными полихроматическими сканированиями и приемлемую скорость сходимости.
[0062] Атомный номер, пропорциональный плотности
[0063] Упрощение полной модели AM на основе допущения Z=cρ приводит к следующему:
И масштабирование принимает вид
что может потребовать дополнительных шагов по сравнению с другими упрощениями, такими как упрощение выбора только фотоэлектрического эффекта. Во-первых, для использования уравнения (8) необходимо рассчитать ρ для каждого пиксела. Затем выполняют проекцию величин ρ и ρ4. Масштабирование по уровням энергии выполняют после проецирования, поскольку масштабирование является линейным по отношению к двум компонентам: ρ и ρ4, но не линейно по отношению к μ. В итоге последовательность выполнения принимает вид
(а) Решение упрощенного уравнения АМ для получения значений р в каждом пикселе, что дает величины ρ и ρ4;
(b) Проецирование ρ;
(c) Проецирование ρ4;
(d) Получение из уравнения (9) проецированной величины ослабления для каждого уровня энергии и проекционных данных P(ρ) и P(ρ4); и
(e) Преобразование проецированной величины ослабления в интенсивность и проведение суммирования после осуществления взвешивания по спектру по уравнению (7).
[0064] Фотоэлектрический эффект как единственный механизм ослабления
[0065] Упрощение выбора фотоэлектрического эффекта в качестве единственного механизма ослабления пренебрегает влиянием компонента модели AM, относящегося к комптоновскому рассеянию, и приводит к следующему
Следует отметить, что масштабирование выполняют после операции проецирования. Это возможно, поскольку значения ослабления для всех материалов находится на одной кривой 1/E3, поэтому для получения проецируемой величины ослабления при другом уровне энергии для каждого пиксела проекции масштабирование проводят аналогично для всех различных материалов. Следовательно, этап проецирования можно обобщить следующим образом:
(a) Проецирование текущей оценки значений коэффициента ослабления; принять, что эта оценка соответствует уровню энергии E0;
(b) Масштабирование проецированного при E0 ослабления по уравнению (11) для получения проецированной величины ослабления для каждого уровня энергии; и
(c) Преобразование проецированной величины ослабления в интенсивность и проведение суммирования после осуществления взвешивания по спектру согласно уравнению (7).
[0066] Комптоновское рассеяние как единственный механизм ослабления
[0067] Это упрощение аналогично вышеописанному упрощению выбора фотоэлектрического эффекта в качестве единственного механизма ослабления. Тем не менее, это упрощение модели AM включает только компонент комптоновского рассеяния. Это упрощение приносит заметно меньше пользы, поскольку фотоэлектрический эффект вызывает наибольшие изменения ослабления и, следовательно, является причиной большинства артефактов увеличения жесткости излучения. Использование этого упрощения приводит к следующей упрощенной модели AM
При масштабировании соотношение принимает вид
Операция проецирования будет иметь такой же вид, как в приведенном выше случае.
[0068] Атомный номер постоянен во всем образце
[0069] Если допустить, что величина эффективного атомного номера Z практически неизменна во всем образце, то плотность ρ становится единственной переменной, и в результате заранее определенное Z является достаточным для всего изображения. Модель AM, основанная на этом упрощении, сводится к следующему виду
Алгоритм масштабирования после операции проецирования путем последующего масштабирования по постоянному значению принимает вид
[0070] Плотность постоянна во всем образце
[0071] Если допустить, что величина плотности практически постоянно во всем образце, возможно, потребуется учитывать плотность объема пространства вне образца. Хотя это может усложнить способ, учесть эту разницу довольно просто. Здесь имеется два значения плотности, которые можно принять: значение ρ внутри образца и 0 г/см3 за пределами образца. Однако отметим, что P(ρZ3)=ρP(Z3) и P(ρ)=ρP(B), где B - бинарная сегментация объема. Z3 может быть получено из объемного значения μ(E0) следующим образом:
B можно рассчитать из этого объема, присваивая ему значение 1 при Z3 > 0 или 0 в противном случае. Проецируемая величина ослабления образца при любом уровне энергии E может быть рассчитана как
Следует отметить, что сегментированный объем B вряд ли сильно изменится с последующими итерациями, и P(B) можно рассчитать на первом шаге итерации, а затем сохранять это значение для повторного использования на всех последующих шагах итерации.
[0072] Возможно, потребуется повторить, что в этом случае операция прямого проецирования не является линейной, и поэтому здесь теряется смысл традиционного понятия сопряженной операции, которое применимо только в контексте умножения линейной матрицы. Тем не менее, при этом сохраняется возможность осуществления псевдосопряженной операции, которая соответствует понятию сопряженных операций и поддерживает совместимость обработки монохроматических данных в поле объемов и полихроматических данных в поле интенсивности/ослабления. Эти методы получены эмпирически, особенно для случаев подбора приемлемых нормировочных множителей, которые сохраняют объемную форму (или сумму значений ослабления для изображения в объемном пространстве) во всех итерациях. Поскольку каждая обратная операция образует пару со своей прямой моделью, в целом требуются пять соответствующих псевдосопряженных операций, т.е. по одной операции для каждого допущения, рассмотренного выше. Для каждой прямой модели доступны две обратные операции: неизмененное обратное проецирование, как в монохроматическом случае; и эмпирически полученное псевдосопряженное обратное проецирование. Хотя по меньшей мере эти две операции доступны, в данном документе рассматривается немодифицированное обратное проецирование, не имеющее сопряженной связи с прямым проецированием.
[0073] Фиг. 5 представляет собой типовой способ 501 оценки значений параметров в соответствии с вариантом реализации изобретения по настоящему описанию. Рассчитанные параметры можно впоследствии использовать в других способах, раскрытых в данном документе, и с их помощью дополнительно выбирать тип упрощения для использования в данной реконструкции. Выбор типа упрощения будет рассмотрен по отношению к фиг. 6. После проведения оценки параметров и выбора упрощения можно выполнять способ 401 с использованием выбранного упрощения. Все способы могут быть выполнены, например системой 100.
[0074] Входные данные ρZ3 и ρ можно получать при выполнении двухэнергетического способа 301 с использованием разреженных низкоэнергетических данных в сочетании с данными сканирования с низким разрешением, данными высокоэнергетического сканирования или вариационными данными (см. Фиг. 7 для способа, основанного на вариационных данных). Промежуточные параметры из блоков 313 и 317 являются начальными входными данными относительно способа 501. Промежуточные значения могут использоваться в способе 501 после получения схождения данных по способу 301 с использованием разреженных низкоэнергетических данных в сочетании с данными сканирования с низким разрешением, данными высокоэнергетического сканирования или вариационными данными.
[0075] Определение значения параметра включает определение средней плотности ρ, среднего атомного номера Z и постоянную пропорциональности c, используемую в соотношении Z=cρ. Значение средней плотности, получаемое из блока обработки 509, может представлять собой среднее значением плотности, найденное для образца. Среднее значение атомного номера Z может быть найдено как среднее значение
в образце. А постоянную пропорциональности c находят с помощью следующего уравнения, также усредненного по всему образцу
[0076] Фиг. 6 представляет собой типовой способ 601 выбора типа упрощения, используемого в способе 301 в соответствии с вариантом реализации изобретения по настоящему описанию. Способ 601 включает выполнение проекций с использованием всех упрощений, основанных на первоначальной оценке и спектре излучения с уровнем энергии E0 и последующее сравнение проецируемых данных с данными измерений при E0. Например, начальную оценку поучают из блока 603, а информация о спектре получают из блока 605, затем оба значения передают в блок обработки 607. В блоке обработки 607 получают проекцию с использованием данных блоков 603 и 605 для каждого из рассмотренных в данном документе упрощений, при этом выходные данные отправляют в соответствующие блоки 609-617. Затем выполняют сравнение выходных данных блоков 609-617 с результатами, полученными в двухэнергетическом способе 301 с использованием данных низкой энергии. Затем проводят выбор допущения, наиболее близкого к результату двухэнергетического способа, для его использования в способе 401.
[0077] Фиг. 7 представляет собой типовой двухэнергетический способ 701 с использованием вариационных данных в соответствии с вариантом реализации изобретения по настоящему описанию. Способ 701 может являться способом аналогичным способу 301 за исключением использования величин E0, E1 и спектра 0. Способ 701 заменяет высокоэнергетические данные сканирования с низким разрешением для выполнения двухэнергетической итерационной схемы, при этом в нем используют вариационные низкоэнергетические данные по спектру 0. Следует отметить, что вариационные данные аналогичны данным ослабления, полученным с использованием спектра, масштабированного по энергии рентгеновского излучения E0, например сдвинутого в область более высоких энергий. В результате этого вариационные данные можно использовать вместо данных высокоэнергетического сканирования.
[0078] Фиг. 8 представляет собой типовую функциональную блок-схему системы КT 800 согласно варианту реализации изобретения по настоящему описанию. Система КТ 800 может представлять собой пример системы КТ 100 и может использоваться для выполнения сканирований КТ образцов и получения реконструкций образца, такого как образец 110, с использованием КЖИ. Система КТ 800 может по меньшей мере содержать контроллер 820, память 822, держатель 824, программы 826, детектор 828 и источник 830. В общем случае система КТ 800 может выполнять сканирование образцов, как описано в данном документе.
[0079] Контроллер 820 может представлять собой микроконтроллер, ЦП, одно или несколько вычислительных ядер или тому подобное. Контроллер 820 может управлять различными характеристиками системы КТ 800 и может быть сопряжен для приема данных от и отправки сигналов управления к иным различным компонентам системы КТ 800. Например, контроллер 820 может выполнять машинные коды программ 826 с подачей команд для системы КТ 800 на выполнение сканирований КT низкого и высокого разрешения. Кроме того, данные сканирования, полученные детектором 828, могут быть направлены в контроллер 820. В ответ контроллер 820 может выполнять набор алгоритмов обработки данных, таких как алгоритмы, проиллюстрированные способами с 201 по 701, для получения реконструкций с проведенной коррекцией увеличения жесткости излучения. Как было отмечено, алгоритмы могут составлять часть программ 826.
[0080] Память 822 может быть энергозависимой или энергонезависимой памятью, установленной для хранения данных и команд. Программы 826, которые могут составлять часть памяти 822, могут содержать управляющие команды и программы для выполнения различных задач системы КТ 800, таких как управление держателем, источником и детектором, а также различными другими вспомогательными компонентами.
[0081] Держатель 824 может поддерживать образец и осуществлять его перемещение в ответ на управляющие сигналы, отправленные контроллером 820. Контроллер 820 может генерировать управляющие сигналы при выполнении команд управления держателем, хранящихся в памяти 822 и/или программах 826. Например, программа управления держателем может выполнять перемещение держателя в соответствии с заданной траекторией, например, такой как спираль или зигзагообразная линия.
[0082] Детектор 828 может обнаруживать рентгеновские лучи, прошедшие через образец, и регистрировать их интенсивность. Интенсивность рентгеновских лучей может быть уменьшена вследствие ослабления при прохождении через образец. Величина ослабления может представлять собой данные, используемые для визуализации образца и для получения реконструкций. Данные изображения или данные сканирования могут быть направлены в контроллер 820 для обработки или сохранены в памяти 822 для последующей обработки. Кроме того, детектор 828 может перемещаться в ответ на управляющие сигналы контроллера 820. Детектор 828 можно настраивать для учета формы, размера образца или для получения требуемого поля зрения (ПЗ).
[0083] Источник 830 направляет рентгеновские лучи на образец и может управляться с помощью контроллера 820. Например, источник 830, в ответ на управляющие сигналы от контроллера 820, может испускать рентгеновские лучи с необходимой интенсивностью и может также перемещаться в различных направлениях относительно образца. Источник 830 можно настраивать для учета формы, размера образца или для получения требуемого поля зрения (ПЗ).
[0084] Многие варианты очевидны для специалистов в данной области техники, не выходя за пределы объема настоящего изобретения.
[0085] Фиг. 9 и 10 иллюстрируют реконструкции с использованием методик, раскрытых в данном документе, вместе с типовым способом реконструкции без использования коррекции увеличения жесткости излучения. Фиг. 9 содержит томографические реконструкции алюминиевого и мраморного стержня, причем алюминиевый стержень находится слева. Изображение A представляет собой реконструкцию без проведения коррекции, на которой видны артефакты, такие как следы от банок (более яркие края) или темные полосы (затененная область между двумя стержнями). Изображение B представляет собой реконструкцию с КЖИ с использованием упрощения выбора фотоэлектрического эффекта, а изображение C представляет собой реконструкцию с КЖИ с использованием допущения плотности, пропорциональной атомному номеру. Как можно увидеть, оба упрощения улучшают реконструкции. Тем не менее, упрощение выбора фотоэлектрического эффекта по меньшей мере визуально создает меньше штриховых полос одновременно с незначительным уменьшением следов от банок. С другой стороны, упрощение пропорциональности, по-видимому, содержит некоторую полосу, но демонстрирует более слабые следы от банок. При этом обе реконструкции с КЖИ демонстрируют улучшенный контраст между двумя стержнями, что приводит к улучшению идентификации материала.
[0086] Фиг. 10 представляет аналогичную последовательность изображений, но при этом изображения содержат тефлоновый стержень (внизу), стеклянный стержень (вверху) и акриловый стержень (вверху справа). Как и в предыдущем случае, оба упрощения улучшают качество изображения и уменьшают артефакты с различной степенью улучшения.
[0087] Варианты реализации изобретения, рассмотренные в данном документе для иллюстрации раскрытых способов, не следует считать ограничивающими, поскольку они являются лишь примерами их использования. Для специалистов в данной области техники понятны и другие безграничные возможности использования раскрытых способов, предусмотренных в данном документе и находящихся в границах объема настоящего изобретения.
Claims (30)
1. Способ реконструкции томографического изображения образца, включающий:
упрощение модели прямой проекции, при этом модель прямой проекции основана на модели ослабления Альвареса-Маковски (AM), причем упрощение модели прямой проекции упрощает модель ослабления AM только для одного из фотоэлектрического эффекта, постоянной плотности, постоянного атомного номера и плотности, пропорциональной атомному номеру; и
выполнение итеративной реконструкции образца с использованием упрощенной модели прямой проекции и обратной проекции, причем измеренные данные изображения образца, полученные при первой энергии, сравнивают с выходным сигналом упрощенной прямой проекции, взвешенным по первому спектру, связанному с первой энергией, и при этом обратная проекция итеративной реконструкции не сопряжена с упрощенной моделью прямой проекции, и при этом обратная проекция предполагает монохроматическое излучение.
2. Способ по п. 1, причем упрощение модели прямой проекции только для фотоэлектрического эффекта включает игнорирование компоненты эффектов комптоновского рассеяния в модели ослабления AM, так что реализуется только компонента фотоэлектрического эффекта.
3. Способ по п. 1, причем упрощение модели прямой проекции для постоянной плотности включает проецирование модели ослабления АМ с плотностью, установленной на постоянной величине.
4. Способ по п. 1, причем упрощение модели прямой проекции для постоянного атомного номера включает проецирование с модели ослабления АМ с атомным номером, установленным на постоянной величине.
5. Способ по п. 1, причем упрощение модели прямой проекции для плотности, пропорциональной атомному номеру, включает установление атомного номера на параметр плотности образца, взвешенный по постоянной пропорциональности и заменяющий атомный номер в модели ослабления AM.
6. Способ по п. 1, причем упрощение модели прямой проекции включает сокращение числа переменных с двух до одной, причем две переменные включают плотность и атомный номер, и при этом сокращение числа переменных включает превращение плотности в функцию атомного номера или превращение атомного номера в функцию плотности.
7. Способ по п. 1, дополнительно включающий:
определение упрощения модели ослабления АМ на основании определения упрощения, которое минимизирует остаточную ошибку при выполнении итеративной реконструкции.
8. Способ по п. 1, причем выполнение итеративной реконструкции образца с использованием упрощенной модели прямой проекции уменьшает артефакты увеличения жесткости излучения в реконструкции.
9. Способ по п. 1, причем высокоэнергетические данные образца, полученные при сканировании с низким разрешением и высокой энергией, используют для определения упрощения и начальных значений плотности и атомного номера, используемых в итеративной реконструкции.
10. Способ по п. 1, причем в итеративной реконструкции используют полученные с высоким разрешением и низкой энергией данные образца.
11. Система для реконструкции томографического изображения образца, содержащая:
держатель для установки образца;
источник рентгеновского излучения, соединенный с возможностью подавать рентгеновское излучение для получения изображений образца;
детектор, соединенный с возможностью принимать рентгеновское излучение после его прохождения через образец; и
контроллер, соединенный по меньшей мере с источником рентгеновского излучения и детектором, причем контроллер содержит или соединен с одним или более процессорами, выполненными с возможностью исполнять сохраненный в памяти код, который, будучи исполняемым, заставляет контроллер:
упрощать модель прямой проекции, при этом модель прямой проекции основана на модели ослабления Альвареса-Маковски (AM), причем упрощение модели прямой проекции упрощает модель ослабления AM только для одного из фотоэлектрического эффекта, постоянной плотности, постоянного атомного номера и плотности, пропорциональной атомному номеру; и
выполнять итеративную реконструкцию образца с использованием упрощенной модели прямой проекции и обратной проекции, причем измеренные данные изображения образца, полученные при первой энергии, сравниваются с выходным сигналом упрощенной прямой проекции, взвешенным по первому спектру, связанному с первой энергией, и при этом обратная проекция итеративной реконструкции не сопряжена с упрощенной моделью прямой проекции, и причем обратная проекция предполагает монохроматическое излучение.
12. Система по п. 11, причем код для исполнения упрощения модели прямой проекции только для фотоэлектрического эффекта дополнительно содержит код, который, будучи исполняемым, заставляет контроллер игнорировать компоненту эффектов комптоновского рассеяния в модели ослабления AM, так что реализуется только компонента фотоэлектрического эффекта.
13. Система по п. 11, причем код для исполнения упрощения модели прямой проекции для постоянной плотности дополнительно содержит код, который, будучи исполняемым, заставляет контроллер проецировать модель ослабления АМ с плотностью, установленной на постоянной величине.
14. Система по п. 11, причем код для исполнения упрощения модели прямой проекции для постоянного атомного номера дополнительно содержит код, который, будучи исполняемым, заставляет контроллер проецировать модель ослабления АМ с атомным номером, установленным на постоянной величине.
15. Система по п. 11, причем код для исполнения упрощения модели прямой проекции для плотности, пропорциональной атомному номеру, дополнительно содержит код, который, будучи исполняемым, заставляет контроллер устанавливать атомный номер на параметр плотности образца, взвешенный по постоянной пропорциональности и заменяющий атомный номер в модели ослабления AM.
16. Система по п. 11, причем код для упрощения модели прямой проекции дополнительно содержит код, который, будучи исполняемым, заставляет контроллер сокращать число переменных с двух до одной, причем две переменные включают плотность и атомный номер, и при этом сокращение числа переменных включает превращение плотности в функцию атомного номера или превращение атомного номера в функцию плотности.
17. Система по п. 11, дополнительно содержащая код, который, будучи исполняемым, заставляет контроллер:
определять упрощение модели ослабления АМ на основании определения упрощения, которое минимизирует остаточную ошибку при выполнении итеративной реконструкции.
18. Система по п. 11, причем для выполнения итеративной реконструкции образца с использованием упрощенной модели прямой проекции уменьшаются артефакты увеличения жесткости излучения в реконструкции.
19. Система по п. 11, причем код, будучи исполняемым одним или более процессорами, заставляет контроллер получать высокоэнергетические данные образца при сканировании с низким разрешением и высокой энергией и дополнительно заставляет использовать высокоэнергетические данные для определения упрощения и начальных значений плотности и атомного номера, используемых в итеративной реконструкции.
20. Система по п. 11, причем в итеративной реконструкции используются полученные с высоким разрешением и низкой энергией данные образца.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201862659887P | 2018-04-19 | 2018-04-19 | |
US62/659,887 | 2018-04-19 |
Publications (3)
Publication Number | Publication Date |
---|---|
RU2019111292A RU2019111292A (ru) | 2020-10-16 |
RU2019111292A3 RU2019111292A3 (ru) | 2021-11-01 |
RU2770544C2 true RU2770544C2 (ru) | 2022-04-18 |
Family
ID=68237003
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU2019111292A RU2770544C2 (ru) | 2018-04-19 | 2019-04-16 | Коррекция увеличения жесткости рентгеновского излучения при томографической реконструкции с использованием модели ослабления альвареса-маковски |
Country Status (5)
Country | Link |
---|---|
US (1) | US10902650B2 (ru) |
JP (1) | JP7181832B2 (ru) |
CN (1) | CN110389138B (ru) |
AU (1) | AU2019202047B2 (ru) |
RU (1) | RU2770544C2 (ru) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10846860B2 (en) * | 2013-03-05 | 2020-11-24 | Nview Medical Inc. | Systems and methods for x-ray tomosynthesis image reconstruction |
JP2020537555A (ja) | 2017-09-22 | 2020-12-24 | エヌビュー メディカル インク | 機械学習レギュラライザを用いた画像再構成 |
US11610297B2 (en) * | 2019-12-02 | 2023-03-21 | Kla Corporation | Tomography based semiconductor measurements using simplified models |
FR3126779B1 (fr) * | 2021-09-09 | 2024-05-31 | Safran | Procédé, dispositif et programme d’ordinateur de contrôle d’une pièce par radiographie aux rayons X |
CN115078214B (zh) * | 2022-06-14 | 2024-01-26 | 西南石油大学 | 基于双能ct的岩心矿物识别及孔隙度分布联合表征方法 |
CN117330592B (zh) * | 2023-11-30 | 2024-03-15 | 赛诺威盛科技(北京)股份有限公司 | 能谱及散射校正方法、装置和设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030156684A1 (en) * | 2002-02-20 | 2003-08-21 | Fessler Jeffrey A. | Method for statistically reconstructing images from a plurality of transmission measurements having energy diversity and image reconstructor apparatus utilizing the method |
US20100008558A1 (en) * | 2006-12-15 | 2010-01-14 | Koninklijke Philips Electronics N. V. | Spectrally resolving x-ray imaging device |
US20100014737A1 (en) * | 2008-06-27 | 2010-01-21 | Ruehrnschopf Ernst-Peter | Method for generating image data relating to a virtually prescribable x-ray tube voltage from first and second CT image data |
US20160267688A1 (en) * | 2015-03-12 | 2016-09-15 | Kabushiki Kaisha Toshiba | Concurrent forward- and back-projection iterative reconstruction system |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4029963A (en) * | 1976-07-30 | 1977-06-14 | The Board Of Trustees Of Leland Stanford Junior University | X-ray spectral decomposition imaging system |
US7356113B2 (en) * | 2003-02-12 | 2008-04-08 | Brandeis University | Tomosynthesis imaging system and method |
US6950492B2 (en) * | 2003-06-25 | 2005-09-27 | Besson Guy M | Dynamic multi-spectral X-ray projection imaging |
DE102011075917B4 (de) * | 2011-05-16 | 2021-10-28 | Siemens Healthcare Gmbh | Verfahren zum Bereitstellen eines 3D-Bilddatensatzes mit unterdrückten Messfeldüberschreitungsartefakten und Computertomograph |
JP6177800B2 (ja) * | 2012-01-12 | 2017-08-09 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | X線システムにおける減衰画像データおよび位相画像データの生成 |
US8885975B2 (en) * | 2012-06-22 | 2014-11-11 | General Electric Company | Method and apparatus for iterative reconstruction |
US9495516B2 (en) * | 2012-11-21 | 2016-11-15 | The Trustees Of Columbia University In The City Of New York | Systems, methods, and devices for image reconstruction using combined PDE-constrained and simplified spherical harmonics algorithm |
US9600924B2 (en) * | 2014-02-05 | 2017-03-21 | Siemens Aktiengesellschaft | Iterative reconstruction of image data in CT |
EP2945123A1 (en) * | 2014-05-12 | 2015-11-18 | Agfa Healthcare | Computerized tomographic image exposure and reconstruction method |
KR20150145633A (ko) * | 2014-06-20 | 2015-12-30 | 삼성전자주식회사 | 엑스선 영상 장치 및 그 제어 방법 |
WO2016028655A1 (en) * | 2014-08-16 | 2016-02-25 | Fei Company | Correction of beam hardening artifacts in microtomography for samples imaged in containers |
-
2018
- 2018-08-23 US US16/110,264 patent/US10902650B2/en active Active
-
2019
- 2019-03-25 AU AU2019202047A patent/AU2019202047B2/en active Active
- 2019-04-15 JP JP2019077112A patent/JP7181832B2/ja active Active
- 2019-04-16 RU RU2019111292A patent/RU2770544C2/ru active
- 2019-04-18 CN CN201910315600.XA patent/CN110389138B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030156684A1 (en) * | 2002-02-20 | 2003-08-21 | Fessler Jeffrey A. | Method for statistically reconstructing images from a plurality of transmission measurements having energy diversity and image reconstructor apparatus utilizing the method |
US20100008558A1 (en) * | 2006-12-15 | 2010-01-14 | Koninklijke Philips Electronics N. V. | Spectrally resolving x-ray imaging device |
US20100014737A1 (en) * | 2008-06-27 | 2010-01-21 | Ruehrnschopf Ernst-Peter | Method for generating image data relating to a virtually prescribable x-ray tube voltage from first and second CT image data |
US20160267688A1 (en) * | 2015-03-12 | 2016-09-15 | Kabushiki Kaisha Toshiba | Concurrent forward- and back-projection iterative reconstruction system |
Also Published As
Publication number | Publication date |
---|---|
JP2019188135A (ja) | 2019-10-31 |
US20190325618A1 (en) | 2019-10-24 |
RU2019111292A (ru) | 2020-10-16 |
RU2019111292A3 (ru) | 2021-11-01 |
US10902650B2 (en) | 2021-01-26 |
JP7181832B2 (ja) | 2022-12-01 |
AU2019202047A1 (en) | 2019-11-07 |
CN110389138A (zh) | 2019-10-29 |
AU2019202047B2 (en) | 2023-10-05 |
CN110389138B (zh) | 2024-06-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
RU2770544C2 (ru) | Коррекция увеличения жесткости рентгеновского излучения при томографической реконструкции с использованием модели ослабления альвареса-маковски | |
US8923583B2 (en) | Methods and systems for performing model-based iterative reconstruction | |
JP5525750B2 (ja) | 二重エネルギーコンピュータ断層撮影における事前再構成分解及び校正を行うための方法及び装置 | |
EP3129957B1 (en) | Spectral x-ray imaging | |
US7760848B2 (en) | Method and system for generating a multi-spectral image of an object | |
US9036879B2 (en) | Multi-material decomposition using dual energy computed tomography | |
JP2018140165A (ja) | 医用画像生成装置 | |
JP2016536032A (ja) | 電子密度画像の連結再構成 | |
JP7043584B2 (ja) | 患者の臓器のx線撮影の方法 | |
JP2011221022A (ja) | 物体の組成を決定するシステム及び方法 | |
EP2702563B1 (en) | Multi-energy imaging | |
JP7079264B2 (ja) | コンピュータ断層撮影における単一エネルギーデータを用いたアーチファクト低減の方法 | |
WO2009102996A2 (en) | System and method for quantitative imaging of chemical composition to decompose more than two materials | |
Dhaene et al. | A realistic projection simulator for laboratory based X-ray micro-CT | |
JP7408664B2 (ja) | X線撮像システム | |
EP2652709B1 (en) | Imaging system for imaging a region of interest | |
EP3316785B1 (en) | Apparatus for multi material decomposition | |
CN107212898B (zh) | 图像重建方法 | |
US20210319600A1 (en) | System and Method for High Fidelity Computed Tomography | |
Lee et al. | A single scatter model for X-ray CT energy spectrum estimation and polychromatic reconstruction | |
WO2022268618A1 (en) | Multi-energy x-ray imaging with anatomical intelligence | |
Chen et al. | Line integral alternating minimization algorithm for dual-energy x-ray CT image reconstruction | |
Six et al. | Gauss-Newton-Krylov for reconstruction of polychromatic X-ray CT images | |
Zhao et al. | A novel iterative reconstruction method for dual-energy computed tomography based on polychromatic forward-projection calibration | |
Alsaffar et al. | Multi-Material Blind Beam Hardening Correction Based on Non-Linearity Adjustment of Projections |