CN114859359B - Time reversal imaging method, system, equipment and storage medium based on ultrasonic subarray - Google Patents
Time reversal imaging method, system, equipment and storage medium based on ultrasonic subarray Download PDFInfo
- Publication number
- CN114859359B CN114859359B CN202210364433.XA CN202210364433A CN114859359B CN 114859359 B CN114859359 B CN 114859359B CN 202210364433 A CN202210364433 A CN 202210364433A CN 114859359 B CN114859359 B CN 114859359B
- Authority
- CN
- China
- Prior art keywords
- ultrasonic
- array
- ultrasound
- subarray
- imaging
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
The invention discloses a time reversal imaging method and system based on an ultrasonic subarray, equipment and a storage medium. Compared with the existing time reversal imaging method, only part of the array elements are adopted as the working array elements, the number of the working array elements is greatly reduced, the number of time domain signals and the data storage quantity are greatly reduced, different imaging results are compared and analyzed through a multi-parameter comparison function, the imaging results are analyzed from multiple angles such as imaging time, imaging quality and data processing quantity, the optimal ultrasonic subarray division scheme is screened out, the imaging quality is guaranteed, the imaging speed is accelerated, and the requirement on the detection speed in actual detection is well met.
Description
Technical Field
The present invention relates to the field of ultrasound imaging technologies, and in particular, to a time reversal imaging method and system based on an ultrasound sub-array, an electronic device, and a computer-readable storage medium.
Background
The ultrasonic phased array imaging detection technology is an advanced nondestructive detection technology, has the advantages of high detection efficiency, flexible detection mode, reliable detection result and the like, and is applied to industrial departments such as aerospace, nuclear industry, infrastructure and the like. The ultrasonic imaging method is the key of the ultrasonic phased array imaging detection technology, and the same ultrasonic time domain signal is processed by different imaging methods, so that ultrasonic images with different qualities can be obtained. The time reversal is an important imaging method, typically represented by a time reversal multi-signal classification method TR-MUSIC, due to the super-resolution characteristic of the time reversal multi-signal classification method, the TR-MUSIC attracts attention in the fields of ultrasonic nondestructive testing, medical imaging and the like, is widely applied to the aspects of structural integrity detection, disease diagnosis and treatment and the like, and has obvious social and economic benefits.
However, the current time reversal imaging method utilizes full matrix acquisition, that is, all array elements of an array element sensor are utilized to acquire ultrasonic array data, and the data is post-processed to acquire an ultrasonic image. All array elements of the array sensor are required to be used for exciting and receiving ultrasonic signals, so that data redundancy is caused, the imaging speed is influenced, and the requirement on the detection speed in actual detection cannot be met. Although some proposals for sparse arrays have been proposed in recent research, the number of schemes to be considered for optimizing the configuration of sparse arrays is too large, which increases the burden of data storage and calculation, and still cannot meet the requirement for detection speed in actual detection.
Disclosure of Invention
The invention provides a time reversal imaging method and system based on an ultrasonic subarray, electronic equipment and a computer readable storage medium, and aims to solve the technical problems of data redundancy and low imaging speed in the existing time reversal imaging method that all array elements of an array element sensor are used for acquiring ultrasonic array data.
According to an aspect of the present invention, there is provided an ultrasound subarray-based time reversal imaging method, comprising:
designing a plurality of ultrasonic subarray division schemes, and selecting one array element from each ultrasonic subarray as a working array element;
for different ultrasonic subarray division schemes, acquiring ultrasonic array data by using a full-matrix acquisition method respectively, and obtaining ultrasonic images by using a time reversal multi-signal classification method;
defining a multi-parameter comparison function for performing comparative analysis on different ultrasonic sub-array division schemes;
and performing contrast analysis on the multiple ultrasonic images by using the defined multi-parameter comparison function, screening out an optimal ultrasonic sub-array division scheme, and outputting the ultrasonic images obtained by using the optimal ultrasonic sub-array division scheme as final ultrasonic images.
Further, the ultrasound sub-array partitioning scheme specifically includes:
setting the number of array elements contained in each ultrasonic subarray as e _ N, judging whether N/e _ N is an integer or not for the ultrasonic array containing N array elements, if so, dividing the N array elements into N ultrasonic subarrays, wherein N = N/e _ N, if not, removing a remainder, taking 1-e _ N as a first ultrasonic subarray, and repeating the steps until the last array element N × e _ N is obtained, and N × e _ N is less than N.
Furthermore, the number of the array elements contained in the ultrasonic subarrays is odd, and the middle array element is selected from each ultrasonic subarray as a working array element.
Further, the entire ultrasound array is taken as one reference ultrasound sub-array.
Further, the multi-parameter comparison function is: y (e, d _ n, t, SA _ API, SA _ w) = f (e _ n), where e represents the number of active array elements, d _ n represents the number of time domain signals included in the ultrasound array data, t represents the imaging time, SA _ API represents the array evaluation factor, SA _ w represents the-6 dB main lobe width, SA _ API = a (-6 dB)/λ 2 A (-6 dB) represents an imaging region with intensity greater than-6 dB in the ultrasonic image, λ represents the ultrasonic wavelength, SA _ w = | w1_ x-w2_ x |, w1 and w2 are the intersection points of the-6 dB threshold line and the transverse intensity curve, and w1_ x and w2_ x are the x coordinate values thereof.
Further, the process of performing comparative analysis on the multiple ultrasound images by using the defined multi-parameter comparison function to screen out the optimal ultrasound sub-array partition scheme specifically comprises the following steps:
and (3) for each ultrasonic image, respectively calculating an array evaluation factor and a-6 dB main lobe width, and screening out the minimum ultrasonic image and the minimum ultrasonic image as an optimal ultrasonic subarray division scheme.
In addition, the invention also provides a time reversal imaging system based on the ultrasonic subarray, which comprises:
the ultrasonic subarray division module is used for designing various ultrasonic subarray division schemes and selecting one array element from each ultrasonic subarray as a working array element;
the ultrasonic imaging module is used for acquiring ultrasonic array data by respectively utilizing a full matrix acquisition method and obtaining ultrasonic images by utilizing a time reversal multi-signal classification method for different ultrasonic subarray division schemes;
the function definition module is used for defining a multi-parameter comparison function for performing comparison analysis on different ultrasonic sub-array division schemes;
and the contrast analysis module is used for performing contrast analysis on the multiple ultrasonic images by using the defined multi-parameter comparison function, screening out an optimal ultrasonic sub-array division scheme, and outputting the ultrasonic images obtained by using the optimal ultrasonic sub-array division scheme as final ultrasonic images.
Further, the multi-parameter comparison function is: y (e, d _ n, t, SA _ API, SA _ w) = f (e _ n), where e represents the number of working array elements, d _ n represents the number of time domain signals included in the ultrasound array data, t represents the imaging time, SA _ API represents the array evaluation factor, SA _ w represents the-6 dB main lobe width, SA _ API = a (-6 dB)/λ 2 A (-6 dB) represents an imaging region with intensity greater than-6 dB in the ultrasonic image, λ represents the ultrasonic wavelength, SA _ w = | w1_ x-w2_ x |, w1 and w2 are the intersection points of the-6 dB threshold line and the transverse intensity curve, and w1_ x and w2_ x are the x coordinate values thereof.
In addition, the present invention also provides an apparatus comprising a processor and a memory, wherein the memory stores a computer program, and the processor is used for executing the steps of the method by calling the computer program stored in the memory.
In addition, the present invention also provides a computer readable storage medium for storing a computer program for time-reversal imaging based on ultrasound sub-arrays, which computer program, when running on a computer, performs the steps of the method as described above.
The invention has the following effects:
the time reversal imaging method based on the ultrasonic subarrays comprises the steps of firstly designing various ultrasonic subarray division schemes, dividing an ultrasonic array of an array element sensor into a plurality of ultrasonic subarrays with the same number of array elements in each ultrasonic subarray division scheme, and selecting one array element from each ultrasonic subarray to serve as a working array element for exciting and receiving ultrasonic signals. And then, respectively carrying out time reversal imaging to obtain ultrasonic images, and defining a multi-parameter comparison function to carry out comparison analysis on different ultrasonic sub-array division schemes. And finally, screening out an optimal ultrasonic subarray division scheme through comparison analysis, and outputting an ultrasonic image corresponding to the optimal ultrasonic subarray division scheme. Compared with the existing time reversal imaging method, the time reversal imaging method based on the ultrasonic subarray only needs to adopt partial array elements as the working array elements, the number of the working array elements is greatly reduced, the number of time domain signals and data storage capacity are greatly reduced, imaging results of different ultrasonic subarray division schemes are compared and analyzed by defining a multi-parameter comparison function, the imaging results can be analyzed from multiple angles such as imaging time, imaging quality and processing data quantity, the optimal ultrasonic subarray division scheme is screened out, the imaging quality is guaranteed, the imaging speed is accelerated, and the requirement for the detection speed in actual detection can be well met.
In addition, the time reversal imaging system, the time reversal imaging device and the time reversal imaging storage medium based on the ultrasonic subarray also have the advantages.
In addition to the objects, features and advantages described above, other objects, features and advantages of the present invention are also provided. The present invention will be described in further detail below with reference to the drawings.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of this application, illustrate embodiments of the invention and, together with the description, serve to explain the invention and not to limit the invention. In the drawings:
fig. 1 is a flow chart diagram of a time reversal imaging method based on an ultrasound sub-array according to a preferred embodiment of the present invention.
Figure 2 is a schematic diagram of the ultrasound sub-array partitioning logic of the preferred embodiment of the present invention.
Figure 3 is a schematic illustration of a preferred embodiment of the present invention for acquiring ultrasound array data using full matrix acquisition.
Fig. 4 is a schematic diagram of ultrasound images corresponding to different ultrasound sub-array division schemes in a preferred embodiment of the present invention.
Figure 5 is a graphical representation of the results of ultrasound imaging and the transverse intensity profile through the image peak for the baseline version of the preferred embodiment of the present invention.
Fig. 6 is a logic diagram of the screening of the optimal ultrasound sub-array partitioning scheme in the preferred embodiment of the present invention.
Fig. 7 is a schematic diagram of a multi-parameter curve drawn based on image construction parameters and image evaluation parameters under different ultrasound sub-array partition schemes in a preferred embodiment of the present invention.
Fig. 8 is a schematic block diagram of an ultrasound sub-array based time reversal imaging system according to another embodiment of the present invention.
Detailed Description
The embodiments of the invention will be described in detail below with reference to the accompanying drawings, but the invention can be embodied in many different forms, which are defined and covered by the following description.
As shown in fig. 1, a preferred embodiment of the present invention provides a time reversal imaging method based on ultrasound sub-arrays, comprising the following:
step S1: designing a plurality of ultrasonic subarray division schemes, and selecting one array element from each ultrasonic subarray as a working array element;
step S2: for different ultrasonic subarray division schemes, acquiring ultrasonic array data by using a full-matrix acquisition method respectively, and obtaining ultrasonic images by using a time reversal multi-signal classification method;
and step S3: defining a multi-parameter comparison function for performing comparative analysis on different ultrasonic sub-array division schemes;
and step S4: and performing contrast analysis on the multiple ultrasonic images by using the defined multi-parameter comparison function, screening out an optimal ultrasonic sub-array division scheme, and outputting the ultrasonic images obtained by using the optimal ultrasonic sub-array division scheme as final ultrasonic images.
It can be understood that, in the time reversal imaging method based on the ultrasound subarrays of this embodiment, a plurality of ultrasound subarray division schemes are first designed, in each ultrasound subarray division scheme, the ultrasound array of the array element sensor is divided into a plurality of ultrasound subarrays including the same number of array elements, and one array element is selected from each ultrasound subarray as a working array element for exciting and receiving an ultrasound signal. And then, respectively carrying out time reversal imaging to obtain ultrasonic images, and defining a multi-parameter comparison function to carry out comparison analysis on different ultrasonic sub-array division schemes. And finally, screening out an optimal ultrasonic subarray division scheme through comparison analysis, and outputting an ultrasonic image corresponding to the optimal ultrasonic subarray division scheme. Compared with the existing time reversal imaging method, the time reversal imaging method based on the ultrasonic subarray only needs to adopt partial array elements as the working array elements, the number of the working array elements is greatly reduced, the number of time domain signals and data storage capacity are greatly reduced, imaging results of different ultrasonic subarray division schemes are compared and analyzed by defining a multi-parameter comparison function, the imaging results can be analyzed from multiple angles such as imaging time, imaging quality and processing data quantity, the optimal ultrasonic subarray division scheme is screened out, the imaging quality is guaranteed, the imaging speed is accelerated, and the requirement for the detection speed in actual detection can be well met.
It can be understood that, in the step S1, the ultrasound sub-array division scheme specifically includes:
setting the number of array elements contained in each ultrasonic subarray as e _ N, judging whether N/e _ N is an integer or not for the ultrasonic array containing N array elements, if so, dividing the N array elements into N ultrasonic subarrays, N = N/e _ N, if not, removing remainder, taking 1-e _ N as a first ultrasonic subarray, and repeating the steps until the last array element N × e _ N is less than N.
Optionally, the number of array elements included in the ultrasound sub-arrays is odd, and a middle array element is selected from each ultrasound sub-array as an operating array element.
Optionally, the entire ultrasound array is taken as one reference ultrasound sub-array.
In an embodiment of the present invention, for a linear array including N array elements, the linear array is divided into N sub-arrays, and e _ N is the number of array elements included in each sub-array. As shown in fig. 2, the ultrasonic sub-array division logic diagram firstly determines whether N/e _ N is an integer, if so, 1-e _ N are sub-arrays 1, and so on until the last array element N, that is, N = N/e _ N; if not, removing the remainder, wherein 1-e _ N are ultrasonic subarrays 1, and repeating the steps until the last array element N × e _ N is less than N. And after the ultrasonic subarrays are divided, selecting the middle array element of each ultrasonic subarray as a working array element for exciting and receiving ultrasonic signals.
Specifically, a 64-element linear array sensor is considered, wherein N =64, and array elements are numbered from left to right as 1-64; e _ n selects 1,3,5,7,9 and 11 respectively, obtains the corresponding ultrasonic sub-array division scheme, and selects the working array element.
The first scheme is as follows: e _ n =1, all array elements are working array elements, and the scheme is a reference scheme and is used for comparing with other schemes;
scheme II: e _ n =3, and the number of the working array elements is 2,5,8, 11, 14, 17, 20, 23, 26, 29, 32, 35, 38, 41, 44, 47, 50, 53, 56, 59, 62;
the third scheme is as follows: e _ n =5, and the number of the working array elements is 3,8, 13, 18, 23, 28, 33, 38, 43, 48, 53, 58;
the scheme four is as follows: e _ n =7, the number of the working array elements is 4, 11, 18, 25, 32, 39, 46, 53, 60;
and a fifth scheme: e _ n =9, the number of the working array elements is 5, 14, 23, 32, 41, 50, 59;
scheme six: e _ n =11, and the number of the working array elements is 6, 17, 28, 39, 50.
Alternatively, in other embodiments of the present invention, the number of array elements included in each ultrasound sub-array may also be an even number, and the array elements in the same position in each ultrasound sub-array are selected as the working array elements, for example, if the number of array elements included in each ultrasound sub-array is 10, the working array elements are numbered as follows: 2. 12, 22, 32, 42, 52. It can be understood that the number of the array elements included in the ultrasonic sub-array and the selection of the working array elements can be adjusted as required, and the selected plurality of working array elements are only required to be uniformly distributed.
It can be understood that, in step S1, the ultrasonic array of the array element sensor is divided into a plurality of ultrasonic sub-arrays, and one array element is selected from each ultrasonic sub-array as a working array element for exciting and receiving ultrasonic signals, and subsequently, only the time domain signals of the working array element need to be stored and processed, so that the number of time sequence signals and the data storage amount are greatly reduced, and the imaging speed is improved.
It can be understood that, in the step S2, as shown in fig. 3, for each ultrasound sub-array division scheme, the ultrasound array data is acquired by using a full matrix acquisition method FMC, and the time domain signal is processed by using a time reversal multi-signal classification method TR-MUSIC, so as to obtain an ultrasound image. The specific FMC method and TR-MUSIC method belong to the prior art, and are not described herein again.
Specifically, the object to be measured is a metal block structure, and a through hole with a diameter of 1mm is processed in the metal block structure as an imaging target, and the size of the through hole is smaller than the wavelength of ultrasonic waves, so that the through hole can be regarded as a point target. The FMC acquires the ultrasonic array data by utilizing any working array element combination, and the number of time domain signals contained in the ultrasonic array data under different ultrasonic sub-array division schemes is as follows:
the first scheme is as follows: e _ n =1, 64 × 64;
scheme two is as follows: e _ n =3, 21 × 21;
and a third scheme is as follows: e _ n =5, 12 × 12;
and the scheme is as follows: e _ n =7,9 x 9;
and a fifth scheme: e _ n =9,7 × 7;
scheme six: e _ n =11,5 × 5.
Meanwhile, the TR-MUSIC method is used to process the ultrasonic array data to obtain ultrasonic images under different ultrasonic sub-array division schemes, as shown in fig. 4.
It is understood that, in the step S3, the multi-parameter comparison function is: y (e, d _ n, t, SA _ API, SA _ w) = f (e _ n), where e represents the number of working array elements, d _ n represents the number of time domain signals included in the ultrasound array data, t represents the imaging time, SA _ API represents the array evaluation factor, SA _ w represents the-6 dB main lobe width, SA _ API = a (-6 dB)/λ 2 Where a (-6 dB) represents the imaging region in the ultrasound image with intensity greater than-6 dB, λ represents the ultrasound wavelength, and SA _ w = | w1_ x-w2_ x |, as shown in fig. 5, w1 and w2 are the intersection points of the-6 dB threshold line and the transverse intensity curve taken from the ultrasound imaging result and passing through the image peak, and w1_ x, w2_ x are their x-coordinate values. Wherein e, d _ n and t are image construction parameters, and SA _ API and SA _ w are image evaluation parameters.
It can be understood that, in the step S4, the process of performing comparative analysis on the plurality of ultrasound images by using the defined multi-parameter comparison function to screen out the optimal ultrasound sub-array partition scheme specifically includes:
and (3) for each ultrasonic image, respectively calculating an array evaluation factor and a-6 dB main lobe width, and screening the minimum value of the array evaluation factor and the-6 dB main lobe width as an optimal ultrasonic subarray division scheme.
Specifically, as shown in fig. 6, for different ultrasound sub-array division schemes, image construction parameters e, d _ n, and t may be obtained, an ultrasound image may be obtained through step S2, image evaluation parameters SA _ API and SA _ w may be obtained through the ultrasound image, and at the same time, e0, d _ n0, t0, SA _ API0, and SA _ w0 are obtained through a reference scheme, and it is determined whether SA _ API/SA _ API0& SA _ w/SA _ w0 is minimum, and if so, e/e0, d _ n/d _ n0, and t/t0 are calculated.
The parameters under different ultrasound subarray division schemes are as follows:
the first scheme comprises the following steps: e _ n =1, e0=64, d_n0 =4096, t0=26.11s, sa _api0=0.235, sa _w0=0.2603mm;
scheme II: e _ n =3, e =21, d _n =441, t =16.33s, sa _api =0.2786, sa _w =0.2903mm;
and a third scheme is as follows: e _ n =5, e =12, d _n =144, t =14.29s, sa _api =0.2502, sa _w =0.2603mm;
the scheme four is as follows: e _ n =7, e =9, d _n =81, t =13.73s, sa _api =0.3791, sa _w =0.3303mm;
and a fifth scheme: e _ n =9, e =7, d _n =49, t =13.17s, sa _api =0.2994, sa _w =0.2903mm;
scheme six: e _ n =11, e =5, d _n =25, t =13.12s, sa _api =0.796, sa _w =0.4404mm.
Then, a multi-parameter graph is drawn based on the obtained parameters, as shown in fig. 7. As can be seen from the multi-parameter graph or the above-mentioned parameter data, when e _ n =5, SA _ API/SA _ API0& SA _ w/SA _ w0 is minimum, and the ultrasound sub-array division scheme at this time is selected as an optimal scheme. Correspondingly, e/e0=0.1875, d _n/d _ n0=0.0352, t/t0=0.5473, that is, under the optimal partitioning scheme of the ultrasound sub-array, compared with the reference scheme, the number of working array elements is 18.75%, which is reduced by 81.25%, the number of time domain signals included in the ultrasound array data is 3.52%, which is reduced by 96.48%, the imaging time is 54.73%, which is reduced by 45.27%, the data storage capacity is greatly reduced, the imaging speed is increased, and the actual detection requirement can be met. It can be understood that the purpose of calculating the image construction parameters of each ultrasound sub-array division scheme is to compare the image construction parameters with a reference scheme so as to embody the advantage of the ultrasound sub-array division scheme in the aspect of improving the detection speed, and the purpose of calculating the image evaluation parameters is to screen out an optimal scheme.
In addition, as shown in fig. 8, another embodiment of the present invention further provides an ultrasound subarray-based time-reversal imaging system, preferably using the imaging method described above, the system comprising:
the ultrasonic subarray division module is used for designing various ultrasonic subarray division schemes and selecting one array element from each ultrasonic subarray as a working array element;
the ultrasonic imaging module is used for acquiring ultrasonic array data by respectively utilizing a full matrix acquisition method and obtaining ultrasonic images by utilizing a time reversal multi-signal classification method for different ultrasonic subarray division schemes;
the function definition module is used for defining a multi-parameter comparison function for performing comparison analysis on different ultrasonic sub-array division schemes;
and the contrast analysis module is used for performing contrast analysis on the multiple ultrasonic images by using the defined multi-parameter comparison function, screening out an optimal ultrasonic sub-array division scheme, and outputting the ultrasonic images obtained by using the optimal ultrasonic sub-array division scheme as final ultrasonic images.
It can be understood that, in the time reversal imaging system based on the ultrasound subarray of this embodiment, a plurality of ultrasound subarray division schemes are first designed, in each ultrasound subarray division scheme, the ultrasound array of the array element sensor is divided into a plurality of ultrasound subarrays containing the same number of array elements, and one array element is selected from each ultrasound subarray as a working array element for exciting and receiving an ultrasound signal. And then, respectively carrying out time reversal imaging to obtain ultrasonic images, and defining a multi-parameter comparison function to carry out comparison analysis on different ultrasonic sub-array division schemes. And finally, screening out an optimal ultrasonic subarray division scheme through comparison analysis, and outputting an ultrasonic image corresponding to the scheme. Compared with the existing time reversal imaging system, the time reversal imaging system based on the ultrasonic subarray only needs to adopt partial array elements as the working array elements, the number of the working array elements is greatly reduced, the number of time domain signals and the data storage quantity are greatly reduced, imaging results of different ultrasonic subarray division schemes are compared and analyzed by defining a multi-parameter comparison function, the imaging results can be analyzed from multiple angles such as imaging time, imaging quality and data processing quantity, the optimal ultrasonic subarray division scheme is screened out, the imaging quality is ensured, the imaging speed is accelerated, and the requirement for the detection speed in actual detection can be well met.
Wherein the multi-parameter comparison function is: y (e, d _ n, t, SA _ API, SA _ w) = f (e _ n), where e represents the number of active array elements, d _ n represents the number of time domain signals included in the ultrasound array data, t represents the imaging time, SA _ API represents the array evaluation factor, SA_w represents-6 dB main lobe width, SA _ API = A (-6 dB)/λ 2 A (-6 dB) represents an imaging region with intensity greater than-6 dB in the ultrasonic image, λ represents the ultrasonic wavelength, SA _ w = | w1_ x-w2_ x |, w1 and w2 are the intersection points of the-6 dB threshold line and the transverse intensity curve, and w1_ x and w2_ x are the x coordinate values thereof.
It can be understood that the specific working principle and working process of each module of this embodiment correspond to each step of the above method embodiment, and therefore, detailed description is not repeated here, and only the above method embodiment is referred to.
In addition, another embodiment of the present invention further provides an apparatus, which includes a processor and a memory, wherein the memory stores a computer program, and the processor is used for executing the steps of the method described above by calling the computer program stored in the memory.
In addition, another embodiment of the present invention also provides a computer readable storage medium for storing a computer program for time-reversal imaging based on ultrasound sub-arrays, which when run on a computer performs the steps of the method as described above.
Typical forms of computer-readable storage media include: floppy disk (floppy disk), flexible disk (flexible disk), hard disk, magnetic tape, any of the other magnetic media, CD-ROM, any of the remaining optical media, punch cards (punch cards), paper tape (paper tape), any of the remaining physical media with patterns of holes, random Access Memory (RAM), programmable Read Only Memory (PROM), erasable Programmable Read Only Memory (EPROM), FLASH erasable programmable read only memory (FLASH-EPROM), any of the remaining memory chips or cartridges, or any of the remaining media readable by a computer. The instructions may further be transmitted or received by a transmission medium. The term transmission medium may include any tangible or intangible medium that is operable to store, encode, or carry instructions for execution by the machine, and includes digital or analog communications signals or intangible medium that facilitates communication of the instructions. Transmission media include coaxial cables, copper wire and fiber optics, including the wires that comprise a bus for transmitting a computer data signal.
The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (8)
1. A time reversal imaging method based on ultrasonic subarrays is characterized by comprising the following steps:
designing a plurality of ultrasonic subarray division schemes, and selecting one array element from each ultrasonic subarray as a working array element;
for different ultrasonic subarray division schemes, acquiring ultrasonic array data by using a full-matrix acquisition method, and obtaining ultrasonic images by using a time reversal multi-signal classification method;
defining a multi-parameter comparison function for performing comparison analysis on different ultrasonic sub-array partition schemes;
carrying out contrast analysis on the multiple ultrasonic images by using a defined multi-parameter comparison function, screening out an optimal ultrasonic sub-array division scheme, and outputting the ultrasonic images obtained by using the optimal ultrasonic sub-array division scheme as final ultrasonic images;
the multi-parameter comparison function is: y (e, d _ n, t, SA _ API, SA _ w) = f (e _ n), where e represents the number of working array elements, d _ n represents the number of time domain signals included in the ultrasound array data, t represents the imaging time, SA _ API represents the array evaluation factor, SA _ w represents the-6 dB main lobe width, SA _ API = a (-6 dB)/λ 2 A (-6 dB) represents an imaging region with an intensity greater than-6 dB in an ultrasonic image, λ represents the wavelength of the ultrasonic wave, SA _ w = | w1_ x-w2_ x |, w1 and w2 are intersections of a-6 dB threshold line and a transverse intensity curve, and w1_ x and w2_ x are x coordinate values thereof.
2. The ultrasound sub-array based time-reversal imaging method of claim 1, wherein the ultrasound sub-array division scheme is specifically:
setting the number of array elements contained in each ultrasonic subarray as e _ N, judging whether N/e _ N is an integer or not for the ultrasonic array containing N array elements, if so, dividing the N array elements into N ultrasonic subarrays, N = N/e _ N, if not, removing remainder, taking 1-e _ N as a first ultrasonic subarray, and repeating the steps until the last array element N × e _ N is less than N.
3. The ultrasound sub-array based time-reversal imaging method of claim 2, wherein the ultrasound sub-arrays include an odd number of array elements, and wherein a middle array element is selected from each ultrasound sub-array as an active array element.
4. The ultrasound sub-array based time-reversal imaging method of claim 3, wherein the entire ultrasound array is taken as one reference ultrasound sub-array.
5. The ultrasound sub-array based time reversal imaging method according to claim 1, wherein the process of performing a comparative analysis on a plurality of ultrasound images by using a defined multi-parameter comparison function to screen out an optimal ultrasound sub-array partition scheme specifically includes:
and (3) for each ultrasonic image, respectively calculating an array evaluation factor and a-6 dB main lobe width, and screening out the minimum ultrasonic image and the minimum ultrasonic image as an optimal ultrasonic subarray division scheme.
6. An ultrasound subarray based time reversed imaging system comprising:
the ultrasonic subarray division module is used for designing various ultrasonic subarray division schemes and selecting one array element from each ultrasonic subarray as a working array element;
the ultrasonic imaging module is used for acquiring ultrasonic array data by respectively utilizing a full matrix acquisition method and obtaining ultrasonic images by utilizing a time reversal multi-signal classification method for different ultrasonic subarray division schemes;
the function definition module is used for defining a multi-parameter comparison function for performing comparison analysis on different ultrasonic sub-array division schemes;
the contrast analysis module is used for performing contrast analysis on the multiple ultrasonic images by using a defined multi-parameter comparison function, screening out an optimal ultrasonic sub-array division scheme, and outputting the ultrasonic images obtained by using the optimal ultrasonic sub-array division scheme as final ultrasonic images;
the multi-parameter comparison function is: y (e, d _ n, t, SA _ API, SA _ w) = f (e _ n), where e represents the number of working array elements, d _ n represents the number of time domain signals included in the ultrasound array data, t represents the imaging time, SA _ API represents the array evaluation factor, SA _ w represents the-6 dB main lobe width, SA _ API = a (-6 dB)/λ 2 A (-6 dB) represents an imaging region with an intensity greater than-6 dB in an ultrasonic image, λ represents the wavelength of the ultrasonic wave, SA _ w = | w1_ x-w2_ x |, w1 and w2 are intersections of a-6 dB threshold line and a transverse intensity curve, and w1_ x and w2_ x are x coordinate values thereof.
7. An apparatus, comprising a processor and a memory, the memory having stored therein a computer program, the processor being configured to perform the steps of the method according to any one of claims 1 to 5 by calling the computer program stored in the memory.
8. A computer-readable storage medium for storing a computer program for time-reversal imaging based on ultrasound sub-arrays, characterized in that the computer program, when run on a computer, performs the steps of the method according to any one of claims 1 to 5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210364433.XA CN114859359B (en) | 2022-04-07 | 2022-04-07 | Time reversal imaging method, system, equipment and storage medium based on ultrasonic subarray |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210364433.XA CN114859359B (en) | 2022-04-07 | 2022-04-07 | Time reversal imaging method, system, equipment and storage medium based on ultrasonic subarray |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114859359A CN114859359A (en) | 2022-08-05 |
| CN114859359B true CN114859359B (en) | 2022-12-02 |
Family
ID=82630310
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210364433.XA Active CN114859359B (en) | 2022-04-07 | 2022-04-07 | Time reversal imaging method, system, equipment and storage medium based on ultrasonic subarray |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114859359B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104280455A (en) * | 2014-09-29 | 2015-01-14 | 北京工业大学 | Ultrasonic scattering coefficient optimal computation method for crack direction recognition |
| CN104730146A (en) * | 2015-03-11 | 2015-06-24 | 北京工业大学 | Ultrasonic array composite imaging method for quantitative evaluation of plate structure cracks |
| CN107802286A (en) * | 2017-11-09 | 2018-03-16 | 中国人民解放军国防科技大学 | Ultrasonic imaging method and system based on multi-frequency time reversal technology |
| CN108061760A (en) * | 2018-01-29 | 2018-05-22 | 中国人民解放军国防科技大学 | Ultrasonic rapid super-resolution imaging method and device, storage medium and imaging system |
| CN109490419A (en) * | 2019-01-04 | 2019-03-19 | 中国航发北京航空材料研究院 | A kind of acoustic beam calibration method of total focus imaging |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9759690B2 (en) * | 2011-02-24 | 2017-09-12 | General Electric Company | Method and system for nondestructive ultrasound testing |
| TW201325556A (en) * | 2011-12-28 | 2013-07-01 | Ind Tech Res Inst | Ultrasound transducer and ultrasound image system and image method |
| US20140155751A1 (en) * | 2012-12-05 | 2014-06-05 | Kabushiki Kaisha Toshiba | Method and system for element-by-element flexible subarray beamforming |
-
2022
- 2022-04-07 CN CN202210364433.XA patent/CN114859359B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104280455A (en) * | 2014-09-29 | 2015-01-14 | 北京工业大学 | Ultrasonic scattering coefficient optimal computation method for crack direction recognition |
| CN104730146A (en) * | 2015-03-11 | 2015-06-24 | 北京工业大学 | Ultrasonic array composite imaging method for quantitative evaluation of plate structure cracks |
| CN107802286A (en) * | 2017-11-09 | 2018-03-16 | 中国人民解放军国防科技大学 | Ultrasonic imaging method and system based on multi-frequency time reversal technology |
| CN108061760A (en) * | 2018-01-29 | 2018-05-22 | 中国人民解放军国防科技大学 | Ultrasonic rapid super-resolution imaging method and device, storage medium and imaging system |
| CN109490419A (en) * | 2019-01-04 | 2019-03-19 | 中国航发北京航空材料研究院 | A kind of acoustic beam calibration method of total focus imaging |
Non-Patent Citations (3)
| Title |
|---|
| 基于时间反转声学理论的超声相控阵超分辨率成像方法;樊程广 等;《国防科技大学学报》;20211214;摘要 * |
| 子阵列差分系数加权的超声成像算法;赵巍 等;《声学学报》;20210115;第111-120页 * |
| 板结构裂纹兰姆波阵列复合成像方法研究;焦敬品 等;《仪器仪表学报》;20160315;第593-601页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114859359A (en) | 2022-08-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20210272339A1 (en) | Systems and Methods for Generating and Estimating Unknown and Unacquired Ultrasound Data | |
| CN114859359B (en) | Time reversal imaging method, system, equipment and storage medium based on ultrasonic subarray | |
| CN113504306A (en) | Steel rail defect detection method based on ultrasonic phased array low-rank matrix recovery | |
| CN111317509B (en) | Method and device for generating blood vessel longitudinal section image, diagnostic equipment and storage medium | |
| CN107802286B (en) | Ultrasonic imaging method and system based on multi-frequency time reversal technology | |
| CN117045327B (en) | Ultrasonic puncture needle developing method and device, ultrasonic equipment and storage medium | |
| US9759690B2 (en) | Method and system for nondestructive ultrasound testing | |
| Ordonez et al. | Evaluation of echosounder data preparation strategies for modern machine learning models | |
| CN111243240B (en) | Landslide early warning method and device | |
| KR101496205B1 (en) | Reconstruction of phased array data | |
| CN115082320B (en) | Ultrasonic image super-resolution processing method, ultrasonic device and storage medium | |
| CN116228753A (en) | Tumor prognosis evaluation method, device, computer equipment and storage medium | |
| CN115774055A (en) | Method, device and system for detecting internal defects of integrated circuit | |
| CN113094933B (en) | Ultrasonic damage detection and analysis method based on attention mechanism and application thereof | |
| US5526473A (en) | Method of and apparatus for reducing the size of a display while substantially maintaining its information content | |
| CN111091144B (en) | Image feature point matching method and device based on depth pseudo-twin network | |
| CN113030891A (en) | Optimization method for target point aggregation processing in radar system | |
| CN111062935B (en) | Mammary gland tumor detection method, storage medium and terminal equipment | |
| EP0340549B1 (en) | Image processing method | |
| CN114926468B (en) | Ultrasonic image quality control method, ultrasonic device, and storage medium | |
| CN117783301A (en) | Abnormality recognition method, apparatus, device, and storage medium | |
| CN116712101B (en) | Ultrasound image generation method, device, computer equipment and storage medium | |
| Mehdinia | Multimodal Imaging of Structural Concrete Using Image Fusion and Deep Learning | |
| Abbate et al. | Ultrasonic gauging and wavelet image processing for wear and erosion mapping | |
| CN117237299A (en) | Ultrasonic image detection method and system |
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 |