 Open Access
 Total Downloads : 1771
 Authors : S.Asif Hussain, K. Lokeswara Reddy
 Paper ID : IJERTV1IS1009
 Volume & Issue : Volume 01, Issue 02 (April 2012)
 Published (First Online): 29042012
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
An Efficient CT Image Reconstruction with Parallel Modeling for Superior Quantitative Measures
S.Asif Hussain1 K. Lokeswara Reddy2
Associate professor M.Tech (E.S) student Department of ECE, A.I.T.S, Department of ECE, Rajampet, Andhra Pradesh, India A.I.T.S Rajampet, A.P, India email: sah.ssk@gmail.com e mail: lokes9@gmail.com
ABSTRACT
Image segmentation algorithms based on ROI typically rely on the homogeneity of image intensities. CT scanner is dedicated as research Scanner which has been developed in view of imaging applications. A key Feature of the work is to use Empirical system to achieve resolution recovery with novel region based method. This method identifies local intensity cluster with local clustering criterion function with respect to neighborhood center. Reconstruction quality is analyzed quantitatively in terms bias field correction for intensity inhomogenity correction. This method is valid on synthetic images of various imaging modalities. A significant improvement in reconstruction quality can be realized by faster and more accurate visual quality quantitative measures where Reconstruction quality is analyzed quantitatively in terms of biasvariance measures (bar phantom) and mean square error (lesion phantom). However, with the inclusion of the empirical kernel, the iterative algorithms provide superior reconstructions compared to FBP, both in terms of visual quality and quantitative measures. Simulated results show improved tumor bias and variance characteristics with the proposed algorithm.
Keywords: Intensity inhomogeneities, Empirical system kernel, Biasvariance, Iterative algorithms

INTRODUCTION
Image segmentation is often an essential step before further imageprocessing of threedimensional medical images canbe done. An object can be segmented based on shape and/or intensity characteristics. The task of image segmentation can be simplified with initialized parameters to guide accurate segmentation. Semiautomated and interactive methods [1] have been relatively successful, but require varying degrees of human input. Segmentation is often an important step in US Bmode image analysis .we consider the problem of correcting for attenuationrelated intensity inhomogenieties i.e., those that cause a slowly changing (lowfrequency) intensity contrast and are not due to specklemode imaging artifacts include speckle noise, attenuation(absorption and scattering), etc. The statistical analysis and reduction of speckle noise has been studied extensively in the literature [1][7]. Other artifacts, particularly those caused by non uniform beam attenuation within the body that are not accounted for by time gain compensation (TGC), also decrease the image signaltonoise ratio (SNR).Existing level set methods for image segmentation can be categorized into two major classes: regionbased models [4], [10], and edge based models [3], [7], [8], [12]. Regionbased models aim to identify each region of interest by using a certain region descriptor to guide the motion of the active contour. However, it is very difficult to define a region descriptor for images with
intensity inhomogeneities. Most of regionbased models [4],
[16] are based on the assumption of intensity homogeneity. Edgebased models use edge information for image segmentation. These models do not assume homogeneity of image intensities, and thus can be applied to images with intensity inhomogeneities. A novel regionbased method for image segmentation. From a generally accepted model of images with intensity inhomogeneities, we derive a local intensity clustering property, and therefore define a local clustering criterion function for the intensities in a neighbourhood of each point. This local clustering criterion is integrated over the neighbourhood center to define an energy functional, which is converted to a level set formulation. Minimization of this energy is achieved by an interleaved process of level set evolution and estimation of the bias field. As an important application, our method can be used for segmentation and bias correction of magnetic resonance(MR) images. 
Methods:
BACKGROUND:
In this section, we review the method proposed by Zhang for estimating the field distortion and simultaneously segmenting an MR image and provide implementation details on how it has been adapted to work with US images. This method essentially estimates the low (spatial)frequency multiplicative degradation field while at the same time identifying regions of similar intensity inhomogeneity using an MRFMAP frame work. As we will explain in Section III, although developed for another imaging modality, under simplified assumptions, we can justify using the same approach on displayed US images.

Model Specification:
Let S be a lattice indexing the pixels in the given image. Further,
let and be the observed and the ideal (that is, without intensity inhomogeneity distortion) intensities of the given image respectively, being the number of pixels in the image. We assume that the distortion at pixel can be expressed by a multiplicative model of the form
Whererepresents the gain of the intensity due to the intensity inhomogeneity at pixel . A logarithmic transformation of this equation yields an addition. Let anddenote, respectively ,the observed and the ideal log transformed intensities, then
Where denotes the logtransformed intensity distortion field. Segmentation can be considered as a problem of statistical classification, which is to assign every pixel a class label from a label set. Let denote the label set. A labeling
of will be denoted by in is the corresponding class label of Given the class label , it is assumed that the intensity value at follows a Gaussian distribution(this assumption will be justified in
Section III) with parameter being the mean and the variance of class respectively
where
With the distortion field taken into account, the above distribution can be written in terms of the observed intensity as
(4)
and, hence, a classindependent intensity distribution
(5)
Thus, the intensity distribution at is modeled as a Gaussian mixture, given the distortion field. Assuming that the pixel intensities are statistically independent, the probability density for the entire image, given the distortion field,
(6)
Bayes rule can be used to obtain the posterior probability of the distortion field, given the observed intensity values
(7)
Where is a normalization constant. The prior probability density of the distortion is modeled as a Gaussian with zero mean to capture its smoothness property. The maximum a posteriori (MAP) principle can be employed to obtain the optimal estimate of the distortion field , given the observed intensity values
. (8)
The optimum satisfies the following condition
. (9)
Solving this equation leads to the update equations (see [12] for detail)
, with I= (1,1,,1 . (11)
Here, is the posterior probability that belongs to class given the distortion field is a lowpass filter, is the mean residual in which for pixel
(12)
And is the mean inverse covariance, in which if otherwise.
(13)

Bayesian criterion for filtering edge information:
It is assumed that interpolated boundaries will partially overlap with the true edges found using edge detection. The probability of edges overlapping with shapeinterpolated boundaries may be modeled using Bayes probability. It is assumed that the probability of overlap at interpolated slces is greater than or equal to that at userinitialized contours. Edges are divided into edge components based on their connectivity. To retain edges with higher saliency, the edge components are sorted in descending order relative to the amount of overlap with the boundary. When the cumulative probability of overlap exceeds that obtained from user initialized contours, the remaining edge components are discarded. The Bayes classification is thus not employed for training, but rather as a guide to how well boundaries can be defined based on edge detection.


Proposed Algorithms:

The 3D Case:
The algorithm can be applied to 3D volumes reconstructed from a sequence of parallel, closely spaced 2D images. We assume that in such a sequence neighbouring slices resemble each other, that is, overlapping pixels in neighbouring slices tend to have the same class labels. Intensity in homogeneity field estimation is performed within each 2D image, while the energy function in the MRF prior model involves a 3D neighbourhood system, which includes, for each pixel in a 2 D scan, the eight nearest neighbours in the same scan and the two direct neighbours in the previous and the next scan. This 3D constraint helps to strengthen ambiguous boundaries that are easily mislocated in2D processing.

Boundaryedge correspondence:
Ideally, the match between boundary and edge should be one toone. However, deviations in the interpolated shape will not initialize Bi well. To prevent manytoone snapping of boundary points, a minimum snappingdistance map is stored for every edge point. Subsequent boundary points will only be allowed to snap to the edge point if the snapping distance is less than or equal to the value in the minimum snapping distance map. Therefore, boundary points will not arbitrarily snap to false edges if there are no edges to be found. During
the first iteration, the search window has not been adaptively altered to match the edge proximity for the image slice. There is a possibility that a false edge will be included in the Bayesian criterion. To prevent this, an inverse weighted distance transform, M, is multiplied to Fi,k, where M is a square matrix. Denoting pq M as an element in M and any two points on the Bi as bp and bq.,pq M is defined in Eq. 1.
(1)

Local Intensity Clustering Property:
Regionbased image segmentation methods typically relies on a specific region descriptor (e.g. intensity mean or a Gaussian distribution) of the intensities in each region to be segmented. However, it is difficult to give such a region descriptor for images with intensity inhomogeneities. Moreover, intensity inhomogeneities often lead to overlap between the distributions of the intensities in the Therefore, it is impossible to segment these regions directly based on the pixel intensities
.Nevertheless, the property of local intensities is simple
,which can be effectively exploited in the formulation of our method for image segmentation with simultaneous estimation of the bias field based on the image model in (3) and the assumptions A1and A2, we are able to derive a useful property of local intensities, which is referred to as a local intensity clustering property as described and justified below. To be specific, we consider a circular neighbourhood with a radius centered at each point defined by .The of the entire
domain induces a partition of the neighbourhood
, i.e., forms a partition of . For as lowly varying bias field the values for all in the circular Neighbourhood are close to i.e. for
(4)
Thus, the intensities in each sub are close to the constant i.e.
(5)
Then, in view of the image model in (3), we have
Where is additive zeromean Gaussian noise. Therefore, the intensities in the set
form a cluster with cluster which can be considered as samples drawn from a Gaussian distribution with mean . Obviously , the clusters ,, , are wellseparated ,with distinct cluster centers
, (because the constants
are distinct and the variance of the Gaussian noise is assumed to be relatively small). This local intensity clustering property is used to formulate the proposed method for image
segmentation and bias field estimation as follows.

Energy Formulation:
The above described local intensity clustering property indicates that the intensities in the neighbourhood can be classified in to clusters with centers
, This allows us to apply the standard Kmeans clustering to classify these local intensities. Specifically, for the in the neighbourhood
, the Kmeans algorithm is an iterative process to minimize the clustering criterion [19], which can be written in a continuous form as
(6)
Where is the cluster center of the cluster, is the
mem and for Since is the membership function of the region , we can rewrite as
dx. (7)
In view of the clustering criterion in (7) and the approximation of the cluster center by
we define a clustering criterion for classifying
the intensities in as (8)
Where is introduced as a nonnegative window
function, also called kernel function, such that for .With the window function, the clustering criterion function can be rewritten as
(9)
This local clustering criterion function is a basic element in the formulation of our method.

Multiphase Level Set Formulation:
For the case of we can use two or more level set functions to define membership functions of the regions , such that
.
For example, in the case of , we use two level set functions to define
And to give a threephase level set formulation of our method. For the fourphase case ,
the definition of can be defined as
functional in our multiphase level set formulation is defined by
. (25)
The minimization of the energy in (25) with respect tto the variable can be performed by solving the following gradient flow equations:
and
For notational simplicity, we denote these level set functions by a vector valued function .
Thus, the membership functions can ) be written as . The energy in (10) can be converted
to a multiphase level set formulation .
With given by .
(16). For the function .
This defines the regularization Terms and ,
Where and are defined by (19) and (20) for
each level set function , respectively. The energy .


IMPLEMENTATION AND SIMULATION RESULTS:
Fig 3. Iterations are Performed to the Input Image
Fig 1. showsInput Of CT Image of bones
Fig 2. gives Histogram of CT Image With Density
50 iterations
Fig 4. CT image with 50 Iterations
Fig 5. Blurring Of Input Image for Checking from Initial Point of the Image
Bias corrected image
Fig 6. Input Image Of Bias Corrected Image
Histogram of Bias corrected image
200
180
160
Author Profile:
Lokeswar Reddy.Kokatam received B.Tech Degree in Electronics& Communication Engg. from JNT University, Annanthapur, India. Presently he is with Annamacharya Institute of Technology & Sciences,
Rajampet, A.P., India in Dept. of ECE and pursuing his M.Tech. His research interests include Signal Processing, Time Series Analysis and Image Processing.
Asif Hussian.Shaik received B.Tech & M.Tech Degree in Electronics & Communication Engg. from JNT University, Hyderabad, India. He is currently working towards PhD Degree in Biomedical Image
Processing at JNTU University, Annantapur, India. Presently he is with Annamacharya Institute of Technology & Sciences, Rajampet, A.P., India. He is working as Assistant Professor in Dept. of ECE. He presented many research papers in National & International Conferences & journals. He is a member of Professional societies like ISTE (India), BMESI (India), IACSIT (Singapore), IAENG (Hongkong) and WASE (Hongkog).His research interests include Signal Processing, Time Series Analysis and Image Processing.
7. REFERENCES
140
120
100
80
60
40
20
0
0 50 100 150 200 250

Olabarriaga, S.D. and Smeulders, A.W.M., Interactionin the Segmentation of Medical Images: A Survey,Med. Image Analysis, 5: 127142, 2001.

Osher, S. and Sethian, J.A., Fronts Propagating withCurvature Dependent Speed: Algorithms Based onHamiltonJacobi Formulations, J. Comp. Physics, 79:12 49, 1988.

P. N. T.Wells and M. Halliwell, Speckle in ultrasonic imaging, Ultrasonics,vol. 19, pp. 225229, 1981.
Fig 7. Image describes Histogram of Bias Corrected Image


Conclusion:
This work presents a variational level set framework for segmentation and bias correction of images with intensity inhomogeneities. Based on a generally accepted model of images with intensity inhomogeneities and a derived local intensity clustering property, the work defines an energy of the level set functions that represent a partition of the image domain and a bias field that accounts for the intensity inhomogeneity. Segmentation and bias field estimation are therefore jointly performed by minimizing the proposed energy functional. The slowly varying property of the bias field derived from the proposed energy is naturally ensured by the data term in our variational framework, without the need to impose an explicit smoothing term on the bias field. The proposed method is much more robust to initialization than the piecewise smooth model. Experimental results have demonstrated superior performance of our method in terms of accuracy, efficiency, and robustness.

ACKNOWLEDGMENTS:
The work was supported by my guide S.Asif Hussain from Annamacharya Institute of Technology & Sciences,,Rajampet ,India under Research grants of R.P.S A.I.C.T.E, New delhi.

A. N. Evans and M. S. Nixon, Biased motionadaptive temporal filteringfor specklereduction in echocardiography, IEEE Trans. Med.Imag., vol. 15, pp. 3950, Feb. 1996.

V. Caselles, R. Kimmel, and G. Sapiro, Geodesic active contours,Int. J. Comput. Vis., vol. 22, no. 1, pp. 6179, Feb. 1997.

T. Chan and L. Vese, Active contours without edges, IEEE Trans.Image.Process., vol. 10, no. 2, pp. 266277, Feb. 2001.

S. Kichenassamy, A. Kumar, P. Olver, A. Tannenbaum, and A. Yezzi,Gradient flows and geometric active contour models, in Proc. 5th Int.Conf. Comput. Vis., 1995, pp. 810 815.

R. Kimmel, A. Amir, and A. Bruckstein, Finding shortest paths onsurfaces using level set propagation, IEEE Trans. Pattern Anal.Mach.Intell., vol. 17, no. 6, pp. 635640, Jun. 1995.

C. Li, C. Kao, J. C. Gore, and Z. Ding, Minimization of regionscalablefitting energy for image segmentation, IEEE Trans. ImageProcess., vol. 17, no. 10, pp. 19401949, Oct. 2008.

R. Malladi, J. A. Sethian, and B. C.Vemuri, Shape modeling with frontpropagation: A level set approach, IEEE Trans. Pattern Anal. Mach.Intell., vol. 17, no. 2, pp. 158175, Feb. 1995.

R. Ronfard, Regionbased strategies for active contour models, Int.J.Comput. Vis., vol. 13, no. 2, pp. 229251, Oct. 1994.

C. Samson, L. BlancFeraud, G. Aubert, and J. Zerubia, A variationalmodel for image classification and restoration, IEEE Trans. PatternAnal. Mach. Intell., vol. 22, no. 5, pp. 460472, May 2000.

S. Theodoridis and K.Koutroumbas, Pattern Recognition. NewYork:Academic, 2003.

A. Tsai, A. Yezzi, and A. S.Willsky, Curve evolution implementationof the MumfordShah functional for image segmentation, denoising,interpolation, and magnification, IEEE Trans. Image Process., vol. 10,no. 8, pp. 11691186, Aug. 2001.

A. Vasilevskiy and K. Siddiqi, Fluxmaximizing geometric flows,IEEE Trans. Pattern Anal. Mach. Intell., vol. 24, no. 12, pp. 15651578,Dec. 2002.

L. Vese and T. Chan, A multiphase level set framework for imagesegmentation using the Mumford and Shah model, Int. J. Comput. Vis.,vol. 50, no. 3, pp. 271293, Dec. 2002.