Effect of adaptive threshold filtering on ultrasonic nakagami parameter to detect variation in scatterer concentration

The Nakagami parameter is associated with the Nakagami distribution estimated from ultrasonic backscattered signals and closely reflects the scatterer concentrations in tissues. There is an interest in exploring the possibility of enhancing the ability of the Nakagami parameter to characterize tissues. In this paper, we explore the effect of adaptive thresholdfiltering based on the noise-assisted empirical mode decomposition of the ultrasonic backscattered signals on the Nakagami parameter as a function of scatterer concentration for improving the Nakagami parameter performance. We carried out phantom experiments using 5 MHz focused and nonfocused transducers. Before filtering, the dynamic ranges of the Nakagami parameter, estimated using focused and nonfocused transducers between the scatterer concentrations of 2 and 32 scatterers/mm3, were 0.44 and 0.1, respectively. After filtering, the dynamic ranges of the Nakagami parameter, using the focused and nonfocused transducers, were 0.71 and 0.79, respectively. The experimental results showed that the adaptive threshold filter makes the Nakagami parameter measured by a focused transducer more sensitive to the variation in the scatterer concentration. The proposed method also endows the Nakagami parameter measured by a nonfocused transducer with the ability to differentiate various scatterer concentrations. However, the Nakagami parameters estimated by focused and nonfocused transducers after adaptive threshold filtering have different physical meanings: the former represents the statistics of signals backscattered from unresolvable scatterers while the latter is associated with stronger resolvable scatterers or local inhomogeneity due to scatterer aggregation.     

The effect of transducer characteristics on the estimation of Nakagami paramater as a function of scatterer concentration

The effect of transducer characteristics on the sensitivity of the Nakagami parameter to detect the variation of scatterer concentrations was studied. The rationale for this study stems from our pilot results which showed that the Nakagami parameters, estimated using a nonfocused transducer were not as sensitive as those of measurements using a commercial ultrasonic scanner in previous reports. This discrepancy may be attributed to the effects of transducer characteristics relative to the size of the resolution cell as verified by measurements of phantoms and 2-D computer simulations. The Nakagami parameter as a function of scatterer concentration was calculated using backscattered signals acquired from the scattering medium of different scatterer concentrations ranging from 2 to 32 scatterers/mm(3) using both 5 MHz nonfocused and focused transducers. Experimental and simulation results obtained from the nonfocused transducer represent that their respective Nakagami parameters increased from 1.17 to 1.31 and from 0.82 to 1.01 corresponding to the increase of scatterer concentrations. For the results obtained from the focused transducer, their average Nakagami parameters increased from 0.27 to 0.72 and from 0.33 to 0.81. These consistent results demonstrated that Nakagami parameter estimated using a focused transducer tends to be more sensitive than that by a nonfocused transducer to detect the variation of low scatterer concentration. This difference is fully due to the effect of transducer characteristics associated with the effective number of scatterers in the resolution cell.     

Imaging local scatterer concentrations by the Nakagami statistical model

The ultrasonic B-mode image is an important clinical tool used to examine the internal structures of the biological tissue. Due to the fact that the conventional B-scans cannot fully reflect the nature of the tissue, some useful quantitative parameters have been applied to quantify the properties of the tissue. Among various possibilities, the Nakagami parameter was demonstrated to have an outstanding ability to detect the variation of the scatterer concentration. This study is aimed to develop a scatterer concentration image based on the Nakagami parameter map to assist in the B-mode image for tissue characterization. In particular, computer simulations are carried out to generate phantoms of different scatterer concentrations and echogenicity coefficients and their B-mode and Nakagami parametric images are compared to evaluate the performance of the Nakagami image in differentiating the properties of the scatterers. The simulated results show that the B-mode image would be affected by the system settings and user operations, whereas the Nakagami parametric image provides a comparatively consistent image result when different diagnosticians use different dynamic ranges and system gains. This is largely because the Nakagami image formation is only based on the backscattered statistics of the ultrasonic signals in local tissues. Such an imaging principle allows the Nakagami image to quantify the local scatterer concentrations in the tissue and to extract the backscattering information from the regions of the weaker echoes that may be lost in the B-mode image. These findings suggest that the Nakagami image can be combined with the use of the B-mode image simultaneously to visualize the tissue structures and the scatterer properties for a better medical diagnosis.     

Statistical variations of ultrasound signals backscattered from flowing blood

The statistical distributions of ultrasonic signals backscattered from blood have recently been used to characterize hemodynamic properties, such as red blood cell (RBC) aggregation and blood coagulation. However, a thorough understanding of the relationship between blood properties and the statistical behavior of signals backscattered from flowing blood is still lacking. This prompted us to use the statistical parameter to characterize signals backscattered from both whole blood and RBC suspensions at different flow velocities (from 10 to 60 cm/s) and hematocrits (from 20% to 50%) under a steady laminar flow condition. The Nakagami parameter, scaling parameter, backscatter amplitude profile and flow velocity profile across a flow tube were acquired using a 10 MHz focused ultrasonic transducer. The backscattered signal peaked approximately at the centerline of the flow tube due to the effects of RBC aggregation, with the peak value increasing as the flow velocity of whole blood decreased. The Nakagami parameter increased from 0.45 to 0.78 as the flow velocity increased from 10 to 60 cm/s. The probability density function (PDF) of signals backscattered from flowing whole blood conformed with a pre-Rayleigh distribution. The Nakagami parameter was close to 1 for signals backscattered from RBC suspensions at all the flow velocities and hematocrits tested, for which the PDF was Rayleigh distributed. These differences in the statistical distributions of backscattered signals between whole blood and RBC suspensions suggest that variations in the size of dynamic scatterers in the flow affect the shape of the backscattered signal envelope, which should be considered in future statistical models used to characterize blood properties. (E-mail: j648816n@ms23.hinet.net and shyhhau@cycu.edu.tw)     

Ultrasonic Nakagami Imaging: A Strategy to Visualize the Scatterer Properties of Benign and Malignant Breast Tumors

Previous studies have demonstrated the usefulness of the Nakagami parameter in characterizing breast tumors by ultrasound. However, physicians or radiologists may need imaging tools in a clinical setting to visually identify the properties of breast tumors. This study proposed the ultrasonic Nakagami image to visualize the scatterer properties of breast tumors and then explored its clinical performance in classifying benign and malignant tumors. Raw data of ultrasonic backscattered signals were collected from 100 patients (50 benign and 50 malignant cases) using a commercial ultrasound scanner with a 7.5 MHz linear array transducer. The backscattered signals were used to form the B-scan and the Nakagami images of breast tumors. For each tumor, the average Nakagami parameter was calculated from the pixel values in the region-of-interest in the Nakagami image. The receiver operating characteristic (ROC) curve was used to evaluate the clinical performance of the Nakagami image. The results showed that the Nakagami image shadings in benign tumors were different from those in malignant cases. The average Nakagami parameters for benign and malignant tumors were 0.69 ± 0.12 and 0.55 ± 0.12, respectively. This means that the backscattered signals received from malignant tumors tend to be more pre-Rayleigh distributed than those from benign tumors, corresponding to a more complex scatterer arrangement or composition. The ROC analysis showed that the area under the ROC curve was 0.81 ± 0.04 and the diagnostic accuracy was 82%, sensitivity was 92% and specificity was 72%. The results showed that the Nakagami image is useful to distinguishing between benign and malignant breast tumors.     

On estimating the number density of random scatterers from backscattered acoustic signals

When an acoustic pulse interacts with an inhomogeneous, attenuating medium, the backscattered signals exhibit random fluctuations which are correlated with the physical properties of the medium. This paper proposes a robust model for characterizing the statistical nature of these backscattered signals. This model takes into account frequency-dependent attenuation, spatially varying media statistics, arbitrary beam geometries, and arbitrary pulse shapes. Based on this model, statistical estimation schemes are proposed for estimating both the attenuation coefficient and scatterer number density of the medium. Using appropriate simplifying assumptions, it is shown that this model is consistent with attenuation estimation algorithms currently used for ultrasonic tissue characterization. A statistical approach for estimating the number density of scatterers is described and its theoretical performance is evaluated. The algorithm for estimating the scatterer number density incorporates measurements of both the statistical moments of the backscattered signals and the point spread function of the acoustic system. The number density algorithm has been applied to simulated waveforms, waveforms obtained from ultrasonic phantoms with known number densities, and in vitro mammalian tissues. There is an excellent agreement between theoretical, simulation, and experimental results. The application of this technique to ultrasonic tissue characterization is also discussed.     

Changes in Backscattered Ultrasonic Envelope Statistics as a Function of Thrombus Age: An in Vitro Study

It is necessary to determine the age of thrombi in planning clinical treatment for thrombolysis. Ultrasound imaging can potentially be used to evaluate thrombus age in real time. The backscattered signals from thrombi may contain useful information regarding their age. On the basis of the randomness of ultrasound backscattering, this study explored changes in backscattered US statistics as a function of thrombus age. Porcine blood samples were used for the in vitro induction of fresh thrombi (day 0) with hematocrits ranging from 0%–40% and aged thrombi (days 0–8) with a hematocrit of 40%. Each thrombus was imaged using a pulse-echo ultrasound scanner equipped with a 7.5-MHz linear array transducer to acquire raw backscattered signals for B-mode and Nakagami imaging, by which the backscattered statistics were visualized. Hematoxylin and eosin staining and scanning electron microscopy were used to observe the histology of fresh and aged thrombi. The results indicated that a decrease in the number of red blood cells in the thrombus caused by the aging effect was observed in the in vitro model, indicating that the proposed model could simulate the structural changes in the thrombus during aging. Compared with fresh thrombi with various hematocrits, the aged thrombi exhibited a trend toward more substantial decreases in the Nakagami parameter with increasing thrombus age (the Nakagami parameter decreased from 1.1 to 0.6 as thrombus age increased from day 0 to day 8), indicating that thrombus aging causes the backscattered statistics to follow a pre-Rayleigh distribution to a high degree. This finding may be applied to the determination of thrombus age using conventional ultrasound imaging in the future.     

Microvascular Flow Estimation by Microbubble-Assisted Nakagami Imaging

The destruction and replenishment of microbubbles has been previously applied to estimating blood flow in the microcirculation. The rate of increase of the time–intensity curve (TIC) due to microbubbles flowing into the region-of-interest (ROI) as measured from the conventional B-mode images reflects the flow velocity. In this study, we monitored microbubble replenishment using a new proposed approach called the time–Nakagami-parameter curve (TNC) obtained from the parametric image based on the Nakagami statistical parameter for quantifying the microvascular flow velocity. The Nakagami parameter is estimated from signal envelope to reflect the backscattered statistics. The feasibility of using the TNC to estimate the microvascular flow was explored by carrying out phantom measurements and in vivo animal experiments. The rates of increase of the TIC and TNC were quantified as the rate constants β_I and β_N of monoexponential fitted curves, respectively. The experimental results showed that β_N behaves similarly to the conventional β_I in quantifying the flow velocity. Moreover, the tolerance to the effects of clutter is greater for the TNC than for the TIC, which makes it possible to use β_N to differentiate various flow velocities even when the ROI contains nonperfused areas. This finding suggests that the TNC-based technique can be used as a complementary tool for the conventional TIC to improve measurement of blood flow in the microcirculation. (E-mail: ckyeh@mx.nthu.edu.tw)     

A simulation model for ultrasonic temperature imaging using change in backscattered energy

Ultrasound backscattered from tissue has previously been shown theoretically and experimentally to change predictably with temperature in the hyperthermia range, i.e., 37 degrees C to 45 degrees C, motivating use of the change in backscattered ultrasonic energy (CBE) for ultrasonic thermometry. Our earlier theoretical model predicts that CBE from an individual scatterer will be monotonic with temperature, with, e.g., positive change for lipid-based scatterers and negative for aqueous-based scatterers. Experimental results have previously confirmed the presence of these positive and negative changes in one-dimensional ultrasonic signals and in two-dimensional images acquired from in vitro bovine, porcine and turkey tissues. In order to investigate CBE for populations of scatterers, we have developed an ultrasonic image simulation model, including temperature dependence for individual scatterers based on predictions from our theoretical model. CBE computed from images simulated for populations of randomly distributed scatterers behaves similarly to experimental results, with monotonic variation for individual pixel measurements and for image regions. Effects on CBE of scatterer type and distribution, size of the image region and signal-to-noise ratio have been examined. This model also provides the basis for future work regarding significant issues relevant to temperature imaging based on ultrasonic CBE such as effects of motion on CBE, limitations of motion-compensation techniques and accuracy of temperature estimation, including tradeoffs between temperature accuracy and available spatial resolution.     

Performance evaluation of ultrasonic Nakagami image in tissue characterization

Conventional ultrasonic B-mode images qualitatively describe tissue structures but are unsuitable for quantitative analyses of scatterer properties. We have recently developed an ultrasonic parametric imaging technique based on the Nakagami statistical distribution that is able to quantify scatterer concentrations. The aim of the present study is to further explore both the behavior of a Nakagami image in characterizing different scatterer structures at different signal-to-noise ratios (SNRs) and the feasibility of Nakagami imaging using a general commercial ultrasound scanner for tissue examinations. Simulations, experiments on a tissue-mimicking phantom and in vitro measurements on a muscle tissue before and after microwave treatment were carried out. The SNR and contrast-to-noise ratio (CNR) were estimated to quantify image performance. The results demonstrate that a Nakagami image can differentiate different scatterer concentrations for single, hypoechoic and hyperechoic targets. Also, a Nakagami image, when combined with an ultrasound scanner, can complement the B-scan to characterize tissue and to identify the region of interest with a larger CNR. However, the noise effect can degrade the performance of a Nakagami image. When the signal SNR decreased to 15 dB in simulations and to 8 dB in experiments, the CNR of the hyperechoic Nakagami image decreased by 4% and 27%, respectively, and that of the hypoechoic one decreased by 42% and 80%, respectively. These results indicate that a Nakagami image behaves well in identifying regions with high scatterer concentrations but does not perform well when both the scatterer concentration and SNR are low.     

Three-dimensional high-frequency backscatter and envelope quantification of cancerous human lymph nodes

Quantitative imaging methods using high-frequency ultrasound (HFU) offer a means of characterizing biological tissue at the microscopic level. Previously, high-frequency, 3-D quantitative ultrasound (QUS) methods were developed to characterize 46 freshly-dissected lymph nodes of colorectal-cancer patients. 3-D ultrasound radiofrequency data were acquired using a 25.6-MHz center-frequency transducer and each node was inked before tissue fixation to recover orientation after sectioning for 3-D histological evaluation. Backscattered echo signals were processed using 3-D cylindrical regions-of-interest (ROIs) to yield four QUS estimates associated with tissue microstructure (i.e., effective scatterer size, acoustic concentration, intercept and slope). These QUS estimates, obtained by parameterizing the backscatter spectrum, showed great potential for cancer detection. In the present study, these QUS methods were applied to 112 lymph nodes from 77 colorectal and gastric cancer patients. Novel QUS methods parameterizing the envelope statistics of the ROIs using Nakagami and homodyned-K distributions were also developed; they yielded four additional QUS estimates. The ability of these eight QUS estimates to classify lymph nodes and detect cancer was evaluated using receiver operating characteristics (ROC) curves. An area under the ROC curve of 0.996 with specificity and sensitivity of 95% were obtained by combining effective scatterer size and one envelope parameter based on the homodyned-K distribution. Therefore, these advanced 3-D QUS methods potentially can be valuable for detecting small metastatic foci in dissected lymph nodes. (E-mail: jmamou@riversideresearch.org)     

Minimum requirement of artificial noise level for using noise-assisted correlation algorithm to suppress artifacts in ultrasonic Nakagami images

The Nakagami image is a complementary imaging mode for pulse-echo ultrasound B-scan to characterize tissues. White noise in anechoic areas induces artifacts in the Nakagami image. Recently, we proposed a noise-assisted correlation algorithm (NCA) for suppressing the Nakagami artifact. In the NCA, artificial white noise is intentionally added twice to backscattered signals to produce two noisy data, which are used to establish a correlation profile for rejecting noise. This study explored the effects of artificial noise level on the NCA to suppress the artifact of the Nakagami image. Simulations were conducted to produce B-mode images of anechoic regions under signal-to-noise ratios (SNRs) of 20, 10 and 5 dB. Various artificial noise levels ranging from 0.1- to 1-fold of the intrinsic noise amplitude were used in the NCA for constructing the Nakagami images. Phantom experiments were conducted to validate the performance of using the optimal artificial noise level suggested by the simulation results to suppress the Nakagami artifacts by the NCA. The simulation results indicated that the artifacts of the Nakagami image in the anechoic regions can be gradually suppressed by increasing the artificial noise level used in the NCA to improve the image contrast-to-noise ratio (CNR). The CNR of the Nakagami image reached 20 dB when the artificial noise level was 0.7-fold of the intrinsic noise amplitude. This criterion was demonstrated by the phantom results to provide the NCA with an excellent ability to obtain artifact-free Nakagami images.     

3-D H-Scan Ultrasound Imaging and Use of a Convolutional Neural Network for Scatterer Size Estimation

H-Scan ultrasound (US) is a new imaging technology that estimates the relative size of acoustic scattering objects and structures. The purpose of this study was to introduce a three-dimensional (3-D) H-scan US imaging approach for scatterer size estimation in volume space. Using a programmable research scanner (Vantage 256, Verasonics Inc, Kirkland, WA, USA) equipped with a custom volumetric imaging transducer (4 DL7, Vermon, Tours, France), raw radiofrequency (RF) data was collected for offline processing to generate H-scan US volumes. A deep convolutional neural network (CNN) was modified and used to achieve voxel mapping from the input H-scan US image to underlying scatterer size. Preliminary studies were conducted using homogeneous gelatin-based tissue-mimicking phantom materials embedded with acoustic scatterers of varying size (15 to 250 μm) and concentrations (0.1 to 1%). Two additional phantoms were embedded with 63 or 125 µm-sized microspheres and used to test CNN estimation accuracy. In vitro results indicate that 3-D H-scan US imaging can visualize the spatial distribution of acoustic scatterers of varying size at different concentrations (R² > 0.85, p < 0.03). The result of scatterer size estimation reveals that a CNN can achieve an average mapping accuracy of 93.3%. Overall, our preliminary in vitro findings reveal that 3-D H-scan US imaging allows the visualization of tissue scatterer patterns and incorporation of a CNN can be used to help estimate size of the acoustic scattering objects.     

Improving the statistics of quantitative ultrasound techniques with deformation compounding: an experimental study

Many quantitative ultrasound (QUS) techniques are based on estimates of the radio-frequency (RF) echo signal power spectrum. Historically, reliable spectral estimates required spatial averaging over large regions-of-interest (ROIs). Spatial compounding techniques have been used to obtain robust spectral estimates for data acquired over small regions of interest. A new technique referred to as “deformation compounding” is another method for providing robust spectral estimates over smaller regions of interest. Motion tracking software is used to follow an ROI while the tissue is deformed (typically by pressing with the transducer). The deformation spatially reorganizes the scatterers so that the resulting echo signal is decorrelated. The RF echo signal power spectrum for the ROI is then averaged over several frames of RF echo data as the tissue is deformed, thus, undergoing deformation compounding. More specifically, averaging spectral estimates among the uncorrelated RF data acquired following small deformations allows reduction in the variance of the power spectral density estimates and, thereby, improves accuracy of spectrum-based tissue property estimation. The viability of deformation compounding has been studied using phantoms with known attenuation and backscatter coefficients. Data from these phantoms demonstrates that a deformation of about 2% frame-to-frame average strain is sufficient to obtain statistically-independent echo signals (with correlations of less than 0.2). Averaging five such frames, where local scatterer reorganization has taken place due to mechanical deformations, reduces the average percent standard deviation among power spectra by 26% and averaging 10 frames reduces the average percent standard deviation by 49%. Deformation compounding is used in this study to improve measurements of backscatter coefficients. These tests show deformation compounding is a promising method to improve the accuracy of spectrum-based quantitative ultrasound for tissue characterization.     

Theoretical estimation of the temperature dependence of backscattered ultrasonic power for noninvasive thermometry

The backscattered signal received from an insonified volume of tissue depends on tissue properties, such as attenuation, velocity, density, and backscatter coefficient and on the characteristics of the transducer at the insonified volume. Analysis of scattering in response to a burst of insonification showed that the temperature dependence of backscattered power was dominated by the effect of temperature on the backscatter coefficient. The temperature dependence of attenuation had a small effect on backscattered power. Backscattered power was independent of effects of temperature on velocity. These results were seen in the analysis of two types of inhomogeneity: 1) an aqueous scatterer in a water-based medium and 2) a lipid-based scatterer in the same water-based medium. The temperature dependence of the backscatter coefficient was inferred assuming that the backscatter coefficient was proportional to the scattering crosssection of a small scatterer. Backscattered power increased nearly logarithmically with temperature over the range from 37° to 50°C. Our model predicted a change of 5 dB for the lipid scatterer and a change of up to 3 db for the aqueous-based scatterer over that temperature range. For situations in which temperature dependence of the backscattered power can be calibrated, it may be possible to use the backscattered power level to track temperature distributions in tissue.     

Characterization of Tissue Microstructure Scatterer Distribution with Spectral Correlation

Characterization of tissue microstructure from the backscattered ultrasound signal using the spectral autocorrelation (SAC) function provides information about the scatterer distribution in biological tissue. This paper demonstrates SAC capabilities in characterizing periodicities in A-scans due to regularity in the scatterer distribution. The A-scan is modelled as a cyclostationary signal, where the statistical parameters of the signal vary in time with single or multiple periodicities. This periodicity manifests itself as spectral peaks both in the power spectral density (PSD) and in the SAC. Periodicity in the PSD will produce a well defined dominant peak in the cepstrum, which has been used to determine the scatterer spacing. The relationship between the scatterer spacing and the spacing of the spectral peaks is established using a stochastic model of the echo-formation process from biological tissue. The distribution of the scatterers within the microstructure is modelled using a Gamma function, which offers a flexible method of simulating parametric regularity in the scatterer spacing. Simulations of the tissue microstructure for lower orders of regularity indicate that the SAC components reveal information about the scatterer spacing that are not seen in the PSD and the cepstrum. The echo-formation process is tested by simulating microstructure of varying regularity and analyzing their effect on the SAC, PSD and cepstrum. Experimental validation of the simulation results are provided using in vivo scans of the breast and liver tissue that show the presence of significant spectral correlation components in the SAC.     

Statistics of Boundaries in Ultrasonic B-Scan Images

The existence of edges and boundaries in regions of interest (ROIs) in B-scan images alters the statistics of the backscattered echo from the ROI. Boundaries are the result of at least two different types of scattering scenarios in tissue, and the Nakagami model, which is being used extensively in ultrasound, is unlikely to fit the statistics of the backscattered echo under these conditions. Furthermore, there are very few other statistical models exist that describe the statistics of the backscattered echo from regions containing boundaries. In this work, the gamma mixture density and the recently proposed McKay density are explored as two viable models to fill this void. Justifications of these models are presented along with methods for estimating their parameters. Random number simulations and studies on tissue-mimicking phantoms indicate that the McKay and gamma mixture densities are the best for the modeling of the backscattered echo intensity when boundaries are present in the regions of interest.     

Model studies of angular behaviour of acoustic backscatter

This note describes some studies of the relationship between filtered 180° backscattered amplitude and rotation angle, which has been suggested as a means for tissue characterisation. The studies include experimental measurements on simple linear and two-dimensional point targets, and theoretical computed results on both regular arrays and arrays with random displacement of the point scatterers.

The studies were carried out on models under favourable conditions with only small amounts of irregularity. They show that in certain situations neither the disposition and spacing of peaks in the curve of backscattered amplitude versus angle, nor a cross-correlation algorithm designed to detect the optimum scatterer-spacing in the model gave a good estimate of average scatterer spacing. The implication of this finding to clinical application is discussed.     

A method for computer simulation of ultrasound doppler color flow images—I. theory and numerical method

Ultrasound imaging systems utilizing the pulsed Doppler principle are capable of providing images of blood flow in real time. We present a useful method for simulating flow images on a computer. Our method assumes that blood and surrounding tissue consist of many point-like scatterers positioned randomly in three dimensions. The position-dependent acoustic response of each scatterer is calculated using the acoustic impulse response method. This method takes into account the spatial effects of the transducer geometry on both the amplitude and temporal response of point-scattering. Details of theory, assumptions made in the simulation, and numerical methods are described fully for a spherically focused transducer, as well as a discussion of signal processing for generation of the flow image. Motion of a single scatterer is investigated to test the performance of the simulation algorithm. This simulation method could potentially be beneficial for detailed study of current and future flow imaging systems.     

Cataract measurement by estimating the ultrasonic statistical parameter using an ultrasound needle transducer: an in vitro study

A cataract is a clouding of the crystalline lens that reduces the amount of incoming light and impairs visual perception. Phacoemulsification is the most common surgical method for treating advanced cataracts, and the optimal phacoemulsification energy is determined by the lens hardness. A previous study proposed using the ultrasonic Nakagami image to complement the B-scan for distinguishing different degrees of lens hardening. However, it is difficult to implement the use of an imaging probe to detect the lens during phacoemulsification surgery in a clinical situation. To resolve this problem, this study applied an ultrasonic needle transducer to estimate the Nakagami parameter as an alternative for characterizing the cataract lens. Cataracts of porcine lenses were artificially induced in vitro, and the Young's modulus, backscattering intensities, and the Nakagami parameters were measured. The results showed that the backscattering intensity was not correlated with Young's modulus. In contrast, the average Nakagami parameter increased from 0.34 to 0.95 with increasing Young's modulus of the lens from 1.71 to 101 kPa. The above findings showed that the Nakagami parameter estimated with a needle transducer may be useful in differentiating different degrees of lens hardening, and implied that determining the optimal ultrasonic energy during clinical cataract surgery is possible if the needle transducer can be combined with the phacoemulsification probe to estimate the Nakagami parameter.