Arxiv:1906.04333V1 [Eess.IV] 11 Jun 2019 Est As an Effective Way for Quantitative Tumor Characteriza- Window Approaches
Total Page:16
File Type:pdf, Size:1020Kb
MULTISCALE NAKAGAMI PARAMETRIC IMAGING FOR IMPROVED LIVER TUMOR LOCALIZATION Omar S. Al-Kadi King Abdullah II School for IT University of Jordan Amman 11942, Jordan ABSTRACT shape of the Nakagami distribution is specified by the µ pa- rameter corresponding to the local concentration of scatterers, Effective ultrasound tissue characterization is usually hin- and the amount of spread (i.e. the local backscattered en- dered by complex tissue structures. The interlacing of speckle ergy) is represented by the scale parameter !. Different con- patterns complicates the correct estimation of backscatter ditions of the RF envelope statistics can be achieved by vary- distribution parameters. Nakagami parametric imaging based ing the µ parameter. Values of µ between 0 and 1 yield pre- on localized shape parameter mapping can model different Rayleigh and Rayleigh distributions. The Rayleigh distribu- backscattering conditions. However, performance of the con- tion case (µ = 1) resembles of having a large number of ran- structed Nakagami image depends on the sensitivity of the domly distributed scatterers, and in the case of high degree of estimation method to the backscattered statistics and scale of variance the distribution conforms to pre-Rayleigh (µ < 1). analysis. Using a fixed focal region of interest in estimating For a mixture of random and periodically located scatterers, the Nakagami parametric image would increase estimation the RF envelope statistics becomes a post-Rayleigh distribu- variance. In this work, localized Nakagami parameters are es- tion (µ > 1). The map of local µ parameter values – that timated adaptively by means of maximum likelihood estima- correspond to tissue properties – is normally considered in tion on a multiscale basis. The varying size kernel integrates constructing the Nakagami parametric image. The estimated the goodness-of-fit of the backscattering distribution param- Nakagami parameters as a function of the backscattered enve- eters at multiple scales for more stable parameter estimation. lope statistics has shown to be a reliable tool for quantitative Results show improved quantitative visualization of changes visualization of tissue structure changes [4, 5, 6]. in tissue specular reflections, suggesting a potential approach Previous work has improved local window-based µ pa- for improving tumor localization in low contrast ultrasound rameter estimation to generate the Nakagami parameter map images. from envelopes of raw ultrasound signals [7, 8, 9]. As using a Index Terms— Nakagami imaging, tumor detection, gamma kernel density estimation to achieve a smooth estima- maximum likelihood estimation, liver tumor, RF envelope tion of the distribution from small fixed-size windows [7], or by using a number of windows having a size 3 times the pulse length of the ultrasound [8], or by summing and averaging 1. INTRODUCTION multiple Nakagami parametric images generated using differ- ent sliding square window sizes (7-10 times the transducer Ultrasound parametric imaging is gaining increased inter- pulse length) [9]. However challenges persist with fixed-size arXiv:1906.04333v1 [eess.IV] 11 Jun 2019 est as an effective way for quantitative tumor characteriza- window approaches. A focal region of interest (i.e. using tion. Changes in properties of soft tissue texture, e.g. liver small windows) enables enhanced resolution of the Nakagami parenchyma, can be reflected in the radio-frequency (RF) parametric image, but large tissue structures require a large backscattered statistics as different Rayleigh distributions spatial scale to achieve stable parameter estimation. There- [1]. However parametric estimation is not a straightforward fore parameter smoothing might not suffice when prominent process and is generally faced with increased estimation vari- parts of the examined tissue structure is truncated or located ance, especially for complex speckle patterns. This may outside the window borders. On the other hand, using large obscure abnormal tissue structures, e.g. tumors and fibrosis, window sizes for summation and averaging may affect the re- which are deemed important for early diagnosis [2]. sults when compounded with windows of smaller-sizes, and The analytical simplicity of the bi-parametric Nakagami hence affecting the reliability of the constructed parametric distribution model, along with its goodness-of-fit with the en- image resolution. velope histogram of the ultrasound-backscattered signal [3], In this work an alternative approach of employing a multi- can be attractive for tissue characterization [1, 4, 5]. The scale kernel-based technique to model the backscattering dis- tribution statistics is proposed. The focal region of inter- est should be large enough to have sufficient tissue variation, while being also as small to avoid inclusion of irrelevant tex- tures from the nearby regions. The backscattered envelope from tissue was estimated voxel-by-voxel via Nakagami dis- tribution and subsequently used to generate optimized local parametric images for improved liver tumor detection. The assumption is based on that tumor regions tend to have differ- ent backscattered distribution than normal tissue [10], and a localized approach based on a varying-size kernel can assist in better identifying the tumor speckle patterns. 2. METHODOLOGY Fig. 1. Samples of simulated ultrasound speckle images rep- resenting: [left-right] fine texture (dense scatterers), coarse The speckle pattern is dependent on the ultrasound wave- texture (sparse scatterers), heterogeneous texture (random length and underlying tissue structure. As the former factor scatterers), and homogeneous texture (periodic scatterers), re- is fixed in this case and the latter varies across the ultrasound ferring to, respectively, phantoms D57, D30, D5, D37 in Ta- image, a multiscale approach that can localize the different ble 1. tissue structures would best suit the modeling of the backscat- tered envelope. A non-linear kernel is applied adaptively to characterize the tissue speckle patterns. The estimation of the been proposed as a general model for ultrasonic backscatter- best parametric Nakagami maps from varying size kernels is ing under different scattering conditions and scatterer densi- described as follows: ties [3]. This distribution has the density function 2.1. Multiscale kernel localization µ Let V be an order set of constructed envelope detected RF µ 1 (2µ−1) − µ x2 N(xjµ, !) = 2 x e ! ; 8x 2 R ≥ 0 images Ii (x; y), where i is a certain slice in the acquired vol- ! Γ(µ) ume, and Pµ,ω(x; y) are the corresponding µ and ! paramet- (2) ric images. The RF images are calculated from the envelope where x is the envelope of the RF signal and Γ(·) is the of the ultrasound-backscattered signal just before performing gamma function. If x has a Nakagami distribution N(x) with any intensity mapping and post processing filtering. This rep- parameters µ and !, then x2 has a gamma distribution Γ with resentation preserves the ultrasound data unaltered while pro- shape µ and scale (energy) parameter !/µ. Although there are viding better quantitative analysis, i.e. without the risk of los- other distributions exist in the literature for modeling ultra- ing information due to RF signal shaping. Then a set of vary- sonic backscattering, the Nakagami probabilistic distribution ing size kernels K can be defined for each Ii (x; y), where was chosen for its simplicity and ability to characterize differ- K = fv1; v2; : : : ; vkg ; vj 2 V and k = 1=8 of the size of ent scattering conditions ranging from pre- to post-Rayleigh Ii (x; y). Different focal regions are investigated in a multi- [6]. scale manner as in (1) by varying the size of two non-negative integer variables a and b, which are used to center each local- Each voxel l in Ii is adaptively transformed via K at different scales to its corresponding parametric Nakagami ized kernel vj(s; t) with a different size of m × n on each parameters by means of maximum likelihood estimation voxel l in Ii (x; y) of size M × N. (MLE) forming a set of parametric vectors for each voxel l. l l a b 2 The MLE θ^(v) for a density function f v ; : : : ; v jθ when X X k 1 k P (x; y) = v (s; t) I (x + s; y + t) θ is a vector of parameters for the Nakagami distribution µ,ω j i j s=−a t=−b family Θ, estimates the most probable parameters θ^(v) = (1) l l argmaxθ D θjv1; : : : ; v , where D (θjv) = f (vjθ) ,θ 2 Θ m+2 n+2 k where a = d 2 e, b = d 2 e, and m, n = 1; 2; : : : ; k. is the score function. Finally the goodness-of-fit is estimated via root mean square error for the calculated Nakagami pa- 2.2. Modeling backscattered statistics rameters θm at different scales, giving the localized paramet- ric Nakagami images Pµ,ω as summarized in Algorithm 1. The Nakagami distribution N(x) is a gamma related distribu- The shape parametric image Pµ is used for subsequent tissue tion which is known for its analytical simplicity [5], and has characterization. Algorithm 1 Multiscale kernel localization Table 1. Mean absolute difference comparison of estimated Input: Set of ultrasound backscattered envelope images I = i Nakagami shape parametric images against ground-truth. f(x1; y1 : : : ; xj; yj)g Output: Localized ultrasound Nakagami shape and scale Nakagami shape parametric images P ;P µ ! estimation methods for all voxels l in I do i Phantom GKF WMC MKL for all localized kernels νl ! νl do 1 k D5 0.29 0.10 0.04 fStep1g // Fit with a Nakagami distribution µ µ 1 (2µ−1) − µ x2 D11 0.18 0.09 0.04 N(xjµ, !) = 2 x e ! ! Γ(µ) D13 0.40