CN102722641A - Method for quickly obtaining dose distribution in concave tumor target region - Google Patents
Method for quickly obtaining dose distribution in concave tumor target region Download PDFInfo
- Publication number
- CN102722641A CN102722641A CN2012101660457A CN201210166045A CN102722641A CN 102722641 A CN102722641 A CN 102722641A CN 2012101660457 A CN2012101660457 A CN 2012101660457A CN 201210166045 A CN201210166045 A CN 201210166045A CN 102722641 A CN102722641 A CN 102722641A
- Authority
- CN
- China
- Prior art keywords
- tumour
- dose
- sub
- target area
- dose distribution
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 206010028980 Neoplasm Diseases 0.000 title claims abstract description 46
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000005457 optimization Methods 0.000 claims description 3
- 238000000926 separation method Methods 0.000 claims description 3
- 238000002721 intensity-modulated radiation therapy Methods 0.000 abstract description 5
- 238000004364 calculation method Methods 0.000 abstract description 3
- 238000005192 partition Methods 0.000 abstract 1
- 238000003325 tomography Methods 0.000 abstract 1
- 238000004422 calculation algorithm Methods 0.000 description 6
- 238000001959 radiotherapy Methods 0.000 description 4
- 238000005094 computer simulation Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000002068 genetic effect Effects 0.000 description 2
- 239000011521 glass Substances 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 238000002922 simulated annealing Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010845 search algorithm Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Landscapes
- Radiation-Therapy Devices (AREA)
Abstract
The invention provides a method for quickly obtaining the dose distribution in a concave tumor target region. The method comprises the following steps of: reading a computerized tomography (CT) image of a patient; drawing a tumor outline; drawing partition auxiliary lines to obtain sub-target regions; calculating and optimizing the dose distribution of each sub-target region; and adjusting the weights of the sub-target regions, and superposing the weights to obtain the overall dose distribution in a tumor target region. The calculation in the overall process is simple, convenient and quick; and a new quick solution method is provided for a treatment plan of intensity modulated radiation therapy (IMRT).
Description
Technical field
The present invention relates to the radiation therapy technology field, relate in particular to a kind of figure means of dividing that adopt and obtain the method that the tumour target dose distributes fast.
Background technology
Radiation therapy is the important means of treatment of cancer, has the advantage that result of treatment is direct, effective, can prolong patient's life for position that is difficult to perform the operation or the difficult patient who performs the operation, thereby by extensive employing.
It is the key that realizes that radiotherapy treatment planning is formulated that radiation dose distribution is calculated.
Existing radiation dose calculation method mainly contains: simulated annealing, genetic algorithm, gradient algorithm etc.Gradient algorithm belongs to deterministic method, and the speed of convergence of this method is very fast, but is absorbed in Local Extremum easily, is difficult to obtain globally optimal solution; And the random search algorithm of simulated annealing, this type of genetic algorithm can be broken away from the puzzlement of local extremum in theory, obtains globally optimal solution, but the method ordinary convergence speed of this type is slower, can not satisfy clinical requirement.So far also there is not a kind of method of Rapid Dose Calculation fast and effectively can directly solve the dose distribution that the tumour target shape is a spill.
Therefore, be necessary to provide a kind of method that obtains the distribution of spill tumour target dose fast to overcome the problems referred to above.
Summary of the invention
The object of the present invention is to provide a kind of quick, easy, effectively obtain the method for spill tumour target dose.
Correspondingly, a kind of method that spill tumour target dose distributes that obtains fast of the present invention may further comprise the steps:
S1 reads CT picture of patient;
S2, the tumour of drawing profile;
S3, the separation boost line of drawing, thus obtain sub-target area;
S4 calculates and optimizes the dose distribution of each sub-target area;
S5 adjusts the weight of sub-target area, and stack obtains tumour target area integral dose and distributes.
As further improvement of the present invention, S3 specifically comprises: confirm that straight line is divided into two parts with tumour, straight line and profile meet at 2 points, calculate the mid point of two intersection points, if mid point is positioned at inside tumor, then final dose is distributed as the addition of two parts dose distribution; If mid point is positioned at the tumour outside, then final dose is distributed as two parts dose distribution and subtracts each other.
As further improvement of the present invention, after the S3 step, judge whether sub-target area is convex, if not convex then repeats the S3 step.
As further improvement of the present invention, the weight of adjusting sub-target area among the S5 is that the sub-launched field of antithetical phrase target area carries out weight optimization.
As further improvement of the present invention, further exporting among the S5 has launched field number, the shape of each launched field and intensity.
The invention has the beneficial effects as follows: this method adopts the spill figure to be divided into the method for convex graphical; The tumour target area of spill is decomposed into the sub-target area of several convexs; Carrying out the stack of pencil beam dosage respectively calculates; At last the dose distribution of each protruding object is superposeed or subtract each other, finally obtain the dose distribution curve of spill tumour target area.Computer simulation results and actual test result are coincide.Whole process is calculated simple and efficient, for the IMRT treatment plan provides new solution fast.
Description of drawings
Fig. 1 obtains the process flow diagram of the method for spill tumour target dose distribution fast for an embodiment of the present invention is a kind of;
Fig. 2 A is the situation of mid point outside the tumour target area of boost line and profile intersection point, and final dose is distributed as sub-target dose distribution and subtracts each other;
Fig. 2 B is the situation of mid point in the tumour target area of boost line and profile intersection point, and final dose is distributed as sub-target dose distribution addition;
Fig. 3 is dose distribution test synoptic diagram.
Embodiment
Below will combine each embodiment shown in the drawings to describe the present invention.But these embodiments do not limit the present invention, and the conversion on the structure that those of ordinary skill in the art makes according to these embodiments, method or the function all is included in protection scope of the present invention.
As shown in Figure 1, a kind of method that obtains the distribution of spill tumour target dose fast of an embodiment of the present invention may further comprise the steps:
S1 reads CT picture of patient;
S2, the tumour of drawing profile;
S3, the separation boost line of drawing, thus obtain sub-target area;
S4 calculates and optimizes the dose distribution of each sub-target area;
S5 adjusts the weight of sub-target area, and stack obtains tumour target area integral dose and distributes.
Particularly, shown in Fig. 1, Fig. 2 A, Fig. 2 B, after patient's CT image obtains from hospital, use Microsoft Visual Studio software, adopt C Plus Plus to programme.Program can read patient's CT image; Operator's tumour recessed profile 1 of drawing in the drawings under the situation that the dosage line shows, manually or is automatically added boost line 2 the spill target area is divided into experimental process target area 3; Further, confirm that straight line is divided into two parts with the tumour target area.Straight line and profile meet at 2 points, calculate the mid point of two intersection points, if mid point is positioned at inside tumor, then final dose is distributed as the addition of two parts dose distribution; If mid point is positioned at the tumour outside, then final dose is distributed as two parts dose distribution and subtracts each other.Judge whether two parts are convex, if not, above segmentation procedure then repeated, till all parts are convex.After figure is cut apart; Form the sub-target area of several convexs, calculate the dose distribution of each sub-target area according to the pencil beam superposition algorithm, the sub-launched field of antithetical phrase target area carries out weight optimization; Make dose distribution reach certain requirement; Final tumour target dose is distributed as the addition of each sub-target dose distribution or subtracts each other, and the weight proportion of each sub-target area of Automatic Optimal makes net result reach some requirements.Launched field number in the dose distribution of acquisition tumour target area and the IMRT treatment plan, the position of each launched field, shape and exposure intensity.
The inventive method adopts the spill figure to be divided into the method for convex graphical; The tumour target area of spill is decomposed into the sub-target area of several convexs; Carrying out the stack of pencil beam dosage respectively calculates; At last the dose distribution of each protruding object is superposeed or subtract each other, finally obtain the dose distribution curve of spill tumour target area.Computer simulation results and actual test result are coincide.Whole process is calculated simple and efficient, for the IMRT treatment plan provides new solution fast.
As shown in Figure 3, in the dose distribution test process, the recessed profile of on patient's CT figure, drawing tumor region.Adopt the figure partitioning algorithm that the target district is decomposed into the stack of the sub-target area of convex, use existing radiotherapy in the treatment planning system program to calculate the strong result's of accent the dose distribution of dose distribution and the final tumour target area of each protruding object.Little tank 4 is adopted in test; The organic glass solid die body 5 with the human body uniform thickness is put in the centre; The centre accompanies tumour film 6; Particularly, employing is made the organic glass die body of 400*400*400 millimeter water die body (water tank) and 250*250*150 millimeter by oneself and is used corresponding K ODAK company tumour film in this embodiment.Export angle and the intensity that the result is provided with accelerator 7 according to program, shine.In this embodiment, accelerator adopts the built-in MLC X2JAW of west door subband accelerator.Obtain actual dosage line on the film, distributing with program dosage line compares.
Be to be understood that; Though this instructions is described according to embodiment; But be not that each embodiment only comprises an independently technical scheme, this narrating mode of instructions only is for clarity sake, and those skilled in the art should make instructions as a whole; Technical scheme in each embodiment also can form other embodiments that it will be appreciated by those skilled in the art that through appropriate combination.
The listed a series of detailed description of preceding text only is specifying to feasibility embodiment of the present invention; They are not in order to restriction protection scope of the present invention, allly do not break away from equivalent embodiment or the change that skill of the present invention spirit done and all should be included within protection scope of the present invention.
Claims (5)
1. one kind obtains the method that spill tumour target dose distributes fast, it is characterized in that, may further comprise the steps:
S1 reads CT picture of patient;
S2, the tumour of drawing profile;
S3, the separation boost line of drawing, thus obtain sub-target area;
S4 calculates and optimizes the dose distribution of each sub-target area;
S5 adjusts the weight of sub-target area, and stack obtains tumour target area integral dose and distributes.
2. a kind of method that spill tumour target dose distributes that obtains fast according to claim 1; It is characterized in that: S3 specifically comprises: confirm that straight line is divided into two parts with tumour; Straight line and profile meet at 2 points; Calculate the mid point of two intersection points, if mid point is positioned at inside tumor, then final dose is distributed as the addition of two parts dose distribution; If mid point is positioned at the tumour outside, then final dose is distributed as two parts dose distribution and subtracts each other.
3. a kind of method that spill tumour target dose distributes that obtains fast according to claim 1 is characterized in that: after the S3 step, judge whether sub-target area is convex, if not convex then repeats the S3 step.
4. a kind of method that spill tumour target dose distributes that obtains fast according to claim 1 is characterized in that: the weight of adjusting sub-target area among the S5 is that the sub-launched field of antithetical phrase target area carries out weight optimization.
5. a kind of method that spill tumour target dose distributes that obtains fast according to claim 1 is characterized in that: further exporting among the S5 has launched field number, the shape of each launched field and intensity.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210166045.7A CN102722641B (en) | 2012-05-25 | 2012-05-25 | Method for quickly obtaining dose distribution in concave tumor target region |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210166045.7A CN102722641B (en) | 2012-05-25 | 2012-05-25 | Method for quickly obtaining dose distribution in concave tumor target region |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102722641A true CN102722641A (en) | 2012-10-10 |
| CN102722641B CN102722641B (en) | 2015-04-08 |
Family
ID=46948398
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210166045.7A Expired - Fee Related CN102722641B (en) | 2012-05-25 | 2012-05-25 | Method for quickly obtaining dose distribution in concave tumor target region |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102722641B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104548372A (en) * | 2015-01-07 | 2015-04-29 | 上海联影医疗科技有限公司 | Radiotherapy planning method and device, radiotherapy dose determining method and device and radiotherapy quality guaranteeing method and device |
| CN105617535A (en) * | 2015-12-24 | 2016-06-01 | 上海联影医疗科技有限公司 | Dose distribution estimation method and sub-field optimization method |
| CN108171243A (en) * | 2017-12-18 | 2018-06-15 | 广州七乐康药业连锁有限公司 | A kind of medical image information recognition methods and system based on deep neural network |
| CN113877072A (en) * | 2020-07-02 | 2022-01-04 | 南京大学 | Dynamic planning algorithm for optimizing particle source distribution in brachytherapy |
| CN114904153A (en) * | 2021-02-09 | 2022-08-16 | 西安大医集团股份有限公司 | Radiotherapy plan generation method, radiotherapy plan system and storage medium |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101091818A (en) * | 2007-05-17 | 2007-12-26 | 大连现代高技术发展有限公司 | Technique for adjusting weight of beam field in real time visually in external irradiation treating tumors |
| CN101295400A (en) * | 2008-06-18 | 2008-10-29 | 山东大学 | Medical image partitioning method for radiation therapy planning system |
| US20110075806A1 (en) * | 2009-09-28 | 2011-03-31 | Janne Ilmari Nord | Treatment Plan Optimization Method For Radiosurgery |
| CN102306239A (en) * | 2011-07-22 | 2012-01-04 | 李宝生 | Method for evaluating and optimizing radiotherapy dose based on cone beam CT (Computer Tomography) image CT value correction technology |
-
2012
- 2012-05-25 CN CN201210166045.7A patent/CN102722641B/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101091818A (en) * | 2007-05-17 | 2007-12-26 | 大连现代高技术发展有限公司 | Technique for adjusting weight of beam field in real time visually in external irradiation treating tumors |
| CN101295400A (en) * | 2008-06-18 | 2008-10-29 | 山东大学 | Medical image partitioning method for radiation therapy planning system |
| US20110075806A1 (en) * | 2009-09-28 | 2011-03-31 | Janne Ilmari Nord | Treatment Plan Optimization Method For Radiosurgery |
| CN102306239A (en) * | 2011-07-22 | 2012-01-04 | 李宝生 | Method for evaluating and optimizing radiotherapy dose based on cone beam CT (Computer Tomography) image CT value correction technology |
Non-Patent Citations (2)
| Title |
|---|
| 王利波等: "sIMRT与IMRT调强放疗技术在临床中的应用", 《吉林医学》 * |
| 钱冯云等: "基于共轭梯度法的调强算法的实现", 《中国医学物理学杂志》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104548372A (en) * | 2015-01-07 | 2015-04-29 | 上海联影医疗科技有限公司 | Radiotherapy planning method and device, radiotherapy dose determining method and device and radiotherapy quality guaranteeing method and device |
| CN104548372B (en) * | 2015-01-07 | 2017-12-22 | 上海联影医疗科技有限公司 | The dosage determining device of radiotherapy |
| CN105617535A (en) * | 2015-12-24 | 2016-06-01 | 上海联影医疗科技有限公司 | Dose distribution estimation method and sub-field optimization method |
| CN108171243A (en) * | 2017-12-18 | 2018-06-15 | 广州七乐康药业连锁有限公司 | A kind of medical image information recognition methods and system based on deep neural network |
| CN113877072A (en) * | 2020-07-02 | 2022-01-04 | 南京大学 | Dynamic planning algorithm for optimizing particle source distribution in brachytherapy |
| CN114904153A (en) * | 2021-02-09 | 2022-08-16 | 西安大医集团股份有限公司 | Radiotherapy plan generation method, radiotherapy plan system and storage medium |
| WO2022170970A1 (en) * | 2021-02-09 | 2022-08-18 | 西安大医集团股份有限公司 | Method for generating radiotherapy plan, and radiotherapy plan system and storage medium |
| CN114904153B (en) * | 2021-02-09 | 2024-01-12 | 西安大医集团股份有限公司 | Method for generating radiotherapy plan, radiotherapy plan system and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102722641B (en) | 2015-04-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Nwankwo et al. | Knowledge-based radiation therapy (KBRT) treatment planning versus planning by experts: validation of a KBRT algorithm for prostate cancer treatment planning | |
| Astreinidou et al. | Adequate margins for random setup uncertainties in head-and-neck IMRT | |
| CN102722641A (en) | Method for quickly obtaining dose distribution in concave tumor target region | |
| CN102184334B (en) | Retrieval-based radiotherapy planning system and retrieval method | |
| Cao et al. | Proton energy optimization and reduction for intensity-modulated proton therapy | |
| CN110124214A (en) | Intensity modulated radiation therapy plan optimization method and application based on predicted dose distribution guidance | |
| Frank et al. | Optimal radiotherapy treatment planning governed by kinetic equations | |
| CN105893772A (en) | Data acquiring method and data acquiring device for radiotherapy plan | |
| Wu et al. | Improving proton dose calculation accuracy by using deep learning | |
| CN110327554A (en) | Intensity modulated radiation therapy plan optimization method and application based on predicted dose distribution guidance | |
| Chanrion et al. | Dosimetric consequences of pencil beam width variations in scanned beam particle therapy | |
| Engwall et al. | Fast robust optimization of proton PBS arc therapy plans using early energy layer selection and spot assignment | |
| Wolf et al. | Dosimetric impact of intrafraction rotations in stereotactic prostate radiotherapy: A subset analysis of the TROG 15.01 SPARK trial | |
| CN110349665A (en) | Carcinoma of the rectum radiotherapy planning the Automation Design method based on deep learning | |
| Carlino et al. | Commissioning of pencil beam and Monte Carlo dose engines for non-isocentric treatments in scanned proton beam therapy | |
| Jacob et al. | Clinical transition to model-based dose calculation algorithm: a retrospective analysis of high-dose-rate tandem and ring brachytherapy of the cervix | |
| CN102024097A (en) | Target spot layout optimization method in treatment plan and treatment plan system | |
| CN115038495A (en) | Method, computer program product and computer system for radiation therapy planning and radiation therapy delivery system | |
| Mackeprang et al. | Independent Monte-Carlo dose calculation for MLC based CyberKnife radiotherapy | |
| Huang et al. | Meta-optimization for fully automated radiation therapy treatment planning | |
| Wang et al. | An integrated solution of deep reinforcement learning for automatic IMRT treatment planning in non-small-cell lung cancer | |
| Li et al. | Applications of artificial intelligence in radiophysics | |
| CN113101548A (en) | Photon intensity modulated radiotherapy control method for reducing skin dose | |
| Li et al. | Deep learning-based dose map prediction for high-dose-rate brachytherapy | |
| Wang et al. | Perspectives on patient access to radiation oncology facilities and services in mainland China |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150408 Termination date: 20160525 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |