《基于光耦探测器X射线三维显微镜重建体素尺寸标定方法.pdf》由会员分享，可在线阅读，更多相关《基于光耦探测器X射线三维显微镜重建体素尺寸标定方法.pdf（9页珍藏版）》请在专利查询网上搜索。
1、10申请公布号CN104048600A43申请公布日20140917CN104048600A21申请号201410263478322申请日20140612G01B11/0020060171申请人天津大学地址300072天津市南开区卫津路92号72发明人胡小唐胡晓东邹晶须颖赵金涛牛鑫74专利代理机构天津市北洋有限责任专利代理事务所12201代理人刘国威54发明名称基于光耦探测器X射线三维显微镜重建体素尺寸标定方法57摘要本发明涉及X射线CT技术领域，为提供重建图像像素尺寸的标定方法。为此，本发明采取的技术方案是，基于光耦探测器X射线三维显微镜重建体素尺寸标定方法，首先利用标准栅格板进行高分辨显微。
2、CT的光学放大标定，然后利用标准的球棒进行三维显微镜整体放大倍数标定，通过计算获取几何放大倍数，进一步建立几何放大倍数与射线源、样品台和探测器相对于各自零位的位置的关系的公式，根据该公式即可简单方便的计算射线源和探测器位置发生变化的几何放大比，利用CCD像素尺寸除以几何放大比和光学放大比的乘积，得到最终的体素尺寸。本发明主要应用于X射线CT设备的设计、制造。51INTCL权利要求书1页说明书5页附图2页19中华人民共和国国家知识产权局12发明专利申请权利要求书1页说明书5页附图2页10申请公布号CN104048600ACN104048600A1/1页21一种基于光耦探测器X射线三维显微镜重建体。
3、素尺寸标定方法，其特征是，首先利用标准栅格板进行高分辨显微CT的光学放大标定，然后利用标准的球棒进行三维显微镜整体放大倍数标定，通过计算获取几何放大倍数，进一步建立几何放大倍数与射线源、样品台和探测器相对于各自零位的位置的关系的公式，根据该公式即可简单方便的计算射线源和探测器位置发生变化的几何放大比，利用CCD像素尺寸除以几何放大比和光学放大比的乘积，得到最终的体素尺寸。2如权利要求1所述的基于光耦探测器X射线三维显微镜重建体素尺寸标定方法，其特征是，利用标准栅格板进行高分辨显微CT的光学放大标定具体为步骤101将光耦探测器前端的闪烁体更换为栅格板，通过调节栅格板的位置，使其在CCD上成清晰像。
4、，此时栅格板所在位置即为原先闪烁体所在位置；步骤102通过CANNY算子对栅格板显微图像进行边缘提取，计算出平行栅格之间的距离；步骤103利用栅格之间的距离进行光学放大倍数的标定利用平行栅格之间的距离除以平行栅格之间的实际距离，即得光学放大倍数；步骤104求取多行多列的光学放大倍数，将结果取平均得最终的光学放大倍数。3如权利要求1所述的基于光耦探测器X射线三维显微镜重建体素尺寸标定方法，其特征是，利用标准的球棒进行三维显微镜整体放大倍数标定具体为步骤201采集标准球棒的显微CT投影图像；步骤202使用CANNY算子提取标准球棒中双球的投影图像的边缘；步骤203使用HOUGH变换确定标准球棒中双。
5、球的圆心的位置，并且计算出球心距离；步骤204利用计算出的球心距离除以标准球棒的实际距离，即得总放大倍数。4如权利要求1所述的基于光耦探测器X射线三维显微镜重建体素尺寸标定方法，其特征是，通过计算获取几何放大倍数具体为步骤301用总放大倍数除以前面标定的光学放大倍数，即得此时的几何放大倍数；步骤302标定出的几何放大倍数和此时射线源、样品台和探测器相对于各自零位的位置，解算出它们处于各自零位时的SDD0和SOD0；步骤303建立几何放大倍数与各位置读数的公式如下其中，LS，LO和LD分别为射线源、样品台和探测器成像时相对于各自零位的坐标，由光栅尺给出；GEO为此时的几何放大倍数；SDD0和SO。
6、D0为射线源、样品台和探测器在各自零位时的射线源到探测器中心距离SDDDISTANCEOFSOURCETODETECTOR和射线源到转台中心距离SODDISTANCEOFSOURCETOOBJECT；步骤304通过读取射线源、样品台和探测器成像时相对于各自零位的坐标，根据公式1即可计算出几何放大比。5如权利要求1所述的基于光耦探测器X射线三维显微镜重建体素尺寸标定方法，其特征是，体素尺寸计算具体为步骤401利用CCD像素尺寸除以几何放大比和光学放大比的乘积。权利要求书CN104048600A1/5页3基于光耦探测器X射线三维显微镜重建体素尺寸标定方法技术领域0001本发明涉及X射线CT技术领域。
7、，具体讲，涉及基于光耦探测器X射线三维显微镜重建体素尺寸标定方法。技术背景0002X射线CT是无损检测领域的重要技术手段，其分辨力一般处于毫米量级。基于光耦探测器的X射线三维显微镜以其高的空间分辨力，在MEMS器件封装和装配误差分析、半导体器件封装和内部缺陷检测、石油地质勘探等方面起到了越来越重要的作用。0003X射线显微镜扫描成像过程是射线源发射出的X射线穿透所述待测样品后投射至所述探测器上，完成一幅投影图像的采集；通过在360度范围内的系列均匀采样，并利用相应的成像算法即可获取三维立体图像。0004分辨率是X射线显微镜的一个非常重要的技术指标。而这个指标是由体素尺寸决定。体素尺寸越小，分辨。
8、率越高。X射线显微镜三维重建结果可以视为一系列代表物质衰减系数的小立方体组成，这些小立方是以像素为单位的，因此无法直接得到体素尺寸。将像素尺寸转换为实际的体素尺寸是X射线显微镜成像中的重要过程。在实际扫描成像过程中，射线源、样品台和探测器的相对位置发生变化时，体素的尺寸也在发生变化，因此自动标定体素尺寸十分必要。0005中国专利申请号2013101094603提出了“锥束3DCT扫描系统重建体素尺寸的自动标定方法”。该专利不需要进行几何放大比的直接计算、校正程序简单、实用。但是该方法是基于平板探测器的X射线CT系统，而且只适用于射线源和探测器位置固定的情况。本专利则是针对基于光耦探测器的高分辨。
9、显微CT系统，关于其标定方法目前国内外的文献资料中并未查找到相关报道。发明内容0006为克服现有技术的不足，提供重建图像像素尺寸的标定方法。为此，本发明采取的技术方案是，基于光耦探测器X射线三维显微镜重建体素尺寸标定方法，首先利用标准栅格板进行高分辨显微CT的光学放大标定，然后利用标准的球棒进行三维显微镜整体放大倍数标定，通过计算获取几何放大倍数，进一步建立几何放大倍数与射线源、样品台和探测器相对于各自零位的位置的关系的公式，根据该公式即可简单方便的计算射线源和探测器位置发生变化的几何放大比，利用CCD像素尺寸除以几何放大比和光学放大比的乘积，得到最终的体素尺寸。0007利用标准栅格板进行高分。
10、辨显微CT的光学放大标定具体为0008步骤101将光耦探测器前端的闪烁体更换为栅格板，通过调节栅格板的位置，使其在CCD上成清晰像，此时栅格板所在位置即为原先闪烁体所在位置；0009步骤102通过CANNY算子对栅格板显微图像进行边缘提取，计算出平行栅格之间的距离；说明书CN104048600A2/5页40010步骤103利用栅格之间的距离进行光学放大倍数的标定利用平行栅格之间的距离除以平行栅格之间的实际距离，即得光学放大倍数；0011步骤104求取多行多列的光学放大倍数，将结果取平均得最终的光学放大倍数。0012利用标准的球棒进行三维显微镜整体放大倍数标定0013步骤201采集标准球棒的显微。
11、CT投影图像；0014步骤202使用CANNY算子提取标准球棒中双球的投影图像的边缘；0015步骤203使用HOUGH变换确定标准球棒中双球的圆心的位置，并且计算出球心距离；0016步骤204利用计算出的球心距离除以标准球棒的实际距离，即得总放大倍数。0017通过计算获取几何放大倍数具体为0018步骤301用总放大倍数除以前面标定的光学放大倍数，即得此时的几何放大倍数；0019步骤302标定出的几何放大倍数和此时射线源、样品台和探测器相对于各自零位的位置，解算出它们处于各自零位时的SDD0和SOD0；0020步骤303建立几何放大倍数与各位置读数的公式如下00210022其中，LS，LO和LD。
12、分别为射线源、样品台和探测器成像时相对于各自零位的坐标，由光栅尺给出；GEO为此时的几何放大倍数；SDD0和SOD0为射线源、样品台和探测器在各自零位时的射线源到探测器中心距离SDDDISTANCEOFSOURCETODETECTOR和射线源到转台中心距离SODDISTANCEOFSOURCETOOBJECT；0023步骤304通过读取射线源、样品台和探测器成像时相对于各自零位的坐标，根据公式1即可计算出几何放大比。0024体素尺寸计算具体为0025步骤401利用CCD像素尺寸除以几何放大比和光学放大比的乘积。0026与已有技术相比，本发明的技术特点与效果0027本发明的优点在于1系统确定后，。
13、光学放大就确定了，因此光学放大倍数只进行一次标定即可。2由于系统整体的放大倍数是根据标准球棒的球心距进行的标定，避免了球边界界定引起的误差，标定的放大比精确度较高。3系统确定后，射线源和探测器位置可以根据实际需要进行位置调节，而且几何放大的标定只需要将射线源、样品台和探测器相对于各自零位的位置代入公式即可。4本文方法简单，不需要进行CT扫描标定，只需要对标准的栅格板、标准球棒进行若干次投影数据采集即可。附图说明0028图1基于光耦探测器的显微CT结构示意。0029图2为光耦探测器部分的示意图。0030图3为栅格板的DR图像。0031图4为球棒的DR图像。0032图5为球棒的轮廓图像。说明书CN。
14、104048600A3/5页5具体实施方式0033本发明的目的在于提供一种基于光耦探测器的高分辨显微CT系统的重建图像像素尺寸的标定方法。利用标准的球棒首先进行光学放大倍数标定，然后进行三维显微镜整体放大倍数标定，通过计算获取几何放大倍数，进一步建立几何放大倍数与射线源、样品台和探测器相对于各自零位的位置的关系的公式。根据该公式即可简单方便的计算射线源和探测器位置发生变化的几何放大比。利用CCD像素尺寸除以几何放大比和光学放大比的乘积，得到最终的体素尺寸。0034本发明的技术方案如下0035步骤一进行光学放大倍数标定0036步骤101将光耦探测器前端的闪烁体更换为栅格板，通过调节栅格板的位置，。
15、使其在CCD上成清晰像。此时栅格板所在位置即为原先闪烁体所在位置，这样可以保证标定得到的光学放大倍数与实际使用时的一致。0037步骤102通过CANNY算子对栅格板显微图像进行边缘提取，计算出平行栅格之间的距离。0038步骤103利用栅格之间的距离进行光学放大倍数的标定。利用平行栅格之间的距离除以平行栅格之间的实际距离，即得光学放大倍数。0039步骤104求取多行多列的光学放大倍数，将结果取平均得最终的光学放大倍数。0040步骤二三维显微镜整体放大倍数标定0041步骤201采集标准球棒的显微CT投影图像。0042步骤202使用CANNY算子提取标准球棒中双球的投影图像的边缘。0043步骤203。
16、使用HOUGH变换确定标准球棒中双球的圆心的位置，并且计算出球心距离。0044步骤204利用计算出的球心距离除以标准球棒的实际距离，即得总放大倍数。0045步骤三几何放大倍数标定0046步骤301用总放大倍数除以前面标定的光学放大倍数，即得此时的几何放大倍数。0047步骤302标定出的几何放大倍数和此时射线源、样品台和探测器相对于各自零位的位置，解算出它们处于各自零位时的SDD0和SOD0。0048步骤303建立几何放大倍数与各位置读数的公式如下00490050其中，LS，LO和LD分别为射线源、样品台和探测器成像时相对于各自零位的坐标，由光栅尺给出；GEO为此时的几何放大倍数；SDD0和SO。
17、D0为射线源、样品台和探测器在各自零位时的SDD射线源到探测器中心距离，DISTANCEOFSOURCETODETECTOR和SOD射线源到转台中心距离，DISTANCEOFSOURCETOOBJECT。0051步骤304通过读取射线源、样品台和探测器成像时相对于各自零位的坐标，根据公式1即可计算出几何放大比。0052步骤四体素尺寸计算说明书CN104048600A4/5页60053步骤401利用CCD像素尺寸除以几何放大比和光学放大比的乘积。0054本发明的特征在于包括以下步骤00551利用标准栅格板进行高分辨显微CT的光学放大标定。00562利用标准球棒标定三维显微镜整体放大倍数标定。3通。
18、过公式直接计算射线源、样品台和探测器位置变化后的几何放大倍数。00574简单方便的获取重建图像像素尺寸。0058本发明的优点在于1系统确定后，光学放大就确定了，因此光学放大倍数只进行一次标定即可。2由于系统整体的放大倍数是根据标准球棒的球心距进行的标定，避免了球边界界定引起的误差，标定的放大比精确度较高。3系统确定后，射线源和探测器位置可以根据实际需要进行位置调节，而且几何放大的标定只需要将射线源、样品台和探测器相对于各自零位的位置代入公式即可。4本文方法简单，不需要进行CT扫描标定，只需要对标准的栅格板、标准球棒进行若干次投影数据采集即可。0059结合实施例对本发明做进一步说明，具体步骤如下。
19、00601将图X中光耦探测器前端的闪烁体更换为如图3所示的标准栅格板黑线中心间距为100M，调整栅格板的位置使其在CCD上成清晰像。00612通过CANNY算子对栅格板显微图像进行边缘提取，获取平行栅格直接的距离。00623根据CCD像素尺寸，计算网格放大后尺寸，再除以网格实际尺寸，即得光学放大倍数。00634将栅格板更换为原来的闪烁片。00645将球棒放置在样品台中心位置，获取球棒的投影图像。00656使用CANNY算子提取图像边缘。00667使用HOUGH变换确定标准球棒中双球的圆心的位置，并且计算出球心距离。利用计算出的球心距除以球棒的两个球心的实际距离，即得总放大倍数00677用总放大。
20、倍数除以前面标定的光学放大倍数，即得此时的几何放大倍数。00688选取了3组射线源和探测器的位置，每组中对样品在光轴上12个不同位置进行投影，采用最小二乘法求解SDD0和SOD0，结果如表2所示。0069记录射线源、样品台和探测器相对于各自零位的位置，解算出它们处于各自零位时的SDD0和SOD0。0070表1SDD0和SOD0标定数据单位MM007100729射线源、样品台和探测器可以进行调整，其几何放大比可以通过如下公式计算得出0073说明书CN104048600A5/5页70074其中，LS，LO和LD分别为射线源、样品台和探测器成像时相对于各自零位的坐标，由光栅尺给出；GEO为此时的几何放大倍数；SDD0和SOD0为射线源、样品台和探测器在各自零位时的SDD和SOD。007510利用CCD像素尺寸除以几何放大比和光学放大比的乘积，即得出像素尺寸。说明书CN104048600A1/2页8图1图2图3说明书附图CN104048600A2/2页9图4图5说明书附图CN104048600A。