3D SMAP Algorithm. April 11, 2012
|
|
- Martin Terry
- 5 years ago
- Views:
Transcription
1 3D SMAP Algorithm April 11, 2012 Baed on the original SMAP paper [1]. Thi report extend the tructure of MSRF into 3D. The prior ditribution i modified to atify the MRF property. In addition, an iterative trategy i propoed to refine the egmentation. 1 Goal Segment 3D image volume into region of ditinct tatitical behavior. 2 Baic Model Y i oberved 3D image volume containing ditinct texture. X i unoberved field containing the cla of each pixel. S i 3D lattice of point. i a member of S. The lower cae letter x, y denote the correponding determinitic realization. 3 Approach Outline 1. Auming the behavior of each oberved pixel in Y i dependent on a correponding unoberved label pixel in X. 1
2 2. The dependence i pecified through p y x (y x), which i modeled by Gauian mixture ditribution (GMM) in our experiment. The parameter of GMM i computed uing cluter algorithm decribed by [2]. 3. Prior knowledge about the region will be pecified by the prior ditribution p(x), which i modeled by multicale random field (MSRF). 4. Solve following optimization problem which minimize the average cot of an erroneou egmentation. ˆx = arg min x E[C(X, x) Y = y] (1) 5. Once the egmentation i obtained, we can re-etimate the GMM parameter for each cla uing only a ubet of the original data, whoe member have the ame egmentation label. Therefore, we can iteratively refine the egmentation reult until the termination condition i atified. 4 Multicale Random Field Multicale Random Field (MSRF) i compoed of a erie of 3D random field at varying cale or reolution. A demontration picture i hown by Fig 1. At each cale, n, the egmentation or labeling i denoted by the random field X, and the et of lattice point i denoted by S. In particular, X (0) i aumed to be the finet cale random field with each point correponding to a ingle image pixel. Each label at the next coarer cale X (1) then correpond to a group of eight point in the original image. Therefore, the number of point in S (1) i 1/8 the number of point in S (0). Aume the equence of random field from coare to fine cale form a Markov chain. Then P (X = x X (l) = x (l), l > n) (2) =P (X = x X (n+1) = x (n+1) ) =P x x (n+1)(x x (n+1) ) Define the cot function a follow: C SMAP (X, x) = 1 L n 1 C n (X, x) (3) n=0
3 Figure 1: 3D random field with three different cale. Each cube repreent a pixel in the 3D image volume where L i the coaret cale in X and C n (X, x) i the penalty function aociate with labeling error at the n th cale. Suppoe K th cale i the unique cale uch that X (K) x (K), and X (i) = x (i) for all i > K. Define the function C n a 1, if n K C n (X, x) = 0, if n > K Uing (3) and (4) We can compute the total cot C SMAP (X, x) = 2 K. Intuitively, C SMAP (X, x) can be interpretated a the width of the larget grouping of miclaified pixel. We can determine the etimator which minimize the propoed cot function by evaluating (1). ˆx = arg min E[C SMAP (X, x) Y = y] (5) x L = arg min 2 n {1 P (X (i) = x (i), i n Y = y)} x = arg max x n=0 L 2 n {P (X (i) = x (i), i n Y = y)} n=0 Uing recurive approach, [1] ha proved the following olution for (6) ˆx = arg max x {log p x x (n+1),y(x ˆx (n+1), y) + ɛ(x )} (6) where ɛ i a econd order term which vanihe to 0 a number of pixel N increae. i.e. lim N ɛ(x ) = 0. (4)
4 Ignoring ɛ when N i large reult in the following recurive equation: ˆx (L) = arg max x (L) log p x (L) y(x (L) y) (7) ˆx = arg max x {log p x x (n+1),y(x ˆx (n+1), y)} (8) Appying the Baye rule, the Markov propertie of X, and auming that X (L) i uniformly ditributed, the SMAP recurion may be rewritten in a form which i more eaily computed. ˆx (L) = arg max x (L) log p y x (L)(y x (L) ) (9) ˆx = arg max x {log p y x (y x ) + log p x x (n+1)(x ˆx (n+1) )} (10) Now we introduce 3D SMAP algorithm to etimate the olution to (9) and (10). 5 3D SMAP algorithm Thi algorithm pecify the conditional denity p y x (y x ) together with the coare to fine cale tranition denitie p x x (n+1)(x ˆx (n+1) ). A hybrid model, which incorporate 3D quadtree and 3D graph i ued to decribe the tructure of MSRF. The algorithm etimate the parameter of MSRF during the egmentation proce. The computation i baed on following aumption and propertie of Markov Random Field (MRF). Firt, the oberved pixel are conditionaly independent given their cla label. Then the conditional denity function for the image ha the form. p y x (0)(y x (0) ) = p y x (0) S (0) (y x(0) ) (11) where p y x ( k) i the conditional denity function for an individual pixel given the cla (0) label k. We ue a multivariate Gauian mixture denity in experiment, but other ditribution can be applied a well. Second, the labeling in X will be conditionally independent, given the labeling in X (n+1).
5 Furthermore, each labeling X will only be dependent on a local neighborhood of labeling at the next coarer cale. That i where p x x (n+1) x (n+1). p x x (n+1)(x ˆx (n+1) ) = p x S i the probability denity for x x (n+1) (x ˆx (n+1) ) (12) given it neighbor at the coarer cale The choice of neighborhood i important, which affect the tructure of multicale pyramid a well a the pecification of probability denity function decribed earlier in thi ection. 1) Compute p y x (y x ) with 3D quadtree The tructure of the quadtree i hown in Fig 2. Each point in the pyramid i only dependent on a ingle point at the coarer cale. According to Hammerley-Clifford theorem, the probability denity function can be characterized by Gibb ditribution. We ue the following function to model the probability that X ha cla m, given that it parent i of cla k. Figure 2: 3D random field with two ucceive cale. Each cube repreent a pixel in the 3D image volume. The red pixel in cale n = 0 have the ame parent in coarer level n = 1, which i alo in red. p x x (n+1) (m k) = 1 Z exp( θ 0 δ m k ) (13) where δ m k = 1 if m k, otherwie δ m k = 0. Z = M m=1 exp{ θ 0 δ m k } i called partition function, where M i the number of poible clae. θ 0 > 0 i the regularization parameter
6 of cale n. Large θ 0 reult higher probability that the label remain the ame from cale n + 1 to n. An important property of the quadtree tructure i that the conditional ditribution of Y given X, ha a product form that may be computed uing a imple fine-to-coare recurion. Let y be the et of leave of the 3D quadtree that are on the branch tarting at S. [1] prove that the conditional denity for Y given X ha the product form Furthermore, the denity function p y x p y (n+1) x (n+1) p y x (y x ) = p y (y x x ) (14) S (y (n+1) k) = r d 1 () m=1 may be computed uing the following recurion, M p y r x r (y r m)p x where d 1 () i a et of point whoe firt ucceive parent i. r x (n+1) r (m k) (15) Conider the dynamic range, we take the log likelihood function defined a follow. l (k) log p y x (y k) (16) Subtitute the tranition ditribution of (13) into (15) and converting to log likelihood function yield the new recurion l (0) (k) = log p y (0) l (n+1) (k) = r d 1 () x (0)(y(0) k) (17) M 1 log{ Z exp(l r (m) θ 0 δ m k )} m=1 Uing (17), we can rewrite the firt term of right-hand-ide of (10) a log p y x (y x ) = l S 2) Compute p x x (n+1)(x ˆx (n+1) ) with 3D graph (x ) (18) A diadvantage of the 3D quadtree model i that pixel that are patially adjacent may not have common neighbor at the next coarer cale. In order to compenate the weakne of dicontinuity of boundary, we introduce 3D graph to model the tranition probability
7 p x x (n+1)(x ˆx (n+1) ). The main idea i to increae the number of coare cale neighbor for each pixel. According to the decription in 1), point ha only one coare cale neighbor, which i the parent node in the 3D quad tree. Now we expand the number of coare cale neighbor from one to four a hown in Fig. 3. In order to expre the poition of the neighbor, we explicitly denote a pixel at cale D a = (t, i, j), where t, i, j are the frame, row and column indice. The four neigbor at cale D + 1 may then be computed uing the function odd(i) = 1 if i i odd and -1 if i i even and the greatet maller integer function. (a) (b) Figure 3: 3D random field with two ucceive cale. Each cube repreent a pixel in the 3D image volume. a) how 32 pixel in level D, b) how 4 pixel in coarer level D + 1 correponding to a). In particular, pixel in level D ha four coarer level neighbor: 0, 1, 2, 3. 0 = ( t/2, i/2, j/2 ) 1 = ( t/2, i/2, j/2 ) + (odd(t), 0, 0) 2 = ( t/2, i/2, j/2 ) + (0, odd(i), 0) 3 = ( t/2, i/2, j/2 ) + (0, 0, odd(j))
8 The tranition function which we chooe for thi pyramid graph i a natural extenion of the tranition function ued for the quadtree baed model. In order to ditinguih (13) from the following formula, we ue p intead. p x x (n+1) =P (X = 1 Z exp( (m m r ), r d() (19) = m X (n+1) r = m r ) r d() β r δ m mr ) Similar to (13), Z = M m=1 exp{ r d() β r δ m mr }, d() = { 0, 1, 2, 3 } i a et of point contain all the neighbor of in the next coarer level. β r can be viewed a a weight parameter aociated with the neighbor point r in cale n + 1. The larger the β r, the larger the probability that the label of will be the ame a r. Beide, we can make the tranition moother by adding neighbor point in the ame cale. Define d 0 () = { 0,..., 5} be a et containing ix neighbor point in the ame cale D a. See Fig. 4 for a viual demontration. The coordinate of are explicitly given by Then (19) can be rewritten a p x 0 = (t, i, j 1), 1 = (t, i, j + 1) 2 = (t, i 1, j), 3 = (t, i + 1, j) 4 = (t 1, i, j), 5 = (t + 1, i, j) x (n+1),x =P (X = 1 Z exp( (m m r, m q ), r d(), q d 0 () (20) = m X (n+1) r r d() = m r, X q = m q ) β r δ m mr q d 0 () β q δ m mq ) where Z = M m=1 exp{ r d() β r δ m mr q d 0 () β q δ m mq }. For implicity, we can ue identical weight for point in the ame cale, i.e. θ 1, q, β q θ 2. r, β r Uing thi hybrid pyramid tructure of MSRF, the conditional likelihood of (8) and (10) can
9 Figure 4: 3D random field of one cale. Each cube repreent a pixel in the 3D image volume. i a pixel in level D with ix neighbor in the ame level: 0, 1, 2, 3, 4, 5. be expreed in the form log p y,x x (n+1)(y, x ˆx (n+1) ) = l S (x ) + log p x x (n+1),x which reult in the following formula for the SMAP etimate of X : (x ˆx (n+1), ˆx ) (21) ˆx (L) ˆx = arg max 1 k M l(l) (k) (22) = arg max 1 k M {l = arg max 1 k M {l (k) + log p x (k) log Z θ x (n+1),x 1 r d() (k ˆx (n+1), ˆx δ m mr θ 2 )} q d 0 () δ m mq } 6 Etimation of Tranition Probability Parameter We introduced three parameter θ = [θ 0, θ 1, θ 2 ] to determine the tranition probability denity function baed on two different model of pyramid tructure. To be pecific, θ 0 are required for fine-to-coare in computing the log likelihood function (17). And the θ 1, θ 2 are required for coare-to-fine operation in the SMAP egmentation (22). The parameter are aigned with ome initial value and can be updated in the proce of egmentation by
10 olving a joint optimization problem pecified by (17) and (22). Thi problem can be olved uing ICM [3]. 7 Etimation of GMM parameter The tatitical behavior of each egment i modeled by a Gauian mixture ditribution (GMM) with parameter η k = [π k, µ k, R k ], k = 1,..., M, where M i the poible number of clae. Let π k R N, µ k R N, R k R N N, where N i the number of Gauian in the mixture denity. And N 1 i=0 π k i = 1, where π ki i the i th entry of π k. Then the probability denity function of one pixel given it label k can be written explicitly a p y x,η(y k, η k ) = N 1 i=0 For impler cae when N = 1, (23) can be reduced to p y x,η(y k, η k ) = π { ki (2π) R N/2 k i 1/2 exp 1 } 2 (y µ ki ) T R 1 k i (y µ ki ) 1 { (2π) R k 1/2 exp 1 } 1/2 2 (y µ k ) T R 1 k (y µ k ) (23) (24) (23) and (24) are ued to compute the log likelihood function defined in (16). We tart the egmentation with ome initially value of η k. After applying the 3D SMAP algorithm to the entire data et Y, we obtain the label of pixel X. Now we can divide the Y into M dijoint ubet. For each ubet, we can ue [2] to etimate new η k, with which we can apply the 3D SMAP algorithm again to refine the egmentation. Thi proce i repeated until the change of GMM parameter between two ucceive iteration i le than certain amount of value. For intance, we can top the iteration when the following criteria i atified. M µ k µ k ɛ 3 (25) k=1 where i the infinity norm and ɛ 3 = 1. The complete SMAP egmentation algorithm with parameter etimation i ummarized below.
11 1. Set the initial parameter value for all n, θ 0 = 5 1 n, θ 1 = 4, and θ 2 = 2, n = 0,..., L. 2. Compute the likelihood function uing (17). 3. Compute ˆx (L) (k) uing (22). 4. For cale n = L 1 to n = 0 (a) Set I = deired number of iteration i. Compute ˆx uing (22). ii. Update ˆθ. iii. Reduce I by one. (b) Proceed if I i negative or no further change happen, otherwie go back to (a) 5. Re-etimate GMM parameter 6. Repeat tep 2 through 5 until (25) i atified. 8 Experiment Reult Data et: NiCrAl 1 kev Simulation, ource from BlueQuartz Software. Data Decription: The data et conit of 257 lice of 513x513 pixel TIFF image. Each image repreent a 15 nm thick lice, and the pixel dimenion are alo 15x15 nm. Microcope parameter i 1 kev. Experiment condition: The algorithm i implemented in C. The experiment for thi particular dataet i performed on a PC with 2.0GHz clock peed, 1.5GB RAM. The running time of one iteration i econd including tep 1 to 4 decribed above, importing and exporting image. Segmentation reult of the 201 t to 205 th image are provided a follow. The three different reult are obtained by applying 2D SMAP, 3D SMAP without GMM parameter re-etimation and 3D SMAP with one iteration of GMM parameter re-etimation.
12 (a) (b) (c) (d) Figure 5: 201(a)Orignial image; (b)smap; (c)3d SMAP; (d)3d SMAP uing re- etimated GMM parameter with number of iteration equal one
13 (a) (b) (c) (d) Figure 6: 202(a)Orignial image; (b)smap; (c)3d SMAP; (d)3d SMAP uing re- etimated GMM parameter with number of iteration equal one
14 (a) (b) (c) (d) Figure 7: 203(a)Orignial image; (b)smap; (c)3d SMAP; (d)3d SMAP uing re- etimated GMM parameter with number of iteration equal one
15 (a) (b) (c) (d) Figure 8: 204(a)Orignial image; (b)smap; (c)3d SMAP; (d)3d SMAP uing re- etimated GMM parameter with number of iteration equal one
16 (a) (b) (c) (d) Figure 9: 205(a)Orignial image; (b)smap; (c)3d SMAP; (d)3d SMAP uing re- etimated GMM parameter with number of iteration equal one
17 9 reference [1] C. A. Bouman and M. Shapiro, A multicale random eld model for Bayeian image egmentation, IEEE Tran. Image Proceing 3: , [2] C. A. Bouman, M. Shapiro. G. W. Cook, C. B. Atkin, H. Cheng, Cluter: An unupervied algorithm for modeling gauian mixture, In cluter [3] C. A. Bouman, B. Liu, Multiple Reolution Segmentation of Textured Image, IEEE Tran. Pattern Analyi and Machine Intelligence, vol. 13, no. 2, pp , Feb
Trainable Context Model for Multiscale Segmentation
Trainable Context Model for Multicale Segmentation Hui Cheng and Charle A. Bouman School of Electrical and Computer Engineering Purdue Univerity Wet Lafayette, IN 47907-1285 {hui, bouman}@ ecn.purdue.edu
More informationMarkov Random Fields in Image Segmentation
Preented at SSIP 2011, Szeged, Hungary Markov Random Field in Image Segmentation Zoltan Kato Image Proceing & Computer Graphic Dept. Univerity of Szeged Hungary Zoltan Kato: Markov Random Field in Image
More informationMAT 155: Describing, Exploring, and Comparing Data Page 1 of NotesCh2-3.doc
MAT 155: Decribing, Exploring, and Comparing Data Page 1 of 8 001-oteCh-3.doc ote for Chapter Summarizing and Graphing Data Chapter 3 Decribing, Exploring, and Comparing Data Frequency Ditribution, Graphic
More informationIMPROVED JPEG DECOMPRESSION OF DOCUMENT IMAGES BASED ON IMAGE SEGMENTATION. Tak-Shing Wong, Charles A. Bouman, and Ilya Pollak
IMPROVED DECOMPRESSION OF DOCUMENT IMAGES BASED ON IMAGE SEGMENTATION Tak-Shing Wong, Charle A. Bouman, and Ilya Pollak School of Electrical and Computer Engineering Purdue Univerity ABSTRACT We propoe
More informationUniversität Augsburg. Institut für Informatik. Approximating Optimal Visual Sensor Placement. E. Hörster, R. Lienhart.
Univerität Augburg à ÊÇÅÍÆ ËÀǼ Approximating Optimal Viual Senor Placement E. Hörter, R. Lienhart Report 2006-01 Januar 2006 Intitut für Informatik D-86135 Augburg Copyright c E. Hörter, R. Lienhart Intitut
More informationMinimum congestion spanning trees in bipartite and random graphs
Minimum congetion panning tree in bipartite and random graph M.I. Otrovkii Department of Mathematic and Computer Science St. John Univerity 8000 Utopia Parkway Queen, NY 11439, USA e-mail: otrovm@tjohn.edu
More informationPerformance of a Robust Filter-based Approach for Contour Detection in Wireless Sensor Networks
Performance of a Robut Filter-baed Approach for Contour Detection in Wirele Senor Network Hadi Alati, William A. Armtrong, Jr., and Ai Naipuri Department of Electrical and Computer Engineering The Univerity
More information1 The secretary problem
Thi i new material: if you ee error, pleae email jtyu at tanford dot edu 1 The ecretary problem We will tart by analyzing the expected runtime of an algorithm, a you will be expected to do on your homework.
More informationOn combining Learning Vector Quantization and the Bayesian classifiers for natural textured images
On combining Learning Vector Quantization and the Bayeian claifier for natural textured image María Guiarro Dept. Ingeniería del Software e Inteligencia Artificial Facultad Informática Univeridad Complutene
More informationRouting Definition 4.1
4 Routing So far, we have only looked at network without dealing with the iue of how to end information in them from one node to another The problem of ending information in a network i known a routing
More informationLecture 14: Minimum Spanning Tree I
COMPSCI 0: Deign and Analyi of Algorithm October 4, 07 Lecture 4: Minimum Spanning Tree I Lecturer: Rong Ge Scribe: Fred Zhang Overview Thi lecture we finih our dicuion of the hortet path problem and introduce
More informationA note on degenerate and spectrally degenerate graphs
A note on degenerate and pectrally degenerate graph Noga Alon Abtract A graph G i called pectrally d-degenerate if the larget eigenvalue of each ubgraph of it with maximum degree D i at mot dd. We prove
More informationOperational Semantics Class notes for a lecture given by Mooly Sagiv Tel Aviv University 24/5/2007 By Roy Ganor and Uri Juhasz
Operational emantic Page Operational emantic Cla note for a lecture given by Mooly agiv Tel Aviv Univerity 4/5/7 By Roy Ganor and Uri Juhaz Reference emantic with Application, H. Nielon and F. Nielon,
More informationThe norm Package. November 15, Title Analysis of multivariate normal datasets with missing values
The norm Package November 15, 2003 Verion 1.0-9 Date 2002/05/06 Title Analyi of multivariate normal dataet with miing value Author Ported to R by Alvaro A. Novo . Original by Joeph
More informationxy-monotone path existence queries in a rectilinear environment
CCCG 2012, Charlottetown, P.E.I., Augut 8 10, 2012 xy-monotone path exitence querie in a rectilinear environment Gregory Bint Anil Mahehwari Michiel Smid Abtract Given a planar environment coniting of
More informationContents. shortest paths. Notation. Shortest path problem. Applications. Algorithms and Networks 2010/2011. In the entire course:
Content Shortet path Algorithm and Network 21/211 The hortet path problem: Statement Verion Application Algorithm (for ingle ource p problem) Reminder: relaxation, Dijktra, Variant of Dijktra, Bellman-Ford,
More informationBuilding a Compact On-line MRF Recognizer for Large Character Set using Structured Dictionary Representation and Vector Quantization Technique
202 International Conference on Frontier in Handwriting Recognition Building a Compact On-line MRF Recognizer for Large Character Set uing Structured Dictionary Repreentation and Vector Quantization Technique
More informationA Novel Feature Line Segment Approach for Pattern Classification
12th International Conference on Information Fuion Seattle, WA, USA, July 6-9, 2009 A Novel Feature Line Segment Approach for Pattern Claification Yi Yang Intitute of Integrated Automation Xi an Jiaotong
More informationLaboratory Exercise 6
Laboratory Exercie 6 Adder, Subtractor, and Multiplier The purpoe of thi exercie i to examine arithmetic circuit that add, ubtract, and multiply number. Each type of circuit will be implemented in two
More informationRepresentations and Transformations. Objectives
Repreentation and Tranformation Objective Derive homogeneou coordinate tranformation matrice Introduce tandard tranformation - Rotation - Tranlation - Scaling - Shear Scalar, Point, Vector Three baic element
More informationShortest Path Routing in Arbitrary Networks
Journal of Algorithm, Vol 31(1), 1999 Shortet Path Routing in Arbitrary Network Friedhelm Meyer auf der Heide and Berthold Vöcking Department of Mathematic and Computer Science and Heinz Nixdorf Intitute,
More informationStress-Blended Eddy Simulation (SBES) - A new Paradigm in hybrid RANS-LES Modeling
Stre-Blended Eddy Simulation (SBES) - A new Paradigm in hybrid RANS-LES Modeling Menter F.R. ANSYS Germany GmbH Introduction It i oberved in many CFD imulation that RANS model how inherent technology limitation
More informationGray-level histogram. Intensity (grey-level) transformation, or mapping. Use of intensity transformations:
Faculty of Informatic Eötvö Loránd Univerity Budapet, Hungary Lecture : Intenity Tranformation Image enhancement by point proceing Spatial domain and frequency domain method Baic Algorithm for Digital
More informationBayesian segmentation for damage image using MRF prior
Bayeian egmentation for damage image uing MRF prior G. Li 1, F.G. Yuan 1, R. Haftka and N. H. Kim 1 Department of Mechanical and Aeropace Engineering, North arolina State Univerity, Raleigh, N, 7695-791,
More informationIntroduction to PET Image Reconstruction. Tomographic Imaging. Projection Imaging. PET Image Reconstruction 11/6/07
Introduction to PET Image Recontruction Adam Aleio Nuclear Medicine Lecture Imaging Reearch Laboratory Diviion of Nuclear Medicine Univerity of Wahington Fall 2007 http://dept.wahington.edu/nucmed/irl/education.html
More informationLecture Outline. Global flow analysis. Global Optimization. Global constant propagation. Liveness analysis. Local Optimization. Global Optimization
Lecture Outline Global flow analyi Global Optimization Global contant propagation Livene analyi Adapted from Lecture by Prof. Alex Aiken and George Necula (UCB) CS781(Praad) L27OP 1 CS781(Praad) L27OP
More informationEdits in Xylia Validity Preserving Editing of XML Documents
dit in Xylia Validity Preerving diting of XML Document Pouria Shaker, Theodore S. Norvell, and Denni K. Peter Faculty of ngineering and Applied Science, Memorial Univerity of Newfoundland, St. John, NFLD,
More informationSLA Adaptation for Service Overlay Networks
SLA Adaptation for Service Overlay Network Con Tran 1, Zbigniew Dziong 1, and Michal Pióro 2 1 Department of Electrical Engineering, École de Technologie Supérieure, Univerity of Quebec, Montréal, Canada
More informationOn successive packing approach to multidimensional (M-D) interleaving
On ucceive packing approach to multidimenional (M-D) interleaving Xi Min Zhang Yun Q. hi ankar Bau Abtract We propoe an interleaving cheme for multidimenional (M-D) interleaving. To achieved by uing a
More informationFactor Graphs and Inference
Factor Graph and Inerence Sargur Srihari rihari@cedar.bualo.edu 1 Topic 1. Factor Graph 1. Factor in probability ditribution. Deriving them rom graphical model. Eact Inerence Algorithm or Tree graph 1.
More informationNew Structural Decomposition Techniques for Constraint Satisfaction Problems
New Structural Decompoition Technique for Contraint Satifaction Problem Yaling Zheng and Berthe Y. Choueiry Contraint Sytem Laboratory Univerity of Nebraka-Lincoln Email: yzheng choueiry@ce.unl.edu Abtract.
More informationStochastic Search and Graph Techniques for MCM Path Planning Christine D. Piatko, Christopher P. Diehl, Paul McNamee, Cheryl Resch and I-Jeng Wang
Stochatic Search and Graph Technique for MCM Path Planning Chritine D. Piatko, Chritopher P. Diehl, Paul McNamee, Cheryl Rech and I-Jeng Wang The John Hopkin Univerity Applied Phyic Laboratory, Laurel,
More informationDelaunay Triangulation: Incremental Construction
Chapter 6 Delaunay Triangulation: Incremental Contruction In the lat lecture, we have learned about the Lawon ip algorithm that compute a Delaunay triangulation of a given n-point et P R 2 with O(n 2 )
More informationDAROS: Distributed User-Server Assignment And Replication For Online Social Networking Applications
DAROS: Ditributed Uer-Server Aignment And Replication For Online Social Networking Application Thuan Duong-Ba School of EECS Oregon State Univerity Corvalli, OR 97330, USA Email: duongba@eec.oregontate.edu
More informationKS3 Maths Assessment Objectives
KS3 Math Aement Objective Tranition Stage 9 Ratio & Proportion Probabilit y & Statitic Appreciate the infinite nature of the et of integer, real and rational number Can interpret fraction and percentage
More informationFocused Video Estimation from Defocused Video Sequences
Focued Video Etimation from Defocued Video Sequence Junlan Yang a, Dan Schonfeld a and Magdi Mohamed b a Multimedia Communication Lab, ECE Dept., Univerity of Illinoi, Chicago, IL b Phyical Realization
More informationCutting Stock by Iterated Matching. Andreas Fritsch, Oliver Vornberger. University of Osnabruck. D Osnabruck.
Cutting Stock by Iterated Matching Andrea Fritch, Oliver Vornberger Univerity of Onabruck Dept of Math/Computer Science D-4909 Onabruck andy@informatikuni-onabrueckde Abtract The combinatorial optimization
More informationMulti-Target Tracking In Clutter
Multi-Target Tracking In Clutter John N. Sander-Reed, Mary Jo Duncan, W.B. Boucher, W. Michael Dimmler, Shawn O Keefe ABSTRACT A high frame rate (0 Hz), multi-target, video tracker ha been developed and
More information/06/$ IEEE 364
006 IEEE International ympoium on ignal Proceing and Information Technology oie Variance Etimation In ignal Proceing David Makovoz IPAC, California Intitute of Technology, MC-0, Paadena, CA, 95 davidm@ipac.caltech.edu;
More informationTexture-Constrained Active Shape Models
107 Texture-Contrained Active Shape Model Shuicheng Yan, Ce Liu Stan Z. Li Hongjiang Zhang Heung-Yeung Shum Qianheng Cheng Microoft Reearch Aia, Beijing Sigma Center, Beijing 100080, China Dept. of Info.
More informationTarget Oriented High Resolution SAR Image Formation via Semantic Information Guided Regularizations
1 Target Oriented High Reolution SAR Image Formation via Semantic Information Guided Regularization Biao Hou, Member, IEEE, Zaidao Wen, Licheng Jiao, Senior Member, IEEE, and Qian Wu arxiv:1704.07082v1
More informationA TOPSIS based Method for Gene Selection for Cancer Classification
Volume 67 No17, April 2013 A TOPSIS baed Method for Gene Selection for Cancer Claification IMAbd-El Fattah,WIKhedr, KMSallam, 1 Department of Statitic, 3 Department of Deciion upport, 2 Department of information
More informationHassan Ghaziri AUB, OSB Beirut, Lebanon Key words Competitive self-organizing maps, Meta-heuristics, Vehicle routing problem,
COMPETITIVE PROBABIISTIC SEF-ORGANIZING MAPS FOR ROUTING PROBEMS Haan Ghaziri AUB, OSB Beirut, ebanon ghaziri@aub.edu.lb Abtract In thi paper, we have applied the concept of the elf-organizing map (SOM)
More informationCERIAS Tech Report EFFICIENT PARALLEL ALGORITHMS FOR PLANAR st-graphs. by Mikhail J. Atallah, Danny Z. Chen, and Ovidiu Daescu
CERIAS Tech Report 2003-15 EFFICIENT PARALLEL ALGORITHMS FOR PLANAR t-graphs by Mikhail J. Atallah, Danny Z. Chen, and Ovidiu Daecu Center for Education and Reearch in Information Aurance and Security,
More informationCENTER-POINT MODEL OF DEFORMABLE SURFACE
CENTER-POINT MODEL OF DEFORMABLE SURFACE Piotr M. Szczypinki Iintitute of Electronic, Technical Univerity of Lodz, Poland Abtract: Key word: Center-point model of deformable urface for egmentation of 3D
More informationImage Morphing Using Deformation Techniques
Image Morphing Uing Deformation Technique Seung-Yong Lee, Kyung-Yong Chwa, Jame Hahn, and Sung Yong Shin Department of Computer Science Korea Advanced Intitute of Science and Technology 373-1 Kuong-dong
More informationChapter S:II (continued)
Chapter S:II (continued) II. Baic Search Algorithm Sytematic Search Graph Theory Baic State Space Search Depth-Firt Search Backtracking Breadth-Firt Search Uniform-Cot Search AND-OR Graph Baic Depth-Firt
More informationLocalized Minimum Spanning Tree Based Multicast Routing with Energy-Efficient Guaranteed Delivery in Ad Hoc and Sensor Networks
Localized Minimum Spanning Tree Baed Multicat Routing with Energy-Efficient Guaranteed Delivery in Ad Hoc and Senor Network Hanne Frey Univerity of Paderborn D-3398 Paderborn hanne.frey@uni-paderborn.de
More informationSee chapter 8 in the textbook. Dr Muhammad Al Salamah, Industrial Engineering, KFUPM
Goal programming Objective of the topic: Indentify indutrial baed ituation where two or more objective function are required. Write a multi objective function model dla a goal LP Ue weighting um and preemptive
More informationNew Structural Decomposition Techniques for Constraint Satisfaction Problems
113 New Structural Decompoition Technique for Contraint Satifaction Problem Yaling Zheng and Berthe Y. Choueiry Contraint Sytem Laboratory, Univerity of Nebraka-Lincoln {yzheng,choueiry}@ce.unl.edu Abtract.
More informationResearch Article Unsupervised SAR Image Segmentation Based on a Hierarchical TMF Model in the Discrete Wavelet Domain for Sea Area Detection
Dicrete Dynamic in Nature and Society, Article ID 354704, 9 page http://dx.doi.org/10.1155/2014/354704 Reearch Article Unupervied SAR Image Segmentation Baed on a Hierarchical TMF Model in the Dicrete
More informationLaboratory Exercise 6
Laboratory Exercie 6 Adder, Subtractor, and Multiplier a a The purpoe of thi exercie i to examine arithmetic circuit that add, ubtract, and multiply number. Each b c circuit will be decribed in Verilog
More informationShortest Paths Problem. CS 362, Lecture 20. Today s Outline. Negative Weights
Shortet Path Problem CS 6, Lecture Jared Saia Univerity of New Mexico Another intereting problem for graph i that of finding hortet path Aume we are given a weighted directed graph G = (V, E) with two
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
Thi article appeared in a journal publihed by Elevier. The attached copy i furnihed to the author for internal non-commercial reearch and education ue, including for intruction at the author intitution
More informationQuadrilaterals. Learning Objectives. Pre-Activity
Section 3.4 Pre-Activity Preparation Quadrilateral Intereting geometric hape and pattern are all around u when we tart looking for them. Examine a row of fencing or the tiling deign at the wimming pool.
More informationTopics. Lecture 37: Global Optimization. Issues. A Simple Example: Copy Propagation X := 3 B > 0 Y := 0 X := 4 Y := Z + W A := 2 * 3X
Lecture 37: Global Optimization [Adapted from note by R. Bodik and G. Necula] Topic Global optimization refer to program optimization that encompa multiple baic block in a function. (I have ued the term
More informationService and Network Management Interworking in Future Wireless Systems
Service and Network Management Interworking in Future Wirele Sytem V. Tountopoulo V. Stavroulaki P. Demeticha N. Mitrou and M. Theologou National Technical Univerity of Athen Department of Electrical Engineering
More informationShortest-Path Routing in Arbitrary Networks
Ž. Journal of Algorithm 31, 105131 1999 Article ID jagm.1998.0980, available online at http:www.idealibrary.com on Shortet-Path Routing in Arbitrary Network Friedhelm Meyer auf der Heide and Berthold Vocking
More informationCSE 250B Assignment 4 Report
CSE 250B Aignment 4 Report March 24, 2012 Yuncong Chen yuncong@c.ucd.edu Pengfei Chen pec008@ucd.edu Yang Liu yal060@c.ucd.edu Abtract In thi project, we implemented the recurive autoencoder (RAE) a decribed
More informationModeling and Analysis of Slow CW Decrease for IEEE WLAN
Modeling and Analyi of Slow CW Decreae for IEEE 82. WLAN Qiang Ni, Imad Aad 2, Chadi Barakat, and Thierry Turletti Planete Group 2 Planete Group INRIA Sophia Antipoli INRIA Rhône-Alpe Sophia Antipoli,
More informationAn Intro to LP and the Simplex Algorithm. Primal Simplex
An Intro to LP and the Simplex Algorithm Primal Simplex Linear programming i contrained minimization of a linear objective over a olution pace defined by linear contraint: min cx Ax b l x u A i an m n
More informationToday s Outline. CS 561, Lecture 23. Negative Weights. Shortest Paths Problem. The presence of a negative cycle might mean that there is
Today Outline CS 56, Lecture Jared Saia Univerity of New Mexico The path that can be trodden i not the enduring and unchanging Path. The name that can be named i not the enduring and unchanging Name. -
More informationThe Data Locality of Work Stealing
The Data Locality of Work Stealing Umut A. Acar School of Computer Science Carnegie Mellon Univerity umut@c.cmu.edu Guy E. Blelloch School of Computer Science Carnegie Mellon Univerity guyb@c.cmu.edu Robert
More informationDigital Image Processing Laboratory: Markov Random Fields and MAP Image Segmentation
Purdue University: Digital Image Processing Laboratories Digital Image Processing Laboratory: Markov Random Fields and MAP Image Segmentation December, 205 Introduction This laboratory explores the use
More informationANALYSIS OF THE FIRST LAYER IN WEIGHTLESS NEURAL NETWORKS FOR 3_DIMENSIONAL PATTERN RECOGNITION
ANALYSIS OF THE FIRST LAYER IN WEIGHTLESS NEURAL NETWORKS FOR 3_DIMENSIONAL PATTERN RECOGNITION A. Váque-Nava * Ecuela de Ingeniería. CENTRO UNIVERSITARIO MEXICO. DIVISION DE ESTUDIOS SUPERIORES J. Figueroa
More informationMulticlass Road Sign Detection using Multiplicative Kernel
Proceeding of the Croatian Computer Viion Workhop, Year 1 Multicla Road Sign Detection uing Multiplicative Kernel Valentina Zadrija Mireo d. d. Zagreb, Croatia valentina.zadrija@mireo.hr Siniša Šegvić
More informationComputer Arithmetic Homework Solutions. 1 An adder for graphics. 2 Partitioned adder. 3 HDL implementation of a partitioned adder
Computer Arithmetic Homework 3 2016 2017 Solution 1 An adder for graphic In a normal ripple carry addition of two poitive number, the carry i the ignal for a reult exceeding the maximum. We ue thi ignal
More informationIMPLEMENTATION OF AREA, VOLUME AND LINE SOURCES
December 01 ADMS 5 P503I1 IMPEMENTATION OF AREA, VOUME AND INE SOURCES The Met. Office (D J Thomon) and CERC 1. INTRODUCTION ADMS model line ource, and area and volume ource with conve polgon bae area.
More informationPlanning of scooping position and approach path for loading operation by wheel loader
22 nd International Sympoium on Automation and Robotic in Contruction ISARC 25 - September 11-14, 25, Ferrara (Italy) 1 Planning of cooping poition and approach path for loading operation by wheel loader
More informationA Hybrid Deployable Dynamic Traffic Assignment Framework for Robust Online Route Guidance
A Hybrid Deployable Dynamic Traffic Aignment Framework for Robut Online Route Guidance Sriniva Peeta School of Civil Engineering, Purdue Univerity Chao Zhou Sabre, Inc. Sriniva Peeta School of Civil Engineering
More informationarxiv: v1 [cs.ds] 27 Feb 2018
Incremental Strong Connectivity and 2-Connectivity in Directed Graph Louka Georgiadi 1, Giueppe F. Italiano 2, and Niko Parotidi 2 arxiv:1802.10189v1 [c.ds] 27 Feb 2018 1 Univerity of Ioannina, Greece.
More informationEVen though remarkable progress has been made in
JOURNAL OF L A T E X CLASS FILES, VOL. 14, NO. 8, AUGUST 2015 1 Dene Semantic 3D Recontruction Chritian Häne, Chritopher Zach, Andrea Cohen, and Marc Pollefey, Fellow, IEEE Abtract Both image egmentation
More informationLaboratory Exercise 6
Laboratory Exercie 6 Adder, Subtractor, and Multiplier The purpoe of thi exercie i to examine arithmetic circuit that add, ubtract, and multiply number. Each circuit will be decribed in Verilog and implemented
More informationA Fast Association Rule Algorithm Based On Bitmap and Granular Computing
A Fat Aociation Rule Algorithm Baed On Bitmap and Granular Computing T.Y.Lin Xiaohua Hu Eric Louie Dept. of Computer Science College of Information Science IBM Almaden Reearch Center San Joe State Univerity
More informationA SIMPLE IMPERATIVE LANGUAGE THE STORE FUNCTION NON-TERMINATING COMMANDS
A SIMPLE IMPERATIVE LANGUAGE Eventually we will preent the emantic of a full-blown language, with declaration, type and looping. However, there are many complication, o we will build up lowly. Our firt
More informationIncreasing Throughput and Reducing Delay in Wireless Sensor Networks Using Interference Alignment
Int. J. Communication, Network and Sytem Science, 0, 5, 90-97 http://dx.doi.org/0.436/ijcn.0.50 Publihed Online February 0 (http://www.scirp.org/journal/ijcn) Increaing Throughput and Reducing Delay in
More informationUsing Mouse Feedback in Computer Assisted Transcription of Handwritten Text Images
2009 10th International Conference on Document Analyi and Recognition Uing Moue Feedback in Computer Aited Trancription of Handwritten Text Image Verónica Romero, Alejandro H. Toelli, Enrique Vidal Intituto
More informationMotion Control (wheeled robots)
3 Motion Control (wheeled robot) Requirement for Motion Control Kinematic / dynamic model of the robot Model of the interaction between the wheel and the ground Definition of required motion -> peed control,
More informationTouring a Sequence of Polygons
Touring a Sequence of Polygon Mohe Dror (1) Alon Efrat (1) Anna Lubiw (2) Joe Mitchell (3) (1) Univerity of Arizona (2) Univerity of Waterloo (3) Stony Brook Univerity Problem: Given a equence of k polygon
More informationSegmentation of ultrasound images multiresolution 2D and 3D algorithm based on global and local statistics
Pattern Recognition Letter 24 (2003) 779 790 www.elevier.com/locate/patrec Segmentation of ultraound image multireolution 2D and 3D algorithm baed on global and local tatitic Djamal Boukerroui a,b, *,
More informationA Practical Model for Minimizing Waiting Time in a Transit Network
A Practical Model for Minimizing Waiting Time in a Tranit Network Leila Dianat, MASc, Department of Civil Engineering, Sharif Univerity of Technology, Tehran, Iran Youef Shafahi, Ph.D. Aociate Profeor,
More informationKey Terms - MinMin, MaxMin, Sufferage, Task Scheduling, Standard Deviation, Load Balancing.
Volume 3, Iue 11, November 2013 ISSN: 2277 128X International Journal of Advanced Reearch in Computer Science and Software Engineering Reearch Paper Available online at: www.ijarce.com Tak Aignment in
More informationAdvanced Encryption Standard and Modes of Operation
Advanced Encryption Standard and Mode of Operation G. Bertoni L. Breveglieri Foundation of Cryptography - AES pp. 1 / 50 AES Advanced Encryption Standard (AES) i a ymmetric cryptographic algorithm AES
More informationThe Set Constraint/CFL Reachability Connection in Practice
The Set Contraint/CFL Reachability Connection in Practice John Kodumal EECS Department Univerity of California, Berkeley jkodumal@c.berkeley.edu Alex Aiken Computer Science Department Stanford Univerity
More informationDifferentially Private Histogram Publication
Differentially Private Hitogram Publication Jia Xu, Zhenjie Zhang, Xiaokui Xiao, Yin Yang, Ge Yu College of Information Science & Engineering, Northeatern Univerity, China xujia,yuge@ie.neu.edu.cn Advanced
More informationLaboratory Exercise 6
Laboratory Exercie 6 Adder, Subtractor, and Multiplier The purpoe of thi exercie i to examine arithmetic circuit that add, ubtract, and multiply number. Each circuit will be decribed in VHL and implemented
More informationParallel MATLAB at FSU: Task Computing
Parallel MATLAB at FSU: Tak John Burkardt Department of Scientific Florida State Univerity... 1:30-2:30 Thurday, 07 April 2011 499 Dirac Science Library... http://people.c.fu.edu/ jburkardt/preentation/...
More informationDrawing Lines in 2 Dimensions
Drawing Line in 2 Dimenion Drawing a traight line (or an arc) between two end point when one i limited to dicrete pixel require a bit of thought. Conider the following line uperimpoed on a 2 dimenional
More informationM 4 CD: A Robust Change Detection Method for Intelligent Visual Surveillance
1 M 4 CD: A Robut Change Detection Method for Intelligent Viual Surveillance Kunfeng Wang, Member, IEEE, Chao Gou, and Fei-Yue Wang, Fellow, IEEE Abtract In thi paper, we propoe a robut change detection
More informationA System Dynamics Model for Transient Availability Modeling of Repairable Redundant Systems
International Journal of Performability Engineering Vol., No. 3, May 05, pp. 03-. RAMS Conultant Printed in India A Sytem Dynamic Model for Tranient Availability Modeling of Repairable Redundant Sytem
More informationPerformance Evaluation of an Advanced Local Search Evolutionary Algorithm
Anne Auger and Nikolau Hanen Performance Evaluation of an Advanced Local Search Evolutionary Algorithm Proceeding of the IEEE Congre on Evolutionary Computation, CEC 2005 c IEEE Performance Evaluation
More informationImproved Inference in Bayesian Segmentation Using Monte Carlo Sampling: Application to Hippocampal Subfield Volumetry
Improved Inference in Bayeian Segmentation Uing Monte Carlo Sampling: Application to Hippocampal Subfield Volumetry Juan Eugenio Igleia a, Mert Rory Sabuncu a, Koen Van Leemput a,b,c, for the Alzheimer
More informationImplementation of a momentum-based distance metric for motion graphs. Student: Alessandro Di Domenico (st.no ), Supervisor: Nicolas Pronost
Implementation of a momentum-baed ditance metric for motion graph Student: Aleandro Di Domenico (t.no 3775682), Supervior: Nicola Pronot April 3, 2014 Abtract Thi report preent the procedure and reult
More informationLaboratory Exercise 2
Laoratory Exercie Numer and Diplay Thi i an exercie in deigning cominational circuit that can perform inary-to-decimal numer converion and inary-coded-decimal (BCD) addition. Part I We wih to diplay on
More informationUsing Partial Evaluation in Distributed Query Evaluation
A X x Z z R r y Y B Uing Partial Evaluation in Ditributed Query Evaluation Peter Buneman Gao Cong Univerity of Edinburgh Wenfei Fan Univerity of Edinburgh & Bell Laboratorie Anataio Kementietidi Univerity
More informationLoad Estimation of Social Networking Web Sites Using Clustering Technique
International Journal of Electronic and Electrical Engineering Vol. 4, No. 6, December 216 Load Etimation of Social Networking Web Site Uing Clutering Technique Deepti Bhagwani and Setu Kumar Chaturvedi
More informationCompressed Sensing Image Processing Based on Stagewise Orthogonal Matching Pursuit
Senor & randucer, Vol. 8, Iue 0, October 204, pp. 34-40 Senor & randucer 204 by IFSA Publihing, S. L. http://www.enorportal.com Compreed Sening Image Proceing Baed on Stagewie Orthogonal Matching Puruit
More informationAN ALGORITHM FOR RESTRICTED NORMAL FORM TO SOLVE DUAL TYPE NON-CANONICAL LINEAR FRACTIONAL PROGRAMMING PROBLEM
RAC Univerity Journal, Vol IV, No, 7, pp 87-9 AN ALGORITHM FOR RESTRICTED NORMAL FORM TO SOLVE DUAL TYPE NON-CANONICAL LINEAR FRACTIONAL PROGRAMMING PROLEM Mozzem Hoain Department of Mathematic Ghior Govt
More informationA METHOD OF REAL-TIME NURBS INTERPOLATION WITH CONFINED CHORD ERROR FOR CNC SYSTEMS
Vietnam Journal of Science and Technology 55 (5) (017) 650-657 DOI: 10.1565/55-518/55/5/906 A METHOD OF REAL-TIME NURBS INTERPOLATION WITH CONFINED CHORD ERROR FOR CNC SYSTEMS Nguyen Huu Quang *, Banh
More informationA Linear Interpolation-Based Algorithm for Path Planning and Replanning on Girds *
Advance in Linear Algebra & Matrix Theory, 2012, 2, 20-24 http://dx.doi.org/10.4236/alamt.2012.22003 Publihed Online June 2012 (http://www.scirp.org/journal/alamt) A Linear Interpolation-Baed Algorithm
More information