Auranuch Lorsakul

Quantitative validation of optical flow based myocardial strain measures using sonomicrometry

Dynamic cardiac metrics, including myocardial strains and displacements, provide a quantitative approach to evaluate cardiac function. However, in current clinical diagnosis, largely 2D strain measures are used despite that cardiac motions are complex 3D volumes over time. Recent advances in 4D ultrasound enable the capability to capture such complex motion in a single image data set. In our previous work, a 4D optical flow based motion tracking algorithm was developed to extract full 4D dynamic cardiac metrics from such 4D ultrasound data. In order to quantitatively evaluate this tracking method, in-vivo coronary artery occlusion experiments at various locations were performed on three canine hearts. Each dog was screened with 4D ultrasound and sonomicrometry data was acquired during each occlusion study. The 4D ultrasound data from these experiments was then analyzed with the tracking method and estimated principal strain measures were directly compared to those recorded by sonomicrometry. Strong agreement was observed independently for the three canine hearts. This is the first validation study of optical flow based strain estimation for 4D ultrasound with a direct comparison with sonomicrometry using in-vivo data.

Coronary Occlusion Detection with 4D Optical Flow Based Strain Estimation on 4D Ultrasound

Real-time three-dimensional echocardiography (RT3DE) offers an efficient way to obtain complete 3D images of the heart over an entire cardiac cycle in just a few seconds. The complex 3D wall motion and temporal information contained in these 4D data sequences has the potential to enhance and supplement other imaging modalities for clinical diagnoses based on cardiac motion analysis. In our previous work, a 4D optical flow based method was developed to estimate dynamic cardiac metrics, including strains anddisplacements, from 4D ultrasound. In this study, in order to evaluate the ability of our method in detecting ischemic regions, coronary artery occlusion experiments at various locations were performed on five dogs. 4D ultrasound data acquired during these experiments were analyzed with our proposed method. Ischemic regions predicted by the outcome of strain measurements were compared to predictions from cardiac physiology with strong agreement. This is the first direct validation study of our image analysis method for biomechanical prediction and in vivo experimental outcome.

Point-cloud-to-point-cloud technique on tool calibration for dental implant surgical path tracking

Dental implant is one of the most popular methods of tooth root replacement used in prosthetic dentistry. Computerize navigation system on a pre-surgical plan is offered to minimize potential risk of damage to critical anatomic structures of patients. Dental tool tip calibrating is basically an important procedure of intraoperative surgery to determine the relation between the hand-piece tool tip and hand-pieces markers. With the transferring coordinates from preoperative CT data to reality, this parameter is a part of components in typical registration problem. It is a part of navigation system which will be developed for further integration. A high accuracy is required, and this relation is arranged by point-cloud-to-point-cloud rigid transformations and singular value decomposition (SVD) for minimizing rigid registration errors. In earlier studies, commercial surgical navigation systems from, such as, BrainLAB and Materialize, have flexibility problem on tool tip calibration. Their systems either require a special tool tip calibration device or are unable to change the different tool. The proposed procedure is to use the pointing device or hand-piece to touch on the pivot and the transformation matrix. This matrix is calculated every time when it moves to the new position while the tool tip stays at the same point. The experiment acquired on the information of tracking device, image acquisition and image processing algorithms. The key success is that point-to-point-cloud requires only 3 post images of tool to be able to converge to the minimum errors 0.77%, and the obtained result is correct in using the tool holder to track the path simulation line displayed in graphic animation.

Toward robot-assisted dental surgery: Path generation and navigation system using optical tracking approach

The main aim of this paper is to develop dental implant surgical navigation system based on homogenous transformation algorithms. This work is a partial section of robot-assisted surgical development. The previous works are presented in numerous basic research. They are methodology design on tool tip calibration, optical marker recognition, and pose determination using neural networks. This paper concerns with tracking path generation system based on fundamental of optical tracking. The intraoperative system is the principal focus area of this study. The homogenous transformation has been calculated in term of kinematics equation among marker relationship. The stereo camera is utilized to retrieve 3D position of different pattern styles of markers. The beneath marker recognition algorithm using rotation-invariant neural network and physical method is performed to identify markers. The fundamental relationship among markers are computed to obtain the orientation and translation between the guided path and the instruments tool tip. The experiment has been demonstrated and performed under prototype model. The method is to work on procedure step by step. They begin with patient information input and continuously perform on marker recognition, tool tip calibrations and marker digitization. The path tracking is executed to observe the accuracy of the system. The result shows that the system can be performed to track path based on beforehand planning.

Lumen Border Detection of Intravascular Ultrasound via Denoising of Directional Wavelet Representations

In this paper, intravascular ultrasound (IVUS) grayscale images, acquired with a single-element mechanically rotating transducer, are processed with wavelet denoising and region-based segmentation to extract various layers of lumen contours and plaques. First, IVUS volumetric data is expanded on complex exponential multi-resolution basis functions, also known as Brushlets, which are well localized in the time and frequency domains. Brushlet denoising has previously demonstrated a great aptitude for denoising ultrasound data and removal of blood speckle. A region-based segmentation framework is then applied for detection of lumen border layers, which remains a challenging problem in IVUS image analysis for images acquired with a single element, mechanically rotating 45 MHz transducer. We evaluated a hard thresholding operator for Brushlet denoising, and compared segmentation results to manually traced lumen borders. We observed good agreement and suggest that the proposed algorithm has a potential to be used as a reliable pre-processing step for accurate lumen border detection.

Spectral CT Using Multiple Balanced K-Edge Filters

Our goal is to validate a spectral computed tomography (CT) system design that uses a conventional X-ray source with multiple balanced K-edge filters. By performing a simultaneously synthetic reconstruction in multiple energy bins, we obtained a good agreement between measurements and model expectations for a reasonably complex phantom. We performed simulation and data acquisition on a phantom containing multiple rods of different materials using a NeuroLogica CT scanner. Five balanced K-edge filters including Molybdenum, Cerium, Dysprosium, Erbium, and Tungsten were used separately proximal to the X-ray tube. For each sinogram bin, measured filtered vector can be defined as a product of a transmission matrix, which is determined by the filters and is independent of the imaging object, and energy-binned intensity vector. The energy-binned sinograms were then obtained by inverting the transmission matrix followed by a multiplication of the filter measurement vector. For each energy bin defined by two consecutive K-edges, a synthesized energy-binned attenuation image was obtained using filtered back-projection reconstruction. The reconstructed attenuation coefficients for each rod obtained from the experiment was in good agreement with the corresponding simulated results. Furthermore, the reconstructed attenuation coefficients for a given energy bin, agreed with National Institute of Standards and Technology reference values when beam hardening within the energy bin is small. The proposed cost-effective system design using multiple balanced K-edge filters can be used to perform spectral CT imaging at clinically relevant flux rates using conventional detectors and integrating electronics.

In-vivo clinical validation of cardiac deformation and strain measurements from 4D ultrasound

An important goal in clinical cardiology is the non-invasive quantification of regional cardiac deformation. While many methods have been proposed for the estimation of 3D left ventricular deformation and strains derived from 4D ultrasound, currently there is a lack of in vivo clinical validation of these algorithms on humans. In this paper, we describe the experiments used in validating cardiac deformation and strain estimates of 4D ultrasound using correlation-based optical flow tracking on two different COPD patients with normal left ventricular function. Validation of the algorithm was done by 1) validation of cardiac volume across multiple scans of the same patient and 2) validation of the repeatability of cardiac displacement and strain results from multiple scan acquisitions of the same patient. The preliminary results are encouraging with our algorithm producing consistent cardiac volume and strain results across multiple acquisitions. Furthermore, our derived 4D cardiac strains showed qualitatively correct results. We also observed particularly interesting results in the radial displacements of the posterior and lateral walls of our COPD patients.

Impact of temporal resolution on LV myocardial regional strain assessment with real-time 3D ultrasound

Non-invasive quantification of regional left ventricular (LV) deformation is crucial for the identification of clinical and subclinical myocardial dysfunction in various conditions. Several software tools now exist to provide regional LV strain estimation for echocardiography images. In this paper, we experimentally investigated the impact of real-time three-dimensional (RT3D) ultrasound temporal resolution on the precision of an integrated speckle-tracking framework. We compared temporal displacement and strain profiles acquired at three different frame rates on five normal volunteers. Results showed that estimated displacement fields and regional strain measurements were more homogeneous and of larger amplitude at higher frame rates.

Numerical observer for atherosclerotic plaque classification in spectral computed tomography

Abstract. Spectral computed tomography (SCT) generates better image quality than conventional computed tomography (CT). It has overcome several limitations for imaging atherosclerotic plaque. However, the literature evaluating the performance of SCT based on objective image assessment is very limited for the task of discriminating plaques. We developed a numerical-observer method and used it to assess performance on discrimination vulnerable-plaque features and compared the performance among multienergy CT (MECT), dual-energy CT (DECT), and conventional CT methods. Our numerical observer was designed to incorporate all spectral information and comprised two-processing stages. First, each energy-window domain was preprocessed by a set of localized channelized Hotelling observers (CHO). In this step, the spectral image in each energy bin was decorrelated using localized prewhitening and matched filtering with a set of Laguerre–Gaussian channel functions. Second, the series of the intermediate scores computed from all the CHOs were integrated by a Hotelling observer with an additional prewhitening and matched filter. The overall signal-to-noise ratio (SNR) and the area under the receiver operating characteristic curve (AUC) were obtained, yielding an overall discrimination performance metric. The performance of our new observer was evaluated for the particular binary classification task of differentiating between alternative plaque characterizations in carotid arteries. A clinically realistic model of signal variability was also included in our simulation of the discrimination tasks. The inclusion of signal variation is a key to applying the proposed observer method to spectral CT data. Hence, the task-based approaches based on the signal-known-exactly/background-known-exactly (SKE/BKE) framework and the clinical-relevant signal-known-statistically/background-known-exactly (SKS/BKE) framework were applied for analytical computation of figures of merit (FOM). Simulated data of a carotid-atherosclerosis patient were used to validate our methods. We used an extended cardiac-torso anthropomorphic digital phantom and three simulated plaque types (i.e., calcified plaque, fatty-mixed plaque, and iodine-mixed blood). The images were reconstructed using a standard filtered backprojection (FBP) algorithm for all the acquisition methods and were applied to perform two different discrimination tasks of: (1) calcified plaque versus fatty-mixed plaque and (2) calcified plaque versus iodine-mixed blood. MECT outperformed DECT and conventional CT systems for all cases of the SKE/BKE and SKS/BKE tasks (all p<0.01). On average of signal variability, MECT yielded the SNR improvements over other acquisition methods in the range of 46.8% to 65.3% (all p<0.01) for FBP-Ramp images and 53.2% to 67.7% (all p<0.01) for FBP-Hanning images for both identification tasks. This proposed numerical observer combined with our signal variability framework is promising for assessing material characterization obtained through the additional energy-dependent attenuation information of SCT. These methods can be further extended to other clinical tasks such as kidney or urinary stone identification applications.

4D numerical observer for lesion detection in respiratory-gated PET

PURPOSE Respiratory-gated positron emission tomography (PET)/computed tomography protocols reduce lesion smearing and improve lesion detection through a synchronized acquisition of emission data. However, an objective assessment of image quality of the improvement gained from respiratory-gated PET is mainly limited to a three-dimensional (3D) approach. This work proposes a 4D numerical observer that incorporates both spatial and temporal informations for detection tasks in pulmonary oncology. METHODS The authors propose a 4D numerical observer constructed with a 3D channelized Hotelling observer for the spatial domain followed by a Hotelling observer for the temporal domain. Realistic (18)F-fluorodeoxyglucose activity distributions were simulated using a 4D extended cardiac torso anthropomorphic phantom including 12 spherical lesions at different anatomical locations (lower, upper, anterior, and posterior) within the lungs. Simulated data based on Monte Carlo simulation were obtained using geant4 application for tomographic emission (GATE). Fifty noise realizations of six respiratory-gated PET frames were simulated by GATE using a model of the Siemens Biograph mMR scanner geometry. PET sinograms of the thorax background and pulmonary lesions that were simulated separately were merged to generate different conditions of the lesions to the background (e.g., lesion contrast and motion). A conventional ordered subset expectation maximization (OSEM) reconstruction (5 iterations and 6 subsets) was used to obtain: (1) gated, (2) nongated, and (3) motion-corrected image volumes (a total of 3200 subimage volumes: 2400 gated, 400 nongated, and 400 motion-corrected). Lesion-detection signal-to-noise ratios (SNRs) were measured in different lesion-to-background contrast levels (3.5, 8.0, 9.0, and 20.0), lesion diameters (10.0, 13.0, and 16.0 mm), and respiratory motion displacements (17.6-31.3 mm). The proposed 4D numerical observer applied on multiple-gated images was compared to the conventional 3D approach applied on the nongated and motion-corrected images. RESULTS On average, the proposed 4D numerical observer improved the detection SNR by 48.6% (p < 0.005), whereas the 3D methods on motion-corrected images improved by 31.0% (p < 0.005) as compared to the nongated method. For all different conditions of the lesions, the relative SNR measurement (Gain = SNRObserved/SNRNongated) of the 4D method was significantly higher than one from the motion-corrected 3D method by 13.8% (p < 0.02), where Gain4D was 1.49 ± 0.21 and Gain3D was 1.31 ± 0.15. For the lesion with the highest amplitude of motion, the 4D numerical observer yielded the highest observer-performance improvement (176%). For the lesion undergoing the smallest motion amplitude, the 4D method provided superior lesion detectability compared with the 3D method, which provided a detection SNR close to the nongated method. The investigation on a structure of the 4D numerical observer showed that a Laguerre-Gaussian channel matrix with a volumetric 3D function yielded higher lesion-detection performance than one with a 2D-stack-channelized function, whereas a different kind of channels that have the ability to mimic the human visual system, i.e., difference-of-Gaussian, showed similar performance in detecting uniform and spherical lesions. The investigation of the detection performance when increasing noise levels yielded decreasing detection SNR by 27.6% and 41.5% for the nongated and gated methods, respectively. The investigation of lesion contrast and diameter showed that the proposed 4D observer preserved the linearity property of an optimal-linear observer while the motion was present. Furthermore, the investigation of the iteration and subset numbers of the OSEM algorithm demonstrated that these parameters had impact on the lesion detectability and the selection of the optimal parameters could provide the maximum lesion-detection performance. The proposed 4D numerical observer outperformed the other observers for the lesion-detection task in various lesion conditions and motions. CONCLUSIONS The 4D numerical observer shows substantial improvement in lesion detectability over the 3D observer method. The proposed 4D approach could potentially provide a more reliable objective assessment of the impact of respiratory-gated PET improvement for lesion-detection tasks. On the other hand, the 4D approach may be used as an upper bound to investigate the performance of the motion correction method. In future work, the authors will validate the proposed 4D approach on clinical data for detection tasks in pulmonary oncology.