how to convert .mat to .dcm

Hi I have a 128x128 matrix I want to convert it to .dcm file What matlab comment I should use

Answer 1

[Walter Roberson]

http://www.mathworks.com/help/images/write-data-to-dicom-files.html

[Meggie Hwang]

thank you!

[Hany]

Hi there I am doing dicomwrite(Reconimg.mat, 'Reconimg.dcm'); but it does not work , could u please help?

[Walter Roberson]

Hany,

You need to read a variable from the .mat file and dicomwrite the variable

[Hany]

Thanks Walter , the mat file includes 320 images , is there any way to convert them all to a series of DCM images? Thanks

[Walter Roberson]

image_dir = 'dicom_images';    %directory to store the dicom images in

datastruct = load('NameOfMat.mat');
image_names = fieldnames(datastruct);

for K = 1 : length(image_names)
  this_image = image_names{K};
  this_data = datastruct.(this_image);
  outfile = fullfile( image_dir, [this_image, '.dcm']);
  dicomwrite( this_data, outfile );
end

[Hany]

Thanks for your answer , i get the following errors:

image_dir'C:\Users\Hany\Dropbox\Hany\PhD\Patients\dicom_images';   
datastruct = load('Reconimg.mat');
image_names = fieldnames(datastruct);
for K = 1 : length(image_names)
this_image = image_names{K};
this_data = datastruct.(this_image);
outfile = fullfile( image_dir, [this_image, '.dcm']);
  dicomwrite( this_data, outfile );
end

Error using dicom_prep_ImagePixel>getPhotometricInterp (line 116)
Cannot determine photometric interpretation.

Error in dicom_prep_ImagePixel (line 10)
metadata.(dicom_name_lookup('0028', '0004', dictionary)) = getPhotometricInterp(X, map,
txfr);

Error in dicom_prep_metadata (line 52)
    metadata = dicom_prep_ImagePixel(metadata, X, map, txfr, dictionary);

Error in dicom_create_IOD (line 27)
metadata = dicom_prep_metadata(IOD_UID, metadata, X, map, options.txfr, dictionary);

Error in dicomwrite>write_message (line 263)
        [attrs, status] = dicom_create_IOD(SOP_UID, X, map, ...

Error in dicomwrite (line 200)
    [status, options] = write_message(X, filename, map, metadata, options);

[Walter Roberson]

It appears that instead of being in the form of one variable per image, that you have a single 3D array that is to be split up. In that case, the approach you used in your newer Question is appropriate.

[Hany]

Thanks Walter fopr your answer :)

[Somayeh Asadi]

Dear Hany

I wanted to do the same that did to convert a matrix to dicom image. I followed the instructer of Walter Roberson and took the same error:

Error using dicom_prep_ImagePixel>getPhotometricInterp (line 116)

Cannot determine photometric interpretation.

Error in dicom_prep_ImagePixel (line 10)

metadata.(dicom_name_lookup('0028', '0004', dictionary)) = getPhotometricInterp(X, map,

txfr);

Error in dicom_prep_metadata (line 52)

metadata = dicom_prep_ImagePixel(metadata, X, map, txfr, dictionary);

Error in dicom_create_IOD (line 27)

metadata = dicom_prep_metadata(IOD_UID, metadata, X, map, options.txfr, dictionary);

Error in dicomwrite>write_message (line 263)

[attrs, status] = dicom_create_IOD(SOP_UID, X, map, ...

Error in dicomwrite (line 200)

[status, options] = write_message(X, filename, map, metadata, options);

Could you please giude me how to solve it?

Thank you in advance