CN110141772A - Radioactive particle source total quantity acquisition methods needed for tumour - Google Patents
Radioactive particle source total quantity acquisition methods needed for tumour Download PDFInfo
- Publication number
- CN110141772A CN110141772A CN201910262824.9A CN201910262824A CN110141772A CN 110141772 A CN110141772 A CN 110141772A CN 201910262824 A CN201910262824 A CN 201910262824A CN 110141772 A CN110141772 A CN 110141772A
- Authority
- CN
- China
- Prior art keywords
- target site
- particle source
- radioactive particle
- tumour
- knurl
- 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/103—Treatment planning systems
- A61N5/1031—Treatment planning systems using a specific method of dose optimization
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 total quantity acquisition methods needed for tumour, comprising steps of S101, a radioactive particle source to be modeled as to the sphere of radius Ucm in three dimensions；S102, it obtains tumor volume, jeopardize organ and tissue volume in knurl, obtaining jeopardizing organ and tissue volume in practical packing volume=tumor volume-knurl；Radioactive particle source total quantity=practical packing volume needed for S103, tumour/(4 π/3) * U3cm3.It is more accurate to calculate for mode of the invention tumour irregular for form, radioactive particle source total activity needed for calculated is more nearly actually required radioactive particle source total activity, the damage that normal tissue is made utmostly is reduced, calculating radioactive particle source total quantity needed for tumour can be used as the switching threshold in later period radioactive particle source cloth source.
Description
Technical field
The invention belongs to radioactive prospecting instrument technical fields, and in particular to radioactive particle source total quantity needed for tumour obtains
Take 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.Finally determine the mode and method of seeds implanted, common particle plantation treatment there are 3 kinds of modes: template
It is planted in plantation, B ultrasound and the lower plantation of CT guidance, art.Such as template plantation, it is to make implant needle template, hand by above- mentioned information
When art, patient is fixed on CT bed, and implant needle template is fixed on skin site of the patient close to tumour, then implant needle is pressed
It is punctured according to step is pre-designed, while implantation pin position is checked by CT scan in real time, then by being arranged on implant needle
Scale provide depth reference.When implant needle reaches designated depth, starts to inject particle, then pull out implant needle, and
Particle is re-injected after reaching designated depth, can pull out implant needle until all particles on the implant needle all have been injected into.
The quantity of above-mentioned particle and position directly affect the degree of injury of therapeutic effect and normal tissue, therefore particle
Quantity and position it is extremely important.
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.
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 sources needed for tumour
Total quantity acquisition methods.
The technical scheme adopted by the invention is as follows:
Radioactive particle source total quantity acquisition methods, include the following steps: needed for tumour
S101, the sphere that a radioactive particle source is modeled as to radius Ucm in three dimensions；
S102, it obtains tumor volume, jeopardize organ and tissue volume in knurl, practical packing volume is obtained by formula (1):
Jeopardize organ and tissue volume (1) in practical packing volume=tumor volume-knurl；
S103, 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) * U3cm3 (2)。
Based on the above technical solution, the U value 1.
Based on the above technical solution, jeopardize the acquisition side of organ and tissue volume in the tumor volume and knurl
Formula are as follows:
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.
Based on the above technical solution, the 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 target site region and target
Position is with exterior domain；
A33, so that target site is connected as entirety with exterior domain, make faulted scanning pattern exposure mask；
A34, only target site region is can be obtained into the corresponding image masks information of initial three-dimensional labeled data dot product
The data of image.
Based on the above technical solution, the differentiation target site region and target site are with the specific reality of exterior domain
Existing mode is as follows:
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.
Based on the above technical solution, the production method of the 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.
Based on the above technical solution, the content of mark includes coordinate information, and coordinate information is based on breaking where mark
Coordinate system in layer scanning figure generates, and is woven in the relative position on faulted scanning pattern, coordinate system three for marked tumor group
Cartesian coordinate system is tieed up, tumor tissues and its hetero-organization in every tension fault scanning figure are indicated using three Cartesian coordinates
Relative position.
Based on the above technical solution, the content of mark further includes identification information, and identification information is used for present bit
The tissue mark set is that tumor tissues or its hetero-organization, identification information match with coordinate information, and current location corresponds to tissue
Identification information is endowed after the coordinate information that current location corresponds to tissue.
Based on the above technical solution, jeopardize the acquisition side of organ and tissue volume in the tumor volume and knurl
Formula is as follows:
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 invention has the benefit that
It is more accurate to calculate for mode of the invention tumour irregular for form, calculated required radioactivity grain
Component total activity is more nearly actually required radioactive particle source total activity, utmostly reduces the damage that normal tissue is made,
Calculating radioactive particle source total quantity needed for tumour can be used as the switching threshold in later period radioactive particle source cloth source.
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 total quantity acquisition methods needed for the tumour of the present embodiment, include the following steps:
S101, 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；
S102, it obtains tumor volume, jeopardize organ and tissue volume in knurl, practical packing volume is obtained by formula (1):
Jeopardize organ and tissue volume (1) in practical packing volume=tumor volume-knurl；
S103, 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 target site region and target
Position is with exterior domain；
A33, so that target site is connected as entirety with exterior domain, make faulted scanning pattern exposure mask；
A34, only target site region is 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.
It is more accurate to calculate for mode of the invention tumour irregular for form, calculated required radioactivity grain
Component total activity is more nearly actually required radioactive particle source total activity, utmostly reduces the damage that normal tissue is made,
Calculating radioactive particle source total quantity needed for tumour can be used as the switching threshold in later period radioactive particle source cloth source.
Radioactive particle source cloth source and needle track paths planning method in knurl is also disclosed in the present embodiment, radiates needed for tumour
It further include radioactive particle source Bu Yuan in knurl on the basis of property particle source total quantity acquisition methods.
Radioactive particle source Bu Yuan includes the following steps: in knurl
Radioactive particle source is filled in S201, knurl；
Radioactive particle source is arranged in S202, knurl.
The specific implementation of above-mentioned S201 is as follows:
S2011, tumour is extracted.
Tumour is the tumour iso-surface patch model divided.
It is realized particular by a1-a4 as above.
Radioactive particle source is filled in S2012, knurl.
It is realized particular by the optional one of following two mode:
The first, random filling:
C1, 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 1cm (effective radial extent) is generated at random in virtual square box, an equivalent sphere is one and puts
Penetrating property particle source；
C2, generation the sphere of radius 1cm (effective radial extent) ray is uniformly discharged by central point around, detection is every
Whether ray intersects with the surface of tumour, if then proving that the sphere is located inside tumor, is effective sphere, retains the ball
Body, if otherwise proving, the sphere is not entirely located in inside tumor, is invalid sphere, deletes the sphere, until number of spheres reaches
To stopping when the total quantity of radioactive particle source needed for tumour.
By central point, left and right uniformly discharges 6 rays to sphere up and down forwards, backwards.
The first has the beneficial effect that
The radioactive particle source filled at random fully considers the optimum position in radioactive particle source, is placed using such position
Radioactive particle source therapeutic effect it is best, and will not normal tissue cause big damage.
Second, linear filling:
D1, 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 1cm (effective radial extent), a ball are equidistantly generated on the cylinder and the center line of the superposition section of tumour
Body is equivalent to a radioactive particle source；
D2, 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, equidistantly generates radius 1cm on the cylinder and the center line of the superposition section of tumour and (effectively puts
Penetrate range) sphere；
D3, 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, radius 1cm (effectively radiation model is equidistantly generated on the cylinder and the center line of the superposition section of tumour
Enclose) sphere；
D4, circulation execute d3, until the virtual centrum formed stops when not intersecting with tumour；
When d5, 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 d2-d5 until radiation needed for number of spheres approach tumour
Terminate when property particle source total quantity.
Above-mentioned I ≠ P, m generally take 5.
Above-mentioned virtual cylinder is cylindrical body or square body.
It has the beneficial effect that for second
The radioactive particle source linearly filled fully considers the execution of automation equipment, compares first way, such line
Property filling radioactive particle source keep the execution convenience of automation equipment and validity higher.
Jeopardize the radioactive particle source processing of organ and tissue in S2013, knurl.
It is realized particular by following manner:
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.
The specific implementation of above-mentioned S202 includes that radioactive particle source level of linearity is arranged, radioactive particle source is dissociated
Jeopardize the radioactive particle source processing of organ and tissue in (non-linear) arrangement and knurl.
The arrangement of radioactive particle source level of linearity:
E1, the square that tumour is divided into multiple X*1*1cm, X are the length and width or high level of tumour, take the center of square respectively
Line is as linear reference line；
E2, the radioactive particle source in square center line range of set value is sucked in the square, and by sucking side
The central point in the radioactive particle source in block is adapted in the linear reference line of the square, and each radioactive particle source only moves one
Secondary, setting value is generally 1cm；
E3, 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.
Free (non-linear) arrangement in radioactive particle source is primarily directed to the trip generated because jeopardizing organ and tissue in knurl
Radioactive particle source from (non-linear) arrangement 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 source of free (non-linear)
Terminate when being 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.
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.
Radioactive particle source cloth source and needle track paths planning method also wrap on the basis of above technical scheme in knurl
It includes needle track path planning and needle track generates, needle track path planning and needle track are generated and generated including the planning of parallel pin path and needle track
It is generated with free angle needle track path planning and needle track.
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.
Radioactive particle source cloth source and needle track paths planning method also wrap on the basis of above technical scheme in knurl
Include coordinate conversion:
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.0031175.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 θ+cln r+d θ2+eln2r+fθln r+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).
2 I125 particle source radiation activity result of table
As can be seen from Table 2, 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 (9)
1. radioactive particle source total quantity acquisition methods needed for tumour, characterized by the following steps:
S101, the sphere that a radioactive particle source is modeled as to radius Ucm in three dimensions；
S102, it obtains tumor volume, jeopardize organ and tissue volume in knurl, practical packing volume is obtained by formula (1):
Jeopardize organ and tissue volume (1) in practical packing volume=tumor volume-knurl；
S103, 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) * U3cm3 (2)。
2. radioactive particle source total quantity acquisition methods needed for tumour according to claim 1, it is characterised in that: the U
Value 1.
3. radioactive particle source total quantity acquisition methods needed for tumour according to claim 1, it is characterised in that: the tumor
The acquisition pattern for jeopardizing organ and tissue volume in body volume and knurl is as follows:
B1, manually delineate or computer delineate in the tumor's profiles and knurl of each faulted scanning pattern of target site jeopardize organ 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, tumor volume is obtained, mesh
It marks and jeopardizes organ and tissue area superposition in the knurl of each faulted scanning pattern in position, then multiplied by thickness, obtain danger in knurl
And organ and tissue area.
4. radioactive particle source total quantity acquisition methods needed for tumour according to claim 1, it is characterised in that: the tumor
Jeopardize the acquisition pattern of organ and tissue volume in body volume and knurl are as follows:
A1, several faulted scanning patterns for obtaining target site, pre-process faulted scanning pattern, faulted scanning pattern are divided
Target site region and target site be out with exterior domain, and weeds out target site with exterior domain, will be in target site region
Tumor tissues and its hetero-organization are labeled, and are marked out and jeopardized organ and tissue in knurl, are obtained multiple mark samples, will be marked
Note sample is stored, and data set is obtained；A2,3D convolutional neural networks model is established, all mark sample data information are defeated
Enter into 3D convolutional Neural model and be trained, exports trained 3D convolutional Neural deep learning model；
A3, the faulted scanning pattern of target site is pre-processed, faulted scanning pattern is marked off into target site region and target
Position weeds out target site with exterior domain with exterior domain；
A4, trained 3D convolutional Neural deep learning model will be input to by pretreated faulted scanning pattern data information
In, and the 3D cubic block data divided are exported, multiple 3D cubic block data divided are merged, obtains having divided
Tumour iso-surface patch model；
A5, it calculates in tumor volume and knurl and jeopardizes organ and tissue volume.
5. radioactive particle source total quantity acquisition methods needed for tumour according to claim 3, it is characterised in that: described pre-
The specific implementation of processing 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 target site region and target site
With exterior domain；
A33, so that target site is connected as entirety with exterior domain, make faulted scanning pattern exposure mask；
A34, only target site area image is can be obtained into the corresponding image masks information of initial three-dimensional labeled data dot product
Data.
6. radioactive particle source total quantity acquisition methods needed for tumour according to claim 5, it is characterised in that: the area
Partial objectives for area and target site are 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 and institute
There is the corresponding color value of regional organization other than target site to be collected arrangement, obtains corresponding color value area, target site regional organization
Between color value corresponding with regional organization other than target site section, in this, as distinguishing target site regional organization and target site
The standard of regional organization in addition.
7. radioactive particle source total quantity acquisition methods needed for tumour according to claim 5, it is characterised in that: described disconnected
The production method of layer scanning figure exposure mask is as follows:
A331, by the corrosion treatment and expansion process in Morphological scale-space, so that target site is connected in one as far as possible with exterior domain
It rises, and eliminates the specific color value part in target site region as far as possible, to complete the production of the exposure mask in target site region.
8. radioactive particle source total quantity acquisition methods needed for tumour according to claim 4, it is characterised in that: mark
Content includes: coordinate information, and coordinate information is generated based on the coordinate system where mark on faulted scanning pattern, and is used for marked tumor
Group is woven in the relative position on faulted scanning pattern, and coordinate system is three Cartesian coordinates, using three Cartesian coordinates come
Indicate the relative position of the tumor tissues and its hetero-organization in every tension fault scanning figure.
9. radioactive particle source total quantity acquisition methods needed for tumour according to claim 8, it is characterised in that: mark
Content further include: identification information, identification information are used to the tissue mark of current location be tumor tissues or its hetero-organization.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910262824.9A CN110141772A (en) | 2019-04-02 | 2019-04-02 | Radioactive particle source total quantity acquisition methods needed for tumour |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910262824.9A CN110141772A (en) | 2019-04-02 | 2019-04-02 | Radioactive particle source total quantity acquisition methods needed for tumour |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110141772A true CN110141772A (en) | 2019-08-20 |
Family
ID=67589406
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910262824.9A Pending CN110141772A (en) | 2019-04-02 | 2019-04-02 | Radioactive particle source total quantity acquisition methods needed for tumour |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110141772A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116740768A (en) * | 2023-08-11 | 2023-09-12 | 南京诺源医疗器械有限公司 | Navigation visualization method, system, equipment and storage medium based on nasoscope |
-
2019
- 2019-04-02 CN CN201910262824.9A patent/CN110141772A/en active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116740768A (en) * | 2023-08-11 | 2023-09-12 | 南京诺源医疗器械有限公司 | Navigation visualization method, system, equipment and storage medium based on nasoscope |
| CN116740768B (en) * | 2023-08-11 | 2023-10-20 | 南京诺源医疗器械有限公司 | Navigation visualization method, system, equipment and storage medium based on nasoscope |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Jen-San Tsai Ph et al. | Dosimetric verification of the dynamic intensity-modulated radiation therapy of 92 patients | |
| CN105658279B (en) | Radiotherapy apparatus, system and method | |
| US8660325B2 (en) | Efficient user interaction with polygonal meshes for medical image segmentation | |
| US6661870B2 (en) | Fluence adjustment for improving delivery to voxels without reoptimization | |
| US9076222B2 (en) | Use of collection of plans to develop new optimization objectives | |
| Moteabbed et al. | Validation of a deformable image registration technique for cone beam CT‐based dose verification | |
| 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 | |
| Fattori et al. | Automated fiducial localization in CT images based on surface processing and geometrical prior knowledge for radiotherapy applications | |
| Hsi et al. | In vivo verification of proton beam path by using post‐treatment PET/CT imaging | |
| 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 | |
| CN110141773A (en) | Radioactive particle source arrangement method in knurl | |
| 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 | |
| CN110141772A (en) | Radioactive particle source total quantity acquisition methods needed for tumour | |
| Ecclestone et al. | Experimental validation of the van Herk margin formula for lung radiation therapy | |
| CN110141775A (en) | Radioactive particle source arrangement method in a kind of knurl | |
| CN110141776A (en) | The linear fill method in radioactive particle source in knurl | |
| CN110227213A (en) | The random fill method in radioactive particle source in knurl | |
| CN110141774A (en) | Radioactive particle source level of linearity arrangement method in knurl | |
| CN105797284A (en) | Radioactive particle radiation dose verification method and system | |
| Su et al. | Dosimetric effect of the intestinal gas of online adaptive stereotactic body radiotherapy on target and critical organs without online electron density correction for pancreatic cancer | |
| Boxwala et al. | Retrospective reconstruction of three-dimensional radiotherapy treatment plans of the thorax from two dimensional planning data | |
| Liu et al. | A fast online replanning algorithm based on intensity field projection for adaptive radiotherapy |
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 |
|
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190820 |