一种基于多波长光学图像特征融合的头部异质体成像方法及系统技术方案

本发明专利技术涉及一种一种基于多波长光学图像特征融合的头部异质体成像方法及系统，该方法包括：利用历史CT图像建立带有异质体的头部分层模型，根据所述头部分层模型异质体位置确定光源探测器放置位置；在所述光源探测器放置位置轮换放置光源与探测器进行光学仿真，获取不同大小异质体在不同光源探测器放置位置下的仿真吸光度数据与头部模型吸收系数数据；改变光源波长，获取不同光源波长下的单波长仿真吸光度数据；根据不同光源波长下的单波长仿真吸光度数据，采用多波长光学图像特征融合方法获取多波长仿真光学特征数据并与所述头部模型吸收系数数据一起建立数据集；将所述数据集划分为训练集与预测集，训练集用于对深度学习模型进行训练来获取训练好的深度学习网络模型；根据所述预测集，采用所述深度学习网络模型确定头部模型的吸收系数分布，根据吸收系数的分布重建头部异质体图像，本发明专利技术有效缓解DOT方法成像过程中存在病态逆问题，提高头部异质体成像的效果，以便精确监测头部异质体。以便精确监测头部异质体。以便精确监测头部异质体。

【技术实现步骤摘要】
一种基于多波长光学图像特征融合的头部异质体成像方法及系统


[0001]本专利技术涉及生物检测
，特别是涉及一种基于多波长光学图像特征融合的头部异质体成像方法及系统。
技术背景
[0002]颅内出血形成异质体是颅脑损伤和心血管疾病中常见且严重的继发性病变，约占闭合性颅脑损伤的10％和重型颅脑损伤的40％～50％，对于此类病人一般需要进行开颅手术，但术后常伴随再次出现颅内出血的情况，这时医生需要根据异质体大小判断是否需要再次通过手术清除异质体或通过药物治疗使异质体自行吸收，因此需要术后实时监测异质体变化情况。
[0003]目前，头部异质体检测的主要方法仍是如CT(电子计算机断层扫描)、MRI(磁共振成像)这类影像学检测方法，由于设备的大小和使用方式的限制，此类方法均无法对头部异质体变化实现连续监控，现有的头部异质体连续监测算法大多采用的是扩散相干层析成像的方法(diffuse optical tomography，DOT)，但传统的DOT成像过程中存在病态逆问题，这导致其由于输入信息不足，成像效果不佳；因此，亟需一种方法或系统可以有效缓解DOT方法成像过程中存在病态逆问题，提高头部异质体成像的效果，以便精确监测头部异质体。

技术实现思路

[0004]本专利技术的目的是提供一种头部异质体成像方法及系统，有效提高异质体成像的效果。
[0005]为了实现上述目的，本专利技术提供了如下方案：
[0006]一种基于多波长光学图像特征融合的头部异质体成像方法，包括：
[0007]S100，利用历史CT图像建立带有异质体的头部分层模型，根据所述头部分层模型异质体位置确定光源探测器放置位置；
[0008]S200，在所述光源探测器放置位置轮换放置光源与探测器进行光学仿真，获取不同大小异质体在不同光源探测器放置位置下的仿真吸光度数据与头部模型吸收系数数据；
[0009]S300，改变光源波长，获取不同光源波长下的单波长仿真吸光度数据；
[0010]S400，根据不同光源波长下的单波长仿真吸光度数据，采用多波长光学图像特征融合方法获取多波长仿真光学特征数据并与所述头部模型吸收系数数据一起建立数据集；
[0011]S500，将所述数据集划分为训练集与预测集，训练集用于对深度学习模型进行训练来获取训练好的深度学习网络模型；
[0012]S600，根据所述预测集，采用所述深度学习网络模型确定头部模型的吸收系数分布，根据吸收系数的分布重建头部异质体图像。
[0013]在S100中，利用历史CT图像建立带有异质体的头部分层模型，根据所述头部分层模型异质体位置确定光源探测器放置位置，具体包括：
[0014]S101，根据所述历史CT图像确定异质体在所述头部模型的位置，将头部模型存在异质体的横断面区域作为重点区域；
[0015]S102，将所述重点区域分为头皮头骨层、脑脊液层、灰质白质层和异质体区域作为头部分层模型，用有限元分割方法将所述头部分层模型分割为有限元网格；
[0016]可选的，S103，将所述异质体区域最为临近的头皮头骨层部分的上方作为光源探测器放置区域，在所述的光源探测器放置区域等间距确定n个点的位置作为光源探测器放置位置。
[0017]在S200中，在所述光源探测器放置位置轮换放置光源与探测器进行光学仿真，获取不同大小异质体在不同光源探测器放置位置下的仿真吸光度数据与头部模型吸收系数数据，具体包括：
[0018]S201，根据所述的光源探测器放置位置将其n个点标记为D1
‑
Dn；
[0019]S202，通过边界提取及放缩方法放缩异质体得到不同大小的仿真异质体组成仿真异质体样本数据集；
[0020]S203，依次在D1
‑
Dn位置放置单独光源，设置光源波长为λ1，在其余n
‑
1个位置放置探测器，在所述不同大小的仿真异质体组成仿真异质体样本数据集上进行光学仿真，得到光源波长为λ1下的不同大小异质体在不同光源探测器放置位置下的仿真吸光度数据与头部模型吸收系数数据。
[0021]在S300中，改变光源波长，获取不同光源波长下的单波长仿真吸光度数据，具体包括：
[0022]S301，将光源波长由λ1变为λ2，根据所述S200的步骤，获取光源波长为λ2下的单波长仿真吸光度数据与头部模型吸收系数数据；
[0023]S302，将光源波长由λ2变为λ3，根据所述S200的步骤，获取光源波长为λ3下的单波长仿真吸光度数据与头部模型吸收系数数据。
[0024]在S400中，根据不同光源波长下的单波长仿真吸光度数据，采用多波长光学图像特征融合方法获取多波长仿真光学特征数据并与所述头部模型吸收系数数据一起建立数据集，具体包括：
[0025]S401，分别对所述不同光源波长下的单波长仿真吸光度数据采用光学图像映射特征提取方法获取不同光源波长下的单波长仿真光学特征图像；
[0026]其中，在S401中，光学图像映射特征提取方法为：
[0027]S4011，将所述每个大小的异质体对应的n组1*n的仿真吸光度数据整合为n*n的光学灰度图像，对n*n的光学灰度图像依次进行卷积、归一化、映射、卷积、归一化、映射、池化、插值的操作使其变为m*m的光学初步特征图像；
[0028]S4012，将所述光学初步特征图像滤波处理后进行二维离散小波变换获取所述单波长仿真光学特征图像；
[0029]S402，对所述不同光源波长下的单波长仿真光学特征图像采用多波长图像融合方法获取多波长仿真光学特征图像；
[0030]其中，在S402中，多波长图像融合方法为：
[0031]S4021，在所述不同光源波长下的单波长仿真光学特征图像中选取2个波长的单波长仿真光学特征图像进行双波长融合得到初步融合仿真光学特征图像；
[0032]其中，在S4021中，双波长图像融合方法为：
[0033]S40211，计算2个波长的单波长仿真光学特征图像的相似性矩阵S；
[0034]S＝[2k*(P
λ1
×
P
λ2
)]÷
[k*P
λ12
+k*P
λ22
][0035]其中P
λ1
、P
λ2
为2个波长的单波长仿真光学特征图像，“*”表示卷积操作，
“×”
表示矩阵对应元素相乘，
“÷”
表示矩阵对应元素相除，K为3*3的卷积核；
[0036]S40212，根据所述相似性矩阵S计算初步融合仿真光学特征图像；
[0037]P
out
＝S
×
P
λ1
+(I
‑
S)
×
P
λ2
[0038]其中I是所有元素都为1的矩阵，Pout是初步融合仿真光学特征图像；
[0039]S4本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于多波长光学图像特征融合的头部异质体成像方法，其特征在于，包括：S100，利用历史CT图像建立带有异质体的头部分层模型，根据所述头部分层模型异质体位置确定光源探测器放置位置；S200，在所述光源探测器放置位置轮换放置光源与探测器进行光学仿真，获取不同大小异质体在不同光源探测器放置位置下的仿真吸光度数据与头部模型吸收系数数据；S300，改变光源波长，获取不同光源波长下的单波长仿真吸光度数据；S400，根据不同光源波长下的单波长仿真吸光度数据，采用多波长光学图像特征融合方法获取多波长仿真光学特征数据并与所述头部模型吸收系数数据一起建立数据集；S500，将所述数据集划分为训练集与预测集，训练集用于对深度学习模型进行训练来获取训练好的深度学习网络模型；S600，根据所述预测集，采用所述深度学习网络模型确定头部模型的吸收系数分布，根据吸收系数的分布重建头部异质体图像。2.根据权力要求1所述的一种基于多波长光学图像特征融合的异质体成像方法，其特征在于，在S100中，利用历史CT图像建立带有异质体的头部分层模型，根据所述头部分层模型异质体位置确定光源探测器放置位置，具体包括：S101，根据所述历史CT图像确定异质体在所述头部模型的位置，将头部模型存在异质体的横断面区域作为重点区域；S102，将所述重点区域分为头皮头骨层、脑脊液层、灰质白质层和异质体区域作为头部分层模型，用有限元分割方法将所述头部分层模型分割为有限元网格；S103，将所述异质体区域最为临近的头皮头骨层部分的上方作为光源探测器放置区域，在所述的光源探测器放置区域等间距确定n个点的位置作为光源探测器放置位置。3.根据权力要求1所述的一种基于多波长光学图像特征融合的异质体成像方法，其特征在于，在S200中，在所述光源探测器放置位置轮换放置光源与探测器进行光学仿真，获取不同大小异质体在不同光源探测器放置位置下的仿真吸光度数据与头部模型吸收系数数据，具体包括：S201，根据所述的光源探测器放置位置将其n个点标记为D1
‑
Dn；S202，通过边界提取及放缩方法放缩异质体得到不同大小的仿真异质体组成仿真异质体样本数据集；S203，依次在D1
‑
Dn位置放置单独光源，设置光源波长为λ1，在其余n
‑
1个位置放置探测器，在所述不同大小的仿真异质体组成仿真异质体样本数据集上进行光学仿真，得到光源波长为λ1下的不同大小异质体在不同光源探测器放置位置下的仿真吸光度数据与头部模型吸收系数数据。4.根据权力要求1所述的一种基于多波长光学图像特征融合的异质体成像方法，其特征在于，在S300中，改变光源波长，获取不同光源波长下的单波长仿真吸光度数据，具体包括：S301，将光源波长由λ1变为λ2，根据所述S200的步骤，获取光源波长为λ2下的单波长仿真吸光度数据与头部模型吸收系数数据；S302，将光源波长由λ2变为λ3，根据所述S200的步骤，获取光源波长为λ3下的单波长仿真吸光度数据与头部模型吸收系数数据。
5.根据权力要求1所述的一种基于多波长光学图像特征融合的异质体成像方法，其特征在于，在S400中，根据不同光源波长下的单波长仿真吸光度数据，采用多波长光学图像特征融合方法获取多波长仿真光学特征数据并与所述头部模型吸收系数数据一起建立数据集，具体包括：S401，分别对所述不同光源波长下的单波长仿真吸光度数据采用光学图像映射特征提取方法获取不同光源波长下的单波长仿真光学特征图像；其中，在S401中，光学图像映射特征提取方法为：S4011，将所述每个大小的异质体对应的n组1*n的仿真吸光度数据整合为n*n的光学灰度图像，对n*n的光学灰度图像依次进行卷积、归一化、映射、卷积、归一化、映射、池化、插值的操作使其变为m*m的光学初步特征图像；S4012，将所述光学初步特征图像滤波处理后进行二维离散小波变换获取所述单波长仿真光学特征图像；S402，对所述不同光源波长下的单波长仿真光学特征图像采用多波长图像融合方法获取多波长仿真光学特征图像；其中，在S402中，多波长图像融合方法为：S4021，在所述不同光源波长下的单波长仿真光学特征图像中选取2个波长的单波长仿真光学特征图像进行双波长融合得到初步融合仿真光学特征图像；其中，在S4021中，双波长图像融合方法为：S40211，计算2个波长的单波长仿真光学特征图像的相似性矩阵S；S＝[2k*(P
λ1
×
P
λ2
)]
÷
[k*P
λ12
+k*P
λ22
]其中P
λ1
、P
λ2
为2个波长的单波长仿真光学特征图像，“*”表示卷积操作，
...

【专利技术属性】
技术研发人员：韩广，潘中华，王慧泉，
申请(专利权)人：天津工业大学，
类型：发明
国别省市：