Diffusion ShoulderDB_ImportSegmentMeasureAnalyze (master) Browse

dicom_create_IOD.m
4b5355e68791master

dicom_create_IOD.m
View Options

			function [attrs, status] = dicom_create_IOD(IOD_UID, X, map, metadata, options)
			%DICOM_CREATE_IOD Create the attributes for a given IOD.
			% [ATTRS, STATUS] = DICOM_CREATE_IOD(UID, X, MAP, METADATA, TXFR)
			% creates a structure array of DICOM attributes for the SOP Class
			% corresponding to UID. The attributes' values are bassed on the image
			% X, colormap MAP, the METADATA struct (as given by DICOMINFO, for
			% example), and the transfer syntax UID (TXFR) used to encode the
			% file.
			%
			% This is the principal function of DICOM information object creation.
			%
			% See also DICOMWRITE, DICOM_PREP_METADATA, DICOM_IODS, DICOM_MODULES.

			% Copyright 1993-2014 The MathWorks, Inc.

			dictionary = dicomdict('get_current');

			% Find modules and other details.
			iod_details = get_iod_details(IOD_UID);

			if (isempty(iod_details))
			error(message('images:dicom_create_IOD:unimplementedIOD', IOD_UID));
			end

			% Set necessary values for this IOD.
			metadata = dicom_prep_metadata(IOD_UID, metadata, X, map, options.txfr, options.usemetadatabitdepths, dictionary);

			% Look for all of the private metadata.
			if (isfield(options, 'writeprivate'))

			if (islogical(options.writeprivate) \|\| ...
			isnumeric(options.writeprivate))

			if (options.writeprivate)

			[private_tags, private_names] = find_private_metadata(metadata, dictionary);

			end

			else

			error(message('images:dicom_create_IOD:invalidWritePrivateValue'))

			end

			end

			% A module definition function must be registered for unsupported or
			% private modules. Private IODs need only use a custom definition
			% function for unimplemented modules, but they must have a definition
			% function listed. (Registration not yet supported.)
			if (isempty(iod_details.Def_fcn))
			error(message('images:dicom_create_IOD:noModuleDefinitionFcn', iod_details.Name, IOD_UID));
			end

			attrs = [];
			status = [];
			for p = 1:numel(iod_details.Modules(:,1))

			% Determine whether to encode this module.
			if (test_module_condition(X, metadata, iod_details.Modules(p,:), dictionary))

			% Find the attributes in the module.
			module_details = dicom_modules(iod_details.Modules{p,1});

			if (isempty(module_details))
			error(message('images:dicom_create_IOD:undefinedModule', iod_details.Modules{ p, 1 }));
			end

			% Process attributes.
			[new_attrs, new_status] = process_modules(module_details.Attrs, ...
			metadata, ...
			dictionary);
			attrs = add_to_IOD(attrs, new_attrs);
			status = add_to_status(status, new_status);

			end

			end

			% Process any private metadata.
			if ((options.writeprivate) && (~isempty(private_tags)))

			new_attrs = process_private(private_tags, private_names, metadata, dictionary, options.txfr);
			attrs = add_to_IOD(attrs, new_attrs);

			end



			function iod_details = get_iod_details(IOD_UID)
			%GET_IOD_DETAILS Find the details of a given IOD.

			% Load IOD definitions.
			iod_directory = dicom_iods;

			% Look for a particular UID.
			idx = strmatch(IOD_UID, {iod_directory.UID}, 'exact');

			if (isempty(idx))
			iod_details = [];
			else
			iod_details = iod_directory(idx);
			end



			function attrs = add_to_IOD(attrs, new_attrs)
			%ADD_TO_IOD Add attributes to an existing IOD.
			attrs = cat(2, attrs, new_attrs);



			function status = add_to_status(status, new_status)
			%ADD_TO_STATUS Add status values to an existing status struct.

			fnames = fieldnames(new_status);
			for p = 1:numel(fnames)
			if (isfield(status, fnames{p}))
			status(1).(fnames{p}) = [status(1).(fnames{p}) ...
			new_status(1).(fnames{p})];
			else
			status(1).(fnames{p}) = new_status(1).(fnames{p});
			end
			end



			function [attrs, status] = process_modules(attr_details, metadata, dictionary)
			%PROCESS_MODULES Create attributes from a module definition.

			attrs = [];
			status.BadAttribute = {};
			status.MissingCondition = {};
			status.MissingData = {};
			status.SuspectAttribute = {};

			p = 0;
			numAttrs = numel(attr_details(:,1));
			% For each attribute in the module...
			while (p < numAttrs)

			p = p + 1;

			level = attr_details{p, 1};
			group = attr_details{p, 2};
			element = attr_details{p, 3};
			attr_type = attr_details{p, 4};
			vr_map = attr_details{p, 5};
			enum_values = attr_details{p, 6};
			condition = attr_details{p, 7};

			% Skip all sequence attributes for now.
			%
			% NOTE: At this point, all sequences in supported modules are
			% type 3 (optional). We can, and will, safely skip all
			% attributes in a sequence (including the top level attribute).

			nextAttr = p + 1;
			if ((p < numAttrs) && ...
			(attr_details{nextAttr, 1} > level))

			% The "level" determines nesting, with 0 as the
			% non-sequence elements
			while ((nextAttr <= numAttrs) && ...
			(attr_details{nextAttr, 1} > level))

			nextAttr = nextAttr + 1;

			end

			% Point to the next attribute after the sequence.
			p = nextAttr - 1;
			continue

			end

			%
			% Look for the attribute from the module.
			%
			try
			name = dicom_name_lookup(group, element, dictionary);
			catch
			status.BadAttribute{end + 1} = sprintf('(%s,%s)', ...
			group, element);
			continue;
			end

			if (isempty(name))
			status.BadAttribute{end + 1} = sprintf('(%s,%s)', ...
			group, element);
			continue;
			end

			%
			% Determine how to handle the attribute.
			%
			switch (attr_type)
			case {'1', '2', '3'}
			% Must process the attribute. No special action.

			case {'1C', '2C'}
			% Conditionally process the attribute.
			if (isempty(condition))

			status.MissingCondition{end + 1} = sprintf('(%s,%s)', ...
			group, element);
			continue;

			else

			if (~check_condition(condition, metadata, dictionary))
			continue;
			end

			end

			otherwise
			status.BadAttribute{end + 1} = sprintf('(%s,%s)', group, element);
			continue;

			end

			% Process special VR mapping.
			if ((isequal(group, '7FE0')) && (isequal(element, '0010')))
			VR = get_pixelData_VR(metadata, dictionary);
			else
			if (numel(vr_map) == 0)
			VR = '';
			else
			VR = remap_vr(group, element, vr_map, metadata, dictionary);
			end
			end

			%
			% Add the attribute.
			%
			if (isfield(metadata, name))

			% Process acceptable values.
			if ((~isempty(enum_values)) && (~isempty(metadata.(name))))
			if (~check_values(metadata.(name), enum_values))
			status.SuspectAttribute{end + 1} = sprintf('(%s,%s)', ...
			group, element);
			end
			end

			% Add the attribute.
			if (isempty(VR))
			attrs = dicom_add_attr(attrs, group, element, dictionary, ...
			metadata.(name));
			else
			attrs = dicom_add_attr(attrs, group, element, dictionary, ...
			metadata.(name), VR);
			end

			else

			switch (attr_type)
			case {'1', '1C'}

			% The attribute should exist.
			status.MissingData{end + 1} = sprintf('(%s,%s)', ...
			group, element);
			continue;

			case {'2', '2C'}

			% Nonexistent attributes have empty data.
			attrs = dicom_add_attr(attrs, group, element, dictionary);

			end

			% Nonexistent type 3 is ignored.

			end

			end



			function tf = check_condition(expr, metadata, dictionary)
			%CHECK_CONDITION Determine whether a particular condition is true.
			%
			% Conditions are LISP-style cell arrays. The first element is a
			% conditional operator, remaining cells are arguments to the operator.
			%
			% Conditions can be nested. Each cell array in expr indicates a new
			% conditional expression.

			%
			% Error checking
			%
			if (~iscell(expr))
			error(message('images:dicom_create_IOD:conditionalsMustBeCellArrays'))
			end

			if (numel(expr) == 1)
			operator = expr{1};
			operands = {};
			else
			operator = expr{1};
			operands = expr(2:end);
			end


			%
			% Process conditional expressions recursively.
			%
			switch (lower(operator))
			case 'and'

			% This AND short circuits.
			for p = 1:numel(operands)
			tf = check_condition(operands{p}, metadata, dictionary);

			if (~tf)
			return
			end
			end

			case 'equal'

			if (numel(operands) ~= 2)
			error(message('images:dicom_create_IOD:presentOpNeeds2Args'))
			end

			% Get the group and element from the first operand.
			group = operands{1}(2:5);
			element = operands{1}(7:10);

			% Look up the metadata value.
			name = dicom_name_lookup(group, element, dictionary);

			if (isfield(metadata, name))
			tf = isequal(metadata.(name), operands{2});
			else
			tf = false;
			end

			case 'false'

			tf = false;

			case 'not'

			if (numel(operands) ~= 1)
			error(message('images:dicom_create_IOD:tooManyArgsToNot'))
			else
			tf = ~check_condition(operands{1}, metadata, dictionary);
			end

			case 'or'

			% This OR short circuits.
			for p = 1:numel(operands)
			tf = check_condition(operands{p}, metadata, dictionary);

			if (tf)
			return
			end
			end

			case 'present'

			if (numel(operands) ~= 1)
			error(message('images:dicom_create_IOD:presentNeeds1Arg'))
			end

			% Get the group and element from the first operand.
			group = operands{1}(2:5);
			element = operands{1}(7:10);

			% Look up the metadata value.
			name = dicom_name_lookup(group, element, dictionary);

			tf = isfield(metadata, name);

			case 'true'

			tf = true;

			otherwise
			error(message('images:dicom_create_IOD:badConditionalOp', operator))
			end



			function tf = test_module_condition(X, metadata, module, dictionary)
			%TEST_MODULE_CONDITIONS Test whether a module should be created.

			moduleType = module{2};
			moduleCondition = module{3};

			switch (moduleType)
			case {'M'}

			tf = true;

			case {'U', 'C'}

			if (isempty(moduleCondition))
			tf = false;
			return
			end

			switch (moduleCondition{1})
			case 'HAS_FILEMETADATA'
			tf = true;

			case 'HAS_BOLUS'
			tf = isfield(metadata, dicom_name_lookup('0018', '0010', dictionary));

			case 'HAS_OVERLAY'
			tf = isfield(metadata, dicom_name_lookup('6000', '0010', dictionary));

			case 'HAS_VOILUT'
			tf = ((isfield(metadata, dicom_name_lookup('0028', '3010', dictionary))) \|\| ...
			(isfield(metadata, dicom_name_lookup('0028', '1050', dictionary))));

			case 'IS_MULTIFRAME'
			tf = size(X,4) > 1;

			case 'NEEDS_SC_CINE'

			framePointerAttr = dicom_name_lookup('0028', '0009', dictionary);
			if (isfield(metadata, framePointerAttr))

			framePointer = metadata.(framePointerAttr);
			tf = (isequal(framePointer(:), sscanf('0018 1063', '%x')) \|\| ...
			isequal(framePointer(:), sscanf('0018 1065', '%x')));

			else

			tf = false;

			end

			case 'TRUE'
			tf = true;

			case 'FALSE'
			tf = false;

			otherwise
			tf = false;

			end

			otherwise
			error(message('images:dicom_create_IOD:unknownModuleType', moduleType))

			end



			function VR = remap_vr(group, element, vr_map, metadata, dictionary)
			%REMAP_VR Translate a VR mapping into a VR value.

			% Default behavior is an empty VR.
			VR = '';

			if (numel(vr_map) > 1)

			% First item in VR map is an attribute tag.
			cond_name = dicom_name_lookup(vr_map{1}(2:5), vr_map{1}(7:10), dictionary);

			% Remaining items map the conditional attribute's value to a VR.
			if (isfield(metadata, cond_name))

			for p = 2:numel(vr_map)

			if (isequal(metadata.(cond_name), vr_map{p}{1}))

			VR = vr_map{p}{2};
			return

			end

			end

			end

			% If the condition isn't met (or the attribute isn't present),
			% use the first entry in the attribute's VR list.
			if (isempty(VR))

			attr = dicom_dict_lookup(group, element);
			VR = attr.VR;

			end

			elseif (numel(vr_map) == 1)

			% Only item in VR map is the VR.
			VR = vr_map{1};

			end



			function VR = get_pixelData_VR(metadata, dictionary)
			%GET_PIXELDATA_VR Get the VR for (7FE0,0010).
			%
			% PS 3.5 Sec. 8 contains the rules for picking the VR of (7FE0,0010).

			uid_details = dicom_uid_decode(metadata.(dicom_name_lookup('0002','0010', dictionary)));

			if (uid_details.Compressed)

			% Compressed pixels are always stored OB.
			VR = 'OB';

			elseif (isequal(uid_details.VR, 'Implicit'))

			% Implicit VR transfer syntaxes are always OW.
			VR = 'OW';

			else

			% The VR of other Explicit VR transfter syntaxes depends on the bit
			% depth of the pixels.
			if (metadata.(dicom_name_lookup('0028','0100', dictionary)) <= 8)
			VR = 'OB';
			else
			VR = 'OW';
			end

			end



			function tf = check_values(attr_data, enum_values)
			%CHECK_VALUE Verify an attribute's data against required values.

			if (iscellstr(enum_values))

			tf = any(strcmp(attr_data, enum_values));

			else

			tmp = 0;
			for p = 1:numel(enum_values)
			tmp = tmp + isequal(attr_data, enum_values{p});
			end

			tf = (tmp > 0);

			end



			function [private_attrs, private_names] = find_private_metadata(metadata, dictionary)
			%FIND_PRIVATE_METADATA Find user-contributed private metadata.

			% Private attributes have odd element values in their tag.
			%
			% By default these values come back from DICOMINFO with field names like
			% "Private_0029_1004" and "Private_0029_10xx_Creator".
			%
			% If a custom data dictionary is provided, the names will look like
			% regular DICOM attributes.

			private_attrs = {};
			private_names = {};

			% Look through metadata for private attributes.
			fields = fieldnames(metadata);

			for p = 1:numel(fields)

			tag = dicom_tag_lookup(fields{p}, dictionary);

			if ((~isempty(tag)) && (rem(tag(1), 2) ~= 0))
			private_attrs{end + 1} = tag; %#ok<AGROW>
			private_names{end + 1} = fields{p}; %#ok<AGROW>
			end

			end



			function [attrs, status] = process_private(private_tags, private_names, metadata, dictionary, txfr)
			%PROCESS_PRIVATE Create attributes from the private data list.

			% Delay judgement of Private attributes.
			status.BadAttribute = {};
			status.MissingCondition = {};
			status.MissingData = {};
			status.SuspectAttribute = {};

			% Preallocate the attribute structure.
			attrs(numel(private_tags)).Group = [];
			attrs(end).Element = [];
			attrs(end).VR = '';
			attrs(end).VM = [];
			attrs(end).Data = [];

			% Convert the private data to a writable attribute.
			for p = 1:numel(private_tags)
			attrs(p) = dicom_convert_meta_to_attr(private_names{p}, metadata, dictionary, txfr);
			end

R8218/CT/private/dicom_create_IOD.m4b5355e68791master

dicom_create_IOD.mView Options

R8218/CT/private/dicom_create_IOD.m
4b5355e68791master

dicom_create_IOD.m
View Options