CN105701814B - A kind of ultra wide band tumor of breast imaging method based on magnetic resonance image compensation - Google Patents
A kind of ultra wide band tumor of breast imaging method based on magnetic resonance image compensation Download PDFInfo
- Publication number
- CN105701814B CN105701814B CN201610019122.4A CN201610019122A CN105701814B CN 105701814 B CN105701814 B CN 105701814B CN 201610019122 A CN201610019122 A CN 201610019122A CN 105701814 B CN105701814 B CN 105701814B
- Authority
- CN
- China
- Prior art keywords
- gland
- breast
- elliptic cone
- wide band
- ultra wide
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- 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/10088—Magnetic resonance imaging [MRI]
-
- 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/30068—Mammography; Breast
Landscapes
- Engineering & Computer Science (AREA)
- Quality & Reliability (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Health & Medical Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
The present invention relates to a kind of ultra wide band tumor of breast imaging methods based on magnetic resonance image compensation, including：The 3 d-dem breast model of MRI slice atlas structure bodies of gland based on breast；Gland structure is modeled, is elliptic cone shape by the structure proximate of gland tissue in breast；Use space coordinate transform, will be in elliptic cone " righting " to a secondary coordinate system；An isosceles triangle region is created in the plane projection region of elliptic cone, and by this region finite element, determines the elliptic cone of the variation tendency of description body of gland rough surface that can be best；Microwave propagation time estimation error is carried out, accurately image is carried out.The present invention can highly accurately carry out ultra wide band tumor of breast imaging.
Description
Technical field
The invention belongs to a kind of tumor of breast imaging methods.
Background technology
Ultra wide band (UWB) detection technique is quickly grown in field of non destructive testing.Very extensive, the technology is applied in UWB detections
It has a extensive future in numerous areas such as Ground Penetrating Radar, material flaw detection, through-wall radar, tumor of breast detections, there is detection device
At low cost, many advantages, such as equipment size is small, injury radiationless to human body.The technology uses the incident object to be detected of ultra-wideband microwave
Body, then by radar reflection principle positioning target or target is imaged.It is swollen to human breast in ultra-wideband microwave
Tumor is detected in application, the labyrinth of breast cause it is ultrafast with microwave propagation when speed it is inhomogenous and be difficult to estimate, cause
Imaging accuracy declines.An important method for promoting imaging accuracy is to carry out breast tissue composition measurement.By estimating breast
The volume of gland tissue and adipose tissue in room carrys out the correct radar return arrival time estimated needed for imaging, and it is accurate to carry out
Imaging.
Magnetic resonance (MRI) image can accurately reflect human breast inner tissue complex distributions situation in gray scale.Due to breast
Interior gland tissue and adipose tissue are woven into a mesh structure, and the two interface is very coarse, and gland tissue volume is difficult to accurately be estimated
Meter.
Invention content
The object of the present invention is to provide a kind of high accuracy ultra wide band tumor of breast imaging methods.Technical scheme of the present invention
It is as follows：
A kind of ultra wide band tumor of breast imaging method based on magnetic resonance image compensation, includes the following steps：
(1) the 3 d-dem breast model of the MRI slice atlas structure bodies of gland based on breast, by MRI slice maps
Gray value interval carries out threshold definitions, internal mammary is divided into a variety of different tissues, by the body of gland part in discrete model
It extracts；
(2) gland structure is modeled, keeps the model of structure approximate with former gland structure, by gland tissue in breast
Structure proximate is elliptic cone shape, chooses the vertex of oval wimble structure first, then chooses the bottom center of oval wimble structure again
Position obtains the central symmetry axis of elliptic cone；
(3) use space coordinate transform keeps the same coordinate in its bottom surface flat by elliptic cone " righting " to a secondary coordinate system
Face is parallel；
(4) an isosceles triangle region is created in the plane projection region of elliptic cone, and by this region finite element, often
A grid is denoted as pi, wherein the half of the bottom edge length of isosceles triangle is the long axis a of elliptic cone, the height h and length of elliptic cone
Axis a is remained unchanged, and adjusts the size of short axle b, elliptic cone is made to be risen or fallen in the both direction up and down of gore, when
When short axle b takes a certain length, p is calculated firstiThe upper circular conical surface of point then is calculated with the absolute value of the difference of highest body of gland pixel
piThe two absolute values are included in a sample, wait for institute by the lower circular conical surface of point with the absolute value of the difference of minimum body of gland pixel
There is piAfter point completes same individual of sample evaluation work, sample variance is calculated, it is most to take b values when sample variance minimum
Suitable minor axis dimension, the variation tendency of the best description body of gland rough surface of elliptic cone energy at this time；
(5) microwave propagation time estimation error is carried out, accurately image is carried out, method is as follows：Give tacit consent to ultra-wideband antenna arrangement
Position is preceding half of breast skin surface, and the vertical line distance on midpoint to the bottom edge of the cut-off arm of angle is that body of gland compensates radius, further according to
Electromagnetic wave calculates the compensation time in the propagation speed differential of mammary fat tissue and gland tissue, carries out accurately image.
This method models the gland tissue of complex distributions in breast by MRI image, is reconstructed a smooth body of gland
Structure, and the surface details by adjusting body of gland enable the overall variation of the coarse gland surface of shape approximate description of reconstruct to become
Gesture, the new gland structure of construction are used for the volume of approximate evaluation body of gland.By the body of gland shapometer evaluation time propagated error of reconstruct,
Delay compensation is carried out in imaging algorithm, is finally accurately positioned the position of tumor of breast.
By delay compensation, image space deviation caused by when can effectively reduce imaging, especially to gland tissue point
The closeer breast of cloth, imaging accuracy raising are of great significance, and can be the direct positioning tumor of Minimally Invasive Surgery.It is now assumed that one suffers from
Knub position in person's body is x=15mm, y=25mm, z=60mm, and such as Fig. 9 (a), (b) shown respectively patient, which does not use, prolongs
The image of slow correction algorithm and the image using Deferred Correction algorithm.As can be seen that tumor image position deviation by it is original-
4mm's changes into+0.5mm, largely has modified the ultrafast image error that tumor of breast is detected with microwave.Figure 10 is to 25
Example patient compares without using before this algorithm and using the imaging positioning after this algorithm, it is seen that imaging deviation is integrally obviously reduced, body
The practicability of method is showed.
Description of the drawings
Fig. 1 breast MRI slice maps
Udder texture figure derived from Fig. 2 MRI
Fig. 3 three-dimensional gland surface structures
Fig. 4 space coordinate transformations
Gland tissue after Fig. 5 space coordinate transformations
Fig. 6 gland surface roughness profiles calculate
Fig. 7 ultra wide bands detect mammary tumor model and antenna configuration
Fig. 8 compensates radium computing method
The image (b) that Fig. 9 (a) is not suitable for backoff algorithm uses the image of backoff algorithm
25 patients of Figure 10 are imaged drift correction effect
Specific implementation method
The present invention will be described with reference to the accompanying drawings and examples.
(1) Fig. 1 show breast MRI slice maps.It is primarily based on the walk-off-mode of the MRI slice atlas structure bodies of gland of breast
Type, by the equidistant heap poststack of MRI pictures, gap between slice and slice into row interpolation, construct one it is complete it is three-dimensional from
Breast model is dissipated, threshold definitions are carried out by the gray value interval to MRI slice maps, internal mammary are divided into a variety of different
Tissue, each tissue are as shown in Figure 2.Directly the body of gland part in discrete model is extracted, as shown in Figure 3.
(2) gland structure is modeled, keeps the model of structure approximate with former gland structure, by comparing a large amount of mammary gland
MRI image figure, is elliptic cone shape by the structure proximate of gland tissue in breast, which can preferably describe absolutely mostly
Gland structure in number female breast.The vertex of oval wimble structure is chosen first, in the bottom surface for then choosing oval wimble structure again
Heart approximate location thus achieves the central symmetry axis of elliptic cone.The symmetry axis of selection is as shown in Figure 4.
(3) due to the uncertainty that body of gland is grown in human breast, the spatial position of elliptic cone and towards random.Therefore by scheming
Shown in 4 and Fig. 5, use space coordinate transform keeps the same coordinate in its bottom surface flat by elliptic cone " righting " to a secondary coordinate system
Face is parallel.
(4) an isosceles triangle region is created in the plane projection region of elliptic cone, as shown in Figures 5 and 6.And by this
Region finite element, each grid are denoted as pi, the half of the wherein bottom edge length of triangle is the long axis a of elliptic cone, elliptic cone
Height h and a remain unchanged, adjust the size of short axle b, make elliptic cone rise in the both direction up and down of gore or
Decline.When b takes a certain length, p is calculated firstiPoint upper circular conical surface with highest body of gland pixel difference absolute value, then
Calculate piThe two absolute values are included in a sample by the lower circular conical surface of point with the absolute value of the difference of minimum body of gland pixel,
Wait for all piPoint completes same individual of sample evaluation work, calculates sample variance, and it is most to take b values when sample variance minimum
Suitable minor axis dimension, the variation tendency of the best description body of gland rough surface of elliptic cone energy at this time.
(5) microwave propagation time estimation error is carried out.As shown in fig. 7, acquiescence ultra-wideband antenna position is first half
Breast skin surface.As shown in figure 8, the vertical line distance on the midpoint of the cut-off arm of angle to bottom edge is that body of gland compensates radius.Further according to electricity
Magnetic wave calculates the compensation time in the propagation speed differential of mammary fat tissue and gland tissue.
Claims (1)
1. a kind of ultra wide band tumor of breast imaging method based on magnetic resonance image compensation, includes the following steps：
(1) the 3 d-dem breast model of the MRI slice atlas structure bodies of gland based on breast, passes through the gray scale to MRI slice maps
It is worth section and carries out threshold definitions, internal mammary is divided into a variety of different tissues, by the body of gland extracting section in discrete model
Out；
(2) gland structure is modeled, keeps the model of structure approximate with former gland structure, by the structure of gland tissue in breast
It is approximately elliptic cone shape, chooses the vertex of oval wimble structure first, then choose the bottom center position of oval wimble structure again,
Obtain the central symmetry axis of elliptic cone；
(3) use space coordinate transform keeps the same coordinate plane in its bottom surface flat by elliptic cone " righting " to a secondary coordinate system
Row；
(4) an isosceles triangle region is created in the plane projection region of elliptic cone, and by this region finite element, each net
Lattice are denoted as pi, wherein the half of the bottom edge length of isosceles triangle is the long axis a of elliptic cone, the height h and long axis a of elliptic cone
It remains unchanged, adjusts the size of short axle b, so that elliptic cone is risen or fallen in the both direction up and down of gore, when short
When axis b takes a certain length, grid p is calculated firstiThe upper circular conical surface of point then is counted with the absolute value of the difference of highest body of gland pixel
Calculate grid piLower circular conical surface with minimum body of gland pixel difference absolute value, the two absolute values are included in a sample,
Wait for all grid piAfter point completes same individual of sample evaluation work, sample variance is calculated, when taking sample variance minimum
B values be most suitable minor axis dimension, the variation tendency for the description body of gland rough surface that elliptic cone at this time can be best；
(5) microwave propagation time estimation error is carried out, accurately image is carried out, method is as follows：Give tacit consent to ultra-wideband antenna position
For preceding half of breast skin surface, determine that body of gland compensates radius, further according to electromagnetic wave in mammary fat tissue and gland tissue
Propagation speed differential calculates the compensation time, carries out accurately image.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610019122.4A CN105701814B (en) | 2016-01-13 | 2016-01-13 | A kind of ultra wide band tumor of breast imaging method based on magnetic resonance image compensation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610019122.4A CN105701814B (en) | 2016-01-13 | 2016-01-13 | A kind of ultra wide band tumor of breast imaging method based on magnetic resonance image compensation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105701814A CN105701814A (en) | 2016-06-22 |
| CN105701814B true CN105701814B (en) | 2018-10-16 |
Family
ID=56226328
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610019122.4A Expired - Fee Related CN105701814B (en) | 2016-01-13 | 2016-01-13 | A kind of ultra wide band tumor of breast imaging method based on magnetic resonance image compensation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105701814B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106959306B (en) * | 2017-05-11 | 2019-11-12 | 天津大学 | A method of mammary tumor imaging is carried out using microwave reflection time domain S21 signal |
| CN107392902A (en) * | 2017-07-24 | 2017-11-24 | 重庆青囊科技有限公司 | A kind of tumour crack analysis method |
| CN116485789B (en) * | 2023-06-16 | 2023-08-25 | 新创碳谷集团有限公司 | Method, equipment and storage medium for detecting carbon fiber splitting defect |
Citations (3)
| 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 |
| CN103549953A (en) * | 2013-10-25 | 2014-02-05 | 天津大学 | Method for extracting microwave detection breast model based on medical magnetic resonance imaging |
| CN104473617A (en) * | 2014-11-10 | 2015-04-01 | 南方科技大学 | Biological tissue detection device, system and method |
-
2016
- 2016-01-13 CN CN201610019122.4A patent/CN105701814B/en not_active Expired - Fee Related
Patent Citations (3)
| 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 |
| CN103549953A (en) * | 2013-10-25 | 2014-02-05 | 天津大学 | Method for extracting microwave detection breast model based on medical magnetic resonance imaging |
| CN104473617A (en) * | 2014-11-10 | 2015-04-01 | 南方科技大学 | Biological tissue detection device, system and method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105701814A (en) | 2016-06-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Liu et al. | A parametric level set-based approach to difference imaging in electrical impedance tomography | |
| US10925582B2 (en) | Method and device for selecting detection area, and elasticity detection system | |
| CN103549953B (en) | Method for extracting microwave detection breast model based on medical magnetic resonance imaging | |
| CN105701814B (en) | A kind of ultra wide band tumor of breast imaging method based on magnetic resonance image compensation | |
| Stasi et al. | Emptying the rectum before treatment delivery limits the variations of rectal dose–volume parameters during 3DCRT of prostate cancer | |
| US10072924B2 (en) | System and method for quantifying deformation, disruption, and development in a sample | |
| CN104408398B (en) | A kind of recognition methods and system of liver boundary | |
| Williams et al. | Laser surface estimation for microwave breast imaging systems | |
| US10746822B2 (en) | System and method for localized processing of quantitative susceptibility maps in magnetic resonance imaging | |
| CN102944872B (en) | Near field-to-near field transformation method of radar scattering cross section | |
| Fontanarosa et al. | A speed of sound aberration correction algorithm for curvilinear ultrasound transducers in ultrasound-based image-guided radiotherapy | |
| O'Shea et al. | Temporal regularization of ultrasound‐based liver motion estimation for image‐guided radiation therapy | |
| Wentz et al. | Accuracy of dynamic patient surface monitoring using a time-of-flight camera and B-spline modeling for respiratory motion characterization | |
| CN103760544B (en) | The scattering centers extraction method and system of radar target | |
| O’Shea et al. | 4D ultrasound speckle tracking of intra-fraction prostate motion: a phantom-based comparison with x-ray fiducial tracking using CyberKnife | |
| Heethuis et al. | Quantification of variations in intra-fraction motion of esophageal tumors over the course of neoadjuvant chemoradiotherapy based on cine-MRI | |
| CN106920242A (en) | A kind of liver volume measuring method based on ultrasonoscopy and threedimensional model | |
| Sato et al. | Registration of temporal sequences of coronal and sagittal MR images through respiratory patterns | |
| Zykov et al. | Vector method of strain estimation in OCT-elastography with adaptive choice of scale for estimating interframe phase-variation gradients | |
| Zhang et al. | Surface extraction with a three-dimensional freehand ultrasound system | |
| KR101307162B1 (en) | Computer Aided System and method for Breast Mass Detection and Analysis in Microwave Tomography | |
| Ma et al. | A level set based algorithm to reconstruct the urinary bladder from multiple views | |
| Umay et al. | Endoscopic capsule localization with unknown signal propagation coefficients | |
| De Boer et al. | Influence of the number of elongated fiducial markers on the localization accuracy of the prostate | |
| Sun et al. | Anatomic surface reconstruction from sampled point cloud data and prior models |
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: 20181016 Termination date: 20220113 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |