UCL / PETPVC

Partial Volume Correction in PET
Apache License 2.0
51 stars 14 forks source link

PETPVC on PET reconstructed with FBP #13

Closed Bahri-M-A closed 7 years ago

Bahri-M-A commented 7 years ago

Dear Thomas,

We have recently tried PETpvc tool on FDG PET images reconstructed with FBP algorithm, but unfortunately we have got a white image for the MG, MTC, IY, and LABBE methods. Would you have an explanation for this kind of issue? Do you think it might due to the negative values of PET images?

Best Regards, Mohamed

bathomas commented 7 years ago

Dear Mohamed,

The negative values are not necessarily a problem. When you say 'white image', do you mean an image of all zeros or something else. Please could you explain what mask image you used? 3D or 4D? What values do the mask regions have?

Kind regards,

Ben

Bahri-M-A commented 7 years ago

Many thanks for your answer. The obtained images were zeros (see attached file "image number 4"). The used mask is a 4D image (WM + GM + CSF). These masks (attached file "1,2,3") were obtained using the spm segmentation of the structural MRI image (values between 0 and 1). results Thank you, Mohamed

bathomas commented 7 years ago

Hi,

The brain looks close to the edge of the image in the transaxial plane. I would recommend adding the background as a region too 1-(GM+WM+CSF). Then try applying MTC or IY. What happens then?

Bahri-M-A commented 7 years ago

Dear Thomas,

With the background as a forth region calculated as (GM+WM+CSF) - 1 the IY method does not work (see the following error: % --------------------------------- ➜ pvc petpvc -i rrnsCOENEN_AverageTW_PET.nii -m Mask_tissues_4D.nii -o PETpvc_IY.nii --pvc IY -x 6.0 -y 6.0 -z 6.0 Performing iterative Yang... sum in region 1 = 617462 sum in region 2 = 461294 sum in region 3 = 426752

[Error] failure applying iY on: rrnsCOENEN_AverageTW_PET.nii

itk::ExceptionObject (0x2fc5640) Location: "void petpvc::FuzzyCorrectionFilter::GenerateData() [with TImage = itk::Image<float, 4u>]" File: /home/bahri/Documents/Codes/PETPVC/lib/petpvcFuzzyCorrectionFilter.txx Line: 126 Description: itk::ERROR: FuzzyCorrFilter(0x2f9e5d0): Region 4contains negative voxels in mask. Remove this region. % ------------------------------------------------------------------------------------------------------------- however, the MTC method gave an image with very small values (see attached file). results_1 I also tried with the background calculated as 1 - (MG+WM+CSF), the IY method seems to give a correct results but MTC provide the same image (very low values : see attached file). results_2

Many Thanks,

Mohamed

bathomas commented 7 years ago

My apologies. I meant 1-(GM+WM+CSF). could you try that?

Ben

Get Outlook for Androidhttps://aka.ms/ghei36


From: Bahri-M-A notifications@github.com Sent: Friday, May 26, 2017 10:43:49 AM To: UCL/PETPVC Cc: Thomas, Ben; Comment Subject: Re: [UCL/PETPVC] PETPVC on PET reconstructed with FBP (#13)

Dear Thomas,

With the background as a forth region calculated as (GM+WM+CSF) - 1 the IY method does not work (see the following error: % --------------------------------- ➜ pvc petpvc -i rrnsCOENEN_AverageTW_PET.nii -m Mask_tissues_4D.nii -o PETpvc_IY.nii --pvc IY -x 6.0 -y 6.0 -z 6.0 Performing iterative Yang... sum in region 1 = 617462 sum in region 2 = 461294 sum in region 3 = 426752

[Error] failure applying iY on: rrnsCOENEN_AverageTW_PET.nii

itk::ExceptionObject (0x2fc5640) Location: "void petpvc::FuzzyCorrectionFilter::GenerateData() [with TImage = itk::Image<float, 4u>]" File: /home/bahri/Documents/Codes/PETPVC/lib/petpvcFuzzyCorrectionFilter.txx Line: 126 Description: itk::ERROR: FuzzyCorrFilter(0x2f9e5d0): Region 4contains negative voxels in mask. Remove this region. % ------------------------------------------------------------------------------------------------------------- however, the MTC method gave an image with very small values (see attached file). [results_1]https://cloud.githubusercontent.com/assets/3977392/26487250/034fed18-41ff-11e7-89c5-cf6d3c356b2d.png I also tried with the background calculated as 1 - (MG+WM+CSF), the IY method seems to give a correct results but MTC provide the same image (very low values : see attached file). [results_2]https://cloud.githubusercontent.com/assets/3977392/26487413/b6a0ac4a-41ff-11e7-934e-577a040594f0.png

Many Thanks,

Mohamed

— You are receiving this because you commented. Reply to this email directly, view it on GitHubhttps://github.com/UCL/PETPVC/issues/13#issuecomment-304226857, or mute the threadhttps://github.com/notifications/unsubscribe-auth/AGX4acwwfIUhzF7jCZ7PLNtaQl4xvxvZks5r9pDFgaJpZM4NlV8k.

Bahri-M-A commented 7 years ago

Dear Thomas, The second figure of my last email shows the results of the results of the IY and MTC methods as sub-figures 5 and 6, respectively. These results were obtained using the 1-(GM+WM+CSF) as background mask. As you can see, the results of the MTC method is not good. Best regards, Mohamed

bathomas commented 7 years ago

Dear Mohamed,

Would you be willing to share the nifti files via Dropbox or similar? If so, my email address is b.a.thomas [at] ucl.ac.uk. I can then take a look. I assume you are using a PSF of 6mm.

Kind regards,

Ben