Biomedical subjects

Parag J Parikh

Publications and source records attributed to Parag J Parikh.

8 recordsLinked to original sources

PubMed2006

A prospective study of differences in duodenum compared to remaining small bowel motion between radiation treatments: implications for radiation dose escalation in carcinoma of the pancreas.

PURPOSE: As a foundation for a dose escalation trial, we sought to characterize duodenal and non-duodenal small bowel organ motion between fractions of pancreatic radiation therapy. PATIENTS AND METHODS: Nine patients (4 women, 5 men) undergoing radiation therapy were enrolled in this prospective study. The patients had up to four weekly CT scans performed during their course of radiation therapy. Pancreas, duodenum and non-duodenal small bowel were then contoured for each CT scan. On the initial scan, a four-field plan was generated to fully cover the pancreas. This plan was registered to each subsequent CT scan. Dose-volume histogram (DVH) analyses were performed for the duodenum, non-duodenal small bowel, large bowel, and pancreas. RESULTS: With significant individual variation, the volume of duodenum receiving at least 80% of the prescribed dose was consistently greater than the remaining small bowel. In the patient with the largest inter-fraction variation, the fractional volume of non-duodenal small bowel irradiated to at least the 80% isodose line ranged from 1% to 20%. In the patient with the largest inter-fraction variation, the fractional volume of duodenum irradiated to at least the 80% isodose line ranged from 30% to 100%. CONCLUSION: The volume of small bowel irradiated during four-field pancreatic radiation therapy changes substantially between fractions. This suggests dose escalation may be possible. However, dose limits to the duodenum should be stricter than for other segments of small bowel.

Dose-Response Relationship, Radiation↗

PubMed2006

A comparison between amplitude sorting and phase-angle sorting using external respiratory measurement for 4D CT.

Respiratory motion can cause significant dose delivery errors in conformal radiation therapy for thoracic and upper abdominal tumors. Four-dimensional computed tomography (4D CT) has been proposed to provide the image data necessary to model tumor motion and consequently reduce these errors. The purpose of this work was to compare 4D CT reconstruction methods using amplitude sorting and phase angle sorting. A 16-slice CT scanner was operated in ciné mode to acquire 25 scans consecutively at each couch position through the thorax. The patient underwent synchronized external respiratory measurements. The scans were sorted into 12 phases based, respectively, on the amplitude and direction (inhalation or exhalation) or on the phase angle (0-360 degrees) of the external respiratory signal. With the assumption that lung motion is largely proportional to the measured respiratory amplitude, the variation in amplitude corresponds to the variation in motion for each phase. A smaller variation in amplitude would associate with an improved reconstructed image. Air content, defined as the amount of air within the lungs, bronchi, and trachea in a 16-slice CT segment and used by our group as a surrogate for internal motion, was correlated to the respiratory amplitude and phase angle throughout the lungs. For the 35 patients who underwent quiet breathing, images (similar to those used for treatment planning) and animations (used to display respiratory motion) generated using amplitude sorting displayed fewer reconstruction artifacts than those generated using phase angle sorting. The variations in respiratory amplitude were significantly smaller (P < 0.001) with amplitude sorting than those with phase angle sorting. The subdivision of the breathing cycle into more (finer) phases improved the consistency in respiratory amplitude for amplitude sorting, but not for phase angle sorting. For 33 of the 35 patients, the air content showed significantly improved (P < 0.001) correlation with the respiratory amplitude than with the phase angle, suggesting a stronger relationship between internal motion and amplitude. Overall, amplitude sorting performed better than phase angle sorting for 33 of the 35 patients and equally well for two patients who were immobilized with a stereotactic body frame and an abdominal compression plate.

Algorithms↗

PubMed2006

A semi-automatic method for peak and valley detection in free-breathing respiratory waveforms.

The existing commercial software often inadequately determines respiratory peaks for patients in respiration correlated computed tomography. A semi-automatic method was developed for peak and valley detection in free-breathing respiratory waveforms. First the waveform is separated into breath cycles by identifying intercepts of a moving average curve with the inspiration and expiration branches of the waveform. Peaks and valleys were then defined, respectively, as the maximum and minimum between pairs of alternating inspiration and expiration intercepts. Finally, automatic corrections and manual user interventions were employed. On average for each of the 20 patients, 99% of 307 peaks and valleys were automatically detected in 2.8 s. This method was robust for bellows waveforms with large variations.

Algorithms↗

PubMed2006

Progress toward a microradiation therapy small animal conformal irradiator.

Microradiation therapy (microRT) systems are being designed to provide conformal radiation therapy to small animals enabling quantitative radiation response evaluation. We used a Monte Carlo approach to estimate the radiation dose distributions from proposed blueprints and developed a beam model to aid in the microRT system design process. This process was applied to a prototype irradiator that uses a small (3 mm long and 3 mm in diameter), cylindrical, high-activity 192Ir source delivering the radiation beam using custom-fabricated tungsten collimators. The BEAMnrc Monte Carlo code was used to simulate dose distributions from these prototype collimators. Simulations were performed at three source-to-surface distances (50, 60, and 70 mm), and with five circular field sizes (5, 7.5, 10, 12.5, and 15 mm). A dose to a 50 X 50 X 50 mm3 water phantom with 1 X 1 X 1 mm3 voxel spacing was computed. A multiparameter dose calculation algorithm was developed to efficiently and accurately calculate doses for treatment planning exercises. The parametrization was selected so that the parameters varied smoothly as a function of depth, source-to-surface distance, and field size, allowing interpolation for geometries that were not simulated using the Monte Carlo simulation. Direct comparison of the model with the Monte Carlo simulations showed that the variations were within 5% error for field sizes larger than 10 mm, and up to 10% for smaller field sizes.

Animals↗

PubMed2005

Novel breathing motion model for radiotherapy.

PURPOSE: An accurate model of breathing motion under quiet respiration is desirable to obtain the most accurate and conformal dose distributions for mobile lung cancer lesions. On the basis of recent lung motion measurements and the physiologic functioning of the lungs, we have determined that the motion of lung and lung tumor tissues can be modeled as a function of five degrees of freedom, the position of the tissues at a user-specified reference breathing phase, tidal volume and its temporal derivative airflow (tidal volume phase space). Time is an implicit variable in this model. METHODS AND MATERIALS: To test this hypothesis, a mathematical model of motion was developed that described the motion of objects p in the lungs as linear functions of tidal volume and airflow. The position of an object was described relative to its position -->P0 at the reference tidal volume and zero airflow, and the motion of the object was referenced to this position. Hysteresis behavior was hypothesized to be caused by pressure imbalances in the lung during breathing and was, in this model, a function of airflow. The motion was modeled as independent tidal volume and airflow displacement vectors, with the position of the object at time t equal to the vector sum -->rP(t) = -->rv(t) + -->rf(t) where -->rv(t) and -->rf(t) were displacement vectors with magnitudes approximated by linear functions of the tidal volume and airflow. To test this model, we analyzed five-dimensional CT scans (CT scans acquired with simultaneous real-time monitoring of the tidal volume) of 4 patients. The scans were acquired throughout the lungs, but the trajectories were analyzed in the couch positions near the diaphragm. A template-matching algorithm was implemented to identify the positions of the points throughout the 15 scans. In total, 76 points throughout the 4 patients were tracked. The lateral motion of these points was minimal; thus, the model was described in two spatial dimensions, with a total of six parameters necessary to describe the 30 degrees of freedom inherent in the 15 positions. RESULTS: For the 76 evaluated points, the average discrepancy (the distance between the measured and prediction positions) of the 15 locations for each tracked point was 0.75 +/- 0.25 mm, with an average maximal discrepancy of 1.55 +/- 0.54 mm. The average discrepancy was also tabulated as a fraction of the breathing motion. Discrepancies of <10% and 15% of the overall motion occurred in 73% and 95% of the tracked points, respectively. CONCLUSION: The motion tracking algorithms are being improved and automated to provide more motion data to test the models. This may allow a measurement of the motion-fitting parameters throughout the lungs. If the parameters vary smoothly, interpolation may be possible, yielding a continuous mathematical model of the breathing motion throughout the lungs. The utility of the model will depend on its stability as a function of time. If the model is only robust during the measurement session, it may be useful for determining lung function. If it is robust for weeks, it may be useful for treatment planning and gating of lung treatments. The use of tidal volume phase space for characterizing breathing motion appears to have provided, for the first time, the potential for a patient-specific mathematical model of breathing motion.

Humans↗

PubMed2005

Quantitation of the reconstruction quality of a four-dimensional computed tomography process for lung cancer patients.

We have developed a four-dimensional computed tomography (4D CT) technique for mapping breathing motion in radiotherapy treatment planning. A multislice CT scanner (1.5 mm slices) operated in ciné mode was used to acquire 12 contiguous slices in each couch position for 15 consecutive scans (0.5 s rotation, 0.25 s between scans) while the patient underwent simultaneous quantitative spirometry measurements to provide a sorting metric. The spirometry-sorted scans were used to reconstruct a 4D data set. A critical factor for 4D CT is quantifying the reconstructed data set quality which we measure by correlating the metric used relative to internal-object motion. For this study, the internal air content within the lung was used as a surrogate for internal motion measurements. Thresholding and image morphological operations were applied to delineate the air-containing tissues (lungs, trachea) from each CT slice. The Hounsfield values were converted to the internal air content (V). The relationship between the air content and spirometer-measured tidal volume (v) was found to be quite linear throughout the lungs and was used to estimate the overall accuracy and precision of tidal volume-sorted 4D CT. Inspection of the CT-scan air content as a function of tidal volume showed excellent correlations (typically r>0.99) throughout the lung volume. Because of the discovered linear relationship, the ratio of internal air content to tidal volume was indicative of the fraction of air change in each couch position. Theoretically, due to air density differences within the lung and in room, the sum of these ratios would equal 1.11. For 12 patients, the mean value was 1.08 +/- 0.06, indicating the high quality of spirometry-based image sorting. The residual of a first-order fit between v and V was used to estimate the process precision. For all patients, the precision was better than 8%, with a mean value of 5.1% +/- 1.9%. This quantitative analysis highlights the value of using spirometry as the metric in sorting CT scans. The 4D reconstruction provides the CT data required to measure the three-dimensional trajectory of tumor and lung tissue during free breathing.

Air↗

PubMed2005

Comparison of spirometry and abdominal height as four-dimensional computed tomography metrics in lung.

An important consideration in four-dimensional CT scanning is the selection of a breathing metric for sorting the CT data and modeling internal motion. This study compared two noninvasive breathing metrics, spirometry and abdominal height, against internal air content, used as a surrogate for internal motion. Both metrics were shown to be accurate, but the spirometry showed a stronger and more reproducible relationship than the abdominal height in the lung. The abdominal height was known to be affected by sensor placement and patient positioning while the spirometer exhibited signal drift. By combining these two, a normalization of the drift-free metric to tidal volume may be generated and the overall metric precision may be improved.

Abdomen↗

PubMed2005

Reduction of motion blurring artifacts using respiratory gated CT in sinogram space: a quantitative evaluation.

Techniques have been developed for reducing motion blurring artifacts by using respiratory gated computed tomography (CT) in sinogram space and quantitatively evaluating the artifact reduction. A synthetic sinogram was built from multiple scans intercepting a respiratory gating window. A gated CT image was then reconstructed using the filtered back-projection algorithm. Wedge phantoms, developed for quantifying the motion artifact reduction, were scanned while being moved using a computer-controlled linear stage. The resulting artifacts appeared between the high and low density regions as an apparent feature with a Hounsfield value that was the average of the two regions. A CT profile through these regions was fit using two error functions, each modeling the partial-volume averaging characteristics for the unmoving phantom. The motion artifact was quantified by determining the apparent distance between the two functions. The blurring artifact had a linear relationship with both the speed and the tangent of the wedge angles. When gating was employed, the blurring artifact was reduced systematically at the air-phantom interface. The gated image of phantoms moving at 20 mm/s showed similar blurring artifacts as the nongated image of phantoms moving at 10 mm/s. Nine patients were also scanned using the synchronized respiratory motion technique. Image artifacts were evaluated in the diaphragm, where high contrast interfaces intercepted the imaging plane. For patients, this respiratory gating technique reduced the blurring artifacts by 9%-41% at the lung-diaphragm interface.

Air↗