CN107212884B - 一种仰卧***压迫***成像方法 - Google Patents
一种仰卧***压迫***成像方法 Download PDFInfo
- Publication number
- CN107212884B CN107212884B CN201710329825.1A CN201710329825A CN107212884B CN 107212884 B CN107212884 B CN 107212884B CN 201710329825 A CN201710329825 A CN 201710329825A CN 107212884 B CN107212884 B CN 107212884B
- Authority
- CN
- China
- Prior art keywords
- breast
- compression
- tumor
- point
- pixels
- 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.)
- Expired - Fee Related
Links
- 210000000481 breast Anatomy 0.000 title claims abstract description 58
- 238000007906 compression Methods 0.000 title claims abstract description 39
- 230000006835 compression Effects 0.000 title claims abstract description 33
- 238000003384 imaging method Methods 0.000 title claims abstract description 24
- 206010028980 Neoplasm Diseases 0.000 claims abstract description 41
- 238000001514 detection method Methods 0.000 claims abstract description 14
- 238000004088 simulation Methods 0.000 claims abstract description 14
- 230000004044 response Effects 0.000 claims abstract description 7
- 238000010226 confocal imaging Methods 0.000 claims abstract description 5
- 238000013507 mapping Methods 0.000 claims abstract description 4
- 238000000034 method Methods 0.000 claims description 4
- 210000000115 thoracic cavity Anatomy 0.000 claims description 3
- 238000006073 displacement reaction Methods 0.000 claims description 2
- 206010006187 Breast cancer Diseases 0.000 description 8
- 208000026310 Breast neoplasm Diseases 0.000 description 8
- 238000011897 real-time detection Methods 0.000 description 3
- 238000001356 surgical procedure Methods 0.000 description 2
- 238000005481 NMR spectroscopy Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 210000002445 nipple Anatomy 0.000 description 1
- 230000002980 postoperative effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
Images
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
- 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/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/70—Means for positioning the patient in relation to the detecting, measuring or recording means
- A61B5/708—Breast positioning means
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Veterinary Medicine (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- High Energy & Nuclear Physics (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明涉及一种仰卧***压迫***成像方法,方法如下:选用***MRI影像图组进行***受压迫建模,对于***边界低于压缩面的列像素进行向上压缩处理,对每列像素进行压缩处理后,对MRI影像图进行适当的横向展宽后完成整个压缩过程;对MRI影像图组的每一张图都进行上述压缩处理,每张图的皮肤压缩终点都选取同一位置C点,将所有压缩后的图按顺序叠放、通过像素点灰度映射导出***受压迫的三维模型;对于***进行仿真成像时,先对不含肿瘤的***模型进行一次仿真检测,再对含肿瘤的***模型进行一次仿真检测。两次检测信号相减获取纯肿瘤响应信号;利用共焦成像算法将纯肿瘤响应信号共焦处理。
Description
技术领域
本发明属于微波成像技术领域,涉及一种微波***内肿瘤成像方法。
背景技术
乳腺癌对女性生命健康威胁很大,手术治疗是国际上根治乳腺癌的通用方法。早发现早手术被认为是提高乳腺癌治愈率、防止乳腺癌复发的有效手段。但是,面对小尺寸肿瘤,医生手术时往往无法确定肿瘤的准确位置,只能扩大手术创口,在一定范围内寻找肿瘤。过大的创口极大地削弱了患者术后恢复力,患者抵抗肿瘤复发的能力大大减弱,导致手术无法达到预期目的,降低了手术效费比。因此,***内肿瘤手术实时检测变得非常重要,现有的***检测手段如核磁共振、X射线、超声等因其设备体积、射线伤害、图像不直观等因素无法作为有效的手术实时检测设备。超宽带微波***内肿瘤成像技术因其设备体积小,价格低廉,长时间使用对人体无伤害等诸多优点有望成为手术实时检测设备。微波***内肿瘤成像主要基于雷达近场成像原理,微波信号经发射天线入射至***内,由于肿瘤同正常组织存在介电特性差异,入射信号遭遇肿瘤后的反射信号可携带差异信息并最终被接收天线接收。这种差异信息经信号处理后,可以显示肿瘤位置。
发明内容
本发明提供一种处于手术仰卧***的***内肿瘤成像方法,技术方案如下:
一种仰卧***压迫***成像方法,人体处于仰卧***时,使用平板阵列天线压置于***表面,在***皮肤近似平面的情况下,寻找位于阵列天线正下方的肿瘤准确位置,方法如下:
1)选用***MRI影像图组进行***受压迫建模,对于***边界低于压缩面的列像素进行向上压缩处理,设某列像素起始点为空气中的A点,终止点为胸腔内的C点,压缩终点置为B点,对该列像素AC段上的所有像素点进行移位处理,使AC段上像素点压缩至BC段内,对每列像素进行压缩处理后,对MRI影像图进行适当的横向展宽后完成整个压缩过程;
2)对MRI影像图组的每一张图都进行上述压缩处理,每张图的皮肤压缩终点都选取同一位置C点,将所有压缩后的图按顺序叠放、通过像素点灰度映射导出***受压迫的三维模型;
3)对于***进行仿真成像时,先对不含肿瘤的***模型进行一次仿真检测,再对含肿瘤的***模型进行一次仿真检测。两次检测信号相减获取纯肿瘤响应信号
4)利用共焦成像算法将纯肿瘤响应信号共焦处理。
附图说明
图1***压迫建模的MRI图像处理示意图
图2受压迫的***MRI图像
图3受压迫的***模型示意图
图4患者仰卧***受压迫***模型切片图
图5患者仰卧******模型内肿瘤成像图
具体实施方式
仰卧时,***因其自身的多脂肪结构在重力的作用下会趋向于平坦。此时,肿瘤距离***皮肤距离减小,距离减小将导致阵列天线聚焦性能大幅度下降,阵列天线适当压迫***即能保障天线聚焦成像的性能达标,又能保证***皮肤被有效压成平面,便于阵列天线同***皮肤紧密接触。具体方案如下:
仿真实验使用***MRI影像图组进行***受压迫建模。如图1所示对于***边界低于压缩面的列像素进行向上压缩处理。设某列像素起始点为A点(空气),终止点为C点(胸腔),压缩终点置为B点。对该列像素AC段上的所有像素点进行移位处理,使AC段上像素点压缩至BC段内。对每列像素进行压缩处理后,对MRI影像图进行适当的横向展宽后完成整个压缩过程,如图2所示。
对MRI影像图组的每一张图都进行上述压缩处理。每张图的皮肤压缩终点都选取同一位置C点。将所有压缩后的图按顺序叠放、通过像素点灰度映射导出***受压迫的三维模型,如图3所示。
图3中皮肤表面圆点表示超宽带阵列天线的位置,本仿真成像共使用32个超宽带天线,天线环形排布在皮肤平面上,中间空出***的位置。每个天线即可作为发射天线,也可作为接收天线使用,信号经某一发射天线发出后,发射天线位置被记录,同时接收反射信号的天线位置也被记录,两个位置代入共焦成像算法的公式(1)中,使信号共焦。***模型经过轻度压缩,压缩比约为26%。经压缩后,患者***模型在横向上有所增大,***皮肤被压缩为一个平面。半径为2mm的球形肿瘤模型被放置在***内x=22mm,y=73mm,z=50mm位置处,肿瘤位置如图4所示。进行成像仿真时,首先使用有肿瘤***模型进行超宽带检测,获得接收信号。然后,使用不带有肿瘤的***模型再进行一次检测获取参考信号。用有肿瘤情况下的接收信号减去无肿瘤情况下的参考信号获取纯肿瘤响应。
使用公式(1)~(3)所示的共焦成像算法将信号共焦处理:
其中I代表成像区域共焦点r处的像素强度,rm代表发射天线的位置是m点,rn代表接收天线的位置是n点,τ代表肿瘤信号从rm点处发射天线发出经共焦点r再到达rn点处接收天线的总传播时间,v代表信号在***模型内平均传播速度,c代表光速,εfat表示***内脂肪的介电常数。
图5给出了仰卧******模型内肿瘤成像仿真结果图。成像结果显示肿瘤位置为x=16mm,y=72mm,z=49mm,同预设的肿瘤位置相符合。可以看出,在受阵列天线轻度压迫下,受压迫的***仍可接受超宽带微波检测并正确成像,成像位置较准确。在俯视***的情况下,在俯视平面的两个方向肿瘤成像偏差(y向和z向,而x向表示肿瘤深度)仅为1mm,非常准确地定位了乳腺肿瘤在体内位置,而手术下刀一般向下碰触肿瘤即可,所以肿瘤深度(x向)偏差为6mm对手术影响非常小。
综上,在***受压迫的情况下,即使肿瘤距离阵列天线较近的情况下,肿瘤成像位置在俯视平面十分准确,对于增强手术下刀准确性,减小创口具有重要意义。
Claims (1)
1.一种仰卧***压迫***成像方法,用于仿真实现,模拟人体处于仰卧***时,使用平板阵列天线压置于***表面,在***皮肤近似平面的情况下,寻找位于阵列天线正下方的肿瘤准确位置,方法如下:
1)选用***MRI影像图组进行***受压迫建模,对于***边界低于压缩面的列像素进行向上压缩处理,设某列像素起始点为空气中的A点,终止点为胸腔内的C点,压缩终点置为B点,对该列像素AC段上的所有像素点进行移位处理,使AC段上像素点压缩至BC段内,对每列像素进行压缩处理后,对MRI影像图进行适当的横向展宽后完成整个压缩过程;
2)对MRI影像图组的每一张图都进行上述压缩处理,每张图的皮肤压缩终点都选取同一位置C点,将所有压缩后的图按顺序叠放、通过像素点灰度映射导出***受压迫的三维模型;
3)对***进行仿真成像时,先对不含肿瘤的***模型进行一次仿真检测,再对含肿瘤的***模型进行一次仿真检测,两次检测信号相减获取纯肿瘤响应信号;
4)利用共焦成像算法将纯肿瘤响应信号共焦处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710329825.1A CN107212884B (zh) | 2017-05-11 | 2017-05-11 | 一种仰卧***压迫***成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710329825.1A CN107212884B (zh) | 2017-05-11 | 2017-05-11 | 一种仰卧***压迫***成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107212884A CN107212884A (zh) | 2017-09-29 |
CN107212884B true CN107212884B (zh) | 2020-07-24 |
Family
ID=59944839
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710329825.1A Expired - Fee Related CN107212884B (zh) | 2017-05-11 | 2017-05-11 | 一种仰卧***压迫***成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107212884B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110288669A (zh) * | 2019-06-11 | 2019-09-27 | 天津大学 | 一种压迫式超宽带微波乳腺肿瘤成像方法 |
CN110664405B (zh) * | 2019-09-27 | 2021-12-03 | 天津大学 | 基于焦点质量度量估计微波***成像平均介电特性的方法 |
CN111067524B (zh) * | 2019-12-26 | 2021-12-03 | 天津大学 | 一种估计微波***成像的平均介电特性的方法 |
WO2023070387A1 (zh) * | 2021-10-27 | 2023-05-04 | 深圳市大疆创新科技有限公司 | 一种图像处理方法、装置、拍摄设备及可移动平台 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103462589A (zh) * | 2013-09-10 | 2013-12-25 | 天津大学 | 一种检测乳腺多发肿瘤的超宽带微波成像方法 |
CN103549953A (zh) * | 2013-10-25 | 2014-02-05 | 天津大学 | 一种基于医学核磁共振影像提取微波检测***模型的方法 |
CN103799982A (zh) * | 2014-01-26 | 2014-05-21 | 天津大学 | 基于希尔伯特黄变换的超宽带微波检测的快速成像方法 |
CN104794749A (zh) * | 2015-04-08 | 2015-07-22 | 天津大学 | 基于临床mri影像建立三维***电磁仿真模型的方法 |
CN105899144A (zh) * | 2014-01-16 | 2016-08-24 | 佳能株式会社 | 图像处理装置、图像诊断***、图像处理方法和存储介质 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006028396A1 (en) * | 2004-09-10 | 2006-03-16 | Industrial Research Limited | Imaging system |
WO2014205023A1 (en) * | 2013-06-18 | 2014-12-24 | Northeastern University | Microwave nearfield radar imaging (nri) using digital breast tomosynthesis (dbt) for non-invasive breast cancer detection |
-
2017
- 2017-05-11 CN CN201710329825.1A patent/CN107212884B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103462589A (zh) * | 2013-09-10 | 2013-12-25 | 天津大学 | 一种检测乳腺多发肿瘤的超宽带微波成像方法 |
CN103549953A (zh) * | 2013-10-25 | 2014-02-05 | 天津大学 | 一种基于医学核磁共振影像提取微波检测***模型的方法 |
CN105899144A (zh) * | 2014-01-16 | 2016-08-24 | 佳能株式会社 | 图像处理装置、图像诊断***、图像处理方法和存储介质 |
CN103799982A (zh) * | 2014-01-26 | 2014-05-21 | 天津大学 | 基于希尔伯特黄变换的超宽带微波检测的快速成像方法 |
CN104794749A (zh) * | 2015-04-08 | 2015-07-22 | 天津大学 | 基于临床mri影像建立三维***电磁仿真模型的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107212884A (zh) | 2017-09-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107212884B (zh) | 一种仰卧***压迫***成像方法 | |
JP6309376B2 (ja) | 医用情報処理システム、医用情報処理プログラム及び超音波診断装置 | |
US10772606B2 (en) | Method and apparatus for displaying ultrasound images | |
JP4888026B2 (ja) | 超音波映像システム及び超音波映像処理方法 | |
CN103445866A (zh) | 外科规划***和导航*** | |
CN104408398B (zh) | 一种肝脏边界的识别方法及*** | |
CN109009110A (zh) | 基于mri影像的腋窝***转移预测*** | |
CN103417301A (zh) | 外科手术导航*** | |
US10368841B2 (en) | Ultrasound diagnostic apparatus | |
EP3193280B1 (en) | Medical imaging device and method of operating the same | |
CN109512464A (zh) | 一种疾病筛查和诊断*** | |
CN109636806A (zh) | 一种基于多步学习的三维核磁共振胰腺图像分割方法 | |
CN107714070B (zh) | 基于数字化断层融合图的乳腺病灶定位方法、***和装置 | |
CN109199381A (zh) | 一种全息微波弹性成像***及其成像方法 | |
KR102433473B1 (ko) | 환자의 증강 현실 기반의 의료 정보를 제공하는 방법, 장치 및 컴퓨터 프로그램 | |
CN104794749A (zh) | 基于临床mri影像建立三维***电磁仿真模型的方法 | |
CN104822326A (zh) | 用于消融治疗的超声体积流量测量 | |
CN114159099A (zh) | 乳腺超声成像方法及设备 | |
Frighetto-Pereira et al. | Semiautomatic classification of benign versus malignant vertebral compression fractures using texture and gray-level features in magnetic resonance images | |
US9449425B2 (en) | Apparatus and method for generating medical image | |
Zamani et al. | Frequency domain method for early stage detection of congestive heart failure | |
Sammali et al. | Experimental setup for objective evaluation of uterine motion analysis by ultrasound speckle tracking | |
Hussain et al. | Strain-initialized robust bone surface detection in 3-D ultrasound | |
CN112641471A (zh) | 一种放疗专用的膀胱容量测定与三维形态评估方法与*** | |
CN105701814B (zh) | 一种基于磁共振图像补偿的超宽带乳腺肿瘤成像方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200724 |