How to convert dataset of .mat to .tif files?

I have 100 hyperspectral images (501x502x38) in .mat format, and I want to convert them to .tif files, but my code does not work. Here is the code. I do appreciate your help.

fileList = dir(fullfile(inputFolder,'*.mat'));
for kk = 1:numel(fileList)
  S = load(fullfile(fileList(kk).folder,fileList(kk).name));
  I = S.cjdata.image;
  fileName = replace(fileList(kk).name,'.mat','.tiff');
  imwrite(I,fullfile(outputFolder,fileName));
end

[ANKUR KUMAR]

If you can attached couple of your matfiles (if size is not more than 50 MB), it would be good for us to help.

[Ronak]

Thank you, but unfortunetaly, the size is about 75 MB.

[Simon Allosserie]

It's difficult to assess without seeing the original files. Can you do

imshow(I)

Just to check if that is an actual image? If that already doesn't work, you first have to fix the way you form I.

Answer 1

[Walter Roberson]

You need to use the Tiff() class https://www.mathworks.com/help/matlab/ref/tiff.html

Use SamplesPerPixel = 38, and use ExtraSamples https://www.awaresystems.be/imaging/tiff/tifftags/extrasamples.html as a vector of (38-3)=35 EXTRASAMPLE_UNSPECIFIED values.

However, you need to know more about what you are going to do with the TIFF files. If you were planning to use them as an "Image Stack" in ImageJ then you need a different arrangement, as ImageJ implements Image Stacks through a private tag. If you were planning to use them as GeoTiff, then you should probably use the geotiff routines in MATLAB.