Grid template position in interventional medicine

Abstract

A processing device 1 and method for determining the position and / or orientation of a multi-hole grid template in a medical interventional procedure are disclosed. An anatomical spatial information processing unit 2 receives and / or processes data representing a target spatial volume within a body, and a grid position sampler 4 generates candidate positions for the grid template. A quality calculator 5 calculates, for each candidate, a quality determination criterion indicative of the suitability of the candidate grid template position for the interventional procedure. A position selector 6 selects a position and / or orientation from the candidates based on the quality determination criterion. For each candidate, a spatial relationship between each grid hole trajectory and the target volume is determined, and the quality determination criterion takes into account at least this spatial relationship.

Description

[Technical field]

[0001] The present invention relates to the field of planning interventional procedures such as brachytherapy and ablation procedures. In particular, the present invention relates to devices and methods for planning procedures in which a grid or template is used to support and guide an interventional tool, such as a needle or catheter, that is inserted into an anatomical target volume through selected holes in the grid. [Background technology]

[0002] In localized therapeutic procedures, one or more focal sources are used near and / or within a treatment volume in the body, such as a cancerous lesion. Such focal sources are incorporated into or positioned by an interventional focal therapy delivery device, such as a needle or catheter. The focal source provides a highly localized energy irradiation to damage cells within the treatment volume, while protecting healthy tissue outside this volume. Generally, but not necessarily strictly limited to, the focal source deposits energy according to an inverse square law to achieve this advantageous highly focused (localized) therapeutic effect. There are many modalities of localized therapy in the art. Perhaps the most well-known modality is brachytherapy, in which one or more radioactive seeds (sealed radioisotopes) are placed in or near the area in need of treatment. Other examples include thermal ablation, microwave ablation, ultrasound ablation, radiofrequency ablation, photothermal therapy, laser therapy, etc. Certain sensitization techniques, such as photodynamic therapy, can also be used. While in the strictest sense cryoablation does not rely on a focal energy source, but instead uses a focal heat sink, it can be considered entirely similar and therefore can be considered yet another example of a focal treatment modality.

[0003] It is known in the art to use a rigid or flexible grid template placed on the patient's body to guide and support a focal therapy delivery device, such as a needle.

[0004] For example, in prostate brachytherapy, needles or catheters are inserted into the prostate and / or surrounding tissue through such an alignment template that is positioned relative to the perineum. In this way, the template guides, supports and (temporarily) fixes the needle or catheter to reach its intended delivery site in the body with good positional accuracy and ensures that the treatment tool remains at the intended site for the required predefined time during the procedure.

[0005] The present disclosure is directed to providing means and methods for planning and executing precise positioning of such grid templates relative to the body and predefined volumes therein that are particularly useful for focal treatments, although embodiments are not necessarily so limited. Other interventional procedures, such as biopsies, in which grid templates can be used to insert tools at precise locations within bodily tissue, would also benefit from the present invention.

[0006] Typically, medical imaging techniques are used to identify and segment regions of interest within the body, such as the primary treatment volume and nearby organs that may be at risk. This information can then be used to perform manual forward planning or automatic inverse planning to discover a set of delivery parameters, including, for example, the selection of the catheter or needle position on a grid template, its insertion depth, dwell (or ablation) time and / or treatment power parameters (e.g., radioisotope flux or ablation power). For example, in a brachytherapy inverse planning method, a grid hole trajectory that intersects the lesion area to be treated defines a target lesion intersection segment that is used to sample a set of dwell positions where the radiation source can dwell for a given optimal time. In a typical plan, multiple such available trajectories are selected and used in a sequential treatment pattern, e.g., a brachytherapy radiation source (or ablation tool) is inserted through each selected hole according to the determined parameters, e.g., according to the determined dwell time and insertion depth.

[0007] However, this planning is typically constrained by the position of a grid template that includes a number of through-holes (e.g., arranged in a grid pattern) that define available positions for the delivery device on this predefined grid, e.g., perpendicular to the grid template and along paths through the number of through-holes.

[0008] For example, especially for small lesions, suboptimal grid locations may limit the best quality of the treatment plan that can be achieved under the planning constraints. In a worst-case scenario, a very small lesion may be placed entirely between grid hole trajectories. Optimal treatment within these constraints would therefore increase the number of local treatment sources, thereby surrounding the lesion as much as possible, potentially resulting in increased undesired damage to the surrounding, possibly healthy tissue, e.g. increasing the risk of undesired side effects. Unfortunately, such scenarios are not uncommon in clinical practice. For example, clinical experts include a margin, e.g. 1-2 centimeters, around the total tumor volume in order to expand the target volume and obtain several intersecting grid trajectories. In such an event, there is no trajectory available for optimization without considering the additional margin, and thus no suitable brachytherapy source or ablator stay position. However, this additional margin also means that healthy tissue near the macroscopic tumor may be adversely affected, e.g. destroyed by the treatment, resulting in potential side effects.

[0009] However, this approach also has advantages. By restricting the possible needle positions to a predefined grid, the complexity of the planning algorithm is reduced, e.g., remaining tractable. For example, the needle positions are constrained to a possibly large, but finite (discrete) number, so direct search and / or integer programming planning strategies can be used. Similarly, the number of dwell positions (with substantially different effects) is also constrained to a discrete number, because the size of the brachytherapy seeds or ablation tools is finite. Although robotic systems that allow substantially continuous positioning of the interventional tool relative to a predefined coordinate system are known in the art, the template-based approach still offers an advantageously cost-effective, simple and fast solution, and is preferred in many scenarios.

[0010] Document WO2016 / 059603A1 discloses an interventional treatment system including at least one catheter configured to be inserted into a target object, and at least one controller configured to acquire a reference image dataset including a plurality of image slices forming a three-dimensional image of the target object, define a restricted area in the reference image dataset, determine location constraints for the at least one catheter according to at least one of the planned catheter intersection points, the peripheral boundary of the target object, and the restricted area defined in the reference dataset, determine at least one of the position and orientation of the distal end of the at least one catheter, and / or determine a planned trajectory of the at least one catheter according to the at least one position and orientation determined for the at least one catheter and the location constraint.

[0011] EP1374949 describes an exemplary approach to brachytherapy planning. This reference discloses a real-time radiation therapy planning system. A three-dimensional image segmentation algorithm is used to determine specific organs in the target anatomical region. A treatment plan for performing radiation therapy is then determined, which defines the number and location of hollow needles in the anatomical region and the radiation dose to be delivered. A single or multi-objective anatomical based genetic optimization algorithm is used to determine in real time the optimal number and location of hollow needles, the location of the energy emitting source in each hollow needle, and the dwell time of the energy emitting source at each location. The needles are then inserted into the anatomical region under the guidance of a template or a guide tool, and the energy emitting source is delivered through the hollow needles. Furthermore, for post-planning purposes, the achieved needle position and dwell time are determined based on the three-dimensional image information. A particular embodiment described in this disclosure relates to a motorized template without holes, which uses a single guide tube through which the needles can be inserted. This guide tube is positioned at each location of a virtual template grid, so that the needle position relative to the template is not limited by a physical grid of holes. Therefore, the virtual grid configuration is limited only to the diameter of the needle being used. Summary of the Invention [Problem to be solved by the invention]

[0012] It is an object of embodiments of the present invention to provide a simple, effective, efficient and / or good means and / or method for determining the position and / or orientation of a grid template for supporting and guiding an interventional tool, such as a needle, catheter, etc., that is inserted into a target anatomical volume, e.g., on or near the surface of a patient's body, through holes in the grid, in an interventional procedure, such as a brachytherapy or ablation procedure. [Means for solving the problem]

[0013] Devices and methods according to embodiments of the present invention achieve the above objectives.

[0014] An advantage of embodiments of the present invention is that the position of the grid is optimized to provide good accessibility (e.g., good physical access and / or coverage) of the target treatment volume within the body by a delivery path (for a needle, catheter, or other interventional tool) through the grid holes of the grid when placed in a predetermined position.

[0015] An advantage of embodiments of the present invention is that the delivery path through the grid holes is optimized to intersect the target treatment volume as much as possible.

[0016] An advantage of embodiments of the present invention is that it is possible to achieve good therapeutic effects of procedures such as ablation or brachytherapy, for example in the treatment of unresectable tumors. It is a further advantage that undesirable side effects can be reduced or avoided, for example damage to healthy tissue can be avoided or reduced.

[0017] A further advantage is that the treatment dose (e.g. radiation dose, ablation energy) can be reduced by enabling delivery of a precise dose within the treatment volume. A further advantage is that a good recovery of the patient after treatment (e.g. fast, e.g. without undesired side effects) can be achieved. For example, to reduce undesired damage to nearby healthy tissue, the total radiation dose delivered or the thermal ablation zone should follow the contour of the tumor as closely as possible (possibly including a margin around the tumor as deemed appropriate by the clinician and with as low an impact as possible due to constraints imposed by the grid).

[0018] An advantage of embodiments of the present invention is that the presently disclosed approach can be applied to a wide variety of interventional procedures where accurate positioning of an interventional tool in a volume of interest inside the body via a delivery path defined by a template grid is advantageous. For example, the approach is applied without requiring knowledge of the specific procedure being applied, but with the sole consideration of providing good access over a given volume of interest inside the body, e.g., via a delivery path provided by a template grid, preferably sampled at as many locations as possible throughout the volume, and preferably distributed as uniformly as possible over the entire volume. Advantageously, further constraints can be taken into account, such as at-risk areas to be avoided, such as tissue that may be damaged by physical insertion of an interventional tool into and / or through the at-risk area, and / or by an energy source, such as an ablation tool or brachytherapy seeds, if placed in or near the at-risk area.

[0019] An advantage of embodiments of the present invention is that the approach is readily available, inexpensive, disposable, and / or can be used with simple template grids, for example as known in the art.

[0020] An advantage of embodiments of the present invention is that the use of a simple template grid avoids compatibility issues with more complex positioning aids, such as compatibility with magnetic resonance imaging, ultrasound and / or X-ray imaging.

[0021] An advantage of embodiments of the present invention is that it is possible to avoid complex robotic systems for substantially free positioning of an interventional tool, such as robotic systems with substantially continuous positioning of the interventional tool with 3 (4, 5 or 6) (rotational) translational degrees of freedom, while providing many advantages of the complex systems in terms of precision. For example, robotic systems offer a significant degree of positioning freedom and allow good optimization of the treatment, but can be expensive, for example in terms of initial investment, maintenance requirements and training of clinical staff. An advantage of embodiments is that it is possible to achieve high quality treatment using inexpensive, simple and readily available treatment systems, such as grid templates and / or treatment planning systems, allowing for better optimization than previously achieved, for example by treatment relying on manually positioned grid templates.

[0022] An advantage of embodiments of the present invention is that a good positioning of the template grid is achieved using an optimization procedure that can be decoupled from the treatment plan based on the determined position, thus reducing the overall complexity of the plan and maintaining a globally tractable treatment plan. The use of a template grid has the advantage of reducing the complexity since it limits the possible positions of the interventional tool considered in the plan, while at the same time providing a set of good available positions to use in the plan. Due to the limited number of available positions, a direct search or similar optimization strategy can be used for the plan, which has the advantage of allowing for example the discovery of a global optimum of the treatment, avoiding suboptimal (locally optimal) solutions.

[0023] An advantage of embodiments of the present invention is that the reliance on the experience of the clinician is reduced. For example, the automated positioning approach disclosed herein increases the efficiency of clinical personnel, avoids problems due to limited experience of physicians in training, and increases the efficacy and / or reproducibility of treatments, for example, avoiding unnecessary variability due to manual positioning of a template grid. By reducing variability between treatments, positioning of the template grid aids in treatment optimization, since the effects of other treatment variability can be more easily detected when analyzing clinical outcomes.

[0024] Furthermore, the available intervention paths can increase the robustness of the treatment against small grid and patient movements by ensuring a good coverage of the target volume: for example, with a good initial position of the template grid, more suitable paths become available, allowing the dose (radiation or ablation) to be distributed over a larger number of treatment delivery points, thus reducing the risk of inaccurate delivery at only a few points.

[0025] An advantage of embodiments of the present invention is that especially small tumors can be effectively treated without increasing the (spatial) tolerance due to inaccurate positioning of the template grid. For example, suboptimal grid positioning can significantly reduce the quality of the delivered treatment plan, especially for very small lesions. In a worst case scenario, a very small lesion may be completely located between the grid hole trajectories, making it inaccessible by the inserted catheter or needle. In such a scenario, the margin around the target needs to be increased to reach the treatment volume, but at the expense of increased dose (or thermal ablation) to nearby healthy tissue.

[0026] It is an advantage of embodiments of the present invention that they provide a simple user interface and / or alignment procedure for positioning a template grid.

[0027] In a first aspect, the present invention relates to a processing device for determining the position and / or orientation of a grid template relative to a human or animal body during a medical intervention procedure, such grid template comprising a plurality of holes defining a corresponding plurality of grid hole trajectories, adapted to support and guide at least one interventional tool along such grid hole trajectory when inserted into the body through at least one of said holes during said intervention procedure.

[0028] The processing device detects at least one target in the body. Target The apparatus includes an anatomical spatial information processing unit for receiving and / or processing data representative of the volume.

[0029] The processing device detects at least one target in the body. Target for generating a plurality of candidate positions and / or orientations of the grid template relative to the volume Generating Unit Includes.

[0030] The processing device includes a quality calculator for calculating, for each candidate position and / or orientation of the plurality of candidate positions and / or orientations, at least one quality determination criterion representative of the suitability of the candidate position of the grid template for the interventional procedure.

[0031] The processing device includes a location selector for selecting a location and / or orientation from a plurality of candidate locations and / or orientations based on at least one quality determination criterion.

[0032] The quality calculator is adapted to determine, for each candidate position and / or orientation, a spatial relationship between each grid hole trajectory of the grid template and at least one said target volume when located according to the candidate position and / or orientation, For example, the quality calculator is adapted to determine an intersection of each grid hole trajectory of the grid template with at least one said target volume when located according to a candidate position and / or orientation.

[0033] The quality calculator is also adapted to calculate, when constrained by the determined spatial relationship, for each candidate position and / or orientation at least one quality determination criterion based on the determined spatial relationship and / or by taking into account at least a value indicative of a geometric overlap and / or proximity of the grid hole trajectory to at least one target volume base, which is indicative of a therapeutic effect measure of the medical intervention procedure.

[0034] In a device according to an embodiment of the invention, this value comprises a value indicative of a therapeutic efficacy measure of a medical intervention procedure when constrained by said determined spatial relationship, said therapeutic efficacy measure being a measure of the therapeutic efficacy of at least one of said targets when one or more radiation sources or ablators are located along said grid hole trajectory. Target In the room R It represents the radiation dose or ablation effect.

[0035] For example, if one or more radiation sources (e.g., emitting a predetermined radiation fluence or having a predetermined intensity in Becquerels) are located along the grid hole trajectory, the treatment efficacy measure may be the percentage of targets that receive at least a predetermined radiation dose (e.g., at least a predetermined value in Gray). Target Similarly, when one or more radiation sources (e.g., emitting a given radiation fluence or having a given intensity in becquerels) are located along the grid hole trajectory, the treatment efficacy measure may be a measure of the target volume or a measure of the relative (e.g., percentage) volume of the target volume. Target The radiation dose received by the volume may include a sum, average, or other summary statistic of the radiation dose received by the volume, or a given volume fraction thereof.

[0036] For example, when one or more ablation sources are positioned along the grid hole trajectory (e.g., according to a predetermined source power and / or dwell time), the treatment efficacy measure may be the percentage of the target receiving a predetermined amount of ablation energy. Target For example, treatment efficacy measures may include the percentage or absolute volume of the target ablated. TargetThis represents the absolute or relative volume of a volume.

[0037] For example, a treatment efficacy measure may be the number of targets that reach a predetermined temperature (e.g., exceed a predetermined temperature threshold) when ablated by one or more ablation probes while located along a grid hole trajectory, e.g., given a predetermined ablation power and / or dwell time and / or other predetermined ablation parameters. Target Similarly, treatment efficacy measures may include the average, maximum, minimum, or other summary statistics of the deposited ablation energy, or the temperature achieved by heating by the ablation source.

[0038] In a device according to an embodiment of the invention, this value comprises a first value indicating the degree of intersection of (e.g. all) of the grid hole trajectories of candidate positions and / or orientations with at least one target volume considered.

[0039] In a processing device according to an embodiment of the invention, the first value comprises a total number of grid hole trajectories that intersect with at least one target volume, and / or a total length of the line segments corresponding to said intersections, and / or an average or other statistical representative, such as for example a median, of the lengths of the line segments corresponding to said intersections.

[0040] In a processing device according to an embodiment of the present invention, the anatomical spatial information processing unit is adapted to receive and / or process data representing at least one spatial volume at risk within the body, such as a volume within the body, to avoid or reduce intersections with the grid hole trajectory and / or to avoid or reduce proximity of the grid hole trajectory to the volume at risk.

[0041] The quality calculator is adapted to determine, for each candidate position and / or orientation, an intersection of each grid hole trajectory of the grid template when located according to the candidate position and / or orientation with at least one said volume at risk. indicating a degree of intersection of at least the first value and (e.g., all) grid hole trajectories for the candidate locations and / or orientations with at least one volume at risk; It is adapted to calculate at least one quality decision criterion by taking the second value into account.

[0042] In a processing device according to an embodiment of the present invention, the second value comprises a total number of grid hole trajectories that intersect with at least one volume at risk, and / or a total length of line segments corresponding to said intersections with at least one volume at risk, and / or an average or other statistical representative value of the lengths of said line segments corresponding to said intersections with at least one volume at risk.

[0043] In a processing device according to an embodiment of the present invention, the quality calculator calculates at least a first said value, optionally The second value , and a third value indicating a minimum distance of the grid hole trajectory to the center of the at least one target volume or to the center of at least one of the at least one target volumes.

[0044] In a processing device according to an embodiment of the invention, the at least one quality decision criterion is a plurality of quality decision criteria, and said quality calculator is adapted to combine said plurality of quality decision criteria into a composite quality decision criterion according to a weighted sum, and a position selector is adapted to select from a plurality of candidate positions and / or orientations a position and / or orientation that reaches an extreme value of the composite quality decision criterion.

[0045] In a processing device according to an embodiment of the present invention, the at least one quality decision criterion is a plurality of quality decision criteria, and the position selector is adapted to select a first subset of the plurality of candidate positions and / or orientations based on a first quality decision criterion among the plurality of quality decision criteria, to select at least a second subset of the first subset based on a second quality decision criterion among the plurality of quality decision criteria that is different from the first quality decision criterion, and to further reduce the second subset, e.g., optionally via one or more nested subsets, using one or more further quality decision criteria.

[0046] A processing device according to an embodiment of the present invention includes a user interface for receiving priority settings from a user to select an ordered set or subset of a plurality of said quality determination criteria, wherein a position selector is adapted to select a first quality determination criterion and a second quality determination criterion according to the ordered set or subset.

[0047] In a processing device according to an embodiment of the present invention, Generating Unit is adapted to generate a plurality of candidate positions and / or orientations by translating and / or rotating a position and / or orientation representing the grid template over a plurality of different translation and / or rotation steps.

[0048] In a processing device according to an embodiment of the present invention, Generating Unit is adapted to translate and / or rotate over a number of different translation and / or rotation steps relative to an initial position and / or orientation.

[0049] A processing device according to an embodiment of the invention includes a user interface for receiving said initial position and / or orientation from a user.

[0050] In a processing device according to an embodiment of the present invention, Generating Unit teeth, - At least one target in the body Target Bo projecting the volume onto a plane tangent to the exterior of the body and / or to the surface of the body and determining an initial position as the center of said projection; - At least one target in the body Target The device is adapted to determine said initial position and / or orientation by determining said initial position by projecting a center of a volume onto a plane tangent to the exterior of the body and / or to a surface of the body.

[0051] In a processing device according to an embodiment of the present invention, the anatomical spatial information processing unit comprises an input port for receiving said data in the form of at least one segmented medical image, wherein at least one target in the body is Target At least one of the spatial volumes at risk within the volume and / or body is represented by a corresponding segmentation label.

[0052] In a processing device according to an embodiment of the invention, the anatomical spatial information processing unit comprises an input port for receiving said data in the form of at least one surface mesh and / or at least one target in the body. Target The information may include a parametric spatial descriptor of the volume and / or at least one spatial volume at risk within the body.

[0053] In the processing device according to the embodiment of the present invention, the anatomical spatial information processing unit is Target an input port for receiving said data in the form of at least one medical image representing at least one said volume at risk within a volume and / or body; Target and an image segmentation unit (3) for segmenting at least one said medical image, such as to determine at least one spatial volume at risk within the body, the volume being examined, and / or the body.

[0054] A processing device according to an embodiment of the invention includes a grid template alignment evaluator for receiving a position signal indicative of a physical position and / or orientation of a physical grid template and for providing a feedback signal indicative of the position and / or orientation selected by a position selector and / or indicative of the physical position and / or orientation and / or the relative position and / or orientation between said physical position and / or orientation and said selected position and / or orientation.

[0055] In a processing device according to an embodiment of the invention, the grid template alignment evaluator is adapted to concomitantly visualize the physical position and / or orientation and the selected position and / or orientation using a user interface, In a processing device according to an embodiment of the invention, the user interface is adapted to display said physical position and / or orientation and said selected position and / or orientation in different display styles.

[0056] In a processing device according to an embodiment of the present invention, the user interface is adapted to display a substantial alignment of the physical position and / or orientation with the selected position and / or orientation in a further different display style (from the style described above).

[0057] In a processing device according to an embodiment of the present invention, the feedback signal comprises an audio signal indicative of a measure of the discrepancy between the selected position and / or orientation and the physical position and / or orientation.

[0058] In a processing device according to an embodiment of the present invention, the feedback signal comprises an actuator signal for controlling one or more actuators adapted to position said physical grid template.

[0059] In a second aspect, the present invention relates to a clinical workstation comprising a device according to an embodiment of the first aspect of the invention.

[0060] In a third aspect, the present invention relates to a computer-implemented method for determining a position and / or orientation of a grid template relative to a human or animal body in a medical intervention procedure, said grid template comprising a plurality of holes defining a corresponding plurality of grid hole trajectories, said grid template adapted to support and guide at least one interventional tool along such grid hole trajectory when inserted into the body through at least one of said holes in said intervention procedure. The method comprises: determining a position and / or orientation of at least one target in the body; Target receiving and / or processing data representative of a volume of the body; and detecting at least one target within the body. Target Bo generating a plurality of candidate positions and / or orientations of the grid template relative to the volume, and for each candidate position and / or orientation, determining a grid hole trajectory of each grid template when located according to the candidate position and / or orientation, e.g., at least one target Target and determining a spatial relationship between at least one of the target volumes, such as an intersection of each grid hole trajectory with a volume of the target volume. The method further comprises calculating, for each candidate location and / or orientation of the plurality of candidate locations and / or orientations, at least one quality determination criterion representative of a suitability of the candidate location of the grid template for the interventional procedure.

[0061] The calculation of the at least one quality determination criterion considers a value indicative of a geometric overlap and / or proximity of the grid hole trajectories to the at least one target volume based on the determined spatial relationship. The calculation of the at least one quality determination criterion considers a value indicative of a therapeutic effect measure of the medical intervention procedure when constrained by the determined spatial relationship.

[0062] For example, at least one quality determination criterion is calculated by considering at least a first value indicative of the degree of intersection of the grid hole trajectory for a candidate position and / or orientation with at least one target volume, and selecting a position and / or orientation from a plurality of candidate positions and / or orientations based on the at least one quality determination criterion.

[0063] In a fourth aspect, the present invention relates to a computer program product comprising executable computer program code for performing a computer-implemented method according to an embodiment of the present invention, when the computer program product (the executable computer program code) is executed on a processor (e.g. a computing device).

[0064] The independent and dependent claims set out certain preferred features of the invention. Features from the dependent claims may be combined with features of the independent claims and with features of other dependent claims as appropriate and not necessarily as explicitly set out in the claims. [Brief description of the drawings]

[0065] [Figure 1] 1 illustrates a device and system according to an embodiment of the present invention; [Diagram 2] To illustrate an embodiment of the present invention, a diagram is shown illustrating intersecting segments of a grid hole trajectory, e.g., a therapy delivery path, with a target volume of interest, e.g., a cancerous lesion. [Diagram 3] FIG. 1 shows two different candidate locations for a grid template and the minimum distance between the grid hole trajectory and the geometric center of the volume of interest to explain an embodiment of the present invention. [Figure 4] FIG. 1 shows a visualization of physical grid locations and selected grid locations (i.e., target locations) to explain an embodiment of the present invention. [Diagram 5] FIG. 1 shows a grid template attached to a stepper device to explain an embodiment of the present invention. [Figure 6] FIG. 1 shows a table mount for supporting a grid template and a stepper device to explain an embodiment of the present invention. [Figure 7] FIG. 1 shows a diagram illustrating the controllable position and orientation of a grid template using a (eg, manually controllable) stepper device to illustrate an embodiment of the present invention. [Figure 8] FIG. 2 illustrates a method according to an embodiment of the present invention. DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0066] The drawings are schematic and non-limiting, elements in the drawings are not necessarily drawn to scale, and the present invention is not necessarily limited to the particular embodiment(s) of the invention shown in the drawings.

[0067] Notwithstanding the exemplary embodiments described below, the present invention is limited only by the scope of the appended claims, which are hereby expressly incorporated into this detailed description, and each claim, and each combination of claims permitted by the substructure defined by the claims, forms a separate embodiment of the present invention.

[0068] The term "comprising" as used in the claims is not limited to the features, elements or steps described below and does not exclude additional features, elements or steps. It therefore specifies the presence of the stated features without excluding the further or additional presence of one or more features.

[0069] In this detailed description, numerous specific details are set forth. Embodiments of the invention may be practiced without these specific details. Additionally, well-known features, elements, and / or steps have not necessarily been described in detail in order to make the disclosure clear and concise.

[0070] In a first aspect, the present invention relates to a processing device for determining the position of a grid template in a medical intervention procedure, where the grid template is used to support and guide at least one interventional tool, such as a needle, an ablation tool, a brachytherapy source or a catheter, into a target anatomical volume through at least one selected hole in the grid template.

[0071] Fig. 1 shows an exemplary device 1 according to an embodiment of the invention. The device includes, for example, a computer specially programmed to implement the device as described. Such a computer includes inputs and outputs, e.g., communication interfaces for receiving data and transmitting data, e.g., via a data carrier and / or a communication network interface. Such inputs and outputs also include user interface hardware, e.g., human interaction devices (e.g., keyboard, mouse, voice interpreter, touch interface, gyroscope, or accelerometer, etc.) for receiving input from a human user, a monitor for presenting information to the user, a speaker, a printer for rendering information on a physical carrier such as paper, a three-dimensional printer for generating a three-dimensional physical model of the data, and other such elements known in the art.

[0072] A computer includes a general purpose processor for executing instructions, e.g., computer code, and a memory for storing such instructions. A computer includes a memory for storing data, e.g., for operating according to instructions. The device is not necessarily limited to a general purpose computer, but includes, e.g., application specific integrated circuits (ASICs) and / or configurable processing hardware, e.g., field programmable gate arrays (FPGAs). Additionally, the device may be included in a single processing device, e.g., a computer, but may be distributed across multiple such devices operably connected to each other, such that, for example, the processes described herein are performed by the cooperation of a server and client devices, or a parallel processing system, e.g., a computing cluster.

[0073] The processing device 1 is adapted to determine the position of the grid template in a medical intervention procedure. Such a medical intervention procedure is a localized treatment, a biopsy, or another medical intervention performed by specifically inserting an intervention tool into a spatial volume. For example, a medical intervention procedure comprises a biopsy or a pathological mapping of an anatomical area, such as a lesion or (part of) an organ or tissue suspected to be caused by a disease or disorder. For example, a localized treatment refers to a percutaneous localized treatment, such as a percutaneous localized cancer treatment, such as a brachytherapy cancer treatment, an ablative procedure, such as cryoablation, thermal ablation, microwave ablation, (focused) ultrasound ablation, radiofrequency ablation, laser ablation, or a treatment in which a focused energy source is used to activate sensitized tissues and / or cells, such as photodynamic therapy. For example, a brachytherapy or an ablation modality is selected by a physician to treat a tumor, depending for example on the available means, the size and location of the lesion, the affected organ, etc. However, in such procedures, it is common to manually place a grid template on (or at least near) the surface of the skin and insert an interventional tool (e.g., biopsy needle, ablator, needle, or catheter preloaded or loadable with brachytherapy seeds, ...) into the body through one or more guide holes in the grid template, typically selected from a number of possible (selectable) holes. Another example of such a procedure is the implantation of a microdevice, e.g., a microstimulator, or other structure, e.g., a time-release capsule for releasing a medicinal compound to a specific location in the body.

[0074] The grid template is referred to as a leading grid frame or treatment grid. The grid template may be a rigid or flexible grid template as known in the art. The grid template is radiopaque, for example, to avoid interference with ionizing radiation used to image the patient when the grid template is in position for treatment. In this disclosure, terms such as "grid", "template", "template grid" and "grid template" are considered equivalent and interchangeable, i.e., refer to the same such leading grid frame for use in an interventional procedure. The grid template includes multiple holes, i.e., through-holes, at different locations on the grid template. For example, the grid template has a major surface with such holes distributed such that each hole protrudes through the grid template from the major surface to a surface opposite the major surface. For example, the grid template has a one-dimensional component that is substantially smaller than the other two (Cartesian) dimensional components, and typically the two "major" surfaces on opposite sides of the template span the larger two dimensions. One such major surface is adapted to contact (or be directed towards) the patient's skin in use, while the opposite major surface allows access to these holes, so that tools, such as needles, pins, catheters and other such elongated elements, can be inserted from this side to the contact side and along this path into the patient's body in a specific position (position of the holes) and direction (obtained, for example, by guiding the tool through the holes). Without being limited thereto, a grid template includes such holes arranged in rows and columns (hence a "grid", e.g. a rectangular grid). This does not exclude other distribution patterns of holes, such as, for example, polar grids, hexagonal grids. The distribution pattern can be regular or irregular, uniformly or non-uniformly distributed. Such a template can be disposable or reusable. The template can be rigid, flexible or have more or less limited flexibility. The grid holes are, for example, spaced apart from each other (i.e. from adjacent holes) at a distance of 1 mm to 20 mm.A typical example is a template with a distance between holes in the range of 2 mm to 10 mm, e.g., 3 mm to 8 mm, e.g., 5 mm. The holes have a diameter that fits a given needle gauge, e.g., in the range of 14 gauge to 18 gauge (but not limited to), e.g., 14g, 17g, or 18g needle gauge. The number of holes is (without limitation) in the range of, e.g., 9 to 1000, e.g., typically in the range of 25 to 225, e.g., 25, 36, 49, 64, 81, 100, 121, 144, 169, or 225 holes. These examples are based on a regular grid with an equal number of rows and columns, but it is to be understood that in embodiments the grid need not be a square grid or have a number of rows equal to the number of columns.

[0075] The grid template includes a number of holes (through holes) that define a number of corresponding grid hole trajectories for a needle, catheter, or other elongated interventional tool (when inserted through the holes). The grid template has a substantial width such that the direction of the grid hole trajectories is constrained to a particular direction, e.g., perpendicular to the face of the grid template (e.g., in the case of a flexible grid template, local surface curvature may or may not be taken into account). However, embodiments are not necessarily limited thereto. For example, other means may be used to control the direction of insertion, or the holes may be provided at an angle other than perpendicular to the surface.

[0076] Those skilled in the art will appreciate that the location of the grid template is important for the success of the procedure. Typically, unless otherwise instructed by a medical professional, only those grid hole trajectories that intersect (e.g., enter) the target volume, e.g., the lesion to be treated, are used for the interventional procedure, e.g., for treatment delivery. For example, in an inverse planning method for brachytherapy, intersecting segments of the target lesion, e.g., as shown in FIG. 2, are used to sample a set of dwell positions where it is determined that the radioactive source may dwell for a given optimal time. Thus, both automatic inverse planning solutions and manual forward planning solutions benefit from optimal leading grid positions. For example, it is preferable that the trajectory of the needle through the grid holes intersects the segmented anatomical structures as much as possible. The success and efficacy of the treatment, and the prevention of post-treatment complications, are particularly sensitive to the initial grid position for very small lesions to be treated.

[0077] When used, the device 1 determines the location of the grid template to be used in such a procedure, e.g., the optimal (or at least good) placement of the grid template to allow the spatial volume to be reached by inserting an interventional tool through one or more holes in the grid template. For example, this location depends on the size and location of the volume, e.g., the location and size of the tumor, and on the parameters of the grid template, e.g., the specifications provided by the grid template manufacturer. Automatically optimizing the grid template location advantageously reduces the dependency on the physician's experience and can increase the physician's efficiency by reducing or avoiding the time spent manually determining the appropriate location. For example, in conventional approaches, the physician needs to place the grid template, select a set of tool positions, and determine a corresponding set of tool parameters. Once the grid template is placed, solutions are known in the prior art for automatically (or in a computationally assisted manner) selecting tool positions on the grid hole trajectory provided by the grid template and determining the appropriate parameters. For example, a manual forward planning solution or an automatic inverse planning solution may be applied to find an optimal set of delivery parameters (eg, brachytherapy catheter positions, thermal ablation probe positions, dwell times, ablation times / power values, etc.).

[0078] For example, the processing device 1 may be configured to detect at least one target in the patient's body, through which at least one interventional tool is to be inserted according to a medical intervention procedure, via at least one hole in the grid template. Target The device is configured to receive spatial data representative of a volume.

[0079] For example, the processing device is also configured to receive spatial data representative of at least one spatial volume that is at risk within the patient's body that the interventional tool should avoid, e.g., by being away by a distance margin and / or by not drilling a hole along an insertion path. TargetFor convenience, the volume, and optionally the spatial volume that is exposed to at least one hazard, are hereinafter collectively referred to as the "volume of interest."

[0080] Thus, the processing device 1 detects at least one target in the patient's body. Target The medical imaging device further comprises an anatomical spatial information processing unit 2 for receiving and / or processing data representative of at least one spatial volume at risk within the body, and optionally an anatomical spatial information processing unit 2 for receiving and / or processing data representative of at least one spatial volume at risk within the body. The spatial information processing unit is or comprises an input port for receiving said data, for example via a digital communication network, a physical data carrier or a dedicated interface, using a medical imaging device, for example a bus interface. The data can be, for example, TargetThe data is received in the form of a medical image, such as a volumetric three-dimensional image of at least a portion of the body, including the volume and / or volume of space at risk. Such medical images include (e.g., three-dimensional or tomographic) images obtained by X-ray computed tomography (CT), magnetic resonance imaging, nuclear medicine imaging, ultrasonography, or other medical imaging modalities (e.g., but not limited to these examples, elastography, optical CT, photoacoustic imaging, magnetic particle imaging). The medical image is processed to define the volume of interest, for example, in an image overlay that annotates the voxels assigned to the volume of interest, or the medical image is composed of the definition of the volume of interest, i.e., a segmented image. Alternatively, the data is received in the form of a descriptor of the volume of interest, such as, for example, a surface mesh that defines the volume of interest, which is a set of vertices, edges, and faces that define the or each volume of interest as a polyhedral object in a reference coordinate space. Such faces include triangles, quadrilaterals, and / or other polygons (e.g., n-gons that are convex or even concave polygons). The descriptor may be defined as a parametric definition of the (each) volume of interest, e.g. as one or more spheres, cubes, Non-Uniform Rational Basis Splines (NURBS), or other parametric shapes, or in another suitable format, e.g. as a volumetric mesh that includes explicit volume information. A descriptor of a volume of interest typically defines the volume (or its enclosing surfaces) directly, although a probabilistic or fuzzy definition of the volume (e.g. being a probability map) is not necessarily excluded.

[0081] The anatomical spatial information processing unit 2 also includes a (volumetric) image segmentation unit 3, which segments the medical image to depict a volume of interest within the medical image, where, for example, the data does not yet include segmentation information or another form of definition of the volume referred to above. Suitable segmentation algorithms are known in the art and include manual segmentation (e.g., using a suitable user interface 7), semi-automatic segmentation, and / or automatic segmentation. Some examples of such segmentation approaches known in the art include, but are not limited to, voxel value thresholding, clustering methods, histogram-based methods, edge detection, region growing methods, partial differential equation-based methods, calculus of variations, graph partitioning methods, watershed transform-based methods, multi-scale methods, machine learning-based methods (e.g., trainable methods), and / or any combination thereof.

[0082] Thus, the processing device receives a definition of a volume of interest or determines such a definition by processing received data, e.g., image data, that already contains such information by implication. Thus, all regions of interest (such as target lesions or organs at risk) are segmented in a pre-processing step and provided to the device, e.g., by processing performed by the device (e.g., based on the received image data) or by a combination of both (e.g., an initial segmentation is provided and refined by an automatic, interactive or manual procedure performed by the device).

[0083] The processing device may, for example, generate a plurality of candidate positions of the grid template relative to the volume of interest in a shared (implicit or explicit) reference coordinate system. Generating Unit 4. Generating the plurality of candidate positions includes translating and / or rotating the grid template through a plurality of different translation and / or rotation steps. For example, Generating Unitsamples the positions (and optionally, orientations) of the grid template over a series of translation steps and / or rotation angles. The initial positions are automatically determined and / or specified by a user, for example, using the user interface 7 of the device. For example, the user provides a best guess for the initial grid position (x0,y0) based on the patient's anatomy and / or prior knowledge and / or experience to initialize the positioning method performed by the device. The user-defined initial positions are verified by the device to see if they can fully encompass all target volumes, such as defined cancer lesions. Although the term "sampling" is used, this sampling is a deterministic determination of candidate positions, for example by direct enumeration, or probabilistic sampling, for example by random sampling from a distribution function of candidate positions. Reference to "position" relates to position coordinates, angular orientation coordinates, and / or combinations thereof. Thus, a "position" should not be construed narrowly as only a point in space, but may also relate to one or more angles or orientations, or a combination of a point in space and an orientation, such as, for example, a localized (e.g., "localized" as opposed to a free or released vector) vector.

[0084] The initial position is determined automatically by using a fixed predefined position or by simple heuristics such as a projection of the target volume of interest onto the surface of the skin of the body, e.g. along a fixed coordinate axis or along a surface normal of the skin surface. For example, the initial grid position is automatically determined by the geometric centroid of the target volume, e.g. a single cancer lesion, projected by a two-dimensional axial projection onto a grid x,y plane (defined, e.g., by a plane parallel to the surface of the skin, or an approximation thereof). Similarly, in case of multiple target volumes, e.g. multiple lesions, the spatial average position of such projected centroids calculated for each target volume is used. Alternatives can easily be envisaged, such as using other measures of centrality or using the centroid of the convex hull of multiple target volumes. Such an automatic determination of the initial position is also performed if the initial position provided by the user is deemed inappropriate, e.g. allowing insufficient access to the target volume, see validation above.

[0085] The transformation steps are determined by a straightforward and deterministic enumeration using integer (e.g., signed) multiples of a predefined step size, for example at one, two, or three coordinates. Such step sizes (or step sizes of different translational and / or rotational components) may be pre-determined or may be configurable by the user using the user interface 7.

[0086] For example, a horizontal step size sx and a vertical step size sy (which may be equal) are used (predetermined or specified by the user using the user interface 7), and candidate positions for the initial position x0,y0 are defined as [x,y]=[x0+i.sx,y0+j.sy] for all combinations of i={-K,-K+1···,0,···,K-1,K} and j={-L,-L+1,-L+2,···,0,···,L-2,L-1,L} to cover a given area range 2K.sx×2L.sy. Consideration of additional steps in the z direction perpendicular to the x,y plane (or, more precisely, steps with a z component that is not constant over all the steps considered) is not necessarily excluded (nor is it necessarily included). Similarly, rotations using, for example, an angle step size sθ for rotations in the x,y plane are also considered. The rotation can also be considered outside the x,y plane, for example using two or even three complementary angular components. Thus, a discrete (finite) set of grid template positions (and / or orientations) spanning 1, 2, 3, 4, 5, or even 6 (rotational) translational components of the transformation is generated. To obtain a good position estimate, at least two components, for example two complementary components of the translation, such as an x ​​component and a y component, are preferred (although in some embodiments optimization in a single direction is useful). A third component, such as a rotation in the x,y plane, is also considered to find positions that allow further optimization. Similarly, including additional degrees of freedom results in a better optimized position estimate. However, the embodiments are designed with a trade-off between optimization quality and complexity in mind. For example, the complexity associated with the number of sampled positions can be scaled as a power O(nk) with the number of degrees of freedom k considered in the transformation step.

[0087] Although this processing step of generating multiple candidate positions is described as a separate step, those skilled in the art will appreciate that the iteration over multiple candidate positions can also be performed in-line, e.g., the positions / orientations do not need to be pre-computed or stored for use in further processing steps, but are determined "on the fly" in those further processing steps.

[0088] While an example of uniform grid sampling of coordinates (possibly including rotational components) has been given, it will be understood that embodiments are not necessarily limited to uniform sampling. For example, polar sampling is used to consider more locations closer to the initial position than those farther away. Alternatively, stochastic sampling of the position space, or a variable step size, for example, is used to increase the sampling density near the initial position. This position ( / orientation) sampling process may include, for example, sampling over horizontal ("x") and vertical ("y") positions that define a plane perpendicular to the grid template at the initial position, e.g., selected or determined to be tangent to the surface of the patient's skin at (or near) the initial position (not necessarily limited to a tangential orientation).

[0089] The displacement is also not necessarily limited to linear displacements (and / or corresponding rotations), but also takes into account the curvature of the body's surface. For example, the sampled positions (e.g. uniform Cartesian sampled positions) are mapped to the body's surface to constrain the sampled positions to those positions where the grid template is placed in contact with the skin. For example, if a flexible grid template is used, positions are sampled over a given area of ​​the skin, possibly including a (curved) in-plane rotational component. However, even in the case of a rigid template, the orientation of the grid template follows the contours of the body, thereby sampling positions on the skin's surface, but with the caveat that the grid template, due to its rigidity and the curvature of the skin, does not contact (nor can it intersect) the skin at all points of the grid template. Generating Unit is adapted to take this into account, for example by sampling positions (and / or relative orientations) on the surface of the skin (determined in a segmentation step based on the imaging data) and adjusting each position, for example in a direction perpendicular to the surface of the skin, to the minimum distance necessary to avoid a collision (i.e. an intersection of the model of the rigid template with the body). Alternatively or additionally, an orientation adjustment is determined at this (each) position to avoid a collision. In case of a flexible grid template, either the grid template is assumed to be flexible enough to fit the surface of the skin, or the degree of flexibility is taken into account to determine a typical displacement (and / or orientation adjustment) that avoids collisions, in order to exclude positions where the curvature of the skin is too strong, or to maximally bend the grid template to fit as well as possible to the local curvature within the physical tolerance for bending.

[0090] The device 1 includes a quality calculator 5 for calculating, for each candidate location of a plurality of candidate locations, at least one quality decision criterion, which represents the suitability of the candidate location of the grid template for the intervention procedure. Reference to the "quality" calculator is made only for ease of disclosure, e.g. to avoid confusion with other functions, and does not mean anything more than what is described as the performed function of this component. Even though a quality decision criterion is referred to that generally increases in value, e.g. when the suitability of the location is high, a cost function that generally decreases, e.g. when the suitability increases, is also fully equivalent and therefore considered as a "quality decision criterion". The quality decision criterion is also referred to as an objective or objective function.

[0091] The quality calculator 5 is adapted to determine, for each candidate position and / or orientation, a spatial relationship between each grid hole trajectory of the grid template when located according to the candidate position and / or orientation and at least one (e.g. each) target volume.

[0092] For example, the quality calculator 5 is adapted to determine the intersection of each grid hole trajectory with each target volume or each target volume, as for example illustrated in Fig. 2. The quality calculator also determines the spatial relationship of each grid hole trajectory with each endangered volume or each endangered volume, such as for example the intersection of each grid hole trajectory with each endangered volume or each endangered volume.

[0093] For example, the quality calculator includes a model of the grid template parameterized to be evaluated as a function of the grid template's position (e.g. each candidate position) and / or orientation as parameters in the coordinate space in which the volume of interest is specified (or applying trivial coordinate transformations to allow the determination of said intersections). Although reference is made to a "model", such a model shall be understood to be particularly simple, for example a set of localized vectors defining each grid hole trajectory with respect to a local coordinate frame, where a transformation of this local coordinate frame depending on the position / orientation under evaluation allows the determination of the grid hole trajectory for each position / orientation and thus the determination of said intersections. It shall also be understood that such a model may be quite complex, for example taking into account deformation properties of the grid template. The quality calculator shall be understood to include multiple models of different grid templates, or a model of a grid template specified in a form that depends on grid template parameters such as, for example, grid hole spacing, grid configuration, grid dimensions, etc. User interface 7 is used to select a model corresponding to a grid template intended for use in the procedure and / or configure such grid template parameters (although embodiments in which only a single fixed grid template model is implemented are not necessarily excluded).

[0094] Each grid hole trajectory corresponds to a line segment passing through a hole in the grid template in a predetermined direction relative to the grid template's orientation, such as, for example, the direction in which the hole penetrates the template perpendicular to the grid template's main surface. The predetermined direction thus corresponds to the direction in which an interventional tool, e.g., a needle or catheter (or other types of generally elongated interventional tools for insertion into the body, including but not limited to needles or catheters), would be guided by the grid template when inserted through the hole. An intersection thus corresponds to such a line segment that coincides with the volume of interest. A calculated intersection may include an explicit identification of such a line segment, but may additionally or alternatively include a value derived therefrom, such as the length of the line segment, or simply a Boolean indicator indicating whether the grid hole trajectory intersects the volume of interest, e.g., whether an intersection exists, e.g., whether the length of the line segment is different from zero.

[0095] Calculating at least one quality determination criterion includes, when constrained by the determined spatial relationship, calculating at least one quality determination criterion for each candidate position and / or orientation based at least on the determined spatial relationship and / or by considering a value indicative of a geometric overlap and / or proximity of the grid hole trajectory to at least one target volume base, which is indicative of a therapeutic effect measure of the medical intervention procedure.

[0096] Calculating the at least one quality determination criterion comprises calculating a value indicative of a treatment effect measure of the medical intervention procedure when constrained by the determined spatial relationship. The treatment effect measure thus comprises a treatment effect measure that is indicative of a treatment effect measure of at least one (or each) target when one or more radiation sources or ablators are located along said grid hole trajectories. Target Bo In the room R It represents the radiation dose or ablation effect.

[0097] For example, if one or more radiation sources (e.g., emitting a predetermined radiation fluence or having a predetermined intensity in Becquerels) are located along the grid hole trajectory, the treatment efficacy measure may be the number of targets that receive at least a predetermined radiation dose (e.g., at least a predetermined value in Gray). Target Similarly, when one or more radiation sources (e.g., emitting a given radiation fluence or having a given intensity in becquerels) are located along the grid hole trajectory, the treatment efficacy measure may be a measure of the target volume or a measure of the relative (e.g., percentage) volume of the target volume. Target The present invention may provide a method for determining the amount of radiation received by a volume fraction, including the sum, mean, maximum, minimum, median, a predetermined percentile (e.g., first and / or third quartile), or other summary statistic, of the radiation dose received by the volume, or a predetermined volume fraction.

[0098] For example, when one or more ablation sources are positioned along the grid hole trajectory (e.g., according to a predetermined source power and / or dwell time), the treatment efficacy measure may be a measure of the number of targets that receive a predetermined amount of ablation energy. Target For example, treatment efficacy measures may include the percentage or absolute volume of the target ablated. Target Bo It represents the absolute or relative volume of a volume.

[0099] For example, a treatment efficacy measure may be the number of targets that reach a predetermined temperature (e.g., exceed a predetermined temperature threshold) when ablated by one or more ablation probes while located along a grid hole trajectory, e.g., given a predetermined ablation power and / or dwell time and / or other predetermined ablation parameters. Target Similarly, treatment efficacy measures may include the average, maximum, minimum, or other summary statistics of the deposited ablation energy, or the temperature achieved by heating by the ablation source.

[0100] For example, for each candidate position / orientation, an optimal position of a radiation source or ablator is determined along each trajectory to determine the treatment effect measure, or multiple such positions along each trajectory are considered, e.g., two or more, e.g., a predetermined number of sources are distributed along the trajectory. The device includes a treatment planning component that determines an appropriate treatment plan constrained by the candidate positions / orientations of the grid to determine the treatment effect measure. As known in the art, forward planning algorithms can be used to determine the optimal position of a radiation source or ablator along each trajectory. R A treatment effect measure is determined by applying a simplified, e.g. less accurate, treatment planning algorithm to each candidate position / orientation. However, if a simplified, e.g. less accurate, treatment planning algorithm is used to determine an approximation of the treatment effect measure and to use a specific position / orientation of the grid template, an efficient estimation of the treatment effect is possible, and after selecting the best position / orientation based on this less accurate measure, a more detailed, e.g. more accurate, treatment planning algorithm is applied to propose the treatment plan to be performed. Furthermore, in embodiments, a simplified approach is applied that, for reasons of efficiency, does not require a detailed treatment plan to evaluate the suitability of the candidate positions / orientations, e.g. based on a measure of trajectory intersection with the target volume. Such an approach advantageously allows to determine an appropriate configuration of the grid template without requiring detailed knowledge of the applied procedure and is applied regardless of the specific procedure to be performed, e.g. by relying only on geometrical considerations. Nevertheless, as already mentioned, a more complex approach that advantageously uses knowledge of the applied procedure is not necessarily excluded in embodiments according to the invention.

[0101] For example, calculating the at least one quality determination criterion includes calculating (at least one) target value indicative of the degree of intersection of all grid hole trajectories with the target volume (or all target volumes) when the grid template is located (and / or oriented) according to the candidate position. The target value includes the total number of grid hole trajectories intersecting the target volume. The target value includes the total length of line segments intersecting the target volume, e.g., the sum of the lengths of all such intersecting line segments. The target value includes the average length of line segments intersecting the target volume, e.g., the sum of the lengths divided by the number of grid hole trajectories intersecting the target volume. The at least one quality determination criterion is calculated such that a higher corresponding value of the target value is more indicative of the suitability of the candidate position of the grid template for the intervention procedure. Thus, the position / orientation of the grid template is preferentially selected based on the quality determination criterion to maximize the total number of grid holes available to reach the target volume and / or maximize the total path length that an interventional tool can be positioned in the target volume via the grid hole trajectories.

[0102] Calculating the at least one quality decision criterion includes calculating a (at least one) risk value indicative of the degree of intersection of all grid hole trajectories with the at-risk volume (or all volumes) when the grid template is located (and / or oriented) according to the candidate position. The risk value includes the total number of grid hole trajectories intersecting the at-risk volume. The risk value includes the total length of line segments intersecting the at-risk volume, e.g., the sum of the lengths of all such intersecting line segments. The risk value can include the average length of line segments intersecting the at-risk volume, e.g., the sum of the lengths divided by the number of grid hole trajectories intersecting the at-risk volume. The at least one quality decision criterion is calculated such that the higher the corresponding value of the risk value, the lower the suitability of the candidate position of the grid template for the intervention procedure. Thus, the position / orientation of the grid template is preferentially selected based on the quality decision criterion that minimizes the total number of grid holes available to reach the at-risk volume and / or minimizes the total path length that the intervention tool is placed in the at-risk volume via the grid hole trajectories.

[0103] Calculating the at least one quality decision criterion comprises calculating (at least one) centrality value indicative of a minimum distance of the grid hole trajectories to the center of the target volume. The distance is a Euclidean distance, but other suitable distance decision criteria of any integer or non-integer (still greater than 1) order p are not necessarily excluded, such as the Manhattan (or "taxi") distance r, the Chebyshev (or "chessboard" or "infinity norm") distance, or in general the Minkowski distance.

[0104] The minimum value is calculated as the minimum over all such distances calculated for all grid hole trajectories of the candidate positions / orientations under consideration. These distances are the distances between the grid hole trajectories (or equivalently, the determined intersections with the target volumes) and the centers of the target volumes, which may be the geometric center, center of gravity, centroid, center point, center of mass, Chebyshev center, center of the convex hull of the target volumes, center of the smallest sphere enclosing the target volumes, or another geometric measure of the center of a spatial volume. If multiple target volumes are defined, the center may refer to the (combined) geometric center (or other measure of centrality) of the collective target volumes, or this minimum distance is calculated for each target volume separately. For example, the centrality value may indicate the mean, median, average, or similar measure of the minimum distances obtained for the target volumes. Alternatively, the centrality value may indicate the maximum of the minimum distances between the grid hole trajectories and the centrality measures of the target volumes over all target volumes. Thus, the grid template positions / orientations are preferentially selected based on a quality decision criterion that minimizes the shortest distance to the center of the target volume (or each target volume) that can be reached via the grid hole trajectories. Therefore, the selection of grid location / orientation that maximizes the probability of having a grid hole trajectory that passes through the center of the target region is recommended. This improves the symmetry of the optimized interventional procedure plan with respect to the shape of the lesion at hand and improves robustness against minute patient / couch movements during the execution of the interventional procedure. For example, Figure 3 shows two different candidate locations for the grid template. It can be seen that the candidate location of the grid template 23, as shown on the left, achieves a smaller average nearest grid distance than the one shown on the right. In Figure 3, the center of gravity 21 is shown for two target volumes 22 (projected onto the plane of the grid template).

[0105] References to "target values," "risk values," and "centrality values" should be viewed as mere designations for the sake of clarity and to avoid confusion, and do not imply any special properties beyond those explained above. Similarly, the values ​​may be referred to as "first values," "second values," and "third values," but no particular ordering attributes are implied by such numerical designations.

[0106] Calculating at least one quality decision criterion includes calculating multiple such quality decision criteria, such as one or more of the mentioned target values ​​(e.g. based on the number of intersection segments, the total length of the intersection segments, and / or the average of the intersection segments), one or more of the mentioned risk values ​​(e.g. based on the number of intersection segments, the total length of the intersection segments, and / or the average of the intersection segments), and / or one or more of the mentioned centrality values ​​(different measures of centrality and / or approaches to consider multiple target volumes). These quality decision criteria are combined into a composite quality decision criterion, for example by calculating a weighted sum of the different quality decision criteria. The weights of such a weighted sum are pre-determined or configurable, for example received as input from a user via the user interface 7. Each quality decision criterion is appropriately scaled and, for example, normalized so that for comparable weight values, the weights have a comparable effect. For example, each value forming a term of the composite quantity decision criterion is divided by its known or assumed maximum value, allowing to select a weight ranging from 0 to 1, indicating "irrelevant" or "not taken into account" or "maximum relevance", respectively. It should also be noted that some of the above quality decision criteria should be considered as penalties, while others as targets, and therefore need to be combined with the appropriate sign, inversion of the other scalings to take this into account. For example, the number of segments intersecting with the at-risk volume is subtracted from the number of segments intersecting with the target volume, or otherwise has its effect inverted before being combined into the composite quality decision criterion (in addition to possible normalization of these terms, as discussed above). Alternatively or additionally, multiple quality decision criteria are evaluated separately, e.g., first selecting a suitable position / orientation based on the extreme value reached for the first quality decision criterion, then using a second (third,...) quality decision criterion to further reduce the number of selected positions / orientations, and finally arriving at one (or a small number) suitable positions / orientations.In particular, some of the exemplary quality decision criteria discussed above render integer values, e.g., number of intersections, allowing multiple candidate positions / orientations to render the same value for such individual quality decision criterion, e.g., maximum (or minimum) of value. Similarly, in case of non-discrete (continuous) quality decision criteria, a (predetermined or configurable) tolerance threshold is considered to select multiple positions / orientations that are deemed sufficiently close to an extreme value of a value to be considered in a further step using another quality decision criterion, in order to further reduce the number of selected positions / orientations.

[0107] The device 1 includes a position selector 6 for selecting a position and / or orientation from a plurality of candidate positions / orientations based on at least one quality decision criterion, e.g. to select a position / orientation that reaches an extreme value (maximum or minimum) for one of the at least one quality decision criterion and / or for a composite quality decision criterion.

[0108] For example, in an embodiment where multiple quality decision criteria are evaluated separately (i.e., when these quality decision criteria are not evaluated together as a composite quality decision criterion), the user interface 7 is adapted to receive a priority setting from the user to select the quality decision criteria to be evaluated and their evaluation order. Thus, to make a first selection of candidate positions / orientations, at least a first quality decision criterion is selected by the user, and to make a second selection, at least a second quality decision criterion is selected by the user from the candidate positions / orientations selected by the first selection. To further reduce the number of selected positions / orientations, a third or more quality decision criterion is selected by the user. Thus, the selector 6 selects positions / orientations by implementing a lexicographical order method, in which a prioritized set of targets (quality decision criteria) is provided by the user. However, in other embodiments, such priorities may be fixed and predefined, or only a single (e.g. composite) quality decision criterion is used.

[0109] In a first illustrative example, the selector implements a lexicographical order method, where a prioritized set of objectives is determined by the user, such as, for example, a highest priority for achieving the maximum number of grid hole trajectories intersecting the target volume, and a lower priority for achieving the maximum average intersection segment of the grid hole trajectories with the target volume. In such an approach, the selector first selects the candidate position / orientation that achieves the maximum number of grid hole trajectories intersecting the target volume. If this maximum is achieved by a single grid hole trajectory, this solution is provided as output without even processing the decision criteria / objectives of the second priority level. If not, as a second step, the best candidate position / orientation is selected that also produces the maximum average intersection segment length within the subset of grid positions / orientations that achieve the maximum number of said intersections. As mentioned above, other decision criteria, possibly including lower priority third, fourth, ... decision criteria, are selected or prioritized in a different manner. After evaluating all quality decision criteria in the assigned order, if multiple candidate positions / orientations remain selected, one configuration is randomly selected for output, or all selected configurations are output and the user is allowed to manually decide on a configuration.

[0110] In a second illustrative example, the composite quality decision criterion is constructed by a weighted sum of all quality decision criteria, such as, for example, the quality decision criteria discussed above. The user selects the quality decision criterion to be used and assigns the corresponding importance. The importance weights are relative values ​​in the interval [0,1], with 0 indicating "not relevant" and 1 indicating "maximum relevant" (although the embodiment need not be limited to this particular value or even to the normalized "relative value"). For example, to avoid bias, the quality decision criterion values ​​used as components of the composite quality decision criterion are normalized, for example, to the interval [0,1], by dividing each decision criterion value by the known maximum value. For example, the quality decision criterion value indicating the number of grid hole trajectory intersections is divided by the total number of grid holes (e.g., 169 for a typical 13x13 high dose rate brachytherapy grid template). The quality decision criterion value indicating the sum (sum) of line segment lengths is divided by the boundary dimension of the target volume, for example, the maximum length of the lesion along the insertion direction of the interventional tool. The quality decision criterion, which indicates the distance of the closest trajectory to the center of the target volume, is normalized by the length of the diagonal of the grid or by the boundary dimension of the target volume, e.g., the maximum diameter or radius of a plane parallel to the grid template. Thus, the selector selects the grid position / orientation that achieves the maximum combined decision criterion / target value and outputs this selected configuration.

[0111] The device 1 includes an output 8 for outputting the position and / or orientation selected by the position selector 6 for further evaluation and / or use by a user, for example via a user interface 7.

[0112] The device 1 includes a grid alignment evaluator 9 for receiving position signals indicative of the position and / or orientation of the physical grid template relative to the patient's body and for providing a feedback signal indicative of the position and / or orientation selected by the position selector 6.

[0113] For example, the grid alignment evaluator may visualize, e.g., in the user interface 7, the position signal and the selected position on a display device, whereby the user can position / orient the physical grid template, e.g., achieve alignment with the selected position / orientation. Thus, real-time visual guidance can be provided to aid in optimal manual grid positioning. For example, both the contour of the physical grid template based on the received position signal and the contour of the grid template corresponding to the selected position and / or orientation are displayed, e.g., using different line styles, colors, or blinking characteristics (e.g., one of them is displayed as a permanent shape, i.e., continuously, while the other is displayed as a blinking shape, i.e., intermittently). For example, the selected position is represented by a blinking contour and / or a dashed contour. For example, the contour of the current (physical) grid position is displayed in a different color than the optimal (selected) grid position. In the example shown in FIG. 4, the optimal grid position is indicated by a dashed line 41, while the current manual position is indicated by a solid line 42. When both contours match, e.g., alignment is achieved, this is indicated by an appropriate change in the presented image. For example, as shown on the right side of FIG. 4, the color of the solid line 42 changes, for example, from red (left) to green (right).

[0114] Additionally or alternatively, the feedback signal comprises an audio signal indicating a measure of the discrepancy between the selected position and the physical position. For example, the frequency of a tone or frequency of audio pulses may increase or decrease as a function of the distance (and / or to take into account differences in orientation) between the selected position and the physical position. Thus, for example, audio pulses of increasing frequency may be used to indicate approach to the optimal position.

[0115] Alternatively or additionally, the feedback signal comprises an actuator signal for controlling one or more actuators adapted to position the physical grid template, such that the device automatically controls the position of the physical grid template, e.g., places the grid template at a selected position / orientation.

[0116] The position signal is received from a tracking sensor 53, which may be, for example, an electromagnetic (EM) tracking sensor integrated into or attachable to the grid template or a carrier mechanism (e.g., a stepper device) for positioning the grid template. For example, the grid template 51 is attached to a translation and / or rotation stage 50, 52, 60, such as an ultrasonic stepper device 50, as shown in Figures 5 and 6. The ultrasonic stepper device 50 provides fine adjustment of the translation and tilt of the grid template, as shown in Figure 7. The table mount 60 (and possibly also the rotation stage 52) may enable coarse positioning and rotation of the ultrasonic stepper device 50, thereby achieving both substantial degrees of freedom in the placement and orientation of the grid template, while at the same time achieving fine movement, at least in the degrees of freedom provided by the ultrasonic stepper device.

[0117] The position signal is received from a tracking sensor 53, e.g. an EM tracking sensor, and embodiments also provide different means of determining the actual position and / or orientation of the physical grid template, for example the grid template is detected in a live imaging stream, such as fluoroscopy, optical imaging, ultrasound imaging, magnetic resonance imaging or other modality allowing acquisition of image information at a sufficiently high frequency, e.g. at least 5 images per minute, preferably at least 20 images per minute, even more preferably at least 1 image per second, ideally at least 20 images per second.

[0118] In a second aspect, the invention relates to a clinical workstation 30 including the device 1 described above. The clinical workstation is adapted to present visual information, for example an image visualization workstation. The workstation includes one or more graphic display devices, for example a monitor, a user interface device such as a keyboard and / or a mouse and / or other human interface devices known in the art, and a processor. The workstation is embodied by a computer, a smartphone, a tablet, a network server computer, and / or combinations thereof. Such a clinical workstation is suitable for or integrated into a radiology suite, a biopsy suite, an operating room suite, a radiation treatment planning and / or execution system, etc.

[0119] In addition to the device 1 and / or the clinical workstation 30, the embodiments also include, e.g. as a kit of parts, a tracking sensor 53 and / or a grid template 51 and / or a stepper device 50. The embodiments further include at least one interventional tool, e.g. a needle, catheter, ablation probe, needle loaded with one or more brachytherapy seeds, adapted for use with the grid template, e.g. having dimensions such that it fits, preferably in a sufficiently tight but slidable manner, into holes provided in the grid template.

[0120] In a third aspect, the present invention relates to a computer-implemented method for determining a position and / or orientation of a grid template relative to a human or animal body in a medical intervention procedure, said grid template comprising a plurality of holes defining a corresponding plurality of grid hole trajectories, said grid template being adapted to support and guide at least one interventional tool along such grid hole trajectory when inserted into the body via at least one of said holes in said intervention procedure.

[0121] FIG. 8 illustrates an exemplary computer-implemented method 100 according to an embodiment of the present invention.

[0122] The computer-implemented method 100 comprises: Target The method includes receiving and / or processing (101) data representative of at least one target volume within the body, and optionally data representative of at least one spatial volume at risk within the body. For example, such data may include data representative of at least one target volume within the body, Target Bo At least one segmented medical image, in which at least one of the spatial volumes at risk within the volume and / or body is represented by a corresponding segmentation label, and at least one target within the body is represented by a corresponding segmentation label. Target Bo At least one surface mesh and / or parametric spatial descriptor of at least one spatial volume at risk within the body and / or at least one medical image. For example, processing the data may include at least one target within the body. Target The method further comprises segmenting at least one medical image (even if the segmentation of the at least one medical image is instead received from an external source according to an embodiment of the present invention), such as to determine a volume and / or at least one spatial volume at risk within the body. This process comprises determining a surface mesh of the segmentation, or such a surface mesh (or an alternative descriptor) is received from an external source.

[0123] The method 100 comprises: Target The method includes generating 102 a plurality of candidate positions and / or orientations of the grid template with respect to the volume.

[0124] The method 100 comprises determining 103, for each candidate position and / or orientation, a spatial relationship between each grid hole trajectory of the grid template when located according to the candidate position and / or orientation and at least one of said target volumes, and optionally at least one of said volumes at risk, e.g., an intersection between each grid hole trajectory and at least one (or each) target volume is determined.

[0125] The method includes calculating (104) for each candidate position and / or orientation of the plurality of candidate positions and / or orientations at least one quality determination criterion representative of a suitability of the candidate position of the grid template for the interventional procedure.

[0126] The calculating (104) includes calculating at least one quality determination criterion by considering a value indicative of a geometric overlap and / or proximity of the grid hole trajectories to the at least one target volume based on the determined spatial relationship, the calculating of the at least one quality determination criterion considering a value indicative of a therapeutic effect measure of the medical intervention procedure when constrained by the determined spatial relationship.

[0127] For example, at least one quality determination criterion is calculated by considering at least a first value indicative of the degree of intersection of the grid hole trajectories for a candidate location and / or orientation with at least one target volume, the first value including a total number of grid hole trajectories intersecting with at least one target volume, and / or a total length of line segments corresponding to said intersections, and / or an average or other statistical representative value of the lengths of the line segments corresponding to said intersections.

[0128] Calculating (104) at least one quality criterion also takes into account a second value indicative of the degree of intersection of the grid hole trajectories for the candidate location and / or orientation with at least one at-risk volume, the second value including a total number of grid hole trajectories intersecting with at least one at-risk target volume, and / or a total length of line segments corresponding to the intersection with at least one at-risk volume, and / or an average or other statistical representative value of the length of the line segments corresponding to the intersection with at least one at-risk volume.

[0129] Calculating (104) the at least one quality determination criterion also takes into account a third value, the third value indicating a minimum distance of the grid hole trajectory to the center of the at least one target volume or to the center of at least one of the at least one target volumes.

[0130] Calculating (104) at least one quality decision criterion comprises calculating multiple decision criteria, such as one or more “first” values ​​(e.g., the number of intersections, total length, and / or average), and / or one or more “second” values ​​(e.g., the number of intersections, total length, and / or average), and / or a third value, and combining the multiple quality decision criteria into a composite quality decision criterion according to a weighted sum.

[0131] The method 100 comprises selecting (105) a position and / or orientation from a plurality of candidate positions and / or orientations based on at least one quality decision criterion, for example selecting (105) a position and / or orientation from a plurality of candidate positions and / or orientations that reach an extreme value of one or more quality decision criteria and / or a composite quality decision criterion. The selecting (105) also comprises a prioritized stepwise selection based on a plurality of said quality decision criteria evaluated in a predefined or configurable order (e.g., the lexicographical order method discussed above).

[0132] Thus, for a finite, discrete set of candidate grid locations and / or orientations, a direct search is performed to find the (substantially) optimal location and / or orientation that maximizes grid hole trajectory intersection with the target lesion, possibly taking into account tissues / organs at risk and / or preferentially considering a more symmetric or at least more central distribution of grid hole trajectories relative to the target lesion.

[0133] The method also includes receiving (107) a position signal indicative of a physical location and / or orientation of the physical grid template.

[0134] The method includes providing (108) a feedback signal indicative of a selected position and / or orientation and / or indicative of a physical position and / or orientation and / or a relative position and / or orientation between the physical position and / or orientation and the selected position and / or orientation.

[0135] For example, providing 108 the feedback signal may include visualizing the physical position and / or orientation and the selected position and / or orientation using a user interface. For example, the current physical position of the reading grid and the determined "optimal" position may be visualized in real time in a suitable graphic user interface (GUI) to guide the user during manual positioning of the reading grid. For example, the calculated optimal grid position may be visualized in real time on the GUI, for example, but not limited to, as a blinking rectangle, to guide the clinical professional during manual positioning of the reading grid template, and preferably the position is displayed in real time using different display styles so that both shapes are easily distinguishable.

[0136] For example, providing a feedback signal (108) may include generating an audio signal to indicate a measure of discrepancy between the selected position and / or orientation and the physical position and / or orientation.

[0137] For example, providing a feedback signal (108) may include generating an actuator signal for controlling one or more actuators adapted to position the physical grid template.

[0138] The method also comprises performing forward or inverse treatment planning of the interventional procedure using the selected position and / or orientation as input, e.g., as predetermined parameters of a planning algorithm. Such planning algorithms are well known in the art and therefore will not be discussed in further detail in the present disclosure. Similarly, a workstation according to an embodiment of the second aspect of the invention includes a treatment planning system for performing said forward or inverse treatment planning.

[0139] An embodiment of the invention also provides a method for obtaining such a (physical) grid template and measuring the position of at least one target in the body. Target Bo The present invention relates to a method for imaging a target area of ​​a body, comprising the steps of acquiring data representative of the volume, using medical imaging techniques, and executing the computer-implemented method discussed above. Thus, it is possible to provide a localized treatment, biopsy or intervention tool, particularly for target areas within the body. Target The position / orientation of the grid template is determined for use in a medical intervention procedure, such as another medical intervention performed by inserting the grid template into a vascular cavity. The grid template, e.g., a leading grid frame or a treatment grid, is a rigid or flexible grid template, as known in the art. The grid template is (without limitation) radiopaque. The grid template includes a plurality of holes, i.e., through-holes, at different locations on the grid template.

[0140] For example, such methods may include imaging a patient's body to provide said data.

[0141] For example, such a method may also include determining a physical position and / or orientation of the physical grid template, eg, using a position sensor, to provide a position signal.

[0142] Such methods may also include actuating an actuator to position the grid template based on a feedback signal, or may include manually positioning the grid template by using a feedback signal, e.g., audio / visual cues, for guidance.

[0143] Other features of the device according to the embodiments of the invention, or details of the features described above, shall be apparent in light of the description provided above with respect to the method according to the embodiments of the invention, and vice versa.

[0144] In a fourth aspect, the present invention relates to a computer program product for performing a computer-implemented method according to an embodiment of the invention, when the computer program product is executed on a suitable processor.

Claims

1. 1. A processing device for determining a position and / or orientation of a grid template relative to a human or animal body in a medical intervention procedure, said grid template including a plurality of holes defining a corresponding plurality of grid hole trajectories, said grid template supporting and guiding at least one interventional tool along said grid hole trajectory when inserted into said body through at least one of said holes in said medical intervention procedure, said processing device comprising: an anatomical spatial information processing unit for receiving and / or processing data representative of at least one target volume within the body; a generation unit for generating a plurality of candidate positions and / or orientations of the grid template relative to at least one of the target volumes within the body; - a quality calculator for calculating, for each candidate position and / or orientation of the plurality of candidate positions and / or orientations, at least one quality criterion representative of a suitability of the candidate position of the grid template for the medical intervention procedure; a position selector for selecting a position and / or orientation from said plurality of candidate positions and / or orientations based on at least one said quality determination criterion; the quality calculator determines, for each candidate position and / or orientation, a spatial relationship between each grid hole trajectory of the grid template and at least one of the target volumes when located according to the candidate position and / or orientation; A processing device, wherein at least one of the quality determination criteria, when constrained by the determined spatial relationship, includes a value indicative of a therapeutic effect measure of the medical intervention procedure, the therapeutic effect measure representing the radiation dose or ablation effect received in at least one of the target volumes when one or more radiation sources or ablators are located along the grid hole trajectory.

2. The processing device of claim 1 , wherein at least one of the quality determination criteria further comprises a value indicative of a geometric overlap and / or proximity of the grid hole trajectory to at least one of the target volumes based on the determined spatial relationship.

3. 3. The processing device of claim 1 or 2, wherein the quality calculator also determines the spatial relationship by determining, for each candidate position and / or orientation, an intersection between each grid hole trajectory of the grid template when located according to the candidate position and / or orientation and at least one of the target volumes, and calculates, for each candidate position and / or orientation, at least one of the quality determination criterion by taking into account a first value indicative of a degree of intersection of the grid hole trajectory of the candidate position and / or orientation with at least one of the target volumes.

4. 4. The processing device of claim 3, wherein the first value comprises a total number of grid hole trajectories that intersect at least one of the target volumes, and / or a total length of line segments corresponding to the intersections, and / or an average or other statistical representative of the lengths of the line segments corresponding to the intersections.

5. The anatomical spatial information processing unit further receives and / or processes data representative of at least one spatial volume within the body that is at risk; said quality calculator determines, for each candidate position and / or orientation, an intersection of each grid hole trajectory of said grid template with at least one said spatial volume at risk when located according to said candidate position and / or orientation; 5. The processing device according to claim 3 or 4, wherein the quality calculator calculates, for each candidate position and / or orientation, at least one of the quality determination criteria by taking into account at least the first value and a second value indicative of a degree of intersection of the grid hole trajectory for the candidate position and / or orientation with at least one of the spatial volumes that is at risk.

6. The processing device of claim 5, wherein the second value comprises a total number of grid hole trajectories that intersect with at least one of the spatial volumes that are exposed to danger, and / or a total length of the line segments that correspond to the intersections with at least one of the spatial volumes that are exposed to danger, and / or an average or other statistical representative value of the lengths of the line segments that correspond to the intersections with at least one of the spatial volumes that are exposed to danger.

7. 7. The processing device according to claim 3, wherein the quality calculator calculates for each candidate position and / or orientation at least one of the quality decision criteria by taking into account at least the first value and a third value, the third value indicating a minimum distance of a grid hole trajectory to a center of at least one target volume or to a center of at least one of the at least one target volumes.

8. 8. A processing device according to claim 3, wherein at least one of the quality decision criterion is a plurality of quality decision criterion, the quality calculator combining the plurality of quality decision criterion into a composite quality decision criterion according to a weighted sum, and the position selector selecting from the plurality of candidate positions and / or orientations the position and / or orientation that reaches an extreme value of the composite quality decision criterion.

9. 9. The processing device of claim 1 , wherein the at least one quality decision criterion is a plurality of quality decision criterion, and the position selector selects a first subset of the plurality of candidate positions and / or orientations based on a first quality decision criterion from the plurality of quality decision criteria, and selects at least a second subset of the first subset based on a second quality decision criterion from the plurality of quality decision criteria, the second quality decision criterion being different from the first quality decision criterion.

10. 10. The processing device of claim 1, wherein the generation unit is configured to generate the plurality of candidate positions and / or orientations by translating and / or rotating a position and / or orientation representing the grid template over a plurality of different translation and / or rotation steps.

11. 11. A processing device according to claim 1, comprising a grid template alignment evaluator for receiving a position signal indicative of a physical position and / or orientation of a physical grid template and for providing a feedback signal indicative of the physical position and / or orientation selected by the position selector and / or indicative of the physical position and / or orientation relative to the physical position and / or orientation and the selected position and / or orientation.

12. The processing device of claim 11 , wherein the grid template alignment evaluator is configured to visualize the physical position and / or orientation and the selected position and / or orientation concomitantly using a user interface.

13. 1. A computer-implemented method for determining a position and / or orientation of a grid template relative to a human or animal body in a medical intervention procedure, said grid template including a plurality of holes defining a corresponding plurality of grid hole trajectories, said grid template supporting and guiding at least one interventional tool along such grid hole trajectory when inserted into said body through at least one of said holes in said medical intervention procedure, said computer-implemented method comprising: receiving and / or processing data representative of at least one target volume within the body; generating a plurality of candidate positions and / or orientations of the grid template relative to at least one of the target volumes within the body; determining, for each candidate location and / or orientation, a spatial relationship between each grid hole trajectory of the grid template when located according to said candidate location and / or orientation and at least one of said target volumes; - calculating, for each candidate position and / or orientation of said plurality of candidate positions and / or orientations, at least one quality determination criterion representative of the suitability of said candidate position of said grid template for said medical intervention procedure when constrained by said determined spatial relationship, by taking into account at least a value indicative of a treatment effect measure of said medical intervention procedure, said treatment effect measure representative of a radiation dose or an ablation effect received in at least one said target volume when one or more radiation sources or ablators are located along said grid hole trajectories; selecting a position and / or orientation from said plurality of candidate positions and / or orientations based on at least one said quality determination criterion.

14. 14. A computer program for performing the computer implemented method of claim 13 when the computer program is run on a processor.

Images