CN105608282B - A kind of ray filtration design method - Google Patents
A kind of ray filtration design method Download PDFInfo
- Publication number
- CN105608282B CN105608282B CN201511031035.2A CN201511031035A CN105608282B CN 105608282 B CN105608282 B CN 105608282B CN 201511031035 A CN201511031035 A CN 201511031035A CN 105608282 B CN105608282 B CN 105608282B
- Authority
- CN
- China
- Prior art keywords
- filtering
- ray
- design method
- kev
- monte carlo
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
A kind of ray filtration design method of the present invention utilizes Monte Carlo software Fast simulation ray filtration, shortens the design cycle, and without processing a series of multiple filter plates, reduce design cost；In addition, directly analyzing filter effect using the Monte Carlo software, it is easier to operate, and improve the accuracy assessed filter effect.
Description
Technical field
The present invention relates to ray filtration pieces, and in particular to a kind of ray filtration design method.
Background technology
The X-ray generated applied to the x-ray source of industry and medical field has a certain amount of Energy distribution, that is, has more
Color.In medical field, low energy ray is small to the contribution of imaging, and will increase the absorption metering of patient；And in industrial circle,
The pleochroism of ray will produce hardening effect, reduce image quality.And ray filtration piece, can under the premise of not influencing imaging,
The content for reducing low energy ray in X-ray, reaches filter effect.
Currently, existing filter plate uses experimental method in design, include the following steps：
A) filtering material is selected, and is processed into a series of filter plates with different thickness；
B) filter effect of different filter plates in the series filter plate is measured, that is, measures ray and passes through agent after these filter plates
The case where amount and power spectrum；
C) filter effect is assessed, selects suitable filter plate.
This filtering design method need to use dedicated instrument when measuring ray energy spectrum, and this instrument cost is higher, and
It is not easy to operate；Meanwhile needing to process a series of multiple filter plates, processing, measurement process are complicated, and the process-cycle is long, causes
Design cycle is long；In addition, this design method is not accurate enough to the assessment of filter effect.
Invention content
In view of this, it is an object of the invention to optimize the design method of ray filtration piece, a kind of ray filtration piece is provided
Design method, with the accuracy for shortening the design cycle, reducing design cost, reducing operation difficulty, improve filter effect assessment.
In order to achieve the above objectives, the present invention provides the following technical solutions：
A kind of ray filtration design method, includes the following steps：
1) Monte Carlo software is utilized to calculate radiographic source in the roentgen dose X D for measuring position1With spectral distribution Ψ1(E),
Middle E is ray energy；
2) Monte Carlo software is utilized to set filtering condition；
3) Monte Carlo software setting filtering material and thickness parameter are utilized；
4) model is established, and filtered roentgen dose X D is calculated using Monte Carlo software2With spectral distribution Ψ2(E),
Middle E is ray energy；
5) utilize Monte Carlo software analytical procedure 4) in obtained roentgen dose X D2With spectral distribution Ψ2(E) whether meet
The filtering condition set in step 2), if satisfied, filter plate design is completed, if not satisfied, adjustment filtering sheet material and thickness ginseng
Number repeats step 3), 4), 5), until meeting filtering condition.
Further, after the completion of step 5), the result calculated according to step 5) simulation processes filter plate, and carries out practical survey
Amount is not required to modify to filtering chip architecture if actual measured results meet filtering requirements；It is adjusted if being unsatisfactory for filtering requirements
Entire filter wave plate thickness, until practical measurement meets filtering requirements.
Further, the filtering condition is filtering completely, that is, ray energy is less than a after filtering1The component of keV is 0, filtering
The ratio between roentgen dose X and filtering predose meet D afterwards2/D1>η, wherein a1KeV is the filtered energy boundary value of setting, and η is setting
D2/D1Ratio boundary.
Further, the filtering condition is that part filters, and ray energy is in a after filtering1KeV~a2Component within the scope of keV
The percentage for accounting for total amount is less than b%, and the ratio between roentgen dose X and filtering predose meet D after filtering2/D1>η, wherein a1KeV~
a2KeV be filtering energy limit range, b% be set filtering after ray energy in a1KeV~a2Accounting within the scope of keV
Boundary value, η are the D of setting2/D1Ratio boundary.
Further, the Monte Carlo software is the software write based on Monte Carlo method.
Further, the Monte Carlo software is Geant4 or EGS.
Further, the filtering material is metal material or alloy material.
Further, the metal material is copper or aluminium or stainless steel.
Further, it in step 5) when adjustment filter plate thickness, under conditions of ensureing that the complete filtering of satisfaction or part filter, selects
Take the ratio between roentgen dose X and filtering predose D after filtering2/D1Thickness when maximum.
The beneficial effects of the present invention are：
A kind of ray filtration design method of the present invention utilizes Monte Carlo software Fast simulation ray filtration, shortens
Design cycle, and without processing a series of multiple filter plates, reduce design cost；In addition, utilizing the Monte Carlo software
Directly analyze filter effect, it is easier to operate, and improve the accuracy assessed filter effect.
Description of the drawings
In order to keep the purpose of the present invention, technical solution and advantageous effect clearer, the present invention provides following attached drawing and carries out
Explanation：
Fig. 1 is existing filtering design method；
Fig. 2 is a kind of step flow diagram of ray filtration design method of the present invention.
Specific implementation mode
Below in conjunction with attached drawing, the preferred embodiment of the present invention is described in detail.
A kind of ray filtration design method, includes the following steps：
1) Monte Carlo software is utilized to calculate radiographic source in the roentgen dose X D for measuring position1With spectral distribution Ψ1(E),
Middle E is ray energy；
2) Monte Carlo software is utilized to set filtering condition；
3) Monte Carlo software setting filtering material and thickness parameter are utilized；
4) model is established, and filtered roentgen dose X D is calculated using Monte Carlo software2With spectral distribution Ψ2(E),
Middle E is ray energy；
5) utilize Monte Carlo software analytical procedure 4) in obtained roentgen dose X D2With spectral distribution Ψ2(E) whether meet
The filtering condition set in step 2), if satisfied, filter plate design is completed, if not satisfied, adjustment filtering sheet material and thickness ginseng
Number repeats step 3), 4), 5), until meeting filtering condition.
In the present embodiment, the Monte Carlo software is the software write based on Monte Carlo method.Specifically, the illiteracy is special
Caro software is Geant4 or EGS.
In the present embodiment, the filtering material is metal material or alloy material.Specifically, the metal material be copper or
Aluminium or stainless steel.The metal material or alloy material have preferable filter effect.
In the present embodiment, the filtering condition is filtering completely, that is, ray energy is less than a after filtering1The component of keV is 0,
The ratio between roentgen dose X and filtering predose meet D after filtering2/D1>η.Certainly, the filtering condition is alternatively part filtering, filtering
Ray energy is in a afterwards1KeV~a2The percentage that component within the scope of keV accounts for total amount is less than b%, roentgen dose X and filtering after filtering
The ratio between predose meets D2/D1>η.In specific implementation process, filtering condition can be set according to actual demand.Wherein,
a1KeV is the filtered energy boundary value of setting, a1KeV~a2KeV is the energy limit range of filtering, and ray energy exists after filtering
a1KeV~a2The percentage that component within the scope of keV accounts for total amount is less than b%, and η is the D of setting2/D1Ratio boundary.
A kind of ray filtration design method of the present invention utilizes Monte Carlo software Fast simulation ray filtration, shortens
Design cycle, and without processing a series of multiple filter plates, reduce design cost；In addition, utilizing the Monte Carlo software
Directly analyze filter effect, it is easier to operate, and improve the accuracy assessed filter effect.
As the improvement of said program, after the completion of step 5), the result calculated according to step 5) simulation processes filter plate, and
Practical measurement is carried out, is not required to modify to filtering chip architecture if actual measured results meet filtering requirements；If being unsatisfactory for filtering
Wave requires then to adjust filter plate thickness, until practical measurement meets filtering requirements.In the present embodiment, using practical measurement, into one
Step demonstrate,proves filter effect, and the filtering reliability of filter plate can be improved.
As the improvement of said program, when adjusting filter plate thickness in step 5), ensures to meet filtering completely or partly filter
Under conditions of wave, the ratio between roentgen dose X and filtering predose D after filtering are chosen2/D1Thickness when maximum, at this point, the thickness is most
The filter plate filter effect of good thickness, the thickness is best.
Finally illustrate, preferred embodiment above is merely illustrative of the technical solution of the present invention and unrestricted, although logical
It crosses above preferred embodiment the present invention is described in detail, however, those skilled in the art should understand that, can be
Various changes are made to it in form and in details, without departing from claims of the present invention limited range.
Claims (9)
1. a kind of ray filtration design method, it is characterised in that include the following steps：
1) Monte Carlo software is utilized to calculate radiographic source in the roentgen dose X D for measuring position1With spectral distribution Ψ1(E), wherein E is
Ray energy；
2) Monte Carlo software is utilized to set filtering condition；
3) Monte Carlo software setting filtering material and thickness parameter are utilized；
4) model is established, and filtered roentgen dose X D is calculated using Monte Carlo software2With spectral distribution Ψ2(E), wherein E
For ray energy；
5) utilize Monte Carlo software analytical procedure 4) in obtained roentgen dose X D2With spectral distribution Ψ2(E) whether meet step
2) filtering condition set in, if satisfied, filter plate design is completed, if not satisfied, adjustment filtering sheet material and thickness parameter, weight
Multiple step 3), 4), 5), until meeting filtering condition.
2. a kind of ray filtration design method according to claim 1, it is characterised in that：After the completion of step 5), according to
The result that step 5) simulation calculates processes filter plate, and carries out practical measurement, if actual measured results meet filtering requirements not
It need to modify to filtering chip architecture；Filter plate thickness is adjusted if being unsatisfactory for filtering requirements, until practical measure meets filtering
It is required that.
3. a kind of design method of ray filtration piece according to claim 1, it is characterised in that：The filtering condition has been
Full filtering, that is, ray energy is less than a after filtering1The component of keV is 0, and the ratio between roentgen dose X and filtering predose meet after filtering
D2/D1>η, wherein a1KeV is the filtered energy boundary value of setting, and η is the D of setting2/D1Ratio boundary.
4. a kind of ray filtration design method according to claim 1, it is characterised in that：The filtering condition is part
Filtering, ray energy is in a after filtering1KeV~a2The percentage that component within the scope of keV accounts for total amount is less than b%, ray after filtering
The ratio between dosage and filtering predose meet D2/D1>η, wherein a1KeV~a2KeV is the energy limit range of filtering, and b% is to set
Ray energy is in a after fixed filtering1KeV~a2Accounting boundary value within the scope of keV, η are the D of setting2/D1Ratio boundary.
5. a kind of ray filtration design method according to claim 1, it is characterised in that：It is characterized in that：The illiteracy
Special Caro software is the software write based on Monte Carlo method.
6. a kind of ray filtration design method according to claim 5, it is characterised in that：The Monte Carlo software is
Geant4 or EGS.
7. a kind of ray filtration design method according to claim 1, it is characterised in that：The filtering material is metal
Material.
8. a kind of ray filtration design method according to claim 7, it is characterised in that：The metal material be copper or
Aluminium or stainless steel.
9. a kind of ray filtration design method according to claim 1, it is characterised in that：Filter plate is adjusted in step 5)
When thickness, under conditions of ensureing that the complete filtering of satisfaction or part filter, the ratio between roentgen dose X and filtering predose after selection filters
D2/D1Thickness when maximum.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201511031035.2A CN105608282B (en) | 2015-12-31 | 2015-12-31 | A kind of ray filtration design method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201511031035.2A CN105608282B (en) | 2015-12-31 | 2015-12-31 | A kind of ray filtration design method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105608282A CN105608282A (en) | 2016-05-25 |
| CN105608282B true CN105608282B (en) | 2018-09-11 |
Family
ID=55988214
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201511031035.2A Active CN105608282B (en) | 2015-12-31 | 2015-12-31 | A kind of ray filtration design method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105608282B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115201227A (en) * | 2022-06-28 | 2022-10-18 | 北京山水云图科技有限公司 | Filter plate design method and device based on X-ray detection and computer readable storage medium |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103349555A (en) * | 2013-04-27 | 2013-10-16 | 中国人民解放军北京军区总医院 | X ray filtration device, X ray filtration system and mobile CT scanner |
| CN104318972A (en) * | 2014-10-29 | 2015-01-28 | 重庆真测科技股份有限公司 | Design method of X-ray flattening filter |
| CN204536229U (en) * | 2014-12-31 | 2015-08-05 | 同方威视技术股份有限公司 | A kind of filter and a kind of substance detecting device |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4611168B2 (en) * | 2005-10-07 | 2011-01-12 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Image reconstruction method and X-ray CT apparatus |
| EP1939147B2 (en) * | 2005-10-19 | 2020-06-24 | Nippon Electric Glass Co., Ltd. | Radiation shielding glass and method for manufacture thereof |
-
2015
- 2015-12-31 CN CN201511031035.2A patent/CN105608282B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103349555A (en) * | 2013-04-27 | 2013-10-16 | 中国人民解放军北京军区总医院 | X ray filtration device, X ray filtration system and mobile CT scanner |
| CN104318972A (en) * | 2014-10-29 | 2015-01-28 | 重庆真测科技股份有限公司 | Design method of X-ray flattening filter |
| CN204536229U (en) * | 2014-12-31 | 2015-08-05 | 同方威视技术股份有限公司 | A kind of filter and a kind of substance detecting device |
Non-Patent Citations (4)
| Title |
|---|
| X光管滤光片效应的蒙特卡罗仿真;卢艳 等;《核技术》;20121031;第35卷(第10期);第751-754页 * |
| X射线吸收光谱技术中的滤光片研制;张静 等;《核技术》;20040731;第27卷(第7期);第497-500页 * |
| 一种基于NP设计X射线均整片的方法;代志力 等;《全国射线数字成像与CT新技术》;20120731;第205-208页 * |
| 微型Ｘ射线管出射谱特征研究及Be窗厚度确定;谷懿 等;《光谱学与光谱分析》;20140131;第34卷(第1期);第252-256页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105608282A (en) | 2016-05-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9025724B2 (en) | Method for analyzing internal density of material by using X-ray computed tomography | |
| JP5975181B2 (en) | X-ray fluorescence analysis method and X-ray fluorescence analyzer | |
| CN104166962B (en) | A kind of cone beam CT scatter correction method using scattering kernel method | |
| Østerås et al. | Average glandular dose in paired digital mammography and digital breast tomosynthesis acquisitions in a population based screening program: effects of measuring breast density, air kerma and beam quality | |
| CN105528288A (en) | Software testing method and apparatus | |
| CN105608282B (en) | A kind of ray filtration design method | |
| US20230110386A1 (en) | System and method to compute a pixel sensitivity map of an imaging device | |
| Gambaccini et al. | Measurement of focal spot size in a 5.5 MeV linac | |
| Hanson et al. | Clinical implementation and rapid commissioning of an EPID based in-vivo dosimetry system | |
| CN110782952A (en) | Quantitative Monte Carlo simulation method and system for absorbed dose and storage medium | |
| CN109173082A (en) | A kind of design method of neutron apparatus for shaping | |
| CN102621576B (en) | Mammary dose distribution detection method and device for special CT (computed tomography) equipment for mammary gland | |
| CN204536229U (en) | A kind of filter and a kind of substance detecting device | |
| Zboray et al. | High-frame rate imaging of two-phase flow in a thin rectangular channel using fast neutrons | |
| CN104677922B (en) | A kind of hot γ irradiation collaboration aging testing method of cable | |
| CN104502947B (en) | The quick apparatus and method for obtaining flushing-free checking film dosimetry response curve | |
| Shu-Quan et al. | Corrections on energy spectrum and scatterings for fast neutron radiography at NECTAR facility | |
| CN106501293A (en) | A kind of method that XRF power spectrum background is deducted based on Wavelet Transform | |
| Dias et al. | Measurements of k0 and Q0 values for 64Zn (n, γ) 65Zn and 68Zn (n, γ) 69mZn reactions with covariance analysis | |
| KR102152639B1 (en) | Quantitative analytic method for lead or arsenic | |
| Hopkins et al. | Analytical corrections for beam-hardening and object scatter in volumetric computed tomography systems | |
| CN104318972A (en) | Design method of X-ray flattening filter | |
| Yu et al. | Heel effect adaptive flat field correction of digital x‐ray detectors | |
| CN106373100B (en) | A kind of hardening artifact bearing calibration of CT image | |
| JP3291253B2 (en) | X-ray fluorescence analyzer |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| PE01 | Entry into force of the registration of the contract for pledge of patent right | ||
| PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: A Design Method for Ray Filters Effective date of registration: 20230605 Granted publication date: 20180911 Pledgee: Chongqing Bank Co.,Ltd. Major Branch Pledgor: CHONGQING ZHENCE SCIENCE AND TECHNOLOGY Co.,Ltd. Registration number: Y2023500000039 |