CN102735169A - System and method for calculating width of dose curve band - Google Patents
System and method for calculating width of dose curve band Download PDFInfo
- Publication number
- CN102735169A CN102735169A CN2011100815808A CN201110081580A CN102735169A CN 102735169 A CN102735169 A CN 102735169A CN 2011100815808 A CN2011100815808 A CN 2011100815808A CN 201110081580 A CN201110081580 A CN 201110081580A CN 102735169 A CN102735169 A CN 102735169A
- Authority
- CN
- China
- Prior art keywords
- dose curve
- image
- curve band
- pixel
- gray
- 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
Images
Landscapes
- Image Processing (AREA)
Abstract
The invention relates to the field of medical devices. The invention provides a system for calculating widths of dose curve bands, a method for calculating the widths of the dose curve bands, a computer program for implementing the method, and a computer readable recording medium for storing the program. The system comprises: a pixel determining module for determining a longitudinal pixel number of more than one dose curve bands in a scan image, wherein the scan image has the more than one dose curve bands; and a width calculating module for calculating the widths of the dose curve bands according to the longitudinal pixel number. The objective of the system and method is to solve the drawback of big personal errors in the prior art, and improve the accuracy of dose curve width measurements.
Description
Technical field
The present invention relates to the imaging of medical field, particularly a kind of system and method that calculates the dose curve bandwidth.
Background technology
X-ray imaging equipment, (Computed Tomography, CT) equipment must carry out the dose curve bandwidth to equipment and measure, so that it meets the requirement of State Food and Drug Administration before dispatching from the factory like computer tomography.In the prior art, normally utilize survey instruments such as vernier caliper directly the width of dose curve band on the film to be carried out manual measurement.Notification number is that the utility model of CN86207064U just discloses a kind of multidimensional rhombus chi, and it is used for measuring or writing down a certain curve, as makes the isodose curve superposition.
But the personal error of said method is bigger.This is to form as follows because of the dose curve band: at first start CT scan, adjustment collimating apparatus opening makes X-ray beam project on a certain position of film through a collimating apparatus opening, on film, to form a dose curve band; Then change the opening of collimating apparatus, and by the technician by the artificial film position of moving of marking light, X-ray beam is projected on another position of film through the collimating apparatus that changes opening, overlapping on film with the zone of avoiding double exposing.Because the artificial film that moves unavoidably can produce deviation, in addition, identification light also has ± error of 12mm, so in said process, just produced personal error.And when measuring the width of dose curve band, the technician utilizes vernier caliper to confirm the width of dose curve through naked eyes usually, is confirming width border and the process that reads vernier caliper numerical value, just inevitably produces personal error.Therefore, prior art is difficult to accurately and targetedly assess the dosage characteristic of CT scanner, also can't information very accurately be provided to the client.
Summary of the invention
The object of the present invention is to provide a kind of system and method that calculates the dose curve bandwidth accurately, with the big defective of personal error in the solution prior art, thus the accuracy of improving agent discharge curve width measure.Simultaneously, the invention allows for corresponding computer programs and computer-readable medium.
In view of this; The present invention proposes a kind of system that is used to calculate the dose curve bandwidth; Said system comprises: the pixel determination module, be used for confirming vertical pixel count of an above dose curve band of one scan image, and wherein said scan image has an above dose curve band; The width computing module is used for the width according to said vertical this dose curve band of pixel count calculating.
According to one embodiment of present invention, said pixel determination module comprises: the dose curve band is confirmed submodule, is used for confirming the longitudinal boundary and the number of said dose curve band; The angle of inclination calculating sub module is used at each dose curve band, vertical coboundary of this dose curve band with vertically choose two frontier points respectively to form a quadrilateral on the lower boundary, calculate this tetragonal angle of inclination; Pixel is confirmed submodule, is used for confirming according to said angle of inclination vertical pixel count of this dose curve band.So that confirm vertical pixel count of each dose curve band in the said scan image, and can be according to the number preliminary identification result's of the present invention of dose curve band accuracy.
According to another embodiment of the invention, said dose curve band confirms that submodule comprises: region of interest is confirmed the unit, is used for calculating the longitudinal extent and the lateral extent of said scan image exposure area pixel, and confirms a region of interest in view of the above; Dose curve band number acquiring unit is used for obtaining the number of said region of interest dose curve band.So that calculate the number of dose curve band.
According to still another embodiment of the invention, said tetragonal center overlaps with the center of said dose curve band.So that calculate tetragonal angle of inclination.
Further, said pixel determination module comprises: the Flame Image Process submodule, if said scan image is coloured image, then convert thereof into bianry image, and with this bianry image as said scan image.So that with discriminating dose curve band and the film that is not made public.
Said according to still another embodiment of the invention Flame Image Process submodule comprises: the gray level image acquiring unit is used for converting said coloured image to gray level image; The bianry image acquiring unit is used to cut apart said gray level image to obtain bianry image.
Further, said Flame Image Process submodule comprises: smooth unit, be used for said gray level image is carried out smothing filtering, and image that will be behind smothing filtering is as said gray level image.
According to still another embodiment of the invention, said bianry image acquiring unit comprises: grey level histogram obtains subelement, is used to obtain the grey level histogram of said gray level image; Level and smooth subelement is used for said grey level histogram is carried out smothing filtering; The frequency spike value is searched subelement, is used for finding out two frequency spike values of the grey level histogram behind smothing filtering; Threshold value is obtained subelement, is used for minimum frequency value corresponding gray between these two frequency spike values as a threshold value; Bianry image obtains subelement, is used for the gray-scale value of more said each pixel of gray level image and the size of this threshold value, if the gray-scale value of a pixel is greater than this threshold value, then this pixel belongs to said dose curve band and its gray-scale value is changed to 1; Otherwise this pixel does not belong to said dose curve band and its gray-scale value is changed to 0, thereby obtains a width of cloth bianry image.
The present invention also provides a kind of method of calculating the dose curve bandwidth, and said method comprises: confirm vertical pixel count of an above dose curve band in the scan image, wherein said scan image has an above dose curve band; Calculate the width of this dose curve band according to said vertical pixel count.
According to one embodiment of present invention, said definite pixel comprises: longitudinal boundary and the number of confirming said dose curve band; In each dose curve band, vertical coboundary of this dose curve band with vertically choose two frontier points respectively to form a quadrilateral on the lower boundary, calculate this tetragonal angle of inclination; Confirm the vertical pixel count in this dose curve band according to said angle of inclination.So that confirm vertical pixel count of each dose curve band in the said scan image, and can be according to the number preliminary identification result's of the present invention of dose curve band accuracy.
According to another embodiment of the invention, the dose curve band in said definite said scan image comprises: calculate the longitudinal extent and the lateral extent of exposure area pixel in institute's scan image, and confirm a region of interest in view of the above; Obtain the number of dose curve band in the said region of interest.
Further, saidly confirm that vertical pixel count comprises: if said scan image is a coloured image, then convert thereof into gray level image earlier, cut apart said gray level image again to obtain a width of cloth bianry image, with said bianry image as said scan image.
Further, said method comprises: said gray level image is carried out smothing filtering, and image that will be behind smothing filtering is as said gray level image.
According to still another embodiment of the invention, the said gray level image of cutting apart comprises to obtain a width of cloth bianry image: the grey level histogram that obtains said gray level image; Said grey level histogram is carried out smothing filtering; Find out two frequency spike values in the grey level histogram behind smothing filtering; With minimum frequency value corresponding gray between these two frequency spike values as a threshold value; The size of the gray-scale value of each pixel and this threshold value in the more said gray level image, if the gray-scale value of a pixel is greater than this threshold value, then this pixel belongs to said dose curve band and its gray-scale value is changed to 1; Otherwise this pixel does not belong to said dose curve band and its gray-scale value is changed to 0, thereby obtains a width of cloth bianry image.
According to still another embodiment of the invention; Calculate the width of dose curve band according to following formula:
wherein; W is the width of a dose curve band; PN is vertical pixel count of this dose curve band, and DPI is the pixel count of per inch.
The present invention also provides a kind of computer program, makes computing machine carry out the method for said calculating dose curve bandwidth.
The present invention also provides a kind of computer-readable recording medium, has stored aforementioned calculation machine program.
From such scheme, can find out; Because system and method for the present invention can calculate the width of dose curve band automatically, thus the personal error in the measuring process avoided, thus improved the degree of accuracy of measuring the dose curve bandwidth; For example; If the resolution of scanner is 600DPI (dot per inch, the pixel count of per inch), the degree of accuracy of then measuring can reach 0.04mm.Simultaneously, system and method for the present invention can discern in the X-ray scanning process with scanner scanning film process in the angle of inclination of dose curve band and film, to carry out slant correction, this has further improved the degree of accuracy of measurement.In addition, system and method for the present invention can be used for various types of x-ray imaging equipment to measure the width of dose curve band.
Description of drawings
To make clearer above-mentioned and other feature and advantage of the present invention of those skilled in the art through describing embodiments of the invention in detail below with reference to accompanying drawing, in the accompanying drawing:
Fig. 1 is the synoptic diagram of the gray level image that converts to from original color image.
Fig. 2 is the grey level histogram of the gray level image behind smothing filtering.
Fig. 3 is for through cutting apart the synoptic diagram of the bianry image that gray-scale map obtains.
Fig. 4 is the synoptic diagram of the various inclinations that occur in the CT scan process, and wherein Fig. 4 a is film position run-off the straight in the CT scan process, the scan image synoptic diagram that obtains; Run-off the straight when Fig. 4 b is the scanner scanning film; The scan image synoptic diagram that obtains; Fig. 4 c is the run-off the straight of CT scan process film, and liquid run-off the straight during the scanner scanning film, the scan image synoptic diagram that obtains; Fig. 4 d is on the film image during porose or hole, the scan image synoptic diagram that obtains.
The tetragonal synoptic diagram that Fig. 5 forms for four frontier points in the dose curve band, getting, wherein Fig. 5 a is the quadrilateral synoptic diagram under the normal condition; Fig. 5 b is in the CT scan process during run-off the straight of film position, the quadrilateral synoptic diagram that obtains; Run-off the straight when Fig. 5 c is scanner scanning film process, the quadrilateral synoptic diagram that obtains; Fig. 5 d is the run-off the straight of CT scan process film, and also run-off the straight during the scanner scanning film, the quadrilateral synoptic diagram that obtains.
Fig. 6 is the width value synoptic diagram of the dose curve band that obtains with system and method for the present invention; Wherein Fig. 6 a is on the film image during porose or hole; The width value synoptic diagram of the dose curve band that obtains; When Fig. 6 b is film position or scanner run-off the straight, the width value synoptic diagram of the dose curve band that obtains; The width value synoptic diagram of the dose curve band under Fig. 6 c normal condition.
Fig. 7 is used to calculate the synoptic diagram of the system of dose curve bandwidth for the present invention.
Fig. 8 is the synoptic diagram of pixel determination module of the present invention.
Fig. 9 is the synoptic diagram of Flame Image Process submodule of the present invention.
Figure 10 is the synoptic diagram of bianry image acquiring unit of the present invention.
Figure 11 confirms the synoptic diagram of submodule for dose curve band of the present invention.
Figure 12 is for measuring the flow chart of steps of dose curve bandwidth based on system and method for the present invention.
Figure 13 calculates the method flow diagram of dose curve bandwidth for the present invention.
Reference numeral
1 the present invention calculates the system of dose curve bandwidth
12 pixel determination modules
13 width computing modules
120 Flame Image Process submodules
121 dose curve bands are confirmed submodule
122 angle of inclination calculating sub module
123 pixels are confirmed submodule
201 gray level image acquiring units
202 bianry image acquiring units
203 smooth unit
It is single that 301 grey level histograms obtain son
302 level and smooth subelements
303 frequency spike values are searched subelement
304 threshold values are obtained subelement
305 bianry images obtain subelement
211 region of interest are confirmed the unit
212 dose curve band number acquiring units
101～103 measure the step of dose curve bandwidth based on system and method for the present invention
201～206 the present invention are used to calculate the step of dose curve bandwidth
Embodiment
For making the object of the invention, technical scheme and advantage clearer, below lift embodiment to further explain of the present invention.
Below be example with CT equipment, specify system and method with calculating dose curve bandwidth of the present invention.
Shown in figure 12, for the step of measuring the dose curve bandwidth based on system and method for the present invention following:
Step 101 places film in the center measurement zone of X ray projection place on the examination couch.
Step 102 starts CT scan, and adjustment collimating apparatus opening projects on the film X-ray beam respectively through the collimating apparatus of different openings size, on film, to form the dose curve band of an above different in width.
Step 103 is handled this film with the system and method for calculating dose curve bandwidth of the present invention, to obtain the width value of each dose curve band.
In the prior art, through behind step 101～102, the technician just utilizes vernier caliper to come the width of various dose curve band on the manual measurement film, and there is bigger personal error in this as previously mentioned.And the present invention utilizes step 103 to substitute the step of manual measurement dose curve bandwidth in the prior art, thereby has realized the robotization and the precision of dose curve bandwidth measuring process, has significantly reduced personal error.
Below the process of calculating the dose curve bandwidth with system and method for the present invention in the step 103 is elaborated.Because dissimilar CT equipment has different scan setting to dose curve; Like tube current, tube voltage, collimating apparatus width etc.; So the present invention confirms the method for the width of each dose curve band and should have ubiquity; Various x-ray imaging equipment being used for, therefore should not be provided with the parameter relevant in the method with scan setting.
Be illustrated in figure 7 as the synoptic diagram that the present invention is used to calculate the system 1 of dose curve bandwidth, it comprises pixel determination module 12 and width computing module 13, particularly:
As shown in Figure 8; Pixel determination module 12 comprises that the dose curve band confirms that submodule 121, angle of inclination calculating sub module 122 and pixel confirm submodule 123; Further pixel determination module 12 also comprises Flame Image Process submodule 120; If said scanner is color scanner, then to convert the coloured image that scanning obtains to bianry image, and with this bianry image as said scan image.If scanner is the black and white scanner, then scan the image that obtains and be black white image, then need not Flame Image Process submodule 120 and handle.
Wherein Flame Image Process submodule 120 comprises: gray level image acquiring unit 201 and bianry image acquiring unit 202, and as shown in Figure 9, particularly:
Gray level image acquiring unit 201 is used for the coloured image I that obtains with color scanner scanning
0Convert gray level image I to
g
Be illustrated in figure 1 as the synoptic diagram of the gray level image that converts to from original color image.8 gray level images with routine are example, and the scope of gray-scale value can be got [0,255] among Fig. 1, can find out gray level image I among the figure
gIn 7 dose curves (shown in the zone that indicates with oblique line among the figure) are arranged, the width of every dose curve band is different, this is because the collimating apparatus of the corresponding different openings size of the width of every dose curve band; The zone of using vertical line to indicate is not by the film of X ray exposure.In the Fig. 1 with the gray-scale value signal, system is difficult to dose curve band and the film strictness that is not made public are distinguished.
Owing to all possibly have some impurity on film and the scanner, also can produce noise in the scanner scanning film process, these all can make image I
0Not " totally ", promptly there is burr in the image, therefore through I
0Convert I to
gAlso can there be burr.
In order to eliminate these unnecessary burrs, further, Flame Image Process submodule 120 also comprises smooth unit 203, is used for said gray level image I
gCarry out smothing filtering, and image I that will be behind smothing filtering
sAs said gray level image I
gThe other types wave filter that can know with Gaussian filter or those skilled in the art is to gray level image I
gCarry out smothing filtering.If use Gaussian filter, available gaussian kernel function is shown in formula (1):
Wherein, (n1 n2) is gaussian kernel function to h, and n1 is the horizontal ordinate of each pixel in the two dimensional image, and n2 is the ordinate of each pixel in the two dimensional image, and σ is the standard deviation of Gaussian distribution.
The gaussian filtering kernel function of available 3*3 is pressed following formula (2) to gray level image I
gCarry out smothing filtering, obtain image I
s
I
s＝I
g*h(n1，n2) (2)
Wherein, I
sBe the gray level image behind smothing filtering, here with the gaussian filtering kernel function of 3*3 in order to remove image I effectively
gIn singular point.
Bianry image acquiring unit 202; Be used to cut apart said gray level image to obtain bianry image; It comprises that grey level histogram obtains subelement 301, level and smooth subelement 302, frequency spike value search subelement 303, threshold value obtains subelement 304 and bianry image obtains subelement 305, and is shown in figure 10.Because the pixel value of binary picture has only two, thereby is easy to the dose curve zone of unexposed film with exposure separated.Wherein:
Grey level histogram obtains subelement 301, is used to obtain gray level image I
sGrey level histogram Hist, gray-scale value scope desirable [0,255].
Level and smooth subelement 302 is used for grey level histogram Hist is carried out smothing filtering, to remove the noise among the Hist, obtains Hist
sCan carry out smothing filtering to gray level image Hist with the other types wave filter that Gaussian filter or those skilled in the art know.If use Gaussian filter, available gaussian kernel function is shown in formula (3):
Wherein, Conv_Gauss is the Gaussian filter that is used for grey level histogram Hist is carried out smothing filtering; Nwide representes full width at half maximum (the Full width at half maxium of Gauss normal distribution curve; FWHM), Nsigma representes the scope of Gaussian distribution, and i representes each gray values of pixel points among the Hist.
So just can use the Gaussian filter in the formula (3) grey level histogram Hist to be carried out smothing filtering, to obtain Hist by following formula (4)
s
Wherein, k is an intermediate variable that carries out convolutional calculation, Hist
sBe the grey level histogram behind smothing filtering, as shown in Figure 2.Horizontal ordinate is a gray-scale value among Fig. 2, and ordinate is a frequency.
The frequency spike value is searched subelement 303, is used for finding out Fig. 2 Hist
sTwo frequency spike values.
At first, find out among Fig. 2 [0,255-Th
1] in frequency spike value PV
1, with PV
1Corresponding gray is made as GV
1Threshold T h wherein
1Be used to remove the histogrammic data of background.PV
1Gray-scale value GV
1Away from Th
1
Then, find out among Fig. 2 [0, GV
1-Th
2] and [GV
1-Th
2, 255-Th
1] in frequency spike value PV
2, with PV
2Corresponding gray is made as GV
2Threshold T h wherein
2Be used to find second kurtosis.Those skilled in the art can confirm Th according to experiment
1And Th
2Value, its result has robustness, does not detail here and confirms Th
1And Th
2Process.
Threshold value is obtained subelement 304, is used to find out PV
1With PV
2Between minimum frequency value (being the concave point value) VV, and obtain VV corresponding gray GV according to formula (5)
3, and with GV
3As a threshold value, wherein GV
3Be positioned at GV
1With GV
2Between.
VV＝Hist
s(GV
3) (5)
With GV
3Be made as the global threshold of Fig. 2, like this, utilize GV
3Just can pixels all in Fig. 1 gray level image be divided into two groups: every gray-scale value is higher than threshold value GV
3Pixel, this gray values of pixel points is made as 0 (black), every gray-scale value is lower than threshold value GV
3Pixel, this gray values of pixel points is made as 255 (whites).So obtain bianry image as shown in Figure 3, system just is easy to and can distinguishes the film that do not made public and the dose curve band of exposure according to Fig. 3 like this.
The coloured image that color scanner is scanned is after the processing of Flame Image Process submodule 120; Just can this scan image be sent into the dose curve band and confirm that submodule 121 handles; And the black and white scanner scanning to image owing to be black white image; Then need not to handle, directly send into the dose curve band and confirm that submodule 121 handles through Flame Image Process submodule 120.
The dose curve band is confirmed submodule 121, is used for confirming the longitudinal boundary and the number of the dose curve band of said scan image, and it comprises that region of interest confirms unit 211 and dose curve band number acquiring unit 212, and is shown in figure 11.Particularly:
Region of interest is confirmed unit 211, is used to calculate said scan image (bianry image I
BwOr black white image) longitudinal extent Rc and the lateral extent Rr of exposure area pixel in, with obtain the definite region of interest of Rc and Rr (Region ofInterest, ROI).
Dose curve band number acquiring unit 212 is used for obtaining the number of ROI dose curve band.
At first, calculate vertical value Vc of intermediate pixel among the Rc.
Then, in ROI, remove the singular point in the image with the window of 3*3; In order to avoid system's mistake is grouped into singular point in the dose curve band of exposure; Thereby make system can confirm the number of dose curve band exactly, and near Vc, longitudinally mark n bar straight line, wherein n is an adjustable parameter.
Angle of inclination calculating sub module 122; Confirm the dose curve band that submodule 121 is confirmed according to the dose curve band; In each dose curve band; Choose two some A, B in vertical coboundary of this dose curve band, vertically choosing two frontier point C, D on the lower boundary, calculate the angle of inclination of quadrilateral ABCD.Said n can determine the width of AB.
Because what the present invention will measure is vertical width of dose curve band, so limit AB and limit CD must get the up-and-down boundary of dose curve band, and limit AC and limit BD then can be on the borders, the left and right sides of dose curve band.
Preferably, the center of quadrilateral ABCD overlaps with the center of the dose curve band at its place, so that calculate tetragonal angle of inclination.
Generally, quadrilateral ABCD is a rectangle, shown in Fig. 5 a.But because the film position maybe run-off the straight in the CT scan process, also maybe run-off the straight during the scanner scanning film, as shown in Figure 4, wherein Fig. 4 a is film position run-off the straight in the CT scan process, the scan image synoptic diagram that obtains; Run-off the straight when Fig. 4 b is the scanner scanning film; The scan image synoptic diagram that obtains; Fig. 4 c is the run-off the straight of CT scan process film, and also run-off the straight during the scanner scanning film, the scan image synoptic diagram that obtains; Fig. 4 d is porose or when hole (pointing out said hole or hole with arrow among Fig. 4 d) on the film image, the scan image synoptic diagram that obtains.Fig. 4 d mesopore or hole are because the technician who has likes plateholder is produced at file, and the skill people personnel that have in addition like on film, doing sign (like literal) etc.System of the present invention can discern these holes, hole, sign etc., and gets rid of them to measuring the influence of dose curve bandwidth.Said n bar ordinate just can be chosen without the ordinate of these holes, hole and sign.
In Fig. 4 a～4c, quadrilateral ABCD is a rectangle no longer just, but parallelogram is as shown in Figure 5.Wherein Fig. 5 b is the quadrilateral synoptic diagram that obtains during the run-off the straight of film position in the CT scan process; The quadrilateral synoptic diagram that run-off the straight obtained when Fig. 5 c was scanner scanning film process; Fig. 5 d is the run-off the straight of CT scan process film, and also run-off the straight during the scanner scanning film, the quadrilateral synoptic diagram that obtains.
The present invention will identify in the CT scan process the whether run-off the straight of film position, and whether run-off the straight during the scanner scanning film, and this can calculate the tilt angle alpha of quadrilateral ABCD through formula (6)～(7), and β realizes.Wherein, α is the angle between AC and the x axle, and β is the angle between AC and the CD.The x direction is a horizontal direction among Fig. 5 a～5d, and the y direction is vertical direction (being aforementioned vertical).
α＝arctan(dist
y(A，C)/dist
x(A，C))(6)
Wherein, dist
x(A C) is A point and the distance of C point on the x direction, dist
y(A C) is A point and the distance of C point on the y direction.
β-α＝arctan(dist
y(A，B)/dist
x(A，B))(7)
Wherein, dist
x(A B) is A point and the distance of B point on the x direction, dist
y(A B) is A point and the distance of B point on the y direction.
Because β=α+(β-α), so just can obtain angle beta according to formula (6)～(7).
Inclination type among Fig. 5 a～5d is as follows shown in (a)～(d):
Pixel is confirmed submodule 123, and according to tilt angle alpha, β confirms the vertical pixel count in each dose curve band.
According to tilt angle alpha, β judges that this dose curve band belongs to any in the above-mentioned inclination type (a)～(d), in view of the above this dose curve band is carried out angularity correction, thereby obtains between AB and the CD in the n bar ordinate pixel count PN on every ordinate
1PN
n, the average value P N (pix number of each vertical line) that gets this n pixel count is as the pixel count on arbitrary ordinate between AB and the CD.
So just just obtained the vertical pixel count PN in each dose curve band through pixel determination module 12.
Width computing module 13 is used for obtaining the width W that vertical pixel count PN calculates each dose curve band according to pixel determination module 12.
At first, calculate the big or small Pix_Size of each pixel according to following formula (8),
Size through the common per inch of the image of scanner scanning (inch) is 25.4mm, if 1DPI=300dot/inch, then the size of each pixel is:
So the width W of each dose curve can be calculated by following formula (9):
W＝Pix_Size×PN(9)。
As shown in Figure 6, be the width value synoptic diagram of the dose curve band that obtains with system and method for the present invention, wherein scheme the width value of the numerical value on right side for each corresponding dose curve band.Fig. 6 a is on the film during porose or hole, when the width value synoptic diagram of the dose curve band that obtains, Fig. 6 b are the run-off the straight of film position, and the width value synoptic diagram of the dose curve band that obtains; The width value synoptic diagram of the dose curve band under Fig. 6 c normal condition.Amplification has been illustrated to calculate the required quadrilateral ABCD of 11.85mm dose curve band on the right side of Fig. 6 c, and Fig. 6 a and Fig. 6 b similarly repeat no more here.
Further; The system that the present invention is used for calculating the dose curve bandwidth also comprises an authentication module (accompanying drawing is not shown); It is used to preserve the width value W of each dose curve band; And together be presented at a visualization interface comparing with said coloured image or black white image, with the accuracy of checking result of calculation of the present invention.
The bar number of dose curve band in bar number through the dose curve band that relatively obtains with system of the present invention and coloured image or the black white image, but the accuracy of preliminary identification measurement result of the present invention.
In addition, the present invention also comprises a kind of method of calculating the dose curve bandwidth, and it comprises the steps 201～203:
Step 202 is confirmed the dose curve band in the said scan image.This confirms existing specifying in the submodule 121 at the dose curve band, repeats no more here.
Further; Also comprise step 203, if said scanner is a color scanner, the coloured image that then earlier scanning is obtained converts gray level image to; Cut apart said gray level image again to obtain a width of cloth bianry image; As said scan image, this is existing specifying in Flame Image Process submodule 120, repeats no more here with said bianry image.
Step 204 in each dose curve band, is chosen two some A, B in vertical coboundary of this dose curve band, vertically choosing two frontier point C, D on the lower boundary, calculates the angle of inclination of quadrilateral ABCD.
Preferably, the center of quadrilateral ABCD overlaps with the center of the dose curve band at its place.
Since possibility run-off the straight in film position in the CT scan process, also run-off the straight during the scanner scanning film, so quadrilateral ABCD just no longer is the rectangle shown in Fig. 5 a, but the parallelogram shown in Fig. 5 b～5d, the present invention will identify said inclination.Can calculate the tilt angle alpha of quadrilateral ABCD through formula (6)～(7), β realizes that this is existing specifying in angle of inclination calculating sub module 122, repeats no more here.
Step 206 is by formula calculated the width W of each dose curve in (8)～(9) according to PN.
Further, the present invention's method of calculating the dose curve bandwidth also comprises:
Step 207 is preserved the width value of each dose curve band in the step 206, and with step 201 in the scan image that obtains of scanner together be presented at a visualization interface to compare, with the accuracy of checking result of calculation of the present invention.
The present invention relates to field of medical.The present invention provides a kind of system that is used to calculate the dose curve bandwidth, calculates the method for dose curve bandwidth, is used to realize the computer program of said method and the computer readable recording medium storing program for performing of storing this program; Wherein said system comprises the pixel determination module; Be used for confirming vertical pixel count of an above dose curve band of one scan image, wherein said scan image has an above dose curve band; The width computing module is used for the width according to said vertical this dose curve band of pixel count calculating.To solve the big defective of personal error in the prior art, the accuracy of improving agent discharge curve width measure.The above is merely preferred embodiment of the present invention, and is in order to restriction the present invention, not all within spirit of the present invention and principle, any modification of being done, is equal to replacement, improvement etc., all should be included within protection scope of the present invention.
Claims (17)
1. system (1) that is used to calculate the dose curve bandwidth, said system comprises:
Pixel determination module (12) is used for confirming vertical pixel count of an above dose curve band of one scan image, and wherein said scan image has an above dose curve band;
Width computing module (13) is used for the width according to said vertical this dose curve band of pixel count calculating.
2. system according to claim 1 is characterized in that, said pixel determination module (12) comprising:
The dose curve band is confirmed submodule (121), is used for confirming the longitudinal boundary and the number of said dose curve band;
Angle of inclination calculating sub module (122) is used at each dose curve band, vertical coboundary of this dose curve band with vertically choose two frontier points respectively to form a quadrilateral on the lower boundary, calculate this tetragonal angle of inclination;
Pixel is confirmed submodule (123), is used for confirming according to said angle of inclination vertical pixel count of this dose curve band.
3. system according to claim 2 is characterized in that, said dose curve band confirms that submodule (121) comprising:
Region of interest is confirmed unit (211), is used for calculating the longitudinal extent and the lateral extent of said scan image exposure area pixel, and confirms a region of interest in view of the above;
Dose curve band number acquiring unit (212) is used for obtaining the number of said region of interest dose curve band.
4. system according to claim 2 is characterized in that, said tetragonal center overlaps with the center of said dose curve band.
5. system according to claim 2 is characterized in that, said pixel determination module (12) further comprises:
Flame Image Process submodule (120) if said scan image is coloured image, then converts thereof into bianry image, and with this bianry image as said scan image.
6. system according to claim 5 is characterized in that, said Flame Image Process submodule (120) comprising:
Gray level image acquiring unit (201) is used for converting said coloured image to gray level image;
Bianry image acquiring unit (202) is used to cut apart said gray level image to obtain bianry image.
7. system according to claim 6, said Flame Image Process submodule (120) further comprises:
Smooth unit (203) is used for said gray level image is carried out smothing filtering, and image that will be behind smothing filtering is as said gray level image.
8. system according to claim 6 is characterized in that, said bianry image acquiring unit (202) comprising:
Grey level histogram obtains subelement (301), is used to obtain the grey level histogram of said gray level image;
Level and smooth subelement (302) is used for said grey level histogram is carried out smothing filtering;
The frequency spike value is searched subelement (303), is used for finding out two frequency spike values of the grey level histogram behind smothing filtering;
Threshold value is obtained subelement (304), is used for minimum frequency value corresponding gray between these two frequency spike values as a threshold value;
Bianry image obtains subelement (305), is used for the gray-scale value of more said each pixel of gray level image and the size of this threshold value, if the gray-scale value of a pixel is greater than this threshold value, then this pixel belongs to said dose curve band and its gray-scale value is changed to 1; Otherwise this pixel does not belong to said dose curve band and its gray-scale value is changed to 0, thereby obtains a width of cloth bianry image.
9. method of calculating the dose curve bandwidth, said method comprises:
Confirm vertical pixel count of an above dose curve band in the one scan image, wherein said scan image has an above dose curve band;
Calculate the width of this dose curve band according to said vertical pixel count.
10. method according to claim 9 is characterized in that, said definite vertical pixel count comprises:
Confirm the longitudinal boundary and the number of said dose curve band;
In each dose curve band, vertical coboundary of this dose curve band with vertically choose two frontier points respectively to form a quadrilateral on the lower boundary, calculate this tetragonal angle of inclination;
Confirm the vertical pixel count in this dose curve band according to said angle of inclination.
11. method according to claim 10 is characterized in that, the dose curve band in said definite said scan image comprises:
Calculate the longitudinal extent and the lateral extent of exposure area pixel in institute's scan image, and confirm a region of interest in view of the above;
Obtain the number of dose curve band in the said region of interest.
12. method according to claim 10 is characterized in that, said definite vertical pixel count further comprises:
If said scan image is coloured image, then with converting thereof into gray level image earlier, cut apart said gray level image again to obtain a width of cloth bianry image, with said bianry image as said scan image.
13. method according to claim 12, said method further comprises:
Said gray level image is carried out smothing filtering, and image that will be behind smothing filtering is as said gray level image.
14. method according to claim 12 is characterized in that, the said gray level image of cutting apart comprises to obtain a width of cloth bianry image:
Obtain the grey level histogram of said gray level image;
Said grey level histogram is carried out smothing filtering;
Find out two frequency spike values in the grey level histogram behind smothing filtering;
With minimum frequency value corresponding gray between these two frequency spike values as a threshold value;
The size of the gray-scale value of each pixel and this threshold value in the more said gray level image, if the gray-scale value of a pixel is greater than this threshold value, then this pixel belongs to said dose curve band and its gray-scale value is changed to 1; Otherwise this pixel does not belong to said dose curve band and its gray-scale value is changed to 0, thereby obtains a width of cloth bianry image.
15. method according to claim 9 is characterized in that, calculates the width of dose curve band according to following formula:
Wherein, W is the width of a dose curve band, and PN is vertical pixel count of this dose curve band, and DPI is the pixel count of per inch.
16. a computer program is carried out like each described method among the claim 9-15 computing machine.
17. a computer-readable recording medium has been stored computer program as claimed in claim 16.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201110081580.8A CN102735169B (en) | 2011-03-31 | 2011-03-31 | System and method for calculating width of dose curve band |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201110081580.8A CN102735169B (en) | 2011-03-31 | 2011-03-31 | System and method for calculating width of dose curve band |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102735169A true CN102735169A (en) | 2012-10-17 |
| CN102735169B CN102735169B (en) | 2015-06-03 |
Family
ID=46991107
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201110081580.8A Active CN102735169B (en) | 2011-03-31 | 2011-03-31 | System and method for calculating width of dose curve band |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102735169B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113533785A (en) * | 2021-07-22 | 2021-10-22 | 中国计量科学研究院 | Scanning electron microscope-based nano line width self-adaptive measurement method and system |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN86207064U (en) * | 1986-09-18 | 1987-09-30 | 肿瘤医院 | Multidimensional rhombus ruler |
| JP2006292922A (en) * | 2005-04-08 | 2006-10-26 | Fuji Photo Film Co Ltd | Method and apparatus for correcting radiation image |
| CN101691994A (en) * | 2009-09-30 | 2010-04-07 | 浙江大学 | Method for automatically positioning and detecting maximum width of tunnel crack |
-
2011
- 2011-03-31 CN CN201110081580.8A patent/CN102735169B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN86207064U (en) * | 1986-09-18 | 1987-09-30 | 肿瘤医院 | Multidimensional rhombus ruler |
| JP2006292922A (en) * | 2005-04-08 | 2006-10-26 | Fuji Photo Film Co Ltd | Method and apparatus for correcting radiation image |
| CN101691994A (en) * | 2009-09-30 | 2010-04-07 | 浙江大学 | Method for automatically positioning and detecting maximum width of tunnel crack |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113533785A (en) * | 2021-07-22 | 2021-10-22 | 中国计量科学研究院 | Scanning electron microscope-based nano line width self-adaptive measurement method and system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102735169B (en) | 2015-06-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10175181B2 (en) | X-ray imaging system allowing the correction of the scatter radiation and precise detection of the distance between the source and the detector | |
| US9177682B2 (en) | Dynamic filter for computed tomography (CT) | |
| US9934603B2 (en) | Three-dimensional resolution gauge for evaluating performance of tomographic imaging systems | |
| KR101670986B1 (en) | System and method for automatic planning of views in 3d images of brain | |
| EP3250909B1 (en) | Gap resolution for linear detector array | |
| CN104771184A (en) | CT scanning imaging method | |
| JPH07181609A (en) | Region of irradiation field extracting device for radiograph | |
| KR102138778B1 (en) | Test chart and method for using thereof | |
| CN104173071A (en) | Method and device for automatically determining scanning range of petrous bone | |
| US6839401B2 (en) | X-ray computed tomography apparatus | |
| US20160210740A1 (en) | Method and system for spine position detection | |
| EP0335204B1 (en) | Method for determining the contour of an irradiation field | |
| CN102735169B (en) | System and method for calculating width of dose curve band | |
| CN106651768B (en) | Image correction method and device and X-ray photographing equipment | |
| US20210212642A1 (en) | Methods and systems for beam attenuation | |
| CN102759538B (en) | Method for measuring and modeling point spread function of cone beam CT system | |
| CN104132950B (en) | CL scanning means projection rotating center scaling method based on original projection information | |
| US7111985B2 (en) | Method and system for measuring table sag | |
| WO2011080808A1 (en) | Radiograph processing device and radiograph processing program | |
| CN111862130B (en) | Body surface contour estimation method for medical imaging | |
| JP2005058315A (en) | Image size reduction factor calculation method | |
| Hasegawa et al. | Quality control of scintillation cameras using a minicomputer | |
| CN103505231A (en) | Method and system for adjusting semi-fan-shaped beam angle of computed tomography system and system thereof | |
| CN111696047B (en) | Imaging quality determining method and system of medical imaging equipment | |
| CN115439353B (en) | CT image ring artifact correction method, system and storage medium |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |