in the image restoration framework , the regularization term is crucial to incorporate knowledge concerning a priori acceptable solutions and to constrain the solution of this ill - posed inverse problem . in this paper , we are concerned with image restoration issues using information from a different modality ( intermodality ) for this fusion - based regularization term .
more precisely , we herein propose an original extension of the nas - rif restoration technique allowing the nas - rif inverse filtering to be constrained from an intermodality registration , that is , a registration between anatomical and the functional medical images to be restored ( from the anatomical structure of the same patient ) . in our application ,
anatomical information is extracted from a high - resolution anatomical procedure such as magnetic resonance imaging ( mri ) or computed tomography ( ct ) and this high spatial resolution modality is exploited to improve the contrast of functional spect images .           
the contrast and signal - to - noise ratio displayed by brain single photon emission computed tomography ( spect ) images is rather limited , when compared with that from anatomical techniques ( mri , ct scanning ) .
for instance , it is not easy to differentiate low tracer uptake due to a functional deficit , where brain tissue still is anatomically intact , from low uptake generated by focal atrophy , where tissue is lost and replaced by cerebrospinal fluid ( csf ) .            up to now , several methods have been proposed to improve directly or indirectly the spatial resolution of spect images .
these methods can be split into two major classes : those using restoration techniques during the reconstruction process from projections , and those where the restoration is performed on the already reconstructed images . in this paper
, we are describing a posttomographic reconstruction process , an approach which has the advantage of being essentially independent from the physical features of the scanner .            in this category , we can cite   where rajabi   et al . compared four widely used filters ( i.e. , hanning , butterworth , metz and wiener ) for myocardial tc - sestamibi spect perfusion studies . in   a nonnegativity   and support constraints recursive inverse filtering ( nas - rif ) algorithm proposed by kundur and hatzinakos , was extended to the 3d spect imaging restoration context .
the nas - rif blind deconvolution technique is relevant to any situations in which a finite object of interest is imaged against a uniformly grey ( or noisy ) background .
this method can thus be efficiently exploited in brain spect imaging since the true undistorted rcbf map of a human brain consists of a finite support imaged against a noisy background ( the background depending of both the poisson noise phenomenon inherent to imaging with radioactive elements and to other , nonpoisson sources of background in the images such as electronic noise from the scanner ) .
the only information required for this deconvolution procedure is the nonnegativity   of the true image and the support of the object to be restored . in ,
this support was accurately determined by an unsupervised 3d markovian segmentation technique applied to the spect volume .         
using the same strategy , but this time during the reconstruction process , we can cite bayesian tomographic reconstruction methods [ 5 , 6 ] which use structural information on the presence and location of important anatomical  landmarks 
( such as local discontinuities or extended homogeneous regions as seen for instance on an mr anatomical image ) to control noise without smoothing edges .
these models usually express that , within a detected and segmented  uniform  anatomical region , neighboring pixels in the functional image should have similar grey level values ( local homogeneity ) or follow a gaussian distribution with a unique mean value ( global homogeneity ) [ 7 , 8 ] .              in this paper , we propose to extend the method presented in   by introducing into the nas - rif algorithm a new spatially - adaptive regularization term for spect image deconvolution .
this regularization term allows to efficiently include anatomical information extracted from a high - resolution anatomical mri image   while stabilizing the solution of the nas - rif inverse filter by preventing noise amplification and ringing artifacts . in our application
, this structural anatomy - based regularization term exploits the result of an unsupervised markovian segmentation obtained after a preliminary registration step between the mri and spect volumes from the same patient .
the proposed regularization term is quadratic and the nas - rif procedure thus involves recursive filtering of the degraded image to minimize a newly convex objective cost function with a conjugate gradient method .
our restoration method was tested on 30 pairs of brain mri and spect images from different patients and on hoffman and jaszczak spect phantoms and compared with a standard supervised deconvolution / restoration approach using a classical metz filter .         
section 2 briefly describes the proposed 3d anatomical constraint version of the nas - rif deconvolution technique .
we then show some of our experimental results on phantom and real brain spect volumes and validate the proposed model in section 5 .
in our application , and as proposed in , we assume that 3d spect images are degraded according to the following , widely - used linear model : 
 ( 1)g(x , y , z)=f(x , y , z)*h(x , y , z)+n(x , y , z ) , 
							 in which g(x , y , z ) , f(x , y , z ) , and h(x , y , z ) , respectively , denote the degraded 3d image , the true image and the point spread function ( psf ) .
n(x , y , z ) represents the additive noise and  designates the 3d discrete linear convolution operator .
the 3d blind deconvolution problem consists then in determining f(x , y , z ) and h(x , y , z ) ( or its inverse ) given the blurred observation g(x , y , z ) .            in the 3d
extended version of the nas - rif deconvolution strategy , the output of the fir filter u(x , y , z ) of dimension nxu  nyu  nzu gives an estimate of the true image f^(x , y , z ) .
each resulting estimation is passed through a nonlinear filter which uses a nonexpansive mapping to project the estimated 3d image into the space representing the known characteristics of the true image ( expressing in fact that the image is assumed to be nonnegative with a known support ) .
the difference between this projected image f^nl and f^ , e(x , y , z ) , is used as the error signal to update the variable filter u(x , y , z ) . in the 3d context , the cost function used in the deconvolution procedure of the 3d image is defined as 
 ( 2)j(u)=j1(u)+j2(u)+  j3(u ) 
							 with 
								 ( 3)j1(u)=(x , y , z)f^2(x , y , z)(1sgn  ( f^(x , y , z))2),j2(u)=(x , y , z)(f^(x , y , z)lb)2,j3(u)=((x , y , z)u(x , y , z)1)2 , 
							 where f^(x , y , z)=g(x , y , z)*u(x , y , z ) , and sgn(f ) =   1 if f < 0 and sgn(f ) = 1 if f  0 .
 is the set of all pixels of g(x , y , z ) inside the region of support , and  is the set of all pixels outside the region of support .             the first term , j1(u ) ,
is used to penalize negative voxels in the support in order to keep the image estimate nonnegative .
the second term j2(u ) penalizes voxels located outside the support with values which deviate significantly from the background average lb . when the background of the true image is black , that is , lb = 0 , the third term , j3(u ) , is used to avoid a trivial all - zero minimum solution (  being a positive constant ) .       
the authors have shown in   that the above equation is convex in the 2d case with respect to u. this property remains true in the 3d case so that convergence of the algorithm to the global minimum is ensured using the conjugate gradient minimization routine .
the major shortcoming of the nas - rif technique is its noise amplification at low snr .
this is due to the high - pass property of the inverse filter u(x , y , z ) which amplifies high - frequency noise . as a result
, the solution at convergence may not be the best estimate of the original object in the presence of noise . in order to solve this problem ,
a solution has been suggested by kundur and hatzinakos , which consists in halting the iterative restoration process through visual inspection . in practice , this requires a strong supervision and , even in this case , it is not so easy to determine which is the optimal iteration for termination ( different parts of the image may converge at different rates ) , making this method unreliable .           in this work
, we propose an alternative regularization approach for the nas - rif algorithm which can also be viewed as a way to incorporate anatomical information extracted from a ( high - resolution ) anatomical mri ( or ct ) image into the spect data .
the proposed regularization term also allows stabilization of the inverse solution by preventing noise amplification , does not require supervision ( parameters tuning or stopping criterion ) and is capable of introducing better constraints on the solution of our restoration problem .
this strategy consists in applying , over predetected and segmented anatomical regions , a piecewise smoothness constraint on the functional spect image to be recovered . to this end
, our regularization term exploits the result of a preliminary registration step between the mri and spect images as well as the result of a segmentation of the mri image into anatomical classes .           in our model , the new cost function related to the deconvolution of the 3d image is now defined as 
 ( 4)j(u)=j1(u)+j2(u)+  j3(u)+  j4(u ) 
							 with : 
								 ( 5)j4(u)=i=13  (x , y , z)ri(f^(x , y , z)ri)2 , 
							 where the first summation is made on the three main  anatomical  types ( tissues ) found in the brain , that is , white matter ( rwm ) , grey matter ( rgm ) , and cerebro - spinal fluid ( rcsf ) and ri designates the mean , in grey levels , of the ith region and  is a weighting factor between this anatomical constraint and the hard constraints of the nas - rif procedure . in this context ,
 = rwm  rgm  rcsf and ri=(1/nri)(x , y , z)rif^(x , y , z ) where nri is the number of voxels in the region ri .    

j4(u ) is proportional to the sum of variance within each anatomical region ( for each transversal slice ) of the spect image .
this term expresses that , within a detected and segmented anatomical region , pixels in the functional image in general tend to have similar grey level values .
this regularization term is edge - preserving since it allows to apply a smoothness constraint , while preserving ( anatomical ) discontinuities .            
furthermore , the introduction of this regularization term j4 does not affect the convexity of the nas - rif cost function , and therefore a unique solution to the problem is still guaranteed .
figure 2 shows the structure of this scheme . a coregistered spect and ( high - resolution ) mri datasets  along with the result of a segmentation of the mri volume into anatomical classes 
the first derivative of the cost function in ( 4 ) is shown in ( 7 ) .
the gradient vector of   j with respect to u is :  
 ( 6)j(u)=(j(u)u(0,0,0 )    j(u)u(i , j , l )    j(u)u(nxu1,nyu1,nzu1))t , 
							 where each entry is expressed as : 
 ( 7)j(u)u(i , j , l )  = 2(x , y , z)f^(x , y , z)(1sgn  ( f^(x , y , z))2 )  g(xi , yj , zl )  + 2(x , y , z)(f^(x , y , z)lb)g(xi , yj , zl )  + 2  ( (x , y , z)u(x , y , z)1 )  + 2n=13  (x , y , z)rn(f^(x , y , z)1nrn(x , y , z)rnf^(x , y , z ) )  (g(xi , yj , zl)1nrn(x , y , z)rng(xi , yj , zl ) ) . 
     a gradient - based iterative restoration algorithm or its conjugate version can be efficiently applied to minimize this convex cost function . besides , since the proposed criterion is quadratic , many other optimization methods can be used .          
the initial inverse fir filter required by the nas - rif algorithm is the kronecker delta function ; the size of this inverse filter is set to 3  3  3 pixels .
furthermore , we have used  = 0 because the background of spect images is not completely  black  .          finally , the convergence criterion of the proposed algorithm is the stability of the cost function to be minimized , that is , 
 ( 8)j(u[l+1])j(u[l])j(u[l]) , 
							 where  is a threshold , typically set in our application to 10 , with the upper - script denoting the iteration number .
figure 3 shows the evolution of the cost function value along the iterations of the gradient descent process for the image restoration presented in figure 8 .          in order to define our anatomically based regularization term j4
, we exploit the result of a 3d registration step between the mri and spect input volumes ( from the same patient ) [ 11 , 12 ] and then an unsupervised markovian segmentation of the ( registered ) mri 3d image into anatomical classes   ( see section 3.2 ) .
in order to define our anatomically based regularization term j4 , we exploit the result of a 3d registration step between the mri and spect input volumes ( from the same patient ) and then an unsupervised markovian segmentation of the ( registered ) mri 3d image into anatomical classes . the 3d registration method used in our application
is based on mutual information ( mi ) and is fully described in . the mi registration criterion c( ) between the input mri and spect volumes
describes the amount of information in the joint histogram of the images ; hence its maximization results in the best match of intensity correspondences between the images for registration .
the optimal set of registration parameters optimal is then found by maximizing c( ) , where the vector  is simply estimated by the powell 's method .
the images are smoothed slightly in order to make the cost function c( ) as smooth as possible to give faster convergence and less chance of finding bad local minima ( related to a wrong registration ) .
the code used to register the mr image to the spect image is mainly inspired from the software package statistical parametric mapping ( spm ) .
( see figure 7 ) . to this end , we consider two random fields ( x , g ) , where g = { gs , s  s } represents the field of observations located on the 3d lattice s of sites s ( voxels ) and x = { xs , s  s } the label field ( related to the class labels xs of a segmented 3d image ) .
each aforementioned label is associated to a specific brain  tissue  category or region on the 3d image .
csf and extra - cerebral tissues are combined in a single class , corresponding to tissue without ( significant ) tracer uptake .
although skin and other structures outside the brain actually have a nonzero tracer presence ( they have blood flow ) , we assume this to be negligible here .
the csf area designates regions within the brain and immediately around it that are actually devoid of activity ( intracerebal ventricles and peri - cerebral cysterns ) .
the white matter and grey matter ( brightest region ) are associated with lower and higher levels of blood flow , respectively , .
each gs takes its value in { 0 ,  , 255 } ( 256 grey levels ) , and each xs in { e1 =  csf 
, e2 =  whitematter  , e3 =  greymatter}. the distribution of ( x , g ) is defined , firstly , by a prior distribution px(x ) , assumed to be markovian and secondly , by the site - wise conditional data likelihoods pgs / xs(gs / xs ) whose shape and parameter vector (xs ) depend on the concerned class label xs ( gs designates the grey level intensity associated to the site s ) .
estimation stepin order to determine  = ( (e1 ) , (e2 ) , (e3 ) ) , we use the iterative conditional estimation ( ice ) algorithm   and three different gaussian laws for the likelihood . in order to determine  = ( (e1 ) , (e2 ) , (e3 ) ) , we use the iterative conditional estimation ( ice ) algorithm   and three different gaussian laws for the likelihood .
segmentation stepbased on the estimates given by the ice procedure , we can compute an unsupervised 3d markovian segmentation of the mr volume . in this framework , the markovian segmentation can be viewed as a statistical labeling problem according to a global bayesian formulation in which the posterior energy has to be maximized : 
									 ( 9)u(x , g)=ssln  pgs    xs(gs    xs)u1(x , g)+s , tst  ( 1(xs , xt))u2(x ) 
								 where u1 expresses the adequacy between observations and labels , and u2 represents the energy of the a priori pott model ( which tends to favor homogeneous regions with no privileged orientation ) .
we use the deterministic iterated conditional modes ( icm ) algorithm   to minimize this global energy function .
for initialization of this algorithm , we use the segmentation map obtained by a maximum likelihood ( ml ) segmentation .
the support  is then determined simply by the set of pixels belonging to csf , white and grey matter classes .
based on the estimates given by the ice procedure , we can compute an unsupervised 3d markovian segmentation of the mr volume . in this framework
, the markovian segmentation can be viewed as a statistical labeling problem according to a global bayesian formulation in which the posterior energy has to be maximized : 
									 ( 9)u(x , g)=ssln  pgs    xs(gs    xs)u1(x , g)+s , tst  ( 1(xs , xt))u2(x ) 
								 where u1 expresses the adequacy between observations and labels , and u2 represents the energy of the a priori pott model ( which tends to favor homogeneous regions with no privileged orientation ) .
we use the deterministic iterated conditional modes ( icm ) algorithm   to minimize this global energy function .
for initialization of this algorithm , we use the segmentation map obtained by a maximum likelihood ( ml ) segmentation .
the support  is then determined simply by the set of pixels belonging to csf , white and grey matter classes .
the spect images were acquired with a triple - head gamma camera ( picker prism , cleveland , oh , usa ) equipped with low - energy , high - resolution parallel - holes collimators .
the spect projections were acquired over 360. 90 projections of 50 seconds each were obtained .
mri images were acquired on a siemens magnetom avanto 1.5 t scanner using a 3d - fisp with a radial trajectory in k - space .
the scanning parameters were tr = 9.2 ms , te = 2.2 ms with n slices of 512      512 voxels with voxel dimensions of 0.5  0.5  1.0 mm , and n  [ 130 , 150 ] . these 3d mri images were further processed to isolate the brain from other tissues , using the brain extraction tool ( bet )   of the mricro software by adjusting bet 's fractional intensity threshold .
two imaging phantoms were considered in order to validate the accuracy of our spect images restoration method for spect imaging , the hoffman 3d brain phantom   was scanned while containing 148 mbq of activity .
the phantom was positioned so that the orientation of the slices within the phantom would match as much as possible those from the mri image .
phantom spect data included a set of 61 slices of 128  128 voxels with voxel dimension of 1.85  1.85  1.85 mm .
the mri phantom data contained 209 slices of 256  256 voxels with 1 mm isotropic voxels .
the other phantom contains 6 spheres of different sizes ( diameters # 1 : 9.5 mm , # 2 : 12.7 mm , # 3 : 15.9 mm , # 4 : 19.1 mm , # 5 : 25.4 mm , # 6 : 31.0 mm ) . in one condition
, the spheres were filled with a radioactive solution and the rest of the cylinder ( background ) with a less concentrated solution of the same radioisotope , giving an activity concentration ratio ( spheres / background ) of approximately 2.7 : 1 ; in another condition , the cylinder was filled with a low - activity solution and the spheres were filled with nonradioactive water .
phantom spect data included a set of 93 slices of 128  128 voxels with voxel dimensions of 1.85  1.85  1.85 mm .
the phantom mri data contained 224 slices of 256  256 voxels with 1 mm isotropic voxels .
simple visual examination is an easy method for evaluation of the restorative power of a technique , but it is obviously an insufficient approach .
a better evaluation approach consists in computing a performance measure based on the improvement in signal - to - noise ratio ( isnr ) , expressed in decibels ( db ) , using both the degraded phantom , the ground truth ( or the original undegraded image given by the mri image ) , and the restored phantom images .
the isnr is defined by , 
 ( 10)isnr=10  log  10(ioriideg  2ioriires2 ) , 
							 where ideg   is a given degraded phantom image , iori is the corresponding original ( ground truth ) phantom image and ires is the restored phantom image . ||
obviously , this metric can only be used with knowledge of the original object ; in our case this came from the mri phantom and knowledge of the radioactivity concentration within each subcompartment of the phantom .          
in addition , restored images were also evaluated by the specific evaluation criteria proposed in [ 22 , 23 ] , based on the estimation of the four following measures : the global contrast   of the image , defined by cg = ( 1  mwm / mgm ) , where mwm and mgm are the means of the pixel value in the white and grey matter areas respectively .
the local contrast of the image , defined by cl = ( ri  bj)/bj , where ri represents the mean grey level value inside the ith sphere and bj represents the mean grey level value outside the ith sphere ( in a circle centered around the sphere and whose radius d is half the distance from one sphere center to the next in the image . )
the image mottle mwm in the white matter region , defined by mwm = wm / mwm , where wm is the standard deviation of pixel values in this area .
the image mottle mgm in the grey matter region , defined by mgm = gm / mgm , where gm is the standard deviation of pixel values in this area .
these two metrics mwm and mgm allow us to measure the amplification of noise and/or to measure the presence of undesirable artifacts that could be created by the restoration procedure in a uniform region of the spect volume . due to the different number of pixels belonging to each brain anatomical tissue category , we considered the total mottle measure given by m = wmmwm + gmmgm , with wm and gm designating the proportion of pixel belonging to the white and gray matter area respectively .
a reliable spect image restoration method is one that enhances image contrast with little increase in the mottle .
inversely , for a given maximal mottle measure , we can measure if the contrast enhancement is significant .
we have compared our blind and unsupervised deconvolution approach to a classical deconvolution technique using a metz filter .
the metz filter is a supervised deconvolution ( restoration ) procedure which assumes knowledge of the point spread function ( psf ) of the imaging system .
this filter allows deconvolution of spect images while attenuating very high frequencies ( i.e. , noise which could be induced by inverse filtering )
. please remember that the metz filter has two parameters to be adjusted , the full width at half maximum ( fwhm ) related to the filter size and p ( which is the order of the filter ) .
restoration of clinical data was performed on pairs of mri and spect images from 30 epileptic patients .
each spect data set contained n slices of 128  128 voxels with voxel dimensions of 1.85  1.85  1.85 mm , and n  [ 69 , 103 ] .
mri data sets contained m slices of 512  512 voxels with voxel dimensions of   0.5  0.5  1.0 mm and m  [ 130 , 150 ] . in this section ,
isolation of the brain ( in mri ) from other tissues was made off - line prior to the registration and restoration steps .
this value was chosen empirically after a series of test runs which were visually evaluated .      

figure 4 shows the variation in isnr of the processed phantom images with respect to the observations for a range of fwhm parameter and for the best value of p. by examining the isnr for different values of fwhm , we noticed that the optimal restoration value was for fwhm = 9.85 mm , p = 40 corresponding to isnr = 0.42 db .             
average contrast and total mottle were first quantified on a set of ( human brain ) degraded and restored spect images ( with and without our anatomical - based regularization term j4(u ) ) and compared to the metz filter .
our proposed algorithm with j4(u ) resulted in an increase of global contrast by 1.87 and of mottle by an acceptable factor of 1.17 .
the metz filter increased global contrast by 1.52 and mottle by a factor of 1.08 .             for spect phantoms ,
for all the tested experiments , we used  = 0 because the background of spect images is not completely  black    and we chose 25 for the weighting factor  for j4 .
this value was chosen empirically after a set of experiments by varying  in the interval [ 10 , 100 ] .         
spect images of the cold spheres phantoms were also restored using partially incorrect anatomic information ( i.e. , cylindrical phantom mri with 5 different sphere sizes instead of 6 ; one sphere was removed in the mri segmented image ) .
table 2 shows that the local contrast measured on this restored spect image ( with our anatomical - based regularization term j4(u ) ) remained similar to the local contrast obtained with correct anatomic information .           
local contrast was quantified on individual spheres of spect degraded and restored images ( with and without our anatomical - based regularization term j4(u ) ) and compared to those obtained with the metz filter .
table 3 shows that our algorithm with j4(u ) resulted in a average value of the increase in local contrast of 2.07 .
the metz filter increased , in average , the local contrast by only 1.53 .          
figures 8   and 9   show , respectively , examples of the hoffman phantom and brain spect volumes obtained with our restoration method .
the reconstructed whole brain volume converged to a very good estimate of the solution without   a priori information about the psf and allowed to noticeably improve the contrast of the original , unrestored spect volume .
this restoration should allow more efficient detection of small , localized singularities associated with different types of lesions ( tumors , epileptogenic foci , etc . ) that often are not clearly visible in the original blurred image .             
the metz filter has two parameters to be adjusted , fwhm and order p. we have used the values of 4.5 , 5 , 5.5 , 6 , 6.5 , 7 , 8 , 9 , 10 and 11 mm for fwhm and 10 , 20 , 30 , 40 and 50 for order p , a total of 50 combinations of both parameters .
all the above values were tested and only the best one ( in terms of isnr ) was selected for comparison with our algorithm .             
in addition , the 30 pairs of brain mri and spect acquisitions from different subjects have been anonymised using jim 's dicom anonymiser tool .
( this tool allows the removal from image files of information that may jeopardise patient 's or physician 's privacy . )
our restoration procedure was not highly sensitive to the bet 's fractional intensity threshold value when this was visually set by an experimented user ( e.g. , the value of contrast was stationary when the threshold went from 0.45 to 0.55 ) .
therefore , this parameter was set to 0.50 for all thirty pairs of mri images and seemed to be optimal in all tested cases .          
notice that the degradation model used by our restoration method assumes that noise is additive while in reality it is a multiplicative poisson process . nevertheless , the results with real spect / mri data are visually impressive .
it seems that at the high count level used in spect , the difference between ( multiplicative ) poisson and ( additive ) gaussian noise does not affect significantly the efficiency of the algorithm
. however it would be meaningful to conduct clinical studies using roc analysis to better assess this restoration procedure .
when compared with a classical restoration approach using a metz filter , our method performed better , in terms of signal - to - noise ratio .
the slight increase in mottle with j4 is a limited price to pay given the gains in global contrast , isnr , and local contrast measures , which from a clinical perspective ensure better detection of focal anomalies , a task at which spect is usually notoriously poor .          
the tests also showed that the restored image is constrained ( or guided ) by our prior knowledge but not wrongly biased by this information while avoiding noise amplification inherent to each deconvolution process .          
the computational time of our technique was approximatively 6.33 minutes as compared to 0.24 minutes for the metz filter on a 2.0 ghz pc workstation running linux .
our method is computationally demanding but remember that it is an unsupervised one which does not require any psf parameter .         
another disadvantage of our method is the need to get an mri scan ; however , many of the subjects submitted to a spect rcbf study actually require one from a clinical perspective .
importantly , the accuracy of the registration procedure is crucial for the final restoration result , fortunately many highly accurate registration procedures are widely available elsewhere .
in this paper , we have presented a robust restoration method using anatomical and functional information fusion with application to spect images .
this method improves the quality of human brain 3d spect images and should thus be helpful to physicians interpreting such studies .
our approach takes advantage of the anatomical information contained in the mri study of each subject ( or any other high - resolution image such as ct ) .
the proposed constraint term allowed both stabilization of the inverse solution of the nas - rif procedure by prevention of noise amplification and the generation of a better constraint on the solution of our restoration problem . in the regularization framework , this term allowed smooth regions to be reconstructed in the spect image , where such homogeneous anatomical regions are found in the high - resolution mri images of the patient , after those were registered to the subject 's spect volume .
this method was tested on a number of spect / mri pairs , demonstrating its efficiency and robustness .
this 3d blind restoration technique is completely data driven and automated , therefore it could be implemented to process large numbers of 3d spect studies .
 image restoration is usually viewed as an ill - posed problem in image processing , since there is no unique solution associated with it . 
 the quality of restored image closely depends on the constraints imposed of the characteristics of the solution . in this paper 
 , we propose an original extension of the nas - rif restoration technique by using information fusion as prior information with application in spect medical imaging . 
 that extension allows the restoration process to be constrained by efficiently incorporating , within the nas - rif method , a regularization term which stabilizes the inverse solution . 
 our restoration method is constrained by anatomical information extracted from a high resolution anatomical procedure such as magnetic resonance imaging ( mri ) . 
 this structural anatomy - based regularization term uses the result of an unsupervised markovian segmentation obtained after a preliminary registration step between the mri and spect data volumes from each patient . 
 this method was successfully tested on 30 pairs of brain mri and spect acquisitions from different subjects and on hoffman and jaszczak spect phantoms . 
 the experiments demonstrated that the method performs better , in terms of signal - to - noise ratio , than a classical supervised restoration approach using a metz filter . 
1. Introduction
2. 3D Anatomical Constraint NAS-RIF Algorithm
3. Registration and Segmentation
4. Validation
5. Experimental Results
6. Discussion
7. Conclusion