Search

4D cone-beam computed tomography (CBCT) is an important imaging modality in image-guided radiation therapy to address the motion-induced artifacts caused by organ movements during the respiratory process. However, due to the extremely sparse projection data for each temporal phase, 4D CBCT reconstructions will suffer from severe streaking artifacts. Therefore, to tackle the streak artifacts and provide high-quality images, we proposed a framework termed Prior-Regularized Iterative Optimization Reconstruction (PRIOR) for 4D CBCT. The PRIOR framework combines the physics-based model and data-driven method simultaneously, with powerful feature extracting capacity, significantly promoting the image quality compared to single model-based or deep learning-based methods. Besides, we designed a specialized deep learning model named PRIOR-Net, which can effectively excavate the static information in the prior image reconstructed from the fully-sampled projections at the encoding stage to improve the reconstruction performance for individual phase-resolved images. Both the simulated and clinical 4D CBCT datasets were performed to evaluate the performance of the PRIOR-Net and the PRIOR framework. Compared with the advanced 4D CBCT reconstruction methods, the proposed methods achieve promising results quantitatively and qualitatively in streak artifact suppression, soft tissue restoration, and tiny detail preservation. [ABSTRACT FROM AUTHOR]

The deep learning-based automatic recognition of the scanning or exposing region in medical imaging automation is a promising new technique, which can decrease the heavy workload of the radiographers, optimize imaging workflow and improve image quality. However, there is little related research and practice in X-ray imaging. In this paper, we focus on two key problems in X-ray imaging automation: automatic recognition of the exposure moment and the exposure region. Consequently, we propose an automatic video analysis framework based on the hybrid model, approaching real-time performance. The framework consists of three interdependent components: Body Structure Detection, Motion State Tracing, and Body Modeling. Body Structure Detection disassembles the patient to obtain the corresponding body keypoints and body Bboxes. Combining and analyzing the two different types of body structure representations is to obtain rich spatial location information about the patient body structure. Motion State Tracing focuses on the motion state analysis of the exposure region to recognize the appropriate exposure moment. The exposure region is calculated by Body Modeling when the exposure moment appears. A large-scale dataset for X-ray examination scene is built to validate the performance of the proposed method. Extensive experiments demonstrate the superiority of the proposed method in automatically recognizing the exposure moment and exposure region. This paradigm provides the first method that can enable automatically and accurately recognize the exposure region in X-ray imaging without the help of the radiographer. [ABSTRACT FROM AUTHOR]

Accurate and robust cephalometric image analysis plays an essential role in orthodontic diagnosis, treatment assessment and surgical planning. This paper proposes a novel landmark localization method for cephalometric analysis using multiscale image patch-based graph convolutional networks. In detail, image patches with the same size are hierarchically sampled from the Gaussian pyramid to well preserve multiscale context information. We combine local appearance and shape information into spatialized features with an attention module to enrich node representations in graph. The spatial relationships of landmarks are built with the incorporation of three-layer graph convolutional networks, and multiple landmarks are simultaneously updated and moved toward the targets in a cascaded coarse-to-fine process. Quantitative results obtained on publicly available cephalometric X-ray images have exhibited superior performance compared with other state-of-the-art methods in terms of mean radial error and successful detection rate within various precision ranges. Our approach performs significantly better especially in the clinically accepted range of 2 mm and this makes it suitable in cephalometric analysis and orthognathic surgery. [ABSTRACT FROM AUTHOR]

The present work describes a Bayesian maximum a posteriori (MAP) method using a statistical multiscale wavelet prior model. Rather than using the orthogonal discrete wavelet transform (DWT), this prior is built on the translation invariant wavelet transform (TIWT). The statistical modeling of wavelet coefficients relies on the generalized Gaussian distribution. Image reconstruction is performed in spatial domain with a fast block sequential iteration algorithm. We study theoretically the TIWT MAP method by analyzing the Hessian of the prior function to provide some insights on noise and resolution properties of image reconstruction. We adapt the key concept of local shift invariance and explore how the TIWT MAP algorithm behaves with different scales. It is also shown that larger support wavelet filters do not offer better performance in contrast recovery studies. These theoretical developments are confirmed through simulation studies. The results show that the proposed method is more attractive than other MAP methods using either the conventional Gibbs prior or the DWT-based wavelet prior. Index Terms--Emission tomography, expectation-maximization, image reconstruction, maximum a posteriori, translation invariant, wavelet transform.

In this paper, we present a PET reconstruction method using the wavelet-based maximum a posteriori (MAP) expectation-maximization (EM) algorithm. The proposed method, namely WV-MAP-EM, shows several advantages over conventional methods. It provides an adaptive way for hyperparameter determination. Since the wavelet transform allows the use of fast algorithms, WV-MAP-EM also does not increase the order of computational complexity. The spatial noise behavior (bias/variance and resolution) of the proposed MAP estimator is analyzed. Quantitative comparisons to MAP methods with Markov random field (MRF) prior models point out that our alternative method, wavelet-base method, offers competitive performance in PET image reconstruction. Index Terms--Expectation-maximization (EM), image reconstruction, maximum a posteriori (MAP), positron emission tomography (PET), wavelet transform.

This paper reports on a method for left ventricle three-dimensional (3-D) reconstruction from two orthogonal ventriculograms. The proposed algorithm is voxel-based and takes into account the conical projection geometry associated with the biplane image acquisition equipment. The reconstruction process starts with an initial ellipsoidal approximation derived from the input ventriculograms. This model is subsequently deformed in such a way as to match the input projections. To this end, the object is modeled as a 3-D Markov-Gibbs random field, and an energy function is defined so that it includes one term that models the projections compatibility and another one that includes the space-time regularity constraints. The performance of this reconstruction method is evaluated by considering the reconstruction of mathematically synthesized phantoms and two 3-D binary databases from two orthogonal synthesized projections. The method is also tested using real biplane ventriculograms. In this case, the performance of the reconstruction is expressed in terms of the projection error, which attains values between 9.50% and 11.78% for two biplane sequences including a total of 55 images. Index Terms--Markov random fields, reconstruction from projections, simulated annealing, three-dimensional (3-D) reconstruction, ventriculographic images.

Clinical interventional hemodynamic studies quantify the ventricular function from two-dimensional (2-D) X-ray projection images without having enough information of the actual three-dimensional (3-D) shape of this cardiac cavity. This paper reports a left ventricle 3-D reconstruction method from two orthogonal angiographic projections. This investigation is motivated by the lack of information about the actual 3-D shape of the cardiac cavity. The proposed algorithm works in 3-D space and considers the oblique projection geometry associated with the biplane image acquisition equipment. The reconstruction process starts by performing an approximate reconstruction based on the Cylindrical Closure Operation and the Dempster-Shafer theory. This approximate reconstruction is appropriately deformed in order to match the given projections. The deformation procedure is carried out by an iterative process that, by means of the Dempster--Shafer and the fuzzy integral theory, combines the information provided by the projection error and the connectivity between voxels. The performance of the proposed reconstruction method is evaluated by considering first the reconstruction of two 3-D binary databases from two orthogonal synthetized projections, obtaining errors as low as 6.48%. The method is then tested on real data, where two orthogonal preprocessed angiographic images are used for reconstruction. The performance of the technique, in this case, is assessed by means of the projection error, whose average for both views is 7.5%. The reconstruction method is also tested by performing the 3-D reconstruction of a ventriculographic sequence throughout an entire cardiac cycle. Index Terms--Angiographic images, evidential combination, fuzzy sets, left ventricle, three-dimensional (3-D) reconstruction.

In this paper, a model-based approach to medical image analysis is presented. It is aimed at understanding the influence of the physiological (related to tissue) and physical (related to image modality) processes underlying the image content. This methodology is exemplified by modeling first, the liver and its vascular network, and second, the standard computed tomography (CT) scan acquisition. After a brief survey on vascular modeling literature, a new method, aimed at the generation of growing three-dimensional vascular structures perfusing the tissue, is described. A solution is proposed in order to avoid intersections among vessels belonging to arterial and/or venous trees, which are physiologically connected. Then it is shown how the propagation of contrast material leads to simulate time-dependent sequences of enhanced liver CT slices. Index Terms--Computed tomography (CT), hepatic enhancement, image simulation, physiologically based modeling, vascular network.

Abstract Background The SYMBIOmatics Specific Support Action (SSA) is "an information gathering and dissemination activity" that seeks "to identify synergies between the bioinformatics and the medical informatics" domain to improve collaborative progress between both domains (ref. to http://www.symbiomatics.org). As part of the project experts in both research fields will be identified and approached through a survey. To provide input to the survey, the scientific literature was analysed to extract topics relevant to both medical informatics and bioinformatics. Results This paper presents results of a systematic analysis of the scientific literature from medical informatics research and bioinformatics research. In the analysis pairs of words (bigrams) from the leading bioinformatics and medical informatics journals have been used as indication of existing and emerging technologies and topics over the period 2000–2005 ("recent") and 1990–1990 ("past"). We identified emerging topics that were equally important to bioinformatics and medical informatics in recent years such as microarray experiments, ontologies, open source, text mining and support vector machines. Emerging topics that evolved only in bioinformatics were system biology, protein interaction networks and statistical methods for microarray analyses, whereas emerging topics in medical informatics were grid technology and tissue microarrays. Conclusion We conclude that although both fields have their own specific domains of interest, they share common technological developments that tend to be initiated by new developments in biotechnology and computer science.

In the field of epilepsy, the analysis of stereoelectroencephalographic (SEEG) signals recorded with depth electrodes provides major information on interactions between brain structures during seizures. A comprehensive methodology of comparing SEEG seizure recordings is presented. It proceeds in three steps: 1) segmentation of SEEG signals; 2) characterization and labeling of segments; and 3) comparison of observations coded as sequences of symbol vectors. The third step reports a vectorial extension of the Wagner and Fischer's algorithm to first, quantify similarities between observations and second, extract invariant sequences of events, referred to as spatio-temporal signatures. The study shows that two observations of nonequal duration can be matched by deforming the first one to optimally fit the second, under cost constraints. Results show that the methodology allows to exhibit signatures occurring during epileptic seizures and to point out different types of seizure patterns. The study brings objective results on reproducible interactions between brain structures during ictal periods and may help in the understanding of epileptogenic networks.

This paper presents a three-dimensional edge operator aimed at the detection of anatomical structures in medical imaging. It uses the spatial moments of gray level surface, and operates in three dimensions with any window size. It allows us to estimate the location and the contrast surface, as well as the surface orientation. The computation of the discrete version is reported. Bias and errors due to the spatial sampling and noise are analyzed both at a theoretical and experimental level. The moment-based operator is compared with other well-known edge operators on simple shaped primitives for which the analytical solution is known. The 3D rendering of real data is then provided by merging the operator in a ray-tracing framework.

Accurate segmentation of uterus, uterine fibroids, and spine from MR images is crucial for high intensity focused ultrasound (HIFU) therapy but remains still difficult to achieve because of 1) the large shape and size variations among individuals, 2) the low contrast between adjacent organs and tissues, and 3) the unknown number of uterine fibroids. To tackle this problem, in this paper, we propose a large kernel Encoder-Decoder Network based on a 2D segmentation model. The use of this large kernel can capture multi-scale contexts by enlarging the valid receptive field. In addition, a deep multiple atrous convolution block is also employed to enlarge the receptive field and extract denser feature maps. Our approach is compared to both conventional and other deep learning methods and the experimental results conducted on a large dataset show its effectiveness. [ABSTRACT FROM AUTHOR]

The wide applications of X-ray computed tomography (CT) bring low-dose CT (LDCT) into a clinical prerequisite, but reducing the radiation exposure in CT often leads to significantly increased noise and artifacts, which might lower the judgment accuracy of radiologists. In this paper, we put forward a domain progressive 3D residual convolution network (DP-ResNet) for the LDCT imaging procedure that contains three stages: sinogram domain network (SD-net), filtered back projection (FBP), and image domain network (ID-net). Though both are based on the residual network structure, the SD-net and ID-net provide complementary effect on improving the final LDCT quality. The experimental results with both simulated and real projection data show that this domain progressive deep-learning network achieves significantly improved performance by combing the network processing in the two domains. [ABSTRACT FROM AUTHOR]

An improved and very fast algorithm dealing with the extraction of vessels in three-dimensional imaging is described. The approach is based on geometrical moments and a local cylindrical approximation. A robust estimation of vessel and background intensity levels, position, orientation, and diameter of the vessels with adaptive control of key parameters, is provided during vessel tracking. Experimental results are presented for lower limb arteries in multidetector computed tomography scanner. Index Terms--CTA, geometrical moment, three-dimensional segmentation.

This work presents a novel method for multimodal medical registration based on histogram estimation of continuous image representation. The proposed method, regarded as “fast continuous histogram estimation,” employs continuous image representation to estimate the joint histogram of two images to be registered. The Jensen–Arimoto (JA) divergence is a similarity measure to measure the statistical dependence between medical images from different modalities. The estimated joint histogram is exploited to calculate the JA divergence in multimodal medical image registration. In addition, to reduce the grid effect caused by the grid-aligning transformations between two images and improve the implementation speed of the registration method, random samples instead of all pixels are extracted from the images to be registered. The goal of the registration is to optimize the JA divergence, which would be maximal when two misregistered images are perfectly aligned using the downhill simplex method, and thus to get the optimal geometric transformation. Experiments are conducted on an affine registration of 2D and 3D medical images. Results demonstrate the superior performance of the proposed method compared to standard histogram, Parzen window estimations, particle filter, and histogram estimation based on continuous image representation without random sampling. [ABSTRACT FROM AUTHOR]

Noise detection accuracy is crucial in suppressing random-valued impulse noise. Both false and miss detections determine the final estimation performance. Deterministic detection methods, which distinctly classify pixels into noisy or uncorrupted pixels, tend to increase the estimation error because some uncorrupted edge points are hard to discriminate from the random-valued impulse noise points. This paper proposes an iterative structure-adaptive fuzzy estimation (SAFE) for random-valued impulse noise suppression. This SAFE method is developed in the framework of Gaussian maximum likelihood estimation. The structure-adaptive fuzziness is reflected by two structure-adaptive metrics based on pixel reliability and patch similarity, respectively. The reliability metric for each pixel (as noise free) is estimated via a novel-minimal-path-based structure propagation to give full consideration of the spatially varying image structures. A robust iteration stopping strategy is also proposed by evaluating the reestimation error of the uncorrupted intensity information. The comparative experimental results show that the proposed structure-adaptive fuzziness can lead to effective restoration. An efficient implementation of this SAFE method is also realized via graphics-processing-unit-based parallelization. [ABSTRACT FROM AUTHOR]

The death of biomedical engineering advocate Dr. Pierre Galletti is announced. He was instrumental for bringing recognition to biomedical engineering from the government and the academe. He received his M.D. degree from the University of Lausanne in 1951 and went on to become the first chairman of the Division of Biological and Medical Sciences at Brown University. At the time of his death, Dr. Galletti was a member of the editorial boards of four research journals on biomedical engineering. He is survived by his wife Sonia and son Marc Henri.

Although the nonlocal means (NLM) algorithm takes a significant step forward in image filtering field, it suffers from a high computational complexity. To deal with this drawback, this paper proposes an acceleration strategy based on a correlation operation. Instead of per-pixel processing, this approach performs a simultaneous calculation of all the image pixels with the help of correlation operators. Complexity analysis and experimental results are reported and show the advantage of the proposed algorithm in terms of computation and time cost. [ABSTRACT FROM AUTHOR]

Purpose: The amount of coronary artery calcification (CAC) is a strong and independent predictor of cardiovascular disease (CVD) events. In clinical practice, CAC is manually identified and automatically quantified in cardiac CT using commercially available software. This is a tedious and time-consuming process in large-scale studies. Therefore, a number of automatic methods that require no interaction and semiautomatic methods that require very limited interaction for the identification of CAC in cardiac CT have been proposed. Thus far, a comparison of their performance has been lacking. The objective of this study was to perform an independent evaluation of (semi)automatic methods for CAC scoring in cardiac CT using a publicly available standardized framework. Methods: Cardiac CT exams of 72 patients distributed over four CVD risk categories were provided for (semi)automatic CAC scoring. Each exam consisted of a noncontrast-enhanced calcium scoring CT (CSCT) and a corresponding coronary CT angiography (CCTA) scan. The exams were acquired in four different hospitals using state-of-the-art equipment from four major CT scanner vendors. The data were divided into 32 training exams and 40 test exams. A reference standard for CAC in CSCT was defined by consensus of two experts following a clinical protocol. The framework organizers evaluated the performance of (semi)automatic methods on test CSCT scans, per lesion, artery, and patient. Results: Five (semi)automatic methods were evaluated. Four methods used both CSCT and CCTA to identify CAC, and one method used only CSCT. The evaluated methods correctly detected between 52% and 94% of CAC lesions with positive predictive values between 65% and 96%. Lesions in distal coronary arteries were most commonly missed and aortic calcifications close to the coronary ostia were the most common false positive errors. The majority (between 88% and 98%) of correctly identified CAC lesions were assigned to the correct artery. Linearly weighted Cohen's kappa for patient CVD risk categorization by the evaluated methods ranged from 0.80 to 1.00. Conclusions: A publicly available standardized framework for the evaluation of (semi)automatic methods for CAC identification in cardiac CT is described. An evaluation of five (semi)automatic methods within this framework shows that automatic per patient CVD risk categorization is feasible. CAC lesions at ambiguous locations such as the coronary ostia remain challenging, but their detection had limited impact on CVD risk determination. [ABSTRACT FROM AUTHOR]

Purpose: Calcium scoring is widely used to assess the risk of coronary heart disease (CHD). Accurate coronary artery calcification detection in noncontrast CT image is a prerequisite step for coronary calcium scoring. Currently, calcified lesions in the coronary arteries are manually identified by radiologists in clinical practice. Thus, in this paper, a fully automatic calcium scoring method was developed to alleviate the work load of the radiologists or cardiologists. Methods: The challenge of automatic coronary calcification detection is to discriminate the calcification in the coronary arteries from the calcification in the other tissues. Since the anatomy of coronary arteries is difficult to be observed in the noncontrast CT images, the contrast CT image of the same patient is used to extract the regions of the aorta, heart, and coronary arteries. Then, a patient-specific region-of-interest (ROI) is generated in the noncontrast CT image according to the segmentation results in the contrast CT image. This patient-specific ROI focuses on the regions in the neighborhood of coronary arteries for calcification detection, which can eliminate the calcifications in the surrounding tissues. A support vector machine classifier is applied finally to refine the results by removing possible image noise. Furthermore, the calcified lesions in the noncontrast images belonging to the different main coronary arteries are identified automatically using the labeling results of the extracted coronary arteries. Results: Forty datasets from four different CT machine vendors were used to evaluate their algorithm, which were provided by the MICCAI 2014 Coronary Calcium Scoring (orCaScore) Challenge. The sensitivity and positive predictive value for the volume of detected calcifications are 0.989 and 0.948. Only one patient out of 40 patients had been assigned to the wrong risk category defined according to Agatston scores (0, 1-100, 101-300, >300) by comparing with the ground truth. Conclusions: The calcified lesions in the noncontrast CT images can be detected automatically by using the segmentation results of the aorta, heart, and coronary arteries obtained in the contrast CT images with a very high accuracy. [ABSTRACT FROM AUTHOR]

Computational load remains a major concern when processing signals by means of sliding transforms. In this paper, we present an efficient algorithm for the fast computation of one-dimensional and two-dimensional sliding discrete Tchebichef moments. To do so, we first establish the relationships that exist between the Tchebichef moments of two neighboring windows taking advantage of Tchebichef polynomials’ properties. We then propose an original way to fast compute the moments of one window by utilizing the moment values of its previous window. We further theoretically establish the complexity of our fast algorithm and illustrate its interest within the framework of digital forensics and more precisely the detection of duplicated regions in an audio signal or an image. Our algorithm is used to extract local features of such a signal tampering. Experimental results show that its complexity is independent of the window size, validating the theory. They also exhibit that our algorithm is suitable to digital forensics and beyond to any applications based on sliding Tchebichef moments. [ABSTRACT FROM PUBLISHER]