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

Arxiv:1906.04333V1 [Eess.IV] 11 Jun 2019 Est As an Effective Way for Quantitative Tumor Characteriza- Window Approaches 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
Recommended publications
  • A New Family of Odd Generalized Nakagami (Nak-G) Distributions
    TURKISH JOURNAL OF SCIENCE http:/dergipark.gov.tr/tjos VOLUME 5, ISSUE 2, 85-101 ISSN: 2587–0971 A New Family of Odd Generalized Nakagami (Nak-G) Distributions Ibrahim Abdullahia, Obalowu Jobb aYobe State University, Department of Mathematics and Statistics bUniversity of Ilorin, Department of Statistics Abstract. In this article, we proposed a new family of generalized Nak-G distributions and study some of its statistical properties, such as moments, moment generating function, quantile function, and prob- ability Weighted Moments. The Renyi entropy, expression of distribution order statistic and parameters of the model are estimated by means of maximum likelihood technique. We prove, by providing three applications to real-life data, that Nakagami Exponential (Nak-E) distribution could give a better fit when compared to its competitors. 1. Introduction There has been recent developments focus on generalized classes of continuous distributions by adding at least one shape parameters to the baseline distribution, studying the properties of these distributions and using these distributions to model data in many applied areas which include engineering, biological studies, environmental sciences and economics. Numerous methods for generating new families of distributions have been proposed [8] many researchers. The beta-generalized family of distribution was developed , Kumaraswamy generated family of distributions [5], Beta-Nakagami distribution [19], Weibull generalized family of distributions [4], Additive weibull generated distributions [12], Kummer beta generalized family of distributions [17], the Exponentiated-G family [6], the Gamma-G (type I) [21], the Gamma-G family (type II) [18], the McDonald-G [1], the Log-Gamma-G [3], A new beta generated Kumaraswamy Marshall-Olkin- G family of distributions with applications [11], Beta Marshall-Olkin-G family [2] and Logistic-G family [20].
    [Show full text]
  • Dictionary-Based Stochastic Expectation– Maximization for SAR
    188 IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING, VOL. 44, NO. 1, JANUARY 2006 Dictionary-Based Stochastic Expectation– Maximization for SAR Amplitude Probability Density Function Estimation Gabriele Moser, Member, IEEE, Josiane Zerubia, Fellow, IEEE, and Sebastiano B. Serpico, Senior Member, IEEE Abstract—In remotely sensed data analysis, a crucial problem problem as a parameter estimation problem. Several strategies is represented by the need to develop accurate models for the have been proposed in the literature to deal with parameter statistics of the pixel intensities. This paper deals with the problem estimation, e.g., the maximum-likelihood methodology [18] of probability density function (pdf) estimation in the context of and the “method of moments” [41], [46], [53]. On the contrary, synthetic aperture radar (SAR) amplitude data analysis. Several theoretical and heuristic models for the pdfs of SAR data have nonparametric pdf estimation approaches do not assume any been proposed in the literature, which have been proved to be specific analytical model for the unknown pdf, thus providing effective for different land-cover typologies, thus making the a higher flexibility, although usually involving internal archi- choice of a single optimal parametric pdf a hard task, especially tecture parameters to be set by the user [18]. In particular, when dealing with heterogeneous SAR data. In this paper, an several nonparametric kernel-based estimation and regression innovative estimation algorithm is described, which faces such a problem by adopting a finite mixture model for the amplitude architectures have been described in the literature, that have pdf, with mixture components belonging to a given dictionary of proved to be effective estimation tools, such as standard Parzen SAR-specific pdfs.
    [Show full text]
  • Estimation and Testing Procedures of the Reliability Functions of Nakagami Distribution
    Austrian Journal of Statistics January 2019, Volume 48, 15{34. AJShttp://www.ajs.or.at/ doi:10.17713/ajs.v48i3.827 Estimation and Testing Procedures of the Reliability Functions of Nakagami Distribution Ajit Chaturvedi Bhagwati Devi Rahul Gupta Dept. of Statistics, Dept. of Statistics, Dept. of Statistics, University of Delhi University of Jammu University of Jammu Abstract A very important distribution called Nakagami distribution is taken into considera- tion. Reliability measures R(t) = P r(X > t) and P = P r(X > Y ) are considered. Point as well as interval procedures are obtained for estimation of parameters. Uniformly Mini- mum Variance Unbiased Estimators (U.M.V.U.Es) and Maximum Likelihood Estimators (M.L.Es) are developed for the said parameters. A new technique of obtaining these esti- mators is introduced. Moment estimators for the parameters of the distribution have been found. Asymptotic confidence intervals of the parameter based on M.L.E and log(M.L.E) are also constructed. Then, testing procedures for various hypotheses are developed. At the end, Monte Carlo simulation is performed for comparing the results obtained. A real data analysis is performed to describe the procedure clearly. Keywords: Nakagami distribution, point estimation, testing procedures, confidence interval, Markov chain Monte Carlo (MCMC) procedure.. 1. Introduction and preliminaries R(t), the reliability function is the probability that a system performs its intended function without any failure at time t under the prescribed conditions. So, if we suppose that life- time of an item or any system is denoted by the random variate X then the reliability is R(t) = P r(X > t).
    [Show full text]
  • Newdistns: an R Package for New Families of Distributions
    JSS Journal of Statistical Software March 2016, Volume 69, Issue 10. doi: 10.18637/jss.v069.i10 Newdistns: An R Package for New Families of Distributions Saralees Nadarajah Ricardo Rocha University of Manchester Universidade Federal de São Carlos Abstract The contributed R package Newdistns written by the authors is introduced. This pack- age computes the probability density function, cumulative distribution function, quantile function, random numbers and some measures of inference for nineteen families of distri- butions. Each family is flexible enough to encompass a large number of structures. The use of the package is illustrated using a real data set. Also robustness of random number generation is checked by simulation. Keywords: cumulative distribution function, probability density function, quantile function, random numbers. 1. Introduction Let G be any valid cumulative distribution function defined on the real line. The last decade or so has seen many approaches proposed for generating new distributions based on G. All of these approaches can be put in the form F (x) = B (G(x)) , (1) where B : [0, 1] → [0, 1] and F is a valid cumulative distribution function. So, for every G one can use (1) to generate a new distribution. The first approach of the kind of (1) proposed in recent years was that due to Marshall and Olkin(1997). In Marshall and Olkin(1997), B was taken to be B(p) = βp/ {1 − (1 − β)p} for β > 0. The distributions generated in this way using (1) will be referred to as Marshall Olkin G distributions. Since Marshall and Olkin(1997), many other approaches have been proposed.
    [Show full text]
  • University of Nevada, Reno Generalized Univariate
    University of Nevada, Reno Generalized Univariate Distributions and a New Asymmetric Laplace Model A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science in Mathematics By Palash Sharma Dr. Tomasz Kozubowski/Thesis Advisor August, 2017 c 2017 Palash Sharma ALL RIGHTS RESERVED THE GRADUATE SCHOOL We recommend that the thesis prepared under our supervision by Palash Sharma entitled Generalized Univariate Distributions and a New Asymmetric Laplace Model be accepted in partial fulfillment of the requirements for the degree of MASTER OF SCIENCE Tomasz J. Kozubowski, Ph.D., Advisor Anna Panorska, Ph.D., Committee Member Minggen Lu, Ph.D., Graduate School Representative David Zeh, Ph.D., Dean, Graduate School August, 2017 i ABSTRACT Generalized Univariate Distributions and a New Asymmetric Laplace Model By Palash Sharma This work provides a survey of general class of distributions generated from a mixture of beta random variables. We provide an extensive review of the literature, concern- ing generating new distributions via the inverse CDF transformation. In particular, we account for beta generated and Kumaraswamy generated families of distributions. We provide a brief summary of each of their families of distributions. We also propose a new asymmetric mixture distribution, which is an alternative to beta generated dis- tributions. We provide basic properties of this new class of distributions generated from the Laplace model. We also address the issue of parameter estimation of this new skew generalized Laplace model. ii ACKNOWLEDGMENTS At first, I would like to thank my honorable thesis advisor, Professor Tomasz J. Kozubowski, who showed me a great interest in the field of statistics and probability theory.
    [Show full text]
  • Field Guide to Continuous Probability Distributions
    Field Guide to Continuous Probability Distributions Gavin E. Crooks v 1.0.0 2019 G. E. Crooks – Field Guide to Probability Distributions v 1.0.0 Copyright © 2010-2019 Gavin E. Crooks ISBN: 978-1-7339381-0-5 http://threeplusone.com/fieldguide Berkeley Institute for Theoretical Sciences (BITS) typeset on 2019-04-10 with XeTeX version 0.99999 fonts: Trump Mediaeval (text), Euler (math) 271828182845904 2 G. E. Crooks – Field Guide to Probability Distributions Preface: The search for GUD A common problem is that of describing the probability distribution of a single, continuous variable. A few distributions, such as the normal and exponential, were discovered in the 1800’s or earlier. But about a century ago the great statistician, Karl Pearson, realized that the known probabil- ity distributions were not sufficient to handle all of the phenomena then under investigation, and set out to create new distributions with useful properties. During the 20th century this process continued with abandon and a vast menagerie of distinct mathematical forms were discovered and invented, investigated, analyzed, rediscovered and renamed, all for the purpose of de- scribing the probability of some interesting variable. There are hundreds of named distributions and synonyms in current usage. The apparent diver- sity is unending and disorienting. Fortunately, the situation is less confused than it might at first appear. Most common, continuous, univariate, unimodal distributions can be orga- nized into a small number of distinct families, which are all special cases of a single Grand Unified Distribution. This compendium details these hun- dred or so simple distributions, their properties and their interrelations.
    [Show full text]
  • Probability Models and Statistical Tests for Extreme Precipitation Based on Generalized Negative Binomial Distributions
    mathematics Article Probability Models and Statistical Tests for Extreme Precipitation Based on Generalized Negative Binomial Distributions Victor Korolev 1,2,3,4 and Andrey Gorshenin 1,2,3,* 1 Moscow Center for Fundamental and Applied Mathematics, Lomonosov Moscow State University, 119991 Moscow, Russia; [email protected] 2 Faculty of Computational Mathematics and Cybernetics, Lomonosov Moscow State University, 119991 Moscow, Russia 3 Federal Research Center “Computer Science and Control” of the Russian Academy of Sciences, 119333 Moscow, Russia 4 Department of Mathematics, School of Science, Hangzhou Dianzi University, Hangzhou 310018, China * Correspondence: [email protected] Received: 4 April 2020; Accepted: 14 April 2020; Published: 16 April 2020 Abstract: Mathematical models are proposed for statistical regularities of maximum daily precipitation within a wet period and total precipitation volume per wet period. The proposed models are based on the generalized negative binomial (GNB) distribution of the duration of a wet period. The GNB distribution is a mixed Poisson distribution, the mixing distribution being generalized gamma (GG). The GNB distribution demonstrates excellent fit with real data of durations of wet periods measured in days. By means of limit theorems for statistics constructed from samples with random sizes having the GNB distribution, asymptotic approximations are proposed for the distributions of maximum daily precipitation volume within a wet period and total precipitation volume for a wet period. It is shown that the exponent power parameter in the mixing GG distribution matches slow global climate trends. The bounds for the accuracy of the proposed approximations are presented. Several tests for daily precipitation, total precipitation volume and precipitation intensities to be abnormally extremal are proposed and compared to the traditional PoT-method.
    [Show full text]
  • Concentration Inequalities from Likelihood Ratio Method
    Concentration Inequalities from Likelihood Ratio Method ∗ Xinjia Chen September 2014 Abstract We explore the applications of our previously established likelihood-ratio method for deriving con- centration inequalities for a wide variety of univariate and multivariate distributions. New concentration inequalities for various distributions are developed without the idea of minimizing moment generating functions. Contents 1 Introduction 3 2 Likelihood Ratio Method 4 2.1 GeneralPrinciple.................................... ...... 4 2.2 Construction of Parameterized Distributions . ............. 5 2.2.1 WeightFunction .................................... .. 5 2.2.2 ParameterRestriction .............................. ..... 6 3 Concentration Inequalities for Univariate Distributions 7 3.1 BetaDistribution.................................... ...... 7 3.2 Beta Negative Binomial Distribution . ........ 7 3.3 Beta-Prime Distribution . ....... 8 arXiv:1409.6276v1 [math.ST] 1 Sep 2014 3.4 BorelDistribution ................................... ...... 8 3.5 ConsulDistribution .................................. ...... 8 3.6 GeetaDistribution ................................... ...... 9 3.7 GumbelDistribution.................................. ...... 9 3.8 InverseGammaDistribution. ........ 9 3.9 Inverse Gaussian Distribution . ......... 10 3.10 Lagrangian Logarithmic Distribution . .......... 10 3.11 Lagrangian Negative Binomial Distribution . .......... 10 3.12 Laplace Distribution . ....... 11 ∗The author is afflicted with the Department of Electrical
    [Show full text]
  • MPS: an R Package for Modelling New Families of Distributions
    MPS: An R package for modelling new families of distributions Mahdi Teimouri Department of Statistics Gonbad Kavous University Gonbad Kavous, IRAN Abstract: We introduce an R package, called MPS, for computing the probability density func- tion, computing the cumulative distribution function, computing the quantile function, simulating random variables, and estimating the parameters of 24 new shifted families of distributions. By considering an extra shift (location) parameter for each family more flexibility yields. Under some situations, since the maximum likelihood estimators may fail to exist, we adopt the well-known maximum product spacings approach to estimate the parameters of shifted 24 new families of distributions. The performance of the MPS package for computing the cdf, pdf, and simulating random samples will be checked by examples. The performance of the maximum product spacings approach is demonstrated by executing MPS package for three sets of real data. As it will be shown, for the first set, the maximum likelihood estimators break down but MPS package find them. For the second set, adding the location parameter leads to acceptance the model while absence of the location parameter makes the model quite inappropriate. For the third set, presence of the location parameter yields a better fit. Keywords: Cumulative distribution function; Maximum likelihood estimation; Method of maxi- mum product spacings; Probability density function; Quantile function; R package; Simulation; 1 Introduction Over the last two decades, generalization of the statistical distributions has attracted much at- tention in the literature. Most of these extensions have been spawned by applications found in arXiv:1809.02959v1 [stat.CO] 9 Sep 2018 analyzing lifetime data.
    [Show full text]
  • The Estimation of the M Parameter of the Nakagami Distribution
    WSEAS TRANSACTIONS on BIOLOGY and BIOMEDICINE Li-Fei Huang, Jen-Jen Lin The estimation of the m parameter of the Nakagami distribution Li-Fei Huang Jen-Jen Lin Ming Chuan University Ming Chuan University Dept. of Applied Stat. and Info. Science Dept. of Applied Stat. and Info. Science 5 Teh-Ming Rd., Gwei-Shan Taoyuan City 5 Teh-Ming Rd., Gwei-Shan Taoyuan City Taiwan Taiwan [email protected] [email protected] Abstract: This paper introduces the Nakagami-m distribution which is usually used to simulate the ultrasound image. The gamma distribution is used to derive the moment estimator and the maximum likelihood estimator because the Nakagami distribution has no moment generation function and too complicate likelihood function. The moment estimator of m is normally distributed with a smaller bias and a larger standard deviation. The second order maximum likelihood estimator of m is also normally distributed with a larger bias and a smaller standard deviation. The confidence interval for the ratio of medians from two independent distributions of Nakagami-m estimators is constructed. The moment estimator provides a quick understanding about the m parameter, while the second order maximum likelihood estimator provides a full understanding about the m parameter. Key–Words: The Nakagami distribution, The moment estimator, The maximum likelihood estimator, The distribu- tion of the Nakagami-m, The ratio of medians. 1 Introduction agami distribution is as follows. 2 m m 2m−1 − m n2 f(n) = ( ) n e Ω The Nakagami distribution is usually used to simu- Γ(m) Ω late the ultrasound image.
    [Show full text]
  • Assessing Probabilistic Modelling for Wind Speed from Numerical Weather Prediction Model and Observation in the Arctic
    www.nature.com/scientificreports OPEN Assessing probabilistic modelling for wind speed from numerical weather prediction model and observation in the Arctic Hao Chen1*, Yngve Birkelund1, Stian Normann Anfnsen2, Reidar Staupe‑Delgado1 & Fuqing Yuan1 Mapping Arctic renewable energy resources, particularly wind, is important to ensure the transition into renewable energy in this environmentally vulnerable region. The statistical characterisation of wind is critical for efectively assessing energy potential and planning wind park sites and is, therefore, an important input for wind power policymaking. In this article, diferent probability density functions are used to model wind speed for fve wind parks in the Norwegian Arctic region. A comparison between wind speed data from numerical weather prediction models and measurements is made, and a probability analysis for the wind speed interval corresponding to the rated power, which is largely absent in the existing literature, is presented. The results of the present study suggest that no single probability function outperforms across all scenarios. However, some diferences emerged from the models when applied to diferent wind parks. The Nakagami and Generalised extreme value distributions were chosen for the numerical weather predicted prediction and the observed wind speed modelling, respectively, due to their superiority and stability compared with other methods. This paper, therefore, provides a novel direction for understanding the numerical weather prediction wind model and shows that its speed statistical features are better captured than those of real wind. With the growing reliance on renewable energy resources in many regions of the world, studying the predictabil- ity of renewable energy is becoming progressively important1. As one of the cleanest renewable energy sources, wind energy has attracted growing attention worldwide 2.
    [Show full text]
  • The Generalized Odd Nakagami-G Family of Distributions: Properties and Applications
    NATURENGS, MTU Journal of Engineering and Natural Sciences 1:2 (2020) 1-16 NATURENGS MTU Journal of Engineering and Natural Sciences https://dergipark.org.tr/tr/pub/naturengs DOI: 10.46572/nat.2020.12 The Generalized Odd Nakagami-G Family of Distributions: Properties and Applications Ibrahim ABDULLAHI1*, Job OBALOWU2 1Department of Mathematics and Statistics, Faculty of Science, Yobe State University, 620, Yobe, Nigeria. 2Department of Statistics, Faculty of Physical Science, University of Ilorin, 240003, Kwara, Nigeria. (Received: 15.08.2020; Accepted: 01.10.2020) ABSTRACT: In this study, the Generalized Odd Nakagami-G distribution has been empirically investigated using Frechet distribution as the baseline distribution. Some mathematical statistics properties viz: quantile function, moments, probability-weighted moment, entropies, and order statistics were derived among others. The method of maximum likelihood estimation is suitable for the derivation of estimators of the distribution parameters. And the Generalized odd Nakagami Frechet (GONak-Fr) distribution gives the best fit vis-a-vis its competitors via application to real life data sets. Keywords: GONak-G, GONak-Fr, Probability-weighted moment, Entropies, Order statistics. 1. INTRODUCTION There have been recent developments focus on generalized odd classes of continuous distributions by adding at least one shape parameters to the baseline distribution, has led mathematical statistician researcher to make of the new flexible distributions and studying the properties of these distributions
    [Show full text]