For various reasons, I needed to have a CT scan of my jaw done. This is a very high resolution 3D scan. The technician needed to review the scan to make sure it was of high quality and I stood behind him and looked over his shoulder. The software was pretty impressive, but the 3D model and resolution were really impressive. And then I left the office and drove home…
… and as I was driving, I thought: wouldn’t it be fun to have a copy of the data?; perhaps I could build a point cloud and shoot it into a crystal (as I’d done with fractals)? So I called back the lab (Park XRay) and asked if I could have a copy of the data. “Sure! It’s your skull.” was the reply and they delivered an extra copy to my dentist.
The files were in DICOM format and were produced or for-use by iCATVision. Fortunately, Python has a DICOM library, so it was fairly easy to parse the files. My code is on GitHub. [The code is not pretty, but it worked.]
I’ve previously “printed” point clouds into crystals using Precision Laser Art, so I needed to convert the 448 16-bit slices of my jaw into 1-bit XYZ point clouds. “visualize.py” provides a simple 2D visualization of the slices. Most importantly, it let me tune the threshold values for the quantizer so that the point cloud would highlight interesting structures in my jaw. Here’s the interface (perhaps it’s obvious, but I’m not a UX expert…):
Once I’d tuned the parameters, I added those parameters to “process.py” and generated the giant XYZ point cloud. The format of the point cloud is just:
X1 Y1 Z1
X2 Y2 Z2
X3 Y3 Z3
[repeat about 3 million times...]
I sent my order to Precision Laser Art and, after 7 days and $100, received this:
Which had a nicely cushioned interior:
And this is the resulting crystal. It’s the C-888 80mm cube.
While it’s not amazingly easy to see in this photo, my vertebrae and hyoid bone are clearly visible in the crystal.
Anyhow, the point is: medical data is cool. You can get it, so get it and play with it!