WO2013175352A1 - Fast scatter estimation in pet reconstruction. - Google Patents
Fast scatter estimation in pet reconstruction. Download PDFInfo
- Publication number
- WO2013175352A1 WO2013175352A1 PCT/IB2013/053965 IB2013053965W WO2013175352A1 WO 2013175352 A1 WO2013175352 A1 WO 2013175352A1 IB 2013053965 W IB2013053965 W IB 2013053965W WO 2013175352 A1 WO2013175352 A1 WO 2013175352A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- scatter
- sinogram
- imaging
- image processing
- scanner
- Prior art date
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
Definitions
- the present application relates to the diagnostic imaging arts. It finds particular application in accelerating the estimation of scatter in a nuclear medicine scanner, and will be described with particular reference thereto. It is to be understood, however, that it also finds application in other usage scenarios, and is not necessarily limited to the aforementioned application.
- SSS single-scatter simulation
- a typical method to estimate the scaling factor for SSS is fitting the "tail" part of a SSS sinogram to the measured sinogram, where tail refers to the portion in the sinogram corresponding to the outside of the imaged object.
- tail refers to the portion in the sinogram corresponding to the outside of the imaged object.
- the tail part in the measured data includes only contributions from scattered events. This assumption is valid for smaller patients when the tail is available and includes enough counts in it.
- the tail part decreases in size or disappears (truncated).
- fitting the tail can have significant error.
- the resulting image 2 typically suffers from scatter over-subtraction near high concentration areas 4 for a large patient.
- Another approach for accurate scatter estimation is to perform a full Monte Carlo simulation to produce shapes of both primary and scatter contributions. It demands significant amount of computations and therefore is too slow for a commercial nuclear medicine imaging system.
- the present application provides new and improved methods and systems which overcome the above-referenced problems and others.
- an image processing apparatus a scatter simulation processor which processes measured sinograms generated from imaging data acquired for an imaging subject by an imaging apparatus to produce a scatter sinogram that represents a shape of scatter contribution.
- a scatter scaling processor utilizes a Monte Carlo simulation to determine a scatter fraction and scales the scatter sinogram to generate a scaled scatter sinogram that matches the scatter contribution in the measured sinogram.
- a reconstruction processor reconstructs the imaging data into an image representation using the scaled scatter sinogram for scatter correction.
- a method of image processing including processing measured sinograms generated from imaging data acquired for an imaging subject by an imaging apparatus, producing a scatter sinogram that represents a shape of scatter contribution, determining a scatter fraction and scaling the scatter sinogram to generate a scaled scatter sinogram that matches the scatter contribution in the measured sinogram, and reconstructing the imaging data into an image representation using the scaled scatter sinogram for scatter correction.
- a method of image processing including with a PET scanner, generating a plurality of events, utilizing SSS and a short Monte Carlo simulation to determine the probability for a detected event pair having encountered scattering, and reconstructing the plurality of generated events into an image representation.
- One advantage resides in more accurate estimation of scatter contribution in nuclear medicine scanners. Another advantage resides in quicker estimation of scatter in nuclear medicine scanners.
- Another advantage resides in improved image quality and quantitative accuracy for nuclear medicine scanners.
- Another advantage resides in eliminating tail fitting in scatter contribution estimation.
- Another advantage resides in quicker image reconstruction.
- Another advantage resides in improved subject throughput.
- the invention may take form in various components and arrangements of components, and in various steps and arrangements of steps.
- the drawings are only for purposes of illustrating the preferred embodiments and are not to be construed as limiting the invention.
- FIGURE 1 is an exemplary prior art PET image reconstruction with SSS scatter estimation utilizing tail fitting.
- FIGURE 2 is an exemplary PET image reconstruction with SSS scatter estimation utilizing Monte Carlo simulation scaling in accordance with the present application.
- FIGURE 3 is a diagrammatic illustration of an imaging system in accordance with the present application.
- FIGURE 4 is a diagrammatic illustration of the comparison of profiles from measures and simulated data in accordance with the present application.
- FIGURE 5 is a flowchart illustration of a method of image processing in accordance with the present application.
- the present invention utilizes Monte Carlo simulations to quickly derive the ratio among primary and scattered photon pairs in measured PET data for each specific patient.
- the Monte Carlo simulation is based on the source distribution and attenuation map.
- the derived ratio is used to scale the sinogram generated from the single-scatter simulation (SSS) so that the scaled sinogram closely approximates the total scatter in the measured sinogram.
- SSS single-scatter simulation
- the resulting image 6 reconstructed using SSS with scaling factor derived from a short Monte Carlo simulation eliminates the scanner over-subtraction problem seen in Figure 1.
- a multi-modality system 10 includes a first imaging system, e.g. a functional modality, preferably, a nuclear imaging system 12, and a second imaging system, e.g. an anatomical modality, such as a computed tomography (CT) scanner 14.
- the CT scanner 14 includes a non-rotating gantry 16.
- An x-ray tube 18 is mounted to a rotating gantry 20.
- a bore 22 defines an examination region 24 of the CT scanner 14.
- An array of radiation detectors 26 is disposed on the rotating gantry 20 to receive radiation from the x-ray tube 18 after the x-rays transverse the examination region 24.
- the array of detectors 26 may be positioned on the non-rotating gantry 16.
- magnetic resonance and other imaging modalities are also contemplated.
- the functional or nuclear imaging system 12 includes a positron emission tomography (PET) scanner 30 which may mounted on tracks 32 to facilitate patient access.
- PET positron emission tomography
- the tracks 32 extend in parallel to a longitudinal axis of a subject support or couch 34, thus enabling the CT scanner 14 and PET scanner 12 to form a closed system.
- a motor and drive 36 is provided to move the PET scanner 12 in and out of the closed position.
- Detectors 38 are arranged around a bore 40 which defines an examination region 42. In the illustrated PET system, the detectors 38 are arranged in a stationery ring, although rotatable heads are also contemplated.
- the detectors 38 are typically incorporated into individual heads, which are mounted for rotational and radial movement relative to the patient.
- a motor and drive 44 or the like provides a longitudinal movement and vertical adjustment of the subject support 34 in the examination regions 24, 42. Mounted CT and PET systems in a single, shared close system with a common examination region is also contemplated.
- the subject support 34 which carries a subject, is moved into the examination region 24 of the CT scanner 14.
- the CT scanner 14 generates radiation attenuated data which is then used by an attenuation reconstruction processor 60 to reconstruct the radiation attenuated data into an attenuation map that is stored in an attenuation map memory 62.
- the patient support 34 moves the subject into the PET scanner 12 in a position that is geometrically and mechanically predicated as being the same as the imaged position in the CT imaging region 24.
- a subject is injected with a radiopharmaceutical.
- PET scanning a pair of gamma rays is produced by a positron annihilation event in the examination region 42 and travel in opposite directions.
- the gamma ray strikes the detectors 38, the location of the struck detector element and the strike time are recorded.
- a triggering processor 52 monitors each detector 38 for an energy spike, e.g., integrated area under the pulse, characteristic of the energy of the gamma rays generated by the radiopharmaceutical.
- each event is identified as an unscattered or scattered event.
- the triggering processor 52 checks a clock 54 and time stamps each detected gamma ray with a time of leading edge receipt stamp.
- the time stamp, energy estimate and detector position estimation are first used by an event verification processor 56 to determine whether there is a coincident event. Accepted pairs of coincident events define lines of response (LORs). Once an event pair is verified by the event verification processor 56, the LOR is passed to an event storage buffer 58 with their time stamps and stored in a list in the event storage buffer 58 as event data, i.e. as list-mode data.
- a sinogram reconstruction processor 64 reconstructs the verified pairs into an image representation of the subject.
- the sinogram reconstruction processor 64 converts the verified pairs into sinograms and accesses the attenuation data stored in the attenuation map memory 62 and reconstructs the sinograms into attenuation corrected source distribution map.
- the attenuation corrected source distribution map is stored in a source map memory 66.
- other reconstruction algorithms may be used including algorithms operating directly with the list-mode data such as list-mode ordered subsets expectation maximization (OSEM), and list-mode reconstruction with time-of- flight (TOF) reconstruction, etc.
- a scatter simulation processor 68 utilizes the determined source distribution map and the attenuation map to generate the shape of the scatter sinogram.
- a scatter scaling processor 72 also utilizes the determined source distribution map and the attenuation map to generate a ratio of scatter contribution with a short Monte Carlo simulation. In one embodiment, the Monte Carlo simulation is performed until the calculated ratio stabilizes. Once the ratio is stable, the Monte Carlo simulation can be terminated to reduce processing time.
- the SSS sinogram is scaled to define a scaled scatter sinogram.
- a reconstruction processor 74 utilizes scaled scatter sinogram and reconstructs the event pairs into a final reconstructed image with attenuation and scatter correction. The final reconstructed image is stored in an image memory 76 and displayed for a user on a display device 78, printed, saved for later use, and the like.
- the measured event data (after random correction) includes primary and scattered coincidence events.
- the single-scatter simulation (SSS) in scatter simulation processor 68 models single scattering from a given source distribution and the corresponding attenuation map.
- the scatter scaling processor 72 determines a scale factor with Monte Carlo simulation using the source distribution and the corresponding attenuation map. In the Monte Carlo simulation plurality of photon pairs are generated according to the source distribution. The trajectory of each photon in the attenuation map is traced until the photon escapes from the imaged object. The escaped photon may hit a detector and therefore is detected. A coincidence event is registered if both photons from a positron annihilation are detected.
- the event is labeled as primary if neither of the photons encountered scattering.
- the event is labeled as scattered if one or both of the photons have encountered one or more Compton scatterings in the attenuation medium.
- the ratio of total detected scattered events and total detected events represents the scatter fraction.
- the scatter fraction changes with number of positron annihilations and eventually stabilizes to a certain value r sc .
- the scatter fraction obtained from the Monte Carlo simulation is a good approximation of the actual scatter fraction in the measured coincidence events.
- a scaling factor for SSS sinogram is then obtained by the following formula: * measured!
- TssS ( 1 ) , where T meaS ured is the total counts in the measured sinogram and T sss is the total counts in the SSS sinogram.
- the SSS sinogram is scaled by k to produce an estimated scatter sinogram.
- the reconstruction processor 74 utilizes the scaled SS S sinogram and reconstructs the final reconstructed image.
- Illustrated in Figure 4 are the line profiles 80 from the sinograms obtained from the measured data (after random correction), Monte Carlo simulated data, Monte Carlo simulated scattering data, and SSS scattering data after scaled using the scaling factor described above. As shown, (1) the Monte Carlo simulation result matches the measured data and (2) the SSS simulation matches with the total scatter derived from Monte Carlo simulation.
- the triggering processor 52, event verification processor 56, attenuation reconstruction processor 60, sinogram reconstruction processor 64, scatter simulation processor 68, and scatter scaling processor 72 include a processor, for example a microprocessor or other software controlled device configured to execute software for performing the operations described above.
- the software is carried on tangible memory or a computer readable medium for execution by the processor.
- Types of computer readable media include memory such as a hard disk drive, CD-ROM, DVD- ROM and the like.
- Other implementations of the processor are also contemplated.
- Display controllers, Application Specific Integrated Circuits (ASICs), FPGAs, and microcontrollers are illustrative examples of other types of component which may be implemented to provide functions of the processor.
- Embodiments may be implemented using software for execution by a processor, hardware, or some combination thereof.
- FIGURE 5 illustrates a method of image processing.
- a step 100 attenuation data is received from an examination region.
- An attenuation map is reconstructed from the received attenuation data in a step 102.
- a step 104 potential radiation events are received from an examination region.
- a time stamp is assigned to the received events in a step 106.
- a verification criterion is applied to detect the coincident events.
- verified pairs of coincident events are defined as LORs.
- a reconstruction processor converts the verified pairs into sinograms.
- a reconstruction processor reconstructs the sinograms into an attenuation corrected source distribution map.
- a scatter simulation processor utilizes the determined source distribution map and the attenuation map to generate the shape of the scatter sinogram in a step 116.
- a scatter scaling processor utilizes the determined source distribution map and the attenuation map to generate a ratio of scatter contribution with a short Monte Carlo simulation, and scales the scatter sinogram to match scatter contribution in the measured sinogram.
- a reconstruction processor reconstructs a final reconstructed image utilizing the scaled scatter sinogram and the measured coincidence events or sinograms.
Landscapes
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Nuclear Medicine (AREA)
- Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Life Sciences & Earth Sciences (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Image Analysis (AREA)
Abstract
Description
Claims
Priority Applications (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2015513315A JP6170549B2 (en) | 2012-05-21 | 2013-05-15 | Rapid scatter estimation in PET reconstruction |
RU2014151735A RU2620862C2 (en) | 2012-05-21 | 2013-05-15 | Quick assessment of scattering at reconstruction by positron-emission tomography |
CN201380026774.5A CN104335247B (en) | 2012-05-21 | 2013-05-15 | Apparatus and method for the quick scattering estimation in being rebuild in PET |
US14/396,464 US9659389B2 (en) | 2012-05-21 | 2013-05-15 | Fast scatter estimation in pet reconstruction |
EP13735422.1A EP2852931B1 (en) | 2012-05-21 | 2013-05-15 | Fast scatter estimation in pet reconstruction |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201261649825P | 2012-05-21 | 2012-05-21 | |
US61/649,825 | 2012-05-21 | ||
US201261658049P | 2012-06-11 | 2012-06-11 | |
US61/658,049 | 2012-06-11 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2013175352A1 true WO2013175352A1 (en) | 2013-11-28 |
Family
ID=48782557
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/IB2013/053965 WO2013175352A1 (en) | 2012-05-21 | 2013-05-15 | Fast scatter estimation in pet reconstruction. |
Country Status (6)
Country | Link |
---|---|
US (1) | US9659389B2 (en) |
EP (1) | EP2852931B1 (en) |
JP (1) | JP6170549B2 (en) |
CN (1) | CN104335247B (en) |
RU (1) | RU2620862C2 (en) |
WO (1) | WO2013175352A1 (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018163362A1 (en) | 2017-03-09 | 2018-09-13 | 株式会社島津製作所 | Scattering estimation method, scattering estimation program, and positron ct device having same installed thereon |
US10223480B2 (en) | 2013-04-11 | 2019-03-05 | Koninklijke Philips N.V. | Method for modeling and accounting for cascade gammas in images |
CN110009707A (en) * | 2019-04-09 | 2019-07-12 | 上海联影医疗科技有限公司 | Image dispersion bearing calibration, device, computer equipment and storage medium |
CN110211198A (en) * | 2019-06-03 | 2019-09-06 | 上海联影医疗科技有限公司 | Image dispersion bearing calibration, device, computer equipment and storage medium |
CN111627084A (en) * | 2020-05-29 | 2020-09-04 | 上海联影医疗科技有限公司 | Image correction method, device and system based on PET scanning |
US11234667B2 (en) | 2017-09-14 | 2022-02-01 | Koninklijke Philips N.V. | Scatter correction using emission image estimate reconstructed from narrow energy window counts in positron emission tomography |
Families Citing this family (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10034640B2 (en) * | 2015-03-11 | 2018-07-31 | David Byron Douglas | Method and apparatus for improvement of spatial resolution in molecular and radiological imaging |
US10036817B2 (en) * | 2015-05-04 | 2018-07-31 | Koninklijke Philips N.V. | Solving outside-field of view scatter correction problem in positron emission tomography via digital experimentation |
EP3297537B1 (en) * | 2015-05-19 | 2019-01-30 | Koninklijke Philips N.V. | Estimation of an attenuation map based on scattered coincidences in a pet system |
JP6440230B2 (en) * | 2015-08-07 | 2018-12-19 | シャンハイ・ユナイテッド・イメージング・ヘルスケア・カンパニー・リミテッド | Multi-modality imaging system and method |
CN105631910B (en) | 2015-12-28 | 2019-02-05 | 沈阳东软医疗系统有限公司 | A kind of method for reconstructing and device of PET image |
DE102016217984A1 (en) * | 2016-09-20 | 2018-04-05 | Siemens Healthcare Gmbh | Sinogram-based scatter correction in computed tomography |
CN106491153B (en) * | 2016-12-29 | 2017-10-27 | 上海联影医疗科技有限公司 | A kind of PET scatter correction methods, PET imaging methods and PET imaging systems |
CN107133996B (en) * | 2017-03-21 | 2020-08-04 | 上海联影医疗科技有限公司 | Method for generating an attenuation map for PET data reconstruction and PET/CT system |
WO2018202648A1 (en) * | 2017-05-01 | 2018-11-08 | Koninklijke Philips N.V. | Generation of accurate hybrid datasets for quantitative molecular imaging |
US11002866B2 (en) * | 2017-07-25 | 2021-05-11 | Shimadzu Corporation | Scattering estimation method and image processor |
JP6984010B2 (en) * | 2017-09-28 | 2021-12-17 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | Deep learning-based scattering correction |
CN111670462B (en) * | 2018-01-31 | 2024-08-27 | 皇家飞利浦有限公司 | Scatter correction for Positron Emission Tomography (PET) |
CN109115809A (en) * | 2018-09-04 | 2019-01-01 | 南京航空航天大学 | A kind of γ photon equilibrium state meets the method for reverting to true coincidence |
CN109124672B (en) * | 2018-10-19 | 2021-08-17 | 上海联影医疗科技股份有限公司 | Method, device, equipment and storage medium for correcting random coincidence events |
JP7534038B2 (en) * | 2018-11-30 | 2024-08-14 | アキュレイ インコーポレイテッド | Asymmetric scatter fitting for optimal panel readout in cone-beam computed tomography |
CN110136076B (en) * | 2019-04-18 | 2021-07-16 | 上海联影医疗科技股份有限公司 | Medical scanning imaging method, device, storage medium and computer equipment |
CN110097611B (en) * | 2019-04-28 | 2023-09-22 | 上海联影智能医疗科技有限公司 | Image reconstruction method, device, equipment and storage medium |
CN110428384B (en) * | 2019-08-08 | 2021-11-16 | 江苏赛诺格兰医疗科技有限公司 | Method for acquiring correction information for attenuation correction of PET images of respiration or heart |
CN110477937B (en) * | 2019-08-26 | 2023-07-25 | 上海联影医疗科技股份有限公司 | Scattering estimation parameter determination method, device, equipment and medium |
CN111568450B (en) * | 2020-05-20 | 2023-04-18 | 上海联影医疗科技股份有限公司 | PET scanning data scattering correction method and device and computer equipment |
CN112598760B (en) * | 2020-12-18 | 2023-07-04 | 上海联影医疗科技股份有限公司 | Image truncation artifact correction method and system |
Family Cites Families (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6163589A (en) | 1998-06-13 | 2000-12-19 | General Electric Company | Monte Carlo scatter correction method for computed tomography of general object geometries |
JP4536212B2 (en) * | 2000-05-24 | 2010-09-01 | 浜松ホトニクス株式会社 | PET equipment |
DE10047720A1 (en) * | 2000-09-27 | 2002-04-11 | Philips Corp Intellectual Pty | Device and method for generating an X-ray computer tomogram with a scattered radiation correction |
US7471813B2 (en) * | 2004-10-01 | 2008-12-30 | Varian Medical Systems International Ag | Systems and methods for correction of scatter in images |
GB0512252D0 (en) * | 2005-06-16 | 2005-07-27 | Hammersmith Imanet Ltd | Method of and software for conducting motion correction in tomographic scanning and system for tomographic scanning using the method |
US7569827B2 (en) * | 2005-08-16 | 2009-08-04 | Chuanyong Bai | Emission-data-based photon scatter correction in computed nuclear imaging technology |
WO2007035775A2 (en) * | 2005-09-19 | 2007-03-29 | Feng Ma | Imaging system and method utilizing primary radiation |
US7397035B2 (en) | 2005-10-14 | 2008-07-08 | Siemens Medical Solutions Usa, Inc. | Scatter correction for time-of-flight positron emission tomography data |
WO2007054843A1 (en) * | 2005-11-10 | 2007-05-18 | Koninklijke Philips Electronics, N.V. | Pet imaging using anatomic list mode mask |
WO2008075037A1 (en) * | 2006-12-18 | 2008-06-26 | Hammersmith Imanet Limited | Method of and software for calculating a scatter estimate for tomographic scanning and system for tomographic scanning |
JP2011509697A (en) * | 2007-12-04 | 2011-03-31 | ジーイー・ヘルスケア・リミテッド | Image analysis |
EP2289048B1 (en) * | 2008-05-15 | 2019-02-27 | Koninklijke Philips N.V. | Using non-attenuation corrected pet emission images to compensate for incomplete anatomic images |
US8193505B2 (en) | 2008-09-29 | 2012-06-05 | Siemens Medical Solutions Usa, Inc. | System and method for scatter normalization of PET images |
CN101499173B (en) * | 2009-03-06 | 2011-03-09 | 刘华锋 | Kalman filtering image reconstruction method in PET imaging |
WO2010133985A2 (en) | 2009-05-20 | 2010-11-25 | Koninklijke Philips Electronics, N.V. | Continuous time-of-flight scatter simulation method |
US20110142367A1 (en) * | 2009-12-15 | 2011-06-16 | Charles William Stearns | Methods and systems for correcting image scatter |
JP5611640B2 (en) * | 2010-04-05 | 2014-10-22 | 株式会社東芝 | Nuclear medicine diagnostic apparatus, medical image processing apparatus, and medical image diagnostic apparatus |
RU2014148796A (en) * | 2012-05-04 | 2016-06-27 | Конинклейке Филипс Н.В. | MISCELLANEOUS MAP OF THE DISTINCTED MATCHES IN POSITRON-EMISSION TOMOGRAPHY |
-
2013
- 2013-05-15 US US14/396,464 patent/US9659389B2/en active Active
- 2013-05-15 CN CN201380026774.5A patent/CN104335247B/en active Active
- 2013-05-15 RU RU2014151735A patent/RU2620862C2/en not_active IP Right Cessation
- 2013-05-15 EP EP13735422.1A patent/EP2852931B1/en active Active
- 2013-05-15 WO PCT/IB2013/053965 patent/WO2013175352A1/en active Application Filing
- 2013-05-15 JP JP2015513315A patent/JP6170549B2/en active Active
Non-Patent Citations (4)
Title |
---|
HOLDSWORTH C H ET AL: "Evaluation of a Monte Carlo scatter correction in clinical 3D PET", 2003 IEEE NUCLEAR SCIENCE SYMPOSIUM CONFERENCE RECORD. / 2003 IEEE NUCLEAR SCIENCE SYMPOSIUM AND MEDICAL IMAGING CONFERENCE. PORTLAND, OR, OCT. 19 - 25, 2003; [IEEE NUCLEAR SCIENCE SYMPOSIUM CONFERENCE RECORD], NEW YORK, NY : IEEE, US, vol. 4, 19 October 2003 (2003-10-19), pages 2540 - 2544, XP010737563, ISBN: 978-0-7803-8257-2 * |
IRENE POLYCARPOU ET AL: "Comparative evaluation of scatter correction in 3D PET using different scatter-level approximations", ANNALS OF NUCLEAR MEDICINE, SPRINGER JAPAN, JAPAN, vol. 25, no. 9, 14 July 2011 (2011-07-14), pages 643 - 649, XP019978843, ISSN: 1864-6433, DOI: 10.1007/S12149-011-0514-Y * |
KYUNG SANG KIM ET AL: "Ultra-fast hybrid CPU-GPU Monte Carlo simulation for scatter correction in 3D PETs", NUCLEAR SCIENCE SYMPOSIUM AND MEDICAL IMAGING CONFERENCE (NSS/MIC), 2011 IEEE, IEEE, 23 October 2011 (2011-10-23), pages 2749 - 2752, XP032121241, ISBN: 978-1-4673-0118-3, DOI: 10.1109/NSSMIC.2011.6152961 * |
READER A J ET AL: "Attenuation and scatter correction of list-mode data driven iterative and analytic image reconstruction algorithms for rotating 3D PET systems", IEEE TRANSACTIONS ON NUCLEAR SCIENCE, IEEE SERVICE CENTER, NEW YORK, NY, US, vol. 46, no. 6, 1 December 1999 (1999-12-01), pages 2218 - 2226, XP011358669, ISSN: 0018-9499, DOI: 10.1109/23.819307 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10223480B2 (en) | 2013-04-11 | 2019-03-05 | Koninklijke Philips N.V. | Method for modeling and accounting for cascade gammas in images |
WO2018163362A1 (en) | 2017-03-09 | 2018-09-13 | 株式会社島津製作所 | Scattering estimation method, scattering estimation program, and positron ct device having same installed thereon |
US11513243B2 (en) | 2017-03-09 | 2022-11-29 | Shimadzu Corporation | Scatter estimation method, scatter estimation program, and positron CT device having same installed thereon |
US11234667B2 (en) | 2017-09-14 | 2022-02-01 | Koninklijke Philips N.V. | Scatter correction using emission image estimate reconstructed from narrow energy window counts in positron emission tomography |
CN110009707A (en) * | 2019-04-09 | 2019-07-12 | 上海联影医疗科技有限公司 | Image dispersion bearing calibration, device, computer equipment and storage medium |
CN110211198A (en) * | 2019-06-03 | 2019-09-06 | 上海联影医疗科技有限公司 | Image dispersion bearing calibration, device, computer equipment and storage medium |
CN110211198B (en) * | 2019-06-03 | 2023-01-03 | 上海联影医疗科技股份有限公司 | Image scattering correction method, device, computer equipment and storage medium |
CN111627084A (en) * | 2020-05-29 | 2020-09-04 | 上海联影医疗科技有限公司 | Image correction method, device and system based on PET scanning |
CN111627084B (en) * | 2020-05-29 | 2023-11-14 | 上海联影医疗科技股份有限公司 | Image correction method, device and system based on PET scanning |
Also Published As
Publication number | Publication date |
---|---|
JP6170549B2 (en) | 2017-07-26 |
CN104335247A (en) | 2015-02-04 |
EP2852931B1 (en) | 2019-05-15 |
RU2620862C2 (en) | 2017-05-30 |
US20150093004A1 (en) | 2015-04-02 |
EP2852931A1 (en) | 2015-04-01 |
US9659389B2 (en) | 2017-05-23 |
RU2014151735A (en) | 2016-07-10 |
JP2015522803A (en) | 2015-08-06 |
CN104335247B (en) | 2018-03-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9659389B2 (en) | Fast scatter estimation in pet reconstruction | |
US10223480B2 (en) | Method for modeling and accounting for cascade gammas in images | |
Zaidi et al. | Scatter compensation techniques in PET | |
US9155514B2 (en) | Reconstruction with partially known attenuation information in time of flight positron emission tomography | |
JP5254810B2 (en) | Local motion compensation based on list mode data | |
EP3494547B1 (en) | Time-of-flight (tof) pet image reconstruction using locally modified tof kernels | |
CN107635469B (en) | Estimation of attenuation map based on scatter coincidence in PET system | |
US20090124900A1 (en) | Three-Dimensional Time-of-Flight Pet With Course Angular and Slice Rebinning | |
US8193505B2 (en) | System and method for scatter normalization of PET images | |
CN103315760A (en) | Systems and methods for attenuation compensation in nuclear medicine imaging based on emission data | |
EP3619686B1 (en) | Generation of accurate hybrid datasets for quantitative molecular imaging | |
US20140107476A1 (en) | Spatial sampling improvement for list-mode pet acquisition using planned table/gantry movement | |
EP3292426A1 (en) | Solving outside-field of view scatter correction problem in positron emission tomography via digital experimentation | |
US11324472B2 (en) | Energy-based scatter correction for PET sinograms | |
HU231327B1 (en) | Multiple emission energies in single photon emission computed tomography | |
US9317942B2 (en) | Fast computation of the sensitivity matrix in iterative algorithms | |
KR20200057450A (en) | Method and system for generating virtual CT(Computed Tomography) image and attenuation-corrected PET(Positron Emission Tomography) image based on PET image | |
Fin et al. | A practical way to improve contrast‐to‐noise ratio and quantitation for statistical‐based iterative reconstruction in whole‐body PET imaging | |
Silva-Rodriguez et al. | A SimSET-STIR hybrid Monte Carlo model for the Philips Vereos digital PET | |
TW201417769A (en) | A method for improving image quality and imaging system using the same | |
Laymon et al. | Anomaly detection and artifact recovery in pet attenuation-correction images using the likelihood function | |
Rankine | The Effects of PET Reconstruction Parameters on Radiotherapy Response Assessment and an Investigation of SUV-peak Sampling Parameters | |
Case et al. | Attenuation and scatter correction in cardiac PET |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 13735422 Country of ref document: EP Kind code of ref document: A1 |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2013735422 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 14396464 Country of ref document: US |
|
ENP | Entry into the national phase |
Ref document number: 2015513315 Country of ref document: JP Kind code of ref document: A |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
ENP | Entry into the national phase |
Ref document number: 2014151735 Country of ref document: RU Kind code of ref document: A |