CN106959306A - A kind of method that mammary tumor imaging is carried out using microwave reflection time domain S21 signals - Google Patents

A kind of method that mammary tumor imaging is carried out using microwave reflection time domain S21 signals Download PDF

Info

Publication number
CN106959306A
CN106959306A CN201710330039.3A CN201710330039A CN106959306A CN 106959306 A CN106959306 A CN 106959306A CN 201710330039 A CN201710330039 A CN 201710330039A CN 106959306 A CN106959306 A CN 106959306A
Authority
CN
China
Prior art keywords
time domain
signal
tumour
antenna
transmitting 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
Application number
CN201710330039.3A
Other languages
Chinese (zh)
Other versions
CN106959306B (en
Inventor
肖夏
王梁
李钦伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
JIANGSU MENSCH AUTO PARTS Co.,Ltd.
Original Assignee
Tianjin University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tianjin University filed Critical Tianjin University
Priority to CN201710330039.3A priority Critical patent/CN106959306B/en
Publication of CN106959306A publication Critical patent/CN106959306A/en
Application granted granted Critical
Publication of CN106959306B publication Critical patent/CN106959306B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N22/00—Investigating or analysing materials by the use of microwaves or radio waves, i.e. electromagnetic waves with a wavelength of one millimetre or more

Abstract

The present invention relates to a kind of method that mammary tumor imaging is carried out using microwave reflection time domain S21 signals, including：Prepare cavity breast model and suspention tumour prosthese, make the breast model with tumour；By the port 1 of transmitting antenna and feeder line connected vector Network Analyzer, by the port 2 of reception antenna and feeder line connected vector Network Analyzer, Network Analyzer time domain S is excited21Measurement pattern；Transmitting antenna and reception antenna are close to placement location, the transmission path time domain S of record now when measuring for the first time21Signal；It is during second of measurement, transmitting antenna and reception antenna is separated by a distance, record new transmission path time domain S21Signal；Time domain S when will have a tumour prosthese21Time domain S when signal is with without tumour prosthese21Signal subtraction, can obtain the time domain S as caused by tumour prosthese21Signal intensity；Using confocal method to subtracting each other after time domain S21Signal carries out tumor imaging.

Description

A kind of method that mammary tumor imaging is carried out using microwave reflection time domain S21 signals
Technical field
The invention belongs to microwave Imaging Technique field, it is related to a kind of mammary tumor imaging method.
Background technology
Breast cancer turns into modern society women fatal disease occurred frequently, and early making a definite diagnosis early treatment to improving disease cured rate has weight Want meaning.It is high, expensive etc. to there is radiation risk, misdiagnosis rate in existing detection method such as molybdenum target X-ray, ultrasound, nuclear-magnetism etc. Shortcoming.Microwave is detected as a kind of emerging lossless detection method, can effectively detect intramammary infantile tumour, and it is clear into Picture.Tumor imaging is based primarily upon Near-field radar imaging principle in microwave breast, and the emitted antenna of microwave signal is incident in breast, Because tumour with normal structure has dielectric property difference, the reflected signal that incoming signal meets with after tumour can carry different information And finally received by reception antenna.This different information after signal processing, can show knub position.Conventional microwave breast into It is imaged as technology using time-domain signal, this method imaging complexity is, it is necessary to be provided simultaneously with microwave signal transmitter module, microwave Signal receiving module, costly, time synchronized is complicated, operating difficulties for price.Therefore, it is necessary to develop a kind of system composition letter Single, imaging method with low cost.
The content of the invention
It is an object of the invention to provide a kind of method of simple and efficient mammary tumor imaging, the present invention uses microwave reflection Time domain S21Signal carries out tumor imaging in breast, and the method for the transmitting time domain fixed waveform signal of abandoning tradition utilizes dual-port Signal meets with the total propagation time returned again to after tumour in the incident breast of time domain model reacting condition in transmission path, and uses Signal processing method had, without subtracting each other time domain S under tumor presence21Signal, most afterwards through confocal algorithm to tumour focal imaging. Technical scheme is as follows：
A kind of method that mammary tumor imaging is carried out using microwave reflection time domain S21 signals, is comprised the following steps：
1) prepare cavity breast model and suspention tumour prosthese, make the breast model with tumour；
2) by the port 1 of transmitting antenna and feeder line connected vector Network Analyzer, by reception antenna and feeder line connected vector The port 2 of Network Analyzer, excites Network Analyzer time domain S21Measurement pattern；
3) transmitting antenna and reception antenna are close to placement location, the transmission path time domain of record now when measuring for the first time S21Signal
4) it is during second of measurement, transmitting antenna and reception antenna is separated by a distance, still keep to configuration state, record New transmission path time domain S21Signal, new time domain S21The more former S of signal21The time location that acute variation occurs in signal is closure The position of air impedance is initially added into transmission path, will appear from change time location be set to 0 initial time, this 0 when Blaze will transmits signals to reception antenna when transmitting antenna and reception antenna are close to placement location, from transmitting antenna and receives signal Required time；
5) transmitting antenna and reception antenna are arranged in the breast model surface with tumour and carry out imaging experiment, two days Port 1 of the line through coaxial feeder and vector network analyzer, port 2 are respectively connected with；
6) 3.5~15GHz scanning bands of excitation vector Network Analyzer, record time domain S21Signal data；
7) tumour prosthese is taken out from breast model, time domain S is recorded again21Signal data；
8) time domain S when will have a tumour prosthese21Time domain S when signal is with without tumour prosthese21Signal subtraction, can obtain by Time domain S caused by tumour prosthese21Signal intensity；
9) using confocal method to subtracting each other after time domain S21Signal carries out tumor imaging.
Brief description of the drawings
Fig. 1 imaging data instrumentation plans
Fig. 2 (a) measures closed path time domain S21 curves twice；(b) aerial position figure is measured for the first time；(c) second of survey Measure aerial position figure
The time domain S21 signal differences that Fig. 3 is measured twice
Fig. 4 subtract each other after time domain S21 signals
The aerial position of Fig. 5 records
Fig. 6 tumour confocal imaging result figures
Embodiment
Measure the moment of transmit antenna port 0 in the dual-port transmission path that antenna, feeder line are constituted.Method as shown in Fig. 2 By the port 1 of transmitting antenna and feeder line connected vector Network Analyzer, by reception antenna and feeder line connected vector Network Analyzer Port 2, excite Network Analyzer time domain S21Measurement pattern.Transmitting antenna and reception antenna are close to when measuring for the first time opposed Place, the transmission path time domain S of record now21Signal.Now, whole antenna, coaxial feeder system default can be closed for one Combining footpath, time domain S21Signal embodies the impedance value in transmission path.Second when measuring, will transmitting, reception antenna it is separately big About 8cm distance, still keeps to configuration state, records new transmission path time domain S21Signal.Now, new two-port network is suitable New closure transmission path is generated in adding one section of air in the antenna of first time measurement, coaxial feeder system-based, newly Time domain S21The more former S of signal21The time location that acute variation occurs in signal is to be initially added into air impedance in closure transmission path Position.The time location that will appear from change is set to 0 initial time, and this 0 moment mark works as transmitting antenna and reception antenna When being close to placement location, transmit signals to reception antenna from transmitting antenna and receive signal required time.
1st, it is by progress imaging experiment as shown in Figure 1, PE breast models frame is good, the fatty edible oil of class is filled into cavity, By the suspention of tumour prosthese wherein.Transmitting antenna and reception antenna are respectively placed in model surface, two antennas through coaxial feeder with The port 1 of vector network analyzer, port 2 are respectively connected with.3.5~15GHz scanning bands of excitation vector Network Analyzer, note Record time domain S21Signal data.Tumour prosthese is taken out from breast model, time domain S is recorded again21Signal data.Time domain twice S21Signal difference is as shown in figure 3, time domain S when will have a tumour prosthese21Time domain S when signal is with without tumour prosthese21Signal phase Subtract, the time domain S as caused by tumour prosthese can be obtained21Signal intensity.As shown in figure 4, peak value corresponding generation time of this change Table microwave signal meets with tumour prosthese and produces reflection, the closed path for returning again to reception antenna is always passed after transmitting antenna is sent Between sowing time.Time domain S after confocal method can be used directly to subtracting each other21Signal carries out tumor imaging.
2nd, after signal is sent through a certain transmitting antenna, position of transmitting antenna is recorded, while receiving the antenna of reflected signal Position is also recorded, and two positions are substituted into the formula (1) of confocal imaging algorithm, make signal confocal.
3rd, transmitting antenna, the position of reception antenna, duplicate measurements, the time domain S after obtaining multigroup subtract each other are changed21Letter Number.Such as Fig. 5, breast surface round dot show transmitting antenna and the reception antenna position of record.Using formula (1)~the (3) Suo Shi Confocal imaging algorithm by the confocal processing of signal：
Wherein I represents the image pixel intensities at the confocal point r of imaging region, rmThe position for representing transmitting antenna is m points, rnRepresent The position of reception antenna is n points, and τ represents tumor signal from rmTransmitting antenna is sent at point reaches r again through confocal point rnReceived at point Total propagation time of antenna, v representation signals the mean propagation velocity in breast model, c represents the light velocity, εfatRepresent that lipoid is eaten The dielectric constant of oil.
Imaging results are as shown in Figure 6, it is seen that use time domain S21Signal can carry out blur-free imaging to tumour.Time domain S21Signal It can reflect and impedance variations when are there occurs in dual-port transmission path, it is simple to operate, only need dual-port vector network analyzer i.e. It can complete, image space is accurate, and focusing effect is good.Compared to the combination waving map pattern using signal generator and oscillograph, The quantity of image-forming module is reduced, is a kind of inexpensive high speed imaging method.

Claims (1)

1. a kind of method that mammary tumor imaging is carried out using microwave reflection time domain S21 signals, is comprised the following steps：
1) prepare cavity breast model and suspention tumour prosthese, make the breast model with tumour；
2) by the port 1 of transmitting antenna and feeder line connected vector Network Analyzer, by reception antenna and feeder line connected vector network The port 2 of analyzer, excites Network Analyzer time domain S21Measurement pattern；
3) transmitting antenna and reception antenna are close to placement location, the transmission path time domain S of record now when measuring for the first time21Letter Number
4) it is during second of measurement, transmitting antenna and reception antenna is separated by a distance, still keep, to configuration state, recording newly Transmission path time domain S21Signal, new time domain S21The more former S of signal21The time location that acute variation occurs in signal is closure transmission The position of air impedance is initially added on path, will appear from change time location be set to 0 initial time, this 0 when blaze Will is transmitted signals to needed for reception antenna receives signal when transmitting antenna and reception antenna are close to placement location from transmitting antenna Time；
5) transmitting antenna and reception antenna are arranged in the breast model surface with tumour and carry out imaging experiment, two antenna warps The port 1 of coaxial feeder and vector network analyzer, port 2 are respectively connected with；
6) 3.5~15GHz scanning bands of excitation vector Network Analyzer, record time domain S21Signal data；
7) tumour prosthese is taken out from breast model, time domain S is recorded again21Signal data；
8) time domain S when will have a tumour prosthese21Time domain S when signal is with without tumour prosthese21Signal subtraction, can obtain false by tumour Time domain S caused by body21Signal intensity；
9) using confocal method to subtracting each other after time domain S21Signal carries out tumor imaging.
CN201710330039.3A 2017-05-11 2017-05-11 A method of mammary tumor imaging is carried out using microwave reflection time domain S21 signal Expired - Fee Related CN106959306B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710330039.3A CN106959306B (en) 2017-05-11 2017-05-11 A method of mammary tumor imaging is carried out using microwave reflection time domain S21 signal

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710330039.3A CN106959306B (en) 2017-05-11 2017-05-11 A method of mammary tumor imaging is carried out using microwave reflection time domain S21 signal

Publications (2)

Publication Number Publication Date
CN106959306A true CN106959306A (en) 2017-07-18
CN106959306B CN106959306B (en) 2019-11-12

Family

ID=59482348

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710330039.3A Expired - Fee Related CN106959306B (en) 2017-05-11 2017-05-11 A method of mammary tumor imaging is carried out using microwave reflection time domain S21 signal

Country Status (1)

Country Link
CN (1) CN106959306B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109350078A (en) * 2018-10-19 2019-02-19 天津大学 Blood sugar concentration detection method based on ultra-wideband microwave inverse Fourier transform
CN109350053A (en) * 2018-10-19 2019-02-19 深圳市太赫兹科技有限公司 A kind of brain imaging method and its system, equipment, storage medium
CN109350054A (en) * 2018-10-19 2019-02-19 深圳市太赫兹科技有限公司 A kind of brain imaging system
CN111081380A (en) * 2019-12-26 2020-04-28 天津大学 Method for optimizing microwave breast image based on image quality index and simulated annealing
CN112569006A (en) * 2020-12-11 2021-03-30 倪超 Microwave-based non-contact intraoperative real-time skin flap thickness monitoring system and measuring and calculating method thereof

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101234022A (en) * 2006-12-19 2008-08-06 华东师范大学 Microwave near-field medicine body detecting method and use thereof
CN103300826A (en) * 2013-06-06 2013-09-18 天津大学 Imaging method for early breast tumor ultra wide band microwave detection
CN105342556A (en) * 2015-09-23 2016-02-24 天津大学 Tumor signal extraction method based on integrated empirical mode decomposition
CN105701814A (en) * 2016-01-13 2016-06-22 天津大学 Ultra-wideband breast tumor imaging method based on magnetic resonance image compensation
CN105816172A (en) * 2016-03-11 2016-08-03 金陵科技学院 Brain tumor microwave detection system

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101234022A (en) * 2006-12-19 2008-08-06 华东师范大学 Microwave near-field medicine body detecting method and use thereof
CN103300826A (en) * 2013-06-06 2013-09-18 天津大学 Imaging method for early breast tumor ultra wide band microwave detection
CN105342556A (en) * 2015-09-23 2016-02-24 天津大学 Tumor signal extraction method based on integrated empirical mode decomposition
CN105701814A (en) * 2016-01-13 2016-06-22 天津大学 Ultra-wideband breast tumor imaging method based on magnetic resonance image compensation
CN105816172A (en) * 2016-03-11 2016-08-03 金陵科技学院 Brain tumor microwave detection system

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
LIANG WANG ET AL: "A novel dual calibration technique for microwave early breast tumour detection in MR images-based model", 《JOURNAL OF MICROWAVE POWER AND ELECTROMAGNETIC ENERGY》 *
LIANG WANG ET AL: "MRI-aided tissues interface characterization: An accurate signal propagation time calculation method for UWB breast tumor imaging", 《APPLIED SURFACE SCIENCE》 *
M. A. ELAHI ET AL: "Hybrid Artifact Removal for Confocal Microwave Breast Imaging", 《ANTENNAS AND WIRELESS PROPAGATION LETTERS》 *
WANG LIANG ET AL: "Contour Detection-Based Realistic Finite-Difference-Time-Domain Models for Microwave Breast Cancer Detection", 《TRANSACTIONS OF TIANJIN UNIVERSITY》 *
XIAO XIA ET AL: "A progressive processing method for breast cancer detection via UWB based on an MRI-derived model", 《CHIN. PHYS. B》 *
肖夏等: "早期乳腺肿瘤的超宽带微波稳健波束形成成像检测系统", 《物理学报》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109350078A (en) * 2018-10-19 2019-02-19 天津大学 Blood sugar concentration detection method based on ultra-wideband microwave inverse Fourier transform
CN109350053A (en) * 2018-10-19 2019-02-19 深圳市太赫兹科技有限公司 A kind of brain imaging method and its system, equipment, storage medium
CN109350054A (en) * 2018-10-19 2019-02-19 深圳市太赫兹科技有限公司 A kind of brain imaging system
CN111081380A (en) * 2019-12-26 2020-04-28 天津大学 Method for optimizing microwave breast image based on image quality index and simulated annealing
CN111081380B (en) * 2019-12-26 2023-04-07 天津大学 Method for optimizing microwave breast image based on image quality index and simulated annealing
CN112569006A (en) * 2020-12-11 2021-03-30 倪超 Microwave-based non-contact intraoperative real-time skin flap thickness monitoring system and measuring and calculating method thereof

Also Published As

Publication number Publication date
CN106959306B (en) 2019-11-12

Similar Documents

Publication Publication Date Title
CN106959306B (en) A method of mammary tumor imaging is carried out using microwave reflection time domain S21 signal
US20230404422A1 (en) Sensors for a portable device
CN105120756A (en) Microwave imaging resilient to background and skin clutter
EP2234539B1 (en) Synthetic aperture radar system
CN106707291A (en) Laser radar system
US20210255314A1 (en) Synthetic ultra-wideband millimeter-wave imaging for tissue diagnostics
US20100315280A1 (en) Multi-Platform Radar with Forced Resonating Antennas for Embedded Detection and Volumetric Imaging
CN104169739A (en) Spread spectrum coded waveforms in ultrasound imaging
CN103176183B (en) Short-range high-resolution detection imaging method and device for implementing same
Zamani et al. Boundary estimation of imaged object in microwave medical imaging using antenna resonant frequency shift
JP2005193017A (en) Method and system for classifying diseased part of mamma
CN105142528B (en) Diagnostic ultrasound equipment
CN106405661A (en) Hand-held safety-check device and safety-check method
Khalesi et al. Free-space operating microwave imaging device for bone lesion detection: A phantom investigation
Sanabria et al. Hand-held sound-speed imaging based on ultrasound reflector delineation
Noritake et al. Boundary extraction enhanced distorted born iterative method for microwave mammography
WO2018083492A1 (en) A breast density meter and method
CN106022369A (en) Bone quality evaluation method based on ultrasonic backscatter signal imaging
Gezimati et al. Curved synthetic aperture radar for near-field terahertz imaging
Jin et al. Time reversal beamforming for microwave breast cancer detection
Sarafianou et al. Breast surface reconstruction algorithm for a multi-static radar-based breast imaging system
Solimene et al. An incoherent radar imaging system for medical applications
Jin et al. Breast cancer detection by time reversal imaging
JP6651012B2 (en) Method, apparatus, and computer readable medium for assessing fit in a system for probing the internal structure of an object
Danjuma et al. Microwave imaging using arrays of Vivaldi antenna for breast cancer applications

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
TR01 Transfer of patent right

Effective date of registration: 20201225

Address after: 224003 No.66, Xinyang Road, Tinghu District, Yancheng City, Jiangsu Province (18)

Patentee after: JIANGSU MENSCH AUTO PARTS Co.,Ltd.

Address before: 300072 Tianjin City, Nankai District Wei Jin Road No. 92

Patentee before: Tianjin University

TR01 Transfer of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20191112

Termination date: 20200511

CF01 Termination of patent right due to non-payment of annual fee