发明名称 一种基于电子倍增CCD的扩散光学层析成像系统及方法
摘要 本发明属于生理参数测量技术领域,涉及一种基于电子倍增CCD的扩散光学层析成像系统,包括稳态激光光源、平移台,耦合透镜,镜头,电子倍增CCD和计算机,由稳态激光光源发射出的稳态激光,经过入射光纤照射到仿体或待测物体上,由平移台带动光纤头移动,使得激光照射到仿体的各个成像点,之后镜头接受经由仿体而出的透射光,使其在电子倍增CCD中成像,成像结果被传送至计算机,由计算机实现图像重建。本发明同时提供一种采用上述系统实现的成像方法。本发明可以提高重建图像的成像质量,最大可能消除成像算法呈现出的欠定的病态特性,提高图像重建的效率和精度。
申请公布号 CN102579011B 申请公布日期 2013.10.16
申请号 CN201210051377.0 申请日期 2012.03.01
申请人 天津大学 发明人 赵会娟;周晓青;侯强;范颖;高峰
分类号 A61B5/00(2006.01)I;G01N21/35(2006.01)I 主分类号 A61B5/00(2006.01)I
代理机构 天津市北洋有限责任专利代理事务所 12201 代理人 程毓英
主权项 1.一种基于电子倍增CCD的扩散光学层析成像系统实现的图像重建方法,所述的基于电子倍增CCD的扩散光学层析成像系统,包括稳态激光光源、平移台、耦合透镜、镜头、电子倍增CCD和计算机,由稳态激光光源发射出的稳态激光,经过入射光纤照射到仿体上,由平移台带动光纤头移动,使得激光照射到仿体的各个成像点,之后镜头接收经由仿体而出的透射光,使其在电子倍增CCD中成像,成像结果被传送至计算机,由计算机实现图像重建;该图像重建方法假设光源的个数即光源随平移台移动扫描的点数为S×S,电子倍增CCD获取的图像有D×D个像素点,即有D×D个探测点,系统结束图像采集时共得到S×S幅均匀平板仿体的图像和S×S幅待测平板仿体的图像,包括下列步骤:1)根据光源位置的不同,将每个光源照射下的均匀平板仿体图像和待测平板仿体图像合并为一组,此时可得到S×S组图像,每组图像中包含有一幅均匀平板仿体的图像和一幅待测平板仿体的图像,由于图像的像素点的幅值代表着该位置的出射光强,所以每幅图像共可得到D×D个出射光强,这里用I<sub>0</sub>(r<sub>s</sub>,r<sub>d</sub>)表示光源在光源位置r<sub>s</sub>处照射、在探测点位置r<sub>d</sub>处探测得到的经过均匀平板仿体的出射光强,I(r<sub>s</sub>,r<sub>d</sub>)表示光源在光源位置r<sub>s</sub>处照射、在探测点位置r<sub>d</sub>处探测得到的经过待测平板仿体的出射光强;2)将在同一个光源位置r<sub>s</sub>处,每组待测平板仿体和均匀平板仿体的图像中对应同一探测点r<sub>d</sub>的出射光强作比,得到相对量<img file="FDA00003621938300011.GIF" wi="441" he="156" />此时,对于每组图像,共可以得到D×D个T(r<sub>s</sub>,r<sub>d</sub>)值,表现为一个D×D的矩阵M<sub>0</sub>,矩阵的每个元素代表对应于(r<sub>s</sub>,r<sub>d</sub>)所测量到的相对量T(r<sub>s</sub>,r<sub>d</sub>)的值,由于共有S×S组图像对,可得到S×S个D×D矩阵,之后将这S×S个D×D矩阵组合成一个以S×S为列数,D×D为行数的矩阵M;3)基于格林函数的微扰理论和一阶利托夫近似,得到中间变量φ(r<sub>s</sub>,r<sub>d</sub>)=-G<sub>0</sub>(r<sub>s</sub>,r<sub>d</sub>)lnM,其中,G<sub>0</sub>(r<sub>s</sub>,r<sub>d</sub>)为模拟计算得到的准确值,代表光源在r<sub>s</sub>处照射、在探测点r<sub>d</sub>处探测得到的均质格林函数解,中间变量φ(r<sub>s</sub>,r<sub>d</sub>)是一个以S×S为列数,D×D为行数的矩阵;4)根据均匀平板仿体,建立三维平面直角坐标系XYZ,将中间变量φ(r<sub>s</sub>,r<sub>d</sub>)中的光源位置r<sub>s</sub>和探测器位置r<sub>d</sub>,通过所建立的坐标系表示为X、Y和Z方向上的分量,在给定的Z坐标下,分别对光源位置r<sub>s</sub>和探测器位置r<sub>d</sub>的X方向和Y方向的分量做傅里叶变换,得到光源位置r<sub>s</sub>和探测器位置r<sub>d</sub>的X方向和Y方向相应的频谱信息;5)然后根据光源位置r<sub>s</sub>和探测器位置r<sub>d</sub>的X方向和Y方向相应的频谱信息,建立一维逆积分方程,求得吸收系数的变化量δμ<sub>a</sub>(r)中X方向和Y方向相应的频谱信息,再对δμ<sub>a</sub>(r)中X方向和Y方向上的分量做傅里叶逆变换,即可求得δμ<sub>a</sub>(r)在X方向和Y方向的分量;6)求出任意Z坐标下δμ<sub>a</sub>(r)在X方向和Y方向的分量,完成对待测平板仿体的图像重建过程。
地址 300072 天津市南开区卫津路92号