CN103995010A - X-射线成像空间分布不均匀性的校正方法 - Google Patents
X-射线成像空间分布不均匀性的校正方法 Download PDFInfo
- Publication number
- CN103995010A CN103995010A CN201410036628.7A CN201410036628A CN103995010A CN 103995010 A CN103995010 A CN 103995010A CN 201410036628 A CN201410036628 A CN 201410036628A CN 103995010 A CN103995010 A CN 103995010A
- Authority
- CN
- China
- Prior art keywords
- flat panel
- image
- panel detector
- theta
- bearing calibration
- 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
Landscapes
- Measurement Of Radiation (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开了X-射线成像空间分布不均匀性的校正方法,首先使平板探测器在位置z1输出一个均匀的图像,然后使平板探测器在位置z2收到一帧强度不均匀分布图像,通过对这一帧强度不均匀分布图像用二次多项式函数进行拟合得到图像校正系数,这样将平板探测器移动到任何一个待矫正的位置z,通过二次多项式函数及图像校正系数即可得到位置z上的均匀图像。因而,就可以将X-射线所成的像的不均匀性进行修正,从而达到均匀成像的效果。
Description
技术领域
本发明涉及X-射线成像方法。
背景技术
通常我们的X-射线成像方法是将X-射线源产生的射线直接照射到被探测的物体、然后透射在探测器表面上而生成的像。如果不考虑量子噪声,理想状态是当一个X-射线源直接照射到远离射线源的平板探测器上时,X-射线强度在平板表面的分布是近似均匀的,因此平板探测器接受到该X-射线所成的像是从一个简单的点光源照射所成的像。但是,实际的X-射线成像,是将射线源放在距离平板探测器相对比较近的范围内,以保证有足够的射线强度穿透物体进行投影成像。这时X-射线投影在平板探测器表面的强度分布就会不是均匀了。这种不均匀性由两方面造成的:一方面是在近距离情况下,由于平板探测器表面接收单元距离射线源不同所形成(射线强度与该距离平方成反比);另一方面是由于X-射线强度本身随着射线离开焦点的角度变化而变化的,这称为Heel效应。Heel效应是由于靠近阳极这一侧的部分射线会被阳极靶材料所吸收,因此靠近阳极这一侧的射线强度要低于靠近阴极一侧的射线强度。例如,对于一个43cmx43cm的平板探测器放置于1米处(如图1所示)。阴极端的X射线相对强度是105%,而阳极端的X射线相对强度只有70%。两边的射线强度差35%。这些射线强度的不均匀性,对实际应用带来困扰。例如,X-射线所成的像就会出现图像亮暗不均匀。
发明内容
本发明所要解决的技术问题就是提供一种X-射线成像空间分布不均匀性的 校正方法,将X-射线所成的像的不均匀性进行修正,从而达到均匀成像的效果。
为解决上述技术问题,本发明采用如下技术方案:X-射线成像空间分布不均匀性的校正方法,至少包括如下步骤:
将平板探测器放在位置z1并正对X-射线源,调节平板探测器每一个像素单元的灵敏度,使得输出的信号形成一个均匀的图像;
将平板探测器移动到位置z2且仍然正对X-射线源,这时平板探测器就会收到一帧强度不均匀分布图像;
将位置z2上这一帧强度不均匀分布图像沿X、Y轴采集信号,并将该信号用二次多项式函数进行拟合得到图像校正系数;
将平板探测器移动到任何一个待矫正的位置z,通过二次多项式函数及图像校正系数得到位置z上的均匀图像。
优选的,位置z1为平板探测器距离X-射线源足够远的位置zmax>>L,这里L是平板的尺度,位置z2为平板探测器所需成像最近的距离位置zmin。
进一步的,位置z上的均匀图像
这里t(x,y)-I(x,y)/I0是相对强度分布,t(0,0)-1,A,B,C,D是二项式拟合系数,θx和θy是该点(x,y)与X‐射线源连线与Z轴分别在Y‐Z面和X‐Z面上投影的夹角,z满足:zmin≤z≤zmax。
进一步的,沿着X轴(θy=0)或Y轴(θy=0)上的强度分布分别是:
依据上述两个公式求出二项式拟合系数A,B,C,D。
优选的,用二次多项式函数进行拟合,求出a0,b0,c0,和d0,
由于a0=A(1/z1-1/z2)、 c0=C(1/z1-1/z2)和 因此可以求出二项式拟合系数A,B,C,D:A=a0×(1/z1-1/z2)-1、 C=c0×(1/z1-1/z2)-1和
进一步的,将平板探测器移到任何一个待矫正的位置z,这时可以得到一组系数:a=A(1/z1-1/z)、 c=C(1/z1-1/z)和 然后将在z处的采集的图像信号乘以以下函数:
Δt(x,y)=(1|ax|bx2)×(1|cy|dy2)
这样,就可以得到矫正之后的图像。
本发明依据Heel效应的原理,首先使平板探测器在位置z1输出一个均匀的图像,然后使平板探测器在位置z2收到一帧强度不均匀分布图像,通过对这一 帧强度不均匀分布图像用二次多项式函数进行拟合得到图像校正系数,这样将平板探测器移动到任何一个待矫正的位置z,通过二次多项式函数及图像校正系数即可得到位置z上的均匀图像。因而,就可以将X-射线所成的像的不均匀性进行修正,从而达到均匀成像的效果。
附图说明
下面结合附图和具体实施方式对本发明作进一步描述:
图1为Heel效应的实例图;
图2为点光源强度与距光源距离的关系图;
图3为实施例1中,将平板探测器移动到所需成像最近的距离位置zmin时,该X-射线源在平板上所成图像示意图;
图4为实施例2的系统方框图。
具体实施方式
首先,结合相关原理对发明的实施1做出具体说明。
在近距离X‐射线成像时,对于点光源而言,我们知道被照射点点光源强度跟距光源距离平方成反比(如图2所示):
这里A是光源强度,z0是点光源离平板表面的最短垂直距离,I0=A/z0 2,r是点光源离平板某一点坐标(x,y)的距离:r2=x2|y2|=z0 2,θx和θy是该点(x,y) 与X‐射线源连线与Z轴分别在Y‐Z面和X‐Z面上投影的夹角:θx=tan-1(x/z0)≈x/z0,θy=tan-1(y/z0)y/z0。
由于Heel效应,实际X‐光源的强度在空间有一个分布,因此照射到探测器平板上是射线强度不是像公式一那样简单。这时X‐射线照射强度是理想点光源照射强度分布函数I1(x,y)乘以一个空间分布函数I2(x,y)。这里空间辐射强度分布函数我们用二次函数来近似表示:
或
这里t(x,y)=I(x,y)/I0是相对强度分布:t(0,0)=1,A,B,C,D,是二项式拟合系数。沿着X轴(θy=0)或Y轴(θx=0)上的强度分布分别是:
和
依据公式(3)和(4),如果求出系数A,B,C,D,我们就可以完全确定X‐射线强度在空间的分布。我们就有可能对射线强度空间分布做修正,从而得到强度均匀的投影图像。
首先将平板探测器正对X‐射线源,并将它移到足够远的位置(zmax>>L),这里L是平板的尺度。这时照射到平板上的射线几乎是均匀的。我们可以调节平板探测器每一个像素单元的灵敏度,使得输出的图像是一个均匀的图像(每一个像素单元输出的信号强度一样)。然后,将平板探测器移动到所需成像最近的距离位置zmin。这时,该X‐射线源在平板上所成都像如图3所示:
从所成的像中,我们可以找到X、Y轴对称,然后将沿X、Y轴的信号采集下来,如图3所示,用二次多项式函数去拟合X,Y轴上信号强度分布,从而得到系数A,B,C,D。这样当点源到平板探测器距离z满足:zmin≤z≤zmax,采用公式(2’),对平板探测器灵敏度进行修正,得到在任何距离处,都是均匀的图像。
进一步的,我们对发明的实施例2做出具体说明。
其实际上是一种更为实用的方法,先将平板探测器正对X‐射线源,先放在位置z1,调节平板探测器每一个像素单元的灵敏度,使得输出的信号能形成一个均匀的图像(每一个像素单元输出的图像信号强度一样);然后将平板探测器仍然正对X‐射线源,并移动到位置z2,保持每个像素单元的灵敏度不变,这时保持平板探测器就会收到一帧强度不均匀分布图像(类似于图3),沿X、Y轴上采集信号,该信号形成一条曲线,然后用二次多项式函数进行拟合,求出a0,b0,c0,和d0,
和
这里,t(x,y)=I(x,y)/I0是相对强度分布,t(0,0)=1,A,B,C,D是二项式拟合系数,θx和θy是该点(x,y)与X‐射线源连线与Z轴分别在Y‐Z面和X‐Z面上投影的夹角,
由于a0=A(1/z1-1/z2)、 c0=C(1/z1-1/z2)和 所以我们就可以求出系数A,B,C,D:A=a0×(1/z1-1/z2)-1、 C=c0×(1/z1-1/z2)-1和
如果将平板探测器移到任何一个待矫正的位置“z”,这时我们可以得到新的一组系数:a=A(1/z1-1/z)、c=C(1/z1-1/z)和 然后将在z处的采集的图像信号乘以以下函数:
Δt(x,y)=(1|ax|bx2)×(1|cy|dy2) (7)
这样,我们就可以得到矫正之后的图像。具体的实现过程,见系统方框图4。
Claims (6)
1.X-射线成像空间分布不均匀性的校正方法,其特征在于至少包括如下步骤:
将平板探测器放在位置z1并正对X-射线源,调节平板探测器每一个像素单元的灵敏度,使得输出的信号形成一个均匀的图像;
将平板探测器移动到位置z2且仍然正对X-射线源,这时平板探测器就会收到一帧强度不均匀分布图像;
将位置z2上这一帧强度不均匀分布图像沿X、Y轴采集信号,并将该信号用二次多项式函数进行拟合得到图像校正系数;
将平板探测器移动到任何一个待矫正的位置z,通过二次多项式函数及图像校正系数得到位置z上的均匀图像。
2.根据权利要求1所述的X-射线成像空间分布不均匀性的校正方法,其特征在于:位置z1为平板探测器距离X-射线源足够远的位置zmax>>L,这里L是平板的尺度,位置z2为平板探测器所需成像最近的距离位置zmin。
3.根据权利要求2所述的X-射线成像空间分布不均匀性的校正方法,其特征在于:位置z上的均匀图像 这里t(x,y)=I(x,y)/I0是相对强度分布,t(0,0)=1,A,B,C,D是二项式拟合系数,θx和θy是该点(x,y)与X‐射线源连线与Z轴分别在Y‐Z面和X‐Z面上投影的夹角,z满足:zmin≤z≤zmax。
4.根据权利要求3所述的X-射线成像空间分布不均匀性的校正方法,其特征在于:沿着X轴(θy=0)或Y轴(θx=0)上的强度分布分别是:
和
依据上述两个公式求出二项式拟合系数A,B,C,D。
5.根据权利要求1所述的X-射线成像空间分布不均匀性的校正方法,其特征在于:将位置z2上这一帧强度不均匀分布图像沿X、Y轴采集信号并用二次多项式函数进行拟合,求出a0,b0,c0,和d0,
和
由于a0=A(1/z1-1/z2)、 c0=C(1/z1-1/z2)和 因此可以求出二项式拟合系数A,B,C,D:A=a0×(1/z1-1/z2)-1、 C=c0×(1/z1-1/z2)-1和
6.根据权利要求5所述的X-射线成像空间分布不均匀性的校正方法,其特征在于:将平板探测器移到任何一个待矫正的位置z,这时可以得到一组系数:a=A(1/z1-1/z)、 c=C(1/z1-1/z)和 然后将在z处的采集的图像信号乘以以下函数:
Δt(x,y)=(1|ax|bx2)×(1|cy|dy2) 这样,就可以得到矫正之后的图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410036628.7A CN103995010B (zh) | 2014-01-26 | 2014-01-26 | X-射线成像空间分布不均匀性的校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410036628.7A CN103995010B (zh) | 2014-01-26 | 2014-01-26 | X-射线成像空间分布不均匀性的校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103995010A true CN103995010A (zh) | 2014-08-20 |
CN103995010B CN103995010B (zh) | 2016-03-23 |
Family
ID=51309233
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410036628.7A Active CN103995010B (zh) | 2014-01-26 | 2014-01-26 | X-射线成像空间分布不均匀性的校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103995010B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1530074A (zh) * | 2003-03-18 | 2004-09-22 | 株式会社岛津制作所 | 放射线图像拾取装置 |
WO2007141689A1 (en) * | 2006-06-02 | 2007-12-13 | Philips Intellectual Property & Standards Gmbh | X-ray image apparatus and device for and method of calibrating an x-ray image apparatus |
US20080123816A1 (en) * | 2004-03-29 | 2008-05-29 | National Institute Of Radiological Sciences | Heel Effect Compensation Filter X-Ray Irradiator, X-Ray Ct Scanner and Method for X-Ray Ct Imaging |
CN101647707A (zh) * | 2008-08-11 | 2010-02-17 | 株式会社东芝 | X射线计算机断层摄影装置 |
US20110013748A1 (en) * | 2009-07-16 | 2011-01-20 | Yokogawa Electric Corporation | Radiation inspection apparatus |
CN102068268A (zh) * | 2010-12-17 | 2011-05-25 | 陈建锋 | 一种利用多能量x射线复合投影数字合成成像的方法及其装置 |
-
2014
- 2014-01-26 CN CN201410036628.7A patent/CN103995010B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1530074A (zh) * | 2003-03-18 | 2004-09-22 | 株式会社岛津制作所 | 放射线图像拾取装置 |
US20080123816A1 (en) * | 2004-03-29 | 2008-05-29 | National Institute Of Radiological Sciences | Heel Effect Compensation Filter X-Ray Irradiator, X-Ray Ct Scanner and Method for X-Ray Ct Imaging |
WO2007141689A1 (en) * | 2006-06-02 | 2007-12-13 | Philips Intellectual Property & Standards Gmbh | X-ray image apparatus and device for and method of calibrating an x-ray image apparatus |
CN101647707A (zh) * | 2008-08-11 | 2010-02-17 | 株式会社东芝 | X射线计算机断层摄影装置 |
US20110013748A1 (en) * | 2009-07-16 | 2011-01-20 | Yokogawa Electric Corporation | Radiation inspection apparatus |
CN102068268A (zh) * | 2010-12-17 | 2011-05-25 | 陈建锋 | 一种利用多能量x射线复合投影数字合成成像的方法及其装置 |
Also Published As
Publication number | Publication date |
---|---|
CN103995010B (zh) | 2016-03-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8077909B2 (en) | Apparatus and method for testing infrared camera | |
US9517045B2 (en) | Radiographic imaging apparatus and a method of correcting threshold energies in a photon-counting radiographic detector | |
US10054555B2 (en) | X-ray transmission inspection apparatus and inspection method using the same | |
EP3342342B1 (en) | Radiation image processing method and radiographic system | |
Huet et al. | Characterization of the gafchromic EBT3 films for dose distribution measurements in stereotactic radiotherapy | |
Blake et al. | Characterization of a novel EPID designed for simultaneous imaging and dose verification in radiotherapy | |
Takabayashi et al. | New method for measuring beam profiles using a parametric X-ray pinhole camera | |
CN103576179A (zh) | 用于温度稳定化的方法、x射线辐射探测器和ct系统 | |
CN102335002A (zh) | Ct机x光发生器及探测器位置校准工具,校准工具的校准方法和ct系统校准方法 | |
CN105717532A (zh) | 闪烁体的空间增益分布的确定 | |
US20150276626A1 (en) | X-Ray Transmission Inspection Apparatus | |
Méndez | Model selection for radiochromic film dosimetry | |
Tamborini et al. | Development and characterization of a 2D scintillation detector for quality assurance in scanned carbon ion beams | |
Bergamaschi et al. | Photon counting microstrip detector for time resolved powder diffraction experiments | |
CN103995010B (zh) | X-射线成像空间分布不均匀性的校正方法 | |
CN105374014A (zh) | 图像校正方法及装置、医学图像生成方法及装置 | |
CN204314473U (zh) | 快速获得免冲洗验证胶片剂量响应曲线的装置 | |
CN106405621B (zh) | 一种焦点尺寸测量方法 | |
CN103091968B (zh) | 具有自动校正功能的x-射线成像方法 | |
US20210165112A1 (en) | Method for estimating a dose rate on the basis of a spectral image | |
Norlin et al. | Characterisation of the charge sharing in pixellated Si detectors with single-photon processing readout | |
KR101629934B1 (ko) | X선 장치의 산란선 분포 측정 시스템, 방법 및 그 방법을 기록한 기록매체 | |
CN104055530A (zh) | 一种医用x射线束光器的预览方法 | |
CN102949193B (zh) | 荧光检查设备及其降低荧光检查剂量的方法 | |
Yokoyama et al. | Ionoluminescence analysis of glass scintillators and application to single-ion-hit real-time detection |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |