CN115736987B - A fast imaging method for breast positioning in a hemispherical ultrasound tomography system - Google Patents
A fast imaging method for breast positioning in a hemispherical ultrasound tomography system Download PDFInfo
- Publication number
- CN115736987B CN115736987B CN202211458517.6A CN202211458517A CN115736987B CN 115736987 B CN115736987 B CN 115736987B CN 202211458517 A CN202211458517 A CN 202211458517A CN 115736987 B CN115736987 B CN 115736987B
- Authority
- CN
- China
- Prior art keywords
- breast
- data
- images
- attenuation
- coronal
- 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.)
- Active
Links
- 210000000481 breast Anatomy 0.000 title claims abstract description 39
- 238000003384 imaging method Methods 0.000 title claims abstract description 27
- 238000003325 tomography Methods 0.000 title claims abstract description 12
- 238000002604 ultrasonography Methods 0.000 title claims description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 12
- 238000000034 method Methods 0.000 claims abstract description 10
- 238000007781 pre-processing Methods 0.000 claims abstract description 4
- 239000011159 matrix material Substances 0.000 claims description 4
- 238000013507 mapping Methods 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000012163 sequencing technique Methods 0.000 claims 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
The invention belongs to the field of ultrasonic tomography, and in particular relates to a rapid imaging method for breast positioning of a hemispherical ultrasonic tomography system, which comprises the steps of acquiring empty water data; the method comprises the steps of collecting breast data, analyzing the collected data, preprocessing the analyzed data, calculating attenuation coefficients, reconstructing images, processing the reconstructed images, storing the processed images in a designated memory position, periodically reading the images through a GUI end, displaying breast position images of a coronal plane and a sagittal plane through the GUI end, and judging whether the breast position of a patient is in an optimal area or not. The invention can image a pair of coronal planes and a pair of sagittal planes in real time, can judge whether the breast is in the optimal imaging area in the XY direction through the coronal planes, can judge whether the patient is prone to be deep enough in the Z direction through the sagittal planes, is integrated on a GUI operation interface, displays images in the GUI in real time, is convenient to operate and judges whether the breast position of the patient is in the optimal area.
Description
Technical Field
The invention belongs to the field of ultrasonic tomography, and particularly relates to a rapid imaging method for breast positioning of a hemispherical ultrasonic tomography system.
Background
The imaging space of ultrasonic tomography is three-dimensional, but the optimal imaging area is generally smaller than the imaging aperture size. When a patient is prone on the device, the naked eye cannot well determine whether the breast is within the optimal imaging region. The existing optical camera mode can not accurately judge the breast position of a patient and has large imaging effect difference.
The Chinese patent with publication number CN113598825A discloses a breast positioning imaging method for an ultrasonic imaging system and application thereof, the scheme comprises the following steps of setting a sequence of ultrasonic transducer array elements, transmitting and collecting ultrasonic signals according to the sequence by each ultrasonic transducer array element, calculating and obtaining a contour of a layer of coronal surface of a breast according to the collected ultrasonic signals, and carrying out space coordinate conversion on the contour to obtain the space distribution of the contour in a cross section of the ultrasonic imaging system.
However, the breast positioning imaging method described above still has the following drawbacks:
the contrast is poor, the noise is obvious, the position of the breast is not easy to be effectively distinguished, only one section is displayed (X, Y directions), and no depth direction information exists.
Disclosure of Invention
In order to make up for the defects of the prior art, the invention provides a technical scheme of a rapid imaging method for breast positioning of a hemispherical ultrasonic tomography system.
A method for rapid imaging for breast positioning in a hemispherical ultrasound tomography system, comprising:
Step 1, acquiring empty water data;
Step 2, breast data acquisition;
step 3, analyzing the collected data;
step 4, preprocessing the analyzed data;
Step 5, calculating an attenuation coefficient;
step 6, reconstructing an image;
and 7, processing the reconstructed image, storing the processed image in a designated memory position, and periodically reading the image through a GUI (graphical user interface) end.
Further, in the step 1 and the step 2, data of the coronal plane and the sagittal plane are acquired simultaneously each time.
Further, the step 3 includes:
And corresponding the transmitting TAS sequence number, the transmitting array element sequence number, the receiving TAS sequence number and a section AScan, and storing the corresponding receiving array element sequence number and the corresponding section AScan in a memory.
Further, the step 3 includes distinguishing sagittal plane data from coronal plane data for reconstruction of sagittal plane and coronal plane images, respectively.
Further, the step 4 includes:
step 4.1, the energy of each segment AScan signal is calculated according to the following formula:
wherein L is the characteristic length of a transmitted signal, and t0 is the arrival time of a direct wave;
Step 4.2, LOR ordering comprises connecting all pairs of transmitting array elements and receiving array elements, solving the result of polar coordinate transformation of the Line under the original rectangular coordinate system, arranging signals into a (R, theta) matrix, wherein R is the distance from a path to a central point, theta is the included angle between the path and a coordinate axis, and LOR ordering is carried out on the data of a sagittal plane and a coronal plane separately.
Further, the step 5 includes:
According to the energy I calculated in the step 4, the ordered empty water and breast data are calculated according to the following formula:
Atten=log(I Empty water /I Object )
wherein Atten is the attenuation coefficient for each transmit-receive pair.
Further, the step 6 includes:
reconstructing the processed Atten according to the following formula:
Where f is the image, k is the number of iterations, n is the pixel number, w mn is the intersection distance of the mth attenuation and the nth pixel, p m is the mth attenuation value, Is the attenuation value of the mth weighted projection, W m+ is the normalization factor corresponding to the mth attenuation, lambda k is the coefficient of the kth iteration,Is the sampling point of the whole imaging system, w m'n is the sum of the pixel n and all the attenuation intersecting distances, and m' is used to calculate the sum of the pixel n and all the attenuation intersecting distances, and is a constant weighting, which can be proposed.
Further, the step 7 includes:
step 7.1, mapping the reconstructed image value to a 0-255 interval;
And 7.2, respectively adding a Mask with a rectangular frame to the images reconstructed from the sagittal plane and the coronal plane according to the theoretical optimal imaging area, reserving the images in the Mask, setting the images outside the Mask to 0, storing the processed images in a designated memory position, and periodically reading through a GUI end.
Further, the method also comprises the step 8 of displaying the breast position images of the coronal plane and the sagittal plane through two windows at the GUI end respectively for judging whether the breast position of the patient is in the optimal area.
Compared with the prior art, the invention has the beneficial effects that:
The invention can image a pair of coronal planes and a pair of sagittal planes in real time, can judge whether the breast is in the optimal imaging area in the XY direction through the coronal planes, can judge whether the patient is prone to be deep enough in the Z direction through the sagittal planes, is integrated on a GUI operation interface, displays images in the GUI in real time, is convenient to operate and judges whether the breast position of the patient is in the optimal area.
Drawings
FIG. 1 is a flow chart of the present invention.
Detailed Description
The invention is further described below with reference to the accompanying drawings.
The hemispherical support adopted by the invention refers to patent CN113598825A, and the ultrasonic transducer array element are provided with serial numbers according to the method in patent CN 113598825A.
Referring to fig. 1, a method for rapid imaging for breast positioning in a hemispherical ultrasound tomography system includes the steps of:
Step 1, acquiring empty water data
The empty water data acquisition is to acquire data only when water does not have an object, and the function is to provide a reference for the invention.
The selection of the transmitting and receiving transducers follows the principle of geometric spatial location by selecting TAS that is intersected by the hemispherical transducer. At each acquisition, coronal and sagittal data are acquired simultaneously.
Step 2, breast data acquisition
And acquiring data of the breast, wherein the adopted ultrasonic transducer and the acquisition mode are the same as those of the breast.
Step 3, data analysis
The method mainly comprises the steps of corresponding a transmitting TAS sequence number, a transmitting array element sequence number, a receiving TAS sequence number, a receiving array element sequence number and a section AScan (4000 points), and storing the corresponding sequences in a memory. Meanwhile, the step distinguishes the sagittal plane and the coronal plane and is respectively used for reconstructing sagittal plane and coronal plane images.
Step 4, data preprocessing
Step 4.1, the energy of each segment AScan signal is calculated and integrated according to the following formula:
where L is the characteristic length of the transmitted signal and t0 is the direct wave arrival time.
Step 4.2, LOR (Line of Response) reorder. The specific way is to reorder the signals according to their geometric position. And connecting all transmitting array elements and receiving array element pairs, solving the polar coordinate transformation of the Line under the original direct coordinate system, and arranging signals into an (R, theta) matrix. R is the distance from the path to the center point, theta is the included angle between the path and the coordinate axis. In this step the data of the sagittal and coronal planes need to be sorted separately for LOR.
Step 5, calculating attenuation coefficient
Based on the energy I calculated in step 4, the ordered empty water and breast data are calculated according to the following formula:
Atten=log(I Water and its preparation method /I Object )
Atten is the attenuation coefficient of each transmitting-receiving pair.
Step 6, image reconstruction
The processed Atten was reconstructed according to the following formula.
Where f is the image, k is the number of iterations, n is the pixel number, w mn is the intersection distance of the mth attenuation and the nth pixel, p m is the mth attenuation value,Is the attenuation value of the mth weighted projection, W m+ is the normalization factor corresponding to the mth attenuation, lambda k is the coefficient of the kth iteration,Is the sampling point of the whole imaging system, m' is used to calculate the sum of the pixel n and all attenuation intersection distances, and is a constant weighting, which can be proposed.
Step 7, image processing
And 7.1, mapping the reconstructed image value to a 0-255 interval.
And 7.2, respectively adding a Mask with a rectangular frame to the images reconstructed from the sagittal plane and the coronal plane according to the theoretical optimal imaging area, wherein the Mask is a two-dimensional matrix with an optimal area of 1 and an external area of 0, and the images in the Mask are reserved and the images are set to 0. The processed image is stored in a designated memory location, and the GUI end periodically reads.
Step 8, image GUI display
The GUI side (graphical user interface) has two windows showing images of breast positions in the coronal and sagittal planes, respectively, for determining whether the patient's breast position is in the optimal region. If the breast is not in the optimal imaging region, the patient is instructed to move and adjust the breast position until the breast is in the optimal region.
It should be noted that the above embodiments are merely for illustrating the technical solution of the present invention and not for limiting the same, and although the present invention has been described in detail with reference to the above embodiments, it should be understood by those skilled in the art that the technical solution described in the above embodiments may be modified or some or all of the technical features may be equivalently replaced, and these modifications or substitutions do not make the essence of the corresponding technical solution deviate from the scope of the technical solution of the embodiments of the present invention.
Claims (2)
1. A method for rapid imaging for breast positioning in a hemispherical ultrasound tomography system, comprising:
step 1, acquiring empty water data, wherein the data of a coronal plane and a sagittal plane are acquired simultaneously when each acquisition is performed;
Step 2, breast data are collected, and data of a coronal plane and a sagittal plane are collected simultaneously when each time of collection is performed;
and 3, analyzing the acquired data, wherein the step comprises the following steps:
The method comprises the steps of storing a transmitting TAS sequence number, a transmitting array element sequence number, a receiving TAS sequence number, a receiving array element sequence number and a section AScan in a memory, distinguishing sagittal plane data from coronal plane data, and respectively reconstructing sagittal plane and coronal plane images;
Step 4, preprocessing the analyzed data, including:
step 4.1, the energy of each segment AScan signal is calculated according to the following formula:
wherein L is the characteristic length of a transmitted signal, and t0 is the arrival time of a direct wave;
Step 4.2, LOR sequencing, which comprises connecting all pairs of transmitting array elements and receiving array elements, solving the result of polar coordinate transformation of the connecting lines under the original rectangular coordinate system, arranging signals into a (R, theta) matrix, wherein R is the distance from a path to a central point, theta is the included angle between the path and a coordinate axis, and the LOR sequencing is carried out on the data of a sagittal plane and a coronal plane separately;
step 5, calculating attenuation coefficient, including:
According to the energy I calculated in the step 4, the ordered empty water and breast data are calculated according to the following formula:
Atten=log(I Empty water /I Object )
wherein Atten is the attenuation coefficient of each transmit-receive pair;
step6, image reconstruction is carried out, which comprises the following steps:
reconstructing the processed Atten according to the following formula:
Where f is the image, k is the number of iterations, n is the pixel number, w mn is the intersection distance of the mth attenuation and the nth pixel, p m is the mth attenuation value, Is the attenuation value of the mth weighted projection, W m+ is the normalization factor corresponding to the mth attenuation, lambda k is the coefficient of the kth iteration,Is the sampling point of the whole imaging system, w m'n is the sum of the pixel n and all attenuation intersecting distances, m' is used for calculating the sum of the pixel n and all attenuation intersecting distances, and is a constant weighting;
and 7, processing the reconstructed image, and storing the processed image in a designated memory position, wherein the method comprises the following steps:
step 7.1, mapping the reconstructed image value to a 0-255 interval;
And 7.2, respectively adding a Mask with a rectangular frame to the images reconstructed from the sagittal plane and the coronal plane according to the theoretical optimal imaging area, reserving the images in the Mask, setting the images outside the Mask to 0, storing the processed images in a designated memory position, and periodically reading through a GUI end.
2. The method for rapid imaging for breast positioning in a hemispherical ultrasound tomography system of claim 1, further comprising step 8, wherein step 8 comprises displaying the breast position images of the coronal and sagittal planes through two windows at the GUI end, respectively, for determining whether the patient's breast position is in the optimal region.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211458517.6A CN115736987B (en) | 2022-11-21 | 2022-11-21 | A fast imaging method for breast positioning in a hemispherical ultrasound tomography system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211458517.6A CN115736987B (en) | 2022-11-21 | 2022-11-21 | A fast imaging method for breast positioning in a hemispherical ultrasound tomography system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115736987A CN115736987A (en) | 2023-03-07 |
| CN115736987B true CN115736987B (en) | 2025-02-18 |
Family
ID=85333954
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211458517.6A Active CN115736987B (en) | 2022-11-21 | 2022-11-21 | A fast imaging method for breast positioning in a hemispherical ultrasound tomography system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115736987B (en) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111194185A (en) * | 2017-10-02 | 2020-05-22 | 百合医疗科技株式会社 | Medical imaging device |
| CN113598825A (en) * | 2021-09-16 | 2021-11-05 | 浙江衡玖医疗器械有限责任公司 | Breast positioning imaging method for ultrasonic imaging system and application thereof |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2010268944A (en) * | 2009-05-21 | 2010-12-02 | Konica Minolta Medical & Graphic Inc | Ultrasonic diagnostic imaging method and apparatus |
| US8663113B2 (en) * | 2011-08-11 | 2014-03-04 | Delphinus Medical Technologies, Inc. | Method for imaging a volume of tissue |
| CN102688071A (en) * | 2012-06-15 | 2012-09-26 | 华东医院 | Ultrasonic superficial tissue and organ volume scanning fracture imaging method |
| JP7085470B2 (en) * | 2015-09-01 | 2022-06-16 | デルフィヌス メディカル テクノロジーズ， インコーポレイテッド | Tissue imaging and analysis using ultrasonic waveform tomography |
| CN107616784B (en) * | 2017-09-19 | 2020-09-22 | 华南师范大学 | A wide-field photoacoustic ultrasound breast coronal scanning imaging device and method based on 1024 linear array detection |
| KR102770978B1 (en) * | 2018-04-04 | 2025-02-19 | 토모웨이브 래버러토리즈, 인코포레이티드 | Quantitative imaging system and its uses |
| CN108830852B (en) * | 2018-07-13 | 2021-08-27 | 上海深博医疗器械有限公司 | Three-dimensional ultrasonic tumor auxiliary measurement system and method |
-
2022
- 2022-11-21 CN CN202211458517.6A patent/CN115736987B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111194185A (en) * | 2017-10-02 | 2020-05-22 | 百合医疗科技株式会社 | Medical imaging device |
| CN113598825A (en) * | 2021-09-16 | 2021-11-05 | 浙江衡玖医疗器械有限责任公司 | Breast positioning imaging method for ultrasonic imaging system and application thereof |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115736987A (en) | 2023-03-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5511892B2 (en) | Ultrasonic image processing system using pixel-oriented processing | |
| US7087022B2 (en) | 3D ultrasound-based instrument for non-invasive measurement of amniotic fluid volume | |
| JP5849048B2 (en) | Three-dimensional (3D) ultrasound imaging system for scoliosis evaluation | |
| US7041059B2 (en) | 3D ultrasound-based instrument for non-invasive measurement of amniotic fluid volume | |
| US7744534B2 (en) | 3D ultrasound-based instrument for non-invasive measurement of amniotic fluid volume | |
| CN100998511A (en) | Real-time, freedom-arm, three-D ultrasonic imaging system and method therewith | |
| US20080242985A1 (en) | 3d ultrasound-based instrument for non-invasive measurement of amniotic fluid volume | |
| WO2020206755A1 (en) | Ray theory-based method and system for ultrasound ct image reconstruction | |
| CN104268412B (en) | Angle gather ray chromatography migration velocity analysis method and device | |
| CA2534287A1 (en) | Image enhancing and segmentation of structures in 3d ultrasound | |
| RU2677191C2 (en) | Rib blockage delineation in anatomically intelligent echocardiography | |
| CN107106128A (en) | Supersonic imaging device and method for splitting anatomical object | |
| CN116058867B (en) | A lightweight imaging ultrasound system with unfixed scanning probe and an ultrasound detection method | |
| Gao et al. | Wireless and sensorless 3D ultrasound imaging | |
| CN103761767A (en) | Quick three-dimensional ultrasound image reconstruction method based on sparse data | |
| Yuan et al. | Optimization of reconstruction time of ultrasound computed tomography with a piecewise homogeneous region-based refract-ray model | |
| CN107928708A (en) | Freely three-dimensional backbone ultrasonic image-forming system and control method | |
| CN115736987B (en) | A fast imaging method for breast positioning in a hemispherical ultrasound tomography system | |
| CN107690312A (en) | Supersonic imaging device | |
| CN111062935B (en) | Mammary gland tumor detection method, storage medium and terminal equipment | |
| CN114521914A (en) | Ultrasonic parameter measuring method and ultrasonic parameter measuring system | |
| CN107076820A (en) | Method for assessing and improving the quality of data in fine-structure distribution data | |
| CN118383798B (en) | Method and system for identifying throat space occupying lesions based on image identification | |
| CN116777879A (en) | Prostate ultrasound image auxiliary diagnosis method and system based on YOLOv5 | |
| CN118236091A (en) | Ultrasound imaging system and method for calculating and displaying probe position adjustment |
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 |