{"id":35,"date":"2018-02-07T12:26:39","date_gmt":"2018-02-07T17:26:39","guid":{"rendered":"https:\/\/pctk.jhu.edu\/?page_id=35"},"modified":"2018-02-09T14:35:56","modified_gmt":"2018-02-09T19:35:56","slug":"some-details-of-pctk","status":"publish","type":"page","link":"https:\/\/pctk.jhu.edu\/index.php\/some-details-of-pctk\/","title":{"rendered":"What PcTK can do"},"content":{"rendered":"<p><strong>How it works<\/strong><\/p>\n<p>An x-ray photon interacts with photon counting detectors (PCDs) via various phenomena and generates an electron charge cloud or multiple clouds.<\/p>\n<figure id=\"attachment_36\" aria-describedby=\"caption-attachment-36\" style=\"width: 300px\" class=\"wp-caption aligncenter\"><a href=\"https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/interactions.png\"><img loading=\"lazy\" decoding=\"async\" class=\"wp-image-36 size-medium\" src=\"https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/interactions-300x138.png\" alt=\"\" width=\"300\" height=\"138\" srcset=\"https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/interactions-300x138.png 300w, https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/interactions-768x353.png 768w, https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/interactions-1024x471.png 1024w, https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/interactions-587x270.png 587w, https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/interactions.png 1239w\" sizes=\"auto, (max-width: 300px) 100vw, 300px\" \/><\/a><figcaption id=\"caption-attachment-36\" class=\"wp-caption-text\">Fig. 1. Various interaction phenomena: (a) No interaction with PCD; (b-d) Photoelectric effects (b) with total absorption, (c) fluorescence x-ray emission and its escape from PCD, and (d) fluorescence x-ray emission and its re-absorption; and (e) Compton scattering and Rayleigh scattering. Compton scattering and Rayleigh scattering are not included in PcTK version 3.2. Figure is from Ref. 2.<\/figcaption><\/figure>\n<p>The clouds may be detected by two or more adjacent PCD pixels and the recorded energy may be lower than the original energy due to charge sharing (Fig. 1b), escape and reabsorption of a K-shell fluorescence x-ray (Figs. 1c\u20131d), etc. This is called double-counting or <em>n<\/em>-tuple-counting. As a result, PCD data are spatially and energetically correlated. The Photon Counting Toolkit (PcTK) models the x-ray spectrum of recorded PCD data, with such correlation.<\/p>\n<p><strong>PcTK output<\/strong><\/p>\n<p>Users define parameters for PCDs such as the pixel size, the pixel depth, the reference size of electron charges, the electronic noise, etc. PcTK then calculates probabilities of the recorded energy and count at neighboring 3&#215;3 PCD pixels when x-rays are incident onto the center of the pixels (pixel 5 in Fig. 2) and outputs a normalized covariance matrix for 1-keV-width energy windows or a given\u00a0energy windows. It will take minutes or hours, depending on parameters.<\/p>\n<figure id=\"attachment_34\" aria-describedby=\"caption-attachment-34\" style=\"width: 300px\" class=\"wp-caption aligncenter\"><a href=\"https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/Echarges.png\"><img loading=\"lazy\" decoding=\"async\" class=\"wp-image-34 size-medium\" src=\"https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/Echarges-300x149.png\" alt=\"\" width=\"300\" height=\"149\" srcset=\"https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/Echarges-300x149.png 300w, https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/Echarges-768x382.png 768w, https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/Echarges-1024x509.png 1024w, https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/Echarges-543x270.png 543w, https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/Echarges.png 1130w\" sizes=\"auto, (max-width: 300px) 100vw, 300px\" \/><\/a><figcaption id=\"caption-attachment-34\" class=\"wp-caption-text\">Fig. 2. An example of electron charge distribution when (a) no fluorescence x-ray is emitted and (b) a fluorescence x-ray is emitted and absorbed by adjacent pixels. Figure is from Ref. 1.<\/figcaption><\/figure>\n<figure id=\"attachment_30\" aria-describedby=\"caption-attachment-30\" style=\"width: 300px\" class=\"wp-caption aligncenter\"><a href=\"https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/nCovE_nCovW.png\"><img loading=\"lazy\" decoding=\"async\" class=\"wp-image-30 size-medium\" src=\"https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/nCovE_nCovW-300x167.png\" alt=\"\" width=\"300\" height=\"167\" srcset=\"https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/nCovE_nCovW-300x167.png 300w, https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/nCovE_nCovW-768x427.png 768w, https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/nCovE_nCovW-485x270.png 485w, https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/nCovE_nCovW.png 992w\" sizes=\"auto, (max-width: 300px) 100vw, 300px\" \/><\/a><figcaption id=\"caption-attachment-30\" class=\"wp-caption-text\">Fig. 3. Normalized covariance matrices for\u00a01-keV-width energy windows (a) and <em>Nl<\/em>(=4) energy windows (b). Figure is from Ref. 1.<\/figcaption><\/figure>\n<p><strong>Use cases<\/strong><\/p>\n<p>From the output of PcTK, users can generate the recorded spectra as shown in Figs. 4\u20136. This process is very fast.<\/p>\n<figure id=\"attachment_31\" aria-describedby=\"caption-attachment-31\" style=\"width: 300px\" class=\"wp-caption aligncenter\"><a href=\"https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/SRE_CXDpixels.png\"><img loading=\"lazy\" decoding=\"async\" class=\"wp-image-31 size-medium\" src=\"https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/SRE_CXDpixels-300x234.png\" alt=\"\" width=\"300\" height=\"234\" srcset=\"https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/SRE_CXDpixels-300x234.png 300w, https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/SRE_CXDpixels-768x600.png 768w, https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/SRE_CXDpixels-346x270.png 346w, https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/SRE_CXDpixels.png 892w\" sizes=\"auto, (max-width: 300px) 100vw, 300px\" \/><\/a><figcaption id=\"caption-attachment-31\" class=\"wp-caption-text\">Fig. 4. The recorded spectra at each of the 3&#215;3 neighboring pixels when only the central pixel receives x-rays (from Ref. 1)<\/figcaption><\/figure>\n<figure id=\"attachment_32\" aria-describedby=\"caption-attachment-32\" style=\"width: 300px\" class=\"wp-caption aligncenter\"><a href=\"https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/SRE_E1s.png\"><img loading=\"lazy\" decoding=\"async\" class=\"wp-image-32 size-medium\" src=\"https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/SRE_E1s-300x239.png\" alt=\"\" width=\"300\" height=\"239\" srcset=\"https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/SRE_E1s-300x239.png 300w, https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/SRE_E1s-338x270.png 338w, https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/SRE_E1s.png 633w\" sizes=\"auto, (max-width: 300px) 100vw, 300px\" \/><\/a><figcaption id=\"caption-attachment-32\" class=\"wp-caption-text\">Fig. 5 The recorded spectra with different incident energies when all of the pixels receive x-rays (from Ref. 1)<\/figcaption><\/figure>\n<figure id=\"attachment_33\" aria-describedby=\"caption-attachment-33\" style=\"width: 161px\" class=\"wp-caption aligncenter\"><a href=\"https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/SRE_PCDs.png\"><img loading=\"lazy\" decoding=\"async\" class=\"wp-image-33 size-medium\" src=\"https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/SRE_PCDs-161x300.png\" alt=\"\" width=\"161\" height=\"300\" srcset=\"https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/SRE_PCDs-161x300.png 161w, https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/SRE_PCDs-549x1024.png 549w, https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/SRE_PCDs-145x270.png 145w, https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/SRE_PCDs.png 666w\" sizes=\"auto, (max-width: 161px) 100vw, 161px\" \/><\/a><figcaption id=\"caption-attachment-33\" class=\"wp-caption-text\">Fig. 6 The recorded spectra with different PCD parameters (from Ref. 1)<\/figcaption><\/figure>\n<p>Using the workflow script, users can generate correlated noisy PCD data (Figs. 7f\u20137i) from &#8220;material-specific sinograms&#8221; (C in Fig. 7)) and other parameters (A\u2013E in Fig. 7).<\/p>\n<figure id=\"attachment_29\" aria-describedby=\"caption-attachment-29\" style=\"width: 563px\" class=\"wp-caption aligncenter\"><a href=\"https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/CTworkflow.png\"><img loading=\"lazy\" decoding=\"async\" class=\"wp-image-29 size-large\" src=\"https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/CTworkflow-563x1024.png\" alt=\"\" width=\"563\" height=\"1024\" srcset=\"https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/CTworkflow-563x1024.png 563w, https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/CTworkflow-165x300.png 165w, https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/CTworkflow-149x270.png 149w, https:\/\/pctk.jhu.edu\/wp-content\/uploads\/2018\/02\/CTworkflow.png 696w\" sizes=\"auto, (max-width: 563px) 100vw, 563px\" \/><\/a><figcaption id=\"caption-attachment-29\" class=\"wp-caption-text\">Fig. 7. PCD-CT workflow using PcTK. Figure is from Ref. 1.<\/figcaption><\/figure>\n<p>There is a caveat. Generating correlated noisy data that are correlated among neighboring pixels, allows for more realistic PCD simulation, but takes longer time (hours). We provide a script that generates <em>uncorrelated<\/em> noisy PCD data with the same expectation. Generating uncorrelated noisy data, where all of PCD data are independent to each other, allows for more efficient computation (minutes) at the expense of the lack of noise correlation. The bottle-neck of the scripts is the way to call a random number generator of Matlab\u2014more specifically, a random number generator is called for each set of 3&#215;3 neighboring PCD pixels. It would be great if users can offer suggestions on more efficient ways or GPU-based programming.<\/p>\n<p>Ref. 1. Taguchi K, Stierstorfer K, Polster C, Lee O, and Kappler S. Spatio-energetic cross-talk in photon counting detectors: Numerical detector model (PcTK) and workflow for CT image quality assessment. Medical Physics. 2018;45(xx):xxx\u2013xxx. (tentatively accepted). This paper outlines PcTK version 3.2.<\/p>\n<p>Ref. 2. Taguchi K, Polster C, Lee O, Stierstorfer K, and Kappler S. Spatio-energetic cross-talk in photon counting detectors: Detector model and correlated Poisson data generator. Medical Physics. 2016;43(12):6386\u20136404. DOI:<a href=\"https:\/\/www.ncbi.nlm.nih.gov\/pubmed\/27908175\">10.1118\/1.4966699<\/a>\u00a0This paper outlines PcTK version 2.1.<\/p>\n<p>&nbsp;<\/p>\n","protected":false},"excerpt":{"rendered":"<p>How it works An x-ray photon interacts with photon counting detectors (PCDs) via various phenomena and generates an electron charge cloud or multiple clouds. The clouds may be detected by two or more adjacent PCD pixels and the recorded energy may be lower than the original energy due to charge sharing (Fig. 1b), escape and &hellip; <\/p>\n<p class=\"read-more\"><a class=\"btn btn-default\" href=\"https:\/\/pctk.jhu.edu\/index.php\/some-details-of-pctk\/\"> Read More<span class=\"screen-reader-text\">  Read More<\/span><\/a><\/p>\n","protected":false},"author":1,"featured_media":0,"parent":0,"menu_order":0,"comment_status":"closed","ping_status":"closed","template":"","meta":{"_exactmetrics_skip_tracking":false,"_exactmetrics_sitenote_active":false,"_exactmetrics_sitenote_note":"","_exactmetrics_sitenote_category":0,"footnotes":""},"class_list":["post-35","page","type-page","status-publish","hentry"],"_links":{"self":[{"href":"https:\/\/pctk.jhu.edu\/index.php\/wp-json\/wp\/v2\/pages\/35","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/pctk.jhu.edu\/index.php\/wp-json\/wp\/v2\/pages"}],"about":[{"href":"https:\/\/pctk.jhu.edu\/index.php\/wp-json\/wp\/v2\/types\/page"}],"author":[{"embeddable":true,"href":"https:\/\/pctk.jhu.edu\/index.php\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/pctk.jhu.edu\/index.php\/wp-json\/wp\/v2\/comments?post=35"}],"version-history":[{"count":7,"href":"https:\/\/pctk.jhu.edu\/index.php\/wp-json\/wp\/v2\/pages\/35\/revisions"}],"predecessor-version":[{"id":67,"href":"https:\/\/pctk.jhu.edu\/index.php\/wp-json\/wp\/v2\/pages\/35\/revisions\/67"}],"wp:attachment":[{"href":"https:\/\/pctk.jhu.edu\/index.php\/wp-json\/wp\/v2\/media?parent=35"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}