Model-based Correlation Measure for Gain and Offset Nonuniformity in Infrared Focal-Plane-Array Sensors
|
|
- Eugene Chandler
- 5 years ago
- Views:
Transcription
1 Model-based Correlation Measure for Gain and Offset Nonuniformity in Infrared Focal-Plane-Array Sensors César San Martin Sergio Torres Abstract In this paper, a model-based correlation measure between gain and offset nonuniformity in infrared focal plane array (FPA) imaging systems is developed. The correlation between the nonuniformity gain and offset is modeled by means of a generalized Gauss-Marov model. The propose model-based correlation is validated by using several real infrared video sequences and three well nown scene-based nonuniformity correction methods. Keywords: Gauss-Marov Model, Image Sequence Processing, Infrared FPA, Signal Processing 1. INTRODUCTION Infrared (IR) cameras use an IR sensor to digitize the information, and due to its high performance, the most employed integrated technology in IR sensors is the Focal Plane Array (FPA). An IR-FPA is a die composed of a group of photodetectors placed in a focal plane forming a matrix of X Y pixels, which gives the sensor the ability to collect the IR radiation. Nonuniformity noise in IR imaging sensors, which is due to pixel-to-pixel variation in the detectors responses, can considerably degrade the quality of IR images since it results in a fixed-pattern-noise (FPN) that is superimposed on the true image. Further, what maes matter worse is that the nonuniformity slowly varies over time, and depending on the FPA technology, this drift can tae from minutes to hours. In order to solve this problem, several scene-based nonuniformity correction (NUC) techniques have been developed 1,, 3, 4. Scene-based techniques perform the NUC Department of Electrical Engineering, Universidad de La Frontera. Casilla 54-D, Temuco, Chile. csmarti@ufro.cl Department of Electrical Engineering, University of Concepción. Casilla 160-C, Concepción, Chile. sertorre@udec.cl using only the video sequences that are being imaged, not requiring any ind of laboratory calibration technique. Our group has been active in the development of novel scene-based algorithms for NUC based on statistical estimation theory. In 5, 6 we have developed a Gauss-Marov model to capture the slow variation in the FPN and have utilized the model to adaptively estimate the nonuniformity in blocs of infrared video sequences using a Kalman Filter. In such wor the gain and offset are assumed uncorrelated. The principal contribution of this wor is to study and model the correlation between the nonuniformity gain and offset of each detector in the focal plane array. The model, called generalized Gauss-Marov model, is based in an experimental measure of the correlation between gain and offset nonuniformity. To estimate the correlation, we employ the detector parameters estimated by three well nown nonuniformity correction methods: a neural networ approach, a non-linear filtering method and a statistical algorithm. With those estimates the correlation function is obtained and the Gauss-Marov model is then generalized. To validate the proposed model experimental testing with several real infrared video sequences are performed. This paper is organized as follows. In Sections the IR-FPA read out model and the Gauss-Marov model are presented. In Section 3 the generalized Gauss-Marov model is developed. In Section 4 the correlation between gain and offset is measure and validate with sequence of real infrared data. In Section 5 the conclusions of the paper are summarized.. Gauss-Marov Model for The Gain and Offset Nonuniformity In this paper, we model the pixel-to-pixel variation in the detectors responses (nonuniformity) using the commonly used linear model for each pixel on the IR FPA. For the () th detector, the measured readout signal Y (n) at a given time n is represented by the approx-
2 imate linear relation Y (n) = X (1) (n)t (n) + X () (n) +V (n), where X (1) (n) and X () (n) are the gain and the offset of the th detector respectively and T (n) is the real incident infrared photon flux collected by the detector. The term V (n) represents the read-out noise associate with the th detector at time n. the subscript is omitted with the understanding that all operations are performed on a pixel by pixel basis. In practice, vectors of observations are captured over which the gain and offset are approximately constants. This is an acceptable assumption in applications where the observation vector is short. At a later time, other vectors of observation are captured for which these parameters may have drift. This slow random variation of the gain and offset between consecutive vectors of observation is modeled by a Gauss-Marov process. It is also considered that the gain and offset at th vector-time are uncorrelated. The driver noise for the gain and offset are assumed gaussian individually and mutually uncorrelated among them. Mathematically, the Gauss-Marov state equation is given by X +1 = Φ X + G W, (1) X (1) where X +1 = X () is the state vector at th vector time. Φ = is called the state transition α 0 0 β matrix between the states at th and + 1th time vector. α and β represent the level of drift in the gain and offset between consecutive vectors of observation, respectively. The parameter α (β ) is chosen according to the magnitude of the drift between the gain (offset) at state and the state + 1. G = is a identity matrix that randomly relates the driving noise vector W to the state vector X. W = W (1) W (), W (1) and W () are the driving noise for the gain and offset process, respectively. Finally, to complete the Gauss-Marov model, the observation model for the detector response is given by Y = H X + V, () where Y is the observation vector at th bloc, and the length of observation vector is l. H = T (1) T (l ) 1 is the irradiance at th bloc. V is the vector of independent, additive temporal noise elements in the th bloc. The main assumption in the observation model is that the input irradiance T (n) is a uniformly distributed random variable in the range T min,t max constituting the range common to all detectors of possible irradiance levels prior to saturation. For a complete development for the gain and offset estimation of each detector in the FPA using an optimized recursion algorithm such that Kalman filter, see Generalized Gauss-Marov Model In this paper, we the model the correlation between the nonuniformity gain and offset as follow: X +1 = Φ X + G W, (3) where the transition matrix between the states at th α ρ and + 1th time vector is modified to Φ = where ρ and δ represent the cross correlation between the gain and offset in consecutive vectors of observations, respectively. The stability of the system (3) is given by λ i (Φ ) < 1, i.e., the product of the eigenvalues must be lower than 1. Note that if we multiplied (3) by X +1 l and applying the expectation value we obtain E X +1 X +1 l T = Φ E X X +1 l T +G E W W l T G T, (4) where l = 0,1,... For a complete nowledge of the generalized Gauss-Marov model we need to determine the G and Φ matrices. In the next sections we present the development of the (4) for two cases: when l = 0 and when l > 0. The first case allow to find the matrix G and the second case, permit to figure out an expression for the transition matrix Φ Solution for G : Case where l = 0 In this case, setting l = 0, using (3) for X +1, and replacing on (4) we obtain E X +1 X T +1 = Φ E X X T Φ T + G E and solve (5) for E W W T we obtain that δ β W W T G T, (5)
3 E W W T { } = G 1 E X +1 X +1 T Φ E X X T Φ T and assuming that X is a stationary random variable and E X +1 X +1 T = E X X T E X +1 X +1 T Φ E X X T σ = a σ ab σ ab σb Φ T = σ 1 σ 0 σ 0 σ where ( σ0 = (1 αβ δρ)σ ab αβσa ρβσb, σ 1 = 1 α ) σa αρσ a b ρ σb and σ = ( 1 β ) σb 1 βδσ a b δ σa g1, and considering G = g 1 1 we found that E W W T σw = a 0 0 σ wb i.e., W is uncorrelated only if g 1 and g 1 satisfies that g 1 g 1σ 1 σ 0 g 1 σ 0 σ with σ w a = σ 1 g 1σ 0 + g 1σ = 0, and σ w b = σ g 1 σ0 + g 1σ1. For example, if g 1 = 0 then g 1 = σ 0 σ and G = 1 σ 0 σ 0 1, (7) which implies that the W is uncorrelated. Similarly, the expectation of W can be obtained from (3) by means E W = G 1 I Φ E X = G 1 I Φ E X Solution for Φ : Case where l > 0 In this case, replacing l > 0 on (4) we obtain E X +1 X T +1 l = Φ E X X +1 l T, (8) and expanding (8) on each elements of X, using correlation function notation and solve for Φ we obtain Φ = R X (l) R 1 (l 1), (9) X G 1 T, (6) where R X (l) = RX (1) X (1) (l) R X () X (1) (l) R X (1) X () (l) R X () X () (l). (10) Note that in 5 they assume that gain and offset are uncorrelated, α 0 = α 1 =... = α = α and β 0 = β 1 =... = β = β. Replacing this assumption on (9) results Φ = R X (1) X (1) (l) R X (1) X (1) (l 1) 0 0 R X () X () (l) R X () X () (l 1) = α 0, 0 β (11) this mean that the gain and offset are signals with exponential correlation function, i.e., following a Gauss- Marov model given in (1). Note that (11) correspond a particular case of (9), called the generalized transition matrix. In the next section we present the methodology to obtain an estimation of Φ using (9), i.e., based on estimate the correlation functions R X (l) and R X (l 1) obtained using some NUC published methods. 4. Estimation and Validation of the Generalized Model In this section we obtain an estimation of Φ from (9) using three NUC published methods: the neural networs approach 1,, 3, 4, nonlinear filtering 7, 8 and statistical algorithm 9. Briefly, the neural networs approach use a hidden layer lie a two-point NUC neurons that models the gain and offset of each detectors. This NUC technique is derived using and adaptive least-square (LMS) approach, maing it a comparison between neighboring pixels as new frames of data are input (the retina-lie processing). The second NUC method, nonlinear filtering, is used to obtain an individual estimation of gain and offset with recursive least-square (RLS) based on minimizing the error between a preliminary true scene estimation and the corresponding frame corrected for each pixel. This preliminary estimation of true irradiance is performed by a spatial nonuniformity filtering of readout data reducing the spatial signal resolution that is restored by the RLS algorithm. Both NUC methods, LMS and RLS based, perform the detector parameters estimation in a frame by frame basis. On the other hand, the statistical algorithm assumes that the irradiance at each detector is an uniformly distributed random variable and that the gain and offset are temporally constant within certain blocs of time. This mean that the gain and offset estimation is performed in a bloc of frame basis, i.e., when a new bloc of raw data arrives, the method perform the NU correction by a finite-impulse response filter and then,
4 the temporal filter is adapted to the changes in the detector parameters and the gain and offset are updated to the next arrivals of raw data. Now, to obtain an estimation for the gain and offset (X (1),X () ) we use real infrared data captured by the InSb FPA cooled camera (Amber Model AE-418) operating in 3 5µm range. The sequences were collected at 1 PM with 3000 frames collected at a rate of 30 frames/s and each pixel was quantized in 16 bit integers. Using this raw data, we are able to estimate for each method the detector parameters correlation function R X (1) and R X (0) to obtain a suitable estimation of Φ. Additionally, in this wor we assume that Φ 0 = Φ 1 =... = Φ = Φ. Now, the correlation function between X i and X j with i = 1, and j = 1,, can be estimated by 10 ˆR n X (i) X ( j) (l ) = 1 n n 1 =0 X (i) ()X ( j) ( + l ) (1) where ˆR n X (i) X ( j)(l ) is the estimation for R X (i) X ( j) (1) and R X (i) X ( j) (0) when i = 1,, j = 1, and l = 0,1. In (1) the nowledge of X (1) and X () at each time is required. Since from the three NUC methods we have an estimation of the gain and offset we can express (1) by the following recursive algorithm ˆR n X (i) X ( j) (l ) = n 1 n ˆR n 1 X (i) X ( j)(l )+ 1 n X (i) (n 1)X ( j) (n 1+l ), (13) and then, using (9) we finally obtain the estimated Φ for each method. In summary, to obtain an estimation for Φ we needed to mae the following steps: i) using a real infrared sequence, estimate the gain and offset by three NUC methods, ii) estimate the correlation function for R X (1) and R X (0) using (13) and finally, iii) obtain the estimated Φ for each method using (11). To validate the parameters of each model, we considered two aspect: the NUC performance and the estimation error. For study the NUC capability of each model we use the roughness parameters defined by ρ( f ) = h f 1 + h T f 1 f 1, (14) where h is a horizontal mas, f 1 is the L 1 norm of frame f and represents discrete convolution. In the other hand, to validate the model we use equation () in order to calculate the mean-square error (MSE) between the raw data Y and the estimation Y ˆ ( Φ j ) for each model, given by 11 Table 1. Validation Results for different models Model λ i (φ j ) ρ( f ) MSE Φ Φ Φ MSE = 1 N N 1 Y ˆ T Y ( Φ ) j Y ˆ Y ( Φ j ), (15) =0 where Y ˆ ( Φ j ) = H X = H Φ X 1 is the a priori estimate of Y ˆ based on the model Φ j for j = 1,,3. Following the previous procedure we obtain three generalized models corresponding to each NUC algorithm and they are given by: Φ 1 = Φ = Φ 3 = (16) (17) (18) where Φ 1, Φ and Φ 3 correspond to the model obtained with the neural networ approach, nonlinear filtering and statistical algorithm respectively. Note that all NUC methods presented are stable (see Table 1) and shows that X (1) and X () are correlated. This mean that the assumption that the gain and offset in consecutive vector of observation are uncorrelated is not valid. Even more, whereas LMS and RLS algorithm exhibit a large correlation, the statistical algorithm shows that the gain strongly is influenced by the gain and offset of the previous state, i.e., we can assume that the offset indirectly is estimate from the gain. This is a very important conclusion because it means that we need the nowledge of only one state to reduce NU. In Table 1, the λ i (Φ j ) for each model, the spatial average of ρ and MSE parameters are presented. We can observe that the best performance is obtained for the statistical algorithm ( Φ 3 ). Additionally, from equations (16),(17) and (18) we can observe that as the magnitude of the drift between the gain (offset) at state and state + 1 is increasingly different (i.e., the correlation between gain and offset is considered), the NU correction and the model estimation are improved. This in addition confirms the assumption of the gain and offset
5 (a) (b) (c) (d) Figure 1. The 0 th frame from the 1 PM sequence captured with the Amber infrared camera: a) the corrupted raw frame, the corrected frame with b) neural networ approach, c) nonlinear filtering, and d) statistical algorithm. It can be seen by using only the naed eyes that a good NUC is obtained from the three NUC methods. are temporally constant in a bloc of time. Finally, for illustration proposed, the LMS-based, RLS-based and the statistical algorithm were applied to the 0-frame image sequence, and the results is shown in Figure 1. In this figure, Fig. 1(a) correspond to the corrupted frames, Fig. 1(b), Fig. 1(c) and Fig. 1(d) shows the corrected frames with neural networ approach, nonlinear filtering, and statistical algorithm respectively. This figure clearly shows that the statistical method present a good reduction of FPN. 5. Conclusions The correlation between gain and offset nonuniformity in an infrared focal plane array has been analyzed. It was shown experimentally using real IR data that the gain and offset for each pixel in an infrared focal plane array are correlated. Indeed, when the asymmetry is more notable in the generalized transition matrix, i.e., the magnitude of the drift between the gain (offset) at consecutive bloc of frames is increasingly different, the method has shown an improved reduction in both nonuniformity and in the error estimation. 6. ACKNOWLEDGMENTS This wor was partially supported by Grants Milenio ICM P0-049, by FONDECYT and Universidad de La Frontera, Proyecto DIUFRO The authors wish to than Ernest E. Armstrong (OptiMetrics Inc., USA) for collecting the data, and the United States Air Force Research Laboratory, Ohio, USA. References 1 D. Scribner, K. Sarady, M. Kruer: Adaptive Nonuniformity Correction for Infrared Focal Plane Arrays using Neural Networs. Proceeding of SPIE (1991) D. Scribner, K. Sarady, M. Kruer: Adaptive Retina-lie Preprocessing for Imaging Detector Arrays. Proceeding of the IEEE International Conference on Neural Networs. 3. (1993)
6 3 S. Torres, E. Vera, R. Reeves, and S. Sobarzo: Adaptive Scene-Based Nonuniformity Correction Method for Infrared Focal Plane Arrays. Proceeding of SPIE (003) E. Vera, S. Torres: Fast Adaptive Nonuniformity Correction for Infrared Focal Plane Arrays. Journal on Applied Signal Processing EURASIP. (005). 5 S. Torres, M. Hayat: Kalman Filtering for Adaptive Nonuniformity Correction in Infrared Focal Plane Arrays. The JOSA-A Opt. Soc. of America. 0 (003) S. Torres, J. Pezoa, and M. Hayat: Scene-based Nonuniformity Correction for Focal Plane Arrays Using the Method of the Inverse Covariance Form. OSA App. Opt. Inf. Proc. 4 (003) F. Torres, S. Torres, and C. San Martn: A Recursive Least Square Adaptive Filter for Nonuniformity Correction of Infrared Image Sequences. LNCS (005) F. Torres, C. San Martin, and S. Torres: A RLS Filter for Nonuniformity and Ghosting Correction of Infrared Image Sequences. Sequences. LNCS. 45 (006) M. Hayat, S. Torres, E. Amstrong, S. Cain, and B. Yasuda: Statistical Algorithm fo Nonuniformity Correction in Focal Plane Arrays. Applied Optics. 38 (1999) Charles W. Therrien : Discrete Random Signal and Statistical Signal Processing, Prentice Hall Press, USA, J. Pezoa, M. Hayat, S. Torres, and Md. Rahman: Multimodel Kalman Filtering for Adaptive Nonuniformity Correction in Infrared Sensors. JOSA A. 3 (006)
Multimodel Kalman filtering for adaptive nonuniformity correction in infrared sensors
1282 J. Opt. Soc. Am. A/ Vol. 23, No. 6/ June 2006 Pezoa et al. Multimodel Kalman filtering for adaptive nonuniformity correction in infrared sensors Jorge E. Pezoa and Majeed M. Hayat Department of Electrical
More informationScene-based nonuniformity correction for focal plane arrays by the method of the inverse covariance form
Scene-based nonuniformity correction for focal plane arrays by the method of the inverse covariance form Sergio N. Torres, Jorge E. Pezoa, and Majeed M. Hayat What is to our knowledge a new scene-based
More informationKalman filtering for adaptive nonuniformity correction in infrared focal-plane arrays
470 J. Opt. Soc. Am. A/ Vol. 0, No. 3/ March 003 S. N. Torres and M. M. Hayat Kalman filtering for adaptive nonuniformity correction in infrared focal-plane arrays Sergio N. Torres Department of Electrical
More informationKalman filtering for drift correction in IR detectors
Proceedings of the 5th WSEAS Int. Conf. on Signal Processing, Robotics and Automation, Madrid, Spain, February 15-17, 006 (pp58-63 Kalman filtering for drift correction in IR detectors PIÑEIRO AE, J. (1,
More informationImproved Kalman Filter Initialisation using Neurofuzzy Estimation
Improved Kalman Filter Initialisation using Neurofuzzy Estimation J. M. Roberts, D. J. Mills, D. Charnley and C. J. Harris Introduction It is traditional to initialise Kalman filters and extended Kalman
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 information12.4 Known Channel (Water-Filling Solution)
ECEn 665: Antennas and Propagation for Wireless Communications 54 2.4 Known Channel (Water-Filling Solution) The channel scenarios we have looed at above represent special cases for which the capacity
More informationLinear Prediction Theory
Linear Prediction Theory Joseph A. O Sullivan ESE 524 Spring 29 March 3, 29 Overview The problem of estimating a value of a random process given other values of the random process is pervasive. Many problems
More informationHadamard speckle contrast reduction
Published in OPTICS LETTERS 29, -3 (2004) Hadamard speckle contrast reduction Jahja I. Trisnadi Silicon Light Machines, 385 Moffett Park Drive, Sunnyvale, California 94089 The conditon for a diffuser to
More informationIntroduction p. 1 Fundamental Problems p. 2 Core of Fundamental Theory and General Mathematical Ideas p. 3 Classical Statistical Decision p.
Preface p. xiii Acknowledgment p. xix Introduction p. 1 Fundamental Problems p. 2 Core of Fundamental Theory and General Mathematical Ideas p. 3 Classical Statistical Decision p. 4 Bayes Decision p. 5
More informationADAPTIVE INVERSE CONTROL BASED ON NONLINEAR ADAPTIVE FILTERING. Information Systems Lab., EE Dep., Stanford University
ADAPTIVE INVERSE CONTROL BASED ON NONLINEAR ADAPTIVE FILTERING Bernard Widrow 1, Gregory Plett, Edson Ferreira 3 and Marcelo Lamego 4 Information Systems Lab., EE Dep., Stanford University Abstract: Many
More informationDesign of FIR Smoother Using Covariance Information for Estimating Signal at Start Time in Linear Continuous Systems
Systems Science and Applied Mathematics Vol. 1 No. 3 2016 pp. 29-37 http://www.aiscience.org/journal/ssam Design of FIR Smoother Using Covariance Information for Estimating Signal at Start Time in Linear
More informationDiffusion LMS Algorithms for Sensor Networks over Non-ideal Inter-sensor Wireless Channels
Diffusion LMS Algorithms for Sensor Networs over Non-ideal Inter-sensor Wireless Channels Reza Abdolee and Benoit Champagne Electrical and Computer Engineering McGill University 3480 University Street
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 informationNeural networks: Associative memory
Neural networs: Associative memory Prof. Sven Lončarić sven.loncaric@fer.hr http://www.fer.hr/ipg 1 Overview of topics l Introduction l Associative memories l Correlation matrix as an associative memory
More informationAperture-Averaging - Theory and Measurements
Aperture-Averaging - Theory and Measurements N. Perlot, D. Fritzsche German Aerospace Center (DLR), D-834 Wessling, Germany ABSTRACT Atmospheric laser communications using direct-detection systems do suffer
More informationHere represents the impulse (or delta) function. is an diagonal matrix of intensities, and is an diagonal matrix of intensities.
19 KALMAN FILTER 19.1 Introduction In the previous section, we derived the linear quadratic regulator as an optimal solution for the fullstate feedback control problem. The inherent assumption was that
More informationMaximum Likelihood Diffusive Source Localization Based on Binary Observations
Maximum Lielihood Diffusive Source Localization Based on Binary Observations Yoav Levinboo and an F. Wong Wireless Information Networing Group, University of Florida Gainesville, Florida 32611-6130, USA
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 informationVEHICLE WHEEL-GROUND CONTACT ANGLE ESTIMATION: WITH APPLICATION TO MOBILE ROBOT TRACTION CONTROL
1/10 IAGNEMMA AND DUBOWSKY VEHICLE WHEEL-GROUND CONTACT ANGLE ESTIMATION: WITH APPLICATION TO MOBILE ROBOT TRACTION CONTROL K. IAGNEMMA S. DUBOWSKY Massachusetts Institute of Technology, Cambridge, MA
More informationRecursive Least Squares for an Entropy Regularized MSE Cost Function
Recursive Least Squares for an Entropy Regularized MSE Cost Function Deniz Erdogmus, Yadunandana N. Rao, Jose C. Principe Oscar Fontenla-Romero, Amparo Alonso-Betanzos Electrical Eng. Dept., University
More informationWhat is Image Deblurring?
What is Image Deblurring? When we use a camera, we want the recorded image to be a faithful representation of the scene that we see but every image is more or less blurry, depending on the circumstances.
More informationA Novel Activity Detection Method
A Novel Activity Detection Method Gismy George P.G. Student, Department of ECE, Ilahia College of,muvattupuzha, Kerala, India ABSTRACT: This paper presents an approach for activity state recognition of
More informationDigital Phase-Locked Loop and its Realization
Proceedings of the 9th WSEAS International Conference on APPLIE INFORMATICS AN COMMUNICATIONS (AIC '9) igital Phase-Loced Loop and its Realiation Tsai-Sheng Kao 1, Sheng-Chih Chen 2, Yuan-Chang Chang 1,
More informationNoise, Image Reconstruction with Noise!
Noise, Image Reconstruction with Noise! EE367/CS448I: Computational Imaging and Display! stanford.edu/class/ee367! Lecture 10! Gordon Wetzstein! Stanford University! What s a Pixel?! photon to electron
More informationLecture Notes 2 Charge-Coupled Devices (CCDs) Part I. Basic CCD Operation CCD Image Sensor Architectures Static and Dynamic Analysis
Lecture Notes 2 Charge-Coupled Devices (CCDs) Part I Basic CCD Operation CCD Image Sensor Architectures Static and Dynamic Analysis Charge Well Capacity Buried channel CCD Transfer Efficiency Readout Speed
More informationExtended Kalman Filtering and Phase Detector Modeling for a Digital Phase-Locked Loop
Extended Kalman Filtering and Phase etector Modeling for a igital Phase-Loced Loop Tsai-Sheng Kao 1, Sheng-Chih Chen 2, Yuan-Chang Chang 1, Sheng-Yun Hou 1, and Chang-Jung Juan 1 1 epartment of Electronic
More informationBobby Hunt, Mariappan S. Nadar, Paul Keller, Eric VonColln, and Anupam Goyal III. ASSOCIATIVE RECALL BY A POLYNOMIAL MAPPING
Synthesis of a Nonrecurrent Associative Memory Model Based on a Nonlinear Transformation in the Spectral Domain p. 1 Bobby Hunt, Mariappan S. Nadar, Paul Keller, Eric VonColln, Anupam Goyal Abstract -
More information18/10/2017. Image Enhancement in the Spatial Domain: Gray-level transforms. Image Enhancement in the Spatial Domain: Gray-level transforms
Gray-level transforms Gray-level transforms Generic, possibly nonlinear, pointwise operator (intensity mapping, gray-level transformation): Basic gray-level transformations: Negative: s L 1 r Generic log:
More informationPattern Recognition with Slow Feature Analysis
Cognitive Sciences EPrint Archive (CogPrints) 404, (2005) http://cogprints.org/404/ Pattern Recognition with Slow Feature Analysis Pietro Beres Institute for Theoretical Biology Humboldt University Berlin
More informationImprovement of Himawari-8 observation data quality
Improvement of Himawari-8 observation data quality 3 July 2017 Meteorological Satellite Center Japan Meteorological Agency The Japan Meteorological Agency (JMA) plans to modify its Himawari-8 ground processing
More informationSliding Window Recursive Quadratic Optimization with Variable Regularization
11 American Control Conference on O'Farrell Street, San Francisco, CA, USA June 29 - July 1, 11 Sliding Window Recursive Quadratic Optimization with Variable Regularization Jesse B. Hoagg, Asad A. Ali,
More informationExpressions for the covariance matrix of covariance data
Expressions for the covariance matrix of covariance data Torsten Söderström Division of Systems and Control, Department of Information Technology, Uppsala University, P O Box 337, SE-7505 Uppsala, Sweden
More information9.4 Enhancing the SNR of Digitized Signals
9.4 Enhancing the SNR of Digitized Signals stepping and averaging compared to ensemble averaging creating and using Fourier transform digital filters removal of Johnson noise and signal distortion using
More informationCity, University of London Institutional Repository
City Research Online City, University of London Institutional Repository Citation: Zhao, S., Shmaliy, Y. S., Khan, S. & Liu, F. (2015. Improving state estimates over finite data using optimal FIR filtering
More informationPERIODIC KALMAN FILTER: STEADY STATE FROM THE BEGINNING
Journal of Mathematical Sciences: Advances and Applications Volume 1, Number 3, 2008, Pages 505-520 PERIODIC KALMAN FILER: SEADY SAE FROM HE BEGINNING MARIA ADAM 1 and NICHOLAS ASSIMAKIS 2 1 Department
More informationMulti-Objective Autonomous Spacecraft Motion Planning around Near-Earth Asteroids using Machine Learning
Multi-Objective Autonomous Spacecraft Motion Planning around Near-Earth Asteroids using Machine Learning CS 229 Final Project - Fall 2018 Category: Physical Sciences Tommaso Guffanti SUNetID: tommaso@stanford.edu
More informationLogarithmic quantisation of wavelet coefficients for improved texture classification performance
Logarithmic quantisation of wavelet coefficients for improved texture classification performance Author Busch, Andrew, W. Boles, Wageeh, Sridharan, Sridha Published 2004 Conference Title 2004 IEEE International
More informationPupil matching of Zernike aberrations
Pupil matching of Zernike aberrations C. E. Leroux, A. Tzschachmann, and J. C. Dainty Applied Optics Group, School of Physics, National University of Ireland, Galway charleleroux@yahoo.fr Abstract: The
More informationFor final project discussion every afternoon Mark and I will be available
Worshop report 1. Daniels report is on website 2. Don t expect to write it based on listening to one project (we had 6 only 2 was sufficient quality) 3. I suggest writing it on one presentation. 4. Include
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 informationQuality & Information Content Of CHRIS Hyper-Spectral Data
Quality & Information Content Of CHRIS Hyper-Spectral Data B. Aiazzi, S. Baronti, P. Marcoionni, I. Pippi, M. Selva Institute of Applied Physics Nello Carrara IFAC-CNR, Florence Italy Information-Theoretic
More informationGround and On-Orbit Characterization and Calibration of the Geosynchronous Imaging Fourier Transform Spectrometer (GIFTS)
Ground and On-Orbit Characterization and Calibration of the Geosynchronous Imaging Fourier Transform Spectrometer (GIFTS) John D. Elwell 1, Deron K. Scott 1 Henry E. Revercomb 2, Fred A. Best 2, Robert
More informationNumerical method of active thermography for the reconstruction of back wall geometry
11 th International Conference on Quantitative InfraRed hermography Numerical method of active thermography for the reconstruction of bac wall geometry by Regina Richter *, Christiane Maierhofer *, Marc
More informationImproved least-squares-based combiners for diffusion networks
Improved least-squares-based combiners for diffusion networs Jesus Fernandez-Bes, Luis A. Azpicueta-Ruiz, Magno T. M. Silva, and Jerónimo Arenas-García Universidad Carlos III de Madrid Universidade de
More informationConvolutional networks. Sebastian Seung
Convolutional networks Sebastian Seung Convolutional network Neural network with spatial organization every neuron has a location usually on a grid Translation invariance synaptic strength depends on locations
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 informationEIGENFILTERS FOR SIGNAL CANCELLATION. Sunil Bharitkar and Chris Kyriakakis
EIGENFILTERS FOR SIGNAL CANCELLATION Sunil Bharitkar and Chris Kyriakakis Immersive Audio Laboratory University of Southern California Los Angeles. CA 9. USA Phone:+1-13-7- Fax:+1-13-7-51, Email:ckyriak@imsc.edu.edu,bharitka@sipi.usc.edu
More informationNew Recursive-Least-Squares Algorithms for Nonlinear Active Control of Sound and Vibration Using Neural Networks
IEEE TRANSACTIONS ON NEURAL NETWORKS, VOL. 12, NO. 1, JANUARY 2001 135 New Recursive-Least-Squares Algorithms for Nonlinear Active Control of Sound and Vibration Using Neural Networks Martin Bouchard,
More informationExtraction of incident irradiance from LWIR hyperspectral imagery
DRDC-RDDC-215-P14 Extraction of incident irradiance from LWIR hyperspectral imagery Pierre Lahaie, DRDC Valcartier 2459 De la Bravoure Road, Quebec, Qc, Canada ABSTRACT The atmospheric correction of thermal
More informationUsing Entropy and 2-D Correlation Coefficient as Measuring Indices for Impulsive Noise Reduction Techniques
Using Entropy and 2-D Correlation Coefficient as Measuring Indices for Impulsive Noise Reduction Techniques Zayed M. Ramadan Department of Electronics and Communications Engineering, Faculty of Engineering,
More informationVISION TRACKING PREDICTION
VISION RACKING PREDICION Eri Cuevas,2, Daniel Zaldivar,2, and Raul Rojas Institut für Informati, Freie Universität Berlin, ausstr 9, D-495 Berlin, German el 0049-30-83852485 2 División de Electrónica Computación,
More informationFATIGUE DAMAGE ASSESMENT USING THERMOELASTIC STRESS ANALYSIS
FATIGUE DAMAGE ASSESMENT USING THERMOELASTIC STRESS ANALYSIS F.A. Díaz, E.A. Patterson and J.R. Yates The University of Sheffield, Department of Mechanical Engineering, Mappin street, Sheffield S10 3JD,
More informationEE-597 Notes Quantization
EE-597 Notes Quantization Phil Schniter June, 4 Quantization Given a continuous-time and continuous-amplitude signal (t, processing and storage by modern digital hardware requires discretization in both
More informationOn the FPA infrared camera transfer function calculation
On the FPA infrared camera transfer function calculation (1) CERTES, Université Paris XII Val de Marne, Créteil, France (2) LTM, Université de Bourgogne, Le Creusot, France by S. Datcu 1, L. Ibos 1,Y.
More informationCLOUD NOWCASTING: MOTION ANALYSIS OF ALL-SKY IMAGES USING VELOCITY FIELDS
CLOUD NOWCASTING: MOTION ANALYSIS OF ALL-SKY IMAGES USING VELOCITY FIELDS Yézer González 1, César López 1, Emilio Cuevas 2 1 Sieltec Canarias S.L. (Santa Cruz de Tenerife, Canary Islands, Spain. Tel. +34-922356013,
More informationShort-Term Electricity Price Forecasting
1 Short-erm Electricity Price Forecasting A. Arabali, E. Chalo, Student Members IEEE, M. Etezadi-Amoli, LSM, IEEE, M. S. Fadali, SM IEEE Abstract Price forecasting has become an important tool in the planning
More information5.1 2D example 59 Figure 5.1: Parabolic velocity field in a straight two-dimensional pipe. Figure 5.2: Concentration on the input boundary of the pipe. The vertical axis corresponds to r 2 -coordinate,
More informationNear-infrared guiding and tip-tilt correction for the UC Berkeley Infrared Spatial Interferometer
Near-infrared guiding and tip-tilt correction for the UC Berkeley Infrared Spatial Interferometer E. A. Lipman, M. Bester, W. C. Danchi, and C. H. Townes Space Sciences Laboratory and Department of Physics
More informationEFFICIENT LINEAR MIMO ADAPTIVE ŜP INVERSE CONTROL. Gregory L. Plett E
PSfrag replacements C COPY P COPY P REV X COPY M M COPY S EFFICIENT LINEAR MIMO ADAPTIVE ŜP INVERSE CONTROL ŜP COPY Ŝ Gregory L. Plett E 1 E COPY F F COPY University of Colorado at Colorado zi Springs,
More informationPerformance Analysis of an Adaptive Algorithm for DOA Estimation
Performance Analysis of an Adaptive Algorithm for DOA Estimation Assimakis K. Leros and Vassilios C. Moussas Abstract This paper presents an adaptive approach to the problem of estimating the direction
More informationUse of Matlab package for external local calibration of IR camera with microbolometric FPA detector
Use of Matlab package for external local calibration of IR camera with microlometric FPA detector More info about this article: http://www.ndt.net/?id=20646 Abstract by T. Kruczek*, W. Adamczyk* and G.
More informationSPATIO-TEMPORAL PREDICTION FOR ADAPTIVE OPTICS WAVEFRONT RECONSTRUCTORS
SPATIO-TEMPORAL PREDICTION FOR ADAPTIVE OPTICS WAVEFRONT RECONSTRUCTORS Michael Lloyd-Hart and Patrick McGuire Center for Astronomical Adaptive Optics, Steward Observatory, University of Arizona, Tucson,
More informationSTOCHASTIC STABILITY OF EXTENDED FILTERING FOR NONLINEAR SYSTEMS WITH MEASUREMENT PACKET LOSSES
Proceedings of the IASTED International Conference Modelling, Identification and Control (AsiaMIC 013) April 10-1, 013 Phuet, Thailand STOCHASTIC STABILITY OF EXTENDED FILTERING FOR NONLINEAR SYSTEMS WITH
More informationPAPER 338 OPTICAL AND INFRARED ASTRONOMICAL TELESCOPES AND INSTRUMENTS
MATHEMATICAL TRIPOS Part III Monday, 12 June, 2017 1:30 pm to 3:30 pm PAPER 338 OPTICAL AND INFRARED ASTRONOMICAL TELESCOPES AND INSTRUMENTS Attempt no more than TWO questions. There are THREE questions
More informationRecursive On-orbit Calibration of Star Sensors
Recursive On-orbit Calibration of Star Sensors D. odd Griffith 1 John L. Junins 2 Abstract Estimation of calibration parameters for a star tracer is investigated. Conventional estimation schemes are evaluated
More informationHand Written Digit Recognition using Kalman Filter
International Journal of Electronics and Communication Engineering. ISSN 0974-2166 Volume 5, Number 4 (2012), pp. 425-434 International Research Publication House http://www.irphouse.com Hand Written Digit
More informationAn Adaptive Sensor Array Using an Affine Combination of Two Filters
An Adaptive Sensor Array Using an Affine Combination of Two Filters Tõnu Trump Tallinn University of Technology Department of Radio and Telecommunication Engineering Ehitajate tee 5, 19086 Tallinn Estonia
More informationSparse Least Mean Square Algorithm for Estimation of Truncated Volterra Kernels
Sparse Least Mean Square Algorithm for Estimation of Truncated Volterra Kernels Bijit Kumar Das 1, Mrityunjoy Chakraborty 2 Department of Electronics and Electrical Communication Engineering Indian Institute
More informationAnalysis of Finite Wordlength Effects
Analysis of Finite Wordlength Effects Ideally, the system parameters along with the signal variables have infinite precision taing any value between and In practice, they can tae only discrete values within
More informationMusical noise reduction in time-frequency-binary-masking-based blind source separation systems
Musical noise reduction in time-frequency-binary-masing-based blind source separation systems, 3, a J. Čermá, 1 S. Arai, 1. Sawada and 1 S. Maino 1 Communication Science Laboratories, Corporation, Kyoto,
More informationMULTI-RESOLUTION SIGNAL DECOMPOSITION WITH TIME-DOMAIN SPECTROGRAM FACTORIZATION. Hirokazu Kameoka
MULTI-RESOLUTION SIGNAL DECOMPOSITION WITH TIME-DOMAIN SPECTROGRAM FACTORIZATION Hiroazu Kameoa The University of Toyo / Nippon Telegraph and Telephone Corporation ABSTRACT This paper proposes a novel
More informationMultimedia communications
Multimedia communications Comunicazione multimediale G. Menegaz gloria.menegaz@univr.it Prologue Context Context Scale Scale Scale Course overview Goal The course is about wavelets and multiresolution
More informationPolarization Processing at LWA Stations
Polarization Processing at LWA Stations Steve Ellingson October 28, 2007 Contents 1 Introduction 2 2 Theory 2 3 Frequency-Dependent Antenna Polarization 3 Bradley Dept. of Electrical & Computer Engineering,
More information1 Kalman Filter Introduction
1 Kalman Filter Introduction You should first read Chapter 1 of Stochastic models, estimation, and control: Volume 1 by Peter S. Maybec (available here). 1.1 Explanation of Equations (1-3) and (1-4) Equation
More informationA NOVEL OPTIMAL PROBABILITY DENSITY FUNCTION TRACKING FILTER DESIGN 1
A NOVEL OPTIMAL PROBABILITY DENSITY FUNCTION TRACKING FILTER DESIGN 1 Jinglin Zhou Hong Wang, Donghua Zhou Department of Automation, Tsinghua University, Beijing 100084, P. R. China Control Systems Centre,
More informationRequirements for the Star Tracker Parallel Science Programme
Requirements for the Star Tracker Parallel Science Programme Rømer System Definition Phase 2000/2001 Document no.: MONS/IFA/PL/RS/0003(1) Date: 22.04.2001 Prepared by: Hans Kjeldsen and Tim Bedding Checked
More informationP. Bruschi - Notes on Mixed Signal Design
Chapter 1. Concepts and definitions about Data Acquisition Systems Electronic systems An electronic systems is a complex electronic networ, which interacts with the physical world through sensors (input
More informationAuxiliaire d enseignement Nicolas Ayotte
2012-02-15 GEL 4203 / GEL 7041 OPTOÉLECTRONIQUE Auxiliaire d enseignement Nicolas Ayotte GEL 4203 / GEL 7041 Optoélectronique VI PN JUNCTION The density of charge sign Fixed charge density remaining 2
More informationMini-Course 07 Kalman Particle Filters. Henrique Massard da Fonseca Cesar Cunha Pacheco Wellington Bettencurte Julio Dutra
Mini-Course 07 Kalman Particle Filters Henrique Massard da Fonseca Cesar Cunha Pacheco Wellington Bettencurte Julio Dutra Agenda State Estimation Problems & Kalman Filter Henrique Massard Steady State
More informationSIMULTANEOUS STATE AND PARAMETER ESTIMATION USING KALMAN FILTERS
ECE5550: Applied Kalman Filtering 9 1 SIMULTANEOUS STATE AND PARAMETER ESTIMATION USING KALMAN FILTERS 9.1: Parameters versus states Until now, we have assumed that the state-space model of the system
More informationImploded Shell Parameter Estimation Based on Radiograph Analysis. George Liu. Pittsford Sutherland High School. LLE Advisor: Reuben Epstein
Imploded Shell Parameter Estimation Based on Radiograph Analysis George Liu Pittsford Sutherland High School LLE Advisor: Reuben Epstein Laboratory for Laser Energetics University of Rochester Summer High
More informationPOLYNOMIAL SINGULAR VALUES FOR NUMBER OF WIDEBAND SOURCES ESTIMATION AND PRINCIPAL COMPONENT ANALYSIS
POLYNOMIAL SINGULAR VALUES FOR NUMBER OF WIDEBAND SOURCES ESTIMATION AND PRINCIPAL COMPONENT ANALYSIS Russell H. Lambert RF and Advanced Mixed Signal Unit Broadcom Pasadena, CA USA russ@broadcom.com Marcel
More informationPhase-Correlation Motion Estimation Yi Liang
EE 392J Final Project Abstract Phase-Correlation Motion Estimation Yi Liang yiliang@stanford.edu Phase-correlation motion estimation is studied and implemented in this work, with its performance, efficiency
More informationARTEFACT DETECTION IN ASTROPHYSICAL IMAGE DATA USING INDEPENDENT COMPONENT ANALYSIS. Maria Funaro, Erkki Oja, and Harri Valpola
ARTEFACT DETECTION IN ASTROPHYSICAL IMAGE DATA USING INDEPENDENT COMPONENT ANALYSIS Maria Funaro, Erkki Oja, and Harri Valpola Neural Networks Research Centre, Helsinki University of Technology P.O.Box
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 informationEM-algorithm for Training of State-space Models with Application to Time Series Prediction
EM-algorithm for Training of State-space Models with Application to Time Series Prediction Elia Liitiäinen, Nima Reyhani and Amaury Lendasse Helsinki University of Technology - Neural Networks Research
More informationCompression methods: the 1 st generation
Compression methods: the 1 st generation 1998-2017 Josef Pelikán CGG MFF UK Praha pepca@cgg.mff.cuni.cz http://cgg.mff.cuni.cz/~pepca/ Still1g 2017 Josef Pelikán, http://cgg.mff.cuni.cz/~pepca 1 / 32 Basic
More informationEUSIPCO
EUSIPCO 3 569736677 FULLY ISTRIBUTE SIGNAL ETECTION: APPLICATION TO COGNITIVE RAIO Franc Iutzeler Philippe Ciblat Telecom ParisTech, 46 rue Barrault 753 Paris, France email: firstnamelastname@telecom-paristechfr
More informationJ. Liang School of Automation & Information Engineering Xi an University of Technology, China
Progress In Electromagnetics Research C, Vol. 18, 245 255, 211 A NOVEL DIAGONAL LOADING METHOD FOR ROBUST ADAPTIVE BEAMFORMING W. Wang and R. Wu Tianjin Key Lab for Advanced Signal Processing Civil Aviation
More informationAnalysis of incremental RLS adaptive networks with noisy links
Analysis of incremental RLS adaptive networs with noisy lins Azam Khalili, Mohammad Ali Tinati, and Amir Rastegarnia a) Faculty of Electrical and Computer Engineering, University of Tabriz Tabriz 51664,
More informationNovel spectrum sensing schemes for Cognitive Radio Networks
Novel spectrum sensing schemes for Cognitive Radio Networks Cantabria University Santander, May, 2015 Supélec, SCEE Rennes, France 1 The Advanced Signal Processing Group http://gtas.unican.es The Advanced
More informationSun Sensor Model. Center for Distributed Robotics Technical Report Number
Sun Sensor Model Nikolas Trawny and Stergios Roumeliotis Department of Computer Science & Engineering University of Minnesota Center for Distributed Robotics Technical Report Number -2005-00 January 2005
More informationLessons in Estimation Theory for Signal Processing, Communications, and Control
Lessons in Estimation Theory for Signal Processing, Communications, and Control Jerry M. Mendel Department of Electrical Engineering University of Southern California Los Angeles, California PRENTICE HALL
More informationA KALMAN FILTERING TUTORIAL FOR UNDERGRADUATE STUDENTS
A KALMAN FILTERING TUTORIAL FOR UNDERGRADUATE STUDENTS Matthew B. Rhudy 1, Roger A. Salguero 1 and Keaton Holappa 2 1 Division of Engineering, Pennsylvania State University, Reading, PA, 19610, USA 2 Bosch
More informationImage Contrast Enhancement and Quantitative Measuring of Information Flow
6th WSEAS International Conference on Information Security and Privacy, Tenerife, Spain, December 14-16, 007 17 Contrast Enhancement and Quantitative Measuring of Information Flow ZHENGMAO YE 1, HABIB
More informationSpitzer Space Telescope
Spitzer Space Telescope (A.K.A. The Space Infrared Telescope Facility) The Infrared Imaging Chain 1/38 The infrared imaging chain Generally similar to the optical imaging chain... 1) Source (different
More informationThe Kalman Filter. Data Assimilation & Inverse Problems from Weather Forecasting to Neuroscience. Sarah Dance
The Kalman Filter Data Assimilation & Inverse Problems from Weather Forecasting to Neuroscience Sarah Dance School of Mathematical and Physical Sciences, University of Reading s.l.dance@reading.ac.uk July
More informationColour. The visible spectrum of light corresponds wavelengths roughly from 400 to 700 nm.
Colour The visible spectrum of light corresponds wavelengths roughly from 4 to 7 nm. The image above is not colour calibrated. But it offers a rough idea of the correspondence between colours and wavelengths.
More informationA STATE-SPACE APPROACH FOR THE ANALYSIS OF WAVE AND DIFFUSION FIELDS
ICASSP 2015 A STATE-SPACE APPROACH FOR THE ANALYSIS OF WAVE AND DIFFUSION FIELDS Stefano Maranò Donat Fäh Hans-Andrea Loeliger ETH Zurich, Swiss Seismological Service, 8092 Zürich ETH Zurich, Dept. Information
More information