EP3460530A1 - One dimensional transmission detector for radiotherapy - Google Patents

One dimensional transmission detector for radiotherapy Download PDF

Info

Publication number
EP3460530A1
EP3460530A1 EP18195642.6A EP18195642A EP3460530A1 EP 3460530 A1 EP3460530 A1 EP 3460530A1 EP 18195642 A EP18195642 A EP 18195642A EP 3460530 A1 EP3460530 A1 EP 3460530A1
Authority
EP
European Patent Office
Prior art keywords
strip
mlc
sensing device
strips
dose
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
EP18195642.6A
Other languages
German (de)
French (fr)
Other versions
EP3460530B1 (en
Inventor
David Menichelli
Lutz MÜLLER
Jochen KRIMMER
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Ion Beam Applications SA
Original Assignee
Ion Beam Applications SA
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Ion Beam Applications SA filed Critical Ion Beam Applications SA
Publication of EP3460530A1 publication Critical patent/EP3460530A1/en
Application granted granted Critical
Publication of EP3460530B1 publication Critical patent/EP3460530B1/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01TMEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/29Measurement performed on radiation beams, e.g. position or section of the beam; Measurement of spatial distribution of radiation
    • G01T1/2914Measurement of spatial distribution of radiation
    • G01T1/2921Static instruments for imaging the distribution of radioactivity in one or two dimensions; Radio-isotope cameras
    • G01T1/2935Static instruments for imaging the distribution of radioactivity in one or two dimensions; Radio-isotope cameras using ionisation detectors
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1048Monitoring, verifying, controlling systems and methods
    • A61N5/1075Monitoring, verifying, controlling systems and methods for testing, calibrating, or quality assurance of the radiation treatment apparatus
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1042X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy with spatial modulation of the radiation beam within the treatment head
    • A61N5/1045X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy with spatial modulation of the radiation beam within the treatment head using a multi-leaf collimator, e.g. for intensity modulated radiation therapy or IMRT
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1048Monitoring, verifying, controlling systems and methods
    • A61N5/1071Monitoring, verifying, controlling systems and methods for verifying the dose delivered by the treatment plan
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1077Beam delivery systems
    • A61N5/1081Rotating beam systems with a specific mechanical construction, e.g. gantries
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N2005/1085X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy characterised by the type of particles applied to the patient
    • A61N2005/1087Ions; Protons

Definitions

  • the present invention relates to radiation therapy equipment, and in particular, to devices adapted for enabling quality assurance of external beam radiation therapy for devices using a multi-leaf collimator.
  • This treatment plan prepared following the prescription of the radiation oncologist, requires precise, accurate, and consistent delivery of an optimal radiation dose to the tumor while simultaneously reducing the side effects on healthy surrounding tissues.
  • state-of-the art radiation therapy includes using a multi-leaf collimator to modulate the intensity of the radiation beam produced by a particle beam accelerator, such as a LINAC (conventional radiation-emitting, linear accelerator), to better shape the radiation dose distribution in order to minimize irradiation of non-tumorous cells.
  • a multi-leaf collimator is a device made up of individual "leaves" of metal, typically high atomic-number materials such as tungsten. The leaves move independently in and out of the beam aperture in order to partially block the beam. MLCs typically include two banks of opposing leaves.
  • each bank has from approximately 50 to 200 leaves, each leaf having a width in the range from about 0.25 cm to about 2 cm (projected at isocenter).
  • a given leaf in a bank and the opposite one in the other bank form a "leaf pair"; they can move one against the other (thus obstructing the beam) or apart (creating an aperture which, in most installations, can be also asymmetric with respect to the central beam axis [leaf overtravel]).
  • IMRT Intensity Modulated Radiation Therapy
  • the delivered dose may vary from the intended (planned) treatment dose due to a number of factors. These factors include the wrong treatment planning system (TPS) calculation, an error or fault in the physical movement of the LINAC, including its collimating system, or an incorrect beam energy of the LINAC, errors in patient positioning, or changes in the patient's anatomy.
  • TPS treatment planning system
  • Risks of errors is particularly high in the case of plans with high modulation (i.e. a highly conformal dose distribution is obtained with a complex dynamics of the MLC, in such a way that many small beams contribute to the overall dose).
  • stereotactic treatments in which the overall treatment dose is delivered in one or a few fractions (i.e. treatment sessions) instead of using the more traditional approach of many fractions (e.g. 20-30) delivered in subsequent days.
  • Stereotactic treatments are appealing since they increase workflow efficiency, reduce problems related to inter-fractional anatomy changes and, in some cases, present a radiobiological advantage.
  • delivering a wrong dose may have severe and irreparable consequences.
  • a phantom-based detector rests on the couch (replacing the patient) and is irradiated with the treatment plan.
  • This test reveals calculation errors and systematic LINAC errors, and sometimes resulting dose distributions are projected into the patient anatomy with the assumption that it conforms to the CT scan taken at the time of planning.
  • this test is not effective against lack of LINAC repeatability, positioning errors, and anatomical changes, as it is done in a special QA fraction and not during patient treatment.
  • pre-treatment QA adds an additional step in the clinical workflow, reducing its efficiency.
  • Active sensors typically organized in an array, provide the spatial resolution necessary to verify the complex radiation patterns generated by an MLC.
  • Known ionization chamber arrays include a device described by Bonin et al., "A pixel chamber to monitor the beam performances in hadron therapy" published in Nuclear Instruments and Methods in Physics Research, No. A 519 (2004) at pages 674-686 . This device includes a 2-D array of 1024 successive ionization chambers arranged in a regular matrix of 32 ⁇ 32 pixels.
  • Ionization chambers generally work by providing a high-voltage difference between two electrodes.
  • a gas or liquid fills the space between the electrodes.
  • this gas is air, nitrogen or an organic gas.
  • an electrically insulating liquid fills this space.
  • the medium gets ionized when exposed to radiation. The ions move towards the electrodes and their charge is collected. An external electrometer measures this charge.
  • the ionization of a molecule of the medium requires a known average energy and, depending on the gas and on the irradiation type, the collected charge is directly proportional to the energy deposited in the gas.
  • United States Pat. No. 4,988,866 describes yet another dosimetry device having a limited number of specifically located sensors for measuring.
  • This device has at least two significant shortcomings. First, due to its configuration and limited number of sensors, it cannot verify the structure of a radiation field of any size. Second, the device can only measure a very narrow energy range due to the single absorber.
  • Still another dosimetry device improves the spatial resolution of measurement without increasing the number of individual sensors.
  • a set of sensors arranged in a first line and positioned on a rotatable support are located at a different radius from a second line of sensors.
  • the support successively rotates along an angle of 1° or 2°.
  • Successive measurements by lines of sensors taken at angular steps result in a set of many measurement points, thus giving higher spatial resolution.
  • this device requires a mechanical drive for the sensor, and the measurement is more time-consuming.
  • one because one takes a set of measurements at successive times, one relies on the stability and constancy of the radiation source. This device cannot measure radiation beam energy.
  • Yang, Y. and Xing, L. describe a method for determining the position of a leaf of a MLC in " Using the volumetric effect of a finite-sized detector for routine quality assurance of multi leaf collimatorleaf positioning" Med. Phys. 30 (2003) 433-441 . They disclose a finite-size detector, such as an ion chamber, positioned on one leaf as projected in the isocenter plane. Leaf-position displacement increases or decreases the irradiated volume of the detector; therefore, measurement of dose relates to a leaf position error.
  • the QA measurement repeats for each treatment fraction while the patient lies on the couch (in-vivo QA).
  • in-vivo QA enables treatment interruption if something starts to go wrong.
  • Fluence detectors estimate dose perturbation based on detected errors and modify the revealed dose distribution inside the patient accordingly using a dose engine combined with the patient's CT images.
  • Updated images of the patient e.g. taken with CBCT
  • a first solution uses transmission detectors.
  • the detector positions between the LINAC's collimating system and the patient.
  • the collection of measured data is accomplished by sensors with high-dosimetric performance, such as ion chambers.
  • Müller in U.S. Pat. No. 8160204 issued 2012-04-17 teaches the determination of photon and electron fluence distributions based on measurements. Comparing the expected detector response (calculated from the patient plan on a detector model) to measured data enables a calculation of the difference between actually delivered fluence and the planned fluence.
  • One particularly well suited detector includes the MatriXX brand detector paired with the COMPASS 3D brand software tool, both available from IBA Dosimetry GmbH (Bruhofstrasse 5, 90592 Schwarzenbruck, Germany). This is a system for 3D treatment plan and delivery verification and patient dose analysis, giving clinically relevant data analysis.
  • the software component of this detector reconstructs dose from measured fluence, compares the patient plan with measurements, and provides 3D dose deposition information inside the patient's anatomy. Accordingly, this arrangement enables visual or statistical plan evaluation (evaluating dose differences/gamma relative to TPS inside patient CT or on a structure-by-structure statistical/quantitative basis via comparison of the TPS generated DVHs to that of this device's independently calculated DVHs).
  • Drawbacks of these known detectors include perturbation of the beam under measurement, reduced space between the patient's head and the couch, and loss of the option of mounting additional accessories at the LINAC. Additionally, radiation scattered from the LINAC's head and collimator, which does not contribute to the dose in the patient, complicates the beam measurement, and modeling with available dose engines and predicting the response becomes difficult.
  • a second solution a system based on portal imagers, offers a practical solution because most LINAC devices are equipped with a system for portal imaging (typically an amorphous-silicon flat-panel with high spatial resolution) dedicated to control patient positioning.
  • the QA system consists simply of an additional software.
  • EPIgray available from DOSIsoft, SA, 47 Avenue Carnot, 94230 Cachon, France
  • PerFraction available from Sun Nuclear Corp., 3275 Suntree Blvd, Melbourne, FL 32940.
  • a third solution to determine the fluence uses protocol files, such as LINAC log files generated during delivery, to calculate dose values and distributions received by the patient.
  • protocol files such as LINAC log files generated during delivery
  • Commercial products using log files in this way include Mobius 3D with FX extension (Mobius medical). Persons having ordinary skill in this art debate whether log files can provide the same degree of QA reliability as measurements, as they are generated by the system to be monitored itself.
  • Options of suitable transmission detectors include a two-dimensional (2D) array with parallel readout (i.e. , each pixel connects to an individual channel of the readout electronics). This enables a true measurement of the 2D fluence distribution.
  • a main limitation of such a 2D detector includes that complexity (e.g.
  • a first commercial solution is the DOLPHIN-branded dosimetry detector, available from IBA Dosimetry GmbH (Bruhofstrasse 5, 90592 Schwarzenbruck, Germany). It measures the dose using a 2D ion-chamber pixel array (about 1600 pixels with 5mm pixel pitch projected to isocenter.
  • the Dolphin-brand detector has limitations.
  • the Dolphin-brand detector is limited in terms of attenuation, spatial resolution and patient clearance.
  • ScandiDos Delta4 Discovery brand (available from Scandidos AB, SE-751 83 Uppsala, Sweden) is based on a 2D transmission detector similar to Dolphin brand detector (available from IBA Dosimetry GmbH (Bruhofstrasse 5, 90592 Schwarzenbruck, Germany).
  • the ScandiDos device uses a larger number (about 4000) of smaller sensors (silicon diodes), which are aligned with MLC leaf pairs. For online use, it requires a pre-treatment measurement with a phantom-based array to correlate transmission array measurements and dose.
  • Two-dimensional (2D) arrays with multiplexed readout using large area solid-state technology associates a set of pixels (e.g. a column) to the same readout channel. Complexity increases linearly in relation to the reduction in pixel pitch. However, the cost of solid-state sensors increase more than linearly with area (one pays for the full area, not only the portion covered by pixels). To date there are no commercial applications, although use of MAPS (monolithic active pixel sensors) has been investigated at research level.
  • An integral ion chamber, covering the full beam area, can be used for QA purposes too (see, e.g., Islam US8119978 ).
  • a commercial implementation is the detector IQM (from iRT Systems GmbH, Koblenz, Germany). Since a large parallel plate ion chamber would be suitable to track the integrated fluence only, one of the electrodes is inclined in order to introduce sensitivity to rigid displacements of the MLC aperture, as anticipated by the teaching of Barthelmes ( US4803368 ).
  • the IQM is actually composed of two integral chambers mounted back to back, in order to increase the overall detector sensitivity to errors.
  • the collecting electrode of the integral chamber can be segmented into two interdigitated electrodes with separate readout, in order to obtain with a parallel plate configuration the same error sensitivity which would be given by an inclined electrode, thus simplifying the construction of the device.
  • One-dimensional (1D) arrays of linear sensors can deliver a complete measurement the number of degrees of freedom of the two-dimensional (2D) fluence distribution generated by the one-dimensional motion of the collimator leaves.
  • one sensor can monitor the motion of each leaf-pair.
  • An early commercial example is provided by the DAVID system (PTW Freiburg, Germany), based on an array of wire ion chambers (one wire is associated to each leaf pair). Measurements taken with this detector during a given fraction are compared with a baseline measured during or before the first fraction (which is assumed to be in agreement with the plan), in order to perform a constancy check.
  • 1D linear detectors can do more than performing a constancy check or measuring the integrated fluence.
  • One example is a one-dimensional (1D) scintillating fiber optic array developed by Beaulieu et al. (see for example US Pat. No. 8183534 ,which describes a scintillating optical fiber placed below each leaf pair, and a readout at both extremities using two photodiodes. Fiber optic arrays, however, are difficult to calibrate, prone to radiation damage, and are less suited for beam monitoring.
  • the 1D array measures the position of a leaf instead of a true dose, it can be relatively insensitive to electron contamination and can work with a minimal buildup. In this way, beam perturbation can be minimized.
  • Such an 1D array detector can also estimate the flux (space integral of fluence), but not the true 2D fluence (e.g. one has to assume a given beam profile e.g. from TPS commissioning.
  • the present invention provides an improvement with respect to the known art and solves limitations known in the art.
  • the present invention includes a one-dimensional (1D) strip detector for measuring and determining radiation dose during treatment, when the patient is in position on the couch.
  • the present invention further overcomes limitations of the current state-of-the-art.
  • the present invention provides a measurement system suitable to determine the fluence generated by a radiotherapy LINAC equipped with a MLC with an optimal balance between performance, cost and reliability.
  • the invention is related to a sensing device for a radiation therapy apparatus, the apparatus comprising an accelerator, a beam-shaping device, the beam shaping device being a multi-leaf collimator (MLC) (2) having a plurality of pairs of leaves, and a rotatable gantry, the sensing device comprising :
  • MLC multi-leaf collimator
  • each strip assembly comprises two collecting electrodes mounted back to back
  • the detector comprises two planar bias electrodes which are conductive plates placed on either side of the planar 1D array of strip assemblies, and oriented in such a way that the distance between every strip on one side of the planar 1D array and the corresponding bias electrode increases monotonically along the length of the strip in a given direction along said length, while the distance between every strip on the other side and the corresponding bias electrode decreases monotonically along the length of the strip, in the same given direction.
  • the bias electrodes may be parallel to each other.
  • each strip assembly comprises two adjacent collecting electrode strips, and wherein the width of the first electrode strip decreases monotonically along the length of the strip assembly in a given direction along said length, while the width of the second electrode strip increases monotonically in the same given direction, and wherein the detector further comprises one bias electrode in the form of a conductive plate that is parallel to the planar array of strip assemblies.
  • the sensing device may further comprise a large area collecting electrode on the back side of the planar array of strip assemblies, and a second bias electrode parallel to the planar array of strip assemblies and facing said large area collecting electrode.
  • the invention is related to a sensing device for a radiation therapy apparatus, the apparatus comprising an accelerator, a beam-shaping device the beam shaping device being a multi-leaf collimator (MLC) having a plurality of pairs of leaves, and a rotatable gantry, the sensing device comprising :
  • MLC multi-leaf collimator
  • a sensing device may further comprise a large area collecting electrode on the back side of the planar array of strips, and a second bias electrode parallel to the planar array of strips and facing said large area collecting electrode.
  • the processor may further be configured to enable real-time signalling of a deviation from an expected behavior and a measured behavior exceeding a user-defined tolerance.
  • the sensing device may further comprise :
  • the processor may further be adapted to compare a calculated delivered dose to a planned dose.
  • the processor may further be configured to extract at least one subset of a 3D dose distribution and a set of quality parameters from the calculated 3D dose distribution.
  • the sensing device may further comprise :
  • the processor may further be adapted to accumulate a plurality of data measurements relating to a corresponding plurality of dose delivered to the patient during consecutive fractions.
  • a beam model is a mathematical description of a radiation therapy apparatus in general and contains a number of parameters. These parameters take into account the characteristics of the accelerator (e.g., energy spectrum and lateral beam quality variations), the shapes and positions of the effective radiation sources, and the geometry and material of the beam shaping devices.
  • a fluence is the number of photons/particles passing through a plane perpendicular to the beam per unit time and unit area.
  • a fluence computation algorithm is a set of mathematical rules that computes the fluence according to the beam model and a given parameter set. The representation of the computed fluence (units, coordinate systems) is such that it is compatible with additional computational procedures for computing deposited dose in tissue and/or detector response. Useful descriptions of basic beam modelling techniques are provided, for example, by Wolfgang A.
  • Barthelmes et al. (U.S. Pat. No. 4803368 issued 1989-02-07 ) generally describes an ion chamber having inclined electrodes. Further, Islam et al. (US Pat. No. 8119978 issued on2012-02-21 ) describes the use of strips with variable width in a single ion chamber.
  • the present invention extends the teaching of the prior art. For example, the present invention includes for every leaf pair of an MLC, at least one ion chamber and two readout channels.
  • the sensing device comprises a detector that is built on the basis of a one-dimensional (1D) array of elongate ion chambers.
  • the detector is configured to be attached to the head of a particle accelerator, usually a linear accelerator (LINAC).
  • LINAC linear accelerator
  • the accelerator head is equipped with a multi-leaf collimator and coupled to a rotatable gantry. Therefore, the detector rotates integrally to (meaning together with) the MLC, i.e. the detector is stationary relative to the rotating MLC, which itself can rotate together with the gantry or about the beam axis.
  • the number of ion chambers is chosen in accordance with the number of leaf pairs of the MLC.
  • the ion chambers are defined by elongate conductive strips on a planar surface, preferably on a PCB, said strips acting as the collector electrodes of the ion chambers, and by a single bias electrode facing the strips, if the strips are present on one side of the planar surface, and by two bias electrodes, if the strips are present on both sides.
  • the strips are arranged in strip assemblies, each strip assembly comprising two strips, yielding two respective output currents.
  • An integer number of strip assemblies is associated with each leaf pair when the detector is attached to the accelerator head.
  • a strip assembly is positioned to receive radiation passing through a particular leaf pair. For instance, the motion of a collimator's two banks of leaves, with 80 leaf-pairs may be tracked with a detector having 80 collector electrode strip assemblies, each strip assembly yielding two output currents.
  • the strips connect to 160 corresponding channels of readout electronics.
  • the present invention's number of 160 channels is remarkably lower than the number of pixels (from 1600 to 4000) of these 2D detectors.
  • the strips are arranged as an array of single elongate strips, each strip yielding two output currents at its opposite outer ends.
  • an integer number of strips is associated with each leaf pair when the detector is attached to the accelerator head.
  • Ion chamber technology offers high reliability and performance with moderate manufacturing costs.
  • the embodiments of the present invention may use, or be based upon, collecting electrodes manufactured using standard printed circuit board technology (PCB), which is widely available and economical.
  • PCB printed circuit board technology
  • Figures 1 to 3 illustrate a first embodiment of the present invention.
  • Figure 1 shows a 3D view of the main components of a detector 1 according to the first embodiment, as it is positioned relative to a MLC 2.
  • the first embodiment configures, or is adapted, to measure the position of MLC leaves with a 1D array 10 of collector electrode strip assemblies.
  • the detector comprises, or consists of, two sub-arrays 10a and 10b, which mount back to back.
  • segmented collecting electrodes in the form of printed PCB strips are combined with bias electrodes 11a and 11b in the form of simple, conductive planes placed at non-zero inclination angle relative to the strips.
  • a 'strip assembly' in this embodiment is thus an assembly of two strips mounted back to back.
  • All the strip assemblies are oriented in the direction of the collimator's leaf-movement.
  • the drawing in Figure 1 is a mere sketch, and does not clearly show the number of collector electrode strip assemblies in the 1D array. This number is preferably equal to the number of leaf pairs of the MLC 2, or it may be an integer multiple thereof.
  • the embodiments described hereafter in more detail are based on the case wherein the number of strip assemblies is equal to the number of leaf pairs, and wherein the width of the strips is configured so that when the detector is attached to the accelerator head, each strip assembly in the array is associated with one corresponding leaf pair in the array of leaf pairs of the MLC (i.e. each strip assembly receives the radiation from one leaf pair).
  • Figure 2 illustrates in a side view the position of two leaves 20 of a leaf pair, relative to the detector 1.
  • a collimated beam 21 passes through the opening defined by the leaves 20 and generates two currents i 1 and i 2 respectively, in the strips of the strip assembly that faces the leaf pair (i.e. that is positioned to receive the irradiation passing through said leaf pair).
  • the first current i 1 is generated in the upper IC, between the bias electrode 11a and the upper collector electrode 10a.
  • the second current i 2 is generated in the lower IC, between the bias electrode 11b and the lower electrode 10b.
  • the detector comprises electrometers known as such in the art, which comprise at least a pair of amplifiers 22, for reading out the currents i 1 and i 2 .
  • the bias electrodes 11a, 11b are mutually parallel.
  • the current collected by one side of the n th strip is given by the following proportionality relation (the symbol ⁇ meaning 'proportional to') i n ⁇ ⁇ ⁇ x y n k x ⁇ ydx
  • is the beam photon fluence
  • k(x) is the distance between the strip and the corresponding bias electrode
  • ⁇ y is the strip width
  • y n is the position of the strip center along the array.
  • each IC measurement is the determination of a dose D, applied in the IC volume by the incident radiation.
  • the proportionality constant K is calculated based on the system architecture or, alternatively, experimentally determined with a calibration procedure.
  • the aperture ⁇ x of a leaf pair and its center of gravity x 0 are determined from the currents i 1 and i 2 by the following proportionality relationships ⁇ x ⁇ i 1 + i 2 x 0 ⁇ i 1 ⁇ i 2 / i 1 + i 2
  • K ⁇ V 2 K ⁇ p s ⁇ ⁇ x ⁇ y min + sin ⁇ ⁇ L ⁇ x 0
  • p s is the pitch of the array of collector electrodes (i.e. the pitch of each subarray). This is a constant value for a given detector geometry. p s is equal to the strip width ⁇ y plus the distance between two adjacent strips.
  • the radiation flux is proportional to the sum of the current in all the strips. Knowing the position of the leaves and flux enables one to prove delivery of the prescribed plan (expressed, e.g., in terms of MU/sampling interval, MU/segment or MU/control point), wherein 'MU' stands for 'Monitor Units'.
  • a second embodiment of the present invention includes the 1D array 10 of collector electrode strips shown in Figures 4-6 .
  • each strip assembly 30 consists of a set of adjacent electrodes 31 and 32 printed on the same side of a PCB 33.
  • single bias electrode 34 is provided.
  • the bias electrode 34 is a conductive surface parallel to the collecting electrodes, i.e. at a constant distance h from the array of collecting electrode strips 30.
  • the width of the collecting electrodes 31 and 32 changes monotonically in the x-direction but at an opposite rate. In other words, the first electrode 31 decreases monotonically from high width to low width in the direction +x, while the second electrode 32 increases monotonically from low width to high width in the same direction +x.
  • the detector is again positioned relative to an MLC (not shown) so that an integer number of strip assemblies 30 is associated with each leaf pair.
  • the area below the leaf-pair is shared by two collecting electrodes 31 and 32, which are read-out by independent electrometers (comprising at least two amplifiers 22 coupled respectively to the first and second electrodes 31 and 32), and which produce i 1 and i 2 for each IC strip 30.
  • the geometry for a single collecting electrode strip is shown in detail in Figure 5 , from which the formulas for ⁇ x and x 0 can be derived.
  • the active area of the IC varies while the height h is constant (see Fig. 4 ). Therefore, the current i in the IC is proportional to the area A.
  • each strip 40 of the array 10 is a single collecting electrode on a common PCB 41.
  • the strips are adjacent and isolated from each other, and arranged at a pitch p s .
  • a single bias electrode 42 is placed at a constant height h above the strips 40 (i.e. parallel to the strip array).
  • the strips 40 are read-out at both ends by separate electrometers comprising at least amplifiers 22. This solution works as long as the resistivity of the strip is high enough to limit the parasitic current generated by the voltage offset of the readout amplifiers, and/or to obtain signals with a sufficiently high signal to noise ratio 22.
  • Figure 8 shows the geometry of one strip 40 relative to the bias electrode 42.
  • Equation 15 16 and 17 can be transformed:
  • i 1 i ⁇ i 2 resp .
  • K is a constant containing the sensitivity of the chamber
  • A is the area of the collecting electrode.
  • the broad beam may be described as the sum of an (infinite) number of pencil beams, each of width dx.
  • the currents i 1 and i 2 created by a broad beam are obtained as the continuous sum of all the partial currents d i 1 respectively d i 2 .
  • i 1 and i 2 K L ⁇ x L 2 + x 0
  • a measurement is required of the dose rate ⁇ in order to calculate the value of K (see equation (4)).
  • the dose rate is typically obtained from ionization chamber detectors included in the standard equipment of most irradiation accelerator types.
  • the dose rate is usually measured by a redundant system in the form of two independent ionization chambers. When the sensing device of the invention is mounted on such a LINAC, the dose rate measurement obtained from this redundant system may be used to determine K.
  • the integration over time of ⁇ yields the dose D.
  • a large area ionization chamber may be added on the back side of the array of strips.
  • the large area chamber may be realized by providing a large collecting electrode on the backside of the PCB carrying the strip assemblies 30 (in the case of the second embodiment), or the strips 40 (in the case of the third embodiment).
  • the surface of the large collecting electrode is large enough to receive a beam shaped by any configuration of the MLC.
  • the surface preferably covers an area corresponding to the combined area of the strip assemblies 30 or of the strips 40.
  • a second bias electrode is then provided, on the opposite side of the first bias electrode 34 or 42, preferably at the same distance h from the large collecting electrode, but this distance may also be different from h.
  • An electrometer is further provided for measuring a current output from the large area collector electrode when an electric field is created between the large area collector electrode and the second bias electrode, and when a beam passes through the large area chamber.
  • the large area ionization chamber thus created on the backside of the strips allows to obtain a measurement of the dose rate ⁇ .
  • a sensing device comprises a detector equipped with an array of ion chambers according to any of the above-described embodiments or equivalents thereof.
  • the sensing device comprises a processor configured to determine the position of at least one pair of leaves of the MLC using the currents i 1 and i 2 obtained from the collecting electrode strips.
  • the processor is configured to calculate ⁇ x and x 0 on the basis of predefined formulae which depend on the treatment setup and the detector type, as described above for three particular embodiments. This functionality of the processor is clear as such.
  • the design and practical implementation of the processor may be executed according to known technology.
  • the sensing device further comprises a gantry sensor configured to determine at least one position associated with a first gantry angle,.
  • the gantry sensor detects the angular position of the gantry, allowing the determineddetected parameters of the MLC leaves ( ⁇ x and x 0 ) to be compared to prescribed parameters included in a treatment program for specific gantry angles.
  • An angle sensor is known as such and therefore not described in detail.
  • Leaf aperture and flux are not sufficient alone to determine fluence ⁇ (x,y). If, however, this data is supplemented by independent measurements, e.g. beam profiles measured during LINAC commissioning or periodic QA, then ⁇ (x,y) can be determined. Fluence data is fed into a dose engine in conjunction with patient images (e.g. CT scan) to determine the actual dose delivered in the patient's anatomy. Lastly, this dose distribution is compared to the planned one.
  • patient images e.g. CT scan
  • log files Many LINACs provide a listing of operational parameters during or after the delivery of a treatment (often called log files). These parameters include the evolution of the collimator configuration, the number of delivered MU and the gantry angle. In principle, these parameters should match the requirements of the plan, although deviations are possible due to malfunctioning, or constraints of real operation that were not accounted in the TPS.
  • the machine determines these parameters with a high resolution (e.g. 0.1mm in leaf position), probably higher than what is achievable by the proposed 1D system:
  • log files are not always suitable for QA purposes because the device under exam also generates the data, and a redundant QA process is not possible.
  • One additional contemplated use of the present invention's proposed system is to verify that measurements are consistent with log files (within experimental uncertainty), and to use the log files to determine the fluence and the dose into the patient.

Abstract

A sensing device for a radiation therapy apparatus, the apparatus comprising an accelerator and a beam-shaping device the beam shaping device being a multi-leaf collimator (MLC) (2) having a plurality of pairs of leaves, and a rotatable gantry, the sensing device comprising: €¢ a transmission electronic detector (1) comprising an array of ionization chambers. The ionization chambers are defined by a bias electrode (11a, 11b, 34,42) on the one hand and by a planar array of conductive strips (40) or strip assemblies (30) on the other hand. The strips or strip assemblies are associated to the leaf pairs of the MLC. The strips are the collecting electrodes of the ionization chambers. Each strip assembly or in the case of one particular embodiment, each strip, yields two currents which allow to determine the position of the leaves of a leaf pair associated with the strip or strip assembly, €¢ a gantry sensor configured to determine at least one position associated with a gantry angle; and €¢ a processor, adapted to determine the position of at least one pair of leaves of the MLC using the currents i 1 and i 2 obtained from the collecting electrode strips.

Description

    TECHNICAL FIELD
  • The present invention relates to radiation therapy equipment, and in particular, to devices adapted for enabling quality assurance of external beam radiation therapy for devices using a multi-leaf collimator.
  • BACKGROUND OF THE INVENTION
  • Minimizing the harm to organs at risk, while maximizing the effectiveness of treatment, requires precise determination of actual dose delivered and a comparison to the treatment plan.
  • This treatment plan, prepared following the prescription of the radiation oncologist, requires precise, accurate, and consistent delivery of an optimal radiation dose to the tumor while simultaneously reducing the side effects on healthy surrounding tissues.
  • At present, state-of-the art radiation therapy includes using a multi-leaf collimator to modulate the intensity of the radiation beam produced by a particle beam accelerator, such as a LINAC (conventional radiation-emitting, linear accelerator), to better shape the radiation dose distribution in order to minimize irradiation of non-tumorous cells. A multi-leaf collimator (MLC) is a device made up of individual "leaves" of metal, typically high atomic-number materials such as tungsten. The leaves move independently in and out of the beam aperture in order to partially block the beam. MLCs typically include two banks of opposing leaves. In typical examples, each bank has from approximately 50 to 200 leaves, each leaf having a width in the range from about 0.25 cm to about 2 cm (projected at isocenter). A given leaf in a bank and the opposite one in the other bank form a "leaf pair"; they can move one against the other (thus obstructing the beam) or apart (creating an aperture which, in most installations, can be also asymmetric with respect to the central beam axis [leaf overtravel]). Accordingly, today's treatment plans rely on a properly functioning MLC for Intensity Modulated Radiation Therapy (IMRT). Ideally, verification, in real time, that the MLC is properly operating would allow adaptions to subsequent treatment plans or even an adaptation during delivery
  • The delivered dose may vary from the intended (planned) treatment dose due to a number of factors. These factors include the wrong treatment planning system (TPS) calculation, an error or fault in the physical movement of the LINAC, including its collimating system, or an incorrect beam energy of the LINAC, errors in patient positioning, or changes in the patient's anatomy.
  • Risks of errors is particularly high in the case of plans with high modulation (i.e. a highly conformal dose distribution is obtained with a complex dynamics of the MLC, in such a way that many small beams contribute to the overall dose).
  • The effect of errors on patient's health is then amplified in the case of stereotactic treatments, in which the overall treatment dose is delivered in one or a few fractions (i.e. treatment sessions) instead of using the more traditional approach of many fractions (e.g. 20-30) delivered in subsequent days. Stereotactic treatments are appealing since they increase workflow efficiency, reduce problems related to inter-fractional anatomy changes and, in some cases, present a radiobiological advantage. However, delivering a wrong dose may have severe and irreparable consequences.
  • In traditional pre-treatment quality assurance (QA), a phantom-based detector rests on the couch (replacing the patient) and is irradiated with the treatment plan. This test reveals calculation errors and systematic LINAC errors, and sometimes resulting dose distributions are projected into the patient anatomy with the assumption that it conforms to the CT scan taken at the time of planning. However, this test is not effective against lack of LINAC repeatability, positioning errors, and anatomical changes, as it is done in a special QA fraction and not during patient treatment. Additionally, such pre-treatment QA adds an additional step in the clinical workflow, reducing its efficiency.
  • Radiation sensors suited to quality assurance categorize as either passive or active sensors. Passive sensors cannot alone produce a signal, thus they need to be coupled to a readout system after irradiation. Examples of passive sensors include radiographic and radiochromic films, gels, sensors based on optically- or thermally stimulated luminescence. Active sensors, permanently connected to an electronic readout system, immediately provide a readout signal during irradiation. Examples of active sensors include semiconductor diodes or ion chambers. Active sensors are preferable since they are easier and faster to use, thus increasing the efficiency of the radiotherapy workflow. Moreover, these detectors allow the instantaneous signal readout (not only integrated signal) and are suited to track dynamical changes.
  • Well known examples of active sensors for radiotherapy include diodes and ionization chambers. United States Pat. No. 6,125,335 describes an example of a diode detector having a beam profiler which is in the form of an array of 46 sensor diodes and an off-axis horizontal sensor diode array. Another example of a typical beam profiler is a multi-sensor radiation detector array, such as the Profiler Model 1170 by Sun Nuclear Corporation of Melbourne, Florida, USA.
  • Active sensors, typically organized in an array, provide the spatial resolution necessary to verify the complex radiation patterns generated by an MLC. Known ionization chamber arrays include a device described by Bonin et al., "A pixel chamber to monitor the beam performances in hadron therapy" published in Nuclear Instruments and Methods in Physics Research, No. A 519 (2004) at pages 674-686. This device includes a 2-D array of 1024 successive ionization chambers arranged in a regular matrix of 32×32 pixels.
  • Ionization chambers generally work by providing a high-voltage difference between two electrodes. A gas or liquid fills the space between the electrodes. Typically, this gas is air, nitrogen or an organic gas. Alternatively, an electrically insulating liquid fills this space. The medium gets ionized when exposed to radiation. The ions move towards the electrodes and their charge is collected. An external electrometer measures this charge. The ionization of a molecule of the medium requires a known average energy and, depending on the gas and on the irradiation type, the collected charge is directly proportional to the energy deposited in the gas.
  • United States Pat. No. 4,988,866 describes yet another dosimetry device having a limited number of specifically located sensors for measuring. This device has at least two significant shortcomings. First, due to its configuration and limited number of sensors, it cannot verify the structure of a radiation field of any size. Second, the device can only measure a very narrow energy range due to the single absorber.
  • Still another dosimetry device, described in DE-101 43 609 , improves the spatial resolution of measurement without increasing the number of individual sensors. A set of sensors arranged in a first line and positioned on a rotatable support are located at a different radius from a second line of sensors. The support successively rotates along an angle of 1° or 2°. Successive measurements by lines of sensors taken at angular steps result in a set of many measurement points, thus giving higher spatial resolution. However, this device requires a mechanical drive for the sensor, and the measurement is more time-consuming. In addition, because one takes a set of measurements at successive times, one relies on the stability and constancy of the radiation source. This device cannot measure radiation beam energy.
  • Although some of the devices of the state of the art measure energy using build-up plates of different thicknesses, the user must enter into the treatment room to perform the measurement several times with the required build-up plate for different energies.
  • Yang, Y. and Xing, L. describe a method for determining the position of a leaf of a MLC in "Using the volumetric effect of a finite-sized detector for routine quality assurance of multi leaf collimatorleaf positioning" Med. Phys. 30 (2003) 433-441 . They disclose a finite-size detector, such as an ion chamber, positioned on one leaf as projected in the isocenter plane. Leaf-position displacement increases or decreases the irradiated volume of the detector; therefore, measurement of dose relates to a leaf position error.
  • Islam et al., in U.S. Pat. No. 8119978, issued 2012-02-21 , describes an area integrating fluence-monitoring sensor for measuring radiation dose. The sensor consists of at least one gradient ion chamber having a volume gradient across its length.
  • To detect LINAC errors, which may not be present during first plan delivery, and to increase the efficiency of the clinical workflow, the QA measurement repeats for each treatment fraction while the patient lies on the couch (in-vivo QA). One important benefit is that in-vivo QA enables treatment interruption if something starts to go wrong.
  • Early attempts for in-vivo QA taught taping a few point detectors (usually diodes) to the patient's skin. This verification, however, was limited to few points and measurement uncertainty was high (conversion from the dose measured at skin level and the one in the target depends on many parameters). In addition, these detectors introduced a perturbation to dose distribution. This method requires positioning of the probes in areas of flat dosimetric gradients, which is usually not given for Intensity Modulated Treatments.
  • More recent known devices incorporate either entrance (transmission) or exit (portal) detectors.
  • Other known devices allow the user to evaluate the beam fluence, and thus to project results into the patient's anatomy. Fluence detectors estimate dose perturbation based on detected errors and modify the revealed dose distribution inside the patient accordingly using a dose engine combined with the patient's CT images. Updated images of the patient (e.g. taken with CBCT) in lieu of the planning CT scan can account for positioning errors and any anatomical changes.
  • Various solutions are known for verifying in vivo the dose distribution in a patient. A first solution uses transmission detectors. The detector positions between the LINAC's collimating system and the patient. The collection of measured data (not perturbed by the patient and the couch) is accomplished by sensors with high-dosimetric performance, such as ion chambers.
  • Müller in U.S. Pat. No. 8160204 issued 2012-04-17 , for example, teaches the determination of photon and electron fluence distributions based on measurements. Comparing the expected detector response (calculated from the patient plan on a detector model) to measured data enables a calculation of the difference between actually delivered fluence and the planned fluence.
  • One particularly well suited detector includes the MatriXX brand detector paired with the COMPASS 3D brand software tool, both available from IBA Dosimetry GmbH (Bahnhofstrasse 5, 90592 Schwarzenbruck, Germany). This is a system for 3D treatment plan and delivery verification and patient dose analysis, giving clinically relevant data analysis. The software component of this detector reconstructs dose from measured fluence, compares the patient plan with measurements, and provides 3D dose deposition information inside the patient's anatomy. Accordingly, this arrangement enables visual or statistical plan evaluation (evaluating dose differences/gamma relative to TPS inside patient CT or on a structure-by-structure statistical/quantitative basis via comparison of the TPS generated DVHs to that of this device's independently calculated DVHs).
  • Drawbacks of these known detectors include perturbation of the beam under measurement, reduced space between the patient's head and the couch, and loss of the option of mounting additional accessories at the LINAC. Additionally, radiation scattered from the LINAC's head and collimator, which does not contribute to the dose in the patient, complicates the beam measurement, and modeling with available dose engines and predicting the response becomes difficult.
  • A second solution, a system based on portal imagers, offers a practical solution because most LINAC devices are equipped with a system for portal imaging (typically an amorphous-silicon flat-panel with high spatial resolution) dedicated to control patient positioning. In this example, the QA system consists simply of an additional software. Although there are many commercial solutions available, only few are suited to evaluate the fluence, and calculate the dose in the patient's anatomy. Examples of such commercial systems include EPIgray (available from DOSIsoft, SA, 47 Avenue Carnot, 94230 Cachon, France), and PerFraction (available from Sun Nuclear Corp., 3275 Suntree Blvd, Melbourne, FL 32940). All these solutions suffer of similar limitations including: Portal imagers are optimized for patient imaging and positioning - not dosimetry - and thus require complex dosimetric calibration. Fluence has to be back-projected through the patient, increasing uncertainty; and the intensive QA use may cause radiation damage to the detector. Moreover, they usually cover only a part of the maximum field opening and are therefore not suited for all treatments.
  • A third solution to determine the fluence uses protocol files, such as LINAC log files generated during delivery, to calculate dose values and distributions received by the patient. Commercial products using log files in this way include Mobius 3D with FX extension (Mobius medical). Persons having ordinary skill in this art debate whether log files can provide the same degree of QA reliability as measurements, as they are generated by the system to be monitored itself.. Options of suitable transmission detectors include a two-dimensional (2D) array with parallel readout (i.e. , each pixel connects to an individual channel of the readout electronics). This enables a true measurement of the 2D fluence distribution. A main limitation of such a 2D detector includes that complexity (e.g. the number of readout channels) increases as spatial resolution is increased.For instance, for a square detector of side length L segmented with a uniform distribution of pixels with pitch h - the number of readout channels is (L/h)2. It is thus difficult to reach a resolution suitable to verify stereotactic treatments, characterized by steep dose gradients. However, two commercial solutions are available: A first commercial solution is the DOLPHIN-branded dosimetry detector, available from IBA Dosimetry GmbH (Bahnhofstrasse 5, 90592 Schwarzenbruck, Germany). It measures the dose using a 2D ion-chamber pixel array (about 1600 pixels with 5mm pixel pitch projected to isocenter. The Dolphin-brand detector, however, has limitations. For example, alignment of pixels' rows with MLC leaf pairs is not possible since pixel pitch is larger than typical leaves width. Also, due to constraints of pixel ion-chamber design, the Dolphin-brand detector is limited in terms of attenuation, spatial resolution and patient clearance.
  • A second commercial solution, the ScandiDos Delta4 Discovery brand (available from Scandidos AB, SE-751 83 Uppsala, Sweden) is based on a 2D transmission detector similar to Dolphin brand detector (available from IBA Dosimetry GmbH (Bahnhofstrasse 5, 90592 Schwarzenbruck, Germany). However the ScandiDos device uses a larger number (about 4000) of smaller sensors (silicon diodes), which are aligned with MLC leaf pairs. For online use, it requires a pre-treatment measurement with a phantom-based array to correlate transmission array measurements and dose.
  • Two-dimensional (2D) arrays with multiplexed readout using large area solid-state technology associates a set of pixels (e.g. a column) to the same readout channel. Complexity increases linearly in relation to the reduction in pixel pitch. However, the cost of solid-state sensors increase more than linearly with area (one pays for the full area, not only the portion covered by pixels). To date there are no commercial applications, although use of MAPS (monolithic active pixel sensors) has been investigated at research level.
  • An integral ion chamber, covering the full beam area, can be used for QA purposes too (see, e.g., Islam US8119978 ). A commercial implementation is the detector IQM (from iRT Systems GmbH, Koblenz, Germany). Since a large parallel plate ion chamber would be suitable to track the integrated fluence only, one of the electrodes is inclined in order to introduce sensitivity to rigid displacements of the MLC aperture, as anticipated by the teaching of Barthelmes ( US4803368 ). The IQM is actually composed of two integral chambers mounted back to back, in order to increase the overall detector sensitivity to errors. According to Islam ( WO2016191883 ) the collecting electrode of the integral chamber can be segmented into two interdigitated electrodes with separate readout, in order to obtain with a parallel plate configuration the same error sensitivity which would be given by an inclined electrode, thus simplifying the construction of the device.
  • One-dimensional (1D) arrays of linear sensors can deliver a complete measurement the number of degrees of freedom of the two-dimensional (2D) fluence distribution generated by the one-dimensional motion of the collimator leaves. Thus, one sensor can monitor the motion of each leaf-pair. An early commercial example is provided by the DAVID system (PTW Freiburg, Germany), based on an array of wire ion chambers (one wire is associated to each leaf pair). Measurements taken with this detector during a given fraction are compared with a baseline measured during or before the first fraction (which is assumed to be in agreement with the plan), in order to perform a constancy check.
  • With a more complex arrangement, 1D linear detectors can do more than performing a constancy check or measuring the integrated fluence. In particular, it is possible to determine the position of a pair of leaves and not the true dose. One example is a one-dimensional (1D) scintillating fiber optic array developed by Beaulieu et al. (see for example US Pat. No. 8183534 ,which describes a scintillating optical fiber placed below each leaf pair, and a readout at both extremities using two photodiodes. Fiber optic arrays, however, are difficult to calibrate, prone to radiation damage, and are less suited for beam monitoring.
  • If the 1D array measures the position of a leaf instead of a true dose, it can be relatively insensitive to electron contamination and can work with a minimal buildup. In this way, beam perturbation can be minimized. Such an 1D array detector can also estimate the flux (space integral of fluence), but not the true 2D fluence (e.g. one has to assume a given beam profile e.g. from TPS commissioning.
  • The above considerations yield the rationale for developing an improved detector for online measurement of radiation dose distribution during treatment, when the patient is in position on the couch.
  • SUMMARY OF THE INVENTION
  • The present invention provides an improvement with respect to the known art and solves limitations known in the art. The present invention includes a one-dimensional (1D) strip detector for measuring and determining radiation dose during treatment, when the patient is in position on the couch. The present invention further overcomes limitations of the current state-of-the-art. Specifically, the present invention provides a measurement system suitable to determine the fluence generated by a radiotherapy LINAC equipped with a MLC with an optimal balance between performance, cost and reliability.
  • According to a first aspect, the invention is related to a sensing device for a radiation therapy apparatus, the apparatus comprising an accelerator, a beam-shaping device, the beam shaping device being a multi-leaf collimator (MLC) (2) having a plurality of pairs of leaves, and a rotatable gantry, the sensing device comprising :
    • a transmission electronic detector comprising one or two bias electrodes and a planar 1D array of electrically conductive strip assemblies, and wherein :
      • ▪ each strip assembly comprises two strips which are oriented in the direction of the leaf pairs, and configured so that each leaf pair of the MLC is associated with one or more strip assemblies, when the detector is positioned between the MLC and a target of interest,
      • ▪ each strip together with one of the bias electrodes defines an ion chamber, with the strip forming the collecting electrode thereof,
      • ▪ the two strips of each assembly are configured to respectively produce two currents i1 and i2 when the detector is irradiated by a beam that is shaped by the MLC, and when an electric field is created between the bias electrode(s) and the strips,
      • ▪ said currents allow to determine the position of the pair of MLC leaves associated with the strip assembly,
    • for each strip assembly, a pair of electrometers comprising respective read-out amplifiers coupled respectively to the strips of the assembly, for measuring the currents i1 and i2,
    • a gantry sensor configured to determine at least one position associated with a gantry angle; and
    • a processor, adapted to determine the position of at least one pair of leaves of the MLC using the currents i1 and i2 obtained from the strip assemblies.
  • According to an embodiment of the first aspect of the invention, each strip assembly comprises two collecting electrodes mounted back to back, and wherein the detector comprises two planar bias electrodes which are conductive plates placed on either side of the planar 1D array of strip assemblies, and oriented in such a way that the distance between every strip on one side of the planar 1D array and the corresponding bias electrode increases monotonically along the length of the strip in a given direction along said length, while the distance between every strip on the other side and the corresponding bias electrode decreases monotonically along the length of the strip, in the same given direction. In the latter embodiment, the bias electrodes may be parallel to each other.
  • According to another embodiment of the first aspect of the invention each strip assembly comprises two adjacent collecting electrode strips, and wherein the width of the first electrode strip decreases monotonically along the length of the strip assembly in a given direction along said length, while the width of the second electrode strip increases monotonically in the same given direction, and wherein the detector further comprises one bias electrode in the form of a conductive plate that is parallel to the planar array of strip assemblies.
  • In the latter embodimenth, the sensing device may further comprise a large area collecting electrode on the back side of the planar array of strip assemblies, and a second bias electrode parallel to the planar array of strip assemblies and facing said large area collecting electrode.
  • According a second aspect, the invention is related to a sensing device for a radiation therapy apparatus, the apparatus comprising an accelerator, a beam-shaping device the beam shaping device being a multi-leaf collimator (MLC) having a plurality of pairs of leaves, and a rotatable gantry, the sensing device comprising :
    • a transmission electronic detector comprising a bias electrode and a planar 1D array of electrically conductive strips oriented parallel to the bias electrode, and wherein :
      • ▪ each strip is oriented in the direction of the leaf pairs, and configured so that each leaf pair of the MLC is associated with one or more strips, when the detector is positioned between the MLC and a target of interest,
      • ▪ each strip together with one of the bias electrodes defines an ion chamber, with the strip forming the collecting electrode thereof,
      • ▪ each strip is configured to respectively produce two currents i1 and i2 at the opposite outer ends of the strip, when the detector is irradiated by a beam that is shaped by the MLC, and when an electric field is created between the bias electrode(s) and the strips,
      • ▪ said currents allow to determine the position of the pair of MLC leaves associated with the strip assembly,
    • for each strip, a pair of electrometers comprising respective read-out amplifiers coupled respectively to the opposite outer ends of the strip, for measuring the currents i1 and i2,
    • a gantry sensor configured to determine at least one position associated with a gantry angle; and
    • a processor, adapted to determine the position of at least one pair of leaves of the MLC using the currents i1 and i2 obtained from the strips.
  • A sensing device according to the second aspect may further comprise a large area collecting electrode on the back side of the planar array of strips, and a second bias electrode parallel to the planar array of strips and facing said large area collecting electrode.
  • According to further embodiments of the invention according to the first or second aspect:
    • the sensing device further comprises a receiver, configured to receive a patient plan (including patient dose, patient images taken before planning and radiation beam configuration); and
    • the processor is adapted to compare expected (from plan) and measured LINAC behavior (in term of leaves position, flux and gantry angle).
  • The processor may further be configured to enable real-time signalling of a deviation from an expected behavior and a measured behavior exceeding a user-defined tolerance.
  • The sensing device may further comprise :
    • a storage means in communication with the processor, the storage means configured to store a beam model of the accelerator;
    • a fluence-reconstructing engine configured to compute a photon fluence based on the position of pairs of the leaves determined by the processor and on the beam model wherein the photon fluence is determined independently from a presence of a patient or phantom in the beam; and
    • a dose calculation engine configured to compute a delivered 3D dose distribution using fluence and patient images.
  • The processor may further be adapted to compare a calculated delivered dose to a planned dose. The processor may further be configured to extract at least one subset of a 3D dose distribution and a set of quality parameters from the calculated 3D dose distribution.
  • The sensing device may further comprise :
    • an image acquisition station configured to acquire updated patient images;
    • an image processor configured to communicate the images to the dose calculation engine whereby the dose calculation engine calculates a dose distribution based on measured fluence and the updated patient images.
  • The processor may further be adapted to accumulate a plurality of data measurements relating to a corresponding plurality of dose delivered to the patient during consecutive fractions.
  • A further summary of characteristic features of the invention is given hereafter.
    1. A. For a radiation therapy apparatus comprising a beam-shaping device comprising an accelerator, a multi-leaf collimator (MLC) and a rotatable gantry, a sensing device comprising:
      1. (a) a transmission electronic detector comprising a 1D array comprising at least two ion chambers, whereby the transmission electronic detector disposes between the beam shaping device and a target of interest, wherein at least one ion chamber of the 1D array is associated with at least one pair of leaves of the MLC;
      2. (b) a gantry sensor configured to determine at least one position associated with a first gantry angle; and
      3. (c) a processor, adapted to determine the position of the at least one pair of leaves of the MLC using data obtained from the at least two ion chambers of the 1D array.
    2. B. The sensing device of paragraph A wherein:
      1. (a) a receiver is configured to receive a patient plan (including patient dose, patient images taken before planning, and radiation beam configuration); and
      2. (b) a processor adapted to compare expected (from plan) and measured LINAC behavior (in term of leaves position, flux and gantry angle).
      The sensing device of paragraph B wherein:
      the processor is further configured to enable real-time signaling a deviation from an expected behavior and a measured behavior exceeding a user-defined tolerance.
    3. C. The sensing device of paragraph B, further comprising:
      1. (a) a storage means in communication with the processor, the storage means configured to store a beam model of the accelerator;
      2. (b) a fluence-reconstructing engine configured to compute a photon fluence based on the position of pairs of the leaves determined by the processor and on the beam model wherein the photon fluence is determined independently from a presence of a patient or phantom in the beam; and
      3. (c) a dose calculation engine configured to compute a delivered 3D dose distribution using fluence and patient images.
    4. D. The sensing device of paragraph C wherein:
      the processor is further adapted to compare a calculated delivered dose to a planned dose.
    5. E. The sensing device of paragraph C wherein:
      the processor is further configured to extract at least one subset of a 3D dose distribution and a set of quality parameters from the calculated 3D dose distribution.
    6. F. The sensing device of paragraph C further comprising:
      1. a. an acquirer configured to acquire updated patient images;
      2. b. an image processor configured to communicate the images to the dose engine whereby the dose engine calculates a dose distribution based on measured fluence and the updated patient images.
    7. G. The sensing device according to paragraph C or paragraph F wherein:
      the processor is further adapted to accumulate a plurality of data measurements relating to a corresponding plurality of dose delivered to the patient during consecutive fractions.
    8. H. The sensing device according to any of previous paragraphs further characterized by:
      the 1D array further comprises two sub arrays of ion chambers mounted back to back and orthogonal to the beam axis, in which:
      1. (a) each collecting electrode comprises at least two conductive strips oriented in the direction of MLC pairs, lying on a plane surface;
      2. (b) an ionization current collected by the conductive strips is measured by an electrometer connected electrically to the respective strip;
      3. (c) at least two bias electrodes comprise a corresponding pair of conductive plates oriented in such a way that the distance between every strip and the bias electrode increases monotonically along the length of the strip;
      4. (d) the two bias electrodes are substantially parallel; and
      5. (e) at least one strip from each sub array is associated to a corresponding MLC leaf.
    9. I. The device according to any of the paragraphs A-H, in which:
      1. (a) an area below each corresponding leaf is shared by two collecting electrodes wherein the electrode comprises an individual thickness that changes monotonically when moving from one side of the leaf to the other; the thickness is measured in the plane of the electrodes. In the remainder of this specification, it will therefore be referred to as the width of the electrodes.
      2. (b) an ionization current is collected by each electrode and the ionization current is measured by an independent electrometer electrically connected to the electrode; and
      3. (c) a distance between the collecting electrodes and the bias electrodes is uniform throughout the array.
    10. J. The device according to any of the paragraphs A-I, in which the 1D array is characterized by an array of ion chambers in which:
      1. (a) collecting electrodes comprise a conductive strip oriented in the direction of MLC pairs, lying on a plane surface; the 'conductive' strip will be referred to in the remainder of this specification as 'resistive' strip.
      2. (b) a bias electrode is a plate with uniform distance from the strip; and
      3. (c) an ionization current collected by each strip is measured via two independent electrometers connected to the opposite extremities of the strip.
    BRIEF DESCRIPTION OF THE DRAWINGS
    • Figure 1 is a 3D view of the main components of a detector in accordance with a first embodiment of the invention, and its position relative to a multi-leaf collimator (MLC).
    • Figure 2 is a side view of the detector in accordance with the first embodiment.
    • Figure 3 is another side view of the collimator in accordance with the first embodiment, including indications of the geometrical parameters needed to determine the position of two leaves of the MLC.
    • Figure 4 is a 3D view of the main components of a detector in accordance with a second embodiment of the invention.
    • Figure 5 is a plane view of one collecting electrode strip in a detector according to the second embodiment.
    • Figure 6 illustrates a number of geometrical considerations made in relation to a detector according to the second embodiment.
    • Figure 7 is a 3D view of the main components of a detector in accordance with a third embodiment of the invention.
    • Figure 8 is a side view of one collecting electrode strip of the detector according to the third embodiment, including a number of geometrical parameters needed to calculate the position of two leaves of the MLC.
    • Figure 9 illustrates the electrical equivalent of one collecting electrode strip in a detector according to the third embodiment.
    DETAILED DESCRIPTION OF THE INVENTION
  • This detailed description sets forth specific details of exemplary embodiments of the present invention to enable a thorough understanding: However, those of ordinary skill in the art understand that the scope of the present invention is not restricted to those embodiments. Other embodiments, whether described herein or omitted, may be practiced without these specific details. In other instances, well-known methods, procedures and components are not described in detail, thus avoiding obscuring the scope of the invention. Accordingly, this description does not limit the invention to the scope of the embodiments described. Although certain features are illustrated and described, many modifications, substitutions, changes, and equivalents will now occur to those of ordinary skill in the art. It is understood that the appended claims are intended to cover such modifications and changes as fall within the true spirit of the invention.
  • A beam model is a mathematical description of a radiation therapy apparatus in general and contains a number of parameters. These parameters take into account the characteristics of the accelerator (e.g., energy spectrum and lateral beam quality variations), the shapes and positions of the effective radiation sources, and the geometry and material of the beam shaping devices. A fluence is the number of photons/particles passing through a plane perpendicular to the beam per unit time and unit area. A fluence computation algorithm is a set of mathematical rules that computes the fluence according to the beam model and a given parameter set. The representation of the computed fluence (units, coordinate systems) is such that it is compatible with additional computational procedures for computing deposited dose in tissue and/or detector response. Useful descriptions of basic beam modelling techniques are provided, for example, by Wolfgang A. Tome, "Beam Modelling for a Convolution/Superposition-Based Treatment Planning System", Medical Dosimetry, Vol. 27, No. 1, pp. 11-19, 2002; or by Nikos Papanikolaou, "Investigation of the convolution method for polyenergetic spectra", Med. Phys. 20(5), 1993.
  • Three technical solutions for the sensor have been identified, which allow the measurement of the aperture of each leaf pair AND the center of gravity of the aperture (it is not simply the measurement of a linear dose integral). These solutions are:
    1. a) A double array of regular strips with inclined electrodes (Fig. 1-3)
    2. b) Strips with variable width (Fig. 4-6)
    3. c) Resistive strips (Fig. 7-9).
  • Barthelmes et al. (U.S. Pat. No. 4803368 issued 1989-02-07 ) generally describes an ion chamber having inclined electrodes. Further, Islam et al. (US Pat. No. 8119978 issued on2012-02-21 ) describes the use of strips with variable width in a single ion chamber. However, the present invention extends the teaching of the prior art. For example, the present invention includes for every leaf pair of an MLC, at least one ion chamber and two readout channels.
  • The sensing device according to the present invention comprises a detector that is built on the basis of a one-dimensional (1D) array of elongate ion chambers. The detector is configured to be attached to the head of a particle accelerator, usually a linear accelerator (LINAC). The accelerator head is equipped with a multi-leaf collimator and coupled to a rotatable gantry. Therefore, the detector rotates integrally to (meaning together with) the MLC, i.e. the detector is stationary relative to the rotating MLC, which itself can rotate together with the gantry or about the beam axis. There is a direct mapping between the ion chambers and the leaf-pairs of the MLC. This means that the number of ion chambers is chosen in accordance with the number of leaf pairs of the MLC. The ion chambers are defined by elongate conductive strips on a planar surface, preferably on a PCB, said strips acting as the collector electrodes of the ion chambers, and by a single bias electrode facing the strips, if the strips are present on one side of the planar surface, and by two bias electrodes, if the strips are present on both sides. According to a first and second embodiment, the strips are arranged in strip assemblies, each strip assembly comprising two strips, yielding two respective output currents. An integer number of strip assemblies is associated with each leaf pair when the detector is attached to the accelerator head. 'Associated with' means that a strip assembly is positioned to receive radiation passing through a particular leaf pair. For instance, the motion of a collimator's two banks of leaves, with 80 leaf-pairs may be tracked with a detector having 80 collector electrode strip assemblies, each strip assembly yielding two output currents. The strips connect to 160 corresponding channels of readout electronics. In contrast to conventional 2D arrays presently commercialized as online dosimeters, the present invention's number of 160 channels is remarkably lower than the number of pixels (from 1600 to 4000) of these 2D detectors.
  • According to a third embodiment, the strips are arranged as an array of single elongate strips, each strip yielding two output currents at its opposite outer ends. In the third embodiment therefore, an integer number of strips is associated with each leaf pair when the detector is attached to the accelerator head.
  • Ion chamber technology offers high reliability and performance with moderate manufacturing costs. The embodiments of the present invention may use, or be based upon, collecting electrodes manufactured using standard printed circuit board technology (PCB), which is widely available and economical.
  • Figures 1 to 3 illustrate a first embodiment of the present invention. Figure 1 shows a 3D view of the main components of a detector 1 according to the first embodiment, as it is positioned relative to a MLC 2. The first embodiment configures, or is adapted, to measure the position of MLC leaves with a 1D array 10 of collector electrode strip assemblies. The detector comprises, or consists of, two sub-arrays 10a and 10b, which mount back to back. In each sub array, segmented collecting electrodes, in the form of printed PCB strips are combined with bias electrodes 11a and 11b in the form of simple, conductive planes placed at non-zero inclination angle relative to the strips. A 'strip assembly' in this embodiment is thus an assembly of two strips mounted back to back. All the strip assemblies are oriented in the direction of the collimator's leaf-movement. The drawing in Figure 1 is a mere sketch, and does not clearly show the number of collector electrode strip assemblies in the 1D array. This number is preferably equal to the number of leaf pairs of the MLC 2, or it may be an integer multiple thereof. The embodiments described hereafter in more detail are based on the case wherein the number of strip assemblies is equal to the number of leaf pairs, and wherein the width of the strips is configured so that when the detector is attached to the accelerator head, each strip assembly in the array is associated with one corresponding leaf pair in the array of leaf pairs of the MLC (i.e. each strip assembly receives the radiation from one leaf pair). Figure 2 illustrates in a side view the position of two leaves 20 of a leaf pair, relative to the detector 1. A collimated beam 21 passes through the opening defined by the leaves 20 and generates two currents i1 and i2 respectively, in the strips of the strip assembly that faces the leaf pair (i.e. that is positioned to receive the irradiation passing through said leaf pair). The first current i1 is generated in the upper IC, between the bias electrode 11a and the upper collector electrode 10a. The second current i2 is generated in the lower IC, between the bias electrode 11b and the lower electrode 10b. The detector comprises electrometers known as such in the art, which comprise at least a pair of amplifiers 22, for reading out the currents i1 and i2. According to a preferred embodiment, the bias electrodes 11a, 11b are mutually parallel.
  • Assuming that scattering and lateral spread of secondary electrons in the chamber volume can be neglected, the current collected by one side of the nth strip is given by the following proportionality relation (the symbol ∝ meaning 'proportional to') i n Φ x y n k x Δ ydx
    Figure imgb0001
  • Where Φ is the beam photon fluence, k(x) is the distance between the strip and the corresponding bias electrode, Δy is the strip width and yn is the position of the strip center along the array. For the detectors described hereafter, the above equation can be written in a simplified way.
  • The aim of each IC measurement is the determination of a dose D, applied in the IC volume by the incident radiation. "Dose" is the energy deposited by radiation in an amount of mass of material: D = E m = E V ρ
    Figure imgb0002
    where
    • D= dose in Gy (1Gy =1J/kg)
    • E = energy
    • m= mass
    • V= volume in m3 or mm3
    • ρ= density in kg/m3
  • What is measured however, is the charge Q produced in the sensitive volume. Q = E W air e = D V ρ e W air
    Figure imgb0003
    where W air is the ionization constant for air (the energy which is required to produce one pair of charges in air).
  • The measured current i is the time derivative of the charge Q collected: i = D ˙ V ρ e W air = K V with K = D ˙ ρ e W air
    Figure imgb0004
    where is the dose rate (in Gy/min or Gy/s). depends on the treatment machine emitting the ionizing radiation and its calibration. Air density is approximated with 1.2kg/m3. W air and e are physical constants with the values W air=33.97eV respectively e = 1.602 · 10-19C. Here and in the following formulae, the proportionality constant K is calculated based on the system architecture or, alternatively, experimentally determined with a calibration procedure.
  • With the assumption that the fluence is uniform within the collimator aperture, the aperture Δx of a leaf pair and its center of gravity x0, are determined from the currents i1 and i2 by the following proportionality relationships Δx i 1 + i 2 x 0 i 1 i 2 / i 1 + i 2
    Figure imgb0005
  • With reference to Figure 3, the above relationships can be derived for the detector according to the first embodiment : i 1 = K V 1 = K p s Δ x y min + sin α x 0
    Figure imgb0006
    i 2 = K V 2 = K p s Δ x y min + sin α L x 0
    Figure imgb0007
  • Herein, p s is the pitch of the array of collector electrodes (i.e. the pitch of each subarray). This is a constant value for a given detector geometry. p s is equal to the strip width Δy plus the distance between two adjacent strips.
  • Δx and x0 can be obtained from building the sum i 1+i 2 and the difference i 1-i 2 and result in: Δ x = i 1 + i 2 K p s 2 y min + sin α L = i 1 + i 2
    Figure imgb0008
    x 0 = L 2 + i 1 i 2 2 K p s Δ x sin α
    Figure imgb0009
  • Additionally, the radiation flux is proportional to the sum of the current in all the strips. Knowing the position of the leaves and flux enables one to prove delivery of the prescribed plan (expressed, e.g., in terms of MU/sampling interval, MU/segment or MU/control point), wherein 'MU' stands for 'Monitor Units'.
  • Those having ordinary skill in the art will appreciate that other solutions, alternative to the above-described inclined electrode ion chambers, can easily be figured out, leading to similar relationship between leave positions and measured current. The key point is that the sensitivity of each strip has to be modulated along its length to give sensitivity to beam position x0.
  • A second embodiment of the present invention includes the 1D array 10 of collector electrode strips shown in Figures 4-6. Here each strip assembly 30 consists of a set of adjacent electrodes 31 and 32 printed on the same side of a PCB 33. single bias electrode 34 is provided. The bias electrode 34 is a conductive surface parallel to the collecting electrodes, i.e. at a constant distance h from the array of collecting electrode strips 30. The width of the collecting electrodes 31 and 32 changes monotonically in the x-direction but at an opposite rate. In other words, the first electrode 31 decreases monotonically from high width to low width in the direction +x, while the second electrode 32 increases monotonically from low width to high width in the same direction +x. According to the preferred embodiment, for a given position x, the sum of the width of both strips at this position is constant. The width of the electrodes thereby modulates sensitivity along the length. The detector is again positioned relative to an MLC (not shown) so that an integer number of strip assemblies 30 is associated with each leaf pair. In detail, when each strip assembly 30 is associated with only one leaf pair, the area below the leaf-pair is shared by two collecting electrodes 31 and 32, which are read-out by independent electrometers (comprising at least two amplifiers 22 coupled respectively to the first and second electrodes 31 and 32), and which produce i1 and i2 for each IC strip 30. The geometry for a single collecting electrode strip is shown in detail in Figure 5, from which the formulas for Δx and x0 can be derived.
  • In case of the triangular strips design shown in Figure 5, the active area of the IC varies while the height h is constant (see Fig. 4). Therefore, the current i in the IC is proportional to the area A. The areas A 1 and A 2 (see Figure 5) can be expressed using geometric parameters of the strip as well as the parameters Δx and x0 as follows: i 1 = K A 1 = K Δ x z = K Δ x Z L x 0
    Figure imgb0010
    i 2 = K A 2 = K Δ x Z z = K Δ x Z 1 x 0 L
    Figure imgb0011
  • To obtain these equations, some geometric considerations shown in Figure 6 were made. On the left hand side, it is illustrated that the respective area Ai can be calculated like a square area, since A = A'. Additionally, the following relationship between z, Z, L, and x 0 (illustrated on the right hand side) is applied: Z L = z x 0 Z = Z L x 0 .
    Figure imgb0012
  • From the two currents i1 and i2, Δx and x 0 can be determined as follows Δ x = i 1 + i 2 K Z
    Figure imgb0013
    x 0 = L 2 i 1 i 2 i 1 + i 2
    Figure imgb0014
  • Figures 7 to 9 illustrate the third embodiment of the present invention. Here, each strip 40 of the array 10 is a single collecting electrode on a common PCB 41. The strips are adjacent and isolated from each other, and arranged at a pitch ps. A single bias electrode 42 is placed at a constant height h above the strips 40 (i.e. parallel to the strip array). The strips 40 are read-out at both ends by separate electrometers comprising at least amplifiers 22. This solution works as long as the resistivity of the strip is high enough to limit the parasitic current generated by the voltage offset of the readout amplifiers, and/or to obtain signals with a sufficiently high signal to noise ratio 22. Figure 8 shows the geometry of one strip 40 relative to the bias electrode 42. The resistance R of the strip 40 depends on the length L (or rather x in the general case as indicated in Figure 8 for the length coordinate) and resistance per unit length λ of the resistive material: λ = R L and R = L λ
    Figure imgb0015
  • When a pencil beam creates free charges in the region above the resistive strip at location x, these currents are collected on the strip and flow towards the electrometer amplifiers through resistors R 1 and R 2 formed by the resistive material. The current causes a voltage drop over those resistors on both sides of the intersection of the pencil beam with the strip. As the inputs to the electrometer amplifiers are at ground potential, this can be described as a parallel circuit with two branches in which a division of currents into partial currents i 1 and i 2 takes place. This principle is visualized in Figure 9. The sum of i 1 and i 2 always equals the total current.
  • For the circuit as sketched in Figure 9, the following general equations can be defined: i = i 1 + i 2
    Figure imgb0016
    u 1 = u 2 or R 1 i 1 = R 2 i 2
    Figure imgb0017
    with R 1 = λ x L 2 and R 2 = λ L 2 x
    Figure imgb0018
  • In order to obtain i 1 and i 2, Equation 15, 16 and 17 can be transformed:
    By applying i 1 = i i 2 resp . i 2 = i i 1
    Figure imgb0019
    and
    R = R 1 + R 2 total resistance = constant
    Figure imgb0020
    the currents i 1 and i 2 can be described as: i 1 = i R 2 R and i 2 = i R 1 R .
    Figure imgb0021
  • As explained above, the current in an IC is proportional to the volume or area of the chamber (depending on the IC design). Therefore, the derivative must be proportional to the derivative of the volume or area, too: d i = K d A
    Figure imgb0022
    where K is a constant containing the sensitivity of the chamber and A is the area of the collecting electrode. In case of the resistive strip design, two dimensions of the chamber volume remain constant: the strip pitch p s as well as the chamber height h (distance between the strip array and the bias electrode, see Figure 8). Therefore, only one dimension of the volume V varies, namely, its length, x. di can be determined as d i = K d x .
    Figure imgb0023
  • The broad beam may be described as the sum of an (infinite) number of pencil beams, each of width dx. The current di (comprised of di 1 and di 2) is created by the charges formed in the area which is hit by the pencil beam of width dx. Therefore, applying Equations 18 and 19, the currents di 1 resp. di 2 induced in the strip can be written as d i 1 = R 2 R d i and d i 2 = R 1 R d i .
    Figure imgb0024
  • With use of Equation 17, it can be concluded that d i 1 = R 2 R d i = λ L 2 x R K d x
    Figure imgb0025
    and analogously d i 2 = R 1 R d i = λ x L 2 R K d x .
    Figure imgb0026
  • The currents i 1 and i 2 created by a broad beam, are obtained as the continuous sum of all the partial currents di 1 respectively di 2. Exemplarily, calculations for i 1 are presented here; calculations for i 2 were performed analogously. i 1 = x 0 Δ x 2 x 0 + Δ x 2 λ L 2 x R K dx = λK 2 R L 2 X | x 0 Δ x 2 2 x 0 + Δ x 2 = λK 2 R L 2 x 0 + Δ x 2 2 L 2 x 0 Δ x 2 2 = λK 2 R 2 Δ x L 2 x 0 = K L Δ x L 2 x 0
    Figure imgb0027
  • After integration, i 1 and i 2 can be written as: i 1 = K L Δ x L 2 x 0 and
    Figure imgb0028
    i 2 = K L Δ x L 2 + x 0
    Figure imgb0029
  • From these two equations Δx and x 0 can be determined: Δ x = i 1 + i 2 K
    Figure imgb0030
    x 0 = L 2 i 1 i 2 i 1 + i 2
    Figure imgb0031
  • In any of the embodiments described above, as in any ionization chamber measurement, a measurement is required of the dose rate in order to calculate the value of K (see equation (4)). The dose rate is typically obtained from ionization chamber detectors included in the standard equipment of most irradiation accelerator types. In a LINAC the dose rate is usually measured by a redundant system in the form of two independent ionization chambers. When the sensing device of the invention is mounted on such a LINAC, the dose rate measurement obtained from this redundant system may be used to determine K. The integration over time of yields the dose D.
  • In special forms of the second and third embodiment, a large area ionization chamber may be added on the back side of the array of strips. The large area chamber may be realized by providing a large collecting electrode on the backside of the PCB carrying the strip assemblies 30 (in the case of the second embodiment), or the strips 40 (in the case of the third embodiment). The surface of the large collecting electrode is large enough to receive a beam shaped by any configuration of the MLC. The surface preferably covers an area corresponding to the combined area of the strip assemblies 30 or of the strips 40. A second bias electrode is then provided, on the opposite side of the first bias electrode 34 or 42, preferably at the same distance h from the large collecting electrode, but this distance may also be different from h. An electrometer is further provided for measuring a current output from the large area collector electrode when an electric field is created between the large area collector electrode and the second bias electrode, and when a beam passes through the large area chamber. The large area ionization chamber thus created on the backside of the strips allows to obtain a measurement of the dose rate .
  • A sensing device according to the invention comprises a detector equipped with an array of ion chambers according to any of the above-described embodiments or equivalents thereof. In addition, the sensing device comprises a processor configured to determine the position of at least one pair of leaves of the MLC using the currents i1 and i2 obtained from the collecting electrode strips. In other words, the processor is configured to calculate Δx and x0 on the basis of predefined formulae which depend on the treatment setup and the detector type, as described above for three particular embodiments. This functionality of the processor is clear as such. The design and practical implementation of the processor may be executed according to known technology. The sensing device further comprises a gantry sensor configured to determine at least one position associated with a first gantry angle,. In other words, the gantry sensor detects the angular position of the gantry, allowing the determineddetected parameters of the MLC leaves (Δx and x0) to be compared to prescribed parameters included in a treatment program for specific gantry angles. An angle sensor is known as such and therefore not described in detail.
  • The procedures outlined above represent a first approximation to obtain the position of the leaves from measured current. This approach is acceptable if scattering and lateral spread of electrons are negligible. If this is not the case, a more refined calculation, using the signal from neighbor strips, is performed. In this case Δx and x0 are calculated using the signal from neighbor strips too.
  • Leaf aperture and flux are not sufficient alone to determine fluence Φ(x,y). If, however, this data is supplemented by independent measurements, e.g. beam profiles measured during LINAC commissioning or periodic QA, then Φ(x,y) can be determined. Fluence data is fed into a dose engine in conjunction with patient images (e.g. CT scan) to determine the actual dose delivered in the patient's anatomy. Lastly, this dose distribution is compared to the planned one.
  • Many LINACs provide a listing of operational parameters during or after the delivery of a treatment (often called log files). These parameters include the evolution of the collimator configuration, the number of delivered MU and the gantry angle. In principle, these parameters should match the requirements of the plan, although deviations are possible due to malfunctioning, or constraints of real operation that were not accounted in the TPS. The machine determines these parameters with a high resolution (e.g. 0.1mm in leaf position), probably higher than what is achievable by the proposed 1D system: However, log files are not always suitable for QA purposes because the device under exam also generates the data, and a redundant QA process is not possible.
  • One additional contemplated use of the present invention's proposed system is to verify that measurements are consistent with log files (within experimental uncertainty), and to use the log files to determine the fluence and the dose into the patient.

Claims (14)

  1. A sensing device for a radiation therapy apparatus, the apparatus comprising an accelerator, a beam-shaping device, the beam shaping device being a multi-leaf collimator (MLC) (2) having a plurality of pairs of leaves, and a rotatable gantry, the sensing device comprising :
    • a transmission electronic detector (1) comprising one or two bias electrodes (11a,11b,34) and a planar 1D array (10) of electrically conductive strip assemblies, and wherein :
    ▪ each strip assembly comprises two strips which are oriented in the direction of the leaf pairs, and configured so that each leaf pair of the MLC (2) is associated with one or more strip assemblies, when the detector is positioned between the MLC and a target of interest,
    ▪ each strip together with one of the bias electrodes defines an ion chamber, with the strip forming the collecting electrode thereof,
    ▪ the two strips of each assembly are configured to respectively produce two currents i1 and i2 when the detector is irradiated by a beam that is shaped by the MLC, and when an electric field is created between the bias electrode(s) and the strips,
    ▪ said currents allow to determine the position of the pair of MLC leaves associated with the strip assembly,
    • for each strip assembly, a pair of electrometers comprising respective read-out amplifiers (22) coupled respectively to the strips of the assembly, for measuring the currents i1 and i2,
    • a gantry sensor configured to determine at least one position associated with a gantry angle; and
    • a processor, adapted to determine the position of at least one pair of leaves of the MLC using the currents i1 and i2 obtained from the strip assemblies.
  2. The sensing device according to claim 1, wherein each strip assembly comprises two collecting electrodes (10a,10b) mounted back to back, and wherein the detector comprises two planar bias electrodes (11 a,11 b) which are conductive plates placed on either side of the planar 1D array of strip assemblies, and oriented in such a way that the distance between every strip on one side of the planar 1D array and the corresponding bias electrode increases monotonically along the length of the strip in a given direction along said length, while the distance between every strip on the other side and the corresponding bias electrode decreases monotonically along the length of the strip, in the same given direction.
  3. The sensing device according to claim 2, wherein the bias electrodes (11a,11b) are parallel to each other.
  4. The sensing device according to claim 1, wherein each strip assembly (30) comprises two adjacent collecting electrode strips (31,32), and wherein the width of the first electrode strip (31) decreases monotonically along the length of the strip assembly in a given direction along said length, while the width of the second electrode strip (32) increases monotonically in the same given direction, and wherein the detector further comprises one bias electrode (34) in the form of a conductive plate that is parallel to the planar array (10) of strip assemblies (30).
  5. The sensing device according to claim 4, further comprising a large area collecting electrode on the back side of the planar array of strip assemblies (30), and a second bias electrode parallel to the planar array of strip assemblies and facing said large area collecting electrode.
  6. A sensing device for a radiation therapy apparatus, the apparatus comprising an accelerator, a beam-shaping device the beam shaping device being a multi-leaf collimator (MLC) (2) having a plurality of pairs of leaves, and a rotatable gantry, the sensing device comprising :
    • a transmission electronic detector (1) comprising a bias electrode (42) and a planar 1D array (10) of electrically conductive strips (40) oriented parallel to the bias electrode, and wherein :
    ▪ each strip is oriented in the direction of the leaf pairs, and configured so that each leaf pair of the MLC (2) is associated with one or more strips, when the detector is positioned between the MLC and a target of interest,
    ▪ each strip together with one of the bias electrodes defines an ion chamber, with the strip forming the collecting electrode thereof,
    ▪ each strip is configured to respectively produce two currents i1 and i2 at the opposite outer ends of the strip, when the detector is irradiated by a beam that is shaped by the MLC, and when an electric field is created between the bias electrode(s) and the strips,
    ▪ said currents allow to determine the position of the pair of MLC leaves associated with the strip assembly,
    • for each strip, a pair of electrometers comprising respective read-out amplifiers (22) coupled respectively to the opposite outer ends of the strip, for measuring the currents i1 and i2,
    • a gantry sensor configured to determine at least one position associated with a gantry angle; and
    • a processor, adapted to determine the position of at least one pair of leaves of the MLC using the currents i1 and i2 obtained from the strips.
  7. The sensing device according to claim 6, further comprising a large area collecting electrode on the back side of the planar array of strips (40), and a second bias electrode parallel to the planar array of strips and facing said large area collecting electrode.
  8. The sensing device of any one of the preceding claims wherein
    • the sensing device further comprises a receiver, configured to receive a patient plan (including patient dose, patient images taken before planning and radiation beam configuration); and
    • the processor is adapted to compare expected (from plan) and measured LINAC behavior (in term of leaves position, flux and gantry angle).
  9. The sensing device of claim 8 wherein the processor is further configured to enable real-time signalling of a deviation from an expected behavior and a measured behavior exceeding a user-defined tolerance.
  10. The sensing device of claim 8, further comprising:
    • a storage means in communication with the processor, the storage means configured to store a beam model of the accelerator;
    • a fluence-reconstructing engine configured to compute a photon fluence based on the position of pairs of the leaves determined by the processor and on the beam model wherein the photon fluence is determined independently from a presence of a patient or phantom in the beam; and
    • a dose calculation engine configured to compute a delivered 3D dose distribution using fluence and patient images.
  11. The sensing device of claim 9 wherein the processor is further adapted to compare a calculated delivered dose to a planned dose.
  12. The sensing device of claim 9 wherein:
    the processor is further configured to extract at least one subset of a 3D dose distribution and a set of quality parameters from the calculated 3D dose distribution.
  13. The sensing device of claim 9 further comprising:
    • an image acquisition station configured to acquire updated patient images;
    • an image processor configured to communicate the images to the dose calculation engine whereby the dose calculation engine calculates a dose distribution based on measured fluence and the updated patient images.
  14. The sensing device according to claim 9 or claim 12 wherein the processor is further adapted to accumulate a plurality of data measurements relating to a corresponding plurality of dose delivered to the patient during consecutive fractions.
EP18195642.6A 2017-09-20 2018-09-20 One dimensional transmission detector for radiotherapy Active EP3460530B1 (en)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
US201762561152P 2017-09-20 2017-09-20

Publications (2)

Publication Number Publication Date
EP3460530A1 true EP3460530A1 (en) 2019-03-27
EP3460530B1 EP3460530B1 (en) 2022-03-30

Family

ID=63667745

Family Applications (1)

Application Number Title Priority Date Filing Date
EP18195642.6A Active EP3460530B1 (en) 2017-09-20 2018-09-20 One dimensional transmission detector for radiotherapy

Country Status (2)

Country Link
US (1) US10702717B2 (en)
EP (1) EP3460530B1 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20220168592A1 (en) * 2019-04-04 2022-06-02 Koninklijke Philips N.V. Fast generation of multi-leaf collimator (mlc) openings using hierarchical multi-resolution matching

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4803368A (en) 1987-12-28 1989-02-07 Siemens Medical Laboratories, Inc. Assembly and method for monitoring the lateral position of a beam of ionizing radiation
US4988866A (en) 1987-07-21 1991-01-29 Westerlund Kjell B Measuring device for checking radiation fields from treatment machines for radiotherapy
US6125335A (en) 1998-04-10 2000-09-26 Sun Nuclear Corporation Wide field calibration of a multi-sensor array
DE10143609A1 (en) 2001-09-06 2003-03-27 Schuster Medizinische Systeme Measurement of the radiation field used in radiation therapy and chemotherapy of patients using a planar radiation sensor support with sensors arranged along equally spaced radial lines to provide high-resolution measurements
US8119978B2 (en) 2006-07-10 2012-02-21 University Health Network Apparatus and methods for real-time verification of radiation therapy
US8160204B2 (en) 2006-11-03 2012-04-17 Ion Beam Applications S.A. Method and device for IMRT verification
US8183534B2 (en) 2007-11-21 2012-05-22 Frederic Lacroix Scintillating fiber dosimeter array
WO2016191883A1 (en) 2015-06-05 2016-12-08 University Health Network Sensors with virtual spatial sensitivity for monitoring a radiation generating device

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1917998A1 (en) * 2006-11-03 2008-05-07 Ion Beam Applications S.A. Method and device for online IMRT verification
JP2008173182A (en) * 2007-01-16 2008-07-31 Mitsubishi Heavy Ind Ltd Radiation irradiation method and radiotherapy apparatus controller
EP2116278A1 (en) * 2008-05-06 2009-11-11 Ion Beam Applications S.A. Device for 3D dose tracking in radiation therapy

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4988866A (en) 1987-07-21 1991-01-29 Westerlund Kjell B Measuring device for checking radiation fields from treatment machines for radiotherapy
US4803368A (en) 1987-12-28 1989-02-07 Siemens Medical Laboratories, Inc. Assembly and method for monitoring the lateral position of a beam of ionizing radiation
US6125335A (en) 1998-04-10 2000-09-26 Sun Nuclear Corporation Wide field calibration of a multi-sensor array
DE10143609A1 (en) 2001-09-06 2003-03-27 Schuster Medizinische Systeme Measurement of the radiation field used in radiation therapy and chemotherapy of patients using a planar radiation sensor support with sensors arranged along equally spaced radial lines to provide high-resolution measurements
US8119978B2 (en) 2006-07-10 2012-02-21 University Health Network Apparatus and methods for real-time verification of radiation therapy
US8160204B2 (en) 2006-11-03 2012-04-17 Ion Beam Applications S.A. Method and device for IMRT verification
US8183534B2 (en) 2007-11-21 2012-05-22 Frederic Lacroix Scintillating fiber dosimeter array
WO2016191883A1 (en) 2015-06-05 2016-12-08 University Health Network Sensors with virtual spatial sensitivity for monitoring a radiation generating device

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
B POPPE ET AL: "DAVID-a translucent multi-wire transmission ionization chamber for in vivo verification of IMRT and conformal irradiation techniques", PHYSICS IN MEDICINE AND BIOLOGY, vol. 51, no. 5, 15 February 2006 (2006-02-15), Bristol GB, pages 1237 - 1248, XP055543750, ISSN: 0031-9155, DOI: 10.1088/0031-9155/51/5/013 *
BONIN ET AL.: "A pixel chamber to monitor the beam performances in hadron therapy", NUCLEAR INSTRUMENTS AND METHODS IN PHYSICS RESEARCH, 2004, pages 674 - 686, XP004492326, DOI: doi:10.1016/j.nima.2003.11.104
ISLAM MOHAMMAD ET AL: "An integral quality monitoring system for real-time verification of intensity modulated radiation therapy", MEDICAL PHYSICS, AIP, MELVILLE, NY, US, vol. 36, no. 12, 5 November 2009 (2009-11-05), pages 5420 - 5428, XP012129815, ISSN: 0094-2405, DOI: 10.1118/1.3250859 *
NIKOS PAPANIKOLAOU: "Investigation of the convolution method for polyenergetic spectra", MED. PHYS., vol. 20, no. 5, 1993
WOLFGANG A.: "Beam Modelling for a Convolution/Superposition-Based Treatment Planning System", , MEDICAL DOSIMETRY, vol. 27, no. 1, 2002, pages 11 - 19
YANG, Y.; XING, L.: "Using the volumetric effect of a finite-sized detector for routine quality assurance of multi leaf collimator leaf positioning", MED. PHYS., vol. 30, 2003, pages 433 - 441, XP012012014, DOI: doi:10.1118/1.1543150

Also Published As

Publication number Publication date
US20190083816A1 (en) 2019-03-21
US10702717B2 (en) 2020-07-07
EP3460530B1 (en) 2022-03-30

Similar Documents

Publication Publication Date Title
Parodi et al. In vivo range verification in particle therapy
JP5580531B2 (en) Dosimetry equipment for measurement of radiotherapy equipment
Chandraraj et al. Comparison of four commercial devices for RapidArc and sliding window IMRT QA
US8874385B2 (en) Radiation detector and verification technique of positioning accuracy for radiation detector
US10471279B2 (en) Proton dose imaging method and apparatus
Amaldi et al. Construction, test and operation of a proton range radiography system
JP5791546B2 (en) Calibration method of radiation measuring apparatus and particle beam therapy apparatus
Boggula et al. Evaluation of a 2D detector array for patient-specific VMAT QA with different setups
US7432510B2 (en) Dosimeter based on a gas electron multiplier for dose measurements of therapeutic radiation
Markovic et al. Characterization of a two‐dimensional liquid‐filled ion chamber detector array used for verification of the treatments in radiotherapy
US8218725B2 (en) Method and device for online IMRT verification
Taylor et al. Proton tracking for medical imaging and dosimetry
Biasi et al. CyberKnife® fixed cone and Iris™ defined small radiation fields: Assessment with a high‐resolution solid‐state detector array
EP3061497B1 (en) Multi-resolution detectors for measuring and controlling a charged particle pencil beam
Guardiola et al. Preliminary microdosimetric measurements with ultra-thin 3D silicon detectors of a 62 MeV proton beam
Alhujaili et al. Quality assurance of Cyberknife robotic stereotactic radiosurgery using an angularly independent silicon detector
Rosenfeld et al. Medipix detectors in radiation therapy for advanced quality-assurance
EP3460530B1 (en) One dimensional transmission detector for radiotherapy
De Martin et al. On the evaluation of edgeless diode detectors for patient-specific QA in high-dose stereotactic radiosurgery
Wroe et al. Initial testing of a pixelated silicon detector prototype in proton therapy
Bucciolini et al. Dosimetric characterization of a bi‐directional micromultileaf collimator for stereotactic applications
Hollebeek et al. A new technology for fast two-dimensional detection of proton therapy beams
US20220390628A1 (en) Particle induced radiography system
Pardo et al. Verification of intensity modulated profiles using a pixel segmented liquid-filled linear array
US20220161061A1 (en) System for charged particle therapy verification

Legal Events

Date Code Title Description
PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE APPLICATION HAS BEEN PUBLISHED

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

AX Request for extension of the european patent

Extension state: BA ME

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE

17P Request for examination filed

Effective date: 20190926

RBV Designated contracting states (corrected)

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: EXAMINATION IS IN PROGRESS

17Q First examination report despatched

Effective date: 20201002

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: EXAMINATION IS IN PROGRESS

GRAP Despatch of communication of intention to grant a patent

Free format text: ORIGINAL CODE: EPIDOSNIGR1

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: GRANT OF PATENT IS INTENDED

INTG Intention to grant announced

Effective date: 20211022

GRAS Grant fee paid

Free format text: ORIGINAL CODE: EPIDOSNIGR3

GRAA (expected) grant

Free format text: ORIGINAL CODE: 0009210

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE PATENT HAS BEEN GRANTED

AK Designated contracting states

Kind code of ref document: B1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

REG Reference to a national code

Ref country code: GB

Ref legal event code: FG4D

REG Reference to a national code

Ref country code: CH

Ref legal event code: EP

REG Reference to a national code

Ref country code: DE

Ref legal event code: R096

Ref document number: 602018032875

Country of ref document: DE

REG Reference to a national code

Ref country code: AT

Ref legal event code: REF

Ref document number: 1479721

Country of ref document: AT

Kind code of ref document: T

Effective date: 20220415

REG Reference to a national code

Ref country code: IE

Ref legal event code: FG4D

REG Reference to a national code

Ref country code: LT

Ref legal event code: MG9D

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: SE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220330

Ref country code: RS

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220330

Ref country code: NO

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220630

Ref country code: LT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220330

Ref country code: HR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220330

Ref country code: BG

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220630

REG Reference to a national code

Ref country code: NL

Ref legal event code: MP

Effective date: 20220330

REG Reference to a national code

Ref country code: AT

Ref legal event code: MK05

Ref document number: 1479721

Country of ref document: AT

Kind code of ref document: T

Effective date: 20220330

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: LV

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220330

Ref country code: GR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220701

Ref country code: FI

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220330

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: NL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220330

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: SM

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220330

Ref country code: SK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220330

Ref country code: RO

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220330

Ref country code: PT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220801

Ref country code: ES

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220330

Ref country code: EE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220330

Ref country code: CZ

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220330

Ref country code: AT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220330

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: PL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220330

Ref country code: IS

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220730

Ref country code: AL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220330

REG Reference to a national code

Ref country code: DE

Ref legal event code: R097

Ref document number: 602018032875

Country of ref document: DE

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: DK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220330

PLBE No opposition filed within time limit

Free format text: ORIGINAL CODE: 0009261

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT

26N No opposition filed

Effective date: 20230103

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: MC

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220330

REG Reference to a national code

Ref country code: CH

Ref legal event code: PL

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: SI

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220330

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: LU

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20220920

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: LI

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20220930

Ref country code: IT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220330

Ref country code: IE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20220920

Ref country code: CH

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20220930

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: GB

Payment date: 20230927

Year of fee payment: 6

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: FR

Payment date: 20230925

Year of fee payment: 6

Ref country code: DE

Payment date: 20230927

Year of fee payment: 6

Ref country code: BE

Payment date: 20230927

Year of fee payment: 6

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: HU

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT; INVALID AB INITIO

Effective date: 20180920