Create a TFORM structure for SPATIAL information given by dicomreadvolume

Anders

2020 11 30

0 个回答

3 次查看（30 天）

0 个投票

I have collected a number of MRI volumes in different orientations cor,sag,tra and also oblique volume. All image are from the same examination and I now like to compare them. Some of the image volumes are on homogenous spherical phantoms so there is not possible to do image based registration to get the TFORM.

Can i calculate a 3D TFORM structure from the SPATIAL information given by dicmreadvolume?

I.e. what I like to do is (for simplicity assume that the volumes have the same FOV and resolution)

    [SagVolume, sagSPATIAL, sagDIM] = dicomreadvolume(sagVolumeDir);
    
    [ObliqueVolume, obliqueSPATIAL, obliqueDIM] = dicomreadvolume(obliqueVolumeDir);  
    
    tformSag2standardSpace = ???? something dependent on sagSPATIAL.
    
    tformOblique2standardSpce  = ???? something dependent on obliqueSPATIAL
    
    SagVolumeInStandradSpcae = affine3D( SagVolume, tformSag2standardSpace);
    ObliqueVolumeInStandradSpcae = affine3D( ObliqueVolume, tformOblique2standardSpace);
    
    DiffVolumes = SagVolumeInStandradSpcae - ObliqueVolumeInStandradSpcae;