An Iteratively Reweighted Norm Algorithm for Total Variation Regularization
|
|
- Damian Grant
- 5 years ago
- Views:
Transcription
1 An Iteratively eweighted Norm Algorithm for Total Variation egularization Paul odríguez and Brendt Wohlberg Abstract Total Variation TV regularization has become a oular method for a wide variety of image restoration roblems, including denoising and deconvolution. ecently, a number of authors have noted the advantages, including suerior erformance with certain non-gaussian noise, of relacing the standard l data fidelity term with an l norm. We roose a simle but very flexible and comutationally efficient method, the Iteratively eweighted Norm algorithm, for minimizing a generalized TV functional which includes both the l -TV and and l -TV roblems. I. INTODUCTION Total Variation TV regularization was first introduced for image denoising [], and has since evolved into a more general tool for solving a wide variety of image restoration roblems, including deconvolution and inainting [], [3]. The l -TV regularized solution of the inverse roblem involving data b and forward linear oerator A the identity in the case of denoising, and a convolution for a deconvolution roblem, for examle is the minimum of the functional Ju = Au b + λ D x u + D y u where D x u + D y u is the total variation of u, and D x, D y denote the horizontal and vertical discrete derivative oerators resectively. While a variety of algorithms [4], [5] have been roosed to solve this otimization roblem, it remains a comutationally exensive task which can be rohibitively costly for large roblems and non-sarse forward oerator A. ecently, the modified TV functional with an l data fidelity term Ju = Au b + λ D x u + D y u has attracted attention [6], [7] due to a number of advantages, including suerior denoising erformance with salt and eer seckle noise [8]. The standard aroaches to solving roblem are not effective for roblem, for Paul odríguez is with T-7 Mathematical Modeling and Analysis, Los Alamos National Laboratory, Los Alamos, NM 87545, USA. rodrig@t7.lanl.gov, Tel: , ax: Brendt Wohlberg is with T-7 Mathematical Modeling and Analysis, Los Alamos National Laboratory, Los Alamos, NM 87545, USA. brendt@t7.lanl.gov, Tel: , ax: This work was carried out under the ausices of the National Nuclear Security Administration of the U.S. Deartment of Energy at Los Alamos National Laboratory under Contract No. DE-AC5-6NA5396 and was artially suorted by the NNSA s Laboratory Directed esearch and Develoment Program. which efficient time-erformance and memory reuirements algorithm develoment is not well advanced [7]. In this aer, we introduce the Iteratively eweighted Norm IN algorithm for solving the generalized TV functional Ju = Au b + λ D x u + D y u 3 for and we have found the algorithm to converge for smaller and values, but cannot rove that the algorithm converges in these cases. The IN algorithm is a simle but comutationally efficient and very flexible method which is cometitive with existing, well-established [4], [5] algorithms for l -TV, and while we have not comared it with all of the l -TV algorithms of which we are aware [7], [8], [9], [], [], it is significantly faster than those with which we have erformed comarisons, and in contrast to many of these other methods, easily allows the inclusion of forward oerator A for a more general inverse roblem than denoising. II. ITEATIVELY EWEIGHTED NOM APPOACH A. Previous elated Work The IN aroach is closely related to the Iteratively eweighted Least Suares ILS [], [3] method, which minimizes the l norm Ju = Au b 4 by aroximating it, within an iterative scheme, by a weighted l norm roblem, for which an algebraic solution is easily derived. It is also worth mentioning similar methods used for constructing sarse signal reresentations [4], [5]. A brief outline of the convergence roof for ILS rovides a useful introduction to the outline of the IN convergence roof resented in Section II-D. At each iteration k, define the functional Q k u = W k/ Au b + Ju k, where W k = diag Au k b, and u k is the current solution i.e., the solution obtained at the end of the revious iteration. It may be shown that Ju Q k u u Ju k = Q k u k Ju k = Q k u k, i.e. the uadratic functional Q k u is tangent to Ju at u = u k, where it is also an uer bound for Ju. As a
2 result [3], the iteration u k+ = A T W k A A T b, 5 which minimizes Q k u using the weights derived from the revious iteration, converges to the minimizer of Ju. When <, the definition W k = diag Au k b must be modified to avoid the ossibility of division by zero. or =, it may be shown see [6] that the choice W k n,n = r k ɛ n if r n k ɛ if r n k < ɛ, where r k = Au k b, and ɛ is a small ositive number, guarantees global convergence to the minimizer of ρ ɛ r n, n where ɛ ρ ɛ r n = rn if r n ɛ r n ɛ if r n > ɛ is the Huber function. We will address a similar uestion for the IN algorithm in the following sections. B. IN: Data idelity Term The data fidelity term of Euation 3 has the form of the ILS functional in Euation 4, and is handled in the same way, reresenting Au b by Au b W / with iteratively udated weights W. Since the choice of W as described in Section II-A gives infinite weights for zerovalued comonents of Au b, when <, we set W = diag f Au b, where f is defined for some small ɛ as x if x > ɛ f x = ɛ if x ɛ, which is a common aroach for ILS algorithms [7]. C. IN: egularization Term It is not uite as obvious how to exress the TV regularization term from Euation 3 as a weighted l norm. Given vectors ξ and χ we have using block-matrix notation so that when we have W / W / ξ χ = k W = diag ξ + χ / W / W / ξ χ w k ξ k + w k χ k = ξ + χ. We therefore set ξ = D x u, χ = D y u, and define the oerator D and weights W D = Dx D y W = W W so that W / Du = W / D xu + W / D yu with weights defined by W = diag Dx u + D y u / gives the desired term. Note that this is not the searable aroximation, as in [8], for examle, which is often used. As in the case of the data fidelity term, care needs to be taken when < and D x u + D y u has zero-valued comonents. We define x / if x > ɛ f x = if x ɛ for some small ɛ, and set W = diag f Dx u + D y u. Note that f sets values smaller than the threshold, ɛ, to zero, as oosed to f, which sets values smaller than the threshold, ɛ, to ɛ. The motivation for this choice is that a region with very small or zero gradient should be allowed to have zero contribution to the regularization term, rather than be clamed to some minimum value. In ractice, however, we have found that this choice does not give significantly different results than the standard ILS aroach reresented by f. D. Convergence of the IN algorithm Here we briefly sketch the roof of global convergence of the IN algorithm. Combining the terms described in Sections II-B and II-C, we define the generalized TV functional Q k u = W k/ Au b u k + λ + λ k/ W Du + u k, 6 where W k, W k, and uk are the current data fidelity weights, regularization weights, and solution resectively. The constant with resect to u terms u k = Au k b u k = D x u k + D y u k do not effect the IN iteration since they have zero gradient with resect to u, but are necessary so that Ju k = Q k u k. Note that Ju k = uk + λ uk, as in Euation 3. It is easily shown that Ju Q k u u Ju k = Q k u k Ju k = Q k u k Q k u = A T W k A + λdt k W D >
3 Algorithm IN algorithm - General case Initialize u = A T A + λd T D A T b Iterate W k W k = diag f Au k b = diag f D x u k + D y u k u k = A T W k +λdy T W k D y A + λdt x W k D x A T W k b ig Lena image. As in the ILS case Section II-A, the uadratic functional Q k u is tangent to Ju at u = u k, where it is also an uer bound for Ju, and it has a ositive definite Hessian. Using these results, it can be shown that the minimizer of Q k u, given by u = A T W k A + λdt W k D A T W k b, 7 converges to the minimizer of 3 as we iterate over k. III. IN ALGOITHMS The IN algorithm can be easily derived see Algorithm from Euation 7. The matrix inversion can be achieved using the Conjugate Gradient CG or Preconditioned CG PCG method such as Jacobi line relaxation JL or symmetric Gauss-Seidel line relaxation SLGS [9]. We have also found that a significant seed imrovement may be achieved by starting with a high CG/PCG tolerance, which is decreased with each main iteration until the final desired value is reached. or the results reorted hereafter, oerators D x and D y are defined by alying the same one-dimensional discrete derivative along image rows and columns resectively. Alied to vector u N, this discrete one-dimensional derivative is comuted as u k u k+ for the derivative at index k,,..., N }, and as u N 3 8u N / for the derivative at index N. Values in the ranges to 4 and 4 to 8 were used for constants ɛ and ɛ resectively. All rogram run times were obtained using the NUMIPAD library [] on a Linux 3.GHz Intel Pentium-4 rocessor. In the remainder of this aer we shall restrict our attention to the l -TV case =, =, but note that this flexible aroach is caable of efficiently solving other cases as well, including the standard l -TV case =, = where we have found it to have similar erformance to the lagged diffusivity algorithm [4]; for a 5 5 inut image igure, normalized between and, corruted with Gaussian noise σ =.5 lagged diffusivity reuired.69s for 5 iterations, giving a denoised image with 6.4 db, while the l -TV IN via PCG-SLGS reuired.69s for 5 iterations, giving a 6. db outut. The number of CG iterations er main loo needed by lagged-diffusivity and several flavors of the IN algorithm are shown in igure. CG Iterations l TV Lagged l TV IN CG l TV IN PCG JL l TV IN PCG SLGS Main Iteration Number ig.. A comarison of CG iterations for lagged-diffusivity and several flavors of the IN algorithm for l -TV denoising. The IN algorithm is slow when couled with a standard CG solver, but the use of reconditioned CG PCG strategies, such as Jacobi line relaxation JL and symmetric Gauss- Seidel line relaxation SLGS, dramatically reduces the number of reuired CG iterations. A. IN l -TV Denoising Direct alication of Euation 7 for l -TV denoising is slow due to the ill-conditioning of the system to be inverted. However, in this case we may aly the substitution ũ = W k/ ũ = u to Euation 7, giving I + λw k / D T k W / DW k / W k/ b. Alying this modification indirect IN algorithm to the general direct algorithm results in a significant reduction in the reuired number of CG iterations; moreover if a PCG strategy is alied e.g. JL or SLGS, the number of CG iterations is further reduced, as shown in igure 4, where a comarison between the number of CG iteration for 5 loos for the direct and indirect CG, PCG-JL and PCG-SLGS IN algorithms is carried out to denoise a 5 5 image igure corruted with 5% seckle noise. The run times were 8.s,.s,.3s
4 a Image with % seckle noise. SN: -.76dB. b Denoised l -TV λ =.5. SN: 7.4dB. Time:.63s c Image with 5% seckle noise. SN: -5.77dB. d Denoised l -TV λ =.5. SN: 9.6dB. Time:.6s ig. 3. Extreme examles for l -TV denoising using the IN algorithm: original image igure was corruted with % and 5% seckle noise. and.3s resectively, with denoised ualities of 7.36 db,.85 db,.35 db, and.36 db resectively. More extreme examles are resented in igure 3, with the original image corruted by % igure 3a and 5% igure 3c seckle noise. The result of l -TV denoising these images, using the IN via PCG-JL algorithm, is surrisingly good: 7.4 db igure 3b in.63s when denoising igure 3a and 9.6 db igure 3d in.6s when denoising igure 3c. B. IN l -TV Deconvolution We aly the IN algorithm to the roblem of deconvolution of an image convolved by a searable smoothing filter having 9 tas, aroximating a Gaussian with standard deviation of.. In this case A is the corresonding linear oerator, and the substitution alied in the revious section is no longer ossible. We constructed a test image by convolving the image in igure by the smoothing kernel, and adding 5% seckle noise see igure 5a, giving an image with an SN of 4.dB. Comaring the erformance of l -TV via the lagged diffusivity algorithm and l -TV via the IN method deconvolution, we obtain a.db reconstruction in 9.s and a 4.4dB reconstruction in 6.8s resectively see igure 5b. IV. CONCLUSIONS The Iteratively eweighted Norm IN aroach rovides a simle but comutationally efficient method for TV regularized otimization roblems, including both denoising and those such as deconvolution having a linear oerator in the data fidelity term. This framework is very flexible, and can be alied to regularized inversions with a wide variety of norms for the data fidelity and regularization terms, including the standard l -TV, and more recently roosed l -TV formulations. This method rovides a significantly faster algorithm for the l -TV formulation than any other algorithm of which we are aware. The NUMIPAD library [], is roving to be a useful tool for solving otimization roblems and can be used to reroduce the results resented in this aer.
5 8 6 4 IN l direct IN l indirect CG IN l indirect PCG JL IN l indirect PCG SLGS CG Iterations Main Iteration Number ig. 4. A comarison of CG iterations for the direct and several flavors of the indirect IN algorithm for l -TV denoising. Both reconditioned CG PCG strategies, Jacobi line relaxation JL and symmetric Gauss-Seidel line relaxation SLGS, dramatically reduced the number of reuired CG iterations. a Blurred image with 5% seckle noise. SN: 4.dB. ACKNOWLEDGMENT The authors thank Markus Berndt for valuable discussion on CG reconditioning methods. EEENCES [] L. udin, S. Osher, and E. atemi, Nonlinear total variation based noise removal algorithms, Physica D, vol. 6, , 99. [] C.. Vogel and M. E. Oman, ast, robust total variation-based reconstruction of noisy, blurred images, IEEE Transactions on Image Processing, vol. 7, , June 998. [3] T. Chan, S. Esedoglu,. Park, and A. Yi, ecent develoments in total variation image restoration, in The Handbook of Mathematical Models in Comuter Vision N. Paragios, Y. Chen, and O. augeras., eds., Sringer, 5. [4] C.. Vogel and M. E. Oman, Iterative methods for total variation denoising, SIAM Journal on Scientific Comuting, vol. 7,. 7 38, Jan [5] A. Chambolle, An algorithm for total variation minimization and alications, Journal of Mathematical Imaging and Vision, vol., , 4. [6] M. Nikolova, Minimizers of cost-functions involving nonsmooth datafidelity terms. alication to the rocessing of outliers, SIAM Journal on Numerical Analysis, vol. 4, no. 3, ,. [7] T.. Chan and S. Esedoglu, Asects of total variation regularized l function aroximation, SIAM Journal on Alied Mathematics, vol. 65, no. 5, , 5. [8] M. Nikolova, A variational aroach to remove outliers and imulse noise, Journal of Mathematical Imaging and Vision, vol., no. -,. 99, 4. 4th Conference on Mathematics and Image Analysis; Setember -3, ; Paris, rance. [9] J. Darbon and M. Sigelle, A fast and exact algorithm for total variation minimization, Lecture Notes in Comuter Science, vol. 35/5, , 5. [] D. Goldfarb and W. Yin, Second-order cone rogramming methods for total variation-based image restoration, SIAM J. Sci. Comut., vol. 7, no., , 5. [] J.-. Aujol, G. Gilboa, T. Chan, and S. Osher, Structure-texture image decomosition - modeling, algorithms, and arameter selection, International Journal of Comuter Vision, vol. 67, no.,. 36, 6. []. Wolke and H. Schwetlick, Iteratively reweighted least suares: Algorithms, convergence analysis, and numerical comarisons, SIAM Journal on Scientific and Statistical Comuting, vol. 9,. 97 9, Set [3] K. P. Bube and. T. Langan, Hybrid l /l minimization with alications to tomograhy, Geohysics, vol. 6, , July-August 997. [4] I.. Gorodnitsky and B. D. ao, A new iterative weighted norm minimization algorithm and its alications, in IEEE Sixth SP Worksho on Statistical Signal and Array Processing, 99, Victoria, BC, Canada, Oct. 99. b Deblurred l -TV λ =.3. SN: 4.4dB. Time: 6.8s ig. 5. l -TV deconvolution using the IN algorithm. The original image igure was convolved by a searable smoothing filter before the addition of 5% seckle noise. [5] B. D. ao and K. Kreutz-Delgado, An affine scaling methodology for best basis selection, IEEE Transactions On Signal Processing, vol. 47,. 87, Jan [6] S. A. uzinsky and E. T. Olsen, L and L minimization via a variant of Karmarkar s algorithm, IEEE Transactions on Acoustics, Seech and Signal Processing, vol. 37, , 989. [7] J. A. Scales and A. Gersztenkorn, obust methods in inverse theory, Inverse Problems, vol. 4,. 7 9, Oct [8] Y. Li and. Santosa, A comutational algorithm for minimizing total variation in image restoration, IEEE Transactions on Image Processing, vol. 5, , June 996. [9] W. L. Briggs, V. E. Henson, and S.. McCormick, A Multigrid Tutorial. SIAM Books, nd ed.,. [] P. odríguez and B. Wohlberg, Numerical methods for inverse roblems and adative decomosition NUMIPAD. Software library available from htt://numiad.sourceforge.net/.
Stereo Disparity Estimation in Moment Space
Stereo Disarity Estimation in oment Sace Angeline Pang Faculty of Information Technology, ultimedia University, 63 Cyberjaya, alaysia. angeline.ang@mmu.edu.my R. ukundan Deartment of Comuter Science, University
More informationAn Efficient VLSI Architecture for Adaptive Rank Order Filter for Image Noise Removal
International Journal of Information and Electronics Engineering, Vol. 1, No. 1, July 011 An Efficient VLSI Architecture for Adative Rank Order Filter for Image Noise Removal M. C Hanumantharaju, M. Ravishankar,
More informationImproved Image Super-Resolution by Support Vector Regression
Proceedings of International Joint Conference on Neural Networks, San Jose, California, USA, July 3 August 5, 0 Imroved Image Suer-Resolution by Suort Vector Regression Le An and Bir Bhanu Abstract Suort
More informationAdaptive Reciprocal Cell based Sparse Representation for Satellite Image Restoration
Vol.49 (SoftTech 04),.8-85 htt://dx.doi.org/0.457/astl.04.49.50 Adative Recirocal Cell based Sarse Reresentation for Satellite Image Restoration Yanfei He and Shunfeng Wang Deartment of Math and Statistic,
More informationA NON-NEGATIVE QUADRATIC PROGRAMMING APPROACH TO MINIMIZE THE GENERALIZED VECTOR-VALUED TOTAL VARIATION FUNCTIONAL
18th European Signal Processing Conference EUSIPCO-010 Aalborg, Denmark, August 3-7, 010 A NON-NEGATIVE QUADATIC POGAMMING APPOACH TO MINIMIZE THE GENEALIZED VECTO-VALUED TOTAL VAIATION UNCTIONAL Paul
More informationWavelet Based Statistical Adapted Local Binary Patterns for Recognizing Avatar Faces
Wavelet Based Statistical Adated Local Binary atterns for Recognizing Avatar Faces Abdallah A. Mohamed 1, 2 and Roman V. Yamolskiy 1 1 Comuter Engineering and Comuter Science, University of Louisville,
More informationA Novel Iris Segmentation Method for Hand-Held Capture Device
A Novel Iris Segmentation Method for Hand-Held Cature Device XiaoFu He and PengFei Shi Institute of Image Processing and Pattern Recognition, Shanghai Jiao Tong University, Shanghai 200030, China {xfhe,
More informationAn Image Inpainting Approach Based on the Poisson Equation
An mage nainting Aroach Based on the Poisson Equation Xiaowei Shao, Zhengkai Liu, Houqiang Li Deartment of Electronic Engineering and nformation Science Universit of Science and Technolog of China sw@mail.ustc.edu.cn,
More informationA CLASS OF STRUCTURED LDPC CODES WITH LARGE GIRTH
A CLASS OF STRUCTURED LDPC CODES WITH LARGE GIRTH Jin Lu, José M. F. Moura, and Urs Niesen Deartment of Electrical and Comuter Engineering Carnegie Mellon University, Pittsburgh, PA 15213 jinlu, moura@ece.cmu.edu
More informationInteractive Image Segmentation
Interactive Image Segmentation Fahim Mannan (260 266 294) Abstract This reort resents the roject work done based on Boykov and Jolly s interactive grah cuts based N-D image segmentation algorithm([1]).
More informationAn Efficient Coding Method for Coding Region-of-Interest Locations in AVS2
An Efficient Coding Method for Coding Region-of-Interest Locations in AVS2 Mingliang Chen 1, Weiyao Lin 1*, Xiaozhen Zheng 2 1 Deartment of Electronic Engineering, Shanghai Jiao Tong University, China
More informationFacial Expression Recognition using Image Processing and Neural Network
Keerti Keshav Kanchi / International Journal of Comuter Science & Engineering Technology (IJCSET) Facial Exression Recognition using Image Processing and eural etwork Keerti Keshav Kanchi PG Student, Deartment
More informationLecture 8: Orthogonal Range Searching
CPS234 Comutational Geometry Setember 22nd, 2005 Lecture 8: Orthogonal Range Searching Lecturer: Pankaj K. Agarwal Scribe: Mason F. Matthews 8.1 Range Searching The general roblem of range searching is
More informationEquality-Based Translation Validator for LLVM
Equality-Based Translation Validator for LLVM Michael Ste, Ross Tate, and Sorin Lerner University of California, San Diego {mste,rtate,lerner@cs.ucsd.edu Abstract. We udated our Peggy tool, reviously resented
More informationWeighted Page Rank Algorithm based on In-Out Weight of Webpages
Indian Journal of Science and Technology, Vol 8(34), DOI: 10.17485/ijst/2015/v8i34/86120, December 2015 ISSN (Print) : 0974-6846 ISSN (Online) : 0974-5645 eighted Page Rank Algorithm based on In-Out eight
More informationDistributed Estimation from Relative Measurements in Sensor Networks
Distributed Estimation from Relative Measurements in Sensor Networks #Prabir Barooah and João P. Hesanha Abstract We consider the roblem of estimating vectorvalued variables from noisy relative measurements.
More informationTotal Variation Denoising with Overlapping Group Sparsity
1 Total Variation Denoising with Overlapping Group Sparsity Ivan W. Selesnick and Po-Yu Chen Polytechnic Institute of New York University Brooklyn, New York selesi@poly.edu 2 Abstract This paper describes
More informationLearning Motion Patterns in Crowded Scenes Using Motion Flow Field
Learning Motion Patterns in Crowded Scenes Using Motion Flow Field Min Hu, Saad Ali and Mubarak Shah Comuter Vision Lab, University of Central Florida {mhu,sali,shah}@eecs.ucf.edu Abstract Learning tyical
More informationEE678 Application Presentation Content Based Image Retrieval Using Wavelets
EE678 Alication Presentation Content Based Image Retrieval Using Wavelets Grou Members: Megha Pandey megha@ee. iitb.ac.in 02d07006 Gaurav Boob gb@ee.iitb.ac.in 02d07008 Abstract: We focus here on an effective
More informationGradient Compared l p -LMS Algorithms for Sparse System Identification
Gradient Comared l -LMS Algorithms for Sarse System Identification Yong Feng,, Jiasong Wu, Rui Zeng, Limin Luo, Huazhong Shu. School of Biological Science and Medical Engineering, Southeast University,
More informationSimultaneous Tracking of Multiple Objects Using Fast Level Set-Like Algorithm
Simultaneous Tracking of Multile Objects Using Fast Level Set-Like Algorithm Martin Maška, Pavel Matula, and Michal Kozubek Centre for Biomedical Image Analysis, Faculty of Informatics Masaryk University,
More informationThe degree constrained k-cardinality minimum spanning tree problem: a lexisearch
Decision Science Letters 7 (2018) 301 310 Contents lists available at GrowingScience Decision Science Letters homeage: www.growingscience.com/dsl The degree constrained k-cardinality minimum sanning tree
More informationAUTOMATIC GENERATION OF HIGH THROUGHPUT ENERGY EFFICIENT STREAMING ARCHITECTURES FOR ARBITRARY FIXED PERMUTATIONS. Ren Chen and Viktor K.
inuts er clock cycle Streaming ermutation oututs er clock cycle AUTOMATIC GENERATION OF HIGH THROUGHPUT ENERGY EFFICIENT STREAMING ARCHITECTURES FOR ARBITRARY FIXED PERMUTATIONS Ren Chen and Viktor K.
More informationA Symmetric FHE Scheme Based on Linear Algebra
A Symmetric FHE Scheme Based on Linear Algebra Iti Sharma University College of Engineering, Comuter Science Deartment. itisharma.uce@gmail.com Abstract FHE is considered to be Holy Grail of cloud comuting.
More informationRandomized algorithms: Two examples and Yao s Minimax Principle
Randomized algorithms: Two examles and Yao s Minimax Princile Maximum Satisfiability Consider the roblem Maximum Satisfiability (MAX-SAT). Bring your knowledge u-to-date on the Satisfiability roblem. Maximum
More informationAN INTEGER LINEAR MODEL FOR GENERAL ARC ROUTING PROBLEMS
AN INTEGER LINEAR MODEL FOR GENERAL ARC ROUTING PROBLEMS Philie LACOMME, Christian PRINS, Wahiba RAMDANE-CHERIF Université de Technologie de Troyes, Laboratoire d Otimisation des Systèmes Industriels (LOSI)
More informationFace Recognition Using Legendre Moments
Face Recognition Using Legendre Moments Dr.S.Annadurai 1 A.Saradha Professor & Head of CSE & IT Research scholar in CSE Government College of Technology, Government College of Technology, Coimbatore, Tamilnadu,
More informationVector sigma filters for noise detection and removal in color images
J. Vis. Commun. Image R. 7 (26) 26 www.elsevier.com/locate/jvci Vector sigma filters for noise detection and removal in color images Rastislav Lukac a, *,, Bogdan Smolka b, Konstantinos N. Plataniotis
More informationImproved heuristics for the single machine scheduling problem with linear early and quadratic tardy penalties
Imroved heuristics for the single machine scheduling roblem with linear early and quadratic tardy enalties Jorge M. S. Valente* LIAAD INESC Porto LA, Faculdade de Economia, Universidade do Porto Postal
More informationarxiv: v1 [cs.mm] 18 Jan 2016
Lossless Intra Coding in with 3-ta Filters Saeed R. Alvar a, Fatih Kamisli a a Deartment of Electrical and Electronics Engineering, Middle East Technical University, Turkey arxiv:1601.04473v1 [cs.mm] 18
More informationLearning Robust Locality Preserving Projection via p-order Minimization
Proceedings of the Twenty-Ninth AAAI Conference on Artificial Intelligence Learning Robust Locality Preserving Projection via -Order Minimization Hua Wang, Feiing Nie, Heng Huang Deartment of Electrical
More informationSensitivity Analysis for an Optimal Routing Policy in an Ad Hoc Wireless Network
1 Sensitivity Analysis for an Otimal Routing Policy in an Ad Hoc Wireless Network Tara Javidi and Demosthenis Teneketzis Deartment of Electrical Engineering and Comuter Science University of Michigan Ann
More informationTexture Mapping with Vector Graphics: A Nested Mipmapping Solution
Texture Maing with Vector Grahics: A Nested Mimaing Solution Wei Zhang Yonggao Yang Song Xing Det. of Comuter Science Det. of Comuter Science Det. of Information Systems Prairie View A&M University Prairie
More informationEfficient Parallel Hierarchical Clustering
Efficient Parallel Hierarchical Clustering Manoranjan Dash 1,SimonaPetrutiu, and Peter Scheuermann 1 Deartment of Information Systems, School of Comuter Engineering, Nanyang Technological University, Singaore
More informationAUTOMATIC EXTRACTION OF BUILDING OUTLINE FROM HIGH RESOLUTION AERIAL IMAGERY
AUTOMATIC EXTRACTION OF BUILDING OUTLINE FROM HIGH RESOLUTION AERIAL IMAGERY Yandong Wang EagleView Technology Cor. 5 Methodist Hill Dr., Rochester, NY 1463, the United States yandong.wang@ictometry.com
More informationFace Recognition Based on Wavelet Transform and Adaptive Local Binary Pattern
Face Recognition Based on Wavelet Transform and Adative Local Binary Pattern Abdallah Mohamed 1,2, and Roman Yamolskiy 1 1 Comuter Engineering and Comuter Science, University of Louisville, Louisville,
More informationA BICRITERION STEINER TREE PROBLEM ON GRAPH. Mirko VUJO[EVI], Milan STANOJEVI] 1. INTRODUCTION
Yugoslav Journal of Oerations Research (00), umber, 5- A BICRITERIO STEIER TREE PROBLEM O GRAPH Mirko VUJO[EVI], Milan STAOJEVI] Laboratory for Oerational Research, Faculty of Organizational Sciences University
More informationApplying the fuzzy preference relation to the software selection
Proceedings of the 007 WSEAS International Conference on Comuter Engineering and Alications, Gold Coast, Australia, January 17-19, 007 83 Alying the fuzzy reference relation to the software selection TIEN-CHIN
More informationA Fast Image Restoration Method Based on an Improved Criminisi Algorithm
A Fast Image Restoration Method Based on an Imroved Algorithm Yue Chi1, Ning He2*, Qi Zhang1 Beijing Key Laboratory of Information Services Engineering, Beijing Union University, Beijing 100101, China.
More informationA New and Efficient Algorithm-Based Fault Tolerance Scheme for A Million Way Parallelism
A New and Efficient Algorithm-Based Fault Tolerance Scheme for A Million Way Parallelism Erlin Yao, Mingyu Chen, Rui Wang, Wenli Zhang, Guangming Tan Key Laboratory of Comuter System and Architecture Institute
More informationPatterned Wafer Segmentation
atterned Wafer Segmentation ierrick Bourgeat ab, Fabrice Meriaudeau b, Kenneth W. Tobin a, atrick Gorria b a Oak Ridge National Laboratory,.O.Box 2008, Oak Ridge, TN 37831-6011, USA b Le2i Laboratory Univ.of
More informationSemi-Supervised Learning Based Object Detection in Aerial Imagery
Semi-Suervised Learning Based Obect Detection in Aerial Imagery Jian Yao Zhongfei (Mark) Zhang Deartment of Comuter Science, State University of ew York at Binghamton, Y 13905, USA yao@binghamton.edu Zhongfei@cs.binghamton.edu
More informationSPARSE SIGNAL REPRESENTATION FOR COMPLEX-VALUED IMAGING Sadegh Samadi 1, M üjdat Çetin 2, Mohammad Ali Masnadi-Shirazi 1
SPARSE SIGNAL REPRESENTATION FOR COMPLEX-VALUED IMAGING Sadegh Samadi 1, M üjdat Çetin, Mohammad Ali Masnadi-Shirazi 1 1. Shiraz University, Shiraz, Iran,. Sabanci University, Istanbul, Turkey ssamadi@shirazu.ac.ir,
More informationAn Efficient Video Program Delivery algorithm in Tree Networks*
3rd International Symosium on Parallel Architectures, Algorithms and Programming An Efficient Video Program Delivery algorithm in Tree Networks* Fenghang Yin 1 Hong Shen 1,2,** 1 Deartment of Comuter Science,
More informationComplexity Issues on Designing Tridiagonal Solvers on 2-Dimensional Mesh Interconnection Networks
Journal of Comuting and Information Technology - CIT 8, 2000, 1, 1 12 1 Comlexity Issues on Designing Tridiagonal Solvers on 2-Dimensional Mesh Interconnection Networks Eunice E. Santos Deartment of Electrical
More informationA DEA-bases Approach for Multi-objective Design of Attribute Acceptance Sampling Plans
Available online at htt://ijdea.srbiau.ac.ir Int. J. Data Enveloment Analysis (ISSN 2345-458X) Vol.5, No.2, Year 2017 Article ID IJDEA-00422, 12 ages Research Article International Journal of Data Enveloment
More informationPREDICTING LINKS IN LARGE COAUTHORSHIP NETWORKS
PREDICTING LINKS IN LARGE COAUTHORSHIP NETWORKS Kevin Miller, Vivian Lin, and Rui Zhang Grou ID: 5 1. INTRODUCTION The roblem we are trying to solve is redicting future links or recovering missing links
More informationEfficient Processing of Top-k Dominating Queries on Multi-Dimensional Data
Efficient Processing of To-k Dominating Queries on Multi-Dimensional Data Man Lung Yiu Deartment of Comuter Science Aalborg University DK-922 Aalborg, Denmark mly@cs.aau.dk Nikos Mamoulis Deartment of
More informationExtracting Optimal Paths from Roadmaps for Motion Planning
Extracting Otimal Paths from Roadmas for Motion Planning Jinsuck Kim Roger A. Pearce Nancy M. Amato Deartment of Comuter Science Texas A&M University College Station, TX 843 jinsuckk,ra231,amato @cs.tamu.edu
More informationLeak Detection Modeling and Simulation for Oil Pipeline with Artificial Intelligence Method
ITB J. Eng. Sci. Vol. 39 B, No. 1, 007, 1-19 1 Leak Detection Modeling and Simulation for Oil Pieline with Artificial Intelligence Method Pudjo Sukarno 1, Kuntjoro Adji Sidarto, Amoranto Trisnobudi 3,
More informationConstrained Path Optimisation for Underground Mine Layout
Constrained Path Otimisation for Underground Mine Layout M. Brazil P.A. Grossman D.H. Lee J.H. Rubinstein D.A. Thomas N.C. Wormald Abstract The major infrastructure comonent reuired to develo an underground
More informationContinuous Visible k Nearest Neighbor Query on Moving Objects
Continuous Visible k Nearest Neighbor Query on Moving Objects Yaniu Wang a, Rui Zhang b, Chuanfei Xu a, Jianzhong Qi b, Yu Gu a, Ge Yu a, a Deartment of Comuter Software and Theory, Northeastern University,
More informationDirected File Transfer Scheduling
Directed File Transfer Scheduling Weizhen Mao Deartment of Comuter Science The College of William and Mary Williamsburg, Virginia 387-8795 wm@cs.wm.edu Abstract The file transfer scheduling roblem was
More informationSensitivity of multi-product two-stage economic lotsizing models and their dependency on change-over and product cost ratio s
Sensitivity two stage EOQ model 1 Sensitivity of multi-roduct two-stage economic lotsizing models and their deendency on change-over and roduct cost ratio s Frank Van den broecke, El-Houssaine Aghezzaf,
More informationGrouping of Patches in Progressive Radiosity
Grouing of Patches in Progressive Radiosity Arjan J.F. Kok * Abstract The radiosity method can be imroved by (adatively) grouing small neighboring atches into grous. Comutations normally done for searate
More informationAn accurate and fast point-to-plane registration technique
Pattern Recognition Letters 24 (23) 2967 2976 www.elsevier.com/locate/atrec An accurate and fast oint-to-lane registration technique Soon-Yong Park *, Murali Subbarao Deartment of Electrical and Comuter
More informationProbabilistic Visual Cryptography Schemes
The Comuter Journal Advance Access ublished December, 25 The Author 25. Published by Oxford University Press on behalf of The British Comuter Society. All rights reserved. For Permissions, lease email:
More informationExperimental Comparison of Shortest Path Approaches for Timetable Information
Exerimental Comarison of Shortest Path roaches for Timetable Information Evangelia Pyrga Frank Schulz Dorothea Wagner Christos Zaroliagis bstract We consider two aroaches that model timetable information
More informationIMS Network Deployment Cost Optimization Based on Flow-Based Traffic Model
IMS Network Deloyment Cost Otimization Based on Flow-Based Traffic Model Jie Xiao, Changcheng Huang and James Yan Deartment of Systems and Comuter Engineering, Carleton University, Ottawa, Canada {jiexiao,
More informationA Generalized Iterated Shrinkage Algorithm for Non-convex Sparse Coding
A Generalized Iterated Shrinkage Algorithm for Non-convex Sarse Coding Wangmeng Zuo,3, Deyu Meng 2, Lei Zhang 3, Xiangchu Feng 4, David Zhang 3 Harbin Institute of Technology 2 Xi an Jiaotong University
More informationVehicle Logo Recognition Using Modest AdaBoost and Radial Tchebichef Moments
Proceedings of 0 4th International Conference on Machine Learning and Comuting IPCSIT vol. 5 (0) (0) IACSIT Press, Singaore Vehicle Logo Recognition Using Modest AdaBoost and Radial Tchebichef Moments
More informationImplementation of Evolvable Fuzzy Hardware for Packet Scheduling Through Online Context Switching
Imlementation of Evolvable Fuzzy Hardware for Packet Scheduling Through Online Context Switching Ju Hui Li, eng Hiot Lim and Qi Cao School of EEE, Block S Nanyang Technological University Singaore 639798
More informationRegularized Kernel Regression for Image Deblurring
Regularized Kernel Regression for Image Deblurring Hiroyuki Takeda, Sina Farsiu, and Peyman Milanfar Department of Electrical Engineering, niversity of California at Santa Cruz {htakeda,farsiu,milanfar}@soeucscedu
More informationImage Restoration using Accelerated Proximal Gradient method
Image Restoration using Accelerated Proximal Gradient method Alluri.Samuyl Department of computer science, KIET Engineering College. D.Srinivas Asso.prof, Department of computer science, KIET Engineering
More informationA Parallel Algorithm for Constructing Obstacle-Avoiding Rectilinear Steiner Minimal Trees on Multi-Core Systems
A Parallel Algorithm for Constructing Obstacle-Avoiding Rectilinear Steiner Minimal Trees on Multi-Core Systems Cheng-Yuan Chang and I-Lun Tseng Deartment of Comuter Science and Engineering Yuan Ze University,
More informationCMSC 425: Lecture 16 Motion Planning: Basic Concepts
: Lecture 16 Motion lanning: Basic Concets eading: Today s material comes from various sources, including AI Game rogramming Wisdom 2 by S. abin and lanning Algorithms by S. M. LaValle (Chats. 4 and 5).
More informationSpaceTwist: Managing the Trade-Offs Among Location Privacy, Query Performance, and Query Accuracy in Mobile Services
SaceTwist: Managing the Trade-Offs Among Location Privacy, Query Performance, and Query Accuracy in Mobile Services Man Lung Yiu, Christian S. Jensen, Xuegang Huang, Hua Lu Deartment of Comuter Science,
More informationHigh Order Accurate Curved-element Implementation in OpenFOAM for Discontinuous Galerkin Method
,.47-52 htt://dx.doi.org/10.14257/astl.2016.142.08 High Order Accurate Curved-element Imlementation in OenFOAM for Discontinuous Galerkin Method Li-Yang Xu 1, Shuai Ye 1, Yong-Quan Feng 1, Hao Li 1, Xin-Hai
More informationBayesian Oil Spill Segmentation of SAR Images via Graph Cuts 1
Bayesian Oil Sill Segmentation of SAR Images via Grah Cuts 1 Sónia Pelizzari and José M. Bioucas-Dias Instituto de Telecomunicações, I.S.T., TULisbon,Lisboa, Portugal sonia@lx.it.t, bioucas@lx.it.t Abstract.
More informationSELF-ADJUSTED ACTIVE CONTOURS USING MULTI-DIRECTIONAL TEXTURE CUES. E.A. Mylona, M.A. Savelonas, D. Maroulis
SELF-ADJUSTED ACTIVE CONTOURS USING MULTI-DIRECTIONAL TEXTURE CUES E.A. Mylona, M.A. Savelonas, D. Maroulis Deartment of Informatics and Telecommunications, University of Athens, Greece {emylona, msavel,
More informationThe Research on Curling Track Empty Value Fill Algorithm Based on Similar Forecast
Research Journal of Alied Sciences, Engineering and Technology 6(8): 1472-1478, 2013 ISSN: 2040-7459; e-issn: 2040-7467 Maxwell Scientific Organization, 2013 Submitted: October 31, 2012 Acceted: January
More informationIdentity-sensitive Points-to Analysis for the Dynamic Behavior of JavaScript Objects
Identity-sensitive Points-to Analysis for the Dynamic Behavior of JavaScrit Objects Shiyi Wei and Barbara G. Ryder Deartment of Comuter Science, Virginia Tech, Blacksburg, VA, USA. {wei,ryder}@cs.vt.edu
More informationOMNI: An Efficient Overlay Multicast. Infrastructure for Real-time Applications
OMNI: An Efficient Overlay Multicast Infrastructure for Real-time Alications Suman Banerjee, Christoher Kommareddy, Koushik Kar, Bobby Bhattacharjee, Samir Khuller Abstract We consider an overlay architecture
More information3D Surface Simplification Based on Extended Shape Operator
3D Surface Simlification Based on Extended Shae Oerator JUI-LIG TSEG, YU-HSUA LI Deartment of Comuter Science and Information Engineering, Deartment and Institute of Electrical Engineering Minghsin University
More informationAPPLICATION OF PARTICLE FILTERS TO MAP-MATCHING ALGORITHM
APPLICATION OF PARTICLE FILTERS TO MAP-MATCHING ALGORITHM Pavel Davidson 1, Jussi Collin 2, and Jarmo Taala 3 Deartment of Comuter Systems, Tamere University of Technology, Finland e-mail: avel.davidson@tut.fi
More informationRandomized Selection on the Hypercube 1
Randomized Selection on the Hyercube 1 Sanguthevar Rajasekaran Det. of Com. and Info. Science and Engg. University of Florida Gainesville, FL 32611 ABSTRACT In this aer we resent randomized algorithms
More informationSubmission. Verifying Properties Using Sequential ATPG
Verifying Proerties Using Sequential ATPG Jacob A. Abraham and Vivekananda M. Vedula Comuter Engineering Research Center The University of Texas at Austin Austin, TX 78712 jaa, vivek @cerc.utexas.edu Daniel
More informationA Petri net-based Approach to QoS-aware Configuration for Web Services
A Petri net-based Aroach to QoS-aware Configuration for Web s PengCheng Xiong, YuShun Fan and MengChu Zhou, Fellow, IEEE Abstract With the develoment of enterrise-wide and cross-enterrise alication integration
More informationFuzzy Logic and ANFIS based Image Centroid Tracking and Filtering Schemes
Fuzzy Logic and ANFIS based Image Centroid Tracking and Filtering Schemes Parimala, A., Asstt. Prof., Det. of Telecommunication, MSRIT; Jain University Research Scholar, Bangalore. Raol, J. R., Prof. Emeritus,
More informationDr Pavan Chakraborty IIIT-Allahabad
GVC-43 Lecture - 5 Ref: Donald Hearn & M. Pauline Baker, Comuter Grahics Foley, van Dam, Feiner & Hughes, Comuter Grahics Princiles & Practice Dr Pavan Chakraborty IIIT-Allahabad Summary of line drawing
More information2-D Fir Filter Design And Its Applications In Removing Impulse Noise In Digital Image
-D Fir Filter Design And Its Alications In Removing Imulse Noise In Digital Image Nguyen Thi Huyen Linh 1, Luong Ngoc Minh, Tran Dinh Dung 3 1 Faculty of Electronic and Electrical Engineering, Hung Yen
More informationPattern Recognition Letters
Pattern Recognition Letters 30 (2009) 114 122 Contents lists available at ScienceDirect Pattern Recognition Letters journal homeage: www.elsevier.com/locate/atrec A stroke filter and its alication to text
More informationMultiscale Flat Norm Signatures for Shapes and Images
Applied Mathematical Sciences, Vol. 4, 21, no. 14, 667-68 Multiscale Flat Norm Signatures for Shapes and Images Kevin R. Vixie, Keith Clawson, Thomas J. Asaki 1 Mathematics Department, Washington State
More informationPower Savings in Embedded Processors through Decode Filter Cache
Power Savings in Embedded Processors through Decode Filter Cache Weiyu Tang Rajesh Guta Alexandru Nicolau Deartment of Information and Comuter Science University of California, Irvine Irvine, CA 92697-3425
More informationAuto-Tuning Distributed-Memory 3-Dimensional Fast Fourier Transforms on the Cray XT4
Auto-Tuning Distributed-Memory 3-Dimensional Fast Fourier Transforms on the Cray XT4 M. Gajbe a A. Canning, b L-W. Wang, b J. Shalf, b H. Wasserman, b and R. Vuduc, a a Georgia Institute of Technology,
More informationIntroduction to Parallel Algorithms
CS 1762 Fall, 2011 1 Introduction to Parallel Algorithms Introduction to Parallel Algorithms ECE 1762 Algorithms and Data Structures Fall Semester, 2011 1 Preliminaries Since the early 1990s, there has
More informationGRAPH CUT OPTIMIZATION FOR THE MUMFORD-SHAH MODEL
GRAPH CUT OPTIMIZATION FOR THE MUMFORD-SHAH MODEL Noha El Zehiry 1,2, Steve Xu 2, Prasanna Sahoo 2 and Adel Elmaghraby 1 1 Deartment of Comuter Engineering and Comuter Science, University of Louisville,
More informationarxiv: v1 [cs.dc] 13 Nov 2018
Task Grah Transformations for Latency Tolerance arxiv:1811.05077v1 [cs.dc] 13 Nov 2018 Victor Eijkhout November 14, 2018 Abstract The Integrative Model for Parallelism (IMP) derives a task grah from a
More informationAccelerating the Hessian-free Gauss-Newton Full-waveform Inversion via Preconditioned Conjugate Gradient Method
Accelerating the Hessian-free Gauss-Newton Full-waveform Inversion via Preconditioned Conjugate Gradient Method Wenyong Pan 1, Kris Innanen 1 and Wenyuan Liao 2 1. CREWES Project, Department of Geoscience,
More informationLecture 2: Fixed-Radius Near Neighbors and Geometric Basics
structure arises in many alications of geometry. The dual structure, called a Delaunay triangulation also has many interesting roerties. Figure 3: Voronoi diagram and Delaunay triangulation. Search: Geometric
More informationTheoretical Analysis of Graphcut Textures
Theoretical Analysis o Grahcut Textures Xuejie Qin Yee-Hong Yang {xu yang}@cs.ualberta.ca Deartment o omuting Science University o Alberta Abstract Since the aer was ublished in SIGGRAPH 2003 the grahcut
More informationGaussian and Laplacian of Gaussian weighting functions for robust feature based tracking q
Pattern Recognition Letters 26 (25) 1995 25 www.elsevier.com/locate/atrec Gaussian and Lalacian of Gaussian weighting functions for robust feature based tracking q Meghna Singh a, Mrinal K. Mandal a, Anu
More informationon continuous k-means monitoring over moving objects,
Continuous k-means Monitoring over Moving Objects Zhenjie Zhang, Yin Yang, Anthony K.H. Tung, and Dimitris Paadias Abstract Given a dataset P, a k-means query returns k oints in sace (called centers),
More informationNON-UNIFORM TO UNIFORM IMAGE RESAMPLING UTILIZING A 2D FARROW STRUCTURE
NON-UNIFORM TO UNIFORM IMAGE RESAMPLING UTILIZING A D FARROW STRUCTURE Harish Essay Sanaran, Mihail Georgiev, Atanas Gotchev and Karen Egiaarian Noia Cororation, Hermianatu, Tamere, FINLAND, harish.essaysanaran@noia.com
More informationCollective communication: theory, practice, and experience
CONCURRENCY AND COMPUTATION: PRACTICE AND EXPERIENCE Concurrency Comutat.: Pract. Exer. 2007; 19:1749 1783 Published online 5 July 2007 in Wiley InterScience (www.interscience.wiley.com)..1206 Collective
More informationMatlab Virtual Reality Simulations for optimizations and rapid prototyping of flexible lines systems
Matlab Virtual Reality Simulations for otimizations and raid rototying of flexible lines systems VAMVU PETRE, BARBU CAMELIA, POP MARIA Deartment of Automation, Comuters, Electrical Engineering and Energetics
More informationLecture 18. Today, we will discuss developing algorithms for a basic model for parallel computing the Parallel Random Access Machine (PRAM) model.
U.C. Berkeley CS273: Parallel and Distributed Theory Lecture 18 Professor Satish Rao Lecturer: Satish Rao Last revised Scribe so far: Satish Rao (following revious lecture notes quite closely. Lecture
More informationBrief Contributions. A Geometric Theorem for Network Design 1 INTRODUCTION
IEEE TRANSACTIONS ON COMPUTERS, VOL. 53, NO., APRIL 00 83 Brief Contributions A Geometric Theorem for Network Design Massimo Franceschetti, Member, IEEE, Matthew Cook, and Jehoshua Bruck, Fellow, IEEE
More informationSwift Template Matching Based on Equivalent Histogram
Swift emlate Matching ased on Equivalent istogram Wangsheng Yu, Xiaohua ian, Zhiqiang ou * elecommunications Engineering Institute Air Force Engineering University Xi an, PR China *corresonding author:
More informationA Study of Protocols for Low-Latency Video Transport over the Internet
A Study of Protocols for Low-Latency Video Transort over the Internet Ciro A. Noronha, Ph.D. Cobalt Digital Santa Clara, CA ciro.noronha@cobaltdigital.com Juliana W. Noronha University of California, Davis
More information