Polychromatic Iterative Statistical Material Image Reconstruction for Photon-Counting Computed Tomography

<div>Normalized deviation of the iodine ((a), (b)) and water ((c), (d)) fraction from the ground truth (<svg height="15.6876pt" id="M249" style="vertical-align:-3.4294pt" version="1.1" viewbox="-0.0498162 -12.2582 89.66 15.6876" width="89.66pt" 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="M600 0V24L339 665L311 657L43 24V0H600ZM497 50H107L303 539L497 50Z" id="g113-133"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="641" vert-adv-y="641"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,7.153,0)"><path d="M619 670C619 686 593 712 555 712S459 686 410 634S335 504 320 430H250L219 400L222 388H312L258 73C223 -133 201 -166 187 -180C175 -191 158 -199 140 -199C123 -199 88 -188 74 -172C68 -166 63 -164 54 -171C38 -185 23 -201 23 -215C23 -236 60 -261 93 -261C122 -261 161 -247 207 -200C268 -138 300 -71 337 94C365 220 376 277 394 387L501 399L521 430H401C432 623 464 665 501 665C524 665 544 651 567 627C577 617 583 618 592 625C601 631 619 651 619 670Z" id="g113-103"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="642" vert-adv-y="642"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,15.505,-5.741)"><path d="M138 535C169 535 191 560 191 589C191 620 168 644 139 644S85 620 85 589C85 560 110 535 138 535ZM257 0V30C192 35 185 43 185 109V453C142 435 94 422 40 414V386C101 378 104 372 104 311V109C104 41 97 35 33 30V0H257Z" id="g58-103"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="275" vert-adv-y="275"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,18.007,-5.741)"><path d="M262 451C169 451 38 375 38 210C38 98 121 -12 261 -12C360 -12 481 64 481 227C481 351 386 451 262 451ZM249 414C337 414 385 320 385 204C385 70 333 25 272 25C189 25 134 117 134 241C134 355 189 414 249 414Z" id="g58-109"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="519" vert-adv-y="519"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,22.812,-5.741)"><path d="M526 56L493 58C457 61 448 68 448 119V710C411 697 344 682 288 675V646C362 641 367 639 367 573V439C344 448 315 451 300 451C164 451 40 342 40 202C40 61 147 -12 227 -12C239 -12 266 -7 305 16L367 53V-12C427 10 505 21 526 26V56ZM367 88C342 69 306 54 271 54C207 54 133 111 133 229C133 373 219 409 264 409C301 409 342 394 367 359V88Z" id="g58-98"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="538" vert-adv-y="538"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,27.708,-5.741)"><path d="M138 535C169 535 191 560 191 589C191 620 168 644 139 644S85 620 85 589C85 560 110 535 138 535ZM257 0V30C192 35 185 43 185 109V453C142 435 94 422 40 414V386C101 378 104 372 104 311V109C104 41 97 35 33 30V0H257Z" id="g58-103"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="275" vert-adv-y="275"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,30.21,-5.741)"><path d="M532 0V30C473 38 467 38 467 107V298C467 395 416 451 336 451C307 451 279 439 251 420C226 404 205 389 185 374V453C142 434 94 422 43 413V386C98 378 104 373 104 309V107C104 41 100 38 31 30V0H251V30C191 35 185 40 185 107V337C216 364 255 389 295 389C360 389 386 347 386 276V110C386 44 380 35 321 30V0H532Z" id="g58-108"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="557" vert-adv-y="557"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,35.252,-5.741)"><path d="M384 109C349 76 310 59 268 59C197 59 123 113 120 245C240 248 367 259 384 262C403 266 407 275 407 296C407 376 338 451 254 451H253C201 451 148 425 105 378C65 334 39 272 39 203C39 88 112 -12 235 -12C266 -12 339 5 402 84L384 109ZM230 411C285 411 319 365 319 313C319 298 313 293 294 293C236 291 180 290 125 290C141 370 185 411 230 411Z" id="g58-99"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="434" vert-adv-y="434"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,39.201,-5.741)"><path d="M167 -163V704H102V-163H167Z" id="g50-9"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="269" vert-adv-y="269"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,41.649,-5.741)"><path d="M708 439H523V409C587 401 594 395 581 344C562 272 534 184 513 118H511C468 225 428 330 391 434H361C318 317 283 217 244 121H242C218 196 193 273 170 349C157 394 160 401 220 409V439H10V409C60 403 67 394 86 343C127 224 165 109 201 -10H233C273 101 314 199 353 299C393 195 432 99 475 -10H507C554 126 596 251 627 334C648 390 657 401 708 409V439Z" id="g58-117"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="698" vert-adv-y="698"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,48.092,-5.741)"><path d="M442 39L431 68C420 63 406 57 394 57C377 57 359 71 359 118V301C359 353 349 393 314 424C291 443 261 451 230 451C171 438 107 401 79 381C59 366 46 355 46 340C46 316 72 295 91 295C106 295 117 303 122 319C131 350 139 371 152 384C163 397 178 404 197 404C246 404 281 365 281 293V278C260 262 186 232 124 213C69 194 41 161 41 112C41 47 92 -12 163 -12C195 -12 243 25 283 51C288 35 293 22 308 7C319 -4 339 -12 354 -12L442 39ZM281 87C262 68 226 51 199 51C169 51 129 74 129 126C129 162 151 181 190 198C211 208 259 229 281 240V87Z" id="g58-95"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="449" vert-adv-y="449"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,52.014,-5.741)"><path d="M303 37L293 66C280 59 256 49 231 49S173 63 173 143V396H281C294 404 298 428 287 439H173V575L158 577L92 509V439H46L17 408L22 396H92V107C92 27 128 -12 192 -12C203 -12 218 -9 237 2L303 37Z" id="g58-114"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="312" vert-adv-y="312"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,54.807,-5.741)"><path d="M384 109C349 76 310 59 268 59C197 59 123 113 120 245C240 248 367 259 384 262C403 266 407 275 407 296C407 376 338 451 254 451H253C201 451 148 425 105 378C65 334 39 272 39 203C39 88 112 -12 235 -12C266 -12 339 5 402 84L384 109ZM230 411C285 411 319 365 319 313C319 298 313 293 294 293C236 291 180 290 125 290C141 370 185 411 230 411Z" id="g58-99"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="434" vert-adv-y="434"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,58.756,-5.741)"><path d="M185 344V453C137 433 92 421 41 413V386C100 378 104 375 104 310V106C104 41 97 35 33 30V0H268V30C190 35 185 41 185 107V283C206 340 240 371 267 371C282 371 296 364 310 351C316 344 325 343 336 348C357 357 369 378 369 399C369 424 345 451 310 451C262 451 219 397 187 344H185Z" id="g58-112"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="380" vert-adv-y="380"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,62.888,0)"><path d="M368 703H309L44 -163H104L368 703Z" id="g113-48"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="412" vert-adv-y="412"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,68.249,0)"><path d="M619 670C619 686 593 712 555 712S459 686 410 634S335 504 320 430H250L219 400L222 388H312L258 73C223 -133 201 -166 187 -180C175 -191 158 -199 140 -199C123 -199 88 -188 74 -172C68 -166 63 -164 54 -171C38 -185 23 -201 23 -215C23 -236 60 -261 93 -261C122 -261 161 -247 207 -200C268 -138 300 -71 337 94C365 220 376 277 394 387L501 399L521 430H401C432 623 464 665 501 665C524 665 544 651 567 627C577 617 583 618 592 625C601 631 619 651 619 670Z" id="g113-103"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="642" vert-adv-y="642"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,76.601,-5.741)"><path d="M700 305H441V272C542 264 548 258 548 186V107C548 64 536 53 516 42C493 30 462 24 429 24C237 24 151 189 151 333C151 517 264 626 415 626C507 626 585 594 609 470L642 476C634 544 629 601 626 636C594 644 513 665 429 665C234 665 46 554 46 321C46 123 194 -15 415 -15C495 -15 578 7 642 22C635 50 635 81 635 116V201C635 259 639 264 700 272V305Z" id="g58-69"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="725" vert-adv-y="725"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,83.199,-5.741)"><path d="M601 496C596 555 591 627 591 675H567C551 656 544 650 516 650H120C91 650 79 651 64 675H41C38 619 34 555 29 492H63C75 537 83 566 96 582C110 604 128 612 207 612H272V126C272 46 264 38 170 33V0H464V33C368 38 359 48 359 126V612H430C502 612 520 605 535 582C547 566 556 541 567 492L601 496Z" id="g58-82"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="628" vert-adv-y="628"></glyph.data></g></svg>), evaluated within the ROIs shown in Figures <a href="../fig4/#a">4(a)</a> and <a href="../fig4/#b">4(b)</a>, respectively.</div>

International Journal of Biomedical Imaging

fig6

Figure 6

Figure 6: Polychromatic Iterative Statistical Material Image Reconstruction for Photon-Counting Computed Tomography 