Statistical construction of a Japanese male liver phantom for internal radionuclide dosimetry

Contour-based lung shape analysis in order to tuberculosis detection: modeling and feature description

Statistical shape analysis of lung is a reliable alternative method for diagnosing pulmonary diseases such as tuberculosis (TB). The 2D contour-based lung shape analysis is investigated and developed using Fourier descriptors (FDs). The proposed 2D lung shape analysis is carried out in threefold: (1) represent the normal and the abnormal (i.e. pulmonary tuberculosis (PTB)) lung shape models using Fourier descriptors modeling (FDM) framework from chest X-ray (CXR) images, (2) estimate and compare the 2D inter-patient lung shape variations for the normal and abnormal lungs by applying principal component analysis (PCA) techniques, and (3) describe the optimal type of contour-based feature vectors to train a classifier in order to detect TB using one publicly available dataset-namely the Montgomery dataset. Since almost all of the previous works in lung shape analysis are content-based analysis, we proposed contour-based lung shape analysis for statistical modeling and feature description of PTB cases. The results show that the proposed approach is able to explain more than 95% of total variations in both of the normal and PTB cases using only 6 and 7 principal component modes for the right and the left lungs, respectively. In case of PTB detection, using 138 lung cases (80 normal and 58 PTB cases), we achieved the accuracy (ACC) and the area under the curve (AUC) of 82.03% and 88.75%, respectively. In comparison with existing state-of-art studies in the same dataset, the proposed approach is a very promising supplement for diagnosis of PTB disease. The method is robust and valuable for application in 2D automatic segmentation, classification, and atlas registration. Moreover, the approach could be used for any kind of pulmonary diseases. Graphical abstract Contour-based lung shape analysis in order to detect tuberculosis: modeling and feature description.

Innovations in Computer Technologies Have Impacted Radiation Dosimetry Through Anatomically Realistic Phantoms and Fast Monte Carlo Simulations

Radiological physics principles have not changed in the past 60 y when computer technologies advanced exponentially. The research field of anatomical modeling for the purpose of radiation dose calculations has experienced an explosion in activity in the past two decades. Such an exciting advancement is due to the feasibility of creating three-dimensional geometric details of the human anatomy from tomographic imaging and of performing Monte Carlo radiation transport simulations on increasingly fast and cheap personal computers. The advent of a new type of high-performance computing hardware in recent years-graphics processing units-has made it feasible to carry out time-consuming Monte Carlo calculations at near real-time speeds. This paper introduces the history of three generations of computational human phantoms (the stylized medical internal radiation dosimetry-type phantoms, the voxelized tomographic phantoms, and the boundary representation deformable phantoms) and new development of the graphics processing unit-based Monte Carlo radiation dose calculations. Examples are given for research projects performed by my students in applying computational phantoms and a new Monte Carlo code, ARCHER, to problems in radiation protection, imaging, and radiotherapy. Finally, the paper discusses challenges and future opportunities for research.

Inter-Patient Modelling of 2D Lung Variations from Chest X-Ray Imaging via Fourier Descriptors

Detailed knowledge of anatomical lung variation is very important in medical image processing. Normal variations of lung consistent with the maintenance of pulmonary health and abnormal lung variations can be as a result of a pulmonary disease. Inter-patient lung variations can be due to the several factors such as sex, age, height, weight and type of disease. This study tries to show the inter-patient lung variations by using one of the shape-based descriptions techniques which is called Fourier descriptors. Shape-based description is an important approach to construct an object according to its parametric values. A different types of techniques are reported in the literature that aim to represent objects based on their shapes; each of these techniques has its cons and pros. Fourier descriptors, a simple yet powerful technique, has interesting properties such as rotational, scale, and translational invariance and these are powerful features for the recognition of two-dimensional connected shapes. In this paper, we use 380 CXR (Chest X-ray) images as a training set to construct the statistical mean model of lung contour. For modelling, the first step is evaluation of lung contour approximation and characterization to get the good spatial and frequency resolution. In the second step, all of the lung contours registered to show the variation and make a mean shape (i.e. lungs). And the final step is calculating the dispersion (i.e. covariance matrix) and analyzing by principle components. The proposed technique used to create the inter-patient statistical model and provide statistical parameters for application in segmentation, classification, 2D atlas based registration, etc. In this paper, we presented an approach for creating 2D modelling of human lungs from CXR image archives and reported some interesting statistical parameters to analysis the left and the right lung shape.

Deformable torso phantoms of Chinese adults for personalized anatomy modelling

In recent years, there has been increasing demand for personalized anatomy modelling for medical and industrial applications, such as ergonomics device development, clinical radiological exposure simulation, biomechanics analysis, and 3D animation character design. In this study, we constructed deformable torso phantoms that can be deformed to match the personal anatomy of Chinese male and female adults. The phantoms were created based on a training set of 79 trunk computed tomography (CT) images (41 males and 38 females) from normal Chinese subjects. Major torso organs were segmented from the CT images, and the statistical shape model (SSM) approach was used to learn the inter-subject anatomical variations. To match the personal anatomy, the phantoms were registered to individual body surface scans or medical images using the active shape model method. The constructed SSM demonstrated anatomical variations in body height, fat quantity, respiratory status, organ geometry, male muscle size, and female breast size. The masses of the deformed phantom organs were consistent with Chinese population organ mass ranges. To validate the performance of personal anatomy modelling, the phantoms were registered to the body surface scan and CT images. The registration accuracy measured from 22 test CT images showed a median Dice coefficient over 0.85, a median volume recovery coefficient (RCvlm ) between 0.85 and 1.1, and a median averaged surface distance (ASD) < 1.5 mm. We hope these phantoms can serve as computational tools for personalized anatomy modelling for the research community.

A feasibility study on the use of phantoms with statistical lung masses for determining the uncertainty in the dose absorbed by the lung from broad beams of incident photons and neutrons

Computational models of the human body have gradually become crucial in the evaluation of doses absorbed by organs. However, individuals may differ considerably in terms of organ size and shape. In this study, the authors sought to determine the energy-dependent standard deviations due to lung size of the dose absorbed by the lung during external photon and neutron beam exposures. One hundred lungs with different masses were prepared and located in an adult male International Commission on Radiological Protection (ICRP) reference phantom. Calculations were performed using the Monte Carlo N-particle code version 5 (MCNP5). Variation in the lung mass caused great uncertainty: ~90% for low-energy broad parallel photon beams. However, for high-energy photons, the lung-absorbed dose dependency on the anatomical variation was reduced to <1%. In addition, the results obtained indicated that the discrepancy in the lung-absorbed dose varied from 0.6% to 8% for neutron beam exposure. Consequently, the relationship between absorbed dose and organ volume was found to be significant for low-energy photon sources, whereas for higher energy photon sources the organ-absorbed dose was independent of the organ volume. In the case of neutron beam exposure, the maximum discrepancy (of 8%) occurred in the energy range between 0.1 and 5 MeV.

Monte Carlo simulation of a whole-body counter using IGOR phantoms

Whole-body counting is a technique of choice for assessing the intake of gamma-emitting radionuclides. An appropriate calibration is necessary, which is done either by experimental measurement or by Monte Carlo (MC) calculation. The aim of this work was to validate a MC model for calibrating whole-body counters (WBCs) by comparing the results of computations with measurements performed on an anthropomorphic phantom and to investigate the effect of a change in phantom's position on the WBC counting sensitivity. GEANT MC code was used for the calculations, and an IGOR phantom loaded with several types of radionuclides was used for the experimental measurements. The results show a reasonable agreement between measurements and MC computation. A 1-cm error in phantom positioning changes the activity estimation by >2%. Considering that a 5-cm deviation of the positioning of the phantom may occur in a realistic counting scenario, this implies that the uncertainty of the activity measured by a WBC is ∼10-20%.

An exponential growth of computational phantom research in radiation protection, imaging, and radiotherapy: a review of the fifty-year history

Radiation dose calculation using models of the human anatomy has been a subject of great interest to radiation protection, medical imaging, and radiotherapy. However, early pioneers of this field did not foresee the exponential growth of research activity as observed today. This review article walks the reader through the history of the research and development in this field of study which started some 50 years ago. This review identifies a clear progression of computational phantom complexity which can be denoted by three distinct generations. The first generation of stylized phantoms, representing a grouping of less than dozen models, was initially developed in the 1960s at Oak Ridge National Laboratory to calculate internal doses from nuclear medicine procedures. Despite their anatomical simplicity, these computational phantoms were the best tools available at the time for internal/external dosimetry, image evaluation, and treatment dose evaluations. A second generation of a large number of voxelized phantoms arose rapidly in the late 1980s as a result of the increased availability of tomographic medical imaging and computers. Surprisingly, the last decade saw the emergence of the third generation of phantoms which are based on advanced geometries called boundary representation (BREP) in the form of Non-Uniform Rational B-Splines (NURBS) or polygonal meshes. This new class of phantoms now consists of over 287 models including those used for non-ionizing radiation applications. This review article aims to provide the reader with a general understanding of how the field of computational phantoms came about and the technical challenges it faced at different times. This goal is achieved by defining basic geometry modeling techniques and by analyzing selected phantoms in terms of geometrical features and dosimetric problems to be solved. The rich historical information is summarized in four tables that are aided by highlights in the text on how some of the most well-known phantoms were developed and used in practice. Some of the information covered in this review has not been previously reported, for example, the CAM and CAF phantoms developed in 1970s for space radiation applications. The author also clarifies confusion about 'population-average' prospective dosimetry needed for radiological protection under the current ICRP radiation protection system and 'individualized' retrospective dosimetry often performed for medical physics studies. To illustrate the impact of computational phantoms, a section of this article is devoted to examples from the author's own research group. Finally the author explains an unexpected finding during the course of preparing for this article that the phantoms from the past 50 years followed a pattern of exponential growth. The review ends on a brief discussion of future research needs (a supplementary file '3DPhantoms.pdf' to figure 15 is available for download that will allow a reader to interactively visualize the phantoms in 3D).

Classification of normal and diseased liver shapes based on Spherical Harmonics coefficients

Liver-shape analysis and quantification is still an open research subject. Quantitative assessment of the liver is of clinical importance in various procedures such as diagnosis, treatment planning, and monitoring. Liver-shape classification is of clinical importance for corresponding intra-subject and inter-subject studies. In this research, we propose a novel technique for the liver-shape classification based on Spherical Harmonics (SH) coefficients. The proposed liver-shape classification algorithm consists of the following steps: (a) Preprocessing, including mesh generation and simplification, point-set matching, and surface to template alignment; (b) Liver-shape parameterization, including surface normalization, SH expansion followed by parameter space registration; (c) Feature selection and classification, including frequency based feature selection, feature space reduction by Principal Component Analysis (PCA), and classification. The above multi-step approach is novel in the sense that registration and feature selection for liver-shape classification is proposed and implemented and validated for the normal and diseases liver in the SH domain. Various groups of SH features after applying conventional PCA and/or ordered by p-value PCA are employed in two classifiers including Support Vector Machine (SVM) and k-Nearest Neighbor (k-NN) in the presence of 101 liver data sets. Results show that the proposed specific features combined with classifiers outperform existing liver-shape classification techniques that employ liver surface information in the spatial domain. In the available data sets, the proposed method can successful classify normal and diseased livers with a correct classification rate of above 90 %. The performed result in average is higher than conventional liver-shape classification method. Several standard metrics such as Leave-one-out cross-validation and Receiver Operating Characteristic (ROC) analysis are employed in the experiments and confirm the effectiveness of the proposed liver-shape classification with respect to conventional techniques.

Comparison of organs' shapes with geometric and Zernike 3D moments

The morphological similarity of organs is studied with feature vectors based on geometric and Zernike 3D moments. It is particularly investigated if outliers and average models can be identified. For this purpose, the relative proximity to the mean feature vector is defined, principal coordinate and clustering analyses are also performed. To study the consistency and usefulness of this approach, 17 livers and 76 hearts voxel models from several sources are considered. In the liver case, models with similar morphological feature are identified. For the limited amount of studied cases, the liver of the ICRP male voxel model is identified as a better surrogate than the female one. For hearts, the clustering analysis shows that three heart shapes represent about 80% of the morphological variations. The relative proximity and clustering analysis rather consistently identify outliers and average models. For the two cases, identification of outliers and surrogate of average models is rather robust. However, deeper classification of morphological feature is subject to caution and can only be performed after cross analysis of at least two kinds of feature vectors. Finally, the Zernike moments contain all the information needed to re-construct the studied objects and thus appear as a promising tool to derive statistical organ shapes.