CN105640582B - 一种深部组织x射线激发多光谱断层成像系统及方法 - Google Patents
一种深部组织x射线激发多光谱断层成像系统及方法 Download PDFInfo
- Publication number
- CN105640582B CN105640582B CN201610118383.1A CN201610118383A CN105640582B CN 105640582 B CN105640582 B CN 105640582B CN 201610118383 A CN201610118383 A CN 201610118383A CN 105640582 B CN105640582 B CN 105640582B
- Authority
- CN
- China
- Prior art keywords
- imaging
- ray
- excitation
- light
- rays
- 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
- 230000005284 excitation Effects 0.000 title claims abstract description 75
- 238000000034 method Methods 0.000 title claims abstract description 22
- 238000002591 computed tomography Methods 0.000 title abstract description 4
- 238000002603 single-photon emission computed tomography Methods 0.000 title abstract description 3
- 238000003384 imaging method Methods 0.000 claims abstract description 123
- 230000003595 spectral effect Effects 0.000 claims abstract description 70
- 239000000463 material Substances 0.000 claims abstract description 52
- 238000004458 analytical method Methods 0.000 claims abstract description 22
- 238000001914 filtration Methods 0.000 claims abstract description 21
- 238000012545 processing Methods 0.000 claims abstract description 8
- 238000000701 chemical imaging Methods 0.000 claims abstract description 3
- 230000003287 optical effect Effects 0.000 claims description 29
- 239000011159 matrix material Substances 0.000 claims description 13
- 230000002902 bimodal effect Effects 0.000 claims description 9
- 238000012937 correction Methods 0.000 claims description 8
- 230000005540 biological transmission Effects 0.000 claims description 7
- 230000011664 signaling Effects 0.000 claims description 6
- 238000001514 detection method Methods 0.000 claims description 5
- 210000000056 organ Anatomy 0.000 claims description 4
- 230000003044 adaptive effect Effects 0.000 claims description 3
- 238000011017 operating method Methods 0.000 claims description 3
- 238000004321 preservation Methods 0.000 claims description 3
- 238000003325 tomography Methods 0.000 claims description 3
- 238000007405 data analysis Methods 0.000 claims description 2
- 230000004927 fusion Effects 0.000 claims description 2
- 238000000354 decomposition reaction Methods 0.000 claims 1
- 230000003760 hair shine Effects 0.000 claims 1
- 230000005855 radiation Effects 0.000 claims 1
- 230000035515 penetration Effects 0.000 description 8
- 229910014323 Lanthanum(III) bromide Inorganic materials 0.000 description 3
- XKUYOJZZLGFZTC-UHFFFAOYSA-K lanthanum(iii) bromide Chemical compound Br[La](Br)Br XKUYOJZZLGFZTC-UHFFFAOYSA-K 0.000 description 3
- 229910052761 rare earth metal Inorganic materials 0.000 description 3
- 238000004875 x-ray luminescence Methods 0.000 description 3
- 229910052691 Erbium Inorganic materials 0.000 description 2
- 229910052769 Ytterbium Inorganic materials 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 238000010253 intravenous injection Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 239000003153 chemical reaction reagent Substances 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 150000004820 halides Chemical class 0.000 description 1
- 238000003331 infrared imaging Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000004020 luminiscence type Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- -1 rare earth halide Chemical class 0.000 description 1
- 150000002910 rare earth metals Chemical class 0.000 description 1
- 230000002269 spontaneous effect Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
- A61B6/035—Mechanical aspects of CT
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/10—Safety means specially adapted therefor
- A61B6/107—Protection against radiation, e.g. shielding
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/44—Constructional features of apparatus for radiation diagnosis
- A61B6/4417—Constructional features of apparatus for radiation diagnosis related to combined acquisition of different diagnostic modalities
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/48—Diagnostic techniques
- A61B6/482—Diagnostic techniques involving multiple energy imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/508—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for non-human patients
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5229—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
- A61B6/5247—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from an ionising-radiation diagnostic technique and a non-ionising radiation diagnostic technique, e.g. X-ray and ultrasound
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Heart & Thoracic Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Pulmonology (AREA)
- Theoretical Computer Science (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种深部组织X射线激发多光谱断层成像系统及方法，包括X射线源、X射线平板探测器、EMCCD相机、电控旋转台、窄带滤波片、铅板、计算机和待成像动物；利用X射线激发纳米发光材料，发射出特定谱段特别是具有较高穿透深度的短波红外光，光子穿过组织到达成像物体表面，经过窄带滤波片的过滤，不同波长谱段的光被EMCCD相机接收，得到每个谱段下的成像和重建结果；采用多变量分析对多谱成像和重建结果进行处理分析，得到纳米发光材料在物体内的准确分布，并与CT断层成像结果融合最终得到深部组织多模态的断层成像图像。本发明可有效提高X射线激发光学断层成像的成像深度，实现活体小动物的深部组织成像。
Description
技术领域
本发明属于多模态分子成像技术领域，尤其涉及一种深部组织X射线激发多光谱断层成像的方法及系统。
背景技术
光学分子成像技术能够在体、非接触、动态地对生物体内细胞分子水平的变化进行观测，为生物医学研究领域若干重要问题的解决提供了有效技术手段。但光在生物组织中的强散射特性极大限制了光学成像的穿透深度及成像分辨率，也限制了光学成像技术在深部组织器官中的应用。
X射线发光断层成像(x-ray luminescence computed tomography,XLCT)是一种新型的融合了CT断层成像和光学分子影像的多模态分子成像技术，于2010年由斯坦福大学的Xing Lei所在的研究小组最早提出。某些纳米发光材料在用X射线激发之后能够发出不同波长的光，激发出的光穿过组织达到成像物体表面被高灵敏度的EMCCD(ElectronMultiplying Charge Coupled Device)相机捕获，经过光学重建最终恢复出成像物体内部的纳米发光材料分布。光学成像与X射线扫描得到的结构成像的融合，最终得到多模态成像。
虽然X射线良好的穿透性以及背景光学信号的消除使得X射线发光断层成像有望实现活体小动物深部组织的断层成像，但是由于目前X射线激发光的发射谱段大部分在可见光以及近红外谱段，穿透深度仍有限。目前XLCT相关专利有CN103110425A，CN103876770A，仍处于公布期，其专利主要针对XLCT的系统结构和成像模式，没有明确深部成像所需材料、发光谱段和针对多光谱数据的分析方法，导致难以实现活体小动物深部组织的成像。短波红外成像(short wavelength infrared imaging,SWIR)利用第二红外窗(1000-2300nm)对组织进行成像，相较传统的可见光或近红外光，SWIR在组织中的光子吸收、散射及自发光现象大大减少，从而在改进空间分辨率和成像质量的同时，能在生物组织中穿透更深。同时，SWIR波段可有效避开生物组织自发光的影响，结合纳米发光材料的研究，基于SWIR谱段的X射线激发发光成像可实现深部组织的断层成像。
发明内容
本发明的目的在于提供一种用于深部组织X射线激发多光谱断层成像的系统及方法，旨在解决目前X射线激发光的发射谱段大部分在可见光以及近红外谱段，XLCT的成像深度有限，重建结果准确性低，难以实现活体小动物深部组织的成像的问题。目前尚未见采用具有更好穿透性的短波红外光XLCT成像相关专利，亦未见与之相配的能够被X射线激发的纳米发光材料和针对深部组织成像的多光谱数据的分析方法。
本发明实现一种深部组织X射线激发多光谱断层成像的系统，该系统包括1、X射线源 2、X射线平板探测器 3、EMCCD相机 4、电控旋转台 5、窄带滤波片 6、计算机 7、铅板 8、成像物体。窄带滤波片置于EMCCD相机的镜头前，二者位于X射线源，成像物体和X射线平板探测器连线的垂直方向；X射线源、X射线平板探测器、EMCCD相机、电控旋转台均与计算机相连，使用计算机控制其运行。
进一步，所述X射线源用来发射X射线。
进一步，所述X射线平板探测器用来接收探测X射线。
进一步，所述EMCCD相机，用来接收短波红外光及其它特定谱段光。
进一步，所述电控旋转台，用来固定支持成像物体并进行旋转以采集多个角度下的X射线投影数据和光学投影数据。
进一步，所述窄带滤波片，用来滤除标示带宽之外谱段的光。
进一步，所述计算机，用来控制X射线源、X射线平板探测器、EMCCD相机、电控旋转台的运行以及参数设置，控制光学数据和X射线投影数据的采集和重建，融合双模态的成像。
进一步，所述铅板，用来滤除多余的X射线以减少射线对EMCCD芯片的损伤。
本发明的另一目的在于提供一种深部组织X射线激发多光谱断层成像的方法，该方法利用X射线激发注射到成像物体体内的特定纳米发光材料，这种材料在X射线激发下能够发出短波红外和其他特定谱段的光；光穿过成像物体到达表面，经过窄带滤波片的过滤，不同波长谱段的光被EMCCD相机接收；在每个谱段下采用迭代重建算法对EMCCD接收到的激发光数据进行重建得到每个谱段下的成像和重建结果；采用多变量分析方法中的PCA分解方法对多个谱段下的成像和重建结果进行处理分析，得到纳米发光材料在成像物体内的准确分布，作为光学断层成像结果；穿过成像物体的X射线被X射线探测器接收保存为X射线投影；通过经典的FDK重建算法对X射线投影进行重建得到CT断层成像结果；对CT断层成像结果与光学断层成像结果根据坐标关系进行配准并融合最终得到深部组织多模态的断层成像图像，该用于深部组织X射线激发多光谱断层成像的方法包括以下步骤：
步骤一，准备成像物体：将成像物体置于旋转台上，向成像物体内部加入X射线激发的纳米发光材料。
步骤二、激发光数据采集：利用X射线激发成像物体深部组织内的激发纳米发光材料产生短波红外及其它特定谱段光，穿过组织的光学信号经过窄带滤波片的过滤被EMCCD接收，每隔一定角度(如15度)采集一副激发光数据并保存，直至成像物体完成360度的旋转；
步骤三、白光数据采集：关闭X射线源，打开外界光源，取下窄带滤波片，每隔一定角度(如5度)采集一幅成像物体的白光图像并保存，直至完成360度的旋转。
步骤四、X射线投影数据的采集：打开X射线源发射X射线，利用X射线平板探测器采集穿过成像物体的X射线，每隔一定角度(如1度)采集一幅投影数据并保存，直至成像物体完成180度或360度的旋转。
步骤五、光学三维断层重建：对步骤一获得的激发光图像进行预处理，采用二维滤波器中的自适应中值滤波去除图像中的孤立亮点、坏点和部分噪点；每个谱段下处理图像，利用步骤三中的白光数据构建物体轮廓并剖分，将各个角度下的光强映射到成像物体轮廓表面，采用迭代重建算法得到该谱段下的断层重建结果。
步骤六、多光谱数据分析：对所有谱段下的激发光图像(步骤二)和重建结果(步骤五)进行多变量分析，利用多光谱信息作为先验，获得短波红外谱段的更准确成像，从而得到纳米发光材料在成像物体深部组织内的准确分布。
步骤七、CT三维断层重建：对步骤四采集到的X射线投影数据进行预处理，包括坏点校正，暗场校正以及亮场校正等，采用经典的FDK重建算法对处理后投影数据进行重建，得到小成像物体的三维结构成像。
步骤八、双模态重建图像的配准融合：对步骤六和步骤七得到的短波红外光学断层图像和CT断层图像根据坐标关系进行配准并融合，实现双模态成像。
进一步，所述步骤二根据本发明的特点和要求采用能够在X射线激发下发射短波红外及其它特定谱段光并具有良好的生物兼容性的纳米发光材料(如稀土卤化物发光材料LaBr3:Ce3+，NaYF4:Er,Yb等)通过静脉注射到成像物体体内。
进一步，所述步骤二关闭外界灯光光源，打开X射线源发射X射线；利用X射线激发成像物体深部组织内的激发纳米发光材料产生短波红外光及其它特定谱段光，穿过组织的光学信号经过窄带滤波片的过滤被EMCCD相机接收，形成激发光数据；计算机控制电控旋转台每隔15度采集一副激发光数据并保存，直至成像物体完成360度的旋转，更换窄带滤波片，重复上述步骤；
进一步，所述步骤四打开X射线源发射X射线，利用X射线平板探测器采集穿过成像物体的X射线，形成一幅投影数据并保存；计算机控制电控旋转台每隔1度采集一副X射线投影数据并保存，直至成像物体完成360度的旋转。
进一步，所述步骤五采用传输方程描述光子在生物组织中的传输过程。对于一个谱段下的激发光数据，将每个角度下的激发光图像映射到物体表面，根据传输方程，利用有限元方法建立表面节点上的光强与内部纳米发光材料分布的关系，采用迭代重建方法恢复出内部纳米发光材料分布，完成一个谱段下的断层重建。
进一步，所述步骤六对所有谱段下的激发光图像和重建结果进行多变量分析，得到纳米发光材料在成像物体深部组织内的分布情况的操作方法为：对于每个谱段下的重建图像Xs，其矩阵大小为N＝W×L×H,W为图像的宽度，L为图像的长度，H为图像的高度。将Xs装入新矩阵X＝{X1,X2,...,Xs,...,XM}作为其中一列，其中M为谱段数目，X为主成分分析方法的输入，大小为N×M。
定义矩阵其中X0为矩阵X每一列减去该列的均值形成的矩阵，通过对角化计算得到矩阵P的特征向量为T＝{T1,T2,...,Tk,...,TM},主成分分析方法的输出为
Q＝X0×T
第k个分量为
Qk＝X0×Tk
输出的各个分量Qk即代表纳米发光材料在不同的深部组织和器官中的分布。
本发明的X射线源发出X射线，X射线与位于成像物体深部组织内部的激发纳米发光材料发生作用，发出具有较好穿透深度的短波红外及特定谱段光，可采用高光产额的LaBr3:Ce3+为代表的卤化物发光材料，或选用文献”X-ray-Induced Shortwave InfraredBiomedical Imaging Using RareEarth Nanoprobes”中提及的已经用于小动物深部组织的稀土掺杂的NaYF4:Er,Yb，这些材料其在X射线激发下能发出峰值≥1500nm的短波红外及其它特定谱段光，穿透深度可达1-2cm，能够满足活体小动物的深部组织成像要求。光穿过组织到达成像物体表面，经过窄带滤波片的过滤，被EMCCD相机接收，得到不同谱段的光学数据保存在计算机，穿过成像物体的X射线被X射线平板探测器探测接收形成X射线投影数据并保存在计算机，由计算机重建得到每个谱段下的光学断层图像和CT断层成像，采用多变量分析方法准确得到纳米发光材料在物体内的分布，与CT断层成像的重建结果融合最终得到多模态的图像。本发明采用可被X射线激发发射出短波红外及其它特定谱段光的纳米发光材料，提高了激发光的穿透深度。采用多变量分析方法，能有效的定位深部组织内的激发纳米发光材料位置，实现活体小动物深部组织的XLCT成像。
附图说明
图1是本发明实施例提供的用于深部组织X射线激发多光谱断层成像的方法流程图。
图2是本发明实施例提供的用于深部组织X射线激发多光谱断层成像的系统结构示意图。
图中：1、X射线源；2、X射线平板探测器 3、EMCCD相机 4、电控旋转台 5、窄带滤波片 6、计算机 7、铅板 8、待成像动物。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
下面结合附图对本发明的应用原理作进一步描述。
如图1所示：一种深部组织X射线激发多光谱断层成像的方法，利用X射线激发纳米发光材料，发射出特定谱段特别是具有较高穿透深度的短波红外光，光子穿过组织到达成像物体表面，经过窄带滤波片的过滤，不同波长谱段的光被EMCCD相机接收，得到每个谱段下的成像和重建结果。采用多变量分析对多谱成像和重建结果进行处理分析，得到纳米发光材料在物体内的准确分布，并与CT断层成像结果融合最终得到深部组织多模态的断层成像图像，该用于深部组织X射线激发多光谱断层成像的方法包括以下步骤：
S101：纳米发光试剂注射：麻醉待成像动物，向其体内注射浓度为1μg/ml的可被X射线激发的纳米发光材料LaBr3:Ce3+。
S102：激发光数据采集：设置X射线源的电压为50kV，电流为1mA,利用X射线激发小动物深部组织内的LaBr3:Ce3+产生波长大于1500nm的短波红外及其它特定谱段光，设置EMCCD的积分时间为1s，穿过组织的光学信号经过窄带滤波片的过滤被EMCCD接收，每隔15度采集一副激发光数据并保存，直至待成像动物完成360度的旋转。
S103:白光数据采集：关闭X射线源，打开外界光源，取下窄带滤波片，每隔5度采集一幅动物的白光图像并保存，直至完成360度的旋转。
S104：X射线投影数据的采集：打开X射线源,设置X射线源的电压为50kV，电流为1mA，积分时间为300ms，利用X射线平板探测器采集穿过动物的X射线，每隔1度采集一幅投影数据并保存，直至完成360度的旋转。
S105：光学三维断层重建：对步骤一获得的激发光图像进行预处理，采用二维滤波器中的自适应中值滤波器去除图像中的孤立亮点、坏点和部分噪点；在每个谱段下处理图像，利用步骤二中的白光数据构建动物轮廓并剖分，将各个角度下的光强映射到动物轮廓表面，采用迭代重建算法得到该谱段下的断层重建结果。
S106：多光谱数据分析：对所有谱段下的激发光图像(步骤二)和重建结果(步骤四)进行多变量分析，利用多光谱信息作为先验，获得短波红外谱段的更准确成像，从而得到纳米发光材料在深部组织内的准确分布。
S107：CT三维断层重建：对步骤三采集到的X射线投影数据进行预处理，包括坏点校正，暗场校正以及亮场校正等，采用经典的FDK重建算法对处理后投影数据进行重建，得到动物的三维结构成像。
S108：双模态重建图像的配准融合：对步骤六和步骤七得到的短波红外光学断层图像和CT断层图像根据坐标关系进行配准并融合，实现双模态成像。
所述S101根据本发明的特点和要求采用能够在X射线激发下发射短波红外及其它特定谱段光并具有良好的生物兼容性的纳米发光材料-稀土卤化物发光材料LaBr3:Ce3+，通过静脉注射到小动物体内。
所述S102关闭外界灯光光源，打开X射线源发射X射线；利用X射线激发小动物深部组织内的激发纳米发光材料产生短波红外及特定谱段光，穿过组织的光学信号经过窄带滤波片的过滤被EMCCD相机接收，形成激发光数据；计算机控制电控旋转台每隔15度采集一副激发光数据并保存，直至小动物完成360度的旋转，更换窄带滤波片，重复上述步骤；
所述S104打开X射线源发射X射线，利用X射线平板探测器采集穿过小动物的X射线，形成一幅投影数据并保存；计算机控制电控旋转台每隔1度采集一副X射线投影数据并保存，直至小动物完成360度的旋转。
所述S105采用传输方程描述光子在生物组织中的传输过程。对于一个谱段下的激发光数据，将每个角度下的激发光图像映射到物体表面，根据传输方程，利用有限元方法建立表面节点上的光强与内部纳米发光材料分布的关系，采用迭代重建方法恢复出内部纳米发光材料分布，完成一个谱段下的断层重建。
所述S106对所有谱段下的激发光图像和重建结果进行多变量分析，得到纳米发光材料在成像物体深部组织内的分布情况的操作方法为：对于每个谱段下的重建图像Xs，其矩阵大小为N＝256×256×400,图像的宽度为256，图像的长度为256，图像的高度为400。将Xs装入新矩阵X＝{X1,X2,X3}作为其中一列，其中谱段数目为3，X为主成分分析方法的输入，大小为N×3。
定义矩阵其中X0为矩阵X每一列减去该列的均值形成的矩阵，通过对角化计算得到矩阵P的特征向量为T＝{T1,T2,T3},主成分分析方法的输出为
Q＝X0×T
第k个分量为
Qk＝X0×Tk
输出的各个分量Qk即代表纳米发光材料在不同的深部组织和器官中的分布。
如图2所示：一种深部组织X射线激发多光谱断层成像的系统，该用于深部组织X射线激发多光谱断层成像的系统包括X射线源1、X射线源 2、X射线平板探测器 3、EMCCD相机4、电控旋转台 5、窄带滤波片 6、计算机 7、铅板 8、待成像动物。窄带滤波片置于EMCCD相机的镜头前，二者位于X射线源，待成像活体小动物和X射线平板探测器连线的垂直方向；X射线源、X射线平板探测器、EMCCD相机、电控旋转台均与计算机相连，使用计算机控制其运行。
进一步，所述X射线源用来发射X射线。
进一步，所述X射线平板探测器用来接收探测X射线。
进一步，所述EMCCD相机，用来接收短波红外及其它特定谱段光。
进一步，所述电控旋转台，用来固定支持小动物并进行旋转以采集多个角度下的X射线投影数据和光学投影数据。
进一步，所述窄带滤波片，用来滤除标示带宽之外谱段的光。进一步，所述计算机，用来控制X射线源、X射线按测器、EMCCD相机、电控旋转台的运行以及参数设置，控制光学数据和X射线投影数据的采集和重建，融合双模态的成像。
进一步，所述铅板，用来滤除多余的X射线以减少射线对EMCCD芯片的损伤。
X射线源1发出X射线，X射线与位于活体小动物8深部组织内部的激发纳米发光材料LaBr3:Ce3+发生作用，发出具有较好穿透深度的短波红外及其它特定谱段光，其在X射线激发下能发出峰值≥1500nm的短波红外及其它特定谱段光，穿透深度可达1-2cm，能够满足活体小动物的深部组织成像要求。光穿过组织到达小动物表面，经过窄带滤波片5的过滤，被EMCCD相机3接收，得到不同谱段的光学数据保存在计算机6，穿过小动物的X射线被X射线平板探测器2探测接收形成X射线投影数据并保存在计算机，由计算机重建得到每个谱段下的光学断层图像和CT断层成像，采用多变量分析方法准确得到纳米发光材料在物体内的分布，与CT断层成像的重建结果融合最终得到多模态的图像。
采用可被X射线激发发射出短波红外及其它特定谱段光的纳米发光材料LaBr3:Ce3 +，提高了激发光的穿透深度。采用多变量分析方法，能有效的定位深部组织内的激发纳米发光材料位置，实现活体小动物深部组织的XLCT成像。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (5)
1.一种深部组织X射线激发多光谱断层成像方法，其特征在于，所述深部组织X射线激发多光谱断层成像方法利用X射线激发成像物体体内的特定纳米发光材料，这种材料在X射线激发下能够发出短波红外和其他特定谱段的光；深部组织内的光子穿过组织到达成像物体表面，经过窄带滤波片的过滤，不同波长谱段的光被EMCCD相机接收；根据光在组织中的传输方程，建立发光材料的分布与表面光强分布的关系，在每个谱段下采用迭代重建算法对EMCCD接收到的激发光数据进行重建，得到每个谱段下的成像和重建结果；采用多变量分析方法中的主成分分析方法对多谱成像和重建结果进行处理分析，得到纳米发光材料在成像物体内的准确分布，作为光学断层成像结果；穿过成像物体的X射线被X射线探测器接收保存为X射线投影；通过经典的FDK重建算法对X射线投影进行重建，得到CT断层成像结果；对CT断层成像结果与光学断层成像结果根据坐标关系进行配准并融合最终得到深部组织多模态的断层成像图像；
所述深部组织X射线激发多光谱断层成像方法使用的深部组织X射线激发多光谱断层成像系统包括：
X射线源，与计算机连接，用于发射X射线；
X射线平板探测器，与计算机连接，用于接收探测X射线；
EMCCD相机，与计算机连接，用于接收短波红外光及其它特定谱段发光；
电控旋转台，与计算机连接，用于固定并支持成像物体进行高精度旋转，以采集多个角度下的X射线投影数据和激发光数据；
窄带滤波片，置于EMCCD相机的镜头前，用于滤除标示带宽之外谱段的光；
铅板，用来滤除多余的X射线，并保护EMCCD相机免于X射线辐射；
计算机，用于控制X射线源、X射线平板探测器、EMCCD相机、电控旋转台的运行以及参数设置，控制光学数据和X射线投影数据的采集，进行多光谱数据和X射线投影数据的处理和重建，融合双模态的成像；
所述窄带滤波片和EMCCD相机位于X射线源和X射线平板探测器连线的垂直方向。
2.如权利要求1所述的深部组织X射线激发多光谱断层成像方法，其特征在于，所述深部组织X射线激发多光谱断层成像方法包括以下步骤：
步骤一，准备成像物体：将内部加入X射线激发的纳米发光材料的成像物体置于旋转台上；
步骤二，激发光数据采集：利用X射线激发成像物体内部纳米发光材料产生短波红外光及其它特定谱段光，穿过组织的光学信号经过窄带滤波片的过滤被EMCCD相机接收，每隔15度采集一幅激发光图像并保存，直至完成360度的旋转；
步骤三，白光数据采集：关闭X射线源，打开外界光源，取下窄带滤波片，每隔5度采集一幅成像物体的白光图像并保存，直至完成360度的旋转；
步骤四，X射线投影数据的采集：打开X射线源发射X射线，利用X射线平板探测器采集穿过成像物体的X射线，每隔1度采集一幅投影数据并保存，直至成像物体完成180度或360度的旋转；
步骤五，光学三维断层重建：对步骤二获得的激发光图像进行预处理，采用二维滤波器中的自适应中值滤波去除图像中的孤立亮点、坏点和部分噪点；在每个谱段下处理图像，利用步骤三中的白光数据构建成像物体轮廓并剖分，将各个角度下的光强映射到成像物体轮廓表面，采用迭代重建算法得到该谱段下的断层重建结果；
步骤六，多光谱数据分析：对所有谱段下的激发光图像和重建结果进行多变量分析，利用多光谱信息作为先验，获得短波红外谱段的更准确成像，从而得到纳米发光材料在成像物体深部组织内的分布；
步骤七，CT三维断层重建：对步骤四采集到的X射线投影数据进行预处理，包括坏点校正，暗场校正以及亮场校正，采用经典的FDK重建算法对处理后投影数据进行重建，得到成像物体的三维结构成像；
步骤八，双模态重建图像的配准融合：对步骤六和步骤七得到的短波红外光学断层图像和CT断层图像根据坐标关系进行配准并融合，实现双模态成像。
3.如权利要求2所述的深部组织X射线激发多光谱断层成像方法，其特征在于，所述步骤二采用能够在X射线激发下发射短波红外及特定谱段光并具有良好的生物兼容性的纳米发光材料，关闭外界灯光光源，打开X射线源发射X射线，利用X射线激发成像物体内的纳米发光材料，产生短波红外及其它特定谱段光；穿过成像物体的光学信号经过窄带滤波片的过滤被EMCCD相机接收，形成多光谱数据；计算机控制电控旋转台每隔15度采集一幅激发光数据并保存，直至成像物体完成360度的旋转，更换窄带滤波片，重复上述步骤。
4.如权利要求2所述的深部组织X射线激发多光谱断层成像方法，其特征在于，所述步骤四X射线源发射X射线，利用X射线平板探测器采集穿过成像物体的X射线，形成一幅投影数据并保存；计算机控制电控旋转台每隔1度采集一幅X射线投影数据并保存，直至成像物体完成360度的旋转。
5.如权利要求2所述的深部组织X射线激发多光谱断层成像方法，其特征在于，所述步骤六对所有谱段下的激发光图像和重建结果进行多变量分析，得到纳米发光材料在成像物体深部组织内的分布情况的操作方法为：对于每个谱段下的重建图像Xs，将Xs组装入新矩阵X＝{X1,X2,...,Xs,...,XM}，矩阵X为主成分分析方法的输入，输出的分解之后的图像矩阵即代表纳米发光材料在不同的深部组织和器官中的分布。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610118383.1A CN105640582B (zh) | 2016-03-02 | 2016-03-02 | 一种深部组织x射线激发多光谱断层成像系统及方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610118383.1A CN105640582B (zh) | 2016-03-02 | 2016-03-02 | 一种深部组织x射线激发多光谱断层成像系统及方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105640582A CN105640582A (zh) | 2016-06-08 |
| CN105640582B true CN105640582B (zh) | 2019-05-07 |
Family
ID=56492961
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610118383.1A Active CN105640582B (zh) | 2016-03-02 | 2016-03-02 | 一种深部组织x射线激发多光谱断层成像系统及方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105640582B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2021168691A1 (en) * | 2020-02-26 | 2021-09-02 | Shenzhen Xpectvision Technology Co., Ltd. | Imaging system using x-ray fluorescence |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106137129B (zh) * | 2016-06-28 | 2020-03-06 | 中国科学院深圳先进技术研究院 | 荧光散射光学断层成像系统及方法 |
| WO2018000186A1 (zh) * | 2016-06-28 | 2018-01-04 | 中国科学院深圳先进技术研究院 | 荧光散射光学断层成像系统及方法 |
| CN106955121B (zh) * | 2017-03-14 | 2020-06-23 | 杭州电子科技大学 | X射线调制的x射线激发发光成像方法 |
| CN106994001B (zh) * | 2017-03-30 | 2023-08-18 | 展谱光电科技（上海）有限公司 | 便携式多光谱摄像装置及其使用方法 |
| CN106959571B (zh) * | 2017-04-07 | 2022-07-08 | 展谱光电科技（上海）有限公司 | 多光谱投影及摄像装置以及多光谱投影方法 |
| CN107411706B (zh) * | 2017-04-07 | 2023-08-18 | 展谱光电科技（上海）有限公司 | 多光谱投影装置和投影方法、多光谱摄像装置和摄像方法 |
| CN107860781A (zh) * | 2017-10-12 | 2018-03-30 | 兰州职业技术学院 | 一种食品微生物定性与定量的检测方法 |
| CN107703165B (zh) * | 2017-10-23 | 2023-08-08 | 中国科学院苏州生物医学工程技术研究所 | 具有精密样品检测转台的ct成像设备 |
| CN109147049B (zh) * | 2018-07-28 | 2022-11-01 | 天津大学 | 用于x射线光动力治疗的影像重建方法 |
| CN109589126B (zh) * | 2018-10-18 | 2022-08-19 | 天津大学 | 基于宽光束小步长扫描方式的x射线发光断层成像方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102429668A (zh) * | 2011-07-28 | 2012-05-02 | 中国科学院自动化研究所 | 一种功能-分子-结构成像系统及方法 |
| CN103110425A (zh) * | 2013-03-11 | 2013-05-22 | 西安电子科技大学 | 一种联合成像的锥束x射线发光断层成像系统及其方法 |
| CN103462623A (zh) * | 2013-09-04 | 2013-12-25 | 江苏美伦影像系统有限公司 | 一种快速采集图像的x射线数字成像系统 |
| CN103810755A (zh) * | 2014-03-04 | 2014-05-21 | 西安电子科技大学 | 基于结构聚类稀疏表示的压缩感知光谱图像重建方法 |
| CN103876770A (zh) * | 2014-03-12 | 2014-06-25 | 西安电子科技大学 | 一种窄束x射线激发的小动物多光谱荧光断层成像系统 |
| CN204649311U (zh) * | 2012-09-14 | 2015-09-16 | 菲力尔系统公司 | 测量装置 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2011133606A2 (en) * | 2010-04-19 | 2011-10-27 | The Regents Of The University Of California | Real-time volumetric image reconstruction and 3d tumor localization based on a single x-ray projection image for lung cancer radiotherapy |
-
2016
- 2016-03-02 CN CN201610118383.1A patent/CN105640582B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102429668A (zh) * | 2011-07-28 | 2012-05-02 | 中国科学院自动化研究所 | 一种功能-分子-结构成像系统及方法 |
| CN204649311U (zh) * | 2012-09-14 | 2015-09-16 | 菲力尔系统公司 | 测量装置 |
| CN103110425A (zh) * | 2013-03-11 | 2013-05-22 | 西安电子科技大学 | 一种联合成像的锥束x射线发光断层成像系统及其方法 |
| CN103462623A (zh) * | 2013-09-04 | 2013-12-25 | 江苏美伦影像系统有限公司 | 一种快速采集图像的x射线数字成像系统 |
| CN103810755A (zh) * | 2014-03-04 | 2014-05-21 | 西安电子科技大学 | 基于结构聚类稀疏表示的压缩感知光谱图像重建方法 |
| CN103876770A (zh) * | 2014-03-12 | 2014-06-25 | 西安电子科技大学 | 一种窄束x射线激发的小动物多光谱荧光断层成像系统 |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2021168691A1 (en) * | 2020-02-26 | 2021-09-02 | Shenzhen Xpectvision Technology Co., Ltd. | Imaging system using x-ray fluorescence |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105640582A (zh) | 2016-06-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105640582B (zh) | 一种深部组织x射线激发多光谱断层成像系统及方法 | |
| CN109938764B (zh) | 一种基于深度学习的自适应多部位扫描成像方法及其系统 | |
| Schulz et al. | Experimental fluorescence tomography of tissues with noncontact measurements | |
| JP2022016591A (ja) | Ｘ線コンピュータ断層撮影装置、スキャン計画設定支援装置、医用画像診断システム、制御方法及び制御プログラム | |
| CN100593389C (zh) | 一种连续动态采集式小动物诱发荧光分子成像系统 | |
| US9635349B2 (en) | Second generation hand held optical imager | |
| CN103070673A (zh) | 一种在体小动物荧光分子断层成像系统与方法 | |
| Shao et al. | Designing a wearable navigation system for image-guided cancer resection surgery | |
| CN106651982B (zh) | 一种基于阵列x射线源和探测器的ct图像重建方法 | |
| CN104168830B (zh) | 内窥式x射线发光断层成像装置及方法 | |
| US8712504B2 (en) | Hand-held optical probe based imaging system with 3D tracking facilities | |
| CN114128243A (zh) | 在缺光环境中利用拓扑激光扫描进行超光谱和荧光成像 | |
| CN105832357A (zh) | 使用多光谱x射线源的全景成像 | |
| ES2800032T3 (es) | Sistema de obtención de imágenes médicas con una red fija de detectores de rayos X y método para usarlo | |
| CN103110425A (zh) | 一种联合成像的锥束x射线发光断层成像系统及其方法 | |
| JP2018515753A (ja) | ハイパースペクトル・イメージングのためのシステムおよび方法 | |
| CN109044277B (zh) | 近红外二区荧光断层成像系统 | |
| CN114245868A (zh) | 使用用于超光谱和荧光成像和形貌激光标测的单色图像传感器的宽动态范围 | |
| CN103239255A (zh) | 一种锥束x射线发光断层成像方法 | |
| CN105796121B (zh) | 一种ct和x射线激发荧光双模同步断层成像方法 | |
| CN106419837A (zh) | 基于螺旋式激发的荧光分子断层成像系统 | |
| EP3531169A1 (en) | Dynamically scanned x-ray detector panel | |
| CN107249465B (zh) | 用于稀疏角度采样的断层摄影成像装置和方法 | |
| CN114025656A (zh) | 脉冲超光谱、荧光和激光标测成像系统中的噪声感知边缘增强 | |
| CN103876770B (zh) | 一种窄束x射线激发的小动物多光谱荧光断层成像系统 |
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 |