Approximation of the Euclidean Distance by Chamfer Distances
|
|
- Hortense Newman
- 5 years ago
- Views:
Transcription
1 Acta Cybernetica 0 ( Aroximation of the Euclidean Distance by Chamfer Distances András Hajdu, Lajos Hajdu, and Robert Tijdeman Abstract Chamfer distances lay an imortant role in the theory of distance transforms. Though the determination of the exact Euclidean distance transform is also a well investigated area, the classical chamfering method based uon small neighborhoods still outerforms it e.g. in terms of comutation time. In this aer we determine the best ossible maximum relative error of chamfer distances under various boundary conditions. In each case some best aroximating sequences are exlicitly given. Further, because of ossible ractical interest, we give all best aroximating sequences in case of small (i.e. and 7 7 neighborhoods. Keywords: Chamfering, Aroximation of the Euclidean distance, Distance transform, Digital image rocessing Introduction Suose we measure distances between grid oints of a two-dimensional grid and we want to aroximate the Euclidean distance by a distance function which can be comuted quickly, without calculating square roots. We may then use the class of chamfer distances. They are obtained by rescribing the lengths of the grid vectors in a so-called mask M := {(x, y Z : max( x, y } (for some ositive integer Research of the Hungarian authors was suorted in art by the OTKA grants F043090, T0498, T04879, K6780, K766, NK0680, by the János Bolyai Research Fellowshi of the Hungarian Academy of Sciences, by the TECH08- roject DRSCREEN - Develoing a comuter based image rocessing system for diabetic retinoathy screening of the National Office for Research and Technology of Hungary (contract no.: OM-0094/008, OM-009/008, OM- 0096/008, and by the TÁMOP 4.../B-09//KONV roject, which is imlemented through the New Hungary Develoment Plan, cofinanced by the Euroean Social Fund and the Euroean Regional Develoment Fund. Faculty of Informatics, University of Debrecen, H-400 Debrecen, P.O. Box, Hungary. hajdu.andras@inf.unideb.hu Number Theory Research Grou of the Hungarian Academy of Sciences and Institute of Mathematics, University of Debrecen, H-400 Debrecen, P.O. Box, Hungary. hajdul@math.unideb.hu Mathematical Institute, Leiden University. Niels Bohrweg, Leiden Postbus 9, 300 RA Leiden, The Netherlands. tijdeman@math.leidenuniv.nl
2 400 András Hajdu, Lajos Hajdu, and Robert Tijdeman such that the values at (±x, ±y and (±y, ±x are all the same, and by defining the length function W as follows: the length W ( v of any vector v Z is defined as the minimal sum of the lengths of those vectors from M, reetitions ermitted, which have sum v. The literature on chamfer distances is very rich. See Borgefors [, 3, 4] for the basics, [7, 8] for lists of ( + ( + neighborhoods for 0, and [7] for an overview of alications. Further, recently many related results have been obtained by several authors, concerning distance transforms and their exlicit calculation using different kinds of neighborhoods in certain (mostly 3D grids. For examle, Strand, Nagy, Fouard and Borgefors [0] gave a sequential algorithm for comuting the distance ma using distances based on neighborhood sequences in the D square grid, and 3D cubic and so-called FCC and BCC cubic grids, resectively. Similar results for other kinds of grids are also known, see e.g. [6] (nd hexagonal grids, [] (diamond grid and [] (general oint grids and the references given there. Classical chamfer distances using 3 3, and 7 7 neighborhoods given by Borgefors [, 3] are generated by the masks , and resectively (with the actual generator entries underlined. For comarison with the Euclidean distance the values of the neighborhoods have to be divided by 3, and, resectively. The aroximations to.4 are therefore 4/3.33, 7/ =.4 and 7/.4, resectively. For alternative neighborhood values see Verwer [, 3], Thiel [], Coquin and Bolon [6], Butt and Maragos [] and Scholtus [7]. More secifically, in [6] the minimization of the error between the Euclidean distance and the local distance was considered over circular trajectories similarly to [, 3] rather than linear ones [3, ]. The aroximation error can also be measured based on area as it is done in [] with calculating the difference between a disk of large size obtained by chamfer metric and a Euclidean disk of the same radius. The determination of the exact Euclidean distance transform is also a well investigated area (see e.g. [, 7, 8, 3, 9], but the classical 3 3 chamfering method still outerforms it in terms of comutation time and simle extendability to other grids. In this aer we determine chamfer distances best aroximating the Euclidean distance in a certain sense. In each neighborhood size some best aroximating sequences are exlicitly given. Further, because of ossible ractical interest, we give all best aroximating sequences in case of small (i.e. and 7 7 neighborhoods. Throughout the aer, as a measure for the quality of a length function W
3 Aroximation of the Euclidean Distance by Chamfer Distances 40 defined on Z we use the so-called maximum relative error (m.r.error for short E := lim su W ( v v v where. denotes the Euclidean length. The M -, M - and M 3 -neighborhoods given above yield rounded E-values 0.07, and 0.038, resectively. Firstly we shall rove that the smallest ossible constant E B for the mask M under the condition that W (x, 0 = x for x Z is given by E B = =. ( + O 4. In articular, E B 0.0, E B and E3 B Comaring these values with the E-values given above, one can see that the E B -values yield aroximately 4%, 6% and 3% imrovement, resectively. The B refers to Borgefors who was the first to consider such neighborhoods. Secondly we consider the case D in which W ( v v for all v Z. (The D refers to the fact that W ( v dominates v. The otimal m.r.error under this restriction equals E D = ( + + = ( 8 + O 4 = 0. ( + O 4. In articular, E D 0.084, E D 0.07 and E3 D Thirdly we shall rove that the otimal E-value without any restriction on the neighborhood defined on M (i.e. droing the condition W (x, 0 = x for x Z equals + + E C = + = ( O 4. In articular, E C , E C and E3 C In 99, on using the symmetry in case C the value of E C was comuted by Verwer [, 3] in terms of trigonometric functions. The C refers to the word central. In 998, because of geometric considerations, Butt and Maragos [] chose to use the error function lim su v v W ( v which of course is small if and only if E C is small. In general it gives different error values, but the values for E C are equal to the values obtained by the above error function (cf. Scholtus [7]. We rove the correctness of the above E C values. In doing so, our motivation is twofold: on the one hand, by a simle reasoning we obtain these values immediately from the values of E D, and on the other hand, our
4 40 András Hajdu, Lajos Hajdu, and Robert Tijdeman roofs are mathematically rigorous while the corresonding arguments of Verwer and Butt and Maragos contain some hidden assumtions. Namely, by certain lausible but not exlicitly verified geometric arguments they restrict their attention and investigations to certain values of the neighborhoods in question, and they erform exact investigations only for these values. We shall further study an auxiliary class of neighborhoods on M, viz. the class of neighborhoods satisfying N c ( v = for all v = (x, y M with either x < or y < 0, N c ( v = for v = (, 0, and N c ( v = c v for v = (, k with 0 < k. Here c is a constant close to and at most equal to. Informally seaking, the use of such neighborhoods means that only such stes (v, v are allowed, where v is a ositive multile of and v is nonnegative. Further, beside N c (, 0 = the weights of the other such neighborhood vectors are their Euclidean lengths, multilied by a factor c. All the other vectors of the neighborhood are forbidden to use, thus they have weights. For examle, the weights for the neighborhood N c with = (i.e. for M are given by c 8 c where the origin is in the middle. We denote the maximum relative error for this class of neighborhoods by E c where we restrict the limsu to vectors v with finite lengths W ( v (i.e. having coordiantes (x, y with 0 y x and x. Our motivation for considering such neighborhoods is that it will turn out that (due to its secial form N c is easier to handle, but yields the same m.r.error as the corresonding neighborhood N c, in which N c (±, 0 = N c (0, ± = and N c (x, y = c x + y otherwise ((x, y M. In Section we introduce some notation and rove some reliminary results. In Sections 3 and 4 we comute the values of E B and E D where E B is the maximum relative error E c for otimal c and E D = E. We give all sequences yielding minimal m.r.error in case of and 7 7 neighborhoods, as well. In Section 4 we rove that E B = E B and E D = E D and further show that E C = E D /( + E D for all. Finally, we draw some conclusions in Section. Definitions and basic roerties Let N be a neighborhood defined on the mask M. Put M = M \ {(0, 0}. We denote the value of N at osition (n, k by w(n, k for (n, k M. Throughout the aer we assume that w(±n, ±k = w(±k, ±n > 0 for all (n, k M and all ossible sign choices. Hence it suffices to consider the values w(n, k with 0 k n. We can measure lengths of vectors and distances between oints using neighborhood sequences. Note that such sequences rovide a flexible and very useful
5 Aroximation of the Euclidean Distance by Chamfer Distances 403 tool in handling several roblems in discrete geometry. For the basics and most imortant facts about such sequences, see e.g. the aers [9, 4, 0,, 4] and the references given there. Here we only give those notions which will be needed for our uroses. Let A = (N i i= be a sequence of neighborhoods defined on M and u, v Z. The sequence u = u 0, u,..., u m = v with u i u i M is called an A-ath from u to v. The A-length of the ath is defined as m w i ( u i u i. The distance W A ( v u between u and v, which is the A-length of v u, is defined as the minimal A-length taken over all A-aths from u to v. If the neighborhood sequence is fixed, then we suress the letter A in the above notation. If N i = N for all i, then the corresonding (constant neighborhood sequence is denoted by A = N. We assume throughout the aer that for such sequences W (n, k = w(n, k holds for (n, k M ; if it would not have been the case, then the function w := W M would have generated W, too. We call W a metric if for all u, v Z W ( u < (W is finite, W ( u = 0 u = 0 (W is ositive definite, W ( u = W ( u (W is symmetric, W ( u + v W ( u + W ( v (W satisfies the triangle inequality. It follows from the above roerties that W ( u 0 for every u Z. By our basic assumtions on w, every induced length function W is ositive definite and symmetric. Furthermore, W satisfies the triangle inequality for u, v with u, v, u+ v M by definition. The first lemma shows that for a constant neighborhood sequence W ( v/ v attains a minimal value which is reached already in M. Lemma. Let N be a neighborhood defined on M which induces the length function W on Z. Then W ( v w( v lim inf = min. v v v M v w( v Proof. Let m = min v M v = w( u u W ( v so that lim inf v v i= ( u M. Then for all n we have W (n u n u = m, m. On the other hand, since w( v v m for every v M, it follows from the definition of shortest ath and the triangle inequality for the Euclidean distance that W ( v i w( v i = i w( v i v i v i m i v i m v for every v Z W ( v not equal to the origin. Thus lim inf v v m.
6 404 András Hajdu, Lajos Hajdu, and Robert Tijdeman The challenge is therefore to comute lim su v W ( v v. 3 The maximum relative error for neighborhoods N c Let c be some ositive real number with < c. We shall study neighborhoods N c on M with N c (n, k = for which either n < or k < 0, N c (, 0 = + and N c (, k = c + k for 0 < k. We are interested in the length function W c induced by A c := N c for oints in the set {(x, y Z : x, 0 y x}. First we secure that under suitable conditions only two distinct stes occur in a shortest A c -ath. Lemma. Let < c. Then a shortest A c-ath from (0, 0 to (m, mr+k + with m, r, k Z, 0 r <, 0 k < m consists only of stes (, r and (, r +. Proof. Suose a shortest ath from (0, 0 to (m, mr + k with m, r, k Z, 0 r <, 0 k < m contains two stes (, t and (, u with t u 0. Relace the two stes with stes (, t and (, u +, and write L and L for the length of the old and new aths, resectively. Then we have L L c + t c + (t + c + u c + (u + = = c(f (t f (u +, where f (x = + x + (x (x Z 0. A simle calculation yields that f (x is strictly monotone increasing in x, which shows that L L > 0. However, this contradicts the minimality of the length of the original ath. Hence a shortest ath may contain stes (, t and (, t + only, for some nonnegative integer t. Since altogether we make m stes, this immediately gives that t = r, and our statement follows. Remark. The latter inequality is the most severe and exlains why we restrict c to values greater than. + Corollary. Let < c Then a shortest A c-ath from (0, 0 to (m, mr + with 0 r consists of m stes (, r. The next theorem gives the value of the aroximation error for general, in case of any neighborhood N c on M.
7 Aroximation of the Euclidean Distance by Chamfer Distances 40 Theorem. Let, distance is given by max( c, < c. Then the m.r.error of A c to the Euclidean + + c + + c c +. Proof. As a general remark we mention that to erform our calculations, we used the rogram ackage Male R. Let be a ositive integer, and fix c with < c. As reviously, it + is sufficient to consider the A c -length of oints of the form (m, k where m is some ositive integer and k is an integer with 0 k m. Write k = mq + r with 0 q and 0 r < m. The ossible stes are (, 0 of length and (, ±i of length W i := c + i (for i. From Lemma and the inequalities = W 0 < W <... < W we see that a ath of minimal length from (0, 0 to a oint (m, mq + r consists of r stes (, q + and m r stes (, q. Hence for the induced length function we get W(m, mq + r = rw q+ + (m rw q. Put t = r/m, and recall that W 0 = and W i = c + i for i =,...,. Set and for q and let H 0 (t = ct + + ( t + t, H q (t = c t + (q + + ( t + q + (q + t, h q (, c = max 0 t H q(t (0 q < and h (, c = H (0. Now we investigate the error functions h q (, c for q =, q = 0, 0 < q <, resectively. Suose first that q =. Then r = 0 and k = m. In this case we trivially have h (, c = c. Assume next that q = 0. Then 0 k <. Put t 0 := (c +. A simle calculation yields that 0 t 0, and that H 0 is monotone increasing on the interval [0, t 0 ] and monotone decreasing on the interval [t 0, ]. Moreover, we have H 0 (0 = 0 and H 0 ( = c, hence H 0 (t 0 0. Thus we have h 0 (, c = max( c, H 0 (t 0 = max( c, + c + + c c +. Male is a registered trademark of Waterloo Male Inc.
8 406 András Hajdu, Lajos Hajdu, and Robert Tijdeman Finally, suose that 0 < q <, that is k < m. Put + q t q := ( ( + q ( + (q + q q (q + + q q. + (q + A simle calculation gives that 0 t q, and that H q is monotone increasing on the interval [0, t q ], while monotone decreasing on the interval [t q, ]. We also have H q (0 = H q ( = c. Hence H q (t q < 0 imlies H q (t q c. Thus we get h q (, c = max( c, H q (t q = max c, c + ( +q ( +(q+ Now we calculate the error function h(, c := lim su W(n, k n + k = max h q(, c. 0 q n, n k 0. Observe first that for fixed and c the function h q (, c is monotone decreasing in q with q. Hence h q (, c h (, c for q =,...,. Further, again by Male, we obtain that for any c with < c + c + ( +( +4 + c + + c c + holds, which imlies h (, c h 0 (, c. Hence h(, c = max( c, + c + + c c + and the theorem follows. The following corollaries rovide the m.r.errors E B (when c = c B and E D (when c =, resectively. Corollary. Let be a ositive integer. Then we have c B = That is, the sequence A = A c B of eriod given by A = N c B yields the smallest m.r.error among all sequences A c of eriod. Moreover, the error is given by E B = c B = =. ( + O ( + O 4.
9 Aroximation of the Euclidean Distance by Chamfer Distances 407 Proof. Put f(c = c and g(c = + c + + c c +. A straightforward comutation shows that f is strictly monotone decreasing, while g is strictly monotone increasing for < c. Hence there is a unique + solution of the equation f(c = g(c in this interval. By Theorem this solution is given by c B = Thus the statement follows. Corollary 3. Let be a ositive integer. Then the sequence A = A of eriod given by A = N (corresonding to the choice c = has m.r.error E D = ( + + = ( 8 + O 4 = 0. ( + O 4. Proof. On substituting c = into the formula of Theorem, the statement follows immediately. Now we give the best aroximating sequences realizing the minimal maximum relative error for matrices ( = in Theorem and for 7 7 matrices ( = 3 in Theorem 3, resectively. Theorem. Let < c. Let A c = N c be the corresonding sequence on M. Then the minimal m.r.error to the Euclidean distance among the neighborhood sequences A c is attained if and only if where c = c B, W = s and u W v, s = , u = s.776 and v = + s Further, the m.r.error is given by E B = c B = s =
10 408 András Hajdu, Lajos Hajdu, and Robert Tijdeman Proof. For any even n with 0 k n the ossible stes are (, 0 of length, (, and (, of length W, and (, and (, of length W. From Lemma and the inequality < W < W we see that the ath from (0, 0 to (n, k of minimal length consists of k stes (, and n k stes (, 0 if 0 k n/ and of k n/ stes (, and n k stes (, if n k n. Hence we have for the induced length function { kw + n k, if k n W(n, k =, (n kw + (k n W, otherwise. Put t = k/n. Then the error function is given by h(w, W := lim su W(n, k n + k = max ( max 0 t t(w + + t n, n k 0, max t ( tw + (t W + t. Our aim is to choose W and W such that h(w, W is minimal. For fixed W, define the function H 0 : R 0 R by H 0 (t = t(w +. + t Put t 0 = W. We observe that H 0 is monotone increasing on [0, t 0 ] and monotone decreasing on [t 0,. Hence, as H 0 (0 =, ( ( max ( H 0 (t = max H 0 (t 0, H 0 0 t ( = max W 4W +, W if W / and ( max 0 (t = H 0 = 0 t ( H W otherwise. Clearly, ( min (h(w, W min max W 4W +, W,W W. ( W A calculation gives that the minimum of the right-hand side is achieved for W = s :=
11 Aroximation of the Euclidean Distance by Chamfer Distances 409 and equals s 4s + = s = Now we fix the value s of W, and show that we can choose W in a way to have equality in (. In fact we comletely describe the set of the aroriate W -s. Consider the maximum over t [/, ]. For fixed W, define the function H : R 0 R by H (t = ( tw + (t W. + t Observe that H attains its maximum at t := (W W W W (which is ositive and further, H is monotone increasing in [0, t ] and monotone decreasing in [t,. Hence ( ( max ( H (t = max H, H (t, H ( = t ( = max W (W W, + 4(W W, W if / t, and max t ( H (t = max ( ( H, H ( = ( = max W, W otherwise. By our choice of W, we have that W = s The values of W (W W and +4(W W do not exceed this value if and only if u W v where u and v are defined in the statement of the theorem. We conclude that h(w, W attains its minimum s if W = s and u W v. The above argument shows that E B = W. Hence the minimum among neighborhoods N c is realized for c = c B = W and for no other value of c. 3 Theorem 3. Let 0 < c. Let A c = N c be the corresonding sequence on M 3. Then the minimal m.r.error to the Euclidean distance among the neighborhood sequences A c is attained if and only if c = c B 3, W = s, u W v, q W 3 r,
12 40 András Hajdu, Lajos Hajdu, and Robert Tijdeman where s = , 9 3 u = s 3.733, 0 v = 43s s q = 3s 4.047, r = 3 3s 0s + 0 0W , + W 3, and in the definition of r, W can be any number with u W v. Further, the m.r.error is given by E3 B = c B 3 = s = Proof. Let 3 n and 0 k n. The ossible stes are (3, 0 of length 3, (3, ± of length W, (3, ± of length W, and (3, ±3 of length W 3. From the inequalities 3 < W < W < W 3 it follows that the ath from (0, 0 to (n, k of minimal length consists of k stes (3, and n 3 k stes (3, 0 if 0 k n 3 ; of k n 3 stes (3, and n 3 k stes (3, if n 3 k n 3 ; of k n 3 stes (3, 3 and n k stes (3, if n 3 k n. Hence we have for the induced length function kw + n 3k, if k n/3, W(n, k = (n/3 kw + (k n/3w, if n/3 < k n/3, (n kw + (k n/3w 3, otherwise. Put t = k/n, and define the functions H i : R 0 R (i = 0,, by H 0 (t = t(w ( 3 + 3, H (t = t W + ( t 3 + t + t W and H (t = ( tw + ( t 3 + t W3. Then for fixed W, W, W 3 the error of aroximation is given by ( h(w, W, W 3 = max max 0 t 3 H 0 (t, max 3 t 3 H (t, max 3 t H (t. Let t 0 = W 3, t = 3(W W W W, t = 3(W 3 W 3W W 3,
13 Aroximation of the Euclidean Distance by Chamfer Distances 4 and observe that all t 0, t and t are ositive. By differentiation and following standard calculus, we get that for i = 0,,, H i is monotone decreasing if t i [i/3, (i + /3], and that H i is monotone increasing in [i/3, t i ] and monotone decreasing in [t i, (i + /3] otherwise. Hence from H 0 (0 = we get that ( ( max ( H 0 (t = max H 0 (t 0, H 0 0 t 3 = 3 ( = max W 6W + 0, W 0. Hence obviously, ( min h(w, W, W 3 min max W 6W + 0, W,W,W 3 W. ( 0 W By a simle calculation we get that the minimum of the right-hand side is achieved for and equals W = s := M := s 6s + 0 = s = Now we fix the value s of W, and show that we can choose W and W 3 in a way to have equality in (. More recisely, we comletely describe the set of the aroriate airs (W, W 3. For this urose, first we consider the maximum of H over t [/3, /3]. In a similar manner as in the roof of Theorem, we obtain that max ( H (t = max 3 t 3 = max W 0, ( ( H 3, H (t, H (W W + 9(W W 3 ( 3 =, W 3. Using our choice for W, a simle calculation gives that the above maximum does not exceed the value of M recisely when u W v, where u and v are defined in the statement of the theorem. So let W be any fixed number from the interval [u, v], and consider the the maximum of H over t [/3, ]. Now we get that ( ( max ( H (t = max H 3 t 3, H (t, H ( =
14 4 András Hajdu, Lajos Hajdu, and Robert Tijdeman (3W W 3 + 9(W 3 W ( = max W 3,, W Using our choice for W and W, a simle calculation yields that the above maximum is not larger than M if and only if q W 3 r, where q and r are given in the statement. (Note that < r < The above argument shows that E B 3 = W 0. Hence the minimum among neighborhoods N c is realized for c = c B 3 = W 0, and the theorem follows. 4 Equivalence of m.r.errors for M neighborhoods In this section we comute the m.r.errors E B, E C and E D. First we introduce neighborhoods N c on M defined by N c (0, 0 =, N c (n, 0 = N c (0, n = n for 0 < n, N c (n, k = c n + k for (n, k M, nk 0. Let W c denote the length function induced by the sequence N c. We show that the corresonding m.r.error E c satisfies E c = E c for every considered value of c. It then follows that E B = E B and E D = E D for every. Lemma 3. Let < c. There is a shortest N c-ath from (0, 0 to (m, k + with 0 k m which consists of stes of the form (, 0 and (,. Proof. Suose a shortest ath from (0, 0 to (m, k contains a ste (g, h with h < 0. Then it also contains a ste (i, j with j. But it is shorter to relace both stes with stes (g, h + and (i, j. A similar argument can be used to exclude stes (g, h with h >. So every shortest ath from (0, 0 to (m, k contains only stes of the forms (g, 0 and (g,. If k = m, then taking only stes (, gives the shortest ath length because of the triangle inequality for the Euclidean distance and the inequality c. Suose that there is a ste (g, with g < in a shortest ath from (0, 0 to (m, k with 0 k < m. Then there is also a ste (h, 0 with h > 0. But we can relace both stes with stes (g +, and (h, 0 and make the ath shorter. Therefore all the stes of the form (g, are of the form (,. The remaining stes can be combined to stes of the form (, 0. Lemma 4. Let be fixed. Let < c. The m.r.error of the neighborhood + sequence N c is equal to E D if c = and equal to E B if c assumes the value c B from Corollary. Proof. Because of symmetry it suffices only to consider oints (n, k with 0 k n. First let c =. By definition N(n, k = n + k for (n, k M. Hence the induced length function satisfies W (n, k (n, k for all (n, k Z. Thus min W (n, k n + k
15 Aroximation of the Euclidean Distance by Chamfer Distances 43 where the minimum is taken over all (n, k Z with (n, k (0, 0. On the other hand, by Lemma 3, the shortest N ath from (0, 0 to (m, k with 0 k m consists of stes of the forms (, 0 and (, which have lengths and +, resectively. Hence W (m, k = W (m, k for 0 k m. If n = m + r with 0 r <, then W (n, k W (m, k <. Note that in view the roof of Theorem (in articular, since h 0 (, c h i (, c for all i there we have W (m, k lim su (m,k (m, k 0 k m Thus on the one hand it follows that W (n, k lim su (n,k (n, k 0 k n = lim su (m,k 0 k m W (m, k. (m, k W (m, k = lim su (m,k (m, k 0 k m W (m, k W (m, k lim su = lim su (m,k (m, k (m,k (m, k 0 k m 0 k m = lim su (m,k 0 k m W (m, k. (m, k On the other hand, by W (m, k W (m, k for all m, and k, we also have that W (m, k lim su (m,k (m, k 0 k m W (m, k lim su (m,k (m, k 0 k m W (n, k = lim su (n,k (n, k. 0 k n Hence and by W (n, k lim su (n,k (n, k 0 k n W (m, k = lim su (m,k (m, k 0 k m W (m, k lim su = + E D, (m,k (m, k the m.r.error of N equals E D. Next let c = c B = E B. Then + < c <, and, by construction, W c (, 0 =, W c (, k = c + k for 0 < k, and W c (n, k = c n + k for 0 < k n. Hence min (n,k M W c (n, k n + k = c = E B. Thus W c B lim inf (n, k = E B. (n,k (n, k
16 44 András Hajdu, Lajos Hajdu, and Robert Tijdeman On the other hand, by Lemma 3, the shortest N c ath from (0, 0 to (m, k with 0 k m consists of stes of the form (, 0 and (,. By a similar reasoning as above we obtain that W c B lim su (n, k = + E B. (n,k (n, k Thus the m.r.error of N c equals E B. Theorem 4. For every we have E B = E B and E D = E D. Proof. We first consider the D-case. Suose the neighborhood N on M induces a length function W : Z R 0 such that W ( v v for all v Z and W has m.r.error E D. It can only imrove the m.r.error if we relace the value N(n, k for some (n, k M with a smaller value (n, k. Therefore we may assume without loss of generality that N = N. Hence E D = E D. Now we turn to the B-case. Suose a neighborhood N on M induces a length function W such that W (n, 0 = W (0, n = n for n Z and ( E B v W ( v ( + E B v for all v Z. Without loss of generality we may relace all values N(n, k for (n, k M with n if k = 0, with k if n = 0, and with ( E B (n, k otherwise. Thus E B equals the m.r.error of the neighborhood sequence N E B. We know from Lemma 4 and Corollary that if c = c B, then the m.r.error of N c equals E B = c B. Hence E B E B. From N(n, k ( E B (n, k c B (n, k for all (n, k M we obtain W ( v W c B ( v for all v Z. Hence + E B W ( v W c B = inf lim su lim su ( v = + E B N v v v v by Lemma 4. Thus E B = E B. Finally, we comute the minimal m.r.error E C for the class of arbitrary neighborhoods N defined on M. Observe that the m.r.error E C is attained by the length function W corresonding to the neighborhood N defined by w( v = ( E C v for v M, since N( v v should not assume a smaller value than E C and the limsu-value cannot increase if we decrease some w( v. Clearly, the length function W corresonding to N is just W E C where W is the length function on N. Recall that N has m.r.error E D. Therefore we have + E C = lim su v W ( v v = ( + E D ( E C. (3 By a simle calculation we get E C = ED. So we have roved +E D
17 Aroximation of the Euclidean Distance by Chamfer Distances 4 Theorem. For every we have E C = ED + E D = = ( O 4. Remark. Observe that E B is about 37% larger than E C. This is the rice to be aid for the restriction W (n, 0 = n for n Z. The value of E D is about twice the error E C. This is due to the fact that the negative and ositive deviations in E C are added to the ositive deviation in E D. Conclusion In this aer, we have determined the smallest ossible maximum relative error of chamfer distances with resect to the Euclidean distance under various conditions. We have dealt with aroximating distances from three main asects: suosing that a horizontal/vertical ste has a weight in the local chamfer neighborhoods, majorating the Euclidean distance, and also without any constraint. We have calculated otimal weights for small ( and 7 7 neighborhoods in a certain case, as well. Our framework is embedded in the theory of neighborhood sequences with ossible generalizations in this field. References [] Bailey, D.G. An efficient euclidean distance transform. Lecture Notes in Comuter Science, 33: , 004. [] Borgefors, G. Distance transformations in arbitrary dimensions. Comuter Vision, Grahics, and Image Processing, 7:3 34, 984. [3] Borgefors, G. Distance transformations in digital images. Comuter Vision, Grahics, and Image Processing, 34:344 37, 986. [4] Borgefors, G. Hierarchical chamfer matching: a arametric edge matching algorithm. IEEE Transactions on Pattern Analysis and Machine Intelligence, 0:849 86, 988. [] Butt, M.A. and Maragos, P. Otimum design of chamfer distance transforms. IEEE Transactions on Image Processing, 7: , 998. [6] Coquin, D. and Bolon, Ph. Discrete distance oerator on rectangular grids. Pattern Recognition Letters, 6:9 93, 99. [7] Cuisenaire, O. Distance Transformation, Fast Algorithms and Alications to Medical Image Processing. PhD thesis, Université Catholique de Louvain, 999.
18 46 András Hajdu, Lajos Hajdu, and Robert Tijdeman [8] Danielsson, P.E. Euclidean distance maing. Comuter Grahics and Image Processing, 4:7 48, 980. [9] Das, P.P., Chakrabarti, P.P., and Chatterji, B.N. Distance functions in digital geometry. Information Sciences, 4:3 36, 987. [0] Fazekas, A., Hajdu, A., and Hajdu, L. Lattice of generalized neighbourhood sequences in nd and d. Publ. Math. Debrecen, 60:40 47, 00. [] Fouard, C., Strand, R., and Borgefors, G. Weighted distance transforms generalized to modules and their comutation on oint lattices. Pattern Recognition, 40:43 474, 007. [] Hajdu, A., Hajdu, L., and Tijdeman, R. General neighborhood sequences in Z n. Discrete Al. Math., :07, 007. [3] Jr, C.R. Maurer, Qi, R., and Raghavan, V. A linear time algorithm for comuting exact euclidean distance transforms of binary images in arbitrary dimensions. IEEE Transactions on Pattern Analysis and Machine Intelligence, :6 70, 003. [4] Nagy, B. Distance with generalized neighbourhood sequences in nd and d. Discrete Al. Math., 6:344 3, 008. [] Nagy, B. and Strand, R. Neighborhood sequences in the diamond grid: algorithms with two and three neighbors. International Journal of Imaging Systems and Technology, 9:46 7, 009. [6] Nagy, B. and Strand, R. Neighborhood sequences on nd hexagonal/facecentered-cubic grids. Lecture Notes in Comuter Science, 8:96 08, 009. [7] Scholtus, S. Chamfer Distances with Integer Neighborhoods, Master Thesis, Leiden University, The Netherlands [8] Scholtus, S. and Tijdeman, R. Chamfer distances with integer neighborhoods. Technical reort, Leiden University, The Netherlands, 006. [9] Shih, F.Y. and Wu, Y.T. Fast euclidean distance transformation in two scans using a 3 3 neighborhood. Comuter Vision and Image Understanding, 93:9 0, 004. [0] Strand, R., Nagy, B., Fouard, C., and Borgefors, G. Generating distance mas with neighbourhood sequences. Lecture Notes in Comuter Science, 44:9 307, 006. [] Thiel, E. Les Distances de Chamfrein en Analyse d Images: Fondements et Alications. PhD thesis, Université Joseh Fourier de Grenoble, 994. [] Verwer, B.H.J. Distance Transforms: Matrics, Algorithms and Alications. PhD thesis, Techn. Univ. Delft, 99.
19 Aroximation of the Euclidean Distance by Chamfer Distances 47 [3] Verwer, B.H.J. Local distances for distance transforms in two and three dimensions. Pattern Recognition Letters, :67 68, 99. [4] Yamashita, M. and Ibaraki, T. Distances defined by neighbourhood sequences. Pattern Recognition, 9:37 46, 986. Received 6th February 00
Approximation of the Euclidean distance by chamfer distances
Approximation of the Euclidean distance by chamfer distances András Hajdu a, Lajos Hajdu b Robert Tijdeman c a Faculty of Informatics, University of Debrecen, H-4010 Debrecen, P.O.Box 1. b Number Theory
More informationImprovement on the Decay of Crossing Numbers
Grahs and Combinatorics 2013) 29:365 371 DOI 10.1007/s00373-012-1137-3 ORIGINAL PAPER Imrovement on the Decay of Crossing Numbers Jakub Černý Jan Kynčl Géza Tóth Received: 24 Aril 2007 / Revised: 1 November
More informationSums of independent random variables
3 Sums of indeendent random variables This lecture collects a number of estimates for sums of indeendent random variables with values in a Banach sace E. We concentrate on sums of the form N γ nx n, where
More informationElementary Analysis in Q p
Elementary Analysis in Q Hannah Hutter, May Szedlák, Phili Wirth November 17, 2011 This reort follows very closely the book of Svetlana Katok 1. 1 Sequences and Series In this section we will see some
More informationResearch Article An iterative Algorithm for Hemicontractive Mappings in Banach Spaces
Abstract and Alied Analysis Volume 2012, Article ID 264103, 11 ages doi:10.1155/2012/264103 Research Article An iterative Algorithm for Hemicontractive Maings in Banach Saces Youli Yu, 1 Zhitao Wu, 2 and
More informationAnalysis of some entrance probabilities for killed birth-death processes
Analysis of some entrance robabilities for killed birth-death rocesses Master s Thesis O.J.G. van der Velde Suervisor: Dr. F.M. Sieksma July 5, 207 Mathematical Institute, Leiden University Contents Introduction
More informationON THE LEAST SIGNIFICANT p ADIC DIGITS OF CERTAIN LUCAS NUMBERS
#A13 INTEGERS 14 (014) ON THE LEAST SIGNIFICANT ADIC DIGITS OF CERTAIN LUCAS NUMBERS Tamás Lengyel Deartment of Mathematics, Occidental College, Los Angeles, California lengyel@oxy.edu Received: 6/13/13,
More informationEstimation of the large covariance matrix with two-step monotone missing data
Estimation of the large covariance matrix with two-ste monotone missing data Masashi Hyodo, Nobumichi Shutoh 2, Takashi Seo, and Tatjana Pavlenko 3 Deartment of Mathematical Information Science, Tokyo
More information4. Score normalization technical details We now discuss the technical details of the score normalization method.
SMT SCORING SYSTEM This document describes the scoring system for the Stanford Math Tournament We begin by giving an overview of the changes to scoring and a non-technical descrition of the scoring rules
More informationApproximating min-max k-clustering
Aroximating min-max k-clustering Asaf Levin July 24, 2007 Abstract We consider the roblems of set artitioning into k clusters with minimum total cost and minimum of the maximum cost of a cluster. The cost
More informationPositive decomposition of transfer functions with multiple poles
Positive decomosition of transfer functions with multile oles Béla Nagy 1, Máté Matolcsi 2, and Márta Szilvási 1 Deartment of Analysis, Technical University of Budaest (BME), H-1111, Budaest, Egry J. u.
More informationMarch 4, :21 WSPC/INSTRUCTION FILE FLSpaper2011
International Journal of Number Theory c World Scientific Publishing Comany SOLVING n(n + d) (n + (k 1)d ) = by 2 WITH P (b) Ck M. Filaseta Deartment of Mathematics, University of South Carolina, Columbia,
More informationMATH 2710: NOTES FOR ANALYSIS
MATH 270: NOTES FOR ANALYSIS The main ideas we will learn from analysis center around the idea of a limit. Limits occurs in several settings. We will start with finite limits of sequences, then cover infinite
More informationHENSEL S LEMMA KEITH CONRAD
HENSEL S LEMMA KEITH CONRAD 1. Introduction In the -adic integers, congruences are aroximations: for a and b in Z, a b mod n is the same as a b 1/ n. Turning information modulo one ower of into similar
More informationOn generalizing happy numbers to fractional base number systems
On generalizing hay numbers to fractional base number systems Enriue Treviño, Mikita Zhylinski October 17, 018 Abstract Let n be a ositive integer and S (n) be the sum of the suares of its digits. It is
More informationGaps in Semigroups. Université Pierre et Marie Curie, Paris 6, Equipe Combinatoire - Case 189, 4 Place Jussieu Paris Cedex 05, France.
Gas in Semigrous J.L. Ramírez Alfonsín Université Pierre et Marie Curie, Paris 6, Equie Combinatoire - Case 189, 4 Place Jussieu Paris 755 Cedex 05, France. Abstract In this aer we investigate the behaviour
More informationAdditive results for the generalized Drazin inverse in a Banach algebra
Additive results for the generalized Drazin inverse in a Banach algebra Dragana S. Cvetković-Ilić Dragan S. Djordjević and Yimin Wei* Abstract In this aer we investigate additive roerties of the generalized
More informationVarious Proofs for the Decrease Monotonicity of the Schatten s Power Norm, Various Families of R n Norms and Some Open Problems
Int. J. Oen Problems Comt. Math., Vol. 3, No. 2, June 2010 ISSN 1998-6262; Coyright c ICSRS Publication, 2010 www.i-csrs.org Various Proofs for the Decrease Monotonicity of the Schatten s Power Norm, Various
More informationFeedback-error control
Chater 4 Feedback-error control 4.1 Introduction This chater exlains the feedback-error (FBE) control scheme originally described by Kawato [, 87, 8]. FBE is a widely used neural network based controller
More informationSolution sheet ξi ξ < ξ i+1 0 otherwise ξ ξ i N i,p 1 (ξ) + where 0 0
Advanced Finite Elements MA5337 - WS7/8 Solution sheet This exercise sheets deals with B-slines and NURBS, which are the basis of isogeometric analysis as they will later relace the olynomial ansatz-functions
More informationLinear diophantine equations for discrete tomography
Journal of X-Ray Science and Technology 10 001 59 66 59 IOS Press Linear diohantine euations for discrete tomograhy Yangbo Ye a,gewang b and Jiehua Zhu a a Deartment of Mathematics, The University of Iowa,
More informationarxiv: v2 [math.na] 6 Apr 2016
Existence and otimality of strong stability reserving linear multiste methods: a duality-based aroach arxiv:504.03930v [math.na] 6 Ar 06 Adrián Németh January 9, 08 Abstract David I. Ketcheson We rove
More informationUniform Law on the Unit Sphere of a Banach Space
Uniform Law on the Unit Shere of a Banach Sace by Bernard Beauzamy Société de Calcul Mathématique SA Faubourg Saint Honoré 75008 Paris France Setember 008 Abstract We investigate the construction of a
More informationMODELING THE RELIABILITY OF C4ISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL
Technical Sciences and Alied Mathematics MODELING THE RELIABILITY OF CISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL Cezar VASILESCU Regional Deartment of Defense Resources Management
More informationOn Wald-Type Optimal Stopping for Brownian Motion
J Al Probab Vol 34, No 1, 1997, (66-73) Prerint Ser No 1, 1994, Math Inst Aarhus On Wald-Tye Otimal Stoing for Brownian Motion S RAVRSN and PSKIR The solution is resented to all otimal stoing roblems of
More informationTranspose of the Weighted Mean Matrix on Weighted Sequence Spaces
Transose of the Weighted Mean Matri on Weighted Sequence Saces Rahmatollah Lashkariour Deartment of Mathematics, Faculty of Sciences, Sistan and Baluchestan University, Zahedan, Iran Lashkari@hamoon.usb.ac.ir,
More informationRepresenting integers as linear combinations of powers
ubl. Math. Debrecen Manuscript (August 15, 2011) Representing integers as linear combinations of powers By Lajos Hajdu and Robert Tijdeman Dedicated to rofessors K. Győry and A. Sárközy on their 70th birthdays
More informationIntroduction to Banach Spaces
CHAPTER 8 Introduction to Banach Saces 1. Uniform and Absolute Convergence As a rearation we begin by reviewing some familiar roerties of Cauchy sequences and uniform limits in the setting of metric saces.
More informationTopic: Lower Bounds on Randomized Algorithms Date: September 22, 2004 Scribe: Srinath Sridhar
15-859(M): Randomized Algorithms Lecturer: Anuam Guta Toic: Lower Bounds on Randomized Algorithms Date: Setember 22, 2004 Scribe: Srinath Sridhar 4.1 Introduction In this lecture, we will first consider
More informationUniformly best wavenumber approximations by spatial central difference operators: An initial investigation
Uniformly best wavenumber aroximations by satial central difference oerators: An initial investigation Vitor Linders and Jan Nordström Abstract A characterisation theorem for best uniform wavenumber aroximations
More information#A64 INTEGERS 18 (2018) APPLYING MODULAR ARITHMETIC TO DIOPHANTINE EQUATIONS
#A64 INTEGERS 18 (2018) APPLYING MODULAR ARITHMETIC TO DIOPHANTINE EQUATIONS Ramy F. Taki ElDin Physics and Engineering Mathematics Deartment, Faculty of Engineering, Ain Shams University, Cairo, Egyt
More informationAsymptotically Optimal Simulation Allocation under Dependent Sampling
Asymtotically Otimal Simulation Allocation under Deendent Samling Xiaoing Xiong The Robert H. Smith School of Business, University of Maryland, College Park, MD 20742-1815, USA, xiaoingx@yahoo.com Sandee
More informationOn the Chvatál-Complexity of Knapsack Problems
R u t c o r Research R e o r t On the Chvatál-Comlexity of Knasack Problems Gergely Kovács a Béla Vizvári b RRR 5-08, October 008 RUTCOR Rutgers Center for Oerations Research Rutgers University 640 Bartholomew
More informationOptimal Design of Truss Structures Using a Neutrosophic Number Optimization Model under an Indeterminate Environment
Neutrosohic Sets and Systems Vol 14 016 93 University of New Mexico Otimal Design of Truss Structures Using a Neutrosohic Number Otimization Model under an Indeterminate Environment Wenzhong Jiang & Jun
More informationIntrinsic Approximation on Cantor-like Sets, a Problem of Mahler
Intrinsic Aroximation on Cantor-like Sets, a Problem of Mahler Ryan Broderick, Lior Fishman, Asaf Reich and Barak Weiss July 200 Abstract In 984, Kurt Mahler osed the following fundamental question: How
More informationSharp gradient estimate and spectral rigidity for p-laplacian
Shar gradient estimate and sectral rigidity for -Lalacian Chiung-Jue Anna Sung and Jiaing Wang To aear in ath. Research Letters. Abstract We derive a shar gradient estimate for ositive eigenfunctions of
More informationThe inverse Goldbach problem
1 The inverse Goldbach roblem by Christian Elsholtz Submission Setember 7, 2000 (this version includes galley corrections). Aeared in Mathematika 2001. Abstract We imrove the uer and lower bounds of the
More informationConvex Optimization methods for Computing Channel Capacity
Convex Otimization methods for Comuting Channel Caacity Abhishek Sinha Laboratory for Information and Decision Systems (LIDS), MIT sinhaa@mit.edu May 15, 2014 We consider a classical comutational roblem
More informationSYMPLECTIC STRUCTURES: AT THE INTERFACE OF ANALYSIS, GEOMETRY, AND TOPOLOGY
SYMPLECTIC STRUCTURES: AT THE INTERFACE OF ANALYSIS, GEOMETRY, AND TOPOLOGY FEDERICA PASQUOTTO 1. Descrition of the roosed research 1.1. Introduction. Symlectic structures made their first aearance in
More informationDiscrete Applied Mathematics. Weighted distances based on neighborhood sequences for point-lattices
Discrete Applied Mathematics 157 009 641 65 Contents lists available at ScienceDirect Discrete Applied Mathematics journal homepage: wwwelseviercom/locate/dam Weighted distances based on neighborhood sequences
More informationAn Estimate For Heilbronn s Exponential Sum
An Estimate For Heilbronn s Exonential Sum D.R. Heath-Brown Magdalen College, Oxford For Heini Halberstam, on his retirement Let be a rime, and set e(x) = ex(2πix). Heilbronn s exonential sum is defined
More information#A47 INTEGERS 15 (2015) QUADRATIC DIOPHANTINE EQUATIONS WITH INFINITELY MANY SOLUTIONS IN POSITIVE INTEGERS
#A47 INTEGERS 15 (015) QUADRATIC DIOPHANTINE EQUATIONS WITH INFINITELY MANY SOLUTIONS IN POSITIVE INTEGERS Mihai Ciu Simion Stoilow Institute of Mathematics of the Romanian Academy, Research Unit No. 5,
More information2 Asymptotic density and Dirichlet density
8.785: Analytic Number Theory, MIT, sring 2007 (K.S. Kedlaya) Primes in arithmetic rogressions In this unit, we first rove Dirichlet s theorem on rimes in arithmetic rogressions. We then rove the rime
More information2 Asymptotic density and Dirichlet density
8.785: Analytic Number Theory, MIT, sring 2007 (K.S. Kedlaya) Primes in arithmetic rogressions In this unit, we first rove Dirichlet s theorem on rimes in arithmetic rogressions. We then rove the rime
More informationMath 4400/6400 Homework #8 solutions. 1. Let P be an odd integer (not necessarily prime). Show that modulo 2,
MATH 4400 roblems. Math 4400/6400 Homework # solutions 1. Let P be an odd integer not necessarily rime. Show that modulo, { P 1 0 if P 1, 7 mod, 1 if P 3, mod. Proof. Suose that P 1 mod. Then we can write
More informationOn Isoperimetric Functions of Probability Measures Having Log-Concave Densities with Respect to the Standard Normal Law
On Isoerimetric Functions of Probability Measures Having Log-Concave Densities with Resect to the Standard Normal Law Sergey G. Bobkov Abstract Isoerimetric inequalities are discussed for one-dimensional
More information#A37 INTEGERS 15 (2015) NOTE ON A RESULT OF CHUNG ON WEIL TYPE SUMS
#A37 INTEGERS 15 (2015) NOTE ON A RESULT OF CHUNG ON WEIL TYPE SUMS Norbert Hegyvári ELTE TTK, Eötvös University, Institute of Mathematics, Budaest, Hungary hegyvari@elte.hu François Hennecart Université
More informationp-adic Measures and Bernoulli Numbers
-Adic Measures and Bernoulli Numbers Adam Bowers Introduction The constants B k in the Taylor series exansion t e t = t k B k k! k=0 are known as the Bernoulli numbers. The first few are,, 6, 0, 30, 0,
More informationON THE NORM OF AN IDEMPOTENT SCHUR MULTIPLIER ON THE SCHATTEN CLASS
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 00, Number 0, Pages 000 000 S 000-9939XX)0000-0 ON THE NORM OF AN IDEMPOTENT SCHUR MULTIPLIER ON THE SCHATTEN CLASS WILLIAM D. BANKS AND ASMA HARCHARRAS
More informationProperties of a Natural Ordering Relation for Octagonal Neighborhood Sequences
Properties of a Natural Ordering Relation for Octagonal Neighborhood Sequences Attila Fazekas Image Processing Group of Debrecen Faculty of Informatics, University of Debrecen P.O.Box 12, H-4010 Debrecen,
More informationOn Erdős and Sárközy s sequences with Property P
Monatsh Math 017 18:565 575 DOI 10.1007/s00605-016-0995-9 On Erdős and Sárközy s sequences with Proerty P Christian Elsholtz 1 Stefan Planitzer 1 Received: 7 November 015 / Acceted: 7 October 016 / Published
More informationMultiplicative group law on the folium of Descartes
Multilicative grou law on the folium of Descartes Steluţa Pricoie and Constantin Udrişte Abstract. The folium of Descartes is still studied and understood today. Not only did it rovide for the roof of
More informationFor q 0; 1; : : : ; `? 1, we have m 0; 1; : : : ; q? 1. The set fh j(x) : j 0; 1; ; : : : ; `? 1g forms a basis for the tness functions dened on the i
Comuting with Haar Functions Sami Khuri Deartment of Mathematics and Comuter Science San Jose State University One Washington Square San Jose, CA 9519-0103, USA khuri@juiter.sjsu.edu Fax: (40)94-500 Keywords:
More informationSome Unitary Space Time Codes From Sphere Packing Theory With Optimal Diversity Product of Code Size
IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 5, NO., DECEMBER 4 336 Some Unitary Sace Time Codes From Shere Packing Theory With Otimal Diversity Product of Code Size Haiquan Wang, Genyuan Wang, and Xiang-Gen
More informationBEST POSSIBLE DENSITIES OF DICKSON m-tuples, AS A CONSEQUENCE OF ZHANG-MAYNARD-TAO
BEST POSSIBLE DENSITIES OF DICKSON m-tuples, AS A CONSEQUENCE OF ZHANG-MAYNARD-TAO ANDREW GRANVILLE, DANIEL M. KANE, DIMITRIS KOUKOULOPOULOS, AND ROBERT J. LEMKE OLIVER Abstract. We determine for what
More informationRadial Basis Function Networks: Algorithms
Radial Basis Function Networks: Algorithms Introduction to Neural Networks : Lecture 13 John A. Bullinaria, 2004 1. The RBF Maing 2. The RBF Network Architecture 3. Comutational Power of RBF Networks 4.
More informationBy Evan Chen OTIS, Internal Use
Solutions Notes for DNY-NTCONSTRUCT Evan Chen January 17, 018 1 Solution Notes to TSTST 015/5 Let ϕ(n) denote the number of ositive integers less than n that are relatively rime to n. Prove that there
More informationMATH 361: NUMBER THEORY EIGHTH LECTURE
MATH 361: NUMBER THEORY EIGHTH LECTURE 1. Quadratic Recirocity: Introduction Quadratic recirocity is the first result of modern number theory. Lagrange conjectured it in the late 1700 s, but it was first
More informationApplications to stochastic PDE
15 Alications to stochastic PE In this final lecture we resent some alications of the theory develoed in this course to stochastic artial differential equations. We concentrate on two secific examles:
More informationNotes on duality in second order and -order cone otimization E. D. Andersen Λ, C. Roos y, and T. Terlaky z Aril 6, 000 Abstract Recently, the so-calle
McMaster University Advanced Otimization Laboratory Title: Notes on duality in second order and -order cone otimization Author: Erling D. Andersen, Cornelis Roos and Tamás Terlaky AdvOl-Reort No. 000/8
More informationLocation of solutions for quasi-linear elliptic equations with general gradient dependence
Electronic Journal of Qualitative Theory of Differential Equations 217, No. 87, 1 1; htts://doi.org/1.14232/ejqtde.217.1.87 www.math.u-szeged.hu/ejqtde/ Location of solutions for quasi-linear ellitic equations
More informationA sharp generalization on cone b-metric space over Banach algebra
Available online at www.isr-ublications.com/jnsa J. Nonlinear Sci. Al., 10 2017), 429 435 Research Article Journal Homeage: www.tjnsa.com - www.isr-ublications.com/jnsa A shar generalization on cone b-metric
More informationOnline Appendix to Accompany AComparisonof Traditional and Open-Access Appointment Scheduling Policies
Online Aendix to Accomany AComarisonof Traditional and Oen-Access Aointment Scheduling Policies Lawrence W. Robinson Johnson Graduate School of Management Cornell University Ithaca, NY 14853-6201 lwr2@cornell.edu
More informationAn Inverse Problem for Two Spectra of Complex Finite Jacobi Matrices
Coyright 202 Tech Science Press CMES, vol.86, no.4,.30-39, 202 An Inverse Problem for Two Sectra of Comlex Finite Jacobi Matrices Gusein Sh. Guseinov Abstract: This aer deals with the inverse sectral roblem
More informationCHAPTER-II Control Charts for Fraction Nonconforming using m-of-m Runs Rules
CHAPTER-II Control Charts for Fraction Nonconforming using m-of-m Runs Rules. Introduction: The is widely used in industry to monitor the number of fraction nonconforming units. A nonconforming unit is
More informationPETER J. GRABNER AND ARNOLD KNOPFMACHER
ARITHMETIC AND METRIC PROPERTIES OF -ADIC ENGEL SERIES EXPANSIONS PETER J. GRABNER AND ARNOLD KNOPFMACHER Abstract. We derive a characterization of rational numbers in terms of their unique -adic Engel
More informationElementary theory of L p spaces
CHAPTER 3 Elementary theory of L saces 3.1 Convexity. Jensen, Hölder, Minkowski inequality. We begin with two definitions. A set A R d is said to be convex if, for any x 0, x 1 2 A x = x 0 + (x 1 x 0 )
More informationCorrespondence Between Fractal-Wavelet. Transforms and Iterated Function Systems. With Grey Level Maps. F. Mendivil and E.R.
1 Corresondence Between Fractal-Wavelet Transforms and Iterated Function Systems With Grey Level Mas F. Mendivil and E.R. Vrscay Deartment of Alied Mathematics Faculty of Mathematics University of Waterloo
More informationGENERICITY OF INFINITE-ORDER ELEMENTS IN HYPERBOLIC GROUPS
GENERICITY OF INFINITE-ORDER ELEMENTS IN HYPERBOLIC GROUPS PALLAVI DANI 1. Introduction Let Γ be a finitely generated grou and let S be a finite set of generators for Γ. This determines a word metric on
More informationOn Doob s Maximal Inequality for Brownian Motion
Stochastic Process. Al. Vol. 69, No., 997, (-5) Research Reort No. 337, 995, Det. Theoret. Statist. Aarhus On Doob s Maximal Inequality for Brownian Motion S. E. GRAVERSEN and G. PESKIR If B = (B t ) t
More informationIDENTIFYING CONGRUENCE SUBGROUPS OF THE MODULAR GROUP
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 24, Number 5, May 996 IDENTIFYING CONGRUENCE SUBGROUPS OF THE MODULAR GROUP TIM HSU (Communicated by Ronald M. Solomon) Abstract. We exhibit a simle
More informationGOOD MODELS FOR CUBIC SURFACES. 1. Introduction
GOOD MODELS FOR CUBIC SURFACES ANDREAS-STEPHAN ELSENHANS Abstract. This article describes an algorithm for finding a model of a hyersurface with small coefficients. It is shown that the aroach works in
More informationA Note on Guaranteed Sparse Recovery via l 1 -Minimization
A Note on Guaranteed Sarse Recovery via l -Minimization Simon Foucart, Université Pierre et Marie Curie Abstract It is roved that every s-sarse vector x C N can be recovered from the measurement vector
More informationThe non-stochastic multi-armed bandit problem
Submitted for journal ublication. The non-stochastic multi-armed bandit roblem Peter Auer Institute for Theoretical Comuter Science Graz University of Technology A-8010 Graz (Austria) auer@igi.tu-graz.ac.at
More informationMultiplicity of weak solutions for a class of nonuniformly elliptic equations of p-laplacian type
Nonlinear Analysis 7 29 536 546 www.elsevier.com/locate/na Multilicity of weak solutions for a class of nonuniformly ellitic equations of -Lalacian tye Hoang Quoc Toan, Quô c-anh Ngô Deartment of Mathematics,
More informationCombining Logistic Regression with Kriging for Mapping the Risk of Occurrence of Unexploded Ordnance (UXO)
Combining Logistic Regression with Kriging for Maing the Risk of Occurrence of Unexloded Ordnance (UXO) H. Saito (), P. Goovaerts (), S. A. McKenna (2) Environmental and Water Resources Engineering, Deartment
More informationRound-off Errors and Computer Arithmetic - (1.2)
Round-off Errors and Comuter Arithmetic - (.). Round-off Errors: Round-off errors is roduced when a calculator or comuter is used to erform real number calculations. That is because the arithmetic erformed
More information16.2. Infinite Series. Introduction. Prerequisites. Learning Outcomes
Infinite Series 6.2 Introduction We extend the concet of a finite series, met in Section 6., to the situation in which the number of terms increase without bound. We define what is meant by an infinite
More informationA Social Welfare Optimal Sequential Allocation Procedure
A Social Welfare Otimal Sequential Allocation Procedure Thomas Kalinowsi Universität Rostoc, Germany Nina Narodytsa and Toby Walsh NICTA and UNSW, Australia May 2, 201 Abstract We consider a simle sequential
More informationOn the irreducibility of a polynomial associated with the Strong Factorial Conjecture
On the irreducibility of a olynomial associated with the Strong Factorial Conecture Michael Filaseta Mathematics Deartment University of South Carolina Columbia, SC 29208 USA E-mail: filaseta@math.sc.edu
More informationTRACES OF SCHUR AND KRONECKER PRODUCTS FOR BLOCK MATRICES
Khayyam J. Math. DOI:10.22034/kjm.2019.84207 TRACES OF SCHUR AND KRONECKER PRODUCTS FOR BLOCK MATRICES ISMAEL GARCÍA-BAYONA Communicated by A.M. Peralta Abstract. In this aer, we define two new Schur and
More informationTHE 2D CASE OF THE BOURGAIN-DEMETER-GUTH ARGUMENT
THE 2D CASE OF THE BOURGAIN-DEMETER-GUTH ARGUMENT ZANE LI Let e(z) := e 2πiz and for g : [0, ] C and J [0, ], define the extension oerator E J g(x) := g(t)e(tx + t 2 x 2 ) dt. J For a ositive weight ν
More information2 K. ENTACHER 2 Generalized Haar function systems In the following we x an arbitrary integer base b 2. For the notations and denitions of generalized
BIT 38 :2 (998), 283{292. QUASI-MONTE CARLO METHODS FOR NUMERICAL INTEGRATION OF MULTIVARIATE HAAR SERIES II KARL ENTACHER y Deartment of Mathematics, University of Salzburg, Hellbrunnerstr. 34 A-52 Salzburg,
More informationCONVOLVED SUBSAMPLING ESTIMATION WITH APPLICATIONS TO BLOCK BOOTSTRAP
Submitted to the Annals of Statistics arxiv: arxiv:1706.07237 CONVOLVED SUBSAMPLING ESTIMATION WITH APPLICATIONS TO BLOCK BOOTSTRAP By Johannes Tewes, Dimitris N. Politis and Daniel J. Nordman Ruhr-Universität
More informationStrong Matching of Points with Geometric Shapes
Strong Matching of Points with Geometric Shaes Ahmad Biniaz Anil Maheshwari Michiel Smid School of Comuter Science, Carleton University, Ottawa, Canada December 9, 05 In memory of Ferran Hurtado. Abstract
More information16.2. Infinite Series. Introduction. Prerequisites. Learning Outcomes
Infinite Series 6. Introduction We extend the concet of a finite series, met in section, to the situation in which the number of terms increase without bound. We define what is meant by an infinite series
More informationANALYTIC NUMBER THEORY AND DIRICHLET S THEOREM
ANALYTIC NUMBER THEORY AND DIRICHLET S THEOREM JOHN BINDER Abstract. In this aer, we rove Dirichlet s theorem that, given any air h, k with h, k) =, there are infinitely many rime numbers congruent to
More informationA generalization of Amdahl's law and relative conditions of parallelism
A generalization of Amdahl's law and relative conditions of arallelism Author: Gianluca Argentini, New Technologies and Models, Riello Grou, Legnago (VR), Italy. E-mail: gianluca.argentini@riellogrou.com
More informationCMSC 425: Lecture 4 Geometry and Geometric Programming
CMSC 425: Lecture 4 Geometry and Geometric Programming Geometry for Game Programming and Grahics: For the next few lectures, we will discuss some of the basic elements of geometry. There are many areas
More informationTHE ERDÖS - MORDELL THEOREM IN THE EXTERIOR DOMAIN
INTERNATIONAL JOURNAL OF GEOMETRY Vol. 5 (2016), No. 1, 31-38 THE ERDÖS - MORDELL THEOREM IN THE EXTERIOR DOMAIN PETER WALKER Abstract. We show that in the Erd½os-Mordell theorem, the art of the region
More informationNEW SUBCLASS OF MULTIVALENT HYPERGEOMETRIC MEROMORPHIC FUNCTIONS
Kragujevac Journal of Mathematics Volume 42(1) (2018), Pages 83 95. NEW SUBCLASS OF MULTIVALENT HYPERGEOMETRIC MEROMORPHIC FUNCTIONS M. ALBEHBAH 1 AND M. DARUS 2 Abstract. In this aer, we introduce a new
More informationState Estimation with ARMarkov Models
Deartment of Mechanical and Aerosace Engineering Technical Reort No. 3046, October 1998. Princeton University, Princeton, NJ. State Estimation with ARMarkov Models Ryoung K. Lim 1 Columbia University,
More informationSmall Zeros of Quadratic Forms Mod P m
International Mathematical Forum, Vol. 8, 2013, no. 8, 357-367 Small Zeros of Quadratic Forms Mod P m Ali H. Hakami Deartment of Mathematics, Faculty of Science, Jazan University P.O. Box 277, Jazan, Postal
More informationSets of Real Numbers
Chater 4 Sets of Real Numbers 4. The Integers Z and their Proerties In our revious discussions about sets and functions the set of integers Z served as a key examle. Its ubiquitousness comes from the fact
More informationIntroduction Consider a set of jobs that are created in an on-line fashion and should be assigned to disks. Each job has a weight which is the frequen
Ancient and new algorithms for load balancing in the L norm Adi Avidor Yossi Azar y Jir Sgall z July 7, 997 Abstract We consider the on-line load balancing roblem where there are m identical machines (servers)
More informationNONLINEAR OPTIMIZATION WITH CONVEX CONSTRAINTS. The Goldstein-Levitin-Polyak algorithm
- (23) NLP - NONLINEAR OPTIMIZATION WITH CONVEX CONSTRAINTS The Goldstein-Levitin-Polya algorithm We consider an algorithm for solving the otimization roblem under convex constraints. Although the convexity
More informationOXFORD UNIVERSITY. MATHEMATICS, JOINT SCHOOLS AND COMPUTER SCIENCE WEDNESDAY 4 NOVEMBER 2009 Time allowed: hours
OXFORD UNIVERSITY MATHEMATICS, JOINT SCHOOLS AND COMPUTER SCIENCE WEDNESDAY 4 NOVEMBER 2009 Time allowed: 2 1 2 hours For candidates alying for Mathematics, Mathematics & Statistics, Comuter Science, Mathematics
More informationMulti-Operation Multi-Machine Scheduling
Multi-Oeration Multi-Machine Scheduling Weizhen Mao he College of William and Mary, Williamsburg VA 3185, USA Abstract. In the multi-oeration scheduling that arises in industrial engineering, each job
More informationEconometrica Supplementary Material
Econometrica Sulementary Material SUPPLEMENT TO WEAKLY BELIEF-FREE EQUILIBRIA IN REPEATED GAMES WITH PRIVATE MONITORING (Econometrica, Vol. 79, No. 3, May 2011, 877 892) BY KANDORI,MICHIHIRO IN THIS SUPPLEMENT,
More information3 Properties of Dedekind domains
18.785 Number theory I Fall 2016 Lecture #3 09/15/2016 3 Proerties of Dedekind domains In the revious lecture we defined a Dedekind domain as a noetherian domain A that satisfies either of the following
More information