CN110141774A - Radioactive particle source level of linearity arrangement method in knurl - Google Patents
Radioactive particle source level of linearity arrangement method in knurl Download PDFInfo
- Publication number
- CN110141774A CN110141774A CN201910262862.4A CN201910262862A CN110141774A CN 110141774 A CN110141774 A CN 110141774A CN 201910262862 A CN201910262862 A CN 201910262862A CN 110141774 A CN110141774 A CN 110141774A
- Authority
- CN
- China
- Prior art keywords
- particle source
- radioactive particle
- tumour
- knurl
- radioactive
- 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
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
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)
- Radiation-Therapy Devices (AREA)
Abstract
The present invention discloses radioactive particle source level of linearity arrangement method in knurl, comprising steps of S101, extraction have been filled with the tumour in radioactive particle source；S102, tumour is divided into multiple squares, takes the center line of square as linear reference line respectively；S103, the radioactive particle source in square center line range of set value is sucked in the square, and the central point in the radioactive particle source sucked in square is adapted in the linear reference line of the square, each radioactive particle source is only mobile primary；S104, mutual distance is carried out to the radioactive particle after movement detect, be overlapped intersection or when distance is less than setting value when multiple radioactive particle sources occur, the radioactive particle source moved after deletion.The convenience that the present invention fully considers that automation equipment executes reduces the degree of injury of normal tissue under the premise of guaranteeing preferable therapeutic effect, reduces the fussy degree of operation.
Description
Technical field
The invention belongs to radioactive prospecting instrument technical fields, and in particular to radioactive particle source linear water flat raft in knurl
Cloth method.
Background technique
Interstitial radiation seeds implanted (also referred to as closely) cure, is by mini-radioactive resource (particle) implantation tumour
It is interior or by tumor-infiltrated tissue, including the tissue that malignant tumour is spread along Lymphatic channel, issued by radioactive particle source
The gamma-rays for continuing low energy, makes tumor tissues by maximum radiation injury and destruction, and normal tissue is injury-free
Or only by slight damage, to reach therapeutic purposes.Permanent implanted method and impermanency implantation can be divided by the seeds implanted time
Method.Expert thinks, compares other oncotherapy technologies, and seeds implantation technology technology content itself is not high, difficult
It spends and little.It but due to being implanted directly into human body, and is radioactive source, so strictly to hold indication.
Have different specific methods to the seeds implanted treatment of various different tumours, due to particle be planted in three-dimensional space into
Row, every kind of radioactive particle physical characteristic is different, needs to formulate a kind of special Three dimensional treatment planning system to every kind of nucleic, into
The spatial distribution of row simulation particle plantation determines the dosage distribution of particle plantation number and target area and surrounding risk organs, guidance
Clinical particle plantation.
Form, position, size and the relationship with adjacent organs, blood vessel for first having to clear tumour, depict the area for the treatment of
Domain.Conventional mode is scanned to the tumor region of patient, can be suffered from by the image checks such as B ultrasound, CT, ECT, MRI
Then the tumor region image of person carries out manually according to image or computer target area is drawn, carried out according to the target area figure drawn
Particle layout.Quantity of its secondary determining implantation particle and position, while confirming needle track position, this depend on tumour size and
The activity intensity of radioactive source, the quantity of particle and position directly affect the degree of injury of therapeutic effect and normal tissue.Most
The mode and method of seeds implanted are determined afterwards, and common particle plantation treatment there are 3 kinds of modes: template plantation, B ultrasound and CT guidance are sowed
It plants, planted in art.Such as template plantation, it is to make implant needle template by above- mentioned information, when operation, patient is fixed on CT
On, and implant needle template is fixed on skin site of the patient close to tumour, then implant needle is worn according to step is pre-designed
Thorn, while checking that implantation pin position, then the scale by being arranged on implant needle provide depth reference by CT scan in real time.When
When implant needle reaches designated depth, starts to inject particle, then pull out implant needle, and infuse again after reaching designated depth
Enter particle, implant needle can be pulled out until all particles on the implant needle all have been injected into.
In the prior art, in particle layout stage, it is all made of TPS system, can be carried out delineating confirmation tumour according to image
Then target area carries out particle layout again.However traditional particle layout type does not fully consider that automation equipment executes just
Benefit, next, which not can guarantee, is preferably minimized the degree of injury of human normal tissue.
Summary of the invention
In order to solve the above problems existing in the present technology, it is an object of that present invention to provide radioactive particle source lines in knurl
Property horizontal arrangement method.
The technical scheme adopted by the invention is as follows:
Radioactive particle source level of linearity arrangement method, includes the following steps: in knurl
S101, extraction have been filled with the tumour in radioactive particle source；
S102, tumour is divided into multiple squares, takes the center line of square as linear reference line respectively；
S103, the radioactive particle source in square center line range of set value is sucked in the square, and will sucking
The central point in the radioactive particle source in square is adapted in the linear reference line of the square, and each radioactive particle source is only moved
Once；
S104, mutual distance detection is carried out to the radioactive particle after movement, when the multiple radioactive particles of generation
Source is overlapped intersection or apart from radioactive particle source when being less than setting value, moved after deletion.
Based on the above technical solution, setting value is 1cm in the S103.
Based on the above technical solution, setting value is 5mm in the S104.
Based on the above technical solution, the size of the square is X*1*1cm, and X is the length and width or high level of tumour.
Based on the above technical solution, the fill method in S101 radioactive particle source includes: radiation in knurl
The property random fill method of particle source or linear fill method.
Based on the above technical solution, the implementation of the random fill method in radioactive particle source is such as in the knurl
Under:
S1011a, tumour is extracted；
S1012a, one and the matched virtual square box of maximum length of tumour are generated, tumour is embedded in the void
In quasi- square box, the sphere of radius Vcm is generated at random in virtual square box, an equivalent sphere is a radioactive particle
Source；
S1013a, generation the sphere of radius Vcm ray is uniformly discharged by central point around, whether detect every ray
Intersect with the surface of tumour, if then proving that the sphere is located inside tumor, is effective sphere, retains the sphere, if otherwise demonstrate,proving
The bright sphere is not entirely located in inside tumor, is invalid sphere, deletes the sphere, until number of spheres reaches and puts needed for tumour
Stopping when penetrating property particle source total quantity.
Based on the above technical solution, the implementation of the linear fill method in radioactive particle source is such as in the knurl
Under:
S1011b, tumour is extracted；
S1012b, an origin is taken in the lower part of tumour, the Z axis to pass through the origin generates a virtual column as center line
Body, equidistantly generates the sphere of radius Vcm on the cylinder and the center line of the superposition section of tumour, and an equivalent sphere is one
Radioactive particle source；
S1013b, m ° is tilted to X-axis with center line and rotates a circle to form virtual centrum, choose uniform I bus and make
It is respectively formed virtual cylinder for center line, equidistantly generates radius Vcm's on the cylinder and the center line of the superposition section of tumour
Sphere；
S1014b, m ° is tilted to X-axis again and rotates a circle to form virtual centrum, choose in uniform P bus conduct
Heart line is respectively formed virtual cylinder, and the sphere of radius Vcm is equidistantly generated on the cylinder and the center line of the superposition section of tumour；
S1015b, circulation execute S1014b, until the virtual centrum formed stops when not intersecting with tumour；
When S1016b, stopping, the number of spheres of generation is counted, the number of spheres and radioactive particle source needed for tumour is total
Quantitative comparison adjusts every time if radioactive particle source total quantity needed for number of spheres and tumour differs by more than setting value to X-axis
The item number of inclined degree and each virtual centrum Up Highway UHW, recycling executes S1013b-S1016b until number of spheres approaches
Terminate when the total quantity of radioactive particle source needed for tumour.
Based on the above technical solution, the tumour is the tumour iso-surface patch model divided.
Based on the above technical solution, the value of the V is 1.
Based on the above technical solution, the virtual cylinder is cylindrical body or square body.
Based on the above technical solution, the fill method in the radioactive particle source further include: in radioactive particle
After source generates, the radioactive particle source processing for jeopardizing organ and tissue in knurl is carried out；
The radioactive particle source processing for jeopardizing organ and tissue in the knurl is accomplished in that
Organ will be jeopardized in the radioactive particle source generated in knurl and knurl and tissue does collision detection and distance exam,
Jeopardize in organ and tissue or certain radioactive particle source central point and knurl in knurl if certain radioactive particle source central point passes through
When jeopardizing the shortest distance of organ and tissue less than setting value, which is considered as invalid radioactive particle source, deletes
The penetrating property particle source；It is on the contrary then be effective radioactive particle source, retain the penetrating property particle source.
The invention has the benefit that
The present invention fully considers the convenience that automation equipment executes, and under the premise of guaranteeing preferable therapeutic effect, reduces
The degree of injury of normal tissue reduces the fussy degree of operation.
Detailed description of the invention
Fig. 1 is the present invention-embodiment method flow diagram.
Fig. 2 is the schematic diagram of the present invention-embodiment I125 dose point.
Specific embodiment
With reference to the accompanying drawing and specific embodiment the present invention is further elaborated.
Embodiment:
As shown in Figure 1, radioactive particle source level of linearity arrangement method in the knurl of the present embodiment, includes the following steps:
S101, extraction have been filled with the tumour in radioactive particle source；
S102, the square that tumour is divided into multiple X*1*1cm, X are the length and width or high level of tumour, are taken in square respectively
Heart line is as linear reference line；
S103, the radioactive particle source in square center line range of set value is sucked in the square, and will sucking
The central point in the radioactive particle source in square is adapted in the linear reference line of the square, and each radioactive particle source is only moved
Once, setting value is generally 1cm；
S104, mutual distance detection is carried out to the radioactive particle after movement, when the multiple radioactive particles of generation
When source is overlapped intersection or distance less than setting value, the radioactive particle source moved after deletion, setting value is generally 5mm.
The present invention fully considers the convenience that automation equipment executes, and under the premise of guaranteeing preferable therapeutic effect, reduces
The degree of injury of normal tissue reduces the fussy degree of operation.
The fill method in above-mentioned S101 radioactive particle source includes: the random fill method in radioactive particle source or line in knurl
Property fill method.
The implementation of the random fill method in radioactive particle source is as follows in knurl:
S1011a, tumour is extracted；
S1012a, one and the matched virtual square box of maximum length of tumour are generated, tumour is embedded in the void
In quasi- square box, the sphere of radius 1cm (effective radial extent) is generated at random in virtual square box, an equivalent sphere is one
A radioactive particle source；
S1013a, generation the sphere of radius 1cm (effective radial extent) ray is uniformly discharged by central point around, examine
It surveys whether every ray intersects with the surface of tumour, is effective sphere if then proving that the sphere is located inside tumor, retaining should
Sphere, if otherwise proving, the sphere is not entirely located in inside tumor, is invalid sphere, deletes the sphere, until number of spheres
Stopping when reaching radioactive particle source total quantity needed for tumour.
The implementation of the linear fill method in radioactive particle source is as follows in knurl:
S1011b, tumour is extracted, tumour is the tumour iso-surface patch model divided；
S1012b, an origin is taken in the lower part of tumour, the Z axis to pass through the origin generates a virtual column as center line
Body, equidistantly generates the sphere of radius 1cm (effective radial extent) on the cylinder and the center line of the superposition section of tumour, and one
Equivalent sphere is a radioactive particle source；
S1013b, m ° is tilted to X-axis with center line and rotates a circle to form virtual centrum, choose uniform I bus and make
It is respectively formed virtual cylinder for center line, equidistantly generating radius 1cm on the cylinder and the center line of the superposition section of tumour (has
Imitate radial extent) sphere；
S1014b, m ° is tilted to X-axis with center line and rotates a circle to form virtual centrum, choose uniform I bus and make
It is respectively formed virtual cylinder for center line, equidistantly generating radius 1cm on the cylinder and the center line of the superposition section of tumour (has
Imitate radial extent) sphere；
S1015b, circulation execute S1014b, until the virtual centrum formed stops when not intersecting with tumour；
When S1016b, stopping, the number of spheres of generation is counted, the number of spheres and radioactive particle source needed for tumour is total
Quantitative comparison adjusts every time if radioactive particle source total quantity needed for number of spheres and tumour differs by more than setting value to X-axis
The item number of inclined degree and each virtual centrum Up Highway UHW, recycling executes S1013b-S1016b until number of spheres approaches
Terminate when the total quantity of radioactive particle source needed for tumour.
Above-mentioned I ≠ P, m generally take 5.
Above-mentioned virtual cylinder is cylindrical body or square body.
Radioactive particle source total activity needed for traditional tumour is by Paris formula (being averaged path length according to knurl) or root
It is obtained according to knurl Mass Calculation.
A, Paris formula (being averaged path length according to knurl) calculates.
The dosage that should be given is calculated using Dimension-averageing method:
A=Da*K
Wherein, A is total activity, unit mCi；Da is that knurl be averaged path length, unit cm, Da=(long+wide+height)/3；K
For fitting parameter, the K value for I125 is 5.
Resulting specific calculating formula such as table 1 is deformed by above-mentioned Paris formula.
Table 1
| Knurl is averaged path length Da (cm) | Particle source total activity A (mCi) |
| Da≤2.4 | A=5 (Da+1) |
| 2.4 Da≤3.24 < | A=3.87 (Da+1)1.293 |
| Da > 3.24 | A=2.76 (Da+1)1.581 |
B, according to knurl Mass Calculation.
Radioactive particle source total activity (mCi) needed for knurl=expectation tissue resorption dosage (cGy) * tumor weight (g)/
182。
Note: the quality of knurl can be calculated by CT.
The calculation pair of radioactive particle source total quantity needed for Paris formula (being averaged path length according to knurl) calculates tumour
Exist for the irregular tumour of form and calculates the biggish defect of error, the operation for needing artificially to adjust.According to knurl quality
The calculation of radioactive particle source total quantity needed for calculating tumour is substantially that the quality of knurl is calculated by CT, but tumour
Density is uncertain and uneven, therefore there is also calculate the biggish defect of error.Above two calculation does not account for
Jeopardize organ and tissue in knurl, for example, in the case where knurl wrapping portion blood vessel, what above two calculation was calculated
Required radioactive particle source total activity includes radioactive particle source radiation activity needed for capped blood vessel, and actually blood vessel is just
It often organizes, does not need radioactive particle source, therefore greatly possible calculated required radioactive particle source total activity is greater than reality
Required radioactive particle source total activity, i.e. overtreatment, this calculating inaccuracy, thereby increases and it is possible to which normal tissue causes big damage.
And the present invention calculates radioactive particle source total quantity needed for tumour in the following way:
S1011, the ball that a radioactive particle source is modeled as to radius 1cm (effective radial extent) in three dimensions
Body, radioactive particle source itself are capsule bodies, and simulation is (the 80% dosage covering half of I125 particle of radioactive particle source
Diameter is 1cm) region of certain dose value covering, it is similar to sphere, is only modeled as the sphere of radius 1cm herein；
S1012, it obtains tumor volume, jeopardize organ and tissue volume in knurl, practical obturator is obtained by formula (1)
Product:
Jeopardize organ and tissue volume (1) in practical packing volume=tumor volume-knurl；
S1013, radioactive particle source total quantity needed for tumour is obtained by formula (2):
Radioactive particle source total quantity=practical packing volume needed for tumour/(4 π/3) cm3 (2)。
Above-mentioned tumor volume is obtained in the following manner:
A1, several faulted scanning patterns for obtaining target site, pre-process faulted scanning pattern, by faulted scanning pattern
Target site region and target site are marked off with exterior domain, and weeds out target site with exterior domain, by target site region
In tumor tissues and its hetero-organization be labeled, and mark out and jeopardize organ and tissue in knurl, obtain multiple mark samples,
Mark sample is stored, data set is obtained；
A2,3D convolutional neural networks model is established, by all mark sample data information inputs to 3D convolutional Neural model
In be trained, export trained 3D convolutional Neural deep learning model；
A3, the faulted scanning pattern of target site is pre-processed, by faulted scanning pattern mark off target site region and
Target site weeds out target site with exterior domain with exterior domain；
A4, trained 3D convolutional Neural deep learning mould will be input to by pretreated faulted scanning pattern data information
In type, and the 3D cubic block data divided are exported, multiple 3D cubic block data divided are merged, obtain dividing
Tumour iso-surface patch model；
A5, it calculates in tumor volume and knurl and jeopardizes organ and tissue volume.
Above-mentioned pretreated specific implementation is as follows:
The pixel value of a31, standardized images, and probability density distribution is done to pixel value；
A32, the boundary between different zones tissue is found according to the distribution of pixel value, distinguish tumor tissue sections and other
Tissue regions；
A33, so that other tissue regions is connected as entirety, make faulted scanning pattern exposure mask；
A34, only other tissue regions are can be obtained into the corresponding image masks information of initial three-dimensional labeled data dot product
The data of image.
It distinguishes target site region and target site is as follows with the specific implementation of exterior domain:
A321, the color value for reading different zones on identification faulted scanning pattern, by the corresponding color value of target site regional organization
Color value corresponding with regional organization other than all target sites is collected arrangement, obtains the corresponding color of target site regional organization
It is worth the corresponding color value section of regional organization other than section and target site, in this, as differentiation target site regional organization and target
The standard of regional organization other than position.
The production method of faulted scanning pattern exposure mask is as follows:
A331, by the corrosion treatment and expansion process in Morphological scale-space, connect target site as far as possible with exterior domain
Together, and as far as possible the specific color value part in target site region is eliminated, to complete the exposure mask in target site region
Production.
The content of mark includes coordinate information, and coordinate information is generated based on the coordinate system where mark on faulted scanning pattern,
And it is woven in the relative position on faulted scanning pattern for marked tumor group, coordinate system is three Cartesian coordinates, utilizes three-dimensional
Cartesian coordinate system indicates the relative positions of tumor tissues and its hetero-organization in every tension fault scanning figure.
The content of mark further includes identification information, identification information be used for be by the tissue mark of current location tumor tissues or
Its hetero-organization, identification information match with coordinate information, and the identification information that current location corresponds to tissue is endowed current location pair
After the coordinate information that should be organized.
3D convolutional neural networks model mainly using U-net divide network+Res-net residual error network integration by the way of structure
It builds.Wherein the shallow-layer network of model mainly saves the location information of pixel, and deep layer network mainly saves the classification letter of pixel
Breath can ensure that image pixel segmentation by shallow-layer information supplement to deep layer using U-net network and Res-net network
Location accuracy and classification accuracy；Meanwhile Res-net network structure changes the target of study, by the defeated of learning model
Become the difference that study is output and input out, it reduce the learning difficulties of network model, restrain faster conducive to model.
Before carrying out 3D convolutional neural networks model training, the three-dimensional data marked in sample data information is cut to picture
The 3D cubic block data of plain 256*256*256 (or 128*128*128,64*64*64) are for training.Because three-dimensional data includes
The information of slice between layers carries out the training of 3D convolutional neural networks using three-dimensional data, the accuracy of result can be improved.
In addition, in the present embodiment by taking 256*256*256 3D cube module unit as an example, 3D convolutional neural networks modelling
The overall number of plies up to 160 layers, specifically: convolution+BN+maxpool is carried out to 256*256*256 cubes of module unit, to extract spy
Characteristic pattern is obtained, centre is used in combination residual error network (Res-net) and adds residual block；It extracts feature and reduces characteristic pattern to 16*
16*16, channel increase to 64；It reuses deconvolution combination residual error method and characteristic pattern is amplified to 256*256*256, in deconvolution
During amplification characteristic figure, with the characteristic pattern of position before being merged using the cat method in U-net network, finally obtain point
The tumour iso-surface patch model cut out.
Above-mentioned tumor volume is obtained by another way:
B1, it manually delineates or computer is delineated in the tumor's profiles and knurl of each faulted scanning pattern of target site and jeopardizes device
Official and tissue contours calculate and jeopardize organ and tissue area in the tumor area and knurl of each faulted scanning pattern of target site；
B2, the tumor area of each faulted scanning pattern of target site is superimposed, then multiplied by thickness, obtains tumor volume,
Jeopardize organ and tissue area superposition in knurl each faulted scanning pattern of target site, then multiplied by thickness, obtains knurl
Inside jeopardize organ and tissue area.
The present invention calculate tumour needed for radioactive particle source total quantity calculation tumour irregular for form and
It is more accurate that speech calculates, and calculated required radioactive particle source total activity is more nearly actually required radioactive particle source and always lives
Degree utmostly reduces the damage made of normal tissue, calculates after radioactive particle source total quantity needed for tumour can be used as
The switching threshold in phase radioactive particle source cloth source.
The fill method in above-mentioned radioactive particle source further include: at the radioactive particle source for jeopardizing organ and tissue in knurl
Reason.
The radioactive particle source processing for jeopardizing organ and tissue in knurl is accomplished in that
Organ will be jeopardized in the radioactive particle source generated in knurl and knurl and tissue does collision detection and distance exam,
Jeopardize in organ and tissue or certain radioactive particle source central point and knurl in knurl if certain radioactive particle source central point passes through
When jeopardizing the shortest distance of organ and tissue less than setting value, which is considered as invalid radioactive particle source, deletes
The penetrating property particle source；It is on the contrary then be effective radioactive particle source, retain the penetrating property particle source.
Above-mentioned setting value is generally 1cm.
In knurl after the arrangement of radioactive particle source level of linearity, the present embodiment further comprises the steps of: radioactive particle source
Free (non-linear) arrangement, free (non-linear) arrangement in radioactive particle source is primarily directed to because jeopardizing organ and tissue in knurl
The radioactive particle source of free (non-linear) arrangement generated, is realized there are two types of mode.
The first:
F11, the point centered on the radioactive particle source of tumour lowermost end, the Z axis to pass through the central point are created as center line
Build the virtual cylinder of a 1*1*Zcm；
F12, virtual cylinder tilt s ° to X-axis with center line and rotate a circle, during rotation, if having detected 2
The radioactive particle source of a or multiple free (non-linear) is then sucked in the virtual column body in virtual column body, and will be inhaled
The central point for entering the intracorporal radioactive particle source of virtual column, which is adapted on the center line of the virtual cylinder, no longer to be moved, and records rotation
Gyration；
F13, s ° is tilted to X-axis again and is rotated a circle, it is virtual that the radioactive particle source of free (non-linear) is sucked this
In cylinder, and the central point for sucking the intracorporal radioactive particle source of virtual column is adapted on the center line of the virtual cylinder no longer
It is mobile, and record rotation angle；
F14, circulation execute f13, until virtual cylinder stops when not intersecting with tumour；
When f15, stopping, whether the radioactive particle source for counting free (non-linear) has all been moved through, if then terminating；
If otherwise adjusting every time to the inclined degree of X-axis, recycling executes f12-f15 up to the radioactive particle of free (non-linear)
Terminate when source is all moved through；
F16, mutual distance detection is carried out to the radioactive particle after movement, when the multiple radioactive particle sources of generation
When being overlapped intersection or distance less than setting value, the radioactive particle source moved after deletion, setting value is generally 5mm.
Above-mentioned virtual cylinder is cylindrical body or square body.
Second:
F21, with the knurl immediately below tumor center's point, any is origin outside, to pass through the Z axis of the origin as center line wound
Build a virtual cylinder；
F22, virtual cylinder tilt t ° to X-axis with center line and rotate a circle, during rotation, if having detected 2
The radioactive particle source of a or multiple free (non-linear) is then sucked in the virtual column body in virtual column body, and will be inhaled
The central point for entering the intracorporal radioactive particle source of virtual column, which is adapted on the center line of the virtual cylinder, no longer to be moved, and records rotation
Gyration；
F23, t ° is tilted to X-axis again and is rotated a circle, it is virtual that the radioactive particle source of free (non-linear) is sucked this
In cylinder, and the central point for sucking the intracorporal radioactive particle source of virtual column is adapted on the center line of the virtual cylinder no longer
It is mobile, and record rotation angle；
F24, circulation execute f23, until virtual cylinder stops when not intersecting with tumour；
When f25, stopping, whether the radioactive particle source for counting free (non-linear) has all been moved through, if then terminating；
If otherwise adjusting every time to the inclined degree of X-axis, recycling executes f22-f25 up to the radioactive particle of free (non-linear)
Terminate when source is all moved through；
F26, mutual distance detection is carried out to the radioactive particle after movement, when the multiple radioactive particle sources of generation
When being overlapped intersection or distance less than setting value, the radioactive particle source moved after deletion, setting value is generally 5mm.
Above-mentioned virtual cylinder is cylindrical body or square body.
After (non-linear) arrangement step is dissociated in radioactive particle source, the present embodiment, which further comprises the steps of: in knurl, to be jeopardized
The processing of the radioactive particle source of organ and tissue.
The radioactive particle source processing for jeopardizing organ and tissue in knurl is as follows:
By after arrangement radioactive particle source and knurl in jeopardize organ and tissue do collision detection and distance exam, if certain
Radioactive particle source central point passes through to jeopardize in knurl to be jeopardized in organ and tissue or certain radioactive particle source central point and knurl
When the shortest distance of organ and tissue is less than setting value, which is considered as invalid radioactive particle source, deletes this and penetrates
Property particle source；It is on the contrary then be effective radioactive particle source, retain the penetrating property particle source.
Needle track path planning and needle track generation method is also disclosed in the present embodiment, the radioactive particle source total quantity needed for tumour
Radioactive particle source fill method in calculation method, knurl, radioactive particle source level of linearity arrangement method and radiation in knurl
Property free (non-linear) arrangement of particle source on the basis of, further comprise the steps of: the planning of parallel pin path and needle track generate and freely
The planning of angle pin path and needle track generate.
It includes two kinds of implementations that the planning of parallel pin path and needle track, which generate:
The first:
S301a, the square that tumour is divided into multiple X*1*1cm, X are the length and width or high level of tumour；
S302a, each square both ends are extended outwardly；
If the Duan Wei that extends in the same direction of S303a, all squares encounters and jeopardizes organ and tissue outside knurl, with each side of tumour
The center line of block generates needle track as linear reference line；If extending in the same direction the sum encountered in section and jeopardize organ and tissue outside knurl
Tumour itself is then carried out rotation until the Duan Wei that extends in the same direction of all squares encounters and jeopardizes organ and group outside knurl by >=setting value
It knits or encounters the square sum for jeopardizing organ and tissue outside knurl less than setting value, with each square of tumour in the rotation angle
Center line as linear reference line generate needle track.
Second:
S301b, the square that tumour is divided into multiple X*1*1cm, X are the length and width or high level of tumour；
S302b, point centered on a point is chosen in three dimensions, line centered on the Z axis to pass through the central point；
S303b, tumour is rotated a circle centered on central point, if the angle of the surfaces of revolution and center line is α；
The angle continuation of S304b, the adjustment surfaces of revolution and center line are rotated a circle centered on central point, which passes through
The normal and center line of central point are always positioned at same plane；
S305b, circulation execute S304b, until the angle of the surfaces of revolution and center line is+180 ° of α；
S306b, adjustment radius of turn, execute S303b-S305b；
S307b, circulation execute S306b, until radius of turn terminates more than setting value；
S308b, in rotary course, take Q in each radius of turn identical rotation angles, and record across knurl
Jeopardize the square sum of organ and tissue outside, degree of taking passes through the minimum corresponding rotation of square sum for jeopardizing organ and tissue outside knurl
Gyration generates needle track using the center line of each square of tumour as linear reference line in the rotation angle.
The implementation that free angle needle track path planning and needle track generate is as follows:
With the virtual cylinder center line under the rotation angle that is recorded in free (non-linear) the arrangement step in radioactive particle source
Needle track is generated, needle track is extended outwardly, retains the needle track if extended segment does not pass through and jeopardizes organ and tissue outside knurl, otherwise deletes
Except the needle track, and re-start free (non-linear) arrangement in radioactive particle source.
Needle track path planning and needle track generation method further include coordinate conversion on the basis of above technical scheme:
According to the isocenter point of automation equipment, corresponding points are found in faulted scanning pattern, are found out in three dimensions every
The coordinate information of a needle track and Δ X, Δ Y value and the rotation angle value and inclination angle value of isocenter point are transmitted to certainly
Dynamicization equipment, making it, needle track is run according to schedule.
Radioactive particle source cloth source and needle track paths planning method also wrap on the basis of above technical scheme in knurl
It includes and dose point in knurl is calculated after the step of radioactive particle source cloth source in knurl, isodose curve is delineated and calculates simultaneously
Draw DVH figure (dose volume histogram).
The specific implementation calculated dose point in knurl is as follows, as shown in Figure 2:
Z1011, dose point in knurl is calculated using formula (3):
Wherein, SkFor air kerma intensity (air kerma strength), meaning is in vacuum apart from long-range
In on the position of source size itself as energy greater than δ photon in infinitely small volume element air caused by air kerma rateWith d squares of distance of product, i.e. formula (4):
Wherein, Gr(r, θ)=r-2, point source fitting.
Wherein,Line source fitting.
Wherein, Λ is dose rate constant, in expression liquid water at reference point caused by the source of unit air Kerma intensity
Dosage rate, i.e. formula (5):
Value: 0.986 or 0.88 or 0.965.
Radial dose function (amendment of inverse ratio square law) is such as formula (6):
gx(r)=a0+a1r+a2r2+a3r3+a4r4+a5r5 (6)
Wherein, a0=1.08；a1=0.9707；a2=0.02957；a3=0.1429；a4=-0.02256； a5=
0.003117.Fitting correlation coefficient is 0.9999, sum of square of deviations 2.78*10-4.The formula scope of application be r=0.1~
10.0cm。
Anisotropy function (influence that description angle, θ is distributed dosage) is such as formula (7):
F (r, θ)=a+b θ+clnr+d θ2+eln2r+fθlnr+gθ3+hln3r+iθln2r+jθ2ln r (7)
In formula, a=0.2635；B=1.4991；C=0.1632；D=-1.005；E=-0.0284；F=-0.3150； g
=0.2224；H=0.0202；I=0.00155；J=0.1322.
It should be noted that the θ in formula (7) is indicated with radian value.The above fitting correlation coefficient R2=0.9926.
Sk value: (0.36-1.27).
3 I125 particle source radiation activity result of table
Seen from table 3, some particle source relative errors 13.4% exceed the limit value of relative deviation ± 5%, if this batch of grain
Component is used for clinical treatment, will be that overdose shines.
The specific implementation for delineating isodose curve is as follows:
Z1021, tumour is divided into multiple 1*1*1mm3Square, calculate wherein radioactive particle source radiation activity be superimposed post dose
Value is averaged, and the average value of calculating is considered as the dose value in the square, square central point is the expression point of dose point；
Z1022, the center progress line of multiple squares for possessing same dose value is formed in same level and is closed
Annular, the isodose curve of composition；
Z1023, the isodose curve line of the dose values such as multiple the dosage body such as formed into three dimensions.
It calculates and draws the specific implementation of DVH figure (dose volume histogram): certain is calculated according to isodose curve
The volume of one magnitude covering accounts for the ratio of knurl total volume, thus draws DVH figure.
The present invention is not limited to above-mentioned optional embodiment, anyone can show that other are various under the inspiration of the present invention
The product of form, however, make any variation in its shape or structure, it is all to fall into the claims in the present invention confining spectrum
Technical solution, be within the scope of the present invention.
Claims (11)
1. radioactive particle source level of linearity arrangement method in knurl, characterized by the following steps:
S101, extraction have been filled with the tumour in radioactive particle source；
S102, tumour is divided into multiple squares, takes the center line of square as linear reference line respectively；
S103, the radioactive particle source in square center line range of set value is sucked in the square, and square will be sucked
The central point in interior radioactive particle source is adapted in the linear reference line of the square, and each radioactive particle source only moves one
It is secondary；
S104, mutual distance detection is carried out to the radioactive particle after movement, it is heavy when multiple radioactive particle sources occur
Close intersection or apart from radioactive particle source when being less than setting value, moved after deletion.
2. radioactive particle source level of linearity arrangement method in knurl according to claim 1, it is characterised in that: described
Setting value is 1cm in S103.
3. radioactive particle source level of linearity arrangement method in knurl according to claim 1, it is characterised in that: described
Setting value is 5mm in S104.
4. radioactive particle source level of linearity arrangement method in knurl according to claim 1, it is characterised in that: the side
The size of block is X*1*1cm, and X is the length and width or high level of tumour.
5. radioactive particle source level of linearity arrangement method in knurl according to claim 1, it is characterised in that: described
The fill method in S101 radioactive particle source includes: the random fill method in radioactive particle source or linear fill method in knurl.
6. radioactive particle source level of linearity arrangement method in knurl according to claim 5, it is characterised in that: the tumor
The implementation of the random fill method of body radioactivity particle source is as follows:
S1011a, tumour is extracted；
S1012a, one and the matched virtual square box of maximum length of tumour are generated, tumour is embedded in the virtual side
In box body, the sphere of radius Vcm is generated at random in virtual square box, an equivalent sphere is a radioactive particle source；
S1013a, generation the sphere of radius Vcm ray is uniformly discharged by central point around, detect every ray whether with it is swollen
The surface of tumor is intersected, if then proving that the sphere is located inside tumor, is effective sphere, is retained the sphere, should if otherwise proving
Sphere is not entirely located in inside tumor, is invalid sphere, deletes the sphere, until number of spheres reaches radioactivity needed for tumour
Stop when particle source total quantity.
7. radioactive particle source level of linearity arrangement method in knurl according to claim 5, it is characterised in that: the tumor
The implementation of the linear fill method of body radioactivity particle source is as follows:
S1011b, tumour is extracted；
S1012b, an origin being taken in the lower part of tumour, the Z axis to pass through the origin generates a virtual cylinder as center line,
The sphere of radius Vcm is equidistantly generated on the cylinder and the center line of the superposition section of tumour, an equivalent sphere is one and puts
Penetrating property particle source；
S1013b, m ° is tilted to X-axis with center line and rotates a circle to form virtual centrum, choose in uniform I bus conduct
Heart line is respectively formed virtual cylinder, and the sphere of radius Vcm is equidistantly generated on the cylinder and the center line of the superposition section of tumour；
S1014b, m ° is tilted to X-axis again and rotates a circle to form virtual centrum, choose line centered on uniform P bus
It is respectively formed virtual cylinder, the sphere of radius Vcm is equidistantly generated on the cylinder and the center line of the superposition section of tumour；
S1015b, circulation execute S1014b, until the virtual centrum formed stops when not intersecting with tumour；
When S1016b, stopping, the number of spheres of generation is counted, by radioactive particle source total quantity needed for the number of spheres and tumour
Comparison is adjusted and is tilted every time to X-axis if radioactive particle source total quantity needed for number of spheres and tumour differs by more than setting value
Degree and each virtual centrum Up Highway UHW item number, recycling executes S1013b-S1016b until number of spheres approaches tumour
Terminate when required radioactive particle source total quantity.
8. radioactive particle source level of linearity arrangement method in knurl according to claim 6 or 7, it is characterised in that: institute
Stating tumour is the tumour iso-surface patch model divided.
9. radioactive particle source level of linearity arrangement method in knurl according to claim 6 or 7, it is characterised in that: institute
The value for stating V is 1.
10. radioactive particle source level of linearity arrangement method in knurl according to claim 7, it is characterised in that: described
Virtual cylinder is cylindrical body or square body.
11. radioactive particle source level of linearity arrangement method in knurl according to claim 6 or 7, it is characterised in that: institute
State the fill method in radioactive particle source further include: after radioactive particle source generates, carry out jeopardizing organ and tissue in knurl
Radioactive particle source processing；
The radioactive particle source processing for jeopardizing organ and tissue in the knurl is accomplished in that
Organ will be jeopardized in the radioactive particle source generated in knurl and knurl and tissue does collision detection and distance exam, if certain
Radioactive particle source central point passes through to jeopardize in knurl to be jeopardized in organ and tissue or certain radioactive particle source central point and knurl
When the shortest distance of organ and tissue is less than setting value, which is considered as invalid radioactive particle source, deletes this and penetrates
Property particle source；It is on the contrary then be effective radioactive particle source, retain the penetrating property particle source.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910262862.4A CN110141774A (en) | 2019-04-02 | 2019-04-02 | Radioactive particle source level of linearity arrangement method in knurl |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910262862.4A CN110141774A (en) | 2019-04-02 | 2019-04-02 | Radioactive particle source level of linearity arrangement method in knurl |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110141774A true CN110141774A (en) | 2019-08-20 |
Family
ID=67589398
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910262862.4A Pending CN110141774A (en) | 2019-04-02 | 2019-04-02 | Radioactive particle source level of linearity arrangement method in knurl |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110141774A (en) |
-
2019
- 2019-04-02 CN CN201910262862.4A patent/CN110141774A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105658279B (en) | Radiotherapy apparatus, system and method | |
| Jen-San Tsai Ph et al. | Dosimetric verification of the dynamic intensity-modulated radiation therapy of 92 patients | |
| US6661870B2 (en) | Fluence adjustment for improving delivery to voxels without reoptimization | |
| Litzenberg et al. | Daily prostate targeting using implanted radiopaque markers | |
| Pickett et al. | Static field intensity modulation to treat a dominant intra-prostatic lesion to 90 Gy compared to seven field 3-dimensional radiotherapy | |
| US7672705B2 (en) | Weighted surface-to-surface mapping | |
| US8660325B2 (en) | Efficient user interaction with polygonal meshes for medical image segmentation | |
| Keall et al. | On the use of EPID‐based implanted marker tracking for 4D radiotherapy: EPID‐based marker tracking | |
| US8929621B2 (en) | Methods and systems for segmentation and surface matching | |
| McCarter et al. | Evaluation of the validity of a convolution method for incorporating tumour movement and set-up variations into the radiotherapy treatment planning system | |
| CN103933672B (en) | Radiation therapy planning device and therapy planning method | |
| CN108815721A (en) | A kind of exposure dose determines method and system | |
| CN110141777A (en) | Radioactive particle source cloth source and needle track paths planning method in knurl | |
| Aubin et al. | Robustness and precision of an automatic marker detection algorithm for online prostate daily targeting using a standard V‐EPID | |
| CN107980008A (en) | Medical Image Processing equipment, treatment system and Medical Image Processing program | |
| Hsi et al. | In vivo verification of proton beam path by using post‐treatment PET/CT imaging | |
| CN110141773A (en) | Radioactive particle source arrangement method in knurl | |
| Lu et al. | Dose‐surface histograms as treatment planning tool for prostate conformal therapy | |
| CN110141771A (en) | The planning of radioactive particle source implant needle path and needle track generation method in knurl | |
| CN110141770A (en) | The planning of radioactive particle source implant needle path and needle track generation method in a kind of knurl | |
| CN110141778A (en) | The acquisition of radioactive particle source total quantity, radioactive particle source cloth source and needle track paths planning method in knurl | |
| Neal et al. | Tangential breast irradiation—rationale and methods for improving dosimetry | |
| Chen et al. | Object‐constrained meshless deformable algorithm for high speed 3D nonrigid registration between CT and CBCT | |
| CN110141775A (en) | Radioactive particle source arrangement method in a kind of knurl | |
| Ecclestone et al. | Experimental validation of the van Herk margin formula for lung radiation therapy |
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 | ||
| CB02 | Change of applicant information |
Address after: 610000 No.804, floor 8, unit 1, building 6, Zone D, West Zhigu, No.166, Wuxing Fourth Road, Wuhou New Town Management Committee, Wuhou District, Chengdu City, Sichuan Province Applicant after: Chengdu real dimension Technology Co., Ltd. Applicant after: Affiliated Zhongshan Hospital of Dalian University Address before: 610000 No. 801, Building 8, Building 6, Unit 1, West Zhigu D District, 166 Wuxing Road, Wuhou New Town Management Committee, Chengdu, Sichuan Province Applicant before: Chengdu real dimension Technology Co., Ltd. Applicant before: Affiliated Zhongshan Hospital of Dalian University |
|
| CB02 | Change of applicant information | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190820 |
|
| RJ01 | Rejection of invention patent application after publication |