MATLAB: Slice thickness information lost when using dicomwrite()

dicom filesdicominfodicomreaddicomwriteImage Processing Toolboxslice thickness

Hi, I have written a bit of code to perform some manipulations on a set of DICOM files whose slice thickness (according to dicominfo) is 1 and then write them using dicomwrite. When I use InfoDcm.Slicethickness=1; in order for the slice thickness to be preserved, it only appears as [] when I apply dicominfo to the output file. Does anyone know how to ensure the slice thickness is unchanged?

Thanks, RH

Best Answer

I think your problem may be the missing capitalization of the "t" in the name SliceThickness. Notice that all the fields in the structure returned by DICOMINFO are in upper camel case

Related Solutions

MATLAB: Efficiently saving multislice tiff using Tifflib

Does this help?

MATLAB: Extract rescale slope value

for K = 135 : -1 : 1
  petname = sprintf('PET_I1001_PT%03d.dcm', K);
  info(K) = dicominfo(petname);
end
rescale_slopes = [info.RescaleSlope];

If, for some reason you need to find the non-zero entries, then

slice_idx = find(rescale_slopes);

Note: the above code has a limitation that the dicominfo returned by each slice must have exactly the same set of fields. If that assumption is violated then the info(K) assignment will give you an error about assignment between dissimilar structures.

Best Answer

Related Solutions

MATLAB: Efficiently saving multislice tiff using Tifflib

MATLAB: Extract rescale slope value

Related Question