Enhancing Surgical Visualization Health And Social Care Essay

Abstract - In recent years, lung tumor diagnosis and the projection of tumor segmentation in 3D has gained significant momentous in the therapeutic field. Establishing the dissimilarity exists in the three dimensional volume representation of tumor cells affords more information, which can sharpen the treatment of a multiplicity of tumors. The volume reconstruction information is indispensable in the case of surgical operations. We adduce a contour based segmentation algorithm to acquire the appropriate differentiation of pixel boundary that scrutinizes the exact difference between tumor and non tumor cells along the tumor boundary. With the aid of aforementioned formulation, extracted tumor part pixels are reconstructed for the entire 2D slices of the patient data set. Proposal work on 3D voxel reconstruction relies on encountering the isosurfaces. Originally, volume data are subjected to the smoothening process which computes the isosurface data from the smoothened volume data. The generated outcome of this process comprises the vertices and faces of the isosurfaces and directly flows to patch the data. Exploit the 3D reconstructed model to enumerate the voxel damaged by tumor. Proposal work associated with the percentage of damaged voxel along with accurate and reliable perception, simplifies the physician task in lung tumor diagnosis and assist the surgical procedure. Experimental evaluations across the wide range of images show the superiority of our proposed work.

Index terms – Image segmentation, Active contours, Level set, 3D voxel reconstruction.

I.INTRODUCTION

For the past few decades, Lung cancer is the leading cause of cancer death. Early diagnosis of lung cancer can maximize the survival rate of tumor influenced person. Therefore, poor long term survivorships of lung cancer increase the obligation for novel methodologies and technologies with respect to lung tumor diagnosis in early stage. When the tumor diagnosing methodologies are hybridized with a utilization of CT (Computerized Tomography) lung screening, then the sensitivity and accuracy of primitive stage tumor detection can be relatively enhanced. Image segmentation and enhancement techniques are exploited to accomplish the early stage tumor diagnosis more accurately. Image segmentation techniques segregate the input image into homogeneous regions, which affiliates some universal properties. This technique is categorized as region based (examine the uniformity within the sub region) and edge based segmentation (examine the discontinuity in the image intensity). Numerous techniques have been introduced to enhance the performance of image segmentation.

To solve the image segmentation problem, we employ contour based image segmentation algorithm in this paper, which is supported by Active contour and Level Set equations. Contour based segmentation algorithm efficiently acquires the exact differentiation of pixel boundary. In this segmentation algorithm, user defines the fundamental assumption for the contour and is moved by image driven forces to the contour of specified objects. Contour based segmentation model accommodate internal and external forces. During the deformation process, internal forces specified within the curve aids to keep the model smooth whereas external forces are specified to precede the model toward an object boundary. These external forces are computed only when the underlying image data is known.

In the implementation of active contours, the avail of level set method produce more convenience and flexibility. Level set method is employed to embed a curve within the surface. It acts as a theoretical and numerical tool for propagating interfaces. The intrinsic idea is to begin with a closed curve in 2D/3D and allow the curve to move perpendicular to itself at a defined speed. This way of propagation of contour by exploiting the image properties makes level set method as a frequently used tool in image segmentation.

Three dimensional (3D) reconstructions of medical images have gained an increasing attention. Furthermore, reconstruction of 3D volumetric data from 2D slices of patient data set is now an integral part in the field of biomedical research. The voluminous 2D slices registration is of extreme significance for the exact morphometric analysis (volume and surface representation) and 3D visualization. The fundamental principle of 3D reconstruction methods from 2D images is classified into assortion as follows: feature based methods using contours, gray level based registration techniques using the whole image intensities, characteristic points or crest lines elicited from the images, fiducial marker based registration methods. Effective 3D volume reconstruction of image is essential to enhance the pathological visualizations. Moreover, it promotes the value added service and assists the physician.

Our work focuses on the lung tumor segmentation and its projection in 3D volume reconstruction. The volume reconstruction information is indispensable in the case of surgical operations. Proposed work introduce contour based segmentation algorithm that acquire the exact dissimilarity between tumor and non tumor cells along the pixel boundary of the tumor. This segmentation algorithm is supported by active contour and level set equations. The formulated curve evolution of active contour and level set method aids to reconstruct the extracted image pixels. The contribution of our proposed work on reconstruction of 3D voxel relies on detection of isosurfaces. The isosurfaces data are computed through smoothening process. It computes a tri-surface geometry of a given volume data at a given isosurfaces value. The outcome of this phase contains the faces and vertices of the isosurfaces that can be passed directly to patch data. The 3D reconstructed model is then manipulated to compute the voxel impaired by tumor. We are focusing towards a system that builds 3D volume reconstruction an accurate 3D visualization of volumetric data. A complete system would comprises the enumerated components: Lung dataset perception, Enhancement of perception, Preprocessing and feature extraction, Classification, Image segmentation, Patch assignment, 3D visualization of segmented tumor portion, Tumor volume calculation. The remainder of the study is composed as follows: Section 2 surveys the related works. Proposed system in 3D reconstruction of volumetric data from 2D slices of CT lung is discussed in section 3. Section 4 analyzes the experimental evaluation of our proposed work. Finally, section 5 sums up the paper with conclusion and path ahead.

II. RELATED WORKS

Medical image segmentation based on the paradigm of deformable contour was illustrated by Domagoj K, et al. [1]. Compare to the traditional segmentation method like region or edge based methods, deformable contour based image segmentation retain many benefits. Deformable contour paradigm is implemented using level set with automatic establishment of active contour original position. This way of image segmentation overcome the limitation occurred in classical region or edge based methods.

M. Airouche, et al. [2] explores the active contour based image segmentation to detect the oil spills in images. Level set method based on partial differential equation denotes the spill surface of the image into propagation interface. Furthermore, speed function is also created by level set method to discover the propagation interface position. Segmentation method exploited in this paper produce more reasonable segmentation results than the classical segmentation. It is also possible to divide the contour into multiplex contours with the aim of accomplishing good flexibility in the use of active contours.

In medical image analysis, preprocessing is the fundamental stage in image segmentation. Computer aided diagnosis of lung CT image has gained an increasing attention in the early detection of lung disease. The preliminary step in this kind of diagnosis involves in segmenting the region of interest like lung, heart, brain etc. The segmentation algorithm proposed in [3] accurately segment the lung parenchyma of lung CT images, which assist the physicians in early lung disease diagnosis. Region growing segmentation algorithm selects the seed pixel from the dark side of the image. Once the dark pixels coordinate are selected, pixels which exist in 4-connected neighborhood of the seed pixel are marked. The original gray scale image is segmented into the lungs and trachea after the analysis of region growing.

Another Contour based Segmentation technique using selective visual attention was proposed by Engin M, et al. [4]. Visual attention is the term of selecting and getting visual information based on knowledge about objects, scenes and their intersection. The Chan Vese active contour model implemented by this author defines the attended location as initial contours. Saliency toolbox employed in this paper extracts the attended locations. This method of image segmentation is effective for accurate detection and extraction of ROI (Region of Interest). The integration of these techniques maximizes the speed of segmentation process.

In [5], region based technique presented for image segmentation integrates the information regarding the curvilinear continuity, color and intensity. The soft contour information is obtained through energy orientation whereas weak contrast gaps are completed through contour propagation.

In [6], novel boundary detection and image segmentation technique based on edge flow was presented by W.-Y.Ma, et al. The direction of color change and texture at each image location is identified by predictive coding model and then edge flow vector is constructed. In this edge flow technique, identification and localization of edges are performed indirectly. Originally, the flow direction at each pixel location, which points to the closest boundary, is identified by iteratively propagating the edge flow. Single scale parameter has been used for the implementation of this image segmentation based on edge flow. This scheme attains general framework for combining various kinds of image information for boundary detection.

Another framework for region based image segmentation methodology using Bayesian Level Sets are formulated by E. Sifakis, et al. [7]. In this framework, statistical approaches are applied for modeling the different regions. This framework is the extension of the level set approach to a multi label framework, which permits the speed of propagation depend on the corresponding region label. Based on the label content description and the region capacity, performance of segmentation is predicted. Propagation speed is maximized by exploiting posterior probability of the corresponding label. Proposed posteriori probability technique yields the optimal segmentation solution, which clearly biases the segmentation result.

C.-L. Huang, et al. [8] proposed a level set method with shape model for image segmentation. The initial shape model is established by manual processing rather than automatically. But depending on the application, shape model is initialized. The initialized shape model is dynamically resized and translated. The distance map of shape model contour is then established by dilation operation. From the obtained shape model, the level set function is updated until the contour converges. This kind of specific shape image segmentation algorithm is widely applicable to temperature image, medical image, etc.

The most widely used image segmentation method is the active contour, which provides sub region with continuous boundaries. The applicability of level set theory in the active contour implementation adds more flexibility in image segmentation. C. P. Lee, et al. [9] proposed a robust image segmentation framework, which segment the image based on statistics of image intensity. Three active contour models are proposed depends on the method to estimate the mixture density functions. Unsupervised multidimensional histogram method, Supervised multivariate Gaussian mixture density method, half supervised multivariate Gaussian mixture density method. Integration of both statistical pattern classification and image segmentation in this framework provides autonomous pattern classification.

A novel level set approach for smooth isosurface extraction was proposed by L. Linsen, et al. [10]. This approach works well for both structured and unstructured volume data. Computation of small neighborhood around the zero isosurface of the level set function is mainly concentrated in this novel framework. Since, it provides gradual smooth evaluation of the isosurface during computations. Significant isosurfaces are substantial for provoking advisory visualizations of volume data.

Numerous volumetric reconstruction methods were developed in literature. Among those, reconstruction volumetric 3D model was developed by P. Eisert, et al. [17]. W.E.Lorensen, et al. [18] proposed a high resolution 3D surface construction algorithm. S. Krinidis, et al. [11] proposed efficient algorithm for 3D volume reconstruction from acquired 2D slices. The alignment problem is considered globally on 3D volume which in turn reduces the global objective function. The algorithm overcomes the main limitation of 3D image alignment.

Similarly, a novel framework for the reconstruction of 3D rat brain volume from 2D histological images was presented by U. Bagci, et al. [12]. This framework works on the basis of three premises. 1. Standard histological images provide accurate 3D volume reconstruction, 2. Reliable feature space only is considered for successful registration in consecutive slices, 3. Proper selection of reference slice prevents high noise slices and distortion. SDM (Standard Deviation Maps) introduced in this framework estimate the smoothness of 3D volume reconstruction.

R. Li, et al. [13] developed a new algorithm based on PCA (Principal Component Analysis) and demonstrates volumetric reconstruction of 3D images. Using this reconstruction algorithm, information about 3D tumor motion is extracted in real time from a single x-ray projection. Such way of 3D volume reconstruction yields more accuracy and efficiency. Likewise, S. Dogan, et al. [14] used software called Medical Image Processing and Analysis System for reconstruction of 3D models of human body by using digital images, CT and MR slices. The fundamental problems in 3D image reconstruction are discussed by V. Cong, et al. [15].

In [16], optical flow based registration method was presented by H. Handels, which enhance the quality of image reconstructed from 4D CT data sets. The surface points analyzed from 4D CT data sets associated with 3D computation of lung tumor center displays the appearance of lung tumor in three orthogonal directions.

This section reviewed the progress and limitations of related works. Examining all those complexities exist in literature, we propose 3D volumetric data reconstruction from 2D slices of a CT lung image, which aids to accurate lung tumor 3D visualization.

III. PROPOSED METHODOLOGY

Constructing early detection of lung tumor is substantial to ensure the most convenient treatment regime is considered. 3D volumetric data reconstruction works on the principle of physical mechanism, which displays points of light within the volume. The reconstructed 3D content uses the display systems as either single rotating display panel or multi-planar display panels. Such 3D volume reconstruction surmounts the problems imposed by the 2D visualization of data. In this paper, we developed 3D reconstruction of volumetric data from 2D slices of a CT lung image, which investigate the tumor slices based on the voxel impaired by tumor. The outer section of proposed system is depicted in figure 1.

Fig 1: Proposed System Design

3.1 Preprocessing Stage

Preprocessing is the most substantial task in medical image processing, aids to enhance the image quality. Preprocessing and enhancement phase in medical image processing, regenerate the input CT lung image into standard format with enhancing image contrast, enhancing the sharpening of edges, noise reduction through background removal, image filtering and elimination of film artifacts. The effectiveness of feature extraction depends highly on the image quality. Denoising, a preprocessing technique employed in this paper denoises the CT lung image by Kernel Based Non-Local Neighborhood denoising method. Turkey-bi-weight is the supreme denoising function exercised in our proposed work, which denoise the lung image and affords noise free image. Turkey bi-weight denoising function furnishes lesser MSE (Mean Square Error) and the higher PSNR (Peak Signal to Noise Ratio). Hence, we accord this kind of denoising function to accomplish denoised image. Originally, we acquire the lung image from the dataset then the labels and marks on the image are removed and finally discharge the high frequency components in the image. Preprocessing stage generates the noise reduced, and filtered edges of given input lung image.

3.2 Feature Extraction

Extracting significantly important features from the image data is the ultimate aim in image processing applications. The denoised image obtained from preprocessing technique is supplied into the second order histogram, a feature extraction technique deployed in this work. In our proposal, we consider some valuable characteristics of images like energy, variance, contrast, entropy and homogeneity. These criteria prevail to select the subset of the most essential robust features, which simplify the classification result and minimizing entire perplexity.

3.3 Lung Image Classification

The robust image features extracted in feature extraction stage is fed up for classification to distinguish the normal and abnormal status of lung image. In our proposal, we classify the lung image by exploiting Multivariate Multinomial Distributed Bayes Classification. Multinomial Distributed Bayes Classification is one of the two Naive Bayes variants used in image classification, which holds the algorithm for multinomially distributed data. For each class x, the distribution framework is represented by vectors, here k denotes the total number of features in lung image classification, implies the probability of feature i present in a sample belonging to class x. is the metrics measured by a smoothened version of maximum likelihood.

(1)

Here, denotes the number of times feature i reflects in the training set (T), implies the total count of all features for class x, is a constant.

In Multivariate Distributed Bayes Classification, distributed data is classified corresponding to multivariate distributions. Precisely, multiple features appear in training data are chosen and make an assumption of each variable as binary valued vector. The classification rule for Multivariate Distributed Bayes Classification follows

(2)

The multivariate decision rule in a class x explicitly chases the non-occurrence of a features i, whereas decision rule based on multinomial classification neglect a non-appearing feature. In our classification algorithm, we evaluate both the models.

Training phase of naive bayes classifier estimates the parameters of a probability distribution using the training samples, by making assumptions of conditionally independent features. Prediction phase of classifier estimates the posterior probability of unseen test sample and then distinguish the test sample based on greatest posterior probability. Adduced work is trained with 150 slices (51 normal, 99 abnormal) for lung tumor diagnosis and segmentation. The ultimate goal of classification is to examine the accuracy rate of generating results. For that purpose, accuracy metrics like Sensitivity, Specificity, Error rate, correct rate, Last error rate, Last correct rate, classified rate, Prevalence, Positive/Negative Likelihood are validated.

3.4 Contour based Segmentation Algorithm

Boundary detection or image segmentation is the most fundamental process in medical image analysis. We employ active contour based segmentation algorithm, ultimate technique for detecting boundary in lung CT images. Our segmentation algorithm incorporates active contour and level set model, which acquire the exact dissimilarity of pixels boundary.

The main premise considered with respect to our proposed work is the energy minimization. Active contour based segmentation technique entrusts on the internal mean (IM) and attraction of contour towards the object (AF) for curve evolution. On the other hand, segmentation through Level set does not entrust on the image gradient. The criterion responsible for curve evolution is formulated with the help of internal and external energy supportive parameters.

Implementation of active contour segmentation in our proposed work provides simple and suitable platform to detect errors and relatively adjust parameters. It is usually deployed to evolve a curve and is subjected to constraints in terms of discovering objects from the given input image. These constraints in active contour modeling are formulated by exploiting Level Set equations in our proposed work. For instance, initiating with the curve around the object to be distinguished, the curve progress towards its interior normal from the image based on few constraints. Once the boundary of the object is reached, the curve needs to be stopped. Traditional snakes and active contour models used edge detector in existing system to halt the curve evolving on the boundaries of the desired object. This parameter can identify only the objects with edge function exemplified by gradient.

Existing approach to stop the curve evolution is formed by two regions namely piecewise constant, intensity of distinct value (. This criterion produced multiple force term and multiple minimizer term. Moreover, initial formulated curve is not accurate, inorder to facilitate accuracy in curve evolution it needs to be adjusted under certain constraints. This in turn paves a way for multiple minimizer term. Our proposed work surmounts the aforementioned minimize problem by exploring one force term (AF) and one minimizer term (ME).

We framed energy minimization equation that accomplishes minimized energy to evolve the contour. The interior and exterior mean calculations drawn in equation (3) (4) and curvature information from the image will support the force of attraction computations.

(3)

(4)

Here, IM and EM are the interior and exterior mean; i = Interior Points; e = Exterior points; ε= 2.2204e-016; I = Image. This force of attraction is necessary to hold the contour in intact position. The force of attraction and curvature value is subjected to compute the energy to be minimized.

(5)

(6)

In equation (5) and (6), AF= Attractive Force; ME= Minimization Energy; α = 0.2; C= Curvature. In curvature dependency, when the differentiation is equal to the divergence in the curvature, minimized energy value is manipulated with the automatic generated mask to initialize the evolution of contour.

(6)

Where, CE= Curve Evolution; M= Mask Image; D= Experimental Constant. Mask image is obtained by acquiring the differentiation between repetitions of pixels (tumor pixels and tissue pixels) from the image. With the aid of above formulation, we solve the energy minimization problem. Evolving contour for segmentation along with acquired level set equation enhance the image segmentation as more accurate, regarding the tumor volume. Precisely, Contour based segmentation algorithm extracts the tumor slices of the lung image. Later, the extracted tumor portion pixels for the entire slices of the patient data set are reconstructed in 3D that makes the visualization of tumor portion more accurately.

3.5 3D Volume Reconstruction

When the segmentation on lung image is completed, then next stage is the 3D volumetric reconstruction of tumor slices from the segmented 2D CT images. Our work on reconstruction of 3D voxel relies on detection of isosurfaces. Identification of isosurfaces from volumetric data is the crucial task to obtain the volume reconstruction information. Extended marching cubes and dual contouring is the recently improved isosurface detection algorithm. To construct isosurface from scattered data sets, our approach builds on adaptive isosurface detection.

In this proposal, formulation of active contour and level set model intends to identify the isosurface and paves a way for 3D voxel reconstruction. Significant isosurfaces are substantial for provoking advisory visualizations of volume data. As we move to higher dimension data the level set contours become level surfaces. These level surfaces are referred as isosurfaces which provokes the initialization of patch data. Patch used on each pixel of the segmented tumor image determines the tumor 3D visualization. Initially, the volume data are subjected to the smoothening process by making use of default convolution kernels. As a result, smoothened volume data is generated. The computation of isosurface data is effectuated from the obtained smoothened volume data. The isosurface value is extracted from the isosurface data for the computation of tri-surface geometry. With the given isosurface value, tri-surface geometry of the volume data is estimated, which form the basis for patch assignment. The resultant output generation contains the faces and vertices of the isosurfaces.

Patches support a variety of coloring options that are useful for visualizing the data superimposed on geometric shapes. The segmented tumor slice is split up into multiple patches according to the faces and vertices of the isosurfaces. Since, the isosurface connects points that have the specified value much the way contour lines connect points of equal elevation. Finally, we compute the normal for the vertices generated. The 3D reconstructed model is then manipulated to compute the voxel impaired by tumor. The percentage of damaged voxel guides the surgical procedure. Finally, tumor volume calculation is percolated to estimate the number of abnormal voxels and percentage of tumor cells.

Representation of tumor size in 3D is termed as tumor volume, which facilitates excellent visualization about the pathological state. The 3D volume construction and visualization adduced in our proposed work allows exploring tumor slices as 3D model. The amount of extracted contour segmented image volume is used to construct the 3D model directly influences its visualization. In our study, greatest tumor volume increases the vulnerability of lung tumor in patient.

IV. PERFORMANCE EVALUATION In this section, validation framework is performed to evaluate the quality of 3D volumetric data reconstruction. As an application, we apply our method to construct a high resolution 3D model of the lung tumor from 2D slices of CT lung image.

To evaluate our proposed work, we conduct experiments on the entire slices of patient data set. Adduced work is trained with 150 slices (51 normal, 99 abnormal) for lung tumor diagnosis and segmentation. In this analysis, preprocessing and feature extraction is the first procedure that outputs the denoised image along with the subset of most essential features. Fig 2. a) signifies the input lung image taken for diagnosis, b) Enhanced lung image after denoising.

Fig 2: a) Input CT lung Image b) Enhanced Lung Image

The features extracted from denoised lung image are fed to classifier to classify the image as normal and abnormal categories. Performance of classification accuracy is predicted in terms of classification measures illustrated in table 1. True Positive/Negative and False Positive/Negative are the balance measures taken into account for classification.

Table 1: Representation of classification measures

Notice the classification measures shown in table 1, TN- True Negative, TP – True Positive, FN – False Negative, FP – False Positive. Moreover, confusion matrix affords knowledge related to actual and predicted cases generated by Multivariate Multinomial Distributed Bayes Classification. It is represented graphically by the obtained classification results. Confusion matrix is defined as the format of table 2.

Table 2: Confusion Matrix Representation

It is obvious from the table 2, TP (51) - Estimate Tumor as Tumor; FP (0) - Estimate Tumor as Normal; TN (98) - Estimate Normal as Normal; FN (1) - Estimate Normal as Tumor. The confusion matrix of our classification result is drawn on fig 3. Our classification system accomplishes the accuracy rate of 99%. It is evidently proved that classification technique exploited in this study act as a proficient classifier for lung image diagnosis, maximize the accuracy rate by distinguishing the lung image as tumor and normal.

Fig 3: Confusion Matrix

Once the diagnosis of image is accomplished through our method, next stage is the segmentation of tumor mass from abnormal images. Active contour and level set model deployed in this work enforce to identify the smooth contour towards the tumor mass. Formulation of level set equations integrated with active contour modeling processed in section 3.4 minimizes the energy by evolving the curve. The tumor mass extracted by the proposed contour based segmentation algorithm shown in fig 4 assists the physician in identifying the appropriate tumor size.

Fig 4: Tumor Segmentation

Our major contribution in this proposal intends in 3D volume reconstruction from the extracted tumor portion pixels. Significant isosurfaces substantial for provoking advisory visualizations of volume data are computed from the smoothened volume data as described in section 3.5.

Table 3: Sample Tumor Slices with Tumor Size in Pixels

The smoothened volume isosurfaces provokes the intialization of patch data as portrayed in fig 5.

Fig 5: Patch conversion of extracted tumor slices

Finally, we reconstruct the 3D tumor visualization (fig 6) and are manipulated to estimate the voxel impaired by the tumor.

Fig 6: Tumor 3D Visualization

Tumor volume calculation is incorporated in this experiment by evaluating the highest axis of the tumor along x, y and z orientation. Table 3 affirms the sample tumor size computed in pixels for the segmented tumor portion of image. Measurements like number of abnormal voxels and the percentage of tumor cells in tumor volume computation expose that the proposed system for 3D volume reconstruction works superior and enhance the tumor diagnosis. With the computed tumor volume visualized in 3D, it is compared with the diagnosis of a physician and then further analysis is conducted.

V. CONCLUSION AND PATH AHEAD

Numerous works have been devoted for 2D image diagnosis and segmentation. In this paper, we adduce a novel 3D reconstruction of volumetric data from 2D slices of a CT lung image. The main contribution of the approach involves lung tumor segmentation and its projection in 3D. A Contour based segmentation algorithm implemented in our proposal, concentrates on acquiring the exact dissimilarity of pixel boundary. With the aid of active contour and level set model, extracted tumor portion pixels are reconstructed. Proposed 3D volume reconstruction detects the isosurface under smoothening process by computing tri-surface geometry of a given volume data. The 3D reconstructed model is then manipulated to compute the voxel impaired by tumor. The volume reconstruction information is indispensable in the case of surgical operations. The 3D volume reconstruction intended in this paper surmounts the problems imposed by the 2D visualization of data.