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 PDF

Info

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
Application number
CN202110725297.8A
Other languages
Chinese (zh)
Other versions
CN113552573A (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.)
Fudan University
Original Assignee
Fudan 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 Fudan University filed Critical Fudan University
Priority to CN202110725297.8A priority Critical patent/CN113552573B/en
Publication of CN113552573A publication Critical patent/CN113552573A/en
Application granted granted Critical
Publication of CN113552573B publication Critical patent/CN113552573B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8997Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using synthetic aperture techniques
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8909Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration
    • G01S15/8915Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array
    • G01S15/8922Short-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

Rapid imaging algorithm based on ultrasonic ring array synthetic aperture receiving
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 button
Figure BDA0003138373080000021
Rotating 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:
Figure BDA0003138373080000022
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:
Figure BDA0003138373080000061
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 button
Figure BDA0003138373080000062
Rotating 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
Figure BDA0003138373080000081
Figure BDA0003138373080000082
The algorithm of the invention has the time complexity of
Figure BDA0003138373080000083
Figure BDA0003138373080000084
Figure BDA0003138373080000085
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 transducer
Figure DEST_PATH_IMAGE001
Each array element in each group
Figure DEST_PATH_IMAGE002
Each adjacent said array element being divided into
Figure 126003DEST_PATH_IMAGE001
An array element group, the calculation comprises the array elements
Figure DEST_PATH_IMAGE003
The 1 st sector area index matrix formed by the array element group and the circle center of the ultrasonic annular array transducer
Figure DEST_PATH_IMAGE004
Step S2, setting circulation value
Figure DEST_PATH_IMAGE005
Step S3, according to the array element
Figure DEST_PATH_IMAGE006
The received ultrasonic echo signal is used for calculating the array element
Figure 96672DEST_PATH_IMAGE006
Each of whichThe array element to index matrix
Figure DEST_PATH_IMAGE007
The 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 algorithm
Figure 610830DEST_PATH_IMAGE007
To correspond to the first
Figure DEST_PATH_IMAGE008
Imaging results of a sector
Figure DEST_PATH_IMAGE009
Step S5, let
Figure DEST_PATH_IMAGE010
And the index matrix is combined
Figure 333935DEST_PATH_IMAGE004
Push button
Figure DEST_PATH_IMAGE011
The index matrix is obtained by rotating and transforming
Figure 882728DEST_PATH_IMAGE007
And then returns to the step S3 until
Figure DEST_PATH_IMAGE012
Step S6, imaging the sector area
Figure DEST_PATH_IMAGE013
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 matrix
Figure 236349DEST_PATH_IMAGE007
Having 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 place
Figure 300120DEST_PATH_IMAGE008
In each sector, the value of the element is 0, which means that the point is located at the second position
Figure 499020DEST_PATH_IMAGE008
Outside the area of the sector of the circle,
in the step S4, for the second step
Figure 902320DEST_PATH_IMAGE008
Within each sector area of the index matrix
Figure 692421DEST_PATH_IMAGE007
At any point where the value of the corresponding element is 1
Figure DEST_PATH_IMAGE014
According to the following formula:
Figure DEST_PATH_IMAGE015
calculate the point
Figure 977909DEST_PATH_IMAGE014
The image luminance value of
Figure 449342DEST_PATH_IMAGE008
The matrix formed by the image luminance values of all the points in the fan-shaped area is the imaging result
Figure 972727DEST_PATH_IMAGE009
In the formula (I), wherein,
Figure DEST_PATH_IMAGE016
as an array element
Figure DEST_PATH_IMAGE017
Transmitting array element
Figure DEST_PATH_IMAGE018
The received signal is at
Figure DEST_PATH_IMAGE019
The magnitude of the amplitude of the time of day,
Figure DEST_PATH_IMAGE020
and
Figure DEST_PATH_IMAGE021
are respectively the array elements
Figure 261626DEST_PATH_IMAGE018
And said array element
Figure 237672DEST_PATH_IMAGE017
To
Figure DEST_PATH_IMAGE022
The 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.
CN202110725297.8A 2021-06-29 2021-06-29 Rapid imaging algorithm based on ultrasonic ring array synthetic aperture receiving Active CN113552573B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (9)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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