CN113171116A - 一种用于检测软组织异物的多模态x射线成像设备及其方法 - Google Patents
一种用于检测软组织异物的多模态x射线成像设备及其方法 Download PDFInfo
- Publication number
- CN113171116A CN113171116A CN202110391180.0A CN202110391180A CN113171116A CN 113171116 A CN113171116 A CN 113171116A CN 202110391180 A CN202110391180 A CN 202110391180A CN 113171116 A CN113171116 A CN 113171116A
- Authority
- CN
- China
- Prior art keywords
- image
- sample
- ray
- grating
- contrast image
- 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.)
- Pending
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 76
- 210000004872 soft tissue Anatomy 0.000 title claims abstract description 29
- 238000000034 method Methods 0.000 title claims abstract description 10
- 238000010521 absorption reaction Methods 0.000 claims abstract description 58
- 238000001514 detection method Methods 0.000 claims abstract description 22
- 238000004458 analytical method Methods 0.000 claims abstract description 16
- 238000012545 processing Methods 0.000 claims abstract description 9
- 230000033228 biological regulation Effects 0.000 claims description 7
- 230000003287 optical effect Effects 0.000 claims description 5
- 238000012216 screening Methods 0.000 claims description 5
- 238000007405 data analysis Methods 0.000 claims description 4
- WFKWXMTUELFFGS-UHFFFAOYSA-N tungsten Chemical compound [W] WFKWXMTUELFFGS-UHFFFAOYSA-N 0.000 claims description 4
- 229910052721 tungsten Inorganic materials 0.000 claims description 4
- 239000010937 tungsten Substances 0.000 claims description 4
- 239000003086 colorant Substances 0.000 claims description 3
- 238000003672 processing method Methods 0.000 claims description 3
- 238000006073 displacement reaction Methods 0.000 claims description 2
- 238000007781 pre-processing Methods 0.000 claims description 2
- 230000003750 conditioning effect Effects 0.000 claims 1
- 238000011282 treatment Methods 0.000 abstract description 8
- 238000005516 engineering process Methods 0.000 abstract description 7
- 230000003902 lesion Effects 0.000 abstract description 6
- 239000000126 substance Substances 0.000 abstract description 6
- 238000007920 subcutaneous administration Methods 0.000 abstract description 4
- 210000003205 muscle Anatomy 0.000 description 21
- 239000004033 plastic Substances 0.000 description 15
- 229920003023 plastic Polymers 0.000 description 15
- 210000000577 adipose tissue Anatomy 0.000 description 12
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 11
- 239000010439 graphite Substances 0.000 description 11
- 229910002804 graphite Inorganic materials 0.000 description 11
- 239000000919 ceramic Substances 0.000 description 10
- 229910052751 metal Inorganic materials 0.000 description 10
- 239000002184 metal Substances 0.000 description 10
- 239000011521 glass Substances 0.000 description 9
- 210000001519 tissue Anatomy 0.000 description 9
- 239000002023 wood Substances 0.000 description 9
- 230000005540 biological transmission Effects 0.000 description 6
- 238000003745 diagnosis Methods 0.000 description 6
- 241000287828 Gallus gallus Species 0.000 description 5
- 238000000338 in vitro Methods 0.000 description 5
- 201000010099 disease Diseases 0.000 description 4
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 4
- 239000000463 material Substances 0.000 description 4
- 239000002985 plastic film Substances 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000000737 periodic effect Effects 0.000 description 3
- 206010006187 Breast cancer Diseases 0.000 description 2
- 208000026310 Breast neoplasm Diseases 0.000 description 2
- 210000000988 bone and bone Anatomy 0.000 description 2
- 210000000481 breast Anatomy 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000007918 intramuscular administration Methods 0.000 description 2
- 230000001788 irregular Effects 0.000 description 2
- 230000035515 penetration Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 235000017166 Bambusa arundinacea Nutrition 0.000 description 1
- 235000017491 Bambusa tulda Nutrition 0.000 description 1
- 241001330002 Bambuseae Species 0.000 description 1
- 208000004434 Calcinosis Diseases 0.000 description 1
- 206010014561 Emphysema Diseases 0.000 description 1
- 235000015334 Phyllostachys viridis Nutrition 0.000 description 1
- 239000004698 Polyethylene Substances 0.000 description 1
- 208000016247 Soft tissue disease Diseases 0.000 description 1
- 208000026062 Tissue disease Diseases 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 206010003246 arthritis Diseases 0.000 description 1
- 210000001188 articular cartilage Anatomy 0.000 description 1
- 239000011425 bamboo Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 201000003149 breast fibroadenoma Diseases 0.000 description 1
- 230000002308 calcification Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 210000000845 cartilage Anatomy 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000007850 degeneration Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 206010020718 hyperplasia Diseases 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 210000005075 mammary gland Anatomy 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 229920006255 plastic film Polymers 0.000 description 1
- -1 polyethylene Polymers 0.000 description 1
- 229920000573 polyethylene Polymers 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/12—Devices for detecting or locating foreign bodies
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. 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/5235—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 the same or different ionising radiation imaging techniques, e.g. PET and CT
Abstract
本发明涉及生物医疗领域及医疗设备技术，公开了一种用于检测软组织异物的多模态X射线成像设备及其方法。成像检测方法包括以下步骤：工作状态调节步骤、采集样本步骤、信号分析步骤、图像增强步骤等。成像设备包括工作状态调节系统、数据采集系统和计算机数据处理系统等。该方法可以获得病变的三个方面的信息：样品的吸收能力、电子密度分布和小角度散射能力，并且通过对样品的吸收能力和小角度散射能力两方面信息的整合和计算，我们可以获得具有一定单角度物质鉴别能力的R值图像，从而使常见异物的检测更加全面。在不同的皮下层次，不同的对比结果对不同类型异物的鉴别有不同的优势。
Description
技术领域
本发明涉及生物医疗领域及医疗设备技术，具体涉及一种用于检测软组织异物的多模态X射线成像设备及其方法。
背景技术
X射线透射成像及其相关的X射线计算机断层成像(CT)技术作为获取医学影像的基础性方法，与核磁共振(MRI)、正负电子湮灭断层成像(PET)等方法相比具有成本低、速度快、效率高的特点，在生物医学研究和疾病诊疗中应用广泛。
目前，医院普遍使用的X射线成像装置均基于点投影成像的原理，利用被检测对象内部元素组分及密度变化导致对X射线吸收能力的差异，获得具有一定衬度(或对比度)的图像。虽然现有X射线成像装置给出的吸收衬度图像包含生物体内部结构的丰富信息，能够有效指导研究和疾病诊疗，然而从实际效果来看，该技术仍存在以下两方面的不足。
首先，由于X射线强度的衰减随着物质原子序数和密度的降低而减小，因此在人体(或其他生物体)的X射线透射图像中，平均原子序数较高、密度较大的骨骼等由于其对X射线吸收较强可以清晰呈现，而由低Z元素组成且密度较低的软组织图像的衬度较低且外形轮廓模糊，制约了X射线成像技术在乳腺、关节(软骨)、肌肉等组织疾病诊疗中的应用。例如，乳腺癌，乳腺纤维腺瘤，乳腺增生，关节软骨退变，关节炎的诊断及鉴别都存在很大困难。
其次，现有X射线成像装置受到使用的毫米焦斑射线源和百微米像元尺寸电子记录装置(如平板探测器、闪烁体阵列等)的限制，所得的点投影图像的空间分辨率较低(0.2-0.5mm)，无法清晰反映一些尺度小但指示性强的病变特征的出现(如乳腺中的钙化点)或者具备多孔结构的组织器官(如肺、骨骼)的形貌信息，不利于相关疾病(如乳腺癌、肺气肿等)的早期发现和诊疗。
由于吸收衬度X射线透射成像技术存在上述不足，因此发展能够高对比度呈现生物体软组织异物形态和尺度病变特征的医用多模态X线成像仪，对于拓展X射线透射成像技术的适用范围、提升软组织疾病的早期发现及诊疗能力具有重要意义和应用价值。并且，鉴于X射线影像在临床医学中的基础性作用，医用多模态X线成像仪提供的生物组织结构及成分的多元信息将会对某些疾病的早期发现和诊断带来深远影响。
发明内容
为了解决上述现有技术存在的不足，本发明的首要目的在于提出一种用于检测软组织异物的多模态X射线成像设备的成像检测方法；
本发明的另一目的在于提供一种用于检测软组织异物的多模态X射线成像设备。
本发明是通过以下技术方案实现的：
一种用于检测软组织异物的多模态X射线成像设备的成像检测方法，使用基于光栅的多模态X射线成像检测设备包括：
产生X射线的光源；
采用X射线吸收光栅从而把从所述光源入射的X射线分离成多个线光源的源光栅；
能使X射线波前产生相位分布从而在光栅下游产生规则的竖状条纹的相位光栅；
采用X射线吸收光栅，配备步进位移控制器，将条纹图样的衰减、扭曲、或模糊信息转化为下游X射线强度信息的分析光栅；
接收X射线并转化为数字信号的X射线探测器；
其光路依次是X射线光源、源光栅、相位光栅、分析光栅和X射线探测器，X射线由X射线光源发出，依次经过光路上的各个组分；样品通常放置在相位光栅前方；
成像检测方法包括以下步骤：
工作状态调节步骤，开启X射线光源和X射线探测器，调节各部件至工作状态，并采集参考序列图像；即，对图像中的每一个像素都获得了无样品情况下的位置-光强曲线，即背景的步进曲线IB；
采集样本步骤，将样品放置于样品台上，校准后，移动样品或者对系统光栅部件进行步进调节，X射线探测器采集到样品序列图像；即，对图像中的每一个像素都获得了有样品情况下的位置-光强曲线，即样品的步进曲线IC；
信号分析步骤，经过校准和预处理之后，对参考序列图像和样品序列图像分别进行正弦拟合分析，然后进行比较和运算，进而获得关于样品的多模态图像；所述多模态图像包括吸收衬度图像imageA、相位衬度图像imagePS、暗场衬度图像imageS；
图像增强步骤，根据多模态图像的局域信息，可以对样品内部的感兴趣区域进行筛选区分，对每个像素的背景、样品进行处理，进而获得样品的R值图像。
数据处理方法包括以下步骤：
a、分别对每一个像素的背景的步进曲线IB、样品的步进曲线IC进行正弦拟合，获得步进曲线的直流分量kB和kC、幅值AB和AC、相角φB和φC；
b、计算-log(kC/kB)，获得吸收衬度图像imageA；
c、计算φC-φB，获得相位衬度图像imagePS；
d、计算-log((ACkB)/(ABkC))，获得暗场衬度图像imageS
e、计算imageA/imageS获得样品的R值图像。
软组织中低密度异物的检测是一个很大的挑战。本发明创新的提出基于光栅的多模态X射线成像，其具有硬件要求低的特点，利用样品的X射线吸收能力、电子密度分布和小角散射能力的差异，能同时获得吸收、相衬和暗场三种图像信息，在皮下软组织中低密度异物检测方面显示出巨大的效果。
本发明意外发现，吸收成像的信号值和暗场的信号值与厚度线性相关，它们的比值(定义为R)可以消除未知参数厚度值。将异物的吸收成像结果除以暗场结果，即可得到体外异物的R值。R值成像能一定程度上消去样品厚度因素，获得X射线路经上样品自身线性衰减系数和线性散射系数的比值，从而具备物质分辨能力。
基于光栅的多模态X射线成像可以一次获得互补的吸收衬度图像、相位衬度图像、暗场衬度图像和R值成像数据，可以清晰地将软组织异物区分开来。
R值计算方法：从源文件获取吸收成像和暗场成像结果中每个像素的信号强度，将每个像素点吸收信号强度值除以该点的暗场信号强度值，可获得每个像素点R值数值。将每个像素点R值数值作为灰度值重新绘制成图像，可得R值图像。
优选地，工作状态调节步骤和采集样本步骤中，在吸收衬度图像和暗场衬度图像中，选择目标像素线段，从源文件中获取该线段每个像素的信号强度，并根据位置序列绘制位置-光强曲线。
优选地，步骤d中将吸收衬度图像imageA、相位衬度图像imagePS、暗场衬度图像imageS分别转换成256级的红、蓝、绿三原色，然后进行图像合并，得到伪彩色图像。
通过数据采集系统光栅部件处于一系列不同特定位置时探测器记录到的图像，获得参考序列图像；通过数据采集系统光栅部件处于一系列不同特定位置时或者样品扫描移动时探测器记录到的样品透射图像，获得样品序列图像。
一种用于检测软组织异物的多模态X射线成像设备包括：工作状态调节系统、数据采集系统和计算机数据处理系统；
所述数据采集系统包括：产生X射线的光源、源光栅、相位光栅、分析光栅、X射线探测器，以及精密调节模块。
优选地，所述工作状态调节系统包含针对各光学硬件的三维角度和位移调节器件，以及样品的宏观调节器件；通过对调试结果或已有数据的回溯，以及对当前X射线探测器采集到的图像结果的分析和判断，调节各光学硬件的位置、角度和样品的相对位置，使得系统满足成像条件。
优选地，所述计算机数据处理系统包括调节控制模块、数据采集模块、数据分析模块、图像增强模块；
1)调节控制模块负责工作状态调节过程中关键数据的记录、调取和反馈，为工作状态调节系统提供依据；
2)数据采集模块负责存储参考序列图像和样品序列图像，并对其进行预处理；
3)数据分析模块负责获取样品的多模态图像，即吸收衬度图像imageA、相位衬度图像imagePS、暗场衬度图像imageS；
4)图像增强模块负责对感兴趣区域进行筛选区分，进而获得样品的R值图像。
优选地，X射线的光源为焦斑尺寸为小于1mm的临床钨靶X射线管，工作电压为小于200kVp。
进一步优选地，X射线的光源为焦斑尺寸为400μm的临床钨靶X射线管，工作电压为55kVp，管电流为8.5mA。
优选地，源光栅与相位光栅之间的距离为2156mm；相位光栅与分析仪光栅之间的距离为308mm。
与现有技术相比，本发明具有如下有益效果：
本发明基于光栅的多模态X射线成像与传统的吸收X射线成像相比，可以获得病变的三个方面的信息：样品的吸收能力、电子密度分布和小角度散射能力，从而使常见异物的检测更加全面。
在不同的皮下层次，不同的对比结果对不同类型异物的鉴别有显著的优势。
本发明多模态X射线成像设备能实现对金属、玻璃、木材、塑料、石墨、陶瓷等异物的检测，采用自制控制系统，用户操作简单，能够精确快速地进行目标部位的检测。使医生和研究者对病变部位进行及时准确的判断和处理，减少检测时间及费用，并降低并发症的发生。这将大大降低患者的医疗费用及国家医保支出，并且提高患者的生存质量，可满足医院医疗等多种应用场所的需求。
附图说明
图1为X射线光栅干涉仪的原理图；
图2为观察面示意图；
图3为材料展示图，其中A、金属、玻璃、木材、塑料、石墨、陶瓷异物；B.猪脂肪组织；C.鸡腿。
图4为基于光栅的体外常见异物X射线多模态成像结果及对应于横断面各位置的原始信号值曲线显示在线段上。A、常见异物照片；B.吸收衬度图像；C.相位衬度图像；D.暗场衬度图像；E.伪彩色图像。
图5脂肪组织和肌肉中异物的基于光栅多模态X射线成像结果。
其中，脂肪组织的吸收衬度图像(A)、相位衬度图像(B)、暗场衬度图像(C)、伪彩色图像(D)；肌肉组织的吸收衬度图像(E)、相位衬度图像(F)、暗场衬度图像(G)、伪彩色图像(H)。(图像右侧的不规则阴影是鸡腿末端组织的边界)(彩色箭头表示异物穿透的位置。蓝色箭头指向金属针的位置，同样，绿色箭头指向玻璃物体，黄色箭头指向木制物体，紫色箭头指向塑料物体，红色箭头指向石墨物体，浅蓝箭头指向陶瓷物体)。
图6为将吸收结果除以暗场结果得到的R值。
其中，A、不同组织和各种异物的R值柱状图；B、肌肉组织中异物的吸收衬度图像；C、肌肉组织中异物的暗场衬度图像；D、将肌肉内异物的R值图像，(红色箭头指向塑料板的位置)。
具体实施方式
下面给出实施例以对本发明进行具体的描述，但不限于此。实施例及应用例所用的原材料均为市购。
实施例1
如图1和图2所示，用于检测软组织异物的多模态X射线成像设备，利用pfx-X射线探测器的三个光栅来实现临床X射线成像。
用于检测软组织异物的多模态X射线成像设备包括：
产生X射线的光源；
采用X射线吸收光栅从而把从所述光源入射的X射线分离成多个线光源的源光栅；
能使X射线波前产生相位分布从而在光栅下游产生规则的竖状条纹的相位光栅；
采用X射线吸收光栅，配备高精度步进位移控制器，将条纹图样的衰减、扭曲、或模糊信息转化为下游X射线强度信息的分析光栅；
接收X射线并转化为数字信号的X射线探测器。
光路上依次是X射线光源、源光栅、相位光栅、分析光栅和探测器，X射线由X射线光源发出，依次经过光路上的各个组分。样品通常放置在相位光栅前方。
样品下游的相位光栅(G1)产生一个周期性的强度模式，该模式可由另一个分析光栅(G2)分析，因为该模式对于X射线探测器的直接观察来说过于密集。利用相位步进法对周期性光强分布进行分析，即平移其中一个光栅并在中间位置获取图像。所获得的位置-光强曲线，对应于作为位置函数的像素中的光强度，通常被称为步进曲线。
该步进曲线的可见度是X射线光栅干涉仪的关键参数。源光栅(G0)位于X射线光源附近，将X射线光斑分割成一系列线源以提高横向相干性。
其中，X射线光源是一种焦斑尺寸为400μm的临床钨靶X射线管，工作电压为55kVp，管电流为8.5mA。加速电压的选择是实现能见度和合理传输率之间的折衷。周期性强度模式分析采用7步相位步进法，每步曝光50秒。光栅之间的距离和光栅的参数如下表所示。
表1光栅间距
| 组件 | 距离 |
| 源光栅(G0)-相位光栅(G1) | 2156mm |
| 相位光栅(G1)–分析仪光栅(G2) | 308mm |
表2光栅参数
| G0 | G1 | G2 | |
| 材料 | Au | Au | Au |
| Period(μm) | 21 | 5.25 | 3 |
| Height of bars(μm) | 85 | 6.29 | 91 |
| Duty cycle | 0.52 | 0.54 | 0.53 |
如图3A所示，本实施例选用异物材料：金属、玻璃、木材、塑料、石墨、陶瓷异物。异物尺寸从左到右依次如下：(i)大小约5mm长的金属针；(ii)大小约为3*5*0.8mm的玻璃片；(iii)大小约为1*5*1mm的竹签；(iv)大小约为3*5*0.2mm的聚乙烯塑料片；(v)大小约为1*4*1mm的石墨笔尖；(vi)大小约为3*5*2mm的陶瓷片。
为了评价这些生活中常见的异物能否在人体软组织中得到有效识别，本实施例用大块猪脂肪组织(如图3B，模拟人体脂肪组织)和鸡腿(如图3C，模拟人体肌肉组织)进行离体组织实验。上述物体依次刺入猪脂肪组织和鸡腿。
将金属、玻璃、木材、塑料、石墨、陶瓷等异物分别刺入猪脂肪组织和鸡腿肌肉组织的软组织中，利用本实验室研制的用于检测软组织异物的多模态X射线成像设备对上述异物进行检测。
图形采集：将鸡腿和猪脂肪组织固定在固定架上，放入用于检测软组织异物的多模态X射线成像设备中进行成像。
成像检测方法包括以下步骤：
工作状态调节步骤，开启X射线光源和X射线探测器，调节各部件至工作状态，并采集参考序列图像；即，对图像中的每一个像素都获得了无样品情况下的位置-光强曲线，即背景的步进曲线IB；
采集样本步骤，将样品放置于样品台上，校准后，移动样品或者对系统光栅部件进行步进调节，X射线探测器采集到样品序列图像；即，对图像中的每一个像素都获得了有样品情况下的位置-光强曲线，即样品的步进曲线IC；
信号分析步骤，经过校准和预处理之后，对参考序列图像和样品序列图像分别进行正弦拟合分析，然后进行比较和运算，进而获得关于样品的多模态图像；所述多模态图像包括吸收衬度图像imageA、相位衬度图像imagePS、暗场衬度图像imageS；
图像增强步骤，根据多模态图像的局域信息，可以对样品内部的感兴趣区域进行筛选区分，对每个像素的背景、样品进行处理，进而获得样品的R值图像。
具体成像操作步骤包括：
1、打开光源、探测器；
2、样品移出光路，三块光栅移入光路；
3、调节三块光栅的位置和角度，使三块光栅的栅线彼此平行；
4、调节分析光栅步进至数个特定位置，并在分析光栅稳定时，用探测器采集对应该位置的图像；
5、对图像中的每一个像素都获得了无样品情况下的位置-光强曲线，即背景的步进曲线IB；
6、样品移入光路；
7、重复4；
8、对图像中的每一个像素都获得了有样品情况下的位置-光强曲线，即样品的步进曲线IC。
吸收和暗场图象的强度：
在吸收衬度图像和暗场衬度图像中，选择一段感兴趣的像素线段，从源文件中获取该线段每个像素的信号强度，并根据位置序列绘制距离-强度曲线。
R值：
本试验意外发现，吸收成像的信号值和暗场的信号值与厚度线性相关，它们的比值(定义为R)可以消除未知参数厚度值。将异物的吸收成像结果除以暗场结果，即可得到体外异物的R值。通过对脂肪和肌肉图像的处理，也可以得到脂肪和肌肉组织的R值。
数据处理方法包括以下步骤：
a、分别对每一个像素的背景的步进曲线IB、样品的步进曲线IC进行正弦拟合，获得步进曲线的直流分量kB和kC、幅值AB和AC、相角φB和φC；
b、计算-log(kC/kB)，获得吸收衬度图像imageA；
c、计算φC-φB，获得相位衬度图像imagePS；
d、计算-log((ACkB)/(ABkC))，获得暗场衬度图像imageS
e、计算imageA/imageS获得样品的R值图像
伪彩色图像：
将吸收衬度图像imageA、相位衬度图像imagePS、暗场衬度图像imageS分别转换成256级的红、蓝、绿三原色，然后进行图像合并，得到伪彩色图像。
检测结果分析：
图4为基于光栅的体外常见异物X射线多模态成像结果及对应于横断面各位置的原始信号值曲线显示在线段上。A、常见异物照片；B.吸收衬度图像；C.相位衬度图像；D.暗场衬度图像；E.伪彩色图像。
体外吸收、相衬和暗场数据清楚地显示了这些异物的形态和大小。通过提取结果线段各像素的信号强度值，绘制曲线。从曲线上可以看出，吸收成像信号的值在异物处较高，相位成像信号的值在异物边界处较高。即吸收成像直接显示物体，而相衬成像主要显示物体的边界。暗场结果可以直接显示金属、木材、石墨和陶瓷，但只能显示玻璃和塑料异物的边界。
图5脂肪组织和肌肉中异物的基于光栅多模态X射线成像结果。(1)脂肪组织的吸收衬度图像(A)、相位衬度图像(B)、暗场衬度图像(C)、伪彩色图像(D)；肌肉组织的吸收衬度图像(E)、相位衬度图像(F)、暗场衬度图像(G)、伪彩色图像(H)。(图像右侧的不规则阴影是鸡腿末端组织的边界)(彩色箭头表示异物穿透的位置。蓝色箭头指向金属针的位置，同样，绿色箭头指向玻璃物体，黄色箭头指向木制物体，紫色箭头指向塑料物体，红色箭头指向石墨物体，浅蓝箭头指向陶瓷物体)。
在脂肪组织中，吸收衬度图像相对清晰地显示除木棍以外的异物(图5A)。相位衬度图像主要显示异物的边界轮廓，仅能清楚地看到金属、玻璃、石墨和陶瓷，但同时脂肪组织中看不到木材和塑料碎片(图5B)。除塑料外，这些异物的暗场衬度图像清晰可见(图5C)。此外，木质异物在暗场衬度图像中表现出明显的高信号，其强度与吸收成像中的金属信号几乎相同。
在肌肉组织中，从吸收衬度图像很难识别木材和石墨的异物(图5E)。在相位衬度图像的结果中，我们能够从肌肉中识别所有六种类型的异物，其中0.2mm的塑料片是最难识别的(图5F)。在肌肉层，暗场衬度图像与脂肪层相似，两者都很难找到塑料薄膜(图5G)。
图6吸收成像的信号值除以暗场的信号值得到的R值。A、不同组织和各种异物的R值直方图；B、肌肉组织中异物的吸收衬度图像；C、肌肉组织中异物的暗场衬度图像。D、将肌肉内异物的暗场值除以衰减值，重建成图像。红色箭头指向塑料板的位置。
通过将吸收成像的信号值中每个像素的信号值除以暗场的信号值中同一像素的信号值，在一定程度上消除了厚度的影响，得到的值(R值)可以代表物质的本质属性。从不同组织和各种异物的R值柱状图(图6A)可以看出，脂肪和肌肉的R值无显著性差异，但异物和软组织的R值有显著性差异(P＜0.001)。其中，木材、塑料、石墨和陶瓷的R值甚至比软组织的R值高出数倍。在目前的硬件条件下，通过检测物体的R值很难直接判断异物的类型。
将肌肉内异物的吸收成像的信号值除以暗场的信号值，再重建为图像，肌肉内塑性异物的成像得到明显改善(图6D)。值得注意的是，塑料和塑料轮廓的R值是肌肉R值的10倍。因此，难以通过吸收衬度图像和暗场衬度图像(图6B和6C)清楚地判断的塑料片可以通过R值重建的图像容易地从肌肉中识别。
基于光栅的多模态X射线成像与传统的吸收X射线成像相比，可以获得病变的三个方面的信息：样品的吸收能力、电子密度分布和小角度散射能力，从而使常见异物的检测更加全面。在不同的皮下层次，不同的对比结果对不同类型异物的鉴别有不同的优势。
基于光栅的多模态X射线成像在组织识别上具有很高的灵敏度。
Claims (8)
1.一种用于检测软组织异物的多模态X射线成像设备的成像检测方法，其特征在于，基于光栅的多模态X射线成像检测设备包括：
产生X射线的光源；
采用X射线吸收光栅从而把从所述光源入射的X射线分离成多个线光源的源光栅；
能使X射线波前产生相位分布从而在光栅下游产生规则的竖状条纹的相位光栅；
采用X射线吸收光栅，配备步进位移控制器，将条纹图样的衰减、扭曲、或模糊信息转化为下游X射线强度信息的分析光栅；
接收X射线并转化为数字信号的X射线探测器；
其光路依次是X射线光源、源光栅、相位光栅、分析光栅和X射线探测器，X射线由X射线光源发出，依次经过光路上的各个组分；样品通常放置在相位光栅前方；
成像检测方法包括以下步骤：
工作状态调节步骤，开启X射线光源和X射线探测器，调节各部件至工作状态，并采集参考序列图像；即，对图像中的每一个像素都获得了无样品情况下的位置-光强曲线，即背景的步进曲线IB；
采集样本步骤，将样品放置于样品台上，校准后，移动样品或者对系统光栅部件进行步进调节，X射线探测器采集到样品序列图像；即，对图像中的每一个像素都获得了有样品情况下的位置-光强曲线，即样品的步进曲线IC；
信号分析步骤，经过校准和预处理之后，对参考序列图像和样品序列图像分别进行正弦拟合分析，然后进行比较和运算，进而获得关于样品的多模态图像；所述多模态图像包括吸收衬度图像imageA、相位衬度图像imagePS、暗场衬度图像imageS；
图像增强步骤，根据多模态图像的局域信息，可以对样品内部的感兴趣区域进行筛选区分，对每个像素的背景、样品进行处理，进而获得样品的R值图像；
数据处理方法包括以下步骤：
a、分别对每一个像素的背景的步进曲线IB、样品的步进曲线IC进行正弦拟合，获得步进曲线的直流分量kB和kC、幅值AB和AC、相角φB和φC；
b、计算-log(kC/kB)，获得吸收衬度图像imageA；
c、计算φC-φB，获得相位衬度图像imagePS；
d、计算-log((ACkB)/(ABkC))，获得暗场衬度图像imageS
e、计算imageA/imageS获得样品的R值图像。
2.根据权利要求1所述用于检测软组织异物的多模态X射线成像设备的成像检测方法，其特征在于，工作状态调节步骤和采集样本步骤中，在吸收衬度图像和暗场衬度图像中，选择目标像素线段，从源文件中获取该线段每个像素的信号强度，并根据位置序列绘制位置-光强曲线。
3.根据权利要求1所述用于检测软组织异物的多模态X射线成像设备的成像检测方法，其特征在于，步骤d中将吸收衬度图像imageA、相位衬度图像imagePS、暗场衬度图像imageS分别转换成256级的红、蓝、绿三原色，然后进行图像合并，得到伪彩色图像。
4.一种用于权利要求1～3任意一项所述用于检测软组织异物的多模态X射线成像设备，其特征在于包括：工作状态调节系统、数据采集系统和计算机数据处理系统；
所述数据采集系统包括：产生X射线的光源、源光栅、相位光栅、分析光栅、X射线探测器，以及精密调节模块。
5.根据权利要求4所述用于检测软组织异物的多模态X射线成像设备，其特征在于，所述工作状态调节系统包含针对各光学硬件的三维角度和位移调节器件，以及样品的宏观调节器件；通过对调试结果或已有数据的回溯，以及对当前X射线探测器采集到的图像结果的分析和判断，调节各光学硬件的位置、角度和样品的相对位置，使得系统满足成像条件。
6.根据权利要求4所述用于检测软组织异物的多模态X射线成像设备，其特征在于，所述计算机数据处理系统包括调节控制模块、数据采集模块、数据分析模块、图像增强模块；
1)调节控制模块负责工作状态调节过程中关键数据的记录、调取和反馈，为工作状态调节系统提供依据；
2)数据采集模块负责存储参考序列图像和样品序列图像，并对其进行预处理；
3)数据分析模块负责获取样品的多模态图像，即吸收衬度图像imageA、相位衬度图像imagePS、暗场衬度图像imageS；
4)图像增强模块负责对感兴趣区域进行筛选区分，进而获得样品的R值图像。
7.根据权利要求4所述用于检测软组织异物的多模态X射线成像设备，其特征在于，X射线的光源为焦斑尺寸为小于1mm的临床钨靶X射线管，工作电压为小于200kVp。
8.根据权利要求4所述用于检测软组织异物的多模态X射线成像设备，其特征在于，源光栅与相位光栅之间的距离为2156mm；
相位光栅与分析仪光栅之间的距离为308mm。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110063422 | 2021-01-18 | ||
| CN2021100634223 | 2021-01-18 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113171116A true CN113171116A (zh) | 2021-07-27 |
Family
ID=76924885
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110391180.0A Pending CN113171116A (zh) | 2021-01-18 | 2021-04-12 | 一种用于检测软组织异物的多模态x射线成像设备及其方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113171116A (zh) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101943668A (zh) * | 2009-07-07 | 2011-01-12 | 清华大学 | X射线暗场成像系统和方法 |
| CN102221565A (zh) * | 2010-04-19 | 2011-10-19 | 清华大学 | X射线源光栅步进成像系统与成像方法 |
| CN102971620A (zh) * | 2010-06-28 | 2013-03-13 | 保罗·谢勒学院 | 使用平面几何结构的光栅设备进行x射线相衬成像和暗场成像的方法 |
| CN104137149A (zh) * | 2012-02-24 | 2014-11-05 | 保罗·谢勒学院 | 用于人体组织中的不同类型微钙化的非侵入性分类的系统 |
| CN104970815A (zh) * | 2014-04-04 | 2015-10-14 | 曹红光 | 基于光栅相位衬度和光子计数的x射线成像系统及方法 |
| US20160113611A1 (en) * | 2013-05-10 | 2016-04-28 | Paul Scherrer Institut | Quantitative x-ray radiology using the absorption and scattering information |
| CN110068585A (zh) * | 2019-03-15 | 2019-07-30 | 中国工程物理研究院流体物理研究所 | 双能x射线光栅干涉成像系统及方法 |
-
2021
- 2021-04-12 CN CN202110391180.0A patent/CN113171116A/zh active Pending
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101943668A (zh) * | 2009-07-07 | 2011-01-12 | 清华大学 | X射线暗场成像系统和方法 |
| CN102221565A (zh) * | 2010-04-19 | 2011-10-19 | 清华大学 | X射线源光栅步进成像系统与成像方法 |
| CN102971620A (zh) * | 2010-06-28 | 2013-03-13 | 保罗·谢勒学院 | 使用平面几何结构的光栅设备进行x射线相衬成像和暗场成像的方法 |
| CN104137149A (zh) * | 2012-02-24 | 2014-11-05 | 保罗·谢勒学院 | 用于人体组织中的不同类型微钙化的非侵入性分类的系统 |
| US20160113611A1 (en) * | 2013-05-10 | 2016-04-28 | Paul Scherrer Institut | Quantitative x-ray radiology using the absorption and scattering information |
| CN104970815A (zh) * | 2014-04-04 | 2015-10-14 | 曹红光 | 基于光栅相位衬度和光子计数的x射线成像系统及方法 |
| CN110068585A (zh) * | 2019-03-15 | 2019-07-30 | 中国工程物理研究院流体物理研究所 | 双能x射线光栅干涉成像系统及方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9420983B2 (en) | System for quantification of neovasculature in CT volumes | |
| US7574249B2 (en) | Device-less gating of physiological movement for improved image detection | |
| US7778380B2 (en) | Data handling and analysis in computed tomography with multiple energy windows | |
| US7292721B2 (en) | Calibration of tissue densities in computerized tomography | |
| EP2782504B1 (en) | Method and device for imaging soft body tissue using x-ray projection and optical tomography | |
| JP4402436B2 (ja) | 気道測定のための方法及びシステム | |
| US8186880B1 (en) | Extended and fixed INTable simultaneously imaged calibration and correction methods and references for 3-D imaging devices | |
| CN1762304A (zh) | 在三维成像尤其是乳腺x射线照片中给出检查结果的方法 | |
| US20140050378A1 (en) | System and method for detecting materials or disease states using multi-energy computed tomography | |
| JP2007152104A (ja) | 悪性腫瘍を自動的に特徴評価する方法及び装置 | |
| CN104408758A (zh) | 一种低剂量能谱ct图像处理方法 | |
| HUE035719T2 (en) | New imaging mode with penetrating radiation | |
| CN111312373B (zh) | 一种pet/ct图像融合自动标注方法 | |
| CN104873213B (zh) | 基于x‑射线数字影像的骨骼病变评估方法和装置 | |
| CN110613471A (zh) | 一种ct系统及ct图像重建方法 | |
| De Ruiter et al. | MARS imaging and reconstruction challenges | |
| CN113171116A (zh) | 一种用于检测软组织异物的多模态x射线成像设备及其方法 | |
| CN109709127A (zh) | 低散射x射线荧光ct成像方法 | |
| Mahnken et al. | Spectral ρZ-Projection Method for Characterization of Body Fluids in Computed Tomography: Ex Vivo Experiments1 | |
| Chacko et al. | Three-dimensional reconstruction of transillumination tomographic images of human breast phantoms by red and infrared lasers | |
| Rajeswari Amma | Study of bone-metal interface in orthopaedic application using spectral CT | |
| Dahal | Characterization of atherosclerotic plaque with MARS Imaging | |
| Uddin | Spectral image quality in MARS scanners. | |
| da Silva Carvalho et al. | Characterization of hydroxyapatite tablet by X-ray computed microtomography aiming at the construction of a low-cost bone tissue phantom | |
| Pacureanu et al. | 3D non-linear enhancement of tubular microscopic bone porosities |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210727 |
|
| RJ01 | Rejection of invention patent application after publication |