Improvement of ultrasound speckle image velocimetry using image enhancement techniques

Ultrasound-based techniques have been developed and widely used in noninvasive measurement of blood velocity. Speckle image velocimetry (SIV), which applies a cross-correlation algorithm to consecutive B-mode images of blood flow has often been employed owing to its better spatial resolution compared with conventional Doppler-based measurement techniques. The SIV technique utilizes speckles backscattered from red blood cell (RBC) aggregates as flow tracers. Hence, the intensity and size of such speckles are highly dependent on hemodynamic conditions. The grayscale intensity of speckle images varies along the radial direction of blood vessels because of the shear rate dependence of RBC aggregation. This inhomogeneous distribution of echo speckles decreases the signal-to-noise ratio (SNR) of a cross-correlation analysis and produces spurious results. In the present study, image-enhancement techniques such as contrast-limited adaptive histogram equalization (CLAHE), min/max technique, and subtraction of background image (SB) method were applied to speckle images to achieve a more accurate SIV measurement. A mechanical sector ultrasound scanner was used to obtain ultrasound speckle images from rat blood under steady and pulsatile flows. The effects of the image-enhancement techniques on SIV analysis were evaluated by comparing image intensities, velocities, and cross-correlation maps. The velocity profiles and wall shear rate (WSR) obtained from RBC suspension images were compared with the analytical solution for validation. In addition, the image-enhancement techniques were applied to in vivo measurement of blood flow in human vein. The experimental results of both in vitro and in vivo SIV measurements show that the intensity gradient in heterogeneous speckles has substantial influence on the cross-correlation analysis. The image-enhancement techniques used in this study can minimize errors encountered in ultrasound SIV measurement in which RBCs are used as flow tracers instead of exogenous contrast agents.     

Carr-Purcell-Meiboom-Gill imaging of prostate cancer: quantitative T2 values for cancer discrimination

Quantitative, apparent T(2) values of suspected prostate cancer and healthy peripheral zone tissue in men with prostate cancer were measured using a Carr-Purcell-Meiboom-Gill (CPMG) imaging sequence in order to assess the cancer discrimination potential of tissue T(2) values. The CPMG imaging sequence was used to image the prostates of 18 men with biopsy-proven prostate cancer. Whole gland coverage with nominal voxel volumes of 0.54 x 1.1 x 4 mm(3) was obtained in 10.7 min, resulting in data sets suitable for generating high-quality images with variable T(2)-weighting and for evaluating quantitative T(2) values on a pixel-by-pixel basis. Region-of-interest analysis of suspected healthy peripheral zone tissue and suspected cancer, identified on the basis of both T(1)- and T(2)-weighted signal intensities and available histopathology reports, yielded significantly (P<.0001) longer apparent T(2) values in suspected healthy tissue (193+/-49 ms) vs. suspected cancer (100+/-26 ms), suggesting potential utility of this method as a tissue specific discrimination index for prostate cancer. We conclude that CPMG imaging of the prostate can be performed in reasonable scan times and can provide advantages over T(2)-weighted fast spin echo (FSE) imaging alone, including quantitative T(2) values for cancer discrimination as well as proton density maps without the point spread function degradation associated with short effective echo time FSE sequences.     

Optical classification of random image fields using spectral synthetic discriminant functions

It is shown that the problem of classification of images that have the perfectly random nature may be solved with the help of synthetic discriminant functions being synthesized by least-squares technique to separate linearly the power spectra of the corresponding random image fields. The realization of the proposed method by means of an optical technique is discussed, and its efficiency is illustrated by two examples of real-life texture classification.     

Improvement of image quality using BLADE sequences in brain MR imaging

The purpose of this study is to compare two types of sequences in brain magnetic resonance (MR) examinations of uncooperative and cooperative patients. For each group of patients, the pairs of sequences that were compared were two T2-weighted (T2-W) fluid attenuated inversion recovery sequences with different k-space trajectories (conventional Cartesian and BLADE) and two T2-TSE weighted with different k-space trajectories (conventional Cartesian and BLADE). Twenty-three consecutive uncooperative patients and 44 cooperative patients, who routinely underwent brain MR imaging examination, participated in the study. Both qualitative and quantitative analyses were performed based on the signal-to-noise ratio, contrast-to-noise ratio (CNR), and relative contrast (ReCon) measures of normal anatomic structures. The qualitative analysis was performed by experienced radiologists. Also, the presence of motion, other (e.g., Gibbs, susceptibility artifacts, phase encoding from vessels) artifacts and pulsatile flow artifacts was evaluated.     

Integration of quantitative DCE-MRI and ADC mapping to monitor treatment response in human breast cancer: initial results

PURPOSE: The objective of this study was to assess changes in the water apparent diffusion coefficient (ADC) and in pharmacokinetic parameters obtained from the fast-exchange regime (FXR) modeling of dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) during neoadjuvant chemotherapy in breast cancer. MATERIALS AND METHODS: Eleven patients with locally advanced breast cancer underwent MRI examination prior to and after chemotherapy but prior to surgery. A 1.5-T scanner was used to obtain T1, ADC and DCE-MRI data. DCE-MRI data were analyzed by the FXR model returning estimates of K(trans) (volume transfer constant), v(e) (extravascular extracellular volume fraction) and tau(i) (average intracellular water lifetime). Histogram and correlation analyses assessed parameter changes post-treatment. RESULTS: Significant (P < .05) changes or trends towards significance (P < .10) were seen in all parameters except tau(i), although there was qualitative reduction in tau(i) values post-treatment. In particular, there was reduction (P < .035) in voxels with K(trans) values in the range 0.2-0.5 min(-1) and a decrease (P < .05) in voxels with ADC values in the range 0.99 x 10(-3) to 1.35 x 10(-3) mm2/s. ADC and v(e) were negatively correlated (r = -.60, P < .02). Parameters sensitive to water distribution and geometry (T(1), v(e), tau(i) and ADC) correlated with a multivariable linear regression model. CONCLUSION: The analysis presented here is sensitive to longitudinal changes in breast tumor status; K(trans) and ADC are most sensitive to these changes. Relationships between parameters provide information on water distribution and geometry in the tumor environment.     

Investigation of prostate cancer using diffusion-weighted intravoxel incoherent motion imaging

Purpose

The objective of this work was to evaluate the diagnostic performance of the intravoxel incoherent motion (IVIM) model to differentiate between healthy and malignant prostate tissue.

Materials and Methods

Regions of interest were drawn in healthy and cancerous tissue of 13 patients with histologically proven prostate carcinoma and fitted to a monoexponential model [yielding the apparent diffusion coefficient (ADC)] and the IVIM signal equation (yielding the perfusion fraction f, the diffusion constant D and the pseudodiffusion coefficient of perfusion D?). Parameter maps were calculated for all parameters.

Results

The ADC, D and f were significantly (P<.005) lowered in cancerous tissue (1.01±0.22 μm²/ms, 0.84±0.19 μm²/ms and 14.27±7.10%, respectively) compared to benign tissue (1.49±0.17 μm²/ms, 1.21±0.22 μm²/ms and 21.25±8.32%, respectively). Parameter maps of D and f allowed for a delineation of the tumor, but showed higher variations compared to the ADC map.

Conclusion

Apparent diffusion coefficient maps provide better diagnostic performance than IVIM maps for tumor detection. However, the results suggest that the reduction of the ADC in prostate cancer stems not only from changes in cellularity but also from perfusion effects. IVIM imaging might hold promise for the diagnosis of other prostatic lesions.     

Feasibility of diffusional kurtosis tensor imaging in prostate MRI for the assessment of prostate cancer: Preliminary results

Purpose

To assess the feasibility of full diffusional kurtosis tensor imaging (DKI) in prostate MRI in clinical routine. Histopathological correlation was achieved by targeted biopsy.

Materials and Methods

Thirty-one men were prospectively included in the study. Twenty-one were referred to our hospital with increased prostate specific antigen (PSA) values (> 4 ng/ml) and suspicion of prostate cancer. The other 10 men were volunteers without any history of prostate disease. DKI applying diffusion gradients in 20 different spatial directions with four b-values (0, 300, 600, 1000 s/mm²) was performed additionally to standard functional prostate MRI. Region of interest (ROI)-based measurements were performed in all histopathologically verified lesions of every patient, as well as in the peripheral zone, and the central gland of each volunteer.

Results

DKI showed a substantially better fit to the diffusion-weighted signal than the monoexponential apparent diffusion coefficient (ADC). Altogether, 29 lesions were biopsied in 14 different patients with the following results: Gleason score 3 + 3 = 6 (n = 1), 3 + 4 = 7 (n = 7), 4 + 3 = 7 (n = 6), 4 + 4 = 8 (n = 1), and 4 + 5 = 9 (n = 2), and prostatitis (n = 12). Values of axial (K_ax) and mean kurtosis (K_mean) were significantly different in the tumor (K_ax 1.78 ± 0.39, K_mean 1.84 ± 0.43) compared with the normal peripheral zone (K_ax 1.09 ± 0.12, K_mean 1.16 ± 0.13; p < 0.001) or the central gland (K_ax 1.40 ± 0.12, K_mean 1.44 ± 0.17; p = 0.01 respectively). There was a minor correlation between axial kurtosis (r = 0.19) and the Gleason score.

Conclusion

Full DKI is feasible to utilize in a routine clinical setting. Although there is some overlap some DKI parameters can significantly distinguish prostate cancer from the central gland or the normal peripheral zone. Nevertheless, the additional value of DKI compared with conventional monoexponential ADC calculation remains questionable and requires further research.     

DWI and PRECISE criteria in men on active surveillance for prostate cancer: A multicentre preliminary experience of different ADC calculations

PurposeThe PRECISE score estimates the likelihood of radiological progression in patients on active surveillance (AS) for prostate cancer (PCa) with serial multiparametric magnetic resonance imaging (mpMRI). A PRECISE score of 1 or 2 denotes radiological regression, PRECISE 3 indicates stability and PRECISE 4 or 5 implies progression.We evaluated the inter-reader reproducibility of different apparent diffusion coefficient (ADC) calculations and their relationship to the PRECISE score.Material and methodsBaseline and follow-up scans (on the same MR systems) of 30 patients with visible lesions from two different institutions (University College London and Sapienza University of Rome) were analysed by two radiologists (one from each site). The PRECISE score was initially assessed in consensus. At least six weeks later, to reduce the likelihood of being influenced by the consensus PRECISE reading, each radiologist independently calculated ADC for the following: lesion, non-cancerous tissue and urine in the bladder. Normalised ADC ratios were calculated with respect to normal prostatic tissue (npADC) and urine. Spearman's correlation (ρ), intraclass correlation coefficients (ICC), differences in ADC and ROC curves were computed.ResultsInterobserver reproducibility was very good (ρ > 0.8; ICC > 0.90). Lesion ADC (0.91 vs 0.73 × 10⁻³ mm²/s; p=0.025) and npADC ratio (0.68 vs 0.53; p=0.012) at follow-up mpMRI were different between patients with radiological regression or stability vs progression. Cut-offs of 0.77 × 10⁻³ mm²/s (lesion ADC) and 0.59 (npADC ratio) could differentiate the two groups (area under the curve: 0.74 and 0.77, respectively).ConclusionThe ADC, npADC ratio and the PRECISE score should be recorded for MRI-based AS.     

A comparative assessment of preclinical chemotherapeutic response of tumors using quantitative non-Gaussian diffusion MRI

BackgroundDiffusion-weighted MRI (DWI) signal attenuation is often not mono-exponential (i.e. non-Gaussian diffusion) with stronger diffusion weighting. Several non-Gaussian diffusion models have been developed and may provide new information or higher sensitivity compared with the conventional apparent diffusion coefficient (ADC) method. However the relative merits of these models to detect tumor therapeutic response is not fully clear.MethodsConventional ADC, and three widely-used non-Gaussian models, (bi-exponential, stretched exponential, and statistical model), were implemented and compared for assessing SW620 human colon cancer xenografts responding to barasertib, an agent known to induce apoptosis via polyploidy. Bayesian Information Criterion (BIC) was used for model selection among all three non-Gaussian models.ResultsAll of tumor volume, histology, conventional ADC, and three non-Gaussian DWI models could show significant differences between control and treatment groups after four days of treatment. However, only the non-Gaussian models detected significant changes after two days of treatment. For any treatment or control group, over 65.7% of tumor voxels indicate the bi-exponential model is strongly or very strongly preferred.ConclusionNon-Gaussian DWI model-derived biomarkers are capable of detecting tumor earlier chemotherapeutic response of tumors compared with conventional ADC and tumor volume. The bi-exponential model provides better fitting compared with statistical and stretched exponential models for the tumor and treatment models used in the current work.     

Motion corrected DWI with integrated T2-mapping for simultaneous estimation of ADC,T2-relaxation and perfusion in prostate cancer

ObjectiveMultiparametric magnetic resonance imaging (MRI) and PI-RADS (Prostate Imaging – Reporting and Data System) has become the standard to determine a probability score for a lesion being a clinically significant prostate cancer. T2-weighted and diffusion-weighted imaging (DWI) are essential in PI-RADS, depending partly on visual assessment of signal intensity, while dynamic-contrast enhanced imaging is less important. To decrease inter-rater variability and further standardize image evaluation, complementary objective measures are in need.MethodsWe here demonstrate a sequence enabling simultaneous quantification of apparent diffusion coefficient (ADC) and T2-relaxation, as well as calculation of the perfusion fraction f from low b-value intravoxel incoherent motion data. Expandable wait pulses were added to a FOCUS DW SE-EPI sequence, allowing the effective echo time to change at run time. To calculate both ADC and f, b-values 200 s/mm² and 600 s/mm² were chosen, and for T2-estimation 6 echo times between 64.9 ms and 114.9 ms were used.ResultsThree patients with prostate cancer were examined and all had significantly decreased ADC and T2-values, while f was significantly increased in 2 of 3 tumors. T2 maps obtained in phantom measurements and in a healthy volunteer were compared to T2 maps from a SE sequence with consecutive scans, showing good agreement. In addition, a motion correction procedure was implemented to reduce the effects of prostate motion, which improved T2-estimation.ConclusionsThis sequence could potentially enable more objective tumor grading, and decrease the inter-rater variability in the PI-RADS classification.     

A quantitative assessment of liver metabolites during jaundice using three dimensional phosphorus chemical shift imaging

Phosphorus metabolites in the jaundiced rat liver were studied by three-dimensional phosphorus chemical shift imaging (CSI). Animals were studied at 1, 2, and 3 weeks post-ligation of the common bile duct. Quantitation of metabolites was performed using an external standard. Metabolite T(1) values were assessed in CSI experiments on normal untreated animals. High-performance liquid chromatography (HPLC) was used to measure adenine nucleotides in a separate group of jaundiced rats. 3D-CSI did not detect significant changes in NTP in jaundiced animals relative to baseline controls. At two and three weeks post bile duct ligation, pH was significantly elevated. HPLC data comparing ATP levels to baseline controls also detected no change except for elevated ATP detected on Day 21. (31)P NMR chemical shift imaging may be used to assess liver metabolites under conditions of stress such as jaundice. However, absolute quantitation requires careful attention to many factors including point spread function, correct T(1) values, and adequate signal-to-noise ratio.     

3D quantitative image analysis of open-cell nickel foams under tension and compression loading using X-ray microtomography

The deformation behaviour and fracture of an open-cell nickel foam were analysed using X-ray microtomography at the ESRF, Grenoble, France. In situ tensile and compression tests were performed at a resolution of 2 and 10 μm. The initial morphology of the foam was studied using 3D image analysis. Parameters such as the cell volume and strut length distributions, number of faces per cell, number of nodes per face and the shape of the most representative cells were determined. The cells are shown to be non-spherical due to the initial geometrical anisotropy of the polyurethane foam template and to the load applied to the nickel foam during processing. This geometrical anisotropy is shown to be related to the observed anisotropy of the elastic properties of the material using a simple beam model. In tension, bending, stretching and alignment of struts are observed. A tensile test in the longitudinal direction is shown to reinforce the privileged orientations of the cells. In contrast, a tensile test in the transverse direction leads to a more isotropic distribution of the cells. These features are illustrated by pole figures of the three axes of equivalent ellipsoids for all cells at different strain levels. Compression tests are associated with strain localization phenomena due to the buckling of struts in a weaker region of the foam. Finally, study of open-cell nickel foam fracture shows that cracks initiate at nodes during tensile tests and that the damaged zone is about five cells wide. Free edge effects on crack initiation are also evidenced.     

High resolution quantitative SIMS analysis of shallow boron implants in silicon using a bevel and image approach

Secondary ion mass spectrometry (SIMS) is frequently used as the preferred tool for dopant profiling due to its sensitivity and depth resolution. However, as dopant profiles become shallower most, if not all of the implant profile lies in the pre-equilibrium or transient region of an SIMS depth profile. In this region sputter yield and ionisation rate vary making accurate quantification of the implant profile very difficult. These problems can be reduced through the use of much lower beam energies or oxygen flooding of the sample. However, most SIMS instruments do not have these capabilities. In this paper an alternative technique for producing an accurate depth profile of a shallow implant, using existing SIMS technology is presented.Through the fabrication of bevels with very small slope angles on a shallow boron implanted silicon via a chemical etch, SIMS ion imaging is performed on the exposed surface. Ion image data is then summed, and in conjunction with accurate measurement of the bevel morphology, a shallow boron implant profile produced. The ‘bevel-image’ profile compares very well with a profile obtained using a 1 keV oxygen beam. To ensure a good dynamic range on the ‘bevel-image’ profile it is important to clean the bevel with a HF etch, prior to imaging.     

CCD技术在塞曼效应实验中应用的改进

利用CCD图像采集和计算机图像程序处理技术改进了塞曼效应实验系统,对比了测微目镜测量干涉圆环直径、利用手动方法处理实验图像以及利用计算机程序处理实验图像3种方法的实验结果.实验结果表明:利用计算机程序自动处理实验图像的方法能够简化实验并减小实验误差.     

Enhancement of image quality by using metamaterial inspired energy harvester

Metamaterial energy harvester based subwavelength imaging structure is presented. Besides, 2×2 array patch antenna is designed to create incident radiated signal towards the proposed harvester structure. Harvested energy is converted to DC voltage signal by using Schottky diodes and each cell is represented by 256 grey levelled pixel value. Therefore, any incident electromagnetic wave placed in the operating frequency range can be monitorized. In the experimental study, experiments are realized with and without antenna to investigate the effects of instantaneous electromagnetic waves and the observed images. Afterwards, another experiment is conducted to observe the effects of metal plate which is located between the antenna and the harvester structure. The last experiment is performed by using a Yagi Uda antenna for sub-wavelength imaging. The proposed harvester structure can be used in various applications such as energy harvesting, incident wave tracing, crack detection, spy device detection, medical imaging, and so on.     

Evaluation of the quality of a speckle pattern in the digital image correlation method by mean subset fluctuation

The quality of the speckle pattern used in digital image correlation is studied in this paper using a parameter called mean subset fluctuation. A numerical translation test is performed on four speckle patterns captured from actual experiments. The translated images are analyzed and the results show that the mean bias error of the calculated displacement is linear with the value of the mean subset fluctuation. The results from the numerical calculation illustrate that speckle size and density have an influence on the quality of speckle pattern.     

Online biospeckle assessment without loss of definition and resolution by motion history image

The application of the dynamic laser speckle as a reliable instrument to achieve maps of activity in biological material is available in literature optics and laser. The application, particularly in live specimens, such as animals and human beings necessitated some approaches to avoid the kinking of the bodies, which creates changes in the patterns undermining the biological activity under monitoring. The adoption of online techniques circumvented the noise generated by the kinking, however, with considerable reduction in the resolution and definition of the activity maps. This work presents a feasible alternative to the routine online methods based on the Motion History Image (MHI) methodology. The adoption of MHI was tested in biological and non-biological samples and compared with online as well as offline procedures of biospeckle image analysis. Tests on paint drying was associated to alcohol volatilization, and tests on a maize seed and on growing of roots confirmed the hypothesis that the MHI would be able to implement an online approach without the reduction of resolution and definition on the resultant images, thereby presenting in some cases results that were comparable to the offline procedures.     

Assessment of robustness and transferability of classification models built for cancer diagnostics using Raman spectroscopy

Over recent years, Raman spectroscopy has been demonstrated as a prospective tool for application in cancer diagnostics. The use of Raman spectroscopy for this purpose relies on pattern recognition methods that have been developed to perform well on data achieved under laboratory conditions. However, the application of Raman spectroscopy as a routine clinical tool is likely to result in imperfect data due to instrument‐to‐instrument variation. Such corruption to the pure tissue spectral data is expected to negatively impact the classification performance of the diagnostic model. In this paper, we present a thorough assessment of the robustness of the Raman approach. This was achieved by perturbing a set of spectra in different ways, including various linear shifts, nonlinear shifts and random noise and using previously optimised classification models to predict the class membership of each spectrum in a testing set. The loss of predictive power with increased corruption was used to calculate a score, which allows an easy comparison of the model robustness. For this approach, three different types of classification models, including linear discriminant analysis (LDA), partial least square discriminant analysis (PLS‐DA) and support vector machine (SVM), built for lymph node diagnostics were the subject of the robustness testing. The results showed that a linear perturbation had the highest impact on the performance of all classification models. Among all linear corruption methods, a gradient y‐shift resulted in the highest performance loss. Thus, the factor most likely to affect the predictive outcome of models when using different systems is a gradient y‐shift. Copyright © 2010 John Wiley & Sons, Ltd.     

利用光斑图像特征确定飞秒激光有效烧蚀焦距

为确定飞秒激光光束对微尺度结构的烧蚀深度,研究了给定功率条件下对应的激光束有效烧蚀焦距。提出采用激光焦点处获得的烧痕阵列图像及在离焦状态下提取烧痕图像特征,通过分析图像特征与离焦距离,获得激光束有效烧蚀焦距范围的方法。在激光束焦点附近的硅晶片表面烧蚀出斑痕阵列,向下逐渐减小焦距,采集硅晶片斑痕图像,提取斑痕平均像素面积及斑痕目标与背景之间的R分量灰度差,获得斑痕像素面积及灰度差随激光束焦距变化的曲线;向上逐渐增大焦距,提取并获得斑痕像素面积及灰度差随激光束焦距变化的曲线。结合激光束向下离焦阈值(633 μm)及向上离焦阈值(993 μm),确定20 mW输出功率条件下,飞秒激光在硅晶片材料表面的有效烧蚀深度为360 μm。采用中位值方法确定了激光束在硅晶片表面聚焦时的焦距为0.823 mm。实验表明,激光烧蚀斑痕像素面积及灰度差与激光束焦距之间的关系能够客观地反映激光束有效烧蚀焦距的变化范围。     

Pre-operative assessment of residual disease in locally advanced breast cancer patients: A sequential study by quantitative diffusion weighted MRI as a function of therapy

PurposeThe potential of diffusion weighted imaging (DWI) in assessing pathologic response and surgical margins in locally advanced breast cancer patients (n = 38) undergoing neoadjuvant chemotherapy was investigated.MethodsDWI was performed at pre-therapy (Tp0), after I (Tp1) and III (Tp3) NACT at 1.5 T. Apparent diffusion coefficient (ADC) of whole tumor (ADC_WT), solid tumor (ADC_ST), intra-tumoral necrosis (ADC_Nec) was determined. Further, ADC of 6 consecutive shells (5 mm thickness each) including tumor margin to outside tumor margins (OM1 to OM5) was calculated and the data analyzed to define surgical margins.ResultsOf 38 patients, 6 were pathological complete responders (pCR), 19 partial responders (pPR) and 13 were non-responders (pNR). Significant increase was observed in ADC_ST and ADC_WT in pCR and pPR following therapy. Pre-therapy ADC was significantly lower in pCR compared to pPR and pNR indicating the heterogeneous nature of tumor which may affect drug perfusion and consequently the response. ADC of outside margins (OM1, OM2, and OM3) was significantly different among pCR, pPR and pNR at Tp3 which may serve as response predictive parameter. Further, at Tp3, ADC of outside margins (OM1, OM2, and OM3) was significantly lower compared to that seen at Tp0 in pCR, indicating the presence of residual disease in these shells.ConclusionPre-surgery information may serve as a guide to define cancer free margins and the extent of residual disease which may be useful in planning breast conservation surgery.