private static void Save(XElement parent, DcmDataset dataset)
{
foreach (XElement attr in parent.Elements("attr"))
{
DicomTag tag = DicomTag.Parse(attr.Attribute("tag").Value);
DicomVR vr = DicomVR.Lookup(attr.Attribute("vr").Value);
int len = int.Parse(attr.Attribute("len").Value, CultureInfo.InvariantCulture);
if (vr == DicomVR.SQ)
{
DcmItemSequence seq = new DcmItemSequence(tag);
foreach (XElement itm in attr.Elements("item"))
{
DcmItemSequenceItem item = new DcmItemSequenceItem();
Save(itm, item.Dataset);
seq.AddSequenceItem(item);
}
dataset.AddItem(seq);
}
else if (len == -1)
{
DcmFragmentSequence seq = new DcmFragmentSequence(tag, vr);
bool first = true;
foreach (XElement itm in attr.Elements("item"))
{
if (first)
{
SaveFragmentOffsetTable(itm, seq);
first = false;
}
else
{
SaveFragmentItem(itm, seq);
}
}
dataset.AddItem(seq);
}
else
{
DcmElement element = DcmElement.Create(tag, vr);
element.SetValueString(attr.FirstText());
dataset.AddItem(element);
}
}
}
/// <summary>
/// Read dataset from stream
/// </summary>
/// <param name="stopAtTag">End parsing at this tag</param>
/// <param name="options">DICOM read options</param>
/// <returns>Status code</returns>
public DicomReadStatus Read(DicomTag stopAtTag, DicomReadOptions options)
{
// Counters:
// _remain - bytes remaining in stream
// _bytes - estimates bytes to end of dataset
// _read - number of bytes read from stream
try {
_need = 0;
_remain = _stream.Length - _stream.Position;
while (_remain > 0)
{
DicomReadStatus status = ParseTag(stopAtTag, options);
if (status == DicomReadStatus.SuccessEndRead)
{
return(DicomReadStatus.Success);
}
if (status != DicomReadStatus.Success)
{
return(status);
}
status = ParseVR(options);
if (status != DicomReadStatus.Success)
{
return(status);
}
status = ParseLength(options);
if (status != DicomReadStatus.Success)
{
return(status);
}
if (_tag.IsPrivate)
{
if (_tag.Element != 0x0000 && _tag.Element <= 0x00ff)
{
// handle UN private creator id
if (_vr != DicomVR.LO && Flags.IsSet(options, DicomReadOptions.ForcePrivateCreatorToLO))
{
Dicom.Debug.Log.Warn("Converting Private Creator VR from '{0}' to 'LO'", _vr.VR);
_vr = DicomVR.LO;
}
}
}
if (_vr == DicomVR.UN && _syntax.IsExplicitVR && Flags.IsSet(options, DicomReadOptions.UseDictionaryForExplicitUN))
{
_vr = _tag.Entry.DefaultVR;
}
if (_fragment != null)
{
status = InsertFragmentItem(options);
if (status != DicomReadStatus.Success)
{
return(status);
}
}
else if (_sqs.Count > 0 &&
(_tag == DicomTags.Item ||
_tag == DicomTags.ItemDelimitationItem ||
_tag == DicomTags.SequenceDelimitationItem))
{
status = InsertSequenceItem(options);
if (status != DicomReadStatus.Success)
{
return(status);
}
}
else
{
if (_sqs.Count > 0)
{
DcmItemSequence sq = _sqs.Peek();
if (sq.StreamLength != UndefinedLength)
{
long end = sq.StreamPosition + 8 + sq.StreamLength;
if (_syntax.IsExplicitVR)
{
end += 2 + 2;
}
if ((_stream.Position - _offset) >= end)
{
if (_sds.Count == _sqs.Count)
{
_sds.Pop();
}
_sqs.Pop();
}
}
}
if (_len == UndefinedLength)
{
if (_vr == DicomVR.SQ)
{
DcmItemSequence sq = new DcmItemSequence(_tag, _pos, _len, _endian);
InsertDatasetItem(sq, options);
_sqs.Push(sq);
}
else
{
_fragment = new DcmFragmentSequence(_tag, _vr, _pos, _endian);
InsertDatasetItem(_fragment, options);
}
}
else
{
if (_vr == DicomVR.SQ)
{
DcmItemSequence sq = new DcmItemSequence(_tag, _pos, _len, _endian);
InsertDatasetItem(sq, options);
_sqs.Push(sq);
}
else
{
if (_len > _remain)
{
return(NeedMoreData(_len));
}
DcmElement elem = DcmElement.Create(_tag, _vr, _pos, _endian, CurrentBuffer(options));
_remain -= _len;
_read += _len;
InsertDatasetItem(elem, options);
}
}
}
_tag = null;
_vr = null;
_len = UndefinedLength;
}
return(DicomReadStatus.Success);
}
catch (EndOfStreamException) {
// should never happen
return(DicomReadStatus.UnknownError);
}
}
public string MakeGreyDicom(byte[] greybytes, ushort imgwidth, ushort imgheight)
{
DcmUID studyUid = DcmUID.Generate();
DcmUID seriesUid = DcmUID.Generate(studyUid, 1);
DcmUID instUid = DcmUID.Generate(seriesUid, 1);
DcmDataset data = new DcmDataset(DcmTS.ExplicitVRBigEndian); //.ImplicitVRLittleEndian ok
data.AddElementWithValue(DcmTags.SOPClassUID, DcmUIDs.CTImageStorage); //ComputedRadiographyImageStorage ok
//data.AddElementWithValue(DcmTags.SOPClassUID, DcmUIDs .SecondaryCapture);
data.AddElementWithValue(DcmTags.StudyInstanceUID, studyUid);
data.AddElementWithValue(DcmTags.SeriesInstanceUID, seriesUid);
data.AddElementWithValue(DcmTags.SOPInstanceUID, instUid);//"1.3.6.1.4.1.30071.6.635719267134010719.1.1"
//data.AddElementWithValue(DcmTags.MediaStorageSOPClassUID, DcmUIDs.ImplicitVRLittleEndian);
//data.AddElementWithValueString(DcmTags.MediaStorageSOPClassUID, DcmUIDs.ComputedRadiographyImageStorage.ToString());
//type 2 attributes
////data.AddElement(DcmTags.PrinterStatus);
if (tags.ContainsKey("0010,0020"))
{
data.AddElementWithValueString(DcmTags.PatientID, tags["0010,0020"].Substring(5));
}
if (tags.ContainsKey("0010,0010"))
{
data.AddElementWithValueString(DcmTags.PatientsName, tags["0010,0010"].Substring(5));
}
if (tags.ContainsKey("0010,0030"))
{
data.AddElementWithValueString(DcmTags.PatientsBirthDate, tags["0010,0030"].Substring(5));
}
if (tags.ContainsKey("0010,0040"))
{
data.AddElementWithValueString(DcmTags.PatientsSex, tags["0010,0040"].Substring(5));
}
if (tags.ContainsKey("0010,1010"))
{
data.AddElementWithValueString(DcmTags.PatientsAge, tags["0010,1010"].Substring(5));
}
if (tags.ContainsKey("0008,0005"))
{
data.AddElementWithValueString(DcmTags.SpecificCharacterSet, tags["0008,0005"].Substring(5));
}
if (tags.ContainsKey("0008,0008"))
{
data.AddElementWithValueString(DcmTags.ImageType, tags["0008,0008"].Substring(5));
}
//if (tags.ContainsKey("0008,0016"))
// data.AddElementWithValueString(DcmTags.ContentTime, DateTime.Now.ToString());
//if (tags.ContainsKey("0008,0018"))
// data.AddElementWithValueString(DcmTags.ContentTime, DateTime.Now.ToString());
if (tags.ContainsKey("0008,0020"))
{
data.AddElementWithValueString(DcmTags.StudyDate, tags["0008,0020"].Substring(5));
}
if (tags.ContainsKey("0008,0021"))
{
data.AddElementWithValueString(DcmTags.SeriesDate, tags["0008,0021"].Substring(5));
}
if (tags.ContainsKey("0008,0022"))
{
data.AddElementWithValueString(DcmTags.AcquisitionDate, tags["0008,0022"].Substring(5));
}
if (tags.ContainsKey("0008,0023"))
{
data.AddElementWithValueString(DcmTags.ContentDate, tags["0008,0023"].Substring(5));
}
if (tags.ContainsKey("0008,002a"))
{
data.AddElementWithValueString(DcmTags.AcquisitionDateTime, tags["0008,002a"].Substring(5));
}
if (tags.ContainsKey("0008,0030"))
{
data.AddElementWithValueString(DcmTags.StudyTime, tags["0008,0030"].Substring(5));
}
if (tags.ContainsKey("0008,0031"))
{
data.AddElementWithValueString(DcmTags.SeriesTime, tags["0008,0031"].Substring(5));
}
if (tags.ContainsKey("0008,0032"))
{
data.AddElementWithValueString(DcmTags.AcquisitionTime, tags["0008,0032"].Substring(5));
}
if (tags.ContainsKey("0008,0033"))
{
data.AddElementWithValueString(DcmTags.ContentTime, tags["0008,0033"].Substring(5));
}
if (tags.ContainsKey("0008,0050"))
{
data.AddElementWithValueString(DcmTags.AcquisitionNumber, tags["0008,0050"].Substring(5));
}
if (tags.ContainsKey("0008,0060"))
{
data.AddElementWithValueString(DcmTags.Modality, tags["0008,0060"].Substring(5));
}
if (tags.ContainsKey("0008,0070"))
{
data.AddElementWithValueString(DcmTags.Manufacturer, tags["0008,0070"].Substring(5));
}
if (tags.ContainsKey("0008,0080"))
{
data.AddElementWithValueString(DcmTags.InstitutionName, tags["0008,0080"].Substring(5));
}
if (tags.ContainsKey("0008,0081"))
{
data.AddElementWithValueString(DcmTags.InstitutionAddress, tags["0008,0081"].Substring(5));
}
if (tags.ContainsKey("0008,0090"))
{
data.AddElementWithValueString(DcmTags.ReferringPhysiciansName, tags["0008,0090"].Substring(5));
}
if (tags.ContainsKey("0008,1010"))
{
data.AddElementWithValueString(DcmTags.StationName, tags["0008,1010"].Substring(5));
}
if (tags.ContainsKey("0008,1030"))
{
data.AddElementWithValueString(DcmTags.StudyDescription, tags["0008,1030"].Substring(5));
}
if (tags.ContainsKey("0008,103e"))
{
data.AddElementWithValueString(DcmTags.SeriesDescription, tags["0008,103e"].Substring(5));
}
if (tags.ContainsKey("0008,1090"))
{
data.AddElementWithValueString(DcmTags.ManufacturersModelName, tags["0008,1090"].Substring(5));
}
if (tags.ContainsKey("0018,0010"))
{
data.AddElementWithValueString(DcmTags.ContrastBolusAgent, tags["0018,0010"].Substring(5));
}
if (tags.ContainsKey("0018,0015"))
{
data.AddElementWithValueString(DcmTags.BodyPartExamined, tags["0018,0015"].Substring(5));
}
if (tags.ContainsKey("0018,0050"))
{
data.AddElementWithValueString(DcmTags.SliceThickness, tags["0018,0050"].Substring(5));
}
if (tags.ContainsKey("0018,0060"))
{
data.AddElementWithValueString(DcmTags.KVP, tags["0018,0060"].Substring(5));
}
if (tags.ContainsKey("0018,0090"))
{
data.AddElementWithValueString(DcmTags.DataCollectionDiameter, tags["0018,0090"].Substring(5));
}
if (tags.ContainsKey("0018,1000"))
{
data.AddElementWithValueString(DcmTags.DeviceSerialNumber, tags["0018,1000"].Substring(5));
}
if (tags.ContainsKey("0018,1020"))
{
data.AddElementWithValueString(DcmTags.SoftwareVersions, tags["0018,1020"].Substring(5));
}
if (tags.ContainsKey("0018,1030"))
{
data.AddElementWithValueString(DcmTags.ProtocolName, tags["0018,1030"].Substring(5));
}
if (tags.ContainsKey("0018,1041"))
{
data.AddElementWithValueString(DcmTags.ContrastBolusVolume, tags["0018,1041"].Substring(5));
}
if (tags.ContainsKey("0018,1042"))
{
data.AddElementWithValueString(DcmTags.ContrastBolusStartTime, tags["0018,1042"].Substring(5));
}
if (tags.ContainsKey("0018,1043"))
{
data.AddElementWithValueString(DcmTags.ContrastBolusStopTime, tags["0018,1043"].Substring(5));
}
if (tags.ContainsKey("0018,1044"))
{
data.AddElementWithValueString(DcmTags.ContrastBolusTotalDose, tags["0018,1044"].Substring(5));
}
if (tags.ContainsKey("0018,1046"))
{
data.AddElementWithValueString(DcmTags.ContrastFlowRate, tags["0018,1046"].Substring(5));
}
if (tags.ContainsKey("0018,1047"))
{
data.AddElementWithValueString(DcmTags.ContrastFlowDuration, tags["0018,1047"].Substring(5));
}
if (tags.ContainsKey("0018,1049"))
{
data.AddElementWithValueString(DcmTags.ContrastBolusIngredientConcentration, tags["0018,1049"].Substring(5));
}
if (tags.ContainsKey("0018,1100"))
{
data.AddElementWithValueString(DcmTags.ReconstructionDiameter, tags["0018,1100"].Substring(5));
}
if (tags.ContainsKey("0018,1110"))
{
data.AddElementWithValueString(DcmTags.DistanceSourceToDetector, tags["0018,1110"].Substring(5));
}
if (tags.ContainsKey("0018,1111"))
{
data.AddElementWithValueString(DcmTags.DistanceSourceToPatient, tags["0018,1111"].Substring(5));
}
if (tags.ContainsKey("0018,1120"))
{
data.AddElementWithValueString(DcmTags.GantryDetectorTilt, tags["0018,1120"].Substring(5));
}
if (tags.ContainsKey("0018,1130"))
{
data.AddElementWithValueString(DcmTags.TableHeight, tags["0018,1130"].Substring(5));
}
if (tags.ContainsKey("0018,1140"))
{
data.AddElementWithValueString(DcmTags.RotationDirection, tags["0018,1140"].Substring(5));
}
if (tags.ContainsKey("0018,1150"))
{
data.AddElementWithValueString(DcmTags.ExposureTime, tags["0018,1150"].Substring(5));
}
if (tags.ContainsKey("0018,1151"))
{
data.AddElementWithValueString(DcmTags.XRayTubeCurrent, tags["0018,1151"].Substring(5));
}
if (tags.ContainsKey("0018,1152"))
{
data.AddElementWithValueString(DcmTags.Exposure, tags["0018,1152"].Substring(5));
}
if (tags.ContainsKey("0018,1160"))
{
data.AddElementWithValueString(DcmTags.FilterType, tags["0018,1160"].Substring(5));
}
if (tags.ContainsKey("0018,1170"))
{
data.AddElementWithValueString(DcmTags.GeneratorPower, tags["0018,1170"].Substring(5));
}
if (tags.ContainsKey("0018,1190"))
{
data.AddElementWithValueString(DcmTags.FocalSpots, tags["0018,1190"].Substring(5));
}
if (tags.ContainsKey("0018,1200"))
{
data.AddElementWithValueString(DcmTags.DateOfLastCalibration, tags["0018,1200"].Substring(5));
}
if (tags.ContainsKey("0018,1201"))
{
data.AddElementWithValueString(DcmTags.TimeOfLastCalibration, tags["0018,1201"].Substring(5));
}
if (tags.ContainsKey("0018,1210"))
{
data.AddElementWithValueString(DcmTags.ConvolutionKernel, tags["0018,1210"].Substring(5));
}
if (tags.ContainsKey("0018,5100"))
{
data.AddElementWithValueString(DcmTags.PatientPosition, tags["0018,5100"].Substring(5));
}
//if (tags.ContainsKey("0020,000D"))
// data.AddElementWithValueString(DcmTags.ContrastBolusStopTime, DateTime.Now.ToString());
//if (tags.ContainsKey("0020,000E"))
// data.AddElementWithValueString(DcmTags.ContrastBolusStopTime, DateTime.Now.ToString());
if (tags.ContainsKey("0020,0010"))
{
data.AddElementWithValueString(DcmTags.StudyID, tags["0020,0010"].Substring(5));
}
if (tags.ContainsKey("0020,0011"))
{
data.AddElementWithValueString(DcmTags.SeriesNumber, tags["0020,0011"].Substring(5));
}
if (tags.ContainsKey("0020,0012"))
{
data.AddElementWithValueString(DcmTags.AccessionNumber, tags["0020,0012"].Substring(5));
}
if (tags.ContainsKey("0020,0013"))
{
data.AddElementWithValueString(DcmTags.InstanceNumber, tags["0020,0013"].Substring(5));
}
if (tags.ContainsKey("0020,0032"))
{
data.AddElementWithValueString(DcmTags.ImagePositionPatient, tags["0020,0032"].Substring(5));
}
if (tags.ContainsKey("0020,0037"))
{
data.AddElementWithValueString(DcmTags.ImageOrientationPatient, tags["0020,0037"].Substring(5));
}
if (tags.ContainsKey("0020,0052"))
{
data.AddElementWithValueString(DcmTags.FrameOfReferenceUID, tags["0020,0052"].Substring(5));
}
if (tags.ContainsKey("0020,1040"))
{
data.AddElementWithValueString(DcmTags.PositionReferenceIndicator, tags["0020,1040"].Substring(5));
}
if (tags.ContainsKey("0020,1041"))
{
data.AddElementWithValueString(DcmTags.SliceLocation, tags["0020,1041"].Substring(5));
}
if (tags.ContainsKey("0020,4000"))
{
data.AddElementWithValueString(DcmTags.ImageComments, tags["0020,4000"].Substring(5));
}
//data.AddElementWithValueString(DcmTags.StudyTime, DateTime.Now.ToString());
//data.AddElementWithValueString(DcmTags.AccessionNumber, "");
//data.AddElementWithValueString(DcmTags.ReferringPhysiciansName, "");
//data.AddElementWithValueString(DcmTags.StudyID, "1");
//data.AddElementWithValueString(DcmTags.SeriesNumber, "1");
//data.AddElementWithValueString(DcmTags.ModalitiesInStudy, "CT");//CR
//data.AddElementWithValueString(DcmTags.Modality, "CT");//CR
//data.AddElementWithValueString(DcmTags.NumberOfStudyRelatedInstances, "1");
//data.AddElementWithValueString(DcmTags.NumberOfStudyRelatedSeries, "1");
//data.AddElementWithValueString(DcmTags.NumberOfSeriesRelatedInstances, "1");
//data.AddElementWithValueString(DcmTags.PatientOrientation, "HFS");//F/A
//data.AddElementWithValueString(DcmTags.ImageLaterality, "U");
if (tags.ContainsKey("0028,1050"))
{
data.AddElementWithValueString(DcmTags.WindowCenter, "1113");
}
if (tags.ContainsKey("0028,1051"))
{
data.AddElementWithValueString(DcmTags.WindowWidth, "749");
}
//data.AddElementWithValueString(DcmTags.WindowCenterWidthExplanation, "WINDOW1\\WINDOW2");
data.AddElementWithValueString(DcmTags.PixelRepresentation, "0");
data.AddElementWithValueString(DcmTags.RescaleIntercept, "0");//0
data.AddElementWithValueString(DcmTags.RescaleSlope, "1");
//data.AddElementWithValueString(DcmTags.RotationDirection, "CW");
//ushort bitdepth = 2;未使用过
DcmPixelData pixelData = new DcmPixelData(DcmTS.ImplicitVRLittleEndian);
pixelData.PixelRepresentation = 0;//ok
pixelData.ImageWidth = imgwidth;
pixelData.ImageHeight = imgheight;
pixelData.SamplesPerPixel = 1; //ok
pixelData.HighBit = 15; //ok
pixelData.BitsStored = 16; //ok
pixelData.BitsAllocated = 16; //ok
//pixelData.SamplesPerPixel = 1;
//pixelData.HighBit = 7;
//pixelData.BitsStored = 8;
//pixelData.BitsAllocated = 8;
pixelData.ImageType = "ORIGINAL\\PRIMARY\\AXIAL";
pixelData.PhotometricInterpretation = "MONOCHROME2";//2 byte gray? //ok
//pixelData.FragmentSize
//pixelData.IsLossy = true;
//pixelData.LossyCompressionMethod = "01";
pixelData.PixelDataElement = DcmElement.Create(DcmTags.PixelData, DcmVR.OW); //OB: Other Byte, OW: Other Word
//pixelData.AddFrame(bmpBytes);
pixelData.AddFrame(greybytes);
pixelData.UpdateDataset(data);
DicomFileFormat ff = new DicomFileFormat(data);
//string fileout = System.IO.Path.Combine(Directory.GetCurrentDirectory(), "greyimg_test.dcm");
ff.Save(fileout, Dicom.DicomWriteOptions.Default);//Default
ff = null;
return(fileout);
}