Homec4science
Diffusion iCAPs public (master) Browse

iCAPs/packages/spm12/spm_dicom_essentials.m
31022b3a13e6master

Last Change
Raw File

spm_dicom_essentials.m
View Options

function hdr1 = spm_dicom_essentials(hdr0)
% Remove unused fields from DICOM header
% FORMAT hdr1 = spm_dicom_essentials(hdr0)
% hdr0 - original DICOM header
% hdr1 - Stripped down DICOM header
%
% With lots of DICOM files, the size of all the headers can become too
% big for all the fields to be saved. The idea here is to strip down
% the headers to their essentials.
%__________________________________________________________________________
% Copyright (C) 2008-2013 Wellcome Trust Centre for Neuroimaging
% John Ashburner
% $Id: spm_dicom_essentials.m 6588 2015-11-03 15:26:38Z john $
used_fields = {...
'AcquisitionDate',...
'AcquisitionNumber',...
'AcquisitionTime',...
'BitsAllocated',...
'BitsStored',...
'CSAImageHeaderInfo',...
'Columns',...
'EchoNumbers',...
'EchoTime',...
'Filename',...
'FlipAngle',...
'HighBit',...
'ImageOrientationPatient',...
'ImagePositionPatient',...
'ImageType',...
'InstanceNumber',...
'MagneticFieldStrength',...
'Manufacturer',...
'Modality',...
'MRAcquisitionType',...
'NumberOfFrames',...
'PatientID',...
'PatientsName',...
'PixelRepresentation',...
'PixelSpacing',...
'Private_0029_1110',...
'Private_0029_1210',...
'ProtocolName',...
'RepetitionTime',...
'RescaleIntercept',...
'RescaleSlope',...
'Rows',...
'SOPClassUID',...
'SamplesPerPixel',...
'ScanningSequence',...
'SequenceName',...
'SeriesDescription',...
'SeriesInstanceUID',...
'SeriesNumber',...
'SizeOfPixelData'...
'SliceNormalVector',...
'SliceThickness',...
'SpacingBetweenSlices',...
'StartOfPixelData',...
'StudyDate',...
'StudyTime',...
'TransferSyntaxUID',...
'VROfPixelData',...
'ScanOptions',...
'GE_ImageType',...
...
'PerFrameFunctionalGroupsSequence',...
'SharedFunctionalGroupsSequence',...
'DimensionIndexSequence',...
};
% Should probably save some memory by stripping out unused fields
% from PerFrameFunctionalGroupsSequence and PerFrameFunctionalGroupsSequence.
% This is not yet a priority though.
fnames = fieldnames(hdr0);
for i=1:numel(used_fields)
if ismember(used_fields{i},fnames)
hdr1.(used_fields{i}) = hdr0.(used_fields{i});
end
end
Private_spectroscopy_fields = {...
'Columns',...
'Rows',...
'ImageOrientationPatient',...
'ImagePositionPatient',...
'SliceThickness',...
'PixelSpacing',...
'VoiPhaseFoV',...
'VoiReadoutFoV',...
'VoiThickness'};
if isfield(hdr1,'Private_0029_1110')
hdr1.Private_0029_1110 = ...
getfields(hdr1.Private_0029_1110,...
Private_spectroscopy_fields);
end
if isfield(hdr1,'Private_0029_1210')
hdr1.Private_0029_1210 = ...
getfields(hdr1.Private_0029_1210,...
Private_spectroscopy_fields);
end
if isfield(hdr1,'CSAImageHeaderInfo')
CSAImageHeaderInfo_fields = {...
'SliceNormalVector',...
'NumberOfImagesInMosaic',...
'AcquisitionMatrixText',...
'ICE_Dims'};
hdr1.CSAImageHeaderInfo = ...
getfields(hdr1.CSAImageHeaderInfo,...
CSAImageHeaderInfo_fields);
end
if isfield(hdr1,'CSASeriesHeaderInfo')
CSASeriesHeaderInfo_fields = {};
hdr1.CSASeriesHeaderInfo = ...
getfields(hdr1.CSASeriesHeaderInfo,...
CSASeriesHeaderInfo_fields);
end
function str1 = getfields(str0,names)
str1 = [];
for i=1:numel(names)
for j=1:numel(str0)
if strcmp(str0(j).name,names{i})
str1 = [str1,str0(j)];
end
end
end
c4science · Help