Checkpoint
Pmc
Racine
Checkpoint
Pmc
Exploration
Accueil
Racine
Racine

Serveur d'exploration sur Mozart

Attention, ce site est en cours de développement !
Attention, site généré par des moyens informatiques à partir de corpus bruts.
Les informations ne sont donc pas validées.

Scattered Radiation Emission Imaging: Principles and Applications

Identifieur interne : 000317 ( Pmc/Checkpoint ); précédent : 000316; suivant : 000318

Scattered Radiation Emission Imaging: Principles and Applications

Auteurs : M. K. Nguyen [France] ; T. T. Truong [France] ; M. Morvidone [Argentine] ; H. Zaidi

Source :

RBID : PMC:3124682

Abstract

Imaging processes built on the Compton scattering effect have been under continuing investigation since it was first suggested in the 50s. However, despite many innovative contributions, there are still formidable theoretical and technical challenges to overcome. In this paper, we review the state-of-the-art principles of the so-called scattered radiation emission imaging. Basically, it consists of using the cleverly collected scattered radiation from a radiating object to reconstruct its inner structure. Image formation is based on the mathematical concept of compounded conical projection. It entails a Radon transform defined on circular cone surfaces in order to express the scattered radiation flux density on a detecting pixel. We discuss in particular invertible cases of such conical Radon transforms which form a mathematical basis for image reconstruction methods. Numerical simulations performed in two and three space dimensions speak in favor of the viability of this imaging principle and its potential applications in various fields.


Url:
DOI: 10.1155/2011/913893
PubMed: 21747823
PubMed Central: 3124682


Affiliations:

Links toward previous steps (curation, corpus...)

Links to Exploration step

PMC:3124682

Le document en format XML

<record>
<TEI>
<teiHeader>
<fileDesc>
<titleStmt>
<title xml:lang="en">Scattered Radiation Emission Imaging: Principles and Applications</title>
<author>
<name sortKey="Nguyen, M K" sort="Nguyen, M K" uniqKey="Nguyen M" first="M. K." last="Nguyen">M. K. Nguyen</name>
<affiliation wicri:level="3">
<nlm:aff id="I1">Laboratoire Equipes Traitement de l'Information et Systèmes, CNRS UMR 8051/ENSEA, Université de Cergy-Pontoise, 95302 Cergy-Pontoise, France</nlm:aff>
<country xml:lang="fr">France</country>
<wicri:regionArea>Laboratoire Equipes Traitement de l'Information et Systèmes, CNRS UMR 8051/ENSEA, Université de Cergy-Pontoise, 95302 Cergy-Pontoise</wicri:regionArea>
<placeName>
<region type="region" nuts="2">Île-de-France</region>
<settlement type="city">Cergy-Pontoise</settlement>
</placeName>
</affiliation>
</author>
<author>
<name sortKey="Truong, T T" sort="Truong, T T" uniqKey="Truong T" first="T. T." last="Truong">T. T. Truong</name>
<affiliation wicri:level="3">
<nlm:aff id="I2">Laboratoire de Physique Théorique et Modélisation, CNRS UMR 8089, Université de Cergy-Pontoise, 95302 Cergy-Pontoise, France</nlm:aff>
<country xml:lang="fr">France</country>
<wicri:regionArea>Laboratoire de Physique Théorique et Modélisation, CNRS UMR 8089, Université de Cergy-Pontoise, 95302 Cergy-Pontoise</wicri:regionArea>
<placeName>
<region type="region" nuts="2">Île-de-France</region>
<settlement type="city">Cergy-Pontoise</settlement>
</placeName>
</affiliation>
</author>
<author>
<name sortKey="Morvidone, M" sort="Morvidone, M" uniqKey="Morvidone M" first="M." last="Morvidone">M. Morvidone</name>
<affiliation wicri:level="1">
<nlm:aff id="I3">Facultad Regional Buenos Aires, Universidad Tecnológica Nacional, Mozart 2300, C1407IVT Buenos Aires, Argentina</nlm:aff>
<country xml:lang="fr">Argentine</country>
<wicri:regionArea>Facultad Regional Buenos Aires, Universidad Tecnológica Nacional, Mozart 2300, C1407IVT Buenos Aires</wicri:regionArea>
<wicri:noRegion>C1407IVT Buenos Aires</wicri:noRegion>
</affiliation>
</author>
<author>
<name sortKey="Zaidi, H" sort="Zaidi, H" uniqKey="Zaidi H" first="H." last="Zaidi">H. Zaidi</name>
<affiliation>
<nlm:aff>NONE</nlm:aff>
</affiliation>
<affiliation>
<nlm:aff>NONE</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt>
<idno type="wicri:source">PMC</idno>
<idno type="pmid">21747823</idno>
<idno type="pmc">3124682</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3124682</idno>
<idno type="RBID">PMC:3124682</idno>
<idno type="doi">10.1155/2011/913893</idno>
<date when="2011">2011</date>
<idno type="wicri:Area/Pmc/Corpus">000845</idno>
<idno type="wicri:Area/Pmc/Curation">000845</idno>
<idno type="wicri:Area/Pmc/Checkpoint">000317</idno>
</publicationStmt>
<sourceDesc>
<biblStruct>
<analytic>
<title xml:lang="en" level="a" type="main">Scattered Radiation Emission Imaging: Principles and Applications</title>
<author>
<name sortKey="Nguyen, M K" sort="Nguyen, M K" uniqKey="Nguyen M" first="M. K." last="Nguyen">M. K. Nguyen</name>
<affiliation wicri:level="3">
<nlm:aff id="I1">Laboratoire Equipes Traitement de l'Information et Systèmes, CNRS UMR 8051/ENSEA, Université de Cergy-Pontoise, 95302 Cergy-Pontoise, France</nlm:aff>
<country xml:lang="fr">France</country>
<wicri:regionArea>Laboratoire Equipes Traitement de l'Information et Systèmes, CNRS UMR 8051/ENSEA, Université de Cergy-Pontoise, 95302 Cergy-Pontoise</wicri:regionArea>
<placeName>
<region type="region" nuts="2">Île-de-France</region>
<settlement type="city">Cergy-Pontoise</settlement>
</placeName>
</affiliation>
</author>
<author>
<name sortKey="Truong, T T" sort="Truong, T T" uniqKey="Truong T" first="T. T." last="Truong">T. T. Truong</name>
<affiliation wicri:level="3">
<nlm:aff id="I2">Laboratoire de Physique Théorique et Modélisation, CNRS UMR 8089, Université de Cergy-Pontoise, 95302 Cergy-Pontoise, France</nlm:aff>
<country xml:lang="fr">France</country>
<wicri:regionArea>Laboratoire de Physique Théorique et Modélisation, CNRS UMR 8089, Université de Cergy-Pontoise, 95302 Cergy-Pontoise</wicri:regionArea>
<placeName>
<region type="region" nuts="2">Île-de-France</region>
<settlement type="city">Cergy-Pontoise</settlement>
</placeName>
</affiliation>
</author>
<author>
<name sortKey="Morvidone, M" sort="Morvidone, M" uniqKey="Morvidone M" first="M." last="Morvidone">M. Morvidone</name>
<affiliation wicri:level="1">
<nlm:aff id="I3">Facultad Regional Buenos Aires, Universidad Tecnológica Nacional, Mozart 2300, C1407IVT Buenos Aires, Argentina</nlm:aff>
<country xml:lang="fr">Argentine</country>
<wicri:regionArea>Facultad Regional Buenos Aires, Universidad Tecnológica Nacional, Mozart 2300, C1407IVT Buenos Aires</wicri:regionArea>
<wicri:noRegion>C1407IVT Buenos Aires</wicri:noRegion>
</affiliation>
</author>
<author>
<name sortKey="Zaidi, H" sort="Zaidi, H" uniqKey="Zaidi H" first="H." last="Zaidi">H. Zaidi</name>
<affiliation>
<nlm:aff>NONE</nlm:aff>
</affiliation>
<affiliation>
<nlm:aff>NONE</nlm:aff>
</affiliation>
</author>
</analytic>
<series>
<title level="j">International Journal of Biomedical Imaging</title>
<idno type="ISSN">1687-4188</idno>
<idno type="e-ISSN">1687-4196</idno>
<imprint>
<date when="2011">2011</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc>
<textClass></textClass>
</profileDesc>
</teiHeader>
<front>
<div type="abstract" xml:lang="en">
<p>Imaging processes built on the Compton scattering effect have been under continuing investigation since it was first suggested in the 50s. However, despite many innovative contributions, there are still formidable theoretical and technical challenges to overcome. In this paper, we review the state-of-the-art principles of the so-called scattered radiation emission imaging. Basically, it consists of using the cleverly collected scattered radiation from a radiating object to reconstruct its inner structure. Image formation is based on the mathematical concept of compounded conical projection. It entails a Radon transform defined on circular cone surfaces in order to express the scattered radiation flux density on a detecting pixel. We discuss in particular invertible cases of such conical Radon transforms which form a mathematical basis for image reconstruction methods. Numerical simulations performed in two and three space dimensions speak in favor of the viability of this imaging principle and its potential applications in various fields.</p>
</div>
</front>
<back>
<div1 type="bibliography">
<listBibl>
<biblStruct>
<analytic>
<author>
<name sortKey="Tuy, H" uniqKey="Tuy H">H Tuy</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Zhao, S" uniqKey="Zhao S">S Zhao</name>
</author>
<author>
<name sortKey="Yu, H" uniqKey="Yu H">H Yu</name>
</author>
<author>
<name sortKey="Wang, G" uniqKey="Wang G">G Wang</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Wernick, M" uniqKey="Wernick M">M Wernick</name>
</author>
<author>
<name sortKey="Aarsvold, J" uniqKey="Aarsvold J">J Aarsvold</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Cherry, Sr" uniqKey="Cherry S">SR Cherry</name>
</author>
<author>
<name sortKey="Sorenson, Ja" uniqKey="Sorenson J">JA Sorenson</name>
</author>
<author>
<name sortKey="Phelps, Me" uniqKey="Phelps M">ME Phelps</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Grangeat, P" uniqKey="Grangeat P">P Grangeat</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Zaidi, H" uniqKey="Zaidi H">H Zaidi</name>
</author>
<author>
<name sortKey="Koral, Kf" uniqKey="Koral K">KF Koral</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Lale, Pg" uniqKey="Lale P">PG Lale</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Lenti, M" uniqKey="Lenti M">M Lenti</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Guzzardi, R" uniqKey="Guzzardi R">R Guzzardi</name>
</author>
<author>
<name sortKey="Licitra, G" uniqKey="Licitra G">G Licitra</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Harding, G" uniqKey="Harding G">G Harding</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Norton, Sj" uniqKey="Norton S">SJ Norton</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Hussein, Ema" uniqKey="Hussein E">EMA Hussein</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Gautam, Sr" uniqKey="Gautam S">SR Gautam</name>
</author>
<author>
<name sortKey="Hopkins, Ff" uniqKey="Hopkins F">FF Hopkins</name>
</author>
<author>
<name sortKey="Klinksiek, R" uniqKey="Klinksiek R">R Klinksiek</name>
</author>
<author>
<name sortKey="Morgan, Il" uniqKey="Morgan I">IL Morgan</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Todd, Rw" uniqKey="Todd R">RW Todd</name>
</author>
<author>
<name sortKey="Nightingale, J" uniqKey="Nightingale J">J Nightingale</name>
</author>
<author>
<name sortKey="Everett, D" uniqKey="Everett D">D Everett</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Singh, M" uniqKey="Singh M">M Singh</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Singh, M" uniqKey="Singh M">M Singh</name>
</author>
<author>
<name sortKey="Doria, D" uniqKey="Doria D">D Doria</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Gerl, J" uniqKey="Gerl J">J Gerl</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Braem, A" uniqKey="Braem A">A Braem</name>
</author>
<author>
<name sortKey="Llatas, Mchamizo" uniqKey="Llatas M">MChamizo Llatas</name>
</author>
<author>
<name sortKey="Chesi, E" uniqKey="Chesi E">E Chesi</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Delarbre, Jl" uniqKey="Delarbre J">JL Delarbre</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Nguyen, Mk" uniqKey="Nguyen M">MK Nguyen</name>
</author>
<author>
<name sortKey="Truong, Tt" uniqKey="Truong T">TT Truong</name>
</author>
<author>
<name sortKey="Driol, C" uniqKey="Driol C">C Driol</name>
</author>
<author>
<name sortKey="Zaidi, H" uniqKey="Zaidi H">H Zaidi</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Zeng, Gl" uniqKey="Zeng G">GL Zeng</name>
</author>
<author>
<name sortKey="Gullberg, Gt" uniqKey="Gullberg G">GT Gullberg</name>
</author>
<author>
<name sortKey="Tsui, Bmw" uniqKey="Tsui B">BMW Tsui</name>
</author>
<author>
<name sortKey="Terry, Ja" uniqKey="Terry J">JA Terry</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Zeng, Gl" uniqKey="Zeng G">GL Zeng</name>
</author>
<author>
<name sortKey="Clack, R" uniqKey="Clack R">R Clack</name>
</author>
<author>
<name sortKey="Gullberg, Gt" uniqKey="Gullberg G">GT Gullberg</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Kadrmas, Dj" uniqKey="Kadrmas D">DJ Kadrmas</name>
</author>
<author>
<name sortKey="Frey, Ec" uniqKey="Frey E">EC Frey</name>
</author>
<author>
<name sortKey="Karimi, Ss" uniqKey="Karimi S">SS Karimi</name>
</author>
<author>
<name sortKey="Tsui, Bmw" uniqKey="Tsui B">BMW Tsui</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Pack, Jd" uniqKey="Pack J">JD Pack</name>
</author>
<author>
<name sortKey="Noo, F" uniqKey="Noo F">F Noo</name>
</author>
<author>
<name sortKey="Clackdoyle, R" uniqKey="Clackdoyle R">R Clackdoyle</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Barrett, Hh" uniqKey="Barrett H">HH Barrett</name>
</author>
<author>
<name sortKey="Swindell, W" uniqKey="Swindell W">W Swindell</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Floyd, Ce" uniqKey="Floyd C">CE Floyd</name>
</author>
<author>
<name sortKey="Jaszczak, Rj" uniqKey="Jaszczak R">RJ Jaszczak</name>
</author>
<author>
<name sortKey="Greer, Kl" uniqKey="Greer K">KL Greer</name>
</author>
<author>
<name sortKey="Coleman, Re" uniqKey="Coleman R">RE Coleman</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Gunter, Dl" uniqKey="Gunter D">DL Gunter</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Kadrmas, Dj" uniqKey="Kadrmas D">DJ Kadrmas</name>
</author>
<author>
<name sortKey="Frey, Ec" uniqKey="Frey E">EC Frey</name>
</author>
<author>
<name sortKey="Tsui, Bmw" uniqKey="Tsui B">BMW Tsui</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Nguyen, Mk" uniqKey="Nguyen M">MK Nguyen</name>
</author>
<author>
<name sortKey="Faye, C" uniqKey="Faye C">C Faye</name>
</author>
<author>
<name sortKey="Eglin, L" uniqKey="Eglin L">L Eglin</name>
</author>
<author>
<name sortKey="Truong, Tt" uniqKey="Truong T">TT Truong</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Eglin, L" uniqKey="Eglin L">L Eglin</name>
</author>
<author>
<name sortKey="Nguyen, Mk" uniqKey="Nguyen M">MK Nguyen</name>
</author>
<author>
<name sortKey="Faye, C" uniqKey="Faye C">C Faye</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Evans, Rd" uniqKey="Evans R">RD Evans</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Nguyen, Mk" uniqKey="Nguyen M">MK Nguyen</name>
</author>
<author>
<name sortKey="Truong, Tt" uniqKey="Truong T">TT Truong</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Hertle, A" uniqKey="Hertle A">A Hertle</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Zaidi, H" uniqKey="Zaidi H">H Zaidi</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Nguyen, Mk" uniqKey="Nguyen M">MK Nguyen</name>
</author>
<author>
<name sortKey="Truong, Tt" uniqKey="Truong T">TT Truong</name>
</author>
<author>
<name sortKey="Bui, Hd" uniqKey="Bui H">HD Bui</name>
</author>
<author>
<name sortKey="Delarbre, Jl" uniqKey="Delarbre J">JL Delarbre</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Siegman, Ae" uniqKey="Siegman A">AE Siegman</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Nguyen, Mk" uniqKey="Nguyen M">MK Nguyen</name>
</author>
<author>
<name sortKey="Truong, Tt" uniqKey="Truong T">TT Truong</name>
</author>
<author>
<name sortKey="Delarbre, Jl" uniqKey="Delarbre J">JL Delarbre</name>
</author>
<author>
<name sortKey="Roux, C" uniqKey="Roux C">C Roux</name>
</author>
<author>
<name sortKey="Zaidi, H" uniqKey="Zaidi H">H Zaidi</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Driol, C" uniqKey="Driol C">C Driol</name>
</author>
<author>
<name sortKey="Nguyen, Mk" uniqKey="Nguyen M">MK Nguyen</name>
</author>
<author>
<name sortKey="Truong, Tt" uniqKey="Truong T">TT Truong</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Nguyen, Mk" uniqKey="Nguyen M">MK Nguyen</name>
</author>
<author>
<name sortKey="Driol, C" uniqKey="Driol C">C Driol</name>
</author>
<author>
<name sortKey="Truong, Tt" uniqKey="Truong T">TT Truong</name>
</author>
<author>
<name sortKey="Zaidi, H" uniqKey="Zaidi H">H Zaidi</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Truong, Tt" uniqKey="Truong T">TT Truong</name>
</author>
<author>
<name sortKey="Nguyen, Mk" uniqKey="Nguyen M">MK Nguyen</name>
</author>
<author>
<name sortKey="Zaidi, H" uniqKey="Zaidi H">H Zaidi</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Truong, Tt" uniqKey="Truong T">TT Truong</name>
</author>
<author>
<name sortKey="Nguyen, Mk" uniqKey="Nguyen M">MK Nguyen</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Morvidone, M" uniqKey="Morvidone M">M Morvidone</name>
</author>
<author>
<name sortKey="Nguyen, Mk" uniqKey="Nguyen M">MK Nguyen</name>
</author>
<author>
<name sortKey="Truong, Tt" uniqKey="Truong T">TT Truong</name>
</author>
<author>
<name sortKey="Zaidi, H" uniqKey="Zaidi H">H Zaidi</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Lavoine, J" uniqKey="Lavoine J">J Lavoine</name>
</author>
</analytic>
</biblStruct>
</listBibl>
</div1>
</back>
</TEI>
<pmc article-type="review-article">
<pmc-dir>properties open_access</pmc-dir>
<front>
<journal-meta>
<journal-id journal-id-type="nlm-ta">Int J Biomed Imaging</journal-id>
<journal-id journal-id-type="publisher-id">IJBI</journal-id>
<journal-title-group>
<journal-title>International Journal of Biomedical Imaging</journal-title>
</journal-title-group>
<issn pub-type="ppub">1687-4188</issn>
<issn pub-type="epub">1687-4196</issn>
<publisher>
<publisher-name>Hindawi Publishing Corporation</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="pmid">21747823</article-id>
<article-id pub-id-type="pmc">3124682</article-id>
<article-id pub-id-type="doi">10.1155/2011/913893</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Review Article</subject>
</subj-group>
</article-categories>
<title-group>
<article-title>Scattered Radiation Emission Imaging: Principles and Applications</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Nguyen</surname>
<given-names>M. K.</given-names>
</name>
<xref ref-type="aff" rid="I1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Truong</surname>
<given-names>T. T.</given-names>
</name>
<xref ref-type="aff" rid="I2">
<sup>2</sup>
</xref>
<xref ref-type="corresp" rid="cor1">*</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Morvidone</surname>
<given-names>M.</given-names>
</name>
<xref ref-type="aff" rid="I3">
<sup>3</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zaidi</surname>
<given-names>H.</given-names>
</name>
<xref ref-type="aff" rid="I4">
<sup>4, 5</sup>
</xref>
<xref ref-type="aff" rid="I5"></xref>
</contrib>
</contrib-group>
<aff id="I1">
<sup>1</sup>
Laboratoire Equipes Traitement de l'Information et Systèmes, CNRS UMR 8051/ENSEA, Université de Cergy-Pontoise, 95302 Cergy-Pontoise, France</aff>
<aff id="I2">
<sup>2</sup>
Laboratoire de Physique Théorique et Modélisation, CNRS UMR 8089, Université de Cergy-Pontoise, 95302 Cergy-Pontoise, France</aff>
<aff id="I3">
<sup>3</sup>
Facultad Regional Buenos Aires, Universidad Tecnológica Nacional, Mozart 2300, C1407IVT Buenos Aires, Argentina</aff>
<aff id="I4">
<sup>4</sup>
Division of Nuclear Medicine, Geneva University Hospital, 1211 Geneva 4, Switzerland</aff>
<aff id="I5">
<sup>5</sup>
Geneva Neuroscience Center, Geneva University, 1211 Geneva 4, Switzerland</aff>
<author-notes>
<corresp id="cor1">*T. T. Truong:
<email>truong@u-cergy.fr</email>
</corresp>
<fn fn-type="other">
<p>Academic Editor: Peter Bruyndonckx</p>
</fn>
</author-notes>
<pub-date pub-type="ppub">
<year>2011</year>
</pub-date>
<pub-date pub-type="epub">
<day>7</day>
<month>6</month>
<year>2011</year>
</pub-date>
<volume>2011</volume>
<elocation-id>913893</elocation-id>
<history>
<date date-type="received">
<day>1</day>
<month>12</month>
<year>2010</year>
</date>
<date date-type="rev-recd">
<day>8</day>
<month>3</month>
<year>2011</year>
</date>
<date date-type="accepted">
<day>10</day>
<month>4</month>
<year>2011</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright © 2011 M. K. Nguyen et al.</copyright-statement>
<copyright-year>2011</copyright-year>
<license license-type="open-access">
<license-p>This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</license-p>
</license>
</permissions>
<abstract>
<p>Imaging processes built on the Compton scattering effect have been under continuing investigation since it was first suggested in the 50s. However, despite many innovative contributions, there are still formidable theoretical and technical challenges to overcome. In this paper, we review the state-of-the-art principles of the so-called scattered radiation emission imaging. Basically, it consists of using the cleverly collected scattered radiation from a radiating object to reconstruct its inner structure. Image formation is based on the mathematical concept of compounded conical projection. It entails a Radon transform defined on circular cone surfaces in order to express the scattered radiation flux density on a detecting pixel. We discuss in particular invertible cases of such conical Radon transforms which form a mathematical basis for image reconstruction methods. Numerical simulations performed in two and three space dimensions speak in favor of the viability of this imaging principle and its potential applications in various fields.</p>
</abstract>
</article-meta>
</front>
<body>
<sec sec-type="section" id="sec1">
<title>1. Introduction</title>
<p>Since the early 50s, ionizing radiations (in particular gamma-rays), because of their penetrating property, have been used to explore the interior of objects. At first, this was done in transmission mode with an external radiation source, which projects a shadow onto a plane detector. Later, it was shown that a three-dimensional image can be reconstructed provided there is a sufficient number of such two-dimensional projections generated by the displacement of the source/detector pair in space. This three-dimensional image reconstruction of the inner object structure relies on the existence of the inverse of the so-called X-ray transform, which correctly models the above process of data acquisition [
<xref ref-type="bibr" rid="B1">1</xref>
,
<xref ref-type="bibr" rid="B2">2</xref>
].</p>
<p>A second modality, called emission imaging, deals with radiation-emitting objects. In nuclear medicine, this modality concerns objects (human organs), which, after injection of a radiotracer, displays its biodistribution in the human body. Nowadays, the image reconstruction of both single-photon emitting and positron-emitting tracer distributions is achieved by single-photon emission computed tomography (SPECT) [
<xref ref-type="bibr" rid="B3">3</xref>
] and positron emission tomography (PET) [
<xref ref-type="bibr" rid="B4">4</xref>
,
<xref ref-type="bibr" rid="B5">5</xref>
], respectively. These two modalities are based on the invertibility of the standard Radon transform.</p>
<p>However, due to the interaction of radiation with matter, gamma-ray imaging is plagued by Compton scatter, which degrades image quality and spatial resolution. Thus, effects due to scattered photons should at best be eliminated or at least be reduced [
<xref ref-type="bibr" rid="B6">6</xref>
]. However, a more astute point of view would be to take advantage of their properties either for improving image quality or for generating new imaging processes. The idea of Compton scatter imaging has been launched many years ago and many ways to exploit Compton scattering for imaging purposes have been introduced.</p>
<p>An early proposal goes back to the 50s [
<xref ref-type="bibr" rid="B7">7</xref>
], but interest in this concept has remained vivid [
<xref ref-type="bibr" rid="B8">8</xref>
], because this idea has many highly desirable features. In the field of diagnostic medical imaging, radiography using scattered radiation could provide a direct and quantitative measurement of the density of the studied object. In nondestructive testing, it offers three advantages.</p>
<list list-type="roman-lower">
<list-item>
<p>It permits to place both the radiation source and the detector on the same side of the object.</p>
</list-item>
<list-item>
<p>It has also greater sensitivity to low-density materials such as gases.</p>
</list-item>
<list-item>
<p>Finally, it allows direct spatial definition with high-contrast resolution.</p>
</list-item>
</list>
<p>With the advent of X-ray computed tomography (CT), interest in Compton scatter imaging has waned for a while. But research in this field has remained very much alive, and a large variety of imaging techniques have been developed [
<xref ref-type="bibr" rid="B3">3</xref>
,
<xref ref-type="bibr" rid="B9">9</xref>
,
<xref ref-type="bibr" rid="B10">10</xref>
]. </p>
<p>Earlier modalities for Compton scatter imaging are classified according to the way measurement of the spatial distribution of scattered radiation is done or the number of simultaneous volume elements being scanned: that is, point by point, line by line, or plane by plane (see reviews [
<xref ref-type="bibr" rid="B9">9</xref>
,
<xref ref-type="bibr" rid="B10">10</xref>
]). Most of the devices work at constant scattering angles (generally at 90 degrees). In the mid-90s, the concept of Compton scatter tomography was introduced by Norton [
<xref ref-type="bibr" rid="B11">11</xref>
] and subsequently developed by many other workers [
<xref ref-type="bibr" rid="B12">12</xref>
,
<xref ref-type="bibr" rid="B13">13</xref>
]. A prominent example in which Compton scattering acts as imaging agent without mechanical collimation is the co-called Compton camera [
<xref ref-type="bibr" rid="B14">14</xref>
–
<xref ref-type="bibr" rid="B16">16</xref>
], as well as gamma-ray tracking imaging or the like. More recently, Compton scatter imaging using annihilation pair photons with coincidence measurements has appeared on the scene as a yet unexploited imaging technique [
<xref ref-type="bibr" rid="B17">17</xref>
]. Related concepts allow enhancing the detection efficiency by reconstructing a significant fraction of events which underwent Compton scattering in crystals [
<xref ref-type="bibr" rid="B18">18</xref>
].</p>
<p>In this work, we review a different approach to Compton-scattered radiation for
<italic> emission</italic>
imaging. We first point out why scattered photons should be used instead of primary photons in gamma-ray imaging. Then, we show how the image formation process is actually performed using the concept of compounded conical projection. The basic mathematical object to be considered is a generalized Radon transform on circular cone surfaces.
<xref ref-type="sec" rid="sec2"> Section 2</xref>
describes the general principles of this imaging modality and its properties. Then, in consecutive subsequent sections, we treat explicitly a special case, involving a standard gamma camera-based SPECT system, in two and three dimensions. Results on numerical simulations are also presented; they speak in favor for the realizability of this imaging method, once real world issues (noise, energy uncertainty, sampling, etc.) are addressed and resolved. A conclusion summarizes this paper and points out future research perspectives. </p>
</sec>
<sec sec-type="section" id="sec2">
<title>2. The Emission Imaging Problem</title>
<p>The aim is to image the interior of radiating (or made radiating) objects, which emit gamma-rays of a given energy
<italic>E</italic>
<sub>0</sub>
resulting from a nuclear transition, (e.g., Technetium 99 m in nuclear medicine with
<italic>E</italic>
<sub>0</sub>
= 140 keV). To intercept the emitted gamma photons, we will use a device called a
<italic>directional detector</italic>
(or DD for short). This is a point-like pixel at a site
<bold>D</bold>
in space, capable of absorbing gamma photons of any energy
<italic>E</italic>
below
<italic>E</italic>
<sub>0</sub>
, in a given direction specified by a unit vector
<bold>n</bold>
. It may be thought as a one-pixel collimated gamma camera. This device is only used for argumentation.</p>
<p>For each pair (
<bold>D</bold>
,
<bold>n</bold>
), the DD records a photon flux density coming directly from the radiating object reaching the pixel site
<bold>D</bold>
in the direction −
<bold>n</bold>
. This measurement is called in the specialized jargon a
<italic>projection</italic>
, or better a
<italic>linear projection</italic>
, since it is done along a straight line passing through
<bold>D</bold>
in the direction
<bold>n</bold>
. As the recorded photons carry the original energy
<italic>E</italic>
<sub>0</sub>
, they are called primary (or nonscattered) photons. In this way, the set of such measurements constitutes the conventional X-ray projection data of the object radioactivity function, which is represented by a nonnegative, well-behaved, and compactly supported function
<italic>f</italic>
(
<bold>r</bold>
) =
<italic>f</italic>
(
<italic>x</italic>
,
<italic>y</italic>
,
<italic>z</italic>
). Mathematically, this data is represented by
<italic>Xf</italic>
(
<bold>D</bold>
,
<bold>n</bold>
), the so-called X-ray transform of
<italic>f</italic>
(
<italic>x</italic>
,
<italic>y</italic>
,
<italic>z</italic>
).
<xref ref-type="fig" rid="fig1"> Figure 1</xref>
illustrates this measurement concept.</p>
<p>As the X-ray transform is invertible, albeit under conditions,
<italic>f</italic>
(
<bold>r</bold>
) can be obtained by various reconstruction algorithms (see, e.g., among the many references [
<xref ref-type="bibr" rid="B19">21</xref>
–
<xref ref-type="bibr" rid="B22">24</xref>
]). </p>
<sec sec-type="subsection" id="sec2.1">
<title>2.1. Scattering</title>
<p>Photon transport through matter suffers from two lossy phenomena: photoelectric absorption and Compton scattering (
<xref ref-type="fig" rid="fig2">Figure 2</xref>
). The net depletion in photon number is described by a macroscopic coefficient of linear attenuation of the traveling photon flux density. However, in the energy range of a few hundred keV, Compton scattering is dominant and photoelectric absorption is negligible in biological tissues. Scattered photons have a disturbing effect in nuclear medicine imaging, see, for example, [
<xref ref-type="bibr" rid="B23">25</xref>
]. It causes blurring, loss of contrast, and false detection of emitting sources in the context of primary (nonscattered) radiation imaging. So it is natural to raise the question: is scattered radiation any good for imaging?</p>
<p>Handling away scatter in gamma-ray imaging has been pioneered by many authors. Algorithms to compensate for Compton scattering in SPECT imaging have been developed, for example, [
<xref ref-type="bibr" rid="B24">26</xref>
], and techniques allowing the determination of source depth via scattered radiation proposed [
<xref ref-type="bibr" rid="B25">27</xref>
]. The idea of using scattered photons to reduce the noise level of SPECT images has emerged in [
<xref ref-type="bibr" rid="B26">28</xref>
], in which data acquired in the photopeak and various scatter energy windows are statistically assembled to improve image quality. In 2001, it was observed that scattered radiation images of an object may be sorted out at a given energy (or at a given wavelength) using standard gamma camera data operating in list mode [
<xref ref-type="bibr" rid="B27">29</xref>
]. A series of apparent images labeled by the photon-scattered energy of the object is then acquired [
<xref ref-type="bibr" rid="B28">30</xref>
]. Subsequently, they are taken into account in the process of image restoration. Small details, unresolved before, emerge clearly separated from each other. As an example in bone scintigraphy,
<xref ref-type="fig" rid="fig3">Figure 3</xref>
shows that small hot spots or nodules, which are invisible on the left image, become perfectly distinguishable on the right image after accounting for scattered radiation. This is, in fact, very valuable for clinical diagnosis, assessment of response to treatment, and radiation therapy treatment planning. Thus, the newly revealed resolving power brought by scattered radiation has appeared very attractive for further development. </p>
<p>Therefore, it is of interest to take a close look at the effect of scattering in photon detection. Let us consider first the case of a point source emitting a monochromatic red light by a clear day. A human eye, placed at a certain distance from this source, would see only a red spot. However, if a fog cloud sets in, then the eye would see a diffuse red cloud much larger than the red spot. The fog cloud has made itself visible because the emitted red light is scattered by the fog droplets and re-emitted as scattered light by the fog droplets acting as scattering centers. This fog cloud has become a kind of secondary radiating object, visible to the human eye. It also implies the existence of a concealed red light source.</p>
<p>In the gamma energy range, a similar observation can be made. A single gamma-ray emitting point source behaves exactly in the same way with respect to our directional detector (DD), which is a kind of gamma-ray sensitive “eye.” But if the gamma-ray emitting point source is embedded in a medium of finite volume, which plays the role of the fog cloud, then light wave scattering by water droplets is replaced by Compton scattering of emitted photons with electric charges of the surrounding medium. If visible light emerges from scattering without changing its wavelength (
<italic>Rayleigh scattering</italic>
), the emerging scattered gamma ray has an energy
<italic>E</italic>
lower than the incident energy
<italic>E</italic>
<sub>0</sub>
, because part of
<italic>E</italic>
<sub>0</sub>
is transferred to electric charges in the traversed medium. As
<italic>E</italic>
is continuously distributed, the gamma-ray sensitive eye (DD) would now “see” a
<italic>red-shifted “polychromatic” radiation</italic>
emanating from the scattering medium volume [
<xref ref-type="bibr" rid="B30">20</xref>
]. The wavelengths
<italic>λ</italic>
of the scattered gamma-rays are longer than the incident wavelength
<italic>λ</italic>
<sub>0</sub>
, as given by the Compton formula [
<xref ref-type="bibr" rid="B29">31</xref>
],</p>
<p>
<disp-formula id="EEq1">
<label>(1)</label>
<mml:math id="M1">
<mml:mi>λ</mml:mi>
<mml:mo>=</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>λ</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mi>  </mml:mi>
<mml:mo>+</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>h</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>m</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mi>  </mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>−</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mi>cos</mml:mi>
<mml:mo> </mml:mo>
</mml:mrow>
<mml:mo> </mml:mo>
<mml:mrow>
<mml:mi>ω</mml:mi>
</mml:mrow>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:math>
</disp-formula>
where
<italic>m</italic>
<sub>
<italic>e</italic>
</sub>
is the electron mass,
<italic>h</italic>
is the Planck constant,
<italic>c</italic>
is the speed of light in vacuum, and
<italic>ω</italic>
is the Compton scattering angle.</p>
<p>Now instead of having a single-point source, consider a nonuniform three-dimensional distribution of gamma-ray emitting point sources embedded in a medium of finite volume. The question is what flux density of photons of energy
<italic>E</italic>
<
<italic>E</italic>
<sub>0</sub>
would a DD record at a site
<bold>D</bold>
in a direction
<bold>n</bold>
? To give an answer, we should follow the backward track taken by the gamma-ray before it reaches the DD.</p>
<p>A gamma photon arriving at
<bold>D</bold>
with energy
<italic>E</italic>
must have gone through a scattering with some electric charge at a site
<bold>N</bold>
situated on the straight line, which starts from
<bold>D</bold>
in the direction of
<bold>n</bold>
, see
<xref ref-type="fig" rid="fig4">Figure 4</xref>
. As the scattering angle
<italic>ω</italic>
is given by the Compton formula </p>
<p>
<disp-formula id="EEq2">
<label>(2)</label>
<mml:math id="M2">
<mml:mi>E</mml:mi>
<mml:mo>=</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>E</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mi>  </mml:mi>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>−</mml:mo>
<mml:mi>ε</mml:mi>
<mml:mi>  </mml:mi>
<mml:mi>cos</mml:mi>
<mml:mo> </mml:mo>
<mml:mi>ω</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>,</mml:mo>
</mml:math>
</disp-formula>
where
<italic>ε</italic>
=
<italic>E</italic>
<sub>0</sub>
/
<italic>m</italic>
<sub>
<italic>e</italic>
</sub>
  
<italic>c</italic>
<sup>2</sup>
, this photon must have originated from a site
<bold>S</bold>
, located on the surface of a circular cone of vertex
<bold>N</bold>
, axis
<italic>DN</italic>
, and opening angle
<italic>ω</italic>
. Following this picture, we can write down the photon flux density detected at
<bold>D</bold>
. </p>
<p>To this end, we need the expression of the Compton differential cross-section, which reads as the product of the electron density
<italic>n</italic>
<sub>
<italic>e</italic>
</sub>
(
<bold>N</bold>
) at the scattering site
<bold>N</bold>
by the so-called Klein-Nishina probability
<italic>𝒫</italic>
(
<italic>ω</italic>
) [
<xref ref-type="bibr" rid="B23">25</xref>
,
<xref ref-type="bibr" rid="B29">31</xref>
], that is,
<italic>n</italic>
<sub>
<italic>e</italic>
</sub>
(
<bold>N</bold>
)
<italic>𝒫</italic>
(
<italic>ω</italic>
), </p>
<p>
<disp-formula id="EEq3">
<label>(3)</label>
<mml:math id="M3">
<mml:mtable>
<mml:mtr>
<mml:mtd columnalign="right">
<mml:mi>𝒫</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>ω</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mtd>
<mml:mtd columnalign="left">
<mml:mo>=</mml:mo>
<mml:mi>π</mml:mi>
<mml:msubsup>
<mml:mrow>
<mml:mi>r</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>e</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mi>  </mml:mi>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>π</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>ε</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>−</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mi>cos</mml:mi>
<mml:mo> </mml:mo>
</mml:mrow>
<mml:mo> </mml:mo>
<mml:mrow>
<mml:mi>ω</mml:mi>
</mml:mrow>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd columnalign="right"></mml:mtd>
<mml:mtd columnalign="left">
<mml:mo> </mml:mo>
<mml:mo>×</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mi>cos</mml:mi>
<mml:mo> </mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mi>ω</mml:mi>
<mml:mo>+</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>−</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mi>cos</mml:mi>
<mml:mo> </mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mi>ω</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>ε</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>−</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mi>cos</mml:mi>
<mml:mo> </mml:mo>
</mml:mrow>
<mml:mo> </mml:mo>
<mml:mrow>
<mml:mi>ω</mml:mi>
</mml:mrow>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</disp-formula>
Radiation emitted at point source
<bold>S</bold>
is assumed to be isotropic. It propagates then to scattering site
<bold>N</bold>
and reaches the detection site
<bold>D</bold>
at the end. The incoming photon flux density on scattering site
<bold>N</bold>
is computed from the emission data at a point source. Let
<italic>f</italic>
(
<bold>S</bold>
)  
<italic>d </italic>
<bold>S</bold>
be the number of gamma photons emitted per unit of time by a volume
<italic>d </italic>
<bold>S</bold>
in the object around site
<bold>S</bold>
. The emission being isotropic, the number of photons emitted in the direction
<inline-formula>
<mml:math id="M4">
<mml:mrow>
<mml:mover accent="true">
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:mi>N</mml:mi>
</mml:mrow>
<mml:mo>⃗</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</inline-formula>
in a solid angle
<italic>dΩ</italic>
<sub>
<italic>S</italic>
</sub>
is </p>
<p>
<disp-formula id="eq1">
<label>(4)</label>
<mml:math id="M5">
<mml:mfrac>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold">S</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi mathvariant="bold">S</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>4</mml:mn>
<mml:mi>π</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mi>  </mml:mi>
<mml:mi>d</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi>Ω</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>.</mml:mo>
</mml:math>
</disp-formula>
Therefore, the incoming photon flux density at scattering site
<bold>N</bold>
is given by</p>
<p>
<disp-formula id="eq2">
<label>(5)</label>
<mml:math id="M6">
<mml:mfrac>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold">S</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi mathvariant="bold">S</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>4</mml:mn>
<mml:mi>π</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mi>  </mml:mi>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mi>e</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>−</mml:mo>
<mml:mi>𝒜</mml:mi>
<mml:mo minsize="0.75em" maxsize="0.75em">(</mml:mo>
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:mi>N</mml:mi>
</mml:mrow>
<mml:mo minsize="0.75em" maxsize="0.75em">)</mml:mo>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:msup>
<mml:mrow>
<mml:mi>N</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mfrac>
<mml:mo>,</mml:mo>
</mml:math>
</disp-formula>
where
<inline-formula>
<mml:math id="M7">
<mml:mi>S</mml:mi>
<mml:mi>N</mml:mi>
<mml:mo>=</mml:mo>
<mml:mo stretchy="false">|</mml:mo>
<mml:mover accent="true">
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:mi>N</mml:mi>
</mml:mrow>
<mml:mo>⃗</mml:mo>
</mml:mover>
<mml:mo stretchy="false">|</mml:mo>
</mml:math>
</inline-formula>
and
<italic>e</italic>
<sup>−
<italic>𝒜</italic>
(
<italic>SN</italic>
)</sup>
is the attenuation factor along the path
<italic>SN</italic>
and is given by the integral </p>
<p>
<disp-formula id="EEq4">
<label>(6)</label>
<mml:math id="M8">
<mml:mi>𝒜</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:mi>N</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>∫</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi mathvariant="bold">S</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="bold">N</mml:mi>
</mml:mrow>
</mml:munderover>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>μ</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold">S</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>s</mml:mi>
<mml:mi mathvariant="bold">k</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:math>
</disp-formula>
where
<inline-formula>
<mml:math id="M9">
<mml:mi mathvariant="bold">k</mml:mi>
<mml:mo>=</mml:mo>
<mml:mover accent="true">
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:mi>N</mml:mi>
</mml:mrow>
<mml:mo>⃗</mml:mo>
</mml:mover>
<mml:mo>/</mml:mo>
<mml:mi>S</mml:mi>
<mml:mi>N</mml:mi>
</mml:math>
</inline-formula>
and
<italic>μ</italic>
(
<bold>M</bold>
) is the linear coefficient of absorption. Recall that, by assumption,
<italic>μ</italic>
is equal to the product of the Compton scattering cross-section and the electron density. </p>
<p>Next, the number of scatterers around site
<bold>N</bold>
in a volume
<italic>d </italic>
<bold>N</bold>
is
<italic>n</italic>
<sub>
<italic>e</italic>
</sub>
(
<bold>N</bold>
)  
<italic>d </italic>
<bold>N</bold>
. The net number of photons emerging from the scattering in an elementary solid angle
<italic>dΩ</italic>
<sub>
<italic>N</italic>
</sub>
is </p>
<p>
<disp-formula id="eq3">
<label>(7)</label>
<mml:math id="M10">
<mml:mfrac>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold">S</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi mathvariant="bold">S</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>4</mml:mn>
<mml:mi>π</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:msup>
<mml:mrow>
<mml:mi>N</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mfrac>
<mml:msub>
<mml:mrow>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold">N</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi mathvariant="bold">N</mml:mi>
<mml:mi>π</mml:mi>
<mml:msubsup>
<mml:mrow>
<mml:mi>r</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>e</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mi>𝒫</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>ω</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>d</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi>Ω</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>N</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>.</mml:mo>
</mml:math>
</disp-formula>
This means in turn that the detected photon flux density at site
<bold>D</bold>
is </p>
<p>
<disp-formula id="eq4">
<label>(8)</label>
<mml:math id="M11">
<mml:mfrac>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold">S</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi mathvariant="bold">S</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>4</mml:mn>
<mml:mi>π</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mi>  </mml:mi>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mi>e</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>−</mml:mo>
<mml:mi>𝒜</mml:mi>
<mml:mo minsize="0.75em" maxsize="0.75em">(</mml:mo>
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:mi>N</mml:mi>
</mml:mrow>
<mml:mo minsize="0.75em" maxsize="0.75em">)</mml:mo>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:msup>
<mml:mrow>
<mml:mi>N</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mfrac>
<mml:msub>
<mml:mrow>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold">N</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi mathvariant="bold">N</mml:mi>
<mml:mi>π</mml:mi>
<mml:msubsup>
<mml:mrow>
<mml:mi>r</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>e</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mi>𝒫</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>ω</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mi>e</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>−</mml:mo>
<mml:mi>𝒜</mml:mi>
<mml:mo minsize="0.75em" maxsize="0.75em">(</mml:mo>
<mml:mrow>
<mml:mi>N</mml:mi>
<mml:mi>D</mml:mi>
</mml:mrow>
<mml:mo minsize="0.75em" maxsize="0.75em">)</mml:mo>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:mi>N</mml:mi>
<mml:msup>
<mml:mrow>
<mml:mi>D</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mfrac>
<mml:mo>.</mml:mo>
</mml:math>
</disp-formula>
</p>
<p>Now, all the contributing point sources
<bold>S</bold>
, for a given scattering center
<bold>N</bold>
, lie on a circular cone sheet of axis identified with
<inline-formula>
<mml:math id="M12">
<mml:mover accent="true">
<mml:mrow>
<mml:mi>D</mml:mi>
<mml:mi>N</mml:mi>
</mml:mrow>
<mml:mo>⃗</mml:mo>
</mml:mover>
<mml:mo>=</mml:mo>
<mml:mi mathvariant="bold">n</mml:mi>
<mml:mi>  </mml:mi>
<mml:mi>D</mml:mi>
<mml:mi>N</mml:mi>
</mml:math>
</inline-formula>
and opening angle
<italic>ω</italic>
. Thus, we must integrate with the measure
<italic>δ</italic>
(cone )
<italic>d </italic>
<bold>S</bold>
first. Next, we must take into account all the scattering sites on the line
<italic>DN</italic>
. Hence, we must perform a second integration with the measure
<italic>δ</italic>
(line)
<italic>d </italic>
<bold>N</bold>
. Consequently, the detected photon flux density at
<bold>D</bold>
in the direction
<bold>n</bold>
for a scattering angle
<italic>ω</italic>
is </p>
<p>
<disp-formula id="EEq5">
<label>(9)</label>
<mml:math id="M13">
<mml:mtable>
<mml:mtr>
<mml:mtd columnalign="right">
<mml:mover accent="true">
<mml:mrow>
<mml:mi>f</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold">D</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold">n</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>ω</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mtd>
<mml:mtd columnalign="left">
<mml:mo>=</mml:mo>
<mml:mstyle displaystyle="true">
<mml:mo>∫</mml:mo>
</mml:mstyle>
<mml:mstyle displaystyle="true">
<mml:mo>∫</mml:mo>
</mml:mstyle>
<mml:mi>δ</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>cone</mml:mi>
<mml:mo> </mml:mo>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold">S</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi mathvariant="bold">S</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>4</mml:mn>
<mml:mi>π</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mi>  </mml:mi>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mi>e</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>−</mml:mo>
<mml:mi>𝒜</mml:mi>
<mml:mo minsize="0.75em" maxsize="0.75em">(</mml:mo>
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:mi>N</mml:mi>
</mml:mrow>
<mml:mo minsize="0.75em" maxsize="0.75em">)</mml:mo>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:msup>
<mml:mrow>
<mml:mi>N</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd columnalign="right"></mml:mtd>
<mml:mtd columnalign="left">
<mml:mo> </mml:mo>
<mml:mo>×</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold">N</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi mathvariant="bold">N</mml:mi>
<mml:mi>π</mml:mi>
<mml:msubsup>
<mml:mrow>
<mml:mi>r</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>e</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>ω</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>δ</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mtext>line</mml:mtext>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mi>e</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>−</mml:mo>
<mml:mi>𝒜</mml:mi>
<mml:mo minsize="0.75em" maxsize="0.75em">(</mml:mo>
<mml:mrow>
<mml:mi>N</mml:mi>
<mml:mi>D</mml:mi>
</mml:mrow>
<mml:mo minsize="0.75em" maxsize="0.75em">)</mml:mo>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:mi>N</mml:mi>
<mml:msup>
<mml:mrow>
<mml:mi>D</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mfrac>
<mml:mo>,</mml:mo>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</disp-formula>
where
<italic>δ</italic>
(cone ) means the delta function concentrated on the cone of vertex
<bold>N</bold>
, opening angle
<italic>ω</italic>
, and axis
<inline-formula>
<mml:math id="M14">
<mml:mi mathvariant="bold">n</mml:mi>
<mml:mo>=</mml:mo>
<mml:mover accent="true">
<mml:mrow>
<mml:mi>D</mml:mi>
<mml:mi>N</mml:mi>
</mml:mrow>
<mml:mo>⃗</mml:mo>
</mml:mover>
<mml:mo>/</mml:mo>
<mml:mi>D</mml:mi>
<mml:mi>N</mml:mi>
</mml:math>
</inline-formula>
;
<italic>δ</italic>
(line) is the delta function concentrated on the line
<italic>DN</italic>
.</p>
<p>
<inline-formula>
<mml:math id="M15">
<mml:mover accent="true">
<mml:mrow>
<mml:mi>f</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi mathvariant="bold">D</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold">n</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>ω</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:math>
</inline-formula>
will be called compounded conical projection data in the direction
<bold>n</bold>
, at site
<bold>D</bold>
, and with given scattering angle
<italic>ω</italic>
[
<xref ref-type="bibr" rid="B30">20</xref>
].
<inline-formula>
<mml:math id="M16">
<mml:mover accent="true">
<mml:mrow>
<mml:mi>f</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi mathvariant="bold">D</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold">n</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>ω</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:math>
</inline-formula>
is also called the compounded conical Radon transform (CCRT) of
<italic>f</italic>
(
<bold>r</bold>
), a generalized Radon transform on circular cone surfaces as originally introduced in [
<xref ref-type="bibr" rid="B31">32</xref>
].</p>
<p> Assuming that exchange of integration is valid, we may view
<inline-formula>
<mml:math id="M17">
<mml:mover accent="true">
<mml:mrow>
<mml:mi>f</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi mathvariant="bold">D</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold">n</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>ω</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:math>
</inline-formula>
as </p>
<p>
<disp-formula id="EEq6">
<label>(10)</label>
<mml:math id="M18">
<mml:mtable>
<mml:mtr>
<mml:mtd columnalign="right">
<mml:mover accent="true">
<mml:mrow>
<mml:mi>f</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold">D</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold">n</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>ω</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mtd>
<mml:mtd columnalign="left">
<mml:mo>=</mml:mo>
<mml:mi>π</mml:mi>
<mml:msubsup>
<mml:mrow>
<mml:mi>r</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>e</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mi>𝒫</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>ω</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mstyle displaystyle="true">
<mml:mo>∫</mml:mo>
</mml:mstyle>
<mml:mi>d</mml:mi>
<mml:mi mathvariant="bold">N</mml:mi>
<mml:mi>δ</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mtext>line</mml:mtext>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold">N</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mi>e</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>−</mml:mo>
<mml:mi>𝒜</mml:mi>
<mml:mo minsize="0.75em" maxsize="0.75em">(</mml:mo>
<mml:mrow>
<mml:mi>N</mml:mi>
<mml:mi>D</mml:mi>
</mml:mrow>
<mml:mo minsize="0.75em" maxsize="0.75em">)</mml:mo>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:mi>N</mml:mi>
<mml:msup>
<mml:mrow>
<mml:mi>D</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd columnalign="right"></mml:mtd>
<mml:mtd columnalign="left">
<mml:mo> </mml:mo>
<mml:mo>×</mml:mo>
<mml:mstyle displaystyle="true">
<mml:mo>∫</mml:mo>
</mml:mstyle>
<mml:mfrac>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold">S</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi mathvariant="bold">S</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>4</mml:mn>
<mml:mi>π</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mi>δ</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>cone</mml:mi>
<mml:mo> </mml:mo>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mi>e</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>−</mml:mo>
<mml:mi>𝒜</mml:mi>
<mml:mo minsize="0.75em" maxsize="0.75em">(</mml:mo>
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:mi>N</mml:mi>
</mml:mrow>
<mml:mo minsize="0.75em" maxsize="0.75em">)</mml:mo>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:msup>
<mml:mrow>
<mml:mi>N</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mfrac>
<mml:mo>.</mml:mo>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</disp-formula>
This result may be regarded as the X-ray transform (on the straight line passing through
<bold>D</bold>
in the direction
<bold>n</bold>
) of the function </p>
<p>
<disp-formula id="EEq7">
<label>(11)</label>
<mml:math id="M19">
<mml:mi>g</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold">N</mml:mi>
<mml:mo> </mml:mo>
<mml:mo>∣</mml:mo>
<mml:mo> </mml:mo>
<mml:mi mathvariant="bold">n</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>ω</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold">N</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mi>e</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>−</mml:mo>
<mml:mi>𝒜</mml:mi>
<mml:mo minsize="0.75em" maxsize="0.75em">(</mml:mo>
<mml:mrow>
<mml:mi>N</mml:mi>
<mml:mi>D</mml:mi>
</mml:mrow>
<mml:mo minsize="0.75em" maxsize="0.75em">)</mml:mo>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:mi>N</mml:mi>
<mml:msup>
<mml:mrow>
<mml:mi>D</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mfrac>
<mml:mstyle displaystyle="true">
<mml:mo>∫</mml:mo>
</mml:mstyle>
<mml:mfrac>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold">S</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi mathvariant="bold">S</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>4</mml:mn>
<mml:mi>π</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mi>δ</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>cone</mml:mi>
<mml:mo> </mml:mo>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mi>e</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>−</mml:mo>
<mml:mi>𝒜</mml:mi>
<mml:mo minsize="0.75em" maxsize="0.75em">(</mml:mo>
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:mi>N</mml:mi>
</mml:mrow>
<mml:mo minsize="0.75em" maxsize="0.75em">)</mml:mo>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:msup>
<mml:mrow>
<mml:mi>N</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mfrac>
<mml:mo>,</mml:mo>
</mml:math>
</disp-formula>
or </p>
<p>
<disp-formula id="EEq8">
<label>(12)</label>
<mml:math id="M20">
<mml:mover accent="true">
<mml:mrow>
<mml:mi>f</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold">D</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold">n</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>ω</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mi>π</mml:mi>
<mml:msubsup>
<mml:mrow>
<mml:mi>r</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>e</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mi>𝒫</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>ω</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>X</mml:mi>
<mml:mi>g</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="bold">N</mml:mi>
<mml:mo> </mml:mo>
<mml:mo>∣</mml:mo>
<mml:mo> </mml:mo>
<mml:mi mathvariant="bold">n</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>ω</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:math>
</disp-formula>
</p>
<p> So by letting the directional detector (DD) take all possible spatial orientations, we generate the totality of possible compounded conical projections, which depends on five parameters: four for a line in ℝ
<sup>3</sup>
and one for the scattering angle (
<italic>ω</italic>
). In principle, the object under study is the support for two functions
<italic>f</italic>
(
<bold>r</bold>
), the object radioactivity distribution and
<italic>n</italic>
<sub>
<italic>e</italic>
</sub>
(
<bold>r</bold>
) its electron density distribution to be simultaneously determined. This problem is akin to the identification problem in the exponential Radon transform [
<xref ref-type="bibr" rid="B32">33</xref>
]. It seems to be difficult to solve at present since, by inspection, the DD data is underdetermined although it would be in principle possible to perform the inversion of the X-ray transform in (
<xref ref-type="disp-formula" rid="EEq8">12</xref>
), to retrieve
<italic>g</italic>
(
<bold>N</bold>
|
<bold>n</bold>
,
<italic>ω</italic>
). Thus, in the coming sections, we will review two major achievements of the compounded conical Radon transforms in ℝ
<sup>
<italic>n</italic>
</sup>
with
<italic>n</italic>
= 2,3 and discuss the properties of the corresponding imaging modalities. </p>
</sec>
</sec>
<sec sec-type="section" id="sec3">
<title>3. Compton-Scattered Radiation Imaging in Three Dimensions</title>
<sec sec-type="subsection" id="sec3.1">
<title>3.1. The Compounded Conical Radon Transform (CCRT)</title>
<p>In this section, we examine a tractable case by which a particular set of compounded conical projections is used for image reconstruction. This is possible under the following conditions: </p>
<list list-type="roman-lower">
<list-item>
<p>First-order scattering events are accounted for since they are vastly dominant and higher-order scattering are neglected [
<xref ref-type="bibr" rid="B6">6</xref>
,
<xref ref-type="bibr" rid="B33">34</xref>
], </p>
</list-item>
<list-item>
<p>the electron density is assumed to be constant. This is a reasonable hypothesis since most human tissues (brain cells, blood, muscles, lung tissues, water, etc.) have an electron density around 3.4 × 10
<sup>23</sup>
 cm
<sup>−3</sup>
. Their density is also around 1.0 g·cm
<sup>−3</sup>
[
<xref ref-type="bibr" rid="B35">19</xref>
]. This means for our purpose, objects containing bones should not be considered, </p>
</list-item>
<list-item>
<p>the set of compounded conical projections has one fixed direction
<bold>n</bold>
, parallel to the
<italic>Oz</italic>
axis direction, </p>
</list-item>
<list-item>
<p>the set of detecting pixels are distributed as array on a two-dimensional area, forming a collimated SPECT gamma camera. </p>
</list-item>
</list>
<p>Using the coordinate system of
<xref ref-type="fig" rid="fig5">Figure 5</xref>
and following the path of a photon from emission to absorption via one scattering at a site
<bold>N</bold>
, the expression of the flux density on the detector at a site
<bold>D</bold>
= (
<italic>x</italic>
<sub>
<italic>D</italic>
</sub>
,
<italic>y</italic>
<sub>
<italic>D</italic>
</sub>
, 0),
<inline-formula>
<mml:math id="M21">
<mml:mover accent="true">
<mml:mrow>
<mml:mi>f</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
<mml:mo stretchy="false">(</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>D</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>y</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>D</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>ω</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:math>
</inline-formula>
has the form of a linear integral transform of the object activity density
<italic>f</italic>
(
<italic>x</italic>
,
<italic>y</italic>
,
<italic>z</italic>
),</p>
<p>
<disp-formula id="EEq9">
<label>(13)</label>
<mml:math id="M22">
<mml:mtable>
<mml:mtr>
<mml:mtd columnalign="right"></mml:mtd>
<mml:mtd columnalign="left">
<mml:mover accent="true">
<mml:mrow>
<mml:mi>f</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>D</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>y</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>D</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>ω</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd columnalign="right"></mml:mtd>
<mml:mtd columnalign="left">
<mml:mo> </mml:mo>
<mml:mo> </mml:mo>
<mml:mo>=</mml:mo>
<mml:msub>
<mml:mstyle displaystyle="true">
<mml:mo>∫</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mi>ℝ</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>3</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:msub>
<mml:mi>d</mml:mi>
<mml:mi>x</mml:mi>
<mml:mo> </mml:mo>
<mml:mi>d</mml:mi>
<mml:mi>y</mml:mi>
<mml:mo> </mml:mo>
<mml:mi>d</mml:mi>
<mml:mi>z</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi>𝒦</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mtext>PSF</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>D</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>y</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>D</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>ω</mml:mi>
<mml:mo> </mml:mo>
<mml:mo>∣</mml:mo>
<mml:mo> </mml:mo>
<mml:mi>x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>y</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>z</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mover>
<mml:mrow>
<mml:mi>f</mml:mi>
</mml:mrow>
<mml:mo>¯</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>y</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>z</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</disp-formula>
where the function
<inline-formula>
<mml:math id="M23">
<mml:mover>
<mml:mrow>
<mml:mi>f</mml:mi>
</mml:mrow>
<mml:mo>¯</mml:mo>
</mml:mover>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>y</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>z</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:math>
</inline-formula>
is defined by </p>
<p>
<disp-formula id="EEq10">
<label>(14)</label>
<mml:math id="M24">
<mml:mover>
<mml:mrow>
<mml:mi>f</mml:mi>
</mml:mrow>
<mml:mo>¯</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>y</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>z</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>∫</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>∞</mml:mi>
</mml:mrow>
</mml:munderover>
<mml:mi>d</mml:mi>
<mml:mi>ζ</mml:mi>
<mml:mi>ν</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>ζ</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>y</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>z</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>ζ</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:math>
</disp-formula>
and the integration kernel is </p>
<p>
<disp-formula id="EEq11">
<label>(15)</label>
<mml:math id="M25">
<mml:mtable>
<mml:mtr>
<mml:mtd columnalign="right"></mml:mtd>
<mml:mtd columnalign="left">
<mml:msub>
<mml:mrow>
<mml:mi>𝒦</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mtext>PSF</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>D</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>y</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>D</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>ω</mml:mi>
<mml:mo> </mml:mo>
<mml:mo>∣</mml:mo>
<mml:mo> </mml:mo>
<mml:mi>x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>y</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>z</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd columnalign="right"></mml:mtd>
<mml:mtd columnalign="left">
<mml:mo> </mml:mo>
<mml:mo>=</mml:mo>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>ω</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>ν</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msqrt>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>−</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>D</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>y</mml:mi>
<mml:mo>−</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>y</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>D</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>z</mml:mi>
<mml:mo>−</mml:mo>
<mml:mi>ζ</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:msqrt>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd columnalign="right"></mml:mtd>
<mml:mtd columnalign="left">
<mml:mo> </mml:mo>
<mml:mo> </mml:mo>
<mml:mo>×</mml:mo>
<mml:mi>δ</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mi>cos</mml:mi>
<mml:mo> </mml:mo>
</mml:mrow>
<mml:mo> </mml:mo>
<mml:mrow>
<mml:mi>ω</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:msqrt>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>−</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>D</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>y</mml:mi>
<mml:mo>−</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>y</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>D</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:msqrt>
<mml:mo>−</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>z</mml:mi>
<mml:mo>−</mml:mo>
<mml:mi>ζ</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mi>sin</mml:mi>
</mml:mrow>
<mml:mo> </mml:mo>
<mml:mrow>
<mml:mi>ω</mml:mi>
</mml:mrow>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</disp-formula>
<italic>K</italic>
(
<italic>ω</italic>
) is the Compton kinematic factor, and
<italic>ν</italic>
(
<italic>d</italic>
) is a function describing a photometric factor for a distance
<italic>d</italic>
, for example,
<italic>ν</italic>
(
<italic>d</italic>
) = 1/
<italic>d</italic>
<sup>2</sup>
. By definition,
<inline-formula>
<mml:math id="M26">
<mml:mover accent="true">
<mml:mrow>
<mml:mi>f</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
<mml:mo stretchy="false">(</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>D</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>y</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>D</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>ω</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:math>
</inline-formula>
is called the compounded conical Radon transform (CCRT) of
<italic>f</italic>
(
<italic>x</italic>
,
<italic>y</italic>
,
<italic>z</italic>
).</p>
<p>The adopted working hypotheses are aimed to avoid unnecessary complications which would mask the main idea.</p>
<p>The inversion of the kernel
<italic>𝒦</italic>
<sub>PSF</sub>
(
<italic>x</italic>
<sub>
<italic>D</italic>
</sub>
,
<italic>y</italic>
<sub>
<italic>D</italic>
</sub>
,
<italic>ω</italic>
|
<italic>x</italic>
,
<italic>y</italic>
,
<italic>z</italic>
) is then obtained via a form of central slice theorem in Fourier space of the detector plane for the function
<inline-formula>
<mml:math id="M27">
<mml:mrow>
<mml:mover>
<mml:mrow>
<mml:mi>f</mml:mi>
</mml:mrow>
<mml:mo>¯</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</inline-formula>
, followed by a deconvolution to get the Fourier transform
<inline-formula>
<mml:math id="M28">
<mml:mover accent="true">
<mml:mrow>
<mml:mi>f</mml:mi>
</mml:mrow>
<mml:mo>˜</mml:mo>
</mml:mover>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>u</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>v</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>w</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:math>
</inline-formula>
of
<italic>f</italic>
(
<italic>x</italic>
,
<italic>y</italic>
,
<italic>z</italic>
), see [
<xref ref-type="bibr" rid="B31">32</xref>
,
<xref ref-type="bibr" rid="B34">35</xref>
],</p>
<p>
<disp-formula id="EEq12">
<label>(16)</label>
<mml:math id="M29">
<mml:mtable>
<mml:mtr>
<mml:mtd columnalign="right"></mml:mtd>
<mml:mtd columnalign="left">
<mml:mover accent="true">
<mml:mrow>
<mml:mi>f</mml:mi>
</mml:mrow>
<mml:mo>˜</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>u</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>v</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>w</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd columnalign="right"></mml:mtd>
<mml:mtd columnalign="left">
<mml:mo>=</mml:mo>
<mml:msub>
<mml:mstyle displaystyle="true">
<mml:mo>∫</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>ℝ</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>d</mml:mi>
<mml:mi>σ</mml:mi>
<mml:mi>exp</mml:mi>
<mml:mo> </mml:mo>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>i</mml:mi>
<mml:mi>π</mml:mi>
<mml:mi>σ</mml:mi>
<mml:mi>w</mml:mi>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:mo>−</mml:mo>
<mml:mrow>
<mml:mo>|</mml:mo>
<mml:mrow>
<mml:mi>z</mml:mi>
</mml:mrow>
<mml:mo>|</mml:mo>
</mml:mrow>
<mml:msqrt>
<mml:msup>
<mml:mrow>
<mml:mi>u</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mi>v</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:msqrt>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi>J</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>w</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:mi>  </mml:mi>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd columnalign="right"></mml:mtd>
<mml:mtd columnalign="left">
<mml:mo> </mml:mo>
<mml:mo>×</mml:mo>
<mml:msub>
<mml:mstyle displaystyle="true">
<mml:mo>∫</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>ℝ</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>+</mml:mo>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mi>t</mml:mi>
<mml:mo> </mml:mo>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi>J</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>π</mml:mi>
<mml:mrow>
<mml:mo>|</mml:mo>
<mml:mrow>
<mml:mi>z</mml:mi>
</mml:mrow>
<mml:mo>|</mml:mo>
</mml:mrow>
<mml:mi>t</mml:mi>
<mml:msqrt>
<mml:msup>
<mml:mrow>
<mml:mi>u</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mi>v</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:msqrt>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd columnalign="right"></mml:mtd>
<mml:mtd columnalign="left">
<mml:mo> </mml:mo>
<mml:mo>×</mml:mo>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>ω</mml:mi>
<mml:mo>−</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>π</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mo>∂</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mo>∂</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mfrac>
<mml:mrow>
<mml:mi>G</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>u</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>v</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:mo>+</mml:mo>
<mml:mi>H</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>π</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:mfrac>
<mml:mo>−</mml:mo>
<mml:mi>ω</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mo>∂</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mo>∂</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mfrac>
<mml:mrow>
<mml:mi>G</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>u</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>v</mml:mi>
<mml:mo>,</mml:mo>
<mml:mo>−</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mo>−</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</disp-formula>
where
<italic>J</italic>
(
<italic>w</italic>
) is the Fourier transform of
<italic>ν</italic>
(
<italic>x</italic>
),
<italic>J</italic>
<sub>1</sub>
(
<italic>x</italic>
) is the Bessel function of order 1,
<italic>H</italic>
(
<italic>x</italic>
) is the Heaviside unit step function,
<italic>t</italic>
= tan
<italic>ω</italic>
,
<italic>G</italic>
(
<italic>u</italic>
,
<italic>v</italic>
,
<italic>t</italic>
) is the Fourier transform of
<italic>g</italic>
(
<italic>x</italic>
<sub>
<italic>D</italic>
</sub>
,
<italic>y</italic>
<sub>
<italic>D</italic>
</sub>
,
<italic>t</italic>
),
<italic>K</italic>
(
<italic>t</italic>
) is the Compton kinematic factor as a function of
<italic>t</italic>
.</p>
<p>Finally,
<italic>f</italic>
(
<italic>x</italic>
,
<italic>y</italic>
,
<italic>z</italic>
) is recovered by three-dimensional inverse Fourier transform. It is observed that the data acquisition can be performed with a
<italic>spatially fixed</italic>
SPECT camera operating at successive scattered energies. As the scattering angle parameterizes series of “images” of the object, one may view it as replacing the spatial rotation angle in a standard SPECT data acquisition. This is a major advantage offered by this approach. </p>
<p>For practical purposes, the effective treatment of (
<xref ref-type="disp-formula" rid="EEq12">16</xref>
) is in itself a numerical challenge. It has been fully done in chapter 4 and appendices
<italic>A</italic>
and
<italic>B</italic>
of [
<xref ref-type="bibr" rid="B35">19</xref>
], where details on samplings in object space and medium space are given. Because the kernel of (
<xref ref-type="disp-formula" rid="EEq12">16</xref>
) is a Bessel function of order zero, to control the oscillations, an exponential discretization step is to be used as suggested in [
<xref ref-type="bibr" rid="B36">36</xref>
]. Along the line perpendicular to the planar detector, the high-frequency components of the activity function information are carried by the weakly scattered (or back-scattered) radiation, the sampling step should then be very small. But for the low-frequency components of the activity function, it is the strongly scattered photons (
<italic>ω</italic>
~
<italic>π</italic>
/2) which carry information; the sampling step should be then much larger. For a primary radiation energy of 140 keV, to obtain the same spatial resolution, the detector should have a very fine energy resolution Δ
<italic>E</italic>
. In fact, simulations show that the activity function reconstruction is satisfactory for Δ
<italic>E</italic>
~ 0.5 keV (even in the presence of a 24 dB white noise). The key point is that a reconstruction of reasonable quality can be achieved using the inverse CCRT. </p>
</sec>
<sec sec-type="subsection" id="sec3.2">
<title>3.2. Point Spread Function of the CCRT and Simulation Results</title>
<p>A way to get an idea of what the CCRT could be or do is to construct its point spread function (PSF), or the response function to a unit point source. It does not have the form of a delta-function as in the usual Radon transform. Because of the integration over all cones standing on top of the detecting site
<bold>D</bold>
, it appears as a function with the shape of a Mexican hat as shown in
<xref ref-type="fig" rid="fig6">Figure 6</xref>
; the point source is located somewhere on the vertical line symmetry axis of the Mexican hat above the planar detector [
<xref ref-type="bibr" rid="B34">35</xref>
].</p>
<p>The gamma detector operates now at a
<italic>fixed position</italic>
. No coincidence detection, as in Compton cameras, is required. Performed numerical simulations are in favor of the feasibility of this new imaging principle [
<xref ref-type="bibr" rid="B34">35</xref>
]. However issues related to higher-order scattering contribution, nonuniform attenuation, Poisson emission noise, detection sensitivity, and collimator efficiency are to be resolved before interesting practical modalities with possible combination with transmission imaging can be proposed [
<xref ref-type="bibr" rid="B37">37</xref>
].</p>
<p>To provide more convincing arguments regarding the viability of this idea, we present numerical simulations which illustrate the reconstruction of a simple cylindrical object using the analytic inversion formula with the following working conditions:</p>
<list list-type="roman-lower">
<list-item>
<p>the used
<italic>γ</italic>
-detector is a conventional SPECT camera. It has discretized dimensions
<italic>N</italic>
length units ×
<italic>N</italic>
length units. The length unit is arbitrary but should remain coherent with reality and in fact is taken equal to 1 mm. We have chosen
<italic>N</italic>
= 16 to keep the calculations required at a reasonable level, </p>
</list-item>
<list-item>
<p>the scattering medium is represented by a cube of dimensions
<italic>N</italic>
×
<italic>N</italic>
×
<italic>N</italic>
 (length unit)
<sup>3</sup>
,</p>
</list-item>
<list-item>
<p>the electron density in biological medium is
<italic>n</italic>
<sub>
<italic>e</italic>
</sub>
= 3.5 × 10
<sup>23</sup>
 electrons/cm
<sup>3</sup>
, </p>
</list-item>
<list-item>
<p>the radionuclide employed is Tc-99 with an activity concentration corresponding to 4.84 × 10
<sup>10</sup>
counts per minutes per cm
<sup>3</sup>
, </p>
</list-item>
<list-item>
<p>the acquisition time per projection is set to 0.1 sec, </p>
</list-item>
<list-item>
<p>the 3D original object (a cylinder of height 6 arbitrary units) is placed at the center of the scattering medium (cube), </p>
</list-item>
<list-item>
<p>the distance from the camera to the upper face of the scattering medium cube is
<italic>l</italic>
= 200 arbitrary length units. </p>
</list-item>
</list>
<p>
<xref ref-type="fig" rid="fig7">Figure 7</xref>
represents the original object.
<xref ref-type="fig" rid="fig8"> Figure 8</xref>
shows the series of images of the object at various scattering angles
<italic>ω</italic>
(5° <
<italic>ω</italic>
< 175°). In
<xref ref-type="fig" rid="fig9">Figure 9</xref>
, the reconstructed object in the absence of noise is illustrated with a relative mean square error (RMSE) = 1.2%, which is perfectly reasonable. We observe a good performance of the CCRT for modeling the new imaging process.</p>
<p>Concerning spatial resolution, the intrinsic resolution depends on the camera design (collimator, crystal, photomultiplier tubes, and measurement electronics). The reconstructed system resolution is further determined by the reconstruction algorithm used. The inclusion of scattered radiation increases considerably the number of detected photons, which might contribute to improve the signal-to-noise ratio (SNR) and the resolution of the imaging system. To evaluate accurately the spatial resolution, it is necessary to use real data and to compare it with conventional methods which do not make use of scattered radiation. At the present time, it is too early to use our preliminary simulation results for this purpose.</p>
<p>Since our main objective in this paper is to show how to exploit advantageously Compton-scattered radiation to propose a new imaging principle, we focus on results illustrating the image formation process as well as image reconstruction from scattered photons.</p>
<p>In real situations, of course, one must take into account other factors such as photon attenuation by the medium, Poisson emission noise, and imperfections of the detector system including the collimator and electronics.</p>
<p>The case of uniform attenuation (often assumed in the literature) was included in [
<xref ref-type="bibr" rid="B34">35</xref>
]. The exact treatment of inhomogeneous attenuation poses enormous mathematical difficulties. Concerning emission noise, several approaches have been suggested to deal with it such as maximum likelihood or wavelets method. They may be used for “denoising” the measured data beforehand or jointly with the inversion process.</p>
<p>As for the imperfections of the detector, the standard way for treating this problem is to make use of a response function usually modeled as a Gaussian defined both in spatial and energy coordinates. These issues are discussed in details in [
<xref ref-type="bibr" rid="B35">19</xref>
,
<xref ref-type="bibr" rid="B27">29</xref>
]. </p>
</sec>
<sec sec-type="subsection" id="sec3.3">
<title>3.3. A Possible Generalization and Its Numerical Test</title>
<p>As mentioned earlier, the presence of a mechanical collimator restricts severely the sensitivity of the imaging process. We have recently advocated a new functional modality following the principle of emission imaging by scattered gamma-rays without mechanical collimation [
<xref ref-type="bibr" rid="B37">37</xref>
,
<xref ref-type="bibr" rid="B38">38</xref>
]. Removing the collimator from the detector allows more gamma-rays to reach a detecting pixel from all directions coming from the upper half-space of this site, therefore increasing the strength of the signal (
<xref ref-type="fig" rid="fig10">Figure 10</xref>
). An introductory study in two dimensions has recently been performed [
<xref ref-type="bibr" rid="B39">39</xref>
] to show convincingly the viability of this idea and to motivate the present work.</p>
<p>The modeling of the image formation process is done by a more general compounded conical Radon transform, whereby one sums over conical projections at one detection pixel over all cone vertices in the upper half-space.
<xref ref-type="fig" rid="fig11"> Figure 11</xref>
shows the position of one conical projection in this generalized gCCRT. </p>
<p>This transform is obtained by summing over all scattering sites for a given site
<bold>D</bold>
on the detector. The mathematical expression of one arbitrary conical projection is quite involved and given in [
<xref ref-type="bibr" rid="B40">40</xref>
]. The summation over all such objects can be still expressed as a linear integral transform of the activity density
<italic>f</italic>
(
<italic>x</italic>
,
<italic>y</italic>
,
<italic>z</italic>
), </p>
<p>
<disp-formula id="EEq13">
<label>(17)</label>
<mml:math id="M30">
<mml:mi>g</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold">D</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>ω</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>∫</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mi>ℝ</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>3</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mrow></mml:mrow>
</mml:munderover>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>x</mml:mi>
<mml:mo> </mml:mo>
<mml:mi>d</mml:mi>
<mml:mi>y</mml:mi>
<mml:mo> </mml:mo>
<mml:mi>d</mml:mi>
<mml:mi>z</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi>𝒦</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mtext>PSF</mml:mtext>
</mml:mrow>
<mml:mrow>
<mml:mi>∗</mml:mi>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="bold">D</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>ω</mml:mi>
<mml:mo> </mml:mo>
<mml:mo>∣</mml:mo>
<mml:mo> </mml:mo>
<mml:mi>x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>y</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>z</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>y</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>z</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:math>
</disp-formula>
</p>
<p>Unfortunately, the explicit form of
<italic>𝒦</italic>
<sub>PSF*</sub>
(
<bold>D</bold>
,
<italic>ω</italic>
|
<italic>x</italic>
,
<italic>y</italic>
,
<italic>z</italic>
) is too complicated to yield a simple interpretation and will not be addressed here. However, this PSF, although no longer a computable function, is in fact an integral of the electronic density over the surface of a torus of revolution whose axis is the line connecting the point source to the detection point (
<xref ref-type="fig" rid="fig12">Figure 12</xref>
). The shape of the PSF is now completely different compared to the collimated geometry (
<xref ref-type="fig" rid="fig13">Figure 13</xref>
).</p>
<p>To compare with the collimated detector geometry,
<xref ref-type="fig" rid="fig14">Figure 14</xref>
shows the computed PSFs for both cases at a scattering angle of 30 degrees. The PSF without collimator is about 12 times “stronger” than that with collimator (approximatively 3200 counts without collimator compared to approximatively 250 counts with collimator; see [
<xref ref-type="bibr" rid="B38">38</xref>
]).</p>
<p>To demonstrate the viability of this idea, we have used data generated for a simple object and applied classical algebraic reconstruction methods. We have taken a simple source immersed in a cubic scattering medium. The source itself consists of two concentric cubes with different activity concentrations (
<xref ref-type="fig" rid="fig15">Figure 15</xref>
).</p>
<p>A 256 × 256 pixels detector is placed on the
<italic>xy</italic>
plane. The pixel size is 0.4 × 0.4 mm
<sup>2</sup>
. The scattering medium is a rectangular box of dimensions 30 cm by 30 cm by 15 cm, which is at a distance of 1 cm above the planar detector. The electronic density inside the scattering medium is
<italic>n</italic>
<sub>
<italic>e</italic>
</sub>
= 3.341023 electrons/cm
<sup>3</sup>
since most biological tissues have an electronic structure close to that of water. The radionuclide used in this simulation is
<sup>99
<italic>m </italic>
</sup>
Tc, which emits photons at an energy of 140.1 keV. The scattering medium is discretized with 13 voxels in
<italic>x</italic>
and
<italic>y</italic>
axis directions and with 9 voxels in
<italic>z</italic>
-axis direction. The detector is reduced to 13 × 13 pixels. We construct the weight matrix of the medium by calculating from our previous models, for each point of the mesh, the PSF of the detector at the different scattering angles. The reconstruction is carried out using the conjugated gradient method with positivity constraint; see
<xref ref-type="fig" rid="fig16">Figure 16</xref>
.</p>
<p>These results are an incentive to pursue the development of this imaging modality. </p>
</sec>
</sec>
<sec sec-type="section" id="sec4">
<title>4. Compton-Scattered Radiation Imaging in Two Dimensions</title>
<sec sec-type="subsection" id="sec4.1">
<title>4.1. Image formation-Compounded V-Line Radon Transform (CVLRT)</title>
<p>In this section, we discuss the transposition of the previous concept into a two-dimensional world. This passage entails a new form of the compounded conical projection, which is called now the compounded V-line projection, since the two-dimensional version of the cone surface is now just a geometric figure made up of two half-lines forming a letter V. This concept is illustrated in
<xref ref-type="fig" rid="fig17">Figure 17</xref>
.</p>
<p> Inspection shows that each compounded V-Line projection depends on three parameters: two parameters for the line, and
<italic>ω</italic>
the scattering angle. Thus, an attempt to determine simultaneously the electron density
<italic>n</italic>
<sub>
<italic>e</italic>
</sub>
(
<italic>x</italic>
,
<italic>y</italic>
) and the radioactivity density of the object with the totality of the compounded V-line projections in the plane will not be successful for lack of sufficient data. So we will examine only the case of constant
<italic>n</italic>
<sub>
<italic>e</italic>
</sub>
(
<italic>x</italic>
,
<italic>y</italic>
) within the hypotheses adopted in the previous section.</p>
<p>This imaging process concerns two-dimensional structures in biomedical imaging, in which a radiotracer has been injected and maintained on a support.
<xref ref-type="fig" rid="fig18"> Figure 18</xref>
shows the corresponding setup with a linear SPECT gamma camera for data acquisition.</p>
<p>In the image formation process, the compounded conical Radon transform is now replaced by the compounded V-line Radon transform. We give now the expression of the photon flux density at detecting site
<bold>D</bold>
for a scattering angle
<italic>ω</italic>
as in the previous section. Physical densities used here are actually derived from their three-dimensional values since we are dealing with real three-dimensional phenomena which are now restricted to a plane. They will appear with a star in formulas, for example,
<italic>n</italic>
<sub>
<italic>e</italic>
</sub>
*(
<italic>x</italic>
,
<italic>y</italic>
) in lieu of
<italic>n</italic>
<sub>
<italic>e</italic>
</sub>
(
<italic>x</italic>
,
<italic>y</italic>
,
<italic>z</italic>
). Let
<inline-formula>
<mml:math id="M31">
<mml:mover accent="true">
<mml:mrow>
<mml:mi>f</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi mathvariant="bold">D</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>τ</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
<mml:mo>=</mml:mo>
<mml:mover accent="true">
<mml:mrow>
<mml:mi>f</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>ζ</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>τ</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:math>
</inline-formula>
be the registered photon flux density, it is given by </p>
<p>
<disp-formula id="EEq14">
<label>(18)</label>
<mml:math id="M32">
<mml:mi>  </mml:mi>
<mml:mi>  </mml:mi>
<mml:mtable>
<mml:mtr>
<mml:mtd columnalign="right">
<mml:mover accent="true">
<mml:mrow>
<mml:mi>f</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>ζ</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>τ</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mtd>
<mml:mtd columnalign="left">
<mml:mo>=</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mi>K</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>∗</mml:mi>
</mml:mrow>
</mml:msup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>τ</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>∫</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>∞</mml:mi>
</mml:mrow>
</mml:munderover>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>η</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>η</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>∫</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>∞</mml:mi>
</mml:mrow>
</mml:munderover>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd columnalign="right"></mml:mtd>
<mml:mtd columnalign="left">
<mml:mo> </mml:mo>
<mml:mo>×</mml:mo>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>ζ</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mi>sin</mml:mi>
</mml:mrow>
<mml:mo> </mml:mo>
<mml:mrow>
<mml:mi>ω</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:mi>η</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mi>cos</mml:mi>
<mml:mo> </mml:mo>
</mml:mrow>
<mml:mo> </mml:mo>
<mml:mrow>
<mml:mi>ω</mml:mi>
</mml:mrow>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd columnalign="right"></mml:mtd>
<mml:mtd columnalign="left">
<mml:mo> </mml:mo>
<mml:mo> </mml:mo>
<mml:mo> </mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mo>+</mml:mo>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>ζ</mml:mi>
<mml:mo>−</mml:mo>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mi>sin</mml:mi>
</mml:mrow>
<mml:mo> </mml:mo>
<mml:mrow>
<mml:mi>ω</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:mi>η</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mi>cos</mml:mi>
<mml:mo> </mml:mo>
</mml:mrow>
<mml:mo> </mml:mo>
<mml:mrow>
<mml:mi>ω</mml:mi>
</mml:mrow>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</disp-formula>
where
<italic>τ</italic>
= tan 
<italic>ω</italic>
,
<italic>K</italic>
*(
<italic>ω</italic>
) =
<italic>πr</italic>
<sub>
<italic>e</italic>
</sub>
<sup>2</sup>
<italic>n</italic>
<sub>
<italic>e</italic>
</sub>
*
<italic>𝒫</italic>
(
<italic>ω</italic>
), the Compton kinematic factor restricted to two dimensions, 1/
<italic>η</italic>
and 1/
<italic>r</italic>
are the photometric propagation factor in two dimensions as opposed to the usual inverse square of the distance rule given in three dimensions. We call
<inline-formula>
<mml:math id="M33">
<mml:mover accent="true">
<mml:mrow>
<mml:mi>f</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi mathvariant="bold">D</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>ω</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:math>
</inline-formula>
the compounded V-line Radon transform (CVLRT) of
<italic>f</italic>
(
<italic>x</italic>
,
<italic>y</italic>
). Inversion of this transform would allow the reconstruction of
<italic>f</italic>
(
<italic>x</italic>
,
<italic>y</italic>
) under the assumptions cited above. </p>
</sec>
<sec sec-type="subsection" id="sec4.2">
<title>4.2. Inversion of the Compounded V-Line Radon Transform</title>
<p>The inversion procedure can be performed in two steps; see [
<xref ref-type="bibr" rid="B41">41</xref>
]. We present it with some details since it appears for the first time.</p>
<p>(i) First, let us define </p>
<p>
<disp-formula id="EEq15">
<label>(19)</label>
<mml:math id="M34">
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>∫</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>∞</mml:mi>
</mml:mrow>
</mml:munderover>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>η</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>η</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mi>  </mml:mi>
<mml:mi>  </mml:mi>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>ξ</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>±</mml:mo>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>η</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mi>cos</mml:mi>
<mml:mo> </mml:mo>
</mml:mrow>
<mml:mo> </mml:mo>
<mml:mrow>
<mml:mi>ω</mml:mi>
</mml:mrow>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mi>h</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>ξ</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>±</mml:mo>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mi>cos</mml:mi>
<mml:mo> </mml:mo>
</mml:mrow>
<mml:mo> </mml:mo>
<mml:mrow>
<mml:mi>ω</mml:mi>
</mml:mrow>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:math>
</disp-formula>
where
<italic>ξ</italic>
<sub>±</sub>
= (
<italic>ζ</italic>
±
<italic>r</italic>
sin
<italic>ω</italic>
). Then, this problem is now shifted to the problem of a simple V-line Radon transform (VLRT) with the integration measure
<italic>dr</italic>
/
<italic>r</italic>
on the function
<italic>h</italic>
(
<italic>ξ</italic>
,
<italic>r</italic>
cos  
<italic>ω</italic>
). Hence, the inversion problem of the CVLRT is that of the inversion of VLRT for
<italic>h</italic>
, followed by the extraction of
<italic>f</italic>
from (
<xref ref-type="disp-formula" rid="EEq15">19</xref>
),</p>
<p>
<disp-formula id="EEq16">
<label>(20)</label>
<mml:math id="M35">
<mml:mtable>
<mml:mtr>
<mml:mtd columnalign="right"></mml:mtd>
<mml:mtd columnalign="left">
<mml:mover accent="true">
<mml:mrow>
<mml:mi>f</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>ζ</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>ω</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd columnalign="right"></mml:mtd>
<mml:mtd columnalign="left">
<mml:mo>=</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mi>K</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>∗</mml:mi>
</mml:mrow>
</mml:msup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>ω</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>∫</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>∞</mml:mi>
</mml:mrow>
</mml:munderover>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd columnalign="right"></mml:mtd>
<mml:mtd columnalign="left">
<mml:mo> </mml:mo>
<mml:mo>×</mml:mo>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:mi>h</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>ζ</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mi>sin</mml:mi>
</mml:mrow>
<mml:mo> </mml:mo>
<mml:mrow>
<mml:mi>ω</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:mi>η</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mi>cos</mml:mi>
<mml:mo> </mml:mo>
</mml:mrow>
<mml:mo> </mml:mo>
<mml:mrow>
<mml:mi>ω</mml:mi>
</mml:mrow>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>+</mml:mo>
<mml:mi>h</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>ζ</mml:mi>
<mml:mo>−</mml:mo>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mi>sin</mml:mi>
</mml:mrow>
<mml:mo> </mml:mo>
<mml:mrow>
<mml:mi>ω</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:mi>η</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mi>cos</mml:mi>
<mml:mo> </mml:mo>
</mml:mrow>
<mml:mo> </mml:mo>
<mml:mrow>
<mml:mi>ω</mml:mi>
</mml:mrow>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</disp-formula>
</p>
<p>The inversion formula for the VLRT can be worked out following [
<xref ref-type="bibr" rid="B42">42</xref>
]. It yields </p>
<p>
<disp-formula id="EEq17">
<label>(21)</label>
<mml:math id="M36">
<mml:mtable>
<mml:mtr>
<mml:mtd columnalign="right"></mml:mtd>
<mml:mtd columnalign="left">
<mml:mi>h</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>y</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd columnalign="right"></mml:mtd>
<mml:mtd columnalign="left">
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>y</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>π</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>∫</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>π</mml:mi>
</mml:mrow>
</mml:munderover>
<mml:mi>d</mml:mi>
<mml:mi>τ</mml:mi>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd columnalign="right"></mml:mtd>
<mml:mtd columnalign="left">
<mml:mo> </mml:mo>
<mml:mo>×</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>  </mml:mi>
<mml:mtext>P</mml:mtext>
<mml:mo>.</mml:mo>
<mml:mtext>V</mml:mtext>
<mml:mo>.</mml:mo>
<mml:mi>  </mml:mi>
<mml:msub>
<mml:mstyle displaystyle="true">
<mml:mo>∫</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>ℝ</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>d</mml:mi>
<mml:mi>ζ</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>ζ</mml:mi>
<mml:mo>−</mml:mo>
<mml:mi>x</mml:mi>
<mml:mo>−</mml:mo>
<mml:mi>y</mml:mi>
<mml:mi>τ</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>+</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>ζ</mml:mi>
<mml:mo>−</mml:mo>
<mml:mi>x</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>y</mml:mi>
<mml:mi>τ</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mo>∂</mml:mo>
<mml:mover accent="true">
<mml:mrow>
<mml:mi>f</mml:mi>
</mml:mrow>
<mml:mo>^</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>ζ</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>τ</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mo>∂</mml:mo>
<mml:mi>ζ</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</disp-formula>
</p>
<p>  (ii) Now, knowing
<italic>h</italic>
(
<italic>x</italic>
,
<italic>y</italic>
),
<italic>f</italic>
is the solution of the convolution equation </p>
<p>
<disp-formula id="EEq18">
<label>(22)</label>
<mml:math id="M37">
<mml:mi>h</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>y</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>∫</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>∞</mml:mi>
</mml:mrow>
</mml:munderover>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>η</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>η</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>η</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>y</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:math>
</disp-formula>
Thus,
<italic>f</italic>
(
<italic>x</italic>
,
<italic>y</italic>
) can be extracted by Fourier techniques [
<xref ref-type="bibr" rid="B43">43</xref>
]. Let </p>
<p>
<disp-formula id="EEq19">
<label>(23)</label>
<mml:math id="M38">
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>y</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:msub>
<mml:mstyle displaystyle="true">
<mml:mo>∫</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>ℝ</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>d</mml:mi>
<mml:mi>k</mml:mi>
<mml:msup>
<mml:mrow>
<mml:mi>e</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>i</mml:mi>
<mml:mi>π</mml:mi>
<mml:mi>k</mml:mi>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:msup>
<mml:mi>  </mml:mi>
<mml:mover>
<mml:mrow>
<mml:mi>f</mml:mi>
</mml:mrow>
<mml:mo>¯</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>∫</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>∞</mml:mi>
</mml:mrow>
</mml:munderover>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>η</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>η</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:msup>
<mml:mrow>
<mml:mi>e</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>−</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>i</mml:mi>
<mml:mi>π</mml:mi>
<mml:mi>k</mml:mi>
<mml:mi>η</mml:mi>
</mml:mrow>
</mml:msup>
<mml:mo>=</mml:mo>
<mml:mo>−</mml:mo>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:mi>log</mml:mi>
<mml:mo> </mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>π</mml:mi>
<mml:mrow>
<mml:mo>|</mml:mo>
<mml:mrow>
<mml:mi>k</mml:mi>
</mml:mrow>
<mml:mo>|</mml:mo>
</mml:mrow>
<mml:mo>+</mml:mo>
<mml:mi>γ</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>i</mml:mi>
<mml:mfrac>
<mml:mrow>
<mml:mi>π</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:mfrac>
<mml:mi>sgn</mml:mi>
<mml:mo> </mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</disp-formula>
be the relevant Fourier transforms and
<italic>γ</italic>
= 0.57721566... the Euler's constant, then (
<xref ref-type="disp-formula" rid="EEq18">22</xref>
) becomes </p>
<p>
<disp-formula id="EEq20">
<label>(24)</label>
<mml:math id="M39">
<mml:mi>h</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>y</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:msub>
<mml:mstyle displaystyle="true">
<mml:mo>∫</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>ℝ</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>d</mml:mi>
<mml:mi>k</mml:mi>
<mml:msup>
<mml:mrow>
<mml:mi>e</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>i</mml:mi>
<mml:mi>π</mml:mi>
<mml:mi>k</mml:mi>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:msup>
<mml:mi>  </mml:mi>
<mml:mover>
<mml:mrow>
<mml:mi>f</mml:mi>
</mml:mrow>
<mml:mo>¯</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mo>−</mml:mo>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:mi>log</mml:mi>
<mml:mo> </mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>π</mml:mi>
<mml:mrow>
<mml:mo>|</mml:mo>
<mml:mrow>
<mml:mi>k</mml:mi>
</mml:mrow>
<mml:mo>|</mml:mo>
</mml:mrow>
<mml:mo>+</mml:mo>
<mml:mi>γ</mml:mi>
<mml:mo>−</mml:mo>
<mml:mi>i</mml:mi>
<mml:mfrac>
<mml:mrow>
<mml:mi>π</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:mfrac>
<mml:mi>sgn</mml:mi>
<mml:mo> </mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:math>
</disp-formula>
The recovery of
<italic>f</italic>
(
<italic>x</italic>
,
<italic>y</italic>
) is achieved by inverse Fourier transform </p>
<p>
<disp-formula id="EEq21">
<label>(25)</label>
<mml:math id="M40">
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>y</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:msub>
<mml:mstyle displaystyle="true">
<mml:mo>∫</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>ℝ</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>d</mml:mi>
<mml:mi>q</mml:mi>
<mml:msup>
<mml:mrow>
<mml:mi>e</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>i</mml:mi>
<mml:mi>π</mml:mi>
<mml:mi>q</mml:mi>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:msup>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mo>∫</mml:mo>
<mml:mrow>
<mml:mi>ℝ</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>d</mml:mi>
<mml:mi>z</mml:mi>
<mml:msup>
<mml:mrow>
<mml:mi>e</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>−</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>i</mml:mi>
<mml:mi>π</mml:mi>
<mml:mi>z</mml:mi>
<mml:mi>q</mml:mi>
</mml:mrow>
</mml:msup>
<mml:mi>h</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>z</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mo>−</mml:mo>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:mi>log</mml:mi>
<mml:mo> </mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>π</mml:mi>
<mml:mrow>
<mml:mo>|</mml:mo>
<mml:mrow>
<mml:mi>q</mml:mi>
</mml:mrow>
<mml:mo>|</mml:mo>
</mml:mrow>
<mml:mo>+</mml:mo>
<mml:mi>γ</mml:mi>
<mml:mo>−</mml:mo>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>π</mml:mi>
<mml:mo>/</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>sgn</mml:mi>
<mml:mo> </mml:mo>
<mml:mi>q</mml:mi>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:mi>  </mml:mi>
<mml:mo>,</mml:mo>
</mml:math>
</disp-formula>
as
<italic>h</italic>
(
<italic>x</italic>
,
<italic>z</italic>
) is known for all
<italic>z</italic>
∈ ℝ.</p>
<p>Formula (
<xref ref-type="disp-formula" rid="EEq17">21</xref>
) lends itself to the derivation of a filtered back-projection (FBP) method for image reconstruction. This is because its structure is basically the one found in the standard Radon transform. Essentially, the procedure consists of the sum of two filtered back-projections on two half-lines forming the letter V. This is at first not obvious at all since in general the inverse of the sum of two operators is not necessarily the sum of their inverses. The advantage of the filtered back-projection inversion formula is that it can be implemented by fast algorithms. Simulation results are presented next. </p>
</sec>
<sec sec-type="subsection" id="sec4.3">
<title>4.3. Numerical Simulations</title>
<p>We present now the results of numerical simulations. The original image (
<xref ref-type="fig" rid="fig19">Figure 19</xref>
) of size 512 × 512 of length units is a thyroid phantom presenting with small nodules.
<xref ref-type="fig" rid="fig20"> Figure 20</xref>
shows the CVLT transform of a thyroid phantom with angular sampling rate
<italic>dω</italic>
= 0.005 rad and 314 projections (
<italic>π</italic>
/2/0.005 = 314) which are the images of Compton-scattered radiation on the camera in terms of the distance
<italic>ξ</italic>
and the scattering angle
<italic>ω</italic>
. The reconstruction using FBP is given in
<xref ref-type="fig" rid="fig21">Figure 21</xref>
. The artifacts are due to the limited number of projections. Moreover, back-projection on V-lines generates more artifacts than back-projection on straight lines, because of more spurious line intersections. As our numerical results are based on the discretization of the inverse formula, a choice of a smaller discretization step
<italic>dω</italic>
would improve image quality. This is indeed a well-established fact and in agreement with the improved sampling resulting from the increase of data collected at more values of the scattering angle
<italic>ω</italic>
. Despite these limitations, the small structures in the object are clearly reconstructed. This result illustrates the feasibility of the new imaging modality, for which the main advantage resides in the use of a
<italic> one-dimensional nonmoving</italic>
Compton camera for
<italic> two-dimensional</italic>
imaging [
<xref ref-type="bibr" rid="B42">42</xref>
]. To show the possible medical application of this imaging modality, we have also presented the simulations on a Shepp-Logan phantom as a second example; see Figures
<xref ref-type="fig" rid="fig22">22</xref>
,
<xref ref-type="fig" rid="fig23">23</xref>
, and
<xref ref-type="fig" rid="fig24">24</xref>
. </p>
</sec>
</sec>
<sec sec-type="section" id="sec5">
<title>5. Conclusion</title>
<p>The aim of this paper on Compton-scattered radiation emission imaging is to review the progress and development steps in this field during the last decade. The theoretical work is based on the inversion of generalized (and compounded) Radon transforms defined on cone surfaces in three dimensions and on V-lines in two dimensions. Numerical simulations appear to support the viability of the corresponding imaging process. The natural next research step to undertake is to study how this imaging principle behave under realistic operating conditions, namely</p>
<list list-type="roman-lower">
<list-item>
<p>how does it react to energy uncertainty measurements, </p>
</list-item>
<list-item>
<p>how the CCRT inversion is modified if the electron recoil is not negligible, </p>
</list-item>
<list-item>
<p>how to include the imperfections of the detection system (collimator response and electronic system response), </p>
</list-item>
<list-item>
<p>how to account for limited statistics and random noise in typical real data acquisition situations. </p>
</list-item>
</list>
<p>In spite of these open questions, we still believe that this imaging technique has a strong potential as a future variant of the SPECT modality as soon as high performance detectors with very fine energy resolution are made available and accessible on a large scale. The future of scattered radiation emission imaging lies perhaps in its foremost advantage; it can provide a three-dimensional image reconstruction without having to move or displace the gamma camera (or planar detector). As discussed earlier, the usual rotation angle of an SPECT camera is now replaced with the Compton scattering angle. Thus, in scattered radiation emission imaging, the gamma camera is motionless and as such it does not require a heavy mechanical apparatus to accurately rotate the gamma camera around an object. It will be therefore less cumbersome and more convenient to use under stringent space conditions. In the end, it has led to a new concept of high energy resolution photon detector, a new development in data acquisition, and new image reconstruction methods derived from the compounded conical Radon transform and compounded V-line Radon transform. Applications in two- or three-dimensional imaging are possible particularly in medicine for clinical diagnosis and treatment planning. Finally, the tantalizing subject of identification via the concept of compounded conical projection looks very exciting and motivates future research in this direction. </p>
</sec>
</body>
<back>
<ack>
<title>Acknowledgment</title>
<p>This work was supported by the Swiss National Science Foundation under Grant no. 3152A0-102143. </p>
</ack>
<ref-list>
<ref id="B1">
<label>1</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tuy</surname>
<given-names>H</given-names>
</name>
</person-group>
<article-title>Inversion formula for cone-beam reconstruction</article-title>
<source>
<italic>SIAM Journal on Applied Mathematics</italic>
</source>
<year>1983</year>
<volume>43</volume>
<issue>3</issue>
<fpage>546</fpage>
<lpage>552</lpage>
</element-citation>
</ref>
<ref id="B2">
<label>2</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhao</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Yu</surname>
<given-names>H</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>G</given-names>
</name>
</person-group>
<article-title>A unified framework for exact cone-beam reconstruction formulas</article-title>
<source>
<italic>Medical Physics</italic>
</source>
<year>2005</year>
<volume>32</volume>
<issue>6</issue>
<fpage>1712</fpage>
<lpage>1721</lpage>
<pub-id pub-id-type="pmid">16013729</pub-id>
</element-citation>
</ref>
<ref id="B3">
<label>3</label>
<element-citation publication-type="book">
<person-group person-group-type="author">
<name>
<surname>Wernick</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Aarsvold</surname>
<given-names>J</given-names>
</name>
</person-group>
<source>
<italic>Emission Tomography: The Fundamentals of PET and SPECT</italic>
</source>
<year>2004</year>
<publisher-loc>San Diego, Calif, USA</publisher-loc>
<publisher-name>Academic Press</publisher-name>
</element-citation>
</ref>
<ref id="B4">
<label>4</label>
<element-citation publication-type="book">
<person-group person-group-type="author">
<name>
<surname>Cherry</surname>
<given-names>SR</given-names>
</name>
<name>
<surname>Sorenson</surname>
<given-names>JA</given-names>
</name>
<name>
<surname>Phelps</surname>
<given-names>ME</given-names>
</name>
</person-group>
<source>
<italic>Physics in Nuclear Medicine</italic>
</source>
<year>2004</year>
<publisher-loc>Philadelphia, Pa, USA</publisher-loc>
<publisher-name>Elsevier Health Sciences</publisher-name>
</element-citation>
</ref>
<ref id="B5">
<label>5</label>
<element-citation publication-type="book">
<person-group person-group-type="author">
<name>
<surname>Grangeat</surname>
<given-names>P</given-names>
</name>
</person-group>
<source>
<italic>La Tomographie</italic>
</source>
<year>2006</year>
<publisher-loc>Paris, France</publisher-loc>
<publisher-name>Hermès Science</publisher-name>
</element-citation>
</ref>
<ref id="B6">
<label>6</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zaidi</surname>
<given-names>H</given-names>
</name>
<name>
<surname>Koral</surname>
<given-names>KF</given-names>
</name>
</person-group>
<article-title>Scatter modelling and compensation in emission tomography</article-title>
<source>
<italic>European Journal of Nuclear Medicine and Molecular Imaging</italic>
</source>
<year>2004</year>
<volume>31</volume>
<issue>5</issue>
<fpage>761</fpage>
<lpage>782</lpage>
<pub-id pub-id-type="pmid">15057488</pub-id>
</element-citation>
</ref>
<ref id="B7">
<label>7</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lale</surname>
<given-names>PG</given-names>
</name>
</person-group>
<article-title>The examination of internal tissues, using gamma-ray scatter with a possible extension to megavoltage radiography</article-title>
<source>
<italic>Physics in Medicine and Biology</italic>
</source>
<year>1959</year>
<volume>4</volume>
<fpage>159</fpage>
<lpage>167</lpage>
<pub-id pub-id-type="pmid">14413525</pub-id>
</element-citation>
</ref>
<ref id="B8">
<label>8</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lenti</surname>
<given-names>M</given-names>
</name>
</person-group>
<article-title>A 3-D imaging device using Compton scattering off the body</article-title>
<source>
<italic>Nuclear Instruments and Methods in Physics Research A</italic>
</source>
<year>2008</year>
<volume>588</volume>
<issue>3</issue>
<fpage>457</fpage>
<lpage>462</lpage>
</element-citation>
</ref>
<ref id="B9">
<label>9</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Guzzardi</surname>
<given-names>R</given-names>
</name>
<name>
<surname>Licitra</surname>
<given-names>G</given-names>
</name>
</person-group>
<article-title>Critical review of Compton imaging</article-title>
<source>
<italic>Critical Reviews in Biomedical Engineering</italic>
</source>
<year>1988</year>
<volume>15</volume>
<issue>3</issue>
<fpage>237</fpage>
<lpage>268</lpage>
<pub-id pub-id-type="pmid">3329596</pub-id>
</element-citation>
</ref>
<ref id="B10">
<label>10</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Harding</surname>
<given-names>G</given-names>
</name>
</person-group>
<article-title>Inelastic photon scattering: effects and applications in biomedical science and industry</article-title>
<source>
<italic>Radiation Physics and Chemistry</italic>
</source>
<year>1997</year>
<volume>50</volume>
<issue>1</issue>
<fpage>91</fpage>
<lpage>111</lpage>
</element-citation>
</ref>
<ref id="B11">
<label>11</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Norton</surname>
<given-names>SJ</given-names>
</name>
</person-group>
<article-title>Compton scattering tomography</article-title>
<source>
<italic>Journal of Applied Physics</italic>
</source>
<year>1994</year>
<volume>76</volume>
<issue>4</issue>
<fpage>2007</fpage>
<lpage>2015</lpage>
</element-citation>
</ref>
<ref id="B12">
<label>12</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hussein</surname>
<given-names>EMA</given-names>
</name>
</person-group>
<article-title>Radiation scattering methods for nondestructive testing and imaging</article-title>
<source>
<italic>International Advances in Nondestructive Testing</italic>
</source>
<year>1989</year>
<volume>14</volume>
<fpage>301</fpage>
<lpage>321</lpage>
</element-citation>
</ref>
<ref id="B13">
<label>13</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gautam</surname>
<given-names>SR</given-names>
</name>
<name>
<surname>Hopkins</surname>
<given-names>FF</given-names>
</name>
<name>
<surname>Klinksiek</surname>
<given-names>R</given-names>
</name>
<name>
<surname>Morgan</surname>
<given-names>IL</given-names>
</name>
</person-group>
<article-title>Compton interaction tomography I: feasibility studies for applications in earthquake engineering</article-title>
<source>
<italic>IEEE Transactions on Nuclear Science</italic>
</source>
<year>1982</year>
<volume>30</volume>
<issue>2</issue>
<fpage>1680</fpage>
<lpage>1684</lpage>
</element-citation>
</ref>
<ref id="B14">
<label>14</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Todd</surname>
<given-names>RW</given-names>
</name>
<name>
<surname>Nightingale</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Everett</surname>
<given-names>D</given-names>
</name>
</person-group>
<article-title>A proposed
<italic>γ</italic>
camera</article-title>
<source>
<italic>Nature</italic>
</source>
<year>1974</year>
<volume>251</volume>
<issue>5471</issue>
<fpage>132</fpage>
<lpage>134</lpage>
</element-citation>
</ref>
<ref id="B15">
<label>15</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Singh</surname>
<given-names>M</given-names>
</name>
</person-group>
<article-title>An electronically collimated gamma camera for single photon emission computed tomography, part I: theoretical considerations and design criteria</article-title>
<source>
<italic>Medical Physics</italic>
</source>
<year>1983</year>
<volume>10</volume>
<issue>4</issue>
<fpage>421</fpage>
<lpage>427</lpage>
<pub-id pub-id-type="pmid">6604216</pub-id>
</element-citation>
</ref>
<ref id="B16">
<label>16</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Singh</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Doria</surname>
<given-names>D</given-names>
</name>
</person-group>
<article-title>Single photon imaging with electronic collimation</article-title>
<source>
<italic>IEEE Transactions on Nuclear Science</italic>
</source>
<year>1984</year>
<volume>32</volume>
<issue>1</issue>
<fpage>843</fpage>
<lpage>847</lpage>
</element-citation>
</ref>
<ref id="B17">
<label>17</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gerl</surname>
<given-names>J</given-names>
</name>
</person-group>
<article-title>Gamma-ray imaging exploiting the Compton effect</article-title>
<source>
<italic>Nuclear Physics A</italic>
</source>
<year>2005</year>
<volume>752</volume>
<issue>1–4</issue>
<fpage>688</fpage>
<lpage>695</lpage>
</element-citation>
</ref>
<ref id="B18">
<label>18</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Braem</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Llatas</surname>
<given-names>MChamizo</given-names>
</name>
<name>
<surname>Chesi</surname>
<given-names>E</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Feasibility of a novel design of high resolution parallax-free compton enhanced PET scanner dedicated to brain research</article-title>
<source>
<italic>Physics in Medicine and Biology</italic>
</source>
<year>2004</year>
<volume>49</volume>
<issue>12</issue>
<fpage>2547</fpage>
<lpage>2562</lpage>
<pub-id pub-id-type="pmid">15272673</pub-id>
</element-citation>
</ref>
<ref id="B35">
<label>19</label>
<element-citation publication-type="book">
<person-group person-group-type="author">
<name>
<surname>Delarbre</surname>
<given-names>JL</given-names>
</name>
</person-group>
<source>
<italic>Imagerie d’mission gamma: nouvelles mthodes d’exploitation du rayonnement diffus</italic>
</source>
<year>2005</year>
<publisher-loc>Paris, France</publisher-loc>
<publisher-name>University of Cergy-Pontoise</publisher-name>
<comment>Ph.D. thesis</comment>
</element-citation>
</ref>
<ref id="B30">
<label>20</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nguyen</surname>
<given-names>MK</given-names>
</name>
<name>
<surname>Truong</surname>
<given-names>TT</given-names>
</name>
<name>
<surname>Driol</surname>
<given-names>C</given-names>
</name>
<name>
<surname>Zaidi</surname>
<given-names>H</given-names>
</name>
</person-group>
<article-title>On a novel approach to compton scattered emission imaging</article-title>
<source>
<italic>IEEE Transactions on Nuclear Science</italic>
</source>
<year>2009</year>
<volume>56</volume>
<issue>3</issue>
<fpage>1430</fpage>
<lpage>1437</lpage>
</element-citation>
</ref>
<ref id="B19">
<label>21</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zeng</surname>
<given-names>GL</given-names>
</name>
<name>
<surname>Gullberg</surname>
<given-names>GT</given-names>
</name>
<name>
<surname>Tsui</surname>
<given-names>BMW</given-names>
</name>
<name>
<surname>Terry</surname>
<given-names>JA</given-names>
</name>
</person-group>
<article-title>Three-dimensional iterative reconstruction algorithms with attenuation and geometric point response correction</article-title>
<source>
<italic>IEEE Transactions on Nuclear Science</italic>
</source>
<year>1991</year>
<volume>38</volume>
<issue>2</issue>
<fpage>693</fpage>
<lpage>702</lpage>
</element-citation>
</ref>
<ref id="B20">
<label>22</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zeng</surname>
<given-names>GL</given-names>
</name>
<name>
<surname>Clack</surname>
<given-names>R</given-names>
</name>
<name>
<surname>Gullberg</surname>
<given-names>GT</given-names>
</name>
</person-group>
<article-title>Implementation of Tuy’s cone-beam inversion formula</article-title>
<source>
<italic>Physics in Medicine and Biology</italic>
</source>
<year>1994</year>
<volume>39</volume>
<issue>3</issue>
<fpage>493</fpage>
<lpage>507</lpage>
<pub-id pub-id-type="pmid">15551594</pub-id>
</element-citation>
</ref>
<ref id="B21">
<label>23</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kadrmas</surname>
<given-names>DJ</given-names>
</name>
<name>
<surname>Frey</surname>
<given-names>EC</given-names>
</name>
<name>
<surname>Karimi</surname>
<given-names>SS</given-names>
</name>
<name>
<surname>Tsui</surname>
<given-names>BMW</given-names>
</name>
</person-group>
<article-title>Fast implementations of reconstruction-based scatter compensation in fully 3D SPECT image reconstruction</article-title>
<source>
<italic>Physics in Medicine and Biology</italic>
</source>
<year>1998</year>
<volume>43</volume>
<issue>4</issue>
<fpage>857</fpage>
<lpage>873</lpage>
<pub-id pub-id-type="pmid">9572510</pub-id>
</element-citation>
</ref>
<ref id="B22">
<label>24</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pack</surname>
<given-names>JD</given-names>
</name>
<name>
<surname>Noo</surname>
<given-names>F</given-names>
</name>
<name>
<surname>Clackdoyle</surname>
<given-names>R</given-names>
</name>
</person-group>
<article-title>Cone-beam reconstruction using the backprojection of locally filtered projections</article-title>
<source>
<italic>IEEE Transactions on Medical Imaging</italic>
</source>
<year>2005</year>
<volume>24</volume>
<issue>1</issue>
<fpage>70</fpage>
<lpage>85</lpage>
<pub-id pub-id-type="pmid">15638187</pub-id>
</element-citation>
</ref>
<ref id="B23">
<label>25</label>
<element-citation publication-type="book">
<person-group person-group-type="author">
<name>
<surname>Barrett</surname>
<given-names>HH</given-names>
</name>
<name>
<surname>Swindell</surname>
<given-names>W</given-names>
</name>
</person-group>
<source>
<italic>Swindell Radiological Imaging</italic>
</source>
<year>1981</year>
<publisher-loc>New York, NY, USA</publisher-loc>
<publisher-name>Academic Press</publisher-name>
</element-citation>
</ref>
<ref id="B24">
<label>26</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Floyd</surname>
<given-names>CE</given-names>
<suffix>Jr.</suffix>
</name>
<name>
<surname>Jaszczak</surname>
<given-names>RJ</given-names>
</name>
<name>
<surname>Greer</surname>
<given-names>KL</given-names>
</name>
<name>
<surname>Coleman</surname>
<given-names>RE</given-names>
</name>
</person-group>
<article-title>Deconvolution of Compton scatter in SPECT</article-title>
<source>
<italic>Journal of Nuclear Medicine</italic>
</source>
<year>1985</year>
<volume>26</volume>
<issue>4</issue>
<fpage>403</fpage>
<lpage>408</lpage>
<pub-id pub-id-type="pmid">3872353</pub-id>
</element-citation>
</ref>
<ref id="B25">
<label>27</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gunter</surname>
<given-names>DL</given-names>
</name>
</person-group>
<article-title>Three-dimensional imaging utilizing energy discrimination II</article-title>
<source>
<italic>Journal of the Optical Society of America A</italic>
</source>
<year>1992</year>
<volume>9</volume>
<issue>6</issue>
<fpage>903</fpage>
<lpage>923</lpage>
</element-citation>
</ref>
<ref id="B26">
<label>28</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kadrmas</surname>
<given-names>DJ</given-names>
</name>
<name>
<surname>Frey</surname>
<given-names>EC</given-names>
</name>
<name>
<surname>Tsui</surname>
<given-names>BMW</given-names>
</name>
</person-group>
<article-title>Analysis of the reconstructibility and noise properties of scattered photons in Tc SPECT</article-title>
<source>
<italic>Physics in Medicine and Biology</italic>
</source>
<year>1997</year>
<volume>42</volume>
<issue>12</issue>
<fpage>2493</fpage>
<lpage>2516</lpage>
<pub-id pub-id-type="pmid">9434303</pub-id>
</element-citation>
</ref>
<ref id="B27">
<label>29</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nguyen</surname>
<given-names>MK</given-names>
</name>
<name>
<surname>Faye</surname>
<given-names>C</given-names>
</name>
<name>
<surname>Eglin</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Truong</surname>
<given-names>TT</given-names>
</name>
</person-group>
<article-title>Apparent image formation by compton-scattered photons in gamma-ray imaging</article-title>
<source>
<italic>IEEE Signal Processing Letters</italic>
</source>
<year>2001</year>
<volume>8</volume>
<issue>9</issue>
<fpage>248</fpage>
<lpage>251</lpage>
</element-citation>
</ref>
<ref id="B28">
<label>30</label>
<element-citation publication-type="book">
<person-group person-group-type="author">
<name>
<surname>Eglin</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Nguyen</surname>
<given-names>MK</given-names>
</name>
<name>
<surname>Faye</surname>
<given-names>C</given-names>
</name>
</person-group>
<person-group person-group-type="editor">
<name>
<surname>Djafari</surname>
<given-names>AM</given-names>
</name>
</person-group>
<article-title>Restoration from a multi-energy scintigraphic image sequence in the Bayesian framework</article-title>
<source>
<italic>Bayesian Inference and Maximum Entropy Methods in Science and Engineering</italic>
</source>
<year>2001</year>
<publisher-name>American Institute of Physics</publisher-name>
<fpage>513</fpage>
<lpage>520</lpage>
</element-citation>
</ref>
<ref id="B29">
<label>31</label>
<element-citation publication-type="book">
<person-group person-group-type="author">
<name>
<surname>Evans</surname>
<given-names>RD</given-names>
</name>
</person-group>
<source>
<italic>The Atomic Nucleus</italic>
</source>
<year>1955</year>
<publisher-loc>New York, NY, USA</publisher-loc>
<publisher-name>McGraw-Hill</publisher-name>
</element-citation>
</ref>
<ref id="B31">
<label>32</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nguyen</surname>
<given-names>MK</given-names>
</name>
<name>
<surname>Truong</surname>
<given-names>TT</given-names>
</name>
</person-group>
<article-title>On an integral transform and its inverse in nuclear imaging</article-title>
<source>
<italic>Inverse Problems</italic>
</source>
<year>2002</year>
<volume>18</volume>
<issue>1</issue>
<fpage>265</fpage>
<lpage>277</lpage>
</element-citation>
</ref>
<ref id="B32">
<label>33</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hertle</surname>
<given-names>A</given-names>
</name>
</person-group>
<article-title>The identification problem for the constantly attenuated radon transform</article-title>
<source>
<italic>Mathematische Zeitschrift</italic>
</source>
<year>1988</year>
<volume>197</volume>
<issue>1</issue>
<fpage>13</fpage>
<lpage>19</lpage>
</element-citation>
</ref>
<ref id="B33">
<label>34</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zaidi</surname>
<given-names>H</given-names>
</name>
</person-group>
<article-title>Relevance of accurate Monte Carlo modeling in nuclear medical imaging</article-title>
<source>
<italic>Medical Physics </italic>
</source>
<year>1999</year>
<volume>26</volume>
<issue>4</issue>
<fpage>574</fpage>
<lpage>608</lpage>
<pub-id pub-id-type="pmid">10227362</pub-id>
</element-citation>
</ref>
<ref id="B34">
<label>35</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nguyen</surname>
<given-names>MK</given-names>
</name>
<name>
<surname>Truong</surname>
<given-names>TT</given-names>
</name>
<name>
<surname>Bui</surname>
<given-names>HD</given-names>
</name>
<name>
<surname>Delarbre</surname>
<given-names>JL</given-names>
</name>
</person-group>
<article-title>A novel inverse problem in
<italic>γ</italic>
-rays emission imaging</article-title>
<source>
<italic>Inverse Problems in Science and Engineering</italic>
</source>
<year>2004</year>
<volume>12</volume>
<issue>2</issue>
<fpage>225</fpage>
<lpage>246</lpage>
</element-citation>
</ref>
<ref id="B36">
<label>36</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Siegman</surname>
<given-names>AE</given-names>
</name>
</person-group>
<article-title>Quasi-fast Hankel transform</article-title>
<source>
<italic>Optics Letters</italic>
</source>
<year>1997</year>
<volume>1</volume>
<issue>1</issue>
<fpage>13</fpage>
<lpage>15</lpage>
<pub-id pub-id-type="pmid">19680315</pub-id>
</element-citation>
</ref>
<ref id="B37">
<label>37</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nguyen</surname>
<given-names>MK</given-names>
</name>
<name>
<surname>Truong</surname>
<given-names>TT</given-names>
</name>
<name>
<surname>Delarbre</surname>
<given-names>JL</given-names>
</name>
<name>
<surname>Roux</surname>
<given-names>C</given-names>
</name>
<name>
<surname>Zaidi</surname>
<given-names>H</given-names>
</name>
</person-group>
<article-title>Novel approach to stationary transmission scanning using Compton scattered radiation</article-title>
<source>
<italic>Physics in Medicine and Biology</italic>
</source>
<year>2007</year>
<volume>52</volume>
<issue>15</issue>
<fpage>4615</fpage>
<lpage>4632</lpage>
<pub-id pub-id-type="pmid">17634654</pub-id>
</element-citation>
</ref>
<ref id="B38">
<label>38</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Driol</surname>
<given-names>C</given-names>
</name>
<name>
<surname>Nguyen</surname>
<given-names>MK</given-names>
</name>
<name>
<surname>Truong</surname>
<given-names>TT</given-names>
</name>
</person-group>
<article-title>Modeling and simulation results on high sensitivity scattered gamma-ray emission imaging</article-title>
<source>
<italic>Simulation Modelling Practice and Theory</italic>
</source>
<year>2008</year>
<volume>16</volume>
<issue>8</issue>
<fpage>1067</fpage>
<lpage>1076</lpage>
</element-citation>
</ref>
<ref id="B39">
<label>39</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nguyen</surname>
<given-names>MK</given-names>
</name>
<name>
<surname>Driol</surname>
<given-names>C</given-names>
</name>
<name>
<surname>Truong</surname>
<given-names>TT</given-names>
</name>
<name>
<surname>Zaidi</surname>
<given-names>H</given-names>
</name>
</person-group>
<article-title>Towards a new concept for high sensitivity compton scatter emission imaging</article-title>
<source>
<italic>Journal of the European Optical Society</italic>
</source>
<year>2008</year>
<volume>3, article 08010</volume>
<comment>Article ID 08010.</comment>
</element-citation>
</ref>
<ref id="B40">
<label>40</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Truong</surname>
<given-names>TT</given-names>
</name>
<name>
<surname>Nguyen</surname>
<given-names>MK</given-names>
</name>
<name>
<surname>Zaidi</surname>
<given-names>H</given-names>
</name>
</person-group>
<article-title>The mathematical foundations of 3D compton scatter emission imaging</article-title>
<source>
<italic>International Journal of Biomedical Imaging</italic>
</source>
<year>2007</year>
<volume>2007</volume>
<comment>Article ID 92780.</comment>
</element-citation>
</ref>
<ref id="B41">
<label>41</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Truong</surname>
<given-names>TT</given-names>
</name>
<name>
<surname>Nguyen</surname>
<given-names>MK</given-names>
</name>
</person-group>
<article-title>On new V-line radon transforms in ℝ
<sup>2</sup>
and their inversion</article-title>
<source>
<italic>Journal of Physics A</italic>
</source>
<year>2011</year>
<volume>44</volume>
<comment>Article ID 075206.</comment>
</element-citation>
</ref>
<ref id="B42">
<label>42</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Morvidone</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Nguyen</surname>
<given-names>MK</given-names>
</name>
<name>
<surname>Truong</surname>
<given-names>TT</given-names>
</name>
<name>
<surname>Zaidi</surname>
<given-names>H</given-names>
</name>
</person-group>
<article-title>On the V-line radon transform and its imaging applications</article-title>
<source>
<italic>International Journal of Biomedical Imaging</italic>
</source>
<year>2010</year>
<volume>2010</volume>
<fpage>6 pages</fpage>
<comment>Article ID 208179.</comment>
</element-citation>
</ref>
<ref id="B43">
<label>43</label>
<element-citation publication-type="book">
<person-group person-group-type="author">
<name>
<surname>Lavoine</surname>
<given-names>J</given-names>
</name>
</person-group>
<source>
<italic>Transformation de Fourier des Pseudo-Fonctions avec Tables de Nouvelles Transformes</italic>
</source>
<year>1963</year>
<publisher-loc>Paris, France</publisher-loc>
<publisher-name>CNRS</publisher-name>
</element-citation>
</ref>
</ref-list>
</back>
<floats-group>
<fig id="fig1" position="float">
<label>Figure 1</label>
<caption>
<p>Linear projection through an object.</p>
</caption>
<graphic xlink:href="IJBI2011-913893.001"></graphic>
</fig>
<fig id="fig2" position="float">
<label>Figure 2</label>
<caption>
<p>Illustration of the problem of photon attenuation resulting from photoelectric absorption and Compton scattering.</p>
</caption>
<graphic xlink:href="IJBI2011-913893.002"></graphic>
</fig>
<fig id="fig3" position="float">
<label>Figure 3</label>
<caption>
<p>Bone scintigraphy: a standard scintigraphic image (left image) is restored with the use of Compton-scattered radiation. Hot spots or nodule are clearly displayed (right image—reprinted with permission from [
<xref ref-type="bibr" rid="B35">19</xref>
]).</p>
</caption>
<graphic xlink:href="IJBI2011-913893.003"></graphic>
</fig>
<fig id="fig4" position="float">
<label>Figure 4</label>
<caption>
<p>Compounded conical projection.</p>
</caption>
<graphic xlink:href="IJBI2011-913893.004"></graphic>
</fig>
<fig id="fig5" position="float">
<label>Figure 5</label>
<caption>
<p>Coordinate system used in CCRT computation.</p>
</caption>
<graphic xlink:href="IJBI2011-913893.005"></graphic>
</fig>
<fig id="fig6" position="float">
<label>Figure 6</label>
<caption>
<p>The PSF at a scattering angle of 50 degrees with collimator (reprinted with permission from [
<xref ref-type="bibr" rid="B30">20</xref>
]).</p>
</caption>
<graphic xlink:href="IJBI2011-913893.006"></graphic>
</fig>
<fig id="fig7" position="float">
<label>Figure 7</label>
<caption>
<p>Original object (cylinder) in a cube consisting of 16 transaxial planes (reprinted with permission from [
<xref ref-type="bibr" rid="B30">20</xref>
]).</p>
</caption>
<graphic xlink:href="IJBI2011-913893.007"></graphic>
</fig>
<fig id="fig8" position="float">
<label>Figure 8</label>
<caption>
<p>Series of images parameterized by the angle of Compton scattering
<italic>ω</italic>
(5° <
<italic>ω</italic>
< 175°) (reprinted with permission from [
<xref ref-type="bibr" rid="B30">20</xref>
]).</p>
</caption>
<graphic xlink:href="IJBI2011-913893.008"></graphic>
</fig>
<fig id="fig9" position="float">
<label>Figure 9</label>
<caption>
<p>Reconstructed object in the absence of noise (RMSE = 1.2%) (reprinted with permission from [
<xref ref-type="bibr" rid="B30">20</xref>
]).</p>
</caption>
<graphic xlink:href="IJBI2011-913893.009"></graphic>
</fig>
<fig id="fig10" position="float">
<label>Figure 10</label>
<caption>
<p>Two imaging modalities using Compton-scattered radiation with and without collimation (reprinted with permission from [
<xref ref-type="bibr" rid="B30">20</xref>
]).</p>
</caption>
<graphic xlink:href="IJBI2011-913893.010"></graphic>
</fig>
<fig id="fig11" position="float">
<label>Figure 11</label>
<caption>
<p>A conical projection in the gCCRT.</p>
</caption>
<graphic xlink:href="IJBI2011-913893.011"></graphic>
</fig>
<fig id="fig12" position="float">
<label>Figure 12</label>
<caption>
<p>Torus surface of scattering sites for the case without collimator (reprinted with permission from [
<xref ref-type="bibr" rid="B30">20</xref>
]).</p>
</caption>
<graphic xlink:href="IJBI2011-913893.012"></graphic>
</fig>
<fig id="fig13" position="float">
<label>Figure 13</label>
<caption>
<p>The PSF at a scattering angle of 50 degrees without collimator (reprinted with permission from [
<xref ref-type="bibr" rid="B30">20</xref>
]).</p>
</caption>
<graphic xlink:href="IJBI2011-913893.013"></graphic>
</fig>
<fig id="fig14" position="float">
<label>Figure 14</label>
<caption>
<p>Comparison of the PSF with (lower blue line) and without collimator (upper red line), (reprinted with permission from [
<xref ref-type="bibr" rid="B30">20</xref>
]).</p>
</caption>
<graphic xlink:href="IJBI2011-913893.014"></graphic>
</fig>
<fig id="fig15" position="float">
<label>Figure 15</label>
<caption>
<p>Two representative slices of the original object illustrating two transaxial slices corresponding to the 6th and 9th planes, respectively (reprinted with permission from [
<xref ref-type="bibr" rid="B30">20</xref>
]).</p>
</caption>
<graphic xlink:href="IJBI2011-913893.015"></graphic>
</fig>
<fig id="fig16" position="float">
<label>Figure 16</label>
<caption>
<p>Two representative slices of the reconstructed images corresponding to the object shown in
<xref ref-type="fig" rid="fig15">Figure 15</xref>
(reprinted with permission from [
<xref ref-type="bibr" rid="B30">20</xref>
]).</p>
</caption>
<graphic xlink:href="IJBI2011-913893.016"></graphic>
</fig>
<fig id="fig17" position="float">
<label>Figure 17</label>
<caption>
<p>Compounded V-Line projection.</p>
</caption>
<graphic xlink:href="IJBI2011-913893.017"></graphic>
</fig>
<fig id="fig18" position="float">
<label>Figure 18</label>
<caption>
<p>Setup and parameters of the compounded V-line Radon transform (CVLT).</p>
</caption>
<graphic xlink:href="IJBI2011-913893.018"></graphic>
</fig>
<fig id="fig19" position="float">
<label>Figure 19</label>
<caption>
<p>Original thyroid phantom.</p>
</caption>
<graphic xlink:href="IJBI2011-913893.019"></graphic>
</fig>
<fig id="fig20" position="float">
<label>Figure 20</label>
<caption>
<p>The
<italic>T𝕍</italic>
transform of the thyroid image shown in
<xref ref-type="fig" rid="fig19">Figure 19</xref>
with
<italic>dω</italic>
= 0.0025 rad.</p>
</caption>
<graphic xlink:href="IJBI2011-913893.020"></graphic>
</fig>
<fig id="fig21" position="float">
<label>Figure 21</label>
<caption>
<p>FBP-IM reconstruction of the thyroid phantom with
<italic>dω</italic>
= 0.0025 rad.</p>
</caption>
<graphic xlink:href="IJBI2011-913893.021"></graphic>
</fig>
<fig id="fig22" position="float">
<label>Figure 22</label>
<caption>
<p>Shepp-Logan phantom.</p>
</caption>
<graphic xlink:href="IJBI2011-913893.022"></graphic>
</fig>
<fig id="fig23" position="float">
<label>Figure 23</label>
<caption>
<p>The
<italic>T𝕍</italic>
transform of the Shepp-Logan image shown in
<xref ref-type="fig" rid="fig19">Figure 19</xref>
with
<italic>dω</italic>
= 0.0025 rad.</p>
</caption>
<graphic xlink:href="IJBI2011-913893.023"></graphic>
</fig>
<fig id="fig24" position="float">
<label>Figure 24</label>
<caption>
<p>FBP-IM reconstruction of the Shepp-Logan phantom with
<italic>dω</italic>
= 0.0025 rad.</p>
</caption>
<graphic xlink:href="IJBI2011-913893.024"></graphic>
</fig>
</floats-group>
</pmc>
<affiliations>
<list>
<country>
<li>Argentine</li>
<li>France</li>
</country>
<region>
<li>Île-de-France</li>
</region>
<settlement>
<li>Cergy-Pontoise</li>
</settlement>
</list>
<tree>
<noCountry>
<name sortKey="Zaidi, H" sort="Zaidi, H" uniqKey="Zaidi H" first="H." last="Zaidi">H. Zaidi</name>
</noCountry>
<country name="France">
<region name="Île-de-France">
<name sortKey="Nguyen, M K" sort="Nguyen, M K" uniqKey="Nguyen M" first="M. K." last="Nguyen">M. K. Nguyen</name>
</region>
<name sortKey="Truong, T T" sort="Truong, T T" uniqKey="Truong T" first="T. T." last="Truong">T. T. Truong</name>
</country>
<country name="Argentine">
<noRegion>
<name sortKey="Morvidone, M" sort="Morvidone, M" uniqKey="Morvidone M" first="M." last="Morvidone">M. Morvidone</name>
</noRegion>
</country>
</tree>
</affiliations>
</record>

Pour manipuler ce document sous Unix (Dilib)

EXPLOR_STEP=$WICRI_ROOT/Wicri/Musique/explor/MozartV1/Data/Pmc/Checkpoint

HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000317 | SxmlIndent | more

Ou

HfdSelect -h $EXPLOR_AREA/Data/Pmc/Checkpoint/biblio.hfd -nk 000317 | SxmlIndent | more

Pour mettre un lien sur cette page dans le réseau Wicri

{{Explor lien
   |wiki=    Wicri/Musique
   |area=    MozartV1
   |flux=    Pmc
   |étape=   Checkpoint
   |type=    RBID
   |clé=     PMC:3124682
   |texte=   Scattered Radiation Emission Imaging: Principles and Applications
}}

Pour générer des pages wiki

HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Checkpoint/RBID.i   -Sk "pubmed:21747823" \
       | HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Checkpoint/biblio.hfd   \
       | NlmPubMed2Wicri -a MozartV1
Wicri

This area was generated with Dilib version V0.6.20.
Data generation: Sun Apr 10 15:06:14 2016. Site generation: Tue Feb 7 15:40:35 2023