WO2013030707A1 - Integration of user inputs and correction of deformation vector field in deformable image registration workflow - Google Patents

Integration of user inputs and correction of deformation vector field in deformable image registration workflow Download PDF

Info

Publication number
WO2013030707A1
WO2013030707A1 PCT/IB2012/054169 IB2012054169W WO2013030707A1 WO 2013030707 A1 WO2013030707 A1 WO 2013030707A1 IB 2012054169 W IB2012054169 W IB 2012054169W WO 2013030707 A1 WO2013030707 A1 WO 2013030707A1
Authority
WO
WIPO (PCT)
Prior art keywords
image
dvf
contour
radiation therapy
adjustment
Prior art date
Application number
PCT/IB2012/054169
Other languages
English (en)
French (fr)
Inventor
Yogisha Mallya
Karl Antonin Bzdusek
Original Assignee
Koninklijke Philips Electronics N.V.
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Koninklijke Philips Electronics N.V. filed Critical Koninklijke Philips Electronics N.V.
Priority to EP12779160.6A priority Critical patent/EP2751778B1/en
Priority to CN201280042377.2A priority patent/CN103782320B/zh
Priority to MX2014002153A priority patent/MX2014002153A/es
Priority to US14/237,910 priority patent/US9336591B2/en
Priority to JP2014527764A priority patent/JP6053792B2/ja
Publication of WO2013030707A1 publication Critical patent/WO2013030707A1/en

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/103Treatment planning systems
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F3/00Input arrangements for transferring data to be processed into a form capable of being handled by the computer; Output arrangements for transferring data from processing unit to output unit, e.g. interface arrangements
    • G06F3/01Input arrangements or combined input and output arrangements for interaction between user and computer
    • G06F3/048Interaction techniques based on graphical user interfaces [GUI]
    • G06F3/0481Interaction techniques based on graphical user interfaces [GUI] based on specific properties of the displayed interaction object or a metaphor-based environment, e.g. interaction with desktop elements like windows or icons, or assisted by a cursor's changing behaviour or appearance
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F3/00Input arrangements for transferring data to be processed into a form capable of being handled by the computer; Output arrangements for transferring data from processing unit to output unit, e.g. interface arrangements
    • G06F3/01Input arrangements or combined input and output arrangements for interaction between user and computer
    • G06F3/048Interaction techniques based on graphical user interfaces [GUI]
    • G06F3/0484Interaction techniques based on graphical user interfaces [GUI] for the control of specific functions or operations, e.g. selecting or manipulating an object, an image or a displayed text element, setting a parameter value or selecting a range
    • G06F3/04845Interaction techniques based on graphical user interfaces [GUI] for the control of specific functions or operations, e.g. selecting or manipulating an object, an image or a displayed text element, setting a parameter value or selecting a range for image manipulation, e.g. dragging, rotation, expansion or change of colour
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/30Determination of transform parameters for the alignment of images, i.e. image registration
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2200/00Indexing scheme for image data processing or generation, in general
    • G06T2200/24Indexing scheme for image data processing or generation, in general involving graphical user interfaces [GUIs]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10081Computed x-ray tomography [CT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10088Magnetic resonance imaging [MRI]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10104Positron emission tomography [PET]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10108Single photon emission computed tomography [SPECT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20004Adaptive image processing
    • G06T2207/20012Locally adaptive
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20092Interactive image processing based on input by user
    • G06T2207/20104Interactive definition of region of interest [ROI]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30096Tumor; Lesion

Definitions

  • the following relates to the medical imaging arts. It particularly relates to image registration of successive images acquired in conjunction with the performance of fractionated radiation therapy, and is described with particular reference thereto. The following is more generally related to registration of images with previously acquired and segmented images generally.
  • the radiation dose is spread out over a series of radiation therapy sessions. Fractionation of the radiation dose provides certain benefits such as allowing the patient to recover in between sessions, and enabling medical personnel to assess the effectiveness of the radiation therapy from session to session, and to make adjustments, for example to accommodate a reduction in size over time of the malignant tumor (perhaps due to effectiveness of the radiation therapy).
  • IMRT intensity modulated radiation therapy
  • the radiation beam is controlled to deliver the radiation dose to the malignant tissue while limiting the radiation exposure of surrounding healthy tissue, especially so-called "critical" organs that may be especially susceptible to radiation damage.
  • IMRT can employ various radiation beam modulation tools such as multi-leaf collimator (MLC) apparatus, a tomographically orbiting radiation source providing irradiation over a large angular range (up to 360°), and so forth.
  • MLC multi-leaf collimator
  • a difficulty in fractionated radiation therapy is that the locations, sizes, orientations, and other aspects of the tumor and surrounding organs or tissue may change over time due to numerous factors such as weight loss or gain, natural movement of organs within the body, and so forth. If these changes are not accommodated, then the IMRT parameters may result in irradiation that is not well targeted to the malignant tumor and instead partially overlaps and irradiates critical structures.
  • the planning image may be a high resolution computed tomography (CT) image while subsequent treatment images may be lower resolution CT images (possibly acquired using a CT scanner integrated with the radiation therapy system) and/or emission images such as positron emission tomography (PET) or single photon emission computed tomography (SPECT) images.
  • PET positron emission tomography
  • SPECT single photon emission computed tomography
  • the initial planning image is manually or semi-automatically segmented to delineate the malignant tumor that is to be the target of the irradiation, and any surrounding critical structures whose irradiation should be limited. These features are delineated by contours.
  • the planning image is also used to generate radiation attenuation map. In the case of a CT planning image, this entails adjusting the CT image, which is essentially an attenuation map for the x-rays used in the CT imaging, to account for the difference in energy of the radiation therapy beam.
  • a treatment image is acquired.
  • the treatment image and the planning image are spatially registered.
  • IMRT non-rigid spatial registration is typically employed in order to precisely account for changes more complex than simple rigid translation or rotation.
  • the tumor and critical structures are contoured in the treatment image and the resulting contours are compared with the contours of the planning image so as to identify any changes.
  • contouring process is manually intensive and potentially prone to human error. Contouring can be partially automated by adapting the planning image contours to the treatment image based on the non-rigid registration deformation vector field (DVF). However, the resulting contours are sometimes not sufficiently accurate, and may require manual correction. Furthermore, the contours are usually not propagated to all aspects of the radiation therapy planning. For example, the contours are typically not used for correcting the radiation attenuation map, or for computing dose accumulation, or so forth.
  • a method comprises: computing a deformation vector field (DVF) relatively spatially registering first and second images; adapting a contour delineating a structure in the first image using the DVF to generate an initial contour for the structure in the second image; receiving a final contour for the structure in the second image; and correcting the DVF based on the initial and final contours for the structure in the second image to generate a corrected DVF; wherein the computing, adapting, and correcting are performed by an electronic processing device.
  • the correction comprises computing an adjustment DVF relating the initial and final contours and combining the DVF and the adjustment DVF to generate the corrected DVF.
  • a non-transitory storage medium stores instructions executable by an electronic processing device to perform a method comprising: adjusting a contour of a second image; and updating a deformation vector field (DVF) mapping between the second image and a first image based on the adjustment of the contour to generate an updated DVF mapping between the second image and the first image.
  • a deformation vector field DVD
  • an apparatus comprises: an electronic processing device configured to compute a deformation vector field (DVF) spatially mapping a first image and a second image; and a user interface configured to display the second image and an overlaid contour comprising an initial contour generated by mapping a contour of the first image to the second image using the DVF.
  • the user interface is further configured to receive and display user adjustments of the overlaid contour wherein the overlaid contour adjusted by said user adjustments defines a final contour.
  • the electronic processing device is further configured to correct the DVF based on the initial contour and the final contour to generate a corrected DVF.
  • the first image is a planning image for planning intensity-modulated radiation therapy (IMRT) and the second image is a treatment image for updating the IMRT
  • the electronic processing device is further configured to update one or more parameters of the IMRT based at least on the corrected DVF.
  • the apparatus further includes an IMRT delivery system configured to perform a session of the IMRT using the one or more updated parameters.
  • One advantage resides in providing automated feedback from manual contour correction back to the non-rigid spatial image registration, so that processing reliant upon accuracy of the image registration benefits from the manual contour correction.
  • Another advantage resides in providing more accurate fractionated radiation therapy by feeding contour corrections back to the non-rigid spatial registration of treatment and planning images.
  • the invention may take form in various components and arrangements of components, and in various process operations and arrangements of process operations.
  • the drawings are only for the purpose of illustrating preferred embodiments and are not to be construed as limiting the invention.
  • FIGURE 1 diagrammatically shows a fractionated intensity modulated radiation therapy (IMRT) system as disclosed herein.
  • IMRT fractionated intensity modulated radiation therapy
  • FIGURE 2 diagrammatically shows the deformation vector field (DVF) correction module in additional detail.
  • FIGURES 3-4 diagrammatically show a user interface display via which a user may input contour adjustments.
  • FIGURE 5 diagrammatically shows a process suitably performed by the contour/DVF correction module of FIGURE 2.
  • IMRT intensity-modulated radiation therapy
  • the intensity-modulated radiation therapy is performed as a fractionated radiation therapy, that is, over a series of radiation therapy sessions.
  • a planning image is acquired of the subject (e.g. oncology patient, animal undergoing veterinary radiation therapy, or so forth).
  • the planning image is typically a transmission computed tomography (CT) image, although the planning image may be acquired by another imaging modality such as magnetic resonance (MR), an emission modality such as positron emission tomography (PET) or single photon emission computed tomography (SPECT), or so forth.
  • CT transmission computed tomography
  • MR magnetic resonance
  • PET positron emission tomography
  • SPECT single photon emission computed tomography
  • the planning imaging may entail acquiring a volumetric (3D) image, a stack of planar (2D) images collectively defining a volumetric image, or so forth, and moreover the planning imaging may employ two or more different imaging modalities.
  • the planning imaging may be performed using (i) CT to acquire anatomical information and (ii) an emission modality such as PET or SPECT to acquire functional information.
  • the planning images are stored in a picture archiving and communications system (PACS) 10 from which the planning images may be accessed by authorized medical personnel such as an oncologist, radiologist, the patient's personal physician, or so forth.
  • PES picture archiving and communications system
  • the planning images are typically analyzed by an oncologist or other qualified medical personnel in order to plan the radiation therapy regimen.
  • the oncologist or other qualified medical personnel manually draw contours on or in one or more the planning images.
  • These contours delineate the malignant tumor (or, more generally, malignant tissue that is to be the target of the radiation therapy) and typically also delineate one or more "critical" structures, which are organs, tissue, or so forth that is susceptible to substantial damage from radiation exposure.
  • the goal of the IMRT is to deliver at least a specified radiation dosage to the target tumor or tissue while keeping the radiation exposure of critical structures below a specified maximum allowable level.
  • These contours are stored with the planning image in the PACS 10.
  • the radiation therapy planning is based on information such as: the contours delineating the tumor and critical structures; radiation dosage thresholds for the tumor (usually a minimum dosage threshold) and for critical structures (usually a maximum allowable exposure); and a radiation attenuation map for the subject.
  • the latter is typically generated from the planning image.
  • the attenuation map can be generated by adjusting the attenuation values of the CT image to account for the difference in photon energy between the x-ray radiation used in the CT imaging and the radiation used in the radiation therapy.
  • the parameters may include multi-leaf collimator (MLC) apparatus settings, rotation speed and/or trajectory of a tomographic radiation source, total beam power, and so forth.
  • MLC multi-leaf collimator
  • Settings or other values for the various parameters of the radiation therapy delivery system that satisfy the dosage thresholds for the tumor and critical organs are determined for each radiation therapy session of the fractionated IMRT using an automated (or "inverse") planning procedure.
  • the planning process can take a relatively long time.
  • the planning images are typically acquired with high resolution, and contouring may be performed in each slice of a stack of slices in order to delineate the tumor and critical structures three-dimensionally.
  • the inverse planning although performed by a computer or other electronic processing device, is complex and can take a substantial amount of computing time.
  • the planning images are typically acquired some time (e.g., hours, days, or longer) before commencement of the first radiation therapy session. During this interval, changes may occur, such as organ movement due to expansion, contraction, formation, or dissolution of air pockets in the subject, changes in patient weight, growth or shrinkage of the malignant tumor, or so forth.
  • a treatment image acquisition system 12 is employed to acquire a treatment image 14 of the subject a short time prior to commencement of the radiation therapy session.
  • the treatment image acquisition system 12 is typically a CT imaging scanner, although an MR or another imaging modality is also contemplated.
  • the treatment image 14 is spatially registered with a corresponding planning image 16 retrieved from the PACS 10.
  • the spatial registration is performed by a deformable image registration module 20 using a deformable image registration algorithm such as the Demons deformable registration technique or another deformable image registration technique such as level sets, B-spline, or so forth.
  • the output of the deformable image registration is a deformation vector field (DVF) 22 relatively spatially registering the planning image 16 and the treatment image 14.
  • the DVF 22 is also known in the art by other, similar nomenclature, such as the displacement field, and the term deformation vector field or DVF as used herein is intended to encompass these various nomenclatures.
  • the DVF 22 relatively spatially registers the planning and treatment images by specifying for pixels or voxels of one image the displacement transformation (including both distance and direction) in order to align with the second image.
  • the DVF 22 indicates the displacement transformations to align the planning image with the treatment image, or indicates the displacement transformations to align the treatment image with the planning image.
  • Non-rigid deformation (or non-rigid spatial registration, or similar phraseology) is obtained because the pixel or voxel displacement transformations can vary across the image to provide non-rigid deformation.
  • the DVF 22 maps one image (e.g., the planning image) to the other image (e.g., the treatment image).
  • mapping in the opposite direction can also be performed using the DVF 22 by suitable reversal of the pixel or voxel vector transformations.
  • the DVF 22 spatially maps a first image (e.g., the planning image 16) and a second image (e.g., the treatment image 14).
  • the DVF 22 can be used to update radiation therapy planning parameters.
  • the DVF 22 can be used to warp the dosage, and/or to warp the attenuation map, in order to account for changes such as movement, shrinkage, or expansion of the tumor and/or various critical organs.
  • the DVF 22 can be used to adjust planning image contours 26 associated with the planning image 16.
  • the planning image 16 along with the contours 26 and other metadata such as acquisition time, imaging parameters, and so forth may be stored in a file 28 denoted diagrammatically in FIGURE 1 by a dotted box surrounding the components 16, 26.
  • the contours 26 are generated by the oncologist by manually delineating the contours 26.
  • the contouring of the planning image 16 may be assisted by an automated contouring algorithm (e.g., semi-automated contouring).
  • Contouring is generally difficult, because the boundaries different soft-tissue organs can be difficult to detect by a trained human observer or by an automated segmentation algorithm.
  • the difficulty is enhanced by poor soft-tissue contrast in CT scans.
  • the oncologist or other medical personnel perform the tedious task of contouring the tumor and critical structures in the planning image. It is preferable to avoid repeating such tedium each time a new treatment image is acquired to track changes in the patient.
  • One approach for simplifying the contouring of the treatment image 14 is to use the planning contours 26 as patient-specific priors, and to use the DVF 22 to map the planning image contours 26 from the planning image 16 to the treatment image 14.
  • This approach is known in the art as region of interest (ROI) warping or propagation.
  • ROI region of interest
  • the warped contours provide automatic delineation of the target tumor and critical structures in the treatment image 14, without performing tedious manual contouring.
  • any mapping error in the DVF 22 (which may be caused, for example, by severe deformation which is difficult to map accurately, or by image artifacts, or by features that appear in only one image, differences in contrast agent distribution between the planning and treatment images, or so forth). Any errors in the DVF 22 that are in the proximity of the contour undergoing mapping will be propagated into the warped contour. Since these contours are used for patient-critical operations such as defining radiation dose distributions that target the malignant tumor while substantially avoiding critical structures, any contour error can have substantial adverse consequences for the subject undergoing IMRT.
  • One way to address errors in the contour mapping is to allow the user (e.g., oncologist or radiologist) to view the contour superimposed on the treatment image 14 and, if errors are observed, to manually edit the contour using a mouse or other input device.
  • this approach corrects only the contour, but not the DVF 22 that was the source of the erroneous contour. This can be a problem because some aspects of the radiation therapy, such as the attenuation map, dose accumulation map, and so forth, depend directly on the DVF 22 rather than on the contours.
  • improved user correction is provided by a contour/DVF correction module 30 which enables the user to correct contours and also propagates any such corrections back into the DVF 22 in order to create a corrected DVF 32.
  • the adjusted or corrected DVF 32 (rather than the original DVF 22) is then input to a radiation therapy session preparation module 34 for use in operations such as dose warping performed by a dose warping sub-module 36, dose accumulation mapping performed by a dose accumulation tracking sub-module 38, or so forth.
  • the resulting IMRT session plan is suitably executed by an IMRT delivery system 40, which may by way of illustrative example comprise a linear electron accelerator (LINAC), configured to perform tomographic irradiation, or a multi-beam radiation therapy apparatus, or so forth.
  • LINAC linear electron accelerator
  • the electronic processing components 20, 30, 34 can be variously embodied by one or more electronic processing devices, such as an illustrative computer 42.
  • the contour/DVF correction module 30 preferably includes user interfacing components such as an illustrative display device 44 and one or more user input devices such as an illustrative keyboard 46 and mouse 48.
  • a non-transitory storage medium stores instructions executable by an electronic processing device (such as the illustrative computer 42) to perform methods performed by the electronic processing components 20, 30, 34.
  • Such a non-transitory storage medium may comprise, by way of illustrative example: a hard disk or other magnetic storage medium; an optical disk or other optical storage medium; an electrostatic memory such as a FLASH memory; a read-only memory (ROM); a random access memory (RAM); or so forth.
  • contour/DVF correction module 30 receives as inputs: (1) the planning image contours 26; the DVF 22; and the treatment image 14.
  • An image contours warping module 50 maps the planning image contours 26 to the treatment image 14 using the DVF 22 to automatically generate initial contours 52 for the treatment image 14.
  • a contours editing module 54 which includes a graphical user interface (GUI), displays the treatment image 14 and an overlaid contour comprising the initial contour 52.
  • GUI graphical user interface
  • FIGURE 3 shows a treatment image I T (corresponding to the treatment image 14 of FIGURES 1 and 2) displayed on the display device 44, along with initial contours Cu, Ci2, Ci3 corresponding to the automatically generated initial contours 52 of FIGURE 2 shown as overlaid contours that are overlaid on the treatment image I T displayed on the display device 44.
  • the contours editing module 54 also provides instructions to the user for editing the contours Cu, Q 2 , Q 3 , e.g. the displayed text "Adjust contours then (CONTINUE)". It is noted that the contours Cu, Q 2 , Q 3 , are close to, but not precisely aligned with, apparent image features. The misalignment is due to some errors in the DVF 22 used to map the planning image contours to the treatment image 14.
  • the user suitably uses an on-screen pointer P controlled by the mouse 48 (or by a trackball or other pointing-type user input device) to adjust the initial contours Cu, Q 2 , Ci3 by making various user adjustments to create final contours C FI , C F2 , C F3 shown in FIGURE 4.
  • a "select, drag-and-drop” process can be used, in which the user points to a portion of a contour, clicks and holds a mouse button to select that portion, drags it to a new position, and releases the mouse button to "drop" the contour portion into its new location.
  • the process can be repeated to make a plurality of user adjustments, and in some embodiments the overlaid contours shown on the display device 44 are updated after each such user adjustment to reflect the updated contour.
  • the overlaid contours shown on the display device 44 are updated after each such user adjustment to reflect the updated contour.
  • the process can be repeated to make a plurality of user adjustments, and in some embodiments the overlaid contours shown on the display device 44 are updated after each such user adjustment to reflect the updated contour.
  • the process can be repeated to make a plurality of user adjustments, and in some embodiments the overlaid contours shown on the display device 44 are updated after each such user adjustment to reflect the updated contour.
  • the final C F i, C F2 , C F3 shown in FIGURE 4 constitute user-adjusted final contours 56 in FIGURE 2.
  • These final contours 56 are corrected for errors in the DVF 22; however, the DVF 22 does not itself include the user adjustments embodied in the final contours 56.
  • a deformable contours registration module 60 receives the initial contours 52 and the final contours 56, and computes an adjustment deformation vector field (adjustment DVF) 62 that accounts for the user adjustments.
  • the deformable contours registration module 60 uses a deformable image registration algorithm such as the Demons deformable registration technique or another deformable image registration technique such as level sets, B-spline, or so forth to generate the adjustment DVF 62 that maps the initial contour 52 to the final contour 56 (or vice versa).
  • a DVF combiner module 64 combines the DVF 22 (which does not have correction) and the adjustment DVF 62 to generate the corrected DVF 32.
  • the combiner module 64 suitably operates by adding the adjustment DVF 62 to the DVF 22 to generate the corrected DVF 32.
  • This addition operation assumes that both the DVF 22 and the adjustment DVF 62 have the same "sign", i.e. map in the same direction (e.g., mapping the planning image to the treatment image and the initial contour to the final contour). If the DVF 22 and the adjustment DVF 62 have opposite "sign" (e.g., the DVF 22 maps the planning image to the treatment image while the adjustment DVF 62 maps the final contour to the initial contour), then the combiner module 64 suitably operates by subtracting the adjustment DVF 62 from the DVF 22 to generate the corrected DVF 32.
  • contour/DVF correction module 30 a process suitably performed by the contour/DVF correction module 30 is described.
  • the contours 26 of the planning image 16 are warped in accordance with the DVF 22 to automatically generate the initial contours 52.
  • a user employs a graphical user interface (GUI) to edit the initial contours 52 in order to generate user-adjusted final contours 56.
  • GUI graphical user interface
  • the adjustment DVF 62 is generated based on differences between the initial contours 52 and the final contours 56.1n an operation 106 performed by the DVF combiner module 64, the DVF 22 is updated with the adjustment DVF 62 to generate the corrected DVF 32.
  • the final contour 56 delineating the structure in the treatment image 14 is generated by user adjustment of the initial contour 52.
  • the final contour may be generated by performing an automated optimization of the mapped initial contour 52 respective to the structure in the treatment image 14 to generate the final second image contour delineating the structure in the treatment image.
  • contour/DVF adjustment has been disclosed herein in the context of IMRT, the contour/DVF adjustment is suitably applied to other applications in which a first image (e.g., the planning image 16 in the illustrative embodiment) is contoured to generate first image contours (e.g., contours 26 in the illustrative embodiment), a second image is acquired (e.g., the treatment image 14 in the illustrative embodiment), and non-rigid image registration is applied to generate a DVF mapping between the first and second images.
  • a first image e.g., the planning image 16 in the illustrative embodiment
  • first image contours e.g., contours 26 in the illustrative embodiment
  • second image e.g., the treatment image 14 in the illustrative embodiment
  • non-rigid image registration is applied to generate a DVF mapping between the first and second images.
  • the first image may be a reference image acquired of a laboratory test animal prior to commencement of a pre-clinical test, and the second image may be subsequent images of the laboratory test animal acquired at various points during the pre-clinical test.
  • the first image may be a reference image acquired of a human subject prior to commencement of a clinical trial, and the second image may be subsequent images of the human subject acquired at various points during the clinical trial.
  • Other contemplated applications include various radiation therapy techniques such as IMRT and specific techniques such as conformal arc therapy or volumetric-modulated arc therapy (VMAT), various pre-operative surgical planning procedures, interventional imaging or image guided surgery, post-operative assessment, and so forth, each of which readily benefit from contour/DVF adjustment as disclosed herein.

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • General Health & Medical Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Human Computer Interaction (AREA)
  • Public Health (AREA)
  • Pathology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • Medical Informatics (AREA)
  • Quality & Reliability (AREA)
  • Radiation-Therapy Devices (AREA)
  • Collating Specific Patterns (AREA)
  • Image Processing (AREA)
PCT/IB2012/054169 2011-08-30 2012-08-16 Integration of user inputs and correction of deformation vector field in deformable image registration workflow WO2013030707A1 (en)

Priority Applications (5)

Application Number Priority Date Filing Date Title
EP12779160.6A EP2751778B1 (en) 2011-08-30 2012-08-16 Integration of user inputs and correction of deformation vector field in deformable image registration workflow
CN201280042377.2A CN103782320B (zh) 2011-08-30 2012-08-16 在可变形图像配准工作流中用户输入和变形矢量场的校正的集成
MX2014002153A MX2014002153A (es) 2011-08-30 2012-08-16 Integracion de entradas de usuario y correcion de campo de vector de deformacion en la dinamica de trabajo del registro deformable de imagenes.
US14/237,910 US9336591B2 (en) 2011-08-30 2012-08-16 Integration of user inputs and correction of deformation vector field in deformable image registration workflow
JP2014527764A JP6053792B2 (ja) 2011-08-30 2012-08-16 可変画像レジストレーションのワークフローにおけるユーザー入力と変形ベクトル場の修正との統合

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US201161528790P 2011-08-30 2011-08-30
US61/528,790 2011-08-30

Publications (1)

Publication Number Publication Date
WO2013030707A1 true WO2013030707A1 (en) 2013-03-07

Family

ID=47089092

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/IB2012/054169 WO2013030707A1 (en) 2011-08-30 2012-08-16 Integration of user inputs and correction of deformation vector field in deformable image registration workflow

Country Status (6)

Country Link
US (1) US9336591B2 (es)
EP (1) EP2751778B1 (es)
JP (1) JP6053792B2 (es)
CN (1) CN103782320B (es)
MX (1) MX2014002153A (es)
WO (1) WO2013030707A1 (es)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3498335A1 (en) * 2017-12-18 2019-06-19 Koninklijke Philips N.V. Evaluation of an anatomic structure with respect to a dose distribution in radiation therapy planning
CN113577577A (zh) * 2017-09-07 2021-11-02 医科达有限公司 自适应放射疗法系统

Families Citing this family (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9401051B2 (en) * 2010-06-23 2016-07-26 Varian Medical Systems International Ag Mechanism for dynamically propagating real-time alterations of medical images
CN103119626B (zh) * 2010-09-17 2016-08-03 皇家飞利浦电子股份有限公司 用于辐射治疗规划的具有实时轮廓片段影响绘制的轮廓勾画
GB2495150B (en) * 2011-09-30 2015-07-01 Mirada Medical Method and system of defining a region of interest on medical scan images
WO2014155346A2 (en) * 2013-03-29 2014-10-02 Koninklijke Philips N.V. Image registration
US10025479B2 (en) * 2013-09-25 2018-07-17 Terarecon, Inc. Advanced medical image processing wizard
US9774838B2 (en) * 2015-06-12 2017-09-26 Accuray Incorporated Ambient light suppression using color space information to derive pixel-wise attenuation factors
CN105260997B (zh) * 2015-09-22 2019-02-01 北京医拍智能科技有限公司 一种自动获取目标图像的方法
JP6164662B2 (ja) * 2015-11-18 2017-07-19 みずほ情報総研株式会社 治療支援システム、治療支援システムの動作方法及び治療支援プログラム
EP3181049B1 (en) * 2015-12-18 2018-02-14 RaySearch Laboratories AB Radiotherapy method, computer program and computer system
JP6732593B2 (ja) * 2016-08-03 2020-07-29 キヤノン株式会社 画像処理装置および画像処理方法
US10713801B2 (en) * 2017-01-06 2020-07-14 Accuray Incorporated Image registration of treatment planning image, intrafraction 3D image, and intrafraction 2D x-ray image
US11449208B2 (en) * 2017-07-06 2022-09-20 Varian Medical Systems International Ag Interactive and intuitive method to shape 3D dose distribution during optimization of IMRT plans
CN107224678B (zh) * 2017-07-17 2020-04-10 上海联影医疗科技有限公司 一种治疗评估系统及存储介质
JP7432329B2 (ja) * 2018-09-13 2024-02-16 キヤノンメディカルシステムズ株式会社 医用画像診断システム
US10918885B2 (en) * 2018-09-27 2021-02-16 Varian Medical Systems International Ag Systems, methods and devices for automated target volume generation
US11307730B2 (en) 2018-10-19 2022-04-19 Wen-Chieh Geoffrey Lee Pervasive 3D graphical user interface configured for machine learning
US11216150B2 (en) 2019-06-28 2022-01-04 Wen-Chieh Geoffrey Lee Pervasive 3D graphical user interface with vector field functionality
EP3789085A1 (en) * 2019-09-05 2021-03-10 Koninklijke Philips N.V. Dose-guided deformable image registration
CN111524081B (zh) * 2020-04-24 2023-10-10 讯飞医疗科技股份有限公司 肺部影像角度矫正方法、装置、电子设备和存储介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007014105A2 (en) * 2005-07-22 2007-02-01 Tomotherapy Incorporated Method and system for adapting a radiation therapy treatment plan based on a biological model
WO2007014092A2 (en) * 2005-07-22 2007-02-01 Tomotherapy Incorporated Method of placing constraints on a deformation map and system for implementing same
US20100232572A1 (en) * 2009-03-11 2010-09-16 Varian Medical Systems International Ag Use of planning atlas in radiation therapy
WO2010148250A2 (en) * 2009-06-17 2010-12-23 Tomotherapy Incorporated System and method of applying anatomically-constrained deformation

Family Cites Families (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6195480A (ja) * 1984-10-17 1986-05-14 Hitachi Ltd 画像間位置合わせ方式
DE69332042T2 (de) * 1992-12-18 2003-01-02 Koninkl Philips Electronics Nv Ortungszurückstellung von relativ elastisch verformten räumlichen Bildern durch übereinstimmende Flächen
CN100586507C (zh) * 2004-07-23 2010-02-03 吴大可 三维适形近距离放射治疗集成系统
WO2006018761A1 (en) * 2004-08-13 2006-02-23 Koninklijke Philips Electronics N.V. Radiotherapeutic treatment plan adaptation
JP4651375B2 (ja) * 2004-12-16 2011-03-16 株式会社日立メディコ 医用画像表示装置及びその方法
US7352370B2 (en) * 2005-06-02 2008-04-01 Accuray Incorporated Four-dimensional volume of interest
US7596207B2 (en) * 2005-07-14 2009-09-29 Koninklijke Philips Electronics N.V. Method of accounting for tumor motion in radiotherapy treatment
JP2009502250A (ja) * 2005-07-22 2009-01-29 トモセラピー・インコーポレーテッド 放射線療法治療計画に関連するデータを処理するための方法およびシステム
US20070116381A1 (en) * 2005-10-19 2007-05-24 Ali Khamene Method for deformable registration of images
EP2018627B1 (en) * 2006-05-11 2017-03-22 Koninklijke Philips N.V. Deformable registration of images for image guided radiation therapy
US7945117B2 (en) * 2006-08-22 2011-05-17 Siemens Medical Solutions Usa, Inc. Methods and systems for registration of images
WO2009016530A2 (en) * 2007-07-27 2009-02-05 Koninklijke Philips Electronics N.V. Interactive atlas to image registration
US8265356B2 (en) 2008-01-30 2012-09-11 Computerized Medical Systems, Inc. Method and apparatus for efficient automated re-contouring of four-dimensional medical imagery using surface displacement fields
WO2010086776A1 (en) 2009-01-30 2010-08-05 Koninklijke Philips Electronics N.V. System for providing lung ventilation information
US9245336B2 (en) 2010-12-15 2016-01-26 Koninklijke Philips N.V. Contour guided deformable image registration
US8526692B2 (en) * 2011-06-30 2013-09-03 Wisconsin Alumni Research Foundation Reduction of transitivity errors in radiotherapy image registration

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007014105A2 (en) * 2005-07-22 2007-02-01 Tomotherapy Incorporated Method and system for adapting a radiation therapy treatment plan based on a biological model
WO2007014092A2 (en) * 2005-07-22 2007-02-01 Tomotherapy Incorporated Method of placing constraints on a deformation map and system for implementing same
US20100232572A1 (en) * 2009-03-11 2010-09-16 Varian Medical Systems International Ag Use of planning atlas in radiation therapy
WO2010148250A2 (en) * 2009-06-17 2010-12-23 Tomotherapy Incorporated System and method of applying anatomically-constrained deformation

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
BURAK EREM ET AL: "Interactive Deformable Registration Visualization and Analysis of 4D Computed Tomography", 4 January 2008, MEDICAL BIOMETRICS; [LECTURE NOTES IN COMPUTER SCIENCE], SPRINGER BERLIN HEIDELBERG, BERLIN, HEIDELBERG, PAGE(S) 232 - 239, ISBN: 978-3-540-77410-5, XP019137069 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113577577A (zh) * 2017-09-07 2021-11-02 医科达有限公司 自适应放射疗法系统
EP3498335A1 (en) * 2017-12-18 2019-06-19 Koninklijke Philips N.V. Evaluation of an anatomic structure with respect to a dose distribution in radiation therapy planning
WO2019121436A1 (en) * 2017-12-18 2019-06-27 Koninklijke Philips N.V. Evaluation of an anatomic structure with respect to a dose distribution in radiation therapy planning
US11383103B2 (en) 2017-12-18 2022-07-12 Koninklijke Philips N.V. Evaluation of an anatomic structure with respect to a dose distribution in radiation therapy planning

Also Published As

Publication number Publication date
CN103782320B (zh) 2017-03-15
EP2751778A1 (en) 2014-07-09
MX2014002153A (es) 2014-05-01
CN103782320A (zh) 2014-05-07
JP6053792B2 (ja) 2016-12-27
US9336591B2 (en) 2016-05-10
US20140201670A1 (en) 2014-07-17
JP2014531634A (ja) 2014-11-27
EP2751778B1 (en) 2017-10-11

Similar Documents

Publication Publication Date Title
US9336591B2 (en) Integration of user inputs and correction of deformation vector field in deformable image registration workflow
AU2018204655B2 (en) Image guidance for radiation therapy
EP1778074B1 (en) Weighted surface-to-surface mapping
US8660325B2 (en) Efficient user interaction with polygonal meshes for medical image segmentation
EP2175931B1 (en) Systems for compensating for changes in anatomy of radiotherapy patients
US8929621B2 (en) Methods and systems for segmentation and surface matching
US10149987B2 (en) Method and system for generating synthetic electron density information for dose calculations based on MRI
US8260013B2 (en) Data representation for rtp
JP5330992B2 (ja) 生物学に導かれた適応的な治療計画
JP2015083068A (ja) 放射線治療装置およびシステムおよび方法
US11273327B2 (en) Virtual beam's-eye view imaging in radiation therapy for patient setup
WO2012069965A1 (en) Interactive deformation map corrections
CN104636618B (zh) 放疗计划系统
Penney Registration of tomographic images to X-ray projections for use in image guided interventions
US9403035B2 (en) Dose-optimised patient positioning for radiotherapy
WO2022009014A1 (en) System for planning and verifying treatment during iort procedures
CN102908725A (zh) 放射治疗中的造影剂强化成像
Hoffmann et al. Correction of prostate misalignment in radiation therapy using US-CT registration

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 12779160

Country of ref document: EP

Kind code of ref document: A1

REEP Request for entry into the european phase

Ref document number: 2012779160

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 14237910

Country of ref document: US

ENP Entry into the national phase

Ref document number: 2014527764

Country of ref document: JP

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: MX/A/2014/002153

Country of ref document: MX

NENP Non-entry into the national phase

Ref country code: DE