CN107990853B - A kind of source image distance detecting method and device - Google Patents

A kind of source image distance detecting method and device Download PDF

Info

Publication number
CN107990853B
CN107990853B CN201810005172.6A CN201810005172A CN107990853B CN 107990853 B CN107990853 B CN 107990853B CN 201810005172 A CN201810005172 A CN 201810005172A CN 107990853 B CN107990853 B CN 107990853B
Authority
CN
China
Prior art keywords
sid
frequency distribution
grid
main body
distributing position
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.)
Active
Application number
CN201810005172.6A
Other languages
Chinese (zh)
Other versions
CN107990853A (en
Inventor
李海春
王柳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Neusoft Medical Systems Co Ltd
Original Assignee
Neusoft Medical Systems Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Neusoft Medical Systems Co Ltd filed Critical Neusoft Medical Systems Co Ltd
Priority to CN201810005172.6A priority Critical patent/CN107990853B/en
Publication of CN107990853A publication Critical patent/CN107990853A/en
Application granted granted Critical
Publication of CN107990853B publication Critical patent/CN107990853B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B15/00Measuring arrangements characterised by the use of electromagnetic waves or particle radiation, e.g. by the use of microwaves, X-rays, gamma rays or electrons
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
    • G01N23/02Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
    • G01N23/04Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01TMEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T7/00Details of radiation-measuring instruments

Abstract

The embodiment of the present application discloses a kind of source image away from SID detection method and device, when needing to detect the SID of the radiation imaging with grid, the grid shadow curve of frequency distribution as caused by grid in available image, the main body distributing position of frequency distribution can be determined according to the grid shadow curve of frequency distribution, relational model corresponding with the grid technical parameter is inquired according to the main body distributing position, since the relational model can embody main body distributing position and the corresponding relationship of SID, therefore pass through the inquiry available SID corresponding with the main body distributing position of the relational model, so as to which the SID to be determined as to the SID of the current desired determination of the radiation equipmen.By grid, the corresponding relationship of main body distributing position and SID can accurately determine current SID in grid shadow curve of frequency distribution on image, without relying upon external measurement devices, reduce the cost of determining SID.

Description

A kind of source image distance detecting method and device
Technical field
This application involves data processing fields, more particularly to a kind of source image distance detecting method and device.
Background technique
In radiation line imaging, source image is away from (Source Image Distance, SID) for embodying radioactive source and image The distance between receiver.SID is a critically important parameter, such as the radiation agent of SID and radioactive source in radiation line imaging There is corresponding relationship, therefore some cases need to determine the required dose of radiation of imaging according to SID between amount.
In radiation imaging, in general SID is not fixed and invariable numerical value, can be adjusted with imaging demand It is whole.As the SID according to needed for imaging demand such as dose of radiation or irradiation with clearly defined objective Polaroid operation, need accurate Radioactive source is adjusted to the distance between panadaptor and meets the SID.For the numerical value of clear SID, traditional approach is by attached Add external measurement devices to measure SID.
However measuring device is additionally arranged in traditional approach needs on radiation imaging, will increase radiation line imaging and sets Standby whole design difficulty and cost.
Summary of the invention
In order to solve the above-mentioned technical problem, this application provides a kind of source image distance detecting methods, not against external equipment Under conditions of source image can be realized away from detection.
The embodiment of the present application discloses following technical solution:
In a first aspect, the embodiment of the present application provides a kind of source image away from SID detection method, applied to putting with grid Radiation imaging apparatus, when determining the SID of the radiation imaging, which comprises
Obtain the grid shadow curve of frequency distribution as caused by the grid in image;
Relational model is inquired according to the main body distributing position of frequency distribution in the grid shadow curve of frequency distribution, is obtained and institute The corresponding SID of main body distributing position is stated, the relational model corresponds to the technical parameter of the grid, for embodying frequency point The main body distributing position of cloth and the corresponding relationship of SID;
The SID is determined as to the SID of the radiation imaging.
Optionally, the relational model is obtained according to following manner:
At different SID, the multiple images obtained by the grid with the technical parameter are obtained;
According to the main body distributing position of frequency distribution in the grid shadow curve of frequency distribution of the multiple image, and To the multiple image SID determine include main body distributing position Yu SID corresponding relationship relational model.
Optionally, in acquired the multiple image, there is image different at same SID.
Optionally, the main body of frequency distribution is distributed in the grid shadow curve of frequency distribution according to the multiple image Position, and obtain the multiple image SID determine include main body distributing position Yu SID corresponding relationship relationship Model, comprising:
For the different images under the same SID, the grid shadow curve of frequency distribution of the different images is closed It and is a grid shadow curve of frequency distribution;
According to the grid shadow curve of frequency distribution after merging, and obtain the multiple image SID determine include The relational model of main body distributing position and SID corresponding relationship.
Optionally, the main body distributing position includes the corresponding frequency of maximum amplitude in the grid shadow curve of frequency distribution, Either determined according to the starting point and ending point of amplitude jump in the grid shadow curve of frequency distribution.
Second aspect, the embodiment of the present application provide a kind of source image away from SID detection device, applied to putting with grid Radiation imaging apparatus, when determining the SID of the radiation imaging, described device includes:
Module is obtained, for obtaining the grid shadow curve of frequency distribution as caused by the grid in image;
Enquiry module, for inquiring relationship according to the main body distributing position of frequency distribution in the grid shadow curve of frequency distribution Model obtains SID corresponding with the main body distributing position, and the relational model corresponds to the technical parameter of the grid, For embodying the main body distributing position of frequency distribution and the corresponding relationship of SID;
Determining module, for the SID to be determined as to the SID of the radiation imaging.
Optionally, the enquiry module includes:
Image acquiring unit, at different SID, acquisition to be obtained by the grid with the technical parameter Multiple images;
Relational model determination unit, for frequency distribution in the grid shadow curve of frequency distribution according to the multiple image Main body distributing position, and obtain the SID of the multiple image and determine to include that main body distributing position is corresponding with SID The relational model of relationship.
Optionally, in acquired the multiple image, there is image different at same SID.
Optionally, the determining relational model unit includes:
Grid shadow curve of frequency distribution merges subelement, the different images for being directed under the same SID, by described in not Grid shadow curve of frequency distribution with image merges into a grid shadow curve of frequency distribution;
Relational model determines subelement, for according to the grid shadow curve of frequency distribution after merging, and obtains the multiple The SID of image determine include main body distributing position Yu SID corresponding relationship relational model.
Optionally, the main body distributing position includes the corresponding frequency of maximum amplitude in the grid shadow curve of frequency distribution, Either determined according to the starting point and ending point of amplitude jump in the grid shadow curve of frequency distribution.
The third aspect, the embodiment of the present application provide a kind of processing equipment detected for source image away from SID, include storage Perhaps more than one program one of them or more than one program is stored in memory by device and one, and is configured It include for performing the following operation to execute the one or more programs by one or more than one processor Instruction:
Obtain the grid shadow curve of frequency distribution as caused by the grid in image;
Relational model is inquired according to the main body distributing position of frequency distribution in the grid shadow curve of frequency distribution, is obtained and institute The corresponding SID of main body distributing position is stated, the relational model corresponds to the technical parameter of the grid, for embodying frequency point The main body distributing position of cloth and the corresponding relationship of SID;
The SID is determined as to the SID of the radiation imaging.
Fourth aspect, the embodiment of the present application provide a kind of machine readable media, are stored thereon with instruction, when by one or When multiple processors execute so that device execute source image as described in one or more in first aspect away from detection method.
It, can it can be seen from above-mentioned technical proposal when it needs to be determined that having the SID of the radiation imaging of grid To obtain the grid shadow curve of frequency distribution as caused by grid in image, can be determined according to the grid shadow curve of frequency distribution The main body distributing position of frequency distribution out inquires relationship mould corresponding with the grid technical parameter according to the main body distributing position Type, since the relational model can embody main body distributing position and the corresponding relationship of SID, therefore can by inquiring the relational model To obtain SID corresponding with the main body distributing position, so as to which the SID is determined as the current desired determination of the radiation equipmen SID.By grid, the corresponding relationship of main body distributing position and SID can be in grid shadow curve of frequency distribution on image It accurately determines current SID, without relying upon external measurement devices, reduces the cost of determining SID.
Detailed description of the invention
In order to illustrate the technical solutions in the embodiments of the present application or in the prior art more clearly, to embodiment or will show below There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this Some embodiments of application without any creative labor, may be used also for those of ordinary skill in the art To obtain other drawings based on these drawings.
Fig. 1 is a kind of flow chart of SID detection method provided by the embodiments of the present application;
Fig. 2 is the grid shadow curve of frequency distribution figure under a certain SID;
Fig. 3 is a kind of flow chart of relational model method for building up provided by the embodiments of the present application;
Fig. 4 is to change a kind of grid shadow frequency corresponding to different images point that conditions of exposure obtains at a certain SID Cloth curve graph;
Fig. 5 is to change the grid shadow frequency corresponding to different images that another conditions of exposure obtains at a certain SID Scatter chart;
Fig. 6 is the grid shadow curve of frequency distribution figure corresponding to different SID;
Fig. 7 is the obtained grid shadow curve of frequency distribution figure corresponding to different SID after smoothed processing;
Fig. 8 is a kind of structure chart of the source image away from SID detection device provided by the embodiments of the present application;
Fig. 9 is a kind of structure drawing of device for source image away from the SID processing equipment detected applying for embodiment and providing.
Specific embodiment
With reference to the accompanying drawing, embodiments herein is described.
In radiation line imaging, for SID for embodying the distance between radioactive source and panadaptor, application is very extensive. For example, when design auto-exposure control scheme, due to dose of radiation and SID square between inversely, can will SID as control dose of radiation foundation, before exposure with obtain SID in exposure process, according to the SID to dose of radiation carry out Control.In addition, usually also needing to be set between radioactive source and panadaptor according to imaging demand in radioactive ray imaging field Distance be a certain sid value, SID is adjusted by adjusting the component in radiation imaging, so that radioactive source and shadow As the distance between receiver accurately attains the required SID, and then meet imaging demand.In conclusion radioactive ray at As in field, SID is a critically important parameter.
In the prior art, SID is generally determined using following three kinds of methods.
First method places master scale ruler, between manual measurement radioactive source and panadaptor at radioactive source Distance SID.
Laser range finder is arranged in second method, emits beam of laser, photoelectric cell to panadaptor from radioactive source The laser beam of panadaptor reflection is received, timer, which obtains the outgoing laser beam from radioactive source and receives to reflect to photoelectric cell, to swash The time expended during light beam calculates SID according to the time.
The third method, using potentiometer measurement SID, specifically, potentiometer is corrected before measurement SID, it is normal to survey When measuring SID, potentiometric value of feedback is obtained, according to the potentiometer value of feedback and correction data, obtains SID.
It can be found through observation, the prior art, which is generally required, obtains SID by the way that additional measuring device is arranged.For Radiation imaging, add obtain SID measuring device, will increase radiation imaging whole design difficulty and at This.
Due to being typically provided with the grid for filtering out scattered rays on radiation imaging, this application provides one kind SID determines method, without additional setting measuring device, grid shadow frequency caused by being directly based upon the grid on image point Cloth curve determines SID.
Specifically, obtaining the frequency of the grid shadow as caused by grid in image in the SID for determining radiation imaging Rate distribution curve determines the main body distributing position of frequency distribution according to the grid shadow curve of frequency distribution, is distributed position according to the main body Set corresponding with the grid technical parameter relational model of inquiry, due to the relational model be able to reflect main body distributing position and Corresponding relationship between SID, therefore, SID corresponding with current topic distributing position can be obtained by inquiring the relational model, by this SID is determined as the SID of the current desired determination of radiation imaging.
Determine that method determines the SID of radiation imaging using above-mentioned SID, without adding on radiation imaging Additional measuring device, directly according to grid on image main body distributing position and SID in grid shadow curve of frequency distribution Corresponding relationship, current SID can be accurately determined, reduce the cost of determining SID.
Embodiment one
It is the flow chart that a kind of SID provided in this embodiment determines method referring to Fig. 1.This method is applied to have filter line The radiation imaging of grid, in the SID for determining radiation imaging, this method comprises:
Step 101: obtaining the grid shadow curve of frequency distribution as caused by grid in image.
Grid in radiation imaging is generally positioned above panadaptor, radioactive source emitting radiation, in shadow As forming image on receiver, grid can filter out the scattered rays in ray, and occur grid in image Grid shadow.When obtaining the grid shadow curve of frequency distribution according to the grid shadow on image, sampling can be obtained on image Data element handles data element, to obtain the grid shadow curve of frequency distribution of the image.
It should be noted that the horizontal axis of the grid shadow curve of frequency distribution typicallys represent sample frequency, the longitudinal axis typicallys represent frequency Rate amplitude.
It is described below by a kind of optional method for obtaining grid shadow curve of frequency distribution provided by the embodiments of the present application:
In image, data element is extracted by row (or column) on the direction vertical with grid grid shadow, it is right respectively The data element of every row (or each column) carries out one-dimensional Fourier transform, obtains the frequency spectrum of the data element, on Jiang Gehang (or column) The frequency spectrum of data element carries out integrated treatment, and grid shadow curve of frequency distribution in present image can be obtained.
In order to retain the raw information in image as much as possible, the flat of the frequency spectrum of the data element on every row (or column) is sought Mean value integrates the spectral average of data element on all rows (or column), obtains grid shadow curve of frequency distribution, and only retain grid shadow The high frequency section of curve of frequency distribution.
Grid shadow curve of frequency distribution, In as caused by grid in image can certainly be obtained using other methods This is not defined the method for obtaining grid shadow curve of frequency distribution.
Step 102: relational model being inquired according to the main body distributing position of frequency distribution in grid shadow curve of frequency distribution, is obtained SID corresponding with main body distributing position.
Step 103: SID is determined as to the SID of radioactive ray query facility.
In the grid shadow curve of frequency distribution caused by grid, in certain sample frequency section, frequency amplitude can occur The jump of significant change, specifically, the Frequency and Amplitude of grid shadow curve of frequency distribution may be from frequency width in the sample frequency section Degree minimum value is increased to Frequency and Amplitude maximum value, then is reduced to Frequency and Amplitude minimum value, grid shadow frequency from Frequency and Amplitude maximum value The part of this amplitude jump can be considered as the main part of the grid shadow curve of frequency distribution in distribution curve, and can embody The frequency of institute's distributing position or frequency separation are the main body of frequency distribution in grid shadow curve of frequency distribution where the main part Distributing position.
At a fixed SID, different radiation imagings pass through radioactive source emitting radiation shape on panadaptor At image, if on the panadaptor of these imaging devices be provided with the identical grid of technical parameter, it is acquired at As image grid shadow curve of frequency distribution in, it is similar for capable of embodying the frequency of the distributing position of its main body.That is, In When having used the grid of same technique parameter, in the image that is obtained under same SID, the main body of grid shadow curve of frequency distribution Distributing position can embody the corresponding relationship between the SID.Moreover, inventors have found that using same technique parameter Under the premise of grid, with SID is changed, the main body distributing position in the grid shadow curve of frequency distribution of image also can be therewith Change, under normal circumstances, when SID increases, main body distributing position can reduce.
It is described below by a kind of optional main body distributing position for obtaining frequency separation form provided by the embodiments of the present application Method:
When obtaining the main body distributing position of frequency distribution in grid shadow curve of frequency distribution, curve of frequency distribution can be first determined The starting point and ending point of middle frequency amplitude jump, and then the sample frequency section among the starting point and ending point is determined as The main body distributing position of frequency distribution.Specifically, a threshold value can be set previously according to actual conditions, compare grid shadow frequency point Whether the absolute value that frequency amplitude jumps between each consecutive points in cloth curve is greater than the threshold value, if certain puts frequency between point adjacent thereto The absolute value of rate amplitude jump is greater than the threshold value, then illustrates that the point is the starting point or terminating point on grid shadow curve of frequency distribution. The starting point and ending point on grid shadow curve of frequency distribution is determined in the method, correspondingly, between the starting point and ending point Sample frequency section be frequency distribution main body distributing position.
The main body distributing position of frequency separation form can certainly be obtained using other methods, the embodiment of the present application is herein Any restriction is not done to the method for the main body distributing position for obtaining frequency separation form.
In addition, the body position of frequency distribution may be a certain frequency in grid shadow curve of frequency distribution in main part Point, if obtaining the Frequency point in main part marginal position as main body distributing position, due to being in the frequency of marginal position The external world that point may be subject to is affected, and therefore, obtains according to this in the main body distributing position of main part marginal position SID may be less accurate.Therefore, it when choosing in a certain Frequency point in main part as main body distributing position, does not select generally The Frequency point being affected in marginal position by the external world in main part is taken, the point of selection is generally to be influenced by the external world Lesser Frequency point.
It is described below by the main body distributing position of two kinds provided by the embodiments of the present application optional acquisition Frequency point forms Method:
First method obtains the corresponding Frequency point of maximum amplitude, as frequency distribution in grid shadow curve of frequency distribution Main body distributing position.When by external interference, the extraneous influence being subject at the maximum amplitude in curve of frequency distribution is relatively Small, correspondingly, the corresponding Frequency point of maximum amplitude is also smaller by being influenced.Therefore, by the corresponding Frequency point of the maximum amplitude As main body distributing position, accurate SID can be obtained according to the main body distributing position.
Second method obtains the intermediate frequency point between starting point and ending point, as in grid shadow curve of frequency distribution The main body distributing position of frequency distribution.The intermediate frequency point is the median of starting point and ending point.Due to the intermediate frequency point Middle part in main part, when by external interference, the intermediate frequency is relatively small by being influenced.It therefore, will be among this Frequency point can obtain accurate SID according to the main body distributing position as main body distributing position.
Certainly main body distributing position of other Frequency points in main part as frequency distribution, this implementation can also be obtained Example does not do any restriction to the position of the Frequency point as main body distributing position herein, is not also used as main body distributing position to acquisition The method of Frequency point do any restriction.
Due to there is certain corresponding pass in grid shadow curve of frequency distribution between the main body distributing position and SID of frequency distribution System, therefore, when it needs to be determined that obtaining the main body point of the frequency distribution of grid shadow curve of frequency distribution in image when current SID Behind cloth position, according to the main body distributing position got, inquiry is able to reflect corresponding relationship between main body distributing position and SID Relational model, can get the corresponding SID of the main body distributing position, and the SID can for it is current it needs to be determined that SID.
It is understood that when main body distributing position is frequency separation, the relational model inquired, to be able to reflect this The relational model of corresponding relationship between frequency separation and SID;When main body distributing position is Frequency point, the relationship mould inquired Type is the relational model that can reflect corresponding relationship between the Frequency point and SID.For example, when the main body distributing position chosen is most When the corresponding Frequency point of amplitude, inquiry is able to reflect the pass of corresponding relationship between the corresponding Frequency point of the maximum amplitude and SID It is model, the corresponding SID of acquisition Frequency point corresponding with current maximum amplitude;When selection main body distributing position be starting point and When intermediate frequency point between terminating point, inquiry is able to reflect the relational model of corresponding relationship between the intermediate frequency point and SID, Obtain SID corresponding with current intermediate frequency point.
It should be noted that above-mentioned relation model corresponds to the technical parameter of grid, the technical parameter of grid is different, Reflect that the relational model of corresponding relationship is different between the main body distributing position and SID of frequency distribution, wherein the technology of grid is joined Number may include: focal length, grid density and grid ratio of grid etc..
It should be noted that above-mentioned relation model can be function model, or for indicating main body distributing position The Table Model of one-to-one relationship between SID, does not do any restriction to the form of relational model herein.
The SID for corresponding to the main body distributing position of frequency distribution in grid shadow curve of frequency distribution is obtained according to the method described above, The SID is determined as to the SID of the radiation imaging.
In order to make it easy to understand, with reference to the accompanying drawing, determining that method is illustrated to SID provided in this embodiment:
The grid that focal length is 1.3m is chosen, the grid density of the grid is 40L/cm, and grid ratio is 10:1.It obtains a certain The grid shadow curve of frequency distribution as caused by grid in image under SID, as shown in Figure 2.
According to the method for above-mentioned acquisition starting point and ending point, starting point A and the end of the grid shadow curve of frequency distribution are obtained Stop B, correspondingly, the sample frequency section between starting point A and terminating point B are the main body distributing position of frequency distribution.It obtains The corresponding sample frequency of the midpoint C of starting point A and terminating point B inquires relational model corresponding with the technical parameter of the grid, What the relational model reflected is the corresponding relationship between midpoint sample frequency and SID, obtains sample frequency phase corresponding with midpoint C Corresponding SID.
SID provided in this embodiment determines method, additional without adding in the SID for determining radiation imaging Measuring device, directly according to grid, main body distributing position, inquiry reflection are somebody's turn to do in grid shadow curve of frequency distribution on image The relational model of main body distributing position and SID corresponding relationship determines current SID, reduces the cost of determining SID.
Above-described embodiment is in getting grid shadow curve of frequency distribution after the main body distributing position of frequency distribution, according to the master Body distributing position obtains the SID of radiation equipmen by inquiry relational model.When determining SID using the above method, need Know the relational model for being able to reflect corresponding relationship between main body distributing position and SID, below to the method for building up of relational model into Row is introduced in detail.
Embodiment two
It is a kind of flow chart of relational model method for building up provided in this embodiment referring to Fig. 3.This method comprises:
Step 301: at different SID, obtaining the multiple images obtained by the grid with technical parameter.
The SID for changing radiation equipmen, at different SID, obtain through grid filter out scattered ray treated at As image, the technical parameter of the grid is certain, correspondingly, is established using relational model method for building up provided in this embodiment Relational model, for the relational model for being directed to the grid with the technical parameter.While obtaining image, record is not With SID corresponding to image.
It should be noted that may be influenced by extraneous factor, directly when due to obtaining image at a certain SID It may inaccuracy according to the grid shadow curve of frequency distribution that the image obtains.It is built by the grid shadow curve of frequency distribution of inaccuracy The accuracy of vertical relational model can also have an impact, therefore in order to guarantee the grid shadow as caused by grid obtained according to image Curve of frequency distribution is more accurate, and different images can be repeatedly obtained at same SID.
Different images is obtained at same SID, can in the case where not changing conditions of exposure, repeatedly obtain at As image;Conditions of exposure can also be varied multiple times at the SID, corresponding to different conditions of exposures, repeatedly obtain image.
Step 302: according to the main body distributing position of frequency distribution in the grid shadow curve of frequency distribution of multiple images, with And obtain multiple images SID determine include main body distributing position Yu SID corresponding relationship relational model.
According to the multiple images for corresponding to different SID, the grid shadow frequency distribution for corresponding to each image is obtained Curve, specifically, can obtain according to the method for the acquisition grid shadow distribution curve provided in embodiment one and correspond to each imaging The grid shadow curve of frequency distribution of image, details are not described herein.
Due to can repeatedly obtain different images at same SID, obtaining the corresponding grid shadow frequency of each SID When rate distribution curve, it is directed to image different under same SID, needs the grid shadow frequency distribution by different images bent Line merges into a grid shadow curve of frequency distribution.
Specifically, can first handle the different images obtained at the SID, obtain with it is each under the SID The corresponding grid shadow curve of frequency distribution of a image is sought in each grid shadow curve of frequency distribution corresponding to same sample frequency The average value of frequency amplitude, the comprehensive frequency amplitude average value for corresponding to each sample frequency, can be obtained corresponding to the SID's Grid shadow curve of frequency distribution.Similarly, the grid shadow curve of frequency distribution corresponding to other SID can be obtained in this way.
At same SID, changes conditions of exposure, repeatedly obtain different images, and then according to multiple under the SID Image obtains curve of frequency distribution corresponding with each image.As shown in Figure 4 and Figure 5, be at same SID, It according to the curve of frequency distribution that the image under different exposure obtains, can be found through observation, change conditions of exposure, Frequency amplitude in each curve of frequency distribution can change, but in curve of frequency distribution frequency distribution main body distributing position It does not change, that is, corresponds to different conditions of exposures, the main body distributing position of frequency distribution does not change in curve of frequency distribution Become.It can be seen that when determining that method determines SID using the SID that embodiment one provides, it, will not be right when conditions of exposure changes It is had an impact according to the SID that the main body distributing position of frequency distribution in curve of frequency distribution is determined.
It, can be according to the determination provided in embodiment one after getting the grid shadow curve of frequency distribution corresponding to each SID The method of the main body distributing position of frequency distribution determines the main body distributing position of frequency distribution in each grid shadow curve of frequency distribution. Since the corresponding SID of each image is different, correspondingly, the main body distributing position of frequency distribution in each grid shadow curve of frequency distribution Corresponding SID is also different, according to the corresponding relationship between each main body distributing position and SID, establishes reflection SID and main body is distributed position The relational model of corresponding relationship between setting.
Specifically, available grid shadow curve of frequency distribution corresponding with a certain SID, and then obtain the grid shadow frequency distribution The main body distributing position of frequency distribution in curve records the corresponding relationship between the main body distributing position and SID.According to same Method obtains main body distributing position corresponding with each SID, and records the corresponding relationship between each main body distributing position and SID. Corresponding relationship between aforementioned body distributing position and SID is arranged, obtain being able to reflect main body distributing position and SID it Between corresponding relationship relational model.
Since the main body distributing position of frequency distribution can be frequency separation, or Frequency point correspondingly works as main body When distributing position is frequency separation, the corresponding relationship between the frequency separation and SID is recorded, according to the corresponding relationship, is established anti- Reflect the relational model of corresponding relationship between the frequency separation and SID;When main body distributing position is Frequency point, the Frequency point is recorded Corresponding relationship between SID establishes the relationship mould for reflecting corresponding relationship between the Frequency point and SID according to the corresponding relationship Type.
Below by taking main body distributing position is the corresponding Frequency point of maximum amplitude as an example, introduces one kind and mentioned by the embodiment of the present application The method of the relational model of confession:
The grid that focal length is 1.3m is chosen, at different SID, obtains the image obtained after grid processing Picture, and then obtain grid shadow curve of frequency distribution corresponding from the image under different SID.
As shown in fig. 6, for the use of focal length is 1.3m grid when, the grid shadow frequency corresponding to different SID of acquisition is divided Cloth curve.After grid shadow curve of frequency distribution corresponding to different SID is smoothed, as shown in Figure 7 correspond to is obtained The grid grid shadow curve of frequency distribution of different SID, in figure 6 and figure 7, the grid shadow frequency that curve 1, which is SID, to be obtained when being 0.8m Distribution curve, the grid shadow curve of frequency distribution that curve 2, which is SID, to be obtained when being 1.0m, the grid that curve 3, which is SID, to be obtained when being 1.3m Shadow curve of frequency distribution, the grid shadow curve of frequency distribution that curve 4, which is SID, to be obtained when being 1.5m, curve 5, which is SID, to be obtained when being 1.8m The grid shadow curve of frequency distribution taken.
Obtain the corresponding frequency of maximum frequency amplitude in each grid shadow curve of frequency distribution, wherein maximum frequency width in curve 1 Corresponding frequency is 0.9001 at value, and corresponding frequency is 0.8948 at maximum frequency amplitude in curve 2, maximum frequency in curve 3 Corresponding frequency is 0.8882 at rate amplitude, and corresponding frequency is 0.8856 at maximum frequency amplitude in curve 4, in curve 5 most Corresponding frequency is 0.8830 at big frequency amplitude.Arrange the corresponding pass at maximum frequency amplitude between corresponding frequency and SID System, available corresponding table as shown in Table 1.
Table 1
Data in table 1 are analyzed, it can be found that it is directed to the grid that the focal length is 1.3m, radiation equipmen Corresponding relationship between SID and frequency corresponding at maximum amplitude meets relational model (1) and relational model (2).
Y=ax2-bx+c (1)
Y=C*x-t (2)
It brings into the frequency in table 1 as x, SID in above-mentioned relation model (1), seeks normal in relational model as y Number system number a, b and c, to obtain the relational model for the corresponding relationship being able to reflect between SID and frequency corresponding at maximum amplitude (1).Similarly, it brings into the frequency in table 1 as x, SID in above-mentioned relation model (2), seeks normal in relational model as y Number system number C and t, to obtain the relational model for the corresponding relationship being able to reflect between SID and frequency corresponding at maximum amplitude (2)。
Further, it is also possible to choose the intermediate frequency point of starting point and ending point in main part as main body distributing position, Establish the relational model for being able to reflect corresponding relationship between the intermediate frequency point and SID.A kind of the embodiment of the present application is described below The optional relational model for establishing corresponding relationship between reflection SID and intermediate frequency point provided is introduced:
The grid that focal length is 1.3m is chosen, at different SID, obtains the image obtained after grid processing Picture, and then obtain grid shadow curve of frequency distribution corresponding from the image under different SID.
According to the method for the starting point and ending point of fixed grid shadow curve of frequency distribution really is provided in embodiment one, obtain each Starting point and ending point in grid shadow curve of frequency distribution calculates the median of starting point and ending point as midpoint frequency, note Record corresponds to the midpoint frequency of different SID.Specifically, midpoint frequency when SID is 0.8m be 0.9001, SID be 1.0m when The midpoint frequency that midpoint frequency when being 1.3m that midpoint frequency is 0.8968, SID is 0.8902, SID when being 1.5m is 0.8869, Midpoint frequency when SID is 1.8m is 0.8849.Arrange the corresponding relationship between above-mentioned SID and midpoint frequency, it is available such as Table is corresponded to shown in table 2.
Table 2
Data in table 2 are analyzed, it can be found that it is directed to the grid that the focal length is 1.3m, radiation equipmen Corresponding relationship between SID and frequency corresponding at maximum amplitude also corresponds to relational model (1) and relational model (2).
Y=ax2-bx+c (1)
Y=C*x-t (2)
It brings into the frequency in table 2 as x, SID in above-mentioned relation model (1), seeks normal in relational model as y Number system number a, b and c, to obtain the relational model for the corresponding relationship being able to reflect between SID and frequency corresponding at maximum amplitude (1).Similarly, it brings into the frequency in table 2 as x, SID in above-mentioned relation model (2), seeks normal in relational model as y Number system number C and t, to obtain the relational model for the corresponding relationship being able to reflect between SID and frequency corresponding at maximum amplitude (2)。
Relational model provided in this embodiment determines method, can obtain at different SID by with technical parameter Multiple images that grid obtains obtain the grid shadow curve of frequency distribution for corresponding to each image, determine the frequency The main body distributing position of frequency distribution in distribution curve, and then according to the corresponding relationship between each main body distributing position and SID, it builds Vertical relational model.Using the relational model, can in known grid shadow curve of frequency distribution frequency distribution main body distributing position In the case where, directly acquire SID.
Embodiment three
It is a kind of source image provided in this embodiment away from SID structure diagram of detection device, which is applied to have referring to Fig. 8 The radiation imaging of grid, when determining the SID of the radiation imaging, described device includes:
Module 801 is obtained, for obtaining the grid shadow curve of frequency distribution as caused by the grid in image;
Enquiry module 802, for being inquired according to the main body distributing position of frequency distribution in the grid shadow curve of frequency distribution Relational model, obtains SID corresponding with the main body distributing position, and the technology that the relational model corresponds to the grid is joined Number, for embodying the main body distributing position of frequency distribution and the corresponding relationship of SID;
Determining module 802, for the SID to be determined as to the SID of the radiation imaging.
Optionally, the enquiry module 802 includes:
Image acquiring unit, at different SID, acquisition to be obtained by the grid with the technical parameter Multiple images;
Relational model determination unit, for frequency distribution in the grid shadow curve of frequency distribution according to the multiple image Main body distributing position, and obtain the SID of the multiple image and determine to include that main body distributing position is corresponding with SID The relational model of relationship.
Optionally, in acquired the multiple image, there is image different at same SID.
Optionally, the determining relational model unit includes:
Grid shadow curve of frequency distribution merges subelement, the different images for being directed under the same SID, by described in not Grid shadow curve of frequency distribution with image merges into a grid shadow curve of frequency distribution;
Relational model determines subelement, for according to the grid shadow curve of frequency distribution after merging, and obtains the multiple The SID of image determine include main body distributing position Yu SID corresponding relationship relational model.
Optionally, the main body distributing position includes the corresponding frequency of maximum amplitude in the grid shadow curve of frequency distribution, Either determined according to the starting point and ending point of amplitude jump in the grid shadow curve of frequency distribution.
SID detection device provided in this embodiment, it is additional without adding in the SID for determining radiation imaging Measuring device, directly according to grid, main body distributing position, inquiry reflection are somebody's turn to do in grid shadow curve of frequency distribution on image The relational model of main body distributing position and SID corresponding relationship determines current SID, reduces the cost of determining SID.
The detection device from the angle of modular functionality entity to the source image in the embodiment of the present application away from SID carries out above Description, the embodiment of the present application also provides a kind of processing equipments detected for source image away from SID, below from the angle of hardware handles Processing equipment in the embodiment of the present application is described.
Fig. 9 is a kind of structure drawing of device for source image away from the SID processing equipment detected provided by the embodiments of the present application, institute Stating processing equipment 900 includes processor 902 and memory 901, in which:
The memory 901 is transferred to the processor 902 for storing program code, and by said program code;
The processor 902, for according to the instruction execution following methods in said program code:
Obtain the grid shadow curve of frequency distribution as caused by the grid in image;
Relational model is inquired according to the main body distributing position of frequency distribution in the grid shadow curve of frequency distribution, is obtained and institute The corresponding SID of main body distributing position is stated, the relational model corresponds to the technical parameter of the grid, for embodying frequency point The main body distributing position of cloth and the corresponding relationship of SID;
The SID is determined as to the SID of the radiation imaging.
A kind of non-transitorycomputer readable storage medium, when the instruction in the storage medium is executed by processor, So as to execute a kind of time detection method of the source image away from SID, which comprises
Obtain the grid shadow curve of frequency distribution as caused by the grid in image;
Relational model is inquired according to the main body distributing position of frequency distribution in the grid shadow curve of frequency distribution, is obtained and institute The corresponding SID of main body distributing position is stated, the relational model corresponds to the technical parameter of the grid, for embodying frequency point The main body distributing position of cloth and the corresponding relationship of SID;
The SID is determined as to the SID of the radiation imaging.
Those of ordinary skill in the art will appreciate that: realize that all or part of the steps of above method embodiment can pass through The relevant hardware of program instruction is completed, and foregoing routine can be stored in a computer readable storage medium, which exists When execution, step including the steps of the foregoing method embodiments is executed;And storage medium above-mentioned can be at least one in following media Kind: read-only memory (English: read-only memory, abbreviation: ROM), RAM, magnetic or disk etc. are various to be can store The medium of program code.
It should be noted that all the embodiments in this specification are described in a progressive manner, each embodiment it Between same and similar part may refer to each other, each embodiment focuses on the differences from other embodiments. For equipment and system embodiment, since it is substantially similar to the method embodiment, so describe fairly simple, The relevent part can refer to the partial explaination of embodiments of method.Equipment and system embodiment described above is only schematic , wherein unit may or may not be physically separated as illustrated by the separation member, it is shown as a unit Component may or may not be physical unit, it can and it is in one place, or may be distributed over multiple networks On unit.Some or all of the modules therein can be selected to achieve the purpose of the solution of this embodiment according to the actual needs. Those of ordinary skill in the art can understand and implement without creative efforts.
The above, only a kind of specific embodiment of the application, but the protection scope of the application is not limited thereto, Within the technical scope of the present application, any changes or substitutions that can be easily thought of by anyone skilled in the art, Should all it cover within the scope of protection of this application.Therefore, the protection scope of the application should be with scope of protection of the claims Subject to.

Claims (12)

1. a kind of source image is away from SID detection method, which is characterized in that applied to the radiation imaging with grid, true When the SID of the fixed radiation imaging, which comprises
Obtain the grid shadow curve of frequency distribution as caused by the grid in image;
Relational model is inquired according to the main body distributing position of frequency distribution in the grid shadow curve of frequency distribution, is obtained and the master The corresponding SID of body distributing position, the relational model corresponds to the technical parameter of the grid, for embodying frequency distribution The corresponding relationship of main body distributing position and SID;
The SID is determined as to the SID of the radiation imaging.
2. the method according to claim 1, wherein the relational model is obtained according to following manner:
At different SID, the multiple images obtained by the grid with the technical parameter are obtained;
According to the main body distributing position of frequency distribution in the grid shadow curve of frequency distribution of the multiple image, and obtain institute The SID for stating multiple images determine include main body distributing position Yu SID corresponding relationship relational model.
3. according to the method described in claim 2, it is characterized in that, having in acquired the multiple image same Different image under one SID.
4. according to the method described in claim 3, it is characterized in that, the grid shadow frequency according to the multiple image point The main body distributing position of frequency distribution in cloth curve, and obtain the SID of the multiple image and determine to include main body The relational model of distributing position and SID corresponding relationship, comprising:
For the different images under the same SID, the grid shadow curve of frequency distribution of the different images is merged into One grid shadow curve of frequency distribution;
According to the grid shadow curve of frequency distribution after merging, and obtain the multiple image SID determine include master The relational model of body distributing position and SID corresponding relationship.
5. the method according to claim 1, which is characterized in that the main body distributing position includes described The corresponding frequency of maximum amplitude in grid shadow curve of frequency distribution, or according to amplitude jump in the grid shadow curve of frequency distribution Starting point and ending point determine.
6. a kind of source image is away from SID detection device, which is characterized in that applied to the radiation imaging with grid, true When the SID of the fixed radiation imaging, described device includes:
Module is obtained, for obtaining the grid shadow curve of frequency distribution as caused by the grid in image;
Enquiry module, for inquiring relationship mould according to the main body distributing position of frequency distribution in the grid shadow curve of frequency distribution Type obtains SID corresponding with the main body distributing position, and the relational model corresponds to the technical parameter of the grid, uses In the corresponding relationship for the main body distributing position and SID for embodying frequency distribution;
Determining module, for the SID to be determined as to the SID of the radiation imaging.
7. device according to claim 6, which is characterized in that the enquiry module includes:
Image acquiring unit, at different SID, acquisition to be obtained more by the grid with the technical parameter A image;
Relational model determination unit, the master for frequency distribution in the grid shadow curve of frequency distribution according to the multiple image Body distributing position, and obtain the SID of the multiple image and determine to include main body distributing position and SID corresponding relationship Relational model.
8. device according to claim 7, which is characterized in that in acquired the multiple image, have same Different image under one SID.
9. device according to claim 8, which is characterized in that the relational model determination unit includes:
Grid shadow curve of frequency distribution merges subelement, for for the different images under the same SID, by it is described it is different at As the grid shadow curve of frequency distribution of image merges into a grid shadow curve of frequency distribution;
Relational model determines subelement, for according to the grid shadow curve of frequency distribution after merging, and obtains the multiple imaging The SID of image determine include main body distributing position Yu SID corresponding relationship relational model.
10. according to device described in claim 6 to 9 any one, which is characterized in that the main body distributing position includes described The corresponding frequency of maximum amplitude in grid shadow curve of frequency distribution, or according to amplitude jump in the grid shadow curve of frequency distribution Starting point and ending point determine.
11. a kind of processing equipment detected for source image away from SID, which is characterized in that include memory and one or one A above program, perhaps more than one program is stored in memory and is configured to by one or one for one of them It includes the instruction for performing the following operation that the above processor, which executes the one or more programs:
Obtain the grid shadow curve of frequency distribution as caused by grid in image;
Relational model is inquired according to the main body distributing position of frequency distribution in the grid shadow curve of frequency distribution, is obtained and the master The corresponding SID of body distributing position, the relational model corresponds to the technical parameter of the grid, for embodying frequency distribution The corresponding relationship of main body distributing position and SID;
The SID is determined as to the SID of radiation imaging.
12. a kind of machine readable media is stored thereon with instruction, when executed by one or more processors, so that processor Execute source image as described in one or more in claim 1 to 5 away from detection method.
CN201810005172.6A 2018-01-03 2018-01-03 A kind of source image distance detecting method and device Active CN107990853B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810005172.6A CN107990853B (en) 2018-01-03 2018-01-03 A kind of source image distance detecting method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810005172.6A CN107990853B (en) 2018-01-03 2018-01-03 A kind of source image distance detecting method and device

Publications (2)

Publication Number Publication Date
CN107990853A CN107990853A (en) 2018-05-04
CN107990853B true CN107990853B (en) 2019-11-22

Family

ID=62040826

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810005172.6A Active CN107990853B (en) 2018-01-03 2018-01-03 A kind of source image distance detecting method and device

Country Status (1)

Country Link
CN (1) CN107990853B (en)

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7806591B2 (en) * 2007-09-27 2010-10-05 Carestream Health, Inc. Alignment apparatus for imaging system using reflective element
JP5438493B2 (en) * 2009-12-22 2014-03-12 富士フイルム株式会社 Radiation imaging system and auxiliary device thereof
US8827554B2 (en) * 2010-04-13 2014-09-09 Carestream Health, Inc. Tube alignment for mobile radiography system
CN102478742B (en) * 2010-11-26 2014-03-05 深圳迈瑞生物医疗电子股份有限公司 Method and system for self-adaptive correction of exposure parameters in digital radiography
CN104068886B (en) * 2013-03-28 2018-02-06 上海联影医疗科技有限公司 Medical photography system and its method for imaging
CN104173066B (en) * 2013-07-31 2015-09-02 上海联影医疗科技有限公司 Detect the method for X-ray shooting system source image distance
CN105832359B (en) * 2016-03-22 2018-10-19 广州七喜医疗设备有限公司 A kind of method of adaptive x-ray grid grid shadow removal

Also Published As

Publication number Publication date
CN107990853A (en) 2018-05-04

Similar Documents

Publication Publication Date Title
US6915008B2 (en) Method and apparatus for multi-nodal, three-dimensional imaging
US8436912B2 (en) Range measurement using multiple coded apertures
EP1683546B1 (en) Relative calibration for dosimetric devices
US20110267485A1 (en) Range measurement using a coded aperture
US9183625B2 (en) Method of processing radiograph and apparatus for processing radiograph using the method in which hough transform and radon transform performed on image
EP2702556A1 (en) A method and apparatus for providing image data for constructing an image of a region of a target object
RU2013136547A (en) COMPUTER TOMOGRAPHY SYSTEM AND BOLUS MONITORING METHOD
CN105361955B (en) Medical imaging component, medical imaging apparatus and the method for acquiring patient motion
US20110115896A1 (en) High-speed and large-scale microscope imaging
CN100466976C (en) X-ray CT system and method with beam-hardening correction
Curry et al. Calibration of an array camera
CN110353714A (en) CT equipment method for detecting abnormality, device, equipment and storage medium
US20070090298A1 (en) Method To Determine The Depth-Of-Interaction Function For PET Detectors
EP2227705B1 (en) Method and apparatus for providing image data
US7433086B2 (en) Edge detection and correcting system and method
CN107990853B (en) A kind of source image distance detecting method and device
EP2953097B1 (en) Biais-free regularization for spectral phase-unwrapping in differential phase contrast imaging
US8442291B2 (en) Method and apparatus for establishing a position of a phantom in an image quality test
Calantropio et al. Evaluation of photogrammetric block orientation using quality descriptors from statistically filtered tie points
US11353599B2 (en) Method for tracking the position of an irradiating source
CN103955610B (en) Medical image computer-aided analysis method
CN105832359A (en) Self-adaptive X-ray grid shadow eliminating method
CN115839667B (en) Height measurement method, device, equipment and storage medium
CN109613462A (en) A kind of scaling method of CT imaging
CN109416276B (en) Vibration detection device, method, computer storage medium and pulse diagnosis instrument

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 110179 No. 177-1 Innovation Road, Hunnan District, Shenyang City, Liaoning Province

Applicant after: DongSoft Medical System Co., Ltd.

Address before: 110179 No. 177-1 Innovation Road, Hunnan District, Shenyang City, Liaoning Province

Applicant before: Dongruan Medical Systems Co., Ltd., Shenyang

GR01 Patent grant
GR01 Patent grant