Rotate Image Orientation(Patient) in DICOM -


i have extracted 3d surface mri acquisition , coordinates of points describing surface (i believe) respect reference system of first image of series (i mean origin corresponds image position(patient) , axes orientation image orientation(patient)). have set of images different image position(patient) , different image orientation(patient); want rotate , translate surface extracted first set in order have match second set of images.

i'm having trouble finding correct 4x4 matrix job, once it, know how apply surface.

any kind of appreciated, thank you.

simon

this page explains how form transformation matrix geometry information in dicom headers. these transformation matrices used transform volume coordinate system (pixel-x, pixel-y, slice number) patient/world coordinate system (in millimeters).

the basic idea transform volume 1 volume 2 tranform volume 1 patient coordinates , patient coordinates volume 2 coordinate system. multiplying both matrices yields matrix transform directly volume 1 volume 2.

caution: obviously, there can no guarantee every coordinate in v1 matches coordinate in v2, i.e. stacks may have different size and/or position.

so have:

m1 - matrix transform volume 1 world coordinate system , m2 - matrix transform volume 2 world coordinate system

then

m1 * (m2^(-1)) matrix transform position vector volume 1 volume 2 (input , output pixel-x, pixel-y, slice number)

and

m2 * (m1^(-1)) matrix transform position vector volume 1 volume 2 (input , output pixel-x, pixel-y, slice number)


Comments

Popular posts from this blog

PHP and MySQL WP -

android - InAppBilling registering BroadcastReceiver in AndroidManifest -

go - golang pprof for c library code -