Input Slices Research Articles

AcquisitionFocus: Joint Optimization of Acquisition Orientation and Cardiac Volume Reconstruction Using Deep Learning.

In cardiac cine imaging, acquiring high-quality data is challenging and time-consuming due to the artifacts generated by the heart's continuous movement. Volumetric, fully isotropic data acquisition with high temporal resolution is, to date, intractable due to MR physics constraints. To assess whole-heart movement under minimal acquisition time, we propose a deep learning model that reconstructs the volumetric shape of multiple cardiac chambers from a limited number of input slices while simultaneously optimizing the slice acquisition orientation for this task. We mimic the current clinical protocols for cardiac imaging and compare the shape reconstruction quality of standard clinical views and optimized views. In our experiments, we show that the jointly trained model achieves accurate high-resolution multi-chamber shape reconstruction with errors of <13 mm HD95 and Dice scores of >80%, indicating its effectiveness in both simulated cardiac cine MRI and clinical cardiac MRI with a wide range of pathological shape variations.

CeLNet: a correlation-enhanced lightweight network for medical image segmentation

Objective. Convolutional neural networks have been widely adopted for medical image segmentation with their outstanding feature representation capabilities. As the segmentation accuracy gets constantly updated, the complexity of networks increases as well. Complex networks can achieve better performance but require more parameters and are hard to train with limited resources, while lightweight models are faster but cannot fully utilize the contextual information of medical images. In this paper, we focus on better balancing the efficiency and accuracy. Approach. We propose a correlation-enhanced lightweight network (CeLNet) for medical image segmentation, which adopts a siamese structure for weight sharing and parameter saving. Through the feature reuse and feature stacking of parallel branches, a point-depth convolution parallel block (PDP Block) is proposed to reduce the model parameters and computational cost while improving the feature extraction capability of encoder. A relation module is also designed to extract feature correlations of input slices, which utilizes global and local attention to enhance feature connections, while reducing feature differences through element subtraction, and finally obtains contextual information of associated slices to improve the segmentation performance. Main results. We conduct extensive experiments on the LiTS2017, MM-WHS and ISIC2018 datasets, and the proposed model consumes merely 5.18M parameters but achieves excellent segmentation performance, specifically, a DSC of 0.9233 in LiTS2017 dataset, an average DSC of 0.7895 on MM-WHS dataset and an average DSC of 0.8401 on ISIC2018 dataset. Significance. CeLNet achieves state-of-the-art performance in multiple datasets while ensuring lightweight.

A multi-kernel and multi-scale learning based deep ensemble model for predicting recurrence of non-small cell lung cancer.

Predicting recurrence in patients with non-small cell lung cancer (NSCLC) before treatment is vital for guiding personalized medicine. Deep learning techniques have revolutionized the application of cancer informatics, including lung cancer time-to-event prediction. Most existing convolutional neural network (CNN) models are based on a single two-dimensional (2D) computational tomography (CT) image or three-dimensional (3D) CT volume. However, studies have shown that using multi-scale input and fusing multiple networks provide promising performance. This study proposes a deep learning-based ensemble network for recurrence prediction using a dataset of 530 patients with NSCLC. This network assembles 2D CNN models of various input slices, scales, and convolutional kernels, using a deep learning-based feature fusion model as an ensemble strategy. The proposed framework is uniquely designed to benefit from (i) multiple 2D in-plane slices to provide more information than a single central slice, (ii) multi-scale networks and multi-kernel networks to capture the local and peritumoral features, (iii) ensemble design to integrate features from various inputs and model architectures for final prediction. The ensemble of five 2D-CNN models, three slices, and two multi-kernel networks, using 5 × 5 and 6 × 6 convolutional kernels, achieved the best performance with an accuracy of 69.62%, area under the curve (AUC) of 72.5%, F1 score of 70.12%, and recall of 70.81%. Furthermore, the proposed method achieved competitive results compared with the 2D and 3D-CNN models for cancer outcome prediction in the benchmark studies. Our model is also a potential adjuvant treatment tool for identifying NSCLC patients with a high risk of recurrence.

DeepAD: A deep learning application for predicting amyloid standardized uptake value ratio through PET for Alzheimer's prognosis.

Amyloid deposition is a vital biomarker in the process of Alzheimer's diagnosis. 18F-florbetapir PET scans can provide valuable imaging data to determine cortical amyloid quantities. However, the process is labor and doctor intensive, requiring extremely specialized education and resources that may not be accessible to everyone, making the amyloid calculation process inefficient. Deep learning is a rising tool in Alzheimer's research which could be used to determine amyloid deposition. Using data from the Alzheimer's Disease Neuroimaging Initiative, we identified 2,980 patients with PET imaging, clinical, and genetic data. We tested various ResNet, EfficientNet, and RegNet convolutional neural networks and later combined the best performing model with Gradient Boosting Decision Tree algorithms to predict standardized uptake value ratio (SUVR) of amyloid in each patient session. We tried several configurations to find the best model tuning for regression-to-SUVR. We found that the RegNet X064 architecture combined with a grid search-tuned Gradient Boosting Decision Tree with 3 axial input slices and clinical and genetic data achieved the lowest loss. Using the mean-absolute-error metric, the loss converged to an MAE of 0.0441, equating to 96.4% accuracy across the 596-patient test set. We showed that this method is more consistent and accessible in comparison to human readers from previous studies, with lower margins of error and substantially faster calculation times. We implemented our deep learning model on to a web application named DeepAD which allows our diagnostic tool to be accessible. DeepAD could be used in hospitals and clinics with resource limitations for amyloid deposition and shows promise for more imaging tasks as well.

Multislice input for 2D and 3D residual convolutional neural network noise reduction in CT.

Deep convolutional neural network (CNN)-based methods are increasingly used for reducing image noise in computed tomography (CT). Current attempts at CNN denoising are based on 2D or 3D CNN models with a single- or multiple-slice input. Our work aims to investigate if the multiple-slice input improves the denoising performance compared with the single-slice input and if a 3D network architecture is better than a 2D version at utilizing the multislice input. Two categories of network architectures can be used for the multislice input. First, multislice images can be stacked channel-wise as the multichannel input to a 2D CNN model. Second, multislice images can be employed as the 3D volumetric input to a 3D CNN model, in which the 3D convolution layers are adopted. We make performance comparisons among 2D CNN models with one, three, and seven input slices and two versions of 3D CNN models with seven input slices and one or three output slices. Evaluation was performed on liver CT images using three quantitative metrics with full-dose images as reference. Visual assessment was made by an experienced radiologist. When the input channels of the 2D CNN model increases from one to three to seven, a trend of improved performance was observed. Comparing the three models with the seven-slice input, the 3D CNN model with a one-slice output outperforms the other models in terms of noise texture and homogeneity in liver parenchyma as well as subjective visualization of vessels. We conclude the that multislice input is an effective strategy for improving performance for 2D deep CNN denoising models. The pure 3D CNN model tends to have a better performance than the other models in terms of continuity across axial slices, but the difference was not significant compared with the 2D CNN model with the same number of slices as the input.

Phase recognition in contrast‐enhanced CT scans based on deep learning and random sampling

A fully automated system for interpreting abdominal computed tomography (CT) scans with multiple phases of contrast enhancement requires an accurate classification of the phases. Current approaches to classify the CT phases are commonly based on three-dimensional (3D) convolutional neural network (CNN) approaches with high computational complexity and high latency. This work aims at developing and validating a precise, fast multiphase classifier to recognize three main types of contrast phases in abdominal CTscans. We propose in this study a novel method that uses a random sampling mechanism on top of deep CNNs for the phase recognition of abdominal CT scans of four different phases: noncontrast, arterial, venous, and others. The CNNs work as a slicewise phase prediction, while random sampling selects input slices for the CNN models. Afterward, majority voting synthesizes the slicewise results of the CNNs to provide the final prediction at the scanlevel. Our classifier was trained on 271426 slices from 830 phase-annotated CT scans, and when combined with majority voting on 30% of slices randomly chosen from each scan, achieved a mean F1 score of 92.09% on our internal test set of 358 scans. The proposed method was also evaluated on two external test sets: CTPAC-CCRCC (N = 242) and LiTS (N = 131), which were annotated by our experts. Although a drop in performance was observed, the model performance remained at a high level of accuracy with a mean F1 scores of 76.79% and 86.94% on CTPAC-CCRCC and LiTS datasets, respectively. Our experimental results also showed that the proposed method significantly outperformed the state-of-the-art 3D approaches while requiring less computation time forinference. In comparison to state-of-the-art classification methods, the proposed approach shows better accuracy with significantly reduced latency. Our study demonstrates the potential of a precise, fast multiphase classifier based on a two-dimensional deep learning approach combined with a random sampling method for contrast phase recognition, providing a valuable tool for extracting multiphase abdomen studies from low veracity, real-worlddata.

An externally validated fully automated deep learning algorithm to classify COVID-19 and other pneumonias on chest computed tomography.

PurposeIn this study, we propose an artificial intelligence (AI) framework based on three-dimensional convolutional neural networks to classify computed tomography (CT) scans of patients with coronavirus disease 2019 (COVID-19), influenza/community-acquired pneumonia (CAP), and no infection, after automatic segmentation of the lungs and lung abnormalities.MethodsThe AI classification model is based on inflated three-dimensional Inception architecture and was trained and validated on retrospective data of CT images of 667 adult patients (no infection n=188, COVID-19 n=230, influenza/CAP n=249) and 210 adult patients (no infection n=70, COVID-19 n=70, influenza/CAP n=70), respectively. The model's performance was independently evaluated on an internal test set of 273 adult patients (no infection n=55, COVID-19 n= 94, influenza/CAP n=124) and an external validation set from a different centre (305 adult patients: COVID-19 n=169, no infection n=76, influenza/CAP n=60).ResultsThe model showed excellent performance in the external validation set with area under the curve of 0.90, 0.92 and 0.92 for COVID-19, influenza/CAP and no infection, respectively. The selection of the input slices based on automatic segmentation of the abnormalities in the lung reduces analysis time (56 s per scan) and computational burden of the model. The Transparent Reporting of a Multivariable Prediction Model for Individual Prognosis or Diagnosis (TRIPOD) score of the proposed model is 47% (15 out of 32 TRIPOD items).ConclusionThis AI solution provides rapid and accurate diagnosis in patients suspected of COVID-19 infection and influenza.

ECSU-Net: An Embedded Clustering Sliced U-Net Coupled With Fusing Strategy for Efficient Intervertebral Disc Segmentation and Classification.

Automatic vertebra segmentation from computed tomography (CT) image is the very first and a decisive stage in vertebra analysis for computer-based spinal diagnosis and therapy support system. However, automatic segmentation of vertebra remains challenging due to several reasons, including anatomic complexity of spine, unclear boundaries of the vertebrae associated with spongy and soft bones. Based on 2D U-Net, we have proposed an Embedded Clustering Sliced U-Net (ECSU-Net). ECSU-Net comprises of three modules named segmentation, intervertebral disc extraction (IDE) and fusion. The segmentation module follows an instance embedding clustering approach, where our three sliced sub-nets use axis of CT images to generate a coarse 2D segmentation along with embedding space with the same size of the input slices. Our IDE module is designed to classify vertebra and find the inter-space between two slices of segmented spine. Our fusion module takes the coarse segmentation (2D) and outputs the refined 3D results of vertebra. A novel adaptive discriminative loss (ADL) function is introduced to train the embedding space for clustering. In the fusion strategy, three modules are integrated via a learnable weight control component, which adaptively sets their contribution. We have evaluated classical and deep learning methods on Spineweb dataset-2. ECSU-Net has provided comparable performance to previous neural network based algorithms achieving the best segmentation dice score of 95.60% and classification accuracy of 96.20%, while taking less time and computation resources.

A myocardial extraction method using deep learning for 99mTc myocardial perfusion SPECT images: A basic study to reduce the effects of extra-myocardial activity

AimThe purpose of this study was to automatically extract myocardial regions from transaxial single-photon emission computed tomography (SPECT) images using deep learning to reduce the effects of extracardiac activity, which has been problematic in cardiac nuclear imaging. Method: Myocardial region extraction was performed using two deep neural network architectures, U-Net and U-Net ++, and 694 myocardial SPECT images manually labeled with myocardial regions were used as the training data. In addition, a multi-slice input method was introduced during the learning session while taking the relationships to adjacent slices into account. Accuracy was assessed using Dice coefficients at both the slice and pixel levels, and the most effective number of input slices was determined. Results: The Dice coefficient was 0.918 at the pixel level, and there were no false positives at the slice level using U-Net++ with 9 input slices. Conclusion: The proposed system based on U-Net++ with multi-slice input provided highly accurate myocardial region extraction and reduced the effects of extracardiac activity in myocardial SPECT images.

Evaluation of multislice inputs to convolutional neural networks for medical image segmentation.

When using convolutional neural networks (CNNs) for segmentation of organs and lesions in medical images, the conventional approach is to work with inputs and outputs either as single slice [two-dimensional (2D)] or whole volumes [three-dimensional (3D)]. One common alternative, in this study denoted as pseudo-3D, is to use a stack of adjacent slices as input and produce a prediction for at least the central slice. This approach gives the network the possibility to capture 3D spatial information, with only a minor additional computational cost. In this study, we systematically evaluate the segmentation performance and computational costs of this pseudo-3D approach as a function of the number of input slices, and compare the results to conventional end-to-end 2D and 3D CNNs, and to triplanar orthogonal 2D CNNs. The standard pseudo-3D method regards the neighboring slices as multiple input image channels. We additionally design and evaluate a novel, simple approach where the input stack is a volumetric input that is repeatably convolved in 3D to obtain a 2D feature map. This 2D map is in turn fed into a standard 2D network. We conducted experiments using two different CNN backbone architectures and on eight diverse data sets covering different anatomical regions, imaging modalities, and segmentation tasks. We found that while both pseudo-3D methods can process a large number of slices at once and still be computationally much more efficient than fully 3D CNNs, a significant improvement over a regular 2D CNN was only observed with two of the eight data sets. triplanar networks had the poorest performance of all the evaluated models. An analysis of the structural properties of the segmentation masks revealed no relations to the segmentation performance with respect to the number of input slices. A post hoc rank sum test which combined all metrics and data sets yielded that only our newly proposed pseudo-3D method with an input size of 13 slices outperformed almost all methods. In the general case, multislice inputs appear not to improve segmentation results over using 2D or 3D CNNs. For the particular case of 13 input slices, the proposed novel pseudo-3D method does appear to have a slight advantage across all data sets compared to all other methods evaluated in this work.

Application of Deep Learning to Predict Standardized Uptake Value Ratio and Amyloid Status on 18F-Florbetapir PET Using ADNI Data.

Cortical amyloid quantification on PET by using the standardized uptake value ratio is valuable for research studies and clinical trials in Alzheimer disease. However, it is resource intensive, requiring co-registered MR imaging data and specialized segmentation software. We investigated the use of deep learning to automatically quantify standardized uptake value ratio and used this for classification. Using the Alzheimer's Disease Neuroimaging Initiative dataset, we identified 2582 18F-florbetapir PET scans, which were separated into positive and negative cases by using a standardized uptake value ratio threshold of 1.1. We trained convolutional neural networks (ResNet-50 and ResNet-152) to predict standardized uptake value ratio and classify amyloid status. We assessed performance based on network depth, number of PET input slices, and use of ImageNet pretraining. We also assessed human performance with 3 readers in a subset of 100 randomly selected cases. We have found that 48% of cases were amyloid positive. The best performance was seen for ResNet-50 by using regression before classification, 3 input PET slices, and pretraining, with a standardized uptake value ratio root-mean-square error of 0.054, corresponding to 95.1% correct amyloid status prediction. Using more than 3 slices did not improve performance, but ImageNet initialization did. The best trained network was more accurate than humans (96% versus a mean of 88%, respectively). Deep learning algorithms can estimate standardized uptake value ratio and use this to classify 18F-florbetapir PET scans. Such methods have promise to automate this laborious calculation, enabling quantitative measurements rapidly and in settings without extensive image processing manpower and expertise.

Ensemble of deep convolutional neural networks based multi‐modality images for Alzheimer's disease diagnosis

Alzheimer's disease (AD) is one of the most common progressive neurodegenerative diseases. Structural magnetic resonance imaging (MRI) would provide abundant information on the anatomical structure of human organs. Fluorodeoxy-glucose positron emission tomography (PET) obtains the metabolic activity of the brain. Previous studies have demonstrated that multi-modality images could contribute to improve diagnosis of AD. However, these methods need to extract the handcrafted features that demand domain specific knowledge and image processing stage is time consuming. In order to tackle these problems, in this study, the authors propose a novel framework that ensembles three state-of-the-art deep convolutional neural networks (DCNNs) with multi-modality images for AD classification. In detail, they extract some slices from each subject of each modality, and every DCNN generates a probabilistic score for the input slices. Furthermore, a `dropout' mechanism is introduced to discard low discrimination slices of the category probabilities. Then average reserved slices of each subject are acquired as a new feature. Finally, they train the Adaboost ensemble classifier based on single decision tree classifier with the MRI and PET probabilistic scores of each DCNN. Evaluations on Alzheimer's Disease Neuroimaging Initiative database show that the proposed algorithm has better performance compared to existing method, the algorithm proposed in this study significantly improved the classification accuracy.

Brain Tumor Detection by Using Stacked Autoencoders in Deep Learning.

Brain tumor detection depicts a tough job because of its shape, size and appearance variations. In this manuscript, a deep learning model is deployed to predict input slices as a tumor (unhealthy)/non-tumor (healthy). This manuscript employs a high pass filter image to prominent the inhomogeneities field effect of the MR slices and fused with the input slices. Moreover, the median filter is applied to the fused slices. The resultant slices quality is improved with smoothen and highlighted edges of the input slices. After that, based on these slices' intensity, a 4-connected seed growing algorithm is applied, where optimal threshold clusters the similar pixels from the input slices. The segmented slices are then supplied to the fine-tuned two layers proposed stacked sparse autoencoder (SSAE) model. The hyperparameters of the model are selected after extensive experiments. At the first layer, 200 hidden units and at the second layer 400 hidden units are utilized. The testing is performed on the softmax layer for the prediction of the images having tumors and no tumors. The suggested model is trained and checked on BRATS datasets i.e., 2012(challenge and synthetic), 2013, and 2013 Leaderboard, 2014, and 2015 datasets. The presented model is evaluated with a number of performance metrics which demonstrates the improved performance.

Synthesis of CT images from digital body phantoms using CycleGAN.

The potential of medical image analysis with neural networks is limited by the restricted availability of extensive data sets. The incorporation of synthetic training data is one approach to bypass this shortcoming, as synthetic data offer accurate annotations and unlimited data size. We evaluated eleven CycleGAN for the synthesis of computed tomography (CT) images based on XCAT body phantoms. The image quality was assessed in terms of anatomical accuracy and realistic noise properties. We performed two studies exploring various network and training configurations as well as a task-based adaption of the corresponding loss function. The CycleGAN using the Res-Net architecture and three XCAT input slices achieved the best overall performance in the configuration study. In the task-based study, the anatomical accuracy of the generated synthetic CTs remained high ([Formula: see text] and [Formula: see text]). At the same time, the generated noise texture was close to real data with a noise power spectrum correlation coefficient of [Formula: see text]. Simultaneously, we observed an improvement in annotation accuracy of 65% when using the dedicated loss function. The feasibility of a combined training on both real and synthetic data was demonstrated in a blood vessel segmentation task (dice similarity coefficient [Formula: see text]). CT synthesis using CycleGAN is a feasible approach to generate realistic images from simulated XCAT phantoms. Synthetic CTs generated with a task-based loss function can be used in addition to real data to improve the performance of segmentation networks.

A model based segmentation approach for lung segmentation from chest computer tomography images

Segmentation of pathological lung regions from the body regions is the most challenging part in any Computer Aided Diagnosis (CAD) system due to the presence of Pathological Bearing Regions (PBR) which lies on the lung’s periphery. Here, a unique pathological lung segmentation method called reference-model based segmentation that uses shape property of human lung is proposed. This method trounces the difficulties in segmentation from traditional approaches by examining the shape knowledge of lung. The proposed segmentation approach constructs a reference lung model from input slices using a novel Sampling Lines Algorithm (SLA) and extracts the shape features. The segmentation work is validated using dataset consisting of Digital Imaging and Communications in Medicines (DICOM) standard chest CT images of seven patients from cancer institute Chennai, nine patients from Gemini Scans, Chennai and fourteen patients from Lung Image Data-base Consortium image collection (LIDC-IDRI). The segmentation method’s performance is analyzed against widely used segmentation methods namely Graph Cuts (GC), Region Growing (RG), Active Contour and Flood Fill in terms of accuracy, specificity, sensitivity, overlap score, Jaccard index, and also Dices similarity coefficients (DSC). The numerical outcomes specify that the proposed work attains an improved result against widely used segmentation techniques.

Brain tumor detection using statistical and machine learning method

Background and ObjectiveBrain tumor occurs because of anomalous development of cells. It is one of the major reasons of death in adults around the globe. Millions of deaths can be prevented through early detection of brain tumor. Earlier brain tumor detection using Magnetic Resonance Imaging (MRI) may increase patient's survival rate. In MRI, tumor is shown more clearly that helps in the process of further treatment. This work aims to detect tumor at an early phase. MethodsIn this manuscript, Weiner filter with different wavelet bands is used to de-noise and enhance the input slices. Subsets of tumor pixels are found with Potential Field (PF) clustering. Furthermore, global threshold and different mathematical morphology operations are used to isolate the tumor region in Fluid Attenuated Inversion Recovery (Flair) and T2 MRI. For accurate classification, Local Binary Pattern (LBP) and Gabor Wavelet Transform (GWT) features are fused. ResultsThe proposed approach is evaluated in terms of peak signal to noise ratio (PSNR), mean squared error (MSE) and structured similarity index (SSIM) yielding results as 76.38, 0.037 and 0.98 on T2 and 76.2, 0.039 and 0.98 on Flair respectively. The segmentation results have been evaluated based on pixels, individual features and fused features. At pixels level, the comparison of proposed approach is done with ground truth slices and also validated in terms of foreground (FG) pixels, background (BG) pixels, error region (ER) and pixel quality (Q). The approach achieved 0.93 FG and 0.98 BG precision and 0.010 ER on a local dataset. On multimodal brain tumor segmentation challenge dataset BRATS 2013, 0.93 FG and 0.99 BG precision and 0.005 ER are acquired. Similarly on BRATS 2015, 0.97 FG and 0.98 BG precision and 0.015 ER are obtained. In terms of quality, the average Q value and deviation are 0.88 and 0.017. At the fused feature based level, specificity, sensitivity, accuracy, area under the curve (AUC) and dice similarity coefficient (DSC) are 1.00, 0.92, 0.93, 0.96 and 0.96 on BRATS 2013, 0.90, 1.00, 0.97, 0.98 and 0.98 on BRATS 2015 and 0.90, 0.91, 0.90, 0.77 and 0.95 on local dataset respectively. ConclusionThe presented approach outperformed as compared to existing approaches.

Super-resolution musculoskeletal MRI using deep learning.

To develop a super-resolution technique using convolutional neural networks for generating thin-slice knee MR images from thicker input slices, and compare this method with alternative through-plane interpolation methods. We implemented a 3D convolutional neural network entitled DeepResolve to learn residual-based transformations between high-resolution thin-slice images and lower-resolution thick-slice images at the same center locations. DeepResolve was trained using 124 double echo in steady-state (DESS) data sets with 0.7-mm slice thickness and tested on 17 patients. Ground-truth images were compared with DeepResolve, clinically used tricubic interpolation, and Fourier interpolation methods, along with state-of-the-art single-image sparse-coding super-resolution. Comparisons were performed using structural similarity, peak SNR, and RMS error image quality metrics for a multitude of thin-slice downsampling factors. Two musculoskeletal radiologists ranked the 3 data sets and reviewed the diagnostic quality of the DeepResolve, tricubic interpolation, and ground-truth images for sharpness, contrast, artifacts, SNR, and overall diagnostic quality. Mann-Whitney U tests evaluated differences among the quantitative image metrics, reader scores, and rankings. Cohen's Kappa (κ) evaluated interreader reliability. DeepResolve had significantly better structural similarity, peak SNR, and RMS error than tricubic interpolation, Fourier interpolation, and sparse-coding super-resolution for all downsampling factors (p < .05, except 4 × and 8 × sparse-coding super-resolution downsampling factors). In the reader study, DeepResolve significantly outperformed (p < .01) tricubic interpolation in all image quality categories and overall image ranking. Both readers had substantial scoring agreement (κ = 0.73). DeepResolve was capable of resolving high-resolution thin-slice knee MRI from lower-resolution thicker slices, achieving superior quantitative and qualitative diagnostic performance to both conventionally used and state-of-the-art methods.

A new feature matching algorithm of CT slices

This paper introduces a novel feature matching algorithm using both Image interpolating method and edge detecting method to make feature matching more precise. The proposed feature matching is composed of four elements: inputting data, interpolating method, features extracting and matching. At first, the tasks in human anatomy and the method to input slices data are introduced. Secondly, we propose a novel slices matching using image interpolating method and edge detecting method. Experimental results demonstrate that the proposed approach owns the properties both quickly and exactly.

A Morphology-Based Approach for Interslice Interpolation of Anatomical Slices From Volumetric Images

This paper proposes a new morphology-based approach for the interslice interpolation of current transformer (CT) and MRI datasets composed of parallel slices. Our approach is object based and accepts as input data binary slices belonging to the same anatomical structure. Such slices may contain one or more regions, since topological changes between two adjacent slices may occur. Our approach handles explicitly interslice topology changes by decomposing a many-to-many correspondence into three fundamental cases: one-to-one, one-to-many, and zero-to-one correspondences. The proposed interpolation process is iterative. One iteration of this process computes a transition sequence between a pair of corresponding input slices, and selects the element located at equal distance from the input slices. This algorithmic design yields a gradual, smooth change of shape between the input slices. Therefore, the main contribution of our approach is its ability to interpolate between two anatomic shapes by creating a smooth, gradual change of shape, and without generating over-smoothed interpolated shapes.

3D visualization of a rock sample, Rotliegend sandstone, Southern Permian Basin: applications for core analysis and petrophysics

A methodology has been developed to visualize and interpret the texture and pore system of rock samples using a virtual reality application. This technique is illustrated using an aeolian sandstone core sample of the Rotliegendes Dethlingen Formation from the Southern Permian basin of northwest Germany. The methodology involves digitizing a series of densely spaced microfocus computer tomography (μCT) images with an appropriate sampling and processing the individual images, which are based on density contrasts, into a volume definition similar to seismic data. Once the volume file has been created the data can be visualized and analysed using virtual reality software. Visualization of scanned core samples has benefits with regards to both the qualitative interpretation of the rock as well as to the quantitative evaluation of the intergranular volumes, cements and pore space. Further advantages exist in terms of the effective communication and collaboration within relevant communities in exploration and production. A key requirement for this technology to be beneficial is calibration of the input slices with rock properties. In terms of the qualitative analysis of the rock, visualization offers the possibility to view slices of the rock samples in any direction. Consequently, bedding planes and fracture planes can be viewed orthogonally even if the physical sample has been cut oblique to these planes. Both slice and transparent volume views are possible, readily revealing the sedimentary and structural features of the rock. Provided that the scanning process has taken place with sufficient resolution, the views can separate grains from cement, matrix and pore space. In addition to the qualitative aspect of visualization, the 3D interpretation capabilities of the software are very valuable. Voxel-based technology allows measurement of grain sizes, cement, matrix and pore space, both in terms of volume as well as geometry. This provides important information in addition, and complementary, to conventional core analysis. In most cases the volume measurements can be expected to be more precise than traditional 2D methods, such as point counting, because they account for the three-dimensional anisotropy of the rock. Stereo-viewing allows pore space to be ‘experienced’ in 3D, offering a better perception and analysis with a view towards hydrocarbon recovery and production strategies. Finally, visualization offers an optimal means to communicate and collaborate, creating a better shared understanding of the relevance of the rock properties within the overall exploration and production evaluation process. Supplementary material: The movie referred to in the article is available at https://doi.org/10.6084/m9.figshare.c.4812969