CN113552573B - Rapid imaging algorithm based on ultrasonic ring array synthetic aperture receiving - Google Patents
Rapid imaging algorithm based on ultrasonic ring array synthetic aperture receiving Download PDFInfo
- Publication number
- CN113552573B CN113552573B CN202110725297.8A CN202110725297A CN113552573B CN 113552573 B CN113552573 B CN 113552573B CN 202110725297 A CN202110725297 A CN 202110725297A CN 113552573 B CN113552573 B CN 113552573B
- Authority
- CN
- China
- Prior art keywords
- imaging
- array
- ultrasonic
- synthetic aperture
- algorithm
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8997—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using synthetic aperture techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8909—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration
- G01S15/8915—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array
- G01S15/8922—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array the array being concentric or annular
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
The invention belongs to the field of ultrasonic detection and imaging, and provides a rapid imaging algorithm based on ultrasonic ring array synthetic aperture receiving, which uses a ring array ultrasonic transducer to transmit and receive ultrasonic signals in a synthetic aperture receiving mode, groups array elements of the ring array ultrasonic transducer according to n adjacent array elements of each group, and calculates an index matrix D of a sector area formed by 1-n array elements of the 1 st group and the circle center of the ring array ultrasonic transducer 1 (ii) a Setting an initial value i to be 1; obtaining corresponding local imaging result I by utilizing synthetic aperture algorithm i (ii) a Obtaining the ith group index matrix D through rotation transformation i Repeatedly executing the first two steps until all local imaging results are obtained; and superposing and fusing all local imaging results according to the positions of the local imaging results to obtain a final imaging result. The algorithm provided by the invention can effectively remove arc artifacts in the traditional synthetic aperture imaging algorithm and greatly improve the performance of the algorithm.
Description
Technical Field
The invention belongs to the field of ultrasonic detection and imaging, and particularly relates to a rapid imaging algorithm based on ultrasonic ring array synthetic aperture receiving.
Background
The synthetic aperture technology is used for radar detection at first, compared with the traditional fixed point focusing technology, the ultrasonic synthetic aperture imaging technology adopts dynamic focusing, a transducer of a small aperture array element and lower central frequency can be used for realizing imaging with high azimuth resolution, and the imaging quality in an area to be imaged is integrally improved. The ultrasonic synthetic aperture imaging technology is applied to intravascular imaging, liver lesion imaging and the like. For a ring array transducer, if a traditional full-aperture synthetic aperture imaging technology is used, the imaging result has the problems of inaccurate inner surface, excessive artifacts, long imaging time consumption and the like.
Disclosure of Invention
The invention is made to solve the above problems, and an object of the invention is to provide a fast imaging algorithm based on ultrasound circular array synthetic aperture reception with less artifacts, short imaging time, and high efficiency.
The invention provides a rapid imaging algorithm based on ultrasonic ring array synthetic aperture receiving, which is characterized by comprising the following steps: step S1, dividing N array elements of the ultrasonic annular array transducer into N array element groups according to each group of N adjacent array elements, and calculating an index matrix D of a 1 st sector area consisting of a 1 st array element group containing array elements 1-N and the circle center of the ultrasonic annular array transducer 1 (ii) a Step S2, setting the cycle value i to 1; step S3, according to the ultrasonic echo signals received by the array elements i to (n + i-1), calculating each array element in the array elements i to (n + i-1) to an index matrix D i The ultrasound propagation time of the corresponding point; step S4, calculating index matrix D by using synthetic aperture algorithm according to ultrasonic echo signal and ultrasonic propagation time i Imaging result I corresponding to sector area i (ii) a Step S5, let i equal i +1, and index matrix D 1 Push buttonRotating and transforming to obtain an index matrix D i Then returns to step S3 until i > N; step S6, imaging result I of the sector area i And N is superposed and fused according to the position of the N, so as to obtain the final imaging result of the sample.
In the fast imaging algorithm based on the ultrasonic circular array synthetic aperture receiving provided by the invention, the fast imaging algorithm can also have the following characteristics: wherein, index matrix D i The imaging device is provided with a plurality of elements which are in one-to-one correspondence with points in the whole to-be-imaged area, when the value of the element is 1, the point is located in the ith sector area, and when the value of the element is 0, the point is located outside the ith sector area.
In the fast imaging algorithm based on the ultrasonic circular array synthetic aperture receiving provided by the invention, the fast imaging algorithm can also have the following characteristics: in step S4, the index matrix D in the i-th sector area i Any point P (x, z) in which the value of the corresponding element in (b) is 1, according to the following formula:
calculating the image brightness value of the point P (x, z), wherein the matrix formed by the image brightness values of all the points in the ith fan-shaped area is the imaging result I i In the formula, s j,k (t) is the amplitude of the signal transmitted by the array element k and received by the array element j at the time t, t j (x, z) and t k (x, z) are the ultrasound propagation times to point P for element j and element k, respectively.
In the fast imaging algorithm based on the ultrasonic circular array synthetic aperture receiving provided by the invention, the fast imaging algorithm can also have the following characteristics: in step S5, a first arrival wave may be extracted from the ultrasonic echo signal by using a akage information criterion algorithm, then a sound velocity distribution model in the region to be imaged is obtained by using an inversion algorithm, and finally the sound velocity distribution model is substituted into an equation of function to solve the ultrasonic propagation time.
Action and Effect of the invention
According to the rapid imaging algorithm based on the ultrasonic ring array synthetic aperture receiving, the region to be imaged is divided into N fan-shaped regions, the imaging result of each fan-shaped region is obtained by utilizing the synthetic aperture algorithm, the imaging result in each fan-shaped region is ensured to be accurate, and finally the N fan-shaped imaging results are overlapped and fused according to the positions of the N fan-shaped imaging results, so that the arc-shaped artifacts of the obtained imaging result are obviously reduced, and the final imaging result is more accurate and clear. In addition, because the method of firstly calculating the index matrix of one sector area and then obtaining the index matrices of the other sector areas through rotation transformation is adopted, compared with the operation of performing complete index value calculation once on each array element group, the method effectively improves the algorithm performance.
Drawings
Fig. 1 is a schematic diagram of a signal transmitting and acquiring system of an ultrasonic circular array transducer in an embodiment of the present invention, and fig. 2 is a flowchart of a fast imaging algorithm based on ultrasonic circular array synthetic aperture receiving in an embodiment of the present invention;
FIG. 3 is a sector area schematic diagram of an array element grouping of an ultrasonic loop array transducer in an embodiment of the invention;
FIG. 4 is a schematic diagram of an imaging result and a final imaging result of an array element group of a sample to be imaged by using a fast imaging algorithm based on ultrasonic ring array synthetic aperture receiving in an embodiment of the present invention;
FIG. 5 is a schematic diagram of an imaging result and a final imaging result of an array element group when a sample to be imaged uses a conventional full aperture synthetic aperture imaging algorithm according to an embodiment of the present invention; and
fig. 6 is a graph illustrating the performance improvement of the fast imaging algorithm based on the ultrasound circular array synthetic aperture reception compared with the conventional synthetic aperture imaging algorithm in the embodiment of the present invention.
Detailed Description
In order to make the technical means, the creation characteristics, the achievement purposes and the effects of the invention easy to understand, the following describes a fast imaging algorithm based on ultrasound circular array synthetic aperture receiving specifically with reference to the embodiments and the accompanying drawings.
< example >
The fast imaging algorithm based on the ultrasonic circular array synthetic aperture receiving is realized based on a signal transmitting and acquiring system of an ultrasonic circular array transducer.
Fig. 1 is a schematic diagram of a signal transmitting and acquiring system of an ultrasonic loop array transducer in the embodiment.
As shown in fig. 1, the system has a water tank 1, an ultrasonic loop array transducer 2, a cortical bone phantom 3, and a signal transmitting and acquiring device 4.
The ultrasonic loop array transducer 2 is fixed in the water tank 1, the area to be imaged of the cortical bone phantom sample 3 is placed in the center of the ultrasonic loop array transducer 2, and water is injected into the water tank to immerse the sample and the ultrasonic loop array transducer 2. The signal transmitting and collecting device 4 controls all array elements of the ultrasonic loop array transducer to transmit ultrasonic pulse signals in sequence at a fixed central frequency, and receives ultrasonic echo signals of the ultrasonic pulse signals through all the array elements.
In this step, the supersound cyclic array transducer that uses is the supersound cyclic array transducer that contains 128 array elements, and the cyclic array diameter is 50mm, and array element center distance is 1.23mm, and the array element interval is 0.2mm, and central frequency is 3.5 MHz. The pulse signal is a Gaussian envelope sine wave with two periods, and the sampling frequency of the system is 25 MHz. After the ultrasonic echo signal acquisition is finished, the ultrasonic propagation time from each array element to any position in the region to be imaged can be calculated by the ultrasonic echo signal.
Fig. 2 is a flowchart of a fast imaging algorithm based on ultrasonic circular array synthetic aperture reception in this embodiment, and fig. 3 is a schematic view of a sector area of a grouping of ultrasonic circular array transducer elements in this embodiment.
As shown in fig. 2 and fig. 3, the fast imaging algorithm based on the ultrasound circular array synthetic aperture receiving includes the following steps:
step S1, dividing 128 array elements of the ultrasonic annular array transducer into 128 array element groups according to 8 adjacent array elements in each group, and calculating an index matrix D of a 1 st sector area consisting of a 1 st group of array element groups containing array elements 1-8 and the circle center of the ultrasonic annular array transducer 1 . The sector area formed by one array element is schematically shown as sector area a in fig. 3 (only three sector areas a, b, c are schematically shown in fig. 3, and the remaining sector areas are not shown).
Index matrix D 1 The imaging device is provided with a plurality of elements which are in one-to-one correspondence with points in the whole to-be-imaged area, when the value of the element is 1, the point is located in the 1 st sector area, and when the value of the element is 0, the point is located outside the 1 st sector area.
In step S2, the cycle value i is set to 1.
Step S3, according to the ultrasonic echo signals received by the array elements i to (8+ i-1), extracting a first arrival wave from the obtained ultrasonic echo signals by using a Chichi information criterion algorithm, then obtaining a sound velocity distribution model in a region to be imaged by using an inversion algorithm, finally substituting the sound velocity distribution model into an equation of function, and solving to obtain an index matrix D from each array element in the array elements i to (8+ i-1) i The ultrasound propagation time of the corresponding point.
In this step, other methods may also be used to calculate the ultrasound propagation time, for example, in the step of obtaining the sound velocity distribution model, the method may be replaced by "reconstructing the sound velocity distribution model of the region to be imaged by using the bayesian estimation method".
Step S4, calculating index matrix D by using synthetic aperture algorithm according to ultrasonic echo signal and ultrasonic propagation time i Imaging result I of the corresponding I-th sector area i The method comprises the following specific steps:
for the index matrix D in the ith sector i Any point P (x, z) in which the value of the corresponding element is 1, according to the following formula:
calculating the image brightness value of the point P (x, z), wherein the matrix formed by the image brightness values of all the points P (x, z) in the ith fan-shaped area is the imaging result I i ,
In the formula, s j,k (t) is the amplitude of the signal transmitted by the array element k and received by the array element j at the time t, t j (x, z) and t k (x, z) are the ultrasound propagation times of element j and element k to point P, respectively.
Step S5, let i equal i +1, and index matrix D 1 Push buttonRotating and transforming to obtain an index matrix D i Then returns to step S3 until i > 128;
step S6, imaging the sector area I i And the position of the sample is superposed and fused by 128 to obtain the final imaging result of the sample.
Fig. 4 is an imaging result of array element group of the sample to be imaged by using a fast imaging algorithm based on ultrasonic ring array synthetic aperture reception in the embodiment.
As shown in fig. 4, for the sector area imaging results as shown in a, b, and c in fig. 4 (only three sector area imaging results a, b, and c are schematically shown in fig. 4, and the remaining sector area imaging results are not shown), the final imaging results obtained by the fusion performed in step S6 are shown as d and e in fig. 4.
< comparative example >
In this comparative example, the sample to be imaged in the example was imaged using a conventional full aperture synthetic aperture imaging algorithm, and the imaging result was compared with that in the example.
Fig. 5 is a schematic diagram of the imaging result and the final imaging result of the array element group of the sample to be imaged in the present comparative example by using the conventional full aperture synthetic aperture imaging algorithm.
The sample to be imaged in the embodiment is imaged by using a conventional full aperture synthetic aperture imaging algorithm, and the imaging result of the array element group and the final imaging result are shown in fig. 5.
As shown in fig. 4 and 5, in the embodiment, the final imaging result of the fast imaging algorithm based on the ultrasound circular array synthetic aperture reception for the sample to be imaged is significantly reduced in arc artifact compared with the final imaging result of the conventional full aperture synthetic aperture imaging algorithm for the sample to be imaged.
Fig. 6 is a graph illustrating the performance improvement of the ultrasound loop array synthetic aperture reception-based fast imaging algorithm of the present invention compared to a comparative example conventional synthetic aperture imaging algorithm.
As shown in fig. 6, when N and N are different and have different values, the efficiency of the fast imaging algorithm based on the ultrasound circular array synthetic aperture reception is superior to that of the conventional synthetic aperture imaging algorithm. In this embodiment, when N is 128 and N is 8, the ratio of the computational efficiency of the fast imaging algorithm based on the ultrasound circular array synthetic aperture reception and the conventional synthetic aperture imaging algorithm is 3.7237.
The fast imaging algorithm based on the ultrasound circular array synthetic aperture receiving provided by this embodiment calculates the index matrix D of the sector area corresponding to the number 1-n array elements 1 Then, an index matrix D of a sector area corresponding to array elements i to (n + i-1) is obtained through rotation transformation i Compared with the operation of performing complete index value calculation once on each array element group in the traditional algorithm, the method effectively improves the performance of the algorithm. The time complexity of the conventional algorithm is The algorithm of the invention has the time complexity of Where Nx × Ny is the size of the region to be imaged, f is the time complexity of signal amplitude addition for one time when a luminance value of a pixel point in the image is calculated, g is the time complexity of judging whether a point in the imaging region is located in the sector region, h is the time complexity of performing rotation transformation on the matrix, and values of f, g, and h are O (1), O (n +1), and O (1), respectively.
Effects and effects of the embodiments
According to the rapid imaging algorithm based on the ultrasonic ring array synthetic aperture receiving provided by the embodiment, the region to be imaged is divided into N fan-shaped regions, the imaging result of each fan-shaped region is obtained by utilizing the synthetic aperture algorithm, the imaging result in each fan-shaped region is ensured to be accurate, and finally the N fan-shaped imaging results are overlapped and fused according to the positions of the N fan-shaped imaging results, so that the arc-shaped artifacts of the obtained imaging result are obviously reduced, and the final imaging result is more accurate and clear. In addition, because the method of calculating the index matrix of one sector area and then obtaining the index matrices of the other sector areas through rotation transformation is adopted, compared with the operation of performing complete index value calculation once on each array element group, the method effectively improves the algorithm performance.
Further, the rapid imaging algorithm based on the ultrasound circular array synthetic aperture receiving provided by the embodiment can adopt parallel computing to a plurality of sector areas to obtain imaging results, and the performance is improved.
In summary, the fast imaging algorithm based on the ultrasound circular array synthetic aperture receiving of the embodiment has a good imaging effect and greatly improved algorithm performance compared with the traditional full-aperture synthetic aperture imaging algorithm.
The above embodiments are preferred examples of the present invention, and are not intended to limit the scope of the present invention.
Claims (2)
1. A fast imaging algorithm based on ultrasonic ring array synthetic aperture receiving is used for imaging a sample by ultrasonic echo signals transmitted and received by an ultrasonic ring array transducer, and is characterized by comprising the following steps:
step S1, the ultrasonic ring array transducerEach array element in each groupEach adjacent said array element being divided intoAn array element group, the calculation comprises the array elementsThe 1 st sector area index matrix formed by the array element group and the circle center of the ultrasonic annular array transducer;
Step S3, according to the array elementThe received ultrasonic echo signal is used for calculating the array elementEach of whichThe array element to index matrixThe ultrasound propagation time of the corresponding point;
step S4, according to the ultrasonic echo signal and the ultrasonic propagation time, the index matrix is calculated by utilizing a synthetic aperture algorithmTo correspond to the firstImaging results of a sector;
Step S5, letAnd the index matrix is combinedPush buttonThe index matrix is obtained by rotating and transformingAnd then returns to the step S3 until;
Step S6, imaging the sector area Superposing and fusing the positions of the two to obtain the final imaging result of the sample,
in step S3, the first step,the index matrixHaving a plurality of elements in one-to-one correspondence with points in the entire region to be imaged, the value of said element being 1 indicating that said point is located in said second placeIn each sector, the value of the element is 0, which means that the point is located at the second positionOutside the area of the sector of the circle,
in the step S4, for the second stepWithin each sector area of the index matrixAt any point where the value of the corresponding element is 1According to the following formula:
calculate the pointThe image luminance value ofThe matrix formed by the image luminance values of all the points in the fan-shaped area is the imaging resultIn the formula (I), wherein,as an array elementTransmitting array elementThe received signal is atThe magnitude of the amplitude of the time of day,andare respectively the array elementsAnd said array elementToThe ultrasound propagation time of a point.
2. The fast imaging algorithm based on ultrasound circular array synthetic aperture reception of claim 1, wherein:
in step S3, a first arrival wave is extracted from the ultrasonic echo signal by using a akage information criterion algorithm, then a sound velocity distribution model in the region to be imaged is obtained by using an inversion algorithm, and finally the sound velocity distribution model is substituted into an engineering function equation to obtain the ultrasonic propagation time.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110725297.8A CN113552573B (en) | 2021-06-29 | 2021-06-29 | Rapid imaging algorithm based on ultrasonic ring array synthetic aperture receiving |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110725297.8A CN113552573B (en) | 2021-06-29 | 2021-06-29 | Rapid imaging algorithm based on ultrasonic ring array synthetic aperture receiving |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113552573A CN113552573A (en) | 2021-10-26 |
CN113552573B true CN113552573B (en) | 2022-07-29 |
Family
ID=78102460
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110725297.8A Active CN113552573B (en) | 2021-06-29 | 2021-06-29 | Rapid imaging algorithm based on ultrasonic ring array synthetic aperture receiving |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113552573B (en) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1998002761A1 (en) * | 1996-07-11 | 1998-01-22 | Science Applications International Corporation | Terrain elevation measurement by interferometric synthetic aperture radar (ifsar) |
CN105411626A (en) * | 2015-12-24 | 2016-03-23 | 华中科技大学 | Ultrasonic CT-based synthetic aperture imaging method and system |
CN109281651A (en) * | 2017-07-19 | 2019-01-29 | 中国科学院声学研究所 | A kind of ultrasonic borehole wall imaging method applied to cylinder supersonic array |
CN110179495A (en) * | 2019-04-23 | 2019-08-30 | 华中科技大学 | Ultrasonic tomography calculation optimization method and system based on distributed system |
CN110441398A (en) * | 2019-07-17 | 2019-11-12 | 复旦大学 | A kind of synthetic aperture ultrasonic imaging method based on multilayer dielectricity velocity of sound model |
CN110974293A (en) * | 2019-12-11 | 2020-04-10 | 华中科技大学 | Synthetic aperture imaging method based on C-type probe |
CN112138970A (en) * | 2020-09-22 | 2020-12-29 | 深圳市赛禾医疗技术有限公司 | Ultrasonic forward loop array transceiver |
WO2021004076A1 (en) * | 2019-07-05 | 2021-01-14 | 山东大学 | Ai chip-based conformal wearable biological information monitoring device and system |
CN112764040A (en) * | 2019-11-01 | 2021-05-07 | 复旦大学 | Synthetic aperture beam forming method based on ray theory phase correction |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1998012667A2 (en) * | 1996-08-29 | 1998-03-26 | Johnson Steven A | Wavefield imaging using inverse scattering techniques |
US6123671A (en) * | 1998-12-31 | 2000-09-26 | General Electric Company | Method and apparatus for distributed, agile calculation of beamforming time delays and apodization values |
ITSV20000027A1 (en) * | 2000-06-22 | 2001-12-22 | Esaote Spa | METHOD AND MACHINE FOR THE ACQUISITION OF ECHOGRAPHIC IMAGES IN PARTICULAR OF THE THREE-DIMENSIONAL TYPE AS WELL AS THE ACQUISITION PROBE |
US7744532B2 (en) * | 2004-03-31 | 2010-06-29 | Siemens Medical Solutions Usa, Inc. | Coherence factor adaptive ultrasound imaging methods and systems |
SG11201706953YA (en) * | 2015-02-25 | 2017-09-28 | Decision Sciences Medical Company Llc | Acoustic signal transmission couplants and coupling mediums |
JP2018011927A (en) * | 2016-07-08 | 2018-01-25 | キヤノン株式会社 | Control device, control method, control system, and program |
JP6730919B2 (en) * | 2016-12-12 | 2020-07-29 | 株式会社日立製作所 | Ultrasonic CT device |
US10705210B2 (en) * | 2017-05-31 | 2020-07-07 | B-K Medical Aps | Three-dimensional (3-D) imaging with a row-column addressed (RCA) transducer array using synthetic aperture sequential beamforming (SASB) |
EP3701877B1 (en) * | 2017-10-24 | 2022-06-22 | Lily Medtech Inc. | Ultrasound diagnostic system and ultrasound diagnostic method |
WO2021016767A1 (en) * | 2019-07-26 | 2021-02-04 | 深圳先进技术研究院 | Ultrasonic endoscope probe and ultrasonic endoscope system |
US11627933B2 (en) * | 2019-10-30 | 2023-04-18 | Worcester Polytechnic Institute | Ring-arrayed forward-viewing ultrasonic imaging system and method with needle guidance and image reconstruction |
-
2021
- 2021-06-29 CN CN202110725297.8A patent/CN113552573B/en active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1998002761A1 (en) * | 1996-07-11 | 1998-01-22 | Science Applications International Corporation | Terrain elevation measurement by interferometric synthetic aperture radar (ifsar) |
CN105411626A (en) * | 2015-12-24 | 2016-03-23 | 华中科技大学 | Ultrasonic CT-based synthetic aperture imaging method and system |
CN109281651A (en) * | 2017-07-19 | 2019-01-29 | 中国科学院声学研究所 | A kind of ultrasonic borehole wall imaging method applied to cylinder supersonic array |
CN110179495A (en) * | 2019-04-23 | 2019-08-30 | 华中科技大学 | Ultrasonic tomography calculation optimization method and system based on distributed system |
WO2021004076A1 (en) * | 2019-07-05 | 2021-01-14 | 山东大学 | Ai chip-based conformal wearable biological information monitoring device and system |
CN110441398A (en) * | 2019-07-17 | 2019-11-12 | 复旦大学 | A kind of synthetic aperture ultrasonic imaging method based on multilayer dielectricity velocity of sound model |
CN112764040A (en) * | 2019-11-01 | 2021-05-07 | 复旦大学 | Synthetic aperture beam forming method based on ray theory phase correction |
CN110974293A (en) * | 2019-12-11 | 2020-04-10 | 华中科技大学 | Synthetic aperture imaging method based on C-type probe |
CN112138970A (en) * | 2020-09-22 | 2020-12-29 | 深圳市赛禾医疗技术有限公司 | Ultrasonic forward loop array transceiver |
Non-Patent Citations (3)
Title |
---|
Open-source Gauss-Newton-based methods for refraction-corrected ultrasound computed tomography;Rehman Ali;《PROCEEDINGS OF SPIE》;20190315;第1-14页 * |
基于内窥超声环阵换能器的快速成像方法研究;谭清源等;《集成技术》;20200515(第03期);第44-45页 * |
超声内窥合成孔径成像技术的研究;郁道银等;《中国激光》;20101110(第11期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113552573A (en) | 2021-10-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104688271B (en) | Ultrasonic imaging method and ultrasonic imaging device by synthetic focusing | |
CN101190134B (en) | Method and device for transmitting and receiving multiple wave beams in ultrasound wave diagnosis system | |
WO2019214127A1 (en) | Transcranial three-dimensional cerebrovascular compound imaging method and system | |
EP3581961A1 (en) | Method and apparatus for ultrasound imaging with improved beamforming | |
CN103181779B (en) | Ultrasonic energy conversion device and ultrasonic imaging system and method | |
CN106680825B (en) | A kind of acoustic array imaging system and method | |
CN110811688B (en) | Ultrafast ultrasonic Doppler blood flow estimation method for multi-angle plane wave repeated compounding | |
JP2000232978A (en) | Ultrasonic image pickup for optimizing image quality in region of interest | |
CN102697524B (en) | Full-focus ultrasonic imaging method and application of method in blood flow imaging | |
US6423004B1 (en) | Real-time ultrasound spatial compounding using multiple angles of view | |
CN104688273A (en) | Ultra high speed ultrasonic imaging device and method based on central processing unit (CPU) + graphic processing unit (GPU) isomeric framework | |
US20220167947A1 (en) | Methods and systems for acquiring composite 3d ultrasound images | |
CN114129185B (en) | Beam forming method, ultrasonic imaging method, device and equipment | |
EP2609866A1 (en) | Providing motion mode image in ultrasound system | |
CN106683083A (en) | Image processing method and device of anal sphincter and ultrasonic device | |
CN105997147A (en) | Ultrasonic pulse Doppler imaging method and device | |
CN113552573B (en) | Rapid imaging algorithm based on ultrasonic ring array synthetic aperture receiving | |
CN103829974A (en) | Composite energy Doppler blood flow imaging method and system | |
US20100191115A1 (en) | Ultrasound imaging system and method | |
CN110037740A (en) | System and method for ultrasonic fluid imaging | |
CN113573645B (en) | Method and system for adjusting field of view of an ultrasound probe | |
CN107997784A (en) | A kind of ultrasonic beam synthetic method and system based on velocity of sound adaptive correction | |
CN114966684A (en) | Fast decomposition projection imaging method for group multi-base synthetic aperture radar | |
CN112826529B (en) | Ultrasonic space compounding method and device based on right trapezoid | |
CN114185047A (en) | Bistatic SAR moving target refocusing method based on optimal polar coordinate transformation |
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 |