CN105023276A - 一种医用微波热疗天线sar图与肿瘤ct图像的匹配方法 - Google Patents
一种医用微波热疗天线sar图与肿瘤ct图像的匹配方法 Download PDFInfo
- Publication number
- CN105023276A CN105023276A CN201510517718.2A CN201510517718A CN105023276A CN 105023276 A CN105023276 A CN 105023276A CN 201510517718 A CN201510517718 A CN 201510517718A CN 105023276 A CN105023276 A CN 105023276A
- Authority
- CN
- China
- Prior art keywords
- image
- sar
- microwave
- tumour
- antenna
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30096—Tumor; Lesion
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
- Radiation-Therapy Devices (AREA)
Abstract
本发明属于微波热疗技术领域,涉及一种医用微波热疗天线SAR图与肿瘤CT图像的匹配方法。本发明通过肿瘤CT图像和微波SAR图像的数字化处理和检索匹配,可以根据不同形状的肿瘤CT图像检索寻找的相似的微波SAR图,并确定所需微波阵列天线的功率、相位、阵列数量、频率的参数。
Description
技术领域
本发明属于医用成像技术领域,特别是涉及一种医用微波热疗天线SAR图与肿瘤CT图像的匹配方法。
背景技术
微波热疗已逐渐成为治疗肿瘤的一种重要的技术,而微波热疗技术中最关键的环节就是在热疗过程中杀死肿瘤细胞的同时降低对正常细胞的损伤。目前微波热疗中常采用微波热疗天线在体外进行热疗,一般针对的是浅表层的肿瘤,位置大约为皮下2~3cm,如皮肤癌、乳腺肿瘤、软组织肉瘤等。微波能量作用于肿瘤组织时,其在肿瘤组织中产生的高热效果能抑制肿瘤细胞的DNA和RNA的复制、蛋白质合成和肿瘤细胞的繁殖,能使肿瘤骨架散乱,细胞的许多重要功能受损,直接导致肿瘤细胞的死亡。
但是在微波热疗中,在长时间均匀加热的基础上,由于靠近辐射器的组织温度较高,所以对于除表浅肿瘤外,极有可能损伤到正常组织。另外一种肿瘤只采用一种辐射器,辐射器产生的热量分布范围与肿瘤并不能完全匹配,导致完全杀死肿瘤的概率降低。
现有的电磁仿真软件HFSS(High Frequeney Strueture Simulator)是Ansoft公司开发的可用于任意三维无源器件的高性能电磁场仿真软件,通过HFSS可以得到天线的微波辐射SAR图。国际上多采用比吸收率SAR描述生物组织对电磁场能量的吸收特性。SAR的定义是单位质量生物组织对电磁场能量的吸收率(单位时间内单位体积吸收的电磁能量)。生物组织吸收电磁波能量后产生热量,温度上升,因此SAR图可以近似的看作为生物组织受到微波辐射后的温度分布图。对于不同微波热疗天线或者天线阵列通过调节辐射频率、天线相位、辐射功率等参数可以产生不同的SAR图。
发明内容
本发明的目的在于提供一种医用微波热疗天线SAR图与肿瘤CT图像的匹配方法,该方法通过对微波热疗天线SAR图和肿瘤CT图像的匹配,可以有效地轮廓覆盖不同形状的肿瘤组织,降低对正常组织损伤。
本发明所采用的技术方案是:
本发明采用的微波热疗天线是微波阵列天线。
一种医用微波热疗天线SAR图与肿瘤CT图像的匹配方法,其特征在于,方法步骤如下:
1)建立图像库:在HFSS中进行微波SAR图像的采集,确定微波阵列天线参数,根据微波阵列天线参数确定SAR图像,通过数据处理构成图像库,所述微波阵列天线参数包括:功率、相位、频率和阵列数量;
2)图像匹配搜索:根据肿瘤CT影像,通过数据处理后,在图像库中找到与肿瘤CT影像完全匹配的SAR图像,再确定微波阵列天线参数;若在图像库中找到SAR图像与肿瘤CT影像低于相似度阈值时,则通过位移填充完成匹配,直到SAR图像与肿瘤CT影像符合设定的相似度后,再来确定微波阵列天线的功率、频率、相位和阵列数量参数值。
所述的相似度阈值为多个参数:1.SAR图和肿瘤CT面积相差的百分数(可设定1%-10%),2.SAR图和肿瘤CT周长相差的百分数(可设定1%-20%),3 SAR图和肿瘤CT长宽比相差的百分数(可设定1%-10%)。
所述微波阵列天线功率为0-100W。
所述微波阵列天线频率433MHz、915MHz或2450 MHz。
所述微波阵列天线的阵列数量为1-5个,阵列排列方式有1、1*1、1*3、2*2或2*3。
所述微波阵列天线的相位为0°-360°。
本发明图像数据处理包括下列步骤:1)图像转换:采用阈值将肿瘤CT图像和SAR图像转换为二值图像;
2)噪声消除:为了消除在边界或边界内部的“空洞”,先用填充函数将“空洞”进行填充,然后进行噪声消除,噪声消除采用均值滤波,用滤波函数滤波之后得到底图;
3)多边形近似:取滤波后底图,根据底图调用函数,把边界提取出来,呈现出一个多边形图片,最后通过这个多边形图片进行匹配,所述函数是根据最小周长多边形近似原理确定。
本发明图像匹配:利用特征匹配中的面特征,对肿瘤图像和SAR图进行提取,在两者的面特征即轮廓提取出来后,计算面特征的二值特征,然后根据二值特征设定检索条件把符合条件的SAR图呈现出来,并显示出它所标注的微波阵列天线类型,所述二值特征包括:面积、偏心率、周长、重心、长宽比。
本发明图像库内的每个SAR图上都标注出所对应的微波阵列天线具体参数。
本发明根据肿瘤CT影像,在SAR图像库中找到和肿瘤CT影像相匹配的SAR图像,并根据SAR图像来确定微波阵列天线的功率、频率、相位和阵列数量。
本发明采用的微波阵列天线的体积小,容易搭配,对不规则的肿瘤通过图像匹配拼接填充,有效地轮廓覆盖不同形状的肿瘤组织,达到参数选取的最佳范围,降低对正常组织损伤。
附图说明
图1 经过数据处理的SAR图;
图2 肿瘤图像;
图3经过处理的肿瘤图像;
图4图库中检索出来的最相似的SAR图像
图5图库中检索出来的较相似的SAR图像
图6是本发明的工作流程图;
图7是本发明的图像处理流程图;
图8是本发明的匹配检索流程图。
具体实施方式
现结合附图和具体实施方式对本方法作进一步描述。
实例1:首先进行的是在HFSS中采集微波SAR图像,在不同的参数(如功率20W、相位(0°、90°、180°、270°、360°)、频率(915MHz或2450MHz)和阵列数量(1或2)的微波阵列天线的SAR图像按一定的顺序进行采集,可采集到大量的SAR图,在采集SAR图的同时要对其进行标注,如果SAR图是由功率20W、相位90°、频率915MHz、阵列数量2的天线产生,则在其经过数据处理形成的图像库中也标注出来,只要满足匹配图片的要求,就可确定微波热疗所需天线的参数为功率20W、相位90°、频率915MHz、阵列数量2。功率20W、相位90°、频率915MHz、阵列数量2的天线产生的SAR图经过数据处理后得到如图1。
采用阈值把肿瘤CT图像图2,通过数据处理成为二值图像,为了消除在边界或边界内部的“空洞”,先用填充函数将“空洞”进行填充,然后进行噪声消除,噪声消除采用均值滤波,用滤波函数滤波之后得到底图首先要做的就是取滤波后图像的大小,然后提取边界得到边界,再以边界为根据,调用根据最小周长多边形近似原理编好的函数得到图像边界的顶点和连接边界的顶点,最后把边界提取出来,呈现在我们面前的就是一个多边形图3。最后通过这个多边形来进行匹配。通过对面特征的计算,利用特征匹配中的面特征,对多边形进行提取,在面特征即轮廓提取出来后,计算面特征的二值特征,所述二值特征包括:面积、偏心率、周长、重心、长宽比。图3的面积为47512、周长为881.09、偏心率为0.6963、宽高比为0.8383、重心(行,列)为【285.79 212.33】,通过对SAR图像库的检索出最相近的SAR图像图4的面积为48956、周长为900.17、偏心率为0.7512、宽高比为0.9031,满足预先设定的相似度范围,即SAR图和肿瘤CT面积相差的3%小于预设的10%,SAR图和肿瘤CT周长相差的2%小于预设的10%, SAR图和肿瘤CT长宽比相差的7.7%小于预设的10%。然后根据图像的标注找到所需微波天线的功率20w、相位90°、阵列数量1*1、频率915MHz,最终确定选用这种参数的天线针对图2所示的肿瘤进行热疗。
如若图像库中最相近的图像为图5,其面积为34968、周长为650、偏心率为0.5042、宽高比0.7321,其天线参数为(20W、相位180°、阵列数为1、频率为2450MHz),经搜索比对没有其他的图像与肿瘤图像更相似,则需要对图片进行覆盖填充处理,覆盖填充的过程需要将图5经过平移、放大或缩小处理,使得SAR图能够完全覆盖肿瘤图像,通过计算的结果显示将图5放大1.3倍再进行平移,平移量(x轴)138.4245、(y轴)0.1961,满足预先设定的相似度范围,则可选定此天线,在天线辐射时,需根据平移量和放大倍数来进行相对的移动和提高辐射功率为26W来满足热疗要求。
本发明通过对肿瘤CT图像和微波SAR图像的数字化处理和检索匹配,可以根据不同形状的肿瘤CT图像检索寻找相似的微波SAR图,并确定所需微波阵列天线的功率、相位、阵列数量、频率的参数,对不规则的肿瘤通过图像匹配拼接填充,达到参数选取的最佳范围,实现最佳匹配。
Claims (8)
1.一种医用微波热疗天线SAR图与肿瘤CT图像的匹配方法,其特征在于,方法步骤如下:
1)建立图像库:在HFSS中进行微波SAR图像的采集,确定微波阵列天线参数,根据微波阵列天线参数确定SAR图像,通过数据处理构成图像库,所述微波阵列天线参数包括:功率、相位、频率和阵列数量;
2)图像匹配搜索:根据肿瘤CT影像,通过数据处理后,在图像库中找到与肿瘤CT影像形状相似度高的SAR图像,再确定微波阵列天线参数;若在图像库中找到SAR图像与肿瘤CT影像低于相似度阈值时,则通过位移填充完成匹配,直到SAR图像与肿瘤CT影像的相似度达到设定值后,再来确定微波阵列天线的功率、频率、相位和阵列数量参数值。
2.根据权利要求1所述的一种医用微波热疗天线SAR图与肿瘤CT图像的匹配方法,其特征在于,微波阵列天线功率为0-100W。
3.根据权利要求1或2所述的一种医用微波热疗天线SAR图与肿瘤CT图像的匹配方法,其特征在于,微波阵列天线频率为433MHz、915MHz或2450 MHz。
4.根据权利要求3所述的一种医用微波热疗天线SAR图与肿瘤CT图像的匹配方法,其特征在于,微波阵列天线的阵列数量为1-5个,阵列排列方式有1、1*1、1*3、2*2或2*3。
5.根据权利要求4所述的一种医用微波热疗天线SAR图与肿瘤CT图像的匹配方法,其特征在于,微波阵列天线的相位为0°-360°。
6.根据权利要求1所述的一种医用微波热疗天线SAR图与肿瘤CT图像的匹配方法,其特征在于,图像数据处理包括下列步骤:1)图像转换:采用阈值将肿瘤CT图像和SAR图像转换为二值图像,
2)噪声消除:为了消除在边界或边界内部的“空洞”,先用填充函数将“空洞”进行填充,然后进行噪声消除,噪声消除采用均值滤波,用滤波函数滤波之后得到底图;
3)多边形近似:取滤波后底图,根据底图调用函数,把边界提取出来,呈现出一个多边形图片,最后通过这个多边形图片进行匹配,所述函数是根据最小周长多边形近似原理确定。
7.根据权利要求1所述的一种医用微波热疗天线SAR图与肿瘤CT图像的匹配方法,其特征在于,图像匹配:利用特征匹配中的面特征,对肿瘤图像和SAR图进行提取,在两者的面特征即轮廓提取出来后,计算面特征的二值特征,然后根据二值特征设定检索条件把符合条件的SAR图呈现出来,并显示出它所标注的微波阵列天线类型,所述二值特征包括:面积、偏心率、周长、重心、长宽比。
8.根据权利要求1所述的一种医用微波热疗天线SAR图与肿瘤CT图像的匹配方法,其特征在于,图像库内的每个SAR图上都标注出所对应的微波阵列天线具体参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510517718.2A CN105023276B (zh) | 2015-08-22 | 2015-08-22 | 一种医用微波热疗天线sar图与肿瘤ct图像的匹配方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510517718.2A CN105023276B (zh) | 2015-08-22 | 2015-08-22 | 一种医用微波热疗天线sar图与肿瘤ct图像的匹配方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105023276A true CN105023276A (zh) | 2015-11-04 |
CN105023276B CN105023276B (zh) | 2018-01-23 |
Family
ID=54413217
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510517718.2A Expired - Fee Related CN105023276B (zh) | 2015-08-22 | 2015-08-22 | 一种医用微波热疗天线sar图与肿瘤ct图像的匹配方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105023276B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112528468A (zh) * | 2020-11-20 | 2021-03-19 | 南京航空航天大学 | 一种基于电磁场sar图像处理的标签天线逆向设计方法 |
CN113902927A (zh) * | 2021-12-09 | 2022-01-07 | 北京车网科技发展有限公司 | 一种融合了图像与点云信息的综合信息处理方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1453052A (zh) * | 2002-04-24 | 2003-11-05 | 董宝玮 | 肿瘤微波介入治疗多天线热场组合方法 |
US20040151380A1 (en) * | 2003-01-30 | 2004-08-05 | Postech Foundation | Multi-layered real-time stereo matching method and system |
CN103714545A (zh) * | 2013-12-26 | 2014-04-09 | 北京工业大学 | 一种热消融治疗肿瘤的适形覆盖方法 |
WO2014201052A2 (en) * | 2013-06-10 | 2014-12-18 | University Of Mississippi Medical Center | Medical image processing method |
-
2015
- 2015-08-22 CN CN201510517718.2A patent/CN105023276B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1453052A (zh) * | 2002-04-24 | 2003-11-05 | 董宝玮 | 肿瘤微波介入治疗多天线热场组合方法 |
US20040151380A1 (en) * | 2003-01-30 | 2004-08-05 | Postech Foundation | Multi-layered real-time stereo matching method and system |
WO2014201052A2 (en) * | 2013-06-10 | 2014-12-18 | University Of Mississippi Medical Center | Medical image processing method |
CN103714545A (zh) * | 2013-12-26 | 2014-04-09 | 北京工业大学 | 一种热消融治疗肿瘤的适形覆盖方法 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112528468A (zh) * | 2020-11-20 | 2021-03-19 | 南京航空航天大学 | 一种基于电磁场sar图像处理的标签天线逆向设计方法 |
CN113902927A (zh) * | 2021-12-09 | 2022-01-07 | 北京车网科技发展有限公司 | 一种融合了图像与点云信息的综合信息处理方法 |
CN113902927B (zh) * | 2021-12-09 | 2022-04-12 | 北京车网科技发展有限公司 | 一种融合了图像与点云信息的综合信息处理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105023276B (zh) | 2018-01-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Nguyen et al. | Three-dimensional microwave hyperthermia for breast cancer treatment in a realistic environment using particle swarm optimization | |
Hameed et al. | A comprehensive survey on image-based computer aided diagnosis systems for skin cancer | |
Lak et al. | Evaluation of SAR distribution in six‐layer human head model | |
CN1380833A (zh) | 利用微波治疗乳房损害的设备和方法 | |
CN105023276A (zh) | 一种医用微波热疗天线sar图与肿瘤ct图像的匹配方法 | |
CN103598915B (zh) | 一种肿瘤射频消融技术中精确定位三球覆盖肿瘤的方法 | |
Karampatzakis et al. | Antenna design and tissue parameters considerations for an improved modelling of microwave ablation in the liver | |
Lodi et al. | Towards the robust and effective design of hyperthermic devices: Improvement of a patch antenna for the case study of abdominal rhabdomyosarcoma with 3D perfusion | |
Meenu et al. | A survey on heating effects of electromagnetic radiation on human body | |
CN104485521B (zh) | 一种适用于微波热疗的平面阵列聚焦天线 | |
Hassan et al. | Sar performance of rectangular microstrip antenna for breast cancer hyperthermia treatment with different period of treatment procedure | |
Polozov et al. | Cylindrical phased dipoles array for hyperthermia of deep-situated tumors | |
CN107715309B (zh) | 可穿戴微波热疗阵列天线 | |
Surita et al. | Finite element analysis for optimizing antenna for microwave coagulation therapy | |
Karli et al. | Early detection of breast tumors using UWB microstrip antenna imaging | |
Razib et al. | Microwave ablation technique (MWA) for cancer treatment: Simulation of single Slot MCA for different slot position | |
Aldhaeebi et al. | New antenna design for hyperthermia treatment of human head | |
Mawaddah et al. | Rhizome Image Classification Using Support Vector Machine | |
Baskaran et al. | Computer simulations of 434 MHz Electromagnetic Phased Array for thermal therapy of locally advanced breast cancer | |
Takook et al. | Evaluation of the 3D time reversal method for hyperthermia treatment planning in head and neck tumors | |
Chauhan et al. | Octuple Ridge Conical Horn Antenna with Hemispherical Dielectric Lens Loading for Biomedical Imaging Applications | |
Garcia-Miquel et al. | Compact microwave applicator for thermal therapy of breast cancer: Comparative assessment of arrays operating at 434 and 915 MHz | |
Wong et al. | Heat distribution improvement with the implementation of polyethylene-covered water bolus into breast cancer hyperthermia | |
CN110288669A (zh) | 一种压迫式超宽带微波乳腺肿瘤成像方法 | |
CN109464755A (zh) | 新型3d精准乳腺、胸壁适形器的制备方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180123 Termination date: 20180822 |
|
CF01 | Termination of patent right due to non-payment of annual fee |