[en] In this paper, we put forward a systematic method to analyze, control and evaluate the safety issues of medical robotics. We created a safety model that consists of three axes to analyze safety factors. Software and hardware are the two material axes. The third axis is the policy that controls all phases of design, production, testing and application of the robot system. The policy was defined as hazard identification and safety insurance control (HISIC) that includes seven principles: definitions and requirements, hazard identification, safety insurance control, safety critical limits, monitoring and control, verification and validation, system log and documentation. HISIC was implemented in the development of a robot for urological applications that was known as URObot. The URObot is a universal robot with different modules adaptable for 3D ultrasound image-guided interstitial laser coagulation, radiation seed implantation, laser resection, and electrical resection of the prostate. Safety was always the key issue in the building of the robot. The HISIC strategies were adopted for safety enhancement in mechanical, electrical and software design. The initial test on URObot showed that HISIC had the potential ability to improve the safety of the system. Further safety experiments are being conducted in our laboratory

[en] A nonrigid B-spline-based point-matching (BPM) method is proposed to match dense surface points. The method solves both the point correspondence and nonrigid transformation without features extraction. The registration method integrates a motion model, which combines a global transformation and a B-spline-based local deformation, into a robust point-matching framework. The point correspondence and deformable transformation are estimated simultaneously by fuzzy correspondence and by a deterministic annealing technique. Prior information about global translation, rotation and scaling is incorporated into the optimization. A local B-spline motion model decreases the degrees of freedom for optimization and thus enables the registration of a larger number of feature points. The performance of the BPM method has been demonstrated and validated using synthesized 2D and 3D data, mouse MRI and micro-CT images. The proposed BPM method can be used to register feature point sets, 2D curves, 3D surfaces and various image data. (paper)

[en] High-frequency ultrasound (HFU) has the ability to image both skeletal and cardiac muscles. The quantitative assessment of these myofiber orientations has a number of applications in both research and clinical examinations; however, difficulties arise due to the severe speckle noise contained in the HFU images. Thus, for the purpose of automatically measuring myofiber orientations from two-dimensional HFU images, we propose a two-step multiscale image decomposition method. It combines a nonlinear anisotropic diffusion filter and a coherence enhancing diffusion filter to extract myofibers. This method has been verified by ultrasound data from simulated phantoms, excised fiber phantoms, specimens of porcine hearts, and human skeletal muscles in vivo. The quantitative evaluations of both phantoms indicated that the myofiber measurements of our proposed method were more accurate than other methods. The myofiber orientations extracted from different layers of the porcine hearts matched the prediction of an established cardiac mode and demonstrated the feasibility of extracting cardiac myofiber orientations from HFU images ex vivo. Moreover, HFU also demonstrated the ability to measure myofiber orientations in vivo. (paper)

[en] An automatic segmentation framework is proposed to segment the right ventricle (RV) in echocardiographic images. The method can automatically segment both epicardial and endocardial boundaries from a continuous echocardiography series by combining sparse matrix transform, a training model, and a localized region-based level set. First, the sparse matrix transform extracts main motion regions of the myocardium as eigen-images by analyzing the statistical information of the images. Second, an RV training model is registered to the eigen-images in order to locate the position of the RV. Third, the training model is adjusted and then serves as an optimized initialization for the segmentation of each image. Finally, based on the initializations, a localized, region-based level set algorithm is applied to segment both epicardial and endocardial boundaries in each echocardiograph. Three evaluation methods were used to validate the performance of the segmentation framework. The Dice coefficient measures the overall agreement between the manual and automatic segmentation. The absolute distance and the Hausdorff distance between the boundaries from manual and automatic segmentation were used to measure the accuracy of the segmentation. Ultrasound images of human subjects were used for validation. For the epicardial and endocardial boundaries, the Dice coefficients were 90.8 ± 1.7% and 87.3 ± 1.9%, the absolute distances were 2.0 ± 0.42 mm and 1.79 ± 0.45 mm, and the Hausdorff distances were 6.86 ± 1.71 mm and 7.02 ± 1.17 mm, respectively. The automatic segmentation method based on a sparse matrix transform and level set can provide a useful tool for quantitative cardiac imaging. (paper)

[en] We are investigating imaging techniques to study the tumor response to photodynamic therapy (PDT). Positron emission tomography (PET) can provide physiological and functional information. High-resolution magnetic resonance imaging (MRI) can provide anatomical and morphological changes. Image registration can combine MRI and PET images for improved tumor monitoring. In this study, we acquired high-resolution MRI and microPET ¹⁸F-fluorodeoxyglucose (FDG) images from C3H mice with RIF-1 tumors that were treated with Pc 4-based PDT. We developed two registration methods for this application. For registration of the whole mouse body, we used an automatic three-dimensional, normalized mutual information algorithm. For tumor registration, we developed a finite element model (FEM)-based deformable registration scheme. To assess the quality of whole body registration, we performed slice-by-slice review of both image volumes; manually segmented feature organs, such as the left and right kidneys and the bladder, in each slice; and computed the distance between corresponding centroids. Over 40 volume registration experiments were performed with MRI and microPET images. The distance between corresponding centroids of organs was 1.5±0.4 mm which is about 2 pixels of microPET images. The mean volume overlap ratios for tumors were 94.7% and 86.3% for the deformable and rigid registration methods, respectively. Registration of high-resolution MRI and microPET images combines anatomical and functional information of the tumors and provides a useful tool for evaluating photodynamic therapy

[en] Based on the Radon transform, a wavelet multiscale denoising method is proposed for MR images. The approach explicitly accounts for the Rician nature of MR data. Based on noise statistics we apply the Radon transform to the original MR images and use the Gaussian noise model to process the MR sinogram image. A translation invariant wavelet transform is employed to decompose the MR 'sinogram' into multiscales in order to effectively denoise the images. Based on the nature of Rician noise we estimate noise variance in different scales. For the final denoised sinogram we apply the inverse Radon transform in order to reconstruct the original MR images. Phantom, simulation brain MR images, and human brain MR images were used to validate our method. The experiment results show the superiority of the proposed scheme over the traditional methods. Our method can reduce Rician noise while preserving the key image details and features. The wavelet denoising method can have wide applications in MRI as well as other imaging modalities

[en] Common CT imaging signs of lung diseases (CISLs) are defined as the imaging signs that frequently appear in lung CT images from patients. CISLs play important roles in the diagnosis of lung diseases. This paper proposes a novel learning method, namely learning with distribution of optimized feature (DOF), to effectively recognize the characteristics of CISLs. We improve the classification performance by learning the optimized features under different distributions. Specifically, we adopt the minimum spanning tree algorithm to capture the relationship between features and discriminant ability of features for selecting the most important features. To overcome the problem of various distributions in one CISL, we propose a hierarchical learning method. First, we use an unsupervised learning method to cluster samples into groups based on their distribution. Second, in each group, we use a supervised learning method to train a model based on their categories of CISLs. Finally, we obtain multiple classification decisions from multiple trained models and use majority voting to achieve the final decision. The proposed approach has been implemented on a set of 511 samples captured from human lung CT images and achieves a classification accuracy of 91.96%. The proposed DOF method is effective and can provide a useful tool for computer-aided diagnosis of lung diseases on CT images. (paper)

[en] A three-dimensional (3D) mutual information registration method was created and used to register MRI volumes of the pelvis and prostate. It had special features to improve robustness. First, it used a multi-resolution approach and performed registration from low to high resolution. Second, it used two similarity measures, correlation coefficient at lower resolutions and mutual information at full resolution, because of their particular advantages. Third, we created a method to avoid local minima by restarting the registration with randomly perturbed parameters. The criterion for restarting was a correlation coefficient below an empirically determined threshold. Experiments determined the accuracy of registration under conditions found in potential applications in prostate cancer diagnosis, staging, treatment and interventional MRI (iMRI) guided therapies. Images were acquired in the diagnostic (supine) and treatment position (supine with legs raised). Images were also acquired as a function of bladder filling and the time interval between imaging sessions. Overall studies on three patients and three healthy volunteers, when both volumes in a pair were obtained in the diagnostic position under comparable conditions, bony landmarks and prostate 3D centroids were aligned within 1.6± 0.2mm and 1.4±0.2 mm, respectively, values only slightly larger than a voxel. Analysis suggests that actual errors are smaller because of the uncertainty in landmark localization and prostate segmentation. Between the diagnostic and treatment positions, bony landmarks continued to register well, but prostate centroids moved towards the posterior 2.8-3.4 mm. Manual cropping to remove voxels in the legs was necessary to register these images. In conclusion, automatic, rigid body registration is probably sufficiently accurate for many applications in prostate cancer. For potential iMRI-guided treatments, the small prostate displacement between the diagnostic and treatment positions can probably be avoided by acquiring volumes in similar positions and by reducing bladder and rectal volumes. (author)

[en] We are investigating three-dimensional (3D) to two-dimensional (2D) registration methods for computed tomography (CT) and dual-energy digital radiography (DEDR). CT is an established tool for the detection of cardiac calcification. DEDR could be a cost-effective alternative screening tool. In order to utilize CT as the ''gold standard'' to evaluate the capability of DEDR images for the detection and localization of calcium, we developed an automatic, intensity-based 3D-to-2D registration method for 3D CT volumes and 2D DEDR images. To generate digitally reconstructed radiography (DRR) from the CT volumes, we developed several projection algorithms using the fast shear-warp method. In particular, we created a Gaussian-weighted projection for this application. We used normalized mutual information (NMI) as the similarity measurement. Simulated projection images from CT values were fused with the corresponding DEDR images to evaluate the localization of cardiac calcification. The registration method was evaluated by digital phantoms, physical phantoms, and clinical data sets. The results from the digital phantoms show that the success rate is 100% with a translation difference of less than 0.8 mm and a rotation difference of less than 0.2 deg. . For physical phantom images, the registration accuracy is 0.43±0.24 mm. Color overlay and 3D visualization of clinical images show that the two images registered well. The NMI values between the DRR and DEDR images improved from 0.21±0.03 before registration to 0.25±0.03 after registration. Registration errors measured from anatomic markers decreased from 27.6±13.6 mm before registration to 2.5±0.5 mm after registration. Our results show that the automatic 3D-to-2D registration is accurate and robust. This technique can provide a useful tool for correlating DEDR with CT images for screening coronary artery calcification

[en] Purpose: To develop and test an automated algorithm to classify different types of tissue in dedicated breast CT images. Methods: Images of a single breast of five different patients were acquired with a dedicated breast CT clinical prototype. The breast CT images were processed by a multiscale bilateral filter to reduce noise while keeping edge information and were corrected to overcome cupping artifacts. As skin and glandular tissue have similar CT values on breast CT images, morphologic processing is used to identify the skin based on its position information. A support vector machine (SVM) is trained and the resulting model used to create a pixelwise classification map of fat and glandular tissue. By combining the results of the skin mask with the SVM results, the breast tissue is classified as skin, fat, and glandular tissue. This map is then used to identify markers for a minimum spanning forest that is grown to segment the image using spatial and intensity information. To evaluate the authors’ classification method, they use DICE overlap ratios to compare the results of the automated classification to those obtained by manual segmentation on five patient images. Results: Comparison between the automatic and the manual segmentation shows that the minimum spanning forest based classification method was able to successfully classify dedicated breast CT image with average DICE ratios of 96.9%, 89.8%, and 89.5% for fat, glandular, and skin tissue, respectively. Conclusions: A 2D minimum spanning forest based classification method was proposed and evaluated for classifying the fat, skin, and glandular tissue in dedicated breast CT images. The classification method can be used for dense breast tissue quantification, radiation dose assessment, and other applications in breast imaging