Estimating the NEMA characteristics of the J-PET tomograph using the GATE package

Feasibility of the J-PET to monitor the range of therapeutic proton beams

PURPOSE

The aim of this work is to investigate the feasibility of the Jagiellonian Positron Emission Tomography (J-PET) scanner for intra-treatment proton beam range monitoring.

METHODS

The Monte Carlo simulation studies with GATE and PET image reconstruction with CASToR were performed in order to compare six J-PET scanner geometries. We simulated proton irradiation of a PMMA phantom with a Single Pencil Beam (SPB) and Spread-Out Bragg Peak (SOBP) of various ranges. The sensitivity and precision of each scanner were calculated, and considering the setup's cost-effectiveness, we indicated potentially optimal geometries for the J-PET scanner prototype dedicated to the proton beam range assessment.

RESULTS

The investigations indicate that the double-layer cylindrical and triple-layer double-head configurations are the most promising for clinical application. We found that the scanner sensitivity is of the order of 10-5 coincidences per primary proton, while the precision of the range assessment for both SPB and SOBP irradiation plans was found below 1 mm. Among the scanners with the same number of detector modules, the best results are found for the triple-layer dual-head geometry. The results indicate that the double-layer cylindrical and triple-layer double-head configurations are the most promising for the clinical application, CONCLUSIONS:: We performed simulation studies demonstrating that the feasibility of the J-PET detector for PET-based proton beam therapy range monitoring is possible with reasonable sensitivity and precision enabling its pre-clinical tests in the clinical proton therapy environment. Considering the sensitivity, precision and cost-effectiveness, the double-layer cylindrical and triple-layer dual-head J-PET geometry configurations seem promising for future clinical application.

Validation of a discovery MI 4-ring model according to the NEMA NU 2-2018 standards: from Monte Carlo simulations to clinical-like reconstructions

BACKGROUND

We propose a comprehensive evaluation of a Discovery MI 4-ring (DMI) model, using a Monte Carlo simulator (GATE) and a clinical reconstruction software package (PET toolbox). The following performance characteristics were compared with actual measurements according to NEMA NU 2-2018 guidelines: system sensitivity, count losses and scatter fraction (SF), coincidence time resolution (CTR), spatial resolution (SR), and image quality (IQ). For SR and IQ tests, reconstruction of time-of-flight (TOF) simulated data was performed using the manufacturer's reconstruction software.

RESULTS

Simulated prompt, random, true, scatter and noise equivalent count rates closely matched the experimental rates with maximum relative differences of 1.6%, 5.3%, 7.8%, 6.6%, and 16.5%, respectively, in a clinical range of less than 10 kBq/mL. A 3.6% maximum relative difference was found between experimental and simulated sensitivities. The simulated spatial resolution was better than the experimental one. Simulated image quality metrics were relatively close to the experimental results.

CONCLUSIONS

The current model is able to reproduce the behaviour of the DMI count rates in the clinical range and generate clinical-like images with a reasonable match in terms of contrast and noise.

Comparative studies of the sensitivities of sparse and full geometries of Total-Body PET scanners built from crystals and plastic scintillators

BACKGROUND

Alongside the benefits of Total-Body imaging modalities, such as higher sensitivity, single-bed position, low dose imaging, etc., their final construction cost prevents worldwide utilization. The main aim of this study is to present a simulation-based comparison of the sensitivities of existing and currently developed tomographs to introduce a cost-efficient solution for constructing a Total-Body PET scanner based on plastic scintillators.

METHODS

For the case of this study, eight tomographs based on the uEXPLORER configuration with different scintillator materials (BGO, LYSO), axial field-of-view (97.4 cm and 194.8 cm), and detector configurations (full and sparse) were simulated. In addition, 8 J-PET scanners with different configurations, such as various axial field-of-view (200 cm and 250 cm), different cross sections of plastic scintillator, and multiple numbers of plastic scintillator layers (2, 3, and 4), based on J-PET technology have been simulated by GATE software. Furthermore, Siemens' Biograph Vision has been simulated to compare the results with standard PET scans. Two types of simulations have been performed. The first one with a centrally located source with a diameter of 1 mm and a length of 250 cm, and the second one with the same source inside a water-filled cylindrical phantom with a diameter of 20 cm and a length of 183 cm.

RESULTS

With regards to sensitivity, among all the proposed scanners, the ones constructed with BGO crystals give the best performance ([Formula: see text] 350 cps/kBq at the center). The utilization of sparse geometry or LYSO crystals significantly lowers the achievable sensitivity of such systems. The J-PET design gives a similar sensitivity to the sparse LYSO crystal-based detectors while having full detector coverage over the body. Moreover, it provides uniform sensitivity over the body with additional gain on its sides and provides the possibility for high-quality brain imaging.

CONCLUSION

Taking into account not only the sensitivity but also the price of Total-Body PET tomographs, which till now was one of the main obstacles in their widespread clinical availability, the J-PET tomography system based on plastic scintillators could be a cost-efficient alternative for Total-Body PET scanners.

Technical opportunities and challenges in developing total-body PET scanners for mice and rats

Positron emission tomography (PET) is the most sensitive in vivo molecular imaging technique available. Small animal PET has been widely used in studying pharmaceutical biodistribution and disease progression over time by imaging a wide range of biological processes. However, it remains true that almost all small animal PET studies using mouse or rat as preclinical models are either limited by the spatial resolution or the sensitivity (especially for dynamic studies), or both, reducing the quantitative accuracy and quantitative precision of the results. Total-body small animal PET scanners, which have axial lengths longer than the nose-to-anus length of the mouse/rat and can provide high sensitivity across the entire body of mouse/rat, can realize new opportunities for small animal PET. This article aims to discuss the technical opportunities and challenges in developing total-body small animal PET scanners for mice and rats.

Multi-photon time-of-flight MLEM application for the positronium imaging in J-PET

We develop a positronium imaging method for the Jagiellonian PET (J-PET) scanners based on the time-of-flight maximum likelihood expectation maximisation (TOF MLEM). The system matrix elements are calculated on-the-fly for the coincidences comprising two annihilation and one de-excitation photons that originate from the ortho-positronium (o-Ps) decay. Using the Geant4 library, a Monte Carlo simulation was conducted for four cylindrical 22Na sources of β+ decay with diverse o-Ps mean lifetimes, placed symmetrically inside the two JPET prototypes. The estimated time differences between the annihilation and the positron emission were aggregated into histograms (one per voxel), updated by the weights of the activities reconstructed by TOF MLEM. The simulations were restricted to include only the o-Ps decays into back-to-back photons, allowing a linear fitting model to be employed for the estimation of the mean lifetime from each histogram built in the log scale. To suppress the noise, the exclusion of voxels with activity below 2% – 10% of the peak was studied. The estimated o-Ps mean lifetimes were consistent with the simulation and distributed quasi -uniformly at high MLEM iterations. The proposed positronium imaging technique can be further upgraded to include various correction factors, as well as be modified according to realistic o-Ps decay models.

ProTheRaMon - a GATE simulation framework for proton therapy range monitoring using PET imaging

OBJECTIVE

This paper reports on the implementation and shows examples of the use of the ProTheRaMon framework for simulating the delivery of proton therapy treatment plans and range monitoring using positron emission tomography (PET). ProTheRaMon offers complete processing of proton therapy treatment plans, patient CT geometries, and intra-treatment PET imaging, taking into account therapy and imaging coordinate systems and activity decay during the PET imaging protocol specific to a given proton therapy facility. We present the ProTheRaMon framework and illustrate its potential use case and data processing steps for a patient treated at the Cyclotron Centre Bronowice (CCB) proton therapy center in Krakow, Poland.

APPROACH

The ProTheRaMon framework is based on GATE Monte Carlo software, the CASToR reconstruction package and in-house developed Python and bash scripts. The framework consists of five separated simulation and data processing steps, that can be further optimized according to the user's needs and specific settings of a given proton therapy facility and PET scanner design.

MAIN RESULTS

ProTheRaMon is presented using example data from a patient treated at CCB and the J-PET scanner to demonstrate the application of the framework for proton therapy range monitoring. The output of each simulation and data processing stage is described and visualized.

SIGNIFICANCE

We demonstrate that the ProTheRaMon simulation platform is a high-performance tool, capable of running on a computational cluster and suitable for multi-parameter studies, with databases consisting of large number of patients, as well as different PET scanner geometries and settings for range monitoring in a clinical environment. Due to its modular structure, the ProTheRaMon framework can be adjusted for different proton therapy centers and/or different PET detector geometries. It is available to the community via github.

Simulating NEMA characteristics of the modular total-body J-PET scanner-an economic total-body PET from plastic scintillators

The purpose of the presented research is estimation of the performance characteristics of the economic total-body Jagiellonian-PET system (TB-J-PET) constructed from plastic scintillators. The characteristics are estimated according to the NEMA NU-2-2018 standards utilizing the GATE package. The simulated detector consists of 24 modules, each built out of 32 plastic scintillator strips (each with cross section of 6 mm times 30 mm and length of 140 or 200 cm) arranged in two layers in regular 24-sided polygon circumscribing a circle with the diameter of 78.6 cm. For the TB-J-PET with an axial field-of-view (AFOV) of 200 cm, a spatial resolutions (SRs) of 3.7 mm (transversal) and 4.9 mm (axial) are achieved. The noise equivalent count rate (NECR) peak of 630 kcps is expected at 30 kBq cc^-1. Activity concentration and the sensitivity at the center amounts to 38 cps kBq^-1. The scatter fraction (SF) is estimated to 36.2 %. The values of SF and SR are comparable to those obtained for the state-of-the-art clinical PET scanners and the first total-body tomographs: uExplorer and PennPET. With respect to the standard PET systems with AFOV in the range from 16 to 26 cm, the TB-J-PET is characterized by an increase in NECR approximately by factor of 4 and by the increase of the whole-body sensitivity by factor of 12.6 to 38. The time-of-flight resolution for the TB-J-PET is expected to be at the level of CRT = 240 ps full width at half maximum. For the TB-J-PET with an AFOV of 140 cm, an image quality of the reconstructed images of a NEMA IEC phantom was presented with a contrast recovery coefficient and a background variability parameters. The increase of the whole-body sensitivity and NECR estimated for the TB-J-PET with respect to current commercial PET systems makes the TB-J-PET a promising cost-effective solution for the broad clinical applications of total-body PET scanners. TB-J-PET may constitute an economic alternative for the crystal TB-PET scanners, since plastic scintillators are much cheaper than BGO or LYSO crystals and axial arrangement of the strips significantly reduces the costs of readout electronics and SiPMs.

Hybrid model of photon propagation based on the analytical and Monte Carlo methods for a dual-head PET system

The construction of photon propagation has a close relationship with the quality of reconstructed images. The classical Monte Carlo (MC) based method can model the photon propagation precisely, but it is time-consuming. The analytical method can often quickly construct a model, but its precision is a problem. How to fully exploit the advantages of the MC simulation and analytical model is an open problem. Inspired by the characteristics of the depth of interaction (DOI) detectors, which can help confirm the deposited position of a photon with DOI-encoding technology, we virtually discretize each crystal into several subcrystals to obtain the statistical distribution by MC-based simulation. Then, the statistical distribution is combined with a spatially variant solid-angle model. This combination strategy provides a hybrid model to describe photon propagation with relatively high accuracy and low computational cost. Three discretization schemes are compared to optimize the constructed photon propagation model. Several experiments are carried out to evaluate the performance of the proposed hybrid method. The metrics of full width at half maximum (FWHM), contrast recovery (CR), and coefficient of variation (COV) are adopted to quantitate the imaging results. The classical MC-based method is compared as a gold-standard reference. When a crystal is divided into two discretized positions, the convergent tendencies of CRs and COVs are consistent with that based on MC simulation method, respectively. In terms of FWHMs, the resolutions of using the MC-based model and the proposed hybrid model are 0.71 mm and 0.68 mm in the direction parallel to the detector head, respectively. This indicates the potential of the proposed method in positron emission tomography imaging.

Optimisation of the event-based TOF filtered back-projection for online imaging in total-body J-PET

We perform a parametric study of the newly developed time-of-flight (TOF) image reconstruction algorithm, proposed for the real-time imaging in total-body Jagiellonian PET (J-PET) scanners. The asymmetric 3D filtering kernel is applied at each most likely position of electron-positron annihilation, estimated from the emissions of back-to-back γ-photons. The optimisation of its parameters is studied using Monte Carlo simulations of a 1-mm spherical source, NEMA IEC and XCAT phantoms inside the ideal J-PET scanner. The combination of high-pass filters which included the TOF filtered back-projection (FBP), resulted in spatial resolution, 1.5 times higher in the axial direction than for the conventional 3D FBP. For realistic 10-minute scans of NEMA IEC and XCAT, which require a trade-off between the noise and spatial resolution, the need for Gaussian TOF kernel components, coupled with median post-filtering, is demonstrated. The best sets of 3D filter parameters were obtained by the Nelder-Mead minimisation of the mean squared error between the resulting and reference images. The approach allows training the reconstruction algorithm for custom scans, using the IEC phantom, when the temporal resolution is below 50 ps. The image quality parameters, estimated for the best outcomes, were systematically better than for the non-TOF FBP.

Advanced Monte Carlo simulations of emission tomography imaging systems with GATE

Built on top of the Geant4 toolkit, GATE is collaboratively developed for more than 15 years to design Monte Carlo simulations of nuclear-based imaging systems. It is, in particular, used by researchers and industrials to design, optimize, understand and create innovative emission tomography systems. In this paper, we reviewed the recent developments that have been proposed to simulate modern detectors and provide a comprehensive report on imaging systems that have been simulated and evaluated in GATE. Additionally, some methodological developments that are not specific for imaging but that can improve detector modeling and provide computation time gains, such as Variance Reduction Techniques and Artificial Intelligence integration, are described and discussed.

Sensitivity of Discrete Symmetry Tests in the Positronium System with the J-PET Detector

Study of certain angular correlations in the three-photon annihilations of the triplet state of positronium, the electron–positron bound state, may be used as a probe of potential CP and CPT-violating effects in the leptonic sector. We present the perspectives of CP and CPT tests using this process recorded with a novel detection system for photons in the positron annihilation energy range, the Jagiellonian Positron Emission Tomography (J-PET). We demonstrate the capability of this system to register three-photon annihilations with an unprecedented range of kinematical configurations and to measure the CPT-odd correlation between positronium spin and annihilation plane orientation with a precision improved by at least an order of magnitude with respect to present results. We also discuss the means to control and reduce detector asymmetries in order to allow J-PET to set the first measurement of the correlation between positronium spin and momentum of the most energetic annihilation photon which has never been studied to date.

Performance assessment of the 2 γpositronium imaging with the total-body PET scanners

Purpose

In living organisms, the positron-electron annihilation (occurring during the PET imaging) proceeds in about 30% via creation of a metastable ortho-positronium atom. In the tissue, due to the pick-off and conversion processes, over 98% of ortho-positronia annihilate into two 511 keV photons. In this article, we assess the feasibility for reconstruction of the mean ortho-positronium lifetime image based on annihilations into two photons. The main objectives of this work include the (i) estimation of the sensitivity of the total-body PET scanners for the ortho-positronium mean lifetime imaging using 2γ annihilations and (ii) estimation of the spatial and time resolution of the ortho-positronium image as a function of the coincidence resolving time (CRT) of the scanner.

Methods

Simulations are conducted assuming that radiopharmaceutical is labeled with ⁴⁴Sc isotope emitting one positron and one prompt gamma. The image is reconstructed on the basis of triple coincidence events. The ortho-positronium lifetime spectrum is determined for each voxel of the image. Calculations were performed for cases of total-body detectors build of (i) LYSO scintillators as used in the EXPLORER PET and (ii) plastic scintillators as anticipated for the cost-effective total-body J-PET scanner. To assess the spatial and time resolution, the four cases were considered assuming that CRT is equal to 500 ps, 140 ps, 50 ps, and 10 ps.

Results

The estimated total-body PET sensitivity for the registration and selection of image forming triple coincidences (2γ+γ_prompt) is larger by a factor of 13.5 (for LYSO PET) and by factor of 5.2 (for plastic PET) with respect to the sensitivity for the standard 2γ imaging by LYSO PET scanners with AFOV = 20 cm. The spatial resolution of the ortho-positronium image is comparable with the resolution achievable when using TOF-FBP algorithms already for CRT = 50 ps. For the 20-min scan, the resolution better than 20 ps is expected for the mean ortho-positronium lifetime image determination.

Conclusions

Ortho-positronium mean lifetime imaging based on the annihilations into two photons and prompt gamma is shown to be feasible with the advent of the high sensitivity total-body PET systems and time resolution of the order of tens of picoseconds.

Estimating relationship between the time over threshold and energy loss by photons in plastic scintillators used in the J-PET scanner

PURPOSE

The time-over-threshold (TOT) technique is being used widely due to itsimplications in developing the multi-channel readouts, mainly when fast signal processing is required. Using the TOT technique, as a measure of energy loss instead of charge integration methods, significantly reduces the signal readout costs by combining the time and energy information. Therefore, this approach can potentially be utilized in J-PET tomograph which is built from plastic scintillators characterized by fast light signals. The drawback in adopting this technique lies in the non-linear correlation between input energy loss and TOT of the signal. The main motivation behind this work is to develop the relationship between TOT and energy loss and validate it by the J-PET tomograph setup.

METHODS

The experiment was performed using a ²²Na beta emitter source placed in the center of the J-PET tomograph. This isotope produces photons of two different energies: 511 keV photons from the positron annihilation (direct annihilation or through the formation of a para-positronium atom or pick-off process of ortho-positronium atoms) and a 1275 keV prompt photon. This allows the study of the correlation between TOT values and energy loss for energy ranges up to 1000 keV. Since the photon interacts predominantly via Compton scattering inside the plastic scintillator, there is no direct information of the energy deposition. However, using the J-PET geometry, one can measure the scattering angle of the interacting photon. Since the ²²Na source emits photons of two different energies, it is necessary to know unambiguously the energy of incident photons and their corresponding scattering angles in order to estimate energy deposition. In summary, this work presents a dedicated algorithm developed to tag photons of different energies and studying their scattering angles to calculate the energy deposition by the interacting photons.

RESULTS

A new method was elaborated to measure the energy loss by photons interacting with plastic scintillators used in the J-PET tomograph. We find the relationship between the energy loss and TOT is non-linear and can be described by the functions TOT = A0 + A1 * ln(E _dep + A2) + A3 * (ln(E _dep + A2))² and TOT = A0 - A1 * A2[Formula: see text]. In addition, we also introduced a theoretical model to calculate the TOT as a function of energy loss in plastic scintillators.

CONCLUSIONS

A relationship between TOT and energy loss by photons interacting inside the plastic scintillators used in J-PET scanner is established for a deposited energy range of 100-1000 keV.

State of the art in total body PET

The idea of a very sensitive positron emission tomography (PET) system covering a large portion of the body of a patient already dates back to the early 1990s. In the period 2000-2010, only some prototypes with long axial field of view (FOV) have been built, which never resulted in systems used for clinical research. One of the reasons was the limitations in the available detector technology, which did not yet have sufficient energy resolution, timing resolution or countrate capabilities for fully exploiting the benefits of a long axial FOV design. PET was also not yet as widespread as it is today: the growth in oncology, which has become the major application of PET, appeared only after the introduction of PET-CT (early 2000).The detector technology used in most clinical PET systems today has a combination of good energy and timing resolution with higher countrate capabilities and has now been used since more than a decade to build time-of-flight (TOF) PET systems with fully 3D acquisitions. Based on this technology, one can construct total body PET systems and the remaining challenges (data handling, fast image reconstruction, detector cooling) are mostly related to engineering. The direct benefits of long axial FOV systems are mostly related to the higher sensitivity. For single organ imaging, the gain is close to the point source sensitivity which increases linearly with the axial length until it is limited by solid angle and attenuation of the body. The gains for single organ (compared to a fully 3D PET 20-cm axial FOV) are limited to a factor 3-4. But for long objects (like body scans), it increases quadratically with scanner length and factors of 10-40 × higher sensitivity are predicted for the long axial FOV scanner. This application of PET has seen a major increase (mostly in oncology) during the last 2 decades and is now the main type of study in a PET centre. As the technology is available and the full body concept also seems to match with existing applications, the old concept of a total body PET scanner is seeing a clear revival. Several research groups are working on this concept and after showing the potential via extensive simulations; construction of these systems has started about 2 years ago. In the first phase, two PET systems with long axial FOV suitable for large animal imaging were constructed to explore the potential in more experimental settings. Recently, the first completed total body PET systems for human use, a 70-cm-long system, called PennPET Explorer, and a 2-m-long system, called uExplorer, have become reality and first clinical studies have been shown. These results illustrate the large potential of this concept with regard to low-dose imaging, faster scanning, whole-body dynamic imaging and follow-up of tracers over longer periods. This large range of possible technical improvements seems to have the potential to change the current clinical routine and to expand the number of clinical applications of molecular imaging. The J-PET prototype is a prototype system with a long axial FOV built from axially arranged plastic scintillator strips.This paper gives an overview of the recent technical developments with regard to PET scanners with a long axial FOV covering at least the majority of the body (so called total body PET systems). After explaining the benefits and challenges of total body PET systems, the different total body PET system designs proposed for large animal and clinical imaging are described in detail. The axial length is one of the major factors determining the total cost of the system, but there are also other options in detector technology, design and processing for reducing the cost these systems. The limitations and advantages of different designs for research and clinical use are discussed taking into account potential applications and the increased cost of these systems.

Feasibility study of the positronium imaging with the J-PET tomograph

A detection system of the conventional PET tomograph is set-up to record data from [Formula: see text] annihilation into two photons with energy of 511 keV, and it gives information on the density distribution of a radiopharmaceutical in the body of the object. In this paper we explore the possibility of performing the three gamma photons imaging based on ortho-positronium annihilation, as well as the possibility of positronium mean lifetime imaging with the J-PET tomograph constructed from plastic scintillators. For this purposes simulations of the ortho-positronium formation and its annihilation into three photons were performed taking into account distributions of photons' momenta as predicted by the theory of quantum electrodynamics and the response of the J-PET tomograph. In order to test the proposed ortho-positronium lifetime image reconstruction method, we concentrate on the decay of the ortho-positronium into three photons and applications of radiopharmaceuticals labeled with isotopes emitting a prompt gamma. The proposed method of imaging is based on the determination of hit-times and hit-positions of registered photons which enables the reconstruction of the time and position of the annihilation point as well as the lifetime of the ortho-positronium on an event-by-event basis. We have simulated the production of the positronium in point-like sources and in a cylindrical phantom composed of a set of different materials in which the ortho-positronium lifetime varied from 2.0 ns to 3.0 ns, as expected for ortho-positronium created in the human body. The presented reconstruction method for total-body J-PET like detector allows to achieve a mean lifetime resolution of  ∼40 ps. Recent positron annihilation lifetime spectroscopy measurements of cancerous and healthy uterine tissues show that this sensitivity may allow to study the morphological changes in cell structures.

The tests of CP and CPT symmetry using the J-PET detector

Symmetries under the parity transformation (P), charge-conjugation (C) and time reversal (T) are of fundamental importance in nuclear and elementary particle physics. Studies of the observables violating the combined CP symmetry constitute precise tests of the Standard Model. However, CP violation was observed to date only for systems involving quarks, raising the importance of searches its manifestations e.g. in purely leptonic systems. The 3γ decay of spin-aligned ortho-positronium atoms (o-Ps) can be used to test CP invariance in such a purely leptonic system. The Jagiellonian Positron Emission Tomograph (J-PET) detection system enables experimental tests of CP and CPT through measurement of the expectation values of angular correlation operators odd under these transformations and constructed from (i) spin vector of the ortho-positronium atom, (ii) co-planar momentum vectors of photons originating from the decay of the positronium atom, and (iii) linear polarization direction of annihilation photons. Precise experimental symmetry tests with J-PET are possible thanks to a dedicated reconstruction technique of 3γ ortho-positronium decays and a positronium production chamber including a highly porous aerogel target, whose setup allows for determining the orthopositronium spin polarization without the use of an external magnetic field.

Time Over Threshold as a measure of energy response of plastic scintillators used in the J-PET detector

The Jagiellonian Positron Emission Tomograph (J-PET) is a multipurpose detector being developed to provide an economical alternative of commercially available PETs as well as to perform the tests on the discrete symmetries and entanglement. It is composed of 192 plastic scintillators axially arranged in three cylindrical layers. In the framework of J-PET detector, Time-Over-Threshold (TOT) approach is adopted for the signal readouts in order to utilize the excellent time resolution of the plastic scintillators. In this paper, we present a method elaborated for establishing a relation between TOT and energy loss.

Feasibility studies of the polarization of photons beyond the optical wavelength regime with the J-PET detector

J-PET is a detector optimized for registration of photons from the electron-positron annihilation via plastic scintillators where photons interact predominantly via Compton scattering. Registration of both primary and scattered photons enables to determinate the linear polarization of the primary photon on the event by event basis with a certain probability. Here we present quantitative results on the feasibility of such polarization measurements of photons from the decay of positronium with the J-PET and explore the physical limitations for the resolution of the polarization determination of 511 keV photons via Compton scattering. For scattering angles of about 82∘ (where the best contrast for polarization measurement is theoretically predicted) we find that the single event resolution for the determination of the polarization is about 40∘ (predominantly due to properties of the Compton effect). However, for samples larger than ten thousand events the J-PET is capable of determining relative average polarization of these photons with the precision of about few degrees. The obtained results open new perspectives for studies of various physics phenomena such as quantum entanglement and tests of discrete symmetries in decays of positronium and extend the energy range of polarization measurements by five orders of magnitude beyond the optical wavelength regime.