Exponential-spline-based Features for Ultrasound Signal Characterization
|
|
- Antony Richards
- 5 years ago
- Views:
Transcription
1 1 / 29 Exponential-spline-based Features for Ultrasound Signal Characterization Simona Maggio Department of Electronics, Computer Science and Systems (DEIS) University of Bologna Scuola di Dottorato: Scienze e Ingegneria dell Informazione Corso di Dottorato: Ingegneria Elettronica, Informatica e delle Telecomunicazioni 19/10/2009
2 2 / 29 Outline 1 Tissue Characterization for Ultrasound Diagnostic 2 Ultrasound Signals Modeling Predictive Deconvolution Continuous/Discrete Approach 3 Exponential Splines 4 From Continuous to Discrete E-spline-based Whitening Model Non Stationary Discrete Whitening 5 Conclusions Analysis Results Further Developments
3 2 / 29 Outline 1 Tissue Characterization for Ultrasound Diagnostic 2 Ultrasound Signals Modeling Predictive Deconvolution Continuous/Discrete Approach 3 Exponential Splines 4 From Continuous to Discrete E-spline-based Whitening Model Non Stationary Discrete Whitening 5 Conclusions Analysis Results Further Developments
4 2 / 29 Outline 1 Tissue Characterization for Ultrasound Diagnostic 2 Ultrasound Signals Modeling Predictive Deconvolution Continuous/Discrete Approach 3 Exponential Splines 4 From Continuous to Discrete E-spline-based Whitening Model Non Stationary Discrete Whitening 5 Conclusions Analysis Results Further Developments
5 2 / 29 Outline 1 Tissue Characterization for Ultrasound Diagnostic 2 Ultrasound Signals Modeling Predictive Deconvolution Continuous/Discrete Approach 3 Exponential Splines 4 From Continuous to Discrete E-spline-based Whitening Model Non Stationary Discrete Whitening 5 Conclusions Analysis Results Further Developments
6 2 / 29 Outline 1 Tissue Characterization for Ultrasound Diagnostic 2 Ultrasound Signals Modeling Predictive Deconvolution Continuous/Discrete Approach 3 Exponential Splines 4 From Continuous to Discrete E-spline-based Whitening Model Non Stationary Discrete Whitening 5 Conclusions Analysis Results Further Developments
7 3 / 29 Ultrasound Imaging Pros: real time, non invasive. Limits: low resolution. Tissue characterization: Highligthing features invisible by visual inspection Cancer detection and staging Avoiding unnecessary biopsy Specific application: Prostate cancer computer-aided detection
8 3 / 29 Ultrasound Imaging Pros: real time, non invasive. Limits: low resolution. Tissue characterization: Highligthing features invisible by visual inspection Cancer detection and staging Avoiding unnecessary biopsy Specific application: Prostate cancer computer-aided detection
9 4 / 29 Procedure for Tissue Characterization Figure: Feature extraction for ultrasound analysis
10 4 / 29 Procedure for Tissue Characterization Figure: Feature extraction for ultrasound analysis
11 4 / 29 Procedure for Tissue Characterization Figure: Feature extraction for ultrasound analysis
12 5 / 29 From Tissue to Echo Discrete model: z = HΣs + n = Hx + n Σ: coherent reflection, macroscopic interactions, mean value of diffused field. s: incoherent reflections, interactions smaller than wavelength, random fluctuations of diffused field. Scattering as generalized gaussian noise: p(s µ, b) = a e s µ b 2
13 5 / 29 From Tissue to Echo Discrete model: z = HΣs + n = Hx + n Σ: coherent reflection, macroscopic interactions, mean value of diffused field. s: incoherent reflections, interactions smaller than wavelength, random fluctuations of diffused field. Scattering as generalized gaussian noise: p(s µ, b) = a e s µ b 2
14 5 / 29 From Tissue to Echo Discrete model: z = HΣs + n = Hx + n Σ: coherent reflection, macroscopic interactions, mean value of diffused field. s: incoherent reflections, interactions smaller than wavelength, random fluctuations of diffused field. Scattering as generalized gaussian noise: p(s µ, b) = a e s µ b 2
15 5 / 29 From Tissue to Echo Discrete model: z = HΣs + n = Hx + n Σ: coherent reflection, macroscopic interactions, mean value of diffused field. s: incoherent reflections, interactions smaller than wavelength, random fluctuations of diffused field. Scattering as generalized gaussian noise: p(s µ, b) = a e s µ b 2
16 Deconvolution to recover Tissue Response Convolutional model: z[k] = x[k] h[k] + n[k] Deconvolution to restore tissue response: x[k] Point Spread Function (PSF) h[n] not known Blind adaptive deconvolution approach 1 Advantages: simplicity, low computational cost, variable PSF. 1 [Ng et al., 2007], [Michailovich and Adam, 2005], [Jensen, 1994],[Rasmussen, 1994] 6 / 29
17 7 / 29 Predictive Deconvolution (PD) W(z) predictive filter to remove predictable features of z[n] Ideal reconstruction if the transducer is an AR system and x[n] is white Gaussian whitening Recursive Least Squares solution for adaptive whitening z[n] as non stationary AR process and x[n] generalized Gaussian Recovering unpredictable part of RF signal: scattering
18 7 / 29 Predictive Deconvolution (PD) W(z) predictive filter to remove predictable features of z[n] Ideal reconstruction if the transducer is an AR system and x[n] is white Gaussian whitening Recursive Least Squares solution for adaptive whitening z[n] as non stationary AR process and x[n] generalized Gaussian Recovering unpredictable part of RF signal: scattering
19 Improvement due to Predictive Deconvolution Prostate Images No Preprocessing RLS Deconvolution SE 0.69 ± ± 0.09 SP 0.94 ± ± 0.01 Acc 0.93 ± ± 0.02 Az 0.92 ± ± / 29
20 Improvement due to Predictive Deconvolution Benignant case No Preprocessing RLS Deconvolution SE 0.69 ± ± 0.09 SP 0.94 ± ± 0.01 Acc 0.93 ± ± 0.02 Az 0.92 ± ± / 29
21 Improvement due to Predictive Deconvolution Malignant case No Preprocessing RLS Deconvolution SE 0.69 ± ± 0.09 SP 0.94 ± ± 0.01 Acc 0.93 ± ± 0.02 Az 0.92 ± ± / 29
22 Continuous/Discrete Approach for PD Continuous whitening g(t) continuous RF signal in a stationary case w(t) continuous uncorrelated unpredictable (scattering) Link with discrete: discretization in a spline basis Exponential splines: multiresolution version of L α Identification of continuous model from sampled data Shift-variant multiresolution description of scattering signal 9 / 29
23 Continuous/Discrete Approach for PD Continuous whitening g(t) continuous RF signal in a stationary case w(t) continuous uncorrelated unpredictable (scattering) Link with discrete: discretization in a spline basis Exponential splines: multiresolution version of L α Identification of continuous model from sampled data Shift-variant multiresolution description of scattering signal 9 / 29
24 Continuous/Discrete Approach for PD Continuous whitening g(t) continuous RF signal in a stationary case w(t) continuous uncorrelated unpredictable (scattering) Link with discrete: discretization in a spline basis Exponential splines: multiresolution version of L α Identification of continuous model from sampled data Shift-variant multiresolution description of scattering signal 9 / 29
25 Exponential Splines E-splines 2 : s(t) = k Z a kρ α (t t k ) + p 0 (t) ρ α green function of L α s(t) reconstructed in exponential B-spline basis: s(t) = k Z c[k]β α(t k) β α (t) = α {ρ α (t)} Scale T for generic sampling time: β αt (t) = αt {ρ α (t)} 2 [Khalidov and Unser, 2006], [Unser and Blu, 2005a], [Unser and Blu, 2005b] 10 / 29
26 Exponential Splines E-splines 2 : s(t) = k Z a kρ α (t t k ) + p 0 (t) ρ α green function of L α s(t) reconstructed in exponential B-spline basis: s(t) = k Z c[k]β α(t k) β α (t) = α {ρ α (t)} Scale T for generic sampling time: β αt (t) = αt {ρ α (t)} 2 [Khalidov and Unser, 2006], [Unser and Blu, 2005a], [Unser and Blu, 2005b] 10 / 29
27 Exponential Splines E-splines 2 : s(t) = k Z a kρ α (t t k ) + p 0 (t) ρ α green function of L α s(t) reconstructed in exponential B-spline basis: s(t) = k Z c[k]β α(t k) β α (t) = α {ρ α (t)} Scale T for generic sampling time: β αt (t) = αt {ρ α (t)} 2 [Khalidov and Unser, 2006], [Unser and Blu, 2005a], [Unser and Blu, 2005b] 10 / 29
28 E-Spline Multi-Resolution Analysis Dyadic scale for multiresolution analysis: T = 2 i Scaling function: ϕ i (t) = β 2 i(t)/ β 2 i L2 At scale i it is not a dilated version of the scaling function at scale 0. Mallat s filter bank algorithms but precomputation of filters for each scale Main property: wavelet coefficients of a signal f(t) are samples of smoothed version of L α (f) 11 / 29
29 E-Spline Multi-Resolution Analysis Dyadic scale for multiresolution analysis: T = 2 i Scaling function: ϕ i (t) = β 2 i(t)/ β 2 i L2 At scale i it is not a dilated version of the scaling function at scale 0. Mallat s filter bank algorithms but precomputation of filters for each scale Main property: wavelet coefficients of a signal f(t) are samples of smoothed version of L α (f) 11 / 29
30 12 / 29 E-spline-based Whitening Model Continuous: R gg (t) = (ρ α1 (t) ρ α1 ( t)) ( ρ αp (t) ρ αp ( t) ) Discrete: R(z) = z p B ( α: α) (z) p i=1 e αi z(1 e αi z)(1 e αi z 1 ) = b[1] p 1 i=1 1 ζ i (1 ζ i z)(1 ζ i z 1 ) p i=1 (1 e e αi αi z)(1 e αi z 1 ) Inter-dependence: ζ i = f( α) Whitening filter: W(z 1 ) = p 1 Y (1 ζ i z 1 ) i=1 py (1 e αi z 1 ) i=1 ARMA(p, p 1) discrete model (stable)
31 12 / 29 E-spline-based Whitening Model Continuous: R gg (t) = (ρ α1 (t) ρ α1 ( t)) ( ρ αp (t) ρ αp ( t) ) Discrete: R(z) = z p B ( α: α) (z) p i=1 e αi z(1 e αi z)(1 e αi z 1 ) = b[1] p 1 i=1 1 ζ i (1 ζ i z)(1 ζ i z 1 ) p i=1 (1 e e αi αi z)(1 e αi z 1 ) Inter-dependence: ζ i = f( α) Whitening filter: W(z 1 ) = p 1 Y (1 ζ i z 1 ) i=1 py (1 e αi z 1 ) i=1 ARMA(p, p 1) discrete model (stable)
32 12 / 29 E-spline-based Whitening Model Continuous: R gg (t) = (ρ α1 (t) ρ α1 ( t)) ( ρ αp (t) ρ αp ( t) ) Discrete: R(z) = z p B ( α: α) (z) p i=1 e αi z(1 e αi z)(1 e αi z 1 ) = b[1] p 1 i=1 1 ζ i (1 ζ i z)(1 ζ i z 1 ) p i=1 (1 e e αi αi z)(1 e αi z 1 ) Inter-dependence: ζ i = f( α) Whitening filter: W(z 1 ) = p 1 Y (1 ζ i z 1 ) i=1 py (1 e αi z 1 ) i=1 ARMA(p, p 1) discrete model (stable)
33 12 / 29 E-spline-based Whitening Model Continuous: R gg (t) = (ρ α1 (t) ρ α1 ( t)) ( ρ αp (t) ρ αp ( t) ) Discrete: R(z) = z p B ( α: α) (z) p i=1 e αi z(1 e αi z)(1 e αi z 1 ) = b[1] p 1 i=1 1 ζ i (1 ζ i z)(1 ζ i z 1 ) p i=1 (1 e e αi αi z)(1 e αi z 1 ) Inter-dependence: ζ i = f( α) Whitening filter: W(z 1 ) = p 1 Y (1 ζ i z 1 ) i=1 py (1 e αi z 1 ) i=1 ARMA(p, p 1) discrete model (stable)
34 12 / 29 E-spline-based Whitening Model Continuous: R gg (t) = (ρ α1 (t) ρ α1 ( t)) ( ρ αp (t) ρ αp ( t) ) Discrete: R(z) = z p B ( α: α) (z) p i=1 e αi z(1 e αi z)(1 e αi z 1 ) = b[1] p 1 i=1 1 ζ i (1 ζ i z)(1 ζ i z 1 ) p i=1 (1 e e αi αi z)(1 e αi z 1 ) Inter-dependence: ζ i = f( α) Whitening filter: W(z 1 ) = p 1 Y (1 ζ i z 1 ) i=1 py (1 e αi z 1 ) i=1 ARMA(p, p 1) discrete model (stable)
35 13 / 29 Exponential Parameters Estimation MA correction: AR(p): Yule-Walker equations. Only poles a i = e α i ARMA(p,p-1): No pole-zero inter-dependence. Annihilating Polynomial: coeff of AP of R gg (k) are functions { of α: p ( i=1 1 ai z 1)} R gg (k) = 0 AR to compute poles pole-zero inter-dependence to estimate zeros 1 MA(z) filtering to cancel zeros new AR computation for better pole estimation possibly iterations
36 13 / 29 Exponential Parameters Estimation MA correction: AR(p): Yule-Walker equations. Only poles a i = e α i ARMA(p,p-1): No pole-zero inter-dependence. Annihilating Polynomial: coeff of AP of R gg (k) are functions { of α: p ( i=1 1 ai z 1)} R gg (k) = 0 AR to compute poles pole-zero inter-dependence to estimate zeros 1 MA(z) filtering to cancel zeros new AR computation for better pole estimation possibly iterations
37 13 / 29 Exponential Parameters Estimation MA correction: AR(p): Yule-Walker equations. Only poles a i = e α i ARMA(p,p-1): No pole-zero inter-dependence. Annihilating Polynomial: coeff of AP of R gg (k) are functions { of α: p ( i=1 1 ai z 1)} R gg (k) = 0 AR to compute poles pole-zero inter-dependence to estimate zeros 1 MA(z) filtering to cancel zeros new AR computation for better pole estimation possibly iterations
38 13 / 29 Exponential Parameters Estimation MA correction: AR(p): Yule-Walker equations. Only poles a i = e α i ARMA(p,p-1): No pole-zero inter-dependence. Annihilating Polynomial: coeff of AP of R gg (k) are functions { of α: p ( i=1 1 ai z 1)} R gg (k) = 0 AR to compute poles pole-zero inter-dependence to estimate zeros 1 MA(z) filtering to cancel zeros new AR computation for better pole estimation possibly iterations
39 13 / 29 Exponential Parameters Estimation MA correction: AR(p): Yule-Walker equations. Only poles a i = e α i ARMA(p,p-1): No pole-zero inter-dependence. Annihilating Polynomial: coeff of AP of R gg (k) are functions { of α: p ( i=1 1 ai z 1)} R gg (k) = 0 AR to compute poles pole-zero inter-dependence to estimate zeros 1 MA(z) filtering to cancel zeros new AR computation for better pole estimation possibly iterations
40 14 / 29 Restoring Non Stationarity Non stationary RF signal Shift-variant whitening to get the non stationary scattering Two possibilities: Sliding window Recursive parameter estimation Recursive Annihilating Polynomial vs LS solution Time dependent parameters: α[n] Shift-variant E-spline wavelet Shift-variant multiresolution analysis (sliding window) Piece-wise multiresolution description of scattering
41 14 / 29 Restoring Non Stationarity Non stationary RF signal Shift-variant whitening to get the non stationary scattering Two possibilities: Sliding window Recursive parameter estimation Recursive Annihilating Polynomial vs LS solution Time dependent parameters: α[n] Shift-variant E-spline wavelet Shift-variant multiresolution analysis (sliding window) Piece-wise multiresolution description of scattering
42 14 / 29 Restoring Non Stationarity Non stationary RF signal Shift-variant whitening to get the non stationary scattering Two possibilities: Sliding window Recursive parameter estimation Recursive Annihilating Polynomial vs LS solution Time dependent parameters: α[n] Shift-variant E-spline wavelet Shift-variant multiresolution analysis (sliding window) Piece-wise multiresolution description of scattering
43 14 / 29 Restoring Non Stationarity Non stationary RF signal Shift-variant whitening to get the non stationary scattering Two possibilities: Sliding window Recursive parameter estimation Recursive Annihilating Polynomial vs LS solution Time dependent parameters: α[n] Shift-variant E-spline wavelet Shift-variant multiresolution analysis (sliding window) Piece-wise multiresolution description of scattering
44 15 / 29 Phantoms Targets Hypo and hyper echoic target detection Linear classifier: learning on ±9 db, testing on ±6 db Comparison with traditional predictive Deconvolution: improvement 40.3% 1 ROC curves: hypoechoic target detection ESW 1 feat PD Sensitivity Specificity
45 15 / 29 Phantoms Targets Hypo and hyper echoic target detection Linear classifier: learning on ±9 db, testing on ±6 db Comparison with traditional predictive Deconvolution: improvement 2.1% 1 ROC curves: hyperechoic target detection ESW 1 feat PD Sensitivity Specificity
46 16 / 29 Real Data Real data cancer detection: prostate trans-rectal Dataset: 15 benignant cases, 22 malignant cases Linear classifier: training 18 cases, testing 19 unknown img Comparison with traditional predictive Deconvolution: improvement 16.9%
47 17 / 29 Conclusions Predictive Deconvolution to restore US tissue response Continuous approach to predictive deconvolution Discretization in Exponential spline basis Identifications of E-splines tuned on US signal Multiresolution description of tissue response Improvement in detection performance: 16.9% Further Developments Improve parameter estitmation Better MA correction algorithm E-spline for texture analysis
48 17 / 29 Conclusions Predictive Deconvolution to restore US tissue response Continuous approach to predictive deconvolution Discretization in Exponential spline basis Identifications of E-splines tuned on US signal Multiresolution description of tissue response Improvement in detection performance: 16.9% Further Developments Improve parameter estitmation Better MA correction algorithm E-spline for texture analysis
49 18 / 29 Thank you for your attention!
50 Bibliography Jensen, J. (1994). Estimation of in vivo pulses in medical ultrasound. Ultrasonic Imaging, 16: Khalidov, I. and Unser, M. (2006). From differential equations to the constructio of new wavelet-like bases. IEEE Transactions on Signal processing, 54(4). Michailovich, O. V. and Adam, D. (2005). A novel approach to the 2-d blind deconvolution problem in medical ultrasound. IEEE Transactions on Medical Imaging, 24(1): Ng, J., Prager, R., Kingsbury, N., Treece, G., and Gee, A. (2007). Wavelet restoration of medical pulse-echo ultrasound images in an em framework. IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, 54(3): Rasmussen, K. (1994). Maximum likelihood estimation of the attenuated uktrasound pulse. IEEE Transactions on Signal Processing, 42: Unser, M. and Blu, T. (2005a). Cardinal exponential splines: Part i - theory and filtering algorithms. IEEE Transactions on Signal Processing, 5373(4). Unser, M. and Blu, T. (2005b). Cardinal exponential splines: Part ii - think analog, act digital. IEEE Transactions on Signal Processing, 53(4). 19 / 29
51 20 / 29 Publications 2009 (to be published in next issue) IEEE Transactions on Medical Imaging S. Maggio, A. Palladini, L. De Marchi, M. Alessandrini, N. Speciale, G. Masetti Predictive deconvolution and hybrid feature Selection for Computer-Aided Detection of prostate cancer 2009 March Proceedings International Symposium on Acoustical Imaging M. Scebran, A. Palladini, S. Maggio, L. De Marchi, N. Speciale Automatic regions of interests segmentation for computer aided classification of prostate TRUS images 2008 November Proceedings IEEE IUS2008 S. Maggio, L. De Marchi, M. Alessandrini, N. Speciale Computer aided detection of prostate cancer based on GDA and predictive deconvolution 2005 November WSEAS Transactions on Systems S. Maggio, N. Testoni, L. De Marchi, N. Speciale, G. Masetti Ultrasound Images Enhancement by means of Deconvolution Algorithms in the Wavelet Domain 2005 September WSEAS Int. Conf. on Signal Processing, Computational Geometry and Artificial Vision (ISCGAV) S. Maggio, N. Testoni, L. De Marchi, N. Speciale, G. Masetti Wavelet-based Deconvolution Algorithms Applied to Ultrasound Images
52 21 / 29 Texture analysis through E-splines Global learning of exponential parameters E-spline wavelet tuned on ultrasound signal E-spline wavelet transform for texture typing Comparison with traditional wavelets Improvement with respect to energy information: 14.9% 1 ROC curves: cancer detection Nakagami Naka + Variance ESW Sensitivity Specificity
53 Comparison with previous works Table: Published methods for ultrasound-based prostate tissue characterization Work Ground Truth Technique Results % # ROIs Features SE SP Acc Az Basset 37 Textural Huynen - Textural Houston 25 Textural Schmitz 3405 Multi Scheipers Multi Feleppa 1019 Spectral Mohamed 96 Textural Llobet 4944 Textural Mohamed 108 Multi Han 2000 Multi Maggio Multi on RLS / 29
54 Autocorrelation Function Discretization 1/4 In the first order case the ACF of g(t) can be obtained as: R gg (τ) = F 1 {Φ g } = F 1 { 1 jω α } F 1 { } 1 jω α = e αt u(t) e αt u( t) = ρ α (t) ρ α ( t) (1) + Green function of L α : ρ α (t) = e αt u(t) = p α [k]β α (t k) R gg (t) = = = + k=0 + k=0 + k=0 p α [k]β α (t k) + p α [k] p α [k] k =0 + k =0 + k =0 k=0 p α [k ]β α ( t k ) p α [k ]β α (t k) β α ( t k ) p α [k ]e α β α (t k) β α (t + k + 1) (2) 23 / 29
55 24 / 29 Autocorrelation Function Discretization 2/4 The expression for ACF simplifies to R gg (t) = + k=0 p α [k] p α [k k 1]e α β ( α:α) (t k ) k Z = (p α p α)[k + 1]e α β ( α:α) (t k ) (3) k Z where p α [k] = p α[ k]. It turns out that the z-transform of the discrete signal ACF is given as: R(z) = e α zp α (z)p α (z)b ( α:α)(z) e α zb ( α:α) (z) = (1 e α z 1 )(1 e α z) (4)
56 25 / 29 Autocorrelation Function Discretization 3/4 The extension to higher-order AR models follows naturally as shown below: { R gg (t) = F 1 1 } { 2 jω α 1 F 1 1 } 2 jω α p = (ρ α1 (t) ρ α1 ( t)) ( ρ αp (t) ρ αp ( t) ) = k1 Z(p α1 p α 1 )[k 1 + 1]e α 1 β ( α1 :α 1 )(t k 1 ) kp Z(p αp p α p )[k p + 1]e αp β ( αp:α p)(t k p ) (5)
57 26 / 29 Autocorrelation Function Discretization 4/4 The final expression of the z-transform of R gg (k) is: R(z) = = p e α i zp αi (z)p α i (z) B ( α: α) (z) i=1 p i=1 e α i z (1 e α i z)(1 e α iz 1 ) zp B ( α: α) (z) (6)
58 27 / 29 Whitening Filter Stability 1/3 The proof that the z-transform of B-spline β ( α: α), B ( α, α) (z), doesn t have zeros on unit circle is shown by the following considerations: ˆβ ( α, α) (ω) = = = F [β α β α (t)] p e αi F [β α β α (t p)] i=1 p e αi e jωp ˆβ α (ω) 2 (7) i=1
59 28 / 29 Whitening Filter Stability 2/3 B ( α, α) (z) is linked to the Fourier transform of samples of β ( α, α) (t), ˆβ s,( α, α) (ω): ˆβ s,( α, α) (ω) = k ˆβ ( α, α) (ω + 2πk) = = p e αi e jωp ˆβ α (ω + 2πk) 2 k p e αi e jωp A α (e jω ) (8) i=1 i=1
60 29 / 29 Whitening Filter Stability 3/3 where A α (e jω ) is the discrete Fourier transform of the Gram sequence of B-splines, and, for the Riesz basis property, it is always grater than zero: 0 < r 2 α < A α(e jω ) < R 2 α < + (9) As a consequence ˆβ s,( α, α) (ω) > 0 for every ω and, since B ( α, α) (z) = ˆβ s,( α, α) (ω) e jω =z, B ( α, α) (z) doesn t have any zeros on the unit circle. This proof and the consideration about reciprocal roots of B ( α, α) (z) guarantee the whitening filter stability.
Parametric Signal Modeling and Linear Prediction Theory 1. Discrete-time Stochastic Processes (cont d)
Parametric Signal Modeling and Linear Prediction Theory 1. Discrete-time Stochastic Processes (cont d) Electrical & Computer Engineering North Carolina State University Acknowledgment: ECE792-41 slides
More informationAN INTRODUCTION TO WAVELET TRANSFORMS
AN INTRODUCTION TO WAVELET TRANSFORMS Luca De Marchi 7/9/005 1 OUTLINE Time-Frequency Analysis Introduction on Wavelet Operators Examples of applications: Radar/Sonar Experimental results Conclusions 7/9/005
More informationSENSITIVITY TO POINT-SPREAD FUNCTION PARAMETERS IN MEDICAL ULTRASOUND IMAGE DECONVOLUTION
SENSITIVITY TO POINT-SPREAD FUNCTION PARAMETERS IN MEDICAL ULTRASOUND IMAGE DECONVOLUTION H.-C. Shin, R. W. Prager, J. K. H. Ng, W. H. Gomersall, N. G. Kingsbury, G. M. Treece and A. H. Gee CUED / F-INFENG
More informationBlind Deconvolution of Ultrasonic Signals Using High-Order Spectral Analysis and Wavelets
Blind Deconvolution of Ultrasonic Signals Using High-Order Spectral Analysis and Wavelets Roberto H. Herrera, Eduardo Moreno, Héctor Calas, and Rubén Orozco 3 University of Cienfuegos, Cuatro Caminos,
More informationA time series is called strictly stationary if the joint distribution of every collection (Y t
5 Time series A time series is a set of observations recorded over time. You can think for example at the GDP of a country over the years (or quarters) or the hourly measurements of temperature over a
More informationOn Moving Average Parameter Estimation
On Moving Average Parameter Estimation Niclas Sandgren and Petre Stoica Contact information: niclas.sandgren@it.uu.se, tel: +46 8 473392 Abstract Estimation of the autoregressive moving average (ARMA)
More information/96$ IEEE
IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL, VOL. 43, NO. 5, SEPTEMBER 1996 979 Correspondence Order Selection Criteria for Detecting Mean Scatterer Spacings with the AR Model
More information31545 Medical Imaging systems
Simulation of ultrasound systems and non-linear imaging 545 Medical Imaging systems Lecture 9: Simulation of ultrasound systems and non-linear imaging Jørgen Arendt Jensen Department of Electrical Engineering
More informationMultiresolution analysis & wavelets (quick tutorial)
Multiresolution analysis & wavelets (quick tutorial) Application : image modeling André Jalobeanu Multiresolution analysis Set of closed nested subspaces of j = scale, resolution = 2 -j (dyadic wavelets)
More informationRepresentation: Fractional Splines, Wavelets and related Basis Function Expansions. Felix Herrmann and Jonathan Kane, ERL-MIT
Representation: Fractional Splines, Wavelets and related Basis Function Expansions Felix Herrmann and Jonathan Kane, ERL-MIT Objective: Build a representation with a regularity that is consistent with
More informationGENERALIZED DEFLATION ALGORITHMS FOR THE BLIND SOURCE-FACTOR SEPARATION OF MIMO-FIR CHANNELS. Mitsuru Kawamoto 1,2 and Yujiro Inouye 1
GENERALIZED DEFLATION ALGORITHMS FOR THE BLIND SOURCE-FACTOR SEPARATION OF MIMO-FIR CHANNELS Mitsuru Kawamoto,2 and Yuiro Inouye. Dept. of Electronic and Control Systems Engineering, Shimane University,
More informationCorrelation analysis of three-dimensional strain imaging using ultrasound two-dimensional array transducers
Correlation analysis of three-dimensional strain imaging using ultrasound two-dimensional array transducers Min Rao a and Tomy Varghese Department of Medical Physics, The University of Wisconsin-Madison,
More informationADAPTIVE FILTER THEORY
ADAPTIVE FILTER THEORY Fourth Edition Simon Haykin Communications Research Laboratory McMaster University Hamilton, Ontario, Canada Front ice Hall PRENTICE HALL Upper Saddle River, New Jersey 07458 Preface
More informationNOISE ROBUST RELATIVE TRANSFER FUNCTION ESTIMATION. M. Schwab, P. Noll, and T. Sikora. Technical University Berlin, Germany Communication System Group
NOISE ROBUST RELATIVE TRANSFER FUNCTION ESTIMATION M. Schwab, P. Noll, and T. Sikora Technical University Berlin, Germany Communication System Group Einsteinufer 17, 1557 Berlin (Germany) {schwab noll
More informationSymmetric Wavelet Tight Frames with Two Generators
Symmetric Wavelet Tight Frames with Two Generators Ivan W. Selesnick Electrical and Computer Engineering Polytechnic University 6 Metrotech Center, Brooklyn, NY 11201, USA tel: 718 260-3416, fax: 718 260-3906
More informationIndex. p, lip, 78 8 function, 107 v, 7-8 w, 7-8 i,7-8 sine, 43 Bo,94-96
p, lip, 78 8 function, 107 v, 7-8 w, 7-8 i,7-8 sine, 43 Bo,94-96 B 1,94-96 M,94-96 B oro!' 94-96 BIro!' 94-96 I/r, 79 2D linear system, 56 2D FFT, 119 2D Fourier transform, 1, 12, 18,91 2D sinc, 107, 112
More informationDiscrete Fourier Transform
Discrete Fourier Transform Virtually all practical signals have finite length (e.g., sensor data, audio records, digital images, stock values, etc). Rather than considering such signals to be zero-padded
More informationarxiv: v1 [cs.cv] 10 Feb 2016
GABOR WAVELETS IN IMAGE PROCESSING David Bařina Doctoral Degree Programme (2), FIT BUT E-mail: xbarin2@stud.fit.vutbr.cz Supervised by: Pavel Zemčík E-mail: zemcik@fit.vutbr.cz arxiv:162.338v1 [cs.cv]
More informationBlind Deconvolution via Maximum Kurtosis Adaptive Filtering
Blind Deconvolution via Maximum Kurtosis Adaptive Filtering Deborah Pereg Doron Benzvi The Jerusalem College of Engineering Jerusalem, Israel doronb@jce.ac.il, deborahpe@post.jce.ac.il ABSTRACT In this
More informationDigital Image Processing Lectures 13 & 14
Lectures 13 & 14, Professor Department of Electrical and Computer Engineering Colorado State University Spring 2013 Properties of KL Transform The KL transform has many desirable properties which makes
More informationBasics on 2-D 2 D Random Signal
Basics on -D D Random Signal Spring 06 Instructor: K. J. Ray Liu ECE Department, Univ. of Maryland, College Park Overview Last Time: Fourier Analysis for -D signals Image enhancement via spatial filtering
More information550 ieee transactions on ultrasonics, ferroelectrics, and frequency control, vol. 54, no. 3, march 2007
55 ieee transactions on ultrasonics, ferroelectrics, and frequency control, vol. 54, no. 3, march 7 Wavelet Restoration of Medical Pulse-Echo Ultrasound Images in an EM Framework James Ng, Richard Prager,
More informationUltrasound is a widely used, nonionizing, and costeffective
340 IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 57, no. 2, February 2010 Trade-Offs in Data Acquisition and Processing Parameters for Backscatter and Scatterer Size Estimations
More informationTime Series Examples Sheet
Lent Term 2001 Richard Weber Time Series Examples Sheet This is the examples sheet for the M. Phil. course in Time Series. A copy can be found at: http://www.statslab.cam.ac.uk/~rrw1/timeseries/ Throughout,
More informationParametric Method Based PSD Estimation using Gaussian Window
International Journal of Engineering Trends and Technology (IJETT) Volume 29 Number 1 - November 215 Parametric Method Based PSD Estimation using Gaussian Window Pragati Sheel 1, Dr. Rajesh Mehra 2, Preeti
More informationAdaptive Filter Theory
0 Adaptive Filter heory Sung Ho Cho Hanyang University Seoul, Korea (Office) +8--0-0390 (Mobile) +8-10-541-5178 dragon@hanyang.ac.kr able of Contents 1 Wiener Filters Gradient Search by Steepest Descent
More informationStatistical and Adaptive Signal Processing
r Statistical and Adaptive Signal Processing Spectral Estimation, Signal Modeling, Adaptive Filtering and Array Processing Dimitris G. Manolakis Massachusetts Institute of Technology Lincoln Laboratory
More informationLesson 9: Autoregressive-Moving Average (ARMA) models
Lesson 9: Autoregressive-Moving Average (ARMA) models Dipartimento di Ingegneria e Scienze dell Informazione e Matematica Università dell Aquila, umberto.triacca@ec.univaq.it Introduction We have seen
More informationLesson 13: Box-Jenkins Modeling Strategy for building ARMA models
Lesson 13: Box-Jenkins Modeling Strategy for building ARMA models Facoltà di Economia Università dell Aquila umberto.triacca@gmail.com Introduction In this lesson we present a method to construct an ARMA(p,
More informationSatellite image deconvolution using complex wavelet packets
Satellite image deconvolution using complex wavelet packets André Jalobeanu, Laure Blanc-Féraud, Josiane Zerubia ARIANA research group INRIA Sophia Antipolis, France CNRS / INRIA / UNSA www.inria.fr/ariana
More informationA SARIMAX coupled modelling applied to individual load curves intraday forecasting
A SARIMAX coupled modelling applied to individual load curves intraday forecasting Frédéric Proïa Workshop EDF Institut Henri Poincaré - Paris 05 avril 2012 INRIA Bordeaux Sud-Ouest Institut de Mathématiques
More informationDigital Signal Processing Lecture 10 - Discrete Fourier Transform
Digital Signal Processing - Discrete Fourier Transform Electrical Engineering and Computer Science University of Tennessee, Knoxville November 12, 2015 Overview 1 2 3 4 Review - 1 Introduction Discrete-time
More informationShift-Variance and Nonstationarity of Generalized Sampling-Reconstruction Processes
Shift-Variance and Nonstationarity of Generalized Sampling-Reconstruction Processes Runyi Yu Eastern Mediterranean University Gazimagusa, North Cyprus Web: faraday.ee.emu.edu.tr/yu Emails: runyi.yu@emu.edu.tr
More informationWhich wavelet bases are the best for image denoising?
Which wavelet bases are the best for image denoising? Florian Luisier a, Thierry Blu a, Brigitte Forster b and Michael Unser a a Biomedical Imaging Group (BIG), Ecole Polytechnique Fédérale de Lausanne
More informationEXPERIMENTAL VALIDATION OF THE SPECTRAL FIT ALGORITHM USING TISSUE MIMICKING PHANTOMS
EXPERIMENTAL VALIDATION OF THE SPECTRAL FIT ALGORITHM USING TISSUE MIMICKING PHANTOMS T.A. Bigelow, W.D. O Brien, Jr. Bioacoustics Research Laboratory, Department of Electrical and Computer Engineering,
More informationAngular Spectrum Decomposition Analysis of Second Harmonic Ultrasound Propagation and its Relation to Tissue Harmonic Imaging
The 4 th International Workshop on Ultrasonic and Advanced Methods for Nondestructive Testing and Material Characterization, June 9, 006 at ABSTRACT Angular Spectrum Decomposition Analysis of Second Harmonic
More informationWavelet Based Image Restoration Using Cross-Band Operators
1 Wavelet Based Image Restoration Using Cross-Band Operators Erez Cohen Electrical Engineering Department Technion - Israel Institute of Technology Supervised by Prof. Israel Cohen 2 Layout Introduction
More informationTime Series: Theory and Methods
Peter J. Brockwell Richard A. Davis Time Series: Theory and Methods Second Edition With 124 Illustrations Springer Contents Preface to the Second Edition Preface to the First Edition vn ix CHAPTER 1 Stationary
More informationSampling Signals of Finite Rate of Innovation*
Sampling Signals of Finite Rate of Innovation* Martin Vetterli http://lcavwww.epfl.ch/~vetterli EPFL & UCBerkeley.15.1.1.5.5.5.5.1.1.15.2 128 128 256 384 512 64 768 896 124 1152.15 128 128 256 384 512
More informationApplied Time. Series Analysis. Wayne A. Woodward. Henry L. Gray. Alan C. Elliott. Dallas, Texas, USA
Applied Time Series Analysis Wayne A. Woodward Southern Methodist University Dallas, Texas, USA Henry L. Gray Southern Methodist University Dallas, Texas, USA Alan C. Elliott University of Texas Southwestern
More informationContents. Digital Signal Processing, Part II: Power Spectrum Estimation
Contents Digital Signal Processing, Part II: Power Spectrum Estimation 5. Application of the FFT for 7. Parametric Spectrum Est. Filtering and Spectrum Estimation 7.1 ARMA-Models 5.1 Fast Convolution 7.2
More informationWavelets and multiresolution representations. Time meets frequency
Wavelets and multiresolution representations Time meets frequency Time-Frequency resolution Depends on the time-frequency spread of the wavelet atoms Assuming that ψ is centred in t=0 Signal domain + t
More informationPower Supply Quality Analysis Using S-Transform and SVM Classifier
Journal of Power and Energy Engineering, 2014, 2, 438-447 Published Online April 2014 in SciRes. http://www.scirp.org/journal/jpee http://dx.doi.org/10.4236/jpee.2014.24059 Power Supply Quality Analysis
More informationMoving Average (MA) representations
Moving Average (MA) representations The moving average representation of order M has the following form v[k] = MX c n e[k n]+e[k] (16) n=1 whose transfer function operator form is MX v[k] =H(q 1 )e[k],
More informationDeep Learning: Approximation of Functions by Composition
Deep Learning: Approximation of Functions by Composition Zuowei Shen Department of Mathematics National University of Singapore Outline 1 A brief introduction of approximation theory 2 Deep learning: approximation
More informationWavelets and differential operators: from fractals to Marr!s primal sketch
Wavelets and differential operators: from fractals to Marrs primal sketch Michael Unser Biomedical Imaging Group EPFL, Lausanne, Switzerland Joint work with Pouya Tafti and Dimitri Van De Ville Plenary
More informationEEG- Signal Processing
Fatemeh Hadaeghi EEG- Signal Processing Lecture Notes for BSP, Chapter 5 Master Program Data Engineering 1 5 Introduction The complex patterns of neural activity, both in presence and absence of external
More informationTopic 7: Filter types and structures
ELEN E481: Digital Signal Processing Topic 7: Filter tpes and structures 1. Some filter tpes 2. Minimum and maximum phase 3. Filter implementation structures 1 1. Some Filter Tpes We have seen the basics
More informationIntroduction to Biomedical Engineering
Introduction to Biomedical Engineering Biosignal processing Kung-Bin Sung 6/11/2007 1 Outline Chapter 10: Biosignal processing Characteristics of biosignals Frequency domain representation and analysis
More informationMachine Learning. A Bayesian and Optimization Perspective. Academic Press, Sergios Theodoridis 1. of Athens, Athens, Greece.
Machine Learning A Bayesian and Optimization Perspective Academic Press, 2015 Sergios Theodoridis 1 1 Dept. of Informatics and Telecommunications, National and Kapodistrian University of Athens, Athens,
More informationECE6604 PERSONAL & MOBILE COMMUNICATIONS. Week 3. Flat Fading Channels Envelope Distribution Autocorrelation of a Random Process
1 ECE6604 PERSONAL & MOBILE COMMUNICATIONS Week 3 Flat Fading Channels Envelope Distribution Autocorrelation of a Random Process 2 Multipath-Fading Mechanism local scatterers mobile subscriber base station
More informationMultirate signal processing
Multirate signal processing Discrete-time systems with different sampling rates at various parts of the system are called multirate systems. The need for such systems arises in many applications, including
More informationComputer Vision & Digital Image Processing
Computer Vision & Digital Image Processing Image Restoration and Reconstruction I Dr. D. J. Jackson Lecture 11-1 Image restoration Restoration is an objective process that attempts to recover an image
More informationComputer Exercise 0 Simulation of ARMA-processes
Lund University Time Series Analysis Mathematical Statistics Fall 2018 Centre for Mathematical Sciences Computer Exercise 0 Simulation of ARMA-processes The purpose of this computer exercise is to illustrate
More informationMultichannel Sampling of Signals with Finite. Rate of Innovation
Multichannel Sampling of Signals with Finite 1 Rate of Innovation Hojjat Akhondi Asl, Pier Luigi Dragotti and Loic Baboulaz EDICS: DSP-TFSR, DSP-SAMP. Abstract In this paper we present a possible extension
More informationA Subspace Approach to Estimation of. Measurements 1. Carlos E. Davila. Electrical Engineering Department, Southern Methodist University
EDICS category SP 1 A Subspace Approach to Estimation of Autoregressive Parameters From Noisy Measurements 1 Carlos E Davila Electrical Engineering Department, Southern Methodist University Dallas, Texas
More informationFast simulation of nonlinear radio frequency ultrasound images in inhomogeneous nonlinear media: CREANUIS
Proceedings of the Acoustics 2012 Nantes Conference Fast simulation of nonlinear radio frequency ultrasound images in inhomogeneous nonlinear media: CREANUIS F. Varray a, M. Toulemonde a,b, O. Basset a
More informationOptimal Beamforming in Ultrasound Using the Ideal Observer. Craig K. Abbey, Nghia Q. Nguyen, and Michael F. Insana, Senior Member, IEEE
178 IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 57, no. 8, August 010 Optimal Beamforming in Ultrasound Using the Ideal Observer Craig K. Abbey, Nghia Q. Nguyen, and Michael
More information64 ieee transactions on ultrasonics, ferroelectrics, and frequency control, vol. 53, no. 1, january Despeckling of Medical Ultrasound Images
64 ieee transactions on ultrasonics, ferroelectrics, and frequency control, vol. 53, no. 1, january 2006 Despeckling of Medical Ultrasound Images Oleg V. Michailovich and Allen Tannenbaum, Member, IEEE
More informationNonparametric and Parametric Defined This text distinguishes between systems and the sequences (processes) that result when a WN input is applied
Linear Signal Models Overview Introduction Linear nonparametric vs. parametric models Equivalent representations Spectral flatness measure PZ vs. ARMA models Wold decomposition Introduction Many researchers
More informationTime Series Examples Sheet
Lent Term 2001 Richard Weber Time Series Examples Sheet This is the examples sheet for the M. Phil. course in Time Series. A copy can be found at: http://www.statslab.cam.ac.uk/~rrw1/timeseries/ Throughout,
More informationMarcel Dettling. Applied Time Series Analysis SS 2013 Week 05. ETH Zürich, March 18, Institute for Data Analysis and Process Design
Marcel Dettling Institute for Data Analysis and Process Design Zurich University of Applied Sciences marcel.dettling@zhaw.ch http://stat.ethz.ch/~dettling ETH Zürich, March 18, 2013 1 Basics of Modeling
More informationForecasting. Simon Shaw 2005/06 Semester II
Forecasting Simon Shaw s.c.shaw@maths.bath.ac.uk 2005/06 Semester II 1 Introduction A critical aspect of managing any business is planning for the future. events is called forecasting. Predicting future
More informationA Hybrid Time-delay Prediction Method for Networked Control System
International Journal of Automation and Computing 11(1), February 2014, 19-24 DOI: 10.1007/s11633-014-0761-1 A Hybrid Time-delay Prediction Method for Networked Control System Zhong-Da Tian Xian-Wen Gao
More informationComputer Exercise 1 Estimation and Model Validation
Lund University Time Series Analysis Mathematical Statistics Fall 2018 Centre for Mathematical Sciences Computer Exercise 1 Estimation and Model Validation This computer exercise treats identification,
More informationESE 531: Digital Signal Processing
ESE 531: Digital Signal Processing Lec 6: January 30, 2018 Inverse z-transform Lecture Outline! z-transform " Tie up loose ends " Regions of convergence properties! Inverse z-transform " Inspection " Partial
More informationCCNY. BME I5100: Biomedical Signal Processing. Stochastic Processes. Lucas C. Parra Biomedical Engineering Department City College of New York
BME I5100: Biomedical Signal Processing Stochastic Processes Lucas C. Parra Biomedical Engineering Department CCNY 1 Schedule Week 1: Introduction Linear, stationary, normal - the stuff biology is not
More informationCentre for Mathematical Sciences HT 2017 Mathematical Statistics. Study chapters 6.1, 6.2 and in the course book.
Lund University Stationary stochastic processes Centre for Mathematical Sciences HT 2017 Mathematical Statistics Computer exercise 2 in Stationary stochastic processes, HT 17. The purpose with this computer
More informationChapter 6: Model Specification for Time Series
Chapter 6: Model Specification for Time Series The ARIMA(p, d, q) class of models as a broad class can describe many real time series. Model specification for ARIMA(p, d, q) models involves 1. Choosing
More informationMedical Image Processing
Medical Image Processing Federica Caselli Department of Civil Engineering University of Rome Tor Vergata Medical Imaging X-Ray CT Ultrasound MRI PET/SPECT Digital Imaging! Medical Image Processing What
More informationDigital Signal Processing. Midterm 2 Solutions
EE 123 University of California, Berkeley Anant Sahai arch 15, 2007 Digital Signal Processing Instructions idterm 2 Solutions Total time allowed for the exam is 80 minutes Please write your name and SID
More informationLinear Prediction 1 / 41
Linear Prediction 1 / 41 A map of speech signal processing Natural signals Models Artificial signals Inference Speech synthesis Hidden Markov Inference Homomorphic processing Dereverberation, Deconvolution
More informationBLIND DECONVOLUTION ALGORITHMS FOR MIMO-FIR SYSTEMS DRIVEN BY FOURTH-ORDER COLORED SIGNALS
BLIND DECONVOLUTION ALGORITHMS FOR MIMO-FIR SYSTEMS DRIVEN BY FOURTH-ORDER COLORED SIGNALS M. Kawamoto 1,2, Y. Inouye 1, A. Mansour 2, and R.-W. Liu 3 1. Department of Electronic and Control Systems Engineering,
More informationDiscrete-time signals and systems
Discrete-time signals and systems 1 DISCRETE-TIME DYNAMICAL SYSTEMS x(t) G y(t) Linear system: Output y(n) is a linear function of the inputs sequence: y(n) = k= h(k)x(n k) h(k): impulse response of the
More information2A1H Time-Frequency Analysis II Bugs/queries to HT 2011 For hints and answers visit dwm/courses/2tf
Time-Frequency Analysis II (HT 20) 2AH 2AH Time-Frequency Analysis II Bugs/queries to david.murray@eng.ox.ac.uk HT 20 For hints and answers visit www.robots.ox.ac.uk/ dwm/courses/2tf David Murray. A periodic
More information446 SCIENCE IN CHINA (Series F) Vol. 46 introduced in refs. [6, ]. Based on this inequality, we add normalization condition, symmetric conditions and
Vol. 46 No. 6 SCIENCE IN CHINA (Series F) December 003 Construction for a class of smooth wavelet tight frames PENG Lizhong (Λ Π) & WANG Haihui (Ξ ) LMAM, School of Mathematical Sciences, Peking University,
More informationPractical Spectral Estimation
Digital Signal Processing/F.G. Meyer Lecture 4 Copyright 2015 François G. Meyer. All Rights Reserved. Practical Spectral Estimation 1 Introduction The goal of spectral estimation is to estimate how the
More informationFilter Banks II. Prof. Dr.-Ing. G. Schuller. Fraunhofer IDMT & Ilmenau University of Technology Ilmenau, Germany
Filter Banks II Prof. Dr.-Ing. G. Schuller Fraunhofer IDMT & Ilmenau University of Technology Ilmenau, Germany Page Modulated Filter Banks Extending the DCT The DCT IV transform can be seen as modulated
More informationCourse content (will be adapted to the background knowledge of the class):
Biomedical Signal Processing and Signal Modeling Lucas C Parra, parra@ccny.cuny.edu Departamento the Fisica, UBA Synopsis This course introduces two fundamental concepts of signal processing: linear systems
More informationSparse Sampling: Theory and Applications
November 24, 29 Outline Problem Statement Signals with Finite Rate of Innovation Sampling Kernels: E-splines and -splines Sparse Sampling: the asic Set-up and Extensions The Noisy Scenario Applications
More informationNEW STEIGLITZ-McBRIDE ADAPTIVE LATTICE NOTCH FILTERS
NEW STEIGLITZ-McBRIDE ADAPTIVE LATTICE NOTCH FILTERS J.E. COUSSEAU, J.P. SCOPPA and P.D. DOÑATE CONICET- Departamento de Ingeniería Eléctrica y Computadoras Universidad Nacional del Sur Av. Alem 253, 8000
More informationCHAPTER 2 RANDOM PROCESSES IN DISCRETE TIME
CHAPTER 2 RANDOM PROCESSES IN DISCRETE TIME Shri Mata Vaishno Devi University, (SMVDU), 2013 Page 13 CHAPTER 2 RANDOM PROCESSES IN DISCRETE TIME When characterizing or modeling a random variable, estimates
More informationOptimized Compact-support Interpolation Kernels
Optimized Compact-support Interpolation Kernels Ramtin Madani, Student Member, IEEE, Ali Ayremlou, Student Member, IEEE, Arash Amini, Farrokh Marvasti, Senior Member, IEEE, Abstract In this paper, we investigate
More information3 Theory of stationary random processes
3 Theory of stationary random processes 3.1 Linear filters and the General linear process A filter is a transformation of one random sequence {U t } into another, {Y t }. A linear filter is a transformation
More informationEE538 Final Exam Fall :20 pm -5:20 pm PHYS 223 Dec. 17, Cover Sheet
EE538 Final Exam Fall 005 3:0 pm -5:0 pm PHYS 3 Dec. 17, 005 Cover Sheet Test Duration: 10 minutes. Open Book but Closed Notes. Calculators ARE allowed!! This test contains five problems. Each of the five
More information1 Random walks and data
Inference, Models and Simulation for Complex Systems CSCI 7-1 Lecture 7 15 September 11 Prof. Aaron Clauset 1 Random walks and data Supposeyou have some time-series data x 1,x,x 3,...,x T and you want
More informationImage Degradation Model (Linear/Additive)
Image Degradation Model (Linear/Additive),,,,,,,, g x y h x y f x y x y G uv H uv F uv N uv 1 Source of noise Image acquisition (digitization) Image transmission Spatial properties of noise Statistical
More informationI. INTRODUCTION J. Acoust. Soc. Am. 115 (6), June /2004/115(6)/3226/9/$ Acoustical Society of America
Defining optimal axial and lateral resolution for estimating scatterer properties from volumes using ultrasound backscatter Michael L. Oelze a) and William D. O Brien, Jr. Bioacoustics Research Laboratory,
More informationIntroduction to Signal Processing
to Signal Processing Davide Bacciu Dipartimento di Informatica Università di Pisa bacciu@di.unipi.it Intelligent Systems for Pattern Recognition Signals = Time series Definitions Motivations A sequence
More informationA Unified Formulation of Gaussian Versus Sparse Stochastic Processes
A Unified Formulation of Gaussian Versus Sparse Stochastic Processes Michael Unser, Pouya Tafti and Qiyu Sun EPFL, UCF Appears in IEEE Trans. on Information Theory, 2014 Presented by Liming Wang M. Unser,
More informationGENERALIZATION OF CAMPBELL S THEOREM TO NONSTATIONARY NOISE. Leon Cohen
GENERALIZATION OF CAMPBELL S THEOREM TO NONSTATIONARY NOISE Leon Cohen City University of New York, Hunter College, 695 Park Ave, New York NY 10056, USA ABSTRACT Campbell s theorem is a fundamental result
More informationLecture 4 - Spectral Estimation
Lecture 4 - Spectral Estimation The Discrete Fourier Transform The Discrete Fourier Transform (DFT) is the equivalent of the continuous Fourier Transform for signals known only at N instants separated
More informationDirect Learning: Linear Classification. Donglin Zeng, Department of Biostatistics, University of North Carolina
Direct Learning: Linear Classification Logistic regression models for classification problem We consider two class problem: Y {0, 1}. The Bayes rule for the classification is I(P(Y = 1 X = x) > 1/2) so
More informationSolutions. Number of Problems: 10
Final Exam February 9th, 2 Signals & Systems (5-575-) Prof. R. D Andrea Solutions Exam Duration: 5 minutes Number of Problems: Permitted aids: One double-sided A4 sheet. Questions can be answered in English
More informationChapter 4: Models for Stationary Time Series
Chapter 4: Models for Stationary Time Series Now we will introduce some useful parametric models for time series that are stationary processes. We begin by defining the General Linear Process. Let {Y t
More informationSpatially adaptive alpha-rooting in BM3D sharpening
Spatially adaptive alpha-rooting in BM3D sharpening Markku Mäkitalo and Alessandro Foi Department of Signal Processing, Tampere University of Technology, P.O. Box FIN-553, 33101, Tampere, Finland e-mail:
More informationADAPTIVE FILTER THEORY
ADAPTIVE FILTER THEORY Fifth Edition Simon Haykin Communications Research Laboratory McMaster University Hamilton, Ontario, Canada International Edition contributions by Telagarapu Prabhakar Department
More information5 Transfer function modelling
MSc Further Time Series Analysis 5 Transfer function modelling 5.1 The model Consider the construction of a model for a time series (Y t ) whose values are influenced by the earlier values of a series
More informationIntroduction to ARMA and GARCH processes
Introduction to ARMA and GARCH processes Fulvio Corsi SNS Pisa 3 March 2010 Fulvio Corsi Introduction to ARMA () and GARCH processes SNS Pisa 3 March 2010 1 / 24 Stationarity Strict stationarity: (X 1,
More informationMultiresolution schemes
Multiresolution schemes Fondamenti di elaborazione del segnale multi-dimensionale Stefano Ferrari Università degli Studi di Milano stefano.ferrari@unimi.it Elaborazione dei Segnali Multi-dimensionali e
More information