An Automatic Sequential Smoothing Method for Processing Biomechanical Kinematic Signals
|
|
- Peter Todd
- 6 years ago
- Views:
Transcription
1 An Automatic Sequential Smoothing Method or Processing Biomechanical Kinematic Signals F.J. ALONSO, F. ROMERO, D.R. SALGADO Department o Mechanical, Energetic and Material Engineering University o Extremadura Escuela de Ingenierías Industriales, Avda. de Elvas s/n, 66, Badajoz SPAIN jas@unex.es Abstract: - he combined application o Singular Spectrum Analysis (SSA) and Cluster Analysis to the automatic smoothing o raw kinematic signals is an alternative to the use o traditional digital iltering and spline based methods. SSA is a non parametric technique that decomposes original time series into a number o additive time series each o which can be easily identiied as being part o the noise present in the acquired signal. Nevertheless, the smoothing automation is not a trivial task. his work presents a heuristic automatic smoothing procedure or processing kinematic biomechanical signals based in sequential SSA. Cluster analysis is used to group the SSA decomposition in order to obtain several independent components in the requency domain. he procedure eliminates iteratively the noise present in the signal in a simple and intuitive way. he new method is applied to several signals to demonstrate its perormance. Key-Words: - Signal processing, Smoothing, Noise removal, Singular Spectrum Analysis, Signal dierentiation, Automatic smoothing, Cluster analysis. Introduction Motion capture systems used in biomechanical analysis introduce introduce systematic and random measurement errors that appear in the orm o highrequency noise in the recorded displacement signals. Raw displacement signals must be iltered or smoothed prior to dierentiation to avoid noise ampliication due to the ill-posed derivative estimation process [-3]. he iltering o displacement signals to obtain noiseless velocities and accelerations has been extensively treated in the literature. raditional iltering techniques include Digital Butterworth ilters, splines, and ilters based on spectral analysis [-5]. In order to ilter non-stationary signals, advanced iltering techniques like Discrete Wavelet ransorms [6], the Wigner Function [7-8], nonstationary Butterworth ilter [9] and Singular Spectrum Analysis (SSA) [] have been used. Nonetheless, the drawback in these cases is the complexity o devising an automatic and systematic procedure. A mother wavelet unction must be selected when using Discrete Wavelet ransorms, the iltering unction parameters must be chosen when using the Wigner Function, window length and grouping strategy must be selected when using SSA []. he goal o this paper is to demonstrate the advantages o automatic smoothing methods based on sequential Singular Spectrum Analysis (SSA) and cluster analysis techniques. SSA is a non-parametric technique that decomposes original time series into a number o additive time series each o which can be easily identiied as being part o the noise present in the acquired signal. his work presents a heuristic automatic smoothing procedure or processing kinematic biomechanical signals based on sequential SSA. Cluster analysis is used to group the SSA decomposition in order to obtain several independent components in the requency domain. he procedure then applies sequential SSA to eliminate iteratively the noise present in the signal in a simple and intuitive way. Singular Spectrum Analysis Singular spectrum analysis is a novel nonparametric technique o time series analysis based on principles o multivariate statistics. It decomposes a given time series into an additive set o independent time series. he set o series resulting rom the decomposition can be interpreted as consisting o a trend representing the signal mean at each instant, a set o periodic series, and an aperiodic noise []. he original application o SSA was to extract trends rom climatic and geophysical time ISBN:
2 series and to identiy periodic motion in complex dynamical systems. he SSA method builds a Hankel matrix, called the trajectory matrix, rom the original time series in a process called embedding. his matrix consists o vectors obtained by means o a sliding window that traverses the series. he trajectory matrix is then subjected to a singular value decomposition (SVD). he SVD decomposes the trajectory matrix into a sum o unit-rank matrices known as elementary matrices. Each o these matrices can be transormed into a reconstructed time series. Elementary matrices are no longer Hankel matrices, but an approximate time series may be recovered by taking the average o the diagonals (diagonal averaging). he resulting time series are called principal components []. he sum o all the principal components is equal to the original time series. he objective is to obtain a requency decomposition o the original signal in which the latent low-requency signal can be detected in a simple ashion. he SSA decomposition algorithm will be described in the ollowing. A more detailed explanation may be ound in Golyandina et al. []. he above description o SSA may be expressed in ormal terms as ollows: Step. Embedding Let F = (,,, ) be the length N time series N representing the original signal. Let L be the window length, with < L < N and L an integer. Each column X j o the Hankel matrix corresponds to the "snapshot" taken by the sliding window: X j = ( j, j,, j+ L ), j =,,, K where K = N L + is the number o columns, i.e., the number o dierent possible positions o the said window. he matrix X = (X, X,, XK ) is a Hankel matrix since all elements on the diagonal i + j = c are equal. his matrix is sometimes reerred to as the trajectory matrix. he orm o this matrix is: = L L L 3 L N L N N N L+ X () Step. Singular value decomposition (SVD) o the trajectory matrix. It can be proven that the trajectory matrix (or any matrix, or that matter) may be expressed as the summation o d rank-one elementary matrices X = E + E + + E d, where d is the number o non-zero eigenvalues o the L L matrix S = X X. he elementary matrices are given by Ei = λ i Ui Vi, i =,,, d, where λ, λ,, λd are the non-zero eigenvalues o S = X X in decreasing order, U, U,, Ud are the corresponding eigenvectors, and the vectors V i are obtained rom Vi = X Ui / λi. he norm o elementary matrix E i equals λ i, so that, the contribution o the irst matrices to the norm o X is much higher than the contribution o the last matrices. hereore, it is likely that these last matrices represent noise in the signal. he plot o the eigenvalues in decreasing order is called the singular spectrum, and gives the method its name. Reconstruction (diagonal averaging) At this step, each elementary matrix E i is transormed into a principal component o length N by applying a linear transormation known as diagonal averaging or Hankelization. he elementary matrices are not themselves Hankel matrices, so that to reconstruct each principal component one calculates the average along the diagonals i + j = c. he diagonal averaging algorithm [] is as ollows: Let Y be any o the elementary matrices E i o dimension L K, the elements o which are y ij, i L, j K. he time series G = g, g,, gn (the principal component) corresponding to this elementary matrix is given by: k + < + ym, k m+ or k L k m= L g k = y, + < m k m or L k K () L m= N K + < ym, k m+ or K k N N k m= k K + Where L = min( L, K), K = max( L, K), and the length N = L + K. It can be shown that the squared norm o each elementary matrix equals the ISBN:
3 corresponding eigenvalue, and that the squared norm o the trajectory matrix is the sum o the squared norms o the elementary matrices []. hus the ratio: d λ / λ (3) i i represents the contribution o the elementary matrix E i in the expansion o the trajectory matrix X = E + E + + E d. he largest eigenvalues in the singular spectrum represent the large amplitude components in the decomposition. Contrariwise, the lowamplitude components o the signal are represented in the singular spectrum by the smallest eigenvalues.. Separability he obtained SSA decomposition is a unction o the window length choice. his choice thereore conditions whether the components obtained will be correlated to a greater or lesser degree in the requency domain. o study whether these components are mutually independent, one deines the ollowing necessary (but not suicient) separability condition []: wo principal components F and F obtained rom the elementary matrices E and E are separable (w-orthogonal) i the inner product o series F and F is null, i.e.: N, = w i = ( F ) w ( i) ( i) F (4) where the weights wi are deined by i + or i L w i = L or L i < K (5) N i or K i N and L = min( L, K), and K = max( L, K) I the original series F is decomposed using SSA into a sum o separable components F, F,, F L that match equation (4) then this sum can be interpreted as an expansion o the original signal F with respect to a certain w-orthogonal basis generated by the original series itsel []. he weights in the inner product have the orm o a trapezium and reduce the inluence o data close to the end-points o the components with respect to the central terms in the series, particularly or large window lengths. Real-lie signals do not match equation (4). In practice, one speaks only about approximate separability. o quantiy the quality o separation between two components, one deines the weighted correlation or w-correlation o two components F and F : w ( ) ( ) (, ) i i i w F F w ρ = = (6) F F N N w w w ( i) ( i) w ( i) ( i) i N i his coeicient is a measure o the deviation o two series, F and F, rom w- orthogonality. I the absolute value o the w- correlation is near zero, then the two series are separable. I it is large (near ), the two components are badly separable. It is desirable to obtain a decomposition such that the principal components are mutually independent in order to extract the noise present in the displacement signal. Figure shows the SSA decomposition o a displacement signal (described in the results section) using a window length L =. Figure shows the decomposition obtained together with the original signal (see igure caption or details). Figure (b) shows a graphical representation o the w- correlation matrix corresponding to the decomposition obtained. Cell ij represents the correlation between components i and j, gray-scale w w coded rom black or ρ = to white or ρ = ij Fig.. (a) Singular spectrum. (b) w-correlation matrix. (c) Original time series and principal components obtained in the time domain. ij ISBN:
4 3 Cluster Analysis As was discussed above, the statistical independence o the principal components depends on the size o the window length chosen. With large values ( L 5) while the separation is good, the number o principal components obtained or the representation o the trend signal is impractically large. In order to reduce the number o principal components with which to extract the trend signal and ensure their statistical independence, we apply a cluster analysis procedure. Namely, we perorm K -means clustering on the w-correlation matrix in order to obtain K groups o independent components. he K -means algorithm classiies a given data set through a certain number K o clusters ixed a priori []. he main idea is to deine K centroids, one or each cluster. Since dierent locations o these centroids cause dierent results, the best choice is to place them as ar as possible rom each other. he next step is to take each point belonging to a given data set and associate it to the nearest centroid. When no point is let, the irst step is completed and a putative grouping is made. At this point, one needs to recalculate K new centroids as barycenters o the clusters resulting rom the previous step. Ater one has these K new centroids, a new binding procedure has to be carried out between the same data set points and the nearest new centroid. A loop has thus been generated. As a result o this loop, the K centroids change their location step by step until no more changes are made, i.e., the centroids do not move any more. Finally, the algorithm aims at minimizing an objective unction, in this case the ollowing squared error unction: k n ( j) J = x i c j (7) j= ( j) Where i c j x, a distance measure ( j) between a data point x i and the cluster centre c j, is an indicator o the distance o the n data points rom their respective cluster centers. his simple version o the K -means procedure can be viewed as a greedy algorithm or partitioning the n samples into K clusters so as to minimize the sum o the squared distances to the cluster centers. Unortunately there is no general theoretical solution to ind the optimal number o clusters or any given data set. A simple and systematic heuristic approach to obtain an optimal number o clusters is to compare the obtained error unction (Equation 7) o multiples runs with dierent K classes and chose the number o clusters that minimize this unction. 4 Smoothing Automation he SSA smoothing procedure presented in Alonso et al. [] is based on the act that raw-displacement acquired signals present a very large signal-to-noise ratio. In this situation, the contribution o the irst matrices to the norm o X is much higher than the contribution o the last matrices, that represent noise. o eliminate the noise present in the displacement signal it is suicient to choose the leading eigenvalues that represent a large percentage o the entire singular spectrum. As it has been pointed out [-], one o the drawbacks o SSA application is the lack o general rules or selecting the values o the parameters L and r that arise in the SSA algorithm. Moreover, certain window lengths and grouping strategy choices produce a poor separation between signal trend and noise. In other words, trend components would be mixed with noise components in the reconstruction o the signal. A way to overcome the uncertainty in the choice o the truncation value r is to apply sequential SSA. his means that we extract some components o the initial series by the standard SSA and then extract the components o interest by applying SSA smoothing to an already smoothed record. Such a recursive SSA application produces a gradual elimination o the noise present in the signal. o ensure that any signiicant part o the noise is eliminated, the number o eigenvalues L r to eliminate in each iteration was chosen to satisy the ollowing criteria: log < log λ L λ λ r λ r L (8) his criterion ensures that the eigenvalues whose λr logarithmic dierence log is lower than the λr average logarithmic range o the entire singular λ spectrum: log are eliminated. his criterion L λ L ensures that we eliminate eigenvalues in a zone where the singular spectrum has suicient latness. he convergence o the sequential procedure may be measured by means o the percentage root mean square (RMS) dierence between the current and previous acceleration signals obtained in each iteration. he algorithm stops when this dierence ISBN:
5 is suiciently small, namely the stop criterion at iteration i is: RMS ( G i Gi ) E = < ε RMS( G ) = (9) i Where G i is the smoothed acceleration reconstructed at iteration i. he automatic iltering procedure is summarized in the ollowing way: Choose an arbitrary window length L. Perorm 4 -means clustering on the w-correlation matrix in order to obtain 3 groups o independent components and extract the trend signal. Apply sequential SSA to the extracted trend signal (truncation value r is ixed to account or the grouping criterion, Equation 8). Calculate the acceleration signal numerically in each iteration. Stop the procedure according to the stop criterion (Equation 9). Fig.. (a) Experimental layout. (b). Noisy vertical displacement signal acquired. Figure 3 shows the w-correlation matrices obtained using L = 5, L =, L = 5 and L =. It is clear rom Figure 3 that the smaller window lengths do not adequately separate the high and low requency components. he greater window lengths, however, yield a good separation o the high and low requencies. 5 Results In order to study the perormance o the iltering procedure three signals were tested. wo reerence acceleration signals were taken rom the literature. A double dierentiation is perormed on each signal in order to quantiy the eect o smoothing. First order central inite dierences are used to calculate the higher derivatives. he sequential SSA algorithm stops when the percentage dierence between the current and previous values o the RMS acceleration is smaller than % in order to prevent excessive smoothing. he irst signal is the motion o a vertical slider moved by hand to obtain a non-stationary mono-dimensional motion o biomechanical origin. A subject was asked to move the slider randomly with ast upward and downward movements. he vertical position (Figure ) was acquired with the use o a marker attached to the slider and the Qualisys camera system (Qualisys Medical AB). he second derivative o this record, or o the record obtained ater smoothing, is compared to the acceleration obtained directly rom an accelerometer attached to the slider. Displacement and acceleration were sampled at Hz during 5.9 seconds, obtaining records o 8 elements each. SSA decomposition was attained using a window length L = (see Figure ). Various window lengths were tested in order to study the inluence o this parameter on the obtained decomposition. Fig. 3. w-correlation matrices. (a) L = 5. (b) L =. (c) L = 5. (d) L =. Ater the SSA decomposition, a 4-means cluster analysis was perormed in order to obtain our requency independent components or each vibration signal. he principal components orming the low-requency signal are I = [,,3]. It is important to stress that our components would also have obtained using a window length o L = 4 in a SSA decomposition o the original series, but without perorming the cluster analysis, these our components would not, however, have been independent in requency, which is the undamental point o this smoothing method. In order to eliminate the noise present in the obtained reconstruction sequential SSA-cluster analysis was applied. Such a recursive SSA application produces a gradual elimination o the ISBN:
6 noise present in the reconstructed signal. Figure 4 shows the results obtained or two iterations o the method using L =. Figure 4 (top) shows the acceleration calculated rom original raw displacement data (dotted line) and acceleration measured by the accelerometer (continuous line). Figure 4 (middle) shows the acceleration calculated ater having been passed through a 4 Hz cut-o requency Butterworth ilter (dotted line) and acceleration measured by the accelerometer (continuous line). his igure shows the large endpoint errors associated to the usage o this ilter. Finally, Figure 4 (bottom) represents the acceleration calculated ater smoothing the raw displacement signal using sequential SSA-cluster analysis ( L = ) and acceleration measured by the accelerometer (continuous line). Eigenvalue 6 Singular spectrum Eigenvalue number Fig. 5. Singular Spectrum evolution or the irst signal able summarizes the results obtained using several window lengths. It can be concluded that the results are robust to variations o the window length L and number o clusters. L RMSE N. iterations able. Summary o the results or several window lengths. Fig. 4. Signal results. See descriptions in the text. Figure 5 shows the singular spectrum evolution in each o the two iterations required to achieve the convergence criterion. It is clear rom the igure that the proposed procedure sequentially reduces the noise amplitude in each iteration. Comparison with traditional and advanced iltering techniques is quantiied by taking the root mean square o the error signal (RMSE) o the acceleration. he error signal is the dierence between the measured reerence signal (accelerometer) and the signal obtained ater iltering and dierentiating the raw data. he errors in terms o RMSE are 9.54 m / s or the Butterworth ilter and.8 m / s using SSA-cluster analysis respectively. he second signal corresponds to the woman ile rom GAILAB (Vaughan et al., 99). he signal measures the lateral displacement o a marker attached on the right tibial tubercle. Data was acquired or.94 seconds using a sampling requency o 5 Hz. he reerence signal is taken as that obtained ater iltering the raw data at 6.5 Hz, but dierent amplitude white noises are added in order to test smoothing methods. An amplitude, time generated, white noise was used or comparisons, see Giakas and Baltzopulos, 997 [5]. For signal, using L = the automatic procedure achieved RMSE = mm / s versus RMSE = 4 mm / s obtained by Giakas and Baltzopoulos, [5] using Power Spectrum Assessment (PSA) and RMSE = mm / s using a SSA non-automatic approach. he accuracy o the obtained acceleration may also be appreciated in Figure 6 (bottom), where the acceleration obtained rom the reerence signal is plotted along with the ISBN:
7 acceleration calculated rom SSA smoothed displacement data. he SSA algorithm decomposes the original signal into independent additive components o decreasing weight. his act allows the method to successully extract the latent trend in the signal rom the random noise inherent to the motion capture system. An automatic heuristic procedure based in sequential SSA and cluster analysis to extract the trend signal has been presented. Fig. 6. (op) Singular Spectrum evolution. (Bottom) Acceleration obtained rom reerence signal (continuous line) and acceleration calculated rom SSA-smoothed displacement data (dotted line). he third signal is a measure o the angular coordinate o a pendulum impacting against a compliant wall [4]. he angular acceleration obtained rom the motion capture system is compared to that obtained directly (ater dividing by pendulum length) rom accelerometers. hree accelerometers were used in order to average their measurements to reduce noise. he average signal, logged at a sampling rate o 5 Hz is used as the acceleration reerence signal. he same procedure was perormed on the third signal. Using L = 5 the automatic procedure achieved RMSE = 4.37 rad / s (Figure 7). he result is similar to the value RMSE = 3.6 rad / s obtained by Giakas et al. [8] with the help o the Wigner distribution. he slight loss o accuracy in the SSA method is compensated by the ease with which the method is applied and the act that one must not extend the ends o the record in order to eliminate end point errors. Moreover, the automatic iltering procedure obtain similar results or a reasonable range o window lengths. 6 Conclusions he present work has studied the applicability o the sequential SSA method and cluster analysis to automatic smoothing o biomechanical kinematic signals. Fig. 7. (op) Singular Spectrum evolution. (Bottom) Acceleration obtained rom reerence signal 3 (continuous line) and acceleration calculated rom SSA-smoothed displacement data (dotted line). he method does not use any inormation rom the reerence acceleration signal (which is not available in practice) to perorm the smoothing. One o the main advantages to the method is the act that the algorithm requires the selection o just one parameter. Namely, the window length L, moreover, the results are robust to variations o the window length. In conclusion, we believe that this new automatic smoothing technique, that has proven its eectiveness with complex signals, will help to improve the accuracy o raw kinematic data processing in biomechanical analysis. Reerences: [] C.L. Vaughan, Smoothing and dierentiation o displacement-time data: an application o splines and digital iltering. International Journal o Bio-Medical Computing, No. 3, 98, pp [] H. Hatze, he undamental problem o myoskeletal inverse dynamics and its implications, Journal o Biomechanics, No. 35,, pp [3] L. Chiari, U. Della Croce, Leardini, A. Cappozzo A., Human movement analysis using stereophotogrammetry Part : Instrumental ISBN:
8 errors, Gait and Posture, No., 5, pp.97-. [4] J. Dowling, A modelling strategy or the smoothing o biomechanical data. In: Johnsson, B. (Ed.), Biomechanics, Vol. XB, Human Kinetics, Champaign, IL, 63 67, 985. [5] G. Giakas, V. Baltzopoulos, A comparison o automatic iltering techniques applied to biomechanical walking data, Journal o Biomechanics, No. 3, 997, pp [6] R.I. Adham, S.A. Shihab, Discrete wavelet transorm: a tool in smoothing kinematic data, Journal o Biomechanics, No.3, 999, pp [7] A. Georgiakis, L.K. Stergioulas, G. Giakas, Wigner iltering with smooth roll-o boundary or dierentiation o noisy non-stationary signals, Signal Processing, No. 8,, pp [8] G. Giakas, L.K. Stergioulas, A. Vourdas, ime-requency analysis and iltering o kinematic signals with impacts using the Wigner unction: accurate estimation o the second derivative, Journal o Biomechanics, No.33,, pp [9] K.S. Erer, Adaptive usage o the Butterworth digital ilter, Journal o Biomechanics, No. 4, 7, pp [] F. J. Alonso, J. M. Del Castillo, P. Pintado, Application o singular spectrum analysis to the smoothing o raw kinematic signals, Journal o Biomechanics, No.38, 5, pp [] N. Golyandina, V. Nekrutkin, A. Zhigljavsky, Analisys o time series structure - SSA and related techniques, Chapman & Hall/CRC, Boca Raton, Florida,. [] J.B. MacQueen, Some methods or classiication and analysis o multivariate observations, in: Proceedings o the Fith Berkeley Symposium on Mathematical Statistics and Probability, vol.. Berkeley, CA, 967, pp ISBN:
MAPI Computer Vision. Multiple View Geometry
MAPI Computer Vision Multiple View Geometry Geometry o Multiple Views 2- and 3- view geometry p p Kpˆ [ K R t]p Geometry o Multiple Views 2- and 3- view geometry Epipolar Geometry The epipolar geometry
More informationSensors & Transducers 2016 by IFSA Publishing, S. L.
Sensors & ransducers 06 by IFSA Publishing, S. L. http://www.sensorsportal.com Polynomial Regression echniques or Environmental Data Recovery in Wireless Sensor Networks Kohei Ohba, Yoshihiro Yoneda, Koji
More informationChapter 3 Image Enhancement in the Spatial Domain
Chapter 3 Image Enhancement in the Spatial Domain Yinghua He School o Computer Science and Technology Tianjin University Image enhancement approaches Spatial domain image plane itsel Spatial domain methods
More informationCompressed Sensing Image Reconstruction Based on Discrete Shearlet Transform
Sensors & Transducers 04 by IFSA Publishing, S. L. http://www.sensorsportal.com Compressed Sensing Image Reconstruction Based on Discrete Shearlet Transorm Shanshan Peng School o Inormation Science and
More informationA Comparison of Parameters Estimation Based on Forecasting Accuracy in Singular Spectrum Analysis (SSA)
A Comparison of Parameters Estimation Based on Forecasting Accuracy in Singular Spectrum Analysis (SSA) Meri Andani Rosmalawati, Iman Bimantara M 2, Gumgum Darmawan, Toni Toharudin Statistics Department,
More informationCS485/685 Computer Vision Spring 2012 Dr. George Bebis Programming Assignment 2 Due Date: 3/27/2012
CS8/68 Computer Vision Spring 0 Dr. George Bebis Programming Assignment Due Date: /7/0 In this assignment, you will implement an algorithm or normalizing ace image using SVD. Face normalization is a required
More informationAutomatic filtering techniques for three-dimensional kinematics data using 3D motion capture system
Automatic filtering techniques for three-dimensional kinematics data using 3D motion capture system Aissaoui R., Husse S., Mecheri H., Parent G., de Guise J.A. Laboratoire de recherche en imagerie et orthopédie,
More information1. Techniques for ill-posed least squares problems. We write throughout this lecture = 2. Consider the following least squares problems:
ILL-POSED LEAST SQUARES PROBLEMS. Techniques for ill-posed least squares problems. We write throughout this lecture = 2. Consider the following least squares problems: where A = min b Ax, min b A ε x,
More informationNeighbourhood Operations
Neighbourhood Operations Neighbourhood operations simply operate on a larger neighbourhood o piels than point operations Origin Neighbourhoods are mostly a rectangle around a central piel Any size rectangle
More informationBinary Morphological Model in Refining Local Fitting Active Contour in Segmenting Weak/Missing Edges
0 International Conerence on Advanced Computer Science Applications and Technologies Binary Morphological Model in Reining Local Fitting Active Contour in Segmenting Weak/Missing Edges Norshaliza Kamaruddin,
More informationResearch on Image Splicing Based on Weighted POISSON Fusion
Research on Image Splicing Based on Weighted POISSO Fusion Dan Li, Ling Yuan*, Song Hu, Zeqi Wang School o Computer Science & Technology HuaZhong University o Science & Technology Wuhan, 430074, China
More informationMethod estimating reflection coefficients of adaptive lattice filter and its application to system identification
Acoust. Sci. & Tech. 28, 2 (27) PAPER #27 The Acoustical Society o Japan Method estimating relection coeicients o adaptive lattice ilter and its application to system identiication Kensaku Fujii 1;, Masaaki
More informationA Proposed Approach for Solving Rough Bi-Level. Programming Problems by Genetic Algorithm
Int J Contemp Math Sciences, Vol 6, 0, no 0, 45 465 A Proposed Approach or Solving Rough Bi-Level Programming Problems by Genetic Algorithm M S Osman Department o Basic Science, Higher Technological Institute
More informationFig. 3.1: Interpolation schemes for forward mapping (left) and inverse mapping (right, Jähne, 1997).
Eicken, GEOS 69 - Geoscience Image Processing Applications, Lecture Notes - 17-3. Spatial transorms 3.1. Geometric operations (Reading: Castleman, 1996, pp. 115-138) - a geometric operation is deined as
More informationECE 6560 Multirate Signal Processing Chapter 8
Multirate Signal Processing Chapter 8 Dr. Bradley J. Bazuin Western Michigan University College o Engineering and Applied Sciences Department o Electrical and Computer Engineering 903 W. Michigan Ave.
More informationMATRIX ALGORITHM OF SOLVING GRAPH CUTTING PROBLEM
UDC 681.3.06 MATRIX ALGORITHM OF SOLVING GRAPH CUTTING PROBLEM V.K. Pogrebnoy TPU Institute «Cybernetic centre» E-mail: vk@ad.cctpu.edu.ru Matrix algorithm o solving graph cutting problem has been suggested.
More informationInternational Journal of Foundations of Computer Science c World Scientic Publishing Company DFT TECHNIQUES FOR SIZE ESTIMATION OF DATABASE JOIN OPERA
International Journal of Foundations of Computer Science c World Scientic Publishing Company DFT TECHNIQUES FOR SIZE ESTIMATION OF DATABASE JOIN OPERATIONS KAM_IL SARAC, OMER E GEC_IO GLU, AMR EL ABBADI
More informationJane Li. Assistant Professor Mechanical Engineering Department, Robotic Engineering Program Worcester Polytechnic Institute
Jane Li Assistant Professor Mechanical Engineering Department, Robotic Engineering Program Worcester Polytechnic Institute (3 pts) Compare the testing methods for testing path segment and finding first
More informationDimension reduction : PCA and Clustering
Dimension reduction : PCA and Clustering By Hanne Jarmer Slides by Christopher Workman Center for Biological Sequence Analysis DTU The DNA Array Analysis Pipeline Array design Probe design Question Experimental
More informationMotion Estimation. There are three main types (or applications) of motion estimation:
Members: D91922016 朱威達 R93922010 林聖凱 R93922044 謝俊瑋 Motion Estimation There are three main types (or applications) of motion estimation: Parametric motion (image alignment) The main idea of parametric motion
More informationAutomatic Singular Spectrum Analysis for Time-Series Decomposition
Automatic Singular Spectrum Analysis for Time-Series Decomposition A.M. Álvarez-Meza and C.D. Acosta-Medina and G. Castellanos-Domínguez Universidad Nacional de Colombia, Signal Processing and Recognition
More informationMultiple attenuation with a modified parabolic Radon transform B. Ursin*, B. Abbad, NTNU, Trondheim, Norway, M. J. Porsani, UFBA, Salvador, Brazil.
Multiple attenuation with a modiied parabolic Radon transorm B. Ursin* B. Abbad TU Trondheim orway M. J. Porsani UFBA Salvador Brazil. Copyright 009 SBG - Sociedade Brasileira de Geoísica This paper was
More informationReducing the Bandwidth of a Sparse Matrix with Tabu Search
Reducing the Bandwidth o a Sparse Matrix with Tabu Search Raael Martí a, Manuel Laguna b, Fred Glover b and Vicente Campos a a b Dpto. de Estadística e Investigación Operativa, Facultad de Matemáticas,
More informationWhat is Clustering? Clustering. Characterizing Cluster Methods. Clusters. Cluster Validity. Basic Clustering Methodology
Clustering Unsupervised learning Generating classes Distance/similarity measures Agglomerative methods Divisive methods Data Clustering 1 What is Clustering? Form o unsupervised learning - no inormation
More informationUnderstanding Signal to Noise Ratio and Noise Spectral Density in high speed data converters
Understanding Signal to Noise Ratio and Noise Spectral Density in high speed data converters TIPL 4703 Presented by Ken Chan Prepared by Ken Chan 1 Table o Contents What is SNR Deinition o SNR Components
More informationFactorization with Missing and Noisy Data
Factorization with Missing and Noisy Data Carme Julià, Angel Sappa, Felipe Lumbreras, Joan Serrat, and Antonio López Computer Vision Center and Computer Science Department, Universitat Autònoma de Barcelona,
More informationRobust Kernel Methods in Clustering and Dimensionality Reduction Problems
Robust Kernel Methods in Clustering and Dimensionality Reduction Problems Jian Guo, Debadyuti Roy, Jing Wang University of Michigan, Department of Statistics Introduction In this report we propose robust
More informationELEC Dr Reji Mathew Electrical Engineering UNSW
ELEC 4622 Dr Reji Mathew Electrical Engineering UNSW Review of Motion Modelling and Estimation Introduction to Motion Modelling & Estimation Forward Motion Backward Motion Block Motion Estimation Motion
More informationMotion based 3D Target Tracking with Interacting Multiple Linear Dynamic Models
Motion based 3D Target Tracking with Interacting Multiple Linear Dynamic Models Zhen Jia and Arjuna Balasuriya School o EEE, Nanyang Technological University, Singapore jiazhen@pmail.ntu.edu.sg, earjuna@ntu.edu.sg
More informationDigital Image Processing. Image Enhancement in the Spatial Domain (Chapter 4)
Digital Image Processing Image Enhancement in the Spatial Domain (Chapter 4) Objective The principal objective o enhancement is to process an images so that the result is more suitable than the original
More informationOptimal Segmentation and Understanding of Motion Capture Data
Optimal Segmentation and Understanding of Motion Capture Data Xiang Huang, M.A.Sc Candidate Department of Electrical and Computer Engineering McMaster University Supervisor: Dr. Xiaolin Wu 7 Apr, 2005
More information6. Dicretization methods 6.1 The purpose of discretization
6. Dicretization methods 6.1 The purpose of discretization Often data are given in the form of continuous values. If their number is huge, model building for such data can be difficult. Moreover, many
More informationA Novel Accurate Genetic Algorithm for Multivariable Systems
World Applied Sciences Journal 5 (): 137-14, 008 ISSN 1818-495 IDOSI Publications, 008 A Novel Accurate Genetic Algorithm or Multivariable Systems Abdorreza Alavi Gharahbagh and Vahid Abolghasemi Department
More informationFX Cadzow / SSA Random Noise Filter: Frequency Extension
FX Cadzow / SSA Random Noise Filter: Frequency Extension Alexander Falkovskiy*, Elvis Floreani, Gerry Schlosser alex@absoluteimaging.ca Absolute Imaging Inc., Calgary, AB, Canada Summary The application
More informationImage Coding with Active Appearance Models
Image Coding with Active Appearance Models Simon Baker, Iain Matthews, and Jeff Schneider CMU-RI-TR-03-13 The Robotics Institute Carnegie Mellon University Abstract Image coding is the task of representing
More informationLab # 2 - ACS I Part I - DATA COMPRESSION in IMAGE PROCESSING using SVD
Lab # 2 - ACS I Part I - DATA COMPRESSION in IMAGE PROCESSING using SVD Goals. The goal of the first part of this lab is to demonstrate how the SVD can be used to remove redundancies in data; in this example
More informationNonlinear Noise Reduction
Chapter 3 Nonlinear Noise Reduction Many natural and engineered systems generate nonlinear deterministic time series that are contaminated by random measurement or dynamic noise. Even clean time series
More informationMobile Robot Static Path Planning Based on Genetic Simulated Annealing Algorithm
Mobile Robot Static Path Planning Based on Genetic Simulated Annealing Algorithm Wang Yan-ping 1, Wubing 2 1. School o Electric and Electronic Engineering, Shandong University o Technology, Zibo 255049,
More informationITU - Telecommunication Standardization Sector. G.fast: Far-end crosstalk in twisted pair cabling; measurements and modelling ABSTRACT
ITU - Telecommunication Standardization Sector STUDY GROUP 15 Temporary Document 11RV-22 Original: English Richmond, VA. - 3-1 Nov. 211 Question: 4/15 SOURCE 1 : TNO TITLE: G.ast: Far-end crosstalk in
More informationModel parametrization strategies for Newton-based acoustic full waveform
Model parametrization strategies for Newton-based acoustic full waveform inversion Amsalu Y. Anagaw, University of Alberta, Edmonton, Canada, aanagaw@ualberta.ca Summary This paper studies the effects
More informationLecture 3: Camera Calibration, DLT, SVD
Computer Vision Lecture 3 23--28 Lecture 3: Camera Calibration, DL, SVD he Inner Parameters In this section we will introduce the inner parameters of the cameras Recall from the camera equations λx = P
More informationGENERATION OF DEM WITH SUB-METRIC VERTICAL ACCURACY FROM 30 ERS-ENVISAT PAIRS
GEERATIO OF DEM WITH SUB-METRIC VERTICAL ACCURACY FROM 30 ERS-EVISAT PAIRS C. Colesanti (), F. De Zan (), A. Ferretti (), C. Prati (), F. Rocca () () Dipartimento di Elettronica e Inormazione, Politecnico
More informationClustering: Classic Methods and Modern Views
Clustering: Classic Methods and Modern Views Marina Meilă University of Washington mmp@stat.washington.edu June 22, 2015 Lorentz Center Workshop on Clusters, Games and Axioms Outline Paradigms for clustering
More information1724. Mobile manipulators collision-free trajectory planning with regard to end-effector vibrations elimination
1724. Mobile manipulators collision-free trajectory planning with regard to end-effector vibrations elimination Iwona Pajak 1, Grzegorz Pajak 2 University of Zielona Gora, Faculty of Mechanical Engineering,
More informationUnsupervised Learning
Unsupervised Learning Learning without Class Labels (or correct outputs) Density Estimation Learn P(X) given training data for X Clustering Partition data into clusters Dimensionality Reduction Discover
More informationRecognition, SVD, and PCA
Recognition, SVD, and PCA Recognition Suppose you want to find a face in an image One possibility: look for something that looks sort of like a face (oval, dark band near top, dark band near bottom) Another
More informationImage Compression with Singular Value Decomposition & Correlation: a Graphical Analysis
ISSN -7X Volume, Issue June 7 Image Compression with Singular Value Decomposition & Correlation: a Graphical Analysis Tamojay Deb, Anjan K Ghosh, Anjan Mukherjee Tripura University (A Central University),
More informationApplication of Genetic Algorithms to the Identification of Website Link Structure
00 International Conerence on Advances in Social Networks Analysis and Mining Application o Genetic Algorithms to the Identiication o Website Link Structure M. R. Martínez Torres, B. Palacios, S. L. Toral,
More informationCS 775: Advanced Computer Graphics. Lecture 3 : Kinematics
CS 775: Advanced Computer Graphics Lecture 3 : Kinematics Traditional Cell Animation, hand drawn, 2D Lead Animator for keyframes http://animation.about.com/od/flashanimationtutorials/ss/flash31detanim2.htm
More informationReview for Exam I, EE552 2/2009
Gonale & Woods Review or Eam I, EE55 /009 Elements o Visual Perception Image Formation in the Ee and relation to a photographic camera). Brightness Adaption and Discrimination. Light and the Electromagnetic
More informationCHAPTER 5 MOTION DETECTION AND ANALYSIS
CHAPTER 5 MOTION DETECTION AND ANALYSIS 5.1. Introduction: Motion processing is gaining an intense attention from the researchers with the progress in motion studies and processing competence. A series
More informationApproximate structural optimization using kriging method and digital modeling technique considering noise in sampling data
Available online at www.sciencedirect.com Computers and Structures 86 (2008) 1477 1485 www.elsevier.com/locate/compstruc Approimate structural optimiation using kriging method and digital modeling technique
More informationADVANCED IMAGE PROCESSING METHODS FOR ULTRASONIC NDE RESEARCH C. H. Chen, University of Massachusetts Dartmouth, N.
ADVANCED IMAGE PROCESSING METHODS FOR ULTRASONIC NDE RESEARCH C. H. Chen, University of Massachusetts Dartmouth, N. Dartmouth, MA USA Abstract: The significant progress in ultrasonic NDE systems has now
More informationUniversity of Florida CISE department Gator Engineering. Clustering Part 2
Clustering Part 2 Dr. Sanjay Ranka Professor Computer and Information Science and Engineering University of Florida, Gainesville Partitional Clustering Original Points A Partitional Clustering Hierarchical
More informationRough Connected Topologized. Approximation Spaces
International Journal o Mathematical Analysis Vol. 8 04 no. 53 69-68 HIARI Ltd www.m-hikari.com http://dx.doi.org/0.988/ijma.04.4038 Rough Connected Topologized Approximation Spaces M. J. Iqelan Department
More information3-D TERRAIN RECONSTRUCTION WITH AERIAL PHOTOGRAPHY
3-D TERRAIN RECONSTRUCTION WITH AERIAL PHOTOGRAPHY Bin-Yih Juang ( 莊斌鎰 ) 1, and Chiou-Shann Fuh ( 傅楸善 ) 3 1 Ph. D candidate o Dept. o Mechanical Engineering National Taiwan University, Taipei, Taiwan Instructor
More informationAn Adaptive Eigenshape Model
An Adaptive Eigenshape Model Adam Baumberg and David Hogg School of Computer Studies University of Leeds, Leeds LS2 9JT, U.K. amb@scs.leeds.ac.uk Abstract There has been a great deal of recent interest
More informationConvex Optimization / Homework 2, due Oct 3
Convex Optimization 0-725/36-725 Homework 2, due Oct 3 Instructions: You must complete Problems 3 and either Problem 4 or Problem 5 (your choice between the two) When you submit the homework, upload a
More informationA Quantitative Comparison of 4 Algorithms for Recovering Dense Accurate Depth
A Quantitative Comparison o 4 Algorithms or Recovering Dense Accurate Depth Baozhong Tian and John L. Barron Dept. o Computer Science University o Western Ontario London, Ontario, Canada {btian,barron}@csd.uwo.ca
More information9.8 Graphing Rational Functions
9. Graphing Rational Functions Lets begin with a deinition. Deinition: Rational Function A rational unction is a unction o the orm P where P and Q are polynomials. Q An eample o a simple rational unction
More informationClassifier Evasion: Models and Open Problems
Classiier Evasion: Models and Open Problems Blaine Nelson 1, Benjamin I. P. Rubinstein 2, Ling Huang 3, Anthony D. Joseph 1,3, and J. D. Tygar 1 1 UC Berkeley 2 Microsot Research 3 Intel Labs Berkeley
More informationCluster Analysis. Ying Shen, SSE, Tongji University
Cluster Analysis Ying Shen, SSE, Tongji University Cluster analysis Cluster analysis groups data objects based only on the attributes in the data. The main objective is that The objects within a group
More informationAnnotated multitree output
Annotated multitree output A simplified version of the two high-threshold (2HT) model, applied to two experimental conditions, is used as an example to illustrate the output provided by multitree (version
More informationA Robust and Efficient Motion Segmentation Based on Orthogonal Projection Matrix of Shape Space
A Robust and Efficient Motion Segmentation Based on Orthogonal Projection Matrix of Shape Space Naoyuki ICHIMURA Electrotechnical Laboratory 1-1-4, Umezono, Tsukuba Ibaraki, 35-8568 Japan ichimura@etl.go.jp
More informationScienceDirect. Using Search Algorithms for Modeling Economic Processes
vailable online at www.sciencedirect.com ScienceDirect Procedia Economics and Finance 6 ( 03 ) 73 737 International Economic onerence o Sibiu 03 Post risis Economy: hallenges and Opportunities, IES 03
More informationClustering and Visualisation of Data
Clustering and Visualisation of Data Hiroshi Shimodaira January-March 28 Cluster analysis aims to partition a data set into meaningful or useful groups, based on distances between data points. In some
More informationClassification Method for Colored Natural Textures Using Gabor Filtering
Classiication Method or Colored Natural Textures Using Gabor Filtering Leena Lepistö 1, Iivari Kunttu 1, Jorma Autio 2, and Ari Visa 1, 1 Tampere University o Technology Institute o Signal Processing P.
More informationGeneral Instructions. Questions
CS246: Mining Massive Data Sets Winter 2018 Problem Set 2 Due 11:59pm February 8, 2018 Only one late period is allowed for this homework (11:59pm 2/13). General Instructions Submission instructions: These
More informationLarger K-maps. So far we have only discussed 2 and 3-variable K-maps. We can now create a 4-variable map in the
EET 3 Chapter 3 7/3/2 PAGE - 23 Larger K-maps The -variable K-map So ar we have only discussed 2 and 3-variable K-maps. We can now create a -variable map in the same way that we created the 3-variable
More informationUnsupervised learning in Vision
Chapter 7 Unsupervised learning in Vision The fields of Computer Vision and Machine Learning complement each other in a very natural way: the aim of the former is to extract useful information from visual
More information5. GENERALIZED INVERSE SOLUTIONS
5. GENERALIZED INVERSE SOLUTIONS The Geometry of Generalized Inverse Solutions The generalized inverse solution to the control allocation problem involves constructing a matrix which satisfies the equation
More informationMotion Interpretation and Synthesis by ICA
Motion Interpretation and Synthesis by ICA Renqiang Min Department of Computer Science, University of Toronto, 1 King s College Road, Toronto, ON M5S3G4, Canada Abstract. It is known that high-dimensional
More informationThe Graph of an Equation Graph the following by using a table of values and plotting points.
Calculus Preparation - Section 1 Graphs and Models Success in math as well as Calculus is to use a multiple perspective -- graphical, analytical, and numerical. Thanks to Rene Descartes we can represent
More informationReddit Recommendation System Daniel Poon, Yu Wu, David (Qifan) Zhang CS229, Stanford University December 11 th, 2011
Reddit Recommendation System Daniel Poon, Yu Wu, David (Qifan) Zhang CS229, Stanford University December 11 th, 2011 1. Introduction Reddit is one of the most popular online social news websites with millions
More informationES 240: Scientific and Engineering Computation. a function f(x) that can be written as a finite series of power functions like
Polynomial Deinition a unction () that can be written as a inite series o power unctions like n is a polynomial o order n n ( ) = A polynomial is represented by coeicient vector rom highest power. p=[3-5
More informationA SEMI-FRAGILE WATERMARKING SCHEME FOR IMAGE AUTHENTICATION AND TAMPER-PROOFING
A SEMI-FRAGILE WATERMARKING SCHEME FOR IMAGE AUTHENTICATION AND TAMPER-PROOFING Nozomi Ishihara and Kôi Abe Department o Computer Science, The University o Electro-Communications, 1-5-1 Chougaoa Chou-shi,
More informationGeneralized trace ratio optimization and applications
Generalized trace ratio optimization and applications Mohammed Bellalij, Saïd Hanafi, Rita Macedo and Raca Todosijevic University of Valenciennes, France PGMO Days, 2-4 October 2013 ENSTA ParisTech PGMO
More informationReflection and Refraction
Relection and Reraction Object To determine ocal lengths o lenses and mirrors and to determine the index o reraction o glass. Apparatus Lenses, optical bench, mirrors, light source, screen, plastic or
More informationROBUST FACE DETECTION UNDER CHALLENGES OF ROTATION, POSE AND OCCLUSION
ROBUST FACE DETECTION UNDER CHALLENGES OF ROTATION, POSE AND OCCLUSION Phuong-Trinh Pham-Ngoc, Quang-Linh Huynh Department o Biomedical Engineering, Faculty o Applied Science, Hochiminh University o Technology,
More informationDimension Reduction CS534
Dimension Reduction CS534 Why dimension reduction? High dimensionality large number of features E.g., documents represented by thousands of words, millions of bigrams Images represented by thousands of
More informationIMAGE ANALYSIS, CLASSIFICATION, and CHANGE DETECTION in REMOTE SENSING
SECOND EDITION IMAGE ANALYSIS, CLASSIFICATION, and CHANGE DETECTION in REMOTE SENSING ith Algorithms for ENVI/IDL Morton J. Canty с*' Q\ CRC Press Taylor &. Francis Group Boca Raton London New York CRC
More informationSUPER RESOLUTION IMAGE BY EDGE-CONSTRAINED CURVE FITTING IN THE THRESHOLD DECOMPOSITION DOMAIN
SUPER RESOLUTION IMAGE BY EDGE-CONSTRAINED CURVE FITTING IN THE THRESHOLD DECOMPOSITION DOMAIN Tsz Chun Ho and Bing Zeng Department o Electronic and Computer Engineering The Hong Kong University o Science
More informationSmoothing Dissimilarities for Cluster Analysis: Binary Data and Functional Data
Smoothing Dissimilarities for Cluster Analysis: Binary Data and unctional Data David B. University of South Carolina Department of Statistics Joint work with Zhimin Chen University of South Carolina Current
More information1498. End-effector vibrations reduction in trajectory tracking for mobile manipulator
1498. End-effector vibrations reduction in trajectory tracking for mobile manipulator G. Pajak University of Zielona Gora, Faculty of Mechanical Engineering, Zielona Góra, Poland E-mail: g.pajak@iizp.uz.zgora.pl
More informationApplication of FX Singular Spectrum Analysis on Structural Data
Application of FX Singular Spectrum Analysis on Structural Data Alexander Falkovskiy*, Elvis Floreani, Gerry Schlosser alex@absoluteimaging.ca Absolute Imaging Inc., Calgary, AB, Canada Summary The application
More informationOptimum design of roll forming process of slide rail using design of experiments
Journal o Mechanical Science and Technology 22 (28) 1537~1543 Journal o Mechanical Science and Technology www.springerlink.com/content/1738-494x DOI 1.17/s1226-8-43-9 Optimum design o roll orming process
More information1. [1 pt] What is the solution to the recurrence T(n) = 2T(n-1) + 1, T(1) = 1
Asymptotics, Recurrence and Basic Algorithms 1. [1 pt] What is the solution to the recurrence T(n) = 2T(n-1) + 1, T(1) = 1 1. O(logn) 2. O(n) 3. O(nlogn) 4. O(n 2 ) 5. O(2 n ) 2. [1 pt] What is the solution
More informationRobotics Programming Laboratory
Chair of Software Engineering Robotics Programming Laboratory Bertrand Meyer Jiwon Shin Lecture 8: Robot Perception Perception http://pascallin.ecs.soton.ac.uk/challenges/voc/databases.html#caltech car
More informationCompressed Sensing Algorithm for Real-Time Doppler Ultrasound Image Reconstruction
Mathematical Modelling and Applications 2017; 2(6): 75-80 http://www.sciencepublishinggroup.com/j/mma doi: 10.11648/j.mma.20170206.14 ISSN: 2575-1786 (Print); ISSN: 2575-1794 (Online) Compressed Sensing
More informationCHAPTER 2 TEXTURE CLASSIFICATION METHODS GRAY LEVEL CO-OCCURRENCE MATRIX AND TEXTURE UNIT
CHAPTER 2 TEXTURE CLASSIFICATION METHODS GRAY LEVEL CO-OCCURRENCE MATRIX AND TEXTURE UNIT 2.1 BRIEF OUTLINE The classification of digital imagery is to extract useful thematic information which is one
More informationImage Transformation Techniques Dr. Rajeev Srivastava Dept. of Computer Engineering, ITBHU, Varanasi
Image Transformation Techniques Dr. Rajeev Srivastava Dept. of Computer Engineering, ITBHU, Varanasi 1. Introduction The choice of a particular transform in a given application depends on the amount of
More informationInterlude: Solving systems of Equations
Interlude: Solving systems of Equations Solving Ax = b What happens to x under Ax? The singular value decomposition Rotation matrices Singular matrices Condition number Null space Solving Ax = 0 under
More informationADAPTIVE LOW RANK AND SPARSE DECOMPOSITION OF VIDEO USING COMPRESSIVE SENSING
ADAPTIVE LOW RANK AND SPARSE DECOMPOSITION OF VIDEO USING COMPRESSIVE SENSING Fei Yang 1 Hong Jiang 2 Zuowei Shen 3 Wei Deng 4 Dimitris Metaxas 1 1 Rutgers University 2 Bell Labs 3 National University
More informationUsing a Projected Subgradient Method to Solve a Constrained Optimization Problem for Separating an Arbitrary Set of Points into Uniform Segments
Using a Projected Subgradient Method to Solve a Constrained Optimization Problem or Separating an Arbitrary Set o Points into Uniorm Segments Michael Johnson May 31, 2011 1 Background Inormation The Airborne
More information1 INTRODUCTION The LMS adaptive algorithm is the most popular algorithm for adaptive ltering because of its simplicity and robustness. However, its ma
MULTIPLE SUBSPACE ULV ALGORITHM AND LMS TRACKING S. HOSUR, A. H. TEWFIK, D. BOLEY University of Minnesota 200 Union St. S.E. Minneapolis, MN 55455 U.S.A fhosur@ee,tewk@ee,boley@csg.umn.edu ABSTRACT. The
More informationAN ALGORITHM FOR BLIND RESTORATION OF BLURRED AND NOISY IMAGES
AN ALGORITHM FOR BLIND RESTORATION OF BLURRED AND NOISY IMAGES Nader Moayeri and Konstantinos Konstantinides Hewlett-Packard Laboratories 1501 Page Mill Road Palo Alto, CA 94304-1120 moayeri,konstant@hpl.hp.com
More informationTwo-view geometry Computer Vision Spring 2018, Lecture 10
Two-view geometry http://www.cs.cmu.edu/~16385/ 16-385 Computer Vision Spring 2018, Lecture 10 Course announcements Homework 2 is due on February 23 rd. - Any questions about the homework? - How many of
More informationMSA220 - Statistical Learning for Big Data
MSA220 - Statistical Learning for Big Data Lecture 13 Rebecka Jörnsten Mathematical Sciences University of Gothenburg and Chalmers University of Technology Clustering Explorative analysis - finding groups
More informationSpectral Clustering and Community Detection in Labeled Graphs
Spectral Clustering and Community Detection in Labeled Graphs Brandon Fain, Stavros Sintos, Nisarg Raval Machine Learning (CompSci 571D / STA 561D) December 7, 2015 {btfain, nisarg, ssintos} at cs.duke.edu
More informationEE 264: Image Processing and Reconstruction. Image Motion Estimation II. EE 264: Image Processing and Reconstruction. Outline
Peman Milanar Image Motion Estimation II Peman Milanar Outline. Introduction to Motion. Wh Estimate Motion? 3. Global s. Local Motion 4. Block Motion Estimation 5. Optical Flow Estimation Basics 6. Optical
More information