1.本发明涉及医疗成像技术领域，具体为一种基于深度学习算法的微波人体实时穿透式成像方法。
背景技术：2.现阶段较为成熟的医疗成像仪器主要有核磁共振成像(magnetic resonance imaging,mri)、断层扫描(computed tomography,ct)、x射线成像、正电子发射断层成像(positron emission tomography,pet)、超声成像(ultrasound,us)，各个方法的成像原理不同，得到的成像信息不同，因此拥有各自的优势和劣势。
3.其中，mri分辨率较高，但价格昂贵，永磁体体积巨大，并且受测患者体内不可以装有金属医疗支架、钉等；ct拥有很高的空间分辨率，可实现组织的灰度成像，是目前医疗成像的主流设备，但基于x射线的成像存在电离辐射，容易引起额外的组织病变，对人体始终存在不可忽视的健康隐患，无法进行实时监测，此外x光成像对比度很低，误诊率较高；mri和ct检测系统设备规模较大，不可携带，均不适合医者现场使用，偏远地区相关医疗建设情况也较少。pet利用放射性元素衰变标记进行肿瘤和癌细胞扩散的追踪诊断，医学原理可靠，诊断精度很高，但对检查医生经验要求较高，且费用十分昂贵，因此不易推广；超声成像对软组织结构分辨率较高，能够实现实时成像和动态连续检测，且无痛苦与电离辐射伤害，但该方法对气体和骨骼传导深度有限，因此对肺、骨骼、肥胖者成像效果欠佳。
4.微波具有穿透性，通过反演组织介电常数判断病灶种类，可对人体进行定量成像。微波器件较为成熟，成本低，且可以集成为便携式设备。微波对人体没有电离辐射，可以对人体进行实时监控。但因传统逆散射成像算法中涉及到正问题迭代，导致对算力要求高，计算速度慢，因此不具备实时性。
5.目前市面上的成像系统，在定性成像判断上通过有经验的医生进行病灶划分，生成训练集，训练集获取困难大，数据获取不易；而基于微波的穿透式实时定量成像方案，可以直接通过计算获得病灶参数，摆脱人为判断，且训练集可通过模拟获得，保护患者数据。
6.为此，我们提出了一种基于深度学习算法的微波人体实时穿透式成像方法。
技术实现要素：7.针对现有技术的不足，本发明提供了一种基于深度学习算法的微波人体实时穿透式成像方法，以解决上述背景技术中提出的问题。
8.为实现上述目的，本发明提供如下技术方案：一种基于深度学习算法的微波人体实时穿透式成像方法，包括以下步骤：
9.s1、首先确定成像系统参数并接收原始测试数据；
10.s2、通过线性算法获得只有线性信息的介电常数粗糙图像；
11.s3、基于注意力机制的机器学习网络将介电常数粗糙图像作为神经网络的输入；
12.s4、加入高次散射信息，得到真实介电常数的成像结果。
13.进一步优化本技术方案，所述s1中，确定的成像系统参数为散射参数，所述散射参数用于描述频域通道的频域特性，原始测试数据为散射场数据，所述散射场数据基于多发多收的阵列天线进行采集，所述散射参数基于成像系统中的矢量网络分析仪进行采集。
14.进一步优化本技术方案，所述阵列天线在目标成像区域四周等间距均匀环绕分布，阵列天线的个数根据成像分辨率需求选择为n，其中n的理论最小值为2，阵列天线中的天线为收发一体天线，当选取1号天线为发射天线，采集1号天线自身的输入反射系数，采集2～n号天线的传输效率参数。
15.进一步优化本技术方案，在电脑端读取成像系统中矢量网络分析仪采集到的数据，根据接收天线的情况进行存储；切换发射天线，采集剩余所有天线测量参数，形成n*n的数据矩阵，如下式(1.1)所示：
[0016][0017]
根据成像系统中的矢量网络分析仪采集到的散射参数，对应转换到系统的电场分布，在天线分布确定的情况下，通过理论计算得到收发天线处的入射场根据矢量网络分析仪测试得到收发天线处的总电场由上述二者得到一个n*n的散射场数据作为系统输入，其中n为阵列天线的个数。
[0018]
进一步优化本技术方案，所述s1中，从输入原始测试数据到得到介电常数分布的过程中，通过理论计算，得到对应位置的格林函数，用于在传统的散射场直接到成像结果的映射中加入物理意义，并将待测目标区域划分为m*m的网格。
[0019]
进一步优化本技术方案，所述s2中，所述线性算法采用反向传播bps成像算法，基于电磁波引起的感应电流是线性的，并忽略高阶散射效应的考虑，线性算法的推导如下所示：
[0020][0021]
其中，算符h表示矩阵的共轭转置，χ是未确定的尺度系数，为了使尽可能准确，将价值函数定义为：
[0022][0023]
为了取到价值函数的最小值，f
obj
(χ)对χ的一阶导数应为0，因此可以得到：
[0024][0025]
于是可以得到而在待测目标区域m*m的每个像素点上，总电场等于入射场和计算所得散射场的和，表示为：
[0026][0027]
对于每一个成像像素点，对比度和感应电流之间的关系可以表示为：
[0028][0029]
对于所有的点，这是一个最小二乘的问题，最终得到介电常数的解：
[0030][0031]
进一步优化本技术方案，所述s3中，机器学习网络需要以计算或者仿真数据作为训练集训练网络，并改变损失函数对成像进行优化。
[0032]
进一步优化本技术方案，所述训练集的获取方式包括以下三种：
[0033]
1)使用自编写矩量法mom程序对模型进行理论计算，建立数据集；
[0034]
2)利用电磁仿真软件，包括但不限于hfss，cst，feko对系统进行仿真，将仿真电场结果作为训练集；
[0035]
3)根据实际模型测试得到的电场值，作为训练集。
[0036]
进一步优化本技术方案，所述s3中，机器学习网络采用基于语义分割的vision transformer模型，即vit模型，该模型分为编码器与解码器两部分，首先将输入图像分割成块，将其映射为一个线性嵌入序列，加入位置编码后用编码器进行编码，再由mask transformer将编码器和类嵌入的输出进行解码，上采样后应用argmax给每个像素一一分好类，输出最终的像素分割图。
[0037]
进一步优化本技术方案，所述机器学习网络基于vit模型可将低分辨介电常数图像的高分辨率信息补全到高分辨率，生成超分辨率图像，机器学习网络一旦完成训练，即可根据输入瞬间得到输出图像，从而基于transformer的超分辨率方法实现实时成像。
[0038]
与现有技术相比，本发明提供了一种基于深度学习算法的微波人体实时穿透式成像方法，具备以下有益效果：
[0039]
1、该基于深度学习算法的微波人体实时穿透式成像方法，通过计算，仿真，实测完成训练集数据收集，利用线性算法得到粗糙的初始粗糙成像结果，基于深度学习，能够实现高精度的实时反演成像，为搭建微波实时定性成像系统提供了有效的技术支持。
[0040]
2、该基于深度学习算法的微波人体实时穿透式成像方法，利用微波成像原理实现了辐射功率小、安全无创、非侵入式检测的目的；系统设计实现了小型化，便携式的目的；使用深度学习算法，提高了成像分辨率和实现了实时成像，可应用于人脑、胳膊、膝盖等部位的医疗成像检测。
附图说明
[0041]
图1为本发明提出的一种基于深度学习算法的微波人体实时穿透式成像方法的流程示意图；
[0042]
图2为本发明提出的一种基于深度学习算法的微波人体实时穿透式成像方法的步骤s1的数据采集示意图；
[0043]
图3为本发明提出的一种基于深度学习算法的微波人体实时穿透式成像方法的机器学习网络的示意图；
[0044]
图4为本发明提出的一种基于深度学习算法的微波人体实时穿透式成像方法的结果示意图。
具体实施方式
[0045]
下面将结合本发明的实施例，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0046]
实施例一：
[0047]
请参阅图1，一种基于深度学习算法的微波人体实时穿透式成像方法，包括以下步骤：
[0048]
s1、首先确定成像系统参数并接收原始测试数据。
[0049]
其中，确定的成像系统参数为散射参数，所述散射参数用于描述频域通道的频域特性，原始测试数据为散射场数据，所述散射场数据基于多发多收的阵列天线进行采集，所述散射参数基于成像系统中的矢量网络分析仪进行采集。如图2所示，所述阵列天线在目标成像区域四周等间距均匀环绕分布，阵列天线的个数根据成像分辨率需求选择为n(n≥2，理论最小值为2)，阵列天线中的天线为收发一体天线，当选取1号天线为发射天线，采集1号天线自身的s11参数(输入反射系数)，采集2～n号天线的s21参数(即天线的传输效率参数)。在电脑端读取成像系统中矢量网络分析仪采集到的数据，根据接收天线的情况进行存储；切换发射天线，采集剩余所有天线测量参数，形成n*n的数据矩阵，如下式(1.1)所示：
[0050][0051]
根据成像系统中的矢量网络分析仪采集到的散射参数，对应转换到系统的电场分布，在天线分布确定的情况下，通过理论计算得到收发天线处的入射场根据矢量网络分析仪测试得到收发天线处的总电场由上述二者得到一个n*n的散射场数据(n为收发天线个数)作为系统输入。从输入原始测试数据到得到介电常数分布的过程中，通过理论计算，得到对应位置的格林函数，用于在传统的散射场直接到成像结果的映射中加入物理意义，并将待测目标区域划分为m*m的网格。
[0052]
s2、通过线性算法获得只有线性信息的介电常数粗糙图像。
[0053]
其中，所述线性算法采用反向传播bps(back-propagation scheme)成像算法，基于电磁波引起的感应电流是线性的，并忽略高阶散射效应的考虑，线性算法的推导如下所示：
[0054][0055]
其中，算符h表示矩阵的共轭转置，χ是未确定的尺度系数，为了使尽可能准确，将价值函数定义为：
[0056][0057]
为了取到价值函数的最小值，f
obj
(χ)对χ的一阶导数应为0，因此可以得到：
[0058][0059]
于是可以得到而在待测目标区域m*m的每个像素点上，总电场等于入射场和计算所得散射场的和，表示为：
[0060][0061]
对于每一个成像像素点，对比度和感应电流之间的关系可以表示为：
[0062][0063]
对于所有的点，这是一个最小二乘的问题，最终得到介电常数的解：
[0064][0065]
s3、基于注意力机制的机器学习网络将介电常数粗糙图像作为神经网络的输入。
[0066]
其中，机器学习网络需要以计算或者仿真数据作为训练集训练网络，并改变损失函数对成像进行优化。
[0067]
更进一步的，所述训练集的获取方式包括以下三种：
[0068]
1)使用自编写矩量法mom程序对模型进行理论计算，建立数据集；
[0069]
2)利用电磁仿真软件，包括但不限于hfss，cst，feko对系统进行仿真，将仿真电场结果作为训练集；
[0070]
3)根据实际模型测试得到的电场值，作为训练集。
[0071]
机器学习网络采用一种基于语义分割的vision transformer(vit)模型。该模型分为编码器(decoder)与解码器(encoder)两部分。首先将输入图像分割成块(patches)，将其映射为一个线性嵌入序列，加入位置编码(position encoding)后用编码器进行编码，再由mask transformer将编码器和类嵌入的输出进行解码，上采样后应用argmax给每个像素一一分好类，输出最终的像素分割图。基于语义分割的vit模型可将低分辨介电常数图像的高分辨率信息补全到高分辨率，生成超分辨率图像。神经网络一旦完成训练，即可根据输入瞬间得到输出图像，因此基于transformer的超分辨率方法可以实现实时成像。
[0072]
s4、网络训练完成后，输入采集到的数据，加入高次散射信息，得到真实介电常数的成像结果，实现高精度实时成像。
[0073]
下表中展示利用此公开中的方式进行成像，较传统s-dbim算法，在准确度和时间上，都有大幅提高。
[0074][0075]
因此，该基于深度学习算法的微波人体实时穿透式成像方法，通过计算，仿真，实测完成训练集数据收集，利用线性算法得到粗糙的初始粗糙成像结果，基于深度学习，能够实现高精度的实时反演成像，为搭建微波实时定性成像系统提供了有效的技术支持。利用微波成像原理实现了辐射功率小、安全无创、非侵入式检测的目的；系统设计实现了小型化，便携式的目的；使用深度学习算法，提高了成像分辨率和实现了实时成像，可应用于人脑、胳膊、膝盖等部位的医疗成像检测。
[0076]
实施例二：
[0077]
线性算法在实施例一中以bps(back-propagation scheme,反向传播方法)算法为例，具体实施中对所出现的线性算法，在本实施例中包括但不局限于ba(born approximation，伯恩近似法)，ra(rytov approximation，rytov近似法),mcc(major current coefficients，主电流系数法),mba(modified born approximation，改进型伯恩近似法),mcc+mba(主电流系数和改进型伯恩近似法)等算法。
[0078]
结合机器学习部分，实施例一中采用了基于transformer的网络架构，在具体实施过程中，在本实施例中包括但不局限于使用基于gan(生成对抗网络)架构的sr-gan(超分辨率图像复原的生成对抗网络)网络，u-net cnn(卷积神经网络)等。
[0079]
本发明的有益效果是：
[0080]
1、该基于深度学习算法的微波人体实时穿透式成像方法，通过计算，仿真，实测完成训练集数据收集，利用线性算法得到粗糙的初始粗糙成像结果，基于深度学习，能够实现高精度的实时反演成像，为搭建微波实时定性成像系统提供了有效的技术支持。
[0081]
2、该基于深度学习算法的微波人体实时穿透式成像方法，利用微波成像原理实现了辐射功率小、安全无创、非侵入式检测的目的；系统设计实现了小型化，便携式的目的；使用深度学习算法，提高了成像分辨率和实现了实时成像，可应用于人脑、胳膊、膝盖等部位的医疗成像检测。
[0082]
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
[0083]
尽管已经示出和描述了本发明的实施例，对于本领域的普通技术人员而言，可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换
和变型，本发明的范围由所附权利要求及其等同物限定。
技术特征：1.一种基于深度学习算法的微波人体实时穿透式成像方法，其特征在于，包括以下步骤：s1、首先确定成像系统参数并接收原始测试数据；s2、通过线性算法获得只有线性信息的介电常数粗糙图像；s3、基于注意力机制的机器学习网络将介电常数粗糙图像作为神经网络的输入；s4、加入高次散射信息，得到真实介电常数的成像结果。2.根据权利要求1所述的一种基于深度学习算法的微波人体实时穿透式成像方法，其特征在于，所述s1中，确定的成像系统参数为散射参数，所述散射参数用于描述频域通道的频域特性，原始测试数据为散射场数据，所述散射场数据基于多发多收的阵列天线进行采集，所述散射参数基于成像系统中的矢量网络分析仪进行采集。3.根据权利要求2所述的一种基于深度学习算法的微波人体实时穿透式成像方法，其特征在于，所述阵列天线在目标成像区域四周等间距均匀环绕分布，阵列天线的个数根据成像分辨率需求选择为n，其中n的理论最小值为2，阵列天线中的天线为收发一体天线，当选取1号天线为发射天线，采集1号天线自身的输入反射系数，采集2～n号天线的传输效率参数。4.根据权利要求3所述的一种基于深度学习算法的微波人体实时穿透式成像方法，其特征在于，在电脑端读取成像系统中矢量网络分析仪采集到的数据，根据接收天线的情况进行存储；切换发射天线，采集剩余所有天线测量参数，形成n*n的数据矩阵，如下式(1.1)所示：根据成像系统中的矢量网络分析仪采集到的散射参数，对应转换到系统的电场分布，在天线分布确定的情况下，通过理论计算得到收发天线处的入射场根据矢量网络分析仪测试得到收发天线处的总电场由上述二者得到一个n*n的散射场数据作为系统输入，其中n为阵列天线的个数。5.根据权利要求1所述的一种基于深度学习算法的微波人体实时穿透式成像方法，其特征在于，所述s1中，从输入原始测试数据到得到介电常数分布的过程中，通过理论计算，得到对应位置的格林函数，用于在传统的散射场直接到成像结果的映射中加入物理意义，并将待测目标区域划分为m*m的网格。6.根据权利要求1所述的一种基于深度学习算法的微波人体实时穿透式成像方法，其特征在于，所述s2中，所述线性算法采用反向传播bps成像算法，基于电磁波引起的感应电流是线性的，并忽略高阶散射效应的考虑，线性算法的推导如下所示：其中，算符h表示矩阵的共轭转置，χ是未确定的尺度系数，为了使尽可能准确，将价值函数定义为：
为了取到价值函数的最小值，f
obj
(χ)对χ的一阶导数应为0，因此可以得到：于是可以得到而在待测目标区域m*m的每个像素点上，总电场等于入射场和计算所得散射场的和，表示为：对于每一个成像像素点，对比度和感应电流之间的关系可以表示为：对于所有的点，这是一个最小二乘的问题，最终得到介电常数的解：7.根据权利要求1所述的一种基于深度学习算法的微波人体实时穿透式成像方法，其特征在于，所述s3中，机器学习网络需要以计算或者仿真数据作为训练集训练网络，并改变损失函数对成像进行优化。8.根据权利要求7所述的一种基于深度学习算法的微波人体实时穿透式成像方法，其特征在于，所述训练集的获取方式包括以下三种：1)使用自编写矩量法mom程序对模型进行理论计算，建立数据集；2)利用电磁仿真软件，包括但不限于hfss，cst，feko对系统进行仿真，将仿真电场结果作为训练集；3)根据实际模型测试得到的电场值，作为训练集。9.根据权利要求1所述的一种基于深度学习算法的微波人体实时穿透式成像方法，其特征在于，所述s3中，机器学习网络采用基于语义分割的vision transformer模型，即vit模型，该模型分为编码器与解码器两部分，首先将输入图像分割成块，将其映射为一个线性嵌入序列，加入位置编码后用编码器进行编码，再由mask transformer将编码器和类嵌入的输出进行解码，上采样后应用argmax给每个像素一一分好类，输出最终的像素分割图。10.根据权利要求9所述的一种基于深度学习算法的微波人体实时穿透式成像方法，其特征在于，所述机器学习网络基于vit模型可将低分辨介电常数图像的高分辨率信息补全到高分辨率，生成超分辨率图像，机器学习网络一旦完成训练，即可根据输入瞬间得到输出图像，从而基于transformer的超分辨率方法实现实时成像。
技术总结本发明公开了一种基于深度学习算法的微波人体实时穿透式成像方法，涉及医疗成像技术领域，包括以下步骤：S1、首先确定成像系统参数并接收原始测试数据；S2、通过线性算法获得只有线性信息的介电常数粗糙图像；S3、基于注意力机制的机器学习网络将介电常数粗糙图像作为神经网络的输入；S4、加入高次散射信息，得到真实介电常数的成像结果。该基于深度学习算法的微波人体实时穿透式成像方法，利用微波成像原理实现了辐射功率小、安全无创、非侵入式检测的目的；系统设计实现了小型化，便携式的目的；使用深度学习算法，提高了成像分辨率和实现了实时成像，可应用于人脑、胳膊、膝盖等部位的医疗成像检测。的医疗成像检测。的医疗成像检测。
技术研发人员：叶修竹 杜耐珂
受保护的技术使用者：北京理工大学
技术研发日：2022.04.18
技术公布日：2022/7/4
