NUFFT-Based Iterative Image Reconstruction via Alternating Direction Total Variation Minimization for Sparse-View CT

<div>An object <svg height="12.6917pt" id="M9" style="vertical-align:-3.40337pt" version="1.1" viewbox="-0.0498162 -9.28833 38.228 12.6917" width="38.228pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M620 667C620 684 593 710 552 710C515 710 459 687 409 634C360 582 334 506 320 433H247L215 400L220 388H310L259 80C224 -129 199 -161 185 -174C175 -183 160 -192 141 -192C123 -192 90 -181 75 -165C69 -158 64 -156 54 -164C39 -178 24 -195 24 -208C24 -231 61 -257 94 -257C124 -257 162 -243 209 -194C267 -134 303 -66 343 107C369 220 378 279 393 385L503 397L523 433H404C432 619 464 660 500 660C522 660 542 646 566 621C573 614 583 615 593 621C603 628 620 648 620 667Z" id="g50-103"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="644" vert-adv-y="644"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,8.378,0)"><path d="M309 -142C211 -61 151 99 151 271S210 601 309 683L288 710C141 611 75 451 75 272V271C75 90 141 -72 288 -169L309 -142Z" id="g50-41"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="358" vert-adv-y="358"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,13.032,0)"><path d="M545 403C545 425 527 451 498 451C448 451 400 404 312 286L290 341C262 412 246 451 219 451C182 451 138 404 92 345L112 323C152 368 169 378 181 378C192 378 201 359 216 321L257 214C184 115 142 69 113 69C102 69 93 73 88 79S78 89 68 89C47 89 24 61 24 40C24 8 49 -12 75 -12C125 -12 175 31 271 175L312 64C330 15 357 -12 382 -12C421 -12 475 35 515 98L496 121C466 88 439 62 420 62C402 62 384 92 363 147L325 247C349 280 374 309 397 333C422 361 446 381 462 381C471 381 481 374 488 366C493 360 499 357 505 357C521 357 545 380 545 403Z" id="g50-121"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="569" vert-adv-y="569"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,20.429,0)"><path d="M98 134C72 134 46 117 46 90C46 73 55 65 60 64C95 55 124 32 124 -4C124 -42 95 -68 44 -89L57 -123C122 -104 194 -60 194 22C194 94 136 134 98 134Z" id="g50-45"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="237" vert-adv-y="237"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,25.716,0)"><path d="M563 395C563 428 543 451 519 451C500 451 480 438 472 423C468 415 465 403 470 394C477 381 482 368 482 347C482 272 386 136 342 74H340C334 163 323 257 312 337C302 410 290 451 257 451C214 451 160 385 127 338L145 310C169 341 203 373 209 373C219 373 225 369 232 325C249 215 269 62 271 -15C220 -80 135 -165 18 -231L26 -257L142 -230C254 -105 279 -74 342 13C394 85 484 211 526 286C553 334 563 365 563 395Z" id="g50-122"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="587" vert-adv-y="587"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,33.349,0)"><path d="M283 271C283 451 220 610 70 710L48 683C146 603 207 443 207 271S149 -59 48 -142L70 -169C220 -67 283 90 283 271Z" id="g50-42"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="358" vert-adv-y="358"></glyph.data></g></svg> and its parallel projection <svg height="12.6134pt" id="M10" style="vertical-align:-3.32507pt" version="1.1" viewbox="-0.0498162 -9.28833 26.9979 12.6134" width="26.9979pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M573 302C573 402 527 451 414 451C386 451 343 446 313 437L330 513L320 522C295 508 261 484 243 463L230 415C194 400 159 383 126 359L131 330C159 344 187 357 222 368L109 -147C96 -204 80 -214 18 -223L13 -255L256 -244L259 -212L236 -210C184 -205 180 -195 191 -141L219 -1C240 -10 268 -12 284 -12C352 4 431 48 484 104C543 166 573 240 573 302ZM481 290C481 165 381 37 305 37C280 37 249 56 235 71L302 395C328 399 353 402 372 402C427 402 481 378 481 290Z" id="g50-113"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="597" vert-adv-y="597"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,7.436,3.001)"><path d="M482 498C482 606 448 710 328 710C138 710 24 416 24 218C24 99 58 -12 184 -12C373 -12 482 294 482 498ZM395 521C395 483 391 449 383 395H128C158 537 226 671 312 671C388 671 395 566 395 521ZM375 342C348 190 284 25 192 25C128 25 109 116 109 199C109 243 112 292 121 342H375Z" id="g50-230"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="506" vert-adv-y="506"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,12.556,0)"><path d="M309 -142C211 -61 151 99 151 271S210 601 309 683L288 710C141 611 75 451 75 272V271C75 90 141 -72 288 -169L309 -142Z" id="g50-41"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="358" vert-adv-y="358"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,17.21,0)"><path d="M357 391C357 416 325 451 270 451C241 451 179 431 146 400C109 365 97 334 97 307C97 248 143 213 196 180C250 146 261 123 261 101C261 75 237 40 190 40C153 40 110 72 86 109C81 116 68 119 55 112C37 102 24 87 24 66C24 30 85 -12 137 -12C229 -12 331 62 331 140C331 189 304 218 236 260C197 284 164 309 164 346C164 381 194 400 220 400C259 400 282 380 304 350C310 342 321 339 330 344C347 353 357 372 357 391Z" id="g50-116"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="381" vert-adv-y="381"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,22.163,0)"><path d="M283 271C283 451 220 610 70 710L48 683C146 603 207 443 207 271S149 -59 48 -142L70 -169C220 -67 283 90 283 271Z" id="g50-42"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="358" vert-adv-y="358"></glyph.data></g></svg>.</div>

Computational and Mathematical Methods in Medicine

fig1

Figure 1

Figure 1: NUFFT-Based Iterative Image Reconstruction via Alternating Direction Total Variation Minimization for Sparse-View CT 