Stablty results for unquely determned sets from two drectons n dscrete tomography Brgt van Dalen arxv:0800316v [mathco] 7 Nov 008 Mathematsch Insttuut, Unverstet Leden, Nels Bohrweg 1, 333 CA Leden, The Netherlands Abstract In ths paper we prove several new stablty results for the reconstructon of bnary mages from two projectons We consder an orgnal mage that s unquely determned by ts projectons and possble reconstructons from slghtly dfferent projectons We show that for a gven dfference n the projectons, the reconstructon can only be dsjont from the orgnal mage f the sze of the mage s not too large We also prove an upper bound for the sze of the mage gven the error n the projectons and the sze of the ntersecton between the mage and the reconstructon Key words: Dscrete tomography, Stablty, Image reconstructon, Bnary mage 1 Introducton Dscrete tomography s concerned wth problems such as reconstructng bnary mages on a lattce from gven projectons n lattce drectons [6] Each pont of a bnary mage has a value equal to zero or one The lne sum of a lne through the mage s the sum of the values of the ponts on ths lne The projecton of the mage n a certan lattce drecton conssts of all the lne sums of the lnes through the mage n ths drecton Several problems related to the reconstructon of bnary mages from two or more projectons have been descrbed n the lterature [6,7] Already n 1957, Ryser gave an algorthm to reconstruct bnary mages from ther horzontal Emal address: dalen@mathledenunvnl (Brgt van Dalen) Preprnt submtted to Elsever 6 January 01
and vertcal projectons and charactersed the set of projectons that correspond to a unque bnary mage [11] For any set of drectons, t s possble to construct mages that are not unquely determned by ther projectons n those drectons [6, Theorem 31] The problem of decdng whether an mage s unquely determned by ts projectons and the problem of reconstructng t are NP-hard for any set of more than two drectons [] Asde from varous nterestng theoretcal problems, dscrete tomography also has applcatons n a wde range of felds The most mportant are electron mcroscopy [8] and medcal magng [5,13], but there are also applcatons n nuclear scence [9,10] and varous other felds [1,15] An nterestng problem n dscrete tomography s the stablty of reconstructons Even f an mage s unquely determned by ts projectons, a very small error n the projectons may lead to a completely dfferent reconstructon [1,3] Alpers et al [1,] showed that n the case of two drectons a total error of at most n the projectons can only cause a small dfference n the reconstructon They also proved a lower bound on the error f the reconstructon s dsjont from the orgnal mage In ths paper we mprove ths bound, and we resolve the open problem of stablty wth a projecton error greater than Notaton and statement of the problems Let F 1 and F be two fnte subsets of Z wth characterstc functons χ 1 and χ (That s, χ h (x, y) = 1 f and only f (x, y) F h, h {1, }) For Z, we defne row as the set {(x, y) Z : x = } We call the ndex of the row For j Z, we defne column j as the set {(x, y) Z : y = j} We call j the ndex of the column Followng matrx notaton, we use row numbers that ncrease when gong downwards and column numbers that ncrease when gong to the rght The row sum r (h) s the number of elements of F h n row, that s r (h) = j Z χ h (, j) The column sum c (h) j of F h s the number of elements of F h n column j, that s c (h) j = Z χ h (, j) We refer to both row and column sums as the lne sums of F h Throughout ths paper, we assume that F 1 s unquely determned by ts row and column sums Such sets were studed by, among others, Ryser [11] and Wang [1] Let a be the number of rows and b the number of columns that
1 3 5 1 3 5 Fg 1 A unquely determned set wth the assumed row and column orderng contan elements of F 1 We renumber the rows and columns such that we have r (1) 1 r (1) r (1) a > 0, c (1) 1 c (1) c (1) b > 0, and such that all elements of F are contaned n rows and columns wth postve ndces By [1, Theorem 3] we have the followng property of F 1 (see Fgure 1): n row the elements of F 1 are precsely the ponts (, 1), (, ),, (, r (1) ), n column j the elements of F 1 are precsely the ponts (1, j), (, j),, (c (1) j, j) We wll refer to ths property as the trangular shape of F 1 Everywhere except n Secton 6 we assume that F 1 = F Note that we do not assume F to be unquely determned As F 1 and F are dfferent and F 1 s unquely determned by ts lne sums, F cannot have exactly the same lne sums as F 1 Defne the dfference or error n the lne sums as c (1) j c () j + r (1) r () j 1 1 As n general t s t + s mod, the above expresson s congruent to j 1 ( (1) c j + c () j ) + 1 ( (1) r + r () ) F1 + F 0 mod, hence the error n the lne sums s always even We wll denote t by α, where α s a postve nteger For notatonal convenence, we wll often wrte p for F 1 F We consder two problems concernng stablty Problem 1 Suppose F 1 F = How large can F 1 be n terms of α? 3
Alpers et al [, Theorem 9] proved that F 1 α They also showed that there s no constant c such that F 1 cα for all F 1 and F In Secton of ths paper we wll prove the new bound F 1 α(1 + log α) and show that ths bound s asymptotcally sharp Problem How small can F 1 F be n terms of F 1 and α, or, equvalently, how large can F 1 be n terms of F 1 F and α? Alpers ([1, Theorem 5118]) showed n the case α = 1 that F 1 F F 1 + 1 F 1 + 1 Ths bound s sharp: f F 1 = 1 n(n + 1) for some postve nteger n, then there exsts an example for whch equalty holds A smlar result s stated n [, Theorem 19] Whle [1,] only deal wth the case α = 1, we wll gve stablty results for general α In Secton 5 we wll gve two dfferent upper bounds for F 1 The bounds have dfferent asymptotc behavour Wrtng p for F 1 F, the second bound reduces to F 1 p + 1 + p + 1 n case α = 1, whch s equvalent to p F 1 F 1 Hence the second new bound can be vewed as a generalsaton of Alpers bound The frst new bound s dfferent and better n the case that α s very large In Secton 6 we wll generalse the results to the case F 1 F 3 Starcases Alpers ntroduced the noton of a starcase to characterse F 1 F n the case α = 1 We wll use a slghtly dfferent defnton and then show that for general α the symmetrc dfference F 1 F conssts of α starcases Defnton 3 A set of ponts (p 1, p,, p n ) n Z s called a starcase f the followng two condtons are satsfed: for each wth 1 n 1 one of the ponts p and p +1 s an element of F 1 \F and the other s an element of F \F 1 ;
Fg A starcase The set F 1 conssts of the whte and the black-and-whte ponts, whle F conssts of the black and the black-and-whte ponts The starcase s ndcated by the dashed lne segments ether for all the ponts p and p +1 are n the same column and the ponts p +1 and p + are n the same row, or for all the ponts p and p +1 are n the same row and the ponts p +1 and p + are n the same column Ths defnton s dfferent from [1,] n the followng way Frstly, the number of ponts does not need to be even Secondly, the ponts p 1 and p n can both be ether n F 1 \F or n F \F 1 So ths defnton s slghtly more general than the one used n [1,] for the case α = 1 Consder a pont p F 1 \F of a starcase (p 1, p,, p n ) Assume p 1 s n the same column as p and p +1 s n the same row as p Because of the trangular shape of F 1, the row ndex of p 1 must be larger than the row ndex of p, and the column ndex of p +1 must be larger than the column ndex of p Therefore, the starcase looks lke a real-world starcase (see Fgure ) From now on, we assume for all starcases that p 1 s the pont wth the largest row ndex and the smallest column ndex, whle p n s the pont wth the smallest row ndex and the largest column ndex We say that the starcase begns wth p 1 and ends wth p n Lemma Let F 1 and F be fnte subsets of Z such that F 1 s unquely determned by ts row and column sums, and F 1 = F Let α be defned as n Secton Then the set F 1 F s the dsjont unon of α starcases PROOF We wll construct the starcases one by one and delete them from F 1 F For a subset A of F 1 F, defne ρ (A) = {j Z : (, j) A F 1 } {j Z : (, j) A F }, Z, σ j (A) = { Z : (, j) A F 1 } { Z : (, j) A F }, j Z, τ(a) = ρ (A) + j σ j (A) We have α = τ(f 1 F ) 5
Assume that the rows and columns are ordered as n Secton Because of the trangular shape of F 1, for any pont (, j) F 1 \F and any pont (k, l) F \F 1 we then have k > or l > j Suppose we have deleted some starcases and are now left wth a non-empty subset A of F 1 F Let (p 1, p,, p n ) be a starcase of maxmal length that s contaned n A Let (x 1, y 1 ) and (x n, y n ) be the coordnates of the ponts p 1 and p n respectvely Each of those two ponts can be ether n A F 1 or n A F, so there are four dfferent cases (If n = 1, so p 1 and p n are the same pont, then there are only two cases) We consder two cases; the other two are smlar Frst suppose p 1 A F 1 and p n A F If (x, y 1 ) s a pont of A F n the same column as p 1, then x > x 1, so we can extend the starcase by addng ths pont That contradcts the maxmal length of the starcase So there are no ponts of A F n column y 1 Therefore σ y1 (A) > 0 Smlarly, snce p n A F, there are no ponts of A F 1 n the same column as p n Therefore σ yn (A) < 0 All rows and all columns that contan ponts of the starcase, except columns y 1 and y n, contan exactly two ponts of the starcase, one n A F 1 and one n A F Let A = A\{p 1, p,, p n } Then ρ (A ) = ρ (A) for all, and σ j (A ) = σ j (A) for all j y 1, y n Furthermore, σ y1 (A ) = σ y1 (A) 1 and σ yn (A ) = σ yn (A) + 1 Snce σ y1 (A) > 0 and σ yn (A) < 0, ths gves τ(a ) = τ(a) Now consder the case p 1 A F 1 and p n A F 1 As above, we have σ y1 (A) > 0 Suppose (x n, y) s a pont of A F n the same row as p n Then y > y n, so we can extend the starcase by addng ths pont That contradcts the maxmal length of the starcase So there are no ponts of A F n row x n Therefore ρ xn (A) > 0 All rows and all columns that contan ponts of the starcase, except column y 1 and row x n, contan exactly two ponts of the starcase, one n A F 1 and one n A F Let A = A\{p 1, p,, p n } Then ρ (A ) = ρ (A) for all x n, and σ j (A ) = σ j (A) for all j y 1 Furthermore, σ y1 (A ) = σ y1 (A) 1 and ρ xn (A ) = ρ xn (A) 1 Snce σ y1 (A) > 0 and ρ xn (A) > 0, ths gves τ(a ) = τ(a) We can contnue deletng starcases n ths way untl all ponts of F 1 F have been deleted Snce τ(a) 0 for all subsets A F 1 F, ths must happen after deletng exactly α starcases Remark 5 Some remarks about the above lemma and ts proof 6
() The α starcases from the prevous lemma have α endponts n total (where we count the same pont twce n case of a starcase consstng of one pont) Each endpont contrbutes a dfference of 1 to the lne sums n one row or column Snce all these dfferences must add up to α, they cannot cancel each other () A starcase consstng of more than one pont can be splt nto two or more starcases So t may be possble to wrte F 1 F as the dsjont unon of more than α starcases However, n that case some of the contrbutons of the endponts of starcases to the dfference n the lne sums cancel each other On the other hand, t s mpossble to decompose F 1 F nto fewer than α starcases () The endponts of a starcase can be n F 1 \F or F \F 1 For a starcase T of whch the two endponts are n dfferent sets, we have T F 1 = T F For a starcase T of whch the two endponts are n the same set, we have T F 1 = 1 + T F or T F = 1 + T F 1 Snce F 1 \F = F \F 1, the number of starcases wth two endponts n F 1 \F must be equal to the number of starcases wth two endponts n F \F 1 Ths mples that of the α endponts, exactly α are n the set F 1 \F and α are n the set F \F 1 Consder a decomposton of F 1 F as n the proof of Lemma We wll now show that for our purposes we may assume that all these starcases begn wth a pont p 1 F 1 \F and end wth a pont p n F \F 1 Suppose there s a starcase begnnng wth a pont (x, y) F \F 1 Then there also exsts a starcase endng wth a pont (x, y ) F 1 \F : otherwse more than half of the α endponts would be n F \F 1, whch s a contradcton to Remark 5() Because of Remark 5() we must have r (1) x < r () x and r (1) x > r() x Let y be such that (x, y ) F 1 F Delete the pont (x, y) from F and add the pont (x, y ) to F Then r x () decreases by 1 and r () x ncreases by 1, so the dfference n the row sums decreases by Meanwhle, the dfference n the column sums ncreases by at most So α does not ncrease, whle F 1, F and F 1 F do not change So the new stuaton s just as good or better than the old one The starcase that began wth (x, y) n the old stuaton now begns wth a pont of F 1 \F The pont that we added becomes the new endpont of the starcase that prevously ended wth (x, y ) Therefore, n our nvestgatons we may assume that all starcases begn wth a pont of F 1 \F and end wth a pont of F \F 1 Ths s an mportant assumpton that we wll use n the proofs throughout the paper An mmedate consequence of t s that r (1) = r () for all The only dfference between correspondng lne sums occurs n the columns 7
A new bound for the dsjont case Usng the concept of starcases, we can prove a new bound for Problem 1 Theorem 6 Let F 1 and F be fnte subsets of Z such that F 1 s unquely determned by ts row and column sums, F 1 = F, and F 1 F = Let α be defned as n Secton Then α α F 1 PROOF Assume that the rows and columns are ordered as n Secton Let a be the number of rows and b the number of columns that contan elements of F 1 Let (k, l) F 1 Then all the ponts n the rectangle {(, j) : 1 k, 1 j l} are elements of F 1 Snce F 1 and F are dsjont, none of the ponts n ths rectangle s an element of F, and all the ponts belong to F 1 F So all of the kl ponts must belong to dfferent starcases, whch mples α kl For all wth 1 a we have (, r (1) ) F 1, hence r (1) an nteger, we have F 1 = a Snce (a, 1) F 1, we have a α, so r (1) α α F 1 a α α Snce r(1) must be Corollary 7 Let F 1, F and α be defned as n Theorem 6 Then F 1 α(1 + log α) PROOF We have α α F 1 α α ( 1 α ) α 1 1 + 1 x dx = α (1 + log α) 8
Fg 3 The constructon from Example 8 wth m = 3 The followng example shows that the upper bound cannot even be mproved by a factor 1 log 07 Example 8 (taken from [1]) Let m 1 be an nteger We construct sets F 1 and F as follows (see also Fgure 3) Row 1: (1, j) F 1 for 1 j m, (1, j) F for m + 1 j m+1 Let 0 l m 1 Row, where l + 1 l+1 : (, j) F 1 for 1 j m l 1, (, j) F for m l 1 + 1 j m l The constructon s almost completely symmetrcal: f (, j) F 1, then (j, ) F 1 ; and f (, j) F wth > 1, then (j, ) F Snce t s clear from the constructon that each row contans exactly as many ponts of F 1 as ponts of F, we conclude that each column j wth j m contans exactly as many ponts of F 1 as ponts of F as well The only dfference n the lne sums occurs n the frst column (whch has m ponts of F 1 and none of F ) and n columns m + 1 up to m+1 (each of whch contans one pont of F and none of F 1 ) So we have Furthermore, F 1 = m + m 1 l=0 α = m Hence for ths famly of examples t holds that l m l 1 = m + m m 1 F 1 = α + 1 α log α, whch s very close to the bound we proved n Corollary 7 9
5 Two bounds for general α In case F 1 and F are not dsjont, we can use an approach very smlar to Secton n order to derve a bound for Problem Theorem 9 Let F 1 and F be fnte subsets of Z such that F 1 s unquely determned by ts row and column sums, and F 1 = F Let α be defned as n Secton, and let p = F 1 F Then F 1 α+p α + p PROOF Assume that the rows and columns are ordered as n Secton Let (k, l) F 1 Then all the ponts n the rectangle {(, j) : 1 k, 1 j l} are elements of F 1 At most p of the ponts n ths rectangle are elements of F, so at least kl p ponts belong to F 1 F None of the ponts n the rectangle s an element of F \F 1, so all of the kl p ponts of F 1 F n the rectangle must belong to dfferent starcases, whch mples α + p kl For all wth 1 a we have (, r (1) ) F 1, hence r (1) an nteger, we have F 1 = a r (1) Snce (a, 1) F 1, we have a α + p, so F 1 α+p a α + p α + p α+p Snce r (1) Corollary 10 Let F 1, F, α and p be defned as n Theorem 9 Then F 1 (α + p)(1 + log(α + p)) must be PROOF Analogous to the proof of Corollary 7 The followng example shows that the upper bound cannot even be mproved 1 p+1 by a factor 07, provded that α > log(p + 1) log log 1 Example 11 10
Fg The constructon from Example 11 wth k = 3 and m = Let k and m be ntegers satsfyng k and m k We construct sets F 1 and F as follows (see also Fgures and 5) Row 1: (1, j) F 1 F for 1 j k 1, (1, j) F 1 for k 1 + 1 j m k 1 + 1, (1, j) F for m k 1 + j m+1 k k 1 + Let 0 l k Row, where l + 1 l+1 : (, 1) F 1 F, (, j) F 1 for j m l 1 k l + 1, (, j) F for m l 1 k l + j m l k l 1 + 1 Let k 1 l m k Row, where l + 1 l+1 : (, j) F 1 for 1 j m l 1, (, j) F for m l 1 + 1 j m l Let m k + 1 l m 1 Row, where l l m+k 1 + l+1 l m+k + 1: (, j) F 1 for 1 j m l 1, (, j) F for m l 1 + 1 j m l The constructon s almost symmetrcal: f (, j) F 1, then (j, ) F 1 ; f (, j) F 1 F, then (j, ) F 1 F ; and f (, j) F wth > 1, then (j, ) F Snce t s clear from the constructon that each row contans exactly as many ponts of F 1 as ponts of F, we conclude that each column j wth j m k 1 + 1 contans exactly as many ponts of F 1 as ponts of F as well The only dfference n the lne sums occurs n the frst column (whch has m k 1 + 1 ponts of F 1 and only k 1 of F ) and n columns m k 1 + up to m+1 k k 1 + (each of whch contans one pont of F and none of F 1 ) So we have 11
Fg 5 The constructon from Example 11 wth k = and m = α = 1 ( ( m k 1 + 1) k 1 + ( m+1 k k 1 + ) ( m k 1 + 1) ) = m k + 1 It s easy to see that p = F 1 F = k 1 Now we count the number of elements of F 1 Row 1 contans m k 1 + 1 elements of F 1 Let 0 l k Rows l + 1 up to l+1 together contan l ( m l 1 k l + 1) = m 1 k + l elements of F 1 Let k 1 l m k Rows l + 1 up to l+1 together contan l m l 1 = m 1 elements of F 1 Let m k + 1 l m 1 Rows l l m+k 1 + up to l+1 l m+k + 1 together contan ( l l m+k 1 )( m l 1 ) = m 1 k elements of F 1 Hence the number of elements of F 1 s F 1 = m k 1 + 1 + (k 1)( m 1 k ) + k l l=0 +(m k + ) m 1 + (k 1)( m 1 k ) = m + m m 1 + k 1 k k 1 For ths famly of examples we now have F 1 = α + p + α + p log (α + p) + p + 1 p + 1 log (p + 1) 1
We wll now prove another bound, whch s better f p = F 1 F s large compared to α Let u be an nteger such that u = F 1 F We wll frst derve an upper bound on u n terms of a, b and α Then we wll derve a lower bound on F 1 n terms of a, b and α By combnng these two, we fnd an upper bound on u n terms of α and p Lemma 1 Let F 1 and F be fnte subsets of Z such that F 1 s unquely determned by ts row and column sums, and F 1 = F Let α, a and b be defned as n Secton Defne u as u = F 1 F Then we have u α (a + b)(a + b + α 1) PROOF Decompose F 1 F nto α starcases as n Lemma, and let T be the set consstng of these starcases Let T T be a starcase and a + 1 a postve nteger Consder the elements of T F n rows, + 1,, a If such elements exst, then let w (T ) be the largest column ndex that occurs among these elements If there are no elements of T F n those rows, then let w (T ) be equal to the smallest column ndex of an element of T F 1 (no longer restrcted to rows,, a) We have w (T ) 1 Defne W = T T w (T ) Let d be the number of elements of F 1 \F n row Let y 1 < < y d be the column ndces of the elements of F 1 \F n row, and let y 1 < < y d be the column ndces of the elements of F \F 1 n row Let T T be the set of starcases wth elements n row The elements n F \F 1 of these starcases are n columns y 1, y,, y d, hence the set {w (T ) : T T } s equal to the set {y 1, y,, y d } The elements n F 1 \F are n columns y 1, y,, y d and are ether the frst element of a starcase or correspond to an element of F \F 1 n the same column but n a row wth ndex at least + 1 In ether case, for a starcase T T we have w +1 (T ) = y j for some j Hence the set {w +1 (T ) : T T } s equal to the set {y 1, y,, y d } We have d d w +1 (T ) = y j (y d j + 1) = d y d 1 T T j=1 j=1 (d 1)d, and Hence d d w (T ) = y j (y d + j) = d y d + 1 T T j=1 j=1 (d + 1)d 13
W = W +1 + (w (T ) w +1 (T )) T T W +1 + 1 (d + 1)d + 1 (d 1)d = W +1 + d Snce W a+1 α, we fnd W 1 α + d 1 + + d a We may assume that f (x, y) s the endpont of a starcase, then (x, y ) s an element of F 1 F for 1 y < y (e there are no gaps between the endponts and other elements of F 1 F on the same row) After all, by movng the endpont of a starcase to another empty poston on the same row, the error n the columns can only become smaller (f the new poston of the endpont happens to be n the same column as the frst pont of another starcase, n whch case the two starcases fuse together to one) but not larger, and u, a and b do not change So on the other hand, as W 1 s the sum of the column ndces of the endponts of the starcases, we have W 1 (b + 1) + (b + ) + + (b + α) = αb + 1 α(α + 1) We conclude a α + d αb + 1 α(α + 1) Note that a d = u By the Cauchy-Schwarz nequalty, we have so From ths t follows that or, equvalently, By symmetry we also have ( a ) ( a ) ( a ) d 1 d = u, a d u a αb + 1 u α(α + 1) α + a, u αab + 1 α(α 1)a u αab + 1 α(α 1)b 1
Hence u αab + 1 α(α 1)(a + b) Usng that ab a+b, we fnd ( ) (a + b) u (α 1)(a + b) α + = α (a + b)(a + b + α 1) Lemma 13 Let F 1 and F be fnte subsets of Z such that F 1 s unquely determned by ts row and column sums, and F 1 = F Let α, a and b be defned as n Secton Then we have F 1 (a + b) (α + 1) PROOF Wthout loss of generalty, we may assume that all rows and columns that contan elements of F 1 also contan at least one pont F 1 F : f a row or column does not contan any ponts of F 1 F, we may delete t By dong so, F 1 F does not change, whle F 1 becomes smaller, so the stuaton becomes better Frst consder the case r (1) +1 < r (1) α for some We wll show that ths s mpossble If a column does not contan an element of F \F 1, then by the assumpton above t contans an element of F 1 \F, whch must then be the frst pont of a starcase Consder all ponts of F \F 1 and all frst ponts of starcases n columns r +1 + 1, r +1 +,, r Snce these are more than α columns, at least two of those ponts must belong to the same starcase On the other hand, f (x, y) F 1 \F s the frst pont of a starcase wth r +1 < y r, then we have x, so the second pont (x, y ) n the starcase, whch s n F \F 1, must satsfy x and therefore y > r So the second pont cannot also be n one of the columns r +1 + 1, r +1 +,, r If two ponts of F \F 1 n columns r +1 + 1, r +1 +,, r belong to the same starcase, then they must be connected by a pont of F 1 \F n the same columns However, by a smlar argument ths forces the next pont to be outsde the mentoned columns, whle we assumed that t was n those columns We conclude that t s mpossble for row sums of two consecutve rows to dffer by more than α By the same argument, column sums of two consecutve columns cannot dffer by more than α Hence we have r (1) +1 r (1) α for all, and c (1) j+1 c (1) j α for all j 15
Fg 6 The number of ponts of F 1 (ndcated by small black dots) s equal to the grey area We now have r (1) b α, r (1) 3 b α, and so on Also, c (1) a α, c (1) 3 a α, and so on Usng ths, we can derve a lower bound on F 1 for fxed a and b Consder Fgure 6 The ponts of F 1 are ndcated by black dots The number of ponts s equal to the grey area n the pcture, whch conssts of all 1 1-squares wth a pont of F 1 n the upper left corner We can estmate ths area from below by drawng a lne wth slope α through the pont (a + 1, 1) and a lne wth slope 1 through the pont (b + 1, 1); the α area closed n by these two lnes and the two axes s less than or equal to the number of ponts of F 1 For α = 1 those lnes do not have a pont of ntersecton Under the assumpton we made at the begnnng of ths proof, we must n ths case have a = b and the number of ponts of F 1 s equal to a(a + 1) a α + 1 = (a + b) (α + 1), so n ths case we are done In order to compute the area for α we swtch to the usual coordnates n R, see Fgure 7 The equaton of the frst lne s y = αx a, and the equaton of the second lne s y = 1 x 1 b We fnd that the pont of ntersecton s α α gven by ( aα b (x, y) = α 1, bα + a ) α 1 The area of the grey part of Fgure 7 s equal to We now have 1 a aα b α 1 + 1 b bα a α 1 = a α + b α ab (α 1) F 1 α(a + b ) ab (α 1) α (a+b) (a+b) (α 1) = (a + b) (α + 1) Theorem 1 Let F 1 and F be fnte subsets of Z such that 16
(0, 0) (b, 0) y = 1 α x 1 α b (0, a) y = αx a Fg 7 Computng the area bounded by the two lnes and the two axes F 1 s unquely determned by ts row and column sums, and F 1 = F Let α be defned as n Secton, and let p = F 1 F Wrte β = α(α + 1) Then ( α F 1 p + β + β(α 1) + (α + 1)p + β + α 1 ) (α 1) α 16 PROOF Wrte s = a + b for convenence of notaton From Lemma 1 we derve ( α u s + α 1 ) We substtute F 1 = u + p n Lemma 13 and use the above bound for u: ( α s + α 1 ) s + p F 1 (α + 1) Solvng for s, we fnd s α(α + 1) + α(α 1) + (α + 1)p + α(α + 1) = β + β(α + 1) + (α + 1)p + β Fnally we substtute ths n Lemma 1: u ( α β + β(α 1) + (α + 1)p + β + α 1 Ths, together wth F 1 = u + p, yelds the clamed result ) (α 1) α 16 Remark 15 By a straghtforward generalsaton of [, Proposton 13 and Lemma 16], we fnd a bound very smlar to the one n Theorem 1: F 1 p + (α + 1)(α 1 ) + (α + 1) p + (α 1) 17
N α N α Fg 8 The constructon from Example 16 wth N = and α = 3 Theorem 1 says that F 1 s asymptotcally bounded by p + α p + α The next example shows that F 1 can be asymptotcally as large as p + αp + α Example 16 Let N be a postve nteger We construct F 1 and F wth total dfference n the lne sums equal to α as follows (see also Fgure 8) Let (, j) F 1 F for 1 N, 1 j N Furthermore, for 1 N: Let (, j), (j, ) F 1 F for N + 1 j N + (N )α Let (, j), (j, ) F 1 for N + (N )α + 1 j N + (N + 1)α Let (, j), (j, ) F for N + (N + 1)α + 1 j N + (N + )α Fnally, for 1 t α, let (, j) F wth = N + t and j = N + α + 1 t The only dfferences n the lne sums occur n the frst column (a dfference of α) and n columns N + Nα + 1 up to N + Nα + α (a dfference of 1 n each column) We have p = N + 1 N(N 1)α = N + N α Nα, and F 1 = N + 1 N(N + 1)α = N + N α + Nα 18
From the frst equalty we derve N = α (α + 1) + p α + 1 + α (α + 1) Hence F 1 = p + Nα = p + α α + 1 + α p α + 1 + α (α + 1) 6 Generalsaton to unequal szes Untl now, we have assumed that F 1 = F However, we can easly generalse all the results to the case F 1 F Suppose F 1 > F Then there must be a row wth r (1) > r () Let j > b be such that (, j) F and defne F 3 = F {(, j)} We have r (3) = r () + 1, so the error n row has decreased by one, whle the error n column j has ncreased by one In ths way, we can keep addng ponts untl F together wth the extra ponts s just as large as F 1, whle the total dfference n the lne sums s stll α Note that p = F 1 F and F 1 have not changed durng ths process, so the results from Theorem 1 and Corollary 10 are stll vald n exactly the same form Suppose on the other hand that F 1 < F Then there must be a row wth r (1) < r () Let j be such that (, j) F \F 1 and defne F 3 = F \{(, j)} The error n row has now decreased by one, whle the error n column j has at most ncreased by one, so the total error n the lne sums has not ncreased We can keep deletng ponts of F untl there are exactly F 1 ponts left, whle the total dfference n the lne sums s at most α By usng F 1 F = ( F 1 p), we can state the results from Theorem 1 and Corollary 10 n a more symmetrc way, not dependng on the sze of F 1 Theorem 17 Let F 1 and F be fnte subsets of Z such that F 1 s unquely determned by ts row and column sums Let α be defned as n Secton, and let p = F 1 F Wrte β = α(α + 1) Then (1) F 1 F α + (α + p) log(α + p) () F 1 F α ( β + β(α 1) + (α + 1)p + β + α 1 ) (α 1) α 19
References [1] A Alpers, Instablty and stablty n dscrete tomography, PhD thess, Technsche Unverstät München, ISBN 3-83-355-X, Shaker Verlag, Aachen (003) [] A Alpers, S Brunett, Stablty results for the reconstructon of bnary pctures from two projectons, Image and Vson Computng 5 (007) 1599-1608 [3] A Alpers, P Grtzmann, L Thorens, Stablty and nstablty n dscrete tomography, Lectures Notes n Computer Scence 3: Dgtal and Image Geometry (001) 175-186 [] RJ Gardner, P Grtzmann, D Prangenberg, On the computatonal complexty of reconstructng lattce sets from ther X-rays, Dscrete Mathematcs 0 (1999) 5-71 [5] GT Herman, A Kuba, Dscrete tomography n medcal magng, Proceedngs of IEEE 91 (003) 161-166 [6] GT Herman, A Kuba, edtors, Dscrete Tomography: Foundatons, Algorthms and Applcatons, Brkhäuser, Boston (1999) [7] GT Herman, A Kuba, edtors, Advances n Dscrete Tomography and Its Applcatons, Brkhäuser, Boston (007) [8] JR Jnschek, KJ Batenburg, HA Calderon, R Klaas, V Radmlovc and C Kselowsk, 3-D reconstructon of the atomc postons n a smulated gold nanocrystal based on dscrete tomography, Ultramcroscopy (007), do:101016/jultramc0071000 [9] A Kuba, L Rodek, Z Kss, L Ruskó, A Nagy, M Balaskó, Dscrete tomography n neutron radography, Nucl Instrum Methods Phys Res, Sect A 5 (005) 376-38 [10] JC Palacos, LC Longora, J Santos, RT Perry, A PC-based dscrete tomography magng software system for assayng radoactve waste contaners, Nucl Instrum Methods Phys Res, Sect A 508 (003) 500-511 [11] HJ Ryser, Combnatoral propertes of matrces of zeros and ones, Can J Math 9 (1957) 371-377 [1] B Schllnger, Proposed combnaton of CAD data and dscrete tomography for the detecton of cokng and lubrcants n turbne blades or engnes, Electronc Notes n Dscrete Mathematcs 0 (005) 93-99 [13] H Slump, JJ Gerbrands, A network flow approach to reconstructon of the left ventrcle from two projectons, Computer graphcs and mage processng 18 (198) 18-36 [1] YR Wang, Characterzaton of bnary patterns and ther projectons, IEEE Trans Comput (1975) 103-1035 0
[15] L Wang, JY Park, Y Fu, Representaton of real partcles for DEM smulaton usng X-ray tomography, Constructon and Buldng Materals 1 (007) 338-36 1