- Open Access
- Total Downloads : 105
- Authors : Dhanya. F. E , Sandhya. S
- Paper ID : IJERTV5IS070407
- Volume & Issue : Volume 05, Issue 07 (July 2016)
- DOI : http://dx.doi.org/10.17577/IJERTV5IS070407
- Published (First Online): 01-08-2016
- ISSN (Online) : 2278-0181
- Publisher Name : IJERT
- License: This work is licensed under a Creative Commons Attribution 4.0 International License
Simulation of Positron Emission Tomography Tumor Images using MRI and CT
Dhanya.F.E Sandhya.S
PG Scholar, Dept of ECE Assistant Professor
Mohandas college of Engineering Dept of Electronics and Communication Thiruvananthapuram Thiruvananthapuram
Abstract:- MRI and CT are commonly used for disease detection. A new technique is introduced in this area Positron Emission Tomography (PET) image analysis. This is able to detect most of the diseases and most appropriate treatment can be provided to the patients. Tumor goes more deeply, PET image is more used in the diagnoses of tumors. Now, some simulation programs such as GATE (Geant4 Application for Tomographic Emission) exist, but too complicated and more time consuming. This work proposes a new method to simulate the PET image of the brain with Monte Carlo simulation in Matlab. . For the simulation, a high resolution, MRI and CT based, segmented image, is used as the original image. In MRI segmentation intensity non – uniformity (INU) artifact present.Hence, Adaptive Spatial Fuzzy Center means segmentation is used. It is based on fuzzy C – means that address both INU artifact and local spatial continuity.
Key word:- INU artifact,FCM,ASFCM, membership values, clusters, fuzziness factor
I. INTRODUCTION
A wide range of different image modalities for medical imaging are available at present, which provide view of internal organs of the human body such as brain, kidney, liver etc. Among these medical image modalities, Magnetic Resonance Imaging (MRI) and Computed Tomography (CT) imaging are commonly used for disease diagnosis. Another technique used in this area is Positron Emission Tomography (PET). Unlike earlier radiotracer techniques PET can quantitatively measure biochemical and physiological processes. PET scans are more accurate in detecting tumors that are smaller or less aggressive than other imaging techniques.
For clinical PET scanning, positron emission radioactive isotopes are injected into the human body. The following are the major drawbacks of this type of scanning:
-
radio active isotopes have side effects on patients and
-
PET scan is expensive. For these two reasons, PET simulators are needed for physical and clinical research. Here a new method is proposed to simulate the PET tumor image of the brain using analytical information from CT and MRI images with Monte Carlo Simulation in Matlab. For the simulation, MRI and CT based segmented image is used as the original image. The main draw back of the MRI scanning is that it cannot provide clear information of hard tissues such as bones and skull while CT scan provides this information. A combination of CT scan and MRI scan is used in this work.
In neuroimaging, Monte Carlo Simulation has been applied in several elds such as the evaluation of reconstruction methods ,statistical toolkits, quantication methods, registration algorithms, segmentation methods, etc. Monte Carlo is a numerical calculation method in which random variables play an important role. Random numbers are used to simulate the transport of particles from its emission site for either its detection in the detector or its absorption in the medium in PET simulation.
-
METHODOLOGY
For the simulation of PET image MRI and CT images are used as initial maps.from these attenuation and activity maps are created .MRI images are segmented to 1) extracted brain image, 2) outskin image, and 3) brain tissue segmented image . Brain tissue segmented image means separate grey matter , white matter and cerebrospinal uid(CSF). Thresholding methods is used for bone tissue extraction from the CT image.From the combination of the bone tissue and outskin image attenuation maps were created. Information from all images are requird for generate the activity maps.MC code was used to simulate PET image.
-
SEGMENTATION
MRI and CT data were used to generate the initial images. MRI studies were brain extracted and segmented to obtain the image set 1) extracted brain image, 2) outskin image, and 3) brain tissue segmented image (grey matter , white matter and cerebrospinal uid(CSF) separately). Bone tissue was extracted from the CT image by thresholding.
Segment the MRI image into extracted brain image and outskin image by linear thresholding. Brain tissue segmention is done by Adaptive Spatial Fuzzy center mean algorithm (FCM)[16]. Adaptive Spatial Fuzzy C-means clustering is a data clustering method .In which each data point belongs to a cluster to a degree specified by a membership value.
-
Adaptive Spatial Fuzzy C-Means Clustering
-
Membership matrix U randomly select according to the size of the image.
U=[Uij], 1 i c,1 j n
It satisfies the following condition:
c
SFCM Segmentation algorithm :
I Uij = 1
i=1
j = 1,2, , n
-
Read the MRI image
Where Uij- membership grade ,n- number of samples c- number of clusters
-
-
Then calculate the Cluster center
-
n Um.xi
-
Assign number of clusters
-
Assume the fuzziness factor
-
For i=1:max_iter
-
Distance between pixel and centroid is calculated
-
Membership values are calculated
ci = j=l ij
U
U
n m
j=l ij
Where m [1,] is the weighting exponent xi is the data point
-
Calculate distance between the pixel and center
d2 = lx c l2
End for
-
-
-
If distance less than threshold
-
Membership values are updated
-
Update cluster centroid
-
-
Do the Segmentation
ij i j
-
If the distance is less than the threshold 1e-5 then update the membership matrix
-
The new membership matrix is
Um = 1
m-l
m-l
k=1 )
k=1 )
ij c (dij 2
dik
-
Terminating condition: Aims at minimizing the
desired function
c n
]ASFCM = I I Um d2
-
-
EXPERIMENTAL RESULTS
MRI studies were brain extracted and segmented to obtain the image set 1) extracted brain image, 2) outskin image, and 3) brain tissue segmented image (grey matter , white matter and cerebrospinal uid(CSF) separately). Bone tissue was extracted from the CT image by thresholding.The attenuation maps
i=1
ij ij
j=1
were generated from a combination of the bone tissue and outskinimage.To generate the activity maps, information from all images was taken into account.
Fig 8.5 Outskin image from MRI
Fig 8.1 Input CT image
-
Using FCM b)Using ASFCM
Fig 8.6 Segmented MRI
Fig 8.2 Input MRI
-
Using FCM b)Using ASFCM
-
-
-
Fig 8.7 Simulated PET
CONCLUSION
Fig 8.3 Bone image from CT
Fig 8.4 Extracted brain image from MRI
In this work ,the efficiency of the PET simulation is demonstated by FCM and ASFCM segmentation.In this method first assign the membership marix and fuzziness factor.From these values update the membership value according to the distance.Hence brain tissue segmentation is obtained correctly.
REFERENCES
-
Berta Marti-Fuster, Oscar Esteban, Kris Thielemans, Xavier Setoain, Andres Santos, Domenec Ros, and Javier Pavia "Including Anatomical and Functional Information in MC Simulation of PET and SPECT Brain Studies". Brain-VISET:A Voxel-Based Iterative Method IEEE trans. On med. Imag.,vol33,no.10,oct 2014.
-
K. Thielemans, C. Tsoumpas, S. Mustafovic, T. Beisel, P. Aguiar,
N. Dikaios, and M.W. Jacobson, "STIR: Software for tomographic image reconstruction release 2," IEEE Trans. on Phys.Med. Biol., vol. 57, no. 4, pp. 867 – 883, 2012.
-
C.Catana, A.vander Kouwe, T.Benner, C.Michel, M.Hamm, M.Fenchel, B.Fischl, B.Rosen, M.Schmand, and A.Sorensen, "Toward implementing an MRI-based PET attenuation correction method for neurologic studies on the MR-PET brain prototype," IEEE J. Nucl. Med. 51, 1431-1438 May 2010.
-
Jun Ma and Malcolm Hudson "Block-Iterative Fisher Scoring Algorithms for Maximum Penalized Likelihood Image
Reconstruction in Emission Tomography" IEEE Trans on med imag, vol. 27, no. 8, aug 2008.
-
Fabiana Crepaldi and Alvaro R. De Pierro, Member, "Activity and Attenuation Reconstruction for Positron Emission Tomography Using Emission Data Only Via Maximum Likelihood and Iterative Data Renement " IEEE Trans on nuclear science, vol. 54, no. 1, feb 2007.
-
J. Cal-Gonzalez, J. L. Herraiz, S. Espana, M. Desco, J. J. Vaquero, J. M. Udias, "Positron Range Effects in High Resolution 3D PET Imaging", IEEE Trans. On Phys. Med. Biol., Vol. 29, pp. 1217-1230, May 2007
-
Kris Thielemans, member IEEE, Sanida Mustafovic, Charalampos Tsoumpas "STIR: Software for Tomographic Image Reconstruction Release 2" Nuclear Science Symposium Conference Record May 2006.
-
Steven Staelens, Daniel Strul, Giovanni Santin, Stefaan Vandenberghe, Michel Koole, Yves D Asseler, Ignace Lemahieu, Rik Van de Walle, "Monte Carlo Simulations of a Scintillation Camera Using GATE: Validation and Application and Modeling", IEEE Trans. on Phys. Med. Bio., vol. 48, pp. 3021-3042, Aug. 2003.
-
A. Reilhac, G. Batan, C. Michel, N. Costes, and A. C. Evans, "PET- SORTEO: a platform for simulating realistic PET studies." IEEE Nucl. Sci. Symp. and Med. Imag. Conf., Oct. 2004.
-
Karine Assie, Vincent Breton, Irene Buvat, Claude Comtat, Sebastien Jan, Magalie Krieguer, Delphine Lazaro, Christian Morel, Martin Rey, Giovanni Santin, Luc Simon, Steven Staelens, Daniel Strul, Jean-Marc Vieira, Rik Van de Walle, "Monte Carlo Simulation in PET and SPECT Instrumentation Using GATE", IEEE Trans on Nuclear Instruments and Methods in Physics Research, vol. 527, pp.180-189 Apr.2004.
-
D. strul, G. Santin, D. Lazaro, V. Breton, C. Morel, "GATE (Geant 4 Application for Tomographic Emission): a PET/SPECT General-Purpose Simulation Platform", IEEE Trans on Nuclear Physics ,vol. 125 ,pp.75-79,Aug 2003.
-
G. Santin, D. strul, D. Lazaro, L. Simon, M. Vieira Martins, V. Breton, C. Morel, "GATE: A Geant4-Based Simulation Platform for PET and SPECT Integrating Movement and Time Management", IEEE Transactions on Nuclear Science, VOL. 50, No. 5, October 2003.
-
V. Frouin, C. Comtat, A. Reilhac, and M.-C. Gregoire, "Correction of partial-volume effect for PET striatal imaging: Fast implementation and study of robustness," J. Nucl. Med., vol. 43, no. 12, pp. 1715-1726, Dec. 2002.
-
D. S. Lalush, E. C. Frey, and B. M. Tsui, "Fast maximum entropy approximation in SPECT using the RBI-MAP algorithm," IEEE Trans. on Med. Imag., vol. 19, pp. 286-294, Apr. 2000.
-
H. Zaidi, "Relevance of accurate Monte Carlo modeling in nuclear medicalimaging,"IEEE on Med.Phys.,vol.26,no.4,pp.574
– 608,Aug 1999
-
Cheng HD, Lui YM, Freimanis RI, "A novel approach to micro calcification detection using fuzzy logic technique", IEEE Trans Med Imaging , Vol. 17, pp. 442-450, Sep 1998
-
C. J. Thompson,J. Moreno-Cantu, Y. Picard, "PETSIM: Monte Carlo Simulation of all Sensitivity and Resolution Parameters of Cylindrical Positron Imaging Systems",IEEE on Phys. Med. Biol., Vol. 37, No. 3, pp.731-739 Aug 1992.
-
C. E. Floyd, R. J. Jaszczak, C. C. Harris, R. E. Coleman, "Energy and Spatial Distribution of Multiple Order Compton Scatter in SPECT: a Monte Carlo Investigation", IEEE Trans. on Phys. Med. Biol., Vol. 29, pp. 1217-1230, May 1984.