CN107990853A - A kind of source image distance detecting method and device - Google Patents
A kind of source image distance detecting method and device Download PDFInfo
- Publication number
- CN107990853A CN107990853A CN201810005172.6A CN201810005172A CN107990853A CN 107990853 A CN107990853 A CN 107990853A CN 201810005172 A CN201810005172 A CN 201810005172A CN 107990853 A CN107990853 A CN 107990853A
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims description 54
- 230000005855 radiation Effects 0.000 claims abstract description 53
- 238000003384 imaging method Methods 0.000 claims abstract description 50
- 238000001514 detection method Methods 0.000 claims abstract description 13
- 238000012545 processing Methods 0.000 claims description 12
- 239000004744 fabric Substances 0.000 claims description 10
- 238000005259 measurement Methods 0.000 abstract description 6
- 230000002285 radioactive effect Effects 0.000 description 15
- 238000000926 separation method Methods 0.000 description 10
- 238000013461 design Methods 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B15/00—Measuring arrangements characterised by the use of electromagnetic waves or particle radiation, e.g. by the use of microwaves, X-rays, gamma rays or electrons
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating 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/02—Investigating 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/04—Investigating 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T7/00—Details of radiation-measuring instruments
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Chemical & Material Sciences (AREA)
- Electromagnetism (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Measurement Of Radiation (AREA)
Abstract
The embodiment of the present application discloses a kind of source image away from SID detection methods 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 can be obtained in 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 about according to the main body distributing position, since the relational model can embody main body distributing position and the correspondence of SID, therefore SID corresponding with the main body distributing position can be obtained by inquiring about the relational model, so as to which the SID is determined as the current desired definite SID of the radiation equipmen.By grid, the correspondence 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 definite SID.
Description
Technical field
This application involves data processing field, more particularly to a kind of source image distance detecting method and device.
Background technology
In line imaging is radiated, source image is used to embody radioactive source and image away from (Source Image Distance, SID)
The distance between receiver.SID is a critically important parameter in line imaging is radiated, such as the radiation agent of SID and radioactive source
There is correspondence, therefore certain situation needs to determine to be imaged required dose of radiation according to SID between amount.
In radiation imaging, in general SID is not changeless 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 the with clearly defined objective Polaroid operation of irradiation, it is necessary to accurate
Radioactive source is adjusted to the distance between panadaptor and meets the SID.For the numerical value of clear and definite SID, traditional approach is by attached
Add external measurement devices to measure SID.
But traditional approach needs additionally to set measuring apparatus on radiation imaging, radiation line imaging can be increased and set
Standby global design difficulty and cost.
The content of the invention
In order to solve the above-mentioned technical problem, this application provides a kind of source image distance detecting method, not against external equipment
Under conditions of i.e. can be achieved source image 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 methods, applied to putting with grid
Radiation imaging apparatus, when determining the SID of the radiation imaging, the described method includes:
Obtain image in as caused by the grid grid shadow curve of frequency distribution;
Relational model is inquired about 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 correspondence of SID;
The SID is determined as to the SID of the radiation imaging.
Optionally, the relational model is obtained according in the following manner:
Under 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
SID to the multiple image determines the relational model for including main body distributing position and SID correspondences.
Optionally, in acquired the multiple image, there is image different under 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 SID of the multiple image and determine the relation for including main body distributing position and SID correspondences
Model, including:
For the different images under same SID, the grid shadow curve of frequency distribution of the different images is closed
And it is a grid shadow curve of frequency distribution;
According to the grid shadow curve of frequency distribution after merging, and obtain the SID of the multiple image and determine to include
The relational model of main body distributing position and SID correspondences.
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 provides a kind of source image away from SID detection devices, applied to putting with grid
Radiation imaging apparatus, when determining the SID of the radiation imaging, described device includes:
Acquisition module, for obtaining in image the grid shadow curve of frequency distribution as caused by the grid;
Enquiry module, for inquiring about relation 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 correspondence 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, under different SID, obtaining and being 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 multiple image SID determine to include main body distributing position it is corresponding with SID
The relational model of relation.
Optionally, in acquired the multiple image, there is image different under same SID.
Optionally, the definite relational model unit includes:
Grid shadow curve of frequency distribution merges subelement, for for the different images under same SID, will described in not
Grid shadow curve of frequency distribution with image merges into a grid shadow curve of frequency distribution;
Relational model determination subelement, for according to the grid shadow curve of frequency distribution after merging, and obtains the multiple
The SID of image determines the relational model for including main body distributing position and SID correspondences.
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
Device, and one either more than one program one of them or more than one program storage in memory, and be configured
It is used for what is operated below to be performed the one or more programs by one or more than one processor and be included
Instruction:
Obtain image in as caused by the grid grid shadow curve of frequency distribution;
Relational model is inquired about 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 correspondence 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 perform so that device perform source image as described in one or more in first aspect away from detection method.
It can be seen from above-mentioned technical proposal it needs to be determined that the radiation imaging with grid SID when, can
To obtain in image the grid shadow curve of frequency distribution as caused by grid, can be determined according to the grid shadow curve of frequency distribution
Go out the main body distributing position of frequency distribution, relation mould corresponding with the grid technical parameter is inquired about according to the main body distributing position
Type, since the relational model can embody main body distributing position and the correspondence of SID, therefore can by inquiring about the relational model
To obtain SID corresponding with the main body distributing position, so as to which the SID is determined as, the radiation equipmen is current desired to be determined
SID.By grid, the correspondence of main body distributing position and SID can be with grid shadow curve of frequency distribution on image
Accurately determine current SID, without relying upon external measurement devices, reduce the cost of definite SID.
Brief description of the drawings
In order to illustrate the technical solutions in the embodiments of the present application or in the prior art more clearly, below will be to embodiment or existing
There is attached drawing needed in technology description to be briefly described, it should be apparent that, drawings in the following description are only this
Some embodiments of application, for those of ordinary skill in the art, without having to pay creative labor, may be used also
To obtain other attached drawings according to these attached drawings.
Fig. 1 is a kind of flow chart of SID detection methods 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 under a certain SID, changes the grid shadow frequency point corresponding to different images that a kind of conditions of exposure obtains
Cloth curve map;
Fig. 5 is under a certain SID, changes the grid shadow frequency corresponding to different images that another conditions of exposure obtains
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 structure chart of a kind of source image away from SID detection devices provided by the embodiments of the present application;
Fig. 9 is that a kind of of application embodiment offer is used for structure drawing of device of the source image away from the SID processing equipments detected.
Embodiment
Below in conjunction with the accompanying drawings, embodiments herein is described.
In line imaging is radiated, SID is used to embody the distance between radioactive source and panadaptor, and application is quite varied.
For example, during 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, in radioactive ray imaging field, usually also need to, according to imaging demand, set between radioactive source and panadaptor
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 into
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.
Second method, sets laser range finder, launches beam of laser, photoelectric cell to panadaptor at radioactive source
The laser beam of panadaptor reflection is received, timer, which obtains the outgoing laser beam at radioactive source and receives to reflect to photoelectric cell, to swash
The time expended during light beam, SID is calculated 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.
By observing it can be found that the prior art is generally required obtains SID by setting extra measuring apparatus.For
Radiation imaging, add obtain SID measuring apparatus, can increase radiation imaging global design difficulty and into
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, sets measuring apparatus without extra, is directly based upon grid grid shadow frequency caused by image point
Cloth curve determines SID.
Specifically, determine radiation imaging SID when, obtain image in as caused by grid grid shadow frequency
Rate distribution curve, the main body distributing position of frequency distribution is determined according to the grid shadow curve of frequency distribution, is distributed position according to the main body
Put corresponding with the grid technical parameter relational model of inquiry, due to the relational model can reflect main body distributing position and
Correspondence between SID, therefore, SID corresponding with current topic distributing position can be obtained by inquiring about the relational model, by this
SID is determined as the current desired definite SID of radiation imaging.
Determine that method determines the SID of radiation imaging using above-mentioned SID, without being added on radiation imaging
Extra measuring apparatus, directly according to grid on image main body distributing position and SID in grid shadow curve of frequency distribution
Correspondence, you can accurately determine current SID, reduce the cost of definite 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, when determining the SID of radiation imaging, this method includes:
Step 101:Obtain image in as caused by grid grid shadow curve of frequency distribution.
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 transverse 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:
It is right respectively by row (or row) extraction data element on the direction vertical with grid grid shadow in image
Often the data element of row (or each column) carries out one-dimensional Fourier transform, obtains the frequency spectrum of the data element, on Jiang Gehang (or row)
The frequency spectrum of data element carries out integrated treatment, you can obtains grid shadow curve of frequency distribution in present image.
In order to retain the raw information in image as much as possible, ask for data element on every row (or row) frequency spectrum it is flat
Average, integrates the spectral average of data element on all rows (or row), obtains grid shadow curve of frequency distribution, and only retain grid shadow
The high frequency section of curve of frequency distribution.
Can certainly utilize other methods obtain image in as caused by grid grid shadow curve of frequency distribution,
This is not defined the method for obtaining grid shadow curve of frequency distribution.
Step 102:Relational model is inquired about 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 grid shadow curve of frequency distribution caused by grid, in certain sample frequency section, frequency amplitude can occur
The saltus step of significant change, specifically, in the sample frequency section, the Frequency and Amplitude of grid shadow curve of frequency distribution may be from frequency width
Degree minimum value is increased to Frequency and Amplitude maximum, then is reduced to Frequency and Amplitude minimum value, grid shadow frequency from Frequency and Amplitude maximum
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.
Under a fixed SID, different radiation imagings shape on panadaptor by radioactive source emitting radiation
Into image, if being provided with the identical grid of technical parameter on the panadaptor of these imaging devices, it is resulting into
As image grid shadow curve of frequency distribution in, the frequency that can embody the distributing position of its main body is similar.That is,
When having used the grid of constructed parameter, in the image obtained under same SID, the main body of grid shadow curve of frequency distribution
Distributing position can embody the correspondence between the SID.Moreover, inventor has found, constructed parameter is being used
On the premise of grid, with change SID, 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 of frequency amplitude saltus step is more than the threshold value between each consecutive points in cloth curve, if certain puts frequency between point adjacent thereto
The absolute value of rate amplitude jump is more than the threshold value, then illustrates the point for 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.
Other methods can certainly be used to obtain the main body distributing position of frequency separation form, the embodiment of the present application is herein
The method of main body distributing position to obtaining frequency separation form does not do any restriction.
In addition, the body position of frequency distribution can also be a certain frequency in grid shadow curve of frequency distribution in main part
Point, if the Frequency point in main part marginal position is obtained as main body distributing position, due to the frequency in marginal position
The ectocine that point may be subject to is larger, therefore, is obtained according to the main body distributing position for being in main part marginal position
SID may be less accurate.Therefore, choose be in main part in a certain Frequency point as main body distributing position when, do not select generally
Take in main part in marginal position by the larger Frequency point of ectocine, the point of selection is generally by ectocine
Less 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 being subject to external interference, the ectocine 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 be subject to 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 being subject to external interference, the intermediate frequency is relatively small be subject to being influenced.Therefore, by among this
Frequency point can obtain accurate SID as main body distributing position according to the 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 during current SID, obtaining the main body point of the frequency distribution of grid shadow curve of frequency distribution in image
Behind cloth position, according to the main body distributing position got, inquiry can reflect correspondence between main body distributing position and SID
Relational model, you can the corresponding SID of the main body distributing position is got, and the SID can be to be currently needed for definite SID.
It is understood that when main body distributing position is frequency separation, the relational model inquired about, for that can reflect this
The relational model of correspondence between frequency separation and SID;When main body distributing position is Frequency point, the relation mould inquired about
Type, is the relational model that can reflect correspondence between the Frequency point and SID.For example, when the main body distributing position chosen is most
During the corresponding Frequency point of amplitude, inquiry can reflect the pass of correspondence 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 for starting point and
During intermediate frequency point between terminating point, inquiry can reflect the relational model of correspondence 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,
The relational model of correspondence is different between the main body distributing position and SID of reflection frequency distribution, wherein, the technology ginseng of grid
Number can 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 representing 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, below in conjunction with the accompanying drawings, method, which is illustrated, to be determined 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.Obtain a certain
Under SID in image as caused by grid grid shadow curve of frequency distribution, 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.Obtain
The corresponding sample frequencys of midpoint C of starting point A and terminating point B, inquire about relational model corresponding with the technical parameter of the grid,
What the relational model reflected is the correspondence between midpoint sample frequency and SID, obtains sample frequency phase corresponding with midpoint C
Corresponding SID.
SID provided in this embodiment determines method, extra without adding when determining the SID of radiation imaging
Measuring apparatus, 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 correspondences, determines current SID, reduces the cost of definite SID.
Above-described embodiment is after the main body distributing position of frequency distribution in getting grid shadow curve of frequency distribution, according to the master
Body distributing position, the SID of radiation equipmen is obtained by inquiring about relational model., it is necessary to when determining SID using the above method
The relational model of correspondence between main body distributing position and SID can be reflected by knowing, 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 includes:
Step 301:Under different SID, the multiple images obtained by the grid with technical parameter are obtained.
Change radiation equipmen SID, under different SID, obtain through grid filter out scattered ray processing after into
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, to be directed to the relational model of the grid with the technical parameter.While image is obtained, record is not
With the SID corresponding to image.
It should be noted that due to that during acquisition image, may be influenced under a certain SID by extraneous factor, directly
The grid shadow curve of frequency distribution obtained according to the image may be inaccurate.Built by inaccurate grid shadow curve of frequency distribution
The accuracy of vertical relational model can also have an impact, therefore in order to ensure 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 under same SID.
Different images is obtained under same SID, can in the case where not changing conditions of exposure, repeatedly obtain into
As image;Conditions of exposure can also be varied multiple times under 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 the SID of multiple images and determine the relational model for including main body distributing position and SID correspondences.
According to multiple images corresponding to different SID, the grid shadow frequency distribution corresponding to each image is obtained
Curve, specifically, each imaging can be corresponded to according to the method for the acquisition grid shadow distribution curve provided in embodiment one, acquisition
The grid shadow curve of frequency distribution of image, details are not described herein.
Due to can repeatedly obtain different images under same SID, obtaining the corresponding grid shadow frequencies of each SID
During rate distribution curve, image different under same SID is directed to, it is necessary to which the grid shadow frequency distribution of different images is bent
Line merges into a grid shadow curve of frequency distribution.
Specifically, can first handle the different images obtained under the SID, obtain with it is each under the SID
The corresponding grid shadow curve of frequency distribution of a image, is asked 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 corresponding to each sample frequency, you can obtain 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.
Under same SID, change 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 under same SID,
The curve of frequency distribution obtained according to the image under different exposure, by observation it can be found that changing 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
Do not change, i.e., corresponding 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 from the above that when determining that method determines SID using the SID that embodiment one provides, will not be right when conditions of exposure changes
Had an impact according to the SID that the main body distributing position of frequency distribution in curve of frequency distribution is determined.
After getting the grid shadow curve of frequency distribution corresponding to each SID, it can be determined according to what is provided in embodiment one
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 correspondence between each main body distributing position and SID, establishes reflection SID and main body distribution position
The relational model of correspondence between putting.
Specifically, grid shadow curve of frequency distribution corresponding with a certain SID can be obtained, and then obtain the grid shadow frequency distribution
The main body distributing position of frequency distribution in curve, records the correspondence between the main body distributing position and SID.According to same
Method, obtains main body distributing position corresponding with each SID, and records the correspondence between each main body distributing position and SID.
Correspondence between aforementioned body distributing position and SID is arranged, obtain reflecting main body distributing position and SID it
Between correspondence 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 correspondence between the frequency separation and SID is recorded, according to the correspondence, is established anti-
Reflect the relational model of correspondence between the frequency separation and SID;When main body distributing position is Frequency point, the Frequency point is recorded
Correspondence between SID, according to the correspondence, establishes the relation mould for reflecting correspondence between the Frequency point and SID
Type.
Below by taking main body distributing position is the corresponding Frequency point of maximum amplitude as an example, introduce one kind and carried 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, under different SID, obtains the image obtained after grid processing
Picture, and then obtain and the corresponding grid shadow curve of frequency distribution of 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 point of acquisition
Cloth curve.After being smoothed corresponding to the grid shadow curve of frequency distribution of different SID, 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, curve 1 are the grid shadow frequency obtained when SID is 0.8m
Distribution curve, curve 2 are the grid shadow curve of frequency distribution obtained when SID is 1.0m, and curve 3 is the grid obtained when SID is 1.3m
Shadow curve of frequency distribution, curve 4 are the grid shadow curve of frequency distribution obtained when SID is 1.5m, and curve 5 is that SID is obtained when being 1.8m
The grid shadow curve of frequency distribution taken.
The corresponding frequency of peak frequency amplitude in each grid shadow curve of frequency distribution is obtained, wherein, peak frequency width in curve 1
Corresponding frequency is 0.9001 at value, and corresponding frequency is 0.8948 at peak 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 peak frequency amplitude in curve 4, in curve 5 most
Corresponding frequency is 0.8830 at big frequency amplitude.Arrange the corresponding pass between corresponding frequency and SID at peak frequency amplitude
System, can obtain corresponding table as shown in Table 1.
Table 1
Data in table 1 are analyzed, it can be found that the grid that the focal length is 1.3m is directed to, radiation equipmen
Correspondence between SID frequencies corresponding with maximum amplitude, meets relational model (1) and relational model (2).
Y=ax2-bx+c (1)
Y=C*x-t (2)
Bring into the frequency in table 1 as x, SID in above-mentioned relation model (1), ask for normal in relational model as y
Number system number a, b and c, to obtain the relational model that can reflect the correspondence between SID frequencies corresponding with maximum amplitude
(1).Similarly, bring into the frequency in table 1 as x, SID in above-mentioned relation model (2), ask for normal in relational model as y
Number system number C and t, to obtain the relational model that can reflect the correspondence between SID frequencies corresponding with maximum amplitude
(2)。
Further, it is also possible to the intermediate frequency point of starting point and ending point in main part is chosen as main body distributing position,
The relational model of correspondence between the intermediate frequency point and SID can be reflected by establishing.A kind of the embodiment of the present application is described below
The optional relational model for establishing correspondence between reflection SID and intermediate frequency point provided is introduced:
The grid that focal length is 1.3m is chosen, under 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 that 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
Midpoint frequency of the record corresponding to different SID.Specifically, midpoint frequency when SID is 0.8m is 0.9001, SID when being 1.0m
Midpoint frequency is that midpoint frequency when 0.8968, SID is 1.3m is that midpoint frequency when 0.8902, SID is 1.5m is 0.8869,
Midpoint frequency when SID is 1.8m is 0.8849.Arrange the correspondence between above-mentioned SID and midpoint frequency, can obtain as
Corresponding table shown in table 2.
Table 2
Data in table 2 are analyzed, it can be found that the grid that the focal length is 1.3m is directed to, radiation equipmen
Correspondence between SID frequencies corresponding with maximum amplitude, also corresponds to relational model (1) and relational model (2).
Y=ax2-bx+c (1)
Y=C*x-t (2)
Bring into the frequency in table 2 as x, SID in above-mentioned relation model (1), ask for normal in relational model as y
Number system number a, b and c, to obtain the relational model that can reflect the correspondence between SID frequencies corresponding with maximum amplitude
(1).Similarly, bring into the frequency in table 2 as x, SID in above-mentioned relation model (2), ask for normal in relational model as y
Number system number C and t, to obtain the relational model that can reflect the correspondence between SID frequencies corresponding with maximum amplitude
(2)。
Relational model provided in this embodiment determines method, can be obtained under different SID by with technical parameter
Multiple images that grid obtains, obtain the grid shadow curve of frequency distribution corresponding to each image, determine the frequency
The main body distributing position of frequency distribution in distribution curve, and then according to the correspondence between each main body distributing position and SID, build
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 of, directly acquire SID.
Embodiment three
Referring to Fig. 8, be a kind of source image provided in this embodiment away from SID structure of the detecting device figures, which is applied to have
The radiation imaging of grid, when determining the SID of the radiation imaging, described device includes:
Acquisition module 801, for obtaining in image the grid shadow curve of frequency distribution as caused by the grid;
Enquiry module 802, for being inquired about 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 correspondence 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, under different SID, obtaining and being 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 multiple image SID determine to include main body distributing position it is corresponding with SID
The relational model of relation.
Optionally, in acquired the multiple image, there is image different under same SID.
Optionally, the definite relational model unit includes:
Grid shadow curve of frequency distribution merges subelement, for for the different images under same SID, will described in not
Grid shadow curve of frequency distribution with image merges into a grid shadow curve of frequency distribution;
Relational model determination subelement, for according to the grid shadow curve of frequency distribution after merging, and obtains the multiple
The SID of image determines the relational model for including main body distributing position and SID correspondences.
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 devices provided in this embodiment, it is extra without adding when determining the SID of radiation imaging
Measuring apparatus, 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 correspondences, determines current SID, reduces the cost of definite SID.
The angle of slave module functional entity carries out detection device of the source image in the embodiment of the present application away from SID above
Description, the embodiment of the present application additionally provides a kind of processing equipment 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 a kind of is used for structure drawing of device of the source image away from the SID processing equipments detected, institute to be provided by the embodiments of the present application
Stating processing equipment 900 includes processor 902 and memory 901, wherein:
The memory 901, the processor 902 is transferred to for store program codes, and by said program code;
The processor 902, following methods are performed for the instruction in said program code:
Obtain image in as caused by the grid grid shadow curve of frequency distribution;
Relational model is inquired about 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 correspondence 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 performed by processor,
So as to a kind of time detection method of the source image away from SID is performed, the described method includes:
Obtain image in as caused by the grid grid shadow curve of frequency distribution;
Relational model is inquired about 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 correspondence of SID;
The SID is determined as to the SID of the radiation imaging.
One of ordinary skill in the art will appreciate that:Realizing all or part of step of above method embodiment can pass through
The relevant hardware of programmed instruction is completed, and foregoing routine can be stored in a computer read/write memory medium, which exists
During execution, execution the step of including above method embodiment;And foregoing storage medium can be in following media at least one
Kind:Read-only storage (English:Read-only memory, abbreviation:ROM), RAM, magnetic disc or CD etc. are various to store
The medium of program code.
It should be noted that each embodiment in this specification is described by the way of progressive, each embodiment it
Between identical similar part mutually referring to what each embodiment stressed is the difference with other embodiment.
For equipment and system embodiment, since it is substantially similar to embodiment of the method, 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 may or may not be as the unit that separating component illustrates physically separate, shown as unit
Component may or may not be physical location, you can with positioned at a place, or can also be distributed to multiple networks
On unit.Some or all of module therein can be selected to realize the purpose of this embodiment scheme according to the actual needs.
Those of ordinary skill in the art are without creative efforts, you can to understand and implement.
The above, is only a kind of embodiment of the application, but the protection domain of the application is not limited thereto,
Any one skilled in the art is in the technical scope that the application discloses, the change or replacement that can readily occur in,
It should all cover within the protection domain of the application.Therefore, the protection domain of the application should be with scope of the claims
Subject to.
Claims (12)
1. a kind of source image is away from SID detection methods, it is characterised in that applied to the radiation imaging with grid, true
During the SID of the fixed radiation imaging, the described method includes:
Obtain image in as caused by the grid grid shadow curve of frequency distribution;
Relational model is inquired about 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 correspondence of main body distributing position and SID;
The SID is determined as to the SID of the radiation imaging.
2. according to the method described in claim 1, it is characterized in that, the relational model is obtained according in the following manner:
Under 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 determines the relational model for including main body distributing position and SID correspondences.
3. according to the method described in claim 2, it is characterized in that, in acquired the multiple image, have same
Different image under one SID.
4. the according to the method described in claim 3, it is characterized in that, grid shadow frequency point according to the multiple image
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 correspondences, including:
For the different images under 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 SID of the multiple image and determine to include master
The relational model of body distributing position and SID correspondences.
5. method according to any one of claims 1 to 4, it is characterised 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 devices, it is characterised in that applied to the radiation imaging with grid, true
During the SID of the fixed radiation imaging, described device includes:
Acquisition module, for obtaining in image the grid shadow curve of frequency distribution as caused by the grid;
Enquiry module, for inquiring about relation 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 main body distributing position and the correspondence of SID that embody frequency distribution;
Determining module, for the SID to be determined as to the SID of the radiation imaging.
7. device according to claim 6, it is characterised in that the enquiry module includes:
Image acquiring unit, under different SID, obtain obtained by the grid with the technical parameter it is more
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 correspondences
Relational model.
8. device according to claim 7, it is characterised in that in acquired the multiple image, have same
Different image under one SID.
9. device according to claim 8, it is characterised in that the definite relational model unit includes:
Grid shadow curve of frequency distribution merges subelement, for for the different images under same SID, by it is described it is different into
As the grid shadow curve of frequency distribution of image merges into a grid shadow curve of frequency distribution;
Relational model determination subelement, for according to the grid shadow curve of frequency distribution after merging, and obtains the multiple imaging
The SID of image determines the relational model for including main body distributing position and SID correspondences.
10. according to the device described in claim 6 to 9 any one, it is characterised 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, it is characterised in that include memory, and one or one
Program more than a, either more than one program storage in memory and is configured to by one or one for one of them
Above processor performs the one or more programs and includes the instruction for being used for being operated below:
Obtain image in as caused by the grid grid shadow curve of frequency distribution;
Relational model is inquired about 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 correspondence of main body distributing position and SID;
The SID is determined as to the SID of the radiation imaging.
12. a kind of machine readable media, is stored thereon with instruction, when executed by one or more processors so that device is held
Source image of the row as described in one or more in claim 1 to 5 away from detection method.
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 true CN107990853A (en) | 2018-05-04 |
CN107990853B 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) |
Citations (7)
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 |
CN102138803A (en) * | 2009-12-22 | 2011-08-03 | 富士胶片株式会社 | Radiation imaging system and assist apparatus for the same |
US20120134475A1 (en) * | 2010-11-26 | 2012-05-31 | Shenzhen Mindray Bio-Medical Electronics Co., Ltd. | Methods and systems for adaptively correcting exposure parameters during digital radiographic imaging |
CN102917646A (en) * | 2010-04-13 | 2013-02-06 | 卡尔斯特里姆保健公司 | Tube alignment for mobile radiography system |
CN104068886A (en) * | 2013-03-28 | 2014-10-01 | 上海联影医疗科技有限公司 | Medical photography system and photography method implemented by same |
CN104173066A (en) * | 2013-07-31 | 2014-12-03 | 上海联影医疗科技有限公司 | Method for detecting source image distance of X-ray shooting system |
CN105832359A (en) * | 2016-03-22 | 2016-08-10 | 广州七喜医疗设备有限公司 | Self-adaptive X-ray grid shadow eliminating method |
-
2018
- 2018-01-03 CN CN201810005172.6A patent/CN107990853B/en active Active
Patent Citations (7)
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 |
CN102138803A (en) * | 2009-12-22 | 2011-08-03 | 富士胶片株式会社 | Radiation imaging system and assist apparatus for the same |
CN102917646A (en) * | 2010-04-13 | 2013-02-06 | 卡尔斯特里姆保健公司 | Tube alignment for mobile radiography system |
US20120134475A1 (en) * | 2010-11-26 | 2012-05-31 | Shenzhen Mindray Bio-Medical Electronics Co., Ltd. | Methods and systems for adaptively correcting exposure parameters during digital radiographic imaging |
CN104068886A (en) * | 2013-03-28 | 2014-10-01 | 上海联影医疗科技有限公司 | Medical photography system and photography method implemented by same |
CN104173066A (en) * | 2013-07-31 | 2014-12-03 | 上海联影医疗科技有限公司 | Method for detecting source image distance of X-ray shooting system |
CN105832359A (en) * | 2016-03-22 | 2016-08-10 | 广州七喜医疗设备有限公司 | Self-adaptive X-ray grid shadow eliminating method |
Also Published As
Publication number | Publication date |
---|---|
CN107990853B (en) | 2019-11-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4366362B2 (en) | Dose response calibration system and dose response calibration method | |
EP2994776B1 (en) | Apparatus and method for the evaluation of gamma radiation events | |
US20020167726A1 (en) | Method and apparatus for multi-nodal, three-dimensional imaging | |
AU2013260650B2 (en) | Rotational phase unwrapping | |
KR20180048680A (en) | Multi-Method Imaging System and Method | |
CN103596502A (en) | Adaptive calibration for tomographic imaging systems | |
US10779776B2 (en) | Apparatus for X-ray imaging an object | |
CN105388478A (en) | Method and arrangement for detecting acoustic and optical information as well as a corresponding computer program and a corresponding computer-readable storage medium | |
CN105361955B (en) | Medical imaging assembly, medical imaging device, and method for capturing patient motion | |
JP2017027357A (en) | Image processing unit, imaging system and image processing method | |
CN108888284A (en) | Image adjusting method, device and equipment, storage medium | |
EP1629438A1 (en) | Automatic optimal view determination for cardiac acquisitions | |
US20130163843A1 (en) | Method of processing radiograph and apparatus for processing radiograph using the method | |
KR20190123865A (en) | Calibration method of x-ray apparatus and calibration apparatus for the same | |
US7433086B2 (en) | Edge detection and correcting system and method | |
CN102119341B (en) | Method of meshing and calculating a volume in an ultrasound imaging system | |
EP2823465A1 (en) | Stereo x-ray tube based suppression of outside body high contrast objects | |
CN107990853B (en) | A kind of source image distance detecting method and device | |
US20220253508A1 (en) | Consecutive approximation calculation method, consecutive approximation calculation device, and program | |
Calantropio et al. | Evaluation of photogrammetric block orientation using quality descriptors from statistically filtered tie points | |
CN113362404A (en) | Scatter correction method, device and storage medium for computer tomography | |
CN105832359B (en) | A kind of method of adaptive x-ray grid grid shadow removal | |
CN109724690A (en) | Automatic reconstruction method of focal spot for far-field measurement by schlieren method based on spectral angle mapping | |
CN109613593B (en) | Gamma photon action position sampling method, gamma photon action position sampling device, electronic equipment and storage medium | |
CN111588997B (en) | Method, device and equipment for determining scattering quantity and storage medium |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB02 | Change of applicant information |
Address after: 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 |
|
CB02 | Change of applicant information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |