Extraction of the respiratory signal from small-animal CT projections for a retrospective gating method

Intrinsic Respiratory Gating for Simultaneous Multi-Mouse μCT Imaging to Assess Liver Tumors

Small animal micro computed tomography (μCT) is an important tool in cancer research and is used to quantify liver and lung tumors. A type of cancer that is intensively investigated with μCT is hepatocellular carcinoma (HCC). μCT scans acquire projections from different angles of the gantry which rotates X-ray source and detector around the animal. Motion of the animal causes inconsistencies between the projections which lead to artifacts in the resulting image. This is problematic in HCC research, where respiratory motion affects the image quality by causing hypodense intensity at the liver edge and smearing out small structures such as tumors. Dealing with respiratory motion is particularly difficult in a high throughput setting when multiple mice are scanned together and projection removal by retrospective respiratory gating may compromise image quality and dose efficiency. In mice, inhalation anesthesia leads to a regular respiration with short gasps and long phases of negligible motion. Using this effect and an iterative reconstruction which can cope with missing angles, we discard the relatively few projections in which the gasping motion occurs. Moreover, since gated acquisition, i.e., acquiring multiple projections from a single gantry angle is not a requirement, this method can be applied to existing scans. We applied our method in a high throughput setting in which four mice with HCC tumors were scanned simultaneously in a multi-mouse bed. To establish a ground truth, we manually selected projections with visible respiratory motion. Our automated intrinsic breathing projection selection achieved an accordance of 97% with manual selection. We reconstructed volumetric images and demonstrated that our intrinsic gating method significantly reduces the hypodense depiction at the cranial liver edge and improves the detectability of small tumors. Furthermore, we show that projection removal in a four mice scan discards only 7.5% more projections than in a single-mouse setting, i.e., four mouse scanning does not substantially compromise dose efficiency or image quality. To the best of our knowledge, no comparable method that combines multi-mouse scans for high throughput, intrinsic respiratory gating, and an available iterative reconstruction has been described for liver tumor imaging before.

Implementation and evaluation of respiratory gating in small animal radiotherapy

Major advance was done in preclinical radiotherapy thanks to the development of image guided micro-irradiator. Nevertheless, some applications still can benefit of improvements, such as the irradiation of mobile tumors. This preclinical radiotherapy case presents increased difficulties compared to clinical practice because of the waveform of small animals breathing cycle, its frequency and amplitude. To answer this issue, we developed a specific beam shutter and implemented respiratory gating on the X-RAD 225Cx preclinical irradiator. In the first step of this study, the shutter was accurately characterized. Opening and closing speed of 1.28 and 0.33 mm ms^-1 were respectively measured, and a transmission of 0.7% of the beam was measured with the shutter fully closed. Beam-on times were also determined for various gating parameters and highlighted a difference of 57 ms between the beam delivery duration and the gate width. This discrepancy was compensated during the respiratory monitoring adjustment. In a second step, a respiratory protocol was evaluated with two vertical beams of 2.5 and 5 mm diameters, for motion amplitudes ranging from 0.5 to 4 mm. This evaluation demonstrated the effectiveness of our set up to perform motion compensation for amplitude as small as 0.5 mm despite a dose gradient of 1.47 cGy mm^-1 observed with the 5 mm irradiation field, due to the shutter opening and closing durations. We also investigated the efficiency of a scintillating fiber dosimeter, adapted to small beams and providing real-time dose rate measurements. This detector showed very good performances to detect motion in small irradiation fields and would be very suitable to monitor the number of delivered gates until the planned delivered dose is achieved. This study presented a new respiratory gating set up and showed that very efficient motion compensation could be achieved in small animal radiotherapy.

ESTRO ACROP: Technology for precision small animal radiotherapy research: Optimal use and challenges

Many radiotherapy research centers have recently installed novel research platforms enabling the investigation of the radiation response of tumors and normal tissues in small animal models, possibly in combination with other treatment modalities. Many more research institutes are expected to follow in the coming years. These novel platforms are capable of mimicking human radiotherapy more closely than older technology. To facilitate the optimal use of these novel integrated precision irradiators and various small animal imaging devices, and to maximize the impact of the associated research, the ESTRO committee on coordinating guidelines ACROP (Advisory Committee in Radiation Oncology Practice) has commissioned a report to review the state of the art of the technology used in this new field of research, and to issue recommendations. This report discusses the combination of precision irradiation systems, small animal imaging (CT, MRI, PET, SPECT, bioluminescence) systems, image registration, treatment planning, and data processing. It also provides guidelines for reporting on studies.

Development of a dynamic phantom and investigation of mobile target imaging and irradiation in preclinical small animal research

OBJECTIVE

Progress made in preclinical radiotherapy makes respiratory gating reachable. Nevertheless, technical means are still needed, as well as accurate investigations of the effect of motion on small animal treatment plans.

METHODS

An animal-scaled dynamic phantom (0.3-11.1-mm motion peak-to-peak amplitude, 30-120 cycles per minute) was developed and characterized. It was used to evaluate respiratory monitoring and high resolution imaging (μPET/CT scans). The width and position variations of a fluorine-18 solution were measured for various motions and gating configurations. The phantom was finally used to measure the impact of motion on dose distribution for vertical irradiation using 2.5- and 5-mm collimations.

RESULTS

Phantom motions accurately reproduced original waveforms with good rate and amplitude linearity (R² = 1 and R² = 0.9995, respectively). µPET/CT acquisitions showed an increase of 92% of the target size caused by a 4.9-mm sine motion and reduced to <12% by gating. Target motion measurements showed consistency better than 18% between modalities. Irradiations showed that motions >0.8 and 1.1 mm (for the 2.5- and 5-mm collimations, respectively) significantly impact dose homogeneity in the target.

CONCLUSION

The phantom allowed studying motion in small animal imaging and irradiation. It showed the important impact of motions >2 mm and provided accurate data to improve the management of mobile tumour irradiation. The implementation of gated irradiation, associated with motion-compensated imaging, is currently under progress. Advances in knowledge: Small animal irradiation gating is not yet used in preclinical studies. As few solutions are under development, tools and accurate studies are highly needed.

A Novel Prior- and Motion-Based Compressed Sensing Method for Small-Animal Respiratory Gated CT

Low-dose protocols for respiratory gating in cardiothoracic small-animal imaging lead to streak artifacts in the images reconstructed with a Feldkamp-Davis-Kress (FDK) method. We propose a novel prior- and motion-based reconstruction (PRIMOR) method, which improves prior-based reconstruction (PBR) by adding a penalty function that includes a model of motion. The prior image is generated as the average of all the respiratory gates, reconstructed with FDK. Motion between respiratory gates is estimated using a nonrigid registration method based on hierarchical B-splines. We compare PRIMOR with an equivalent PBR method without motion estimation using as reference the reconstruction of high dose data. From these data acquired with a micro-CT scanner, different scenarios were simulated by changing photon flux and number of projections. Methods were evaluated in terms of contrast-to-noise-ratio (CNR), mean square error (MSE), streak artefact indicator (SAI), solution error norm (SEN), and correction of respiratory motion. Also, to evaluate the effect of each method on lung studies quantification, we have computed the Jaccard similarity index of the mask obtained from segmenting each image as compared to those obtained from the high dose reconstruction. Both iterative methods greatly improved FDK reconstruction in all cases. PBR was prone to streak artifacts and presented blurring effects in bone and lung tissues when using both a low number of projections and low dose. Adopting PBR as a reference, PRIMOR increased CNR up to 33% and decreased MSE, SAI and SEN up to 20%, 4% and 13%, respectively. PRIMOR also presented better compensation for respiratory motion and higher Jaccard similarity index. In conclusion, the new method proposed for low-dose respiratory gating in small-animal scanners shows an improvement in image quality and allows a reduction of dose or a reduction of the number of projections between two and three times with respect to previous PBR approaches.

Investigation of different sparsity transforms for the PICCS algorithm in small-animal respiratory gated CT

Respiratory gating helps to overcome the problem of breathing motion in cardiothoracic small-animal imaging by acquiring multiple images for each projection angle and then assigning projections to different phases. When this approach is used with a dose similar to that of a static acquisition, a low number of noisy projections are available for the reconstruction of each respiratory phase, thus leading to streak artifacts in the reconstructed images. This problem can be alleviated using a prior image constrained compressed sensing (PICCS) algorithm, which enables accurate reconstruction of highly undersampled data when a prior image is available. We compared variants of the PICCS algorithm with different transforms in the prior penalty function: gradient, unitary, and wavelet transform. In all cases the problem was solved using the Split Bregman approach, which is efficient for convex constrained optimization. The algorithms were evaluated using simulations generated from data previously acquired on a micro-CT scanner following a high-dose protocol (four times the dose of a standard static protocol). The resulting data were used to simulate scenarios with different dose levels and numbers of projections. All compressed sensing methods performed very similarly in terms of noise, spatiotemporal resolution, and streak reduction, and filtered back-projection was greatly improved. Nevertheless, the wavelet domain was found to be less prone to patchy cartoon-like artifacts than the commonly used gradient domain.

Evaluation of intrinsic respiratory signal determination methods for 4D CBCT adapted for mice

PURPOSE

4D CT imaging in mice is important in a variety of areas including studies of lung function and tumor motion. A necessary step in 4D imaging is obtaining a respiratory signal, which can be done through an external system or intrinsically through the projection images. A number of methods have been developed that can successfully determine the respiratory signal from cone-beam projection images of humans, however only a few have been utilized in a preclinical setting and most of these rely on step-and-shoot style imaging. The purpose of this work is to assess and make adaptions of several successful methods developed for humans for an image-guided preclinical radiation therapy system.

METHODS

Respiratory signals were determined from the projection images of free-breathing mice scanned on the X-RAD system using four methods: the so-called Amsterdam shroud method, a method based on the phase of the Fourier transform, a pixel intensity method, and a center of mass method. The Amsterdam shroud method was modified so the sharp inspiration peaks associated with anesthetized mouse breathing could be detected. Respiratory signals were used to sort projections into phase bins and 4D images were reconstructed. Error and standard deviation in the assignment of phase bins for the four methods compared to a manual method considered to be ground truth were calculated for a range of region of interest (ROI) sizes. Qualitative comparisons were additionally made between the 4D images obtained using each of the methods and the manual method.

RESULTS

4D images were successfully created for all mice with each of the respiratory signal extraction methods. Only minimal qualitative differences were noted between each of the methods and the manual method. The average error (and standard deviation) in phase bin assignment was 0.24 ± 0.08 (0.49 ± 0.11) phase bins for the Fourier transform method, 0.09 ± 0.03 (0.31 ± 0.08) phase bins for the modified Amsterdam shroud method, 0.09 ± 0.02 (0.33 ± 0.07) phase bins for the intensity method, and 0.37 ± 0.10 (0.57 ± 0.08) phase bins for the center of mass method. Little dependence on ROI size was noted for the modified Amsterdam shroud and intensity methods while the Fourier transform and center of mass methods showed a noticeable dependence on the ROI size.

CONCLUSIONS

The modified Amsterdam shroud, Fourier transform, and intensity respiratory signal methods are sufficiently accurate to be used for 4D imaging on the X-RAD system and show improvement over the existing center of mass method. The intensity and modified Amsterdam shroud methods are recommended due to their high accuracy and low dependence on ROI size.

Fully automated intrinsic respiratory and cardiac gating for small animal CT

A fully automated, intrinsic gating algorithm for small animal cone-beam CT is described and evaluated. A parameter representing the organ motion, derived from the raw projection images, is used for both cardiac and respiratory gating. The proposed algorithm makes it possible to reconstruct motion-corrected still images as well as to generate four-dimensional (4D) datasets representing the cardiac and pulmonary anatomy of free-breathing animals without the use of electrocardiogram (ECG) or respiratory sensors. Variation analysis of projections from several rotations is used to place a region of interest (ROI) on the diaphragm. The ROI is cranially extended to include the heart. The centre of mass (COM) variation within this ROI, the filtered frequency response and the local maxima are used to derive a binary motion-gating parameter for phase-sensitive gated reconstruction. This algorithm was implemented on a flat-panel-based cone-beam CT scanner and evaluated using a moving phantom and animal scans (seven rats and eight mice). Volumes were determined using a semiautomatic segmentation. In all cases robust gating signals could be obtained. The maximum volume error in phantom studies was less than 6%. By utilizing extrinsic gating via externally placed cardiac and respiratory sensors, the functional parameters (e.g. cardiac ejection fraction) and image quality were equivalent to this current gold standard. This algorithm obviates the necessity of both gating hardware and user interaction. The simplicity of the proposed algorithm enables adoption in a wide range of small animal cone-beam CT scanners.