public async Task TestImageOpening()
{
var prostateFolder = TestData.GetFullImagesPath("sample_dicom");
var files = Directory.EnumerateFiles(prostateFolder)
.Where(x => !x.EndsWith("rtstruct_mod.dcm")).ToList();
var dicomFiles = new DicomFile[files.Count];
for (var i = 0; i < files.Count; i++)
{
dicomFiles[i] = await DicomFile.OpenAsync(files[i]);
}
var volume = DicomSeriesReader.BuildVolume(
dicomFiles.Select(x => x.Dataset),
new ModerateGeometricAcceptanceTest(string.Empty, string.Empty), true);
Assert.AreEqual(512, volume.DimX);
Assert.AreEqual(262144, volume.DimXY);
Assert.AreEqual(512, volume.DimY);
Assert.AreEqual(2, volume.DimZ);
Assert.AreEqual(3, volume.Dimensions);
Assert.AreEqual(524288, volume.Length);
Assert.AreEqual(-250d, volume.Origin.X);
Assert.AreEqual(-250d, volume.Origin.Y);
Assert.AreEqual(125.5d, volume.Origin.Z);
Assert.AreEqual(0.9765625, volume.SpacingX);
Assert.AreEqual(0.9765625, volume.SpacingY);
Assert.AreEqual(2.86102294921875, volume.VoxelVolume);
Assert.AreEqual(-1000, volume.Array[3453]);
Assert.AreEqual(-1000, volume.Array[8453]);
Assert.AreEqual(-1000, volume.Array[10453]);
Assert.AreEqual(-1000, volume.Array[100453]);
}
/// <summary>
/// Attempt to load a volume from the given SeriesUID for the given DicomFolderContents
/// </summary>
/// <param name="dfc">A pre-built description of DICOM contents within a particular folder</param>
/// <param name="seriesUID">The DICOM seriesUID you wish to construct a volume for</param>
/// <param name="acceptanceTests">An implementation of IVolumeGeometricAcceptanceTest defining the geometric constraints of your application</param>
/// <param name="loadStructuresIfExists">True if rt-structures identified in the folder and referencing seriesUID should be loaded</param>
/// <param name="supportLossyCodecs">If you wish to accept lossy encodings of image pixel data</param>
/// <returns></returns>
private static VolumeLoaderResult LoadDicomSeries(
DicomFolderContents dfc, DicomUID seriesUID, IVolumeGeometricAcceptanceTest acceptanceTests, bool loadStructuresIfExists, bool supportLossyCodecs)
{
try
{
var dicomSeriesContent = dfc.Series.FirstOrDefault((s) => s.SeriesUID == seriesUID);
var warnings = new List <string>();
RadiotherapyStruct rtStruct = null;
if (dicomSeriesContent != null)
{
var volumeData = DicomSeriesReader.BuildVolume(dicomSeriesContent.Content.Select(x => x.File.Dataset), acceptanceTests, supportLossyCodecs);
if (volumeData != null && loadStructuresIfExists)
{
var rtStructData = dfc.RTStructs.FirstOrDefault(rt => rt.SeriesUID == seriesUID);
if (rtStructData != null)
{
if (rtStructData.Content.Count == 1)
{
var rtStructAndWarnings = RtStructReader.LoadContours(
rtStructData.Content.First().File.Dataset,
volumeData.Transform.DicomToData,
seriesUID.UID,
null,
false);
rtStruct = rtStructAndWarnings.Item1;
var warning = rtStructAndWarnings.Item2;
if (!string.IsNullOrEmpty(warning))
{
warnings.Add(warning);
}
}
else if (rtStructData.Content.Count > 1)
{
warnings.Add("There is more than 1 RT STRUCT referencing this series - skipping structure set load");
}
}
}
var dicomIdentifiers = dicomSeriesContent.Content.Select((v) => DicomIdentifiers.ReadDicomIdentifiers(v.File.Dataset)).ToArray();
if (rtStruct == null)
{
rtStruct = RadiotherapyStruct.CreateDefault(dicomIdentifiers);
}
var result = new MedicalVolume(
volumeData,
dicomIdentifiers,
dicomSeriesContent.Content.Select((d) => d.Path).ToArray(),
rtStruct);
return(new VolumeLoaderResult(seriesUID.UID, result, null, warnings));
}
throw new Exception("Could not find that series");
}
catch (Exception oops)
{
return(new VolumeLoaderResult(seriesUID.UID, null, oops, new List <string>()));
}
}
