A New Framework of Quantifying Differences Between Images by Matching Gradient Fields and Its Application to Image Blending
|
|
- Samson Leslie Barnett
- 5 years ago
- Views:
Transcription
1 A New Framework o Quantiying Dierences Between Images by Matching Gradient Fields and Its Application to Image Blending Wei-Hsun Liao Chin-Lung Yu Marvin Bergsneider Luminita Vese and Sung-Cheng Huang Abstract--This paper presents a new method or quantiying the dierences between images The proposed model is based on matching the gradient ields o two images We irst deine new image spaces in which images are considered equivalent under a similarity group actions and the dierence between two image classes is then deined by employing the Cauchy-Schwarz inequality to the gradient ields The advantage o our approach is that images are identiied by their relative contrasts and thus is scale ree Using this approach we are able to achieve image blending in a novel way By modiying the group actions we extend our basic model to more general equivalence classes The variational problems and the corresponding Euler-Lagrange equations associated to these models are proposed and the gradient descent time dependent partial dierential equations are derived Fast and eicient solvers employing the Additive Operator Splitting scheme are also presented We tested our models on simulation images as well as real brain MRI and PET images rom normal control subjects I INTRODUCTION HE mathematical ormulation and theories o objects Tand shapes have been studied extensively in computer vision statistics geology biomedical science and image processing literature A lot o eorts have been made to mathematically quantiy the dierence and to deine a distance between shapes or objects In this paper we turn to the problem o deining dierences between images By doing so we could achieve image blending in a novel way by minimizing the distance between two images Manuscript received November 8 00 This work was supported in part by NIH Grant No NS30308 and DoE Grant No DE-FC0387ER60615 Wei-Hsun Liao ( euillet@uclaedu) and Sung-Cheng Huang ( hhuang@mednetuclaedu) are with the Department o Biomathematics and Department o Molecular and Medical Pharmacology at UCLA Los Angeles CA USA Chin-Lung Yu is with the Department o Molecular and Medical Pharmacology at UCLA Los Angeles CA USA ( clyu@uclaedu) Luminita Vese is with the Department o Mathematics at UCLA Los Angeles CA USA ( lvese@mathuclaedu) Marvin Bergsneider is with the Brain Injury Research Center o David Geen School o Medicine at UCLA Los Angeles CA USA ( mbergsneider@mednetuclaedu) Image blending [1] and image morphing have wide applications in computer vision image processing and biomedical imaging For a more complete treatment o this subject we reer the reader to [ 3] Image morphing usually consists mainly o two steps The irst one is warping which involves a coordinate transormation to align user-deined landmarks The next step is to apply image blending which is usually done by linear interpolation (or cross dissolving) Linear interpolation is based on the sum o least squares o two images and g: 1 dist( g) ( g) dx Minimizing this distance with respect to simply gives us a linear interpolation between images and g Although it is the most commonly used it is not necessarily the best First it is not scale invariant i we rescale one o the images then we obtain a dierent distance Moreover two images that only dier in scaling are considered dierent under this distance measure In this paper we try to interpret images in a dierent way and avoid these disadvantages In order to accomplish this we irst redeine the meaning o identical or equivalent images We turn to the idea o shape spaces to look or solutions II RELATED WORK AND MOTIVATIONS In this section we closely ollow the discussion in [4] The study o shapes has been proposed in [5-9] as a novel environment in statistics or comparing shapes where a metric and a probabilistic measure could be deined It has been successully used in comparing shapes where distinct landmarks could be ound In this case shapes are deined as the equivalence classes o M points in R N where the equivalence class is deined by a similarity group including rotation translation and scaling The space o all shapes could thus be viewed as a collection o iber bundles where motions in the space could be along the iber or across the iber Motions along the iber are simply motions between equivalent classes while motions across the iber correspond to deormations rom one class to another class Another approach that also has motivated our work is the success o variational methods and the use o PDEs to model /03/$ IEEE 109
2 shapes and images For example in [10-13] a general rame work o deormations and warping was proposed and the application o this ramework to object comparison has shown to be very successul In this ramework deormations could be rigorously deined in terms o group actions and the notion o distance could be made precise using these theories For contributions rom other authors see or example [14-16] Our work is also closely related to nonlinear image smoothing techniques Anisotropic diusion [17 18] is one o the irst nonlinear smoothing techniques The original idea o anisotropic diusion is to smooth selectively within more homogeneous areas while preserving edges when a large gradient is detected The PDE that governs anisotropic diusion can be written in the ollowing orm: I t div( c( I ) I) where I is an image and c is the diusion coeicient It has been noted in [19] that this PDE could also be derived by minimizing the ollowing energy: min E( I) p( I ) dx I Ω The unction p is the penalty unction The PDE that minimizes the energy has the orm: I t div ( p' ( I ) I I ) Thus the diusion coeicient c in the original model can be written as a unction o p: c( ξ) p'( ξ) ξ In developing our models we were originally looking or ways to extend the anisotropic diusion which could be viewed as a prior-less smoothing as opposed to a priorbased smoothing The connection between anisotropic diusion and our work resides in the act that by orming a suitable variational problem our work could be incorporated into a prior-based smoothing when two images being compared contain similar inormation with one (the template) o the two images having a better image quality For more details please reer to the results and discussion sections III DESCRIPTION OF THE MODEL A A natural way o deining equivalent images Inspired by the idea o shape spaces we want to deine dierences between images in a way that the dierence is deined in terms o equivalence classes o images Dierence between two images will only be nonzero when they are not equivalent But what similarity group action should be chosen to deine equivalent classes? Human perception o an image depends mostly by the contrast o the image Areas with high contrast are perceived as boundaries o objects while areas with more homogeneous intensities are oten viewed as inside an object So our goal is to construct a new distance measure that only takes into account the relative change in intensity o the images This means we should consider two images I 1 and I to be equivalent i there exist real numbers k 1 > 0 and k such that the ollowing holds: I ( x) k I ( x) + k x Ω (1) 1 1 So our similarity group should include translation (k ) and rescaling (k 1 ) o the intensities Later on we will relax the restrictions on k 1 by allowing negativity o k 1 This concept shares many eatures with the shape spaces For example images like shape spaces are organized in iber bundles where a iber bundle is just the collection o images that are equivalent under the similarity group actions we just deined Also motions rom one image to another also all into two categories: across the bundles and within the bundles by deining them in a similar way Thus we have a new space on which images must be identiied by their gradient ields regardless o the scaling Now we need a proper distance measure which is always nonnegative and is zero only between equivalent images We will derive distance measure between images and g based on the ollowing inequality which is derived rom the Cauchy- Schwarz inequality dx g d x ( g)dx 0 B The models We will present three dierent distance measures based on the above inequality The irst model is related to measuring dierences o images where equivalent images are deined as in (1) without allowing negative k 1 while the second model allows any real value o k 1 The last one is a urther generalization Each one o them has slightly dierent properties and leads to its own PDE by employing gradient descent to solve the Euler-Lagrange equations For images and g we can look at the ollowing minimization problems to match the gradient ield o to the gradient ield o g: model a min E( ) dx g d x ( g)dx model b min E( ) d x g d x ( gd x) Model B is dierent rom model a in that the negative image o g (up to a scalar) is also considered to be equivalent to g model c We propose a third model that is even more nonlinear by applying the absolute unction to the integrand o the last term on the right hand side o model b It also turns out to be the most interesting o the three models (see results) By equivalence class in this space we mean that the ollowing energy between two images is 0: /03/$ IEEE 1093
3 min E( ) d x g d x ( g d x) C Modiied models An undesirable property o these models is that since they are not properly normalized images are collapsed to constant images when time goes to ininity In order to avoid this problem we could normalize the energy to be minimized so that it is always between 0 and 1 Thus we look at the ollowing modiied versions o our models: model a g dx min E( ) 1 1/ ( d x g d x ) g dx g t d x g model b min E( ) 1 ( g d x) d d x g x ( g d x) ( d ) g x g t d x g model c g dx min E( ) 1 ( d x g d x ) 1/ g dx div ( sgn( g) g) t d x g IV IMPLEMENTATION Since the PDE we need to solve is an inhomogeneous heat equation let us write down the PDE in the ollowing orm: F + α t The semi-implicit scheme is used to solve this PDE which now reads: r r r d 1 n + 1 n n n I tα A tf l + l 1 ( ) Furthermore by using Additive Operator Splitting (AOS) scheme [0 1] we update in the ollowing way: r 1 1 d n n 1 r r + n n ( I d tα Al ) ( + tf ) d l 1 Here d is the dimension o the problem and A is the onedimensional discretized Laplacian operator along lth axis The splitting is o order one in time and order two in space Moreover it is unconditionally stable By splitting the operator into a coordinate-by-coordinate ashion we now only need to invert a tri-diagonal matrix along each coordinate and this allows an O(m ) implementation by the Thomas algorithm V RESULTS For results presented in this section all calculation were perormed on 18 by 18 grid points and the spatial step is 1 A Test o the basic model We only present the numerical results or model a The goal is to blend a circular shape (Fig 1(a)) to an oval shape The intermediate images during blending process are shown in Fig 1(b)-(d) We notice that along the blending process the image intensities outside the circles also change with time which does not happen in simple interpolation technique This phenomenon generates a low-like motion that gives novel visual eects These eects might be closely related the resemblance o the governing PDE to the heat equation B Test o modiied model on medical images Simple modiication o model c allows us to combine eatures o two images and to generate a hybrid or usion image The idea is similar to [18 ] We consider denoising an image 0 by solving in the ollowing problem: 1 min E( ) ( 0 ) dx 1/ + c 1 g d x ( d x g d x ) The PDE that governs this minimization problem is: 1 ( 0 ) + c( dx g d x) t () g dx div[ sgn ( g) g] The above variational problem could be viewed as a priorbased image denoising The idea is that 0 is a corrupted image that we would like to recover Furthermore we know that g contains the same inormation but we do not know exactly how the intensity values relate in these two images By applying the above model we seek to ind a recovered image that is close to g while constraining the result to be not ar away rom the initial image /03/$ IEEE 1094
4 Figure 1 Image blending process described in A o the results section (a): The initial image containing an circular object; (b) (c) and (d): Intermediate images The template image is visually similar to (d) Figure Image blending in brain MRI imaging (a): SPGR MRI image as the initial image; (b): FDG-PET image as the template; (c): hybrid o image (a) and (b); (d): T-star MRI image as another template; (e): hybrid image o (a) and (d) In order to test this modiied model we turn to biomedical imaging or test images In biomedicine images o dierent modalities such as CT MRI and PET are oten obtained or the same subject This is especially true in brain imaging as dierent images could provide inormation on brain structures and unctions We chose to test our model on MRI and PET images FDG-PET imaging and MRI imaging o the brain were obtained rom a normal subject in the Brain Trauma Project at UCLA The available sequences in MRI images include SPGR T and T-star All MRI images were irst co-registered to PET using mutual inormation [3 4] as the similarity measure The mutual inormation between MRI and PET images were optimized using Powell s multi-dimensional search algorithm We looked at two speciic issues First we combined the eatures in MRI and PET images by solving () Second we examined how dierent sequences in MRI images inluenced each other under this model For all computations the intensity values o all images were normalized to 0 to 56 beore computation To investigate the irst issue we chose SPGR MRI images as the test images since SPGR images had the best spatial resolution and image quality Ater MRI-PET coregistration was perormed an SPGR image slice (Fig (a)) and a PET image slice (Fig (b)) were extracted We then used the SPGR image as the initial image and the corresponding PET image as the template The usion image was calculated by solving () The weight c o the regularizer is determined empirically The 6 inal c used in this test was 10 The corresponding PDE was solved at t to get the inal usion image The result is shown in Fig (c) Note that the overall contrast o the usion images is very much like that o the initial SPGR image while the shape o the gray and white matters resemble that o the corresponding PET image For the second issue we chose T-star and SPGR images T-star images usually have lower resolution and in this subject the extracted slice was noted to have magnetic ield inhomogeniety that resulted in shading o the image (rontal lobe has a much higher value than the occipital lobe) as shown in Fig (d) On the other hand the corresponding SPGR (Fig (a)) was o good quality without shading We thus used the SPGR image as the template and the T-star as the initial image and see i we could improve the quality o the T-star image The result is shown in Fig (e) The weight c used in 6 this test was 6 10 and again the corresponding PDE was solved at t We notice that although the resolution o the usion image did not improve signiicantly the intensity inhomogeneity almost completely disappeared VI DISCUSSIONS AND CONCLUSIONS In this paper we propose a new method o deining the dierence o two images The idea is to look at the gradient ield instead o the image itsel This results in a novel approach o image blending and a method o combining or hybridizing two images The advantage o our approach is that it is scale ree; the scale o the image is determined by the PDE that governs the minimization process The drawback o our approach is that since we are working on the gradient ields we have to restrict ourselves to the space H 1 which in theory does not allow discontinuous solutions Though we do point out that numerically images are discretized signals and the gradient could always be viewed as inite Experimental results showed that our approach generated low-like transition between images This low-like blending phenomenon is not observed in simple linear interpolation blending and provides a novel perception to human eyes Although our model generates better results than simple interpolation method it still has the undesirable overlapping eect o objects in the transition process Thus we have to emphasize that our model is not intended to improve or replace /03/$ IEEE 1095
5 image warping but instead oers a dierent view o image blending and o quantiying dierences between images We also applied our approach to biomedical images We generated a hybrid MRI-PET image with eatures rom both modalities We also showed that by incorporating the inormation rom SPGR image we were able to correct the magnetic ield inhomogeneity o the corresponding T-star images Another issue is that in model c we seek to minimize the ollowing energy: g dx min E( ) 1 1/ ( d x g d x ) We can urther generalize the model by replacing the absolute value unction in the numerator with a continuous unction Ψ : Ψ( g)dx min E( ) 1 1/ ( d x g d x ) I we urther require that this unction is always non-negative and symmetric and its graph lies below the absolute value unction then our energy is still deined and between 0 and 1 Moreover i Ψ is dierentiable then the associated gradient descent PDE becomes: Ψ ( g)dx div ( Ψ'( g) g) t d x g As a concrete example the ollowing choice o Ψ satisies all the requirements: Ψ ( ξ ) ε + log(1 + ξ ) where ε is a small number By modiying the model this way one can achieve an anisotropic orm o the original model Work to explore how dierent choices o Ψ will inluence the results is still ongoing in our laboratory VII ACKNOWLEDGMENTS This study was partially supported by NIH Grant NS30308 and DoE Grant DE-FC0387-ER60615 We would like to thank Dr Richard Tsai or the valuable discussions [6] K V Mardia and I L Dryden "Shape Distributions or Landmark Data" Advances in Applied Probability vol 1 pp [7] D G Kendall "Shape Maniolds Procrustean Metrics and Complex Projective Spaces" Bulletin o the London Mathematical Society vol 16 pp [8] D G Kendall Shape and shape theory New York: Wiley 1999 [9] H L Le and D G Kendall "The Riemannian Structure o Euclidean Shape Spaces - a Novel Environment or Statistics" Annals o Statistics vol 1 pp [10] L Younes "Computable elastic distances between shapes" Siam Journal on Applied Mathematics vol 58 pp [11] L Younes "A distance or elastic matching in object recognition" Comptes Rendus De L Academie Des Sciences Serie I-Mathematique vol 3 pp [1] L Younes "Optimal matching between shapes via elastic deormations" Image and Vision Computing vol 17 pp [13] M I Miller and L Younes "Group actions homeomorphisms and matching: A general ramework" International Journal o Computer Vision vol 41 pp [14] B B Kimia A R Tannenbaum and S W Zucker "Shapes Shocks and Deormations 1 The Components o -Dimensional Shape and the Reaction-Diusion Space" International Journal o Computer Vision vol 15 pp [15] B B Kimia A Tannenbaum and S W Zucker "On the Evolution o Curves Via a Function o Curvature 1 The Classical Case" Journal o Mathematical Analysis and Applications vol 163 pp [16] B B Kimia A Tannenbaum and S W Zucker "Toward a Computational Theory o Shape - an Overview" Lecture Notes in Computer Science vol 47 pp [17] P Perona and J Malik "Scale-space and edge detection using anisotropic diusion" IEEE Transactions on Pattern Analysis and Machine Intelligence vol 1 pp [18] L I Rudin S Osher and E Fatemi "Nonlinear total variation based noise removal algorithms" 199 [19] Y Yu-Li X Wenyuan A Tannenbaum and M Kaveh "Behavioral analysis o anisotropic diusion in image processing" IEEE Transactions on Image Processing vol 5 pp [0] J Weickert "Recursive separable schemes or nonlinear diusion ilters" 1997 [1] J Weickert B M T H Romeny and M A Viergever "Eicient and reliable schemes or nonlinear diusion iltering" IEEE Transactions on Image Processing vol 7 pp [] G Aubert and L Vese "A variational method in image recovery" SIAM Journal on Numerical Analysis vol 34 pp [3] P Viola and W M Wells III "Alignment by maximization o mutual inormation" 1995 [4] F Maes A Collignon D Vandermeulen G Marchal and P Suetens "Multimodality image registration by maximization o mutual inormation" IEEE Transactions on Medical Imaging vol 16 pp VIII REFERENCES [1] R T Whitaker "A level-set approach to image blending" IEEE Transactions on Image Processing vol 9 pp [] C Watkins A Sadun and S Marenka Modern image processing: warping morphing and classical techniques Boston: Academic Press 1993 [3] G Wolberg Digital image warping Los Alamitos Cali: IEEE Computer Society Press 199 [4] S Soatto and A J Yezzi "Deormotion: Deorming Motion Shape Average and the Joint Registration and Segmentation o Images" UCLA Los Angeles UCLA CAM report 01-3 December 001 [5] T K Carne "The Geometry o Shape Spaces" Proceedings o the London Mathematical Society vol 61 pp /03/$ IEEE 1096
Neighbourhood 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 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 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 informationMAPI 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 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 informationGeometric Registration for Deformable Shapes 2.2 Deformable Registration
Geometric Registration or Deormable Shapes 2.2 Deormable Registration Variational Model Deormable ICP Variational Model What is deormable shape matching? Example? What are the Correspondences? Eurographics
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 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 informationBrain Warping Via Landmark Points and Curves with a Level Set Representation
Brain Warping Via Landmark Points and Curves with a Level Set Representation Andrew Y. Wang, Alex D. Leow, 2 Hillary D. Protas, Arthur W. Toga, Paul M. Thompson UCLA Laboratory of Neuro Imaging, Los Angeles,
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 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 informationNIH Public Access Author Manuscript Proc Soc Photo Opt Instrum Eng. Author manuscript; available in PMC 2014 October 07.
NIH Public Access Author Manuscript Published in final edited form as: Proc Soc Photo Opt Instrum Eng. 2014 March 21; 9034: 903442. doi:10.1117/12.2042915. MRI Brain Tumor Segmentation and Necrosis Detection
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 informationDenoising using Curvelet Transform with Improved Thresholds for Application on Ultrasound Images
7 Denoising using Curvelet Transorm with Improved Thresholds or Application on Ultrasound Images Er. Tainder Kaur, Assistant Proessor, SBBSIET Padhiana, alandhar Er. Harpreet Kaur, Assistant Proessor,
More informationYunyun Yang, Chunming Li, Chiu-Yen Kao and Stanley Osher. Speaker: Chiu-Yen Kao (Math Department, The Ohio State University) BIRS, Banff, Canada
Yunyun Yang, Chunming Li, Chiu-Yen Kao and Stanley Osher Speaker: Chiu-Yen Kao (Math Department, The Ohio State University) BIRS, Banff, Canada Outline Review of Region-based Active Contour Models Mumford
More informationImage denoising using TV-Stokes equation with an orientation-matching minimization
Image denoising using TV-Stokes equation with an orientation-matching minimization Xue-Cheng Tai 1,2, Sofia Borok 1, and Jooyoung Hahn 1 1 Division of Mathematical Sciences, School of Physical Mathematical
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 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 informationPiecewise polynomial interpolation
Chapter 2 Piecewise polynomial interpolation In ection.6., and in Lab, we learned that it is not a good idea to interpolate unctions by a highorder polynomials at equally spaced points. However, it transpires
More informationA SAR IMAGE REGISTRATION METHOD BASED ON SIFT ALGORITHM
A SAR IMAGE REGISTRATION METHOD BASED ON SIFT ALGORITHM W. Lu a,b, X. Yue b,c, Y. Zhao b,c, C. Han b,c, * a College o Resources and Environment, University o Chinese Academy o Sciences, Beijing, 100149,
More informationSTATISTICS AND ANALYSIS OF SHAPE
Control and Cybernetics vol. 36 (2007) No. 2 Book review: STATISTICS AND ANALYSIS OF SHAPE by H. Krim, A. Yezzi, Jr., eds. There are numerous definitions of a notion of shape of an object. These definitions
More informationAutomated Segmentation Using a Fast Implementation of the Chan-Vese Models
Automated Segmentation Using a Fast Implementation of the Chan-Vese Models Huan Xu, and Xiao-Feng Wang,,3 Intelligent Computation Lab, Hefei Institute of Intelligent Machines, Chinese Academy of Science,
More information3D Surface Reconstruction of the Brain based on Level Set Method
3D Surface Reconstruction of the Brain based on Level Set Method Shijun Tang, Bill P. Buckles, and Kamesh Namuduri Department of Computer Science & Engineering Department of Electrical Engineering University
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 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 informationStudy and Analysis of Edge Detection and Implementation of Fuzzy Set. Theory Based Edge Detection Technique in Digital Images
Study and Analysis o Edge Detection and Implementation o Fuzzy Set Theory Based Edge Detection Technique in Digital Images Anju K S Assistant Proessor, Department o Computer Science Baselios Mathews II
More informationDenoising an Image by Denoising its Components in a Moving Frame
Denoising an Image by Denoising its Components in a Moving Frame Gabriela Ghimpețeanu 1, Thomas Batard 1, Marcelo Bertalmío 1, and Stacey Levine 2 1 Universitat Pompeu Fabra, Spain 2 Duquesne University,
More informationFace Detection for Automatic Avatar Creation by using Deformable Template and GA
Face Detection or Automatic Avatar Creation by using Deormable Template and GA Tae-Young Park*, Ja-Yong Lee **, and Hoon Kang *** * School o lectrical and lectronics ngineering, Chung-Ang University, Seoul,
More informationNotes 9: Optical Flow
Course 049064: Variational Methods in Image Processing Notes 9: Optical Flow Guy Gilboa 1 Basic Model 1.1 Background Optical flow is a fundamental problem in computer vision. The general goal is to find
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 informationDistribution Fields with Adaptive Kernels for Large Displacement Image Alignment
MEARS et al.: DISTRIBUTION FIELDS WITH ADAPTIVE KERNELS 1 Distribution Fields with Adaptive Kernels or Large Displacement Image Alignment Benjamin Mears bmears@cs.umass.edu Laura Sevilla Lara bmears@cs.umass.edu
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 informationAutomatic Video Segmentation for Czech TV Broadcast Transcription
Automatic Video Segmentation or Czech TV Broadcast Transcription Jose Chaloupka Laboratory o Computer Speech Processing, Institute o Inormation Technology and Electronics Technical University o Liberec
More informationGlobal Minimization of the Active Contour Model with TV-Inpainting and Two-Phase Denoising
Global Minimization of the Active Contour Model with TV-Inpainting and Two-Phase Denoising Shingyu Leung and Stanley Osher Department of Mathematics, UCLA, Los Angeles, CA 90095, USA {syleung, sjo}@math.ucla.edu
More informationEXPERIMENTS ON GEOMETRIC IMAGE ENHANCEMENT. aco at
EXPERIMENTS ON GEOMETRIC IMAGE ENHANCEMENT Guillermo Sapiro Hewlett-Packard Labs 1501 Page Mill Road Palo Alto, CA 94304 guille@lids.mit.edu Allen Tannenbaum, Yu-La You, Mos Kaueh Department of Electrical
More informationOutlines. Medical Image Processing Using Transforms. 4. Transform in image space
Medical Image Processing Using Transforms Hongmei Zhu, Ph.D Department of Mathematics & Statistics York University hmzhu@yorku.ca Outlines Image Quality Gray value transforms Histogram processing Transforms
More informationSolution of the Hyperbolic Partial Differential Equation on Graphs and Digital Spaces: A Klein Bottle a Projective Plane and a 4D Sphere
International Journal o Discrete Mathematics 2017; 2(3): 88-94 http://wwwsciencepublishinggroupcom/j/dmath doi: 1011648/jdmath2017020315 olution o the Hyperbolic Partial Dierential Equation on Graphs and
More informationIterative Estimation of 3D Transformations for Object Alignment
Iterative Estimation of 3D Transformations for Object Alignment Tao Wang and Anup Basu Department of Computing Science, Univ. of Alberta, Edmonton, AB T6G 2E8, Canada Abstract. An Iterative Estimation
More informationA MULTI-LEVEL IMAGE DESCRIPTION MODEL SCHEME BASED ON DIGITAL TOPOLOGY
In: Stilla U et al (Eds) PIA7. International Archives o Photogrammetry, Remote Sensing and Spatial Inormation Sciences, 36 (3/W49B) A MULTI-LEVEL IMAGE DESCRIPTION MODEL SCHEME BASED ON DIGITAL TOPOLOG
More informationA Novel Field-source Reverse Transform for Image Structure Representation and Analysis
A Novel Field-source Reverse Transform for Image Structure Representation and Analysis X. D. ZHUANG 1,2 and N. E. MASTORAKIS 1,3 1. WSEAS Headquarters, Agiou Ioannou Theologou 17-23, 15773, Zografou, Athens,
More informationOptical flow estimation on image sequences with differently exposed frames. Tomas Bengtsson Tomas McKelvey Konstantin Lindström
Optical low estimation on image sequences with dierently exposed rames Tomas Bengtsson Tomas McKelvey Konstantin Lindström Optical Engineering 54(9), 093103 (September 015) Optical low estimation on image
More informationMULTISCALE IMAGE REGISTRATION. Dana Paquin. Doron Levy. Eduard Schreibmann. Lei Xing. (Communicated by Yang Kuang)
MATHEMATICAL BIOSCIENCES http://www.mbejournal.org/ AND ENGINEERING Volume 3, Number 2, April 2006 pp. 389 418 MULTISCALE IMAGE REGISTRATION Dana Paquin Department of Mathematics, Stanford University Stanford,
More informationImage Registration. Prof. Dr. Lucas Ferrari de Oliveira UFPR Informatics Department
Image Registration Prof. Dr. Lucas Ferrari de Oliveira UFPR Informatics Department Introduction Visualize objects inside the human body Advances in CS methods to diagnosis, treatment planning and medical
More informationAutomated Modelization of Dynamic Systems
Automated Modelization o Dynamic Systems Ivan Perl, Aleksandr Penskoi ITMO University Saint-Petersburg, Russia ivan.perl, aleksandr.penskoi@corp.imo.ru Abstract Nowadays, dierent kinds o modelling settled
More informationA Novel Image Transform Based on Potential field Source Reverse for Image Analysis
A Novel Image Transform Based on Potential field Source Reverse for Image Analysis X. D. ZHUANG 1,2 and N. E. MASTORAKIS 1,3 1. WSEAS Headquarters, Agiou Ioannou Theologou 17-23, 15773, Zografou, Athens,
More informationFast Image Registration via Joint Gradient Maximization: Application to Multi-Modal Data
MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.merl.com Fast Image Registration via Joint Gradient Maximization: Application to Multi-Modal Data Xue Mei, Fatih Porikli TR-19 September Abstract We
More informationKANGAL REPORT
Individual Penalty Based Constraint handling Using a Hybrid Bi-Objective and Penalty Function Approach Rituparna Datta Kalyanmoy Deb Mechanical Engineering IIT Kanpur, India KANGAL REPORT 2013005 Abstract
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 informationA Hierarchical Quality-Dependent Approach Toward Establishing A Seamless Nationwide Topographic Database
A Hierarchical Quality-Dependent Approach Toward Establishing A Seamless Nationwide Topographic Database Sagi Dalyot, Arieal Gershkovich, Yerach Doytsher Mapping and Geo-Inormation Engineering Technion,
More informationsmooth coefficients H. Köstler, U. Rüde
A robust multigrid solver for the optical flow problem with non- smooth coefficients H. Köstler, U. Rüde Overview Optical Flow Problem Data term and various regularizers A Robust Multigrid Solver Galerkin
More information2 DETERMINING THE VANISHING POINT LOCA- TIONS
IEEE TRANSACTIONS ON PATTERN ANALYSIS AND MACHINE INTELLIGENCE, VOL.??, NO.??, DATE 1 Equidistant Fish-Eye Calibration and Rectiication by Vanishing Point Extraction Abstract In this paper we describe
More informationLevel lines based disocclusion
Level lines based disocclusion Simon Masnou Jean-Michel Morel CEREMADE CMLA Université Paris-IX Dauphine Ecole Normale Supérieure de Cachan 75775 Paris Cedex 16, France 94235 Cachan Cedex, France Abstract
More informationAn Introduc+on to Mathema+cal Image Processing IAS, Park City Mathema2cs Ins2tute, Utah Undergraduate Summer School 2010
An Introduc+on to Mathema+cal Image Processing IAS, Park City Mathema2cs Ins2tute, Utah Undergraduate Summer School 2010 Luminita Vese Todd WiCman Department of Mathema2cs, UCLA lvese@math.ucla.edu wicman@math.ucla.edu
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 informationA Method of Sign Language Gesture Recognition Based on Contour Feature
Proceedings o the World Congress on Engineering and Computer Science 014 Vol I WCECS 014, -4 October, 014, San Francisco, USA A Method o Sign Language Gesture Recognition Based on Contour Feature Jingzhong
More informationGesture Recognition using a Probabilistic Framework for Pose Matching
Gesture Recognition using a Probabilistic Framework or Pose Matching Ahmed Elgammal Vinay Shet Yaser Yacoob Larry S. Davis Computer Vision Laboratory University o Maryland College Park MD 20742 USA elgammalvinayyaserlsd
More informationIsophote-Based Interpolation
Isophote-Based Interpolation Bryan S. Morse and Duane Schwartzwald Department of Computer Science, Brigham Young University 3361 TMCB, Provo, UT 84602 {morse,duane}@cs.byu.edu Abstract Standard methods
More informationYing-Chuan Road, Tamsui, Taipei County, Taiwan, 251, R.O.C. I-Lan, Taiwan, 260, R.O.C.
Materials Science Forum Online: 006-0-5 ISSN: 66-975, Vols. 505-507, pp 697-70 doi:0.408/www.scientiic.net/msf.505-507.697 006 Trans Tech Publications, Switzerland Application o Streamline Method and ANFIS
More informationISSN: X International Journal of Advanced Research in Electronics and Communication Engineering (IJARECE) Volume 6, Issue 8, August 2017
ENTROPY BASED CONSTRAINT METHOD FOR IMAGE SEGMENTATION USING ACTIVE CONTOUR MODEL M.Nirmala Department of ECE JNTUA college of engineering, Anantapuramu Andhra Pradesh,India Abstract: Over the past existing
More informationImage Segmentation and Registration
Image Segmentation and Registration Dr. Christine Tanner (tanner@vision.ee.ethz.ch) Computer Vision Laboratory, ETH Zürich Dr. Verena Kaynig, Machine Learning Laboratory, ETH Zürich Outline Segmentation
More informationMonocular 3D human pose estimation with a semi-supervised graph-based method
Monocular 3D human pose estimation with a semi-supervised graph-based method Mahdieh Abbasi Computer Vision and Systems Laboratory Department o Electrical Engineering and Computer Engineering Université
More informationImage Processing and related PDEs Lecture 1: Introduction to image processing
Image Processing and related PDEs Lecture 1: Introduction to image processing Yves van Gennip School of Mathematical Sciences, University of Nottingham Minicourse on Image Processing and related PDEs University
More informationFuzzy Inference System based Edge Detection in Images
Fuzzy Inference System based Edge Detection in Images Anjali Datyal 1 and Satnam Singh 2 1 M.Tech Scholar, ECE Department, SSCET, Badhani, Punjab, India 2 AP, ECE Department, SSCET, Badhani, Punjab, India
More informationAn Active Contour Model without Edges
An Active Contour Model without Edges Tony Chan and Luminita Vese Department of Mathematics, University of California, Los Angeles, 520 Portola Plaza, Los Angeles, CA 90095-1555 chan,lvese@math.ucla.edu
More informationUnstructured Mesh Generation for Implicit Moving Geometries and Level Set Applications
Unstructured Mesh Generation for Implicit Moving Geometries and Level Set Applications Per-Olof Persson (persson@mit.edu) Department of Mathematics Massachusetts Institute of Technology http://www.mit.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 informationAC : DEVELOPMENT OF A ROBOTIC PLATFORM FOR TEACH- ING MODEL-BASED DESIGN TECHNIQUES IN DYNAMICS AND CON- TROL PROGRAM
AC 011-714: DEVELOPMENT OF A ROBOTIC PLATFORM FOR TEACH- ING MODEL-BASED DESIGN TECHNIQUES IN DYNAMICS AND CON- TROL PROGRAM Bingen Yang, University o Southern Caliornia Dr. Bingen Yang is Proessor o Aerospace
More informationRobust Realignment of fmri Time Series Data
Robust Realignment of fmri Time Series Data Ben Dodson bjdodson@stanford.edu Olafur Gudmundsson olafurg@stanford.edu December 12, 2008 Abstract FMRI data has become an increasingly popular source for exploring
More informationRoad Sign Analysis Using Multisensory Data
Road Sign Analysis Using Multisensory Data R.J. López-Sastre, S. Lauente-Arroyo, P. Gil-Jiménez, P. Siegmann, and S. Maldonado-Bascón University o Alcalá, Department o Signal Theory and Communications
More informationA Cartesian Grid Method with Transient Anisotropic Adaptation
Journal o Computational Physics 179, 469 494 (2002) doi:10.1006/jcph.2002.7067 A Cartesian Grid Method with Transient Anisotropic Adaptation F. E. Ham, F. S. Lien, and A. B. Strong Department o Mechanical
More informationGRID WARPING IN TOTAL VARIATION IMAGE ENHANCEMENT METHODS. Andrey Nasonov, and Andrey Krylov
GRID WARPING IN TOTAL VARIATION IMAGE ENHANCEMENT METHODS Andrey Nasonov, and Andrey Krylov Lomonosov Moscow State University, Moscow, Department of Computational Mathematics and Cybernetics, e-mail: nasonov@cs.msu.ru,
More informationHybrid Spline-based Multimodal Registration using a Local Measure for Mutual Information
Hybrid Spline-based Multimodal Registration using a Local Measure for Mutual Information Andreas Biesdorf 1, Stefan Wörz 1, Hans-Jürgen Kaiser 2, Karl Rohr 1 1 University of Heidelberg, BIOQUANT, IPMB,
More informationTHE FINANCIAL CALCULATOR
Starter Kit CHAPTER 3 Stalla Seminars THE FINANCIAL CALCULATOR In accordance with the AIMR calculator policy in eect at the time o this writing, CFA candidates are permitted to use one o two approved calculators
More information9. Reviewing Printed Circuit Board Schematics with the Quartus II Software
November 2012 QII52019-12.1.0 9. Reviewing Printed Circuit Board Schematics with the Quartus II Sotware QII52019-12.1.0 This chapter provides guidelines or reviewing printed circuit board (PCB) schematics
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 informationWhole Body MRI Intensity Standardization
Whole Body MRI Intensity Standardization Florian Jäger 1, László Nyúl 1, Bernd Frericks 2, Frank Wacker 2 and Joachim Hornegger 1 1 Institute of Pattern Recognition, University of Erlangen, {jaeger,nyul,hornegger}@informatik.uni-erlangen.de
More informationAnnales UMCS Informatica AI 1 (2003) UMCS. Registration of CT and MRI brain images. Karol Kuczyński, Paweł Mikołajczak
Annales Informatica AI 1 (2003) 149-156 Registration of CT and MRI brain images Karol Kuczyński, Paweł Mikołajczak Annales Informatica Lublin-Polonia Sectio AI http://www.annales.umcs.lublin.pl/ Laboratory
More informationFormalizing Cardinality-based Feature Models and their Staged Configuration
Formalizing Cardinality-based Feature Models and their Staged Coniguration Krzyszto Czarnecki, Simon Helsen, and Ulrich Eisenecker 2 University o Waterloo, Canada 2 University o Applied Sciences Kaiserslautern,
More informationXavier: A Robot Navigation Architecture Based on Partially Observable Markov Decision Process Models
Xavier: A Robot Navigation Architecture Based on Partially Observable Markov Decision Process Models Sven Koenig and Reid G. Simmons Carnegie Mellon University School o Computer Science Pittsburgh, PA
More informationImage Processing with Nonparametric Neighborhood Statistics
Image Processing with Nonparametric Neighborhood Statistics Ross T. Whitaker Scientific Computing and Imaging Institute School of Computing University of Utah PhD Suyash P. Awate University of Pennsylvania,
More informationExpectation Maximization and Total Variation Based Model for Computed Tomography Reconstruction from Undersampled Data
Expectation Maximization and Total Variation Based Model for Computed Tomography Reconstruction from Undersampled Data Ming Yan and Luminita A. Vese Department of Mathematics, University of California,
More informationFast 3D Mean Shift Filter for CT Images
Fast 3D Mean Shift Filter for CT Images Gustavo Fernández Domínguez, Horst Bischof, and Reinhard Beichel Institute for Computer Graphics and Vision, Graz University of Technology Inffeldgasse 16/2, A-8010,
More informationA non-rigid registration method for mouse whole body skeleton registration
A non-rigid registration method or mouse whole body skeleton registration Di Xiao* a, David Zahra b, Pierrick Bourgeat a, Paula Berghoer b, Oscar Acosta amayo a, Catriona Wimberley b, Marie Claude Gregoire
More informationLecture 2 September 3
EE 381V: Large Scale Optimization Fall 2012 Lecture 2 September 3 Lecturer: Caramanis & Sanghavi Scribe: Hongbo Si, Qiaoyang Ye 2.1 Overview of the last Lecture The focus of the last lecture was to give
More informationPATH PLANNING OF UNMANNED AERIAL VEHICLE USING DUBINS GEOMETRY WITH AN OBSTACLE
Proceeding o International Conerence On Research, Implementation And Education O Mathematics And Sciences 2015, Yogyakarta State University, 17-19 May 2015 PATH PLANNING OF UNMANNED AERIAL VEHICLE USING
More informationAnalysis of Image and Video Using Color, Texture and Shape Features for Object Identification
IOSR Journal of Computer Engineering (IOSR-JCE) e-issn: 2278-0661,p-ISSN: 2278-8727, Volume 16, Issue 6, Ver. VI (Nov Dec. 2014), PP 29-33 Analysis of Image and Video Using Color, Texture and Shape Features
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 informationComparison between Various Edge Detection Methods on Satellite Image
Comparison between Various Edge Detection Methods on Satellite Image H.S. Bhadauria 1, Annapurna Singh 2, Anuj Kumar 3 Govind Ballabh Pant Engineering College ( Pauri garhwal),computer Science and Engineering
More informationIntensity gradient based registration and fusion of multi-modal images
Intensity gradient based registration and fusion of multi-modal images Eldad Haber 1 and Jan Modersitzki 2 1 Mathematics and Computer Science, Emory University, Atlanta, GA, USA, haber@mathcs.emory,edu
More informationImproving Alignment of Faces for Recognition
Improving Alignment o Faces or Recognition Md. Kamrul Hasan Département de génie inormatique et génie logiciel École Polytechnique de Montréal, Québec, Canada md-kamrul.hasan@polymtl.ca Christopher J.
More informationAn Automatic Sequential Smoothing Method for Processing Biomechanical Kinematic Signals
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
More informationKernel Density Estimation and Intrinsic Alignment for Knowledge-driven Segmentation: Teaching Level Sets to Walk
C. Rasmussen et al. (Eds.), Pattern Recognition, Tübingen, Sept. 2004. LNCS Vol. 3175, pp. 36 44. c Springer Kernel Density Estimation and Intrinsic Alignment for Knowledge-driven Segmentation: Teaching
More informationMethodological progress in image registration for ventilation estimation, segmentation propagation and multi-modal fusion
Methodological progress in image registration for ventilation estimation, segmentation propagation and multi-modal fusion Mattias P. Heinrich Julia A. Schnabel, Mark Jenkinson, Sir Michael Brady 2 Clinical
More informationLevel-set and ALE Based Topology Optimization Using Nonlinear Programming
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, USA Level-set and ALE Based Topology Optimization Using Nonlinear Programming Shintaro Yamasaki
More informationMultimodality Imaging for Tumor Volume Definition in Radiation Oncology
81 There are several commercial and academic software tools that support different segmentation algorithms. In general, commercial software packages have better implementation (with a user-friendly interface
More informationDeformable Registration Using Scale Space Keypoints
Deformable Registration Using Scale Space Keypoints Mehdi Moradi a, Purang Abolmaesoumi a,b and Parvin Mousavi a a School of Computing, Queen s University, Kingston, Ontario, Canada K7L 3N6; b Department
More informationContent Based Image Retrieval Algorithm Based On The Dual-Tree Complex Wavelet Transform: Efficiency Analysis
Content Based Image etrieval Algorithm Based On The Dual-Tree Complex Wavelet Transorm: Eiciency Analysis STELLA VETOVA, IVAN IVANOV 2 Institute o Inormation and Communication Technologies, 2 Telecommunication
More informationAn Efficient, Geometric Multigrid Solver for the Anisotropic Diffusion Equation in Two and Three Dimensions
1 n Efficient, Geometric Multigrid Solver for the nisotropic Diffusion Equation in Two and Three Dimensions Tolga Tasdizen, Ross Whitaker UUSCI-2004-002 Scientific Computing and Imaging Institute 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 informationECSE 626 Project Report Multimodality Image Registration by Maximization of Mutual Information
ECSE 626 Project Report Multimodality Image Registration by Maximization of Mutual Information Emmanuel Piuze McGill University Montreal, Qc, Canada. epiuze@cim.mcgill.ca Abstract In 1997, Maes et al.
More information