# Scatter rate estimation for LM-TOF-MLEM **URL:** \ **Category:** General Discussions\ **Created:** [August 17, 2021, 10:52pm UTC](https://castor-project.discourse.group/t/scatter-rate-estimation-for-lm-tof-mlem/229 "2021-08-17T22:52:29Z")\ **Posts on this page:** 5\ **Page:** 1
**Author:** ![jbaran](https://avatars.discourse-cdn.com/v4/letter/j/90db22/32.png) [@jbaran](https://castor-project.discourse.group/u/jbaran)\ **Post date:** [August 17, 2021, 10:52pm UTC](https://castor-project.discourse.group/t/scatter-rate-estimation-for-lm-tof-mlem/229/1 "2021-08-17T22:52:29Z")
Dear Castor Developers and Users, I am trying to incorporate the information about the scatter correction in LM-TOF-MLEM reconstruction. I have calculated the scatter fraction for each TOF bin. The values of scatter rate are in range 1-2. I have applied them for the reconstruction, however the reconstructed images seem to be not corrected for scatter. Please find the attached screenshot. On the top panel, the reconstructed image with true+scatter coincidences is presented and on the bottom the same coincidences were used but with scatter correction applied. On the right corresponding line profiles are presented (ots - no scatter correction; corrected1 - scatter correction applied) FYI: The total number of coincidences used for the reconstruction is 1 order of magnitude lower than the number of possible LORs. Do you have any idea what kind of values should be applied for the scatter rate? Did you have any similar case and found the solution how the scatter rate should be calculated? All the best, Jakub ![Screenshot from 2021-08-18 00-36-18.jpg](https://global.discourse-cdn.com/free1/uploads/castor_project/original/1X/400efd47c2e11ea58a9f237b38c80594eff9bf7b.jpeg) ---
**Author:** ![mfilipovic](https://avatars.discourse-cdn.com/v4/letter/m/858c86/32.png) [@mfilipovic](https://castor-project.discourse.group/u/mfilipovic)\ **Post date:** [August 18, 2021, 2:54pm UTC](https://castor-project.discourse.group/t/scatter-rate-estimation-for-lm-tof-mlem/229/2 "2021-08-18T14:54:04Z")
Hi Jacub, Have you divided the estimated TOF bin scatter rates with the spatial width of the TOF bin, directly in the list-mode datafile, as explained in the documentation (general and TOF-specific) ? Marina ---
**Author:** ![jbaran](https://avatars.discourse-cdn.com/v4/letter/j/90db22/32.png) [@jbaran](https://castor-project.discourse.group/u/jbaran)\ **Post date:** [August 25, 2021, 6:52am UTC](https://castor-project.discourse.group/t/scatter-rate-estimation-for-lm-tof-mlem/229/3 "2021-08-25T06:52:29Z")
Hi Marina, many thanks for your answer. could you please clarify what do you mean to "divide the estimated TOF bin scatter rates with the spatial width of the TOF bin, directly in the list-mode datafile". Once I estimate the scatter fraction I presume I should divide it by the TOF bin (e.g. 200ps)? However, I do not specify the TOF bin size in my header. Unless the bin size is then defined as a TOF resolution? Example of the header which I am using: Data filename: data\_df.Cdf Number of events: 298175691 Data mode: list-mode Data type: PET Start time (s): 0 Duration (s): 500 Scanner name: scanner1 Calibration factor: 1 Isotope: unknown TOF information flag: 1 TOF resolution (ps): 200 List TOF measurement range (ps): 3000 Scatter correction flag: 1 What is the way to proceed in that case? All the best, Jakub ---
**Author:** ![mfilipovic](https://avatars.discourse-cdn.com/v4/letter/m/858c86/32.png) [@mfilipovic](https://castor-project.discourse.group/u/mfilipovic)\ **Post date:** [August 26, 2021, 7:08am UTC](https://castor-project.discourse.group/t/scatter-rate-estimation-for-lm-tof-mlem/229/4 "2021-08-26T07:08:28Z")
Hi Jakub, You should divide your scatter rate with the size of the TOF bin converted to mm. That’s if you’re using the “continuous” LM-TOF recon, which is run if you don’t specify the size of the TOF bin in the header. If you specify the size of the TOF bin in the header, then you don’t have to change your scatter rates per bin, that’s easier, I should have mentioned this option first. The TOF resolution field is only for TOF resolution. The documentation is maybe not so clear about these options? Let me know if this helps. Cheers Marina ---
**Author:** ![jbaran](https://avatars.discourse-cdn.com/v4/letter/j/90db22/32.png) [@jbaran](https://castor-project.discourse.group/u/jbaran)\ **Post date:** [September 20, 2021, 11:18am UTC](https://castor-project.discourse.group/t/scatter-rate-estimation-for-lm-tof-mlem/229/5 "2021-09-20T11:18:22Z")
Hi Marina, many thanks for your comments. Unfortunately, the solution you provided didn't change the results much. The greater background in the middle of the imaged object, characteristic for the lack of scatter correction is still visible. We are suspecting that probably we need greater coincidence statistics to sort it out. I will let you know if we found anything interesting. All the best, Jakub