RADIATION MODULATOR AND METHODS OF USE AND PRODUCTION THEREOF
This application claims the benefit of U.S. Provisional Application Ser. No. 62/502,092, filed May 5, 2017, entitled RADIATION MODULATOR AND METHODS OF USE AND PRODUCTION THEREOF, which is hereby incorporated in its entirety herein. The present disclosure generally relates to a brachytherapy applicator, methods of producing the applicator, and methods of treatment using the applicator. High Dose Rate (HDR) brachytherapy, either paired with external beam radiation therapy (EBRT) or delivered alone, is a known treatment modality for cervical cancer at any stage. Like traditional EBRT, the dose delivered to the tumor in a brachytherapy treatment is limited by the presence of surrounding organs at risk (OARs). Existing brachytherapy tandem applicators typically include a single, central lumen and a source that is characterized by a highly isotropic dose profile. These characteristics limit the ability of these existing applicators to satisfy OAR dose constraints while simultaneously delivering the prescribed dose to the tumor. The limitations of these existing applicators are especially apparent in cases where the tumor is large, laterally extended, and/or anisotropically distributed. Existing approaches for treating these extended and/or asymmetric tumors aim to create a more conformal dose profile by supplementing the intracavitary tandem with interstitial brachytherapy. Such approaches include using a modified tandem and ring applicators in which the ring also acts as a template for interstitial needles. Existing approaches further include intensity-modulated brachytherapy (IMBT) that enables anisotropic modulation of the source distribution, dynamic modulated brachytherapy in which the source is encapsulated by a cylindrical shield having a delivery window for radiation, rotating shield brachytherapy that makes uses of an applicator with an electronic brachytherapy source housed in a tandem applicator with an external rotating shield (e.g., rotating window), and direction modulated brachytherapy (DMBT) that makes use of an applicator with a source positioned on the periphery of the applicator as opposed to within a central lumen. Despite the improvements offered by these systems, many of these systems are complicated and may further increase the invasiveness of the treatment. Moreover, these approaches are ultimately still limited by the isotropic dose distribution of the source. Further, none of the existing approaches provide an applicator that is based on a patient's individual anatomy (e.g., the tumor size and position in relation to the position of the surrounding OARs), and is capable of delivering continuous and more conformal radiation dose profiles to extended and/or asymmetric tumors without harming the OARs. Rather, existing approaches utilize standard applicators that deliver doses of radiation that compromise coverage of the tumor (e.g., by covering only portions of the tumor) and/or cover the OARs. Therefore, a need exists for a patient specific intensity-modulated HDR brachytherapy tandem applicator that yields conformal dose distributions, leads to improved target coverage compared to existing brachytherapy treatments, minimizes damage to the OARs, and improves radiation dose delivery time compared to the delivery time of existing brachytherapy treatments. In one aspect, a patient-specific intensity-modulated high dose rate (HDR) brachytherapy applicator for administering an HDR brachytherapy treatment to a patient is provided. The applicator includes a plurality of shielding segments distributed along a central longitudinal axis. Each shielding segment corresponds to one dwell position and includes a shielding wall. Each shielding wall includes a plurality of equiangular shielding sections of varying thickness distributed circumferentially about the central longitudinal axis. Each equilangular shielding section has a shielding thickness. Each shield thickness of each equiangular shielding section at each shielding segment is configured to transmit radiation from an HDR source positioned within each shielding segment into the patient at a predetermined dose rate distribution to administer the HDR brachytherapy treatment. In another aspect, a computer-implemented method for designing a patient-specific intensity-modulated high dose rate (HDR) brachytherapy applicator for administering an HDR brachytherapy treatment to a patient is provided. The applicator includes a plurality of shielding segments distributed along a central longitudinal axis, and each shielding segment includes a plurality of equiangular shielding sections distributed circumferentially about the central longitudinal axis. The method is implemented using at least one processor in communication with at least one memory. The method includes receiving, by a computing device, a radiation treatment plan for administering the HDR brachytherapy treatment. The radiation treatment plan includes a prescribed radiation dosage to be delivered to a region of interest and patient anatomical data representative of the region of interest to be treated. The method also includes determining, by the computing device, an optimal shielding thickness profile and a plurality of optimal dwell times using an inverse planning optimization model constrained by the radiation treatment plan. Each optimal dwell time corresponds to one dwell position, and each dwell position corresponds to one shielding segment. The optimal thickness profile includes a plurality of shield thicknesses, and each shield thickness corresponds to one equiangular shielding section of one shielding segment. The method further includes generating a dwell position-dependent shielding thickness profile that includes the positions of the plurality of the shielding segments and each shield thickness of each equiangular shielding section at each shielding segment. The method additionally includes transmitting, by the computing device, design instructions to a three dimensional (3D) printer for fabrication of the applicator. The design instructions include at least the dwell position-dependent shielding thickness profile. In an additional aspect, a computing device for designing a patient-specific intensity-modulated high dose rate (HDR) brachytherapy applicator for administering an HDR brachytherapy treatment to a patient is provided. The applicator includes a plurality of shielding segments distributed along a central longitudinal axis, and each shielding segment includes a plurality of equiangular shielding sections distributed circumferentially about the central longitudinal axis. The computing device includes at least one processor in communication with at least one memory device, and the at least one processor is programmed to receive a radiation treatment plan for administering the HDR brachytherapy treatment. The radiation treatment plan includes a prescribed radiation dosage to be delivered to a region of interest and patient anatomical data representative of the region of interest to be treated. The at least one processor is also programmed to determine an optimal shielding thickness profile and a plurality of optimal dwell times using an inverse planning optimization model constrained by the radiation treatment plan. Each optimal dwell time corresponds to one dwell position, each dwell position corresponds to one shielding segment, and the optimal thickness profile includes a plurality of shield thicknesses. Each shield thickness corresponds to one equiangular shielding section of one shielding segment. The at least one processor is further programmed to generate a dwell position-dependent shielding thickness profile that includes the positions of the plurality of the shielding segments and each shield thickness of each equiangular shielding section at each shielding segment. The at least one processor is additionally programmed to transmit design instructions to a three dimensional (3D) printer for fabrication of the applicator. The design instructions include at least the dwell position-dependent thickness profile. In another additional aspect, a high-dose radiation (HDR) modulating system configured to improve target coverage of tumor volume during an HDR treatment is provided. The HDR modulating system includes a patient-specific intensity-modulated high dose rate (HDR) brachytherapy applicator that includes a plurality of shielding segments distributed along a central longitudinal axis. Each shielding segment includes a plurality of equiangular shielding sections distributed circumferentially about the central longitudinal axis. The plurality of shielding segments define a central lumen extending along the central longitudinal axis. Each shielding segment further defines a dwell position within the central lumen. The HDR modulating system also includes an HDR source that is movably insertable into the central lumen during an HDR treatment. The HDR source is configured to reside at each dwell position within each shielding segment for a corresponding dwell time. Each corresponding dwell time is based on a radiation therapy plan. Each equiangular shielding section at each shielding segment includes a shield thickness configured to transmit radiation from the HDR source residing at each dwell position at a predetermined dose rate distribution. The patent or application file contains at least one drawing executed in color. Copies of this patent or patent application publication with color drawing(s) will be provided by the Office upon request and payment of the necessary fee. The following drawings illustrate various aspects of the disclosure. The present disclosure relates to radiotherapy systems, devices, and methods for modulating the intensity of x-rays and/or gamma-rays emanating from a radiation source utilized to treat cancerous tumors. Such techniques can enable treatments that provide a non-invasive alternative to existing brachytherapy (CBT) treatments to treat cancerous tumors, such as cervical tumors. In various aspects, mathematical modeling is used for the optimization of a 3D printed patient-specific intensity-modulated brachytherapy (IMBT) applicator configured for use in high dose rate (HDR) treatments to successfully modulate radiation intensity to deliver focused radiation to a pathology site (e.g., gross tumor volume) while minimizing radiation exposure to surrounding organs at risk (OARs). As described herein, the present disclosure is directed to patient-specific IMBT applicators for use in HDR treatments. In one aspect, the external shape of the applicator resembles an existing brachytherapy (BT) tandem applicator, as shown in In contrast, existing BT tandem applicators consist of a single central lumen with a uniform cross-sectional profile, as shown in In various aspects, a patient's anatomical information includes information in regards to the patient's tumor. The patient's anatomical information can be provided by imaging modalities typically used for radiation therapy planning, such as, but not limited to, computerized tomography (CT) scans, ultrasonography scans, and magnetic resonance imaging (MRI) scans of the tumor and surrounding OARs. More specifically, the patient's anatomical information may include information defining the gross tumor volume (GTV), defined herein as the known tumor volume that can be seen, measured, and/or palpated. The patient's anatomical information may also include information defining the clinical target volume (CTV), defined herein as the volume of suspected tumor infiltration surrounding the GTV. The CTV can include the volume of suspected tumor extensions that may or may not be fully imaged and/or accurately defined. In various aspects, the radiation dosage distributions may be generated by HDR sources such as isotopes including, but not limited to,192Ir,131Cs,125I,103Pd,198Au,187W,169Yb,145Sm,137Cs,109Cd,65Zn,153Gd,57Co,56Co, and58Co. In another aspect, radiation dosage distributions can be generated by an HDR source, such as an electronic brachytherapy (eBT) source contained within a novel modulator comprising of high-Z material (e.g., an atomic number “Z” that is greater than or equal to 22). Such isotopes can be referred to as, for example, non-electronic brachytherapy (BT) sources. In one aspect, an individualized applicator, specifically the interior of the individualized applicator, is designed by optimizing wall thickness of the applicator and the dwell time based on a patient's anatomy and the prescribed treatment dosage. In this aspect, the interior of the individualized applicator is sub-divided into a plurality of equiangular shielding sections. Each section of the plurality of equiangular shielding sections may independently vary in shielding thickness at each dwell position of the applicator to enable anisotropic modulation by defining multiple emission windows at each dwell position, similar to paddle-based rotating shield brachytherapy (P-RSBT), shown illustrated in In one aspect, the thickness of the inner tandem applicator varies such that maximum radiation is delivered to tumor regions and minimal and/or zero radiation is delivered to nearby OARs. In another aspect, the disclosed applicator is optimized to modulate radiation intensity such that focused radiation is delivered to the GTV and the CTV. In various aspects, the number of equiangular shielding sections provided at each dwell position of the individualized applicator may depend upon a specific patient's anatomical information. By way of non-limiting example, the number of equiangular shielding sections can be increased to create a more spatially tailored dose distribution for a patient having a tumor that is laterally extended and/or anisotropically distributed. Increasing the number of equiangular shielding sections at each dwell position increases the computational cost of design optimization by the disclosed computer-implemented methods. However, the increased computation time is less than the computation time typically used for existing BT approaches. In another non-limiting example, the number of sections can be reduced for a patient depending on tumor position and/or volume with respect to OARs. Restricting the number of equiangular shielding sections at each dwell position reduces the computational cost of design optimization by the disclosed computer-implemented methods. In another additional aspect, the number of equiangular shielding sections at each dwell position of the individualized applicator may be selected based on the resolution of the patient's anatomical information, as well as the estimated degree of movement of anatomical features of the patient during acquisition of the patient's anatomical information and/or during treatment. Without being limited to any particular theory, internal organs, tissues, and other anatomical landmarks are subject to a limited degree of movement. Consequently, the location of the tumor (i.e. GTV and CTV) and OARs at any given time in each patient may be subject to some degree of uncertainty. If a high number of equiangular shielding sections are included at each dwell position of the individualized applicator, the spatial resolution of the resulting radiotherapy dosage maps may be needlessly high given the lower degree of precision at which the patient's anatomical data is known. In one aspect, the inner tandem wall of the applicator is divided into six equiangular shielding sections at each dwell position of the HDR source, as shown in Optimization Model for Designing Patient-SpecificIntensity-Modulated HDR BT Applicator In various aspects, an optimization model is used to determine the thicknesses of the inner tandem wall (e.g., the shielding wall) of the applicator for each equiangular segment of each dwell position such that focused radiation (e.g., maximum radiation at the prescribed dosage) is delivered to the GTV and CTV, and minimum or no radiation is delivered to surrounding OARs. Applicator design parameters determined by the optimization model including, but not limited to, the inner wall thicknesses at each equiangular segment of each dwell position are input into 3D modeling software for 3D printing of the patient-specific intensity-modulated HDR BT tandem applicator. In various aspects, the 3D printed applicator is formed form a tungsten material (shown in In one aspect, to compute the radiation dose rate, a 1D isotropic point source dose rate calculation formulation suggested by AAPM Report TG-43 is utilized. Specifically, the dose rate {dot over (D)}() at a voxel at position from a point source at position is calculated as expressed in Equation 1: where SKdenotes the air-kerma strength (units μGym2h−1), Λ denotes the dose rate constant in water, gP( ) denotes the radial dose function of the point source, and ϕan( ) denotes the 1D anisotropy function. For the purpose of calculation, the volume of interest is discretized into voxels with resolution [rxmm×rymm×rzmm] and index i. The dose rate in voxel i is subsequently induced by the jthsource, denoted by {dot over (D)}iJ. Assuming there are Nssource positions located along the tandem separated by an equal distance dsand the tandem consists of Ntpieces of shielding, the total dose received by the ithvoxel is given by In Equation 2, Tg(i,j)α(i,j)is the transmission rate of a given piece of shield indexed by g(i,j) ranging from 1 to Nt, which is a function that determines which shield a ray originating at the jthsource will cross in travelling to the ithvoxel. In Equation 2, α(i, j) is a constant representing the ratio which reflects the effect of the thickness of the shield on the ray passing through it, as shown in In one aspect, to determine the dwell time for each dwell position and the transmission rate for each shielding portion, the following optimization model is utilized. In this optimization model, Npdwell positions are assumed and t∈ represents the vector of dwell times, and T∈ represents the vector of transmission rates, where ={t∈N In this optimization model, l, u∈(0,1) are the lower and upper bounds of the transmission rate, which correspond to the thickness of each shielding portion. Considering the design limitations of the shield, l=0.38 and u=0.90 were set corresponding to thicknesses of 4.61 mm and 0.50 mm respectively. Subsequently, ={1,2, . . . , No} represents the set of indexed OARs to be considered in the optimization method. o⊂ represents the set of voxels making up the o-th OAR, while t⊂ represents the set of voxels in the HR-CTV. The cost function for this optimization model is defined according to Equation 3 below: The first term in Equation 3 corresponds to the costs for OARs with the given configuration of dwell times t and the transmission rates T while the second term corresponds to the cost for the tumor. Additionally, the cost function at voxel level for the OARs is defined according to Equation 4 as follows: The cost function at voxel level for the tumor is defined according to Equation 5 as follows: where Di(t,T) is the dose calculated using Equation 2, and Dt∈Nrepresents the target dose for each voxel, which is defined as Here, Dois the maximum dose that could be delivered to the o-th OAR while Dtis the minimum dose that should be delivered to the tumor. Additionally, in Equations 4 and 5, C∈ is a constant that scales down the cost, while Soand Stare constants that control the relative importance for the o-th OAR or the tumor respectively. Unlike the general multi-objective optimization approach, in which a set of weighting parameters are utilized in presenting the total cost function as a convex combination of all the individual terms, horizontal shifting constants Soand Stare used to balance the relative importance of the individual cost functions. These cost functions are illustrated in The cost function presented in Equation 3 above may be expressed in a denser form to provide the optimization model given by Equation 7 below: Equation 7 includes two blocks of variables and is convex for dwell time t when fixing transmission rates T, but not for the alternate case for transmission rates T when fixing dwell time t. To show this, it can be verified that fio(t,T) is convex for one block of variables when fixing the other and fit(t,T) is convex for t but concave for T. Specifically, D, is a linear function of t so it is both convex and concave, but it is convex for T as the Hessian is positive semidefinite. On the other hand, fiois a convex non-decreasing function for t and T while fitis convex and non-increasing. Since the summation of convex functions is convex, it follows that F is convex for t but the convexity with respect to T is not clear. In one aspect, to solve Equation 7 an alternating minimization scheme is utilized to search for t and T in turns. This alternating minimization algorithm starts at an initial point (t0, T0)∈× and solves the two subproblems with respect to t and T while fixing the other by gradient descent with back-tracking line search. Specifically, the partial gradients of fi(t,T) with respect to the jthdwell time is expressed as Equation 9 as follows: and partial gradients of fi(t,T) for the kthtransmission rate Tkis expressed as Equation 10 below: The gradient of F(t,T) with respect to t and T can be obtained by taking summations of Equations 9 and 10. Without being limited to any particular theory, the gradient descent method with constant step size bounded by is unlikely to converge to a strict saddle point, where L is the Lipschitz constant for a function ƒ∈C2. This conclusion does not necessarily hold for a gradient descent method with line search. In some aspects, a random perturbation to the iterate Timay be made when the step size is sufficiently small, or approximately less than the Lipschitz constant, and then the constant step size may be used to continue the optimization. The scheme of the algorithm in this aspect is summarized as shown below in Table 1. In one aspect, for user-specified set of constants Soand St, Algorithm 1 will generate a plan with dwell time t and transmission rates T. However, the estimated dose volume histogram associated with this generated plan may not satisfy the clinical goal for tumor coverage and OAR dose. In another aspect, an automatic mechanism to tune the control constants Soand Stmay be introduced to generate a satisfactory plan. In this other aspect, the algorithm assumes a relatively large Stinitially to ensure that the OAR sparing fails for at least some OARs, indicating that the current weighting of the cost function favors tumor coverage. Subsequently, the algorithm gradually increases Soif the o-th OAR received excessive dose, and perform an additional iteration with the updated So. This procedure terminates when the OAR doses satisfy the prescribed criteria. This process of tuning the control constants Soand Stis summarized as Algorithm 2, shown below in Table 2. In various additional aspects, an initial guess of Stmay be used to tune Soand Stautomatically in a manner similar to Algorithm 2. In various aspects, the optimal thickness of the shielding wall and the dwell time (e.g., the amount of time at which the HDR source delivers radiotherapy at each dwell position) are determined using the optimization model as described above. The optimal thickness and dwell time at each dwell position are determined according to the specific patient's anatomical information and the dosage prescribed by the patient's physician. As described above, the optimization model considers the transmission rates of the shielding wall of the tandem applicator, which depend on shielding thickness, and the dwell time of the HDR source as variables to be calculated in order to achieve the best possible target coverage. In various aspects, the best possible target coverage is where maximum radiation is delivered to the GTV and the CTV while minimum to no radiation is delivered to the OARs. By using alternating minimization to solve the optimization model described above, optimal thickness of the shielding wall and dwell time at each dwell position are calculated. In various aspects, these calculated parameters are subsequently input into 3D modeling software and used to 3D print an individualized applicator that, on the exterior, appears similar to existing applicators (e.g., cylindrical shape), but on the interior, differs in wall thickness based on the individual patient's anatomy and treatment needs. The determined thickness of the shielding wall at each dwell position is based on factors including, but not limited to, the thickness of the applicator inserted in the patient, the patient's anatomy, the patient's tolerance level, and whether or not certain degrees of thickness are suitable for continuous delivery of radiation. In various aspects, tungsten material, as shown in In some aspects, the above described methods and processes may be implemented using a computing system, including one or more computing devices. The methods and processes described herein may be implemented as computer applications, computer services, computer APIs, computer libraries, and/or any other computer program product without limitation. In some aspects, the computing device 1002 may be communicably coupled to at least one of a medical scanner 1010, a patient records server 1012, and a 3D printing device 1014 via a network 1016. In various aspects, the medical scanner 1010 can be any medical imaging system configured to obtain suitable images of a tumor (i.e. GTV and CTV) and organs (i.e. OARs) for analysis by the computing device 1002 including, but not limited to, an MRI scanner, a CT scanner, and any other suitable medical imaging device without limitation. In various aspects, the computing device 1002 receives imaging data from medical scanner 1010 as well as a radiation treatment plan from the database 1008, and applies an optimization model to the imaging data and the radiation treatment plan to determine optimal thickness parameters for the applicator and dwell times at each dwell position, and transmits instructions to the 3D printing device 1014 for automatically creating a patient-specific applicator. The network 1016 may be any network that allows local area or wide area communication between the devices. For example, the network 1016 may allow communicative coupling to the Internet through at least one of many interfaces including, but not limited to, at least one of a network, such as the Internet, a local area network (LAN), a wide area network (WAN), an integrated services digital network (ISDN), a dial-up-connection, a digital subscriber line (DSL), a cellular phone connection, and a cable modem. The processor 1105 may be operatively coupled to a communication interface 1115 such that the server system 1102 may be capable of communicating with a remote device such as a medical imaging device 1010, a patient records server 1012, a 3D printing device 1014, a treatment device 1018 (all shown in The processor 1105 may also be operatively coupled to a storage device 1125. The storage device 1125 may be any computer-operated hardware suitable for storing and/or retrieving data. In some aspects, the storage device 1125 may be integrated within the server system 1102. By way of non-limiting example, the server system 1102 may include storage device 1125 in the form of one or more hard disk drives. In other aspects, the storage device 1125 may be external to the server system 1102 and may be accessed by a plurality of server systems in addition to the server system 1102. By way of non-limiting example, the storage device 1125 may include multiple storage units including, but not limited to, hard disks or solid state disks in a redundant array of inexpensive disks (RAID) configuration. Other non-limiting examples of suitable storage devices 1125 include storage area networks (SAN) and/or network attached storage (NAS) systems. In some aspects, the processor 1105 may be operatively coupled to the storage device 1125 via a storage interface 1120. The storage interface 1120 may be any component capable of providing the processor 1105 with access to the storage device 1125. Non-limiting examples of suitable storage interfaces 1120 include Advanced Technology Attachment (ATA) adapters, Serial ATA (SATA) adapters, Small Computer System Interface (SCSI) adapters, RAID controllers, SAN adapters, network adapters, and/or any suitable components providing the processor 1105 with access to the storage device 1125. The memory 1110 may include, but is not limited to, random access memory (RAM) such as dynamic RAM (DRAM) or static RAM (SRAM), read-only memory (ROM), erasable programmable read-only memory (EPROM), electrically erasable programmable read-only memory (EEPROM), and non-volatile RAM (NVRAM). The above memory types are example only, and are thus not limiting as to the types of memory suitable for storage of a computer program. Referring again to The applicator design parameters 1226 may include a plurality of parameters defining the geometry of the patient-specific applicator as calculated by the optimization model component 1260 using the optimization model 1222, including, but not limited to, a position of each dwell position along the applicator and thicknesses for each angular segment of each dwell position. The 3D printing instructions may include instructions generated by the fabrication component 1250 and based on the applicator design parameters 1226 to be used by a 3D printing device 1014 to fabricate the patient-specific applicator. The treatment device instructions 1232 may include instructions used by a treatment component 1270 to operate a treatment device 1018 to administer a treatment using the patient-specific applicator. The computing device 1202 also includes at least several components configured to perform specific tasks associated with designing and producing a patient-specific applicator and to administer a treatment using the patient-specific applicator as disclosed herein. In one aspect, the computing device 1202 includes a data storage device 1240, an optimization model component 1250, a fabrication component 1260, and a treatment component 1270. The data storage device 1240 is configured to store data received or generated by computing device 1202, such as any of the data stored in database 1220 or any outputs of processes implemented by any component of computing device 1202. The optimization model component 1250 is configured to receive a radiation treatment plan 1224 from the database 1220 for treating a region of interest. More specifically, the fabrication component 1260 is configured to receive a radiation treatment plan 1224 that includes, but is not limited to, a prescribed radiation dosage to be delivered to the region of interest and patient anatomical data of the region of interest to be treated for use in the design of the patient-specific applicator using the optimization model 1222. The optimization model component 1260 is configured to apply an inverse planning optimization model using the received prescribed radiation dosage and the received patient anatomical data (the radiation treatment plan 1224 and the imaging data 1230) to determine an optimal thickness of the interior shielding within the patient-specific applicator at a plurality of dwell positions within the region of interest. In addition, the optimization model component 1250 is configured to transmit the calculated plurality of dwell positions and associated shield thickness profiles for each dwell position to the fabrication component 1260. The optimization model component 1250 is further configured to transmit the plurality of dwell positions and associated dwell times to the treatment component 1270. The fabrication component 1260 is configured to generate instructions used to operate the 3D printing device 1014 for fabrication of the patient-specific applicator. More specifically, the fabrication component 1260 is configured to receive design instructions from the optimization model component 1250 that include at least the plurality of dwell position and the associated shield thickness profiles. The treatment component 1270 is configured to generate instructions that may be used to operate a treatment device 1018 to administer a treatment to a patient using the patient-specific applicator produced by the fabrication component 1260 using the 3D printing device 1014. In one aspect, the treatment component 1270 is configured to receive a plurality of dwell positions and associated schedule of dwell times for the patient-specific applicator from the optimization model component. The treatment fabrication component 1270 is further configured to modify each dwell time from the schedule of dwell times based on the age/condition of the radiation source using any method known in the art. In one aspect, the treatment component 1270 may control the operation of a treatment device 1018 including, but not limited to, a radiotherapy device, to administer a radiotherapy treatment to the patient using the patient-specific applicator. The computer systems and computer-implemented methods discussed herein may include additional, less, or alternate actions and/or functionalities, including those discussed elsewhere herein. The computer systems may include or be implemented via computer-executable instructions stored on non-transitory computer-readable media. The methods may be implemented via one or more local or remote processors, transceivers, servers, and/or sensors (such as processors, transceivers, servers, and/or sensors mounted on vehicle or mobile devices, or associated with smart infrastructure or remote servers), and/or via computer executable instructions stored on non-transitory computer-readable media or medium. In one embodiment, a computer program is provided, and the program is embodied on a computer-readable medium. In an example embodiment, the system is executed on a single computer system, without requiring a connection to a server computer. In a further example embodiment, the system is being run in a Windows® environment (Windows is a registered trademark of Microsoft Corporation, Redmond, Wash.). In yet another embodiment, the system is run on a mainframe environment and a UNIX® server environment (UNIX is a registered trademark of X/Open Company Limited located in Reading, Berkshire, United Kingdom). In a further embodiment, the system is run on an iOS® environment (iOS is a registered trademark of Cisco Systems, Inc. located in San Jose, Calif.). In yet a further embodiment, the system is run on a Mac OS® environment (Mac OS is a registered trademark of Apple Inc. located in Cupertino, Calif.). In still yet a further embodiment, the system is run on Android® OS (Android is a registered trademark of Google, Inc. of Mountain View, Calif.). In another embodiment, the system is run on Linux® OS (Linux is a registered trademark of Linus Torvalds of Boston, Mass.). The application is flexible and designed to run in various different environments without compromising any major functionality. The following detailed description illustrates embodiments of the disclosure by way of example and not by way of limitation. It is contemplated that the disclosure has general application to providing an on-demand ecosystem in industrial, commercial, and residential applications. The methods and systems described herein can be used to treat any disease, disorder, or condition that can be treated with traditional brachytherapy. For example, diseases, disorders, and/or conditions can include pathology, tumor, and cancer such as, but not limited to, prostate cancer, breast cancer, lung cancer, esophageal cancer, gynecologic cancer, anal/rectal tumor, sarcoma; and head or neck cancer. The tumor can be cancerous or non-cancerous lesions. The present devices and methods enable treatment of lateral tumor extensions by delivering targeted radiation dosages to these areas. Lateral tumor extensions are difficult to treat with existing intracavitary brachytherapy (BT) due to radiation dose limitations imposed by the presence of nearby healthy tissues and organs, such as the bladder, rectum, and/or sigmoid in the case of cervical cancer treatment. In another aspect, the present devices and methods can enable increased dosage conformity for non-symmetric tumors by utilizing a device that can shield radiation emanated from an electronic brachytherapy (eBT) source or non-electronic brachytherapy (BT) source. In one example, the device includes a radiation modulator that includes a material having a position-dependent thickness that is based at least on (i) a radiation therapy plan specific to a patient and (ii) a geometry of a patient region to be treated (e.g., tumor region). In an additional or alternative aspect, the device includes an HDR source that is movably inserted into an enclosure coupled to the radiation modulator. In various aspects, the methods as described herein can include, the HDR source residing at a plurality of locations within the radiation modulator during a respective plurality of dwell times based on a patient's radiation therapy plan. Also provided is a process of treating a pathology, cancer, or tumor in a subject in need administration of a therapeutically effective amount of radiation, so as to destroy pathologic cells and shrink tumors. Methods described herein are generally performed on a subject in need thereof. A subject in need of the therapeutic methods described herein can be a subject having, diagnosed with, suspected of having, or at risk for developing pathologic cells, tumors, or cancer. A determination of the need for treatment will typically be assessed by a history and physical exam consistent with the disease or condition at issue. Diagnosis of the various conditions treatable by the methods described herein is within the skill of the art. The subject can be an animal subject, including a mammal, such as horses, cows, dogs, cats, sheep, pigs, mice, rats, monkeys, hamsters, guinea pigs, and chickens, and humans. For example, the subject can be a human subject. Generally, a safe and effective amount of radiation is, for example, that amount that would cause the desired therapeutic effect in a subject while minimizing undesired side effects. In various embodiments, an effective amount of radiation described herein can substantially inhibit tumor or pathologic cell growth, slow the progress of tumor or pathologic cell growth, or limit the development of tumor or pathologic cell growth. When used in the treatments described herein, a therapeutically effective amount of radiation can be any amount as prescribed by a radiologist. Again, each of the states, diseases, disorders, and conditions, described herein, as well as others, can benefit from the treatment methods described herein. Generally, treating a state, disease, disorder, or condition includes preventing or delaying the appearance of clinical symptoms in a mammal that may be afflicted with or predisposed to the state, disease, disorder, or condition but does not yet experience or display clinical or subclinical symptoms thereof. Treating can also include inhibiting the state, disease, disorder, or condition, e.g., arresting or reducing the development of the disease or at least one clinical or subclinical symptom thereof. Furthermore, treating can include relieving the disease, e.g., causing regression of the state, disease, disorder, or condition or at least one of its clinical or subclinical symptoms. A benefit to a subject to be treated can be either statistically significant or at least perceptible to the subject or to a physician. Administration of radiation can occur as a single event or over a time course of treatment. For example, radiation can be administered daily, weekly, bi-weekly, or monthly. For treatment of acute conditions, the time course of treatment will usually be at least several days. Certain conditions could extend treatment from several days to several weeks. For example, treatment could extend over one week, two weeks, or three weeks. For more chronic conditions, treatment could extend from several weeks to several months or even a year or more. As another example, a radiation delivery device can be implanted. Treatment in accord with the methods described herein can be performed prior to, concurrent with, or after existing treatment modalities for tumor or pathologic cell (e.g., cancer) growth. Radiation treatment, as described herein, can be administered according to methods described herein and in a variety of means known to the art (see e.g., U.S. Patent Application Publication No. 2014/0249406; U.S. Patent Application Publication No. 2015/0367144; and U.S. Patent Application Publication No. 2016/0271379, incorporated by reference in their entireties herein). As discussed above, radiation therapy can be administered in a dose or a plurality of doses or the radiation can be delivered via an implant. Also provided are kits. Such kits can include an agent or composition described herein and, in certain embodiments, instructions for administration. Such kits can facilitate performance of the methods described herein. When supplied as a kit, the different components of the composition can be packaged in separate containers and admixed immediately before use. Components include, but are not limited to software, 3D printing materials, or a 3D printer. Such packaging of the components separately can, if desired, be presented in a pack or dispenser device which may contain one or more unit dosage forms containing the composition. The pack may, for example, comprise metal or plastic foil such as a blister pack. In certain embodiments, kits can be supplied with instructional materials. Instructions may be printed on paper or other substrate, and/or may be supplied as an electronic-readable medium, such as a floppy disc, mini-CD-ROM, CD-ROM, DVD-ROM, Zip disc, videotape, audio tape, and the like. Detailed instructions may not be physically associated with the kit; instead, a user may be directed to an Internet web site specified by the manufacturer or distributor of the kit. Definitions and methods described herein are provided to better define the present disclosure and to guide those of ordinary skill in the art in the practice of the present disclosure. Unless otherwise noted, terms are to be understood according to conventional usage by those of ordinary skill in the relevant art. As employed in this specification and annexed drawings, the terms “unit,” “component,” “interface,” “system,” “platform,” “stage,” and the like are intended to include a computer-related entity or an entity related to an operational apparatus with one or more specific functionalities, wherein the computer-related entity or the entity related to the operational apparatus can be either hardware, a combination of hardware and software, software, or software in execution. One or more of such entities are also referred to as “functional elements.” As an example, a unit may be, but is not limited to being, a process running on a processor, a processor, an object, an executable computer program, a thread of execution, a program, a memory (e.g., a hard disc drive), and/or a computer. As another example, a unit can be an apparatus with specific functionality provided by mechanical parts operated by electric or electronic circuitry which is operated by a software or a firmware application executed by a processor, wherein the processor can be internal or external to the apparatus and executes at least a part of the software or firmware application. In addition or in the alternative, a unit can provide specific functionality based on physical structure or specific arrangement of hardware elements. As yet another example, a unit can be an apparatus that provides specific functionality through electronic functional elements without mechanical parts, the electronic functional elements can include a processor therein to execute software or firmware that provides at least in part the functionality of the electronic functional elements. An illustration of such apparatus can be control circuitry, such as a programmable logic controller. The foregoing example and related illustrations are but a few examples and are not intended to be limiting. Moreover, while such illustrations are presented for a unit, the foregoing examples also apply to a component, a system, a platform, and the like. It is noted that in certain embodiments, or in connection with certain aspects or features thereof, the terms “unit,” “component,” “system,” “interface,” “platform” can be utilized interchangeably. In some embodiments, numbers expressing quantities of ingredients, properties such as molecular weight, reaction conditions, and so forth, used to describe and claim certain embodiments of the present disclosure are to be understood as being modified in some instances by the term “about.” In some embodiments, the term “about” is used to indicate that a value includes the standard deviation of the mean for the device or method being employed to determine the value. In some embodiments, the numerical parameters set forth in the written description and attached claims are approximations that can vary depending upon the desired properties sought to be obtained by a particular embodiment. In some embodiments, the numerical parameters should be construed in light of the number of reported significant digits and by applying ordinary rounding techniques. Notwithstanding that the numerical ranges and parameters setting forth the broad scope of some embodiments of the present disclosure are approximations, the numerical values set forth in the specific examples are reported as precisely as practicable. The numerical values presented in some embodiments of the present disclosure may contain certain errors necessarily resulting from the standard deviation found in their respective testing measurements. The recitation of ranges of values herein is merely intended to serve as a shorthand method of referring individually to each separate value falling within the range. Unless otherwise indicated herein, each individual value is incorporated into the specification as if it were individually recited herein. In some embodiments, the terms “a” and “an” and “the” and similar references used in the context of describing a particular embodiment (especially in the context of certain of the following claims) can be construed to cover both the singular and the plural, unless specifically noted otherwise. In some embodiments, the term “or” as used herein, including the claims, is used to mean “and/or” unless explicitly indicated to refer to alternatives only or the alternatives are mutually exclusive. The terms “comprise,” “have” and “include” are open-ended linking verbs. Any forms or tenses of one or more of these verbs, such as “comprises,” “comprising,” “has,” “having,” “includes” and “including,” are also open-ended. For example, any method that “comprises,” “has” or “includes” one or more steps is not limited to possessing only those one or more steps and can also cover other unlisted steps. Similarly, any composition or device that “comprises,” “has” or “includes” one or more features is not limited to possessing only those one or more features and can cover other unlisted features. All methods described herein can be performed in any suitable order unless otherwise indicated herein or otherwise clearly contradicted by context. The use of any and all examples, or exemplary language (e.g. “such as”) provided with respect to certain embodiments herein is intended merely to better illuminate the present disclosure and does not pose a limitation on the scope of the present disclosure otherwise claimed. No language in the specification should be construed as indicating any non-claimed element essential to the practice of the present disclosure. Groupings of alternative elements or embodiments of the present disclosure disclosed herein are not to be construed as limitations. Each group member can be referred to and claimed individually or in any combination with other members of the group or other elements found herein. One or more members of a group can be included in, or deleted from, a group for reasons of convenience or patentability. When any such inclusion or deletion occurs, the specification is herein deemed to contain the group as modified thus fulfilling the written description of all Markush groups used in the appended claims. All publications, patents, patent applications, and other references cited in this application are incorporated herein by reference in their entirety for all purposes to the same extent as if each individual publication, patent, patent application or other reference was specifically and individually indicated to be incorporated by reference in its entirety for all purposes. Citation of a reference herein shall not be construed as an admission that such is prior art to the present disclosure. Having described the present disclosure in detail, it will be apparent that modifications, variations, and equivalent embodiments are possible without departing the scope of the present disclosure defined in the appended claims. Furthermore, it should be appreciated that all examples in the present disclosure are provided as non-limiting examples. The following non-limiting examples are provided to further illustrate the present disclosure. It should be appreciated by those of skill in the art that the techniques disclosed in the examples that follow represent approaches that may function well in the practice of the present disclosure, and thus can be considered to constitute examples of modes for its practice. However, those of skill in the art should, in light of the present disclosure, appreciate that many changes can be made in the specific embodiments that are disclosed and still obtain a like or similar result without departing from the spirit and scope of the present disclosure. To validate the methods described above, the following experiments were conducted. A tandem applicator was designed such that the external shape of the patient-specific tandem applicator resembles the existing tandem applicator's external shape (e.g., a cylindrical shape). The wall thickness inside the patient-specific tandem applicator and the dwell time at each dwell position p were simultaneously optimized to provide varying degrees of thickness around each circumference at each dwell position based on the specific patient's anatomy and the prescribed radiation dosage. Tungsten material was used to 3D print the patient-specific tandem applicator. To generate a model for 3D printing, an optimization model as described above was solved to generate the design parameters of the patient-specific tandem applicator. More specifically, the optimization model was solved to generate the transmission rate relating to tungsten thickness for a given schedule of dwell times. The following cost function was used to calculate the optimization model parameters: where Xr(i)=ΣpKtpDp(i)Tg(i,p)α(i,g(i,p))−Do,r(i), D=dose, p=dwell position, t is the dwell time, and T=transmission factor at each specific angular section. After optimizing the cost function, the model parameters were inputted into 3D modeling software (Cinema 4D R17, Maxon) to generate a stereolithographic (STL) file for use by the 3D printing device. At each dwell position of the HDR source, the surrounding tandem wall was divided into sections or segments with varying thicknesses. The thickness of each section varied from 0.12 cm to 0.48 cm inwards, such that the exterior surface of the tandem applicator resembled the traditional applicator (diameter=1.2 cm). The inner wall thickness profiles were generated using the ‘inner extrude’ function in the 3D modeling software for all pre-defined dwell positions. The resulting 3D model was then exported to the 3D printing software (Simplify3D) and converted from 3D volumetric data into an STL file. A 3D metal printer (ProX DMP 320) was utilized to fabricate the tandem applicator using tungsten, as shown in 3D printing the patient-specific tandem applicator using design parameters calculated from the optimization model described above was achieved without any major discrepancies between the digital and physical models. A comparison between the model parameters and measurements from the 3D printed model of the applicator indicated an accuracy within about 0.1 mm. The disclosed method of utilizing an inverse planning optimization model to design a patient-specific tandem applicator was demonstrated to be a feasible alternative to existing tandem applicators, especially when surrounding OARS significantly constrain the tumor-dose coverage during HDR. Further, the disclosed method may be adapted to other HDR sites such as rectal (shown in To validate the method described above, the following experiments were conducted. Numerical experiments were performed on a 2D phantom, shown in Numerical experiments were also performed using 2D patient data from a clinically-treated cervical cancer patient for further validation. Numerical experiments were also performed using 3D patient data from a clinically treated cervical cancer patient for further validation. The following dose rate calculation formulation, based on the formulation described in the AAPM TG-43 Report, the contents of which are incorporated by reference herein in its entirety, was used to calculate the dose rate for the region of interest, as expressed in Equation 12: The dose at one dwell position was subsequently calculated using the following: In Equation 13, tiis the HDR source dwell time for each dwell position in seconds, Tiis the transmission factor (e.g., the transmission rate) of the HDR source, and 0.25 (Maximum Thickness)≤Ti≤0.8 (Minimum thickness). The transmission factor Tiof the HDR source is related to the thickness of the shield wall of the tandem applicator. IMRT optimization was performed using Equation 14 below: In Equation 14, Dcis the calculated dose of a 2D or 3D matrix depending on the dwell time tiand the transmission rate Ti. The exponential cost function of Equation 15 was chosen for this optimization: An objective function was utilized as expressed in Equations 16 and For the objective functions, tpwas dwell time at location p, Tjwas the transmission rate at index j, and both tpand Tjwere varied during optimization. Dpis the dose rate matrix at dwell location p. Additionally, c and srare constants that control the shape of the cost function. Further, g(i,p) is the index of the transmission rate, which depends on the dwell location p and the pixel location, and α(i,g(i,p)) is a constant that depends on the pixel location and the index of transmission rate. Both the dwell time and the thickness of the equiangular shielding sections of the tandem applicator were optimized using an alternating minimization scheme to search for tpand Tj. More specifically, alternative minimization with gradient descent and back-tracking line search was utilized using the method as follows: STEP 1. Update the dwell time: STEP 2. Update the transmission rate: STEP 3. Check convergence, if there is no convergence, go to STEP 1; otherwise go to STEP 4. STEP 4. Check the constraints as expressed in Equation 20 below. If the constraints are not satisfied, modify the model by changing srand proceed to STEP 1. If the constraints are satisfied, terminate. Equation 20 provided partial derivatives with respect to the dwell time tiand the transmission rate Ti. In Equation 20, I(p,i,j) was an indicator function that defined whether, at dwell location p, the ithpixel would be affected by the jthtransmission rate. The cost function was then alternately minimized with respect to one variable (while fixing the other), using the following: where St={t|t} and ST={T|0.25T1}. For each sub-problem, the gradient descent with back-tracking line search was used. Upon convergence, the constraints for the OARs were checked. If the constraints were not met, the cost function was automatically modified by changing srsuch that the OARs were favored. The compare the disclosed patient-specific tandem applicator method with the existing HDR method, the same optimization model, as described above, was used with the exception of fixing the transmission rates as a constant (T=1.0). For both the experiments performed on the 2D phantom and the 2D patient data, the dose constraints outlined in Table 3 below were assigned to the CTV, bladder, rectum, and sigmoid. More specifically, Table 3 provided dose constraints as prescribed by the patient's physician. For the experiments performed on the 3D patient data, the dose constraint for the CTV was 560 cGy. The dose constraints for the bladder, rectum, and sigmoid were the same as those outlined in Table 3 for the 2D phantom and the 2D patient data. For the patient data experiments, clinically treated cervical HDR patient data were used. A radioactive Ir-192 source was utilized, and data was collected for the twelve dwell positions of the Ir-192 source that were monitored. A configuration was calculated for an existing HDR method, as shown in The advantages offered by the disclosed patient-specific tandem applicator design were demonstrated in the most realistic 2D patient case as well. A configuration was calculated using the 2D patient model for an existing HDR method, as shown in The comparison of treatments administered using the existing HDR method and the patient-specific tandem applicator method was repeated using 3D patient data. The 3D patient data included dimensions of 332×502×118 cm3and had an image resolution of 0.29×0.29×0.9 cm3. Experiments on the 3D patient data were implemented using CUDA C++ to enable parallel computation. The computation time was under one minute using a system with the following specifications: a CPU of Intel i7-6700K 4.00 GHz, a GPU of NVidia GTX 1080, and a memory of 32 GB DDR4 3200 MHz. The patient-specific tandem applicator design yielded benefits when applied to the 3D patient data. The results of these experiments validated the patient-specific tandem applicator design and treatment method. The patient-specific tandem applicator design improved dose coverage of the CTV by 9-44% without compromising the surrounding OARs. The patient-specific tandem applicator design yielded benefits when applied to the 2D patient case by covering 99.92% of the CTV in comparison to the existing HDR method only covering 56.21% of the CTV. For both the 2D phantom case and the 2D patient case, the complete tumor coverage was achieved while simultaneously satisfying the OAR constraints. The patient-specific tandem applicator method significantly improved the coverage by approximately 70% in the 2D phantom case and 78% in the 2D patient case. The patient-specific tandem applicator design also yielded benefits when applied to 3D patient data. The patient-specific tandem applicator design improved coverage of the CTV with respect to the existing HDR method without exceeding any OAR dose constraints when applied to the 2D phantom, the 2D patient data, or the 3D patient data. The present disclosure is directed to a computer-implemented method for designing a patient-specific brachytherapy (BT) tandem applicator. The method is implemented using at least one processor in communication with at least one memory. The method includes receiving a radiation treatment plan for treating a region of interest. The radiation treatment plan includes a prescribed radiation dosage and patient anatomical data of the region of interested to be treated. The method also includes applying an inverse planning optimization model to determine an optimal thickness of an interior surface of the tandem applicator at a plurality of dwell positions within the region of interest. The method also includes generating a schedule of dwell times for the tandem applicator based on the generated position-dependent thickness profile. The method also includes transmitting design instructions to a 3D printer for fabrication of the tandem applicator. 1. A patient-specific intensity-modulated high dose rate (HDR) brachytherapy applicator for administering an HDR brachytherapy treatment to a patient, the applicator comprising a plurality of shielding segments distributed along a central longitudinal axis, each shielding segment corresponding to one dwell position and comprising a shielding wall, each shielding wall comprising a plurality of equiangular shielding sections of varying thickness distributed circumferentially about the central longitudinal axis, each equilangular shielding section comprising a shielding thickness, wherein each shield thickness of each equiangular shielding section at each shielding segment is configured to transmit radiation from an HDR source positioned within each shielding segment into the patient at a predetermined dose rate distribution to administer the HDR brachytherapy treatment. 2. The applicator of 3. The applicator of 4. The applicator of 5. The applicator of 6. A computer-implemented method for designing a patient-specific intensity-modulated high dose rate (HDR) brachytherapy applicator for administering an HDR brachytherapy treatment to a patient, the applicator comprising a plurality of shielding segments distributed along a central longitudinal axis, each shielding segment comprising a plurality of equiangular shielding sections distributed circumferentially about the central longitudinal axis, the method implemented using at least one processor in communication with at least one memory, the method comprising:
receiving, by a computing device, a radiation treatment plan for administering the HDR brachytherapy treatment, the radiation treatment plan comprising a prescribed radiation dosage to be delivered to a region of interest and patient anatomical data representative of the region of interest to be treated; determining, by the computing device, an optimal shielding thickness profile and a plurality of optimal dwell times using an inverse planning optimization model constrained by the radiation treatment plan, each optimal dwell time corresponding to one dwell position, each dwell position corresponding to one shielding segment, and the optimal thickness profile comprising a plurality of shield thicknesses, each shield thickness corresponding to one equiangular shielding section of one shielding segment; generating a dwell position-dependent shielding thickness profile comprising the positions of the plurality of the shielding segments and each shield thickness of each equiangular shielding section at each shielding segment; and transmitting, by the computing device, design instructions to a three dimensional (3D) printer for fabrication of the applicator, wherein the design instructions include at least the dwell position-dependent shielding thickness profile. 7. The computer-implemented method of calculating, by the computing device, a plurality of radiation dose rate maps and a plurality of transmission rate maps, each radiation dose rate map and each transmission rate map corresponding to one dwell position of the plurality of dwell positions; calculating, by the computing device, a radiation dose distribution based on the plurality of radiation dose rate maps, the plurality of transmission rate maps, and the plurality of dwell times, the radiation dose distribution comprising a spatial map of a cumulative amount of radiation delivered from a HDR source positioned at each dwell position for each corresponding dwell time; minimizing a cost function by alternately varying the plurality of dwell times with the plurality of transmission rate maps held constant and varying the plurality of transmission rate maps with the plurality of dwell times held constant; and calculating the optimal shielding thickness profile and the plurality of optimal dwell times based on the plurality of transmission rate maps and the plurality of dwell times determined to minimize the cost function. 8. The computer-implemented method of 9. The computer-implemented method of generating a dwell position-dependent dwell time schedule for the applicator comprising the plurality of dwell positions and a corresponding plurality of optimal dwell times; and transmit the dwell position-dependent dwell time schedule to a treatment device for administering the HDR brachytherapy treatment to the patient using the applicator. 10. The computer-implemented method of 11. The computer-implemented method of 12. A computing device for designing a patient-specific intensity-modulated high dose rate (HDR) brachytherapy applicator for administering an HDR brachytherapy treatment to a patient, the applicator comprising a plurality of shielding segments distributed along a central longitudinal axis, each shielding segment comprising a plurality of equiangular shielding sections distributed circumferentially about the central longitudinal axis, the computing device including at least one processor in communication with at least one memory device, the at least one processor programmed to:
receive a radiation treatment plan for administering the HDR brachytherapy treatment, the radiation treatment plan comprising a prescribed radiation dosage to be delivered to a region of interest and patient anatomical data representative of the region of interest to be treated; determine an optimal shielding thickness profile and a plurality of optimal dwell times using an inverse planning optimization model constrained by the radiation treatment plan, each optimal dwell time corresponding to one dwell position, each dwell position corresponding to one shielding segment, and the optimal thickness profile comprising a plurality of shield thicknesses, each shield thickness corresponding to one equiangular shielding section of one shielding segment; generate a dwell position-dependent shielding thickness profile comprising the positions of the plurality of the shielding segments and each shield thickness of each equiangular shielding section at each shielding segment; and transmit design instructions to a three dimensional (3D) printer for fabrication of the applicator, wherein the design instructions include at least the dwell position-dependent thickness profile. 13. The computing device of calculating, by the computing device, a plurality of radiation dose rate maps and a plurality of transmission rate maps, each radiation dose rate map and each transmission rate map corresponding to one dwell position of the plurality of dwell positions; calculating, by the computing device, a radiation dose distribution based on the plurality of radiation dose rate maps, the plurality of transmission rate maps, and the plurality of dwell times, the radiation dose distribution comprising a spatial map of a cumulative amount of radiation delivered from a HDR source positioned at each dwell position for each corresponding dwell time; minimizing a cost function by alternately varying the plurality of dwell times with the plurality of transmission rate maps held constant and varying the plurality of transmission rate maps with the plurality of dwell times held constant; and calculating the optimal shielding thickness profile and the plurality of optimal dwell times based on the plurality of transmission rate maps and the plurality of dwell times determined to minimize the cost function. 14. The computing device of 15. The computing device of generate a dwell position-dependent dwell time schedule for the applicator comprising the plurality of dwell positions and a corresponding plurality of optimal dwell times; and transmit the dwell position-dependent dwell time schedule to a treatment device for administering the HDR brachytherapy treatment to the patient using the applicator. 16. The computing device of 17. The computing device of 18. A high-dose radiation (HDR) modulating system configured to improve target coverage of tumor volume during an HDR treatment, the HDR modulating system including:
a patient-specific intensity-modulated high dose rate (HDR) brachytherapy applicator comprising a plurality of shielding segments distributed along a central longitudinal axis, each shielding segment comprising a plurality of equiangular shielding sections distributed circumferentially about the central longitudinal axis, the plurality of shielding segments defining a central lumen extending along the central longitudinal axis, each shielding segment further defining a dwell position within the central lumen; and an HDR source movably insertable into the central lumen during an HDR treatment, the HDR source configured to reside at each dwell position within each shielding segment for a corresponding dwell time, wherein each corresponding dwell time is based on a radiation therapy plan; wherein each equiangular shielding section at each shielding segment comprises a shield thickness configured to transmit radiation from the HDR source residing at each dwell position at a predetermined dose rate distribution. 19. The system of 20. The system of CROSS-REFERENCE TO RELATED APPLICATION
FIELD OF THE INVENTION
BACKGROUND
BRIEF DESCRIPTION
BRIEF DESCRIPTION OF THE DRAWINGS
DETAILED DESCRIPTION
Input: z0= (t0, T0), (So)oϵO, St, ST, C, λt0, λT0, tolin, tolout, {circumflex over (L)}. While 1. Solve t subproblem: Set nt= 1, and do (a) λt= linesearch(F, tn 2. tn= tn 3. Solve T subproblem: Set nT= 1, and do (a) if λt≥ {circumflex over (L)}: λT= linesearch(F, tn, Tn else: Perturb Tn 4. Tn= Tn 5. zn= (tn, Tn), n = n + 1. Output: zn= (tn, Tn). Input: t0, T0, (So)o∈O, St, δS, C, λt0, λT0, tolin, tolout. While the calculated dose volume D does not satisfy the treatment target do 1. (tn,Tn) = Algorithm 1(tn−1,Tn−1,So,St,C,λt0,λT0). 2. Calculate the dose volume based on t and T. 3. for o = 1, ..., No, if D does not satisfy the condition for the o-th OAR So= So+ δS. Output: tn,Tn. Using the Calculated Parameters to Create a Patient-Specific Applicator
Computing System
Therapeutic Methods
Administration
Kits
EXAMPLES
Example 1
Validation of a Method for Fabricating a Patient-specific IMBT Tandem HDR Applicator for Cervical Cancer Using a 3D Printer
Example 2
Validation of a Method for Designing a Patient-specific IMBT Tandem HDR Applicator Using an Inverse Planning Optimization Model
Dose Rate Calculation (Based on AAPM TG-43 Report)
β=Linesearch(
γ=Linesearch(Dose Constraints and Parameters
Dose Constraints Assigned to Structures in the 2D Phantom and 2D Patient Cases Structure Dose Constraint CTV ≥550 cGy Bladder ≤460 cGy Rectum ≤420 cGy Sigmoid ≤420 cGy 2D Phantom Criterions checking Bladder Rectum Sigmoid CTV Existing HDR 100% 100% 100% 58.32% Patient-Specific 100% 100% 100% 99.18% Tandem Applicator Method 2D Patient Data Criterions checking Bladder Rectum Sigmoid CTV Existing HDR 100% 100% 100% 56.21% Patient-Specific 100% 100% 100% 99.92% Tandem Applicator Method




































