CN102722641B - 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
- CN102722641B CN102722641B CN201210166045.7A CN201210166045A CN102722641B CN 102722641 B CN102722641 B CN 102722641B CN 201210166045 A CN201210166045 A CN 201210166045A CN 102722641 B CN102722641 B CN 102722641B
- Authority
- CN
- China
- Prior art keywords
- tumor
- profiles
- sub
- dosage
- 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.)
- Expired - Fee Related
Links
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 radiation therapy technology field, particularly relate to a kind of method adopting figure programming project to obtain the distribution of tumor target dosage fast.
Background technology
Radiation therapy is the important means for the treatment of of cancer, has result for the treatment of directly, effectively, can extend the advantage of patients ' lives, be thus widely adopted for the position being difficult to perform the operation or the patient that not easily performs the operation.
It is the key realizing radiotherapy treatment planning formulation that radiation dose distribution calculates.
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 the method is very fast, but is easily absorbed in Local Extremum, is difficult to obtain globally optimal solution; And the random search algorithm that simulated annealing, genetic algorithm are this kind of, the puzzlement of local extremum can be broken away from theory, obtain globally optimal solution, but this kind of method ordinary convergence speed is comparatively slow, clinical requirement can not be met.So far also do not have a kind of method of Rapid Dose Calculation fast and effectively can directly solve tumor target shape be spill dosage distribution.
Therefore, be necessary to provide a kind of method obtaining the distribution of concave tumor target region dosage 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, method of effectively obtaining concave tumor target region dosage.
Correspondingly, a kind of method obtaining the distribution of concave tumor target region dosage fast of the present invention, comprises the following steps:
S1, reads CT picture of patient;
S2, draws tumor's profiles;
S3, determine that one meets at the separation boost line of 2 with tumor's profiles, described separation boost line and tumor's profiles form two parts, judge the point midway of described separation boost line and tumor's profiles two intersection points, if mid point is positioned at tumor's profiles inside, the dosage that then accumulated dose is distributed as two parts that this separation boost line is formed with tumor's profiles distributes and is added, if mid point is positioned at tumor's profiles outside, then accumulated dose is distributed as two parts dosage that this separation boost line and tumor's profiles formed and distributes and subtract each other, judge whether the shape of two parts is convex, if not convex then to above-mentioned be not that the part of convex adopts new separation boost line to be convex graphical to the part obtained again, thus obtain the sub-target area of convex,
S4, calculate and optimize each sub-target area dosage distribution;
S5, adjusts the weight of sub-target area, and superposition obtains the distribution of tumor target integral dose.
As a further improvement on the present invention, the sub-launched field of the weight and antithetical phrase target area that adjust sub-target area in S5 carries out weight optimization.
As a further improvement on the present invention, in S5, output has launched field number, the shape of each launched field and intensity further.
The invention has the beneficial effects as follows: the method adopts concave pattern to be divided into the method for convex graphical, the tumor target of spill is decomposed into the sub-target area of several convexs, carry out pencil beam dosage superposition calculation respectively, finally the distribution of the dosage of the sub-target area of each convex is carried out superposing or subtracting each other, finally obtain the dose profile of concave tumor target region.Computer simulation results and actual test result are coincide.Whole process computation is simple and efficient, for IMRT treatment plan provides new solution fast.
Accompanying drawing explanation
Fig. 1 is a kind of process flow diagram obtaining the method for concave tumor target region dosage distribution fast of an embodiment of the present invention;
Fig. 2 A is the situation of mid point outside tumor target of boost line and profile intersection point, and final dose is distributed as the distribution of sub-target dose and subtracts each other;
Fig. 2 B is the situation of mid point in tumor target of boost line and profile intersection point, and final dose is distributed as the distribution of sub-target dose and is added;
Fig. 3 is dosage distribution tests schematic diagram.
Embodiment
Describe the present invention below with reference to each embodiment shown in the drawings.But these embodiments do not limit the present invention, the structure that those of ordinary skill in the art makes according to these embodiments, method or conversion functionally are all included in protection scope of the present invention.
As shown in Figure 1, a kind of method obtaining the distribution of concave tumor target region dosage fast of an embodiment of the present invention, comprises the following steps:
S1, reads CT picture of patient;
S2, draws tumor's profiles;
S3, draws separation boost line, thus obtains sub-target area;
S4, calculate and optimize each sub-target area dosage distribution;
S5, adjusts the weight of sub-target area, and superposition obtains the distribution of tumor target integral dose.
Particularly, as shown in Figure 1, Figure 2 shown in A, Fig. 2 B, the CT image of patient, after hospital's acquisition, uses Microsoft Visual Studio software, adopts C Plus Plus to programme.Program can read the CT image of patient, operator draws tumour recessed profile 1 in the drawings, when dosage line shows, add boost line 2 manually or automatically and spill target area is divided into some sub-target areas 3, further, determine that one meets at the separation boost line of 2 with tumor's profiles, described separation boost line and tumor's profiles form two parts, judge the point midway of described separation boost line and tumor's profiles two intersection points, if mid point is positioned at tumor's profiles inside, the dosage that then accumulated dose is distributed as two parts that this separation boost line is formed with tumor's profiles distributes and is added, if mid point is positioned at tumor's profiles outside, then accumulated dose is distributed as two parts dosage that this separation boost line and tumor's profiles formed and distributes and subtract each other, judge whether the shape of two parts is convex, if not convex then to above-mentioned be not that the part of convex adopts new separation boost line to be convex graphical to the part obtained again.After figure segmentation, form the sub-target area of several convexs, the dosage distribution of each sub-target area is calculated according to pencil beam superposition algorithm, the sub-launched field of antithetical phrase target area carries out weight optimization, dosage is distributed and reaches certain requirement, final tumor target dosage is distributed as the addition of each sub-target dose distribution or subtracts each other, the weight proportion of Automatic Optimal each sub-target area, makes net result reach some requirements.Obtain launched field number in the dosage distribution of tumor target and IMRT treatment plan, the position of each launched field, shape and exposure intensity.
The inventive method adopts concave pattern to be divided into the method for convex graphical, the tumor target of spill is decomposed into the sub-target area of several convexs, carry out pencil beam dosage superposition calculation respectively, finally the distribution of the dosage of each convex object is carried out superposing or subtracting each other, finally obtain the dose profile of concave tumor target region.Computer simulation results and actual test result are coincide.Whole process computation is simple and efficient, for IMRT treatment plan provides new solution fast.
As shown in Figure 3, in dosage distribution tests process, the CT figure of patient draws the recessed profile of tumor region.Adopt figure partitioning algorithm target district to be decomposed into the superposition of the sub-target area of convex, use existing radiotherapy in the treatment planning system program computation to obtain the dosage distribution of the dosage distribution of each convex object and the strong result of tune of final tumor target.Test adopts little tank 4, the organic glass solid model body 5 with human body uniform thickness is put in centre, centre accompanies tumour film 6, particularly, adopt the organic glass die body of self-control 400*400*400 millimeter Water ball (water tank) and 250*250*150 millimeter in present embodiment and use corresponding K ODAK company tumour film.Angle and the intensity of accelerator 7 are set according to program Output rusults, irradiate.In present embodiment, accelerator adopts west gate subband built-in MLC X2JAW accelerator.Film obtains actual dosage line, distribute with program dosage line and compare.
Be to be understood that, although this instructions is described according to embodiment, but not each embodiment only comprises an independently technical scheme, this narrating mode of instructions is only for clarity sake, those skilled in the art should by instructions integrally, technical scheme in each embodiment also through appropriately combined, can form other embodiments that it will be appreciated by those skilled in the art that.
A series of detailed description listed is above only illustrating for feasibility embodiment of the present invention; they are also not used to limit the scope of the invention, all do not depart from the skill of the present invention equivalent implementations done of spirit or change all should be included within protection scope of the present invention.
Claims (3)
1. obtain a method for concave tumor target region dosage distribution fast, it is characterized in that, comprise the following steps:
S1, reads CT picture of patient;
S2, draws tumor's profiles;
S3, determine that one meets at the separation boost line of 2 with tumor's profiles, described separation boost line and tumor's profiles form two parts, judge the point midway of described separation boost line and tumor's profiles two intersection points, if mid point is positioned at tumor's profiles inside, the dosage that then accumulated dose is distributed as two parts that this separation boost line is formed with tumor's profiles distributes and is added, if mid point is positioned at tumor's profiles outside, then accumulated dose is distributed as two parts dosage that this separation boost line and tumor's profiles formed and distributes and subtract each other, judge whether the shape of two parts is convex, if not convex then to above-mentioned be not that the part of convex adopts new separation boost line to be convex graphical to the part obtained again, thus obtain the sub-target area of convex,
S4, calculate and optimize each sub-target area dosage distribution;
S5, adjusts the weight of sub-target area, and superposition obtains the distribution of tumor target integral dose.
2. a kind of method obtaining the distribution of concave tumor target region dosage fast according to claim 1, is characterized in that: the sub-launched field of the weight and antithetical phrase target area that adjust sub-target area in S5 carries out weight optimization.
3. a kind of method obtaining the distribution of concave tumor target region dosage fast according to claim 1, is characterized in that: exporting further in 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 CN102722641A (en) | 2012-10-10 |
| CN102722641B true 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) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
| CN108171243B (en) * | 2017-12-18 | 2021-07-30 | 广州七乐康药业连锁有限公司 | Medical image information identification method 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 |
| CN114904153B (en) * | 2021-02-09 | 2024-01-12 | 西安大医集团股份有限公司 | Method for generating radiotherapy plan, radiotherapy plan system and storage medium |
Citations (3)
| 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 |
| 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 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8009803B2 (en) * | 2009-09-28 | 2011-08-30 | Varian Medical Systems International Ag | Treatment plan optimization method for radiosurgery |
-
2012
- 2012-05-25 CN CN201210166045.7A patent/CN102722641B/en not_active Expired - Fee Related
Patent Citations (3)
| 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 |
| 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调强放疗技术在临床中的应用;王利波等;《吉林医学》;20091105;第30卷(第21期);2705～2706 * |
| 基于共轭梯度法的调强算法的实现;钱冯云等;《中国医学物理学杂志》;20081115;第25卷(第6期);867～869 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102722641A (en) | 2012-10-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106730411B (en) | Radiotherapy treatment planning optimization method and system | |
| Astreinidou et al. | Adequate margins for random setup uncertainties in head-and-neck IMRT | |
| CN102722641B (en) | Method for quickly obtaining dose distribution in concave tumor target region | |
| Gupta et al. | Image-guided radiation therapy: Physician's perspectives | |
| Mutanga et al. | Margin evaluation in the presence of deformation, rotation, and translation in prostate and entire seminal vesicle irradiation with daily marker-based setup corrections | |
| CN104888364B (en) | Ion beam radiation therapy dosage verifying method based on water equivalent coefficient | |
| CN104117151A (en) | Optimization method of online self-adaption radiotherapy plan | |
| van Asselen et al. | Intensity-modulated radiotherapy of breast cancer using direct aperture optimization | |
| CN104107062A (en) | Radiotherapeutic effect evaluation method and system | |
| CN110124214A (en) | Intensity modulated radiation therapy plan optimization method and application based on predicted dose distribution guidance | |
| CN104338240A (en) | Automatic optimization method for on-line self-adaption radiotherapy plan and device | |
| CN110415785A (en) | The method and system of artificial intelligence guidance radiotherapy planning | |
| CN105893772A (en) | Data acquiring method and data acquiring device for radiotherapy plan | |
| CN104225809A (en) | Implementation method and equipment for 4D radiotherapy plan with respiratory compensation | |
| Wolf et al. | Dosimetric impact of intrafraction rotations in stereotactic prostate radiotherapy: A subset analysis of the TROG 15.01 SPARK trial | |
| Frank et al. | Optimal radiotherapy treatment planning governed by kinetic equations | |
| CN110349665A (en) | Carcinoma of the rectum radiotherapy planning the Automation Design method based on deep learning | |
| Song et al. | Fully automatic volumetric modulated arc therapy plan generation for rectal cancer | |
| CN109513121A (en) | A kind of dosage guidance adaptive radiation therapy plan re-optimization system and method | |
| Ravkilde et al. | Fast motion-including dose error reconstruction for VMAT with and without MLC tracking | |
| 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 | |
| Biston et al. | Time of PTV is ending, robust optimization comes next | |
| Tilly et al. | Fast dose algorithm for generation of dose coverage probability for robustness analysis of fractionated radiotherapy | |
| CN110021399A (en) | A kind of automatic design method of radiotherapy treatment planning | |
| CN113101548A (en) | Photon intensity modulated radiotherapy control method for reducing skin dose |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150408 Termination date: 20160525 |