The embodiments described herein relate to the field of x-ray and particle radiation dosimetry, and in particular to two-dimensional (planar) and three-dimensional (volumetric) dosimetry using tomographic reconstruction algorithms.
Radiation therapy, or “radiotherapy”, is the medical use of ionizing radiation, for example to control malignant cells in cancer treatments. Most frequently, radiation therapy makes use of x-ray or particle beams originating from outside a patient's body (i.e. external radiotherapy).
Because of their inherent complexity, some more advanced modes of external beam radiation treatment, such as Intensity Modulated Radiation Therapy (IMRT), helical tomotherapy, radiosurgery and Intensity Modulated Arc Therapy (IMAT), are verified experimentally before the treatment is administered. In particular, the dose distribution may at least be measured in a plurality of positions in space, such as in a two-dimensional plane.
Most often, such dosimetry measurements are performed by an array of periodic detectors. Such detectors can be stacked in a 2D assembly and are usually composed of ionization chambers, semi-conductor diodes (e.g. the Sun Nuclear Mapcheck™), scintillation detectors, or thermoluminescent dosimeters (TLD).
Ionization chamber assemblies are usually limited in spatial resolution because of their significant detector size (usually many millimeters squared). Moreover, the materials used in the fabrication of ionization chambers and diodes (such as silicon and air) are normally not water-equivalent (neither by chemical composition or physical density), and as such their presence in the radiation field tends to perturb the radiation fluence and the resulting dose measurement.
Scintillation detector assemblies on the other hand are usually limited in the number of detectors that can be conveniently stacked in a single array, thus reducing the area of detection or the resolution of the array.
Finally, each TLD normally needs to be read individually after irradiation, which limits the use of a large number of TLDs in a 2D array.
In other cases, direct two-dimensional dose measurement can be conducted by using high resolution radiographic or radiochromic films. However such techniques are subject to several drawbacks. For example, the need to develop radiographic films before they are readable limits their speed and makes their use for online dosimetry generally impractical. Moreover, the development process can affect the film response. Radiographic films are also prone to over-responding to low-energy photons, are not water-equivalent, and their precision is often limited to ±5%.
Radiochromic films on the other hand may be closer to water-equivalence in the energy range used in external beam radiotherapy and do not require developing. However, radiochromic films are temperature dependant and sensitive to ultraviolet light, and achieving better than 5% reproducibility in a routine fashion with radiochromic film tends to be difficult.
Direct two-dimensional dose measurement can also be conducted by a scintillator plate as described for example in U.S. Patent Application Publication Nos. 2004/0178361 (to Fontbonne) and 2009/0218495 (to Leroux et al). However, such methods require extensive correction of the system readings in order to correct for inhomogeneous light dispersion and {hacek over (C)}erenkov radiation. Accordingly, the precision of such assemblies tends to be limited, and errors of ±10% on the dose reading are typical.
The above-mentioned detectors allow, at most, measurements in two-dimensional planes. Accordingly, 3D dosimetry with these detectors usually involves summing many 1D or 2D measurements of the dose distribution.
Currently, known detectors that can be used for three-dimensional measurements are dosimetric gels. Dosimetric gels are either based on the behavior of ferrous ions or on the polymerization of a monomer. They can be produced using a large variety of chemical formulas and each has its own set of advantages and disadvantages. However, most gels are made using a delicate fabrication process and require time-consuming development processes that makes them unsuitable for online measurements.
The embodiments herein will now be described, by way of example only, with reference to the following drawings, in which:
It will be appreciated that numerous specific details are set forth in order to provide a thorough understanding of the exemplary embodiments described herein. However, it will be understood by those of ordinary skill in the art that the embodiments described herein may be practiced without these specific details. In particular, in some instances, well-known methods, procedures and components have not been described in detail so as not to obscure the embodiments described herein.
Furthermore, this description is not to be considered as limiting the scope of the embodiments described herein in any way, but rather as merely describing the implementation of various embodiments as described herein. In addition, while certain features have been illustrated and described herein, many modifications, substitutions, changes, and equivalents will now occur to those of ordinary skill in the art. It is therefore to be understood that the appended claims are intended to cover such modifications and changes.
According to one aspect, there is provided a planar and volumetric dosimeter for use with a radiotherapy machine having a radiation source, the planar and volumetric dosimeter comprising a scintillating assembly adapted to generate a light output in response to an incident radiation distribution, and a photo-detector adapted to convert the transmitted optical energy to electrical signals for calculating a two-dimensional (2D) or three-dimensional (3D) incident dose distribution on the scintillating assembly by using tomographic reconstruction algorithms.
The scintillating assembly may include a plurality of scintillating optical fibers, usually in a sizeable length (e.g. at least a few centimeters long), each scintillating optical fiber acting as a weighted dose line integrator. These scintillating optical fibers may be embedded in a phantom medium. The phantom medium may be made of a material with properties relative to radiation that are similar to the material used in the scintillating optical fibers.
In some embodiments, a plurality of collection optical fibers may be coupled to one or both of the opposing ends of the scintillating optical fibers. In such cases one or more photo-detectors may be coupled to the collection optical fibers and be operable to convert optical energy transmitted by the collection optical fibers into electrical signals.
In some embodiments, the scintillation optical fibers of the scintillating assembly may be disposed in one or more imaging planes. These imaging planes may be disposed and oriented following different geometries in the dosimeter volume (e.g. in a series of parallel disks or in a series of concentric cylinders or “rings”). The light intensities collected in one imaging plane may be configured to represent weighted dose line integrals of the incident 2D dose distribution on that imaging plane. By using tomographic reconstruction algorithms on these weighted dose line integrals, the incident 2D dose distribution on the imaging plane may be calculated.
In some embodiments, the 3D dose distribution incident on at least a part of the planar and volumetric dosimeter may be determined by the summation or interpolation (or both) of a plurality of 2D dose distributions obtained by one or more imaging plane(s).
In some embodiments, the scintillating optical fibers may be conditioned with UV light in order to increase their optical attenuation.
In some embodiments, both the scintillating assembly and the photo-detector may be rotatable about a given axis of rotation. This axis of rotation may be parallel to the radiation beam axis, perpendicular (orthogonal) to the radiation beam axis, or in some embodiments generally oriented in any spatial direction.
The scintillating assembly may be made of a water-equivalent material with properties relative to radiation that are equivalent to water.
In some embodiments, the scintillating assembly may be located on a treatment table at or near a height of treatment delivery.
In some examples, the light collected by each end of each scintillating optical fiber may be determined by a summation of the contribution of each element of the dose pattern volume according to the following equations:
where Ij,± represents the light collected by a single end (+ or −) of the jth scintillating optical fiber, Di represents the dose deposited in the ith volume element, Fj,i represents the percentage of the ith volume element covered by the jth scintillating optical fiber, κj represents the scintillation efficiency of the jth scintillating optical fiber (i.e. the number of scintillation photons emitted per unit of deposited dose), λj accounts for the differential light attenuation along the jth scintillating optical fiber, Cj,± represent the light losses due to the optical coupling of the jth scintillating optical fiber to the photo-detector, xi, yi and zi represent the spatial coordinate of the ith volume element, ui represents the position on the scintillating optical fiber with respect to the reference point on the fiber and f represents the function relating ui to the coordinates xi, yi and zi. The preceding equation can also be written in the following vectorial form:
Ij,±=Aj,±·{right arrow over (D)}
where {right arrow over (D)} represents the dose distribution in vectorial form and Aj,± represents the row of the projection matrix associated with the + or − side of the jth scintillating optical fiber.
According to another aspect there is provided a method for determining the weighted dose line integral along a single scintillating optical fiber from the light collected by a single end of the scintillating optical fiber according to the following equations:
where pj represents the weighted dose line integral for the jth scintillating optical fiber, Ij,± represents the light collected by a single end (+ or −) of the jth scintillating optical fiber when irradiated by a given dose distribution, INorm,j,± represents the light collected by the same end (+ or −) of the jth scintillating optical fiber when irradiated by a known normalization dose distribution, {right arrow over (D)}norm represents the known normalization dose distribution in vectorial form,
where pj represents the weighted dose line integral for the jth scintillating optical fiber, Ij,± represents the light collected by the + or − end of the jth scintillating optical fiber when irradiated by a given dose distribution, INorm,j,± represents the light collected by the + or − end of the jth scintillating optical fiber when irradiated by a known normalization dose distribution, {right arrow over (D)}norm represents the known dose distribution in vectorial form, {right arrow over (D)} represents the measured dose distribution in vectorial form, Bj represents the row of the modified projection matrix associated with the jth scintillating optical fiber, Fj,i represents the percentage of the ith volume element covered by the jth scintillating optical fiber, κj represents the scintillation efficiency of the jth scintillating optical fiber (i.e. the number of scintillation photons emitted per unit of deposited dose), λj accounts for the differential light attenuation along the jth scintillating optical fiber, ui represents the position on the scintillating optical fiber with respect to the reference point on the fiber and Di represents the dose deposited in the ith volume element of the measured dose distribution.
According to another aspect there is provided a method for determining the incident 2D dose distribution incident on an imaging plane composed of a plurality of scintillating fibers for use on a radiotherapy machine having a radiation source, the method comprising using iterative reconstruction techniques on the measured weighted dose line integrals (such as an algebraic reconstruction algorithm or an expectation-maximization algorithm) in conjunction with a total variation (TV) minimization of the resultant dose pattern, the total variation (TV) of the 2D dose distribution being defined as:
where Di represents the dose measured in the ith volume element and m is the number of such volume elements in each row of the 2D dose distribution.
According to another aspect there is provided a method for determining the incident 2D dose distribution on an imaging plane composed of a plurality of scintillating fibers for use on a radiotherapy machine having a radiation source, the method comprising using linear programming algorithms to determine the 2D dose distribution under the constraints that the measured weighted dose line integrals should be respected and that the reconstructed dose distribution should be physically possible and noiseless.
According to some embodiments, a planar and volumetric dosimeter as described herein includes a scintillating assembly and a photo-detector configured to determine the 2D or 3D dose distribution incident on the scintillating assembly. Generally, the scintillating assembly is situated under the primary beam of a medical linear accelerator (also known as a “linac”) operable to generate x-rays or particle beams.
In some embodiments, the scintillating assembly includes a plurality of scintillating optical fibers assembled in a phantom volume (i.e. a non-scintillating material).
In some embodiments, the scintillating optical fibers may have a length of a few centimeters or more.
In some embodiments, the scintillating optical fibers are conditioned using UV light in order to increase their optical attenuation.
In some embodiments, the phantom volume is made of plastic or other material that is equivalent to the composition of the scintillating optical fibers.
In some embodiments, the phantom volume is made of a material with at least some properties relative to radiation that are equivalent to water (e.g. a water-equivalent phantom material).
In some embodiments, the scintillating optical fibers of the scintillating assembly are grouped in imaging planes in which a plurality of scintillating optical fibers are aligned parallel to each other in a single planar region within the dosimeter volume. The light intensities collected by all scintillating optical fibers in a given imaging plane are configured to represent weighted dose line integrals of the 2D dose distribution generated in that imaging plane.
In some embodiments, the imaging plane(s) may rotate around an axis of rotation that is perpendicular or normal to the plane(s), obtaining a plurality of weighted dose line integrals of the 2D dose distribution generated in the imaging plane at different angles between the scintillating optical fibers and the dose distribution. Accordingly, by using tomographic reconstruction algorithms on the acquired weighted dose line integrals of the 2D dose distribution, the 2D dose distribution incident on the imaging plane can be determined.
Furthermore, by summing a plurality of such imaging planes, or by translating one or more imaging planes in the dosimeter area, one can obtain a plurality of 2D dose distributions incident on the imaging plane simultaneously. By summing and/or interpolating these 2D dose distributions, one can calculate the 3D dose distribution incident on at least a part of the dosimeter volume.
In some embodiments, the scintillating optical fibers may be grouped in circular imaging planes or “rings”, in which a plurality of scintillating optical fibers are grouped in a single cylindrical region or ring in the scintillating assembly. Each scintillating optical fiber may be oriented in the circular imaging plane following a distinct angle in the circular imaging plane. For example, the angle of a particular scintillating optical fiber in the circular imaging plane may be distinct from the angles of the other scintillating optical fibers in the same circular imaging plane.
In some embodiments, the circular imaging planes may be adapted to rotate around their central axis of symmetry. The sum of the weighted dose line integrals of a single scintillating optical fiber after a complete rotation (i.e. 360 degrees) or at least a portion of a complete rotation (i.e. 180 degrees) of the circular imaging plane around its central axis of symmetry may represent a line (1D) projection of the 2D dose distribution generated in the circular imaging plane at an angle defined by the orientation of the scintillating optical fiber in the circular imaging plane. By using tomographic reconstruction algorithms on the acquired weighted dose line integrals of the 2D dose distribution, the 2D dose distribution incident on the circular imaging plane can be determined. Furthermore, by summing a plurality of such circular imaging planes, one can obtain a plurality of 2D dose distributions incident on the imaging plane simultaneously. By summing and interpolating these 2D dose distributions, one can calculate the 3D dose distribution incident on at least a part of the dosimeter volume.
Also generally described herein is a theoretical model for predicting the optical energy produced by scintillating optical fibers and collected at each end thereof as a function of the measured dose distribution.
Also generally described herein are methods of determining the weighted dose line integral along scintillating optical fibers as a function of the optical energy produced by scintillating optical fibers and collected at each end thereof.
Also generally described herein are methods of determining the incident 2D dose distribution on an imaging plane of the planar and volumetric dosimeter. A first determination method uses algebraic reconstruction techniques on the calculated weighted dose line integrals and a total variation (TV) minimization of the resultant dose pattern. A second determination method uses linear programming optimization of the calculated weighted dose line integrals in parallel with physical assumptions on the 2D dose distribution geometry.
Turning now to
The phantom material 21 and the imaging plane 19 are positioned under the radiation source (indicated schematically as 14) so that the axis of rotation 15 of the planar and volumetric dosimeter 10 is generally aligned with or parallel to the radiation beam axis 22.
In some embodiments, at least one of the imaging plane 19, the collection optical fibers 20 and the photo-detector 12 can undergo a vertical translation 18 (e.g. generally along a direction aligned with the radiation beam axis 22) while the phantom material 21 remains static (i.e. is held in place).
In some embodiments, the light fluxes 24+ and 24− emitted by the scintillating fibers 23 in the imaging plane 19 and converted to electrical signals by the photo-detector 12 are configured to represent weighted dose line integrals of the 2D dose distribution incident on the imaging plane 19.
When moving along this circular trajectory 16, each imaging plane 19 is able to take a plurality of weighted dose line integrals of the 2D dose distribution incident on the imaging plane 19 at different orientations between the scintillating fibers 23 and the deposited dose distribution. By using tomographic reconstruction algorithms on this plurality of weighted dose line integrals, the 2D dose distribution incident on the imaging plane 19 can be determined.
Furthermore, by applying a vertical translation 18 to each of the imaging planes 19 in addition to the rotation of the planar and volumetric dosimeter 10 around the axis of rotation 15, a plurality of 2D dose distributions at different depths in the phantom material 21 can be acquired simultaneously. By summing or interpolating (or both) between the obtained 2D dose distributions at different depths in the phantom material 21, the 3D dose distribution incident on at least a part of the planar and volumetric dosimeter 10 can be determined.
Turning now to
As shown, the planar and volumetric dosimeter 110 includes a scintillating assembly composed of a plurality of circular imaging planes 125 (as shown on
Each scintillating optical fiber 123 may have one end exposed to the photo-detector 112, or both ends exposed to the photo-detectors 112a and 112b. The circular imaging planes 125, the phantom material 121 and the photo-detector(s) 112 can all rotate together about a central axis of rotation 115. For clarity, a perspective view of the planar and volumetric dosimeter 10 is shown in
As shown, the circular imaging plane 125 includes the plurality of scintillating optical fibers 123, with each scintillating fiber 123 inclined at a different angle 127 with respect to the faces 128a and 128b of the phantom material 121. Thus, after a complete rotation (e.g. 360 degrees) of the planar and volumetric dosimeter 110 around its central axis of rotation 115, the sum of all the light fluxes emitted by a single scintillating optical fiber 123 and converted to a weighted dose line integral by the photo-detector(s) 112 represents a line (1D) projection of the 2D dose distribution incident on the circular imaging plane 125 in which the scintillating optical fiber belongs.
With a plurality of line (1D) projections obtained from the plurality of scintillating optical fibers 123 and using tomographic reconstruction algorithms, the 2D dose distribution incident on each circular imaging plane 125 may be determined. Furthermore, by using a plurality of circular imaging planes 125, a plurality of 2D dose distributions at different radius from the central axis of rotation 115 can be acquired simultaneously. By summing or interpolating (or both) between the obtained 2D dose distributions at different radii in the phantom material 121, the 3D dose distribution incident on at least a part of the planar and volumetric dosimeter 110 can be determined.
Turning now to
By using such a layout of scintillating optical fibers 123 in the circular imaging plane 125, only half a rotation (e.g. 180 degrees) of the dosimeter 110 around its central axis of rotation 115 may be required to obtain 1D projections (e.g. line projections) of the 2D dose distribution incident on the circular imaging plane 125, for a total of one 1D projection per pair of scintillating optical fibers sharing the same unique angle 127a.
Also described herein is a model that predicts the optical energy produced by scintillating optical fibers (e.g. scintillating optical fibers 123) when irradiated by a given dose distribution 129. From this model and following the deposition of a dose distribution 129, as shown in
where Di represents the dose deposited in the ith volume element 130, Fj,i represents the percentage of the ith volume element 130 covered by the jth scintillating optical fiber 123, κj represents the scintillation efficiency of the jth scintillating optical fiber 123 (i.e. the number of scintillation photons emitted per unit of deposited dose), λj accounts for the differential light attenuation along the scintillating optical fiber 123, Cj,± represent the light losses due to the optical coupling of the jth scintillating optical fiber 123 to the photo-detector, xi, yi and zi represent the spatial coordinate of the ith volume element 130, ui represents the position on the scintillating optical fiber 123 with respect to a reference point on the fiber 132 (e.g. the position on the scintillating optical fiber 123 closest to the axis of rotation 131) and f represents the function relating ui to the coordinates xi, yi and zi. In some embodiments, this reference point on the fiber 132 can be defined as the center of the scintillating optical fiber 123. Equation 1 can also be written in the following vectorial form:
Ij,±=Aj,±·{right arrow over (D)} (1a)
where
One method for calibrating a planar and volumetric dosimeter is described in U.S. Provisional Patent Application No. 61/375,125 filed on Aug. 19, 2010 (now U.S. Non-Provisional application Ser. No. 13/207,624) entitled FLUENCE MONITORING DEVICES WITH SCINTILLATING FIBERS FOR X-RAY RADIOTHERAPY TREATMENT AND METHODS FOR CALIBRATION AND VALIDATION OF SAME, the entire contents of which are hereby incorporated herein by reference for all purposes.
Also described herein is a method for determining the weighted dose line integral along a single scintillating optical fiber from the light collected by a single end of the scintillating optical fiber according to the following equation:
where pj represents the weighted dose line integral for the jth scintillating optical fiber, Ij,± represents the light collected by a single end (+ or −) of the jth scintillating optical fiber when irradiated by a given dose distribution, INorm,j,± represents the light collected by the same end (+ or −) of the jth scintillating optical fiber when irradiated by a known normalization dose distribution, {right arrow over (D)}norm represents the known dose distribution in vectorial form, {right arrow over (D)} represents the measured dose distribution in vectorial form, Bj,± represents the row of the modified projection matrix associated with the + or − side of the jth scintillating optical fiber, represents the percentage of the ith volume element covered by the jth scintillating optical fiber, κj represents the scintillation efficiency of the jth scintillating optical fiber (i.e. the number of scintillation photons emitted per unit of deposited dose), λj accounts for the differential light attenuation along the jth scintillating optical fiber, ui represents the position on the scintillating optical fiber with respect to the reference point on the fiber and Di represents the dose deposited in the ith volume element of the measured dose distribution. It will be understood that the projection calculated using equation (3) and (3a) will normally be dependant on the light attenuation in the scintillating optical fiber, so that pj will be different if calculated from the + or the − side of the scintillating optical fiber.
Also described herein is a method for determining the weighted dose line integral along a single scintillating optical fiber from the light collected by the two opposed ends of the scintillating optical fiber according to the following equation:
where pj represents the weighted dose line integral for the jth scintillating optical fiber, Ij,± represents the light collected by the + or − end of the jth scintillating optical fiber when irradiated by a given dose distribution, INorm,j,± represents the light collected by the + or − end of the jth scintillating optical fiber when irradiated by a known normalization dose distribution, {right arrow over (D)}norm represents the known dose distribution in vectorial form, {right arrow over (D)} represents the measured dose distribution in vectorial form, Bj,± represents the row of the modified projection matrix associated with the + or − side of the jth scintillating optical fiber, Fj,i represents the percentage of the ith volume element covered by the jth scintillating optical fiber, κj represents the scintillation efficiency of the jth scintillating optical fiber (i.e. the number of scintillation photons emitted per unit of deposited dose), λ accounts for the differential light attenuation along the jth scintillating optical fiber, ui represents the position on the scintillating optical fiber with respect to the reference point on the fiber and Di,norm represents the dose deposited in the ith volume element of the known normalization dose distribution (in vectorial form). It will be understood that equation (4) and (4a) represent an approximation that is normally valid only when the λj term in equation (1), accounting for differential light attenuation along the scintillating optical fiber, is linear in u, and that under this approximation the projection calculated using equation (4) will be approximately independent of the light attenuation in the scintillating optical fiber.
Also provided herein are two exemplary methods of determination of the incident dose distribution on an array of scintillating optical fibers (e.g. the imaging planes 19, or 125) rotating around a pre-determined axis of rotation (e.g. central axis of rotations 15, 115 or 131). A first determination method uses iterative reconstruction techniques on the measured weighted dose line integrals in conjunction with a total variation (TV) minimization of the resultant dose pattern, with the total variation (TV) of the 2D dose distribution being defined as:
where Di represents the dose measured in the ith volume element 130 as shown in
The iterative reconstruction technique used in this first method of determination can be an algebraic reconstruction technique (ART). For example, suppose the incident dose distribution {right arrow over (D)} (in vectorial form) can be described using:
pj=Bj·D (6)
where Bj represents the row of the modified projection matrix associated with pj, which represents the jth weighted dose line integral of the incident dose distribution measured as described herein. Using this technique the nth iteration of the reconstructed dose distribution ({right arrow over (D)}(n)) can be obtained from the previous iteration ({right arrow over (D)}(n−1)) by the following formula:
where βn is the control parameter of the nth iteration, Bj is the jth row of the projection matrix associated with pj, the jth weighted dose line integral measured as described herein.
The iterative reconstruction technique used in this first method of determination can also be an expectation-maximization technique (EM). For example, the nth iteration of the reconstructed dose distribution ({right arrow over (D)}(n)) can be obtained from the previous iteration ({right arrow over (D)}(n−1)) by the following formula:
where Di(n) represents the dose measured in the ith volume element 130 in the nth iteration of the reconstructed dose distribution, bj,i is the element of the matrix B associated with the ith volume element and the jth weighted dose line integral measured as described herein, Bj is the jth row of the projection matrix associated with pj, the jth weighted dose line integral measured as described herein.
The total variation (TV) minimization technique using this first method of determination can be a steepest-descent method. For example, the total variation of the reconstructed dose distribution ({right arrow over (D)}(n)) can be diminished by the following equation:
where Di(n′) represents the dose measured in the ith volume element 130 in the nth iteration of the reconstructed dose pattern after TV minimization, αn is the control parameter of the nth iteration of the TV minimization and TV is calculated using Equation 5.
The two control parameters αn or βn (or both) can be adjusted so that the dose distribution modification caused by the iterative reconstruction technique is of the same magnitude as the dose pattern modification caused by the TV minimization. It will be understood that, by performing a plurality of iterations in the first determination method described above, the reconstructed dose pattern will tend to more closely approach the actual delivered dose pattern, under the assumption that the incident dose pattern has a low total variation.
A second determination method uses linear programming algorithms to determine the dose in each volume element 130 by linear optimization. For example and using the representations of Equations 6 and 6a, the optimization problem used in this determination method may include the following:
where Rj is the jth element of the residual vector that represents the difference between the projection Bj·{right arrow over (D)} calculated from the reconstructed incident dose distribution and the measured weighted dose line integral pj, and r generally represents the maximum allowed difference between the projection Bj·{right arrow over (D)} calculated from the reconstructed incident dose pattern and the measured weighted dose line integral pj.
The physical dose distribution constraints generally ensure that the reconstructed dose pattern is physically acceptable. Such constraints can ensure that the dose distribution is smooth (i.e. not noisy) and/or that the dose gradient between two adjacent volume elements is inferior to a certain threshold. In some cases, such constraints can be represented mathematically for example by:
where the indices k applies to the N neighboring volume elements of the ith volume element Di, and hmax and η are constants regulating the physical dose distribution constraints. Such constants hmax and η can be adjusted based on a previous reconstruction of the dose pattern or a theoretical dose pattern, in order for the constraints to be more permissive in regions of high dose gradient and more severe in regions of dose plateau.
While certain features have been illustrated and described herein, many modifications, substitutions, changes, and equivalents will now occur to those of ordinary skill in the art. It is, therefore, to be understood that the appended claims are intended to cover all such modifications and changes.
This application claims the benefit of U.S. Provisional Patent Application Ser. No. 61/470,817 filed Apr. 1, 2012 and entitled PLANAR AND VOLUMETRIC DOSIMETER WITH SCINTILLATING MATERIAL FOR RADIOTHERAPY TREATMENT USING TOMOGRAPHIC RECONSTRUCTION, the entire contents of which are hereby incorporated by reference herein for all purposes.
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/CA2012/000283 | 3/29/2012 | WO | 00 | 2/28/2014 |
Publishing Document | Publishing Date | Country | Kind |
---|---|---|---|
WO2012/129661 | 10/4/2012 | WO | A |
Number | Name | Date | Kind |
---|---|---|---|
5856673 | Ikegami | Jan 1999 | A |
20020131552 | Nishizawa et al. | Sep 2002 | A1 |
20060153341 | Guyonnet et al. | Jul 2006 | A1 |
Number | Date | Country | |
---|---|---|---|
20140217295 A1 | Aug 2014 | US |
Number | Date | Country | |
---|---|---|---|
61470817 | Apr 2011 | US |