Lung Segmentation in 4D CT Volumes Based on Robust Active Shape Model Matching

<div>Comparison of results of methods <svg height="11.8174pt" id="M219" style="vertical-align:-3.1815pt" version="1.1" viewbox="-0.0498162 -8.6359 24.1077 11.8174" width="24.1077pt" 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="M962 650H795L470 145L347 650H176L170 622C268 613 275 606 244 503L190 322C150 188 132 126 118 91C102 50 80 33 18 28L12 0H245L251 28C175 35 170 48 174 93C177 128 191 180 220 284L292 542H294C331 392 383 150 409 4H432L774 555H776L714 137C700 40 694 34 612 28L606 0H868L874 28C793 34 784 37 797 137L849 533C859 612 863 616 956 622L962 650Z" id="g113-78"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="985" vert-adv-y="985"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,12.22,3.132)"><path d="M462 177V227H365V632H320C217 496 116 350 21 208V177H284V109C284 43 280 38 190 31V0H451V31C369 38 365 43 365 108V177H462ZM284 227H88C155 336 218 430 282 519H284V227Z" id="g50-53"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,16.652,3.132)"><path d="M43 650V617C122 612 130 606 130 522V128C130 44 122 37 34 32V0H275C400 0 497 25 573 82C652 142 698 229 698 350C698 445 662 518 613 565C549 625 456 650 327 650H43ZM217 544C217 580 220 591 229 599C240 607 264 613 311 613C380 613 437 599 482 570C561 523 598 436 598 336C598 177 511 40 324 40C248 40 217 61 217 135V544Z" id="g58-66"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="744" vert-adv-y="744"></glyph.data></g></svg>, <svg height="11.927pt" id="M220" style="vertical-align:-3.291101pt" version="1.1" viewbox="-0.0498162 -8.6359 24.1077 11.927" width="24.1077pt" 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="M962 650H795L470 145L347 650H176L170 622C268 613 275 606 244 503L190 322C150 188 132 126 118 91C102 50 80 33 18 28L12 0H245L251 28C175 35 170 48 174 93C177 128 191 180 220 284L292 542H294C331 392 383 150 409 4H432L774 555H776L714 137C700 40 694 34 612 28L606 0H868L874 28C793 34 784 37 797 137L849 533C859 612 863 616 956 622L962 650Z" id="g113-78"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="985" vert-adv-y="985"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,12.22,3.132)"><path d="M290 377C321 398 342 415 358 430C378 450 389 473 389 502C389 578 329 635 238 635H237C184 635 137 610 109 578L64 515L88 493C112 529 154 573 208 573S303 542 303 482C303 409 233 370 141 341L149 308C165 313 190 319 215 319C272 319 341 283 341 193C342 98 292 43 222 43C163 43 122 72 96 94C88 101 79 100 70 94C61 87 47 73 46 60C44 47 48 37 62 23C76 10 118 -12 165 -12C238 -12 430 62 430 223C430 297 379 359 290 375V377Z" id="g50-52"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,16.652,3.132)"><path d="M43 650V617C122 612 130 606 130 522V128C130 44 122 37 34 32V0H275C400 0 497 25 573 82C652 142 698 229 698 350C698 445 662 518 613 565C549 625 456 650 327 650H43ZM217 544C217 580 220 591 229 599C240 607 264 613 311 613C380 613 437 599 482 570C561 523 598 436 598 336C598 177 511 40 324 40C248 40 217 61 217 135V544Z" id="g58-66"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="744" vert-adv-y="744"></glyph.data></g></svg>, <svg height="14.1238pt" id="M221" style="vertical-align:-5.487901pt" version="1.1" viewbox="-0.0498162 -8.6359 47.6918 14.1238" width="47.6918pt" 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="M962 650H795L470 145L347 650H176L170 622C268 613 275 606 244 503L190 322C150 188 132 126 118 91C102 50 80 33 18 28L12 0H245L251 28C175 35 170 48 174 93C177 128 191 180 220 284L292 542H294C331 392 383 150 409 4H432L774 555H776L714 137C700 40 694 34 612 28L606 0H868L874 28C793 34 784 37 797 137L849 533C859 612 863 616 956 622L962 650Z" id="g113-78"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="985" vert-adv-y="985"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,12.22,3.132)"><path d="M462 177V227H365V632H320C217 496 116 350 21 208V177H284V109C284 43 280 38 190 31V0H451V31C369 38 365 43 365 108V177H462ZM284 227H88C155 336 218 430 282 519H284V227Z" id="g50-53"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,16.652,3.132)"><path d="M43 650V617C122 612 130 606 130 522V128C130 44 122 37 34 32V0H275C400 0 497 25 573 82C652 142 698 229 698 350C698 445 662 518 613 565C549 625 456 650 327 650H43ZM217 544C217 580 220 591 229 599C240 607 264 613 311 613C380 613 437 599 482 570C561 523 598 436 598 336C598 177 511 40 324 40C248 40 217 61 217 135V544Z" id="g58-66"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="744" vert-adv-y="744"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,23.422,3.132)"><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(.0091,0,0,-0.0091,26.771,3.132)"><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,30.72,3.132)"><path d="M472 439C432 432 382 428 334 425C303 442 271 451 236 451H235C158 451 54 398 54 285C54 216 97 171 140 152C124 133 93 105 53 90C52 81 55 63 64 48C77 27 102 5 137 -6C115 -22 77 -53 57 -73C40 -88 32 -107 32 -130C32 -189 96 -253 208 -253C340 -253 461 -159 461 -56C461 44 378 63 314 63C281 63 246 62 207 62C163 62 145 80 145 99C145 114 160 130 175 141C191 138 210 136 226 136C314 136 405 188 405 294C405 329 393 359 376 382L433 378C448 387 469 415 477 430L472 439ZM226 417C282 417 320 364 320 288S281 172 236 171C180 171 142 224 142 300S182 417 226 417ZM248 -9C291 -9 319 -11 346 -22S390 -57 390 -90C390 -151 339 -199 246 -199C171 -199 114 -159 114 -106C114 -79 132 -53 159 -30C178 -15 200 -9 248 -9Z" id="g58-101"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="479" vert-adv-y="479"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,35.079,3.132)"><path d="M682 0V33C615 38 603 49 571 134C507 305 441 494 378 665L343 655L136 136C101 48 86 41 22 33V0H246V33C170 41 161 50 183 112C196 152 209 192 224 232H441C461 176 478 126 493 87C504 49 499 41 433 33V0H682ZM425 279H242C273 363 303 448 335 533H337L425 279Z" id="g58-63"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="703" vert-adv-y="703"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,41.476,3.132)"><path d="M384 354C435 370 515 408 515 502C515 549 494 593 459 616C426 639 384 650 308 650H45V618C127 613 130 606 130 524V127C130 45 125 38 37 33V0H262C342 0 409 12 461 40C521 73 558 124 558 192C558 295 481 336 384 352V354ZM217 365V556C217 587 221 597 229 604C238 611 261 615 288 615C382 615 421 558 421 490C421 406 374 365 266 365H217ZM217 327H264C384 327 457 284 457 183C457 80 383 38 301 38C235 38 217 55 217 127V327Z" id="g58-64"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="598" vert-adv-y="598"></glyph.data></g></svg>, and <svg height="14.1238pt" id="M222" style="vertical-align:-5.487901pt" version="1.1" viewbox="-0.0498162 -8.6359 47.4726 14.1238" width="47.4726pt" 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="M962 650H795L470 145L347 650H176L170 622C268 613 275 606 244 503L190 322C150 188 132 126 118 91C102 50 80 33 18 28L12 0H245L251 28C175 35 170 48 174 93C177 128 191 180 220 284L292 542H294C331 392 383 150 409 4H432L774 555H776L714 137C700 40 694 34 612 28L606 0H868L874 28C793 34 784 37 797 137L849 533C859 612 863 616 956 622L962 650Z" id="g113-78"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="985" vert-adv-y="985"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,12.22,3.132)"><path d="M462 177V227H365V632H320C217 496 116 350 21 208V177H284V109C284 43 280 38 190 31V0H451V31C369 38 365 43 365 108V177H462ZM284 227H88C155 336 218 430 282 519H284V227Z" id="g50-53"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,16.652,3.132)"><path d="M43 650V617C122 612 130 606 130 522V128C130 44 122 37 34 32V0H275C400 0 497 25 573 82C652 142 698 229 698 350C698 445 662 518 613 565C549 625 456 650 327 650H43ZM217 544C217 580 220 591 229 599C240 607 264 613 311 613C380 613 437 599 482 570C561 523 598 436 598 336C598 177 511 40 324 40C248 40 217 61 217 135V544Z" id="g58-66"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="744" vert-adv-y="744"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,23.422,3.132)"><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(.0091,0,0,-0.0091,26.771,3.132)"><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,30.72,3.132)"><path d="M472 439C432 432 382 428 334 425C303 442 271 451 236 451H235C158 451 54 398 54 285C54 216 97 171 140 152C124 133 93 105 53 90C52 81 55 63 64 48C77 27 102 5 137 -6C115 -22 77 -53 57 -73C40 -88 32 -107 32 -130C32 -189 96 -253 208 -253C340 -253 461 -159 461 -56C461 44 378 63 314 63C281 63 246 62 207 62C163 62 145 80 145 99C145 114 160 130 175 141C191 138 210 136 226 136C314 136 405 188 405 294C405 329 393 359 376 382L433 378C448 387 469 415 477 430L472 439ZM226 417C282 417 320 364 320 288S281 172 236 171C180 171 142 224 142 300S182 417 226 417ZM248 -9C291 -9 319 -11 346 -22S390 -57 390 -90C390 -151 339 -199 246 -199C171 -199 114 -159 114 -106C114 -79 132 -53 159 -30C178 -15 200 -9 248 -9Z" id="g58-101"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="479" vert-adv-y="479"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,35.079,3.132)"><path d="M384 354C435 370 515 408 515 502C515 549 494 593 459 616C426 639 384 650 308 650H45V618C127 613 130 606 130 524V127C130 45 125 38 37 33V0H262C342 0 409 12 461 40C521 73 558 124 558 192C558 295 481 336 384 352V354ZM217 365V556C217 587 221 597 229 604C238 611 261 615 288 615C382 615 421 558 421 490C421 406 374 365 266 365H217ZM217 327H264C384 327 457 284 457 183C457 80 383 38 301 38C235 38 217 55 217 127V327Z" id="g58-64"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="598" vert-adv-y="598"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,40.302,3.132)"><path d="M682 0V33C615 38 603 49 571 134C507 305 441 494 378 665L343 655L136 136C101 48 86 41 22 33V0H246V33C170 41 161 50 183 112C196 152 209 192 224 232H441C461 176 478 126 493 87C504 49 499 41 433 33V0H682ZM425 279H242C273 363 303 448 335 533H337L425 279Z" id="g58-63"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="703" vert-adv-y="703"></glyph.data></g></svg> on a longitudinal CT volume pair. The time between scans A and B was 19 months, and the imaging protocol was quite different.</div>

International Journal of Biomedical Imaging

fig8

Figure 8

Figure 8: Lung Segmentation in 4D CT Volumes Based on Robust Active Shape Model Matching 