CN109199381A - 一种全息微波弹性成像系统及其成像方法 - Google Patents
一种全息微波弹性成像系统及其成像方法 Download PDFInfo
- Publication number
- CN109199381A CN109199381A CN201811061290.5A CN201811061290A CN109199381A CN 109199381 A CN109199381 A CN 109199381A CN 201811061290 A CN201811061290 A CN 201811061290A CN 109199381 A CN109199381 A CN 109199381A
- Authority
- CN
- China
- Prior art keywords
- microwave
- signal
- ultrasonic
- electric field
- target object
- 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.)
- Granted
Links
- 238000003384 imaging method Methods 0.000 title claims description 48
- 238000001093 holography Methods 0.000 title 1
- 230000005684 electric field Effects 0.000 claims abstract description 74
- 238000012545 processing Methods 0.000 claims abstract description 64
- 238000000034 method Methods 0.000 claims abstract description 36
- 238000002604 ultrasonography Methods 0.000 claims abstract description 6
- 230000006870 function Effects 0.000 claims description 39
- 239000013598 vector Substances 0.000 claims description 24
- 238000006073 displacement reaction Methods 0.000 claims description 19
- 230000008569 process Effects 0.000 claims description 17
- 230000005669 field effect Effects 0.000 claims description 16
- 238000002091 elastography Methods 0.000 claims description 15
- 230000008859 change Effects 0.000 claims description 6
- 230000004660 morphological change Effects 0.000 claims description 6
- 238000010521 absorption reaction Methods 0.000 claims description 4
- 230000005855 radiation Effects 0.000 claims description 4
- 125000004432 carbon atom Chemical group C* 0.000 claims description 3
- 230000021615 conjugation Effects 0.000 claims description 3
- 230000010287 polarization Effects 0.000 claims description 3
- 210000001519 tissue Anatomy 0.000 description 44
- 238000012937 correction Methods 0.000 description 19
- 238000004364 calculation method Methods 0.000 description 6
- 206010028980 Neoplasm Diseases 0.000 description 5
- 235000019197 fats Nutrition 0.000 description 5
- 238000007689 inspection Methods 0.000 description 5
- 230000003902 lesion Effects 0.000 description 5
- 210000000481 breast Anatomy 0.000 description 4
- 238000001514 detection method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 201000010099 disease Diseases 0.000 description 4
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000036548 skin texture Effects 0.000 description 4
- 239000003981 vehicle Substances 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 206010006187 Breast cancer Diseases 0.000 description 3
- 208000026310 Breast neoplasm Diseases 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 238000012984 biological imaging Methods 0.000 description 3
- 210000000056 organ Anatomy 0.000 description 3
- 230000036285 pathological change Effects 0.000 description 3
- 231100000915 pathological change Toxicity 0.000 description 3
- 208000006011 Stroke Diseases 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 210000004369 blood Anatomy 0.000 description 2
- 239000008280 blood Substances 0.000 description 2
- 210000000988 bone and bone Anatomy 0.000 description 2
- 201000011510 cancer Diseases 0.000 description 2
- 238000012790 confirmation Methods 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 230000004927 fusion Effects 0.000 description 2
- 210000004907 gland Anatomy 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000001575 pathological effect Effects 0.000 description 2
- 230000003252 repetitive effect Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 239000013589 supplement Substances 0.000 description 2
- 206010008190 Cerebrovascular accident Diseases 0.000 description 1
- 206010030113 Oedema Diseases 0.000 description 1
- 208000001132 Osteoporosis Diseases 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 239000012620 biological material Substances 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 230000037237 body shape Effects 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 230000002490 cerebral effect Effects 0.000 description 1
- 235000019864 coconut oil Nutrition 0.000 description 1
- 239000003240 coconut oil Substances 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000013399 early diagnosis Methods 0.000 description 1
- 230000005672 electromagnetic field Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000007499 fusion processing Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000003702 image correction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 208000032839 leukemia Diseases 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 238000002610 neuroimaging Methods 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 239000013535 sea water Substances 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012285 ultrasound imaging Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/0507—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves using microwaves or terahertz waves
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0825—Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of the breast, e.g. mammography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
- A61B8/4416—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device related to combined acquisition of different diagnostic modalities, e.g. combination of ultrasound and X-ray acquisitions
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
- A61B8/5238—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image
- A61B8/5261—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image combining images from different diagnostic modalities, e.g. ultrasound and X-ray
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Veterinary Medicine (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Radar Systems Or Details Thereof (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明涉及一种全息微波弹性成像系统和方法，系统包括微波信号收发模块、超声波信号收发模块和信号处理模块，微波信号收发模块不间断地朝目标区域发射单一频率的微波信号作为入射电场，同时超声波信号收发模块不间断地朝目标区域发射单一频率的超声波信号作为入射超声波场；微波信号穿透位于目标区域的目标物后至少部分微波信号被目标物内部的不同部分反射而形成散射电场并由微波接收天线对散射电场进行探测得到散射电场回波信号；超声波信号收发模块发射超声波信号到目标物从而施加一个外力并接收反射的超声波回波信号；信号处理模块对超声波回波信号和散射电场回波信号进行处理以构建目标物的二维图像。本发明的成像速度快且成像质量高。
Description
技术领域
本发明涉及全息成像和弹性成像技术领域，尤其涉及一种全息微波弹性成像系统及其成像方法。
背景技术
微波成像是指以微波作为信息载体的一种成像手段，实质属于电磁逆散射问题。由于它既用被成像目标散射的幅度信息，也用它的相位信息，因此也称为微波全息成像或者全息微波成像。
目前，微波生物成像因其无创、无有害辐射、价格低廉等优势，为疾病如乳腺癌、脑中风等疾病早期诊断带来新机遇。生物组织的电特性与生理系统的健康状况息息相关，当乳腺癌、白血病等疾病发生时，生物组织的介电性发生显著变化，介电性常用介电常数表示。不同类型生物组织的介电性差异明显，例如骨骼、脂肪、肺以及外层皮肤这些低含水量的组织比例如肌肉、血液和大脑这些高含水量的组织的介电常数低，而典型病变组织，例如恶性肿瘤与其正常组织之间的介电常数差别大于为10％～20％，这些差异为磁感应生物成像检测活体生物组织的生理病理状态提供了可行的物理基础。
弹性成像就是把获取的生物体材料的弹性信息转换成医生习惯的可见光图像，从而让医生能够通过可见光图像判别组织的材料力学特性，进而根据组织的软硬情况判断相应组织或器官可能发生的病理改变以及其位置、形状和大小。具体地，根据各种不同生物组织的弹性系数不同，在加外力或交变振动后其应变或者说形态改变亦不同。在相同外力作用下，弹性系数大的，引起的应变比较小；反之，弹性系数较小的，相应的应变比较大。也就是比较柔软的正常组织变形超过坚硬的肿瘤组织。全息微波弹性成像即利用肿瘤或其他病变区域与周围正常组织间弹性系数的不同和介电属性的不同，在微波和超声波激励作用下，生物组织产生应变大小的不同，以彩色编码显示，来判别病变组织的弹性大小、介电常数分布、电导率分布、温度分布和血液含氧量中的至少一种重要特征，从而推断某些病变的可能性。
过去的二十年，大量的工作主要围绕微波成像应技术在生物成像和诊断方面的应用进行，如脑部成像检测脑卒中和脑水肿、乳房成像检测乳腺癌和骨骼成像检测骨质疏松等。
近年来，有将微波和超声波进行结合成像的技术出现，例如，公开号为CN105708492A的中国专利文献公开了一种融合B超成像和微波成像的方法及系统，通过将所述B超图像和所述微波图像中一幅图像的像素点映射到另一幅图像中，使两幅图像的相关像素点在空间位置上达到一致完成该两幅图像的融合。在诊断早期癌症患者的过程中可以通过使用不同的成像模态进行优势互补，图像融合的潜力在于综合应用这些成像设备所得信息，可以准确地确定病变体的空间位置大小、几何形状及它与周围生物组织之间的空间关系，从而及时高效地诊断疾病，也可以用在手术计划的制定、病理变化的跟踪、治疗效果的评价等方面。但是，该专利进行图像融合所需的计算量大、成像速度慢且成像质量欠佳。
发明内容
针对现有技术之不足，本发明提供了一种全息微波弹性成像方法主要是利用生物组织的介电属性和弹性信息帮助疾病的诊断。其基本原理为：根据各种不同生物组织的弹性系数不同，在加外力或交变振动后其应变亦不同，主要为形态改变不同，同时根据各种不同生物组织的介电属性不同，在微波照射下，目标生物组织发射或散射不同，通过分析处理形态改变和散射场从而重构目标的外观像或介质目标内部的结构成像，也可以是空间电磁场分布的直观显示，相比于融合成像，本发明的成像速度更快，成像质量更好，本发明可用于生物医学成像、疾病诊断、异物探测等领域。
根据一个优选实施方式，一种全息微波弹性成像系统，包括微波信号收发模块、超声波信号收发模块和信号处理模块，所述微波信号收发模块包括至少一个微波发射天线和布置在同一平面上的至少三个微波接收天线，其中，所述微波信号收发模块不间断地朝目标区域发射单一频率的微波信号作为入射电场，同时所述超声波信号收发模块不间断地朝目标区域发射单一频率的超声波信号作为入射超声波场；所述微波信号穿透位于所述目标区域的目标物后至少部分微波信号被所述目标物内部的具有不同介电常数的若干部分反射而形成散射电场并由所述至少三个微波接收天线对所述散射电场进行探测得到散射电场回波信号；所述超声波信号收发模块发射超声波信号到所述目标区域的目标物从而对目标物施加一个外力并接收目标物反射的超声波回波信号；所述信号处理模块对超声波回波信号和散射电场回波信号进行处理以构建目标物的二维图像。
根据一个优选实施方式，所述信号处理模块对超声波回波信号和散射电场回波信号进行处理以构建目标物的二维图像的处理包括：信号处理模块对所述至少三个微波接收天线探测到的散射电场回波信号进行对比获得的可见度函数结合超声波信号收发模块接收到的超声波回波信号通过信号处理方式用于构建目标物的二维图像。
根据一个优选实施方式，所述信号处理模块对所述至少三个微波接收天线探测到的散射电场回波信号进行对比获得的可见度函数结合超声波信号收发模块接收到的超声波回波信号通过信号处理方式用于构建目标物的二维图像的处理包括：对所有微波接收天线中任意两个微波接收天线探测到的散射电场信号进行对比得到可见度函数；根据对比得到的可见度函数获得能够反映目标物的电磁属性分布信息；对所述超声传感器接收的施加外力前和施加外力后目标物的超声波回波信号进行对比得到形态变化信号；根据对比得到的形态变化信号获得目标物内部位移以及振动强度分布信息；基于电磁属性分布信息、内部位移和振动强度分布信息通过傅里叶逆变换处理构建目标物的二维图像。
根据一个优选实施方式，所述信号处理模块对超声波回波信号和散射电场回波信号进行处理以构建目标物的二维图像的处理包括：基于微波接收天线的分布排列建立目标物的电磁属性和散射电场之间的非线性观测模型；至少部分地基于所述非线性观测模型实现对散射回波信号进行两两对比。
根据一个优选实施方式，所述非线性观测模型包括内部场效应模型和外部场效应模型；其中，内部场效应模型描述为： 其中，表示入射电场，表示入射电场到目标物的位置矢量，为场源点到目标物内一点的位置矢量，f1表示微波工作频率，t表示时间，c1表示光速，为TE模式波振幅，AN和BB分别为微波发射天线的窄孔径尺寸和宽孔径尺寸，表示天线辐射场，表示极化矢量；外部场效应模型描述为： 其中，表示散射电场，表示从场源点到场域内任一点的位置矢量，k0表示自由空间的波数，εr表示目标物的介电常数，εb表示背景和/或媒介质的介电常数，表示总场强，G表示格林函数， dV表示体积分，V表示目标物体积,I表示目标物可见度强度；结合内部场效应模型和外部场效应模型得到用于定义目标物的可见度强度的非线性观测模型描述为：
上式中，表示单位矢量，k0表示自由空间的波数，εr表示目标物的介电常数，εb表示背景和/或媒介质的介电常数，表示总场强， *表示复数共轭。
根据一个优选实施方式，所述至少部分地基于所述非线性观测模型实现对散射回波信号进行两两对比的处理包括：通过计算同一平面上所述任意两个微波接收天线的散射电场信号以计算任意两个微波接收天线接收到的可见度函数：其中,表示位于的两个微波接收天线的可见度函数包含相位延迟和/或振幅差异信息，表示目标区域隐藏物中任意点到第i个微波接收天线的距离矢量,表示目标区域隐藏物中任意点到第j个微波接收天线的距离矢量，表示位于处的微波接收天线探测到的散射电场，表示位于处的微波接收天线探测到的散射电场的共轭，*表示复数共轭，＜＞表示平均时间；依次计算任意两个微波接收天线的可见度函数之和，得到所有微波接收天线的总可见度函数：
Υ表示所有微波接收天线的总可见度函数，N表示微波接收天线的总数量，表示位于的两个微波接收天线的可见度函数，其包含相位延迟、振幅差异信息、内部位移和内部组织的振动强度。
根据一个优选实施方式，所述基于介电属性分布的幅值、相位和内部组织的振动强度的信息通过傅里叶逆变换处理方式构建目标物的二维图像的处理包括：
通过对任意两个微波接收天线探测到的可见度函数信号进行二维傅里叶逆变换，得到目标物的二维重构图像：
其中,其中,表示位于的两个微波接收天线的可见度函数，其包含相位延迟、振幅差异信息、内部位移和内部组织的振动强度，l＝sinθcosφ,m＝sinθsinφ,分别为沿x,y,z轴正时空方向的单位矢量,分别为任一微波接收天线Ai在直角坐标系中沿x,y,z轴的位置,分别为任一微波接收天线Aj在直角坐标系中沿x,y,z轴的位置,θ是原点o和空间任一点P的连线与正向z轴的夹角，φ为xoz平面与通过空间任一点P的半平面之间的夹角，若P点在z轴上则φ角是不确定的。表示背景/媒介质的波长。
根据一个优选实施方式，超声波信号收发模块对目标物施加的力表示为：其中，α表示吸收常数，c表示超声波在目标物组织中的速度，ρ表示目标物密度，f为超声波信号收发模块的工作频率，t表示时间，P0表示超声波信号收发模块的压力峰值P为超声波信号收发模块的功率，a和b分别表示超声波信号收发模块的宽带和窄带；超声波信号收发模块探测到的目标物内部组织的内部位移表示为：X(t)＝X0cos(2πft+Phase)；
其中，X0表示最大位移，其取决于局部组织的机械特性和施加的超声波力的强度，Phase表示相位变化，其取决于总路径长度；病灶组织内部的振动强度表示为：其中，εr表示目标物的介电常数，εb表示背景和/或媒介质的介电常数，ε0表示自由空间的介电常数，表示入射场。
根据一个优选实施方式，所述同一平面上的至少三个微波接收天线以非规则排列的方式布置在该平面上。
根据一个优选实施方式，一种全息微波弹性成像方法，包括：配置和/或使用一系统用于成像，所述系统包括微波信号收发模块、超声波信号收发模块和信号处理模块，所述微波信号收发模块包括至少一个微波发射天线和布置在同一平面上的至少三个微波接收天线，其中，所述微波信号收发模块不间断地朝目标区域发射单一频率的微波信号作为入射电场，同时所述超声波信号收发模块不间断地朝目标区域发射单一频率的超声波信号作为入射超声波场；所述微波信号穿透位于所述目标区域的目标物后至少部分微波信号被所述目标物内部的具有不同介电常数的若干部分反射而形成散射电场并由所述至少三个微波接收天线对所述散射电场进行探测得到散射电场回波信号；所述超声波信号收发模块发射超声波信号到所述目标区域的目标物从而对目标物施加一个外力并接收目标物反射的超声波回波信号；所述信号处理模块对超声波回波信号和散射电场回波信号进行处理以构建目标物的二维图像。
附图说明
图1是本发明的系统的一个优选实施方式的简化结构图；
图2是本发明的系统的一个优选实施方式的模块连接图；
图3是系统的其中两个微波接收天线(12)的几何排列示意图；
图4是乳房模型图；
图5是乳房模型的第一重构图，表示针对目标物的介电属性探测到的信号强度；和
图6是乳房模型的第二重构图，表示针对目标物的电导率属性探测到的信号强度。
附图标记列表
10：微波信号收发模块 11：微波发射天线
12：微波接收天线 13：微波控制电路
20：超声波信号收发模块 21：超声传感器
22：超声波控制电路 30：计算设备
31：信号处理模块 32：控制模块
40：图像显示模块 51：纤维淋巴腺
52：人体组织 53：肿瘤
54：皮肤 55：脂肪
具体实施方式
下面结合附图1、2、3、4、5和6进行详细说明。
实施例1
本实施例公开了一种系统，该系统适于执行本发明记载的各个方法步骤，以达到预期的技术效果。
根据一个优选实施方式，一种全息微波弹性成像系统，可以包括微波信号收发模块10、超声波信号收发模块20和信号处理模块31中的至少一个。微波信号收发模块10可以包括至少一个微波发射天线11和在同一平面上布置的至少三个微波接收天线12。微波信号收发模块10可以不间断地朝目标区域发射单一频率的微波信号作为入射电场，同时超声波信号收发模块20可以不间断地朝目标区域发射单一频率的超声波信号作为入射超声波场。微波信号穿透位于目标区域的目标物后至少部分微波信号被目标物内部的具有不同介电常数的若干部分反射而形成散射电场并由至少三个微波接收天线12对散射电场进行探测得到散射电场回波信号。超声波信号收发模块20可以发射超声波信号到目标区域的目标物从而施加一个外力并接收目标物反射的超声波回波信号。信号处理模块31可以对超声波回波信号和至少三个微波接收天线12中任意两个微波接收天线12探测到的散射电场回波信号进行处理以构建目标物的二维图像。优选地，微波信号收发模块10和超声波信号收发模块20设置于目标物的同侧或异侧。优选地，在同一平面上布置的至少三个微波接收天线12使得各个微波接收天线12距离与目标物的垂直距离均相等，由此能降低计算难度。优选地，至少三个微波接收天线12在同一平面上呈非均匀排布，由此提高成像质量。
根据一个优选实施方式，表1和表2给出了一些参考的参数。
表1模拟参数
表2超声波参数
| 衰减 | 速度 | |
| 脂肪 | 0.34dB/cm/MHz | 1479m/s |
| 肿瘤 | 0.79dB/cm/MHz | 1550m/s |
| 腺体 | 1.5dB/cm/MHz | 1553m/s |
优选地，信号处理模块31可以对至少三个微波接收天线12探测到的散射电场回波信号进行对比获得的可见度函数结合超声波信号收发模块20接收到的超声波回波信号通过信号处理方式用于构建目标物的二维图像。优选地，超声波信号收发模块20包括至少一个超声传感器21。超声传感器21可以设于超声波探头内。
根据一个优选实施方式，信号处理模块31可以对超声波回波信号和至少三个微波接收天线12中任意两个微波接收天线12探测到的散射电场回波信号进行处理以构建目标物的二维图像的处理可以包括：对所有微波接收天线12中任意两个微波接收天线12探测到的散射电场信号进行对比得到可见度函数；根据对比得到的可见度函数获得能够反映目标物的电磁属性分布信息；对超声传感器21接收的施加外力前和施加外力后目标物的超声波回波信号进行对比得到形态变化信号；根据对比得到的形态变化信号获得目标物内部位移以及振动强度分布信息；和/或基于电磁属性分布信息、内部位移和振动强度分布信息通过傅里叶逆变换处理构建目标物的二维图像。优选地，电磁属性分布信息可以包括相位延迟和/或振幅差异信息。
根据一个优选实施方式，信号处理模块31对超声波回波信号和至少三个微波接收天线12中任意两个微波接收天线12探测到的散射电场回波信号进行处理以构建目标物的二维图像的处理可以包括：基于微波接收天线12的分布排列建立目标物的电磁属性和散射电场之间的非线性观测模型；和/或至少部分地基于非线性观测模型实现对散射回波信号和/或散射电场进行两两对比。
根据一个优选实施方式，用于定义目标物的可见度强度的非线性观测模型可以描述为：
上式中，表示单位矢量，k0表示自由空间的波数，εr表示目标物的介电常数，εb表示背景和/或媒介质的介电常数，表示总场强， *表示复数共轭。
非线性观测模型包括内部场效应模型和外部场效应模型；
其中，内部场效应模型描述为：
其中，表示入射电场，表示入射电场到目标物的位置矢量，为场源点到目标物内一点的位置矢量，f1表示微波工作频率，t表示时间，c1表示光速，为TE模式波振幅，TE模式波优选为TE10模式波，AN和BB分别为微波发射天线11的窄孔径尺寸和宽孔径尺寸，表示天线辐射场，表示极化矢量。
外部场效应模型可以描述为：
其中，表示散射电场，表示从场源点到场域内任一点的位置矢量，k0表示自由空间的波数，εr表示目标物的介电常数，εb表示背景和/或媒介质的介电常数，表示总场强， G表示格林函数， dV表示体积分，V表示目标物体积,I表示目标物可见度强度。
优选地，超声波衰减主要是由吸收引起的，超声传感器21对目标物施加的力可以表示为：
其中，α表示吸收常数，c表示超声波在目标物组织中的速度，ρ表示目标物密度，f为超声传感器21的工作频率，P0表示超声传感器21的压力峰值P为超声传感器21的功率，a和b分别表示超声传感器21的宽带和窄带。
优选地，超声传感器21探测到的目标物内部组织的内部位移可表示为：
X(t)＝X0cos(2πft+Phase)
其中，X0表示最大位移，其取决于局部组织的机械特性和施加的超声波力的强度，Phase表示相位变化，其取决于总路径长度。
优选地，如果目标物足够小，病灶组织内部的振动强度可以表示为：
其中，εr表示目标物的介电常数，εb表示背景和/或媒介质的介电常数，ε0表示自由空间的介电常数，表示入射场。
根据一个优选实施方式，至少部分地基于非线性观测模型实现对散射电场和/或散射回波信号进行两两对比的处理可以包括：
通过计算同一平面上任意两个微波接收天线12的散射电场信号以计算任意两个微波接收天线12接收到的可见度函数：
其中,表示位于的两个微波接收天线12的可见度函数，其包含相位延迟、振幅差异信息、内部位移和内部组织的振动强度，表示目标区域隐藏物中任意点到第i个微波接收天线12的距离矢量,表示目标区域隐藏物中任意点到第j个微波接收天线12的距离矢量，表示位于处的微波接收天线12探测到的散射电场，表示位于处的微波接收天线12探测到的散射电场的共轭，*表示复数共轭，＜＞表示平均时间。
优选地，依次计算任意两个微波接收天线12的可见度函数之和，得到所有微波接收天线12的总可见度函数：
γ表示所有微波接收天线12的总可见度函数，N表示微波接收天线12的总数量，表示位于的两个微波接收天线12的可见度函数，其包含相位延迟、振幅差异信息、内部位移和内部组织的振动强度。
根据一个优选实施方式，基于介电属性分布的幅值和相位以及内部组织的振动强度的信息通过傅里叶逆变换处理方式构建目标物的二维图像的处理可以包括：
通过对任意两个微波接收天线12探测到的可见度函数信号进行二维傅里叶逆变换，得到目标物的二维重构图像：
其中,其中,表示位于的两个微波接收天线12的可见度函数，其包含相位延迟、振幅差异信息、内部位移和内部组织的振动强度，l＝sinθcosφ,m＝sinθsinφ,分别为沿x,y,z轴正时空方向的单位矢量,分别为任一微波接收天线12Ai在直角坐标系中沿x,y,z轴的位置,分别为任一微波接收天线12Aj在直角坐标系中沿x,y,z轴的位置,θ是原点o和空间任一点P的连线与正向z轴的夹角，φ为xoz平面与通过空间任一点P的半平面之间的夹角，若P点在z轴上则φ角是不确定的。λb表示背景/媒介质的波长。
根据一个优选实施方式，本发明提供了一种全息微波弹性成像系统，其可以包括微波信号收发模块10、超声波信号收发模块20、控制模块32、信号处理模块31和图像显示模块40中的至少一个。优选地，控制模块32和信号处理模块31可以设置在一个设备中，例如设置在同一个计算设备30中。计算设备30可以是服务器、电脑、工作站和手机中的至少一种。
优选地，微波信号收发模块10可以包括微波发射天线11、微波接收天线12和微波控制电路13。微波控制电路13可以连接于微波发射天线11、微波接收天线12和微波控制电路13。超声波信号收发模块可以包括彼此连接的超声传感器21和超声波控制电路22。
优选地，微波信号收发模块10可以包括微波信号发生器、微波信号发射器、微波信号接收器和多通道控制开关电路板，超声波信号收发模块20包括超声信号发生器、超声信号接收器和单通道控制开关电路板，控制模块32和微波信号收发模块10中的多通道控制开关电路板以及超声波信号收发模块20中的单通道控制开关电路板相连接。优选地，多通道控制开关电路板分别与微波信号发生器和微波信号接收器连接，单通道控制开关电路板分别与超声波信号发生器和超声波信号接收器连接。
优选地，控制模块32可以控制微波信号收发模块10中的微波信号发生器产生单频率的微波信号。微波信号以交变电流的形式施加在微波信号发射器中，交变电流产生交变磁场，微波信号接收器对交变磁场进行接收。微波信号接收器将接收到的信号通过控制模块32传输至信号处理模块31；同时，控制模块32控制超声波信号收发模块20中的超声信号发生器产生单频率的超声信号并传输到目标物表面对其施加外力，超声信号接收器接收目标物对外力的回波信号并将接收到的信号通过控制模块32传输至信号处理模块31。信号处理模块31对接收到的信号进行数据处理分析通过傅里叶逆变换获得目标物的二维重构图像，并将二维重构图像传输至图像显示模块40进行显示。
优选地，上述实施例中，微波信号发生器可采用网络分析仪。微波信号发射器和微波信号接收器可以由微波天线实现。超声波信号发射器和超声波信号接收器可以由超声传感器21实现；
优选地，微波信号发射器采用作为激励天线的微波发射天线11，作为微波发射天线11的数量为NT，NT为自然数且NT≥1。作为微波发射天线11的天线可以排列成任意形状且各天线距目标物的垂直距离均相同。微波信号接收器采用作为信号接收天线的微波接收天线12，作为微波接收天线12的数量为NR，NR为自然数且NR≥3。作为微波接收天线12可以排列成均匀和/或非均匀形状且各微波接收天线12距目标物的距离或高度均相同。尤其优选地，微波接收天线12排列成非均匀形状。
优选地，例如，微波天线的数量为16个，以目标物为中心，16个微波天线围绕着目标物均匀排列成正方形。每个微波天线既作为微波发射天线11传输微波信号，又作为微波接收天线12探测目标物内部及其周围的电场变化以及介电属性和电导率的分布状态。
优选地，为减少信号耦合，提高检测灵敏度，目标生物与微波天线之间以及各微波天线之间的间隙内都填充介电属性与正常脂肪相似的媒介物质，例如海水、盐水和椰子油中的至少一种。
优选地，微波信号发生器可以采用网络分析仪，其可以产生频率范围为1GHz-10GHz的微波信号。
优选地，本发明的微波信号收发模块10的工作频率为单频率，其最佳工作频率范围为1GHz-10GHz。比如，可以采用2GHz、3GHz、4GHz、5GHz、6GHz、7GHz、8GHz或9GHz，尤其优选8～9GHz。
优选地，本发明的超声波信号收发模块20的工作频率为单频率，其最佳工作频率范围为1MHz-10MHz。比如，可以采用2MHz、3MHz、4MHz、5MHz、6MHz、7MHz、8MHz或9MHz，尤其优选6～7Mhz。
优选地，微波天线例如可采用开口波导和/或贴片天线，微波天线的大小受工作频率的限制。
优选地，作为微波发射天线11的发射天线与作为测量天线的微波接收天线12之间的位置关系为：作为微波发射天线11的发射天线与作为测量天线的微波接收天线12均位于目标物的同侧，位于相同高度或不同高度；或作为微波发射天线11的发射天线与作为测量天线的微波接收天线12相互平行；或作为微波发射天线11的发射天线与作为测量天线的微波接收天线12阵列重合；或作为微波发射天线11的发射天线与作为测量天线的微波接收天线12阵列均与目标物体一定的角度。
实施例2
本实施例还公开了一种系统，该系统适于执行本发明记载的各个方法步骤，以达到预期的技术效果。本实施例可以是对实施例1的进一步改进和/或补充，重复的内容不再赘述。在不造成冲突或者矛盾的情况下，其他实施例的优选实施方式的整体和/或部分内容可以作为本实施例的补充。
根据一个优选实施方式，超声波信号收发模块20还接收由微波信号收发模块10发射的微波信号照射到目标物上产生的热声信号，在信号处理模块31对至少三个微波接收天线12探测到的散射电场回波信号进行对比获得的可见度函数结合超声波信号收发模块20接收到的超声波回波信号通过信号处理方式用于构建目标物的二维图像后，信号处理模块31至少部分地基于热声信号对构建的目标物的二维图像进行修正，并将经过修正的二维图像发送到图像显示模块40进行显示。通常，热声信号是利用图像融合处理方式进行成像，其成像质量欠佳且速度慢，而本发明是利用主动发射和探测的超声波信号和微波信号来成像，然后利用微波信号照射到目标物产生的热声信号来对图像进行修正，其成像速度更快且成像质量更高。
根据一个优选实施方式，信号处理模块31基于热声信号、超声波回波信号和散射电场回波信号中的至少一个分析目标物的运动并由此对构建的目标物的二维图像进行修正。检测对象内部组织的运动会对成像质量造成影响，例如，心跳、血液流动或者呼吸，会对弹性成像的结果造成影响，由此，本发明利用采集到的数据分析目标物的运动，并对图像进行运动修正，可以提高成像的准确性和质量。
优选地，在信号处理模块31基于热声信号对构建的目标物的二维图像进行修正前，控制模块32对目标物的二维图像的成像质量进行评估，仅在成像质量低于质量阈值时才对构建的目标物的二维图像进行修正。而在成像质量高于或等于质量阈值时，直接将目标物的二维图像发送到图像显示模块40进行显示。
根据一个优选实施方式，信号处理模块31根据热声信号分析目标物从第一时间到第二时间的运动以确定第一修正参数，信号处理模块31利用第一修正参数对构建的目标物的二维图像进行修正。
根据一个优选实施方式，信号处理模块31根据超声波回波信号分析目标物从第一时间到第二时间的运动以确定第二修正参数，信号处理模块31基于第一修正参数和第二修正参数采用加权法计算第三修正参数，信号处理模块31响应于用户请求而基于第三修正参数对构建的目标物的二维图像进行修正，并将经过修正的二维图像发送到图像显示模块40进行显示。
根据一个优选实施方式，信号处理模块31根据散射电场回波信号分析目标物从第一时间到第二时间的运动以确定第四修正参数，信号处理模块31基于第三修正参数和第四修正参数采用加权法计算第五修正参数，信号处理模块31响应于用户请求而基于第五修正参数对构建的目标物的二维图像进行修正，并将经过修正的二维图像发送到图像显示模块40进行显示。在对图像进行修正时，基于第五修正参数的计算量大于基于第三修正参数的计算量，基于第三修正参数的计算量大于基于第一修正参数的计算量，计算量越大导致成像速度越慢，因此，实际操作过程中，基于第五修正参数或者第三修正参数的图像修正参数由用户决定是否启用，比如，设置两个专用按键，更适于操作者根据实际的成像情况做出快速地调整，以更好地满足实际成像需求。
根据一个优选实施方式，微波信号收发模块10包括在同一平面上呈非均匀排布的至少四个微波接收天线12；
在成像前，控制模块32获取待检查对象的身份标识并根据身份标识查询待检测对象的历史检查数据；
在未查询到历史检查数据的情况下，控制模块32从至少四个微波接收天线12中以非重复随机方式选择至少三个微波接收天线12同时作为成像过程所需的微波发射天线11，并在预设时间内至少部分基于各个选择的硬件条件实现对相应的目标区域的成像，然后控制模块32基于对成像质量的排序输出成像质量最高的图像和确定相应的目标区域的推荐硬件条件；
在查询到历史检查数据的情况下，控制模块32根据当前的目标区域的推荐硬件条件配置微波发射天线11，基于配置实现对相应的目标区域的成像。本发明采用此方式至少能够实现以下有益技术效果：第一，每个待检测对象的自身特性存在，例如身体形态、脂肪层厚度、组织含水量差异和器官大小等等，由此，其最佳的成像硬件条件有差异，因此，本发明在检测时，通过非重复随机方式选择最优的配置方式，实现最优的成像质量；第二，保存对象的历史检查数据和相应的目标区域的推荐硬件条件，再该对象再次检查时，快速采用推荐硬件条件进行配置，减少配置时间，提高成像效率。优选地，同时作为微波成像过程中所需的微波发射天线11是指同时作为微波发射天线11和微波接收天线12，而未被选择的微波接收天线12仅作为微波接收天线12。例如，假设系统中共有四个微波发射天线11，以非重复随机方式选择至少三个微波接收天线12则共有七种选择，在两秒内用这七种硬件条件实现对相应的目标区域的七组成像，然后对七组成像的成像质量进行排序，选择成像质量最高的一组图像和其对于的硬件条件作为该目标区域的推荐硬件条件，作为以后同一个对象再次检查时配置系统的依据。优选地，身份标识例如可以是身份证号、姓名、就诊卡号、指纹和虹膜中的至少一种。优选地，控制模块32基于峰值信噪比、结构相似度和模糊系数中的至少两种采用加权法对成像质量进行评估并由此进行排序。
根据一个优选实施方式，系统还包括图像采集装置，图像采集装置与微波信号收发模块10同步运动采集皮肤纹理，并将皮肤纹理与对象的相应的目标区域关联地保存到历史检查数据；控制模块32基于对目标物的二维图像的图像识别初步判断当前的目标区域，随后控制模块32至少部分地基于图像采集装置采集到的实时皮肤纹理从历史检查数据中选取与初步判断的目标区域相对应的皮肤纹理对当前的目标区域进行二次确认，并在二次确认通过的情况下，控制模块32才根据当前的目标区域的推荐硬件条件配置微波发射天线11。本发明采用此方式至少能够实现以下有益技术效果：大幅减少确定目标区域计算量，减少系统配置时间。
实施例3
本实施例还公开了一种方法，该方法可以由本发明的系统和/或其他可替代的零部件实现。比如，通过使用本发明的系统中的各个零部件实现本发明的方法。
根据一个优选实施方式，一种全息微波弹性成像方法，可以包括：配置和/或使用一系统用于成像。
根据一个优选实施方式，一种全息微波弹性成像方法，该方法可以包括：利用将微波信号收发模块10向目标物发送微波信号并探测目标物反射微波信号形成的散射电场获得散射电场回波信号，利用超声波信号收发模块20向目标物发送入射超声波并探测目标物反射微波信号形成的散射超声波场获得超声波回波信号，将接收到的散射电场回波信号和超声波回波信号之和进行两两对比得到的可见信号通过傅里叶信号处理方式以构建目标物的二维图像。该方法能够利用生物组织的机械特性和电磁属性获得更清晰的目标生物组织图像，提高成像质量和目标物检测灵敏度。本发明可应用于无损检测、医学成像和目标探测等领域。
如本文所用的词语“模块”描述任一种硬件、软件或软硬件组合，其能够执行与“模块”相关联的功能。
需要注意的是，上述具体实施例是示例性的，本领域技术人员可以在本发明公开内容的启发下想出各种解决方案，而这些解决方案也都属于本发明的公开范围并落入本发明的保护范围之内。本领域技术人员应该明白，本发明说明书及其附图均为说明性而并非构成对权利要求的限制。本发明的保护范围由权利要求及其等同物限定。
Claims (10)
1.一种全息微波弹性成像系统，其特征在于，包括微波信号收发模块(10)、超声波信号收发模块(20)和信号处理模块(31)，所述微波信号收发模块(10)包括至少一个微波发射天线(11)和布置在同一平面上的至少三个微波接收天线(12)，其中，所述微波信号收发模块(10)不间断地朝目标区域发射单一频率的微波信号作为入射电场，同时所述超声波信号收发模块(20)不间断地朝目标区域发射单一频率的超声波信号作为入射超声波场；
所述微波信号穿透位于所述目标区域的目标物后至少部分微波信号被所述目标物内部的具有不同介电常数的若干部分反射而形成散射电场并由所述至少三个微波接收天线(12)对所述散射电场进行探测得到散射电场回波信号；
所述超声波信号收发模块(20)发射超声波信号到所述目标区域的目标物从而对目标物施加一个外力并接收目标物反射的超声波回波信号；
所述信号处理模块(31)对超声波回波信号和散射电场回波信号进行处理以构建目标物的二维图像。
2.如权利要求1所述的系统，其特征在于，所述信号处理模块(31)对超声波回波信号和散射电场回波信号进行处理以构建目标物的二维图像的处理包括：
信号处理模块(31)对所述至少三个微波接收天线(12)探测到的散射电场回波信号进行对比获得的可见度函数结合超声波信号收发模块(20)接收到的超声波回波信号通过信号处理方式用于构建目标物的二维图像。
3.如前述权利要求之一所述的系统，其特征在于，所述信号处理模块(31)对所述至少三个微波接收天线(12)探测到的散射电场回波信号进行对比获得的可见度函数结合超声波信号收发模块(20)接收到的超声波回波信号通过信号处理方式用于构建目标物的二维图像的处理包括：
对所有微波接收天线(12)中任意两个微波接收天线(12)探测到的散射电场信号进行对比得到可见度函数；
根据对比得到的可见度函数获得能够反映目标物的电磁属性分布信息；
对所述超声传感器(21)接收的施加外力前和施加外力后目标物的超声波回波信号进行对比得到形态变化信号；
根据对比得到的形态变化信号获得目标物内部位移以及振动强度分布信息；
基于电磁属性分布信息、内部位移和振动强度分布信息通过傅里叶逆变换处理构建目标物的二维图像。
4.如前述权利要求之一所述的系统，其特征在于，所述信号处理模块(31)对超声波回波信号和散射电场回波信号进行处理以构建目标物的二维图像的处理包括：
基于微波接收天线(12)的分布排列建立目标物的电磁属性和散射电场之间的非线性观测模型；
至少部分地基于所述非线性观测模型实现对散射回波信号进行两两对比。
5.如前述权利要求之一所述的系统，其特征在于，
所述非线性观测模型包括内部场效应模型和外部场效应模型；
其中，内部场效应模型描述为：
其中，表示入射电场，表示入射电场到目标物的位置矢量，为场源点到目标物内一点的位置矢量，f1表示微波工作频率，t表示时间，c1表示光速，为TE模式波振幅，AN和BB分别为微波发射天线(11)的窄孔径尺寸和宽孔径尺寸，表示天线辐射场，表示极化矢量；
外部场效应模型描述为：
其中，表示散射电场，表示从场源点到场域内任一点的位置矢量，k0表示自由空间的波数，εr表示目标物的介电常数，εb表示背景和/或媒介质的介电常数，表示总场强， G表示格林函数， dV表示体积分，V表示目标物体积,I表示目标物可见度强度；
结合内部场效应模型和外部场效应模型得到用于定义目标物的可见度强度的非线性观测模型描述为：
上式中，表示单位矢量，k0表示自由空间的波数，εr表示目标物的介电常数，εb表示背景和/或媒介质的介电常数，表示总场强， *表示复数共轭。
6.如前述权利要求之一所述的系统，其特征在于，所述至少部分地基于所述非线性观测模型实现对散射回波信号进行两两对比的处理包括：
通过计算同一平面上所述任意两个微波接收天线(12)的散射电场信号以计算任意两个微波接收天线(12)接收到的可见度函数：
其中,表示位于的两个微波接收天线(12)的可见度函数，其包含相位延迟、振幅差异信息、内部位移和内部组织的振动强度，表示目标区域隐藏物中任意点到第i个微波接收天线(12)的距离矢量,表示目标区域隐藏物中任意点到第j个微波接收天线(12)的距离矢量，表示位于处的微波接收天线(12)探测到的散射电场，表示位于处的微波接收天线(12)探测到的散射电场的共轭，*表示复数共轭，<>表示平均时间；
依次计算任意两个微波接收天线(12)的可见度函数之和，得到所有微波接收天线(12)的总可见度函数：
Υ表示所有微波接收天线(12)的总可见度函数，N表示微波接收天线(12)的总数量，表示位于的两个微波接收天线(12)的可见度函数。
7.如前述权利要求之一所述的系统，其特征在于，所述基于介电属性分布的幅值和相位的信息通过傅里叶逆变换处理方式构建目标物的二维图像的处理包括：
通过对任意两个微波接收天线(12)探测到的可见度函数信号进行二维傅里叶逆变换，得到目标物的二维重构图像：
其中,其中,表示位于的两个微波接收天线(12)的可见度函数，其包含相位延迟、振幅差异信息、内部位移和内部组织的振动强度，l＝sinθcosφ,m＝sinθsinφ, 分别为沿x,y,z轴正时空方向的单位矢量,分别为任一微波接收天线(12)Ai在直角坐标系中沿x,y,z轴的位置,分别为任一微波接收天线(12)Aj在直角坐标系中沿x,y,z轴的位置,θ是原点o和空间任一点P的连线与正向z轴的夹角，φ为xoz平面与通过空间任一点P的半平面之间的夹角，若P点在z轴上则φ角是不确定的；λb表示背景和/或媒介质的波长。
8.如前述权利要求之一所述的系统，其特征在于，所述超声波信号收发模块(20)对目标物施加的力表示为：
其中，α表示吸收常数，c表示超声波在目标物组织中的速度，ρ表示目标物密度，f为超声波信号收发模块(20)的工作频率，t表示时间，P0表示超声波信号收发模块(20)的压力峰值P为超声波信号收发模块(20)的功率，a和b分别表示超声波信号收发模块(20)的宽带和窄带；
超声波信号收发模块(20)探测到的目标物内部组织的内部位移表示为：
X(t)＝X0cos(2πft+Phase)；
其中，X0表示最大位移，其取决于局部组织的机械特性和施加的超声波力的强度，Phase表示相位变化，其取决于总路径长度；
内部组织的振动强度表示为：
其中，εr表示目标物的介电常数，εb表示背景和/或媒介质的介电常数，ε0表示自由空间的介电常数，表示入射场。
9.如前述权利要求之一所述的系统，其特征在于，所述同一平面上的至少三个微波接收天线(12)以非规则排列的方式布置在该平面上。
10.一种全息微波弹性成像方法，其特征在于，包括：配置和/或使用一系统用于成像，所述系统包括微波信号收发模块(10)、超声波信号收发模块(20)和信号处理模块(31)，所述微波信号收发模块(10)包括至少一个微波发射天线(11)和布置在同一平面上的至少三个微波接收天线(12)，其中，所述微波信号收发模块(10)不间断地朝目标区域发射单一频率的微波信号作为入射电场，同时所述超声波信号收发模块(20)不间断地朝目标区域发射单一频率的超声波信号作为入射超声波场；
所述微波信号穿透位于所述目标区域的目标物后至少部分微波信号被所述目标物内部的具有不同介电常数的若干部分反射而形成散射电场并由所述至少三个微波接收天线(12)对所述散射电场进行探测得到散射电场回波信号；
所述超声波信号收发模块(20)发射超声波信号到所述目标区域的目标物从而对目标物施加一个外力并接收目标物反射的超声波回波信号；
所述信号处理模块(31)对超声波回波信号和散射电场回波信号进行处理以构建目标物的二维图像。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811061290.5A CN109199381B (zh) | 2018-09-11 | 2018-09-11 | 一种全息微波弹性成像系统及其成像方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811061290.5A CN109199381B (zh) | 2018-09-11 | 2018-09-11 | 一种全息微波弹性成像系统及其成像方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109199381A true CN109199381A (zh) | 2019-01-15 |
| CN109199381B CN109199381B (zh) | 2021-11-02 |
Family
ID=64983549
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811061290.5A Active CN109199381B (zh) | 2018-09-11 | 2018-09-11 | 一种全息微波弹性成像系统及其成像方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109199381B (zh) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110720914A (zh) * | 2019-10-25 | 2020-01-24 | 深圳技术大学 | 基于稀疏采样的全息磁感应胸腔成像方法及成像系统 |
| CN110960216A (zh) * | 2019-10-25 | 2020-04-07 | 深圳技术大学 | 多频全息微波脑成像系统及成像方法 |
| WO2021067265A1 (en) | 2019-10-03 | 2021-04-08 | Endra Life Sciences Inc. | Method and system for determining a parameter of a material of interest |
| WO2021077547A1 (zh) * | 2019-10-25 | 2021-04-29 | 深圳技术大学 | 用于生物组织检测的磁感应分子成像方法及系统 |
| CN112914540A (zh) * | 2021-01-29 | 2021-06-08 | 深圳技术大学 | 基于PXIe总线的全息微波成像系统及其成像方法 |
| CN113820398A (zh) * | 2021-09-30 | 2021-12-21 | 电子科技大学 | 一种极化微波热声成像装置及方法 |
| CN114396864A (zh) * | 2021-08-20 | 2022-04-26 | 上海交通大学 | 微波多维形变及振动测量方法与系统 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090316854A1 (en) * | 2008-06-23 | 2009-12-24 | Ismail Aly M | Multi-modality system for screening, imaging and diagnosis in dense compressive media and method of use thereof |
| CN102283679A (zh) * | 2011-08-04 | 2011-12-21 | 中国科学院深圳先进技术研究院 | 弹性测量的超声成像系统及测量生物组织弹性的方法 |
| CN105708492A (zh) * | 2015-12-31 | 2016-06-29 | 深圳市一体医疗科技有限公司 | 一种融合b超成像和微波成像的方法及系统 |
| CN105997070A (zh) * | 2016-06-15 | 2016-10-12 | 合肥工业大学 | 一种非接触式磁感应成像系统及其成像方法 |
| CN106037816A (zh) * | 2016-05-19 | 2016-10-26 | 深圳大学 | 一种超声粘弹性测量方法及系统 |
| WO2018078315A1 (en) * | 2016-10-27 | 2018-05-03 | Micrima Limited | System and method for combined microwave and ultrasound imaging |
-
2018
- 2018-09-11 CN CN201811061290.5A patent/CN109199381B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090316854A1 (en) * | 2008-06-23 | 2009-12-24 | Ismail Aly M | Multi-modality system for screening, imaging and diagnosis in dense compressive media and method of use thereof |
| CN102283679A (zh) * | 2011-08-04 | 2011-12-21 | 中国科学院深圳先进技术研究院 | 弹性测量的超声成像系统及测量生物组织弹性的方法 |
| CN105708492A (zh) * | 2015-12-31 | 2016-06-29 | 深圳市一体医疗科技有限公司 | 一种融合b超成像和微波成像的方法及系统 |
| CN106037816A (zh) * | 2016-05-19 | 2016-10-26 | 深圳大学 | 一种超声粘弹性测量方法及系统 |
| CN105997070A (zh) * | 2016-06-15 | 2016-10-12 | 合肥工业大学 | 一种非接触式磁感应成像系统及其成像方法 |
| WO2018078315A1 (en) * | 2016-10-27 | 2018-05-03 | Micrima Limited | System and method for combined microwave and ultrasound imaging |
Non-Patent Citations (1)
| Title |
|---|
| LULU WANG ET AL: "Holographic microwave imaging for medical applications", 《J. BIOMEDICAL SCIENCE AND ENGINEERING》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2021067265A1 (en) | 2019-10-03 | 2021-04-08 | Endra Life Sciences Inc. | Method and system for determining a parameter of a material of interest |
| CN110720914A (zh) * | 2019-10-25 | 2020-01-24 | 深圳技术大学 | 基于稀疏采样的全息磁感应胸腔成像方法及成像系统 |
| CN110960216A (zh) * | 2019-10-25 | 2020-04-07 | 深圳技术大学 | 多频全息微波脑成像系统及成像方法 |
| WO2021077547A1 (zh) * | 2019-10-25 | 2021-04-29 | 深圳技术大学 | 用于生物组织检测的磁感应分子成像方法及系统 |
| CN110960216B (zh) * | 2019-10-25 | 2023-04-07 | 深圳技术大学 | 多频全息微波脑成像系统及成像方法 |
| CN112914540A (zh) * | 2021-01-29 | 2021-06-08 | 深圳技术大学 | 基于PXIe总线的全息微波成像系统及其成像方法 |
| CN114396864A (zh) * | 2021-08-20 | 2022-04-26 | 上海交通大学 | 微波多维形变及振动测量方法与系统 |
| CN114396864B (zh) * | 2021-08-20 | 2023-02-28 | 上海交通大学 | 微波多维形变及振动测量方法与系统 |
| CN113820398A (zh) * | 2021-09-30 | 2021-12-21 | 电子科技大学 | 一种极化微波热声成像装置及方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109199381B (zh) | 2021-11-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109199381B (zh) | 一种全息微波弹性成像系统及其成像方法 | |
| JP7018024B2 (ja) | マイクロ波トモグラフィシステム | |
| Nikolova | Microwave imaging for breast cancer | |
| Fear et al. | Microwaves for breast cancer detection? | |
| Ireland et al. | Microwave head imaging for stroke detection | |
| O'Halloran et al. | Quasi-multistatic MIST beamforming for the early detection of breast cancer | |
| CN107569256B (zh) | 基于热膨胀和门控算法测量生物组织温度变化的超声方法 | |
| EP1788940A1 (en) | Synthetic focusing method | |
| Jalilvand et al. | UWB synthetic aperture-based radar system for hemorrhagic head-stroke detection | |
| Godinho et al. | Experimental evaluation of an axillary microwave imaging system to aid breast cancer staging | |
| EP3556278A1 (en) | Scanning device for living objects | |
| Ojaroudi et al. | A novel approach of brain tumor detection using miniaturized high-fidelity UWB slot antenna array | |
| Elahi | Confocal microwave imaging and artifact removal algorithms for the early detection of breast cancer | |
| Bicer et al. | A matching-pursuit based approach for detecting and imaging breast cancer tumor | |
| CN110960216B (zh) | 多频全息微波脑成像系统及成像方法 | |
| Malla et al. | Investigation of breast tumor detection using microwave imaging technique | |
| Solimene et al. | An incoherent radar imaging system for medical applications | |
| Ricci et al. | Beamforming algorithms for UWB radar-based stroke detection: Trade-off performance-complexity | |
| Ricci et al. | Artifact removal algorithms for stroke detection using a multistatic MIST beamforming algorithm | |
| Töpfer et al. | Microwave cancer diagnosis | |
| Fogel et al. | An integrated microwave-ultrasound breast imaging system: Initial phantom results | |
| Cheng et al. | 3-D microwave imaging for breast cancer | |
| Bicer et al. | An experimental study on microwave imaging of breast cancer with the use of tumor phantom | |
| Nath | Simulation of a portable microwave breast cancer detection system and derivation of a mathematical model for the relative electromagnetic fields | |
| CN205729400U (zh) | 一种融合b超成像和雷达成像的装置 |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |