CN102488971B - The image imaging method of real-time dynamic proton imaging and radiotherapy and the imaging system of use the method - Google Patents
The image imaging method of real-time dynamic proton imaging and radiotherapy and the imaging system of use the method Download PDFInfo
- Publication number
- CN102488971B CN102488971B CN201110328222.2A CN201110328222A CN102488971B CN 102488971 B CN102488971 B CN 102488971B CN 201110328222 A CN201110328222 A CN 201110328222A CN 102488971 B CN102488971 B CN 102488971B
- Authority
- CN
- China
- Prior art keywords
- proton
- imaging
- image
- proton radiation
- radiotherapy
- 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
Landscapes
- Radiation-Therapy Devices (AREA)
Abstract
Real-time dynamic proton imaging designed by the present invention and the image imaging method of radiotherapy, it is the energy being gathered proton radiation by the TFT Amorphous silicon flat-panel detectors of band cesium iodide coating, transformed by digital-to-analogue again and obtain image, it is characterized in that the energy acquisition carrying out twice proton radiation, it is before proton radiation is injected and is detected object that first time gathers, second time collection be at proton radiation after being detected object; Then the data collected for twice by process form image.The feature of this method is, by the deviate of comparison twice image data, and then the impact that reduction proton radiation brings imaging through testee generation scattering.
Description
Technical field
The present invention relates to a kind of formation method, the imaging system of the image imaging method of particularly a kind of real-time dynamic proton imaging and radiotherapy and image reconstruction optimization method and use the method.
Background technology
Proton radiation therapy technology is a kind of new anti-cancer therapies.Concentrate due to energy major part and be deposited on Bragg peak, compared with proton radiotherapy is treated with traditional photon mode, have and significantly treat accuracy benefits: tumor exposure dose is high; Low to the proximal tissue radiation dose between incident top layer to tumor; The radiationless dosage of remote organization of contrast tumor deep layer.High accuracy in order to ensure proton radiotherapy changes into high curative effect, needs the Two dimensional Distribution of Real-Time Monitoring proton beam energies whether to meet the requirement for the treatment of planning systems to precision.Along with proton radiation therapy cancer reaches its maturity, in clinical practice, the demand of proton detector is also increased gradually.Proton is charged particle, therefore to penetrate the uncharged X-ray of scattering ratio in human body and detector process remarkable.The resolution of the diffuse transmission influence image of proton, therefore needs to revise to be applicable to clinical practice to scattering track in image reconstruction process, as proton images guides radiotherapy etc.
The methods such as existing proton imaging technology mainly uses monolayer film, monolayer ionization chamber matrix-scanning.Film Dosimeter program is complicated, cannot be applied to real time image collection and process.Ionization chamber scan method uses method comparatively widely in current clinical practice, and its defect is that picture position resolution is poor.And above-mentioned single layer image system can only measure the positional information after proton injection object under test, therefore effectively cannot revise scattering of proton and also cannot optimize raising resolution.
Summary of the invention
The object of the invention is the deficiency in order to solve above-mentioned technology and provide one can implement image acquisition and processing, and effectively revise scattering of proton, obtain image resolution ratio high, the imaging system of the image imaging method of clear reliable real-time dynamic proton imaging and radiotherapy and image reconstruction optimization method and use the method.
In order to achieve the above object, real-time dynamic proton imaging designed by the present invention and the image imaging method of radiotherapy, it is the energy being gathered proton radiation by the TFT Amorphous silicon flat-panel detectors of band cesium iodide coating, transformed by digital-to-analogue again and obtain image, it is characterized in that the energy acquisition carrying out twice proton radiation, it is before proton radiation is injected and is detected object that first time gathers, second time collection be at proton radiation after being detected object; Then the data collected for twice by process form image.The feature of this method is, by the deviate of comparison twice image data, and then the impact that reduction proton radiation brings imaging through testee generation scattering.
Further method is, described proton radiation source is the proton radiation source of movement, proton radiation is injected from three different angles and is detected object, and the described position of TFT Amorphous silicon flat-panel detectors coordinates with the exposure pathways of proton radiation, and gather the proton radiation energy that different angles inject, finally obtain the image information of testee different depth layer.
The image reconstruction optimization method being applicable to the image imaging method of real-time dynamic proton imaging and radiotherapy designed by the present invention, utilize the proton radiation that collects to inject energy attenuation rate before and after testee and positional information, reducing scattering of proton to the impact of imaging by calculating, finally obtaining two-dimentional water equivalence range (thickness) distributed image that proton beam injects Energy distribution two dimensional image before and after testee and object under test; It is calculated the maximum of probability track of proton by Monte Carlo simulation, and the maximum of probability track calculated the most at last is converted into two dimensional image that described calculating reduces the impact of scattering of proton on imaging.Calculation procedure and the formula of wherein said maximum of probability track are as follows:
The degree of depth (x
1) maximum of probability track point:
Wherein:
Wherein: wherein: assuming that incidence point coordinate is (0,0), x is eye point abscissa, and y is eye point vertical coordinate; θ is shooting angle, y
1with θ
1be respectively particle to be asked at degree of depth x
1the vertical coordinate at place and scattering angle; X
0for the radiant intensity constant relevant to passing material; β is the ratio of Particles Moving speed and the light velocity defined in the Theory of Relativity; P is particle momentum.X ' is integration parameter; Θ
0=13.6MeV/c is a constant.Other parameters are calculating intermediate variable, without actual physics meaning.
Proton transport process in material is a complicated lasting scattering and energy attenuation process, and has the nuclear physics reaction of small probability to occur.The target atom of proton and penetrating material interacts and energy is decayed gradually, square being inversely proportional to of rate of energy loss and speed.Angular distribution probability all Gaussian distributed of test card phaneroplasm scattering each time, therefore for the proton of known incident and Exit positions, its probability of scattering in object under test inside can calculate.And the scattering of proton in each depth layer has a maximum probability point, the track connecting each degree of depth maximum of probability point is exactly maximum of probability track.According to above principle, utilize the data collected for twice, use monte carlo method analog computation can obtain maximum of probability track.And resolution and the precision of image effectively can be improved by the maximum of probability track imaging calculated.
The imaging system of the use real-time dynamic proton imaging designed by the present invention and the image imaging method of radiotherapy, it mainly comprises TFT Amorphous silicon flat-panel detectors and the analog-digital converter of band cesium iodide coating, this system has spaced two pieces of TFT Amorphous silicon flat-panel detectors, detection thing placement space is set between two pieces of TFT Amorphous silicon flat-panel detectors, described analog-digital converter is connected with signal synchronous processing device, and signal synchronous processing device is connected with picture system again.
The imaging system of the image imaging method of the real-time dynamic proton imaging that the present invention obtains and radiotherapy and image reconstruction optimization method and use the method, that innovates devises twice acquisition testing method, utilize the data residual quantity gathered, correction image, can be applicable to various field.But also devise and inject testee from different perspectives, gather the method for image, the proton images (proton energy scattergram) of different angles, by image layer analysis system, can produce the image information of object under test different depth layer.In proton radiation therapy field, from patient's body surface to the image of tumor depth and range information to guarantee with to improve therapeutic quality particularly important.In sum, the imaging system of the real-time dynamic proton imaging designed by the present invention and the image imaging method of radiotherapy and image reconstruction optimization method and use the method, there is output image resolution high, definition is high, image exports in real time, in field of radiation therapy potential applicability in clinical practice feature widely.
Accompanying drawing explanation
Fig. 1 is the diagram that embodiment 1 comprises the embodiment of the proton radiation detector assembly of the TFT Amorphous silicon flat-panel detectors of preposition and rearmounted band cesium iodide coating, signal synchronous processing device, analog-digital converter and picture system.
Fig. 2 is the computed tomography scanning application example of embodiment 2 proton imaging system.
Fig. 3 is proton maximum of probability track method schematic diagram.
In figure: 1. proton radiation source, 2. cesium iodide coating, 3.TFT Amorphous silicon flat-panel detectors, 4. analog-digital converter, 5. signal synchronous processing device, 6. picture system, 7. testee
Detailed description of the invention
Below by embodiment, the invention will be further described by reference to the accompanying drawings.
Embodiment 1:
As shown in Figure 1 and Figure 2, the imaging system of the image imaging method of the use real-time dynamic proton imaging that the present embodiment describes and radiotherapy, it mainly comprises TFT Amorphous silicon flat-panel detectors 3 and the analog-digital converter 4 of band cesium iodide coating 2, this system has spaced two pieces of TFT Amorphous silicon flat-panel detectors 3, arrange between two pieces of TFT Amorphous silicon flat-panel detectors 3 and be detected thing 7 placement space, described analog-digital converter 4 is connected with signal synchronous processing device 5, and signal synchronous processing device 5 is connected with picture system 6 again.
Real-time dynamic proton imaging described by the present embodiment and the image imaging method of radiotherapy, it is the energy being gathered proton radiation by the TFT Amorphous silicon flat-panel detectors of band cesium iodide coating, transformed by digital-to-analogue again and obtain image, it is characterized in that the energy acquisition carrying out twice proton radiation, it is before proton radiation is injected and is detected object that first time gathers, second time collection be at proton radiation after being detected object; Then the data collected for twice by process form image.
The image reconstruction optimization method being applicable to the image imaging method of real-time dynamic proton imaging and radiotherapy that the present embodiment describes, utilize the proton radiation that collects to inject energy attenuation rate before and after testee and positional information, reducing scattering of proton to the impact of imaging by calculating, finally obtaining two-dimentional water equivalence range (thickness) distributed image that proton beam injects Energy distribution two dimensional image before and after testee and object under test; It is calculated the maximum of probability track of proton by Monte Carlo simulation, and the maximum of probability track calculated the most at last is converted into two dimensional image that described calculating reduces the impact of scattering of proton on imaging.Calculation procedure and the formula of wherein said maximum of probability track are as follows:
The degree of depth (x
1) maximum of probability track point:
Wherein:
formula (1)
formula (2)
formula (3)
formula (4)
Wherein: assuming that incidence point coordinate is (0,0), x is eye point abscissa, and y is eye point vertical coordinate; θ is shooting angle, y
1with θ
1be respectively particle to be asked at degree of depth x
1the vertical coordinate at place and scattering angle; X
0for the radiant intensity constant relevant to passing material; β is the ratio of Particles Moving speed and the light velocity defined in the Theory of Relativity; P is particle momentum.X ' is integration parameter; Θ
0=13.6MeV/c is a constant.Other parameters are calculating intermediate variable, without actual physics meaning.(such as
respectively can by formula (1) above, (2) and (3) calculate, and then use formula (4) to calculate det (Δ again
1) and det (Δ
2))
As shown in Figure 3, proton is from (0,0) point along in incident 20 cm of water of x-axis positive direction, and wherein solid line is actual scattering track.And dotted line is the maximum of probability track by calculating.Actual scattering track and maximum of probability track error are 1.1 millimeters to the maximum as shown in Figure 3, it can thus be appreciated that the image resolution ratio obtained by this method also can reach 1.1 millimeters.
Embodiment 2
As shown in Figure 2, the image imaging of the use real-time dynamic proton imaging that the present embodiment describes and radiotherapy, moveable proton radiation source 1 and TFT Amorphous silicon flat-panel detectors 3 gather the image information of testee 7 from different perspectives.The proton images (proton energy scattergram) of different angles, by image layer analysis system, can produce the image information of object under test different depth layer.Be applied in proton radiation therapy field, the image from patient's body surface to tumor depth and range information are provided, ensure and improve therapeutic quality.
Claims (3)
1. the image imaging method of a real-time dynamic proton imaging and radiotherapy, it is the energy being gathered proton radiation by the TFT Amorphous silicon flat-panel detectors of band cesium iodide coating, transformed by digital-to-analogue again and obtain image, it is characterized in that the energy acquisition carrying out twice proton radiation, it is before proton radiation is injected and is detected object that first time gathers, second time collection be at proton radiation after being detected object; Then the data collected for twice by process form image; The proton radiation that collects is utilized to inject energy attenuation rate before and after testee and positional information, reducing scattering of proton to the impact of imaging by calculating, finally obtaining two-dimentional water equivalence range and thickness distribution image that proton beam injects Energy distribution two dimensional image before and after testee and object under test; It is calculated the maximum of probability track of proton by Monte Carlo simulation, and the maximum of probability track calculated the most at last is converted into two dimensional image that described calculating reduces the impact of scattering of proton on imaging.
2. the image imaging method of real-time dynamic proton imaging according to claim 1 and radiotherapy, it is characterized in that described proton radiation is penetrated by proton radiation source, proton radiation source is the proton radiation source of movement, proton radiation is injected from three different angles and is detected object, and the described position of TFT Amorphous silicon flat-panel detectors coordinates with the exposure pathways of proton radiation, and gather the proton radiation energy that different angles inject, finally obtain the image information of testee different depth layer.
3. one kind uses the imaging system of the image imaging method of real-time dynamic proton imaging as claimed in claim 1 and radiotherapy, it mainly comprises TFT Amorphous silicon flat-panel detectors and the analog-digital converter of band cesium iodide coating, this system has spaced two pieces of TFT Amorphous silicon flat-panel detectors, detection thing placement space is set between two pieces of TFT Amorphous silicon flat-panel detectors, described analog-digital converter is connected with signal synchronous processing device, and signal synchronous processing device is connected with picture system again.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201110328222.2A CN102488971B (en) | 2011-10-25 | 2011-10-25 | The image imaging method of real-time dynamic proton imaging and radiotherapy and the imaging system of use the method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201110328222.2A CN102488971B (en) | 2011-10-25 | 2011-10-25 | The image imaging method of real-time dynamic proton imaging and radiotherapy and the imaging system of use the method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102488971A CN102488971A (en) | 2012-06-13 |
| CN102488971B true CN102488971B (en) | 2016-02-17 |
Family
ID=46180859
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201110328222.2A Expired - Fee Related CN102488971B (en) | 2011-10-25 | 2011-10-25 | The image imaging method of real-time dynamic proton imaging and radiotherapy and the imaging system of use the method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102488971B (en) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2573178C1 (en) * | 2014-09-15 | 2016-01-20 | Российская Федерация, от имени которой выступает Государственная корпорация по атомной энергии "Росатом" - Госкорпорация "Росатом" | Method of imaging high-speed process using proton radiation |
| CN108478936B (en) * | 2018-02-28 | 2019-10-01 | 武汉大学 | The method and apparatus of proton therapeutic dosage and range is determined by proton-induced thermoacoustic signal |
| RU2687840C1 (en) * | 2018-08-17 | 2019-05-16 | Российская Федерация, от имени которой выступает Государственная корпорация по атомной энергии "Росатом" (Госкорпорация "Росатом") | Method of investigating behavior of materials during impact-wave loading using proton radiography |
| RU2690713C1 (en) * | 2018-09-03 | 2019-06-05 | Российская Федерация, от имени которой выступает Государственная корпорация по атомной энергии "Росатом" (Госкорпорация "Росатом") | Method of obtaining and processing images formed using proton radiation |
| CN112569478B (en) * | 2019-09-29 | 2023-05-16 | 王宏凯 | Method and device for synchronously reading digital proton CT imaging |
| RU2757865C1 (en) * | 2021-01-11 | 2021-10-21 | Российская Федерация, от имени которой выступает Государственная корпорация по атомной энергии "Росатом" (Госкорпорация "Росатом") | Method for obtaining images formed using proton radiation |
| RU2767592C1 (en) * | 2021-06-11 | 2022-03-17 | Российская Федерация, от имени которой выступает Государственная корпорация по атомной энергии "Росатом" (Госкорпорация "Росатом") | Method for obtaining and processing images formed using proton radiation |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101900826A (en) * | 2010-06-13 | 2010-12-01 | 中国科学院近代物理研究所 | Heavy ion beam current transverse dosage distribution measuring detector and two-dimensional imaging method thereof |
| CN102119586A (en) * | 2008-05-22 | 2011-07-06 | 弗拉迪米尔·叶戈罗维奇·巴拉金 | Multi-field charged particle cancer therapy method and apparatus |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070031337A1 (en) * | 2004-06-22 | 2007-02-08 | Reinhard Schulte | Nanoparticle enhanced proton computed tomography and proton therapy |
| US8278633B2 (en) * | 2009-01-20 | 2012-10-02 | Varian Medical Systems International Ag | Gated radiation procedure using packages |
-
2011
- 2011-10-25 CN CN201110328222.2A patent/CN102488971B/en not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102119586A (en) * | 2008-05-22 | 2011-07-06 | 弗拉迪米尔·叶戈罗维奇·巴拉金 | Multi-field charged particle cancer therapy method and apparatus |
| CN101900826A (en) * | 2010-06-13 | 2010-12-01 | 中国科学院近代物理研究所 | Heavy ion beam current transverse dosage distribution measuring detector and two-dimensional imaging method thereof |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102488971A (en) | 2012-06-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11819308B2 (en) | Systems and methods for real-time target validation for image-guided radiation therapy | |
| CN102488971B (en) | The image imaging method of real-time dynamic proton imaging and radiotherapy and the imaging system of use the method | |
| Richter et al. | First clinical application of a prompt gamma based in vivo proton range verification system | |
| Zhu et al. | Proton therapy verification with PET imaging | |
| CN104548372B (en) | The dosage determining device of radiotherapy | |
| Kurosawa et al. | Prompt gamma detection for range verification in proton therapy | |
| Testa et al. | Proton radiography and proton computed tomography based on time-resolved dose measurements | |
| US8515011B2 (en) | System and method for dose verification radiotherapy | |
| Rinaldi et al. | Experimental investigations on carbon ion scanning radiography using a range telescope | |
| McCurdy et al. | In vivo dosimetry for lung radiotherapy including SBRT | |
| KR101948800B1 (en) | 3d scattering radiation imager, radiation medical apparatus having the same and method for placing the 3d scattering radiation imager | |
| Cilla et al. | An in-vivo dosimetry procedure for Elekta step and shoot IMRT | |
| CN103729868A (en) | Dual-energy CT (Computer Tomography) scan data based detection method for reconstructing projected image | |
| Fuangrod et al. | A cine-EPID based method for jaw detection and quality assurance for tracking jaw in IMRT/VMAT treatments | |
| Bert et al. | Scanned carbon beam irradiation of moving films: comparison of measured and calculated response | |
| Han et al. | Proton radiography and fluoroscopy of lung tumors: a Monte Carlo study using patient‐specific 4DCT phantoms | |
| CN110392550B (en) | System and method for radiation beam range verification using sonic measurements | |
| Li et al. | Using Cherenkov imaging to monitor the match line between photon and electron radiation therapy fields on biological tissue phantoms | |
| CN107102349A (en) | Cerenkov's ray beam scanning measuring system | |
| Sellner et al. | Real-time imaging of antiprotons stopping in biological targets–Novel uses of solid state detectors | |
| US11975219B2 (en) | Systems and methods for particle portal imaging | |
| Parodi | In vivo treatment verification | |
| US12029922B2 (en) | X-ray CT-or MRI-based quantitative correction of Cherenkov light emission in radiation dose imaging | |
| US20220257982A1 (en) | X-ray ct calibration for quantitative correction of cherenkov light emission in radiation dose imaging | |
| Ting et al. | The feasibility of an approximate irregular field dose distribution simulation program applied to a respiratory motion compensation system |
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: 20160217 Termination date: 20191025 |