CN112461872A - 基于双微球模体的显微ct焦斑漂移测量系统及补偿方法 - Google Patents
基于双微球模体的显微ct焦斑漂移测量系统及补偿方法 Download PDFInfo
- Publication number
- CN112461872A CN112461872A CN202011110009.XA CN202011110009A CN112461872A CN 112461872 A CN112461872 A CN 112461872A CN 202011110009 A CN202011110009 A CN 202011110009A CN 112461872 A CN112461872 A CN 112461872A
- Authority
- CN
- China
- Prior art keywords
- focal spot
- image
- double
- microsphere
- ray source
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- 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
- G01N23/046—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 using tomography, e.g. computed tomography [CT]
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B21/00—Microscopes
- G02B21/36—Microscopes arranged for photographic purposes or projection purposes or digital imaging or video purposes including associated control and data processing arrangements
- G02B21/365—Control or image processing arrangements for digital or video microscopes
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Multimedia (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- General Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biochemistry (AREA)
- Pulmonology (AREA)
- Theoretical Computer Science (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Radiology & Medical Imaging (AREA)
- Optics & Photonics (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明涉及一种基于双微球模体的显微CT焦斑漂移测量系统,包括:射线源、样品旋转台、由光学显微镜、CCD、置于光学显微镜前的闪烁体共同组成的CCD光学耦合探测器,其特征在于,还包括置于射线源前端的双微球模体单元和两个位移传感器单元,双微球模体单元的投影位于测量系统视场内,两个位移传感器单元,一个用于测量射线源的位移,一个用于测量CCD光学耦合探测器的位移。本发明同时给出采用此种测量系统实现的显微CT焦斑漂移测量和补偿方法。
Description
技术领域:
本发明涉及一种高分辨显微CT系统的焦斑漂移测量及补偿方法。
背景技术:
CT(Computed Tomography,计算机断层成像)技术是利用X射线穿透样品,穿透后的X射线携带样品信息被探测器接收成像,并利用重构算法实现样品内部三维结构检测的一种技术,空间分辨率是系统的关键指标之一。根据空间分辨率的不同,通常将CT分为工业CT、显微CT、纳米CT等,其中,显微CT是指空间分辨率介于0.5~10μm的CT系统。
其中,基于光学耦合探测器的显微CT技术是一种主流的高分辨显微CT技术,其基本原理是利用闪烁体将X射线转换成可见光,然后利用显微镜组将可见光放大成像至CCD芯片上,从而实现高放大比高分辨率的X射线成像。此类系统采用几何—光学的两级放大结构,可以有效提高系统空间分辨率,使系统空间分辨率达到微米/亚微米量级。但这放大了由于射线源焦斑漂移引起的图像偏移,严重影响系统成像质量。因此,对射线源焦斑漂移进行测量及补偿引起的图像偏移就尤为重要。
目前,用于对射线源焦斑漂移进行检测的方法有单球法、双球法、漂移向量法等,然而,以上方法均只对焦斑漂移进行定性检测,无法实现定量测量,且不能有效补偿焦斑漂移引起的图像偏移。因此,希望有一种技术方案来克服或至少减轻现有技术中的缺陷,以实现焦斑漂移的测量及补偿。
发明内容:
本发明的目的在于提供一种高分辨显微CT焦斑漂移测量及补偿方法来克服或至少减轻现有技术中的缺陷。技术方案如下:
一种基于双微球模体的显微CT焦斑漂移测量系统,包括:射线源、样品旋转台、由光学显微镜、CCD、置于光学显微镜前的闪烁体共同组成的CCD光学耦合探测器,其特征在于,还包括置于射线源前端的双微球模体单元和两个位移传感器单元,双微球模体单元的投影位于测量系统视场内,两个位移传感器单元,一个用于测量射线源的位移,一个用于测量CCD光学耦合探测器的位移。
采用所述的测量系统实现的显微CT焦斑漂移测量及补偿方法,该方法包括如下步骤:
(1)开启射线源,以不大于0.5°的步长旋转样品台,采集多帧高分辨率图像,称之为投影图像Pk,同时获取采集每幅图像时的两个位移传感器数据,分别记为DSk与DDk,其中DSk为射线源位移传感器数据,DDk为探测器位移传感器数据;;
(2)利用获得的高分辨率图像中的双微球模体图像及对应的位移传感器数据实现射线源焦斑漂移的定量化测量,包括:
①利用光学显微镜对双微球模体的球心距进行离线测量,记为Dcal;
②利用霍夫变换对投影图像Pk中双微球投影的圆心进行亚像素提取,并结合CCD像素尺寸获得双微球投影图像球心距Dmea;
③计算射线源焦斑到双微球模体的距离SPD,
其中,SDD为射线源焦斑到闪烁体的距离,Mopt为CCD光学耦合探测器的放大倍数;
④计算焦斑漂移量在图像上的放大倍数Mshift,
⑤以第一幅高分辨率图像为基准,计算投影图像中双微球中心的偏移量,记为(U‘k,V’k);
⑥结合两个位移传感器数据,对因机械偏移引起的图像偏移进行剔除,获取因射线源焦斑偏移引起的图像偏移,记为(Uk,Vk),其中,
⑦射线源焦斑漂移的定量化测量,其在Y轴和Z轴的偏移量记为(uk,vk)
(3)利用步骤(2)所获得的射线源焦斑漂移定量化测量数据,对获得的高分辨率图像中的样品图像位置进行修正,实现高分辨显微CT焦斑漂移的补偿,再对补偿后的高分辨率图像进行三维重建。
步骤(3)具体包括:
①计算因射线源焦斑漂移引起的样品图像偏移量,记为(Xk,Yk)
其中,SOD为射线源焦斑到样品旋转台中心旋转轴的距离;
②利用所得到样品图像偏移量对投影图像中的样品图像位置进行修正,实现高分辨显微CT焦斑漂移的补偿;
③对补偿后的投影图像进行三维重建。
本发明的技术特点和创新:利用双微球模体和电容式位移传感器对射线源焦斑漂移进行了定量化测量,并实现了射线源焦斑漂移引起图像偏移的补偿,从而提高系统的空间分辨率。
附图说明:
图1高分辨显微CT焦斑漂移测量系统及补偿流程图。
图2高分辨显微CT焦斑漂移测量系统结构图;具体包括:射线源1;双微球模体单元2;旋转台3;CCD光学耦合探测器4,包括闪烁体5,物镜6,管镜组7,反光镜8,CCD 9;电容式位移传感器单元10。
图3双微球模体单元示意图
具体实施方式:
本发明的目的在于提供一种高分辨显微CT焦斑漂移测量及补偿方法,从而提高显微CT重建结果的三维空间分辨率,所述的测量方法简单、实用、可靠。
本发明是通过下述技术方案加以实现的,一种高分辨显微CT焦斑漂移测量及补偿方法,其流程如图1所示,具体可通过如下方案实施:
步骤1:在图2所示的高分辨显微CT系统射线源1前端安装双微球模体单元2,并通过设计加工在碳纤维板上的四个顶丝调整其位置,使双微球在CCD 9上的投影位于视场内靠上区域;利用磁性表座安装固定电容式位移传感器单元10,使其正对射线源1、CCD光学耦合探测器4侧面,以监测射线源1、CCD光学耦合探测器4发生的机械偏移;如图2所示;
步骤2:开启射线源1,调整射线源1电压、电流、CCD光学耦合探测器4采集时间等相关参数,使其满足成像需求;通过闪烁体5将穿透双微球模体单元2及样品后的射线转化为可见光,利用CCD光学耦合探测器4进行光学放大并成像,获得高分辨扫描数据;以不大于0.5°的步长旋转旋转台3,带动样品旋转,获得不同角度下的高分辨扫描数据,记为Pk,k=1,……K,其中,k表示第k个角度;同时获取每帧图像对应下的电容式位移传感器10的数据,分别记为DSk与DDk,其中DSk为射线源位移传感器数据,DDk为探测器位移传感器数据;
步骤3:利用所述步骤2中获得的所述投影图像Pk中的双微球模体图像及对应的电容式位移传感器数据实现射线源焦斑漂移的定量化测量;“定量化测量”具体实施如下:
步骤31:利用光学共聚焦显微镜对双微球模体的球心距进行离线测量,记为Dcal;
步骤32:利用霍夫变换对投影图像Pk中双微球投影的圆心进行亚像素提取,并结合CCD像素尺寸获得双微球投影图像球心距Dmea;
步骤33:计算射线源焦斑到双微球模体的距离SPD,
其中,SDD为射线源焦斑到闪烁体的距离,Mopt为光学显微镜组的放大倍数;
步骤34:计算焦斑漂移量在图像上的放大倍数Mshift,
步骤35:以第一幅图像为基准,计算投影图像中双微球中心的偏移量,记为(U‘k,V’k);
步骤36:结合电容式位移传感器数据,对因机械偏移引起的图像偏移进行剔除,获取因射线源焦斑漂移引起的图像偏移,记为(Uk,Vk),其中,
步骤37:射线源焦斑漂移的定量化测量,其在Y轴和Z轴的偏移量记为(uk,vk)
步骤4:利用所述步骤3中获得的射线源焦斑漂移定量化测量数据,对所述步骤2中获得的所述高分辨扫描数据Pk中的样品图像位置进行修正,实现高分辨显微CT焦斑漂移的补偿;再对补偿后的投影图像进行三维重建。“补偿”具体实施如下:
步骤41:计算因焦斑漂移引起的样品图像偏移量,记为(Xk,Yk)
步骤42:利用步骤41得出的偏移量对投影图像中的样品图像位置进行修正,实现高分辨显微CT焦斑漂移的补偿。
步骤43:对补偿后的投影图进行三维重建。
Claims (3)
1.一种基于双微球模体的显微CT焦斑漂移测量系统,包括:射线源、样品旋转台、由光学显微镜、CCD、置于光学显微镜前的闪烁体共同组成的CCD光学耦合探测器,其特征在于,还包括置于射线源前端的双微球模体单元和两个位移传感器单元,双微球模体单元的投影位于测量系统视场内,两个位移传感器单元,一个用于测量射线源的位移,一个用于测量CCD光学耦合探测器的位移。
2.采用权利要求1所述的测量系统实现的显微CT焦斑漂移测量及补偿方法,该方法包括如下步骤:
(1)开启射线源,以不大于0.5°的步长旋转样品台,采集多帧高分辨率图像,称之为投影图像Pk,同时获取采集每幅图像时的两个位移传感器数据,分别记为DSk与DDk,其中DSk为射线源位移传感器数据,DDk为探测器位移传感器数据;;
(2)利用获得的高分辨率图像中的双微球模体图像及对应的位移传感器数据实现射线源焦斑漂移的定量化测量,包括:
①利用光学显微镜对双微球模体的球心距进行离线测量,记为Dcal;
②利用霍夫变换对投影图像Pk中双微球投影的圆心进行亚像素提取,并结合CCD像素尺寸获得双微球投影图像球心距Dmea;
③计算射线源焦斑到双微球模体的距离SPD,
其中,SDD为射线源焦斑到闪烁体的距离,Mopt为CCD光学耦合探测器的放大倍数;
④计算焦斑漂移量在图像上的放大倍数Mshift,
⑤以第一幅高分辨率图像为基准,计算投影图像中双微球中心的偏移量,记为(U‘k,V’k);
⑥结合两个位移传感器数据,对因机械偏移引起的图像偏移进行剔除,获取因射线源焦斑偏移引起的图像偏移,记为(Uk,Vk),其中,
⑦射线源焦斑漂移的定量化测量,其在Y轴和Z轴的偏移量记为(uk,vk)
(3)利用步骤(2)所获得的射线源焦斑漂移定量化测量数据,对获得的高分辨率图像中的样品图像位置进行修正,实现高分辨显微CT焦斑漂移的补偿,再对补偿后的高分辨率图像进行三维重建。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011110009.XA CN112461872B (zh) | 2020-10-16 | 2020-10-16 | 基于双微球模体的显微ct焦斑漂移测量系统及补偿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011110009.XA CN112461872B (zh) | 2020-10-16 | 2020-10-16 | 基于双微球模体的显微ct焦斑漂移测量系统及补偿方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112461872A true CN112461872A (zh) | 2021-03-09 |
CN112461872B CN112461872B (zh) | 2023-01-24 |
Family
ID=74833581
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011110009.XA Active CN112461872B (zh) | 2020-10-16 | 2020-10-16 | 基于双微球模体的显微ct焦斑漂移测量系统及补偿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112461872B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6185275B1 (en) * | 1998-08-25 | 2001-02-06 | General Electric Company | Systems and methods for correcting focal spot thermal drift |
CN101515370A (zh) * | 2009-03-06 | 2009-08-26 | 北京航空航天大学 | 三维显微ct扫描系统中射线源焦点的投影坐标的标定方法 |
US20120328076A1 (en) * | 2011-06-23 | 2012-12-27 | Abdelaziz Ikhlef | Systems and methods for focal spot motion correction |
CN103759681A (zh) * | 2014-01-10 | 2014-04-30 | 天津大学 | 显微ct转轴运动误差校正方法 |
CN104048600A (zh) * | 2014-06-12 | 2014-09-17 | 天津大学 | 基于光耦探测器x射线三维显微镜重建体素尺寸标定方法 |
CN106403811A (zh) * | 2016-05-07 | 2017-02-15 | 天津大学 | 基于光耦探测器x射线三维显微镜系统参数标定方法 |
CN107684435A (zh) * | 2017-08-16 | 2018-02-13 | 深圳先进技术研究院 | 锥束ct系统几何校准方法及其校准装置 |
CN109685877A (zh) * | 2018-12-27 | 2019-04-26 | 重庆大学 | 一种基于自适应投影图像特征区域匹配的微纳ct焦点漂移校正方法 |
CN110068918A (zh) * | 2019-03-26 | 2019-07-30 | 北京航空航天大学 | 基于叠加双微球透镜的光学超分辨率成像系统和方法 |
-
2020
- 2020-10-16 CN CN202011110009.XA patent/CN112461872B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6185275B1 (en) * | 1998-08-25 | 2001-02-06 | General Electric Company | Systems and methods for correcting focal spot thermal drift |
CN101515370A (zh) * | 2009-03-06 | 2009-08-26 | 北京航空航天大学 | 三维显微ct扫描系统中射线源焦点的投影坐标的标定方法 |
US20120328076A1 (en) * | 2011-06-23 | 2012-12-27 | Abdelaziz Ikhlef | Systems and methods for focal spot motion correction |
CN103759681A (zh) * | 2014-01-10 | 2014-04-30 | 天津大学 | 显微ct转轴运动误差校正方法 |
CN104048600A (zh) * | 2014-06-12 | 2014-09-17 | 天津大学 | 基于光耦探测器x射线三维显微镜重建体素尺寸标定方法 |
CN106403811A (zh) * | 2016-05-07 | 2017-02-15 | 天津大学 | 基于光耦探测器x射线三维显微镜系统参数标定方法 |
CN107684435A (zh) * | 2017-08-16 | 2018-02-13 | 深圳先进技术研究院 | 锥束ct系统几何校准方法及其校准装置 |
CN109685877A (zh) * | 2018-12-27 | 2019-04-26 | 重庆大学 | 一种基于自适应投影图像特征区域匹配的微纳ct焦点漂移校正方法 |
CN110068918A (zh) * | 2019-03-26 | 2019-07-30 | 北京航空航天大学 | 基于叠加双微球透镜的光学超分辨率成像系统和方法 |
Non-Patent Citations (3)
Title |
---|
JINTAO ZHAO等: "Method for correction of rotation errors in Micro-CT System", 《NUCLEAR INSTRUMENTS AND METHODS IN PHYSICS RESEARCH A》 * |
蔡玉芳等: "基于投影图像特征匹配的微纳CT焦点漂移校正", 《仪器仪表学报》 * |
赵金涛: "显微CT系统标定关键技术的研究", 《中国博士学位论文全文数据库 (信息科技辑)》 * |
Also Published As
Publication number | Publication date |
---|---|
CN112461872B (zh) | 2023-01-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7974379B1 (en) | Metrology and registration system and method for laminography and tomography | |
CN101487809B (zh) | 零点定标方法及在光学微扫描显微热成像系统中的应用 | |
KR101850320B1 (ko) | 단층 촬영 장치 및 오차 제거 방법 | |
CN106500965A (zh) | 基于ccd探测器的龙虾眼x射线光学元件聚焦性能测试装置与方法 | |
EP1313065A1 (de) | Verfahren und Vorrichtung zur Kalibrierung von und zur Bilderzeugung mit einer schwerkraftempfindlichen Röntgenaufnahmeeinrichtung | |
CN106679595B (zh) | 楔角球面透镜的中心偏差和楔角的测量方法 | |
US20100208969A1 (en) | Rotation center identifying method, ring artifact correction method, rotation center identifying apparatus, x-ray diagnostic apparatus, recording medium on which program for executing rotation center identification is recorded and recording medium on which program for executing ring artifact correction is recorded | |
CN107505340A (zh) | 一种陶瓷艺术品鉴定方法 | |
CN103759681B (zh) | 显微ct转轴运动误差校正方法 | |
CN101936720A (zh) | 一种适用于锥束xct系统的探测器扭转角的标定方法 | |
CN103654835A (zh) | 一种评估spect针孔准直器性能的测试装置及其测试方法 | |
CN104048600A (zh) | 基于光耦探测器x射线三维显微镜重建体素尺寸标定方法 | |
CN112461872B (zh) | 基于双微球模体的显微ct焦斑漂移测量系统及补偿方法 | |
CN106447729A (zh) | 一种基于坐标变换的二维数字图像相关补偿方法及二维光学数字图像相关引伸计 | |
JP2011145160A (ja) | マルチフォーカス検査装置及びマルチフォーカス検査方法 | |
US10133052B2 (en) | Image acquiring method and image acquiring apparatus using the same | |
CN201277864Y (zh) | 内方位元素及畸变测试仪 | |
CN106403811A (zh) | 基于光耦探测器x射线三维显微镜系统参数标定方法 | |
CN109709127A (zh) | 低散射x射线荧光ct成像方法 | |
US6118843A (en) | Quantitative stereoscopic radiography method | |
JPH10137234A (ja) | X線断層撮像装置 | |
JPH04307035A (ja) | X線撮影装置 | |
JP5572521B2 (ja) | X線ct装置、x線ct装置の画像再構成方法 | |
Gao et al. | Optical microscanning X-ray real-time imaging system | |
CN112446248A (zh) | 面部表情采集系统 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |