Abstract:
A tomography system based on Cerenkov tomography, comprising: a detector of Cerenkov fluorescence for acquiring optical plane images; a structural imaging system for acquiring three-dimensional structural images; 
     a bed device for supporting an object to be imaged; a computer for forming an optical image, a structural image and a CLT image. The invention adopts the SP3 model and the semi threshold iterator to implement the global reconstruction of the CLT, and obtains the three-dimensional tomography to image of the distribution of the radiopharmaceutical and the molecular probe in vivo within a short time. Since ordinary CCD camera is used, the cost of the imaging system has been sharply reduced for the equipment&#39;s construct and maintenance compared with PET/SPECT or γ camera. Therefore the present invention expands the options of the molecular probe, and application of the is medicine Imaging.

Description:
FIELD OF THE INVENTION 
       [0001]    The present invention generally relates to an imaging system. More particularly, the present invention relates to method and system for Cerenkov Luminescence Tomography. 
       DESCRIPTION OF THE RELATED ART 
       [0002]    The propositions of Cerenkov Luminescence Imaging (CLI) establish the foundation of clinical application of molecular imaging. In 2010, A. Ruggiero (J. Nucl. Med. 2010,51:1123-1130) proved that the high-speed movement of high-energy α or β particles cause Cerenkov Emission, which can result in visible photons. Those photons can penetrate through the object to be imaged and get to the surface. With those photons, the optical imaging system can provide a vivo nuclear medicine image. It means the optical imaging becomes a new mode besides PET (Positron Emission Tomography) and SPECT (Single Photon Emission Computed Tomography). 
         [0003]    While CLI can obtain the two-dimensional distribution of the nuclide without the depth information of the radiation source, Spinelli ( Phys. Med. Biol.  2010,55:483-495) adopt multispectral technology to estimate the depth of the radiation source, but the essence of the method is still single perspective two-dimensional image. Li ( Opt. Lett.  2010,35:1109-1111) provides a new method named CLT (Cerenkov Luminescence Tomography) and achieves three-dimensional image based on the assumption of uniformity model. Both of the methods listed above admit the assumption of Diffusion Approximation theory. But actually, the spectrum of Cerenkov Effect is not limited in 300-700 nm, the distribution of the energy is inversely proportional to the square of the wavelength. A complex biological tissue presents characters as below when penetrated by a blue dominated light: the ratio of the optical scattering coefficient and the absorption coefficient is usually less than 10. As a result, the DA model is no long fit to describe the dissemination of the light in body. To improve the quality of the image got by CLT, we have already built and proved the RTE (Radiation Transport Equation) and SP3 (The Third-order Simplified Spherical Harmonics Approximation) model to describe the phenomenon of the Cerenkov photons&#39; dissemination in the body. However, distribution of the radiopharmaceutical in body is partial centralized and global sparse. However it is difficult to have an inverse model for the multiple-peak and sparse global method based on the CLT 
       SUMMARY OF THE INVENTION 
       [0004]    The subject matter of the present invention is to provide a tomography system and a method based on Cerenkov Effect. 
         [0005]    According to one aspect of the present invention, a tomography method based on the Cerenkov Tomography Imaging comprising: 
         [0006]    Acquiring an optical plane image and a three-dimensional structural image; 
         [0007]    Obtaining a space registration of the optical plane image collected by different perspective and the three-dimensional structural image based on fixed marking points, and constructing a geometrical mapping table of the optical and a tomography image; 
         [0008]    Converting every pixel&#39;s gray value of the optical plane image into an illumination intensity of three-dimensional object based on the geometrical mapping table, and obtaining a distribution B of the illumination intensity on a surface of the object; 
         [0009]    Segmenting the three-dimensional tomography image in order to obtain each organ, and the organs being assigned by different pixel values; 
         [0010]    The segmented body data being discratized by a finite element method, and outputting a volume gridding; 
         [0011]    Based on the SP3 forward model, constructing the linear mapping table between the distribution of the surface light intensity B and the unknown distribution of the isotope in vivo S, B=AS, where A stands for the system parameter matrix of the forward model, B stands for a m-dimensional column vector, S stands for an n-dimensional column vector, m is a signless integrals less than n; 
         [0012]    Based on the theory of the L 1/2  regularization, objective function F(S) being calculated for reconstructing tomography image using iterator as follows: F(S)=min{∥B−AS∥ 2 +p∥s∥ 1/2 }, where p stands for the parameter of the regularization, p=(1−ε)/∥A∥ 2 ; 
         [0013]    In the K+1 St  iterate, the iterator update the distribution of the isotope in vivo. S k+1 =H pq  (T k (S k )), T k (S k )=S k +q k A T (B−AS k ), where k is a signless integral, q k =4|T k | 3/2 ∥A∥ 2 /3, H=(h(S 1 ),h(S 2 ) . . . h(S n )) T , 
         [0000]    
       
         
           
             
               
                 h 
                  
                 
                   ( 
                   S 
                   ) 
                 
               
               = 
               
                 2 
                  
                 
                   
                     S 
                      
                     
                       ( 
                       
                         1 
                         + 
                         
                           cos 
                            
                           
                             
                               
                                 2 
                                  
                                 π 
                               
                               - 
                               
                                 2 
                                  
                                 φ 
                               
                             
                             3 
                           
                         
                       
                       ) 
                     
                   
                   / 
                   3 
                 
               
             
             , 
           
         
       
     
         [0000]    when S&gt;η(pq) 3/2 , h(S)=0 and S≦η(pq) 3/2 , 
         [0000]    
       
         
           
             
               φ 
               = 
               
                 ar 
                  
                 
                     
                 
                  
                 
                   cos 
                   ( 
                   
                     0.125 
                      
                     
                         
                     
                      
                     
                       
                         pq 
                          
                         
                           ( 
                           
                             
                                
                               S 
                                
                             
                             3 
                           
                           ) 
                         
                       
                       
                         
                           - 
                           3 
                         
                         / 
                         2 
                       
                     
                   
                   ) 
                 
               
             
             , 
           
         
       
     
         [0000]    η stands for a real constant; 
         [0014]    If the S got by the kst makes the objective function F k &lt;V, stop the iterator, where V is a threshold constant. 
         [0015]    According to other aspect of the present invention, a tomography system based on Cerenkov tomography, comprising: 
         [0016]    a detector of Cerenkov fluorescence for acquiring optical plane images; 
         [0017]    a structural imaging system for acquiring three-dimensional structural images; 
         [0018]    a bed device for supporting an object to be imaged; 
         [0019]    a computer for forming an optical image, a structural image and a CLT image. 
         [0020]    The invention adopts the SP3 model and the semi threshold iterator to implement the global reconstruction of the CLT, and obtains the three-dimensional tomography image of the distribution of the radiopharmaceutical and the molecular probe in vivo within a short time. Since ordinary CCD camera is used, the cost of the imaging system has been sharply reduced for the equipment&#39;s construct and maintenance compared with PET/SPECT or γ camera. Therefore the present invention expands the options of the molecular probe, and application of the medicine Imaging. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0021]    to  FIG. 1  shows the three-dimensional structure of the system; 
           [0022]      FIG. 2  shows the constitution of the entire invention; 
           [0023]      FIG. 3  shows the flow chart of the method to reconstruct the Cerenkov tomography image; 
           [0024]      FIG. 4  shows the schematic diagram of the process of the reconstructing the Cerenkov tomography image; 
           [0025]      FIG. 5  shows the registering relational graph of the optical and the structural image; 
           [0026]      FIG. 6  shows the distribution of the vexels which have the most isotope according to the algorithm of the reconstruction. 
           [0027]      FIG. 7  shows the Cerenkov tomography and the three orthogonal sectional drawing which go through the maximum vexel. 
       
    
    
     DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 
       [0028]    As shown in  FIG. 1 , a Tomography system of the present invention comprises a detector ( 101 ) for Cerenkov fluorescence, a computer ( 109 ), a camera obscura ( 110 ), a structural imaging system ( 102 , 103 ), a bed device ( 104 , 105 , 106 , 107 , 108 ) which supports a object to be imaged (i.e patient, animal or others), the bed is called a small animal bed. The detector for Cerenkov fluorescence ( 101 ) is the key point of the present invention, comprising a CCD camera that works at −110 degree, a band-pass filter bank with a bandwidth ranged from 10 nm to 20 nm and a centre wavelength ranged from 500 nm to 650 nm. The band-pass filter bank is installed on the axis of the stepper motor. There are six optical holes to collect photons of different wave bands. One of them does not have an optical filter plate. The small animal stents ( 104 ), which can be detached from the imaging system, has 24 marking point with diameter of 2 mm, which can be recognized by the optical camera and the micro-CT. The bottom of the stents is fixed on the 360-degree rotating platform ( 105 ), which has three irrelevant stepper motor. The small animal can change position by the rotating platform and the orthogonal translation guide ( 106 , 107 ). The structural imaging device is a high-resolution micro-CT system, which is orthogonal to the Cerenkov photon detector on the air-cushion optical platform. The front side of the camera obscura ( 110 ) is made of lead crystals, the other sides are made of lead plate. The inner face of the camera obscura is cover by blank dope. There is a 1 cm diameter hole on the bottom of the camera obscura, which can let the cable, which connect the controller and the computer, and the pipe of anesthesia gas get through. 
         [0029]    The  FIG. 2  shows the frame diagram of the tomography system based on Cerenkov Luminescence. The Cerenkov photon detector ( 202 ) comprises a rotating band-pass filter bank ( 205 ), a CCD device to convert the optical signal into the electrical signal ( 206 ), and an optical imaging controller ( 207 ). The optical imaging controller can adjust the bandwidth and the centre wavelength of the band-pass filter bank, the working temperature of the CCD camera, the frequency of the data acquisition, the time of the exposure, and the openness of the aperture by the command of the computer. The small animal bed ( 203 ) comprise a small animal stents ( 208 ) for fixing the object to be imaged, a 360-degree rotating platform ( 209 ), a translation guide ( 210 , 211 ), and a Controller ( 212 ). The controller for the bed device ( 212 ) can change the position of the object to be imaged based on the command of the computer. The structural imaging system micro-CT ( 204 ) comprises a X-ray source ( 213 ), a X-ray detector ( 214 ), and a micro-CT controller ( 215 ). The micro-CT controller can adjust the parameters of the micro-CT system based on the command of the computer. The optical imaging controller ( 207 ) delivers the image acquired by the CCD camera to the computer, and obtains the optical plane image ( 216 ). The micro-CT controller ( 215 ) delivers the three-dimensional data acquired by the X-ray detector to the computer ( 219 ), and obtains the three-dimensional structural image ( 218 ). The computer calculates the Cerenkov tomography image ( 217 ) according to the optical plane image, the three-dimensional structural image, and the reconstruction algorithm. Therefore, the computer outside the camera obscura not only control the whole system, but also process and output the optical plane image ( 216 ), three-dimensional structural image ( 218 ), and Cerenkov tomography image ( 217 ). 
         [0030]    As shown in the  FIG. 3 , a computer ( 301 ) calculates the Cerenkov is tomography image ( 312 ) according to an optical plane image ( 302 ), a three-dimensional structural image ( 303 ), and a reconstruction algorithm (the arrows in  FIG. 3  show the message passing direction.) 
         [0000]    Step  304 : obtaining a space registration image of the optical plane image ( 302 ) and the three-dimensional structural image ( 303 ) based on a fixed marking point, And constructing a geometrical mapping table of the optical and the CT tomography image.
 
Step  305 : Converting every pixel&#39;s gray value of the optical plane image into illumination intensity of the three-dimensional object based on the geometrical mapping table acquired by step  304 , and obtaining the distribution of the illumination intensity on the surface of the object (B).
 
Step  306 : Segmenting the three-dimensional tomography image acquired by the micro-CT, in order to obtain each organ, and the organs are assigned by different pixel values. For instance, after the image is segmented, we obtain the images of the heart, the kidney, the liver, the skeleton, the muscle, and the bladder, etc.
 
         [0031]    Step  307 : obtaining the volume gridding of the body is dispersed by using the finite element method. For example, a tetrahedron gridding is used. 
         [0000]    Step  308 : Based on the SP3 forward model, constructing a linear mapping table between the distribution of the surface light intensity B and the unknown distribution of the isotope in vivo S , i.e B=AS, where A stands for the system parameter matrix of the forward model, B stands for a m-dimensional column vector, S stands for an n-dimensional column vector. m and n are signless integral, and m&lt;n.
 
Step  309 : Based on the theory of the L 1/2  regularization, objective function F(S) is calculated for reconstructing tomography image using iterator as follows: F(S)=min {∥B−AS∥ 2 +p∥s∥ 1/2 }, where p stands for the parameter of the regularization, p=(1−ε)/∥A∥ 2 .
 
Step  310 : In the K+1 st  iterate, the iterator update the distribution of the isotope in vivo. S k+1 =H pq (T k (S k )), T k (S k )=S k +q k A T  (B−AS k ), where k is a signless integral, q k =4|T k | 3/2 ∥A∥ 2 /3, H=(h(S 1 ), h(S 2 ) . . . h(S n )) T ,
 
         [0000]    
       
         
           
             
               
                 h 
                  
                 
                   ( 
                   S 
                   ) 
                 
               
               = 
               
                 2 
                  
                 
                   
                     S 
                      
                     
                       ( 
                       
                         1 
                         + 
                         
                           cos 
                            
                           
                             
                               
                                 2 
                                  
                                 π 
                               
                               - 
                               
                                 2 
                                  
                                 φ 
                               
                             
                             3 
                           
                         
                       
                       ) 
                     
                   
                   / 
                   3 
                 
               
             
             , 
           
         
       
     
         [0000]    when S&gt;η(pq) 3/2 , h(S)=0 and S≦η(pq) 3/2 , 
         [0000]    
       
         
           
             
               φ 
               = 
               
                 ar 
                  
                 
                     
                 
                  
                 
                   cos 
                   ( 
                   
                     0.125 
                      
                     
                         
                     
                      
                     
                       
                         pq 
                          
                         
                           ( 
                           
                             
                                
                               S 
                                
                             
                             3 
                           
                           ) 
                         
                       
                       
                         
                           - 
                           3 
                         
                         / 
                         2 
                       
                     
                   
                   ) 
                 
               
             
             , 
           
         
       
     
         [0000]    η stands for a real constant.
 
Step  311 : If the S got by the kst makes the objective function F k &lt;V stop the iterator, where V is a threshold constant. Output the final reconstruction image, namely the Cerenkov tomography image ( 312 ).
 
         [0032]    As shown in the  FIG. 4 , the present invention utilizes 2 fluorine-2 deoxidation-D glucose ( 18 F-FDG) to implement tomography image as follows: anaesthetize a femina nu/nu nude mouse, which has a weight of 23 g, with the isoflurane. Inject 0.2 ml of small animal CT contrast medium Fenestra LC, and 1.11 MBp  18 F-FDG into the mouse, and set the parameter of the micro-CT system. Scanning voltage of the X-ray emission source is 50V, the power is 50 W, the integral time of the detector is 0.467 s, the rotational speed of the platform is 1.0°/s, the size of the projected image is 1120×2344, the imaging time is 3.0 s, and the number of the projection is 360. The CT scanning is conducted 30 min after the injection is completed, and acquire the three-dimensional structural image. After setting the parameter of the Cerenkov optical imaging system, CCD camera&#39;s time of exposure is set to be 3 min, the F number of the aperture is set to be 2.8, the focal length is 55 mm, the distance between the object to be imaged and the lens is 15 cm, the binning number is 2. Conduct the optical scanning without an optical filter 60 min after injection, acquire the optical image every 90° rotation. After the registering and fusion of the CT tomography image and the optical image, construct the distribution of the illumination intensity on the surface of the object (B), in which, there are 1934 nodes, 3899 dough sheet. After the segment of the volume data, obtain the kidney, the skeleton, the bladder, and the muscle. The geometric centre of the bladder is 34.7 mm, 14.5 mm, and 5.0 mm. Discretize the surface gridding using the finite element method, acquire the tetrahedron gridding with 3952 nodes, 25578 borders, 42471 triangle patch, and 20844 tetrahedrons. After acquiring the distribution of the illumination intensity on the surface and the volume gridding, set the weight of mixture spectrum as the optical parameters of different biological tissues. According to the reconstruction method, reconstruct the CLT image, get the distribution of the medical isotopes in vivo. 
         [0033]    The  FIG. 5  shows the distance errors of the registering of the micro-CT image and the optical image based on the fixed marking point. 
         [0034]    The  FIG. 6  shows the result of the image segment, while  601  points to the kidney,  602  points to skeleton,  603  points to bladder,  604  points to the vexel that has the most illumination intensity of the fluorescence, 605 points to the muscle. The geometric centre of  604  is (34.4 mm, 13.2 mm, 4.7 mm), in the bladder. 
         [0035]    The  FIG. 7  shows the distribution of the  18 F-FDG in vivo in the CLT reconstruction image. While  306  is the sectional drawing of X=34.4 mm,  307  is the sectional drawing of Y=13.2 mm,  308  is the sectional drawing of Z=4.7 mm,  309  shows the position of the three sectional drawing listed above in the CLT reconstruction image. The parameter of the CLT reconstruction is η=125. In the sectional drawing of the CLT reconstruction image, the illumination intensity is normalize into [0,1]. The reconstruction time of the whole image is less than 0.65s. Based on the assumption of FDG is uniform in bladder, the reconstruction errors is 1.4 mm. 
         [0036]    The present invention can achieve the three-dimensional imaging of radiopharmaceutical in vivo, provide a low cost molecular imaging method for clinical application. 
         [0037]    The foregoing exemplary embodiments and advantages are merely exemplary and are not to be construed as limiting the present invention. The present teaching can be readily applied to other type of apparatus. Also the description of the exemplary embodiments of the present invention is intended to be illustrative, and alternatives, modifications, and variations will be apparent to those skilled in the art.