Raw to Dicom conversion

Hello

I have written a code to convert .raw medical images to .dcm format. But I am having some error while trying to add some metadata information.

Following is my code:

fid3 = fopen( 'F:\patient\same patientCT.raw','r');

pCT = fread(fid3, [512*512*160,1], 'int16');

pCT_p = reshape(pCT, [512 512 160]);

fclose(fid3);

new_pCT1=permute(pCT_p,[2 1 3]);

for j =1:size(new_pCT1,3)

      File=int16(new_pCT1(:,:,j)+1000);

      m = -54.7500 -(j-1)*1.25;

      metadata.ImagePositionPatient=[-325 ; -325; m];

      metadata.SliceLocation=m;

      metadata.InStackPositionNumber=j;

      metadata.InstanceNumber=j;

      dicomwrite(File,['CT_' int2str(j) '.dcm'], metadata, 'CreateMode', 'copy');

     end

Please help me to resolve this issue. Thankyou

[Amy]

Hi Hamza, could you please share the exact error messages you are receiving when trying to run this code?

[Hamza Naeem]

Following are the errors:

Error using dicom_copy_IOD>getIOD (line 138)
Missing required attribute (0008,0016) "SOPClassUID"

Error in dicom_copy_IOD (line 26)
IOD_UID = getIOD(metadata, options, dictionary);

Error in dicomwrite>write_message (line 277)
        [attrs, status] = dicom_copy_IOD(X, map, ...

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

Error in Untitled (line 18)
    dicomwrite(File,['pCT' int2str(j) '.dcm'], metadata, 'CreateMode' , 'copy');

[Walter Roberson]

Try passing the appropriate ObjectType parameter to the dicomwrite() call.

[Hamza Naeem]

Can you please give some example how to do?

[Walter Roberson]

      dicomwrite(File,['pCT' int2str(j) '.dcm'], metadata, 'CreateMode' , 'copy', 'ObjectType', 'CT Image Storage')

https://www.mathworks.com/help/images/ref/dicomwrite.html

[Hamza Naeem]

Still not working, and now I am having this kind of error:

Warning: The 'SOPClassUID' and 'ObjectType' parameters are ignored when 'CreateMode' is 'copy'. 
> In dicomwrite>checkOptionConsistency (line 947)
  In dicomwrite (line 193)
  In Untitled (line 19) 
Error using dicom_copy_IOD>getIOD (line 138)
Missing required attribute (0008,0016) "SOPClassUID"

Error in dicom_copy_IOD (line 26)
IOD_UID = getIOD(metadata, options, dictionary);

Error in dicomwrite>write_message (line 277)
        [attrs, status] = dicom_copy_IOD(X, map, ...

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

Error in Untitled (line 19)
   dicomwrite(File,['pCT' int2str(j) '.dcm'], metadata, 'CreateMode' , 'copy', 'ObjectType', 'CT Image Storage')

[Walter Roberson]

As you are creating a new DICOM object, you will need to set a number of attributes in the metadata. See http://dicomlookup.com/lookup.asp?sw=Ttable&q=C.12-1 and http://dicom.nema.org/dicom/2013/output/chtml/part04/sect_B.5.html and https://dicom.innolitics.com/ciods/rt-plan/sop-common/00080016 and http://dicomiseasy.blogspot.ca/2011/12/chapter-4-dicom-objects-in-chapter-3.html

For example for plain CT image storage you might perhaps use SOPClassUID of 1.2.840.10008.5.1.4.1.1.2

There is an example at http://dicom.nema.org/DICOM/2013/output/chtml/part17/sect_LLL.2.html

... In short, if you are going to create your own metadata then you have a lot to create. It is usually far far easier to take another existing DICOM image of the exact same type that you want to create, and copy its metadata (setting different patient information and dates and a unique instance identifier).