예제 #1
0
        private void LoadStructureExecute()
        {
            Debug.WriteLine("LoadStructure pressed!");
            InitialState = false;

            StructureSet.Clear();
            string sFile = @"DicomSamples\RSAnonymized.dcm";

            DicomLib.RtStructureDicom rtStructureDicom = RtStructureDicom.ImportFromFile(sFile);
            var ss = rtStructureDicom.StructureSet;

            foreach (Structure item in ss.StructureCollection)
            {
                StructureSet.Add(item);
                Trace.WriteLine(item);
            }

            var structure = ss.StructureCollection.FirstOrDefault(s => s.Name == "Prostate" || s.Name == "Prostata");

            EventAggregator.GetEvent <StrucutureSetUpdatedEvent>().Publish(new StrucutureSetUpdatedEventArg(StructureSet));
        }
예제 #2
0
        private static RtStructureDicom ReadFile(string sDicomFile)
        {
            DicomFile df      = DicomFile.Open(sDicomFile);
            var       dataset = df.Dataset;

            var modality = dataset.Get <string>(DicomTag.Modality);

            if (!modality.ToUpper().Contains("RTSTR"))
            {
                throw new ArgumentException("No RT Structure file provided.");
            }

            string ssName = String.Empty;

            try
            {
                ssName = dataset.Get <string>(DicomTag.StructureSetName);
            }
            catch (Exception e)
            {
                Trace.WriteLine(e);
            }

            var ssObject = new StructureSet(ssName);

            RtStructureDicom rtStructureDicom = new RtStructureDicom(ssObject);

            var ssFrameOfReferenceSequence = dataset.Get <DicomSequence>(DicomTag.ReferencedFrameOfReferenceSequence);

            foreach (DicomDataset item in ssFrameOfReferenceSequence.Items)
            {
                try
                {
                    var id = item.Get <string>(DicomTag.FrameOfReferenceUID);
                    ssObject.AddFrameOfReferenceUid(id);
                }
                catch (Exception e)
                {
                    Trace.WriteLine(e);
                }
            }

            var structureSetRoiSequence = dataset.Get <DicomSequence>(DicomTag.StructureSetROISequence);

            foreach (DicomDataset item in structureSetRoiSequence.Items)
            {
                var name  = item.Get <string>(DicomTag.ROIName);
                var roiId = item.Get <int>(DicomTag.ROINumber);
                //Trace.WriteLine(item.Get<string>(DicomTag.ROIName));
                Structure structure = new Structure(name, roiId);
                ssObject.Add(structure);
            }

            var roiContourSequence = dataset.Get <DicomSequence>(DicomTag.ROIContourSequence);

            foreach (DicomDataset item in roiContourSequence.Items)
            {
                try
                {
                    var       roiColor        = item.Get <int[]>(DicomTag.ROIDisplayColor);
                    var       contourSequence = item.Get <DicomSequence>(DicomTag.ContourSequence);
                    var       roiId           = item.Get <int>(DicomTag.ReferencedROINumber);
                    Structure structure       = ssObject.GetStructureById(roiId);
                    structure.Color = Color.FromRgb((byte)roiColor[0], (byte)roiColor[1], (byte)roiColor[2]);
                    foreach (DicomDataset contourItem in contourSequence.Items)
                    {
                        Contour contourOnAPlane      = new Contour();
                        int     numerOfContourPoints = contourItem.Get <int>(DicomTag.NumberOfContourPoints);
                        var     contourData          = contourItem.Get <float[]>(DicomTag.ContourData);
                        Debug.Assert(contourData.Length == numerOfContourPoints * 3);

                        contourOnAPlane.SetZPosition(contourData[2]);
                        for (int i = 0; i < contourData.Length; i = i + 3)
                        {
                            float x = contourData[i];
                            float y = contourData[i + 1];
                            float z = contourData[i + 2];

                            contourOnAPlane.AddPoint(x, y, z);
                        }
                        structure.AddContour(contourOnAPlane);
                    }
                }
                catch (Exception e)
                {
                    Trace.WriteLine(e);
                }
            }
            return(rtStructureDicom);
        }
예제 #3
0
 public static RtStructureDicom ImportFromFile(string sDicomFile)
 {
     return(RtStructureDicom.ReadFile(sDicomFile));
 }