CN107812325A - The method being tracked using motion of the magnetic resonance imaging to internal target - Google Patents

The method being tracked using motion of the magnetic resonance imaging to internal target Download PDF

Info

Publication number
CN107812325A
CN107812325A CN201710578671.XA CN201710578671A CN107812325A CN 107812325 A CN107812325 A CN 107812325A CN 201710578671 A CN201710578671 A CN 201710578671A CN 107812325 A CN107812325 A CN 107812325A
Authority
CN
China
Prior art keywords
target
image
accordance
imaging
motion
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201710578671.XA
Other languages
Chinese (zh)
Inventor
张坤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wu Quanliang
Original Assignee
Wu Quanliang
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 Wu Quanliang filed Critical Wu Quanliang
Priority to CN201710578671.XA priority Critical patent/CN107812325A/en
Publication of CN107812325A publication Critical patent/CN107812325A/en
Pending legal-status Critical Current

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00—Radiation therapy
    • A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1001—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy using radiation sources introduced into or applied onto the body; brachytherapy
    • A61N5/1007—Arrangements or means for the introduction of sources into the body
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00—Radiation therapy
    • A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/103—Treatment planning systems
    • A61N5/1039—Treatment planning systems using functional images, e.g. PET or MRI
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00—Radiation therapy
    • A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1048—Monitoring, verifying, controlling systems and methods
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00—Radiation therapy
    • A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1048—Monitoring, verifying, controlling systems and methods
    • A61N5/1049—Monitoring, verifying, controlling systems and methods for verifying the position of the patient with respect to the radiation beam
    • A61N2005/1052—Monitoring, verifying, controlling systems and methods for verifying the position of the patient with respect to the radiation beam using positron emission tomography [PET] single photon emission computer tomography [SPECT] imaging
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00—Radiation therapy
    • A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1048—Monitoring, verifying, controlling systems and methods
    • A61N5/1049—Monitoring, verifying, controlling systems and methods for verifying the position of the patient with respect to the radiation beam
    • A61N2005/1055—Monitoring, verifying, controlling systems and methods for verifying the position of the patient with respect to the radiation beam using magnetic resonance imaging [MRI]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00—Radiation therapy
    • A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N2005/1092—Details

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Pathology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Radiation-Therapy Devices (AREA)

Abstract

A kind of method that real-time tracking is carried out using motion of the magnetic resonance imaging to internal target.This method can apply in mri guided radiotherapy equipment.The flow (Fig. 5) that the flow (Fig. 1) and a realtime graphic that the present invention is tracked including a motion realized to internal target compare with texture map collection progress similarity, and the method used in two flows, including, the method (102) of target three-dimensional is partitioned into from image, the method for calculating the center (3) of objective (2), establish the method (103) of target fault texture map collection, the method (104) of real time imagery is carried out to area-of-interest, and the method (105) that realtime graphic compares with texture map collection progress similarity.

Description

The method being tracked using motion of the magnetic resonance imaging to internal target
Technical field
The present invention relates to a kind of method that motion to internal target is tracked, this method can apply to magnetic resonance It is imaged in the radiotherapy equipment of guiding.
Background technology
In three dimensions, the motion of an object (rigid body) has six-freedom degree, that is, vertically moves, move left and right, preceding After move, dive and rotate, spiral and rotation and sidewinder rotation.
Under study for action, the compound movement of object is often counted as the combination of the motion in some frees degree.To target in three-dimensional Motion in space, which is tracked, can also be reduced to be tracked its motion on six-freedom degree respectively.
Motion to internal target is tracked to be made great sense in terms of accurate radiotherapy.
At present, the radiotherapy equipment (MRI-Linac) of guiding in real time is carried out just in development using magnetic resonance imaging. MRI-Linac research and using will largely reduce the damage of normal tissue in radiation therapy process, while can be with The exposure dose to knub position is improved, to reach to the more preferable inhibition of cancer cell.
In a MRI-Linac system, realize to tumour in patient body be from motion tracking a unusual core work( Energy.There are several different research approaches at present.Such as：T.Bjerre et al are in periodical Physics in Medicine and Biology, the article on 58 (2013), texts of the L.Brix et al on periodical Medical Physics, 41 (4) (2014) Chapter, and articles of the B.Stemkens et al on periodical Physics in Medicine and Biology, 61 (2016), Each different implementation methods is elaborated respectively.
The common drawback of these methods is the movement that can only detect internal target, and can not detect the rotation of target, i.e., only Motion (vertically move, move left and right and move forward and backward) of the target on three degree of freedom in human body can be tracked.
However, the motion of the tumour in actual patient body is the motion on six-freedom degree, existing movement has rotation again, For example lung cancer patient, during breathing, tumour can carry out complicated movement and rotation with the diastole and contraction of lung.
The present invention is a kind of method that motion to internal target on six-freedom degree is tracked.
The present invention has higher tracking accuracy using unique method compared based on target fault texture.
In other correlative studys, a kind of method (US7260426) is using X-ray realtime graphics and is obtained ahead of time more Individual target area image is compared, to find the real time position of target area.The shortcomings that this method, has, and first, use is whole Individual target region is compared, and higher error probability be present, and second, can be to where target using X-ray real time imageries The patient body in region causes multiple additional radiation.
A kind of method (US9248316) is constantly carried out using in radiation therapy process to target region Ultrasound ultrasonic imagings, the position of target is then obtained with ultrasonoscopy comparison.The shortcomings that this method equally have with Track precision is low, and the rotation of target can not be detected.
The content of the invention
Disclosed by the invention is a kind of method for carrying out real-time tracking using motion of the magnetic resonance imaging to internal target.Bag Include flow (Fig. 1) and a realtime graphic and texture map collection progress phase that a motion realized to internal target is tracked The flow (Fig. 5) compared like degree, and the method used in two flows, including, target three-dimensional is partitioned into from image The method (102) of model, the method for calculating the center (3) of objective (2), the method for establishing target fault texture map collection (103) method (104) of real time imagery, is carried out to area-of-interest, and realtime graphic carries out similarity with texture map collection The method (105) of comparison.
Application claims carry out 3 D medical imaging (101) to the area-of-interest where target.Imaging method can be MRI is imaged or CT imagings, or the fusion of the imaging method imaging such as MRI, CT, PET.
After the 3-D view of the area-of-interest where obtaining target, using a kind of image partition method target Threedimensional model (2) extract (102).This dividing method can be full automatic image segmentation algorithm or hand The dynamic method for describing border, the automanual algorithm increased based on region can also be used.
Application claims use the center of gravity of objective as its center (3).
Application claims establish the atlas (103) of texture of the objective on different tomographies.Atlas it is every One figure is image of the target on a certain tomography.In the atlas established, should have target with three coordinate plane phases Faultage image in parallel plane, also to have target with the faultage image in three nonparallel planes of coordinate plane.This Invention is required in atlas, records position and angle of each image in three-dimensional system of coordinate.
Application claims are carried out continuously two-dimensional imaging (104) to a fixed position plane of target region.This The position of plane can be a plane through magnetic resonance system ISO centers (7), or adjoining imaging plane.
In image (4) formed by each frame, the faultage image (5) of patient body and the faultage image of target should be included (6)。
Similarity is carried out with the target fault image (6) in the first two field picture (4) and atlas to compare (201).According to most Similar map A position and angle information, position and target of the target's center (3) in three-dimensional system of coordinate can be extrapolated and turned Dynamic angle (202).
Similarity is carried out with the target fault image (6) in later each two field picture (4) with atlas one by one to compare (203). According to most like map B and the most like map A of previous frame relation (204), the movement locus of target is calculated.
Application claims similarity, which compares, uses " absolute error with " can also use other similarities as Evaluation criterion Evaluation criterion.
Brief description of the drawings
Fig. 1 is the flow chart for the method that the motion in the present invention to internal target is tracked.
Fig. 2 is the threedimensional model example of internal target in the present invention.
Fig. 3 A and Fig. 3 B are the preparation method examples of target section texture atlas in the present invention.Fig. 3 A are in an angle On three parallel planes on atlas, Fig. 3 B are the atlases on three parallel planes in another angle.
Fig. 4 is a kind of embodiment that real-time magnetic resonance is imaged in the present invention.
Fig. 5 is the flow chart that realtime graphic compares with texture map collection similarity in the present invention.
Embodiment
In the flow (Fig. 1) of the present invention, step 101 requires to carry out three-dimensional imaging to the area-of-interest where target. A kind of embodiment is to use magnetic resonance imaging.The parameter of magnetic resonance imaging is set so that imaging has higher space point Resolution.Meanwhile suppressed in scanning process using respiration gate control or ecg-gating due to artifact caused by physiological movement.
In the case of negligible target self-deformation (being rigid body i.e. depending on target), three-dimensional imaging only need to be at one of patient Into once as finally giving a threedimensional model of target in the physiological movement cycle.
In the case where needing to consider target self-deformation, three-dimensional imaging was needed within the physiological movement cycle of patient Repeatedly, image segmentation is then carried out, it is multiple in different time points within one physiological movement cycle of patient to obtain target Threedimensional model.For example, carrying out three-dimensional imaging to the lung of lung cancer patient, can be imaged respectively in the diastole of lung and systole phase, To obtain two kinds of different three-dimensional images, and then it is partitioned into two different object modules.
In the flow (Fig. 1) of the present invention, step 102 requires the threedimensional model of the segmentation object from image.Specific implementation In, can be with the threedimensional model as the reference for formulating radiotherapy treatment planning.
In the flow (Fig. 1) of the present invention, step 103 requires to extract faultage image from the threedimensional model of target, and builds Found an atlas.Selection to fault plane can be carried out according to the angle of radiation exposure in treatment plan.Should be tried one's best selection Image specified in irradiation plan in the plane of approximate angle.
Fig. 3 A are the examples that faultage image extraction is carried out to target.Plane 301,302 and 303 is three parallel to each other Plane, their thickness is identical, and at intervals of 0, they are angled with coordinate plane.Image in these three planes is mutual not It is identical, it is faultage image of the target in these three planes respectively.
Fig. 3 B are the examples that target fault image is extracted on another group of parallel plane 304,305 and 306.
The faultage image extracted by the method in Fig. 3 A and 3B, each has a defined location and angle Degree.In atlas, position and the angle information of each map are recorded.
, be in step 103 respectively not if having obtained multiple threedimensional models of target from step 101 and 102 Faultage image is extracted on same model.
Fig. 4 is a kind of embodiment that real-time magnetic resonance scanning (104) is carried out in the present invention.It is flat by scan patient three Piece, the image of the area-of-interest where target can be obtained.Then sweeping comprising target is selected on three plain film images Plane is retouched to be scanned.
In one embodiment, the pulse train used in real-time magnetic resonance scanning (104) is steady state free precession (SSFP), the temporal resolution of imaging is that 4 frames are per second.
In the flow (Fig. 1) of the present invention, step 105 will carry out the comparison of realtime graphic and target texture atlas.For Reduction amount of calculation, in one embodiment, alignment algorithm use the map in atlas, only with including target in image Sub-fraction region carries out registration.Because the motion of target is gradual, the map of energy Optimum Matching is in the time in atlas There is correlation in sequence, in one embodiment, by optimizing the way of search to atlas, can equally reduce algorithm Amount of calculation.
In the flow (Fig. 5) that realtime graphic and texture map collection carry out that similarity compares, when map B and map A not phases Meanwhile position and the angle information (205) of respective map are searched from atlas, the movement locus of target can be extrapolated (206)。

Claims (10)

  1. A kind of 1. method that real-time tracking is carried out using motion of the magnetic resonance imaging to internal target.Realized including one to internal The flow that the flow (Fig. 1) and a realtime graphic that the motion of target is tracked compare with texture map collection progress similarity (Fig. 5), and the method used in two flows, including, the method that target three-dimensional is partitioned into from image (102) method for, calculating the center (3) of objective (2), the method (103) of target fault texture map collection is established, it is emerging to feeling Interesting region carries out the method (104) of real time imagery, and the method that realtime graphic compares with texture map collection progress similarity (105)。
  2. 2. in accordance with the method for claim 1, the area-of-interest where internal target is imaged twice.For the first time Imaging is high-resolution three-dimensional imaging, and imaging mode can be MRI, or the medical imaging such as MRI, CT, PET mode institute Into picture fusion.Second of imaging is MRI two dimension dynamic images.
  3. 3. in accordance with the method for claim 1, to be partitioned into the threedimensional model of target.The quantity of target three-dimensional can be One, or multiple.
  4. 4. in accordance with the method for claim 1, the center of target three-dimensional is its center of gravity.
  5. 5. in accordance with the method for claim 1, to establish the atlas of image of the target on different tomographies.
  6. 6. in accordance with the method for claim 2, be the two dimension dynamics of MRI in one plane to second of target imaging into Picture.
  7. 7. in accordance with the method for claim 1, the image in target texture atlas will be with the target in two-dimentional dynamic image Carry out similarity comparison.
  8. " absolute error and " can also make as an Evaluation criterion of similarity 8. in accordance with the method for claim 7, use With other Evaluation criterions.
  9. 9. in accordance with the method for claim 7, position target in each frame two dimension dynamic image using image registration algorithm. The hunting zone that target is searched in two dimensional image can be whole image, or subregion.
  10. 10. in accordance with the method for claim 7, the method for " exhaustive search " can be used to search atlas, can also use The lookup method of optimization.
CN201710578671.XA 2017-07-17 2017-07-17 The method being tracked using motion of the magnetic resonance imaging to internal target Pending CN107812325A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710578671.XA CN107812325A (en) 2017-07-17 2017-07-17 The method being tracked using motion of the magnetic resonance imaging to internal target

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710578671.XA CN107812325A (en) 2017-07-17 2017-07-17 The method being tracked using motion of the magnetic resonance imaging to internal target

Publications (1)

Publication Number Publication Date
CN107812325A true CN107812325A (en) 2018-03-20

Family

ID=61600838

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710578671.XA Pending CN107812325A (en) 2017-07-17 2017-07-17 The method being tracked using motion of the magnetic resonance imaging to internal target

Country Status (1)

Country Link
CN (1) CN107812325A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113195051A (en) * 2018-07-28 2021-07-30 瓦里安医疗系统公司 Radiation therapy system with multiple X-ray imagers for near real time positioning

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1639675A (en) * 2001-07-06 2005-07-13 皇家菲利浦电子有限公司 Image processing method for interacting with a 3-d surface represented in a 3-d image
CN101057785A (en) * 2006-03-31 2007-10-24 美国西门子医疗解决公司 Cross reference measurement for diagnostic medical imaging
CN101120871A (en) * 2006-12-29 2008-02-13 成都川大奇林科技有限责任公司 Precise radiotherapy planning system
CN103876763A (en) * 2012-12-21 2014-06-25 重庆伟渡医疗设备股份有限公司 Image guide method implemented by aid of two-dimensional images

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1639675A (en) * 2001-07-06 2005-07-13 皇家菲利浦电子有限公司 Image processing method for interacting with a 3-d surface represented in a 3-d image
CN101057785A (en) * 2006-03-31 2007-10-24 美国西门子医疗解决公司 Cross reference measurement for diagnostic medical imaging
CN101120871A (en) * 2006-12-29 2008-02-13 成都川大奇林科技有限责任公司 Precise radiotherapy planning system
CN103876763A (en) * 2012-12-21 2014-06-25 重庆伟渡医疗设备股份有限公司 Image guide method implemented by aid of two-dimensional images

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113195051A (en) * 2018-07-28 2021-07-30 瓦里安医疗系统公司 Radiation therapy system with multiple X-ray imagers for near real time positioning
CN113195051B (en) * 2018-07-28 2023-09-12 瓦里安医疗系统公司 Radiation therapy system with multiple X-ray imagers for near real-time positioning

Similar Documents

Publication Publication Date Title
CN111161326B (en) System and method for unsupervised deep learning of deformable image registration
EP2904583B1 (en) Method for tracking three-dimensional object
CN103402453B (en) Auto-initiation and the system and method for registration for navigation system
CN105025803B (en) Segmentation from multiple 3-D views to blob
Livyatan et al. Gradient-based 2-D/3-D rigid registration of fluoroscopic X-ray to CT
JP5134957B2 (en) Dynamic tracking of moving targets
EP2646979B1 (en) Image registration apparatus
EP1516290A2 (en) Physiological model based non-rigid image registration
Kim et al. 3D reconstruction of leg bones from X-ray images using CNN-based feature analysis
Brunie et al. Pre-and intra-irradiation multimodal image registration: principles and first experiments
Chen et al. Development of Automatic Assessment Framework for Spine Deformity using Freehand 3D Ultrasound Imaging System
JP6692817B2 (en) Method and system for calculating displacement of target object
Yang et al. Ultrasound fusion image error correction using subject-specific liver motion model and automatic image registration
CN107812325A (en) The method being tracked using motion of the magnetic resonance imaging to internal target
Oulbacha et al. MRI to C‐arm spine registration through Pseudo‐3D CycleGANs with differentiable histograms
Zhang et al. An automatic ICP-based 2D-3D registration method for a high-speed biplanar videoradiography imaging system
Pan et al. A fast registration from 3D CT images to 2D X-ray images
Shanmuganathan et al. Two-Step Rigid and Non-Rigid Image Registration for the Alignment of Three-Dimensional Echocardiography Sequences From Multiple Views
KR20140054783A (en) 4d fusion display technique of motion compensated pet, and ct and/or mr
Handels et al. BVM 2009 Advances and recent developments in medical image computing
Maris et al. Deformable surface registration for breast tumors tracking: a phantom study
Feng et al. CT Image Segmentation for Preoperative Tracker Registration of Robot-Assisted Surgery
US11430203B2 (en) Computer-implemented method for registering low dimensional images with a high dimensional image, a method for training an aritificial neural network useful in finding landmarks in low dimensional images, a computer program and a system for registering low dimensional images with a high dimensional image
US11816831B2 (en) System and method for navigating and illustrating a procedure
US20240277417A1 (en) System And Method For Navigating And Illustrating A Procedure

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20180320