Effect of loading direction and anatomical location on the ultimate tensile stress, fracture toughness, and failure patterns of knee meniscus.

BACKGROUND: Rupture of the knee menisci is a common injury that can have implications for other conditions, such as osteoarthritis. The fracture toughness of soft tissue (Jc) is a mechanical property that characterizes its resistance to tear extension. To date, Jc of the meniscus has not been quantified. METHODS: Cyclic tensile tests were conducted on meniscus samples to determine Jc and explore its characteristics. Initially, the study investigated the impact of an initial notch on the ultimate tensile stress. This allowed for an understanding of how the presence of a notch affects its structural integrity. Subsequently, Jc was measured in both the radial and circumferential directions to assess its loading direction dependency. Furthermore, the study assessed the effect of anatomical location by comparing samples collected from the femoral and tibial layers. RESULTS: Defect tolerance of the meniscus is influenced by the loading direction. In the circumferential direction, the presence of an initial notch did not affect the ultimate stress, and no crack expansion was observed. In radial samples with a notch length of 40% or more of the total width, crack propagation occurred, leading to a decrease in the ultimate stress (p< 0.01). Additionally, Jc was found to be higher in the femoral layer compared to the tibial layer (p= 0.017). CONCLUSION: The study also examined the failure patterns of the meniscus to enhance our understanding of its pathology. These insights contribute to a better comprehension of meniscus injuries and can aid in the development of more effective treatment strategies.

Dosimetric Evaluation and Reproducibility of Breath-hold Plans in Intensity Modulated Proton Therapy: An Initial Clinical Experience.

Purpose: Breath-hold (BH) technique can mitigate target motion, minimize target margins, reduce normal tissue doses, and lower the effect of interplay effects with intensity-modulated proton therapy (IMPT). This study presents dosimetric comparisons between BH and nonbreath-hold (non-BH) IMPT plans and investigates the reproducibility of BH plans using frequent quality assurance (QA) computed tomography scans (CT). Methods and Materials: Data from 77 consecutive patients with liver (n = 32), mediastinal/lung (n = 21), nonliver upper abdomen (n = 20), and malignancies in the gastroesophageal junction (n = 4), that were treated with a BH spirometry system (SDX) were evaluated. All patients underwent both BH CT and 4-dimensional CT simulations. Clinically acceptable BH and non-BH plans were generated on each scan, and dose-volume histograms of the 2 plans were compared. Reproducibility of the BH plans for 30 consecutive patients was assessed using 1 to 3 QA CTs per patient and variations in dose-volume histograms for deformed target and organs at risk (OARs) volumes were compared with the initial CT plan. Results: Use of BH scans reduced initial and boost target volumes to 72% ± 20% and 70% ± 17% of non-BH volumes, respectively. Additionally, mean dose to liver, stomach, kidney, esophagus, heart, and lung V20 were each reduced to 71% to 79% with the BH technique. Similarly, small and large bowels, heart, and spinal cord maximum doses were each lowered to 68% to 84%. Analysis of 62 QA CT scans demonstrated that mean target and OAR doses using BH scans were reproducible to within 5% of their nominal plan values. Conclusions: The BH technique reduces the irradiated volume, leading to clinically significant reductions in OAR doses. By mitigating tumor motion, the BH technique leads to reproducible target coverage and OAR doses. Its use can reduce motion-related uncertainties that are normally associated with the treatment of thoracic and abdominal tumors and, therefore, optimize IMPT delivery.

Overview and Recommendations for Prospective Multi-institutional Spatially Fractionated Radiation Therapy Clinical Trials.

PURPOSE: The highly heterogeneous dose delivery of spatially fractionated radiation therapy (SFRT) is a profound departure from standard radiation planning and reporting approaches. Early SFRT studies have shown excellent clinical outcomes. However, prospective multi-institutional clinical trials of SFRT are still lacking. This NRG Oncology/American Association of Physicists in Medicine working group consensus aimed to develop recommendations on dosimetric planning, delivery, and SFRT dose reporting to address this current obstacle toward the design of SFRT clinical trials. METHODS AND MATERIALS: Working groups consisting of radiation oncologists, radiobiologists, and medical physicists with expertise in SFRT were formed in NRG Oncology and the American Association of Physicists in Medicine to investigate the needs and barriers in SFRT clinical trials. RESULTS: Upon reviewing the SFRT technologies and methods, this group identified challenges in several areas, including the availability of SFRT, the lack of treatment planning system support for SFRT, the lack of guidance in the physics and dosimetry of SFRT, the approximated radiobiological modeling of SFRT, and the prescription and combination of SFRT with conventional radiation therapy. CONCLUSIONS: Recognizing these challenges, the group further recommended several areas of improvement for the application of SFRT in cancer treatment, including the creation of clinical practice guidance documents, the improvement of treatment planning system support, the generation of treatment planning and dosimetric index reporting templates, and the development of better radiobiological models through preclinical studies and through conducting multi-institution clinical trials.

Proton Therapy for Spinal Tumors: A Consensus Statement From the Particle Therapy Cooperative Group.

PURPOSE: Proton beam therapy (PBT) plays an important role in the management of primary spine tumors. The purpose of this consensus statement was to summarize safe and optimal delivery of PBT for spinal tumors. METHODS AND MATERIALS: The Particle Therapy Cooperative Group Skull Base/Central nervous system/Sarcoma Subcommittee consisting of radiation oncologists and medical physicists with specific expertise in spinal irradiation developed expert recommendations discussing treatment planning considerations and current approaches in the treatment of primary spinal tumors. RESULTS: Computed tomography simulation: factors that require significant consideration include (1) patient comfort, (2) setup reproducibility and stability, and (3) accessibility of appropriate beam angles. SPINE STABILIZATION HARDWARE: If present, hardware should be placed with cross-links well above/below the level of the primary tumor to reduce the metal burden at the level of the tumor bed. New materials that can reduce uncertainties include polyether-ether-ketone and composite polyether-ether-ketone-carbon fiber implants. FIELD ARRANGEMENT: Appropriate beam selection is required to ensure robust target coverage and organ at risk sparing. Commonly, 2 to 4 treatment fields, typically from posterior and/or posterior-oblique directions, are used. TREATMENT PLANNING METHODOLOGY: Robust optimization is recommended for all pencil beam scanning plans (the preferred treatment modality) and should consider setup uncertainty (between 3 and 7 mm) and range uncertainty (3%-3.5%). In the presence of metal hardware, use of an increased range uncertainty up to 5% is recommended. CONCLUSIONS: The Particle Therapy Cooperative Group Skull Base/Central nervous system/Sarcoma Subcommittee has developed recommendations to enable centers to deliver PBT safely and effectively for the management of primary spinal tumors.

Design and performance evaluation of a slit-slat camera for 2D prompt gamma imaging in proton therapy monitoring: A Monte Carlo simulation study.

PURPOSE: We investigated the design of a prompt gamma camera for real-time dose delivery verification and the partial mitigation of range uncertainties. METHODS: A slit slat (SS) camera was optimized using the trade-off between the signal-to-noise ratio and spatial resolution. Then, using the GATE Monte Carlo package, the camera performances were estimated by means of target shifts, beam position quantification, changing the camera distance from the beam, and air cavity inserting. A homogeneous PMMA phantom and the air gaps induced PMMA phantom were used. The air gaps ranged from 5 mm to 30 mm by 5 mm increments were positioned in the middle of the beam range. To reduce the simulation time, phase space scoring was used. The batch method with five realizations was used for stochastic error calculations. RESULTS: The system's detection efficiency was 1.1 × 10 - 4 PGs Emitted PGs ( 1.8 × 10 - 5 $1.1 \times {10}^{-4}\frac{{\rm PGs}}{{\rm Emitted}\ {\rm PGs}}\ (1.8 \times {10}^{-5}$ PGs/proton) for a 10 × 20 cm2 detector (source-to-collimator distance = 15.0 cm). Axial and transaxial resolutions were 23 mm and 18 mm, respectively. The SS camera estimated the range as 69.0 ± 3.4 (relative stochastic error 1-sigma is 5%) and 67.6 ± 1.8 mm (2.6%) for the real range of 67.0 mm for 107 and 108 protons of 100 MeV, respectively. Considering 160 MeV, these values are 155.5 ± 3.1 (2%) and 152.2 ± 2.0 mm (1.3%) for the real range of 152.0 mm for 107 and 108 protons, respectively. Considering phantom shift, for a 100 MeV beam, the precision of the quantification (1-sigma) in the axial and lateral phantom shift estimation is 2.6 mm and 1 mm, respectively. Accordingly, the axial and lateral quantification precisions were 1.3 mm and 1 mm for a 160 MeV beam, respectively. Furthermore, the quantification of an air gap formulated as gap d e t = 0.98 × gap real ${{\rm gap}}_{det}=0.98 \times {{\rm gap}}_{{\rm real}}$ , where gap d e t ${{\rm gap}}_{det}$ and gapreal are the estimated and real air gap, respectively. The precision of the air gap quantification is 1.6 mm (1 sigma). Moreover, 2D PG images show the trajectory of the proton beam through the phantom. CONCLUSION: The proposed slit-slat imaging systems can potentially provide a real-time, in-vivo, and non-invasive treatment monitoring method for proton therapy.

Cleaning the dose falloff with low modulation in SBRT lung plans.

Purpose.This dosimetric study is intended to lower the modulation factor in lung SBRT plans generated in the Eclipse TPS that could replace highly modulated plans that are prone to the interplay effect.Materials and methods.Twenty clinical lung SBRT plans with high modulation factors (≥4) were replanned in Varian Eclipse TPS version 15.5 utilizing 2 mm craniocaudal and 1 mm axial block margins followed by light optimization in order to reduce modulation. A unique plan optimization methodology, which utilizes a novel shell structure (OptiForR50) for R50%optimization in addition to five consecutive concentric 5 mm shells, was utilized to control dose falloff according to RTOG 0813 and 0915 recommendations. The prescription varied from 34-54 Gy in 1-4 fractions, and the dose objectives were PTV D95%= Rx, PTV Dmax< 140% of Rx, and minimizing the modulation factor. Plan evaluation metrics included modulation factor, CIRTOG, homogeneity index (HI), R50%, D2cm, V105%, and lung V8-12.8Gy(Timmerman Constraint). A random-intercept linear mixed effects model was used with a p ≤ 0.05 threshold to test for statistical significance.Results.The retrospectively generated plans had significantly lower modulation factors (3.65 ± 0.35 versus 4.59 ± 0.54; p < 0.001), lower CIRTOG(0.97 ± 0.02 versus 1.02 ± 0.06; p = 0.001), higher HI (1.35 ± 0.06 versus 1.14 ± 0.04; p < 0.001), lower R50%(4.09 ± 0.45 versus 4.56 ± 0.56; p < 0.001), and lower lungs V8-12.8Gy(Timmerman) (4.61% ± 3.18% versus 4.92% ± 3.37%; p < 0.001). The high dose spillage V105%was borderline significantly lower (0.44% ± 0.49% versus 1.10% ± 1.64%; p = 0.051). The D2cmwas not statistically different (46.06% ± 4.01% versus 46.19% ± 2.80%; p = 0.835).Conclusion.Lung SBRT plans with significantly lower modulation factors can be generated that meet the RTOG constraints, using our planning strategy.

A phantom-based study and clinical implementation of brainlab's treatment planning system for radiosurgical treatments of arteriovenous malformations.

PURPOSE: Development of a simple, phantom-based methodology allowing for pilot applications for the Elements TPS cranio-vascular module and clinical implementation prior to AVM treatments. METHODS: A customized phantom was developed to be visible in MRI and CT images. High resolution digital subtraction angiograms (DSAs) and CT images of the phantom were acquired and imported into the Brainlab Elements treatment planning system. A clinical treatment plan with 5 arcs was generated in cranial vascular planning module and delivered to the phantom using a Varian TrueBeam STx Linac equipped with HD-MLCs and Brainlab ExacTrac imaging system for non-coplanar setup verification. The delivered dose was verified using a calibrated ionization chamber placed in the phantom. Upon verification of the TPS workflow, three patients with AVM who have been treated to date at our center using the Brainlab's cranial vascular module for AVM are presented here for retrospective review. RESULTS: The difference between the planed and measured dose by the ionization chamber was found to be less than 1%. Following a successful dose verification study, a clinical workflow was created. Currently, three AVM patients have been treated successfully. Clinical aspects of imaging and treatment planning consideration are presented in retrospective setting. CONCLUSIONS: Dose verification of the Brainlab Elements cranial vascular planning module for intracranial SRS treatments of AVM on Varian TrueBeam was successfully implemented using a custom-made phantom with <1% discrepancy. The Brainlab Elements' cranial vascular module was successfully implemented in clinical workflow to treat patients with AVM. This manuscript provides a guideline for clinical implementation of frameless Linac-based AVM treatment using the Brainlab Elements TPS.

Influence of crack length and anatomical location on the fracture toughness of annulus fibrosus.

Fracture toughness (Jc) of a soft biological tissue is an important mechanical property that characterizes its resistance to crack or tear extension. To date, no information is available on fracture toughness of annulus fibrosus (AF); therefore, its defect tolerance is not known. The present study modified a previously introduced method to determine Jc of ovine AF. Then, the effect of the notch length on the failure pattern and Jc was investigated. Also, the test samples of anterior and lateral regions were collected to determine the effect of the location on Jc. Results showed that for a notch length of less than 45% of total width, no crack extension occurred, but for a notch length above 45% of the width, crack propagation and ultimately the failure of the AF were observed. However, statistical analysis indicated no significant difference on Jc (p = 0.5) for the initial notch length of 50% and 70% of total width. The fracture toughness was significantly higher for the samples extracted from the lateral site than those from the anterior site (p < 0.05). Dissimilar failure patterns were observed for different initial notch lengths. Among the parameters studied, the defect tolerance of AF was dependent on the initial tear size.

Simulation of an HDR "Boost" with Stereotactic Proton versus Photon Therapy in Prostate Cancer: A Dosimetric Feasibility Study.

PURPOSE/OBJECTIVES: To compare the dose escalation potential of stereotactic body proton therapy (SBPT) versus stereotactic body photon therapy (SBXT) using high-dose rate prostate brachytherapy (HDR-B) dose-prescription metrics. PATIENTS AND METHODS: Twenty-five patients previously treated with radiation for prostate cancer were identified and stratified by prostate size (≤ 50cc; n = 13, > 50cc; n = 12). Initial CT simulation scans were re-planned using SBXT and SBPT modalities using a prescription dose of 19Gy in 2 fractions. Target coverage goals were designed to mimic the dose distributions of HDR-B and maximized to the upper limit constraint for the rectum and urethra. Dosimetric parameters between SBPT and SBXT were compared using the signed-rank test and again after stratification for prostate size (≤ 50cm3 and >50cm3) using the Wilcoxon rank test. RESULTS: Prostate volume receiving 100% of the dose (V100) was significantly greater for SBXT (99%) versus SBPT (96%) (P ≤ 0.01), whereas the median V125 (82% vs. 73%, P < 0.01) and V200 (12% vs. 2%, P < 0.01) was significantly greater for SBPT compared to SBXT. Median V150 was 49% for both cohorts (P = 0.92). V125 and V200 were significantly correlated with prostate size. For prostates > 50cm3, V200 was significantly greater with SBPT compared to SBXT (14.5% vs. 1%, P = 0.005), but not for prostates 50cm3 (9% vs 4%, P = 0.11). Median dose to 2cm3 of the bladder neck was significantly lower with SBPT versus SBXT (9.6 Gy vs. 14 Gy, P < 0.01). CONCLUSION: SBPT and SBXT can be used to simulate an HDR-B boost for locally advanced prostate cancer. SBPT demonstrated greater dose escalation potential than SBXT. These results are relevant for future trial design, particularly in patients with high risk prostate cancer who are not amenable to brachytherapy.

Plan quality effects of maximum monitor unit constraints in pencil beam scanning proton therapy for central nervous system and skull base tumors.

PURPOSE/OBJECTIVE(S): With reports of CNS toxicity in patients treated with proton therapy at doses lower than would be expected based on photon data, it has been proposed that heavy monitor unit (MU) weighting of pencil beam scanning (PBS) proton therapy spots may potentially increase the risk of toxicity. We evaluated the impact of maximum MU weighting per spot (maxMU/spot) restrictions on PBS plan quality, prior to implementing clinic-wide maxMU/spot restrictions. MATERIALS/METHODS: PBS plans of 11 patients, of which 3 plans included boosts, for a total of 14 PBS sample cases were included. Per sample case, a single dosimetrist created 4 test plans, gradually reducing the maxMU/spot in the plan. Test Plan 1, unrestricted in maxMU/spot, was the reference for all restricted plan comparisons (comparison sets 2 vs. 1; 3 vs. 1; and 4 vs. 1). The impact of MU/spot restrictions on plan quality metrics were analyzed with Wilcoxon signed rank test analyses. Treatment delivery time was modeled for a representative case. RESULTS: A total of 14 PBS sample cases, 7 (50%) single-field optimized, 7 (50%) multi-field optimized, 9 (64%) delivering > 3500 cGy, 9 (64%) with 3 beams, and 7 (50%) without a range shifter were included. There were no differences in plan quality metrics of target coverage (V95% and V100% prescription), conformality and gradient indices, hot spot volume (V105% prescription), and dose to normal brain (V10%/30%/50%/70%/90%/100% prescription) with reductions of allowable maxMU/spot across all comparison sets (p > 0.05). Max MU/spot restrictions did not increase treatment delivery time when analyzed for a representative case. CONCLUSION: MaxMU/spot restrictions within the thresholds evaluated in this study did not degrade overall plan quality metrics. Future studies should evaluate spot weighting with linear energy transfer/relative biologic effectiveness-informed planning to determine if spot weighting manipulation impacts clinical outcomes and mitigates toxicity.

Initial Validation of Proton Dose Calculations on SPR Images from DECT in Treatment Planning System.

PURPOSE: To investigate and quantify the potential benefits associated with the use of stopping-power-ratio (SPR) images created from dual-energy computed tomography (DECT) images for proton dose calculation in a clinical proton treatment planning system (TPS). MATERIALS AND METHODS: The DECT and single-energy computed tomography (SECT) scans obtained for 26 plastic tissue surrogate plugs were placed individually in a tissue-equivalent plastic phantom. Relative-electron density (ρe) and effective atomic number (Z eff) images were reconstructed from the DECT scans and used to create an SPR image set for each plug. Next, the SPR for each plug was measured in a clinical proton beam for comparison of the calculated values in the SPR images. The SPR images and SECTs were then imported into a clinical TPS, and treatment plans were developed consisting of a single field delivering a 10 × 10 × 10-cm3 spread-out Bragg peak to a clinical target volume that contained the plugs. To verify the accuracy of the TPS dose calculated from the SPR images and SECTs, treatment plans were delivered to the phantom containing each plug, and comparisons of point-dose measurements and 2-dimensional γ-analysis were performed. RESULTS: For all 26 plugs considered in this study, SPR values for each plug from the SPR images were within 2% agreement with measurements. Additionally, treatment plans developed with the SPR images agreed with the measured point dose to within 2%, whereas a 3% agreement was observed for SECT-based plans. γ-Index pass rates were > 90% for all SECT plans and > 97% for all SPR image-based plans. CONCLUSION: Treatment plans created in a TPS with SPR images obtained from DECT scans are accurate to within guidelines set for validation of clinical treatment plans at our center. The calculated doses from the SPR image-based treatment plans showed better agreement to measured doses than identical plans created with standard SECT scans.

A novel deformable lung phantom with programably variable external and internal correlation.

PURPOSE: Lung motion phantoms used to validate radiotherapy motion management strategies have fairly simplistic designs that do not adequately capture complex phenomena observed in human respiration such as external and internal deformation, variable hysteresis and variable correlation between different parts of the thoracic anatomy. These limitations make reliable evaluation of sophisticated motion management techniques quite challenging. In this work, we present the design and implementation of a programmable, externally and internally deformable lung motion phantom that allows for a reproducible change in external-internal and internal-internal correlation of embedded markers. METHODS: An in-house-designed lung module, made from natural latex foam was inserted inside the outer shell of a commercially available lung phantom (RSD, Long Beach, CA, USA). Radiopaque markers were placed on the external surface and embedded into the lung module. Two independently programmable high-precision linear motion actuators were used to generate primarily anterior-posterior (AP) and primarily superior-inferior (SI) motion in a reproducible fashion in order to enable (a) variable correlation between the displacement of interior volume and the exterior surface, (b) independent changes in the amplitude of the AP and SI motions, and (c) variable hysteresis. The ability of the phantom to produce complex and variable motion accurately and reproducibly was evaluated by programming the two actuators with mathematical and patient-recorded lung tumor motion traces, and recording the trajectories of various markers using kV fluoroscopy. As an example application, the phantom was used to evaluate the performance of lung motion models constructed from kV fluoroscopy and 4DCT images. RESULTS: The phantom exhibited a high degree of reproducibility and marker motion ranges were reproducible to within 0.5 mm. Variable correlation was observed between the displacements of internal-internal and internal-external markers. The SI and AP components of motion of a specific marker had a correlation parameter that varied from -11 to 17. Monitoring a region of interest on the phantom's surface to estimate internal marker motion led to considerably lower uncertainties than when a single point was monitored. CONCLUSIONS: We successfully designed and implemented a programmable, externally and internally deformable lung motion phantom that allows for a reproducible change in external-internal and internal-internal correlation of embedded markers.

Commissioning and acceptance guide for the GammaPod.

The GammaPod breast treatment device has been introduced to provide stereotactic radiation therapy to the breast to patients in the prone position. The GammaPod, using a stereotactic coordinate system, dynamically delivers dose to the target by rotating 25 non-overlapping Co-60 beams while the patient's breast is translated continuously in three axes on the couch during delivery. From simulation to treatment, the patient's breast is immobilized using mild negative pressure (150 mm Hg below atmospheric pressure) through a device-specific dual-cup system with stereotactic fiducials. This technology can be used for boost, multi-fraction partial-breast steterotactic body radiotherapy (SBRT), or single-fraction stereotactic radiosurgery (SRS). This paper reports the commissioning of the system for clinical use. The GammaPod device has four major subsystems: mechanical, dosimetric, radiation safety, and safety interlocks. Detailed methods for testing each subsystem have been identified and quantified. Mechanical systems include couch motion and accuracy along with couch sag. Dosimetric tests include absolute dose calibration, dose profiles, timer error, and plan verifications. Radiation safety includes room and wall surveys, along with device leakage measurements. Safety interlocks deal with power systems, immobilization, and treatment interrupts. The absolute dose rate of the 25 mm collimator was determined using TG-21 dosimetry protocol. The relative output factor for the 15 mm collimator was 0.94. The difference of the full-width-at-half-maximum of the single shot of the 25 mm collimator between the treatment planning system and the measurement was 0.2 mm. All interlocks were found to perform correctly, and the shield was within state and Nuclear Regulatory Commission limits. The items and techniques for commissioning the GammaPod have been developed and tested using the methods reported here.

Inversed-Planned Respiratory Phase Gating in Lung Conformal Radiation Therapy.

PURPOSE: To assess whether the optimal gating window for each beam during lung radiation therapy with respiratory gating will be dependent on a variety of patient-specific factors, such as tumor size and location and the extent of relative tumor and organ motion. METHODS AND MATERIALS: To create optimal gating treatment plans, we started from an optimized clinical plan, created a plan per respiratory phase using the same beam arrangements, and used an inverse planning optimization approach to determine the optimal gating window for each beam and optimal beam weights (ie, monitor units). Two pieces of information were used for optimization: (1) the state of the anatomy at each phase, extracted from 4-dimensional computed tomography scans; and (2) the time spent in each state, estimated from a 2-minute monitoring of the patient's breathing motion. We retrospectively studied 15 lung cancer patients clinically treated by hypofractionated conformal radiation therapy, for whom 45 to 60 Gy was administered over 3 to 15 fractions using 7 to 13 beams. Mean gross tumor volume and respiratory-induced tumor motion were 82.5 cm3 and 1.0 cm, respectively. RESULTS: Although patients spent most of their respiratory cycle in end-exhalation (EE), our optimal gating plans used EE for only 34% of the beams. Using optimal gating, maximum and mean doses to the esophagus, heart, and spinal cord were reduced by an average of 15% to 26%, and the beam-on times were reduced by an average of 23% compared with equivalent single-phase EE gated plans (P<.034, paired 2-tailed t test). CONCLUSIONS: We introduce a personalized respiratory-gating technique in which inverse planning optimization is used to determine patient- and beam-specific gating phases toward enhancing dosimetric quality of radiation therapy treatment plans.

A continuous surface reconstruction method on point cloud captured from a 3D surface photogrammetry system.

PURPOSE: To accurately and efficiently reconstruct a continuous surface from noisy point clouds captured by a surface photogrammetry system (VisionRT). METHODS: The authors have developed a level-set based surface reconstruction method on point clouds captured by a surface photogrammetry system (VisionRT). The proposed method reconstructs an implicit and continuous representation of the underlying patient surface by optimizing a regularized fitting energy, offering extra robustness to noise and missing measurements. By contrast to explicit/discrete meshing-type schemes, their continuous representation is particularly advantageous for subsequent surface registration and motion tracking by eliminating the need for maintaining explicit point correspondences as in discrete models. The authors solve the proposed method with an efficient narrowband evolving scheme. The authors evaluated the proposed method on both phantom and human subject data with two sets of complementary experiments. In the first set of experiment, the authors generated a series of surfaces each with different black patches placed on one chest phantom. The resulting VisionRT measurements from the patched area had different degree of noise and missing levels, since VisionRT has difficulties in detecting dark surfaces. The authors applied the proposed method to point clouds acquired under these different configurations, and quantitatively evaluated reconstructed surfaces by comparing against a high-quality reference surface with respect to root mean squared error (RMSE). In the second set of experiment, the authors applied their method to 100 clinical point clouds acquired from one human subject. In the absence of ground-truth, the authors qualitatively validated reconstructed surfaces by comparing the local geometry, specifically mean curvature distributions, against that of the surface extracted from a high-quality CT obtained from the same patient. RESULTS: On phantom point clouds, their method achieved submillimeter reconstruction RMSE under different configurations, demonstrating quantitatively the faith of the proposed method in preserving local structural properties of the underlying surface in the presence of noise and missing measurements, and its robustness toward variations of such characteristics. On point clouds from the human subject, the proposed method successfully reconstructed all patient surfaces, filling regions where raw point coordinate readings were missing. Within two comparable regions of interest in the chest area, similar mean curvature distributions were acquired from both their reconstructed surface and CT surface, with mean and standard deviation of (µrecon=-2.7×10(-3) mm(-1), σrecon=7.0×10(-3) mm(-1)) and (µCT=-2.5×10(-3) mm(-1), σCT=5.3×10(-3) mm(-1)), respectively. The agreement of local geometry properties between the reconstructed surfaces and the CT surface demonstrated the ability of the proposed method in faithfully representing the underlying patient surface. CONCLUSIONS: The authors have integrated and developed an accurate level-set based continuous surface reconstruction method on point clouds acquired by a 3D surface photogrammetry system. The proposed method has generated a continuous representation of the underlying phantom and patient surfaces with good robustness against noise and missing measurements. It serves as an important first step for further development of motion tracking methods during radiotherapy.