public byte[] DICOMSeriesToVolume(DicomDirectoryRecord series, out Int3 size)
{
size = GetSizeForRecord(series);
var voxels = new VolumeBuffer <Color32>(size);
var tex = new Texture2D(2, 2);
var instanceNumbers = series.LowerLevelDirectoryRecordCollection.Select(x => x.Get <int>(DicomTag.InstanceNumber)).OrderBy(x => x).ToArray();
for (var z = 0; z < instanceNumbers.Length; z++)
{
var frame = instanceNumbers[z];
tex = loadDicomInstance.GetTexture2DForRecord(series, frame);
var fromPixels = tex.GetPixels32();
for (var y = 0; y < size.y; y++)
{
for (var x = 0; x < size.x; x++)
{
var from = fromPixels[x * resolution + ((size.y - 1 - y) * size.x) * resolution * 2];
voxels.SetVoxel(new Int3(x, y, z), from);
}
}
}
voxels.ClearEdges(new Color32(0, 0, 0, 0));
return(VolumeTextureUtils.Color32ArrayToByteArray(voxels.DataArray));
}
public void ButtonPush(string button)
{
if (button == "3D")
{
transform.Find("zstack slider").gameObject.SetActive(false);
var vc = transform.Find("Volume Cube");
vc.gameObject.SetActive(true);
if (vc.GetComponent <Renderer>().material.GetTexture("_Volume") == null)
{
var id = record.Get <string>(DicomTag.SeriesInstanceUID);
var path = Path.Combine(Application.persistentDataPath, "Volumes", id);
byte[] vol;
Int3 size;
var startTime = Time.realtimeSinceStartup;
if (File.Exists(path) && useCache) // a cached volume for this series exists
{
vol = File.ReadAllBytes(path);
size = GetSizeForRecord(record);
}
else
{
vol = DICOMSeriesToVolume(record, out size);
File.WriteAllBytes(path, vol);
}
var volumeSizePow2 = MathExtensions.PowerOfTwoGreaterThanOrEqualTo(size);
var tex3D = VolumeTextureUtils.BuildTexture(vol, size, volumeSizePow2);
Debug.Log("created volume in " + (Time.realtimeSinceStartup - startTime) + "s");
vcRenderer.material.SetTexture("_Volume", tex3D);
var first = loadDicomInstance.GetImageForRecord(record, 0);
var last = loadDicomInstance.GetImageForRecord(record, -2);
var sliceLocFirst = first.Dataset.Get <float>(DicomTag.SliceLocation);
var sliceLocLast = last.Dataset.Get <float>(DicomTag.SliceLocation);
var pos = first.Dataset.Get <string[]>(DicomTag.ImagePositionPatient);
var spacing = first.Dataset.Get <float[]>(DicomTag.PixelSpacing);
var zDepth = sliceLocLast - sliceLocFirst;
var xWidth = spacing[0] * first.Dataset.Get <int>(DicomTag.Rows);
var yHeight = spacing[1] * first.Dataset.Get <int>(DicomTag.Columns);
Debug.Log("volume is " + xWidth + "mm x " + yHeight + "mm x " + zDepth + "mm, and located at " + pos[0] + "," + pos[1]);
vc.transform.localScale = new Vector3(xWidth / realWorldScale, yHeight / realWorldScale, zDepth / realWorldScale);
}
renderer.enabled = false;
is3D = true;
var rootDir = loadDicomInstance.rootDirectoryMap[record];
if (rootDir.Contains("ChestSeries"))
{
rootDir = "ChestSeries";
}
else if (rootDir.Contains("HeadNeckChestBody"))
{
rootDir = "HeadNeckChestBody";
}
var studyFID = record.LowerLevelDirectoryRecord.Get <string>(DicomTag.ReferencedFileID, 1);
var seriesFID = record.LowerLevelDirectoryRecord.Get <string>(DicomTag.ReferencedFileID, 2);
foreach (var s in loadDicomInstance.meshMarkers)
{
if (s.Contains(rootDir) && s.Contains(studyFID) && s.Contains(seriesFID))
{
Debug.Log(s + " seems to be a mesh marker for " + gameObject.name);
var meshes = Parabox.STL.pb_Stl_Importer.Import(s);
var marker = Instantiate(meshMarkerPrefab, meshMarkers.transform);
var meshFilter = marker.GetComponent <MeshFilter>();
meshFilter.mesh = meshes[0];
marker.name = s;
if (s.Contains("Sampling Site"))
{
marker.GetComponent <Renderer>().material.color = Color.green;
}
}
}
}
else if (button == "2D")
{
renderer.enabled = true;
transform.Find("zstack slider").gameObject.SetActive(true);
transform.Find("Volume Cube").gameObject.SetActive(false);
is3D = false;
}
}
