CN108074636A - A kind of surface incident dose computational methods, equipment and storage medium - Google Patents
A kind of surface incident dose computational methods, equipment and storage medium Download PDFInfo
- Publication number
- CN108074636A CN108074636A CN201711436101.3A CN201711436101A CN108074636A CN 108074636 A CN108074636 A CN 108074636A CN 201711436101 A CN201711436101 A CN 201711436101A CN 108074636 A CN108074636 A CN 108074636A
- Authority
- CN
- China
- Prior art keywords
- value
- incident dose
- apart
- kilovolt
- bulb
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/17—Function evaluation by approximation methods, e.g. inter- or extrapolation, smoothing, least mean square method
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Measurement Of Radiation (AREA)
Abstract
The invention discloses a kind of surface incident dose computational methods, include the following steps:The kilovolt value of high pressure generator and milliampere seconds value are set, changes kilovolt value, obtains the incident dose apart from bulb focus different position;According to N number of incident dose apart from bulb focus different position, first matched curve of the N items apart from bulb focal position and incident dose is obtained;Wherein, N is positive integer;Assuming that the milliampere seconds value of high pressure generator is fixed, fitting coefficient during arbitrary kilovolt value is determined, obtain the second matched curve apart from bulb focal position and incident dose;According to the milliampere seconds value of high pressure generator and the relation of incident dose, the 3rd matched curve when obtaining arbitrary kilovolt value, arbitrary milliampere seconds value apart from bulb focal position and incident dose;When carrying out image check, patient is inputted into the 3rd matched curve apart from the position of bulb focus, obtains the incident dose of X-ray.This method can prevent incident dose overshoot risk.
Description
Technical field
The present invention relates to a kind of surface incident dose computational methods, equipment and storage mediums, belong to medical measurement technology neck
Domain.
Background technology
X ray image equipment is the present hospitals master equipment to be used that medical imaging is provided for patient, main
Principle is X-ray generates different X-ray attenuations through different institutional frameworks after human body, and flat panel detector is according to by human body
The X-ray generated afterwards be distributed to be formed with the corresponding intensity profile image of human tissue structure, pass through this image so as to doctor
Complete the observation and diagnosis to human tissue structure.
In the clinical practice of X ray image equipment, the surface incident dose control of patient is that each side all paid close attention in recent years
The problem of, because excessive X-ray radiation can increase the risk of human body radiation damage.EU countries has all formulated statutory standard will
It asks and establishes patient radiation dose's filing system, suffered X-ray radiation dosage is remembered when a patient per is received X-ray examination
It records in patient files, medical institutions can also carry out about that X is penetrated several times according to personal intergal dose record value to assess in the current year
Ray examination or assess whether this patient has received excessive X-ray examination.This provides the European Union for having been written to new edition
In statutory standard.It is every export to Europe various X ray image equipment must have patient surface incident dose acquisition and
Writing function.So doctor, when carrying out X ray image inspection to patient, used equipment is measurable or calculating patient enters
Penetrating dosage just becomes a function that must have.
Existing X ray image equipment, as angiography machine, Direct Digital flat panel X-ray imaging system or digital GI are made
Shadow machine etc. mainly measures the dose value of X-ray by being mounted on the dose sensor before X-ray tube.Measure obtained dosage
Value can be converted to a voltage signal, be transmitted on process circuit plate and be AD converted, AD conversion value is sent to image workstation
In computer.While radioscopic image is gathered, the dose value measured at this time is also stored in image format file and case
In file.In this way, the surface incident dose that patient when shooting this image receives, disease can be seen when checking this image by doctor
Also this value can be recorded in example text part, facilitates medical institutions' statistics.
But sensor in the path of X-ray must be installed using sensor measurement patient surface's incident dose and needed
Circuit design is AD converted, can so dramatically increase the complexity and cost of product design.
The content of the invention
In view of the deficiencies of the prior art, primary technical problem to be solved by this invention is to provide a kind of surface incidence agent
Measure computational methods.
Another technical problem to be solved by this invention provides a kind of surface incident dose computing device.
Another technical problem to be solved by this invention provides a kind of storage medium.
For achieving the above object, the present invention uses following technical solutions:
It is according to embodiments of the present invention in a first aspect, provide a kind of surface incident dose computational methods, include the following steps:
The kilovolt value of high pressure generator and milliampere seconds value are set, changes kilovolt value, obtains apart from bulb focus different position
Incident dose;
According to N number of incident dose apart from bulb focus different position, N items are obtained apart from bulb focal position and incident agent
First matched curve of amount;Wherein, N is positive integer;
Assuming that the milliampere seconds value of high pressure generator is fixed, determine fitting coefficient during arbitrary kilovolt value, obtain apart from bulb
Focal position and the second matched curve of incident dose;
According to the milliampere seconds value of high pressure generator and the relation of incident dose, arbitrary kilovolt value, arbitrary milliampere seconds value are obtained
When the 3rd matched curve apart from bulb focal position and incident dose;
When carrying out image check, patient is inputted into the 3rd matched curve apart from the position of bulb focus, obtains X-ray
Incident dose.
Wherein more preferably, the kilovolt value of high pressure generator and milliampere seconds value are set, change kilovolt value, is obtained burnt apart from bulb
The incident dose of point different position, including following sub-step:
S111 obtains the kilovolt maximum of the kilovolt value of high pressure generator and kilovolt minimum value, in kilovolt maximum and thousand
N number of kilovolt of value is chosen between volt minimum value, the N that is set to kVj, j=1,2 ...;
S112 sets value R1 on the basis of the milliampere seconds value of high pressure generator;
S113, according to patient support to the distance of bulb focus, determine it is N number of apart from Along ent, apart from bulb focus
Position is respectively that FIDi, i=1,2 ... N, N are positive integer;
S114 sets the kilovolt value kVj of high pressure generator, obtains the incidence apart from bulb focus different position FIDi one by one
Dosage;
S115 repeats step S114, and j is made to distinguish value 1,2 ... N, is obtained under each kilovolt value, apart from bulb focus not
With the incident dose of position.
Wherein more preferably, according to N number of incident dose apart from bulb focus different position, N items are obtained apart from bulb focus
Position and the first matched curve of incident dose, including following sub-step:
S121, using, as abscissa, the incident dose value of each point is ordinate, establishes coordinate apart from the position of bulb focus
System;
S122 when the kilovolt value of high pressure generator is kVj, will be represented apart from the position of bulb focus and incident dose one by one
Point distribution in a coordinate system;
S123, the point in coordinate system, carries out curve fitting, and obtains one apart from bulb focal position and incident dose
The first matched curve;
S124 repeats step S122~S123, and j difference value 1,2 ... N is made to carry out curve fitting, obtains N item distances
Bulb focal position and the first matched curve of incident dose.
Wherein more preferably, it is assumed that the milliampere seconds value of high pressure generator is fixed, and determines fitting coefficient during arbitrary kilovolt value, bag
Include following sub-step:
Using kilovolt value as abscissa, fitting coefficient is ordinate, establishes coordinate system;
The point distribution of kilovolt value and fitting coefficient will be represented in a coordinate system one by one;
Point in coordinate system, carries out curve fitting, and obtains the functional relation of fitting coefficient and kilovolt value, passes through kilovolt
It is worth to fitting coefficient.
Wherein more preferably, the surface incident dose computational methods before the incident dose of X-ray is calculated, further include
Following sub-step:
According to the milliampere seconds value of application mesohigh generator and the relationship error of incident dose, the 3rd matched curve is carried out
It corrects, obtains final matched curve.
Wherein more preferably, according to the milliampere seconds value of application mesohigh generator and the relationship error of incident dose, to the 3rd
Matched curve is corrected, and obtains final matched curve, including following sub-step:
In fixed FID points position, the kilovolt value of high pressure generator is arranged to kV1, and milliampere seconds value is arranged to R2;Make
The incident dose RB of this position is measured with dosemeter;
In the fixed FID points position, the kilovolt value of high pressure generator is arranged to kV1, and milliampere seconds value is arranged to R1;It uses
Dosemeter measures the incident dose RA of this position;
It calculates in FID points position, milliampere seconds value is arranged to incident dose RBS values during R2;In RBS=(R2/R1) × RA;
According to the incident dose that dosemeter measures and the incident dose for calculating acquisition, parameter dosage multiplying power factor is obtained, it will
The parameter dosage multiplying power factor brings the 3rd matched curve into, obtains final matched curve.
Wherein more preferably, apart from bulb focal position and incident dose when the arbitrary kilovolt value, arbitrary milliampere seconds value
Final matched curve is represented using equation below:
Y=(R/R1) × M × f (kV) × x-u;
Wherein, y is incident dose, and R is the milliampere seconds value under conditions present, and R1 is a reference value of milliampere seconds value, and f (kV) is
Fitting coefficient;X is current location apart from the position of bulb focus;U is function coefficients, and M is parameter dosage multiplying power system.
Second aspect according to embodiments of the present invention provides a kind of surface incident dose computing device, including processor and
Memory;The available computer program run on the processor is stored on the memory, when the computer program
Following steps are realized when being performed by the processor:
The kilovolt value of high pressure generator and milliampere seconds value are set, changes kilovolt value, obtains apart from bulb focus different position
Incident dose;
According to N number of incident dose apart from bulb focus different position, N items are obtained apart from bulb focal position and incident agent
First matched curve of amount;Wherein, N is positive integer;
Assuming that the milliampere seconds value of high pressure generator is fixed, determine fitting coefficient during arbitrary kilovolt value, obtain apart from bulb
Focal position and the second matched curve of incident dose;
According to the milliampere seconds value of high pressure generator and the relation of incident dose, arbitrary kilovolt value, arbitrary milliampere seconds value are obtained
When the 3rd matched curve apart from bulb focal position and incident dose;
When carrying out image check, patient is inputted into the 3rd matched curve apart from the position of bulb focus, obtains X-ray
Incident dose.
Wherein more preferably, when the computer program is performed by the processor, following steps are also realized;
According to the milliampere seconds value of application mesohigh generator and the relationship error of incident dose, the 3rd matched curve is carried out
It corrects, obtains final matched curve.
The third aspect according to embodiments of the present invention provides a kind of computer readable storage medium, described computer-readable
Storage medium is stored with one or more program, and one or more of programs can be held by one or more processor
Row, to realize the step in above-mentioned surface incident dose computational methods.
Incident dose computational methods in surface provided by the present invention, during by obtaining arbitrary kilovolt value, arbitrary milliampere seconds value
The 3rd matched curve apart from bulb focal position and incident dose;When carrying out image check, by patient apart from bulb focus
Position input the 3rd matched curve, obtain the incident dose of X-ray.So that in the case of without using dose sensor, it can
The surface incidence agent of patient's receiving is calculated apart from bulb focal position according to kV setting values and mAs setting values and patient surface
Amount, escapable cost, and can patient surface's acceptable dose be prompted by calculating at any time, it can prevent excess radiation risk.
Description of the drawings
Fig. 1 is the flow chart of incident dose computational methods in surface provided by the present invention;
Fig. 2 (a) is in embodiment provided by the present invention, and the structure of the first matched curve under 31kV, 20mAs setting is shown
It is intended to;
Fig. 2 (b) is in embodiment provided by the present invention, 35kV set under, mAs values when taking 20,50,100 respectively the
The structure diagram of one matched curve;
Fig. 2 (c) is in embodiment provided by the present invention, and the structure of the first matched curve under 38kV, 20mAs setting is shown
It is intended to;
Fig. 3 is the structure diagram of the surface incident dose computing device provided in the embodiment of the present invention.
Specific embodiment
Detailed specific description is carried out to the technology contents of the present invention in the following with reference to the drawings and specific embodiments.
As shown in Figure 1, incident dose computational methods in surface provided by the present invention, include the following steps:First, set high
The kV values of generator and mAs values are pressed, changes kV values, obtains N number of incident dose apart from bulb focus different position;Secondly, root
According to N number of incident dose apart from bulb focus different position, N items are obtained apart from bulb focal position and the first plan of incident dose
Close curve;Then, it is assumed that the mAs values of high pressure generator are fixed, and determine fitting coefficient during arbitrary kV values, are obtained burnt apart from bulb
Point position and the second matched curve of incident dose.Finally, according to the mAs values of high pressure generator and the relation of incident dose, obtain
To the 3rd matched curve apart from the position of bulb focus and incident dose when arbitrary kV values, arbitrary mAs values;When progress image inspection
When looking into, patient is inputted into the 3rd matched curve apart from the position of bulb focus (i.e. from bulb focal position), obtains entering for X-ray
Penetrate dosage.This processing procedure is described in detail below.
It should be noted that the kV values and mAs values of high pressure generator are to characterize the two indices of high pressure generator performance,
In, kV values are also known as kilovolt value, are the voltage values that high pressure generator bulb both ends add.MAs values are also known as milliampere seconds value, are sent out in X-ray
The total amount of X-ray is usually described inside raw device with milliampere seconds value.
S11 sets the kV values of high pressure generator and mAs values, changes kV values, obtains entering apart from bulb focus different position
Penetrate dosage.
Above-mentioned steps specifically include following sub-step:
S111 obtains the kV maximums of the kV values of high pressure generator and kV minimum values, kV maximums and kV minimum values it
Between choose N number of kV values, the N that is set to kVj, j=1,2 ..., N is positive integer.
All there are one adjustable ranges for the kV values of medical high pressure generator, it is first determined the kV of currently used high pressure generator
It is worth adjustable range, that is, determines kV maximums and kV minimum values.In embodiment of the method provided by the present invention, by what is be adjusted out
KV minimum values are defined as kV1.The value of kV1 is recorded.
Choose N number of kV values point.According to the kV maximums of kV values and kV minimum values, several kV points are chosen, selection rule is general
It is to choose several kV points according to certain interval.In embodiment of the method provided by the present invention, kV2, kV3 are defined as ...
kVN.These kV values that record has been chosen well.
S112 sets value R1 on the basis of the mAs values of high pressure generator.
The mAs values of medical high pressure generator choose a clinically used mAs value and define this there is also an adjustable range
It is worth for R1.Record the mAs values R1 chosen.
S113 according to patient support to the distance of bulb focus, determines N number of apart from Along ent, is apart from bulb focus
Position, the N that is respectively set to FIDi, i=1,2 ..., N is positive integer.
It according to actual patient supporting surface to the distance of bulb focus, determines N number of apart from Along ent, distinguishes away from focal length
The FIND that is defined as FID1, FID2, FID3 ..., unit use mm.
S114 sets the kV value kVj of high pressure generator, obtains the incident agent apart from bulb focus different position FIDi one by one
Amount.
The kV values of high pressure generator are set, in the case where kV1 and R1mAs is set, dosemeter is individually positioned in each FIDi positions
Point measures the incident dose value of FID1, FID2, FID3, FIDN each point one by one, and the measurement incident dose value of each position point is recorded
Into form.
S115 repeats step S114, and j is made to distinguish value 1,2 ... N, is obtained under each kV values, different apart from bulb focus
The incident dose of position.
Step S114 is repeated, measures the incident dose of each FIDi location points under kV2 and R1mAs.According to the kV values of selection
Point, when measuring kV2, the incident dose value of FID1, FID2, FID3, FIDN each point remembers the measurement incident dose value of each position point
It records in form.
The incident dose value that step carries out kV3 each FIDi points when the R1mAs of the kVN each kV points chosen is set is repeated,
The incident dose value of each FID point of the retest step by each kV points in R1mAs tests to obtain, and form is recorded
In.
S12 according to N number of incident dose apart from bulb focus different position, obtains N items apart from bulb focal position with entering
Penetrate the first matched curve of dosage.
Above-mentioned steps specifically include following sub-step::
S121, using, as abscissa, the incident dose value of each point is ordinate, establishes coordinate apart from the position of bulb focus
System;
Data derived above are the incident dose values of each FID points under each kV points R1mAs, are horizontal stroke with FIDi points distance
Coordinate, the incident dose value of each point is ordinate, establishes coordinate system.
S122 when the kV values of high pressure generator are kVj, will be represented apart from the position of bulb focus and incident dose one by one
Point distribution is in a coordinate system.
When high pressure generator kV values be kV1 when, can obtain in step sl it is N number of apart from the position of bulb focus with entering
The respective value of dosage is penetrated, according to coordinate system, the point represented apart from the position of bulb focus and incident dose is distributed in seat one by one
In mark system.
S123, the point in coordinate system, carries out curve fitting, and obtains one apart from bulb focal position and incident dose
The first matched curve.
Data derived above are the incident dose values of each FID points under each kV points R1mAs, for each kV value,
It can be fitted and generate function y=A × x-u.Wherein, the method for obtaining the first matched curve can be any existing by existing
The curve-fitting method apart from bulb focal position and incident dose incidence relation is found a little.It is not just repeating herein.
S124 repeats step S122~S123, and j difference value 1,2 ... N is made to carry out curve fitting, obtains N item distances
Bulb focal position and the first matched curve of incident dose.
Step S122~S123 is repeated, j difference value 1,2 ... N is made to carry out curve fitting, it is burnt apart from bulb to obtain N items
Point position and the first matched curve of incident dose, will obtain N number of A values and u values.
S13, the mAs values of high pressure generator fix, determine fitting coefficient during arbitrary kV values, obtain apart from bulb focus position
Put the second matched curve with incident dose.
Preceding to have addressed, first matched curve of the N items apart from bulb focal position and incident dose will obtain N number of A values and u
Value.
The AN that N number of A values is defined as A1, A2 ....A1 is corresponding with kV1, and A2 is corresponding with kV2, and AN and kVN is pair
It answers.This fitting, which obtains x in formula, can represent arbitrary FID values, and u is the coefficient that fitting obtains.
On fitting coefficient u, the basic phase of function coefficients u values fitted through overtesting discovery in each kV1 ... kVN gears
Together, in embodiment of the method provided by the present invention, only provide part kV values and mAs values setting under using FIDi point distances as
Abscissa, the incident dose value of each point is the data point of ordinate and the first matched curve being fitted to.Such as table 1, be 31kV,
Under 20mAs is set, using FIDi point distances as abscissa, the incident dose value of each point is the data point of ordinate and is fitted to
First matched curve;Shown in its corresponding matched curve such as Fig. 2 (a).
Under table 1 31kV, 20mAs are set, FIDi points and incident dose value contrast relationship table
Such as table 2, under being set for 35kV, mAs values are when taking 20,50,100 respectively, using FIDi point distances as abscissa, each point
Incident dose value is the data point of ordinate and the first matched curve being fitted to;Its corresponding matched curve such as Fig. 2 (b) institutes
Show.Wherein, the line that diamond shape is worn is represented under 35kV, 20mAs setting, using FIDi point distances as abscissa, the incident dose of each point
It is worth the first matched curve being fitted to for ordinate;The line that square is worn represent 35kV, 50mAs set under, with FIDi points away from
The first matched curve being fitted to from the incident dose value for abscissa, each point for ordinate;The line that triangle is worn represents
Under 35kV, 100mAs are set, using FIDi point distances as abscissa, the incident dose value of each point is intended for ordinate is fitted to first
Close curve.
Under 2 35kV of table is set, during different mA s values, FIDi points and incident dose value contrast relationship table
Such as table 3, under being set for 38kV, 20mAs, using FIDi point distances as abscissa, the incident dose value of each point is sat to be vertical
Target data point and the first matched curve being fitted to;Shown in its corresponding matched curve such as Fig. 2 (c).
Under table 3 38kV, 20mAs are set, FIDi points and incident dose value contrast relationship table
According to tracing analysis, the value of u is essentially identical.In embodiment of the method provided by the present invention, each u values of N are taken
Average for formula calculate.
So at this time the R1mAs under these selected kV shelves of this formula calculating kV1, kV2, kV3, kVN can be utilized to set
When putting, the incident dose value apart from bulb focus FID values position.But actual not just these kV gears of high pressure generator,
So it needs to be determined that a fitting coefficient A is determined under arbitrary kV values, when mAs values are R1, the incident dose values of each FID points.
Assuming that the mAs values of high pressure generator are fixed, determine fitting coefficient during arbitrary kV values, specifically comprise the following steps:
S131, using kV values as abscissa, fitting coefficient is ordinate, establishes coordinate system.
It is abscissa by kV1 ... kVN, N number of fitting coefficient A1 ... the AN ordinates obtained in above-mentioned steps are established
Coordinate system.
S132 will represent kV values one by one and the point of fitting coefficient be distributed in a coordinate system.
Each point of the kV values obtained in step S2 and fitting coefficient A compositions is brought into coordinate system.
S133, the point in coordinate system, carries out curve fitting, and obtains the functional relation of fitting coefficient A and kV value, passes through
KV is worth to fitting coefficient A.
Point in coordinate system, carries out curve fitting, and obtains the functional relation of fitting coefficient A and kV value, is expressed as A=
f(kV).At this time using the formula fitted, it may be determined that A coefficients during arbitrary kV values.
After this step, the second matched curve apart from bulb focal position and incident dose, formula y=f can be obtained
(kV)×x-u.KV is the kV values of high pressure generator setting, and x is apart from the position of bulb focus.It can be calculated by this formula
Arbitrary kV under the conditions of R1mAs (mAs values be R1) is in the incident dose value of FIDi points.
S14, according to the relation of the mAs values of high pressure generator and incident dose, when obtaining arbitrary kV values, arbitrary mAs values away from
The 3rd matched curve from bulb focal position and incident dose.
In actual use, the mAs values of high pressure generator are also adjustable within the scope of one that can be counted after above step
Incident dose when calculating arbitrary FID, arbitrary kV values under R1mAs values (mAs values are R1), can not use whole mAs values.So
It needs y=f (kV) × x-uIt is extended.Because incident dose and mAs values are directly proportional, assume according to this principle real
Border is R using adjusting to arbitrary mAs values, at this point, when obtained arbitrary kV values, arbitrary mAs values, apart from bulb focal position with entering
The 3rd matched curve for penetrating dosage is:Y=(R/R1) × f (kV) × x-u, arbitrary kV can be calculated by the 3rd matched curve
When value, arbitrary mAs values condition, the incident dose apart from bulb focus any position.
Patient when carrying out image check, is inputted the 3rd matched curve apart from the position of bulb focus, obtains X and penetrate by S15
The incident dose of line.
S16 according to the mAs values of high pressure generator and the relationship error of incident dose, corrects the 3rd matched curve,
Obtain final matched curve.
It is found in actual tests, there are certain error, this error should for mAs values and incident dose proportional relation
Caused by being exactly the mAs value accuracy differences of different generators.In order to correct this error, M parameter dosage multiplying power factors are introduced.
According to the mAs values of application mesohigh generator and the relationship error of incident dose, the 3rd matched curve is corrected, is obtained most
Whole matched curve, specifically comprises the following steps:
S161, in fixed FID points position, the kV values of high pressure generator are arranged to kV1, and mAs values are arranged to R2;Make
The incident dose RB of this position is measured with dosemeter;
In fixed FID points position, the kV values of high pressure generator are arranged to kV1, and it is R2 to select a mAs value, is made
The incident dose of this FID point is measured with dosemeter, records this incident dose value as RB.
S162, in the fixed FID points position, the kV values of high pressure generator are arranged to kV1, and mAs values are arranged to R1;It uses
Dosemeter measures the incident dose RA of this position;
R1mAs values are RA in the incident dose value of this FID point when step has measured kV1 before setting.
S163 is calculated in FID points position, and mAs values are arranged to incident dose RBS values during R2;In RBS=(R2/R1) ×
RA;
(R2/R1) should be equal to according to incident dose RBS value of this FID point of theoretical calculation under R2mAs (mAs values are R2)
×RA。
S164 according to the incident dose that dosemeter measures with calculating the incident dose obtained, obtains parameter dosage multiplying power system
Number, brings the parameter dosage multiplying power factor into the 3rd matched curve, obtains final matched curve.
Actually measured RB values are at this time certainly M dosage multiplying power factors with RB/RBS, by this not equal to calculated value RBS
Value add in the 3rd fit curve equation can the calculation error brought of the accurate sex differernce of modifying factor mAs values, obtain arbitrary kV values, times
During mAs values of anticipating, the final matched curve apart from bulb focal position and incident dose:
Y=(R/R1) × M × f (kV) × x-u。
Wherein, y is incident dose, and y is calculated apart from bulb focal position x points, and mAs values are R, in this kV value
When incident dose value, the dosage unit that selects when dose value unit is the test of dosage instrument.R is the mAs under conditions present
Value, R1 are a reference value of mAs values, are a fixed values.F (kV) is fitting coefficient, according to measurement fitting determine using kV values as
The A coefficient formulas of variable.X is distance of the current location apart from bulb focus;U is function coefficients, and what when fitting formula determined consolidates
Definite value.M is to determine dosage multiplying power factor according to experiment, corrects the calculation error that mAs othernesses are brought.
It is the incident dose value being calculated using incident dose computational methods in surface provided by the present invention and to make below
The value tested with dosemeter compares, and calculation error meets requirement.
The incident dose value and the comparison sheet of dosage instrument test value that 4 part of table is calculated
The incident dose value and the comparison sheet of dosage instrument test value that 5 part of table is calculated pass through actual test and meter
Verification is calculated, can extremely accurate be calculated in ray path and appointed using incident dose computational methods in surface provided by the present invention
One apart from the corresponding incident dose in the position of bulb focus.If by the method in X ray image equipment, it can be without using
Existing dose sensor measures, you can obtains incident dose value, while can also provide incident dose value to doctor and carry
Show, so as to more effectively protect patient from receiving excess X-radiation radiation.
In conclusion incident dose computational methods in surface provided by the present invention, by the kV values for setting high pressure generator
With mAs values, change kV values, obtain N number of incident dose apart from bulb focus different position;According to it is N number of apart from bulb focus not
With the incident dose of position, first matched curve of the N items apart from bulb focal position and incident dose is obtained;Then, high pressure
The mAs values of raw device are fixed, and are determined fitting coefficient during arbitrary kV values, are obtained second apart from bulb focal position and incident dose
Matched curve.Finally, according to the relation of the mAs values of high pressure generator and incident dose, when obtaining arbitrary kV values, arbitrary mAs values,
The 3rd matched curve apart from bulb focal position and incident dose;When carrying out image check, by patient apart from bulb focus
Position input the 3rd matched curve, obtain the incident dose of X-ray.It can be used according to this method in X ray image equipment
When the surface incidence agent that patient received calculated according to distances of the kVmAs and patient's incidence surface of setting apart from bulb focus
Amount.This method can show the incident dose being calculated when patient puts position and dosage is set at any time not against dose sensor
Value, to doctor to prompt, can prevent incident dose overshoot risk.
The present invention also provides a kind of surface incident dose computing device, the surface incident dose to realize above-mentioned calculates
Method.As shown in figure 3, the equipment includes processor 22 and is stored with the memory 21 of 22 executable instruction of processor;
Wherein, processor 22 can be general processor, such as central processing unit (CPU), can also be at digital signal
Reason device (DSP), application-specific integrated circuit (ASIC) or be arranged to implement the embodiment of the present invention one or more integrate electricity
Road.
Wherein, which for storing program code, and is transferred to CPU by memory 21.Memory 21 can wrap
Include volatile memory, such as random access memory (RAM);Memory 21 can also include nonvolatile memory, such as
Read-only memory, flash memory, hard disk or solid state disk;Memory 21 can also include the group of the memory of mentioned kind
It closes.
Specifically, the surface incident dose computing device that the embodiment of the present invention is provided, including processor 22 and memory
21;The computer program run on processor 22 can be used by being stored on memory 21, when computer program is held by processor 32
Following steps are realized during row:
The kilovolt value of high pressure generator and milliampere seconds value are set, changes kilovolt value, obtains apart from bulb focus different position
Incident dose;
According to N number of incident dose apart from bulb focus different position, N items are obtained apart from bulb focal position and incident agent
First matched curve of amount;Wherein, N is positive integer;
Assuming that the milliampere seconds value of high pressure generator is fixed, determine fitting coefficient during arbitrary kilovolt value, obtain apart from bulb
Focal position and the second matched curve of incident dose;
According to the milliampere seconds value of high pressure generator and the relation of incident dose, arbitrary kilovolt value, arbitrary milliampere seconds value are obtained
When the 3rd matched curve apart from bulb focal position and incident dose;
When carrying out image check, patient is inputted into the 3rd matched curve apart from the position of bulb focus, obtains X-ray
Incident dose.
Wherein, the kV values of high pressure generator and mAs values are set, changes kV values, obtains apart from bulb focus different position
Incident dose;Following steps are realized when computer program is performed by processor 22;
S211 obtains the kV maximums of the kV values of high pressure generator and kV minimum values, kV maximums and kV minimum values it
Between choose N number of kV values, the N that is set to kVj, j=1,2 ..., N is positive integer.
S212 sets value R1 on the basis of the mAs values of high pressure generator.
S213, according to patient support to the distance of bulb focus, determine it is N number of apart from Along ent, apart from bulb focus
Position is respectively that FIDi, i=1,2 ... N, N are positive integer.
S214 sets the kV value kVj of high pressure generator, obtains the incident agent apart from bulb focus different position FIDi one by one
Amount.
S215 repeats step S214, and j is made to distinguish value 1,2 ... N, is obtained under each kV values, different apart from bulb focus
The incident dose of position.
Wherein, according to N number of incident dose apart from bulb focus different position, obtain N items apart from bulb focal position with
First matched curve of incident dose;Following steps are realized when computer program is performed by processor 22;
S221, using, as abscissa, the incident dose value of each point is ordinate, establishes coordinate apart from the position of bulb focus
System;
S222 when the kV values of high pressure generator are kVj, will be represented apart from the position of bulb focus and incident dose one by one
Point distribution is in a coordinate system;
S223, the point in coordinate system, carries out curve fitting, and obtains one apart from bulb focal position and incident dose
The first matched curve;
S224 repeats step S222~S223, and j difference value 1,2 ... N is made to carry out curve fitting, obtains N item distances
Bulb focal position and the first matched curve of incident dose.
Wherein, the mAs values of high pressure generator are fixed, and determine fitting coefficient during arbitrary kV values;When computer program is located
Reason device 22 realizes following steps when performing;
Using kV values as abscissa, fitting coefficient is ordinate, establishes coordinate system;
KV values will be represented one by one and the point of fitting coefficient is distributed in a coordinate system;
Point in coordinate system, carries out curve fitting, and obtains the functional relation of fitting coefficient A and kV value, passes through kV values
Obtain fitting coefficient A.
Wherein, following steps are also realized when computer program is performed by processor 22;
According to the mAs values of application mesohigh generator and the relationship error of incident dose, the 3rd matched curve is entangled
Just, final matched curve is obtained.
Wherein, according to application mesohigh generator mAs values and incident dose relationship error, to the 3rd matched curve into
Row is corrected, and obtains final matched curve;Following steps are realized when computer program is performed by processor 22;
In fixed FID points position, the kV values of high pressure generator are arranged to kV1, and mAs values are arranged to R2;Use agent
Amount instrument measures the incident dose RB of this position;
In the fixed FID points position, the kV values of high pressure generator are arranged to kV1, and mAs values are arranged to R1;Dosage
Instrument measures the incident dose RA of this position;
It calculates in FID points position, mAs values are arranged to incident dose RBS values during R2;In RBS=(R2/R1) × RA;
According to the incident dose that dosemeter measures and the incident dose for calculating acquisition, parameter dosage multiplying power factor is obtained, it will
The parameter dosage multiplying power factor brings the 3rd matched curve into, obtains final matched curve.
The present invention also provides a kind of computer readable storage mediums.Here computer-readable recording medium storage has one
A or multiple programs.Wherein, computer readable storage medium can include volatile memory, such as random access memory
Device;Memory can also include nonvolatile memory, such as read-only memory, flash memory, hard disk or solid state disk;It deposits
Reservoir can also include the combination of the memory of mentioned kind.When one or more program in computer readable storage medium can
It is performed by one or more processor, to realize the part steps in above-mentioned surface incident dose computational methods or all walk
Suddenly.
Incident dose computational methods in surface provided by the present invention, equipment and storage medium have been carried out specifically above
It is bright.For those of ordinary skill in the art, it is done on the premise of without departing substantially from true spirit any aobvious
And the change being clear to, it will all form to infringement of patent right of the present invention, corresponding legal liabilities will be undertaken.
Claims (10)
1. a kind of surface incident dose computational methods, it is characterised in that include the following steps:
The kilovolt value of high pressure generator and milliampere seconds value are set, changes kilovolt value, obtains entering apart from bulb focus different position
Penetrate dosage;
According to N number of incident dose apart from bulb focus different position, N items are obtained apart from bulb focal position and incident dose
First matched curve;Wherein, N is positive integer;
Assuming that the milliampere seconds value of high pressure generator is fixed, fitting coefficient during arbitrary kilovolt value is determined, obtain apart from bulb focus
Position and the second matched curve of incident dose;
According to the milliampere seconds value of high pressure generator and the relation of incident dose, when obtaining arbitrary kilovolt value, arbitrary milliampere seconds value away from
The 3rd matched curve from bulb focal position and incident dose;
When carrying out image check, patient is inputted into the 3rd matched curve apart from the position of bulb focus, obtains the incidence of X-ray
Dosage.
2. incident dose computational methods in surface as described in claim 1, it is characterised in that the kilovolt value of high pressure generator is set
With milliampere seconds value, change kilovolt value, the incident dose apart from bulb focus different position is obtained, including following sub-step:
S111 obtains the kilovolt maximum of the kilovolt value of high pressure generator and kilovolt minimum value, in kilovolt maximum and kilovolt most
N number of kilovolt of value is chosen between small value, the N that is set to kVj, j=1,2 ...;
S112 sets value R1 on the basis of the milliampere seconds value of high pressure generator;
S113, according to patient support to the distance of bulb focus, determine it is N number of apart from Along ent, apart from the position of bulb focus
Respectively FIDi, i=1,2 ... N, N are positive integer;
S114 sets the kilovolt value kVj of high pressure generator, obtains the incident agent apart from bulb focus different position FIDi one by one
Amount;
S115 repeats step S114, and j is made to distinguish value 1,2 ... N, is obtained under each kilovolt value, apart from bulb focus difference position
The incident dose put.
3. incident dose computational methods in surface as described in claim 1, it is characterised in that according to N number of different apart from bulb focus
The incident dose of position obtains first matched curve of the N items apart from bulb focal position and incident dose, including following sub-step
Suddenly:
S121, using, as abscissa, the incident dose value of each point is ordinate, establishes coordinate system apart from the position of bulb focus;
S122 when the kilovolt value of high pressure generator is kVj, will be represented apart from the position of bulb focus and the point of incident dose one by one
Distribution is in a coordinate system;
S123, the point in coordinate system, carries out curve fitting, and obtains one article apart from the of bulb focal position and incident dose
One matched curve;
S124 repeats step S122~S123, and j difference value 1,2 ... N is made to carry out curve fitting, obtains N items apart from bulb
Focal position and the first matched curve of incident dose.
4. incident dose computational methods in surface as described in claim 1, it is characterised in that assuming that the milliampere-second of high pressure generator
Value is fixed, and fitting coefficient during arbitrary kilovolt value is determined, including following sub-step:
Using kilovolt value as abscissa, fitting coefficient is ordinate, establishes coordinate system;
The point distribution of kilovolt value and fitting coefficient will be represented in a coordinate system one by one;
Point in coordinate system, carries out curve fitting, and obtains the functional relation of fitting coefficient and kilovolt value, is worth by kilovolt
To fitting coefficient.
5. incident dose computational methods in surface as described in claim 1, it is characterised in that calculate X-ray incident dose it
Before, further include following sub-step:
According to the milliampere seconds value of application mesohigh generator and the relationship error of incident dose, the 3rd matched curve is entangled
Just, final matched curve is obtained.
6. incident dose computational methods in surface as claimed in claim 5, it is characterised in that according to application mesohigh generator
The relationship error of milliampere seconds value and incident dose, corrects the 3rd matched curve, final matched curve is obtained, including as follows
Sub-step:
In fixed FID points position, the kilovolt value of high pressure generator is arranged to kV1, and milliampere seconds value is arranged to R2;Use agent
Amount instrument measures the incident dose RB of this position;
In the fixed FID points position, the kilovolt value of high pressure generator is arranged to kV1, and milliampere seconds value is arranged to R1;Dosage
Instrument measures the incident dose RA of this position;
It calculates in FID points position, milliampere seconds value is arranged to incident dose RBS values during R2;In RBS=(R2/R1) × RA;
According to the incident dose that dosemeter measures with calculating the incident dose obtained, parameter dosage multiplying power factor is obtained, this is joined
Number dosage multiplying power factor brings the 3rd matched curve into, obtains final matched curve.
7. incident dose computational methods in surface as claimed in claim 5, it is characterised in that the arbitrary kilovolt value, arbitrary milliampere
Final matched curve during seconds value apart from bulb focal position and incident dose is represented using equation below:
Y=(R/R1) × M × f (kV) × x-u;
Wherein, y is incident dose, and R is the milliampere seconds value under conditions present, and R1 is a reference value of milliampere seconds value, and f (kV) is fitting
Coefficient;X is current location apart from the position of bulb focus;U is function coefficients, and M is parameter dosage multiplying power system.
8. a kind of surface incident dose computing device, it is characterised in that including processor and memory;It is stored on the memory
There is the available computer program run on the processor, realized such as when the computer program is performed by the processor
Lower step:
The kilovolt value of high pressure generator and milliampere seconds value are set, changes kilovolt value, obtains entering apart from bulb focus different position
Penetrate dosage;
According to N number of incident dose apart from bulb focus different position, N items are obtained apart from bulb focal position and incident dose
First matched curve;Wherein, N is positive integer;
Assuming that the milliampere seconds value of high pressure generator is fixed, fitting coefficient during arbitrary kilovolt value is determined, obtain apart from bulb focus
Position and the second matched curve of incident dose;
According to the milliampere seconds value of high pressure generator and the relation of incident dose, when obtaining arbitrary kilovolt value, arbitrary milliampere seconds value away from
The 3rd matched curve from bulb focal position and incident dose;
When carrying out image check, patient is inputted into the 3rd matched curve apart from the position of bulb focus, obtains the incidence of X-ray
Dosage.
9. incident dose computing device in surface as claimed in claim 8, it is characterised in that when the computer program is described
When processor performs, following steps are also realized;
According to the milliampere seconds value of application mesohigh generator and the relationship error of incident dose, the 3rd matched curve is entangled
Just, final matched curve is obtained.
10. a kind of computer readable storage medium, which is characterized in that the computer-readable recording medium storage there are one or
Multiple programs, one or more of programs can be performed by one or more processor, to realize in claim 1~7
The step in the incident dose computational methods of surface described in any one.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711436101.3A CN108074636B (en) | 2017-12-26 | 2017-12-26 | Surface incident dose calculation method, device and storage medium |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711436101.3A CN108074636B (en) | 2017-12-26 | 2017-12-26 | Surface incident dose calculation method, device and storage medium |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108074636A true CN108074636A (en) | 2018-05-25 |
CN108074636B CN108074636B (en) | 2021-10-26 |
Family
ID=62155804
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711436101.3A Active CN108074636B (en) | 2017-12-26 | 2017-12-26 | Surface incident dose calculation method, device and storage medium |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108074636B (en) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005074000A (en) * | 2003-09-01 | 2005-03-24 | Ge Medical Systems Global Technology Co Llc | X-ray ct apparatus and exposure dose distribution calculating method |
US7309615B2 (en) * | 2002-09-27 | 2007-12-18 | Siemens Medical Solutions Diagnostic | High quantum yield acridinium compounds and their uses in improving assay sensitivity |
CN101985739A (en) * | 2010-10-21 | 2011-03-16 | 复旦大学 | Dosage control system of plasma immersion implanting device |
WO2013049455A1 (en) * | 2011-09-29 | 2013-04-04 | William Blakely | Biodosimetry panels and methods |
CN103340646A (en) * | 2013-07-23 | 2013-10-09 | 奚岩 | Method for calculating, displaying and storing radiation dosages of CT image formation |
CN105030266A (en) * | 2014-04-21 | 2015-11-11 | 株式会社东芝 | X-ray computer tomographic apparatus and scan plan setting supporting apparatus |
CN105664376A (en) * | 2015-12-31 | 2016-06-15 | 上海联影医疗科技有限公司 | Method and device of determining dose distribution and optimizing radiotherapy treatment plan |
CN106232007A (en) * | 2014-03-14 | 2016-12-14 | 株式会社日立制作所 | X ray CT device and processing means |
CN106662656A (en) * | 2014-08-26 | 2017-05-10 | 三菱电机株式会社 | Dose rate measurement device |
-
2017
- 2017-12-26 CN CN201711436101.3A patent/CN108074636B/en active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7309615B2 (en) * | 2002-09-27 | 2007-12-18 | Siemens Medical Solutions Diagnostic | High quantum yield acridinium compounds and their uses in improving assay sensitivity |
JP2005074000A (en) * | 2003-09-01 | 2005-03-24 | Ge Medical Systems Global Technology Co Llc | X-ray ct apparatus and exposure dose distribution calculating method |
CN101985739A (en) * | 2010-10-21 | 2011-03-16 | 复旦大学 | Dosage control system of plasma immersion implanting device |
WO2013049455A1 (en) * | 2011-09-29 | 2013-04-04 | William Blakely | Biodosimetry panels and methods |
CN103340646A (en) * | 2013-07-23 | 2013-10-09 | 奚岩 | Method for calculating, displaying and storing radiation dosages of CT image formation |
CN106232007A (en) * | 2014-03-14 | 2016-12-14 | 株式会社日立制作所 | X ray CT device and processing means |
CN105030266A (en) * | 2014-04-21 | 2015-11-11 | 株式会社东芝 | X-ray computer tomographic apparatus and scan plan setting supporting apparatus |
CN106662656A (en) * | 2014-08-26 | 2017-05-10 | 三菱电机株式会社 | Dose rate measurement device |
CN105664376A (en) * | 2015-12-31 | 2016-06-15 | 上海联影医疗科技有限公司 | Method and device of determining dose distribution and optimizing radiotherapy treatment plan |
Non-Patent Citations (1)
Title |
---|
陈晓飞 等: ""数字X 线摄影焦点至探测器距离"", 《中国临床研究》 * |
Also Published As
Publication number | Publication date |
---|---|
CN108074636B (en) | 2021-10-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Sabol | A Monte Carlo estimation of effective dose in chest tomosynthesis | |
Andria et al. | Towards patient dose optimization in digital radiography | |
Zarb et al. | Image quality assessment tools for optimization of CT images | |
Balter | Methods for measuring fluoroscopic skin dose | |
McKenney et al. | Experimental validation of a method characterizing bow tie filters in CT scanners using a real‐time dose probe | |
Lin et al. | Accuracy and calibration of integrated radiation output indicators in diagnostic radiology: a report of the AAPM Imaging Physics Committee Task Group 190 | |
Gislason et al. | Dose optimization in pediatric cardiac x‐ray imaging | |
CN105138813B (en) | Patient data is mapped as another system and method from a kind of physiological status | |
Suliman et al. | Survey of computed tomography technique and radiation dose in Sudanese hospitals | |
Andria et al. | Dose optimization in chest radiography: System and model characterization via experimental investigation | |
Richard et al. | Quantitative breast tomosynthesis: From detectability to estimability | |
Dong et al. | Relationship between x‐ray illumination field size and flat field intensity and its impacts on x‐ray imaging | |
JP6339684B2 (en) | Photon counting CT apparatus and estimated dose calculation method | |
Lee et al. | Determination of the reliability and repeatability of a quantitative occlusal analyzer by using a piezoelectric film sensor: An in vitro study | |
Hernandez et al. | Average glandular dose coefficients for pendant‐geometry breast CT using realistic breast phantoms | |
Borrego et al. | A hybrid phantom system for patient skin and organ dosimetry in fluoroscopically guided interventions | |
CN102621576B (en) | Mammary dose distribution detection method and device for special CT (computed tomography) equipment for mammary gland | |
Gatt et al. | Optimisation of the AP abdomen projection for larger patient body thicknesses | |
Fajrin et al. | Kilovolt peak meter design as a calibrator of X-ray machine. | |
Shim et al. | Radiation dose estimates based on Monte Carlo simulation for spiral breast computed tomography imaging in a large cohort of patients | |
CN104207794B (en) | Mammography multi-parameter sensor array and measuring method thereof | |
CN108074636A (en) | A kind of surface incident dose computational methods, equipment and storage medium | |
Figl et al. | Factors for conversion between human and automatic read‐outs of CDMAM images | |
Asada et al. | Proposal for local diagnostic reference levels in general radiography in Japan | |
Li et al. | Advancing size‐specific dose estimates in CT examinations: Dose estimates at longitudinal positions of scans |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |