Table of Contents

CalibrateDSX Overview

CalibrateDSX prepares the X-ray images and calculates the configuration of the X-ray hardware. It performs three important tasks:

     1. X-ray Correction
X-ray images typically suffer from uniformity and distortion errors due to cross-scatter, non-uniform beam intensity, and objects (other than the subject) in the field of view. CalibrateDSX corrects these errors using images taken with an empty field of view (“white” images) and images taken of a uniform grid of holes or beads at known locations (“grid” images). Information on collecting white and grid images can be found here. X-ray images of the calibration object and of the subject are corrected by first pre-processing the white and grid images as follows:
Stacks of white images are processed via these steps:

- remove outlier images (based on the mean intensity)
- average the remaining images
- add an optional border to the image, based on the value of the Image Padding parameter
- resize the image to the size specified by the Image Size parameter
- save the resulting image in the TIFF format

Stacks of grid images are processed via these steps:

- remove outlier images (based on the mean intensity)
- add an optional border to the images, based on the value of the Image Padding parameter
- resize the images to the size specified by the Image Size parameter
- uniformity-correct the remaining images: output image = white image - grid image
- save the resulting image stack in the TIFF format
- calculate a displacement map for the first image
- calculate initial pixel sizes from the displacements in the center of the map

Once the white and grid images have been processed, all other X-ray images (calibration object, static subject, motion subject, etc.) can be corrected with these steps:

- remove outlier images (calibration images only, based on the mean intensity)
- optionally smooth the remaining images with a Gaussian filter
- add an optional border to the images, based on the value of the Image Padding parameter
- resize the images to the size specified by the Image Size parameter
- uniformity-correct the remaining images: output image = white image - input image
- distortion-correct the remaining images: use the average displacement map to move each pixel to its proper location
- saving the resulting image stack in the TIFF format

Note: uniformity correction and distortion correction are both optional. If no white images are specified in the calibration data, uniformity correction will be skipped and the images to be corrected will simply be inverted. Similarly, if no grid images are specified, distortion correction will be skipped.
     2. 3D X-ray Configuration
CalibrateDSX calculates the 3D configuration of the X-ray hardware (the Pose of the X-ray sources and image planes) from images of a rigid calibration object with precisely placed embedded radio opaque beads. Given corrected images of the calibration object, the 3D configuration of the X-ray hardware can be computed. The first step is to digitize the beads of the object in each X-ray view. Identification of the beads can be difficult if there are many of them in the object. A commonly used object contains 64 beads, but only 40 - 50 are visible in each view. To aid identification, uniquely shaped trinkets are sometimes added to the object. CalibrateDSX enables digitizing of the trinkets, or, if there aren’t any, 4 or 5 beads that are easy to identify. From these 4 - 5 trinket/bead locations, the software can calculate where the other beads are expected to be, and locate and identify their centroids on the images. The user can check these automatically digitized beads and manually digitize any that were misidentified. With this step complete, CalibrateDSX now has, for each X-ray view, a set of 2D bead locations in the image plane and a set of 3D bead locations in the lab frame. A simulated annealing algorithm is used for the DLT (direct linear transform) calculations, outputting the following parameters:

the XYZ position of the X-ray source in the lab frame

the XYZ position of the geometric center of the X-ray image (not the perpendicular point of the X-ray source) in the lab frame

the XYZ Euler angles (or 3×3 transform) expressing the orientation of the X-ray image in the lab frame

the XY pixel size of the X-ray image
     3. Motion Capture Transform
If traditional motion capture data (e.g., surface markers, inertial sensors) are collected in conjunction with X-ray data, they can be used in X4D to speed the tracking process by determining initial poses of the bones. CalibrateDSX calculates the transform between the X-ray (lab) frame and the frame of the motion capture system (mocap frame) so that the mocap-based bone poses can be converted to the lab frame. It does this using surface markers fixed to the calibration object. The 3D positions of these markers are known in the lab frame, which is embedded in the object. The motion capture system records their 3D positions in the mocap frame, and thus the transform between the lab and mocap frames can be computed.

Tutorials

How To: Correct X-ray Images and Perform 3D Calibration
How To: Collect Calibration Images
How To: Digitize Points of Interest
The Importance of Image Smoothing and Downsampling
[Video Tutorial on Calibration]
File Menu

View Menu

Tools Menu

Mocap Menu

Options Menu

Help Menu

Widgets

DSX Configuration
Calibration Object
3D Calibration
Image Correction
Image Preview
2D X-ray Windows
Settings Dialog