metadatapydicomniftimedical-imaging

DICOM to Nifti metadata not transfering


I am trying to take a number of DICOM stacks and convert them to Nifti files. When I do the conversion and open the new Nifti file in a 3D viewer the volume is smashed together in the z direction. The Nifti files do not know what the spacing is between slices. To my understanding imageio.volread() does not read the metadata. I tried using pydicom.filereader.dcmread() but that only reads one file. How can I copy the metadata from the DICOM stack to the Nifti file when converting formats?

import nibabel as nib
import imageio
import numpy as np
import os, sys

DIR = '\\all scans\\'
savefold = '\\nifti\\'

for root, dirs, files in os.walk(DIR):
    for directory in dirs:
        vol = imageio.volread(DIR + directory).astype(int)
        vol = np.transpose(vol, (2,1,0)).astype(int)
        niftisave = nib.Nifti1Image(vol, affine=np.eye(4))
        nib.save(niftisave, os.path.join(savefold + directory) + '.nii')

UPDATE:

I am using Nifti1Header and setting my voxel spacing but the voxel spacing is still 1x1x1 when I save and open the file in other programs. When I print the header right before saving the pixdim shows [1. 0.09 0.09 0.09 1. 1. 1. 1. ].

header = nib.Nifti1Header()
OM = np.eye(4)
header.set_data_shape((224,352,224))
voxel_spacing = ((.09,.09,.09))
header.set_zooms(voxel_spacing)
header.set_sform(OM)
header.set_dim_info(slice = 2)

vol=imageio.volread(source) 

ROI_save = nib.Nifti1Image(vol, OM, header=header)
print(ROI_save.header)

HEADER:

<class 'nibabel.nifti1.Nifti1Header'> object, endian='<'
sizeof_hdr      : 348
data_type       : b''
db_name         : b''
extents         : 0
session_error   : 0
regular         : b''
dim_info        : 48
dim             : [  3 224 352 224   1   1   1   1]
intent_p1       : 0.0
intent_p2       : 0.0
intent_p3       : 0.0
intent_code     : none
datatype        : float32
bitpix          : 32
slice_start     : 0
pixdim          : [1.   0.09 0.09 0.09 1.   1.   1.   1.  ]
vox_offset      : 0.0
scl_slope       : nan
scl_inter       : nan
slice_end       : 0
slice_code      : unknown
xyzt_units      : 0
cal_max         : 0.0
cal_min         : 0.0
slice_duration  : 0.0
toffset         : 0.0
glmax           : 0
glmin           : 0
descrip         : b''
aux_file        : b''
qform_code      : unknown
sform_code      : aligned
quatern_b       : 0.0
quatern_c       : 0.0
quatern_d       : 0.0
qoffset_x       : 0.0
qoffset_y       : 0.0
qoffset_z       : 0.0
srow_x          : [1. 0. 0. 0.]
srow_y          : [0. 1. 0. 0.]
srow_z          : [0. 0. 1. 0.]
intent_name     : b''
magic           : b'n+1'

AFFINE:

np.eye(4)
--->[[1. 0. 0. 0.]
     [0. 1. 0. 0.]
     [0. 0. 1. 0.]
     [0. 0. 0. 1.]]

DESIRED AFFINE:

[[-0.09  0.    0.   -0.  ]
 [ 0.   -0.09  0.   -0.  ]
 [ 0.    0.    0.09  0.  ]
 [ 0.    0.    0.    1.  ]]

Solution

  • You need to directly specify pixel spacing and array shape, consider you have 512x512x128 3D volume, with 0.5 x 0.5 x 2.5 mm voxel spacing, and identity orientation matrix see example below:

    from nibabel import Nifti1Header, Nifti1Image
    
    img_array = np.zeros((512, 512, 128))
    
    voxel_spacing = [0.5, 0.5, 2.5, 1]
    OM = np.eye(4)
    OM = OM * np.diag(voxel_spacing)
    
    header = Nifti1Header()
    header.set_data_shape((512, 512, 128))
    header.set_dim_info(slice=2)
    header.set_xyzt_units('mm')
    
    nifti =  Nifti1Image(img_array, OM, header=header)
    

    upd. Save file using nibabel.save (or img.to_filename) and open it in MRIcron https://people.cas.sc.edu/rorden/mricron/index.html , gives the following result:

    enter image description here