public void LoadAndSaveMedicalVolumeTest()
{
var directory = TestData.GetFullImagesPath("sample_dicom");
var acceptanceTests = new StrictGeometricAcceptanceTest(string.Empty, string.Empty);
var medicalVolume = MedIO.LoadAllDicomSeriesInFolderAsync(directory, acceptanceTests).Result.First().Volume;
Directory.CreateDirectory(_tempFolder);
Console.WriteLine($"Directory created {_tempFolder}");
Volume3D <byte>[] contourVolumes = new Volume3D <byte> [medicalVolume.Struct.Contours.Count];
for (int i = 0; i < medicalVolume.Struct.Contours.Count; i++)
{
contourVolumes[i] = new Volume3D <byte>(
medicalVolume.Volume.DimX,
medicalVolume.Volume.DimY,
medicalVolume.Volume.DimZ,
medicalVolume.Volume.SpacingX,
medicalVolume.Volume.SpacingY,
medicalVolume.Volume.SpacingZ,
medicalVolume.Volume.Origin,
medicalVolume.Volume.Direction);
contourVolumes[i].Fill(medicalVolume.Struct.Contours[i].Contours, (byte)1);
}
// Calculate contours based on masks
var rtContours = new List <RadiotherapyContour>();
for (int i = 0; i < medicalVolume.Struct.Contours.Count; i++)
{
var contour = medicalVolume.Struct.Contours[i];
var contourForAllSlices = GetContoursForAllSlices(contourVolumes[i]);
var rtcontour = new DicomRTContour(
contour.DicomRtContour.ReferencedRoiNumber,
contour.DicomRtContour.RGBColor,
contourForAllSlices);
DicomRTStructureSetROI rtROIstructure = new DicomRTStructureSetROI(
contour.StructureSetRoi.RoiNumber,
contour.StructureSetRoi.RoiName,
string.Empty,
ERoiGenerationAlgorithm.Semiautomatic);
DicomRTObservation observation = new DicomRTObservation(
contour.DicomRtObservation.ReferencedRoiNumber, new DicomPersonNameConverter("Left^Richard^^Dr"), ROIInterpretedType.EXTERNAL);
rtContours.Add(new RadiotherapyContour(rtcontour, rtROIstructure, observation));
}
var rtStructureSet = new RadiotherapyStruct(
medicalVolume.Struct.StructureSet,
medicalVolume.Struct.Patient,
medicalVolume.Struct.Equipment,
medicalVolume.Struct.Study,
medicalVolume.Struct.RTSeries,
rtContours);
MedIO.SaveMedicalVolumeAsync(_tempFolder, new MedicalVolume(
medicalVolume.Volume,
medicalVolume.Identifiers,
medicalVolume.FilePaths,
rtStructureSet)).Wait();
var medicalVolume2 = MedIO.LoadAllDicomSeriesInFolderAsync(_tempFolder, acceptanceTests).Result.First().Volume;
foreach (var radiotherapyContour in medicalVolume2.Struct.Contours.Where(x => x.DicomRtContour.DicomRtContourItems.First().GeometricType == "CLOSED_PLANAR"))
{
var savedContour =
medicalVolume2.Struct.Contours.First(x => x.StructureSetRoi.RoiName == radiotherapyContour.StructureSetRoi.RoiName);
foreach (var contour in radiotherapyContour.Contours)
{
Assert.AreEqual(radiotherapyContour.DicomRtObservation.ROIInterpretedType, ROIInterpretedType.EXTERNAL);
for (int i = 0; i < contour.Value.Count; i++)
{
if (!contour.Value[i].Equals(savedContour.Contours.ContoursForSlice(contour.Key)[i]))
{
Console.WriteLine(radiotherapyContour.StructureSetRoi.RoiName);
Assert.Fail();
}
}
}
}
}
